Comprehensive Evolutionary Analysis Research Articles

Characterization, Evolution, Expression and Functional Divergence of the DMP Gene Family in Plants.

The DMP (DOMAIN OF UNKNOWN FUNCTION 679 membrane protein) domain, containing a family of membrane proteins specific to green plants, is involved in numerous biological functions including physiological processes, reproductive development and senescence in Arabidopsis, but their evolutionary relationship and biological function in most crops remains unknown. In this study, we scrutinized phylogenetic relationships, gene structure, conserved domains and motifs, promoter regions, gene loss/duplication events and expression patterns. Overall, 240 DMPs were identified and analyzed in 24 plant species selected from lower plants to angiosperms. Comprehensive evolutionary analysis revealed that these DMPs underwent purifying selection and could be divided into five groups (I-V). DMP gene structure showed that it may have undergone an intron loss event during evolution. The five DMP groups had the same domains, which were distinct from each other in terms of the number of DMPs; group III was the largest, closely followed by group V. The DMP promotor region with various cis-regulatory elements was predicted to have a potential role in development, hormone induction and abiotic stresses. Based on transcriptomic data, expression profiling revealed that DMPs were primarily expressed in reproductive organs and were moderately expressed in other tissues. Evolutionary analysis suggested that gene loss events occurred more frequently than gene duplication events among all groups. Overall, this genome-wide study elucidates the potential function of the DMP gene family in selected plant species, but further research is needed in many crops to validate their biological roles.

Evolution and post-transcriptional regulation insights of m6A writers, erasers, and readers in plant epitranscriptome.

As a dynamic and reversible post-transcriptional marker, N6-methyladenosine (m6A) plays an important role in the regulation of biological functions, which are mediated by m6A pathway components including writers (MT-A70, FIP37, VIR and HAKAI family), erasers (ALKBH family) and readers (YTH family). There is an urgent need for a comprehensive analysis of m6A pathway components across species at evolutionary levels. In this study, we identified 4062 m6A pathway components from 154 plant species including green algae, utilizing large-scale phylogenetic to explore their origin and evolution. We discovered that the copy number of writers was conserved among different plant lineages, with notable expansions in the ALKBH and YTH families. Synteny network analysis revealed conserved genomic contexts and lineage-specific transpositions. Furthermore, we used Direct RNA Sequencing (DRS) to reveal the Poly(A) length (PAL) and m6A ratio profiles in six angiosperms species, with a particular focus on the m6A pathway components. The ECT1/2-Poeaece4 sub-branches (YTH family) with unique genomic contexts exhibited significantly higher expression level than genes of other ECT1/2 poeaece sub-branches (ECT1/2-Poeaece1-3), accompanied by lower m6A modification and PAL. Besides, conserved m6A sites distributed in CDS and 3'UTR were detected in the ECT1/2-Poaceae4, and the dual-luciferase assay further demonstrated that these conserved m6A sites in the 3'UTR negatively regulated the expression of Firefly luciferase (LUC) gene. Finally, we developed transcription factor regulatory networks for m6A pathway components, using yeast one-hybrid assay demonstrated that PheBPC1 could interact with the PheECT1/2-5 promoter. Overall, this study presents a comprehensive evolutionary and functional analysis of m6A pathway components and their modifications in plants, providing a valuable resource for future functional analysis in this field.

Comprehensive Genomic and Evolutionary Analysis of Biofilm Matrix Clusters and Proteins in the Vibrio Genus.

Vibrio cholerae pathogens cause cholera, an acute diarrheal disease resulting in significant morbidity and mortality worldwide. Biofilm formation by V. cholerae enhances its survival in natural ecosystems and facilitates transmission during cholera outbreaks. Critical components of the biofilm matrix are the Vibrio polysaccharide (VPS) produced by the vps-1 and vps-2 gene clusters, and biofilm matrix proteins encoded in the rbm cluster. However, the biofilm matrix clusters and associated matrix proteins in other Vibrio species remain under investigated, and their evolutionary patterns are largely unknown. In this study, we systematically annotated the biofilm matrix clusters across 6,121 Vibrio genomes, revealing their distribution, diversity, and evolution. We found that biofilm matrix clusters not only exist in V. cholerae but also in phylogenetically distant Vibrio species. Additionally, vps-1 clusters tend to co-locate with rbmABC genes, while vps-2 clusters are often adjacent to rbmDEF genes in various Vibrio species, which helps explain the separation of these clusters by the rbm cluster in well-characterized V. cholerae strains. Evolutionary analysis of RbmC and Bap1 reveals that these two major biofilm matrix proteins are sequentially and structurally related and have undergone domain/modular alterations during their evolution. RbmC genes are more prevalent, while bap1 likely resulted from an ancient duplication event of rbmC and is only present in a major clade of species containing rbmC counterparts. Notably, a novel loop-less Bap1 variant, identified in two subspecies clades of V. cholerae, was found to be associated with altered biofilm formation and the loss of antibiotic efflux pumps and chemotaxis. Another rbm cluster gene, rbmB, involved in biofilm dispersal, was found to share a common ancestor with Vibrio prophage pectin lyase-like tail proteins, indicating its functional and evolutionary linkages to Vibriophage proteins. In summary, our findings establish a foundational understanding of the proteins and gene clusters that contribute to Vibrio biofilm formation from an evolutionary perspective across a broad taxonomic scale. This knowledge paves the way for future strategies aimed at engineering and controlling biofilms through genetic modification.

Computational evolution of social norms in well-mixed and group-structured populations

Models of indirect reciprocity study how social norms promote cooperation. In these models, cooperative individuals build up a positive reputation, which in turn helps them in their future interactions. The exact reputational benefits of cooperation depend on the norm in place, which may change over time. Previous research focused on the stability of social norms. Much less is known about how social norms initially evolve when competing with many others. A comprehensive evolutionary analysis, however, has been difficult. Even among the comparably simple space of so-called third-order norms, there are thousands of possibilities, each one inducing its own reputation dynamics. To address this challenge, we use large-scale computer simulations. We study the reputation dynamics of each third-order norm and all evolutionary transitions between them. In contrast to established work with only a handful of norms, we find that cooperation is hard to maintain in well-mixed populations. However, within group-structured populations, cooperation can emerge. The most successful norm in our simulations is particularly simple. It regards cooperation as universally positive, and defection as usually negative-unless defection takes the form of justified punishment. This research sheds light on the complex interplay of social norms, their induced reputation dynamics, and population structure.

A novel deep-benthic sea cucumber species of Benthodytes (Holothuroidea, Elasipodida, Psychropotidae) and its comprehensive mitochondrial genome sequencing and evolutionary analysis

BackgroundThe holothurians, commonly known as sea cucumbers, are marine organisms that possess significant dietary, nutritional, and medicinal value. However, the National Center for Biotechnology Information (NCBI) currently possesses only approximately 70 complete mitochondrial genome datasets of Holothurioidea, which poses limitations on conducting comprehensive research on their genetic resources and evolutionary patterns. In this study, a novel species of sea cucumber belonging to the genus Benthodytes, was discovered in the western Pacific Ocean. The genomic DNA of the novel sea cucumber was extracted, sequenced, assembled and subjected to thorough analysis.ResultsThe mtDNA of Benthodytes sp. Gxx-2023 (GenBank No. OR992091) exhibits a circular structure spanning 17,386 bp, comprising of 13 protein-coding genes (PCGs), 24 non-coding RNAs (2 rRNA genes and 22 tRNA genes), along with two putative control regions measuring 882 bp and 1153 bp, respectively. It exhibits a high AT% content and negative AT-skew, which distinguishing it from the majority of sea cucumbers in terms of environmental adaptability evolution. The mitochondrial gene homology between Gxx-2023 and other sea cucumbers is significantly low, with less than 91% similarity to Benthodytes marianensis, which exhibits the highest level of homology. Additionally, its homology with other sea cucumbers is below 80%. The mitogenome of this species exhibits a unique pattern in terms of start and stop codons, featuring only two types of start codons (ATG and ATT) and three types of stop codons including the incomplete T. Notably, the abundance of AT in the Second position of the codons surpasses that of the First and Third position. The gene arrangement of PCGs exhibits a relatively conserved pattern, while there exists substantial variability in tRNA. Evolutionary analysis revealed that it formed a distinct cluster with B. marianensis and exhibited relatively distant phylogenetic relationships with other sea cucumbers.ConclusionsThese findings contribute to the taxonomic diversity of sea cucumbers in the Elasipodida order, thereby holding significant implications for the conservation of biological genetic resources, evolutionary advancements, and the exploration of novel sea cucumber resources.

Genome-Wide Identification and Evolutionary and Expression Analyses of the Cyclin B Gene Family in Brassica napus.

Cyclin B (CYCB) is a regulatory subunit of cyclin-dependent kinase (CDK), the concentration of which fluctuates to regulate cell cycle progression. Extensive studies have been performed on cyclins in numerous species, yet the evolutionary relationships and biological functions of the CYCB family genes in Brassica napus remain unclear. In this study, we identified 299 CYCB genes in 11 B. napus accessions. Phylogenetic analysis suggests that CYCB genes could be divided into three subfamilies in angiosperms and that the CYCB3 subfamily members may be a newer group that evolved in eudicots. The expansion of BnaCYCB genes underwent segmental duplication and purifying selection in genomes, and a number of drought-responsive and light-responsive cis-elements were found in their promoter regions. Additionally, expression analysis revealed that BnaCYCBs were strongly expressed in the developing seed and silique pericarp, as confirmed by the obviously reduced seed size of the mutant cycb3;1 in Arabidopsis thaliana compared with Col-0. This study provides a comprehensive evolutionary analysis of CYCB genes as well as insight into the biological function of CYCB genes in B. napus.

Comprehensive molecular evolutionary analysis of small heat shock proteins in five diploid Gossypium species.

The small heat shock proteins (sHSPs) are important components in plant growth and development, and stress response. However, a systematical understanding of the sHSP family is yet to be reported in five diploid Gossypium species. In this study, 34 GlsHSPs, 36 GrsHSPs, 37 GtsHSPs, 37 GasHSPs, and 38 GhesHSPs were identified in Gossypium longicalyx, Gossypium raimondii, Gossypium turneri, Gossypium arboreum, and Gossypium herbaceum, respectively. These sHSP members can be clustered into 10 subfamilies. Different subfamilies had different member numbers, motif distributions, gene structures, gene duplication events, gene loss numbers, and cis-regulatory elements. Besides, the paleohexaploidization event in cotton ancestor led to expanding the sHSP members and it was also inherited by five diploid Gossypium species. After the cotton ancestor divergence, the sHSP members had the relatively conserved evolution in five diploid Gossypium species. The comprehensive evolutionary history of the sHSP family was revealed in five diploid Gossypium species. Furthermore, several GasHSPs and GhesHSPs were important candidates in plant growth and development, and stress response. These current findings can provide valuable information for the molecular evolution and further functional research of the sHSP family in cotton.

Lineage-specific gene duplication and expansion of DUF1216 gene family in Brassicaceae.

Proteins containing domain of unknown function (DUF) are prevalent in eukaryotic genome. The DUF1216 proteins possess a conserved DUF1216 domain resembling to the mediator protein of Arabidopsis RNA polymerase II transcriptional subunit-like protein. The DUF1216 family are specifically existed in Brassicaceae, however, no comprehensive evolutionary analysis of DUF1216 genes have been performed. We performed a first comprehensive genome-wide analysis of DUF1216 proteins in Brassicaceae. Totally 284 DUF1216 genes were identified in 27 Brassicaceae species and classified into four subfamilies on the basis of phylogenetic analysis. The analysis of gene structure and conserved motifs revealed that DUF1216 genes within the same subfamily exhibited similar intron/exon patterns and motif composition. The majority members of DUF1216 genes contain a signal peptide in the N-terminal, and the ninth position of the signal peptide in most DUF1216 is cysteine. Synteny analysis revealed that segmental duplication is a major mechanism for expanding of DUF1216 genes in Brassica oleracea, Brassica juncea, Brassica napus, Lepidium meyneii, and Brassica carinata, while in Arabidopsis thaliana and Capsella rubella, tandem duplication plays a major role in the expansion of the DUF1216 gene family. The analysis of Ka/Ks (non-synonymous substitution rate/synonymous substitution rate) ratios for DUF1216 paralogous indicated that most of gene pairs underwent purifying selection. DUF1216 genes displayed a specifically high expression in reproductive tissues in most Brassicaceae species, while its expression in Brassica juncea was specifically high in root. Our studies offered new insights into the phylogenetic relationships, gene structures and expressional patterns of DUF1216 members in Brassicaceae, which provides a foundation for future functional analysis.

Genome-wide identification, characterization and evolutionary analysis of betaine aldehyde dehydrogenase (BADH), mitogen-activated protein kinase (MAPK) and sodium/hydrogen exchanger (NHX) genes in maize (Zea mays) under salt stress

Betaine aldehyde dehydrogenase (BADH), mitogen-activated protein kinase (MAPK) and sodium/hydrogen exchanger (NHX) play important roles in the response to salt stress. This is the first study to identify the BADH and NHX genes in maize (Zea mays) via genome-wide analysis. The qRT‒PCR results indicated that ZmNHX was upregulated by 4.38-fold, while a significant difference was not observed in ZmBADH or ZmMAPK, with fold changes of 0.96 and 1.06, respectively, under salinity stress. Genome-wide analysis revealed 8 ZmBADH, 19 ZmMAPK and 11 ZmNHX proteins in Z. mays. Domain analysis confirmed the presence of the aldehyde dehydrogenase superfamily (ALDH-SF), protein kinase and Na_H_Exchanger domains in the ZmBADH, ZmMAPK and ZmNHX proteins, respectively. Motif analysis indicated that the phylogenetic relationships were similar to the conserved motif distributions within the clade. The Ka/Ks ratio indicated that the ZmBADH, ZmMAPK and ZmNHX genes were influenced primarily by purifying selection. This study provides comprehensive identification, characterization, and evolutionary analysis for a better understanding of the ZmBADH, ZmMAPK and ZmNHX genes in maize.

Caspase-1 activates gasdermin A in non-mammals.

Gasdermins oligomerize to form pores in the cell membrane, causing regulated lytic cell death called pyroptosis. Mammals encode five gasdermins that can trigger pyroptosis: GSDMA, B, C, D, and E. Caspase and granzyme proteases cleave the linker regions of and activate GSDMB, C, D, and E, but no endogenous activation pathways are yet known for GSDMA. Here, we perform a comprehensive evolutionary analysis of the gasdermin family. A gene duplication of GSDMA in the common ancestor of caecilian amphibians, reptiles, and birds gave rise to GSDMA-D in mammals. Uniquely in our tree, amphibian, reptile, and bird GSDMA group in a separate clade than mammal GSDMA. Remarkably, GSDMA in numerous bird species contain caspase-1 cleavage sites like YVAD or FASD in the linker. We show that GSDMA from birds, amphibians, and reptiles are all cleaved by caspase-1. Thus, GSDMA was originally cleaved by the host-encoded protease caspase-1. In mammals the caspase-1 cleavage site in GSDMA is disrupted; instead, a new protein, GSDMD, is the target of caspase-1. Mammal caspase-1 uses exosite interactions with the GSDMD C-terminal domain to confer the specificity of this interaction, whereas we show that bird caspase-1 uses a stereotypical tetrapeptide sequence to confer specificity for bird GSDMA. Our results reveal an evolutionarily stable association between caspase-1 and the gasdermin family, albeit a shifting one. Caspase-1 repeatedly changes its target gasdermin over evolutionary time at speciation junctures, initially cleaving GSDME in fish, then GSDMA in amphibians/reptiles/birds, and finally GSDMD in mammals.

Caspase-1 activates gasdermin A in non-mammals

Gasdermins oligomerize to form pores in the cell membrane, causing regulated lytic cell death called pyroptosis. Mammals encode five gasdermins that can trigger pyroptosis: GSDMA, B, C, D, and E. Caspase and granzyme proteases cleave the linker regions of and activate GSDMB, C, D, and E, but no endogenous activation pathways are yet known for GSDMA. Here, we perform a comprehensive evolutionary analysis of the gasdermin family. A gene duplication of GSDMA in the common ancestor of caecilian amphibians, reptiles, and birds gave rise to GSDMA–D in mammals. Uniquely in our tree, amphibian, reptile, and bird GSDMA group in a separate clade than mammal GSDMA. Remarkably, GSDMA in numerous bird species contain caspase-1 cleavage sites like YVAD or FASD in the linker. We show that GSDMA from birds, amphibians, and reptiles are all cleaved by caspase-1. Thus, GSDMA was originally cleaved by the host-encoded protease caspase-1. In mammals the caspase-1 cleavage site in GSDMA is disrupted; instead, a new protein, GSDMD, is the target of caspase-1. Mammal caspase-1 uses exosite interactions with the GSDMD C-terminal domain to confer the specificity of this interaction, whereas we show that bird caspase-1 uses a stereotypical tetrapeptide sequence to confer specificity for bird GSDMA. Our results reveal an evolutionarily stable association between caspase-1 and the gasdermin family, albeit a shifting one. Caspase-1 repeatedly changes its target gasdermin over evolutionary time at speciation junctures, initially cleaving GSDME in fish, then GSDMA in amphibians/reptiles/birds, and finally GSDMD in mammals.

Comprehensive virulence profiling and evolutionary analysis of specificity determinants in Staphylococcus aureus two-component systems.

In the Staphylococcus aureus genome, a set of highly conserved two-component systems (TCSs) composed of histidine kinases (HKs) and their cognate response regulators (RRs) sense and respond to environmental stimuli, which drive the adaptation of the bacteria. This study investigates the complex interplay between TCSs in S. aureus USA300, a predominant methicillin-resistant S. aureus strain, revealing shared and unique virulence regulatory pathways and genetic variations mediating signal specificity within TCSs. Using TCS-related mutants from the Nebraska Transposon Mutant Library, we analyzed the effects of inactivated TCS HKs and RRs on the production of various virulence factors, in vitro infection abilities, and adhesion assays. We found that the TCSs' influence on virulence determinants was not associated with their phylogenetic relationship, indicating divergent functional evolution. Using the co-crystallized structure of the DesK-DesR from Bacillus subtilis and the modeled structures of the four NarL TCSs in S. aureus, we identified interacting residues, revealing specificity determinants and conservation within the same TCS, even from different strain backgrounds. The interacting residues were highly conserved within strains but varied between species due to selection pressures and the coevolution of cognate pairs. This study unveils the complex interplay and divergent functional evolution of TCSs, highlighting their potential for future experimental exploration of phosphotransfer between cognate and non-cognate recombinant HK and RRs.IMPORTANCEGiven the widespread conservation of two-component systems (TCSs) in bacteria and their pivotal role in regulating metabolic and virulence pathways, they present a compelling target for anti-microbial agents, especially in the face of rising multi-drug-resistant infections. Harnessing TCSs therapeutically necessitates a profound understanding of their evolutionary trajectory in signal transduction, as this underlies their unique or shared virulence regulatory pathways. Such insights are critical for effectively targeting TCS components, ensuring an optimized impact on bacterial virulence, and mitigating the risk of resistance emergence via the evolution of alternative pathways. Our research offers an in-depth exploration of virulence determinants controlled by TCSs in S. aureus, shedding light on the evolving specificity determinants that orchestrate interactions between their cognate pairs.

Comprehensive Evolutionary Analysis of the SMXL Gene Family in Rosaceae: Further Insights into Its Origin, Expansion, Diversification, and Role in Regulating Pear Branching.

SMXL genes constitute a conserved gene family that is ubiquitous in angiosperms and involved in regulating various plant processes, including branching, leaf elongation, and anthocyanin biosynthesis, but little is known about their molecular functions in pear branching. Here, we performed genome-wide identification and investigation of the SMXL genes in 16 angiosperms and analyzed their phylogenetics, structural features, conserved motifs, and expression patterns. In total, 121 SMXLs genes were identified and were classified into four groups. The number of non-redundant SMXL genes in each species varied from 3 (Amborella trichopoda Baill.) to 18 (Glycine max Merr.) and revealed clear gene expansion events over evolutionary history. All the SMXL genes showed conserved structures, containing no more than two introns. Three-dimensional protein structure prediction revealed distinct structures between but similar structures within groups. A quantitative real-time PCR analysis revealed different expressions of 10 SMXL genes from pear branching induced by fruit-thinning treatment. Overall, our study provides a comprehensive investigation of SMXL genes in the Rosaceae family, especially pear. The results offer a reference for understanding the evolutionary history of SMXL genes and provide excellent candidates for studying fruit tree branching regulation, and in facilitating pear pruning and planting strategies.

Genome sequencing and functional analysis of a multipurpose medicinal herb Tinospora cordifolia (Giloy).

Tinospora cordifolia (Willd.) Hook.f. & Thomson, also known as Giloy, is among the most important medicinal plants that have numerous therapeutic applications in human health due to the production of a diverse array of secondary metabolites. To gain genomic insights into the medicinal properties of T. cordifolia, the genome sequencing was carried out using 10× Genomics linked read and Nanopore long-read technologies. The draft genome assembly of T. cordifolia was comprised of 1.01 Gbp, which is the genome sequenced from the plant family Menispermaceae. We also performed the genome size estimation for T. cordifolia, which was found to be 1.13 Gbp. The deep sequencing of transcriptome from the leaf tissue was also performed. The genome and transcriptome assemblies were used to construct the gene set, resulting in 17,245 coding gene sequences. Further, the phylogenetic position of T. cordifolia was also positioned as basal eudicot by constructing a genome-wide phylogenetic tree using multiple species. Further, a comprehensive comparative evolutionary analysis of gene families contraction/expansion and multiple signatures of adaptive evolution was performed. The genes involved in benzyl iso-quinoline alkaloid, terpenoid, lignin and flavonoid biosynthesis pathways were found with signatures of adaptive evolution. These evolutionary adaptations in genes provide genomic insights into the presence of diverse medicinal properties of this plant. The genes involved in the common symbiosis signalling pathway associated with endosymbiosis (Arbuscular Mycorrhiza) were found to be adaptively evolved. The genes involved in adventitious root formation, peroxisome biogenesis, biosynthesis of phytohormones, and tolerance against abiotic and biotic stresses were also found to be adaptively evolved in T. cordifolia.

FEZF2 and AIRE1: An Evolutionary Trade-off in the Elimination of Auto-reactive T Cells in the Thymus.

Autoimmune Regulator 1 (AIRE1) and Forebrain Embryonic Zinc Finger-Like Protein 2 (FEZF2) play pivotal roles in orchestrating the expression of tissue-restricted antigens (TRA) to facilitate the elimination of autoreactive T cells. AIRE1's presence in the gonads of various vertebrates has raised questions about its potential involvement in gene expression control for germline cell selection. Nevertheless, the evolutionary history of these genes has remained enigmatic, as has the rationale behind their apparent redundancy in vertebrates. Furthermore, the origin of the elimination process itself has remained elusive. To shed light on these mysteries, we conducted a comprehensive evolutionary analysis employing a range of tools, including multiple sequence alignment, phylogenetic tree construction, ancestral sequence reconstruction, and positive selection assessment. Our investigations revealed intriguing insights. AIRE1 homologs emerged during the divergence of T cells in higher vertebrates, signifying its role in this context. Conversely, FEZF2 exhibited multiple homologs spanning invertebrates, lampreys, and higher vertebrates. Ancestral sequence reconstruction demonstrated distinct origins for AIRE1 and FEZF2, underscoring that their roles in regulating TRA have evolved through disparate pathways. Furthermore, it became evident that both FEZF2 and AIRE1 govern a diverse repertoire of genes, encompassing ancient and more recently diverged targets. Notably, FEZF2 demonstrates expression in both vertebrate and invertebrate embryos and germlines, accentuating its widespread role. Intriguingly, FEZF2 harbors motifs associated with autophagy, such as DKFPHP, SYSELWKSSL, and SYSEL, a process integral to cell selection in invertebrates. Our findings suggest that FEZF2 initially emerged to regulate self-elimination in the gonads of invertebrates. As organisms evolved toward greater complexity, AIRE1 likely emerged to complement FEZF2's role, participating in the regulation of cell selection for elimination in both gonads and the thymus. This dynamic interplay between AIRE1 and FEZF2 underscores their multifaceted contributions to TRA expression regulation across diverse evolutionary contexts.

Macroevolution of NLR genes in family Fabaceae provides evidence of clade specific expansion and contraction of NLRome in Vicioid clade

Whole genome duplication plays a significant role in plant genome evolution by providing raw materials that can be modified by natural or artificial selection. Nucleotide binding site leucine rich repeat receptor (NLRs) like other gene families are clusters of genes created by duplication and their size reflects the number of duplicated genes. NLR genes encode immune receptors that facilitate identification and binding of effector compounds produced by pathogen as a part of effector triggered immunity (ETI). The accurate identification and characterization of NLR genes substantially contributes to the repertoire of resistance and improves production. The ancestors of Fabaceae family have underwent whole genome duplication (WGD) approximately 58.5 million years ago. In this study, we focused on the subsequent effects of WGD on the evolution of NLR genes within the Vicioid clade, which consists of various legume crops such as chickpea, clover, alfalfa, and pea. The Vicioid clade is divided into three major tribes: Cicereae, Fabeae, and Trifolieae. The analysis of 22 species from the Vicioid clade revealed an overall contraction of the NLRome (the complete set of NLR genes) in members of the Cicereae and Fabeae tribes. This contraction aligns with previous observations that WGD events are often followed by diploidization, leading to a reduced number of duplicated genes. Contrary to this contraction trend, tribe Trifolieae have shown large scale expansion of NLRome irrespective to their genome size. Additionally, the primary diversification of relatively conserved NLR gene subclasses, including helper genes (CCR-NLR) and CCG10-NLR, was reported. Comprehensive evolutionary analysis suggests that NLRome expansion have occurred in during recent 1-6 Mya probably due to their higher substitutions per site per year in Trifolieae. We further hypothesized that this higher rate directed accelerated gene duplications after speciation from their common ancestor and later gene conversion and asymmetric recombination played an active role in subgroup diversification.

Genome of Phyllanthus emblica: the medicinal plant Amla with super antioxidant properties.

Phyllanthus emblica or Indian gooseberry, commonly known as amla, is an important medicinal horticultural plant used in traditional and modern medicines. It bears stone fruits with immense antioxidant properties due to being one of the richest natural sources of vitamin C and numerous flavonoids. This study presents the first genome sequencing of this species performed using 10x Genomics and Oxford Nanopore Technology. The draft genome assembly was 519 Mbp in size and consisted of 4,384 contigs, N50 of 597 Kbp, 98.4% BUSCO score, and 37,858 coding sequences. This study also reports the genome-wide phylogeny of this species with 26 other plant species that resolved the phylogenetic position of P. emblica. The presence of three ascorbate biosynthesis pathways including L-galactose, galacturonate, and myo-inositol pathways was confirmed in this genome. A comprehensive comparative evolutionary genomic analysis including gene family expansion/contraction and identification of multiple signatures of adaptive evolution provided evolutionary insights into ascorbate and flavonoid biosynthesis pathways and stone fruit formation through lignin biosynthesis. The availability of this genome will be beneficial for its horticultural, medicinal, dietary, and cosmetic applications and will also help in comparative genomics analysis studies.

Comprehensive genetic diversity and molecular evolutionary analysis of Theileria annulata isolates based on TAMS 1 gene

Molecular epidemiological studies related to the phylogenetic characterization of Theileria annulata are important in delineating the evolutionary history of the parasite. In the current study, the Theileria annulata (T. annulata) merozoite surface antigen 1 (TAMS 1) gene from 14 bovine isolates of T. annulata originating from semi-arid zone of northern India were amplified and sequenced. TAMS 1 gene sequences (n= 337) reported from 16 countries were subsequently analyzed for haplotype network along with genetic diversity. A total of five haplotypes out of the 14 sequenced isolates and 92 haplotypes out of 337 worldwide sequences are documented in this study. Phylogenetic and molecular evolutionary analyses based on TAMS 1 gene sequences showed that T. annulata is dissipated across different countries and numerous strains are closely linked, even though they belong to different geographical locations. The nucleotide homology between 14 isolates from northern India varied between 91.3 and 100%, whereas it was between 31.5 and 100% when sequences across the globe were compared. Haplotype 14 was recognized as most widely distributed haplotype, with 46 isolates circulating in 10 countries. Globally, negligible genetic distance (FST˂0.15) and very high gene flow (Nm˃1) was found in the five populations of the world (South Asia, East Asia, West Asia, Europe and Africa), supporting the absence of clearly defined subgroups in the phylogenetic analysis. Significant negative values of neutrality tests; Tajima's D (D) and Fu and Li's F (F) provided evidence for recent population expansion through positive selection of advantageous variations.

Isolation and genome characterization of Lloviu virus from Italian Schreibers’s bats

Lloviu cuevavirus (LLOV) was the first identified member of Filoviridae family outside the Ebola and Marburgvirus genera. A massive die-off of Schreibers’s bats (Miniopterus schreibersii) in the Iberian Peninsula in 2002 led to its initial discovery. Recent studies with recombinant and wild-type LLOV isolates confirmed the zoonotic nature of the virus in vitro. We examined bat samples from Italy for the presence of LLOV in an area outside of the currently known distribution range of the virus. We detected one positive sample from 2020, sequenced the complete coding region of the viral genome and established an infectious isolate of the virus. In addition, we performed the first comprehensive evolutionary analysis of the virus, using the Spanish, Hungarian and the Italian sequences. The most important achievement of this study is the establishment of an additional infectious LLOV isolate from a bat sample using the SuBK12-08 cells, demonstrating that this cell line is highly susceptible to LLOV infection and confirming the previous observation that these bats are effective hosts of the virus in nature. This result further strengthens the role of bats as the natural hosts for zoonotic filoviruses.

Whole-genome and dispersed duplication, including transposed duplication, jointly advance the evolution of TLP genes in seven representative Poaceae lineages

BackgroundIn the evolutionary study of gene families, exploring the duplication mechanisms of gene families helps researchers understand their evolutionary history. The tubby-like protein (TLP) family is essential for growth and development in plants and animals. Much research has been done on its function; however, limited information is available with regard to the evolution of the TLP gene family. Herein, we systematically investigated the evolution of TLP genes in seven representative Poaceae lineages.ResultsOur research showed that the evolution of TLP genes was influenced not only by whole-genome duplication (WGD) and dispersed duplication (DSD) but also by transposed duplication (TRD), which has been neglected in previous research. For TLP family size, we found an evolutionary pattern of progressive shrinking in the grass family. Furthermore, the evolution of the TLP gene family was at least affected by evolutionary driving forces such as duplication, purifying selection, and base mutations.ConclusionsThis study presents the first comprehensive evolutionary analysis of the TLP gene family in grasses. We demonstrated that the TLP gene family is also influenced by a transposed duplication mechanism. Several new insights into the evolution of the TLP gene family are presented. This work provides a good reference for studying gene evolution and the origin of duplication.