Social genomics, cognition, and well-being during the COVID-19 pandemic.

INTRODUCTION: Adverse psychosocial exposure is associated with increased proinflammatory gene expression and reduced type-1 interferon gene expression, a profile known as the conserved transcriptional response to adversity (CTRA). Little is known about CTRA activity in the context of cognitive impairment, although chronic inflammatory activation has been posited as one mechanism contributing to late-life cognitive decline. METHODS: We studied 171 community-dwelling older adults from the Wake Forest Alzheimer's Disease Research Center who answered questions via a telephone questionnaire battery about their perceived stress, loneliness, well-being, and impact of COVID-19 on their life, and who provided a self-collected dried blood spot sample. Of those, 148 had adequate samples for mRNA analysis, and 143 were included in the final analysis, which including participants adjudicated as having normal cognition (NC, n = 91) or mild cognitive impairment (MCI, n = 52) were included in the analysis. Mixed effect linear models were used to quantify associations between psychosocial variables and CTRA gene expression. RESULTS: In both NC and MCI groups, eudaimonic well-being (typically associated with a sense of purpose) was inversely associated with CTRA gene expression whereas hedonic well-being (typically associated with pleasure seeking) was positively associated. In participants with NC, coping through social support was associated with lower CTRA gene expression, whereas coping by distraction and reframing was associated with higher CTRA gene expression. CTRA gene expression was not related to coping strategies for participants with MCI, or to either loneliness or perceived stress in either group. DISCUSSION: Eudaimonic and hedonic well-being remain important correlates of molecular markers of stress, even in people with MCI. However, prodromal cognitive decline appears to moderate the significance of coping strategies as a correlate of CTRA gene expression. These results suggest that MCI can selectively alter biobehavioral interactions in ways that could potentially affect the rate of future cognitive decline and may serve as targets for future intervention efforts.

Resilience, aging, and response to radiation exposure (RARRE) in nonhuman primates: a resource review.

The Wake Forest nonhuman primate (NHP) Radiation Late Effects Cohort (RLEC) is a unique and irreplaceable population of aging NHP radiation survivors which serves the nation's need to understand the late effects of radiation exposure. Over the past 16 years, Wake Forest has evaluated > 250 previously irradiated rhesus macaques (Macaca mulatta) that were exposed to single total body irradiation (IR) doses of 1.14-8.5 Gy or to partial body exposures of up to 10 Gy (5% bone marrow sparing) or 10.75 Gy (whole thorax). Though primarily used to examine IR effects on disease-specific processes or to develop radiation countermeasures, this resource provides insights on resilience across physiologic systems and its relationship with biological aging. Exposure to IR has well documented deleterious effects on health, but the late effects of IR are highly variable. Some animals exhibit multimorbidity and accumulated health deficits, whereas others remain relatively resilient years after exposure to total body IR. This provides an opportunity to evaluate biological aging at the nexus of resilient/vulnerable responses to a stressor. Consideration of inter-individual differences in response to this stressor can inform individualized strategies to manage late effects of radiation exposure, and provide insight into mechanisms underlying systemic resilience and aging. The utility of this cohort for age-related research questions was summarized at the 2022 Trans-NIH Geroscience Interest Group's Workshop on Animal Models for Geroscience. We present a brief review of radiation injury and its relationship to aging and resilience in NHPs with a focus on the RLEC.

A New Method of Bone Stromal Cell Characterization by Flow Cytometry.

The bone microenvironment cellular composition plays an essential role in bone health and is disrupted in bone pathologies, such as osteoporosis, osteoarthritis, and cancer. Flow cytometry protocols for hematopoietic stem cell lineages are well defined and well established. Additionally, a consensus for mesenchymal stem cell flow markers has been developed. However, flow cytometry markers for bone-residing cells-osteoblasts, osteoclasts, and osteocytes-have not been proposed. Here, we describe a novel partial digestion method to separate these cells from the bone matrix and present new markers for enumerating these cells by flow cytometry. We optimized bone digestion and analyzed markers across murine, nonhuman primate, and human bone. The isolation and staining protocols can be used with either cell sorting or flow cytometry. Our method allows for the enumeration and collection of hematopoietic and mesenchymal lineage cells in the bone microenvironment combined with bone-residing stromal cells. Thus, we have established a multi-fluorochrome bone marrow cell-typing methodology. © 2022 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Partial digestion for murine long bone stromal cell isolation Alternate Protocol 1: Partial digestion for primate vertebrae stromal cell isolation Alternate Protocol 2: Murine vertebrae crushing for bone stromal cell isolation Basic Protocol 2: Staining of bone stromal cells Support Protocol 1: Fluorescence minus one control, isotype control, and antibody titration Basic Protocol 3: Cell sorting of bone stromal cells Alternate Protocol 3: Flow cytometry analysis of bone stromal cells Support Protocol 2: Preparing compensation beads.

Circulating miRNAs associated with bone mineral density in healthy adult baboons.

MicroRNAs (miRNAs) regulate gene expression post-transcriptionally and circulate in the blood, making them attractive biomarkers of disease state for tissues like bone that are challenging to interrogate directly. Here, we report on five miRNAs-miR-197-3p, miR-320a, miR-320b, miR-331-5p, and miR-423-5p-associated with bone mineral density (BMD) in 147 healthy adult baboons. These baboons ranged in age from 15 to 25 years (45-75 human equivalent years) and 65% were female with a broad range of BMD values including a minority of osteopenic animals. miRNAs were generated via RNA sequencing from buffy coats collected at necropsy and areal BMD (aBMD) measured postmortem via dual-energy X-ray absorptiometry (DXA) of the lumbar vertebrae. Differential expression analysis controlled for the underlying pedigree structure of these animals to account for genetic variation which may drive miRNA abundance and aBMD values. While many of these miRNAs have been associated with the risk of osteoporosis in humans, this finding is of interest because the cohort represents a model of normal aging and bone metabolism rather than a disease cohort. The replication of miRNA associations with osteoporosis or other bone metabolic disorders in animals with healthy aBMD suggests an overlap in normal variation and disease states. We suggest that these miRNAs are involved in the regulation of cellular proliferation, apoptosis, and protein composition in the extracellular matrix throughout life; and age-related dysregulation of these systems may lead to disease. These miRNAs may be early indicators of progression to disease in advance of clinically detectible osteoporosis.

Integrated omics analysis reveals sirtuin signaling is central to hepatic response to a high fructose diet.

BACKGROUND: Dietary high fructose (HFr) is a known metabolic disruptor contributing to development of obesity and diabetes in Western societies. Initial molecular changes from exposure to HFr on liver metabolism may be essential to understand the perturbations leading to insulin resistance and abnormalities in lipid and carbohydrate metabolism. We studied vervet monkeys (Clorocebus aethiops sabaeus) fed a HFr (n=5) or chow diet (n=5) for 6 weeks, and obtained clinical measures of liver function, blood insulin, cholesterol and triglycerides. In addition, we performed untargeted global transcriptomics, proteomics, and metabolomics analyses on liver biopsies to determine the molecular impact of a HFr diet on coordinated pathways and networks that differed by diet. RESULTS: We show that integration of omics data sets improved statistical significance for some pathways and networks, and decreased significance for others, suggesting that multiple omics datasets enhance confidence in relevant pathway and network identification. Specifically, we found that sirtuin signaling and a peroxisome proliferator activated receptor alpha (PPARA) regulatory network were significantly altered in hepatic response to HFr. Integration of metabolomics and miRNAs data further strengthened our findings. CONCLUSIONS: Our integrated analysis of three types of omics data with pathway and regulatory network analysis demonstrates the usefulness of this approach for discovery of molecular networks central to a biological response. In addition, metabolites aspartic acid and docosahexaenoic acid (DHA), protein ATG3, and genes ATG7, and HMGCS2 link sirtuin signaling and the PPARA network suggesting molecular mechanisms for altered hepatic gluconeogenesis from consumption of a HFr diet.

Heterozygosity of the major histocompatibility complex predicts later self-reported pubertal maturation in men.

Individual variation in the age of pubertal onset is linked to physical and mental health, yet the factors underlying this variation are poorly understood. Life history theory predicts that individuals at higher risk of mortality due to extrinsic causes such as infectious disease should sexually mature and reproduce earlier, whereas those at lower risk can delay puberty and continue to invest resources in somatic growth. We examined relationships between a genetic predictor of infectious disease resistance, heterozygosity of the major histocompatibility complex (MHC), referred to as the human leukocyte antigen (HLA) gene in humans, and self-reported pubertal timing. In a combined sample of men from Canada (n = 137) and the United States (n = 43), MHC heterozygosity predicted later self-reported pubertal development. These findings suggest a genetic trade-off between immunocompetence and sexual maturation in human males.

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.

Dichloroacetate reverses sepsis-induced hepatic metabolic dysfunction.

Metabolic reprogramming between resistance and tolerance occurs within the immune system in response to sepsis. While metabolic tissues such as the liver are subjected to damage during sepsis, how their metabolic and energy reprogramming ensures survival is unclear. Employing comprehensive metabolomic, lipidomic, and transcriptional profiling in a mouse model of sepsis, we show that hepatocyte lipid metabolism, mitochondrial tricarboxylic acid (TCA) energetics, and redox balance are significantly reprogrammed after cecal ligation and puncture (CLP). We identify increases in TCA cycle metabolites citrate, cis-aconitate, and itaconate with reduced fumarate and triglyceride accumulation in septic hepatocytes. Transcriptomic analysis of liver tissue supports and extends the hepatocyte findings. Strikingly, the administration of the pyruvate dehydrogenase kinase (PDK) inhibitor dichloroacetate reverses dysregulated hepatocyte metabolism and mitochondrial dysfunction. In summary, our data indicate that sepsis promotes hepatic metabolic dysfunction and that targeting the mitochondrial PDC/PDK energy homeostat rebalances transcriptional and metabolic manifestations of sepsis within the liver.

An evolutionary perspective of DNA methylation patterns in skeletal tissues using a baboon model of osteoarthritis.

Epigenetic factors, such as DNA methylation, play an influential role in the development of the degenerative joint disease osteoarthritis (OA). These molecular mechanisms have been heavily studied in humans, and although OA affects several other animals in addition to humans, few efforts have taken an evolutionary perspective. This study explores the evolution of OA epigenetics by assessing the relationship between DNA methylation variation and knee OA development in baboons (Papio spp.) and by comparing these findings to human OA epigenetic associations. Genome-wide DNA methylation patterns were identified in bone and cartilage of the right distal femora from 56 pedigreed, adult baboons (28 with and 28 without knee OA) using the Illumina Infinium MethylationEPIC BeadChip. Several significantly differentially methylated positions (DMPs) and regions were found between tissue types. Substantial OA-related differential methylation was also identified in cartilage, but not in bone, suggesting that cartilage epigenetics may be more influential in OA than bone epigenetics. Additionally, some genes containing OA-related DMPs overlap with and display methylation patterns similar to those previously identified in human OA, revealing a mixture of evolutionarily conserved and divergent OA-related methylation patterns in primates. Overall, these findings reinforce the current etiological perspectives of OA and enhance our evolutionary understanding of epigenetic mechanisms associated with OA. This study further establishes baboons as a valuable nonhuman primate model of OA, and continued investigations in baboons will help to disentangle the molecular mechanisms contributing to OA and their evolutionary histories.

Insights into the genetic architecture of the human face.

The human face is complex and multipartite, and characterization of its genetic architecture remains challenging. Using a multivariate genome-wide association study meta-analysis of 8,246 European individuals, we identified 203 genome-wide-significant signals (120 also study-wide significant) associated with normal-range facial variation. Follow-up analyses indicate that the regions surrounding these signals are enriched for enhancer activity in cranial neural crest cells and craniofacial tissues, several regions harbor multiple signals with associations to different facial phenotypes, and there is evidence for potential coordinated actions of variants. In summary, our analyses provide insights into the understanding of how complex morphological traits are shaped by both individual and coordinated genetic actions.

Use of Metabolomic Profiling to Understand Variability in Adiposity Changes Following an Intentional Weight Loss Intervention in Older Adults.

Inter-individual response to dietary interventions remains a major challenge to successful weight loss among older adults. This study applied metabolomics technology to identify small molecule signatures associated with a loss of fat mass and overall weight in a cohort of older adults on a nutritionally complete, high-protein diet. A total of 102 unique metabolites were measured using liquid chromatography-mass spectrometry (LC-MS) for 38 adults aged 65-80 years randomized to dietary intervention and 36 controls. Metabolite values were analyzed in both baseline plasma samples and samples collected following the six-month dietary intervention to consider both metabolites that could predict the response to diet and those that changed in response to diet or weight loss.Eight metabolites changed over the intervention at a nominally significant level: D-pantothenic acid, L-methionine, nicotinate, aniline, melatonin, deoxycarnitine, 6-deoxy-L-galactose, and 10-hydroxydecanoate. Within the intervention group, there was broad variation in the achieved weight-loss and dual-energy x-ray absorptiometry (DXA)-defined changes in total fat and visceral adipose tissue (VAT) mass. Change in the VAT mass was significantly associated with the baseline abundance of α-aminoadipate (p = 0.0007) and an additional mass spectrometry peak that may represent D-fructose, myo-inositol, mannose, α-D-glucose, allose, D-galactose, D-tagatose, or L-sorbose (p = 0.0001). This hypothesis-generating study reflects the potential of metabolomic biomarkers for the development of personalized dietary interventions.

Differential DNA methylation of vocal and facial anatomy genes in modern humans.

Changes in potential regulatory elements are thought to be key drivers of phenotypic divergence. However, identifying changes to regulatory elements that underlie human-specific traits has proven very challenging. Here, we use 63 reconstructed and experimentally measured DNA methylation maps of ancient and present-day humans, as well as of six chimpanzees, to detect differentially methylated regions that likely emerged in modern humans after the split from Neanderthals and Denisovans. We show that genes associated with face and vocal tract anatomy went through particularly extensive methylation changes. Specifically, we identify widespread hypermethylation in a network of face- and voice-associated genes (SOX9, ACAN, COL2A1, NFIX and XYLT1). We propose that these repression patterns appeared after the split from Neanderthals and Denisovans, and that they might have played a key role in shaping the modern human face and vocal tract.

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.

Assessment of DNA Methylation Patterns in the Bone and Cartilage of a Nonhuman Primate Model of Osteoarthritis.

OBJECTIVE: Osteoarthritis (OA) affects humans and several other animals. Thus, the mechanisms underlying this disorder, such as specific skeletal tissue DNA methylation patterns, may be evolutionary conserved. However, associations between methylation and OA have not been readily studied in nonhuman animals. Baboons serve as important models of disease and develop OA at rates similar to those in humans. Therefore, this study investigated the associations between methylation and OA in baboons to advance the evolutionary understanding of OA. DESIGN: Trabecular bone and cartilage was collected from the medial condyles of adult female baboon femora, 5 with and 5 without knee OA. The Infinium HumanMethylation450 BeadChip (450K array) was used to identify DNA methylation patterns in these tissues. RESULTS: Approximately 44% of the 450K array probes reliably align to the baboon genome, contain a CpG site of interest, and maintain a wide distribution throughout the genome. Of the 2 filtering methods tested, both identified significantly differentially methylated positions (DMPs) between healthy and OA individuals in cartilage tissues, and some of these patterns overlap with those previously identified in humans. Conversely, no DMPs were found between tissue types or between disease states in bone tissues. CONCLUSIONS: Overall, the 450K array can be used to measure genome-wide DNA methylation in baboon tissues and identify significant associations with complex traits. The results of this study indicate that some DNA methylation patterns associated with OA are evolutionarily conserved, while others are not. This warrants further investigation in a larger and more phylogenetically diverse sample set.

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.

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.

A variance component method for integrated pathway analysis of gene expression data.

BACKGROUND: The application of pathway and gene-set based analyses to high-throughput data is increasingly common and represents an effort to understand underlying biology where single-gene or single-marker analyses have failed. Many such analyses rely on the a priori identification of genes associated with the trait of interest. In contrast, this variance-component-based approach creates a similarity matrix of individuals based on the expression of genes in each pathway. METHODS: We compared 16 methods of calculating similarity for positive control matrices based on probes for the genes used to model the simulated Genetic Analysis Workshop phenotypes. RESULTS: A simple correlation matrix outperforms the other methods by identifying pathways associated with the simulated phenotypes at nearly twice the rate expected based on the associations of the component transcripts and an approximate false-positive rate of 0.05. CONCLUSIONS: This method has a number of additional advantages compared to single-transcript and pathway overrepresentation analyses, including the ability to estimate the proportion of variation explained by each pathway and the logistical advantage of only calculating the distance matrices once for each messenger RNA data set regardless of the number of phenotypes. Additionally, it offers a significant reduction in the multiple testing burden over individual consideration of each probe.

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.

Evaluation of estimated genetic values and their application to genome-wide investigation of systolic blood pressure.

The concept of breeding values, an individual's phenotypic deviation from the population mean as a result of the sum of the average effects of the genes they carry, is of great importance in livestock, aquaculture, and cash crop industries where emphasis is placed on an individual's potential to pass desirable phenotypes on to the next generation. As breeding or genetic values (as referred to here) cannot be measured directly, estimated genetic values (EGVs) are based on an individual's own phenotype, phenotype information from relatives, and, increasingly, genetic data. Because EGVs represent additive genetic variation, calculating EGVs in an extended human pedigree is expected to provide a more refined phenotype for genetic analyses. To test the utility of EGVs in genome-wide association, EGVs were calculated for 847 members of 20 extended Mexican American families based on 100 replicates of simulated systolic blood pressure. Calculations were performed in GAUSS to solve a variation on the standard Best Linear Unbiased Predictor (BLUP) mixed model equation with age, sex, and the first 3 principal components of sample-wide genetic variability as fixed effects and the EGV as a random effect distributed around the relationship matrix. Three methods of calculating kinship were considered: expected kinship from pedigree relationships, empirical kinship from common variants, and empirical kinship from both rare and common variants. Genome-wide association analysis was conducted on simulated phenotypes and EGVs using the additive measured genotype approach in the SOLAR software package. The EGV-based approach showed only minimal improvement in power to detect causative loci.

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.