Evolutionary history and adaptation from high-coverage whole-genome sequences of diverse african hunter-gatherers

To reconstruct modern human evolutionary history and identify loci that have shaped hunter-gatherer adaptation, we sequenced the whole genomes of five individuals in each of three different hunter-gatherer populations at >60× coverage: Pygmies from Cameroon and Khoesan-speaking Hadza and Sandawe from Tanzania. We identify 13.4 million variants, substantially increasing the set of known human variation. We found evidence of archaic introgression in all three populations, and the distribution of time to most recent common ancestors from these regions is similar to that observed for introgressed regions in Europeans. Additionally, we identify numerous loci that harbor signatures of local adaptation, including genes involved in immunity, metabolism, olfactory and taste perception, reproduction, and wound healing. Within the Pygmy population, we identify multiple highly differentiated loci that play a role in growth and anterior pituitary function and are associated with height.     

Effect of Heat Stress on Expression of DILP2 and DILP3 Insulin-Like Peptide Genes in <Emphasis Type="Italic">Drosophila melanogaster</Emphasis> Adults

In a number of works, it was demonstrated that insulin/insulin-like growth factor signaling pathway in the Drosophila melanogaster can be involved in the control of the organism reaction to stress. However, it remains unclear which of eight insulin-like peptides (ligands of insulin/insulin-like growth factor signaling pathway) known in the D. melanogaster are involved in the response to different types of stress. We conducted immunohistochemical analysis of the expression of two insulin-like peptide genes (DILP2 and DILP3) in insulin-producing cells of the brain in adult D. melanogaster females after heat stress. We for the first time found that the DILP3 is one of the components of the response to heat stress, while the DILP2 is apparently not involved in the organism response to heat stress.     

Adaptive evolution of loci covarying with the human African Pygmy phenotype

African Pygmies are hunter-gatherer populations from the equatorial rainforest that present the lowest height averages among humans. The biological basis and the putative adaptive role of the short stature of Pygmy populations has been one of the most intriguing topics for human biologists in the last century, which still remains elusive. Worldwide convergent evolution of the Pygmy size suggests the presence of strong selective pressures on the phenotype. We developed a novel approach to survey the genetic architecture of phenotypes and applied it to study the genomic covariation between allele frequencies and height measurements among Pygmy and non-Pygmy populations. Among the regions that were most associated with the phenotype, we identified a significant excess of genes with pivotal roles in bone homeostasis, such as PPPT3B and the height associated SUPT3H-RUNX2. We hypothesize that skeletal remodeling could be a key biological process underlying the Pygmy phenotype. In addition, we showed that these regions have most likely evolved under positive selection. These results constitute the first genetic hint of adaptive evolution in the African Pygmy phenotype, which is consistent with the independent emergence of the Pygmy height in other continents with similar environments.     

Four loci explain 83% of size variation in the horse

Horse body size varies greatly due to intense selection within each breed. American Miniatures are less than one meter tall at the withers while Shires and Percherons can exceed two meters. The genetic basis for this variation is not known. We hypothesize that the breed population structure of the horse should simplify efforts to identify genes controlling size. In support of this, here we show with genome-wide association scans (GWAS) that genetic variation at just four loci can explain the great majority of horse size variation. Unlike humans, which are naturally reproducing and possess many genetic variants with weak effects on size, we show that horses, like other domestic mammals, carry just a small number of size loci with alleles of large effect. Furthermore, three of our horse size loci contain the LCORL, HMGA2 and ZFAT genes that have previously been found to control human height. The LCORL/NCAPG locus is also implicated in cattle growth and HMGA2 is associated with dog size. Extreme size diversification is a hallmark of domestication. Our results in the horse, complemented by the prior work in cattle and dog, serve to pinpoint those very few genes that have played major roles in the rapid evolution of size during domestication.     

Identification of a genetic network for an ecologically relevant behavioural phenotype in Drosophila melanogaster

Pupation site choice of Drosophila third‐instar larvae is critical for the survival of individuals, as pupae are exposed to various biotic and abiotic dangers while immobilized during the 3–4 days of metamorphosis. This singular behavioural choice is sensitive to both environmental and genetic factors. Here, we developed a high‐throughput phenotyping approach to assay the variation in pupation height in Drosophila melanogaster, while controlling for possibly confounding factors. We find substantial variation of mean pupation height among sampled natural stocks and we show that the Drosophila Genetic Reference Panel (DGRP) reflects this variation. Using the DGRP stocks for genome‐wide association (GWA) mapping, 16 loci involved in determining pupation height could be resolved. The candidate genes in these loci are enriched for high expression in the larval central nervous system. A genetic network could be constructed from the candidate loci, which places scribble (scrib) at the centre, plus other genes known to be involved in nervous system development, such as Epidermal growth factor receptor (Egfr) and p53. Using gene disruption lines, we could functionally validate several of the initially identified loci, as well as additional loci predicted from network analysis. Our study shows that the combination of high‐throughput phenotyping with a genetic analysis of variation captured from the wild can be used to approach the genetic dissection of an environmentally relevant behavioural phenotype.     

Identification of IGF1, SLC4A4, WWOX, and SFMBT1 as hypertension susceptibility genes in Han Chinese with a genome-wide gene-based association study

Hypertension is a complex disorder with high prevalence rates all over the world. We conducted the first genome-wide gene-based association scan for hypertension in a Han Chinese population. By analyzing genome-wide single-nucleotide-polymorphism data of 400 matched pairs of young-onset hypertensive patients and normotensive controls genotyped with the Illumina HumanHap550-Duo BeadChip, 100 susceptibility genes for hypertension were identified and also validated with permutation tests. Seventeen of the 100 genes exhibited differential allelic and expression distributions between patient and control groups. These genes provided a good molecular signature for classifying hypertensive patients and normotensive controls. Among the 17 genes, IGF1, SLC4A4, WWOX, and SFMBT1 were not only identified by our gene-based association scan and gene expression analysis but were also replicated by a gene-based association analysis of the Hong Kong Hypertension Study. Moreover, cis-acting expression quantitative trait loci associated with the differentially expressed genes were found and linked to hypertension. IGF1, which encodes insulin-like growth factor 1, is associated with cardiovascular disorders, metabolic syndrome, decreased body weight/size, and changes of insulin levels in mice. SLC4A4, which encodes the electrogenic sodium bicarbonate cotransporter 1, is associated with decreased body weight/size and abnormal ion homeostasis in mice. WWOX, which encodes the WW domain-containing protein, is related to hypoglycemia and hyperphosphatemia. SFMBT1, which encodes the scm-like with four MBT domains protein 1, is a novel hypertension gene. GRB14, TMEM56 and KIAA1797 exhibited highly significant differential allelic and expressed distributions between hypertensive patients and normotensive controls. GRB14 was also found relevant to blood pressure in a previous genetic association study in East Asian populations. TMEM56 and KIAA1797 may be specific to Taiwanese populations, because they were not validated by the two replication studies. Identification of these genes enriches the collection of hypertension susceptibility genes, thereby shedding light on the etiology of hypertension in Han Chinese populations.     

Juvenile hormone regulates vitellogenin gene expression through insulin-like peptide signaling pathway in the red flour beetle, Tribolium castaneum

Our recent studies identified juvenile hormone (JH) and nutrition as the two key signals that regulate vitellogenin (Vg) gene expression in the red flour beetle, Tribolium castaneum. Juvenile hormone regulation of Vg synthesis has been known for a long time in several insects, but the mechanism of JH action is not known. Experiments were conducted to determine the mechanism of action of these two signals in regulation of Vg gene expression. Injection of bovine insulin or FOXO double-stranded RNA into the previtellogenic, starved, or JH-deficient female adults increased Vg mRNA and protein levels, thereby implicating the pivotal role for insulin-like peptide signaling in the regulation of Vg gene expression and possible cross-talk between JH and insulin-like peptide signaling pathways. Reduction in JH synthesis or its action by RNAi-mediated silencing of genes coding for acid methyltransferase or methoprene-tolerant decreased expression of genes coding for insulin-like peptides (ILPs) and influenced FOXO subcellular localization, resulting in the down-regulation of Vg gene expression. Furthermore, JH application to previtellogenic female beetles induced the expression of genes coding for ILP2 and ILP3, and induced Vg gene expression. FOXO protein expressed in baculovirus system binds to FOXO response element present in the Vg gene promoter. These data suggest that JH functions through insulin-like peptide signaling pathway to regulate Vg gene expression.     

Beyond pathways: genetic dissection of tocopherol content in maize kernels by combining linkage and association analyses

Although tocopherols play an important role in plants and animals, the genetic architecture of tocopherol content in maize kernels has remained largely unknown. In this study, linkage and association analyses were conducted to examine the genetic architecture of tocopherol content in maize kernels. Forty‐one unique quantitative trait loci (QTLs) were identified by linkage mapping in six populations of recombinant inbred lines (RILs). In addition, 32 significant loci were detected via genome‐wide association study (GWAS), 18 of which colocalized with the QTLs identified by linkage mapping. Fine mapping of a major QTL validated the accuracy of GWAS and QTL mapping results and suggested a role for nontocopherol pathway genes in the modulation of natural tocopherol variation. We provided genome‐wide evidence that genes involved in fatty acid metabolism, chlorophyll metabolism and chloroplast function may affect natural variation in tocopherols. These findings were confirmed through mutant analysis of a particular gene from the fatty acid pathway. In addition, the favourable alleles for many of the significant SNPs/QTLs represented rare alleles in natural populations. Together, our results revealed many novel genes that are potentially involved in the variation of tocopherol content in maize kernels. Pyramiding of the favourable alleles of the newly elucidated genes and the well‐known tocopherol pathway genes would greatly improve tocopherol content in maize.     

Indirect evidence for the genetic determination of short stature in African Pygmies

Central African Pygmy populations are known to be the shortest human populations worldwide. Many evolutionary hypotheses have been proposed to explain this short stature: adaptation to food limitations, climate, forest density, or high mortality rates. However, such hypotheses are difficult to test given the lack of long-term surveys and demographic data. Whether the short stature observed nowadays in African Pygmy populations as compared to their Non-Pygmy neighbors is determined by genetic factors remains widely unknown. Here, we study a uniquely large new anthropometrical dataset comprising more than 1,000 individuals from 10 Central African Pygmy and neighboring Non-Pygmy populations, categorized as such based on cultural criteria rather than height. We show that climate, or forest density may not play a major role in the difference in adult stature between existing Pygmies and Non-Pygmies, without ruling out the hypothesis that such factors played an important evolutionary role in the past. Furthermore, we analyzed the relationship between stature and neutral genetic variation in a subset of 213 individuals and found that the Pygmy individuals' stature was significantly positively correlated with levels of genetic similarity with the Non-Pygmy gene-pool for both men and women. Overall, we show that a Pygmy individual exhibiting a high level of genetic admixture with the neighboring Non-Pygmies is likely to be taller. These results show for the first time that the major morphological difference in stature found between Central African Pygmy and Non-Pygmy populations is likely determined by genetic factors.     

Candidate genes for colour and vision exhibit signals of selection across the pied flycatcher (Ficedula hypoleuca) breeding range

The role of natural selection in shaping adaptive trait differentiation in natural populations has long been recognized. Determining its molecular basis, however, remains a challenge. Here, we search for signals of selection in candidate genes for colour and its perception in a passerine bird. Pied flycatcher plumage varies geographically in both its structural and pigment-based properties. Both characteristics appear to be shaped by selection. A single-locus outlier test revealed 2 of 14 loci to show significantly elevated signals of divergence. The first of these, the follistatin gene, is expressed in the developing feather bud and is found in pathways with genes that determine the structure of feathers and may thus be important in generating variation in structural colouration. The second is a gene potentially underlying the ability to detect this variation: SWS1 opsin. These two loci were most differentiated in two Spanish pied flycatcher populations, which are also among the populations that have the highest UV reflectance. The follistatin and SWS1 opsin genes thus provide strong candidates for future investigations on the molecular basis of adaptively significant traits and their co-evolution.     

Serotonin receptors antagonistically modulate Caenorhabditis elegans longevity

The neurotransmitter serotonin has been implicated in affecting the variation of longevity in natural Drosophila populations and age-related diseases in mammals. Based on these observations, it has been predicted that serotonin signal, perhaps at levels of serotonin biosynthesis, may control lifespan. Here, we investigated a variety of mutations in serotonin-signal genes, including serotonin biosynthesis genes, a serotonin transporter gene, and serotonin receptor genes. Despite this prediction, mutations in the serotonin biosynthesis genes had little or modest effects on lifespan, while the mod-5 mutation with increased availability of serotonin caused a modest life-shortening effect. In contrast, a deletion mutation of the ser-1 serotonin receptor gene increased longevity by up to 46%, likely through the insulin/insulin-like growth factor 1 pathway. This result suggests an interaction between the serotonin pathway and the insulin/insulin-like growth factor 1 pathway. A deletion mutation of another serotonin receptor gene, ser-4 , shortened early to mid lifespan. The results suggest that serotonin signal antagonistically modulates longevity through different serotonin receptors. This study may indicate serotonin receptors as a potential target for antigeric interventions.     

Latitudinal clines in<Emphasis Type="Italic">Drosophila melanogaster</Emphasis>: Body size,allozyme frequencies,inversion frequencies,and the insulin-signalling pathway

Many latitudinal clines exist inDrosophila melanogaster: in adult body size, in allele frequency at allozyme loci, and in frequencies of common cosmopolitan inversions. The question is raised whether these latitudinal clines are causally related. This review aims to connect data from two very different fields of study, evolutionary biology and cell biology, in explaining such natural genetic variation inD. melanogaster body size and development time. It is argued that adult body size clines, inversion frequency clines, and clines in allele frequency at loci involved in glycolysis and glycogen storage are part of the same adaptive strategy. Selection pressure is expected to differ at opposite ends of the clines. At high latitudes, selection onD. melanogaster would favour high larval growth rate at low temperatures, and resource storage in adults to survive winter. At low latitudes selection would favour lower larval critical size to survive crowding, and increased male activity leading to high male reproductive success. Studies of the insulin-signalling pathway inD. melanogaster point to the involvement of this pathway in metabolism and adult body size. The genes involved in the insulin-signalling pathway are associated with common cosmopolitan inversions that show latitudinal clines. Each chromosome region connected with a large common cosmopolitan inversion possesses a gene of the insulin transmembrane complex, a gene of the intermediate pathway and a gene of the TOR branch. The hypothesis is presented that temperateD. melanogaster populations have a higher frequency of a ’thrifty’ genotype corresponding to high insulin level or high signal level, while tropical populations possess a more ’spendthrift’ genotype corresponding to low insulin or low signal level.     

Genome-wide polymorphisms show unexpected targets of natural selection

Natural selection can act on all the expressed genes of an individual, leaving signatures of genetic differentiation or diversity at many loci across the genome. New power to assay these genome-wide effects of selection comes from associating multi-locus patterns of polymorphism with gene expression and function. Here, we performed one of the first genome-wide surveys in a marine species, comparing purple sea urchins, Strongylocentrotus purpuratus, from two distant locations along the species' wide latitudinal range. We examined 9112 polymorphic loci from upstream non-coding and coding regions of genes for signatures of selection with respect to gene function and tissue- and ontogenetic gene expression. We found that genetic differentiation (F(ST)) varied significantly across functional gene classes. The strongest enrichment occurred in the upstream regions of E3 ligase genes, enzymes known to regulate protein abundance during development and environmental stress. We found enrichment for high heterozygosity in genes directly involved in immune response, particularly NALP genes, which mediate pro-inflammatory signals during bacterial infection. We also found higher heterozygosity in immune genes in the southern population, where disease incidence and pathogen diversity are greater. Similar to the major histocompatibility complex in mammals, balancing selection may enhance genetic diversity in the innate immune system genes of this invertebrate. Overall, our results show that how genome-wide polymorphism data coupled with growing databases on gene function and expression can combine to detect otherwise hidden signals of selection in natural populations.     

The functions of insulin signaling: size isn't everything,even in Drosophila

Mammalian insulin and insulin-like growth factors (IGFs) signal through several receptors with different ligand specificities to regulate metabolism and growth. This regulation is defective in diabetes and in a wide variety of human tumors. Recent analysis in Drosophila melanogaster has revealed that insulin-like molecules (known as DILPs in flies) also control growth and metabolism, but probably do so by signaling through a single insulin receptor (InR). The intracellular signaling molecules regulated by this receptor are highly evolutionarily conserved. Work in flies has helped to dissect the network of InR-regulated intracellular signaling pathways and identify some of the critical players in these pathways and in interacting signaling cascades. Surprisingly, these studies have shown that DILPs control tissue and body growth primarily by regulating cell growth and cell size. Changes in cell growth produced by these molecules may subsequently modulate the rate of cell proliferation in a cell type-specific fashion. At least part of this growth effect is mediated by two small groups of neurons in the Drosophila brain, which secrete DILPs into the circulatory system at levels that are modulated by nutrition. This signaling center is also involved in DILP-dependent control of the fly's rate of development, fertility, and life span. These surprisingly diverse functions of InR signaling, which appear to be conserved in all higher animals, reflect a central role for this pathway in coordinating development, physiology, and properly proportioned growth of the organism in response to its nutritional state. Studies in flies are providing important new insights into the biology of this system, and the identification of novel components in the InR-regulated signaling cascade is already beginning to inform the development of new therapeutic strategies for insulin-linked diseases in the clinic.     

Genomics,environment and balancing selection in behaviourally bimodal populations: The caribou case

Selection forces that favour different phenotypes in different environments can change frequencies of genes between populations along environmental clines. Clines are also compatible with balancing forces, such as negative frequency‐dependent selection (NFDS), which maintains phenotypic polymorphisms within populations. For example, NFDS is hypothesized to maintain partial migration, a dimorphic behavioural trait prominent in species where only a fraction of the population seasonally migrates. Overall, NFDS is believed to be a common phenomenon in nature, yet a scarcity of studies were published linking naturally occurring allelic variation with bimodal or multimodal phenotypes and balancing selection. We applied a Pool‐seq approach and detected selection on alleles associated with environmental variables along a North–South gradient in western North American caribou, a species displaying partially migratory behaviour. On 51 loci, we found a signature of balancing selection, which could be related to NFDS and ultimately the maintenance of the phenotypic polymorphisms known within these populations. Yet, remarkably, we detected directional selection on a locus when our sample was divided into two behaviourally distinctive groups regardless of geographic provenance (a subset of GPS‐collared migratory or sedentary individuals), indicating that, within populations, phenotypically homogeneous groups were genetically distinctive. Loci under selection were linked to functional genes involved in oxidative stress response, body development and taste perception. Overall, results indicated genetic differentiation along an environmental gradient of caribou populations, which we found characterized by genes potentially undergoing balancing selection. We suggest that the underlining balancing force, NFDS, plays a strong role within populations harbouring multiple haplotypes and phenotypes, as it is the norm in animals, plants and humans too.     

A genome-wide scan for signatures of directional selection in domesticated pigs

Background

Animal domestication involved drastic phenotypic changes driven by strong artificial selection and also resulted in new populations of breeds, established by humans. This study aims to identify genes that show evidence of recent artificial selection during pig domestication.

Results

Whole-genome resequencing of 30 individual pigs from domesticated breeds, Landrace and Yorkshire, and 10 Asian wild boars at ~16-fold coverage was performed resulting in over 4.3 million SNPs for 19,990 genes. We constructed a comprehensive genome map of directional selection by detecting selective sweeps using an F_ST-based approach that detects directional selection in lineages leading to the domesticated breeds and using a haplotype-based test that detects ongoing selective sweeps within the breeds. We show that candidate genes under selection are significantly enriched for loci implicated in quantitative traits important to pig reproduction and production. The candidate gene with the strongest signals of directional selection belongs to group III of the metabolomics glutamate receptors, known to affect brain functions associated with eating behavior, suggesting that loci under strong selection include loci involved in behaviorial traits in domesticated pigs including tameness.

Conclusions

We show that a significant proportion of selection signatures coincide with loci that were previously inferred to affect phenotypic variation in pigs. We further identify functional enrichment related to behavior, such as signal transduction and neuronal activities, for those targets of selection during domestication in pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1330-x) contains supplementary material, which is available to authorized users.     

A scan of molecular variation leads to the narrow localization of a selective sweep affecting both Afrotropical and cosmopolitan populations of Drosophila melanogaster

Drosophila melanogaster originated in tropical Africa but has achieved a cosmopolitan distribution in association with human habitation. Cosmopolitan populations of D. melanogaster are known to have reduced genetic variation, particularly on the X chromosome. However, the relative importance of population bottlenecks and selective sweeps in explaining this reduction is uncertain. We surveyed variation at 31 microsatellites across a 330-kb section of the X chromosome located between the white and kirre genes. Two linked clusters of loci were observed with reduced variation and a skew toward rare alleles in both an Ecuador and a Zimbabwe population sample. Examining Zimbabwe DNA sequence polymorphism within one of these regions allowed us to localize a selective sweep to a 361-bp window within the 5' regulatory region of the roughest gene, with one nucleotide substitution representing the best candidate for the target of selection. Estimates of sweep age suggested that this fixation event occurred prior to the expansion of D. melanogaster from sub-Saharan Africa. For both putative sweep regions in our data set, cosmopolitan populations showed wider footprints of selection compared to those in Zimbabwe. This pattern appears consistent with the demographic amplification of preexisting sweep signals due to one or more population bottlenecks.     

Phosphorylation of receptors for insulin and insulin-like growth factor I. Effects of hormones and phorbol esters

The phosphorylation of receptors for insulin and insulin-like growth factor I was studied by phosphoamino acid analysis and tryptic phosphopeptide maps in an attempt to determine if protein kinase C is involved in their phosphorylation in response to insulin and insulin-like growth factor I, respectively. Two cell lines were utilized, Hep G2 and IM-9 cells. sn-1,2-Dioctanoylglycerol and 12-O-tetradecanoylphorbol 13-acetate (TPA), agents known to activate protein kinase C, stimulated the phosphorylation of the beta subunits of both receptors, as did their hormones. In unstimulated cells, phosphorylation of the insulin receptor occurred on seryl and to a lesser extent on threonyl residues. TPA stimulated seryl and threonyl phosphorylation that resulted in the appearance of four major phosphoserine-containing phosphopeptides which were not detected in the basal state and an increase in phosphorylation of a phosphothreonine-containing peptide which was present in the basal state. Insulin treatment resulted in the appearance of three major phosphotyrosine-containing tryptic peptides. In IM-9 cells, insulin also increased the phosphoserine and possibly the phosphothreonine content of the beta subunit. In both cells, the major phosphoserine-containing peptides that were stimulated by TPA were not detected following treatment with insulin. Very similar results, including similar peptide maps, were obtained for the insulin-like growth factor I receptor from cells treated with TPA and insulin-like growth factor I. Although not entirely conclusive, these results suggest that the insulin- and insulin-like growth factor I-stimulated phosphorylation of their receptors does not result from activation of protein kinase C.     

Genetic Variation of Microsatellite Loci in the Major Histocompatibility Complex (MHC) Region in the Southeast Asian House Mouse (<Emphasis Type="Italic">Mus musculus castaneus</Emphasis>)

Major histocompatibility complex (MHC) genes are the most polymorphic loci known for vertebrates. Here we employed five microsatellite loci closely linked to the MHC region in an attempt to study the amount of genetic variation in 19 populations of the southeast Asian house mouse (Mus musculus castaneus) in Taiwan. The overall polymorphism at the five loci was high (H_e = 0.713), and the level of polymorphism varied from locus to locus. Furthermore, in order to investigate if selection is operating on MHC genes in natural mouse populations, we compared the extent and pattern of genetic variation for the MHC-linked microsatellite loci (the MHC loci) with those for the microsatellite loci located outside the MHC region (the non-MHC loci). The number of alleles and the logarithm of variance in repeat number were significantly higher for the MHC loci than for the non-MHC loci, presumably reflecting linkage to a locus under balancing selection. Although three statistical tests used do not provide support for selection, their lack of support may be due to low statistical power of the tests, to weakness of selection, or to a profound effect of genetic drift reducing the signature of balancing selection. Our results also suggested that the populations in the central and the southwestern regions of Taiwan might be one part of a metapopulation structure.