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.     

Family‐based genome scan for age at onset of late‐onset Alzheimer's disease in whole exome sequencing data

Alzheimer's disease (AD) is a common and complex neurodegenerative disease. Age at onset (AAO) of AD is an important component phenotype with a genetic basis, and identification of genes in which variation affects AAO would contribute to identification of factors that affect timing of onset. Increase in AAO through prevention or therapeutic measures would have enormous benefits by delaying AD and its associated morbidities. In this paper, we performed a family‐based genome‐wide association study for AAO of late‐onset AD in whole exome sequence data generated in multigenerational families with multiple AD cases. We conducted single marker and gene‐based burden tests for common and rare variants, respectively. We combined association analyses with variance component linkage analysis, and with reference to prior studies, in order to enhance evidence of the identified genes. For variants and genes implicated by the association study, we performed a gene‐set enrichment analysis to identify potential novel pathways associated with AAO of AD. We found statistically significant association with AAO for three genes (WRN, NTN4 and LAMC3) with common associated variants, and for four genes (SLC8A3, SLC19A3, MADD and LRRK2) with multiple rare‐associated variants that have a plausible biological function related to AD. The genes we have identified are in pathways that are strong candidates for involvement in the development of AD pathology and may lead to a better understanding of AD pathogenesis.     

Multi‐genic pattern found in rare type of hypopituitarism: a whole‐exome sequencing study of Han Chinese with pituitary stalk interruption syndrome

Pituitary stalk interruption syndrome (PSIS) is a rare type of hypopituitarism manifesting various degrees of pituitary hormone deficiency. Although mutations have been identified in some familial cases, the underpinning mechanisms of sporadic patients with PSIS who are in a vast majority remain elusive, necessitating a comprehensive study using systemic approaches. We postulate that other genetic mechanisms may be responsible for the sporadic PSIS. To test this hypothesis, we conducted a study in 24 patients with PSIS of Han Chinese with no family history using whole‐exome sequencing (WES) and bioinformatic analysis. We identified a group of heterozygous mutations in 92% (22 of 24) of the patients, and these genes are mostly associated with Notch, Shh, Wnt signalling pathways. Importantly, 83% (20 of 24) of the patients had more than one mutation in those pathways suggesting synergy of compound mutations underpin the pathogenesis of sporadic PSIS.     

Analytic power and sample size calculation for the genotypic transmission/disequilibrium test in case‐parent trio studies

Case‐parent trio studies considering genotype data from children affected by a disease and their parents are frequently used to detect single nucleotide polymorphisms (SNPs) associated with disease. The most popular statistical tests for this study design are transmission/disequilibrium tests (TDTs). Several types of these tests have been developed, for example, procedures based on alleles or genotypes. Therefore, it is of great interest to examine which of these tests have the highest statistical power to detect SNPs associated with disease. Comparisons of the allelic and the genotypic TDT for individual SNPs have so far been conducted based on simulation studies, since the test statistic of the genotypic TDT was determined numerically. Recently, however, it has been shown that this test statistic can be presented in closed form. In this article, we employ this analytic solution to derive equations for calculating the statistical power and the required sample size for different types of the genotypic TDT. The power of this test is then compared with the one of the corresponding score test assuming the same mode of inheritance as well as the allelic TDT based on a multiplicative mode of inheritance, which is equivalent to the score test assuming an additive mode of inheritance. This is, thus, the first time the power of these tests are compared based on equations, yielding instant results and omitting the need for time‐consuming simulation studies. This comparison reveals that these tests have almost the same power, with the score test being slightly more powerful.     

Targeted exome sequencing of unselected heavy‐ion beam‐irradiated populations reveals less‐biased mutation characteristics in the rice genome

Heavy‐ion beams have been widely utilized as a novel and effective mutagen for mutation breeding in diverse plant species, but the induced mutation spectrum is not fully understood at the genome scale. We describe the development of a multiplexed and cost‐efficient whole‐exome sequencing procedure in rice, and its application to characterize an unselected population of heavy‐ion beam‐induced mutations. The bioinformatics pipeline identified single‐nucleotide mutations as well as small and large (>63 kb) insertions and deletions, and showed good agreement with the results obtained with conventional polymerase chain reaction (PCR) and sequencing analyses. We applied the procedure to analyze the mutation spectrum induced by heavy‐ion beams at the population level. In total, 165 individual M₂ lines derived from six irradiation conditions as well as eight pools from non‐irradiated ‘Nipponbare’ controls were sequenced using the newly established target exome sequencing procedure. The characteristics and distribution of carbon‐ion beam‐induced mutations were analyzed in the absence of bias introduced by visual mutant selections. The average (±SE) number of mutations within the target exon regions was 9.06 ± 0.37 induced by 150 Gy irradiation of dry seeds. The mutation frequency changed in parallel to the irradiation dose when dry seeds were irradiated. The total number of mutations detected by sequencing unselected M₂ lines was correlated with the conventional mutation frequency determined by the occurrence of morphological mutants. Therefore, mutation frequency may be a good indicator for sequencing‐based determination of the optimal irradiation condition for induction of mutations.     

Genome‐wide association studies on the phyllosphere microbiome: Embracing complexity in host–microbe interactions

Environmental sequencing shows that plants harbor complex communities of microbes that vary across environments. However, many approaches for mapping plant genetic variation to microbe‐related traits were developed in the relatively simple context of binary host–microbe interactions under controlled conditions. Recent advances in sequencing and statistics make genome‐wide association studies (GWAS) an increasingly promising approach for identifying the plant genetic variation associated with microbes in a community context. This review discusses early efforts on GWAS of the plant phyllosphere microbiome and the outlook for future studies based on human microbiome GWAS. A workflow for GWAS of the phyllosphere microbiome is then presented, with particular attention to how perspectives on the mechanisms, evolution and environmental dependence of plant–microbe interactions will influence the choice of traits to be mapped.     

Combined analysis of whole‐exon sequencing and lncRNA sequencing in type 2 diabetes mellitus patients with obesity

This study sought to find more exon mutation sites and lncRNA candidates associated with type 2 diabetes mellitus (T2DM) patients with obesity (O‐T2DM). We used O‐T2DM patients and healthy individuals to detect mutations in their peripheral blood by whole‐exon sequencing. And changes in lncRNA expression caused by mutation sites were studied at the RNA level. Then, we performed GO analysis and KEGG pathway analysis. We found a total of 277 377 mutation sites between O‐T2DM and healthy individuals. Then, we performed a DNA‐RNA joint analysis. Based on the screening of harmful sites, 30 mutant genes shared in O‐T2DM patients were screened. At the RNA level, mutations of 106 differentially expressed genes were displayed. Finally, a consensus mutation site and differential expression consensus gene screening were performed. In the current study, the results revealed significant differences in exon sites in peripheral blood between O‐T2DM and healthy individuals, which may play an important role in the pathogenesis of O‐T2DM by affecting the expression of the corresponding lncRNA. This study provides clues to the molecular mechanisms of metabolic disorders in O‐T2DM patients at the DNA and RNA levels, as well as biomarkers of the risk of these disorders.     

Genotyping‐in‐Thousands by sequencing (GT‐seq) panel development and application to minimally invasive DNA samples to support studies in molecular ecology

Minimally invasive sampling (MIS) is widespread in wildlife studies; however, its utility for massively parallel DNA sequencing (MPS) is limited. Poor sample quality and contamination by exogenous DNA can make MIS challenging to use with modern genotyping‐by‐sequencing approaches, which have been traditionally developed for high‐quality DNA sources. Given that MIS is often more appropriate in many contexts, there is a need to make such samples practical for harnessing MPS. Here, we test the ability for Genotyping‐in‐Thousands by sequencing (GT‐seq), a multiplex amplicon sequencing approach, to effectively genotype minimally invasive cloacal DNA samples collected from the Western Rattlesnake (Crotalus oreganus), a threatened species in British Columbia, Canada. As there was no previous genetic information for this species, an optimized panel of 362 SNPs was selected for use with GT‐seq from a de novo restriction site‐associated DNA sequencing (RADseq) assembly. Comparisons of genotypes generated within and among RADseq and GT‐seq for the same individuals found low rates of genotyping error (GT‐seq: 0.50%; RADseq: 0.80%) and discordance (2.57%), the latter likely due to the different genotype calling models employed. GT‐seq mean genotype discordance between blood and cloacal swab samples collected from the same individuals was also minimal (1.37%). Estimates of population diversity parameters were similar across GT‐seq and RADseq data sets, as were inferred patterns of population structure. Overall, GT‐seq can be effectively applied to low‐quality DNA samples, minimizing the inefficiencies presented by exogenous DNA typically found in minimally invasive samples and continuing the expansion of molecular ecology and conservation genetics in the genomics era.     

Disease management at the wildlife‐livestock interface: Using whole‐genome sequencing to study the role of elk in Mycobacterium bovis transmission in Michigan,USA

The role of wildlife in the persistence and spread of livestock diseases is difficult to quantify and control. These difficulties are exacerbated when several wildlife species are potentially involved. Bovine tuberculosis (bTB), caused by Mycobacterium bovis, has experienced an ecological shift in Michigan, with spillover from cattle leading to an endemically infected white‐tailed deer (deer) population. It has potentially substantial implications for the health and well‐being of both wildlife and livestock and incurs a significant economic cost to industry and government. Deer are known to act as a reservoir of infection, with evidence of M. bovis transmission to sympatric elk and cattle populations. However, the role of elk in the circulation of M. bovis is uncertain; they are few in number, but range further than deer, so may enable long distance spread. Combining Whole Genome Sequences (WGS) for M. bovis isolates from exceptionally well‐observed populations of elk, deer and cattle with spatiotemporal locations, we use spatial and Bayesian phylogenetic analyses to show strong spatiotemporal admixture of M. bovis isolates. Clustering of bTB in elk and cattle suggests either intraspecies transmission within the two populations, or exposure to a common source. However, there is no support for significant pathogen transfer amongst elk and cattle, and our data are in accordance with existing evidence that interspecies transmission in Michigan is likely only maintained by deer. This study demonstrates the value of whole genome population studies of M. bovis transmission at the wildlife‐livestock interface, providing insights into bTB management in an endemic system.     

LMO1 polymorphisms and the risk of neuroblastoma: Assessment of meta‐analysis of case‐control studies

Neuroblastoma (NB), a neuroendocrine tumour, is one of the most prevalent cancers in children. The link between LMO1 polymorphisms and NB has been investigated by several groups, rendering inconclusive results. Here, with this comprehensive systematic review and up‐to‐date meta‐analysis, we aim to distinctively elucidate the possible correlation between LMO1 polymorphisms and NB susceptibility. Eligible studies were systematically researched and identified using PubMed, Web of Science and Scopus databases up to 10 February 2019. Odds ratios (ORs) and 95% confidence intervals (CIs) were calculated to assess the strength of the associations. Our findings revealed that rs110419 and rs2168101 polymorphisms were significantly associated with a decreased risk of NB in all genetic models. In addition, the rs4758051 variant appeared protective against NB in homozygous, dominant and allele genetic models, whereas the rs10840002 variant markedly decreased the risk of NB in the allele model. In contrast, the rs204938 polymorphism showed a positive association with NB susceptibility in allele genetic models. In summary, our meta‐analysis is the first to provide clear evidence of an association between specific polymorphisms of LMO1 and susceptibility to NB. Of note, additional larger well‐designed studies would be helpful to further evaluate and confirm this association.     

A survey of single nucleotide polymorphisms identified from whole‐genome sequencing and their functional effect in the porcine genome,

Genetic variants detected from sequence have been used to successfully identify causal variants and map complex traits in several organisms. High and moderate impact variants, those expected to alter or disrupt the protein coded by a gene and those that regulate protein production, likely have a more significant effect on phenotypic variation than do other types of genetic variants. Hence, a comprehensive list of these functional variants would be of considerable interest in swine genomic studies, particularly those targeting fertility and production traits. Whole‐genome sequence was obtained from 72 of the founders of an intensely phenotyped experimental swine herd at the U.S. Meat Animal Research Center (USMARC). These animals included all 24 of the founding boars (12 Duroc and 12 Landrace) and 48 Yorkshire–Landrace composite sows. Sequence reads were mapped to the Sscrofa10.2 genome build, resulting in a mean of 6.1 fold (×) coverage per genome. A total of 22 342 915 high confidence SNPs were identified from the sequenced genomes. These included 21 million previously reported SNPs and 79% of the 62 163 SNPs on the PorcineSNP60 BeadChip assay. Variation was detected in the coding sequence or untranslated regions (UTRs) of 87.8% of the genes in the porcine genome: loss‐of‐function variants were predicted in 504 genes, 10 202 genes contained nonsynonymous variants, 10 773 had variation in UTRs and 13 010 genes contained synonymous variants. Approximately 139 000 SNPs were classified as loss‐of‐function, nonsynonymous or regulatory, which suggests that over 99% of the variation detected in our pigs could potentially be ignored, allowing us to focus on a much smaller number of functional SNPs during future analyses.     

Whole exome sequencing identified a novel DAG1 mutation in a patient with rare,mild and late age of onset muscular dystrophy‐dystroglycanopathy

Muscular dystrophy‐dystroglycanopathy (limb‐girdle), type C, 9 (MDDGC9) is the rarest type of autosomal recessive muscular dystrophies. MDDGC9 is manifested with an early onset in childhood. Patients with MDDGC9 usually identified with defective glycosylation of DAG1, hence it is known as “dystroglycanopathies”. Here, we report a Chinese pedigree presented with mild MDDGC9. The proband is a 64 years old Chinese man. In this family, both the proband and proband's younger brother have been suffering from mild and late onset MDDGC9. Muscle biopsy showed that the left deltoid muscle with an advanced stage of dystrophic change. Immunohistochemistry staining of dystrophin, α‐sarcoglycan, β‐sarcoglycan and dysferlin are normal. Molecular genetic analysis of the proband has been done with whole exome sequencing. A homozygous novel missense mutation (c.2326C>T; p.R776C) in the exon 3 of the DAG1 gene has been identified in the proband. Sanger sequencing revealed that this missense mutation is co‐segregated well among the affected and unaffected (carrier) family members. This mutation is not detected in 200 normal healthy control individuals. This novel homozygous missense mutation (c.2326C>T) causes substitution of arginine by cystine at the position of 776 (p.R776C) which is evolutionarily highly conserved. Immunoblotting studies revealed that a significant reduction of α‐dystroglycan expression in the muscle tissue. The novelty of our study is that it is a first report of DAG1 associated muscular dystrophy‐dystroglycanopathy (limb‐girdle), type C, 9 (MDDGC9) with mild and late age of onset. In Chinese population this is the first report of DAG1 associated MDDGC9.     

Whole‐exome sequencing identifies a novel mutation of GPD1L (R189X) associated with familial conduction disease and sudden death

Cardiac conduction disease (CCD) is a serious disorder and the leading cause of mortality worldwide. It is characterized by arrhythmia, syncope or even sudden cardiac death caused by the dysfunction of cardiac voltage‐gated channel. Previous study has demonstrated that mutations in genes encoding voltage‐gated channel and related proteins were the crucial genetic lesion of CCD. In this study, we employed whole‐exome sequencing to explore the potential causative genes in a Chinese family with ventricular tachycardia and syncope. A novel nonsense mutation (c.565C>T/p.R189X) of glycerol‐3‐phosphate dehydrogenase‐like (GPD1L) was identified and co‐segregated with the affected family members. GPD1L is a crucial interacting protein of SCN5A, a gene encoded sodium channel α‐subunit Na_v1.5 and mainly associated with Brugada syndrome (BrS). The novel mutation (c.565C>T/p.R189X) may result in a premature stop codon at position 189 in exon 4 of the GPD1L gene and lead to functional haploinsufficiency of GPD1L due to mRNA carrying this mutation will be degraded by nonsense‐mediated mRNA decay, which has been confirmed by Western blot in HEK293 cells transfected HIS‐GPD1L plasmid. The levels of GPD1L decreasing may disturb the function of Na_v1.5 and induce arrhythmia and syncope in the end. In conclusion, our study not only further supported the important role of GPD1L in CCD, but also expanded the spectrum of GPD1L mutations and will contribute to the genetic diagnosis and counselling of families with CCD.