Nucleotide Substitution Rates Research Articles

Comparing the molecular evolution and recombination patterns of predominant PRRSV-2 lineages co-circulating in China.

Porcine reproductive and respiratory syndrome virus (PRRSV) poses widespread epidemics in swine herds, yet the drivers underlying lineage replacements/fitness dynamics remain unclear. To delineate the evolutionary trajectories of PRRSV-2 lineages prevalent in China, we performed a comprehensive longitudinal phylodynamic analysis of 822 viral sequences spanning 1991-2022. The objectives encompassed evaluating lineage dynamics, genetic diversity, recombination patterns and glycosylation profiles. A significant shift in the dominance of PRRSV-2 sub-lineages has been observed over the past 3 decades, transitioning from sub-lineage 8.7 to sub-lineage 1.8, followed by extensive diversification. The analysis revealed discordant recombination patterns between the two dominant viral sub-lineages 1.8 and 8.7, underscoring that modular genetic exchanges contribute significantly to their evolutionary shaping. Additionally, a strong association was found between recombination breakpoint locations and transcriptional regulatory sequences (TRSs). Glycosylation patterns also demonstrated considerable variability across sub-lineages and temporally, providing evidence for immune-driven viral evolution. Furthermore, we quantified different evolutionary rates across sub-lineages, with sub-lineage 1.8 uniquely displaying the highest nucleotide substitution rates. Taken together, these findings provide refined insight into the evolutionary mechanisms underpinning cyclic shifts in dominance among regionally circulating PRRSV sub-lineages.

A phylogenetic method linking nucleotide substitution rates to rates of continuous trait evolution.

Genomes contain conserved non-coding sequences that perform important biological functions, such as gene regulation. We present a phylogenetic method, PhyloAcc-C, that associates nucleotide substitution rates with changes in a continuous trait of interest. The method takes as input a multiple sequence alignment of conserved elements, continuous trait data observed in extant species, and a background phylogeny and substitution process. Gibbs sampling is used to assign rate categories (background, conserved, accelerated) to lineages and explore whether the assigned rate categories are associated with increases or decreases in the rate of trait evolution. We test our method using simulations and then illustrate its application using mammalian body size and lifespan data previously analyzed with respect to protein coding genes. Like other studies, we find processes such as tumor suppression, telomere maintenance, and p53 regulation to be related to changes in longevity and body size. In addition, we also find that skeletal genes, and developmental processes, such as sprouting angiogenesis, are relevant.

Morphological Characteristics and Taxonomical Positions of Annonaceae in Indonesia

Annonaceae is one of the most important and economically valuable plant families in Malesia region. Taxonomical studies could contribute to further utilization and conservation of the studied taxa. Due to their importance, an update in the morphological characters and taxanomical position of Annonaceae is required. This study aims to characterize the taxonomical position and relationships of 24 species belongs to the family Annonaceae found in the Purwodadi Botanic Gardens. In this work, we reported the taxonomical position and the phenetic clustering of 24 Annonaceae species collected from the Purwodadi Botanic Garden, East Java, Indonesia. The taxonomical study was performed based on morphological and phenetic approach where 81 characters (leaves, stem, and flower) were used in the analysis. A dendogram was produced following the scoring of the character using statistical program PAST ver. 3.0. (Hierarchial Cluster Analysis using Jukes Cantor Similarity approach). The dendogram showed clear separation between members of sub-family Annonoideae and Malmeoideae; however Alphonsea javanica and Sageraea lanceolata were not grouped into their respective tribe or sub-family. Further studies are required to evaluate these taxonomical separations. Overall, from the 24 studied species, we found that they grouped into 2 sub-family (Annonoideae and Malmeoideae), 3 tribes (Annoneae, Miliuseae, and Uvarieae), and 13 genus (Alphonsea, Annona, Desmos, Fissistigma, Miliusa, Mitrephora, Orophea, Polyalthia, Popowia, Pseuduvaria, Sageraea, Stelechocarpus, dan Uvaria).

Assembly and comparative analysis of the first complete mitochondrial genome of the invasive water hyacinth, Eichhornia crassipes

Eichhornia crassipes is an aquatic plant in tropical and subtropical regions, renowned for its notorious invasive tendencies. In this study, we assembled the complete mitogenome of E. crassipes into a single circle molecule of 397,361 bp. The mitogenome has 58 unique genes, including 37 protein-coding genes (PCGs), 18 tRNA genes, three rRNA genes, and 47 % GC content. Sixteen (6.93 %) homologous fragments, ranging from 31 bp to 8548 bp, were identified, indicating the transfer of genetic material from chloroplasts to mitochondria. In addition, we detected positive selection in six PCGs (ccmB, ccmC, ccmFC, nad3, nad4 and sdh4), along with the identification of 782 RNA editing sites across 37 mt-PCGs. These findings suggest a potential contribution to the robust adaptation of this invasive plant to the stressful environment. Lastly, we inferred that phylogenetic conflicts of E. crassipes between the plastome and mitogenome may be attributed to the difference in nucleotide substitution rates between the two organelle genomes. In conclusion, our study provided vital genomic resources for further understanding the invasive mechanism of this species and exploring the dynamic evolution of mitogenomes within the monocot clade.

Comparative analysis of the organelle genomes of Aconitum carmichaelii revealed structural and sequence differences and phylogenetic relationships

In this study, we conducted an assembly and analysis of the organelle genomes of Aconitum carmichaelii. Our investigation encompassed the examination of organelle genome structures, gene transfer events, and the environmental selection pressures affecting A. carmichaelii. The results revealed distinct evolutionary patterns in the organelle genomes of A. carmichaelii. Especially, the plastome exhibited a more conserved structure but a higher nucleotide substitution rate (NSR), while the mitogenome displayed a more complex structure with a slower NSR. Through homology analysis, we identified several instances of unidirectional protein-coding genes (PCGs) transferring from the plastome to the mitogenome. However, we did not observe any events which genes moved from the mitogenome to the plastome. Additionally, we observed multiple transposable element (TE) fragments in the organelle genomes, with both organelles showing different preferences for the type of nuclear TE insertion. Divergence time estimation suggested that rapid differentiation occurred in Aconitum species approximately 7.96 million years ago (Mya). This divergence might be associated with the reduction in CO2 levels and the significant uplift of the Qinghai-Tibet Plateau (QTP) during the late Miocene. Selection pressure analysis indicated that the dN/dS values of both organelles were less than 1, suggested that organelle PCGs were subject to purification selection. However, we did not detect any positively selected genes (PSGs) in Subg. Aconitum and Subg. Lycoctonum. This observation further supports the idea that stronger negative selection pressure on organelle genes in Aconitum results in a more conserved amino acid sequence. In conclusion, this study contributes to a deeper understanding of organelle evolution in Aconitum species and provides a foundation for future research on the genetic mechanisms underlying the structure and function of the Aconitum plastome and mitogenome.

The genomic signatures of evolutionary stasis.

Evolutionary stasis characterizes lineages that seldom speciate and show little phenotypic change over long stretches of geological time. Although lineages that appear to exhibit evolutionary stasis are often called living fossils, no single mechanism is thought to be responsible for their slow rates of morphological evolution and low species diversity. Some analyses of molecular evolutionary rates in a handful of living fossil lineages have indicated that these clades exhibit slow rates of genomic change. Here, we investigate mechanisms of evolutionary stasis using a dataset of 1,105 exons for 481 vertebrate species. We demonstrate that two ancient clades of ray-finned fishes classically called living fossils, gars and sturgeons, exhibit the lowest rates of molecular substitution in protein-coding genes among all jawed vertebrates. Comparably low rates of evolution are observed at fourfold degenerate sites in gars and sturgeons, implying a mechanism of stasis decoupled from selection that we speculate is linked to a highly effective DNA repair apparatus. We show that two gar species last sharing common ancestry over 100 million years ago produce morphologically intermediate and fertile hybrids in the wild. This makes gars the oldest naturally hybridizing divergence among eukaryotes and supports a theoretical prediction that slow rates of nucleotide substitution across the genome slow the accumulation of genetic incompatibilities, enabling hybridization across deeply divergent lineages and slowing the rate of speciation over geological timescales. Our results help establish molecular stasis as a barrier to speciation and phenotypic innovation and provide a mechanism to explain the low species diversity in living fossil lineages.

Genotype F of Echovirus 25 with multiple recombination pattern have been persistently and extensively circulating in Chinese mainland

Echovirus 25 (E25), a member of the Enterovirus B (EV-B) species, can cause aseptic meningitis (AM), viral meningitis (VM), and acute flaccid paralysis (AFP). However, systematic studies on the molecular epidemiology of E25, especially those concerning its evolution and recombination, are lacking. In this study, 18 strains of E25, isolated from seven provinces of China between 2009 and 2018, were collected based on the Chinese hand, foot, and mouth disease (HFMD) surveillance network, and 95 sequences downloaded from GenBank were also screened. Based on the phylogenetic analysis of 113 full-length VP1 sequences worldwide, globally occurring E25 strains were classified into 9 genotypes (A–I), and genotype F was the dominant genotype in the Chinese mainland. The average nucleotide substitution rate of E25 was 6.08 × 10–3 substitutions/site/year, and six important transmission routes were identified worldwide. Seventeen recombination patterns were determined, of which genotype F can be divided into 9 recombination patterns. A positive selector site was found in the capsid protein region of genotype F. Recombination analysis and pressure selection analysis for genotype F showed multiple recombination patterns and evolution characteristics, which may be responsible for it being the dominant genotype in the Chinese mainland. This study provides a theoretical basis for the subsequent prevention and control of E25.

Dating in the Dark: Elevated Substitution Rates in Cave Cockroaches (Blattodea: Nocticolidae) Have Negative Impacts on Molecular Date Estimates.

Rates of nucleotide substitution vary substantially across the Tree of Life, with potentially confounding effects on phylogenetic and evolutionary analyses. A large acceleration in mitochondrial substitution rate occurs in the cockroach family Nocticolidae, which predominantly inhabit subterranean environments. To evaluate the impacts of this among-lineage rate heterogeneity on estimates of phylogenetic relationships and evolutionary timescales, we analysed nuclear ultraconserved elements (UCEs) and mitochondrial genomes from nocticolids and other cockroaches. Substitution rates were substantially elevated in nocticolid lineages compared with other cockroaches, especially in mitochondrial protein-coding genes. This disparity in evolutionary rates is likely to have led to different evolutionary relationships being supported by phylogenetic analyses of mitochondrial genomes and UCE loci. Furthermore, Bayesian dating analyses using relaxed-clock models inferred much deeper divergence times compared with a flexible local clock. Our phylogenetic analysis of UCEs, which is the first genome-scale study to include all thirteen major cockroach families, unites Corydiidae and Nocticolidae and places Anaplectidae as the sister lineage to the rest of Blattoidea. We uncover an extraordinary level of genetic divergence in Nocticolidae, including two highly distinct clades that separated ~115 million years ago despite both containing representatives of the genus Nocticola. The results of our study highlight the potential impacts of high among-lineage rate variation on estimates of phylogenetic relationships and evolutionary timescales.

Increasing evolution, prevalence, and outbreaks for rift valley fever virus in the process of breaking geographical barriers

BackgroundRift valley fever (RVF) is listed as one of prioritized diseases by WHO. This study aims to describe RVF virus' landscape distribution globally, and to insight dynamics change of its evolution, prevalence, and outbreaks in the process of breaking geographical barriers. MethodsA systematic literature review and meta-analyses was conducted to estimate RVF prevalence by hosts using a random-effect model. Molecular clock-based phylogenetic analyses were performed to estimate RVF virus nucleotide substitution rates using nucleotide sequences in NCBI database. RVF virus prevalence, nucleotide substitution rates, and outbreaks were compared before and after breaking geographical barriers twice, respectively. ResultsRVF virus was reported from 26 kinds of hosts covering 48 countries from 1930 to 2022. Since RVF broke geographical barriers, (1) nucleotide substitution rates significantly increased after firstly spreading out of Africa in 2000, (2) prevalence in humans significantly increased from 1.92 % (95 % CI: 0.86–3.25 %) to 3.03 % (95 % CI: 2.09–4.12 %) after it broke Sahara Desert geographical barriers in 1977, and to 5.24 % (95 % CI: 3.81–6.82 %) after 2000, (3) RVF outbreaks in humans and the number of wildlife hosts presented increasing trends. RVF virus spillover may exist between bats and humans, and accelerate viral substitution rates in humans. During outbreaks, the RVF virus substitution rates accelerated in humans. 60.00 % RVF outbreaks occurred 0–2 months after floods and (or) heavy rainfall. ConclusionRVF has the increasing risk to cause pandemics, and global collaboration on “One Health” is needed to prevent potential pandemics.

Extraordinary preservation of gene collinearity over three hundred million years revealed in homosporous lycophytes

Homosporous lycophytes (Lycopodiaceae) are a deeply diverged lineage in the plant tree of life, having split from heterosporous lycophytes (Selaginella and Isoetes) ~400 Mya. Compared to the heterosporous lineage, Lycopodiaceae has markedly larger genome sizes and remains the last major plant clade for which no chromosome-level assembly has been available. Here, we present chromosomal genome assemblies for two homosporous lycophyte species, the allotetraploid Huperzia asiatica and the diploid Diphasiastrum complanatum. Remarkably, despite that the two species diverged ~350 Mya, around 30% of the genes are still in syntenic blocks. Furthermore, both genomes had undergone independent whole genome duplications, and the resulting intragenomic syntenies have likewise been preserved relatively well. Such slow genome evolution over deep time is in stark contrast to heterosporous lycophytes and is correlated with a decelerated rate of nucleotide substitution. Together, the genomes of H. asiatica and D. complanatum not only fill a crucial gap in the plant genomic landscape but also highlight a potentially meaningful genomic contrast between homosporous and heterosporous species.

Genomic Analyses Uncover Evolutionary Features of Influenza A/H3N2 Viruses in Yunnan Province, China, from 2017 to 2022.

Influenza A viruses evolve at a high rate of nucleotide substitution, thereby requiring continuous monitoring to determine the efficacy of vaccines and antiviral drugs. In the current study, we performed whole-genome sequencing analyses of 253 influenza A/H3N2 strains from Yunnan Province, China, during 2017-2022. The hemagglutinin (HA) segments of Yunnan A/H3N2 strains isolated during 2017-2018 harbored a high genetic diversity due to heterogeneous distribution across branches. The mutation regularity of the predominant antigenic epitopes of HA segments in Yunnan was inconsistent in different years. Some important functional mutations in gene segments associated with viral adaptation and drug tolerance were revealed. The rapid genomic evolution of Yunnan A/H3N2 strains from 2017 to 2022 mainly concentrated on segments, i.e., matrix protein 2 (M2), non-structural protein 1 (NS1), neuraminidase (NA), NS2, and HA, with a high overall non-synonymous/synonymous substitution ratio (dN/dS). Our results highlighted a decline in vaccine efficacy against the A/H3N2 circulating strains, particularly against the Yunnan 2021-2022 A/H3N2 strains. These findings aid our understanding of evolutionary characteristics and epidemiological monitoring of the A/H3N2 viruses and provide in-depth insights into the protective efficacy of influenza vaccines.

Genome reduction occurred in early Prochlorococcus with an unusually low effective population size.

In the oligotrophic sunlit ocean, the most abundant free-living planktonic bacterial lineages evolve convergently through genome reduction. The cyanobacterium Prochlorococcus responsible for 10% global oxygen production is a prominent example. The dominant theory known as "genome streamlining" posits that they have extremely large effective population sizes (Ne) such that selection for metabolic efficiency acts to drive genome reduction. Because genome reduction largely took place anciently, this theory builds on the assumption that their ancestors' Ne was similarly large. Constraining Ne for ancient ancestors is challenging because experimental measurements of extinct organisms are impossible and alternatively reconstructing ancestral Ne with phylogenetic models gives large uncertainties. Here, we develop a new strategy that leverages agent-based modeling to simulate the changes in the genome-wide ratio of radical to conservative nonsynonymous nucleotide substitution rate (dR/dC) in a possible range of Ne in ancestral populations. This proxy shows expected increases with decreases of Ne only when Ne falls to about 10k - 100k or lower, magnitudes characteristic of Ne of obligate endosymbiont species where drift drives genome reduction. Our simulations therefore strongly support a scenario where the primary force of Prochlorococcus genome reduction is drift rather than selection.

Molecular characterization of mungbean yellow mosaic India virus infecting Vigna radiata in Oman

AbstractThe full‐genome sequences of a bipartite begomovirus were identified in a mungbean (Vigna radiata) crop displaying symptoms such as yellow mosaics and stunting in the Al‐Batinah North region (23.6867° N 57.9058° E), Oman. Initial detection and confirmation of the virus genome were carried out through polymerase chain reaction, followed by obtaining complete genomes using rolling circle amplification. Bioinformatics analysis on both genomic components, DNA‐A and DNA‐B, revealed over 99% nucleotide sequence identity to the mungbean yellow mosaic India virus (MYMIV), previously known to infect tomato plants. Pairwise sequence analysis using the STD tool and subsequent phylogenetic analysis unveiled that the DNA‐A of the MYMIV isolate formed a cluster with closely related DNA‐A sequences of MYMIV from GenBank, while the DNA‐B clustered with the corresponding DNA‐B of the MYMIV isolate. Additionally, no evidence of recombination events was observed in the recombination analysis. Similarly, the nucleotide substitution rates between the DNA‐A and DNA‐B segments of MYMIV did not show significant differences. This finding represents the first documented instance of bipartite MYMIV infecting V. radiata in Oman.

Phylogenetic Determination of Chenopodium quinoaBased on the Chloroplast Genes rbcL and matK

Chenopodium quinoa is an Andean species of great interest because of its excellent nutritional quality and great adaptability to different environmental conditions. In addition, the high phenotypic diversity has caused difficulties in the correct taxonomic identification, and there are few studies on the phylogenetic relationships of quinoa in Colombia. Therefore, the objective of this research was to determine the phylogenetic relationships of quinoa with the matK and rcbL chloroplastid genes to characterize the genetic diversity in Colombian quinoa. Evolutionary analyses were performed using nucleotide substitution rates, pattern, base composition, and phylogeny construction. The rbcL gene presented approximately 1344 bp, and matK had 646 bp, which were translated into 434 and 215 amino acids, respectively. The nucleotide composition of the genes showed high percentages of similarity and identity with the Chenopodium quinoa sequences registered in GenBank and BOLD. Similar phylogenetic trees were obtained with the rbcL and matK genes, and both concatenated sequences grouped the accessions into clades. The results showed that Colombian quinoa has low rates of genetic differentiation that may be due to the domestication processes of the species, the lack of certified seeds, and the constant exchange of seeds between farmers in the principal producing areas of the Andean region.

Genetic Diversity, Evolutionary Dynamics, and Ongoing Spread of Pedilanthus Leaf Curl Virus.

Pedilanthus leaf curl virus (PeLCV) is a monopartite begomovirus (family Geminiviridae) discovered just a few decades ago. Since then, it has become a widely encountered virus, with reports from ca. 25 plant species across Pakistan and India, indicative of its notable evolutionary success. Viruses mutate at such a swift rate that their ecological and evolutionary behaviors are inextricably linked, and all of these behaviors are imprinted on their genomes as genetic diversity. So, all these imprints can be mapped by computational methods. This study was designed to map the sequence variation dynamics, genetic heterogeneity, regional diversity, phylogeny, and recombination events imprinted on the PeLCV genome. Phylogenetic and network analysis grouped the full-length genome sequences of 52 PeLCV isolates into 7 major clades, displaying some regional delineation but lacking host-specific demarcation. The progenitor of PeLCV was found to have originated in Multan, Pakistan, in 1977, from where it spread concurrently to India and various regions of Pakistan. A high proportion of recombination events, distributed unevenly throughout the genome and involving both inter- and intraspecies recombinants, were inferred. The findings of this study highlight that the PeLCV population is expanding under a high degree of genetic diversity (π = 0.073%), a high rate of mean nucleotide substitution (1.54 × 10-3), demographic selection, and a high rate of recombination. This sets PeLCV apart as a distinctive begomovirus among other begomoviruses. These factors could further exacerbate the PeLCV divergence and adaptation to new hosts. The insights of this study that pinpoint the emergence of PeLCV are outlined.

Buried treasure in a public repository: Mining mitochondrial genes of 32 annelid species from sequence reads deposited in the Sequence Read Archive (SRA).

The mitochondrial genomes (mitogenomes) of metazoans generally include the same set of protein-coding genes, which ensures the homology of mitochondrial genes between species. The mitochondrial genes are often used as reference data for species identification based on genetic data (DNA barcoding). The need for such reference data has been increasing due to the application of environmental DNA (eDNA) analysis for environmental assessments. Recently, the number of publicly available sequence reads obtained with next-generation sequencing (NGS) has been increasing in the public database (the NCBI Sequence Read Archive, SRA). Such freely available NGS reads would be promising sources for assembling mitochondrial protein-coding genes (mPCGs) of organisms whose mitochondrial genes are not available in GenBank. The present study aimed to assemble annelid mPCGs from raw data deposited in the SRA. The recent progress in the classification of Annelida was briefly introduced. In the present study, the mPCGs of 32 annelid species of 19 families in clitellates and allies in Sedentaria (echiurans and polychaetes) were newly assembled from the reads deposited in the SRA. Assembly was performed with a recently published pipeline mitoRNA, which includes cycles of Bowtie2 mapping and Trinity assembly. Assembled mPCGs were deposited in GenBank as Third Party Data (TPA) data. A phylogenetic tree was reconstructed with maximum likelihood (ML) analysis, together with other mPCGs deposited in GenBank. mPCG assembly was largely successful except for Travisia forbesii; only four genes were detected from the assembled contigs of the species probably due to the reads targeting its parasite. Most genes were largely successfully obtained, whereas atp8, nad2, and nad4l were only successful in 22-24 species. The high nucleotide substitution rates of these genes might be relevant to the failure in the assembly although nad6, which showed a similarly high substitution rate, was successfully assembled. Although the phylogenetic positions of several lineages were not resolved in the present study, the phylogenetic relationships of some polychaetes and leeches that were not inferred by transcriptomes were well resolved probably due to a more dense taxon sampling than previous phylogenetic analyses based on transcriptomes. Although NGS data are generally better sources for resolving phylogenetic relationships of both higher and lower classifications, there are ensuring needs for specific loci of the mitochondrial genes for analyses that do not require high resolutions, such as DNA barcoding, eDNA, and phylogenetic analysis among lower taxa. Assembly from publicly available NGS reads would help design specific primers for the mitochondrial gene sequences of species, whose mitochondrial genes are hard to amplify by Sanger sequencing using universal primers.

Influence of the Rate of Changes in the COX1 Gene on Body Size and Sexual Selection in Carp Hybridization

The influence of mtDNA cytochrome c-oxidase I gene fragment variability on body length was studied in twelve species of cyprinids, which may have hybrids with Rutilus rutilus L. and Abramis brama L., and in reciprocal hybrids (RA, AR) and alloplasmatic backcrosses (ARR, RAA) of roach (R) and bream (A). It has been established that the rate of nucleotide substitutions in COX1 is negatively related not only to body size but also to fish life span, which differentiates them into two groups: group I – species with a high rate of COX1 changes and a relatively small body size and group II – species with low sequence variability and relatively large body size. The boundary for the distinguished groups runs between species the same genus Leuciscus leuciscus and L. idus: with a twofold decrease in the rate of substitutions in ide, a twofold increase in body size and lifespan occurs, which indicates a decrease in the rate of cellular respiration and free radical leak, and the exact mitonuclear match respiratory complexes. Presumably, the decrease in the rate of COX1 changes in species of group II and in bleak Alburnus alburnus is associated with an increase in the size of genome, which provides additional protection of genes from chemical mutagens and, regardless of body size, reduces the rate of aerobic metabolism. It has been experimentally shown that mtDNA affects body length. When bream mtDNA is included in the roach nuclear genome, ARR backcrosses have the body length of a bream and high viability, while RAA backcrosses with roach mtDNA and the bream nuclear genome inherit the roach body length and reduce viability. Species of group II are not able to effectively use the highly polymorphic mtDNA of species of group I, which is also manifested by a violation of the inheritance of a longer bream body length in RA hybrids and leads to reproductive isolation. Group I species, such as Rutilus rutilus, can include mtDNA of both groups in their genome, which underlies sexual selection in hybridization. Accordingly, sexual size dimorphism has a genetic origin, and body size for a potential partner can be a signal for determining the mitonuclear compatibility of genomes in respiratory complexes.

Genomic estimates of mutation and substitution rates contradict the evolutionary speed hypothesis of the latitudinal diversity gradient.

The latitudinal diversity gradient (LDG) refers to a decrease in biodiversity from the equator to the poles. The evolutionary speed hypothesis, backed by the metabolic theory of ecology, asserts that nucleotide mutation and substitution rates per site per year are higher and thereby speciation rates are higher at higher temperatures, generating the LDG. However, prior empirical investigations of the relationship between the temperature and mutation or substitution rate were based on a few genes and the results were mixed. We here revisit this relationship using genomic data. No significant correlation between the temperature and mutation rate is found in 13 prokaryotes or in 107 eukaryotes. An analysis of 234 diverse trios of bacterial taxa indicates that the synonymous substitution rate is not significantly associated with the growth temperature. The same data, however, reveal a significant negative association between the nonsynonymous substitution rate and temperature, which is explainable by a larger fraction of detrimental nonsynonymous mutations at higher temperatures due to a stronger demand for protein stability. We conclude that the evolutionary speed hypothesis of the LDG is unsupported by genomic data and advise that future mechanistic studies of the LDG should focus on other hypotheses.

Genomic diversity and natural recombination of equid gammaherpesvirus 5 isolates

BackgroundEquid gammaherpesvirus 5 (EHV5) is closely related to equid gammaherpesvirus 2 (EHV2). Detection of EHV5 is frequent in horse populations worldwide, but it is often without a clear and significant clinical impact. Infection in horses can often present as subclinical disease; however, it has been associated with respiratory disease, including equine multinodular pulmonary fibrosis (EMPF). Genetic heterogeneity within small regions of the EHV5 glycoprotein B (gB) sequences have been reported and multiple genotypes of this virus have been identified within individual horses, but full genome sequence data for these viruses is limited. The primary focus of this study was to assess the genomic diversity and natural recombination among EHV5 isolates. ResultsThe genome size of EHV5 prototype strain and the five EHV5 isolates cultured for this study, including four isolates from the same horse, ranged from 181,929 to 183,428 base pairs (bp), with the sizes of terminal repeat regions varying from 0 to 10 bp. The nucleotide sequence identity between the six EHV5 genomes ranged from 95.5 to 99.1%, and the estimated average nucleotide diversity between isolates was 1%. Individual genes displayed varying levels of nucleotide diversity that ranged from 0 to 19%. The analysis of nonsynonymous substitution (Ka > 0.025) revealed high diversity in eight genes. Genome analysis using RDP4 and SplitsTree programs detected evidence of past recombination events between EHV5 isolates. ConclusionGenomic diversity and recombination hotspots were identified among EHV5 strains. Recombination can drive genetic diversity, particularly in viruses that have a low rate of nucleotide substitutions. Therefore, the results from this study suggest that recombination is an important contributing factor to EHV5 genomic diversity. The findings from this study provide additional insights into the genetic heterogeneity of the EHV5 genome.

Natural recombination of the torque teno canis virus within the ORF1, -2, and -3 genes and the untranslated region

The torque teno canis virus (TTCaV) was first reported in 2001 and it shares similarities with the known Torque teno virus (TTV) in terms of genomic organization and putative transcriptional features. It is a single-stranded DNA virus characterized by high rates of recombination and nucleotide substitution, like RNA viruses. Studies reported recombination events in torque teno virus; however, there is limited reporting of TTCaV reorganization events. This study screened fecal samples from domestic dogs in Henan Province. There was a positivity rate of 16.5% (19/115) for TTCaV. Four nearly complete TTCaV genomes, namely Canine/HeNan/4, 5, 6, and 13/2019, were obtained from the 19 positive fecal samples, whose genome sequence was obtained using gap-filling PCR. Sequence analysis revealed two unique amino acid mutation sites in the TTCaV strains, K278Q (compared with the first isolate Cf-TTV10 in Japan) and V/L268I (compared with the TTCaV strain from southern China). Subsequently, 17 near full-length TTCaV genome sequences were subjected to phylogenetic and recombination detection program analyzes. We obtained evidence supporting recombination events in the Chinese TTCaV strains. These findings suggest that mutation and recombination occurred in the three individual gene segments (ORF1, ORF2, ORF3) and the untranslated region, an area of major recombination in the Chinese TTCaV strain GX265 genome. Interestingly, the TTCaV strain (Canine/HeNan/6/2019) was a major parent involved in the genetic recombination of the GX265 strain. This study provides insights into the genetic variability and evolution of TTCaV.