Accelerated Laboratory Evolution Reveals the Influence of Replication on the GC Skew in Escherichia coli

Terrestrial Adaptation in Chelonoidis vicina as Revealed Based on Analysis of the Complete Mitochondrial Genome

BACKGROUND/OBJECTIVES

Mitochondrial genomes are widely used in phylogenetics and evolutionary and ecological research.

METHODS

In this study, the newest mitochondrial genome of Chelonoidis vicina was assembled and annotated. The comparative mitochondrial genome and selection pressure analyses were used to examine the terrestrial adaptive evolution characteristics of C. vicina and other terrestrial reptiles.

RESULTS

The results reveal that the mitochondrial genome of the tortoise C. vicina is consistent with that of other tortoise species, comprising 13 protein-coding genes (PCGs), 2 rRNAs, 22 tRNAs, and 1 noncoding control region (CR). The analysis of selection pressure reveals the presence of positive selection sites in the COX2, COX3, Cytb, ND3, ND4, ND4L, ND5, and ND6 genes of terrestrial reptiles. Of these, the COX2 and ND3 genes exhibited faster evolutionary rates. The mitochondrial genome structure of C. vicina is consistent with that of different terrestrial reptiles. The positive selection sites of COX2 and ND3 in terrestrial reptiles are closely related to a change in mitochondrial energy metabolism, which is possibly related to terrestrial adaptability.

CONCLUSIONS

The results of this study provide new insights into the adaptive evolution of C. vicina to terrestrial niches from a mitogenomic perspective, as well as genetic resources for the protection of C. vicina.

The complete mitochondrial genome of Brachytarsina amboinensis (Diptera: Hippoboscoidea: Streblidae) provides new insights into phylogenetic relationships of Hippoboscoidea

The family Streblidae is a significant grouping of dipteran insects within the superfamily Hippoboscoidea, which parasitizes the body surface of bats. With the global spread of bat-related pathogens in recent years, Streblidae has gained increasing attention due to its potential for pathogen transmission. A sample of Brachytarsina amboinensis was sequenced on the B. amboinensis were obtained, compared with available Streblidae mitogenomes, and the phylogeny of Hippoboscoidea was reconstructed. The results indicate that the mitochondrial genome of B. amboinensis exhibits a relatively high degree of conservation, with an identical gene count, arrangement, and orientation as the ancestral insect's genome. Base composition analysis revealed a strong bias towards A and T in the base composition. Selection pressure analysis indicated strong purifying selection acting on cox1. Pairwise genetic distance analysis showed that cox1 evolved at a relatively slow rate. Regarding phylogenetic relationships, the constructed phylogenetic trees using Bayesian inference and Maximum Likelihood methods supported the monophyly of the Hippoboscoidea, Glossinidae, Hippoboscidae, and Nycteribiidae clades, with high nodal support values. Our research confirmed the paraphyly of the families Streblidae. In the familial relations between Nycteribiidae and Streblidae, New World Streblidae share a closer kinship with Nycteribiidae. This contrasts with prior findings which indicated that Old World Streblidae share a closer kinship with Nycteribiidae. This study not only enhances the molecular database for bat flies but also provides a valuable reference for the identification and phylogenetic analysis of Streblidae.

Analysis of the Complete Mitochondrial Genome of Pteronura brasiliensis and Lontra canadensis

P. brasiliensis and L. canadensis are two otter species, which successfully occupied semi-aquatic habitats and diverged from other Mustelidae. Herein, the full-length mitochondrial genome sequences were constructed for these two otter species for the first time. Comparative mitochondrial genome, selection pressure, and phylogenetic independent contrasts (PICs) analyses were conducted to determine the structure and evolutionary characteristics of their mitochondrial genomes. Phylogenetic analyses were also conducted to confirm these two otter species' phylogenetic position. The results demonstrated that the mitochondrial genome structure of P. brasiliensis and L. canadensis were consistent across Mustelidae. However, selection pressure analyses demonstrated that the evolutionary rates of mitochondrial genome protein-coding genes (PCGs) ND1, ND4, and ND4L were higher in otters than in terrestrial Mustelidae, whereas the evolutionary rates of ND2, ND6, and COX1 were lower in otters. Additionally, PIC analysis demonstrated that the evolutionary rates of ND2, ND4, and ND4L markedly correlated with a niche type. Phylogenetic analysis showed that P. brasiliensis is situated at the base of the evolutionary tree of otters, and then L. canadensis diverged from it. This study suggests a divergent evolutionary pattern of Mustelidae mitochondrial genome PCGs, prompting the otters' adaptation to semi-aquatic habitats.

Shared properties of gene transfer agent and core genes revealed by comparative genomics of Alphaproteobacteria

Gene transfer agents (GTAs) are phage-like particles that transfer pieces of cellular genomic DNA to other cells. Homologues of the Rhodobacter capsulatus GTA (RcGTA) structural genes are widely distributed in the Alphaproteobacteria and particularly well conserved in the order Rhodobacterales. Possible reasons for their widespread conservation are still being discussed. It has been suggested that these alphaproteobacterial elements originate from a prophage that was present in an ancestral bacterium and subsequently evolved into a GTA that is now widely maintained in extant descendant lineages. Here, we analysed genomic properties that might relate to the conservation of these alphaproteobacterial GTAs. This revealed that the chromosomal locations of the GTA gene clusters are biased. They primarily occur on the leading strand of DNA replication, at large distances from long repetitive elements, and thus are in regions of lower plasticity, and in areas of extreme GC skew, which also accumulate core genes. These extreme GC skew regions arise from the preferential use of codons with an excess of G over C, a distinct phenomenon from the elevated GC content that has previously been found to be associated with GTA genes. The observed properties, along with their high level of conservation, show that GTA genes share multiple features with core genes in the examined lineages of the Alphaproteobacteria.

Analysis of Complete Mitochondrial Genome of Bohadschia argus (Jaeger, 1833) (Aspidochirotida, Holothuriidae)

Bohadschia argu is a kind of sea cucumber with high economic value; it is the only undisputed species in the genus Bohadschia. In this study, the complete mitochondrial genome (mitogenome) of B. argus was acquired through high-throughput sequencing. The mitochondrial genome of B. argus was 15,656 bp in total length and contained a putative control region (CR) and 37 typical genes of animal mitochondrial genomes, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rrnS and rrnL) and 22 transfer RNA genes (tRNA). The sizes of the PCGs ranged from 168 bp to 1833 bp, and all PCGs except nad6 were encoded on the heavy chain (H). Both rrnS and rrnL were also encoded on the H chain. Twenty-two tRNA genes had positive AT skew and GC skew. All tRNAs had a typical cloverleaf secondary structure except for trnI, in which an arm of dihydrouridine was missing. B. argus shared the same gene arrangement order (the echinoderm ground pattern) as other species in Aspidochirotida. Phylogenetic analysis clearly revealed that B. argus belongs as a member of the Holothuriidae, and it is closely related to members of Actinopyga and Holothuria.

A Mitochondrial Genome Phylogeny of Cleridae (Coleoptera, Cleroidea)

The predaceous beetle family Cleridae includes a large and widely distributed rapid radiation, which is vital for the ecosystem. Despite its important role, a number of problems remain to be solved regarding the phylogenetic inter-relationships, the timing of divergence, and the mitochondrial biology. Mitochondrial genomes have been widely used to reconstruct phylogenies of various insect groups, but never introduced to Cleridae until now. Here, we generated 18 mitochondrial genomes to address these issues, which are all novel to the family. In addition to phylogenomic analysis, we have leveraged our new sources to study the mitochondrial biology in terms of nucleotide composition, codon usage and substitutional rate, to understand how these vital cellular components may have contributed to the divergence of the Cleridae. Our results recovered Korynetinae sister to the remaining clerids, and the calde of Clerinae+Hydnocerinae is indicated more related to Tillinae. A time-calibrated phylogeny estimated the earliest divergence time of Cleridae was soon after the origin of the family, not later than 160.18 Mya (95% HPD: 158.18–162.07 Mya) during the mid-Jurassic. This is the first mitochondrial genome-based phylogenetic study of the Cleridae that covers nearly all subfamily members, which provides an alternative evidence for reconstructing the phylogenetic relationships.

Delineation of the Ancestral Tus-Dependent Replication Fork Trap

In Escherichia coli, DNA replication termination is orchestrated by two clusters of Ter sites forming a DNA replication fork trap when bound by Tus proteins. The formation of a ‘locked’ Tus–Ter complex is essential for halting incoming DNA replication forks. However, the absence of replication fork arrest at some Ter sites raised questions about their significance. In this study, we examined the genome-wide distribution of Tus and found that only the six innermost Ter sites (TerA–E and G) were significantly bound by Tus. We also found that a single ectopic insertion of TerB in its non-permissive orientation could not be achieved, advocating against a need for ‘back-up’ Ter sites. Finally, examination of the genomes of a variety of Enterobacterales revealed a new replication fork trap architecture mostly found outside the Enterobacteriaceae family. Taken together, our data enabled the delineation of a narrow ancestral Tus-dependent DNA replication fork trap consisting of only two Ter sites.

Mitochondrial mutations in Caenorhabditis elegans show signatures of oxidative damage and an AT-bias

Rapid mutation rates are typical of mitochondrial genomes (mtDNAs) in animals, but it is not clear why. The difficulty of obtaining measurements of mtDNA mutation that are not biased by natural selection has stymied efforts to distinguish between competing hypotheses about the causes of high mtDNA mutation rates. Several studies which have measured mtDNA mutations in nematodes have yielded small datasets with conflicting conclusions about the relative abundance of different substitution classes (i.e., the mutation spectrum). We therefore leveraged Duplex Sequencing, a high-fidelity DNA sequencing technique, to characterize de novo mtDNA mutations in Caenorhabditis elegans. This approach detected nearly an order of magnitude more mtDNA mutations than documented in any previous nematode mutation study. Despite an existing extreme AT bias in the C. elegans mtDNA (75.6% AT), we found that a significant majority of mutations increase genomic AT content. Compared to some prior studies in nematodes and other animals, the mutation spectrum reported here contains an abundance of CG→AT transversions, supporting the hypothesis that oxidative damage may be a driver of mtDNA mutations in nematodes. Furthermore, we found an excess of G→T and C→T changes on the coding DNA strand relative to the template strand, consistent with increased exposure to oxidative damage. Analysis of the distribution of mutations across the mtDNA revealed significant variation among protein-coding genes and as well as among neighboring nucleotides. This high-resolution view of mitochondrial mutations in C. elegans highlights the value of this system for understanding relationships among oxidative damage, replication error, and mtDNA mutation.

Characterization of the First Complete Mitochondrial Genome of Cyphonocerinae (Coleoptera: Lampyridae) with Implications for Phylogeny and Evolution of Fireflies

Complete mitochondrial genomes are valuable resources for phylogenetics in insects. The Cyphonoceridae represents an important lineage of fireflies. However, no complete mitogenome is available until now. Here, the first complete mitochondrial genome from this subfamily was reported, with Cyphonocerus sanguineus klapperichi as a representative. The mitogenome of C. sanguineus klapperichi was conserved in the structure and comparable to that of others in size and A+T content. Nucleotide composition was A+T-biased, and all genes exhibited a positive AT-skew and negative GC-skew. Two types of tandem repeat sequence units were present in the control region (136 bp × 2; 171 bp × 2 + 9 bp). For reconstruction of Lampyridae's phylogeny, three different datasets were analyzed by both maximum likelihood (ML) and Bayesian inference (BI) methods. As a result, the same topology was produced by both ML analysis of 13 protein-coding genes and 2rRNA and BI analysis of 37 genes. The results indicated that Lampyridae, Lampyrinae, Luciolinae (excluding Emeia) were monophyletic, but Ototretinae was paraphyletic, of which Stenocladius was recovered as the sister taxon to all others, while Drilaster was more closely related to Cyphonocerinae; Phturinae + Emeia were included in a monophyletic clade, which comprised sister groups with Lampyridae. Vesta was deeply rooted in the Luciolinae.

Evolutionary Trajectory of the Replication Mode of Bacterial Replicons

Chromosome replication is an essential process for cell division. The mode of chromosome replication has important impacts on the structure of the chromosome and replication speed.

As typical bacterial replicons, circular chromosomes replicate bidirectionally and circular plasmids replicate either bidirectionally or unidirectionally. Whereas the finding of chromids (plasmid-derived chromosomes) in multiple bacterial lineages provides circumstantial evidence that chromosomes likely evolved from plasmids, all experimentally assayed chromids were shown to use bidirectional replication. Here, we employed a model system, the marine bacterial genus Pseudoalteromonas, members of which consistently carry a chromosome and a chromid. We provide experimental and bioinformatic evidence that while chromids in a few strains replicate bidirectionally, most replicate unidirectionally. This is the first experimental demonstration of the unidirectional replication mode in bacterial chromids. Phylogenomic and comparative genomic analyses showed that the bidirectional replication evolved only once from a unidirectional ancestor and that this transition was associated with insertions of exogenous DNA and relocation of the replication terminus region (ter2) from near the origin site (ori2) to a position roughly opposite it. This process enables a plasmid-derived chromosome to increase its size and expand the bacterium’s metabolic versatility while keeping its replication synchronized with that of the main chromosome. A major implication of our study is that the uni- and bidirectionally replicating chromids may represent two stages on the evolutionary trajectory from unidirectionally replicating plasmids to bidirectionally replicating chromosomes in bacteria. Further bioinformatic analyses predicted unidirectionally replicating chromids in several unrelated bacterial phyla, suggesting that evolution from unidirectionally to bidirectionally replicating replicons occurred multiple times in bacteria.

Potential causes and consequences of rapid mitochondrial genome evolution in thermoacidophilic Galdieria (Rhodophyta)

BACKGROUND

The Cyanidiophyceae is an early-diverged red algal class that thrives in extreme conditions around acidic hot springs. Although this lineage has been highlighted as a model for understanding the biology of extremophilic eukaryotes, little is known about the molecular evolution of their mitochondrial genomes (mitogenomes).

RESULTS

To fill this knowledge gap, we sequenced five mitogenomes from representative clades of Cyanidiophyceae and identified two major groups, here referred to as Galdieria-type (G-type) and Cyanidium-type (C-type). G-type mitogenomes exhibit the following three features: (i) reduction in genome size and gene inventory, (ii) evolution of unique protein properties including charge, hydropathy, stability, amino acid composition, and protein size, and (iii) distinctive GC-content and skewness of nucleotides. Based on GC-skew-associated characteristics, we postulate that unidirectional DNA replication may have resulted in the rapid evolution of G-type mitogenomes.

CONCLUSIONS

The high divergence of G-type mitogenomes was likely driven by natural selection in the multiple extreme environments that Galdieria species inhabit combined with their highly flexible heterotrophic metabolism. We speculate that the interplay between mitogenome divergence and adaptation may help explain the dominance of Galdieria species in diverse extreme habitats.

Evolutionary advantage of anti-parallel strand orientation of duplex DNA

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.

Strand-specific Single-stranded DNA Sequencing (4S-seq) of E. coli genomes

Most bacterial genomes have biased nucleotide composition, and the asymmetry is considered to be caused by a single-stranded DNA (ssDNA) deamination arising from the bacterial replication machinery. In order to evaluate the relationship experimentally, the position and frequency of ssDNA formed during replication must be verified clearly. Although many ssDNA detection technologies exist, almost all methods have been developed for eukaryotic genomes. To apply these to bacterial genomes, which harbor a smaller amount of DNA than those of eukaryotes, more efficient, new methods are required. Therefore, we developed a novel strand-specific ssDNA sequencing method, called 4S-seq, for the bacterial genome. The 4S-seq method enriches ssDNA in the extracted genomic DNA by a dsDNA-specific nuclease and implements a strand-specific library using a biotin label with a customized tag. As a result, the 4S-seq is able to calculate the ssDNA content in each strand (Watson/Crick) at each position of the genome efficiently.