Genome-wide association analysis with selective genotyping identifies candidate loci for adult height at 8q21.13 and 15q22.33-q23 in Mongolians

We performed a genome-wide association study with 23,465 microsatellite markers to identify genes related to adult height. Selective genotyping was applied to extremely tall and extremely short individuals from the Khalkh-Mongolian population. Two loci, 8q21.13 and 15q22.33, which showed the strongest association with microsatellites were subjected to further analyses of SNPs in 782 tall and 773 short individuals. The most significant association was observed with SNP rs2220456 at 8q21.13 (P = 0.000016). In the LD block at 15q22.32, SNP rs8038652 located in intron 1 of IQCH was strongly associated (P = 0.0003), especially the AA genotype of the SNP under a recessive model was strongly associated with adult height (P = 0.000046).     

Genetic Polymorphisms Associated with Hearing Threshold Shift in Subjects during First Encounter with Occupational Impulse Noise

Noise-induced hearing loss (NIHL) is the most significant occupational health issue worldwide. We conducted a genome-wide association study to identify single-nucleotide polymorphisms (SNPs) associated with hearing threshold shift in young males undergoing their first encounter with occupational impulse noise. We report a significant association of SNP rs7598759 (p < 5 x 10^-7; p = 0.01 after permutation and correction; Odds Ratio = 12.75) in the gene coding for nucleolin, a multifunctional phosphoprotein involved in the control of senescence and protection against apoptosis. Interestingly, nucleolin has been shown to mediate the anti-apoptotic effect of HSP70, a protein found to prevent ototoxicity and whose polymorphisms have been associated with susceptibility to NIHL. Increase in nucleolin expression has also been associated with the prevention of apoptosis in cells undergoing oxidative stress, a well-known metabolic sequela of noise exposure. To assess the potential role of nucleolin in hearing loss, we tested down-regulation of nucleolin in cochlear sensory cells HEI-OC1 under oxidative stress conditions and report increased sensitivity to cisplatin, a chemotherapeutic drug with ototoxic side effects. Additional SNPs were found with suggestive association (p < 5 x 10^-4), of which 7 SNPs were located in genes previously reported to be related to NIHL and 43 of them were observed in 36 other genes previously not reported to be associated with NIHL. Taken together, our GWAS data and in vitro studies reported herein suggest that nucleolin is a potential candidate associated with NIHL in this population.     

Multiple Functional Risk Variants in a SMAD7 Enhancer Implicate a Colorectal Cancer Risk Haplotype

Genome-wide association studies (GWAS) of colorectal cancer (CRC) have led to the identification of a number of common variants associated with modest risk. Several risk variants map within the vicinity of TGFβ/BMP signaling pathway genes, including rs4939827 within an intron of SMAD7 at 18q21.1. A previous study implicated a novel SNP (novel 1 or rs58920878) as a functional variant within an enhancer element in SMAD7 intron 4. In this study, we show that four SNPs including novel 1 (rs6507874, rs6507875, rs8085824, and rs58920878) in linkage disequilibrium (LD) with the index SNP rs4939827 demonstrate allele-specific enhancer effects in a large, multi-component enhancer of SMAD7. All four SNPs demonstrate allele-specific protein binding to nuclear extracts of CRC cell lines. Furthermore, some of the risk-associated alleles correlate with increased expression of SMAD7 in normal colon tissues. Finally, we show that the enhancer is responsive to BMP4 stimulation. Taken together, we propose that the associated CRC risk at 18q21.1 is due to four functional variants that regulate SMAD7 expression and potentially perturb a BMP negative feedback loop in TGFβ/BMP signaling pathways.     

A genome-wide association study of nephrolithiasis in the Japanese population identifies novel susceptible Loci at 5q35.3, 7p14.3, and 13q14.1

Nephrolithiasis is a common nephrologic disorder with complex etiology. To identify the genetic factor(s) for nephrolithiasis, we conducted a three-stage genome-wide association study (GWAS) using a total of 5,892 nephrolithiasis cases and 17,809 controls of Japanese origin. Here we found three novel loci for nephrolithiasis: RGS14-SLC34A1-PFN3-F12 on 5q35.3 (rs11746443; P = 8.51×10⁻¹², odds ratio (OR) = 1.19), INMT-FAM188B-AQP1 on 7p14.3 (rs1000597; P = 2.16×10⁻¹⁴, OR = 1.22), and DGKH on 13q14.1 (rs4142110; P = 4.62×10⁻⁹, OR = 1.14). Subsequent analyses in 21,842 Japanese subjects revealed the association of SNP rs11746443 with the reduction of estimated glomerular filtration rate (eGFR) (P = 6.54×10⁻⁸), suggesting a crucial role for this variation in renal function. Our findings elucidated the significance of genetic variations for the pathogenesis of nephrolithiasis.     

Genome-Wide Association Study Implicates Chromosome 9q21.31 as a Susceptibility Locus for Asthma in Mexican Children

Many candidate genes have been studied for asthma, but replication has varied. Novel candidate genes have been identified for various complex diseases using genome-wide association studies (GWASs). We conducted a GWAS in 492 Mexican children with asthma, predominantly atopic by skin prick test, and their parents using the Illumina HumanHap 550 K BeadChip to identify novel genetic variation for childhood asthma. The 520,767 autosomal single nucleotide polymorphisms (SNPs) passing quality control were tested for association with childhood asthma using log-linear regression with a log-additive risk model. Eleven of the most significantly associated GWAS SNPs were tested for replication in an independent study of 177 Mexican case–parent trios with childhood-onset asthma and atopy using log-linear analysis. The chromosome 9q21.31 SNP rs2378383 (p = 7.10×10⁻⁶ in the GWAS), located upstream of transducin-like enhancer of split 4 (TLE4), gave a p-value of 0.03 and the same direction and magnitude of association in the replication study (combined p = 6.79×10⁻⁷). Ancestry analysis on chromosome 9q supported an inverse association between the rs2378383 minor allele (G) and childhood asthma. This work identifies chromosome 9q21.31 as a novel susceptibility locus for childhood asthma in Mexicans. Further, analysis of genome-wide expression data in 51 human tissues from the Novartis Research Foundation showed that median GWAS significance levels for SNPs in genes expressed in the lung differed most significantly from genes not expressed in the lung when compared to 50 other tissues, supporting the biological plausibility of our overall GWAS findings and the multigenic etiology of childhood asthma.     

Genetic Variants at 12p11 and 12q24 Are Associated with Breast Cancer Risk in a Chinese Population

Background

A recent genome-wide association study (GWAS) has identified three new breast cancer susceptibility loci at 12p11, 12q24 and 21q21 in populations of European descent. However, because of the genetic heterogeneity, it is largely unknown for the role of these loci in the breast cancer susceptibility in the populations of non-European descent.

Methodology/Principal Findings

Here, we genotyped three variants (rs10771399 at 12p11, rs1292011 at 12q24 and rs2823093 at 21q21) in an independent case–control study with a total of 1792 breast cancer cases and 1867 cancer-free controls in a Chinese population. We found that rs10771399 and rs1292011 were significantly associated with risk of breast cancer with per-allele odds ratios (ORs) of 0.85 (95% confidence interval (CI): 0.76–0.96; P = 0.010) and 0.84 (95% CI: 0.76–0.95; P = 4.50×10⁻³), respectively, which was consistent with those reported in populations of European descent. Similar effects were observed between ER/PR positive and negative breast cancer for both loci. However, we did not found significant association between rs2823093 and breast cancer risk (OR = 0.97, 95%CI = 0.76–1.24; P  = 0.795).

Conclusions/Significance

Our results indicate that genetic variants at 12p11 and 12q24 may also play an important role in breast cancer development in Chinese women.     

Dense‐map genome scan for dyslexia supports loci at 4q13, 16p12, 17q22; suggests novel locus at 7q36

Analysis of genetic linkage to dyslexia was performed using 133,165 array‐based SNPs genotyped in 718 persons from 101 dyslexia‐affected families. Results showed five linkage peaks with lod scores >2.3 (4q13.1, 7q36.1‐q36.2, 7q36.3, 16p12.1, and 17q22). Of these five regions, three have been previously implicated in dyslexia (4q13.1, 16p12.1, and 17q22), three have been implicated in attention‐deficit hyperactivity disorder (ADHD, which highly co‐occurs with dyslexia; 4q13.1, 7q36.3, 16p12.1) and four have been implicated in autism (a condition characterized by language deficits; 7q36.1‐q36.2, 7q36.3, 16p12.1, and 17q22). These results highlight the reproducibility of dyslexia linkage signals, even without formally significant lod scores, and suggest dyslexia predisposing genes with relatively major effects and locus heterogeneity. The largest lod score (2.80) occurred at 17q22 within the MSI2 gene, involved in neuronal stem cell lineage proliferation. Interestingly, the 4q13.1 linkage peak (lod 2.34) occurred immediately upstream of the LPHN3 gene, recently reported both linked and associated with ADHD. Separate analyses of larger pedigrees revealed lods >2.3 at 1–3 regions per family; one family showed strong linkage (lod 2.9) to a known dyslexia locus (18p11) not detected in our overall data, demonstrating the value of analyzing single large pedigrees. Association analysis identified no SNPs with genome‐wide significance, although a borderline significant SNP (P = 6 × 10^–7) occurred at 5q35.1 near FGF18, involved in laminar positioning of cortical neurons during development. We conclude that dyslexia genes with relatively major effects exist, are detectable by linkage analysis despite genetic heterogeneity, and show substantial overlapping predisposition with ADHD and autism.     

Genetic Variants at 6p21.1 and 7p15.3 Are Associated with Risk of Multiple Cancers in Han Chinese

Cancer susceptibility loci identified in reported genome-wide association studies (GWAS) are often tumor-specific; however, evidence of pleiotropy of some genes/loci has also been observed and biologically plausible. We hypothesized that there are important regions in the genome harboring genetic variants associated with risk of multiple types of cancer. In the current study, we attempted to map genetic variants that have consistent effects on risk of multiple cancers using our existing genome-wide scan data of lung cancer, noncardia gastric cancer, and esophageal squamous-cell carcinoma with overall 5,368 cases and 4,006 controls (GWAS stage), followed by a further evaluation in additional 9,001 cases with one of these cancer types and 11,436 controls (replication stage). Five variants satisfying the criteria of pleiotropy with p values from 1.10 × 10⁻⁸ to 8.96 × 10⁻⁶ for genome-wide scans of three cancer types were further evaluated in the replication stage. We found consistent associations of rs2494938 at 6p21.1 and rs2285947 at 7p15.3 with these three cancers in both GWAS and replication stages. In combined samples of GWAS and replication stages, the minor alleles of rs2494938 and rs2285947 were significantly associated with an increased risk of the cancers (odds ratio [OR] = 1.15, 95% confidence interval [CI], 1.10–1.19 and OR = 1.17, 95% CI, 1.12–1.21), with the p values being 1.20 × 10⁻¹² and 1.26 × 10⁻¹⁶, respectively, which are at a genome-wide significance level. Our findings highlight the potential importance of variants at 6p21.1 and 7p15.3 in the susceptibility to multiple cancers.     

Copy Number Variation Analysis on a Non-Hodgkin Lymphoma Case-Control Study Identifies an 11q25 Duplication Associated with Diffuse Large B-Cell Lymphoma

Recent GWAS have identified several susceptibility loci for NHL. Despite these successes, much of the heritable variation in NHL risk remains to be explained. Common copy-number variants are important genomic sources of variability, and hence a potential source to explain part of this missing heritability. In this study, we carried out a CNV analysis using GWAS data from 681 NHL cases and 749 controls to explore the relationship between common structural variation and lymphoma susceptibility. Here we found a novel association with diffuse large B-cell lymphoma (DLBCL) risk involving a partial duplication of the C-terminus region of the LOC283177 long non-coding RNA that was further confirmed by quantitative PCR. For chronic lymphocytic leukemia/small lymphocytic lymphoma (CLL/SLL), known somatic deletions were identified on chromosomes 13q14, 11q22-23, 14q32 and 22q11.22. Our study shows that GWAS data can be used to identify germline CNVs associated with disease risk for DLBCL and somatic CNVs for CLL/SLL.     

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.     

Common Variation at 1q24.1 (ALDH9A1) Is a Potential Risk Factor for Renal Cancer

So far six susceptibility loci for renal cell carcinoma (RCC) have been discovered by genome-wide association studies (GWAS). To identify additional RCC common risk loci, we performed a meta-analysis of published GWAS (totalling 2,215 cases and 8,566 controls of Western-European background) with imputation using 1000 Genomes Project and UK10K Project data as reference panels and followed up the most significant association signals [22 single nucleotide polymorphisms (SNPs) and 3 indels in eight genomic regions] in 383 cases and 2,189 controls from The Cancer Genome Atlas (TCGA). A combined analysis identified a promising susceptibility locus mapping to 1q24.1 marked by the imputed SNP rs3845536 (P _combined =2.30x10^-8). Specifically, the signal maps to intron 4 of the ALDH9A1 gene (aldehyde dehydrogenase 9 family, member A1). We further evaluated this potential signal in 2,461 cases and 5,081 controls from the International Agency for Research on Cancer (IARC) GWAS of RCC cases and controls from multiple European regions. In contrast to earlier findings no association was shown in the IARC series (P=0.94; P _combined =2.73x10^-5). While variation at 1q24.1 represents a potential risk locus for RCC, future replication analyses are required to substantiate our observation.     

Identification of novel loci influencing refractive error in East Asian populations using an extreme phenotype design

The global “myopia boom” has raised significant international concerns. Despite a higher myopia prevalence in Asia, previous large-scale genome-wide association studies (GWASs) were mostly based on European descendants. Here, we report a GWAS of spherical equivalent (SE) in 1852 Chinese Han individuals with extreme SE from Guangzhou (631 < ?6.00D and 574 > 0.00D) and Wenzhou (593 < ?6.00D and 54 > ?1.75D), followed by a replication study in two independent cohorts with totaling 3538 East Asian individuals. The discovery GWAS and meta-analysis identify three novel loci, which show genome-wide significant associations with SE, including 1q25.2 FAM163A, 10p11.22 NRP1/PRAD3, and 10p11.21 ANKRD30A/MTRNR2L7, together explaining 3.34% of SE variance. 10p11.21 is successfully replicated. The allele frequencies of all three loci show significant differences between major continental groups (P < 0.001). The SE reducing (more myopic) allele of rs10913877 (1q25.2 FAM163A) demonstrates the highest frequency in East Asians and much lower frequencies in Europeans and Africans (EAS = 0.60, EUR = 0.20, and AFR = 0.18). The gene-based analysis additionally identifies three novel genes associated with SE, including EI24, LHX5, and ARPP19. These results provide new insights into myopia pathogenesis and indicate the role of genetic heterogeneity in myopia epidemiology among different ethnicities.     

Association analysis of single nucleotide polymorphisms at five loci: comparison between atopic dermatitis and asthma in the Chinese Han population

Atopic diseases, such as atopic dermatitis (AD) and asthma, are closely related to clinical phenotypes with hypersensitivity, and often share some similar genetic and pathogenic bases. Our recent GWAS identified three susceptibility gene/loci FLG (rs11204971 and rs3126085), 5q22.1 (rs10067777, rs7701890, rs13360927 and rs13361382) and 20q13.33 (rs6010620) to AD. The effect of these AD associated polymorphisms in asthma is so far unknown. To investigate whether AD relevant genetic variants is identical to asthma and reveal the differences in genetic factors between AD and asthma in Chinese Han population, seven AD associated single nucleotide polymorphisms (SNPs) as well as 3 other SNPs (rs7936562 and rs7124842 at 11q13.5 and rs4982958 at 14q11.2) from our previous AD GWAS were genotyped in 463 asthma patients and 985 controls using Sequenom MassArray system. We found rs4982958 at 14q11.2 was significantly associated with asthma (P = 3.04×10⁻⁴, OR = 0.73). We also detected one significant risk haplotype GGGA from the 4 SNPs (rs10067777, rs7701890, rs13360927 and rs13361382) at 5q22.1 in AD cases (P _correction = 3.60×10⁻¹⁰, OR = 1.26), and the haplotype was suggestive of risk in asthma cases in this study (P = 0.014, P _correction = 0.084, OR = 1.38). These SNPs (rs11204971, rs3126085, rs7936562, rs712484 and rs6010620) at AD susceptibility genes/loci FLG, 11q13.5 and 20q13.33 were not associated with asthma in this study. Our results further comfirmed that 14q11.2 was an important candidate locus for asthma and demonstrated that 5q22.1 might be shared by AD and asthma in Chinese Han 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.     

Identification of potential functional variants and genes at 18q21.1 associated with the carcinogenesis of colorectal cancer

Genome-wide association studies (GWAS) have identified more than 160 susceptibility loci for colorectal cancer (CRC). The effects of these variants, particularly their mechanisms, however, remain unclear. In this study, a comprehensive functional annotation of CRC-related GWAS signals was firstly conducted to identify the potential causal variants. We found that the SNP rs7229639 in intron 3 of SMAD7 at 18q21.1 might serve as a putative functional variant in CRC. The SNP rs7229639 is located in a region with evidence of regulatory potential. Dual-luciferase reporter assays revealed that three other SNPs (rs77544449, rs60385309 and rs72917785), in strong linkage disequilibrium (LD) with rs7229639, exhibited allele-specific enhancer activity, of which one of the target genes may conceivably be LIPG, as suggested by eQTL association data and Hi-C data. We also verified that LIPG promoted malignancy of CRC cells in vitro, with supporting clinical data indicating that LIPG is upregulated and correlated with a poor prognosis in CRC. Finally, pitavastatin was observed to exhibit an anti-CRC activity and modest inhibition of LIPG mRNA levels. Collectively, our data suggest that these functional variants at 18q21.1 are involved in the pathogenesis of CRC by modulating enhancer activity, and possibly LIPG expression, thus indicating a promising therapeutic target for CRC. The results of functional annotation in our investigation could also serve as an inventory for CRC susceptibility SNPs and offer guides for post-GWAS downstream functional studies.     

Association of chromosome 8q24 variants with prostate cancer risk in the Siberian region of Russia and meta-analysis

Compelling evidence demonstrates the importance of chromosome 8q24 as a locus of susceptibility to prostate cancer. In this work, the association of common 8q24 variants, rs6983267 and rs1447295, with a sporadic risk of prostate cancer was analyzed in the Russian population of Siberia. For this purpose, the above polymorphisms were genotyped in 393 cases and 384 control individuals. The A allele of rs1447295 was significantly associated with prostate cancer risk (OR[CI 95%] = 1.74 [1.26–2.4], p = 7.8 × 10⁻⁴). The common G-A haplotype of rs6983267-rs1447295 also showed association with prostate cancer risk in Russians (OR[CI 95%] = 2.03 [1.1–3.75], p = 0.02). A meta-analysis combining our data with previously published results was performed to better evaluate the association between the SNPs studied and prostate cancer risk; its results strongly supported the association for both loci (p < 10⁻⁶). Thus, our study has confirmed the association of chromosome 8q24 with a risk of prostate cancer.     

Associations between common variants in <Emphasis Type="Italic">GC</Emphasis> and <Emphasis Type="Italic">DHCR7</Emphasis>/<Emphasis Type="Italic">NADSYN1</Emphasis> and vitamin D concentration in Chinese Hans

Recent studies have identified common variants in or near GC, CYP2R1 and NADSYN1/DHCR7 to be associated with 25-hydroxyvitamin D [25(OH)D] levels in European populations. We aimed to examine whether these variants also influence 25(OH)D levels in Chinese. Seven common variants were successfully genotyped and tested for associations with plasma 25(OH)D levels in a population-based cohort of 3,210 Chinese Hans from Beijing and Shanghai. Six common variants at GC (rs4588, rs7041, rs2282679 and rs1155563) and NADSYN1/DHCR7 (rs3829251 and rs1790349) loci were all significantly associated with lower plasma 25(OH)D levels (−0.036 ≤ β ≤ −0.076 per risk-allele, P ≤ 5.7 × 10⁻⁵), while CYP2R1-rs2060793 showed a trend toward association with 25(OH)D levels in the Shanghai subpopulation (P = 0.08), but not in the Beijing subpopulation (P = 0.82). Haplotype-based association analyses of the four GC variants showed that only the haplotype that contained all risk-alleles (TACC) was significantly associated with lower plasma 25(OH)D levels (β = −0.085, P = 2.3 × 10⁻⁹), while the haplotype containing the risk-alleles of rs4588 and rs2282679 (TATC) was marginally associated with lower 25(OH)D levels (β = −0.054, P = 0.0562) when compared with GCTA haplotype carrying the four protective alleles. Most notably, conditional analyses showed that only GC-rs4588 and GC-rs2282679 (r ² = 0.97) remained significantly associated with 25(OH)D concentrations (P ≤ 1.9 × 10⁻⁵) after adjusting for the other two SNPs in GC. In conclusion, GC and NADSYN1/DHCR7 loci individually and collectively contribute to variation in plasma vitamin D levels in Chinese Hans.     

Lung eQTLs to Help Reveal the Molecular Underpinnings of Asthma

Genome-wide association studies (GWAS) have identified loci reproducibly associated with pulmonary diseases; however, the molecular mechanism underlying these associations are largely unknown. The objectives of this study were to discover genetic variants affecting gene expression in human lung tissue, to refine susceptibility loci for asthma identified in GWAS studies, and to use the genetics of gene expression and network analyses to find key molecular drivers of asthma. We performed a genome-wide search for expression quantitative trait loci (eQTL) in 1,111 human lung samples. The lung eQTL dataset was then used to inform asthma genetic studies reported in the literature. The top ranked lung eQTLs were integrated with the GWAS on asthma reported by the GABRIEL consortium to generate a Bayesian gene expression network for discovery of novel molecular pathways underpinning asthma. We detected 17,178 cis- and 593 trans- lung eQTLs, which can be used to explore the functional consequences of loci associated with lung diseases and traits. Some strong eQTLs are also asthma susceptibility loci. For example, rs3859192 on chr17q21 is robustly associated with the mRNA levels of GSDMA (P = 3.55×10⁻¹⁵¹). The genetic-gene expression network identified the SOCS3 pathway as one of the key drivers of asthma. The eQTLs and gene networks identified in this study are powerful tools for elucidating the causal mechanisms underlying pulmonary disease. This data resource offers much-needed support to pinpoint the causal genes and characterize the molecular function of gene variants associated with lung diseases.     

Genetic variants at chromosomes 2q35, 5p12, 6q25.1, 10q26.13, and 16q12.1 influence the risk of breast cancer in men

Male breast cancer accounts for approximately 1% of all breast cancer. To date, risk factors for male breast cancer are poorly defined, but certain risk factors and genetic features appear common to both male and female breast cancer. Genome-wide association studies (GWAS) have recently identified common single nucleotide polymorphisms (SNPs) that influence female breast cancer risk; 12 of these have been independently replicated. To examine if these variants contribute to male breast cancer risk, we genotyped 433 male breast cancer cases and 1,569 controls. Five SNPs showed a statistically significant association with male breast cancer: rs13387042 (2q35) (odds ratio (OR)  = 1.30, p = 7.98×10⁻⁴), rs10941679 (5p12) (OR = 1.26, p = 0.007), rs9383938 (6q25.1) (OR = 1.39, p = 0.004), rs2981579 (FGFR2) (OR = 1.18, p = 0.03), and rs3803662 (TOX3) (OR = 1.48, p = 4.04×10⁻⁶). Comparing the ORs for male breast cancer with the published ORs for female breast cancer, three SNPs—rs13387042 (2q35), rs3803662 (TOX3), and rs6504950 (COX11)—showed significant differences in ORs (p<0.05) between sexes. Breast cancer is a heterogeneous disease; the relative risks associated with loci identified to date show subtype and, based on these data, gender specificity. Additional studies of well-defined patient subgroups could provide further insight into the biological basis of breast cancer development.