Mapping eQTLs in the Norfolk Island genetic isolate identifies candidate genes for CVD risk traits

Enhancing selection of alcohol consumption-associated genes by random forest

Machine learning methods have been used in identifying omics markers for a variety of phenotypes. We aimed to examine whether a supervised machine learning algorithm can improve identification of alcohol-associated transcriptomic markers. In this study, we analysed array-based, whole-blood derived expression data for 17 873 gene transcripts in 5508 Framingham Heart Study participants. By using the Boruta algorithm, a supervised random forest (RF)-based feature selection method, we selected twenty-five alcohol-associated transcripts. In a testing set (30 % of entire study participants), AUC (area under the receiver operating characteristics curve) of these twenty-five transcripts were 0·73, 0·69 and 0·66 for non-drinkers v. moderate drinkers, non-drinkers v. heavy drinkers and moderate drinkers v. heavy drinkers, respectively. The AUC of the selected transcripts by the Boruta method were comparable to those identified using conventional linear regression models, for example, AUC of 1958 transcripts identified by conventional linear regression models (false discovery rate < 0·2) were 0·74, 0·66 and 0·65, respectively. With Bonferroni correction for the twenty-five Boruta method-selected transcripts and three CVD risk factors (i.e. at P < 6·7e-4), we observed thirteen transcripts were associated with obesity, three transcripts with type 2 diabetes and one transcript with hypertension. For example, we observed that alcohol consumption was inversely associated with the expression of DOCK4, IL4R, and SORT1, and DOCK4 and SORT1 were positively associated with obesity, and IL4R was inversely associated with hypertension. In conclusion, using a supervised machine learning method, the RF-based Boruta algorithm, we identified novel alcohol-associated gene transcripts.

Exploring the Functional Basis of Epigenetic Aging in Relation to Body Fat Phenotypes in the Norfolk Island Cohort

DNA methylation is an epigenetic factor that is modifiable and can change over a lifespan. While many studies have identified methylation sites (CpGs) related to aging, the relationship of these to gene function and age-related disease phenotypes remains unclear. This research explores this question by testing for the conjoint association of age-related CpGs with gene expression and the relation of these to body fat phenotypes. The study included blood-based gene transcripts and intragenic CpG methylation data from Illumina 450 K arrays in 74 healthy adults from the Norfolk Island population. First, a series of regression analyses were performed to detect associations between gene transcript level and intragenic CpGs and their conjoint relationship with age. Second, we explored how these age-related expression CpGs (eCpGs) correlated with obesity-related phenotypes, including body fat percentage, body mass index, and waist-to-hip ratio. We identified 35 age-related eCpGs associated with age. Of these, ten eCpGs were associated with at least one body fat phenotype. Collagen Type XI Alpha 2 Chain (COL11A2), Complement C1s (C1s), and four and a half LIM domains 2 (FHL2) genes were among the most significant genes with multiple eCpGs associated with both age and multiple body fat phenotypes. The COL11A2 gene contributes to the correct assembly of the extracellular matrix in maintaining the healthy structural arrangement of various components, with the C1s gene part of complement systems functioning in inflammation. Moreover, FHL2 expression was upregulated under hypermethylation in both blood and adipose tissue with aging. These results suggest new targets for future studies and require further validation to confirm the specific function of these genes on body fat regulation.

Blood transcriptomic biomarkers of alcohol consumption and cardiovascular disease risk factors: the Framingham Heart Study

BACKGROUND

The relations of alcohol consumption and gene expression remain to be elucidated.

MATERIALS AND METHODS

We examined cross-sectional associations between alcohol consumption and whole blood derived gene expression levels and between alcohol-associated genes and obesity, hypertension, and diabetes in 5531 Framingham Heart Study (FHS) participants.

RESULTS

We identified 25 alcohol-associated genes. We further showed cross-sectional associations of 16 alcohol-associated genes with obesity, nine genes with hypertension, and eight genes with diabetes at P < 0.002. For example, we observed decreased expression of PROK2 (β = -0.0018; 95%CI: -0.0021, -0.0007; P = 6.5e - 5) and PAX5 (β = -0.0014; 95%CI: -0.0021, -0.0007; P = 6.5e - 5) per 1 g/day increase in alcohol consumption. Consistent with our previous observation on the inverse association of alcohol consumption with obesity and positive association of alcohol consumption with hypertension, we found that PROK2 was positively associated with obesity (OR = 1.42; 95%CI: 1.17, 1.72; P = 4.5e - 4) and PAX5 was negatively associated with hypertension (OR = 0.73; 95%CI: 0.59, 0.89; P = 1.6e - 3). We also observed that alcohol consumption was positively associated with expression of ABCA13 (β = 0.0012; 95%CI: 0.0007, 0.0017; P = 1.3e - 6) and ABCA13 was positively associated with diabetes (OR = 2.57; 95%CI: 1.73, 3.84; P = 3.5e - 06); this finding, however, was inconsistent with our observation of an inverse association between alcohol consumption and diabetes.

CONCLUSIONS

We showed strong cross-sectional associations between alcohol consumption and expression levels of 25 genes in FHS participants. Nonetheless, complex relationships exist between alcohol-associated genes and CVD risk factors.

Evaluating the suitability of current mitochondrial DNA interpretation guidelines for multigenerational whole mitochondrial genome comparisons

Sanger sequencing of the mitochondrial DNA (mtDNA) control region was previously the only method available for forensic casework involving degraded samples from skeletal remains. The introduction of Next Generation Sequencing (NGS) has transformed genetic data generation and human identification using mtDNA. Whole mitochondrial genome (mtGenome) analysis is now being introduced into forensic laboratories around the world to analyze historical remains. Research into large pedigrees using the mtGenome is critical to evaluate currently available interpretation guidelines for mtDNA analysis, which were developed for comparisons using the control region. This study included mtGenomes from 225 individuals from the last four generations of the Norfolk Island (NI) genetic isolate pedigree consisting of 49 distinct maternal lineages. The data from these individuals were arranged into 2339 maternally related pairs separated by up to 18 meioses. Our results show that 97.3% of maternally related pairs were concordant at all nucleotide positions, resulting in the correct interpretation of “Cannot Exclude”; 2.7% of pairs produced an “Inconclusive” result, and there were no instances of false exclusion. While these results indicate that existing guidelines are suitable for multigenerational whole mtGenome analysis, we recommend caution be taken when classifying heteroplasmic changes as differences for human identification. Our data showed the classification of heteroplasmic changes as differences increases the prevalence of inconclusive identification by 6%, with false exclusions observed in 0.34% of pairs examined. Further studies of multigenerational pedigrees, however, are needed to validate mtGenome interpretation guidelines for historical case work to more fully utilize emerging advancements.

Pedigree derived mutation rate across the entire mitochondrial genome of the Norfolk Island population

Estimates of mutation rates for various regions of the human mitochondrial genome (mtGenome) vary widely, depending on whether they are inferred using a phylogenetic approach or obtained directly from pedigrees. Traditionally, only the control region, or small portions of the coding region have been targeted for analysis due to the cost and effort required to produce whole mtGenome Sanger profiles. Here, we report one of the first pedigree derived mutation rates for the entire human mtGenome. The entire mtGenome from 225 individuals originating from Norfolk Island was analysed to estimate the pedigree derived mutation rate and compared against published mutation rates. These individuals were from 45 maternal lineages spanning 345 generational events. Mutation rates for various portions of the mtGenome were calculated. Nine mutations (including two transitions and seven cases of heteroplasmy) were observed, resulting in a rate of 0.058 mutations/site/million years (95% CI 0.031-0.108). These mutation rates are approximately 16 times higher than estimates derived from phylogenetic analysis with heteroplasmy detected in 13 samples (n = 225, 5.8% individuals). Providing one of the first pedigree derived estimates for the entire mtGenome, this study provides a better understanding of human mtGenome evolution and has relevance to many research fields, including medicine, anthropology and forensics.

Multi-phenotype genome-wide association studies of the Norfolk Island isolate implicate pleiotropic loci involved in chronic kidney disease

Chronic kidney disease (CKD) is a persistent impairment of kidney function. Genome-wide association studies (GWAS) have revealed multiple genetic loci associated with CKD susceptibility but the complete genetic basis is not yet clear. Since CKD shares risk factors with cardiovascular diseases and diabetes, there may be pleiotropic loci at play but may go undetected when using single phenotype GWAS. Here, we used multi-phenotype GWAS in the Norfolk Island isolate (n = 380) to identify new loci associated with CKD. We performed a principal components analysis on different combinations of 29 quantitative traits to extract principal components (PCs) representative of multiple correlated phenotypes. GWAS of a PC derived from glomerular filtration rate, serum creatinine, and serum urea identified a suggestive peak (pmin = 1.67 × 10-7) that mapped to KCNIP4. Inclusion of other secondary CKD measurements with these three kidney function traits identified the KCNIP4 locus with GWAS significance (pmin = 1.59 × 10-9). Finally, we identified a group of two SNPs with increased minor allele frequencies as potential functional variants. With the use of genetic isolate and the PCA-based multi-phenotype GWAS approach, we have revealed a potential pleotropic effect locus for CKD. Further studies are required to assess functional relevance of this locus.

Identification of Distant Regulatory Elements Using Expression Quantitative Trait Loci Mapping for Heat-Responsive Genes in Oysters

Many marine ectotherms, especially those inhabiting highly variable intertidal zones, develop high phenotypic plasticity in response to rapid climate change by modulating gene expression levels. Herein, we examined the regulatory architecture of heat-responsive gene expression plasticity in oysters using expression quantitative trait loci (eQTL) analysis. Using a backcross family of Crassostrea gigas and its sister species Crassostrea angulata under acute stress, 56 distant regulatory regions accounting for 6–26.6% of the gene expression variation were identified for 19 heat-responsive genes. In total, 831 genes and 164 single nucleotide polymorphisms (SNPs) that could potentially regulate expression of the target genes were screened in the eQTL region. The association between three SNPs and the corresponding target genes was verified in an independent family. Specifically, Marker13973 was identified for heat shock protein (HSP) family A member 9 (HspA9). Ribosomal protein L10a (RPL10A) was detected approximately 2 kb downstream of the distant regulatory SNP. Further, Marker14346-48 and Marker14346-85 were in complete linkage disequilibrium and identified for autophagy-related gene 7 (ATG7). Nuclear respiratory factor 1 (NRF1) was detected approximately 3 kb upstream of the two SNPs. These results suggested regulatory relationships between RPL10A and HSPA9 and between NRF1 and ATG7. Our findings indicate that distant regulatory mutations play an important role in the regulation of gene expression plasticity by altering upstream regulatory factors in response to heat stress. The identified eQTLs provide candidate biomarkers for predicting the persistence of oysters under future climate change scenarios.

A genome-wide methylation study of body fat traits in the Norfolk Island isolate

BACKGROUND AND AIMS

Natural variation in body fat is explained by both genetic and environmental effects. Epigenetic mechanisms such as DNA methylation can mediate these effects causing changes in gene expression leading to onset of obesity. Studies of genetic isolates have the potential to provide new epigenetic insights with advantages such as reduced genetic diversity and environmental exposures.

METHODS AND RESULTS

This was an exploratory study of genome-wide DNA methylation in relation to body fat traits in 47 healthy adults from the genetic isolate of Norfolk Island. Quantitative body fat traits (body fat percentage, body mass index, hip circumference, waist circumference, waist-hip-ratio and weight) were carefully measured. DNA methylation data was obtained from peripheral blood using Illumina 450K arrays. Multi-trait analysis was performed using Principal Component Analysis (PCA). CpG by trait association testing was performed using stepwise linear regressions. Two components were identified that explained approximately 89% of the phenotypic variance. In total, 5 differential methylated positions (DMPs) were identified at genome-wide significance (P≤ 2.4 × 10^-7), which mapped to GOT2-CDH8, LYSMD3, HIBADH, ADGRD1 and EBF4 genes. Gene set enrichment analysis of 848 genes containing suggestive DMPs (P≤ 1.0 × 10^-4) implicated the Cadherin (28 genes, P_adj = 6.76 × 10^-7) and Wnt signaling pathways (38 genes, P_adj = 7.78 × 10^-6).

CONCLUSION

This study provides new insights into the epigenetically influenced genes and pathways underlying body fat variation in a healthy cohort and provides targets for consideration in future studies of obesity risk.

AKR1C2 and AKR1C3 expression in adipose tissue: Association with body fat distribution and regulatory variants

BACKGROUND

Changes in androgen dynamics within adipose tissue have been proposed as modulators of body fat accumulation. In this context, AKR1C2 likely plays a significant role by inactivating 5α-dihydrotestosterone.

AIM

To characterize AKR1C2 expression patterns across adipose depots and cell populations and to provide insight into the link with body fat distribution and genetic regulation.

METHODS

We used RNA sequencing data from severely obese patients to assess patterns of AKR1C2 and AKR1C3 expression in abdominal adipose tissue depots and cell fractions. We additionally used data from 856 women to assess AKR1C2 heritability and to link its expression in adipose tissue with body fat distribution. Further, we used public resources to study AKR1C2 genetic regulation as well as reference epigenome data for regulatory element profiling and functional interpretation of genetic data.

RESULTS

We found that mature adipocytes and adipocyte-committed adipocyte progenitor cells (APCs) had enriched expression of AKR1C2. We found adipose tissue AKR1C2 and AKR1C3 expression to be significantly and positively associated with percentage trunk fat mass in women. We identified strong genetic regulation of AKR1C2 by rs28571848 and rs34477787 located on the binding sites of two nuclear transcription factors, namely retinoid acid-related orphan receptor alpha and the glucocorticoid receptor.

CONCLUSION

We confirm the link between AKR1C2, adipogenic differentiation and adipose tissue distribution. We provide insight into genetic regulation of AKR1C2 by identifying regulatory variants mapping to binding sites for the glucocorticoid receptor and retinoid acid-related orphan receptor alpha which may in part mediate the effect of AKR1C2 expression on body fat distribution.

FN3K expression in COPD: a potential comorbidity factor for cardiovascular disease

INTRODUCTION

Cigarette smoking and oxidative stress are common risk factors for the multi-morbidities associated with chronic obstructive pulmonary disease (COPD). Elevated levels of advanced glycation endproducts (AGE) increase the risk of cardiovascular disease (CVD) comorbidity and mortality. The enzyme fructosamine-3-kinase (FN3K) reduces this risk by lowering AGE levels.

METHODS

The distribution and expression of FN3K protein in lung tissues from stable COPD and control subjects, as well as an animal model of COPD, was assessed by immunohistochemistry. Serum FN3K protein and AGE levels were assessed by ELISA in patients with COPD exacerbations receiving metformin. Genetic variants within the FN3K and FN3K-RP genes were evaluated for associations with cardiorespiratory function in the Subpopulations and Intermediate Outcome Measures in COPD Study cohort.

RESULTS

This pilot study demonstrates that FN3K expression in the blood and human lung epithelium is distributed at either high or low levels irrespective of disease status. The percentage of lung epithelial cells expressing FN3K was higher in control smokers with normal lung function, but this induction was not observed in COPD patients nor in a smoking model of COPD. The top five nominal FN3K polymorphisms with possible association to decreased cardiorespiratory function (p<0.008-0.02), all failed to reach the threshold (p<0.0028) to be considered highly significant following multi-comparison analysis. Metformin enhanced systemic levels of FN3K in COPD subjects independent of their high-expression or low-expression status.

DISCUSSION

The data highlight that low and high FN3K expressors exist within our study cohort and metformin induces FN3K levels, highlighting a potential mechanism to reduce the risk of CVD comorbidity and mortality.

Transcriptome-wide Profiling of Cerebral Cavernous Malformations Patients Reveal Important Long noncoding RNA molecular signatures

Cerebral cavernous malformations (CCMs) are low-flow vascular malformations in the brain associated with recurrent hemorrhage and seizures. The current treatment of CCMs relies solely on surgical intervention. Henceforth, alternative non-invasive therapies are urgently needed to help prevent subsequent hemorrhagic episodes. Long non-coding RNAs (lncRNAs) belong to the class of non-coding RNAs and are known to regulate gene transcription and involved in chromatin remodeling via various mechanism. Despite accumulating evidence demonstrating the role of lncRNAs in cerebrovascular disorders, their identification in CCMs pathology remains unknown. The objective of the current study was to identify lncRNAs associated with CCMs pathogenesis using patient cohorts having 10 CCM patients and 4 controls from brain. Executing next generation sequencing, we performed whole transcriptome sequencing (RNA-seq) analysis and identified 1,967 lncRNAs and 4,928 protein coding genes (PCGs) to be differentially expressed in CCMs patients. Among these, we selected top 6 differentially expressed lncRNAs each having significant correlative expression with more than 100 differentially expressed PCGs. The differential expression status of the top lncRNAs, SMIM25 and LBX2-AS1 in CCMs was further confirmed by qRT-PCR analysis. Additionally, gene set enrichment analysis of correlated PCGs revealed critical pathways related to vascular signaling and important biological processes relevant to CCMs pathophysiology. Here, by transcriptome-wide approach we demonstrate that lncRNAs are prevalent in CCMs disease and are likely to play critical roles in regulating important signaling pathways involved in the disease progression. We believe, that detailed future investigations on this set of identified lncRNAs can provide useful insights into the biology and, ultimately, contribute in preventing this debilitating disease.

Ion torrent high throughput mitochondrial genome sequencing (HTMGS)

The implementation and popularity of next generation sequencing (NGS) has led to the development of various rapid whole mitochondrial genome sequencing techniques. We summarise an efficient and cost-effective NGS approach for mitochondrial genomic DNA in humans using the Ion Torrent platform, and further discuss our bioinformatics pipeline for streamlined variant calling. Ion 316 chips were utilised with the Ion Torrent semi-conductor platform Personal Genome Machine (PGM) to perform tandem sequencing of mitochondrial genomes from the core pedigree (n = 315) of the Norfolk Island Health Study. Key improvements from commercial methods focus on the initial PCR step, which currently requires extensive optimisation to ensure the accurate and reproducible elongation of each section of the complete mitochondrial genome. Dual-platform barcodes were incorporated into our protocol thereby extending its potential application onto Illumina-based systems. Our bioinformatics pipeline consists of a modified version of GATK best practices tailored for mitochondrial genomic data. When compared with current commercial methods, our method, termed high throughput mitochondrial genome sequencing (HTMGS), allows high multiplexing of samples and the use of alternate library preparation reagents at a lower cost per sample (~1.7 times) when compared to current commercial methodologies. Our HTMGS methodology also provides robust mitochondrial sequencing data (>450X average coverage) that can be applied and modified to suit various study designs. On average, we were able to identify ~30 variants per sample with 572 variants observed across 315 samples. We have developed a high throughput sequencing and analysis method targeting complete mitochondrial genomes; with the potential to be platform agnostic with analysis options that adhere to current best practices.

Genome-wide allele-specific methylation is enriched at gene regulatory regions in a multi-generation pedigree from the Norfolk Island isolate

Background

Allele-specific methylation (ASM) occurs when DNA methylation patterns exhibit asymmetry among alleles. ASM occurs at imprinted loci, but its presence elsewhere across the human genome is indicative of wider importance in terms of gene regulation and disease risk. Here, we studied ASM by focusing on blood-based DNA collected from 24 subjects comprising a 3-generation pedigree from the Norfolk Island genetic isolate. We applied a genome-wide bisulphite sequencing approach with a genotype-independent ASM calling method to map ASM across the genome. Regions of ASM were then tested for enrichment at gene regulatory regions using Genomic Association Test (GAT) tool.

Results

In total, we identified 1.12 M CpGs of which 147,170 (13%) exhibited ASM (P ≤ 0.05). When including contiguous ASM signal spanning ≥ 2 CpGs, this condensed to 12,761 ASM regions (AMRs). These AMRs tagged 79% of known imprinting regions and most (98.1%) co-localised with known single nucleotide variants. Notably, miRNA and lncRNA showed a 3.3- and 1.8-fold enrichment of AMRs, respectively (P < 0.005). Also, the 5′ UTR and start codons each showed a 3.5-fold enrichment of AMRs (P < 0.005). There was also enrichment of AMRs observed at subtelomeric regions of many chromosomes. Five out of 11 large AMRs localised to the protocadherin cluster on chromosome 5.

Conclusions

This study shows ASM extends far beyond genomic imprinting in humans and that gene regulatory regions are hotspots for ASM. Future studies of ASM in pedigrees should help to clarify transgenerational inheritance patterns in relation to genotype and disease phenotypes.

Influence of tissue context on gene prioritization for predicted transcriptome-wide association studies

Transcriptome-wide association studies (TWAS) have recently gained great attention due to their ability to prioritize complex trait-associated genes and promote potential therapeutics development for complex human diseases. TWAS integrates genotypic data with expression quantitative trait loci (eQTLs) to predict genetically regulated gene expression components and associates predictions with a trait of interest. As such, TWAS can prioritize genes whose differential expressions contribute to the trait of interest and provide mechanistic explanation of complex trait(s). Tissue-specific eQTL information grants TWAS the ability to perform association analysis on tissues whose gene expression profiles are otherwise hard to obtain, such as liver and heart. However, as eQTLs are tissue context-dependent, whether and how the tissue-specificity of eQTLs influences TWAS gene prioritization has not been fully investigated. In this study, we addressed this question by adopting two distinct TWAS methods, PrediXcan and UTMOST, which assume single tissue and integrative tissue effects of eQTLs, respectively. Thirty-eight baseline laboratory traits in 4,360 antiretroviral treatment-naïve individuals from the AIDS Clinical Trials Group (ACTG) studies comprised the input dataset for TWAS. We performed TWAS in a tissue-specific manner and obtained a total of 430 significant gene-trait associations (q-value < 0.05) across multiple tissues. Single tissue-based analysis by PrediXcan contributed 116 of the 430 associations including 64 unique gene-trait pairs in 28 tissues. Integrative tissue-based analysis by UTMOST found the other 314 significant associations that include 50 unique gene-trait pairs across all 44 tissues. Both analyses were able to replicate some associations identified in past variant-based genome-wide association studies (GWAS), such as high-density lipoprotein (HDL) and CETP (PrediXcan, q-value = 3.2e-16). Both analyses also identified novel associations. Moreover, single tissue-based and integrative tissuebased analysis shared 11 of 103 unique gene-trait pairs, for example, PSRC1-low-density lipoprotein (PrediXcan's lowest q-value = 8.5e-06; UTMOST's lowest q-value = 1.8e-05). This study suggests that single tissue-based analysis may have performed better at discovering gene-trait associations when combining results from all tissues. Integrative tissue-based analysis was better at prioritizing genes in multiple tissues and in trait-related tissue. Additional exploration is needed to confirm this conclusion. Finally, although single tissue-based and integrative tissue-based analysis shared significant novel discoveries, tissue context-dependency of eQTLs impacted TWAS gene prioritization. This study provides preliminary data to support continued work on tissue contextdependency of eQTL studies and TWAS.

Expression QTL analysis of glaucoma endophenotypes in the Norfolk Island isolate provides evidence that immune-related genes are associated with optic disc size

Primary open-angle glaucoma (POAG) is influenced by both genetic and environmental factors. Despite significant progress in identifying genetic variants associated with POAG, there remains a substantial amount of unexplained heritability. Study design features that may enhance knowledge of the genetic architecture include focusing on multiple quantitative traits related to ocular disorders (i.e. endophenotypes), targeting genetic variants that directly influence gene expression (i.e. cis-eQTLs) and utilising genetically isolated populations to reduce genetic and environmental noise and thus enhance association signals. In this study we performed heritability and blood-based eQTL association analysis of five key POAG endophenotypes in 330 individuals from the Norfolk Island (NI) isolate. Results showed evidence of heritability for all five traits, with H² estimates ranging from 0.35 for intraocular pressure (IOP) to 0.82 for central corneal thickness (CCT) (P < 0.05). The primary finding was for BTN3A2, whereby both cis-SNP and transcript were significantly associated with disc size within a conditional regression model. Specifically, this model included rs853676 (β = 0.23,P = 0.008) and transcript (β = 0.23, P = 0.03). We also observed a cis-SNP association between optic disc size and LPCAT2 independent of transcript (P = 0.0004). These genes have specific functions in immune system pathways and suggest a role for an inherited immune component of POAG risk. This study also demonstrates an alternate approach to understanding the functional genetic basis of POAG and ocular health more generally.

Genome-wide linkage and association analysis of primary open-angle glaucoma endophenotypes in the Norfolk Island isolate

PURPOSE

Primary open-angle glaucoma (POAG) refers to a group of heterogeneous diseases involving optic nerve damage. Two well-established risk factors for POAG are elevated intraocular pressure (IOP) and a thinner central corneal thickness (CCT). These endophenotypes exhibit a high degree of heritability across populations. Large-scale genome-wide association studies (GWASs) of outbred populations have robustly implicated several susceptibility gene variants for both IOP and CCT. Despite this progress, a substantial amount of genetic variance remains unexplained. Population-specific variants that might be rare in outbred populations may also influence POAG endophenotypes. The Norfolk Island population is a founder-effect genetic isolate that has been well characterized for POAG endophenotypes. This population is therefore a suitable candidate for mapping new variants that influence these complex traits.

METHODS

Three hundred and thirty participants from the Norfolk Island Eye Study (NIES) core pedigree provided DNA. Ocular measurements of CCT and IOP were also taken for analysis. Heritability analyses and genome-wide linkage analyses of short tandem repeats (STRs) were conducted using SOLAR. Pedigree-based GWASs of single-nucleotide polymorphisms (SNPs) were performed using the GenABEL software.

RESULTS

CCT was the most heritable endophenotype in this cohort (h² = 0.77, p = 6×10^-6), while IOP showed a heritability of 0.39 (p = 0.008). A genome-wide linkage analysis of these POAG phenotypes identified a maximum logarithm of the odds (LOD) score of 1.9 for CCT on chromosome 20 (p = 0.0016) and 1.3 for IOP on chromosome 15 (p = 0.0072). The GWAS results revealed a study-wise significant association for IOP at rs790357, which is located within DLG2 on chr11q14.1 (p = 1.02×10^-7). DLG2 is involved in neuronal signaling and development, and while it has not previously been associated with IOP, it has been associated with myopia. An analysis of 12 known SNPs for IOP showed that rs12419342 in RAPSN on chromosome 11 was nominally associated in Norfolk Island (NI; p = 0.0021). For CCT, an analysis of 26 known SNPs showed rs9938149 in BANP-ZNF469 on chromosome 16 was nominally associated in NI (p = 0.002).

CONCLUSIONS

These study results indicate that CCT and IOP exhibit a substantial degree of heritability in the NI pedigree, indicating a genetic component. A genome-wide linkage analysis of POAG endophenotypes did not reveal any major effect loci, but the GWASs did implicate several known loci, as well as a potential new locus in DLG2, suggesting a role for neuronal signaling in development in IOP and perhaps POAG. These results also highlight the need to target rarer variants via whole genome sequencing in this genetic isolate.

Methylome-wide association study of whole blood DNA in the Norfolk Island isolate identifies robust loci associated with age

Epigenetic regulation of various genomic functions, including gene expression, provide mechanisms whereby an organism can dynamically respond to changes in its environment and modify gene expression accordingly. One epigenetic mechanism implicated in human aging and age-related disorders is DNA methylation. Isolated populations such as Norfolk Island (NI) should be advantageous for the identification of epigenetic factors related to aging due to reduced genetic and environmental variation. Here we conducted a methylome-wide association study of age using whole blood DNA in 24 healthy female individuals from the NI genetic isolate (aged 24-47 years). We analysed 450K methylation array data using a machine learning approach (GLMnet) to identify age-associated CpGs. We identified 497 CpG sites, mapping to 422 genes, associated with age, with 11 sites previously associated with age. The strongest associations identified were for a single CpG site in MYOF and an extended region within the promoter of DDO. These hits were validated in curated public data from 2316 blood samples (MARMAL-AID). This study is the first to report robust age associations for MYOF and DDO, both of which have plausible functional roles in aging. This study also illustrates the value of genetic isolates to reveal new associations with epigenome-level data.

Four Susceptibility Loci for Gallstone Disease Identified in a Meta-analysis of Genome-Wide Association Studies

BACKGROUND & AIMS

A genome-wide association study (GWAS) of 280 cases identified the hepatic cholesterol transporter ABCG8 as a locus associated with risk for gallstone disease, but findings have not been reported from any other GWAS of this phenotype. We performed a large-scale, meta-analysis of GWASs of individuals of European ancestry with available prior genotype data, to identify additional genetic risk factors for gallstone disease.

METHODS

We obtained per-allele odds ratio (OR) and standard error estimates using age- and sex-adjusted logistic regression models within each of the 10 discovery studies (8720 cases and 55,152 controls). We performed an inverse variance weighted, fixed-effects meta-analysis of study-specific estimates to identify single-nucleotide polymorphisms that were associated independently with gallstone disease. Associations were replicated in 6489 cases and 62,797 controls.

RESULTS

We observed independent associations for 2 single-nucleotide polymorphisms at the ABCG8 locus: rs11887534 (OR, 1.69; 95% confidence interval [CI], 1.54-1.86; P = 2.44 × 10(-60)) and rs4245791 (OR, 1.27; P = 1.90 × 10(-34)). We also identified and/or replicated associations for rs9843304 in TM4SF4 (OR, 1.12; 95% CI, 1.08-1.16; P = 6.09 × 10(-11)), rs2547231 in SULT2A1 (encodes a sulfoconjugation enzyme that acts on hydroxysteroids and cholesterol-derived sterol bile acids) (OR, 1.17; 95% CI, 1.12-1.21; P = 2.24 × 10(-10)), rs1260326 in glucokinase regulatory protein (OR, 1.12; 95% CI, 1.07-1.17; P = 2.55 × 10(-10)), and rs6471717 near CYP7A1 (encodes an enzyme that catalyzes conversion of cholesterol to primary bile acids) (OR, 1.11; 95% CI, 1.08-1.15; P = 8.84 × 10(-9)). Among individuals of African American and Hispanic American ancestry, rs11887534 and rs4245791 were associated positively with gallstone disease risk, whereas the association for the rs1260326 variant was inverse.

CONCLUSIONS

In this large-scale GWAS of gallstone disease, we identified 4 loci in genes that have putative functions in cholesterol metabolism and transport, and sulfonylation of bile acids or hydroxysteroids.

Choline metabolites: gene by diet interactions

PURPOSE OF REVIEW

The review highlights recent advances in our understanding of the interactions between genetic polymorphisms in genes that metabolize choline and the dietary requirements of choline and how these interactions relate to human health and disease.

RECENT FINDINGS

The importance of choline as an essential nutrient has been well established, but our appreciation of the interaction between our underlying genetic architecture and dietary choline requirements is only beginning. It has been shown in both human and animal studies that choline deficiencies contribute to diseases such as nonalcoholic fatty liver disease and various neurodegenerative diseases. An adequate supply of dietary choline is important for optimum development, highlighted by the increased maternal requirements during fetal development and in breast-fed infants. We discuss recent studies investigating variants in PEMT and MTHFR1 that are associated with a variety of birth defects. In addition to genetic interactions, we discuss several recent studies that uncover changes in fetal global methylation patterns in response to maternal dietary choline intake that result in changes in gene expression in the offspring. In contrast to the developmental role of adequate choline, there is now an appreciation of the role choline has in cardiovascular disease through the gut microbiota-mediated metabolite trimethylamine N-oxide. This pathway highlights some of our understanding of how the microbiome affects nutrient processing and bioavailability. Finally, to better characterize the genetic architecture regulating choline requirements, we discuss recent results focused on identifying polymorphisms that regulate choline and its derivative products.

SUMMARY

Here we discuss recent studies that have advanced our understanding of how specific alleles in key choline metabolism genes are related to dietary choline requirements and human disease.

Serum bilirubin concentration is modified by UGT1A1 haplotypes and influences risk of type-2 diabetes in the Norfolk Island genetic isolate

Background

Located in the Pacific Ocean between Australia and New Zealand, the unique population isolate of Norfolk Island has been shown to exhibit increased prevalence of metabolic disorders (type-2 diabetes, cardiovascular disease) compared to mainland Australia. We investigated this well-established genetic isolate, utilising its unique genomic structure to increase the ability to detect related genetic markers. A pedigree-based genome-wide association study of 16 routinely collected blood-based clinical traits in 382 Norfolk Island individuals was performed.

Results

A striking association peak was located at chromosome 2q37.1 for both total bilirubin and direct bilirubin, with 29 SNPs reaching statistical significance (P < 1.84 × 10⁻⁷). Strong linkage disequilibrium was observed across a 200 kb region spanning the UDP-glucuronosyltransferase family, including UGT1A1, an enzyme known to metabolise bilirubin. Given the epidemiological literature suggesting negative association between CVD-risk and serum bilirubin we further explored potential associations using stepwise multivariate regression, revealing significant association between direct bilirubin concentration and type-2 diabetes risk. In the Norfolk Island cohort increased direct bilirubin was associated with a 28 % reduction in type-2 diabetes risk (OR: 0.72, 95 % CI: 0.57-0.91, P = 0.005). When adjusted for genotypic effects the overall model was validated, with the adjusted model predicting a 30 % reduction in type-2 diabetes risk with increasing direct bilirubin concentrations (OR: 0.70, 95 % CI: 0.53-0.89, P = 0.0001).

Conclusions

In summary, a pedigree-based GWAS of blood-based clinical traits in the Norfolk Island population has identified variants within the UDPGT family directly associated with serum bilirubin levels, which is in turn implicated with reduced risk of developing type-2 diabetes within this population.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0291-z) contains supplementary material, which is available to authorized users.

A Phenomic Scan of the Norfolk Island Genetic Isolate Identifies a Major Pleiotropic Effect Locus Associated with Metabolic and Renal Disorder Markers

Multiphenotype genome-wide association studies (GWAS) may reveal pleiotropic genes, which would remain undetected using single phenotype analyses. Analysis of large pedigrees offers the added advantage of more accurately assessing trait heritability, which can help prioritise genetically influenced phenotypes for GWAS analysis. In this study we performed a principal component analysis (PCA), heritability (h²) estimation and pedigree-based GWAS of 37 cardiovascular disease -related phenotypes in 330 related individuals forming a large pedigree from the Norfolk Island genetic isolate. PCA revealed 13 components explaining >75% of the total variance. Nine components yielded statistically significant h² values ranging from 0.22 to 0.54 (P<0.05). The most heritable component was loaded with 7 phenotypic measures reflecting metabolic and renal dysfunction. A GWAS of this composite phenotype revealed statistically significant associations for 3 adjacent SNPs on chromosome 1p22.2 (P<1x10^-8). These SNPs form a 42kb haplotype block and explain 11% of the genetic variance for this renal function phenotype. Replication analysis of the tagging SNP (rs1396315) in an independent US cohort supports the association (P = 0.000011). Blood transcript analysis showed 35 genes were associated with rs1396315 (P<0.05). Gene set enrichment analysis of these genes revealed the most enriched pathway was purine metabolism (P = 0.0015). Overall, our findings provide convincing evidence for a major pleiotropic effect locus on chromosome 1p22.2 influencing risk of renal dysfunction via purine metabolism pathways in the Norfolk Island population. Further studies are now warranted to interrogate the functional relevance of this locus in terms of renal pathology and cardiovascular disease risk.

While many large genetic association studies have identified genes playing a role in complex disorders, there is still concern over the amount of missing genetic heritability. With this in mind, we have used a data reduction approach alongside pedigree-based association to obtain highly heritable components which explain 'hidden' variance of multiphenotypes within a large pedigree from the Norfolk Island genetic isolate. The most heritable of these components involved 7 traits reflecting metabolic and renal functionality, association of which locates to an intergenic region on chromosome 1p22.2. By integrating gene expression information, we identified enrichment of a purine metabolism pathway, further strengthening the metabolic nature of the observed association. Adding additional support to our approach, we show association of the tagging SNP (rs1396315) in an independent US population. The findings presented here are of particular interest as they implicate pleiotropic effect loci and newly associated biological pathways underlying cardiovascular disease risk.

'Mutiny on the Bounty': the genetic history of Norfolk Island reveals extreme gender-biased admixture

Background

The Pacific Oceania region was one of the last regions of the world to be settled via human migration. Here we outline a settlement of this region that has given rise to a uniquely admixed population. The current Norfolk Island population has arisen from a small number of founders with mixed Caucasian and Polynesian ancestry, descendants of a famous historical event. The ‘Mutiny on the Bounty’ has been told in history books, songs and the big screen, but recently this story can be portrayed through comprehensive molecular genetics. Written history details betrayal and murder leading to the founding of Pitcairn Island by European mutineers and the Polynesian women who left Tahiti with them. Investigation of detailed genealogical records supports historical accounts.

Findings

Using genetics, we show distinct maternal Polynesian mitochondrial lineages in the present day population, as well as a European centric Y-chromosome phylogeny. These results comprehensively characterise the unique gender-biased admixture of this genetic isolate and further support the historical records relating to Norfolk Island.

Conclusions

Our results significantly refine previous population genetic studies investigating Polynesian versus Caucasian diversity in the Norfolk Island population and add information that is beneficial to future disease and gene mapping studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s13323-015-0028-9) contains supplementary material, which is available to authorized users.

Common polygenic variation contributes to risk of migraine in the Norfolk Island population

Migraine has been defined as a common disabling primary headache disorder. Epidemiology studies have provided with the undeniable evidence of genetic components as active players in the development of the disease under a polygenic model in which multiple risk alleles exert modest individual effects. Our objective was to test the contribution of a polygenic effect to migraine risk in the Norfolk Island population using a panel of SNPs reported to be disease associated in published migraine GWAS. We also investigated whether individual SNPs were associated with gene expression levels measured in whole blood. Polygenic scores were calculated in a total of 285 related individuals (74 cases, 211 controls) from the Norfolk Island using 51 SNPs previously reported to be associated with migraine in published GWAS. The association between polygenic score and migraine case-control status was tested using logistic regression. Results indicate that a migraine polygenic risk score was associated with migraine case-control status in this population (P = 0.016). This supports the hypothesis that multiple SNPs with weak effects collectively contribute to migraine risk in this population. Amongst the SNPs included in the polygenic model, four were associated with the expression of the USMG5 gene, including rs171251 (P = 0.012). Results from this study provide evidence for a polygenic contribution to migraine risk in an isolated population and highlight specific SNPs that regulate the expression of USMG5, a gene critical for mitochondrial function.

A genome-wide expression quantitative trait loci analysis of proprotein convertase subtilisin/kexin enzymes identifies a novel regulatory gene variant for FURIN expression and blood pressure

Proprotein convertase subtilisin/kexin (PCSK) enzymes cleave and convert their immature substrates into biologically active forms. Polymorphisms in the PCSK genes have been reported to associate with human diseases and phenotypes, including hypercholesterolemia and blood pressure (BP), and targeting PCSKs is considered a promising future form of drug therapy. PCSK processing is readily induced upon upregulation of the enzyme, but the genetic factors contributing to PCSK expression have not been thoroughly characterized. To gain a comprehensive understanding of the genetic regulation of PCSK expression, we performed, for the first time, a genome-wide expression quantitative trait loci (eQTL) analysis using mRNA expression in >1400 human peripheral blood samples from the Cardiovascular Risk in Young Finns Study and ca. ten million single-nucleotide polymorphisms (SNPs). The expression data showed clear expression for FURIN, PCSK5, PCSK7 and MBTPS1 (membrane-bound transcription factor peptidase, site 1) mRNAs in virtually all tested samples. A discovery analysis demonstrated a genome-wide significant (p < 5 × 10(-8)) association with the selected PCSK probes for 1024 variants, which were located at ten independent loci. Of these loci, 5/10 could be confirmed to regulate PCSK expression in two additional and independent sample sets. Finally, a phenotypic analysis demonstrated that a novel cis-eQTL SNP rs4702 for FURIN is strongly associated with both diastolic (p = 0.012) and systolic (p = 0.035) BP levels, as well as peripheral vascular resistance (p = 0.003). These findings indicate that the expression of the PCSK enzymes is regulated by genetic factors, which have biological roles in health and disease.

Effects of long-term averaging of quantitative blood pressure traits on the detection of genetic associations

Blood pressure (BP) is a heritable, quantitative trait with intraindividual variability and susceptibility to measurement error. Genetic studies of BP generally use single-visit measurements and thus cannot remove variability occurring over months or years. We leveraged the idea that averaging BP measured across time would improve phenotypic accuracy and thereby increase statistical power to detect genetic associations. We studied systolic BP (SBP), diastolic BP (DBP), mean arterial pressure (MAP), and pulse pressure (PP) averaged over multiple years in 46,629 individuals of European ancestry. We identified 39 trait-variant associations across 19 independent loci (p < 5 × 10(-8)); five associations (in four loci) uniquely identified by our LTA analyses included those of SBP and MAP at 2p23 (rs1275988, near KCNK3), DBP at 2q11.2 (rs7599598, in FER1L5), and PP at 6p21 (rs10948071, near CRIP3) and 7p13 (rs2949837, near IGFBP3). Replication analyses conducted in cohorts with single-visit BP data showed positive replication of associations and a nominal association (p < 0.05). We estimated a 20% gain in statistical power with long-term average (LTA) as compared to single-visit BP association studies. Using LTA analysis, we identified genetic loci influencing BP. LTA might be one way of increasing the power of genetic associations for continuous traits in extant samples for other phenotypes that are measured serially over time.