1
|
Krishnakant Kushwaha S, Wu Y, Leonardo Avila H, Anand A, Sicheritz-Pontén T, Millard A, Amol Marathe S, Nobrega FL. Comprehensive blueprint of Salmonella genomic plasticity identifies hotspots for pathogenicity genes. PLoS Biol 2024; 22:e3002746. [PMID: 39110680 PMCID: PMC11305592 DOI: 10.1371/journal.pbio.3002746] [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: 01/15/2024] [Accepted: 07/10/2024] [Indexed: 08/10/2024] Open
Abstract
Understanding the dynamic evolution of Salmonella is vital for effective bacterial infection management. This study explores the role of the flexible genome, organised in regions of genomic plasticity (RGP), in shaping the pathogenicity of Salmonella lineages. Through comprehensive genomic analysis of 12,244 Salmonella spp. genomes covering 2 species, 6 subspecies, and 46 serovars, we uncover distinct integration patterns of pathogenicity-related gene clusters into RGP, challenging traditional views of gene distribution. These RGP exhibit distinct preferences for specific genomic spots, and the presence or absence of such spots across Salmonella lineages profoundly shapes strain pathogenicity. RGP preferences are guided by conserved flanking genes surrounding integration spots, implicating their involvement in regulatory networks and functional synergies with integrated gene clusters. Additionally, we emphasise the multifaceted contributions of plasmids and prophages to the pathogenicity of diverse Salmonella lineages. Overall, this study provides a comprehensive blueprint of the pathogenicity potential of Salmonella. This unique insight identifies genomic spots in nonpathogenic lineages that hold the potential for harbouring pathogenicity genes, providing a foundation for predicting future adaptations and developing targeted strategies against emerging human pathogenic strains.
Collapse
Affiliation(s)
- Simran Krishnakant Kushwaha
- Department of Biological Sciences, Birla Institute of Technology & Science (BITS), Pilani, Rajasthan, India
- School of Biological Sciences, University of Southampton, Southampton, United Kingdom
| | - Yi Wu
- School of Biological Sciences, University of Southampton, Southampton, United Kingdom
| | - Hugo Leonardo Avila
- Laboratory for Applied Science and Technology in Health, Instituto Carlos Chagas, FIOCRUZ Paraná, Brazil
| | - Abhirath Anand
- Department of Computer Sciences and Information Systems, Birla Institute of Technology & Science (BITS), Pilani, Rajasthan, India
| | - Thomas Sicheritz-Pontén
- Center for Evolutionary Hologenomics, Globe Institute, University of Copenhagen, Copenhagen, Denmark
- Centre of Excellence for Omics-Driven Computational Biodiscovery (COMBio), AIMST University, Bedong, Kedah, Malaysia
| | - Andrew Millard
- Centre for Phage Research, Department of Genetics and Genome Biology, University of Leicester, Leicester, United Kingdom
| | - Sandhya Amol Marathe
- Department of Biological Sciences, Birla Institute of Technology & Science (BITS), Pilani, Rajasthan, India
| | - Franklin L. Nobrega
- School of Biological Sciences, University of Southampton, Southampton, United Kingdom
| |
Collapse
|
2
|
Guo Y, Wang T, Lu X, Li W, Lv X, Peng Q, Zhang J, Gao J, Hu M. Comparative genome-wide analysis of circular RNAs in Brassica napus L.: target-site versus non-target-site resistance to herbicide stress. TAG. THEORETICAL AND APPLIED GENETICS. THEORETISCHE UND ANGEWANDTE GENETIK 2024; 137:176. [PMID: 38969812 DOI: 10.1007/s00122-024-04678-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/31/2023] [Accepted: 06/15/2024] [Indexed: 07/07/2024]
Abstract
Circular RNAs (circRNAs), a class of non-coding RNA molecules, are recognized for their unique functions; however, their responses to herbicide stress in Brassica napus remain unclear. In this study, the role of circRNAs in response to herbicide treatment was investigated in two rapeseed cultivars: MH33, which confers non-target-site resistance (NTSR), and EM28, which exhibits target-site resistance (TSR). The genome-wide circRNA profiles of herbicide-stressed and non-stressed seedlings were analyzed. The findings indicate that NTSR seedlings exhibited a greater abundance of circRNAs, shorter lengths of circRNAs and their parent genes, and more diverse functions of parent genes compared with TSR seedlings. Compared to normal-growth plants, the herbicide-stressed group exhibited similar trends in the number of circRNAs, functions of parent genes, and differentially expressed circRNAs as observed in NTSR seedlings. In addition, a greater number of circRNAs that function as competing microRNA (miRNA) sponges were identified in the herbicide stress and NTSR groups compared to the normal-growth and TSR groups, respectively. The differentially expressed circRNAs were validated by qPCR. The differntially expressed circRNA-miRNA networks were predicted, and the mRNAs targeted by these miRNAs were annotated. Our results suggest that circRNAs play a crucial role in responding to herbicide stress, exhibiting distinct responses between NTSR and TSR in rapeseed. These findings offer valuable insights into the mechanisms underlying herbicide resistance in rapeseed.
Collapse
Affiliation(s)
- Yue Guo
- Institute of Industrial Crops, Jiangsu Academy of Agricultural Sciences, Nanjing, China
| | - Ting Wang
- Institute of Industrial Crops, Jiangsu Academy of Agricultural Sciences, Nanjing, China
| | - Xinyu Lu
- Institute of Industrial Crops, Jiangsu Academy of Agricultural Sciences, Nanjing, China
| | - Weilong Li
- Institute of Industrial Crops, Jiangsu Academy of Agricultural Sciences, Nanjing, China
| | - Xinlei Lv
- Institute of Industrial Crops, Jiangsu Academy of Agricultural Sciences, Nanjing, China
| | - Qi Peng
- Institute of Industrial Crops, Jiangsu Academy of Agricultural Sciences, Nanjing, China
| | - Jiefu Zhang
- Institute of Industrial Crops, Jiangsu Academy of Agricultural Sciences, Nanjing, China
| | - Jianqin Gao
- Institute of Industrial Crops, Jiangsu Academy of Agricultural Sciences, Nanjing, China
| | - Maolong Hu
- Institute of Industrial Crops, Jiangsu Academy of Agricultural Sciences, Nanjing, China.
| |
Collapse
|
3
|
McCoy MJ, Fire AZ. Parallel gene size and isoform expansion of ancient neuronal genes. Curr Biol 2024; 34:1635-1645.e3. [PMID: 38460513 PMCID: PMC11043017 DOI: 10.1016/j.cub.2024.02.021] [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/22/2023] [Revised: 12/16/2023] [Accepted: 02/11/2024] [Indexed: 03/11/2024]
Abstract
How nervous systems evolved is a central question in biology. A diversity of synaptic proteins is thought to play a central role in the formation of specific synapses leading to nervous system complexity. The largest animal genes, often spanning hundreds of thousands of base pairs, are known to be enriched for expression in neurons at synapses and are frequently mutated or misregulated in neurological disorders and diseases. Although many of these genes have been studied independently in the context of nervous system evolution and disease, general principles underlying their parallel evolution remain unknown. To investigate this, we directly compared orthologous gene sizes across eukaryotes. By comparing relative gene sizes within organisms, we identified a distinct class of large genes with origins predating the diversification of animals and, in many cases, the emergence of neurons as dedicated cell types. We traced this class of ancient large genes through evolution and found orthologs of the large synaptic genes potentially driving the immense complexity of metazoan nervous systems, including in humans and cephalopods. Moreover, we found that while these genes are evolving under strong purifying selection, as demonstrated by low dN/dS ratios, they have simultaneously grown larger and gained the most isoforms in animals. This work provides a new lens through which to view this distinctive class of large and multi-isoform genes and demonstrates how intrinsic genomic properties, such as gene length, can provide flexibility in molecular evolution and allow groups of genes and their host organisms to evolve toward complexity.
Collapse
Affiliation(s)
- Matthew J McCoy
- Department of Pathology, Stanford University School of Medicine, 300 Pasteur Drive, Stanford, CA 94305, USA.
| | - Andrew Z Fire
- Department of Pathology, Stanford University School of Medicine, 300 Pasteur Drive, Stanford, CA 94305, USA; Department of Genetics, Stanford University School of Medicine, 300 Pasteur Drive, Stanford, CA 94305, USA.
| |
Collapse
|
4
|
Gupta S, Fernandes R, Natarajan S, Jose NP, Giri J, Dahal S. Comparative evaluation of arch form among the Nepalese population: A morphological study. J Oral Maxillofac Pathol 2024; 28:111-118. [PMID: 38800435 PMCID: PMC11126270 DOI: 10.4103/jomfp.jomfp_280_23] [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: 06/25/2023] [Revised: 09/20/2023] [Accepted: 10/09/2023] [Indexed: 05/29/2024] Open
Abstract
Aims The study aims to identify sexual dimorphic features in the arch patterns based on tooth arrangement patterns and the maxillary and mandibular arches using Euclidean Distance Matrix Analysis (EDMA). Settings and Design A total of 96 Nepalese subjects, aged 18 to 25 were assessed using casts and photographs. Materials and Methods Thirteen landmarks representing the most facial portions of the proximal contact areas on the maxillary and mandibular casts were digitised. Seventy-eight possible, Euclidean distances between the 13 landmarks were calculated using the Analysis ToolPak of Microsoft Excel®. The male-to-female ratios of the corresponding distances were computed and ratios were compared to evaluate the arch form for variation in the genders, among the Nepalese population. Statistical Analysis Used Microsoft Excel Analysis ToolPak and SPSS 20.0 (IBM Chicago) were used to perform EDMA and an independent t-test to compare the significant differences between the two genders. Results The maxillary arch's largest ratio (1.008179001) was discovered near the location of the right and left lateral incisors, indicating that the anterior region may have experienced the greatest change. The posterior-molar region is where the smallest ratio was discovered, suggesting less variation. At the intercanine region, female arches were wider than male ones; however, at the interpremolar and intermolar sections, they were similar in width. Females' maxillary arches were discovered to be bigger antero-posteriorly than those of males. The highest ratio (1.014336113) in the mandibular arch was discovered at the intermolar area, suggesting that males had a larger mandibular posterior arch morphology. At the intercanine area, the breadth of the arch form was greater in males and nearly the same in females at the interpremolar and intermolar regions. Female mandibular arch forms were also discovered to be longer than those of males from the anterior to the posterior. Conclusions The male and female arches in the Nepalese population were inferred to be different in size and shape. With references to the landmarks demonstrating such a shift, the EDMA established objectively the presence of square arch forms in Nepali males and tapering arch forms in Nepalese females.
Collapse
Affiliation(s)
- Simran Gupta
- Intern, Manipal College of Dental Sciences, Mangalore, Manipal Academy of Higher Education, Manipal, Karnataka, India
| | - Rhea Fernandes
- Intern, Manipal College of Dental Sciences, Mangalore, Manipal Academy of Higher Education, Manipal, Karnataka, India
| | - Srikant Natarajan
- Department of Forensic Odontology, Manipal College of Dental Sciences, Mangalore, Manipal Academy of Higher Education, Manipal, Karnataka, India
| | - Nidhin P. Jose
- Department of Orthodontics and Dentofacial Orthopaedics, Manipal College of Dental Sciences, Mangalore, Manipal Academy of Higher Education, Manipal, Karnataka, India
| | - Jamal Giri
- Department of Orthodontics, B.P. Koirala Institute of Health Sciences, Dharan, Nepal, India
| | - Samarika Dahal
- Department of Oral Pathology and Forensic Dentistry, Maharajgunj Medical Campus, Institute of Medicine, Nepal, India
| |
Collapse
|
5
|
Zanin O, Eastham M, Winczura K, Ashe M, Martinez-Nunez RT, Hebenstreit D, Grzechnik P. Ceg1 depletion reveals mechanisms governing degradation of non-capped RNAs in Saccharomyces cerevisiae. Commun Biol 2023; 6:1112. [PMID: 37919390 PMCID: PMC10622555 DOI: 10.1038/s42003-023-05495-6] [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: 01/20/2023] [Accepted: 10/20/2023] [Indexed: 11/04/2023] Open
Abstract
Most functional eukaryotic mRNAs contain a 5' 7-methylguanosine (m7G) cap. Although capping is essential for many biological processes including mRNA processing, export and translation, the fate of uncapped transcripts has not been studied extensively. Here, we employed fast nuclear depletion of the capping enzymes in Saccharomyces cerevisiae to uncover the turnover of the transcripts that failed to be capped. We show that although the degradation of cap-deficient mRNA is dominant, the levels of hundreds of non-capped mRNAs increase upon depletion of the capping enzymes. Overall, the abundance of non-capped mRNAs is inversely correlated to the expression levels, altogether resembling the effects observed in cells lacking the cytoplasmic 5'-3' exonuclease Xrn1 and indicating differential degradation fates of non-capped mRNAs. The inactivation of the nuclear 5'-3' exonuclease Rat1 does not rescue the non-capped mRNA levels indicating that Rat1 is not involved in their degradation and consequently, the lack of the capping does not affect the distribution of RNA Polymerase II on the chromatin. Our data indicate that the cap presence is essential to initiate the Xrn1-dependent degradation of mRNAs underpinning the role of 5' cap in the Xrn1-dependent buffering of the cellular mRNA levels.
Collapse
Affiliation(s)
- Onofrio Zanin
- School of Biosciences, University of Birmingham, Edgbaston, Birmingham, B15 2TT, UK
- School of Immunology & Microbial Sciences, King's College London, Guy's Campus, London, SE1 9RT, UK
| | - Matthew Eastham
- Division of Molecular and Cellular Function, School of Biological Sciences, University of Manchester, Oxford Road, Manchester, M13 9PT, UK
| | - Kinga Winczura
- Division of Molecular and Cellular Function, School of Biological Sciences, University of Manchester, Oxford Road, Manchester, M13 9PT, UK
| | - Mark Ashe
- Division of Molecular and Cellular Function, School of Biological Sciences, University of Manchester, Oxford Road, Manchester, M13 9PT, UK
| | - Rocio T Martinez-Nunez
- School of Immunology & Microbial Sciences, King's College London, Guy's Campus, London, SE1 9RT, UK
| | | | - Pawel Grzechnik
- Division of Molecular and Cellular Function, School of Biological Sciences, University of Manchester, Oxford Road, Manchester, M13 9PT, UK.
| |
Collapse
|
6
|
Belcher LJ, Dewar AE, Hao C, Ghoul M, West SA. Signatures of kin selection in a natural population of the bacteria Bacillus subtilis. Evol Lett 2023; 7:315-330. [PMID: 37829498 PMCID: PMC10565896 DOI: 10.1093/evlett/qrad029] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2023] [Revised: 06/14/2023] [Accepted: 07/07/2023] [Indexed: 10/14/2023] Open
Abstract
Laboratory experiments have suggested that bacteria perform a range of cooperative behaviors, which are favored because they are directed toward relatives (kin selection). However, there is a lack of evidence for cooperation and kin selection in natural bacterial populations. Molecular population genetics offers a promising method to study natural populations because the theory predicts that kin selection will lead to relaxed selection, which will result in increased polymorphism and divergence at cooperative genes. Examining a natural population of Bacillus subtilis, we found consistent evidence that putatively cooperative traits have higher polymorphism and greater divergence than putatively private traits expressed at the same rate. In addition, we were able to eliminate alternative explanations for these patterns and found more deleterious mutations in genes controlling putatively cooperative traits. Overall, our results suggest that cooperation is favored by kin selection, with an average relatedness of r = .79 between interacting individuals.
Collapse
Affiliation(s)
| | - Anna E Dewar
- Department of Biology, University of Oxford, Oxford, United Kingdom
| | - Chunhui Hao
- Department of Biology, University of Oxford, Oxford, United Kingdom
| | - Melanie Ghoul
- Department of Biology, University of Oxford, Oxford, United Kingdom
| | - Stuart A West
- Department of Biology, University of Oxford, Oxford, United Kingdom
| |
Collapse
|
7
|
Huang Q, Szklarczyk D, Wang M, Simonovic M, von Mering C. PaxDb 5.0: Curated Protein Quantification Data Suggests Adaptive Proteome Changes in Yeasts. Mol Cell Proteomics 2023; 22:100640. [PMID: 37659604 PMCID: PMC10551891 DOI: 10.1016/j.mcpro.2023.100640] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2023] [Revised: 08/25/2023] [Accepted: 08/30/2023] [Indexed: 09/04/2023] Open
Abstract
The "Protein Abundances Across Organisms" database (PaxDb) is an integrative metaresource dedicated to protein abundance levels, in tissue-specific or whole-organism proteomes. PaxDb focuses on computing best-estimate abundances for proteins in normal/healthy contexts and expresses abundance values for each protein in "parts per million" in relation to all other protein molecules in the cell. The uniform data reprocessing, quality scoring, and integrated orthology relations have made PaxDb one of the preferred tools for comparisons between individual datasets, tissues, or organisms. In describing the latest version 5.0 of PaxDb, we particularly emphasize the data integration from various types of raw data and how we expanded the number of organisms and tissue groups as well as the proteome coverage. The current collection of PaxDb includes 831 original datasets from 170 species, including 22 Archaea, 81 Bacteria, and 67 Eukaryota. Apart from detailing the data update, we also present a comparative analysis of the human proteome subset of PaxDb against the two most widely used human proteome data resources: Human Protein Atlas and Genotype-Tissue Expression. Lastly, through our protein abundance data, we reveal an evolutionary trend in the usage of sulfur-containing amino acids in the proteomes of Fungi.
Collapse
Affiliation(s)
- Qingyao Huang
- Department of Molecular Life Sciences and Swiss Institute of Bioinformatics, University of Zurich, Zurich, Switzerland
| | - Damian Szklarczyk
- Department of Molecular Life Sciences and Swiss Institute of Bioinformatics, University of Zurich, Zurich, Switzerland
| | - Mingcong Wang
- Department of Molecular Life Sciences and Swiss Institute of Bioinformatics, University of Zurich, Zurich, Switzerland
| | - Milan Simonovic
- Department of Molecular Life Sciences and Swiss Institute of Bioinformatics, University of Zurich, Zurich, Switzerland
| | - Christian von Mering
- Department of Molecular Life Sciences and Swiss Institute of Bioinformatics, University of Zurich, Zurich, Switzerland.
| |
Collapse
|
8
|
Cheng Y, Chi Y, Sun L, Wang GZ. Dominant constraints on the evolution of rhythmic gene expression. Comput Struct Biotechnol J 2023; 21:4301-4311. [PMID: 37692081 PMCID: PMC10492206 DOI: 10.1016/j.csbj.2023.08.035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2023] [Revised: 08/21/2023] [Accepted: 08/31/2023] [Indexed: 09/12/2023] Open
Abstract
Although the individual transcriptional regulators of the core circadian clock are distinct among different organisms, the autoregulatory feedback loops they form are conserved. This unified design principle explains how daily physiological activities oscillate across species. However, it is unknown whether analogous design principles govern the gene expression output of circadian clocks. In this study, we performed a comparative analysis of rhythmic gene expression in eight diverse species and identified four common distribution patterns of cycling gene expression across these species. We hypothesized that the maintenance of reduced energetic costs constrains the evolution of rhythmic gene expression. Our large-scale computational simulations support this hypothesis by showing that selection against high-energy expenditure completely regenerates all cycling gene patterns. Moreover, we find that the peaks of rhythmic expression have been subjected to this type of selective pressure. The results suggest that selective pressure from circadian regulation efficiently removes unnecessary gene products from the transcriptome, thereby significantly impacting its evolutionary path.
Collapse
Affiliation(s)
| | | | | | - Guang-Zhong Wang
- CAS Key Laboratory of Computational Biology, Shanghai Institute of Nutrition and Health, University of Chinese Academy of Sciences, Chinese Academy of Sciences, Shanghai 200031, China
| |
Collapse
|
9
|
McCoy MJ, Fire AZ. Ancient origins of complex neuronal genes. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.03.28.534655. [PMID: 37034725 PMCID: PMC10081198 DOI: 10.1101/2023.03.28.534655] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/19/2023]
Abstract
How nervous systems evolved is a central question in biology. An increasing diversity of synaptic proteins is thought to play a central role in the formation of specific synapses leading to nervous system complexity. The largest animal genes, often spanning millions of base pairs, are known to be enriched for expression in neurons at synapses and are frequently mutated or misregulated in neurological disorders and diseases. While many of these genes have been studied independently in the context of nervous system evolution and disease, general principles underlying their parallel evolution remain unknown. To investigate this, we directly compared orthologous gene sizes across eukaryotes. By comparing relative gene sizes within organisms, we identified a distinct class of large genes with origins predating the diversification of animals and in many cases the emergence of dedicated neuronal cell types. We traced this class of ancient large genes through evolution and found orthologs of the large synaptic genes driving the immense complexity of metazoan nervous systems, including in humans and cephalopods. Moreover, we found that while these genes are evolving under strong purifying selection as demonstrated by low dN/dS scores, they have simultaneously grown larger and gained the most isoforms in animals. This work provides a new lens through which to view this distinctive class of large and multi-isoform genes and demonstrates how intrinsic genomic properties, such as gene length, can provide flexibility in molecular evolution and allow groups of genes and their host organisms to evolve toward complexity.
Collapse
Affiliation(s)
- Matthew J. McCoy
- Department of Pathology, Stanford University School of Medicine, Stanford, CA 94305, USA
- Whitman Center, Marine Biological Laboratory, Woods Hole, MA 02543, USA
| | - Andrew Z. Fire
- Department of Pathology, Stanford University School of Medicine, Stanford, CA 94305, USA
- Department of Genetics, Stanford University School of Medicine, Stanford, CA 94305, USA
| |
Collapse
|
10
|
Nevers Y, Glover NM, Dessimoz C, Lecompte O. Protein length distribution is remarkably uniform across the tree of life. Genome Biol 2023; 24:135. [PMID: 37291671 PMCID: PMC10251718 DOI: 10.1186/s13059-023-02973-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2022] [Accepted: 05/16/2023] [Indexed: 06/10/2023] Open
Abstract
BACKGROUND In every living species, the function of a protein depends on its organization of structural domains, and the length of a protein is a direct reflection of this. Because every species evolved under different evolutionary pressures, the protein length distribution, much like other genomic features, is expected to vary across species but has so far been scarcely studied. RESULTS Here we evaluate this diversity by comparing protein length distribution across 2326 species (1688 bacteria, 153 archaea, and 485 eukaryotes). We find that proteins tend to be on average slightly longer in eukaryotes than in bacteria or archaea, but that the variation of length distribution across species is low, especially compared to the variation of other genomic features (genome size, number of proteins, gene length, GC content, isoelectric points of proteins). Moreover, most cases of atypical protein length distribution appear to be due to artifactual gene annotation, suggesting the actual variation of protein length distribution across species is even smaller. CONCLUSIONS These results open the way for developing a genome annotation quality metric based on protein length distribution to complement conventional quality measures. Overall, our findings show that protein length distribution between living species is more uniform than previously thought. Furthermore, we also provide evidence for a universal selection on protein length, yet its mechanism and fitness effect remain intriguing open questions.
Collapse
Affiliation(s)
- Yannis Nevers
- Department of Computational Biology, University of Lausanne, Lausanne, Switzerland.
- Swiss Institute for Bioinformatics, University of Lausanne, Lausanne, Switzerland.
| | - Natasha M Glover
- Department of Computational Biology, University of Lausanne, Lausanne, Switzerland
- Swiss Institute for Bioinformatics, University of Lausanne, Lausanne, Switzerland
| | - Christophe Dessimoz
- Department of Computational Biology, University of Lausanne, Lausanne, Switzerland
- Swiss Institute for Bioinformatics, University of Lausanne, Lausanne, Switzerland
- Department of Computer Science, University College London, London, UK
- Centre for Life's Origins and Evolution, Department of Genetics, Evolution and Environment, University College London, London, UK
| | - Odile Lecompte
- Department of Computer Science, Centre de Recherche en Biomédecine de Strasbourg, ICube, UMR 7357, University of Strasbourg, CNRS, Strasbourg, France
| |
Collapse
|
11
|
Roberts M, Josephs EB. Weaker selection on genes with treatment-specific expression consistent with a limit on plasticity evolution in Arabidopsis thaliana. Genetics 2023; 224:iyad074. [PMID: 37094602 PMCID: PMC10484170 DOI: 10.1093/genetics/iyad074] [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: 03/06/2023] [Revised: 03/06/2023] [Accepted: 04/07/2023] [Indexed: 04/26/2023] Open
Abstract
Differential gene expression between environments often underlies phenotypic plasticity. However, environment-specific expression patterns are hypothesized to relax selection on genes, and thus limit plasticity evolution. We collated over 27 terabases of RNA-sequencing data on Arabidopsis thaliana from over 300 peer-reviewed studies and 200 treatment conditions to investigate this hypothesis. Consistent with relaxed selection, genes with more treatment-specific expression have higher levels of nucleotide diversity and divergence at nonsynonymous sites but lack stronger signals of positive selection. This result persisted even after controlling for expression level, gene length, GC content, the tissue specificity of expression, and technical variation between studies. Overall, our investigation supports the existence of a hypothesized trade-off between the environment specificity of a gene's expression and the strength of selection on said gene in A. thaliana. Future studies should leverage multiple genome-scale datasets to tease apart the contributions of many variables in limiting plasticity evolution.
Collapse
Affiliation(s)
- Miles Roberts
- Genetics and Genome Sciences Program, Michigan State University, East Lansing, MI 48824, USA
| | - Emily B Josephs
- Department of Plant Biology, Michigan State University, East Lansing, MI 48824, USA
- Ecology, Evolution, and Behavior Program, Michigan State University, East Lansing, MI 48824, USA
| |
Collapse
|
12
|
Mimoso CA, Adelman K. U1 snRNP increases RNA Pol II elongation rate to enable synthesis of long genes. Mol Cell 2023; 83:1264-1279.e10. [PMID: 36965480 PMCID: PMC10135401 DOI: 10.1016/j.molcel.2023.03.002] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2022] [Revised: 02/06/2023] [Accepted: 02/28/2023] [Indexed: 03/27/2023]
Abstract
The expansion of introns within mammalian genomes poses a challenge for the production of full-length messenger RNAs (mRNAs), with increasing evidence that these long AT-rich sequences present obstacles to transcription. Here, we investigate RNA polymerase II (RNAPII) elongation at high resolution in mammalian cells and demonstrate that RNAPII transcribes faster across introns. Moreover, we find that this acceleration requires the association of U1 snRNP (U1) with the elongation complex at 5' splice sites. The role of U1 to stimulate elongation rate through introns reduces the frequency of both premature termination and transcriptional arrest, thereby dramatically increasing RNA production. We further show that changes in RNAPII elongation rate due to AT content and U1 binding explain previous reports of pausing or termination at splice junctions and the edge of CpG islands. We propose that U1-mediated acceleration of elongation has evolved to mitigate the risks that long AT-rich introns pose to transcript completion.
Collapse
Affiliation(s)
- Claudia A Mimoso
- Department of Biological Chemistry and Molecular Pharmacology, Blavatnik Institute, Harvard Medical School, Boston, MA 02115, USA
| | - Karen Adelman
- Department of Biological Chemistry and Molecular Pharmacology, Blavatnik Institute, Harvard Medical School, Boston, MA 02115, USA; Ludwig Center at Harvard, Boston, MA 02115, USA; Broad Institute of MIT and Harvard, Cambridge, MA 02142, USA.
| |
Collapse
|
13
|
Mielczarek M, Frąszczak M, Zielak-Steciwko AE, Nowak B, Hofman B, Pierścińska J, Kruszyński W, Szyda J. An effect of large-scale deletions and duplications on transcript expression. Funct Integr Genomics 2022; 23:19. [PMID: 36564645 PMCID: PMC9789009 DOI: 10.1007/s10142-022-00946-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2022] [Revised: 12/14/2022] [Accepted: 12/15/2022] [Indexed: 12/25/2022]
Abstract
Since copy number variants (CNVs) have been recognized as an important source of genetic and transcriptomic variation, we aimed to characterize the impact of CNVs located within coding, intergenic, upstream, and downstream gene regions on the expression of transcripts. Regions in which deletions occurred most often were introns, while duplications in coding regions. The transcript expression was lower for deleted coding (P = 0.008) and intronic regions (P = 1.355 × 10-10), but it was not changed in the case of upstream and downstream gene regions (P = 0.085). Moreover, the expression was decreased if duplication occurred in the coding region (P = 8.318 × 10-5). Furthermore, a negative correlation (r = - 0.27) between transcript length and its expression was observed. The correlation between the percent of deleted/duplicated transcript and transcript expression level was not significant for all concerned genomic regions in five out of six animals. The exceptions were deletions in coding regions (P = 0.004) and duplications in introns (P = 0.01) in one individual. CNVs in coding (deletions, duplications) and intronic (deletions) regions are important modulators of transcripts by reducing their expression level. We hypothesize that deletions imply severe consequences by interrupting genes. The negative correlation between the size of the transcript and its expression level found in this study is consistent with the hypothesis that selection favours shorter introns and a moderate number of exons in highly expressed genes. This may explain the transcript expression reduction by duplications. We did not find the correlation between the size of deletions/duplications and transcript expression level suggesting that expression is modulated by CNVs regardless of their size.
Collapse
Affiliation(s)
- Magda Mielczarek
- Wroclaw University of Environmental and Life Sciences, Kozuchowska 7, 51-631, Wroclaw, Poland.
| | - Magdalena Frąszczak
- Wroclaw University of Environmental and Life Sciences, Kozuchowska 7, 51-631, Wroclaw, Poland
| | - Anna E Zielak-Steciwko
- Wroclaw University of Environmental and Life Sciences, Kozuchowska 7, 51-631, Wroclaw, Poland
| | - Błażej Nowak
- Wroclaw University of Environmental and Life Sciences, Kozuchowska 7, 51-631, Wroclaw, Poland
| | - Bartłomiej Hofman
- Wroclaw University of Environmental and Life Sciences, Kozuchowska 7, 51-631, Wroclaw, Poland
| | - Jagoda Pierścińska
- Wroclaw University of Environmental and Life Sciences, Kozuchowska 7, 51-631, Wroclaw, Poland
| | - Wojciech Kruszyński
- Wroclaw University of Environmental and Life Sciences, Kozuchowska 7, 51-631, Wroclaw, Poland
| | - Joanna Szyda
- Wroclaw University of Environmental and Life Sciences, Kozuchowska 7, 51-631, Wroclaw, Poland
| |
Collapse
|
14
|
Wen R, Yang D, Wang K, Zan X. Characterization of two full-length Araneus ventricosus major ampullate silk protein genes. Int J Biol Macromol 2022; 213:297-304. [PMID: 35654219 DOI: 10.1016/j.ijbiomac.2022.05.178] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2022] [Revised: 05/25/2022] [Accepted: 05/26/2022] [Indexed: 11/24/2022]
Abstract
Major ampullate silk is noted for its great tensile strength and extensibility. The impressive material properties of major ampullate silk result from their component proteins that encoded by members of the spidroin (spider fibroin) gene family. Although the major ampullate spidroin type has evolved multiple variants within specific-species, most sequences are fragmented. Here, we present two complete major ampullate spidroin genes from the orb-weaving spider Araneus ventricosus. Due to the abundant GPG motifs in their repetitive region, the two MaSp genes were grouped in MaSp2 subclass and named MaSp2C and MaSp2D, respectively. Analysis of the full-length gene sequences reveals that both of them include a single enormous exon (10,851 bp for MaSp2C and 8640 bp for MaSp2D) that mainly translates into a central repetitive region containing multiple amino acid motifs that can be organized into five ensemble types. We use gene-specific PCR primers to search the cDNA from major ampullate glands and find evidence for alternative splicing of MaSp2D transcripts into a minor spliceoform lacking the entire repetitive domain as well as the partial terminal regions. Our results not only provide new templates for protein-based materials with tailored properties, but suggest gene and transcriptional diversity of major ampullate silk.
Collapse
Affiliation(s)
- Rui Wen
- School of Ophthalmology and Optometry, Eye Hospital, Wenzhou Medical University, Wenzhou, Zhejiang Province 325035, China; Oujiang Laboratory, Wenzhou Key Laboratory of Perioperative Medicine, Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325000, China; School of Biomedical Engineering, Wenzhou Medical University, Wenzhou, Zhejiang Province 325035, China.
| | - Dong Yang
- School of Ophthalmology and Optometry, Eye Hospital, Wenzhou Medical University, Wenzhou, Zhejiang Province 325035, China
| | - Kangkang Wang
- Oujiang Laboratory, Wenzhou Key Laboratory of Perioperative Medicine, Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325000, China
| | - Xingjie Zan
- School of Ophthalmology and Optometry, Eye Hospital, Wenzhou Medical University, Wenzhou, Zhejiang Province 325035, China; Oujiang Laboratory, Wenzhou Key Laboratory of Perioperative Medicine, Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325000, China; School of Biomedical Engineering, Wenzhou Medical University, Wenzhou, Zhejiang Province 325035, China.
| |
Collapse
|
15
|
Khandia R, Saeed M, Alharbi AM, Ashraf GM, Greig NH, Kamal MA. Codon Usage Bias Correlates With Gene Length in Neurodegeneration Associated Genes. Front Neurosci 2022; 16:895607. [PMID: 35860292 PMCID: PMC9289476 DOI: 10.3389/fnins.2022.895607] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2022] [Accepted: 06/08/2022] [Indexed: 11/13/2022] Open
Abstract
Codon usage analysis is a crucial part of molecular characterization and is used to determine the factors affecting the evolution of a gene. The length of a gene is an important parameter that affects the characteristics of the gene, such as codon usage, compositional parameters, and sometimes, its functions. In the present study, we investigated the association of various parameters related to codon usage with the length of genes. Gene expression is affected by nucleotide disproportion. In sixty genes related to neurodegenerative disorders, the G nucleotide was the most abundant and the T nucleotide was the least. The nucleotide T exhibited a significant association with the length of the gene at both the overall compositional level and the first and second codon positions. Codon usage bias (CUB) of these genes was affected by pyrimidine and keto skews. Gene length was found to be significantly correlated with codon bias in neurodegeneration associated genes. In gene segments with lengths below 1,200 bp and above 2,400 bp, CUB was positively associated with length. Relative synonymous CUB, which is another measure of CUB, showed that codons TTA, GTT, GTC, TCA, GGT, and GGA exhibited a positive association with length, whereas codons GTA, AGC, CGT, CGA, and GGG showed a negative association. GC-ending codons were preferred over AT-ending codons. Overall analysis indicated that the association between CUB and length varies depending on the segment size; however, CUB of 1,200–2,000 bp gene segments appeared not affected by gene length. In synopsis, analysis suggests that length of the genes correlates with various imperative molecular signatures including A/T nucleotide disproportion and codon choices. In the present study we additionally evaluated various molecular features and their correlation with different indices of codon usage, like the Codon Adaptation Index (CAI) and Relative Dynonymous Codon Usage (RSCU) of codons. We also considered the impact of gene fragment size on different molecular features in genes related to neurodegeneration. This analysis will aid our understanding of and in potentially modulating gene expression in cases of defective gene functioning in clinical settings.
Collapse
Affiliation(s)
- Rekha Khandia
- Department of Biochemistry and Genetics, Barkatullah University, Bhopal, India
- *Correspondence: Rekha Khandia, ;
| | - Mohd. Saeed
- Department of Biology, College of Sciences, University of Hail, Hail, Saudi Arabia
| | - Ahmed M. Alharbi
- Department of Clinical Laboratory Sciences, College of Applied Medical Sciences, University of Hail, Hail, Saudi Arabia
| | - Ghulam Md. Ashraf
- Pre-clinical Research Unit, King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia
- Department of Medical Laboratory Sciences, Faculty of Applied Medical Sciences, King Abdulaziz University, Jeddah, Saudi Arabia
| | - Nigel H. Greig
- Drug Design and Development Section, Translational Gerontology Branch, Intramural Research Program National Institute on Aging, NIH, Baltimore, MD, United States
| | - Mohammad Amjad Kamal
- Institutes for Systems Genetics, Frontiers Science Center for Disease-Related Molecular Network, West China Hospital, Sichuan University, Chengdu, China
- King Fahd Medical Research Center, King Abdulaziz University, Jeddah, Saudi Arabia
- Department of Pharmacy, Faculty of Allied Health Sciences, Daffodil International University, Dhaka, Bangladesh
- Enzymoics, Novel Global Community Educational Foundation, Hebersham, NSW, Australia
| |
Collapse
|
16
|
Tao W, Bian J, Tang M, Zeng Y, Luo R, Ke Q, Li T, Li Y, Cui L. Genomic insights into positive selection during barley domestication. BMC PLANT BIOLOGY 2022; 22:267. [PMID: 35641942 PMCID: PMC9158214 DOI: 10.1186/s12870-022-03655-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/19/2022] [Accepted: 05/23/2022] [Indexed: 06/15/2023]
Abstract
BACKGROUND Cultivated barley (Hordeum vulgare) is widely used in animal feed, beverages, and foods and has become a model crop for molecular evolutionary studies. Few studies have examined the evolutionary fates of different types of genes in barley during the domestication process. RESULTS The rates of nonsynonymous substitution (Ka) to synonymous substitution (Ks) were calculated by comparing orthologous genes in different barley groups (wild vs. landrace and landrace vs. improved cultivar). The rates of evolution, properties, expression patterns, and diversity of positively selected genes (PSGs) and negatively selected genes (NSGs) were compared. PSGs evolved more rapidly, possessed fewer exons, and had lower GC content than NSGs; they were also shorter and had shorter intron, exon, and first exon lengths. Expression levels were lower, the tissue specificity of expression was higher, and codon usage bias was weaker for PSGs than for NSGs. Nucleotide diversity analysis revealed that PSGs have undergone a more severe genetic bottleneck than NSGs. Several candidate PSGs were involved in plant growth and development, which might make them as excellent targets for the molecular breeding of barley. CONCLUSIONS Our comprehensive analysis of the evolutionary, structural, and functional divergence between PSGs and NSGs in barley provides new insight into the evolutionary trajectory of barley during domestication. Our findings also aid future functional studies of PSGs in barley.
Collapse
Affiliation(s)
- Wenjing Tao
- College of Bioscience and Engineering, Jiangxi Agricultural University, Nanchang, Jiangxi, 330045 China
| | - Jianxin Bian
- Peking University Institute of Advanced Agricultural Sciences, Weifang, Shandong, 261325 China
| | - Minqiang Tang
- College of Forestry, Hainan University, Haikou, Hainan, 570228 China
| | - Yan Zeng
- College of Bioscience and Engineering, Jiangxi Agricultural University, Nanchang, Jiangxi, 330045 China
| | - Ruihan Luo
- College of Bioscience and Engineering, Jiangxi Agricultural University, Nanchang, Jiangxi, 330045 China
| | - Qinglin Ke
- College of Bioscience and Engineering, Jiangxi Agricultural University, Nanchang, Jiangxi, 330045 China
| | - Tingting Li
- College of Bioscience and Engineering, Jiangxi Agricultural University, Nanchang, Jiangxi, 330045 China
| | - Yihan Li
- College of Bioscience and Engineering, Jiangxi Agricultural University, Nanchang, Jiangxi, 330045 China
| | - Licao Cui
- College of Bioscience and Engineering, Jiangxi Agricultural University, Nanchang, Jiangxi, 330045 China
| |
Collapse
|
17
|
Ramezaniaghdam M, Nahdi ND, Reski R. Recombinant Spider Silk: Promises and Bottlenecks. Front Bioeng Biotechnol 2022; 10:835637. [PMID: 35350182 PMCID: PMC8957953 DOI: 10.3389/fbioe.2022.835637] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2021] [Accepted: 02/01/2022] [Indexed: 02/02/2023] Open
Abstract
Spider silk threads have exceptional mechanical properties such as toughness, elasticity and low density, which reach maximum values compared to other fibre materials. They are superior even compared to Kevlar and steel. These extraordinary properties stem from long length and specific protein structures. Spider silk proteins can consist of more than 20,000 amino acids. Polypeptide stretches account for more than 90% of the whole protein, and these domains can be repeated more than a hundred times. Each repeat unit has a specific function resulting in the final properties of the silk. These properties make them attractive for innovative material development for medical or technical products as well as cosmetics. However, with livestock breeding of spiders it is not possible to reach high volumes of silk due to the cannibalistic behaviour of these animals. In order to obtain spider silk proteins (spidroins) on a large scale, recombinant production is attempted in various expression systems such as plants, bacteria, yeasts, insects, silkworms, mammalian cells and animals. For viable large-scale production, cost-effective and efficient production systems are needed. This review describes the different types of spider silk, their proteins and structures and discusses the production of these difficult-to-express proteins in different host organisms with an emphasis on plant systems.
Collapse
Affiliation(s)
- Maryam Ramezaniaghdam
- Plant Biotechnology, Faculty of Biology, University of Freiburg, Freiburg, Germany
- Cluster of Excellence livMatS at FIT – Freiburg Center for Interactive Materials and Bioinspired Technologies, University of Freiburg, Freiburg, Germany
| | - Nadia D. Nahdi
- Plant Biotechnology, Faculty of Biology, University of Freiburg, Freiburg, Germany
| | - Ralf Reski
- Plant Biotechnology, Faculty of Biology, University of Freiburg, Freiburg, Germany
- Cluster of Excellence livMatS at FIT – Freiburg Center for Interactive Materials and Bioinspired Technologies, University of Freiburg, Freiburg, Germany
| |
Collapse
|
18
|
Gaunt ER, Digard P. Compositional biases in RNA viruses: Causes, consequences and applications. WILEY INTERDISCIPLINARY REVIEWS. RNA 2022; 13:e1679. [PMID: 34155814 PMCID: PMC8420353 DOI: 10.1002/wrna.1679] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/08/2021] [Revised: 05/29/2021] [Accepted: 05/31/2021] [Indexed: 01/05/2023]
Abstract
If each of the four nucleotides were represented equally in the genomes of viruses and the hosts they infect, each base would occur at a frequency of 25%. However, this is not observed in nature. Similarly, the order of nucleotides is not random (e.g., in the human genome, guanine follows cytosine at a frequency of ~0.0125, or a quarter the number of times predicted by random representation). Codon usage and codon order are also nonrandom. Furthermore, nucleotide and codon biases vary between species. Such biases have various drivers, including cellular proteins that recognize specific patterns in nucleic acids, that once triggered, induce mutations or invoke intrinsic or innate immune responses. In this review we examine the types of compositional biases identified in viral genomes and current understanding of the evolutionary mechanisms underpinning these trends. Finally, we consider the potential for large scale synonymous recoding strategies to engineer RNA virus vaccines, including those with pandemic potential, such as influenza A virus and Severe Acute Respiratory Syndrome Coronavirus Virus 2. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Evolution and Genomics > Computational Analyses of RNA RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition.
Collapse
Affiliation(s)
- Eleanor R. Gaunt
- Department of Infection and ImmunityThe Roslin Institute, The University of EdinburghEdinburghUK
| | - Paul Digard
- Department of Infection and ImmunityThe Roslin Institute, The University of EdinburghEdinburghUK
| |
Collapse
|
19
|
Kin selection for cooperation in natural bacterial populations. Proc Natl Acad Sci U S A 2022; 119:2119070119. [PMID: 35193981 PMCID: PMC8892524 DOI: 10.1073/pnas.2119070119] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/06/2021] [Indexed: 01/03/2023] Open
Abstract
Bacteria secrete many molecules outside the cell, where they provide benefits to other cells. One potential reason for producing these “public goods” is that they benefit closely related cells that share the gene for cooperation (kin selection). While many laboratory studies have supported this hypothesis, there is a lack of evidence that kin selection favors cooperation in natural populations. We examined bacterial genomes from the environment and used population genetics theory to analyze the DNA sequences. Our analyses suggest that public goods cooperation has indeed been favored by kin selection in natural populations. Bacteria produce a range of molecules that are secreted from the cell and can provide a benefit to the local population of cells. Laboratory experiments have suggested that these “public goods” molecules represent a form of cooperation, favored because they benefit closely related cells (kin selection). However, there is a relative lack of data demonstrating kin selection for cooperation in natural populations of bacteria. We used molecular population genetics to test for signatures of kin selection at the genomic level in natural populations of the opportunistic pathogen Pseudomonas aeruginosa. We found consistent evidence from multiple traits that genes controlling putatively cooperative traits have higher polymorphism and greater divergence and are more likely to harbor deleterious mutations relative to genes controlling putatively private traits, which are expressed at similar rates. These patterns suggest that cooperative traits are controlled by kin selection, and we estimate that the relatedness for social interactions in P. aeruginosa is r = 0.84. More generally, our results demonstrate how molecular population genetics can be used to study the evolution of cooperation in natural populations.
Collapse
|
20
|
Khandia R, Ali Khan A, Alexiou A, Povetkin SN, Nikolaevna VM. Codon Usage Analysis of Pro-Apoptotic Bim Gene Isoforms. J Alzheimers Dis 2022; 86:1711-1725. [DOI: 10.3233/jad-215691] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Background: Bim is a Bcl-2 homology 3 (BH3)-only proteins, a group of pro-apoptotic proteins involved in physiological and pathological conditions. Both the overexpression and under-expression of Bim protein are associated with the diseased condition, and various isoforms of Bim protein are present with differential apoptotic potential. Objective: The present study attempted to envisage the association of various molecular signatures with the codon choices of Bim isoforms. Methods: Molecular signatures like composition, codon usage, nucleotide skews, the free energy of mRNA transcript, physical properties of proteins, codon adaptation index, relative synonymous codon usage, and dinucleotide odds ratio were determined and analyzed for their associations with codon choices of Bim gene. Results: Skew analysis of the Bim gene indicated the preference of C nucleotide over G, A, and T and preference of G over T and A nucleotides was observed. An increase in C content at the first and third codon position increased gene expression while it decreased at the second codon position. Compositional constraints on nucleotide C at all three codon positions affected gene expression. The analysis revealed an exceptionally high usage of CpC dinucleotide in all the envisaged 31 isoforms of Bim. We correlated it with the requirement of rapid demethylation machinery to fine-tune the Bimgene expression. Also, mutational pressure played a dominant role in shaping codon usage bias in Bim isoforms. Conclusion: An exceptionally high usage of CpC dinucleotide in all the envisaged 31 isoforms of Bim indicates a high order selectional force to fine tune Bim gene expression.
Collapse
Affiliation(s)
- Rekha Khandia
- Department of Biochemistry and Genetics, Barkatullah University, Bhopal, India
| | - Azmat Ali Khan
- Pharmaceutical Biotechnology Laboratory, Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, Riyadh, Saudi Arabia
| | - Athanasios Alexiou
- Novel Global Community Educational Foundation, Australia & AFNP Med, Austria
| | | | | |
Collapse
|
21
|
Auboeuf D. The Physics-Biology continuum challenges darwinism: Evolution is directed by the homeostasis-dependent bidirectional relation between genome and phenotype. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2021; 167:121-139. [PMID: 34097984 DOI: 10.1016/j.pbiomolbio.2021.05.008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/19/2021] [Revised: 05/19/2021] [Accepted: 05/31/2021] [Indexed: 10/21/2022]
Abstract
The physics-biology continuum relies on the fact that life emerged from prebiotic molecules. Here, I argue that life emerged from the coupling between nucleic acid and protein synthesis during which proteins (or proto-phenotypes) maintained the physicochemical parameter equilibria (or proto-homeostasis) in the proximity of their encoding nucleic acids (or proto-genomes). This protected the proto-genome physicochemical integrity (i.e., atomic composition) from environmental physicochemical constraints, and therefore increased the probability of reproducing the proto-genome without variation. From there, genomes evolved depending on the biological activities they generated in response to environmental fluctuations. Thus, a genome maintaining homeostasis (i.e., internal physicochemical parameter equilibria), despite and in response to environmental fluctuations, maintains its physicochemical integrity and has therefore a higher probability to be reproduced without variation. Consequently, descendants have a higher probability to share the same phenotype than their parents. Otherwise, the genome is modified during replication as a consequence of the imbalance of the internal physicochemical parameters it generates, until new mutation-deriving biological activities maintain homeostasis in offspring. In summary, evolution depends on feedforward and feedback loops between genome and phenotype, as the internal physicochemical conditions that a genome generates ─ through its derived phenotype in response to environmental fluctuations ─ in turn either guarantee its stability or direct its variation. Evolution may not be explained by the Darwinism-derived, unidirectional principle (random mutations-phenotypes-natural selection) but rather by the bidirectional relationship between genome and phenotype, in which the phenotype in interaction with the environment directs the evolution of the genome it derives from.
Collapse
Affiliation(s)
- Didier Auboeuf
- ENS de Lyon, Univ Claude Bernard, CNRS UMR 5239, INSERM U1210, Laboratory of Biology and Modelling of the Cell, 46 Allée D'Italie, Site Jacques Monod, F-69007, Lyon, France.
| |
Collapse
|
22
|
Analysis of synonymous codon usage bias in human monocytes, B, and T lymphocytes based on transcriptome data. GENE REPORTS 2021. [DOI: 10.1016/j.genrep.2021.101034] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
23
|
Dyson CJ, Goodisman MAD. Gene Duplication in the Honeybee: Patterns of DNA Methylation, Gene Expression, and Genomic Environment. Mol Biol Evol 2021; 37:2322-2331. [PMID: 32243528 DOI: 10.1093/molbev/msaa088] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2019] [Accepted: 03/27/2020] [Indexed: 12/12/2022] Open
Abstract
Gene duplication serves a critical role in evolutionary adaptation by providing genetic raw material to the genome. The evolution of duplicated genes may be influenced by epigenetic processes such as DNA methylation, which affects gene function in some taxa. However, the manner in which DNA methylation affects duplicated genes is not well understood. We studied duplicated genes in the honeybee Apis mellifera, an insect with a highly sophisticated social structure, to investigate whether DNA methylation was associated with gene duplication and genic evolution. We found that levels of gene body methylation were significantly lower in duplicate genes than in single-copy genes, implicating a possible role of DNA methylation in postduplication gene maintenance. Additionally, we discovered associations of gene body methylation with the location, length, and time since divergence of paralogous genes. We also found that divergence in DNA methylation was associated with divergence in gene expression in paralogs, although the relationship was not completely consistent with a direct link between DNA methylation and gene expression. Overall, our results provide further insight into genic methylation and how its association with duplicate genes might facilitate evolutionary processes and adaptation.
Collapse
Affiliation(s)
- Carl J Dyson
- School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA
| | | |
Collapse
|
24
|
Mordvinov VA, Ershov NI, Zaparina OG, Pakharukova MY. Genomics and proteomics of the liver fluke Opisthorchis felineus. Vavilovskii Zhurnal Genet Selektsii 2021; 24:383-390. [PMID: 33659821 PMCID: PMC7716572 DOI: 10.18699/vj20.44-o] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022] Open
Abstract
The causative agent of opisthorchiasis, the liver fluke Opisthorchis felineus (Rivolta, 1884) is one of the helminths of humans and animals in Russia. Together with closely related species of trematodes O. viverrini (Poirier, 1886) and Clonorchis sinensis (Loos, 1907), O. felineus is a part of a triad of epidemiologically important trematodes in the family Opisthorchiidae. Adult O. felineus worms infest the hepatobiliary system of warm-blooded animals and might provoke the development of severe pathologies, including malignancy of bile duct epithelium. The high medical importance of O. felineus attracts the attention of researchers. This review briefly summarizes the data about O. felineus genomics and proteomics. The review provides a comparative analysis of the number of genes and sizes of nuclear genomes of a number of flatworms, the distribution of intron lengths, as well as results of synteny between the O. felineus, O. viverrini and C. sinensis genomes. Special attention is paid to a particular form of RNA processing known as trans-splicing, widely presented in the opisthorchiid genomes. We also provide the results of a comparative analysis of the xenobiotic metabolizing system between parasitic and free-living flatworms. Moreover, data on parasitic granulins, which are potential promoters of cholangiocyte neoplasia, are also presented. Data on the O. felineus genomics and proteomics provide first insights into the structural and functional organization of the genome of this parasitic flatworm with a complex life cycle as well as provide a significant contribution to our understanding of “host-parasite” interaction and evolution of this group of parasitic flatworms.
Collapse
Affiliation(s)
- V A Mordvinov
- Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia
| | - N I Ershov
- Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia
| | - O G Zaparina
- Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia
| | - M Y Pakharukova
- Institute of Cytology and Genetics of Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russia
| |
Collapse
|
25
|
Li L, Yang W, Shen Y, Xu X, Li J. The evolutionary analysis of complement component C5 and the gene co-expression network and putative interaction between C5a and C5a anaphylatoxin receptor (C5AR/CD88) in human and two Cyprinid fish. DEVELOPMENTAL AND COMPARATIVE IMMUNOLOGY 2021; 116:103958. [PMID: 33290783 DOI: 10.1016/j.dci.2020.103958] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/16/2020] [Revised: 12/02/2020] [Accepted: 12/02/2020] [Indexed: 06/12/2023]
Abstract
The complement system is a complex network of soluble and membrane-associated serum proteins that regulate immune response. Activation of the complement C5 generates C5a and C5b which generate chemoattractive effect on myeloid cells and initiate the membrane attack complex (MAC) assembly. However, the study of evolutionary process and systematic function of C5 are still limited. In this study, we performed an evolutionary analysis of C5. Phylogeny analysis indicated that C5 sequences underwent complete divergence in fish and non-fish vertebrate. It was found that codon usage bias improved and provided evolution evidence of C5 in species. Notably, the codon usage bias of grass carp was evolutionarily closer to the zebrafish genome compared with humans and stickleback. This suggested that the zebrafish cell line may provide an alternative environment for heterologous protein expression of grass carp. Sequence comparison showed a higher similarity between human and mouse, grass carp, and zebrafish. Moreover, selective pressure analysis revealed that the C5 genes in fish and non-fish vertebrates exhibited different evolutionary patterns. To study the function of C5, gene co-expression networks of human and zebrafish were built which revealed the complexity of C5 function networks in different species. The protein structure simulation of C5 indicated that grass carp and zebrafish are more similar than to human, however, differences between species in C5a proteins are extremely smaller. Spatial conformations of C5a-C5AR (CD88) protein complex were constructed, which showed that possible interaction may exist between C5a and CD88 proteins. Furthermore, the protein docking sites/residues were measured and calculated according to the minimum distance for all atoms from C5a and CD88 proteins. In summary, this study provides insights into the evolutionary history, function and potential regulatory mechanism of C5 in fish immune responses.
Collapse
Affiliation(s)
- Lisen Li
- Key Laboratory of Freshwater Aquatic Genetic Resources, Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, 201306, China
| | - Weining Yang
- Key Laboratory of Freshwater Aquatic Genetic Resources, Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, 201306, China
| | - Yubang Shen
- Key Laboratory of Freshwater Aquatic Genetic Resources, Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, 201306, China.
| | - Xiaoyan Xu
- Key Laboratory of Freshwater Aquatic Genetic Resources, Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, 201306, China
| | - Jiale Li
- Key Laboratory of Freshwater Aquatic Genetic Resources, Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, 201306, China; Shanghai Engineering Research Center of Aquaculture, Shanghai Ocean University, Shanghai, 201306, China.
| |
Collapse
|
26
|
Blumberg A, Zhao Y, Huang YF, Dukler N, Rice EJ, Chivu AG, Krumholz K, Danko CG, Siepel A. Characterizing RNA stability genome-wide through combined analysis of PRO-seq and RNA-seq data. BMC Biol 2021; 19:30. [PMID: 33588838 PMCID: PMC7885420 DOI: 10.1186/s12915-021-00949-x] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Accepted: 01/05/2021] [Indexed: 01/03/2023] Open
Abstract
BACKGROUND The concentrations of distinct types of RNA in cells result from a dynamic equilibrium between RNA synthesis and decay. Despite the critical importance of RNA decay rates, current approaches for measuring them are generally labor-intensive, limited in sensitivity, and/or disruptive to normal cellular processes. Here, we introduce a simple method for estimating relative RNA half-lives that is based on two standard and widely available high-throughput assays: Precision Run-On sequencing (PRO-seq) and RNA sequencing (RNA-seq). RESULTS Our method treats PRO-seq as a measure of transcription rate and RNA-seq as a measure of RNA concentration, and estimates the rate of RNA decay required for a steady-state equilibrium. We show that this approach can be used to assay relative RNA half-lives genome-wide, with good accuracy and sensitivity for both coding and noncoding transcription units. Using a structural equation model (SEM), we test several features of transcription units, nearby DNA sequences, and nearby epigenomic marks for associations with RNA stability after controlling for their effects on transcription. We find that RNA splicing-related features are positively correlated with RNA stability, whereas features related to miRNA binding and DNA methylation are negatively correlated with RNA stability. Furthermore, we find that a measure based on U1 binding and polyadenylation sites distinguishes between unstable noncoding and stable coding transcripts but is not predictive of relative stability within the mRNA or lincRNA classes. We also identify several histone modifications that are associated with RNA stability. CONCLUSION We introduce an approach for estimating the relative half-lives of individual RNAs. Together, our estimation method and systematic analysis shed light on the pervasive impacts of RNA stability on cellular RNA concentrations.
Collapse
Affiliation(s)
- Amit Blumberg
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | - Yixin Zhao
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | - Yi-Fei Huang
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
- Present Address: Department of Biology and Huck Institutes of the Life Sciences, Pennsylvania State University, University Park, PA, USA
| | - Noah Dukler
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | - Edward J Rice
- Baker Institute for Animal Health, College of Veterinary Medicine, Cornell University, Ithaca, NY, USA
| | - Alexandra G Chivu
- Baker Institute for Animal Health, College of Veterinary Medicine, Cornell University, Ithaca, NY, USA
| | - Katie Krumholz
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA
| | - Charles G Danko
- Baker Institute for Animal Health, College of Veterinary Medicine, Cornell University, Ithaca, NY, USA
| | - Adam Siepel
- Simons Center for Quantitative Biology, Cold Spring Harbor Laboratory, Cold Spring Harbor, NY, USA.
| |
Collapse
|
27
|
Role of Gene Length in Control of Human Gene Expression: Chromosome-Specific and Tissue-Specific Effects. Int J Genomics 2021; 2021:8902428. [PMID: 33688492 PMCID: PMC7911607 DOI: 10.1155/2021/8902428] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2020] [Revised: 01/12/2021] [Accepted: 02/03/2021] [Indexed: 11/19/2022] Open
Abstract
This study was carried out to pursue the observation that the level of gene expression is affected by gene length in the human genome. As transcription is a time-dependent process, it is expected that gene expression will be inversely related to gene length, and this is found to be the case. Here, I describe the results of studies performed to test whether the gene length/gene expression linkage is affected by two factors, the chromosome where the gene is located and the tissue where it is expressed. Studies were performed with a database of 3538 human genes that were divided into short, midlength, and long groups. Chromosome groups were then compared in the expression level of genes with the same length. A similar analysis was performed with 19 human tissues. Tissue-specific groups were compared in the expression level of genes with the same length. Both chromosome and tissue studies revealed new information about the role of gene length in control of gene expression. Chromosome studies led to the identification of two chromosome populations that differ in the expression level of short genes. A high level of expression was observed in chromosomes 2-10, 12-15, and 18 and a low level in 1, 11, 16-17, 19-20, 22, and 24. Studies with tissue-specific genes led to the identification of two tissues, brain and liver, which differ in the expression level of short genes. The results are interpreted to support the view that the level of a gene's expression can be affected by the chromosome and the tissue where the gene is transcribed.
Collapse
|
28
|
Lopes I, Altab G, Raina P, de Magalhães JP. Gene Size Matters: An Analysis of Gene Length in the Human Genome. Front Genet 2021; 12:559998. [PMID: 33643374 PMCID: PMC7905317 DOI: 10.3389/fgene.2021.559998] [Citation(s) in RCA: 61] [Impact Index Per Article: 20.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2020] [Accepted: 01/06/2021] [Indexed: 12/23/2022] Open
Abstract
While it is expected for gene length to be associated with factors such as intron number and evolutionary conservation, we are yet to understand the connections between gene length and function in the human genome. In this study, we show that, as expected, there is a strong positive correlation between gene length, transcript length, and protein size as well as a correlation with the number of genetic variants and introns. Among tissue-specific genes, we find that the longest transcripts tend to be expressed in the blood vessels, nerves, thyroid, cervix uteri, and the brain, while the smallest transcripts tend to be expressed in the pancreas, skin, stomach, vagina, and testis. We report, as shown previously, that natural selection suppresses changes for genes with longer transcripts and promotes changes for genes with smaller transcripts. We also observe that genes with longer transcripts tend to have a higher number of co-expressed genes and protein-protein interactions, as well as more associated publications. In the functional analysis, we show that bigger transcripts are often associated with neuronal development, while smaller transcripts tend to play roles in skin development and in the immune system. Furthermore, pathways related to cancer, neurons, and heart diseases tend to have genes with longer transcripts, with smaller transcripts being present in pathways related to immune responses and neurodegenerative diseases. Based on our results, we hypothesize that longer genes tend to be associated with functions that are important in the early development stages, while smaller genes tend to play a role in functions that are important throughout the whole life, like the immune system, which requires fast responses.
Collapse
Affiliation(s)
| | | | | | - João Pedro de Magalhães
- Integrative Genomics of Ageing Group, Institute of Ageing and Chronic Disease, University of Liverpool, Liverpool, United Kingdom
| |
Collapse
|
29
|
Chamani Mohasses F, Solouki M, Ghareyazie B, Fahmideh L, Mohsenpour M. Correlation between gene expression levels under drought stress and synonymous codon usage in rice plant by in-silico study. PLoS One 2020; 15:e0237334. [PMID: 32776991 PMCID: PMC7416939 DOI: 10.1371/journal.pone.0237334] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2020] [Accepted: 07/23/2020] [Indexed: 11/24/2022] Open
Abstract
We studied the correlation of synonymous codon usage (SCU) on gene expression levels under drought stress in rice. Sixty genes related to drought stress (with high, intermediate and low expression) were selected from rice meta-analysis data and various codon usage indices such as the effective number of codon usage (ENC), codon adaptation index (CAI) and relative synonymous codon usage (RSCU) were calculated. We found that in genes highly expressing under drought 1) GC content was higher, 2) ENC value was lower, 3) the preferred codons of some amino acids changed and 4) the RSCU ratio of GC-end codons relative to AT-end codons for 18 amino acids increased significantly compared with those in other genes. We introduce ARSCU as the Average ratio of RSCUs of GC-end codons to AT-end codons in each gene that could significantly separate high-expression genes under drought from low-expression genes. ARSCU is calculated using the program ARSCU-Calculator developed by our group to help predicting expression level of rice genes under drought. An index above ARSCU threshold is expected to indicate that the gene under study may belong to the "high expression group under drought". This information may be applied for codon optimization of genes for rice genetic engineering. To validate these findings, we further used 60 other genes (randomly selected subset of 43233 genes studied for their response to drought stress). ARSCU value was able to predict the level of expression at 88.33% of the cases. Using third set of 60 genes selected amongst high expressing genes not related to drought, only 31.65% of the genes showed ARSCU value of higher than the set threshold. This indicates that the phenomenon we described in this report may be unique for drought related genes. To justify the observed correlation between CUB and high expressing genes under drought, possible role of tRNA post transcriptional modification and tRFs was hypothesized as possible underlying biological mechanism.
Collapse
Affiliation(s)
- Fatemeh Chamani Mohasses
- Department of Plant Breeding and Biotechnology (PBB), Faculty of Agriculture, University of Zabol, Zabol, Iran
| | - Mahmood Solouki
- Department of Plant Breeding and Biotechnology (PBB), Faculty of Agriculture, University of Zabol, Zabol, Iran
| | - Behzad Ghareyazie
- Agricultural Biotechnology Research Institute of Iran (ABRII), Agricultural Research Education and Extension Organization (AREEO), Karaj, Iran
| | - Leila Fahmideh
- Department of Plant Breeding and Biotechnology (PBB), Faculty of Agriculture, University of Zabol, Zabol, Iran
| | - Motahhareh Mohsenpour
- Agricultural Biotechnology Research Institute of Iran (ABRII), Agricultural Research Education and Extension Organization (AREEO), Karaj, Iran
| |
Collapse
|
30
|
Abstract
Darwin's theory of evolution emphasized that positive selection of functional proficiency provides the fitness that ultimately determines the structure of life, a view that has dominated biochemical thinking of enzymes as perfectly optimized for their specific functions. The 20th-century modern synthesis, structural biology, and the central dogma explained the machinery of evolution, and nearly neutral theory explained how selection competes with random fixation dynamics that produce molecular clocks essential e.g. for dating evolutionary histories. However, quantitative proteomics revealed that selection pressures not relating to optimal function play much larger roles than previously thought, acting perhaps most importantly via protein expression levels. This paper first summarizes recent progress in the 21st century toward recovering this universal selection pressure. Then, the paper argues that proteome cost minimization is the dominant, underlying 'non-function' selection pressure controlling most of the evolution of already functionally adapted living systems. A theory of proteome cost minimization is described and argued to have consequences for understanding evolutionary trade-offs, aging, cancer, and neurodegenerative protein-misfolding diseases.
Collapse
|
31
|
Auboeuf D. Physicochemical Foundations of Life that Direct Evolution: Chance and Natural Selection are not Evolutionary Driving Forces. Life (Basel) 2020; 10:life10020007. [PMID: 31973071 PMCID: PMC7175370 DOI: 10.3390/life10020007] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2019] [Revised: 01/15/2020] [Accepted: 01/16/2020] [Indexed: 12/11/2022] Open
Abstract
The current framework of evolutionary theory postulates that evolution relies on random mutations generating a diversity of phenotypes on which natural selection acts. This framework was established using a top-down approach as it originated from Darwinism, which is based on observations made of complex multicellular organisms and, then, modified to fit a DNA-centric view. In this article, it is argued that based on a bottom-up approach starting from the physicochemical properties of nucleic and amino acid polymers, we should reject the facts that (i) natural selection plays a dominant role in evolution and (ii) the probability of mutations is independent of the generated phenotype. It is shown that the adaptation of a phenotype to an environment does not correspond to organism fitness, but rather corresponds to maintaining the genome stability and integrity. In a stable environment, the phenotype maintains the stability of its originating genome and both (genome and phenotype) are reproduced identically. In an unstable environment (i.e., corresponding to variations in physicochemical parameters above a physiological range), the phenotype no longer maintains the stability of its originating genome, but instead influences its variations. Indeed, environment- and cellular-dependent physicochemical parameters define the probability of mutations in terms of frequency, nature, and location in a genome. Evolution is non-deterministic because it relies on probabilistic physicochemical rules, and evolution is driven by a bidirectional interplay between genome and phenotype in which the phenotype ensures the stability of its originating genome in a cellular and environmental physicochemical parameter-depending manner.
Collapse
Affiliation(s)
- Didier Auboeuf
- Laboratory of Biology and Modelling of the Cell, Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS UMR 5239, INSERM U1210, 46 Allée d'Italie, Site Jacques Monod, F-69007, Lyon, France
| |
Collapse
|
32
|
Evolutionary Forces and Codon Bias in Different Flavors of Intrinsic Disorder in the Human Proteome. J Mol Evol 2019; 88:164-178. [DOI: 10.1007/s00239-019-09921-4] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2019] [Accepted: 11/26/2019] [Indexed: 12/22/2022]
|
33
|
Lemaire S, Fontrodona N, Aubé F, Claude JB, Polvèche H, Modolo L, Bourgeois CF, Mortreux F, Auboeuf D. Characterizing the interplay between gene nucleotide composition bias and splicing. Genome Biol 2019; 20:259. [PMID: 31783898 PMCID: PMC6883713 DOI: 10.1186/s13059-019-1869-y] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2019] [Accepted: 10/28/2019] [Indexed: 12/18/2022] Open
Abstract
BACKGROUND Nucleotide composition bias plays an important role in the 1D and 3D organization of the human genome. Here, we investigate the potential interplay between nucleotide composition bias and the regulation of exon recognition during splicing. RESULTS By analyzing dozens of RNA-seq datasets, we identify two groups of splicing factors that activate either about 3200 GC-rich exons or about 4000 AT-rich exons. We show that splicing factor-dependent GC-rich exons have predicted RNA secondary structures at 5' ss and are dependent on U1 snRNP-associated proteins. In contrast, splicing factor-dependent AT-rich exons have a large number of decoy branch points, SF1- or U2AF2-binding sites and are dependent on U2 snRNP-associated proteins. Nucleotide composition bias also influences local chromatin organization, with consequences for exon recognition during splicing. Interestingly, the GC content of exons correlates with that of their hosting genes, isochores, and topologically associated domains. CONCLUSIONS We propose that regional nucleotide composition bias over several dozens of kilobase pairs leaves a local footprint at the exon level and induces constraints during splicing that can be alleviated by local chromatin organization at the DNA level and recruitment of specific splicing factors at the RNA level. Therefore, nucleotide composition bias establishes a direct link between genome organization and local regulatory processes, like alternative splicing.
Collapse
Affiliation(s)
- Sébastien Lemaire
- Laboratory of Biology and Modelling of the Cell, Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS UMR 5239, INSERM U1210, 46 Allée d'Italie Site Jacques Monod, F-69007, Lyon, France
| | - Nicolas Fontrodona
- Laboratory of Biology and Modelling of the Cell, Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS UMR 5239, INSERM U1210, 46 Allée d'Italie Site Jacques Monod, F-69007, Lyon, France
| | - Fabien Aubé
- Laboratory of Biology and Modelling of the Cell, Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS UMR 5239, INSERM U1210, 46 Allée d'Italie Site Jacques Monod, F-69007, Lyon, France
| | - Jean-Baptiste Claude
- Laboratory of Biology and Modelling of the Cell, Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS UMR 5239, INSERM U1210, 46 Allée d'Italie Site Jacques Monod, F-69007, Lyon, France
| | | | - Laurent Modolo
- LBMC Biocomputing Center, CNRS UMR 5239, INSERM U1210, 46 Allée d'Italie Site Jacques Monod, F-69007, Lyon, France
| | - Cyril F Bourgeois
- Laboratory of Biology and Modelling of the Cell, Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS UMR 5239, INSERM U1210, 46 Allée d'Italie Site Jacques Monod, F-69007, Lyon, France
| | - Franck Mortreux
- Laboratory of Biology and Modelling of the Cell, Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS UMR 5239, INSERM U1210, 46 Allée d'Italie Site Jacques Monod, F-69007, Lyon, France
| | - Didier Auboeuf
- Laboratory of Biology and Modelling of the Cell, Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS UMR 5239, INSERM U1210, 46 Allée d'Italie Site Jacques Monod, F-69007, Lyon, France.
| |
Collapse
|
34
|
Hunt R, Hettiarachchi G, Katneni U, Hernandez N, Holcomb D, Kames J, Alnifaidy R, Lin B, Hamasaki-Katagiri N, Wesley A, Kafri T, Morris C, Bouché L, Panico M, Schiller T, Ibla J, Bar H, Ismail A, Morris H, Komar A, Kimchi-Sarfaty C. A Single Synonymous Variant (c.354G>A [p.P118P]) in ADAMTS13 Confers Enhanced Specific Activity. Int J Mol Sci 2019; 20:ijms20225734. [PMID: 31731663 PMCID: PMC6888508 DOI: 10.3390/ijms20225734] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2019] [Revised: 11/09/2019] [Accepted: 11/13/2019] [Indexed: 12/18/2022] Open
Abstract
Synonymous variants within coding regions may influence protein expression and function. We have previously reported increased protein expression levels ex vivo (~120% in comparison to wild-type) from a synonymous polymorphism variant, c.354G>A [p.P118P], of the ADAMTS13 gene, encoding a plasma protease responsible for von Willebrand Factor (VWF) degradation. In the current study, we investigated the potential mechanism(s) behind the increased protein expression levels from this variant and its effect on ADAMTS13 physico-chemical properties. Cell-free assays showed enhanced translation of the c.354G>A variant and the analysis of codon usage characteristics suggested that introduction of the frequently used codon/codon pair(s) may have been potentially responsible for this effect. Limited proteolysis, however, showed no substantial influence of altered translation on protein conformation. Analysis of post-translational modifications also showed no notable differences but identified three previously unreported glycosylation markers. Despite these similarities, p.P118P variant unexpectedly showed higher specific activity. Structural analysis using modeled interactions indicated that subtle conformational changes arising from altered translation kinetics could affect interactions between an exosite of ADAMTS13 and VWF resulting in altered specific activity. This report highlights how a single synonymous nucleotide variation can impact cellular expression and specific activity in the absence of measurable impact on protein structure.
Collapse
Affiliation(s)
- Ryan Hunt
- Hemostasis Branch, Division of Plasma Protein Therapeutics, Office of Tissues and Advanced Therapies, Center for Biologics Evaluation & Research, US FDA, Silver Spring, MD 20993, USA
| | - Gaya Hettiarachchi
- Hemostasis Branch, Division of Plasma Protein Therapeutics, Office of Tissues and Advanced Therapies, Center for Biologics Evaluation & Research, US FDA, Silver Spring, MD 20993, USA
| | - Upendra Katneni
- Hemostasis Branch, Division of Plasma Protein Therapeutics, Office of Tissues and Advanced Therapies, Center for Biologics Evaluation & Research, US FDA, Silver Spring, MD 20993, USA
| | - Nancy Hernandez
- Hemostasis Branch, Division of Plasma Protein Therapeutics, Office of Tissues and Advanced Therapies, Center for Biologics Evaluation & Research, US FDA, Silver Spring, MD 20993, USA
| | - David Holcomb
- Hemostasis Branch, Division of Plasma Protein Therapeutics, Office of Tissues and Advanced Therapies, Center for Biologics Evaluation & Research, US FDA, Silver Spring, MD 20993, USA
| | - Jacob Kames
- Hemostasis Branch, Division of Plasma Protein Therapeutics, Office of Tissues and Advanced Therapies, Center for Biologics Evaluation & Research, US FDA, Silver Spring, MD 20993, USA
| | - Redab Alnifaidy
- Hemostasis Branch, Division of Plasma Protein Therapeutics, Office of Tissues and Advanced Therapies, Center for Biologics Evaluation & Research, US FDA, Silver Spring, MD 20993, USA
| | - Brian Lin
- Hemostasis Branch, Division of Plasma Protein Therapeutics, Office of Tissues and Advanced Therapies, Center for Biologics Evaluation & Research, US FDA, Silver Spring, MD 20993, USA
| | - Nobuko Hamasaki-Katagiri
- Hemostasis Branch, Division of Plasma Protein Therapeutics, Office of Tissues and Advanced Therapies, Center for Biologics Evaluation & Research, US FDA, Silver Spring, MD 20993, USA
| | - Aaron Wesley
- Gene Therapy Center, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
- Present Address: Department of Emergency Medicine, Banner University Medical Center, The University of Arizona, Tucson, AZ 85724, USA
| | - Tal Kafri
- Gene Therapy Center, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | | | - Laura Bouché
- Department of Life Sciences, Imperial College London, South Kensington Campus, London SW7 2AZ, UK
- Present Address: Antikor Biopharma Ltd., Stevenage Bioscience Catalyst, Gunnels Wood Road, Stevenage SG1 2FX, UK
| | - Maria Panico
- BioPharmaSpec Ltd., St. Saviour JE2 7LA, UK or or
- Department of Life Sciences, Imperial College London, South Kensington Campus, London SW7 2AZ, UK
| | - Tal Schiller
- Hemostasis Branch, Division of Plasma Protein Therapeutics, Office of Tissues and Advanced Therapies, Center for Biologics Evaluation & Research, US FDA, Silver Spring, MD 20993, USA
| | - Juan Ibla
- Departments of Cardiac Surgery and Anesthesiology, Perioperative and Pain Medicine, Boston Children’s Hospital and Harvard Medical School, Boston, MA 02115, USA
| | - Haim Bar
- Department of Statistics, University of Connecticut, Storrs, CT 06269, USA
| | - Amra Ismail
- Center for Gene Regulation in Health and Disease, Department of Biological, Geological & Environmental Sciences, Cleveland State University, Cleveland, OH 44115, USA
| | - Howard Morris
- BioPharmaSpec Ltd., St. Saviour JE2 7LA, UK or or
- Department of Life Sciences, Imperial College London, South Kensington Campus, London SW7 2AZ, UK
| | - Anton Komar
- Center for Gene Regulation in Health and Disease, Department of Biological, Geological & Environmental Sciences, Cleveland State University, Cleveland, OH 44115, USA
| | - Chava Kimchi-Sarfaty
- Hemostasis Branch, Division of Plasma Protein Therapeutics, Office of Tissues and Advanced Therapies, Center for Biologics Evaluation & Research, US FDA, Silver Spring, MD 20993, USA
- Correspondence: ; Tel.: +1-(240)-402-8203
| |
Collapse
|
35
|
Cao Y, Jiang L, Wang L, Cai Y. Evolutionary Rate Heterogeneity and Functional Divergence of Orthologous Genes in Pyrus. Biomolecules 2019; 9:biom9090490. [PMID: 31527450 PMCID: PMC6770726 DOI: 10.3390/biom9090490] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2019] [Revised: 09/09/2019] [Accepted: 09/12/2019] [Indexed: 11/21/2022] Open
Abstract
Negatively selected genes (NSGs) and positively selected genes (PSGs) are the two types of most nuclear protein-coding genes in organisms. However, the evolutionary rates and characteristics of different types of genes have been rarely understood. In the present study, we investigate the rates of synonymous substitution (Ks) and the rates of non-synonymous substitution (Ka) by comparing the orthologous genes of two sequenced Pyrus species, Pyrus bretschneideri and Pyrus communis. Subsequently, we compared the evolutionary rates, gene structures, and expression profiles during different fruit development between PSGs and NSGs. Compared with the NSGs, the PSGs have fewer exons, shorter gene length, lower synonymous substitution rates and have higher evolutionary rates. Remarkably, gene expression patterns between two Pyrus species fruit indicated functional divergence for most of the orthologous genes derived from a common ancestor, and subfunctionalization for some of them. Overall, the present study shows that PSGs differs from NSGs not only under environmental selective pressure (Ka/Ks), but also in their structural, functional, and evolutionary properties. Additionally, our resulting data provides important insights for the evolution and highlights the diversification of orthologous genes in two Pyrus species.
Collapse
Affiliation(s)
- Yunpeng Cao
- Key Laboratory of Cultivation and Protection for Non-Wood Forest Trees, Ministry of Education, Central South University of Forestry and Technology, Changsha 410004, China.
- School of Life Sciences, Anhui Agricultural University, Hefei 230036, China.
| | - Lan Jiang
- Key Laboratory of Cultivation and Protection for Non-Wood Forest Trees, Ministry of Education, Central South University of Forestry and Technology, Changsha 410004, China.
| | - Lihu Wang
- College of Landscape and Ecological Engineering, Hebei University of Engineering, Handan 056038, China.
| | - Yongping Cai
- School of Life Sciences, Anhui Agricultural University, Hefei 230036, China.
| |
Collapse
|
36
|
Molecular cloning and analysis of the full-length aciniform spidroin gene from Araneus ventricosus. Int J Biol Macromol 2018; 117:1352-1360. [DOI: 10.1016/j.ijbiomac.2017.12.090] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2017] [Revised: 12/16/2017] [Accepted: 12/16/2017] [Indexed: 11/19/2022]
|
37
|
Wei X, Zhang J. On the Origin of Compositional Features of Ribosomes. Genome Biol Evol 2018; 10:2010-2016. [PMID: 30059996 PMCID: PMC6097593 DOI: 10.1093/gbe/evy169] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/28/2018] [Indexed: 01/08/2023] Open
Abstract
Ribosomes are highly abundant in cells and comprise, besides RNAs of varying lengths, 55–80 similarly sized, short proteins. This seemingly unusual composition is thought to have resulted from selection for rapid autocatalytic ribosome production. Here, we demonstrate that ribosomal protein-splitting mutations cannot accelerate ribosome production. The autocatalytic explanation is also unnecessary, because protein lengths generally decline with expression levels. Although ribosomal proteins are shorter than expected from their expression levels, they are not outliers among members of large protein complexes in mean protein length or coefficient of variation. These observations are explainable because 1) shortening proteins lowers their synthetic cost and reduces the waste from mistranslation-induced protein dysfunction and degradation, 2) such benefits rise with expression levels, and 3) members of large complexes participate in more protein–protein interactions so are less tolerant to mistranslation. These and other considerations suggest that the compositional features of ribosomes originate from cellular energy economics.
Collapse
Affiliation(s)
- Xinzhu Wei
- 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
|
38
|
Abrahams L, Hurst LD. Adenine Enrichment at the Fourth CDS Residue in Bacterial Genes Is Consistent with Error Proofing for +1 Frameshifts. Mol Biol Evol 2018; 34:3064-3080. [PMID: 28961919 PMCID: PMC5850271 DOI: 10.1093/molbev/msx223] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Beyond selection for optimal protein functioning, coding sequences (CDSs) are under selection at the RNA and DNA levels. Here, we identify a possible signature of “dual-coding,” namely extensive adenine (A) enrichment at bacterial CDS fourth sites. In 99.07% of studied bacterial genomes, fourth site A use is greater than expected given genomic A-starting codon use. Arguing for nucleotide level selection, A-starting serine and arginine second codons are heavily utilized when compared with their non-A starting synonyms. Several models have the ability to explain some of this trend. In part, A-enrichment likely reduces 5′ mRNA stability, promoting translation initiation. However T/U, which may also reduce stability, is avoided. Further, +1 frameshifts on the initiating ATG encode a stop codon (TGA) provided A is the fourth residue, acting either as a frameshift “catch and destroy” or a frameshift stop and adjust mechanism and hence implicated in translation initiation. Consistent with both, genomes lacking TGA stop codons exhibit weaker fourth site A-enrichment. Sequences lacking a Shine–Dalgarno sequence and those without upstream leader genes, that may be more error prone during initiation, have greater utilization of A, again suggesting a role in initiation. The frameshift correction model is consistent with the notion that many genomic features are error-mitigation factors and provides the first evidence for site-specific out of frame stop codon selection. We conjecture that the NTG universal start codon may have evolved as a consequence of TGA being a stop codon and the ability of NTGA to rapidly terminate or adjust a ribosome.
Collapse
Affiliation(s)
- Liam Abrahams
- Department of Biology and Biochemistry, The Milner Centre for Evolution, University of Bath, Bath, United Kingdom
| | - Laurence D Hurst
- Department of Biology and Biochemistry, The Milner Centre for Evolution, University of Bath, Bath, United Kingdom
| |
Collapse
|
39
|
Panda A, Acharya D, Chandra Ghosh T. Insights into human intrinsically disordered proteins from their gene expression profile. MOLECULAR BIOSYSTEMS 2018; 13:2521-2530. [PMID: 29051952 DOI: 10.1039/c7mb00311k] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Expression level provides important clues about gene function. Previously, various efforts have been undertaken to profile human genes according to their expression level. Intrinsically disordered proteins (IDPs) do not adopt any rigid conformation under physiological conditions, however, are considered as an important functional class in all domains of life. Based on a human tissue-averaged gene expression level, previous studies showed that IDPs are expressed at a lower level than ordered globular proteins. Here, we examined the gene expression pattern of human ordered and disordered proteins in 32 normal tissues. We noticed that in most of the tissues, ordered and disordered proteins are expressed at a similar level. Moreover, in a number of tissues IDPs were found to be expressed at a higher level than ordered proteins. Rigorous statistical analyses suggested that the lower tissue-averaged gene expression level of IDPs (reported earlier) may be the consequence of their biased gene expression in some specific tissues and higher protein length. When we considered the gene repertory of each tissue we noticed that a number of human tissues (brain, testes, etc.) selectively express a higher fraction of disordered proteins, which help them to maintain higher protein connectivity by forming disordered binding motifs and to sustain their functional specificities. Our results demonstrated that the disordered proteins are indispensable in these tissues for their functional advantages.
Collapse
Affiliation(s)
- Arup Panda
- Bioinformatics Centre, Bose Institute, P 1/12, C.I.T. Scheme VII M, Kolkata 700 054, West Bengal, India.
| | | | | |
Collapse
|
40
|
Wen R, Liu X, Meng Q. Characterization of full-length tubuliform spidroin gene from Araneus ventricosus. Int J Biol Macromol 2017; 105:702-710. [DOI: 10.1016/j.ijbiomac.2017.07.086] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2017] [Revised: 06/12/2017] [Accepted: 07/13/2017] [Indexed: 10/19/2022]
|
41
|
Goswami AM. Codon usage patterns of 3β-hydroxysteroid dehydrogenase type 2 gene across mammalian species and the influence of mutation and selection pressure. GENE REPORTS 2017. [DOI: 10.1016/j.genrep.2017.08.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
42
|
Analysis of the codon usage pattern in Middle East Respiratory Syndrome Coronavirus. Oncotarget 2017; 8:110337-110349. [PMID: 29299151 PMCID: PMC5746386 DOI: 10.18632/oncotarget.22738] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2017] [Accepted: 10/27/2017] [Indexed: 11/25/2022] Open
Abstract
Middle East Respiratory Syndrome Coronavirus (MERS-CoV), which first broken out in Jeddah in 2012, causes a severe acute respiratory illness with a high mortality rate. To better understand the molecular characteristics of isolated MERS-CoV genomes, we first analysed the codon usage pattern of the zoonotic MERS-CoV strains comprehensively to gain an insight into the mechanism of cross-species transmission. We found that MERS human/camel isolates showed a low codon usage bias. Both mutation and nature selection pressure have contributed to this low codon usage bias, with the former being the main determining factor. We also observed that gene function, evolution time and the different host species of the virus all contributed to the bias of MERS-CoV, to some extent. Additionally, the codon usage pattern of MERS-CoV isolates is different from other related Nidovirales viruses isolated from bats and hedgehogs. In the future, more epidemiological surveys are required to examine the factors that resulted in the emergence and outbreak of this virus.
Collapse
|
43
|
Yi S, Li Y, Wang W. Selection shapes the patterns of codon usage in three closely related species of genus Misgurnus. Genomics 2017; 110:134-142. [PMID: 28911975 DOI: 10.1016/j.ygeno.2017.09.004] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2017] [Revised: 08/29/2017] [Accepted: 09/07/2017] [Indexed: 11/17/2022]
Abstract
Neutrality plots revealed that selection probably dominates codon bias, whereas mutation plays only a minor role, in shaping the codon bias in three loaches, Misgurnus anguillicaudatus, M. mohoity, and M. bipartitus. These three species also clearly showed similar tendencies in the preferential usage of codons. Nineteen, nine, and 14 preferred codon pairs and 179, 182, and 173 avoided codon pairs were also detected in M. anguillicaudatus, M. bipartitus, and M. mohoity, respectively, and the most frequently avoided type of cP3-cA1 dinucleotide in these species was nnUAnn. The expression-linked patterns of codon usage revealed that higher expression was associated with higher GC3, lower ENC, and a smaller proportion of amino acids with high size/complexity (S/C) scores in these three species. These results elucidate selectively driven codon bias in Misgurnus species, and reveal the potential importance of expression-mediated selection in shaping the genome evolution of fish.
Collapse
Affiliation(s)
- Shaokui Yi
- College of Fisheries, Key Lab of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education/Key Lab of Freshwater Animal Breeding, Ministry of Agriculture, Huazhong Agricultural University, Wuhan 430070, China
| | - Yanhe Li
- College of Fisheries, Key Lab of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education/Key Lab of Freshwater Animal Breeding, Ministry of Agriculture, Huazhong Agricultural University, Wuhan 430070, China.
| | - Weimin Wang
- College of Fisheries, Key Lab of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education/Key Lab of Freshwater Animal Breeding, Ministry of Agriculture, Huazhong Agricultural University, Wuhan 430070, China.
| |
Collapse
|
44
|
Guo Y, Liu J, Zhang J, Liu S, Du J. Selective modes determine evolutionary rates, gene compactness and expression patterns in Brassica. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2017; 91:34-44. [PMID: 28332757 DOI: 10.1111/tpj.13541] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/28/2016] [Revised: 02/28/2017] [Accepted: 03/15/2017] [Indexed: 05/18/2023]
Abstract
It has been well documented that most nuclear protein-coding genes in organisms can be classified into two categories: positively selected genes (PSGs) and negatively selected genes (NSGs). The characteristics and evolutionary fates of different types of genes, however, have been poorly understood. In this study, the rates of nonsynonymous substitution (Ka ) and the rates of synonymous substitution (Ks ) were investigated by comparing the orthologs between the two sequenced Brassica species, Brassica rapa and Brassica oleracea, and the evolutionary rates, gene structures, expression patterns, and codon bias were compared between PSGs and NSGs. The resulting data show that PSGs have higher protein evolutionary rates, lower synonymous substitution rates, shorter gene length, fewer exons, higher functional specificity, lower expression level, higher tissue-specific expression and stronger codon bias than NSGs. Although the quantities and values are different, the relative features of PSGs and NSGs have been largely verified in the model species Arabidopsis. These data suggest that PSGs and NSGs differ not only under selective pressure (Ka /Ks ), but also in their evolutionary, structural and functional properties, indicating that selective modes may serve as a determinant factor for measuring evolutionary rates, gene compactness and expression patterns in Brassica.
Collapse
Affiliation(s)
- Yue Guo
- Provincial Key Laboratory of Agrobiology, Institute of Biotechnology, Jiangsu Academy of Agricultural Sciences, Nanjing, 210014, China
| | - Jing Liu
- Provincial Key Laboratory of Agrobiology, Institute of Biotechnology, Jiangsu Academy of Agricultural Sciences, Nanjing, 210014, China
| | - Jiefu Zhang
- Key Laboratory of Cotton and Rapeseed, Ministry of Agriculture of People's Republic of China, Institute of Industrial Crops, Jiangsu Academy of Agricultural Sciences, Nanjing, 210014, China
| | - Shengyi Liu
- Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture of People's Republic of China, Oil Crops Research Institute, Chinese Academy of Agricultural Sciences, Wuhan, 430062, China
| | - Jianchang Du
- Provincial Key Laboratory of Agrobiology, Institute of Biotechnology, Jiangsu Academy of Agricultural Sciences, Nanjing, 210014, China
- Key Laboratory of Cotton and Rapeseed, Ministry of Agriculture of People's Republic of China, Institute of Industrial Crops, Jiangsu Academy of Agricultural Sciences, Nanjing, 210014, China
- Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture of People's Republic of China, Oil Crops Research Institute, Chinese Academy of Agricultural Sciences, Wuhan, 430062, China
| |
Collapse
|
45
|
Abstract
Males and females exhibit highly dimorphic phenotypes, particularly in their gonads, which is believed to be driven largely by differential gene expression. Typically, the protein sequences of genes upregulated in males, or male-biased genes, evolve rapidly as compared to female-biased and unbiased genes. To date, the specific study of gonad-biased genes remains uncommon in metazoans. Here, we identified and studied a total of 2927, 2013, and 4449 coding sequences (CDS) with ovary-biased, testis-biased, and unbiased expression, respectively, in the yellow fever mosquito Aedes aegypti The results showed that ovary-biased and unbiased CDS had higher nonsynonymous to synonymous substitution rates (dN/dS) and lower optimal codon usage (those codons that promote efficient translation) than testis-biased genes. Further, we observed higher dN/dS in ovary-biased genes than in testis-biased genes, even for genes coexpressed in nonsexual (embryo) tissues. Ovary-specific genes evolved exceptionally fast, as compared to testis- or embryo-specific genes, and exhibited higher frequency of positive selection. Genes with ovary expression were preferentially involved in olfactory binding and reception. We hypothesize that at least two potential mechanisms could explain rapid evolution of ovary-biased genes in this mosquito: (1) the evolutionary rate of ovary-biased genes may be accelerated by sexual selection (including female-female competition or male-mate choice) affecting olfactory genes during female swarming by males, and/or by adaptive evolution of olfactory signaling within the female reproductive system (e.g., sperm-ovary signaling); and/or (2) testis-biased genes may exhibit decelerated evolutionary rates due to the formation of mating plugs in the female after copulation, which limits male-male sperm competition.
Collapse
|
46
|
Chang CH, Larracuente AM. Genomic changes following the reversal of a Y chromosome to an autosome in Drosophila pseudoobscura. Evolution 2017; 71:1285-1296. [PMID: 28322435 PMCID: PMC5485016 DOI: 10.1111/evo.13229] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2017] [Accepted: 03/10/2017] [Indexed: 01/10/2023]
Abstract
Robertsonian translocations resulting in fusions between sex chromosomes and autosomes shape karyotype evolution by creating new sex chromosomes from autosomes. These translocations can also reverse sex chromosomes back into autosomes, which is especially intriguing given the dramatic differences between autosomes and sex chromosomes. To study the genomic events following a Y chromosome reversal, we investigated an autosome‐Y translocation in Drosophila pseudoobscura. The ancestral Y chromosome fused to a small autosome (the dot chromosome) approximately 10–15 Mya. We used single molecule real‐time sequencing reads to assemble the D. pseudoobscura dot chromosome, including this Y‐to‐dot translocation. We find that the intervening sequence between the ancestral Y and the rest of the dot chromosome is only ∼78 Kb and is not repeat‐dense, suggesting that the centromere now falls outside, rather than between, the fused chromosomes. The Y‐to‐dot region is 100 times smaller than the D. melanogaster Y chromosome, owing to changes in repeat landscape. However, we do not find a consistent reduction in intron sizes across the Y‐to‐dot region. Instead, deletions in intergenic regions and possibly a small ancestral Y chromosome size may explain the compact size of the Y‐to‐dot translocation.
Collapse
Affiliation(s)
- Ching-Ho Chang
- Department of Biology, University of Rochester, Rochester, New York, 14627
| | | |
Collapse
|
47
|
Deb A, Johnson WA, Kline AP, Scott BJ, Meador LR, Srinivas D, Martin-Garcia JM, Dörner K, Borges CR, Misra R, Hogue BG, Fromme P, Mor TS. Bacterial expression, correct membrane targeting and functional folding of the HIV-1 membrane protein Vpu using a periplasmic signal peptide. PLoS One 2017; 12:e0172529. [PMID: 28225803 PMCID: PMC5321405 DOI: 10.1371/journal.pone.0172529] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2016] [Accepted: 02/06/2017] [Indexed: 12/04/2022] Open
Abstract
Viral protein U (Vpu) is a type-III integral membrane protein encoded by Human Immunodeficiency Virus-1 (HIV- 1). It is expressed in infected host cells and plays several roles in viral progeny escape from infected cells, including down-regulation of CD4 receptors. But key structure/function questions remain regarding the mechanisms by which the Vpu protein contributes to HIV-1 pathogenesis. Here we describe expression of Vpu in bacteria, its purification and characterization. We report the successful expression of PelB-Vpu in Escherichia coli using the leader peptide pectate lyase B (PelB) from Erwinia carotovora. The protein was detergent extractable and could be isolated in a very pure form. We demonstrate that the PelB signal peptide successfully targets Vpu to the cell membranes and inserts it as a type I membrane protein. PelB-Vpu was biophysically characterized by circular dichroism and dynamic light scattering experiments and was shown to be an excellent candidate for elucidating structural models.
Collapse
Affiliation(s)
- Arpan Deb
- School of Life Sciences, Arizona State University, Tempe, Arizona, United States of America
- The Biodesign Center for Immunotherapy, Vaccines, and Virotherapy, The Biodesign Institute, Arizona State University, Tempe, Arizona, United States of America
| | - William A. Johnson
- School of Life Sciences, Arizona State University, Tempe, Arizona, United States of America
- The Biodesign Center for Immunotherapy, Vaccines, and Virotherapy, The Biodesign Institute, Arizona State University, Tempe, Arizona, United States of America
| | - Alexander P. Kline
- School of Life Sciences, Arizona State University, Tempe, Arizona, United States of America
- The Biodesign Center for Immunotherapy, Vaccines, and Virotherapy, The Biodesign Institute, Arizona State University, Tempe, Arizona, United States of America
| | - Boston J. Scott
- School of Life Sciences, Arizona State University, Tempe, Arizona, United States of America
- The Biodesign Center for Immunotherapy, Vaccines, and Virotherapy, The Biodesign Institute, Arizona State University, Tempe, Arizona, United States of America
| | - Lydia R. Meador
- School of Life Sciences, Arizona State University, Tempe, Arizona, United States of America
- The Biodesign Center for Immunotherapy, Vaccines, and Virotherapy, The Biodesign Institute, Arizona State University, Tempe, Arizona, United States of America
| | - Dustin Srinivas
- School of Life Sciences, Arizona State University, Tempe, Arizona, United States of America
- The Biodesign Center for Immunotherapy, Vaccines, and Virotherapy, The Biodesign Institute, Arizona State University, Tempe, Arizona, United States of America
| | - Jose M. Martin-Garcia
- School of Molecular Sciences, Arizona State University, Tempe, Arizona, United States of America
- The Biodesign Center for Applied Structural Discovery, The Biodesign Institute, Arizona State University, Tempe, Arizona, United States of America
| | - Katerina Dörner
- School of Molecular Sciences, Arizona State University, Tempe, Arizona, United States of America
| | - Chad R. Borges
- School of Molecular Sciences, Arizona State University, Tempe, Arizona, United States of America
- The Biodesign Center for Personal Diagnostics, The Biodesign Institute, Arizona State University, Tempe, Arizona, United States of America
| | - Rajeev Misra
- School of Life Sciences, Arizona State University, Tempe, Arizona, United States of America
| | - Brenda G. Hogue
- School of Life Sciences, Arizona State University, Tempe, Arizona, United States of America
- The Biodesign Center for Immunotherapy, Vaccines, and Virotherapy, The Biodesign Institute, Arizona State University, Tempe, Arizona, United States of America
- The Biodesign Center for Applied Structural Discovery, The Biodesign Institute, Arizona State University, Tempe, Arizona, United States of America
| | - Petra Fromme
- School of Molecular Sciences, Arizona State University, Tempe, Arizona, United States of America
- The Biodesign Center for Applied Structural Discovery, The Biodesign Institute, Arizona State University, Tempe, Arizona, United States of America
| | - Tsafrir S. Mor
- School of Life Sciences, Arizona State University, Tempe, Arizona, United States of America
- The Biodesign Center for Immunotherapy, Vaccines, and Virotherapy, The Biodesign Institute, Arizona State University, Tempe, Arizona, United States of America
| |
Collapse
|
48
|
Romiguier J, Roux C. Analytical Biases Associated with GC-Content in Molecular Evolution. Front Genet 2017; 8:16. [PMID: 28261263 PMCID: PMC5309256 DOI: 10.3389/fgene.2017.00016] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2016] [Accepted: 02/06/2017] [Indexed: 12/19/2022] Open
Abstract
Molecular evolution is being revolutionized by high-throughput sequencing allowing an increased amount of genome-wide data available for multiple species. While base composition summarized by GC-content is one of the first metrics measured in genomes, its genomic distribution is a frequently neglected feature in downstream analyses based on DNA sequence comparisons. Here, we show how base composition heterogeneity among loci and taxa can bias common molecular evolution analyses such as phylogenetic tree reconstruction, detection of natural selection and estimation of codon usage. We then discuss the biological, technical and methodological causes of these GC-associated biases and suggest approaches to overcome them.
Collapse
Affiliation(s)
- Jonathan Romiguier
- Department of Ecology and Evolution, University of Lausanne Lausanne, Switzerland
| | - Camille Roux
- Department of Ecology and Evolution, University of Lausanne Lausanne, Switzerland
| |
Collapse
|
49
|
Pastro L, Smircich P, Di Paolo A, Becco L, Duhagon MA, Sotelo-Silveira J, Garat B. Nuclear Compartmentalization Contributes to Stage-Specific Gene Expression Control in Trypanosoma cruzi. Front Cell Dev Biol 2017; 5:8. [PMID: 28243589 PMCID: PMC5303743 DOI: 10.3389/fcell.2017.00008] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2016] [Accepted: 01/25/2017] [Indexed: 12/24/2022] Open
Abstract
In the protozoan parasite Trypanosoma cruzi, as in other trypanosomatids, transcription of protein coding genes occurs in a constitutive fashion, producing large polycistronic transcription units. These units are composed of non-functionally related genes which are pervasively processed to yield each mRNA. Therefore, post-transcriptional processes are crucial to regulate gene expression. Considering that nuclear compartmentalization could contribute to gene expression regulation, we comparatively studied the nuclear, cytoplasmic and whole cell transcriptomes of the non-infective epimastigote stage of T. cruzi, using RNA-Seq. We found that the cytoplasmic transcriptome tightly correlates with the whole cell transcriptome and both equally correlate with the proteome. Nonetheless, 1,200 transcripts showed differential abundance between the nuclear and cytoplasmic fractions. For the genes with transcript content augmented in the nucleus, significant structural and compositional differences were found. The analysis of the reported epimastigote translatome and proteome, revealed scarce ribosome footprints and encoded proteins for them. Ontology analyses unveiled that many of these genes are distinctive of other parasite life-cycle stages. Finally, the relocalization of transcript abundance in the metacyclic trypomastigote infective stage was confirmed for specific genes. While gene expression is strongly dependent on transcript steady-state level, we here highlight the importance of the distribution of transcripts abundance between compartments in T. cruzi. Particularly, we show that nuclear compartmentation is playing an active role in the developmental stage determination preventing off-stage expression.
Collapse
Affiliation(s)
- Lucía Pastro
- Laboratorio de Interacciones Moleculares, Facultad de Ciencias, Universidad de la RepúblicaMontevideo, Uruguay; Departamento de Genética, Facultad de Medicina, Universidad de la RepúblicaMontevideo, Uruguay
| | - Pablo Smircich
- Laboratorio de Interacciones Moleculares, Facultad de Ciencias, Universidad de la RepúblicaMontevideo, Uruguay; Departamento de Genética, Facultad de Medicina, Universidad de la RepúblicaMontevideo, Uruguay
| | - Andrés Di Paolo
- Departamento de Genómica, Instituto de Investigaciones Biológicas Clemente Estable Montevideo, Uruguay
| | - Lorena Becco
- Laboratorio de Interacciones Moleculares, Facultad de Ciencias, Universidad de la República Montevideo, Uruguay
| | - María A Duhagon
- Laboratorio de Interacciones Moleculares, Facultad de Ciencias, Universidad de la RepúblicaMontevideo, Uruguay; Departamento de Genética, Facultad de Medicina, Universidad de la RepúblicaMontevideo, Uruguay
| | - José Sotelo-Silveira
- Departamento de Genómica, Instituto de Investigaciones Biológicas Clemente Estable Montevideo, Uruguay
| | - Beatriz Garat
- Laboratorio de Interacciones Moleculares, Facultad de Ciencias, Universidad de la República Montevideo, Uruguay
| |
Collapse
|
50
|
The transcriptional architecture of phenotypic dimorphism. Nat Ecol Evol 2017; 1:6. [PMID: 28812569 DOI: 10.1038/s41559-016-0006] [Citation(s) in RCA: 94] [Impact Index Per Article: 13.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2016] [Accepted: 09/06/2016] [Indexed: 12/11/2022]
Abstract
The profound differences in gene expression between the sexes are increasingly used to study the molecular basis of sexual dimorphism, sexual selection and sexual conflict. Studies of transcriptional architecture, based on comparisons of gene expression, have also been implemented for a wide variety of other intra-specific polymorphisms. These efforts are based on key assumptions regarding the relationship between transcriptional architecture, phenotypic variation and the target of selection. Some of these assumptions are better supported by available evidence than others. In all cases, the evidence is largely circumstantial, leaving considerable gaps in our understanding of the relationship between transcriptional and phenotypic dimorphism.
Collapse
|