Genetic variation in the HLA region is associated with susceptibility to herpes zoster.

Herpes zoster, commonly referred to as shingles, is caused by the varicella zoster virus (VZV). VZV initially manifests as chicken pox, most commonly in childhood, can remain asymptomatically latent in nerve tissues for many years and often re-emerges as shingles. Although reactivation may be related to immune suppression, aging and female sex, most inter-individual variability in re-emergence risk has not been explained to date. We performed a genome-wide association analyses in 22,981 participants (2280 shingles cases) from the electronic Medical Records and Genomics Network. Using Cox survival and logistic regression, we identified a genomic region in the combined and European ancestry groups that has an age of onset effect reaching genome-wide significance (P>1.0 × 10(-8)). This region tags the non-coding gene HCP5 (HLA Complex P5) in the major histocompatibility complex. This gene is an endogenous retrovirus and likely influences viral activity through regulatory functions. Variants in this genetic region are known to be associated with delay in development of AIDS in people infected by HIV. Our study provides further suggestion that this region may have a critical role in viral suppression and could potentially harbor a clinically actionable variant for the shingles vaccine.

Genome-wide association study of recurrent early-onset major depressive disorder.

A genome-wide association study was carried out in 1020 case subjects with recurrent early-onset major depressive disorder (MDD) (onset before age 31) and 1636 control subjects screened to exclude lifetime MDD. Subjects were genotyped with the Affymetrix 6.0 platform. After extensive quality control procedures, 671 424 autosomal single nucleotide polymorphisms (SNPs) and 25 068 X chromosome SNPs with minor allele frequency greater than 1% were available for analysis. An additional 1 892 186 HapMap II SNPs were analyzed based on imputed genotypic data. Single-SNP logistic regression trend tests were computed, with correction for ancestry-informative principal component scores. No genome-wide significant evidence for association was observed, assuming that nominal P<5 × 10(-8) approximates a 5% genome-wide significance threshold. The strongest evidence for association was observed on chromosome 18q22.1 (rs17077540, P=1.83 × 10(-7)) in a region that has produced some evidence for linkage to bipolar-I or -II disorder in several studies, within an mRNA detected in human brain tissue (BC053410) and approximately 75 kb upstream of DSEL. Comparing these results with those of a meta-analysis of three MDD GWAS data sets reported in a companion article, we note that among the strongest signals observed in the GenRED sample, the meta-analysis provided the greatest support (although not at a genome-wide significant level) for association of MDD to SNPs within SP4, a brain-specific transcription factor. Larger samples will be required to confirm the hypothesis of association between MDD (and particularly the recurrent early-onset subtype) and common SNPs.

Genome-wide association study of bipolar disorder in European American and African American individuals.

To identify bipolar disorder (BD) genetic susceptibility factors, we conducted two genome-wide association (GWA) studies: one involving a sample of individuals of European ancestry (EA; n=1001 cases; n=1033 controls), and one involving a sample of individuals of African ancestry (AA; n=345 cases; n=670 controls). For the EA sample, single-nucleotide polymorphisms (SNPs) with the strongest statistical evidence for association included rs5907577 in an intergenic region at Xq27.1 (P=1.6 x 10(-6)) and rs10193871 in NAP5 at 2q21.2 (P=9.8 x 10(-6)). For the AA sample, SNPs with the strongest statistical evidence for association included rs2111504 in DPY19L3 at 19q13.11 (P=1.5 x 10(-6)) and rs2769605 in NTRK2 at 9q21.33 (P=4.5 x 10(-5)). We also investigated whether we could provide support for three regions previously associated with BD, and we showed that the ANK3 region replicates in our sample, along with some support for C15Orf53; other evidence implicates BD candidate genes such as SLITRK2. We also tested the hypothesis that BD susceptibility variants exhibit genetic background-dependent effects. SNPs with the strongest statistical evidence for genetic background effects included rs11208285 in ROR1 at 1p31.3 (P=1.4 x 10(-6)), rs4657247 in RGS5 at 1q23.3 (P=4.1 x 10(-6)), and rs7078071 in BTBD16 at 10q26.13 (P=4.5 x 10(-6)). This study is the first to conduct GWA of BD in individuals of AA and suggests that genetic variations that contribute to BD may vary as a function of ancestry.

Evidence for association of the TCF7 locus with type I diabetes.

The Type I Diabetes Genetics Consortium (T1DGC) has collected thousands of multiplex and simplex families with type I diabetes (T1D) with the goal of identifying genes involved in T1D susceptibility. These families have been genotyped for the HLA class I and class II loci and, recently, for a genome-wide panel of single-nucleotide polymorphisms (SNPs). In addition, multiple SNPs in specific candidate genes have been genotyped in these families in an attempt to evaluate previously reported T1D associations, including the C883A (Pro-Thr) polymorphism in exon 2 of TCF7, a T-cell transcription factor. The TCF7 883A allele was associated with T1D in subjects with T1D not carrying the high-risk HLA genotype DR3/DR4. A panel of 11 SNPs in TCF7 was genotyped in 2092 families from 9 cohorts of the T1DGC. SNPs at two positions in TCF7 were associated with T1D. One associated SNP, C883A (rs5742913), was reported earlier to have a T1D association. A second SNP, rs17653687, represents a novel T1D susceptibility allele in TCF7. After stratification on the high T1D risk DR3/DR4 genotype, the variant (A) allele of C883A was significantly associated with T1D among non-DR3/DR4 cases (transmission=55.8%, P=0.004; OR=1.26) but was not significantly associated in the DR3/DR4 patient subgroup, replicating the earlier report. The reference A allele of intronic SNP rs17653687 was modestly associated with T1D in both DR3/DR4 strata (transmission=54.4% in DR3/DR4; P=0.03; transmission=52.9% in non-DR3/DR4; P=0.03). These results support the previously reported association of the non-synonymous Pro-Thr SNP in TCF7 with T1D, and suggest that other alleles at this locus may also confer risk.

Association analysis of SNPs in the IL4R locus with type I diabetes.

The Type I Diabetes Genetics Consortium (T1DGC) has collected thousands of multiplex and simplex families with type I diabetes (T1D) with the goal of identifying genes involved in T1D susceptibility. These families have all been genotyped for the HLA class I and class II loci and a subset of samples has been typed for an major histocompatibility complex (MHC) single-nucleotide polymorphism (SNP) panel. In addition, the T1DGC has genotyped SNPs in candidate genes to evaluate earlier reported T1D associations. Individual SNPs and SNP haplotypes in IL4R, which encodes the alpha-chain of the IL4 and IL13 receptors, have been associated with T1D in some reports, but not in others. In this study, 38 SNPs in IL4R were genotyped using the Sequenom iPLEX Gold MassARRAY technology in 2042 multiplex families from nine cohorts. Association analyses (transmission-disequilibrium test and parental-disequilibrium test) were performed on individual SNPs and on three-SNP haplotypes. Analyses were also stratified on the high-risk HLA DR3/DR4-DQB1*0302 genotype. A modest T1D association in HBDI families (n=282) was confirmed in this larger collection of HBDI families (n=424). The variant alleles at the non-synonymous SNPs (rs1805011 (E400A), rs1805012 (C431R), and rs1801275 (Q576R)), which are in strong linkage disequilibrium, were negatively associated with T1D risk. These SNPs were more associated with T1D among non-DR3/DR4-DQB1*0302 genotypes than DR3/DR4-DQB1*0302 genotypes. This association was stronger, both in terms of odds ratio and P-values, than the initial report of the smaller collection of HBDI families. However, the IL4R SNPs and the three-SNP haplotype containing the variant alleles were not associated with T1D in the total data. Thus, in the overall families, these results do not show evidence for an association of SNPs in IL4R with T1D.

Overview of the Rapid Response data.

The Type I Diabetes Genetics Consortium (T1DGC) Rapid Response Workshop was established to evaluate published candidate gene associations in a large collection of affected sib-pair (ASP) families. We report on our quality control (QC) and preliminary family-based association analyses. A random sample of blind duplicates was analyzed for QC. Quality checks, including examination of plate-panel yield, marker yield, Hardy-Weinberg equilibrium, mismatch error rate, Mendelian error rate, and allele distribution across plates, were performed. Genotypes from 2324 families within nine cohorts were obtained from a panel of 21 candidate genes, including 384 single-nucleotide polymorphisms on two genotyping platforms performed at the Broad Institute Center for Genotyping and Analysis (Cambridge, MA, USA). The T1DGC Rapid Response project, following rigorous QC procedures, resulted in a 2297 family, 9688 genotyped individual database on a single-candidate gene panel. The available data include 9005 individuals with genotype data from both platforms and 683 individuals genotyped (276 in Illumina; 407 in Sequenom) on only one platform.

Cloning and characterization of the rat and human phosducin-like protein genes: structure, expression and chromosomal localization.

We isolated and characterized the rat gene encoding phosducin-like protein (PhLP), a putative heterotrimeric G protein modulator. The transcription start site was mapped by primer extension. The putative promoter region lacked a TATA sequence but contained a potential initiator element. Two splice variants were identified by RT-PCR of rat brain RNA, potentially generating either the full length or an amino-truncated protein. Only the full-length protein was immunodetected in all mouse tissues surveyed. Comparison of the conceptual translation product of the rat PhLP gene with those from human and Drosophila clones shows a striking conservation in the amino-terminal region of PhLP from these species. This contrasts with the relatively low degree of homology between PhLP and phosducin in this region, suggesting a functional role for this portion of the PhLP protein. Finally, we mapped the human PhLP gene by PCR analysis of somatic cell hybrids and the Stanford G3 radiation hybrid panel. The human PhLP gene (PDCL) is located on chromosome 9, linked to the polymorphic markers D9S1876 and D9S1674 (66-71 cM).

Characterization of the human mitochondrial aconitase gene (ACO2).

We have cloned and characterized the ACO2 gene on human chromosome 22q13, which encodes the essential iron-dependent metabolic enzyme mitochondrial aconitase. We determined that the ACO2 gene comprises 18 translated exons distributed over approximately 35 kilobasepairs (kbp) of DNA. We have shown that the ACO2 mRNA is 2.7kb in length and is expressed ubiquitously, and we can detect multiple isoforms of the ACO2 protein. As others had reported the existence of biochemically active electrophoretic variants of mitochondrial aconitase, we wished to find common ACO2 gene allozymes, functional polymorphisms that might be associated with susceptibility to human genetic diseases. We looked for ACO2 allozymes by DNA sequencing and genotyping in a population of 217 subjects, many of which had idiopathic Parkinson's disease (IPD). We studied patients with IPD because this movement disorder is thought to arise from defects in neuronal iron and energy metabolism, two properties with which aconitase is involved. Furthermore, reports of associations between alleles of the CYP2D6 locus (nearby on 22q13) and IPD, although inconsistent, indicated that an IPD susceptibility locus might be in strong linkage disequilibrium with CYP2D6. We found three functionally silent single nucleotide polymorphisms (SNPs) present in transcribed sequences that exist in similar frequencies in IPD patients and healthy controls. These ACO2 SNPs are in linkage disequilibrium with each other, providing evidence for distinct ACO2 haplotypes. We have, as yet, not detected polymorphisms that would lead to ACO2 allozymes, nor have we observed differences in ACO2 isoform prevalence or distribution in our population of IPD patients and controls. We conclude it is unlikely that polymorphism in the ACO2 gene or post-translational modification of the enzyme predispose to IPD.

The Bacillus subtilis flagellin gene (hag) is transcribed by the sigma 28 form of RNA polymerase.

The Bacillus subtilis gene hag, which codes for the flagellin structural protein, was identified by DNA sequence analysis in a collection of DNA fragments bearing in vitro promoters for the sigma 28 form of RNA polymerase. The hag gene and adjacent regions of the B. subtilis chromosome were restriction mapped, and the nucleotide sequence was determined. The hag gene was transcribed at all stages of growth from a single promoter that had sequences in the promoter recognition region characteristic of the consensus sequence for the sigma 28 holoenzyme. Transcription of hag was eliminated by insertion mutations that blocked synthesis of the sigma 28 protein. These findings provide strong support for the previous proposal that the sigma 28 form of RNA polymerase controls transcription of a regulon specifying flagellar, chemotaxis, and motility functions in B. subtilis (J. D. Helmann and M. J. Chamberlin, Proc. Natl. Acad. Sci. USA 84:6422-6424, 1987). The steady-state levels of hag mRNA increased during exponential growth and peaked as the B. subtilis cells entered the stationary phase. The transcript levels then decreased to zero within 4 h after the onset of sporulation. Hence, sigma 28 RNA polymerase function is temporally regulated.

Identification of flagellar synthesis regulatory and structural genes in a sigma D-dependent operon of Bacillus subtilis.

The sigma D form of RNA polymerase from Bacillus subtilis has been shown previously to direct the synthesis of several transcription units bearing genes for flagellin, motility proteins, and autolysins. In this report, we describe an operon of genes transcribed from the sigma D-dependent promoter PD-1. We have identified three complete open reading frames and one partial one downstream of this promoter; immediately upstream is the previously identified comF locus. The PD-1 operon encodes the presumptive B. subtilis homologs of two Salmonella typhimurium late flagellar genes, flgM and flgK. Also present in this operon are two genes of unknown function, orf139 and orf160, whose products show similarities to the eukaryotic cytoskeletal proteins myosin and vimentin, respectively. orf139 and orf160 may encode proteins that form extended alpha-helical secondary structures and coiled-coil quaternary structures which may be filamentous components of the gram-positive bacterial flagellum. We have characterized the B. subtilis flgM gene further by constructing an in-frame deletion mutation, flgM delta 80, and creating strains of B. subtilis in which this allele has replaced the wild-type copy. By primer extension analysis of cellular RNA, we have shown that the flgM delta 80 mutation relieves the block to transcription of two other sigma D-dependent operons imposed by an unlinked mutation in a gene directing early flagellar synthesis. We conclude that, as in the case of S. typhimurium, early flagellar synthesis in B. subtilis is coupled to late flagellar synthesis through repression of sigma D-dependent transcription by the flgM gene product.

An operon of Bacillus subtilis motility genes transcribed by the sigma D form of RNA polymerase.

Two genes controlling motility functions in Bacillus subtilis were identified by DNA sequence analysis of a chromosomal fragment containing a strong promoter for sigma D RNA polymerase. Previous studies had shown that this sigma D-dependent promoter controls synthesis of a 1.6-kb transcript in vivo and in vitro. Sequence analysis revealed that the 1.6-kb transcript contains two open reading frames coding for protein sequences homologous to the Escherichia coli motA and motB gene products, respectively, and ends in a rho-independent termination site. Direct evidence linking these genes to motility functions in B. subtilis was obtained by precise localization by polymerase chain reaction of Tn917 transposon insertion mutations of Mot- strains, isolated by Zuberi et al. (A. R. Zuberi, C. Ying, H. M. Parker, and G. W. Ordal, J. Bacteriol. 172:6841-6848, 1990), to within this mot. operon. Replacement of each wild-type gene by in-frame deletion mutations yielded strains possessing paralyzed flagella and confirmed that both motA and motB are required for the motility of B. subtilis. These current findings support our earlier suggestions that sigma D in B. subtilis plays a central role in the control of gene expression for flagellar assembly, chemotaxis, and motility functions. Sigma F, the enteric homolog of sigma D, controls similar functions in E. coli and Salmonella typhimurium, and these factors appear to be representative of a family of factors implicated in flagellar synthesis in many bacterial species, which we propose to designate the sigma 28 family.

Regulation of sigma D expression and activity by spo0, abrB, and sin gene products in Bacillus subtilis.

Expression of sigma D protein and of the hag gene, which is transcribed by the sigma D holoenzyme, is not dependent on spo0, abrB, or sin gene products in Bacillus subtilis. Preliminary results, however, suggest that a signal mediated by the spo0K locus may be responsible for the inhibition of sigma D activity during the stationary phase.

Environmental regulation of Bacillus subtilis sigma(D)-dependent gene expression.

The sigma(D) regulon of Bacillus subtilis is composed of genes encoding proteins for flagellar synthesis, motility, and chemotaxis. Concurrent analyses of sigma(D) protein levels and flagellin mRNA demonstrate that sigD expression and sigma(D) activity are tightly coupled during growth in both complex and minimal media, although they exhibit different patterns of expression. We therefore used the sigma(D)-dependent flagellin gene (hag) as a model gene to study the effects of different nutritional environments on sigma(D)-dependent gene expression. In complex medium, the level of expression of a hag-lacZ fusion increased exponentially during the exponential growth phase and peaked early in the transition state. In contrast, the level of expression of this reporter remained constant and high throughout growth in minimal medium. These results suggest the existence of a nutritional signal(s) that affects sigD expression and/or sigma(D) activity. This signal(s) allows for nutritional repression early in growth and, based on reconstitution studies, resides in the complex components of sporulation medium, as well as in a mixture of mono-amino acids. However, the addition of Casamino Acids to minimal medium results in a dose-dependent decrease in hag-lacZ expression throughout growth and the postexponential growth phase. In work by others, CodY has been implicated in the nutritional repression of several genes. Analysis of a codY mutant bearing a hag-lacZ reporter revealed that flagellin expression is released from nutritional repression in this strain, whereas mutations in the transition state preventor genes abrB, hpr, and sinR failed to elicit a similar effect during growth in complex medium. Therefore, the CodY protein appears to be the physiologically relevant regulator of hag nutritional repression in B. subtilis.

Analysis of IL4R haplotypes in predisposition to multiple sclerosis.

We have investigated the association of multiple sclerosis (MS) with polymorphisms in the IL4R gene in 332 single-case MS families. IL4R encodes a subunit of the interleukin-4 receptor, a molecule important for T-cell development and differentiation, and is a gene shown to be associated with immune-related diseases such as asthma and type I diabetes. By genotyping two promoter and eight coding IL4R SNPs and identifying haplotypes (complex alleles) in the MS families, stratified for HLA genotype, we have observed evidence of the association of the IL4R gene to MS. In particular, we have identified a specific susceptibility haplotype, and observe that the risk is conferred primarily to individuals not carrying the high MS-risk HLA DR2 (DRB1(*)1501-DQB1(*)0602) haplotype (nominal P=0.009). These findings suggest a potentially important role for the IL4R gene in predisposition to MS, and provide further evidence of its relevance as a candidate gene for immune-related diseases.

Is there a genetic susceptibility locus for Parkinson's disease on chromosome 22q13?

The cytochrome P450 mono-oxygenase gene, CYP2D6 on chromosome 22q13 (ch22q13), has been inconsistently associated with Parkinson's disease. Associations with CYP2D6 have either been absent altogether or have involved more than one polymorphism, many of which have the same metabolic effect on gene expression. We examined the association between CYP2D6 polymorphisms and Parkinson's disease in a case-control study and included 10 polymorphic dinucleotide repeat markers linked to CYP2D6 to determine whether the association was present or due to linkage disequilibrium. There was no association between any polymorphism of CYP2D6 and Parkinson's disease, but two of 10 dinucleotide repeat markers linked to CYP2D6 were associated with the disease. These results provide evidence to suggest that there may be an unidentified locus for susceptibility to Parkinson's disease that is in linkage disequilibrium with dinucleotide repeat markers mapping near CYP2D6 on ch22q13.