Identification of Non-HLA Genes Associated with Celiac Disease and Country-Specific Differences in a Large,International Pediatric Cohort

Objectives

There are significant geographical differences in the prevalence and incidence of celiac disease that cannot be explained by HLA alone. More than 40 loci outside of the HLA region have been associated with celiac disease. We investigated the roles of these non-HLA genes in the development of tissue transglutaminase autoantibodies (tTGA) and celiac disease in a large international prospective cohort study.

Methods

A total of 424,788 newborns from the US and European general populations and first-degree relatives with type 1 diabetes were screened for specific HLA genotypes. Of these, 21,589 carried 1 of the 9 HLA genotypes associated with increased risk for type 1 diabetes and celiac disease; we followed 8676 of the children in a 15 y prospective follow-up study. Genotype analyses were performed on 6010 children using the Illumina ImmunoChip. Levels of tTGA were measured in serum samples using radio-ligand binding assays; diagnoses of celiac disease were made based on persistent detection of tTGA and biopsy analysis. Data were analyzed using Cox proportional hazards analyses.

Results

We found 54 single-nucleotide polymorphisms (SNPs) in 5 genes associated with celiac disease (TAGAP, IL18R1, RGS21, PLEK, and CCR9) in time to celiac disease analyses (10⁻⁴>P>5.8x10⁻⁶). The hazard ratios (HR) for the SNPs with the smallest P values in each region were 1.59, 1.45, 2.23, 2.64, and 1.40, respectively. Outside of regions previously associated with celiac disease, we identified 10 SNPs in 8 regions that could also be associated with the disease (P<10⁻⁴). A SNP near PKIA (rs117128341, P = 6.5x10⁻⁸, HR = 2.8) and a SNP near PFKFB3 (rs117139146, P<2.8x10⁻⁷, HR = 4.9) reached the genome-wide association threshold in subjects from Sweden. Analyses of time to detection of tTGA identified 29 SNPs in 2 regions previously associated with celiac disease (CTLA4, P = 1.3x10⁻⁶, HR = 0.76 and LPP, P = 2.8x10⁻⁵, HR = .80) and 6 SNPs in 5 regions not previously associated with celiac disease (P<10⁻⁴); non-HLA genes are therefore involved in development of tTGA.

Conclusions

In conclusion, using a genetic analysis of a large international cohort of children, we associated celiac disease development with 5 non-HLA regions previously associated with the disease and 8 regions not previously associated with celiac disease. We identified 5 regions associated with development of tTGA. Two loci associated with celiac disease progression reached a genome-wide association threshold in subjects from Sweden.     

Genome-Wide Association Study of Pancreatic Cancer in Japanese Population

Pancreatic cancer shows very poor prognosis and is the fifth leading cause of cancer death in Japan. Previous studies indicated some genetic factors contributing to the development and progression of pancreatic cancer; however, there are limited reports for common genetic variants to be associated with this disease, especially in the Asian population. We have conducted a genome-wide association study (GWAS) using 991 invasive pancreatic ductal adenocarcinoma cases and 5,209 controls, and identified three loci showing significant association (P-value<5×10⁻⁷) with susceptibility to pancreatic cancer. The SNPs that showed significant association carried estimated odds ratios of 1.29, 1.32, and 3.73 with 95% confidence intervals of 1.17–1.43, 1.19–1.47, and 2.24–6.21; P-value of 3.30×10⁻⁷, 3.30×10⁻⁷, and 4.41×10⁻⁷; located on chromosomes 6p25.3, 12p11.21 and 7q36.2, respectively. These associated SNPs are located within linkage disequilibrium blocks containing genes that have been implicated some roles in the oncogenesis of pancreatic cancer.     

Genome-Wide Association Study for Incident Myocardial Infarction and Coronary Heart Disease in Prospective Cohort Studies: The CHARGE Consortium

Background

Data are limited on genome-wide association studies (GWAS) for incident coronary heart disease (CHD). Moreover, it is not known whether genetic variants identified to date also associate with risk of CHD in a prospective setting.

Methods

We performed a two-stage GWAS analysis of incident myocardial infarction (MI) and CHD in a total of 64,297 individuals (including 3898 MI cases, 5465 CHD cases). SNPs that passed an arbitrary threshold of 5×10⁻⁶ in Stage I were taken to Stage II for further discovery. Furthermore, in an analysis of prognosis, we studied whether known SNPs from former GWAS were associated with total mortality in individuals who experienced MI during follow-up.

Results

In Stage I 15 loci passed the threshold of 5×10⁻⁶; 8 loci for MI and 8 loci for CHD, for which one locus overlapped and none were reported in previous GWAS meta-analyses. We took 60 SNPs representing these 15 loci to Stage II of discovery. Four SNPs near QKI showed nominally significant association with MI (p-value<8.8×10⁻³) and three exceeded the genome-wide significance threshold when Stage I and Stage II results were combined (top SNP rs6941513: p = 6.2×10⁻⁹). Despite excellent power, the 9p21 locus SNP (rs1333049) was only modestly associated with MI (HR = 1.09, p-value = 0.02) and marginally with CHD (HR = 1.06, p-value = 0.08). Among an inception cohort of those who experienced MI during follow-up, the risk allele of rs1333049 was associated with a decreased risk of subsequent mortality (HR = 0.90, p-value = 3.2×10⁻³).

Conclusions

QKI represents a novel locus that may serve as a predictor of incident CHD in prospective studies. The association of the 9p21 locus both with increased risk of first myocardial infarction and longer survival after MI highlights the importance of study design in investigating genetic determinants of complex disorders.     

Association of Genetic Variants at 3q22 with Nephropathy in Patients with Type 1 Diabetes Mellitus

Diabetic nephropathy (DN) is the primary cause of morbidity and mortality in patients with type 1 diabetes mellitus (T1DM) and affects about 30% of these patients. We have previously localized a DN locus on chromosome 3q with suggestive linkage in Finnish individuals. Linkage to this region has also been reported earlier by several other groups. To fine map this locus, we conducted a multistage case-control association study in T1DM patients, comprising 1822 cases with nephropathy and 1874 T1DM patients free of nephropathy, from Finland, Iceland, and the British Isles. At the screening stage, we genotyped 3072 tag SNPs, spanning a 28 Mb region, in 234 patients and 215 controls from Finland. SNPs that met the significance threshold of p < 0.01 at this stage were followed up by a series of sample sets. A genetic variant, rs1866813, in the noncoding region at 3q22 was associated with increased risk of DN (overall p = 7.07 × 10⁻⁶, combined odds ratio [OR] of the allele = 1.33). The estimated genotypic ORs of this variant in all Finnish samples suggested a codominant effect, resulting in significant association, with a p value of 4.7 × 10⁻⁵ (OR = 1.38; 95% confidence interval = 1.18–1.62). Additionally, an 11 kb segment flanked by rs62408925 and rs1866813, two strongly correlated variants (r² = 0.95), contains three elements highly conserved across multiple species. Independent replication will clarify the role of the associated variants at 3q22 in influencing the risk of DN.     

Susceptibility Loci Associated with Specific and Shared Subtypes of Lymphoid Malignancies

The genetics of lymphoma susceptibility reflect the marked heterogeneity of diseases that comprise this broad phenotype. However, multiple subtypes of lymphoma are observed in some families, suggesting shared pathways of genetic predisposition to these pathologically distinct entities. Using a two-stage GWAS, we tested 530,583 SNPs in 944 cases of lymphoma, including 282 familial cases, and 4,044 public shared controls, followed by genotyping of 50 SNPs in 1,245 cases and 2,596 controls. A novel region on 11q12.1 showed association with combined lymphoma (LYM) subtypes. SNPs in this region included rs12289961 near LPXN, (P_LYM = 3.89×10⁻⁸, OR = 1.29) and rs948562 (P_LYM = 5.85×10⁻⁷, OR = 1.29). A SNP in a novel non-HLA region on 6p23 (rs707824, P_NHL = 5.72×10⁻⁷) was suggestive of an association conferring susceptibility to lymphoma. Four SNPs, all in a previously reported HLA region, 6p21.32, showed genome-wide significant associations with follicular lymphoma. The most significant association with follicular lymphoma was for rs4530903 (P_FL = 2.69×10⁻¹², OR = 1.93). Three novel SNPs near the HLA locus, rs9268853, rs2647046, and rs2621416, demonstrated additional variation contributing toward genetic susceptibility to FL associated with this region. Genes implicated by GWAS were also found to be cis-eQTLs in lymphoblastoid cell lines; candidate genes in these regions have been implicated in hematopoiesis and immune function. These results, showing novel susceptibility regions and allelic heterogeneity, point to the existence of pathways of susceptibility to both shared as well as specific subtypes of lymphoid malignancy.     

A Genome Wide Association Study of Mathematical Ability Reveals an Association at Chromosome 3q29, a Locus Associated with Autism and Learning Difficulties: A Preliminary Study

Mathematical ability is heritable, but few studies have directly investigated its molecular genetic basis. Here we aimed to identify specific genetic contributions to variation in mathematical ability. We carried out a genome wide association scan using pooled DNA in two groups of U.K. samples, based on end of secondary/high school national academic exam achievement: high (n = 419) versus low (n = 183) mathematical ability while controlling for their verbal ability. Significant differences in allele frequencies between these groups were searched for in 906,600 SNPs using the Affymetrix GeneChip Human Mapping version 6.0 array. After meeting a threshold of p<1.5×10⁻⁵, 12 SNPs from the pooled association analysis were individually genotyped in 542 of the participants and analyzed to validate the initial associations (lowest p-value 1.14 ×10⁻⁶). In this analysis, one of the SNPs (rs789859) showed significant association after Bonferroni correction, and four (rs10873824, rs4144887, rs12130910 rs2809115) were nominally significant (lowest p-value 3.278 × 10⁻⁴). Three of the SNPs of interest are located within, or near to, known genes (FAM43A, SFT2D1, C14orf64). The SNP that showed the strongest association, rs789859, is located in a region on chromosome 3q29 that has been previously linked to learning difficulties and autism. rs789859 lies 1.3 kbp downstream of LSG1, and 700 bp upstream of FAM43A, mapping within the potential promoter/regulatory region of the latter. To our knowledge, this is only the second study to investigate the association of genetic variants with mathematical ability, and it highlights a number of interesting markers for future study.     

Genome-Wide Association Study of Maternal and Inherited Loci for Conotruncal Heart Defects

Conotruncal and related heart defects (CTDs) are a group of serious and relatively common birth defects. Although both maternal and inherited genotypes are thought to play a role in the etiology of CTDs, few specific genetic risk factors have been identified. To determine whether common variants acting through the genotype of the mother (e.g. via an in utero effect) or the case are associated with CTDs, we conducted a genome-wide association study of 750 CTD case-parent triads, with follow-up analyses in 358 independent triads. Log-linear analyses were used to assess the association of CTDs with the genotypes of both the mother and case. No association achieved genomewide significance in either the discovery or combined (discovery+follow-up) samples. However, three loci with p-values suggestive of association (p<10⁻⁵) in the discovery sample had p-values <0.05 in the follow-up sample and p-values in the combined data that were lower than in the discovery sample. These included suggestive association with an inherited intergenic variant at 20p12.3 (rs6140038, combined p = 1.0×10⁻⁵) and an inherited intronic variant in KCNJ4 at 22q13.1 (rs2267386, combined p = 9.8×10⁻⁶), as well as with a maternal variant in SLC22A24 at 11q12.3 (rs11231379, combined p = 4.2×10⁻⁶). These observations suggest novel candidate loci for CTDs, including loci that appear to be associated with the risk of CTDs via the maternal genotype, but further studies are needed to confirm these associations.     

Fine Mapping on Chromosome 10q22-q23 Implicates Neuregulin 3 in Schizophrenia

Linkage studies have implicated 10q22-q23 as a schizophrenia (SZ) susceptibility locus in Ashkenazi Jewish (AJ) and Han Chinese from Taiwan populations. To further explore our previous linkage signal in the AJ population (NPL score: 4.27, empirical p = 2 × 10⁻⁵), we performed a peakwide association fine mapping study by using 1414 SNPs across ~12.5 Mb in 10q22-q23. We genotyped 1515 AJ individuals, including 285 parent-child trios, 173 unrelated cases, and 487 unrelated controls. We analyzed the binary diagnostic phenotype of SZ and 9 heritable quantitative traits derived from a principal components factor analysis of 73 items from our consensus diagnostic ratings and direct assessment interviews. Although no marker withstood multiple test correction for association with the binary SZ phenotype, we found strong evidence of association by using the “delusion” factor as the quantitative trait at three SNPs (rs10883866, rs10748842, and rs6584400) located in a 13 kb interval in intron 1 of Neuregulin 3 (NRG3). Our best p value from family-based association analysis was 7.26 × 10⁻⁷. We replicated this association in the collection of 173 unrelated AJ cases (p = 1.55 × 10⁻²), with a combined p value of 2.30 × 10⁻⁷. After performing 10,000 permutations of each of the phenotypes, we estimated the empirical study-wide significance across all 9 factors (90,000 permutations) to be p = 2.7 × 10⁻³. NRG3 is primarily expressed in the central nervous system and is one of three paralogs of NRG1, a gene strongly implicated in SZ. These biological properties together with our linkage and association results strongly support NRG3 as a gene involved in SZ.     

Novel Crohn disease locus identified by genome-wide association maps to a gene desert on 5p13.1 and modulates expression of PTGER4

To identify novel susceptibility loci for Crohn disease (CD), we undertook a genome-wide association study with more than 300,000 SNPs characterized in 547 patients and 928 controls. We found three chromosome regions that provided evidence of disease association with p-values between 10⁻⁶ and 10⁻⁹. Two of these (IL23R on Chromosome 1 and CARD15 on Chromosome 16) correspond to genes previously reported to be associated with CD. In addition, a 250-kb region of Chromosome 5p13.1 was found to contain multiple markers with strongly suggestive evidence of disease association (including four markers with p < 10⁻⁷). We replicated the results for 5p13.1 by studying 1,266 additional CD patients, 559 additional controls, and 428 trios. Significant evidence of association (p < 4 × 10⁻⁴) was found in case/control comparisons with the replication data, while associated alleles were over-transmitted to affected offspring (p < 0.05), thus confirming that the 5p13.1 locus contributes to CD susceptibility. The CD-associated 250-kb region was saturated with 111 SNP markers. Haplotype analysis supports a complex locus architecture with multiple variants contributing to disease susceptibility. The novel 5p13.1 CD locus is contained within a 1.25-Mb gene desert. We present evidence that disease-associated alleles correlate with quantitative expression levels of the prostaglandin receptor EP4, PTGER4, the gene that resides closest to the associated region. Our results identify a major new susceptibility locus for CD, and suggest that genetic variants associated with disease risk at this locus could modulate cis-acting regulatory elements of PTGER4.     

Variants in TF and HFE Explain ∼40% of Genetic Variation in Serum-Transferrin Levels

Only a small proportion of genetic variation in complex traits has been explained by SNPs from genome-wide association studies (GWASs). We report the results from two GWASs for serum markers of iron status (serum iron, serum transferrin, transferrin saturation with iron, and serum ferritin), which are important in iron overload (e.g., hemochromatosis) and deficiency (e.g., anemia) conditions. We performed two GWASs on samples of Australians of European descent. In the first GWAS, 411 adolescent twins and their siblings were genotyped with 100K SNPs. rs1830084, 10.8 kb 3′ of TF, was significantly associated with serum transferrin (p total association test = 1.0 × 10⁻⁹; p within-family test = 2.2 × 10⁻⁵). In the second GWAS on an independent sample of 459 female monozygotic (MZ) twin pairs genotyped with 300K SNPs, we found rs3811647 (within intron 11 of TF, HapMap CEU r² with rs1830084 = 0.86) was significantly associated with serum transferrin (p = 3.0 × 10⁻¹⁵). In the second GWAS, we found two additional and independent SNPs on TF (rs1799852 and rs2280673) and confirmed the known C282Y mutation in HFE to be independently associated with serum transferrin. The three variants in TF (rs3811647, rs1799852 and rs2280673) plus the HFE C282Y mutation explained ~40% of genetic variation in serum transferrin (p = 7.8 × 10⁻²⁵). These findings are potentially important for our understanding of iron metabolism and of regulation of hepatic protein secretion, and also strongly support the hypothesis that the genetic architecture of some endophenotypes may be simpler than that of disease.     

Identification of New Genetic Risk Variants for Type 2 Diabetes

Although more than 20 genetic susceptibility loci have been reported for type 2 diabetes (T2D), most reported variants have small to moderate effects and account for only a small proportion of the heritability of T2D, suggesting that the majority of inter-person genetic variation in this disease remains to be determined. We conducted a multistage, genome-wide association study (GWAS) within the Asian Consortium of Diabetes to search for T2D susceptibility markers. From 590,887 SNPs genotyped in 1,019 T2D cases and 1,710 controls selected from Chinese women in Shanghai, we selected the top 2,100 SNPs that were not in linkage disequilibrium (r²<0.2) with known T2D loci for in silico replication in three T2D GWAS conducted among European Americans, Koreans, and Singapore Chinese. The 5 most promising SNPs were genotyped in an independent set of 1,645 cases and 1,649 controls from Shanghai, and 4 of them were further genotyped in 1,487 cases and 3,316 controls from 2 additional Chinese studies. Consistent associations across all studies were found for rs1359790 (13q31.1), rs10906115 (10p13), and rs1436955 (15q22.2) with P-values (per allele OR, 95%CI) of 6.49×10⁻⁹ (1.15, 1.10–1.20), 1.45×10⁻⁸ (1.13, 1.08–1.18), and 7.14×10⁻⁷ (1.13, 1.08–1.19), respectively, in combined analyses of 9,794 cases and 14,615 controls. Our study provides strong evidence for a novel T2D susceptibility locus at 13q31.1 and the presence of new independent risk variants near regions (10p13 and 15q22.2) reported by previous GWAS.     

Replication of putative susceptibility loci from genome-wide association studies associated with coronary atherosclerosis in Chinese Han population

Background

Coronary atherosclerosis, the main cause of cardiovascular disease, is a progressive disease. Recent Genome Wide Association Studies (GWASs) discovered several novel loci associated with coronary artery disease (CAD) or its main complication myocardial infarction (MI). In this study, we investigated the associations between previously reported CAD- and MI-associated variants and coronary atherosclerosis in Chinese Han population.

Methodology/Principal Findings

We performed a case-control association study with 2,335 coronary atherosclerosis patients and 1,078 controls undergoing coronary angiography of Chinese Han from China. Fourteen single nucleotide polymorphisms (SNPs), located at 1p13.3, 1q41, 2q36.3, 6q25.1, 9p21.3, 10q11.21 and 15q22.33, were genotyped in our sample collection. Six SNPs at 9p21 were associated with coronary atherosclerosis susceptibility (P_trend<0.05) and rs10757274 showed the most significant association (P = 2.38×10⁻⁰⁸, OR = 1.34). These associations remained significant after adjustment for multiple comparisons. Rs17465637 at 1q41 (P_trend = 6.83×10⁻⁰³, OR = 0.86) also showed significant association with coronary atherosclerosis, but the association was not significant after multiple comparisons. Additionally, rs501120 (P = 8.36×10⁻⁰³, OR = 0.80) at 10q11.21 was associated with coronary atherosclerosis in females, but did not show association in males and all participants. Variants at 1p13.3, 2q36.3, 6q25.1 and 15q22.33 showed no associations with coronary atherosclerosis and main cardiovascular risk factors in our data.

Conclusions/Significance

Our findings indicated variants at 9p21 were significantly associated with coronary atherosclerosis in Han Chinese. Variants at 1q41 showed suggestive evidence of association and variants at 10q11.21 showed suggestive evidence of association in females, which warrant further study in a larger sample.     

Genome-Wide Search for Gene-Gene Interactions in Colorectal Cancer

Genome-wide association studies (GWAS) have successfully identified a number of single-nucleotide polymorphisms (SNPs) associated with colorectal cancer (CRC) risk. However, these susceptibility loci known today explain only a small fraction of the genetic risk. Gene-gene interaction (GxG) is considered to be one source of the missing heritability. To address this, we performed a genome-wide search for pair-wise GxG associated with CRC risk using 8,380 cases and 10,558 controls in the discovery phase and 2,527 cases and 2,658 controls in the replication phase. We developed a simple, but powerful method for testing interaction, which we term the Average Risk Due to Interaction (ARDI). With this method, we conducted a genome-wide search to identify SNPs showing evidence for GxG with previously identified CRC susceptibility loci from 14 independent regions. We also conducted a genome-wide search for GxG using the marginal association screening and examining interaction among SNPs that pass the screening threshold (p<10⁻⁴). For the known locus rs10795668 (10p14), we found an interacting SNP rs367615 (5q21) with replication p = 0.01 and combined p = 4.19×10⁻⁸. Among the top marginal SNPs after LD pruning (n = 163), we identified an interaction between rs1571218 (20p12.3) and rs10879357 (12q21.1) (nominal combined p = 2.51×10⁻⁶; Bonferroni adjusted p = 0.03). Our study represents the first comprehensive search for GxG in CRC, and our results may provide new insight into the genetic etiology of CRC.     

Genome-wide association scan shows genetic variants in the FTO gene are associated with obesity-related traits

The obesity epidemic is responsible for a substantial economic burden in developed countries and is a major risk factor for type 2 diabetes and cardiovascular disease. The disease is the result not only of several environmental risk factors, but also of genetic predisposition. To take advantage of recent advances in gene-mapping technology, we executed a genome-wide association scan to identify genetic variants associated with obesity-related quantitative traits in the genetically isolated population of Sardinia. Initial analysis suggested that several SNPs in the FTO and PFKP genes were associated with increased BMI, hip circumference, and weight. Within the FTO gene, rs9930506 showed the strongest association with BMI (p = 8.6 ×10⁻⁷), hip circumference (p = 3.4 × 10⁻⁸), and weight (p = 9.1 × 10⁻⁷). In Sardinia, homozygotes for the rare “G” allele of this SNP (minor allele frequency = 0.46) were 1.3 BMI units heavier than homozygotes for the common “A” allele. Within the PFKP gene, rs6602024 showed very strong association with BMI (p = 4.9 × 10⁻⁶). Homozygotes for the rare “A” allele of this SNP (minor allele frequency = 0.12) were 1.8 BMI units heavier than homozygotes for the common “G” allele. To replicate our findings, we genotyped these two SNPs in the GenNet study. In European Americans (N = 1,496) and in Hispanic Americans (N = 839), we replicated significant association between rs9930506 in the FTO gene and BMI (p-value for meta-analysis of European American and Hispanic American follow-up samples, p = 0.001), weight (p = 0.001), and hip circumference (p = 0.0005). We did not replicate association between rs6602024 and obesity-related traits in the GenNet sample, although we found that in European Americans, Hispanic Americans, and African Americans, homozygotes for the rare “A” allele were, on average, 1.0–3.0 BMI units heavier than homozygotes for the more common “G” allele. In summary, we have completed a whole genome–association scan for three obesity-related quantitative traits and report that common genetic variants in the FTO gene are associated with substantial changes in BMI, hip circumference, and body weight. These changes could have a significant impact on the risk of obesity-related morbidity in the general population.     

Genome-wide Association Study Identifies Five Susceptibility Loci for Follicular Lymphoma outside the HLA Region

Genome-wide association studies (GWASs) of follicular lymphoma (FL) have previously identified human leukocyte antigen (HLA) gene variants. To identify additional FL susceptibility loci, we conducted a large-scale two-stage GWAS in 4,523 case subjects and 13,344 control subjects of European ancestry. Five non-HLA loci were associated with FL risk: 11q23.3 (rs4938573, p = 5.79 × 10⁻²⁰) near CXCR5; 11q24.3 (rs4937362, p = 6.76 × 10⁻¹¹) near ETS1; 3q28 (rs6444305, p = 1.10 × 10⁻¹⁰) in LPP; 18q21.33 (rs17749561, p = 8.28 × 10⁻¹⁰) near BCL2; and 8q24.21 (rs13254990, p = 1.06 × 10⁻⁸) near PVT1. In an analysis of the HLA region, we identified four linked HLA-DRβ1 multiallelic amino acids at positions 11, 13, 28, and 30 that were associated with FL risk (p_omnibus = 4.20 × 10⁻⁶⁷ to 2.67 × 10⁻⁷⁰). Additional independent signals included rs17203612 in HLA class II (odds ratio [OR_per-allele] = 1.44; p = 4.59 × 10⁻¹⁶) and rs3130437 in HLA class I (OR_per-allele = 1.23; p = 8.23 × 10⁻⁹). Our findings further expand the number of loci associated with FL and provide evidence that multiple common variants outside the HLA region make a significant contribution to FL risk.     

Genome Wide Linkage Study,Using a 250K SNP Map,of Plasmodium falciparum Infection and Mild Malaria Attack in a Senegalese Population

Multiple factors are involved in the variability of host''s response to P. falciparum infection, like the intensity and seasonality of malaria transmission, the virulence of parasite and host characteristics like age or genetic make-up. Although admitted nowadays, the involvement of host genetic factors remains unclear. Discordant results exist, even concerning the best-known malaria resistance genes that determine the structure or function of red blood cells. Here we report on a genome-wide linkage and association study for P. falciparum infection intensity and mild malaria attack among a Senegalese population of children and young adults from 2 to 18 years old. A high density single nucleotide polymorphisms (SNP) genome scan (Affimetrix GeneChip Human Mapping 250K-nsp) was performed for 626 individuals: i.e. 249 parents and 377 children out of the 504 ones included in the follow-up. The population belongs to a unique ethnic group and was closely followed-up during 3 years. Genome-wide linkage analyses were performed on four clinical and parasitological phenotypes and association analyses using the family based association tests (FBAT) method were carried out in regions previously linked to malaria phenotypes in literature and in the regions for which we identified a linkage peak. Analyses revealed three strongly suggestive evidences for linkage: between mild malaria attack and both the 6p25.1 and the 12q22 regions (empirical p-value = 5×10⁻⁵ and 9×10⁻⁵ respectively), and between the 20p11q11 region and the prevalence of parasite density in asymptomatic children (empirical p-value = 1.5×10⁻⁴). Family based association analysis pointed out one significant association between the intensity of plasmodial infection and a polymorphism located in ARHGAP26 gene in the 5q31–q33 region (p-value = 3.7×10⁻⁵). This study identified three candidate regions, two of them containing genes that could point out new pathways implicated in the response to malaria infection. Furthermore, we detected one gene associated with malaria infection in the 5q31–q33 region.     

Meta-analysis of Genome-wide Association Studies Identifies 1q22 as a Susceptibility Locus for Intracerebral Hemorrhage

Intracerebral hemorrhage (ICH) is the stroke subtype with the worst prognosis and has no established acute treatment. ICH is classified as lobar or nonlobar based on the location of ruptured blood vessels within the brain. These different locations also signal different underlying vascular pathologies. Heritability estimates indicate a substantial genetic contribution to risk of ICH in both locations. We report a genome-wide association study of this condition that meta-analyzed data from six studies that enrolled individuals of European ancestry. Case subjects were ascertained by neurologists blinded to genotype data and classified as lobar or nonlobar based on brain computed tomography. ICH-free control subjects were sampled from ambulatory clinics or random digit dialing. Replication of signals identified in the discovery cohort with p < 1 × 10⁻⁶ was pursued in an independent multiethnic sample utilizing both direct and genome-wide genotyping. The discovery phase included a case cohort of 1,545 individuals (664 lobar and 881 nonlobar cases) and a control cohort of 1,481 individuals and identified two susceptibility loci: for lobar ICH, chromosomal region 12q21.1 (rs11179580, odds ratio [OR] = 1.56, p = 7.0 × 10⁻⁸); and for nonlobar ICH, chromosomal region 1q22 (rs2984613, OR = 1.44, p = 1.6 × 10⁻⁸). The replication included a case cohort of 1,681 individuals (484 lobar and 1,194 nonlobar cases) and a control cohort of 2,261 individuals and corroborated the association for 1q22 (p = 6.5 × 10⁻⁴; meta-analysis p = 2.2 × 10⁻¹⁰) but not for 12q21.1 (p = 0.55; meta-analysis p = 2.6 × 10⁻⁵). These results demonstrate biological heterogeneity across ICH subtypes and highlight the importance of ascertaining ICH cases accordingly.     

Comprehensive Linkage and Association Analyses Identify Haplotype, Near to the TNFSF15 Gene, Significantly Associated with Spondyloarthritis

Spondyloarthritis (SpA) is a chronic inflammatory disorder with a strong genetic predisposition dominated by the role of HLA-B27. However, the contribution of other genes to the disease susceptibility has been clearly demonstrated. We previously reported significant evidence of linkage of SpA to chromosome 9q31–34. The current study aimed to characterize this locus, named SPA2. First, we performed a fine linkage mapping of SPA2 (24 cM) with 28 microsatellite markers in 149 multiplex families, which allowed us to reduce the area of investigation to an 18 cM (13 Mb) locus delimited by the markers D9S279 and D9S112. Second, we constructed a linkage disequilibrium (LD) map of this region with 1,536 tag single-nucleotide polymorphisms (SNPs) in 136 families (263 patients). The association was assessed using a transmission disequilibrium test. One tag SNP, rs4979459, yielded a significant P-value (4.9×10⁻⁵). Third, we performed an extension association study with rs4979459 and 30 surrounding SNPs in LD with it, in 287 families (668 patients), and in a sample of 139 cases and 163 controls. Strong association was observed in both familial and case/control datasets for several SNPs. In the replication study, carried with 8 SNPs in an independent sample of 232 cases and 149 controls, one SNP, rs6478105, yielded a nominal P-value<3×10⁻². Pooled case/control study (371 cases and 312 controls) as well as combined analysis of extension and replication data showed very significant association (P<5×10⁻⁴) for 6 of the 8 latter markers (rs7849556, rs10817669, rs10759734, rs6478105, rs10982396, and rs10733612). Finally, haplotype association investigations identified a strongly associated haplotype (P<8.8×10⁻⁵) consisting of these 6 SNPs and located in the direct vicinity of the TNFSF15 gene. In conclusion, we have identified within the SPA2 locus a haplotype strongly associated with predisposition to SpA which is located near to TNFSF15, one of the major candidate genes in this region.     

A Genome-Wide Association Study Uncovers a Genetic Locus Associated with Thoracic-to-Hip Ratio in Koreans

The thoracic-to-hip circumference ratio (THR) is an anthropometric marker recently described as a predictor of type 2 diabetes. In this study, we performed a genome-wide association study (GWAS) followed by confirmatory analyses to identify genetic markers associated with THR. A total of 7,240 Korean subjects (4,988 for the discovery stage and 2,252 for the confirmatory analyses) were recruited for this study, and genome-wide single nucleotide polymorphism (SNP) genotyping of the initial 4,988 individuals was performed using Affymetrix Human SNP array 5.0. Linear regression analysis was then performed to adjust for the effects of age, sex, and current diabetes medication status on the THR of the study subjects. In the initial discovery stage, there was a statistically nominal association between minor alleles of SNP markers on chromosomes 4, 8, 10, and 12, and THR changes (p < 5.0 × 10⁻⁶). The subsequent confirmatory analyses of these markers, however, only detected a significant association between two SNPs in the HECTD4 gene and decreased THRs. Notably, this association was detected in male (rs11066280: p = 1.14 × 10⁻²; rs2074356: p = 1.10 × 10⁻²), but not in female subjects. Meanwhile, the combined results from the two analyses (initial and confirmatory) indicated that minor alleles of these two intronic variants exhibited a significant genome-wide association with decreased THR in the male subjects (n = 3,155; rs11066280: effect size = −0.008624, p = 6.19 × 10⁻⁹; rs2074356: effect size = −0.008762, p = 1.89 × 10⁻⁸). Furthermore, minor alleles of these two SNPs exhibited protective effects on patients’ risks for developing type 2 diabetes. In conclusion, we have identified two genetic variations in HECTD4 that are associated with THR, particularly in men.