Functional Evolution of a cis-Regulatory Module

Lack of knowledge about how regulatory regions evolve in relation to their structure–function may limit the utility of comparative sequence analysis in deciphering cis-regulatory sequences. To address this we applied reverse genetics to carry out a functional genetic complementation analysis of a eukaryotic cis-regulatory module—the even-skipped stripe 2 enhancer—from four Drosophila species. The evolution of this enhancer is non-clock-like, with important functional differences between closely related species and functional convergence between distantly related species. Functional divergence is attributable to differences in activation levels rather than spatiotemporal control of gene expression. Our findings have implications for understanding enhancer structure–function, mechanisms of speciation and computational identification of regulatory modules.     

cis-Regulatory Complexity within a Large Non-Coding Region in the Drosophila Genome

Analysis of cis-regulatory enhancers has revealed that they consist of clustered blocks of highly conserved sequences. Although most characterized enhancers reside near their target genes, a growing number of studies have shown that enhancers located over 50 kb from their minimal promoter(s) are required for appropriate gene expression and many of these ‘long-range’ enhancers are found in genomic regions that are devoid of identified exons. To gain insight into the complexity of Drosophila cis-regulatory sequences within exon-poor regions, we have undertaken an evolutionary analysis of 39 of these regions located throughout the genome. This survey revealed that within these genomic expanses, clusters of conserved sequence blocks (CSBs) are positioned once every 1.1 kb, on average, and that a typical cluster contains multiple (5 to 30 or more) CSBs that have been maintained for at least 190 My of evolutionary divergence. As an initial step toward assessing the cis-regulatory activity of conserved clusters within gene-free genomic expanses, we have tested the in-vivo enhancer activity of 19 consecutive CSB clusters located in the middle of a 115 kb gene-poor region on the 3^rd chromosome. Our studies revealed that each cluster functions independently as a specific spatial/temporal enhancer. In total, the enhancers possess a diversity of regulatory functions, including dynamically activating expression in defined patterns within subsets of cells in discrete regions of the embryo, larvae and/or adult. We also observed that many of the enhancers are multifunctional–that is, they activate expression during multiple developmental stages. By extending these results to the rest of the Drosophila genome, which contains over 70,000 non-coding CSB clusters, we suggest that most function as enhancers.     

Interdomain connecting loop and J loop structures determine cross-species compatibility of PCNA

Eukaryotic proliferating cell nuclear antigen (PCNA) plays an essential role in orchestrating the assembly of the replisome complex, stimulating processive DNA synthesis, and recruiting other regulatory proteins during the DNA damage response. PCNA and its binding partner network are relatively conserved in eukaryotes, and it exhibits extraordinary structural similarity across species. However, despite this structural similarity, the PCNA of a given species is rarely functional in heterologous systems. In this report, we determined the X-ray crystal structure of Neurospora crassa PCNA (NcPCNA) and compared its structure–function relationship with other available PCNA studies to understand this cross-species incompatibility. We found two regions, the interdomain connecting loop (IDCL) and J loop structures, vary significantly among PCNAs. In particular, the J loop deviates in NcPCNA from that in Saccharomyces cerevisiae PCNA (ScPCNA) by 7 Å. Differences in the IDCL structures result in varied binding affinities of PCNAs for the subunit Pol32 of DNA polymerase delta and for T2-amino alcohol, a small-molecule inhibitor of human PCNA. To validate that these structural differences are accountable for functional incompatibility in S. cerevisiae, we generated NcPCNA mutants mimicking IDCL and J loop structures of ScPCNA. Our genetic analyses suggested that NcPCNA mutants are fully functional in S. cerevisiae. The susceptibility of the strains harboring ScPCNA mimics of NcPCNA to various genotoxic agents was similar to that in yeast cells expressing ScPCNA. Taken together, we conclude that in addition to the overall architecture of PCNA, structures of the IDCL and J loop of PCNA are critical determinants of interspecies functional compatibility.     

Asexuality Associated with Marked Genomic Expansion of Tandemly Repeated rRNA and Histone Genes

How does asexual reproduction influence genome evolution? Although is it clear that genomic structural variation is common and important in natural populations, we know very little about how one of the most fundamental of eukaryotic traits—mode of genomic inheritance—influences genome structure. We address this question with the New Zealand freshwater snail Potamopyrgus antipodarum, which features multiple separately derived obligately asexual lineages that coexist and compete with otherwise similar sexual lineages. We used whole-genome sequencing reads from a diverse set of sexual and asexual individuals to analyze genomic abundance of a critically important gene family, rDNA (the genes encoding rRNAs), that is notable for dynamic and variable copy number. Our genomic survey of rDNA in P. antipodarum revealed two striking results. First, the core histone and 5S rRNA genes occur between tandem copies of the 18S–5.8S–28S gene cluster, a unique architecture for these crucial gene families. Second, asexual P. antipodarum harbor dramatically more rDNA–histone copies than sexuals, which we validated through molecular and cytogenetic analysis. The repeated expansion of this genomic region in asexual P. antipodarum lineages following distinct transitions to asexuality represents a dramatic genome structural change associated with asexual reproduction—with potential functional consequences related to the loss of sexual reproduction.     

Red flag on the white reporter: a versatile insulator abuts the white gene in Drosophila and is omnipresent in mini-white constructs

Much of the research on insulators in Drosophila has been done with transgenic constructs using the white gene (mini-white) as reporter. Hereby we report that the sequence between the white and CG32795 genes in Drosophila melanogaster contains an insulator of a novel kind. Its functional core is within a 368 bp segment almost contiguous to the white 3′UTR, hence we name it as Wari (white-abutting resident insulator). Though Wari contains no binding sites for known insulator proteins and does not require Su(Hw) or Mod(mdg4) for its activity, it can equally well interact with another copy of Wari and with unrelated Su(Hw)-dependent insulators, gypsy or 1A2. In its natural downstream position, Wari reinforces enhancer blocking by any of the three insulators placed between the enhancer and the promoter; again, Wari–Wari, Wari–gypsy or 1A2–Wari pairing results in mutual neutralization (insulator bypass) when they precede the promoter. The distressing issue is that this element hides in all mini-white constructs employed worldwide to study various insulators and other regulatory elements as well as long-range genomic interactions, and its versatile effects could have seriously influenced the results and conclusions of many works.     

Genetic Analysis of the Dumpy Complex Locus in DROSOPHILA MELANOGASTER: Complementation, Fine Structure and Function

An extensive genetic analysis of the dumpy locus is presented. This study includes complementation, fine structure mapping and allelic interaction. A number of complementing recessive lethals of the dp complex have been genetically mapped. Two alleles of the ol(v) type that complement l alleles map to the left portion of the locus. A number of olv alleles that complement both l and lv lethals map within the right portion of the locus.——Fine-structure analysis demonstrated that both olv and o alleles are distributed among various subloci. Evidence for spacer regions between subloci is presented.——An extensive discussion of the data considers whether the locus is unicistronic or multicistronic. The conclusion reached is that the locus is not a single functional cistron. The possibility of a single cistron encoding a multifunctional polypeptide is discussed.——The hypothesis is proposed that the left portion of the map and the l mutations function as regulatory sequences and that the right portion of the map encodes structural sequences.     

A Comprehensive Map of Insulator Elements for the Drosophila Genome

Insulators are DNA sequences that control the interactions among genomic regulatory elements and act as chromatin boundaries. A thorough understanding of their location and function is necessary to address the complexities of metazoan gene regulation. We studied by ChIP–chip the genome-wide binding sites of 6 insulator-associated proteins—dCTCF, CP190, BEAF-32, Su(Hw), Mod(mdg4), and GAF—to obtain the first comprehensive map of insulator elements in Drosophila embryos. We identify over 14,000 putative insulators, including all classically defined insulators. We find two major classes of insulators defined by dCTCF/CP190/BEAF-32 and Su(Hw), respectively. Distributional analyses of insulators revealed that particular sub-classes of insulator elements are excluded between cis-regulatory elements and their target promoters; divide differentially expressed, alternative, and divergent promoters; act as chromatin boundaries; are associated with chromosomal breakpoints among species; and are embedded within active chromatin domains. Together, these results provide a map demarcating the boundaries of gene regulatory units and a framework for understanding insulator function during the development and evolution of Drosophila.     

<Emphasis Type="Italic">cis</Emphasis>-Decoder discovers constellations of conserved DNA sequences shared among tissue-specific enhancers

A systematic approach is described for analysis of evolutionarily conserved cis-regulatory DNA using cis-Decoder, a tool for discovery of conserved sequence elements that are shared between similarly regulated enhancers. Analysis of 2,086 conserved sequence blocks (CSBs), identified from 135 characterized enhancers, reveals most CSBs consist of shorter overlapping/adjacent elements that are either enhancer type-specific or common to enhancers with divergent regulatory behaviors. Our findings suggest that enhancers employ overlapping repertoires of highly conserved core elements.     

A novel method for tissue-specific RNAi rescue in Drosophila

Targeted gene silencing by RNA interference allows the study of gene function in plants and animals. In cell culture and small animal models, genetic screens can be performed—even tissue-specifically in Drosophila—with genome-wide RNAi libraries. However, a major problem with the use of RNAi approaches is the unavoidable false-positive error caused by off-target effects. Until now, this is minimized by computational RNAi design, comparing RNAi to the mutant phenotype if known, and rescue with a presumed ortholog. The ultimate proof of specificity would be to restore expression of the same gene product in vivo. Here, we present a simple and efficient method to rescue the RNAi-mediated knockdown of two independent genes in Drosophila. By exploiting the degenerate genetic code, we generated Drosophila RNAi Escape Strategy Construct (RESC) rescue proteins containing frequent silent mismatches in the complete RNAi target sequence. RESC products were no longer efficiently silenced by RNAi in cell culture and in vivo. As a proof of principle, we rescue the RNAi-induced loss of function phenotype of the eye color gene white and tracheal defects caused by the knockdown of the heparan sulfate proteoglycan syndecan. Our data suggest that RESC is widely applicable to rescue and validate ubiquitous or tissue-specific RNAi and to perform protein structure–function analysis.     

Effector caspase Dcp-1 and IAP protein Bruce regulate starvation-induced autophagy during Drosophila melanogaster oogenesis

A complex relationship exists between autophagy and apoptosis, but the regulatory mechanisms underlying their interactions are largely unknown. We conducted a systematic study of Drosophila melanogaster cell death–related genes to determine their requirement in the regulation of starvation-induced autophagy. We discovered that six cell death genes—death caspase-1 (Dcp-1), hid, Bruce, Buffy, debcl, and p53—as well as Ras–Raf–mitogen activated protein kinase signaling pathway components had a role in autophagy regulation in D. melanogaster cultured cells. During D. melanogaster oogenesis, we found that autophagy is induced at two nutrient status checkpoints: germarium and mid-oogenesis. At these two stages, the effector caspase Dcp-1 and the inhibitor of apoptosis protein Bruce function to regulate both autophagy and starvation-induced cell death. Mutations in Atg1 and Atg7 resulted in reduced DNA fragmentation in degenerating midstage egg chambers but did not appear to affect nuclear condensation, which indicates that autophagy contributes in part to cell death in the ovary. Our study provides new insights into the molecular mechanisms that coordinately regulate autophagic and apoptotic events in vivo.     

Identification of a Gene Regulating the Tissue Expression of a Phosphoglucomutase Locus in Rainbow Trout

Nine percent of the rainbow trout (Salmo gairdneri) from a hatchery source have a greater than 100-fold increase in expression of a phosphoglucomutase (PGM) locus, Pgm1, in the liver but have normal expression of this locus in other tissues. The results of genetic crosses are consistent with a single regulatory gene with additive inheritance being responsible for the differences in the amount of PGM activity in the liver.—The allele responsible for the expression of Pgm1 in the liver is apparently a recent mutation. This is supported by its restricted distribution in rainbow trout and the absence of liver Pgm1 expression in closely related species. This genetic system is valuable for future analysis of the control of gene expression and in determining the relative evolutionary importance of genetic variation at structural and regulatory genes.     

Single Nucleotide Polymorphisms with Cis-Regulatory Effects on Long Non-Coding Transcripts in Human Primary Monocytes

We applied genome-wide allele-specific expression analysis of monocytes from 188 samples. Monocytes were purified from white blood cells of healthy blood donors to detect cis-acting genetic variation that regulates the expression of long non-coding RNAs. We analysed 8929 regions harboring genes for potential long non-coding RNA that were retrieved from data from the ENCODE project. Of these regions, 60% were annotated as intergenic, which implies that they do not overlap with protein-coding genes. Focusing on the intergenic regions, and using stringent analysis of the allele-specific expression data, we detected robust cis-regulatory SNPs in 258 out of 489 informative intergenic regions included in the analysis. The cis-regulatory SNPs that were significantly associated with allele-specific expression of long non-coding RNAs were enriched to enhancer regions marked for active or bivalent, poised chromatin by histone modifications. Out of the lncRNA regions regulated by cis-acting regulatory SNPs, 20% (n = 52) were co-regulated with the closest protein coding gene. We compared the identified cis-regulatory SNPs with those in the catalog of SNPs identified by genome-wide association studies of human diseases and traits. This comparison identified 32 SNPs in loci from genome-wide association studies that displayed a strong association signal with allele-specific expression of non-coding RNAs in monocytes, with p-values ranging from 6.7×10⁻⁷ to 9.5×10⁻⁸⁹. The identified cis-regulatory SNPs are associated with diseases of the immune system, like multiple sclerosis and rheumatoid arthritis.     

Discovery and analysis of evolutionarily conserved intronic splicing regulatory elements

Knowledge of the functional cis-regulatory elements that regulate constitutive and alternative pre-mRNA splicing is fundamental for biology and medicine. Here we undertook a genome-wide comparative genomics approach using available mammalian genomes to identify conserved intronic splicing regulatory elements (ISREs). Our approach yielded 314 ISREs, and insertions of ~70 ISREs between competing splice sites demonstrated that 84% of ISREs altered 5′ and 94% altered 3′ splice site choice in human cells. Consistent with our experiments, comparisons of ISREs to known splicing regulatory elements revealed that 40%–45% of ISREs might have dual roles as exonic splicing silencers. Supporting a role for ISREs in alternative splicing, we found that 30%–50% of ISREs were enriched near alternatively spliced (AS) exons, and included almost all known binding sites of tissue-specific alternative splicing factors. Further, we observed that genes harboring ISRE-proximal exons have biases for tissue expression and molecular functions that are ISRE-specific. Finally, we discovered that for Nova1, neuronal PTB, hnRNP C, and FOX1, the most frequently occurring ISRE proximal to an alternative conserved exon in the splicing factor strongly resembled its own known RNA binding site, suggesting a novel application of ISRE density and the propensity for splicing factors to auto-regulate to associate RNA binding sites to splicing factors. Our results demonstrate that ISREs are crucial building blocks in understanding general and tissue-specific AS regulation and the biological pathways and functions regulated by these AS events.     

Function and Regulation of AUTS2, a Gene Implicated in Autism and Human Evolution

Nucleotide changes in the AUTS2 locus, some of which affect only noncoding regions, are associated with autism and other neurological disorders, including attention deficit hyperactivity disorder, epilepsy, dyslexia, motor delay, language delay, visual impairment, microcephaly, and alcohol consumption. In addition, AUTS2 contains the most significantly accelerated genomic region differentiating humans from Neanderthals, which is primarily composed of noncoding variants. However, the function and regulation of this gene remain largely unknown. To characterize auts2 function, we knocked it down in zebrafish, leading to a smaller head size, neuronal reduction, and decreased mobility. To characterize AUTS2 regulatory elements, we tested sequences for enhancer activity in zebrafish and mice. We identified 23 functional zebrafish enhancers, 10 of which were active in the brain. Our mouse enhancer assays characterized three mouse brain enhancers that overlap an ASD–associated deletion and four mouse enhancers that reside in regions implicated in human evolution, two of which are active in the brain. Combined, our results show that AUTS2 is important for neurodevelopment and expose candidate enhancer sequences in which nucleotide variation could lead to neurological disease and human-specific traits.