Search

Background: The genus Robinsonia DC. (tribe Senecioneae, Asteraceae) endemic to the Juan Fernández Islands in Chile is one of the most conspicuous insular plant groups in the world. Unlike typical herbaceous Asteraceae plants, these plants demonstrate spectacular and unusual rosette tree growth forms as shown by the alpine giant senecios (genus Dendrosenecio, tribe Senecioneae) endemic to the East African mountains. However, monophyly of the genus and phylogenetic relationships among species of Robinsonia as well as their plastome evolution remain elusive. This study aims to explore their phylogeny, species diversification, and molecular evolution based on the complete plastome sequences in the context of adaptive radiation on oceanic islands. Results: The insular Robinsonia plastomes are highly conserved in their structures and organization of contents. Five divergence hotspots as potential chloroplast markers and five positively selected coding genes (accD, ndhF, rpoA, ycf1, and ycf2) are identified. Robinsonia plastomes has an overall nucleotide diversity higher than that of the sky island Dendrosenecio, but much lower than herbaceous Senecio. Phylogenetic analysis demonstrates the monophyly of Robinsonia and identifies two major infrageneric lineages. Both Robinsonia and Dendrosenecio are deeply nested within large genus Senecio. Conclusions: While plastid genomes of Robinsonia are highly conserved, their sequences strongly demonstrated the monophyly of the genus and inferred robust interspecific relationships, including herbaceous Senecio and woody Dendrosenecio. Different sets of positively selected chloroplast genes, five for Robinsonia and two for Dendrosenecio, may play an important role in the adaptation strategies of these fascinating woody species in insular and continental sky island habitats. Overall phylogenetic positions and sister lineages of Robinsonia and Dendrosenecio require additional study based on broader sampling of Senecio. [ABSTRACT FROM AUTHOR]

Background: Lithospermeae is the largest tribe within Boraginaceae. The tribe has been the focus of multiple phylogenetic studies over the last 15 years, with most focused on one genus or a few genera. In the present study, we newly sequenced 69 species of Lithospermeae and relatives to analyze the phylogenomic relationships among its members as well as the evolution of the plastid genome. Results: The phylogeny of Lithospermeae resolved from the plastid genome and nrDNA cistron is generally congruent with prior studies, but is better resolved and supported. Increasing character sampling across the plastid genome results in gradually more similar trees to that from the entire plastid genome. Overall, plastid genome structure was quite consistent across Lithospermeae. Codon Usage Bias (CUB) analyses demonstrate that across Lithospermeae plastid genomes were rich in AT and poor in GC. Mutation may play a greater role than selection across the plastid genome of Lithospermeae. The present study is the first to highlight the CUB characteristics of Lithospermeae species, which can help elucidate the mechanisms underlying patterns of molecular evolution and improve the expression levels of exogenous genes by codon optimization. Conclusions: This study provides a comprehensive phylogenomic analysis of Lithospermeae, significantly enhancing our understanding of the phylogenetic relationships and plastid genome evolution within this largest tribe of Boraginaceae. By utilizing an expanded genomic sampling approach, we have achieved increased resolution and support among the evolutionary relationships of the tribe, in line with but improving upon previous studies. The analyses of plastid genome structure revealed consistency across Lithospermeae, with a notable CUB. This study marks the first investigation into the CUB of Lithospermeae species and sets the stage for further research on the molecular evolution of plastid genomes across Boraginaceae. [ABSTRACT FROM AUTHOR]

Background: Mitochondrial genes and nuclear genes cooperate closely to maintain the functions of mitochondria, especially in the oxidative phosphorylation (OXPHOS) pathway. However, mitochondrial genes among arthropod lineages have dramatic evolutionary rate differences. Haplodiploid arthropods often show fast-evolving mitochondrial genes. One hypothesis predicts that the small effective population size of haplodiploid species could enhance the effect of genetic drift leading to higher substitution rates in mitochondrial and nuclear genes. Alternatively, positive selection or compensatory changes in nuclear OXPHOS genes could lead to the fast-evolving mitochondrial genes. However, due to the limited number of arthropod genomes, the rates of evolution for nuclear genes in haplodiploid species, besides hymenopterans, are largely unknown. To test these hypotheses, we used data from 76 arthropod genomes, including 5 independently evolved haplodiploid lineages, to estimate the evolutionary rates and patterns of gene family turnover of mitochondrial and nuclear genes. Results: We show that five haplodiploid lineages tested here have fast-evolving mitochondrial genes and fast-evolving nuclear genes related to mitochondrial functions, while nuclear genes not related to mitochondrion showed no significant evolutionary rate differences. Among hymenopterans, bees and ants show faster rates of molecular evolution in mitochondrial genes and mitochondrion-related nuclear genes than sawflies and wasps. With genome data, we also find gene family expansions and contractions in mitochondrion-related genes of bees and ants. Conclusions: Our results reject the small population size hypothesis in haplodiploid species. A combination of positive selection and compensatory changes could lead to the observed patterns in haplodiploid species. The elevated evolutionary rates in OXPHOS complex 2 genes of bees and ants suggest a unique evolutionary history of social hymenopterans. [ABSTRACT FROM AUTHOR]

Background: The tribe Ampelopsideae plants are important garden plants with both medicinal and ornamental values. The study of codon usage bias (CUB) facilitates a deeper comprehension of the molecular genetic evolution of species and their adaptive strategies. The joint analysis of CUB in chloroplast genomes (cpDNA) offers valuable insights for in-depth research on molecular genetic evolution, biological resource conservation, and elite breeding within this plant family. Results: The base composition and codon usage preferences of the eighteen chloroplast genomes were highly similar, with the GC content of bases at all positions of their codons being less than 50%. This indicates that they preferred A/T bases. Their effective codon numbers were all in the range of 35–61, which indicates that the codon preferences of the chloroplast genomes of the 18 Ampelopsideae plants were relatively weak. A series of analyses indicated that the codon preference of the chloroplast genomes of the 18 Ampelopsideae plants was influenced by a combination of multiple factors, with natural selection being the primary influence. The clustering tree generated based on the relative usage of synonymous codons is consistent with some of the results obtained from the phylogenetic tree of chloroplast genomes, which indicates that the clustering tree based on the relative usage of synonymous codons can be an important supplement to the results of the sequence-based phylogenetic analysis. Eventually, 10 shared best codons were screened on the basis of the chloroplast genomes of 18 species. Conclusion: The codon preferences of the chloroplast genome in Ampelopsideae plants are relatively weak and are primarily influenced by natural selection. The codon composition of the chloroplast genomes of the eighteen Ampelopsideae plants and their usage preferences were sufficiently similar to demonstrate that the chloroplast genomes of Ampelopsideae plants are highly conserved. This study provides a scientific basis for the genetic evolution of chloroplast genes in Ampelopsideae species and their suitable strategies. [ABSTRACT FROM AUTHOR]

To reveal the molecular function of elongation family of very long chain fatty acids(ELO) protein in Cyrtotrachelus buqueti, we have identified 15 ELO proteins from C.buqueti genome. 15 CbuELO proteins were located on four chromosomes. Their isoelectric points ranged from 9.22 to 9.68, and they were alkaline. These CbuELO proteins were stable and hydrophobic. CbuELO proteins had transmembrane movement, and had multiple phosphorylation sites. The secondary structure of CbuELO proteins was mainly α-helix. A total of 10 conserved motifs were identified in CbuELO protein family. Phylogenetic analysis showed that molecular evolutionary relationships of ELO protein family between C. buqueti and Tribolium castaneum was the closest. Developmental transcriptome analysis indicated that CbuELO10, CbuELO13 and CbuELO02 genes were key enzyme genes that determine the synthesis of very long chain fatty acids in pupae and eggs, CbuELO6 and CbuELO7 were that in the male, and CbuELO8 and CbuELO11 were that in the larva. Transcriptome analysis under different temperature conditions indicated that CbuELO1, CbuELO5, CbuELO12 and CbuELO14 participated in regulating temperature stress responses. Transcriptome analysis at different feeding times showed CbuELO12 gene expression level in all feeding time periods was significant downregulation. The qRT-PCR experiment verified expression level changes of CbuELO gene family under different temperature and feeding time conditions. Protein-protein interaction analysis showed that 9 CbuELO proteins were related to each other, CbuELO1, CbuELO4 and CbuELO12 had more than one interaction relationship. These results lay a theoretical foundation for further studying its molecular function during growth and development of C. buqueti. [ABSTRACT FROM AUTHOR]

DNA binding with one finger(Dof) gene family is a class of transcription factors which play an important role on plant growth and development. Genome-wide identification results indicated that there were 45 Dof genes(ColDof) in C.oleifera genome. All 45 ColDof proteins were non-transmembrane and non-secretory proteins. Phosphorylation site analysis showed that biological function of ColDof proteins were mainly realized by phosphorylation at serine (Ser) site. The secondary structure of 44 ColDof proteins was dominated by random coil, and only one ColDof protein was dominated by α-helix. ColDof genes' promoter region contained a variety of cis-acting elements, including light responsive regulators, gibberellin responsive regulators, abscisic acid responsive regulators, auxin responsive regulators and drought induction responsive regulators. The SSR sites analysis showed that the proportion of single nucleotide repeats and the frequency of A/T in ColDof genes were the largest. Non-coding RNA analysis showed that 45 ColDof genes contained 232 miRNAs. Transcription factor binding sites of ColDof genes showed that ColDof genes had 5793 ERF binding sites, 4381 Dof binding sites, 2206 MYB binding sites, 3702 BCR-BPC binding sites. ColDof9, ColDof39 and ColDof44 were expected to have the most TFBSs. The collinearity analysis showed that there were 40 colinear locis between ColDof proteins and AtDof proteins. Phylogenetic analysis showed that ColDof gene family was most closely related to that of Camellia sinensis var. sinensis cv.Biyun and Camellia lanceoleosa. Protein-protein interaction analysis showed that ColDof34, ColDof20, ColDof28, ColDof35, ColDof42 and ColDof26 had the most protein interactions. The transcriptome analysis of C. oleifera seeds showed that 21 ColDof genes were involved in the growth and development process of C. oleifera seeds, and were expressed in 221 C. oleifera varieties. The results of qRT-PCR experiments treated with different concentrations NaCl and PEG6000 solutions indicated that ColDof1, ColDof2, ColDof14 and ColDof36 not only had significant molecular mechanisms for salt stress tolerance, but also significant molecular functions for drought stress tolerance in C. oleifera. The results of this study provide a reference for further understanding of the function of ColDof genes in C.oleifera. [ABSTRACT FROM AUTHOR]

Dof transcription factor family in Cyperus esculentus genome was identified and analyzed using bioinformatics. The analysis results revealed that C.esculentus genome contains 29 Dof genes (CesDof), all of which are located in the nucleus according to subcellular localization prediction. CesDof proteinrs have a range of 124 to 512 amino acids, with most being basic proteins. Their secondary structure was mainly irregular curl. The promoter sequence of CesDof genes contains cis-acting elements that respond to light, drought, hormones, low temperature, and circadian rhythm. Codon preference analysis showed that CesDof genes' codon preference ends in T/A. Collinearity analysis revealed that C.esculentus had three pairs of collinear CesDof genes. Additionally, there were 15 pairs of collinear genes between C.esculentus and Arabidopsis thaliana. The genetic relationship between C.esculentus and Rhynchospora pubera was found to be the closest. Phylogenetic tree analysis revealed that 29 CesDof genes of C.esculentus can be classified into 4 subgroups. Additionally, 144 miRNAs were predicted to target these CesDof genes. Furthermore, protein interaction analysis indicated that 15 Dof proteins in C.esculentus had interactions. The qRT-PCR verification results of drought stress and salt stress treatment experiments showed that most CesDof genes were involved in drought stress and salt stress responses, and the gene expression trends under drought stress and salt stress conditions were consistent. These results lay a theoretical foundation for further studying the molecular functions of Dof gene family in C.esculentus and its molecular mechanisms in regulating the life activities of C.esculentus. [ABSTRACT FROM AUTHOR]

Background: The C-type lectin family 18 (CLEC18) with lipid and glycan binding capabilities is important to metabolic regulation and innate immune responses against viral infection. However, human CLEC18 comprises three paralogous genes with highly similar sequences, making it challenging to distinguish genetic variations, expression patterns, and biological functions of individual CLEC18 paralogs. Additionally, the evolutionary relationship between human CLEC18 and its counterparts in other species remains unclear. Methods: To identify the sequence variation and evolutionary divergence of human CLEC18 paralogs, we conducted a comprehensive analysis using various resources, including human and non-human primate reference genome assemblies, human pangenome assemblies, and long-read-based whole-genome and -transcriptome sequencing datasets. Results: We uncovered paralogous sequence variants (PSVs) and polymorphic variants (PVs) of human CLEC18 proteins, and identified distinct signatures specific to each CLEC18 paralog. Furthermore, we unveiled a novel segmental duplication for human CLEC18A gene. By comparing CLEC18 across human and non-human primates, our research showed that the CLEC18 paralogy probably occurred in the common ancestor of human and closely related non-human primates, and the lipid-binding CAP/SCP/TAPS domain of CLEC18 is more diverse than its glycan-binding CTLD. Moreover, we found that certain amino acids alterations at variant positions are exclusive to human CLEC18 paralogs. Conclusions: Our findings offer a comprehensive profiling of the intricate variations and evolutionary characteristics of human CLEC18. [ABSTRACT FROM AUTHOR]

Background: PC (phytocyanin) is a class of copper-containing electron transfer proteins closely related to plant photosynthesis, abiotic stress responses growth and development in plants, and regulation of the expression of some flavonoids and phenylpropanoids, etc., however, compared with other plants, the PC gene family has not been systematically characterized in apple. Results: A total of 59 MdPC gene members unevenly distributed across 12 chromosomes were identified at the genome-wide level. The proteins of the MdPC family were classified into four subfamilies based on differences in copper binding sites and glycosylation sites: Apple Early nodulin-like proteins (MdENODLs), Apple Uclacyanin-like proteins (MdUCLs), Apple Stellacyanin-like proteins (MdSCLs), and Apple Plantacyanin-like proteins (MdPLCLs). Some MdPC members with similar gene structures and conserved motifs belong to the same group or subfamily. The internal collinearity analysis revealed 14 collinearity gene pairs among members of the apple MdPC gene. Interspecific collinearity analysis showed that apple had 31 and 35 homologous gene pairs with strawberry and grape, respectively. Selection pressure analysis indicated that the MdPC gene was under purifying selection. Prediction of protein interactions showed that MdPC family members interacted strongly with the Nad3 protein. GO annotation results indicated that the MdPC gene also regulated the biosynthesis of phenylpropanoids. Chip data analysis showed that (MdSCL3, MdSCL7 and MdENODL27) were highly expressed in mature fruits and peels. Many cis-regulatory elements related to light response, phytohormones, abiotic stresses and flavonoid biosynthetic genes regulation were identified 2000 bp upstream of the promoter of the MdPC gene, and qRT-PCR results showed that gene members in Group IV (MdSCL1/3, MdENODL27) were up-regulated at all five stages of apple coloring, but the highest expression was observed at the DAF13 (day after fruit bag removal) stage. The gene members in Group II (MdUCL9, MdPLCL3) showed down-regulated or lower expression in the first four stages of apple coloring but up-regulated and highest expression in the DAF 21 stage. Conclusion: Herein, one objective of these findings is to provide valuable information for understanding the structure, molecular evolution, and expression pattern of the MdPC gene, another major objective in this study was designed to lay the groundwork for further research on the molecular mechanism of PC gene regulation of apple fruit coloration. [ABSTRACT FROM AUTHOR]

Background: Insect Cytochrome P450 monooxygenase (CYPs or P450s) plays an important role in detoxifying insecticides, causing insect populations to develop resistance. However, the molecular functions of P450 gene family in Cyrtotrachelus buqueti genome are still lacking. Results: In this study, 71 CbuP450 genes have been identified. The amino acids length of CbuP450 proteins was between 183 aa ~ 1041 aa. They are proteins with transmembrane domains. The main component of their secondary structure is α-helix and random coils. Phylogenetic analysis showed that C. buqueti and Rhynchophorus ferrugineus were the most closely related. This gene family has 29 high-frequency codons, which tend to use A/T bases and A/T ending codons. Gene expression analysis showed that CbuP450_23 in the female adult may play an important role on high temperature resistance, and CbuP450_17 in the larval may play an important role on low temperature tolerance. CbuP450_10, CbuP450_17, CbuP450_23, CbuP450_10, CbuP450_16, CbuP450_20, CbuP450_23 and CbuP450_ 29 may be related to the regulation of bamboo fiber degradation genes in C. buqueti. Protein interaction analysis indicates that most CbuP450 proteins are mainly divided into three aspects: encoding the biosynthesis of ecdysteroids, participating in the decomposition of synthetic insecticides, metabolizing insect hormones, and participating in the detoxification of compounds. Conclusions: We systematically analyzed the gene and protein characteristics, gene expression, and protein interactions of CbuP450 gene family, revealing the key genes involved in the stress response of CbuP450 gene family in the resistance of C. buqueti to high or low temperature stress, and identified the key CbuP450 proteins involved in important life activity metabolism. These results provided a reference for further research on the function of P450 gene family in C. buqueti. [ABSTRACT FROM AUTHOR]

Background: Myxozoa is a class of cnidarian parasites that encompasses over 2,400 species. Phylogenetic relationships among myxozoans remain highly debated, owing to both a lack of informative morphological characters and a shortage of molecular markers. Mitochondrial (mt) genomes are a common marker in phylogeny and biogeography. However, only five complete myxozoan mt genomes have been sequenced: four belonging to two closely related genera, Enteromyxum and Kudoa, and one from the genus Myxobolus. Interestingly, while cytochrome oxidase genes could be identified in Enteromyxum and Kudoa, no such genes were found in Myxobolus squamalis, and another member of the Myxobolidae (Henneguya salminicola) was found to have lost its entire mt genome. To evaluate the utility of mt genomes to reconstruct myxozoan relationships and to understand if the loss of cytochrome oxidase genes is a characteristic of myxobolids, we sequenced the mt genome of five myxozoans (Myxobolus wulii, M. honghuensis, M. shantungensis, Thelohanellus kitauei and, Sphaeromyxa zaharoni) using Illumina and Oxford Nanopore platforms. Results: Unlike Enteromyxum, which possesses a partitioned mt genome, the five mt genomes were encoded on single circular chromosomes. An mt plasmid was found in M. wulii, as described previously in Kudoa iwatai. In all new myxozoan genomes, five protein-coding genes (cob, cox1, cox2, nad1, and nad5) and two rRNAs (rnl and rns) were recognized, but no tRNA. We found that Myxobolus and Thelohanellus species shared unidentified reading frames, supporting the view that these mt open reading frames are functional. Our phylogenetic reconstructions based on the five conserved mt genes agree with previously published trees based on the 18S rRNA gene. Conclusions: Our results suggest that the loss of cytochrome oxidase genes is not a characteristic of all myxobolids, the ancestral myxozoan mt genome was likely encoded on a single circular chromosome, and mt plasmids exist in a few lineages. Our findings indicate that myxozoan mt sequences are poor markers for reconstructing myxozoan phylogenetic relationships because of their fast-evolutionary rates and the abundance of repeated elements, which complicates assembly. [ABSTRACT FROM AUTHOR]

Background: Rhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism. Results: We assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening. Conclusion: Overall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies. [ABSTRACT FROM AUTHOR]

Mitochondrial genomes play important roles in studying genome evolution, phylogenetic analyses, and species identification. Amphipods (Class Malacostraca, Order Amphipoda) are one of the most ecologically diverse crustacean groups occurring in a diverse array of aquatic and terrestrial environments globally, from freshwater streams and lakes to groundwater aquifers and the deep sea, but we have a limited understanding of how habitat influences the molecular evolution of mitochondrial energy metabolism. Subterranean amphipods likely experience different evolutionary pressures on energy management compared to surface-dwelling taxa that generally encounter higher levels of predation and energy resources and live in more variable environments. In this study, we compared the mitogenomes, including the 13 protein-coding genes involved in the oxidative phosphorylation (OXPHOS) pathway, of surface and subterranean amphipods to uncover potentially different molecular signals of energy metabolism between surface and subterranean environments in this diverse crustacean group. We compared base composition, codon usage, gene order rearrangement, conducted comparative mitogenomic and phylogenomic analyses, and examined evolutionary signals of 35 amphipod mitogenomes representing 13 families, with an emphasis on Crangonyctidae. Mitogenome size, AT content, GC-skew, gene order, uncommon start codons, location of putative control region (CR), length of rrnL and intergenic spacers differed between surface and subterranean amphipods. Among crangonyctid amphipods, the spring-dwelling Crangonyx forbesi exhibited a unique gene order, a long nad5 locus, longer rrnL and rrnS loci, and unconventional start codons. Evidence of directional selection was detected in several protein-encoding genes of the OXPHOS pathway in the mitogenomes of surface amphipods, while a signal of purifying selection was more prominent in subterranean species, which is consistent with the hypothesis that the mitogenome of surface-adapted species has evolved in response to a more energy demanding environment compared to subterranean amphipods. Overall, gene order, locations of non-coding regions, and base-substitution rates points to habitat as an important factor influencing the evolution of amphipod mitogenomes. [ABSTRACT FROM AUTHOR]

Background: Evolutionarily conserved in plants, the enzyme D-myo-inositol-3-phosphate synthase (MIPS; EC 5.5.1.4) regulates the initial, rate-limiting reaction in the phytic acid biosynthetic pathway. They are reported to be transcriptional regulators involved in various physiological functions in the plants, growth, and biotic/abiotic stress responses. Even though the genomes of most legumes are fully sequenced and available, an all-inclusive study of the MIPS family members in legumes is still ongoing. Results: We found 24 MIPS genes in ten legumes: Arachis hypogea, Cicer arietinum, Cajanus cajan, Glycine max, Lablab purpureus, Medicago truncatula, Pisum sativum, Phaseolus vulgaris, Trifolium pratense and Vigna unguiculata. The total number of MIPS genes found in each species ranged from two to three. The MIPS genes were classified into five clades based on their evolutionary relationships with Arabidopsis genes. The structural patterns of intron/exon and the protein motifs that were conserved in each gene were highly group-specific. In legumes, MIPS genes were inconsistently distributed across their genomes. A comparison of genomes and gene sequences showed that this family was subjected to purifying selection and the gene expansion in MIPS family in legumes was mainly caused by segmental duplication. Through quantitative PCR, expression patterns of MIPS in response to various abiotic stresses, in the vegetative tissues of various legumes were studied. Expression pattern shows that MIPS genes control the development and differentiation of various organs, and have significant responses to salinity and drought stress. Conclusion: The MIPS genes in the genomes of legumes have been identified, characterized and their expression was analysed. The findings pave way for understanding their molecular functions and evolution, and lead to identify the putative MIPS genes associated with different cell and tissue development. [ABSTRACT FROM AUTHOR]

Background: MYB transcription factors are splay a vital role in plant biology, with previous research highlighting the significant impact of the R2R3-MYB-like transcription factor MYB5 on seed mucilage biosynthesis, trichome branching, and seed coat development. However, there is a dearth of studies investigating its role in the regulation of proanthocyanidin (PA) biosynthesis. Results: In this study, a total of 51 MYB5 homologous genes were identified across 31 species belonging to the Brassicaceae family, with particular emphasis on Brassica napus for subsequent investigation. Through phylogenetic analysis, these genes were categorized into four distinct subclasses. Protein sequence similarity and identity analysis demonstrated a high degree of conservation of MYB5 among species within the Brassicaceae family. Additionally, the examination of selection pressure revealed that MYB5 predominantly underwent purifying selection during its evolutionary history, as indicated by the Ka/Ks values of all MYB5 homologous gene pairs being less than one. Notably, we observed a higher rate of non-synonymous mutations in orthologous genes compared to paralogous genes, and the Ka/Ks value displayed a stronger correlation with Ka. In B. napus, an examination of expression patterns in five tissues revealed that MYB5 exhibited particularly high expression in the black seed coat. The findings from the WGCNA demonstrated a robust correlation between MYB5 and BAN(ANR) associated with PA biosynthesis in the black seed coat, providing further evidence of their close association and co-expression. Furthermore, the results obtained from of the analysis of protein interaction networks offer supplementary support for the proposition that MYB5 possesses the capability to interact with transcriptional regulatory proteins, specifically TT8 and TT2, alongside catalytic enzymes implicated in the synthesis of PAs, thereby making a contribution to the biosynthesis of PAs. These findings imply a plausible and significant correlation between the nuique expression pattern of MYB5 and the pigmentation of rapeseed coats. Nevertheless, additional research endeavors are imperative to authenticate and substantiate these findings. Conclusions: This study offers valuable insights into the genetic evolution of Brassicaceae plants, thereby serving as a significant reference for the genetic enhancement of Brassicaceae seed coat color. [ABSTRACT FROM AUTHOR]

Background: Costaceae, commonly known as the spiral ginger family, consists of approximately 120 species distributed in the tropical regions of South America, Africa, and Southeast Asia, of which some species have important ornamental, medicinal and ecological values. Previous studies on the phylogenetic and taxonomic of Costaceae by using nuclear internal transcribed spacer (ITS) and chloroplast genome fragments data had low resolutions. Additionally, the structures, variations and molecular evolution of complete chloroplast genomes in Costaceae still remain unclear. Herein, a total of 13 complete chloroplast genomes of Costaceae including 8 newly sequenced and 5 from the NCBI GenBank database, representing all three distribution regions of this family, were comprehensively analyzed for comparative genomics and phylogenetic relationships. Result: The 13 complete chloroplast genomes of Costaceae possessed typical quadripartite structures with lengths from 166,360 to 168,966 bp, comprising a large single copy (LSC, 90,802 − 92,189 bp), a small single copy (SSC, 18,363 − 20,124 bp) and a pair of inverted repeats (IRs, 27,982 − 29,203 bp). These genomes coded 111 − 113 different genes, including 79 protein-coding genes, 4 rRNA genes and 28 − 30 tRNAs genes. The gene orders, gene contents, amino acid frequencies and codon usage within Costaceae were highly conservative, but several variations in intron loss, long repeats, simple sequence repeats (SSRs) and gene expansion on the IR/SC boundaries were also found among these 13 genomes. Comparative genomics within Costaceae identified five highly divergent regions including ndhF, ycf1-D2, ccsA-ndhD, rps15-ycf1-D2 and rpl16-exon2-rpl16-exon1. Five combined DNA regions (ycf1-D2 + ndhF, ccsA-ndhD + rps15-ycf1-D2, rps15-ycf1-D2 + rpl16-exon2-rpl16-exon1, ccsA-ndhD + rpl16-exon2-rpl16-exon1, and ccsA-ndhD + rps15-ycf1-D2 + rpl16-exon2-rpl16-exon1) could be used as potential markers for future phylogenetic analyses and species identification in Costaceae. Positive selection was found in eight protein-coding genes, including cemA, clpP, ndhA, ndhF, petB, psbD, rps12 and ycf1. Maximum likelihood and Bayesian phylogenetic trees using chloroplast genome sequences consistently revealed identical tree topologies with high supports between species of Costaceae. Three clades were divided within Costaceae, including the Asian clade, Costus clade and South American clade. Tapeinochilos was a sister of Hellenia, and Parahellenia was a sister to the cluster of Tapeinochilos + Hellenia with strong support in the Asian clade. The results of molecular dating showed that the crown age of Costaceae was about 30.5 Mya (95% HPD: 14.9 − 49.3 Mya), and then started to diverge into the Costus clade and Asian clade around 23.8 Mya (95% HPD: 10.1 − 41.5 Mya). The Asian clade diverged into Hellenia and Parahellenia at approximately 10.7 Mya (95% HPD: 3.5 − 25.1 Mya). Conclusion: The complete chloroplast genomes can resolve the phylogenetic relationships of Costaceae and provide new insights into genome structures, variations and evolution. The identified DNA divergent regions would be useful for species identification and phylogenetic inference in Costaceae. [ABSTRACT FROM AUTHOR]

Background: GPR65 is a pH-sensing G-protein-coupled receptor that acts as a key innate immune checkpoint in the human tumor microenvironment, inhibiting the release of inflammatory factors and inducing significant upregulation of tissue repair genes. However, the expression pattern and function of GPR65 in osteosarcoma (OS) remain unclear. The purpose of this study was to investigate and elucidate the role of GPR65 in the microenvironment, proliferation and migration of OS. Methods: Retrospective RNA-seq data analysis was conducted in a cohort of 97 patients with OS data in the TAEGET database. In addition, single-cell sequencing data from six surgical specimens of human OS patients was used to analyze the molecular evolution process during OS genesis. Tissues chips and bioinformatics results were used to verify GPR65 expression level in OS. MTT, colony formation, EdU assay, wound healing, transwell assay and F-actin assay were utilized to analyze cell proliferation and invasion of OS cancer cells. RNA-seq was used to explore the potential mechanism of GPR65's role in OS. Results: GPR65 expression was significantly low in OS, and subgroup analysis found that younger OS patients, OS patients in metastatic status, and overall survival and progression free survival OS patients had lower GPR65 expression. From ScRNA-seq data of GSE162454, we found the expression of GPR65 is significantly positively correlated with CD4 + T cells CD8 + T cells and OS related macrophage infiltration. Verification experiment found that silencing the expression of GPR65 in osteosarcoma cells U2OS and HOS could promote the proliferation and invasion process, RNA-seq results showed that the role of GPR65 in OS cells was related to immune system, metabolism and signal transduction. Conclusion: The low expression of GPR65 in OS leads to high metastasis rate and poor prognosis in OS patients. The suppression of immune escape and inhibition of proliferation may be a key pathway for GPR65 to participate in the progression of OS. The current study strengthens the role of GPR65 in OS development and provides a potential biomarker for the prognosis of OS patients. [ABSTRACT FROM AUTHOR]

Background: Human mitochondrial heteroplasmy is an extensively investigated phenomenon in the context of medical diagnostics, forensic identification and molecular evolution. However, technical limitations of high-throughput sequencing hinder reliable determination of point heteroplasmies (PHPs) with minor allele frequencies (MAFs) within the noise threshold. Results: To investigate the PHP landscape at an MAF threshold down to 0.1%, we sequenced whole mitochondrial genomes at approximately 7.700x coverage, in multiple technical and biological replicates of longitudinal blood and buccal swab samples from 11 human donors (159 libraries in total). The results obtained by two independent sequencing platforms and bioinformatics pipelines indicate distinctive PHP patterns below and above the 1% MAF cut-off. We found a high inter-individual prevalence of low-level PHPs (MAF < 1%) at polymorphic positions of the mitochondrial DNA control region (CR), their tissue preference, and a tissue-specific minor allele linkage. We also established the position-dependent potential of minor allele expansion in PHPs, and short-term PHP instability in a mitotically active tissue. We demonstrate that the increase in sensitivity of PHP detection to minor allele frequencies below 1% within a robust experimental and analytical pipeline, provides new information with potential applicative value. Conclusions: Our findings reliably show different mutational loads between tissues at sub-1% allele frequencies, which may serve as an informative medical biomarker of time-dependent, tissue-specific mutational burden, or help discriminate forensically relevant tissues in a single person, close maternal relatives or unrelated individuals of similar phylogenetic background. [ABSTRACT FROM AUTHOR]

Background: The chloroplast genome (cp genome) is directly related to the study and analysis of molecular phylogeny and evolution of plants in the phylogenomics era. The cp genome, whereas, is highly plastic and exists as a heterogeneous mixture of sizes and physical conformations. It is advantageous to purify/enrich the circular chloroplast DNA (cpDNA) to reduce sequence complexity in cp genome research. Large-insert, ordered DNA libraries are more practical for genomics research than conventional, unordered ones. From this, a technique of constructing the ordered BAC library with the goal-insert cpDNA fragment is developed in this paper. Results: This novel in-situ-process technique will efficiently extract circular cpDNA from crops and construct a high-quality cpDNA library. The protocol combines the in-situ chloroplast lysis for the high-purity circular cpDNA with the in-situ substitute/ligation for the high-quality cpDNA library. Individually, a series of original buffers/solutions and optimized procedures for chloroplast lysis in-situ is different than bacterial lysis in-situ; the in-situ substitute/ligation that reacts on the MCE membrane is suitable for constructing the goal-insert, ordered cpDNA library while preventing the large-insert cpDNA fragment breakage. The goal-insert, ordered cpDNA library is arrayed on the microtiter plate by three colonies with the definite cpDNA fragment that are homologous-corresponds to the whole circular cpDNA of the chloroplast. Conclusion: The novel in-situ-process technique amply furtherance of research in genome-wide functional analysis and characterization of chloroplasts, such as genome sequencing, bioinformatics analysis, cloning, physical mapping, molecular phylogeny and evolution. [ABSTRACT FROM AUTHOR]

Background: The evolution of protein residues depends on the mutation rates of their encoding nucleotides, but it may also be affected by co-evolution with other residues. Chloroplasts function as environmental sensors, transforming fluctuating environmental signals into different physiological responses. We reasoned that habitat diversity may affect their rate and mode of evolution, which might be evidenced in the chloroplast genome. The Pteridaceae family of ferns occupy an unusually broad range of ecological niches, which provides an ideal system for analysis. Results: We conducted adaptive evolution and intra-molecular co-evolution analyses of Pteridaceae chloroplast DNAs (cpDNAs). The results indicate that the residues undergoing adaptive evolution and co-evolution were mostly independent, with only a few residues being simultaneously involved in both processes, and these overlapping residues tend to exhibit high mutations. Additionally, our data showed that Pteridaceae chloroplast genes are under purifying selection. Regardless of whether we grouped species by lineage (which corresponded with ecological niches), we determined that positively selected residues mainly target photosynthetic genes. Conclusions: Our work provides evidence for the adaptive evolution of Pteridaceae cpDNAs, especially photosynthetic genes, to different habitats and sheds light on the adaptive evolution and co-evolution of proteins. [ABSTRACT FROM AUTHOR]

Mechanosensitive (MS) ion channels provide efficient molecular mechanism for transducing mechanical forces into intracellular ion fluxes in all kingdoms of life. The mechanosensitive channel of small conductance (MscS) was one of the best-studied MS channels and its homologs (MSL, MscS-like) were widely distributed in cell-walled organisms. However, the origin, evolution and expansion of MSL proteins in plants are still not clear. Here, we identified more than 2100 MSL proteins from 176 plants and conducted a broad-scale phylogenetic analysis. The phylogenetic tree showed that plant MSL proteins were divided into three groups (I, II and III) prior to the emergence of chlorophytae algae, consistent with their specific subcellular localization. MSL proteins were distributed unevenly into each of plant species, and four parallel expansion was identified in angiosperms. In Brassicaceae, most MSL duplicates were derived by whole-genome duplication (WGD)/segmental duplications. Finally, a hypothetical evolutionary model of MSL proteins in plants was proposed based on phylogeny. Our studies illustrate the evolutionary history of the MSL proteins and provide a guide for future functional diversity analyses of these proteins in plants. [ABSTRACT FROM AUTHOR]

Background: Molluscan shell, composed of a diverse range of architectures and microstructures, is a classic model system to study the relationships between molecular evolution and biomineralized structure formation. The shells of oysters differ from those of other molluscs by possessing a novel microstructure, chalky calcite, which facilitates adaptation to the sessile lifestyle. However, the genetic basis and evolutionary origin of this adaptive innovation remain largely unexplored. Results: We report the first whole-genome assembly and shell proteomes of the Iwagaki oyster Crassostrea nippona. Multi-omic integrative analyses revealed that independently expanded and co-opted tyrosinase, peroxidase, TIMP genes may contribute to the chalky layer formation in oysters. Comparisons with other molluscan shell proteomes imply that von Willebrand factor type A and chitin-binding domains are basic members of molluscan biomineralization toolkit. Genome-wide identification and analyses of these two domains in 19 metazoans enabled us to propose that the well-known Pif may share a common origin in the last common ancestor of Bilateria. Furthermore, Pif and LamG3 genes acquire new genetic function for shell mineralization in bivalves and the chalky calcite formation in oysters likely through a combination of gene duplication and domain reorganization. Conclusions: The spatial expression of SMP genes in the mantle and molecular evolution of Pif are potentially involved in regulation of the chalky calcite deposition, thereby shaping the high plasticity of the oyster shell to adapt to a sessile lifestyle. This study further highlights neo-functionalization as a crucial mechanism for the diversification of shell mineralization and microstructures in molluscs, which may be applied more widely for studies on the evolution of metazoan biomineralization. [ABSTRACT FROM AUTHOR]

Background: Alfalfa mosaic virus (AMV) is an important virus affecting many vegetable crops in Egypt. In this study, virus isolates were collected from naturally infected potato, tomato, alfalfa and clover plants that showed suspected symptoms of AMV in different locations of Beheira and Alexandria governorates during the 2019–2020 growing season. The relative incidence of the virus ranged from 11–25% based on visual observations of symptoms and ELISA testing. A total of 41 samples were tested by ELISA using polyclonal antisera for AMV. Four AMV isolates collected from different host plants, named AM1 from potato, AM2 from tomato, AM3 from alfalfa and AM4 from alfalfa, were maintained on Nicotiana glutinosa plants for further characterization of AMV. Results: Electron micrographs of the purified viral preparation showed spheroidal particles with a diameter of 18 nm and three bacilliform particles with lengths of roughly 55, 68, and 110 nm and diameters identical to those of the spheroidal particles. The CP gene sequence comparisons of four AMV isolates (AM1, AM2, AM3 and AM4) showed the highest nucleotide identity of 99.7% with the Gomchi isolate from South Korea infecting Gomchi (Ligularia fischeri) plants. Phylogenetic analysis showed that the present isolates were grouped together into a distinct separate clade (GPI) along with the Gomchi isolate from South Korea. Similarly, the deduced amino acid sequence comparisons of Egyptian AMV isolates revealed that amino acids Q29, S30, T34, V92 and V175 were conserved among the Egyptian isolates in GPI. Conclusion: The present study found strong evolutionary evidence for the genetic diversity of AMV isolates by the identification of potential recombination events involving parents from GPI and GPII lineages. Additionally, the study found that Egyptian AMV isolates are genetically stable with low nucleotide diversity. Genetic analysis of the AMV population suggested that the AMV populations differ geographically, and AMV CP gene is under mild purifying selection. Furthermore, the study proposed that the Egyptian AMV population had common evolutionary ancestors with the Asian AMV population. Antioxidant enzymes activity was assessed on N. glutinosa plants in response to infection with each AMV isolate studied, and the results revealed that the enzyme activity varied. [ABSTRACT FROM AUTHOR]

Background: Paenarthrobacter nicotinovorans ATCC 49919 uses the pyridine-pathway to degrade nicotine and could provide a renewable source of precursors from nicotine-containing waste as well as a model for studying the molecular evolution of catabolic pathways and their spread by horizontal gene transfer via soil bacterial plasmids. Results: In the present study, the strain was sequenced using the Illumina NovaSeq 6000 and Oxford Nanopore Technology (ONT) MinION platforms. Following hybrid assembly with Unicycler, the complete genome sequence of the strain was obtained and used as reference for whole-genome-based phylogeny analyses. A total of 64 related genomes were analysed; five Arthrobacter strains showed both digital DNA-DNA hybridization and average nucleotide identity values over the species threshold when compared to P. nicotinovorans ATCC 49919. Five plasmids and two contigs belonging to Arthrobacter and Paenarthrobacter strains were shown to be virtually identical with the pAO1 plasmid of Paenarthrobacter nicotinovorans ATCC 49919. Moreover, a highly syntenic nic-genes cluster was identified on five plasmids, one contig and three chromosomes. The nic-genes cluster contains two major locally collinear blocks that appear to form a putative catabolic transposon. Although the origins of the nic-genes cluster and the putative transposon still elude us, we hypothesise here that the ATCC 49919 strain most probably evolved from Paenarthrobacter sp. YJN-D or a very closely related strain by acquiring the pAO1 megaplasmid and the nicotine degradation pathway. Conclusions: The data presented here offers another snapshot into the evolution of plasmids harboured by Arthrobacter and Paenarthrobacter species and their role in the spread of metabolic traits by horizontal gene transfer among related soil bacteria. [ABSTRACT FROM AUTHOR]

Background: The glyoxalase system includes glyoxalase I (GLXI), glyoxalase II (GLXII) and glyoxalase III (GLXIII), which are responsible for methylglyoxal (MG) detoxification and involved in abiotic stress responses such as drought, salinity and heavy metal. Results: In this study, a total of 620 GLX family genes were identified from 21 different plant species. The results of evolutionary analysis showed that GLX genes exist in all species from lower plants to higher plants, inferring that GLX genes might be important for plants, and GLXI and GLXII account for the majority. In addition, motif showed an expanding trend in the process of evolution. The analysis of cis-acting elements in 21 different plant species showed that the promoter region of the GLX genes were rich in phytohormones and biotic and abiotic stress-related elements, indicating that GLX genes can participate in a variety of life processes. In cotton, GLXs could be divided into two groups and most GLXIs distributed in group I, GLXIIs and GLXIIIs mainly belonged to group II, indicating that there are more similarities between GLXII and GLXIII in cotton evolution. The transcriptome data analysis and quantitative real-time PCR analysis (qRT-PCR) show that some members of GLX family would respond to high temperature treatment in G.hirsutum. The protein interaction network of GLXs in G.hirsutum implied that most members can participate in various life processes through protein interactions. Conclusions: The results elucidated the evolutionary history of GLX family genes in plants and lay the foundation for their functions analysis in cotton. [ABSTRACT FROM AUTHOR]

Background: Tea geometrid Ectropis grisescens (Geometridae: Lepidoptera), is one of the most destructive defoliators in tea plantations in China. The MAPK cascade is known to be an evolutionarily conserved signaling module, acting as pivotal cores of host–pathogen interactions. Although the chromosome-level reference genome of E. grisescens was published, the whole MAPK cascade gene family has not been fully identified yet, especially the expression patterns of MAPK cascade gene family members upon an ecological biopesticide, Metarhizium anisopliae, remains to be understood. Results: In this study, we have identified 19 MAPK cascade gene family members in E. grisescens, including 5 MAPKs, 4 MAP2Ks, 8 MAP3Ks, and 2 MAP4Ks. The molecular evolution characteristics of the whole Eg-MAPK cascade gene family, including gene structures, protein structural organization, chromosomal localization, orthologs construction and gene duplication, were systematically investigated. Our results showed that the members of Eg-MAPK cascade gene family were unevenly distributed in 13 chromosomes, and the clustered members in each group shared similar structures of the genes and proteins. Gene expression data revealed that MAPK cascade genes were expressed in all four developmental stages of E. grisescens and were fairly and evenly distributed in four different larva tissues. Importantly, most of the MAPK cascade genes were induced or constitutively expressed upon M. anisopliae infection. Conclusions: In summary, the present study was one of few studies on MAPK cascade gene in E. grisescens. The characterization and expression profiles of Eg-MAPK cascades genes might help develop new ecofriendly biological insecticides to protect tea trees. [ABSTRACT FROM AUTHOR]

Barley yellow dwarf virus (BYDV) has caused considerable losses in the global production of grain crops such as wheat, barley and maize. We investigated the phylodynamics of the virus by analysing 379 and 485 nucleotide sequences of the genes encoding the coat protein and movement protein, respectively. The maximum clade credibility tree indicated that BYDV-GAV and BYDV-MAV, BYDV-PAV and BYDV-PAS share the same evolutionary lineage, respectively. The diversification of BYDV arises from its adaptability to vector insects and geography. Bayesian phylogenetic analyses showed that the mean substitution rates of the coat and movement proteins of BYDV ranged from 8.327 × 10− 4 (4.700 × 10− 4–1.228 × 10− 3) and 8.671 × 10− 4 (6.143 × 10− 4–1.130 × 10− 3) substitutions/site/year, respectively. The time since the most recent common BYDV ancestor was 1434 (1040–1766) CE (Common Era). The Bayesian skyline plot (BSP) showed that the BYDV population experienced dramatic expansions approximately 8 years into the 21st century, followed by a dramatic decline in less than 15 years. Our phylogeographic analysis showed that the BYDV population originating in the United States was subsequently introduced to Europe, South America, Australia and Asia. The migration pathways of BYDV suggest that the global spread of BYDV is associated with human activities. [ABSTRACT FROM AUTHOR]

Background: Simulating DNA evolution has been done through coevolution-agnostic probabilistic frameworks for the past 3 decades. The most common implementation is by using the converse of the probabilistic approach used to infer phylogenies which, in the simplest form, simulates a single sequence at a time. However, biological systems are multi-genic, and gene products can affect each other's evolutionary paths through coevolution. These crucial evolutionary dynamics still remain to be simulated, and we believe that modelling them can lead to profound insights for comparative genomics. Results: Here we present CastNet, a genome evolution simulator that assumes each genome is a collection of genes with constantly evolving regulatory interactions in between them. The regulatory interactions produce a phenotype in the form of gene expression profiles, upon which fitness is calculated. A genetic algorithm is then used to evolve a population of such entities through a user-defined phylogeny. Importantly, the regulatory mutations are a response to sequence mutations, thus making a 1–1 relationship between the rate of evolution of sequences and of regulatory parameters. This is, to our knowledge, the first time the evolution of sequences and regulation have been explicitly linked in a simulation, despite there being a multitude of sequence evolution simulators, and a handful of models to simulate Gene Regulatory Network (GRN) evolution. In our test runs, we see a coevolutionary signal among genes that are active in the GRN, and neutral evolution in genes that are not included in the network, showing that selective pressures imposed on the regulatory output of the genes are reflected in their sequences. Conclusion: We believe that CastNet represents a substantial step for developing new tools to study genome evolution, and more broadly, coevolutionary webs and complex evolving systems. This simulator also provides a new framework to study molecular evolution where sequence coevolution has a leading role. [ABSTRACT FROM AUTHOR]

Background: Mycoheterotrophs, acquiring organic carbon and other nutrients from mycorrhizal fungi, have evolved repeatedly with substantial plastid genome (plastome) variations. To date, the fine-scale evolution of mycoheterotrophic plastomes at the intraspecific level is not well-characterized. A few studies have revealed unexpected plastome divergence among species complex members, possibly driven by various biotic/abiotic factors. To illustrate evolutionary mechanisms underlying such divergence, we analyzed plastome features and molecular evolution of 15 plastomes of Neottia listeroides complex from different forest habitats. Results: These 15 samples of Neottia listeroides complex split into three clades according to their habitats approximately 6 million years ago: Pine Clade, including ten samples from pine-broadleaf mixed forests, Fir Clade, including four samples from alpine fir forests and Fir-willow Clade with one sample. Compared with those of Pine Clade members, plastomes of Fir Clade members show smaller size and higher substitution rates. Plastome size, substitution rates, loss and retention of plastid-encoded genes are clade-specific. We propose to recognized six species in N. listeroides complex and slightly modify the path of plastome degradation. Conclusions: Our results provide insight into the evolutionary dynamics and discrepancy of closely related mycoheterotrophic orchid lineages at a high phylogenetic resolution. [ABSTRACT FROM AUTHOR]

Background: Influenza A viruses have undergone rapid evolution with virulent; however, complete and comprehensive data on gene evolution and amino acid variation of HA and NA in immunosuppressed patients was few. In this study, we analysed molecular epidemiology and evolution of influenza A viruses in immunosuppressed population, and immunocompetent population were used as controls. Methods: Full sequences of HA and NA of A(H1N1)pdm09 and A(H3N2) were acquired through reverse transcription-polymerase chain reaction (RT-PCR). HA and NA genes were sequenced using the Sanger method and phylogenetically analysed using ClustalW 2.10 and MEGA software version 11.0. Results: During the 2018–2020 influenza seasons, 54 immunosuppressed and 46 immunocompetent inpatients screened positive for influenza A viruses by using the quantitative real-time PCR (qRT-PCR) were enrolled. 27 immunosuppressed and 23 immunocompetent nasal swab or bronchoalveolar lavage fluid samples were randomly selected and sequenced using the Sanger method. A(H1N1)pdm09 were detected in 15 samples and the remaining 35 samples were A(H3N2) positive. By analyzing the HA and NA gene sequences of these virus strains, we found that all A(H1N1)pdm09 viruses shared high similarities to each other and the HA and NA genes of these viruses exclusively belonged to subclade 6B.1A.1. Some NA genes of A(H3N2) viruses were not in the same clade as those of A/Singapore/INFIMH-16-0019/2016 and A/Kansas/14/2017, which may have led to A(H3N2) being the dominant strain in the 2019–2020 influenza season. Both A(H1N1)pdm09 and A(H3N2) viruses showed similar evolutionary lineages patterns of HA and NA between immunosuppressed and immunocompetent patients. Compared with the vaccine strains, there were no statistically significant of HA and NA genes and amino acid sequences of influenza A viruses in immunosuppressed and immunocompetent patients. However, the oseltamivir resistance substitution of NA-H275Y and R292K have been observed in immunosuppressed patients. Conclusions: A(H1N1)pdm09 and A(H3N2) viruses showed similar evolutionary lineages patterns of HA and NA between immunosuppressed and immunocompetent patients. Both immunocompetent and immunosuppressed patients have some key substitutions, which should be of note monitored, especially those with potential to affect the viral antigen. [ABSTRACT FROM AUTHOR]

Background: The Salmonella enterica serovar Newport red onion outbreak of 2020 was the largest foodborne outbreak of Salmonella in over a decade. The epidemiological investigation suggested two farms as the likely source of contamination. However, single nucleotide polymorphism (SNP) analysis of the whole genome sequencing data showed that none of the Salmonella isolates collected from the farm regions were linked to the clinical isolates—preventing the use of phylogenetics in source identification. Here, we explored an alternative method for analyzing the whole genome sequencing data driven by the hypothesis that if the outbreak strain had come from the farm regions, then the clinical isolates would disproportionately contain plasmids found in isolates from the farm regions due to horizontal transfer. Results: SNP analysis confirmed that the clinical isolates formed a single, nearly-clonal clade with evidence for ancestry in California going back a decade. The clinical clade had a large core genome (4,399 genes) and a large and sparsely distributed accessory genome (2,577 genes, at least 64% on plasmids). At least 20 plasmid types occurred in the clinical clade, more than were found in the literature for Salmonella Newport. A small number of plasmids, 14 from 13 clinical isolates and 17 from 8 farm isolates, were found to be highly similar (> 95% identical)—indicating they might be related by horizontal transfer. Phylogenetic analysis was unable to determine the geographic origin, isolation source, or time of transfer of the plasmids, likely due to their promiscuous and transient nature. However, our resampling analysis suggested that observing a similar number and combination of highly similar plasmids in random samples of environmental Salmonella enterica within the NCBI Pathogen Detection database was unlikely, supporting a connection between the outbreak strain and the farms implicated by the epidemiological investigation. Conclusion: Horizontally transferred plasmids provided evidence for a connection between clinical isolates and the farms implicated as the source of the outbreak. Our case study suggests that such analyses might add a new dimension to source tracking investigations, but highlights the need for detailed and accurate metadata, more extensive environmental sampling, and a better understanding of plasmid molecular evolution. [ABSTRACT FROM AUTHOR]

Hand-foot-and-mouth disease (HFMD) is a viral illness commonly seen in young children under 5 years of age, characterized by typical manifestations such as oral herpes and rashes on the hands and feet. These symptoms typically resolve spontaneously within a few days without complications. Over the past two decades, our understanding of HFMD has greatly improved and it has received significant attention. A variety of research studies, including epidemiological, animal, and in vitro studies, suggest that the disease may be associated with potentially fatal neurological complications. These findings reveal clinical, epidemiological, pathological, and etiological characteristics that are quite different from initial understandings of the illness. It is important to note that HFMD has been linked to severe cardiopulmonary complications, as well as severe neurological sequelae that can be observed during follow-up. At present, there is no specific pharmaceutical intervention for HFMD. An inactivated Enterovirus A71 (EV-A71) vaccine that has been approved by the China Food and Drug Administration (CFDA) has been shown to provide a high level of protection against EV-A71-related HFMD. However, the simultaneous circulation of multiple pathogens and the evolution of the molecular epidemiology of infectious agents make interventions based solely on a single agent comparatively inadequate. Enteroviruses are highly contagious and have a predilection for the nervous system, particularly in child populations, which contributes to the ongoing outbreak. Given the substantial impact of HFMD around the world, this Review synthesizes the current knowledge of the virology, epidemiology, pathogenesis, therapy, sequelae, and vaccine development of HFMD to improve clinical practices and public health efforts. [ABSTRACT FROM AUTHOR]

Background: Thymic malignancies are a heterogeneous group of rare cancers for which systemic chemotherapy is the standard treatment in the setting of advanced, recurrent or refractory diseases. Both environmental and genetic risk factors have not been fully clarified and few target-specific drugs have been developed for thymic epithelial tumors. A major challenge in studying thymic epithelial tumors is the lack of preclinical models for translational studies. Main body: Starting from bioptic material of two consecutive recurrences of the same patient, we generated two patient-derived xenografts. The patient-derived xenografts models were characterized for histology by immunohistochemistry and mutations using next-generation sequencing. When compared to the original tumors resected from the patient, the two patient-derived xenografts had preserved morphology after the stain with hematoxylin and eosin, although there was a moderate degree of de-differentiation. From a molecular point of view, the two patient-derived xenografts maintained 74.3 and 61.8% of the mutations present in the human tumor of origin. Short conclusion: The newly generated patient-derived xenografts recapitulate both the molecular characteristics and the evolution of the thymoma it derives from well, allowing to address open questions for this rare cancer. [ABSTRACT FROM AUTHOR]

Background: Dicranostigma leptopodum (Maxim.) Fedde is a perennial herb with bright yellow flowers, well known as "Hongmao Cao" for its medicinal properties, and is an excellent early spring flower used in urban greening. However, its molecular genomic information remains largely unknown. Here, we sequenced and analyzed the chloroplast genome of D. leptopodum to discover its genome structure, organization, and phylogenomic position within the subfamily Papaveroideae. Results: The chloroplast genome size of D. leptopodum was 162,942 bp, and D. leptopodum exhibited a characteristic circular quadripartite structure, with a large single-copy (LSC) region (87,565 bp), a small single-copy (SSC) region (18,759 bp) and a pair of inverted repeat (IR) regions (28,309 bp). The D. leptopodum chloroplast genome encoded 113 genes, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The dynamics of the genome structures, genes, IR contraction and expansion, long repeats, and single sequence repeats exhibited similarities, with slight differences observed among the eight Papaveroideae species. In addition, seven interspace regions and three coding genes displayed highly variable divergence, signifying their potential to serve as molecular markers for phylogenetic and species identification studies. Molecular evolution analyses indicated that most of the genes were undergoing purifying selection. Phylogenetic analyses revealed that D. leptopodum formed a clade with the tribe Chelidonieae. Conclusions: Our study provides detailed information on the D. leptopodum chloroplast genome, expanding the available genomic resources that may be used for future evolution and genetic diversity studies. [ABSTRACT FROM AUTHOR]

Background: Foxtail millet (Setaria italica L.) is a millet species with high tolerance to stressful environments. Plant non-specific lipid transfer proteins (nsLTPs) are a kind of small, basic proteins involved in many biological processes. So far, the genome of S. italica has been fully sequenced, and a comprehensive understanding of the evolution and expression of the nsLTP family is still lacking in foxtail millet. Results: Forty-five nsLTP genes were identified in S. italica and clustered into 5 subfamilies except three single genes (SinsLTP38, SinsLTP7, and SinsLTP44). The proportion of SinsLTPs was different in each subfamily, and members within the same subgroup shared conserved exon–intron structures. Besides, 5 SinsLTP duplication events were investigated. Both tandem and segmental duplication contributed to nsLTP expansion in S. italica, and the duplicated SinsLTPs had mainly undergone purifying selection pressure, which suggested that the function of the duplicated SinsLTPs might not diverge much. Moreover, we identified the nsLTP members in 5 other monocots, and 41, 13, 10, 4, and 1 orthologous gene pairs were identified between S. italica and S. viridis, S. bicolor, Z. mays, O. sativa, and B. distachyon, respectively. The functional divergence within the nsLTP orthologous genes might be limited. In addition, the tissue-specific expression patterns of the SinsLTPs were investigated, and the expression profiles of the SinsLTPs in response to abiotic stress were analyzed, all the 10 selected SinsLTPs were responsive to drought, salt, and cold stress. Among the selected SinsLTPs, 2 paired duplicated genes shared almost equivalent expression profiles, suggesting that these duplicated genes might retain some essential functions during subsequent evolution. Conclusions: The present study provided the first systematic analysis for the phylogenetic classification, conserved domain and gene structure, expansion pattern, and expression profile of the nsLTP family in S. italica. These findings could pave a way for further comparative genomic and evolution analysis of nsLTP family in foxtail millet and related monocots, and lay the foundation for the functional analysis of the nsLTPs in S. italica. [ABSTRACT FROM AUTHOR]

Background: The study presents results of research on the evolution of plastid genomes in Stipa L. which is a large genus of the Poaceae family, comprising species diverse in terms of geographic distribution, growing under highly variated habitat conditions. Complete plastome sequences of 43 taxa from Stipeae and Ampelodesmae tribes were analyzed for the variability of the coding regions against the background of phylogenetic relationships within the genus Stipa. The research hypothesis put forward in our research was that some of coding regions are affected by a selection pressure differentiated between individual phylogenetic lines of Stipa, potentially reducing the phylogenetic informativeness of these CDS. The study aimed to answer the question, which genes evolve in Stipa most rapidly and what kind of changes in the properties of encoded amino acids this entails. Another goal of this research was to find out whether individual genes are affected by positive selection and finally, whether selective pressure is uniform within the genus or does it vary between particular evolutionary lines within the genus. Results: Results of our study proved the presence of selective pressure in 11 genes: ccsA, matK, ndhC, ndhF, ndhK, rbcL, rpoA rpoC1, rpoC2, rps8 and rps11. For the first time the effect of positive selection on the rps8, rps11, and ndhK genes was documented in grasses. The varied pace of evolution, different intensity and effects of selective pressure have been demonstrated between particular phylogenetic lines of the genus tested. Conclusions: Positive selection in plastid genome in Stipa mostly affects photosynthetic genes. The potential strongest adaptive pressure was observed in the rbcL gene, especially in the oldest evolutionary group comprising Central Asian high-mountain species: S. basiplumosa, S. klimesii, S. penicillata and S. purpurea, where adaptive pressure probably affected the amino acids directly related to the efficiency of CO2 assimilation. [ABSTRACT FROM AUTHOR]

Background: Functional complexity of the eukaryotic mitochondrial proteome is augmented by independent gene acquisition from bacteria since its endosymbiotic origins. Mammalian homologs of many ancestral mitochondrial proteins have uncharacterized catalytic activities. Recent forward genetic approaches attributed functions to proteins in established metabolic pathways, thereby limiting the possibility of identifying novel biology relevant to human disease. We undertook a bottom-up biochemistry approach to discern evolutionarily conserved mitochondrial proteins with catalytic potential. Results: Here, we identify a Parkinson-associated DJ-1/PARK7-like protein—glutamine amidotransferase-like class 1 domain-containing 3A (GATD3A), with bacterial evolutionary affinities although not from alphaproteobacteria. We demonstrate that GATD3A localizes to the mitochondrial matrix and functions as a deglycase. Through its amidolysis domain, GATD3A removes non-enzymatic chemical modifications produced during the Maillard reaction between dicarbonyls and amines of nucleotides and amino acids. GATD3A interacts with factors involved in mitochondrial mRNA processing and translation, suggestive of a role in maintaining integrity of important biomolecules through its deglycase activity. The loss of GATD3A in mice is associated with accumulation of advanced glycation end products (AGEs) and altered mitochondrial dynamics. Conclusions: An evolutionary perspective helped us prioritize a previously uncharacterized but predicted mitochondrial protein GATD3A, which mediates the removal of early glycation intermediates. GATD3A restricts the formation of AGEs in mitochondria and is a relevant target for diseases where AGE deposition is a pathological hallmark. [ABSTRACT FROM AUTHOR]

Background: After domestication, the evolution of phenotypically-varied sheep breeds has generated rich biodiversity. This wide phenotypic variation arises as a result of hidden genomic changes that range from a single nucleotide to several thousands of nucleotides. Thus, it is of interest and significance to reveal and understand the genomic changes underlying the phenotypic variation of sheep breeds in order to drive selection towards economically important traits. Review: Various traits contribute to the emergence of variation in sheep phenotypic characteristics, including coat color, horns, tail, wool, ears, udder, vertebrae, among others. The genes that determine most of these phenotypic traits have been investigated, which has generated knowledge regarding the genetic determinism of several agriculturally-relevant traits in sheep. In this review, we discuss the genomic knowledge that has emerged in the past few decades regarding the phenotypic traits in sheep, and our ultimate aim is to encourage its practical application in sheep breeding. In addition, in order to expand the current understanding of the sheep genome, we shed light on research gaps that require further investigation. Conclusions: Although significant research efforts have been conducted in the past few decades, several aspects of the sheep genome remain unexplored. For the full utilization of the current knowledge of the sheep genome, a wide practical application is still required in order to boost sheep productive performance and contribute to the generation of improved sheep breeds. The accumulated knowledge on the sheep genome will help advance and strengthen sheep breeding programs to face future challenges in the sector, such as climate change, global human population growth, and the increasing demand for products of animal origin. [ABSTRACT FROM AUTHOR]

Background: Saposhnikovia divaricata (Turcz.) Schischk. is a perennial herb whose dried roots are commonly used as a source of traditional medicines. To elucidate the organelle-genome-based phylogeny of Saposhnikovia species and the transfer of DNA between organelle genomes, we sequenced and characterised the mitochondrial genome (mitogenome) of S. divaricata. Results: The mitogenome of S. divaricata is a circular molecule of 293,897 bp. The nucleotide composition of the mitogenome is as follows: A, 27.73%; T, 27.03%; C, 22.39%; and G, 22.85. The entire gene content is 45.24%. A total of 31 protein-coding genes, 20 tRNAs and 4 rRNAs, including one pseudogene (rpl16), were annotated in the mitogenome. Phylogenetic analysis of the organelle genomes from S. divaricata and 10 related species produced congruent phylogenetic trees. Selection pressure analysis revealed that most of the mitochondrial genes of related species are highly conserved. Moreover, 2 and 46 RNA-editing sites were found in the chloroplast genome (cpgenome) and mitogenome protein-coding regions, respectively. Finally, a comparison of the cpgenome and the mitogenome assembled from the same dataset revealed 10 mitochondrial DNA fragments with sequences similar to those in the repeat regions of the cpgenome, suggesting that the repeat regions might be transferred into the mitogenome. Conclusions: In this study, we assembled and annotated the mitogenome of S. divaricata. This study provides valuable information on the taxonomic classification and molecular evolution of members of the family Apiaceae. [ABSTRACT FROM AUTHOR]

Objectives: The most recent research projects in scientific field of eutherian comparative genomics included intentions to sequence every extant eutherian species genome in foreseeable future, so that future revisions and updates of eutherian gene data sets were expected. Data description: Using 35 public eutherian reference genomic sequence assemblies and free available software, the eutherian comparative genomic analysis protocol RRID:SCR_014401 was published as guidance against potential genomic sequence errors. The protocol curated 14 eutherian third-party data gene data sets, including, in aggregate, 2615 complete coding sequences that were deposited in European Nucleotide Archive. The published eutherian gene collections were used in revisions and updates of eutherian gene data set classifications and nomenclatures that included gene annotations, phylogenetic analyses and protein molecular evolution analyses. [ABSTRACT FROM AUTHOR]

Background: Both the protein domains and transcript structures influence protein functional variation. The genomic location of both protein domains and transcript structural features can be described using the genomic coordinates of their encoded sequences. However, the coordinates of protein domains and transcriptional features often differ greatly, and it is difficult to view them in combination at the genome-wide level. In this paper, we describe the development of a new tool that allows users to visualize domains and transcript features together, using either built-in or uploaded genome datasets, and export publication-ready figures. Results: We developed a user-friendly, independent R package and Shiny web application named "VisProDom". VisProDom consists of a genome-wide database containing entire annotated transcripts merged with annotated protein domains from the Pfam database. The built-in dataset includes 82 files, which merge genome general feature format (GFF) annotations with rpsblast tabular outputs from protein sequence searches in the Pfam database. Multiple genomes can be simultaneously screened for protein domains or transcript names. VisProDom includes step-by-step introductions and clickable elements for ease of use. Conclusion: VisProDom can display hundreds of transcripts alongside protein domains and export figures in a publication-ready format. This makes it a valuable tool for molecular evolution and comparative genomics. [ABSTRACT FROM AUTHOR]

The World Health Organization (WHO) recommends surveillance of molecular markers of resistance to anti-malarial drugs. This is particularly important in the case of mass drug administration (MDA), which is endorsed by the WHO in some settings to combat malaria. Dihydroartemisinin-piperaquine (DHA-PPQ) is an artemisinin-based combination therapy which has been used in MDA. This review analyses the impact of MDA with DHA-PPQ on the evolution of molecular markers of drug resistance. The review is split into two parts. Section I reviews the current evidence for different molecular markers of resistance to DHA-PPQ. This includes an overview of the prevalence of these molecular markers in Plasmodium falciparum Whole Genome Sequence data from the MalariaGEN Pf3k project. Section II is a systematic literature review of the impact that MDA with DHA-PPQ has had on the evolution of molecular markers of resistance. This systematic review followed PRISMA guidelines. This review found that despite being a recognised surveillance tool by the WHO, the surveillance of molecular markers of resistance following MDA with DHA-PPQ was not commonly performed. Of the total 96 papers screened for eligibility in this review, only 20 analysed molecular markers of drug resistance. The molecular markers published were also not standardized. Overall, this warrants greater reporting of molecular marker prevalence following MDA implementation. This should include putative pfcrt mutations which have been found to convey resistance to DHA-PPQ in vitro. [ABSTRACT FROM AUTHOR]

Background: As a new epi-center of COVID-19 in Asia and a densely populated developing country, Indonesia is facing unprecedented challenges in public health. SARS-CoV-2 lineage B.1.466.2 was reported to be an indigenous dominant strain in Indonesia (once second only to the Delta variant). However, it remains unclear how this variant evolved and spread within such an archipelagic nation. Methods: For statistical description, the spatiotemporal distributions of the B.1.466.2 variant were plotted using the publicly accessible metadata in GISAID. A total of 1302 complete genome sequences of Indonesian B.1.466.2 strains with high coverage were downloaded from the GISAID's EpiCoV database on 28 August 2021. To determine the molecular evolutionary characteristics, we performed a time-scaled phylogenetic analysis using the maximum likelihood algorithm and called the single nucleotide variants taking the Wuhan-Hu-1 sequence as reference. To investigate the spatiotemporal transmission patterns, we estimated two dynamic parameters (effective population size and effective reproduction number) and reconstructed the phylogeography among different islands. Results: As of the end of August 2021, nearly 85% of the global SARS-CoV-2 lineage B.1.466.2 sequences (including the first one) were obtained from Indonesia. This variant was estimated to account for over 50% of Indonesia's daily infections during the period of March–May 2021. The time-scaled phylogeny suggested that SARS-CoV-2 lineage B.1.466.2 circulating in Indonesia might have originated from Java Island in mid-June 2020 and had evolved into two disproportional and distinct sub-lineages. High-frequency non-synonymous mutations were mostly found in the spike and NSP3; the S-D614G/N439K/P681R co-mutations were identified in its larger sub-lineage. The demographic history was inferred to have experienced four phases, with an exponential growth from October 2020 to February 2021. The effective reproduction number was estimated to have reached its peak (11.18) in late December 2020 and dropped to be less than one after early May 2021. The relevant phylogeography showed that Java and Sumatra might successively act as epi-centers and form a stable transmission loop. Additionally, several long-distance transmission links across seas were revealed. Conclusions: SARS-CoV-2 variants circulating in the tropical archipelago may follow unique patterns of evolution and transmission. Continuous, extensive and targeted genomic surveillance is essential. [ABSTRACT FROM AUTHOR]

Background: Components of the antioxidant defense system in Trypanosoma cruzi are potential targets for new drug development. Superoxide dismutases (SODs) constitute key components of antioxidant defense systems, removing excess superoxide anions by converting them into oxygen and hydrogen peroxide. The main goal of the present study was to investigate the genes coding for iron superoxide dismutase (FeSOD) in T. cruzi strains from an evolutionary perspective. Methods: In this study, molecular biology methods and phylogenetic studies were combined with drug assays. The FeSOD-A and FeSOD-B genes of 35 T. cruzi strains, belonging to six discrete typing units (Tcl–TcVI), from different hosts and geographical regions were amplified by PCR and sequenced using the Sanger method. Evolutionary trees were reconstructed based on Bayesian inference and maximum likelihood methods. Drugs that potentially interacted with T. cruzi FeSODs were identified and tested against the parasites. Results: Our results suggest that T. cruzi FeSOD types are members of distinct families. Gene copies of FeSOD-A (n = 2), FeSOD-B (n = 4) and FeSOD-C (n = 4) were identified in the genome of the T. cruzi reference clone CL Brener. Phylogenetic inference supported the presence of two functional variants of each FeSOD type across the T. cruzi strains. Phylogenetic trees revealed a monophyletic group of FeSOD genes of T. cruzi TcIV strains in both distinct genes. Altogether, our results support the hypothesis that gene duplication followed by divergence shaped the evolution of T. cruzi FeSODs. Two drugs, mangafodipir and polaprezinc, that potentially interact with T. cruzi FeSODs were identified and tested in vitro against amastigotes and trypomastigotes: mangafodipir had a low trypanocidal effect and polaprezinc was inactive. Conclusions: Our study contributes to a better understanding of the molecular biodiversity of T. cruzi FeSODs. Herein we provide a successful approach to the study of gene/protein families as potential drug targets. [ABSTRACT FROM AUTHOR]

Frameshift mutations have been considered of significant importance for the molecular evolution of proteins and their coding genes, while frameshift protein sequences encoded in the alternative reading frames of coding genes have been considered to be meaningless. However, functional frameshifts have been found widely existing. It was puzzling how a frameshift protein kept its structure and functionality while substantial changes occurred in its primary amino-acid sequence. This study shows that the similarities among frameshifts and wild types are higher than random similarities and are determined at different levels. Frameshift substitutions are more conservative than random substitutions in the standard genetic code (SGC). The frameshift substitutions score of SGC ranks in the top 2.0-3.5% of alternative genetic codes, showing that SGC is nearly optimal for frameshift tolerance. In many genes and certain genomes, frameshift-resistant codons and codon pairs appear more frequently than expected, suggesting that frameshift tolerance is achieved through not only the optimality of the genetic code but, more importantly, the further optimization of a specific gene or genome through the usages of codons/codon pairs, which sheds light on the role of frameshift mutations in molecular and genomic evolution. [ABSTRACT FROM AUTHOR]

Background: Diplozoidae are monogenean (Monogenea: Polyopisthocotylea) fish parasites characterised by a unique life history: two larvae permanently fuse into an X-shaped "Siamese" organism. Taxonomy and phylogeny of Diplozoidae and Polyopisthocotylea remain unresolved due to the unavailability of molecular markers with sufficiently high resolution. Mitogenomes may be a suitable candidate, but there are currently only 12 available for the Polyopisthocotylea (three for Diplozoidae). The only available study of diplozoid mitogenomes found unique base composition patterns and elevated evolution rates in comparison with other Monogenean mitogenomes. Methods: To further explore their evolution and generate molecular data for evolutionary studies, we sequenced the complete mitogenomes of two Diplozoidae species, Paradiplozoon homoion and Paradiplozoon yarkandense, and conducted a number of comparative mitogenomic analyses with other polyopisthocotyleans. Results: We found further evidence that mitogenomes of Diplozoidae evolve at a unique, elevated rate, which was reflected in their exceptionally long branches, large sizes, unique base composition, skews, and very low gene sequence similarity levels between the two newly sequenced species. They also exhibited remarkably large overlaps between some genes. Phylogenetic analysis of Polyopisthocotylea resolved all major taxa as monophyletic, and Mazocraeidea was split into two major clades: (Diplozoidae) + (all four remaining families: Diclidophoridae, Chauhaneidae, Mazocraeidae and Microcotylidae). It also provided further confirmation that the genus Paradiplozoon is paraphyletic and requires a taxonomic revision, so the two species may have to be renamed Indodiplozoon homoion and Diplozoon yarkandense comb. nov. Conclusions: Although our findings indicate that mitogenomes may be a promising tool for resolving the phylogeny of Polyopisthocotylea, elevated evolutionary rates of Diplozoidae may cause phylogenetic artefacts, so future studies should pay caution to this problem. Furthermore, as the reason for their elevated evolution remains unknown, Diplozoidae are a remarkably interesting lineage for other types of evolutionary mitogenomic studies. [ABSTRACT FROM AUTHOR]

Background: Paspalum notatum 'Flugge' is a diploid with 20 chromosomes (2n = 20) multi-purpose subtropical herb native to South America and has a high ecological significance. It is currently widely planted in tropical and subtropical regions. Despite the gene pool of P. notatum 'Flugge' being unearthed to a large extent in the past decade, no details about the genomic information of relevant species in Paspalum have been reported. In this study, the complete genome information of P. notatum was established and annotated through sequencing and de novo assembly of its genome. Results: The latest PacBio third-generation HiFi assembly and sequencing revealed that the genome size of P. notatum 'Flugge' is 541 M. The assembly result is the higher index among the genomes of the gramineous family published so far, with a contig N50 = 52Mbp, scaffold N50 = 49Mbp, and BUSCOs = 98.1%, accounting for 98.5% of the estimated genome. Genome annotation revealed 36,511 high-confidence gene models, thus providing an important resource for future molecular breeding and evolutionary research. A comparison of the genome annotation results of P. notatum 'Flugge' with other closely related species revealed that it had a close relationship with Zea mays but not close compared to Brachypodium distachyon, Setaria viridis, Oryza sativa, Puccinellia tenuiflora, Echinochloa crusgalli. An analysis of the expansion and contraction of gene families suggested that P. notatum 'Flugge' contains gene families associated with environmental resistance, increased reproductive ability, and molecular evolution, which explained its excellent agronomic traits. Conclusion: This study is the first to report the high-quality chromosome-scale-based genome of P. notatum 'Flugge' assembled using the latest PacBio third-generation HiFi sequencing reads. The study provides an excellent genetic resource bank for gramineous crops and invaluable perspectives regarding the evolution of gramineous plants. [ABSTRACT FROM AUTHOR]

The origin of SARS-CoV-2 is uncertain. Findings support a "bat origin" but results are not highly convincing. Studies found evidence that SARS-CoV-2 was around for many years before the pandemic outbreak. Evidence has been published that the progenitor of SARS-CoV-2 already had the capability to bind strongly to the human ACE2 receptor. This may be an indication that many other animal viruses are capable to jump to humans, having already affinity for a human receptor. This is quite worrying since current ecosystems' collapse brings people to high proximity with animals, increasing probabilities for random viral transitions. On the other hand, future adaptation of SARS-CoV-2 is of great concern. Virus-host interactions are complicated and unfortunately, we still do not have accurate tools for predicting viruses' future evolution. Viral adaptation is a multifactorial process and probably SARS-CoV-2 will not become soon, as we wish, a harmless infection. However, humanity is currently under the largest vaccination program and it's of great interest to see if vaccinations will change the evolutionary game against the virus. [ABSTRACT FROM AUTHOR]