Dara G. Torgerson

Meta-analysis of genome-wide association studies of asthma in ethnically diverse North American populations

Asthma is a common disease with a complex risk architecture including both genetic and environmental factors. We performed a meta-analysis of North American genome-wide association studies of asthma in 5,416 individuals with asthma (cases) including individuals of European American, African American or African Caribbean, and Latino ancestry, with replication in an additional 12,649 individuals from the same ethnic groups. We identified five susceptibility loci. Four were at previously reported loci on 17q21, near IL1RL1, TSLP and IL33, but we report for the first time, to our knowledge, that these loci are associated with asthma risk in three ethnic groups. In addition, we identified a new asthma susceptibility locus at PYHIN1, with the association being specific to individuals of African descent (P = 3.9 × 10−9). These results suggest that some asthma susceptibility loci are robust to differences in ancestry when sufficiently large samples sizes are investigated, and that ancestry-specific associations also contribute to the complex genetic architecture of asthma.

Targets of Balancing Selection in the Human Genome

Balancing selection is potentially an important biological force for maintaining advantageous genetic diversity in populations, including variation that is responsible for long-term adaptation to the environment. By serving as a means to maintain genetic variation, it may be particularly relevant to maintaining phenotypic variation in natural populations. Nevertheless, its prevalence and specific targets in the human genome remain largely unknown. We have analyzed the patterns of diversity and divergence of 13,400 genes in two human populations using an unbiased single-nucleotide polymorphism data set, a genome-wide approach, and a method that incorporates demography in neutrality tests. We identified an unbiased catalog of genes with signatures of long-term balancing selection, which includes immunity genes as well as genes encoding keratins and membrane channels; the catalog also shows enrichment in functional categories involved in cellular structure. Patterns are mostly concordant in the two populations, with a small fraction of genes showing population-specific signatures of selection. Power considerations indicate that our findings represent a subset of all targets in the genome, suggesting that although balancing selection may not have an obvious impact on a large proportion of human genes, it is a key force affecting the evolution of a number of genes in humans.

Novel Susceptibility Variants at 10p12.31-12.2 for Childhood Acute Lymphoblastic Leukemia in Ethnically Diverse Populations

BACKGROUND Acute lymphoblastic leukemia (ALL) is the most common cancer in children and the incidence of ALL varies by ethnicity. Although accumulating evidence indicates inherited predisposition to ALL, the genetic basis of ALL susceptibility in diverse ancestry has not been comprehensively examined. METHODS We performed a multiethnic genome-wide association study in 1605 children with ALL and 6661 control subjects after adjusting for population structure, with validation in three replication series of 845 case subjects and 4316 control subjects. Association was tested by two-sided logistic regression. RESULTS A novel ALL susceptibility locus at 10p12.31-12.2 (BMI1-PIP4K2A, rs7088318, P = 1.1 × 10(-11)) was identified in the genome-wide association study, with independent replication in European Americans, African Americans, and Hispanic Americans (P = .001, .009, and .04, respectively). Association was also validated at four known ALL susceptibility loci: ARID5B, IKZF1, CEBPE, and CDKN2A/2B. Associations at ARID5B, IKZF1, and BMI1-PIP4K2A variants were consistent across ethnicity, with multiple independent signals at IKZF1 and BMI1-PIP4K2A loci. The frequency of ARID5B and BMI1-PIP4K2A variants differed by ethnicity, in parallel with ethnic differences in ALL incidence. Suggestive evidence for modifying effects of age on genetic predisposition to ALL was also observed. ARID5B, IKZF1, CEBPE, and BMI1-PIP4K2A variants cumulatively conferred strong predisposition to ALL, with children carrying six to eight copies of risk alleles at a ninefold (95% confidence interval = 6.9 to 11.8) higher ALL risk relative to those carrying zero to one risk allele at these four single nucleotide polymorphisms. CONCLUSIONS These findings indicate strong associations between inherited genetic variation and ALL susceptibility in children and shed new light on ALL molecular etiology in diverse ancestry.

Darwinian and demographic forces affecting human protein coding genes

Past demographic changes can produce distortions in patterns of genetic variation that can mimic the appearance of natural selection unless the demographic effects are explicitly removed. Here we fit a detailed model of human demography that incorporates divergence, migration, admixture, and changes in population size to directly sequenced data from 13,400 protein coding genes from 20 European-American and 19 African-American individuals. Based on this demographic model, we use several new and established statistical methods for identifying genes with extreme patterns of polymorphism likely to be caused by Darwinian selection, providing the first genome-wide analysis of allele frequency distributions in humans based on directly sequenced data. The tests are based on observations of excesses of high frequency-derived alleles, excesses of low frequency-derived alleles, and excesses of differences in allele frequencies between populations. We detect numerous new genes with strong evidence of selection, including a number of genes related to psychiatric and other diseases. We also show that microRNA controlled genes evolve under extremely high constraints and are more likely to undergo negative selection than other genes. Furthermore, we show that genes involved in muscle development have been subject to positive selection during recent human history. In accordance with previous studies, we find evidence for negative selection against mutations in genes associated with Mendelian disease and positive selection acting on genes associated with several complex diseases.

Fast and accurate inference of local ancestry in Latino populations

MOTIVATION It is becoming increasingly evident that the analysis of genotype data from recently admixed populations is providing important insights into medical genetics and population history. Such analyses have been used to identify novel disease loci, to understand recombination rate variation and to detect recent selection events. The utility of such studies crucially depends on accurate and unbiased estimation of the ancestry at every genomic locus in recently admixed populations. Although various methods have been proposed and shown to be extremely accurate in two-way admixtures (e.g. African Americans), only a few approaches have been proposed and thoroughly benchmarked on multi-way admixtures (e.g. Latino populations of the Americas). RESULTS To address these challenges we introduce here methods for local ancestry inference which leverage the structure of linkage disequilibrium in the ancestral population (LAMP-LD), and incorporate the constraint of Mendelian segregation when inferring local ancestry in nuclear family trios (LAMP-HAP). Our algorithms uniquely combine hidden Markov models (HMMs) of haplotype diversity within a novel window-based framework to achieve superior accuracy as compared with published methods. Further, unlike previous methods, the structure of our HMM does not depend on the number of reference haplotypes but on a fixed constant, and it is thereby capable of utilizing large datasets while remaining highly efficient and robust to over-fitting. Through simulations and analysis of real data from 489 nuclear trio families from the mainland US, Puerto Rico and Mexico, we demonstrate that our methods achieve superior accuracy compared with published methods for local ancestry inference in Latinos.

Inherited GATA3 variants are associated with Ph-like childhood acute lymphoblastic leukemia and risk of relapse

Recent genomic profiling of childhood acute lymphoblastic leukemia (ALL) identified a high-risk subtype with an expression signature resembling that of Philadelphia chromosome–positive ALL and poor prognosis (Ph-like ALL). However, the role of inherited genetic variation in Ph-like ALL pathogenesis remains unknown. In a genome-wide association study (GWAS) of 511 ALL cases and 6,661 non-ALL controls, we identified a susceptibility locus for Ph-like ALL (GATA3, rs3824662; P = 2.17 × 10−14, odds ratio (OR) = 3.85 for Ph-like ALL versus non-ALL; P = 1.05 × 10−8, OR = 3.25 for Ph-like ALL versus non-Ph-like ALL), with independent validation. The rs3824662 risk allele was associated with somatic lesions underlying Ph-like ALL (CRLF2 rearrangement, JAK gene mutation and IKZF1 deletion) and with variation in GATA3 expression. Finally, genotype at the GATA3 SNP was also associated with early treatment response and risk of ALL relapse. Our results provide insights into interactions between inherited and somatic variants and their role in ALL pathogenesis and prognosis.

ARID5B Genetic Polymorphisms Contribute to Racial Disparities in the Incidence and Treatment Outcome of Childhood Acute Lymphoblastic Leukemia

PURPOSE Recent genome-wide screens have identified genetic variations in ARID5B associated with susceptibility to childhood acute lymphoblastic leukemia (ALL). We sought to determine the contribution of ARID5B single nucleotide polymorphisms (SNPs) to racial disparities in ALL susceptibility and treatment outcome. PATIENTS AND METHODS We compared the association between ARID5B SNP genotype and ALL susceptibility in whites (> 95% European genetic ancestry; 978 cases and 1,046 controls) versus in Hispanics (> 10% Native American ancestry; 330 cases and 541 controls). We determined the relationships between ARID5B SNP genotype and ALL relapse risk in 1,605 children treated on the Childrens Oncology Group (COG) P9904/9905 clinical trials. RESULTS Among 49 ARID5B SNPs interrogated, 10 were significantly associated with ALL susceptibility in both whites and Hispanics (P < .05), with risk alleles consistently more frequent in Hispanics than in whites. rs10821936 exhibited the most significant association in both races (P = 8.4 × 10(-20) in whites; P = 1 × 10(-6) in Hispanics), and genotype at this SNP was highly correlated with local Native American genetic ancestry (P = 1.8 × 10(-8)). Multivariate analyses in Hispanics identified an additional SNP associated with ALL susceptibility independent of rs10821936. Eight ARID5B SNPs were associated with both ALL susceptibility and relapse hazard; the alleles related to higher ALL incidence were always linked to poorer treatment outcome and were more frequent in Hispanics. CONCLUSION ARID5B polymorphisms are important determinants of childhood ALL susceptibility and treatment outcome, and they contribute to racial disparities in this disease.

Evolutionary processes acting on candidate cis-regulatory regions in humans inferred from patterns of polymorphism and divergence.

Analysis of polymorphism and divergence in the non-coding portion of the human genome yields crucial information about factors driving the evolution of gene regulation. Candidate cis-regulatory regions spanning more than 15,000 genes in 15 African Americans and 20 European Americans were re-sequenced and aligned to the chimpanzee genome in order to identify potentially functional polymorphism and to characterize and quantify departures from neutral evolution. Distortions of the site frequency spectra suggest a general pattern of selective constraint on conserved non-coding sites in the flanking regions of genes (CNCs). Moreover, there is an excess of fixed differences that cannot be explained by a Gamma model of deleterious fitness effects, suggesting the presence of positive selection on CNCs. Extensions of the McDonald-Kreitman test identified candidate cis-regulatory regions with high probabilities of positive and negative selection near many known human genes, the biological characteristics of which exhibit genome-wide trends that differ from patterns observed in protein-coding regions. Notably, there is a higher probability of positive selection in candidate cis-regulatory regions near genes expressed in the fetal brain, suggesting that a larger portion of adaptive regulatory changes has occurred in genes expressed during brain development. Overall we find that natural selection has played an important role in the evolution of candidate cis-regulatory regions throughout hominid evolution.

Genome-wide association studies of asthma indicate opposite immunopathogenesis direction from autoimmune diseases

BACKGROUND Genome-wide association studies (GWASs) of asthma have consistently implicated the ORM1-like 3 and gasdermin B (ORMDL3-GSDMB), IL33, IL-1 receptor-like 1 and IL-18 receptor 1 (IL1RL1-IL18R1), RAD50-IL13, thymic stromal lymphopoietin and WD repeat domain 36 region (TSLP-WDR36), and HLA-DR/DQ regions. OBJECTIVE A GWAS of asthma was performed in a non-Hispanic white population. METHODS A GWAS was performed in 813 Severe Asthma Research Program/Collaborative Studies on the Genetics of Asthma/Chicago Asthma Genetics Study cases and 1564 control subjects. The GWAS results were compared with those of the published GWASs of autoimmune diseases. RESULTS Multiple single nucleotide polymorphisms in the TNFAIP3 interacting protein 1 (TNIP1) gene, which interacts with TNFAIP3 and inhibits the TNF-α-induced nuclear factor κB inflammation pathway, were associated with asthma: rs1422673 (P = 3.44 × 10(-7)) and rs10036748 (P = 1.41 × 10(-6), r(2) = 0.67). rs1422673 was also associated with asthma in the published GABRIEL (P = .018) and EVE (P = 1.31 × 10(-5)) studies. The minor allele T of rs20541 in IL13 is the risk allele for asthma but the protective allele for psoriasis. The minor allele T of rs2395185 in HLA-DRA is the risk allele for asthma but the protective allele for ulcerative colitis. The minor allele A of rs2872507 in GSDMB is the protective allele for asthma but the risk allele for rheumatoid arthritis, Crohn disease, and ulcerative colitis. The T allele of rs10036748 in the TNIP1 gene is the minor protective allele for asthma but the minor or major risk allele for systemic lupus erythematosus and systemic sclerosis in non-Hispanic white or Chinese subjects, respectively. CONCLUSIONS Our study suggests that single nucleotide polymorphisms associated with both asthma and autoimmune diseases might have opposite effects on immunopathogenesis.

Recombination rates in admixed individuals identified by ancestry-based inference

Studies of recombination and how it varies depend crucially on accurate recombination maps. We propose a new approach for constructing high-resolution maps of relative recombination rates based on the observation of ancestry switch points among admixed individuals. We show the utility of this approach using simulations and by applying it to SNP genotype data from a sample of 2,565 African Americans and 299 African Caribbeans and detecting several hundred thousand recombination events. Comparison of the inferred map with high-resolution maps from non-admixed populations provides evidence of fine-scale differentiation in recombination rates between populations. Overall, the admixed map is well predicted by the average proportion of admixture and the recombination rate estimates from the source populations. The exceptions to this are in areas surrounding known large chromosomal structural variants, specifically inversions. These results suggest that outside of structurally variable regions, admixture does not substantially disrupt the factors controlling recombination rates in humans.