Emma E. Thompson

CYP3A Variation and the Evolution of Salt-Sensitivity Variants

Members of the cytochrome P450 3A subfamily catalyze the metabolism of endogenous substrates, environmental carcinogens, and clinically important exogenous compounds, such as prescription drugs and therapeutic agents. In particular, the CYP3A4 and CYP3A5 genes play an especially important role in pharmacogenetics, since they metabolize >50% of the drugs on the market. However, known genetic variants at these two loci are not sufficient to account for the observed phenotypic variability in drug response. We used a comparative genomics approach to identify conserved coding and noncoding regions at these genes and resequenced them in three ethnically diverse human populations. We show that remarkable interpopulation differences exist with regard to frequency spectrum and haplotype structure. The non-African samples are characterized by a marked excess of rare variants and the presence of a homogeneous group of long-range haplotypes at high frequency. The CYP3A5*1/*3 polymorphism, which is likely to influence salt and water retention and risk for salt-sensitive hypertension, was genotyped in >1,000 individuals from 52 worldwide population samples. The results reveal an unusual geographic pattern whereby the CYP3A5*3 frequency shows extreme variation across human populations and is significantly correlated with distance from the equator. Furthermore, we show that an unlinked variant, AGT M235T, previously implicated in hypertension and pre-eclampsia, exhibits a similar geographic distribution and is significantly correlated in frequency with CYP3A5*1/*3. Taken together, these results suggest that variants that influence salt homeostasis were the targets of a shared selective pressure that resulted from an environmental variable correlated with latitude.

Complex signatures of natural selection at the Duffy blood group locus.

The Duffy blood group locus (FY) has long been considered a likely target of natural selection, because of the extreme pattern of geographic differentiation of its three major alleles (FY*B, FY*A, and FY*O). In the present study, we resequenced the FY region in samples of Hausa from Cameroon (fixed for FY*O), Han Chinese (fixed for FY*A), Italians, and Pakistanis. Our goals were to characterize the signature of directional selection on FY*O in sub-Saharan Africa and to understand the extent to which natural selection has also played a role in the extreme geographic differentiation of the other derived allele at this locus, FY*A. The data from the FY region are compared with the patterns of variation observed at 10 unlinked, putatively neutral loci from the same populations, as well as with theoretical expectations from the neutral-equilibrium model. The FY region in the Hausa shows evidence of directional selection in two independent properties of the data (i.e., level of sequence variation and frequency spectrum), observations that are consistent with the FY*O mutation being the target. The Italian and Chinese FY data show patterns of variation that are very unusual, particularly with regard to frequency spectrum and linkage disequilibrium, but do not fit the predictions of any simple model of selection. These patterns may represent a more complex and previously unrecognized signature of positive selection.

Human spermatogenic failure purges deleterious mutation load from the autosomes and both sex chromosomes, including the gene DMRT1.

Gonadal failure, along with early pregnancy loss and perinatal death, may be an important filter that limits the propagation of harmful mutations in the human population. We hypothesized that men with spermatogenic impairment, a disease with unknown genetic architecture and a common cause of male infertility, are enriched for rare deleterious mutations compared to men with normal spermatogenesis. After assaying genomewide SNPs and CNVs in 323 Caucasian men with idiopathic spermatogenic impairment and more than 1,100 controls, we estimate that each rare autosomal deletion detected in our study multiplicatively changes a mans risk of disease by 10% (OR 1.10 [1.04–1.16], p<2×10−3), rare X-linked CNVs by 29%, (OR 1.29 [1.11–1.50], p<1×10−3), and rare Y-linked duplications by 88% (OR 1.88 [1.13–3.13], p<0.03). By contrasting the properties of our case-specific CNVs with those of CNV callsets from cases of autism, schizophrenia, bipolar disorder, and intellectual disability, we propose that the CNV burden in spermatogenic impairment is distinct from the burden of large, dominant mutations described for neurodevelopmental disorders. We identified two patients with deletions of DMRT1, a gene on chromosome 9p24.3 orthologous to the putative sex determination locus of the avian ZW chromosome system. In an independent sample of Han Chinese men, we identified 3 more DMRT1 deletions in 979 cases of idiopathic azoospermia and none in 1,734 controls, and found none in an additional 4,519 controls from public databases. The combined results indicate that DMRT1 loss-of-function mutations are a risk factor and potential genetic cause of human spermatogenic failure (frequency of 0.38% in 1306 cases and 0% in 7,754 controls, p = 6.2×10−5). Our study identifies other recurrent CNVs as potential causes of idiopathic azoospermia and generates hypotheses for directing future studies on the genetic basis of male infertility and IVF outcomes.

The ABO blood group is a trans-species polymorphism in primates

The ABO histo-blood group, the critical determinant of transfusion incompatibility, was the first genetic polymorphism discovered in humans. Remarkably, ABO antigens are also polymorphic in many other primates, with the same two amino acid changes responsible for A and B specificity in all species sequenced to date. Whether this recurrence of A and B antigens is the result of an ancient polymorphism maintained across species or due to numerous, more recent instances of convergent evolution has been debated for decades, with a current consensus in support of convergent evolution. We show instead that genetic variation data in humans and gibbons as well as in Old World monkeys are inconsistent with a model of convergent evolution and support the hypothesis of an ancient, multiallelic polymorphism of which some alleles are shared by descent among species. These results demonstrate that the A and B blood groups result from a trans-species polymorphism among distantly related species and has remained under balancing selection for tens of millions of years—to date, the only such example in hominoids and Old World monkeys outside of the major histocompatibility complex.

Novel Functional Germline Variants in the VEGF Receptor 2 Gene and Their Effect on Gene Expression and Microvessel Density in Lung Cancer

Purpose: VEGF receptor 2 (VEGFR-2) plays a crucial role in mediating angiogenic endothelial cell responses via the VEGF pathway, and angiogenesis inhibitors targeting VEGFR-2 are in clinical use. As angiogenesis is a host-driven process, functional heritable variation in KDR, the gene encoding VEGFR-2, may affect VEGFR-2 function and, ultimately, the extent of tumor angiogenesis. Experimental Design: We resequenced KDR using 24 DNAs each from healthy Caucasian, African American, and Asian groups. Nonsynonymous genetic variants were assessed for function by phosphorylation assays. Luciferase reporter gene assays were used to examine effects of variants on gene expression. KDR mRNA and protein expression and microvessel density (MVD) were measured in non–small cell lung cancer (NSCLC) tumor samples, and matching patient DNA samples were genotyped to test for associations with variants of interest. Results: KDR resequencing led to the discovery of 120 genetic variants, of which 25 had not been previously reported. Q472H had increased VEGFR-2 protein phosphorylation and associated with increased MVD in NSCLC tumor samples. −2854C and −2455A increased luciferase expression and associated with higher KDR mRNA levels in NSCLC samples. −271A reduced luciferase expression and associated with lower VEGFR-2 levels in NSCLC samples. −906C and 23408G associated with higher KDR mRNA levels in NSCLC samples. Conclusions: This study has defined KDR genetic variation in 3 populations and identified common variants that impact on tumoral KDR expression and vascularization. These findings may have important implications for understanding the molecular basis of genetic associations between KDR variation and clinical phenotypes related to VEGFR-2 function. Clin Cancer Res; 17(16); 5257–67. ©2011 AACR.

Sequencing the IL4 locus in African Americans implicates rare noncoding variants in asthma susceptibility.

BACKGROUND Common genetic variations in the IL4 gene have been associated with asthma and atopy in European and Asian populations, but not in African Americans. OBJECTIVE Because populations of African descent have increased levels of genetic variation compared with other populations, particularly with respect to low frequency or rare variants, we hypothesized that rare variants in the IL4 gene contribute to the development of asthma in African Americans. METHODS To test this hypothesis, we sequenced the IL4 locus in 72 African Americans with asthma and 70 African American controls without asthma to identify novel and rare polymorphisms in the IL4 gene that may be contributing to asthma susceptibility. RESULTS We report an excess of private noncoding single nucleotide polymorphisms (SNPs) in the subjects with asthma compared with control subjects without asthma (P = .031). Tajimas D is significantly more negative in subjects with asthma (-0.375) than controls (-0.073; P = .04), reflecting an excess of rare variants in the subjects with asthma. CONCLUSION Our findings indicate that SNPs at the IL4 locus that are potentially exclusive to African Americans are associated with susceptibility to asthma. Only 3 of the 26 private SNPs (ie, SNPs present only in the subjects with asthma or only in the controls) are tagged by single SNPs on one of the common genotyping platforms used in genome-wide association studies. We also find that most of the private SNPs cannot be reliably imputed, highlighting the importance of sequencing to identify genetic variants contributing to common diseases in African Americans.

Sequence diversity and haplotype structure at the human CYP3A cluster

The four members of the human CYP3A subfamily play important roles in the clearance of xenobiotics, hormones, and environmental compounds. Many SNPs at the CYP3A locus have been characterized, with several showing large allele frequency differences across populations. In addition to the effects of CYP3A SNPs on drug metabolism, recent studies have highlighted the potential for CYP3A variation in susceptibility to several common phenotypes, including hypertension and cancer. We previously showed that the CYP3A4 and CYP3A5 genes have a strong haplotype structure at varying frequencies across ethnic groups. Here, we extend our re-sequencing survey to the remaining CYP3A genes in the same cluster, CYP3A7 and CYP3A43. Our study identified a large number of SNPs in coding and conserved noncoding sequences, several of which are common. The combined data set allows us to investigate patterns of sequence variation and linkage disequilibrium at the entire CYP3A locus for use in future association studies.

Genome-wide association study of lung function phenotypes in a founder population.

BACKGROUND Lung function is a long-term predictor of mortality and morbidity. OBJECTIVE We sought to identify single nucleotide polymorphisms (SNPs) associated with lung function. METHODS We performed a genome-wide association study (GWAS) of FEV1, forced vital capacity (FVC), and FEV1/FVC in 1144 Hutterites aged 6 to 89 years, who are members of a founder population of European descent. We performed least absolute shrinkage and selection operation regression to select the minimum set of SNPs that best predict FEV1/FVC in the Hutterites and used the GRAIL algorithm to mine the Gene Ontology database for evidence of functional connections between genes near the predictive SNPs. RESULTS Our GWAS identified significant associations between FEV1/FVC and SNPs at the THSD4-UACA-TLE3 locus on chromosome 15q23 (P = 5.7 × 10(-8) to 3.4 × 10(-9)). Nine SNPs at or near 4 additional loci had P < 10(-5) with FEV1/FVC. Only 2 SNPs were found with P < 10(-5) for FEV1 or FVC. We found nominal levels of significance with SNPs at 9 of the 27 previously reported loci associated with lung function measures. Among a predictive set of 80 SNPs, 6 loci were identified that had a significant degree of functional connectivity (GRAIL P < .05), including 3 clusters of β-defensin genes, 2 chemokine genes (CCL18 and CXCL12), and TNFRSF13B. CONCLUSION This study identifies genome-wide significant associations and replicates results of previous GWASs. Multimarker modeling implicated for the first time common variation in genes involved in antimicrobial immunity in airway mucosa that influences lung function.

Fraction of exhaled nitric oxide values in childhood are associated with 17q11.2-q12 and 17q12-q21 variants

BACKGROUND The fraction of exhaled nitric oxide (Feno) value is a biomarker of eosinophilic airway inflammation and is associated with childhood asthma. Identification of common genetic variants associated with childhood Feno values might help to define biological mechanisms related to specific asthma phenotypes. OBJECTIVE We sought to identify the genetic variants associated with childhood Feno values and their relation with asthma. METHODS Feno values were measured in children age 5 to 15 years. In 14 genome-wide association studies (N = 8,858), we examined the associations of approximately 2.5 million single nucleotide polymorphisms (SNPs) with Feno values. Subsequently, we assessed whether significant SNPs were expression quantitative trait loci in genome-wide expression data sets of lymphoblastoid cell lines (n = 1,830) and were related to asthma in a previously published genome-wide association data set (cases, n = 10,365; control subjects: n = 16,110). RESULTS We identified 3 SNPs associated with Feno values: rs3751972 in LYR motif containing 9 (LYRM9; P = 1.97 × 10(-10)) and rs944722 in inducible nitric oxide synthase 2 (NOS2; P = 1.28 × 10(-9)), both of which are located at 17q11.2-q12, and rs8069176 near gasdermin B (GSDMB; P = 1.88 × 10(-8)) at 17q12-q21. We found a cis expression quantitative trait locus for the transcript soluble galactoside-binding lectin 9 (LGALS9) that is in linkage disequilibrium with rs944722. rs8069176 was associated with GSDMB and ORM1-like 3 (ORMDL3) expression. rs8069176 at 17q12-q21, but not rs3751972 and rs944722 at 17q11.2-q12, were associated with physician-diagnosed asthma. CONCLUSION This study identified 3 variants associated with Feno values, explaining 0.95% of the variance. Identification of functional SNPs and haplotypes in these regions might provide novel insight into the regulation of Feno values. This study highlights that both shared and distinct genetic factors affect Feno values and childhood asthma.

Shades of gray: a comparison of linkage disequilibrium between Hutterites and Europeans

Founder or isolated populations have advantages for genetic studies due to decreased genetic and environmental heterogeneity. However, whereas longer‐range linkage disequilibrium (LD) in these populations is expected to facilitate gene localization, extensive LD may actually limit the ability for gene discovery. The North American Hutterite population is one of the best characterized young founder populations and members of this isolate have been the subjects of our studies of complex traits, including fertility, asthma and cardiovascular disease, for >20 years. Here, we directly assess the patterns and extent of global LD using single nucleotide polymorphism genotypes with minor allele frequencies (MAFs) ≥5% from the Affymetrix GeneChip® Mapping 500 K array in 60 relatively unrelated Hutterites and 60 unrelated Europeans (HapMap CEU). Although LD among some marker pairs extends further in the Hutterites than in Europeans, the pattern of LD and MAF are surprisingly similar. These results indicate that (1) identifying disease genes should be no more difficult in the Hutterites than in outbred European populations, (2) the same common susceptibility alleles for complex diseases should be present in the Hutterites and outbred European populations, and (3) imputation algorithms based on HapMap CEU should be applicable to the Hutterites. Genet. Epidemiol. 34: 133–139, 2010.