Replication of the 4p16 Susceptibility Locus in Congenital Heart Disease in Han Chinese Populations

Congenital heart disease (CHD) is the most common form of congenital human birth anomalies and a leading cause of perinatal and infant mortality. Some studies including our published genome-wide association study (GWAS) of CHD have indicated that genetic variants may contribute to the risk of CHD. Recently, Cordell et al. published a GWAS of multiple CHD phenotypes in European Caucasians and identified 3 susceptibility loci (rs870142, rs16835979 and rs6824295) for ostium secundum atrial septal defect (ASD) at chromosome 4p16. However, whether these loci at 4p16 confer the predisposition to CHD in Chinese population is unclear. In the current study, we first analyzed the associations between these 3 single nucleotide polymorphisms (SNPs) at 4p16 and CHD risk by using our existing genome-wide scan data and found all of the 3 SNPs showed significant associations with ASD in the same direction as that observed in Cordell’s study, but not with other subtypes- ventricular septal defect (VSD) and ASD combined VSD. As these 3 SNPs were in high linkage disequilibrium (LD) in Chinese population, we selected one SNP with the lowest P value in our GWAS scan (rs16835979) to perform a replication study with additional 1,709 CHD cases with multiple phenotypes and 1,962 controls. The significant association was also observed only within the ASD subgroup, which was heterogeneous from other disease groups. In combined GWAS and replication samples, the minor allele of rs16835979 remained significant association with the risk of ASD (OR = 1.22, 95% CI = 1.08–1.38, P = 0.001). Our findings suggest that susceptibility loci of ASD identified from Cordell’s European GWAS are generalizable to Chinese population, and such investigation may provide new insights into the roles of genetic variants in the etiology of different CHD phenotypes.     

Genome-Wide Association Study for Levels of Total Serum IgE Identifies HLA-C in a Japanese Population

Most of the previously reported loci for total immunoglobulin E (IgE) levels are related to Th2 cell-dependent pathways. We undertook a genome-wide association study (GWAS) to identify genetic loci responsible for IgE regulation. A total of 479,940 single nucleotide polymorphisms (SNPs) were tested for association with total serum IgE levels in 1180 Japanese adults. Fine-mapping with SNP imputation demonstrated 6 candidate regions: the PYHIN1/IFI16, MHC classes I and II, LEMD2, GRAMD1B, and chr13∶60576338 regions. Replication of these candidate loci in each region was assessed in 2 independent Japanese cohorts (n = 1110 and 1364, respectively). SNP rs3130941 in the HLA-C region was consistently associated with total IgE levels in 3 independent populations, and the meta-analysis yielded genome-wide significance (P = 1.07×10⁻¹⁰). Using our GWAS results, we also assessed the reproducibility of previously reported gene associations with total IgE levels. Nine of 32 candidate genes identified by a literature search were associated with total IgE levels after correction for multiple testing. Our findings demonstrate that SNPs in the HLA-C region are strongly associated with total serum IgE levels in the Japanese population and that some of the previously reported genetic associations are replicated across ethnic groups.     

Genome-Wide Association Analyses for Fatty Acid Composition in Porcine Muscle and Abdominal Fat Tissues

Fatty acid composition is an important phenotypic trait in pigs as it affects nutritional, technical and sensory quality of pork. Here, we reported a genome-wide association study (GWAS) for fatty acid composition in the longissimus muscle and abdominal fat tissues of 591 White Duroc×Erhualian F₂ animals and in muscle samples of 282 Chinese Sutai pigs. A total of 46 loci surpassing the suggestive significance level were identified on 15 pig chromosomes (SSC) for 12 fatty acids, revealing the complex genetic architecture of fatty acid composition in pigs. Of the 46 loci, 15 on SSC5, 7, 14 and 16 reached the genome-wide significance level. The two most significant SNPs were ss131535508 (P = 2.48×10⁻²⁵) at 41.39 Mb on SSC16 for C20∶0 in abdominal fat and ss478935891 (P = 3.29×10⁻¹³) at 121.31 Mb on SSC14 for muscle C18∶0. A meta-analysis of GWAS identified 4 novel loci and enhanced the association strength at 6 loci compared to those evidenced in a single population, suggesting the presence of common underlying variants. The longissimus muscle and abdominal fat showed consistent association profiles at most of the identified loci and distinct association signals at several loci. All loci have specific effects on fatty acid composition, except for two loci on SSC4 and SSC7 affecting multiple fatness traits. Several promising candidate genes were found in the neighboring regions of the lead SNPs at the genome-wide significant loci, such as SCD for C18∶0 and C16∶1 on SSC14 and ELOVL7 for C20∶0 on SSC16. The findings provide insights into the molecular basis of fatty acid composition in pigs, and would benefit the final identification of the underlying mutations.     

A Single Nucleotide Polymorphism within DUSP9 Is Associated with Susceptibility to Type 2 Diabetes in a Japanese Population

Aims

The DUSP9 locus on chromosome X was identified as a susceptibility locus for type 2 diabetes in a meta-analysis of European genome-wide association studies (GWAS), and GWAS in South Asian populations identified 6 additional single nucleotide polymorphism (SNP) loci for type 2 diabetes. However, the association of these loci with type 2 diabetes have not been examined in the Japanese. We performed a replication study to investigate the association of these 7 susceptibility loci with type 2 diabetes in the Japanese population.

Methods

We genotyped 11,319 Japanese participants (8,318 with type 2 diabetes and 3,001 controls) for each of the 7 SNPs–rs5945326 near DUSP9, rs3923113 near GRB14, rs16861329 in ST6GAL1, rs1802295 in VPS26A, rs7178572 in HMG20A, rs2028299 near AP3S2, and rs4812829 in HNF4A–and examined the association of each of these 7 SNPs with type 2 diabetes by using logistic regression analysis.

Results

All SNPs had the same direction of effect (odds ratio [OR]>1.0) as in the original reports. One SNP, rs5945326 near DUSP9, was significantly associated with type 2 diabetes at a genome-wide significance level (p = 2.21×10⁻⁸; OR 1.39, 95% confidence interval [CI]: 1.24−1.56). The 6 SNPs derived from South Asian GWAS were not significantly associated with type 2 diabetes in the Japanese population by themselves (p≥0.007). However, a genetic risk score constructed from 6 South Asian GWAS derived SNPs was significantly associated with Japanese type 2 diabetes (p = 8.69×10⁻⁴, OR  = 1.06. 95% CI; 1.03−1.10).

Conclusions/interpretation

These results indicate that the DUSP9 locus is a common susceptibility locus for type 2 diabetes across different ethnicities, and 6 loci identified in South Asian GWAS also have significant effect on susceptibility to Japanese type 2 diabetes.     

Association of Genetic Variants with Primary Angle Closure Glaucoma in Two Different Populations

Purpose

A recent large genome-wide association study (GWAS) identified multiple variants associated with primary angle-closure glaucoma (PACG). The present study investigated the role of these variants in two cohorts with PACG recruited from Australia and Nepal.

Method

Patients with PACG and appropriate controls were recruited from eye clinics in Australia (n = 232 cases and n = 288 controls) and Nepal (n = 106 cases and 204 controls). Single nucleotide polymorphisms (SNPs) rs3753841 (COL11A1), rs1015213 (located between PCMTD1 and ST18), rs11024102 (PLEKHA7), and rs3788317 (TXNRD2) were selected and genotyped on the Sequenom. Analyses were conducted using PLINK and METAL.

Results

After adjustment for age and sex, SNP rs3753841 was found to be significantly associated with PACG in the Australian cohort (p = 0.017; OR = 1.34). SNPs rs1015213 (p = 0.014; OR 2.35) and rs11024102 (p = 0.039; OR 1.43) were significantly associated with the disease development in the Nepalese cohort. None of these SNPs survived Bonferroni correction (p = 0.05/4 = 0.013). However, in the combined analysis, of both cohorts, rs3753841 and rs1015213 showed significant association with p-values of 0.009 and 0.004, respectively both surviving Bonferroni correction. SNP rs11024102 showed suggestive association with PACG (p-value 0.035) and no association was found with rs3788317.

Conclusion

The present results support the initial GWAS findings, and confirm the SNP’s contribution to PACG. This is the first study to investigate these loci in both Australian Caucasian and Nepalese populations.     

Association of CHRDL1 Mutations and Variants with X-linked Megalocornea,Neuh?user Syndrome and Central Corneal Thickness

We describe novel CHRDL1 mutations in ten families with X-linked megalocornea (MGC1). Our mutation-positive cohort enabled us to establish ultrasonography as a reliable clinical diagnostic tool to distinguish between MGC1 and primary congenital glaucoma (PCG). Megalocornea is also a feature of Neuhäuser or megalocornea-mental retardation (MMR) syndrome, a rare condition of unknown etiology. In a male patient diagnosed with MMR, we performed targeted and whole exome sequencing (WES) and identified a novel missense mutation in CHRDL1 that accounts for his MGC1 phenotype but not his non-ocular features. This finding suggests that MMR syndrome, in some cases, may be di- or multigenic. MGC1 patients have reduced central corneal thickness (CCT); however no X-linked loci have been associated with CCT, possibly because the majority of genome-wide association studies (GWAS) overlook the X-chromosome. We therefore explored whether variants on the X-chromosome are associated with CCT. We found rs149956316, in intron 6 of CHRDL1, to be the most significantly associated single nucleotide polymorphism (SNP) (p = 6.81×10⁻⁶) on the X-chromosome. However, this association was not replicated in a smaller subset of whole genome sequenced samples. This study highlights the importance of including X-chromosome SNP data in GWAS to identify potential loci associated with quantitative traits or disease risk.     

Common Genetic Variation Near the Phospholamban Gene Is Associated with Cardiac Repolarisation: Meta-Analysis of Three Genome-Wide Association Studies

To identify loci affecting the electrocardiographic QT interval, a measure of cardiac repolarisation associated with risk of ventricular arrhythmias and sudden cardiac death, we conducted a meta-analysis of three genome-wide association studies (GWAS) including 3,558 subjects from the TwinsUK and BRIGHT cohorts in the UK and the DCCT/EDIC cohort from North America. Five loci were significantly associated with QT interval at P<1×10⁻⁶. To validate these findings we performed an in silico comparison with data from two QT consortia: QTSCD (n = 15,842) and QTGEN (n = 13,685). Analysis confirmed the association between common variants near NOS1AP (P = 1.4×10⁻⁸³) and the phospholamban (PLN) gene (P = 1.9×10⁻²⁹). The most associated SNP near NOS1AP (rs12143842) explains 0.82% variance; the SNP near PLN (rs11153730) explains 0.74% variance of QT interval duration. We found no evidence for interaction between these two SNPs (P = 0.99). PLN is a key regulator of cardiac diastolic function and is involved in regulating intracellular calcium cycling, it has only recently been identified as a susceptibility locus for QT interval. These data offer further mechanistic insights into genetic influence on the QT interval which may predispose to life threatening arrhythmias and sudden cardiac death.     

Genome-Wide Association of Body Fat Distribution in African Ancestry Populations Suggests New Loci

Central obesity, measured by waist circumference (WC) or waist-hip ratio (WHR), is a marker of body fat distribution. Although obesity disproportionately affects minority populations, few studies have conducted genome-wide association study (GWAS) of fat distribution among those of predominantly African ancestry (AA). We performed GWAS of WC and WHR, adjusted and unadjusted for BMI, in up to 33,591 and 27,350 AA individuals, respectively. We identified loci associated with fat distribution in AA individuals using meta-analyses of GWA results for WC and WHR (stage 1). Overall, 25 SNPs with single genomic control (GC)-corrected p-values<5.0×10⁻⁶ were followed-up (stage 2) in AA with WC and with WHR. Additionally, we interrogated genomic regions of previously identified European ancestry (EA) WHR loci among AA. In joint analysis of association results including both Stage 1 and 2 cohorts, 2 SNPs demonstrated association, rs2075064 at LHX2, p = 2.24×10⁻⁸ for WC-adjusted-for-BMI, and rs6931262 at RREB1, p = 2.48×10⁻⁸ for WHR-adjusted-for-BMI. However, neither signal was genome-wide significant after double GC-correction (LHX2: p = 6.5×10⁻⁸; RREB1: p = 5.7×10⁻⁸). Six of fourteen previously reported loci for waist in EA populations were significant (p<0.05 divided by the number of independent SNPs within the region) in AA studied here (TBX15-WARS2, GRB14, ADAMTS9, LY86, RSPO3, ITPR2-SSPN). Further, we observed associations with metabolic traits: rs13389219 at GRB14 associated with HDL-cholesterol, triglycerides, and fasting insulin, and rs13060013 at ADAMTS9 with HDL-cholesterol and fasting insulin. Finally, we observed nominal evidence for sexual dimorphism, with stronger results in AA women at the GRB14 locus (p for interaction = 0.02). In conclusion, we identified two suggestive loci associated with fat distribution in AA populations in addition to confirming 6 loci previously identified in populations of EA. These findings reinforce the concept that there are fat distribution loci that are independent of generalized adiposity.     

Sequence Analysis of Six Blood Pressure Candidate Regions in 4,178 Individuals: The Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Targeted Sequencing Study

Background

Genome-wide association studies (GWAS) identified multiple loci for blood pressure (BP) and hypertension. Six genes – ATP2B1, CACNB2, CYP17A1, JAG1, PLEKHA7, and SH2B3 – were evaluated for sequence variation with large effects on systolic blood pressure (SBP), diastolic blood pressure (DBP), pulse pressure (PP), and mean arterial pressure (MAP).

Methods and Results

Targeted genomic sequence was determined in 4,178 European ancestry participants from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium. Common variants (≥50 minor allele copies) were evaluated individually and rare variants (minor allele frequency, MAF≤1%) were aggregated by locus. 464 common variants were identified across the 6 genes. An upstream CYP17A1 variant, rs11191416 (MAF = 0.09), was the most significant association for SBP (P = 0.0005); however the association was attenuated (P = 0.0469) after conditioning on the GWAS index single nucleotide polymorphism (SNP). A PLEKHA7 intronic variant was the strongest DBP association (rs12806040, MAF = 0.007, P = 0.0006) and was not in LD (r² = 0.01) with the GWAS SNP. A CACNB2 intronic SNP, rs1571787, was the most significant association with PP (MAF = 0.27, P = 0.0003), but was not independent from the GWAS SNP (r² = 0.34). Three variants (rs6163 and rs743572 in the CYP17A1 region and rs112467382 in PLEKHA7) were associated with BP traits (P<0.001). Rare variation, aggregately assessed in the 6 regions, was not significantly associated with BP measures.

Conclusion

Six targeted gene regions, previously identified by GWAS, did not harbor novel variation with large effects on BP in this sample.     

A Divergent Artiodactyl MYADM-like Repeat Is Associated with Erythrocyte Traits and Weight of Lamb Weaned in Domestic Sheep

A genome-wide association study (GWAS) was performed to investigate seven red blood cell (RBC) phenotypes in over 500 domestic sheep (Ovis aries) from three breeds (Columbia, Polypay, and Rambouillet). A single nucleotide polymorphism (SNP) showed genome-wide significant association with increased mean corpuscular hemoglobin concentration (MCHC, P = 6.2×10⁻¹⁴) and genome-wide suggestive association with decreased mean corpuscular volume (MCV, P = 2.5×10⁻⁶). The ovine HapMap project found the same genomic region and the same peak SNP has been under extreme historical selective pressure, demonstrating the importance of this region for survival, reproduction, and/or artificially selected traits. We observed a large (>50 kb) variant haplotype sequence containing a full-length divergent artiodactyl MYADM-like repeat in strong linkage disequilibrium with the associated SNP. MYADM gene family members play roles in membrane organization and formation in myeloid cells. However, to our knowledge, no member of the MYADM gene family has been identified in development of morphologically variant RBCs. The specific RBC differences may be indicative of alterations in morphology. Additionally, erythrocytes with altered morphological structure often exhibit increased structural fragility, leading to increased RBC turnover and energy expenditure. The divergent artiodactyl MYADM-like repeat was also associated with increased ewe lifetime kilograms of lamb weaned (P = 2×10⁻⁴). This suggests selection for normal RBCs might increase lamb weights, although further validation is required before implementation in marker-assisted selection. These results provide clues to explain the strong selection on the artiodactyl MYADM-like repeat locus in sheep, and suggest MYADM family members may be important for RBC morphology in other mammals.     

Oligoclonal Band Status in Scandinavian Multiple Sclerosis Patients Is Associated with Specific Genetic Risk Alleles

The presence of oligoclonal bands (OCB) in cerebrospinal fluid (CSF) is a typical finding in multiple sclerosis (MS). We applied data from Norwegian, Swedish and Danish (i.e. Scandinavian) MS patients from a genome-wide association study (GWAS) to search for genetic differences in MS relating to OCB status. GWAS data was compared in 1367 OCB positive and 161 OCB negative Scandinavian MS patients, and nine of the most associated SNPs were genotyped for replication in 3403 Scandinavian MS patients. HLA-DRB1 genotypes were analyzed in a subset of the OCB positive (n = 2781) and OCB negative (n = 292) MS patients and compared to 890 healthy controls. Results from the genome-wide analyses showed that single nucleotide polymorphisms (SNPs) from the HLA complex and six other loci were associated to OCB status. In SNPs selected for replication, combined analyses showed genome-wide significant association for two SNPs in the HLA complex; rs3129871 (p = 5.7×10⁻¹⁵) and rs3817963 (p = 5.7×10⁻¹⁰) correlating with the HLA-DRB1*15 and the HLA-DRB1*04 alleles, respectively. We also found suggestive association to one SNP in the Calsyntenin-2 gene (p = 8.83×10⁻⁷). In HLA-DRB1 analyses HLA-DRB1*15∶01 was a stronger risk factor for OCB positive than OCB negative MS, whereas HLA-DRB1*04∶04 was associated with increased risk of OCB negative MS and reduced risk of OCB positive MS. Protective effects of HLA-DRB1*01∶01 and HLA-DRB1*07∶01 were detected in both groups. The groups were different with regard to age at onset (AAO), MS outcome measures and gender. This study confirms both shared and distinct genetic risk for MS subtypes in the Scandinavian population defined by OCB status and indicates different clinical characteristics between the groups. This suggests differences in disease mechanisms between OCB negative and OCB positive MS with implications for patient management, which need to be further studied.     

Replication Study Confirms Link between TSPAN18 Mutation and Schizophrenia in Han Chinese

Schizophrenia (SCZ) is a severe psychiatric disorder associated with many different risk factors, both genetic and environmental. A recent genome-wide association study (GWAS) of Han Chinese identified three single-nucleotide polymorphisms (SNPs rs11038167, rs11038172, and rs835784) in the tetraspanins gene TSPAN18 as possible susceptibility loci for schizophrenia. Hoping to validate these findings, we conducted a case-control study of Han Chinese with 1093 schizophrenia cases and 1022 healthy controls. Using the LDR-PCR method to genotype polymorphisms in TSPAN18, we found no significant differences (P>0.05) between patients and controls in either the allele or genotype frequency of the SNPs rs11038167 and rs11038172. We did find, however, that the frequency of the ‘A’ allele of SNP rs835784 is significantly higher in patients than in controls. We further observed a significant association (OR  = 1.197, 95%CI  = 1.047–1.369) between risk for SCZ and this ‘A’ allele. These results confirm the significant association, in Han Chinese populations, of increased SCZ risk and the variant of the TSPAN18 gene containing the ‘A’ allele of SNP rs835784.     

Blood cis-eQTL Analysis Fails to Identify Novel Association Signals among Sub-Threshold Candidates from Genome-Wide Association Studies in Restless Legs Syndrome

Restless legs syndrome (RLS) is a common neurologic disorder characterized by nightly dysesthesias affecting the legs primarily during periods of rest and relieved by movement. RLS is a complex genetic disease and susceptibility factors in six genomic regions have been identified by means of genome-wide association studies (GWAS). For some complex genetic traits, expression quantitative trait loci (eQTLs) are enriched among trait-associated single nucleotide polymorphisms (SNPs). With the aim of identifying new genetic susceptibility factors for RLS, we assessed the 332 best-associated SNPs from the genome-wide phase of the to date largest RLS GWAS for cis-eQTL effects in peripheral blood from individuals of European descent. In 740 individuals belonging to the KORA general population cohort, 52 cis-eQTLs with p_nominal<10⁻³ were identified, while in 976 individuals belonging to the SHIP-TREND general population study 53 cis-eQTLs with p_nominal<10⁻³ were present. 23 of these cis-eQTLs overlapped between the two cohorts. Subsequently, the twelve of the 23 cis-eQTL SNPs, which were not located at an already published RLS-associated locus, were tested for association in 2449 RLS cases and 1462 controls. The top SNP, located in the DET1 gene, was nominally significant (p<0.05) but did not withstand correction for multiple testing (p = 0.42). Although a similar approach has been used successfully with regard to other complex diseases, we were unable to identify new genetic susceptibility factor for RLS by adding this novel level of functional assessment to RLS GWAS data.     

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.     

Evaluation of Replication of Variants Associated with Genetic Risk of Otitis Media

The first Genome Wide Association Study (GWAS) of otitis media (OM) found evidence of association in the Western Australian Pregnancy Cohort (Raine) study, but lacked replication in an independent OM population. The aim of this study was to investigate association at these loci in our family-based sample of chronic otitis media with effusion and recurrent otitis media (COME/ROM). Autosomal SNPs were selected from the Raine OM GWAS results. SNPs from the Raine cohort GWAS genotyped in our GWAS of COME/ROM had P-values ranging from P = 0.06–0.80. After removal of SNPs previously genotyped in our GWAS of COME/ROM (N = 21) and those that failed Fluidigm assay design (N = 1), 26 SNPs were successfully genotyped in 716 individuals from our COME/ROM family population. None of the SNP associations replicated in our family-based population (unadjusted P = 0.03–0.93). Replication in an independent sample would confirm that these represent novel OM loci, and that further investigation is warranted.     

Large-Scale Genome-Wide Association Studies and Meta-Analyses of Longitudinal Change in Adult Lung Function

Background

Genome-wide association studies (GWAS) have identified numerous loci influencing cross-sectional lung function, but less is known about genes influencing longitudinal change in lung function.

Methods

We performed GWAS of the rate of change in forced expiratory volume in the first second (FEV₁) in 14 longitudinal, population-based cohort studies comprising 27,249 adults of European ancestry using linear mixed effects model and combined cohort-specific results using fixed effect meta-analysis to identify novel genetic loci associated with longitudinal change in lung function. Gene expression analyses were subsequently performed for identified genetic loci. As a secondary aim, we estimated the mean rate of decline in FEV₁ by smoking pattern, irrespective of genotypes, across these 14 studies using meta-analysis.

Results

The overall meta-analysis produced suggestive evidence for association at the novel IL16/STARD5/TMC3 locus on chromosome 15 (P  =  5.71 × 10^-7). In addition, meta-analysis using the five cohorts with ≥3 FEV₁ measurements per participant identified the novel ME3 locus on chromosome 11 (P  =  2.18 × 10^-8) at genome-wide significance. Neither locus was associated with FEV₁ decline in two additional cohort studies. We confirmed gene expression of IL16, STARD5, and ME3 in multiple lung tissues. Publicly available microarray data confirmed differential expression of all three genes in lung samples from COPD patients compared with controls. Irrespective of genotypes, the combined estimate for FEV₁ decline was 26.9, 29.2 and 35.7 mL/year in never, former, and persistent smokers, respectively.

Conclusions

In this large-scale GWAS, we identified two novel genetic loci in association with the rate of change in FEV₁ that harbor candidate genes with biologically plausible functional links to lung function.     

Genome-Wide Scan Identifies Variant in TNFSF13 Associated with Serum IgM in a Healthy Chinese Male Population

IgM provides a first line of defense during microbial infections. Serum IgM levels are detected routinely in clinical practice. And IgM is a genetically complex trait. We conducted a two-stage genome-wide association study (GWAS) to identify genetic variants affecting serum IgM levels in a Chinese population of 3495, including 1999 unrelated subjects in the first stage and 1496 independent individuals in the second stage. Our data show that a common single nucleotide polymorphism (SNP), rs11552708 located in the TNFSF13 gene was significantly associated with IgM levels (p = 5.00×10⁻⁷ in first stage, p = 1.34×10⁻³ in second stage, and p = 4.22×10⁻⁹ when combined). Besides, smoking was identified to be associated with IgM levels in both stages (P<0.05), but there was no significant interaction between smoking and the identified SNP (P>0.05). It is suggested that TNFSF13 may be a susceptibility gene affecting serum IgM levels in Chinese male population.     

Genome-Wide Contribution of Genotype by Environment Interaction to Variation of Diabetes-Related Traits

While genome-wide association studies (GWAS) and candidate gene approaches have identified many genetic variants that contribute to disease risk as main effects, the impact of genotype by environment (GxE) interactions remains rather under-surveyed. To explore the importance of GxE interactions for diabetes-related traits, a tool for Genome-wide Complex Trait Analysis (GCTA) was used to examine GxE variance contribution of 15 macronutrients and lifestyle to the total phenotypic variance of diabetes-related traits at the genome-wide level in a European American population. GCTA identified two key environmental factors making significant contributions to the GxE variance for diabetes-related traits: carbohydrate for fasting insulin (25.1% of total variance, P-nominal = 0.032) and homeostasis model assessment of insulin resistance (HOMA-IR) (24.2% of total variance, P-nominal = 0.035), n-6 polyunsaturated fatty acid (PUFA) for HOMA-β-cell-function (39.0% of total variance, P-nominal = 0.005). To demonstrate and support the results from GCTA, a GxE GWAS was conducted with each of the significant dietary factors and a control E factor (dietary protein), which contributed a non-significant GxE variance. We observed that GxE GWAS for the environmental factor contributing a significant GxE variance yielded more significant SNPs than the control factor. For each trait, we selected all significant SNPs produced from GxE GWAS, and conducted anew the GCTA to estimate the variance they contributed. We noted the variance contributed by these SNPs is higher than that of the control. In conclusion, we utilized a novel method that demonstrates the importance of genome-wide GxE interactions in explaining the variance of diabetes-related traits.     

A Copy Number Variant at the KITLG Locus Likely Confers Risk for Canine Squamous Cell Carcinoma of the Digit

The domestic dog is a robust model for studying the genetics of complex disease susceptibility. The strategies used to develop and propagate modern breeds have resulted in an elevated risk for specific diseases in particular breeds. One example is that of Standard Poodles (STPOs), who have increased risk for squamous cell carcinoma of the digit (SCCD), a locally aggressive cancer that causes lytic bone lesions, sometimes with multiple toe recurrence. However, only STPOs of dark coat color are at high risk; light colored STPOs are almost entirely unaffected, suggesting that interactions between multiple pathways are necessary for oncogenesis. We performed a genome-wide association study (GWAS) on STPOs, comparing 31 SCCD cases to 34 unrelated black STPO controls. The peak SNP on canine chromosome 15 was statistically significant at the genome-wide level (P_raw = 1.60×10⁻⁷; P_genome = 0.0066). Additional mapping resolved the region to the KIT Ligand (KITLG) locus. Comparison of STPO cases to other at-risk breeds narrowed the locus to a 144.9-Kb region. Haplotype mapping among 84 STPO cases identified a minimal region of 28.3 Kb. A copy number variant (CNV) containing predicted enhancer elements was found to be strongly associated with SCCD in STPOs (P = 1.72×10⁻⁸). Light colored STPOs carry the CNV risk alleles at the same frequency as black STPOs, but are not susceptible to SCCD. A GWAS comparing 24 black and 24 light colored STPOs highlighted only the MC1R locus as significantly different between the two datasets, suggesting that a compensatory mutation within the MC1R locus likely protects light colored STPOs from disease. Our findings highlight a role for KITLG in SCCD susceptibility, as well as demonstrate that interactions between the KITLG and MC1R loci are potentially required for SCCD oncogenesis. These findings highlight how studies of breed-limited diseases are useful for disentangling multigene disorders.