Detection of Clinically Relevant Genetic Variants in Autism Spectrum Disorder by Whole-Genome Sequencing

Autism Spectrum Disorder (ASD) demonstrates high heritability and familial clustering, yet the genetic causes remain only partially understood as a result of extensive clinical and genomic heterogeneity. Whole-genome sequencing (WGS) shows promise as a tool for identifying ASD risk genes as well as unreported mutations in known loci, but an assessment of its full utility in an ASD group has not been performed. We used WGS to examine 32 families with ASD to detect de novo or rare inherited genetic variants predicted to be deleterious (loss-of-function and damaging missense mutations). Among ASD probands, we identified deleterious de novo mutations in six of 32 (19%) families and X-linked or autosomal inherited alterations in ten of 32 (31%) families (some had combinations of mutations). The proportion of families identified with such putative mutations was larger than has been previously reported; this yield was in part due to the comprehensive and uniform coverage afforded by WGS. Deleterious variants were found in four unrecognized, nine known, and eight candidate ASD risk genes. Examples include CAPRIN1 and AFF2 (both linked to FMR1, which is involved in fragile X syndrome), VIP (involved in social-cognitive deficits), and other genes such as SCN2A and KCNQ2 (linked to epilepsy), NRXN1, and CHD7, which causes ASD-associated CHARGE syndrome. Taken together, these results suggest that WGS and thorough bioinformatic analyses for de novo and rare inherited mutations will improve the detection of genetic variants likely to be associated with ASD or its accompanying clinical symptoms.     

Klinik und Genetik der spastischen Spinalparalysen

Even before the advent of next generation sequencing (NGS), multiple loci for hereditary spastic paraplegias (HSPs) had already been identified. In the last 2 years, dozens of new disease genes have been added, accounting for a total of 52 established HSP loci and 35 known HSP disease genes. With overlapping phenotypes for distinct genetic entities, the clinical diagnosis is often demanding and high-throughput genetic testing has to parallel a diagnostic workflow. Notwithstanding this aspect, spastin (SPAST) mutations evidently constitute the most important genetic cause in autosomal dominant paraplegia 4 (SPG4). Recently, large studies established that SPAST mutations are even causative in roughly 10?% of sporadic HSP patients. In this review, we suggest a diagnostic routine for HSP and elaborate on how detailed phenotyping and extensive genotyping will assist in the diagnosis of many more HSP subtypes. This ultimately will set a basis for selective clinical observations and therapy development.     

Non-Syndromic Hearing Impairment in India: High Allelic Heterogeneity among Mutations in TMPRSS3, TMC1, USHIC,CDH23 and TMIE

Mutations in the autosomal genes TMPRSS3, TMC1, USHIC, CDH23 and TMIE are known to cause hereditary hearing loss. To study the contribution of these genes to autosomal recessive, non-syndromic hearing loss (ARNSHL) in India, we examined 374 families with the disorder to identify potential mutations. We found four mutations in TMPRSS3, eight in TMC1, ten in USHIC, eight in CDH23 and three in TMIE. Of the 33 potentially pathogenic variants identified in these genes, 23 were new and the remaining have been previously reported. Collectively, mutations in these five genes contribute to about one-tenth of ARNSHL among the families examined. New mutations detected in this study extend the allelic heterogeneity of the genes and provide several additional variants for structure-function correlation studies. These findings have implications for early DNA-based detection of deafness and genetic counseling of affected families in the Indian subcontinent.     

Germline Variation in Cancer-Susceptibility Genes in a Healthy,Ancestrally Diverse Cohort: Implications for Individual Genome Sequencing

Technological advances coupled with decreasing costs are bringing whole genome and whole exome sequencing closer to routine clinical use. One of the hurdles to clinical implementation is the high number of variants of unknown significance. For cancer-susceptibility genes, the difficulty in interpreting the clinical relevance of the genomic variants is compounded by the fact that most of what is known about these variants comes from the study of highly selected populations, such as cancer patients or individuals with a family history of cancer. The genetic variation in known cancer-susceptibility genes in the general population has not been well characterized to date. To address this gap, we profiled the nonsynonymous genomic variation in 158 genes causally implicated in carcinogenesis using high-quality whole genome sequences from an ancestrally diverse cohort of 681 healthy individuals. We found that all individuals carry multiple variants that may impact cancer susceptibility, with an average of 68 variants per individual. Of the 2,688 allelic variants identified within the cohort, most are very rare, with 75% found in only 1 or 2 individuals in our population. Allele frequencies vary between ancestral groups, and there are 21 variants for which the minor allele in one population is the major allele in another. Detailed analysis of a selected subset of 5 clinically important cancer genes, BRCA1, BRCA2, KRAS, TP53, and PTEN, highlights differences between germline variants and reported somatic mutations. The dataset can serve a resource of genetic variation in cancer-susceptibility genes in 6 ancestry groups, an important foundation for the interpretation of cancer risk from personal genome sequences.     

Comprehensive Analysis of Deafness Genes in Families with Autosomal Recessive Nonsyndromic Hearing Loss

Comprehensive genetic testing has the potential to become the standard of care for individuals with hearing loss. In this study, we investigated the genetic etiology of autosomal recessive nonsyndromic hearing loss (ARNSHL) in a Turkish cohort including individuals with cochlear implant, who had a pedigree suggestive of an autosomal recessive inheritance. A workflow including prescreening of GJB2 and a targeted next generation sequencing panel (Illumına TruSight^TM Exome) covering 2761 genes that we briefly called as mendelian exome sequencing was used. This panel includes 102 deafness genes and a number of genes causing Mendelian disorders. Using this approach, we identified causative variants in 21 of 29 families. Three different GJB2 variants were present in seven families. Remaining 14 families had 15 different variants in other known NSHL genes (MYO7A, MYO15A, MARVELD2, TMIE, DFNB31, LOXHD1, GPSM2, TMC1, USH1G, CDH23). Of these variants, eight are novel. Mutation detection rate of our workflow is 72.4%, confirming the usefulness of targeted sequencing approach in NSHL.     

Mutations of 60 known causative genes in 157 families with retinitis pigmentosa based on exome sequencing

Retinitis pigmentosa (RP) is the most common and highly heterogeneous form of hereditary retinal degeneration. This study was to identify mutations in the 60 genes that were known to be associated with RP in 157 unrelated Chinese families with RP. Genomic DNA from probands was initially analyzed by whole exome sequencing. Sanger sequencing was used to confirm potential candidate variants affecting the encoded residues in the 60 genes, including heterozygous variants from genes that are related to autosomal dominant RP, homozygous or compound heterozygous variants from genes that are related to autosomal recessive RP, and hemizygous variants from genes that are related to X-linked RP. Synonymous and intronic variants were also examined to confirm whether they could affect splicing. A total of 244 candidate variants were detected by exome sequencing. Sanger sequencing confirmed 240 variants out of the 244 candidates. Informatics and segregation analyses suggested 110 potential pathogenic mutations in 28 out of the 60 genes involving 79 of the 157 (50 %) families, including 31 (39 %, 31/79) families with heterozygous mutations in autosomal dominant genes, 37 (47 %, 37/79) families with homozygous (9) or compound heterozygous (28) mutations in autosomal recessive genes, and 11 (14 %, 11/79) families with hemizygous mutations in X-linked genes. Of the 110 identified variants, 74 (67 %) were novel. The genetic defects in approximately half of the 157 studies families were detected by exome sequencing. A comprehensive analysis of the 60 known genes not only expanded the mutation spectrum and frequency of the 60 genes in Chinese patients with RP, but also provided an overview of the molecular etiology of RP in Chinese patients. The analysis of the known genes also supplied the foundation and clues for discovering novel causative RP genes.     

Whole-Exome Sequencing Efficiently Detects Rare Mutations in Autosomal Recessive Nonsyndromic Hearing Loss

Identification of the pathogenic mutations underlying autosomal recessive nonsyndromic hearing loss (ARNSHL) is difficult, since causative mutations in 39 different genes have so far been reported. After excluding mutations in the most common ARNSHL gene, GJB2, via Sanger sequencing, we performed whole-exome sequencing (WES) in 30 individuals from 20 unrelated multiplex consanguineous families with ARNSHL. Agilent SureSelect Human All Exon 50 Mb kits and an Illumina Hiseq2000 instrument were used. An average of 93%, 84% and 73% of bases were covered to 1X, 10X and 20X within the ARNSHL-related coding RefSeq exons, respectively. Uncovered regions with WES included those that are not targeted by the exome capture kit and regions with high GC content. Twelve homozygous mutations in known deafness genes, of which eight are novel, were identified in 12 families: MYO15A-p.Q1425X, -p.S1481P, -p.A1551D; LOXHD1-p.R1494X, -p.E955X; GIPC3-p.H170N; ILDR1-p.Q274X; MYO7A-p.G2163S; TECTA-p.Y1737C; TMC1-p.S530X; TMPRSS3-p.F13Lfs*10; TRIOBP-p.R785Sfs*50. Each mutation was within a homozygous run documented via WES. Sanger sequencing confirmed co-segregation of the mutation with deafness in each family. Four rare heterozygous variants, predicted to be pathogenic, in known deafness genes were detected in 12 families where homozygous causative variants were already identified. Six heterozygous variants that had similar characteristics to those abovementioned variants were present in 15 ethnically-matched individuals with normal hearing. Our results show that rare causative mutations in known ARNSHL genes can be reliably identified via WES. The excess of heterozygous variants should be considered during search for causative mutations in ARNSHL genes, especially in small-sized families.     

Evaluation of the Genetic Basis of Familial Aggregation of Pacemaker Implantation by a Large Next Generation Sequencing Panel

Background

The etiology of conduction disturbances necessitating permanent pacemaker (PPM) implantation is often unknown, although familial aggregation of PPM (faPPM) suggests a possible genetic basis. We developed a pan-cardiovascular next generation sequencing (NGS) panel to genetically characterize a selected cohort of faPPM.

Materials and Methods

We designed and validated a custom NGS panel targeting the coding and splicing regions of 246 genes with involvement in cardiac pathogenicity. We enrolled 112 PPM patients and selected nine (8%) with faPPM to be analyzed by NGS.

Results

Our NGS panel covers 95% of the intended target with an average of 229x read depth at a minimum of 15-fold depth, reaching a SNP true positive rate of 98%. The faPPM patients presented with isolated cardiac conduction disease (ICCD) or sick sinus syndrome (SSS) without overt structural heart disease or identifiable secondary etiology. Three patients (33.3%) had heterozygous deleterious variants previously reported in autosomal dominant cardiac diseases including CCD: LDB3 (p.D117N) and TRPM4 (p.G844D) variants in patient 4; TRPM4 (p.G844D) and ABCC9 (p.V734I) variants in patient 6; and SCN5A (p.T220I) and APOB (p.R3527Q) variants in patient 7.

Conclusion

FaPPM occurred in 8% of our PPM clinic population. The employment of massive parallel sequencing for a large selected panel of cardiovascular genes identified a high percentage (33.3%) of the faPPM patients with deleterious variants previously reported in autosomal dominant cardiac diseases, suggesting that genetic variants may play a role in faPPM.     

Morbus Parkinson

Monogenic Parkinson’s disease (PD), i.?e. parkinsonism caused by mutations in single genes, represents ~5% of all PD cases. Over the past 20 years, three autosomal dominantly (SNCA, LRRK2, VPS35) and three autosomal recessively (Parkin, PINK1, DJ-1) inherited causal PD genes have been identified and validated. Although pathogenic changes in SNCA are very rare, begin early, and may be associated with the development of dementia, pathogenic variants in LRRK2-linked PD are most common among monogenic PS and patients are clinically indistinguishable from those with idiopathic PD. In patients with onset before the age of 40 years, pathogenic variants in the Parkin and PINK1 genes should be suspected, and in patients with a positive family history, genetic counseling should be carried out. Recently, dynamic developments in the area of Parkinson’s genetics have led to new therapeutic approaches and the first gene-specific therapies have entered the early testing phase. Besides the established monogenic PD genes, candidate genes have been identified, but not yet conclusively validated. In addition to established monogenic PD, as yet unvalidated Parkinson’s candidate genes and well-characterized genetic risk exist at this time. As monogenic PD represents a “model disease” for idiopathic PD too, further progress toward more personalized medicine may be expected for both monogenic and idiopathic PD.     

The Characteristics of Heterozygous Protein Truncating Variants in the Human Genome

Sequencing projects have identified large numbers of rare stop-gain and frameshift variants in the human genome. As most of these are observed in the heterozygous state, they test a gene’s tolerance to haploinsufficiency and dominant loss of function. We analyzed the distribution of truncating variants across 16,260 autosomal protein coding genes in 11,546 individuals. We observed 39,893 truncating variants affecting 12,062 genes, which significantly differed from an expectation of 12,916 genes under a model of neutral de novo mutation (p<10⁻⁴). Extrapolating this to increasing numbers of sequenced individuals, we estimate that 10.8% of human genes do not tolerate heterozygous truncating variants. An additional 10 to 15% of truncated genes may be rescued by incomplete penetrance or compensatory mutations, or because the truncating variants are of limited functional impact. The study of protein truncating variants delineates the essential genome and, more generally, identifies rare heterozygous variants as an unexplored source of diversity of phenotypic traits and diseases.     

Segregation of a novel p.(Ser270Tyr) <Emphasis Type="Italic">MAF</Emphasis> mutation and p.(Tyr56∗) <Emphasis Type="Italic">CRYGD</Emphasis> variant in a family with dominantly inherited congenital cataracts

A bilaterally blind woman, with a three generation family history of autosomal dominant congenital cataracts, variably associated with iris colobomata and microcornea, sought preconception genetic consultation. Whole-exome sequencing was performed in three affected family members, one unaffected first degree relative, and one spouse. The sequence variant c.168C>G; p.(Tyr56?) in CRYGD, previously reported as pathogenic, and a novel mutation c.809C>A; p.(Ser270Tyr) in MAF, were identified in two affected family members; the grandmother, and half-brother of the proband. The proband inherited only the MAF mutation, whereas her clinically unaffected sister had the CRYGD change. In silico analysis supported a pathogenic role of p.(Ser270Tyr) in MAF, which was absent from publicly available whole-exome datasets, and 1161 Czech individuals. The frequency of CRYGD p.(Tyr56?) in the ExAC dataset was higher than the estimated incidence of congenital cataract in the general population. Our study highlights that patients with genetically heterogeneous conditions may exhibit rare variants in more than one disease-associated gene, warranting caution with data interpretation, and supporting parallel screening of all genes known to harbour pathogenic mutations for a given phenotype. The pathogenicity of sequence variants previously reported as cataract-causing may require re-assessment in light of recently released datasets of human genomic variation.     

Genetic and functional analyses of SHANK2 mutations suggest a multiple hit model of autism spectrum disorders

Autism spectrum disorders (ASD) are a heterogeneous group of neurodevelopmental disorders with a complex inheritance pattern. While many rare variants in synaptic proteins have been identified in patients with ASD, little is known about their effects at the synapse and their interactions with other genetic variations. Here, following the discovery of two de novo SHANK2 deletions by the Autism Genome Project, we identified a novel 421 kb de novo SHANK2 deletion in a patient with autism. We then sequenced SHANK2 in 455 patients with ASD and 431 controls and integrated these results with those reported by Berkel et al. 2010 (n = 396 patients and n = 659 controls). We observed a significant enrichment of variants affecting conserved amino acids in 29 of 851 (3.4%) patients and in 16 of 1,090 (1.5%) controls (P = 0.004, OR = 2.37, 95% CI = 1.23–4.70). In neuronal cell cultures, the variants identified in patients were associated with a reduced synaptic density at dendrites compared to the variants only detected in controls (P = 0.0013). Interestingly, the three patients with de novo SHANK2 deletions also carried inherited CNVs at 15q11–q13 previously associated with neuropsychiatric disorders. In two cases, the nicotinic receptor CHRNA7 was duplicated and in one case the synaptic translation repressor CYFIP1 was deleted. These results strengthen the role of synaptic gene dysfunction in ASD but also highlight the presence of putative modifier genes, which is in keeping with the “multiple hit model” for ASD. A better knowledge of these genetic interactions will be necessary to understand the complex inheritance pattern of ASD.     

Genomic landscapes of Chinese sporadic autism spectrum disorders revealed by whole-genome sequencing

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with considerable clinical and genetic heterogeneity.In this study,we identified all classes of genomic variants from whole-genome sequencing (WGS) dataset of 32 Chinese trios with ASD,including de novo mutations,inherited variants,copy number variants (CNVs) and genomic structural variants.A higher mutation rate (Poisson test,P2.2×10~(-16)) in exonic (1.37×10~(-8)) and 3'-UTR regions (1.42×10~(-8)) was revealed in comparison with that of whole genome (1.05×10~(-8)).Using an integrated model,we identified 87 potentially risk genes (P0.01) from 4832 genes harboring various rare deleterious variants,including CHD8 and NRXN2,implying that the disorders may be in favor to multiple-hit.In particular,frequent rare inherited mutations of several microcephaly-associated genes (ASPM,WDR62,and ZNF335)were found in ASD.In chromosomal structure analyses,we found four de novo CNVs and one de novo chromosomal rearrangement event,including a de novo duplication of UBE3A-containing region at 15q11.2-q13.1,which causes Angelman syndrome and microcephaly,and a disrupted TNR due to de novo chromosomal translocation t (1;5) (q25.1;q33.2).Taken together,our results suggest that abnormalities of centrosomal function and chromatin remodeling of the microcephaly-associated genes may be implicated in pathogenesis of ASD.Adoption of WGS as a new yet efficient technique to illustrate the full genetic spectrum in complex disorders,such as ASD,could provide novel insights into pathogenesis,diagnosis and treatment.     

Investigation of Rare Single-Nucleotide PCDH15 Variants in Schizophrenia and Autism Spectrum Disorders

Both schizophrenia (SCZ) and autism spectrum disorders (ASD) are neuropsychiatric disorders with overlapping genetic etiology. Protocadherin 15 (PCDH15), which encodes a member of the cadherin super family that contributes to neural development and function, has been cited as a risk gene for neuropsychiatric disorders. Recently, rare variants of large effect have been paid attention to understand the etiopathology of these complex disorders. Thus, we evaluated the impacts of rare, single-nucleotide variants (SNVs) in PCDH15 on SCZ or ASD. First, we conducted coding exon-targeted resequencing of PCDH15 with next-generation sequencing technology in 562 Japanese patients (370 SCZ and 192 ASD) and detected 16 heterozygous SNVs. We then performed association analyses on 2,096 cases (1,714 SCZ and 382 ASD) and 1,917 controls with six novel variants of these 16 SNVs. Of these six variants, four (p.R219K, p.T281A, p.D642N, c.3010-1G>C) were ultra-rare variants (minor allele frequency < 0.0005) that may increase disease susceptibility. Finally, no statistically significant association between any of these rare, heterozygous PCDH15 point variants and SCZ or ASD was found. Our results suggest that a larger sample size of resequencing subjects is necessary to detect associations between rare PCDH15 variants and neuropsychiatric disorders.     

Genome-Wide Association Study for Autism Spectrum Disorder in Taiwanese Han Population

Background

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with strong genetic components. Several recent genome-wide association (GWA) studies in Caucasian samples have reported a number of gene regions and loci correlated with the risk of ASD—albeit with very little consensus across studies.

Methods

A two-stage GWA study was employed to identify common genetic variants for ASD in the Taiwanese Han population. The discovery stage included 315 patients with ASD and 1,115 healthy controls, using the Affymetrix SNP array 6.0 platform for genotyping. Several gene regions were then selected for fine-mapping and top markers were examined in extended samples. Single marker, haplotype, gene-based, and pathway analyses were conducted for associations.

Results

Seven SNPs had p-values ranging from 3.4~9.9*10⁻⁶, but none reached the genome-wide significant level. Five of them were mapped to three known genes (OR2M4, STYK1, and MNT) with significant empirical gene-based p-values in OR2M4 (p = 3.4*10⁻⁵) and MNT (p = 0.0008). Results of the fine-mapping study showed single-marker associations in the GLIS1 (rs12082358 and rs12080993) and NAALADL2 (rs3914502 and rs2222447) genes, and gene-based associations for the OR2M3-OR2T5 (olfactory receptor genes, p = 0.02), and GLIPR1/KRR1 gene regions (p = 0.015). Pathway analyses revealed important pathways for ASD, such as olfactory and G protein–coupled receptors signaling pathways.

Conclusions

We reported Taiwanese Han specific susceptibility genes and variants for ASD. However, further replication in other Asian populations is warranted to validate our findings. Investigation in the biological functions of our reported genetic variants might also allow for better understanding on the underlying pathogenesis of autism.     

RET and NRG1 interplay in Hirschsprung disease

Hirschsprung disease (HSCR, aganglionic megacolon) is a complex genetic disorder of the enteric nervous system (ENS) characterized by the absence of enteric neurons along a variable length of the intestine. While rare variants (RVs) in the coding sequence (CDS) of several genes involved in ENS development lead to disease, the association of common variants (CVs) with HSCR has only been reported for RET (the major HSCR gene) and NRG1. Importantly, RVs in the CDS of these two genes are also associated with the disorder. To assess independent and joint effects between the different types of RET and NRG1 variants identified in HSCR patients, we used 254 Chinese sporadic HSCR patients and 143 ethnically matched controls for whom the RET and/or NRG1 variants genotypes (rare and common) were available. Four genetic risk factors were defined and interaction effects were modeled using conditional logistic regression analyses and pair-wise Kendall correlations. Our analysis revealed a joint effect of RET CVs with RET RVs, NRG1 CVs or NRG1 RVs. To assess whether the genetic interaction translated into functional interaction, mouse neural crest cells (NCCs; enteric neuron precursors) isolated from embryonic guts were treated with NRG1 (ErbB2 ligand) or/and GDNF (Ret ligand) and monitored during the subsequent neural differentiation process. Nrg1 inhibited the Gdnf-induced neuronal differentiation and Gdnf negatively regulated Nrg1-signaling by down-regulating the expression of its receptor, ErbB2. This preliminary data suggest that the balance neurogenesis/gliogenesis is critical for ENS development.     

De Novo Occurrence of a Variant in ARL3 and Apparent Autosomal Dominant Transmission of Retinitis Pigmentosa

Background

Retinitis pigmentosa is a phenotype with diverse genetic causes. Due to this genetic heterogeneity, genome-wide identification and analysis of protein-altering DNA variants by exome sequencing is a powerful tool for novel variant and disease gene discovery. In this study, exome sequencing analysis was used to search for potentially causal DNA variants in a two-generation pedigree with apparent dominant retinitis pigmentosa.

Methods

Variant identification and analysis of three affected members (mother and two affected offspring) was performed via exome sequencing. Parental samples of the index case were used to establish inheritance. Follow-up testing of 94 additional retinitis pigmentosa pedigrees was performed via retrospective analysis or Sanger sequencing.

Results and Conclusions

A total of 136 high quality coding variants in 123 genes were identified which are consistent with autosomal dominant disease. Of these, one of the strongest genetic and functional candidates is a c.269A>G (p.Tyr90Cys) variant in ARL3. Follow-up testing established that this variant occurred de novo in the index case. No additional putative causal variants in ARL3 were identified in the follow-up cohort, suggesting that if ARL3 variants can cause adRP it is an extremely rare phenomenon.     

Whole Exome Sequencing in Patients with the Cuticular Drusen Subtype of Age-Related Macular Degeneration

Age-related macular degeneration (AMD) is the leading cause of irreversible blindness in elderly people worldwide. Cuticular drusen (CD) is a clinical subtype of AMD, which typically displays an earlier age at onset, and has a strong genetic component. Genetic studies support a role for rare sequence variants in CD susceptibility, and rare sequence variants in the CFH gene have been identified in 8.8% of CD cases. To further explore the role of rare variants in CD, we performed whole exome sequencing (WES) in 14 affected members of six families and 12 sporadic cases with CD. We detected rare sequence variants in CFH and FBLN5, which previously were shown to harbor rare variants in patients with CD. In addition, we detected heterozygous rare sequence variants in several genes encoding components of the extracellular matrix (ECM), including FBLN1, FBLN3/EFEMP1, FBLN5, FBLN6/HMCN1, FBN2, and COL15A1. Two rare pathogenic variants were identified in the COL15A1 gene: one in a sporadic case and another was found to segregate in a family with six affected individuals with CD. In addition, two rare pathogenic variants were identified in the FGL1 gene in three unrelated CD cases. These findings suggest that alterations in the ECM and in the coagulation pathway may play a role in the pathogenesis of CD. The identified candidate genes require further analyses in larger cohorts to confirm their role in the CD subtype of AMD. No evidence was found of rare sequence variants in a single gene that segregate with CD in the six families, suggesting that the disease is genetically heterogeneous.