Sliced microRNA targets and precise loop-first processing of MIR319 hairpins revealed by analysis of the Physcomitrella patens degradome

Expression profiling of the 5′ ends of uncapped mRNAs (“degradome” sequencing) can be used to empirically catalog microRNA (miRNA) targets, to probe patterns of miRNA hairpin processing, to examine mRNA decay, and to analyze accumulation of endogenous short interfering RNA (siRNA) precursors. We sequenced and analyzed the degradome of the moss Physcomitrella patens, an important model system for functional genomic analyses in plant evolution. A total of 52 target mRNAs of 27 different Physcomitrella miRNA families were identified. Many targets of both more conserved and less conserved miRNA families encoded putative regulatory proteins. Remnants of MIRNA hairpin processing also populated the degradome data and indicated an unusual “loop-first” mode of precise processing for the MIR319 gene family. Precise loop-first processing was confirmed for native Physcomitrella, rice, and Arabidopsis MIR319 hairpins, as well as an Arabidopsis artificial MIRNA (aMIRNA) based upon a MIR319 backbone. MIR319 is thus a conserved exception to the general rule of loop-last processing of MIRNA hairpins. Loop-first MIR319 processing may contribute to the high efficacy of a widely used MIR319-based strategy for aMIRNA production in plants.     

Expansion of MIR482/2118 by a class‐II transposable element in cotton

Some plant microRNA (miRNA) families contain multiple members generating identical or highly similar mature miRNA variants. Mechanisms underlying the expansion of miRNA families remain elusive, although tandem and/or segmental duplications have been proposed. In this study of two tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, and their extant diploid progenitors, Gossypium arboreum and Gossypium raimondii, we investigated the gain and loss of members of the miR482/2118 superfamily, which modulates the expression of nucleotide‐binding site leucine‐rich repeat (NBS‐LRR) disease resistance genes. We found significant expansion of MIR482/2118d in G. barbadense, G. hirsutum and G. raimondii, but not in G. arboreum. Several newly expanded MIR482/2118d loci have mutated to produce different miR482/2118 variants with altered target‐gene specificity. Based on detailed analysis of sequences flanking these MIR482/2118 loci, we found that this expansion of MIR482/2118d and its derivatives resulted from an initial capture of an MIR482/2118d by a class‐II DNA transposable element (TE) in G. raimondii prior to the tetraploidization event, followed by transposition to new genomic locations in G. barbadense, G. hirsutum and G. raimondii. The ‘GosTE’ involved in the capture and proliferation of MIR482/2118d and its derivatives belongs to the PIF/Harbinger superfamily, generating a 3‐bp target site duplication upon insertion at new locations. All orthologous MIR482/2118 loci in the two diploids were retained in the two tetraploids, but mutation(s) in miR482/2118 were observed across all four species as well as in different cultivars of both G. barbadense and G. hirsutum, suggesting a dynamic co‐evolution of miR482/2118 and its NBS‐LRR targets. Our results provide fresh insights into the mechanisms contributing to MIRNA proliferation and enrich our knowledge on TEs.     

Arabidopsis lyrata Small RNAs: Transient MIRNA and Small Interfering RNA Loci within the Arabidopsis Genus

Twenty-one-nucleotide microRNAs (miRNAs) and 24-nucleotide Pol IV-dependent small interfering RNAs (p4-siRNAs) are the most abundant types of small RNAs in angiosperms. Some miRNAs are well conserved among different plant lineages, whereas others are less conserved, and it is not clear whether less-conserved miRNAs have the same functionality as the well conserved ones. p4-siRNAs are broadly produced in the Arabidopsis genome, sometimes from active hot spot loci, but it is unknown whether individual p4-siRNA hot spots are retained as hot spots between plant species. In this study, we compare small RNAs in two closely related species (Arabidopsis thaliana and Arabidopsis lyrata) and find that less-conserved miRNAs have high rates of divergence in MIRNA hairpin structures, mature miRNA sequences, and target-complementary sites in the other species. The fidelity of miRNA biogenesis from many less-conserved MIRNA hairpins frequently deteriorates in the sister species relative to the species of first discovery. We also observe that p4-siRNA occupied loci have a slight tendency to be retained as p4-siRNA loci between species, but the most active A. lyrata p4-siRNA hot spots are generally not syntenic to the most active p4-siRNA hot spots of A. thaliana. Altogether, our findings indicate that many MIRNAs and most p4-siRNA hot spots are rapidly changing and evolutionarily transient within the Arabidopsis genus.     

A novel method for the multiplexed target enrichment of MinION next generation sequencing libraries using PCR-generated baits

The enrichment of targeted regions within complex next generation sequencing libraries commonly uses biotinylated baits to capture the desired sequences. This method results in high read coverage over the targets and their flanking regions. Oxford Nanopore Technologies recently released an USB3.0-interfaced sequencer, the MinION. To date no particular method for enriching MinION libraries has been standardized. Here, using biotinylated PCR-generated baits in a novel approach, we describe a simple and efficient way for multiplexed enrichment of MinION libraries, overcoming technical limitations related with the chemistry of the sequencing-adapters and the length of the DNA fragments. Using Phage Lambda and Escherichia coli as models we selectively enrich for specific targets, significantly increasing the corresponding read-coverage, eliminating unwanted regions. We show that by capturing genomic fragments, which contain the target sequences, we recover reads extending targeted regions and thus can be used for the determination of potentially unknown flanking sequences. By pooling enriched libraries derived from two distinct E. coli strains and analyzing them in parallel, we demonstrate the efficiency of this method in multiplexed format. Crucially we evaluated the optimal bait size for large fragment libraries and we describe for the first time a standardized method for target enrichment in MinION platform.     

Natural genetic variation in MIR172 isolated from Brassica species

The present study reports a natural variation in microRNA172 (MIR172) family members isolated from six species of genus Brassica. The analysis of nucleotide polymorphism across 44 Brassica MIR172 homologs revealed a higher conservation in the predicted precursors relative to flanking regions. Single nucleotide polymorphisms (SNPs) were detected in miRNA and miRNA*. The 21-nt miRNA sequence was conserved in all MIR172 members except MIR172a. However, the miRNA* sequence was conserved only in MIR172a compared to A. thaliana. Non-canonical Brassica variants of precursor miR172a were detected wherein SNP at 5′ terminal in mature miR172a resulted in a sequence identical to mature miR172e. SNPs and indels in precursors resulted in varied stem-loop structures of differing stabilities (ΔG) implying a differential efficiency of miRNA biogenesis. A sequence based phylogram revealed ortholog specific groupings of MIR172 irrespective of genetic background. A Northern analysis in Brassica juncea displayed the cumulative expression of miR172 isoforms in all tissues representing different developmental stages with levels gradually increasing from vegetative to reproductive stages. Detection of high content of miR172 in roots indicates the possibility of additional roles of Brassica miR172 in root development.     

Ligation-mediated suppression-PCR as a powerful tool to analyse nuclear gene sequences in the green alga Chlamydomonas reinhardtii

To improve the analysis of unknown flanking DNA sequences adjacent to known sequences in nuclear genomes of photoautotrophic eukaryotic organisms, we established the technique of ligation-mediated suppression-PCR (LMS-PCR) in the green alga Chlamydomonas reinhardtii for (1) walking from a specific nuclear insertion fragment of random knockout mutants into the unknown flanking DNA sequence to identify and analyse disrupted genomic DNA regions and for (2) walking from highly conserved DNA regions derived from known gene iso-forms into flanking DNA sequences to identify new members of protein families. The feasibility of LMS-PCR for these applications was successfully demonstrated in two different approaches. The first resulted in the identification of a genomic DNA fragment flanking a nuclear insertion vector in a random knockout mutant whose phenotype was characterised by its inability to perform functional LHC state transitions. The second approach targeted the cab gene family. An oligonucleotide of a cabII gene, derived from a highly conserved region, was used to identify potential cab gene regions in the nuclear genome of Chlamydomonas. LMS-PCR combined with 3′ rapid amplification of cDNA ends (3′ RACE) and a PCR-based screening of a cDNA library resulted in the identification of the new cabII gene lhcb4. Both results clearly indicate that LMS-PCR is a powerful tool for the identification of flanking DNA sequences in the nuclear genome of Chlamydomonas reinhardtii. This revised version was published online in June 2006 with corrections to the Cover Date.     

RNA polymerase V-dependent small RNAs in Arabidopsis originate from small,intergenic loci including most SINE repeats

In plants, heterochromatin is maintained by a small RNA-based gene silencing mechanism known as RNA-directed DNA methylation (RdDM). RdDM requires the non-redundant functions of two plant-specific DNA-dependent RNA polymerases (RNAP), RNAP IV and RNAP V. RNAP IV plays a major role in siRNA biogenesis, while RNAP V may recruit DNA methylation machinery to target endogenous loci for silencing. Although small RNA-generating regions that are dependent on both RNAP IV and RNAP V have been identified previously, the genomic loci targeted by RNAP V for siRNA accumulation and silencing have not been described extensively. To characterize the RNAP V-dependent, heterochromatic siRNA-generating regions in the Arabidopsis genome, we deeply sequenced the small RNA populations of wild-type and RNAP V null mutant (nrpe1) plants. Our results showed that RNAP V-dependent siRNA-generating loci are associated predominately with short repetitive sequences in intergenic regions. Suppression of small RNA production from short repetitive sequences was also prominent in RdDM mutants including dms4, drd1, dms3 and rdm1, reflecting the known association of these RdDM effectors with RNAP V. The genomic regions targeted by RNAP V were small, with an estimated average length of 238 bp. Our results suggest that RNAP V affects siRNA production from genomic loci with features dissimilar to known RNAP IV-dependent loci. RNAP V, along with RNAP IV and DRM1/2, may target and silence a set of small, intergenic transposable elements located in dispersed genomic regions for silencing. Silencing at these loci may be actively reinforced by RdDM.     

Linkage of two distinct AT-rich minisatellites at multiple loci in the genome of Theileria parva

Minisatellite tandem repeat elements are well known components of vertebrate genomes, but have not yet been extensively characterized in lower eukaryotes. We describe two unusual, AT-rich minisatellites of the protozoan parasite Theileria parva whose sequences are unrelated to the G/C-rich `chi minisatellite superfamily' of vertebrate and plant genomes. The T. parva tandem repeats, one with a conserved sequence T_2-5ACACA (6–17 copies), and the other with a 6-bp core sequence of either ACTATA or TATACT associated with additional variable sequences in repeats of 10–17 bp (3–7 copies), were closely linked at more than 20 sites in the T. parva genome, separated by 390, 510 and 660 bp at three loci analysed in detail. Such linkage is without precedent in minisatellites so far analysed in other organisms. The minisatellite loci were widely dispersed on 13 out of 33 genomic SfiI fragments, on all four T. parva chromosomes and did not exhibit a telomeric bias in their distribution. Analysis of flanking sequences revealed no obvious conserved sequences between the five loci, or other multicopy repeat sequences outside the minisatellite regions. The T_2-5 ACACA minisatellite was highly effective as a multilocus fingerprinting probe for discrimination of T. parva isolates. Analysis of two individual minisatellite loci revealed variation between the genomic DNAs of two T. parva isolates in the copy number of the constituent repeats within the array, similar to that typical of vertebrate minisatellites.     

The mechanism of gene targeting in Physcomitrella patens: homologous recombination, concatenation and multiple integration

The model bryophyte Physcomitrella patens exhibits high frequencies of gene targeting when transformed with DNA constructs containing sequences homologous with genomic loci. ‘Targeted gene replacement’ (TGR) resulting from homologous recombination (HR) between each end of a targeting construct and the targeted locus occurs when either single or multiple targeting vectors are delivered. In the latter instance simultaneous, multiple, independent integration of different transgenes occurs at the targeted loci. In both single gene and ‘batch’ transformations, DNA can also be found to undergo ‘targeted insertion’ (TI), integrating at one end of the targeted locus by HR with one flanking sequence of the vector accompanied by an apparent non-homologous end-joining (NHEJ) event at the other. Untargeted integration at nonhomologous sites also occurs, but at a lower frequency. Molecular analysis of TI at a single locus shows that this occurs as a consequence of concatenation of the transforming DNA, in planta, prior to integration, followed by HR between a single site in the genomic target and two of its repeated homologues in the concatenated vector. This reinforces the view that HR is the major pathway by which transforming DNA is integrated in Physcomitrella.     

Nucleotide sequence of a γ gliadin gene: Comparisons with other γ gliadin sequences show the structure of γ gliadin genes and the general primary structure of γ gliadins

A low stringency screening of a wheat (Triticum aestivum L.) genomic library produced three types of γ gliadin clones. The sequence of one clone, λ10–20, encoded a γ gliadin of 34.3 kDa. Comparisons of this protein with the proteins encoded by other γ gliadin DNA sequences revealed a general γ gliadin structure: a 19-residue signal peptide; 12-residue mature amino terminus; 12–16 copies of a proline- and glutamine-rich heptapeptide repeat; a 76-residue region high in glutamine containing most of the cyysteines and charged residues: a 6–16-residue polyglutamine region; and the 41-residue carboxyl terminus. Comparisons of the 5′ and 3′ flanking regions of several γ gliadins reveals the high homology and general structure of γ gliadin genes. A further comparison of the 5′ flanking regions with the 5′ flanking regions of other prolamin genes showed that γ gliadin genes contain three copies of a conserved sequence seen within approx. 600 b.p. upstream of the translation start sites of prolamin genes.     

A Phylogenomic Perspective on the Radiation of Ray-Finned Fishes Based upon Targeted Sequencing of Ultraconserved Elements (UCEs)

Ray-finned fishes constitute the dominant radiation of vertebrates with over 32,000 species. Although molecular phylogenetics has begun to disentangle major evolutionary relationships within this vast section of the Tree of Life, there is no widely available approach for efficiently collecting phylogenomic data within fishes, leaving much of the enormous potential of massively parallel sequencing technologies for resolving major radiations in ray-finned fishes unrealized. Here, we provide a genomic perspective on longstanding questions regarding the diversification of major groups of ray-finned fishes through targeted enrichment of ultraconserved nuclear DNA elements (UCEs) and their flanking sequence. Our workflow efficiently and economically generates data sets that are orders of magnitude larger than those produced by traditional approaches and is well-suited to working with museum specimens. Analysis of the UCE data set recovers a well-supported phylogeny at both shallow and deep time-scales that supports a monophyletic relationship between Amia and Lepisosteus (Holostei) and reveals elopomorphs and then osteoglossomorphs to be the earliest diverging teleost lineages. Our approach additionally reveals that sequence capture of UCE regions and their flanking sequence offers enormous potential for resolving phylogenetic relationships within ray-finned fishes.     

The impact of age,biogenesis, and genomic clustering on Drosophila microRNA evolution

The molecular evolutionary signatures of miRNAs inform our understanding of their emergence, biogenesis, and function. The known signatures of miRNA evolution have derived mostly from the analysis of deeply conserved, canonical loci. In this study, we examine the impact of age, biogenesis pathway, and genomic arrangement on the evolutionary properties of Drosophila miRNAs. Crucial to the accuracy of our results was our curation of high-quality miRNA alignments, which included nearly 150 corrections to ortholog calls and nucleotide sequences of the global 12-way Drosophilid alignments currently available. Using these data, we studied primary sequence conservation, normalized free-energy values, and types of structure-preserving substitutions. We expand upon common miRNA evolutionary patterns that reflect fundamental features of miRNAs that are under functional selection. We observe that melanogaster-subgroup-specific miRNAs, although recently emerged and rapidly evolving, nonetheless exhibit evolutionary signatures that are similar to well-conserved miRNAs and distinct from other structured noncoding RNAs and bulk conserved non-miRNA hairpins. This provides evidence that even young miRNAs may be selected for regulatory activities. More strikingly, we observe that mirtrons and clustered miRNAs both exhibit distinct evolutionary properties relative to solo, well-conserved miRNAs, even after controlling for sequence depth. These studies highlight the previously unappreciated impact of biogenesis strategy and genomic location on the evolutionary dynamics of miRNAs, and affirm that miRNAs do not evolve as a unitary class.     

Development of a versatile cassette for directional genome walking using cassette ligation-mediated PCR and its application in the cloning of complete lipolytic genes from Bacillus species

Since the invention of the PCR technology, adaptation techniques to clone DNA fragments flanking the known sequence continue to be developed. We describe a perfectly annealed cassette available in almost unlimited quantities with variable sticky-and blunt-end restriction enzyme recognition sites for efficient restriction and ligation with the restricted target genomic DNA. The cassette provides a 200-bp sequence, which is used to design a variety of cassette-specific primers. The dephosphorylation prevents cassette self-ligation and creates a nick at the cassette: target genome DNA ligation site suppressing unspecific PCR amplifications. We introduce the single-strand amplification PCR (SSA-PCR) technique where a lone known locus-specific primer is firstly used to enrich the targeted template DNA strand resulting in significant PCR product specificity during the second round conventional nested PCR. The distance between the known locus-specific primer and the nearest location of the restriction enzyme used determined the length of the obtained PCR product. We used this technique to walk downstream into the isochorismatase and upstream into the hypothetical conserved genes flanking the mature extracellular lipase gene from Bacillus licheniformis. We further demonstrated the potential of the technique as a cost-effective method during PCR-based prospecting for novel genes by designing "universal" degenerate primers that detected homologues of Family VII bacterial lipolytic genes in Bacillus species. The cassette ligation-mediated PCR was used to clone complete nucleotide sequences encoding functional lipolytic genes from B. licheniformis and Bacillus pumilus.     

Microfluidic droplet enrichment for targeted sequencing

Targeted sequence enrichment enables better identification of genetic variation by providing increased sequencing coverage for genomic regions of interest. Here, we report the development of a new target enrichment technology that is highly differentiated from other approaches currently in use. Our method, MESA (Microfluidic droplet Enrichment for Sequence Analysis), isolates genomic DNA fragments in microfluidic droplets and performs TaqMan PCR reactions to identify droplets containing a desired target sequence. The TaqMan positive droplets are subsequently recovered via dielectrophoretic sorting, and the TaqMan amplicons are removed enzymatically prior to sequencing. We demonstrated the utility of this approach by generating an average 31.6-fold sequence enrichment across 250 kb of targeted genomic DNA from five unique genomic loci. Significantly, this enrichment enabled a more comprehensive identification of genetic polymorphisms within the targeted loci. MESA requires low amounts of input DNA, minimal prior locus sequence information and enriches the target region without PCR bias or artifacts. These features make it well suited for the study of genetic variation in a number of research and diagnostic applications.