A combination of targeted enrichment methodologies for whole-exome sequencing reveals novel pathogenic mutations

Biallelic structural variants in three patients with ERCC8-related Cockayne syndrome and a potential pitfall of copy number variation analysis

Cockayne syndrome (CS) is a rare autosomal recessive disorder caused by mutations in ERCC8 or ERCC6. Most pathogenic variants in ERCC8 are single nucleotide substitutions. Structural variants (SVs) have been reported in patients with ERCC8-related CS. However, comprehensive molecular detection, including SVs of ERCC8, in CS patients remains problematic. Herein, we present three Japanese patients with ERCC8-related CS in whom causative SVs were identified using whole-exome-based copy number variation (CNV) detection tools. One patient showed compound heterozygosity for a 259-kb deletion and a deletion of exon 4 which has previously been reported as an Asia-specific variant. The other two patients were homozygous for the same exon 4 deletion. The exon 4 deletion was detected only by the ExomeDepth software. Intrigued by the discrepancy in the detection capability of various tools for the SVs, we evaluated the analytic performance of four whole-exome-based CNV detection tools using an exome data set from 337 healthy individuals. A total of 1,278,141 exons were predicted as being affected by the 4 CNV tools. Interestingly 95.1% of these affected exons were detected by one tool alone. Thus, we expect that the use of multiple tools may improve the detection rate of SVs from aligned exome data.

Clinical diagnosis of genetic disorders at both single-nucleotide and chromosomal levels based on BGISEQ-500 platform

Most variations in the human genome refer to single-nucleotide variation (SNV), small fragment insertions and deletions, and genomic copy number variation (CNV). Many human diseases including genetic disorders are associated with variations in the genome. These disorders are often difficult to be diagnosed because of their complex clinical conditions, therefore, an effective detection method is needed to facilitate clinical diagnosis and prevent birth defects. With the development of high-throughput sequencing technology, the method of targeted sequence capture chip has been extensively used owing to its high throughput, high accuracy, fast speed, and low cost. In this study, we designed a chip that potentially captured the coding region of 3043 genes associated with 4013 monogenic diseases, with an addition of 148 chromosomal abnormalities that can be identified by targeting specific regions. To assess the efficiency, a strategy of combining the BGISEQ500 sequencing platform with the designed chip was utilized to screen variants in 63 patients. Eventually, 67 disease-associated variants were found, 31 of which were novel. The results of the evaluation test also show that this combined strategy complies with the requirements of clinical testing and has proper clinical application value.

Phylogenomics and evolutionary history of Oreobates (Anura: Craugastoridae) Neotropical frogs along elevational gradients

Mountain ranges offer opportunities for understanding how species evolved and diversified across different environmental conditions. Neotropical frogs of the genus Oreobates (Anura: Craugastoridae) are adapted to highland and lowland habitats along the Andes, but many aspects of their evolution remain unknown. We studied their evolutionary history using ~18,000 exons enriched by targeted sequence-capture. Since capture success was very variable across samples, we evaluated to what degree differing data filtering produced robust inferences. The inferred evolutionary framework evidenced phylogenetic discordances among lowland species that can be explained by taxonomic misidentification or admixture of ancestral lineages. Highland species showed smaller effective populations than lowland frogs, probably due to greater habitat fragmentation in montane environments. Stronger genetic drift likely decreased the power of purifying selection and led to an increased proportion of nonsynonymous mutations in highland populations that could play an important role in their adaptation. Overall, our work sheds light on the evolutionary history and diversification of this group of Neotropical frogs along elevational gradients in the Andes as well as on their patterns of intraspecific diversity.

Whole exome sequencing of a family revealed a novel variant in the CHM gene, c.22delG p.(Glu8Serfs*4), which co-segregated with choroideremia

Choroideremia is a complex form of blindness-causing retinal degeneration. The aim of the present study was to investigate the pathogenic variant and molecular etiology associated with choroideremia in a Chinese family. All available family members underwent detailed ophthalmological examinations. Whole exome sequencing, bioinformatics analysis, Sanger sequencing, and co-segregation analysis of family members were used to validate sequencing data and confirm the presence of the disease-causing gene variant. The proband was diagnosed with choroideremia on the basis of clinical manifestations. Whole exome sequencing showed that the proband had a hemizygous variant in the CHM gene, c.22delG p. (Glu8Serfs*4), which was confirmed by Sanger sequencing and found to co-segregate with choroideremia. The variant was classified as likely pathogenic and has not previously been described. These results expand the spectrum of variants in the CHM gene, thus potentially enriching the understanding of the molecular basis of choroideremia. Moreover, they may provide insight for future choroideremia diagnosis and gene therapy.

IMSindel: An accurate intermediate-size indel detection tool incorporating de novo assembly and gapped global-local alignment with split read analysis

Insertions and deletions (indels) have been implicated in dozens of human diseases through the radical alteration of gene function by short frameshift indels as well as long indels. However, the accurate detection of these indels from next-generation sequencing data is still challenging. This is particularly true for intermediate-size indels (≥50 bp), due to the short DNA sequencing reads. Here, we developed a new method that predicts intermediate-size indels using BWA soft-clipped fragments (unmatched fragments in partially mapped reads) and unmapped reads. We report the performance comparison of our method, GATK, PINDEL and ScanIndel, using whole exome sequencing data from the same samples. False positive and false negative counts were determined through Sanger sequencing of all predicted indels across these four methods. The harmonic mean of the recall and precision, F-measure, was used to measure the performance of each method. Our method achieved the highest F-measure of 0.84 in one sample, compared to 0.56 for GATK, 0.52 for PINDEL and 0.46 for ScanIndel. Similar results were obtained in additional samples, demonstrating that our method was superior to the other methods for detecting intermediate-size indels. We believe that this methodology will contribute to the discovery of intermediate-size indels associated with human disease.

Targeted Gene Panel Sequencing for Early-onset Inflammatory Bowel Disease and Chronic Diarrhea

BACKGROUND

In contrast to adult-onset inflammatory bowel disease (IBD), where many genetic loci have been shown to be involved in complex disease etiology, early-onset IBD (eoIBD) and associated syndromes can sometimes present as monogenic conditions. As a result, the clinical phenotype and ideal disease management in these patients often differ from those in adult-onset IBD. However, due to high costs and the complexity of data analysis, high-throughput screening for genetic causes has not yet become a standard part of the diagnostic work-up of eoIBD patients.

METHODS

We selected 28 genes of interest associated with monogenic IBD and performed targeted panel sequencing in 71 patients diagnosed with eoIBD or early-onset chronic diarrhea to detect causative variants. We compared these results to whole-exome sequencing (WES) data available for 25 of these patients.

RESULTS

Target coverage was significantly higher in the targeted gene panel approach compared with WES, whereas the cost of the panel was considerably lower (approximately 25% of WES). Disease-causing variants affecting protein function were identified in 5 patients (7%), located in genes of the IL10 signaling pathway (3), WAS (1), and DKC1 (1). The functional effects of 8 candidate variants in 5 additional patients (7%) are under further investigation. WES did not identify additional causative mutations in 25 patients.

CONCLUSIONS

Targeted gene panel sequencing is a fast and effective screening method for monogenic causes of eoIBD that should be routinely established in national referral centers.

Deactivated CRISPR Associated Protein 9 for Minor-Allele Enrichment in Cell-Free DNA

BACKGROUND

Cell-free DNA (cfDNA) diagnostics are emerging as a new paradigm of disease monitoring and therapy management. The clinical utility of these diagnostics is relatively limited by a low signal-to-noise ratio, such as with low allele frequency (AF) mutations in cancer. While enriching for rare alleles to increase their AF before sample analysis is one strategy that can greatly improve detection capability, current methods are limited in their generalizability, ease of use, and applicability to point mutations.

METHODS

Leveraging the robust single-base-pair specificity and generalizability of the CRISPR associated protein 9 (Cas9) system, we developed a deactivated Cas9 (dCas9)-based method of minor-allele enrichment capable of efficient single-target and multiplexed enrichment. The dCas9 protein was complexed with single guide RNAs targeted to mutations of interest and incubated with cfDNA samples containing mutant strands at low abundance. Mutation-bound dCas9 complexes were isolated, dissociated, and the captured DNA purified for downstream use.

RESULTS

Targeting the 3 most common epidermal growth factor receptor mutations (exon 19 deletion, T790M, L858R) found in non-small cell lung cancer (NSCLC), we achieved >20-fold increases in AF and detected mutations by use of qPCR at an AF of 0.1%. In a cohort of 18 NSCLC patient-derived cfDNA samples, our method enabled detection of 8 out of 13 mutations that were otherwise undetected by qPCR.

CONCLUSIONS

The dCas9 method provides an important application of the CRISPR/Cas9 system outside the realm of genome editing and can provide a step forward for the detection capability of cfDNA diagnostics.

Targeted exome sequencing of Korean triple-negative breast cancer reveals homozygous deletions associated with poor prognosis of adjuvant chemotherapy-treated patients

Triple-negative breast cancer is characterized by the absence of estrogen and progesterone receptors and human epidermal growth factor receptor 2, and is associated with a poorer outcome than other subtypes of breast cancer. Moreover, there are no accurate prognostic genes or effective therapeutic targets, thereby necessitating continued intensive investigation. This study analyzed the genetic mutation landscape in 70 patients with triple-negative breast cancer by targeted exome sequencing of tumor and matched normal samples. Sequencing showed that more than 50% of these patients had deleterious mutations and homozygous deletions of DNA repair genes, such as ATM, BRCA1, BRCA2, WRN, and CHEK2. These findings suggested that a large number of patients with triple-negative breast cancer have impaired DNA repair function and that therefore a poly ADP-ribose polymerase inhibitor may be an effective drug in the treatment of this disease. Notably, homozygous deletion of three genes, EPHA5, MITF, and ACSL3, was significantly associated with an increased risk of recurrence or distant metastasis in adjuvant chemotherapy-treated patients.

Role of a heterotrimeric G-protein, Gi2, in the corticogenesis: possible involvement in periventricular nodular heterotopia and intellectual disability

We analyzed the role of a heterotrimeric G-protein, Gi2, in the development of the cerebral cortex. Acute knockdown of the α-subunit (Gαi2) with in utero electroporation caused delayed radial migration of excitatory neurons during corticogenesis, perhaps because of impaired morphology. The migration phenotype was rescued by an RNAi-resistant version of Gαi2. On the other hand, silencing of Gαi2 did not affect axon elongation, dendritic arbor formation or neurogenesis at ventricular zone in vivo. When behavior analyses were conducted with acute Gαi2-knockdown mice, they showed defects in social interaction, novelty recognition and active avoidance learning as well as increased anxiety. Subsequently, using whole-exome sequencing analysis, we identified a de novo heterozygous missense mutation (c.680C>T; p.Ala227Val) in the GNAI2 gene encoding Gαi2 in an individual with periventricular nodular heterotopia and intellectual disability. Collectively, the phenotypes in the knockdown experiments suggest a role of Gαi2 in the brain development, and impairment of its function might cause defects in neuronal functions which lead to neurodevelopmental disorders.

Diagnostics based on nucleic acid sequence variant profiling: PCR, hybridization, and NGS approaches

Nucleic acid sequence variations have been implicated in many diseases, and reliable detection and quantitation of DNA/RNA biomarkers can inform effective therapeutic action, enabling precision medicine. Nucleic acid analysis technologies being translated into the clinic can broadly be classified into hybridization, PCR, and sequencing, as well as their combinations. Here we review the molecular mechanisms of popular commercial assays, and their progress in translation into in vitro diagnostics.

Advances in clinical next-generation sequencing: target enrichment and sequencing technologies

The huge parallel sequencing capabilities of next generation sequencing technologies have made them the tools of choice to characterize genomic aberrations for research and diagnostic purposes. For clinical applications, screening the whole genome or exome is challenging owing to the large genomic area to be sequenced, associated costs, complexity of data, and lack of known clinical significance of all genes. Consequently, routine screening involves limited markers with established clinical relevance. This process, referred to as targeted genome sequencing, requires selective enrichment of the genomic areas comprising these markers via one of several primer or probe-based enrichment strategies, followed by sequencing of the enriched genomic areas. Here, the authors review current target enrichment approaches and next generation sequencing platforms, focusing on the underlying principles, capabilities, and limitations of each technology along with validation and implementation for clinical testing.

What Is the Best NGS Enrichment Method for the Molecular Diagnosis of Monogenic Diabetes and Obesity?

Molecular diagnosis of monogenic diabetes and obesity is of paramount importance for both the patient and society, as it can result in personalized medicine associated with a better life and it eventually saves health care spending. Genetic clinical laboratories are currently switching from Sanger sequencing to next-generation sequencing (NGS) approaches but choosing the optimal protocols is not easy. Here, we compared the sequencing coverage of 43 genes involved in monogenic forms of diabetes and obesity, and variant detection rates, resulting from four enrichment methods based on the sonication of DNA (Agilent SureSelect, RainDance technologies), or using enzymes for DNA fragmentation (Illumina Nextera, Agilent HaloPlex). We analyzed coding exons and untranslated regions of the 43 genes involved in monogenic diabetes and obesity. We found that none of the methods achieves yet full sequencing of the gene targets. Nonetheless, the RainDance, SureSelect and HaloPlex enrichment methods led to the best sequencing coverage of the targets; while the Nextera method resulted in the poorest sequencing coverage. Although the sequencing coverage was high, we unexpectedly found that the HaloPlex method missed 20% of variants detected by the three other methods and Nextera missed 10%. The question of which NGS technique for genetic diagnosis yields the highest diagnosis rate is frequently discussed in the literature and the response is still unclear. Here, we showed that the RainDance enrichment method as well as SureSelect, which are both based on the sonication of DNA, resulted in a good sequencing quality and variant detection, while the use of enzymes to fragment DNA (HaloPlex or Nextera) might not be the best strategy to get an accurate sequencing.