Evolutionary Genetic Signatures of Selection on Bone-Related Variation within Human and Chimpanzee Populations.

Bone strength and the incidence and severity of skeletal disorders vary significantly among human populations, due in part to underlying genetic differentiation. While clinical models predict that this variation is largely deleterious, natural population variation unrelated to disease can go unnoticed, altering our perception of how natural selection has shaped bone morphologies over deep and recent time periods. Here, we conduct the first comparative population-based genetic analysis of the main bone structural protein gene, collagen type I α 1 (COL1A1), in clinical and 1000 Genomes Project datasets in humans, and in natural populations of chimpanzees. Contrary to predictions from clinical studies, we reveal abundant COL1A1 amino acid variation, predicted to have little association with disease in the natural population. We also find signatures of positive selection associated with intron haplotype structure, linkage disequilibrium, and population differentiation in regions of known gene expression regulation in humans and chimpanzees. These results recall how recent and deep evolutionary regimes can be linked, in that bone morphology differences that developed among vertebrates over 450 million years of evolution are the result of positive selection on subtle type I collagen functional variation segregating within populations over time.

Positive selection on human gamete-recognition genes.

Coevolution of genes that encode interacting proteins expressed on the surfaces of sperm and eggs can lead to variation in reproductive compatibility between mates and reproductive isolation between members of different species. Previous studies in mice and other mammals have focused in particular on evidence for positive or diversifying selection that shapes the evolution of genes that encode sperm-binding proteins expressed in the egg coat or zona pellucida (ZP). By fitting phylogenetic models of codon evolution to data from the 1000 Genomes Project, we identified candidate sites evolving under diversifying selection in the human genes ZP3 and ZP2. We also identified one candidate site under positive selection in C4BPA, which encodes a repetitive protein similar to the mouse protein ZP3R that is expressed in the sperm head and binds to the ZP at fertilization. Results from several additional analyses that applied population genetic models to the same data were consistent with the hypothesis of selection on those candidate sites leading to coevolution of sperm- and egg-expressed genes. By contrast, we found no candidate sites under selection in a fourth gene (ZP1) that encodes an egg coat structural protein not directly involved in sperm binding. Finally, we found that two of the candidate sites (in C4BPA and ZP2) were correlated with variation in family size and birth rate among Hutterite couples, and those two candidate sites were also in linkage disequilibrium in the same Hutterite study population. All of these lines of evidence are consistent with predictions from a previously proposed hypothesis of balancing selection on epistatic interactions between C4BPA and ZP3 at fertilization that lead to the evolution of co-adapted allele pairs. Such patterns also suggest specific molecular traits that may be associated with both natural reproductive variation and clinical infertility.

Predicting childhood effortful control from interactions between early parenting quality and children's dopamine transporter gene haplotypes.

Children's observed effortful control (EC) at 30, 42, and 54 months (n = 145) was predicted from the interaction between mothers' observed parenting with their 30-month-olds and three variants of the solute carrier family C6, member 3 (SLC6A3) dopamine transporter gene (single nucleotide polymorphisms in intron8 and intron13, and a 40 base pair variable number tandem repeat [VNTR] in the 3'-untranslated region [UTR]), as well as haplotypes of these variants. Significant moderating effects were found. Children without the intron8-A/intron13-G, intron8-A/3'-UTR VNTR-10, or intron13-G/3'-UTR VNTR-10 haplotypes (i.e., haplotypes associated with the reduced SLC6A3 gene expression and thus lower dopamine functioning) appeared to demonstrate altered levels of EC as a function of maternal parenting quality, whereas children with these haplotypes demonstrated a similar EC level regardless of the parenting quality. Children with these haplotypes demonstrated a trade-off, such that they showed higher EC, relative to their counterparts without these haplotypes, when exposed to less supportive maternal parenting. The findings revealed a diathesis-stress pattern and suggested that different SLC6A3 haplotypes, but not single variants, might represent different levels of young children's sensitivity/responsivity to early parenting.

Interactions among catechol-O-methyltransferase genotype, parenting, and sex predict children's internalizing symptoms and inhibitory control: Evidence for differential susceptibility.

We used sex, observed parenting quality at 18 months, and three variants of the catechol-O-methyltransferase gene (Val158Met [rs4680], intron1 [rs737865], and 3'-untranslated region [rs165599]) to predict mothers' reports of inhibitory and attentional control (assessed at 42, 54, 72, and 84 months) and internalizing symptoms (assessed at 24, 30, 42, 48, and 54 months) in a sample of 146 children (79 male). Although the pattern for all three variants was very similar, Val158Met explained more variance in both outcomes than did intron1, the 3'-untranslated region, or a haplotype that combined all three catechol-O-methyltransferase variants. In separate models, there were significant three-way interactions among each of the variants, parenting, and sex, predicting the intercepts of inhibitory control and internalizing symptoms. Results suggested that Val158Met indexes plasticity, although this effect was moderated by sex. Parenting was positively associated with inhibitory control for methionine-methionine boys and for valine-valine/valine-methionine girls, and was negatively associated with internalizing symptoms for methionine-methionine boys. Using the "regions of significance" technique, genetic differences in inhibitory control were found for children exposed to high-quality parenting, whereas genetic differences in internalizing were found for children exposed to low-quality parenting. These findings provide evidence in support of testing for differential susceptibility across multiple outcomes.

Development of Ego-Resiliency: Relations to Observed Parenting and Polymorphisms in the Serotonin Transporter Gene During Early Childhood.

We used observed parenting behaviors, along with genetic variants and haplotypes of the serotonin transporter gene (SLC6A4), as predictors of children's ego-resiliency during early childhood (N =153). Quality of mothers' parenting was observed at 18 months of age and mothers' reports of ego-resiliency were collected at six time points from 18 to 84 months. Genetic data were collected at 72 months. Observed parenting was positively associated with initial levels of children's ego-resiliency. Furthermore, although individual genetic variants of the serotonin transporter gene (LPR, STin2) were not associated with ego-resiliency, the S10 haplotype (that combines information from these two variants) was negatively associated with initial levels of ego-resiliency. Both parenting and serotonin genetic variation uniquely predicted children's ego-resiliency, suggesting an additive effect of genetic and parental factors.

Interactions between serotonin transporter gene haplotypes and quality of mothers' parenting predict the development of children's noncompliance.

The LPR and STin2 polymorphisms of the serotonin transporter gene (SLC6A4) were combined into haplotypes that, together with quality of maternal parenting, were used to predict initial levels and linear change in children's (N = 138) noncompliance and aggression from age 18-54 months. Quality of mothers' parenting behavior was observed when children were 18 months old, and nonparental caregivers' reports of noncompliance and aggression were collected annually from 18 to 54 months of age. Quality of early parenting was negatively related to the slope of noncompliance only for children with the LPR-S/STin2-10 haplotype and to 18-month noncompliance only for children with haplotypes that did not include LPR-S. The findings support the notion that SLC6A4 haplotypes index differential susceptibility to variability in parenting quality, with certain haplotypes showing greater reactivity to both supportive and unsupportive environments. These different genetic backgrounds likely reflect an evolutionary response to variation in the parenting environment.

Comparative vertebrate evolutionary analyses of type I collagen: potential of COL1a1 gene structure and intron variation for common bone-related diseases.

Collagen type I alpha 1 (COL1a1), which encodes the primary subunit of type I collagen, the main structural and most abundant protein in vertebrates, harbors hundreds of mutations linked to human diseases like osteoporosis and osteogenesis imperfecta. Previous studies have attempted to predict the phenotypic severity associated with type I collagen mutations, yet an evolutionary analysis that compares historical and recent selective pressures, including across noncoding regions, has never been conducted. Here, we use a comparative genomic and species evolutionary analysis representing ∼450 My of vertebrate history to investigate functional constraints associated with both exons and introns of the >17-kb COL1a1 gene. We find that although the COL1a1 amino acid sequence is highly conserved, there are both spatial and temporal signatures of varying selective constraint across protein domains. Furthermore, sites of high evolutionary constraint significantly correlate with the location of disease-associated mutations, the latter of which also cluster with respect to specific severity classes typically categorized in clinical studies. Finally, we find that COL1a1 introns are significantly short in length with high GC content, patterns that are shared across highly diverged vertebrates, and which may be a signature of strong stabilizing selection for high COL1a1 gene expression. In conclusion, although previous studies focused on COL1a1 coding regions, the current results implicate introns as areas of high selective constraint and targets of bone-related phenotypic variation. From a broader perspective, our comparative evolutionary approach provides further resolution to models predicting mutations associated with bone-related function and disease severity.

The beta -globin recombinational hotspot reduces the effects of strong selection around HbC, a recently arisen mutation providing resistance to malaria.

Recombination is expected to reduce the effect of selection on the extent of linkage disequilibrium (LD), but the impact that recombinational hotspots have on sites linked to selected mutations has not been investigated. We empirically determine chromosomal linkage phase for 5.2 kb spanning the beta -globin gene and hotspot. We estimate that the HbC mutation, which is positively selected because of malaria, originated <5,000 years ago and that selection coefficients are 0.04-0.09. Despite strong selection and the recent origin of the HbC allele, recombination (crossing-over or gene conversion) is observed within 1 kb 5' of the selected site on more than one-third of the HbC chromosomes sampled. The rapid decay in LD upstream of the HbC allele demonstrates the large effect the ss-globin hotspot has in mitigating the effects of positive selection on linked variation.

Contrasting patterns of Y chromosome and mtDNA variation in Africa: evidence for sex-biased demographic processes.

To investigate associations between genetic, linguistic, and geographic variation in Africa, we type 50 Y chromosome SNPs in 1122 individuals from 40 populations representing African geographic and linguistic diversity. We compare these patterns of variation with those that emerge from a similar analysis of published mtDNA HVS1 sequences from 1918 individuals from 39 African populations. For the Y chromosome, Mantel tests reveal a strong partial correlation between genetic and linguistic distances (r=0.33, P=0.001) and no correlation between genetic and geographic distances (r=-0.08, P>0.10). In contrast, mtDNA variation is weakly correlated with both language (r=0.16, P=0.046) and geography (r=0.17, P=0.035). AMOVA indicates that the amount of paternal among-group variation is much higher when populations are grouped by linguistics (Phi(CT)=0.21) than by geography (Phi(CT)=0.06). Levels of maternal genetic among-group variation are low for both linguistics and geography (Phi(CT)=0.03 and 0.04, respectively). When Bantu speakers are removed from these analyses, the correlation with linguistic variation disappears for the Y chromosome and strengthens for mtDNA. These data suggest that patterns of differentiation and gene flow in Africa have differed for men and women in the recent evolutionary past. We infer that sex-biased rates of admixture and/or language borrowing between expanding Bantu farmers and local hunter-gatherers played an important role in influencing patterns of genetic variation during the spread of African agriculture in the last 4000 years.