Finding genes for complex diseases has been the goal of many genetic studies. Most of these studies have been successful by searching for genes and mutations in rare familial cases, by screening ...candidate genes and by performing genome wide association studies. However, only a small fraction of the total genetic risk for these complex genetic diseases can be explained by the identified mutations and associated genetic loci. In this review we focus on Hirschsprung disease (HSCR) as an example of a complex genetic disorder. We describe the genes identified in this congenital malformation and postulate that both common ‘low penetrant’ variants in combination with rare or private ‘high penetrant’ variants determine the risk on HSCR, and likely, on other complex diseases. We also discuss how new technological advances can be used to gain further insights in the genetic background of complex diseases. Finally, we outline a few steps to develop functional assays in order to determine the involvement of these variants in disease development.
•The genetics of Hirschsprung disease resembles other complex diseases.•The risk on complex diseases is likely determined by common and rare variants.•Statistical selection with systems biology will identify potential disease variants.•Functional analyses of selected variants are needed to infer disease involvement.
Recently, the American College of Medical Genetics (ACMG) recommended the return of actionable secondary findings detected from clinical sequencing. The reported frequency of secondary findings in ...Asian populations were highly variable and it is unclear whether the uniformity in coverage offered by whole-genome sequencing (WGS) may impact the estimate. In this analysis, we aimed to refine the rate of secondary findings on East Asians through a large-scale WGS study. We classified 1256 protein-altering or splicing variants of the 59 actionable genes detected from WGS of 954 East Asians in strict accordance with the ACMG and the Association for Molecular Pathology guidelines. A total of 21 pathogenic or likely pathogenic variants were detected in 24 of the 954 East Asian genomes with an estimate of 2.5% of East Asians carrying actionable variants. Although the overall estimate of secondary findings was consistent with those reported for non-East Asian ethnicities, genetic and allelic heterogeneity was observed. WGS offers a wider breadth of coverage over WES, which highlights the need to further investigate the variable sensitivity of WES and WGS in the detection of secondary findings. Identifying secondary findings in populations underrepresented in previous genetic literature might improve variant interpretation and has a profound impact on local decision-making with regard to the cost-effectiveness of returning the secondary findings from clinical sequencing.
Hirschsprung disease, or congenital aganglionosis, is believed to be oligogenic—that is, caused by multiple genetic factors. We performed whole-genome sequence analyses of patients with Hirschsprung ...disease to identify genetic factors that contribute to disease development and analyzed the functional effects of these variants.
We performed whole-genome sequence analyses of 443 patients with short-segment disease, recruited from hospitals in China and Vietnam, and 493 ethnically matched individuals without Hirschsprung disease (controls). We performed genome-wide association analyses and gene-based rare-variant burden tests to identify rare and common disease–associated variants and study their interactions. We obtained induced pluripotent stem cell (iPSC) lines from 4 patients with Hirschsprung disease and 2 control individuals, and we used these to generate enteric neural crest cells for transcriptomic analyses. We assessed the neuronal lineage differentiation capability of iPSC-derived enteric neural crest cells using an in vitro differentiation assay.
We identified 4 susceptibility loci, including 1 in the phospholipase D1 gene (PLD1) (P = 7.4 × 10–7). The patients had a significant excess of rare protein-altering variants in genes previously associated with Hirschsprung disease and in the β-secretase 2 gene (BACE2) (P = 2.9 × 10–6). The epistatic effects of common and rare variants across these loci provided a sensitized background that increased risk for the disease. In studies of the iPSCs, we observed common and distinct pathways associated with variants in RET that affect risk. In functional assays, we found variants in BACE2 to protect enteric neurons from apoptosis. We propose that alterations in BACE1 signaling via amyloid β precursor protein and BACE2 contribute to pathogenesis of Hirschsprung disease.
In whole-genome sequence analyses of patients with Hirschsprung disease, we identified rare and common variants associated with disease risk. Using iPSC cells, we discovered some functional effects of these variants.
Display omitted
It is widely recognized that noncoding genetic variants play important roles in many human diseases, but there are multiple challenges that hinder the identification of functional disease-associated ...noncoding variants. The number of noncoding variants can be many times that of coding variants; many of them are not functional but in linkage disequilibrium with the functional ones; different variants can have epistatic effects; different variants can affect the same genes or pathways in different individuals; and some variants are related to each other not by affecting the same gene but by affecting the binding of the same upstream regulator. To overcome these difficulties, we propose a novel analysis framework that considers convergent impacts of different genetic variants on protein binding, which provides multiscale information about disease-associated perturbations of regulatory elements, genes, and pathways. Applying it to our whole-genome sequencing data of 918 short-segment Hirschsprung disease patients and matched controls, we identify various novel genes not detected by standard single-variant and region-based tests, functionally centering on neural crest migration and development. Our framework also identifies upstream regulators whose binding is influenced by the noncoding variants. Using human neural crest cells, we confirm cell stage-specific regulatory roles of three top novel regulatory elements on our list, respectively in the
,
, and
loci. In the
regulatory element, we further show that a noncoding variant found only in the patients affects the binding of the gliogenesis regulator NFIA, with a corresponding up-regulation of multiple genes in the same topologically associating domain.
Schizophrenia is a highly heritable, severe psychiatric disorder affecting approximately 1% of the world population. A substantial portion of heritability is still unexplained and the pathophysiology ...of schizophrenia remains to be elucidated. To identify more schizophrenia susceptibility loci, we performed a genome-wide association study (GWAS) on 498 patients with schizophrenia and 2025 controls from the Han Chinese population, and a follow-up study on 1027 cases and 1005 controls. In the follow-up study, we included 384 single nucleotide polymorphisms (SNPs) which were selected from the top hits in our GWAS (130 SNPs) and from previously implicated loci for schizophrenia based on the SZGene database, NHGRI GWAS Catalog, copy number variation studies, GWAS meta-analysis results from the international Psychiatric Genomics Consortium (PGC) and candidate genes from plausible biological pathways (254 SNPs). Within the chromosomal region Xq28, SNP rs2269372 in RENBP achieved genome-wide significance with a combined P value of 3.98 × 10(-8) (OR of allele A = 1.31). SNPs with suggestive P values were identified within 2 genes that have been previously implicated in schizophrenia, MECP2 (rs2734647, P combined = 8.78 × 10(-7), OR = 1.28; rs2239464, P combined = 6.71 × 10(-6), OR = 1.26) and ARHGAP4 (rs2269368, P combined = 4.74 × 10(-7), OR = 1.25). In addition, the patient sample in our follow-up study showed a significantly greater burden for pre-defined risk alleles based on the SNPs selected than the controls. This indicates the existence of schizophrenia susceptibility loci among the SNPs we selected. This also further supports multigenic inheritance in schizophrenia. Our findings identified a new schizophrenia susceptibility locus on Xq28, which harbor the genes RENBP, MECP2, and ARHGAP4.
Hirschsprung disease (HSCR) is a complex genetic disease characterized by absence of ganglia in the intestine. HSCR etiology can be explained by a unique combination of genetic alterations: rare ...coding variants, predisposing haplotypes and Copy Number Variation (CNV). Approximately 18% of patients have additional anatomical malformations or neurological symptoms (HSCR-AAM). Pinpointing the responsible culprits within a CNV is challenging as often many genes are affected. Therefore, we selected candidate genes based on gene enrichment strategies using mouse enteric nervous system transcriptomes and constraint metrics. Next, we used a zebrafish model to investigate whether loss of these genes affects enteric neuron development in vivo.
Celotno besedilo
Dostopno za:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Exome sequencing studies have facilitated the detection of causal genetic variants in yet-unsolved Mendelian diseases. However, the identification of disease causal genes among a list of candidates ...in an exome sequencing study is still not fully settled, and it is often difficult to prioritize candidate genes for follow-up studies. The inheritance mode provides crucial information for understanding Mendelian diseases, but none of the existing gene prioritization tools fully utilize this information.
We examined the characteristics of Mendelian disease genes under different inheritance modes. The results suggest that Mendelian disease genes with autosomal dominant (AD) inheritance mode are more haploinsufficiency and de novo mutation sensitive, whereas those autosomal recessive (AR) genes have significantly more non-synonymous variants and regulatory transcript isoforms. In addition, the X-linked (XL) Mendelian disease genes have fewer non-synonymous and synonymous variants. As a result, we derived a new scoring system for prioritizing candidate genes for Mendelian diseases according to the inheritance mode. Our scoring system assigned to each annotated protein-coding gene (N = 18 859) three pathogenic scores according to the inheritance mode (AD, AR and XL). This inheritance mode-specific framework achieved higher accuracy (area under curve = 0.84) in XL mode.
The inheritance-mode specific pathogenicity prioritization (ISPP) outperformed other well-known methods including Haploinsufficiency, Recessive, Network centrality, Genic Intolerance, Gene Damage Index and Gene Constraint scores. This systematic study suggests that genes manifesting disease inheritance modes tend to have unique characteristics.
ISPP is included in KGGSeq v1.0 (http://grass.cgs.hku.hk/limx/kggseq/), and source code is available from (https://github.com/jacobhsu35/ISPP.git).
mxli@hku.hkSupplementary information: Supplementary data are available at Bioinformatics online.
Objective
Previous genome‐wide association studies (GWAS), which were mainly based on single‐variant analysis, have identified many systemic lupus erythematosus (SLE) susceptibility loci. However, ...the genetic architecture of this complex disease is far from being understood. The aim of this study was to investigate whether using a gene‐based analysis may help to identify novel loci, by considering global evidence of association from a gene or a genomic region rather than focusing on evidence for individual variants.
Methods
Based on the results of a meta‐analysis of 2 GWAS of SLE conducted in 2 Asian cohorts, we performed an in‐depth gene‐based analysis followed by replication in a total of 4,626 patients and 7,466 control subjects of Asian ancestry. Differential allelic expression was measured by pyrosequencing.
Results
More than one‐half of the reported SLE susceptibility loci showed evidence of independent effects, and this finding is important for understanding the mechanisms of association and explaining disease heritability. ANXA6 was detected as a novel SLE susceptibility gene, with several single‐nucleotide polymorphisms (SNPs) contributing independently to the association with disease. The risk allele of rs11960458 correlated significantly with increased expression of ANXA6 in peripheral blood mononuclear cells from heterozygous healthy control subjects. Several other associated SNPs may also regulate ANXA6 expression, according to data obtained from public databases. Higher expression of ANXA6 in patients with SLE was also reported previously.
Conclusion
Our study demonstrated the merit of using gene‐based analysis to identify novel susceptibility loci, especially those with independent effects, and also demonstrated the widespread presence of loci with independent effects in SLE susceptibility genes.
Tetralogy of Fallot (TOF) is the most common cyanotic heart defect, yet the underlying genetic mechanisms remain poorly understood. Here, we performed whole-genome sequencing analysis on 146 ...nonsyndromic TOF parent-offspring trios of Chinese ethnicity. Comparison of de novo variants and recessive genotypes of this data set with data from a European cohort identified both overlapping and potentially novel gene loci and revealed differential functional enrichment between cohorts. To assess the impact of these mutations on early cardiac development, we integrated single-cell and spatial transcriptomics of early human heart development with our genetic findings. We discovered that the candidate gene expression was enriched in the myogenic progenitors of the cardiac outflow tract. Moreover, subsets of the candidate genes were found in specific gene coexpression modules along the cardiomyocyte differentiation trajectory. These integrative functional analyses help dissect the pathogenesis of TOF, revealing cellular hotspots in early heart development resulting in cardiac malformations.
In mammals, urorectal development starts at early embryonic stage, defective urorectal development results in anorectal malformations, which are common congenital developmental defects of the anus ...and the urethra in newborns. The etiology and embryology of the defects are still largely unknown. Platelet-derived growth factor receptor alpha (Pdgfra) is a cell surface receptor tyrosine kinase, upon binding to its ligands (Pdgfa-d), mediates intracellular signaling and regulates embryonic development. The expression of Pdgfra is tightly regulated in the developing urorectal mesenchyme, and its dysregulation is associated with urorectal defects in animals with urorectal defects. Knockout of Pdgfra induces early embryo lethality which precludes investigation of Pdgfra in urorectal development. To address the temporal requirement of Pdgfra in urorectal development, we conditionally deleted Pdgfra in Pdgfra-expressing tissues using a tamoxifen inducible Cre-loxP approach in mice, examined the urorectal development in Pdgfra conditional knockout (Pdgfra-cKO) embryos. We showed that conditional deletion of Pdgfra in Pdgfra-expressing tissues at E10-E11 caused cloaca septation defect, anteriorly displaced anus, defective urogenital folds development and abnormal urethra tubularization in both male and female mice. Furthermore, we showed that Pdgfra was required for the survival of urorectal mesenchyme, deletion of Pdgfra caused apoptosis in the peri-cloacal, the peri-urethra and the urorectal septum mesenchyme of Pdgfra-cKO mutants, associated with an induction of p53, Ndrg1 and activation of caspase-3 in Pdgfra-cKO embryos. In conclusion, Pdgfra is required for the development and survival of the urorectal mesenchyme in embryo, dysregulated Pdgfra signaling induced urorectal defects in mice resembling human congenital diseases of anorectal malformations and hypospadias. Perturbation of PDGFRA signaling may contribute to anorectal malformations and hypospadias in human.