Search

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.

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 to be 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 being 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.

Sequence discrepancy between genomic DNA and transcribed RNAs, namely RNA editing (RE), was reported to be extensive in metazoan nuclear genomes and plant organellar genomes, but largely unexplored in plant nuclear genomes. By comparing the tomato microTom RNA-seq data to its newly assembled genome, we, for the first time, detected extensive RE events (37,984 sites in 5,873 genes) on in the plant nuclear genome, which is supported by evidence from both experimental validation and proteomic data. Other than the prevalence of RE in the microTom genome, the biased pattern-transition overwhelmingly outnumbering transversion-rivals that in animals, and may imply an ancient and shared mechanism underlying RE across eukaryotes. We observed conserved RE sites that were always edited under all conditions, and could be considered as "fixed mutations" at the RNA level; also we detected conditional RE sites occurring only in specific organ/developmental stage/treatment, possibly suggesting distinct functional roles. Evolutionarily, RE results in sequence diversity at the RNA level and may lead to even functional innovation of proteins, hence may be considered as an additional mechanism to expand gene diversity other than gene duplications and alternative splicing.

The early diversification of angiosperms is thought to have been a rapid process, which may complicate phylogenetic analyses of early angiosperm relationships. Plastid and nuclear phylogenomic studies have raised several conflicting hypotheses regarding overall angiosperm phylogeny, but mitochondrial genomes have been largely ignored as a relevant source of information. Here we sequenced mitochondrial genomes from 18 angiosperms to fill taxon-sampling gaps in Austrobaileyales, magnoliids, Chloranthales, Ceratophyllales, and major lineages of eudicots and monocots. We assembled a data matrix of 38 mitochondrial genes from 107 taxa to assess how well mitochondrial genomic data address current uncertainties in angiosperm relationships. Although we recovered conflicting phylogenies based on different data sets and analytical methods, we also observed congruence regarding deep relationships of several major angiosperm lineages: Chloranthales were always inferred to be the sister group of Ceratophyllales, Austrobaileyales to mesangiosperms, and the unplaced Dilleniales was consistently resolved as the sister to superasterids. Substitutional saturation, GC compositional heterogeneity, and codon-usage bias are possible reasons for the noise/conflict that may impact phylogenetic reconstruction; and angiosperm mitochondrial genes may not be substantially affected by these factors. The third codon positions of the mitochondrial genes appear to contain more parsimony-informative sites than the first and second codon positions, and therefore produced better resolved phylogenetic relationships with generally strong support. The relationships among these major lineages remain incompletely resolved, perhaps as a result of the rapidity of early radiations. Nevertheless, data from mitochondrial genomes provide additional evidence and alternative hypotheses for exploring the early evolution and diversification of the angiosperms. Jiangsu Key Laboratory for the Research and Utilization of Plant Resources (Institute of Botany, Jiangsu Province) [KSPKLB201835]; Jiangsu Key Laboratory for the Research and Utilization of Plant Resources (Chinese Academy of Sciences) [KSPKLB201835]; Shenzhen Key Laboratory of Southern Subtropical Plant Diversity [SLPD-2018-3]; National Natural Science Foundation of ChinaNational Natural Science Foundation of China (NSFC) [31470314, 31471200, 31470448]; Public Welfare Forestry Industry Project of State Forestry Administration, China [201504322]; Special Research Fund of Ghent UniversityGhent University [BOFPDO2018001701]; Strategic Priority Research Program of Chinese Academy of SciencesChinese Academy of Sciences [XDA13020603]; Shenzhen Fairy Lake Botanical Garden Science Foundation [FLSF-2019-02]; Ghent UniversityGhent University; FWOFWO; Flemish Government-department EWI info:eu-repo/semantics/publishedVersion

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.

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. The study assembled a chromosome-level genome of Cycas panzhihuaensis, the last major lineage of seed plants for which a high-quality genome assembly was lacking. The study closes an important gap in our understanding of genome structure and evolution in seed plants.

STT3 is a catalytic subunit of hetero‐oligomeric oligosaccharyltransferase (OST), which is important for asparagine‐linked glycosylation. In mammals and plants, OSTs with different STT3 isoforms exhibit distinct levels of enzymatic efficiency or different responses to stressors. Although two different STT3 isoforms have been identified in both plants and animals, it remains unclear whether these isoforms result from gene duplication in an ancestral eukaryote. Furthermore, the molecular mechanisms underlying the functional divergences between the two STT3 isoforms in plant have not been well elucidated. Here, we conducted phylogenetic analysis of the major evolutionary node species and suggested that gene duplications of STT3 may have occurred independently in animals and plants. Across land plants, the exon–intron structure differed between the two STT3 isoforms, but was highly conserved for each isoform. Most angiosperm STT3a genes had 23 exons with intron phase 0, while STT3b genes had 6 exons with intron phase 2. Characteristic motifs (motif 18 and 19) of STT3s were mapped to different structure domains in the plant STT3 proteins. These two motifs overlap with regions of high nonsynonymous‐to‐synonymous substitution rates, suggesting the regions may be related to functional difference between STT3a and STT3b. In addition, promoter elements and gene expression profiles were different between the two isoforms, indicating expression pattern divergence of the two genes. Collectively, the identified differences may result in the functional divergence of plant STT3s., STT3 is a catalytic subunit of oligosaccharyltransferase (OST) which is important for N‐glycosylation. Two subtypes of STT3 exist in animal and plant through independent duplications. An ancient separation of STT3 in the common ancestor of plants was discovered. Different elements in the promoters, distinct motifs, and auxiliary subunits may confer divergent catalytic ability on OST containing STT3a or STT3b.

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.

Nucleotide-binding site-leucine-rich repeat receptor (NLR) genes comprise the largest family of plant disease resistance genes. 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 plant lifestyle, or between NLR content and signal transduction components, are poorly characterized due to limited genome availability. Here, we established an angiosperm NLRatlas (ANNA, http://compbio.nju.edu.cn/app/ANNA/), which 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 due 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 long-term evolutionary silence of NLR genes in 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 recovered a conserved TNL lineage that may function independently of the RNL pathway. Our findings provide new insights into the evolution of NLR genes in the context of plant lifestyles and genome content variation.

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.

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.

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.

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.

Soybean mosaic virus (SMV) is widely recognized as a highly damaging pathogen of soybean, and various strains/isolates have been reported to date. However, the pathogenic differences and phylogenetic relationships of these SMV strains/isolates have not been extensively studied. In the present work, by first obtaining 18 new genomic sequences of Chinese SMV strains/isolates and further compiling these with available data, we have explored the evolution of SMV from multiple aspects. First, as in other potyviruses, recombination has occurred frequently during SMV evolution, and a total of 32 independent events were detected. Second, using a maximum-likelihood method and removing recombinant fragments, a phylogeny covering 83 SMV sequences sampled from all over the world was reconstructed and the results showed four separate SMV clades, with clade I and II recovered for the first time. Third, the population structure analysis of SMV revealed significant genetic differentiations between China and two other countries (Korea and U.S.A.). Fourth, certain SMV-encoded genes, such as P1, HC-Pro and P3, exhibited higher non-synonymous substitution rate (dN) than synonymous substitution rate (dS), indicating that positive selection has influenced these genes. Finally, four Chinese SMV strains/isolates were selected for inoculation of both USA and Chinese differential soybean cultivars, and their pathogenic phenotypes were significantly different from that of the American strains. Overall, these findings have further broadened our understanding on SMV evolution, which would assist researchers to better deal with this harmful virus.

Plant genomes harbor dozens to hundreds of nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes; however, the long-term evolutionary history of these resistance genes has not been fully understood. This study focuses on five Brassicaceae genomes and the Carica papaya genome to explore changes in NBS-LRR genes that have taken place in this Rosid II lineage during the past 72 million years. Various numbers of NBS-LRR genes were identified from Arabidopsis lyrata (198), A. thaliana (165), Brassica rapa (204), Capsella rubella (127), Thellungiella salsuginea (88), and C. papaya (51). In each genome, the identified NBS-LRR genes were found to be unevenly distributed among chromosomes and most of them were clustered together. Phylogenetic analysis revealed that, before and after Brassicaceae speciation events, both toll/interleukin-1 receptor-NBS-LRR (TNL) genes and non-toll/interleukin-1 receptor-NBS-LRR (nTNL) genes exhibited a pattern of first expansion and then contraction, suggesting that both subclasses of NBS-LRR genes were responding to pathogen pressures synchronically. Further, by examining the gain/loss of TNL and nTNL genes at different evolutionary nodes, this study revealed that both events often occurred more drastically in TNL genes. Finally, the phylogeny of nTNL genes suggested that this NBS-LRR subclass is composed of two separate ancient gene types: RPW8-NBS-LRR and Coiled-coil-NBS-LRR.

Streptococcal histidine triad proteins HTPs are widely distributed within the Streptococcus genus. Based on the phylogenetic relationship and domain composition, HTPs are classified into type I and type II subfamilies. Previous studies revealed that several pathogenic streptococci contain more than one htp gene. We found that the highly virulent strain of Streptococcus suis 2 (S. suis 2), 05ZYH33 encodes 3 HTPs, designated HtpsA (previously described as HtpS), HtpsB, and HtpsC. Among them, HtpsC is the only member that contains leucine-rich repeat (LRR) domains at the C-terminal. In this study, we demonstrated that the recombinant HtpsC could bind to 2 different components of human ECM complex laminin and fibronectin in vitro, suggesting that it is a novel adhesin of S. suis 2. Having constructed an htpsC mutant, we evaluated its role in the pathogenesis of the highly virulent S. suis 2 strain 05ZYH33. Our data showed that inactivation of htpsC significantly affected adherence of S. suis 2 to Hep-2 cells and shortened the survival of the bacteria in whole blood. Furthermore, deletion of htpsC significantly attenuated the virulence of S. suis 2 in mice. These results demonstrated that htpsC was involved in the pathogenesis of the highly virulent S. suis 2 strain 05ZYH33. In line with the observation, immunization with HtpsC significantly prolonged mice's survival after S. suis 05ZYH33 challenge, indicating its potential use in the vaccine development against S. suis.

Widely known as a severe pathogen of bean plants, the bean common mosaic virus (BCMV) has been reported to infect soybeans only sporadically and the involved strains were all found in China regions. To explore variations among soybean-infecting BCMV strains, hundreds of soybean mosaic leave samples were collected throughout China, with a total of 30 BCMV isolates detected and their genomes sequenced. These newly obtained genomes, together with 16 other BCMV genomes available in GenBank were examined from multiple aspects to characterize BCMV evolutionary processes. Phylogenetic analysis showed that both soybean-infecting BCMVs (group I) and peanut-infecting BCMVs (group II) are distantly related to other BCMVs, suggesting ancestral differentiation and host adaptation. Genetic variation analysis showed that P1, P3 and 6K2 genes and the beginning portion of CP gene showed higher levels of variation relative to other genes. Moreover, selection analyses further confirmed that a number of sites within the P1 and P3 genes have suffered positive selection. These obtained BCMV sequences also exhibit high recombination frequencies, indicating a more dynamic evolutionary history. Finally, 12 different soybean cultivars were challenged with two BCMV isolates (DXH015 and HZZB011), with most of the cultivars successfully infected. These findings suggest that BCMV is indeed a potential threat to soybean production.

Proper utilization of plant disease resistance genes requires a good understanding of their short- and long-term evolution. Here we present a comprehensive study of the long-term evolutionary history of nucleotide-binding site (NBS)-leucine-rich repeat (LRR) genes within and beyond the legume family. The small group of NBS-LRR genes with an amino-terminal RESISTANCE TO POWDERY MILDEW8 (RPW8)-like domain (referred to as RNL) was first revealed as a basal clade sister to both coiled-coil-NBS-LRR (CNL) and Toll/Interleukin1 receptor-NBS-LRR (TNL) clades. Using Arabidopsis (Arabidopsis thaliana) as an outgroup, this study explicitly recovered 31 ancestral NBS lineages (two RNL, 21 CNL, and eight TNL) that had existed in the rosid common ancestor and 119 ancestral lineages (nine RNL, 55 CNL, and 55 TNL) that had diverged in the legume common ancestor. It was shown that, during their evolution in the past 54 million years, approximately 94% (112 of 119) of the ancestral legume NBS lineages experienced deletions or significant expansions, while seven original lineages were maintained in a conservative manner. The NBS gene duplication pattern was further examined. The local tandem duplications dominated NBS gene gains in the total number of genes (more than 75%), which was not surprising. However, it was interesting from our study that ectopic duplications had created many novel NBS gene loci in individual legume genomes, which occurred at a significant frequency of 8% to 20% in different legume lineages. Finally, by surveying the legume microRNAs that can potentially regulate NBS genes, we found that the microRNA-NBS gene interaction also exhibited a gain-and-loss pattern during the legume evolution.

Streptococcus suis (S. suis) is an important zoonotic pathogen that causes multiple diseases in both pigs and humans. Many studies suggest that Streptococcus utilizes host extracellular matrix proteins, including laminin, for adhesion and invasion of host cells. Recently, we identified a putative Lmb protein (CDS 0330) of a highly virulent strain of S. suis (serotype 2). In this study, we characterized the ability of CDS 0330 to bind human laminin, and evaluated the protective efficacy of a recombinant protein vaccine. Bioinformatic analysis revealed that both the amino acid sequence and tertiary structure of CDS 0330 were similar to Lmb proteins in other Streptococcus. In addition, the sequence of CDS 0330 was present in the genomes of 26 of the 38 sequenced streptococci species, indicating an early origin and conservation of this gene. Particularly, all 17 sequenced S. suis genomes, regardless of serotype or geographic origin, contained CDS 0330 gene in their genome with a minimum pair-wise amino acid identity of 92%. PCR amplification revealed that CDS 0330 gene is distributed throughout 35 S. suis serotypes in the lmb-htp format. Flow cytometry analysis confirmed that CDS 0330 was expressed on the cell surface of S. suis, and ELISA revealed the recombinant CDS 0330 protein could bind laminin in vitro. Finally, vaccinating mice with recombinant CDS 0330 protein significantly prolonged survival after S. suis infection. Together, these data reveal that CDS 0330 is a laminin binding protein of S. suis 2, and open new avenues for preventing S. suis 2 infection.

Nucleotide-Binding Site-Leucine-Rich Repeat (NBS-LRR) genes are the largest plant disease resistance (R) gene family, accounting for ∼80% of more than 140 cloned R genes. Recently, we systematically investigated NBS-LRR genes in 22 angiosperm genomes. By performing phylogenetic analysis of these genes in major angiosperm clades separately and as a whole, we gained strong evidence supporting that angiosperm NBS-LRR genes are derived from 3 anciently separated NBS-LRR classes: RPW8-NBS-LRR (RNL), TIR-NBS-LRR (TNL) and CC-NBS-LRR (CNL). A total of 23 ancestral NBS-LRR lineages gave rise to the current NBS-LRR diversity of angiosperm through dynamic expansions. Comparative analysis of RNL, TNL, and CNL classes revealed that while RNL genes evolved conservatively to maintain its role in defense signal transduction, the latter 2 classes underwent convergent recent expansions in various plant genomes. The revealed evolutionary pattern of angiosperm NBS-LRR genes reflects a long history of competition between plant and pathogen.

The Rsv1 - h gene in cultivar Suweon 97, which confers resistance to SMVs, was mapped to a 97.5-kb location (29,815,195–29,912,667 bp on chromosome 13) in the Rsv1 locus, thereby providing additional insights into the molecular nature underlying variations in resistance alleles in this particular locus. Soybean mosaic virus (SMV) is a well-known devastating pathogen of soybean (Glycine max (L.) Merrill.) causing significant yield losses and seed quality deterioration. A single dominant allele, Rsv1-h, which confers resistance to multiple SMV strains, was previously reported in the cultivar Suweon 97, but its exact location is unknown. In the present study, Suweon 97 was crossed with a SMV-sensitive cultivar, Williams 82. Inoculating 267 F 2 individuals with two Chinese SMV strains (SC6-N and SC7-N) demonstrated that one single dominant gene confers SMV resistance. Another 1,150 F 2 individuals were then screened for two simple sequence repeat (SSR) markers (BARCSOYSSR_13_1103 and BARCSOYSSR_13_1187) that flank the Rsv1 locus. Seventy-four recombinants were identified and 20 additional polymorphic SSR markers within the Rsv1 region were then employed in genotyping these recombinants. F 2:3 and F 3:4 recombinant lines were also inoculated with SC6-N and SC7-N to determine their phenotypes. The final data revealed that in Suweon 97, the Rsv1-h gene that confers resistance to SC6-N and SC7-N was flanked by BARCSOYSSR_13_1114 and BARCSOYSSR_13_1115, two markers that delimit a 97.5-kb region in the reference Williams 82 genome. In such region, eight genes were present, of which two, Glyma13g184800 and Glyma13g184900, encode the characteristic CC-NBS-LRR type of resistance gene and were considered potential candidates for Rsv1-h.

A majority of land plants can form symbiosis with arbuscular mycorrhizal (AM) fungi. MicroRNAs (miRNAs) have been implicated to regulate this process in legumes, but their involvement in non-legume species is largely unknown. In this study, by performing deep sequencing of sRNA libraries in tomato roots and comparing with tomato genome, a total of 700 potential miRNAs were predicted, among them, 187 are known plant miRNAs that have been previously deposited in miRBase. Unlike the profiles in other plants such as rice and Arabidopsis, a large proportion of predicted tomato miRNAs was 24 nt in length. A similar pattern was observed in the potato genome but not in tobacco, indicating a Solanum genus-specific expansion of 24-nt miRNAs. About 40% identified tomato miRNAs showed significantly altered expressions upon Rhizophagus irregularis inoculation, suggesting the potential roles of these novel miRNAs in AM symbiosis. The differential expression of five known and six novel miRNAs were further validated using qPCR analysis. Interestingly, three up-regulated known tomato miRNAs belong to a known miR171 family, a member of which has been reported in Medicago truncatula to regulate AM symbiosis. Thus, the miR171 family likely regulates AM symbiosis conservatively across different plant lineages. More than 1000 genes targeted by potential AM-responsive miRNAs were provided and their roles in AM symbiosis are worth further exploring.

Streptococcal histidine triad protein was identified recently as a cell surface-associated protein family. Five members of this family (PhtA, PhtB, PhtD, PhtE and HtpA), derived from Streptococcus pneumoniae and Streptococcus pyogenes, have been shown as antigens that confer protection to the host on infection. In this report, a gene sequence highly homologous to htpA and phtD (designated htpS, the histidine triad protein of Streptococcus suis) was identified from S. suis 2 Chinese strain 05ZYH33. Our data revealed that htpS is extremely conserved in S. suis 2 and widely distributed in 83% (29/35) of 35 S. suis serotypes. It was also demonstrated by Western blot and flow cytometry that HtpS is a cell surface-associated protein that was expressed during S. suis 2 infection. An antibody against HtpS could increase the deposition of human complement 3 on S. suis 2 and also enhance the clearance of S. suis 2 in whole blood. In addition, mice could be immunized against S. suis 2 infection and were well protected after immunization with recombinant HtpS.

Plant genomes harbor dozens to hundreds of nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes; however, the long-term evolutionary history of these resistance genes has not been fully understood. This study focuses on five Brassicaceae genomes and the Carica papaya genome to explore changes in NBS-LRR genes that have taken place in this Rosid II lineage during the past 72 million years. Various numbers of NBS-LRR genes were identified from Arabidopsis lyrata (198), A. thaliana (165), Brassica rapa (204), Capsella rubella (127), Thellungiella salsuginea (88), and C. papaya (51). In each genome, the identified NBS-LRR genes were found to be unevenly distributed among chromosomes and most of them were clustered together. Phylogenetic analysis revealed that, before and after Brassicaceae speciation events, both toll/interleukin-1 receptor-NBS-LRR (TNL) genes and non-toll/interleukin-1 receptor-NBS-LRR (nTNL) genes exhibited a pattern of first expansion and then contraction, suggesting that both subclasses of NBS-LRR genes were responding to pathogen pressures synchronically. Further, by examining the gain/loss of TNL and nTNL genes at different evolutionary nodes, this study revealed that both events often occurred more drastically in TNL genes. Finally, the phylogeny of nTNL genes suggested that this NBS-LRR subclass is composed of two separate ancient gene types: RPW8-NBS-LRR and Coiled-coil-NBS-LRR.

A genome triplication took place in the ancestor of Brassiceae species after the split of the Arabidopsis lineage. The postfragmentation and shuffling of the genome turned the ancestral hexaploid back to diploids and caused the radiation of Brassiceae species. The course of speciation was accompanied by the loss of duplicate genes and also influenced the evolution of retained genes. Of all the genes, those encoding NBS domains are typical R genes that confer resistance to invading pathogens. In this study, using the genome of Arabidopsis thaliana as a reference, we examined the loss/retention of orthologous NBS-encoding loci in the tripled Brassica rapa genome and discovered differential loss/retention frequencies. Further analysis indicated that loci of different retention ratios showed different evolutionary patterns. The loci of classesII and III (maintaining two and three syntenic loci, respectively, multi-loci) show sharper expansions by tandem duplications, have faster evolutionary rates and have more potential to be associated with novel gene functions. On the other hand, the loci that are retained at the minimal rate (keeping only one locus, class I, single locus) showed opposite patterns. Phylogenetic analysis indicated that recombination and translocation events were common among multi-loci in B. rapa, and differential evolutionary patterns between multi- and single-loci are likely the consequence of recombination. Investigations towards other gene families demonstrated different evolutionary characteristics between different gene families. The evolution of genes is more likely determined by the property of each gene family, and the whole genome triplication provided only a specific condition.

Non-random arrangement of synonymous codons in coding sequences has been recently reported in eukaryotic and bacterial genomes, but the case in archaeal genomes is largely undetermined. Here, we systematically investigated 122 archaeal genomes for their synonymous codon co-occurrence patterns. We found that in most archaeal coding sequences, the order of synonymous codons is not arranged randomly, but rather some successive codon pairs appear significantly more often than expected. Importantly, such codon pairing bias (CPB) pattern in archaea does not seem to completely follow the co-tRNA codon pairing (CCP) rule previously reported for eukaryotes, but largely obeys an identical codon pairing (ICP) rule. Further, synonymous codon permutation test demonstrated that in many archaeal genomes, random mutation alone is unable to cause the observed high level of ICP bias, which strongly indicates that selection force has been involved to shape synonymous codon orders, potentially meeting a global requirement to optimize translation rate.