Whole-exome sequencing reveals genetic risks of early-onset sporadic dilated cardiomyopathy in the Chinese Han population

To reveal genetic risks of early-onset sporadic dilated cardiomyopathy (DCM) patients in the Chinese Han population, we enlisted 363 DCM cases and 414 healthy controls. Whole-exome sequencing and phenotypic characterization were conducted. In total, we identified 26 loss-of-function (LOF) candidates and 66 pathogenic variants from 33 genes, most of which were novel. The deleterious variants can account for 25.07% (91/363) of all patients. Furthermore, rare missense variants in 21 genes were found to be significantly associated with DCM in burden tests. Other than rare variants, twelve common SNPs were significantly associated with an increased risk of DCM in allele-based genetic model association analysis. Of note, in the cumulative risk model, high-risk subjects had a 3.113-fold higher risk of developing DCM than low-risk subjects. Also, DCM in the high-risk group had a younger age of onset than that in the low-risk group. In terms of cardiac function, the mean left ventricular ejection fraction of patients with the deleterious variants was lower than those without (27.73%±10.02% vs. 30.61%±10.85%, P=0.026). To conclude, we mapped a comprehensive atlas of genetic risks in Chinese patients with DCM that might lead to new insights into the mechanisms and risk stratification for DCM.

     

Actionable,Pathogenic Incidental Findings in 1,000 Participants’ Exomes

The incorporation of genomics into medicine is stimulating interest on the return of incidental findings (IFs) from exome and genome sequencing. However, no large-scale study has yet estimated the number of expected actionable findings per individual; therefore, we classified actionable pathogenic single-nucleotide variants in 500 European- and 500 African-descent participants randomly selected from the National Heart, Lung, and Blood Institute Exome Sequencing Project. The 1,000 individuals were screened for variants in 114 genes selected by an expert panel for their association with medically actionable genetic conditions possibly undiagnosed in adults. Among the 1,000 participants, 585 instances of 239 unique variants were identified as disease causing in the Human Gene Mutation Database (HGMD). The primary literature supporting the variants’ pathogenicity was reviewed. Of the identified IFs, only 16 unique autosomal-dominant variants in 17 individuals were assessed to be pathogenic or likely pathogenic, and one participant had two pathogenic variants for an autosomal-recessive disease. Furthermore, one pathogenic and four likely pathogenic variants not listed as disease causing in HGMD were identified. These data can provide an estimate of the frequency (∼3.4% for European descent and ∼1.2% for African descent) of the high-penetrance actionable pathogenic or likely pathogenic variants in adults. The 23 participants with pathogenic or likely pathogenic variants were disproportionately of European (17) versus African (6) descent. The process of classifying these variants underscores the need for a more comprehensive and diverse centralized resource to provide curated information on pathogenicity for clinical use to minimize health disparities in genomic medicine.     

Targeted Exon Capture and Sequencing in Sporadic Amyotrophic Lateral Sclerosis

Amyotrophic lateral sclerosis (ALS) is a devastating neurodegenerative disease that results in progressive degeneration of motor neurons, ultimately leading to paralysis and death. Approximately 10% of ALS cases are familial, with the remaining 90% of cases being sporadic. Genetic studies in familial cases of ALS have been extremely informative in determining the causative mutations behind ALS, especially as the same mutations identified in familial ALS can also cause sporadic disease. However, the cause of ALS in approximately 30% of familial cases and in the majority of sporadic cases remains unknown. Sporadic ALS cases represent an underutilized resource for genetic information about ALS; therefore, we undertook a targeted sequencing approach of 169 known and candidate ALS disease genes in 242 sporadic ALS cases and 129 matched controls to try to identify novel variants linked to ALS. We found a significant enrichment in novel and rare variants in cases versus controls, indicating that we are likely identifying disease associated mutations. This study highlights the utility of next generation sequencing techniques combined with functional studies and rare variant analysis tools to provide insight into the genetic etiology of a heterogeneous sporadic disease.     

Identification of Medically Actionable Secondary Findings in the 1000 Genomes

The American College of Medical Genetics and Genomics (ACMG) recommends that clinical sequencing laboratories return secondary findings in 56 genes associated with medically actionable conditions. Our goal was to apply a systematic, stringent approach consistent with clinical standards to estimate the prevalence of pathogenic variants associated with such conditions using a diverse sequencing reference sample. Candidate variants in the 56 ACMG genes were selected from Phase 1 of the 1000 Genomes dataset, which contains sequencing information on 1,092 unrelated individuals from across the world. These variants were filtered using the Human Gene Mutation Database (HGMD) Professional version and defined parameters, appraised through literature review, and examined by a clinical laboratory specialist and expert physician. Over 70,000 genetic variants were extracted from the 56 genes, and filtering identified 237 variants annotated as disease causing by HGMD Professional. Literature review and expert evaluation determined that 7 of these variants were pathogenic or likely pathogenic. Furthermore, 5 additional truncating variants not listed as disease causing in HGMD Professional were identified as likely pathogenic. These 12 secondary findings are associated with diseases that could inform medical follow-up, including cancer predisposition syndromes, cardiac conditions, and familial hypercholesterolemia. The majority of the identified medically actionable findings were in individuals from the European (5/379) and Americas (4/181) ancestry groups, with fewer findings in Asian (2/286) and African (1/246) ancestry groups. Our results suggest that medically relevant secondary findings can be identified in approximately 1% (12/1092) of individuals in a diverse reference sample. As clinical sequencing laboratories continue to implement the ACMG recommendations, our results highlight that at least a small number of potentially important secondary findings can be selected for return. Our results also confirm that understudied populations will not reap proportionate benefits of genomic medicine, highlighting the need for continued research efforts on genetic diseases in these populations.     

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.     

Rapid detection of gene mutations responsible for non-syndromic aortic aneurysm and dissection using two different methods: resequencing microarray technology and next-generation sequencing

Aortic aneurysm and/or dissection (AAD) is a life-threatening condition, and several syndromes are known to be related to AAD. In this study, two new technologies, resequencing array technology (ResAT) and next-generation sequencing (NGS), were used to analyze eight genes associated with syndromic AAD in 70 patients with non-syndromic AAD. Eighteen sequence variants were detected using both ResAT and NGS. In addition one of these sequence variants was detected by ResAT only and two additional variants by NGS only. Three of the 18 variants are likely to be pathogenic (in 4.3% of AAD patients and in 8.6% of a subset of patients with thoracic AAD), highlighting the importance of genetic analysis in non-syndromic AAD. ResAT and NGS similarly detected most, but not all, of the variants. Resequencing array technology was a rapid and efficient method for detecting most nucleotide substitutions, but was unable to detect short insertions/deletions, and it is impractical to update custom arrays frequently. Next-generation sequencing was able to detect almost all types of mutation, but requires improved informatics methods.     

Whole-exome sequencing reveals insights into genetic susceptibility to Congenital Zika Syndrome

Congenital Zika Syndrome (CZS) is a critical illness with a wide range of severity caused by Zika virus (ZIKV) infection during pregnancy. Life-threatening neurodevelopmental dysfunctions are among the most common phenotypes observed in affected newborns. Risk factors that contribute to susceptibility and response to ZIKV infection may be related to the virus itself, the environment, and maternal genetic background. Nevertheless, the newborn’s genetic contribution to the critical illness is still not elucidated. Here, we aimed to identify possible genetic variants as well as relevant biological pathways that might be associated with CZS phenotypes. For this purpose, we performed a whole-exome sequencing in 40 children born to women with confirmed exposure to ZIKV during pregnancy. We investigated the occurrence of rare harmful single-nucleotide variants (SNVs) possibly associated with inborn errors in genes ontologically related to CZS phenotypes. Moreover, an exome-wide association analysis was also performed using a case-control design (29 CZS cases and 11 controls), for both common and rare variants. Five out of the 29 CZS patients harbored known pathogenic variants likely to contribute to mild to severe manifestations observed. Approximately, 30% of affected individuals carried at least one pathogenic or likely pathogenic SNV in genes candidates to play a role in CZS. Our common variant association analysis detected a suggestive protective effect of the rs2076469 in DISP3 gene (p-value: 1.39 x 10⁻⁵). The IL12RB2 gene (p-value: 2.18x10^-11) also showed an unusual distribution of nonsynonymous rare SNVs in control samples. Finally, genes harboring harmful variants are involved in processes related to CZS phenotypes such as neurological development and immunity. Therefore, both rare and common variations may be likely to contribute as the underlying genetic cause of CZS susceptibility. The variations and pathways identified in this study may also have implications for the development of therapeutic strategies in the future.     

Exome Sequencing Identified NRG3 as a Novel Susceptible Gene of Hirschsprung’s Disease in a Chinese Population

Hirschsprung’s disease (HSCR) is a complex developmental defect characterized by the absence of enteric ganglia in the gastrointestinal tract. Although the genetic defect of enteric nervous system (ENS) was identified to play a critical role in the progress of HSCR, the systemic genetic dissection of HSCR still needs to be clarified. In this study, we firstly performed exome sequencing of two HSCR patients from a Han Chinese family, including the affected mother and son. After the initial quality filtering (coverage?≥?5X and SNP quality score?≥?40) of the raw data, we identified 13,948 and 13,856 single nucleotide variants (SNVs), respectively. We subsequently compared the SNVs against public databases (dbSNP130, HapMap, and 1000 Genome Project) and obtained a total of 15 novel nonsynonymous SNVs in 15 genes, which were shared between these two patients. Follow-up Sanger sequencing and bioinformatics analysis highlighted variant c.853G>A (p.E285K) in NRG3, a gene involved in the development of ENS. In the validation phase, we sequenced all nine exons of NRG3 in 96 additional sporadic HSCR cases and 110 healthy individuals and identified another nonsynonymous variant c.1329G>A (p.M443I) and two synonymous variants c.828G>A (p.T276T) and c.1365T>A (p.P455P) only in the cases. Our results indicated that NRG3 may be a susceptibility gene for HSCR in a Chinese population.     

Germline loss‐of‐function variants in MBD4 are rare in Finnish patients with uveal melanoma

Uveal melanoma (UM) is a rare intraocular cancer with the highest incidence in northern latitudes. Metastases develop in approximately 50% of patients, whereafter the median survival is 13 months. Generally, the mutation burden of these tumors is low. Germline variants predisposing to UM have been previously described in BRCA1‐associated protein 1 (BAP1). Recently, germline and somatic loss‐of‐function (LOF) variants in the methyl‐CpG‐binding domain 4 (MBD4) gene have been found to cause a hypermutated UM, and MBD4 also has been put forward as a gene predisposing to UM. We sequenced for MBD4 germline variants in 440 Finnish patients with UM and identified seven rare exonic missense variants in 16 (3.6%) patients, of which one likely alters MBD4 function. The frequency of likely pathogenic variants in our cohort is 0.23% (1/432; 95% CI, 0.01–1.28). We identified no LOF variants though their frequency in the Finnish population is 0.052%. Thus, our data do not support the suggestion that MBD4 germline variants predispose to UM. Somatic loss of MBD4 might modify the mutational burden in UM and change its response to immune checkpoint inhibitors.     

Examination of Candidate Exonic Variants for Association to Alzheimer Disease in the Amish

Alzheimer disease (AD) is the most common cause of dementia. As with many complex diseases, the identified variants do not explain the total expected genetic risk that is based on heritability estimates for AD. Isolated founder populations, such as the Amish, are advantageous for genetic studies as they overcome heterogeneity limitations associated with complex population studies. We determined that Amish AD cases harbored a significantly higher burden of the known risk alleles compared to Amish cognitively normal controls, but a significantly lower burden when compared to cases from a dataset of unrelated individuals. Whole-exome sequencing of a selected subset of the overall study population was used as a screening tool to identify variants located in the regions of the genome that are most likely to contribute risk. By then genotyping the top candidate variants from the known AD genes and from linkage regions implicated previous studies in the full dataset, new associations could be confirmed. The most significant result (p = 0.0012) was for rs73938538, a synonymous variant in LAMA1 within the previously identified linkage peak on chromosome 18. However, this association is specific to the Amish and did not generalize when tested in a dataset of unrelated individuals. These results suggest that additional risk variation in the Amish remains to be identified and likely resides outside of the classical protein coding gene regions.     

ORAI1 genetic polymorphisms associated with the susceptibility of atopic dermatitis in Japanese and Taiwanese populations

Atopic dermatitis is a chronic inflammatory skin disease. Multiple genetic and environmental factors are thought to be responsible for susceptibility to AD. In this study, we collected 2,478 DNA samples including 209 AD patients and 729 control subjects from Taiwanese population and 513 AD patients and 1027 control subject from Japanese population for sequencing and genotyping ORAI1. A total of 14 genetic variants including 3 novel single-nucleotide polymorphisms (SNPs) in the ORAI1 gene were identified. Our results indicated that a non-synonymous SNP (rs3741596, Ser218Gly) associated with the susceptibility of AD in the Japanese population but not in the Taiwanese population. However, there is another SNP of ORAI1 (rs3741595) associated with the risk of AD in the Taiwanese population but not in the Japanese population. Taken together, our results indicated that genetic polymorphisms of ORAI1 are very likely to be involved in the susceptibility of AD.     

Variants of OTOF and PJVK genes in Chinese patients with auditory neuropathy spectrum disorder

Background

Mutations in OTOF and PJVK genes cause DFNB9 and DFNB59 types of hearing loss, respectively. The patients carrying pathogenic mutations in either of these genes may show the typical phenotype of auditory neuropathy spectrum disorder (ANSD). The aim of the present study was to identify OTOF and PJVK mutations in sporadic ANSD patients.

Methods and Findings

A total of 76 unrelated Chinese non-syndromic ANSD patients were sequenced on the gene OTOF and PJVK exon by exon. Variants were valued in 105 controls with normal hearing to verify the carrying rate. We identified one pathogenic mutation (c.1194T>A) and three novel, possibly pathogenic, variants (c.3570+2T>C, c.4023+1 G>A, and c.1102G>A) in the OTOF gene, and one novel, possibly pathogenic, variant (c.548G>A) in PJVK. Moreover, we found three novel missense mutations within the exons of OTOF.

Conclusions

As we identified 4 and 1 possible pathogenic variants of the OTOF gene and the PJVK gene, respectively, we believe that screening in these genes are important in sporadic ANSD patients. The pathogenicity of these novel mutations needs further study because of their single heterozygous nature. Knowledge on the mutation spectra of these genes in Chinese would be beneficial in understanding the genetic character of this worldwide disease.     

Whole-genome sequencing of monozygotic twins discordant for schizophrenia indicates multiple genetic risk factors for schizophrenia

Schizophrenia is a common disorder with a high heritability, but its genetic architecture is still elusive.We implemented whole-genome sequencing(WGS) analysis of 8 families with monozygotic(MZ) twin pairs discordant for schizophrenia to assess potential association of de novo mutations(DNMs) or inherited variants with susceptibility to schizophrenia. Eight non-synonymous DNMs(including one splicing site) were identified and shared by twins, which were either located in previously reported schizophrenia risk genes(p.V24689 I mutation in TTN, p.S2506 T mutation in GCN1L1, IVS3+1G T in DOCK1) or had a benign to damaging effect according to in silico prediction analysis. By searching the inherited rare damaging or loss-of-function(LOF) variants and common susceptible alleles from three classes of schizophrenia candidate genes, we were able to distill genetic alterations in several schizophrenia risk genes, including GAD1, PLXNA2, RELN and FEZ1. Four inherited copy number variations(CNVs; including a large deletion at 16p13.11) implicated for schizophrenia were identified in four families, respectively. Most of families carried both missense DNMs and inherited risk variants, which might suggest that DNMs, inherited rare damaging variants and common risk alleles together conferred to schizophrenia susceptibility. Our results support that schizophrenia is caused by a combination of multiple genetic factors, with each DNM/variant showing a relatively small effect size.     

Germline variants in oculocutaneous albinism genes and predisposition to familial cutaneous melanoma

Approximately 1%–2% of cutaneous melanoma (CM) is classified as strongly familial. We sought to investigate unexplained CM predisposition in families negative for the known susceptibility genes using next‐generation sequencing of affected individuals. Segregation of germline variants of interest within families was assessed by Sanger sequencing. Several heterozygous variants in oculocutaneous albinism (OCA) genes: TYR, OCA2, TYRP1 and SLC45A2, were present in our CM cohort. OCA is a group of autosomal recessive genetic disorders, resulting in pigmentation defects of the eyes, hair and skin. Missense variants classified as pathogenic for OCA were present in multiple families and some fully segregated with CM. The functionally compromised TYR p.T373K variant was present in three unrelated families. In OCA2, known pathogenic variants: p.V443I and p.N489D, were present in three families and one family, respectively. We identified a likely pathogenic SLC45A2 frameshift variant that fully segregated with CM in a family of four cases. Another four‐case family harboured cosegregating variants (p.A24T and p.R153C) of uncertain functional significance in TYRP1. We conclude that rare, heterozygous variants in OCA genes confer moderate risk for CM.     

Rare variants in the endocytic pathway are associated with Alzheimer’s disease,its related phenotypes,and functional consequences

Late-onset Alzheimer’s disease (LOAD) is the most common type of dementia causing irreversible brain damage to the elderly and presents a major public health challenge. Clinical research and genome-wide association studies have suggested a potential contribution of the endocytic pathway to AD, with an emphasis on common loci. However, the contribution of rare variants in this pathway to AD has not been thoroughly investigated. In this study, we focused on the effect of rare variants on AD by first applying a rare-variant gene-set burden analysis using genes in the endocytic pathway on over 3,000 individuals with European ancestry from three large whole-genome sequencing (WGS) studies. We identified significant associations of rare-variant burden within the endocytic pathway with AD, which were successfully replicated in independent datasets. We further demonstrated that this endocytic rare-variant enrichment is associated with neurofibrillary tangles (NFTs) and age-related phenotypes, increasing the risk of obtaining severer brain damage, earlier age-at-onset, and earlier age-of-death. Next, by aggregating rare variants within each gene, we sought to identify single endocytic genes associated with AD and NFTs. Careful examination using NFTs revealed one significantly associated gene, ANKRD13D. To identify functional associations, we integrated bulk RNA-Seq data from over 600 brain tissues and found two endocytic expression genes (eGenes), HLA-A and SLC26A7, that displayed significant influences on their gene expressions. Differential expressions between AD patients and controls of these three identified genes were further examined by incorporating scRNA-Seq data from 48 post-mortem brain samples and demonstrated distinct expression patterns across cell types. Taken together, our results demonstrated strong rare-variant effect in the endocytic pathway on AD risk and progression and functional effect of gene expression alteration in both bulk and single-cell resolution, which may bring more insight and serve as valuable resources for future AD genetic studies, clinical research, and therapeutic targeting.     

Exome-Based Mapping and Variant Prioritization for Inherited Mendelian Disorders

Exome sequencing in families affected by rare genetic disorders has the potential to rapidly identify new disease genes (genes in which mutations cause disease), but the identification of a single causal mutation among thousands of variants remains a significant challenge. We developed a scoring algorithm to prioritize potential causal variants within a family according to segregation with the phenotype, population frequency, predicted effect, and gene expression in the tissue(s) of interest. To narrow the search space in families with multiple affected individuals, we also developed two complementary approaches to exome-based mapping of autosomal-dominant disorders. One approach identifies segments of maximum identity by descent among affected individuals; the other nominates regions on the basis of shared rare variants and the absence of homozygous differences between affected individuals. We showcase our methods by using exome sequence data from families affected by autosomal-dominant retinitis pigmentosa (adRP), a rare disorder characterized by night blindness and progressive vision loss. We performed exome capture and sequencing on 91 samples representing 24 families affected by probable adRP but lacking common disease-causing mutations. Eight of 24 families (33%) were revealed to harbor high-scoring, most likely pathogenic (by clinical assessment) mutations affecting known RP genes. Analysis of the remaining 17 families identified candidate variants in a number of interesting genes, some of which have withstood further segregation testing in extended pedigrees. To empower the search for Mendelian-disease genes in family-based sequencing studies, we implemented them in a cross-platform-compatible software package, MendelScan, which is freely available to the research community.     

Mutation screening of multiple genes in Spanish patients with autosomal recessive retinitis pigmentosa by targeted resequencing

Retinitis Pigmentosa (RP) is a heterogeneous group of inherited retinal dystrophies characterised ultimately by the loss of photoreceptor cells. RP is the leading cause of visual loss in individuals younger than 60 years, with a prevalence of about 1 in 4000. The molecular genetic diagnosis of autosomal recessive RP (arRP) is challenging due to the large genetic and clinical heterogeneity. Traditional methods for sequencing arRP genes are often laborious and not easily available and a screening technique that enables the rapid detection of the genetic cause would be very helpful in the clinical practice. The goal of this study was to develop and apply microarray-based resequencing technology capable of detecting both known and novel mutations on a single high-throughput platform. Hence, the coding regions and exon/intron boundaries of 16 arRP genes were resequenced using microarrays in 102 Spanish patients with clinical diagnosis of arRP. All the detected variations were confirmed by direct sequencing and potential pathogenicity was assessed by functional predictions and frequency in controls. For validation purposes 4 positive controls for variants consisting of previously identified changes were hybridized on the array. As a result of the screening, we detected 44 variants, of which 15 are very likely pathogenic detected in 14 arRP families (14%). Finally, the design of this array can easily be transformed in an equivalent diagnostic system based on targeted enrichment followed by next generation sequencing.     

Exome Sequencing of Germline DNA from Non-BRCA1/2 Familial Breast Cancer Cases Selected on the Basis of aCGH Tumor Profiling

The bulk of familial breast cancer risk (∼70%) cannot be explained by mutations in the known predisposition genes, primarily BRCA1 and BRCA2. Underlying genetic heterogeneity in these cases is the probable explanation for the failure of all attempts to identify further high-risk alleles. While exome sequencing of non-BRCA1/2 breast cancer cases is a promising strategy to detect new high-risk genes, rational approaches to the rigorous pre-selection of cases are needed to reduce heterogeneity. We selected six families in which the tumours of multiple cases showed a specific genomic profile on array comparative genomic hybridization (aCGH). Linkage analysis in these families revealed a region on chromosome 4 with a LOD score of 2.49 under homogeneity. We then analysed the germline DNA of two patients from each family using exome sequencing. Initially focusing on the linkage region, no potentially pathogenic variants could be identified in more than one family. Variants outside the linkage region were then analysed, and we detected multiple possibly pathogenic variants in genes that encode DNA integrity maintenance proteins. However, further analysis led to the rejection of all variants due to poor co-segregation or a relatively high allele frequency in a control population. We concluded that using CGH results to focus on a sub-set of families for sequencing analysis did not enable us to identify a common genetic change responsible for the aggregation of breast cancer in these families. Our data also support the emerging view that non-BRCA1/2 hereditary breast cancer families have a very heterogeneous genetic basis.