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.     

Genetic Analysis of ‘PAX6-Negative’ Individuals with Aniridia or Gillespie Syndrome

We report molecular genetic analysis of 42 affected individuals referred with a diagnosis of aniridia who previously screened as negative for intragenic PAX6 mutations. Of these 42, the diagnoses were 31 individuals with aniridia and 11 individuals referred with a diagnosis of Gillespie syndrome (iris hypoplasia, ataxia and mild to moderate developmental delay). Array-based comparative genomic hybridization identified six whole gene deletions: four encompassing PAX6 and two encompassing FOXC1. Six deletions with plausible cis-regulatory effects were identified: five that were 3ʹ (telomeric) to PAX6 and one within a gene desert 5ʹ (telomeric) to PITX2. Sequence analysis of the FOXC1 and PITX2 coding regions identified two plausibly pathogenic de novo FOXC1 missense mutations (p.Pro79Thr and p.Leu101Pro). No intragenic mutations were detected in PITX2. FISH mapping in an individual with Gillespie-like syndrome with an apparently balanced X;11 reciprocal translocation revealed disruption of a gene at each breakpoint: ARHGAP6 on the X chromosome and PHF21A on chromosome 11. In the other individuals with Gillespie syndrome no mutations were identified in either of these genes, or in HCCS which lies close to the Xp breakpoint. Disruption of PHF21A has previously been implicated in the causation of intellectual disability (but not aniridia). Plausibly causative mutations were identified in 15 out of 42 individuals (12/32 aniridia; 3/11 Gillespie syndrome). Fourteen of these mutations presented in the known aniridia genes; PAX6, FOXC1 and PITX2. The large number of individuals in the cohort with no mutation identified suggests greater locus heterogeneity may exist in both isolated and syndromic aniridia than was previously appreciated.     

Functional dissection of a napin gene promoter: identification of promoter elements required for embryo and endosperm-specific transcription

The promoter region (?309 to +44) of the Brassica napus storage protein gene napA was studied in transgenic tobacco by successive 5′ as well as internal deletions fused to the reporter gene GUS (β-glucuronidase). The expression in the two main tissues of the seed, the endosperm and the embryo, was shown to be differentially regulated. This tissue-specific regulation within the seed was found to affect the developmental expression during seed development. The region between ?309 to ?152, which has a large effect on quantitative expression, was shown to harbour four elements regulating embryo and one regulating endosperm expression. This region also displayed enhancer activity. Deletion of eight bp from position ?152 to position ?144 totally abolished the activity of the napA promoter. This deletion disrupted a cis element with similarity to an ABA-responsive element (ABRE) overlapping with an E-box, demonstrating its crucial importance for quantitative expression. An internal deletion of the region ?133 to ?120, resulted in increased activity in both leaves and endosperm and a decreased activity in the embryo. Within this region, a cis element similar to the (CA)_n element, found in other storage protein promoters, was identified. This suggest that the (CA)_n element is important for conferring seed specificity by serving both as an activator and a repressor element.     

Evolved Bmp6 enhancer alleles drive spatial shifts in gene expression during tooth development in sticklebacks

Mutations in enhancers have been shown to often underlie natural variation but the evolved differences in enhancer activity can be difficult to identify in vivo. Threespine sticklebacks (Gasterosteus aculeatus) are a robust system for studying enhancer evolution due to abundant natural genetic variation, a diversity of evolved phenotypes between ancestral marine and derived freshwater forms, and the tractability of transgenic techniques. Previous work identified a series of polymorphisms within an intronic enhancer of the Bone morphogenetic protein 6 (Bmp6) gene that are associated with evolved tooth gain, a derived increase in freshwater tooth number that arises late in development. Here, we use a bicistronic reporter construct containing a genetic insulator and a pair of reciprocal two-color transgenic reporter lines to compare enhancer activity of marine and freshwater alleles of this enhancer. In older fish, the two alleles drive partially overlapping expression in both mesenchyme and epithelium of developing teeth, but the freshwater enhancer drives a reduced mesenchymal domain and a larger epithelial domain relative to the marine enhancer. In younger fish, these spatial shifts in enhancer activity are less pronounced. Comparing Bmp6 expression by in situ hybridization in developing teeth of marine and freshwater fish reveals similar evolved spatial shifts in gene expression. Together, these data support a model in which the polymorphisms within this enhancer underlie evolved tooth gain by shifting the spatial expression of Bmp6 during tooth development, and provide a general strategy to identify spatial differences in enhancer activity in vivo.     

Multiple control elements are required for expression of the human CD34 gene

Two cis regulatory elements of the human CD34 gene, the promoter and a 3′ enhancer, have previously been described. In transient transfection assays, the promoter was not sufficient to direct cell type specific expression. In contrast, the 3′ enhancer was active only in CD34⁺ cell lines, suggesting that this element might be responsible for stem cell-restricted expression of the CD34 gene. In the current work, through deletion and transient transfection experiments, we delineated the core enhancer sequence. We examined the role of this element upon stable integration. Our data suggested the presence of additional control elements. In order to identify them, using DNaseI hypersensitivity and methylation studies, we determined the chromatin structure of the entire CD34 locus. Amongst a number of DNaseI hypersensitive sites, we detected a strong CD34⁺ cell type-specific site in intron 4. This region, however, did not work as an enhancer by itself. By analyzing stable transfectants and transgenic animals, we demonstrated that the 3′ enhancer and intron 4 hypersensitive regions, either alone or together, did not function as a locus control region upon chromosomal integration. In contrast, a 160 kb genomic fragment encompassing the entire CD34 gene contained regulatory elements sufficient for high-level CD34 mRNA expression in murine stable lines. Our data indicate that combinatorial action of multiple, proximal and long-range, cis elements is necessary for proper regulation of CD34 expression.     

A 12 kb multi-allelic copy number variation encompassing a GC gene enhancer is associated with mastitis resistance in dairy cattle

Clinical mastitis (CM) is an inflammatory disease occurring in the mammary glands of lactating cows. CM is under genetic control, and a prominent CM resistance QTL located on chromosome 6 was reported in various dairy cattle breeds. Nevertheless, the biological mechanism underpinning this QTL has been lacking. Herein, we mapped, fine-mapped, and discovered the putative causal variant underlying this CM resistance QTL in the Dutch dairy cattle population. We identified a ~12 kb multi-allelic copy number variant (CNV), that is in perfect linkage disequilibrium with a lead SNP, as a promising candidate variant. By implementing a fine-mapping and through expression QTL mapping, we showed that the group-specific component gene (GC), a gene encoding a vitamin D binding protein, is an excellent candidate causal gene for the QTL. The multiplicated alleles are associated with increased GC expression and low CM resistance. Ample evidence from functional genomics data supports the presence of an enhancer within this CNV, which would exert cis-regulatory effect on GC. We observed that strong positive selection swept the region near the CNV, and haplotypes associated with the multiplicated allele were strongly selected for. Moreover, the multiplicated allele showed pleiotropic effects for increased milk yield and reduced fertility, hinting that a shared underlying biology for these effects may revolve around the vitamin D pathway. These findings together suggest a putative causal variant of a CM resistance QTL, where a cis-regulatory element located within a CNV can alter gene expression and affect multiple economically important traits.     

Interspersed Homologous DNA of Autographa californica Nuclear Polyhedrosis Virus Enhances Delayed-Early Gene Expression

The five regions of homologous DNA which are interspersed in the genome of the baculovirus Autographa californica nuclear polyhedrosis virus increased the expression of a delayed-early gene of this virus. Although this activity was first observed as a 10-fold trans effect, the homologous region 5 (hr5) enhanced the expression of linked genes 1,000-fold. The hr5 enhancer also exhibited the other characteristics associated with viral enhancer elements, including orientation independence and the abilities to function at a distance from the linked promoter, to regulate heterologous promoters, and to increase the number of RNA polymerase molecules transcribing the linked genes. The expression of the immediate-early regulatory gene was not enhanced by cis-linked hr5, although the enhancer function may require the immediate-early regulatory gene product. The hr5 enhancer was relatively insensitive to competition by an excess of enhancer molecules. The nucleotide sequence of hr5 revealed two different conserved repeats separated by nonhomologous DNA. Deletion analysis of the hr5 enhancer indicated that a 30-base-pair inverted repeat was essential for enhancer function.     

PAX6 Regulates Melanogenesis in the Retinal Pigmented Epithelium through Feed-Forward Regulatory Interactions with MITF

During organogenesis, PAX6 is required for establishment of various progenitor subtypes within the central nervous system, eye and pancreas. PAX6 expression is maintained in a variety of cell types within each organ, although its role in each lineage and how it acquires cell-specific activity remain elusive. Herein, we aimed to determine the roles and the hierarchical organization of the PAX6-dependent gene regulatory network during the differentiation of the retinal pigmented epithelium (RPE). Somatic mutagenesis of Pax6 in the differentiating RPE revealed that PAX6 functions in a feed-forward regulatory loop with MITF during onset of melanogenesis. PAX6 both controls the expression of an RPE isoform of Mitf and synergizes with MITF to activate expression of genes involved in pigment biogenesis. This study exemplifies how one kernel gene pivotal in organ formation accomplishes a lineage-specific role during terminal differentiation of a single lineage.     

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.