Recent studies have demonstrated the important role of plant microRNAs (miRNAs) under nutrient deficiencies. In this study, deep sequencing of Arabidopsis (Arabidopsis thaliana) small RNAs was conducted to reveal miRNAs and other small RNAs that were differentially expressed in response to phosphate (Pi) deficiency. About 3.5 million sequence reads corresponding to 0.6 to 1.2 million unique sequence tags from each Pi-sufficient or Pi-deficient root or shoot sample were mapped to the Arabidopsis genome. We showed that upon Pi deprivation, the expression of miR156, miR399, miR778, miR827, and miR2111 was induced, whereas the expression of miR169, miR395, and miR398 was repressed. We found cross talk coordinated by these miRNAs under different nutrient deficiencies. In addition to miRNAs, we identified one Pi starvation-induced DICER-LIKE1-dependent small RNA derived from the long terminal repeat of a retrotransposon and a group of 19-nucleotide small RNAs corresponding to the 5′ end of tRNA and expressed at a high level in Pi-starved roots. Importantly, we observed an increased abundance of TAS4-derived trans-acting small interfering RNAs (ta-siRNAs) in Pi-deficient shoots and uncovered an autoregulatory mechanism of PAP1/MYB75 via miR828 and TAS4-siR81(−) that regulates the biosynthesis of anthocyanin. This finding sheds light on the regulatory network between miRNA/ta-siRNA and its target gene. Of note, a substantial amount of miR399* accumulated under Pi deficiency. Like miR399, miR399* can move across the graft junction, implying a potential biological role for miR399*. This study represents a comprehensive expression profiling of Pi-responsive small RNAs and advances our understanding of the regulation of Pi homeostasis mediated by small RNAs.

Gene expression is regulated both by cis elements, which are DNA segments closely linked to the genes they regulate, and by trans factors, which are usually proteins capable of diffusing to unlinked genes. Understanding the patterns and sources of regulatory variation is crucial for understanding phenotypic and genome evolution. Here, we measure genome-wide allele-specific expression by deep sequencing to investigate the patterns of cis and trans expression variation between two strains of Saccharomyces cerevisiae. We propose a statistical modeling framework based on the binomial distribution that simultaneously addresses normalization of read counts derived from different parents and estimating the cis and trans expression variation parameters. We find that expression polymorphism in yeast is common for both cis and trans, though trans variation is more common. Constraint in expression evolution is correlated with other hallmarks of constraint, including gene essentiality, number of protein interaction partners, and constraint in amino acid substitution, indicating that both cis and trans polymorphism are clearly under purifying selection, though trans variation appears to be more sensitive to selective constraint. Comparing interspecific expression divergence between S. cerevisiae and S. paradoxus to our intraspecific variation suggests a significant departure from a neutral model of molecular evolution. A further examination of correlation between polymorphism and divergence within each category suggests that cis divergence is more frequently mediated by positive Darwinian selection than is trans divergence.

Gene expression evolution can be caused by changes in cis- or trans-regulatory elements or both. As cis and trans regulation operate through different molecular mechanisms, cis and trans mutations may show different inheritance patterns and may be subjected to different selective constraints. To investigate these issues, we obtained and analyzed gene expression data from two Saccharomyces cerevisiae strains and their hybrid, using high-throughput sequencing. Our data indicate that compared with other types of genes, those with antagonistic cis–trans interactions are more likely to exhibit over- or underdominant inheritance of expression level. Moreover, in accordance with previous studies, genes with trans variants tend to have a dominant inheritance pattern, whereas cis variants are enriched for additive inheritance. In addition, cis regulatory differences contribute more to expression differences between species than within species, whereas trans regulatory differences show a stronger association between divergence and polymorphism. Our data indicate that in the trans component of gene expression differences genes subjected to weaker selective constraints tend to have an excess of polymorphism over divergence compared with those subjected to stronger selective constraints. In contrast, in the cis component, this difference between genes under stronger and weaker selective constraint is mostly absent. To explain these observations, we propose that purifying selection more strongly shapes trans changes than cis changes and that positive selection may have significantly contributed to cis regulatory divergence.

Animal microRNAs (miRNAs) are short RNAs that function as posttranscriptional regulators of gene expression by binding to the target mRNAs. Noting that some miRNAs are highly conserved in evolution, we explored the possibility of evolutionary conservation of their targets. We identified human orthologues of experimentally verified let-7 miRNA target genes in Caenorhabditis elegans and used the luciferase reporter system to examine whether these human genes are still the targets of let-7 miRNA. We found that in some cases, the miRNA–target relationship has indeed been conserved in human. Interestingly, human TRIM71, an orthologue of C. elegans let-7–target lin-41 gene, can be repressed by hsa-let-7a and hsa-let-7c. This repression was abolished when both predicted let-7 target sites of TRIM71 were mutated. Moreover, the zebrafish lin-41 orthologue was also repressed by let-7 to a similar degree as was TRIM71. When the expression of zebrafish lin-41 orthologue was silenced by microinjection of RNA interference or morpholino into zebrafish zygotes, retarded embryonic development was observed, providing direct evidence for an essential role of lin-41 in zebrafish development. Taken together, our results suggest that the regulation of TRIM71 expression by let-7 has been evolutionarily conserved and that TRIM71 likely plays an important role in development.

Phytoplasmas are bacterial phytopathogens that release virulence effectors into sieve cells and act systemically to affect the physiological and morphological state of host plants to promote successful pathogenesis. We show here that transgenic Nicotiana benthamiana lines expressing the secreted effector SAP11 from Candidatus Phytoplasma mali exhibit an altered aroma phenotype. This phenomenon is correlated with defects in the development of glandular trichomes and the biosynthesis of 3-isobutyl-2-methoxypyrazine (IBMP). IBMP is a volatile organic compound (VOC) synthesized by an O-methyltransferase, via a methylation step, from a non-volatile precursor, 3-isobutyl-2-hydroxypyrazine (IBHP). Based on comparative and functional genomics analyses, NbOMT1, which encodes an O-methyltransferase, was found to be highly suppressed in SAP11-transgenic plants. We further silenced NbOMT1 through virus-induced gene silencing and demonstrated that this enzyme influenced the accumulation of IBMP in N. benthamiana. In vitro biochemical analyses also showed that NbOMT1 can catalyse IBHP O-methylation in the presence of S-adenosyl-L-methionine. Our study suggests that the phytoplasma effector SAP11 has the ability to modulate host VOC emissions. In addition, we also demonstrated that SAP11 destabilized TCP transcription factors and suppressed jasmonic acid responses in N. benthamiana. These findings provide valuable insights into understanding how phytoplasma effectors influence plant volatiles.

Several research advances have indicated an important role of transfer RNA (tRNA)-derived small RNAs in modulating developmental processes or stress responses. Recently, from the deep sequencing of small RNAs in Arabidopsis (Arabidopsis thaliana), we identified a new class of 19-nucleotide (nt) small RNAs corresponding to the 5′ end of tRNA accumulated at high levels in phosphate-starved roots. In two very recent studies, 19-nt tRNA fragments were also observed in human cells, suggesting their widespread nature. In our study, tRNA halves cleaved at the anticodon loop, the most common tRNA fragments found, were predominant in roots. These results showed a spatial and temporal expression pattern of small RNAs derived from specific cleavage of tRNA molecules. Although the function of these tRNA-derived small RNAs under phosphate deficiency remains unknown, their diversity, biogenesis and potential function are henceforth summarized and discussed. Certainly, they will emerge as a novel class of regulatory small RNAs.