Conserved piRNA Expression from a Distinct Set of piRNA Cluster Loci in Eutherian Mammals.

The Piwi pathway is deeply conserved amongst animals because one of its essential functions is to repress transposons. However, many Piwi-interacting RNAs (piRNAs) do not base-pair to transposons and remain mysterious in their targeting function. The sheer number of piRNA cluster (piC) loci in animal genomes and infrequent piRNA sequence conservation also present challenges in determining which piC loci are most important for development. To address this question, we determined the piRNA expression patterns of piC loci across a wide phylogenetic spectrum of animals, and reveal that most genic and intergenic piC loci evolve rapidly in their capacity to generate piRNAs, regardless of known transposon silencing function. Surprisingly, we also uncovered a distinct set of piC loci with piRNA expression conserved deeply in Eutherian mammals. We name these loci Eutherian-Conserved piRNA cluster (ECpiC) loci. Supporting the hypothesis that conservation of piRNA expression across ~100 million years of Eutherian evolution implies function, we determined that one ECpiC locus generates abundant piRNAs antisense to the STOX1 transcript, a gene clinically associated with preeclampsia. Furthermore, we confirmed reduced piRNAs in existing mouse mutations at ECpiC-Asb1 and -Cbl, which also display spermatogenic defects. The Asb1 mutant testes with strongly reduced Asb1 piRNAs also exhibit up-regulated gene expression profiles. These data indicate ECpiC loci may be specially adapted to support Eutherian reproduction.

Solution structure of the MID1 B-box2 CHC(D/C)C(2)H(2) zinc-binding domain: insights into an evolutionarily conserved RING fold.

The B-box type 2 domain is a prominent feature of a large and growing family of RING, B-box, coiled-coil (RBCC) domain-containing proteins and is also present in more than 1500 additional proteins. Most proteins usually contain a single B-box2 domain, although some proteins contain tandem domains consisting of both type 1 and type 2 B-boxes, which actually share little sequence similarity. Recently, we determined the solution structure of B-box1 from MID1, a putative E3 ubiquitin ligase that is mutated in X-linked Opitz G/BBB syndrome, and showed that it adopted a betabetaalpha RING-like fold. Here, we report the tertiary structure of the B-box2 (CHC(D/C)C(2)H(2)) domain from MID1 using multidimensional NMR spectroscopy. This MID1 B-box2 domain consists of a short alpha-helix and a structured loop with two short anti-parallel beta-strands and adopts a tertiary structure similar to the B-box1 and RING structures, even though there is minimal primary sequence similarity between these domains. By mutagenesis, ESI-FTICR and ICP mass spectrometry, we show that the B-box2 domain coordinates two zinc atoms with a 'cross-brace' pattern: one by Cys175, His178, Cys195 and Cys198 and the other by Cys187, Asp190, His204, and His207. Interestingly, this is the first case that an aspartic acid is involved in zinc atom coordination in a zinc-finger domain, although aspartic acid has been shown to coordinate non-catalytic zinc in matrix metalloproteinases. In addition, the finding of a Cys195Phe substitution identified in a patient with X-linked Opitz GBBB syndrome supports the importance of proper zinc coordination for the function of the MID1 B-box2 domain. Notably, however, our structure differs from the only other published B-box2 structure, that from XNF7, which was shown to coordinate one zinc atom. Finally, the similarity in tertiary structures of the B-box2, B-box1 and RING domains suggests these domains have evolved from a common ancestor.

Interdomain interactions regulate the activation of the heme-regulated eIF 2 alpha kinase.

The heme-regulated inhibitor of protein synthesis (HRI) regulates translation through the phosphorylation of the alpha-subunit of eukaryotic initiation factor-2 (eIF 2). While HRI is best known for its activation in response to heme-deficiency, we recently showed that the binding of NO and CO to the N-terminal heme-binding domain (NT-HBD) of HRI activated and suppressed its activity, respectively. Here, we examined the effect of hemin, NO, and CO on the interaction between the NT-HBD and the catalytic domain of HRI (HRI/Delta HBD). Hemin stabilized the interaction of NT-HBD with HRI/Delta HBD, and NO and CO disrupted and stabilized this interaction, respectively. Mutant HRI (Delta H-HRI), lacking amino acids 116-158 from the NT-HBD, was less sensitive to heme-induced inhibition, and mutant NT-HBD lacking these residues did not bind to HRI/Delta HBD. HRI/Delta HBD and Delta H-HRI also activated more readily than HRI in response to heme-deficiency. Thus, HRI's activity is regulated through the modulation of the interaction between its NT-HBD and catalytic domain.

Small RNA library construction from minute biological samples.

Increasingly, the discovery and characterization of small regulatory RNAs from a variety of organisms have all required deep-sequencing methodologies. However, the crux to successful deep-sequencing analysis depends upon optimal construction of a cDNA library compatible for the high-throughput sequencing platform. Challenges to small RNA library constructions arise when dealing with minute tissue samples because certain structural RNA fragments can dominate and mask the desired characterization of regulatory small RNAs like microRNAs (miRNAs), endogenous small interfering RNAs (endo-siRNAs), and Piwi-interacting RNAs (piRNAs). Here, we describe methods that improve the chances of constructing a successful library from small RNAs isolated from minute tissues such as enriched follicle cells from the Drosophila ovarium. Because the ribosomal RNA (rRNA) fragments are frequently the major contaminants in small RNA preparations from minute amounts of tissue, we demonstrate the utility of antisense oligonucleotide depletion and an acryloylaminophenylboronic acid (APB) polyacrylamide gel system for separating the abundant 2S rRNA in Drosophila from endo-siRNAs and piRNAs. Finally, our methodology generates libraries amenable to multiplex sequencing on the Illumina Hi-Seq platform.

Quantitative genome-wide enhancer activity maps for five Drosophila species show functional enhancer conservation and turnover during cis-regulatory evolution.

Phenotypic differences between closely related species are thought to arise primarily from changes in gene expression due to mutations in cis-regulatory sequences (enhancers). However, it has remained unclear how frequently mutations alter enhancer activity or create functional enhancers de novo. Here we use STARR-seq, a recently developed quantitative enhancer assay, to determine genome-wide enhancer activity profiles for five Drosophila species in the constant trans-regulatory environment of Drosophila melanogaster S2 cells. We find that the functions of a large fraction of D. melanogaster enhancers are conserved for their orthologous sequences owing to selection and stabilizing turnover of transcription factor motifs. Moreover, hundreds of enhancers have been gained since the D. melanogaster-Drosophila yakuba split about 11 million years ago without apparent adaptive selection and can contribute to changes in gene expression in vivo. Our finding that enhancer activity is often deeply conserved and frequently gained provides functional insights into regulatory evolution.