Relationship between BCL2 mutations and follicular lymphoma outcome in the chemoimmunotherapy era.

How to identify follicular lymphoma (FL) patients with low disease burden but high risk for early progression is unclear. Building on a prior study demonstrating the early transformation of FLs with high variant allele frequency (VAF) BCL2 mutations at activation-induced cytidine deaminase (AICDA) sites, we examined 11 AICDA mutational targets, including BCL2, BCL6, PAX5, PIM1, RHOH, SOCS, and MYC, in 199 newly diagnosed grade 1 and 2 FLs. BCL2 mutations with VAF ≥20% occurred in 52% of cases. Among 97 FL patients who did not initially receive rituximab-containing therapy, nonsynonymous BCL2 mutations at VAF ≥20% were associated with increased transformation risk (HR 3.01, 95% CI 1.04-8.78, p = 0.043) and a trend toward shorter event-free survival (EFS, median 20 months with mutations versus 54 months without, p = 0.052). Other sequenced genes were less frequently mutated and did not increase the prognostic value of the panel. Across the entire population, nonsynonymous BCL2 mutations at VAF ≥20% were associated with decreased EFS (HR 1.55, 95% CI 1.02-2.35, p = 0.043 after correction for FLIPI and treatment) and decreased overall survival after median 14-year follow-up (HR 1.82, 95% CI 1.05-3.17, p = 0.034). Thus, high VAF nonsynonymous BCL2 mutations remain prognostic even in the chemoimmunotherapy era.

Non-random genetic alterations in the cyanobacterium Nostoc sp. exposed to space conditions.

Understanding the impact of long-term exposure of microorganisms to space is critical in understanding how these exposures impact the evolution and adaptation of microbial life under space conditions. In this work we subjected Nostoc sp. CCCryo 231-06, a cyanobacterium capable of living under many different ecological conditions, and also surviving in extreme ones, to a 23-month stay at the International Space Station (the Biology and Mars Experiment, BIOMEX, on the EXPOSE-R2 platform) and returned it to Earth for single-cell genome analysis. We used microfluidic technology and single cell sequencing to identify the changes that occurred in the whole genome of single Nostoc cells. The variant profile showed that biofilm and photosystem associated loci were the most altered, with an increased variant rate of synonymous base pair substitutions. The cause(s) of these non-random alterations and their implications to the evolutionary potential of single bacterial cells under long-term cosmic exposure warrants further investigation.

Whole genome sequencing of cyanobacterium Nostoc sp. CCCryo 231-06 using microfluidic single cell technology.

The Nostoc sp. strain CCCryo 231-06 is a cyanobacterial strain capable of surviving under extreme conditions and thus is of great interest for the astrobiology community. The knowledge of its complete genome sequence would serve as a guide for further studies. However, a major concern has been placed on the effects of contamination on the quality of sequencing data without a reference genome. Here, we report the use of microfluidic technology combined with single cell sequencing and de novo assembly to minimize the contamination and recover the complete genome of the Nostoc strain CCCryo 231-06 with high quality. 100% of the whole genome was recovered with all contaminants removed and a strongly supported phylogenetic tree. The data reported can be useful for comparative genomics for phylogenetic and taxonomic studies. The method used in this work can be applied to studies that require high-quality assemblies of genomes of unknown microorganisms.

RNA-Seq Reveals Differences in Expressed Tumor Mutation Burden in Colorectal and Endometrial Cancers with and without Defective DNA-Mismatch Repair.

Tumor mutation burden (TMB) is an emerging biomarker of immunotherapy response. RNA sequencing in FFPE tissue samples was used for determining TMB in microsatellite-stable (MSS) and microsatellite instability-high (MSI-H) tumors in patients with colorectal or endometrial cancer. Tissue from tumors and paired normal tissue from 46 MSI-H and 12 MSS cases were included. Of the MSI-H tumors, 29 had defective DNA mismatch-repair mutations, and 17 had MLH1 promoter hypermethylation. TMB was measured using the expressed somatic nucleotide variants (eTMB). A method of accurate measurement of eTMB was developed that removes FFPE-derived artifacts by leveraging mutation signatures. There was a significant difference in the median eTMB values observed between MSI-H and MSS cases: 27.3 versus 6.7 mutations/megabase (mut/Mb) (P = 3.5 × 10-9). Among tumors with defective DNA-mismatch repair, those with mismatch-repair mutations had a significantly higher median eTMB than those with hypermethylation: 28.1 versus 17.5 mut/Mb (P = 0.037). Multivariate analysis showed that MSI status, tumor type (endometrial or colorectal), and age were significantly associated with eTMB. Additionally, using whole-exome sequencing in a subset of these patients, it was determined that DNA TMB correlated well with eTMB (Spearman correlation coefficient, 0.83). These results demonstrate that RNA sequencing can be used for measuring eTMB in FFPE tumor specimens.

Chromosomal Junction Detection from Whole-Genome Sequencing on Formalin-Fixed, Paraffin-Embedded Tumors.

DNA junctions (DNAJs) frequently impact clinically relevant genes in tumors and are important for diagnostic and therapeutic purposes. Although routinely screened through fluorescence in situ hybridization assays, such testing only allows the interrogation of single-gene regions or known fusion partners. Comprehensive assessment of DNAJs present across the entire genome can only be determined from whole-genome sequencing. Structural variance analysis from whole-genome paired-end sequencing data is, however, frequently restricted to copy number changes without DNAJ detection. Through optimized whole-genome sequencing and specialized bioinformatics algorithms, complete structural variance analysis is reported, including DNAJs, from formalin-fixed DNA. Selective library assembly from larger fragments (>500 bp) and economical sequencing depths (300 to 400 million reads) provide representative genomic coverage profiles and increased allelic coverage to levels compatible with DNAJ calling (40× to 60×). Although consistently fragmented, more recently formalin-fixed, specimens (<2 years' storage) revealed consistent populations of larger DNA fragments. Optimized bioinformatics efficiently detected >90% of DNAJs in two prostate tumors (approximately 60% tumor) previously analyzed by mate-pair sequencing on fresh frozen tissue, with evidence of at least one spanning-read in 99% of DNAJs. Rigorous masking with data from unrelated formalin-fixed tissue progressively eliminated many false-positive DNAJs, without loss of true positives, resulting in low numbers of false-positive passing current filters. This methodology enables more comprehensive clinical genomics testing on formalin-fixed clinical specimens.

Next-Generation Sequencing of CYP2C19 in Stent Thrombosis: Implications for Clopidogrel Pharmacogenomics.

PURPOSE: Describe CYP2C19 sequencing results in the largest series of clopidogrel-treated cases with stent thrombosis (ST), the closest clinical phenotype to clopidogrel resistance. Evaluate the impact of CYP2C19 genetic variation detected by next-generation sequencing (NGS) with comprehensive annotation and functional studies. METHODS: Seventy ST cases on clopidogrel identified from the PLATO trial (n = 58) and Mayo Clinic biorepository (n = 12) were matched 1:1 with controls for age, race, sex, diabetes mellitus, presentation, and stent type. NGS was performed to cover the entire CYP2C19 gene. Assessment of exonic variants involved measuring in vitro protein expression levels. Intronic variants were evaluated for potential splicing motif variations. RESULTS: Poor metabolizers (n = 4) and rare CYP2C19*8, CYP2C19*15, and CYP2C19*11 alleles were identified only in ST cases. CYP2C19*17 heterozygote carriers were observed more frequently in cases (n = 29) than controls (n = 18). Functional studies of CYP2C19 exonic variants (n = 11) revealed 3 cases and only 1 control carrying a deleterious variant as determined by in vitro protein expression studies. Greater intronic variation unique to ST cases (n = 169) compared with controls (n = 84) was observed with predictions revealing 13 allele candidates that may lead to a potential disruption of splicing and a loss-of-function effect of CYP2C19 in ST cases. CONCLUSION: NGS detected CYP2C19 poor metabolizers and paradoxically greater number of so-called rapid metabolizers in ST cases. Rare deleterious exonic variation occurs in 4%, and potentially disruptive intronic alleles occur in 16% of ST cases. Additional studies are required to evaluate the role of these variants in platelet aggregation and clopidogrel metabolism.

Comparative evaluation for the globin gene depletion methods for mRNA sequencing using the whole blood-derived total RNAs.

BACKGROUND: There are challenges in generating mRNA-Seq data from whole-blood derived RNA as globin gene and rRNA are frequent contaminants. Given the abundance of erythrocytes in whole blood, globin genes comprise some 80% or more of the total RNA. Therefore, depletion of globin gene RNA and rRNA are critical steps required to have adequate coverage of reads mapping to the reference transcripts and thus reduce the total cost of sequencing. In this study, we directly compared the performance of probe hybridization (GLOBINClear Kit and Globin-Zero Gold rRNA Removal Kit) and RNAse-H enzymatic depletion (NEBNext® Globin & rRNA Depletion Kit and Ribo-Zero Plus rRNA Depletion Kit) methods from 1 µg of whole blood-derived RNA on mRNA-Seq profiling. All RNA samples were treated with DNaseI for additional cleanup before the depletion step and were processed for poly-A selection for library generation. RESULTS: Probe hybridization revealed a better overall performance than the RNAse-H enzymatic depletion method, detecting a higher number of genes and transcripts without 3' region bias. After depletion, samples treated with probe hybridization showed globin genes at 0.5% (±0.6%) of the total mapped reads; the RNAse-H enzymatic depletion had 3.2% (±3.8%). Probe hybridization showed more junction reads and transcripts compared with RNAse-H enzymatic depletion and also had a higher correlation (R > 0.9) than RNAse-H enzymatic depletion (R > 0.85). CONCLUSION: In this study, our results showed that 1 µg of high-quality RNA from whole blood could be routinely used for transcriptional profiling analysis studies with globin gene and rRNA depletion pre-processing. We also demonstrated that the probe hybridization depletion method is better suited to mRNA sequencing analysis with minimal effect on RNA quality during depletion procedures.

Bacterial Single Cell Whole Transcriptome Amplification in Microfluidic Platform Shows Putative Gene Expression Heterogeneity.

Single cell RNA sequencing is a technology that provides the capability of analyzing the transcriptome of a single cell from a population. So far, single cell RNA sequencing has been focused mostly on human cells due to the larger starting amount of RNA template for subsequent amplification. One of the major challenges of applying single cell RNA sequencing to microbial cells is to amplify the femtograms of the RNA template to obtain sufficient material for downstream sequencing with minimal contamination. To achieve this goal, efforts have been focused on multiround RNA amplification, but would introduce additional contamination and bias. In this work, we for the first time coupled a microfluidic platform with multiple displacement amplification technology to perform single cell whole transcriptome amplification and sequencing of Porphyromonas somerae, a microbe of interest in endometrial cancer, as a proof-of-concept demonstration of using single cell RNA sequencing tool to unveil gene expression heterogeneity in single microbial cells. Our results show that the bacterial single-cell gene expression regulation is distinct across different cells, supporting widespread heterogeneity.

Recurrent MSC E116K mutations in ALK-negative anaplastic large cell lymphoma.

Anaplastic large cell lymphomas (ALCLs) represent a relatively common group of T-cell non-Hodgkin lymphomas (T-NHLs) that are unified by similar pathologic features but demonstrate marked genetic heterogeneity. ALCLs are broadly classified as being anaplastic lymphoma kinase (ALK)+ or ALK-, based on the presence or absence of ALK rearrangements. Exome sequencing of 62 T-NHLs identified a previously unreported recurrent mutation in the musculin gene, MSC E116K, exclusively in ALK- ALCLs. Additional sequencing for a total of 238 T-NHLs confirmed the specificity of MSC E116K for ALK- ALCL and further demonstrated that 14 of 15 mutated cases (93%) had coexisting DUSP22 rearrangements. Musculin is a basic helix-loop-helix (bHLH) transcription factor that heterodimerizes with other bHLH proteins to regulate lymphocyte development. The E116K mutation localized to the DNA binding domain of musculin and permitted formation of musculin-bHLH heterodimers but prevented their binding to authentic target sequence. Functional analysis showed MSCE116K acted in a dominant-negative fashion, reversing wild-type musculin-induced repression of MYC and cell cycle inhibition. Chromatin immunoprecipitation-sequencing and transcriptome analysis identified the cell cycle regulatory gene E2F2 as a direct transcriptional target of musculin. MSCE116K reversed E2F2-induced cell cycle arrest and promoted expression of the CD30-IRF4-MYC axis, whereas its expression was reciprocally induced by binding of IRF4 to the MSC promoter. Finally, ALCL cells expressing MSC E116K were preferentially targeted by the BET inhibitor JQ1. These findings identify a novel recurrent MSC mutation as a key driver of the CD30-IRF4-MYC axis and cell cycle progression in a unique subset of ALCLs.

A Prospective Correlation of Tissue Histopathology With Nucleic Acid Yield in Metastatic Castration-Resistant Prostate Cancer Biopsy Specimens.

OBJECTIVE: To determine histopathologic, exome, and transcriptome nucleic acid material yield from prospectively collected metastatic tissue biopsy specimens in patients with metastatic castration-resistant prostate cancer (mCRPC). PATIENTS AND METHODS: Patients with mCRPC initiating abiraterone acetate therapy underwent 2 serial metastatic site core needle biopsies after study activation on May 17, 2013. Multiple cores were obtained, and from each core, 1- to 2-mm segments were separated and formalin fixed for histopathologic examination. Tumor purity was determined for DNA and RNA from the rest of the biopsy specimen. RNA quality was assessed by calculation of an RNA integrity number and a DV200 score. RESULTS: A total of 89 patients underwent 172 uniformly processed core needle biopsies (89 on visit 1 and 83 on visit 2) between May 30, 2013, and September 10, 2015. Metastatic sites biopsied included bone (131), lymph nodes (31), liver (5), lung (3), and pelvic soft tissues (2). Of the 172 biopsy specimens, 85 (49%) had at least one of the multiple cores positive for tumor on histopathologic examination (53 of 88 [60%] from visit 1 and 32 of 83 [39%] from visit 2; P=.006). Metastatic carcinoma was observed in 50 of 130 bone lesion specimens (38%), compared to 35 of 41 nonbone specimens (85%) (P<.001). More than 10% tumoral DNA purity was observed in 89% and 80% of visit 1 and visit 2 biopsy specimens, respectively. Similarly, more than 10% tumor RNA purity was observed in 79% of visit 1 vs 59% for visit 2 (P=.008). In all, 134 of 172 procedures (78%) yielded tumor material either by histopathologic or nucleic acid purity analysis. CONCLUSION: This study found that biopsy specimens from mCRPC sites yield adequate histopathologic, exome, and transcriptome material in most, but not all, cases. This finding has relevance for future genome sequencing studies on the introduction of targeted therapeutic agents. TRIAL REGISTRATION: clinicaltrials.gov Identifier: 01953640.

NUP155 insufficiency recalibrates a pluripotent transcriptome with network remodeling of a cardiogenic signaling module.

BACKGROUND: Atrial fibrillation is a cardiac disease driven by numerous idiopathic etiologies. NUP155 is a nuclear pore complex protein that has been identified as a clinical driver of atrial fibrillation, yet the precise mechanism is unknown. The present study employs a systems biology algorithm to identify effects of NUP155 disruption on cardiogenicity in a model of stem cell-derived differentiation. METHODS: Embryonic stem (ES) cell lines (n = 5) with truncated NUP155 were cultured in parallel with wild type (WT) ES cells (n = 5), and then harvested for RNAseq. Samples were run on an Illumina HiSeq 2000. Reads were analyzed using Strand NGS, Cytoscape, DAVID and Ingenuity Pathways Analysis to deconvolute the NUP155-disrupted transcriptome. Network topological analysis identified key features that controlled framework architecture and functional enrichment. RESULTS: In NUP155 truncated ES cells, significant expression changes were detected in 326 genes compared to WT. These genes segregated into clusters that enriched for specific gene ontologies. Deconvolution of the collective framework into discrete sub-networks identified a module with the highest score that enriched for Cardiovascular System Development, and revealed NTRK1/TRKA and SRSF2/SC35 as critical hubs within this cardiogenic module. CONCLUSIONS: The strategy of pluripotent transcriptome deconvolution used in the current study identified a novel association of NUP155 with potential drivers of arrhythmogenic AF. Here, NUP155 regulates cardioplasticity of a sub-network embedded within a larger framework of genome integrity, and exemplifies how transcriptome cardiogenicity in an embryonic stem cell genome is recalibrated by nucleoporin dysfunction.

Targeted sequencing of 36 known or putative colorectal cancer susceptibility genes.

BACKGROUND: Mutations in several genes predispose to colorectal cancer. Genetic testing for hereditary colorectal cancer syndromes was previously limited to single gene tests; thus, only a very limited number of genes were tested, and rarely those infrequently mutated in colorectal cancer. Next-generation sequencing technologies have made it possible to sequencing panels of genes known and suspected to influence colorectal cancer susceptibility. METHODS: Targeted sequencing of 36 known or putative CRC susceptibility genes was conducted for 1231 CRC cases from five subsets: (1) Familial Colorectal Cancer Type X (n = 153); (2) CRC unselected by tumor immunohistochemical or microsatellite stability testing (n = 548); (3) young onset (age <50 years) (n = 333); (4) proficient mismatch repair (MMR) in cases diagnosed at ≥50 years (n = 68); and (5) deficient MMR CRCs with no germline mutations in MLH1, MSH2, MSH6, or PMS2 (n = 129). Ninety-three unaffected controls were also sequenced. RESULTS: Overall, 29 nonsense, 43 frame-shift, 13 splice site, six initiator codon variants, one stop codon, 12 exonic deletions, 658 missense, and 17 indels were identified. Missense variants were reviewed by genetic counselors to determine pathogenicity; 13 were pathogenic, 61 were not pathogenic, and 584 were variants of uncertain significance. Overall, we identified 92 cases with pathogenic mutations in APC,MLH1,MSH2,MSH6, or multiple pathogenic MUTYH mutations (7.5%). Four cases with intact MMR protein expression by immunohistochemistry carried pathogenic MMR mutations. CONCLUSIONS: Results across case subsets may help prioritize genes for inclusion in clinical gene panel tests and underscore the issue of variants of uncertain significance both in well-characterized genes and those for which limited experience has accumulated.

Androgen Receptor Variant AR-V9 Is Coexpressed with AR-V7 in Prostate Cancer Metastases and Predicts Abiraterone Resistance.

Purpose: Androgen receptor (AR) variant AR-V7 is a ligand-independent transcription factor that promotes prostate cancer resistance to AR-targeted therapies. Accordingly, efforts are under way to develop strategies for monitoring and inhibiting AR-V7 in castration-resistant prostate cancer (CRPC). The purpose of this study was to understand whether other AR variants may be coexpressed with AR-V7 and promote resistance to AR-targeted therapies.Experimental Design: We utilized complementary short- and long-read sequencing of intact AR mRNA isoforms to characterize AR expression in CRPC models. Coexpression of AR-V7 and AR-V9 mRNA in CRPC metastases and circulating tumor cells was assessed by RNA-seq and RT-PCR, respectively. Expression of AR-V9 protein in CRPC models was evaluated with polyclonal antisera. Multivariate analysis was performed to test whether AR variant mRNA expression in metastatic tissues was associated with a 12-week progression-free survival endpoint in a prospective clinical trial of 78 CRPC-stage patients initiating therapy with the androgen synthesis inhibitor, abiraterone acetate.Results: AR-V9 was frequently coexpressed with AR-V7. Both AR variant species were found to share a common 3' terminal cryptic exon, which rendered AR-V9 susceptible to experimental manipulations that were previously thought to target AR-V7 uniquely. AR-V9 promoted ligand-independent growth of prostate cancer cells. High AR-V9 mRNA expression in CRPC metastases was predictive of primary resistance to abiraterone acetate (HR = 4.0; 95% confidence interval, 1.31-12.2; P = 0.02).Conclusions: AR-V9 may be an important component of therapeutic resistance in CRPC. Clin Cancer Res; 23(16); 4704-15. ©2017 AACR.

Integrated mate-pair and RNA sequencing identifies novel, targetable gene fusions in peripheral T-cell lymphoma.

Peripheral T-cell lymphomas (PTCLs) represent a heterogeneous group of T-cell malignancies that generally demonstrate aggressive clinical behavior, often are refractory to standard therapy, and remain significantly understudied. The most common World Health Organization subtype is PTCL, not otherwise specified (NOS), essentially a "wastebasket" category because of inadequate understanding to assign cases to a more specific diagnostic entity. Identification of novel fusion genes has contributed significantly to improving the classification, biologic understanding, and therapeutic targeting of PTCLs. Here, we integrated mate-pair DNA and RNA next-generation sequencing to identify chromosomal rearrangements encoding expressed fusion transcripts in PTCL, NOS. Two of 11 cases had novel fusions involving VAV1, encoding a truncated form of the VAV1 guanine nucleotide exchange factor important in T-cell receptor signaling. Fluorescence in situ hybridization studies identified VAV1 rearrangements in 10 of 148 PTCLs (7%). These were observed exclusively in PTCL, NOS (11%) and anaplastic large cell lymphoma (11%). In vitro, ectopic expression of a VAV1 fusion promoted cell growth and migration in a RAC1-dependent manner. This growth was inhibited by azathioprine, a clinically available RAC1 inhibitor. We also identified novel kinase gene fusions, ITK-FER and IKZF2-ERBB4, as candidate therapeutic targets that show similarities to known recurrent oncogenic ITK-SYK fusions and ERBB4 transcript variants in PTCLs, respectively. Additional novel and potentially clinically relevant fusions also were discovered. Together, these findings identify VAV1 fusions as recurrent and targetable events in PTCLs and highlight the potential for clinical sequencing to guide individualized therapy approaches for this group of aggressive malignancies.

Target-enrichment sequencing and copy number evaluation in inherited polyneuropathy.

OBJECTIVE: To assess the efficiency of target-enrichment next-generation sequencing (NGS) with copy number assessment in inherited neuropathy diagnosis. METHODS: A 197 polyneuropathy gene panel was designed to assess for mutations in 93 patients with inherited or idiopathic neuropathy without known genetic cause. We applied our novel copy number variation algorithm on NGS data, and validated the identified copy number mutations using CytoScan (Affymetrix). Cost and efficacy of this targeted NGS approach was compared to earlier evaluations. RESULTS: Average coverage depth was ∼760× (median = 600, 99.4% > 100×). Among 93 patients, 18 mutations were identified in 17 cases (18%), including 3 copy number mutations: 2 PMP22 duplications and 1 MPZ duplication. The 2 patients with PMP22 duplication presented with bulbar and respiratory involvement and had absent extremity nerve conductions, leading to axonal diagnosis. Average onset age of these 17 patients was 25 years (2-61 years), vs 45 years for those without genetic discovery. Among those with onset age less than 40 years, the diagnostic yield of targeted NGS approach is high (27%) and cost savings is significant (∼20%). However, the cost savings for patients with late onset age and without family history is not demonstrated. CONCLUSIONS: Incorporating copy number analysis in target-enrichment NGS approach improved the efficiency of mutation discovery for chronic, inherited, progressive length-dependent polyneuropathy diagnosis. The new technology is facilitating a simplified genetic diagnostic algorithm utilizing targeted NGS, clinical phenotypes, age at onset, and family history to improve diagnosis efficiency. Our findings prompt a need for updating the current practice parameters and payer guidelines.

Prevalence of Pathogenic Mutations in Cancer Predisposition Genes among Pancreatic Cancer Patients.

The prevalence of germline pathogenic mutations in a comprehensive panel of cancer predisposition genes is not well-defined for patients with pancreatic ductal adenocarcinoma (PDAC). To estimate the frequency of mutations in a panel of 22 cancer predisposition genes, 96 patients unselected for a family history of cancer who were recruited to the Mayo Clinic Pancreatic Cancer patient registry over a 12-month period were screened by next-generation sequencing. Fourteen pathogenic mutations in 13 patients (13.5%) were identified in eight genes: four in ATM, two in BRCA2, CHEK2, and MSH6, and one in BARD1, BRCA1, FANCM, and NBN. These included nine mutations (9.4%) in established pancreatic cancer genes. Three mutations were found in patients with a first-degree relative with PDAC, and 10 mutations were found in patients with first- or second-degree relatives with breast, pancreas, colorectal, ovarian, or endometrial cancers. These results suggest that a substantial proportion of patients with PDAC carry germline mutations in predisposition genes associated with other cancers and that a better understanding of pancreatic cancer risk will depend on evaluation of families with broad constellations of tumors. These findings highlight the need for recommendations governing germline gene-panel testing of patients with pancreatic cancer.

Draft Genome Sequences of Nine Pseudomonas aeruginosa Strains, Including Eight Clinical Isolates.

We report on nine draft genomes of Pseudomonas aeruginosa isolates, assembled using a hybrid paired-end and Nextera mate-pair library approach. Eight are of clinical origin, and one is the ATCC 27853 strain. We also report their multilocus sequence types.

Recessive MYH6 Mutations in Hypoplastic Left Heart With Reduced Ejection Fraction.

BACKGROUND: The molecular underpinnings of hypoplastic left heart are poorly understood. Staged surgical palliation has dramatically improved survival, yet eventual failure of the systemic right ventricle necessitates cardiac transplantation in a subset of patients. We sought to identify genetic determinants of hypoplastic left heart with latent right ventricular dysfunction in individuals with a Fontan circulation. METHODS AND RESULTS: Evaluation of cardiac structure and function by echocardiography in patients with hypoplastic left heart and their first-degree relatives identified 5 individuals with right ventricular ejection fraction ≤40% after Fontan operation. Whole genome sequencing was performed on DNA from 21 family members, filtering for genetic variants with allele frequency <1% predicted to alter protein structure or expression. Secondary family-based filtering for de novo and recessive variants revealed rare inherited missense mutations on both paternal and maternal alleles of MYH6, encoding myosin heavy chain 6, in 2 patients who developed right ventricular dysfunction 3 to 11 years postoperatively. Parents and siblings who were heterozygous carriers had normal echocardiograms. Protein modeling of the 4 highly conserved amino acid substitutions, residing in both head and tail domains, predicted perturbation of protein structure and function. CONCLUSIONS: In contrast to dominant MYH6 mutations with variable penetrance identified in other congenital heart defects and dilated cardiomyopathy, this study reveals compound heterozygosity for recessive MYH6 mutations in patients with hypoplastic left heart and reduced systemic right ventricular ejection fraction. These findings implicate a shared molecular basis for the developmental arrest and latent myopathy of left and right ventricles, respectively.