Search

Bladder exstrophy (BE) is a rare, lower ventral midline defect with the bladder and part of the urethra exposed. The etiology of BE is unknown but thought to be influenced by genetic variation with more recent studies suggesting a role for rare variants. As such, we conducted paired-end exome sequencing in 26 child/mother/father trios. Three children had rare (allele frequency ≤ 0.0001 in several public databases) inherited variants in TSPAN4, one with a loss-of-function variant and two with missense variants. Two children had loss-of-function variants in TUBE1. Four children had rare missense or nonsense variants (one per child) in WNT3, CRKL, MYH9, or LZTR1, genes previously associated with BE. We detected 17 de novo missense variants in 13 children and three de novo loss-of-function variants (AKR1C2, PRRX1, PPM1D) in three children (one per child). We also detected rare compound heterozygous loss-of-function variants in PLCH2 and CLEC4M and rare inherited missense or loss-of-function variants in additional genes applying autosomal recessive (three genes) and X-linked recessive inheritance models (13 genes). Variants in two genes identified may implicate disruption in cell migration (TUBE1) and adhesion (TSPAN4) processes, mechanisms proposed for BE, and provide additional evidence for rare variants in the development of this defect.

Objective Our goal was to perform detailed clinical and genomic analysis of a large multigenerational Chinese family with 21 individuals showing symptoms of Familial Cortical Myoclonic Tremor with Epilepsy (FCMTE) that we have followed for over 20 years. Methods Patients were subjected to clinical evaluation, routine EEG, and structural magnetic resonance imaging. Whole exome sequencing, repeat‐primed PCR, long‐range PCR, and PacBio sequencing were performed to characterize the disease‐causing mutation in this family. Results All evaluated patients manifested adult‐onset seizures and presented with progressive myoclonic postural tremors starting after the third or fourth decade of life. Seizures typically diminished markedly in frequency with implementation of antiseizure medications but did not completely cease. The electroencephalogram of affected individuals showed generalized or multifocal spikes and slow wave complexes. An expansion of TTTTA motifs with addition of TTTCA motifs in intron 4 of SAMD12 was identified to segregate with the disease phenotype in this family. Furthermore, we found that the mutant allele is unstable and can undergo both contraction and expansion by changes in the number of repeat motifs each time it is passed to the next generation. The size of mutant allele varied from 5 to 5.5 kb with 549‐603 copies of TTTTA and 287‐343 copies of TTTCA repeat motifs in this family. Significance Our study provides a detailed description of clinical progression of FCMTE symptoms and its management with antiseizure medications. Our method of repeat analysis by PacBio sequencing of long‐range PCR products does not require high‐quality DNA and hence can be easily applied to other families to elucidate any correlation between the repeat size and phenotypic variables, such as, age of onset, and severity of symptoms.

Significance HLA alleles have previously been implicated with scleroderma risk, but, in this study, using a European American ancestral cohort and a newly recruited large cohort of African Americans, we comprehensively define the HLA alleles and amino acid residues associated with scleroderma. Scleroderma is characterized by mutually exclusive and specific autoantibodies. We demonstrated ancestry-predominant HLA alleles that were much more strongly associated with autoantibody subsets of scleroderma than with the overall risk of SSc. We bioinformatically predicted immunodominant peptides of self-antigens and demonstrated homology of these peptides with viral protein sequences from Mimiviridae and Phycodnaviridae families. Our findings suggest the hypothesis that scleroderma-specific autoantibodies may arise through molecular mimicry, driven by the interaction of specific viral antigens with corresponding HLA α/β heterodimers., Systemic sclerosis (SSc) is a clinically heterogeneous autoimmune disease characterized by mutually exclusive autoantibodies directed against distinct nuclear antigens. We examined HLA associations in SSc and its autoantibody subsets in a large, newly recruited African American (AA) cohort and among European Americans (EA). In the AA population, the African ancestry-predominant HLA-DRB1*08:04 and HLA-DRB1*11:02 alleles were associated with overall SSc risk, and the HLA-DRB1*08:04 allele was strongly associated with the severe antifibrillarin (AFA) antibody subset of SSc (odds ratio = 7.4). These African ancestry-predominant alleles may help explain the increased frequency and severity of SSc among the AA population. In the EA population, the HLA-DPB1*13:01 and HLA-DRB1*07:01 alleles were more strongly associated with antitopoisomerase (ATA) and anticentromere antibody-positive subsets of SSc, respectively, than with overall SSc risk, emphasizing the importance of HLA in defining autoantibody subtypes. The association of the HLA-DPB1*13:01 allele with the ATA+ subset of SSc in both AA and EA patients demonstrated a transancestry effect. A direct correlation between SSc prevalence and HLA-DPB1*13:01 allele frequency in multiple populations was observed (r = 0.98, P = 3 × 10−6). Conditional analysis in the autoantibody subsets of SSc revealed several associated amino acid residues, mostly in the peptide-binding groove of the class II HLA molecules. Using HLA α/β allelic heterodimers, we bioinformatically predicted immunodominant peptides of topoisomerase 1, fibrillarin, and centromere protein A and discovered that they are homologous to viral protein sequences from the Mimiviridae and Phycodnaviridae families. Taken together, these data suggest a possible link between HLA alleles, autoantibodies, and environmental triggers in the pathogenesis of SSc.

Receptor Interacting Protein Kinase 1 (RIPK1) is a key regulator of innate immune signalling pathways. To ensure an optimal inflammatory response, RIPK1 is post-translationally regulated by well characterised ubiquitylation and phosphorylation events, as well as caspase-8 mediated cleavage1–7. The physiological relevance of this cleavage remains unclear, though it is believed to inhibit activation of RIPK3 and necroptosis8. Here we show that heterozygous missense mutations p.D324N, p.D324H and p.D324Y prevent caspase cleavage of RIPK1 in humans and result in early-onset periodic fever episodes and severe intermittent lymphadenopathy, a condition we designate ‘Cleavage-resistant RIPK1-Induced Autoinflammatory’ (CRIA) syndrome. To define the mechanism for this disease we generated a cleavage-resistant Ripk1D325A mutant mouse strain. While Ripk1-/- mice die postnatally from systemic inflammation, Ripk1D325A/D325A mice died during embryogenesis. Embryonic lethality was completely prevented by combined loss of Casp8 and Ripk3 but not by loss of Ripk3 or Mlkl alone. Loss of RIPK1 kinase activity also prevented Ripk1D325A/D325A embryonic lethality, however the mice died before weaning from multi organ inflammation in a RIPK3 dependent manner. Consistently, Ripk1D325A/D325A and Ripk1D325A/+ cells were hypersensitive to RIPK3 dependent TNF-induced apoptosis and necroptosis. Heterozygous Ripk1D325A/+ mice were viable and grossly normal, but were hyper-responsive to inflammatory stimuli in vivo. Our results demonstrate the importance of caspase-mediated RIPK1 cleavage during embryonic development and show that caspase cleavage of RIPK1 not only inhibits necroptosis but maintains inflammatory homeostasis throughout life.

In 2001, Celera Genomics and the International Human Genome Sequencing Consortium published their initial drafts of the human genome, which revolutionized the field of genomics. While these drafts and the updates that followed effectively covered the euchromatic fraction of the genome, the heterochromatin and many other complex regions were left unfinished or erroneous. Addressing this remaining 8% of the genome, the Telomere-to-Telomere (T2T) Consortium has finished the first truly complete 3.055 billion base pair (bp) sequence of a human genome, representing the largest improvement to the human reference genome since its initial release. The new T2T-CHM13 reference includes gapless assemblies for all 22 autosomes plus Chromosome X, corrects numerous errors, and introduces nearly 200 million bp of novel sequence containing 2,226 paralogous gene copies, 115 of which are predicted to be protein coding. The newly completed regions include all centromeric satellite arrays and the short arms of all five acrocentric chromosomes, unlocking these complex regions of the genome to variational and functional studies for the first time.

Endometrioid endometrial carcinomas (EECs) are the most common histological subtype of uterine cancer. Late-stage disease is an adverse prognosticator for EEC. The purpose of this study was to analyze EEC exome mutation data to identify late-stage-specific statistically significantly mutated genes (SMGs), which represent candidate driver genes potentially associated with disease progression. We exome sequenced 15 late-stage (stage III or IV) non-ultramutated EECs and paired non-tumor DNAs; somatic variants were called using Strelka, Shimmer, Somatic Sniper and MuTect. Additionally, somatic mutation calls were extracted from The Cancer Genome Atlas (TCGA) data for 66 late-stage and 270 early-stage (stage I or II) non-ultramutated EECs. MutSigCV (v1.4) was used to annotate SMGs in the two late-stage cohorts and to derive p-values for all mutated genes in the early-stage cohort. To test whether late-stage SMGs are statistically significantly mutated in early-stage tumors, q-values for late-stage SMGs were re-calculated from the MutSigCV (v1.4) early-stage p-values, adjusting for the number of late-stage SMGs tested. We identified 14 SMGs in the combined late-stage EEC cohorts. When the 14 late-stage SMGs were examined in the TCGA early-stage data, only KLF3 and PAX6 failed to reach significance as early-stage SMGs, despite the inclusion of enough early-stage cases to ensure adequate statistical power. Within TCGA, nonsynonymous mutations in KLF3 and PAX6 were, respectively, exclusive or nearly exclusive to the microsatellite instability (MSI)-hypermutated molecular subgroup and were dominated by insertions-deletions at homopolymer tracts. In conclusion, our findings are hypothesis-generating and suggest that KLF3 and PAX6, which encode transcription factors, are MSI target genes and late-stage-specific SMGs in EEC.

(1) Background: Biomarkers might play a significant role in predicting the clinical outcomes of patients with ovarian cancer. By analyzing lipid metabolism genes, future perspectives may be uncovered, (2) Methods: RNA-seq data for serous ovarian cancer were downloaded from The Cancer Genome Atlas and Gene Expression Omnibus databases. The non-negative matrix factorization package in programming language R was used to classify molecular subtypes of lipid metabolism genes and the limma package in R was performed for functional enrichment analysis. Through lasso regression, we constructed a multi-gene prognosis model, (3) Results: Two molecular subtypes were obtained and an 11-gene signature was constructed (PI3, RGS, ADORA3, CH25H, CCDC80, PTGER3, MATK, KLRB1, CCL19, CXCL9 and CXCL10). Our prognostic model shows a good independent prognostic ability in ovarian cancer. In a nomogram, the predictive efficiency was notably superior to that of traditional clinical features. Related to known models in ovarian cancer with a comparable amount of genes, ours has the highest concordance index, (4) Conclusions: We propose an 11-gene signature prognosis prediction model based on lipid metabolism genes in serous ovarian cancer.

Identifying the molecular mechanisms by which genome-wide association study (GWAS) loci influence traits remains challenging. Chromatin accessibility quantitative trait loci (caQTLs) help identify GWAS loci that may alter GWAS traits by modulating chromatin structure, but caQTLs have been identified in a limited set of human tissues. Here we mapped caQTLs in human liver tissue in 20 liver samples and identified 3,123 caQTLs. The caQTL variants are enriched in liver tissue promoter and enhancer states and frequently disrupt binding motifs of transcription factors expressed in liver. We predicted target genes for 861 caQTL peaks using proximity, chromatin interactions, correlation with promoter accessibility or gene expression, and colocalization with expression QTLs. Using GWAS signals for 19 liver function and/or cardiometabolic traits, we identified 110 colocalized caQTLs and GWAS signals, 56 of which contained a predicted caPeak target gene. At the LITAF LDL-cholesterol GWAS locus, we validated that a caQTL variant showed allelic differences in protein binding and transcriptional activity. These caQTLs contribute to the epigenomic characterization of human liver and help identify molecular mechanisms and genes at GWAS loci.

After two decades of improvements, the current human reference genome (GRCh38) is the most accurate and complete vertebrate genome ever produced. However, no single chromosome has been finished end to end, and hundreds of unresolved gaps persist1,2. Here we present a human genome assembly that surpasses the continuity of GRCh382, along with a gapless, telomere-to-telomere assembly of a human chromosome. This was enabled by high-coverage, ultra-long-read nanopore sequencing of the complete hydatidiform mole CHM13 genome, combined with complementary technologies for quality improvement and validation. Focusing our efforts on the human X chromosome3, we reconstructed the centromeric satellite DNA array (approximately 3.1 Mb) and closed the 29 remaining gaps in the current reference, including new sequences from the human pseudoautosomal regions and from cancer-testis ampliconic gene families (CT-X and GAGE). These sequences will be integrated into future human reference genome releases. In addition, the complete chromosome X, combined with the ultra-long nanopore data, allowed us to map methylation patterns across complex tandem repeats and satellite arrays. Our results demonstrate that finishing the entire human genome is now within reach, and the data presented here will facilitate ongoing efforts to complete the other human chromosomes., High-coverage, ultra-long-read nanopore sequencing is used to create a new human genome assembly that improves on the coverage and accuracy of the current reference (GRCh38) and includes the gap-free, telomere-to-telomere sequence of the X chromosome.

Methotrexate, a competitive inhibitor of dihydrofolate reductase, is an integral part of treatment for childhood leukemia. High-dose methotrexate (HDMTX) significantly improves survival in NCI-Rome...

Hypertension occurs at a higher rate in African Americans than in European Americans. Based on the assumption that causal variants are more frequently found on DNA segments inherited from the ancestral population with higher disease risk, we employed admixture mapping to identify genetic loci with excess local African ancestry associated with blood pressure. Chromosomal regions 1q21.2-21.3, 4p15.1, 19q12 and 20p13 were significantly associated with diastolic blood pressure (β = 5.28, -7.94, -6.82 and 5.89, P-value = 6.39E-04, 2.07E-04, 6.56E-05 and 5.04E-04, respectively); 1q21.2-21.3 and 19q12 were also significantly associated with mean arterial pressure (β = 5.86 and -6.40, P-value = 5.32E-04 and 6.37E-04, respectively). We further selected SNPs that had large allele frequency differences within these regions and tested their association with blood pressure. SNP rs4815428 was significantly associated with diastolic blood pressure after Bonferroni correction (β = -2.42, P-value = 9.57E-04), and it partially explained the admixture mapping signal at 20p13. SNPs rs771205 (β = -1.99, P-value = 3.37E-03), rs3126067, rs2184953 and rs58001094 (the latter three exhibit strong linkage disequilibrium, β = -2.3, P-value = 1.4E-03) were identified to be significantly associated with mean arterial pressure, and together they fully explained the admixture signal at 1q21.2-21.3. Although no SNP at 4p15.1 showed large ancestral allele frequency differences in our dataset, we detected association at low-frequency African-specific variants that mapped predominantly to the gene PCDH7, which is most highly expressed in aorta. Our results suggest that these regions may harbor genetic variants that contribute to the different prevalence of hypertension.

We show that BDNF regulates the timing of neurodevelopment via a novel mechanism of extranuclear sequestration of NFATc4 in Golgi. This leads to accelerated derepression of an NFI temporal occupancy gene program in cerebellar granule cells that includes Bdnf itself, revealing an autoregulatory loop within the program driven by BDNF and NFATc4., How intrinsic and extrinsic signals are coordinated to regulate synaptic maturation and its timing is an important question for neurodevelopment and its disorders. We investigated the influence of the neurotrophin BDNF on the developmental timing of a dendrite/synapse-related gene program controlled by nuclear factor I (NFI) in maturing cerebellar granule neurons (CGNs). BDNF accelerated the onset of NFI-regulated late-gene expression and NFI temporal occupancy in CGN cultures in a MEK5/ERK5-dependent manner. BDNF and NFI occupancy were mutually regulating, with BDNF enhancing the temporal binding of NFI to the Bdnf4 promoter itself. Moreover, BDNF induced phosphorylation and accelerated the departure of the trans-repressor NFATc4 from NFI late-gene promoters, including Bdnf4, which is permissive for NFI binding. BDNF dismissal of NFATc4 from late genes was linked to MEK5/ERK5-dependent sequestration of NFATc4 in the cis–Golgi, an event mirrored in CGNs developing in vivo. These studies reveal an expanded autoregulatory gene network for NFI temporal occupancy involving BDNF and NFATc4 extranuclear sequestration. Based on these and earlier findings, NFATc4 integrates intrinsic developmental signaling from membrane potential/calcineurin and autocrine/paracrine BDNF/TrkB to control initiation of NFI occupancy in maturing CGNs. We also identify a local Bdnf/Etv1 gene circuit within the larger NFI autoregulatory network.

Background Dermatophagoides pteronyssinus (DP) and Dermatophagoides farinae (DF) are highly similar disease-associated mites with frequently overlapping geographic distributions. A draft genome of DP was assembled to identify candidate allergens in DP homologous to those in DF, investigate allergen isoforms, and facilitate comparisons with related Acari. Methods PacBio and Illumina whole genome sequencing was performed on DP. Assembly and reconstruction of the genomes were optimized for isoform identification in a heterogeneous population. Bioinformatic analyses of Acari genomes were peformed. Results The predicted size of the DP nuclear genome is 52.5 Mb. A predicted protein set of 19,368 proteins was identified, including all 19 currently recognized allergens from this species. Orthologs for 12 allergens established for DF were found. The population of DP mites showed a high level of heterozygosity that allowed the identification of 43 new isoforms for both established and candidate allergens in DP, including a new isoform for the major allergen Der p 23. Reanalyzing the previous DF data assuming heterozygosity, 14 new allergen isoforms could be indentified. Some new isoforms were observed in both species suggesting that these isoforms pre-dated speciation. The high quality of both genomes allowed an examination of synteny which showed many allergen orthologs are physically clustered but with species specific exon/intron structures. Comparative genomic analyses with other Acariformes mites showed that most of the allergen homologs are widely conserved within this Superorder. Conclusions Candidate allergens in DP were identified to facilitate future serological studies. While DP and DF are highly similar genetically, species-specific allergen isoforms exist to facilitate molecular differentiation.

BACKGROUND Uterine carcinosarcomas (UCSs) are a rare but clinically aggressive form of cancer. They are biphasic tumors consisting of both epithelial and sarcomatous components. The majority of uterine carcinosarcomas are clonal, with the carcinomatous cells undergoing metaplasia to give rise to the sarcomatous component. The objective of the current study was to identify novel somatically mutated genes in UCSs. METHODS We whole exome sequenced paired tumor and nontumor DNAs from 14 UCSs and orthogonally validated 464 somatic variants using Sanger sequencing. Fifteen genes that were somatically mutated in at least 2 tumor exomes were Sanger sequenced in another 39 primary UCSs. RESULTS Overall, among 53 UCSs in the current study, the most frequently mutated of these 15 genes were tumor protein p53 (TP53) (75.5%), phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA) (34.0%), protein phosphatase 2, regulatory subunit A, alpha (PPP2R1A) (18.9%), F-box and WD repeat domain containing 7 (FBXW7) (18.9%), chromodomain helicase DNA binding protein 4 (CHD4) (17.0%), and forkhead box A2 (FOXA2) (15.1%). FOXA2 has not previously been implicated in UCSs and was predominated by frameshift and nonsense mutations. One UCS with a FOXA2 frameshift mutation expressed truncated FOXA2 protein by immunoblotting. Sequencing of FOXA2 in 160 primary endometrial carcinomas revealed somatic mutations in 5.7% of serous, 22.7% of clear cell, 9% of endometrioid, and 11.1% of mixed endometrial carcinomas, the majority of which were frameshift mutations. CONCLUSIONS Collectively, the findings of the current study provide compelling genetic evidence that FOXA2 is a pathogenic driver gene in the etiology of primary uterine cancers, including UCSs. Cancer 2017. © 2017 American Cancer Society.

BACKGROUND The molecular pathogenesis of clear cell endometrial cancer (CCEC), a tumor type with a relatively unfavorable prognosis, is not well defined. We searched exome-wide for novel somatically mutated genes in CCEC and assessed the mutational spectrum of known and candidate driver genes in a large cohort of cases. METHODS We conducted whole exome sequencing of paired tumor-normal DNAs from 16 cases of CCEC (12 CCECs and the CCEC components of 4 mixed histology tumors). Twenty-two genes-of-interest were Sanger-sequenced from another 47 cases of CCEC. Microsatellite instability (MSI) and microsatellite stability (MSS) were determined by genotyping 5 mononucleotide repeats. Results Two tumor exomes had relatively high mutational loads and MSI. The other 14 tumor exomes were MSS and had 236 validated nonsynonymous or splice junction somatic mutations among 222 protein-encoding genes. Among the 63 cases of CCEC in this study, we identified frequent somatic mutations in TP53 (39.7%), PIK3CA (23.8%), PIK3R1 (15.9%), ARID1A (15.9%), PPP2R1A (15.9%), SPOP (14.3%), and TAF1 (9.5%), as well as MSI (11.3%). Five of 8 mutations in TAF1, a gene with no known role in CCEC, localized to the putative histone acetyltransferase domain and included 2 recurrently mutated residues. Based on patterns of MSI and mutations in 7 genes, CCEC subsets molecularly resembled serous endometrial cancer (SEC) or endometrioid endometrial cancer (EEC). CONCLUSION Our findings demonstrate molecular similarities between CCEC and SEC and EEC and implicate TAF1 as a novel candidate CCEC driver gene. Cancer 2017. © 2017 American Cancer Society.

Purpose Hermansky-Pudlak syndrome (HPS) is a rare inherited disorder with ten reported genetic types; each type has defects in subunits of either Adaptor Protein-3 complex or Biogenesis of Lysosome-related Organelles Complex (BLOC)-1, -2, or -3. Very few patients with BLOC-1 deficiency (HPS-7, -8, and -9 types) have been diagnosed. We report results of comprehensive clinical testing and molecular analyses of primary fibroblasts from a new case of HPS-7. Results A 6-year old Paraguayan male presented with hypopigmentation, ocular albinism, nystagmus, reduced visual acuity, and easy bruising. He also experienced delayed motor and language development as a very young child; head and chest trauma resulted in intracranial hemorrhage with subsequent right hemiparesis and lung scarring. There was no clinical evidence of immunodeficiency or colitis. Whole mount transmission electron microscopy revealed absent platelet delta granules; platelet aggregation testing was abnormal. Exome sequencing revealed a homozygous nonsense mutation in the Dystrobrevin binding protein 1 (DTNBP1) gene [ NM_032122.4 : c.307C > T; p.Gln103*], previously reported in a Portuguese adult. The gene encodes the dysbindin subunit of BLOC-1. Dysbindin protein expression was negligible in our patient's dermal fibroblasts, while his DTNBP1 mRNA level was similar to that of a normal control. Conclusions Comprehensive clinical evaluation of the first pediatric case reported with HPS-7 reveals oculocutaneous albinism and platelet storage pool deficiency; his phenotype is consistent with findings in other patients with BLOC-1 disorders. This patient's markedly reduced Dysbindin protein expression in HPS-7 resulted from a mechanism other than nonsense mediated decay.

Joubert syndrome and related disorders (JSRD) are a heterogeneous group of ciliopathies defined based on the mid-hindbrain abnormalities that result in the characteristic “molar tooth sign” on brain imaging. The core clinical findings of JSRD are hypotonia, developmental delay, abnormal eye movements and breathing abnormalities. To date, more than 30 JSRD genes that encode proteins important for structure and/or function of cilia have been identified. Here, we present 2 siblings with Joubert syndrome associated with growth hormone deficiency. Whole exome sequencing of the family identified compound heterozygous mutations in KIAA0753, i.e., a missense mutation (p.Arg257Gly) and an intronic mutation (c.2359-1G>C). The intronic mutation alters normal splicing by activating a cryptic acceptor splice site in exon 16. The novel acceptor site skips nine nucleotides, deleting three amino acids from the protein coding frame. KIAA0753 (OFIP) is a centrosome and pericentriolar satellite protein, previously not known to cause Joubert syndrome. We present comprehensive clinical descriptions of the Joubert syndrome patients as well as the cellular phenotype of defective ciliogenesis in the patients’ fibroblasts.

Neurofibromatosis type 1 (NF1) is a common tumor-predisposition disorder due to germline mutations in the tumor suppressor gene NF1 A virtually pathognomonic finding of NF1 is the plexiform neurofibroma (PN), a benign, likely congenital tumor that arises from bi-allelic inactivation of NF1 PN can undergo transformation to a malignant peripheral nerve sheath tumor, an aggressive soft-tissue sarcoma. To better understand the non-NF1 genetic contributions to PN pathogenesis, we performed whole-exome sequencing, RNASeq profiling and genome-wide copy-number determination for 23 low-passage Schwann cell cultures established from surgical PN material with matching germline DNA. All resected tumors were derived from routine debulking surgeries. None of the tumors were considered at risk for malignant transformation at the time;for example, there was no pain or rapid growth. Deep (~500X) NF1 exon sequencing was also conducted on tumor DNA. Non-NF1 somatic mutation verification was performed using the Ampliseq/IonTorrent platform. We identified 100% of the germline NF1 mutations and found somatic NF1 inactivation in 74% of the PN. One individual with three PNs had different NF1 somatic mutations in each tumor. The median number of somatic mutations per sample, including NF1, was one (range 0–8). NF1 was the only gene that was recurrently somatically inactivated in multiple tumors. Gene Set Enrichment Analysis of transcriptome-wide tumor RNA sequencing identified five significant (FDR < 0.01) and seven trending (0.01 ≤ FDR < 0.02) gene sets related to DNA replication, telomere maintenance and elongation, cell cycle progression, signal transduction and cell proliferation. We found no recurrent non-NF1 locus copy-number variation in PN. This is the first multi-sample whole-exome and whole-transcriptome sequencing study of NF1-associated PN. Taken together with concurrent copy-number data, our comprehensive genetic analysis reveals the primacy of NF1 loss as the driver of PN tumorigenesis.

The analysis of pineal cell biology has undergone remarkable development as techniques have become available which allow for sequencing of entire transcriptomes and, most recently, the sequencing of the transcriptome of individual cells. Identification of at least nine distinct cell types in the rat pineal gland has been made possible, allowing identification of the precise cells of origin and expression of transcripts for the first time. Here the history and current state of knowledge generated by these transcriptomic efforts is reviewed, with emphasis on the insights suggested by the findings.

Background Neurofibromatosis type 1 (NF1) is a tumor-predisposition disorder caused by germline mutations in NF1. NF1 patients have an 8–16% lifetime risk of developing a malignant peripheral nerve sheath tumor (MPNST), a highly aggressive soft-tissue sarcoma, often arising from preexisting benign plexiform neurofibromas (PNs) and atypical neurofibromas (ANFs). ANFs are distinct from both PN and MPNST, representing an intermediate step in malignant transformation. Methods In the first comprehensive genomic analysis of ANF originating from multiple patients, we performed tumor/normal whole-exome sequencing (WES) of 16 ANFs. In addition, we conducted WES of 3 MPNSTs, copy-number meta-analysis of 26 ANFs and 28 MPNSTs, and whole transcriptome sequencing analysis of 5 ANFs and 5 MPNSTs. Results We identified a low number of mutations (median 1, range 0–5) in the exomes of ANFs (only NF1 somatic mutations were recurrent), and frequent deletions of CDKN2A/B (69%) and SMARCA2 (42%). We determined that polycomb repressor complex 2 (PRC2) genes EED and SUZ12 were frequently mutated, deleted, or downregulated in MPNSTs but not in ANFs. Our pilot gene expression study revealed upregulated NRAS, MDM2, CCND1/2/3, and CDK4/6 in ANFs and MPNSTs, and overexpression of EZH2 in MPNSTs only. Conclusions The PN-ANF transition is primarily driven by the deletion of CDKN2A/B. Further progression from ANF to MPNST likely involves broad chromosomal rearrangements and frequent inactivation of the PRC2 genes, loss of the DNA repair genes, and copy-number increase of signal transduction and cell-cycle and pluripotency self-renewal genes.

For over a thousand years, the common goldfish (Carassius auratus) was raised throughout Asia for food and as an ornamental pet. As a very close relative of the common carp (Cyprinus carpio), goldfish share the recent genome duplication that occurred approximately 14 million years ago in their common ancestor. The combination of centuries of breeding and a wide array of interesting body morphologies provides an exciting opportunity to link genotype to phenotype and to understand the dynamics of genome evolution and speciation. We generated a high-quality draft sequence and gene annotations of a “Wakin” goldfish using 71X PacBio long reads. The two subgenomes in goldfish retained extensive synteny and collinearity between goldfish and zebrafish. However, genes were lost quickly after the carp whole-genome duplication, and the expression of 30% of the retained duplicated gene diverged substantially across seven tissues sampled. Loss of sequence identity and/or exons determined the divergence of the expression levels across all tissues, while loss of conserved noncoding elements determined expression variance between different tissues. This assembly provides an important resource for comparative genomics and understanding the causes of goldfish variants.

Joubert and Meckel-Gruber syndromes (JS and MGS) are ciliopathies with overlapping features. JS patients manifest the “molar tooth sign” on brain imaging and variable eye, kidney and liver disease. MGS presents with polycystic kidneys, occipital encephalocele and polydactyly; it is typically perinatally fatal. Both syndromes are genetically heterogeneous; some genes cause either syndrome. Here, we report two brothers married to unrelated women. The first brother had three daughters with JS and a son with polycystic kidneys who died at birth. The second brother's wife had a fetal demise due to MGS. Whole exome sequencing identified TMEM231 NM_001077416.2: c.784G>A; p.(Asp262Asn) in all children and the wife of the first brother; the second brother's wife had a c.406T>G;p.(Trp136Gly) change. In-depth analysis uncovered a rare gene conversion event in TMEM231, leading to loss of exon 4, in all the affected children of first brother. We believe that the combination of this gene conversion with different missense mutations led to a spectrum of phenotypes that span JS and MGS. This article is protected by copyright. All rights reserved

Antibodies capable of neutralizing HIV-1 often target variable regions 1 and 2 (V1V2) of the HIV-1 envelope, but the mechanism of their elicitation has been unclear. Here we define the developmental pathway by which such antibodies are generated and acquire the requisite molecular characteristics for neutralization. Twelve somatically related neutralizing antibodies (CAP256-VRC26.01–12) were isolated from donor CAP256 (from the Centre for the AIDS Programme of Research in South Africa (CAPRISA)); each antibody contained the protruding tyrosine-sulphated, anionic antigen-binding loop (complementarity-determining region (CDR) H3) characteristic of this category of antibodies. Their unmutated ancestor emerged between weeks 30–38 post-infection with a 35-residue CDR H3, and neutralized the virus that superinfected this individual 15 weeks after initial infection. Improved neutralization breadth and potency occurred by week 59 with modest affinity maturation, and was preceded by extensive diversification of the virus population. HIV-1 V1V2-directed neutralizing antibodies can thus develop relatively rapidly through initial selection of B cells with a long CDR H3, and limited subsequent somatic hypermutation. These data provide important insights relevant to HIV-1 vaccine development. A longitudinal study of an individual patient developing neutralizing antibodies against HIV-1 (targeting the V1V2 region of gp120) reveals how such neutralizing antibodies develop and evolve over time, providing important insights relevant to vaccine development. A better understanding of how HIV-1-neutralizing antibodies are generated could be a useful contribution to the design of improved AIDS vaccines. John Mascola and colleagues have now elucidated the immunological pathway of an important category of HIV-1-neutralizing antibody — those that target the variable V1V2 region of the viral envelope. These antibodies are more frequently elicited than CD4-binding site antibodies in the early stages of HIV infection and feature modest affinity maturation, a process that favours mutations in antibody variable domains that enhance antigen binding.

Antibodies capable of neutralizing divergent influenza A viruses could form the basis of a universal vaccine. Here, from subjects enrolled in an H5N1 DNA/MIV-prime-boost influenza vaccine trial, we sorted hemagglutinin cross-reactive memory B cells and identified three antibody classes, each capable of neutralizing diverse subtypes of group 1 and group 2 influenza A viruses. Co-crystal structures with hemagglutinin revealed that each class utilized characteristic germline genes and convergent sequence motifs to recognize overlapping epitopes in the hemagglutinin stem. All six analyzed subjects had sequences from at least one multidonor class, and—in half the subjects—multidonor-class sequences were recovered from >40% of cross-reactive B cells. By contrast, these multidonor-class sequences were rare in published antibody datasets. Vaccination with a divergent hemagglutinin can thus increase the frequency of B cells encoding broad influenza A-neutralizing antibodies. We propose the sequence signature-quantified prevalence of these B cells as a metric to guide universal influenza A immunization strategies.

Temporal control of dendritogenesis is poorly understood. Mutual feedback between NFIA temporal occupancy and ETV1 drives the timing of gene expression associated with dendrite formation in maturing neurons. A sequential timing model is proposed in which ETV1 autoregulation precedes activation of downstream NFIA/ETV1 coregulated genes., Nuclear Factor One (NFI) transcription factors regulate temporal gene expression required for dendritogenesis and synaptogenesis via delayed occupancy of target promoters in developing cerebellar granule neurons (CGNs). Mechanisms that promote NFI temporal occupancy have not been previously defined. We show here that the transcription factor ETV1 directly binds to and is required for expression and NFI occupancy of a cohort of NFI-dependent genes in CGNs maturing in vivo. Expression of ETV1 is low in early postnatal cerebellum and increases with maturation, mirroring NFI temporal occupancy of coregulated target genes. Precocious expression of ETV1 in mouse CGNs accelerated onset of expression and NFI temporal occupancy of late target genes and enhanced Map2(+) neurite outgrowth. ETV1 also activated expression and NFI occupancy of the Etv1 gene itself, and this autoregulatory loop preceded ETV1 binding and activation of other coregulated target genes in vivo. These findings suggest a potential model in which ETV1 activates NFI temporal binding to a subset of late-expressed genes in a stepwise manner by initial positive feedback regulation of the Etv1 gene itself followed by activation of downstream coregulated targets as ETV1 expression increases. Sequential transcription factor autoregulation and subsequent binding to downstream promoters may provide an intrinsic developmental timer for dendrite/synapse gene expression.

BACKGROUND Next-generation sequencing (NGS) data are used for both clinical care and clinical research. DNA sequence variants identified using NGS are often returned to patients/participants as part of clinical or research protocols. The current standard of care is to validate NGS variants using Sanger sequencing, which is costly and time-consuming. METHODS We performed a large-scale, systematic evaluation of Sanger-based validation of NGS variants using data from the ClinSeq® project. We first used NGS data from 19 genes in 5 participants, comparing them to high-throughput Sanger sequencing results on the same samples, and found no discrepancies among 234 NGS variants. We then compared NGS variants in 5 genes from 684 participants against data from Sanger sequencing. RESULTS Of over 5800 NGS-derived variants, 19 were not validated by Sanger data. Using newly designed sequencing primers, Sanger sequencing confirmed 17 of the NGS variants, and the remaining 2 variants had low quality scores from exome sequencing. Overall, we measured a validation rate of 99.965% for NGS variants using Sanger sequencing, which was higher than many existing medical tests that do not necessitate orthogonal validation. CONCLUSIONS A single round of Sanger sequencing is more likely to incorrectly refute a true-positive variant from NGS than to correctly identify a false-positive variant from NGS. Validation of NGS-derived variants using Sanger sequencing has limited utility, and best practice standards should not include routine orthogonal Sanger validation of NGS variants.

Most expression quantitative trait locus (eQTL) studies to date have been performed in heterogeneous tissues as opposed to specific cell types. To better understand the cell-type–specific regulatory landscape of human melanocytes, which give rise to melanoma but account for cis-eQTL SNPs prior to linkage disequilibrium (LD) pruning and 4997 eGenes (FDR < 0.05). Melanocyte eQTLs differed considerably from those identified in the 44 GTEx tissue types, including skin. Over a third of melanocyte eGenes, including key genes in melanin synthesis pathways, were unique to melanocytes compared to those of GTEx skin tissues or TCGA melanomas. The melanocyte data set also identified trans-eQTLs, including those connecting a pigmentation-associated functional SNP with four genes, likely through cis-regulation of IRF4. Melanocyte eQTLs are enriched in cis-regulatory signatures found in melanocytes as well as in melanoma-associated variants identified through genome-wide association studies. Melanocyte eQTLs also colocalized with melanoma GWAS variants in five known loci. Finally, a transcriptome-wide association study using melanocyte eQTLs uncovered four novel susceptibility loci, where imputed expression levels of five genes (ZFP90, HEBP1, MSC, CBWD1, and RP11-383H13.1) were associated with melanoma at genome-wide significant P-values. Our data highlight the utility of lineage-specific eQTL resources for annotating GWAS findings, and present a robust database for genomic research of melanoma risk and melanocyte biology.

SummaryFor over a thousand years throughout Asia, the common goldfish (Carassius auratus) was raised for both food and as an ornamental pet. Selective breeding over more than 500 years has created a wide array of body and pigmentation variation particularly valued by ornamental fish enthusiasts. As a very close relative of the common carp (Cyprinus carpio), goldfish shares the recent genome duplication that occurred approximately 14-16 million years ago (mya) in their common ancestor. The combination of centuries of breeding and a wide array of interesting body morphologies is an exciting opportunity to link genotype to phenotype as well as understanding the dynamics of genome evolution and speciation. Here we generated a high-quality draft sequence of a “Wakin” goldfish using 71X PacBio long-reads. We identified 70,324 coding genes and more than 11,000 non-coding transcripts. We found that the two sub-genomes in goldfish retained extensive synteny and collinearity between goldfish and zebrafish. However, “ohnologous” genes were lost quickly after the carp whole-genome duplication, and the expression of 30% of the retained duplicated gene diverged significantly across seven tissues sampled. Loss of sequence identity and/or exons determined the divergence of the expression across all tissues, while loss of conserved, non-coding elements determined expression variance between different tissues. This draft assembly also provides an important resource for comparative genomics with the very commonly used zebrafish model (Danio rerio), and for understanding the underlying genetic causes of goldfish variants.

The prevalence of chronic kidney disease (CKD) is rising worldwide and 10-15% of the global population currently suffers from CKD and its complications. Given the increasing prevalence of CKD there is an urgent need to find novel treatment options. The American black bear (Ursus americanus) copes with months of lowered kidney function and metabolism during hibernation without the devastating effects on metabolism and other consequences observed in humans. In a biomimetic approach to better understand kidney adaptations and physiology in hibernating black bears, we established a high-quality genome assembly. Subsequent RNA-Seq analysis of kidneys comparing gene expression profiles in black bears entering (late fall) and emerging (early spring) from hibernation identified 169 protein-coding genes that were differentially expressed. Of these, 101 genes were downregulated and 68 genes were upregulated after hibernation. Fold changes ranged from 1.8-fold downregulation (RTN4RL2) to 2.4-fold upregulation (CISH). Most notable was the upregulation of cytokine suppression genes (SOCS2, CISH, and SERPINC1) and the lack of increased expression of cytokines and genes involved in inflammation. The identification of these differences in gene expression in the black bear kidney may provide new insights in the prevention and treatment of CKD.

Here, we applied targeted capture to examine 153 genes representative of all the major vertebrate developmental pathways among 333 probands to rank their relative significance as causes for holoprosencephaly (HPE). We now show that comparisons of variant transmission versus nontransmission among 136 HPE Trios indicates some reported genes now lack confirmation, while novel genes are implicated. Furthermore, we demonstrate that variation of modest intrinsic effect can synergize with these driver mutations as gene modifiers.

Background Mutations of TCF4, which encodes a basic helix-loop-helix transcription factor, cause Pitt-Hopkins syndrome (PTHS) via multiple genetic mechanisms. TCF4 is a complex locus expressing multiple transcripts by alternative splicing and use of multiple promoters. To address the relationship between mutation of these transcripts and phenotype, we report a three-generation family segregating mild intellectual disability with a chromosomal translocation disrupting TCF4. Results Using whole genome sequencing, we detected a complex unbalanced karyotype disrupting TCF4 (46,XY,del(14)(q23.3q23.3)del(18)(q21.2q21.2)del(18)(q21.2q21.2)inv(18)(q21.2q21.2)t(14;18)(q23.3;q21.2)(14pter®14q23.3::18q21.2®18q21.2::18q21.1®18qter;18pter®18q21.2::14q23.3®14qter). Subsequent transcriptome sequencing, qRT-PCR and nCounter analyses revealed that cultured skin fibroblasts and peripheral blood had normal expression of genes along chromosomes 14 or 18 and no marked changes in expression of genes other than TCF4. Affected individuals had 12–33 fold higher mRNA levels of TCF4 than did unaffected controls or individuals with PTHS. Although the derivative chromosome generated a PLEKHG3-TCF4 fusion transcript, the increased levels of TCF4 mRNA arose from transcript variants originating distal to the translocation breakpoint, not from the fusion transcript. Conclusions Although validation in additional patients is required, our findings suggest that the dysmorphic features and severe intellectual disability characteristic of PTHS are partially rescued by overexpression of those short TCF4 transcripts encoding a nuclear localization signal, a transcription activation domain, and the basic helix-loop-helix domain. Electronic supplementary material The online version of this article (doi:10.1186/s13023-016-0439-6) contains supplementary material, which is available to authorized users.

Broadly neutralizing antibodies (bNAbs) against HIV-1 Env V1V2 arise in multiple donors. However, atomic-level interactions had previously been determined only with antibodies from a single donor, thus making commonalities in recognition uncertain. Here we report the cocrystal structure of V1V2 with antibody CH03 from a second donor and model Env interactions of antibody CAP256-VRC26 from a third donor. These V1V2-directed bNAbs used strand-strand interactions between a protruding antibody loop and a V1V2 strand but differed in their N-glycan recognition. Ontogeny analysis indicated that protruding loops develop early, and glycan interactions mature over time. Altogether, the multidonor information suggested that V1V2-directed bNAbs form an 'extended class', for which we engineered ontogeny-specific antigens: Env trimers with chimeric V1V2s that interacted with inferred ancestor and intermediate antibodies. The ontogeny-based design of vaccine antigens described here may provide a general means for eliciting antibodies of a desired class.

Somatic mutational landscape of AML with inv(16) or t(8;21) identifies patterns of clonal evolution in relapse leukemia

Next-generation sequencing provides the opportunity to practice predictive medicine based on identified variants. Putative loss-of-function (pLOF) variants are common in genomes and understanding their contribution to disease is critical for predictive medicine. To this end, we characterized the consequences of pLOF variants in an exome cohort by iterative phenotyping. Exome data were generated on 951 participants from the ClinSeq cohort and filtered for pLOF variants in genes likely to cause a phenotype in heterozygotes. 103 of 951 exomes had such a pLOF variant and 79 participants were evaluated. Of those 79, 34 had findings or family histories that could be attributed to the variant (28 variants in 18 genes), 2 had indeterminate findings (2 variants in 2 genes), and 43 had no findings or a negative family history for the trait (34 variants in 28 genes). The presence of a phenotype was correlated with two mutation attributes: prior report of pathogenicity for the variant (p = 0.0001) and prior report of other mutations in the same exon (p = 0.0001). We conclude that 1/30 unselected individuals harbor a pLOF mutation associated with a phenotype either in themselves or their family. This is more common than has been assumed and has implications for the setting of prior probabilities of affection status for predictive medicine.

// Bo Yu 1 , Valya R. Russanova 2 , Silvia Gravina 3 , Stephen Hartley 4 , James C. Mullikin 4, 5 , Alice Ignezweski 6 , James Graham 6 , James H. Segars 7 , Alan H. DeCherney 7 , Bruce H. Howard 2 1 Department of Obstetrics and Gynecology & Women's Health, Albert Einstein College of Medicine, Bronx, New York 10461, USA 2 Program in Genomics of Differentiation, Eunice Kennedy Shriver National Institute for Child Health and Human Development, National Institutes of Health, Bethesda, Maryland 20892, USA 3 Department of Genetics, Albert Einstein College of Medicine, Bronx, New York 10461, USA 4 Comparative Genomics Unit, Genome Technology Branch, National Human Genome Research Institute, National Institutes of Health, Bethesda, Maryland 20892, USA 5 NIH Intramural Sequencing Center, National Human Genome Research Institute, National Institutes of Health, Bethesda, Maryland 20892, USA 6 Shady Grove Fertility Reproductive Science Center, Rockville, Maryland 20850, USA 7 Program in Reproductive and Adult Endocrinology, Eunice Kennedy Shriver National Institute for Child Health and Human Development, National Institutes of Health, Bethesda, Maryland 20892, USA Correspondence to: Bo Yu, e-mail: boyumich2@gmail.com Keywords: DNA methylation, transcription end site, fertility, ovarian granulosa cell, transcriptome Received: November 14, 2014 Accepted: December 08, 2014 Published: February 17, 2015 ABSTRACT Diminished ovarian function occurs early and is a primary cause for age-related decline in female fertility; however, its underlying mechanism remains unclear. This study investigated the roles that genome and epigenome structure play in age-related changes in gene expression and ovarian function, using human ovarian granulosa cells as an experimental system. DNA methylomes were compared between two groups of women with distinct age-related differences in ovarian functions, using both Methylated DNA Capture followed by Next Generation Sequencing (MethylCap-seq) and Reduced Representation Bisulfite Sequencing (RRBS); their transcriptomes were investigated using mRNA-seq. Significant, non-random changes in transcriptome and DNA methylome features are observed in human ovarian granulosa cells as women age and their ovarian functions deteriorate. The strongest correlations between methylation and the age-related changes in gene expression are not confined to the promoter region; rather, high densities of hypomethylated CpG-rich regions spanning the gene body are preferentially associated with gene down-regulation. This association is further enhanced where CpG regions are localized near the 3'-end of the gene. Such features characterize several genes crucial in age-related decline in ovarian function, most notably the AMH (Anti-Mullerian Hormone) gene. The genome-wide correlation between the density of hypomethylated intragenic and 3'-end regions and gene expression suggests previously unexplored mechanisms linking epigenome structure to age-related physiology and pathology.

A single monoclonal protein typically characterizes monoclonal gammopathies, but a small proportion may have more than one M protein identifiable. In the setting of symptomatic multiple myeloma (MM), the development of a new monoclonal protein following therapy is associated with better outcomes. As for the precursor conditions, monoclonal gammopathy undetermined significance (MGUS) and smoldering multiple myeloma (SMM), there is limited information on the impact of a second monoclonal protein on the disease course, including progression and response to treatment. The outcomes of patients with MGUS and SMM with more than one monoclonal protein, after identifying 539 patients with biclonal proteins on electrophoresis and/or immunofixation, were reported. About 22 of 393 patients with MGUS/biclonal gammopathy of undetermined significance (BGUS) progressed to SMM (6), MM (11), AL (3), or WM (2), and 5 of 16 patients with biclonal SMM progressed to MM. The rate of progression for BGUS was approximately 1% per year, which is similar to MGUS with one monoclonal protein. The median estimated time of progression of biclonal SMM was 2.6 years; similar to monoclonal SMM. For patients with biclonal MM, both M spikes responded to treatment and, upon relapse, the original dominant M protein remained dominant as the disease progressed. In conclusion, the presence of a second monoclonal protein does not appear to affect the progression of precursor states and suggests multiple monoclonal proteins do not clinically impact one another in the course of the disease.

Fanconi anemia (FA) is a rare recessive DNA repair deficiency resulting from mutations in one of at least 22 genes. Two-thirds of FA families harbor mutations in FANCA. To genotype patients in the International Fanconi Anemia Registry (IFAR) we employed multiple methodologies, screening 216 families for FANCA mutations. We describe identification of 57 large deletions and 261 sequence variants, in 159 families. All but seven families harbored distinct combinations of two mutations demonstrating high heterogeneity. Pathogenicity of the 18 novel missense variants was analyzed functionally by determining the ability of the mutant cDNA to improve the survival of a FANCA-null cell line when treated with MMC. Overexpressed pathogenic missense variants were found to reside in the cytoplasm, and nonpathogenic in the nucleus. RNA analysis demonstrated that two variants (c.522G > C and c.1565A > G), predicted to encode missense variants, which were determined to be nonpathogenic by a functional assay, caused skipping of exons 5 and 16 respectively, and are most likely pathogenic. We report 48 novel FANCA sequence variants. Defining both variants in a large patient cohort is a major step towards cataloging all FANCA variants, and permitting studies of genotype-phenotype correlations. This article is protected by copyright. All rights reserved

DATA ACCESSThe cDNA sequences reported in this paper have been deposited in the GenBank database (accession numbers): KF175165-KF175225 and BACs accession numbers that are used in this study: AC148621, AC190226, AC097327, AC097332, AC190226, AC097333, AC145401, AC187943, AC166855, AC166597, AC167295, AC235773, AC202644, and AC234805.ABSTRACTThe burst of segmental duplications during human and great ape evolution focuses on a set of “core” duplicons encoding great-ape-specific gene families. Characterization of these gene families is complicated by their high copy number, incomplete sequence, and polymorphic nature. We investigate the structure, transcriptional diversity, and protein localization of the nuclear pore complex-interacting protein (NPIP) or Morpheus gene family. The corresponding core, LCRA, encodes one of the most rapidly evolving genes in the human genome; LCRA has expanded to ~20 copies from a single ancestral locus in Old World monkey and is associated with most of the recurrent chromosome 16 microdeletions implicated in autism and mental retardation. Phylogenetic analysis and cDNA sequencing suggest two distinct subfamilies or subtypes, NPIPA and NPIPB. The latter expanded recently within the great apes due to a series of structural changes within the canonical gene structure. Among Old World monkey, we observe a testis-specific pattern of expression that contrasts with the ubiquitous pattern observed among human tissues. This change in the expression profile coincides with the structural events that reshaped the structure and organization of the gene family. Most of the expressed human copies are capable of producing an open reading frame. Immunofluorescence analyses of the morpheus genes showed a primary localization to both the nucleus and its periphery. We show that morpheus genes may be upregulated upon pI:C treatment and find evidence of human autoantibodies produced against the NPIPB protein, raising the possibility that morpheus genes may be related to immune- or autoimmune-related function.

Somatic hypermutation (SHM) plays a critical role in the maturation of antibodies, optimizing recognition initiated by recombination of V(D)J genes. Previous studies have shown that the propensity to mutate is modulated by the context of surrounding nucleotides and that SHM machinery generates biased substitutions. To investigate the intrinsic mutation frequency and substitution bias of SHMs at the amino acid level, we analyzed functional human antibody repertoires and developed mGSSP (method for Gene-Specific Substitution Profile), a method to construct amino acid substitution profiles from next-generation sequencing-acquired B cell transcripts. We demonstrated that these gene-specific substitution profiles (GSSP) are unique to each V gene and highly consistent between donors. We also showed that the GSSPs constructed from functional antibody repertoires are highly similar to those constructed from antibody sequences amplified from nonproductively rearranged passenger alleles, which do not undergo functional selection. This suggests that the types and frequencies, or mutational space, of a majority of amino acid changes sampled by the SHM machinery is well-captured by GSSPs. We further observed that the rates of mutational exchange between some amino acids are asymmetric and context-dependent, and correlate weakly with their biochemical properties. GSSPs provide an improved, position-dependent alternative to standard substitution matrices, for developing software for accurately modeling the SHM process. GSSPs can also be used for predicting the amino acid mutational space available for antigen-driven selection, and better understanding factors modulating the maturation pathways of antibody lineages in a gene-specific context. The mGSSP method can be used to build, compare, and plot GSSPs (https://github.com/scharch/SONAR/tree/master/sonar/mGSSP) and we report the GSSPs constructed for 69 common human V genes (DOI: 10.6084/m9.figshare.3511083) and provide high-resolution logo plots for each (DOI: 10.6084/m9.figshare.3511085).

Smith-Magenis syndrome (SMS), a neurodevelopmental disorder characterized by dysmorphic features, intellectual disability (ID), and sleep disturbances, results from a 17p11.2 microdeletion or a mutation in the RAI1 gene. We performed exome sequencing on 6 patients with SMS-like phenotypes but without chromosomal abnormalities or RAI1 variants. We identified pathogenic de novo variants in two cases, a nonsense variant in IQSEC2 and a missense variant in the SAND domain of DEAF1, and candidate de novo missense variants in an additional two cases. One candidate variant was located in an alpha helix of Necdin (NDN), phased to the paternally inherited allele. NDN is maternally imprinted within the 15q11.2 Prader-Willi Syndrome (PWS) region. This can help clarify NDN's role in the PWS phenotype. No definitive pathogenic gene variants were detected in the remaining SMS-like cases, but we report our findings for future comparison. This study provides information about the inheritance pattern and recurrence risk for patients with identified variants and demonstrates clinical and genetic overlap of neurodevelopmental disorders. Identification and characterization of ID-related genes that assist in development of common developmental pathways and/or gene-networks, may inform disease mechanism and treatment strategies.

Genome integrity of induced pluripotent stem cells (iPSCs) has been extensively studied in recent years, but it is still unclear whether iPSCs contain more genomic variations than cultured somatic cells. One important question is the origin of genomic variations detected in iPSCs–whether iPSC reprogramming induces such variations. Here, we undertook a unique approach by deriving fibroblast subclones and clonal iPSC lines from the same fibroblast population and applied next-generation sequencing to compare genomic variations in these lines. Targeted deep sequencing of parental fibroblasts revealed that most variants detected in clonal iPSCs and fibroblast subclones were rare variants inherited from the parental fibroblasts. Only a small number of variants remained undetectable in the parental fibroblasts, which were thus likely to be de novo. Importantly, the clonal iPSCs and fibroblast subclones contained comparable numbers of de novo variants. Collectively, our data suggest that iPSC reprogramming is not mutagenic.

Over the past 5 years, a new generation of highly potent and broadly neutralizing HIV-1 antibodies has been identified. These antibodies can protect against lentiviral infection in nonhuman primates (NHPs), suggesting that passive antibody transfer would prevent HIV-1 transmission in humans. To increase the protective efficacy of such monoclonal antibodies, we employed next-generation sequencing, computational bioinformatics, and structure-guided design to enhance the neutralization potency and breadth of VRC01, an antibody that targets the CD4 binding site of the HIV-1 envelope. One variant, VRC07-523, was 5- to 8-fold more potent than VRC01, neutralized 96% of viruses tested, and displayed minimal autoreactivity. To compare its protective efficacy to that of VRC01 in vivo , we performed a series of simian-human immunodeficiency virus (SHIV) challenge experiments in nonhuman primates and calculated the doses of VRC07-523 and VRC01 that provide 50% protection (EC 50 ). VRC07-523 prevented infection in NHPs at a 5-fold lower concentration than VRC01. These results suggest that increased neutralization potency in vitro correlates with improved protection against infection in vivo , documenting the improved functional efficacy of VRC07-523 and its potential clinical relevance for protecting against HIV-1 infection in humans. IMPORTANCE In the absence of an effective HIV-1 vaccine, alternative strategies are needed to block HIV-1 transmission. Direct administration of HIV-1-neutralizing antibodies may be able to prevent HIV-1 infections in humans. This approach could be especially useful in individuals at high risk for contracting HIV-1 and could be used together with antiretroviral drugs to prevent infection. To optimize the chance of success, such antibodies can be modified to improve their potency, breadth, and in vivo half-life. Here, knowledge of the structure of a potent neutralizing antibody, VRC01, that targets the CD4-binding site of the HIV-1 envelope protein was used to engineer a next-generation antibody with 5- to 8-fold increased potency in vitro . When administered to nonhuman primates, this antibody conferred protection at a 5-fold lower concentration than the original antibody. Our studies demonstrate an important correlation between in vitro assays used to evaluate the therapeutic potential of antibodies and their in vivo effectiveness.

In this commentary I have only highlighted a few dimensions that comparative genomics has reached into. Looking at a PubMed search of publications with the exact combination and order of the words “comparative genomics” in the title or abstract identifies 3752 articles as of the date of this writing. The chart in Figure ​Figure22 shows the growth of this field, which at first lagged in growth relative to the same search for “genomics,” but overall tracks this more general field of research. Other dimensions of comparative genomics, beyond the three areas I touched on above, include intraindividual comparative genomics (Cheng et al. 2012; Biesecker and Spinner 2013; Watson et al. 2013), human microbiome comparative genomics (Human Microbiome Project Consortium 2012) and how comparative genomics can shed light on a multidrug-resistant bacteria spread through a hospital (Snitkin et al. 2012). Clearly, as the field of genomics continues to expand, comparative genomics will always be an essential and central enabling component to help us discover and better understand the complexities, intricacies, and interrelatedness of the genomics of life. Figure 2 A PubMed search of publications with the exact combination and order of the words “comparative genomics” in the title or abstract identifies 3752 articles. This chart shows the growth of publications in this field year-by-year, and for ... Acknowledgments This research was supported by the Intramural Research Program of the National Human Genome Research Institute, National Institutes of Health.

Substantial intrastrain variation at the nucleotide level complicates molecular and genetic studies in zebrafish, such as the use of CRISPRs or morpholinos to inactivate genes. In the absence of robust inbred zebrafish lines, we generated NHGRI-1, a healthy and fecund strain derived from founder parents we sequenced to a depth of ∼50×. Within this strain, we have identified the majority of the genome that matches the reference sequence and documented most of the variants. This strain has utility for many reasons, but in particular it will be useful for any researcher who needs to know the exact sequence (with all variants) of a particular genomic region or who wants to be able to robustly map sequences back to a genome with all possible variants defined.

Coronary artery calcification (CAC) is a heritable and definitive morphologic marker of atherosclerosis that strongly predicts risk for future cardiovascular events. To search for genes involved in CAC, we used an integrative transcriptomic, genomic, and protein expression strategy by using next-generation DNA sequencing in the discovery phase with follow-up studies using traditional molecular biology and histopathology techniques. RNA sequencing of peripheral blood from a discovery set of CAC cases and controls was used to identify dysregulated genes, which were validated by ClinSeq and Framingham Heart Study data. Only a single gene, TREML4, was upregulated in CAC cases in both studies. Further examination showed that rs2803496 was a TREML4 cis-eQTL and that the minor allele at this locus conferred up to a 6.5-fold increased relative risk of CAC. We characterized human TREML4 and demonstrated by immunohistochemical techniques that it is localized in macrophages surrounding the necrotic core of coronary plaques complicated by calcification (but not in arteries with less advanced disease). Finally, we determined by von Kossa staining that TREML4 colocalizes with areas of microcalcification within coronary plaques. Overall, we present integrative RNA, DNA, and protein evidence implicating TREML4 in coronary artery calcification. Our findings connect multimodal genomics data with a commonly used clinical marker of cardiovascular disease.

Objective Early-onset epileptic encephalopathies have been associated with de novo mutations of numerous ion channel genes. We employed techniques of modern translational medicine to identify a disease-causing mutation, analyze its altered behavior, and screen for therapeutic compounds to treat the proband. Methods Three modern translational medicine tools were utilized: (1) high-throughput sequencing technology to identify a novel de novo mutation; (2) in vitro expression and electrophysiology assays to confirm the variant protein's dysfunction; and (3) screening of existing drug libraries to identify potential therapeutic compounds. Results A de novo GRIN2A missense mutation (c.2434C>A; p.L812M) increased the charge transfer mediated by N-methyl-D-aspartate receptors (NMDAs) containing the mutant GluN2A-L812M subunit. In vitro analysis with NMDA receptor blockers indicated that GLuN2A-L812M-containing NMDARs retained their sensitivity to the use-dependent channel blocker memantine; while screening of a previously reported GRIN2A mutation (N615K) with these compounds produced contrasting results. Consistent with these data, adjunct memantine therapy reduced our proband's seizure burden. Interpretation This case exemplifies the potential for personalized genomics and therapeutics to be utilized for the early diagnosis and treatment of infantile-onset neurological disease.

Background: Malignant hyperthermia susceptibility (MHS) is a life-threatening, inherited disorder of muscle calcium metabolism, triggered by anesthetics and depolarizing muscle relaxants. An unselected cohort was screened for MHS mutations using exome sequencing. The aim of this study was to pilot a strategy for the RYR1 and CACNA1S genes. Methods: Exome sequencing was performed on 870 volunteers not ascertained for MHS. Variants in RYR1 and CACNA1S were annotated using an algorithm that filtered results based on mutation type, frequency, and information in mutation databases. Variants were scored on a six-point pathogenicity scale. Medical histories and pedigrees were reviewed for malignant hyperthermia and related disorders. Results: The authors identified 70 RYR1 and 53 CACNA1S variants among 870 exomes. Sixty-three RYR1 and 41 CACNA1S variants passed the quality and frequency metrics but the authors excluded synonymous variants. In RYR1, the authors identified 65 missense mutations, one nonsense, two that affected splicing, and one non-frameshift indel. In CACNA1S, 48 missense, one frameshift deletion, one splicing, and one non-frameshift indel were identified. RYR1 variants predicted to be pathogenic for MHS were found in three participants without medical or family histories of MHS. Numerous variants, previously described as pathogenic in mutation databases, were reclassified by the authors as being of unknown pathogenicity. Conclusions: Exome sequencing can identify asymptomatic patients at risk for MHS, although the interpretation of exome variants can be challenging. The use of exome sequencing in unselected cohorts is an important tool to understand the prevalence and penetrance of MHS, a critical challenge for the field. hyperkalemia, as well some or all of the following symptoms: tachycardia, tachypnea, arrhythmias, skeletal muscle rigidity, and lethal hyperthermia. It is inherited in a predominately autosomal dominant pattern and associated with RYR1 or CACNA1S mutations, with other mapped loci. Seventy to 86% of patients with MHS have RYR1 mutations1-5 and 1% have CACNA1S mutations.6 The prevalence and penetrance of MHS mutations are difficult to determine because the pharmacologic exposure rate is low and it is an inconsistently manifesting gene-environment interaction; that is, when a susceptible patient is exposed to a triggering agent, the probability of malignant hyperthermia (MH) is less than 100%. Most MHS gene and variant studies have been performed on families with multiple generations affected with typical MHS. Studying these families made possible the discovery of the two implicated genes. However, these studies had ascertainment biases for those with severe reactions to the drugs. This has complicated efforts to establish the true prevalence and penetrance of MHS mutations. In addition, assigning pathogenicity to RYR1 and CACNA1S variants is challenging for several reasons. First is the issue of locus heterogeneity. With several mapped loci without identified genes, some RYR1 and CACNA1S variants may have been erroneously determined to be pathogenic when there was a causative variant in another (untested) gene. In addition, RYR1 and CACNA1S are large genes with 106 and 44 exons, respectively, making mutation screening challenging. Thus, some RYR1 and CACNA1S variants previously determined to be pathogenic may be benign, as has been shown for other genes.7 New sequencing technologies, including exome sequencing (ES), have made sequencing of the human exome (exons of known genes) feasible. This provides the opportunity to detect mutations in MHS genes in a less biased manner. Using this approach, we can improve our understanding of the mutational spectra of the RYR1 and CACNA1S genes, and estimate their penetrance. Our objective was to identify mutations in RYR1 and CACNA1S in a population not ascertained for MHS, as a pilot for the use of exome data for predictive medicine.

Chromatin-based functional genomic analyses and genomewide association studies (GWASs) together implicate enhancers as critical elements influencing gene expression and risk for common diseases. Here, we performed systematic chromatin and transcriptome profiling in human pancreatic islets. Integrated analysis of islet data with those from nine cell types identified specific and significant enrichment of type 2 diabetes and related quantitative trait GWAS variants in islet enhancers. Our integrated chromatin maps reveal that most enhancers are short (median = 0.8 kb). Each cell type also contains a substantial number of more extended (≥ 3 kb) enhancers. Interestingly, these stretch enhancers are often tissue-specific and overlap locus control regions, suggesting that they are important chromatin regulatory beacons. Indeed, we show that (i) tissue specificity of enhancers and nearby gene expression increase with enhancer length; (ii) neighborhoods containing stretch enhancers are enriched for important cell type-specific genes; and (iii) GWAS variants associated with traits relevant to a particular cell type are more enriched in stretch enhancers compared with short enhancers. Reporter constructs containing stretch enhancer sequences exhibited tissue-specific activity in cell culture experiments and in transgenic mice. These results suggest that stretch enhancers are critical chromatin elements for coordinating cell type-specific regulatory programs and that sequence variation in stretch enhancers affects risk of major common human diseases.

SummaryAntibodies of the VRC01 class neutralize HIV-1, arise in diverse HIV-1-infected donors, and are potential templates for an effective HIV-1 vaccine. However, the stochastic processes that generate repertoires in each individual of >1012 antibodies make elicitation of specific antibodies uncertain. Here we determine the ontogeny of the VRC01 class by crystallography and next-generation sequencing. Despite antibody-sequence differences exceeding 50%, antibody-gp120 cocrystal structures reveal VRC01-class recognition to be remarkably similar. B cell transcripts indicate that VRC01-class antibodies require few specific genetic elements, suggesting that naive-B cells with VRC01-class features are generated regularly by recombination. Virtually all of these fail to mature, however, with only a few—likely one—ancestor B cell expanding to form a VRC01-class lineage in each donor. Developmental similarities in multiple donors thus reveal the generation of VRC01-class antibodies to be reproducible in principle, thereby providing a framework for attempts to elicit similar antibodies in the general population.

Synonymous mutations, which do not alter the protein sequence, have been shown to affect protein function [Sauna ZE, Kimchi-Sarfaty C (2011) Nat Rev Genet 12(10):683–691]. However, synonymous mutations are rarely investigated in the cancer genomics field. We used whole-genome and -exome sequencing to identify somatic mutations in 29 melanoma samples. Validation of one synonymous somatic mutation in BCL2L12 in 285 samples identified 12 cases that harbored the recurrent F17F mutation. This mutation led to increased BCL2L12 mRNA and protein levels because of differential targeting of WT and mutant BCL2L12 by hsa-miR-671–5p. Protein made from mutant BCL2L12 transcript bound p53, inhibited UV-induced apoptosis more efficiently than WT BCL2L12 , and reduced endogenous p53 target gene transcription. This report shows selection of a recurrent somatic synonymous mutation in cancer. Our data indicate that silent alterations have a role to play in human cancer, emphasizing the importance of their investigation in future cancer genome studies.