1
|
Gendron EMS, Qing X, Sevigny JL, Li H, Liu Z, Blaxter M, Powers TO, Thomas WK, Porazinska DL. Comparative mitochondrial genomics in Nematoda reveal astonishing variation in compositional biases and substitution rates indicative of multi-level selection. BMC Genomics 2024; 25:615. [PMID: 38890582 PMCID: PMC11184840 DOI: 10.1186/s12864-024-10500-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2024] [Accepted: 06/05/2024] [Indexed: 06/20/2024] Open
Abstract
BACKGROUND Nematodes are the most abundant and diverse metazoans on Earth, and are known to significantly affect ecosystem functioning. A better understanding of their biology and ecology, including potential adaptations to diverse habitats and lifestyles, is key to understanding their response to global change scenarios. Mitochondrial genomes offer high species level characterization, low cost of sequencing, and an ease of data handling that can provide insights into nematode evolutionary pressures. RESULTS Generally, nematode mitochondrial genomes exhibited similar structural characteristics (e.g., gene size and GC content), but displayed remarkable variability around these general patterns. Compositional strand biases showed strong codon position specific G skews and relationships with nematode life traits (especially parasitic feeding habits) equal to or greater than with predicted phylogeny. On average, nematode mitochondrial genomes showed low non-synonymous substitution rates, but also high clade specific deviations from these means. Despite the presence of significant mutational saturation, non-synonymous (dN) and synonymous (dS) substitution rates could still be significantly explained by feeding habit and/or habitat. Low ratios of dN:dS rates, particularly associated with the parasitic lifestyles, suggested the presence of strong purifying selection. CONCLUSIONS Nematode mitochondrial genomes demonstrated a capacity to accumulate diversity in composition, structure, and content while still maintaining functional genes. Moreover, they demonstrated a capacity for rapid evolutionary change pointing to a potential interaction between multi-level selection pressures and rapid evolution. In conclusion, this study helps establish a background for our understanding of the potential evolutionary pressures shaping nematode mitochondrial genomes, while outlining likely routes of future inquiry.
Collapse
Affiliation(s)
- Eli M S Gendron
- Department of Entomology and Nematology, University of Florida, Gainesville, FL, USA.
| | - Xue Qing
- Department of Plant Pathology, Nanjing Agricultural University, Nanjing, China.
| | - Joseph L Sevigny
- Molecular, Cellular, and Biomedical Sciences, University of New Hampshire, Durham, NH, USA
- Hubbard Center for Genome Studies, University of New Hampshire, Durham, NH, USA
| | - Hongmei Li
- Department of Plant Pathology, Nanjing Agricultural University, Nanjing, China
| | - Zhiyin Liu
- Department of Plant Pathology, Nanjing Agricultural University, Nanjing, China
| | | | - Thomas O Powers
- Department of Plant Pathology, University of Nebraska, Lincoln, NE, USA
| | - W Kelly Thomas
- Molecular, Cellular, and Biomedical Sciences, University of New Hampshire, Durham, NH, USA
- Hubbard Center for Genome Studies, University of New Hampshire, Durham, NH, USA
| | - Dorota L Porazinska
- Department of Entomology and Nematology, University of Florida, Gainesville, FL, USA
| |
Collapse
|
2
|
Moeckel C, Zaravinos A, Georgakopoulos-Soares I. Strand Asymmetries Across Genomic Processes. Comput Struct Biotechnol J 2023; 21:2036-2047. [PMID: 36968020 PMCID: PMC10030826 DOI: 10.1016/j.csbj.2023.03.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2023] [Revised: 03/08/2023] [Accepted: 03/08/2023] [Indexed: 03/12/2023] Open
Abstract
Across biological systems, a number of genomic processes, including transcription, replication, DNA repair, and transcription factor binding, display intrinsic directionalities. These directionalities are reflected in the asymmetric distribution of nucleotides, motifs, genes, transposon integration sites, and other functional elements across the two complementary strands. Strand asymmetries, including GC skews and mutational biases, have shaped the nucleotide composition of diverse organisms. The investigation of strand asymmetries often serves as a method to understand underlying biological mechanisms, including protein binding preferences, transcription factor interactions, retrotransposition, DNA damage and repair preferences, transcription-replication collisions, and mutagenesis mechanisms. Research into this subject also enables the identification of functional genomic sites, such as replication origins and transcription start sites. Improvements in our ability to detect and quantify DNA strand asymmetries will provide insights into diverse functionalities of the genome, the contribution of different mutational mechanisms in germline and somatic mutagenesis, and our knowledge of genome instability and evolution, which all have significant clinical implications in human disease, including cancer. In this review, we describe key developments that have been made across the field of genomic strand asymmetries, as well as the discovery of associated mechanisms.
Collapse
Affiliation(s)
- Camille Moeckel
- Institute for Personalized Medicine, Department of Biochemistry and Molecular Biology, The Pennsylvania State University College of Medicine, Hershey, PA, USA
| | - Apostolos Zaravinos
- Department of Life Sciences, European University Cyprus, Diogenis Str., 6, Nicosia 2404, Cyprus
- Cancer Genetics, Genomics and Systems Biology laboratory, Basic and Translational Cancer Research Center (BTCRC), Nicosia 1516, Cyprus
- Corresponding author at: Department of Life Sciences, European University Cyprus, Diogenis Str., 6, Nicosia 2404, Cyprus.
| | - Ilias Georgakopoulos-Soares
- Institute for Personalized Medicine, Department of Biochemistry and Molecular Biology, The Pennsylvania State University College of Medicine, Hershey, PA, USA
- Corresponding author.
| |
Collapse
|
3
|
Inverted base composition skews and discontinuous mitochondrial genome architecture evolution in the Enoplea (Nematoda). BMC Genomics 2022; 23:376. [PMID: 35585506 PMCID: PMC9115964 DOI: 10.1186/s12864-022-08607-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2021] [Accepted: 04/25/2022] [Indexed: 11/26/2022] Open
Abstract
Background Within the class Enoplea, the earliest-branching lineages in the phylum Nematoda, the relatively highly conserved ancestral mitochondrial architecture of Trichinellida is in stark contrast to the rapidly evolving architecture of Dorylaimida and Mermithida. To better understand the evolution of mitogenomic architecture in this lineage, we sequenced the mitogenome of a fish parasite Pseudocapillaria tomentosa (Trichinellida: Capillariidae) and compared it to all available enoplean mitogenomes. Results P. tomentosa exhibited highly reduced noncoding regions (the largest was 98 bp), and a unique base composition among the Enoplea. We attributed the latter to the inverted GC skew (0.08) in comparison to the ancestral skew in Trichinellidae (-0.43 to -0.37). Capillariidae, Trichuridae and Longidoridae (Dorylaimida) generally exhibited low negative or low positive skews (-0.1 to 0.1), whereas Mermithidae exhibited fully inverted low skews (0 to 0.05). This is indicative of inversions in the strand replication order or otherwise disrupted replication mechanism in the lineages with reduced/inverted skews. Among the Trichinellida, Trichinellidae and Trichuridae have almost perfectly conserved architecture, whereas Capillariidae exhibit multiple rearrangements of tRNA genes. In contrast, Mermithidae (Mermithida) and Longidoridae (Dorylaimida) exhibit almost no similarity to the ancestral architecture. Conclusions Longidoridae exhibited more rearranged mitogenomic architecture than the hypervariable Mermithidae. Similar to the Chromadorea, the evolution of mitochondrial architecture in enoplean nematodes exhibits a strong discontinuity: lineages possessing a mostly conserved architecture over tens of millions of years are interspersed with lineages exhibiting architectural hypervariability. As Longidoridae also have some of the smallest metazoan mitochondrial genomes, they contradict the prediction that compact mitogenomes should be structurally stable. Lineages exhibiting inverted skews appear to represent the intermediate phase between the Trichinellidae (ancestral) and fully derived skews in Chromadorean mitogenomes (GC skews = 0.18 to 0.64). Multiple lines of evidence (CAT-GTR analysis in our study, a majority of previous mitogenomic results, and skew disruption scenarios) support the Dorylaimia split into two sister-clades: Dorylaimida + Mermithida and Trichinellida. However, skew inversions produce strong base composition biases, which can hamper phylogenetic and other evolutionary studies, so enoplean mitogenomes have to be used with utmost care in evolutionary studies. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08607-4.
Collapse
|
4
|
Li Z, Ma B, Li X, Lv Y, Jiang X, Ren C, Hu C, Luo P. The Complete Mitochondrial Genome of Stichopus naso (Aspidochirotida: Stichopodidae: Stichopus) and Its Phylogenetic Position. Genes (Basel) 2022; 13:genes13050825. [PMID: 35627210 PMCID: PMC9141342 DOI: 10.3390/genes13050825] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2022] [Revised: 05/02/2022] [Accepted: 05/03/2022] [Indexed: 01/21/2023] Open
Abstract
The mitochondrial genome is widely used to study the molecular evolution of and perform phylogenetic analyses on animals. In this study, the complete mitochondrial genome (mitogenome) of Stichopus naso was sequenced. The mitogenome was 16,239 bp in length and contained 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and 2 ribosomal RNA genes (rRNAs). The genome composition showed positive AT-skew (0.023) and negative GC-skew (−0.158). The order of the mitochondrial genes was consistent with those from the Stichopus and Isostichopus species, whereas it was different from those of other species of Aspidochirotida. The phylogenetic analysis, based on the nucleotide sequences of 13 PCGs through the methods of Bayesian inference (BI) and maximum likelihood (ML), indicated that S. naso has close relationships with S. horrens and S. monotuberculatus, and belongs to a member of Stichopodidae. Our study provides a reference mitogenome for further molecular evolution studies and phylogenetic research on sea cucumbers.
Collapse
Affiliation(s)
- Zhuobo Li
- CAS Key Laboratory of Tropical Marine Bio-Resources and Ecology (LMB), Guangdong Provincial Key Laboratory of Applied Marine Biology (LAMB), South China Sea Institute of Oceanology Chinese Academy of Sciences, Guangzhou 510301, China; (Z.L.); (B.M.); (X.L.); (X.J.); (C.R.); (C.H.)
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Bo Ma
- CAS Key Laboratory of Tropical Marine Bio-Resources and Ecology (LMB), Guangdong Provincial Key Laboratory of Applied Marine Biology (LAMB), South China Sea Institute of Oceanology Chinese Academy of Sciences, Guangzhou 510301, China; (Z.L.); (B.M.); (X.L.); (X.J.); (C.R.); (C.H.)
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Xiaomin Li
- CAS Key Laboratory of Tropical Marine Bio-Resources and Ecology (LMB), Guangdong Provincial Key Laboratory of Applied Marine Biology (LAMB), South China Sea Institute of Oceanology Chinese Academy of Sciences, Guangzhou 510301, China; (Z.L.); (B.M.); (X.L.); (X.J.); (C.R.); (C.H.)
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Ying Lv
- Marin College, Beibu Gulf University, Qinzhou 535011, China;
| | - Xiao Jiang
- CAS Key Laboratory of Tropical Marine Bio-Resources and Ecology (LMB), Guangdong Provincial Key Laboratory of Applied Marine Biology (LAMB), South China Sea Institute of Oceanology Chinese Academy of Sciences, Guangzhou 510301, China; (Z.L.); (B.M.); (X.L.); (X.J.); (C.R.); (C.H.)
- Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou 510301, China
| | - Chunhua Ren
- CAS Key Laboratory of Tropical Marine Bio-Resources and Ecology (LMB), Guangdong Provincial Key Laboratory of Applied Marine Biology (LAMB), South China Sea Institute of Oceanology Chinese Academy of Sciences, Guangzhou 510301, China; (Z.L.); (B.M.); (X.L.); (X.J.); (C.R.); (C.H.)
- Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou 510301, China
| | - Chaoqun Hu
- CAS Key Laboratory of Tropical Marine Bio-Resources and Ecology (LMB), Guangdong Provincial Key Laboratory of Applied Marine Biology (LAMB), South China Sea Institute of Oceanology Chinese Academy of Sciences, Guangzhou 510301, China; (Z.L.); (B.M.); (X.L.); (X.J.); (C.R.); (C.H.)
- Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou 510301, China
| | - Peng Luo
- CAS Key Laboratory of Tropical Marine Bio-Resources and Ecology (LMB), Guangdong Provincial Key Laboratory of Applied Marine Biology (LAMB), South China Sea Institute of Oceanology Chinese Academy of Sciences, Guangzhou 510301, China; (Z.L.); (B.M.); (X.L.); (X.J.); (C.R.); (C.H.)
- Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou), Guangzhou 510301, China
- Correspondence:
| |
Collapse
|
5
|
Jakovlić I, Zou H, Zhao XM, Zhang J, Wang GT, Zhang D. Evolutionary History of Inversions in Directional Mutational Pressures in Crustacean Mitochondrial Genomes: Implications for Evolutionary Studies. Mol Phylogenet Evol 2021; 164:107288. [PMID: 34365015 DOI: 10.1016/j.ympev.2021.107288] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2021] [Revised: 07/31/2021] [Accepted: 08/03/2021] [Indexed: 01/29/2023]
Abstract
Inversions of the origin of replication (ORI) in mitochondrial genomes produce asymmetrical mutational pressures that can cause strong base composition skews. Due to skews often being overlooked, the total number of crustacean lineages that underwent ORI events remains unknown. We analysed skews, cumulative skew plots, conserved sequence motifs, and mitochondrial architecture of all 965 available crustacean mitogenomes (699 unique species). We found indications of an ORI in 159 (22.7%) species, and mapped these to 23 ORI events: 16 identified with confidence and 7 putative (13 newly proposed, and for 5 we improved the resolution). Two ORIs occurred at or above the order level: Isopoda and Copepoda. Shifts in skew plots are not a precise tool for identifying the replication mechanism. We discuss how ORIs can produce mutational bursts in mitogenomes and show how these can interfere with various types of evolutionary studies. Phylogenetic analyses were plagued by artefactual clustering, and ORI lineages exhibited longer branches, a higher number of synonymous substitutions, higher mutational saturation, and higher compositional heterogeneity. ORI events also affected codon usage and protein properties. We discuss how this may have caused erroneous interpretation of data in previous studies that did not account for skew patterns.
Collapse
Affiliation(s)
- Ivan Jakovlić
- State Key Laboratory of Grassland Agro-Ecosystem, Institute of Innovation Ecology, Lanzhou University, 730000 Lanzhou, China
| | - Hong Zou
- Key Laboratory of Aquaculture Disease Control, Ministry of Agriculture, and State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, China
| | - Xu-Mao Zhao
- State Key Laboratory of Grassland Agro-Ecosystem, Institute of Innovation Ecology, Lanzhou University, 730000 Lanzhou, China
| | - Jin Zhang
- Bio-Transduction Lab, Wuhan 430075, China
| | - Gui-Tang Wang
- Key Laboratory of Aquaculture Disease Control, Ministry of Agriculture, and State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan 430072, China; University of Chinese Academy of Sciences, Beijing 100049, China
| | - Dong Zhang
- State Key Laboratory of Grassland Agro-Ecosystem, Institute of Innovation Ecology, Lanzhou University, 730000 Lanzhou, China.
| |
Collapse
|
6
|
Mitochondrial Genomic Landscape: A Portrait of the Mitochondrial Genome 40 Years after the First Complete Sequence. Life (Basel) 2021; 11:life11070663. [PMID: 34357035 PMCID: PMC8303319 DOI: 10.3390/life11070663] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2021] [Revised: 07/02/2021] [Accepted: 07/03/2021] [Indexed: 12/11/2022] Open
Abstract
Notwithstanding the initial claims of general conservation, mitochondrial genomes are a largely heterogeneous set of organellar chromosomes which displays a bewildering diversity in terms of structure, architecture, gene content, and functionality. The mitochondrial genome is typically described as a single chromosome, yet many examples of multipartite genomes have been found (for example, among sponges and diplonemeans); the mitochondrial genome is typically depicted as circular, yet many linear genomes are known (for example, among jellyfish, alveolates, and apicomplexans); the chromosome is normally said to be “small”, yet there is a huge variation between the smallest and the largest known genomes (found, for example, in ctenophores and vascular plants, respectively); even the gene content is highly unconserved, ranging from the 13 oxidative phosphorylation-related enzymatic subunits encoded by animal mitochondria to the wider set of mitochondrial genes found in jakobids. In the present paper, we compile and describe a large database of 27,873 mitochondrial genomes currently available in GenBank, encompassing the whole eukaryotic domain. We discuss the major features of mitochondrial molecular diversity, with special reference to nucleotide composition and compositional biases; moreover, the database is made publicly available for future analyses on the MoZoo Lab GitHub page.
Collapse
|
7
|
Xia X. Improving Phylogenetic Signals of Mitochondrial Genes Using a New Method of Codon Degeneration. Life (Basel) 2020; 10:life10090171. [PMID: 32872619 PMCID: PMC7555918 DOI: 10.3390/life10090171] [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] [Received: 07/08/2020] [Revised: 08/04/2020] [Accepted: 08/27/2020] [Indexed: 12/31/2022] Open
Abstract
Recovering deep phylogeny is challenging with animal mitochondrial genes because of their rapid evolution. Codon degeneration decreases the phylogenetic noise and bias by aiming to achieve two objectives: (1) alleviate the bias associated with nucleotide composition, which may lead to homoplasy and long-branch attraction, and (2) reduce differences in the phylogenetic results between nucleotide-based and amino acid (AA)-based analyses. The discrepancy between nucleotide-based analysis and AA-based analysis is partially caused by some synonymous codons that differ more from each other at the nucleotide level than from some nonsynonymous codons, e.g., Leu codon TTR in the standard genetic code is more similar to Phe codon TTY than to synonymous CTN codons. Thus, nucleotide similarity conflicts with AA similarity. There are many such examples involving other codon families in various mitochondrial genetic codes. Proper codon degeneration will make synonymous codons more similar to each other at the nucleotide level than they are to nonsynonymous codons. Here, I illustrate a "principled" codon degeneration method that achieves these objectives. The method was applied to resolving the mammalian basal lineage and phylogenetic position of rheas among ratites. The codon degeneration method was implemented in the user-friendly and freely available DAMBE software for all known genetic codes (genetic codes 1 to 33).
Collapse
Affiliation(s)
- Xuhua Xia
- Department of Biology, University of Ottawa, 30 Marie Curie, Ottawa, ON K1N 6N5, Canada;
- Ottawa Institute of Systems Biology, 451 Smyth Road, Ottawa, ON K1H 8M5, Canada
| |
Collapse
|
8
|
Subramanian H, Gatenby RA. Evolutionary advantage of anti-parallel strand orientation of duplex DNA. Sci Rep 2020; 10:9883. [PMID: 32555277 PMCID: PMC7303137 DOI: 10.1038/s41598-020-66705-3] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/01/2020] [Accepted: 05/22/2020] [Indexed: 11/09/2022] Open
Abstract
DNA in all living systems shares common properties that are remarkably well suited to its function, suggesting refinement by evolution. However, DNA also shares some counter-intuitive properties which confer no obvious benefit, such as strand directionality and anti-parallel strand orientation, which together result in the complicated lagging strand replication. The evolutionary dynamics that led to these properties of DNA remain unknown but their universality suggests that they confer as yet unknown selective advantage to DNA. In this article, we identify an evolutionary advantage of anti-parallel strand orientation of duplex DNA, within a given set of plausible premises. The advantage stems from the increased rate of replication, achieved by dividing the DNA into predictable, independently and simultaneously replicating segments, as opposed to sequentially replicating the entire DNA, thereby parallelizing the replication process. We show that anti-parallel strand orientation is essential for such a replicative organization of DNA, given our premises, the most important of which is the assumption of the presence of sequence-dependent asymmetric cooperativity in DNA.
Collapse
Affiliation(s)
| | - Robert A Gatenby
- Integrated Mathematical Oncology Department, Cancer Biology and Evolution Program, H. Lee Moffitt Cancer Center and Research Institute, 12902, USF Magnolia Dr, Tampa, Florida, USA
| |
Collapse
|
9
|
Demongeot J, Seligmann H. Deamination gradients within codons after 1<->2 position swap predict amino acid hydrophobicity and parallel β-sheet conformational preference. Biosystems 2020; 191-192:104116. [PMID: 32081715 DOI: 10.1016/j.biosystems.2020.104116] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2019] [Revised: 12/04/2019] [Accepted: 02/10/2020] [Indexed: 12/16/2022]
Abstract
Deaminations C->T and A->G are frequent mutations producing nucleotide content gradients across genomes proportional to singlestrandedness during replication/transcription. Hence, within single codons, deamination risks increase from first to third codon positions, while second codon positions are functionally most crucial. Here genetic codes are analyzed assuming that after anticodons protected codons from deaminations, first and second codon positions swapped (N2N1N3->N1N2N3), with lowest deamination risks for N2 in presumed primitive N2N1N3 codons. N2N1N3, not standard N1N2N3, codon structure minimizes deaminations inversely proportionally to cognate amino acid hydrophobicity and parallel betasheet conformational preference. For N1N2N3, deamination minimization increases with genetic code integration order of cognate amino acids: during the presumed N2N1N3->N1N2N3 codon structure transition, protein synthesis combined direct codon-amino acid interactions for late amino acids and tRNA-based translation for early amino acids. Hence N2N1N3 codons would correspond to tRNA-free translation by spontaneous codon-amino acid affinities, and tRNA-mediated translation presumably caused N2N1N3->N1N2N3 swaps. Results show that rational, not arbitrary rules link codon and amino acid structures. Some analyses detect mitochondrial RNAs and peptides in public data corresponding to systematic position swaps, suggesting occasional swapping polymerase activity.
Collapse
Affiliation(s)
- Jacques Demongeot
- Université Grenoble Alpes, Faculty of Medicine, Laboratory AGEIS EA 7407, Team Tools for e-Gnosis Medical, F-38700, La Tronche, France.
| | - Hervé Seligmann
- Université Grenoble Alpes, Faculty of Medicine, Laboratory AGEIS EA 7407, Team Tools for e-Gnosis Medical, F-38700, La Tronche, France; The National Natural History Collections, The Hebrew University of Jerusalem, 91404, Jerusalem, Israel.
| |
Collapse
|
10
|
Noah KE, Hao J, Li L, Sun X, Foley B, Yang Q, Xia X. Major Revisions in Arthropod Phylogeny Through Improved Supermatrix, With Support for Two Possible Waves of Land Invasion by Chelicerates. Evol Bioinform Online 2020; 16:1176934320903735. [PMID: 32076367 PMCID: PMC7003163 DOI: 10.1177/1176934320903735] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2019] [Accepted: 01/02/2020] [Indexed: 01/04/2023] Open
Abstract
Deep phylogeny involving arthropod lineages is difficult to recover because the erosion of phylogenetic signals over time leads to unreliable multiple sequence alignment (MSA) and subsequent phylogenetic reconstruction. One way to alleviate the problem is to assemble a large number of gene sequences to compensate for the weakness in each individual gene. Such an approach has led to many robustly supported but contradictory phylogenies. A close examination shows that the supermatrix approach often suffers from two shortcomings. The first is that MSA is rarely checked for reliability and, as will be illustrated, can be poor. The second is that, to alleviate the problem of homoplasy at the third codon position of protein-coding genes due to convergent evolution of nucleotide frequencies, phylogeneticists may remove or degenerate the third codon position but may do it improperly and introduce new biases. We performed extensive reanalysis of one of such "big data" sets to highlight these two problems, and demonstrated the power and benefits of correcting or alleviating these problems. Our results support a new group with Xiphosura and Arachnopulmonata (Tetrapulmonata + Scorpiones) as sister taxa. This favors a new hypothesis in which the ancestor of Xiphosura and the extinct Eurypterida (sea scorpions, of which many later forms lived in brackish or freshwater) returned to the sea after the initial chelicerate invasion of land. Our phylogeny is supported even with the original data but processed with a new "principled" codon degeneration. We also show that removing the 1673 codon sites with both AGN and UCN codons (encoding serine) in our alignment can partially reconcile discrepancies between nucleotide-based and AA-based tree, partly because two sequences, one with AGN and the other with UCN, would be identical at the amino acid level but quite different at the nucleotide level.
Collapse
Affiliation(s)
| | - Jiasheng Hao
- College of Life Sciences, Anhui Normal University, Wuhu, China
| | - Luyan Li
- Nanjing Institute of Geology and Paleontology, Chinese Academy of Sciences, Nanjing, China
| | - Xiaoyan Sun
- Nanjing Institute of Geology and Paleontology, Chinese Academy of Sciences, Nanjing, China
| | - Brian Foley
- Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, Los Alamos, NM, USA
| | - Qun Yang
- Nanjing Institute of Geology and Paleontology, Chinese Academy of Sciences, Nanjing, China
| | - Xuhua Xia
- Department of Biology, University of Ottawa, Ottawa, ON, Canada
- Ottawa Institute of Systems Biology, University of Ottawa, Ottawa, ON, Canada
| |
Collapse
|
11
|
Nieuwenhuis M, van de Peppel LJJ, Bakker FT, Zwaan BJ, Aanen DK. Enrichment of G4DNA and a Large Inverted Repeat Coincide in the Mitochondrial Genomes of Termitomyces. Genome Biol Evol 2020; 11:1857-1869. [PMID: 31209489 PMCID: PMC6609731 DOI: 10.1093/gbe/evz122] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/11/2019] [Indexed: 12/20/2022] Open
Abstract
Mitochondria retain their own genome, a hallmark of their bacterial ancestry. Mitochondrial genomes (mtDNA) are highly diverse in size, shape, and structure, despite their conserved function across most eukaryotes. Exploring extreme cases of mtDNA architecture can yield important information on fundamental aspects of genome biology. We discovered that the mitochondrial genomes of a basidiomycete fungus (Termitomyces spp.) contain an inverted repeat (IR), a duplicated region half the size of the complete genome. In addition, we found an abundance of sequences capable of forming G-quadruplexes (G4DNA); structures that can disrupt the double helical formation of DNA. G4DNA is implicated in replication fork stalling, double-stranded breaks, altered gene expression, recombination, and other effects. To determine whether this occurrence of IR and G4DNA was correlated within the genus Termitomyces, we reconstructed the mitochondrial genomes of 11 additional species including representatives of several closely related genera. We show that the mtDNA of all sampled species of Termitomyces and its sister group, represented by the species Tephrocybe rancida and Blastosporella zonata, are characterized by a large IR and enrichment of G4DNA. To determine whether high mitochondrial G4DNA content is common in fungi, we conducted the first broad survey of G4DNA content in fungal mtDNA, revealing it to be a highly variable trait. The results of this study provide important direction for future research on the function and evolution of G4DNA and organellar IRs.
Collapse
Affiliation(s)
| | | | - Freek T Bakker
- Biosystematics Group, Wageningen University & Research, The Netherlands
| | - Bas J Zwaan
- Laboratory of Genetics, Wageningen University & Research, The Netherlands
| | - Duur K Aanen
- Laboratory of Genetics, Wageningen University & Research, The Netherlands
| |
Collapse
|
12
|
Pentamers with Non-redundant Frames: Bias for Natural Circular Code Codons. J Mol Evol 2020; 88:194-201. [DOI: 10.1007/s00239-019-09925-0] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2019] [Accepted: 12/17/2019] [Indexed: 02/06/2023]
|
13
|
Ma C, Wang Y, Zhang L, Li J. Mitochondrial genome characterization of the family Trigonidiidae (Orthoptera) reveals novel structural features and nad1 transcript ends. Sci Rep 2019; 9:19092. [PMID: 31836821 PMCID: PMC6911046 DOI: 10.1038/s41598-019-55740-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2019] [Accepted: 12/02/2019] [Indexed: 11/13/2022] Open
Abstract
The Trigonidiidae, a family of crickets, comprises 981 valid species with only one mitochondrial genome (mitogenome) sequenced to date. To explore mitogenome features of Trigonidiidae, six mitogenomes from its two subfamilies (Nemobiinae and Trigonidiinae) were determined. Two types of gene rearrangements involving a trnN-trnS1-trnE inversion and a trnV shuffling were shared by Trigonidiidae. A long intergenic spacer was observed between trnQ and trnM in Trigonidiinae (210-369 bp) and Nemobiinae (80-216 bp), which was capable of forming extensive stem-loop secondary structures in Trigonidiinae but not in Nemobiinae. The anticodon of trnS1 was TCT in Trigonidiinae, rather than GCT in Nemobiinae and other related subfamilies. There was no overlap between nad4 and nad4l in Dianemobius, as opposed to a conserved 7-bp overlap commonly found in insects. Furthermore, combined comparative analysis and transcript verification revealed that nad1 transcripts ended with a U, corresponding to the T immediately preceding a conserved motif GAGAC in the superfamily Grylloidea, plus poly-A tails. The resultant UAA served as a stop codon for species lacking full stop codons upstream of the motif. Our findings gain novel understanding of mitogenome structural diversity and provide insight into accurate mitogenome annotation.
Collapse
Affiliation(s)
- Chuan Ma
- Institute of Apicultural Research/Key Laboratory of Pollinating Insect Biology, Ministry of Agriculture, Chinese Academy of Agricultural Sciences, Beijing, 100093, China
| | - Yeying Wang
- Key Laboratory of State Forestry Administration on Biodiversity Conservation in Karst Mountainous Areas of Southwestern China, Guizhou Normal University, Guiyang, 550025, China
| | - Licui Zhang
- Institute of Apicultural Research/Key Laboratory of Pollinating Insect Biology, Ministry of Agriculture, Chinese Academy of Agricultural Sciences, Beijing, 100093, China
| | - Jianke Li
- Institute of Apicultural Research/Key Laboratory of Pollinating Insect Biology, Ministry of Agriculture, Chinese Academy of Agricultural Sciences, Beijing, 100093, China.
| |
Collapse
|
14
|
Li Q, Xiang D, Wan Y, Wu Q, Wu X, Ma C, Song Y, Zhao G, Huang W. The complete mitochondrial genomes of five important medicinal Ganoderma species: Features, evolution, and phylogeny. Int J Biol Macromol 2019; 139:397-408. [DOI: 10.1016/j.ijbiomac.2019.08.003] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2019] [Revised: 07/10/2019] [Accepted: 08/01/2019] [Indexed: 12/31/2022]
|
15
|
Demongeot J, Seligmann H. Theoretical minimal RNA rings designed according to coding constraints mimic deamination gradients. THE SCIENCE OF NATURE - NATURWISSENSCHAFTEN 2019; 106:44. [DOI: 10.1007/s00114-019-1638-5] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/15/2018] [Revised: 06/18/2019] [Accepted: 06/19/2019] [Indexed: 11/27/2022]
|
16
|
Sequeira-Mendes J, Vergara Z, Peiró R, Morata J, Aragüez I, Costas C, Mendez-Giraldez R, Casacuberta JM, Bastolla U, Gutierrez C. Differences in firing efficiency, chromatin, and transcription underlie the developmental plasticity of the Arabidopsis DNA replication origins. Genome Res 2019; 29:784-797. [PMID: 30846531 PMCID: PMC6499314 DOI: 10.1101/gr.240986.118] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2018] [Accepted: 02/25/2019] [Indexed: 12/20/2022]
Abstract
Eukaryotic genome replication depends on thousands of DNA replication origins (ORIs). A major challenge is to learn ORI biology in multicellular organisms in the context of growing organs to understand their developmental plasticity. We have identified a set of ORIs of Arabidopsis thaliana and their chromatin landscape at two stages of post-embryonic development. ORIs associate with multiple chromatin signatures including transcription start sites (TSS) but also proximal and distal regulatory regions and heterochromatin, where ORIs colocalize with retrotransposons. In addition, quantitative analysis of ORI activity led us to conclude that strong ORIs have high GC content and clusters of GGN trinucleotides. Development primarily influences ORI firing strength rather than ORI location. ORIs that preferentially fire at early developmental stages colocalize with GC-rich heterochromatin, but at later stages with transcribed genes, perhaps as a consequence of changes in chromatin features associated with developmental processes. Our study provides the set of ORIs active in an organism at the post-embryo stage that should allow us to study ORI biology in response to development, environment, and mutations with a quantitative approach. In a wider scope, the computational strategies developed here can be transferred to other eukaryotic systems.
Collapse
Affiliation(s)
- Joana Sequeira-Mendes
- Centro de Biologia Molecular Severo Ochoa, CSIC-UAM, Nicolas Cabrera 1, Cantoblanco, 28049 Madrid, Spain
| | - Zaida Vergara
- Centro de Biologia Molecular Severo Ochoa, CSIC-UAM, Nicolas Cabrera 1, Cantoblanco, 28049 Madrid, Spain
| | - Ramon Peiró
- Centro de Biologia Molecular Severo Ochoa, CSIC-UAM, Nicolas Cabrera 1, Cantoblanco, 28049 Madrid, Spain
| | - Jordi Morata
- Center for Research in Agricultural Genomics, CRAG (CSIC-IRTA-UAB-UB), Campus Universitat Autónoma de Barcelona, Bellaterra, Cerdanyola del Valles, 08193 Barcelona, Spain
| | - Irene Aragüez
- Centro de Biologia Molecular Severo Ochoa, CSIC-UAM, Nicolas Cabrera 1, Cantoblanco, 28049 Madrid, Spain
| | - Celina Costas
- Centro de Biologia Molecular Severo Ochoa, CSIC-UAM, Nicolas Cabrera 1, Cantoblanco, 28049 Madrid, Spain
| | - Raul Mendez-Giraldez
- Centro de Biologia Molecular Severo Ochoa, CSIC-UAM, Nicolas Cabrera 1, Cantoblanco, 28049 Madrid, Spain
| | - Josep M Casacuberta
- Center for Research in Agricultural Genomics, CRAG (CSIC-IRTA-UAB-UB), Campus Universitat Autónoma de Barcelona, Bellaterra, Cerdanyola del Valles, 08193 Barcelona, Spain
| | - Ugo Bastolla
- Centro de Biologia Molecular Severo Ochoa, CSIC-UAM, Nicolas Cabrera 1, Cantoblanco, 28049 Madrid, Spain
| | - Crisanto Gutierrez
- Centro de Biologia Molecular Severo Ochoa, CSIC-UAM, Nicolas Cabrera 1, Cantoblanco, 28049 Madrid, Spain
| |
Collapse
|
17
|
Le Ret M, Belcher S, Graindorge S, Wallet C, Koechler S, Erhardt M, Williams-Carrier R, Barkan A, Gualberto JM. Efficient Replication of the Plastid Genome Requires an Organellar Thymidine Kinase. PLANT PHYSIOLOGY 2018; 178:1643-1656. [PMID: 30305373 PMCID: PMC6288739 DOI: 10.1104/pp.18.00976] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/18/2018] [Accepted: 10/01/2018] [Indexed: 05/17/2023]
Abstract
Thymidine kinase (TK) is a key enzyme of the salvage pathway that recycles thymidine nucleosides to produce deoxythymidine triphosphate. Here, we identified the single TK of maize (Zea mays), denoted CPTK1, as necessary in the replication of the plastidial genome (cpDNA), demonstrating the essential function of the salvage pathway during chloroplast biogenesis. CPTK1 localized to both plastids and mitochondria, and its absence resulted in an albino phenotype, reduced cpDNA copy number and a severe deficiency in plastidial ribosomes. Mitochondria were not affected, indicating they are less reliant on the salvage pathway. Arabidopsis (Arabidopsis thaliana) TKs, TK1A and TK1B, apparently resulted from a gene duplication after the divergence of monocots and dicots. Similar but less-severe effects were observed for Arabidopsis tk1a tk1b double mutants in comparison to those in maize cptk1 TK1B was important for cpDNA replication and repair in conditions of replicative stress but had little impact on the mitochondrial phenotype. In the maize cptk1 mutant, the DNA from the small single-copy region of the plastidial genome was reduced to a greater extent than other regions, suggesting preferential abortion of replication in this region. This was accompanied by the accumulation of truncated genomes that resulted, at least in part, from unfaithful microhomology-mediated repair. These and other results suggest that the loss of normal cpDNA replication elicits the mobilization of new replication origins around the rpoB (beta subunit of plastid-encoded RNA polymerase) transcription unit and imply that increased transcription at rpoB is associated with the initiation of cpDNA replication.
Collapse
Affiliation(s)
- Monique Le Ret
- Institut de Biologie Moléculaire des Plantes, CNRS-UPR2357, Université de Strasbourg, 67084 Strasbourg, France
| | - Susan Belcher
- Institute of Molecular Biology, University of Oregon, Eugene, Oregon 97403
| | - Stéfanie Graindorge
- Institut de Biologie Moléculaire des Plantes, CNRS-UPR2357, Université de Strasbourg, 67084 Strasbourg, France
| | - Clémentine Wallet
- Institut de Biologie Moléculaire des Plantes, CNRS-UPR2357, Université de Strasbourg, 67084 Strasbourg, France
| | - Sandrine Koechler
- Institut de Biologie Moléculaire des Plantes, CNRS-UPR2357, Université de Strasbourg, 67084 Strasbourg, France
| | - Mathieu Erhardt
- Institut de Biologie Moléculaire des Plantes, CNRS-UPR2357, Université de Strasbourg, 67084 Strasbourg, France
| | | | - Alice Barkan
- Institute of Molecular Biology, University of Oregon, Eugene, Oregon 97403
| | - José M Gualberto
- Institut de Biologie Moléculaire des Plantes, CNRS-UPR2357, Université de Strasbourg, 67084 Strasbourg, France
| |
Collapse
|
18
|
Li Q, Liao M, Yang M, Xiong C, Jin X, Chen Z, Huang W. Characterization of the mitochondrial genomes of three species in the ectomycorrhizal genus Cantharellus and phylogeny of Agaricomycetes. Int J Biol Macromol 2018; 118:756-769. [DOI: 10.1016/j.ijbiomac.2018.06.129] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2018] [Revised: 06/23/2018] [Accepted: 06/26/2018] [Indexed: 12/15/2022]
|
19
|
Abstract
Codon usage depends on mutation bias, tRNA-mediated selection, and the need for high efficiency and accuracy in translation. One codon in a synonymous codon family is often strongly over-used, especially in highly expressed genes, which often leads to a high dN/dS ratio because dS is very small. Many different codon usage indices have been proposed to measure codon usage and codon adaptation. Sense codon could be misread by release factors and stop codons misread by tRNAs, which also contribute to codon usage in rare cases. This chapter outlines the conceptual framework on codon evolution, illustrates codon-specific and gene-specific codon usage indices, and presents their applications. A new index for codon adaptation that accounts for background mutation bias (Index of Translation Elongation) is presented and contrasted with codon adaptation index (CAI) which does not consider background mutation bias. They are used to re-analyze data from a recent paper claiming that translation elongation efficiency matters little in protein production. The reanalysis disproves the claim.
Collapse
|
20
|
Luo H, Quan CL, Peng C, Gao F. Recent development of Ori-Finder system and DoriC database for microbial replication origins. Brief Bioinform 2018; 20:1114-1124. [DOI: 10.1093/bib/bbx174] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2017] [Revised: 12/04/2017] [Indexed: 01/28/2023] Open
Abstract
Abstract
DNA replication begins at replication origins in all three domains of life. Identification and characterization of replication origins are important not only in providing insights into the structure and function of the replication origins but also in understanding the regulatory mechanisms of the initiation step in DNA replication. The Z-curve method has been used in the identification of replication origins in archaeal genomes successfully since 2002. Furthermore, the Web servers of Ori-Finder and Ori-Finder 2 have been developed to predict replication origins in both bacterial and archaeal genomes based on the Z-curve method, and the replication origins with manual curation have been collected into an online database, DoriC. Ori-Finder system and DoriC database are currently used in the research field of DNA replication origins in prokaryotes, including: (i) identification of oriC regions in bacterial and archaeal genomes; (ii) discovery and analysis of the conserved sequences within oriC regions; and (iii) strand-biased analysis of bacterial genomes.
Up to now, more and more predicted results by Ori-Finder system were supported by subsequent experiments, and Ori-Finder system has been used to identify the replication origins in > 100 newly sequenced prokaryotes in their genome reports. In addition, the data in DoriC database have been widely used in the large-scale analyses of replication origins and strand bias in prokaryotic genomes. Here, we review the development of Ori-Finder system and DoriC database as well as their applications. Some future directions and aspects for extending the application of Ori-Finder and DoriC are also presented.
Collapse
|
21
|
Sadhasivam A, Vetrivel U. Genome-wide codon usage profiling of ocular infective Chlamydia trachomatis serovars and drug target identification. J Biomol Struct Dyn 2017. [PMID: 28627970 DOI: 10.1080/07391102.2017.1343685] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Chlamydia trachomatis (C.t) is a Gram-negative obligate intracellular bacteria and is a major causative of infectious blindness and sexually transmitted diseases. Among the varied serovars of this organism, A, B and C are reported as prominent ocular pathogens. Genomic studies of these strains shall aid in deciphering potential drug targets and genomic influence on pathogenesis. Hence, in this study we performed deep statistical profiling of codon usage in these serovars. The overall base composition analysis reveals that these serovars are over biased to AU than GC. Similarly, relative synonymous codon usage also showed preference towards A/U ending codons. Parity Rule 2 analysis inferred unequal distribution of AT and GC, indicative of other unknown factors acting along with mutational pressure to influence codon usage bias (CUB). Moreover, absolute quantification of CUB also revealed lower bias across these serovars. The effect of natural selection on CUB was also confirmed by neutrality plot, reinforcing natural selection under mutational pressure turned to be a pivotal role in shaping the CUB in the strains studied. Correspondence analysis (COA) clarified that, C.t C/TW-3 to show a unique trend in codon usage variation. Host influence analysis on shaping the codon usage pattern also inferred some speculative relativity. In a nutshell, our finding suggests that mutational pressure is the dominating factor in shaping CUB in the strains studied, followed by natural selection. We also propose potential drug targets based on cumulative analysis of strand bias, CUB and human non-homologue screening.
Collapse
Affiliation(s)
- Anupriya Sadhasivam
- a Centre for Bioinformatics , Kamalnayan Bajaj Institute for Research in Vision and Ophthalmology, Vision Research Foundation, Sankara Nethralaya , Chennai 600 006 , Tamil Nadu , India
| | - Umashankar Vetrivel
- a Centre for Bioinformatics , Kamalnayan Bajaj Institute for Research in Vision and Ophthalmology, Vision Research Foundation, Sankara Nethralaya , Chennai 600 006 , Tamil Nadu , India
| |
Collapse
|
22
|
The linear plastid chromosomes of maize: terminal sequences, structures, and implications for DNA replication. Curr Genet 2015; 62:431-42. [PMID: 26650613 DOI: 10.1007/s00294-015-0548-0] [Citation(s) in RCA: 48] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2015] [Revised: 11/15/2015] [Accepted: 11/22/2015] [Indexed: 01/13/2023]
Abstract
The structure of a chromosomal DNA molecule may influence the way in which it is replicated and inherited. For decades plastid DNA (ptDNA) was believed to be circular, with breakage invoked to explain linear forms found upon extraction from the cell. Recent evidence indicates that ptDNA in vivo consists of linear molecules with discrete termini, although these ends were not characterized. We report the sequences of two terminal regions, End1 and End2, for maize (Zea mays L.) ptDNA. We describe structural features of these terminal regions and similarities found in other plant ptDNAs. The terminal sequences are within inverted repeat regions (leading to four genomic isomers) and adjacent to origins of replication. Conceptually, stem-loop structures may be formed following melting of the double-stranded DNA ends. Exonuclease digestion indicates that the ends in maize are unobstructed, but tobacco (Nicotiana tabacum L.) ends may have a 5'-protein. If the terminal structure of ptDNA molecules influences the retention of ptDNA, the unprotected molecular ends in mature leaves of maize may be more susceptible to degradation in vivo than the protected ends in tobacco. The terminal sequences and cumulative GC skew profiles are nearly identical for maize, wheat (Triticum aestivum L.) and rice (Oryza sativa L.), with less similarity among other plants. The linear structure is now confirmed for maize ptDNA and inferred for other plants and suggests a virus-like recombination-dependent replication mechanism for ptDNA. Plastid transformation vectors containing the terminal sequences may increase the chances of success in generating transplastomic cereals.
Collapse
|
23
|
Multiple Factors Drive Replicating Strand Composition Bias in Bacterial Genomes. Int J Mol Sci 2015; 16:23111-26. [PMID: 26404268 PMCID: PMC4613354 DOI: 10.3390/ijms160923111] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2015] [Revised: 08/18/2015] [Accepted: 09/18/2015] [Indexed: 11/18/2022] Open
Abstract
Composition bias from Chargaff’s second parity rule (PR2) has long been found in sequenced genomes, and is believed to relate strongly with the replication process in microbial genomes. However, some disagreement on the underlying reason for strand composition bias remains. We performed an integrative analysis of various genomic features that might influence composition bias using a large-scale dataset of 1111 genomes. Our results indicate (1) the bias was stronger in obligate intracellular bacteria than in other free-living species (p-value = 0.0305); (2) Fusobacteria and Firmicutes had the highest average bias among the 24 microbial phyla analyzed; (3) the strength of selected codon usage bias and generation times were not observably related to strand composition bias (p-value = 0.3247); (4) significant negative relationships were found between GC content, genome size, rearrangement frequency, Clusters of Orthologous Groups (COG) functional subcategories A, C, I, Q, and composition bias (p-values < 1.0 × 10−8); (5) gene density and COG functional subcategories D, F, J, L, and V were positively related with composition bias (p-value < 2.2 × 10−16); and (6) gene density made the most important contribution to composition bias, indicating transcriptional bias was associated strongly with strand composition bias. Therefore, strand composition bias was found to be influenced by multiple factors with varying weights.
Collapse
|
24
|
Gao F. Bacteria may have multiple replication origins. Front Microbiol 2015; 6:324. [PMID: 25941523 PMCID: PMC4403523 DOI: 10.3389/fmicb.2015.00324] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2015] [Accepted: 03/31/2015] [Indexed: 01/15/2023] Open
Affiliation(s)
- Feng Gao
- Department of Physics, Tianjin University Tianjin, China ; Key Laboratory of Systems Bioengineering (Ministry of Education), Tianjin University Tianjin, China ; SynBio Research Platform, Collaborative Innovation Center of Chemical Science and Engineering Tianjin, China
| |
Collapse
|
25
|
Maternal ancestry and population history from whole mitochondrial genomes. INVESTIGATIVE GENETICS 2015; 6:3. [PMID: 25798216 PMCID: PMC4367903 DOI: 10.1186/s13323-015-0022-2] [Citation(s) in RCA: 69] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/26/2014] [Accepted: 02/04/2015] [Indexed: 01/12/2023]
Abstract
MtDNA has been a widely used tool in human evolutionary and population genetic studies over the past three decades. Its maternal inheritance and lack of recombination have offered the opportunity to explore genealogical relationships among individuals and to study the frequency differences of matrilineal clades among human populations at continental and regional scales. The whole mtDNA genome sequencing delivers molecular resolution that is sufficient to distinguish patterns that have arisen over thousands of years. However, mutation rate is highly variable among the functional and non-coding domains of mtDNA which makes it challenging to obtain accurate split dates of the mitochondrial clades. Due to the shallow coalescent time of mitochondrial TMRCA at approximately 100 to 200 thousand years (ky), mtDNA data have only limited power to inform us about the more distant past and the early stages of human evolutionary history. The variation shared by mitochondrial genomes of individuals drawn from different continents outside Africa has been used to illuminate the details of the colonization process of the Old World, whereas regional patterns of variation have been at the focus of studies addressing questions of a more recent time scale. In the era of whole nuclear genome sequencing, mitochondrial genomes are continuing to be informative as a unique tool for the assessment of female-specific aspects of the demographic history of human populations.
Collapse
|
26
|
Prabhakaran R, Chithambaram S, Xia X. Escherichia coli and Staphylococcus phages: effect of translation initiation efficiency on differential codon adaptation mediated by virulent and temperate lifestyles. J Gen Virol 2015; 96:1169-1179. [PMID: 25614589 PMCID: PMC4631060 DOI: 10.1099/vir.0.000050] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2014] [Accepted: 01/11/2015] [Indexed: 12/19/2022] Open
Abstract
Rapid biosynthesis is key to the success of bacteria and viruses. Highly expressed genes in bacteria exhibit a strong codon bias corresponding to the differential availability of tRNAs. However, a large clade of lambdoid coliphages exhibits relatively poor codon adaptation to the host translation machinery, in contrast to other coliphages that exhibit strong codon adaptation to the host. Three possible explanations were previously proposed but dismissed: (1) the phage-borne tRNA genes that reduce the dependence of phage translation on host tRNAs, (2) lack of time needed for evolving codon adaptation due to recent host switching, and (3) strong strand asymmetry with biased mutation disrupting codon adaptation. Here, we examined the possibility that phages with relatively poor codon adaptation have poor translation initiation which would weaken the selection on codon adaptation. We measured translation initiation by: (1) the strength and position of the Shine-Dalgarno (SD) sequence, and (2) the stability of the secondary structure of sequences flanking the SD and start codon known to affect accessibility of the SD sequence and start codon. Phage genes with strong codon adaptation had significantly stronger SD sequences than those with poor codon adaptation. The former also had significantly weaker secondary structure in sequences flanking the SD sequence and start codon than the latter. Thus, lambdoid phages do not exhibit strong codon adaptation because they have relatively inefficient translation initiation and would benefit little from increased elongation efficiency. We also provided evidence suggesting that phage lifestyle (virulent versus temperate) affected selection intensity on the efficiency of translation initiation and elongation.
Collapse
Affiliation(s)
- Ramanandan Prabhakaran
- Department of Biology and Center for Advanced Research in Environmental Genomics, University of Ottawa, 30 Marie Curie, PO Box 450, Station A, Ottawa, Ontario K1N 6N5, Canada
| | - Shivapriya Chithambaram
- Department of Biology and Center for Advanced Research in Environmental Genomics, University of Ottawa, 30 Marie Curie, PO Box 450, Station A, Ottawa, Ontario K1N 6N5, Canada
| | - Xuhua Xia
- Department of Biology and Center for Advanced Research in Environmental Genomics, University of Ottawa, 30 Marie Curie, PO Box 450, Station A, Ottawa, Ontario K1N 6N5, Canada
- Correspondence Xuhua Xia
| |
Collapse
|
27
|
Fonseca MM, Harris DJ, Posada D. The inversion of the Control Region in three mitogenomes provides further evidence for an asymmetric model of vertebrate mtDNA replication. PLoS One 2014; 9:e106654. [PMID: 25268704 PMCID: PMC4182315 DOI: 10.1371/journal.pone.0106654] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2014] [Accepted: 08/04/2014] [Indexed: 11/29/2022] Open
Abstract
Mitochondrial genomes are known to have a strong strand-specific compositional bias that is more pronounced at fourfold redundant sites of mtDNA protein-coding genes. This observation suggests that strand asymmetries, to a large extent, are caused by mutational asymmetric mechanisms. In vertebrate mitogenomes, replication and not transcription seems to play a major role in shaping compositional bias. Hence, one can better understand how mtDNA is replicated – a debated issue – through a detailed picture of mitochondrial genome evolution. Here, we analyzed the compositional bias (AT and GC skews) in protein-coding genes of almost 2,500 complete vertebrate mitogenomes. We were able to identify three fish mitogenomes with inverted AT/GC skew coupled with an inversion of the Control Region. These findings suggest that the vertebrate mitochondrial replication mechanism is asymmetric and may invert its polarity, with the leading-strand becoming the lagging-strand and vice-versa, without compromising mtDNA maintenance and expression. The inversion of the strand-specific compositional bias through the inversion of the Control Region is in agreement with the strand-displacement model but it is also compatible with the RITOLS model of mtDNA replication.
Collapse
Affiliation(s)
- Miguel M. Fonseca
- Department of Biochemistry, Genetics and Immunology, University of Vigo, Vigo, Spain
- CIBIO/InBIO, Research Center in Biodiversity and Genetic Resources, University of Porto, Vairão, Portugal
- * E-mail:
| | - D. James Harris
- CIBIO/InBIO, Research Center in Biodiversity and Genetic Resources, University of Porto, Vairão, Portugal
| | - David Posada
- Department of Biochemistry, Genetics and Immunology, University of Vigo, Vigo, Spain
| |
Collapse
|
28
|
Nolan JR, Bergthorsson U, Adema CM. Physella acuta: atypical mitochondrial gene order among panpulmonates (Gastropoda). ACTA ACUST UNITED AC 2014; 80:388-399. [PMID: 25368439 PMCID: PMC4214460 DOI: 10.1093/mollus/eyu025] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2013] [Accepted: 03/14/2014] [Indexed: 11/12/2022]
Abstract
Mitochondrial (mt) sequences are frequently used for phylogenetic reconstruction and for identification of species of molluscs. This study expands the phylogenetic range of Hygrophila (Panpulmonata) for which such sequence data are available by characterizing the full mt genome of the invasive freshwater snail Physella acuta (Physidae). The mt genome sequences of two P. acuta isolates from Stubblefield Lake, New Mexico, USA, differed in length (14,490 vs 14,314 bp) and showed 11.49% sequence divergence, whereas ITS1 and ITS2 sequences from the nuclear genome differed by 1.75%. The mt gene order of P. acuta (cox1, P, nad6, nad5, nad1, D, F, cox2, Y, W, nad4L, C, Q, atp6, R, E, rrnS, M, T, cox3, I, nad2, K, V, rrnL, L1, A, cytb, G, H, L2, atp8, N, nad2, S1, S2, nad4) differs considerably from the relatively conserved gene order within Panpulmonata. Phylogenetic trees show that the 13 protein-encoding mt gene sequences (equivalent codons) of P. acuta group according to gastropod phylogeny, yet branch lengths and dN/dS ratios for P. acuta indicate elevated amino acid substitutions relative to other gastropods. This study indicates that mt sequences of P. acuta are phylogenetically informative despite a considerable intraspecific divergence and the atypical gene order in its mt genome.
Collapse
Affiliation(s)
- Journey R Nolan
- Center for Evolutionary and Theoretical Immunology (CETI), Department of Biology MSC03 2020 , University of New Mexico , 1 University Blvd NE, Albuquerque, NM 87131 , USA
| | - Ulfar Bergthorsson
- Center for Evolutionary and Theoretical Immunology (CETI), Department of Biology MSC03 2020 , University of New Mexico , 1 University Blvd NE, Albuquerque, NM 87131 , USA
| | - Coen M Adema
- Center for Evolutionary and Theoretical Immunology (CETI), Department of Biology MSC03 2020 , University of New Mexico , 1 University Blvd NE, Albuquerque, NM 87131 , USA
| |
Collapse
|
29
|
Kim T, Kim J, Cho S, Min GS, Park C, Carreno RA, Nadler SA, Park JK. Phylogeny of Rhigonematomorpha based on the complete mitochondrial genome ofRhigonema thysanophora(Nematoda: Chromadorea). ZOOL SCR 2014. [DOI: 10.1111/zsc.12047] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Taeho Kim
- Program in Cell Biology and Genetics; College of Medicine; Chungbuk National University; Cheongju 361-763 Korea
| | - Jiyeon Kim
- Program in Cell Biology and Genetics; College of Medicine; Chungbuk National University; Cheongju 361-763 Korea
| | - Soowon Cho
- Department of Plant Medicine; Chungbuk National University; Cheongju 361-763 Korea
| | - Gi-Sik Min
- Department of Biological Sciences; Inha University; Incheon 402-751 Korea
| | - Chungoo Park
- School of Biological Sciences and Technology; Chonnam National University; GwangJu 500-757 Korea
| | - Ramon A. Carreno
- Department of Zoology; Ohio Wesleyan University; Delaware OH 43015 USA
| | - Steven A. Nadler
- Department of Entomology and Nematology; University of California; Davis CA 95616 USA
| | - Joong-Ki Park
- Program in Cell Biology and Genetics; College of Medicine; Chungbuk National University; Cheongju 361-763 Korea
| |
Collapse
|
30
|
Abstract
Studying phage codon adaptation is important not only for understanding the process of translation elongation, but also for reengineering phages for medical and industrial purposes. To evaluate the effect of mutation and selection on phage codon usage, we developed an index to measure selection imposed by host translation machinery, based on the difference in codon usage between all host genes and highly expressed host genes. We developed linear and nonlinear models to estimate the C→T mutation bias in different phage lineages and to evaluate the relative effect of mutation and host selection on phage codon usage. C→T-biased mutations occur more frequently in single-stranded DNA (ssDNA) phages than in double-stranded DNA (dsDNA) phages and affect not only synonymous codon usage, but also nonsynonymous substitutions at second codon positions, especially in ssDNA phages. The host translation machinery affects codon adaptation in both dsDNA and ssDNA phages, with a stronger effect on dsDNA phages than on ssDNA phages. Strand asymmetry with the associated local variation in mutation bias can significantly interfere with codon adaptation in both dsDNA and ssDNA phages.
Collapse
|
31
|
Xia X. Phylogenetic Bias in the Likelihood Method Caused by Missing Data Coupled with Among-Site Rate Variation: An Analytical Approach. BIOINFORMATICS RESEARCH AND APPLICATIONS 2014. [DOI: 10.1007/978-3-319-08171-7_2] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
|
32
|
Kennedy SR, Salk JJ, Schmitt MW, Loeb LA. Ultra-sensitive sequencing reveals an age-related increase in somatic mitochondrial mutations that are inconsistent with oxidative damage. PLoS Genet 2013; 9:e1003794. [PMID: 24086148 PMCID: PMC3784509 DOI: 10.1371/journal.pgen.1003794] [Citation(s) in RCA: 248] [Impact Index Per Article: 22.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2013] [Accepted: 07/29/2013] [Indexed: 01/12/2023] Open
Abstract
Mitochondrial DNA (mtDNA) is believed to be highly vulnerable to age-associated damage and mutagenesis by reactive oxygen species (ROS). However, somatic mtDNA mutations have historically been difficult to study because of technical limitations in accurately quantifying rare mtDNA mutations. We have applied the highly sensitive Duplex Sequencing methodology, which can detect a single mutation among >107 wild type molecules, to sequence mtDNA purified from human brain tissue from both young and old individuals with unprecedented accuracy. We find that the frequency of point mutations increases ∼5-fold over the course of 80 years of life. Overall, the mutation spectra of both groups are comprised predominantly of transition mutations, consistent with misincorporation by DNA polymerase γ or deamination of cytidine and adenosine as the primary mutagenic events in mtDNA. Surprisingly, G→T mutations, considered the hallmark of oxidative damage to DNA, do not significantly increase with age. We observe a non-uniform, age-independent distribution of mutations in mtDNA, with the D-loop exhibiting a significantly higher mutation frequency than the rest of the genome. The coding regions, but not the D-loop, exhibit a pronounced asymmetric accumulation of mutations between the two strands, with G→A and T→C mutations occurring more often on the light strand than the heavy strand. The patterns and biases we observe in our data closely mirror the mutational spectrum which has been reported in studies of human populations and closely related species. Overall our results argue against oxidative damage being a major driver of aging and suggest that replication errors by DNA polymerase γ and/or spontaneous base hydrolysis are responsible for the bulk of accumulating point mutations in mtDNA. Owing to their evolutionary history, mitochondria harbor independently replicating genomes. Failure to faithfully transmit the genetic information of mtDNA during replication can lead to the production of dysfunctional electron transport proteins and a subsequent decline in energy production. Cellularly-derived reactive oxygen species (ROS) and environmental agents preferentially damage mtDNA compared to nuclear DNA. However, little is known about the consequences of mtDNA damage for mutagenesis. This lack of knowledge stems, in part, from an absence of methods capable of accurately detecting these mutations throughout the mitochondrial genome. Using a new, highly sensitive DNA sequencing strategy, we find that the frequency of point mutations is 10–100-fold lower than what has been previously reported using less precise means. Moreover, the frequency increases 5-fold over an 80 year lifespan. We also find that it is predominantly transition mutations, rather than mutations commonly associated with oxidative damage to mtDNA, that increase with age. This finding is inconsistent with free radical theories of aging. Finally, the mutagenic patterns and biases we observe in our data are similar to what is seen in population studies of mitochondrial polymorphisms and suggest a common mechanism by which somatic and germline mtDNA mutations arise.
Collapse
Affiliation(s)
- Scott R. Kennedy
- Department of Pathology, University of Washington, Seattle, Washington, United States of America
| | - Jesse J. Salk
- Department of Pathology, University of Washington, Seattle, Washington, United States of America
- Department of Medicine, University of Washington, Seattle, Washington, United States of America
| | - Michael W. Schmitt
- Department of Pathology, University of Washington, Seattle, Washington, United States of America
- Department of Medicine, University of Washington, Seattle, Washington, United States of America
| | - Lawrence A. Loeb
- Department of Pathology, University of Washington, Seattle, Washington, United States of America
- Department of Biochemistry, University of Washington, Seattle, Washington, United States of America
- * E-mail:
| |
Collapse
|
33
|
Xia X. DAMBE5: a comprehensive software package for data analysis in molecular biology and evolution. Mol Biol Evol 2013; 30:1720-8. [PMID: 23564938 PMCID: PMC3684854 DOI: 10.1093/molbev/mst064] [Citation(s) in RCA: 739] [Impact Index Per Article: 67.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Since its first release in 2001 as mainly a software package for phylogenetic analysis, data analysis for molecular biology and evolution (DAMBE) has gained many new functions that may be classified into six categories: 1) sequence retrieval, editing, manipulation, and conversion among more than 20 standard sequence formats including MEGA, NEXUS, PHYLIP, GenBank, and the new NeXML format for interoperability, 2) motif characterization and discovery functions such as position weight matrix and Gibbs sampler, 3) descriptive genomic analysis tools with improved versions of codon adaptation index, effective number of codons, protein isoelectric point profiling, RNA and protein secondary structure prediction and calculation of minimum folding energy, and genomic skew plots with optimized window size, 4) molecular phylogenetics including sequence alignment, testing substitution saturation, distance-based, maximum parsimony, and maximum-likelihood methods for tree reconstructions, testing the molecular clock hypothesis with either a phylogeny or with relative-rate tests, dating gene duplication and speciation events, choosing the best-fit substitution models, and estimating rate heterogeneity over sites, 5) phylogeny-based comparative methods for continuous and discrete variables, and 6) graphic functions including secondary structure display, optimized skew plot, hydrophobicity plot, and many other plots of amino acid properties along a protein sequence, tree display and drawing by dragging nodes to each other, and visual searching of the maximum parsimony tree. DAMBE features a graphic, user-friendly, and intuitive interface and is freely available from http://dambe.bio.uottawa.ca (last accessed April 16, 2013).
Collapse
Affiliation(s)
- Xuhua Xia
- Department of Biology and Center for Advanced Research in Environmental Genomics, University of Ottawa, Ottawa, Ontario, Canada
| |
Collapse
|