ABSTRACT
A large number of putative risk genes for autism spectrum disorder (ASD) have been reported. The functions of most of these susceptibility genes in developing brains remain unknown, and causal relationships between their variation and autism traits have not been established. The aim of this study was to predict putative risk genes at the whole-genome level based on the analysis of gene co-expression with a group of high-confidence ASD risk genes (hcASDs). The results showed that three gene features - gene size, mRNA abundance, and guanine-cytosine content - affect the genome-wide co-expression profiles of hcASDs. To circumvent the interference of these features in gene co-expression analysis, we developed a method to determine whether a gene is significantly co-expressed with hcASDs by statistically comparing the co-expression profile of this gene with hcASDs to that of this gene with permuted gene sets of feature-matched genes. This method is referred to as "matched-gene co-expression analysis" (MGCA). With MGCA, we demonstrated the convergence in developmental expression profiles of hcASDs and improved the efficacy of risk gene prediction. The results of analysis of two recently-reported ASD candidate genes, CDH11 and CDH9, suggested the involvement of CDH11, but not CDH9, in ASD. Consistent with this prediction, behavioral studies showed that Cdh11-null mice, but not Cdh9-null mice, have multiple autism-like behavioral alterations. This study highlights the power of MGCA in revealing ASD-associated genes and the potential role of CDH11 in ASD.
Subject(s)
Animals , Mice , Autism Spectrum Disorder/genetics , Brain , Cadherins/genetics , Gene Expression , Mice, KnockoutABSTRACT
Three prevalent SARS-CoV-2 Variants of Concern (VOCs) were emerged and caused epidemic waves. It is essential to uncover the key genetic changes that cause the high transmissibility of VOCs. However, different viral mutations are generally tightly linked so traditional population genetic methods may not reliably detect beneficial mutation. In this study, we proposed a new pandemic-scale phylogenomic approach to detect mutations crucial to transmissibility. We analyzed 3,646,973 high-quality SARS-CoV-2 genomic sequences and the epidemiology metadata. Based on the sequential occurrence order of mutations and the instantaneously accelerated furcation rate, the analysis revealed that two non-coding mutations at the position of 28271 (g.a28271-/t) might be crucial for the high transmissibility of Alpha, Delta and Omicron VOCs. Both two mutations cause an A-to-T change at the core Kozak site of the N gene. The analysis also revealed that the non-coding mutations (g.a28271-/t) alone are unlikely to cause high viral transmissibility, indicating epistasis or multilocus interaction in viral transmissibility. A convergent evolutionary analysis revealed that g.a28271-/t, S:P681H/R and N:R203K/M occur independently in the three-VOC lineages, suggesting a potential interaction among these mutations. Therefore, this study unveils that non-synonymous and non-coding mutations could affect the transmissibility synergistically.
ABSTRACT
Genomic epidemiology is important to study the COVID-19 pandemic and more than two million SARS-CoV-2 genomic sequences were deposited into public databases. However, the exponential increase of sequences invokes unprecedented bioinformatic challenges. Here, we present the Coronavirus GenBrowser (CGB) based on a highly efficient analysis framework and a movie maker strategy. In total, 1,002,739 high quality genomic sequences with the transmission-related metadata were analyzed and visualized. The size of the core data file is only 12.20 MB, efficient for clean data sharing. Quick visualization modules and rich interactive operations are provided to explore the annotated SARS-CoV-2 evolutionary tree. CGB binary nomenclature is proposed to name each internal lineage. The pre-analyzed data can be filtered out according to the user-defined criteria to explore the transmission of SARS-CoV-2. Different evolutionary analyses can also be easily performed, such as the detection of accelerated evolution and on-going positive selection. Moreover, the 75 genomic spots conserved in SARS-CoV-2 but non-conserved in other coronaviruses were identified, which may indicate the functional elements specifically important for SARS-CoV-2. The CGB not only enables users who have no programming skills to analyze millions of genomic sequences, but also offers a panoramic vision of the transmission and evolution of SARS-CoV-2.
