Li-Ching Hsieh

Uncovering Small RNA-Mediated Responses to Phosphate Deficiency in Arabidopsis by Deep Sequencing

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.

Natural selection on cis and trans regulation in yeasts

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.

Minimal Model for Genome Evolution and Growth

Textual analysis of typical microbial genomes reveals that they have the statistical characteristics of a DNA sequence of a much shorter length. This peculiar property supports an evolutionary model in which a genome evolves by random mutation but primarily grows by random segmental duplication. That genomes grew mostly by duplication is consistent with the observation that repeat sequences in all genomes are widespread and intragenomic and intergenomic homologous genes are preponderant across all life forms.

Inheritance of Gene Expression Level and Selective Constraints on Trans- and Cis-Regulatory Changes in Yeast

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.

Abundance of tRNA-derived small RNAs in phosphate-starved Arabidopsis roots

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.

Divergence and Shannon Information in Genomes

Shannon information (SI) and its special case, divergence, are defined for a DNA sequence in terms of probabilities of chemical words in the sequence and are computed for a set of complete genomes highly diverse in length and composition. We find the following: SI (but not divergence) is inversely proportional to sequence length for a random sequence but is length independent for genomes; the genomic SI is always greater and, for shorter words and longer sequences, hundreds to thousands times greater than the SI in a random sequence whose length and composition match those of the genome; genomic SIs appear to have word-length dependent universal values. The universality is inferred to be an evolution footprint of a universal mode for genome growth.

SHANNON INFORMATION IN COMPLETE GENOMES

Shannon information in the genomes of all completely sequenced prokaryotes and eukaryotes are measured in word lengths of two to ten letters. It is found that in a scale-dependent way, the Shannon information in complete genomes are much greater than that in matching random sequences - thousands of times greater in the case of short words. Furthermore, with the exception of the 14 chromosomes of Plasmodium falciparum, the Shannon information in all available complete genomes belong to a universality class given by an extremely simple formula. The data are consistent with a model for genome growth composed of two main ingredients: random segmental duplications that increase the Shannon information in a scale-independent way, and random point mutations that preferentially reduces the larger-scale Shannon information. The inference drawn from the present study is that the large-scale and coarse-grained growth of genomes was selectively neutral and this suggests an independent corroboration of Kimuras neutral theory of evolution.

Short segmental duplication: parsimony in growth of microbial genomes

We compare the distributions of occurrence frequencies of oligonucleotides two to ten bases long (2 to 10-mers) in microbial complete genomes with corresponding distributions obtained from random sequences and find that the genomic distributions are uniformly many times wider in a universal manner, that is, the same for all microbial complete genomes. The difference increases with decreasing word length, with the genomic spectral width about 40 times wider for 2-mers. We show that the observed genomic properties are characteristic of sequences generated in a simple growth model, where a very short initial random sequence (less than 1 kb) grows mainly by maximally stochastic duplication of short segments (of about 25 b). We discuss a number issues related to the findings and the model, including the proposition that life began in an RNA world before the birth of proteins.

Phytoplasma SAP11 alters 3-isobutyl-2-methoxypyrazine biosynthesis in Nicotiana benthamiana by suppressing NbOMT1

Highlight Phytoplasma effector SAP11 modulates plant volatile organic compound emissions by suppressing the expression of NbOMT1, which encodes an O-methyltransferase required for the biosynthesis of 3-isobutyl-2-methoxypyrazine.

Parity-dependent hairpin configurations of repetitive DNA sequence promote slippage associated with DNA expansion

Significance We found that TGGAA DNA repeats, which are involved in the neurological disease spinocerebellar ataxia 31, are capable of assuming two different hairpin structures depending on repeat number parity. We determined the interconversion kinetics by single-molecule spectroscopy and probed the interconversion mechanism through elucidation of the TGGAA repeat stem structure. Our results suggest that the two hairpin structures interconvert through motion slippage, and the process can be explained by the overall stem stability and local destabilization of the kinked GGA motif. Divalent cations and stem length affected the equilibrium and kinetics of slippage. Our findings suggest a mechanism by which a binary dynamic property of DNA repeats may affect repeat expansion and may be applicable to other repetitive DNA systems. Repetitive DNA sequences are ubiquitous in life, and changes in the number of repeats often have various physiological and pathological implications. DNA repeats are capable of interchanging between different noncanonical and canonical conformations in a dynamic fashion, causing configurational slippage that often leads to repeat expansion associated with neurological diseases. In this report, we used single-molecule spectroscopy together with biophysical analyses to demonstrate the parity-dependent hairpin structural polymorphism of TGGAA repeat DNA. We found that the DNA adopted two configurations depending on the repeat number parity (even or odd). Transitions between these two configurations were also observed for longer repeats. In addition, the ability to modulate this transition was found to be enhanced by divalent ions. Based on the atomic structure, we propose a local seeding model where the kinked GGA motifs in the stem region of TGGAA repeat DNA act as hot spots to facilitate the transition between the two configurations, which may give rise to disease-associated repeat expansion.