Arylsulphatase A Pseudodeficiency (ARSA-PD), hypertension and chronic renal disease in Aboriginal Australians

Chronic renal disease (CRD) associated with cardiovascular disease (CVD) and/or type 2 diabetes (T2D) is a significant problem in Aboriginal Australians. Whole exome sequencing data (N = 72) showed enrichment for ClinVar pathogenic variants in gene sets/pathways linking lipoprotein, lipid and glucose metabolism. The top Ingenuity Pathway Analysis canonical pathways were Farsenoid X Receptor and Retinoid Receptor (FXR/RXR; (P = 1.86 × 10⁻⁷), Liver X Receptor and Retinoid Receptor (LXR/RXR; P = 2.88 × 10⁻⁶), and atherosclerosis signalling (P = 3.80 × 10⁻⁶). Top pathways/processes identified using Enrichr included: Reactome 2016 chylomicron-mediated lipid transport (P = 3.55 × 10⁻⁷); Wiki 2016 statin (P = 8.29 × 10⁻⁸); GO Biological Processes 2017 chylomicron remodelling (P = 1.92 × 10⁻⁸). ClinVar arylsulfatase A pseudodeficiency (ARSA-PD) pathogenic variants were common, including the missense variant c.511 G > A (p.Asp171Asn; rs74315466; frequency 0.44) only reported in Polynesians. This variant is in cis with known ARSA-PD 3′ regulatory c.*96 A > G (rs6151429; frequency 0.47) and missense c.1055 A > G (p.Asn352Ser; rs2071421; frequency 0.47) variants. These latter two variants are associated with T2D (risk haplotype GG; odds ratio 2.67; 95% CI 2.32–3.08; P = 2.43 × 10⁻⁴) in genome-wide association data (N = 402), but are more strongly associated with quantitative traits (DBP, SBP, ACR, eGFR) for hypertension and renal function in non-diabetic than diabetic subgroups. Traits associated with CVD, CRD and T2D in Aboriginal Australians provide novel insight into function of ARSA-PD variants.

First genome-wide association study in an Australian aboriginal population provides insights into genetic risk factors for body mass index and type 2 diabetes

A body mass index (BMI) >22kg/m² is a risk factor for type 2 diabetes (T2D) in Aboriginal Australians. To identify loci associated with BMI and T2D we undertook a genome-wide association study using 1,075,436 quality-controlled single nucleotide polymorphisms (SNPs) genotyped (Illumina 2.5M Duo Beadchip) in 402 individuals in extended pedigrees from a Western Australian Aboriginal community. Imputation using the thousand genomes (1000G) reference panel extended the analysis to 6,724,284 post quality-control autosomal SNPs. No associations achieved genome-wide significance, commonly accepted as P<5x10^-8. Nevertheless, genes/pathways in common with other ethnicities were identified despite the arrival of Aboriginal people in Australia >45,000 years ago. The top hit (rs10868204 P_genotyped = 1.50x10^-6; rs11140653 P_{imputed_1000G} = 2.90x10^-7) for BMI lies 5’ of NTRK2, the type 2 neurotrophic tyrosine kinase receptor for brain-derived neurotrophic factor (BDNF) that regulates energy balance downstream of melanocortin-4 receptor (MC4R). PIK3C2G (rs12816270 P_genotyped = 8.06x10^-6; rs10841048 P_{imputed_1000G} = 6.28x10^-7) was associated with BMI, but not with T2D as reported elsewhere. BMI also associated with CNTNAP2 (rs6960319 P_genotyped = 4.65x10^-5; rs13225016 P_{imputed_1000G} = 6.57x10^-5), previously identified as the strongest gene-by-environment interaction for BMI in African-Americans. The top hit (rs11240074 P_genotyped = 5.59x10^-6, P_{imputed_1000G} = 5.73x10^-6) for T2D lies 5’ of BCL9 that, along with TCF7L2, promotes beta-catenin’s transcriptional activity in the WNT signaling pathway. Additional hits occurred in genes affecting pancreatic (KCNJ6, KCNA1) and/or GABA (GABRR1, KCNA1) functions. Notable associations observed for genes previously identified at genome-wide significance in other populations included MC4R (P_genotyped = 4.49x10^-4) for BMI and IGF2BP2 P_{imputed_1000G} = 2.55x10^-6) for T2D. Our results may provide novel functional leads in understanding disease pathogenesis in this Australian Aboriginal population.

Genetic and functional evidence for a locus controlling otitis media at chromosome 10q26.3

BACKGROUND

Otitis media (OM) is a common childhood disease characterised by middle ear effusion and inflammation. Susceptibility to recurrent acute OM and chronic OM with effusion is 40-70% heritable. Linkage studies provide evidence for multiple putative OM susceptibility loci. This study attempts to replicate these linkages in a Western Australian (WA) population, and to identify the etiological gene(s) in a replicated region.

METHODS

Microsatellites were genotyped in 468 individuals from 101 multicase families (208 OM cases) from the WA Family Study of OM (WAFSOM) and non-parametric linkage analysis carried out in ALLEGRO. Association mapping utilized dense single nucleotide polymorphism (SNP) data extracted from Illumina 660 W-Quad analysis of 256 OM cases and 575 controls from the WA Pregnancy Cohort (Raine) Study. Logistic regression analysis was undertaken in ProbABEL. RT-PCR was used to compare gene expression in paired adenoid and tonsil samples, and in epithelial and macrophage cell lines. Comparative genomics methods were used to identify putative regulatory elements and transcription factor binding sites potentially affected by associated SNPs.

RESULTS

Evidence for linkage was observed at 10q26.3 (Zlr = 2.69; P = 0.0036; D10S1770) with borderline evidence for linkage at 10q22.3 (Zlr = 1.64; P = 0.05; D10S206). No evidence for linkage was seen at 3p25.3, 17q12, or 19q13.43. Peak association at 10q26.3 was in the intergenic region between TCERG1L and PPP2R2D (rs7922424; P = 9.47 × 10-6), immediately under the peak of linkage. Independent associations were observed at DOCK1 (rs9418832; P = 7.48 × 10-5) and ADAM12 (rs7902734; P = 8.04 × 10-4). RT-PCR analysis confirmed expression of all 4 genes in adenoid samples. ADAM12, DOCK1 and PPP2R2D, but not TCERG1L, were expressed in respiratory epithelial and macrophage cell lines. A significantly associated polymorphism (rs7087384) in strong LD with the top SNP (rs7922424; r2 = 0.97) alters a transcription factor binding site (CREB/CREBP) in the intergenic region between TCERG1L and PPP2R2D.

CONCLUSIONS

OM linkage was replicated at 10q26.3. Whilst multiple genes could contribute to this linkage, the weight of evidence supports PPP2R2D, a TGF-β/Activin/Nodal pathway modulator, as the more likely functional candidate lying immediately under the linkage peak for OM susceptibility at chromosome 10q26.3.

Alternative recruitment strategies influence saliva sample return rates in community-based genetic association studies

Collection of saliva for DNA extraction has created new opportunities to recruit participants from the community for genetic association studies. However, sample return rates are variable. No prior study has specifically addressed how study design impacts sample return. Using data from three large-scale genetic association studies we compared recruitment strategy and sample return rates. We found highly significant differences in sample return rates between the studies. In studies that recruited retrospectively, overall returns were much lower from families with a self-limiting condition who provided samples at a research centre or home visit, than adult elderly individuals with a chronic disease who provided samples by post (59% vs. 84%). Prospective recruitment was associated with high agreement to participate (72%), but subsequent low return of actual saliva samples (42%). A telephone call had marginal effect on recruitment in a retrospective family study, but significantly improved returns in a prospective family study. We found no effect upon DNA yield comparing observed versus unobserved sample collection, or between male and female adult participants. Overall, study design significantly impacts upon response rates for genetic association studies recruiting from the community. Our findings will help researchers in constructing and costing a recruitment protocol.

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.

FBXO11, a regulator of the TGFβ pathway, is associated with severe otitis media in Western Australian children

Otitis media (OM) is a common childhood disease characterised by middle ear inflammation following infection. Susceptibility to recurrent acute OM (rAOM) and chronic OM with effusion (COME) is highly heritable. Two murine mutants, Junbo and Jeff, spontaneously develop severe OM with similar phenotypes to human disease. Fine-mapping of these mutants identified two genes (Evi1 and Fbxo11) that interact with the transforming growth factor β (TGFβ) signalling pathway. We investigated these genes, as well as four Sma- and Mad-related (SMAD) genes of the TGFβ pathway, as candidate rAOM/COME susceptibility genes in two predominantly Caucasian populations. Single-nucleotide polymorphisms (SNPs) within FBXO11 (family-based association testing Z-Score=2.61; P(best)=0.009) were associated with severe OM in family-based analysis of 434 families (561 affected individuals) from the Western Australian Family Study of OM. The FBXO11 association was replicated by directed analysis of Illumina 660W-Quad Beadchip data available for 253 cases and 866 controls (OR=1.55 (95% CI 1.28-1.89); P(best)=6.9 × 10(-6)) available within the Western Australian Pregnancy Cohort (Raine) Study. Combined primary and replication results show P(combined)=2.98 × 10(-6). Neither cohort showed an association with EVI1 variants. Family-based associations at SMAD2 (P=0.038) and SMAD4 (P=0.048) were not replicated. Together, these data provide strong evidence for FBXO11 as a susceptibility gene for severe OM.