1
|
Gonzalez-Ferrer J, Lehrer J, O'Farrell A, Paten B, Teodorescu M, Haussler D, Jonsson VD, Mostajo-Radji MA. SIMS: A deep-learning label transfer tool for single-cell RNA sequencing analysis. CELL GENOMICS 2024; 4:100581. [PMID: 38823397 PMCID: PMC11228957 DOI: 10.1016/j.xgen.2024.100581] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/17/2023] [Revised: 04/02/2024] [Accepted: 05/09/2024] [Indexed: 06/03/2024]
Abstract
Cell atlases serve as vital references for automating cell labeling in new samples, yet existing classification algorithms struggle with accuracy. Here we introduce SIMS (scalable, interpretable machine learning for single cell), a low-code data-efficient pipeline for single-cell RNA classification. We benchmark SIMS against datasets from different tissues and species. We demonstrate SIMS's efficacy in classifying cells in the brain, achieving high accuracy even with small training sets (<3,500 cells) and across different samples. SIMS accurately predicts neuronal subtypes in the developing brain, shedding light on genetic changes during neuronal differentiation and postmitotic fate refinement. Finally, we apply SIMS to single-cell RNA datasets of cortical organoids to predict cell identities and uncover genetic variations between cell lines. SIMS identifies cell-line differences and misannotated cell lineages in human cortical organoids derived from different pluripotent stem cell lines. Altogether, we show that SIMS is a versatile and robust tool for cell-type classification from single-cell datasets.
Collapse
Affiliation(s)
- Jesus Gonzalez-Ferrer
- Genomics Institute, University of California, Santa Cruz, Santa Cruz, CA 95060, USA; Live Cell Biotechnology Discovery Lab, University of California, Santa Cruz, Santa Cruz, CA 95060, USA; Department of Biomolecular Engineering, University of California, Santa Cruz, Santa Cruz, CA 95060, USA
| | - Julian Lehrer
- Genomics Institute, University of California, Santa Cruz, Santa Cruz, CA 95060, USA; Live Cell Biotechnology Discovery Lab, University of California, Santa Cruz, Santa Cruz, CA 95060, USA; Department of Applied Mathematics, University of California, Santa Cruz, Santa Cruz, CA 95060, USA
| | - Ash O'Farrell
- Genomics Institute, University of California, Santa Cruz, Santa Cruz, CA 95060, USA
| | - Benedict Paten
- Genomics Institute, University of California, Santa Cruz, Santa Cruz, CA 95060, USA; Department of Biomolecular Engineering, University of California, Santa Cruz, Santa Cruz, CA 95060, USA
| | - Mircea Teodorescu
- Genomics Institute, University of California, Santa Cruz, Santa Cruz, CA 95060, USA; Department of Biomolecular Engineering, University of California, Santa Cruz, Santa Cruz, CA 95060, USA; Department of Electrical and Computer Engineering, University of California, Santa Cruz, Santa Cruz, CA 95060, USA
| | - David Haussler
- Genomics Institute, University of California, Santa Cruz, Santa Cruz, CA 95060, USA; Department of Biomolecular Engineering, University of California, Santa Cruz, Santa Cruz, CA 95060, USA
| | - Vanessa D Jonsson
- Department of Biomolecular Engineering, University of California, Santa Cruz, Santa Cruz, CA 95060, USA; Department of Applied Mathematics, University of California, Santa Cruz, Santa Cruz, CA 95060, USA.
| | - Mohammed A Mostajo-Radji
- Genomics Institute, University of California, Santa Cruz, Santa Cruz, CA 95060, USA; Live Cell Biotechnology Discovery Lab, University of California, Santa Cruz, Santa Cruz, CA 95060, USA.
| |
Collapse
|
2
|
Majane AC, Cridland JM, Blair LK, Begun DJ. Evolution and genetics of accessory gland transcriptome divergence between Drosophila melanogaster and D. simulans. Genetics 2024; 227:iyae039. [PMID: 38518250 PMCID: PMC11151936 DOI: 10.1093/genetics/iyae039] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2023] [Revised: 08/27/2023] [Accepted: 02/15/2024] [Indexed: 03/24/2024] Open
Abstract
Studies of allele-specific expression in interspecific hybrids have provided important insights into gene-regulatory divergence and hybrid incompatibilities. Many such investigations in Drosophila have used transcriptome data from complex mixtures of many tissues or from gonads, however, regulatory divergence may vary widely among species, sexes, and tissues. Thus, we lack sufficiently broad sampling to be confident about the general biological principles of regulatory divergence. Here, we seek to fill some of these gaps in the literature by characterizing regulatory evolution and hybrid misexpression in a somatic male sex organ, the accessory gland, in F1 hybrids between Drosophila melanogaster and D. simulans. The accessory gland produces seminal fluid proteins, which play an important role in male and female fertility and may be subject to adaptive divergence due to male-male or male-female interactions. We find that trans differences are relatively more abundant than cis, in contrast to most of the interspecific hybrid literature, though large effect-size trans differences are rare. Seminal fluid protein genes have significantly elevated levels of expression divergence and tend to be regulated through both cis and trans divergence. We find limited misexpression (over- or underexpression relative to both parents) in this organ compared to most other Drosophila studies. As in previous studies, male-biased genes are overrepresented among misexpressed genes and are much more likely to be underexpressed. ATAC-Seq data show that chromatin accessibility is correlated with expression differences among species and hybrid allele-specific expression. This work identifies unique regulatory evolution and hybrid misexpression properties of the accessory gland and suggests the importance of tissue-specific allele-specific expression studies.
Collapse
Affiliation(s)
- Alex C Majane
- Department of Evolution and Ecology, University of California, Davis, CA 95616, USA
| | - Julie M Cridland
- Department of Evolution and Ecology, University of California, Davis, CA 95616, USA
| | - Logan K Blair
- Department of Evolution and Ecology, University of California, Davis, CA 95616, USA
| | - David J Begun
- Department of Evolution and Ecology, University of California, Davis, CA 95616, USA
| |
Collapse
|
3
|
Peng J, Gu Y, Liu J, Yi H, Ruan D, Huang H, Shu Y, Zong Z, Wu R, Li H. Identification of SOCS3 and PTGS2 as new biomarkers for the diagnosis of gout by cross-species comprehensive analysis. Heliyon 2024; 10:e30020. [PMID: 38707281 PMCID: PMC11066387 DOI: 10.1016/j.heliyon.2024.e30020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2023] [Revised: 04/18/2024] [Accepted: 04/18/2024] [Indexed: 05/07/2024] Open
Abstract
Background Gout is the most common inflammatory arthritis in adults. Gout is an arthritic disease caused by the deposition of monosodium urate crystal (MSU) in the joints, which can lead to acute inflammation and damage adjacent tissue. Hyperuricemia is the main risk factor for MSU crystal deposition and gout. With the increasing burden of gout disease, the identification of potential biomarkers and novel targets for diagnosis is urgently needed. Methods For the analysis of this subject paper, we downloaded the human gout data set GSE160170 and the gout mouse model data set GSE190138 from the GEO database. To obtain the differentially expressed genes (DEGs), we intersected the two data sets. Using the cytohubba algorithm, we identified the key genes and enriched them through GO and KEGG. The gene expression trends of three subgroups (normal control group, intermittent gout group and acute gout attack group) were analyzed by Series Test of Cluster (STC) analysis, and the key genes were screened out, and the diagnostic effect was verified by ROC curve. The expression of key genes in dorsal root nerve and spinal cord of gout mice was analyzed. Finally, the clinical samples of normal control group, hyperuricemia group, intermittent gout group and acute gout attack group were collected, and the expression of key genes at protein level was verified by ELISA. Result We obtained 59 co-upregulated and 28 co-downregulated genes by comparing the DEGs between gout mouse model data set and human gout data set. 7 hub DEGs(IL1B, IL10, NLRP3, SOCS3, PTGS2) were screened out via Cytohubba algorithm. The results of both GO and KEGG enrichment analyses indicate that 7 hub genes play a significant role in regulating the inflammatory response, cytokine production in immune response, and the TNF signaling pathway. The most representative hub genes SOCS3 and PTGS2 were screened out by Series Test of Cluster, and ROC analysis results showed the AUC values were both up to 1.000. In addition, we found that PTGS2 expression was significantly elevated in the dorsal root ganglia and spinal cord in monosodium urate(MSU)-induced gout mouse model. The ELISA results revealed that the expression of SOCS3 and PTGS2 was notably higher in the acute gout attack and intermittent gout groups compared to the normal control group. This difference was statistically significant, indicating a clear distinction between the groups. Conclusion Through cross-species comprehensive analysis and experimental verification, SOCS3 and PTGS2 were proved to be new biomarkers for diagnosing gout and predicting disease progression.
Collapse
Affiliation(s)
- Jie Peng
- Department of Rheumatology and Immunology, the First Affiliated Hospital, Jiangxi Medical College, Nanchang University, 330006, Nanchang, China
- The Second Clinical Medical College, Jiangxi Medical College, Nanchang University, 330006, Nanchang, China
- Department of Sports Medicine, Huashan Hospital, Fudan University, 200040, Shanghai, China
| | - Yawen Gu
- Department of Rheumatology and Immunology, the First Affiliated Hospital, Jiangxi Medical College, Nanchang University, 330006, Nanchang, China
- The Second Clinical Medical College, Jiangxi Medical College, Nanchang University, 330006, Nanchang, China
| | - Jiang Liu
- Department of Gastrointestinal Surgery, the Second Affiliated Hospital, Jiangxi Medical College, Nanchang University, 1 MinDe Road, 330006, Nanchang, China
| | - Hao Yi
- Department of Gastrointestinal Surgery, the Second Affiliated Hospital, Jiangxi Medical College, Nanchang University, 1 MinDe Road, 330006, Nanchang, China
| | - Dong Ruan
- Department of Rheumatology and Immunology, the First Affiliated Hospital, Jiangxi Medical College, Nanchang University, 330006, Nanchang, China
- The Second Clinical Medical College, Jiangxi Medical College, Nanchang University, 330006, Nanchang, China
- Department of Rehabilitation Medicine, the Second Affiliated Hospital, Jiangxi Medical College, Nanchang University, 1 Minde Road, 330006, Nanchang, China
| | - Haoyu Huang
- Department of Gastrointestinal Surgery, the Second Affiliated Hospital, Jiangxi Medical College, Nanchang University, 1 MinDe Road, 330006, Nanchang, China
| | - Yuan Shu
- The Second Clinical Medical College, Jiangxi Medical College, Nanchang University, 330006, Nanchang, China
| | - Zhen Zong
- Department of Gastrointestinal Surgery, the Second Affiliated Hospital, Jiangxi Medical College, Nanchang University, 1 MinDe Road, 330006, Nanchang, China
| | - Rui Wu
- Department of Rheumatology and Immunology, the First Affiliated Hospital, Jiangxi Medical College, Nanchang University, 330006, Nanchang, China
| | - Hui Li
- Department of Rheumatology and Immunology, the First Affiliated Hospital, Jiangxi Medical College, Nanchang University, 330006, Nanchang, China
| |
Collapse
|
4
|
Powell G, Simon MM, Pulit S, Mallon AM, Lindgren CM. Genic constraint against nonsynonymous variation across the mouse genome. BMC Genomics 2023; 24:562. [PMID: 37736706 PMCID: PMC10514939 DOI: 10.1186/s12864-023-09637-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2023] [Accepted: 08/30/2023] [Indexed: 09/23/2023] Open
Abstract
BACKGROUND Selective constraint, the depletion of variation due to negative selection, provides insights into the functional impact of variants and disease mechanisms. However, its characterization in mice, the most commonly used mammalian model, remains limited. This study aims to quantify mouse gene constraint using a new metric called the nonsynonymous observed expected ratio (NOER) and investigate its relationship with gene function. RESULTS NOER was calculated using whole-genome sequencing data from wild mouse populations (Mus musculus sp and Mus spretus). Positive correlations were observed between mouse gene constraint and the number of associated knockout phenotypes, indicating stronger constraint on pleiotropic genes. Furthermore, mouse gene constraint showed a positive correlation with the number of pathogenic variant sites in their human orthologues, supporting the relevance of mouse models in studying human disease variants. CONCLUSIONS NOER provides a resource for assessing the fitness consequences of genetic variants in mouse genes and understanding the relationship between gene constraint and function. The study's findings highlight the importance of pleiotropy in selective constraint and support the utility of mouse models in investigating human disease variants. Further research with larger sample sizes can refine constraint estimates in mice and enable more comprehensive comparisons of constraint between mouse and human orthologues.
Collapse
Affiliation(s)
- George Powell
- Li Ka Shing Centre for Health Information and Discovery, Big Data Institute, University of Oxford, Oxford, UK.
- Mammalian Genetics Unit, MRC Harwell Institute, Oxfordshire, OX11 0RD, UK.
| | - Michelle M Simon
- Mammalian Genetics Unit, MRC Harwell Institute, Oxfordshire, OX11 0RD, UK
| | - Sara Pulit
- Li Ka Shing Centre for Health Information and Discovery, Big Data Institute, University of Oxford, Oxford, UK
| | - Ann-Marie Mallon
- Mammalian Genetics Unit, MRC Harwell Institute, Oxfordshire, OX11 0RD, UK
| | - Cecilia M Lindgren
- Li Ka Shing Centre for Health Information and Discovery, Big Data Institute, University of Oxford, Oxford, UK
- Wellcome Centre for Human Genetics, University of Oxford, Oxford, UK
- Medical and Population Genetics Program, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| |
Collapse
|
5
|
Niu H, Pang Y, Xie L, Yu Q, Shen Y, Li J, Xu X. Clustering pattern and evolution characteristic of microRNAs in grass carp (Ctenopharyngodon idella). BMC Genomics 2023; 24:73. [PMID: 36782132 PMCID: PMC9926789 DOI: 10.1186/s12864-023-09159-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2022] [Accepted: 01/31/2023] [Indexed: 02/15/2023] Open
Abstract
BACKGROUND A considerable fraction of microRNAs (miRNAs) are highly conserved, and certain miRNAs correspond to genomic clusters. The clustering of miRNAs can be advantageous, possibly by allowing coordinated expression. However, little is known about the evolutionary forces responsible for the loss and acquisition of miRNA and miRNA clusters. RESULTS The results demonstrated that several novel miRNAs arose throughout grass carp evolution. Duplication and de novo production were critical strategies for miRNA cluster formation. Duplicates accounted for a smaller fraction of the expansion in the grass carp miRNA than de novo creation. Clustered miRNAs are more conserved and change slower, whereas unique miRNAs usually have high evolution rates and low expression levels. The expression level of miRNA expression in clusters is strongly correlated. CONCLUSIONS This study examines the genomic distribution, evolutionary background, and expression regulation of grass carp miRNAs. Our findings provide novel insights into the genesis and development of miRNA clusters in teleost.
Collapse
Affiliation(s)
- Huiqin Niu
- grid.412514.70000 0000 9833 2433Key Laboratory of Freshwater Aquatic Genetic Resources Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, China ,grid.412514.70000 0000 9833 2433National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China ,grid.412514.70000 0000 9833 2433Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, China
| | - Yifan Pang
- grid.412514.70000 0000 9833 2433Key Laboratory of Freshwater Aquatic Genetic Resources Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, China ,grid.412514.70000 0000 9833 2433National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China ,grid.412514.70000 0000 9833 2433Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, China
| | - Lingli Xie
- grid.412514.70000 0000 9833 2433Key Laboratory of Freshwater Aquatic Genetic Resources Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, China ,grid.412514.70000 0000 9833 2433National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China ,grid.412514.70000 0000 9833 2433Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, China
| | - Qiaozhen Yu
- grid.412514.70000 0000 9833 2433Key Laboratory of Freshwater Aquatic Genetic Resources Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, China ,grid.412514.70000 0000 9833 2433National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China ,grid.412514.70000 0000 9833 2433Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, China
| | - Yubang Shen
- grid.412514.70000 0000 9833 2433Key Laboratory of Freshwater Aquatic Genetic Resources Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, China ,grid.412514.70000 0000 9833 2433National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China ,grid.412514.70000 0000 9833 2433Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, China
| | - Jiale Li
- Key Laboratory of Freshwater Aquatic Genetic Resources Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, China. .,National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China. .,Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, China.
| | - Xiaoyan Xu
- Key Laboratory of Freshwater Aquatic Genetic Resources Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, China. .,National Demonstration Center for Experimental Fisheries Science Education, Shanghai Ocean University, Shanghai, China. .,Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, China.
| |
Collapse
|
6
|
Kasianov AS, Klepikova AV, Mayorov AV, Buzanov GS, Logacheva MD, Penin AA. Interspecific comparison of gene expression profiles using machine learning. PLoS Comput Biol 2023; 19:e1010743. [PMID: 36626392 PMCID: PMC9879537 DOI: 10.1371/journal.pcbi.1010743] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2021] [Revised: 01/26/2023] [Accepted: 11/16/2022] [Indexed: 01/11/2023] Open
Abstract
Interspecific gene comparisons are the keystones for many areas of biological research and are especially important for the translation of knowledge from model organisms to economically important species. Currently they are hampered by the low resolution of methods based on sequence analysis and by the complex evolutionary history of eukaryotic genes. This is especially critical for plants, whose genomes are shaped by multiple whole genome duplications and subsequent gene loss. This requires the development of new methods for comparing the functions of genes in different species. Here, we report ISEEML (Interspecific Similarity of Expression Evaluated using Machine Learning)-a novel machine learning-based algorithm for interspecific gene classification. In contrast to previous studies focused on sequence similarity, our algorithm focuses on functional similarity inferred from the comparison of gene expression profiles. We propose novel metrics for expression pattern similarity-expression score (ES)-that is suitable for species with differing morphologies. As a proof of concept, we compare detailed transcriptome maps of Arabidopsis thaliana, the model species, Zea mays (maize) and Fagopyrum esculentum (common buckwheat), which are species that represent distant clades within flowering plants. The classifier resulted in an AUC of 0.91; under the ES threshold of 0.5, the specificity was 94%, and sensitivity was 72%.
Collapse
Affiliation(s)
- Artem S. Kasianov
- Institute for Information Transmission Problems of the Russian Academy of Sciences, Moscow, Russia
| | - Anna V. Klepikova
- Institute for Information Transmission Problems of the Russian Academy of Sciences, Moscow, Russia
| | - Alexey V. Mayorov
- Institute for Information Transmission Problems of the Russian Academy of Sciences, Moscow, Russia
| | | | - Maria D. Logacheva
- Institute for Information Transmission Problems of the Russian Academy of Sciences, Moscow, Russia
- Skolkovo Institute of Science and Technology, Moscow, Russia
| | - Aleksey A. Penin
- Institute for Information Transmission Problems of the Russian Academy of Sciences, Moscow, Russia
- * E-mail:
| |
Collapse
|
7
|
Church SH, Munro C, Dunn CW, Extavour CG. The evolution of ovary-biased gene expression in Hawaiian Drosophila. PLoS Genet 2023; 19:e1010607. [PMID: 36689550 PMCID: PMC9894553 DOI: 10.1371/journal.pgen.1010607] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2022] [Revised: 02/02/2023] [Accepted: 01/09/2023] [Indexed: 01/24/2023] Open
Abstract
With detailed data on gene expression accessible from an increasingly broad array of species, we can test the extent to which our developmental genetic knowledge from model organisms predicts expression patterns and variation across species. But to know when differences in gene expression across species are significant, we first need to know how much evolutionary variation in gene expression we expect to observe. Here we provide an answer by analyzing RNAseq data across twelve species of Hawaiian Drosophilidae flies, focusing on gene expression differences between the ovary and other tissues. We show that over evolutionary time, there exists a cohort of ovary specific genes that is stable and that largely corresponds to described expression patterns from laboratory model Drosophila species. Our results also provide a demonstration of the prediction that, as phylogenetic distance increases, variation between species overwhelms variation between tissue types. Using ancestral state reconstruction of expression, we describe the distribution of evolutionary changes in tissue-biased expression, and use this to identify gains and losses of ovary-biased expression across these twelve species. We then use this distribution to calculate the evolutionary correlation in expression changes between genes, and demonstrate that genes with known interactions in D. melanogaster are significantly more correlated in their evolution than genes with no or unknown interactions. Finally, we use this correlation matrix to infer new networks of genes that share evolutionary trajectories, and we present these results as a dataset of new testable hypotheses about genetic roles and interactions in the function and evolution of the Drosophila ovary.
Collapse
Affiliation(s)
- Samuel H Church
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, United States of America
- Current address: Department of Ecology and Evolutionary Biology, Yale University, New Haven, Connecticut, United States of America
| | - Catriona Munro
- Collège de France, PSL Research University, CNRS, Inserm, Center for Interdisciplinary Research in Biology, Paris, France
| | - Casey W Dunn
- Current address: Department of Ecology and Evolutionary Biology, Yale University, New Haven, Connecticut, United States of America
| | - Cassandra G Extavour
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, United States of America
- Department of Molecular and Cellular Biology, Harvard University, Cambridge, Massachusetts, United States of America
- Howard Hughes Medical Institute, Chevy Chase, Maryland, United States of America
| |
Collapse
|
8
|
Xu W, Ren H, Qi X, Zhang S, Yu Z, Xie J. Conserved hierarchical gene regulatory networks for drought and cold stress response in Myrica rubra. FRONTIERS IN PLANT SCIENCE 2023; 14:1155504. [PMID: 37123838 PMCID: PMC10140524 DOI: 10.3389/fpls.2023.1155504] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/31/2023] [Accepted: 03/20/2023] [Indexed: 05/03/2023]
Abstract
Stress response in plant is regulated by a large number of genes co-operating in diverse networks that serve multiple adaptive process. To understand how gene regulatory networks (GRNs) modulating abiotic stress responses, we compare the GRNs underlying drought and cold stresses using samples collected at 4 or 6 h intervals within 48 h in Chinese bayberry (Myrica rubra). We detected 7,583 and 8,840 differentially expressed genes (DEGs) under drought and cold stress respectively, which might be responsive to environmental stresses. Drought- and cold-responsive GRNs, which have been built according to the timing of transcription under both abiotic stresses, have a conserved trans-regulator and a common regulatory network. In both GRNs, basic helix-loop-helix family transcription factor (bHLH) serve as central nodes. MrbHLHp10 transcripts exhibited continuous increase in the two abiotic stresses and acts upstream regulator of ASCORBATE PEROXIDASE (APX) gene. To examine the potential biological functions of MrbHLH10, we generated a transgenic Arabidopsis plant that constitutively overexpresses the MrbHLH10 gene. Compared to wild-type (WT) plants, overexpressing transgenic Arabidopsis plants maintained higher APX activity and biomass accumulation under drought and cold stress. Consistently, RNAi plants had elevated susceptibility to both stresses. Taken together, these results suggested that MrbHLH10 mitigates abiotic stresses through the modulation of ROS scavenging.
Collapse
Affiliation(s)
- Weijie Xu
- State Key Laboratory for Managing Biotic and Chemical Treats to the Quality and Safety of Agro-Products, Institute of Horticulture, Institute of Agro-product Safety and Nutrition, Zhejiang Academy of Agricultural Sciences, Hangzhou, China
- State Key Laboratory of Tree Genetics and Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
| | - Haiying Ren
- State Key Laboratory for Managing Biotic and Chemical Treats to the Quality and Safety of Agro-Products, Institute of Horticulture, Institute of Agro-product Safety and Nutrition, Zhejiang Academy of Agricultural Sciences, Hangzhou, China
- State Key Laboratory for Managing Biotic and Chemical Treats to the Quality and Safety of Agro-products, Hangzhou, China
- Xianghu Lab., Hangzhou, China
- *Correspondence: Haiying Ren, ; Jianbo Xie,
| | - Xingjiang Qi
- State Key Laboratory for Managing Biotic and Chemical Treats to the Quality and Safety of Agro-Products, Institute of Horticulture, Institute of Agro-product Safety and Nutrition, Zhejiang Academy of Agricultural Sciences, Hangzhou, China
- State Key Laboratory for Managing Biotic and Chemical Treats to the Quality and Safety of Agro-products, Hangzhou, China
- Xianghu Lab., Hangzhou, China
| | - Shuwen Zhang
- State Key Laboratory for Managing Biotic and Chemical Treats to the Quality and Safety of Agro-Products, Institute of Horticulture, Institute of Agro-product Safety and Nutrition, Zhejiang Academy of Agricultural Sciences, Hangzhou, China
- State Key Laboratory for Managing Biotic and Chemical Treats to the Quality and Safety of Agro-products, Hangzhou, China
| | - Zheping Yu
- State Key Laboratory for Managing Biotic and Chemical Treats to the Quality and Safety of Agro-Products, Institute of Horticulture, Institute of Agro-product Safety and Nutrition, Zhejiang Academy of Agricultural Sciences, Hangzhou, China
- State Key Laboratory for Managing Biotic and Chemical Treats to the Quality and Safety of Agro-products, Hangzhou, China
| | - Jianbo Xie
- State Key Laboratory of Tree Genetics and Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- The Tree and Ornamental Plant Breeding and Biotechnology Laboratory of National Forestry and Grassland Administration, Beijing Forestry University, Beijing, China
- *Correspondence: Haiying Ren, ; Jianbo Xie,
| |
Collapse
|
9
|
Duan G, Wu G, Chen X, Tian D, Li Z, Sun Y, Du Z, Hao L, Song S, Gao Y, Xiao J, Zhang Z, Bao Y, Tang B, Zhao W. HGD: an integrated homologous gene database across multiple species. Nucleic Acids Res 2022; 51:D994-D1002. [PMID: 36318261 PMCID: PMC9825607 DOI: 10.1093/nar/gkac970] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2022] [Revised: 09/28/2022] [Accepted: 10/17/2022] [Indexed: 11/06/2022] Open
Abstract
Homology is fundamental to infer genes' evolutionary processes and relationships with shared ancestry. Existing homolog gene resources vary in terms of inferring methods, homologous relationship and identifiers, posing inevitable difficulties for choosing and mapping homology results from one to another. Here, we present HGD (Homologous Gene Database, https://ngdc.cncb.ac.cn/hgd), a comprehensive homologs resource integrating multi-species, multi-resources and multi-omics, as a complement to existing resources providing public and one-stop data service. Currently, HGD houses a total of 112 383 644 homologous pairs for 37 species, including 19 animals, 16 plants and 2 microorganisms. Meanwhile, HGD integrates various annotations from public resources, including 16 909 homologs with traits, 276 670 homologs with variants, 398 573 homologs with expression and 536 852 homologs with gene ontology (GO) annotations. HGD provides a wide range of omics gene function annotations to help users gain a deeper understanding of gene function.
Collapse
Affiliation(s)
| | | | - Xiaoning Chen
- National Genomics Data Center, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China,University of Chinese Academy of Sciences, Beijing 100049, China
| | - Dongmei Tian
- National Genomics Data Center, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China,CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China
| | - Zhaohua Li
- National Genomics Data Center, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China,University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yanling Sun
- National Genomics Data Center, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China,CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China
| | - Zhenglin Du
- National Genomics Data Center, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China,CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China
| | - Lili Hao
- National Genomics Data Center, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China,CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China
| | - Shuhui Song
- National Genomics Data Center, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China,CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China,University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yuan Gao
- National Genomics Data Center, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China,CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China,University of Chinese Academy of Sciences, Beijing 100049, China
| | - Jingfa Xiao
- National Genomics Data Center, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China,CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China,University of Chinese Academy of Sciences, Beijing 100049, China
| | - Zhang Zhang
- National Genomics Data Center, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China,CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China,University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yiming Bao
- National Genomics Data Center, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China,CAS Key Laboratory of Genome Sciences and Information, Beijing Institute of Genomics, Chinese Academy of Sciences and China National Center for Bioinformation, Beijing 100101, China,University of Chinese Academy of Sciences, Beijing 100049, China
| | - Bixia Tang
- Correspondence may also be addressed to Bixia Tang.
| | - Wenming Zhao
- To whom correspondence should be addressed. Tel: +86 1084097636; Fax: +86 1084097720;
| |
Collapse
|
10
|
Jiang W, Chen L. Tissue Specificity of Gene Expression Evolves Across Mammal Species. J Comput Biol 2022; 29:880-891. [PMID: 35776510 PMCID: PMC9464367 DOI: 10.1089/cmb.2021.0592] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Tissue specificity of gene expression sheds light on the tissue-selective manifestation of hereditary disease despite the same DNA across all tissues. The evolutionary path of such tissue specificity provides essential information about the tissue-specific function of genes and the validity of disease animal models. With recent improvements of the sequencing technology, more and more large-scale transcriptomics studies have been conducted among different species across multiple tissues. In this study, we exploit existing transcriptomics resources of humans, cynomolgus macaques, rats, mice, and dogs across 13 tissues. We find that although tissue specificity of homologous gene expression is largely well conserved across species, a total of 380 genes shift or are in the process of shifting their tissue specificity. The tissue-specificity-shifting genes are less conserved than those preserving their tissue specificity or housekeeping genes. Interestingly, tissue-specificity-shifting genes tend to be less conserved at the third codon positions, likely due to their relaxed synonymous codon usage bias. Moreover, compared with genes, cassette exons are more likely to shift their tissue specificity of splicing across the five species.
Collapse
Affiliation(s)
- Wei Jiang
- Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, California, USA
| | - Liang Chen
- Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, California, USA.,Address correspondence to: Dr. Liang Chen, Department of Quantitative and Computational Biology, University of Southern California, 1050 Childs Way, Los Angeles, CA 90089, USA
| |
Collapse
|
11
|
Abstract
Gene expression divergence through evolutionary processes is thought to be important for achieving programmed development in multicellular organisms. To test this premise in filamentous fungi, we investigated transcriptional profiles of 3,942 single-copy orthologous genes (SCOGs) in five related sordariomycete species that have morphologically diverged in the formation of their flask-shaped perithecia. We compared expression of the SCOGs to inferred gene expression levels of the most recent common ancestor of the five species, ranking genes from their largest increases to smallest increases in expression during perithecial development in each of the five species. We found that a large proportion of the genes that exhibited evolved increases in gene expression were important for normal perithecial development in Fusarium graminearum. Many of these genes were previously uncharacterized, encoding hypothetical proteins without any known functional protein domains. Interestingly, the developmental stages during which aberrant knockout phenotypes appeared largely coincided with the elevated expression of the deleted genes. In addition, we identified novel genes that affected normal perithecial development in Magnaporthe oryzae and Neurospora crassa, which were functionally and transcriptionally diverged from the orthologous counterparts in F. graminearum. Furthermore, comparative analysis of developmental transcriptomes and phylostratigraphic analysis suggested that genes encoding hypothetical proteins are generally young and transcriptionally divergent between related species. This study provides tangible evidence of shifts in gene expression that led to acquisition of novel function of orthologous genes in each lineage and demonstrates that several genes with hypothetical function are crucial for shaping multicellular fruiting bodies. IMPORTANCE The fungal class Sordariomycetes includes numerous important plant and animal pathogens. It also provides model systems for studying fungal fruiting body development, as its members develop fruiting bodies with a few well-characterized tissue types on common growth media and have rich genomic resources that enable comparative and functional analyses. To understand transcriptional divergence of key developmental genes between five related sordariomycete fungi, we performed targeted knockouts of genes inferred to have evolved significant upward shifts in expression. We found that many previously uncharacterized genes play indispensable roles at different stages of fruiting body development, which have undergone transcriptional activation in specific lineages. These novel genes are predicted to be phylogenetically young and tend to be involved in lineage- or species-specific function. Transcriptional activation of genes with unknown function seems to be more frequent than ever thought, which may be crucial for rapid adaption to changing environments for successful sexual reproduction.
Collapse
|
12
|
Wang S, García-Seisdedos D, Prakash A, Kundu DJ, Collins A, George N, Fexova S, Moreno P, Papatheodorou I, Jones AR, Vizcaíno JA. Integrated view and comparative analysis of baseline protein expression in mouse and rat tissues. PLoS Comput Biol 2022; 18:e1010174. [PMID: 35714157 PMCID: PMC9246241 DOI: 10.1371/journal.pcbi.1010174] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2021] [Revised: 06/30/2022] [Accepted: 05/05/2022] [Indexed: 11/19/2022] Open
Abstract
The increasingly large amount of proteomics data in the public domain enables, among other applications, the combined analyses of datasets to create comparative protein expression maps covering different organisms and different biological conditions. Here we have reanalysed public proteomics datasets from mouse and rat tissues (14 and 9 datasets, respectively), to assess baseline protein abundance. Overall, the aggregated dataset contained 23 individual datasets, including a total of 211 samples coming from 34 different tissues across 14 organs, comprising 9 mouse and 3 rat strains, respectively. In all cases, we studied the distribution of canonical proteins between the different organs. The number of canonical proteins per dataset ranged from 273 (tendon) and 9,715 (liver) in mouse, and from 101 (tendon) and 6,130 (kidney) in rat. Then, we studied how protein abundances compared across different datasets and organs for both species. As a key point we carried out a comparative analysis of protein expression between mouse, rat and human tissues. We observed a high level of correlation of protein expression among orthologs between all three species in brain, kidney, heart and liver samples, whereas the correlation of protein expression was generally slightly lower between organs within the same species. Protein expression results have been integrated into the resource Expression Atlas for widespread dissemination.
Collapse
Affiliation(s)
- Shengbo Wang
- European Molecular Biology Laboratory - European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
| | - David García-Seisdedos
- European Molecular Biology Laboratory - European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
- Open Targets, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
| | - Ananth Prakash
- European Molecular Biology Laboratory - European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
- Open Targets, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
| | - Deepti Jaiswal Kundu
- European Molecular Biology Laboratory - European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
| | - Andrew Collins
- Institute of Systems, Molecular and Integrative Biology, University of Liverpool, Liverpool, United Kingdom
| | - Nancy George
- European Molecular Biology Laboratory - European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
| | - Silvie Fexova
- European Molecular Biology Laboratory - European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
| | - Pablo Moreno
- European Molecular Biology Laboratory - European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
| | - Irene Papatheodorou
- European Molecular Biology Laboratory - European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
- Open Targets, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
| | - Andrew R. Jones
- Institute of Systems, Molecular and Integrative Biology, University of Liverpool, Liverpool, United Kingdom
| | - Juan Antonio Vizcaíno
- European Molecular Biology Laboratory - European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
- Open Targets, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom
| |
Collapse
|
13
|
Assis R. No Expression Divergence despite Transcriptional Interference between Nested Protein-Coding Genes in Mammals. Genes (Basel) 2021; 12:genes12091381. [PMID: 34573363 PMCID: PMC8467205 DOI: 10.3390/genes12091381] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2021] [Revised: 08/23/2021] [Accepted: 08/24/2021] [Indexed: 01/05/2023] Open
Abstract
Nested protein-coding genes accumulated throughout metazoan evolution, with early analyses of human and Drosophila microarray data indicating that this phenomenon was simply due to the presence of large introns. However, a recent study employing RNA-seq data uncovered evidence of transcriptional interference driving rapid expression divergence between Drosophila nested genes, illustrating that accurate expression estimation of overlapping genes can enhance detection of their relationships. Hence, here I apply an analogous approach to strand-specific RNA-seq data from human and mouse to revisit the role of transcriptional interference in the evolution of mammalian nested genes. A genomic survey reveals that whereas mammalian nested genes indeed accrued over evolutionary time, they are retained at lower frequencies than in Drosophila. Though several properties of mammalian nested genes align with observations in Drosophila and with expectations under transcriptional interference, contrary to both, their expression divergence is not statistically different from that between unnested genes, and also does not increase after nesting. Together, these results support the hypothesis that lower selection efficiencies limit rates of gene expression evolution in mammals, leading to their reliance on immediate eradication of deleterious nested genes to avoid transcriptional interference.
Collapse
Affiliation(s)
- Raquel Assis
- Department of Electrical Engineering and Computer Science, Institute for Human Health and Disease Intervention, Florida Atlantic University, Boca Raton, FL 33431, USA
| |
Collapse
|
14
|
Razban RM, Dasmeh P, Serohijos AWR, Shakhnovich EI. Avoidance of protein unfolding constrains protein stability in long-term evolution. Biophys J 2021; 120:2413-2424. [PMID: 33932438 PMCID: PMC8390877 DOI: 10.1016/j.bpj.2021.03.042] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2020] [Revised: 02/24/2021] [Accepted: 03/17/2021] [Indexed: 11/28/2022] Open
Abstract
Every amino acid residue can influence a protein's overall stability, making stability highly susceptible to change throughout evolution. We consider the distribution of protein stabilities evolutionarily permittable under two previously reported protein fitness functions: flux dynamics and misfolding avoidance. We develop an evolutionary dynamics theory and find that it agrees better with an extensive protein stability data set for dihydrofolate reductase orthologs under the misfolding avoidance fitness function rather than the flux dynamics fitness function. Further investigation with ribonuclease H data demonstrates that not any misfolded state is avoided; rather, it is only the unfolded state. At the end, we discuss how our work pertains to the universal protein abundance-evolutionary rate correlation seen across organisms' proteomes. We derive a closed-form expression relating protein abundance to evolutionary rate that captures Escherichia coli, Saccharomyces cerevisiae, and Homo sapiens experimental trends without fitted parameters.
Collapse
Affiliation(s)
- Rostam M Razban
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts
| | - Pouria Dasmeh
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts; Departement de Biochimie, Université de Montréal, Montreal, Quebec, Canada
| | | | - Eugene I Shakhnovich
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts.
| |
Collapse
|
15
|
Sivaprakasam B, Sadagopan P. Development of shiny dashboard application for “genome-wide association study on analysis of SNPs injected in Homo sapiens genome (snips-HsG)”. GENE REPORTS 2021. [DOI: 10.1016/j.genrep.2021.101033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
16
|
Fodoulian L, Tuberosa J, Rossier D, Boillat M, Kan C, Pauli V, Egervari K, Lobrinus JA, Landis BN, Carleton A, Rodriguez I. SARS-CoV-2 Receptors and Entry Genes Are Expressed in the Human Olfactory Neuroepithelium and Brain. iScience 2020; 23:101839. [PMID: 33251489 PMCID: PMC7685946 DOI: 10.1016/j.isci.2020.101839] [Citation(s) in RCA: 120] [Impact Index Per Article: 30.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2020] [Revised: 09/25/2020] [Accepted: 11/18/2020] [Indexed: 12/21/2022] Open
Abstract
Reports indicate an association between COVID-19 and anosmia, as well as the presence of SARS-CoV-2 virions in the olfactory bulb. To test whether the olfactory neuroepithelium may represent a target of the virus, we generated RNA-seq libraries from human olfactory neuroepithelia, in which we found substantial expression of the genes coding for the virus receptor angiotensin-converting enzyme-2 (ACE2) and for the virus internalization enhancer TMPRSS2. We analyzed a human olfactory single-cell RNA-seq dataset and determined that sustentacular cells, which maintain the integrity of olfactory sensory neurons, express ACE2 and TMPRSS2. ACE2 protein was highly expressed in a subset of sustentacular cells in human and mouse olfactory tissues. Finally, we found ACE2 transcripts in specific brain cell types, both in mice and humans. Sustentacular cells thus represent a potential entry door for SARS-CoV-2 in a neuronal sensory system that is in direct connection with the brain.
Collapse
Affiliation(s)
- Leon Fodoulian
- Department of Genetics and Evolution, Faculty of Sciences, University of Geneva, quai Ernest-Ansermet 30, 1211 Geneva, Switzerland
- Department of Basic Neurosciences, Faculty of Medicine, University of Geneva, 1 rue Michel-Servet, 1211 Geneva, Switzerland
| | - Joël Tuberosa
- Department of Genetics and Evolution, Faculty of Sciences, University of Geneva, quai Ernest-Ansermet 30, 1211 Geneva, Switzerland
| | - Daniel Rossier
- Department of Genetics and Evolution, Faculty of Sciences, University of Geneva, quai Ernest-Ansermet 30, 1211 Geneva, Switzerland
| | - Madlaina Boillat
- Department of Genetics and Evolution, Faculty of Sciences, University of Geneva, quai Ernest-Ansermet 30, 1211 Geneva, Switzerland
| | - Chenda Kan
- Department of Genetics and Evolution, Faculty of Sciences, University of Geneva, quai Ernest-Ansermet 30, 1211 Geneva, Switzerland
| | - Véronique Pauli
- Department of Genetics and Evolution, Faculty of Sciences, University of Geneva, quai Ernest-Ansermet 30, 1211 Geneva, Switzerland
| | - Kristof Egervari
- Service of Clinical Pathology, Department of Genetic Medicine, Geneva University Hospitals, Geneva, Switzerland
- Department of Pathology and Immunology, Faculty of Medicine, University of Geneva, 1 rue Michel-Servet, 1211 Geneva, Switzerland
| | - Johannes A. Lobrinus
- Service of Clinical Pathology, Department of Genetic Medicine, Geneva University Hospitals, Geneva, Switzerland
| | - Basile N. Landis
- Rhinology-Olfactology Unit, Department of Otorhinolaryngology, Head and Neck Surgery, Geneva University Hospitals, Geneva, Switzerland
| | - Alan Carleton
- Department of Basic Neurosciences, Faculty of Medicine, University of Geneva, 1 rue Michel-Servet, 1211 Geneva, Switzerland
| | - Ivan Rodriguez
- Department of Genetics and Evolution, Faculty of Sciences, University of Geneva, quai Ernest-Ansermet 30, 1211 Geneva, Switzerland
| |
Collapse
|
17
|
Fodoulian L, Tuberosa J, Rossier D, Boillat M, Kan C, Pauli V, Egervari K, Lobrinus JA, Landis BN, Carleton A, Rodriguez I. SARS-CoV-2 Receptors and Entry Genes Are Expressed in the Human Olfactory Neuroepithelium and Brain. iScience 2020; 23:101839. [PMID: 33251489 DOI: 10.1101/2020.03.31.013268] [Citation(s) in RCA: 48] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2020] [Revised: 09/25/2020] [Accepted: 11/18/2020] [Indexed: 05/23/2023] Open
Abstract
Reports indicate an association between COVID-19 and anosmia, as well as the presence of SARS-CoV-2 virions in the olfactory bulb. To test whether the olfactory neuroepithelium may represent a target of the virus, we generated RNA-seq libraries from human olfactory neuroepithelia, in which we found substantial expression of the genes coding for the virus receptor angiotensin-converting enzyme-2 (ACE2) and for the virus internalization enhancer TMPRSS2. We analyzed a human olfactory single-cell RNA-seq dataset and determined that sustentacular cells, which maintain the integrity of olfactory sensory neurons, express ACE2 and TMPRSS2. ACE2 protein was highly expressed in a subset of sustentacular cells in human and mouse olfactory tissues. Finally, we found ACE2 transcripts in specific brain cell types, both in mice and humans. Sustentacular cells thus represent a potential entry door for SARS-CoV-2 in a neuronal sensory system that is in direct connection with the brain.
Collapse
Affiliation(s)
- Leon Fodoulian
- Department of Genetics and Evolution, Faculty of Sciences, University of Geneva, quai Ernest-Ansermet 30, 1211 Geneva, Switzerland
- Department of Basic Neurosciences, Faculty of Medicine, University of Geneva, 1 rue Michel-Servet, 1211 Geneva, Switzerland
| | - Joël Tuberosa
- Department of Genetics and Evolution, Faculty of Sciences, University of Geneva, quai Ernest-Ansermet 30, 1211 Geneva, Switzerland
| | - Daniel Rossier
- Department of Genetics and Evolution, Faculty of Sciences, University of Geneva, quai Ernest-Ansermet 30, 1211 Geneva, Switzerland
| | - Madlaina Boillat
- Department of Genetics and Evolution, Faculty of Sciences, University of Geneva, quai Ernest-Ansermet 30, 1211 Geneva, Switzerland
| | - Chenda Kan
- Department of Genetics and Evolution, Faculty of Sciences, University of Geneva, quai Ernest-Ansermet 30, 1211 Geneva, Switzerland
| | - Véronique Pauli
- Department of Genetics and Evolution, Faculty of Sciences, University of Geneva, quai Ernest-Ansermet 30, 1211 Geneva, Switzerland
| | - Kristof Egervari
- Service of Clinical Pathology, Department of Genetic Medicine, Geneva University Hospitals, Geneva, Switzerland
- Department of Pathology and Immunology, Faculty of Medicine, University of Geneva, 1 rue Michel-Servet, 1211 Geneva, Switzerland
| | - Johannes A Lobrinus
- Service of Clinical Pathology, Department of Genetic Medicine, Geneva University Hospitals, Geneva, Switzerland
| | - Basile N Landis
- Rhinology-Olfactology Unit, Department of Otorhinolaryngology, Head and Neck Surgery, Geneva University Hospitals, Geneva, Switzerland
| | - Alan Carleton
- Department of Basic Neurosciences, Faculty of Medicine, University of Geneva, 1 rue Michel-Servet, 1211 Geneva, Switzerland
| | - Ivan Rodriguez
- Department of Genetics and Evolution, Faculty of Sciences, University of Geneva, quai Ernest-Ansermet 30, 1211 Geneva, Switzerland
| |
Collapse
|
18
|
Abstract
New species arise as the genomes of populations diverge. The developmental 'alarm clock' of speciation sounds off when sufficient divergence in genetic control of development leads hybrid individuals to infertility or inviability, the world awoken to the dawn of new species with intrinsic post-zygotic reproductive isolation. Some developmental stages will be more prone to hybrid dysfunction due to how molecular evolution interacts with the ontogenetic timing of gene expression. Considering the ontogeny of hybrid incompatibilities provides a profitable connection between 'evo-devo' and speciation genetics to better link macroevolutionary pattern, microevolutionary process, and molecular mechanisms. Here, we explore speciation alongside development, emphasizing their mutual dependence on genetic network features, fitness landscapes, and developmental system drift. We assess models for how ontogenetic timing of reproductive isolation can be predictable. Experiments and theory within this synthetic perspective can help identify new rules of speciation as well as rules in the molecular evolution of development.
Collapse
Affiliation(s)
- Asher D Cutter
- Department of Ecology & Evolutionary Biology, University of TorontoTorontoCanada
| | - Joanna D Bundus
- Department of Integrative Biology, University of Wisconsin – MadisonMadisonUnited States
| |
Collapse
|
19
|
Xu W, Zhao Y, Chen S, Xie J, Zhang D. Evolution and Functional Divergence of the Fructokinase Gene Family in Populus. FRONTIERS IN PLANT SCIENCE 2020; 11:484. [PMID: 32499793 PMCID: PMC7243158 DOI: 10.3389/fpls.2020.00484] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/27/2019] [Accepted: 03/31/2020] [Indexed: 05/06/2023]
Abstract
New kinase has emerged throughout evolution, but how new kinase evolve while maintaining their functions and acquiring new functions remains unclear. Fructokinase (FRK), the gateway kinase to fructose metabolism, plays essential roles in plant development, and stress tolerance. Here, we explored the evolution of FRK gene family in 20 plant species (from green algae to angiosperms) and their functional roles in Populus. We identified 125 putative FRK genes in the 20 plant species with an average of 6 members per species. Phylogenetic analysis separated these 125 genes into 8 clades including 3 conserved clades and 5 specific clades, the 5 of which only exist in green algae or angiosperms. Evolutionary analysis revealed that FRK genes in ancient land plants have the largest number of functional domains with the longest amino acid sequences, and the length of FRK genes became shorter during the transition to vascular plants. This was accompanied by loss, acquisition, and diversification of functional domains. In Populus, segmental duplication appears to be the main mechanism for the expansion of FRK genes. Specially, most FRK genes duplicated in salicoids are regulated by Populus-specific microRNAs. Furthermore, compared with common FRKs, Populus-specific FRKs have showed higher expression specificity and are associated with fewer growth and wood property traits, which suggests that these FRKs may have undergone functional divergence. Our study explores the specific roles of FRKs in the Populus genome and provides new insights for functional investigation of this gene family.
Collapse
Affiliation(s)
- Weijie Xu
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, College of Biological Sciences and Technology, Beijing Forestry University, Ministry of Education, Beijing, China
| | - Yiyang Zhao
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, College of Biological Sciences and Technology, Beijing Forestry University, Ministry of Education, Beijing, China
| | - Sisi Chen
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, College of Biological Sciences and Technology, Beijing Forestry University, Ministry of Education, Beijing, China
| | - Jianbo Xie
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, College of Biological Sciences and Technology, Beijing Forestry University, Ministry of Education, Beijing, China
| | - Deqiang Zhang
- Beijing Advanced Innovation Center for Tree Breeding by Molecular Design, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, Beijing, China
- Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, College of Biological Sciences and Technology, Beijing Forestry University, Ministry of Education, Beijing, China
| |
Collapse
|
20
|
David KT, Oaks JR, Halanych KM. Patterns of gene evolution following duplications and speciations in vertebrates. PeerJ 2020; 8:e8813. [PMID: 32266119 PMCID: PMC7120047 DOI: 10.7717/peerj.8813] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2019] [Accepted: 02/27/2020] [Indexed: 11/24/2022] Open
Abstract
BACKGROUND Eukaryotic genes typically form independent evolutionary lineages through either speciation or gene duplication events. Generally, gene copies resulting from speciation events (orthologs) are expected to maintain similarity over time with regard to sequence, structure and function. After a duplication event, however, resulting gene copies (paralogs) may experience a broader set of possible fates, including partial (subfunctionalization) or complete loss of function, as well as gain of new function (neofunctionalization). This assumption, known as the Ortholog Conjecture, is prevalent throughout molecular biology and notably plays an important role in many functional annotation methods. Unfortunately, studies that explicitly compare evolutionary processes between speciation and duplication events are rare and conflicting. METHODS To provide an empirical assessment of ortholog/paralog evolution, we estimated ratios of nonsynonymous to synonymous substitutions (ω = dN/dS) for 251,044 lineages in 6,244 gene trees across 77 vertebrate taxa. RESULTS Overall, we found ω to be more similar between lineages descended from speciation events (p < 0.001) than lineages descended from duplication events, providing strong support for the Ortholog Conjecture. The asymmetry in ω following duplication events appears to be largely driven by an increase along one of the paralogous lineages, while the other remains similar to the parent. This trend is commonly associated with neofunctionalization, suggesting that gene duplication is a significant mechanism for generating novel gene functions.
Collapse
Affiliation(s)
- Kyle T. David
- Department of Biological Sciences, Auburn University, Auburn, AL, USA
| | - Jamie R. Oaks
- Department of Biological Sciences, Auburn University, Auburn, AL, USA
| | | |
Collapse
|
21
|
Youlten SE, Baldock PA. Using mouse genetics to understand human skeletal disease. Bone 2019; 126:27-36. [PMID: 30776501 DOI: 10.1016/j.bone.2019.02.015] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/03/2018] [Revised: 01/25/2019] [Accepted: 02/13/2019] [Indexed: 02/06/2023]
Abstract
Technological advances have enabled the study of the human genome in incredible detail with relative ease. However, our ability to interpret the functional significance of the millions of genetic variants present within each individual is limited. As a result, the confident assignment of disease-causing variant calls remains a significant challenge. Here we explore how mouse genetics can help address this deficit in functional genomic understanding. Underpinned by marked genetic correspondence, skeletal biology shows inter-species similarities which provide important opportunities to use data from mouse models to direct research into the genetic basis of skeletal pathophysiology. In this article we outline critical resources that may be used to establish genotype/phenotype relationships in skeletal tissue, identify genes with established skeletal effects and define the transcriptome of critical skeletal cell types. Finally, we outline how these mouse resources might be utilized to progress from a list of human sequence variants toward plausible gene candidates that contribute to skeletal disease.
Collapse
Affiliation(s)
- Scott E Youlten
- Division of Bone Biology, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia; St Vincent's Clinical School, Faculty of Medicine, UNSW Australia, Sydney, NSW, 2010, Australia.
| | - Paul A Baldock
- Division of Bone Biology, Garvan Institute of Medical Research, Sydney, NSW 2010, Australia; St Vincent's Clinical School, Faculty of Medicine, UNSW Australia, Sydney, NSW, 2010, Australia; University of Notre Dame Australia, Sydney, NSW, 2010, Australia
| |
Collapse
|
22
|
Cutter AD, Garrett RH, Mark S, Wang W, Sun L. Molecular evolution across developmental time reveals rapid divergence in early embryogenesis. Evol Lett 2019; 3:359-373. [PMID: 31388446 PMCID: PMC6675142 DOI: 10.1002/evl3.122] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2019] [Accepted: 05/30/2019] [Indexed: 12/16/2022] Open
Abstract
Ontogenetic development hinges on the changes in gene expression in time and space within an organism, suggesting that the demands of ontogenetic growth can impose or reveal predictable pattern in the molecular evolution of genes expressed dynamically across development. Here, we characterize coexpression modules of the Caenorhabditis elegans transcriptome, using a time series of 30 points from early embryo to adult. By capturing the functional form of expression profiles with quantitative metrics, we find fastest evolution in the distinctive set of genes with transcript abundance that declines through development from a peak in young embryos. These genes are highly enriched for oogenic function and transient early zygotic expression, are nonrandomly distributed in the genome, and correspond to a life stage especially prone to inviability in interspecies hybrids. These observations conflict with the "early conservation model" for the evolution of development, although expression-weighted sequence divergence analysis provides some support for the "hourglass model." Genes in coexpression modules that peak toward adulthood also evolve fast, being hyper-enriched for roles in spermatogenesis, implicating a history of sexual selection and relaxation of selection on sperm as key factors driving rapid change to ontogenetically distinguishable coexpression modules of genes. We propose that these predictable trends of molecular evolution for dynamically expressed genes across ontogeny predispose particular life stages, early embryogenesis in particular, to hybrid dysfunction in the speciation process.
Collapse
Affiliation(s)
- Asher D. Cutter
- Department of Ecology and Evolutionary BiologyUniversity of TorontoTorontoONM6G1W3Canada
| | - Rose H. Garrett
- Department of Ecology and Evolutionary BiologyUniversity of TorontoTorontoONM6G1W3Canada
- Division of Biostatistics, Dalla Lana School of Public HealthUniversity of TorontoTorontoONM6G1W3Canada
- Department of Statistical SciencesUniversity of TorontoTorontoONM6G1W3Canada
| | - Stephanie Mark
- Department of Ecology and Evolutionary BiologyUniversity of TorontoTorontoONM6G1W3Canada
| | - Wei Wang
- Department of Ecology and Evolutionary BiologyUniversity of TorontoTorontoONM6G1W3Canada
| | - Lei Sun
- Division of Biostatistics, Dalla Lana School of Public HealthUniversity of TorontoTorontoONM6G1W3Canada
- Department of Statistical SciencesUniversity of TorontoTorontoONM6G1W3Canada
| |
Collapse
|
23
|
Mark S, Weiss J, Sharma E, Liu T, Wang W, Claycomb JM, Cutter AD. Genome structure predicts modular transcriptome responses to genetic and environmental conditions. Mol Ecol 2019; 28:3681-3697. [PMID: 31325381 DOI: 10.1111/mec.15185] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2019] [Revised: 07/05/2019] [Accepted: 07/10/2019] [Indexed: 12/13/2022]
Abstract
Understanding the plasticity, robustness and modularity of transcriptome expression to genetic and environmental conditions is crucial to deciphering how organisms adapt in nature. To test how genome architecture influences transcriptome profiles, we quantified expression responses for distinct temperature-adapted genotypes of the nematode Caenorhabditis briggsae when exposed to chronic temperature stresses throughout development. We found that 56% of the 8,795 differentially expressed genes show genotype-specific changes in expression in response to temperature (genotype-by-environment interactions, GxE). Most genotype-specific responses occur under heat stress, indicating that cold vs. heat stress responses involve distinct genomic architectures. The 22 co-expression modules that we identified differ in their enrichment of genes with genetic vs. environmental vs. interaction effects, as well as their genomic spatial distributions, functional attributes and rates of molecular evolution at the sequence level. Genes in modules enriched for simple effects of either genotype or temperature alone tend to evolve especially rapidly, consistent with disproportionate influence of adaptation or weaker constraint on these subsets of loci. Chromosome-scale heterogeneity in nucleotide polymorphism, however, rather than the scale of individual genes predominates as the source of genetic differences among expression profiles, and natural selection regimes are largely decoupled between coding sequences and noncoding flanking sequences that contain cis-regulatory elements. These results illustrate how the form of transcriptome modularity and genome structure contribute to predictable profiles of evolutionary change.
Collapse
Affiliation(s)
- Stephanie Mark
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, ON, Canada
| | - Joerg Weiss
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, ON, Canada
| | - Eesha Sharma
- Department of Molecular Genetics, University of Toronto, Toronto, ON, Canada
| | - Ting Liu
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, ON, Canada
| | - Wei Wang
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, ON, Canada
| | - Julie M Claycomb
- Department of Molecular Genetics, University of Toronto, Toronto, ON, Canada
| | - Asher D Cutter
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, ON, Canada
| |
Collapse
|
24
|
Wang Y, Gao S, Zhao Y, Chen WH, Shao JJ, Wang NN, Li M, Zhou GX, Wang L, Shen WJ, Xu JT, Deng WD, Wang W, Chen YL, Jiang Y. Allele-specific expression and alternative splicing in horse×donkey and cattle×yak hybrids. Zool Res 2019; 40:293-304. [PMID: 31271004 PMCID: PMC6680129 DOI: 10.24272/j.issn.2095-8137.2019.042] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
Divergence of gene expression and alternative splicing is a crucial driving force in the evolution of species; to date, however the molecular mechanism remains unclear. Hybrids of closely related species provide a suitable model to analyze allele-specific expression (ASE) and allele-specific alternative splicing (ASS). Analysis of ASE and ASS can uncover the differences in cis-regulatory elements between closely related species, while eliminating interference of trans-regulatory elements. Here, we provide a detailed characterization of ASE and ASS from 19 and 10 transcriptome datasets across five tissues from reciprocal-cross hybrids of horse×donkey (mule/hinny) and cattle×yak (dzo), respectively. Results showed that 4.8%-8.7% and 10.8%-16.7% of genes exhibited ASE and ASS, respectively. Notably, lncRNAs and pseudogenes were more likely to show ASE than protein-coding genes. In addition, genes showing ASE and ASS in mule/hinny were found to be involved in the regulation of muscle strength, whereas those of dzo were involved in high-altitude adaptation. In conclusion, our study demonstrated that exploration of genes showing ASE and ASS in hybrids of closely related species is feasible for species evolution research.
Collapse
Affiliation(s)
- Yu Wang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling Shaanxi 712100, China
| | - Shan Gao
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling Shaanxi 712100, China
| | - Yue Zhao
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling Shaanxi 712100, China
| | - Wei-Huang Chen
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling Shaanxi 712100, China
| | - Jun-Jie Shao
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling Shaanxi 712100, China
| | - Ni-Ni Wang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling Shaanxi 712100, China
| | - Ming Li
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling Shaanxi 712100, China
| | - Guang-Xian Zhou
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling Shaanxi 712100, China
| | - Lei Wang
- Stake Key Laboratory of Plateau Ecology and Agriculture, Qinghai Academy of Animal Science and Veterinary Medicine, Qinghai University, Xining Qinghai 810016, China
| | - Wen-Jing Shen
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming Yunnan 650223, China
| | - Jing-Tao Xu
- Stake Key Laboratory of Plateau Ecology and Agriculture, Qinghai Academy of Animal Science and Veterinary Medicine, Qinghai University, Xining Qinghai 810016, China
| | - Wei-Dong Deng
- Faculty of Animal Science and Technology, Yunnan Agricultural University, Kunming Yunnan 650223, China
| | - Wen Wang
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming Yunnan 650223, China
| | - Yu-Lin Chen
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling Shaanxi 712100, China
| | - Yu Jiang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling Shaanxi 712100, China
| |
Collapse
|
25
|
Chen W, Wu X, Yan X, Xu A, Yang A, You H. Multitranscriptome analyses reveal prioritized genes specifically associated with liver fibrosis progression independent of etiology. Am J Physiol Gastrointest Liver Physiol 2019; 316:G744-G754. [PMID: 30920297 DOI: 10.1152/ajpgi.00339.2018] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Elimination or suppression of causative factors can raise the possibility of liver fibrosis regression. However, different injurious stimuli will give fibrosis from somewhat different etiologies, which, in turn, may hamper the discovery of liver fibrosis-specific therapeutic drugs. Therefore, the analogical cellular and molecular events shared by various etiology-evoked liver fibrosis should be clarified. Our present study systematically integrated five publicly available transcriptomic data sets regarding liver fibrosis with different etiologies from the Gene Expression Omnibus database and performed a series of bioinformatics analyses and experimental verifications. A total of 111 significantly upregulated and 16 downregulated genes were identified specific to liver fibrosis independent of any etiology. These genes were predominately enriched in some Kyoto Encyclopedia of Genes and Genomes pathways, including the "PI3K-AKT signaling pathway," "Focal adhesion," and "ECM-receptor interaction." Subsequently, five prioritized liver fibrosis-specific genes, including COL4A2, THBS2, ITGAV, LAMB1, and PDGFRA, were screened. These genes were positively associated with each other and liver fibrosis progression. In addition, they could robustly separate all stages of samples in both training and validation data sets with diverse etiologies when they were regarded as observed variables applied to principal component analysis plots. Expressions of all five genes were confirmed in activated primary mouse hepatic stellate cells (HSCs) and transforming growth factor β1-treated LX-2 cells. Moreover, THBS2 protein was enhanced in liver fibrosis rodent models, which could promote HSC activation and proliferation and facilitate NOTCH1/JAG1 expression in HSCs. Overall, our current study may provide potential targets for liver fibrosis therapy and aid to a deeper understanding of the molecular underpinnings of liver fibrosis. NEW & NOTEWORTHY Prioritized liver fibrosis-specific genes THBS2, COL4A2, ITGAV, LAMB1, and PDGFRA were identified and significantly associated with liver fibrosis progression and could be combined to discriminate liver fibrosis stages regardless of any etiology. Among the identified prioritized liver fibrosis-specific targets, THBS2 protein was confirmed to be enhanced in liver fibrosis rodent models, which could promote hepatic stellate cell (HSC) activation and proliferation and facilitate NOTCH1/JAG1 expression in HSCs.
Collapse
Affiliation(s)
- Wei Chen
- Experimental and Translational Research Center, Beijing Friendship Hospital, Capital Medical University , Beijing , China.,Beijing Key Laboratory of Tolerance Induction and Organ Protection in Transplantation, Beijing Friendship Hospital, Capital Medical University , Beijing , China
| | - Xiaoning Wu
- Liver Research Center, Beijing Friendship Hospital, Capital Medical University, Beijing Key Laboratory of Translational Medicine in Liver Cirrhosis, National Clinical Research Center of Digestive Diseases , Beijing , China
| | - Xuzhen Yan
- Liver Research Center, Beijing Friendship Hospital, Capital Medical University, Beijing Key Laboratory of Translational Medicine in Liver Cirrhosis, National Clinical Research Center of Digestive Diseases , Beijing , China
| | - Anjian Xu
- Experimental and Translational Research Center, Beijing Friendship Hospital, Capital Medical University , Beijing , China.,Beijing Key Laboratory of Tolerance Induction and Organ Protection in Transplantation, Beijing Friendship Hospital, Capital Medical University , Beijing , China
| | - Aiting Yang
- Experimental and Translational Research Center, Beijing Friendship Hospital, Capital Medical University , Beijing , China.,Beijing Key Laboratory of Tolerance Induction and Organ Protection in Transplantation, Beijing Friendship Hospital, Capital Medical University , Beijing , China
| | - Hong You
- Experimental and Translational Research Center, Beijing Friendship Hospital, Capital Medical University , Beijing , China.,Beijing Key Laboratory of Tolerance Induction and Organ Protection in Transplantation, Beijing Friendship Hospital, Capital Medical University , Beijing , China.,Liver Research Center, Beijing Friendship Hospital, Capital Medical University, Beijing Key Laboratory of Translational Medicine in Liver Cirrhosis, National Clinical Research Center of Digestive Diseases , Beijing , China
| |
Collapse
|
26
|
Assis R. Lineage-Specific Expression Divergence in Grasses Is Associated with Male Reproduction, Host-Pathogen Defense, and Domestication. Genome Biol Evol 2019; 11:207-219. [PMID: 30398650 PMCID: PMC6331041 DOI: 10.1093/gbe/evy245] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/03/2018] [Indexed: 02/02/2023] Open
Abstract
Poaceae (grasses) is an agriculturally important and widely distributed family of plants with extraordinary phenotypic diversity, much of which was generated under recent lineage-specific evolution. Yet, little is known about the genes and functional modules involved in the lineage-specific divergence of grasses. Here, I address this question on a genome-wide scale by applying a novel branch-based statistic of lineage-specific expression divergence, LED, to RNA-seq data from nine tissues of the wild grass Brachypodium distachyon and its domesticated relatives Oryza sativa japonica (rice) and Sorghum bicolor (sorghum). I find that LED is generally smallest in B. distachyon and largest in O. sativa japonica, which underwent domestication earlier than S. bicolor, supporting the hypothesis that domestication may increase the rate of lineage-specific expression divergence in grasses. Moreover, in all three species, LED is positively correlated with protein-coding sequence divergence and tissue specificity, and negatively correlated with network connectivity. Further analysis reveals that genes with large LED are often primarily expressed in anther, implicating lineage-specific expression divergence in the evolution of male reproductive phenotypes. Gene ontology enrichment analysis also identifies an overrepresentation of terms related to male reproduction in the two domesticated grasses, as well as to those involved in host-pathogen defense in all three species. Last, examinations of genes with the largest LED reveal that their lineage-specific expression divergence may have contributed to antimicrobial functions in B. distachyon, to enhanced adaptation and yield during domestication in O. sativa japonica, and to defense against a widespread and devastating fungal pathogen in S. bicolor. Together, these findings suggest that lineage-specific expression divergence in grasses may increase under domestication and preferentially target rapidly evolving genes involved in male reproduction, host-pathogen defense, and the origin of domesticated phenotypes.
Collapse
Affiliation(s)
- Raquel Assis
- Department of Biology, Pennsylvania State University, University Park
| |
Collapse
|
27
|
Ghiselli F, Iannello M, Puccio G, Chang PL, Plazzi F, Nuzhdin SV, Passamonti M. Comparative Transcriptomics in Two Bivalve Species Offers Different Perspectives on the Evolution of Sex-Biased Genes. Genome Biol Evol 2018; 10:1389-1402. [PMID: 29897459 PMCID: PMC6007409 DOI: 10.1093/gbe/evy082] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 04/19/2018] [Indexed: 12/13/2022] Open
Abstract
Comparative genomics has become a central tool for evolutionary biology, and a better knowledge of understudied taxa represents the foundation for future work. In this study, we characterized the transcriptome of male and female mature gonads in the European clam Ruditapes decussatus, compared with that in the Manila clam Ruditapes philippinarum providing, for the first time in bivalves, information about transcription dynamics and sequence evolution of sex-biased genes. In both the species, we found a relatively low number of sex-biased genes (1,284, corresponding to 41.3% of the orthologous genes between the two species), probably due to the absence of sexual dimorphism, and the transcriptional bias is maintained in only 33% of the orthologs. The dN/dS is generally low, indicating purifying selection, with genes where the female-biased transcription is maintained between the two species showing a significantly higher dN/dS. Genes involved in embryo development, cell proliferation, and maintenance of genome stability show a faster sequence evolution. Finally, we report a lack of clear correlation between transcription level and evolutionary rate in these species, in contrast with studies that reported a negative correlation. We discuss such discrepancy and call into question some methodological approaches and rationales generally used in this type of comparative studies.
Collapse
Affiliation(s)
- Fabrizio Ghiselli
- Department of Biological, Geological, and Environmental Sciences, University of Bologna, Italy
| | - Mariangela Iannello
- Department of Biological, Geological, and Environmental Sciences, University of Bologna, Italy
| | - Guglielmo Puccio
- Department of Biological, Geological, and Environmental Sciences, University of Bologna, Italy
| | - Peter L Chang
- Program in Molecular and Computational Biology, Department of Biological Sciences, University of Southern California, Los Angeles, USA
| | - Federico Plazzi
- Department of Biological, Geological, and Environmental Sciences, University of Bologna, Italy
| | - Sergey V Nuzhdin
- Program in Molecular and Computational Biology, Department of Biological Sciences, University of Southern California, Los Angeles, USA
| | - Marco Passamonti
- Department of Biological, Geological, and Environmental Sciences, University of Bologna, Italy
| |
Collapse
|
28
|
Conserved noncoding sequences conserve biological networks and influence genome evolution. Heredity (Edinb) 2018; 120:437-451. [PMID: 29396421 DOI: 10.1038/s41437-018-0055-4] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2017] [Revised: 12/14/2017] [Accepted: 01/08/2018] [Indexed: 01/24/2023] Open
Abstract
Comparative genomics approaches have identified numerous conserved cis-regulatory sequences near genes in plant genomes. Despite the identification of these conserved noncoding sequences (CNSs), our knowledge of their functional importance and selection remains limited. Here, we used a combination of DNA methylome analysis, microarray expression analyses, and functional annotation to study these sequences in the model tree Populus trichocarpa. Methylation in CG contexts and non-CG contexts was lower in CNSs, particularly CNSs in the 5'-upstream regions of genes, compared with other sites in the genome. We observed that CNSs are enriched in genes with transcription and binding functions, and this also associated with syntenic genes and those from whole-genome duplications, suggesting that cis-regulatory sequences play a key role in genome evolution. We detected a significant positive correlation between CNS number and protein interactions, suggesting that CNSs may have roles in the evolution and maintenance of biological networks. The divergence of CNSs indicates that duplication-degeneration-complementation drives the subfunctionalization of a proportion of duplicated genes from whole-genome duplication. Furthermore, population genomics confirmed that most CNSs are under strong purifying selection and only a small subset of CNSs shows evidence of adaptive evolution. These findings provide a foundation for future studies exploring these key genomic features in the maintenance of biological networks, local adaptation, and transcription.
Collapse
|
29
|
Celedon JM, Yuen MMS, Chiang A, Henderson H, Reid KE, Bohlmann J. Cell-type- and tissue-specific transcriptomes of the white spruce (Picea glauca) bark unmask fine-scale spatial patterns of constitutive and induced conifer defense. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2017; 92:710-726. [PMID: 28857307 DOI: 10.1111/tpj.13673] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/24/2017] [Revised: 08/07/2017] [Accepted: 08/22/2017] [Indexed: 05/09/2023]
Abstract
Plant defenses often involve specialized cells and tissues. In conifers, specialized cells of the bark are important for defense against insects and pathogens. Using laser microdissection, we characterized the transcriptomes of cortical resin duct cells, phenolic cells and phloem of white spruce (Picea glauca) bark under constitutive and methyl jasmonate (MeJa)-induced conditions, and we compared these transcriptomes with the transcriptome of the bark tissue complex. Overall, ~3700 bark transcripts were differentially expressed in response to MeJa. Approximately 25% of transcripts were expressed in only one cell type, revealing cell specialization at the transcriptome level. MeJa caused cell-type-specific transcriptome responses and changed the overall patterns of cell-type-specific transcript accumulation. Comparison of transcriptomes of the conifer bark tissue complex and specialized cells resolved a masking effect inherent to transcriptome analysis of complex tissues, and showed the actual cell-type-specific transcriptome signatures. Characterization of cell-type-specific transcriptomes is critical to reveal the dynamic patterns of spatial and temporal display of constitutive and induced defense systems in a complex plant tissue or organ. This was demonstrated with the improved resolution of spatially restricted expression of sets of genes of secondary metabolism in the specialized cell types.
Collapse
Affiliation(s)
- Jose M Celedon
- Michael Smith Laboratories, University of British Columbia, 2185 East Mall, Vancouver, BC, V6T 1Z4, Canada
| | - Macaire M S Yuen
- Michael Smith Laboratories, University of British Columbia, 2185 East Mall, Vancouver, BC, V6T 1Z4, Canada
| | - Angela Chiang
- Michael Smith Laboratories, University of British Columbia, 2185 East Mall, Vancouver, BC, V6T 1Z4, Canada
| | - Hannah Henderson
- Michael Smith Laboratories, University of British Columbia, 2185 East Mall, Vancouver, BC, V6T 1Z4, Canada
| | - Karen E Reid
- Michael Smith Laboratories, University of British Columbia, 2185 East Mall, Vancouver, BC, V6T 1Z4, Canada
| | - Jörg Bohlmann
- Michael Smith Laboratories, University of British Columbia, 2185 East Mall, Vancouver, BC, V6T 1Z4, Canada
| |
Collapse
|
30
|
Begum T, Ghosh TC, Basak S. Systematic Analyses and Prediction of Human Drug Side Effect Associated Proteins from the Perspective of Protein Evolution. Genome Biol Evol 2017; 9:337-350. [PMID: 28391292 PMCID: PMC5499873 DOI: 10.1093/gbe/evw301] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/16/2017] [Indexed: 12/20/2022] Open
Abstract
Identification of various factors involved in adverse drug reactions in target proteins to develop therapeutic drugs with minimal/no side effect is very important. In this context, we have performed a comparative evolutionary rate analyses between the genes exhibiting drug side-effect(s) (SET) and genes showing no side effect (NSET) with an aim to increase the prediction accuracy of SET/NSET proteins using evolutionary rate determinants. We found that SET proteins are more conserved than the NSET proteins. The rates of evolution between SET and NSET protein primarily depend upon their noncomplex (protein complex association number = 0) forming nature, phylogenetic age, multifunctionality, membrane localization, and transmembrane helix content irrespective of their essentiality, total druggability (total number of drugs/target), m-RNA expression level, and tissue expression breadth. We also introduced two novel terms—killer druggability (number of drugs with killing side effect(s)/target), essential druggability (number of drugs targeting essential proteins/target) to explain the evolutionary rate variation between SET and NSET proteins. Interestingly, we noticed that SET proteins are younger than NSET proteins and multifunctional younger SET proteins are candidates of acquiring killing side effects. We provide evidence that higher killer druggability, multifunctionality, and transmembrane helices support the conservation of SET proteins over NSET proteins in spite of their recent origin. By employing all these entities, our Support Vector Machine model predicts human SET/NSET proteins to a high degree of accuracy (∼86%).
Collapse
Affiliation(s)
- Tina Begum
- Bioinformatics Centre, Tripura University, Suryamaninagar, Tripura, India
| | | | - Surajit Basak
- Bioinformatics Centre, Tripura University, Suryamaninagar, Tripura, India.,Department of Molecular Biology & Bioinformatics, Tripura University, Suryamaninagar, Tripura, India
| |
Collapse
|
31
|
Chen W, Zhao W, Yang A, Xu A, Wang H, Cong M, Liu T, Wang P, You H. Integrated analysis of microRNA and gene expression profiles reveals a functional regulatory module associated with liver fibrosis. Gene 2017; 636:87-95. [PMID: 28919164 DOI: 10.1016/j.gene.2017.09.027] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2017] [Revised: 07/13/2017] [Accepted: 09/13/2017] [Indexed: 02/06/2023]
Abstract
BACKGROUND Liver fibrosis, characterized with the excessive accumulation of extracellular matrix (ECM) proteins, represents the final common pathway of chronic liver inflammation. Ever-increasing evidence indicates microRNAs (miRNAs) dysregulation has important implications in the different stages of liver fibrosis. However, our knowledge of miRNA-gene regulation details pertaining to such disease remains unclear. METHODS The publicly available Gene Expression Omnibus (GEO) datasets of patients suffered from cirrhosis were extracted for integrated analysis. Differentially expressed miRNAs (DEMs) and genes (DEGs) were identified using GEO2R web tool. Putative target gene prediction of DEMs was carried out using the intersection of five major algorithms: DIANA-microT, TargetScan, miRanda, PICTAR5 and miRWalk. Functional miRNA-gene regulatory network (FMGRN) was constructed based on the computational target predictions at the sequence level and the inverse expression relationships between DEMs and DEGs. DAVID web server was selected to perform KEGG pathway enrichment analysis. Functional miRNA-gene regulatory module was generated based on the biological interpretation. Internal connections among genes in liver fibrosis-related module were determined using String database. MiRNA-gene regulatory modules related to liver fibrosis were experimentally verified in recombinant human TGFβ1 stimulated and specific miRNA inhibitor treated LX-2 cells. RESULTS We totally identified 85 and 923 dysregulated miRNAs and genes in liver cirrhosis biopsy samples compared to their normal controls. All evident miRNA-gene pairs were identified and assembled into FMGRN which consisted of 990 regulations between 51 miRNAs and 275 genes, forming two big sub-networks that were defined as down-network and up-network, respectively. KEGG pathway enrichment analysis revealed that up-network was prominently involved in several KEGG pathways, in which "Focal adhesion", "PI3K-Akt signaling pathway" and "ECM-receptor interaction" were remarked significant (adjusted p<0.001). Genes enriched in these pathways coupled with their regulatory miRNAs formed a functional miRNA-gene regulatory module that contains 7 miRNAs, 22 genes and 42 miRNA-gene connections. Gene interaction analysis based on String database revealed that 8 out of 22 genes were highly clustered. Finally, we experimentally confirmed a functional regulatory module containing 5 miRNAs (miR-130b-3p, miR-148a-3p, miR-345-5p, miR-378a-3p, and miR-422a) and 6 genes (COL6A1, COL6A2, COL6A3, PIK3R3, COL1A1, CCND2) associated with liver fibrosis. CONCLUSIONS Our integrated analysis of miRNA and gene expression profiles highlighted a functional miRNA-gene regulatory module associated with liver fibrosis, which, to some extent, may provide important clues to better understand the underlying pathogenesis of liver fibrosis.
Collapse
Affiliation(s)
- Wei Chen
- Experimental and Translational Research Center, Beijing Friendship Hospital, Capital Medical University, Beijing, China; Beijing Key Laboratory of Tolerance Induction and Organ Protection in Transplantation, Beijing Friendship Hospital, Capital Medical University, Beijing, China
| | - Wenshan Zhao
- Liver Research Center, Beijing Friendship Hospital, Capital Medical University, Beijing Key Laboratory of Translational Medicine in Liver Cirrhosis, National Clinical Research Center of Digestive Diseases, Beijing, China
| | - Aiting Yang
- Experimental and Translational Research Center, Beijing Friendship Hospital, Capital Medical University, Beijing, China; Beijing Key Laboratory of Tolerance Induction and Organ Protection in Transplantation, Beijing Friendship Hospital, Capital Medical University, Beijing, China
| | - Anjian Xu
- Experimental and Translational Research Center, Beijing Friendship Hospital, Capital Medical University, Beijing, China; Beijing Key Laboratory of Tolerance Induction and Organ Protection in Transplantation, Beijing Friendship Hospital, Capital Medical University, Beijing, China
| | - Huan Wang
- Experimental and Translational Research Center, Beijing Friendship Hospital, Capital Medical University, Beijing, China; Beijing Key Laboratory of Tolerance Induction and Organ Protection in Transplantation, Beijing Friendship Hospital, Capital Medical University, Beijing, China
| | - Min Cong
- Liver Research Center, Beijing Friendship Hospital, Capital Medical University, Beijing Key Laboratory of Translational Medicine in Liver Cirrhosis, National Clinical Research Center of Digestive Diseases, Beijing, China
| | - Tianhui Liu
- Liver Research Center, Beijing Friendship Hospital, Capital Medical University, Beijing Key Laboratory of Translational Medicine in Liver Cirrhosis, National Clinical Research Center of Digestive Diseases, Beijing, China
| | - Ping Wang
- Liver Research Center, Beijing Friendship Hospital, Capital Medical University, Beijing Key Laboratory of Translational Medicine in Liver Cirrhosis, National Clinical Research Center of Digestive Diseases, Beijing, China
| | - Hong You
- Experimental and Translational Research Center, Beijing Friendship Hospital, Capital Medical University, Beijing, China; Beijing Key Laboratory of Tolerance Induction and Organ Protection in Transplantation, Beijing Friendship Hospital, Capital Medical University, Beijing, China; Liver Research Center, Beijing Friendship Hospital, Capital Medical University, Beijing Key Laboratory of Translational Medicine in Liver Cirrhosis, National Clinical Research Center of Digestive Diseases, Beijing, China.
| |
Collapse
|
32
|
Yang JR, Maclean CJ, Park C, Zhao H, Zhang J. Intra and Interspecific Variations of Gene Expression Levels in Yeast Are Largely Neutral: (Nei Lecture, SMBE 2016, Gold Coast). Mol Biol Evol 2017; 34:2125-2139. [PMID: 28575451 PMCID: PMC5850415 DOI: 10.1093/molbev/msx171] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
It is commonly, although not universally, accepted that most intra and interspecific genome sequence variations are more or less neutral, whereas a large fraction of organism-level phenotypic variations are adaptive. Gene expression levels are molecular phenotypes that bridge the gap between genotypes and corresponding organism-level phenotypes. Yet, it is unknown whether natural variations in gene expression levels are mostly neutral or adaptive. Here we address this fundamental question by genome-wide profiling and comparison of gene expression levels in nine yeast strains belonging to three closely related Saccharomyces species and originating from five different ecological environments. We find that the transcriptome-based clustering of the nine strains approximates the genome sequence-based phylogeny irrespective of their ecological environments. Remarkably, only ∼0.5% of genes exhibit similar expression levels among strains from a common ecological environment, no greater than that among strains with comparable phylogenetic relationships but different environments. These and other observations strongly suggest that most intra and interspecific variations in yeast gene expression levels result from the accumulation of random mutations rather than environmental adaptations. This finding has profound implications for understanding the driving force of gene expression evolution, genetic basis of phenotypic adaptation, and general role of stochasticity in evolution.
Collapse
Affiliation(s)
- Jian-Rong Yang
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI
| | - Calum J. Maclean
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI
| | - Chungoo Park
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI
| | - Huabin Zhao
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI
| | - Jianzhi Zhang
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI
| |
Collapse
|
33
|
Abstract
Cross-species comparisons of genomes, transcriptomes and gene regulation are now feasible at unprecedented resolution and throughput, enabling the comparison of human and mouse biology at the molecular level. Insights have been gained into the degree of conservation between human and mouse at the level of not only gene expression but also epigenetics and inter-individual variation. However, a number of limitations exist, including incomplete transcriptome characterization and difficulties in identifying orthologous phenotypes and cell types, which are beginning to be addressed by emerging technologies. Ultimately, these comparisons will help to identify the conditions under which the mouse is a suitable model of human physiology and disease, and optimize the use of animal models.
Collapse
Affiliation(s)
- Alessandra Breschi
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), 08002 Barcelona, Spain
| | - Thomas R Gingeras
- Cold Spring Harbor Laboratory, Cold Spring Harbor, New York 11742, USA
| | - Roderic Guigó
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), 08002 Barcelona, Spain
| |
Collapse
|
34
|
Du Y, Deng W, Wang Z, Ning M, Zhang W, Zhou Y, Lo EH, Xing C. Differential subnetwork of chemokines/cytokines in human, mouse, and rat brain cells after oxygen-glucose deprivation. J Cereb Blood Flow Metab 2017; 37:1425-1434. [PMID: 27328691 PMCID: PMC5453462 DOI: 10.1177/0271678x16656199] [Citation(s) in RCA: 47] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Mice and rats are the most commonly used animals for preclinical stroke studies, but it is unclear whether targets and mechanisms are always the same across different species. Here, we mapped the baseline expression of a chemokine/cytokine subnetwork and compared responses after oxygen-glucose deprivation in primary neurons, astrocytes, and microglia from mouse, rat, and human. Baseline profiles of chemokines (CX3CL1, CXCL12, CCL2, CCL3, and CXCL10) and cytokines (IL-1α, IL-1β, IL-6, IL-10, and TNFα) showed significant differences between human and rodents. The response of chemokines/cytokines to oxygen-glucose deprivation was also significantly different between species. After 4 h oxygen-glucose deprivation and 4 h reoxygenation, human and rat neurons showed similar changes with a downregulation in many chemokines, whereas mouse neurons showed a mixed response with up- and down-regulated genes. For astrocytes, subnetwork response patterns were more similar in rats and mice compared to humans. For microglia, rat cells showed an upregulation in all chemokines/cytokines, mouse cells had many down-regulated genes, and human cells showed a mixed response with up- and down-regulated genes. This study provides proof-of-concept that species differences exist in chemokine/cytokine subnetworks in brain cells that may be relevant to stroke pathophysiology. Further investigation of differential gene pathways across species is warranted.
Collapse
Affiliation(s)
- Yang Du
- 1 Department of Neurology, Beijing Tiantan Hospital, Capital Medical University, Beijing, China.,2 Departments of Radiology and Neurology, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, USA.,3 Department of Geriatrics, Beijing Tiantan Hospital, Capital Medical University, Beijing, China
| | - Wenjun Deng
- 4 Department of Neurology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | - Zixing Wang
- 5 Departments of Bioinformatics and Computational Biology, The University of Texas MD Anderson Cancer Center, Houston, TX, USA
| | - MingMing Ning
- 4 Department of Neurology, Massachusetts General Hospital, Harvard Medical School, Boston, MA, USA
| | - Wei Zhang
- 1 Department of Neurology, Beijing Tiantan Hospital, Capital Medical University, Beijing, China.,3 Department of Geriatrics, Beijing Tiantan Hospital, Capital Medical University, Beijing, China
| | - Yiming Zhou
- 2 Departments of Radiology and Neurology, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, USA
| | - Eng H Lo
- 2 Departments of Radiology and Neurology, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, USA
| | - Changhong Xing
- 2 Departments of Radiology and Neurology, Massachusetts General Hospital, Harvard Medical School, Charlestown, MA, USA
| |
Collapse
|
35
|
The Evolution and Expression Pattern of Human Overlapping lncRNA and Protein-coding Gene Pairs. Sci Rep 2017; 7:42775. [PMID: 28344339 PMCID: PMC5366806 DOI: 10.1038/srep42775] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2016] [Accepted: 01/13/2017] [Indexed: 12/27/2022] Open
Abstract
Long non-coding RNA overlapping with protein-coding gene (lncRNA-coding pair) is a special type of overlapping genes. Protein-coding overlapping genes have been well studied and increasing attention has been paid to lncRNAs. By studying lncRNA-coding pairs in human genome, we showed that lncRNA-coding pairs were more likely to be generated by overprinting and retaining genes in lncRNA-coding pairs were given higher priority than non-overlapping genes. Besides, the preference of overlapping configurations preserved during evolution was based on the origin of lncRNA-coding pairs. Further investigations showed that lncRNAs promoting the splicing of their embedded protein-coding partners was a unilateral interaction, but the existence of overlapping partners improving the gene expression was bidirectional and the effect was decreased with the increased evolutionary age of genes. Additionally, the expression of lncRNA-coding pairs showed an overall positive correlation and the expression correlation was associated with their overlapping configurations, local genomic environment and evolutionary age of genes. Comparison of the expression correlation of lncRNA-coding pairs between normal and cancer samples found that the lineage-specific pairs including old protein-coding genes may play an important role in tumorigenesis. This work presents a systematically comprehensive understanding of the evolution and the expression pattern of human lncRNA-coding pairs.
Collapse
|
36
|
Wada N, Kazuki Y, Kazuki K, Inoue T, Fukui K, Oshimura M. Maintenance and Function of a Plant Chromosome in Human Cells. ACS Synth Biol 2017; 6:301-310. [PMID: 27696824 DOI: 10.1021/acssynbio.6b00180] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Replication, segregation, gene expression, and inheritance are essential features of all eukaryotic chromosomes. To delineate the extent of conservation of chromosome functions between humans and plants during evolutionary history, we have generated the first human cell line containing an Arabidopsis chromosome. The Arabidopsis chromosome was mitotically stable in hybrid cells following cell division, and initially existed as a translocated chromosome. During culture, the translocated chromosomes then converted to two types of independent plant chromosomes without human DNA sequences, with reproducibility. One pair of localization signals of CENP-A, a marker of functional centromeres was detected in the Arabidopsis genomic region in independent plant chromosomes. These results suggest that the chromosome maintenance system was conserved between human and plants. Furthermore, the expression of plant endogenous genes was observed in the hybrid cells, implicating that the plant chromosomal region existed as euchromatin in a human cell background and the gene expression system is conserved between two organisms. The present study suggests that the essential chromosome functions are conserved between evolutionarily distinct organisms such as humans and plants. Systematic analyses of hybrid cells may lead to the production of a shuttle vector between animal and plant, and a platform for the genome writing.
Collapse
Affiliation(s)
| | | | | | | | - Kiichi Fukui
- Department
of Biotechnology, Graduate School of Engineering, Osaka University, 565-0871, Osaka, Japan
| | | |
Collapse
|
37
|
Xie J, Yang X, Song Y, Du Q, Li Y, Chen J, Zhang D. Adaptive evolution and functional innovation of Populus-specific recently evolved microRNAs. THE NEW PHYTOLOGIST 2017; 213:206-219. [PMID: 27277139 DOI: 10.1111/nph.14046] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/11/2016] [Accepted: 05/03/2016] [Indexed: 05/23/2023]
Abstract
Lineage-specific microRNAs (miRNAs) undergo rapid turnover during evolution; however, their origin and functional importance have remained controversial. Here, we examine the origin, evolution, and potential roles in local adaptation of Populus-specific miRNAs, which originated after the recent salicoid-specific, whole-genome duplication. RNA sequencing was used to generate extensive, comparable miRNA and gene expression data for six tissues. A natural population of Populus trichocarpa and closely related species were used to study the divergence rates, evolution, and adaptive variation of miRNAs. MiRNAs that originated in 5' untranslated regions had higher expression levels and their expression showed high correlation with their host genes. Compared with conserved miRNAs, a significantly higher proportion of Populus-specific miRNAs appear to target genes that were duplicated in salicoids. Examination of single nucleotide polymorphisms in Populus-specific miRNA precursors showed high amounts of population differentiation. We also characterized the newly emerged MIR6445 family, which could trigger the production of phased small interfering RNAs from NAC mRNAs, which encode a transcription factor with primary roles in a variety of plant developmental processes. Together, these observations provide evolutionary insights into the birth and potential roles of Populus-specific miRNAs in genome maintenance, local adaptation, and functional innovation.
Collapse
Affiliation(s)
- Jianbo Xie
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, No. 35 Qinghua East Road, Beijing, 100083, China
- Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education, College of Biological Sciences and Technology, Beijing Forestry University, No. 35 Qinghua East Road, Beijing, 100083, China
| | - Xiaohui Yang
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, No. 35 Qinghua East Road, Beijing, 100083, China
- Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education, College of Biological Sciences and Technology, Beijing Forestry University, No. 35 Qinghua East Road, Beijing, 100083, China
| | - Yuepeng Song
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, No. 35 Qinghua East Road, Beijing, 100083, China
- Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education, College of Biological Sciences and Technology, Beijing Forestry University, No. 35 Qinghua East Road, Beijing, 100083, China
| | - Qingzhang Du
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, No. 35 Qinghua East Road, Beijing, 100083, China
- Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education, College of Biological Sciences and Technology, Beijing Forestry University, No. 35 Qinghua East Road, Beijing, 100083, China
| | - Ying Li
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, No. 35 Qinghua East Road, Beijing, 100083, China
- Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education, College of Biological Sciences and Technology, Beijing Forestry University, No. 35 Qinghua East Road, Beijing, 100083, China
| | - Jinhui Chen
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, No. 35 Qinghua East Road, Beijing, 100083, China
- Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education, College of Biological Sciences and Technology, Beijing Forestry University, No. 35 Qinghua East Road, Beijing, 100083, China
| | - Deqiang Zhang
- National Engineering Laboratory for Tree Breeding, College of Biological Sciences and Technology, Beijing Forestry University, No. 35 Qinghua East Road, Beijing, 100083, China
- Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education, College of Biological Sciences and Technology, Beijing Forestry University, No. 35 Qinghua East Road, Beijing, 100083, China
| |
Collapse
|
38
|
Kryuchkova-Mostacci N, Robinson-Rechavi M. Tissue-Specificity of Gene Expression Diverges Slowly between Orthologs, and Rapidly between Paralogs. PLoS Comput Biol 2016; 12:e1005274. [PMID: 28030541 PMCID: PMC5193323 DOI: 10.1371/journal.pcbi.1005274] [Citation(s) in RCA: 42] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2016] [Accepted: 11/26/2016] [Indexed: 11/18/2022] Open
Abstract
The ortholog conjecture implies that functional similarity between orthologous genes is higher than between paralogs. It has been supported using levels of expression and Gene Ontology term analysis, although the evidence was rather weak and there were also conflicting reports. In this study on 12 species we provide strong evidence of high conservation in tissue-specificity between orthologs, in contrast to low conservation between within-species paralogs. This allows us to shed a new light on the evolution of gene expression patterns. While there have been several studies of the correlation of expression between species, little is known about the evolution of tissue-specificity itself. Ortholog tissue-specificity is strongly conserved between all tetrapod species, with the lowest Pearson correlation between mouse and frog at r = 0.66. Tissue-specificity correlation decreases strongly with divergence time. Paralogs in human show much lower conservation, even for recent Primate-specific paralogs. When both paralogs from ancient whole genome duplication tissue-specific paralogs are tissue-specific, it is often to different tissues, while other tissue-specific paralogs are mostly specific to the same tissue. The same patterns are observed using human or mouse as focal species, and are robust to choices of datasets and of thresholds. Our results support the following model of evolution: in the absence of duplication, tissue-specificity evolves slowly, and tissue-specific genes do not change their main tissue of expression; after small-scale duplication the less expressed paralog loses the ancestral specificity, leading to an immediate difference between paralogs; over time, both paralogs become more broadly expressed, but remain poorly correlated. Finally, there is a small number of paralog pairs which stay tissue-specific with the same main tissue of expression, for at least 300 million years.
Collapse
Affiliation(s)
- Nadezda Kryuchkova-Mostacci
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Marc Robinson-Rechavi
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
- Swiss Institute of Bioinformatics, Lausanne, Switzerland
| |
Collapse
|
39
|
Liu L, Chang Y, Yang T, Noren DP, Long B, Kornblau S, Qutub A, Ye J. Evolution-informed modeling improves outcome prediction for cancers. Evol Appl 2016; 10:68-76. [PMID: 28035236 PMCID: PMC5192825 DOI: 10.1111/eva.12417] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2016] [Accepted: 08/17/2016] [Indexed: 12/19/2022] Open
Abstract
Despite wide applications of high-throughput biotechnologies in cancer research, many biomarkers discovered by exploring large-scale omics data do not provide satisfactory performance when used to predict cancer treatment outcomes. This problem is partly due to the overlooking of functional implications of molecular markers. Here, we present a novel computational method that uses evolutionary conservation as prior knowledge to discover bona fide biomarkers. Evolutionary selection at the molecular level is nature's test on functional consequences of genetic elements. By prioritizing genes that show significant statistical association and high functional impact, our new method reduces the chances of including spurious markers in the predictive model. When applied to predicting therapeutic responses for patients with acute myeloid leukemia and to predicting metastasis for patients with prostate cancers, the new method gave rise to evolution-informed models that enjoyed low complexity and high accuracy. The identified genetic markers also have significant implications in tumor progression and embrace potential drug targets. Because evolutionary conservation can be estimated as a gene-specific, position-specific, or allele-specific parameter on the nucleotide level and on the protein level, this new method can be extended to apply to miscellaneous "omics" data to accelerate biomarker discoveries.
Collapse
Affiliation(s)
- Li Liu
- Department of Biomedical Informatics Arizona State University Tempe AZ USA
| | - Yung Chang
- School of Life Science Arizona State University Tempe AZ USA
| | - Tao Yang
- Department of Computer Science and Engineering Arizona State University Tempe AZ USA
| | - David P Noren
- Department of Bioengineering Rice University Houston TX USA
| | - Byron Long
- Department of Bioengineering Rice University Houston TX USA
| | - Steven Kornblau
- The University of Texas MD Anderson Cancer Center Houston TX USA
| | - Amina Qutub
- Department of Bioengineering Rice University Houston TX USA
| | - Jieping Ye
- Department of Computational Medicine and Bioinformatics University of Michigan Ann Arbor MI USA
| |
Collapse
|
40
|
Chen W, Xia X, Song N, Wang Y, Zhu H, Deng W, Kong Q, Pan X, Qin C. Cross-Species Analysis of Gene Expression and Function in Prefrontal Cortex, Hippocampus and Striatum. PLoS One 2016; 11:e0164295. [PMID: 27716781 PMCID: PMC5055290 DOI: 10.1371/journal.pone.0164295] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2016] [Accepted: 09/22/2016] [Indexed: 01/08/2023] Open
Abstract
Background Mouse has been extensively used as a tool for investigating the onset and development of human neurological disorders. As a first step to construct a transgenic mouse model of human brain lesions, it is of fundamental importance to clarify the similarity and divergence of genetic background between non-diseased human and mouse brain tissues. Methods We systematically compared, based on large scale integrated microarray data, the transcriptomes of three anatomically distinct brain regions; prefrontal cortex (PFC), hippocampus (HIP) and striatum (STR), across human and mouse. The widely used DAVID web server was used to decipher the biological functions of the highly expressed genes that were identified using a previously reported approach. Venn analysis was used to depict the overlapping ratios of the notably enriched biological process (BP) terms (one-tailed Fisher’s exact test and Benjamini correction; adjusted p < 0.01) between two brain tissues. GOSemSim, an R package, was selected to perform GO semantic similarity analysis. Next, we adjusted signal intensities of orthologous genes by the total signals in all samples within species, and used one minus Pearson’s correlation coefficient to assess the expression distance. Hierarchical clustering and principal component analysis (PCA) were selected for expression pattern analysis. Lineage specific expressed orthologous genes were identified by comparison of the most extreme sub-datasets across species and further verified using reverse transcription PCR (RT-PCR) and quantitative real-time PCR (qRT-PCR). Results We found that the number of the significantly enriched BP terms of the highly expressed genes in human brain regions is larger than that in mouse corresponding brain regions. The mainly involved BP terms in human brain tissues associated with protein-membrane targeting and selenium metabolism are species-specific. The overlapping ratios of all the significantly enriched BP terms between any two brain tissues across species are lower than that within species, but the pairwise semantic similarities are very high between any two brain tissues from either human or mouse. Hierarchical clustering analysis shows the biological functions of the highly expressed genes in brain tissues are more consistent within species than interspecies; whereas it shows the expression patterns of orthologous genes are evidently conserved between human and mouse equivalent brain tissues. In addition, we identified four orthologous genes (COX5B, WIF1, SLC4A10 and PLA2G7) that are species-specific, which have been widely studied and confirmed to be closely linked with neuro- physiological and pathological functions. Conclusion Our study highlights the similarities and divergences in gene function and expression between human and mouse corresponding brain regions, including PFC, HIP and STR.
Collapse
Affiliation(s)
- Wei Chen
- Institute of Laboratory Animal Science, Chinese Academy of Medical Sciences (CAMS) & Comparative Medicine Centre, Peking Union Medical Collage (PUMC), Beijing, P.R. China
| | - Xiayu Xia
- Institute of Laboratory Animal Science, Chinese Academy of Medical Sciences (CAMS) & Comparative Medicine Centre, Peking Union Medical Collage (PUMC), Beijing, P.R. China
| | - Nan Song
- Institute of Laboratory Animal Science, Chinese Academy of Medical Sciences (CAMS) & Comparative Medicine Centre, Peking Union Medical Collage (PUMC), Beijing, P.R. China
| | - Ying Wang
- Institute of Laboratory Animal Science, Chinese Academy of Medical Sciences (CAMS) & Comparative Medicine Centre, Peking Union Medical Collage (PUMC), Beijing, P.R. China
| | - Hua Zhu
- Institute of Laboratory Animal Science, Chinese Academy of Medical Sciences (CAMS) & Comparative Medicine Centre, Peking Union Medical Collage (PUMC), Beijing, P.R. China
| | - Wei Deng
- Institute of Laboratory Animal Science, Chinese Academy of Medical Sciences (CAMS) & Comparative Medicine Centre, Peking Union Medical Collage (PUMC), Beijing, P.R. China
| | - Qi Kong
- Institute of Laboratory Animal Science, Chinese Academy of Medical Sciences (CAMS) & Comparative Medicine Centre, Peking Union Medical Collage (PUMC), Beijing, P.R. China
| | - Xianmin Pan
- Institute of Laboratory Animal Science, Chinese Academy of Medical Sciences (CAMS) & Comparative Medicine Centre, Peking Union Medical Collage (PUMC), Beijing, P.R. China
- Ministry of Education, The Key Laboratory of Bioinformatics, School of Life Sciences, Tsinghua University, Beijing, P.R. China
| | - Chuan Qin
- Institute of Laboratory Animal Science, Chinese Academy of Medical Sciences (CAMS) & Comparative Medicine Centre, Peking Union Medical Collage (PUMC), Beijing, P.R. China
- * E-mail:
| |
Collapse
|
41
|
Breschi A, Djebali S, Gillis J, Pervouchine DD, Dobin A, Davis CA, Gingeras TR, Guigó R. Gene-specific patterns of expression variation across organs and species. Genome Biol 2016; 17:151. [PMID: 27391956 PMCID: PMC4937605 DOI: 10.1186/s13059-016-1008-y] [Citation(s) in RCA: 57] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2016] [Accepted: 06/14/2016] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND A comparison of transcriptional profiles derived from different tissues in a given species or among different species assumes that commonalities reflect evolutionarily conserved programs and that differences reflect species or tissue responses to environmental conditions or developmental program staging. Apparently conflicting results have been published regarding whether organ-specific transcriptional patterns dominate over species-specific patterns, or vice versa, making it unclear to what extent the biology of a given organism can be extrapolated to another. These studies have in common that they treat the transcriptomes monolithically, implicitly ignoring that each gene is likely to have a specific pattern of transcriptional variation across organs and species. RESULTS We use linear models to quantify this pattern. We find a continuum in the spectrum of expression variation: the expression of some genes varies considerably across species and little across organs, and simply reflects evolutionary distance. At the other extreme are genes whose expression varies considerably across organs and little across species; these genes are much more likely to be associated with diseases than are genes whose expression varies predominantly across species. CONCLUSIONS Whether transcriptomes, when considered globally, cluster preferentially according to one component or the other may not be a property of the transcriptomes, but rather a consequence of the dominant behavior of a subset of genes. Therefore, the values of the components of the variance of expression for each gene could become a useful resource when planning, interpreting, and extrapolating experimental data from mouse to humans.
Collapse
Affiliation(s)
- Alessandra Breschi
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, Barcelona, 08003, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Sarah Djebali
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, Barcelona, 08003, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
- GenPhySE, Université de Toulouse, INRA, INPT, INP-ENVT, Castanet Tolosan, France
| | - Jesse Gillis
- Cold Spring Harbor LaboratoryCold Spring Harbor, NY, 11742, USA
| | - Dmitri D Pervouchine
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, Barcelona, 08003, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Alex Dobin
- Cold Spring Harbor LaboratoryCold Spring Harbor, NY, 11742, USA
| | - Carrie A Davis
- Cold Spring Harbor LaboratoryCold Spring Harbor, NY, 11742, USA
| | | | - Roderic Guigó
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, Barcelona, 08003, Spain.
- Universitat Pompeu Fabra (UPF), Barcelona, Spain.
| |
Collapse
|
42
|
Chen G, Chen J, Yang J, Chen L, Qu X, Shi C, Ning B, Shi L, Tong W, Zhao Y, Zhang M, Shi T. Significant variations in alternative splicing patterns and expression profiles between human-mouse orthologs in early embryos. SCIENCE CHINA-LIFE SCIENCES 2016; 60:178-188. [PMID: 27378339 DOI: 10.1007/s11427-015-0348-5] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/18/2015] [Accepted: 02/11/2016] [Indexed: 02/05/2023]
Abstract
Human and mouse orthologs are expected to have similar biological functions; however, many discrepancies have also been reported. We systematically compared human and mouse orthologs in terms of alternative splicing patterns and expression profiles. Human-mouse orthologs are divergent in alternative splicing, as human orthologs could generally encode more isoforms than their mouse orthologs. In early embryos, exon skipping is far more common with human orthologs, whereas constitutive exons are more prevalent with mouse orthologs. This may correlate with divergence in expression of splicing regulators. Orthologous expression similarities are different in distinct embryonic stages, with the highest in morula. Expression differences for orthologous transcription factor genes could play an important role in orthologous expression discordance. We further detected largely orthologous divergence in differential expression between distinct embryonic stages. Collectively, our study uncovers significant orthologous divergence from multiple aspects, which may result in functional differences and dynamics between human-mouse orthologs during embryonic development.
Collapse
Affiliation(s)
- Geng Chen
- The Center for Bioinformatics and Computational Biology, Shanghai Key Laboratory of Regulatory Biology, the Institute of Biomedical Sciences and School of Life Sciences, East China Normal University, Shanghai, 200241, China.,Center for Pharmacogenomics, School of Pharmacy, Fudan University, Shanghai, 201203, China
| | - Jiwei Chen
- The Center for Bioinformatics and Computational Biology, Shanghai Key Laboratory of Regulatory Biology, the Institute of Biomedical Sciences and School of Life Sciences, East China Normal University, Shanghai, 200241, China
| | - Jianmin Yang
- The Center for Bioinformatics and Computational Biology, Shanghai Key Laboratory of Regulatory Biology, the Institute of Biomedical Sciences and School of Life Sciences, East China Normal University, Shanghai, 200241, China
| | - Long Chen
- The Center for Bioinformatics and Computational Biology, Shanghai Key Laboratory of Regulatory Biology, the Institute of Biomedical Sciences and School of Life Sciences, East China Normal University, Shanghai, 200241, China
| | - Xiongfei Qu
- The Center for Bioinformatics and Computational Biology, Shanghai Key Laboratory of Regulatory Biology, the Institute of Biomedical Sciences and School of Life Sciences, East China Normal University, Shanghai, 200241, China
| | - Caiping Shi
- The Center for Bioinformatics and Computational Biology, Shanghai Key Laboratory of Regulatory Biology, the Institute of Biomedical Sciences and School of Life Sciences, East China Normal University, Shanghai, 200241, China
| | - Baitang Ning
- National Center for Toxicological Research, US Food and Drug Administration, Jefferson, AR, 72079, USA
| | - Leming Shi
- Center for Pharmacogenomics, School of Pharmacy, Fudan University, Shanghai, 201203, China.,National Center for Toxicological Research, US Food and Drug Administration, Jefferson, AR, 72079, USA
| | - Weida Tong
- National Center for Toxicological Research, US Food and Drug Administration, Jefferson, AR, 72079, USA
| | - Yongxiang Zhao
- Biological Targeting Diagnosis and Therapy Research Center, Guangxi Medical University, Nanning, 530021, China.
| | - Meixia Zhang
- Department of Ophthalmology, West China Hospital, Sichuan University, Chengdu, 610041, China.
| | - Tieliu Shi
- The Center for Bioinformatics and Computational Biology, Shanghai Key Laboratory of Regulatory Biology, the Institute of Biomedical Sciences and School of Life Sciences, East China Normal University, Shanghai, 200241, China.
| |
Collapse
|
43
|
Contina A, Bridge ES, Kelly JF. Exploring novel candidate genes from the Mouse Genome Informatics database: Potential implications for avian migration research. Integr Zool 2016; 11:240-9. [PMID: 27061206 DOI: 10.1111/1749-4877.12199] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
To search for genes associated with migratory phenotypes in songbirds, we selected candidate genes through annotations from the Mouse Genome Informatics database and assembled an extensive candidate-gene library. Then, we implemented a next-generation sequencing approach to obtain DNA sequences from the Painted Bunting genome. We focused on those sequences that were conserved across avian species and that aligned with candidate genes in our mouse library. We genotyped short sequence repeats from the following candidate genes: ADRA1d, ANKRD17, CISH and MYH7. We studied the possible correlations between allelic variations occurring in these novel candidate migration genes and avian migratory phenotypes available from the published literature. We found that allele variation at MYH7 correlated with a calculated index of speed of migration (km/day) across 11 species of songbirds. We highlight the potential of the Mouse Genome Informatics database in providing new candidate genes that might play a crucial role in regulating migration in birds and possibly in other taxa. Our research effort shows the benefits and limitations of working with extensive genomic datasets and offers a snapshot of the challenges related to cross-species validation in behavioral and molecular ecology studies.
Collapse
Affiliation(s)
- Andrea Contina
- Oklahoma Biological Survey, University of Oklahoma, Norman, Oklahoma, USA
| | - Eli S Bridge
- Oklahoma Biological Survey, University of Oklahoma, Norman, Oklahoma, USA
| | - Jeffrey F Kelly
- Oklahoma Biological Survey, University of Oklahoma, Norman, Oklahoma, USA.,Department of Biology, University of Oklahoma, Norman, Oklahoma, USA
| |
Collapse
|
44
|
A comparison of human and mouse gene co-expression networks reveals conservation and divergence at the tissue, pathway and disease levels. BMC Evol Biol 2015; 15:259. [PMID: 26589719 PMCID: PMC4654840 DOI: 10.1186/s12862-015-0534-7] [Citation(s) in RCA: 68] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2015] [Accepted: 11/09/2015] [Indexed: 12/25/2022] Open
Abstract
Background A deeper understanding of differences and similarities in transcriptional regulation between species can uncover important information about gene functions and the role of genes in disease. Deciphering such patterns between mice and humans is especially important since mice play an essential role in biomedical research. Results Here, in order to characterize evolutionary changes between humans and mice, we compared gene co-expression maps to evaluate the conservation of co-expression. We show that the conservation of co-expression connectivity of homologous genes is negatively correlated with molecular evolution rates, as expected. Then we investigated evolutionary aspects of gene sets related to functions, tissues, pathways and diseases. Genes expressed in the testis, eye and skin, and those associated with regulation of transcription, olfaction, PI3K signalling, response to virus and bacteria were more divergent between mice and humans in terms of co-expression connectivity. Surprisingly, a deeper investigation of the PI3K signalling cascade revealed that its divergence is caused by the most crucial genes of this pathway, such as mTOR and AKT2. On the other hand, our analysis revealed that genes expressed in the brain and in the bone, and those associated with cell adhesion, cell cycle, DNA replication and DNA repair are most strongly conserved in terms of co-expression network connectivity as well as having a lower rate of duplication events. Genes involved in lipid metabolism and genes specific to blood showed a signature of increased co-expression connectivity in the mouse. In terms of diseases, co-expression connectivity of genes related to metabolic disorders is the most strongly conserved between mice and humans and tumor-related genes the most divergent. Conclusions This work contributes to discerning evolutionary patterns between mice and humans in terms of gene interactions. Conservation of co-expression is a powerful approach to identify gene targets and processes with potential similarity and divergence between mice and humans, which has implications for drug testing and other studies employing the mouse as a model organism. Electronic supplementary material The online version of this article (doi:10.1186/s12862-015-0534-7) contains supplementary material, which is available to authorized users.
Collapse
|
45
|
Delhomme N, Sundström G, Zamani N, Lantz H, Lin YC, Hvidsten TR, Höppner MP, Jern P, Van de Peer Y, Lundeberg J, Grabherr MG, Street NR. Serendipitous Meta-Transcriptomics: The Fungal Community of Norway Spruce (Picea abies). PLoS One 2015; 10:e0139080. [PMID: 26413905 PMCID: PMC4586145 DOI: 10.1371/journal.pone.0139080] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2015] [Accepted: 09/09/2015] [Indexed: 11/18/2022] Open
Abstract
After performing de novo transcript assembly of >1 billion RNA-Sequencing reads obtained from 22 samples of different Norway spruce (Picea abies) tissues that were not surface sterilized, we found that assembled sequences captured a mix of plant, lichen, and fungal transcripts. The latter were likely expressed by endophytic and epiphytic symbionts, indicating that these organisms were present, alive, and metabolically active. Here, we show that these serendipitously sequenced transcripts need not be considered merely as contamination, as is common, but that they provide insight into the plant’s phyllosphere. Notably, we could classify these transcripts as originating predominantly from Dothideomycetes and Leotiomycetes species, with functional annotation of gene families indicating active growth and metabolism, with particular regards to glucose intake and processing, as well as gene regulation.
Collapse
Affiliation(s)
- Nicolas Delhomme
- Department of Plant Physiology, Umeå Plant Science Centre, Umeå University, Umeå, Sweden
| | - Görel Sundström
- Department of Plant Physiology, Umeå Plant Science Centre, Umeå University, Umeå, Sweden
- Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
| | - Neda Zamani
- Department of Plant Physiology, Umeå Plant Science Centre, Umeå University, Umeå, Sweden
- Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
| | - Henrik Lantz
- Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
- Bioinformatics Infrastructure for Life Sciences (BILS), Uppsala, Sweden
| | - Yao-Cheng Lin
- Department of Plant Systems Biology (VIB) and Department of Plant Biotechnology and Bioinformatics (Ghent University), Ghent, Belgium
| | - Torgeir R. Hvidsten
- Department of Plant Physiology, Umeå Plant Science Centre, Umeå University, Umeå, Sweden
- Department of Chemistry, Biotechnology and Food Science, Norwegian University of Life Sciences, Ås, Norway
| | - Marc P. Höppner
- Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
- Bioinformatics Infrastructure for Life Sciences (BILS), Uppsala, Sweden
| | - Patric Jern
- Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
| | - Yves Van de Peer
- Department of Plant Systems Biology (VIB) and Department of Plant Biotechnology and Bioinformatics (Ghent University), Ghent, Belgium
- Genomics Research Institute, University of Pretoria, Pretoria, South Africa
| | - Joakim Lundeberg
- School of Biotechnology, Science for Life Laboratory, KTH Royal Institute of Technology, Solna, Sweden
| | - Manfred G. Grabherr
- Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
- Bioinformatics Infrastructure for Life Sciences (BILS), Uppsala, Sweden
| | - Nathaniel R. Street
- Department of Plant Physiology, Umeå Plant Science Centre, Umeå University, Umeå, Sweden
- * E-mail:
| |
Collapse
|
46
|
Raherison ESM, Giguère I, Caron S, Lamara M, MacKay JJ. Modular organization of the white spruce (Picea glauca) transcriptome reveals functional organization and evolutionary signatures. THE NEW PHYTOLOGIST 2015; 207:172-187. [PMID: 25728802 PMCID: PMC5024012 DOI: 10.1111/nph.13343] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/27/2014] [Accepted: 01/18/2015] [Indexed: 05/13/2023]
Abstract
Transcript profiling has shown the molecular bases of several biological processes in plants but few studies have developed an understanding of overall transcriptome variation. We investigated transcriptome structure in white spruce (Picea glauca), aiming to delineate its modular organization and associated functional and evolutionary attributes. Microarray analyses were used to: identify and functionally characterize groups of co-expressed genes; investigate expressional and functional diversity of vascular tissue preferential genes which were conserved among Picea species, and identify expression networks underlying wood formation. We classified 22 857 genes as variable (79%; 22 coexpression groups) or invariant (21%) by profiling across several vegetative tissues. Modular organization and complex transcriptome restructuring among vascular tissue preferential genes was revealed by their assignment to coexpression groups with partially overlapping profiles and partially distinct functions. Integrated analyses of tissue-based and temporally variable profiles identified secondary xylem gene networks, showed their remodelling over a growing season and identified PgNAC-7 (no apical meristerm (NAM), Arabidopsis transcription activation factor (ATAF) and cup-shaped cotyledon (CUC) transcription factor 007 in Picea glauca) as a major hub gene specific to earlywood formation. Reference profiling identified comprehensive, statistically robust coexpressed groups, revealing that modular organization underpins the evolutionary conservation of the transcriptome structure.
Collapse
Affiliation(s)
- Elie S M Raherison
- Center for Forest Research and Institute for Integrative and Systems Biology, Université Laval, Québec, QC, G1V 0A6, Canada
| | - Isabelle Giguère
- Center for Forest Research and Institute for Integrative and Systems Biology, Université Laval, Québec, QC, G1V 0A6, Canada
| | - Sébastien Caron
- Center for Forest Research and Institute for Integrative and Systems Biology, Université Laval, Québec, QC, G1V 0A6, Canada
| | - Mebarek Lamara
- Center for Forest Research and Institute for Integrative and Systems Biology, Université Laval, Québec, QC, G1V 0A6, Canada
| | - John J MacKay
- Department of Plant Sciences, University of Oxford, OX1 3RB, Oxford, UK
| |
Collapse
|
47
|
Roux J, Rosikiewicz M, Robinson-Rechavi M. What to compare and how: Comparative transcriptomics for Evo-Devo. JOURNAL OF EXPERIMENTAL ZOOLOGY PART B-MOLECULAR AND DEVELOPMENTAL EVOLUTION 2015; 324:372-82. [PMID: 25864439 PMCID: PMC4949521 DOI: 10.1002/jez.b.22618] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Received: 11/07/2014] [Accepted: 02/19/2015] [Indexed: 12/30/2022]
Abstract
Evolutionary developmental biology has grown historically from the capacity to relate patterns of evolution in anatomy to patterns of evolution of expression of specific genes, whether between very distantly related species, or very closely related species or populations. Scaling up such studies by taking advantage of modern transcriptomics brings promising improvements, allowing us to estimate the overall impact and molecular mechanisms of convergence, constraint or innovation in anatomy and development. But it also presents major challenges, including the computational definitions of anatomical homology and of organ function, the criteria for the comparison of developmental stages, the annotation of transcriptomics data to proper anatomical and developmental terms, and the statistical methods to compare transcriptomic data between species to highlight significant conservation or changes. In this article, we review these challenges, and the ongoing efforts to address them, which are emerging from bioinformatics work on ontologies, evolutionary statistics, and data curation, with a focus on their implementation in the context of the development of our database Bgee (http://bgee.org). J. Exp. Zool. (Mol. Dev. Evol.) 324B: 372–382, 2015. © 2015 The Authors. J. Exp. Zool. (Mol. Dev. Evol.) published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Julien Roux
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland.,Swiss Institute of Bioinformatics, Lausanne, Switzerland.,Department of Human Genetics, University of Chicago, Chicago, Illinois
| | - Marta Rosikiewicz
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland.,Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Marc Robinson-Rechavi
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland.,Swiss Institute of Bioinformatics, Lausanne, Switzerland
| |
Collapse
|
48
|
Unraveling the association between mRNA expressions and mutant phenotypes in a genome-wide assessment of mice. Proc Natl Acad Sci U S A 2015; 112:4707-12. [PMID: 25825715 DOI: 10.1073/pnas.1415046112] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
High-throughput gene expression profiling has revealed substantial leaky and extraneous transcription of eukaryotic genes, challenging the perceptions that transcription is strictly regulated and that changes in transcription have phenotypic consequences. To assess the functional implications of mRNA transcription directly, we analyzed mRNA expression data derived from microarrays, RNA-sequencing, and in situ hybridization, together with phenotype data of mouse mutants as a proxy of gene function at the tissue level. The results indicated that despite the presence of widespread ectopic transcription, mRNA expression and mutant phenotypes of mammalian genes or tissues remain associated. The expression-phenotype association at the gene level was particularly strong for tissue-specific genes, and the association could be underestimated due to data insufficiency and incomprehensive phenotyping of mouse mutants; the strength of expression-phenotype association at the tissue level depended on tissue functions. Mutations on genes expressed at higher levels or expressed at earlier embryonic stages more often result in abnormal phenotypes in the tissues where they are expressed. The mRNA expression profiles that have stronger associations with their phenotype profiles tend to be more evolutionarily conserved, indicating that the evolution of transcriptome and the evolution of phenome are coupled. Therefore, mutations resulting in phenotypic aberrations in expressed tissues are more likely to occur in highly transcribed genes, tissue-specific genes, genes expressed during early embryonic stages, or genes with evolutionarily conserved mRNA expression profiles.
Collapse
|
49
|
Characterization of the humanTCAM1Ppseudogene and its activation by a potential dual promoter-enhancer: Comparison with a protein-coding mouse orthologue. FEBS Lett 2015; 589:540-7. [DOI: 10.1016/j.febslet.2015.01.023] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2014] [Revised: 01/14/2015] [Accepted: 01/15/2015] [Indexed: 11/17/2022]
|
50
|
Diament A, Pinter RY, Tuller T. Three-dimensional eukaryotic genomic organization is strongly correlated with codon usage expression and function. Nat Commun 2014; 5:5876. [PMID: 25510862 DOI: 10.1038/ncomms6876] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2014] [Accepted: 11/17/2014] [Indexed: 01/08/2023] Open
Abstract
It has been shown that the distribution of genes in eukaryotic genomes is not random; however, formerly reported relations between gene function and genomic organization were relatively weak. Previous studies have demonstrated that codon usage bias is related to all stages of gene expression and to protein function. Here we apply a novel tool for assessing functional relatedness, codon usage frequency similarity (CUFS), which measures similarity between genes in terms of codon and amino acid usage. By analyzing chromosome conformation capture data, describing the three-dimensional (3D) conformation of the DNA, we show that the functional similarity between genes captured by CUFS is directly and very strongly correlated with their 3D distance in Saccharomyces cerevisiae, Schizosaccharomyces pombe, Arabidopsis thaliana, mouse and human. This emphasizes the importance of three-dimensional genomic localization in eukaryotes and indicates that codon usage is tightly linked to genome architecture.
Collapse
Affiliation(s)
- Alon Diament
- Department of Biomedical Engineering, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Ron Y Pinter
- Department of Computer Science, Technion-Israel Institute of Technology, Haifa 32000, Israel
| | - Tamir Tuller
- 1] Department of Biomedical Engineering, Tel Aviv University, Tel Aviv 6997801, Israel [2] The Sagol School of Neuroscience, Tel Aviv University, Tel Aviv 6997801, Israel
| |
Collapse
|