Identification of Novel Genetic Risk Loci in Maltese Dogs with Necrotizing Meningoencephalitis and Evidence of a Shared Genetic Risk across Toy Dog Breeds

Necrotizing meningoencephalitis (NME) affects toy and small breed dogs causing progressive, often fatal, inflammation and necrosis in the brain. Genetic risk loci for NME previously were identified in pug dogs, particularly associated with the dog leukocyte antigen (DLA) class II complex on chromosome 12, but have not been investigated in other susceptible breeds. We sought to evaluate Maltese and Chihuahua dogs, in addition to pug dogs, to identify novel or shared genetic risk factors for NME development. Genome-wide association testing of single nucleotide polymorphisms (SNPs) in Maltese dogs with NME identified 2 regions of genome-wide significance on chromosomes 4 (chr4:74522353T>A, p = 8.1×10⁻⁷) and 15 (chr15:53338796A>G, p = 1.5×10⁻⁷). Haplotype analysis and fine-mapping suggests that ILR7 and FBXW7, respectively, both important for regulation of immune system function, could be the underlying associated genes. Further evaluation of these regions and the previously identified DLA II locus across all three breeds, revealed an enrichment of nominal significant SNPs associated with chromosome 15 in pug dogs and DLA II in Maltese and Chihuahua dogs. Meta-analysis confirmed effect sizes the same direction in all three breeds for both the chromosome 15 and DLA II loci (p = 8.6×10^–11 and p = 2.5×10⁻⁷, respectively). This suggests a shared genetic background exists between all breeds and confers susceptibility to NME, but effect sizes might be different among breeds. In conclusion, we identified the first genetic risk factors for NME development in the Maltese, chromosome 4 and chromosome 15, and provide evidence for a shared genetic risk between breeds associated with chromosome 15 and DLA II. Last, DLA II and IL7R both have been implicated in human inflammatory diseases of the central nervous system such as multiple sclerosis, suggesting that similar pharmacotherapeutic targets across species should be investigated.     

Two New Loci for Body-Weight Regulation Identified in a Joint Analysis of Genome-Wide Association Studies for Early-Onset Extreme Obesity in French and German Study Groups

Meta-analyses of population-based genome-wide association studies (GWAS) in adults have recently led to the detection of new genetic loci for obesity. Here we aimed to discover additional obesity loci in extremely obese children and adolescents. We also investigated if these results generalize by estimating the effects of these obesity loci in adults and in population-based samples including both children and adults. We jointly analysed two GWAS of 2,258 individuals and followed-up the best, according to lowest p-values, 44 single nucleotide polymorphisms (SNP) from 21 genomic regions in 3,141 individuals. After this DISCOVERY step, we explored if the findings derived from the extremely obese children and adolescents (10 SNPs from 5 genomic regions) generalized to (i) the population level and (ii) to adults by genotyping another 31,182 individuals (GENERALIZATION step). Apart from previously identified FTO, MC4R, and TMEM18, we detected two new loci for obesity: one in SDCCAG8 (serologically defined colon cancer antigen 8 gene; p = 1.85×10⁻⁸ in the DISCOVERY step) and one between TNKS (tankyrase, TRF1-interacting ankyrin-related ADP-ribose polymerase gene) and MSRA (methionine sulfoxide reductase A gene; p = 4.84×10⁻⁷), the latter finding being limited to children and adolescents as demonstrated in the GENERALIZATION step. The odds ratios for early-onset obesity were estimated at ∼1.10 per risk allele for both loci. Interestingly, the TNKS/MSRA locus has recently been found to be associated with adult waist circumference. In summary, we have completed a meta-analysis of two GWAS which both focus on extremely obese children and adolescents and replicated our findings in a large followed-up data set. We observed that genetic variants in or near FTO, MC4R, TMEM18, SDCCAG8, and TNKS/MSRA were robustly associated with early-onset obesity. We conclude that the currently known major common variants related to obesity overlap to a substantial degree between children and adults.     

Replication and Predictive Value of SNPs Associated with Melanoma and Pigmentation Traits in a Southern European Case-Control Study

Background

Genetic association studies have revealed numerous polymorphisms conferring susceptibility to melanoma. We aimed to replicate previously discovered melanoma-associated single-nucleotide polymorphisms (SNPs) in a Greek case-control population, and examine their predictive value.

Methods

Based on a field synopsis of genetic variants of melanoma (MelGene), we genotyped 284 patients and 284 controls at 34 melanoma-associated SNPs of which 19 derived from GWAS. We tested each one of the 33 SNPs passing quality control for association with melanoma both with and without accounting for the presence of well-established phenotypic risk factors. We compared the risk allele frequencies between the Greek population and the HapMap CEU sample. Finally, we evaluated the predictive ability of the replicated SNPs.

Results

Risk allele frequencies were significantly lower compared to the HapMap CEU for eight SNPs (rs16891982 – SLC45A2, rs12203592 – IRF4, rs258322 – CDK10, rs1805007 – MC1R, rs1805008 - MC1R, rs910873 - PIGU, rs17305573- PIGU, and rs1885120 - MTAP) and higher for one SNP (rs6001027 – PLA2G6) indicating a different profile of genetic susceptibility in the studied population. Previously identified effect estimates modestly correlated with those found in our population (r = 0.72, P<0.0001). The strongest associations were observed for rs401681-T in CLPTM1L (odds ratio [OR] 1.60, 95% CI 1.22–2.10; P = 0.001), rs16891982-C in SCL45A2 (OR 0.51, 95% CI 0.34–0.76; P = 0.001), and rs1805007-T in MC1R (OR 4.38, 95% CI 2.03–9.43; P = 2×10⁻⁵). Nominally statistically significant associations were seen also for another 5 variants (rs258322-T in CDK10, rs1805005-T in MC1R, rs1885120-C in MYH7B, rs2218220-T in MTAP and rs4911442-G in the ASIP region). The addition of all SNPs with nominal significance to a clinical non-genetic model did not substantially improve melanoma risk prediction (AUC for clinical model 83.3% versus 83.9%, p = 0.66).

Conclusion

Overall, our study has validated genetic variants that are likely to contribute to melanoma susceptibility in the Greek population.     

Insights into the Genetic Architecture of Early Stage Age-Related Macular Degeneration: A Genome-Wide Association Study Meta-Analysis

Genetic factors explain a majority of risk variance for age-related macular degeneration (AMD). While genome-wide association studies (GWAS) for late AMD implicate genes in complement, inflammatory and lipid pathways, the genetic architecture of early AMD has been relatively under studied. We conducted a GWAS meta-analysis of early AMD, including 4,089 individuals with prevalent signs of early AMD (soft drusen and/or retinal pigment epithelial changes) and 20,453 individuals without these signs. For various published late AMD risk loci, we also compared effect sizes between early and late AMD using an additional 484 individuals with prevalent late AMD. GWAS meta-analysis confirmed previously reported association of variants at the complement factor H (CFH) (peak P = 1.5×10⁻³¹) and age-related maculopathy susceptibility 2 (ARMS2) (P = 4.3×10⁻²⁴) loci, and suggested Apolipoprotein E (ApoE) polymorphisms (rs2075650; P = 1.1×10⁻⁶) associated with early AMD. Other possible loci that did not reach GWAS significance included variants in the zinc finger protein gene GLI3 (rs2049622; P = 8.9×10⁻⁶) and upstream of GLI2 (rs6721654; P = 6.5×10⁻⁶), encoding retinal Sonic hedgehog signalling regulators, and in the tyrosinase (TYR) gene (rs621313; P = 3.5×10⁻⁶), involved in melanin biosynthesis. For a range of published, late AMD risk loci, estimated effect sizes were significantly lower for early than late AMD. This study confirms the involvement of multiple established AMD risk variants in early AMD, but suggests weaker genetic effects on the risk of early AMD relative to late AMD. Several biological processes were suggested to be potentially specific for early AMD, including pathways regulating RPE cell melanin content and signalling pathways potentially involved in retinal regeneration, generating hypotheses for further investigation.     

A Genome-Wide Association Study Identifies Potential Susceptibility Loci for Hirschsprung Disease

Hirschsprung disease (HSCR) is a congenital and heterogeneous disorder characterized by the absence of intramural nervous plexuses along variable lengths of the hindgut. Although RET is a well-established risk factor, a recent genome-wide association study (GWAS) of HSCR has identified NRG1 as an additional susceptibility locus. To discover additional risk loci, we performed a GWAS of 123 sporadic HSCR patients and 432 unaffected controls using a large-scale platform with coverage of over 1 million polymorphic markers. The result was that our study replicated the findings of RET-CSGALNACT2-RASGEF1A genomic region (_rawP = 5.69×10⁻¹⁹ before a Bonferroni correction; _corrP = 4.31×10⁻¹³ after a Bonferroni correction) and NRG1 as susceptibility loci. In addition, this study identified SLC6A20 (_adjP = 2.71×10⁻⁶), RORA (_adjP = 1.26×10⁻⁵), and ABCC9 (_adjP = 1.86×10⁻⁵) as new potential susceptibility loci under adjusting the already known loci on the RET-CSGALNACT2-RASGEF1A and NRG1 regions, although none of the SNPs in these genes passed the Bonferroni correction. In further subgroup analysis, the RET-CSGALNACT2-RASGEF1A genomic region was observed to have different significance levels among subgroups: short-segment (S-HSCR, _corrP = 1.71×10⁻⁵), long-segment (L-HSCR, _corrP = 6.66×10⁻⁴), and total colonic aganglionosis (TCA, _corrP>0.05). This differential pattern in the significance level suggests that other genomic loci or mechanisms may affect the length of aganglionosis in HSCR subgroups during enteric nervous system (ENS) development. Although functional evaluations are needed, our findings might facilitate improved understanding of the mechanisms of HSCR pathogenesis.     

Skin Transcriptome Profiles Associated with Skin Color in Chickens

Nutritional and medicinal benefits have been attributed to the consumption of tissues from the black-boned chickens in oriental countries. Lueyang black-boned chicken is one of the native chicken breeds. However, some birds may instead have white or lighter skin, which directly causes economic losses every year. Previous studies of pigmentation have focused on a number of genes that may play important roles in coat color regulation. Illumina2000 sequencing technology was used to catalog the global gene expression profiles in the skin of the Lueyang chicken with white versus black skin. A total of 18,608 unigenes were assembled from the reads obtained from the skin of the white and black chickens. A total of 649 known genes were differentially expressed in the black versus white chickens, with 314 genes that were up regulated and 335 genes that were down-regulated, and a total of 162 novel genes were differentially expressed in the black versus white chickens, consisting of 73 genes that were up-regulated (including 4 highly expressed genes that were expressed exclusively in the skin of the black chickens) and 89 genes that were down-regulated. There were also a total of 8 known coat-color genes expressed in previous studies (ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R). In this study, 4 of which showed greater expression in the black chickens, and several were up-regulated, such as KIT, ASIP, TYR and OCA2. To our surprise, KITLG, MITF and MC1R showed no significant difference in expression between the black- and white-skinned chickens, and the expression of TYRP1 was not detected in either skin color. The expression of ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R was validated by real-time quantitative polymerase chain reaction (qPCR), and the results of the qPCR were consistent with the RNA-seq. This study provides several candidate genes that may be associated with the development of black versus white skin. More importantly, the fact that the MC1R gene showed no significant difference in expression between the black and white chickens is of particular interest for future studies that aim to elucidate its functional role in the regulation of skin color.     

Genome-Wide Association Study in BRCA1 Mutation Carriers Identifies Novel Loci Associated with Breast and Ovarian Cancer Risk

BRCA1-associated breast and ovarian cancer risks can be modified by common genetic variants. To identify further cancer risk-modifying loci, we performed a multi-stage GWAS of 11,705 BRCA1 carriers (of whom 5,920 were diagnosed with breast and 1,839 were diagnosed with ovarian cancer), with a further replication in an additional sample of 2,646 BRCA1 carriers. We identified a novel breast cancer risk modifier locus at 1q32 for BRCA1 carriers (rs2290854, P = 2.7×10⁻⁸, HR = 1.14, 95% CI: 1.09–1.20). In addition, we identified two novel ovarian cancer risk modifier loci: 17q21.31 (rs17631303, P = 1.4×10⁻⁸, HR = 1.27, 95% CI: 1.17–1.38) and 4q32.3 (rs4691139, P = 3.4×10⁻⁸, HR = 1.20, 95% CI: 1.17–1.38). The 4q32.3 locus was not associated with ovarian cancer risk in the general population or BRCA2 carriers, suggesting a BRCA1-specific association. The 17q21.31 locus was also associated with ovarian cancer risk in 8,211 BRCA2 carriers (P = 2×10⁻⁴). These loci may lead to an improved understanding of the etiology of breast and ovarian tumors in BRCA1 carriers. Based on the joint distribution of the known BRCA1 breast cancer risk-modifying loci, we estimated that the breast cancer lifetime risks for the 5% of BRCA1 carriers at lowest risk are 28%–50% compared to 81%–100% for the 5% at highest risk. Similarly, based on the known ovarian cancer risk-modifying loci, the 5% of BRCA1 carriers at lowest risk have an estimated lifetime risk of developing ovarian cancer of 28% or lower, whereas the 5% at highest risk will have a risk of 63% or higher. Such differences in risk may have important implications for risk prediction and clinical management for BRCA1 carriers.     

Genome-Wide Association Studies Identify Two Novel BMP15 Mutations Responsible for an Atypical Hyperprolificacy Phenotype in Sheep

Some sheep breeds are naturally prolific, and they are very informative for the studies of reproductive genetics and physiology. Major genes increasing litter size (LS) and ovulation rate (OR) were suspected in the French Grivette and the Polish Olkuska sheep populations, respectively. To identify genetic variants responsible for the highly prolific phenotype in these two breeds, genome-wide association studies (GWAS) followed by complementary genetic and functional analyses were performed. Highly prolific ewes (cases) and normal prolific ewes (controls) from each breed were genotyped using the Illumina OvineSNP50 Genotyping Beadchip. In both populations, an X chromosome region, close to the BMP15 gene, harbored clusters of markers with suggestive evidence of association at significance levels between 1E⁻⁰⁵ and 1E⁻⁰⁷. The BMP15 candidate gene was then sequenced, and two novel non-conservative mutations called FecX^Gr and FecX^O were identified in the Grivette and Olkuska breeds, respectively. The two mutations were associated with the highly prolific phenotype (p_FecX^Gr = 5.98E⁻⁰⁶ and p_FecX^O = 2.55E⁻⁰⁸). Homozygous ewes for the mutated allele showed a significantly increased prolificacy (FecX^Gr/FecX^Gr, LS = 2.50±0.65 versus FecX⁺/FecX^Gr, LS = 1.93±0.42, p<1E⁻⁰³ and FecX^O/FecX^O, OR = 3.28±0.85 versus FecX⁺/FecX^O, OR = 2.02±0.47, p<1E⁻⁰³). Both mutations are located in very well conserved motifs of the protein and altered the BMP15 signaling activity in vitro using a BMP-responsive luciferase test in COV434 granulosa cells. Thus, we have identified two novel mutations in the BMP15 gene associated with increased LS and OR. Notably, homozygous FecX^Gr/FecX^Gr Grivette and homozygous FecX^O/FecX^O Olkuska ewes are hyperprolific in striking contrast with the sterility exhibited by all other known homozygous BMP15 mutations. Our results bring new insights into the key role played by the BMP15 protein in ovarian function and could contribute to a better understanding of the pathogenesis of women′s fertility disorders.     

Insight in Genome-Wide Association of Metabolite Quantitative Traits by Exome Sequence Analyses

Metabolite quantitative traits carry great promise for epidemiological studies, and their genetic background has been addressed using Genome-Wide Association Studies (GWAS). Thus far, the role of less common variants has not been exhaustively studied. Here, we set out a GWAS for metabolite quantitative traits in serum, followed by exome sequence analysis to zoom in on putative causal variants in the associated genes. ¹H Nuclear Magnetic Resonance (¹H-NMR) spectroscopy experiments yielded successful quantification of 42 unique metabolites in 2,482 individuals from The Erasmus Rucphen Family (ERF) study. Heritability of metabolites were estimated by SOLAR. GWAS was performed by linear mixed models, using HapMap imputations. Based on physical vicinity and pathway analyses, candidate genes were screened for coding region variation using exome sequence data. Heritability estimates for metabolites ranged between 10% and 52%. GWAS replicated three known loci in the metabolome wide significance: CPS1 with glycine (P-value  = 1.27×10⁻³²), PRODH with proline (P-value  = 1.11×10⁻¹⁹), SLC16A9 with carnitine level (P-value  = 4.81×10⁻¹⁴) and uncovered a novel association between DMGDH and dimethyl-glycine (P-value  = 1.65×10⁻¹⁹) level. In addition, we found three novel, suggestively significant loci: TNP1 with pyruvate (P-value  = 1.26×10⁻⁸), KCNJ16 with 3-hydroxybutyrate (P-value  = 1.65×10⁻⁸) and 2p12 locus with valine (P-value  = 3.49×10⁻⁸). Exome sequence analysis identified potentially causal coding and regulatory variants located in the genes CPS1, KCNJ2 and PRODH, and revealed allelic heterogeneity for CPS1 and PRODH. Combined GWAS and exome analyses of metabolites detected by high-resolution ¹H-NMR is a robust approach to uncover metabolite quantitative trait loci (mQTL), and the likely causative variants in these loci. It is anticipated that insight in the genetics of intermediate phenotypes will provide additional insight into the genetics of complex traits.     

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.     

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.     

Obesity-Related Genomic Loci Are Associated with Type 2 Diabetes in a Han Chinese Population

Background and Aims

Obesity is a well-known risk factor for type 2 diabetes. Genome-wide association studies have identified a number of genetic loci associated with obesity. The aim of this study is to examine the contribution of obesity-related genomic loci to type 2 diabetes in a Chinese population.

Methods

We successfully genotyped 18 obesity-related single nucleotide polymorphisms among 5338 type 2 diabetic patients and 4663 controls. Both individual and joint effects of these single nucleotide polymorphisms on type 2 diabetes and quantitative glycemic traits (assessing β-cell function and insulin resistance) were analyzed using logistic and linear regression models, respectively.

Results

Two single nucleotide polymorphisms near MC4R and GNPDA2 genes were significantly associated with type 2 diabetes before adjusting for body mass index and waist circumference (OR (95% CI) = 1.14 (1.06, 1.22) for the A allele of rs12970134, P = 4.75×10⁻⁴; OR (95% CI) = 1.10 (1.03, 1.17) for the G allele of rs10938397, P = 4.54×10⁻³). When body mass index and waist circumference were further adjusted, the association of MC4R with type 2 diabetes remained significant (P = 1.81×10⁻²) and that of GNPDA2 was attenuated (P = 1.26×10⁻¹), suggesting the effect of the locus including GNPDA2 on type 2 diabetes may be mediated through obesity. Single nucleotide polymorphism rs2260000 within BAT2 was significantly associated with type 2 diabetes after adjusting for body mass index and waist circumference (P = 1.04×10⁻²). In addition, four single nucleotide polymorphisms (near or within SEC16B, BDNF, MAF and PRL genes) showed significant associations with quantitative glycemic traits in controls even after adjusting for body mass index and waist circumference (all P values<0.05).

Conclusions

This study indicates that obesity-related genomic loci were associated with type 2 diabetes and glycemic traits in the Han Chinese population.     

Modeling of Environmental Effects in Genome-Wide Association Studies Identifies SLC2A2 and HP as Novel Loci Influencing Serum Cholesterol Levels

Genome-wide association studies (GWAS) have identified 38 larger genetic regions affecting classical blood lipid levels without adjusting for important environmental influences. We modeled diet and physical activity in a GWAS in order to identify novel loci affecting total cholesterol, LDL cholesterol, HDL cholesterol, and triglyceride levels. The Swedish (SE) EUROSPAN cohort (N _SE = 656) was screened for candidate genes and the non-Swedish (NS) EUROSPAN cohorts (N _NS = 3,282) were used for replication. In total, 3 SNPs were associated in the Swedish sample and were replicated in the non-Swedish cohorts. While SNP rs1532624 was a replication of the previously published association between CETP and HDL cholesterol, the other two were novel findings. For the latter SNPs, the p-value for association was substantially improved by inclusion of environmental covariates: SNP rs5400 (p _{SE,unadjusted} = 3.6×10⁻⁵, p _SE,adjusted = 2.2×10⁻⁶, p _{NS,unadjusted} = 0.047) in the SLC2A2 (Glucose transporter type 2) and rs2000999 (p _{SE,unadjusted} = 1.1×10⁻³, p _SE,adjusted = 3.8×10⁻⁴, p _{NS,unadjusted} = 0.035) in the HP gene (Haptoglobin-related protein precursor). Both showed evidence of association with total cholesterol. These results demonstrate that inclusion of important environmental factors in the analysis model can reveal new genetic susceptibility loci.     

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.     

Genome-Wide Association Study to Identify the Genetic Determinants of Otitis Media Susceptibility in Childhood

Background

Otitis media (OM) is a common childhood disease characterised by middle ear inflammation and effusion. Susceptibility to recurrent acute OM (rAOM; ≥3 episodes of AOM in 6 months) and chronic OM with effusion (COME; MEE ≥3 months) is 40–70% heritable. Few underlying genes have been identified to date, and no genome-wide association study (GWAS) of OM has been reported.

Methods and Findings

Data for 2,524,817 single nucleotide polymorphisms (SNPs; 535,544 quality-controlled SNPs genotyped by Illumina 660W-Quad; 1,989,273 by imputation) were analysed for association with OM in 416 cases and 1,075 controls from the Western Australian Pregnancy Cohort (Raine) Study. Logistic regression analyses under an additive model undertaken in GenABEL/ProbABEL adjusting for population substructure using principal components identified SNPs at CAPN14 (rs6755194: OR = 1.90; 95%CI 1.47–2.45; P_adj-PCA = 8.3×10⁻⁷) on chromosome 2p23.1 as the top hit, with independent effects (rs1862981: OR = 1.60; 95%CI 1.29–1.99; P_adj-PCA = 2.2×10⁻⁵) observed at the adjacent GALNT14 gene. In a gene-based analysis in VEGAS, BPIFA3 (P_Gene = 2×10⁻⁵) and BPIFA1 (P_Gene = 1.07×10⁻⁴) in the BPIFA gene cluster on chromosome 20q11.21 were the top hits. In all, 32 genomic regions show evidence of association (P_adj-PCA<10⁻⁵) in this GWAS, with pathway analysis showing a connection between top candidates and the TGFβ pathway. However, top and tag-SNP analysis for seven selected candidate genes in this pathway did not replicate in 645 families (793 affected individuals) from the Western Australian Family Study of Otitis Media (WAFSOM). Lack of replication may be explained by sample size, difference in OM disease severity between primary and replication cohorts or due to type I error in the primary GWAS.

Conclusions

This first discovery GWAS for an OM phenotype has identified CAPN14 and GALNT14 on chromosome 2p23.1 and the BPIFA gene cluster on chromosome 20q11.21 as novel candidate genes which warrant further analysis in cohorts matched more precisely for clinical phenotypes.     

First Autonomous Bio-Optical Profiling Float in the Gulf of Mexico Reveals Dynamic Biogeochemistry in Deep Waters

Profiling floats equipped with bio-optical sensors well complement ship-based and satellite ocean color measurements by providing highly-resolved time-series data on the vertical structure of biogeochemical processes in oceanic waters. This is the first study to employ an autonomous profiling (APEX) float in the Gulf of Mexico for measuring spatiotemporal variability in bio-optics and hydrography. During the 17-month deployment (July 2011 to December 2012), the float mission collected profiles of temperature, salinity, chlorophyll fluorescence, particulate backscattering (b_bp), and colored dissolved organic matter (CDOM) fluorescence from the ocean surface to a depth of 1,500 m. Biogeochemical variability was characterized by distinct depth trends and local “hot spots”, including impacts from mesoscale processes associated with each of the water masses sampled, from ambient deep waters over the Florida Plain, into the Loop Current, up the Florida Canyon, and eventually into the Florida Straits. A deep chlorophyll maximum (DCM) occurred between 30 and 120 m, with the DCM depth significantly related to the unique density layer ρ = 1023.6 (R² = 0.62). Particulate backscattering, b_bp, demonstrated multiple peaks throughout the water column, including from phytoplankton, deep scattering layers, and resuspension. The bio-optical relationship developed between b_bp and chlorophyll (R² = 0.49) was compared to a global relationship and could significantly improve regional ocean-color algorithms. Photooxidation and autochthonous production contributed to CDOM distributions in the upper water column, whereas in deep water, CDOM behaved as a semi-conservative tracer of water masses, demonstrating a tight relationship with density (R² = 0.87). In the wake of the Deepwater Horizon oil spill, this research lends support to the use of autonomous drifting profilers as a powerful tool for consideration in the design of an expanded and integrated observing network for the Gulf of Mexico.     

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.     

Replication of Association between Schizophrenia and Chromosome 6p21-6p22.1 Polymorphisms in Chinese Han Population

Chromosome 6p21-p22.1, spanning the extended major histocompatibility complex (MHC) region, is a highly polymorphic, gene-dense region. It has been identified as a susceptibility locus of schizophrenia in Europeans, Japanese, and Chinese. In our previous two-stage genome-wide association study (GWAS), polymorphisms of zinc finger with KRAB and SCAN domains 4 (ZKSCAN4), nuclear factor-κB-activating protein-like (NKAPL), and piggyBac transposable element derived 1 (PGBD1), localized to chromosome 6p21-p22.1, were strongly associated with schizophrenia. To further investigate the association between polymorphisms at this locus and schizophrenia in the Chinese Han population, we selected eight other single-nucleotide polymorphisms (SNPs) distributed in or near these genes for a case-control association study in an independent sample of 902 cases and 1,091 healthy controls in an attempt to replicate the GWAS results. Four of these eight SNPs (rs12214383, rs1150724, rs3800324, and rs1997660) displayed a nominal difference in allele frequencies between the case and control groups. The association between two of these SNPs and schizophrenia were significant even after Bonferroni correction (rs12000: allele A>G, P = 2.50E-04, odds ratio [OR] = 1.27, 95% confidence interval [CI] = 1.12–1.45; rs1150722: allele C>T, P = 4.28E-05, OR = 0.55, 95% CI = 0.41–0.73). Haplotype ATTGACGC, comprising these eight SNPs (rs2235359, rs2185955, rs12214383, rs12000, rs1150724, rs1150722, rs3800324, and rs1997660), was significantly associated with schizophrenia (P = 6.60E-05). We also performed a combined study of this replication sample and the first-stage GWAS sample. The combined study revealed that rs12000 and rs1150722 were still strongly associated with schizophrenia (rs12000: allele G>A, P _combined  = 0.0019, OR = 0.81; rs1150722: allele G>A, P _combined  = 3.00E-04, OR = 0.61). These results support our findings that locus 6p21-p22.1 is significantly associated with schizophrenia in the Chinese Han population and encourage further studies of the functions of these genetic factors.     

Thrombotic Antiphospholipid Syndrome Shows Strong Haplotypic Association with SH2B3-ATXN2 Locus

Background

Thrombotic antiphospholipid syndrome is defined as a complex form of thrombophilia that is developed by a fraction of antiphospholipid antibody (aPLA) carriers. Little is known about the genetic risk factors involved in thrombosis development among aPLA carriers.

Methods

To identify new loci conferring susceptibility to thrombotic antiphospholipid syndrome, a two-stage genotyping strategy was performed. In stage one, 19,000 CNV loci were genotyped in 14 thrombotic aPLA+ patients and 14 healthy controls by array-CGH. In stage two, significant CNV loci were fine-mapped in a larger cohort (85 thrombotic aPLA+, 100 non-thrombotic aPLA+ and 569 healthy controls).

Results

Array-CGH and fine-mapping analysis led to the identification of 12q24.12 locus as a new susceptibility locus for thrombotic APS. Within this region, a TAC risk haplotype comprising one SNP in SH2B3 gene (rs3184504) and two SNPs in ATXN2 gene (rs10774625 and rs653178) exhibited the strongest association with thrombotic antiphospholipid syndrome (p-value = 5,9 × 10⁻⁴ OR 95% CI 1.84 (1.32–2.55)).

Conclusion

The presence of a TAC risk haplotype in ATXN2-SH2B3 locus may contribute to increased thrombotic risk in aPLA carriers.     

A Genome-Wide Association Study of Circulating Galectin-3

Galectin-3 is a lectin involved in fibrosis, inflammation and proliferation. Increased circulating levels of galectin-3 have been associated with various diseases, including cancer, immunological disorders, and cardiovascular disease. To enhance our knowledge on galectin-3 biology we performed the first genome-wide association study (GWAS) using the Illumina HumanCytoSNP-12 array imputed with the HapMap 2 CEU panel on plasma galectin-3 levels in 3,776 subjects and follow-up genotyping in an additional 3,516 subjects. We identified 2 genome wide significant loci associated with plasma galectin-3 levels. One locus harbours the LGALS3 gene (rs2274273; P = 2.35×10⁻¹⁸⁸) and the other locus the ABO gene (rs644234; P = 3.65×10⁻⁴⁷). The variance explained by the LGALS3 locus was 25.6% and by the ABO locus 3.8% and jointly they explained 29.2%. Rs2274273 lies in high linkage disequilibrium with two non-synonymous SNPs (rs4644; r² = 1.0, and rs4652; r² = 0.91) and wet lab follow-up genotyping revealed that both are strongly associated with galectin-3 levels (rs4644; P = 4.97×10⁻⁴⁶⁵ and rs4652 P = 1.50×10⁻⁴²¹) and were also associated with LGALS3 gene-expression. The origins of our associations should be further validated by means of functional experiments.