The evolutionary history of small RNAs in Solanaceae.

The Solanaceae or "nightshade" family is an economically important group with remarkable diversity. To gain a better understanding of how the unique biology of the Solanaceae relates to the family's small RNA (sRNA) genomic landscape, we downloaded over 255 publicly available sRNA data sets that comprise over 2.6 billion reads of sequence data. We applied a suite of computational tools to predict and annotate two major sRNA classes: (1) microRNAs (miRNAs), typically 20- to 22-nucleotide (nt) RNAs generated from a hairpin precursor and functioning in gene silencing and (2) short interfering RNAs (siRNAs), including 24-nt heterochromatic siRNAs typically functioning to repress repetitive regions of the genome via RNA-directed DNA methylation, as well as secondary phased siRNAs and trans-acting siRNAs generated via miRNA-directed cleavage of a polymerase II-derived RNA precursor. Our analyses described thousands of sRNA loci, including poorly understood clusters of 22-nt siRNAs that accumulate during viral infection. The birth, death, expansion, and contraction of these sRNA loci are dynamic evolutionary processes that characterize the Solanaceae family. These analyses indicate that individuals within the same genus share similar sRNA landscapes, whereas comparisons between distinct genera within the Solanaceae reveal relatively few commonalities.

Soybean DICER-LIKE2 Regulates Seed Coat Color via Production of Primary 22-Nucleotide Small Interfering RNAs from Long Inverted Repeats.

In plants, 22-nucleotide small RNAs trigger the production of secondary small interfering RNAs (siRNAs) and enhance silencing. DICER-LIKE2 (DCL2)-dependent 22-nucleotide siRNAs are rare in Arabidopsis (Arabidopsis thaliana) and are thought to function mainly during viral infection; by contrast, these siRNAs are abundant in many crops such as soybean (Glycine max) and maize (Zea mays). Here, we studied soybean 22-nucleotide siRNAs by applying CRISPR-Cas9 to simultaneously knock out the two copies of soybean DCL2, GmDCL2a and GmDCL2b, in the Tianlong1 cultivar. Small RNA sequencing revealed that most 22-nucleotide siRNAs are derived from long inverted repeats (LIRs) and disappeared in the Gmdcl2a/2b double mutant. De novo assembly of a Tianlong1 reference genome and transcriptome profiling identified an intronic LIR formed by the chalcone synthase (CHS) genes CHS1 and CHS3 This LIR is the source of primary 22-nucleotide siRNAs that target other CHS genes and trigger the production of secondary 21-nucleotide siRNAs. Disruption of this process in Gmdcl2a/2b mutants substantially increased CHS mRNA levels in the seed coat, thus changing the coat color from yellow to brown. Our results demonstrated that endogenous LIR-derived transcripts in soybean are predominantly processed by GmDCL2 into 22-nucleotide siRNAs and uncovered a role for DCL2 in regulating natural traits.

Next-Generation Sequence Databases: RNA and Genomic Informatics Resources for Plants.

We developed public web sites and resources for data access, display, and analysis of plant small RNAs. These web sites are interconnected with related data types. The current generation of these informatics tools was developed for Illumina data, evolving over more than 15 years of improvements. Our online databases have customized web interfaces to uniquely handle and display RNA-derived data from diverse plant species, ranging from Arabidopsis (Arabidopsis thaliana) to wheat (Triticum spp.), including many crop and model species. The web interface displays the abundance and genomic context of data for small RNAs, parallel analysis of RNA ends/degradome reads, RNA sequencing, and even chromatin immunoprecipitation sequencing data; it also provides information about potentially novel transcripts (antisense transcripts, alternative splice isoforms, and regulatory intergenic transcripts). Numerous options are included for downloading data as tables or via web services. Interpretation of these data is facilitated by the inclusion of extensive repeat or transposon data in our genome viewer. We have developed graphical and analytical tools, including a new viewer and a query page for the analysis of phased small RNAs; these are particularly useful for understanding the complex small RNA pathways of plants. These public databases are accessible at https://mpss.danforthcenter.org.

The asparagus genome sheds light on the origin and evolution of a young Y chromosome.

Sex chromosomes evolved from autosomes many times across the eukaryote phylogeny. Several models have been proposed to explain this transition, some involving male and female sterility mutations linked in a region of suppressed recombination between X and Y (or Z/W, U/V) chromosomes. Comparative and experimental analysis of a reference genome assembly for a double haploid YY male garden asparagus (Asparagus officinalis L.) individual implicates separate but linked genes as responsible for sex determination. Dioecy has evolved recently within Asparagus and sex chromosomes are cytogenetically identical with the Y, harboring a megabase segment that is missing from the X. We show that deletion of this entire region results in a male-to-female conversion, whereas loss of a single suppressor of female development drives male-to-hermaphrodite conversion. A single copy anther-specific gene with a male sterile Arabidopsis knockout phenotype is also in the Y-specific region, supporting a two-gene model for sex chromosome evolution.

Small RNA Functions Are Required for Growth and Development of Magnaporthe oryzae.

RNA interference (RNAi) is conserved in eukaryotic organisms, and it has been well studied in many animal and plant species and some fungal species, yet it is not well studied in fungal plant pathogens. In the rice blast fungus Magnaporthe oryzae, we examined small RNA (sRNA) and their biogenesis in the context of growth and pathogenicity. Through genetic and genomic analyses, we demonstrate that loss of a single gene encoding Dicer, RNA-dependent RNA polymerase, or Argonaute reduces sRNA levels. These three proteins are required for the biogenesis of sRNA-matching genome-wide regions (coding regions, repeats, and intergenic regions). The loss of one Argonaute reduced both sRNA and fungal virulence on barley leaves. Transcriptome analysis of multiple mutants revealed that sRNA play an important role in transcriptional regulation of repeats and intergenic regions in M. oryzae. Together, these data support that M. oryzae sRNA regulate developmental processes including, fungal growth and virulence.

miTRATA: a web-based tool for microRNA Truncation and Tailing Analysis.

SUMMARY: We describe miTRATA, the first web-based tool for microRNA Truncation and Tailing Analysis--the analysis of 3' modifications of microRNAs including the loss or gain of nucleotides relative to the canonical sequence. miTRATA is implemented in Python (version 3) and employs parallel processing modules to enhance its scalability when analyzing multiple small RNA (sRNA) sequencing datasets. It utilizes miRBase, currently version 21, as a source of known microRNAs for analysis. miTRATA notifies user(s) via email to download as well as visualize the results online. miTRATA's strengths lie in (i) its biologist-focused web interface, (ii) improved scalability via parallel processing and (iii) its uniqueness as a webtool to perform microRNA truncation and tailing analysis. AVAILABILITY AND IMPLEMENTATION: miTRATA is developed in Python and PHP. It is available as a web-based application from https://wasabi.dbi.udel.edu/∼apps/ta/. CONTACT: meyers@dbi.udel.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

sPARTA: a parallelized pipeline for integrated analysis of plant miRNA and cleaved mRNA data sets, including new miRNA target-identification software.

Parallel analysis of RNA ends (PARE) is a technique utilizing high-throughput sequencing to profile uncapped, mRNA cleavage or decay products on a genome-wide basis. Tools currently available to validate miRNA targets using PARE data employ only annotated genes, whereas important targets may be found in unannotated genomic regions. To handle such cases and to scale to the growing availability of PARE data and genomes, we developed a new tool, 'sPARTA' (small RNA-PARE target analyzer) that utilizes a built-in, plant-focused target prediction module (aka 'miRferno'). sPARTA not only exhibits an unprecedented gain in speed but also it shows greater predictive power by validating more targets, compared to a popular alternative. In addition, the novel 'seed-free' mode, optimized to find targets irrespective of complementarity in the seed-region, identifies novel intergenic targets. To fully capitalize on the novelty and strengths of sPARTA, we developed a web resource, 'comPARE', for plant miRNA target analysis; this facilitates the systematic identification and analysis of miRNA-target interactions across multiple species, integrated with visualization tools. This collation of high-throughput small RNA and PARE datasets from different genomes further facilitates re-evaluation of existing miRNA annotations, resulting in a 'cleaner' set of microRNAs.

An integrated workflow for DNA methylation analysis.

The analysis of cytosine methylation provides a new way to assess and describe epigenetic regulation at a whole-genome level in many eukaryotes. DNA methylation has a demonstrated role in the genome stability and protection, regulation of gene expression and many other aspects of genome function and maintenance. BS-seq is a relatively unbiased method for profiling the DNA methylation, with a resolution capable of measuring methylation at individual cytosines. Here we describe, as an example, a workflow to handle DNA methylation analysis, from BS-seq library preparation to the data visualization. We describe some applications for the analysis and interpretation of these data. Our laboratory provides public access to plant DNA methylation data via visualization tools available at our "Next-Gen Sequence" websites (http://mpss.udel.edu), along with small RNA, RNA-seq and other data types.

Roles of DCL4 and DCL3b in rice phased small RNA biogenesis.

Higher plants have evolved multiple proteins in the RNase III family to produce and regulate different classes of small RNAs with specialized molecular functions. In rice (Oryza sativa), numerous genomic clusters are targeted by one of two microRNAs (miRNAs), miR2118 and miR2275, to produce secondary small interfering RNAs (siRNAs) of either 21 or 24 nucleotides in a phased manner. The biogenesis requirements or the functions of the phased small RNAs are completely unknown. Here we examine the rice Dicer-Like (DCL) family, including OsDCL1, -3a, -3b and -4. By deep sequencing of small RNAs from different tissues of the wild type and osdcl4-1, we revealed that the processing of 21-nucleotide siRNAs, including trans-acting siRNAs (tasiRNA) and over 1000 phased small RNA loci, was largely dependent on OsDCL4. Surprisingly, the processing of 24-nucleotide phased small RNA requires the DCL3 homolog OsDCL3b rather than OsDCL3a, suggesting functional divergence within DCL3 family. RNA ligase-mediated 5' rapid amplification of cDNA ends and parallel analysis of RNA ends (PARE)/degradome analysis confirmed that most of the 21- and 24-nucleotide phased small RNA clusters were initiated from the target sites of miR2118 and miR2275, respectively. Furthermore, the accumulation of the two triggering miRNAs requires OsDCL1 activity. Finally, we show that phased small RNAs are preferentially produced in the male reproductive organs and are likely to be conserved in monocots. Our results revealed significant roles of OsDCL4, OsDCL3b and OsDCL1 in the 21- and 24-nucleotide phased small RNA biogenesis pathway in rice.

A comparative study on testicular microstructure and relative sperm production in gorillas, chimpanzees, and orangutans.

We performed histological analyses for comparing testicular microstructure between the gorilla, chimpanzee, and orangutan. Testicular samples were obtained by autopsy or biopsy from 10 gorillas, 11 chimpanzees, and 7 orangutans from several zoos and institutes. The seminiferous epithelia were thick in the chimpanzee and orangutan but thin in the gorilla. Leydig cells in the interstitial tissue were abundant in the gorilla. The acrosomic system was extremely well developed in the orangutans. Our study reveals that the cycle of seminiferous epithelium in orangutan testis can be divided into ten stages, whereas that in human, chimpanzee, and gorilla testes can be divided into only six stages. Phylogenetic analyses of the number of divisions may indicate that the seminiferous epithelium of our common ancestor has changed since the orangutan diverged from it. Furthermore, we performed comparative analyses of testicular microstructure to estimate relative sperm production among these three animals, and proposed a new indicator (namely the spermatogenic index, SI) closely related to sperm production. The SI indicated that a chimpanzee usually produces about 223 times more sperm than a gorilla and about 14 times more than an orangutan. Our data demonstrate the significance of the SI for estimating sperm production, thus aiding our understanding of the reproductive strategy as well as testis weight and relative testis size in investigated primates.

The influence of excessive IL-6 production in vivo on the development and function of Foxp3+ regulatory T cells.

IL-6 is a proinflammatory cytokine and its overproduction is implicated in a variety of inflammatory disorders. Recent in vitro analyses suggest that IL-6 is a key cytokine that determines the balance between Foxp3(+) regulatory T cells (Tregs) and Th17 cells. However, it remains unclear whether excessive IL-6 production in vivo alters the development and function of Foxp3(+) Tregs. In this study, we analyzed IL-6 transgenic (Tg) mice in which serum IL-6 levels are constitutively elevated. Interestingly, in IL-6 Tg mice, whereas peripheral lymphoid organs were enlarged, and T cells exhibited activated phenotype, Tregs were not reduced but rather increased compared with wild-type mice. In addition, Tregs from Tg mice normally suppressed proliferation of naive T cells in vitro. Furthermore, Tregs cotransferred with naive CD4 T cells into SCID-IL-6 Tg mice inhibited colitis as successfully as those transferred into control SCID mice. These results indicate that overproduction of IL-6 does not inhibit development or function of Foxp3(+) Tregs in vivo. However, when naive CD4 T cells alone were transferred, Foxp3(+) Tregs retrieved from SCID-IL-6 Tg mice were reduced compared with SCID mice. Moreover, the Helios(-) subpopulation of Foxp3(+) Tregs, recently defined as extrathymic Tregs, was significantly reduced in IL-6 Tg mice compared with wild-type mice. Collectively, these results suggest that IL-6 overproduced in vivo inhibits inducible Treg generation from naive T cells, but does not affect the development and function of natural Tregs.

Bioinformatics analysis of small RNAs in plants using next generation sequencing technologies.

Next-generation sequencing technologies have a substantial impact on a broad range of biological applications. Like many other groups, we use these new technologies, especially SBS (Sequence-By-Synthesis), for deep profiling of small RNA molecules in plants. Small RNAs are 21-24 nucleotides in length and are known to play a major role in the activation of mRNAs and genomic DNAs. We have generated numerous SBS small RNA libraries; each can consist of more than three million signatures of more than 33 nucleotides in length. Here, we describe the challenges and our strategies to handle the very large quantity of small RNA data generated by these next-generation sequencing technologies.

Topography-guided (NIDEK customized aspheric treatment zone) photorefractive keratectomy with mitomycin C after penetrating keratoplasty for keratoconus: case report.

PURPOSE: To report topography-guided photorefractive keratectomy (PRK) with mitomycin C (MMC) after penetrating keratoplasty for keratoconus. METHODS: A 34-year-old woman with irregular astigmatism after penetrating keratoplasty in the right eye underwent PRK. Topography-guided surface ablation using the customized aspheric treatment zone ablation (CATz) was programmed for a 5.00-mm optical zone and an 8.50-mm transition zone. Mitomycin C was applied to corneal stroma for 20 seconds immediately after excimer laser ablation. RESULTS: Two months after surgery, the patient was very satisfied with the refractive outcome. Best spectacle-corrected visual acuity was 20/20 with a manifest refraction of -0.50 diopters. No haze formation was detected at 6-month follow-up. A significant improvement in the central 3 mm of the axial topography was noted. CONCLUSIONS: Topography-guided ablation using MMC may be a reasonable alternative for the management of refractive error after penetrating keratoplasty.

Plant MPSS databases: signature-based transcriptional resources for analyses of mRNA and small RNA.

MPSS (massively parallel signature sequencing) is a sequencing-based technology that uses a unique method to quantify gene expression level, generating millions of short sequence tags per library. We have created a series of databases for four species (Arabidopsis, rice, grape and Magnaporthe grisea, the rice blast fungus). Our MPSS databases measure the expression level of most genes under defined conditions and provide information about potentially novel transcripts (antisense transcripts, alternative splice isoforms and regulatory intergenic transcripts). A modified version of MPSS has been used to perform deep profiling of small RNAs from Arabidopsis, and we have recently adapted our database to display these data. Interpretation of the small RNA MPSS data is facilitated by the inclusion of extensive repeat data in our genome viewer. All the data and the tools introduced in this article are available at http://mpss.udel.edu.

Testicular histological examination of spermatogenetic activity in captive gorillas (Gorilla gorilla).

To clarify the reproductive state of male gorillas, we performed histological examinations on the testicles of 10 male gorillas (Gorilla gorilla). The testicular samples were obtained by autopsy, and ordinal histological preparations were made for light microscopy. The poor spermatogenesis of this species was characterized by the following findings: First, spermatogenesis was evident in only four samples. Meiosis progressed in two samples, but they lacked spermatogenesis. In the remaining four specimens, seminiferous tubules hyalinized without any sign of spermatogenesis. Second, seminiferous epithelia were thin even in the males in which spermatogenesis was observed. Third, degenerated seminiferous tubules were found in all specimens. Fourth, abnormally large syncytial cells were found in the tubules. Six stages in the epithelial cycle of the seminiferous tubules were identified. Testosterone staining made it clear that there were many Leydig cells with spherical or fusiform nuclei in an abundance of interstitial tissue. The relevance of the testicular architecture of gorillas to the mating system is discussed.