Population genomics reveal recent speciation and rapid evolutionary adaptation in polar bears.

Polar bears are uniquely adapted to life in the High Arctic and have undergone drastic physiological changes in response to Arctic climates and a hyper-lipid diet of primarily marine mammal prey. We analyzed 89 complete genomes of polar bear and brown bear using population genomic modeling and show that the species diverged only 479-343 thousand years BP. We find that genes on the polar bear lineage have been under stronger positive selection than in brown bears; nine of the top 16 genes under strong positive selection are associated with cardiomyopathy and vascular disease, implying important reorganization of the cardiovascular system. One of the genes showing the strongest evidence of selection, APOB, encodes the primary lipoprotein component of low-density lipoprotein (LDL); functional mutations in APOB may explain how polar bears are able to cope with life-long elevated LDL levels that are associated with high risk of heart disease in humans.

Paternally Expressed Imprinted Genes under Positive Darwinian Selection in Arabidopsis thaliana.

Genomic imprinting is an epigenetic phenomenon where autosomal genes display uniparental expression depending on whether they are maternally or paternally inherited. Genomic imprinting can arise from parental conflicts over resource allocation to the offspring, which could drive imprinted loci to evolve by positive selection. We investigate whether positive selection is associated with genomic imprinting in the inbreeding species Arabidopsis thaliana. Our analysis of 140 genes regulated by genomic imprinting in the A. thaliana seed endosperm demonstrates they are evolving more rapidly than expected. To investigate whether positive selection drives this evolutionary acceleration, we identified orthologs of each imprinted gene across 34 plant species and elucidated their evolutionary trajectories. Increased positive selection was sought by comparing its incidence among imprinted genes with nonimprinted controls. Strikingly, we find a statistically significant enrichment of imprinted paternally expressed genes (iPEGs) evolving under positive selection, 50.6% of the total, but no such enrichment for positive selection among imprinted maternally expressed genes (iMEGs). This suggests that maternally- and paternally expressed imprinted genes are subject to different selective pressures. Almost all positively selected amino acids were fixed across 80 sequenced A. thaliana accessions, suggestive of selective sweeps in the A. thaliana lineage. The imprinted genes under positive selection are involved in processes important for seed development including auxin biosynthesis and epigenetic regulation. Our findings support a genomic imprinting model for plants where positive selection can affect paternally expressed genes due to continued conflict with maternal sporophyte tissues, even when parental conflict is reduced in predominantly inbreeding species.

Molecular evolution of mammalian genes with epistatic interactions in fertilization.

BACKGROUND: Genes that encode proteins associated with sperm competition, fertilization, and sexual conflicts of interest are often among the most rapidly evolving parts of animal genomes. One family of sperm-expressed genes (Zp3r, C4bpa) in the mammalian gene cluster called the regulator of complement activation (RCA) encodes proteins that bind eggs and mediate reproductive success, and are therefore expected to show high relative rates of nonsynonymous nucleotide substitution in response to sexual selection in comparison to other genes not involved in gamete binding at fertilization. We tested that working hypothesis by using phylogenetic models of codon evolution to identify episodes of diversifying positive selection. We used a comparative approach to quantify the evidence for episodic diversifying selection acting on RCA genes with known functions in fertilization (and sensitivity to sexual selection), and contrast them with other RCA genes in the same gene family that function in innate immunity (and are not sensitive to sexual selection). RESULTS: We expected but did not find evidence for more episodes of positive selection on Zp3r in Glires (the rodents and lagomorphs) or on C4BPA in Primates, in comparison to other paralogous RCA genes in the same taxon, or in comparison to the same orthologous RCA gene in the other taxon. That result was not unique to RCA genes: we also found little evidence for more episodes of diversifying selection on genes that encode selective sperm-binding molecules in the egg coat or zona pellucida (Zp2, Zp3) in comparison to members of the same gene family that encode structural elements of the egg coat (Zp1, Zp4). Similarly, we found little evidence for episodic diversifying selection acting on two other recently discovered genes (Juno, Izumo1) that encode essential molecules for sperm-egg fusion. CONCLUSIONS: These negative results help to illustrate the importance of a comparative context for this type of codon model analysis. The results may also point to other phylogenetic contexts in which the effects of selection acting on these fertilization proteins might be more readily discovered and documented in mammals and other taxa.

Adaptive Evolution as a Predictor of Species-Specific Innate Immune Response.

It has been proposed that positive selection may be associated with protein functional change. For example, human and macaque have different outcomes to HIV infection and it has been shown that residues under positive selection in the macaque TRIM5α receptor locate to the region known to influence species-specific response to HIV. In general, however, the relationship between sequence and function has proven difficult to fully elucidate, and it is the role of large-scale studies to help bridge this gap in our understanding by revealing major patterns in the data that correlate genotype with function or phenotype. In this study, we investigate the level of species-specific positive selection in innate immune genes from human and mouse. In total, we analyzed 456 innate immune genes using codon-based models of evolution, comparing human, mouse, and 19 other vertebrate species to identify putative species-specific positive selection. Then we used population genomic data from the recently completed Neanderthal genome project, the 1000 human genomes project, and the 17 laboratory mouse genomes project to determine whether the residues that were putatively positively selected are fixed or variable in these populations. We find evidence of species-specific positive selection on both the human and the mouse branches and we show that the classes of genes under positive selection cluster by function and by interaction. Data from this study provide us with targets to test the relationship between positive selection and protein function and ultimately to test the relationship between positive selection and discordant phenotypes.

Heterogeneous models place the root of the placental mammal phylogeny.

Heterogeneity among life traits in mammals has resulted in considerable phylogenetic conflict, particularly concerning the position of the placental root. Layered upon this are gene- and lineage-specific variation in amino acid substitution rates and compositional biases. Life trait variations that may impact upon mutational rates are longevity, metabolic rate, body size, and germ line generation time. Over the past 12 years, three main conflicting hypotheses have emerged for the placement of the placental root. These hypotheses place the Atlantogenata (common ancestor of Xenarthra plus Afrotheria), the Afrotheria, or the Xenarthra as the sister group to all other placental mammals. Model adequacy is critical for accurate tree reconstruction and by failing to account for these compositional and character exchange heterogeneities across the tree and data set, previous studies have not provided a strongly supported hypothesis for the placental root. For the first time, models that accommodate both tree and data set heterogeneity have been applied to mammal data. Here, we show the impact of accurate model assignment and the importance of data sets in accommodating model parameters while maintaining the power to reject competing hypotheses. Through these sophisticated methods, we demonstrate the importance of model adequacy, data set power and provide strong support for the Atlantogenata over other competing hypotheses for the position of the placental root.

Mitochondrial data are not suitable for resolving placental mammal phylogeny.

Mitochondrial data have traditionally been used in reconstructing a variety of species phylogenies. The low rates of recombination and thorough characterization of mitochondrial data across vertebrate species make it a particularly attractive phylogenetic marker. The relatively low number of fully sequenced mammal genomes and the lack of extensive sampling within Superorders have posed a serious problem for reaching agreement on the placement mammal species. The use of mitochondrial data sequences from large numbers of mammals could serve to circumvent the taxon-sampling deficit. Here we assess the suitability of mitochondrial data as a phylogenetic marker in mammal phylogenetics. MtDNA datasets of mammal origin have been filtered as follows: (i) we have sampled sparsely across the phylogenetic tree, (ii) we have constrained our sampling to genes with high taxon coverage, (iii) we have categorised rates across sites in a phylogeny independent manner and have removed fast evolving sites, and (iv), we have sampled from very shallow divergence times to reduce phylogenetic conflict. However, topologies obtained using these filters are not consistent with previous studies and are discordant across different genes. Individual mitochondrial genes, and indeed all mitochondrial genes analysed as a supermatrix, resulted in poor resolution of the species phylogeny. Overall, our study highlights the limitations of mitochondrial data, not only for resolving deep divergences and but also for shallow divergences in the mammal phylogeny.

Molecular adaptation of telomere associated genes in mammals.

BACKGROUND: Placental mammals display a huge range of life history traits, including size, longevity, metabolic rate and germ line generation time. Although a number of general trends have been proposed between these traits, there are exceptions that warrant further investigation. Species such as naked mole rat, human and certain bat species all exhibit extreme longevity with respect to body size. It has long been established that telomeres and telomere maintenance have a clear role in ageing but it has not yet been established whether there is evidence for adaptation in telomere maintenance proteins that could account for increased longevity in these species. RESULTS: Here we carry out a molecular investigation of selective pressure variation, specifically focusing on telomere associated genes across placental mammals. In general we observe a large number of instances of positive selection acting on telomere genes. Although these signatures of selection overall are not significantly correlated with either longevity or body size we do identify positive selection in the microbat species Myotis lucifugus in functionally important regions of the telomere maintenance genes DKC1 and TERT, and in naked mole rat in the DNA repair gene BRCA1. CONCLUSION: These results demonstrate the multifarious selective pressures acting across the mammal phylogeny driving lineage-specific adaptations of telomere associated genes. Our results show that regardless of the longevity of a species, these proteins have evolved under positive selection thereby removing increased longevity as the single selective force driving this rapid rate of evolution. However, evidence of molecular adaptations specific to naked mole rat and Myotis lucifugus highlight functionally significant regions in genes that may alter the way in which telomeres are regulated and maintained in these longer-lived species.

Colon cancer associated genes exhibit signatures of positive selection at functionally significant positions.

BACKGROUND: Cancer, much like most human disease, is routinely studied by utilizing model organisms. Of these model organisms, mice are often dominant. However, our assumptions of functional equivalence fail to consider the opportunity for divergence conferred by ~180 Million Years (MY) of independent evolution between these species. For a given set of human disease related genes, it is therefore important to determine if functional equivalency has been retained between species. In this study we test the hypothesis that cancer associated genes have different patterns of substitution akin to adaptive evolution in different mammal lineages. RESULTS: Our analysis of the current literature and colon cancer databases identified 22 genes exhibiting colon cancer associated germline mutations. We identified orthologs for these 22 genes across a set of high coverage (>6X) vertebrate genomes. Analysis of these orthologous datasets revealed significant levels of positive selection. Evidence of lineage-specific positive selection was identified in 14 genes in both ancestral and extant lineages. Lineage-specific positive selection was detected in the ancestral Euarchontoglires and Hominidae lineages for STK11, in the ancestral primate lineage for CDH1, in the ancestral Murinae lineage for both SDHC and MSH6 genes and the ancestral Muridae lineage for TSC1. CONCLUSION: Identifying positive selection in the Primate, Hominidae, Muridae and Murinae lineages suggests an ancestral functional shift in these genes between the rodent and primate lineages. Analyses such as this, combining evolutionary theory and predictions - along with medically relevant data, can thus provide us with important clues for modeling human diseases.

Positive selection neighboring functionally essential sites and disease-implicated regions of mammalian reproductive proteins.

BACKGROUND: Reproductive proteins are central to the continuation of all mammalian species. The evolution of these proteins has been greatly influenced by environmental pressures induced by pathogens, rival sperm, sexual selection and sexual conflict. Positive selection has been demonstrated in many of these proteins with particular focus on primate lineages. However, the mammalia are a diverse group in terms of mating habits, population sizes and germ line generation times. We have examined the selective pressures at work on a number of novel reproductive proteins across a wide variety of mammalia. RESULTS: We show that selective pressures on reproductive proteins are highly varied. Of the 10 genes analyzed in detail, all contain signatures of positive selection either across specific sites or in specific lineages or a combination of both. Our analysis of SP56 and Col1a1 are entirely novel and the results show positively selected sites present in each gene. Our findings for the Col1a1 gene are suggestive of a link between positive selection and severe disease type. We find evidence in our dataset to suggest that interacting proteins are evolving in symphony: most likely to maintain interacting functionality. CONCLUSION: Our in silico analyses show positively selected sites are occurring near catalytically important regions suggesting selective pressure to maximize efficient fertilization. In those cases where a mechanism of protein function is not fully understood, the sites presented here represent ideal candidates for mutational study. This work has highlighted the widespread rate heterogeneity in mutational rates across the mammalia and specifically has shown that the evolution of reproductive proteins is highly varied depending on the species and interacting partners. We have shown that positive selection and disease are closely linked in the Col1a1 gene.

Therapeutic potential of KLF2-induced exosomal microRNAs in pulmonary hypertension.

Pulmonary arterial hypertension (PAH) is a severe disorder of lung vasculature that causes right heart failure. Homoeostatic effects of flow-activated transcription factor Krüppel-like factor 2 (KLF2) are compromised in PAH. Here, we show that KLF2-induced exosomal microRNAs, miR-181a-5p and miR-324-5p act together to attenuate pulmonary vascular remodelling and that their actions are mediated by Notch4 and ETS1 and other key regulators of vascular homoeostasis. Expressions of KLF2, miR-181a-5p and miR-324-5p are reduced, while levels of their target genes are elevated in pre-clinical PAH, idiopathic PAH and heritable PAH with missense p.H288Y KLF2 mutation. Therapeutic supplementation of miR-181a-5p and miR-324-5p reduces proliferative and angiogenic responses in patient-derived cells and attenuates disease progression in PAH mice. This study shows that reduced KLF2 signalling is a common feature of human PAH and highlights the potential therapeutic role of KLF2-regulated exosomal miRNAs in PAH and other diseases associated with vascular remodelling.

Author Correction: Therapeutic potential of KLF2-induced exosomal microRNAs in pulmonary hypertension.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Human pancreatic islet three-dimensional chromatin architecture provides insights into the genetics of type 2 diabetes.

Genetic studies promise to provide insight into the molecular mechanisms underlying type 2 diabetes (T2D). Variants associated with T2D are often located in tissue-specific enhancer clusters or super-enhancers. So far, such domains have been defined through clustering of enhancers in linear genome maps rather than in three-dimensional (3D) space. Furthermore, their target genes are often unknown. We have created promoter capture Hi-C maps in human pancreatic islets. This linked diabetes-associated enhancers to their target genes, often located hundreds of kilobases away. It also revealed >1,300 groups of islet enhancers, super-enhancers and active promoters that form 3D hubs, some of which show coordinated glucose-dependent activity. We demonstrate that genetic variation in hubs impacts insulin secretion heritability, and show that hub annotations can be used for polygenic scores that predict T2D risk driven by islet regulatory variants. Human islet 3D chromatin architecture, therefore, provides a framework for interpretation of T2D genome-wide association study (GWAS) signals.

Functionally Conserved Noncoding Regulators of Cardiomyocyte Proliferation and Regeneration in Mouse and Human.

BACKGROUND: The adult mammalian heart has little regenerative capacity after myocardial infarction (MI), whereas neonatal mouse heart regenerates without scarring or dysfunction. However, the underlying pathways are poorly defined. We sought to derive insights into the pathways regulating neonatal development of the mouse heart and cardiac regeneration post-MI. METHODS AND RESULTS: Total RNA-seq of mouse heart through the first 10 days of postnatal life (referred to as P3, P5, P10) revealed a previously unobserved transition in microRNA (miRNA) expression between P3 and P5 associated specifically with altered expression of protein-coding genes on the focal adhesion pathway and cessation of cardiomyocyte cell division. We found profound changes in the coding and noncoding transcriptome after neonatal MI, with evidence of essentially complete healing by P10. Over two-thirds of each of the messenger RNAs, long noncoding RNAs, and miRNAs that were differentially expressed in the post-MI heart were differentially expressed during normal postnatal development, suggesting a common regulatory pathway for normal cardiac development and post-MI cardiac regeneration. We selected exemplars of miRNAs implicated in our data set as regulators of cardiomyocyte proliferation. Several of these showed evidence of a functional influence on mouse cardiomyocyte cell division. In addition, a subset of these miRNAs, miR-144-3p, miR-195a-5p, miR-451a, and miR-6240 showed evidence of functional conservation in human cardiomyocytes. CONCLUSIONS: The sets of messenger RNAs, miRNAs, and long noncoding RNAs that we report here merit further investigation as gatekeepers of cell division in the postnatal heart and as targets for extension of the period of cardiac regeneration beyond the neonatal period.

Re-analysis of public genetic data reveals a rare X-chromosomal variant associated with type 2 diabetes.

The reanalysis of existing GWAS data represents a powerful and cost-effective opportunity to gain insights into the genetics of complex diseases. By reanalyzing publicly available type 2 diabetes (T2D) genome-wide association studies (GWAS) data for 70,127 subjects, we identify seven novel associated regions, five driven by common variants (LYPLAL1, NEUROG3, CAMKK2, ABO, and GIP genes), one by a low-frequency (EHMT2), and one driven by a rare variant in chromosome Xq23, rs146662057, associated with a twofold increased risk for T2D in males. rs146662057 is located within an active enhancer associated with the expression of Angiotensin II Receptor type 2 gene (AGTR2), a modulator of insulin sensitivity, and exhibits allelic specific activity in muscle cells. Beyond providing insights into the genetics and pathophysiology of T2D, these results also underscore the value of reanalyzing publicly available data using novel genetic resources and analytical approaches.

Publisher Correction: Re-analysis of public genetic data reveals a rare X-chromosomal variant associated with type 2 diabetes.

In the originally published version of this Article, the affiliation details for Santi González, Jian'an Luan and Claudia Langenberg were inadvertently omitted. Santi González should have been affiliated with 'Barcelona Supercomputing Center (BSC), Joint BSC-CRG-IRB Research Program in Computational Biology, 08034 Barcelona, Spain', and Jian'an Luan and Claudia Langenberg should have been affiliated with 'MRC Epidemiology Unit, University of Cambridge School of Clinical Medicine, Cambridge Biomedical Campus, Cambridge CB2 0QQ, UK'. Furthermore, the abstract contained an error in the SNP ID for the rare variant in chromosome Xq23, which was incorrectly given as rs146662057 and should have been rs146662075. These errors have now been corrected in both the PDF and HTML versions of the Article.

Changes in the Coding and Non-coding Transcriptome and DNA Methylome that Define the Schwann Cell Repair Phenotype after Nerve Injury.

Repair Schwann cells play a critical role in orchestrating nerve repair after injury, but the cellular and molecular processes that generate them are poorly understood. Here, we perform a combined whole-genome, coding and non-coding RNA and CpG methylation study following nerve injury. We show that genes involved in the epithelial-mesenchymal transition are enriched in repair cells, and we identify several long non-coding RNAs in Schwann cells. We demonstrate that the AP-1 transcription factor C-JUN regulates the expression of certain micro RNAs in repair Schwann cells, in particular miR-21 and miR-34. Surprisingly, unlike during development, changes in CpG methylation are limited in injury, restricted to specific locations, such as enhancer regions of Schwann cell-specific genes (e.g., Nedd4l), and close to local enrichment of AP-1 motifs. These genetic and epigenomic changes broaden our mechanistic understanding of the formation of repair Schwann cell during peripheral nervous system tissue repair.