Principal component and linkage analysis of cardiovascular risk traits in the Norfolk isolate

Determining independence and associations among various cardiovascular disease risk factors in 9-12 years old school-children: a cross sectional study

BACKGROUND

Cardiovascular disease (CVD) risk assessment of children typically includes evaluating multiple CVD risk factors some of which tend to correlate each other. However, in older children and young adolescents, there are little data on the level of independence of CVD risk factors. The purpose of this study was to examine the relationships among various CVD risk factors to determine the level of independence of each risk factor in a sample of 5th-grade public school students.

METHOD

A cross-sectional analysis of 1525 children (856 girls and 669 boys; age: 9-12 years) who participated in baseline CVD risk assessment for the (S)Partners for Heart Health program from 2010 - 2018. Thirteen CVD risk factor variables were used in the analysis and included blood lipids [low-density lipoprotein (LDL), high-density lipoprotein (HDL), total cholesterol (TC), and triglycerides], resting systolic and diastolic blood pressure (BP); anthropometrics [height, weight, body mass index (BMI), % body fat, waist circumference (WC)]. Additionally, acanthosis nigricans (a marker insulin resistance and diabetes), and cardiorespiratory fitness (VO2 ml/kg) was estimated using the PACER. Descriptive statistics, bivariate Pearson correlations, and principal component analysis were used to determine the relationships among these variables and the independence.

RESULTS

Parallel analysis indicated two components should be extracted. Among the two components extracted, WC, % body fat, and BMI loaded highest on component 1, which explained 34% of the total variance. Systolic BP and diastolic BP loaded predominantly on component 2 and accounted for 17% of the variance. Cardiorespiratory fitness, acanthosis nigricans, HDL, and triglycerides loaded highest on the first component (loadings between 0.42 and 0.57) but still suggest some non-shared variance with this component. Low-density lipoprotein had low loadings on each component. Factor loadings were stable across sex.

CONCLUSION

Among the various CVD risk indicators, measures of adiposity loaded highest on the component that explained the largest proportion of variability in the data reinforcing the importance of assessing adiposity in CVD risk assessment. In addition, blood pressure loaded highest on the second component, suggesting their relative independence when assessing CVD risk. The data also provide support and rationale for determining what CVD risk factors to include- based on resource needs. For example, researchers or public health programs may choose to assess WC instead of lipid profile for cardiovascular related problems if ease of assessment and cost are considerations.

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.

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 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.

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.

Principal component analysis of cardiovascular risk traits in three generations cohort among Indian Punjabi population

The current study focused to determine significant cardiovascular risk factors through principal component factor analysis (PCFA) among three generations on 1827 individuals in three generations including 911 males (378 from offspring, 439 from parental and 94 from grand-parental generations) and 916 females (261 from offspring, 515 from parental and 140 from grandparental generations). The study performed PCFA with orthogonal rotation to reduce 12 inter-correlated variables into groups of independent factors. The factors have been identified as 2 for male grandparents, 3 for male offspring, female parents and female grandparents each, 4 for male parents and 5 for female offspring. This data reduction method identified these factors that explained 72%, 84%, 79%, 69%, 70% and 73% for male and female offspring, male and female parents and male and female grandparents respectively, of the variations in original quantitative traits. The factor 1 accounting for the largest portion of variations was strongly loaded with factors related to obesity (body mass index (BMI), waist circumference (WC), waist to hip ratio (WHR), and thickness of skinfolds) among all generations with both sexes, which has been known to be an independent predictor for cardiovascular morbidity and mortality. The second largest components, factor 2 and factor 3 for almost all generations reflected traits of blood pressure phenotypes loaded, however, in male offspring generation it was observed that factor 2 was loaded with blood pressure phenotypes as well as obesity. This study not only confirmed but also extended prior work by developing a cumulative risk scale from factor scores. Till today, such a cumulative and extensive scale has not been used in any Indian studies with individuals of three generations. These findings and study highlight the importance of global approach for assessing the risk and need for studies that elucidate how these different cardiovascular risk factors interact with each other over the time to create clinical disease. The findings also added depth to the negligible amount of literature of factor analysis of cardiovascular risk in any Indian ethnic population.

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

Cardiovascular disease (CVD) affects millions of people worldwide and is influenced by numerous factors, including lifestyle and genetics. Expression quantitative trait loci (eQTLs) influence gene expression and are good candidates for CVD risk. Founder-effect pedigrees can provide additional power to map genes associated with disease risk. Therefore, we identified eQTLs in the genetic isolate of Norfolk Island (NI) and tested for associations between these and CVD risk factors. We measured genome-wide transcript levels of blood lymphocytes in 330 individuals and used pedigree-based heritability analysis to identify heritable transcripts. eQTLs were identified by genome-wide association testing of these transcripts. Testing for association between CVD risk factors (i.e., blood lipids, blood pressure, and body fat indices) and eQTLs revealed 1,712 heritable transcripts (p < 0.05) with heritability values ranging from 0.18 to 0.84. From these, we identified 200 cis-acting and 70 trans-acting eQTLs (p < 1.84 × 10(-7)) An eQTL-centric analysis of CVD risk traits revealed multiple associations, including 12 previously associated with CVD-related traits. Trait versus eQTL regression modeling identified four CVD risk candidates (NAAA, PAPSS1, NME1, and PRDX1), all of which have known biological roles in disease. In addition, we implicated several genes previously associated with CVD risk traits, including MTHFR and FN3KRP. We have successfully identified a panel of eQTLs in the NI pedigree and used this to implicate several genes in CVD risk. Future studies are required for further assessing the functional importance of these eQTLs and whether the findings here also relate to outbred populations.

The Norfolk Island Eye Study (NIES): rationale, methodology and distribution of ocular biometry (biometry of the bounty)

AIM

To describe the recruitment, ophthalmic examination methods and distribution of ocular biometry of participants in the Norfolk Island Eye Study, who were individuals descended from the English Bounty mutineers and their Polynesian wives.

METHODS

All 1,275 permanent residents of Norfolk Island aged over 15 years were invited to participate, including 602 individuals involved in a 2001 cardiovascular disease study. Participants completed a detailed questionnaire and underwent a comprehensive eye assessment including stereo disc and retinal photography, ocular coherence topography and conjunctival autofluorescence assessment. Additionally, blood or saliva was taken for DNA testing.

RESULTS

781 participants aged over 15 years were seen (54% female), comprising 61% of the permanent Island population. 343 people (43.9%) could trace their family history to the Pitcairn Islanders (Norfolk Island Pitcairn Pedigree). Mean anterior chamber depth was 3.32mm, mean axial length (AL) was 23.5mm, and mean central corneal thickness was 546 microns. There were no statistically significant differences in these characteristics between persons with and without Pitcairn Island ancestry. Mean intra-ocular pressure was lower in people with Pitcairn Island ancestry: 15.89mmHg compared to those without Pitcairn Island ancestry 16.49mmHg (P = .007). The mean keratometry value was lower in people with Pitcairn Island ancestry (43.22 vs. 43.52, P = .007). The corneas were flatter in people of Pitcairn ancestry but there was no corresponding difference in AL or refraction.

CONCLUSION

Our study population is highly representative of the permanent population of Norfolk Island. Ocular biometry was similar to that of other white populations. Heritability estimates, linkage analysis and genome-wide studies will further elucidate the genetic determinants of chronic ocular diseases in this genetic isolate.

European and Polynesian admixture in the Norfolk Island population

The Norfolk Island population in the South Pacific is primarily the product of recent admixture between a small number of British male and Polynesian female founders. We identified and genotyped 128 Ancestry Informative Markers (AIMs) spread across the autosomes, X/Y chromosomes and mitochondrial DNA genome, to explore and quantify the current levels of genetic admixture in the Norfolk Islanders. On the basis of autosomal AIMs, the population shows mean European and Polynesian ancestry proportions of 88 and 12%, respectively. However, there is a substantial variation between individuals ranging from total European ancestry to near total Polynesian origin. There is a strong correlation between individual genetic estimates of Polynesian ancestry and those derived from the extensive pedigree and genealogical records of Islanders. Also in line with historical accounts, there is a substantial asymmetry in the maternal and paternal origins of the Islanders with almost all Y-chromosomes of European origin whereas at least 25% of mtDNAs appear to have a Polynesian origin. Accurate knowledge of ancestry will be important in future attempts to use the Island population in admixture mapping approaches to find the genes that underlie differences in the risk to some diseases between Europeans and Polynesians.