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.

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.

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 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.

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.

Divergence and Conservative Evolution of XTNX Genes in Land Plants

The Toll-interleukin-1 receptor (TIR) and Nucleotide-binding site (NBS) domains are two major components of the TIR-NBS-leucine-rich repeat family plant disease resistance genes. Extensive functional and evolutionary studies have been performed on these genes; however, the characterization of a small group of genes that are composed of atypical TIR and NBS domains, namely XTNX genes, is limited. The present study investigated this specific gene family by conducting genome-wide analyses of 59 green plant genomes. A total of 143 XTNX genes were identified in 51 of the 52 land plant genomes, whereas no XTNX gene was detected in any green algae genomes, which indicated that XTNX genes originated upon emergence of land plants. Phylogenetic analysis revealed that the ancestral XTNX gene underwent two rounds of ancient duplications in land plants, which resulted in the formation of clades I/II and clades IIa/IIb successively. Although clades I and IIb have evolved conservatively in angiosperms, the motif composition difference and sequence divergence at the amino acid level suggest that functional divergence may have occurred since the separation of the two clades. In contrast, several features of the clade IIa genes, including the absence in the majority of dicots, the long branches in the tree, the frequent loss of ancestral motifs, and the loss of expression in all detected tissues of Zea mays, all suggest that the genes in this lineage might have undergone pseudogenization. This study highlights that XTNX genes are a gene family originated anciently in land plants and underwent specific conservative pattern in evolution.

Fine mapping of the Rsv1-h gene in the soybean cultivar Suweon 97 that confers resistance to two Chinese strains of the soybean mosaic virus

KEY MESSAGE

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.

Tracking ancestral lineages and recent expansions of NBS-LRR genes in angiosperms

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 type II histidine triad protein HtpsC is a novel adhesion with the involvement of Streptococcus suis virulence

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.

Large-Scale Analyses of Angiosperm Nucleotide-Binding Site-Leucine-Rich Repeat Genes Reveal Three Anciently Diverged Classes with Distinct Evolutionary Patterns

Nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes make up the largest plant disease resistance gene family (R genes), with hundreds of copies occurring in individual angiosperm genomes. However, the expansion history of NBS-LRR genes during angiosperm evolution is largely unknown. By identifying more than 6,000 NBS-LRR genes in 22 representative angiosperms and reconstructing their phylogenies, we present a potential framework of NBS-LRR gene evolution in the angiosperm. Three anciently diverged NBS-LRR classes (TNLs, CNLs, and RNLs) were distinguished with unique exon-intron structures and DNA motif sequences. A total of seven ancient TNL, 14 CNL, and two RNL lineages were discovered in the ancestral angiosperm, from which all current NBS-LRR gene repertoires were evolved. A pattern of gradual expansion during the first 100 million years of evolution of the angiosperm clade was observed for CNLs. TNL numbers remained stable during this period but were eventually deleted in three divergent angiosperm lineages. We inferred that an intense expansion of both TNL and CNL genes started from the Cretaceous-Paleogene boundary. Because dramatic environmental changes and an explosion in fungal diversity occurred during this period, the observed expansions of R genes probably reflect convergent adaptive responses of various angiosperm families. An ancient whole-genome duplication event that occurred in an angiosperm ancestor resulted in two RNL lineages, which were conservatively evolved and acted as scaffold proteins for defense signal transduction. Overall, the reconstructed framework of angiosperm NBS-LRR gene evolution in this study may serve as a fundamental reference for better understanding angiosperm NBS-LRR genes.

Identification of Arbuscular Mycorrhiza (AM)-Responsive microRNAs in Tomato

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.

Uncovering the dynamic evolution of nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes in Brassicaceae

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.

The evolution of soybean mosaic virus: An updated analysis by obtaining 18 new genomic sequences of Chinese strains/isolates

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.

A genomic survey of thirty soybean-infecting bean common mosaic virus (BCMV) isolates from China pointed BCMV as a potential threat to soybean production

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.

Long-term evolution of nucleotide-binding site-leucine-rich repeat genes: understanding gained from and beyond the legume family

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.

Non-random arrangement of synonymous codons in archaea coding sequences

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.

Insight into the evolution of the histidine triad protein (HTP) family in Streptococcus

The Histidine Triad Proteins (HTPs), also known as Pht proteins in Streptococcus pneumoniae, constitute a family of surface-exposed proteins that exist in many pathogenic streptococcal species. Although many studies have revealed the importance of HTPs in streptococcal physiology and pathogenicity, little is known about their origin and evolution. In this study, after identifying all htp homologs from 105 streptococcal genomes representing 38 different species/subspecies, we analyzed their domain structures, positions in genome, and most importantly, their evolutionary histories. By further projecting this information onto the streptococcal phylogeny, we made several major findings. First, htp genes originated earlier than the Streptococcus genus and gene-loss events have occurred among three streptococcal groups, resulting in the absence of the htp gene in the Bovis, Mutans and Salivarius groups. Second, the copy number of htp genes in other groups of Streptococcus is variable, ranging from one to four functional copies. Third, both phylogenetic evidence and domain structure analyses support the division of two htp subfamilies, designated as htp I and htp II. Although present mainly in the pyogenic group and in Streptococcus suis, htp II members are distinct from htp I due to the presence of an additional leucine-rich-repeat domain at the C-terminus. Finally, htp genes exhibit a faster nucleotide substitution rate than do housekeeping genes. Specifically, the regions outside the HTP domains are under strong positive selection. This distinct evolutionary pattern likely helped Streptococcus to easily escape from recognition by host immunity.

Synonymous codon ordering: a subtle but prevalent strategy of bacteria to improve translational efficiency

BACKGROUND

In yeast coding sequences, once a particular codon has been used, subsequent occurrence of the same amino acid tends to use codons sharing the same tRNA. Such a phenomenon of co-tRNA codons pairing bias (CTCPB) is also found in some other eukaryotes but it is not known whether it occurs in prokaryotes.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we focused on a total of 773 bacterial genomes to investigate their synonymous codon pairing preferences. After calculating the actual frequencies of synonymous codon pairs and comparing them with their expected values, we detected an obvious pairing bias towards identical codon pairs. This seems consistent with the previously reported CTCPB phenomenon, since identical codons are certainly read by the same tRNA. However, among co-tRNA but non-identical codon pairs, only 22 were often found overrepresented, suggesting that many co-tRNA codons actually do not preferentially pair together in prokaryotes. Therefore, the previously reported co-tRNA codons pairing rule needs to be more rigorously defined. The affinity differences between a tRNA anticodon and its readable codons should be taken into account. Moreover, both within-gene-shuffling tests and phylogenetic analyses support the idea that translational selection played an important role in shaping the observed synonymous codon pairing pattern in prokaryotes.

CONCLUSIONS

Overall, a high level of synonymous codon pairing bias was detected in 73% investigated bacterial species, suggesting the synonymous codon ordering strategy has been prevalently adopted by prokaryotes to improve their translational efficiencies. The findings in this study also provide important clues to better understand the complex dynamics of translational process.

Optimal codon identities in bacteria: implications from the conflicting results of two different methods

A correlation method was recently adopted to identify selection-favored 'optimal' codons from 675 bacterial genomes. Surprisingly, the identities of these optimal codons were found to track the bacterial GC content, leading to a conclusion that selection would generally shape the codon usages to the same direction as the overall mutation does. Raising several concerns, here we report a thorough comparative study on 203 well-selected bacterial species, which strongly suggest that the previous conclusion is likely an illusion. Firstly, the previous study did not preclude species that are suffering weak or no selection pressures on their codon usages. For these species, as showed in this study, the optimal codon identities are prone to be incorrect and follow GC content. Secondly, the previous study only adopted the correlation method, without considering another method to test the reliability of inferred optimal codons. Actually by definition, optimal codons can also be identified by simply comparing codon usages between high- and low-expression genes. After using both methods to identify optimal codons for the selected species, we obtained highly conflicting results, suggesting at least one method is misleading. Further we found a critical problem of correlation method at the step of calculating gene bias level. Due to a failure of accurately defining the background mutation, the problem would result in wrong optimal codon identities. In other words, partial mutational effects on codon choices were mistakenly regarded as selective influences, leading to incorrect and biased optimal codon identities. Finally, considering the translational dynamics, optimal codons identified by comparison method can be well-explained by tRNA compositions, whereas optimal codons identified by correlation method can not be. For all above reasons, we conclude that real optimal codons actually do not track the genomic GC content, and correlation method is misleading in identifying optimal codons and better be avoided.

Metabolic regulation of lrp gene expression in Escherichia coli K-12

Expression of the lrp gene is regulated in part by the nutrients available to the cell, and is decreased in rich medium, in glucose minimal media enriched with amino acids, and in minimal medium with alternative carbon sources, such as acetate and succinate. When Lrp production is increased in a given medium, expression of its target genes is also increased. However, when the medium is changed from glucose to acetate, the response of the target genes is governed by many factors.

Characterization of the expression of the thcB gene, coding for a pesticide-degrading cytochrome P-450 in Rhodococcus strains

A cytochrome P-450 system in Rhodococcus strains, encoded by thcB, thcC, and thcD, participates in the degradation of thiocarbamates and several other pesticides. The regulation of the system was investigated by fusing a truncated lacZ in frame to thcB, the structural gene for the cytochrome P-450 monooxygenase. Analysis of the thcB-lacZ fusion showed that the expression of thcB was 10-fold higher in the presence of the herbicide EPTC (s-ethyl dipropylthiocarbamate). Similar enhancement of the thcB-lacZ expression was found with other thiocarbamate pesticides. Atrazine, simazine, or carbofuran, although metabolized by the system, had no effect on the thcB-lacZ expression. The presence of glucose slightly increased the expression of thcB-lacZ, indicating no catabolic repression of the thcB-lacZ expression. The expression of thcB-lacZ was decreased more than twofold in Luria-Bertani medium. This was due in part to cysteine, which repressed thcB-lacZ expression. It was confirmed that the thcR gene, which is transcribed divergently from thcB, codes for a positive regulatory protein which is essential for the thcB-lacZ expression. Studies of the thcR-lacZ protein fusion showed that the thcR gene is expressed constitutively.

Cloning and expression of the s-triazine hydrolase gene (trzA) from Rhodococcus corallinus and development of Rhodococcus recombinant strains capable of dealkylating and dechlorinating the herbicide atrazine

We used degenerate oligodeoxyribonucleotides derived from the N-terminal sequence of the s-triazine hydrolase from Rhodococcus corallinus NRRL B-15444R in an amplification reaction to isolate a DNA segment containing a 57-bp fragment from the trzA gene. By using the nucleotide sequence of this fragment, a nondegenerate oligodeoxyribonucleotide was synthesized and used to screen a genomic library of R. corallinus DNA for fragments containing trzA. A 5.3-kb PstI fragment containing trzA was cloned, and the nucleotide sequence of a 2,450-bp region containing trzA was determined. No trzA expression was detected in Escherichia coli or several other gram-negative bacteria. The trzA gene was subcloned into a Rhodococcus-E. coli shuttle vector, pBS305, and transformed into several Rhodococcus strains. Expression of trzA was demonstrated in all Rhodococcus transformants. Rhodococcus sp. strain TE1, which possesses the catabolic gene (atrA) for the N-dealkylation of the herbicides atrazine and simazine, was able to dechlorinate the dealkylated metabolites of atrazine and simazine when carrying the trzA gene on a plasmid. A plasmid carrying both atrA and trzA was constructed and transformed into three atrA- and trzA-deficient Rhodococcus strains. Both genes were expressed in the transformants. The s-triazine hydrolase activity of the recombinant strains carrying the trzA plasmid were compared with that of the R. corallinus strain from which it was derived.

Cloning of the genes for degradation of the herbicides EPTC (S-ethyl dipropylthiocarbamate) and atrazine from Rhodococcus sp. strain TE1

The degradation of the herbicides EPTC (S-ethyl dipropylthiocarbamate) and atrazine (2-chloro-4-ethyl-amino-6-isopropylamino-1,3,5-triazine) is associated with an indigenous plasmid in Rhodococcus sp. strain TE1. Plasmid DNA libraries of Rhodococcus sp. strain TE1 were constructed in a Rhodococcus-Escherichia coli shuttle vector, pBS305, and transferred into Rhodococcus sp. strain TE3, a derivative of Rhodococcus sp. strain TE1 lacking herbicide degradation activity, to select transformants capable of growing on EPTC as the sole source of carbon (EPTC+). Analysis of plasmids from the EPTC+ transformants indicated that the eptA gene, which codes for the enzyme required for EPTC degradation, residues on a 6.2-kb KpnI fragment. The cloned fragment also harbored the gene required for atrazine N dealkylation (atrA). The plasmid carrying the cloned fragment could be electroporated into a number of other Rhodococcus strains in which both eptA and atrA were fully expressed. No expression of the cloned genes was evident in E. coli strains. Subcloning of the 6.2-kb fragment to distinguish between EPTC- and atrazine-degrading genes was not successful.