Arslan A Zaidi

Modeling 3D Facial Shape from DNA

Human facial diversity is substantial, complex, and largely scientifically unexplained. We used spatially dense quasi-landmarks to measure face shape in population samples with mixed West African and European ancestry from three locations (United States, Brazil, and Cape Verde). Using bootstrapped response-based imputation modeling (BRIM), we uncover the relationships between facial variation and the effects of sex, genomic ancestry, and a subset of craniofacial candidate genes. The facial effects of these variables are summarized as response-based imputed predictor (RIP) variables, which are validated using self-reported sex, genomic ancestry, and observer-based facial ratings (femininity and proportional ancestry) and judgments (sex and population group). By jointly modeling sex, genomic ancestry, and genotype, the independent effects of particular alleles on facial features can be uncovered. Results on a set of 20 genes showing significant effects on facial features provide support for this approach as a novel means to identify genes affecting normal-range facial features and for approximating the appearance of a face from genetic markers.

The conditioning of intervention effects on early adolescent alcohol use by maternal involvement and dopamine receptor D4 (DRD4) and serotonin transporter linked polymorphic region (5-HTTLPR) genetic variants.

Data drawn from the in-home subsample of the PROSPER intervention dissemination trial were used to investigate the moderation of intervention effects on underage alcohol use by maternal involvement and candidate genes. The primary gene examined was dopamine receptor D4 (DRD4). Variation in this gene and maternal involvement were hypothesized to moderate the influence of intervention status on alcohol use. The PROSPER data used were drawn from 28 communities randomly assigned to intervention or comparison conditions. Participating youth were assessed in five in-home interviews from sixth to ninth grades. A main effect of sixth-grade pretest maternal involvement on ninth-grade alcohol use was found. Neither intervention status nor DRD4 variation was unconditionally linked to ninth-grade drinking. However, moderation analyses revealed a significant three-way interaction among DRD4 status, maternal involvement, and intervention condition. Follow-up analyses revealed that prevention reduced drinking risk, but only for youth with at least one DRD4 seven-repeat allele who reported average or greater pretest levels of maternal involvement. To determine if this conditional pattern was limited to the DRD4 gene, we repeated analyses using the serotonin transporter linked polymorphic region site near the serotonin transporter gene. The results for this supplemental analysis revealed a significant three-way interaction similar but not identical to that found for DRD4.

Genome-wide mapping of global-to-local genetic effects on human facial shape.

Genome-wide association scans of complex multipartite traits like the human face typically use preselected phenotypic measures. Here we report a data-driven approach to phenotyping facial shape at multiple levels of organization, allowing for an open-ended description of facial variation while preserving statistical power. In a sample of 2,329 persons of European ancestry, we identified 38 loci, 15 of which replicated in an independent European sample (n = 1,719). Four loci were completely new. For the others, additional support (n = 9) or pleiotropic effects (n = 2) were found in the literature, but the results reported here were further refined. All 15 replicated loci highlighted distinctive patterns of global-to-local genetic effects on facial shape and showed enrichment for active chromatin elements in human cranial neural crest cells, suggesting an early developmental origin of the facial variation captured. These results have implications for studies of facial genetics and other complex morphological traits.The authors report a data-driven approach to phenotyping 3D facial shape. They apply their methodology to 2,329 individuals of European ancestry and identify 38 loci that associate with specific facial morphologies, some of which overlap with neural-crest-specific regulatory regions.

Investigating the case of human nose shape and climate adaptation

The evolutionary reasons for variation in nose shape across human populations have been subject to continuing debate. An import function of the nose and nasal cavity is to condition inspired air before it reaches the lower respiratory tract. For this reason, it is thought the observed differences in nose shape among populations are not simply the result of genetic drift, but may be adaptations to climate. To address the question of whether local adaptation to climate is responsible for nose shape divergence across populations, we use Qst–Fst comparisons to show that nares width and alar base width are more differentiated across populations than expected under genetic drift alone. To test whether this differentiation is due to climate adaptation, we compared the spatial distribution of these variables with the global distribution of temperature, absolute humidity, and relative humidity. We find that width of the nares is correlated with temperature and absolute humidity, but not with relative humidity. We conclude that some aspects of nose shape may indeed have been driven by local adaptation to climate. However, we think that this is a simplified explanation of a very complex evolutionary history, which possibly also involved other non-neutral forces such as sexual selection.

High Levels of Copy Number Variation of Ampliconic Genes across Major Human Y Haplogroups

Abstract Because of its highly repetitive nature, the human male-specific Y chromosome remains understudied. It is important to investigate variation on the Y chromosome to understand its evolution and contribution to phenotypic variation, including infertility. Approximately 20% of the human Y chromosome consists of ampliconic regions which include nine multi-copy gene families. These gene families are expressed exclusively in testes and usually implicated in spermatogenesis. Here, to gain a better understanding of the role of the Y chromosome in human evolution and in determining sexually dimorphic traits, we studied ampliconic gene copy number variation in 100 males representing ten major Y haplogroups world-wide. Copy number was estimated with droplet digital PCR. In contrast to low nucleotide diversity observed on the Y in previous studies, here we show that ampliconic gene copy number diversity is very high. A total of 98 copy-number-based haplotypes were observed among 100 individuals, and haplotypes were sometimes shared by males from very different haplogroups, suggesting homoplasies. The resulting haplotypes did not cluster according to major Y haplogroups. Overall, only two gene families (RBMY and TSPY) showed significant differences in copy number among major Y haplogroups, and the haplogroup of a male could not be predicted based on his ampliconic gene copy numbers. Finally, we did not find significant correlations either between copy number variation and individual’s height, or between the former and facial masculinity/femininity. Our results suggest rapid evolution of ampliconic gene copy numbers on the human Y, and we discuss its causes.

A Population Phylogenetic View of Mitochondrial Heteroplasmy

The mitochondria contained within the human body are genetically diverse. This type of variation, called heteroplasmy, is emerging as an important factor... The mitochondrion has recently emerged as an active player in myriad cellular processes. Additionally, it was recently shown that >200 diseases are known to be linked to variants in mitochondrial DNA or in nuclear genes interacting with mitochondria. This has reinvigorated interest in its biology and population genetics. Mitochondrial heteroplasmy, or genotypic variation of mitochondria within an individual, is now understood to be common in humans and important in human health. However, it is still not possible to make quantitative predictions about the inheritance of heteroplasmy and its proliferation within the body, partly due to the lack of an appropriate model. Here, we present a population-genetic framework for modeling mitochondrial heteroplasmy as a process that occurs on an ontogenetic phylogeny, with genetic drift and mutation changing heteroplasmy frequencies during the various developmental processes represented in the phylogeny. Using this framework, we develop a Bayesian inference method for inferring rates of mitochondrial genetic drift and mutation at different stages of human life. Applying the method to previously published heteroplasmy frequency data, we demonstrate a severe effective germline bottleneck comprised of the cumulative genetic drift occurring between the divergence of germline and somatic cells in the mother, and the separation of germ layers in the offspring. Additionally, we find that the two somatic tissues we analyze here undergo tissue-specific bottlenecks during embryogenesis, less severe than the effective germline bottleneck, and that these somatic tissues experience little additional genetic drift during adulthood. We conclude with a discussion of possible extensions of the ontogenetic phylogeny framework and its possible applications to other ontogenetic processes in addition to mitochondrial heteroplasmy.

Mito-nuclear effects uncovered in admixed populations

To function properly, mitochondria utilize products of 37 and >1,000 genes encoded by the mitochondrial and nuclear genomes, respectively, which should be compatible with each other. Discordance between mitochondrial and nuclear genetic ancestry could contribute to phenotypic variation in admixed populations. Here we explored potential mito-nuclear incompatibility in six admixed human populations from the Americas: African Americans, African Caribbeans, Colombians, Mexicans, Peruvians, and Puerto Ricans. For individuals in these populations, we determined nuclear genome proportions derived from Africans, Europeans, and Native Americans, the geographic origins of the mitochondrial DNA (mtDNA), as well as mtDNA copy number in lymphoblastoid cell lines. By comparing nuclear vs. mitochondrial ancestry in admixed populations, we show that, first, mtDNA copy number decreases with increasing discordance between nuclear and mitochondrial DNA ancestry, in agreement with mito-nuclear incompatibility. The direction of this effect is consistent across mtDNA haplogroups of different geographic origins. This observation suggests suboptimal regulation of mtDNA replication when its components are encoded by nuclear and mtDNA genes with different ancestry. Second, while most populations analyzed exhibit no such trend, in Puerto Ricans and African Americans we find a significant enrichment of ancestry at nuclear-encoded mitochondrial genes towards the source populations contributing the most prevalent mtDNA haplogroups (Native American and African, respectively). This likely reflects compensatory effects of selection in recovering mito-nuclear interactions optimized in the source populations. Our results provide the first evidence of mito-nuclear effects in human admixed populations and we discuss its implications for human health and disease.

Facial masculinity does not appear to be a condition-dependent male ornament in humans and does not reflect MHC heterozygosity

Facial masculinity is thought to be a condition-dependent male ornament, reflecting immunocompetence in humans. To test this hypothesis, we calculated an objective measure of facial masculinity/femininity using three-dimensional images in a large sample (N = 1,233) of people of European ancestry. We show that facial masculinity is positively correlated with adult height in both males and females. This suggests that variation in growth contributes, at least in part, to variation in facial masculinity, which is characteristic of condition-dependent traits. However, facial masculinity scales with growth similarly in males and females, suggesting that facial masculinity is not specifically a male ornament. Additionally, we measured immunocompetence via heterozygosity at the major histocompatibility complex (MHC), a well known genetic marker of immunity. We show that while height is positively correlated with MHC heterozygosity, facial masculinity is not. Thus, facial masculinity does not reflect immunocompetence measured by MHC heterozygosity in humans as thought previously. Overall, we find no support for the idea that facial masculinity is a condition-dependent male ornament that has evolved to indicate immunocompetence.

Correction: Investigating the case of human nose shape and climate adaptation

[This corrects the article DOI: 10.1371/journal.pgen.1006616.].

Copy number variation of ampliconic genes across major human Y haplogroups

Due to its highly repetitive nature, the human male-specific Y chromosome remains understudied. It is important to investigate variation on the Y chromosome to understand its evolution and contribution to phenotypic variation, including infertility. Approximately 20% of the human Y chromosome consists of ampliconic regions which include nine multi-copy gene families. These gene families are expressed exclusively in testes and usually implicated in spermatogenesis. Here, to gain a better understanding of the role of the Y chromosome in human evolution and in determining sexually dimorphic traits, we studied ampliconic gene copy number variation in 100 males representing ten major Y haplogroups world-wide. Copy number was estimated with droplet digital PCR. In contrast to low nucleotide diversity observed on the Y in previous studies, here we show that ampliconic gene copy number diversity is very high. A total of 98 copy-number-based haplotypes were observed among 100 individuals, and haplotypes were sometimes shared by males from very different haplogroups, suggesting homoplasies. The resulting haplotypes did not cluster according to major Y haplogroups. Overall, only three gene families (DAZ, RBMY, and TSPY) showed significant differences in copy number among major Y haplogroups, and the haplogroup of an individual could not be predicted based on his ampliconic gene copy numbers. Our results suggest rapid evolution of ampliconic gene copy numbers on the human Y. Finally, we tested whether ampliconic gene copy number variation is correlated with variation in height and facial masculinity/femininity, but found no significant associations.