Evolutionary Processes Acting on Candidate cis-Regulatory Regions in Humans Inferred from Patterns of Polymorphism and Divergence

Analysis of polymorphism and divergence in the non-coding portion of the human genome yields crucial information about factors driving the evolution of gene regulation. Candidate cis-regulatory regions spanning more than 15,000 genes in 15 African Americans and 20 European Americans were re-sequenced and aligned to the chimpanzee genome in order to identify potentially functional polymorphism and to characterize and quantify departures from neutral evolution. Distortions of the site frequency spectra suggest a general pattern of selective constraint on conserved non-coding sites in the flanking regions of genes (CNCs). Moreover, there is an excess of fixed differences that cannot be explained by a Gamma model of deleterious fitness effects, suggesting the presence of positive selection on CNCs. Extensions of the McDonald-Kreitman test identified candidate cis-regulatory regions with high probabilities of positive and negative selection near many known human genes, the biological characteristics of which exhibit genome-wide trends that differ from patterns observed in protein-coding regions. Notably, there is a higher probability of positive selection in candidate cis-regulatory regions near genes expressed in the fetal brain, suggesting that a larger portion of adaptive regulatory changes has occurred in genes expressed during brain development. Overall we find that natural selection has played an important role in the evolution of candidate cis-regulatory regions throughout hominid evolution.     

Developmentally programmed DNA rearrangement in Tetrahymena thermophila: Isolation and sequence characterization of three new alternative deletion systems

Extensive developmentally programmed DNA rearrangements, including thousands of internal deletions, occur in the differentiating somatic macronucleus in Tetrahymena thermophila. Some deletion systems involve the use of multiple alternative deletion sites. We report here the cloning and the sequences of three new alternative deletion systems (RR, RP and B) obtained using genomic subtraction. The RP and RR deletion systems are 2 kb apart on chromosome IR, and both involve the removal of < 2 kb of micronuclear sequences. The B deletion system is on chromosome 5 and involves a deletion of > 5 kb. All three deleted regions are very AT rich (∼ 80%) and do not appear to encode any protein. Sequences of the regions flanking the deletion junctions of all three systems revealed no sequence similarity among them nor with any previously reported deletion systems, suggesting that different cis-acting elements are involved for rearrangement. Unlike other deletion systems in ciliates, the B deletion system lacks short terminal direct repeats. Our results suggest an average of at least one alternative deletion system per 134 kb of micronuclear DNA and lead to an estimate that at least 25% of all deletion systems in Tetrahymena utilize alternative ends. The genomic subtraction method employed in this study could prove useful for the isolation of alternatively deleted DNA in special-purpose cases in Tetrahymena and other ciliates. The hybridization parameters for genomic subtraction worked out here for highly AT-rich DNA may have wider usefulness.     

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.     

A survey of ancient conserved non-coding elements in the PAX6 locus reveals a landscape of interdigitated cis-regulatory archipelagos

Biological differences between cell types and developmental processes are characterised by differences in gene expression profiles. Gene-distal enhancers are key components of the regulatory networks that specify the tissue-specific expression patterns driving embryonic development and cell fate decisions, and variations in their sequences are a major contributor to genetic disease and disease susceptibility. Despite advances in the methods for discovery of putative cis-regulatory sequences, characterisation of their spatio-temporal enhancer activities in a mammalian model system remains a major bottle-neck. We employed a strategy that combines gnathostome sequence conservation with transgenic mouse and zebrafish reporter assays to survey the genomic locus of the developmental control gene PAX6 for the presence of novel cis-regulatory elements. Sequence comparison between human and the cartilaginous elephant shark (Callorhinchus milii) revealed several ancient gnathostome conserved non-coding elements (agCNEs) dispersed widely throughout the PAX6 locus, extending the range of the known PAX6 cis-regulatory landscape to contain the full upstream PAX6-RCN1 intergenic region. Our data indicates that ancient conserved regulatory sequences can be tested effectively in transgenic zebrafish even when not conserved in zebrafish themselves. The strategy also allows efficient dissection of compound regulatory regions previously assessed in transgenic mice. Remarkable overlap in expression patterns driven by sets of agCNEs indicates that PAX6 resides in a landscape of multiple tissue-specific regulatory archipelagos.     

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.     

A developmentally regulated deletion element with long terminal repeats has cis-acting sequences in the flanking DNA

Approximately 6000 specific DNA deletion events occur during development of the somatic macronucleus of the ciliate Tetrahymena. The eliminated Tlr1 element is 13 kb or more in length and has an 825 bp inverted repeat near the rearrangement junctions. A functional analysis of the cis-acting sequences required for Tlr1 rearrangement was performed. A construct consisting of the entire inverted repeat and several hundred base pairs of flanking DNA on each side was rearranged accurately in vivo and displayed junctional variability similar to the chromosomal Tlr1 rearrangement. Thus, 11 kb or more of internal element DNA is not required in cis for DNA rearrangement. A second construct with only 51 bp of Tetrahymena DNA flanking the right junction underwent aberrant rearrangement. Thus, a signal for determination of the Tlr1 junction is located in the flanking DNA, 51 bp or more from the right junction. Within the Tlr1 inverted repeat are 19 bp tandem repeats. A construct with the 19mer repeat region deleted from the right half of the inverted repeat utilized normal rearrangement junctions. Thus, despite its transposon-like structure, Tlr1 is similar to other DNA rearrangements in Tetrahymena in possessing cis-acting sequences outside the deleted DNA.     

Disruption of Autoregulatory Feedback by a Mutation in a Remote,Ultraconserved PAX6 Enhancer Causes Aniridia

The strictly regulated expression of most pleiotropic developmental control genes is critically dependent on the activity of long-range cis-regulatory elements. This was revealed by the identification of individuals with a genetic condition lacking coding-region mutations in the gene commonly associated with the disease but having a variety of nearby chromosomal abnormalities, collectively described as cis-ruption disease cases. The congenital eye malformation aniridia is caused by haploinsufficiency of the developmental regulator PAX6. We discovered a de novo point mutation in an ultraconserved cis-element located 150 kb downstream from PAX6 in an affected individual with intact coding region and chromosomal locus. The element SIMO acts as a strong enhancer in developing ocular structures. The mutation disrupts an autoregulatory PAX6 binding site, causing loss of enhancer activity, resulting in defective maintenance of PAX6 expression. These findings reveal a distinct regulatory mechanism for genetic disease by disruption of an autoregulatory feedback loop critical for maintenance of gene expression through development.