Massively parallel discovery of human-specific substitutions that alter enhancer activity.

Genetic changes that altered the function of gene regulatory elements have been implicated in the evolution of human traits such as the expansion of the cerebral cortex. However, identifying the particular changes that modified regulatory activity during human evolution remain challenging. Here we used massively parallel enhancer assays in neural stem cells to quantify the functional impact of >32,000 human-specific substitutions in >4,300 human accelerated regions (HARs) and human gain enhancers (HGEs), which include enhancers with novel activities in humans. We found that >30% of active HARs and HGEs exhibited differential activity between human and chimpanzee. We isolated the effects of human-specific substitutions from background genetic variation to identify the effects of genetic changes most relevant to human evolution. We found that substitutions interacted in both additive and nonadditive ways to modify enhancer function. Substitutions within HARs, which are highly constrained compared to HGEs, showed smaller effects on enhancer activity, suggesting that the impact of human-specific substitutions is buffered in enhancers with constrained ancestral functions. Our findings yield insight into how human-specific genetic changes altered enhancer function and provide a rich set of candidates for studies of regulatory evolution in humans.

Disrupting the three-dimensional regulatory topology of the Pitx1 locus results in overtly normal development.

Developmental gene expression patterns are orchestrated by thousands of distant-acting transcriptional enhancers. However, identifying enhancers essential for the expression of their target genes has proven challenging. Maps of long-range regulatory interactions may provide the means to identify enhancers crucial for developmental gene expression. To investigate this hypothesis, we used circular chromosome conformation capture coupled with interaction maps in the mouse limb to characterize the regulatory topology of Pitx1, which is essential for hindlimb development. We identified a robust hindlimb-specific interaction between Pitx1 and a putative hindlimb-specific enhancer. To interrogate the role of this interaction in Pitx1 regulation, we used genome editing to delete this enhancer in mouse. Although deletion of the enhancer completely disrupts the interaction, Pitx1 expression in the hindlimb is only mildly affected, without any detectable compensatory interactions between the Pitx1 promoter and potentially redundant enhancers. Pitx1 enhancer null mice did not exhibit any of the characteristic morphological defects of the Pitx1-/- mutant. Our results suggest that robust, tissue-specific physical interactions at essential developmental genes have limited predictive value for identifying enhancer mutations with strong loss-of-function phenotypes.

The genomic landscape of species divergence in Ficedula flycatchers.

Unravelling the genomic landscape of divergence between lineages is key to understanding speciation. The naturally hybridizing collared flycatcher and pied flycatcher are important avian speciation models that show pre- as well as postzygotic isolation. We sequenced and assembled the 1.1-Gb flycatcher genome, physically mapped the assembly to chromosomes using a low-density linkage map and re-sequenced population samples of each species. Here we show that the genomic landscape of species differentiation is highly heterogeneous with approximately 50 'divergence islands' showing up to 50-fold higher sequence divergence than the genomic background. These non-randomly distributed islands, with between one and three regions of elevated divergence per chromosome irrespective of chromosome size, are characterized by reduced levels of nucleotide diversity, skewed allele-frequency spectra, elevated levels of linkage disequilibrium and reduced proportions of shared polymorphisms in both species, indicative of parallel episodes of selection. Proximity of divergence peaks to genomic regions resistant to sequence assembly, potentially including centromeres and telomeres, indicate that complex repeat structures may drive species divergence. A much higher background level of species divergence of the Z chromosome, and a lower proportion of shared polymorphisms, indicate that sex chromosomes and autosomes are at different stages of speciation. This study provides a roadmap to the emerging field of speciation genomics.

Quantitative Mass Spectrometry Reveals Partial Translational Regulation for Dosage Compensation in Chicken.

There is increasing evidence that dosage compensation is not a ubiquitous feature following sex chromosome evolution, especially not in organisms where females are the heterogametic sex, like in birds. Even when it occurs, compensation can be incomplete and limited to dosage-sensitive genes. However, previous work has mainly studied transcriptional regulation of sex-linked genes, which may not reflect expression at the protein level. Here, we used liquid chromatography-tandem mass spectrometry to detect and quantify expressed levels of more than 2,400 proteins in ten different tissues of male and female chicken embryos. For comparison, transcriptome sequencing was performed in the same individuals, five of each sex. The proteomic analysis revealed that dosage compensation was incomplete, with a mean male-to-female (M:F) expression ratio of Z-linked genes of 1.32 across tissues, similar to that at the RNA level (1.29). The mean Z chromosome-to-autosome expression ratio was close to 1 in males and lower than 1 in females, consistent with partly reduced Z chromosome expression in females. Although our results exclude a general mechanism for chromosome-wide dosage compensation at translation, 30% of all proteins encoded from Z-linked genes showed a significant change in the M:F ratio compared with the corresponding ratio at the RNA level. This resulted in a pattern where some genes showed balanced expression between sexes and some close to 2-fold higher expression in males. This suggests that proteomic analyses will be necessary to reveal a more complete picture of gene regulation and sex chromosome evolution.

Divergence in gene expression within and between two closely related flycatcher species.

Relatively little is known about the character of gene expression evolution as species diverge. It is for instance unclear if gene expression generally evolves in a clock-like manner (by stabilizing selection or neutral evolution) or if there are frequent episodes of directional selection. To gain insights into the evolutionary divergence of gene expression, we sequenced and compared the transcriptomes of multiple organs from population samples of collared (Ficedula albicollis) and pied flycatchers (F. hypoleuca), two species which diverged less than one million years ago. Ordination analysis separated samples by organ rather than by species. Organs differed in their degrees of expression variance within species and expression divergence between species. Variance was negatively correlated with expression breadth and protein interactivity, suggesting that pleiotropic constraints reduce gene expression variance within species. Variance was correlated with between-species divergence, consistent with a pattern expected from stabilizing selection and neutral evolution. Using an expression PST approach, we identified genes differentially expressed between species and found 16 genes uniquely expressed in one of the species. For one of these, DPP7, uniquely expressed in collared flycatcher, the absence of expression in pied flycatcher could be associated with a ≈20-kb deletion including 11 of 13 exons. This study of a young vertebrate speciation model system expands our knowledge of how gene expression evolves as natural populations become reproductively isolated.

CpG island turnover events predict evolutionary changes in enhancer activity.

BACKGROUND: Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. CpG islands (CGIs) have recently been shown to influence enhancer activity, and here we test how their turnover across species contributes to enhancer evolution. RESULTS: We integrate maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and find that CGI content in enhancers is strongly associated with increased histone modification levels. CGIs show widespread turnover across species and species-specific CGIs are strongly enriched for enhancers exhibiting species-specific activity across all tissues and species. Genes associated with enhancers with species-specific CGIs show concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon. CONCLUSIONS: Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

Massively parallel disruption of enhancers active in human neural stem cells.

Changes in gene regulation have been linked to the expansion of the human cerebral cortex and to neurodevelopmental disorders, potentially by altering neural progenitor proliferation. However, the effects of genetic variation within regulatory elements on neural progenitors remain obscure. We use sgRNA-Cas9 screens in human neural stem cells (hNSCs) to disrupt 10,674 genes and 26,385 conserved regions in 2,227 enhancers active in the developing human cortex and determine effects on proliferation. Genes with proliferation phenotypes are associated with neurodevelopmental disorders and show biased expression in specific fetal human brain neural progenitor populations. Although enhancer disruptions overall have weaker effects than gene disruptions, we identify enhancer disruptions that severely alter hNSC self-renewal. Disruptions in human accelerated regions, implicated in human brain evolution, also alter proliferation. Integrating proliferation phenotypes with chromatin interactions reveals regulatory relationships between enhancers and their target genes contributing to neurogenesis and potentially to human cortical evolution.

Evolutionary innovation in conserved regulatory elements across the mammalian tree of life.

Transcriptional enhancers orchestrate cell type- and time point-specific gene expression programs. Evolution of enhancer sequences can alter target gene expression without causing detrimental misexpression in other contexts. It has long been thought that this modularity allows evolutionary changes in enhancers to escape pleiotropic constraints, which is especially important for evolutionary constrained developmental patterning genes. However, there is still little data supporting this hypothesis. Here we identified signatures of accelerated evolution in conserved enhancer elements across the mammalian phylogeny. We found that pleiotropic genes involved in gene regulatory and developmental processes were enriched for accelerated sequence evolution within their enhancer elements. These genes were associated with an excess number of enhancers compared to other genes, and due to this they exhibit a substantial degree of sequence acceleration over all their enhancers combined. We provide evidence that sequence acceleration is associated with turnover of regulatory function. We studied one acceleration event in depth and found that its sequence evolution led to the emergence of a new enhancer activity domain that may be involved in the evolution of digit reduction in hoofed mammals. Our results provide tangible evidence that enhancer evolution has been a frequent contributor to modifications involving constrained developmental signaling genes in mammals.

CpG island turnover events predict evolutionary changes in enhancer activity.

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. Here we show that turnover of CpG islands (CGIs), which contribute to enhancer activation, is broadly associated with changes in enhancer activity across mammals, including humans. We integrated maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and found that CGI content in enhancers was strongly associated with increased histone modification levels. CGIs showed widespread turnover across species and species-specific CGIs were strongly enriched for enhancers exhibiting species-specific activity across all tissues and species we examined. Genes associated with enhancers with species-specific CGIs showed concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon. Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

A systems biology approach identifies the role of dysregulated PRDM6 in the development of hypertension.

Genetic variants in the third intron of the PRDM6 gene have been associated with BP traits in multiple GWAS. By combining fine mapping, massively parallel reporter assays, and gene editing, we identified super enhancers that drive the expression of PRDM6 and are partly regulated by STAT1 as the causal variants for hypertension. The heterozygous disruption of Prdm6 in mice expressing Cre recombinase under the control of mouse smooth muscle cell protein 22-α promoter (Prdm6fl/+ SM22-Cre) exhibited a markedly higher number of renin-producing cells in the kidneys at E18.5 compared with WT littermates and developed salt-induced systemic hypertension that was completely responsive to the renin inhibitor aliskiren. Strikingly, RNA-Seq analysis of the mouse aortas identified a network of PRDM6-regulated genes that are located in GWAS-associated loci for blood pressure, most notably Sox6, which modulates renin expression in the kidney. Accordingly, the smooth muscle cell-specific disruption of Sox6 in Prdm6fl/+ SM22-Cre mice resulted in a dramatic reduction of renin. Fate mapping and histological studies also showed increased numbers of neural crest-derived cells accompanied by increased collagen deposition in the kidneys of Prdm6fl/+ Wnt1Cre-ZsGreen1Cre mice compared with WT mice. These findings establish the role of PRDM6 as a regulator of renin-producing cell differentiation into smooth muscle cells and as an attractive target for the development of antihypertensive drugs.

Sequencing of the needle transcriptome from Norway spruce (Picea abies Karst L.) reveals lower substitution rates, but similar selective constraints in gymnosperms and angiosperms.

BACKGROUND: A detailed knowledge about spatial and temporal gene expression is important for understanding both the function of genes and their evolution. For the vast majority of species, transcriptomes are still largely uncharacterized and even in those where substantial information is available it is often in the form of partially sequenced transcriptomes. With the development of next generation sequencing, a single experiment can now simultaneously identify the transcribed part of a species genome and estimate levels of gene expression. RESULTS: mRNA from actively growing needles of Norway spruce (Picea abies) was sequenced using next generation sequencing technology. In total, close to 70 million fragments with a length of 76 bp were sequenced resulting in 5 Gbp of raw data. A de novo assembly of these reads, together with publicly available expressed sequence tag (EST) data from Norway spruce, was used to create a reference transcriptome. Of the 38,419 PUTs (putative unique transcripts) longer than 150 bp in this reference assembly, 83.5% show similarity to ESTs from other spruce species and of the remaining PUTs, 3,704 show similarity to protein sequences from other plant species, leaving 4,167 PUTs with limited similarity to currently available plant proteins. By predicting coding frames and comparing not only the Norway spruce PUTs, but also PUTs from the close relatives Picea glauca and Picea sitchensis to both Pinus taeda and Taxus mairei, we obtained estimates of synonymous and non-synonymous divergence among conifer species. In addition, we detected close to 15,000 SNPs of high quality and estimated gene expression differences between samples collected under dark and light conditions. CONCLUSIONS: Our study yielded a large number of single nucleotide polymorphisms as well as estimates of gene expression on transcriptome scale. In agreement with a recent study we find that the synonymous substitution rate per year (0.6 × 10-09 and 1.1 × 10-09) is an order of magnitude smaller than values reported for angiosperm herbs. However, if one takes generation time into account, most of this difference disappears. The estimates of the dN/dS ratio (non-synonymous over synonymous divergence) reported here are in general much lower than 1 and only a few genes showed a ratio larger than 1.

Modeling uniquely human gene regulatory function via targeted humanization of the mouse genome.

The evolution of uniquely human traits likely entailed changes in developmental gene regulation. Human Accelerated Regions (HARs), which include transcriptional enhancers harboring a significant excess of human-specific sequence changes, are leading candidates for driving gene regulatory modifications in human development. However, insight into whether HARs alter the level, distribution, and timing of endogenous gene expression remains limited. We examined the role of the HAR HACNS1 (HAR2) in human evolution by interrogating its molecular functions in a genetically humanized mouse model. We find that HACNS1 maintains its human-specific enhancer activity in the mouse embryo and modifies expression of Gbx2, which encodes a transcription factor, during limb development. Using single-cell RNA-sequencing, we demonstrate that Gbx2 is upregulated in the limb chondrogenic mesenchyme of HACNS1 homozygous embryos, supporting that HACNS1 alters gene expression in cell types involved in skeletal patterning. Our findings illustrate that humanized mouse models provide mechanistic insight into how HARs modified gene expression in human evolution.

Transcriptome-Wide Patterns of the Genetic and Expression Variations in Two Sympatric Schizothoracine Fishes in a Tibetan Plateau Glacier Lake.

Sympatric speciation remains a central focus of evolutionary biology. Although some evidence shows speciation occurring in this way, little is known about the gene expression evolution and the characteristics of population genetics as species diverge. Two closely related Gymnocypris fish (Gymnocypris chui and Gymnocypris scleracanthus), which come from a small glacier lake in the Tibetan Plateau, Lake Langcuo, exist a possible incipient sympatric adaptive ecological speciation. We generated large amounts of RNA-Seq data from multiple individuals and tissues from each of the two species and compared gene expression patterns and genetic polymorphisms between them. Ordination analysis separated samples by organ rather than by species. The degree of expression difference between organs within and between species was different. Phylogenetic analyses indicated that the two closely related taxa formed a monophyletic complex. Population structure analysis displayed two distinctly divergent clusters of G. chui and G. scleracanthus populations. By contrast, G. scleracanthus population genetic diversity is higher than that of G. chui. Considerable sites of the two populations were differentiated with a coefficient of FST = 0.25-0.50, implying that a small proportion of loci nevertheless exhibited deep divergence in two comparisons. Concomitantly, putatively selected genes during speciation revealed functional categories are enriched in bone morphogenesis, cell growth, neurogenetics, enzyme activity, and binding activity in G. chui population. In contrast, nutrition and localization were highlighted in G. scleracanthus. Collectively, morphological traits and dietary preference combine with genetic variation and expression variation, probably contributed to the incipient speciation of two sympatric populations.

RPASE: Individual-based allele-specific expression detection without prior knowledge of haplotype phase.

Variation in gene expression is believed to make a significant contribution to phenotypic diversity and divergence. The analysis of allele-specific expression (ASE) can reveal important insights into gene expression regulation. We developed a novel method called RPASE (Read-backed Phasing-based ASE detection) to test for genes that show ASE. With mapped RNA-seq data from a single individual and a list of SNPs from the same individual as the only input, RPASE is capable of aggregating information across multiple dependent SNPs and producing individual-based gene-level tests for ASE. RPASE performs well in simulations and comparisons. We applied RPASE to multiple bird species and found a potentially rich landscape of ASE.

Bayesian Inference of Allele-Specific Gene Expression Indicates Abundant Cis-Regulatory Variation in Natural Flycatcher Populations.

Polymorphism in cis-regulatory sequences can lead to different levels of expression for the two alleles of a gene, providing a starting point for the evolution of gene expression. Little is known about the genome-wide abundance of genetic variation in gene regulation in natural populations but analysis of allele-specific expression (ASE) provides a means for investigating such variation. We performed RNA-seq of multiple tissues from population samples of two closely related flycatcher species and developed a Bayesian algorithm that maximizes data usage by borrowing information from the whole data set and combines several SNPs per transcript to detect ASE. Of 2,576 transcripts analyzed in collared flycatcher, ASE was detected in 185 (7.2%) and a similar frequency was seen in the pied flycatcher. Transcripts with statistically significant ASE commonly showed the major allele in >90% of the reads, reflecting that power was highest when expression was heavily biased toward one of the alleles. This would suggest that the observed frequencies of ASE likely are underestimates. The proportion of ASE transcripts varied among tissues, being lowest in testis and highest in muscle. Individuals often showed ASE of particular transcripts in more than one tissue (73.4%), consistent with a genetic basis for regulation of gene expression. The results suggest that genetic variation in regulatory sequences commonly affects gene expression in natural populations and that it provides a seedbed for phenotypic evolution via divergence in gene expression.

Evolutionary analysis of the female-specific avian W chromosome.

The typically repetitive nature of the sex-limited chromosome means that it is often excluded from or poorly covered in genome assemblies, hindering studies of evolutionary and population genomic processes in non-recombining chromosomes. Here, we present a draft assembly of the non-recombining region of the collared flycatcher W chromosome, containing 46 genes without evidence of female-specific functional differentiation. Survival of genes during W chromosome degeneration has been highly non-random and expression data suggest that this can be attributed to selection for maintaining gene dose and ancestral expression levels of essential genes. Re-sequencing of large population samples revealed dramatically reduced levels of within-species diversity and elevated rates of between-species differentiation (lineage sorting), consistent with low effective population size. Concordance between W chromosome and mitochondrial DNA phylogenetic trees demonstrates evolutionary stable matrilineal inheritance of this nuclear-cytonuclear pair of chromosomes. Our results show both commonalities and differences between W chromosome and Y chromosome evolution.

Genomic identification and characterization of the pseudoautosomal region in highly differentiated avian sex chromosomes.

The molecular characteristics of the pseudoautosomal region (PAR) of sex chromosomes remain elusive. Despite significant genome-sequencing efforts, the PAR of highly differentiated avian sex chromosomes remains to be identified. Here we use linkage analysis together with whole-genome re-sequencing to uncover the 630-kb PAR of an ecological model species, the collared flycatcher. The PAR contains 22 protein-coding genes and is GC rich. The genetic length is 64 cM in female meiosis, consistent with an obligate crossing-over event. Recombination is concentrated to a hotspot region, with an extreme rate of >700 cM/Mb in a 67-kb segment. We find no signatures of sexual antagonism and propose that sexual antagonism may have limited influence on PAR sequences when sex chromosomes are nearly fully differentiated and when a recombination hotspot region is located close to the PAR boundary. Our results demonstrate that a very small PAR suffices to ensure homologous recombination and proper segregation of sex chromosomes during meiosis.

Transcriptome sequencing reveals the character of incomplete dosage compensation across multiple tissues in flycatchers.

Sex chromosome divergence, which follows the cessation of recombination and degeneration of the sex-limited chromosome, can cause a reduction in expression level for sex-linked genes in the heterozygous sex, unless some mechanisms of dosage compensation develops to counter the reduction in gene dose. Because large-scale perturbations in expression levels arising from changes in gene dose might have strong deleterious effects, the evolutionary response should be strong. However, in birds and in at least some other female heterogametic organisms, wholesale sex chromosome dosage compensation does not seem to occur. Using RNA-seq of multiple tissues and individuals, we investigated male and female expression levels of Z-linked and autosomal genes in the collared flycatcher, a bird for which a draft genome sequence recently has been reported. We found that male expression of Z-linked genes was on average 50% higher than female expression, although there was considerable variation in the male-to-female ratio among genes. The ratio for individual genes was well correlated among tissues and there was also a correlation in the extent of compensation between flycatcher and chicken orthologs. The relative excess of male expression was positively correlated with expression breadth, expression level, and number of interacting proteins (protein connectivity), and negatively correlated with variance in expression. These observations lead to a model of compensation occurring on a gene-by-gene basis, supported by an absence of clustering of genes on the Z chromosome with respect to the extent of compensation. Equal mean expression level of autosomal and Z-linked genes in males, and 50% higher expression of autosomal than Z-linked genes in females, is compatible with that partial compensation is achieved by hypertranscription from females' single Z chromosome. A comparison with male-to-female expression ratios in orthologous Z-linked genes of ostriches, where Z-W recombination still occurs, suggests that male-biased expression of Z-linked genes is a derived trait after avian sex chromosome divergence.