Disentangling Signatures of Selection Before and After European Colonization in Latin Americans.

Throughout human evolutionary history, large-scale migrations have led to intermixing (i.e., admixture) between previously separated human groups. Although classical and recent work have shown that studying admixture can yield novel historical insights, the extent to which this process contributed to adaptation remains underexplored. Here, we introduce a novel statistical model, specific to admixed populations, that identifies loci under selection while determining whether the selection likely occurred post-admixture or prior to admixture in one of the ancestral source populations. Through extensive simulations, we show that this method is able to detect selection, even in recently formed admixed populations, and to accurately differentiate between selection occurring in the ancestral or admixed population. We apply this method to genome-wide SNP data of ∼4,000 individuals in five admixed Latin American cohorts from Brazil, Chile, Colombia, Mexico, and Peru. Our approach replicates previous reports of selection in the human leukocyte antigen region that are consistent with selection post-admixture. We also report novel signals of selection in genomic regions spanning 47 genes, reinforcing many of these signals with an alternative, commonly used local-ancestry-inference approach. These signals include several genes involved in immunity, which may reflect responses to endemic pathogens of the Americas and to the challenge of infectious disease brought by European contact. In addition, some of the strongest signals inferred to be under selection in the Native American ancestral groups of modern Latin Americans overlap with genes implicated in energy metabolism phenotypes, plausibly reflecting adaptations to novel dietary sources available in the Americas.

Population-Based Surveillance for Birth Defects Potentially Related to Zika Virus Infection - 22 States and Territories, January 2016-June 2017.

Zika virus infection during pregnancy can cause congenital brain and eye abnormalities and is associated with neurodevelopmental abnormalities (1-3). In areas of the United States that experienced local Zika virus transmission, the prevalence of birth defects potentially related to Zika virus infection during pregnancy increased in the second half of 2016 compared with the first half (4). To update the previous report, CDC analyzed population-based surveillance data from 22 states and territories to estimate the prevalence of birth defects potentially related to Zika virus infection, regardless of laboratory evidence of or exposure to Zika virus, among pregnancies completed during January 1, 2016-June 30, 2017. Jurisdictions were categorized as those 1) with widespread local transmission of Zika virus; 2) with limited local transmission of Zika virus; and 3) without local transmission of Zika virus. Among 2,004,630 live births, 3,359 infants and fetuses with birth defects potentially related to Zika virus infection during pregnancy were identified (1.7 per 1,000 live births, 95% confidence interval [CI] = 1.6-1.7). In areas with widespread local Zika virus transmission, the prevalence of birth defects potentially related to Zika virus infection during pregnancy was significantly higher during the quarters comprising July 2016-March 2017 (July-September 2016 = 3.0; October-December 2016 = 4.0; and January-March 2017 = 5.6 per 1,000 live births) compared with the reference period (January-March 2016) (1.3 per 1,000). These findings suggest a fourfold increase (prevalence ratio [PR] = 4.1, 95% CI = 2.1-8.4) in birth defects potentially related to Zika virus in widespread local transmission areas during January-March 2017 compared with that during January-March 2016, with the highest prevalence (7.0 per 1,000 live births) in February 2017. Population-based birth defects surveillance is critical for identifying infants and fetuses with birth defects potentially related to Zika virus regardless of whether Zika virus testing was conducted, especially given the high prevalence of asymptomatic disease. These data can be used to inform follow-up care and services as well as strengthen surveillance.

Population-Based Surveillance of Birth Defects Potentially Related to Zika Virus Infection - 15 States and U.S. Territories, 2016.

Zika virus infection during pregnancy can cause serious birth defects, including microcephaly and brain abnormalities (1). Population-based birth defects surveillance systems are critical to monitor all infants and fetuses with birth defects potentially related to Zika virus infection, regardless of known exposure or laboratory evidence of Zika virus infection during pregnancy. CDC analyzed data from 15 U.S. jurisdictions conducting population-based surveillance for birth defects potentially related to Zika virus infection.* Jurisdictions were stratified into the following three groups: those with 1) documented local transmission of Zika virus during 2016; 2) one or more cases of confirmed, symptomatic, travel-associated Zika virus disease reported to CDC per 100,000 residents; and 3) less than one case of confirmed, symptomatic, travel-associated Zika virus disease reported to CDC per 100,000 residents. A total of 2,962 infants and fetuses (3.0 per 1,000 live births; 95% confidence interval [CI] = 2.9-3.2) (2) met the case definition.† In areas with local transmission there was a non-statistically significant increase in total birth defects potentially related to Zika virus infection from 2.8 cases per 1,000 live births in the first half of 2016 to 3.0 cases in the second half (p = 0.10). However, when neural tube defects and other early brain malformations (NTDs)§ were excluded, the prevalence of birth defects strongly linked to congenital Zika virus infection increased significantly, from 2.0 cases per 1,000 live births in the first half of 2016 to 2.4 cases in the second half, an increase of 29 more cases than expected (p = 0.009). These findings underscore the importance of surveillance for birth defects potentially related to Zika virus infection and the need for continued monitoring in areas at risk for Zika.

Language Deficits as a Preclinical Window into Parkinson's Disease: Evidence from Asymptomatic Parkin and Dardarin Mutation Carriers.

OBJECTIVES: The worldwide spread of Parkinson's disease (PD) calls for sensitive and specific measures enabling its early (or, ideally, preclinical) detection. Here, we use language measures revealing deficits in PD to explore whether similar disturbances are present in asymptomatic individuals at risk for the disease. METHODS: We administered executive, semantic, verb-production, and syntactic tasks to sporadic PD patients, genetic PD patients with PARK2 (parkin) or LRRK2 (dardarin) mutation, asymptomatic first-degree relatives of the latter with similar mutations, and socio-demographically matched controls. Moreover, to detect sui generis language disturbances, we ran analysis of covariance tests using executive functions as covariate. RESULTS: The two clinical groups showed impairments in all measures, most of which survived covariation with executive functions. However, the key finding concerned asymptomatic mutation carriers. While these subjects showed intact executive, semantic, and action-verb production skills, they evinced deficits in a syntactic test with minimal working memory load. CONCLUSIONS: We propose that this sui generis disturbance may constitute a prodromal sign anticipating eventual development of PD. Moreover, our results suggest that mutations on specific genes (PARK2 and LRRK2) compromising basal ganglia functioning may be subtly related to language-processing mechanisms. (JINS, 2017, 23, 150-158).

Admixture in Latin America: geographic structure, phenotypic diversity and self-perception of ancestry based on 7,342 individuals.

The current genetic makeup of Latin America has been shaped by a history of extensive admixture between Africans, Europeans and Native Americans, a process taking place within the context of extensive geographic and social stratification. We estimated individual ancestry proportions in a sample of 7,342 subjects ascertained in five countries (Brazil, Chile, Colombia, México and Perú). These individuals were also characterized for a range of physical appearance traits and for self-perception of ancestry. The geographic distribution of admixture proportions in this sample reveals extensive population structure, illustrating the continuing impact of demographic history on the genetic diversity of Latin America. Significant ancestry effects were detected for most phenotypes studied. However, ancestry generally explains only a modest proportion of total phenotypic variation. Genetically estimated and self-perceived ancestry correlate significantly, but certain physical attributes have a strong impact on self-perception and bias self-perception of ancestry relative to genetically estimated ancestry.

Increasing Prevalence of Gastroschisis--14 States, 1995-2012.

Gastroschisis is a serious congenital defect in which the intestines protrude through an opening in the abdominal wall. Gastroschisis requires surgical repair soon after birth and is associated with an increased risk for medical complications and mortality during infancy. Reports from multiple surveillance systems worldwide have documented increasing prevalence of gastroschisis since the 1980s, particularly among younger mothers; however, since publication of a multistate U.S. report that included data through 2005, it is not known whether prevalence has continued to increase. Data on gastroschisis from 14 population-based state surveillance programs were pooled and analyzed to assess the average annual percent change (AAPC) in prevalence and to compare the prevalence during 2006-2012 with that during 1995-2005, stratified by maternal age and race/ethnicity. The pooled data included approximately 29% of U.S. births for the period 1995-2012. During 1995-2012, gastroschisis prevalence increased in every category of maternal age and race/ethnicity, and the AAPC ranged from 3.1% in non-Hispanic white (white) mothers aged <20 years to 7.9% in non-Hispanic black (black) mothers aged <20 years. These corresponded to overall percentage increases during 1995-2012 that ranged from 68% in white mothers aged <20 years to 263% in black mothers aged <20 years. Gastroschisis prevalence increased 30% between the two periods, from 3.6 per 10,000 births during 1995-2005 to 4.9 per 10,000 births during 2006-2012 (prevalence ratio = 1.3, 95% confidence interval [CI]: 1.3-1.4), with the largest increase among black mothers aged <20 years (prevalence ratio = 2.0, 95% CI: 1.6-2.5). Public health research is urgently needed to identify factors contributing to this increase.

Neanderthal introgression in SCN9A impacts mechanical pain sensitivity.

The Nav1.7 voltage-gated sodium channel plays a key role in nociception. Three functional variants in the SCN9A gene (encoding M932L, V991L, and D1908G in Nav1.7), have recently been identified as stemming from Neanderthal introgression and to associate with pain symptomatology in UK BioBank data. In 1000 genomes data, these variants are absent in Europeans but common in Latin Americans. Analysing high-density genotype data from 7594 Latin Americans, we characterized Neanderthal introgression in SCN9A. We find that tracts of introgression occur on a Native American genomic background, have an average length of ~123 kb and overlap the M932L, V991L, and D1908G coding positions. Furthermore, we measured experimentally six pain thresholds in 1623 healthy Colombians. We found that Neanderthal ancestry in SCN9A is significantly associated with a lower mechanical pain threshold after sensitization with mustard oil and evidence of additivity of effects across Nav1.7 variants. Our findings support the reported association of Neanderthal Nav1.7 variants with clinical pain, define a specific sensory modality affected by archaic introgression in SCN9A and are consistent with independent effects of the Neanderthal variants on Nav1.7 function.

Automatic landmarking identifies new loci associated with face morphology and implicates Neanderthal introgression in human nasal shape.

We report a genome-wide association study of facial features in >6000 Latin Americans based on automatic landmarking of 2D portraits and testing for association with inter-landmark distances. We detected significant associations (P-value <5 × 10-8) at 42 genome regions, nine of which have been previously reported. In follow-up analyses, 26 of the 33 novel regions replicate in East Asians, Europeans, or Africans, and one mouse homologous region influences craniofacial morphology in mice. The novel region in 1q32.3 shows introgression from Neanderthals and we find that the introgressed tract increases nasal height (consistent with the differentiation between Neanderthals and modern humans). Novel regions include candidate genes and genome regulatory elements previously implicated in craniofacial development, and show preferential transcription in cranial neural crest cells. The automated approach used here should simplify the collection of large study samples from across the world, facilitating a cosmopolitan characterization of the genetics of facial features.

Prediction of eye, hair and skin colour in Latin Americans.

Here we evaluate the accuracy of prediction for eye, hair and skin pigmentation in a dataset of > 6500 individuals from Mexico, Colombia, Peru, Chile and Brazil (including genome-wide SNP data and quantitative/categorical pigmentation phenotypes - the CANDELA dataset CAN). We evaluated accuracy in relation to different analytical methods and various phenotypic predictors. As expected from statistical principles, we observe that quantitative traits are more sensitive to changes in the prediction models than categorical traits. We find that Random Forest or Linear Regression are generally the best performing methods. We also compare the prediction accuracy of SNP sets defined in the CAN dataset (including 56, 101 and 120 SNPs for eye, hair and skin colour prediction, respectively) to the well-established HIrisPlex-S SNP set (including 6, 22 and 36 SNPs for eye, hair and skin colour prediction respectively). When training prediction models on the CAN data, we observe remarkably similar performances for HIrisPlex-S and the larger CAN SNP sets for the prediction of hair (categorical) and eye (both categorical and quantitative), while the CAN sets outperform HIrisPlex-S for quantitative, but not for categorical skin pigmentation prediction. The performance of HIrisPlex-S, when models are trained in a world-wide sample (although consisting of 80% Europeans, https://hirisplex.erasmusmc.nl), is lower relative to training in the CAN data (particularly for hair and skin colour). Altogether, our observations are consistent with common variation of eye and hair colour having a relatively simple genetic architecture, which is well captured by HIrisPlex-S, even in admixed Latin Americans (with partial European ancestry). By contrast, since skin pigmentation is a more polygenic trait, accuracy is more sensitive to prediction SNP set size, although here this effect was only apparent for a quantitative measure of skin pigmentation. Our results support the use of HIrisPlex-S in the prediction of categorical pigmentation traits for forensic purposes in Latin America, while illustrating the impact of training datasets on its accuracy.

A GWAS in Latin Americans identifies novel face shape loci, implicating VPS13B and a Denisovan introgressed region in facial variation.

To characterize the genetic basis of facial features in Latin Americans, we performed a genome-wide association study (GWAS) of more than 6000 individuals using 59 landmark-based measurements from two-dimensional profile photographs and ~9,000,000 genotyped or imputed single-nucleotide polymorphisms. We detected significant association of 32 traits with at least 1 (and up to 6) of 32 different genomic regions, more than doubling the number of robustly associated face morphology loci reported until now (from 11 to 23). These GWAS hits are strongly enriched in regulatory sequences active specifically during craniofacial development. The associated region in 1p12 includes a tract of archaic adaptive introgression, with a Denisovan haplotype common in Native Americans affecting particularly lip thickness. Among the nine previously unidentified face morphology loci we identified is the VPS13B gene region, and we show that variants in this region also affect midfacial morphology in mice.

Genetic make up and structure of Colombian populations by means of uniparental and biparental DNA markers.

Colombia is a country with great geographic heterogeneity and marked regional differences in pre-Columbian native population density and in the extent of past African and European immigration. As a result, Colombia has one of the most diverse populations in Latin America. Here we evaluated ancestry in over 1,700 individuals from 24 Colombian populations using biparental (autosomal and X-Chromosome), maternal (mtDNA), and paternal (Y-chromosome) markers. Autosomal ancestry varies markedly both within and between regions, confirming the great genetic diversity of the Colombian population. The X-chromosome, mtDNA, and Y-chromosome data indicate that there is a pattern across regions indicative of admixture involving predominantly Native American women and European and African men.

A GWAS in Latin Americans highlights the convergent evolution of lighter skin pigmentation in Eurasia.

We report a genome-wide association scan in >6,000 Latin Americans for pigmentation of skin and eyes. We found eighteen signals of association at twelve genomic regions. These include one novel locus for skin pigmentation (in 10q26) and three novel loci for eye pigmentation (in 1q32, 20q13 and 22q12). We demonstrate the presence of multiple independent signals of association in the 11q14 and 15q13 regions (comprising the GRM5/TYR and HERC2/OCA2 genes, respectively) and several epistatic interactions among independently associated alleles. Strongest association with skin pigmentation at 19p13 was observed for an Y182H missense variant (common only in East Asians and Native Americans) in MFSD12, a gene recently associated with skin pigmentation in Africans. We show that the frequency of the derived allele at Y182H is significantly correlated with lower solar radiation intensity in East Asia and infer that MFSD12 was under selection in East Asians, probably after their split from Europeans.

Non-smoking pregnant women and their fetuses are exposed to environmental tobacco smoke as a result of living in multiunit housing.

This study investigates whether pregnant women and their fetuses are exposed to environmental tobacco smoke (ETS) as a result of living in apartments. We measured cotinine concentrations in serum, a biomarker of exposure to ETS, in non-smoking women's umbilical cord blood collected at delivery and in maternal blood drawn shortly after delivering a baby. Concurrently, information was collected regarding the women's housing situation, whether family members or co-workers smoked, and other potential exposure factors. Newborns whose non-smoking mothers lived in an apartment during pregnancy were more than three times (OR 3.17, 95% CI 1.62-6.21) more likely to have detectable levels of cotinine in their cord blood serum than babies whose mothers lived in a detached house. There is a strong association between detectable concentrations of cotinine in cord blood serum and living in an apartment, even after adjusting for confounders, such as exposure at home or at work. A similar association was observed between the detectable levels of cotinine in maternal serum and living in an apartment (OR 1.95, 95% CI 1.03-3.71).

A genome-wide association scan in admixed Latin Americans identifies loci influencing facial and scalp hair features.

We report a genome-wide association scan in over 6,000 Latin Americans for features of scalp hair (shape, colour, greying, balding) and facial hair (beard thickness, monobrow, eyebrow thickness). We found 18 signals of association reaching genome-wide significance (P values 5 × 10(-8) to 3 × 10(-119)), including 10 novel associations. These include novel loci for scalp hair shape and balding, and the first reported loci for hair greying, monobrow, eyebrow and beard thickness. A newly identified locus influencing hair shape includes a Q30R substitution in the Protease Serine S1 family member 53 (PRSS53). We demonstrate that this enzyme is highly expressed in the hair follicle, especially the inner root sheath, and that the Q30R substitution affects enzyme processing and secretion. The genome regions associated with hair features are enriched for signals of selection, consistent with proposals regarding the evolution of human hair.

A genome-wide association scan implicates DCHS2, RUNX2, GLI3, PAX1 and EDAR in human facial variation.

We report a genome-wide association scan for facial features in ∼6,000 Latin Americans. We evaluated 14 traits on an ordinal scale and found significant association (P values<5 × 10(-8)) at single-nucleotide polymorphisms (SNPs) in four genomic regions for three nose-related traits: columella inclination (4q31), nose bridge breadth (6p21) and nose wing breadth (7p13 and 20p11). In a subsample of ∼3,000 individuals we obtained quantitative traits related to 9 of the ordinal phenotypes and, also, a measure of nasion position. Quantitative analyses confirmed the ordinal-based associations, identified SNPs in 2q12 associated to chin protrusion, and replicated the reported association of nasion position with SNPs in PAX3. Strongest association in 2q12, 4q31, 6p21 and 7p13 was observed for SNPs in the EDAR, DCHS2, RUNX2 and GLI3 genes, respectively. Associated SNPs in 20p11 extend to PAX1. Consistent with the effect of EDAR on chin protrusion, we documented alterations of mandible length in mice with modified Edar funtion.

A genome-wide association study identifies multiple loci for variation in human ear morphology.

Here we report a genome-wide association study for non-pathological pinna morphology in over 5,000 Latin Americans. We find genome-wide significant association at seven genomic regions affecting: lobe size and attachment, folding of antihelix, helix rolling, ear protrusion and antitragus size (linear regression P values 2 × 10(-8) to 3 × 10(-14)). Four traits are associated with a functional variant in the Ectodysplasin A receptor (EDAR) gene, a key regulator of embryonic skin appendage development. We confirm expression of Edar in the developing mouse ear and that Edar-deficient mice have an abnormally shaped pinna. Two traits are associated with SNPs in a region overlapping the T-Box Protein 15 (TBX15) gene, a major determinant of mouse skeletal development. Strongest association in this region is observed for SNP rs17023457 located in an evolutionarily conserved binding site for the transcription factor Cartilage paired-class homeoprotein 1 (CART1), and we confirm that rs17023457 alters in vitro binding of CART1.