LTR_FINDER: an efficient tool for the prediction of full-length LTR retrotransposons

Long terminal repeat retrotransposons (LTR elements) are ubiquitous eukaryotic transposable elements. They play important roles in the evolution of genes and genomes. Ever-growing amount of genomic sequences of many organisms present a great challenge to fast identifying them. That is the first and indispensable step to study their structure, distribution, functions and other biological impacts. However, until today, tools for efficient LTR retrotransposon discovery are very limited. Thus, we developed LTR_FINDER web server. Given DNA sequences, it predicts locations and structure of full-length LTR retrotransposons accurately by considering common structural features. LTR_FINDER is a system capable of scanning large-scale sequences rapidly and the first web server for ab initio LTR retrotransposon finding. We illustrate its usage and performance on the genome of Saccharomyces cerevisiae. The web server is freely accessible at http://tlife.fudan.edu.cn/ltr_finder/.

Genome sequence of the Brown Norway rat yields insights into mammalian evolution

The laboratory rat (Rattus norvegicus) is an indispensable tool in experimental medicine and drug development, having made inestimable contributions to human health. We report here the genome sequence of the Brown Norway (BN) rat strain. The sequence represents a high-quality 'draft' covering over 90% of the genome. The BN rat sequence is the third complete mammalian genome to be deciphered, and three-way comparisons with the human and mouse genomes resolve details of mammalian evolution. This first comprehensive analysis includes genes and proteins and their relation to human disease, repeated sequences, comparative genome-wide studies of mammalian orthologous chromosomal regions and rearrangement breakpoints, reconstruction of ancestral karyotypes and the events leading to existing species, rates of variation, and lineage-specific and lineage-independent evolutionary events such as expansion of gene families, orthology relations and protein evolution.

Trex1 prevents cell-intrinsic initiation of autoimmunity

Detection of nucleic acids and induction of type I interferons (IFNs) are principal elements of antiviral defense but can cause autoimmunity if misregulated. Cytosolic DNA detection activates a potent, cell-intrinsic antiviral response through a poorly defined pathway. In a screen for proteins relevant to this IFN-stimulatory DNA (ISD) response, we identify 3' repair exonuclease 1 (Trex1). Mutations in the human trex1 gene cause Aicardi-Goutieres syndrome (AGS) and chilblain lupus, but the molecular basis of these diseases is unknown. We define Trex1 as an essential negative regulator of the ISD response and delineate the genetic pathway linking Trex1 deficiency to lethal autoimmunity. We show that single-stranded DNA derived from endogenous retroelements accumulates in Trex1-deficient cells, and that Trex1 can metabolize reverse-transcribed DNA. These findings reveal a cell-intrinsic mechanism for initiation of autoimmunity, implicate the ISD pathway as the cause of AGS, and suggest an unanticipated contribution of endogenous retroelements to autoimmunity.

Alu repeats and human genomic diversity

During the past 65 million years, Alu elements have propagated to more than one million copies in primate genomes, which has resulted in the generation of a series of Alu subfamilies of different ages. Alu elements affect the genome in several ways, causing insertion mutations, recombination between elements, gene conversion and alterations in gene expression. Alu-insertion polymorphisms are a boon for the study of human population genetics and primate comparative genomics because they are neutral genetic markers of identical descent with known ancestral states.

Epigenetic inheritance at the agouti locus in the mouse

Epigenetic modifications have effects on phenotype, but they are generally considered to be cleared on passage through the germ line in mammals, so that only genetic traits are inherited. Here we describe the inheritance of an epigenetic modification at the agouti locus in mice. In viable yellow ( A(vy)/a) mice, transcription originating in an intra-cisternal A particle (IAP) retrotransposon inserted upstream of the agouti gene (A) causes ectopic expression of agouti protein, resulting in yellow fur, obesity, diabetes and increased susceptibility to tumours. The pleiotropic effects of ectopic agouti expression are presumably due to effects of the paracrine signal on other tissues. Avy mice display variable expressivity because they are epigenetic mosaics for activity of the retrotransposon: isogenic Avy mice have coats that vary in a continuous spectrum from full yellow, through variegated yellow/agouti, to full agouti (pseudoagouti). The distribution of phenotypes among offspring is related to the phenotype of the dam; when an A(vy) dam has the agouti phenotype, her offspring are more likely to be agouti. We demonstrate here that this maternal epigenetic effect is not the result of a maternally contributed environment. Rather, our data show that it results from incomplete erasure of an epigenetic modification when a silenced Avy allele is passed through the female germ line, with consequent inheritance of the epigenetic modification. Because retrotransposons are abundant in mammalian genomes, this type of inheritance may be common.

Tumour microvesicles contain retrotransposon elements and amplified oncogene sequences

Tumour cells release an abundance of microvesicles containing a selected set of proteins and RNAs. Here, we show that tumour microvesicles also carry DNA, which reflects the genetic status of the tumour, including amplification of the oncogene c-Myc. We also find amplified c-Myc in serum microvesicles from tumour-bearing mice. Further, we find remarkably high levels of retrotransposon RNA transcripts, especially for some human endogenous retroviruses, such as LINE-1 and Alu retrotransposon elements, in tumour microvesicles and these transposable elements could be transferred to normal cells. These findings expand the nucleic acid content of tumour microvesicles to include: elevated levels of specific coding and non-coding RNA and DNA, mutated and amplified oncogene sequences and transposable elements. Thus, tumour microvesicles contain a repertoire of genetic information available for horizontal gene transfer and potential use as blood biomarkers for cancer.

Meiotic catastrophe and retrotransposon reactivation in male germ cells lacking Dnmt3L

Mammalian genomes employ heritable cytosine methylation in the long-term silencing of retrotransposons and genes subject to genomic imprinting and X chromosome inactivation. Little is known of the mechanisms that direct cytosine methylation to specific sequences. Here we show that DNA methyltransferase 3-like (Dnmt3L (ref. 1)) is expressed in testes during a brief perinatal period in the non-dividing precursors of spermatogonial stem cells at a stage where retrotransposons undergo de novo methylation. Deletion of the Dnmt3L gene prevented the de novo methylation of both long-terminal-repeat (LTR) and non-LTR retrotransposons, which were transcribed at high levels in spermatogonia and spermatocytes. Loss of Dnmt3L from early germ cells also caused meiotic failure in spermatocytes, which do not express Dnmt3L. Whereas dispersed repeated sequences were demethylated in mutant germ cells, tandem repeats in pericentric regions were methylated normally. This result indicates that the Dnmt3L protein might have a function in the de novo methylation of dispersed repeated sequences in a premeiotic genome scanning process that occurs in male germ cells at about the time of birth.

Human L1 retrotransposon encodes a conserved endonuclease required for retrotransposition

Human L1 elements are highly abundant poly(A) (non-LTR) retrotransposons whose second open reading frame (ORF2) encodes a reverse transcriptase (RT). We have identified an endonuclease (EN) domain at the L1 ORF2 N-terminus that is highly conserved among poly(A) retrotransposons and resembles the apurinic/apyrimidinic (AP) endonucleases. Purified L1 EN protein (L1 ENp) makes 5'-PO4, 3'-OH nicks in supercoiled plasmids, shows no preference for AP sites, and preferentially cleaves sequences resembling L1 in vivo target sequences. Mutations in conserved amino acid residues of L1 EN abolish its nicking activity and eliminate L1 retrotransposition. We propose that L1 EN cleaves the target site for L1 insertion and primes reverse transcription.

High frequency retrotransposition in cultured mammalian cells

We previously isolated two human L1 elements (L1.2 and LRE2) as the progenitors of disease-producing insertions. Here, we show these elements can actively retrotranspose in cultured mammalian cells. When stably expressed from an episome in HeLa cells, both elements retrotransposed into a variety of chromosomal locations at a high frequency. The retrotransposed products resembled endogenous L1 insertions, since they were variably 5' truncated, ended in poly(A) tracts, and were flanked by target-site duplications or short deletions. Point mutations in conserved domains of the L1.2-encoded proteins reduced retrotransposition by 100- to 1000-fold. Remarkably, L1.2 also retrotransposed in a mouse cell line, suggesting a potential role for L1-based vectors in random insertional mutagenesis.

Somatic mosaicism in neuronal precursor cells mediated by L1 retrotransposition

Revealing the mechanisms for neuronal somatic diversification remains a central challenge for understanding individual differences in brain organization and function. Here we show that an engineered human LINE-1 (for long interspersed nuclear element-1; also known as L1) element can retrotranspose in neuronal precursors derived from rat hippocampus neural stem cells. The resulting retrotransposition events can alter the expression of neuronal genes, which, in turn, can influence neuronal cell fate in vitro. We further show that retrotransposition of a human L1 in transgenic mice results in neuronal somatic mosaicism. The molecular mechanism of action is probably mediated through Sox2, because a decrease in Sox2 expression during the early stages of neuronal differentiation is correlated with increases in both L1 transcription and retrotransposition. Our data therefore indicate that neuronal genomes might not be static, but some might be mosaic because of de novo L1 retrotransposition events.

Two classes of short interfering RNA in RNA silencing

RNA silencing is a eukaryotic genome defence system that involves processing of double-stranded RNA (dsRNA) into 21-26 nt, short interfering RNA (siRNA). The siRNA mediates suppression of genes corresponding to the dsRNA through targeted RNA degradation. In some plant systems there are additional silencing processes, involving systemic spread of silencing and RNA-directed methylation/transcriptional suppression of homologous genomic DNA. We show here that siRNAs produced in plants from a green fluorescent protein (GFP) transgene are in short (21-22 nt) and long (24-26 nt) size classes, whereas those from endogenous retroelements are only in the long class. Viral suppressors of RNA silencing and mutations in Arabidopsis indicate that these classes of siRNA have different roles. The long siRNA is dispensable for sequence-specific mRNA degradation, but correlates with systemic silencing and methylation of homologous DNA. Conversely, the short siRNA class correlates with mRNA degradation but not with systemic signalling or methylation. These findings reveal an unexpected level of complexity in the RNA silencing pathway in plants that may also apply in animals.

Annotation, submission and screening of repetitive elements in Repbase: RepbaseSubmitter and Censor

Background

Repbase is a reference database of eukaryotic repetitive DNA, which includes prototypic sequences of repeats and basic information described in annotations. Updating and maintenance of the database requires specialized tools, which we have created and made available for use with Repbase, and which may be useful as a template for other curated databases.

Results

We describe the software tools RepbaseSubmitter and Censor, which are designed to facilitate updating and screening the content of Repbase. RepbaseSubmitter is a java-based interface for formatting and annotating Repbase entries. It eliminates many common formatting errors, and automates actions such as calculation of sequence lengths and composition, thus facilitating curation of Repbase sequences. In addition, it has several features for predicting protein coding regions in sequences; searching and including Pubmed references in Repbase entries; and searching the NCBI taxonomy database for correct inclusion of species information and taxonomic position. Censor is a tool to rapidly identify repetitive elements by comparison to known repeats. It uses WU-BLAST for speed and sensitivity, and can conduct DNA-DNA, DNA-protein, or translated DNA-translated DNA searches of genomic sequence. Defragmented output includes a map of repeats present in the query sequence, with the options to report masked query sequence(s), repeat sequences found in the query, and alignments.

Conclusion

Censor and RepbaseSubmitter are available as both web-based services and downloadable versions. They can be found at (RepbaseSubmitter) and (Censor).

The paleontology of intergene retrotransposons of maize

Retrotransposons, transposable elements related to animal retroviruses, are found in all eukaryotes investigated and make up the majority of many plant genomes. Their ubiquity points to their importance, especially in their contribution to the large-scale structure of complex genomes. The nature and frequency of retro-element appearance, activation and amplification are poorly understood in all higher eukaryotes. Here we employ a novel approach to determine the insertion dates for 17 of 23 retrotransposons found near the maize adh1 gene, and two others from unlinked sites in the maize genome, by comparison of long terminal repeat (LTR) divergences with the sequence divergence between adh1 in maize and sorghum. All retrotransposons examined have inserted within the last six million years, most in the last three million years. The structure of the adh1 region appears to be standard relative to the other gene-containing regions of the maize genome, thus suggesting that retrotransposon insertions have increased the size of the maize genome from approximately 1200 Mb to 2400 Mb in the last three million years. Furthermore, the results indicate an increased mutation rate in retrotransposons compared with genes.

L1 retrotransposition in human neural progenitor cells

Long interspersed element 1 (LINE-1 or L1) retrotransposons have markedly affected the human genome. L1s must retrotranspose in the germ line or during early development to ensure their evolutionary success, yet the extent to which this process affects somatic cells is poorly understood. We previously demonstrated that engineered human L1s can retrotranspose in adult rat hippocampus progenitor cells in vitro and in the mouse brain in vivo. Here we demonstrate that neural progenitor cells isolated from human fetal brain and derived from human embryonic stem cells support the retrotransposition of engineered human L1s in vitro. Furthermore, we developed a quantitative multiplex polymerase chain reaction that detected an increase in the copy number of endogenous L1s in the hippocampus, and in several regions of adult human brains, when compared to the copy number of endogenous L1s in heart or liver genomic DNAs from the same donor. These data suggest that de novo L1 retrotransposition events may occur in the human brain and, in principle, have the potential to contribute to individual somatic mosaicism.

The genome of a songbird

The zebra finch is an important model organism in several fields with unique relevance to human neuroscience. Like other songbirds, the zebra finch communicates through learned vocalizations, an ability otherwise documented only in humans and a few other animals and lacking in the chicken-the only bird with a sequenced genome until now. Here we present a structural, functional and comparative analysis of the genome sequence of the zebra finch (Taeniopygia guttata), which is a songbird belonging to the large avian order Passeriformes. We find that the overall structures of the genomes are similar in zebra finch and chicken, but they differ in many intrachromosomal rearrangements, lineage-specific gene family expansions, the number of long-terminal-repeat-based retrotransposons, and mechanisms of sex chromosome dosage compensation. We show that song behaviour engages gene regulatory networks in the zebra finch brain, altering the expression of long non-coding RNAs, microRNAs, transcription factors and their targets. We also show evidence for rapid molecular evolution in the songbird lineage of genes that are regulated during song experience. These results indicate an active involvement of the genome in neural processes underlying vocal communication and identify potential genetic substrates for the evolution and regulation of this behaviour.

The wheat and barley vernalization gene VRN3 is an orthologue of FT

Winter wheat and barley varieties require an extended exposure to low temperatures to accelerate flowering (vernalization), whereas spring varieties do not have this requirement. In this study, we show that in these species, the vernalization gene VRN3 is linked completely to a gene similar to Arabidopsis FLOWERING LOCUS T (FT). FT induction in the leaves results in a transmissible signal that promotes flowering. Transcript levels of the barley and wheat orthologues, designated as HvFT and TaFT, respectively, are significantly higher in plants homozygous for the dominant Vrn3 alleles (early flowering) than in plants homozygous for the recessive vrn3 alleles (late flowering). In wheat, the dominant Vrn3 allele is associated with the insertion of a retroelement in the TaFT promoter, whereas in barley, mutations in the HvFT first intron differentiate plants with dominant and recessive VRN3 alleles. Winter wheat plants transformed with the TaFT allele carrying the promoter retroelement insertion flowered significantly earlier than nontransgenic plants, supporting the identity between TaFT and VRN-B3. Statistical analyses of flowering times confirmed the presence of significant interactions between vernalization and FT allelic classes in both wheat and barley (P < 0.0001). These interactions were supported further by the observed up-regulation of HvFT transcript levels by vernalization in barley winter plants (P = 0.002). These results confirmed that the wheat and barley FT genes are responsible for natural allelic variation in vernalization requirement, providing additional sources of adaptive diversity to these economically important crops.

Human LINE retrotransposons generate processed pseudogenes

Long interspersed elements (LINEs) are endogenous mobile genetic elements that have dispersed and accumulated in the genomes of higher eukaryotes via germline transposition, with up to 100,000 copies in mammalian genomes. In humans, LINEs are the major source of insertional mutagenesis, being involved in both germinal and somatic mutant phenotypes. Here we show that the human LINE retrotransposons, which transpose through the reverse transcription of their own transcript, can also mobilize transcribed DNA not associated with a LINE sequence by a process involving the diversion of the LINE enzymatic machinery by the corresponding mRNA transcripts. This results in the 'retroposition' of the transcribed gene and the formation of new copies that disclose features characteristic of the widespread and naturally occurring processed pseudogenes: loss of intron and promoter, acquisition of a poly(A) 3' end and presence of target-site duplications of varying length. We further show-by introducing deletions within either coding sequence of the human LINE-that both ORFs are necessary for the formation of the processed pseudogenes, and that retroviral-like elements are not able to produce similar structures in the same assay. Our results strengthen the unique versatility of LINEs as genome modellers.

The regulated retrotransposon transcriptome of mammalian cells

Although repetitive elements pervade mammalian genomes, their overall contribution to transcriptional activity is poorly defined. Here, as part of the FANTOM4 project, we report that 6-30% of cap-selected mouse and human RNA transcripts initiate within repetitive elements. Analysis of approximately 250,000 retrotransposon-derived transcription start sites shows that the associated transcripts are generally tissue specific, coincide with gene-dense regions and form pronounced clusters when aligned to full-length retrotransposon sequences. Retrotransposons located immediately 5' of protein-coding loci frequently function as alternative promoters and/or express noncoding RNAs. More than a quarter of RefSeqs possess a retrotransposon in their 3' UTR, with strong evidence for the reduced expression of these transcripts relative to retrotransposon-free transcripts. Finally, a genome-wide screen identifies 23,000 candidate regulatory regions derived from retrotransposons, in addition to more than 2,000 examples of bidirectional transcription. We conclude that retrotransposon transcription has a key influence upon the transcriptional output of the mammalian genome.

Plant transposable elements: where genetics meets genomics

Transposable elements are the single largest component of the genetic material of most eukaryotes. The recent availability of large quantities of genomic sequence has led to a shift from the genetic characterization of single elements to genome-wide analysis of enormous transposable-element populations. Nowhere is this shift more evident than in plants, in which transposable elements were first discovered and where they are still actively reshaping genomes.

Genome sequence of the lignocellulose degrading fungus Phanerochaete chrysosporium strain RP78

White rot fungi efficiently degrade lignin, a complex aromatic polymer in wood that is among the most abundant natural materials on earth. These fungi use extracellular oxidative enzymes that are also able to transform related aromatic compounds found in explosive contaminants, pesticides and toxic waste. We have sequenced the 30-million base-pair genome of Phanerochaete chrysosporium strain RP78 using a whole genome shotgun approach. The P. chrysosporium genome reveals an impressive array of genes encoding secreted oxidases, peroxidases and hydrolytic enzymes that cooperate in wood decay. Analysis of the genome data will enhance our understanding of lignocellulose degradation, a pivotal process in the global carbon cycle, and provide a framework for further development of bioprocesses for biomass utilization, organopollutant degradation and fiber bleaching. This genome provides a high quality draft sequence of a basidiomycete, a major fungal phylum that includes important plant and animal pathogens.

Epigenetics and aging

Over the past decade, a growing number of studies have revealed that progressive changes to epigenetic information accompany aging in both dividing and nondividing cells. Functional studies in model organisms and humans indicate that epigenetic changes have a huge influence on the aging process. These epigenetic changes occur at various levels, including reduced bulk levels of the core histones, altered patterns of histone posttranslational modifications and DNA methylation, replacement of canonical histones with histone variants, and altered noncoding RNA expression, during both organismal aging and replicative senescence. The end result of epigenetic changes during aging is altered local accessibility to the genetic material, leading to aberrant gene expression, reactivation of transposable elements, and genomic instability. Strikingly, certain types of epigenetic information can function in a transgenerational manner to influence the life span of the offspring. Several important conclusions emerge from these studies: rather than being genetically predetermined, our life span is largely epigenetically determined; diet and other environmental influences can influence our life span by changing the epigenetic information; and inhibitors of epigenetic enzymes can influence life span of model organisms. These new findings provide better understanding of the mechanisms involved in aging. Given the reversible nature of epigenetic information, these studies highlight exciting avenues for therapeutic intervention in aging and age-associated diseases, including cancer.

Cyclophilin A retrotransposition into TRIM5 explains owl monkey resistance to HIV-1

In Old World primates, TRIM5-alpha confers a potent block to human immunodeficiency virus type 1 (HIV-1) infection that acts after virus entry into cells. Cyclophilin A (CypA) binding to viral capsid protects HIV-1 from a similar activity in human cells. Among New World primates, only owl monkeys exhibit post-entry restriction of HIV-1 (ref. 1). Paradoxically, the barrier to HIV-1 in owl monkey cells is released by capsid mutants or drugs that disrupt capsid interaction with CypA. Here we show that knockdown of owl monkey CypA by RNA interference (RNAi) correlates with suppression of anti-HIV-1 activity. However, reintroduction of CypA protein to RNAi-treated cells did not restore antiviral activity. A search for additional RNAi targets unearthed TRIMCyp, an RNAi-responsive messenger RNA encoding a TRIM5-CypA fusion protein. TRIMCyp accounts for post-entry restriction of HIV-1 in owl monkeys and blocks HIV-1 infection when transferred to otherwise infectable human or rat cells. It seems that TRIMCyp arose after the divergence of New and Old World primates when a LINE-1 retrotransposon catalysed the insertion of a CypA complementary DNA into the TRIM5 locus. This is the first vertebrate example of a chimaeric gene generated by this mechanism of exon shuffling.

An endogenous small interfering RNA pathway in Drosophila

Drosophila endogenous small RNAs are categorized according to their mechanisms of biogenesis and the Argonaute protein to which they bind. MicroRNAs are a class of ubiquitously expressed RNAs of approximately 22 nucleotides in length, which arise from structured precursors through the action of Drosha-Pasha and Dicer-1-Loquacious complexes. These join Argonaute-1 to regulate gene expression. A second endogenous small RNA class, the Piwi-interacting RNAs, bind Piwi proteins and suppress transposons. Piwi-interacting RNAs are restricted to the gonad, and at least a subset of these arises by Piwi-catalysed cleavage of single-stranded RNAs. Here we show that Drosophila generates a third small RNA class, endogenous small interfering RNAs, in both gonadal and somatic tissues. Production of these RNAs requires Dicer-2, but a subset depends preferentially on Loquacious rather than the canonical Dicer-2 partner, R2D2 (ref. 14). Endogenous small interfering RNAs arise both from convergent transcription units and from structured genomic loci in a tissue-specific fashion. They predominantly join Argonaute-2 and have the capacity, as a class, to target both protein-coding genes and mobile elements. These observations expand the repertoire of small RNAs in Drosophila, adding a class that blurs distinctions based on known biogenesis mechanisms and functional roles.

Somatic retrotransposition alters the genetic landscape of the human brain

Retrotransposons are mobile genetic elements that use a germline 'copy-and-paste' mechanism to spread throughout metazoan genomes. At least 50 per cent of the human genome is derived from retrotransposons, with three active families (L1, Alu and SVA) associated with insertional mutagenesis and disease. Epigenetic and post-transcriptional suppression block retrotransposition in somatic cells, excluding early embryo development and some malignancies. Recent reports of L1 expression and copy number variation in the human brain suggest that L1 mobilization may also occur during later development. However, the corresponding integration sites have not been mapped. Here we apply a high-throughput method to identify numerous L1, Alu and SVA germline mutations, as well as 7,743 putative somatic L1 insertions, in the hippocampus and caudate nucleus of three individuals. Surprisingly, we also found 13,692 somatic Alu insertions and 1,350 SVA insertions. Our results demonstrate that retrotransposons mobilize to protein-coding genes differentially expressed and active in the brain. Thus, somatic genome mosaicism driven by retrotransposition may reshape the genetic circuitry that underpins normal and abnormal neurobiological processes.

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.