text:
- © 2013 American Society of Plant Biologists. All Rights Reserved.
The Arabidopsis (Arabidopsis thaliana) genome is the most well-annotated plant genome. However, transcriptome sequencing in Arabidopsis continues to suggest the presence of polyadenylated (polyA) transcripts originating from presumed intergenic regions. It is not clear whether these transcripts represent novel noncoding or protein-coding genes. To understand the nature of intergenic polyA transcription, we first assessed its abundance using multiple messenger RNA sequencing data sets. We found 6,545 intergenic transcribed fragments (ITFs) occupying 3.6% of Arabidopsis intergenic space. In contrast to transcribed fragments that map to protein-coding and RNA genes, most ITFs are significantly shorter, are expressed at significantly lower levels, and tend to be more data set specific. A surprisingly large number of ITFs (32.1%) may be protein coding based on evidence of translation. However, our results indicate that these "translated" ITFs tend to be close to and are likely associated with known genes. To investigate if ITFs are under selection and are functional, we assessed ITF conservation through cross-species as well as within-species comparisons. Our analysis reveals that 237 ITFs, including 49 with translation evidence, are under strong selective constraint and relatively distant from annotated features. These ITFs are likely parts of novel genes. However, the selective pressure imposed on most ITFs is similar to that of randomly selected, untranscribed intergenic sequences. Our findings indicate that despite the prevalence of ITFs, apart from the possibility of genomic contamination, many may be background or noisy transcripts derived from "junk" DNA, whose production may be inherent to the process of transcription and which, on rare occasions, may act as catalysts for the creation of novel genes.
The advent of tiling arrays and high-throughput sequencing has led to the discovery of a complex transcriptional landscape in eukaryotic genomes. Studies in yeast (Saccharomyces cerevisiae; David et al., 2006), animals (Bertone et al., 2004; Carninci et al., 2005), and plants (Yamada et al., 2003; Li et al., 2007; Matsui et al., 2008) have revealed the presence of a large number of unannotated, novel transcripts. These novel transcripts may represent alternatively spliced forms of known genes (Filichkin et al., 2010), products of antisense (Yamada et al., 2003) or bidirectional transcription (Xu et al., 2009), retained introns (Ner-Gaon et al., 2004; Filichkin et al., 2010), transcript fusions (Ruan et al., 2007), or intergenic transcriptional units (referred to hereafter as intergenic transcribed fragments [ITFs]). Among these novel transcripts, ITFs are unique in that they do not overlap with known genomic features and may represent novel genic sequences. The prevalence of intergenic transcription raises the possibility that there are many more functional genes yet to be discovered. However, there are two outstanding questions regarding ITFs. First, it is not clear what proportion of ITFs code for proteins. Second, whether most ITFs are functional is under debate (Mattick, 2009; Ponting and Belgard, 2010).
After ITFs are identified with whole-genome tiling arrays or high-throughput sequencing, computational methods are used to determine if they display characteristics of noncoding RNA (ncRNA; Fahlgren et al., 2007; Li et al., 2007; Gregory et al., 2008). These methods rely on secondary structure prediction, similarity to known ncRNAs, and conservation between species. The protein-coding potential of ITFs, on the other hand, is determined based on ab initio gene prediction, open reading frame (ORF) length, evolutionary conservation measures, pairwise alignment scores, predicted secondary structure, and entropy (Nekrutenko et al., 2002; Liu et al., 2006; Dinger et al., 2008). For example, in a global gene expression study in Arabidopsis (Arabidopsis thaliana), a 50-amino acid length threshold was used to define potential protein-coding intergenic transcripts (Stolc et al., 2005). Similarly, the Functional Annotation of the Mammalian Genome consortium defined putative protein-coding mRNAs using an ORF length cutoff of 300 nucleotides (Okazaki et al., 2002). Reliance on length cutoffs can result in longer random ORFs being falsely annotated as protein coding and can also lead to the exclusion of true small ORFs such as those that have been identified in yeast, humans, and Arabidopsis (Basrai et al., 1997; Hanada et al., 2007; Pruitt et al., 2007). Proteomics and polyribosome immunoprecipitation (Zanetti et al., 2005; Sparkes et al., 2006) allow more direct identification of potentially protein-coding ITFs than computational approaches. Currently, there has yet to be a systematic assessment of ITF protein-coding potential based on a combination of computational and experimental approaches.
In addition to the question of whether ITFs code for proteins, the functional relevance of intergenic transcription is not well understood. One hypothesis is that most transcripts simply represent transcriptional noise. For example, based on the genome-wide distribution of RNA polymerase II and TATA-binding protein in yeast, approximately 90% of RNA polymerase II transcriptional initiation events were estimated to be the result of low polymerase fidelity and may represent transcriptional noise (Struhl, 2007). Consistent with the "noise" hypothesis, several studies have shown that ITFs tend to have significantly higher evolutionary rates than known genes. For example, the Encyclopedia of DNA Elements (ENCODE) consortium found that 93% of the unannotated transcribed regions in the human genome show no clear evidence of evolutionary constraint (Birney et al., 2007). The alternative hypothesis is that most ITFs are functional (Dinger et al., 2009). Differential expression, alternative splicing, and/or association with chromatin modification marks have been cited as evidence for ITF functionality (Guttman et al., 2009; Hiller et al., 2009). In addition, the functions of a growing number of novel transcripts have been experimentally determined. Examples include Xist, RepA, Air, and Hotair, which regulate the recruitment of Polycomb proteins onto DNA (Mercer et al., 2009), as well as a recently discovered long noncoding RNA called COLDAIR shown to be important in regulating vernalization responses in Arabidopsis (Heo and Sung, 2011). Based on these studies, it is clear that some ITFs are functional. The main question concerns the abundance of functional ITFs relative to those derived from noisy transcription.
To date, most studies of intergenic transcription have focused on the presumably noncoding fraction of the transcriptome. In addition, currently, there is no published study assessing the evolutionary significance of plant intergenic transcription. In this study, we focused on intergenic polyadenylated (polyA) RNA transcripts to gain more insight into the nature of plant intergenic transcription by RNA polymerase II. We first analyzed eight different Arabidopsis mRNA sequencing (mRNA-seq) data sets from this study and two other sources (Filichkin et al., 2010; Jiao and Meyerowitz, 2010) to determine the extent of intergenic polyA transcription. We then investigated whether ITFs are likely protein coding using (1) ribosome immunoprecipitation data generated in this study as well as public data sets (Jiao and Meyerowitz, 2010), (2) proteomics data (Baerenfaller et al., 2008; Castellana et al., 2008), and (3) fusion protein expression studies on selected targets. Finally, making use of the polymorphism data from 80 different Arabidopsis accessions (Cao et al., 2011) and protein-coding genes and genome sequences of other plants, we explored whether ITFs, especially those that may code for proteins, are likely functional based on within- and cross-species conservation.
Defining Transcribed Regions in the Arabidopsis Genome
To explore the functional significance of intergenic transcription further, a rigorous definition of transcribed regions within the Arabidopsis genome is necessary. To this end, we analyzed mRNA-seq data from three different sources: (1) 7-d-old seedlings generated in this study, (2) whole flower (Jiao and Meyerowitz, 2010), and (3) 12-d-old seedlings grown under six environmental conditions (Filichkin et al., 2010; Supplemental Table S1). We assembled transcript fragments (TxFrags) using two approaches (see "Materials and Methods"). In the first approach, contiguous regions in Arabidopsis occupied by mapped mRNA-seq reads were defined as expressed (Set 1 TxFrags). In the second, more stringent approach, we assembled TxFrags using the transcript assembly program Cufflinks (Set 2 TxFrags).
We first compared the characteristics of Set 1 TxFrags among annotated features including protein-coding genes, RNA genes, pseudogenes, and transposons. Regardless of the genomic feature and data set, the Set 1 TxFrag length distributions are bimodal with the first peaks located near the mRNA-seq single read length, indicating that most Set 1 TxFrags consist of a single read (Fig. 1A; Supplemental Fig. S1). Next, the fragments per kilobase of exon model per million mapped reads (FPKM) measure was used to assess Set 1 TxFrag expression level. Similar to length distributions, the Set 1 TxFrag FPKM distributions are bimodal with the first peaks at very low FPKM, mostly consisting of single-read TxFrags (Fig. 1B; Supplemental Fig. S2). The likely sources of low-FPKM TxFrags are (1) genes with very low-level or highly specific expression, (2) "transcriptional noise" representing background genome transcription (Struhl, 2007), or (3) low-level genomic DNA contamination in the sequenced mRNA sample. If the presence of one or more Set 1 TxFrags is considered evidence of expression, 78% to 94% of protein-coding genes are expressed. However, 19% to 68% of pseudogenes and 5% to 69% of transposons would be considered expressed based on the same criterion (Supplemental File S1A). Given that Arabidopsis transposons have been documented to be underexpressed (Schmid et al., 2005) and subject to strong posttranscriptional silencing through DNA methylation (Zhang et al., 2006; Zilberman et al., 2007), transposon expression was used as a conservative error estimate of expression calls.
To stringently control for false positives arising from background transcription and/or low-level genomic contamination, we applied multiple FPKM thresholds defined according to the percentage of transposon TxFrags considered expressed (Fig. 1C; Supplemental File S1A). Comparing percentage of transposons expressed, we found that the FPKM thresholds have significantly different impacts on data sets (Fig. 1C). For example, an FPKM threshold based on the 90th percentile of the transposon expression distribution results in a 57.4% reduction of transposon expression in the 12-d seedling drought stress data set compared with no FPKM threshold but causes no reduction in the 7-d data (Fig. 1C). This difference in the degree of transposon expression due to FPKM threshold choice is not simply due to differences in sequencing depth, as the numbers of mapped reads are both approximately 4.8 × 106 (Supplemental Table S1). In addition, this difference cannot be attributed to stress treatments, as degrees of transposon expression in the stress treatment and control samples are similar regardless of FPKM thresholds (Fig. 1C, green). We note that only 22% to 38% of reads from the 12-d data sets can be mapped to the Arabidopsis genome compared with 71% and 75% for 7-d and flower data, respectively. Thus, data quality may significantly impact gene expression calls, even after quality filtering and mapping the reads to the genome.
For comparison, we applied a second, more stringent transcript assembly approach using Cufflinks with bias corrections of transcript models based on sequences, positions, and abundance (Trapnell et al., 2010) to generate Set 2 TxFrags. Compared with Set 1 TxFrags, Set 2 TxFrags are significantly longer (Fig. 1D; Kolmogrov-Smirnov [KS] test, P < 2.2e-16) and have significantly higher FPKM values (Fig. 1E; KS test, P < 2.2e-16). In addition, Set 2 TxFrags length and coverage distributions overlap with the right tails of Set 1 TxFrags (Fig. 1, A, B, D, and E), indicating that the main difference between these two sets is enrichment for longer and more abundant transcripts in Set 2 TxFrags. Increasingly stringent FPKM thresholds still have a significant effect on the numbers of transposons considered expressed for several 12-d data sets (Fig. 1F). Nonetheless, the second approach allows for better control in calling transposon expression, which we considered to be mostly false positive, than the first, simpler approach.
Pervasiveness of Transcription in Arabidopsis Intergenic Regions
Previous microarray-based studies in Arabidopsis have shown that a large number of polyA transcripts are produced from the intergenic regions of the genome (Yamada et al., 2003; Matsui et al., 2008). Considering the advantages of RNA-seq over microarrays for expression studies (Agarwal et al., 2010), we reassessed the preponderance of intergenic transcription using RNA-seq data sets. Here, TxFrags located within intergenic regions are referred to as ITFs. We found that the analysis method (Set 1 versus Set 2), FPKM threshold, and data set significantly influence estimates of ITF abundance (Fig. 1, G and H; Supplemental File S1). For example, 9.2% of Set 1 TxFrags from the flower data are considered ITFs when no FPKM threshold is applied, but this proportion drops to 3.7% with an FPKM threshold of 1.33, which corresponds to a 10% false-positive rate (Fig. 1G). Comparing between data sets by allowing a 10% false-positive rate, ITF estimates differ by 7-fold (2.3%–16.4%) and 73-fold (0.2%–14.6%) based on Set 1 TxFrags and Set 2 TxFrags, respectively (Supplemental File S1). Despite these differences, there are two consistent characteristics among data sets that separate Set 1 and Set 2 ITFs. Set 1 ITFs tend to be significantly shorter than Set 2 ITFs (Fig. 1, A and D; KS test, P < 2.2e-16). In addition, Set 1 ITF expression levels are not significantly different from those of Set 1 TxFrags (Fig. 1B; KS test, P = 0.28) but are significantly lower than protein-coding gene TxFrags (KS test, P < 2.2e-16). Set 2 ITFs have significantly lower expression levels than protein-coding gene TxFrags as well (Fig. 1E; KS test, P < 1e-2), although the pattern is not as pronounced as for Set 1 ITFs, presumably due to the bias corrections applied on the data set by Cufflinks. Our findings are consistent with earlier studies in Arabidopsis (Hanada et al., 2007; Matsui et al., 2008) and mammals (Wang et al., 2004; van Bakel et al., 2010), which found that intergenic sequences tend to be lowly expressed.
We next focused on Set 2 TxFrags, which represent a more stringently defined set of transcripts. Across data sets, 0.2% to 14.6% of TxFrags are potentially derived from intergenic transcription based on a 5% false-positive rate (Fig. 1H; Supplemental File S1B). This proportion corresponds to 10,511 ITFs across eight RNA-seq data sets, together representing 6,545 nonoverlapping intergenic transcribed genomic regions and spanning 3.6% of the assembled intergenic region in Arabidopsis. Our ITF estimate is comparable to an earlier tiling array-based study in Arabidopsis, where 7,719 unannotated transcriptional units were defined as novel, non-protein-coding RNAs (Matsui et al., 2008). Other studies have provided more conservative estimates of Arabidopsis intergenic expressed regions, from 104 (Stolc et al., 2005) to 2,397 (Yamada et al., 2003). In mammals, however, the ENCODE project as well as other studies have reported significantly more pervasive intergenic transcription (Bertone et al., 2004; Birney et al., 2007; Kapranov et al., 2007). The ENCODE project reported that 488,906 (22.6%) TxFrags lie in intergenic regions and that 93% of the ENCODE bases have transcription evidence (Birney et al., 2007). Compared with Arabidopsis (3.6%), a significantly larger proportion of the ENCODE region is transcribed, even if we consider Set 1 TxFrags (13.7% at a 5% false-positive rate) that are not as rigorously defined as Set 2 TxFrags.
There are several possible explanations for the differences in ITF pervasiveness between plants and humans. First, the ENCODE study analyzed transcripts obtained from 31 different cell lines and tissues, which represents a much broader sampling of the transcriptome than our study. Second, known issues with tiling arrays used in the ENCODE study, particularly cross hybridization (Agarwal et al., 2010; van Bakel et al., 2010), may lead to an overestimation of ITFs. Consistent with this possibility, a previous RNA-seq study of human 293T cell total RNA found that only approximately 4% of reads were intergenic (van Bakel et al., 2011), similar to our Arabidopsis estimate. Third, the intergenic space in Arabidopsis constitutes only approximately 40% of the genome, compared with approximately 99% in the human genome. If intergenic transcripts are largely derived from noisy transcription or genomic contamination, species with larger genomes may have more RNA-seq reads from intergenic space. The fourth reason may be that larger genomes have more functional elements. However, variation in genome size can be due to extreme proliferation of transposable elements (Hawkins et al., 2006; Piegu et al., 2006). Thus larger genomes do not necessarily contain more genes. Finally, elements of our experimental design, such as the use of tissue samples with multiple cell types or insufficient coverage, may lead to an underestimate of ITFs. To address some of the issues concerning our study design, we analyzed cell type-specific transcriptome data obtained using directional Illumina sequencing.
Factors Affecting ITF Estimates
The data sets we analyzed have the following limitations that may affect estimates of ITF abundance (Clark et al., 2011). First, all data sets were generated using complex tissue samples that may render cell type-specific ITFs undetectable. Second, the sequencing was performed using single reads without directionality information, which may result in misassembly of ITFs. Third, the read length and coverage may be insufficient for detecting ITFs expressed at low levels. To address these issues, we directionally sequenced polyA-selected RNA from T87 suspension culture cells with longer reads (72 bp) and greater depth (two to nine times more sequenced bases; approximately 2.3 Gb, approximately 3 × 107 reads). We found that 0.9% of reads and 4.2% of TxFrags (identified using the same criteria as Set 2 TxFrags) from the suspension culture data are intergenic (Supplemental Fig. S3), consistent with the proportions of intergenic reads and TxFrags identified from more complicated tissues (Supplemental Table S1). In addition, in the suspension cell data set, 3,052 (9.8%) TxFrags and 170 (13.1%) ITFs overlap with one or more other TxFrags and ITFs, respectively, that are in the opposite orientation. Thus, lack of read directionality information and misassembly can lead to an approximately 13% underestimate of ITFs that overlap in opposite orientations.
Another factor affecting the estimate of ITFs is that the data sets we have analyzed so far are derived from polyA RNA. Nonpolyadenylated (nonpolyA) RNA may constitute the bulk of the transcriptome and significantly contribute to intergenic expression (Cheng et al., 2005; Armour et al., 2009; Xu et al., 2010). However, an earlier study focusing on both polyA and nonpolyA RNAs in Arabidopsis found that 3.5% of reads are intergenic (Lister et al., 2008), which is comparable to the 0.4% to 8.2% reads that are intergenic in the mRNA-seq data sets we analyzed (Supplemental Table S1). In addition, a study of human 293T cell ribosomal RNA-depleted total RNA revealed that approximately 4% of reads were intergenic (van Bakel et al., 2011). This suggests that our estimates of intergenic transcription in Arabidopsis based on polyA RNA sequencing are reasonable. Nonetheless, detailed studies of nonpolyA ITFs will be necessary to estimate the contribution of nonpolyA transcripts to intergenic transcription.
Taken together, we have identified 6,545 ITFs (5% false-positive rate) that are likely novel transcriptional units not previously defined in the Arabidopsis genome. Two outstanding questions remain. First, because these ITFs are derived from polyA RNAs, are they parts of novel protein-coding or noncoding RNA genes? Second, do some of these ITFs have clear evidence of selection, therefore suggesting their functionality? To address the first question, we assessed the protein-coding potential of ITFs by analyzing ribosome-associated transcripts and shotgun proteomics data sets.
Distinguishing Coding from Noncoding Intergenic Transcripts Based on Ribosome Association
Translation initiation is the rate-limiting step in protein translation; therefore, transcripts associated with the ribosome are more likely to be translated (Kawaguchi and Bailey-Serres, 2002; Bailey-Serres et al., 2009). Studies in Arabidopsis (Branco-Price et al., 2008; Jiao and Meyerowitz, 2010), mouse (Doyle et al., 2008), and yeast (Ingolia et al., 2009) have taken advantage of this property to globally investigate translational regulation. To assess whether ribosome association of intergenic transcripts is a good measure of their translation potential, we first sequenced ribosome-associated transcripts from 7-d-old seedlings. After identifying the ribosome-associated transcript fragments (R-TxFrags), we selected eight genomic regions with evidence of ribosome association and seven without for in vivo translation studies (Supplemental Table S2). These regions overlap with putative small open reading frame (sORF) genes that were originally computationally predicted from intergenic regions (Hanada et al., 2007). Several of these regions have since been annotated solely based on computational predictions and/or complementary DNA (cDNA) evidence. The 5′ untranslated regions (UTRs) and coding sequences of the sORFs were fused in frame to a yellow fluorescent protein (YFP) reporter that lacks a translational start codon, and the translation of these sequences in transiently transformed tobacco (Nicotiana tabacum) leaf epidermal cells was evaluated (see "Materials and Methods"; Fig. 2).
Of the eight genomic regions with R-TxFrag evidence, five were translated in tobacco, while only one of the seven regions without R-TxFrag support was translated (Fig. 2; Supplemental Table S2). Thus, there was a significant enrichment of sORFs with R-TxFrag evidence among those translated in vivo (Fisher's exact test, P < 0.05). The observed localization patterns of the protein fusions were largely consistent with signal peptide predictions (Fig. 2; Supplemental Table S2), indicating that the fusion proteins were likely correctly translated and targeted in tobacco. However, three sORFs with ribosome association evidence do not appear to be translated in the transient expression assay. These sORFs may not be translated, or this may be an artifact due to the use of a heterologous system (tobacco). One translated sORF is an annotated "other RNA" gene (At1g31935; Fig. 2I; Supplemental Table S2). In addition, a number of annotated other RNA genes have either ribosome association or proteomics evidence, which highlights the importance of experimentally evaluating protein-coding potential. Overall, based on the findings of our in vivo translation assays, we conclude that features with evidence of ribosome association are more likely to be translated than those without.
Translation Evidence for ITFs
Given that ribosome association is a good indicator of the translation potential of intergenic sequences, we further analyzed R-TxFrags from the 7-d-old seedling data to estimate the proportion of ITFs likely to be parts of coding genes. To address potential issues due to sequencing coverage or tissue-specific expression and translation, R-TxFrags were also identified using ribosome-associated transcript data of whole flowers and specific floral domains expressing three homeotic genes (Jiao and Meyerowitz, 2010). For comparison, we also incorporated shotgun proteomics data from two studies examining protein expression in multiple tissues and developmental stages (Baerenfaller et al., 2008; Castellana et al., 2008). These data are collectively referred to as "translation data sets" and are summarized in Supplemental Table S3.
As with mRNA-seq, most R-TxFrags and proteomics tags mapped to previously annotated regions in the genome, particularly protein-coding genes. Among ribosome immunoprecipitation data sets, 63% to 73% of protein-coding genes have one or more R-TxFrags (Fig. 3A). Similarly, 74% of protein-coding genes have one or more proteomics tags (Fig. 3A). In addition, 62% to 67% of the R-TxFrags overlap with one or more proteomic tags. These findings demonstrate that ribosome-associated transcripts tend to be translated, consistent with our in vivo translation studies (Fig. 2). On the other hand, 5% to 23% of annotated ncRNA genes and 7% to 15% of pseudogenes have uniquely mapped R-TxFrags and/or proteomics tags. If all annotated RNA genes are truly noncoding, calling a feature translated based on a corresponding R-TxFrag and/or proteomics tag can have a 5% to 23% false-positive rate depending on the data set (Fig. 3A, RNA). One anomaly is that 34.0% of transposons have proteomics tags, although only 1.9% to 3.4% have R-TxFrags (Fig. 3A). This discrepancy is in sharp contrast to our finding that the proportions of protein-coding genes possessing R-TxFrags and proteomics tags are both approximately 70% (Fig. 3A). This observation, also noted in the original study (Castellana et al., 2008), is inconsistent with studies demonstrating reduced transcription of transposons (Schmid et al., 2005) and their extensive methylation (Zhang et al., 2006; Zilberman et al., 2007). Using the number of proteomics tags as a proxy for protein expression level, transposons with proteomics evidence tend to have significantly fewer tags than protein-coding genes (Supplemental Fig. S5; KS test, P < 2.2e-16). In addition, 67.6% of transposons with proteomics evidence have only one tag compared with 26.6% of protein-coding genes, suggesting that if the transposons are expressed and translated, it happens at significantly lower levels than protein-coding genes.
How many ITFs have evidence of translation? Among the 6,545 nonoverlapping ITFs identified from eight mRNA-seq data sets, 2,107 (32.2%) have one or more R-TxFrags from one or more of the translation data sets analyzed (Fig. 3A; Supplemental Fig. S4). Unlike protein-coding genes, there is substantially stronger support for ITF translation from ribosome association than from proteomics data (Fig. 3A). Some of the ribosome-associated ITFs may contain protein-coding regions even though there is no proteomics support, partly due to the fact that proteomics data tend to be biased toward more abundantly translated proteins (Baerenfaller et al., 2008). It is also likely that a significant number of ribosome-associated ITFs are derived from the UTRs of protein-coding transcripts. Taken together, even if the false-positive rate is 23%, approximately 1,622 ITFs are likely parts of transcripts destined to be translated after eliminating potential false positives. Thus, a significant number of intergenic transcripts may be part of larger protein-coding genes, either as coding sequences or as UTRs. Our finding highlights the importance of assessing the translational potential of polyA intergenic transcripts before defining them as sequences that function solely at the RNA level.
Relationship between ITFs and Neighboring, Annotated Genes
Based on an analysis of mRNA sequencing data, we uncovered thousands of short, low-abundance transcripts from intergenic regions. In addition, many of these ITFs are supported by translation evidence. One immediate question is whether these ITFs, translated or not, are extensions of previously annotated or novel protein-coding genes. To address this question, we assessed whether there is a significant bias in where ITFs are located within the Arabidopsis genome by calculating the distance between each ITF and its closest annotated protein-coding gene. We found that although a substantial number of ITFs are closer to genes, they are not any closer than intergenic sequences sampled randomly based on ITF number and size (Fig. 3C). This is contrary to the expectation that ITFs are predominantly extensions of existing genes.
Given that ITFs in general are not closer to neighboring genes than randomly selected intergenic sequences, do ITFs with translation evidence behave similarly? First, we found that translation evidence (proteomics tags and R-TxFrags) tends to lie farther away from protein-coding genes than random expectation (Supplemental Fig. S6A). However, ITFs with translation evidence tend to lie closer to genes than ITFs without translation evidence (Supplemental Fig. S6B), suggesting that most ITFs with translation evidence may be parts of neighboring protein-coding genes. If translated ITFs are indeed missing parts of annotated genes, ITFs closer to genes with transcription evidence should be enriched in the translated set compared with ITFs closer to nontranscribed genes. Consistent with this expectation, among the 4,942 ITFs closest to a transcribed annotated protein, 37.9% have translation evidence, while among the 563 ITFs closest to a nontranscribed annotated protein, only 14.2% have translation evidence (Fisher's exact test, P < 2.2e-16; Fig. 3D). These observations suggest that most ITFs with translation evidence that are close to annotated genes may be missing parts of those genes or associated with the transcription of those genes via an unknown mechanism.
Taken together, we have demonstrated the presence of ITFs from more than 6,000 intergenic regions in Arabidopsis from multiple RNA sequencing data sets. More than 20% of these ITFs are likely translated or are part of protein-coding transcripts. Among the 6,545 ITFs, 59.4% are located more than 300 bp away from an annotated gene. Of these, 847 (21.7%) have translation evidence. Considering that 300 bp is approximately the 90th percentile of both Arabidopsis intron and UTR lengths, these relatively distant ITFs may be parts of novel transcriptional units. However, ITFs, in general, tend to be significantly shorter and expressed at lower levels than protein-coding genes. We also find that ITFs, in general, tend to be expressed narrowly, in a data set-specific manner, while TxFrags corresponding to annotated features are present in multiple data sets (Fig. 3B). The translation of approximately 32.2% of ITFs is supported by one or more ribosome immunoprecipitation and/or proteomics data sets, compared with 88.0%, 44.6%, and 36.9% for protein-coding genes, pseudogenes, and transposons, respectively. In terms of translation, ITFs behave similarly to pseudogenes and transposons. Previous studies have suggested that the breadth of expression as well as the level of expression can be considered as proxy indicators of functionality (Nuzhdin et al., 2004; Subramanian and Kumar, 2004; Movahedi et al., 2011). However, genes can have highly specific expression and/or low expression levels. Thus, one remaining question is whether these ITFs are parts of functional sequences with clear evidence of selection.
Evidence of Natural Selection on ITFs at the Nucleotide Level
Intergenic transcripts that are independent transcriptional units may be derived from noisy, background transcription or unannotated genes that are functional (Struhl, 2007; Dinger et al., 2009; van Bakel et al., 2010, 2011; Clark et al., 2011). Transcripts not important to cellular function are expected to accumulate mutations much like neutrally evolving sequences. In contrast, functional ITFs should be selected for and show signs of nonneutral evolution. To assess whether there is a clear signature of natural selection that is indicative of functionality, the ITF nucleotide substitution rates were estimated using syntenic genomic regions of Arabidopsis and Arabidopsis lyrata, which diverged from their common ancestor approximately 10 million years ago (Hu et al., 2011). We also estimated the 4-fold degenerate site substitution rates of protein-coding orthologs as proxies for neutral evolution rates. Substitution rates for protein-coding genes, RNA genes, and randomly chosen intergenic regions not overlapping with ITFs were also estimated for comparison.
Among 6,545 ITFs, only 1,238 (18.9%) have identifiable syntenic regions for substitution rate estimation between the two Arabidopsis species. The proportion of syntenic ITFs is significantly lower than those of protein-coding genes (90.9%; Fisher's exact test, P < 2.2e-16), RNA genes (35.7%; P < 9.1e-10), and pseudogenes (33.4%; P < 2.2e-16) but significantly higher than transposons (10.9%; P < 1.5e-16). Thus, many "orphan" ITFs without putative orthologs likely evolved rapidly with little or no selective constraint. For ITFs found within syntenic regions, substitution rates are significantly higher than those of annotated protein-coding genes (KS test, P < 2.2e-16; Fig. 4A). On the other hand, ITF substitution rates in general are significantly lower than those of 4-fold degenerate sites (KS test, P < 2.2e-16; Fig. 4A). These observations suggest that ITFs may constitute a mixed population, with the first population under strong selective constraint and the second one evolving neutrally. Using the 5th percentile of the 4-fold degenerate site rate distribution (rate = 0.07) as a threshold, only 6.4% of the 6,545 ITFs are likely under strong purifying selection. The remaining 93.6% are likely under little or no purifying selection. To control for local rate variation, we compared the rate of each ITF with the 4-fold site rates of neighboring genes. Based on this approach, a much smaller percentage, 2.7%, of ITFs were found to be under selection (Fig. 4B).
One issue in comparing any sequence feature with 4-fold sites is that there can be significant alignment bias. This is because a sequence feature (e.g. an ITF) is aligned to its putative ortholog at the nucleotide level, whereas 4-fold sites are identified from nucleotide sites originally aligned based on protein sequences. Given that the alignment process involves finding an alignment with the best score, regardless of whether the sites are homologous or not, it will tend to make a nucleotide-based alignment look more similar than it really is. Thus, the lower substitution rate among sequence features compared with 4-fold sites can simply be due to this artifact. To account for this, we selected random intergenic regions that have no evidence of expression and calculated their substitution rates. Similar to ITFs (6.4%), 7.5% of random intergenic region samples are under strong purifying selection. More importantly, there is little, but statistically significant, difference in the substitution rate distributions between ITFs and random intergenic sequence (median rates, 0.09 and 0.08, respectively; KS test, P < 2e-05). Therefore, after accounting for potential alignment bias, potentially even fewer than 6.4% of ITFs are under significant selective constraints. The implication is that, based on cross-species comparison, the majority of ITFs appear to evolve in a way similar to presumably nonfunctional, nonexpressed random intergenic regions.
To assess the possibility that some ITFs may have a species-specific function in Arabidopsis, making the signature of selection only obvious at the intraspecific level, we analyzed genomic sequences of 80 accessions of Arabidopsis (Cao et al., 2011) and estimated the nucleotide diversity (π) for ITFs, annotated sequence features, and randomly selected intergenic sequences not overlapping with ITFs. The π value allows us to assess the genetic variability of different genomic features among Arabidopsis populations (Li, 1997). Our findings suggest that ITFs have π values significantly higher than coding sequences and RNA genes (KS tests, both P < 2.2e-16) but similar to random intergenic sequences not overlapping with ITFs (Fig. 5). We also estimated Tajima's D, Fu and Li's D, and Fay and Wu's H statistics based on site-frequency spectrum to assess if ITFs are under selection. The distributions of all three statistics were comparable between ITFs and randomly sampled, unexpressed intergenic sequences (KS tests, all P > 0.1) but significantly different from those of protein-coding genes and RNA genes (KS tests, all P < 0.001; Supplemental Fig. S7). These findings suggest that there is much more relaxed selection within species on ITFs compared with protein-coding genes. Furthermore, the intensity of selection on ITFs is similar to that on random intergenic regions, which are likely largely nonfunctional and evolve neutrally.
Selection on ITFs with Translation Evidence
Our findings indicate that some ITFs are under strong selective constraint and may be functional. However, a much larger number of ITFs do not have clear signatures of selection. One immediate question is whether ITFs under strong selective constraint tend to be those that are translated, given that ITFs with translation evidence tend to be located closer to neighboring genes. We performed a similarity search between ITF sequences and the genomes of 15 land plants ranging from a bryophyte to angiosperms. For comparison and to address potential annotation issues, we also analyzed TxFrags mapping to protein-coding genes and RNA genes. ITFs with evidence of translation were slightly more conserved over ITFs with no evidence of translation (compare Fig. 4, C and D), but most ITFs have significant similarities only between Arabidopsis and A. lyrata, with sequence similarity rapidly declining beyond the Arabidopsis genus. On the other hand, there is a significantly higher degree of cross-species similarity between protein-coding genes: 6,895 of the 10,000 randomly selected protein sequences had E values of less than 1e-5 in more than one species (Fig. 4E). Even RNA genes, which are not expected to be translated, have higher sequence similarities than ITFs (Fig. 4F). Thus, at both the nucleotide and amino acid sequence levels, relatively few ITFs are under selection based on cross-species comparisons.
Of the 847 ITFs with translation evidence that are located more than 300 bp away from genes, 799 did not show significant conservation. Considering that these ITFs tend to be expressed at low levels, such sequences may represent translational noise. But we cannot rule out the possibility that they are lineage-specific coding sequences. Of the 49 ITFs that did show conservation, 16 had similar sequences present in more than one species at the amino acid level (E value < 1e-5), and 10 showed overlap with computationally predicted sORFs with high protein-coding potential (Hanada et al., 2010). The π value distribution of these 49 ITFs is statistically indistinguishable from that of TxFrags mapping to protein-coding sequences (KS test, P = 0.2; Fig. 5). These ITFs may thus represent novel functional genes. Nonetheless, only approximately 5% of ITFs with translational evidence are subject to strong purifying selection among Arabidopsis accessions, reinforcing the notion that most of them are products of noisy transcription.
In this study, we analyzed the intergenic polyA transcriptome of Arabidopsis to address the issues of abundance, coding/noncoding nature, and functional relevance of intergenic polyA transcripts. Our results indicate that approximately 5% of the TxFrags in the Arabidopsis transcriptome can be reliably called intergenic. One limitation of our analyses, as we have noted before, is our focus on the polyA fraction of the transcriptome. It is likely that the nonpolyA fraction of the transcriptome may harbor additional novel noncoding genes that need to be further investigated. Another limitation is that the read lengths of the RNA-seq data we used are short. It is possible that some ITFs belong to the same transcriptional units, making the number of ITFs an overestimate.
Our results indicate that approximately 3.6% of the intergenic space in Arabidopsis is transcribed by RNA polymerase II, and approximately 40% of what is transcribed tends to lie within 300 bp of annotated genes. Around one-third of ITFs have translation evidence, and we find a significant bias in their distribution; they tend to be closer to transcribed protein-coding genes, raising the possibility that some ITFs may in fact be unannotated extensions of known genes. Our primary sequence-level evolutionary analysis indicates that a relatively low fraction (approximately 5%) of the ITFs have experienced strong purifying selection either within species or between species. We should emphasize that our criteria for evaluating selection is stringent. Furthermore, some ITFs may be more strongly constrained at the secondary structural level, similar to noncoding RNA genes (Washietl et al., 2005). In addition, some long noncoding RNAs such as Air and Xist are poorly conserved (Pang et al., 2006; Ponting et al., 2009), indicating that a lack of conservation may not always mean lack of function. Thus, it is likely we will miss some ITFs that are under selection or are functional. However, most ITFs are short, and unlike the long noncoding RNAs, tend to be data set specific and expressed at a very low level compared with annotated genes. In addition, most ITFs have characteristics more similar to pseudogenes and transposons than to protein-coding and RNA genes. Taken together, most ITFs bear the hallmarks of neutrally evolving sequences, suggesting that they are products of noisy transcription, as proposed earlier (Struhl, 2007).
The idea of transcriptional noise has been intensely debated over the past few years. Some studies support the theory that the transcriptional machinery might be error prone and that many transcripts may be the result of false starts and/or stops (Li et al., 2007; Struhl, 2007; Xu et al., 2009; van Bakel et al., 2010). Such errors may occur because it is not possible to regulate any biological process to the point that there is no error; noisy transcription may exist simply because it incurs little fitness cost. Considering that the vast majority of mutations are neutral or nearly neutral (Ohta, 1992), this paradigm for gene evolution may also apply to other molecular events, including transcription. As has been postulated before, the target of natural selection may be the effects of error-prone transcription rather than the transcriptional process itself (Hurst, 2009). Another possibility is that the effect of genetic drift, particularly on organisms with smaller effective population sizes, may render selection against erroneous transcript production ineffective. In either case, the hypothesis is that transcriptional errors may not always be subjected to purifying selection. Based on our findings, it would seem that much of the intergenic transcription falls into this category. We should emphasize that, for most ITFs, there is little or no evidence to reject the null hypothesis that ITFs are nonfunctional. The reason to consider nonfunctionality as null is simply because only functionality can be experimentally tested (van Bakel et al., 2010).
In a recently published series of papers by the ENCODE consortium, 80.4% of the human genome was found to have evidence of functionality (Dunham et al., 2012). Functionality of a sequence, in this case, was defined at the biochemical level. That is, a functional sequence is presumed to have at least one RNA- and/or chromatin-associated event, such as transcription factor binding, nucleosome binding, or DNA methylation, in at least one cell type. However, it remains unclear to what extent these biochemically functional sequences may have physiological function, given that these biochemical events can also be due to noise. As an example, among the novel long intergenic polyA TxFrags obtained in the ENCODE study, only 4% of the bases were conserved between humans and macaques. Among the 96% of bases not conserved, only 6% to 11% showed evidence of lineage-specific constraint in humans, comparable to what we found in Arabidopsis (Ward and Kellis, 2012). In addition, the ITFs found in this study were found to be present at less than 0.1 copy per cell (Djebali et al., 2012), consistent with our finding that most of the Arabidopsis ITFs have a very low expression level.
Considering the dramatically increased sensitivity and throughput in sequencing, noisy transcription and even contaminating sequences, such as trace amounts of genomic DNA, can be readily detected. Thus, a transcription event as detected by sequencing may not be considered functional by default. Evolutionary constraint acting on novel sequences or other evidence should be demonstrated prior to their annotation. The increasing availability of population-wide polymorphism data sets and genome sequences of related species provide more robust tools for such evolutionary studies, especially those focusing on lineage-specific selection (Ward and Kellis, 2012). In addition to the question of functionality, we show that a significant number of ITFs are associated with ribosomes and a smaller fraction of them have proteomics tags. Thus, novel transcripts should not be regarded as noncoding by default without rigorous experimental analysis of their coding potential. Our results also suggest the need to have a clearer understanding of the mechanistic aspects of RNA polymerase action on how noisy transcription may arise. Use of an integrated approach to validate novel RNA predictions and their functionality would be important in this regard. For example, a previous study in mouse used an array of approaches, including the identification of conserved histone modification marks, evolutionary analyses of promoter regions, gene set enrichment analysis, transcription factor chromatin immunoprecipitation, and RNA interference assays, to identify putative functional long ncRNAs (Guttman et al., 2009). We surmise that such an approach will allow us to explore more deeply the mechanistic and evolutionary aspects of the transcriptional process in plants.
Plant Material and RNA Isolation
For transcriptome and ribosome immunoprecipitation studies, transgenic Arabidopsis (Arabidopsis thaliana) Columbia seeds expressing a His6FLAG-tagged version of the ribosomal large subunit protein L18B (35S:HF-RPL18B) were surface sterilized, stratified for 3 d, and sown on 0.5× Murashige and Skoog medium containing 1% (w/v) Suc and 0.4% phytagel. Seedlings were grown vertically under a 16-h-day (125 μE m−2 s−1 photosynthetically active radiation)/8-h-night cycle for 7 d as described previously (Branco-Price et al., 2008). Seven-day-old seedlings were harvested at the end of the light period. Total RNA extraction and ribosome immunoprecipitation were done as described previously (Branco-Price et al., 2008) for three biological replicates, except that the RNeasy Plant Mini Kit purification step was omitted. Total RNA and ribosome-immunoprecipitated RNA were quantified using a Nanodrop spectrophotometer (Nanodrop Technologies), and RNA quality was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies).
Illumina RNA-seq and Data Analysis
For each in-house (7-d seedlings) RNA sample, cDNA libraries were constructed by first isolating polyA RNA from 3 µg of total or ribosome-associated RNA. Libraries for RNA-seq were prepared using the Illumina mRNA-seq Sample Prep Kit. Briefly, polyA RNA was fragmented and reverse transcribed using random primers. Adapters were ligated to double-stranded cDNA, and fragments from 175 to 225 bp were gel purified. After PCR amplification, the cDNA libraries were sequenced on an Illumina Genome Analyzer. Each library was loaded onto at least two lanes; however, usable sequence was only obtained for one polyA RNA library and two ribosome-associated RNA libraries (four lanes total). Three lanes of ribosome-associated RNA sequencing, corresponding to two biological replicates, were combined to give 19,818,643 reads, and one lane of polyA RNA yielding 7,028,772 36-bp reads was obtained. The original sequencing reads were deposited in the National Center for Biotechnology Information Short Read Archive under accession number SRA053376 (http://www.ncbi.nlm.nih.gov/Traces/sra). All public data sets (Supplemental Table S1; Supplemental Fig. S2) were downloaded from the National Center for Biotechnology Information Short Read Archive. The following procedure was commonly performed on both the in-house and public data sets.
The short reads, after quality trimming, were mapped to The Arabidopsis Information Resource Arabidopsis genome release 10 using Bowtie version 0.12.7 (Langmead et al., 2009) and TopHat version 1.2.0 (Trapnell et al., 2009). The default settings were used except that the maximum combined intron size was set at 5,000 bp. The mapped reads were assembled with two approaches. In the first approach, reads with overlapping genomic locations were merged into TxFrags (Set 1 TxFrags) without considering the possibility that neighboring TxFrags may be derived from the same transcriptional units. In the second approach, Cufflinks 0.9.3 (Trapnell et al., 2010) was used with default parameters, except that a maximum combined intron size was set at 5,000 bp (Set 2 TxFrags). All TxFrags overlapping with annotated features by 1 bp or more, including those in introns or UTRs, were flagged as genic transcripts.
Estimating Expression Level and Breadth of Expression of Features
For estimating expression level, the FPKM measure was used. Since Set 1 TxFrags represent a set of unique, nonoverlapping TxFrags, the entire TxFrag was considered an exon for the purpose of FPKM estimation. For Set 2 TxFrags, FPKM values were estimated by Cufflinks. The breadth of expression was calculated for the 6,545 merged ITFs. For comparison, we also measured the expression breadth of TxFrags mapping to annotated features. We used the number of data sets in which a particular feature had expression evidence (one or more overlapping TxFrag) as a measure of the breadth of expression of that feature.
5′ RACE and Transient Expression of YFP Fusion Proteins in Tobacco
The 5′ UTRs of putative coding sORF sequences were identified from publicly available cDNA sequences (Aubourg et al., 2007) or were amplified by 5′ RACE (RLM-RACE kit [Ambion] or SMART RACE cDNA amplification kit [Clontech]). The 5′ UTRs and coding sequences of each sORF were amplified from genomic DNA and cloned into the TOPO-TA entry vector (Invitrogen). The sequences were then transferred by recombination mediated by LR Clonase (Invitrogen) into a modified pMDC83 destination vector (Curtis and Grossniklaus, 2003), containing the enhanced YFP sequence (Clontech) and lacking a translational start codon, under the control of the 35S promoter. Constructs containing sORFs fused in frame with YFP were transformed into Agrobacterium tumefaciens GV3101. Transient transformation was performed to express sORF-YFP fusions in tobacco (Nicotiana tabacum) cells (Sparkes et al., 2006). Transgenic A. tumefaciens cells were cultivated overnight, and 200 μL of the culture (optical density A600 approximately 1–2) was pelleted and resuspended with sterile water to 0.1 optical density. A. tumefaciens cells were infiltrated into tobacco leaves, and the infiltrated tobacco was kept under constant light for 72 h. Infiltrated areas of tobacco leaves were detached and observed with an inverted laser scanning confocal microscope (Olympus Spectral FV 1000). YFP signals were detected with the 514-nm argon laser excitation line with a band-pass emission filter of 517.5 to 542.5 nm. For visualization of AT_1|-|2|5786755-5786853 (Hanada et al., 2007) and Endoplasmic Reticulum Retention Defective2 (ERD2) colocalization, equal volumes of A. tumefaciens cultures were mixed prior to infiltration. Fluorescence was visualized after 3 d with a Meta Zeiss confocal microscope using the argon laser excitation lines of 458 and 514 nm and band-pass emission filters of 475 to 525 nm and 530 to 600 nm for blue-shifted GFP and YFP, respectively.
Evolutionary Conservation Analyses
To identify ITFs under selection at the nucleotide level, the Arabidopsis ITFs were first mapped to the Arabidopsis lyrata genome using GMAP version 2007-09-28 (Wu and Watanabe, 2005) with default settings. Putative ITF orthologs were defined as pairs of similar sequences (80% or greater coverage, 80% or greater identity, 40 bp or greater match length) between Arabidopsis and A. lyrata flanked by one or more putative orthologous genes among 10 protein-coding genes on either side of the ITF. Putative orthologs between these two species were identified based on reciprocal best-match and synteny information. The orthologous ITFs were aligned using Clustal 2.1 (Thompson et al., 1994), and the nucleotide substitution rate was calculated using baseml with the HKY substitution model in PAML (Yang, 2007). ITFs with a substitution rate lower than the 95th percentile (HKY distance ≤ 0.07) of the 4-fold degenerate site substitution rates of all protein-coding orthologs were deemed to be evolving under strong purifying selection. To control for genome-wide variation in local substitution rates, the 4-fold degenerate site substitution rates of up to 60 protein-coding genes in the vicinity of the ITFs were used to determine the 5% significance level using a z test. We did not conduct a nonsynonymous substitutions per nonsynonymous site (Ka)/synonymous substitutions per synonymous site (Ks) analysis for ITFs because (1) most ITFs are short, so the variance of Ka and Ks estimates for short sequences tend to be high, and (2) it is not clear what the correct reading frame is, if these ITFs are translated. Instead, to compare the levels of conservation at the coding sequence level, we performed a translated BLAST search between ITF/TxFrag sequences and the draft assemblies of 14 plant species in Phytozome 5.0 (http://www.phytozome.org/). The negative logarithm of the E value of the top match in each species was used to plot a heat map. All negative log values of 10 or greater or 0 or less were set to 10 and 0, respectively.
For conservation analyses within species, we used polymorphism data in the form of a genome matrix file from 80 different Arabidopsis accessions (Cao et al., 2011). For each genomic feature type, we reconstructed the aligned sequences based on the genome matrix file. The aligned sequences were analyzed for π, Tajima's D, and Fu and Li's D using Variscan (Vilella et al., 2005) with the following parameters: RefPos = 1, Outgroup = none, RunMode = 12, UseMuts = 0, CompleteDeletion = 0, FixNum = 1, NumNuc = 60. For Fay and Wu's H, we used the orthologs in A. lyrata as outgroups with RunMode = 22. For comparison, π values for 10,000 randomly chosen protein-coding genes, RNA genes, transposons, and pseudogenes were also calculated. For features with fewer than 10,000 sequences, bootstrap samples were used. To determine the background π values, 10,000 random intergenic sequences were sampled based on the size distribution of ITFs. Only those intergenic sequences not overlapping with any TxFrags were used for analysis. For each sequence in each feature type, a π value was estimated. The π distributions were then compared statistically.
The presence of ambiguous nucleotides or the short size of the ITFs can affect the error margins associated with π estimates. To assess whether these factors influenced our findings, we conducted additional analysis by changing the minimum number of sites analyzed (MinLength) and the proportion of the aligned length with nonambiguous bases (coverage). We sampled a range of MinLength (0, 50, 100, 150, and 200) at no coverage threshold (Supplemental Fig. S8) and a range of coverage (0, 0.25, 0.50, 0.75, and 1) at no MinLength threshold (Supplemental Fig. S9). Our analyses suggested that the trend observed in Figure 5 is not affected by the presence of ambiguous nucleotides or the short length of the ITF.
Sequence data from this article can be found in the GenBank/EMBL data libraries under accession number SRA053376.
Supplemental Data
The following materials are available in the online version of this article.
Supplemental Figure S1. Size distributions of TxFrags in different data sets.
Supplemental Figure S2. Expression level distributions of TxFrags in different data sets.
Supplemental Figure S3. Directional sequencing of mRNA from T87 cells.
Supplemental Figure S4. Pair-wise proportions of ITFs with translation evidence across different data sets.
Supplemental Figure S5. Number of peptides per feature.
Supplemental Figure S6. Distances of translation evidence from annotated genes.
Supplemental Figure S7. Population-level tests of neutrality for different features.
Supplemental Figure S8. Changing MinLength thresholds does not impact π distributions.
Supplemental Figure S9. Changing coverage thresholds does not impact π distributions.
Supplemental Table S1. Description of input data sets.
Supplemental Table S2. Sequences tested for coding potential.
Supplemental Table S3. Description of data sets used for translation analyses.
Supplemental File S1. Percentages of features considered expressed at different false-positive rate thresholds.
We are grateful to Kevin Carr, Shari Tjugum-Holland, and Jeff Landgraff for help with sequencing and data preprocessing, Melinda Frame for help with confocal laser scanning microscopy, John Johnston for computational support, Kelian Sun for help with cultivating T87 cells, and Andy Lin for his contributions toward preliminary studies.
The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantphysiol.org) is: Shin-Han Shiu (shius{at}msu.edu).
↵1 This work was supported by the National Science Foundation (grant nos. DEB–0919452, IOS–1126998, and MCB–0929100 to S.-H.S.).
↵[OA] Open Access articles can be viewed online without a subscription.
↵[W] The online version of this article contains Web-only data.
- polyA
- polyadenylated
- nonpolyA
- nonpolyadenylated
- intergenic transcribed fragment
- ncRNA
- noncoding RNA
- open reading frame
- mRNA-seq
- mRNA sequencing
- TxFrags
- transcript fragments
- fragments per kilobase of exon model per million mapped reads
- Kolmogrov-Smirnov
- R-TxFrags
- ribosome-associated transcript fragments
- sORF
- small open reading frame
- yellow fluorescent protein
- untranslated region
- π
- nucleotide diversity
- cDNA
- complementary DNA
- Encyclopedia of DNA Elements
- Received August 8, 2012.
- Accepted October 21, 2012.
- Published November 6, 2012.


text:
Beastie Bloomz is a pH balanced blend of extra strength nutrients specially formulated to produce heavy weight fruit and maximize blossom development during the critical last weeks of flowering. Beastie Bloomz is designed to create unbelievable internal bud density, which means that plants create fat, tightly packed blossoms and plenty of them. Use in soilless mixes, organic soils, hydroponics and aeroponic systems.
Beastie Bloomz is perfect for gardeners and growers who are ready to take their plants to the next level. This is a high-phosphorus soluble fertilizer that's geared towards building amazing buds and flowers. Our revolutionary micronutrient pack is designed to create unbelievable internal bud density, which means that plants create fat, dense, tightly-packed blossoms, and plenty of them. You'll also see exceptionally brilliant color and a long, dazzling blooming season.
Garden tip: Start with Open Sesame in the fifth or sixth week of a plant's life. Once bud set begins, feed your plants with Beastie Bloomz every other watering according to package directions and our FoxFarm Feeding Schedule. For high-yield results, start indoor plants in our Light Warrior Agro-Ponic Grow Medium, and use our premium Ocean Forest Potting Soil for container plants.
In the eighth to ninth week of the plants life cycle, switch to Cha Ching for an incredible harvest.
Item Number: 142649
Universal Product Code: 752289501616
Manufacturer: FoxFarm Soil & Fertilizer Company
Manufacturer's Part Number: 732764
Brand: FoxFarm Soil & Fertilizer Company
|MSDS Sheet

subset: commoncrawl


text:
18ct white gold halo design diamond ring set with one fancy yellow round brilliant cut diamond in the centre and grain set round brilliant cut diamonds flanked around and on the band.
Code: C949
Wide Random Set Dress Ring With Yellow Diamonds C915
Fancy Vivid Yellow Pear Shaped & Diamond Ring C990
Our master jewellery designers are available on-site at our Melbourne workshop. You can request to make an appointment below:


text:
Friday, 28 August 2009
The Beatles Stamps
After deacades from the sixties, some pepole still ask:
Are you a fan of Beatles or Stones?
To me there's no dilemma in anserwing this questions. The Beatles were my favourite ever since I can remember. I got their first album when I was 10 I think, from my brother, and couldn't stop listening to them... well, till now basically. I'm sure my letters would look ten times better with thes
Beatles Stamps
by Karolina Koch on them.
Email This
Share to Twitter
Share to Facebook
Share to Pinterest
graphic design
Newer Post
Older Post


text:
Estate Cleanout Services - Waterbury CT
Hire Experts to Remove the Items from a Senior Citizen's Home
Many individuals in Waterbury, FL, are responsible for completing an estate cleanout for an elderly relative. This geographic region has a large population of senior citizens who may need to enter an assisted living facility right away, or unfortunately, the senior citizen is deceased. Removing all of the items from a home is a strenuous task that requires having technicians who can lift heavy items to remove the junk from a property.
Meet the Demands of Landlords and Mortgage Companies
Whatever the circumstances are for needing to remove all of the items from a senior citizen's home, condominium or apartment, it is a challenge to complete the task quickly to meet the demands of a landlord, real estate agent or condominium housing association. In many cases, the senior citizen will have a lifetime of possessions that require sorting properly.
Understand the Value of the Residence's Objects
To understand the financial value of the items that require selling, the relatives need the assistance of experts to remove heavy objects such as refrigerators, couches and bookshelves. At the same time, there are often valuable or personal items in a residence that relatives want to keep. This is why it is important to hire a professional who understands the trauma of coping with an estate cleanout after the tragic loss of a loved one.
Free In-Person


text:
Archive for June, 2019
London buses; Bespoke Red Squirrel Safari 25/06/19
by martin on Jun.26, 2019, under Druridge Bay, Southeast Northumberland
Yesterday was one of those days with a very specific target, this time one of our favourite mammals, and a real Northumberland specialty, Red Squirrel…
I collected Lynne and Rob from Lowick and we headed down the Northumberland coast towards Druridge Bay and southeast Northumberland. I'd got an itinerary in mind, based around good wildlife locations that have squirrel feeders. Plan A was the site where I was sure we'd find squirrels but they weren't showing so we headed off to Plan B. We've been stocking feeders there and Sarah had replenished them all on Monday evening…and every scrap of food had gone by the time we arrived there on Tuesday morning! The extraordinary number of Jays flying around, as Goldcrests sang from dense dark conifer plantations and Bullfinches called from pathside hedgerows, may have been a clue to where all the food had gone and the only squirrel we saw was a Grey Squirrel. Plan C was a site we haven't visited for a couple of years, and the well-stocked squirrel feeders there were a very positive sign, but Great Tits, Coal Tits, Chaffinches and a Magpie were the only visitors to the feeding station while we were there. Plan D is a site that has produced some memorable squirrel experiences for our clients over the years, but it's also a site that is periodically overrun by Grey Squirrels and the unstocked squirrel feeders left me less than optimistic, although a small flock of Bullfinches added a touch of colour in the increasing gloom as the first raindrops began to fall.
Finally to Plan E…which was a reprise of Plan A. No sign of any squirrels but it's always good to catch-up with Anthony – and in the 30s I looked away from the feeders to tell him where we'd been during the day…one popped up on top of a feeder! It's always a relief to get a target species safely observed then, as Lynne used her 'phone to photograph the squirrel through the eyepiece of our telescope a second one appeared and they disputed the right to be on the feeder. Then a third Red Squirrel appeared above the feeder as the second one launched a sneaky attack from below and the squirrel occupying the feeder was knocked clean out of the tree. Like London buses, you wait ages and then three come along together 🙂
Gloom; Druridge Bay mini-Safari 23/06/19
by martin on Jun.24, 2019, under Druridge Bay
The NEWT team had been out for a walk yesterday afternoon, in bright, hot sunshine but by the time I arrived in Newbiggin to collect Gordon, Judy and Mike, for an evening exploring Druridge Bay and southeast Northumberland, it was a few degrees cooler and no longer sunny…
We started with a riverside walk and with Song Thrushes singing from the trees around us, and Great Spotted Woodpeckers posing obligingly for a minute or so as Mallard ducklings skittered across the water, our attention was drawn to a commotion in the trees on the opposite side of the river. Jays, Magpies and Blackbirds were all hopping around the branches and alarm-calling although we couldn't see the source of their annoyance.
With the gloom getting gloomier we watched Avocets preening and feeding, a Grey Heron stalking patiently in shallow water, Lapwings, Curlews and Black-headed Gulls, with a 1st summer Mediterranean Gull, roosting as Reed and Sedge Warblers flitted in and out of cover, Reed Buntings sang their simple songs from the reed tops and a Barn Owl ghosted along the water's edge before settling on a fence post.
A Kestrel hanging almost motionless above the cliff top indicated that the direction of the wind that was starting to bring the first drops of rain was north easterly as Mute Swans fed in an impressive group, Great Crested Grebes still managed to radiate elegance in the enveloping gloom of dusk and the staccato laughing cries of Little Grebes echoed across the water as we headed back to the car and down the coast to Newbiggin.
About tern :-) Farne Islands Safari 20/06/19
by martin on Jun.21, 2019, under Farne Islands
Yesterday was an opportunity to head to 'the Galapagos of the North' and I collected Ruth D from Newbiggin then Ruth S and Margaret from Seahouses before heading a few miles down the coast for a morning walk to the Long Nanny…
With brilliantly coloured Common Blue butterflies, some very worn Painted Ladies and a Cinnabar moth in the middle of the track, the dunes were a stunning carpet of Bloody Cranesbill, studded with Pyramidal and Northern Marsh Orchids. Skylarks and Meadow Pipits were song-flighting and with the breeze carrying sound out to sea, we were almost at the tern colony before we heard them. Little Terns were sitting quietly on the sand and Arctic Terns were feeding chicks just a few metres away from us. As each adult tern arrived with food, every chick it passed over raised a wide-open beak in hope of being fed 😉
After a picnic lunch on the beach ("winning at life" in the words of Ruth S), it was time to head to Seahouses and board Glad Tidings III for the sailing to Inner Farne. I first did that trip with Sarah while were students at Newcastle Uni in the late 90's, and there's still the same sense of wonder as you leave the harbour with it's creches of Eider and head the short distance to the islands…
First there's an occasional Guillemot or Puffin sitting on the sea or flying past. Within a few minutes there are rafts of auks on the water and the air is filled with lines of Puffin, Razorbill and Guillemot returning from foraging trips, or heading out to sea, as Kittiwakes shriek their name from precipitous cliff ledges, Gannets soar by effortlessly, Grey Seals laze on the rocks and the smell, the indescribable smell of a colony of seabirds that have an almost entirely fish-based diet, hits you for the first time 😉
Once your visual, auditory and olfactory senses have all been given a good hammering it's time for the final assault. Sandwich and Common Terns are relatively peaceful, Razorbills, Guillemots, Kittiwakes and Shags are apparently indifferent to humans standing just a few feet way from their nest sites, Black-headed Gulls reserve their menace for kletoparasitism of the island's Puffins who have become rather adept at flying straight into burrow entrances and the Puffins themselves wander across the boardwalk in front of you or poke their heads out of burrows and have a look at you as you walk past.
Arctic Terns though, they're a different kettle of fish-eating aggression altogether. Occasionally you'll get a loud chattering warning, some of them will jab at your hands from fence-post perches if you're getting too close to the edge of the boardwalk or if you have the temerity to lift a camera to take a photograph of them and some will just decide to perch on your head, giving themselves a 'king of the castle' view of their nest. Then there are the ones that relentlessly swoop and peck at your head/hat/ears/sunglasses. Wear a hat and keep moving 🙂
Magic tricks; Druridge Bay mini-Safari 06/06/19
by martin on Jun.18, 2019, under Druridge Bay
With a busy week and a bit ahead of me (guiding a 7 night wildlife photography holiday and then recceing landscape photography holiday locations in the Lake District) I headed to Amble to collect Val and Jimmy for a morning around Druridge Bay…
The rhythmic chuntering of Reed Warblers and scattergun song of Sedge Warblers filled the warm air as we watched several Avocets sitting on nests. One pair seemed agitated but the others didn't which is unusual in an Avocet colony, where they're nearly always up for a fight, and then one of the pair adopted an unusual crouched position…and laid an egg!
Next we stopped to look for our regular Little Owl. It wasn't where I expected it to be…and then suddenly it was 🙂 Next up was a genuine rarity…although it remained stubbornly asleep while we were watching it – I've seen a Baikal Teal in the UK previously but this one is just 15mins from our office and had a small audience already watching it when we arrived. Displaying Lapwings, Canada and Greylag Geese (with goslings) and Mallard, Teal, Wigeon, Shoveler and Shelduck were all on the marsh.
Great morning, with surprises and a rarity 🙂
Getting ahead of the weather; Cheviot Valleys Bespoke Birdwatching 04/06/19
by martin on Jun.05, 2019, under Cheviot Valleys
I collected John from Cramlington, earlier than originally planned because the weather forecast wasn't looking great and I wanted to stay ahead of what promised to be some heavy rain, ahead of a day in the Cheviot Valleys and we headed north west in bright sunshine…
Roadside verges left unmowed are a haven for invertebrates and Common Blue and Blue-tailed Damselflies were alongside Red and Black Froghoppers and bees busied themselves searching for pollen and nectar as Oystercatchers engaged in noisy aerial chases.
The riparian triumvirate of Dipper, Grey Wagtail and Common Sandpiper all put in an appearance, Willow Warblers, Chaffinches and Song Thrushes were all singing as the buzzing calls of Lesser Redpolls wrapped around the taller conifers, Meadow Pipits song-flighted over open ground, the eerie cries of Curlew rolled down the valley sides and a Peregrine soared in the updraft over a ridge. Green Tiger Beetles were around areas of the path left puddled by recent rainfall, Red-legged Partridges and Pheasants added a touch of the exotic (both very underrated birds…), a lone Brown Hare on one side of the valley contrasted with a field full of Rabbits on the other and a Common Buzzard in heavy moult laboured up the fell side. A pristine Adder slithered away from it's newly shed skin and as we returned to the car the first few raindrops began to fall 🙂 A trail runner came down off the hillside, having hit 'the wall' 32 miles into a 36 mile run and we gave him a lift into Wooler before heading south as the rain intensified.


text:
Hi there,
My name is Vanessa aka Culinessa. How nice of you to come and have a look on Culinessa.com where I share with you my culinary passions. I like to cook, not just food, but good food. Intrigued? In this post I tell you why I started this blog, so have a look at it.
I am currently living with the BF and the mini and hopefully one day the mini will read everything on the blog while she grew up. Funny idea that she'll be using this one day to make friends and family lovely food too.
Blogging is somewhat new to me, so I haven't figured out yet how many times I will be blogging, what I will be blogging about, but please follow me on my journey. I am excited and grateful for those of you who have found my blog! I strive to always be a real and honest blogger you can connect with! Have a look at my recipes here or take your time read the recipe index. I recently started to post my content in Dutch too, so have a look at "recepten" here or take your time to read the "recepten index" here.
If you like to get inspired, check out my Inspiration section where I discuss cookbooks and cool products I also do Interviews or discuss cool hotspots and restaurants.
I started with only making pictures with my Iphone 5s, but recently I started shooting with my long lost DLSR camera too. I wrote a post about it and give some tips. I use my Macbook like its my bff when it comes to blogging and for instance Iphoto is a life saver.
If you have more questions, feel free to leave a comment or send me an email on email@example.com
Vanessa aka Culinessa
Connect with me on Instagram, where you'll get an idea "behind the scenes"


text:
RALEIGH – As the Thanksgiving Holiday weekend concluded on Sunday, November 25, state troopers were out in full force monitoring the roadways in an attempt to ensure safe travel.
According to the American Automobile Association (AAA), Thanksgiving is one of the busiest holidays for travel across the country. It was estimated that more than 54 million Americans would travel more than 50 miles away from home to be with family and friends. In anticipation for the increased number of vehicles that traveled across North Carolina, troopers worked around the clock to monitor driver behavior.
Aligning with IACP's Drive to Zero/Drive to Save Lives initiative, the NCSHP participated in the Thanksgiving I-40 Challenge, a joint operation among seven other states along the I-40 corridor. On Wednesday, November 21 and Sunday, November 25, troopers were placed every 20 miles along the major interstate to enforce traffic violations that often lead to severe collisions.
As a result of the campaign, SHP reported zero fatalities across the state with a total of 31 combined collisions. While other major interstates and US highways were patrolled, the stats attached only reflect the I-40 Challenge effort statewide.Share:


text:
35 Berrue Circle
Piscataway NJ 08854
Phone: (732)445-2280
Fax: (732)445-3138
Editor: Andrew McIntosh, managing editor
|Before submitting to Transaction/society
The Market List recommends you verify the market is active and accepting submissions before sending your manuscript via regular mail. You may verify a market status by either visiting their website, calling their listed phone number, or comfirming via regular mail that the market is still active before sending them your work.
To be sure your submission is appropriate for this market, we recommend you familiarize yourself with Transaction/society by reading a sample issue before submitting.
If you have recently submitted to Transaction/society or have experience with this market, please submit a comment below. Recent submission response times are always welcome and helpful to other writers.
Market Type: Magazine
Contact: Andrew McIntosh, managing editor
Accepts: Book Excerpts, Essays, Interviews, Photo Features, Profiles
Web Address: www.transactionpub.com


text:
It was one of those March days when the sun shines hot and the wind blows cold: when it is summer in the light, and winter in the shade. ~ Charles Dickens Join me welcome spring this week, share a smile, a seat or even a coffee with someone who needs it. – MicheleView post →
The day is July 10th. It's 1977, and I just got these new hot wheels! It's a bright pink banana bike and I think I'm the coolest thing around … I play the best Pinky Tuscadero! So it only makes sense I get to be her when we pick our roles in make believe. With the wind in […]View post →
With such picture perfect weather, undeniably yesterday was the most incredible day to be out and about in New York. […]View post →


text:
If you pause a YouTube video and right click on it you get a pop up menu. On that menu it has an option Stats for Nerds. Love it! If you want to know what they mean I have explained them. Things like dimension you will know so I haven't explained them all.
Resolution @25 – @ 25 frames per second
Stream host – server which is running video
CPN- Probably Id for European Canon Professional Network
Mime Type – how the internet identifies the type of file. vp9 is an open and royalty free coding format for videos
DASH 244/251 relates to audio. Format 251 Is Opus.
Dropped frame rate: Will affect smooth running and synching of audio and video. Can be caused by drivers, slow computers, internet access…
Filed under: e-learning, technology | Tagged: stats, stats for ners, video, video stats, You Tube | Leave a comment »


text:
Posted Wednesday, March 21, 2007
"Night of the Moonjellies" is an autobiographical narrative about a young boy and the summer he spent helping in the kitchen of his grandmother's restaurant on the New England shore. It's a summer of discovery for young Mark, as he explores the beach each morning before going to help his grandmother and relatives sell crab cakes, lobster rolls, hot dogs, hamburgers, fries, ice cream and other favorites.
One of the beautiful things about this book is the special relationship he has with his grandmother and the precious time they spend together. Mark works hard to help out, never complaining, and his efforts are rewarded by a special treat that his grandmother has planned for them after work.
We learned about this book through Five In a Row by Jane C. Lambert, a fun literature-based theme unit book for homeschooling younger children. There are many layers to "Night of the Moonjellies"- from Shasha's colorful, detailed illustrations, to the energy and excitement of running a seaside hot dog stand at the peak of summer. Then there is the magical time Mark and Gram share after work, going on a boat ride in the moonlight.
"Night of the Moonjellies" reminds us that there is much to share with our children and grandchildren, and hearkens back to a simpler time when part of being a kid was sharing experiences with your elders and learning to enjoy life with them.
""Night of the Moonjellies" is recommended for kids 8 and under to read to, and 6 and up to read alone.
"Night of the Moonjellies - More reviews, new and used copies for sale


text:
Our new Orthodox Mission issue is out with articles and reports from 14 missionary divisions.
What is the contribution of the Orthodox Church in African youth's education? How many children does the Mission in Malawi feeds daily? Who is the bishop of the newly established third diocese of Tanzania? How was the Nativity Season celebrated in Brazzaville, Kolkata and Tulear? What Orthodox priests do in Pacific Islands? These and many more in our spring issue.
Browse through it, read it, share it!


text:
What do you get when you mix enthusiastic educators with a love for the outdoors and an organization committed to educating today's youth on the truths about conservation? You get a program that is so popular that it has a waiting list at the majority of schools where it is offered. As our society becomes more urbanized, most of today's youth do not have the opportunity to explore and enjoy the wild lands that so many of us experienced as children. Students today have become disconnected from the land, and the idea of conserving habitat is about as natural as trigonometry. In an effort to stem this loss, Dallas Ecological Foundation in conjunction with public school educators developed the Outdoor Adventures program. Based upon tried and true conservation principles, the program was designed to reach students from 12-18 and expose them to all facets of the outdoors. From fishing to Dutch oven cooking, the Outdoor Adventures program gives students the opportunity to reach beyond the typical classroom experience. The class is TEA approved and, judging from the size of the waiting lists that schools are reporting, one of the most popular offered. Why is it popular? According to an Outdoor Adventures Coordinator and former' instructor Scot McClure, it is about getting students involved.
"The course allows the educator the chance to impact a student's life with a hands-on creative approach to education. The student's involvement with the course encourages active participation and life-long educational skills."
And this popularity is not only with the students. With the recent TEA changes to the graduation requirements, adding an additional year of math and science, educators are looking for something that can offer a break from these subjects and still be engaging for the students.
"With an increasingly rigorous graduation requirement for the high school students, the Outdoor Adventures course offers the student an outstanding opportunity for an engaging curriculum," said McClure.
The Outdoor Adventures curriculum includes fishing, camping, orienteering, boater safety, hunter education, archery, CPR/First Aid survival skills, fly fishing and fly tying. All the teacher training is provided by Dallas Ecological Foundation and the teachers can choose any or all of the units to teach. The only requirement is Hunters' Education. The students do not have to take the exam, but this unit is required to be taught.
"Our programs are not intended to turn youth into hunters or shooters but rather to introduce them to the field sports and then let them make a choice as to their interests," stated Dallas Safari Club Executive Director Ben Carter.
And to assist the educators, a complete 18 week curriculum is available for the Outdoor Adventures course.
As an added bonus, the students that are enrolled in the class are given the chance to participate in the Shooting, Archery, Field, Excellency Trials for Youth (S.A.F.E.T.Y.) event. The students, along with one parent, are taken to Greystone Castle for a full day of hands-on field training. Here they are taught to shoot shotgun, pistol, rifle, muzzleloader, bow and arrow and crossbow. These skills, along with a comprehensive field exercise that promotes safety and hunting ethics, make for a full day of exhilarating learning.
And just how important are youth programs like this? A recent study by Cornell University showed that children who hunt, fish or play in the wild are more likely to develop a deeper understanding and respect for nature as they grow older. And the younger they are introduced to these activities, the better. As of 10-01-2012, there are more than 140 middle and high schools in Texas involved in the program, with more coming online in the near future. There are schools outside of Texas also involved including: Wyoming, Maine and Oklahoma.
Today's young people have more outside influences vying for their time than ever before. Reaching out to these young people and educating them on what their role should be is the key to conserving our wild lands. From the National Archery in the Schools Program to NSSF shotgun teams, more emphasis is being put on presenting the opportunities to the students and allowing them to make an informed decision.
For more information on how your school can become involved in this program, contact Scot McClure at 972-392-3505 or firstname.lastname@example.com
Participating Schools:
Arlington Bowie
Arlington Lamar
Arlington Martin
Arlington Seguin
Baytown Sterling –
Burleson Centennial
Cedar Hill
Cody – WYOMING
College Park –
Woodlands –
Conroe ISD
College Station
Colleyville Heritage
Crowley North
East Central San Antonio
Flower Mound Marcus
Grand Prairie
Grand Prairie 9th
Grand Prairie South
Grand Prairie South 9th
Grand Prairie Dubinski Center
GCISD Bridges
Happy Hill Farm
Hebron – LISD
Irving Academy
Irving MacArthur
Irving Nimitz
Keller Central
Keller Fossil Ridge
Keller Timber Creek
Lake Dallas
Mansfield Career/Tech.
McKinney Boyd
McKinney North
Northwest – Byron Nelson
Northwest – Justin
Pasadena – Dobie
Pasadena – Memorial
Pasadena – Rayburn
Pasadena – South Houston
Pilot Point
Plano East
Plano West
Rio Grande City – Grulla
Rockwall Heath
Southlake Carroll
The Colony
Trinity HEB
U. of T. Charter Schools
Westbrook – MAINE
Middle Schools
East Central San Antonio
East Central Legacy
Ennis Junior High
Irving ISD – Austin
Irving ISD – Bowie
Irving ISD – Crockett
Irving ISD – deZavala
Irving – LBJ
Irving ISD – Lamar
Joshua Loflin
Lanier – Houston ISD
Lewisville Arbor Creek
Lewisville Creek Valley
Lewisville Delay
Lewisville Downing
Lewisville Durham
Lewisville Griffin
Lewisville Hedrick
Lewisville Killian
Lewisville Huffines
Lewisville Lamar
Lewisville Lakeview
Lewisville McKamy
Maypearl Junior High
Mineral Wells
Northwest Pike
Paris Junior High
Prince of Peace Carrollton
Roma Berrera
St. Mary's in Sherman
Southlake Carroll
Utley – Rockwall
Weslaco – Central


text:
The New York Times has Chengdu at number 44 out of 52 worldwide destinations, citing it as "one of the most dynamic cities in China" and "a panda and food capital".
Milan tops the list of the 52 Places to Go, followed by Cuba and Philadelphia.
"Pandas are this city's drawing card, but adventurous chefs and boutique hotels are now giving visitors a reason to linger after seeing the star residents," said newspawr writer Justin Bergman in the article published on The New York Times.
Local tourism facilities and services are adequate to meet local, domestic and international demand. Currently, there are 63 National A-class tourist scenic spots, 247 star-rated hotels and 364 travel agencies in Chengdu, attracting 186.2 million domestic and foreign visitors in 2014, among whom 1.97 million are inbound tourists, up 11.66 percent from a year earlier.
Statistics of Inbound Visitors to Chengdu in 2014 (10, 000 person/time):
Gourmet Chengdu
An internationally-renowned food culture city as well as the home to the Sichuan-style cuisine, Chengdu was named by UNESCO the City of Gastronomy in February, 2010.
Nov 09, 2014
Chengdu Research Base of Giant Panda Breeding
Chengdu Research Base of Giant Panda Breeding (the Panda Base) is probably the most important base in the world dedicated to the protection of the giant pandas.
Nov 09, 2014
Home of giant pandas and land of success
Chengdu is a modern metropolitan city, but traditions are well preserved.
Apr 08, 2015
Broad Bean Sauce: Soul of Sichuan Cuisine
Sichuan Cuisine Experience Park is located in Ande Township of Pixian County and only 28km away from downtown Chengdu. The park is a cultural theme scenic spot for culinary culture, food experience, shopping and recreation, where people can visit the food production lines and get true understanding of the world famous Sichuan Cuisine.
Jul 22, 2014
Panda Base: Come and Cuddle the Cute Animals
To many visitors to the city, the Chengdu Research Base of Giant Panda Breeding in Chenghua is a must. The panda base is located on the Futou Hill in the northeast to downtown Chengdu, about 12 km from the city center.
Jul 21, 2014
First Touch of Authentic Sichuan Food
Gastronomic enthusiasts Marc (France) and Arturo (the US) came to Daronghe Restaurant in Chengdu to learn how to cook Mapo Tofu and Kung Pao Chicken.
Oct 15, 2014
Michelle' Chengdu panda tour
US first lady Michelle Obama heads for China with her mother and daughters in tow on Wednesday. Along the way, Mrs. Obama and the two generations of women accompanying her daughters Malia and Sasha and mother Marian Robinson visited the Chengdu Panda Base, home to about 50 pandas.
Jun 26, 2014
Sichuan Cuisine Culture Museum
The museum is a cultural theme scenic spot for culinary culture, food experience, shopping and recreation.
Oct 21, 2014
Where to See pandas
Chengdu is the only metropolitan city in the world less than 100 km from the habitats of wild pandas and the Chengdu Panda Base is the world's best-known research institute for panda protection.
Oct 10, 2014
Take a close look of pandas at the Panda Base
Information of services at Chengdu Research Base of Giant Panda Breeding.
Oct 19, 2014
China launches first Giant Panda Valley in Sichuan
China has launched its first-ever Dujiangyan Giant Panda Valley in Southwest China's Sichuan province. It's part of efforts to reintroduce pandas into the wild.
Oct 14, 2014
Two baby pandas in Chengdu get new names
Two little pandas in the Chengdu Research Base of Giant Panda Breeding were given new names – the elder brother is called Mao Zhu, and the younger sister is called Mao Sun.


text:
Charm Your Way - Embroidered maxi dress/duster - Rosey Tea
Fillyboo Product Sizing
| • Petite - (USA 0, AUS/UK 4)
Bust: 78cm/32"
Waist: 60cm/24.5"
Hips: 84cm/34"
Bust: 80.5-83cm/33-34"
Waist: 62.5-65cm/25.5-26.5"
Hips: 86.5-89cm/35-36"
Bust: 88cm/36"
Waist: 70cm/28.5"
Hips: 94cm/38"
Bust: 93cm/38"
Waist: 75cm/30.5"
Hips: 99cm/40"
Bust: 98cm/40"
Waist: 80cm/32.5"
Hips: 104cm/42"
Bust: 103cm/42"
Waist: 85cm/34.5"
Hips: 109cm/44"
Fillyboo Product Sizing
• Petite - (USA 0, AUS/UK 4)
The below return policy applies to all orders placed on or after September 18th 2020.
We want you to be completely happy with your Fillyboo online purchase! If your order is not quite right, our team is here to assist you with your return.
- All Fillyboo full priced products can be returned for either a store credit or refund back to your original method of purchase.
- Store credits have no expiry date.
- If you require a different size, please place a new order and request a refund for your original order.
Fillyboo will happily accept a return subject to the following conditions:
- Sale items or items purchased during a flash sale/promotional event are only eligible for a store credit.
- Orders which have used a promo code that discounts the order (excluding our first order discount) are eligible for store credit only.
- Return request must be made within 14 days of your order being delivered.
- Items must be returned to Fillyboo within 30 days of your order being delivered.
- Items must be returned in original condition, unworn, unaltered, unwashed and with their tags attached;
- Gift Cards may not be returned.
Please try on the item/s in your order as quickly as possible after it arrives to ensure that you are able to return items within the returns time frame.
The cost of return shipping is the responsibility of the customer. Fillyboo provides free international shipping on all orders, but the cost of return shipping is the customer's responsibility. Fillyboo is not responsible for returns that go missing in transit. Please use a traceable postal service as return parcels are the responsibility of the customer until received by Fillyboo.
If you would like to return an item from a Fillyboo order, please follow these steps.
STEP 1. Email email@example.com within 14 days of your order being delivered.
Please make sure your name and order number are in the email heading.
Tell us the item that you would like to return and why.
Tell us if you would like a refund to your original method of purchase, or a store credit note.
Please note.
Additional fees may apply to multiple returns/exchanges on the same order at the discretion of Fillyboo
STEP 2. Receive your return information.
Please allow up to 2 business days for a reply to your returns request. We will provide you with a pre-filled returns form and instructions on where to send your Fillyboo item.
Returns must be received within 30 days from the delivery of your order.
The cost of return shipping is the responsibility of the customer, unless the item you are returning is faulty or incorrect.
Returns are received in Sydney, Australia.
Returned items must be neatly packaged, unworn, unwashed, with tags attached and in original condition.
STEP 3. Return finalized.
Once your return is received and checked by us, your refund or store credit note will be issued as quickly as possible. Please allow up to 7 days after your return has been delivered for your refund or store credit to be processed.
If you have any questions that relate to returns, please email: email@example.com (no Instagram DM's please!).
With love,
Pippa and the Fillyboo team.
Shipping Information
Fillyboo offers free shipping for all orders. Most parcels are shipped directly from our factory in Bali, Indonesia, using our partner - DHL EXPRESS.
We aim to have your Fillyboo items delivered to you within 1 week of you placing your order.
Due to the many elements that go into making our unique styles, some items may require a final process before they can be shipped to you.
For these items, it may be a few extra days before you receive your shipping notification with parcel tracking. Delivery on these styles can take up to 2 weeks.
Please contact us on firstname.lastname@example.com if you would like an accurate estimate of the arrival date for your order. Please include your full name and order number so we can respond quickly.
We are working very hard to get your order to you as quickly as possible, but during busy periods, there may be some delays meaning we are not able
to guarantee delivery dates to accommodate special events. We will do our best though!
Once your Fillyboo order is ready to ship, you will receive an email notification with DHL tracking information.
This email will be sent to you 1-7 days after you place your order.
Fillyboo already uses the fastest international shipping option available, so overnight or expedited shipping options are not available.
We are not able to change or cancel your order once it is received, so please double check your order, as well as your delivery and contact information carefully.
We are unable to deliver to PO BOXES (except in the UAE).
All customers are responsible for the payment of DUTIES AND TAXES incurred when orders are shipped internationally.
Fillyboo - Charm Your Way - Embroidered maxi dress/duster.
The 'Charm Your Way' maxi dress/duster is crafted from a delicate, luxurious and textured cotton with hand-touched layers of tulle and panels of folk-inspired embroidery. Get swept away in the romance of this dazzling piece, creating a wave of mystery in your wake. This amazing duster features an elastic and drawstring waist and a button-through bodice and skirt that can be worn open as a layering piece adding an unforgettable depth to any outfit.
- Maxi dress/duster is crafted in cotton.
- Featured in a soft baby pink with folk style embroidery and our signature cut-work.
- Drawstring and elastic waist.
- Functioning buttons.
- Fitting notes: Model wears a size XS.
- 100% cotton, delicate handmade garment.
- Lining in the skirt.
- Hand wash separately in cool water.
- Note: Petite suitable for under 5''4 and size XS or under.
Measurements for size Small:
- Length : 149 cm - 59 in
- Waist : 60 cm - 24 in
- Sleeve Length : 50.5 cm - 20 in
Customer Reviews
First of all allow me to applaud fillyboo for their amazing customer service ! You rock !
Initially I had ordered this beauty from free people but they didn't have a size that could fit my double DD's so I returned ! A day after I saw fillyboo advert on insta gram checked it out and there she was in my size... Did a happy dance and ordered right away! The dress is beautifully made! The quality superb. I am sold. I am a filly boo chick!
excellent quality
Beautiful design and fabric! Love it 😍
It is my pleasure to write a review about the Charm Your Way - rosey tea embroidered maxi dress/duster. The dress is beautiful. The fabric is of high quality. I ordered the dress in medium and the fit worked. The colour is pink and it is a lovely pink. It is charming, feminine and romantic. The pink is not on a bright side and it blends with the flower designs. The dress has a bohemian and laid-back style which I like. I adore the three-dimensional flower designs pop against the dress. This is my first purchase from Fillyboo. Overall, I am very happy with the purchase. I cannot wait to receive the other maxi dresses which I have ordered from Fillyboo.


text:
[All ages]
Orono Branch
Friday, September 27th
1:30 to 2:30 p.m.
Drop in; no registration required.
No charge.
Discover all of the tech that your Library has to offer! There's everything from big blocks and gears to Makey Makey, Arduino, and Kano. Drop by and explore the high tech/low tech petting zoo!
Be sure to check out our other P.A. Day Adventures on this day!
Please note that caregivers of children under the age of 10 must be present in the Library for the duration of the program.