Endemic does not mean constant as SARS-CoV-2 continues to evolve

COVID-19 has become endemic, with dynamics that reflect the waning of immunity and re-exposure, by contrast to the epidemic phase driven by exposure in immunologically naïve populations. Endemic does not, however, mean constant. Further evolution of SARS-CoV-2, as well as changes in behavior and public health policy, continue to play a major role in the endemic load of disease and mortality. In this article, we analyze evolutionary models to explore the impact that a newly arising variant can have on the short-term and longer-term endemic load, characterizing how these impacts depend on the transmission and immunological properties of the variants. We describe how evolutionary changes in the virus will increase the endemic load most for a persistently immune-escape variant, by an intermediate amount for a more transmissible variant, and least for a transiently immune-escape variant. Balancing the tendency for evolution to favor variants that increase the endemic load, we explore the impact of vaccination strategies and non-pharmaceutical interventions that can counter these increases in the impact of disease. We end with some open questions about the future of COVID-19 as an endemic disease.

Codon switching of conserved Ser residues in coagulation and fibrinolytic proteases

BACKGROUND

Unique among all amino acids, Ser is encoded by two sets of codons, TCN and AGY (N=any nucleotide, Y=pyrimidine), that cannot interconvert through single nucleotide substitutions. Both codons are documented at the essential residues S195 and S214 within the active site of serine proteases. However, it is not known how the codons interconverted during evolution because replacement of S195 or S214 by other amino acids typically results in loss of activity.

OBJECTIVE

Characterize the prevalence of codon switching among essential and non-essential Ser residues in coagulation and fibrinolytic proteases from different vertebrate lineages.

METHODS

TCN and AGY codon usage was analyzed in >550 sequences.

RESULTS

Evolutionary pressure to preserve the codon of S195 is absolute, with no evidence of interconversion. Pressure to preserve the codon of S214 is also strong, but an AGY ↔ TCN interconversion is observed in factor VII-inactive and protein C from ray-finned fish. In both cases, the interconversion occurred in genes that were rapidly evolving. In contrast, codon switching at non-essential Ser residues in the kringle domains of coagulation and fibrinolytic proteases are quite common and could be identified in half of the kringles analyzed.

CONCLUSION

Codon interconversion of essential Ser residues of coagulation and fibrinolytic proteases only occurred in genes rapidly evolving and which-at least in some cases-evolved following genome duplication. Interconversion is common at non-essential Ser residues as found in kringle domains.

Evolution and virulence of porcine epidemic diarrhea virus following in vitro and in vivo propagation

Practice of inoculating porcine epidemic diarrhea virus (PEDV) in piglets generating feedback material might influence the genetic evolution and attenuation of PEDV. The study was conducted to evaluate evolutionary rate and attenuation following serial in vitro and in vivo propagation. In the study, PED-JPFP0-PJ, Passage 0 (P0), was isolated from infected pigs and serially passaged in Vero cells for 5 consecutive times, P1-P5. P0, P2 and P5 were then subjected to orally inoculate 3-day-old piglets. At 24 h post inoculation, intestines of each passage (F1), were collected, and subsequently sub-passaged in piglets for 2 additional passages (F2-F3). Virus titration, PEDV genomic copies number, VH:CD ratios, and immunohistochemistry were evaluated. S and ORF3 genes were characterized. The results of the study demonstrated that virus titer and virulence were negatively correlated with increased passages, both in vitro and in vivo. Increased substitution rate was observed in higher passages. The evolutionary rate of S gene was higher than that of ORF3. Seven aa changes at positions 223, 291, 317, 607, 694, 1114 and 1199, with reduced N-linked glycan were observed in P5F3. In conclusion, serial passage of PEDV, both in vitro and in vivo, influence the genetic development and the attenuation of PEDV.

Functional fixedness in chimpanzees

Differences in the tool use of non-human primates and humans are subject of ongoing debate. In humans, representations of object functions underpin efficient tool use. Yet, representations of object functions can lead to functional fixedness, which describes the fixation on a familiar tool function leading to less efficient problem solving when the problem requires using the tool for a new function. In the current study, we examined whether chimpanzees exhibit functional fixedness. After solving a problem with a tool, chimpanzees were less efficient in solving another problem which required using the same tool with a different function compared to a control group. This fixation effect was still apparent after a period of nine months and when chimpanzees had learned about the function of a tool by observation of a conspecific. These results suggest that functional fixedness in our closest living relatives likely exists and cast doubt on the notion that stable function representations are uniquely human.

The diverse liver viromes of Australian geckos and skinks are dominated by hepaciviruses and picornaviruses and reflect host taxonomy and habitat

Lizards have diverse ecologies and evolutionary histories, and represent a promising group to explore how hosts shape virome structure and virus evolution. Yet, little is known about the viromes of these animals. In Australia, squamates (lizards and snakes) comprise the most diverse order of vertebrates, and Australia hosts the highest diversity of lizards globally, with the greatest breadth of habitat use. We used meta-transcriptomic sequencing to determine the virome of nine co-distributed, tropical lizard species from three taxonomic families in Australia and analyzed these data to identify host traits associated with viral abundance and diversity. We show that lizards carry a large diversity of viruses, identifying more than thirty novel, highly divergent vertebrate-associated viruses. These viruses were from nine viral families, including several that contain well known pathogens, such as the Flaviviridae, Picornaviridae, Bornaviridae, Iridoviridae, and Rhabdoviridae. Members of the Flaviviridae were particularly abundant across species sampled here, largely belonging to the genus Hepacivirus: fourteen novel hepaciviruses were identified, broadening the known diversity of this group and better defining its evolution by uncovering new reptilian clades. The evolutionary histories of the viruses studied here frequently aligned with the biogeographic and phylogenetic histories of the hosts, indicating that exogenous viruses may help infer host evolutionary history if sampling is strategic and sampling density high enough. Notably, analysis of alpha and beta diversity revealed that virome composition and richness in the animals sampled here was shaped by host taxonomy and habitat. In sum, we identified a diverse range of reptile viruses that broadly contributes to our understanding of virus-host ecology and evolution.

The modern human aryl hydrocarbon receptor is more active when ancestralized by genome editing

The aryl hydrocarbon receptor (AHR) is a transcription factor that has many functions in mammals. Its best known function is that it binds aromatic hydrocarbons and induces the expression of cytochrome P450 genes, which encode enzymes that metabolize aromatic hydrocarbons and other substrates. All present-day humans carry an amino acid substitution at position 381 in the AHR that occurred after the divergence of modern humans from Neandertals and Denisovans. Previous studies that have expressed the ancestral and modern versions of AHR from expression vectors have yielded conflicting results with regard to their activities. Here, we use genome editing to modify the endogenous AHR gene so that it encodes to the ancestral, Neandertal-like AHR protein in human cells. In the absence of exogenous ligands, the expression of AHR target genes is higher in cells expressing the ancestral AHR than in cells expressing the modern AHR, and similar to the expression in chimpanzee cells. Furthermore, the modern human AHR needs higher doses of three ligands than the ancestral AHR to induce the expression of target genes. Thus, the ability of AHR to induce the expression of many of its target genes is reduced in modern humans.

Gene transfer agents: structural and functional properties of domesticated viruses

Horizontal exchange of DNA between bacteria and archaea is prevalent and has major potential implications for genome evolution, plasticity, and population fitness. Several transfer mechanisms have been identified, including gene transfer agents (GTAs). GTAs are intricately regulated domesticated viruses that package host DNA into virus-like capsids and transfer this DNA throughout the bacterial community. Several important advances have recently been made in our understanding of these unusual particles. In this review, we highlight some of these findings, primarily for the model GTA produced by Rhodobacter capsulatus but also for newly identified GTA producers. We provide key insights into these important genetic elements, including the differences between GTAs from their ancestral bacteriophages, their regulation and control, and their elusive evolutionary function.

Transgenerational epigenetic inheritance during plant evolution and breeding

Plants can program and reprogram their genomes to create genetic variation and epigenetic modifications, leading to phenotypic plasticity. Although consequences of genetic changes are comprehensible, the basis for transgenerational inheritance of epigenetic variation is elusive. This review addresses contributions of external (environmental) and internal (genomic) factors to the establishment and maintenance of epigenetic memory during plant evolution, crop domestication, and modern breeding. Dynamic and pervasive changes in DNA methylation and chromatin modifications provide a diverse repertoire of epigenetic variation potentially for transgenerational inheritance. Elucidating and harnessing epigenetic inheritance will help us develop innovative breeding strategies and biotechnological tools to improve crop yield and resilience in the face of environmental challenges. Beyond plants, epigenetic principles are shared across sexually reproducing organisms including humans with relevance to medicine and public health.

Multidimensional primate niche space sheds light on interspecific competition in primate evolution

Characterising how the totality of primate diversity is distributed across the order, and how it evolved, is challenging because diversity in individual traits often show opposing phylogenetic patterns. A species' combination of traits can be conceptualised as its 'niche'. Here, we describe and analyse seven-dimensional niche space, comprising 11 traits, for 191 primate species. Multifaceted diversity is distributed unequally among taxonomic groups. Cercopithecoidea and Hominidae occupy the largest areas of niche space, and are the most diverse families; platyrrhine families occupy small areas, and this space overlaps with strepsirrhines. The evolution of species' locations in niche space is regulated by selection for adaptive optima in trait combinations. Given that niche similarity results in interspecific competition, we quantify two measures of species' niche locations relative to others. We find that omnivores, frugivores, and species tolerating higher temperatures experience stronger interspecific competition. Hominidae occupation of niche space suggests competitive exclusion from niches by Cercopithecoidea over evolutionary time; but living great apes experience the lowest levels of interspecific competition. Callitrichids experience the highest levels of interspecific competition. Our results provide a standardised measure of primate niches that sheds light on the partitioning and evolution of primate diversity, and how this is driven by interspecific competition.

A new and widespread group of fish apicomplexan parasites

Apicomplexans are obligate intracellular parasites that have evolved from a free-living, phototrophic ancestor. They have been reported from marine environmental samples in high numbers,1 with several clades of apicomplexan-related lineages (ARLs) having been described from environmental sequencing data (16S rRNA gene metabarcoding).2 The most notable of these are the corallicolids (previously ARL-V), which possess chlorophyll-biosynthesis genes in their relic chloroplast (apicoplast) and are geographically widespread and abundant symbionts of anthozoans.3 Corallicolids are related to the Eimeriorina, a suborder of apicomplexan coccidians that include other notable members such as Toxoplasma gondii.4Ophioblennius macclurei, the redlip blenny, along with other tropical reef fishes, is known to be infected by Haemogregarina-like and Haemohormidium-like parasites5 supposedly belonging to the Adeleorina; however, phylogenetics shows that these parasites are instead related to the Eimeriorina.6,7 Hybrid genomic sequencing of apicomplexan-infected O. macclurei blood recovered the entire rRNA operon of this apicomplexan parasite along with the complete mitochondrion and apicoplast genomes. Phylogenetic analyses using this new genomic information consistently place these fish-infecting apicomplexans, hereby informally named ichthyocolids, sister to the corallicolids within Coccidia. The apicoplast genome did not contain chlorophyll biosynthesis genes, providing evidence for another independent loss of this pathway within Apicomplexa. Based on the 16S rRNA gene found in the apicoplast, this group corresponds to the previously described ARL-VI. Screening of fish microbiome studies using the plastid 16S rRNA gene shows these parasites to be geographically and taxonomically widespread in fish species across the globe with implications for commercial fisheries and oceanic food webs.

Characterisation of defensins across the marsupial family tree

Defensins are antimicrobial peptides involved in innate immunity, and gene number differs amongst eutherian mammals. Few studies have investigated defensins in marsupials, despite their potential involvement in immunological protection of altricial young. Here we use recently sequenced marsupial genomes and transcriptomes to annotate defensins in nine species across the marsupial family tree. We characterised 35 alpha and 286 beta defensins; gene number differed between species, although Dasyuromorphs had the largest repertoire. Defensins were encoded in three gene clusters within the genome, syntenic to eutherians, and were expressed in the pouch and mammary gland. Marsupial beta defensins were closely related to eutherians, however marsupial alpha defensins were more divergent. We identified marsupial orthologs of human DEFB3 and 6, and several marsupial-specific beta defensin lineages which may have novel functions. Marsupial predicted mature peptides were highly variable in length and sequence composition. We propose candidate peptides for future testing to elucidate the function of marsupial defensins.

The utilization of an unconventional approach to introduce basic bacteriology in a medical school bridge program

Bacteria form an intense portion of reading and learning for students enrolled in microbiology education. As a part of the foundational course outline of bacteriology, bacterial classification is a significant topic of discussion. The purpose of our study was to analyze whether bacterial classification can be taught with a phylogenetic tree approach that might be more engaging and beneficial to student learners of microbiology. This methodology is unique compared to the conventional approach applied in introductory lectures of bacteriology that relies on morphology and Gram-staining to classify bacteria. The participants of this study were students enrolled in a two-semester medical school bridge program that offers a Master's degree in Pre-clinical Sciences. We presented bacterial origin and classification in the light of evolution and used a phylogenetic tree to signify clinically relevant groups of bacteria. Students were also taught the traditional bacterial classification using Gram stains and morphology. Both methods of classification were delivered in a didactic classroom session considering equal time spent and utilizing the same format. An online survey was distributed to the students after the session to collect their feedback. The results from the survey showed that 74% of participants would prefer learning bacterial classification using a combined approach that includes both Gram-staining and morphology as well as the phylogenetic tree. When asked if the study of bacterial classification through an evolutionary tree diagram is a clear and concise way of understanding bacteria, 79% of the students either agreed or strongly agreed with this statement. Interestingly, the alternative phylogenetic tree approach was considered more engaging and regarded as a means to expand the clinical knowledge of bacteria by 78% and 71% of the students, respectively. Overall, our study strongly supports the use of tree-based classification as an additional method to improve the learning of medically important groups of bacteria at varying levels of education.

An evolutionary view of CC-type glutaredoxins functions in land plants development and adaptation to environment

Land plants have to face an oxidizing, heterogeneous and fast changing environment. Redox-dependent post-translational modifications emerge as a critical component of plant responses to stresses. Among thiols oxidoreductases superfamily, class III CC-type glutaredoxins (called ROXYs) are land plant specific, and their evolutionary history is highly dynamic. Angiosperms encode many isoforms classified into five subgroups (Aα, Aβ, Bα, Bβ, Bγ) that probably evolved from five common ancestral ROXYs, with higher evolutionary dynamics in Bγ compared to other subgroups. ROXYs can modulate the transcriptional activity of TGAs transcription factors target genes, although their biochemical function is still debated. ROXYs participate in the control of proper plant development and reproduction, and are mainly negative regulators of plant responses to biotic and abiotic stresses. This suggests that most ROXYs could play essential and conserved functions in resetting redox-dependent changes in transcriptional activity upon stress signaling to ensure the responsiveness of the system and/or avoid exaggerated responses that could lead to major defects in plant growth and reproduction. Bγ members in Arabidopsis acquired important functions in responses to nitrogen availability and endogenous status, but the rapid and independent evolution of this subclass could suggest that this function results from neofunctionalization, specifically observed in core Eudicots.

Genome evolution in intracellular parasites: Microsporidia and Apicomplexa

Microsporidia and Apicomplexa are eukaryotic, single-celled, intracellular parasites with huge public health and economic importance. Typically, these parasites are studied separately, emphasizing their uniqueness and diversity. In this review, we explore the huge amount of genomic data that has recently become available for the two groups. We compare and contrast their genome evolution and discuss how their transitions to intracellular life may have shaped it. In particular, we explore genome reduction and compaction, genome expansion and ploidy, gene shuffling and rearrangements, and the evolution of centromeres and telomeres.

Diversity, evolution, and emergence of fish viruses

The production of aquatic animals has more than doubled over the last 50 years and is anticipated to continually increase. While fish are recognized as a valuable and sustainable source of nutrition, particularly in the context of human population growth and climate change, the rapid expansion of aquaculture coincides with the emergence of highly pathogenic viruses that often spread globally through aquacultural practices. Here, we provide an overview of the fish virome and its relevance for disease emergence, with a focus on the insights gained through metagenomic sequencing, noting potential areas for future study. In particular, we describe the diversity and evolution of fish viruses, for which the majority have no known disease associations, and demonstrate how viruses emerge in fish populations, most notably at an expanding domestic-wild interface. We also show how wild fish are a powerful and tractable model system to study virus ecology and evolution more broadly and can be used to identify the major factors that shape vertebrate viromes. Central to this is a process of virus-host co-divergence that proceeds over many millions of years, combined with ongoing cross-species virus transmission.

SMC-mediated dosage compensation in C. elegans evolved in the presence of an ancestral nematode mechanism

Mechanisms of X chromosome dosage compensation have been studied extensively in a few model species representing clades of shared sex chromosome ancestry. However, the diversity within each clade as a function of sex chromosome evolution is largely unknown. Here, we anchor ourselves to the nematode Caenorhabditis elegans, for which a well-studied mechanism of dosage compensation occurs through a specialized structural maintenance of chromosomes (SMC) complex, and explore the diversity of dosage compensation in the surrounding phylogeny of nematodes. Through phylogenetic analysis of the C. elegans dosage compensation complex and a survey of its epigenetic signatures, including X-specific topologically associating domains (TADs) and X-enrichment of H4K20me1, we found that the condensin-mediated mechanism evolved recently in the lineage leading to Caenorhabditis through an SMC-4 duplication. Intriguingly, an independent duplication of SMC-4 and the presence of X-specific TADs in Pristionchus pacificus suggest that condensin-mediated dosage compensation arose more than once. mRNA-seq analyses of gene expression in several nematode species indicate that dosage compensation itself is ancestral, as expected from the ancient XO sex determination system. Indicative of the ancestral mechanism, H4K20me1 is enriched on the X chromosomes in Oscheius tipulae, which does not contain X-specific TADs or SMC-4 paralogs. Together, our results indicate that the dosage compensation system in C. elegans is surprisingly new, and condensin may have been co-opted repeatedly in nematodes, suggesting that the process of evolving a chromosome-wide gene regulatory mechanism for dosage compensation is constrained.

Significance statement

X chromosome dosage compensation mechanisms evolved in response to Y chromosome degeneration during sex chromosome evolution. However, establishment of dosage compensation is not an endpoint. As sex chromosomes change, dosage compensation strategies may have also changed. In this study, we performed phylogenetic and epigenomic analyses surrounding Caenorhabditis elegans and found that the condensin-mediated dosage compensation mechanism in C. elegans is surprisingly new, and has evolved in the presence of an ancestral mechanism. Intriguingly, condensin-based dosage compensation may have evolved more than once in the nematode lineage, the other time in Pristionchus. Together, our work highlights a previously unappreciated diversity of dosage compensation mechanisms within a clade, and suggests constraints in evolving new mechanisms in the presence of an existing one.

Illuminating the function of the orphan transporter, SLC22A10, in humans and other primates

SLC22A10 is an orphan transporter with unknown substrates and function. The goal of this study is to elucidate its substrate specificity and functional characteristics. In contrast to orthologs from great apes, human SLC22A10, tagged with green fluorescent protein, is not expressed on the plasma membrane. Cells expressing great ape SLC22A10 orthologs exhibit significant accumulation of estradiol-17β-glucuronide, unlike those expressing human SLC22A10. Sequence alignments reveal a proline at position 220 in humans, which is a leucine in great apes. Replacing proline with leucine in SLC22A10-P220L restores plasma membrane localization and uptake function. Neanderthal and Denisovan genomes show proline at position 220, akin to modern humans, indicating functional loss during hominin evolution. Human SLC22A10 is a unitary pseudogene due to a fixed missense mutation, P220, while in great apes, its orthologs transport sex steroid conjugates. Characterizing SLC22A10 across species sheds light on its biological role, influencing organism development and steroid homeostasis.

Ultrasound description of the coelomic cavity of the axolotl (Ambystoma mexicanum) in a clinically healthy population: a pilot study

Axolotls (Ambystoma mexicanum) are extensively studied for their relevance in human medical research. Despite being critically endangered in the wild, they have gained popularity as household pets. Although they have been kept in captivity for over a century, detailed descriptions of their coelomic organ anatomy remain limited. Also, this species exhibits significant variations compared to other amphibians. Ultrasound is a non-invasive and painless medical imaging technique, ideally suited for investigating internal organs or structures. This study focused on describing the ultrasound appearance of the axolotl coelomic cavity. It details the identification, localization and parenchymal description of major organs in 28 neotenic axolotls using ultrasound frequencies ranging from 7 to 15 MHz. The accuracy of the results was validated by comparing ultrasound findings with necropsy results from one male and one female axolotl. The heart, lung surface, liver and reproductive tracts were visualized. Measurements, along with confidence intervals, were calculated for the spleen, kidneys, testicles, gastric wall, gallbladder, and pylorus. Occasional detection of hyperechoic millimetric particles in the gallbladder or ascites was noted. However, visualization of the pancreas and bladder was not possible. This research outcomes involve the development of a comprehensive atlas comprising images obtained throughout the study. Additionally, the experiment established a reproducible and readily accessible protocol for conducting anatomy-morphological assessments in axolotl medicine. This protocol stands as a crucial preliminary stage before advancing to lesion identification.

Satiety controls behavior in Hydra through an interplay of pre-enteric and central nervous system-like neuron populations

Hunger and satiety can have an influence on decision-making, sensory processing, and motor behavior by altering the internal state of the brain. This process necessitates the integration of peripheral sensory stimuli into the central nervous system. Here, we show how animals without a central nervous system such as the cnidarian Hydra measure and integrate satiety into neuronal circuits and which specific neuronal populations are involved. We demonstrate that this simple nervous system, previously referred to as diffuse, has an endodermal subpopulation (N4) similar to the enteric nervous system (feeding-associated behavior) and an ectodermal population (N3) that performs central nervous system-like functions (physiology/motor). This view of a supposedly simple nervous system could open an important window into the origin of more complex nervous systems.

DNA methylation reprogramming in teleosts

Early embryonic development is crucially important but also remarkably diverse among animal taxa. Axis formation and cell lineage specification occur due to both spatial and temporal control of gene expression. This complex system involves various signaling pathways and developmental genes such as transcription factors as well as other molecular interactants that maintain cellular states, including several types of epigenetic marks. 5mC DNA methylation, the chemical modification of cytosines in eukaryotes, represents one such mark. By influencing the compaction of chromatin (a high-order DNA structure), DNA methylation can either repress or induce transcriptional activity. Mammals exhibit a reprogramming of DNA methylation from the parental genomes in the zygote following fertilization, and later in primordial germ cells (PGCs). Whether these periods of methylation reprogramming are evolutionarily conserved, or an innovation in mammals, is an emerging question. Looking into these processes in other vertebrate lineages is thus important, and teleost fish, with their extensive species richness, phenotypic diversity, and multiple rounds of whole genome duplication, provide the perfect research playground for answering such a question. This review aims to present a concise state of the art of DNA methylation reprogramming in early development in fish by summarizing findings from different research groups investigating methylation reprogramming patterns in teleosts, while keeping in mind the ramifications of the methodology used, then comparing those patterns to reprogramming patterns in mammals.

Arginines of the CGN codon family are Achilles' heels of cancer genes

Recent studies have revealed that arginine is the most favorable target of amino acid alteration in most cancer types and it has been suggested that the high preference for arginine mutations reflects the critical roles of this amino acid in the function of proteins. High rates of mutations of arginine residues in cancer, however, might also be due to increased mutability of arginine codons of the CGN family as the CpG dinucleotides of these codons may be methylated. In the present work we have analyzed spectra of single base substitutions of cancer genes (oncogenes, tumor suppressor genes) and passenger genes in cancer tissues to assess the contributions of CpG hypermutability and selection to arginine mutations. Our studies have shown that arginines encoded by the CGN codon family display higher rates of mutation in both cancer genes and passenger genes than arginine codons AGA and AGG that are devoid of CpG dinucleotide, suggesting that the predominance of arginine mutations in cancer is primarily due to CpG hypermutability, rather than selection for arginine replacement. Nevertheless, our results also suggest that CGN codons for arginines may serve as Achilles' heels of cancer genes. CpG hypermutability of key arginines of proto-oncogenes, leading to high rates of recurrence of driver mutations, contributes significantly to carcinogenesis. Similarly, our results indicate that hypermutability of the CpG dinucleotide of CGA codons (converting them to TGA stop codons) contributes significantly to recurrent truncation and inactivation of tumor suppressor genes.

Using structure prediction of negative sense RNA virus nucleoproteins to assess evolutionary relationships

Negative sense RNA viruses (NSV) include some of the most detrimental human pathogens, including the influenza, Ebola and measles viruses. NSV genomes consist of one or multiple single-stranded RNA molecules that are encapsidated into one or more ribonucleoprotein (RNP) complexes. These RNPs consist of viral RNA, a viral RNA polymerase, and many copies of the viral nucleoprotein (NP). Current evolutionary relationships within the NSV phylum are based on alignment of conserved RNA-directed RNA polymerase (RdRp) domain amino acid sequences. However, the RdRp domain-based phylogeny does not address whether NP, the other core protein in the NSV genome, evolved along the same trajectory or whether several RdRp-NP pairs evolved through convergent evolution in the segmented and non-segmented NSV genomes architectures. Addressing how NP and the RdRp domain evolved may help us better understand NSV diversity. Since NP sequences are too short to infer robust phylogenetic relationships, we here used experimentally-obtained and AlphaFold 2.0-predicted NP structures to probe whether evolutionary relationships can be estimated using NSV NP sequences. Following flexible structure alignments of modeled structures, we find that the structural homology of the NSV NPs reveals phylogenetic clusters that are consistent with RdRp-based clustering. In addition, we were able to assign viruses for which RdRp sequences are currently missing to phylogenetic clusters based on the available NP sequence. Both our RdRp-based and NP-based relationships deviate from the current NSV classification of the segmented Naedrevirales, which cluster with the other segmented NSVs in our analysis. Overall, our results suggest that the NSV RdRp and NP genes largely evolved along similar trajectories and that even short pieces of genetic, protein-coding information can be used to infer evolutionary relationships, potentially making metagenomic analyses more valuable.

Genome-Wide Identification, Expression Analysis and Functional Study of DELLA Genes in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

BACKGROUND

DELLA protein is a crucial factor which played pivotal roles in regulating numerous intriguing biological processes in plant development and abiotic stress responses. However, little is known about the function and information of DELLA protein in Chinese cabbage.

METHODS

Using 5 DELLA gene sequences in Arabidopsis Thaliana as probes, 5 DELLA genes in Chinese cabbage were identified by Blast search in Chinese cabbage database (Brassica database (BRAD)). The National Center for Biotechnology Information (NCBI), ExPaSy, SWISS-MODEL, DNAMAN, MEGA 11, PlantCARE were used to identify and analyze the DELLA gene family of Chinese cabbage. Gene expression was analyzed by quantitative real-time polymerase chain reaction (qRT-PCR). The function of BraA10gRGL3 was verified by overexpression and phenotypic analysis of BraA10gRGL3 and yeast hybrid.

RESULTS

In this study, 5 BraDELLAs homologous to Arabidopsis thaliana were identified and cloned based on the Brassica database, namely, BraA02gRGL1, BraA05gRGL2, BraA10gRGL3, BraA06gRGA and BraA09gRGA. All BraDELLAs contain the DELLA, TVHYNP, and GRAS conserved domains. Cis-element analysis revealed that the promoter regions of these 5 DELLA genes all contain light-responsive elements, TCT motif, I-box, G-box, and box 4, which are associated with GA signaling. Transcriptome analysis results proved that the expression of BraA02gRGL1, BraA05gRGL2, and BraA10gRGL3 in Y2 at different growth stages were lower than them in Y7, which is consistent with the phenotype that Y7 exhibited stronger stress tolerance than Y2. It is worth emphasizing that even through the overexpression of BraA10gRGL3-Y7 in Arabidopsis resulted in smaller leaf size and lower fresh weight compared to the wild type (WT) Arabidopsis: Columbia, a stronger response to abiotic stresses was observed in BraA10gRGL3-Y7. It indicated that BraA10gRGL3-Y7 can improve the stress resistance of plants by inhibiting their growth. Moreover, the yeast two-hybrid experiment confirmed that BraA10gRGL3-Y7 can interacted with BraA05gGID1a-Y7, BraA04gGID1b1, BraA09gGID1b3-Y2, and BraA06gGID1c, whereas BraA10gRGL3-Y2 cannot interact with any BraGID1.

CONCLUSIONS

Collectively, BraDELLAs play important role in plant development and response to abiotic stress. The differences in amino acid sequences between BraA10gRGL3-Y2 and BraA10gRGL3-Y7 may result in variations in their protein binding sites, thus affecting their interaction with the BraGID1 family proteins. This systematic analysis lays the foundation for further study of the functional characteristics of DELLA genes of Chinese cabbage.

Evolution of endosymbiosis-mediated nuclear calcium signaling in land plants

The ability of fungi to establish mycorrhizal associations with plants and enhance the acquisition of mineral nutrients stands out as a key feature of terrestrial life. Evidence indicates that arbuscular mycorrhizal (AM) association is a trait present in the common ancestor of land plants,1,2,3,4 suggesting that AM symbiosis was an important adaptation for plants in terrestrial environments.5 The activation of nuclear calcium signaling in roots is essential for AM within flowering plants.6 Given that the earliest land plants lacked roots, whether nuclear calcium signals are required for AM in non-flowering plants is unknown. To address this question, we explored the functional conservation of symbiont-induced nuclear calcium signals between the liverwort Marchantia paleacea and the legume Medicago truncatula. In M. paleacea, AM fungi penetrate the rhizoids and form arbuscules in the thalli.7 Here, we demonstrate that AM germinating spore exudate (GSE) activates nuclear calcium signals in the rhizoids of M. paleacea and that this activation is dependent on the nuclear-localized ion channel DOES NOT MAKE INFECTIONS 1 (MpaDMI1). However, unlike flowering plants, MpaDMI1-mediated calcium signaling is only required for the thalli colonization but not for the AM penetration within rhizoids. We further demonstrate that the mechanism of regulation of DMI1 has diverged between M. paleacea and M. truncatula, including a key amino acid residue essential to sustain DMI1 in an inactive state. Our study reveals functional evolution of nuclear calcium signaling between liverworts and flowering plants and opens new avenues of research into the mechanism of endosymbiosis signaling.