Admixture Mapping of 15,280 African Americans Identifies Obesity Susceptibility Loci on Chromosomes 5 and X

The prevalence of obesity (body mass index (BMI) ≥30 kg/m²) is higher in African Americans than in European Americans, even after adjustment for socioeconomic factors, suggesting that genetic factors may explain some of the difference. To identify genetic loci influencing BMI, we carried out a pooled analysis of genome-wide admixture mapping scans in 15,280 African Americans from 14 epidemiologic studies. Samples were genotyped at a median of 1,411 ancestry-informative markers. After adjusting for age, sex, and study, BMI was analyzed both as a dichotomized (top 20% versus bottom 20%) and a continuous trait. We found that a higher percentage of European ancestry was significantly correlated with lower BMI (ρ=−0.042, P=1.6×10⁻⁷). In the dichotomized analysis, we detected two loci on chromosome X as associated with increased African ancestry: the first at Xq25 (locus-specific LOD=5.94; genome-wide score=3.22; case-control Z=−3.94); and the second at Xq13.1 (locus-specific LOD=2.22; case-control Z=−4.62). Quantitative analysis identified a third locus at 5q13.3 where higher BMI was highly significantly associated with greater European ancestry (locus-specific LOD=6.27; genome-wide score=3.46). Further mapping studies with dense sets of markers will be necessary to identify the alleles in these regions of chromosomes X and 5 that may be associated with variation in BMI.     

NRXN3 Is a Novel Locus for Waist Circumference: A Genome-Wide Association Study from the CHARGE Consortium

Central abdominal fat is a strong risk factor for diabetes and cardiovascular disease. To identify common variants influencing central abdominal fat, we conducted a two-stage genome-wide association analysis for waist circumference (WC). In total, three loci reached genome-wide significance. In stage 1, 31,373 individuals of Caucasian descent from eight cohort studies confirmed the role of FTO and MC4R and identified one novel locus associated with WC in the neurexin 3 gene [NRXN3 (rs10146997, p=6.4×10⁻⁷)]. The association with NRXN3 was confirmed in stage 2 by combining stage 1 results with those from 38,641 participants in the GIANT consortium (p=0.009 in GIANT only, p=5.3×10⁻⁸ for combined analysis, n=70,014). Mean WC increase per copy of the G allele was 0.0498 z-score units (0.65 cm). This SNP was also associated with body mass index (BMI) [p=7.4×10⁻⁶, 0.024 z-score units (0.10 kg/m²) per copy of the G allele] and the risk of obesity (odds ratio 1.13, 95% CI 1.07–1.19; p=3.2×10⁻⁵ per copy of the G allele). The NRXN3 gene has been previously implicated in addiction and reward behavior, lending further evidence that common forms of obesity may be a central nervous system-mediated disorder. Our findings establish that common variants in NRXN3 are associated with WC, BMI, and obesity.     

Meta-Analysis of Genome-Wide Scans for Human Adult Stature Identifies Novel Loci and Associations with Measures of Skeletal Frame Size

Recent genome-wide (GW) scans have identified several independent loci affecting human stature, but their contribution through the different skeletal components of height is still poorly understood. We carried out a genome-wide scan in 12,611 participants, followed by replication in an additional 7,187 individuals, and identified 17 genomic regions with GW-significant association with height. Of these, two are entirely novel (rs11809207 in CATSPER4, combined P-value=6.1×10⁻⁸ and rs910316 in TMED10, P-value=1.4×10⁻⁷) and two had previously been described with weak statistical support (rs10472828 in NPR3, P-value=3×10⁻⁷ and rs849141 in JAZF1, P-value=3.2×10⁻¹¹). One locus (rs1182188 at GNA12) identifies the first height eQTL. We also assessed the contribution of height loci to the upper- (trunk) and lower-body (hip axis and femur) skeletal components of height. We find evidence for several loci associated with trunk length (including rs6570507 in GPR126, P-value=4×10⁻⁵ and rs6817306 in LCORL, P-value=4×10⁻⁴), hip axis length (including rs6830062 at LCORL, P-value=4.8×10⁻⁴ and rs4911494 at UQCC, P-value=1.9×10⁻⁴), and femur length (including rs710841 at PRKG2, P-value=2.4×10⁻⁵ and rs10946808 at HIST1H1D, P-value=6.4×10⁻⁶). Finally, we used conditional analyses to explore a possible differential contribution of the height loci to these different skeletal size measurements. In addition to validating four novel loci controlling adult stature, our study represents the first effort to assess the contribution of genetic loci to three skeletal components of height. Further statistical tests in larger numbers of individuals will be required to verify if the height loci affect height preferentially through these subcomponents of height.     

Reduced Neutrophil Count in People of African Descent Is Due To a Regulatory Variant in the Duffy Antigen Receptor for Chemokines Gene

Persistently low white blood cell count (WBC) and neutrophil count is a well-described phenomenon in persons of African ancestry, whose etiology remains unknown. We recently used admixture mapping to identify an approximately 1-megabase region on chromosome 1, where ancestry status (African or European) almost entirely accounted for the difference in WBC between African Americans and European Americans. To identify the specific genetic change responsible for this association, we analyzed genotype and phenotype data from 6,005 African Americans from the Jackson Heart Study (JHS), the Health, Aging and Body Composition (Health ABC) Study, and the Atherosclerosis Risk in Communities (ARIC) Study. We demonstrate that the causal variant must be at least 91% different in frequency between West Africans and European Americans. An excellent candidate is the Duffy Null polymorphism (SNP rs2814778 at chromosome 1q23.2), which is the only polymorphism in the region known to be so differentiated in frequency and is already known to protect against Plasmodium vivax malaria. We confirm that rs2814778 is predictive of WBC and neutrophil count in African Americans above beyond the previously described admixture association (P=3.8×10⁻⁵), establishing a novel phenotype for this genetic variant.     

A genome-wide association study identifies novel alleles associated with hair color and skin pigmentation

We conducted a multi-stage genome-wide association study of natural hair color in more than 10,000 men and women of European ancestry from the United States and Australia. An initial analysis of 528,173 single nucleotide polymorphisms (SNPs) genotyped on 2,287 women identified IRF4 and SLC24A4 as loci highly associated with hair color, along with three other regions encompassing known pigmentation genes. We confirmed these associations in 7,028 individuals from three additional studies. Across these four studies, SLC24A4 rs12896399 and IRF4 rs12203592 showed strong associations with hair color, with p=6.0×10⁻⁶² and p=7.46×10⁻¹²⁷, respectively. The IRF4 SNP was also associated with skin color (p=6.2×10⁻¹⁴), eye color (p=6.1×10⁻¹³), and skin tanning response to sunlight (p=3.9×10⁻⁸⁹). A multivariable analysis pooling data from the initial GWAS and an additional 1,440 individuals suggested that the association between rs12203592 and hair color was independent of rs1540771, a SNP between the IRF4 and EXOC2 genes previously found to be associated with hair color. After adjustment for rs12203592, the association between rs1540771 and hair color was not significant (p=0.52). One variant in the MATP gene was associated with hair color. A variant in the HERC2 gene upstream of the OCA2 gene showed the strongest and independent association with hair color compared with other SNPs in this region, including three previously reported SNPs. The signals detected in a region around the MC1R gene were explained by MC1R red hair color alleles. Our results suggest that the IRF4 and SLC24A4 loci are associated with human hair color and skin pigmentation.     

Meta-Analysis of 28,141 Individuals Identifies Common Variants within Five New Loci That Influence Uric Acid Concentrations

Elevated serum uric acid levels cause gout and are a risk factor for cardiovascular disease and diabetes. To investigate the polygenetic basis of serum uric acid levels, we conducted a meta-analysis of genome-wide association scans from 14 studies totalling 28,141 participants of European descent, resulting in identification of 954 SNPs distributed across nine loci that exceeded the threshold of genome-wide significance, five of which are novel. Overall, the common variants associated with serum uric acid levels fall in the following nine regions: SLC2A9 (p=5.2×10⁻²⁰¹), ABCG2 (p=3.1×10⁻²⁶), SLC17A1 (p=3.0×10⁻¹⁴), SLC22A11 (p=6.7×10⁻¹⁴), SLC22A12 (p=2.0×10⁻⁹), SLC16A9 (p=1.1×10⁻⁸), GCKR (p=1.4×10⁻⁹), LRRC16A (p=8.5×10⁻⁹), and near PDZK1 (p=2.7×10⁻⁹). Identified variants were analyzed for gender differences. We found that the minor allele for rs734553 in SLC2A9 has greater influence in lowering uric acid levels in women and the minor allele of rs2231142 in ABCG2 elevates uric acid levels more strongly in men compared to women. To further characterize the identified variants, we analyzed their association with a panel of metabolites. rs12356193 within SLC16A9 was associated with DL-carnitine (p=4.0×10⁻²⁶) and propionyl-L-carnitine (p=5.0×10⁻⁸) concentrations, which in turn were associated with serum UA levels (p=1.4×10⁻⁵⁷ and p=8.1×10⁻⁵⁴, respectively), forming a triangle between SNP, metabolites, and UA levels. Taken together, these associations highlight additional pathways that are important in the regulation of serum uric acid levels and point toward novel potential targets for pharmacological intervention to prevent or treat hyperuricemia. In addition, these findings strongly support the hypothesis that transport proteins are key in regulating serum uric acid levels.     

Specificity of the STAT4 genetic association for severe disease manifestations of systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a genetically complex disease with heterogeneous clinical manifestations. A polymorphism in the STAT4 gene has recently been established as a risk factor for SLE, but the relationship with specific SLE subphenotypes has not been studied. We studied 137 SNPs in the STAT4 region genotyped in 4 independent SLE case series (total n=1398) and 2560 healthy controls, along with clinical data for the cases. Using conditional testing, we confirmed the most significant STAT4 haplotype for SLE risk. We then studied a SNP marking this haplotype for association with specific SLE subphenotypes, including autoantibody production, nephritis, arthritis, mucocutaneous manifestations, and age at diagnosis. To prevent possible type-I errors from population stratification, we reanalyzed the data using a subset of subjects determined to be most homogeneous based on principal components analysis of genome-wide data. We confirmed that four SNPs in very high LD (r²=0.94 to 0.99) were most strongly associated with SLE, and there was no compelling evidence for additional SLE risk loci in the STAT4 region. SNP rs7574865 marking this haplotype had a minor allele frequency (MAF)=31.1% in SLE cases compared with 22.5% in controls (OR=1.56, p=10⁻¹⁶). This SNP was more strongly associated with SLE characterized by double-stranded DNA autoantibodies (MAF=35.1%, OR=1.86, p<10⁻¹⁹), nephritis (MAF=34.3%, OR=1.80, p<10⁻¹¹), and age at diagnosis<30 years (MAF=33.8%, OR=1.77, p<10⁻¹³). An association with severe nephritis was even more striking (MAF=39.2%, OR=2.35, p<10⁻⁴ in the homogeneous subset of subjects). In contrast, STAT4 was less strongly associated with oral ulcers, a manifestation associated with milder disease. We conclude that this common polymorphism of STAT4 contributes to the phenotypic heterogeneity of SLE, predisposing specifically to more severe disease.     

Heterogeneity of breast cancer associations with five susceptibility loci by clinical and pathological characteristics

A three-stage genome-wide association study recently identified single nucleotide polymorphisms (SNPs) in five loci (fibroblast growth receptor 2 (FGFR2), trinucleotide repeat containing 9 (TNRC9), mitogen-activated protein kinase 3 K1 (MAP3K1), 8q24, and lymphocyte-specific protein 1 (LSP1)) associated with breast cancer risk. We investigated whether the associations between these SNPs and breast cancer risk varied by clinically important tumor characteristics in up to 23,039 invasive breast cancer cases and 26,273 controls from 20 studies. We also evaluated their influence on overall survival in 13,527 cases from 13 studies. All participants were of European or Asian origin. rs2981582 in FGFR2 was more strongly related to ER-positive (per-allele OR (95%CI)=1.31 (1.27–1.36)) than ER-negative (1.08 (1.03–1.14)) disease (P for heterogeneity=10⁻¹³). This SNP was also more strongly related to PR-positive, low grade and node positive tumors (P=10⁻⁵, 10⁻⁸, 0.013, respectively). The association for rs13281615 in 8q24 was stronger for ER-positive, PR-positive, and low grade tumors (P=0.001, 0.011 and 10⁻⁴, respectively). The differences in the associations between SNPs in FGFR2 and 8q24 and risk by ER and grade remained significant after permutation adjustment for multiple comparisons and after adjustment for other tumor characteristics. Three SNPs (rs2981582, rs3803662, and rs889312) showed weak but significant associations with ER-negative disease, the strongest association being for rs3803662 in TNRC9 (1.14 (1.09–1.21)). rs13281615 in 8q24 was associated with an improvement in survival after diagnosis (per-allele HR=0.90 (0.83–0.97). The association was attenuated and non-significant after adjusting for known prognostic factors. Our findings show that common genetic variants influence the pathological subtype of breast cancer and provide further support for the hypothesis that ER-positive and ER-negative disease are biologically distinct. Understanding the etiologic heterogeneity of breast cancer may ultimately result in improvements in prevention, early detection, and treatment.     

Genome-Wide Association Study of Plasma Polyunsaturated Fatty Acids in the InCHIANTI Study

Polyunsaturated fatty acids (PUFA) have a role in many physiological processes, including energy production, modulation of inflammation, and maintenance of cell membrane integrity. High plasma PUFA concentrations have been shown to have beneficial effects on cardiovascular disease and mortality. To identify genetic contributors of plasma PUFA concentrations, we conducted a genome-wide association study of plasma levels of six omega-3 and omega-6 fatty acids in 1,075 participants in the InCHIANTI study on aging. The strongest evidence for association was observed in a region of chromosome 11 that encodes three fatty acid desaturases (FADS1, FADS2, FADS3). The SNP with the most significant association was rs174537 near FADS1 in the analysis of arachidonic acid (AA; p=5.95×10⁻⁴⁶). Minor allele homozygotes had lower AA compared to the major allele homozygotes and rs174537 accounted for 18.6% of the additive variance in AA concentrations. This SNP was also associated with levels of eicosadienoic acid (EDA; p=6.78×10⁻⁹) and eicosapentanoic acid (EPA; p=1.07×10⁻¹⁴). Participants carrying the allele associated with higher AA, EDA, and EPA also had higher low-density lipoprotein (LDL-C) and total cholesterol levels. Outside the FADS gene cluster, the strongest region of association mapped to chromosome 6 in the region encoding an elongase of very long fatty acids 2 (ELOVL2). In this region, association was observed with EPA (rs953413; p=1.1×10⁻⁶). The effects of rs174537 were confirmed in an independent sample of 1,076 subjects participating in the GOLDN study. The ELOVL2 SNP was associated with docosapentanoic and DHA but not with EPA in GOLDN. These findings show that polymorphisms of genes encoding enzymes in the metabolism of PUFA contribute to plasma concentrations of fatty acids.     

Multiple Loci within the Major Histocompatibility Complex Confer Risk of Psoriasis

Psoriasis is a common inflammatory skin disease characterized by thickened scaly red plaques. Previously we have performed a genome-wide association study (GWAS) on psoriasis with 1,359 cases and 1,400 controls, which were genotyped for 447,249 SNPs. The most significant finding was for SNP rs12191877, which is in tight linkage disequilibrium with HLA-Cw*0602, the consensus risk allele for psoriasis. However, it is not known whether there are other psoriasis loci within the MHC in addition to HLA-C. In the present study, we searched for additional susceptibility loci within the human leukocyte antigen (HLA) region through in-depth analyses of the GWAS data; then, we followed up our findings in an independent Han Chinese 1,139 psoriasis cases and 1,132 controls. Using the phased CEPH dataset as a reference, we imputed the HLA-Cw*0602 in all samples with high accuracy. The association of the imputed HLA-Cw*0602 dosage with disease was much stronger than that of the most significantly associated SNP, rs12191877. Adjusting for HLA-Cw*0602, there were two remaining association signals: one demonstrated by rs2073048 (p=2×10⁻⁶, OR=0.66), located within c6orf10, a potential downstream effecter of TNF-alpha, and one indicated by rs13437088 (p=9×10⁻⁶, OR=1.3), located 30 kb centromeric of HLA-B and 16 kb telomeric of MICA. When HLA-Cw*0602, rs2073048, and rs13437088 were all included in a logistic regression model, each of them was significantly associated with disease (p=3×10⁻⁴⁷, 6×10⁻⁸, and 3×10⁻⁷, respectively). Both putative loci were also significantly associated in the Han Chinese samples after controlling for the imputed HLA-Cw*0602. A detailed analysis of HLA-B in both populations demonstrated that HLA-B*57 was associated with an increased risk of psoriasis and HLA-B*40 a decreased risk, independently of HLA-Cw*0602 and the C6orf10 locus, suggesting the potential pathogenic involvement of HLA-B. These results demonstrate that there are at least two additional loci within the MHC conferring risk of psoriasis.     

A candidate gene approach identifies the CHRNA5-A3-B4 region as a risk factor for age-dependent nicotine addiction

People who begin daily smoking at an early age are at greater risk of long-term nicotine addiction. We tested the hypothesis that associations between nicotinic acetylcholine receptor (nAChR) genetic variants and nicotine dependence assessed in adulthood will be stronger among smokers who began daily nicotine exposure during adolescence. We compared nicotine addiction—measured by the Fagerstrom Test of Nicotine Dependence—in three cohorts of long-term smokers recruited in Utah, Wisconsin, and by the NHLBI Lung Health Study, using a candidate-gene approach with the neuronal nAChR subunit genes. This SNP panel included common coding variants and haplotypes detected in eight α and three β nAChR subunit genes found in European American populations. In the 2,827 long-term smokers examined, common susceptibility and protective haplotypes at the CHRNA5-A3-B4 locus were associated with nicotine dependence severity (p=2.0×10⁻⁵; odds ratio=1.82; 95% confidence interval 1.39–2.39) in subjects who began daily smoking at or before the age of 16, an exposure period that results in a more severe form of adult nicotine dependence. A substantial shift in susceptibility versus protective diplotype frequency (AA versus BC=17%, AA versus CC=27%) was observed in the group that began smoking by age 16. This genetic effect was not observed in subjects who began daily nicotine use after the age of 16. These results establish a strong mechanistic link among early nicotine exposure, common CHRNA5-A3-B4 haplotypes, and adult nicotine addiction in three independent populations of European origins. The identification of an age-dependent susceptibility haplotype reinforces the importance of preventing early exposure to tobacco through public health policies.     

Genetic Differences between the Determinants of Lipid Profile Phenotypes in African and European Americans: The Jackson Heart Study

Genome-wide association analysis in populations of European descent has recently found more than a hundred genetic variants affecting risk for common disease. An open question, however, is how relevant the variants discovered in Europeans are to other populations. To address this problem for cardiovascular phenotypes, we studied a cohort of 4,464 African Americans from the Jackson Heart Study (JHS), in whom we genotyped both a panel of 12 recently discovered genetic variants known to predict lipid profile levels in Europeans and a panel of up to 1,447 ancestry informative markers allowing us to determine the African ancestry proportion of each individual at each position in the genome. Focusing on lipid profiles—HDL-cholesterol (HDL-C), LDL-cholesterol (LDL-C), and triglycerides (TG)—we identified the lipoprotein lipase (LPL) locus as harboring variants that account for interethnic variation in HDL-C and TG. In particular, we identified a novel common variant within LPL that is strongly associated with TG (p=2.7×10⁻⁶) and explains nearly 1% of the variability in this phenotype, the most of any variant in African Americans to date. Strikingly, the extensively studied “gain-of-function” S447X mutation at LPL, which has been hypothesized to be the major determinant of the LPL-TG genetic association and is in trials for human gene therapy, has a significantly diminished strength of biological effect when it is found on a background of African rather than European ancestry. These results suggest that there are other, yet undiscovered variants at the locus that are truly causal (and are in linkage disequilibrium with S447X) or that work synergistically with S447X to modulate TG levels. Finally, we find systematically lower effect sizes for the 12 risk variants discovered in European populations on the African local ancestry background in JHS, highlighting the need for caution in the use of genetic variants for risk assessment across different populations.     

Neuropeptide Y Gene Polymorphisms Confer Risk of Early-Onset Atherosclerosis

Neuropeptide Y (NPY) is a strong candidate gene for coronary artery disease (CAD). We have previously identified genetic linkage to familial CAD in the genomic region of NPY. We performed follow-up genetic, biostatistical, and functional analysis of NPY in early-onset CAD. In familial CAD (GENECARD, N=420 families), we found increased microsatellite linkage to chromosome 7p14 (OSA LOD=4.2, p=0.004) in 97 earliest age-of-onset families. Tagged NPY SNPs demonstrated linkage to CAD of a 6-SNP block (LOD=1.58–2.72), family-based association of this block with CAD (p=0.02), and stronger linkage to CAD in the earliest age-of-onset families. Association of this 6-SNP block with CAD was validated in: (a) 556 non-familial early-onset CAD cases and 256 controls (OR 1.46–1.65, p=0.01–0.05), showing stronger association in youngest cases (OR 1.84–2.20, p=0.0004–0.09); and (b) GENECARD probands versus non-familial controls (OR 1.79–2.06, p=0.003–0.02). A promoter SNP (rs16147) within this 6-SNP block was associated with higher plasma NPY levels (p=0.04). To assess a causal role of NPY in atherosclerosis, we applied the NPY1-receptor–antagonist BIBP-3226 adventitially to endothelium-denuded carotid arteries of apolipoprotein E-deficient mice; treatment reduced atherosclerotic neointimal area by 50% (p=0.03). Thus, NPY variants associate with atherosclerosis in two independent datasets (with strong age-of-onset effects) and show allele-specific expression with NPY levels, while NPY receptor antagonism reduces atherosclerosis in mice. We conclude that NPY contributes to atherosclerosis pathogenesis.     

Genome-Wide Association Scan Meta-Analysis Identifies Three Loci Influencing Adiposity and Fat Distribution

To identify genetic loci influencing central obesity and fat distribution, we performed a meta-analysis of 16 genome-wide association studies (GWAS, N=38,580) informative for adult waist circumference (WC) and waist–hip ratio (WHR). We selected 26 SNPs for follow-up, for which the evidence of association with measures of central adiposity (WC and/or WHR) was strong and disproportionate to that for overall adiposity or height. Follow-up studies in a maximum of 70,689 individuals identified two loci strongly associated with measures of central adiposity; these map near TFAP2B (WC, P=1.9×10⁻¹¹) and MSRA (WC, P=8.9×10⁻⁹). A third locus, near LYPLAL1, was associated with WHR in women only (P=2.6×10⁻⁸). The variants near TFAP2B appear to influence central adiposity through an effect on overall obesity/fat-mass, whereas LYPLAL1 displays a strong female-only association with fat distribution. By focusing on anthropometric measures of central obesity and fat distribution, we have identified three loci implicated in the regulation of human adiposity.     

Identification and replication of three novel myopia common susceptibility gene loci on chromosome 3q26 using linkage and linkage disequilibrium mapping

Refractive error is a highly heritable quantitative trait responsible for considerable morbidity. Following an initial genome-wide linkage study using microsatellite markers, we confirmed evidence for linkage to chromosome 3q26 and then conducted fine-scale association mapping using high-resolution linkage disequilibrium unit (LDU) maps. We used a preliminary discovery marker set across the 30-Mb region with an average SNP density of 1 SNP/15 kb (Map 1). Map 1 was divided into 51 LDU windows and additional SNPs were genotyped for six regions (Map 2) that showed preliminary evidence of multi-marker association using composite likelihood. A total of 575 cases and controls selected from the tails of the trait distribution were genotyped for the discovery sample. Malecot model estimates indicate three loci with putative common functional variants centred on MFN1 (180,566 kb; 95% confidence interval 180,505–180, 655 kb), approximately 156 kb upstream from alternate-splicing SOX2OT (182,595 kb; 95% CI 182,533–182,688 kb) and PSARL (184,386 kb; 95% CI 184,356–184,411 kb), with the loci showing modest to strong evidence of association for the Map 2 discovery samples (p<10⁻⁷, p<10⁻¹⁰, and p=0.01, respectively). Using an unselected independent sample of 1,430 individuals, results replicated for the MFN1 (p=0.006), SOX2OT (p=0.0002), and PSARL (p=0.0005) gene regions. MFN1 and PSARL both interact with OPA1 to regulate mitochondrial fusion and the inhibition of mitochondrial-led apoptosis, respectively. That two mitochondrial regulatory processes in the retina are implicated in the aetiology of myopia is surprising and is likely to provide novel insight into the molecular genetic basis of common myopia.     

Genome-Wide Analyses of Exonic Copy Number Variants in a Family-Based Study Point to Novel Autism Susceptibility Genes

The genetics underlying the autism spectrum disorders (ASDs) is complex and remains poorly understood. Previous work has demonstrated an important role for structural variation in a subset of cases, but has lacked the resolution necessary to move beyond detection of large regions of potential interest to identification of individual genes. To pinpoint genes likely to contribute to ASD etiology, we performed high density genotyping in 912 multiplex families from the Autism Genetics Resource Exchange (AGRE) collection and contrasted results to those obtained for 1,488 healthy controls. Through prioritization of exonic deletions (eDels), exonic duplications (eDups), and whole gene duplication events (gDups), we identified more than 150 loci harboring rare variants in multiple unrelated probands, but no controls. Importantly, 27 of these were confirmed on examination of an independent replication cohort comprised of 859 cases and an additional 1,051 controls. Rare variants at known loci, including exonic deletions at NRXN1 and whole gene duplications encompassing UBE3A and several other genes in the 15q11–q13 region, were observed in the course of these analyses. Strong support was likewise observed for previously unreported genes such as BZRAP1, an adaptor molecule known to regulate synaptic transmission, with eDels or eDups observed in twelve unrelated cases but no controls (p=2.3×10⁻⁵). Less is known about MDGA2, likewise observed to be case-specific (p=1.3×10⁻⁴). But, it is notable that the encoded protein shows an unexpectedly high similarity to Contactin 4 (BLAST E-value=3×10⁻³⁹), which has also been linked to disease. That hundreds of distinct rare variants were each seen only once further highlights complexity in the ASDs and points to the continued need for larger cohorts.     

In Vivo Protein Interactions within the Escherichia coli DNA Polymerase III Core

The mechanisms that control the fidelity of DNA replication are being investigated by a number of approaches, including detailed kinetic and structural studies. Important tools in these studies are mutant versions of DNA polymerases that affect the fidelity of DNA replication. It has been suggested that proper interactions within the core of DNA polymerase III (Pol III) of Escherichia coli could be essential for maintaining the optimal fidelity of DNA replication (H. Maki and A. Kornberg, Proc. Natl. Acad. Sci. USA 84:4389–4392, 1987). We have been particularly interested in elucidating the physiological role of the interactions between the DnaE (α subunit [possessing DNA polymerase activity]) and DnaQ ( subunit [possessing 3′→5′ exonucleolytic proofreading activity]) proteins. In an attempt to achieve this goal, we have used the Saccharomyces cerevisiae two-hybrid system to analyze specific in vivo protein interactions. In this report, we demonstrate interactions between the DnaE and DnaQ proteins and between the DnaQ and HolE (θ subunit) proteins. We also tested the interactions of the wild-type DnaE and HolE proteins with three well-known mutant forms of DnaQ (MutD5, DnaQ926, and DnaQ49), each of which leads to a strong mutator phenotype. Our results show that the mutD5 and dnaQ926 mutations do not affect the subunit-α subunit and subunit-θ subunit interactions. However, the dnaQ49 mutation greatly reduces the strength of interaction of the subunit with both the α and the θ subunits. Thus, the mutator phenotype of dnaQ49 may be the result of an altered conformation of the protein, which leads to altered interactions within the Pol III core.     

Exocyclic amino groups of flanking guanines govern sequence-dependent adduct conformations and local structural distortions for minor groove-aligned benzo[a]pyrenyl-guanine lesions in a GG mutation hotspot context

The environmental carcinogen benzo[a]pyrene (BP) is metabolized to reactive diol epoxides that bind to cellular DNA by predominantly forming N²-guanine adducts (G*). Mutation hotspots for these adducts are frequently found in 5′-···GG··· dinucleotide sequences, but their origins are poorly understood. Here we used high resolution NMR and molecular dynamics simulations to investigate differences in G* adduct conformations in 5′-···CG*GC··· and 5′-···CGG*C··· sequence contexts in otherwise identical 12-mer duplexes. The BP rings are positioned 5′ along the modified strand in the minor groove in both cases. However, subtle orientational differences cause strong distinctions in structural distortions of the DNA duplexes, because the exocyclic amino groups of flanking guanines on both strands compete for space with the BP rings in the minor groove, acting as guideposts for placement of the BP. In the 5′-···CGG*C··· case, the 5′-flanking G · C base pair is severely untwisted, concomitant with a bend deduced from electrophoretic mobility. In the 5′-···CG*GC··· context, there is no untwisting, but there is significant destabilization of the 5′-flanking Watson–Crick base pair. The minor groove width opens near the lesion in both cases, but more for 5′-···CGG*C···. Differential sequence-dependent removal rates of this lesion result and may contribute to the mutation hotspot phenomenon.     

The Individual Blood Cell Telomere Attrition Rate Is Telomere Length Dependent

Age-associated telomere shortening is a well documented feature of peripheral blood cells in human population studies, but it is not known to what extent these data can be transferred to the individual level. Telomere length (TL) in two blood samples taken at ~10 years interval from 959 individuals was investigated using real-time PCR. TL was also measured in 13 families from a multigenerational cohort. As expected, we found an age-related decline in TL over time (r=–0.164, P<0.001, n=959). However, approximately one-third of the individuals exhibited a stable or increased TL over a decade. The individual telomere attrition rate was inversely correlated with initial TL at a highly significant level (r=–0.752, P<0.001), indicating that the attrition rate was most pronounced in individuals with long telomeres at baseline. In accordance, the age-associated telomere attrition rate was more prominent in families with members displaying longer telomeres at a young age (r=–0.691, P<0.001). Abnormal blood TL has been reported at diagnosis of various malignancies, but in the present study there was no association between individual telomere attrition rate or prediagnostic TL and later tumor development. The collected data strongly suggest a TL maintenance mechanism acting in vivo, providing protection of short telomeres as previously demonstrated in vitro. Our findings might challenge the hypothesis that individual TL can predict possible life span or later tumor development.     

Saturation Mutagenesis of Toluene ortho-Monooxygenase of Burkholderia cepacia G4 for Enhanced 1-Naphthol Synthesis and Chloroform Degradation

Directed evolution of toluene ortho-monooxygenase (TOM) of Burkholderia cepacia G4 previously created the hydroxylase α-subunit (TomA3) V106A variant (TOM-Green) with increased activity for both trichloroethylene degradation (twofold enhancement) and naphthalene oxidation (six-times-higher activity). In the present study, saturation mutagenesis was performed at position A106 with Escherichia coli TG1/pBS(Kan)TOMV106A to improve TOM activity for both chloroform degradation and naphthalene oxidation. Whole cells expressing the A106E variant had two times better naphthalene-to-1-naphthol activity than the wild-type cells (V_max of 9.3 versus 4.5 nmol·min⁻¹·mg of protein⁻¹ and unchanged K_m), and the regiospecificity of the A106E variant was unchanged, with 98% 1-naphthol formed, as was confirmed with high-pressure liquid chromatography. The A106E variant degrades its natural substrate toluene 63% faster than wild-type TOM does (2.12 ± 0.07 versus 1.30 ± 0.06 nmol·min⁻¹·mg of protein⁻¹ [mean ± standard deviation]) at 91 μM and has a substantial decrease in regiospecificity, since o-cresol (50%), m-cresol (25%), and p-cresol (25%) are formed, in contrast to the 98% o-cresol formed by wild-type TOM. The A106E variant also has an elevated expression level compared to that of wild-type TOM, as evidenced by sodium dodecyl sulfate-polyacrylamide gel electrophoresis. Another variant, the A106F variant, has 2.8-times-better chloroform degradation activity based on gas chromatography (V_max of 2.61 versus 0.95 nmol·min⁻¹·mg of protein⁻¹ and unchanged K_m) and chloride release (0.034 ± 0.002 versus 0.012 ± 0.001 nmol·min⁻¹·mg of protein⁻¹). The A106F variant also was expressed at levels similar to those of wild-type TOM and 62%-better toluene oxidation activity than wild-type TOM (2.11 ± 0.3 versus 1.30 ± 0.06 nmol·min⁻¹·mg of protein⁻¹). A shift in regiospecificity of toluene hydroxylation was also observed for the A106F variant, with o-cresol (28%), m-cresol (18%), and p-cresol (54%) being formed. Statistical analysis was used to estimate that 292 colonies must be screened for a 99% probability that all 64 codons were sampled during saturation mutagenesis.