Positive Selection Sites Research Articles

Comparative chloroplast genomics reveals the phylogeny and the adaptive evolution of Begonia in China

BackgroundThe Begonia species are common shade plants that are mostly found in southwest China. They have not been well studied despite their medicinal and decorative uses because gene penetration, decreased adaptability, and restricted availability are all caused by frequent interspecific hybridization.ResultTo understand the patterns of mutation in the chloroplast genomes of different species of Begonia, as well as their evolutionary relationships, we collected seven Begonia species in China and sequenced their chloroplast genomes. Begonia species exhibit a quadripartite structure of chloroplast genomes (157,634 − 169,694 bp), consisting of two pairs of inverted repeats (IR: 26,529 − 37,674 bp), a large single copy (LSC: 75,477 − 86,500 bp), and a small single copy (SSC: 17,861 − 18,367 bp). 128–143 genes (comprising 82–93 protein-coding genes, 8 ribosomal RNAs, and 36–43 transfer RNAs) are found in the chloroplast genomes. Based on comparative analyses, this taxon has a relatively similar genome structure. A total of six substantially divergent DNA regions (trnT-UGU-trnL-UAA, atpF-atpH, ycf4-cemA, psbC-trnS-UGA, rpl32-trnL-UAG, and ccsA-ndhD) are found in the seventeen chloroplast genomes. These regions are suitable for species identification and phylogeographic analysis. Phylogenetic analysis shows that Begonia species that were suited to comparable environments grouped in a small clade and that all Begonia species formed one big clade in the phylogenetic tree, supporting the genus’ monophyly. In addition, positive selection sites were discovered in eight genes (rpoC1, rpoB, psbE, psbK, petA, rps12, rpl2, and rpl22), the majority of which are involved in protein production and photosynthesis.ConclusionUsing these genome resources, we can resolve deep-level phylogenetic relationships between Begonia species and their families, leading to a better understanding of evolutionary processes. In addition to enhancing species identification and phylogenetic resolution, these results demonstrate the utility of complete chloroplast genomes in phylogenetically and taxonomically challenging plant groupings.

A Comprehensive Analysis of Chloroplast Genome Provides New Insights into the Evolution of the Genus Chrysosplenium.

Chrysosplenium, a perennial herb in the family Saxifragaceae, prefers to grow in low light and moist environments and is divided into two sections of Alternifolia and Oppositifolia based on phyllotaxy. Although there has been some progress in the phylogeny of Chrysosplenium over the years, the phylogenetic position of some species is still controversial. In this study, we assembled chloroplast genomes (cp genomes) of 34 Chrysosplenium species and performed comparative genomic and phylogenetic analyses in combination with other cp genomes of previously known Chrysosplenium species, for a total of 44 Chrysosplenium species. The comparative analyses revealed that cp genomes of Chrysosplenium species were more conserved in terms of genome structure, gene content and arrangement, SSRs, and codon preference, but differ in genome size and SC/IR boundaries. Phylogenetic analysis showed that cp genomes effectively improved the phylogenetic support and resolution of Chrysosplenium species and strongly supported Chrysosplenium species as a monophyletic taxon and divided into three branches. The results also showed that the sections of Alternifolia and Oppositifolia were not monophyletic with each other, and that C. microspermum was not clustered with other Chrysosplenium species with alternate leaves, but with C. sedakowii into separate branches. In addition, we identified 10 mutational hotspot regions that could serve as potential DNA barcodes for Chrysosplenium species identification. In contrast to Peltoboykinia, the clpP and ycf2 genes of Chrysosplenium were subjected to positive selection and had multiple significant positive selection sites. We further detected a significant positive selection site on the petG gene between the two sections of Chrysosplenium. These evolutionary characteristics may be related to the growth environment of Chrysosplenium species. This study enriches the cp genomes of Chrysosplenium species and provides a reference for future studies on its evolution and origin.

Molecular Mechanism of Response and Adaptation of Antioxidant Enzyme System to Salt Stress in Leaves of Gymnocarpos przewalskii.

The antioxidant enzyme system is the main defense system responsible for maintaining cellular reactive oxygen species (ROS) homeostasis and normal plant growth and development after saline stress. In this study, we identified and characterized the members of the SOD, APX and CAT gene families of the antioxidant enzyme system in Gymnocarpos przewalskii, using plant physiology and molecular biology methods, and analyzed the pattern of enzyme activity in response to NaCl stress. It was found that seven, six and two genes of SOD, APX and CAT gene families, respectively, were expressed in the leaf tissue of G. przewalskii, in which most of the genes were significantly upregulated under NaCl stress, and the enzymatic activities were in accordance with the gene expression. Three positive selection sites in the GpCAT1 gene can increase the hydrophilicity of the GpCAT1 protein, increase the volume of the active site and increase the affinity for H2O2, thus improving the catalytic efficiency of GpCAT1. The results of the present study provide new insights for further investigations of the evolution and function of the SOD, APX and CAT gene families in G. przewalskii and their essential roles under salt stress, and the findings will be useful for revealing the molecular mechanism of salt tolerance and breeding of salt-tolerant plants.

The epidemiological and genetic characteristics of human parvovirus B19 in patients with febrile rash illnesses in China

To understand the epidemiological and genetic characteristics of B19V, a multiple-province surveillance of patients with febrile rash illnesses (FRIs) were conducted in China during 2009 ~ 2021. The clinical specimens of 3,820 FRI patients were collected and tested for B19V DNA. A total of 99 (2.59%) patients were positive for B19V, and 49 (49.49%) were children under 5 years old. B19V infections occurred throughout the year without obvious seasonal pattern. Ten NS1-VP1u sequences and seven genome sequences were obtained in this study, identified as subgenotype 1a. Combined with the globally representative genome sequences, no temporal and geographic clustering trends of B19V were observed, and there was no significant correlation between B19V sequences and clinical manifestations. The evolutionary rate of the B19V genome was 2.30 × 10–4 substitutions/site/year. The number of negative selection sites was higher than that of positive selection sites. It was the first to comprehensively describe the prevalence patterns and evolutionary characteristics of B19V in FRI patients in China. B19V played the role in FRI patients. Children under 5 years old were the main population of B19V infection. Subgenotype 1a was prevalent in FRI patients in China. B19V showed a high mutation rate, while negative selection acted on the genome.

Complete genome analysis of echovirus 30 strains isolated from hand-foot-and-mouth disease in Yunnan province, China

BackgroundEchovirus 30 is prone to cause hand-foot-and-mouth disease in infants and children. However, molecular epidemiologic information on the spread of E30 in southwestern China remains limited. In this study, we determined and analyzed the whole genomic sequences of E30 strains isolated from the stools of patients with hand-foot-and-mouth disease in Yunnan Province, China, in 2019.MethodsE30 isolates were obtained from fecal samples of HFMD patients. The whole genomes were sequenced by segmented PCR and analyzed for phylogeny, mutation and recombination. MEGA and DNAStar were used to align the present isolates with the reference strains. The VP1 sequence of the isolates were analyzed for selection pressure using datamonkey server.ResultsThe complete genome sequences of four E30 were obtained from this virus isolation. Significant homologous recombination signals in the P2-3’UTR region were found in all four isolates with other serotypes. Phylogenetic analysis showed that the four E30 isolates belonged to lineage H. Comparison of the VP1 sequences of these four isolates with other E30 reference strains using three selection pressure analysis models FUBAR, FEL, and MEME, revealed a positive selection site at 133rd position.ConclusionsThis study extends the whole genome sequence of E30 in GenBank, in which mutations and recombinations have driven the evolution of E30 and further improved and enriched the genetic characteristics of E30, providing fundamental data for the prevention and control of diseases caused by E30. Furthermore, we demonstrated the value of continuous and extensive surveillance of enterovirus serotypes other than the major HFMD-causing viruses.

Molecular Evolutionary Analyses of the Spike Protein Gene and Spike Protein in the SARS-CoV-2 Omicron Subvariants.

To better understand the evolution of the SARS-CoV-2 Omicron subvariants, we performed molecular evolutionary analyses of the spike (S) protein gene/S protein using advanced bioinformatics technologies. First, time-scaled phylogenetic analysis estimated that a common ancestor of the Wuhan, Alpha, Beta, Delta variants, and Omicron variants/subvariants diverged in May 2020. After that, a common ancestor of the Omicron variant generated various Omicron subvariants over one year. Furthermore, a chimeric virus between the BM.1.1.1 and BJ.1 subvariants, known as XBB, diverged in July 2021, leading to the emergence of the prevalent subvariants XBB.1.5 and XBB.1.16. Next, similarity plot (SimPlot) data estimated that the recombination point (breakpoint) corresponded to nucleotide position 1373. As a result, XBB.1.5 subvariants had the 5' nucleotide side from the breakpoint as a strain with a BJ.1 sequence and the 3' nucleotide side as a strain with a BM.1.1.1 sequence. Genome network data showed that Omicron subvariants were genetically linked with the common ancestors of the Wuhan and Delta variants, resulting in many amino acid mutations. Selective pressure analysis estimated that the prevalent subvariants, XBB.1.5 and XBB.1.16, had specific amino acid mutations, such as V445P, G446S, N460K, and F486P, located in the RBD when compared with the BA.4 and BA.5 subvariants. Moreover, some representative immunogenicity-associated amino acid mutations, including L452R, F486V, R493Q, and V490S, were also found in these subvariants. These substitutions were involved in the conformational epitopes, implying that these mutations affect immunogenicity and vaccine evasion. Furthermore, these mutations were identified as positive selection sites. These results suggest that the S gene/S protein Omicron subvariants rapidly evolved, and mutations observed in the conformational epitopes may reduce the effectiveness of the current vaccine, including bivalent vaccines such as mRNA vaccines containing the BA.4/BA.5 subvariants.

Comprehensive transcripts analysis based on single-molecule real-time sequencing and Illumina sequencing provides insights into the mining of Toll-like receptor family in Schizothorax lissolabiatus.

Schizothorax lissolabiatus is an economically important cold-water fish species in southwestern China. Because of water pollution and habitat destruction, the number of wild populations has dramatically decreased. In this study, we used PacBio single-molecule real-time (SMRT) sequencing and Illumina sequencing to generate the first full-length transcriptome and transcriptome, respectively. A total of 19310 polished consensus reads (PC) were obtained, with an average length of 1379 bp and an N50 length of 1485 bp. Meanwhile, 12253 transcripts were successfully annotated as known homologous genes. The pathway annotation indicated that the enrichment and expression of most genes were mainly related to membrane, signal transduction and binding, and immune response. Furthermore, we identified 16 Toll-like receptors (TLRs) by mining the data from the transcripts. Phylogeny analysis showed that S. lissolabiatus TLR genes (slTLRs) supported the classification of TLRs into six families as in other vertebrates. Selection pressure analyses showed that 16 slTLRs revealed purification selection at the overall evolutionary selection. Further, positive selection signals were still detected in eight slTLRs, and most of the positive selection sites were located in the leucine-rich repeat region (LRR domain) associated with the recognition of pathogenic microorganisms, indicating that the function of these slTLR genes may be affected. Tissue specific expression analysis showed all slTLRs are present in kidney, spleen and liver but the relative expression varied among tissues. In conclusion, this study not only provided a valuable resource of transcripts for further research on S. lissolabiatus, but also contributed to improve the current understanding of the evolutionary history of immune-related genes and the TLR gene family in S. lissolabiatus.

Evolutionary dynamics and adaptive analysis of Seneca Valley virus.

Over the past 20years, the Seneca Valley virus (SVV) has emerged in various countries and regions around the world. Infected pigs display symptoms similar to foot-and-mouth disease and other vesicular diseases, causing severe economic losses to affected countries. In recent years, the number of SVV infections has been increasing in Brazil, China, and the United States. In this study, we comprehensively analyzed SVV genomic sequence data from the perspectives of evolutionary dynamics, phylogeography, and codon usage bias. We aimed to gain further insights into SVV's genetic diversity, spatiotemporal distribution patterns, and evolutionary adaptations. Phylogenetic analysis revealed that SVV has evolved into eight distinct lineages. Based on the results of phylogeographic analysis, it is speculated that the United States might have been the source of SVV, from where it subsequently spread to different countries and regions. Moreover, our analysis of positive selection sites in SVV capsid proteins suggests their potential importance in the process of receptor recognition. Finally, codon preference analysis indicates that natural selection has been a primary evolutionary driver influencing SVV codon usage bias. In conclusion, our in-depth investigation into SVV's origin, dissemination, evolution, and adaptation emphasizes the significance of SVV surveillance and control measures.

Genomic diversity and evolution analysis of severe fever with thrombocytopenia syndrome in East Asia from 2010 to 2022.

Conducting an up-to-date analysis on the genomic diversity and evolution patterns of severe fever with thrombocytopenia syndrome virus (SFTSV) is crucial for elucidating the underlying mechanisms of its emergency and pathogenicity, as well as assessing the extent of its threat to public health. Complete genome sequences of SFTSV were obtained from GenBank until December 19, 2022. A thorough phylogenetic analysis was conducted using comprehensive bioinformatics methods to estimate the genomic diversity and evolution. The phylogenetic classification of SFTSV strains yielded seven lineages (A-G) for each genome segment. SFTSV displayed notable variations in evolutionary patterns among different regions and segments, without a linear accumulation of nucleotide substitutions within segments and regions. The comprehensive analysis revealed 54 recombination events and 17 reassortment strains, including the first discovery of recombination events involving sea-crossing and species-crossing. Selection analysis identified three positive sites (2, 671, 1353) in RNA-dependent RNA polymerase, three positive sites (22, 298, 404) in glycoprotein, and two positive sites (9, 289) in nonstructural protein. No positive selection sites were found in nucleoprotein. Our study unveiled the existence of multiple evolutionary forces influencing SFTSV, contributing to its increasing genetic diversity, which had the potential to modify its antigenicity and pathogenicity. Furthermore, our study highlights the importance of tracking the spread of SFTSV across regions and species.

Evolutionary analysis of actinins(ACTN) gene family for positive selective sites by finding the dn/ds (ω) value

Identifying positively selected amino acid sites is an important approach for making a conclusion about the function of proteins; an amino acid site that is undergoing positive selection is likely to play a key role in the function of the protein. Here, we used codon-based substitution models and maximum-likelihood (ML) methods to identify positively selected sites that are likely to influence the cohesiveness and mechanics of the cytoskeleton by cross-linking actinins filaments and other cytoskeleton components to create a scaffold that imparts stability and forms a bridge between the cytoskeleton and signaling pathways. We use MEGA 7 and PAML for polygenetic Analysis and positive selective Sites. Alpha-actinins from skeletal muscle as a protein factor promoting the super precipitation of actomyosin and inducing the gelation of F-actin solutions. The protein is located in Z-bands of sarcomeres in skeletal and cardiac muscles as well as near the fascia adherents of intercalated discs in cardiac muscle. In smooth muscle, it exists in cytoplasmic dense bodies and membrane-associated dense plaques. We used PAMLX to find out specific sites in these DNA sequences. These results show the location and identification of positively selected sites in the ACTININS gene family. ACTN1, ACTN2, ACTN3 and ACTN4 had (248G,250L), (234I,236N ,238P, 239K, 252F, 260E, 261Q, 314Q, 590N, 607L, 869S), (NO) and (248G, 250L) positive selective sites respectively. The aims of our work are the prediction of these sites to find the specific location of these sites in 3D structure and to replace it with suitable amino acids for the cure of obesity, muscles, cancer, and other chronic diseases

Molecular evolution of Toll-like receptors in rodents.

Toll-like receptors (TLRs), the key sensor molecules in vertebrates, trigger the innate immunity and prime the adaptive immune system. The TLR family of rodents, the largest order of mammals, typically contains 13 TLR genes. However, a clear picture of the evolution of the rodent TLR family has not yet emerged and the TLR evolutionary patterns are unclear in rodent clades. Here, we analyzed the natural variation and the evolutionary processes acting on the TLR family in rodents at both the interspecific and population levels. Our results showed that rodent TLRs were dominated by purifying selection, but a series of positively selected sites (PSSs) primarily located in the ligand-binding domain was also identified. The numbers of PSSs differed among TLRs, and nonviral-sensing TLRs had more PSSs than those in viral-sensing TLRs. Gene-conversion events were found between TLR1 and TLR6 in most rodent species. Population genetic analyses showed that TLR2, TLR8, and TLR12 were under positive selection in Rattus norvegicus and R. tanezumi, whereas positive selection also acted on TLR5 and TLR9 in the former species, as well as TLR1 and TLR7 in the latter species. Moreover, we found that the proportion of polymorphisms with potentially functional change was much lower in viral-sensing TLRs than in nonviral-sensing TLRs in both of these rat species. Our findings revealed the first thorough insight into the evolution of the rodent TLR genetic variability and provided important novel insights into the evolutionary history of TLRs over long and short timescales.

Identification and Characterization of the Very-Low-Density Lipoprotein Receptor Gene from Branchiostoma belcheri: Insights into the Origin and Evolution of the Low-Density Lipoprotein Receptor Gene Family.

Low-density lipoprotein receptors (LDLRs) are a class of cell-surface endocytosis receptors that are mainly involved in cholesterol homeostasis and cellular signal transduction. Very-low-density lipoprotein receptors (VLDLRs), which are members of the LDLR family, have been regarded as multi-function receptors that fulfill diverse physiological functions. However, no VLDLR gene has been identified in protochordates to date. As a representative protochordate species, amphioxi are the best available example of vertebrate ancestors. Identifying and characterizing the VLDLR gene in amphioxi has high importance for exploring the evolutionary process of the LDLR family. With this study, a new amphioxus VLDLR gene (designated AmphiVLDLR) was cloned and characterized using RACE-PCR. The 3217 nt transcript of the AmphiVLDLR had a 2577 nt ORF, and the deduced 858 amino acids were highly conserved within vertebrate VLDLRs according to their primary structure and three-dimensional structure, both of which contained five characteristic domains. In contrast to other vertebrate VLDLRs, which had a conserved genomic structure organization with 19 exons and 18 introns, the AmphiVLDLR had 13 exons and 12 introns. The results of a selective pressure analysis showed that the AmphiVLDLR had numerous positive selection sites. Furthermore, the tissue expression of AmphiVLDLR using RT-qPCR showed that AmphiVLDLR RNA expression levels were highest in the gills and muscles, moderate in the hepatic cecum and gonads, and lowest in the intestines. The results of the evolutionary analysis demonstrated that the AmphiVLDLR gene is a new member of the VLDLR family whose family members have experienced duplications and deletions over the evolutionary process. These results imply that the functions of LDLR family members have also undergone differentiation. In summary, we found a new VLDLR gene homolog (AmphiVLDLR) in amphioxi. Our results provide insight into the function and evolution of the LDLR gene family.

Molecular Evolutionary Analyses of the RNA-Dependent RNA Polymerase (RdRp) Region and VP1 Gene in Human Norovirus Genotypes GII.P6-GII.6 and GII.P7-GII.6.

To understand the evolution of GII.P6-GII.6 and GII.P7-GII.6 strains, the prevalent human norovirus genotypes, we analysed both the RdRp region and VP1 gene in globally collected strains using authentic bioinformatics technologies. A common ancestor of the P6- and P7-type RdRp region emerged approximately 50 years ago and a common ancestor of the P6- and P7-type VP1 gene emerged approximately 110 years ago. Subsequently, the RdRp region and VP1 gene evolved. Moreover, the evolutionary rates were significantly faster for the P6-type RdRp region and VP1 gene than for the P7-type RdRp region and VP1 genes. Large genetic divergence was observed in the P7-type RdRp region and VP1 gene compared with the P6-type RdRp region and VP1 gene. The phylodynamics of the RdRp region and VP1 gene fluctuated after the year 2000. Positive selection sites in VP1 proteins were located in the antigenicity-related protruding 2 domain, and these sites overlapped with conformational epitopes. These results suggest that the GII.6 VP1 gene and VP1 proteins evolved uniquely due to recombination between the P6- and P7-type RdRp regions in the HuNoV GII.P6-GII.6 and GII.P7-GII.6 virus strains.

Evolutionary Study of DGAT1 Gene Indicates the Presence of Positive Selection in Mammals

This study was focused on the evolving role of DGAT1 using phylogenetic analysis. DGAT1 was retrieved from available databases of 37 mammalian genomes to detect signatures of positive selection in individual codons. Evolutionary analyses likelihood (REL), fast unconstrained Bayesian approximation (FUBAR), and Mixed-effects model evolution (MEME). REL detected 3 positive selection sites (at 109, 195, and 211 positions), and MEME detected 20 sites. Among them, 211 sites were confirmed as a strong positive selection of the DGAT1 gene. The codon model evolution (CME) model clustered amino acid substitution rates at three class levels (0.06, 0.20, and 0.40). Multiple hit tests showed that 20 of the DGAT1 triplet codes for amino acid formation were changed into new codes. In which, the GAC code for aspartate and GTC for valine had the maximum counts of 14 and 16. Strong evidence of positive selection of the DGAT1 gene in mammals was found.

Molecular evolutionary analyses of the fusion protein gene in human respirovirus 1

Few evolutionary studies of the human respiratory virus (HRV) have been conducted, but most of them have focused on HRV3. In this study, the full-length fusion (F) genes in HRV1 strains collected from various countries were subjected to time-scaled phylogenetic, genome population size, and selective pressure analyses. Antigenicity analysis was performed on the F protein. The time-scaled phylogenetic tree using the Bayesian Markov Chain Monte Carlo method estimated that the common ancestor of the HRV1 F gene diverged in 1957 and eventually formed three lineages. Phylodynamic analyses showed that the genome population size of the F gene has doubled over approximately 80 years. Phylogenetic distances between the strains were short (< 0.02). No positive selection sites were detected for the F protein, whereas many negative selection sites were identified. Almost all conformational epitopes of the F protein, except one in each monomer, did not correspond to the neutralising antibody (NT-Ab) binding sites. These results suggest that the HRV1 F gene has constantly evolved over many years, infecting humans, while the gene may be relatively conserved. Mismatches between computationally predicted epitopes and NT-Ab binding sites may be partially responsible for HRV1 reinfection and other viruses such as HRV3 and respiratory syncytial virus.

Genome-wide characterization of plant CTP:phosphocholine cytidylyltransferases through evolutionary, biochemical and structural analyses.

Phosphatidylcholine has essential functions in many eukaryotic cells, and its de novo biosynthesis is rate-limited by cytidine triphosphate:phosphocholine cytidylyltransferase (CCT). Although the biological and biochemical functions of CCT have been reported in mammals and several plants, this key enzyme has yet to be examined at a genome-wide level. As such, certain fundamental questions remain unanswered, including the evolutionary history, genetic and functional relationships, and structural variations among CCTs in the green lineage. In the current study, in-depth phylogenetic analysis, as well as the conservation and diversification in CCT gene structure and motif patterns, indicated that CCTs exist broadly in chlorophytes, bryophytes, lycophytes, monilophytes, gymnosperms, early-diverging angiosperms, monocots, and eudicots, and form eight relatively conserved clades. To further explore the potential function of selection pressure, we conducted extensive selection pressure analysis with a representative CCT gene, CCT1 from the model plant Arabidopsis thaliana (AthCCT1), and identified two positive selection sites, L59 and Q156. Site-directed mutagenesis and in vitro enzyme assays demonstrated that these positively selected sites were indeed important for the activity and substrate affinity of AthCCT1, and subsequent 3D structure analyses explained the potential biochemical mechanism. Taken together, our results unraveled the evolution and diversity of CCTs in the green lineage, as well as their association with the enzyme's biochemical and structural properties, and expanded our understanding of this important enzyme at the genome-wide level.

Structural Characterization of Trivalvaria costata Chloroplast Genome and Molecular Evolution of rps12 Gene in Magnoliids

Chloroplasts (cp) are important organelles in plant cells that have been widely used in phylogenetic, molecular evolution, and gene expression studies due to their conserved molecular structure. In this study, we obtained the complete cp genome of Trivalvaria costata (Annonaceae) and analyzed its structural characteristics. Additionally, we analyzed the rps12 gene in the phylogenetic framework of magnoliids. The T. costata cp genome comprises 162,002 bp and contains 132 genes. We detected 48 simple sequence repeats (SSRs) and identified 29 high-frequency codons as well as 8 optimal codons. Our multiple analyses show that codon usage bias is mainly influenced by natural selection. For the first time, we found the rps12 gene to be entirely located in the IR region (in Annona). In groups with exon 1 located in the single-copy (SC) region and exons 2–3 located in the inverted repeat (IR) region, the transition rate and synonymous substitution rate of exon 1 were higher than those of exons 2–3. Adaptive evolution identified a positive selection site (116) located in the 310-helix region, suggesting that the rps12 gene may undergo adaptive changes during the evolutionary history of magnoliids. This study enhances our knowledge regarding genetic information on T. costata and provides support for reduced substitution rates in the IR region.

Characterization of Two Mitogenomes of Hyla sanchiangensis (Anura: Hylidae), with Phylogenetic Relationships and Selection Pressure Analyses of Hylidae.

Hyla sanchiangensis (Anura: Hylidae) is endemic to China and is distributed across Anhui, Zhejiang, Fujian, Guangdong, Guangxi, Hunan, and Guizhou provinces. The mitogenomes of H. sanchiangensis from two different sites (Jinxiu, Guangxi, and Wencheng, Zhejiang) were sequenced. Phylogenetic analyses were conducted, including 38 mitogenomes of Hylidae from the NCBI database, and assessed the phylogenetic relationship of H. sanchiangensis within the analyzed dataset. Two mitogenomes of H. sanchiangensis showed the typical mitochondrial gene arrangement with 13 protein-coding genes (PCGs), two ribosomal RNA genes (12S rRNA and 16S rRNA), 22 transfer RNA (tRNA) genes, and one non-coding control region (D-loop). The lengths of the 12S rRNA and 16S rRNA genes from both samples (Jinxiu and Wencheng) were 933 bp and 1604 bp, respectively. The genetic distance (p-distance transformed into percent) on the basis of the mitogenomes (excluding the control region) of the two samples was calculated as 4.4%. Hyla sanchiangensis showed a close phylogenetic relationship with the clade of (H. annectans + H. tsinlingensis), which was supported by ML and BI analyses. In the branch-site model, five positive selection sites were found in the clade of Hyla and Dryophytes: Cytb protein (at position 316), ND3 protein (at position 85), and ND5 protein (at position 400) have one site, respectively, and two sites in ND4 protein (at positions 47 and 200). Based on the results, we hypothesized that the positive selection of Hyla and Dryophytes was due to their experience of cold stress in historical events, but more evidence is needed to support this conclusion.

DNA repair-related genes are the key for stony coral ancestors to survive under elevated levels of UVR

DNA damage response (DDR) is a complicated network to defend against physical or chemical changes in DNA in all animals. Elevated levels of ultraviolet radiation (UVR) caused DNA damage, which was a reason for the mass extinction that occurred at the Devonian/Carboniferous (D/C) boundary approximately 359 million years ago (Ma). However, the molecular adaptation of the stony coral ancestors that strangely survived the D/C boundary mass extinction is not well understood. In the present study, the molecular clock analysis using fourfold degenerate sites of 1,463 homologous genes of different stony coral species (a representative group of marine organisms with calcareous skeletons) suggested that their common ancestors originated 384.24 Ma, i.e., slightly earlier than the D/C transition. We identified 21 rapidly evolving genes (REGs) and 49 positive selection genes (PSGs) that were significantly enriched in diverse pathways, including the mitotic cell cycle process, intracellular protein transport, and DNA synthesis involved in DNA repair. Interestingly, four REGs and 21 PSGs were significantly enriched in DDR pathways, including the mitotic cell cycle process, DNA synthesis involved in DNA repair, and cellular response to DNA damage stimulus pathways. We hypothesize that enriched DDR genes are likely involved in the enhanced ability of ancient stony corals to detect and repair DNA damage. For example, the DNA polymerase epsilon catalytic subunit (POLE) gene, which encodes the DNA polymerase epsilon catalytic subunit A, mediates interaction with the other three catalytic subunits through its nonenzymatic carboxy-terminal domain. POLE may potentially enhance the binding ability to other subunits to strengthen its function because eight positive selection sites were distributed in the C-terminal of POLE and on the surface of the simulated 3D protein model. Therefore, our results demonstrated for the first time that the precise transfer of DNA information may help stony coral ancestors survive in elevated levels of ultraviolet radiation, suggesting that DDR levels may be critical to the environmental adaptation of calcareous skeletal organisms during climate change.

Molecular characterization of human respiratory syncytial virus in Seoul, South Korea, during 10 consecutive years, 2010-2019.

Respiratory syncytial virus (RSV) is the leading cause of lower respiratory tract infections and hospitalization in infants and young children. Here, we analyzed the genetic diversity of RSV using partial G gene sequences in 84 RSV-A and 78 RSV- B positive samples collected in Seoul, South Korea, for 10 consecutive years, from 2010 to 2019. Our phylogenetic analysis revealed that RSV-A strains were classified into either the ON1 (80.9%) or NA1 (19.0%) genotypes. On the other hand, RSV-B strains demonstrated diversified clusters within the BA genotype. Notably, some sequences designated as BA-SE, BA-SE1, and BA-DIS did not cluster with previously identified BA genotypes in the phylogenetic trees. Despite this, they did not meet the criteria for the assignment of a new genotype based on recent classification methods. Selection pressure analysis identified three positive selection sites (amino acid positions 273, 274, and 298) in RSV-A, and one possible positive selection site (amino acid position 296) in RSV-B, respectively. The mean evolutionary rates of Korean RSV-A from 1999 to 2019 and RSV-B strains from 1991 and 2019 were estimated at 3.51 × 10-3 nucleotides (nt) substitutions/site/year and 3.32 × 10-3 nt substitutions/site/year, respectively. The population dynamics in the Bayesian skyline plot revealed fluctuations corresponding to the emergence of dominant strains, including a switch of the dominant genotype from NA1 to ON1. Our study on time-scaled cumulative evolutionary analysis contributes to a better understanding of RSV epidemiology at the local level in South Korea.