Scalable detection of technically challenging variants through modified next-generation sequencing.

BACKGROUND: Some clinically important genetic variants are not easily evaluated with next-generation sequencing (NGS) methods due to technical challenges arising from high- similarity copies (e.g., PMS2, SMN1/SMN2, GBA1, HBA1/HBA2, CYP21A2), repetitive short sequences (e.g., ARX polyalanine repeats, FMR1 AGG interruptions in CGG repeats, CFTR poly-T/TG repeats), and other complexities (e.g., MSH2 Boland inversions). METHODS: We customized our NGS processes to detect the technically challenging variants mentioned above with adaptations including target enrichment and bioinformatic masking of similar sequences. Adaptations were validated with samples of known genotypes. RESULTS: Our adaptations provided high-sensitivity and high-specificity detection for most of the variants and provided a high-sensitivity primary assay to be followed with orthogonal disambiguation for the others. The sensitivity of the NGS adaptations was 100% for all of the technically challenging variants. Specificity was 100% for those in PMS2, GBA1, SMN1/SMN2, and HBA1/HBA2, and for the MSH2 Boland inversion; 97.8%-100% for CYP21A2 variants; and 85.7% for ARX polyalanine repeats. CONCLUSIONS: NGS assays can detect technically challenging variants when chemistries and bioinformatics are jointly refined. The adaptations described support a scalable, cost-effective path to identifying all clinically relevant variants within a single sample.

Correction: Sequencing as a first-line methodology for cystic fibrosis carrier screening.

The original version of this Article contained an error in Figure 3. Specifically, the result "3 (67%) TOP" should read "2 (67%) TOP." This has now been corrected in both the PDF and HTML versions of the Article.

Sequencing as a first-line methodology for cystic fibrosis carrier screening.

PURPOSE: Medical society guidelines recommend offering genotyping-based cystic fibrosis (CF) carrier screening to pregnant women or women considering pregnancy. We assessed the performance of sequencing-based CF screening relative to genotyping, in terms of analytical validity, clinical validity, clinical impact, and clinical utility. METHODS: Analytical validity was assessed using orthogonal confirmation and reference samples. Clinical validity was evaluated using the CFTR2 database. Clinical impact was assessed using ~100,000 screened patients. Three screening strategies were compared: genotyping 23 guideline-recommended variants ("CF23"), sequencing all coding bases in CFTR ("NGS"), and sequencing with large copy-number variant (CNV) identification ("NGS + CNV"). Clinical utility was determined via self-reported actions of at-risk couples (ARCs). RESULTS: Analytical accuracy of NGS + CNV was 100% for SNVs, indels, and CNVs; interpretive clinical specificity relative to CFTR2 was 99.5%. NGS + CNV detected 58 ARCs, 18 of whom would have gone undetected with CF23 alone. Most ARCs (89% screened preconceptionally, 56% prenatally) altered pregnancy management, and no significant differences were observed between ARCs with or without at least one non-CF23 variant. CONCLUSION: Modern NGS and variant interpretation enable accurate sequencing-based CF screening. Limiting screening to 23 variants does not improve analytical validity, clinical validity, or clinical utility, but does fail to detect approximately 30% (18/58) of ARCs.

Detecting clinically actionable variants in the 3' exons of PMS2 via a reflex workflow based on equivalent hybrid capture of the gene and its pseudogene.

BACKGROUND: Hereditary cancer screening (HCS) for germline variants in the 3' exons of PMS2, a mismatch repair gene implicated in Lynch syndrome, is technically challenging due to homology with its pseudogene PMS2CL. Sequences of PMS2 and PMS2CL are so similar that next-generation sequencing (NGS) of short fragments-common practice in multigene HCS panels-may identify the presence of a variant but fail to disambiguate whether its origin is the gene or the pseudogene. Molecular approaches utilizing longer DNA fragments, such as long-range PCR (LR-PCR), can definitively localize variants in PMS2, yet applying such testing to all samples can have logistical and economic drawbacks. METHODS: To address these drawbacks, we propose and characterize a reflex workflow for variant discovery in the 3' exons of PMS2. We cataloged the natural variation in PMS2 and PMS2CL in 707 samples and designed hybrid-capture probes to enrich the gene and pseudogene with equal efficiency. For PMS2 exon 11, NGS reads were aligned, filtered using gene-specific variants, and subject to standard diploid variant calling. For PMS2 exons 12-15, the NGS reads were permissively aligned to PMS2, and variant calling was performed with the expectation of observing four alleles (i.e., tetraploid calling). In this reflex workflow, short-read NGS identifies potentially reportable variants that are then subject to disambiguation via LR-PCR-based testing. RESULTS: Applying short-read NGS screening to 299 HCS samples and cell lines demonstrated >99% analytical sensitivity and >99% analytical specificity for single-nucleotide variants (SNVs) and short insertions and deletions (indels), as well as >96% analytical sensitivity and >99% analytical specificity for copy-number variants. Importantly, 92% of samples had resolved genotypes from short-read NGS alone, with the remaining 8% requiring LR-PCR reflex. CONCLUSION: Our reflex workflow mitigates the challenges of screening in PMS2 and serves as a guide for clinical laboratories performing multigene HCS. To facilitate future exploration and testing of PMS2 variants, we share the raw and processed LR-PCR data from commercially available cell lines, as well as variant frequencies from a diverse patient cohort.

Glucose deprivation elicits phenotypic plasticity via ZEB1-mediated expression of NNMT.

Glucose is considered the primary energy source for all cells, and some cancers are addicted to glucose. Here, we investigated the functional consequences of chronic glucose deprivation in serous ovarian cancer cells. We found that cells resistant to glucose starvation (glucose-restricted cells) demonstrated increased metabolic plasticity that was dependent on NNMT (Nicotinamide N-methyltransferase) expression. We further show that ZEB1 induced NNMT, rendered cells resistant to glucose deprivation and recapitulated metabolic adaptations and mesenchymal gene expression observed in glucose-restricted cells. NNMT depletion reversed metabolic plasticity in glucose-restricted cells and prevented de novo formation of glucose-restricted colonies. In addition to its role in glucose independence, we found that NNMT was required for other ZEB1-induced phenotypes, such as increased migration. NNMT protein levels were also elevated in metastatic and recurrent tumors compared to matched primary carcinomas, while normal ovary and fallopian tube tissue had no detectable NNMT expression. Our studies define a novel ZEB1/NNMT signaling axis, which elicits mesenchymal gene expression, as well as phenotypic and metabolic plasticity in ovarian cancer cells upon chronic glucose starvation. Understanding the causes of cancer cell plasticity is crucial for the development of therapeutic strategies to counter intratumoral heterogeneity, acquired drug resistance and recurrence in high-grade serous ovarian cancer (HGSC).

Development and validation of a 36-gene sequencing assay for hereditary cancer risk assessment.

The past two decades have brought many important advances in our understanding of the hereditary susceptibility to cancer. Numerous studies have provided convincing evidence that identification of germline mutations associated with hereditary cancer syndromes can lead to reductions in morbidity and mortality through targeted risk management options. Additionally, advances in gene sequencing technology now permit the development of multigene hereditary cancer testing panels. Here, we describe the 2016 revision of the Counsyl Inherited Cancer Screen for detecting single-nucleotide variants (SNVs), short insertions and deletions (indels), and copy number variants (CNVs) in 36 genes associated with an elevated risk for breast, ovarian, colorectal, gastric, endometrial, pancreatic, thyroid, prostate, melanoma, and neuroendocrine cancers. To determine test accuracy and reproducibility, we performed a rigorous analytical validation across 341 samples, including 118 cell lines and 223 patient samples. The screen achieved 100% test sensitivity across different mutation types, with high specificity and 100% concordance with conventional Sanger sequencing and multiplex ligation-dependent probe amplification (MLPA). We also demonstrated the screen's high intra-run and inter-run reproducibility and robust performance on blood and saliva specimens. Furthermore, we showed that pathogenic Alu element insertions can be accurately detected by our test. Overall, the validation in our clinical laboratory demonstrated the analytical performance required for collecting and reporting genetic information related to risk of developing hereditary cancers.

Profile of Inflammation-associated genes during Hepatic Differentiation of Human Pluripotent Stem Cells.

Expression of genes associated with inflammation was analyzed during differentiation of human pluripotent stem cells (PSCs) to hepatic cells. Messenger RNA transcript profiles of differentiated endoderm (day 5), hepatoblast (day 15) and hepatocyte-like cells (day 21) were obtained by RNA sequencing analysis. When compared to endoderm cells an immature cell type, the hepatic cells (days 15 and 21) had significantly higher expression of acute phase protein genes including complement factors, coagulation factors, serum amyloid A and serpins. Furthermore, hepatic phase of cells expressed proinflammatory cytokines IL18 and IL32 as well as cytokine receptors IL18R1, IL1R1, IL1RAP, IL2RG, IL6R, IL6ST and IL10RB. These cells also produced CCL14, CCL15, and CXCL- 1, 2, 3, 16 and 17 chemokines. Endoderm cells had higher levels of chemokine receptors, CXCR4 and CXCR7, than that of hepatic cells. Sirtuin family of genes involved in aging, inflammation and metabolism were differentially regulated in endoderm and hepatic phase cells. Ligands and receptors of the tumor necrosis factor (TNF) family as well as downstream signaling factors TRAF2, TRAF4, FADD, NFKB1 and NFKBIB were differentially expressed during hepatic differentiation.

Characterization of type I interferon pathway during hepatic differentiation of human pluripotent stem cells and hepatitis C virus infection.

Pluripotent stem cells are being actively studied as a cell source for regenerating damaged liver. For long-term survival of engrafting cells in the body, not only do the cells have to execute liver-specific function but also withstand the physical strains and invading pathogens. The cellular innate immune system orchestrated by the interferon (IFN) pathway provides the first line of defense against pathogens. The objective of this study is to assess the innate immune function as well as to systematically profile the IFN-induced genes during hepatic differentiation of pluripotent stem cells. To address this objective, we derived endodermal cells (day 5 post-differentiation), hepatoblast (day 15) and hepatocyte-like cells (day 21) from human embryonic stem cells (hESCs). Day 5, 15 and 21 cells were stimulated with IFN-α and subjected to IFN pathway analysis. Transcriptome analysis was carried out by RNA sequencing. The results showed that the IFN-α treatment activated STAT-JAK pathway in differentiating cells. Transcriptome analysis indicated stage specific expression of classical and non-classical IFN-stimulated genes (ISGs). Subsequent validation confirmed the expression of novel ISGs including RASGRP3, CLMP and TRANK1 by differentiated hepatic cells upon IFN treatment. Hepatitis C virus replication in hESC-derived hepatic cells induced the expression of ISGs--LAMP3, ETV7, RASGRP3, and TRANK1. The hESC-derived hepatic cells contain intact innate system and can recognize invading pathogens. Besides assessing the tissue-specific functions for cell therapy applications, it may also be important to test the innate immune function of engrafting cells to ensure adequate defense against infections and improve graft survival.

The development of innate lymphoid cells requires TOX-dependent generation of a common innate lymphoid cell progenitor.

Diverse innate lymphoid cell (ILC) subtypes have been defined on the basis of effector function and transcription factor expression. ILCs derive from common lymphoid progenitors, although the transcriptional pathways that lead to ILC-lineage specification remain poorly characterized. Here we found that the transcriptional regulator TOX was required for the in vivo differentiation of common lymphoid progenitors into ILC lineage-restricted cells. In vitro modeling demonstrated that TOX deficiency resulted in early defects in the survival or proliferation of progenitor cells, as well as ILC differentiation at a later stage. In addition, comparative transcriptome analysis of bone marrow progenitors revealed that TOX-deficient cells failed to upregulate many genes of the ILC program, including genes that are targets of Notch, which indicated that TOX is a key determinant of early specification to the ILC lineage.

Cyclin E1 and RTK/RAS signaling drive CDK inhibitor resistance via activation of E2F and ETS.

High-grade serous ovarian cancers (HGSOC) are genomically complex, heterogeneous cancers with a high mortality rate, due to acquired chemoresistance and lack of targeted therapy options. Cyclin-dependent kinase inhibitors (CDKi) target the retinoblastoma (RB) signaling network, and have been successfully incorporated into treatment regimens for breast and other cancers. Here, we have compared mechanisms of response and resistance to three CDKi that target either CDK4/6 or CDK2 and abrogate E2F target gene expression. We identify CCNE1 gain and RB1 loss as mechanisms of resistance to CDK4/6 inhibition, whereas receptor tyrosine kinase (RTK) and RAS signaling is associated with CDK2 inhibitor resistance. Mechanistically, we show that ETS factors are mediators of RTK/RAS signaling that cooperate with E2F in cell cycle progression. Consequently, CDK2 inhibition sensitizes cyclin E1-driven but not RAS-driven ovarian cancer cells to platinum-based chemotherapy. In summary, this study outlines a rational approach for incorporating CDKi into treatment regimens for HGSOC.

Differentiation of human limbal-derived induced pluripotent stem cells into limbal-like epithelium.

Limbal epithelial stem cell (LESC) deficiency (LSCD) leads to corneal abnormalities resulting in compromised vision and blindness. LSCD can be potentially treated by transplantation of appropriate cells, which should be easily expandable and bankable. Induced pluripotent stem cells (iPSCs) are a promising source of transplantable LESCs. The purpose of this study was to generate human iPSCs and direct them to limbal differentiation by maintaining them on natural substrata mimicking the native LESC niche, including feederless denuded human amniotic membrane (HAM) and de-epithelialized corneas. These iPSCs were generated with nonintegrating vectors from human primary limbal epithelial cells. This choice of parent cells was supposed to enhance limbal cell differentiation from iPSCs by partial retention of parental epigenetic signatures in iPSCs. When the gene methylation patterns were compared in iPSCs to parental LESCs using Illumina global methylation arrays, limbal-derived iPSCs had fewer unique methylation changes than fibroblast-derived iPSCs, suggesting retention of epigenetic memory during reprogramming. Limbal iPSCs cultured for 2 weeks on HAM developed markedly higher expression of putative LESC markers ABCG2, ΔNp63α, keratins 14, 15, and 17, N-cadherin, and TrkA than did fibroblast iPSCs. On HAM culture, the methylation profiles of select limbal iPSC genes (including NTRK1, coding for TrkA protein) became closer to the parental cells, but fibroblast iPSCs remained closer to parental fibroblasts. On denuded air-lifted corneas, limbal iPSCs even upregulated differentiated corneal keratins 3 and 12. These data emphasize the importance of the natural niche and limbal tissue of origin in generating iPSCs as a LESC source with translational potential for LSCD treatment.