Profiling of short RNAs during fleshy fruit development reveals stage-specific sRNAome expression patterns.

Plants feature a particularly diverse population of short (s)RNAs, the central component of all RNA silencing pathways. Next generation sequencing techniques enable deeper insights into this complex and highly conserved mechanism and allow identification and quantification of sRNAs. We employed deep sequencing to monitor the sRNAome of developing tomato fruits covering the period between closed flowers and ripened fruits by profiling sRNAs at 10 time-points. It is known that microRNAs (miRNAs) play an important role in development but very little information is available about the majority of sRNAs that are not miRNAs. Here we show distinctive patterns of sRNA expression that often coincide with stages of the developmental process such as flowering, early and late fruit maturation. Moreover, thousands of non-miRNA sRNAs are differentially expressed during fruit development and ripening. Some of these differentially expressed sRNAs derived from transposons but many derive from protein coding genes or regions that show homology to protein coding genes, several of which are known to play a role in flower and fruit development. These findings raise the possibility of a regulative role of these sRNAs during fruit onset and maturation in a crop species. We also identified six new miRNAs and experimentally validated two target mRNAs. These two mRNAs are targeted by the same miRNA but do not belong to the same gene family, which is rare for plant miRNAs. Expression pattern and putative function of these targets indicate a possible role in glutamate accumulation, which contributes to establishing the taste of the fruit.

The genetic diversity and evolution of field pea (Pisum) studied by high throughput retrotransposon based insertion polymorphism (RBIP) marker analysis.

BACKGROUND: The genetic diversity of crop species is the result of natural selection on the wild progenitor and human intervention by ancient and modern farmers and breeders. The genomes of modern cultivars, old cultivated landraces, ecotypes and wild relatives reflect the effects of these forces and provide insights into germplasm structural diversity, the geographical dimension to species diversity and the process of domestication of wild organisms. This issue is also of great practical importance for crop improvement because wild germplasm represents a rich potential source of useful under-exploited alleles or allele combinations. The aim of the present study was to analyse a major Pisum germplasm collection to gain a broad understanding of the diversity and evolution of Pisum and provide a new rational framework for designing germplasm core collections of the genus. RESULTS: 3020 Pisum germplasm samples from the John Innes Pisum germplasm collection were genotyped for 45 retrotransposon based insertion polymorphism (RBIP) markers by the Tagged Array Marker (TAM) method. The data set was stored in a purpose-built Germinate relational database and analysed by both principal coordinate analysis and a nested application of the Structure program which yielded substantially similar but complementary views of the diversity of the genus Pisum. Structure revealed three Groups (1-3) corresponding approximately to landrace, cultivar and wild Pisum respectively, which were resolved by nested Structure analysis into 14 Sub-Groups, many of which correlate with taxonomic sub-divisions of Pisum, domestication related phenotypic traits and/or restricted geographical locations. Genetic distances calculated between these Sub-Groups are broadly supported by principal coordinate analysis and these, together with the trait and geographical data, were used to infer a detailed model for the domestication of Pisum. CONCLUSIONS: These data provide a clear picture of the major distinct gene pools into which the genus Pisum is partitioned and their geographical distribution. The data strongly support the model of independent domestications for P. sativum ssp abyssinicum and P. sativum. The relationships between these two cultivated germplasms and the various sub-divisions of wild Pisum have been clarified and the most likely ancestral wild gene pools for domesticated P. sativum identified. Lastly, this study provides a framework for defining global Pisum germplasm which will be useful for designing core collections.

High-throughput sequencing of Medicago truncatula short RNAs identifies eight new miRNA families.

BACKGROUND: High-throughput sequencing technology is capable to identify novel short RNAs in plant species. We used Solexa sequencing to find new microRNAs in one of the model legume species, barrel medic (Medicago truncatula). RESULTS: 3,948,871 reads were obtained from two separate short RNA libraries generated from total RNA extracted from M. truncatula leaves, representing 1,563,959 distinct sequences. 2,168,937 reads were mapped to the available M. truncatula genome corresponding to 619,175 distinct sequences. 174,504 reads representing 25 conserved miRNA families showed perfect matches to known miRNAs. We also identified 26 novel miRNA candidates that were potentially generated from 32 loci. Nine of these loci produced eight distinct sequences, for which the miRNA* sequences were also sequenced. These sequences were not described in other plant species and accumulation of these eight novel miRNAs was confirmed by Northern blot analysis. Potential target genes were predicted for most conserved and novel miRNAs. CONCLUSION: Deep sequencing of short RNAs from M. truncatula leaves identified eight new miRNAs indicating that specific miRNAs exist in legume species.

Gene-based sequence diversity analysis of field pea (Pisum).

Sequence diversity of 39 dispersed gene loci was analyzed in 48 diverse individuals representative of the genus Pisum. The different genes show large variation in diversity parameters, suggesting widely differing levels of selection and a high overall diversity level for the species. The data set yields a genetic diversity tree whose deep branches, involving wild samples, are preserved in a tree derived from a polymorphic retrotransposon insertions in an identical sample set. Thus, gene regions and intergenic "junk DNA" share a consistent picture for the genomic diversity of Pisum, despite low linkage disequilibrium in wild and landrace germplasm, which might be expected to allow independent evolution of these very different DNA classes. Additional lines of evidence indicate that recombination has shuffled gene haplotypes efficiently within Pisum, despite its high level of inbreeding and widespread geographic distribution. Trees derived from individual gene loci show marked differences from each other, and genetic distance values between sample pairs show high standard deviations. Sequence mosaic analysis of aligned sequences identifies nine loci showing evidence for intragenic recombination. Lastly, phylogenetic network analysis confirms the non-treelike structure of Pisum diversity and indicates the major germplasm classes involved. Overall, these data emphasize the artificiality of simple tree structures for representing genomic sequence variation within Pisum and emphasize the need for fine structure haplotype analysis to accurately define the genetic structure of the species.

Variation in Splicing Efficiency Underlies Morphological Evolution in Capsella.

Understanding the molecular basis of morphological change remains a central challenge in evolutionary-developmental biology. The transition from outbreeding to selfing is often associated with a dramatic reduction in reproductive structures and functions, such as the loss of attractive pheromones in hermaphroditic Caenorhabditis elegans and a reduced flower size in plants. Here, we demonstrate that variation in the level of the brassinosteroid-biosynthesis enzyme CYP724A1 contributes to the reduced flower size of selfing Capsella rubella compared with its outbreeding ancestor Capsella grandiflora. The primary transcript of the C. rubella allele is spliced more efficiently than that of C. grandiflora, resulting in higher brassinosteroid levels. These restrict organ growth by limiting cell proliferation. More efficient splicing of the C. rubella allele results from two de novo mutations in the selfing lineage. Thus, our results highlight the potentially widespread importance of differential splicing efficiency and higher-than-optimal hormone levels in generating phenotypic variation.

Insertional polymorphism and antiquity of PDR1 retrotransposon insertions in pisum species.

Sequences flanking 73 insertions of the retrotransposon PDR1 have been characterized, together with an additional 270 flanking regions from one side alone, from a diverse collection of Pisum germ plasm. Most of the identified flanking sequences are repetitious DNAs but more than expected (7%) lie within nuclear gene protein-coding regions. The approximate age of 52 of the PDR1 insertions has been determined by measuring sequence divergence among LTR pairs. These data show that PDR1 transpositions occurred within the last 5 MY, with a peak at 1-2.5 MYA. The insertional polymorphism of 68 insertions has been assessed across 47 selected Pisum accessions, representing the diversity of the genus. None of the insertions are fixed, showing that PDR1 insertions can persist in a polymorphic state for millions of years in Pisum. The insertional polymorphism data have been compared with the age estimations to ask what rules control the proliferation of PDR1 insertions in Pisum. Relatively recent insertions (< approximately 1.5 MYA) tend to be found in small subsets of the Pisum accessions set, "middle-aged" insertions (between approximately 1.5 and 2.5 MYA) vary greatly in their occurrence, and older insertions (> approximately 2.5 MYA) are mostly found in small subsets of Pisum. Finally, the average age estimate for PDR1 insertions, together with an existing data set for PDR1 retrotransposon SSAP markers, has been used to derive an estimate of the effective population size for Pisum of approximately 7.5 x 10(5).

A microarray-based high throughput molecular marker genotyping method: the tagged microarray marker (TAM) approach.

A microarray-based method has been developed for scoring thousands of DNAs for a co-dominant molecular marker on a glass slide. The approach was developed to detect insertional polymorphism of transposons and works well with single nucleotide polymorphism (SNP) markers. Biotin- terminated allele-specific PCR products are spotted unpurified onto streptavidin-coated glass slides and visualised by hybridisation of fluorescent detector oligonucleotides to tags attached to the allele- specific PCR primers. Two tagged primer oligonucleotides are used per locus and each tag is detected by hybridisation to a concatameric DNA probe labelled with multiple fluorochromes.

Deep sequencing of tomato short RNAs identifies microRNAs targeting genes involved in fruit ripening.

In plants there are several classes of 21-24-nt short RNAs that regulate gene expression. The most conserved class is the microRNAs (miRNAs), although some miRNAs are found only in specific species. We used high-throughput pyrosequencing to identify conserved and nonconserved miRNAs and other short RNAs in tomato fruit and leaf. Several conserved miRNAs showed tissue-specific expression, which, combined with target gene validation results, suggests that miRNAs may play a role in fleshy fruit development. We also identified four new nonconserved miRNAs. One of the validated targets of a novel miRNA is a member of the CTR family involved in fruit ripening. However, 62 predicted targets showing near perfect complementarity to potential new miRNAs did not validate experimentally. This suggests that target prediction of plant short RNAs could have a high false-positive rate and must therefore be validated experimentally. We also found short RNAs from a Solanaceae-specific foldback transposon, which showed a miRNA/miRNA*-like distribution, suggesting that this element may function as a miRNA gene progenitor. The other Solanaceae-specific class of short RNA was derived from an endogenous pararetrovirus sequence inserted into the tomato chromosomes. This study opens a new avenue in the field of fleshy fruit biology by raising the possibility that fruit development and ripening may be under miRNA regulation.

The tagged microarray marker (TAM) method for high-throughput detection of single nucleotide and indel polymorphisms.

The tagged microarray marker (TAM) method allows high-throughput differentiation between predicted alternative PCR products. Typically, the method is used as a molecular marker approach to determining the allelic states of single nucleotide polymorphisms (SNPs) or insertion-deletion (indel) alleles at genomic loci in multiple individuals. Biotin-labeled PCR products are spotted, unpurified, onto a streptavidin-coated glass slide and the alternative products are differentiated by hybridization to fluorescent detector oligonucleotides that recognize corresponding allele-specific tags on the PCR primers. The main attractions of this method are its high throughput (thousands of PCRs are analyzed per slide), flexibility of scoring (any combination, from a single marker in thousands of samples to thousands of markers in a single sample, can be analyzed) and flexibility of scale (any experimental scale, from a small lab setting up to a large project). This protocol describes an experiment involving 3,072 PCRs scored on a slide. The whole process from the start of PCR setup to receiving the data spreadsheet takes 2 d.