Next-Generation Sequencing of Connective Tissue Genes in Patients with Classical Ehlers-Danlos Syndrome

Background: Ehlers-Danlos syndrome (EDS) is a common non-inflammatory, congenital connective tissue disorder. Classical type (cEDS) EDS is one of the more common forms, typically caused by mutations in the COL5A1 and COL5A2 genes, though causative mutations in the COL1A1 gene have also been described. Material and methods: The study group included 59 patients of Polish origin, diagnosed with cEDS. The analysis was performed on genomic DNA (gDNA) with NGS technology, using an Illumina sequencer. Thirty-five genes related to connective tissue were investigated. The pathogenicity of the detected variants was assessed by VarSome. Results: The NGS of 35 genes revealed variants within the COL5A1, COL5A2, COL1A1, and COL1A2 genes for 30 of the 59 patients investigated. Our panel detected no sequence variations for the remaining 29 patients. Discussion: Next-generation sequencing, with an appropriate multigene panel, showed great potential to assist in the diagnosis of EDS and other connective tissue disorders. Our data also show that not all causative genes giving rise to cEDS have been elucidated yet.     

VarWalker: Personalized Mutation Network Analysis of Putative Cancer Genes from Next-Generation Sequencing Data

A major challenge in interpreting the large volume of mutation data identified by next-generation sequencing (NGS) is to distinguish driver mutations from neutral passenger mutations to facilitate the identification of targetable genes and new drugs. Current approaches are primarily based on mutation frequencies of single-genes, which lack the power to detect infrequently mutated driver genes and ignore functional interconnection and regulation among cancer genes. We propose a novel mutation network method, VarWalker, to prioritize driver genes in large scale cancer mutation data. VarWalker fits generalized additive models for each sample based on sample-specific mutation profiles and builds on the joint frequency of both mutation genes and their close interactors. These interactors are selected and optimized using the Random Walk with Restart algorithm in a protein-protein interaction network. We applied the method in >300 tumor genomes in two large-scale NGS benchmark datasets: 183 lung adenocarcinoma samples and 121 melanoma samples. In each cancer, we derived a consensus mutation subnetwork containing significantly enriched consensus cancer genes and cancer-related functional pathways. These cancer-specific mutation networks were then validated using independent datasets for each cancer. Importantly, VarWalker prioritizes well-known, infrequently mutated genes, which are shown to interact with highly recurrently mutated genes yet have been ignored by conventional single-gene-based approaches. Utilizing VarWalker, we demonstrated that network-assisted approaches can be effectively adapted to facilitate the detection of cancer driver genes in NGS data.     

Mutation Scanning in Wheat by Exon Capture and Next-Generation Sequencing

Targeted Induced Local Lesions in Genomes (TILLING) is a reverse genetics approach to identify novel sequence variation in genomes, with the aims of investigating gene function and/or developing useful alleles for breeding. Despite recent advances in wheat genomics, most current TILLING methods are low to medium in throughput, being based on PCR amplification of the target genes. We performed a pilot-scale evaluation of TILLING in wheat by next-generation sequencing through exon capture. An oligonucleotide-based enrichment array covering ~2 Mbp of wheat coding sequence was used to carry out exon capture and sequencing on three mutagenised lines of wheat containing previously-identified mutations in the TaGA20ox1 homoeologous genes. After testing different mapping algorithms and settings, candidate SNPs were identified by mapping to the IWGSC wheat Chromosome Survey Sequences. Where sequence data for all three homoeologues were found in the reference, mutant calls were unambiguous; however, where the reference lacked one or two of the homoeologues, captured reads from these genes were mis-mapped to other homoeologues, resulting either in dilution of the variant allele frequency or assignment of mutations to the wrong homoeologue. Competitive PCR assays were used to validate the putative SNPs and estimate cut-off levels for SNP filtering. At least 464 high-confidence SNPs were detected across the three mutagenized lines, including the three known alleles in TaGA20ox1, indicating a mutation rate of ~35 SNPs per Mb, similar to that estimated by PCR-based TILLING. This demonstrates the feasibility of using exon capture for genome re-sequencing as a method of mutation detection in polyploid wheat, but accurate mutation calling will require an improved genomic reference with more comprehensive coverage of homoeologues.     

Extended RAS and BRAF Mutation Analysis Using Next-Generation Sequencing

Somatic mutations in KRAS, NRAS, and BRAF genes are related to resistance to anti-EGFR antibodies in colorectal cancer. We have established an extended RAS and BRAF mutation assay using a next-generation sequencer to analyze these mutations. Multiplexed deep sequencing was performed to detect somatic mutations within KRAS, NRAS, and BRAF, including minor mutated components. We first validated the technical performance of the multiplexed deep sequencing using 10 normal DNA and 20 formalin-fixed, paraffin-embedded (FFPE) tumor samples. To demonstrate the potential clinical utility of our assay, we profiled 100 FFPE tumor samples and 15 plasma samples obtained from colorectal cancer patients. We used a variant calling approach based on a Poisson distribution. The distribution of the mutation-positive population was hypothesized to follow a Poisson distribution, and a mutation-positive status was defined as a value greater than the significance level of the error rate (α = 2 x 10^-5). The cut-off value was determined to be the average error rate plus 7 standard deviations. Mutation analysis of 100 clinical FFPE tumor specimens was performed without any invalid cases. Mutations were detected at a frequency of 59% (59/100). KRAS mutation concordance between this assay and Scorpion-ARMS was 92% (92/100). DNA obtained from 15 plasma samples was also analyzed. KRAS and BRAF mutations were identified in both the plasma and tissue samples of 6 patients. The genetic screening assay using next-generation sequencer was validated for the detection of clinically relevant RAS and BRAF mutations using FFPE and liquid samples.     

MuA-based Molecular Indexing for Rare Mutation Detection by Next-Generation Sequencing

Detection of low-frequency mutations in cancer genomes or other heterogeneous cell populations requires high-fidelity sequencing. Molecular barcoding is one of the key technologies that enables the differentiation of true mutations from errors, which can be caused by sequencing or library preparation processes. However, current approaches where barcodes are introduced via primer extension or adaptor ligation do not utilize the full power of barcoding, due to complicated library preparation workflows and biases. Here we demonstrate the remarkable tolerance of MuA transposase to the presence of multiple replacements in transposon sequence, and explore this unique feature to engineer the MuA transposome complex with randomised nucleotides in 12 transposon positions, which can be introduced as a barcode into the target molecule after transposition event. We applied the approach of Unique MuA-based Molecular Indexing (UMAMI) to assess the power of rare mutation detection by shortgun sequencing on the Illumina platform. Our results show that UMAMI allows detection of rare mutations readily and reliably, and in this paper we report error rate values for the number of thermophilic DNA polymerases measured by using UMAMI.     

High-Throughput Discovery of Mutations in Tef Semi-Dwarfing Genes by Next-Generation Sequencing Analysis

Tef (Eragrostis tef) is a major cereal crop in Ethiopia. Lodging is the primary constraint to increasing productivity in this allotetraploid species, accounting for losses of ∼15–45% in yield each year. As a first step toward identifying semi-dwarf varieties that might have improved lodging resistance, an ∼6× fosmid library was constructed and used to identify both homeologues of the dw3 semi-dwarfing gene of Sorghum bicolor. An EMS mutagenized population, consisting of ∼21,210 tef plants, was planted and leaf materials were collected into 23 superpools. Two dwarfing candidate genes, homeologues of dw3 of sorghum and rht1 of wheat, were sequenced directly from each superpool with 454 technology, and 120 candidate mutations were identified. Out of 10 candidates tested, six independent mutations were validated by Sanger sequencing, including two predicted detrimental mutations in both dw3 homeologues with a potential to improve lodging resistance in tef through further breeding. This study demonstrates that high-throughput sequencing can identify potentially valuable mutations in under-studied plant species like tef and has provided mutant lines that can now be combined and tested in breeding programs for improved lodging resistance.     

CHOPER Filters Enable Rare Mutation Detection in Complex Mutagenesis Populations by Next-Generation Sequencing

Next-generation sequencing (NGS) has revolutionized genetics and enabled the accurate identification of many genetic variants across many genomes. However, detection of biologically important low-frequency variants within genetically heterogeneous populations remains challenging, because they are difficult to distinguish from intrinsic NGS sequencing error rates. Approaches to overcome these limitations are essential to detect rare mutations in large cohorts, virus or microbial populations, mitochondria heteroplasmy, and other heterogeneous mixtures such as tumors. Modifications in library preparation can overcome some of these limitations, but are experimentally challenging and restricted to skilled biologists. This paper describes a novel quality filtering and base pruning pipeline, called Complex Heterogeneous Overlapped Paired-End Reads (CHOPER), designed to detect sequence variants in a complex population with high sequence similarity derived from All-Codon-Scanning (ACS) mutagenesis. A novel fast alignment algorithm, designed for the specified application, has O(n) time complexity. CHOPER was applied to a p53 cancer mutant reactivation study derived from ACS mutagenesis. Relative to error filtering based on Phred quality scores, CHOPER improved accuracy by about 13% while discarding only half as many bases. These results are a step toward extending the power of NGS to the analysis of genetically heterogeneous populations.     

中国人群中16个已知肺癌基因突变的综合性分析

目的:对过去已知的肺癌基因在中国人群中的突变分布进行综合性的分析,指导肺癌的临床治疗。方法:通过文献查阅,挑选出16个已知的肺癌基因。在112例肺癌样本中,对这16个基因进行大样本的靶基因测序并用Sanger测序来验证。同时,对突变在不同亚组中的分布差异进行分析。结果:16个已知肺癌基因突变可评价60.4%肺癌样本。同时,这些基因的在不同的样本亚组中表现出不同的突变特点;通过功能域的分析及蛋白空间结构的模拟,发现9个可能的突变热点既位于蛋白的功能域内又能导致蛋白空间结构或者表面电荷分布的异常。结论:通过靶基因深度测序,全面分析了16个已知肺癌基因在肺癌不同亚组中的突变分布差异,发现并初步验证了9个可能的肺癌突变热点。     

Identification of a Comprehensive Spectrum of Genetic Factors for Hereditary Breast Cancer in a Chinese Population by Next-Generation Sequencing

The genetic etiology of hereditary breast cancer has not been fully elucidated. Although germline mutations of high-penetrance genes such as BRCA1/2 are implicated in development of hereditary breast cancers, at least half of all breast cancer families are not linked to these genes. To identify a comprehensive spectrum of genetic factors for hereditary breast cancer in a Chinese population, we performed an analysis of germline mutations in 2,165 coding exons of 152 genes associated with hereditary cancer using next-generation sequencing (NGS) in 99 breast cancer patients from families of cancer patients regardless of cancer types. Forty-two deleterious germline mutations were identified in 21 genes of 34 patients, including 18 (18.2%) BRCA1 or BRCA2 mutations, 3 (3%) TP53 mutations, 5 (5.1%) DNA mismatch repair gene mutations, 1 (1%) CDH1 mutation, 6 (6.1%) Fanconi anemia pathway gene mutations, and 9 (9.1%) mutations in other genes. Of seven patients who carried mutations in more than one gene, 4 were BRCA1/2 mutation carriers, and their average onset age was much younger than patients with only BRCA1/2 mutations. Almost all identified high-penetrance gene mutations in those families fulfill the typical phenotypes of hereditary cancer syndromes listed in the National Comprehensive Cancer Network (NCCN) guidelines, except two TP53 and three mismatch repair gene mutations. Furthermore, functional studies of MSH3 germline mutations confirmed the association between MSH3 mutation and tumorigenesis, and segregation analysis suggested antagonism between BRCA1 and MSH3. We also identified a lot of low-penetrance gene mutations. Although the clinical significance of those newly identified low-penetrance gene mutations has not been fully appreciated yet, these new findings do provide valuable epidemiological information for the future studies. Together, these findings highlight the importance of genetic testing based on NCCN guidelines and a multi-gene analysis using NGS may be a supplement to traditional genetic counseling.     

Targeted Next-Generation Sequencing Revealed Novel Mutations in Chinese Ataxia Telangiectasia Patients: A Precision Medicine Perspective

Ataxia telangiectasia (AT) is an autosomal recessive disease characterized by progressive cerebellar ataxia, oculocutaneous telangiectasia and immunodeficiency due to mutations in the ATM gene. We performed targeted next-generation sequencing (NGS) on three unrelated patients and identified five disease-causing variants in three probands, including two pairs of heterozygous variants (FAT–1:c.4396C>T/p.R1466X, c.1608-2A>G; FAT–2:c.4412_4413insT/p.L1472Ffs*19, c.8824C>T/p.Q2942X) and one pair of homozygous variants (FAT–3: c.8110T>G/p.C2704G, Hom). With regard to precision medicine for rare genetic diseases, targeted NGS currently enables the rapid and cost-effective identification of causative mutations and is an updated molecular diagnostic tool that merits further optimization. This high-throughput data-based strategy would propel the development of precision diagnostic methods and establish a foundation for precision medicine.     

Mutation Spectrum of Common Deafness-Causing Genes in Patients with Non-Syndromic Deafness in the Xiamen Area,China

In China, approximately 30,000 babies are born with hearing impairment each year. However, the molecular factors causing congenital hearing impairment in the Xiamen area of Fujian province have not been evaluated. To provide accurate genetic testing and counseling in the Xiamen area, we investigated the molecular etiology of non-syndromic deafness in a deaf population from Xiamen. Unrelated students with hearing impairment (n = 155) who attended Xiamen Special Education School in Fujian Province were recruited for this study. Three common deafness-related genes, GJB2, SLC26A4, and mtDNA12SrRNA, were analyzed using all-exon sequencing. GJB2 mutations were detected in 27.1% (42/155) of the entire cohort. The non-syndromic hearing loss (NSHL) hotspot mutations c.109G>A (p.V37I) and c.235delC were found in this population, whereas the Caucasian hotspot mutation c.35delG was not. The allelic frequency of the c.109G>A mutation was 9.03% (28/310), slightly higher than that of c.235delC (8.39%, 26/310), which is the most common GJB2 mutation in most areas of China. The allelic frequency of the c.109G>A mutation was significantly higher in this Xiamen’s deaf population than that in previously reported cohorts (P = 0.00). The SLC26A4 mutations were found in 16.77% (26/155) of this cohort. The most common pathogenic allele was c.IVS7-2A>G (6.13%, 19/310), and the second most common was the c.1079C>T (p.A360V) mutation (1.94%, 6/310) which has rarely been reported as a hotspot mutation in other studies. The mutation rate of mtDNA12SrRNA in this group was 3.87% (6/155), all being the m.A1555G mutation. These findings show the specificity of the common deaf gene-mutation spectrum in this area. According to this study, there were specific hotspot mutations in Xiamen deaf patients. Comprehensive sequencing analysis of the three common deaf genes can help portray the mutation spectrum and develop optimal testing strategies for deaf patients in this area.     

Mitochondrial DNA Variant Discovery and Evaluation in Human Cardiomyopathies through Next-Generation Sequencing

Mutations in mitochondrial DNA (mtDNA) may cause maternally-inherited cardiomyopathy and heart failure. In homoplasmy all mtDNA copies contain the mutation. In heteroplasmy there is a mixture of normal and mutant copies of mtDNA. The clinical phenotype of an affected individual depends on the type of genetic defect and the ratios of mutant and normal mtDNA in affected tissues. We aimed at determining the sensitivity of next-generation sequencing compared to Sanger sequencing for mutation detection in patients with mitochondrial cardiomyopathy. We studied 18 patients with mitochondrial cardiomyopathy and two with suspected mitochondrial disease. We “shotgun” sequenced PCR-amplified mtDNA and multiplexed using a single run on Roche''s 454 Genome Sequencer. By mapping to the reference sequence, we obtained 1,300× average coverage per case and identified high-confidence variants. By comparing these to >400 mtDNA substitution variants detected by Sanger, we found 98% concordance in variant detection. Simulation studies showed that >95% of the homoplasmic variants were detected at a minimum sequence coverage of 20× while heteroplasmic variants required >200× coverage. Several Sanger “misses” were detected by 454 sequencing. These included the novel heteroplasmic 7501T>C in tRNA serine 1 in a patient with sudden cardiac death. These results support a potential role of next-generation sequencing in the discovery of novel mtDNA variants with heteroplasmy below the level reliably detected with Sanger sequencing. We hope that this will assist in the identification of mtDNA mutations and key genetic determinants for cardiomyopathy and mitochondrial disease.     

Next-Generation Sequencing in the Mycology Lab

New state-of-the-art techniques in sequencing offer valuable tools in both detection of mycobiota and in understanding of the molecular mechanisms of resistance against antifungal compounds and virulence. Introduction of new sequencing platform with enhanced capacity and a reduction in costs for sequence analysis provides a potential powerful tool in mycological diagnosis and research. In this review, we summarize the applications of next-generation sequencing techniques in mycology.     

Targeted Next-Generation Sequencing Identifies a Recurrent Mutation in MCPH1 Associating with Hereditary Breast Cancer Susceptibility

Breast cancer is strongly influenced by hereditary risk factors, a majority of which still remain unknown. Here, we performed a targeted next-generation sequencing of 796 genes implicated in DNA repair in 189 Finnish breast cancer cases with indication of hereditary disease susceptibility and focused the analysis on protein truncating mutations. A recurrent heterozygous mutation (c.904_916del, p.Arg304ValfsTer3) was identified in early DNA damage response gene, MCPH1, significantly associating with breast cancer susceptibility both in familial (5/145, 3.4%, P = 0.003, OR 8.3) and unselected cases (16/1150, 1.4%, P = 0.016, OR 3.3). A total of 21 mutation positive families were identified, of which one-third exhibited also brain tumors and/or sarcomas (P = 0.0007). Mutation carriers exhibited significant increase in genomic instability assessed by cytogenetic analysis for spontaneous chromosomal rearrangements in peripheral blood lymphocytes (P = 0.0007), suggesting an effect for MCPH1 haploinsufficiency on cancer susceptibility. Furthermore, 40% of the mutation carrier tumors exhibited loss of the wild-type allele. These findings collectively provide strong evidence for MCHP1 being a novel breast cancer susceptibility gene, which warrants further investigations in other populations.     

Comprehensive Mutation Analysis for Congenital Muscular Dystrophy: A Clinical PCR-Based Enrichment and Next-Generation Sequencing Panel

The congenital muscular dystrophies (CMDs) comprise a heterogeneous group of heritable muscle disorders with often difficult to interpret muscle pathology, making them challenging to diagnose. Serial Sanger sequencing of suspected CMD genes, while the current molecular diagnostic method of choice, can be slow and expensive. A comprehensive panel test for simultaneous screening of mutations in all known CMD-associated genes would be a more effective diagnostic strategy. Thus, the CMDs are a model disorder group for development and validation of next-generation sequencing (NGS) strategies for diagnostic and clinical care applications. Using a highly multiplexed PCR-based target enrichment method (RainDance) in conjunction with NGS, we performed mutation detection in all CMD genes of 26 samples and compared the results with Sanger sequencing. The RainDance NGS panel showed great consistency in coverage depth, on-target efficiency, versatility of mutation detection, and genotype concordance with Sanger sequencing, demonstrating the test''s appropriateness for clinical use. Compared to single tests, a higher diagnostic yield was observed by panel implementation. The panel''s limitation is the amplification failure of select gene-specific exons which require Sanger sequencing for test completion. Successful validation and application of the CMD NGS panel to improve the diagnostic yield in a clinical laboratory was shown.     

Mutation Prevalence of Cerebral Cavernous Malformation Genes in Spanish Patients

Objective

To study the molecular genetic and clinical features of cerebral cavernous malformations (CCM) in a cohort of Spanish patients.

Methods

We analyzed the CCM1, CCM2, and CCM3 genes by MLPA and direct sequencing of exons and intronic boundaries in 94 familial forms and 41 sporadic cases of CCM patients of Spanish extraction. When available, RNA studies were performed seeking for alternative or cryptic splicing.

Results

A total of 26 pathogenic mutations, 22 of which predict truncated proteins, were identified in 29 familial forms and in three sporadic cases. The repertoire includes six novel non-sense and frameshift mutations in CCM1 and CCM3. We also found four missense mutations, one of them located at the third NPXY motif of CCM1 and another one that leads to cryptic splicing of CCM1 exon 6. We found four genomic deletions with the loss of the whole CCM2 gene in one patient and a partial loss of CCM1and CCM2 genes in three other patients. Four families had mutations in CCM3. The results include a high frequency of intronic variants, although most of them localize out of consensus splicing sequences. The main symptoms associated to clinical debut consisted of cerebral haemorrhage, migraines and epileptic seizures. The rare co-occurrence of CCM with Noonan and Chiari syndromes and delayed menarche is reported.

Conclusions

Analysis of CCM genes by sequencing and MLPA has detected mutations in almost 35% of a Spanish cohort (36% of familial cases and 10% of sporadic patients). The results include 13 new mutations of CCM genes and the main clinical symptoms that deserves consideration in molecular diagnosis and genetic counselling of cerebral cavernous malformations.     

基于高通量测序技术的微生物检测数据分析方法

高通量测序技术的发展正在逐渐改变诸多生物学领域的研究方法.为应对突发疫情以及新发未知微生物威胁的需求,微生物鉴定技术逐渐从传统的物理化学方法及核酸杂交等分子水平方法进一步走向利用无需培养的测序数据进行快速分析检测.随之而来的是对高通量数据分析在精度及速度的要求.基于高通量测序数据的微生物检测数据分析方法在近些年得到了快速的发展.本文分析了目前基于高通量测序数据的微生物检测数据分析方法,对其数据分析的处理流程和计算方法进行了研究,比较了各个微生物检测数据分析方法的特点及适用场景.最后结合本实验室工作总结微生物检测数据分析方法在实际应用中可能遇到的问题,希望对该应用领域的研究有一定的参考意义.     

Screening of Rice Genes from the cDNA Catalog Using the Data Obtained by Protein Sequencing

The partial amino acid sequences of 121 rice proteins separated by two-dimensional gel electrophoresis (2D-PAGE), were determined for a protein sequence data file. In the Rice Genome Research Program (RGP), more than 20,000 cDNA clones randomly selected from rice cDNA libraries have been sequenced to construct a cDNA catalog. Complimentary DNAs encoding about 30% of proteins in the protein sequence data file could be identified in the catalog by computer search. It was deduced that 20,000–40,000 genes are present in the rice genome. Only half of about 20,000 cDNAs sequenced in the RGP, corresponding to 1/4–1/2 of genes present in the entire rice genome, should have unique sequences after considering gene redundancy. This is consistent with the fact that the cDNAs encoding about 30% of the sequenced proteins could be identified in the catalog. If the size of the cDNA catalog is enlarged further, cDNAs encoding all proteins separated by 2D-PAGE could be easily identified from the catalog by using the protein sequence data.     

Next-Generation Sequencing of Duplication CNVs Reveals that Most Are Tandem and Some Create Fusion Genes at Breakpoints

Interpreting the genomic and phenotypic consequences of copy-number variation (CNV) is essential to understanding the etiology of genetic disorders. Whereas deletion CNVs lead obviously to haploinsufficiency, duplications might cause disease through triplosensitivity, gene disruption, or gene fusion at breakpoints. The mutational spectrum of duplications has been studied at certain loci, and in some cases these copy-number gains are complex chromosome rearrangements involving triplications and/or inversions. However, the organization of clinically relevant duplications throughout the genome has yet to be investigated on a large scale. Here we fine-mapped 184 germline duplications (14.7 kb–25.3 Mb; median 532 kb) ascertained from individuals referred for diagnostic cytogenetics testing. We performed next-generation sequencing (NGS) and whole-genome sequencing (WGS) to sequence 130 breakpoints from 112 subjects with 119 CNVs and found that most (83%) were tandem duplications in direct orientation. The remainder were triplications embedded within duplications (8.4%), adjacent duplications (4.2%), insertional translocations (2.5%), or other complex rearrangements (1.7%). Moreover, we predicted six in-frame fusion genes at sequenced duplication breakpoints; four gene fusions were formed by tandem duplications, one by two interconnected duplications, and one by duplication inserted at another locus. These unique fusion genes could be related to clinical phenotypes and warrant further study. Although most duplications are positioned head-to-tail adjacent to the original locus, those that are inverted, triplicated, or inserted can disrupt or fuse genes in a manner that might not be predicted by conventional copy-number assays. Therefore, interpreting the genetic consequences of duplication CNVs requires breakpoint-level analysis.