Improving rigor and reproducibility in nonhuman primate research.

Nonhuman primates (NHPs) are a critical component of translational/preclinical biomedical research due to the strong similarities between NHP and human physiology and disease pathology. In some cases, NHPs represent the most appropriate, or even the only, animal model for complex metabolic, neurological, and infectious diseases. The increased demand for and limited availability of these valuable research subjects requires that rigor and reproducibility be a prime consideration to ensure the maximal utility of this scarce resource. Here, we discuss a number of approaches that collectively can contribute to enhanced rigor and reproducibility in NHP research.

Intraspecific and interspecific investigations of skeletal DNA methylation and femur morphology in primates.

OBJECTIVES: Epigenetic mechanisms influence the development and maintenance of complex phenotypes and may also contribute to the evolution of species-specific phenotypes. With respect to skeletal traits, little is known about the gene regulation underlying these hard tissues or how tissue-specific patterns are associated with bone morphology or vary among species. To begin exploring these topics, this study evaluates one epigenetic mechanism, DNA methylation, in skeletal tissues from five nonhuman primate species which display anatomical and locomotor differences representative of their phylogenetic groups. MATERIALS AND METHODS: First, we test whether intraspecific variation in skeletal DNA methylation is associated with intraspecific variation in femur morphology. Second, we identify interspecific differences in DNA methylation and assess whether these lineage-specific patterns may have contributed to species-specific morphologies. Specifically, we use the Illumina Infinium MethylationEPIC BeadChip to identify DNA methylation patterns in femur trabecular bone from baboons (n = 28), macaques (n = 10), vervets (n = 10), chimpanzees (n = 4), and marmosets (n = 6). RESULTS: Significant differentially methylated positions (DMPs) were associated with a subset of morphological variants, but these likely have small biological effects and may be confounded by other variables associated with morphological variation. Conversely, several species-specific DMPs were identified, and these are found in genes enriched for functions associated with complex skeletal traits. DISCUSSION: Overall, these findings reveal that while intraspecific epigenetic variation is not readily associated with skeletal morphology differences, some interspecific epigenetic differences in skeletal tissues exist and may contribute to evolutionarily distinct phenotypes. This work forms a foundation for future explorations of gene regulation and skeletal trait evolution in primates.

Shades of complexity: New perspectives on the evolution and genetic architecture of human skin.

Like many highly variable human traits, more than a dozen genes are known to contribute to the full range of skin color. However, the historical bias in favor of genetic studies in European and European-derived populations has blinded us to the magnitude of pigmentation's complexity. As deliberate efforts are being made to better characterize diverse global populations and new sequencing technologies, better measurement tools, functional assessments, predictive modeling, and ancient DNA analyses become more widely accessible, we are beginning to appreciate how limited our understanding of the genetic bases of human skin color have been. Novel variants in genes not previously linked to pigmentation have been identified and evidence is mounting that there are hundreds more variants yet to be found. Even for genes that have been exhaustively characterized in European populations like MC1R, OCA2, and SLC24A5, research in previously understudied groups is leading to a new appreciation of the degree to which genetic diversity, epistatic interactions, pleiotropy, admixture, global and local adaptation, and cultural practices operate in population-specific ways to shape the genetic architecture of skin color. Furthermore, we are coming to terms with how factors like tanning response and barrier function may also have influenced selection on skin throughout human history. By examining how our knowledge of pigmentation genetics has shifted in the last decade, we can better appreciate how far we have come in understanding human diversity and the still long road ahead for understanding many complex human traits.

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.

Exome sequences of multiplex, multigenerational families reveal schizophrenia risk loci with potential implications for neurocognitive performance.

Schizophrenia is a serious mental illness, involving disruptions in thought and behavior, with a worldwide prevalence of about one percent. Although highly heritable, much of the genetic liability of schizophrenia is yet to be explained. We searched for susceptibility loci in multiplex, multigenerational families affected by schizophrenia, targeting protein-altering variation with in silico predicted functional effects. Exome sequencing was performed on 136 samples from eight European-American families, including 23 individuals diagnosed with schizophrenia or schizoaffective disorder. In total, 11,878 non-synonymous variants from 6,396 genes were tested for their association with schizophrenia spectrum disorders. Pathway enrichment analyses were conducted on gene-based test results, protein-protein interaction (PPI) networks, and epistatic effects. Using a significance threshold of FDR < 0.1, association was detected for rs10941112 (p = 2.1 × 10-5 ; q-value = 0.073) in AMACR, a gene involved in fatty acid metabolism and previously implicated in schizophrenia, with significant cis effects on gene expression (p = 5.5 × 10-4 ), including brain tissue data from the Genotype-Tissue Expression project (minimum p = 6.0 × 10-5 ). A second SNP, rs10378 located in TMEM176A, also shows risk effects in the exome data (p = 2.8 × 10-5 ; q-value = 0.073). PPIs among our top gene-based association results (p < 0.05; n = 359 genes) reveal significant enrichment of genes involved in NCAM-mediated neurite outgrowth (p = 3.0 × 10-5 ), while exome-wide SNP-SNP interaction effects for rs10941112 and rs10378 indicate a potential role for kinase-mediated signaling involved in memory and learning. In conclusion, these association results implicate AMACR and TMEM176A in schizophrenia risk, whose effects may be modulated by genes involved in synaptic plasticity and neurocognitive performance.

The Evolution of Tanning Needs Its Day in the Sun.

Variation in human pigmentation has long been an area of interest in biological anthropology, with the advent of genetic technologies allowing deeper plumbing of its evolutionary history. Genome-wide scans of selection show that pigmentation genes have undergone some of the strongest selection in many geographically distant populations. A variety of hypotheses for the photoprotective effects of melanin have been developed, but these hypotheses, as well as genetic studies, focus nearly exclusively on constitutive (basal) pigmentation levels. Failing to consider the contribution of the ultraviolet radiation (UVR) environment neglects the true interface between humans and our environment. Data drawn largely from dermatology demonstrate that constitutive pigmentation and tanning response are weakly coupled in populations from East Asia and the Americas. This suggests a possible role for persistent, UVR-induced pigmentation as a convergent adaptation akin to the protective effect of constitutive pigmentation. The adaptive potential of tanned skin, particularly in the Americas, where constitutive pigmentation is lower than expected, may fill in an important gap in our understanding of the evolution of skin color.

Natural selection for the Duffy-null allele in the recently admixed people of Madagascar.

While gene flow between distantly related populations is increasingly recognized as a potentially important source of adaptive genetic variation for humans, fully characterized examples are rare. In addition, the role that natural selection for resistance to vivax malaria may have played in the extreme distribution of the protective Duffy-null allele, which is nearly completely fixed in mainland sub-Saharan Africa and absent elsewhere, is controversial. We address both these issues by investigating the evolution of the Duffy-null allele in the Malagasy, a recently admixed population with major ancestry components from both East Asia and mainland sub-Saharan Africa. We used genome-wide genetic data and extensive computer simulations to show that the high frequency of the Duffy-null allele in Madagascar can only be explained in the absence of positive natural selection under extreme demographic scenarios involving high genetic drift. However, the observed genomic single nucleotide polymorphism diversity in the Malagasy is incompatible with such extreme demographic scenarios, indicating that positive selection for the Duffy-null allele best explains the high frequency of the allele in Madagascar. We estimate the selection coefficient to be 0.066. Because vivax malaria is endemic to Madagascar, this result supports the hypothesis that malaria resistance drove fixation of the Duffy-null allele in mainland sub-Saharan Africa.

ALDH2 is associated to alcohol dependence and is the major genetic determinant of "daily maximum drinks" in a GWAS study of an isolated rural Chinese sample.

Alcohol dependence (AD) is a moderately heritable phenotype with a small number of known risk genes mapped via linkage or candidate gene studies. We considered 313 males from among 595 members of documented, extended pedigrees in which AD segregates collected in Northern Hunan Province, China. A joint analysis of both males and females could not be performed as the difference in alcohol consumption variance was too large. Genome-wide association analyses were performed for approximately 300,000 single nucleotide polymorphisms (SNPs). Significant associations found in the ALDH2 region for AD (minimum P = 4.73 × 10(-8)) and two AD-related phenotypes: flushing response (minimum P = 4.75 × 10(-26)) and maximum drinks in a 24-hr period (minimum P = 1.54 × 10(-16)). Association of previous candidate SNP, rs10774610 in CCDC63, was confirmed but resulted from linkage disequilibrium with ALDH2. ALDH2 is strongly associated with flushing response, AD, and maximum drinks in males, with nonsynonymous SNP rs671 explaining 29.2%, 7.9%, and 22.9% of phenotypic variation, respectively, in this sample. When rs671 was considered as a candidate SNP in females, it explained 23.6% of the variation in flushing response, but alcohol consumption rates were too low among females-despite familial enrichment for AD-for an adequate test of association for either AD or maximum drinks. These results support a mediating effect of aldehyde dehydrogenase deficiency on alcohol consumption in males and a secondary, culturally mediated limitation on alcohol consumption by females that should be appropriately modeled in future studies of alcohol consumption in populations where this may be a factor.

Novel QTL at chromosome 6p22 for alcohol consumption: Implications for the genetic liability of alcohol use disorders.

Linkage studies of alcoholism have implicated several chromosome regions, leading to the successful identification of susceptibility genes, including ADH4 and GABRA2 on chromosome 4. Quantitative endophenotypes that are potentially closer to gene action than clinical endpoints offer a means of obtaining more refined linkage signals of genes that predispose alcohol use disorders (AUD). In this study we examine a self-reported measure of the maximum number of drinks consumed in a 24-hr period (abbreviated Max Drinks), a significantly heritable phenotype (h(2) = 0.32 ± 0.05; P = 4.61 × 10(-14)) with a strong genetic correlation with AUD (ρg = 0.99 ± 0.13) for the San Antonio Family Study (n = 1,203). Genome-wide SNPs were analyzed using variance components linkage methods in the program SOLAR, revealing a novel, genome-wide significant QTL (LOD = 4.17; P = 5.85 × 10(-6)) for Max Drinks at chromosome 6p22.3, a region with a number of compelling candidate genes implicated in neuronal function and psychiatric illness. Joint analysis of Max Drinks and AUD status shows that the QTL has a significant non-zero effect on diagnosis (P = 4.04 × 10(-3)), accounting for 8.6% of the total variation. Significant SNP associations for Max Drinks were also identified at the linkage region, including one, rs7761213 (P = 2.14 × 10(-4)), obtained for an independent sample of Chinese families. Thus, our study identifies a potential risk locus for AUD at 6p22.3, with significant pleiotropic effects on the heaviness of alcohol consumption that may not be population specific.

Technical note: comparing von Luschan skin color tiles and modern spectrophotometry for measuring human skin pigmentation.

Prior to the introduction of reflectance spectrophotometry into anthropological field research during the 1950s, human skin color was most commonly classified by visual skin color matching using the von Luschan tiles, a set of 36 standardized, opaque glass tiles arranged in a chromatic scale. Our goal was to establish a conversion formula between the tile-based color matching method and modern reflectance spectrophotometry to make historical and contemporary data comparable. Skin pigmentation measurements were taken on the forehead, inner upper arms, and backs of the hands using both the tiles and a spectrophotometer on 246 participants showing a broad range of skin pigmentation. From these data, a second-order polynomial conversion formula was derived by jackknife analysis to estimate melanin index (M-index) based on tile values. This conversion formula provides a means for comparing modern data to von Luschan tile measurements recorded in historical reports. This is particularly important for populations now extinct, extirpated, or admixed for which tile-based measures of skin pigmentation are the only data available.