Ancestral components of admixed genomes in a Mexican cohort

Fast, accurate local ancestry inference with FLARE

Local ancestry is the source ancestry at each point in the genome of an admixed individual. Inferred local ancestry is used for admixture mapping and population genetic analyses. We present FLARE (fast local ancestry estimation), a method for local ancestry inference. FLARE achieves high accuracy through the use of an extended Li and Stephens model, and it achieves exceptional computational performance through incorporation of computational techniques developed for genotype imputation. Memory requirements are reduced through on-the-fly compression of reference haplotypes and stored checkpoints. Computation time is reduced through the use of composite reference haplotypes. These techniques allow FLARE to scale to datasets with hundreds of thousands of sequenced individuals and to provide superior accuracy on large-scale data. FLARE is open source and available at https://github.com/browning-lab/flare.

Ancestry-dependent genetic structure of the Xq28 risk haplotype in the Mexican population and its association with childhood-onset systemic lupus erythematosus

Objective

Here we aimed to investigate the association of the Xq28 risk haplotype (H1) with susceptibility to childhood-onset systemic lupus erythematosus (SLE), and to compare its frequency and genetic structure in the Mexican population with those in other continental populations.

Methods

We genotyped 15 single-nucleotide variants (SNVs) that form the H1 haplotype, using TaqMan real-time PCR. The association analysis [case-control and transmission disequilibrium test (TDT)] included 376 cases and 400 adult controls, all of whom were mestizos (MEZ). To identify risk alleles in Mexican Indigenous individuals, SNVs were imputed from whole-exome sequencing data of 1,074 individuals. The allelic frequencies determined in MEZ and Indigenous individuals were compared with those of the continental populations from the 1,000 Genomes database phase 3. Linkage disequilibrium (LD) analysis of risk alleles was performed on all populations. Interleukin-1 receptor associated kinase 1 (IRAK1) and methyl CpG binding protein 2 (MECP2) mRNA levels were determined using real-time PCR.

Results

Case-control analysis revealed genetic association with childhood-onset SLE for all 15 SNVs (OR = 1.49-1.75; p = 0.0095 to 1.81 × 10^-4) and for the Xq28 risk haplotype (OR = 1.97, p = 4 × 10^-6). Comparing with individuals of European ancestry (0.14-0.16), the frequencies of the risk alleles were significantly higher in the MEZ individuals (0.55-0.68) and even higher in Indigenous individuals (0.57-0.83). LD analysis indicated a differential haplotype structure within the Indigenous groups, which was inherited to the MEZ population as a result of genetic admixture. Individuals homozygous for the Xq28 risk haplotype exhibited decreased levels of both MECP2A and B transcripts.

Conclusion

We found that the H1 risk haplotype differs in its conformation in the Mexican population. This difference could be attributed to positive selection within the Indigenous population, with its inheritance now having an autoimmune health impact in both the Mexican Indigenous and MEZ populations.

Persistence and expansion of cryptic endangered red wolf genomic ancestry along the American Gulf coast

Admixture and introgression play a critical role in adaptation and genetic rescue that has only recently gained a deeper appreciation. Here, we explored the geographical and genomic landscape of cryptic ancestry of the endangered red wolf that persists within the genome of a ubiquitous sister taxon, the coyote, all while the red wolf has been extinct in the wild since the early 1980s. We assessed admixture across 120,621 single nucleotiode polymorphism (SNP) loci genotyped in 293 canid genomes. We found support for increased red wolf ancestry along a west-to-east gradient across the southern United States associated with historical admixture in the past 100 years. Southwestern Louisiana and southeastern Texas, the geographical zone where the last red wolves were known prior to extinction in the wild, contained the highest and oldest levels of red wolf ancestry. Further, given the paucity of inferences based on chromosome types, we compared patterns of ancestry on the X chromosome and autosomes. We additionally aimed to explore the relationship between admixture timing and recombination rate variation to investigate gene flow events. We found that X-linked regions of low recombination rates were depleted of introgression, relative to the autosomes, consistent with the large X effect and enrichment with loci involved in maintaining reproductive isolation. Recombination rate was positively correlated with red wolf ancestry across coyote genomes, consistent with theoretical predictions. The geographical and genomic extent of cryptic red wolf ancestry can provide novel genomic resources for recovery plans targeting the conservation of the endangered red wolf.

Simultaneous inference of parental admixture proportions and admixture times from unphased local ancestry calls

Population genetic analyses of local ancestry tracts routinely assume that the ancestral admixture process is identical for both parents of an individual, an assumption that may be invalid when considering recent admixture. Here, we present Parental Admixture Proportion Inference (PAPI), a Bayesian tool for inferring the admixture proportions and admixture times for each parent of a single admixed individual. PAPI analyzes unphased local ancestry tracts and has two components: a binomial model that leverages genome-wide ancestry fractions to infer parental admixture proportions and a hidden Markov model (HMM) that infers admixture times from tract lengths. Crucially, the HMM accounts for unobserved within-ancestry recombination by approximating the pedigree crossover dynamics, enabling inference of parental admixture times. In simulations, we find that PAPI's admixture proportion estimates deviate from the truth by 0.047 on average, outperforming ANCESTOR and PedMix by 46.0% and 57.6%, respectively. Moreover, PAPI's admixture time estimates were strongly correlated with the truth (R=0.76) but have an average downward bias of 1.01 generations that is partly attributable to inaccuracies in local ancestry inference. As an illustration of its utility, we ran PAPI on African American genotypes from the PAGE study (N = 5,786) and found strong evidence of assortative mating by ancestry proportion: couples' ancestry proportions are highly correlated (R = 0.87) and are closer to each other than expected under random mating (p < 10-6). We anticipate that PAPI will be useful in studying the population dynamics of admixture and will also be of interest to individuals seeking to learn about their personal genealogies.

Reviving ghost alleles: Genetically admixed coyotes along the American Gulf Coast are critical for saving the endangered red wolf

The last known red wolves were captured in southwestern Louisiana and eastern Texas in 1980 to establish a captive breeding population. Before their extirpation, gene flow with coyotes resulted in the persistence of endangered red wolf genetic variation in local coyote populations. We assessed genomic ancestry and morphology of coyotes in southwestern Louisiana. We detected that 38 to 62% of the coyote genomes contained red wolf ancestry acquired in the past 30 years and have an admixture profile similar to that of the canids captured before the extirpation of red wolves. We further documented a positive correlation between ancestry and weight. Our findings highlight the importance of hybrids and admixed genomes as a reservoir of endangered species ancestry for innovative conservation efforts. Together, this work presents an unprecedented system that conservation can leverage to enrich the recovery program of an endangered species.

Genomic variation in baboons from central Mozambique unveils complex evolutionary relationships with other Papio species

Background

Gorongosa National Park in Mozambique hosts a large population of baboons, numbering over 200 troops. Gorongosa baboons have been tentatively identified as part of Papio ursinus on the basis of previous limited morphological analysis and a handful of mitochondrial DNA sequences. However, a recent morphological and morphometric analysis of Gorongosa baboons pinpointed the occurrence of several traits intermediate between P. ursinus and P. cynocephalus, leaving open the possibility of past and/or ongoing gene flow in the baboon population of Gorongosa National Park. In order to investigate the evolutionary history of baboons in Gorongosa, we generated high and low coverage whole genome sequence data of Gorongosa baboons and compared it to available Papio genomes.

Results

We confirmed that P. ursinus is the species closest to Gorongosa baboons. However, the Gorongosa baboon genomes share more derived alleles with P. cynocephalus than P. ursinus does, but no recent gene flow between P. ursinus and P. cynocephalus was detected when available Papio genomes were analyzed. Our results, based on the analysis of autosomal, mitochondrial and Y chromosome data, suggest complex, possibly male-biased, gene flow between Gorongosa baboons and P. cynocephalus, hinting to direct or indirect contributions from baboons belonging to the “northern” Papio clade, and signal the presence of population structure within P. ursinus.

Conclusions

The analysis of genome data generated from baboon samples collected in central Mozambique highlighted a complex set of evolutionary relationships with other baboons. Our results provided new insights in the population dynamics that have shaped baboon diversity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01999-7.

The genomic signatures of natural selection in admixed human populations

Admixture has been a pervasive phenomenon in human history, extensively shaping the patterns of population genetic diversity. There is increasing evidence to suggest that admixture can also facilitate genetic adaptation to local environments, i.e., admixed populations acquire beneficial mutations from source populations, a process that we refer to as "adaptive admixture." However, the role of adaptive admixture in human evolution and the power to detect it remain poorly characterized. Here, we use extensive computer simulations to evaluate the power of several neutrality statistics to detect natural selection in the admixed population, assuming multiple admixture scenarios. We show that statistics based on admixture proportions, F_adm and LAD, show high power to detect mutations that are beneficial in the admixed population, whereas other statistics, including iHS and F_ST, falsely detect neutral mutations that have been selected in the source populations only. By combining F_adm and LAD into a single, powerful statistic, we scanned the genomes of 15 worldwide, admixed populations for signatures of adaptive admixture. We confirm that lactase persistence and resistance to malaria have been under adaptive admixture in West Africans and in Malagasy, North Africans, and South Asians, respectively. Our approach also uncovers other cases of adaptive admixture, including APOL1 in Fulani nomads and PKN2 in East Indonesians, involved in resistance to infection and metabolism, respectively. Collectively, our study provides evidence that adaptive admixture has occurred in human populations whose genetic history is characterized by periods of isolation and spatial expansions resulting in increased gene flow.

Deciphering the genetic architecture and ethnographic distribution of IRD in three ethnic populations by whole genome sequence analysis

Patients with inherited retinal dystrophies (IRDs) were recruited from two understudied populations: Mexico and Pakistan as well as a third well-studied population of European Americans to define the genetic architecture of IRD by performing whole-genome sequencing (WGS). Whole-genome analysis was performed on 409 individuals from 108 unrelated pedigrees with IRDs. All patients underwent an ophthalmic evaluation to establish the retinal phenotype. Although the 108 pedigrees in this study had previously been examined for mutations in known IRD genes using a wide range of methodologies including targeted gene(s) or mutation(s) screening, linkage analysis and exome sequencing, the gene mutations responsible for IRD in these 108 pedigrees were not determined. WGS was performed on these pedigrees using Illumina X10 at a minimum of 30X depth. The sequence reads were mapped against hg19 followed by variant calling using GATK. The genome variants were annotated using SnpEff, PolyPhen2, and CADD score; the structural variants (SVs) were called using GenomeSTRiP and LUMPY. We identified potential causative sequence alterations in 61 pedigrees (57%), including 39 novel and 54 reported variants in IRD genes. For 57 of these pedigrees the observed genotype was consistent with the initial clinical diagnosis, the remaining 4 had the clinical diagnosis reclassified based on our findings. In seven pedigrees (12%) we observed atypical causal variants, i.e. unexpected genotype(s), including 4 pedigrees with causal variants in more than one IRD gene within all affected family members, one pedigree with intrafamilial genetic heterogeneity (different affected family members carrying causal variants in different IRD genes), one pedigree carrying a dominant causative variant present in pseudo-recessive form due to consanguinity and one pedigree with a de-novo variant in the affected family member. Combined atypical and large structural variants contributed to about 20% of cases. Among the novel mutations, 75% were detected in Mexican and 50% found in European American pedigrees and have not been reported in any other population while only 20% were detected in Pakistani pedigrees and were not previously reported. The remaining novel IRD causative variants were listed in gnomAD but were found to be very rare and population specific. Mutations in known IRD associated genes contributed to pathology in 63% Mexican, 60% Pakistani and 45% European American pedigrees analyzed. Overall, contribution of known IRD gene variants to disease pathology in these three populations was similar to that observed in other populations worldwide. This study revealed a spectrum of mutations contributing to IRD in three populations, identified a large proportion of novel potentially causative variants that are specific to the corresponding population or not reported in gnomAD and shed light on the genetic architecture of IRD in these diverse global populations.

Admixture with indigenous people helps local adaptation: admixture-enabled selection in Polynesians

Background

Homo sapiens have experienced admixture many times in the last few thousand years. To examine how admixture affects local adaptation, we investigated genomes of modern Polynesians, who are shaped through admixture between Austronesian-speaking people from Southeast Asia (Asian-related ancestors) and indigenous people in Near Oceania (Papuan-related ancestors).

Methods

In this study local ancestry was estimated across the genome in Polynesians (23 Tongan subjects) to find the candidate regions of admixture-enabled selection contributed by Papuan-related ancestors.

Results

The mean proportion of Papuan-related ancestry across the Polynesian genome was estimated as 24.6% (SD = 8.63%), and two genomic regions, the extended major histocompatibility complex (xMHC) region on chromosome 6 and the ATP-binding cassette transporter sub-family C member 11 (ABCC11) gene on chromosome 16, showed proportions of Papuan-related ancestry more than 5 SD greater than the mean (> 67.8%). The coalescent simulation under the assumption of selective neutrality suggested that such signals of Papuan-related ancestry enrichment were caused by positive selection after admixture (false discovery rate = 0.045). The ABCC11 harbors a nonsynonymous SNP, rs17822931, which affects apocrine secretory cell function. The approximate Bayesian computation indicated that, in Polynesian ancestors, a strong positive selection (s = 0.0217) acted on the ancestral allele of rs17822931 derived from Papuan-related ancestors.

Conclusions

Our results suggest that admixture with Papuan-related ancestors contributed to the rapid local adaptation of Polynesian ancestors. Considering frequent admixture events in human evolution history, the acceleration of local adaptation through admixture should be a common event in humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01900-y.

Interface of Human/Wildlife Interactions: An Example of a Bold Coyote (Canis latrans) in Atlanta, GA, USA

There is arguably no other North American species that better illustrates the complexities of the human-wildlife interface than the coyote. In this study, a melanistic coyote in metropolitan Atlanta, Georgia was exhibiting unusually bold behaviors that included encounters with humans, domestic dogs, and attempts to enter homes. After tracking this coyote (nicknamed Carmine) across a highly urbanized landscape with participatory science, including at least 80 publicly reported sightings, he was captured and relocated to a wildlife sanctuary. Genome-wide analyses revealed 92.8% coyote ancestry, 1.7% gray wolf ancestry, and 5.5% domestic dog ancestry. The dog alleles in Carmine’s genome were estimated to have been acquired by his ancestors 14–29 years ago. Despite his bold behavior, Carmine did not carry any mutations known to shape hypersociability in canines. He did, however, carry a single copy of the dominant mutation responsible for his melanistic coat color. This detailed study of Carmine dispels common assumptions about the reticent coyote personality and the origins of behavior. His unusual bold behavior created a higher level of human-coyote interaction. He now serves as a public ambassador for human-wildlife coexistence, urging the global community to reconsider mythologies about wildlife and promote coexistence with them in landscapes significantly altered by human activity in our rapidly changing world.

Genetic Ancestry Inference and Its Application for the Genetic Mapping of Human Diseases

Admixed populations arise when two or more ancestral populations interbreed. As a result of this admixture, the genome of admixed populations is defined by tracts of variable size inherited from these parental groups and has particular genetic features that provide valuable information about their demographic history. Diverse methods can be used to derive the ancestry apportionment of admixed individuals, and such inferences can be leveraged for the discovery of genetic loci associated with diseases and traits, therefore having important biomedical implications. In this review article, we summarize the most common methods of global and local genetic ancestry estimation and discuss the use of admixture mapping studies in human diseases.

Genetic admixture in the culturally unique Peranakan Chinese population in Southeast Asia

The Peranakan Chinese are culturally unique descendants of immigrants from China who settled in the Malay Archipelago ∼300-500 years ago. Today, among large communities in Southeast Asia, the Peranakans have preserved Chinese traditions with strong influence from the local indigenous Malays. Yet, whether or to what extent genetic admixture co-occurred with the cultural mixture has been a topic of ongoing debate. We performed whole-genome sequencing (WGS) on 177 Singapore (SG) Peranakans and analyzed the data jointly with WGS data of Asian and European populations. We estimated that Peranakan Chinese inherited ∼5.62% (95% confidence interval [CI]: 4.75-6.46%) Malay ancestry, much higher than that in SG Chinese (1.08%, 0.69-1.53%), southern Chinese (0.86%, 0.57-1.31%), and northern Chinese (0.25%, 0.18-0.33%). A sex-biased admixture history, in which the Malay ancestry was contributed primarily by females, was supported by X chromosomal variants, and mitochondrial (MT) and Y haplogroups. Finally, we identified an ancient admixture event shared by Peranakan Chinese and SG Chinese ∼1,612 (95% CI: 1,345-1,923) years ago, coinciding with the settlement history of Han Chinese in southern China, apart from the recent admixture event with Malays unique to Peranakan Chinese ∼190 (159-213) years ago. These findings greatly advance our understanding of the dispersal history of Chinese and their interaction with indigenous populations in Southeast Asia.

Women in Neurosurgery: First Neurosurgeon in Latin America and Current Mexican Leaders

BACKGROUND

History has taught us that Mexican culture has been largely supported by women, despite gender prejudice from the society. Neurosurgery has not been the exception. Therefore, we investigated the challenges and influence of female neurosurgeons in Mexico.

METHODS

We conducted a review of the literature and an analysis of the internal database of the Mexican Society of Neurological Surgery focusing on 3 topics: 1) the historical presence of women and gender inequality in Mexico; 2) the life and legacy of the woman who became the first neurosurgeon in Mexico and in Latin America; and 3) the participation of women in neurosurgery in the past 3 decades.

RESULTS

In Latin America, the first woman in neurosurgery was María Cristina García-Sancho, who completed her neurosurgical training in 1951. Currently, women represent 6.2% of the total members of the Mexican Society of Neurological Surgery (MSNS). This percentage is still low, although data collected in this study suggest that it might increase in the next few years because 16.7% of Board Directors of the MSNS are women, the next elected president is a female neurosurgeon, and 14.5% of neurosurgery residents are women.

CONCLUSIONS

Although a steady increase has occurred of women in neurosurgery in Mexico, there is still work to do, especially to overcome the barriers related to the old assumptions of the cultural and social roles of women.

Human genetic admixture

Throughout human history, large-scale migrations have facilitated the formation of populations with ancestry from multiple previously separated populations. This process leads to subsequent shuffling of genetic ancestry through recombination, producing variation in ancestry between populations, among individuals in a population, and along the genome within an individual. Recent methodological and empirical developments have elucidated the genomic signatures of this admixture process, bringing previously understudied admixed populations to the forefront of population and medical genetics. Under this theme, we present a collection of recent PLOS Genetics publications that exemplify recent progress in human genetic admixture studies, and we discuss potential areas for future work.

Monocyte Low-Density Lipoprotein Receptor-Related Protein 1 (LRP1) Expression Correlates with cIMT in Mexican Hypertensive Patients

Fundamento

A hipertensão arterial (HTA) representa um grande fator de risco de morbidade e mortalidade cardiovascular. Ainda não se sabe que mecanismos moleculares específicos estão associados ao desenvolvimento de hipertensão essencial.

Objetivo

Neste trabalho, analisamos a associação entre expressão mRNA de monócito LRP1, expressão de proteína LRP1, e espessura íntima-média de carótida (EIMC) de pacientes com hipertensão essencial.

Métodos

A expressão mRNA de monócito LRP1 e os níveis de proteína e EIMC foram quantificados em 200 indivíduos mexicanos, sendo 91 normotensos (NT) e 109 hipertensos (HT) A significância estatística foi definida em p < 0,05.

Resultados

O grupo de pacientes HT tinha EIMC maior altamente significativa em comparação com os pacientes NT (p = 0,002), e isso está relacionado ao aumento na expressão mRNA de LRP1 (6,54 versus. 2,87) (p = 0,002) e expressão de proteína LRP1 (17,83 versus 6,25), respectivamente (p = 0,001). Essas diferenças foram mantidas mesmo quando dividimos nossos grupos de estudo, levando em consideração apenas aqueles que apresentavam dislipidemia na expressão de mRNA (p = 0,041) e de proteínas (p < 0,001). Também se identificou que a indução de LRP1 mediada por LRP1 em monócitos em de maneira dependente de dose e tempo, com diferença significativa em NT versus HT (0,195 ± 0,09 versus 0,226 ± 0,12, p = 0,046).

Conclusão

Foi encontrado um aumento em EIMC em indivíduos com hipertensão, associada a expressões de proteína LRP1 e mRNA mais altas em monócitos, independente da presença de dislipidemia em pacientes HT. Esses resultados que a upregulation de LRP1 em monócitos de pacientes hipertensos mexicanos poderia estar envolvida na diminuição da EIMC. (Arq Bras Cardiol. 2021; 116(1):56-65)

Rare Functional Variants Associated with Antidepressant Remission in Mexican-Americans: Short title: Antidepressant remission and pharmacogenetics in Mexican-Americans

INTRODUCTION

Rare genetic functional variants can contribute to 30-40% of functional variability in genes relevant to drug action. Therefore, we investigated the role of rare functional variants in antidepressant response.

METHOD

Mexican-American individuals meeting the Diagnostic and Statistical Manual-IV criteria for major depressive disorder (MDD) participated in a prospective randomized, double-blind study with desipramine or fluoxetine. The rare variant analysis was performed using whole-exome genotyping data. Network and pathway analyses were carried out with the list of significant genes.

RESULTS

The Kernel-Based Adaptive Cluster method identified functional rare variants in 35 genes significantly associated with treatment remission (False discovery rate, FDR <0.01). Pathway analysis of these genes supports the involvement of the following gene ontology processes: olfactory/sensory transduction, regulation of response to cytokine stimulus, and meiotic cell cycleprocess.

LIMITATIONS

Our study did not have a placebo arm. We were not able to use antidepressant blood level as a covariate. Our study is based on a small sample size of only 65 Mexican-American individuals. Further studies using larger cohorts are warranted.

CONCLUSION

Our data identified several rare functional variants in antidepressant drug response in MDD patients. These have the potential to serve as genetic markers for predicting drug response.

TRIAL REGISTRATION

ClinicalTrials.gov NCT00265291.

Adaptation and co-adaptation of skin pigmentation and vitamin D genes in native Americans

We carried out an exhaustive review regarding human skin color variation and how much it may be related to vitamin D metabolism and other photosensitive molecules. We discuss evolutionary contexts that modulate this variability and hypotheses postulated to explain them; for example, a small amount of melanin in the skin facilitates vitamin D production, making it advantageous to have fair skin in an environment with little radiation incidence. In contrast, more melanin protects folate from degradation in an environment with a high incidence of radiation. Some Native American populations have a skin color at odds with what would be expected for the amount of radiation in the environment in which they live, a finding challenging the so-called "vitamin D-folate hypothesis." Since food is also a source of vitamin D, dietary habits should also be considered. Here we argue that a gene network approach provides tools to explain this phenomenon since it indicates potential alleles co-evolving in a compensatory way. We identified alleles of the vitamin D metabolism and pigmentation pathways segregated together, but in different proportions, in agriculturalists and hunter-gatherers. Finally, we highlight how an evolutionary approach can be useful to understand current topics of medical interest.

Reconstructed Lost Native American Populations from Eastern Brazil Are Shaped by Differential Jê/Tupi Ancestry

After the colonization of the Americas by Europeans and the consequent Trans-Atlantic Slave Trade, most Native American populations in eastern Brazil disappeared or went through an admixture process that configured a population composed of three main genetic components: the European, the sub-Saharan African, and the Native American. The study of the Native American genetic history is challenged by the lack of availability of genome-wide samples from Native American populations, the technical difficulties to develop ancient DNA studies, and the low proportions of the Native American component in the admixed Brazilian populations (on average 7%). We analyzed genome-wide data of 5,825 individuals from three locations of eastern Brazil: Salvador (North-East), Bambui (South-East), and Pelotas (South) and we reconstructed populations that emulate the Native American groups that were living in the 16th century around the sampling locations. This genetic reconstruction was performed after local ancestry analysis of the admixed Brazilian populations, through the rearrangement of the Native American haplotypes into reconstructed individuals with full Native American ancestry (51 reconstructed individuals in Salvador, 45 in Bambui, and 197 in Pelotas). We compared the reconstructed populations with nonadmixed Native American populations from other regions of Brazil through haplotype-based methods. Our results reveal a population structure shaped by the dichotomy of Tupi-/Jê-speaking ancestry related groups. We also show evidence of a decrease of the diversity of nonadmixed Native American groups after the European contact, in contrast with the reconstructed populations, suggesting a reservoir of the Native American genetic diversity within the admixed Brazilian population.

Ancestry-Specific Analyses Reveal Differential Demographic Histories and Opposite Selective Pressures in Modern South Asian Populations

Genetic variation in contemporary South Asian populations follows a northwest to southeast decreasing cline of shared West Eurasian ancestry. A growing body of ancient DNA evidence is being used to build increasingly more realistic models of demographic changes in the last few thousand years. Through high-quality modern genomes, these models can be tested for gene and genome level deviations. Using local ancestry deconvolution and masking, we reconstructed population-specific surrogates of the two main ancestral components for more than 500 samples from 25 South Asian populations and showed our approach to be robust via coalescent simulations.

Our f3 and f4 statistics–based estimates reveal that the reconstructed haplotypes are good proxies for the source populations that admixed in the area and point to complex interpopulation relationships within the West Eurasian component, compatible with multiple waves of arrival, as opposed to a simpler one wave scenario. Our approach also provides reliable local haplotypes for future downstream analyses. As one such example, the local ancestry deconvolution in South Asians reveals opposite selective pressures on two pigmentation genes (SLC45A2 and SLC24A5) that are common or fixed in West Eurasians, suggesting post-admixture purifying and positive selection signals, respectively.

Genetic variability of five ADRB2 polymorphisms among Mexican Amerindian ethnicities and the Mestizo population

The Mexican population is characterized by high and particular admixture, and the picture of variants associated with disease remains unclear. Here we investigated the distribution of single nucleotide polymorphisms (SNPs) in the Mexican population. We focused on two non-synonymous and three synonymous SNPs in the beta-2 adrenergic receptor gene (ADRB2), which plays key roles in energy balance regulation. These SNPs were genotyped in 2,011 Mexican Amerindians (MAs) belonging to 62 ethnic groups and in 1,980 geographically matched Mexican Mestizos (MEZs). The frequency distribution of all five ADRB2 variants significantly differed between MAs, MEZs, and other continental populations (CPs) from the 1000 Genomes database. Allele frequencies of the three synonymous SNPs rs1042717A, rs1042718A, and rs1042719C were significantly higher in Mexican individuals, particularly among MAs, compared to in the other analyzed populations (P<0.05). The non-synonymous ADRB2 Glu27 allele (rs1042714G), which is associated with several common conditions, showed the lowest frequency in MAs (0.03) compared to other populations worldwide. Among MEZs, this allele showed a frequency of 0.15, intermediate between that in MAs and in Iberians (0.43). Moreover, Glu27 was the only SNP exhibiting a geographic gradient within the MEZ population (from 0.22 to 0.11), reflecting admixed mestizo ancestry across the country. Population differentiation analysis demonstrated that Glu27 had the highest F_ST value in MAs compared with Europeans (CEU) (0.71), and the lowest between MAs and Japanese (JPT) (0.01), even lower than that observed between MAs and MEZs (0.08). This analysis demonstrated the genetic diversity among Amerindian ethnicities, with the most extreme F_ST value (0.34) found between the Nahuatls from Morelos and the Seris. This is the first study of ADRB2 genetic variants among MA ethnicities. Our findings add to our understanding of the genetic contribution to variability in disease susceptibility in admixed populations.

Inter-individual genomic heterogeneity within European population isolates

A number of studies carried out since the early ‘70s has investigated the effects of isolation on genetic variation within and among human populations in diverse geographical contexts. However, no extensive analysis has been carried out on the heterogeneity among genomes within isolated populations. This issue is worth exploring since events of recent admixture and/or subdivision could potentially disrupt the genetic homogeneity which is to be expected when isolation is prolonged and constant over time. Here, we analyze literature data relative to 87,815 autosomal single-nucleotide polymorphisms, which were obtained from a total of 28 European populations. Our results challenge the traditional paradigm of population isolates as structured as genetically (and genomically) uniform entities. In fact, focusing on the distribution of variance of intra-population diversity measures across individuals, we show that the inter-individual heterogeneity of isolated populations is at least comparable to the open ones. More in particular, three small and highly inbred isolates (Sappada, Sauris and Timau in Northeastern Italy) were found to be characterized by levels of inter-individual heterogeneity largely exceeding that of all other populations, possibly due to relatively recent events of genetic introgression. Finally, we propose a way to monitor the effects of inter-individual heterogeneity in disease-gene association studies.

Catalytically Impaired TYK2 Variants are Protective Against Childhood- and Adult-Onset Systemic Lupus Erythematosus in Mexicans

Type I interferon (IFN-I) pathway plays a central role in the systemic lupus erythematosus (SLE) pathogenesis. Recent data suggest that SLE is associated with variants in IFN-I genes, such as tyrosine kinase 2 (TYK2), which is crucial in anti-viral immunity. Here, five TYK2 single nucleotide polymorphisms (SNPs) were genotyped in 368 childhood-onset SLE Mexican patients and 516 sex-matched healthy controls. Allele frequencies were also estimated in four indigenous groups. SLE protection was associated with TYK2 risk infection variants affecting residually its catalytic domain, rs12720356 (OR = 0.308; p = 0.041) and rs34536443 (OR = 0.370; p = 0.034), but not with rs2304256, rs12720270, and rs280500. This association was replicated in a 506 adult-onset SLE patients sample (OR = 0.250; p = 0.005, and OR = 0.277; p = 0.008, respectively). The minor alleles of both associated SNPs had a lower frequency in Mestizos than in Spaniards and were absent or rare in indigenous, suggesting that the presence of these alleles in the Mexican Mestizo population was derived from the Spaniards. For the first time, we report genetic variants with a protective effect in childhood- and adult-onset SLE Mexican population. Our results suggest that the frequency of IFN-I alleles associated with SLE, may have been shaped in populations exposed to infectious diseases for long periods, and this could be an explanation why Native American ancestry is associated with a higher SLE prevalence and an earlier onset.

Genomic approaches to identify hybrids and estimate admixture times in European wildcat populations

The survival of indigenous European wildcat (Felis silvestris silvestris) populations can be locally threatened by introgressive hybridization with free-ranging domestic cats. Identifying pure wildcats and investigating the ancestry of admixed individuals becomes thus a conservation priority. We analyzed 63k cat Single Nucleotide Polymorphisms (SNPs) with multivariate, Bayesian and gene-search tools to better evaluate admixture levels between domestic and wild cats collected in Europe, timing and ancestry proportions of their hybrids and backcrosses, and track the origin (wild or domestic) of the genomic blocks carried by admixed cats, also looking for possible deviations from neutrality in their inheritance patterns. Small domestic ancestry blocks were detected in the genomes of most admixed cats, which likely originated from hybridization events occurring from 6 to 22 generations in the past. We identified about 1,900 outlier coding genes with excess of wild or domestic ancestry compared to random expectations in the admixed individuals. More than 600 outlier genes were significantly enriched for Gene Ontology (GO) categories mainly related to social behavior, functional and metabolic adaptive processes (wild-like genes), involved in cognition and neural crest development (domestic-like genes), or associated with immune system functions and lipid metabolism (parental-like genes). These kinds of genomic ancestry analyses could be reliably applied to unravel the admixture dynamics in European wildcats, as well as in other hybridizing populations, in order to design more efficient conservation plans.

Isolated postaxial polydactyly: Epidemiologic characteristics from a multicenter birth defects study

Isolated postaxial polydactyly (I-PAP), as a single defect, is a frequent malformation, characterized by an extra digit placed on the ulnar or fibular side of the limbs. Worldwide prevalence varies from as high as 225/10,000 in Nigerians to so low as 6.08/10,000 in Argentinians. Genetic-ethnic background significantly affects worldwide prevalence and type of I-PAP. Herein we describe the epidemiological characteristics of I-PAP in 697 newborns, 383 males and 314 females identified in 1,178,993 examined live births from a multicenter case-control hospital-based population study, the Mexican program of Registry and Epidemiological Surveillance of Congenital Malformations (RYVEMCE). The main characteristics analyzed included total I-PAP, stratified in Types A and B, defined as complete or incomplete extra-digit formation, respectively, sex prevalence, affected limb, laterality, parity, prematurity, delivery-type, twinning, consanguinity, and parental age. Males (6.35/10,000) are significantly more frequently affected than females (5.45/10,000), hands more than feet, left more than right limbs, and Type B (74.50%) more than A (25.50%). Prematurity and forceps use were significantly more frequent in cases than controls. An evident decreasing time-trend prevalence was present. Similar findings with other studies were males, upper and left limbs more frequently affected. Findings that were not previously reported include prematurity, forceps use, a significant decreasing time trend and an inverse ethnic prevalence for Types A (75%) and B (25%) in the Mayan population in contrast to other worldwide ethnic groups.

Glioma risk associated with extent of estimated European genetic ancestry in African Americans and Hispanics

Glioma incidence is highest in non-Hispanic Whites, and to date, glioma genome-wide association studies (GWAS) to date have only included European ancestry (EA) populations. African Americans and Hispanics in the US have varying proportions of EA, African (AA) and Native American ancestries (NAA). It is unknown if identified GWAS loci or increased EA is associated with increased glioma risk. We assessed whether EA was associated with glioma in African Americans and Hispanics. Data were obtained for 832 cases and 675 controls from the Glioma International Case-Control Study and GliomaSE Case-Control Study previously estimated to have <80% EA, or self-identify as non-White. We estimated global and local ancestry using fastStructure and RFMix, respectively, using 1,000 genomes project reference populations. Within groups with ≥40% AA (AFR_≥0.4 ), and ≥15% NAA (AMR_≥0.15 ), genome-wide association between local EA and glioma was evaluated using logistic regression conditioned on global EA for all gliomas. We identified two regions (7q21.11, p = 6.36 × 10^-4 ; 11p11.12, p = 7.0 × 10^-4 ) associated with increased EA, and one associated with decreased EA (20p12.13, p = 0.0026) in AFR_≥0.4 . In addition, we identified a peak at rs1620291 (p = 4.36 × 10^-6 ) in 7q21.3. Among AMR_≥0.15 , we found an association between increased EA in one region (12q24.21, p = 8.38 × 10^-4 ), and decreased EA in two regions (8q24.21, p = 0. 0010; 20q13.33, p = 6.36 × 10^-4 ). No other significant associations were identified. This analysis identified an association between glioma and two regions previously identified in EA populations (8q24.21, 20q13.33) and four novel regions (7q21.11, 11p11.12, 12q24.21 and 20p12.13). The identifications of novel association with EA suggest regions to target for future genetic association studies.

Contribution of MSMB promoter region gene polymorphism to early-onset prostate cancer risk in Mexican males

Sexually transmitted infections and its contribution to prostate cancer (PC) development have been relevant in different populations. MSMB gene polymorphism (rs10993994) has exhibited an association both with PC as well as the susceptibility to sexually transmitted infections. Hitherto, these conditions have been not studied in Mexico yet, neither if sexually transmitted infections could modify the MSMB and PC association. Herein, socio-demographic features, sexually transmitted infections records, the reproductive backgrounds, and the genetic characterisation were analysed in 322 incident PC cases and 628 population healthy controls from Mexico City. Whole PC, early-onset PC (PC at < 60 years old), late-onset PC (≥ 60 years old), and PC aggressiveness were used to evaluate the genetic variants contribution to PC risk using unconditional logistic regression models. Overall, none associations between the allelic variants of rs10993994 polymorphisms with whole and PC aggressiveness were found. Howbeit, the TT genotype carriers presented the highest susceptibility to develop early-onset PC (OR = 2.66; 95% CI = 1.41, 5.04; p = 0.03) than CC+CT carriers, both with codominant and recessive models. Although none association between whole PC and MSMB gene polymorphism was found, our results were reinforced by prior studies in European descendent populations, suggesting a contribution between rs10993994 and early-onset PC development.

Native ancestry is associated with optic neuritis and age of onset in hispanics with multiple sclerosis

Background and Objective

Hispanics with multiple sclerosis (MS) present younger and more often with optic neuritis (ON) as compared to Whites in the western United States. Regional differences related to Hispanic genetic admixture could be responsible. We investigated the association between global genetic ancestry and ON and age at onset of MS in Hispanics.

Methods

Data were obtained for 1033 self‐identified Hispanics with MS from four MS‐based registries from four academic institutions across the United States January 2016–April 2017. Multivariate regression models, utilizing genetic ancestry estimates for Native American (NA), African, and European ancestry, were used to assess the relationship between genetic ancestry and ON presentation and age of MS onset, defined as age at first symptom.

Results

Genetic ancestry and ON proportions varied by region where NA ancestry and ON proportions were highest among Hispanics in the southwestern United States (40% vs. 19% overall for NA and 38% vs. 25% overall for ON). A strong inverse correlation was observed between NA and European ancestry (r = −0.83). ON presentation was associated with younger age of onset (OR: 0.98; 95% CI: 0.96–0.99; P = 7.80 × 10⁻⁰³) and increased NA ancestry (OR: 2.35 for the highest versus the lowest quartile of NA ancestry; 95% CI: 1.35–4.10; P = 2.60 × 10⁻⁰³). Younger age of onset was found to be associated with a higher proportion NA (Beta: −5.58; P = 3.49 × 10⁻⁰²) and African ancestry (Beta: −10.07; P = 1.39 × 10⁻⁰³).

Interpretation

Ethnic differences associated with genetic admixture could influence clinical presentation in Hispanics with MS; underscoring the importance of considering genetic substructure in future clinical, genetic, and epigenetic studies in Hispanics.

Wolf outside, dog inside? The genomic make-up of the Czechoslovakian Wolfdog

BACKGROUND

Genomic methods can provide extraordinary tools to explore the genetic background of wild species and domestic breeds, optimize breeding practices, monitor and limit the spread of recessive diseases, and discourage illegal crossings. In this study we analysed a panel of 170k Single Nucleotide Polymorphisms with a combination of multivariate, Bayesian and outlier gene approaches to examine the genome-wide diversity and inbreeding levels in a recent wolf x dog cross-breed, the Czechoslovakian Wolfdog, which is becoming increasingly popular across Europe.

RESULTS

Pairwise FST values, multivariate and assignment procedures indicated that the Czechoslovakian Wolfdog was significantly differentiated from all the other analysed breeds and also well-distinguished from both parental populations (Carpathian wolves and German Shepherds). Coherently with the low number of founders involved in the breed selection, the individual inbreeding levels calculated from homozygosity regions were relatively high and comparable with those derived from the pedigree data. In contrast, the coefficient of relatedness between individuals estimated from the pedigrees often underestimated the identity-by-descent scores determined using genetic profiles. The timing of the admixture and the effective population size trends estimated from the LD patterns reflected the documented history of the breed. Ancestry reconstruction methods identified more than 300 genes with excess of wolf ancestry compared to random expectations, mainly related to key morphological features, and more than 2000 genes with excess of dog ancestry, playing important roles in lipid metabolism, in the regulation of circadian rhythms, in learning and memory processes, and in sociability, such as the COMT gene, which has been described as a candidate gene for the latter trait in dogs.

CONCLUSIONS

In this study we successfully applied genome-wide procedures to reconstruct the history of the Czechoslovakian Wolfdog, assess individual wolf ancestry proportions and, thanks to the availability of a well-annotated reference genome, identify possible candidate genes for wolf-like and dog-like phenotypic traits typical of this breed, including commonly inherited disorders. Moreover, through the identification of ancestry-informative markers, these genomic approaches could provide tools for forensic applications to unmask illegal crossings with wolves and uncontrolled trades of recent and undeclared wolfdog hybrids.

Modeling SNP array ascertainment with Approximate Bayesian Computation for demographic inference

Single nucleotide polymorphisms (SNPs) in commercial arrays have often been discovered in a small number of samples from selected populations. This ascertainment skews patterns of nucleotide diversity and affects population genetic inferences. We propose a demographic inference pipeline that explicitly models the SNP discovery protocol in an Approximate Bayesian Computation (ABC) framework. We simulated genomic regions according to a demographic model incorporating parameters for the divergence of three well-characterized HapMap populations and recreated the SNP distribution of a commercial array by varying the number of haploid samples and the allele frequency cut-off in the given regions. We then calculated summary statistics obtained from both the ascertained and genomic data and inferred ascertainment and demographic parameters. We implemented our pipeline to study the admixture process that gave rise to the present-day Mexican population. Our estimate of the time of admixture is closer to the historical dates than those in previous works which did not consider ascertainment bias. Although the use of whole genome sequences for demographic inference is becoming the norm, there are still underrepresented areas of the world from where only SNP array data are available. Our inference framework is applicable to those cases and will help with the demographic inference.

Investigation of short tandem repeats in major depression using whole-genome sequencing data

BACKGROUND

Major depressive disorder (MDD) is a leading contributor to global disease burden. Recent studies have shown that genetic factors play significant roles in the susceptibility to this condition; however, the underlying genetic basis currently remains largely unknown. Short tandem repeat (STR) has been proposed as an explanatory factor in the "missing heritability" of complex diseases or traits.

METHODS

We investigated STR variations from 15 MDD patients and 10 ethnically matched healthy controls based on their deep whole-genome sequencing (WGS) data. The lobSTR software was used to computationally determine STRs.

RESULTS

The results of the Mexican-American sample showed that STRs are significantly richer in healthy controls than in MDD cases on each of the 23 chromosomes (all false discovery rates, FDR P-values < 0.0062); while for the Australian of European-ancestry sample, there was no statistically significant STRs difference between MDD cases and controls.

LIMITATIONS

High quality WGS costs limited obtaining larger datasets.

CONCLUSIONS

This preliminary work is the first study that STR variations are applied to investigate MDD based on WGS data. The results on Mexican-American population may imply that within the same ancestry, targeted sequencing on a specific chromosome or region of genome would be sufficient for examining the relationship between STR and MDD. Further studies should examine larger sequencing datasets on other ethnic groups.

Properties of global- and local-ancestry adjustments in genetic association tests in admixed populations

Population substructure can lead to confounding in tests for genetic association, and failure to adjust properly can result in spurious findings. Here we address this issue of confounding by considering the impact of global ancestry (average ancestry across the genome) and local ancestry (ancestry at a specific chromosomal location) on regression parameters and relative power in ancestry-adjusted and -unadjusted models. We examine theoretical expectations under different scenarios for population substructure; applying different regression models, verifying and generalizing using simulations, and exploring the findings in real-world admixed populations. We show that admixture does not lead to confounding when the trait locus is tested directly in a single admixed population. However, if there is more complex population structure or a marker locus in linkage disequilibrium (LD) with the trait locus is tested, both global and local ancestry can be confounders. Additionally, we show the genotype parameters of adjusted and unadjusted models all provide tests for LD between the marker and trait locus, but in different contexts. The local ancestry adjusted model tests for LD in the ancestral populations, while tests using the unadjusted and the global ancestry adjusted models depend on LD in the admixed population(s), which may be enriched due to different ancestral allele frequencies. Practically, this implies that global-ancestry adjustment should be used for screening, but local-ancestry adjustment may better inform fine mapping and provide better effect estimates at trait loci.

A genomic perspective on HLA evolution

Several decades of research have convincingly shown that classical human leukocyte antigen (HLA) loci bear signatures of natural selection. Despite this conclusion, many questions remain regarding the type of selective regime acting on these loci, the time frame at which selection acts, and the functional connections between genetic variability and natural selection. In this review, we argue that genomic datasets, in particular those generated by next-generation sequencing (NGS) at the population scale, are transforming our understanding of HLA evolution. We show that genomewide data can be used to perform robust and powerful tests for selection, capable of identifying both positive and balancing selection at HLA genes. Importantly, these tests have shown that natural selection can be identified at both recent and ancient timescales. We discuss how findings from genomewide association studies impact the evolutionary study of HLA genes, and how genomic data can be used to survey adaptive change involving interaction at multiple loci. We discuss the methodological developments which are necessary to correctly interpret genomic analyses involving the HLA region. These developments include adapting the NGS analysis framework so as to deal with the highly polymorphic HLA data, as well as developing tools and theory to search for signatures of selection, quantify differentiation, and measure admixture within the HLA region. Finally, we show that high throughput analysis of molecular phenotypes for HLA genes-namely transcription levels-is now a feasible approach and can add another dimension to the study of genetic variation.

Genomic regions associated with susceptibility to Barrett's esophagus and esophageal adenocarcinoma in African Americans: The cross BETRNet admixture study

Background

Barrett’s esophagus (BE) and esophageal adenocarcinoma (EAC) are far more prevalent in European Americans than in African Americans. Hypothesizing that this racial disparity in prevalence might represent a genetic susceptibility, we used an admixture mapping approach to interrogate disease association with genomic differences between European and African ancestry.

Methods

Formalin fixed paraffin embedded samples were identified from 54 African Americans with BE or EAC through review of surgical pathology databases at participating Barrett’s Esophagus Translational Research Network (BETRNet) institutions. DNA was extracted from normal tissue, and genotyped on the Illumina OmniQuad SNP chip. Case-only admixture mapping analysis was performed on the data from both all 54 cases and also on a subset of 28 cases with high genotyping quality. Haplotype phases were inferred with Beagle 3.3.2, and local African and European ancestries were inferred with SABER plus. Disease association was tested by estimating and testing excess European ancestry and contrasting it to excess African ancestry.

Results

Both datasets, the 54 cases and the 28 cases, identified two admixture regions. An association of excess European ancestry on chromosome 11p reached a 5% genome-wide significance threshold, corresponding to -log₁₀(P) = 4.28. A second peak on chromosome 8q reached -log₁₀(P) = 2.73. The converse analysis examining excess African ancestry found no genetic regions with significant excess African ancestry associated with BE and EAC. On average, the regions on chromosomes 8q and 11p showed excess European ancestry of 15% and 20%, respectively.

Conclusions

Chromosomal regions on 11p15 and 8q22-24 are associated with excess European ancestry in African Americans with BE and EAC. Because GWAS have not reported any variants in these two regions, low frequency and/or rare disease associated variants that confer susceptibility to developing BE and EAC may be driving the observed European ancestry association evidence.

Population Stratification in Genetic Association Studies

Population stratification (PS) is a primary consideration in studies of genetic determinants of human traits. Failure to control for PS may lead to confounding, causing a study to fail for lack of significant results, or resources to be wasted following false-positive signals. Here, historical and current approaches for addressing PS when performing genetic association studies in human populations are reviewed. Methods for detecting the presence of PS, including global and local ancestry methods, are described. Also described are approaches for accounting for PS when calculating association statistics, such that measures of association are not confounded. Many traits are being examined for the first time in minority populations, which may inherently feature PS. © 2017 by John Wiley & Sons, Inc.

Investigation of copy number variation in subjects with major depression based on whole-genome sequencing data

BACKGROUND

Despite recent intensive research using genome-wide association studies, the underlying biological basis of major depressive disorder (MDD) still remains unknown. In contrast to genotyping platforms which identify specific variations, whole-genome sequencing (WGS) allows us to detect all private genetic variations within an individual. So far there have been no studies investigating copy number variations (CNVs) in subjects with MDD using WGS data.

METHODS

We obtained complete WGS paired-end reads data of 15 MDD patients and 10 ethnically matched healthy controls. We performed alignments for the sequencing reads and used GASV package to call CNVs including deletion, inversion, translocation and divergence for those subjects.

RESULTS

Our results show that, in the Mexican-American sample, deletion CNVs were significantly richer in MDD cases than healthy controls on each of 23 chromosomes. However, other types of CNVs failed to reach any significance. In the Australian sample, there was no statistically significant difference of CNVs between MDD cases and controls. Furthermore, we found that the Australian group had significantly more deletion CNVs than the Mexican-American group.

LIMITATIONS

High quality WGS costs limited obtaining larger datasets. The GASV package does not currently support duplication or insertion CNVs.

CONCLUSIONS

To our knowledge this is the first time that CNVs detected by WGS data are used to study major depression. The conclusion that deletion CNVs are significantly richer in MDD cases than healthy controls is consistent with the previous finding about recurrent depressive disorder by genome-wide association analysis of CNVs on a large genotyping microarray data.

Evidence of epigenetic admixture in the Colombian population

DNA methylation (DNAm) measured in lymphoblastoid cell lines has been repeatedly demonstrated to differ between various human populations. Due to the role that DNAm plays in controlling gene expression, these differences could significantly contribute to ethnic phenotypic differences. However, because previous studies have compared distinct ethnic groups where genetic and environmental context are confounded, their relative contribution to phenotypic differences between ethnicities remains unclear. Using DNAm assayed in whole blood and colorectal tissue of 132 admixed individuals from Colombia, we identified sites where differential DNAm levels were associated with the local ancestral genetic context. Our results are consistent with population specific DNAm being primarily driven by between population genetic differences in cis, with little environmental contribution, and with consistent effects across tissues. The findings offer new insights into a possible mechanism driving phenotypic differences among different ethnic groups, and could help explain ethnic differences in colorectal cancer incidence.

Disentangling Timing of Admixture, Patterns of Introgression, and Phenotypic Indicators in a Hybridizing Wolf Population

Hybridization is a natural or anthropogenic process that can deeply affect the genetic make-up of populations, possibly decreasing individual fitness but sometimes favoring local adaptations. The population of Italian wolves (Canis lupus), after protracted demographic declines and isolation, is currently expanding in anthropic areas, with documented cases of hybridization with stray domestic dogs. However, identifying admixture patterns in deeply introgressed populations is far from trivial. In this study, we used a panel of 170,000 SNPs analyzed with multivariate, Bayesian and local ancestry reconstruction methods to identify hybrids, estimate their ancestry proportions and timing since admixture. Moreover, we carried out preliminary genotype-phenotype association analyses to identify the genetic bases of three phenotypic traits (black coat, white claws, and spur on the hind legs) putative indicators of hybridization. Results showed no sharp subdivisions between nonadmixed wolves and hybrids, indicating that recurrent hybridization and deep introgression might have started mostly at the beginning of the population reexpansion. In hybrids, we identified a number of genomic regions with excess of ancestry in one of the parental populations, and regions with excess or resistance to introgression compared with neutral expectations. The three morphological traits showed significant genotype-phenotype associations, with a single genomic region for black coats and white claws, and with multiple genomic regions for the spur. In all cases the associated haplotypes were likely derived from dogs. In conclusion, we show that the use of multiple genome-wide ancestry reconstructions allows clarifying the admixture dynamics even in highly introgressed populations, and supports their conservation management.

Genome-Wide Analysis in Brazilians Reveals Highly Differentiated Native American Genome Regions

Despite its population, geographic size, and emerging economic importance, disproportionately little genome-scale research exists into genetic factors that predispose Brazilians to disease, or the population genetics of risk. After identification of suitable proxy populations and careful analysis of tri-continental admixture in 1,538 North-Eastern Brazilians to estimate individual ancestry and ancestral allele frequencies, we computed 400,000 genome-wide locus-specific branch length (LSBL) Fst statistics of Brazilian Amerindian ancestry compared to European and African; and a similar set of differentiation statistics for their Amerindian component compared with the closest Asian 1000 Genomes population (surprisingly, Bengalis in Bangladesh). After ranking SNPs by these statistics, we identified the top 10 highly differentiated SNPs in five genome regions in the LSBL tests of Brazilian Amerindian ancestry compared to European and African; and the top 10 SNPs in eight regions comparing their Amerindian component to the closest Asian 1000 Genomes population. We found SNPs within or proximal to the genes CIITA (rs6498115), SMC6 (rs1834619), and KLHL29 (rs2288697) were most differentiated in the Amerindian-specific branch, while SNPs in the genes ADAMTS9 (rs7631391), DOCK2 (rs77594147), SLC28A1 (rs28649017), ARHGAP5 (rs7151991), and CIITA (rs45601437) were most highly differentiated in the Asian comparison. These genes are known to influence immune function, metabolic and anthropometry traits, and embryonic development. These analyses have identified candidate genes for selection within Amerindian ancestry, and by comparison of the two analyses, those for which the differentiation may have arisen during the migration from Asia to the Americas.

National Institute of Environmental Health Sciences: 50 Years of Advancing Science and Improving Lung Health

The American Thoracic Society celebrates the 50th anniversary of the National Institute of Environmental Health Sciences (NIEHS). The NIEHS has had enormous impact through its focus on research, training, and translational science on lung health. It has been an advocate for clean air both in the United States and across the world. The cutting-edge science funded by the NIEHS has led to major discoveries that have broadened our understanding of the pathogenesis and treatment for lung disease. Importantly, the NIEHS has developed and fostered mechanisms that require cross-cutting science across the spectrum of areas of inquiry, bringing together environmental and social scientists with clinicians to bring their expertise on specific areas of investigation. The intramural program of the NIEHS nurtures cutting-edge science, and the extramural program encourages investigator-initiated research while at the same time providing broader direction through important initiatives. Under the umbrella of the NIEHS and guided by Dr. Linda Birnbaum, the director of the NIEHS, important collaborative programs, such as the Superfund Program and the National Toxicology Program, work to discover mechanisms to protect from environmental toxins. The American Thoracic Society has overlapping goals with the NIEHS, and the strategic plans of both august bodies converge to synergize on population lung health. These bonds must be tightened and highlighted as we work toward our common goals.

The PHF21B gene is associated with major depression and modulates the stress response

Major depressive disorder (MDD) affects around 350 million people worldwide; however, the underlying genetic basis remains largely unknown. In this study, we took into account that MDD is a gene-environment disorder, in which stress is a critical component, and used whole-genome screening of functional variants to investigate the 'missing heritability' in MDD. Genome-wide association studies (GWAS) using single- and multi-locus linear mixed-effect models were performed in a Los Angeles Mexican-American cohort (196 controls, 203 MDD) and in a replication European-ancestry cohort (499 controls, 473 MDD). Our analyses took into consideration the stress levels in the control populations. The Mexican-American controls, comprised primarily of recent immigrants, had high levels of stress due to acculturation issues and the European-ancestry controls with high stress levels were given higher weights in our analysis. We identified 44 common and rare functional variants associated with mild to moderate MDD in the Mexican-American cohort (genome-wide false discovery rate, FDR, <0.05), and their pathway analysis revealed that the three top overrepresented Gene Ontology (GO) processes were innate immune response, glutamate receptor signaling and detection of chemical stimulus in smell sensory perception. Rare variant analysis replicated the association of the PHF21B gene in the ethnically unrelated European-ancestry cohort. The TRPM2 gene, previously implicated in mood disorders, may also be considered replicated by our analyses. Whole-genome sequencing analyses of a subset of the cohorts revealed that European-ancestry individuals have a significantly reduced (50%) number of single nucleotide variants compared with Mexican-American individuals, and for this reason the role of rare variants may vary across populations. PHF21b variants contribute significantly to differences in the levels of expression of this gene in several brain areas, including the hippocampus. Furthermore, using an animal model of stress, we found that Phf21b hippocampal gene expression is significantly decreased in animals resilient to chronic restraint stress when compared with non-chronically stressed animals. Together, our results reveal that including stress level data enables the identification of novel rare functional variants associated with MDD.

Joint genotype- and ancestry-based genome-wide association studies in admixed populations

In genome-wide association studies (GWAS) genetic loci that influence complex traits are localized by inspecting associations between genotypes of genetic markers and the values of the trait of interest. On the other hand, admixture mapping, which is performed in case of populations consisting of a recent mix of two ancestral groups, relies on the ancestry information at each locus (locus-specific ancestry). Recently it has been proposed to jointly model genotype and locus-specific ancestry within the framework of single marker tests. Here, we extend this approach for population-based GWAS in the direction of multimarker models. A modified version of the Bayesian information criterion is developed for building a multilocus model that accounts for the differential correlation structure due to linkage disequilibrium (LD) and admixture LD. Simulation studies and a real data example illustrate the advantages of this new approach compared to single-marker analysis or modern model selection strategies based on separately analyzing genotype and ancestry data, as well as to single-marker analysis combining genotypic and ancestry information. Depending on the signal strength, our procedure automatically chooses whether genotypic or locus-specific ancestry markers are added to the model. This results in a good compromise between the power to detect causal mutations and the precision of their localization. The proposed method has been implemented in R and is available at http://www.math.uni.wroc.pl/~mbogdan/admixtures/.

Whole-genome single nucleotide variant distribution on genomic regions and its relationship to major depression

Recent advances in DNA technologies have provided unprecedented opportunities for biological and medical research. In contrast to current popular genotyping platforms which identify specific variations, whole-genome sequencing (WGS) allows for the detection of all private mutations within an individual. Major depressive disorder (MDD) is a chronic condition with enormous medical, social and economic impacts. Genetic analysis, by identifying risk variants and thereby increasing our understanding of how MDD arises, could lead to improved prevention and the development of new and more effective treatments. Here we investigated the distributions of whole-genome single nucleotide variants (SNVs) on 12 different genomic regions for 25 human subjects using the symmetrised Kullback-Leibler divergence to measure the similarity between their SNV distributions. We performed cluster analysis for MDD patients and ethnically matched healthy controls. The results showed that Mexican-American controls grouped closer; in contrast depressed Mexican-American participants grouped away from their ethnically matched controls. This implies that whole-genome SNV distribution on the genomic regions may be related to major depression.

A latent genetic subtype of major depression identified by whole-exome genotyping data in a Mexican-American cohort

Identifying data-driven subtypes of major depressive disorder (MDD) is an important topic of psychiatric research. Currently, MDD subtypes are based on clinically defined depression symptom patterns. Although a few data-driven attempts have been made to identify more homogenous subgroups within MDD, other studies have not focused on using human genetic data for MDD subtyping. Here we used a computational strategy to identify MDD subtypes based on single-nucleotide polymorphism genotyping data from MDD cases and controls using Hamming distance and cluster analysis. We examined a cohort of Mexican-American participants from Los Angeles, including MDD patients (n=203) and healthy controls (n=196). The results in cluster trees indicate that a significant latent subtype exists in the Mexican-American MDD group. The individuals in this hidden subtype have increased common genetic substrates related to major depression and they also have more anxiety and less middle insomnia, depersonalization and derealisation, and paranoid symptoms. Advances in this line of research to validate this strategy in other patient groups of different ethnicities will have the potential to eventually be translated to clinical practice, with the tantalising possibility that in the future it may be possible to refine MDD diagnosis based on genetic data.

The time and place of European admixture in Ashkenazi Jewish history

The Ashkenazi Jewish (AJ) population is important in genetics due to its high rate of Mendelian disorders. AJ appeared in Europe in the 10^th century, and their ancestry is thought to comprise European (EU) and Middle-Eastern (ME) components. However, both the time and place of admixture are subject to debate. Here, we attempt to characterize the AJ admixture history using a careful application of new and existing methods on a large AJ sample. Our main approach was based on local ancestry inference, in which we first classified each AJ genomic segment as EU or ME, and then compared allele frequencies along the EU segments to those of different EU populations. The contribution of each EU source was also estimated using GLOBETROTTER and haplotype sharing. The time of admixture was inferred based on multiple statistics, including ME segment lengths, the total EU ancestry per chromosome, and the correlation of ancestries along the chromosome. The major source of EU ancestry in AJ was found to be Southern Europe (≈60–80% of EU ancestry), with the rest being likely Eastern European. The inferred admixture time was ≈30 generations ago, but multiple lines of evidence suggest that it represents an average over two or more events, pre- and post-dating the founder event experienced by AJ in late medieval times. The time of the pre-bottleneck admixture event, which was likely Southern European, was estimated to ≈25–50 generations ago.

The Ashkenazi Jewish population has resided in Europe for much of its 1000-year existence. However, its ethnic and geographic origins are controversial, due to the scarcity of reliable historical records. Previous genetic studies have found links to Middle-Eastern and European ancestries, but the admixture history has not been studied in detail yet, partly due to technical difficulties in disentangling signals from multiple admixture events. Here, we present an in-depth analysis of the sources of European gene flow and the time of admixture events by using multiple new and existing methods and extensive simulations. Our results suggest a model of at least two events of European admixture. One event slightly pre-dated a late medieval founder event and was likely from a Southern European source. Another event post-dated the founder event and likely occurred in Eastern Europe. These results, as well as the methods introduced, will be highly valuable for geneticists and other researchers interested in Ashkenazi Jewish origins.

A novel strategy for clustering major depression individuals using whole-genome sequencing variant data

Major depressive disorder (MDD) is highly prevalent, resulting in an exceedingly high disease burden. The identification of generic risk factors could lead to advance prevention and therapeutics. Current approaches examine genotyping data to identify specific variations between cases and controls. Compared to genotyping, whole-genome sequencing (WGS) allows for the detection of private mutations. In this proof-of-concept study, we establish a conceptually novel computational approach that clusters subjects based on the entirety of their WGS. Those clusters predicted MDD diagnosis. This strategy yielded encouraging results, showing that depressed Mexican-American participants were grouped closer; in contrast ethnically-matched controls grouped away from MDD patients. This implies that within the same ancestry, the WGS data of an individual can be used to check whether this individual is within or closer to MDD subjects or to controls. We propose a novel strategy to apply WGS data to clinical medicine by facilitating diagnosis through genetic clustering. Further studies utilising our method should examine larger WGS datasets on other ethnical groups.