Next Article in Journal
Gut Microbial Composition and Predicted Functions Are Not Associated with Feather Pecking and Antagonistic Behavior in Laying Hens
Next Article in Special Issue
Cross-Validation of Next-Generation Sequencing Technologies for Diagnosis of Chromosomal Mosaicism and Segmental Aneuploidies in Preimplantation Embryos Model
Previous Article in Journal
Classification of the Biogenicity of Complex Organic Mixtures for the Detection of Extraterrestrial Life
Previous Article in Special Issue
Whole-Genome Sequencing in Diagnostics of Selected Slovenian Undiagnosed Patients with Rare Disorders
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Life
Volume 11
Issue 3
10.3390/life11030233
life-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cross-Disorder Analysis of De Novo Variants Increases the Power of Prioritising Candidate Genes

by
Kuokuo Li
1,2,3,†,
Zhengbao Ling
3,†,
Tengfei Luo
3,
Guihu Zhao
1,
Qiao Zhou
1,
Xiaomeng Wang
3,
Kun Xia
3,
Jinchen Li
1,3,4,* and
Bin Li
1,5,*
1
National Clinical Research Center for Geriatric Disorders, Department of Geriatrics, Xiangya Hospital, Central South University, Changsha 410008, China
2
Department of Obstetrics and Gynecology, The First Affiliated Hospital of Anhui Medical University, No 218 Jixi Road, Hefei 230022, China
3
Center for Medical Genetics & Hunan Key Laboratory of Medical Genetics, School of Life Sciences, Central South University, Changsha 410008, China
4
Department of Neurology, Xiangya Hospital, Central South University, Changsha 410008, China
5
Mobile Health Ministry of Education—China Mobile Joint Laboratory, Xiangya Hospital, Central South University, Changsha 410008, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to the paper as first authors.
Life 2021, 11(3), 233; https://doi.org/10.3390/life11030233
Submission received: 29 January 2021 / Revised: 1 March 2021 / Accepted: 6 March 2021 / Published: 12 March 2021
(This article belongs to the Special Issue The Increasing Role of Next Generation Sequencing Methods in Mutation Analysis)
Browse Figures
Versions Notes

Abstract

:
De novo variants (DNVs) are critical to the treatment of neurodevelopmental disorders (NDDs). However, effectively identifying candidate genes in small cohorts is challenging in most NDDs because of high genetic heterogeneity. We hypothesised that integrating DNVs from multiple NDDs with genetic similarity can significantly increase the possibility of prioritising the candidate gene. We catalogued 66,186 coding DNVs in 50,028 individuals with nine types of NDDs in cohorts with sizes spanning from 118 to 31,260 from Gene4Denovo database to validate this hypothesis. Interestingly, we found that integrated DNVs can effectively increase the number of prioritised candidate genes for each disorder. We identified 654 candidate genes including 481 shared candidate genes carrying putative functional variants in at least two disorders. Notably, 13.51% (65/481) of shared candidate genes were prioritised only via integrated analysis including 44.62% (29/65) genes validated in recent large cohort studies. Moreover, we estimated that more novel candidate genes will be prioritised with the increase in cohort size, in particular for some disorders with high putative functional DNVs per individual. In conclusion, integrated DNVs may increase the power of prioritising candidate genes, which is important for NDDs with small cohort size.

1. Introduction

Neurodevelopmental disorders (NDDs) are disorders with high clinical heterogeneity, leading to considerable personal suffering, morbidity, and disability, which also increase the burden of global healthcare [1]. NDDs, including autism spectrum disorder (ASD), congenital heart disease (CHD), developmental disorders/intellectual disability (DD/ID), epileptic encephalopathy (EE), and schizophrenia (SCZ), are diagnosed following clinical practise guidelines created by practiced experts based on symptoms and signs. However, the diagnosis does not follow a uniform biological framework, which leads to the clinical heterogeneity and overlap between different kinds of NDDs. According to previous studies, almost 26% of patients with ASD, DD, or ID, presented with the clinical phenotypes of EE [2]. Additionally, patients with CHD were reported to share significant clinical features [3,4] and genetic components with those with other NDDs [5]; 10% of CHD cases and 50% of severe CHD cases have a similar clinical phenotype to those with other NDDs [6].
With the development of next-generation sequencing technologies, genetic disruptions have been identified as the major reasons for NDDs, especially for de novo variants (DNVs) with important functional effects contributing to early NDDs. Candidate genes in ASD [7], CHD [8], DD [9], EE [10], ID [11], and SCZ [12] were detected and prioritised successfully using whole-exome sequencing (WES) or whole-genome sequencing (WGS) in recent studies. For instance, due to strong functional effects of DNVs, CHD8 [13] was found to be associated with ASD, ID, sleeping problems, macrocephaly, and gastrointestinal symptoms, while DYRK1A [14] is associated with ID, microcephaly, and febrile seizures infancy. However, for most NDDs, only several candidate genes have been identified based on DNVs due to genetic heterogeneity, rarity DNVs, limited cohort size, and small gene-level relative risks. This genotype-phenotype association method was validated in our studies based on GeneMatcher [15,16,17]. To prioritise additional candidate genes, detect DNVs in specific genes, and perform genotype–phenotype analysis, as well as statistical burden analysis, large cohort size is required. However, sequencing new patients with NDDs requires time and funds, which is not the optimal choice. Another method is to integrate DNVs from multiple NDDs with phenotypic similarities to increase the statistical power of candidate gene discovery, which has been validated recently [18,19,20].
In this study, we aimed to validate the performance of cross-disorder analysis in prioritising candidate genes. We catalogued DNVs in 50,028 individuals with nine types of NDDs that had clinical and genetic similarities. Our results demonstrated that integrating DNVs of different disorders effectively increased the number of prioritised candidate genes. All novel candidate genes shared putative functional variants with more than one disorder. Moreover, we found that both gene-level relative risks and cohort size were the major contributors to candidate gene prioritisation.

2. Materials and Methods

2.1. Data Collection and Annotation

We collected DNVs detected using WGS or WES in 48 published studies (Table S1). The redundant DNVs were removed based on the description of these studies. We focused on DNVs in nine kinds of NDDs with cohort sizes spanning from 118 to 31,260, including ASD, SCZ, EE, DD/ID, CHD, Tourette disorder (TD), bipolar disorder (BP), obsessive-compulsive disorder (OCD), and complex motor stereotypies (CMS). ANNOVAR [21] and VarCards [22] were used to annotate DNVs in the human reference genome (hg19). We catalogued DNVs into five classes as follows: (1) Loss-of-function (LoF) variant including splicing (≤2 bp), stopgain, and stoploss SNVs, and frameshift indels; (2) deleterious missense (Dmis) variant; (3) tolerant missense (Tmis) variant; (4) synonymous (Syn) variant; (5) non-frameshift indels (NF) variant. We used ReVe [23] to predict Dmis/Tmis variants. Both LoF and Dmis variants were defined as putative functional (Pfun) variants. All these variants are available in our Gene4Denovo database [24].

2.2. Overlap of Genes across NDDs Based on De Novo Variants

To test genetic similarities among different NDDs, we used DNENRICH [12,25] software, taking gene size, structure, and local trinucleotide variant rate into consideration to test whether one specific kind of DNVs was significantly shared between any two NDDs. We focused on LoF, Dmis, and Pfun, which increase the genetic risk of NDDs. For this test, we counted the variant number of each gene in each disorder and then calculated the number of overlapping variants between any two disorder. Based on the observed number of variants for each disorder, we randomly selected the matching gene number from all human genes and calculated the expected genetic overlap between any two disorders. Permutation tests were performed to estimate genetic similarity. Details of this method can be found at https://fromem03.u.hpc.mssm.edu/dnenrich/ (accessed on 25 September 2018).

2.3. Candidate Genes Prioritization Based on TADA

We next performed transmitted and de novo association (TADA) analysis [26,27] to prioritize candidate genes. In this study, we used TADA-Denovo which identified candidate genes only based on DNMs. TADA-Denovo is a Bayesian model which used observed Pfun DNVs including LoF and deleterious missense variant and expected de novo mutation rate to prioritised candidate gene of disease. Genes carrying significant more Pfun DNVs than expectation were defined as candidate genes. In the first strategy, we counted LoF and Dmis DNVs number in each disorder and performed TADA to calculate the false discovery rate (FDR) for each gene. In the second strategy, LoF and Dmis DNVs in each gene of all NDDs were counted to perform TADA analysis based on the shared genetic components of NDDs. Genes with FDR < 0.05 in these two strategies were defined as candidate genes. Genes carrying Pfun DNVs in more than one disorder were defined as shared genes, and those carrying DNVs in only one disorder were defined as unique genes.

2.4. Predicted Gene Discovery Rate

To determine gene discovery rate of each disorder in increased sample size, we sampled (with replacement) populations of 500, 1000, 2000, 4000, 8000, 16,000, and 32,000 cases, and performed TADA analysis to prioritise candidate genes. As TADA only considers Pfun (LoF and Dmis), we used Pfun per individual of each disorder and the sampling number to calculate Pfun DNVs and then performed TADA analysis. The number of genes with FDR < 0.05 of each prediction was counted.

2.5. Statistical Analysis

We performed statistical analyses by using R software (v3.5.0) and Linux system (vCentOS 7.1). The R code are available in supplementary file and related input files are available from the corresponding author on reasonable request. The genetic similarity between any two NDDs were performed by DNENRICH software (v1.0). TADA software was used to prioritize candidate gene. p-value or FDR in method less than 0.05 were defined as statistically significant. The detail information of two statistical methods were as follows:
DNENRICH (https://fromem03.u.hpc.mssm.edu/dnenrich/) (accessed on 25 September 2018).
TADA (http://www.compgen.pitt.edu/TADA/TADA_guide.html) (accessed on 1 September 2020).

3. Results

We curated a total of 348,812 DNVs from 50,028 patients with nine kinds of NDDs with varying sample sizes as reported in 48 published studies (Table 1; Table S1). The DD/ID (n = 31,260) accounted for the highest proportion of NDDs in this study and included the highest number of trios based on a WES study to detect coding DNVs (n = 44,825) in NDDs. ASD (n = 10,318), SCZ (n = 3402), CHD (n = 2645), EE (n = 973), and TD (n = 909) were also major NDDs in the identification of candidate genes based on DNVs. However, there was only a small number of patients with BP (n = 219), OCD (n = 118), and CMS (n = 184). In addition, we also found significant Pfun enrichment in NDDs, which were not associated with the sample size (Table 1). DD/ID, OCD, and EE exhibited the most DNVs enrichment, revealing 0.60, 0.58, and 0.55 Pfun DNVs per patient, respectively. Whereas ASD, CHD, SCZ, and TD showed 0.40, 0.39, 0.32, and 0.31 Pfun DNVs per patient, respectively, and exhibited lower DNV enrichment followed by BP (0.25) and CMS (0.21) (Table 1).
As LoF, Dmis and Pfun contributed to the formation of NDDs, we explored genetic similarities between any two NDDs using these two kinds of DNVs based on DNENRICH. As nine kinds of NDDs were involved in the genetic similarity analysis, we performed 36 (9 × 8/2) tests. We found that most NDDs showed similar genetic components with others in Dmis, LoF, and Pfun, in particular for ASD and DD/ID. ASD exhibited significant genetic similarity (p < 0.05) with 7/8, 7/8, and 8/8 of other NDDs in Dmis, LoF, and Pfun, respectively (Figure 1). DD/ID exhibited significant genetic similarity (p < 0.05) with 6/8 of other NDDs in Dmis, LoF, and Pfun, respectively (Figure 1). We did not find statistically significant genetic overlapping among NDDs in 44.44% (16/36), 55.56% (20/36), and 41.67% (15/36) of Dmis-, LoF-, and Pfun-based comparisons, respectively, which might be due to genetic heterogeneity, limited cohort size, or small gene-level relative risks, but 94.4% (34/36) of observed genetic overlapping was high than expected (OE > 1) (Figure 1).
Based on the TADA analysis of Pfun DNVs in nine NDDs, we prioritised 589 candidate genes (FDR < 0.05), containing 104, 8, 24, 527, 22, and 1 genes in ASD, SCZ, EE, DD/ID, CHD, and CMS, respectively (Table 2; Table S2). Due to the smaller number of sample size and low contribution of DNVs, we did not prioritise any candidate gene in TD, BP, and OCD based on DNVs in single disorder. Since most NDDs exhibited more genetic similarities than expected (OE > 1), we integrated Pfun DNVs to all NDDs and performed TADA analysis. We prioritised a total of 523 candidate genes with FDR < 0.05, including 65 novel genes that were not included in the above 589 genes (Table S2). Genes carrying putative functional DNVs in a specific disorder, which passed each FDR threshold in the integration analysis was defined as a candidate gene of this disorder. We found that integrated DNVs increased the number of prioritised candidate genes for each disorder from 5 to 258, in particular for disorders that exhibited more genetic overlapping with integrated DNVs or those with large sample size. This was observed in ASD (n = 10,318, p-value = 1.00 × 10−4, OE = 3.72), SCZ (n = 3402, p-value = 1.00 × 10−4, OE = 1.67), EE (n = 973, p-value = 1.00 × 10−4, OE = 6.54), and CHD (n = 2645, p-value = 1.00 × 10−4, OE = 2.84). Moreover, we prioritised putative candidate genes of NDDs with a small cohort size, including TD, OCD, CMS, and BP. For TD and OCD, we prioritised 28 and 14 novel candidate genes, respectively, compared with zero in a single disorder-based analysis. These two disorders exhibited relatively more genetic similarities with other NDDs and allowed us to prioritise more candidate genes. For CMS and BP, we prioritised nine and five novel candidate genes, respectively, compared with one and zero in a single disorder-based analysis.
We prioritised a total of 654 candidate genes with FDR < 0.05 (Table 3; Table S2). Based on the strength of the statistical evidence, we ranked candidate genes into four ranks as follows: Rank 1 (FDR ≤ 0.0001, n = 316); rank 2 (0.0001 < FDR < 0.001, n = 60); rank 3 (0.001 < FDR < 0.01, n = 101); and rank 4 (0.01 < FDR < 0.05, n = 177). Moreover, based on the number of disorders carrying Pfun DNVs of a specific gene, we identified six groups of candidate genes. Precisely 26.45% (173/654) of candidate genes showed Pfun DNMs in only one disorder and 36.54% (239/654), 25.54% (167/654), 9.02% (59/654), 1.99% (13/654), and 0.46% (3/654) of candidate genes showed Pfun DNMs in two, three, four, five, and six disorders, respectively (Table 3). For example, CACNA1E (FDR = 1.00 × 10−11), KMT2C (FDR = 6.31 × 10−15), and KDM5B (FDR < 2.00 × 10−18) showed Pfun in six NDDs. Integrated analysis prioritised 65 novel candidate genes compared to a single disorder-based analysis and all showed Pfun DNVs in more than one disorder (Table S2). To validate these novel candidate genes, we compared them with previously identified candidate genes of NDDs [20,28,29,30,31,32]. In addition, the gene with an expression value > 1 read per kilobase per million map reads in > 50% human brain samples or in >50% human foetal brain samples were defined as a gene expressed in the brain using the BrainSpan database. We found that 44.62% (29/65) of novel candidate genes were reported previously [20,28,29,31,32,33]. Moreover, we used the probability of loss-of-function intolerance (pLI) > 0.9 to filter candidate genes and 55.18% (16/29) passed this threshold, including SPRY2, PSMD12, RALA, CIC, ATP1A1, ZMYND8, BHLHE40, NR6A1, RYR2, GGNBP2, EIF4A2, RAB11A, CTR9, RAB2A, UPF3B, and KCNC1. We did not identify additional 55.38% (36/65) novel genes in previous studies [20,28,29,31,32,33]; however, we identified 30.56% (11/36) of the genes with a pLI > 0.9, including RNF220, AP1G1, TRIM8, LHX2, CRIM1, UBR3, RPSA, WDR20, SUFU, PHEX, and KPNA1.
To strengthen the evidence for candidate genes, we sourced them on a priority basis from genome wide association studies (GWAS) of NDDs based on a threshold of p < 10–5 in GWAS Catalog database (https://www.ebi.ac.uk/gwas/ (accessed on 6 March 2021)). We noted that 14.68% (96/654) of candidate genes that were associated with NDDs by common SNPs (Table S2) which was significant than random expectation (Fisher’s exact test, p = 0.023, OR = 1.29, 95% CI 1.03–1.61). This result provided bidirectional genetic evidence for these genes.
Based on the existing DNVs, we projected the gene discovery rate in an increased sample size. We sampled cohorts with 500, 1000, 2000, 4000, 8000, 16,000, and 32,000 sample sizes for each disorder and performed TADA analysis to prioritise candidate genes. We found that the gene discovery rate rapidly increased with the increasing sample size and then gradually reached a plateau for most NDDs (Figure 2). In addition, we found that putative functional DNVs enrichment was positively correlated with candidate gene discovery rate and negatively correlated with candidate gene discovery plateau. DD/ID, OCD, and EE revealed 0.60, 0.58, and 0.55 putative functional DNVs per patient, respectively, and prioritised more candidate genes with the same cohort size than other disorders. We did not find a candidate gene discovery plateau in the cohort with 32,000 samples. ASD, CHD, SCZ, TD, BP, and CMS exhibited low putative functional DNV enrichment, revealing 0.40, 0.39, 0.32, 0.31, 0.25, and 0.21 Pfun DNVs per patient, respectively. These six NDDs gradually reached a plateau with reduced putative functional DNVs enrichment. For example, BP and CMS plateau in candidate genes were prioritised in the cohort with 16,000 samples.

4. Discussion

In this study, we tried to explore a new method to improve the effectiveness of NDD candidate genes by integrating NDDs with similar clinical features. We observed that DD/ID, OCD, and EE exhibited the most DNVs enrichment, revealing 0.60, 0.58, and 0.55 Pfun DNVs per patient, while BP exhibited 0.25 Pfun DNVs per patient, and CMS exhibited 0.21 Pfun DNVs per patient.
NDDs cause suffering, morbidity, and disability, with challenges in the diagnosis and treatment, due to the high clinical heterogeneity within individual disorders [34]. Next-generation sequencing technologies have revealed that DNVs play an important role, functionally contributing to the development of NDDs. However, the genetic heterogeneity, rarity of DNVs, limited cohort size, and small gene-level relative risks, are major challenges for identifying novel candidate genes, which might result in reduced effectiveness in screening for disease-candidate genes [7,8,9,10,11,12]. Due to the limited cohort size, for CMS (n = 184), only one candidate gene was prioritised and for BP (n = 219), OCD (n = 118), and TD (n = 909) no candidate genes were prioritised in the single disease-based analysis. It is laborious and costly to collect sufficient cohort information and genomic data, especially for the NDDs with low incidence rates. Therefore, it is a great challenge to study NDDs in-depth with small sample size.
Further analyses of genetic components in the present study showed that NDDs presented similar genetic components in Dmis, LoF, and Pfun. We prioritised a total of 523 candidate genes, including 65 novel genes, which could not be screened by analysing a single NDD. This suggested that integrated analyses could reveal potentially useful data that single disorder-based analysis might have missed. By comparing our findings with previously reported NDD candidate genes and analysing the expression of novel genes using the BrainSpan database, we observed that 29 novel genes were expressed in the brain, overlapping with known candidate genes. Exactly 16 out of 29 genes, including SPRY2, PSMD12, RALA, CIC, ATP1A1, ZMYND8, BHLHE40, NR6A1, RYR2, GGNBP2, EIF4A2, RAB11A, CTR9, RAB2A, UPF3B, and KCNC1 were more intolerant to LoF with a pLI > 0.9. This suggested that the results of integrated analyses were highly reliable. In other 36 genes, we detected that 11 genes with LoF, such as RNF220, AP1G1, TRIM8, LHX2, CRIM1, UBR3, RPSA, WDR20, SUFU, PHEX, and KPNA1, were more likely to cause disease.
Compared with previous studies [20,28,29,31,32,33], the 11 genes have not been screened out as the candidate genes of NDDs. Through literature research, we found that most of the 11 novel genes were associated with NDDs. For instance, RNF220 contributes to noradrenergic neuron development [35] and specifies spinal progenitor domains [36]. AP1G1 plays an important role in the PI3K/AKT pathway, which is not only associated with cancers but also with megalencephaly [37], ASD [37,38], neurodevelopmental delay [39], and other NDDs [40]. Both TRIM8 [41] and UBR3 can act on CTNNB1, which can promote the Wnt signalling pathway. Notably, the Wnt signalling pathway is a classical neuron development signalling pathway [42]. Meanwhile, TRIM8 [43] and RPSA [44] are also an important regulators of the PI3K-AKT-mTOR signalling pathway, which is a developmental disease-related signalling pathway [40]. LHX2 is a regulator of neural differentiation [45,46]. Furthermore, LHX2 controls cortical size by regulating the balance between proliferation and differentiation in cortical progenitors [45]. Other novel genes identified in this study are associated with cancer, such as WDR20 with medulloblastoma [47] and SUFU with renal cell carcinoma [48,49]. We previously found that NDDs share a common genomic basis with cancers [42]. Therefore, we hypothesised that these genes contribute effectively to the development of NDDs.
Based on projection estimates, we found that gene discovery was affected by sample size and DNV per individual for each disorder. These results were consistent with previous studies. In our future studies, we will prioritise more novel candidate genes and strengthen the genetic evidence of previous putative candidate genes with the increased sample size. For example, Deciphering Developmental Disorders Study identified 94 candidate genes in 4293 families [9] and 285 genes involved in developmental disorders in 31,058 parent-offspring trios [29]. DNVs are also significant influences for candidate gene discovery. For example, based on the sampling method and hypothetical 3000 trios, about 200 high confidence candidate gene (FDR < 0.1) and 300 probable candidate gene (FDR < 0.3) were proposed for the OCD [50] but for the TD only about 25 high confidence candidate genes (FDR < 0.1) and 80 probable candidate genes (FDR < 0.3) were discovered [51]. Moreover, the degree of functional disruption of the variant can also influence the gene discovery rate. Coe et al. found that gene discovery based on LoF and Dmis variants reaches a plateau and may identify a few novel candidate genes. However, increased cohort size will probably identify more novel candidate genes based on fewer severe de novo missense variants, in particular, missense variant cluster in specific hotspot regions, which was not studied previously [20].
DNVs contribute significant to NDDs and one single DNV might result in the formation of NDDs. Previous studies also found than DNV exhibited potential cumulative effect to NDDs which was defined as “oligogenic model” [52,53,54]. Du et al. found that the number of patients with ASD carrying multiple extreme DNVs are significant more than controls. In addition, they also found that patients with ASD carrying more than 2 DNVs exhibited lower IQ than patients carrying 1 or 0 DNV [55]. This was consistent with another study which found that patients with DD/ID carrying more DNVs than ASD [20]. Gifford et al. found three missense variants contribute to heart disease [56]. These results indicated that genes with deleterious variants in one patient might participate in common biological pathway or one variant work as genetic modifier to other genes.
There are still some limitations to this study. First, the large difference in sample sizes of different NDDs may lead to genetic statistical bias, although we statistically corrected these. Second, only classical LoF including classical splicing, stopgain, and stoploss SNVs, and frameshift indels were involved into conventional TADA analysis. Other kinds of variants such as de novo cryptic splice variants predicted by SpliceAI [57], inframe indels variants and small de novo CNV deletions (SmallDel) could also participate in the formation of NDDs. For example, Ruzzo take SmallDel into adjusted TADA analysis [31]. Third, we used in silico tool to predicted missense variant and not all predicted deleterious missense are pathogenic. The combined of multiple tools are useful to identify truly positive candidate genes [58]. Fourth, the novel candidate genes screened in this study by using bioinformatic tools and further functional experimental verification were necessary to validate whether genes were involved in NDDs. Fifth, environmental factors are also involved in the aetiology of NDDs. Combining environmental factors with genetic factors in further research will improve the accuracy and efficiency of DNVs research.

5. Conclusions

Taken together, it is inferred from the existing evidence that 11 new candidate genes are relatively reliable for further research. Integrated analysis can effectively improve the candidate gene discovery rate in NDDs. This study provides a new idea for genetic research of NDDs with insufficient samples.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-1729/11/3/233/s1, Table S1: WES or WGS studies regarding to DNVs, Table S2: Candidate gene with FDR < 0.05 in this study.

Author Contributions

Conceptualization, K.L., Z.L., T.L., G.Z., Q.Z., X.W., K.X., J.L., and B.L.; data curation, K.L., Z.L., T.L., G.Z., Q.Z., X.W., K.X., J.L., and B.L.; investigation, Z.L., T.L., G.Z., X.W., and J.L.; methodology, K.L., Z.L., T.L., J.L., and B.L.; software, G.Z. and Q.Z.; supervision, K.X., J.L., and B.L.; validation, Z.L.; writing—original draft, K.L., Z.L., J.L., and B.L.; writing—review and editing, Z.L., T.L., G.Z., Q.Z., X.W., K.X., J.L., and B.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China, grant number 82001362; Hunan Natural Science Foundation Outstanding Youth Fund, grant number 2020JJ3059; Young Elite Scientist Sponsorship Program by CAST, grant number 2018QNRC001; and Innovation-Driven Project of Central South University, grant number 20180033040004.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data in this study can be download from Gene4Denovo database (http://genemed.tech/gene4denovo/home (accessed on 6 March 2021)).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Willsey, A.J.; Morris, M.T.; Wang, S.; Willsey, H.R.; Sun, N.; Teerikorpi, N.; Baum, T.B.; Cagney, G.; Bender, K.J.; Desai, T.A.; et al. The Psychiatric Cell Map Initiative: A Convergent Systems Biological Approach to Illuminating Key Molecular Pathways in Neuropsychiatric Disorders. Cell 2018, 174, 505–520. [Google Scholar] [CrossRef] [Green Version]
  2. Heyne, H.O.; Singh, T.; Stamberger, H.; Abou Jamra, R.; Caglayan, H.; Craiu, D.; De Jonghe, P.; Guerrini, R.; Helbig, K.L.; Koeleman, B.P.C.; et al. De novo variants in neurodevelopmental disorders with epilepsy. Nat. Genet. 2018, 50, 1048–1053. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Brosig, C.L.; Bear, L.; Allen, S.; Hoffmann, R.G.; Pan, A.; Frommelt, M.; Mussatto, K.A. Preschool Neurodevelopmental Outcomes in Children with Congenital Heart Disease. J. Pediatr. 2017, 183, 80–86.e1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Latal, B. Neurodevelopmental Outcomes of the Child with Congenital Heart Disease. Clin. Perinatol. 2016, 43, 173–185. [Google Scholar] [CrossRef]
  5. Homsy, J.; Zaidi, S.; Shen, Y.; Ware, J.S.; Samocha, K.E.; Karczewski, K.J.; DePalma, S.R.; McKean, D.; Wakimoto, H.; Gorham, J.; et al. De novo mutations in congenital heart disease with neurodevelopmental and other congenital anomalies. Science 2015, 350, 1262–1266. [Google Scholar] [CrossRef] [Green Version]
  6. Marino, B.S.; Lipkin, P.H.; Newburger, J.W.; Peacock, G.; Gerdes, M.; Gaynor, J.W.; Mussatto, K.A.; Uzark, K.; Goldberg, C.S.; Johnson, W.H., Jr.; et al. Neurodevelopmental outcomes in children with congenital heart disease: Evaluation and management: A scientific statement from the American Heart Association. Circulation 2012, 126, 1143–1172. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Iossifov, I.; O’Roak, B.J.; Sanders, S.J.; Ronemus, M.; Krumm, N.; Levy, D.; Stessman, H.A.; Witherspoon, K.T.; Vives, L.; Patterson, K.E.; et al. The contribution of de novo coding mutations to autism spectrum disorder. Nature 2014, 515, 216–221. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Jin, S.C.; Homsy, J.; Zaidi, S.; Lu, Q.; Morton, S.; DePalma, S.R.; Zeng, X.; Qi, H.; Chang, W.; Sierant, M.C.; et al. Contribution of rare inherited and de novo variants in 2871 congenital heart disease probands. Nat. Genet. 2017, 49, 1593–1601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Deciphering Developmental Disorders, S. Prevalence and architecture of de novo mutations in developmental disorders. Nature 2017, 542, 433–438. [Google Scholar] [CrossRef]
  10. Epi, K.C.; Epilepsy Phenome/Genome, P.; Allen, A.S.; Berkovic, S.F.; Cossette, P.; Delanty, N.; Dlugos, D.; Eichler, E.E.; Epstein, M.P.; Glauser, T.; et al. De novo mutations in epileptic encephalopathies. Nature 2013, 501, 217–221. [Google Scholar]
  11. Lelieveld, S.H.; Reijnders, M.R.; Pfundt, R.; Yntema, H.G.; Kamsteeg, E.J.; de Vries, P.; de Vries, B.B.; Willemsen, M.H.; Kleefstra, T.; Lohner, K.; et al. Meta-analysis of 2,104 trios provides support for 10 new genes for intellectual disability. Nat. Neurosci. 2016, 19, 1194–1196. [Google Scholar] [CrossRef] [PubMed]
  12. Fromer, M.; Pocklington, A.J.; Kavanagh, D.H.; Williams, H.J.; Dwyer, S.; Gormley, P.; Georgieva, L.; Rees, E.; Palta, P.; Ruderfer, D.M.; et al. De novo mutations in schizophrenia implicate synaptic networks. Nature 2014, 506, 179–184. [Google Scholar] [CrossRef] [Green Version]
  13. Bernier, R.; Golzio, C.; Xiong, B.; Stessman, H.A.; Coe, B.P.; Penn, O.; Witherspoon, K.; Gerdts, J.; Baker, C.; Vulto-van Silfhout, A.T.; et al. Disruptive CHD8 mutations define a subtype of autism early in development. Cell 2014, 158, 263–276. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Van Bon, B.W.; Coe, B.P.; Bernier, R.; Green, C.; Gerdts, J.; Witherspoon, K.; Kleefstra, T.; Willemsen, M.H.; Kumar, R.; Bosco, P.; et al. Disruptive de novo mutations of DYRK1A lead to a syndromic form of autism and ID. Mol. Psychiatry 2016, 21, 126–132. [Google Scholar] [CrossRef] [PubMed]
  15. Sobreira, N.; Schiettecatte, F.; Valle, D.; Hamosh, A. GeneMatcher: A matching tool for connecting investigators with an interest in the same gene. Hum. Mutat. 2015, 36, 928–930. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Guo, H.; Bettella, E.; Marcogliese, P.C.; Zhao, R.; Andrews, J.C.; Nowakowski, T.J.; Gillentine, M.A.; Hoekzema, K.; Wang, T.; Wu, H.; et al. Disruptive mutations in TANC2 define a neurodevelopmental syndrome associated with psychiatric disorders. Nat. Commun. 2019, 10, 4679. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Guo, H.; Li, Y.; Shen, L.; Wang, T.; Jia, X.; Liu, L.; Xu, T.; Ou, M.; Hoekzema, K.; Wu, H.; et al. Disruptive variants of CSDE1 associate with autism and interfere with neuronal development and synaptic transmission. Sci. Adv. 2019, 5, eaax2166. [Google Scholar] [CrossRef] [Green Version]
  18. Li, J.; Cai, T.; Jiang, Y.; Chen, H.; He, X.; Chen, C.; Li, X.; Shao, Q.; Ran, X.; Li, Z.; et al. Genes with de novo mutations are shared by four neuropsychiatric disorders discovered from NPdenovo database. Mol. Psychiatry 2016, 21, 290–297. [Google Scholar] [CrossRef] [Green Version]
  19. Gonzalez-Mantilla, A.J.; Moreno-De-Luca, A.; Ledbetter, D.H.; Martin, C.L. A Cross-Disorder Method to Identify Novel Candidate Genes for Developmental Brain Disorders. JAMA Psychiatry 2016, 73, 275–283. [Google Scholar] [CrossRef] [Green Version]
  20. Coe, B.P.; Stessman, H.A.F.; Sulovari, A.; Geisheker, M.R.; Bakken, T.E.; Lake, A.M.; Dougherty, J.D.; Lein, E.S.; Hormozdiari, F.; Bernier, R.A.; et al. Neurodevelopmental disease genes implicated by de novo mutation and copy number variation morbidity. Nat. Genet. 2019, 51, 106–116. [Google Scholar] [CrossRef]
  21. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef]
  22. Li, J.; Shi, L.; Zhang, K.; Zhang, Y.; Hu, S.; Zhao, T.; Teng, H.; Li, X.; Jiang, Y.; Ji, L.; et al. VarCards: An integrated genetic and clinical database for coding variants in the human genome. Nucleic Acids Res. 2018, 46, D1039–D1048. [Google Scholar] [CrossRef] [Green Version]
  23. Li, J.; Zhao, T.; Zhang, Y.; Zhang, K.; Shi, L.; Chen, Y.; Wang, X.; Sun, Z. Performance evaluation of pathogenicity-computation methods for missense variants. Nucleic Acids Res. 2018, 46, 7793–7804. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zhao, G.; Li, K.; Li, B.; Wang, Z.; Fang, Z.; Wang, X.; Zhang, Y.; Luo, T.; Zhou, Q.; Wang, L.; et al. Gene4Denovo: An integrated database and analytic platform for de novo mutations in humans. Nucleic Acids Res. 2020, 48, D913–D926. [Google Scholar] [CrossRef]
  25. Shohat, S.; Ben-David, E.; Shifman, S. Varying Intolerance of Gene Pathways to Mutational Classes Explain Genetic Convergence across Neuropsychiatric Disorders. Cell Rep. 2017, 18, 2217–2227. [Google Scholar] [CrossRef] [Green Version]
  26. He, X.; Sanders, S.J.; Liu, L.; De Rubeis, S.; Lim, E.T.; Sutcliffe, J.S.; Schellenberg, G.D.; Gibbs, R.A.; Daly, M.J.; Buxbaum, J.D.; et al. Integrated model of de novo and inherited genetic variants yields greater power to identify risk genes. PLoS Genet. 2013, 9, e1003671. [Google Scholar] [CrossRef] [Green Version]
  27. De Rubeis, S.; He, X.; Goldberg, A.P.; Poultney, C.S.; Samocha, K.; Cicek, A.E.; Kou, Y.; Liu, L.; Fromer, M.; Walker, S.; et al. Synaptic, transcriptional and chromatin genes disrupted in autism. Nature 2014, 515, 209–215. [Google Scholar] [CrossRef]
  28. Zhao, W.; Quan, Y.; Wu, H.; Han, L.; Bai, T.; Ma, L.; Li, B.; Xun, G.; Ou, J.; Zhao, J.; et al. POGZ de novo missense variants in neuropsychiatric disorders. Mol. Genet. Genomic Med. 2019, 7, e900. [Google Scholar] [CrossRef] [Green Version]
  29. Kaplanis, J.; Samocha, K.E.; Wiel, L.; Zhang, Z.; Arvai, K.J.; Eberhardt, R.Y.; Gallone, G.; Lelieveld, S.H.; Martin, H.C.; McRae, J.F.; et al. Evidence for 28 genetic disorders discovered by combining healthcare and research data. Nature 2020, 586, 757–762. [Google Scholar] [CrossRef]
  30. Satterstrom, F.K.; Kosmicki, J.A.; Wang, J.; Breen, M.S.; De Rubeis, S.; An, J.-Y.; Peng, M.; Collins, R.L.; Grove, J.; Klei, L.; et al. Novel genes for autism implicate both excitatory and inhibitory cell lineages in risk. bioRxiv 2018. [Google Scholar] [CrossRef]
  31. Ruzzo, E.K.; Perez-Cano, L.; Jung, J.Y.; Wang, L.K.; Kashef-Haghighi, D.; Hartl, C.; Singh, C.; Xu, J.; Hoekstra, J.N.; Leventhal, O.; et al. Inherited and De Novo Genetic Risk for Autism Impacts Shared Networks. Cell 2019, 178, 850–866.e26. [Google Scholar] [CrossRef] [Green Version]
  32. Abrahams, B.S.; Arking, D.E.; Campbell, D.B.; Mefford, H.C.; Morrow, E.M.; Weiss, L.A.; Menashe, I.; Wadkins, T.; Banerjee-Basu, S.; Packer, A. SFARI Gene 2.0: A community-driven knowledgebase for the autism spectrum disorders (ASDs). Mol. Autism 2013, 4, 36. [Google Scholar] [CrossRef] [Green Version]
  33. Satterstrom, F.K.; Kosmicki, J.A.; Wang, J.; Breen, M.S.; De Rubeis, S.; An, J.Y.; Peng, M.; Collins, R.; Grove, J.; Klei, L.; et al. Large-Scale Exome Sequencing Study Implicates Both Developmental and Functional Changes in the Neurobiology of Autism. Cell 2020, 180, 568–584.e23. [Google Scholar] [CrossRef]
  34. Doherty, J.L.; Owen, M.J. Genomic insights into the overlap between psychiatric disorders: Implications for research and clinical practice. Genome Med. 2014, 6, 29. [Google Scholar] [CrossRef] [Green Version]
  35. Song, N.N.; Ma, P.; Zhang, Q.; Zhang, L.; Wang, H.; Zhang, L.; Zhu, L.; He, C.H.; Mao, B.; Ding, Y.Q. Rnf220/Zc4h2-mediated monoubiquitylation of Phox2 is required for noradrenergic neuron development. Development 2020, 147, 6. [Google Scholar] [CrossRef]
  36. Kim, J.; Choi, T.I.; Park, S.; Kim, M.H.; Kim, C.H.; Lee, S. Rnf220 cooperates with Zc4h2 to specify spinal progenitor domains. Development 2018, 145, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Crawley, J.N.; Heyer, W.D.; LaSalle, J.M. Autism and Cancer Share Risk Genes, Pathways, and Drug Targets. Trends Genet. 2016, 32, 139–146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Butler, M.G.; Dasouki, M.J.; Zhou, X.P.; Talebizadeh, Z.; Brown, M.; Takahashi, T.N.; Miles, J.H.; Wang, C.H.; Stratton, R.; Pilarski, R.; et al. Subset of individuals with autism spectrum disorders and extreme macrocephaly associated with germline PTEN tumour suppressor gene mutations. J. Med. Genet. 2005, 42, 318–321. [Google Scholar] [CrossRef] [Green Version]
  39. Wang, L.; Zhou, K.; Fu, Z.; Yu, D.; Huang, H.; Zang, X.; Mo, X. Brain Development and Akt Signaling: The Crossroads of Signaling Pathway and Neurodevelopmental Diseases. J. Mol. Neurosci. 2017, 61, 379–384. [Google Scholar] [CrossRef] [Green Version]
  40. Madsen, R.R.; Vanhaesebroeck, B.; Semple, R.K. Cancer-Associated PIK3CA Mutations in Overgrowth Disorders. Trends Mol. Med. 2018, 24, 856–870. [Google Scholar] [CrossRef] [Green Version]
  41. Ma, P.; Yang, X.; Kong, Q.; Li, C.; Yang, S.; Li, Y.; Mao, B. The ubiquitin ligase RNF220 enhances canonical Wnt signaling through USP7-mediated deubiquitination of beta-catenin. Mol. Cell Biol. 2014, 34, 4355–4366. [Google Scholar] [CrossRef] [Green Version]
  42. Li, B.; Li, K.; Tian, D.; Zhou, Q.; Xie, Y.; Fang, Z.; Wang, X.; Luo, T.; Wang, Z.; Zhang, Y.; et al. De novo mutation of cancer-related genes associates with particular neurodevelopmental disorders. J. Mol. Med. 2020, 98, 1701–1712. [Google Scholar] [CrossRef] [PubMed]
  43. Dang, X.; Qin, Y.; Gu, C.; Sun, J.; Zhang, R.; Peng, Z. Knockdown of Tripartite Motif 8 Protects H9C2 Cells Against Hypoxia/Reoxygenation-Induced Injury Through the Activation of PI3K/Akt Signaling Pathway. Cell Transpl. 2020, 29, 963689720949247. [Google Scholar] [CrossRef]
  44. Wu, Y.; Tan, X.; Liu, P.; Yang, Y.; Huang, Y.; Liu, X.; Meng, X.; Yu, B.; Wu, M.; Jin, H. ITGA6 and RPSA synergistically promote pancreatic cancer invasion and metastasis via PI3K and MAPK signaling pathways. Exp. Cell Res. 2019, 379, 30–47. [Google Scholar] [CrossRef]
  45. Chou, S.J.; O’Leary, D.D. Role for Lhx2 in corticogenesis through regulation of progenitor differentiation. Mol. Cell. Neurosci. 2013, 56, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Hou, P.S.; Chuang, C.Y.; Kao, C.F.; Chou, S.J.; Stone, L.; Ho, H.N.; Chien, C.L.; Kuo, H.C. LHX2 regulates the neural differentiation of human embryonic stem cells via transcriptional modulation of PAX6 and CER1. Nucleic Acids Res. 2013, 41, 7753–7770. [Google Scholar] [CrossRef] [Green Version]
  47. Takahashi, M.; Tsukamoto, Y.; Kai, T.; Tokunaga, A.; Nakada, C.; Hijiya, N.; Uchida, T.; Daa, T.; Nomura, T.; Sato, F.; et al. Downregulation of WDR20 due to loss of 14q is involved in the malignant transformation of clear cell renal cell carcinoma. Cancer Sci. 2016, 107, 417–423. [Google Scholar] [CrossRef] [Green Version]
  48. Guerrini-Rousseau, L.; Dufour, C.; Varlet, P.; Masliah-Planchon, J.; Bourdeaut, F.; Guillaud-Bataille, M.; Abbas, R.; Bertozzi, A.I.; Fouyssac, F.; Huybrechts, S.; et al. Germline SUFU mutation carriers and medulloblastoma: Clinical characteristics, cancer risk, and prognosis. Neuro Oncol. 2018, 20, 1122–1132. [Google Scholar] [CrossRef] [Green Version]
  49. Korshunov, A.; Sahm, F.; Stichel, D.; Schrimpf, D.; Ryzhova, M.; Zheludkova, O.; Golanov, A.; Lichter, P.; Jones, D.T.W.; von Deimling, A.; et al. Molecular characterization of medulloblastomas with extensive nodularity (MBEN). Acta Neuropathol. 2018, 136, 303–313. [Google Scholar] [CrossRef]
  50. Cappi, C.; Oliphant, M.E.; Peter, Z.; Zai, G.; Conceicao do Rosario, M.; Sullivan, C.A.W.; Gupta, A.R.; Hoffman, E.J.; Virdee, M.; Olfson, E.; et al. De Novo Damaging DNA Coding Mutations Are Associated With Obsessive-Compulsive Disorder and Overlap With Tourette’s Disorder and Autism. Biol. Psychiatry 2020, 87, 1035–1044. [Google Scholar] [CrossRef]
  51. Willsey, A.J.; Fernandez, T.V.; Yu, D.; King, R.A.; Dietrich, A.; Xing, J.; Sanders, S.J.; Mandell, J.D.; Huang, A.Y.; Richer, P.; et al. De Novo Coding Variants Are Strongly Associated with Tourette Disorder. Neuron 2017, 94, 486–499.e9. [Google Scholar] [CrossRef] [Green Version]
  52. Guo, H.; Wang, T.; Wu, H.; Long, M.; Coe, B.P.; Li, H.; Xun, G.; Ou, J.; Chen, B.; Duan, G.; et al. Inherited and multiple de novo mutations in autism/developmental delay risk genes suggest a multifactorial model. Mol. Autism 2018, 9, 64. [Google Scholar] [CrossRef]
  53. Girirajan, S.; Rosenfeld, J.A.; Cooper, G.M.; Antonacci, F.; Siswara, P.; Itsara, A.; Vives, L.; Walsh, T.; McCarthy, S.E.; Baker, C.; et al. A recurrent 16p12.1 microdeletion supports a two-hit model for severe developmental delay. Nat. Genet. 2010, 42, 203–209. [Google Scholar] [CrossRef] [PubMed]
  54. Schaaf, C.P.; Sabo, A.; Sakai, Y.; Crosby, J.; Muzny, D.; Hawes, A.; Lewis, L.; Akbar, H.; Varghese, R.; Boerwinkle, E.; et al. Oligogenic heterozygosity in individuals with high-functioning autism spectrum disorders. Hum. Mol. Genet. 2011, 20, 3366–3375. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Du, Y.; Li, Z.; Liu, Z.; Zhang, N.; Wang, R.; Li, F.; Zhang, T.; Jiang, Y.; Zhi, X.; Wang, Z.; et al. Nonrandom occurrence of multiple de novo coding variants in a proband indicates the existence of an oligogenic model in autism. Genet Med. 2020, 22, 170–180. [Google Scholar] [CrossRef] [PubMed]
  56. Gifford, C.A.; Ranade, S.S.; Samarakoon, R.; Salunga, H.T.; de Soysa, T.Y.; Huang, Y.; Zhou, P.; Elfenbein, A.; Wyman, S.K.; Bui, Y.K.; et al. Oligogenic inheritance of a human heart disease involving a genetic modifier. Science 2019, 364, 865–870. [Google Scholar] [CrossRef] [PubMed]
  57. Jaganathan, K.; Kyriazopoulou Panagiotopoulou, S.; McRae, J.F.; Darbandi, S.F.; Knowles, D.; Li, Y.I.; Kosmicki, J.A.; Arbelaez, J.; Cui, W.; Schwartz, G.B.; et al. Predicting Splicing from Primary Sequence with Deep Learning. Cell 2019, 176, 535–548.e24. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Nguyen, H.T.; Bryois, J.; Kim, A.; Dobbyn, A.; Huckins, L.M.; Munoz-Manchado, A.B.; Ruderfer, D.M.; Genovese, G.; Fromer, M.; Xu, X.; et al. Integrated Bayesian analysis of rare exonic variants to identify risk genes for schizophrenia and neurodevelopmental disorders. Genome Med. 2017, 9, 114. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Life 11 00233 g001 550
Figure 1. Genetic similarity between different neurodevelopmental disorders. Genetic similarity among disorders were performed based on three classes of variants include LoF, Dmis and Pfun. OE, ratio of observed to expected numbers of shared genes. Solid and coloured circle indicate OE greater than 1 and p value less than 0.05. Solid circle with no colour indicate OE greater than 1 but p value great than 0.05. Solid circle with no colour indicate OE greater than 1. Hollow circle indicate that OE less than 1. Dmis, Deleterious missense variants; LoF, loss of function. LoF include frameshift, stoploss and stopgain, splicing variants. Pfun, Putative functional variants, including Dmis and LoF variants. p value was calculated by using DNENRICH software (v1.0). ASD, autism spectrum disorder; SCZ, schizophrenia; EE, epileptic encephalopathy; DD/ID, developmental disorders/intellectual disability; CHD, congenital heart disease, TD, Tourette disorder; BP, bipolar disorder; OCD, obsessive-compulsive disorder; CMS, complex motor stereotypies.
Figure 1. Genetic similarity between different neurodevelopmental disorders. Genetic similarity among disorders were performed based on three classes of variants include LoF, Dmis and Pfun. OE, ratio of observed to expected numbers of shared genes. Solid and coloured circle indicate OE greater than 1 and p value less than 0.05. Solid circle with no colour indicate OE greater than 1 but p value great than 0.05. Solid circle with no colour indicate OE greater than 1. Hollow circle indicate that OE less than 1. Dmis, Deleterious missense variants; LoF, loss of function. LoF include frameshift, stoploss and stopgain, splicing variants. Pfun, Putative functional variants, including Dmis and LoF variants. p value was calculated by using DNENRICH software (v1.0). ASD, autism spectrum disorder; SCZ, schizophrenia; EE, epileptic encephalopathy; DD/ID, developmental disorders/intellectual disability; CHD, congenital heart disease, TD, Tourette disorder; BP, bipolar disorder; OCD, obsessive-compulsive disorder; CMS, complex motor stereotypies.
Life 11 00233 g001
Life 11 00233 g002 550
Figure 2. Projected gene discovery in larger cohort size. We assume the sample size were 500, 1000, 2000, 4000, 8000, 16,000, and 32,000 for each disorder and sampling de novo variant from exist based on putative functional de novo variant rate per individual. We then estimate the number of candidate gene for each disorders with FDR < 0.05 by transmitted and de novo association (TADA) analysis. ASD, autism spectrum disorder; SCZ, schizophrenia; EE, epileptic encephalopathy; DD/ID, developmental disorders/intellectual disability; CHD, congenital heart disease, TD, Tourette disorder; BP, bipolar disorder; OCD, obsessive-compulsive disorder; CMS, complex motor stereotypies.
Figure 2. Projected gene discovery in larger cohort size. We assume the sample size were 500, 1000, 2000, 4000, 8000, 16,000, and 32,000 for each disorder and sampling de novo variant from exist based on putative functional de novo variant rate per individual. We then estimate the number of candidate gene for each disorders with FDR < 0.05 by transmitted and de novo association (TADA) analysis. ASD, autism spectrum disorder; SCZ, schizophrenia; EE, epileptic encephalopathy; DD/ID, developmental disorders/intellectual disability; CHD, congenital heart disease, TD, Tourette disorder; BP, bipolar disorder; OCD, obsessive-compulsive disorder; CMS, complex motor stereotypies.
Life 11 00233 g002
Table
Table 1. Summary of collected DNVs in neurodevelopmental disorders.
Table 1. Summary of collected DNVs in neurodevelopmental disorders.
PhenotypesStudyTriosDNVsCoding DNVsPTVsDmisPfunPfun per Individual
ASD1410,318287,44412,1411580250740870.40
SCZa1134023422335735871610740.32
EE9973124811911703645340.55
DD/ID631,26045,54144,825707811,68318,7610.60
CHD126452990298136965410230.39
TD3909842818851992840.31
BPa321969951993421550.25
OCD11181341284820680.58
CMS11842051982712390.21
ASD, autism spectrum disorder; SCZ, schizophrenia; EE, epileptic encephalopathy; DD/ID, developmental disorders/intellectual disability; CHD, congenital heart disease, TD, tourette disorder; BP, bipolar disorder; OCD, obsessive-compulsive disorder; CMS, complex motor stereotypies; DNVs, de novo variants; PTVs, protein-truncating variants; Dmis, deleterious missense variant; Pfun, putative functional variant, combining PTVs and Dmis. a, several patients with SCZ/BP come from one study.
Table
Table 2. Comparison of prioritised candidate gene number by integrated analysis based on mutation type.
Table 2. Comparison of prioritised candidate gene number by integrated analysis based on mutation type.
Disorders (N)Genetic SimilarityCategoryTypeFDR < 0.00010.0001 < FDR < 0.0010.001 < FDR < 0.010.01 < FDR < 0.05
p-ValueOE
ASD (10,318)1.00 × 10−43.72Before 2472350
AfterPfun229314755
LoF141162430
Dmis175213133
SCZ (3402)1.00 × 10−41.67Before 0035
AfterPfun6891718
LoF291813
Dmis479116
EE (973)1.00 × 10−46.54Before 7458
AfterPfun876106
LoF38192
Dmis58515
DD/ID (31,260)1.00 × 10−46.80Before 2785381115
AfterPfun287567996
LoF237466465
Dmis267507073
CHD (2645)1.00 × 10−42.84Before 33412
AfterPfun78141620
LoF456814
Dmis4681110
TD (909)1.00 × 10−42.02Before 0000
AfterPfun21160
LoF7020
Dmis14140
BP (219)2.89 × 10−21.90Before 0000
AfterPfun3110
LoF2000
Dmis2110
OCD (118)2.00 × 10−43.20Before 0000
AfterPfun10130
LoF2000
Dmis9130
CMS (184)1.00 × 10−22.49Before 0001
AfterPfun4113
LoF3000
Dmis1113
ASD, autism spectrum disorder; SCZ, schizophrenia; DD/ID, developmental disorders/intellectual disability; CHD, congenital heart disease; TD, Tourette disorder; BP, bipolar disorder; OCD, obsessive-compulsive disorder; CMS, complex motor stereotypies; Pfun, putative functional variant; LoF, loss of function variant; Dmis, deleterious missense variant; Before, prioritised candidate gene base on putative functional DNVs of specific disorder with FDR < 0.05; After, prioritised candidate gene base on the integration of DNVs in all disorders. Gene carrying Pfun, LoF and Dmis in specific disorder and pass each FDR threshold in integration analysis was defined as candidate gene of this disorder. OE, ratio of observed to expected numbers of shared genes with putative functional de novo variants.
Table
Table 3. Candidate gene carrying putative functional variants in different number of disorders (FDR < 0.05).
Table 3. Candidate gene carrying putative functional variants in different number of disorders (FDR < 0.05).
Rank (FDR)Unique Disorders
n = 173, 26.45%		Two Disorders
n = 239, 36.54%		Three Disorders
n = 167, 25.54%		Four Disorders
n = 59, 9.02%		Five Disorders
n = 13, 1.99%		Six Disorders
n = 3, 0.46%
[0, 0.0001) (48.32%)421139850103
[0.0001, 0.001) (9.17%)142618200
[0.001, 0.01) (15.44%)314121620
[0.01, 0.05) (27.06%)865930110
Candidate genes are split into six parts based on the number of disorders with putative functional DNMs in specific gene. Unique genes means gene only carry putative functional DNMs in one disorder. We ranked all candidate genes into four tiers based on the strength of FDR.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, K.; Ling, Z.; Luo, T.; Zhao, G.; Zhou, Q.; Wang, X.; Xia, K.; Li, J.; Li, B. Cross-Disorder Analysis of De Novo Variants Increases the Power of Prioritising Candidate Genes. Life 2021, 11, 233. https://doi.org/10.3390/life11030233

AMA Style

Li K, Ling Z, Luo T, Zhao G, Zhou Q, Wang X, Xia K, Li J, Li B. Cross-Disorder Analysis of De Novo Variants Increases the Power of Prioritising Candidate Genes. Life. 2021; 11(3):233. https://doi.org/10.3390/life11030233

Chicago/Turabian Style

Li, Kuokuo, Zhengbao Ling, Tengfei Luo, Guihu Zhao, Qiao Zhou, Xiaomeng Wang, Kun Xia, Jinchen Li, and Bin Li. 2021. "Cross-Disorder Analysis of De Novo Variants Increases the Power of Prioritising Candidate Genes" Life 11, no. 3: 233. https://doi.org/10.3390/life11030233

APA Style

Li, K., Ling, Z., Luo, T., Zhao, G., Zhou, Q., Wang, X., Xia, K., Li, J., & Li, B. (2021). Cross-Disorder Analysis of De Novo Variants Increases the Power of Prioritising Candidate Genes. Life, 11(3), 233. https://doi.org/10.3390/life11030233

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop