Analysis of Small RNAs in Streptococcus mutans under Acid Stress—A New Insight for Caries Research

Streptococcus mutans (S. mutans) is the major clinical pathogen responsible for dental caries. Its acid tolerance has been identified as a significant virulence factor for its survival and cariogenicity in acidic conditions. Small RNAs (sRNAs) are recognized as key regulators of virulence and stress adaptation. Here, we constructed three libraries of sRNAs with small size exposed to acidic conditions for the first time, followed by verification using qRT-PCR. The levels of two sRNAs and target genes predicted to be bioinformatically related to acid tolerance were further evaluated under different acid stress conditions (pH 7.5, 6.5, 5.5, and 4.5) at three time points (0.5, 1, and 2 h). Meanwhile, bacterial growth characteristics and vitality were assessed. We obtained 1879 sRNAs with read counts of at least 100. One hundred and ten sRNAs were perfectly mapped to reported msRNAs in S. mutans. Ten out of 18 sRNAs were validated by qRT-PCR. The survival of bacteria declined as the acid was increased from pH 7.5 to 4.5 at each time point. The bacteria can proliferate under each pH except pH 4.5 with time. The levels of sRNAs gradually decreased from pH 7.5 to 5.5, and slightly increased in pH 4.5; however, the expression levels of target mRNAs were up-regulated in acidic conditions than in pH 7.5. These results indicate that some sRNAs are specially induced at acid stress conditions, involving acid adaptation, and provide a new insight into exploring the complex acid tolerance for S. mutans.     

Glucosinolate and Trichome Defenses in a Natural Arabidopsis lyrata Population

Glucosinolates (GS) and trichomes contribute to plant resistance against insect herbivores in the model Arabidopsis thaliana. The functional and genetic characteristics of herbivore defense, however, can differ even between closely related species. In a quantitative genetic experiment with the out-crossing perennial Arabidopsis lyrata spp. petraea, we measured constitutive GS composition, trichome density, leaf thickness, and plant resistance in four different herbivore interactions. In a single population of A. lyrata, we found heritable variation for trichome density as well as GS amount and carbon side-chain elongation ratios associated with activity in methylthioalkylmalate synthase (MAM). Unexpectedly, heritabilities for indole GS in A. lyrata were high and less affected by differences in plant age and environment than aliphatic GS. We found significant heritability in plant resistance to the specialist Plutella xylostella and generalist Trichoplusia ni, but not to the specialists Pieris brassicae and Phyllotreta cruciferae. Analyses of phenotypic and genetic correlations between candidate defense traits and insect resistance suggested that A. lyrata resistance was conferred by a combination of indole GS amount and trichome density, and, to a lesser extent, aliphatic GS ratios and leaf thickness. Variation in the most abundant compound, the aliphatic 3-hydroxypropyl GS, had little impact on A. lyrata herbivore resistance. The contribution of defense traits to resistance depended on the experimental herbivory context, and resistances were weakly correlated. A diversified defense strategy is likely to be important for long-lived individuals of A. lyrata that are subject to attack by many different herbivores in nature.     

The QseEF Two-Component System-GlmY Small RNA Regulatory Pathway Controls Swarming in Uropathogenic Proteus mirabilis

Bacterial sensing of environmental signals through the two-component system (TCS) plays a key role in modulating virulence. In the search for the host hormone-sensing TCS, we identified a conserved qseEGF locus following glmY, a small RNA (sRNA) gene in uropathogenic Proteus mirabilis. Genes of glmY-qseE-qseG-qseF constitute an operon, and QseF binding sites were found in the glmY promoter region. Deletion of glmY or qseF resulted in reduced swarming motility and swarming-related phenotypes relative to the wild-type and the respective complemented strains. The qseF mutant had decreased glmYqseEGF promoter activity. Both glmY and qseF mutants exhibited decreased flhDC promoter activity and mRNA level, while increased rcsB mRNA level was observed in both mutants. Prediction by TargetRNA2 revealed cheA as the target of GlmY. Then, construction of the translational fusions containing various lengths of cheA 5′UTR for reporter assay and site-directed mutagenesis were performed to investigate the cheA-GlmY interaction in cheA activation. Notably, loss of glmY reduced the cheA mRNA level, and urea could inhibit swarming in a QseF-dependent manner. Altogether, this is the first report elucidating the underlying mechanisms for modulation of swarming motility by a QseEF-regulated sRNA GlmY, involving expression of cheA, rcsB and flhDC in uropathogenic P. mirabilis.     

Identification and Evolution of Functional Alleles of the Previously Described Pollen Specific Myrosinase Pseudogene AtTGG6 in Arabidopsis thaliana

Myrosinases are β-thioglucoside glucohydrolases and serve as defense mechanisms against insect pests and pathogens by producing toxic compounds. AtTGG6 in Arabidopsis thaliana was previously reported to be a myrosinase pseudogene but specifically expressed in pollen. However, we found that AlTGG6, an ortholog to AtTGG6 in A. lyrata (an outcrossing relative of A. thaliana) was functional, suggesting that functional AtTGG6 alleles may still exist in A. thaliana. AtTGG6 alleles in 29 A. thaliana ecotypes were cloned and sequenced. Results indicate that ten alleles were functional and encoded Myr II type myrosinase of 512 amino acids, and myrosinase activity was confirmed by overexpressing AtTGG6 in Pichia pastoris. However, the 19 other ecotypes had disabled alleles with highly polymorphic frame-shift mutations and diversified sequences. Thirteen frame-shift mutation types were identified, which occurred independently many times in the evolutionary history within a few thousand years. The functional allele was expressed specifically in pollen similar to the disabled alleles but at a higher expression level, suggesting its role in defense of pollen against insect pests such as pollen beetles. However, the defense function may have become less critical after A. thaliana evolved to self-fertilization, and thus resulted in loss of function in most ecotypes.     

smartPARE: An R Package for Efficient Identification of True mRNA Cleavage Sites

Degradome sequencing is commonly used to generate high-throughput information on mRNA cleavage sites mediated by small RNAs (sRNA). In our datasets of potato (Solanum tuberosum, St) and Phytophthora infestans (Pi), initial predictions generated high numbers of cleavage site predictions, which highlighted the need of improved analytic tools. Here, we present an R package based on a deep learning convolutional neural network (CNN) in a machine learning environment to optimize discrimination of false from true cleavage sites. When applying smartPARE to our datasets on potato during the infection process by the late blight pathogen, 7.3% of all cleavage windows represented true cleavages distributed on 214 sites in P. infestans and 444 sites in potato. The sRNA landscape of the two organisms is complex with uneven sRNA production and cleavage regions widespread in the two genomes. Multiple targets and several cases of complex regulatory cascades, particularly in potato, was revealed. We conclude that our new analytic approach is useful for anyone working on complex biological systems and with the interest of identifying cleavage sites particularly inferred by sRNA classes beyond miRNAs.     

Conserved Opposite Functions in Plant Resistance to Biotrophic and Necrotrophic Pathogens of the Immune Regulator SRFR1

Plant immunity is mediated in large part by specific interactions between a host resistance protein and a pathogen effector protein, named effector-triggered immunity (ETI). ETI needs to be tightly controlled both positively and negatively to enable normal plant growth because constitutively activated defense responses are detrimental to the host. In previous work, we reported that mutations in SUPPRESSOR OF rps4-RLD1 (SRFR1), identified in a suppressor screen, reactivated EDS1-dependent ETI to Pseudomonas syringae pv. tomato (Pto) DC3000. Besides, mutations in SRFR1 boosted defense responses to the generalist chewing insect Spodoptera exigua and the sugar beet cyst nematode Heterodera schachtii. Here, we show that mutations in SRFR1 enhance susceptibility to the fungal necrotrophs Fusarium oxysporum f. sp. lycopersici (FOL) and Botrytis cinerea in Arabidopsis. To translate knowledge obtained in AtSRFR1 research to crops, we generated SlSRFR1 alleles in tomato using a CRISPR/Cas9 system. Interestingly, slsrfr1 mutants increased expression of SA-pathway defense genes and enhanced resistance to Pto DC3000. In contrast, slsrfr1 mutants elevated susceptibility to FOL. Together, these data suggest that SRFR1 is functionally conserved in both Arabidopsis and tomato and functions antagonistically as a negative regulator to (hemi-) biotrophic pathogens and a positive regulator to necrotrophic pathogens.     

Overlapping Functions of the Paralogous Proteins AtPAP2 and AtPAP9 in Arabidopsis thaliana

Arabidopsis thaliana purple acid phosphatase 2 (AtPAP2), which is anchored to the outer membranes of chloroplasts and mitochondria, affects carbon metabolism by modulating the import of some preproteins into chloroplasts and mitochondria. AtPAP9 bears a 72% amino acid sequence identity with AtPAP2, and both proteins carry a hydrophobic motif at their C-termini. Here, we show that AtPAP9 is a tail-anchored protein targeted to the outer membrane of chloroplasts. Yeast two-hybrid and bimolecular fluorescence complementation experiments demonstrated that both AtPAP9 and AtPAP2 bind to a small subunit of rubisco 1B (AtSSU1B) and a number of chloroplast proteins. Chloroplast import assays using [³⁵S]-labeled AtSSU1B showed that like AtPAP2, AtPAP9 also plays a role in AtSSU1B import into chloroplasts. Based on these data, we propose that AtPAP9 and AtPAP2 perform overlapping roles in modulating the import of specific proteins into chloroplasts. Most plant genomes contain only one PAP-like sequence encoding a protein with a hydrophobic motif at the C-terminus. The presence of both AtPAP2 and AtPAP9 in the Arabidopsis genome may have arisen from genome duplication in Brassicaceae. Unlike AtPAP2 overexpression lines, the AtPAP9 overexpression lines did not exhibit early-bolting or high-seed-yield phenotypes. Their differential growth phenotypes could be due to the inability of AtPAP9 to be targeted to mitochondria, as the overexpression of AtPAP2 on mitochondria enhances the capacity of mitochondria to consume reducing equivalents.     

Differential Effects of Indole and Aliphatic Glucosinolates on Lepidopteran Herbivores

Glucosinolates are a diverse group of defensive secondary metabolites that is characteristic of the Brassicales. Arabidopsis thaliana (L.) Heynh. (Brassicaceae) lines with mutations that greatly reduce abundance of indole glucosinolates (cyp79B2 cyp79B3), aliphatic glucosinolates (myb28 myb29), or both (cyp79B2 cyp79B3 myb28 myb29) make it possible to test the in vivo defensive function of these two major glucosinolate classes. In experiments with Lepidoptera that are not crucifer-feeding specialists, aliphatic and indole glucosinolates had an additive effect on Spodoptera exigua (Hübner) (Lepidoptera: Noctuidae) larval growth, whereas Trichoplusia ni (Hübner) (Lepidoptera: Noctuidae) and Manduca sexta (L.) (Lepidoptera: Sphingidae) were affected only by the absence of aliphatic glucosinolates. In the case of two crucifer-feeding specialists, Pieris rapae (L.) (Lepidoptera: Pieridae) and Plutella xylostella (L.) (Lepidoptera: Plutellidae), there were no major changes in larval performance due to decreased aliphatic and/or indole glucosinolate content. Nevertheless, choice tests show that aliphatic and indole glucosinolates act in an additive manner to promote larval feeding of both species and P. rapae oviposition. Together, these results support the hypothesis that a diversity of glucosinolates is required to limit the growth of multiple insect herbivores.     

Characterization of a Crabs Claw Gene in Basal Eudicot Species Epimedium sagittatum (Berberidaceae)

The Crabs Claw (CRC) YABBY gene is required for regulating carpel development in angiosperms and has played an important role in nectary evolution during core eudicot speciation. The function or expression of CRC-like genes has been explored in two basal eudicots, Eschscholzia californica and Aquilegia formosa. To further investigate the function of CRC orthologous genes related to evolution of carpel and nectary development in basal eudicots, a CRC ortholog, EsCRC, was isolated and characterized from Epimedium sagittatum (Sieb. and Zucc.) Maxim. A phylogenetic analysis of EsCRC and previously identified CRC-like genes placed EsCRC within the basal eudicot lineage. Gene expression results suggest that EsCRC is involved in the development of sepals and carpels, but not nectaries. Phenotypic complementation of the Arabidopsis mutant crc-1 was achieved by constitutive expression of EsCRC. In addition, over-expression of EsCRC in Arabidopsis and tobacco gave rise to abaxially curled leaves. Transgenic results together with the gene expression analysis suggest that EsCRC may maintain a conserved function in carpel development and also play a novel role related to sepal formation. Absence of EsCRC and ElCRC expression in nectaries further indicates that nectary development in non-core eudicots is unrelated to expression of CRC-like genes.     

Over-Expression of GmGIa-Regulated Soybean miR172a Confers Early Flowering in Transgenic Arabidopsis thaliana

Flowering is a pivotal event in the life cycle of plants. miR172 has been widely confirmed to play critical roles in flowering time control by regulating its target gene expression in Arabidopsis. However, the role of its counterpart in soybean remains largely unclear. In the present study, we found that the gma-miR172a was regulated by a GIGANTEA ortholog, GmGIa, in soybean through miRNA metabolism. The expression analysis revealed that gma-miR172a has a pattern of diurnal rhythm expression and its abundance increased rapidly as plants grew until the initiation of flowering phase in soybean. One target gene of gma-miR172a, Glyma03g33470, was predicted and verified using a modified RLM 5′-RACE (RNA ligase-mediated rapid amplification of 5′ cDNA ends) assay. Overexpression of gma-miR172a exhibited an early flowering phenotype and the expression of FT, AP1 and LFY were simultaneously increased in gma-miR172a-transgenic Arabidopsis plants, suggesting that the early flowering phenotype was associated with up-regulation of these genes. The overexpression of the gma-miR172a-resistant version of Glyma03g33470 weakened early flowering phenotype in the toe1 mutant of Arabidopsis. Taken together, our results suggested that gma-miR172a played an important role in GmGIa-mediated flowering by repressing Glyma03g33470, which in turn increased the expression of FT, AP1 and LFY to promote flowering in soybean.