Ancient exaptation of a CORE-SINE retroposon into a highly conserved mammalian neuronal enhancer of the proopiomelanocortin gene

The proopiomelanocortin gene (POMC) is expressed in the pituitary gland and the ventral hypothalamus of all jawed vertebrates, producing several bioactive peptides that function as peripheral hormones or central neuropeptides, respectively. We have recently determined that mouse and human POMC expression in the hypothalamus is conferred by the action of two 5′ distal and unrelated enhancers, nPE1 and nPE2. To investigate the evolutionary origin of the neuronal enhancer nPE2, we searched available vertebrate genome databases and determined that nPE2 is a highly conserved element in placentals, marsupials, and monotremes, whereas it is absent in nonmammalian vertebrates. Following an in silico paleogenomic strategy based on genome-wide searches for paralog sequences, we discovered that opossum and wallaby nPE2 sequences are highly similar to members of the superfamily of CORE-short interspersed nucleotide element (SINE) retroposons, in particular to MAR1 retroposons that are widely present in marsupial genomes. Thus, the neuronal enhancer nPE2 originated from the exaptation of a CORE-SINE retroposon in the lineage leading to mammals and remained under purifying selection in all mammalian orders for the last 170 million years. Expression studies performed in transgenic mice showed that two nonadjacent nPE2 subregions are essential to drive reporter gene expression into POMC hypothalamic neurons, providing the first functional example of an exapted enhancer derived from an ancient CORE-SINE retroposon. In addition, we found that this CORE-SINE family of retroposons is likely to still be active in American and Australian marsupial genomes and that several highly conserved exonic, intronic and intergenic sequences in the human genome originated from the exaptation of CORE-SINE retroposons. Together, our results provide clear evidence of the functional novelties that transposed elements contributed to their host genomes throughout evolution.     

Cis-regulatory characterization of sequence conservation surrounding the Hox4 genes

Hox genes are key regulators of anterior-posterior axis patterning and have a major role in hindbrain development. The zebrafish Hox4 paralogs have strong overlapping activities in hindbrain rhombomeres 7 and 8, in the spinal cord and in the pharyngeal arches. With the aim to predict enhancers that act on the hoxa4a, hoxb4a, hoxc4a and hoxd4a genes, we used sequence conservation around the Hox4 genes to analyze all fish:human conserved non-coding sequences by reporter assays in stable zebrafish transgenesis. Thirty-four elements were functionally tested in GFP reporter gene constructs and more than 100 F1 lines were analyzed to establish a correlation between sequence conservation and cis-regulatory function, constituting a catalog of Hox4 CNEs. Sixteen tissue-specific enhancers could be identified. Multiple alignments of the CNEs revealed paralogous cis-regulatory sequences, however, the CNE sequence similarities were found not to correlate with tissue specificity. To identify ancestral enhancers that direct Hox4 gene activity, genome sequence alignments of mammals, teleosts, horn shark and the cephalochordate amphioxus, which is the most basal extant chordate possessing a single prototypical Hox cluster, were performed. Three elements were identified and two of them exhibited regulatory activity in transgenic zebrafish, however revealing no specificity. Our data show that the approach to identify cis-regulatory sequences by genome sequence alignments and subsequent testing in zebrafish transgenesis can be used to define enhancers within the Hox clusters and that these have significantly diverged in their function during evolution.     

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.     

Hypothalamic‐specific proopiomelanocortin deficiency reduces alcohol drinking in male and female mice

Opioid receptor antagonist naltrexone reduces alcohol consumption and relapse in both humans and rodents. This study investigated whether hypothalamic proopiomelanocortin (POMC) neurons (producing beta‐endorphin and melanocortins) play a role in alcohol drinking behaviors. Both male and female mice with targeted deletion of two neuronal Pomc enhancers nPE1 and nPE2 (nPE?/?), resulting in hypothalamic‐specific POMC deficiency, were studied in short‐access (4‐h/day) drinking‐in‐the‐dark (DID, alcohol in one bottle, intermittent access (IA, 24‐h cycles of alcohol access every other day, alcohol vs. water in a two‐bottle choice) and alcohol deprivation effect (ADE) models. Wild‐type nPE+/+ exposed to 1‐week DID rapidly established stable alcohol drinking behavior with more intake in females, whereas nPE?/? mice of both sexes had less intake and less preference. Although nPE?/? showed less saccharin intake and preference than nPE+/+, there was no genotype difference in sucrose intake or preference in the DID paradigm. After 3‐week IA, nPE+/+ gradually escalated to high alcohol intake and preference, with more intake in females, whereas nPE?/? showed less escalation. Pharmacological blockade of mu‐opioid receptors with naltrexone reduced intake in nPE+/+ in a dose‐dependent manner, but had blunted effects in nPE?/? of both sexes. When alcohol was presented again after 1‐week abstinence from IA, nPE+/+ of both sexes displayed significant increases in alcohol intake (ADE or relapse‐like drinking), with more pronounced ADE in females, whereas nPE?/? did not show ADE in either sex. Our results suggest that neuronal POMC is involved in modulation of alcohol ‘binge’ drinking, escalation and ‘relapse’, probably via hypothalamic‐mediated mechanisms, with sex differences.     

Evolution of the vertebrate glucose-dependent insulinotropic polypeptide (GIP) gene

The glucose-dependent insulinotropic polypeptide (GIP) gene is believed to have originated from a gene duplication event very early in vertebrate evolution that also produced the proglucagon gene, yet so far GIP has only been described within mammals. Here we report the identification of GIP genes in chicken, frogs, and zebrafish. The chicken and frog genes are organized in a similar fashion to mammalian GIP genes and contain 6 exons and 5 introns in homologous locations. These genes can also potentially be proteolytically processed in identical patterns as observed in the mammalian sequences that would yield a GIP hormone that is only one amino shorter than the mammalian sequences due to the removal of an extra basic residue by carboxypeptidase E. The zebrafish GIP gene and precursor protein is shorter than other vertebrate GIP genes and is missing exon 5. The predicted zebrafish GIP hormone is also shorter, being only 31 amino acids in length. The zebrafish GIP hormone is similar in length to the proglucagon-derived peptide hormones, peptides encoded from the gene most closely related to GIP. We suggest that the structure of zebrafish GIP is more similar to the ancestral gene, and that tetrapod GIP has been extended. The mammalian GIP hormone has also undergone a period of rapid sequence evolution early in mammalian evolution. The discovery of a conserved GIP in diverse vertebrate suggests that it has an essential role in physiology in diverse vertebrates, although it may have only recently evolved a role as an incretin hormone.     

Discovery of a Novel Prolactin in Non-Mammalian Vertebrates: Evolutionary Perspectives and Its Involvement in Teleost Retina Development

Background

The three pituitary hormones, viz. prolactin (PRL), growth hormone (GH) and somatolactin (SL), together with the mammalian placental lactogen (PL), constitute a gene family of hormones with similar gene structure and encoded protein sequences. These hormones are believed to have evolved from a common ancestral gene through several rounds of gene duplication and subsequent divergence.

Principal Findings

In this study, we have identified a new PRL-like gene in non-mammalian vertebrates through bioinformatics and molecular cloning means. Phylogenetic analyses showed that this novel protein is homologous to the previously identified PRL. A receptor transactivation assay further showed that this novel protein could bind to PRL receptor to trigger the downstream post-receptor event, indicating that it is biologically active. In view of its close phylogenetic relationship with PRL and also its ability to activate PRL receptor, we name it as PRL2 and the previously identified PRL as PRL1. All the newly discovered PRL2 sequences possess three conserved disulfide linkages with the exception of the shark PRL2 which has only two. In sharp contrast to the classical PRL1 which is predominantly expressed in the pituitary, PRL2 was found to be mainly expressed in the eye and brain of the zebrafish but not in the pituitary. A largely reduced inner nuclear layer of the retina was observed after morpholino knockdown of zebrafish PRL2, indicating its role on retina development in teleost.

Significance

The discovery of this novel PRL has revitalized our understanding on the evolution of the GH/PRL/SL/PL gene family. Its unique expression and functions in the zebrafish eye also provide a new avenue of research on the neuroendocrine control of retina development in vertebrates.     

Conserved synteny of mammalian imprinted genes in chicken, frog, and fish genomes

Conservation of synteny of mammalian imprinted genes between chicken and human suggested that highly conserved gene clusters were selected long before these genes were recruited for genomic imprinting in mammals. Here we have applied in silico mapping of orthologous genes in pipid frog, zebrafish, spotted green and Japanese pufferfish to show considerable conservation of synteny in lower vertebrates. More than 400 million years ago in a common ancestor of teleost fish and tetrapods, 'preimprinted' chromosome regions homologous to human 6q25, 7q21, 7q32, 11p15, and 15q11-->q12 already contained most present-day mammalian imprinted genes. Interestingly, some imprinted gene orthologues which are isolated from imprinted clusters in mouse and human could be linked to preimprinted regions in lower vertebrates, indicating that separation occurred during mammalian evolution. On the contrary, newly arisen genes by segmental duplication in the mammalian lineage, i.e. SNRPN and FRAT3, were transposed or translocated to imprinted clusters and recruited for parent-specific activity. By analysis of currently available sequences of non-mammalian vertebrates, the imprinted gene clusters homologous to human chromosomes 14q32 and 19q12 are only poorly conserved in chicken, frog, and fish and, therefore, may not have evolved from ancestral preimprinted gene arrays. Evidently, evolution of imprinted gene clusters is an ongoing and dynamic process in mammals. In general, imprinted gene orthologues do not show a higher degree of synteny conservation in vertebrates than non-imprinted genes interspersed with or adjacent to an imprinted cluster.     

Identification of a dopaminergic enhancer indicates complexity in vertebrate dopamine neuron phenotype specification

The dopaminergic neurons of the basal ganglia play critical roles in CNS function and human disease, but specification of dopamine neuron phenotype is poorly understood in vertebrates. We performed an in vivo screen in zebrafish to identify dopaminergic neuron enhancers, in order to facilitate studies on the specification of neuronal identity, connectivity, and function in the basal ganglia. Based primarily on identification of conserved non-coding elements, we tested 54 DNA elements from four species (zebrafish, pufferfish, mouse, and rat), that included 21 genes with known or putative roles in dopaminergic neuron specification or function. Most elements failed to drive CNS expression or did not express specifically in dopaminergic neurons. However, we did isolate a discrete enhancer from the otpb gene that drove specific expression in diencephalic dopaminergic neurons, although it did not share sequence conservation with regulatory regions of otpa or other dopamine-specific genes. For the otpb enhancer, regulation of expression in dopamine neurons requires multiple elements spread across a large genomic area. In addition, we compared our in vivo testing with in silico analysis of genomic regions for genes involved in dopamine neuron function, but failed to find conserved regions that functioned as enhancers. We conclude that regulation of dopaminergic neuron phenotype in vertebrates is regulated by dispersed regulatory elements.     

Regulatory divergence of the duplicated chromosomal loci sox11a/b by subpartitioning and sequence evolution of enhancers in zebrafish

We used the classic example of the duplicated zebrafish sox11a and -b loci to test the duplication, degeneration, complementation (DDC) model of genome evolution through whole genome duplication. While recent reports have demonstrated sub-partitioning of regulatory sequences in duplicated regions, a comparison of the regulatory capabilities of extant regulatory sequences derived from ancient ancestral elements has been scarce. Consistent with the DDC model, we find that ancestral regulatory elements distributed over several hundred kb were lost in either one or the other zebrafish duplicate, leading to subpartitioning. However, regulatory sequences kept as duplicates near both sox11 co-orthologs diverged in sequence from each other and from human elements and in the regulatory patterns they drive in transgenic zebrafish. Evolutionary analysis of the loci suggested that both zebrafish protein coding sox11 orthologs have been maintained by purifying selection, and have evolved at comparable rates, indicative of non-diverged protein functions. The duplicated regulatory elements, conversely, evolved with different divergence rates and degrees of subfunctionalization. These data show that regulatory evolution of gene expression patterns occurred both through differential loss as well as through regulatory sequence evolution in zebrafish versus human genomes.     

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.     

Multiple Chromosomal Rearrangements Structured the Ancestral Vertebrate Hox-Bearing Protochromosomes

While the proposal that large-scale genome expansions occurred early in vertebrate evolution is widely accepted, the exact mechanisms of the expansion—such as a single or multiple rounds of whole genome duplication, bloc chromosome duplications, large-scale individual gene duplications, or some combination of these—is unclear. Gene families with a single invertebrate member but four vertebrate members, such as the Hox clusters, provided early support for Ohno's hypothesis that two rounds of genome duplication (the 2R-model) occurred in the stem lineage of extant vertebrates. However, despite extensive study, the duplication history of the Hox clusters has remained unclear, calling into question its usefulness in resolving the role of large-scale gene or genome duplications in early vertebrates. Here, we present a phylogenetic analysis of the vertebrate Hox clusters and several linked genes (the Hox “paralogon”) and show that different phylogenies are obtained for Dlx and Col genes than for Hox and ErbB genes. We show that these results are robust to errors in phylogenetic inference and suggest that these competing phylogenies can be resolved if two chromosomal crossover events occurred in the ancestral vertebrate. These results resolve conflicting data on the order of Hox gene duplications and the role of genome duplication in vertebrate evolution and suggest that a period of genome reorganization occurred after genome duplications in early vertebrates.     

Cooperation of sonic hedgehog enhancers in midline expression

In zebrafish, as in other vertebrates, the secreted signalling molecule Sonic hedgehog (Shh) is expressed in organiser regions such as the embryonic midline and the zona limitans intrathalamica (zli). To investigate the regulatory mechanisms underlying the pattern of shh expression, we carried out a systematic analysis of the intronic regulatory sequences of zebrafish shh using stable transgenesis. Deletion analysis identified the modules responsible for expression in the embryonic shield, the hypothalamus and the zli and confirmed the activities of previously identified notochord and floor plate enhancers. We detected a strong synergism between regulatory regions. The degree of synergy varied over time in the hypothalamus suggesting different mechanisms for initiation and maintenance of expression. Our data show that the pattern of shh expression in the embryonic central nervous system involves an intricate crosstalk of at least 4 different regulatory regions. When compared to the enhancer activities of the mouse Shh gene, we observed a remarkable divergence of function of structurally conserved enhancer sequences. The activating region ar-C (61% identical to SFPE2 in mouse Shh), for example, mediates floor plate expression in the mouse embryo while it directs expression in the forebrain and the notochord and only weakly in the floor plate in the zebrafish embryo. This raises doubts on the predictive power of phylogenetic footprinting and indicates a stunning divergence of function of structurally conserved regulatory modules during vertebrate evolution.     

BAC transgenic zebrafish reveal hypothalamic enhancer activity around obesity associated SNP rs9939609 within the human FTO gene

Single Nucleotide Polymorphisms in FTO intron 1 have been associated with obesity risk, leading to the hypothesis that FTO is the obesity‐related gene. However, other studies have shown that the FTO gene is part of the regulatory domain of the neighboring IRX3 gene and that enhancers in FTO intron 1 regulate IRX3. While Irx3 activity was shown to be necessary in the hypothalamus for the metabolic function of Irx3 in mouse, no enhancers with hypothalamic activity have been demonstrated in the risk‐associated region within FTO. In order to identify potential enhancers at the human FTO locus in vivo, we tested regulatory activity in FTO intron 1 using BAC transgenesis in zebrafish. A minimal gata2 promoter‐GFP cassette was inserted 1.3 kb upstream of the obesity associated SNP rs9939609 in a human FTO BAC plasmid. In addition to the previously identified expression domains in notochord and kidney, human FTO BAC:GFP transgenic zebrafish larvae expressed GFP in the ventral posterior tuberculum, the posterior hypothalamus and the anterior brainstem, which are also expression domains of zebrafish irx3a. In contrast, an in‐frame insertion of a GFP cassette at the FTO start codon resulted in weak ubiquitous GFP expression indicating that the promoter of FTO does likely not react to enhancers located in the obesity risk‐associated region. genesis 53:640–651, 2015. © 2015 The Authors. genesis Published by Wiley Periodicals, Inc.