Evolutionary and immune-activating character analyses of NLR genes in algae suggest the ancient origin of plant intracellular immune receptors.

Nucleotide-binding leucine-rich repeat (NLR) proteins are crucial intracellular immune receptors in plants, responsible for detecting invading pathogens and initiating defense responses. While previous studies on the evolution and function of NLR genes were mainly limited to land plants, the evolutionary trajectory and immune-activating character of NLR genes in algae remain less explored. In this study, genome-wide NLR gene analysis was conducted on 44 chlorophyte species across seven classes and seven charophyte species across five classes. A few but variable number of NLR genes, ranging from one to 20, were identified in five chlorophytes and three charophytes, whereas no NLR gene was identified from the remaining algal genomes. Compared with land plants, algal genomes possess fewer or usually no NLR genes, implying that the expansion of NLR genes in land plants can be attributed to their adaptation to the more complex terrestrial pathogen environments. Through phylogenetic analysis, domain composition analysis, and conserved motifs profiling of the NBS domain, we detected shared and lineage-specific features between NLR genes in algae and land plants, supporting the common origin and continuous evolution of green plant NLR genes. Immune-activation assays revealed that both TNL and RNL proteins from green algae can elicit hypersensitive responses in Nicotiana benthamiana, indicating the molecular basis for immune activation has emerged in the early evolutionary stage of different types of NLR proteins. In summary, the results from this study suggest that NLR proteins may have taken a role as intracellular immune receptors in the common ancestor of green plants.

Genome-Wide Identification and Evolutionary Analysis of Receptor-like Kinase Family Genes Provides Insights into Anthracnose Resistance of Dioscorea alata.

Dioscorea alata, commonly known as "greater yam", is a vital crop in tropical and subtropical regions of the world, yet it faces significant threats from anthracnose disease, mainly caused by Colletotrichum gloeosporioides. However, exploring disease resistance genes in this species has been challenging due to the difficulty of genetic mapping resulting from the loss of the flowering trait in many varieties. The receptor-like kinase (RLK) gene family represents essential immune receptors in plants. In this study, genomic analysis revealed 467 RLK genes in D. alata. The identified RLKs were distributed unevenly across chromosomes, likely due to tandem duplication events. However, a considerable number of ancient whole-genome or segmental duplications dating back over 100 million years contributed to the diversity of RLK genes. Phylogenetic analysis unveiled at least 356 ancient RLK lineages in the common ancestor of Dioscoreaceae, which differentially inherited and expanded to form the current RLK profiles of D. alata and its relatives. The analysis of cis-regulatory elements indicated the involvement of RLK genes in diverse stress responses. Transcriptome analysis identified RLKs that were up-regulated in response to C. gloeosporioides infection, suggesting their potential role in resisting anthracnose disease. These findings provide novel insights into the evolution of RLK genes in D. alata and their potential contribution to disease resistance.

Comprehensive regulatory networks for tomato organ development based on the genome and RNAome of MicroTom tomato.

MicroTom has a short growth cycle and high transformation efficiency, and is a prospective model plant for studying organ development, metabolism, and plant-microbe interactions. Here, with a newly assembled reference genome for this tomato cultivar and abundant RNA-seq data derived from tissues of different organs/developmental stages/treatments, we constructed multiple gene co-expression networks, which will provide valuable clues for the identification of important genes involved in diverse regulatory pathways during plant growth, e.g. arbuscular mycorrhizal symbiosis and fruit development. Additionally, non-coding RNAs, including miRNAs, lncRNAs, and circRNAs were also identified, together with their potential targets. Interacting networks between different types of non-coding RNAs (miRNA-lncRNA), and non-coding RNAs and genes (miRNA-mRNA and lncRNA-mRNA) were constructed as well. Our results and data will provide valuable information for the study of organ differentiation and development of this important fruit. Lastly, we established a database (http://eplant.njau.edu.cn/microTomBase/) with genomic and transcriptomic data, as well as details of gene co-expression and interacting networks on MicroTom, and this database should be of great value to those who want to adopt MicroTom as a model plant for research.

The evolution of plant NLR immune receptors and downstream signal components.

Along with the emergence of green plants on this planet one billion years ago, the nucleotide binding site leucine-rich repeat (NLR) gene family originated and diverged into at least three subclasses. Two of them, with either characterized N-terminal toll/interleukin-1 receptor (TIR) or coiled-coil (CC) domain, serve as major types of immune receptor of effector-triggered immunity (ETI) in plants, whereas the one having a N-terminal Resistance to powdery mildew8 (RPW8) domain, functions as signal transfer component to them. In this review, we briefly summarized the history of identification of diverse NLR subclasses across Viridiplantae lineages during the establishment of NLR category, and highlighted recent advances on the evolution of NLR genes and several key downstream signal components under the background of ecological adaption.

The Saururus chinensis genome provides insights into the evolution of pollination strategies and herbaceousness in magnoliids.

Saururus chinensis, an herbaceous magnoliid without perianth, represents a clade of early-diverging angiosperms that have gone through woodiness-herbaceousness transition and pollination obstacles: the characteristic white leaves underneath inflorescence during flowering time are considered a substitute for perianth to attract insect pollinators. Here, using the newly sequenced S. chinensis genome, we revisited the phylogenetic position of magnoliids within mesangiosperms, and recovered a sister relationship for magnoliids and Chloranthales. By considering differentially expressed genes, we identified candidate genes that are involved in the morphogenesis of the white leaves in S. chinensis. Among those genes, we verified - in a transgenic experiment with Arabidopsis - that increasing the expression of the "pseudo-etiolation in light" gene (ScPEL) can inhibit the biosynthesis of chlorophyll. ScPEL is thus likely responsible for the switches between green and white leaves, suggesting that changes in gene expression may underlie the evolution of pollination strategies. Despite being an herbaceous plant, S. chinensis still has vascular cambium and maintains the potential for secondary growth as a woody plant, because the necessary machinery, i.e., the entire gene set involved in lignin biosynthesis, is well preserved. However, similar expression levels of two key genes (CCR and CAD) between the stem and other tissues in the lignin biosynthesis pathway are possibly associated with the herbaceous nature of S. chinensis. In conclusion, the S. chinensis genome provides valuable insights into the adaptive evolution of pollination in Saururaceae and reveals a possible mechanism for the evolution of herbaceousness in magnoliids.

The RNAome landscape of tomato during arbuscular mycorrhizal symbiosis reveals an evolving RNA layer symbiotic regulatory network.

Arbuscular mycorrhizal symbiosis (AMS) is an ancient plant-fungus relationship that is widely distributed in terrestrial plants. The formation of symbiotic structures and bidirectional nutrient exchange requires the regulation of numerous genes. However, the landscape of RNAome during plant AMS involving different types of regulatory RNA is poorly understood. In this study, a combinatorial strategy utilizing multiple sequencing approaches was used to decipher the landscape of RNAome in tomato, an emerging AMS model. The annotation of the tomato genome was improved by a multiple-platform sequencing strategy. A total of 3,174 protein-coding genes were upregulated during AMS, 42% of which were alternatively spliced. Comparative-transcriptome analysis revealed that genes from 24 orthogroups were consistently induced by AMS in eight phylogenetically distant angiosperms. Seven additional orthogroups were specifically induced by AMS in all surveyed dicot AMS host plants. However, these orthogroups were absent or not induced in monocots and/or non-AMS hosts, suggesting a continuously evolving AMS-responsive network in addition to a conserved core regulatory module. Additionally, we detected 587 lncRNAs, ten miRNAs, and 146 circRNAs that responded to AMS, which were incorporated to establish a tomato AMS-responsive, competing RNA-responsive endogenous RNA (ceRNA) network. Finally, a tomato symbiotic transcriptome database (TSTD, https://efg.nju.edu.cn/TSTD) was constructed to serve as a resource for deep deciphering of the AMS regulatory network. These results help elucidate the reconfiguration of the tomato RNAome during AMS and suggest a sophisticated and evolving RNA layer responsive network during AMS processes.

The Cycas genome and the early evolution of seed plants.

Cycads represent one of the most ancient lineages of living seed plants. Identifying genomic features uniquely shared by cycads and other extant seed plants, but not non-seed-producing plants, may shed light on the origin of key innovations, as well as the early diversification of seed plants. Here, we report the 10.5-Gb reference genome of Cycas panzhihuaensis, complemented by the transcriptomes of 339 cycad species. Nuclear and plastid phylogenomic analyses strongly suggest that cycads and Ginkgo form a clade sister to all other living gymnosperms, in contrast to mitochondrial data, which place cycads alone in this position. We found evidence for an ancient whole-genome duplication in the common ancestor of extant gymnosperms. The Cycas genome contains four homologues of the fitD gene family that were likely acquired via horizontal gene transfer from fungi, and these genes confer herbivore resistance in cycads. The male-specific region of the Y chromosome of C. panzhihuaensis contains a MADS-box transcription factor expressed exclusively in male cones that is similar to a system reported in Ginkgo, suggesting that a sex determination mechanism controlled by MADS-box genes may have originated in the common ancestor of cycads and Ginkgo. The C. panzhihuaensis genome provides an important new resource of broad utility for biologists.

Comparative Analysis of HSF Genes From Secale cereale and its Triticeae Relatives Reveal Ancient and Recent Gene Expansions.

Plants have evolved sophisticated systems to cope with the environmental stresses, with the heat shock factor (HSF) family proteins composing an integral part of the transcriptional regulation system. Understanding the evolutionary history and functional diversity of HSFs will facilitate improving tolerance of crops to adverse environmental conditions. In this study, genome-wide analysis of Secale cereale identified 31 HSF genes. The total number of HSF genes in S. cereale is larger than that in barley and the three subgenomes of wheat, suggesting it is a valuable resource for mining functional HSFs. Chromosome analysis revealed an uneven distribution of HSF genes among the 7 S. cereale chromosomes, with no HSF gene was detected on chromosome 4. Further interspecies synteny analysis revealed that chromosome reorganization during species-speciation may lead to the escape of HSF genes from the S. cereale chromosome 4. Phylogenetic analysis revealed that S. cereale experienced more HSF gene duplications than barley and the three wheat subgenomes. Expression analysis demonstrated that S. cereale HSF genes showed diverse expression patterns across plant developmental stages and upon drought and freezing treatment, suggesting functional diversity of the gene family. Notably, we detected distinct expression patterns for a recently duplicated HSF gene pair, indicating functional divergence may have occurred between the two genes. The study presents the genome organization, evolutionary features and expression patterns of the S. cereale HSF genes. These results provide new insights into the evolution of HSF genes in Triticeae and may serve as a resource for Triticeae molecular breeding.

An angiosperm NLR Atlas reveals that NLR gene reduction is associated with ecological specialization and signal transduction component deletion.

Nucleotide-binding leucine-rich-repeat (NLR) genes comprise the largest family of plant disease-resistance genes. Angiosperm NLR genes are phylogenetically divided into the TNL, CNL, and RNL subclasses. NLR copy numbers and subclass composition vary tremendously across angiosperm genomes. However, the evolutionary associations between genomic NLR content and ecological adaptation, or between NLR content and signal transduction components, are poorly characterized because of limited genome availability. In this study, we established an angiosperm NLR atlas (ANNA, https://biobigdata.nju.edu.cn/ANNA/) that includes NLR genes from over 300 angiosperm genomes. Using ANNA, we revealed that NLR copy numbers differ up to 66-fold among closely related species owing to rapid gene loss and gain. Interestingly, NLR contraction was associated with adaptations to aquatic, parasitic, and carnivorous lifestyles. The convergent NLR reduction in aquatic plants resembles the lack of NLR expansion during the long-term evolution of green algae before the colonization of land. A co-evolutionary pattern between NLR subclasses and plant immune pathway components was also identified, suggesting that immune pathway deficiencies may drive TNL loss. Finally, we identified a conserved TNL lineage that may function independently of the EDS1-SAG101-NRG1 module. Collectively, these findings provide new insights into the evolution of NLR genes in the context of ecological adaptation and genome content variation.

A chromosome-level genome assembly of rugged rose (Rosa rugosa) provides insights into its evolution, ecology, and floral characteristics.

Rosa rugosa, commonly known as rugged rose, is a perennial ornamental shrub. It produces beautiful flowers with a mild fragrance and colorful seed pods. Unlike many other cultivated roses, R. rugosa adapts to a wide range of habitat types and harsh environmental conditions such as salinity, alkaline, shade, drought, high humidity, and frigid temperatures. Here, we produced and analyzed a high-quality genome sequence for R. rugosa to understand its ecology, floral characteristics and evolution. PacBio HiFi reads were initially used to construct the draft genome of R. rugosa, and then Hi-C sequencing was applied to assemble the contigs into 7 chromosomes. We obtained a 382.6 Mb genome encoding 39,704 protein-coding genes. The genome of R. rugosa appears to be conserved with no additional whole-genome duplication after the gamma whole-genome triplication (WGT), which occurred ~100 million years ago in the ancestor of core eudicots. Based on a comparative analysis of the high-quality genome assembly of R. rugosa and other high-quality Rosaceae genomes, we found a unique large inverted segment in the Chinese rose R. chinensis and a retroposition in strawberry caused by post-WGT events. We also found that floral development- and stress response signaling-related gene modules were retained after the WGT. Two MADS-box genes involved in floral development and the stress-related transcription factors DREB2A-INTERACTING PROTEIN 2 (DRIP2) and PEPTIDE TRANSPORTER 3 (PTR3) were found to be positively selected in evolution, which may have contributed to the unique ability of this plant to adapt to harsh environments. In summary, the high-quality genome sequence of R. rugosa provides a map for genetic studies and molecular breeding of this plant and enables comparative genomic studies of Rosa in the near future.

Genome-Wide Analysis of NLR Disease Resistance Genes in an Updated Reference Genome of Barley.

Barley is one of the top 10 crop plants in the world. During its whole lifespan, barley is frequently infected by various pathogens. In this study, we performed genome-wide analysis of the largest group of plant disease resistance (R) genes, the nucleotide binding site-leucine-rich repeat receptor (NLR) gene, in an updated barley genome. A total of 468 NLR genes were identified from the improved barley genome, including one RNL subclass and 467 CNL subclass genes. Proteins of 43 barley CNL genes were shown to contain 25 different integrated domains, including WRKY and BED. The NLR gene number identified in this study is much larger than previously reported results in earlier versions of barley genomes, and only slightly fewer than that in the diploid wheat Triticum urartu. Barley Chromosome 7 contains the largest number of 112 NLR genes, which equals to seven times of the number of NLR genes on Chromosome 4. The majority of NLR genes (68%) are located in multigene clusters. Phylogenetic analysis revealed that at least 18 ancestral CNL lineages were presented in the common ancestor of barley, T. urartu and Arabidopsis thaliana. Among them fifteen lineages expanded to 533 sub-lineages prior to the divergence of barley and T. urartu. The barley genome inherited 356 of these sub-lineages and duplicated to the 467 CNL genes detected in this study. Overall, our study provides an updated profile of barley NLR genes, which should serve as a fundamental resource for functional gene mining and molecular breeding of barley.

Inactivation of the htpsA gene affects capsule development and pathogenicity of Streptococcus suis.

STREPTOCOCCUS SUIS: serotype 2 (S. suis 2) is an important swine pathogen and also an emerging zoonotic agent. HtpsA has been reported as an immunogenic cell surface protein on the bacterium. In the present study, we constructed an isogenic mutant strain of htpsA, namely ΔhtpsA, to study its role in the development and virulence of S. suis 2. Our results showed that the mutant strain lost its typical encapsulated structure with decreased concentrations of sialic acid. Furthermore, the survival rate in whole blood, the anti-phagocytosis by RAW264.7 murine macrophage, and the adherence ability to HEp-2 cells were all significantly affected in the ΔhtpsA. In addition, the deletion of htpsA sharply attenuated the virulence of S. suis 2 in an infection model of mouse. RNA-seq analysis revealed that 126 genes were differentially expressed between the ΔhtpsA and the wild-type strains, including 28 upregulated and 98 downregulated genes. Among the downregulated genes, many were involved in carbohydrate metabolism and synthesis of virulence-associated factors. Taken together, htpsA was demonstrated to play a role in the morphological development and pathogenesis of the highly virulent S. suis 2 05ZYH33 strain.

Maternal Inheritance of U's Triangle and Evolutionary Process of Brassica Mitochondrial Genomes.

The sequences and genomic structures of plant mitochondrial (mt) genomes provide unique material for phylogenetic studies. The nature of uniparental inheritance renders an advantage when utilizing mt genomes for determining the parental sources of hybridized taxa. In this study, a concatenated matrix of mt genes was used to infer the phylogenetic relationships of six cultivated Brassica taxa and explore the maternal origins of three allotetraploids. The well-resolved sister relationships between two pairs of diploid and allotetraploid taxa suggest that Brassica carinata (car) possessed a maternal origin from Brassica nigra, while Brassica juncea (jun) was maternally derived from Brassica rapa (cam). Another allotetraploid taxon, Brassica napus (cv. Wester) may have been maternally derived from the common ancestor of B. rapa and Brassica oleracea (ole), and/or have undergone (an) extra hybridization event(s) along its evolutionary history. The characteristics of Brassica mt genomic structures also supported the phylogenetic results. Sinapis arvensis was nested inside the Brassica species, sister to the B. nigra-B. carinata lineage, and possessed an mt genome structure that mostly resembled B. nigra. Collectively, the evidence supported a systematic revision that placed S. arvensis within Brassica. Finally, ancestral mt genomes at each evolutionary node of Brassica were reconstructed, and the detailed and dynamic evolution of Brassica mt genomes was successfully reproduced. The mt genome of B. nigra structurally resembled that of the Brassica ancestor the most, with only one reversion of a block, and the Brassica oleracea underwent the most drastic changes. These findings suggested that repeat-mediated recombinations were largely responsible for the observed structural variations in the evolutionary history of Brassica mt genomes.

Revisiting the Origin of Plant NBS-LRR Genes.

The NBS-LRR genes are functionally responsible for plant resistance to alien pathogens. Here, we show that NBS-LRR genes originated in the common ancestor of the whole green lineage, and have rapidly diverged into three subclasses with different domain combinations (TNL, CNL, and RNL) before the split of green algae.

Genome- Wide Analysis of the Nucleotide Binding Site Leucine-Rich Repeat Genes of Four Orchids Revealed Extremely Low Numbers of Disease Resistance Genes.

Orchids are one of the most diverse flowering plant families, yet possibly maintain the smallest number of the nucleotide-binding site-leucine-rich repeat (NBS-LRR) type plant resistance (R) genes among the angiosperms. In this study, a genome-wide search in four orchid taxa identified 186 NBS-LRR genes. Furthermore, 214 NBS-LRR genes were identified from seven orchid transcriptomes. A phylogenetic analysis recovered 30 ancestral lineages (29 CNL and one RNL), far fewer than other angiosperm families. From the genetics aspect, the relatively low number of ancestral R genes is unlikely to explain the low number of R genes in orchids alone, as historical gene loss and scarce gene duplication has continuously occurred, which also contributes to the low number of R genes. Due to recent sharp expansions, Phalaenopsis equestris and Dendrobium catenatum having 52 and 115 genes, respectively, and exhibited an "early shrinking to recent expanding" evolutionary pattern, while Gastrodia elata and Apostasia shenzhenica both exhibit a "consistently shrinking" evolutionary pattern and have retained only five and 14 NBS-LRR genes, respectively. RNL genes remain in extremely low numbers with only one or two copies per genome. Notably, all of the orchid RNL genes belong to the ADR1 lineage. A separate lineage, NRG1, was entirely absent and was likely lost in the common ancestor of all monocots. All of the TNL genes were absent as well, coincident with the RNL NRG1 lineage, which supports the previously proposed notion that a potential functional association between the TNL and RNL NRG1 genes.

Genome-Wide Analysis Reveals Ancestral Lack of Seventeen Different tRNAs and Clade-Specific Loss of tRNA-CNNs in Archaea.

Transfer RNA (tRNA) is a category of RNAs that specifically decode messenger RNAs (mRNAs) into proteins by recognizing a set of 61 codons commonly adopted by different life domains. The composition and abundance of tRNAs play critical roles in shaping codon usage and pairing bias, which subsequently modulate mRNA translation efficiency and accuracy. Over the past few decades, effort has been concentrated on evaluating the specificity and redundancy of different tRNA families. However, the mechanism and processes underlying tRNA evolution have only rarely been investigated. In this study, by surveying tRNA genes in 167 completely sequenced genomes, we systematically investigated the composition and evolution of tRNAs in Archaea from a phylogenetic perspective. Our data revealed that archaeal genomes are compact in both tRNA types and copy number. Generally, no more than 44 different types of tRNA are present in archaeal genomes to decode the 61 canonical codons, and most of them have only one gene copy per genome. Among them, tRNA-Met was significantly overrepresented, with an average of three copies per genome. In contrast, the tRNA-UAU and 16 tRNAs with A-starting anticodons (tRNA-ANNs) were rarely detected in all archaeal genomes. The conspicuous absence of these tRNAs across the archaeal phylogeny suggests they might have not been evolved in the common ancestor of Archaea, rather than have lost independently from different clades. Furthermore, widespread absence of tRNA-CNNs in the Methanococcales and Methanobacteriales genomes indicates convergent loss of these tRNAs in the two clades. This clade-specific tRNA loss may be attributing to the reductive evolution of their genomes. Our data suggest that the current tRNA profiles in Archaea are contributed not only by the ancestral tRNA composition, but also by differential maintenance and loss of redundant tRNAs.

Identification of oligopeptide-binding protein (OppA) and its role in the virulence of Streptococcus suis serotype 2.

The oligopeptide permease (Opp) cassette, an oligopeptide transport system belongs to the superfamily of ATP-binding cassette (ABC) transporter, is widely distributed in bacteria, including Streptococcus suis (S. suis). It is encoded by the opp operon containing oppA, oppB, oppC, oppD, and oppF. In addition to the uptake of peptide, the oppA gene also plays an important role in virulence of many pathogens. In this study, an oppA homologue from the highly virulent S. suis serotype 2 (S. suis 2) strain 05ZYH33 was identified. Flow cytometry and Western blot confirmed that OppA is a surface immunogenic protein and is expressed during S. suis 2 infection. To explore the role of oppA in S. suis 2 growth and pathogenicity, an isogenic 05ZYH33 mutant of oppA (△oppA) was obtained by homologous recombination. Although the complementary strain was not obtained due to the △oppA strain is not transformable, the current data revealed that deletion of the oppA gene in S. suis 2 has greatly affected its growth and virulence. Our data revealed that the growth rate is significantly slow for the △oppA. Adherence of the △oppA strain to human epithelial cells is greatly reduced comparing to the wild strain. Mouse infection experiment showed that inactivation of oppA greatly attenuated the high pathogenicity of S. suis 2. The observed results suggest that OppA is a surface-exposed protein and plays important roles in the growth and pathogenicity of S. suis 2.

Streptococcus suis DivIVA Protein Is a Substrate of Ser/Thr Kinase STK and Involved in Cell Division Regulation.

Streptococcus suis serotype 2 is an important swine pathogen and an emerging zoonotic agent that causes severe infections. Recent studies have reported a eukaryotic-like Ser/Thr protein kinase (STK) gene and characterized its role in the growth and virulence of different S. suis 2 strains. In the present study, phosphoproteomic analysis was adopted to identify substrates of the STK protein. Seven proteins that were annotated to participate in different cell processes were identified as potential substrates, which suggests the pleiotropic effects of stk on S. suis 2 by targeting multiple pathways. Among them, a protein characterized as cell division initiation protein (DivIVA) was further investigated. In vitro analysis demonstrated that the recombinant STK protein directly phosphorylates threonine at amino acid position 199 (Thr-199) of DivIVA. This effect could be completely abolished by the T199A mutation. To determine the specific role of DivIVA in growth and division, a divIVA mutant was constructed. The ΔdivIVA strain exhibited impaired growth and division, including lower viability, enlarged cell mass, asymmetrical division caused by aberrant septum, and extremely weak pathogenicity in a mouse infection model. Collectively, our results reveal that STK regulates the cell growth and virulence of S. suis 2 by targeting substrates that are involved in different biological pathways. The inactivation of DivIVA leads to severe defects in cell division and strongly attenuates pathogenicity, thereby indicating its potential as a molecular drug target against S. suis.