Nucleotide Substitution Rates Research Articles

Genomic surveillance reveals low-level circulation of two subtypes of genogroup C coxsackievirus A10 in Nanchang, Jiangxi Province, China, 2015-2023.

In recent years, coxsackievirus (CV) A10 has been associated with increasing sporadic hand, foot, and mouth disease (HFMD) cases and outbreaks globally. In addition to mild symptoms such as pharyngitis and herpangina, CVA10-related complications or even fatality can occur. Currently, systematic phylogenetic studies of CVA10 are limited. In this study, we first explored the epidemiological and genetic characteristics of CVA10 in Nanchang, an inland southeastern city of China, based on the HFMD surveillance network from 2015-2023. Among 3429 enterovirus-positive cases, 110 (3.04%) were associated with CVA10, with a male-to-female ratio of 1.62. The median age of the CVA10 patients was 2.3 years (interquartile range, IQR 1.0-4.0), with 94.55% (104/110) of the patients aged less than 5 years. Phylogenetic analyses using the full-length VP1, 5'UTR, P1, P2, P3 sequences and near full-length genomes indicated that CVA10 strains (n = 57) isolated in Nanchang belonged to genogroup C; two strains identified in 2017 belonged to C1 subtypes clustered with strains from Vietnam, Madagascar, France and Spain; and the others belonged to C2 subtypes interdigitating with CVA10 isolates from mainland China, the United States and Australia. Through extensive analysis, we identified a rare F168Y mutation in epitope 4 of VP1 in a Madagascar strain of genogroup F and a Chinese strain of genogroup C. Based on Bayesian evolutionary analyses, the average nucleotide substitution rate for the VP1 gene of CV10 strains was 3.07×10-3 substitutions/site/year. The most recent common ancestor (tMRCA) of genogroup C was dated 1990.84, and the tMRCA of CVA10 strains from Nanchang was dated approximately 2003.16, similar to strains circulating in other regions of China, suggesting that the viruses were likely introduced and cryptically circulated in China before the establishment of the HFMD surveillance network. Recombination analysis indicated intertypic recombination of the Nanchang strain with the genogroup G strain in the 3D region. Given the shifting dominance of viral genotypes and frequent recombination events, the existing surveillance system needs to be regulated to enhance genomic surveillance efforts on a more diverse spectrum of genotypes in the future.

Genome-wide identification of nuclear factor -Y (NF–Y) transcription factor family in finger millet reveals structural and functional diversity

The Nuclear Factor Y (NF–Y) is one of the widely explored transcription factors (TFs) family for its potential role in regulating molecular mechanisms related to stress response and developmental processes.Finger millet (Eleusine coracana (L.) Gaertn) is a hardy and stress-tolerant crop where partial efforts have been made to characterize a few transcription factors. However, the NF–Y TF is still poorly explored and not well documented. The present study aims to identify and characterize NF–Y genes of finger millet using a bioinformatics approach. Genome mining revealed 57 EcNF-Y (Eleusine coracana Nuclear Factor-Y) genes in finger millet, comprising 18 NF-YA, 23 NF-YB, and 16 NF-YC genes. The gene organization, conserved motif, cis-regulatory elements, miRNA target sites, and three-dimensional structures of these NF-Ys were analyzed. The nucleotide substitution rate and gene duplication analysis showed the presence of 7 EcNF-YA, 10 EcNF-YB, and 8 EcNF-YC paralogous genes and revealed the possibilities of synonymous substitution and stabilizing selection during evolution. The role of NF-Ys of finger millet in abiotic stress tolerance was evident by the presence of relevant cis-elements such as ABRE (abscisic acid-responsive elements), DRE (dehydration-responsive element), MYB (myeloblastosis) or MYC (myelocytomatosis). Twenty-three isoforms of miR169, mainly targeting a single NF–Y gene, i.e., the EcNF-YA13 gene, were observed. This interaction could be targeted for finger millet improvement against Magnaporthe oryzae (blast fungus). Therefore, by this study, the putative functions related to biotic and abiotic stress tolerance for many of the EcNF-Y genes could be explored in finger millet.

The GC-content at the 5′ ends of human protein-coding genes is undergoing mutational decay

BackgroundIn vertebrates, most protein-coding genes have a peak of GC-content near their 5′ transcriptional start site (TSS). This feature promotes both the efficient nuclear export and translation of mRNAs. Despite the importance of GC-content for RNA metabolism, its general features, origin, and maintenance remain mysterious. We investigate the evolutionary forces shaping GC-content at the transcriptional start site (TSS) of genes through both comparative genomic analysis of nucleotide substitution rates between different species and by examining human de novo mutations.ResultsOur data suggests that GC-peaks at TSSs were present in the last common ancestor of amniotes, and likely that of vertebrates. We observe that in apes and rodents, where recombination is directed away from TSSs by PRDM9, GC-content at the 5′ end of protein-coding gene is currently undergoing mutational decay. In canids, which lack PRDM9 and perform recombination at TSSs, GC-content at the 5′ end of protein-coding is increasing. We show that these patterns extend into the 5′ end of the open reading frame, thus impacting synonymous codon position choices.ConclusionsOur results indicate that the dynamics of this GC-peak in amniotes is largely shaped by historic patterns of recombination. Since decay of GC-content towards the mutation rate equilibrium is the default state for non-functional DNA, the observed decrease in GC-content at TSSs in apes and rodents indicates that the GC-peak is not being maintained by selection on most protein-coding genes in those species.

Seek, and you will find: Cryptic diversity of the cardiopulmonary nematode Angiostrongylus vasorum in the Americas

Angiostrongylus vasorum is a metastrongylid parasite infecting wild canids and domestic dogs. Its patchy distribution, high pathogenicity and taxonomical classification makes the evolutionary history of A. vasorum intriguing and important to study. First larval stages of A. vasorum were recovered from feces of two grey foxes, Urocyon cinereoargenteus, from Costa Rica. Sequencing and phylogenetic and haplotypic analyses of the ITS2, 18S and cytochrome oxidase subunit 1 (cox1) fragments were performed. Then p- and Nei´s genetic distance, nucleotide substitution rates and species delimitation analyses were conducted with cox1 data of the specimens collected herein and other Angiostrongylus spp. Cophylogenetic congruence and coevolutionary events of Angiostrongylus spp. and their hosts were evaluated using patristic and phenetic distances and maximum parsimony reconciliations. Specimens from Costa Rica clustered in a separate branch from European and Brazilian A. vasorum sequences in the phylogenetic and haplotype network analyses using the ITS2 and cox1 data. In addition, cox1 p-distance of the sequences derived from Costa Rica were up to 8.6 % different to the ones from Europe and Brazil, a finding mirrored in Nei´s genetic distance PCoA. Species delimitation analysis supported a separate group with the sequences from Costa Rica, suggesting that these worms may represent cryptic variants of A. vasorum, a new undescribed taxon or Angiocaulus raillieti, a synonym species of A. vasorum described in Brazil. Moreover, nucleotide substitution rates in A. vasorum were up to six times higher than in the congener Angiostrongylus cantonensis. This finding and the long time elapsed since the last common ancestor between both species may explain the larger diversity in A. vasorum. Finally, cophylogenetic congruence was observed between Angiostrongylus spp. and their hosts, with cospeciation events occurring at deeper taxonomic branching of host order. Altogether, our data suggest that the diversity of the genus Angiostrongylus is larger than expected, since additional species may be circulating in wild canids from the Americas.

Accumulation of Large Lineage-Specific Repeats Coincides with Sequence Acceleration and Structural Rearrangement in Plantago Plastomes.

Repeats can mediate rearrangements and recombination in plant mitochondrial genomes and plastid genomes. While repeat accumulations are linked to heightened evolutionary rates and complex structures in specific lineages, debates persist regarding the extent of their influence on sequence and structural evolution. In this study, 75 Plantago plastomes were analyzed to investigate the relationships between repeats, nucleotide substitution rates, and structural variations. Extensive repeat accumulations were associated with significant rearrangements and inversions in the large inverted repeats (IRs), suggesting that repeats contribute to rearrangement hotspots. Repeats caused infrequent recombination that potentially led to substoichiometric shifting, supported by long-read sequencing. Repeats were implicated in elevating evolutionary rates by facilitating localized hypermutation, likely through DNA damage and repair processes. This study also observed a decrease in nucleotide substitution rates for loci translocating into IRs, supporting the role of biased gene conversion in maintaining lower substitution rates. Combined with known parallel changes in mitogenomes, it is proposed that potential dysfunction in nuclear-encoded genes associated with DNA replication, recombination, and repair may drive the evolution of Plantago organellar genomes. These findings contribute to understanding how repeats impact organellar evolution and stability, particularly in rapidly evolving plant lineages.

Draft Genome of Akame (Lates Japonicus) Reveals Possible Genetic Mechanisms for Long-Term Persistence and Adaptive Evolution with Low Genetic Diversity.

It is known that some endangered species have persisted for thousands of years despite their very small effective population sizes and low levels of genetic polymorphisms. To understand the genetic mechanisms of long-term persistence in threatened species, we determined the whole genome sequences of akame (Lates japonicus), which has survived for a long time with extremely low genetic variations. Genome-wide heterozygosity in akame was estimated to be 3.3 to 3.4 × 10-4/bp, one of the smallest values in teleost fishes. Analysis of demographic history revealed that the effective population size in akame was around 1,000 from 30,000 years ago to the recent past. The relatively high ratio of nonsynonymous to synonymous heterozygosity in akame indicated an increased genetic load. However, a detailed analysis of genetic diversity in the akame genome revealed that multiple genomic regions, including genes involved in immunity, synaptic development, and olfactory sensory systems, have retained relatively high nucleotide polymorphisms. This implies that the akame genome has preserved the functional genetic variations by balancing selection, to avoid a reduction in viability and loss of adaptive potential. Analysis of synonymous and nonsynonymous nucleotide substitution rates has detected signs of positive selection in many akame genes, suggesting adaptive evolution to temperate waters after the speciation of akame and its close relative, barramundi (Lates calcarifer). Our results indicate that the functional genetic diversity likely contributed to the long-term persistence of this species by avoiding the harmful effects of the population size reduction.

Purifying selection drove the adaptation of mitochondrial genes along with correlation of gene rearrangements and evolutionary rates in two subfamilies of Whitefly (Insecta: Hemiptera).

Insectmitochondrial genomes(mitogenomes) are usually represented by a conserved gene order. Whiteflies exhibit gene rearrangement in their mitogenomes; however, understanding how nucleotide substitution rates shape gene rearrangement in whiteflies is unclear due to the limited number of mitogenomes. Additionally, the mechanisms by which selection pressure drives adaptations in mitochondrial genes in the two subfamilies of whiteflies are not yet known. Here, we analyzed 18 whitefly mitogenomes, including one newly generated mitogenome, to compare nucleotide substitution rates, selection pressure, and gene arrangements. The newly generated mitogenome is reported along with reannotation of Pealius mori and comparisons to other whitefly mitogenomes. Comparative studies on nucleotide composition of 18 whiteflies revealed the positive GC skewness, confirming the reversal of strand asymmetry. We found 11 rearranged gene orders within two subfamilies of whiteflies with 8-18 breakpoints of gene rearrangements. Members of the subfamily Aleyrodinae exhibit more complex pathways in the evolution of gene order as compared to the subfamily Aleurodicinae. Our findings also revealed that the increase or reduction of nucleotide substitution rates does not have an impact on any of the gene rearrangement scenarios depicting neutral correlation. Selection pressure analysis revealed that the mitogenomes from members of both the subfamilies Aleurodicinaeand Aleyrodinae are characterized by intense purifying selection pressure.

Mitochondrial genome complexity in Stemona sessilifolia: nanopore sequencing reveals chloroplast gene transfer and DNA rearrangements.

Mitochondria are semi-autonomous organelles in eukaryotic cells with their own genome. Plant mitogenomes differ from animal mitogenomes in size, structure, and repetitive DNA sequences. Despite larger sizes, plant mitogenomes do not have significantly more genes. They exhibit diverse structures due to variations in size, repetitive DNA, recombination frequencies, low gene densities, and reduced nucleotide substitution rates. In this study, we analyzed the mitochondrial genome of Stemona sessilifolia using Nanopore and Illumina sequencing. De-novo assembly and annotation were conducted using Unicycler, Geseq, tRNAscan-SE and BLASTN, followed by codon usage, repeat sequence, RNA-editing, synteny, and phylogenetic analyses. S. sessilifolia's mitogenome consisted of one linear contig and six circular contigs totaling 724,751bp. It had 39 protein-coding genes, 27 tRNA genes, and 3 rRNA genes. Transfer of chloroplast sequences accounted for 13.14% of the mitogenome. Various analyses provided insights into genetic characteristics, evolutionary dynamics, and phylogenetic placement. Further investigations can explore transferred genes' functions and RNA-editing's role in mitochondrial gene expression in S. sessilifolia.

Practical guidance and workflows for identifying fast evolving non-coding genomic elements using PhyloAcc.

Comparative genomics provides ample ways to study genome evolution and its relationship to phenotypic traits. By developing and testing alternate models of evolution throughout a phylogeny, one can estimate rates of molecular evolution along different lineages in a phylogeny and link these rates with observations in extant species, such as convergent phenotypes. Pipelines for such work can help identify when and where genomic changes may be associated with, or possibly influence, phenotypic traits. We recently developed a set of models called PhyloAcc, using a Bayesian framework to estimate rates of nucleotide substitution on different branches a phylogenetic tree and evaluate their association with pre-defined or estimated phenotypic traits PhyloAcc-ST and PhyloAcc-GT both allow users to define a priori a set of target lineages and then compare different models to identify loci accelerating in one or more target lineages. Whereas ST considers only one species tree across all input loci, GT considers alternate topologies for every locus. PhyloAcc-C simultaneously models molecular rates and rates of continuous trait evolution,allowing the user to ask whether the two are associated. Here we describe these models and provide tips and workflows on how to prepare the input data and run PhyloAcc.

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 analyzed 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 13 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.