外显子组测序技术在人类疾病研究中的应用

外显子组测序是针对基因组中的蛋白质编码区,靶向富集外显子区域测序,以发现疾病相关遗传变异的技术。该技术近年越来越多地应用于发现人类基因组低频变异、鉴定单基因遗传病致病基因和肿瘤等复杂疾病易感基因研究,成为人类疾病相关变异研究的重要工具。综述了外显子组测序技术的基本原理及其在人类疾病相关基因研究中的应用。     

高通量测序技术在遗传性耳聋研究中的应用及研究进展

耳聋是一种常见的严重出生缺陷,阐明遗传性耳聋的致病机理不仅能够在临床上辅助诊断,为遗传咨询及耳聋预防提供依据,而且能促进人们更深入地了解耳聋的致病机制,开发新的治疗方法。随着基因组研究技术不断创新,以全基因组测序、全外显子组测序、目标区域测序为代表的高通量测序技术在遗传性耳聋研究中已得到广泛应用。本文总结了近5年全外显子组测序和目标区域测序在遗传性耳聋致病基因研究及临床分子诊断中应用及研究进展,希望能够有助于我国临床耳聋基因诊断技术的发展及诊断水平的提升。     

Whole genome sequencing and its applications in medical genetics

Fundamental improvement was made for genome sequencing since the next-generation sequencing (NGS) came out in the 2000s. The newer technologies make use of the power of massively-parallel short-read DNA sequencing, genome alignment and assembly methods to digitally and rapidly search the genomes on a revolutionary scale, which enable large-scale whole genome sequencing (WGS) accessible and practical for researchers. Nowadays, whole genome sequencing is more and more prevalent in detecting the genetics of diseases, studying causative relations with cancers, making genome-level comparative analysis, reconstruction of human population history, and giving clinical implications and instructions. In this review, we first give a typical pipeline of whole genome sequencing, including the lab template preparation, sequencing, genome assembling and quality control, variants calling and annotations. We compare the difference between whole genome and whole exome sequencing (WES), and explore a wide range of applications of whole genome sequencing for both mendelian diseases and complex diseases in medical genetics. We highlight the impact of whole genome sequencing in cancer studies, regulatory variant analysis, predictive medicine and precision medicine, as well as discuss the challenges of the whole genome sequencing.      

Identification of de novo germline mutations and causal genes for sporadic diseases using trio‐based whole‐exome/genome sequencing

Whole‐genome or whole‐exome sequencing (WGS/WES) of the affected proband together with normal parents (trio) is commonly adopted to identify de novo germline mutations (DNMs) underlying sporadic cases of various genetic disorders. However, our current knowledge of the occurrence and functional effects of DNMs remains limited and accurately identifying the disease‐causing DNM from a group of irrelevant DNMs is complicated. Herein, we provide a general‐purpose discussion of important issues related to pathogenic gene identification based on trio‐based WGS/WES data. Specifically, the relevance of DNMs to human sporadic diseases, current knowledge of DNM biogenesis mechanisms, and common strategies or software tools used for DNM detection are reviewed, followed by a discussion of pathogenic gene prioritization. In addition, several key factors that may affect DNM identification accuracy and causal gene prioritization are reviewed. Based on recent major advances, this review both sheds light on how trio‐based WGS/WES technologies can play a significant role in the identification of DNMs and causal genes for sporadic diseases, and also discusses existing challenges.     

GinMicrosatDb: a genome-wide microsatellite markers database for sesame (<Emphasis Type="Italic">Sesamum indicum</Emphasis> L.)

Molecular breeding in sesame is still at infancy due to limited number of microsatellite markers available and the low level of polymorphism exhibited by them. Therefore, whole genome sequencing was used for development of microsatellite markers so as to ensure availability of substantial number of polymorphic markers for use in marker assisted breeding programs. Whole genome sequencing of sesame variety ‘Swetha’ was done using Illumina paired-end sequencing and Roche 454 shotgun sequencing technologies (GCA_000975565.1 in GenBank). ‘GinMicrosatDb’, a genome-wide microsatellite marker database has been developed using the whole genome sequence data of sesame variety ‘Swetha’. The database consists of microsatellites localized on both linkage groups and scaffolds with their genomic co-ordinates. It provides five sets of forward and reverse primers for each of the microsatellite loci along with the flanking sequences, primer GC content, product size and melting temperature etc. The distribution of microsatellites can be viewed and selected through a genome browser as well as through a physical map. The newly identified microsatellite markers are expected to help sesame breeders in developing marker tags for traits of economic importance thereby bringing about greater efficiency in marker-assisted selection programs.     

New insights into the performance of human whole-exome capture platforms

Whole exome sequencing (WES) is increasingly used in research and diagnostics. WES users expect coverage of the entire coding region of known genes as well as sufficient read depth for the covered regions. It is, however, unknown which recent WES platform is most suitable to meet these expectations. We present insights into the performance of the most recent standard exome enrichment platforms from Agilent, NimbleGen and Illumina applied to six different DNA samples by two sequencing vendors per platform. Our results suggest that both Agilent and NimbleGen overall perform better than Illumina and that the high enrichment performance of Agilent is stable among samples and between vendors, whereas NimbleGen is only able to achieve vendor- and sample-specific best exome coverage. Moreover, the recent Agilent platform overall captures more coding exons with sufficient read depth than NimbleGen and Illumina. Due to considerable gaps in effective exome coverage, however, the three platforms cannot capture all known coding exons alone or in combination, requiring improvement. Our data emphasize the importance of evaluation of updated platform versions and suggest that enrichment-free whole genome sequencing can overcome the limitations of WES in sufficiently covering coding exons, especially GC-rich regions, and in characterizing structural variants.     

One Size Doesn't Fit All - RefEditor: Building Personalized Diploid Reference Genome to Improve Read Mapping and Genotype Calling in Next Generation Sequencing Studies

With rapid decline of the sequencing cost, researchers today rush to embrace whole genome sequencing (WGS), or whole exome sequencing (WES) approach as the next powerful tool for relating genetic variants to human diseases and phenotypes. A fundamental step in analyzing WGS and WES data is mapping short sequencing reads back to the reference genome. This is an important issue because incorrectly mapped reads affect the downstream variant discovery, genotype calling and association analysis. Although many read mapping algorithms have been developed, the majority of them uses the universal reference genome and do not take sequence variants into consideration. Given that genetic variants are ubiquitous, it is highly desirable if they can be factored into the read mapping procedure. In this work, we developed a novel strategy that utilizes genotypes obtained a priori to customize the universal haploid reference genome into a personalized diploid reference genome. The new strategy is implemented in a program named RefEditor. When applying RefEditor to real data, we achieved encouraging improvements in read mapping, variant discovery and genotype calling. Compared to standard approaches, RefEditor can significantly increase genotype calling consistency (from 43% to 61% at 4X coverage; from 82% to 92% at 20X coverage) and reduce Mendelian inconsistency across various sequencing depths. Because many WGS and WES studies are conducted on cohorts that have been genotyped using array-based genotyping platforms previously or concurrently, we believe the proposed strategy will be of high value in practice, which can also be applied to the scenario where multiple NGS experiments are conducted on the same cohort. The RefEditor sources are available at https://github.com/superyuan/refeditor.

This is a PLOS Computational Biology Software Article.

     

Exome Sequencing and Functional Validation in Zebrafish Identify GTDC2 Mutations as a Cause of Walker-Warburg Syndrome

Whole-exome sequencing (WES), which analyzes the coding sequence of most annotated genes in the human genome, is an ideal approach to studying fully penetrant autosomal-recessive diseases, and it has been very powerful in identifying disease-causing mutations even when enrollment of affected individuals is limited by reduced survival. In this study, we combined WES with homozygosity analysis of consanguineous pedigrees, which are informative even when a single affected individual is available, to identify genetic mutations responsible for Walker-Warburg syndrome (WWS), a genetically heterogeneous autosomal-recessive disorder that severely affects the development of the brain, eyes, and muscle. Mutations in seven genes are known to cause WWS and explain 50%-60% of cases, but multiple additional genes are expected to be mutated because unexplained cases show suggestive linkage to diverse loci. Using WES in consanguineous WWS-affected families, we found multiple deleterious mutations in GTDC2 (also known as AGO61). GTDC2's predicted role as an uncharacterized glycosyltransferase is consistent with the function of other genes that are known to be mutated in WWS and that are involved in the glycosylation of the transmembrane receptor dystroglycan. Therefore, to explore the role of GTDC2 loss of function during development, we used morpholino-mediated knockdown of its zebrafish ortholog, gtdc2. We found that gtdc2 knockdown in zebrafish replicates all WWS features (hydrocephalus, ocular defects, and muscular dystrophy), strongly suggesting that GTDC2 mutations cause WWS.     

Secondary findings in 421 whole exome-sequenced Chinese children

Background

Variants with known or possible pathogenicity located in genes that are unrelated to primary disease conditions are defined as secondary findings. Secondary findings are not the primary targets of whole exome and genome sequencing (WES/WGS) assay but can be of great practical value in early disease prevention and intervention. The driving force for this study was to investigate the impact of racial difference and disease background on secondary findings. Here, we analyzed secondary findings frequencies in 421 whole exome-sequenced Chinese children who are phenotypically normal or bear congenital heart diseases/juvenile obesity. In total, 421 WES datasets were processed for potential deleterious variant screening. A reference gene list was defined according to the American College of Medical Genetics and Genomics (ACMG) recommendations for reporting secondary findings v2.0 (ACMG SF v2.0). The variant classification was performed according to the evidence-based guidelines recommended by the joint consensus of the ACMG and the Association for Molecular Pathology (AMP).

Results

Among the 421 WES datasets, we identified 11 known/expected pathogenic variants in 12 individuals, accounting for 2.85% of our samples, which is much higher than the reported frequency in a Caucasian population. In conclusion, secondary findings are not so rare in Chinese children, which means that we should pay more attention to the clinical interpretation of sequencing results.

     

Exomsequenzierung zur Identifizierung von Krankheitsgenen für seltene Syndrome

Background

Whole exome sequencing (WES) is the state-of-the-art method for identification of pathogenic mutations in patients with a Mendelian disorder. WES comprehensively covers the coding sequence of the genome and is a fast and cost-effective technique.

Purpose

As most of the technical difficulties have been overcome for WES, the major issue is data processing and analysis to find the pathogenic sequence variation among tens of thousands of sequence changes. Bioinformatic analysis pipelines for filtering sequence variants have to be adapted according to the patients and family members examined by WES and the most likely inheritance pattern underlying the disease.

Possible approaches

Based on 4 cases, different variant prioritization strategies which led to identification of the most likely causative changes in the index patients are described.     

Sequencing and assembly of low copy and genic regions of isolated Triticum aestivum chromosome arm 7DS

The genome of bread wheat (Triticum aestivum) is predicted to be greater than 16 Gbp in size and consist predominantly of repetitive elements, making the sequencing and assembly of this genome a major challenge. We have reduced genome sequence complexity by isolating chromosome arm 7DS and applied second‐generation technology and appropriate algorithmic analysis to sequence and assemble low copy and genic regions of this chromosome arm. The assembly represents approximately 40% of the chromosome arm and all known 7DS genes. Comparison of the 7DS assembly with the sequenced genomes of rice (Oryza sativa) and Brachypodium distachyon identified large regions of conservation. The syntenic relationship between wheat, B. distachyon and O. sativa, along with available genetic mapping data, has been used to produce an annotated draft 7DS syntenic build, which is publicly available at http://www.wheatgenome.info . Our results suggest that the sequencing of isolated chromosome arms can provide valuable information of the gene content of wheat and is a step towards whole‐genome sequencing and variation discovery in this important crop.     

BioNano genome mapping of individual chromosomes supports physical mapping and sequence assembly in complex plant genomes

The assembly of a reference genome sequence of bread wheat is challenging due to its specific features such as the genome size of 17 Gbp, polyploid nature and prevalence of repetitive sequences. BAC‐by‐BAC sequencing based on chromosomal physical maps, adopted by the International Wheat Genome Sequencing Consortium as the key strategy, reduces problems caused by the genome complexity and polyploidy, but the repeat content still hampers the sequence assembly. Availability of a high‐resolution genomic map to guide sequence scaffolding and validate physical map and sequence assemblies would be highly beneficial to obtaining an accurate and complete genome sequence. Here, we chose the short arm of chromosome 7D (7DS) as a model to demonstrate for the first time that it is possible to couple chromosome flow sorting with genome mapping in nanochannel arrays and create a de novo genome map of a wheat chromosome. We constructed a high‐resolution chromosome map composed of 371 contigs with an N50 of 1.3 Mb. Long DNA molecules achieved by our approach facilitated chromosome‐scale analysis of repetitive sequences and revealed a ~800‐kb array of tandem repeats intractable to current DNA sequencing technologies. Anchoring 7DS sequence assemblies obtained by clone‐by‐clone sequencing to the 7DS genome map provided a valuable tool to improve the BAC‐contig physical map and validate sequence assembly on a chromosome‐arm scale. Our results indicate that creating genome maps for the whole wheat genome in a chromosome‐by‐chromosome manner is feasible and that they will be an affordable tool to support the production of improved pseudomolecules.     

Network- and attribute-based classifiers can prioritize genes and pathways for autism spectrum disorders and intellectual disability

Autism spectrum disorders (ASD) are a group of related neurodevelopmental disorders with significant combined prevalence (～1%) and high heritability. Dozens of individually rare genes and loci associated with high-risk for ASD have been identified, which overlap extensively with genes for intellectual disability (ID). However, studies indicate that there may be hundreds of genes that remain to be identified. The advent of inexpensive massively parallel nucleotide sequencing can reveal the genetic underpinnings of heritable complex diseases, including ASD and ID. However, whole exome sequencing (WES) and whole genome sequencing (WGS) provides an embarrassment of riches, where many candidate variants emerge. It has been argued that genetic variation for ASD and ID will cluster in genes involved in distinct pathways and protein complexes. For this reason, computational methods that prioritize candidate genes based on additional functional information such as protein-protein interactions or association with specific canonical or empirical pathways, or other attributes, can be useful. In this study we applied several supervised learning approaches to prioritize ASD or ID disease gene candidates based on curated lists of known ASD and ID disease genes. We implemented two network-based classifiers and one attribute-based classifier to show that we can rank and classify known, and predict new, genes for these neurodevelopmental disorders. We also show that ID and ASD share common pathways that perturb an overlapping synaptic regulatory subnetwork. We also show that features relating to neuronal phenotypes in mouse knockouts can help in classifying neurodevelopmental genes. Our methods can be applied broadly to other diseases helping in prioritizing newly identified genetic variation that emerge from disease gene discovery based on WES and WGS.     

Down syndrome as a genetic model to evaluate the role of oxidative stress and transsulfuration abnormalities in autism spectrum disorder: A 10‐year longitudinal cohort study

Autism spectrum disorder (ASD) is a neurodevelopmental disorder in which evidence reveals oxidative stress and transsulfuration pathway abnormalities. Down syndrome (DS) is a genetic disorder characterized by similar oxidative stress and transsulfuration pathway abnormalities. This hypothesis‐testing longitudinal cohort study determined whether transsulfuration abnormalities and oxidative stress are important susceptibility factors in ASD etiology by evaluating the rate of ASD diagnoses in DS as compared to the general population. The Independent Healthcare Research Database was analyzed for healthcare records prospectively generated in Florida Medicaid. A cohort of 101,736 persons (born: 1990–1999) with ≥10 outpatient office visits and continuously followed for 120 months after birth was examined. There were 942 children in the DS cohort (ICD‐9 code: 758.0) and 100,749 children in the undiagnosed cohort (no DS diagnosis). ASD diagnoses were defined as autistic disorder (ICD‐9 code: 299.00) or Asperger's disorder/pervasive developmental disorder—not otherwise specified (ICD‐9 code: 299.80). ASDs were diagnosed in 5.31% of the DS cohort and 1.34% of the undiagnosed cohort. The risk ratio of being diagnosed with an ASD in the DS cohort as compared to the undiagnosed cohort was 3.97‐fold significantly increased with a risk difference of 3.97%. Among children diagnosed with DS, less than 6% were also diagnosed with an ASD. Among children diagnosed with an ASD, less than 5% were also diagnosed with DS. Children diagnosed with DS are apparently more susceptible to ASD diagnosis relative to the general population suggesting oxidative stress and transsulfuration pathway abnormalities are important susceptibility factors in ASD.     

Chromosome 21-Encoded microRNAs (mRNAs): Impact on Down’s Syndrome and Trisomy-21 Linked Disease

Down’s syndrome (DS; also known as trisomy 21; T21) is caused by a triplication of all or part of human chromosome 21 (chr21). DS is the most common genetic cause of intellectual disability attributable to a naturally-occurring imbalance in gene dosage. DS incurs huge medical, healthcare, and socioeconomic costs, and there are as yet no effective treatments for this incapacitating human neurogenetic disorder. There is a remarkably wide variability in the ‘phenotypic spectrum’ associated with DS; the progression of symptoms and the age of DS onset fluctuate, and there is further variability in the biophysical nature of the chr21 duplication. Besides the cognitive disruptions and dementia in DS patients other serious health problems such as atherosclerosis, altered lipogenesis, Alzheimer’s disease, amyotrophic lateral sclerosis (Lou Gehrig’s disease), autoimmune disease, various cancers including lymphoma, leukemia, glioma and glioblastoma, status epilepticus, congenital heart disease, hypotonia, manic depression, prostate cancer, Usher syndrome, motor disorders, Hirschsprung disease, and various physical anomalies such as early aging occur at elevated frequencies, and all are part of the DS ‘phenotypic spectrum.’ This communication will review the genetic link between these fore-mentioned diseases and a small group of just five stress-associated microRNAs (miRNAs)—that include let-7c, miRNA-99a, miRNA-125b, miRNA-155, and miRNA-802—encoded and clustered on the long arm of human chr21 and spanning the chr21q21.1-chr21q21.3 region.     

Genome-wide genetic diversity of 'Nici', the DNA source for the CHORI-260 turkey BAC library and candidate for whole genome sequencing

Vertebrate whole genome sequence assembly can benefit from a priori knowledge of variability in the target genome, with researchers often selecting highly inbred individuals for sequencing. However, for most species highly inbred research lines are lacking, requiring the use of an outbred individual(s). Here we examined the source DNA [Nicholas inbred (Nici)] of the CHORI-260 turkey bacterial artificial chromosome (BAC) library through analysis of microsatellites and BAC sequences. Heterozygosity of Nici was compared with that of individuals from several breeder lines. Seventy-eight microsatellites were screened for polymorphism in a total of 43 birds, identifying an average individual heterozygosity of 0.39, with Nici at 0.35. Additional loci (total of 147) were examined on a subset of individuals to obtain better genome coverage. The mean heterozygosity for this subset was 0.33 with Nici at 0.31. Examination of approximately 200 kb of genome sequence identified SNPs in the order of one per 200 bp in Nici. These data suggest that the heterozygosity of Nici is comparable to other birds of selected breeder lines and that whole genome sequencing would result in an abundant resource of genome-wide polymorphisms.     

一个非综合征性耳聋家系致病基因的研究

非综合征性耳聋(nonsyndromic hearing impairment, NSHI)是一种十分常见的人类神经系统疾病, 约有1/1000的新生儿患有语前聋。GJB2基因编码间隙连接蛋白Cx26, 是最常见的NSHI致病基因, 大约50%的常染色体隐性遗传NSHI是由GJB2基因突变引起的。在本研究中, 收集了江苏省一个复杂的非综合征性耳聋家系, 并对其进行了分子遗传学研究。对所有已知常染色体隐性遗传的NSHI致病基因, 选用其侧翼的微卫星标记进行连锁分析, 发现该家系的致病基因与D13S175连锁。对GJB2基因进行整个编码区域的测序, 发现235碱基处发生了碱基C的纯合缺失, 这一突变可能是该家系中绝大多数患者致病的遗传基础。     

Whole genome sequencing identifies <Emphasis Type="Italic">ANXA3</Emphasis> and <Emphasis Type="Italic">MTHFR</Emphasis> mutations in a large family with an unknown equinus deformity associated genetic disorder

The aim of this study was to characterize a previously uncharacterized genetic disorder associated with equinus deformity in a large Chinese family at the genetic level. Blood samples were obtained and whole genome sequencing was performed. Differential gene variants were identified and potential impacts on protein structure were predicted. Based on the control sample, several diseases associated variants were identified and selected for further validation. One of the potential variants identified was a ANXA3 gene [chr4, c.C820T(p.R274*)] variant. Further bioinformatic analysis showed that the observed mutation could lead to a three-dimensional conformational change. Moreover, a MTHFR variant that is different from variants associated with clubfoot was also identified. Bioinformatic analysis showed that this mutation could alter the protein binding region. These findings imply that this uncharacterized genetic disorder is not clubfoot, despite sharing some similar symptoms. Furthermore, specific CNV profiles were identified in association with the diseased samples, thus further speaking to the complexity of this multigenerational disorder. This study examined a previously uncharacterized genetic disorder appearing similar to clubfoot and yet having distinct features. Following whole genome sequencing and comparative analysis, several differential gene variants were identified to enable a further distinction from clubfoot. It is hoped that these findings will provide further insight into this disorder and other similar disorders.