ViroFind: A novel target-enrichment deep-sequencing platform reveals a complex JC virus population in the brain of PML patients

Deep nucleotide sequencing enables the unbiased, broad-spectrum detection of viruses in clinical samples without requiring an a priori hypothesis for the source of infection. However, its use in clinical research applications is limited by low cost-effectiveness given that most of the sequencing information from clinical samples is related to the human genome, which renders the analysis of viral genomes challenging. To overcome this limitation we developed ViroFind, an in-solution target-enrichment platform for virus detection and discovery in clinical samples. ViroFind comprises 165,433 viral probes that cover the genomes of 535 selected DNA and RNA viruses that infect humans or could cause zoonosis. The ViroFind probes are used in a hybridization reaction to enrich viral sequences and therefore enhance the detection of viral genomes via deep sequencing. We used ViroFind to detect and analyze all viral populations in the brain of 5 patients with progressive multifocal leukoencephalopathy (PML) and of 18 control subjects with no known neurological disease. Compared to direct deep sequencing, by using ViroFind we enriched viral sequences present in the clinical samples up to 127-fold. We discovered highly complex polyoma virus JC populations in the PML brain samples with a remarkable degree of genetic divergence among the JC virus variants of each PML brain sample. Specifically for the viral capsid protein VP1 gene, we identified 24 single nucleotide substitutions, 12 of which were associated with amino acid changes. The most frequent (4 of 5 samples, 80%) amino acid change was D66H, which is associated with enhanced tissue tropism, and hence likely a viral fitness advantage, compared to other variants. Lastly, we also detected sparse JC virus sequences in 10 of 18 (55.5%) of control samples and sparse human herpes virus 6B (HHV6B) sequences in the brain of 11 of 18 (61.1%) control subjects. In sum, ViroFind enabled the in-depth analysis of all viral genomes in PML and control brain samples and allowed us to demonstrate a high degree of JC virus genetic divergence in vivo that has been previously underappreciated. ViroFind can be used to investigate the structure of the virome with unprecedented depth in health and disease state

Loss of FHL1 induces an age-dependent skeletal muscle myopathy associated with myofibrillar and intermyofibrillar disorganization in mice

Recent human genetic studies have provided evidences that sporadic or inherited missense mutations in four-and-a-half LIM domain protein 1 (FHL1), resulting in alterations in FHL1 protein expression, are associated with rare congenital myopathies, including reducing body myopathy and Emery–Dreifuss muscular dystrophy. However, it remains to be clarified whether mutations in FHL1 cause skeletal muscle remodeling owing to gain- or loss of FHL1 function. In this study, we used FHL1-null mice lacking global FHL1 expression to evaluate loss-of-function effects on skeletal muscle homeostasis. Histological and functional analyses of soleus, tibialis anterior and sternohyoideus muscles demonstrated that FHL1-null mice develop an age-dependent myopathy associated with myofibrillar and intermyofibrillar (mitochondrial and sarcoplasmic reticulum) disorganization, impaired muscle oxidative capacity and increased autophagic activity. A longitudinal study established decreased survival rates in FHL1-null mice, associated with age-dependent impairment of muscle contractile function and a significantly lower exercise capacity. Analysis of primary myoblasts isolated from FHL1-null muscles demonstrated early muscle fiber differentiation and maturation defects, which could be rescued by re-expression of the FHL1A isoform, highlighting that FHL1A is necessary for proper muscle fiber differentiation and maturation in vitro. Overall, our data show that loss of FHL1 function leads to myopathy in vivo and suggest that loss of function of FHL1 may be one of the mechanisms underlying muscle dystrophy in patients with FHL1 mutations

Burden of Rare Sarcomere Gene Variants in the Framingham and Jackson Heart Study Cohorts

Rare sarcomere protein variants cause dominant hypertrophic and dilated cardiomyopathies. To evaluate whether allelic variants in eight sarcomere genes are associated with cardiac morphology and function in the community, we sequenced 3,600 individuals from the Framingham Heart Study (FHS) and Jackson Heart Study (JHS) cohorts. Out of the total, 11.2% of individuals had one or more rare nonsynonymous sarcomere variants. The prevalence of likely pathogenic sarcomere variants was 0.6%, twice the previous estimates; however, only four of the 22 individuals had clinical manifestations of hypertrophic cardiomyopathy. Rare sarcomere variants were associated with an increased risk for adverse cardiovascular events (hazard ratio: 2.3) in the FHS cohort, suggesting that cardiovascular risk assessment in the general population can benefit from rare variant analysis

Cardiomyocyte Proliferative Capacity Is Restricted in Mice With Lmna Mutation

LMNA is one of the leading causative genes of genetically inherited dilated cardiomyopathy (DCM). Unlike most DCM-causative genes, which encode sarcomeric or sarcomere-related proteins, LMNA encodes nuclear envelope proteins, lamin A and C, and does not directly associate with contractile function. However, a mutation in this gene could lead to the development of DCM. The molecular mechanism of how LMNA mutation contributes to DCM development remains largely unclear and yet to be elucidated. The objective of this study was to clarify the mechanism of developing DCM caused by LMNA mutation.Methods and Results: We assessed cardiomyocyte phenotypes and characteristics focusing on cell cycle activity in mice with Lmna mutation. Both cell number and cell size were reduced, cardiomyocytes were immature, and cell cycle activity was retarded in Lmna mutant mice at both 5 weeks and 2 years of age. RNA-sequencing and pathway analysis revealed “proliferation of cells” had the most substantial impact on Lmna mutant mice. Cdkn1a, which encodes the cell cycle regulating protein p21, was strongly upregulated in Lmna mutants, and upregulation of p21 was confirmed by Western blot and immunostaining. DNA damage, which is known to upregulate Cdkn1a, was more abundantly detected in Lmna mutant mice. To assess the proliferative capacity of cardiomyocytes, the apex of the neonate mouse heart was resected, and recovery from the insult was observed. A restricted cardiomyocyte proliferating capacity after resecting the apex of the heart was observed in Lmna mutant mice.Conclusions: Our results strongly suggest that loss of lamin function contributes to impaired cell proliferation through cell cycle defects. The inadequate inborn or responsive cell proliferation capacity plays an essential role in developing DCM with LMNA mutation

De novo mutations in histone modifying genes in congenital heart disease

Congenital heart disease (CHD) is the most frequent birth defect, affecting 0.8% of live births1. Many cases occur sporadically and impair reproductive fitness, suggesting a role for de novo mutations. By analysis of exome sequencing of parent-offspring trios, we compared the incidence of de novo mutations in 362 severe CHD cases and 264 controls. CHD cases showed a significant excess of protein-altering de novo mutations in genes expressed in the developing heart, with an odds ratio of 7.5 for damaging mutations. Similar odds ratios were seen across major classes of severe CHD. We found a marked excess of de novo mutations in genes involved in production, removal or reading of H3K4 methylation (H3K4me), or ubiquitination of H2BK120, which is required for H3K4 methylation2–4. There were also two de novo mutations in SMAD2; SMAD2 signaling in the embryonic left-right organizer induces demethylation of H3K27me5. H3K4me and H3K27me mark `poised' promoters and enhancers that regulate expression of key developmental genes6. These findings implicate de novo point mutations in several hundred genes that collectively contribute to ~10% of severe CHD

Integrated genomic characterization of oesophageal carcinoma

Oesophageal cancers are prominent worldwide; however, there are few targeted therapies and survival rates for these cancers remain dismal. Here we performed a comprehensive molecular analysis of 164 carcinomas of the oesophagus derived from Western and Eastern populations. Beyond known histopathological and epidemiologic distinctions, molecular features differentiated oesophageal squamous cell carcinomas from oesophageal adenocarcinomas. Oesophageal squamous cell carcinomas resembled squamous carcinomas of other organs more than they did oesophageal adenocarcinomas. Our analyses identified three molecular subclasses of oesophageal squamous cell carcinomas, but none showed evidence for an aetiological role of human papillomavirus. Squamous cell carcinomas showed frequent genomic amplifications of CCND1 and SOX2 and/or TP63, whereas ERBB2, VEGFA and GATA4 and GATA6 were more commonly amplified in adenocarcinomas. Oesophageal adenocarcinomas strongly resembled the chromosomally unstable variant of gastric adenocarcinoma, suggesting that these cancers could be considered a single disease entity. However, some molecular features, including DNA hypermethylation, occurred disproportionally in oesophageal adenocarcinomas. These data provide a framework to facilitate more rational categorization of these tumours and a foundation for new therapies

Comprehensive molecular characterization of gastric adenocarcinoma

Gastric cancer is a leading cause of cancer deaths, but analysis of its molecular and clinical characteristics has been complicated by histological and aetiological heterogeneity. Here we describe a comprehensive molecular evaluation of 295 primary gastric adenocarcinomas as part of The Cancer Genome Atlas (TCGA) project. We propose a molecular classification dividing gastric cancer into four subtypes: tumours positive for Epstein–Barr virus, which display recurrent PIK3CA mutations, extreme DNA hypermethylation, and amplification of JAK2, CD274 (also known as PD-L1) and PDCD1LG2 (also knownasPD-L2); microsatellite unstable tumours, which show elevated mutation rates, including mutations of genes encoding targetable oncogenic signalling proteins; genomically stable tumours, which are enriched for the diffuse histological variant and mutations of RHOA or fusions involving RHO-family GTPase-activating proteins; and tumours with chromosomal instability, which show marked aneuploidy and focal amplification of receptor tyrosine kinases. Identification of these subtypes provides a roadmap for patient stratification and trials of targeted therapies

A functional genomic approach to actionable gene fusions for precision oncology

Fusion genes represent a class of attractive therapeutic targets. Thousands of fusion genes have been identified in patients with cancer, but the functional consequences and therapeutic implications of most of these remain largely unknown. Here, we develop a functional genomic approach that consists of efficient fusion reconstruction and sensitive cell viability and drug response assays. Applying this approach, we characterize similar to 100 fusion genes detected in patient samples of The Cancer Genome Atlas, revealing a notable fraction of low-frequency fusions with activating effects on tumor growth. Focusing on those in the RTK-RAS pathway, we identify a number of activating fusions that can markedly affect sensitivity to relevant drugs. Last, we propose an integrated, level-of-evidence classification system to prioritize gene fusions systematically. Our study reiterates the urgent clinical need to incorporate similar functional genomic approaches to characterize gene fusions, thereby maximizing the utility of gene fusions for precision oncology

The Somatic Genomic Landscape of Chromophobe Renal Cell Carcinoma

We describe the landscape of somatic genomic alterations of 66 chromophobe renal cell carcinomas (ChRCCs) based on multidimensional and comprehensive characterization, including mitochondrial DNA (mtDNA) and whole genome sequencing. The result is consistent that ChRCC originates from the distal nephron compared to other kidney cancers with more proximal origins. Combined mtDNA and gene expression analysis implicates changes in mitochondrial function as a component of the disease biology, while suggesting alternative roles for mtDNA mutations in cancers relying on oxidative phosphorylation. Genomic rearrangements lead to recurrent structural breakpoints within TERT promoter region, which correlates with highly elevated TERT expression and manifestation of kataegis, representing a mechanism of TERT up-regulation in cancer distinct from previously-observed amplifications and point mutations

Multiplatform Analysis of 12 Cancer Types Reveals Molecular Classification within and across Tissues of Origin

Recent genomic analyses of pathologically-defined tumor types identify “within-a-tissue” disease subtypes. However, the extent to which genomic signatures are shared across tissues is still unclear. We performed an integrative analysis using five genome-wide platforms and one proteomic platform on 3,527 specimens from 12 cancer types, revealing a unified classification into 11 major subtypes. Five subtypes were nearly identical to their tissue-of-origin counterparts, but several distinct cancer types were found to converge into common subtypes. Lung squamous, head & neck, and a subset of bladder cancers coalesced into one subtype typified by TP53 alterations, TP63 amplifications, and high expression of immune and proliferation pathway genes. Of note, bladder cancers split into three pan-cancer subtypes. The multi-platform classification, while correlated with tissue-of-origin, provides independent information for predicting clinical outcomes. All datasets are available for data-mining from a unified resource to support further biological discoveries and insights into novel therapeutic strategies