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 RNAs are ubiquitous, versatile repressors and include (1) microRNAs (miRNAs), processed from mRNA forming stem-loops; and (2) small interfering RNAs (siRNAs), the latter derived in plants by a process typically requiring an RNAdependent RNA polymerase. We constructed and analyzed an expression atlas of soybean (Glycine max) small RNAs, identifying over 500 loci generating 21-nucleotide phased siRNAs (phasiRNAs; from PHAS loci), of which 483 overlapped annotated protein-coding genes. Via the integration of miRNAs with parallel analysis of RNA end (PARE) data, 20 miRNA triggers of 127 PHAS loci were detected. The primary class of PHAS loci (208 or 41% of the total) corresponded to NB-LRR genes; some of these small RNAs preferentially accumulate in nodules. Among the PHAS loci, novel representatives of TAS3 and noncanonical phasing patterns were also observed. A noncoding PHAS locus, triggered by miR4392, accumulated preferentially in anthers; the phasiRNAs are predicted to target transposable elements, with their peak abundance during soybean reproductive development. Thus, phasiRNAs show tremendous diversity in dicots. We identified novel miRNAs and assessed the veracity of soybean miRNAs registered in miRBase, substantially improving the soybean miRNA annotation, facilitating an improvement of miRBase annotations and identifying at high stringency novel miRNAs and their targets.

Although successful in identifying new cataract-linked genes, the previous version of the database iSyTE (integrated Systems Tool for Eye gene discovery) was based on expression information on just three mouse lens stages and was functionally limited to visualization by only UCSC-Genome Browser tracks. To increase its efficacy, here we provide an enhanced iSyTE version 2.0 (URL: http://research.bioinformatics.udel.edu/iSyTE) based on well-curated, comprehensive genome-level lens expression data as a one-stop portal for the effective visualization and analysis of candidate genes in lens development and disease. iSyTE 2.0 includes all publicly available lens Affymetrix and Illumina microarray datasets representing a broad range of embryonic and postnatal stages from wild-type and specific gene-perturbation mouse mutants with eye defects. Further, we developed a new user-friendly web interface for direct access and cogent visualization of the curated expression data, which supports convenient searches and a range of downstream analyses. The utility of these new iSyTE 2.0 features is illustrated through examples of established genes associated with lens development and pathobiology, which serve as tutorials for its application by the end-user. iSyTE 2.0 will facilitate the prioritization of eye development and disease-linked candidate genes in studies involving transcriptomics or next-generation sequencing data, linkage analysis and GWAS approaches.

Opacification of the ocular lens, termed cataract, is a common cause of blindness. To become transparent, lens fiber cells undergo degradation of their organelles, including their nuclei, presenting a fundamental question: does signaling/transcription sufficiently explain differentiation of cells progressing toward compromised transcriptional potential? We report that a conserved RNA-binding protein Celf1 post-transcriptionally controls key genes to regulate lens fiber cell differentiation. Celf1-targeted knockout mice and celf1-knockdown zebrafish and Xenopus morphants have severe eye defects/cataract. Celf1 spatiotemporally down-regulates the cyclin-dependent kinase (Cdk) inhibitor p27Kip1 by interacting with its 5’ UTR and mediating translation inhibition. Celf1 deficiency causes ectopic up-regulation of p21Cip1. Further, Celf1 directly binds to the mRNA of the nuclease Dnase2b to maintain its high levels. Together these events are necessary for Cdk1-mediated lamin A/C phosphorylation to initiate nuclear envelope breakdown and DNA degradation in fiber cells. Moreover, Celf1 controls alternative splicing of the membrane-organization factor beta-spectrin and regulates F-actin-crosslinking factor Actn2 mRNA levels, thereby controlling fiber cell morphology. Thus, we illustrate new Celf1-regulated molecular mechanisms in lens development, suggesting that post-transcriptional regulatory RNA-binding proteins have evolved conserved functions to control vertebrate oculogenesis.

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.

Although majority of the genes linked to early-onset cataract exhibit lens fiber cell-enriched expression, our understanding of gene regulation in these cells is limited to function of just eight transcription factors and largely in the context of crystallins. We report on small Maf transcription factors Mafg and Mafk as regulators of several non-crystallin human cataract-associated genes in fiber cells and establish their significance to this disease. We applied a bioinformatics tool for cataract gene discovery iSyTE to identify Mafg and its co-regulators in the lens, and generated various null-allelic combinations of Mafg:Mafk mouse mutants for phenotypic and molecular analysis. By age 4-months, Mafg−/−:Mafk+/− mutants exhibit lens defects that progressively develop into cataract. High-resolution phenotypic characterization of Mafg−/−:Mafk+/− mouse lens reveals severely disorganized fiber cells, while microarrays-based expression profiling identifies 97 differentially regulated genes (DRGs). Integrative analysis of Mafg−/−:Mafk+/− lens-DRGs with 1) binding-motifs and genomic targets of small Mafs and their regulatory partners, 2) iSyTE lens-expression data, and 3) interactions between DRGs in the String database, unravels a detailed small Maf regulatory network in the lens, several nodes of which are linked to cataract. This approach identifies 36 high-priority candidates from the original 97 DRGs. Significantly, 8/36 (22%) DRGs are associated with cataracts in human (GSTO1, MGST1, SC4MOL, UCHL1) or mouse (Aldh3a1, Crygf, Hspb1, Pcbd1), suggesting a multifactorial etiology that includes oxidative stress and mis-regulation of sterol synthesis. These data identify Mafg and Mafk as new cataract-associated candidates and define their function in regulating largely non-crystallin genes linked to human cataract.

In grasses, two pathways that generate diverse and numerous 21-nt (premeiotic) and 24-nt (meiotic) phased siRNAs are highly enriched in anthers, the male reproductive organs. These "phasiRNAs" are analogous to mammalian piRNAs, yet their functions and evolutionary origins remain largely unknown. The 24-nt meiotic phasiRNAs have only been described in grasses, wherein their biogenesis is dependent on a specialized Dicer (DCL5). To assess how evolution gave rise to this pathway, we examined reproductive phasiRNA pathways in nongrass monocots: garden asparagus, daylily, and lily. The common ancestors of these species diverged approximately 115-117 million years ago (MYA). We found that premeiotic 21-nt and meiotic 24-nt phasiRNAs were abundant in all three species and displayed spatial localization and temporal dynamics similar to grasses. The miR2275-triggered pathway was also present, yielding 24-nt reproductive phasiRNAs, and thus originated more than 117 MYA. In asparagus, unlike in grasses, these siRNAs are largely derived from inverted repeats (IRs); analyses in lily identified thousands of precursor loci, and many were also predicted to form foldback substrates for Dicer processing. Additionally, reproductive phasiRNAs were present in female reproductive organs and thus may function in both male and female germinal development. These data describe several distinct mechanisms of production for 24-nt meiotic phasiRNAs and provide new insights into the evolution of reproductive phasiRNA pathways in monocots.

Isolated or syndromic congenital cataract are heterogeneous developmental defects, making the identification of their associated genes challenging. In the past, mouse lens expression microarrays have been successfully applied in bioinformatics tools (e.g. iSyTE) to facilitate human cataract-associated gene discovery. To develop a new resource for geneticists, we report high-throughput RNA-sequencing (RNA-seq) profiles of mouse lens at key embryonic stages (E)10.5 (lens pit), E12.5 (primary fiber cell differentiation), E14.5 and E16.5 (secondary fiber cell differentiation). These stages capture important events as the lens develops from an invaginating placode into a transparent tissue. Previously, in silico whole-embryo body (WB)-subtraction-based “lens-enriched” expression has been effective in prioritizing cataract-linked genes. To apply an analogous approach, we generated new mouse WB RNA-seq datasets and show that in silico WB-subtraction of lens RNA-seq datasets successfully identifies key genes based on lens-enriched expression. At ≥2 counts-per-million expression, ≥1.5 log2 fold-enrichment (p<0.05) cut-off, E10.5 lens exhibits 1401 enriched-genes (17% lens-expressed genes), E12.5 lens exhibits 1937 enriched-genes (22% lens-expressed genes), E14.5 lens exhibits 2514 enriched-genes (31% lens-expressed genes), and E16.5 lens exhibits 2745 enriched-genes (34% lens-expressed genes). Biological pathway analysis identified genes associated with lens development, transcription regulation and signaling pathways, among other functional groups. Furthermore, these new RNA-seq data confirmed high expression of established cataract-linked genes and identified new potential regulators in the lens. Finally, we developed new lens stage-specific UCSC Genome Brower annotation-tracks and made these publicly accessible through iSyTE (https://research.bioinformatics.udel.edu/iSyTE/) for user-friendly visualization of lens gene expression/enrichment to prioritize genes from high-throughput data from cataract cases.