Positive Selection Pressure Research Articles

Pan-Genome-Wide Investigation and Co-Expression Network Analysis of HSP20 Gene Family in Maize.

Heat shock protein 20 (HSP20) is a diverse and functionally important protein family that plays a crucial role in plants' tolerance to various abiotic stresses. In this study, we systematically analyzed the structural and functional characteristics of the HSP20 gene family within the Zea pan-genome. By identifying 56 HSP20 pan-genes, we revealed the variation in the number of these genes across different maize inbreds or relatives. Among those 56 genes, only 31 are present in more than 52 inbreds or relatives. Further phylogenetic analysis classified these genes into four major groups (Class A, B, C, D) and explored their diversity in subcellular localization, physicochemical properties, and the terminal structures of those HSP20s. Through collinearity analysis and Ka/Ks ratio calculations, we found that most HSP20 genes underwent purifying selection during maize domestication, although a few genes showed signs of positive selection pressure. Additionally, expression analysis showed that several HSP20 genes were significantly upregulated under high temperatures, particularly in tassels and leaves. Co-expression network analysis revealed that HSP20 genes were significantly enriched in GO terms related to environmental stress responses, suggesting that HSP20 genes not only play key roles in heat stress but may also be involved in regulating various other biological processes, such as secondary metabolism and developmental processes. These findings expand our understanding of the functions of the maize HSP20 family and provide new insights for further research into maize's response mechanisms to environmental stresses.

Characterization, expressional and evolutionary analysis of five fish-specific CCRs (CCR4La, CCR4Lc, CCR12a1, CCR12a2, and CCR12b) in largemouth bass (Micropterus salmoides)

CC chemokine receptors (CCRs), the numbers of the G protein-coupled receptor (GPCR) superfamily, had crucial roles in treating infection, inflammation, and tissue damage by binding to their ligands. In this study, five fish-specific CCRs, namely CCR4La, CCR4Lc, CCR12a1, CCR12a2, and CCR12b, were identified in largemouth bass (Micropterus salmoides). The correction of nomenclatures of these CCRs were confirmed by phylogenetic analysis, structural analysis and genomic synteny analysis. Following 1×106 CFU/mL and 1×107 CFU/mL Edwardsiella piscicida infection, these five CCRs were significantly induced in spleen of largemouth bass, indicating their important roles in the immune response against bacterial infection. Selection pressure analysis revealed that CCR4La, CCR4Lc, CCR12a1, and CCR12a2 underwent negative selection pressure, whereas CCR12b experienced positive selection pressure. Robust selection site detection methods identified that positive selected sites of CCR4La, CCR4Lc, CCR12a1, and CCR12a2 mainly distributed in their extracellular regions, which involved in ligand binding and pathogen interaction. Similarly, the positive selected sites of CCR12b were also located in its extracellular regions. The accuracy of the pressure selected sites were also validated by molecular docking analysis. The potential ligands for these five CCRs were identified by molecular docking analysis, with finding that CCL3 and CCL5 might be the ligands of largemouth bass CCR4La/Lc, and CCL5, CCL8, CCL7, CCL13 and CCL26 might be that of largemouth bass CCR12a1/a2/b. Our results provided basis for elucidating the functions of chemokine-receptor complex in largemouth bass.

The mitochondrial genomes of the reef-dwelling spiny lobsters Panulirus echinatus and Panulirus interruptus with insights into the phylogeny and adaptive evolution of protein-coding genes in the Achelata

AbstractTemperature and oxygen levels drive the evolution of morphological, behavioral, and physiological traits in marine invertebrates, including crustaceans. Environmental conditions are also expected to prompt the adaptive evolution of mitochondrial protein-coding genes (PCGs), which are vital for energy production via the oxidative phosphorylation pathway. We formally tested for adaptive evolution in mitochondrial protein-coding genes in representatives of the decapod infraorder Achelata, including two spiny lobsters, Panulirus echinatus and P. interruptus, for which we sequenced complete mitochondrial genomes (15,644 and 15,659 bp long, respectively). A phylomitogenomic analysis supported the monophyly of the genus Panulirus, the families Palinuridae and Scyllaridae, and the infraorder Achelata. Over the strong negative selection background observed for mitochondrial PCGs in the Achelata, signatures of positive selective pressure were detected within PCGs in equatorial Panulirus spp. and deepwater Scyllaridae. In Panulirus spp. inhabiting equatorial latitudes with consistently high temperatures, the Datamonkey analysis RELAX suggested intensified purifying selection strength in 9 of the 13 PCGs and relaxation in purifying selection strength in atp6, while aBSREL, BUSTED, and MEME recovered signatures of positive selection on PCGs within Complex I, III, and IV PCGs. Likewise, in Scyllaridae species inhabiting depths with low-oxygen levels, RELAX indicated relaxed selection strength in 6 of the 13 PCGs, while aBSREL, BUSTED, and MEME recovered signatures of positive selection on PCGs within Complexes I, III, IV, and V. The newly assembled mitochondrial genomes of P. echinatus and P. interruptus represent new genomic resources to aid with the conservation and management of lobsters targeted by major fisheries and contribute to our understanding of how environmental conditions drive adaptive evolution in spiny and slipper lobster mitochondrial PCGs.

Evolutionary origins and innovations sculpting the mammalian PRPS enzyme complex.

The phosphoribosyl pyrophosphate synthetase (PRPS) enzyme conducts a chokepoint reaction connecting central carbon metabolism and nucleotide production pathways, making it essential for life 1,2 . Here, we show that the presence of multiple PRPS-encoding genes is a hallmark trait of eukaryotes, and we trace the evolutionary origins and define the individual functions of each of the five mammalian PRPS homologs - three isozymes (one testis-restricted) 3,4 and two non-enzymatic associated proteins (APs) 5,6 - which we demonstrate operate together as a large molecular weight complex capable of attaining a heterogeneous array of functional multimeric configurations. Employing a repertoire of isogenic fibroblast clones in all viable individual or combinatorial assembly states, we define preferential interactions between subunits, and we show that cells lacking PRPS2, PRPSAP1, and PRPSAP2 render PRPS1 into aberrant homo-oligomeric assemblies with diminished metabolic flux and impaired proliferative capacity. We demonstrate how numerous evolutionary innovations in the duplicated genes have created specialized roles for individual complex members and identify translational control mechanisms that enable fine-tuned regulation of PRPS assembly and activity, which provide clues into the positive and negative selective pressures that facilitate metabolic flexibility and tissue specialization in advanced lifeforms. Collectively, our study demonstrates how evolution has transformed a single PRPS gene into a multimeric complex endowed with functional and regulatory features that govern cellular biochemistry.

Whole genome characteristics of hedgehog coronaviruses from Poland and analysis of the evolution of the Spike protein for its interspecies transmission potential

BackgroundThe hedgehogs have been recently identified as possible reservoir of Middle East respiratory syndrome coronavirus like (MERS-CoV-like). These viruses were classified as a distinct Betacoronavirus erinacei (BCoV-Eri) species within the MerBCoV-Eriirus subgenus. As coronaviruses are known for their ability to jump between different hosts, including humans, this can pose a particular threat to people in direct contact with hedgehogs, such as those working at animal asylums. Our previous studies have shown the presence of BCoV-Eri strains in animals collected in the wildlife rehabilitation centre. This study aimed to investigate the presence of CoV in subsequent hedgehogs collected from the urban area of Poland and their molecular characteristics.ResultsMonitoring for the presence of coronavirus infection in hedgehogs revealed five positive individuals. The presence of BCoV-Eri was found in a total of 20% of animals tested. Our analyses revealed no correlation between CoVs positivity and animal health conditions but a higher probability of such infection in juveniles and females. The whole genome of two Polish Hedgehog coronavirus 1 strains were sequenced and compared with available counterparts from European and Asian countries. Phylogenetic analysis showed that both CoV strains formed common cluster with other similar MerBCoV-Eriirus, but they were also found to be genetically variable and most changes in the S protein were identified. Our analysis revealed that some S protein sites of the Hedgehog coronavirus 1 strains evolved under positive selection pressure and of five such sites, three are in the S1 region while the other two in the S2 region of the Spike.ConclusionsBCoV-Eri is to some extent prevalent in wildlife asylums in Poland. Given that the S protein of BCoVs-Eri is highly variable and that some sites of this protein evolve under positive selection pressure, these strains could potentially acquire a favourable feature for cross-species transmission. Consequently, the threat to humans working in such asylums is particularly high. Adequate biosecurity safeguards, but also human awareness of such risks, are therefore essential.

Cooperative folding as a molecular switch in an evolved antibody binder

Designing proteins with tunable activities from easily accessible external cues remains a biotechnological challenge. Here, we set out to create a small antibody-binding domain equipped with a molecular switch inspired by the allosteric response to calcium seen in naturally derived proteins like calmodulin. We have focused on one of the three domains of Protein G that show inherent affinity to antibodies. By combining a semi-rational protein design with directed evolution, we engineered novel variants containing a calcium-binding loop rendering the inherent antibody affinity calcium-dependent. The evolved variants resulted from a designed selection strategy subjecting them to negative and positive selection pressures focused on conditional antibody-binding. Hence, these variants contained molecular “on/off” switches, controlling the target affinity towards antibody fragments simply by the presence or absence of calcium. From NMR spectroscopy we found that the molecular mechanism underlying the evolved switching behavior was a coupled calcium-binding and folding event where the target binding surface was intact and functional only in the presence of bound calcium. Notably, it was observed that the response to the employed selection pressures gave rise to the evolution of a cooperative folding mechanism. This observation illustrates why the cooperative folding reaction is an effective solution seen repeatedly in the natural evolution of fine-tuned macromolecular recognition. Engineering binding moieties to confer conditional target interaction has great potential due to the exquisite interaction control that is tunable to application requirements. Improved understanding of the molecular mechanisms behind regulated interactions is crucial to unlock how to engineer switchable proteins useful in a variety of biotechnological applications.

Genetic diversity and evolutionary significance of Cahi-DRB1 gene in Indian goat from semi-arid regions

Unlocking the genetic diversity of the Cahi-DRB1 gene in Indian goat breeds sheds light on their immune system's resilience. While extensive studies have elucidated MHC diversity in cattle and sheep, research on goats remains limited. Thus, this study aimed to investigate the genetic diversity of cahi-DRB1 across five goat breeds in India, encompassing a total of 233 individual goats. Using a PCR sequence-based typing approach, DRB1 fragments (285 bp) were amplified and sequenced. The PCR products were purified and sequenced. Amino acid variability within breeds was evaluated using the Wu-Kabat index at peptide-binding pockets. Selection analysis included assessing positive selection through relative rates of nonsynonymous and synonymous mutations, Z-test, and Tajima's test of neutrality. The findings revealed notable heterozygosity across breeds (mean: 100 %), with nucleotide content averaging 22.51 % for A, 18.62 % for T, 24.33 % for C, and 34.35 % for G. The GC content was rich (59.26 %) across breeds for the DRB1 region. The mean distance for DRB1 was highest in Barbari (0.20) and lowest in Osmanabadi goat (0.12), with a higher rate of transversions observed. The Wu-Kabat index assessments underscored substantial genetic variation within the peptide binding sites (PBS) of DRB1, with specific residues (12, 14, 33, 58, and 87) showing polymorphism relevant to antigen presentation. Positive selection pressure was evident across all breeds for the DRB1 locus (dN>dS), highlighting its evolutionary significance. The rich diversity of the DRB1 gene underscores the fitness of the population and the significance of this locus for future selection programs. These results highlight the highly polymorphic nature of the DRB1 locus in Indian goat breeds, emphasizing its crucial role in modulating immune responses. Consequently, our study offers valuable insights for developing targeted breeding strategies to enhance disease resistance in goats, thereby providing practical implications for sustainable livestock management.

Enriched G4 forming repeats in the human genome are associated with robust well-coordinated transcription and reduced cancer transcriptome variation

Non-B DNA G-quadruplex (G4) structures with guanine (G) runs of 2-4 repeats can trigger opposing experimental transcriptional impacts. Here, we employed bioinformatic algorithms to comprehensively assess correlations of steady-state RNA transcript levels with all putative G4 sequence (pG4) locations genome-wide in three mammalian genomes and in normal and tumor human tissues. The human pG4-containing gene set displays higher expression levels than the set without pG4, supporting and extending some prior observations. pG4 enrichment at transcription start sites (TSS) in human, but not chimpanzee and mouse genomes, suggests possible positive selection pressure for pG4 at human TSS, potentially driving genome rewiring and gene expression divergence between human and chimpanzee. Comprehensive bioinformatic analyses revealed lower pG4-containing gene set variability in humans and among different pG4 genes in tumors. As G4 stabilizers are under therapeutic consideration for cancer and pathogens, such distinctions between human normal and tumor G4s along with other species merit attention. Furthermore, in germline and cancer sequences, the most mutagenic pG4 mapped to regions promoting alternative DNA structures. Overall findings establish high pG4 at TSS as a human genome attribute statistically associated with robust well-coordinated transcription and reduced cancer transcriptome variation with implications for biology, model organisms, and medicine.

The Chrysosplenium sinicum genome provides insights into adaptive evolution of shade plants

Chrysosplenium sinicum, a traditional Tibetan medicinal plant, can successfully thrive in low-light environments for long periods of time. To investigate the adaptive evolution of shade plants in low-light environments, we generated a chromosome-scale genome assembly (~320 Mb) for C. sinicum by combining PacBio sequencing and Hi-C technologies. Based on our results, gene families related to photosynthesis and cell respiration greatly expanded and evolved in C. sinicum genome due to intracellular DNA transfer from organelle genome to nuclear genome. Under positive selective pressure, adaptive evolution of light-harvesting complex II (LHCII) component protein CsLhcb1s resulted in the expansion of threonine residues at the phosphorylation site of STN7 kinase, potentially establishing a crucial genomic foundation for enhancing C. sinicum’s adaptability in low-light environments. Through transcriptome and metabolome analysis, we identified chrysosplenol and chrysosplenoside as predominant flavonoid metabolites of C. sinicum and predicted their synthesis pathways. In addition, analysis of alternative splicing (AS) revealed that AS events help regulate state transition and flavonoid biosynthesis. The present study provides new insights into the genomes of shade plants exposed to low-light conditions and adaptive evolution of these genomes; in addition, the results improve our current knowledge on the biosynthetic and regulatory processes of chrysosplenol and chrysosplenoside.

Intercompatibility of eukaryotic and Asgard archaea ribosome-translocon machineries

In all domains of life, the ribosome-translocon complex inserts nascent transmembrane proteins into, and processes and transports signal peptide-containing proteins across, membranes. Eukaryotic translocons are anchored in the endoplasmic reticulum, while the prokaryotic complexes reside in cell membranes. Phylogenetic analyses indicate the inheritance of eukaryotic Sec61/oligosaccharyltransferase/translocon-associated protein translocon subunits from an Asgard archaea ancestor. However, the mechanism for translocon migration from a peripheral membrane to an internal cellular compartment (the proto-endoplasmic reticulum) during eukaryogenesis is unknown. Here we show compatibility between the eukaryotic ribosome-translocon complex and Asgard signal peptides and transmembrane proteins. We find that Asgard translocon proteins from Candidatus Prometheoarchaeum syntrophicum strain Candidatus Prometheoarchaeum syntrophicum strain MK-D1, a Lokiarchaeon confirmed to contain no internal cellular membranes, are targeted to the eukaryotic endoplasmic reticulum on ectopic expression. Furthermore, we show that the cytoplasmic domain of Candidatus Prometheoarchaeum syntrophicum strain MK-D1 oligosaccharyltransferase 1 (ribophorin I) can interact with eukaryotic ribosomes. Our data indicate that the location of existing ribosome-translocon complexes, at the protein level, determines the future placement of yet-to-be-translated translocon subunits. This principle predicts that during eukaryogenesis, under positive selection pressure, the relocation of a few translocon complexes to the proto-endoplasmic reticulum will have contributed to propagating the new translocon location, leading to their loss from the cell membrane.

Bebtelovimab-bound SARS-CoV-2 RBD mutants: resistance profiling and validation with escape mutations, clinical results, and viral genome sequences.

The dynamic evolution of SARS-CoV-2 variants necessitates ongoing advancements in therapeutic strategies. Despite the promise of monoclonal antibody (mAb) therapies like bebtelovimab, concerns persist regarding resistance mutations, particularly single-to-multipoint mutations in the receptor-binding domain (RBD). Our study addresses this by employing interface-guided computational protein design to predict potential bebtelovimab-resistance mutations. Through extensive physicochemical analysis, mutational preferences, precision-recall metrics, protein-protein docking, and energetic analyses, combined with all-atom, and coarse-grained molecular dynamics (MD) simulations, we elucidated the structural-dynamics-binding features of the bebtelovimab-RBD complexes. Identification of susceptible RBD residues under positive selection pressure, coupled with validation against bebtelovimab-escape mutations, clinically reported resistance mutations, and viral genomic sequences enhances the translational significance of our findings and contributes to a better understanding of the resistance mechanisms of SARS-CoV-2.

Genetic diversity and natural selection analysis of VAR2CSA and vir genes: implication for vaccine development

Variable surface antigens (VSAs) encoded by var and vir genes in Plasmodium falciparum and Plasmodium vivax, respectively, are known to be involved in malaria pathogenesis and host immune escape through antigenic variations. Knowledge of the genetic diversity of these antigens is essential for malaria control and effective vaccine development. In this study, we analysed the genetic diversity and evolutionary patterns of two fragments (DBL2X and DBL3X) of VAR2CSA gene and four vir genes (vir 4, vir 12, vir 21 and vir 27) from different endemic regions, including Southeast Asia and sub-Saharan Africa. High levels of segregating sites (S) and haplotype diversity (Hd) were observed in both var and vir genes. Among vir genes, vir 12 (S = 131, Hd = 0.996) and vir 21 (S = 171, Hd = 892) were found to be more diverse as compared to vir 4 (S = 11, Hd = 0.748) and vir 27 (S = 23, Hd = 0.814). DBL2X (S = 99, Hd = 0.996) and DBL3X (S = 307, Hd = 0.999) fragments showed higher genetic diversity. Our analysis indicates that var and vir genes are highly diverse and follow the similar evolutionary pattern globally. Some codons showed signatures of positive or negative selection pressure, but vir and var genes are likely to be under balancing selection. This study highlights the high variability of var and vir genes and underlines the need of functional experimental studies to determine the most relevant allelic forms for effective progress towards vaccine formulation and testing.

Exploring selection signatures in the divergence and evolution of lipid droplet (LD) associated genes in major oilseed crops

BackgroundOil bodies or lipid droplets (LDs) in the cytosol are the subcellular storage compartments of seeds and the sites of lipid metabolism providing energy to the germinating seeds. Major LD-associated proteins are lipoxygenases, phospholipaseD, oleosins, TAG-lipases, steroleosins, caleosins and SEIPINs; involved in facilitating germination and enhancing peroxidation resulting in off-flavours. However, how natural selection is balancing contradictory processes in lipid-rich seeds remains evasive. The present study was aimed at the prediction of selection signatures among orthologous clades in major oilseeds and the correlation of selection effect with gene expression.ResultsThe LD-associated genes from the major oil-bearing crops were analyzed to predict natural selection signatures in phylogenetically close-knit ortholog clusters to understand adaptive evolution. Positive selection was the major force driving the evolution and diversification of orthologs in a lineage-specific manner. Significant positive selection effects were found in 94 genes particularly in oleosin and TAG-lipases, purifying with excess of non-synonymous substitution in 44 genes while 35 genes were neutral to selection effects. No significant selection impact was noticed in Brassicaceae as against LOX genes of oil palm. A heavy load of deleterious mutations affecting selection signatures was detected in T-lineage oleosins and LOX genes of Arachis hypogaea. The T-lineage oleosin genes were involved in mainly anther, tapetum and anther wall morphogenesis. In Ricinus communis and Sesamum indicum > 85% of PLD genes were under selection whereas selection pressures were low in Brassica juncea and Helianthus annuus. Steroleosin, caleosin and SEIPINs with large roles in lipid droplet organization expressed mostly in seeds and were under considerable positive selection pressures. Expression divergence was evident among paralogs and homeologs with one gene attaining functional superiority compared to the other. The LOX gene Glyma.13g347500 associated with off-flavor was not expressed during germination, rather its paralog Glyma.13g347600 showed expression in Glycine max. PLD-α genes were expressed on all the tissues except the seed,δ genes in seed and meristem while β and γ genes expressed in the leaf.ConclusionsThe genes involved in seed germination and lipid metabolism were under strong positive selection, although species differences were discernable. The present study identifies suitable candidate genes enhancing seed oil content and germination wherein directional selection can become more fruitful.

Assembly and Comparative Analyses of the Chloroplast Genomes of the Threatened Plant Rosa anemoniflora

Due to insufficient molecular biology and genetic research on Rosa anemoniflora, this endangered plant has not yet received effective protection. Therefore, the complete chloroplast genome sequence of R. anemoniflora, along with comparative analysis of the chloroplast genomes of related species, is necessary and crucial for reconstructing phylogenetic relationships and developing genetic markers to conserve these species. A series of analyses, including genome structure, GC content, gene number, selection pressure, and nucleotide diversity, were performed by comparing the chloroplast genomes of R. anemoniflora and its relatives. The results indicate that the chloroplast genomes of R. anemoniflora and its close relatives are highly conserved in all genome characteristics, and all protein-coding genes in R. anemoniflora have not experienced significant positive or negative selection pressures. Comparative analysis revealed several variation hotspots, such as the atpH-atpI region, which can serve as a DNA barcode for distinguishing R. anemoniflora from its close relatives. Finally, the results confirmed that R. anemoniflora belongs to Rosa section Synstylae and that R. anemoniflora and its close relatives likely originated from the Fujian and Taiwan regions of China and diverged approximately 3.24 million years ago. This study provides crucial information for future biodiversity conservation and genetic resource management.

Antarctica's hidden mycoviral treasures in fungi isolated from mosses: A first genomic approach.

This study investigates the presence of mycoviruses in Antarctic fungi and elucidates their evolutionary relationships. To achieve this, we aligned mycoviral gene sequences with genomes of previously sequenced Antarctic endophytic fungi, made available by our research group and accessible via Joint Genome Institute. Our findings reveal that the most prevalent genetic regions in all endophytic fungi are homologous to Partitiviruses, Baculoviridae, and Phycodnaviridae. These regions display evidence of positive selection pressure, suggesting genetic diversity and the accumulation of nonsynonymous mutations. This phenomenon implies a crucial role for these regions in the adaptation and survival of these fungi in the challenging Antarctic ecosystems. The presence of mycoviruses in Antarctic endophytic fungi may indicate shared survival strategies between the virus and its host, shedding light on their evolutionary dynamics. This study underscores the significance of exploring mycoviruses within endophytic fungi and their contributions to genetic diversity. Future research avenues could delve into the functional implications of these conserved mycoviral genetic regions in Antarctic endophytic fungi, providing a comprehensive understanding of this intriguing association and genomic retention of viral region in fungi.

Molecular characterization of positively selected genes contributing aquatic adaptation in marine mammals.

Marine mammals, which have evolved independently into three distinct lineages, share common physiological features that contribute to their adaptation to the marine environment. To identify positively selected genes (PSGs) for adaptation to the marine environment using available genomic data from three taxonomic orders: cetaceans, pinnipeds, and sirenians. Based on the genomes within each group of Artiodactyla, Carnivora and Afrotheria, we performed selection analysis using the branch-site model in CODEML. Based on the branch-site model, 460, 614, and 359 PSGs were predicted for the cetaceans, pinnipeds, and sirenians, respectively. Functional enrichment analysis indicated that genes associated with hemostasis were positively selected across all lineages of marine mammals. We observed positive selection signals for the hemostasis and coagulation-related genes plasminogen activator, urokinase (PLAU), multimerin 1 (MMRN1), gamma-glutamyl carboxylase (GGCX), and platelet endothelial aggregation receptor 1 (PEAR1). Additionally, we found out that the sodium voltage-gated channel alpha subunit 9 (SCN9A), serine/arginine repetitive matrix 4 (SRRM4), and Ki-ras-induced actin-interacting protein (KRAP) are under positive selection pressure and are associated with cognition, neurite outgrowth, and IP3-mediated Ca2 + release, respectively. This study will contribute to our understanding of the adaptive evolution of marine mammals by providing information on a group of candidate genes that are predicted to influence adaptation to aquatic environments, as well as their functional characteristics.

Whole genome sequence analysis of shallot virus X from India reveals it to be a natural recombinant with positive selection pressure

BackgroundShallots are infected by various viruses like Onion yellow dwarf virus (OYDV), Leek yellow stripe virus (LYSV), Shallot latent virus (SLV) and Shallot virus X (ShVX). In India, they have been found to be persistently infected by ShVX. ShVX also infects onion and garlic in combination with other carlaviruses and potyviruses. ShVX is a member of genus Allexivirus of family Alphaflexiviridae. ShVX has a monopartite genome, which is represented by positive sense single-stranded RNA. Globally, only six complete and 3 nearly complete genome sequences of ShV X are reported to date. This number is insufficient to measure a taxon’s true molecular diversity. Moreover, the complete genome sequence of ShVX from Asia has not been reported as yet. Therefore, this study was undertaken to generate a complete genome sequence of ShVX from India.ResultsShallot virus X (ShVX) is one of the significant threats to Allium crop production. In this study, we report the first complete genome sequence of the ShVX from India through Next-generation sequencing (NGS). The complete genome of the ShVX (Accession No. OK104171), from this study comprised 8911 nucleotides. In-silico analysis of the sequence revealed variability between this isolate and isolates from other countries. The dissimilarities are spread all over the genome specifically some non-coding intergenic regions. Statistical analysis of individual genes for site-specific selection indicates a positive selection in NABP region. The presence of a recombination event was detected in coat protein region. The sequence similarity percentage and phylogenetic analysis indicate ShVX Indian isolate is a distinctly different isolate. Recombination and site-specific selection may have a function in the evolution of this isolate. This is the first detailed study of the ShVX complete genome sequence from Southeast Asia.ConclusionThis study presents the first report of the entire genome sequence of an Indian isolate of ShVX along with an in-depth exploration of its evolutionary traits. The findings highlight the Indian variant as a naturally occurring recombinant, emphasizing the substantial role of recombination in the evolution of this viral species. This insight into the molecular diversity of strains within a specific geographical region holds immense significance for comprehending and forecasting potential epidemics. Consequently, the insights garnered from this research hold practical value for shaping ShVX management strategies and providing a foundation for forthcoming studies delving into its evolutionary trajectory.

Evolutionary dynamics and comparative pathogenicity of clade 2.3.4.4b H5 subtype avian influenza viruses, China, 2021–2022

The recent concurrent emergence of H5N1, H5N6, and H5N8 avian influenza viruses (AIVs) has led to significant avian mortality globally. Since 2020, frequent human-animal interactions have been documented. To gain insight into the novel H5 subtype AIVs (i.e., H5N1, H5N6 and H5N8), we collected 6102 samples from various regions of China between January 2021 and September 2022, and identified 41 H5Nx strains. Comparative analyses on the evolution and biological properties of these isolates were conducted. Phylogenetic analysis revealed that the 41 H5Nx strains belonged to clade 2.3.4.4b, with 13 related to H5N1, 19 to H5N6, and 9 to H5N8. Analysis based on global 2.3.4.4b viruses showed that all the viruses described in this study were likely originated from H5N8, exhibiting a heterogeneous evolutionary history between H5N1 and H5N6 during 2015–2022 worldwide. H5N1 showed a higher rate of evolution in 2021–2022 and more sites under positive selection pressure in 2015–2022. The antigenic profiles of the novel H5N1 and H5N6 exhibited notable variations. Further hemagglutination inhibition assay suggested that some A(H5N1) viruses may be antigenically distinct from the circulating H5N6 and H5N8 strains. Mammalian challenge assays demonstrated that the H5N8 virus (21GD001_H5N8) displayed the highest pathogenicity in mice, followed by the H5N1 virus (B1557_H5N1) and then the H5N6 virus (220086_H5N6), suggesting a heterogeneous virulence profile of H5 AIVs in the mammalian hosts. Based on the above results, we speculate that A(H5N1) viruses have a higher risk of emergence in the future. Collectively, these findings unveil a new landscape of different evolutionary history and biological characteristics of novel H5 AIVs in clade 2.3.4.4b, contributing to a better understanding of designing more effective strategies for the prevention and control of novel H5 AIVs.

Complete Chloroplast Genome of Alternanthera sessilis and Comparative Analysis with Its Congeneric Invasive Weed Alternanthera philoxeroides.

Alternanthera sessilis is considered the closest relative to the invasive weed Alternanthera philoxeroides in China, making it an important native species for studying the invasive mechanisms and adaptations of A. philoxeroides. Chloroplasts play a crucial role in a plant's environmental adaptation, with their genomes being pivotal in the evolution and adaptation of both invasive and related species. However, the chloroplast genome of A. sessilis has remained unknown until now. In this study, we sequenced and assembled the complete chloroplast genome of A. sessilis using high-throughput sequencing. The A. sessilis chloroplast genome is 151,935 base pairs long, comprising two inverted repeat regions, a large single copy region, and a small single copy region. This chloroplast genome contains 128 genes, including 8 rRNA-coding genes, 37 tRNA-coding genes, 4 pseudogenes, and 83 protein-coding genes. When compared to the chloroplast genome of the invasive weed A. philoxeroides and other Amaranthaceae species, we observed significant variations in the ccsA, ycf1, and ycf2 regions in the A. sessilis chloroplast genome. Moreover, two genes, ccsA and accD, were found to be undergoing rapid evolution due to positive selection pressure. The phylogenetic trees were constructed for the Amaranthaceae family, estimating the time of independent species formation between A. philoxeroides and A. sessilis to be approximately 3.5186-8.8242 million years ago. These findings provide a foundation for understanding the population variation within invasive species among the Alternanthera genus.

Comparative genomic studies on the TGF-β superfamily in blue whale.

TGF-β supergene family has a wide range of physiological functions including cell adhesion, motility, proliferation, apoptosis, and differentiation. We systematically analyzed and characterized the TGF-β gene superfamily from the whole blue whale (Balaenoptera musculus) genome, using comparative genomic and evolutionary analysis. We identified 30 TGF-β genes and were split into two subgroups, BMP-like and TGF-like. All TGF-β proteins demonstrating a basic nature, with the exception of BMP1, BMP2, BMP10, GDF2, MSTN, and NODAL modulator, had acidic characteristics. All the blue whale (B. musculus) TGF-β proteins, excluding BMP1, are thermostable based on aliphatic index. The instability index showed all proteins except the NODAL modulator was unstable. TGF-β proteins showed a hydrophilic character, with the exception of GDF1 and INHBC. Moreover, all the detected TGF-β genes showed evolutionary conserved nature. A segmental duplication was indicated by TGF-β gene family, and the Ka/Ks ratio showed that the duplicated gene pairs were subjected to selection pressure, indicating both purifying and positive selection pressure. Two possible recombination breakpoints were also predicted. This study provides insights into the genetic characterization and evolutionary aspects of the TGF-β superfamily in blue whales (B. musculus).