Precision phenotyping, panomics, and system-level bioinformatics to delineate complex biologies of atherosclerosis: rationale and design of the "Genetic Loci and the Burden of Atherosclerotic Lesions" study

BACKGROUND

Complex biological networks of atherosclerosis are largely unknown.

OBJECTIVE

The main objective of the Genetic Loci and the Burden of Atherosclerotic Lesions study is to assemble comprehensive biological networks of atherosclerosis using advanced cardiovascular imaging for phenotyping, a panomic approach to identify underlying genomic, proteomic, metabolomic, and lipidomic underpinnings, analyzed by systems biology-driven bioinformatics.

METHODS

By design, this is a hypothesis-free unbiased discovery study collecting a large number of biologically related factors to examine biological associations between genomic, proteomic, metabolomic, lipidomic, and phenotypic factors of atherosclerosis. The Genetic Loci and the Burden of Atherosclerotic Lesions study (NCT01738828) is a prospective, multicenter, international observational study of atherosclerotic coronary artery disease. Approximately 7500 patients are enrolled and undergo non-contrast-enhanced coronary calcium scanning by CT for the detection and quantification of coronary artery calcium, as well as coronary artery CT angiography for the detection and quantification of plaque, stenosis, and overall coronary artery disease burden. In addition, patients undergo whole genome sequencing, DNA methylation, whole blood-based transcriptome sequencing, unbiased proteomics based on mass spectrometry, as well as metabolomics and lipidomics on a mass spectrometry platform. The study is analyzed in 3 subsequent phases, and each phase consists of a discovery cohort and an independent validation cohort. For the primary analysis, the primary phenotype will be the presence of any atherosclerotic plaque, as detected by cardiac CT. Additional phenotypic analyses will include per patient maximal luminal stenosis defined as 50% and 70% diameter stenosis. Single-omic and multi-omic associations will be examined for each phenotype; putative biomarkers will be assessed for association, calibration, discrimination, and reclassification.

Gene-centric meta-analyses for central adiposity traits in up to 57 412 individuals of European descent confirm known loci and reveal several novel associations

Waist circumference (WC) and waist-to-hip ratio (WHR) are surrogate measures of central adiposity that are associated with adverse cardiovascular events, type 2 diabetes and cancer independent of body mass index (BMI). WC and WHR are highly heritable with multiple susceptibility loci identified to date. We assessed the association between SNPs and BMI-adjusted WC and WHR and unadjusted WC in up to 57 412 individuals of European descent from 22 cohorts collaborating with the NHLBI's Candidate Gene Association Resource (CARe) project. The study population consisted of women and men aged 20-80 years. Study participants were genotyped using the ITMAT/Broad/CARE array, which includes ∼50 000 cosmopolitan tagged SNPs across ∼2100 cardiovascular-related genes. Each trait was modeled as a function of age, study site and principal components to control for population stratification, and we conducted a fixed-effects meta-analysis. No new loci for WC were observed. For WHR analyses, three novel loci were significantly associated (P < 2.4 × 10(-6)). Previously unreported rs2811337-G near TMCC1 was associated with increased WHR (β ± SE, 0.048 ± 0.008, P = 7.7 × 10(-9)) as was rs7302703-G in HOXC10 (β = 0.044 ± 0.008, P = 2.9 × 10(-7)) and rs936108-C in PEMT (β = 0.035 ± 0.007, P = 1.9 × 10(-6)). Sex-stratified analyses revealed two additional novel signals among females only, rs12076073-A in SHC1 (β = 0.10 ± 0.02, P = 1.9 × 10(-6)) and rs1037575-A in ATBDB4 (β = 0.046 ± 0.01, P = 2.2 × 10(-6)), supporting an already established sexual dimorphism of central adiposity-related genetic variants. Functional analysis using ENCODE and eQTL databases revealed that several of these loci are in regulatory regions or regions with differential expression in adipose tissue.

Gene-centric meta-analysis in 87,736 individuals of European ancestry identifies multiple blood-pressure-related loci

Blood pressure (BP) is a heritable risk factor for cardiovascular disease. To investigate genetic associations with systolic BP (SBP), diastolic BP (DBP), mean arterial pressure (MAP), and pulse pressure (PP), we genotyped ~50,000 SNPs in up to 87,736 individuals of European ancestry and combined these in a meta-analysis. We replicated findings in an independent set of 68,368 individuals of European ancestry. Our analyses identified 11 previously undescribed associations in independent loci containing 31 genes including PDE1A, HLA-DQB1, CDK6, PRKAG2, VCL, H19, NUCB2, RELA, HOXC@ complex, FBN1, and NFAT5 at the Bonferroni-corrected array-wide significance threshold (p < 6 × 10(-7)) and confirmed 27 previously reported associations. Bioinformatic analysis of the 11 loci provided support for a putative role in hypertension of several genes, such as CDK6 and NUCB2. Analysis of potential pharmacological targets in databases of small molecules showed that ten of the genes are predicted to be a target for small molecules. In summary, we identified previously unknown loci associated with BP. Our findings extend our understanding of genes involved in BP regulation, which may provide new targets for therapeutic intervention or drug response stratification.

The NIH Toolbox Cognition Battery: results from a large normative developmental sample (PING)

OBJECTIVE

The NIH Toolbox Cognition Battery (NTCB) was designed to provide a brief, efficient computerized test of key neuropsychological functions appropriate for use in children as young as 3 years of age. This report describes the performance of a large group of typically developing children and adolescents and examines the impact of age and sociocultural variables on test performance.

METHOD

The NTCB was administered to a sample of 1,020 typically developing males and females ranging in age from 3 to 20 years, diverse in terms of socioeconomic status (SES) and race/ethnicity, as part of the new publicly accessible Pediatric Imaging, Neurocognition, and Genetics (PING) data resource, at 9 sites across the United States.

RESULTS

General additive models of nonlinear age-functions were estimated from age-differences in test performance on the 8 NTCB subtests while controlling for family SES and genetic ancestry factors (GAFs). Age accounted for the majority of the variance across all NTCB scores, with additional significant contributions of gender on some measures, and of SES and race/ethnicity (GAFs) on all. After adjusting for age and gender, SES and GAFs explained a substantial proportion of the remaining unexplained variance in Picture Vocabulary scores.

CONCLUSIONS

The results highlight the sensitivity to developmental effects and efficiency of this new computerized assessment battery for neurodevelopmental research. Limitations are observed in the form of some ceiling effects in older children, some floor effects, particularly on executive function tests in the youngest participants, and evidence for variable measurement sensitivity to cultural/socioeconomic factors.

Genome-wide association study of shared components of reading disability and language impairment

Written and verbal languages are neurobehavioral traits vital to the development of communication skills. Unfortunately, disorders involving these traits-specifically reading disability (RD) and language impairment (LI)-are common and prevent affected individuals from developing adequate communication skills, leaving them at risk for adverse academic, socioeconomic and psychiatric outcomes. Both RD and LI are complex traits that frequently co-occur, leading us to hypothesize that these disorders share genetic etiologies. To test this, we performed a genome-wide association study on individuals affected with both RD and LI in the Avon Longitudinal Study of Parents and Children. The strongest associations were seen with markers in ZNF385D (OR = 1.81, P = 5.45 × 10(-7) ) and COL4A2 (OR = 1.71, P = 7.59 × 10(-7) ). Markers within NDST4 showed the strongest associations with LI individually (OR = 1.827, P = 1.40 × 10(-7) ). We replicated association of ZNF385D using receptive vocabulary measures in the Pediatric Imaging Neurocognitive Genetics study (P = 0.00245). We then used diffusion tensor imaging fiber tract volume data on 16 fiber tracts to examine the implications of replicated markers. ZNF385D was a predictor of overall fiber tract volumes in both hemispheres, as well as global brain volume. Here, we present evidence for ZNF385D as a candidate gene for RD and LI. The implication of transcription factor ZNF385D in RD and LI underscores the importance of transcriptional regulation in the development of higher order neurocognitive traits. Further study is necessary to discern target genes of ZNF385D and how it functions within neural development of fluent language.

Genome-wide association analysis of blood-pressure traits in African-ancestry individuals reveals common associated genes in African and non-African populations

High blood pressure (BP) is more prevalent and contributes to more severe manifestations of cardiovascular disease (CVD) in African Americans than in any other United States ethnic group. Several small African-ancestry (AA) BP genome-wide association studies (GWASs) have been published, but their findings have failed to replicate to date. We report on a large AA BP GWAS meta-analysis that includes 29,378 individuals from 19 discovery cohorts and subsequent replication in additional samples of AA (n = 10,386), European ancestry (EA) (n = 69,395), and East Asian ancestry (n = 19,601). Five loci (EVX1-HOXA, ULK4, RSPO3, PLEKHG1, and SOX6) reached genome-wide significance (p < 1.0 × 10(-8)) for either systolic or diastolic BP in a transethnic meta-analysis after correction for multiple testing. Three of these BP loci (EVX1-HOXA, RSPO3, and PLEKHG1) lack previous associations with BP. We also identified one independent signal in a known BP locus (SOX6) and provide evidence for fine mapping in four additional validated BP loci. We also demonstrate that validated EA BP GWAS loci, considered jointly, show significant effects in AA samples. Consequently, these findings suggest that BP loci might have universal effects across studied populations, demonstrating that multiethnic samples are an essential component in identifying, fine mapping, and understanding their trait variability.

Genome-wide association study of age at menarche in African-American women

African-American (AA) women have earlier menarche on average than women of European ancestry (EA), and earlier menarche is a risk factor for obesity and type 2 diabetes among other chronic diseases. Identification of common genetic variants associated with age at menarche has a potential value in pointing to the genetic pathways underlying chronic disease risk, yet comprehensive genome-wide studies of age at menarche are lacking for AA women. In this study, we tested the genome-wide association of self-reported age at menarche with common single-nucleotide polymorphisms (SNPs) in a total of 18 089 AA women in 15 studies using an additive genetic linear regression model, adjusting for year of birth and population stratification, followed by inverse-variance weighted meta-analysis (Stage 1). Top meta-analysis results were then tested in an independent sample of 2850 women (Stage 2). First, while no SNP passed the pre-specified P < 5 × 10(-8) threshold for significance in Stage 1, suggestive associations were found for variants near FLRT2 and PIK3R1, and conditional analysis identified two independent SNPs (rs339978 and rs980000) in or near RORA, strengthening the support for this suggestive locus identified in EA women. Secondly, an investigation of SNPs in 42 previously identified menarche loci in EA women demonstrated that 25 (60%) of them contained variants significantly associated with menarche in AA women. The findings provide the first evidence of cross-ethnic generalization of menarche loci identified to date, and suggest a number of novel biological links to menarche timing in AA women.

Loci influencing blood pressure identified using a cardiovascular gene-centric array

Blood pressure (BP) is a heritable determinant of risk for cardiovascular disease (CVD). To investigate genetic associations with systolic BP (SBP), diastolic BP (DBP), mean arterial pressure (MAP) and pulse pressure (PP), we genotyped ∼50 000 single-nucleotide polymorphisms (SNPs) that capture variation in ∼2100 candidate genes for cardiovascular phenotypes in 61 619 individuals of European ancestry from cohort studies in the USA and Europe. We identified novel associations between rs347591 and SBP (chromosome 3p25.3, in an intron of HRH1) and between rs2169137 and DBP (chromosome1q32.1 in an intron of MDM4) and between rs2014408 and SBP (chromosome 11p15 in an intron of SOX6), previously reported to be associated with MAP. We also confirmed 10 previously known loci associated with SBP, DBP, MAP or PP (ADRB1, ATP2B1, SH2B3/ATXN2, CSK, CYP17A1, FURIN, HFE, LSP1, MTHFR, SOX6) at array-wide significance (P < 2.4 × 10(-6)). We then replicated these associations in an independent set of 65 886 individuals of European ancestry. The findings from expression QTL (eQTL) analysis showed associations of SNPs in the MDM4 region with MDM4 expression. We did not find any evidence of association of the two novel SNPs in MDM4 and HRH1 with sequelae of high BP including coronary artery disease (CAD), left ventricular hypertrophy (LVH) or stroke. In summary, we identified two novel loci associated with BP and confirmed multiple previously reported associations. Our findings extend our understanding of genes involved in BP regulation, some of which may eventually provide new targets for therapeutic intervention.

FMR1, circadian genes and depression: suggestive associations or false discovery?

Background

There are several indications that malfunctions of the circadian clock contribute to depression. To search for particular circadian gene polymorphisms associated with depression, diverse polymorphisms were genotyped in two samples covering a range of depressed volunteers and participants with normal mood.

Methods

Depression mood self-ratings and DNA were collected independently from a sample of patients presenting to a sleep disorders center (1086 of European origin) and from a separate sample consisting of 399 participants claiming delayed sleep phase symptoms and 406 partly-matched controls. A custom Illumina Golden Gate array of 768 selected single nucleotide polymorphisms (SNPs) was assayed in both samples, supplemented by additional SNPlex and Taqman assays, including assay of 41 ancestry-associated markers (AIMs) to control stratification.

Results

In the Sleep Clinic sample, these assays yielded Bonferroni-significant association with depressed mood in three linked SNPs of the gene FMR1: rs25702 (nominal P=1.77E-05), rs25714 (P=1.83E-05), and rs28900 (P=5.24E-05). This FMR1 association was supported by 8 SNPs with nominal significance and a nominally-significant gene-wise set test. There was no association of depressed mood with FMR1 in the delayed sleep phase case–control sample or in downloaded GWAS data from the GenRED 2 sample contrasting an early-onset recurrent depression sample with controls. No replication was located in other GWAS studies of depression. Our data did weakly replicate a previously-reported association of depression with PPARGC1B rs7732671 (P=0.0235). Suggestive associations not meeting strict criteria for multiple testing and replication were found with GSK3B, NPAS2, RORA, PER3, CRY1, MTNR1A and NR1D1. Notably, 16 SNPs nominally associated with depressed mood (14 in GSK3B) were also nominally associated with delayed sleep phase syndrome (P=3E10^-6).

Conclusions

Considering the inconsistencies between samples and the likelihood that the significant three FMR1 SNPs might be linked to complex polymorphisms more functionally related to depression, large gene resequencing studies may be needed to clarify the import for depression of these circadian genes.

Common genetic variation near the connexin-43 gene is associated with resting heart rate in African Americans: a genome-wide association study of 13,372 participants

BACKGROUND

Genome-wide association studies have identified several genetic loci associated with variation in resting heart rate in European and Asian populations. No study has evaluated genetic variants associated with heart rate in African Americans.

OBJECTIVE

To identify novel genetic variants associated with resting heart rate in African Americans.

METHODS

Ten cohort studies participating in the Candidate-gene Association Resource and Continental Origins and Genetic Epidemiology Network consortia performed genome-wide genotyping of single nucleotide polymorphisms (SNPs) and imputed 2,954,965 SNPs using HapMap YRI and CEU panels in 13,372 participants of African ancestry. Each study measured the RR interval (ms) from 10-second resting 12-lead electrocardiograms and estimated RR-SNP associations using covariate-adjusted linear regression. Random-effects meta-analysis was used to combine cohort-specific measures of association and identify genome-wide significant loci (P≤2.5×10(-8)).

RESULTS

Fourteen SNPs on chromosome 6q22 exceeded the genome-wide significance threshold. The most significant association was for rs9320841 (+13 ms per minor allele; P = 4.98×10(-15)). This SNP was approximately 350 kb downstream of GJA1, a locus previously identified as harboring SNPs associated with heart rate in Europeans. Adjustment for rs9320841 also attenuated the association between the remaining 13 SNPs in this region and heart rate. In addition, SNPs in MYH6, which have been identified in European genome-wide association study, were associated with similar changes in the resting heart rate as this population of African Americans.

CONCLUSIONS

An intergenic region downstream of GJA1 (the gene encoding connexin 43, the major protein of the human myocardial gap junction) and an intragenic region within MYH6 are associated with variation in resting heart rate in African Americans as well as in populations of European and Asian origin.

Novel loci associated with PR interval in a genome-wide association study of 10 African American cohorts

BACKGROUND

The PR interval, as measured by the resting, standard 12-lead ECG, reflects the duration of atrial/atrioventricular nodal depolarization. Substantial evidence exists for a genetic contribution to PR, including genome-wide association studies that have identified common genetic variants at 9 loci influencing PR in populations of European and Asian descent. However, few studies have examined loci associated with PR in African Americans.

METHODS AND RESULTS

We present results from the largest genome-wide association study to date of PR in 13 415 adults of African descent from 10 cohorts. We tested for association between PR (ms) and ≈2.8 million genotyped and imputed single-nucleotide polymorphisms. Imputation was performed using HapMap 2 YRI and CEU panels. Study-specific results, adjusted for global ancestry and clinical correlates of PR, were meta-analyzed using the inverse variance method. Variation in genome-wide test statistic distributions was noted within studies (λ range: 0.9-1.1), although not after genomic control correction was applied to the overall meta-analysis (λ: 1.008). In addition to generalizing previously reported associations with MEIS1, SCN5A, ARHGAP24, CAV1, and TBX5 to African American populations at the genome-wide significance level (P<5.0 × 10(-8)), we also identified a novel locus: ITGA9, located in a region previously implicated in SCN5A expression. The 3p21 region harboring SCN5A also contained 2 additional independent secondary signals influencing PR (P<5.0 × 10(-8)).

CONCLUSIONS

This study demonstrates the ability to map novel loci in African Americans as well as the generalizability of loci associated with PR across populations of African, European, and Asian descent.

Impact of ancestry and common genetic variants on QT interval in African Americans

BACKGROUND

Ethnic differences in cardiac arrhythmia incidence have been reported, with a particularly high incidence of sudden cardiac death and low incidence of atrial fibrillation in individuals of African ancestry. We tested the hypotheses that African ancestry and common genetic variants are associated with prolonged duration of cardiac repolarization, a central pathophysiological determinant of arrhythmia, as measured by the electrocardiographic QT interval.

METHODS AND RESULTS

First, individual estimates of African and European ancestry were inferred from genome-wide single-nucleotide polymorphism (SNP) data in 7 population-based cohorts of African Americans (n=12,097) and regressed on measured QT interval from ECGs. Second, imputation was performed for 2.8 million SNPs, and a genome-wide association study of QT interval was performed in 10 cohorts (n=13,105). There was no evidence of association between genetic ancestry and QT interval (P=0.94). Genome-wide significant associations (P<2.5 × 10(-8)) were identified with SNPs at 2 loci, upstream of the genes NOS1AP (rs12143842, P=2 × 10(-15)) and ATP1B1 (rs1320976, P=2 × 10(-10)). The most significant SNP in NOS1AP was the same as the strongest SNP previously associated with QT interval in individuals of European ancestry. Low probability values (P<10(-5)) were observed for SNPs at several other loci previously identified in genome-wide association studies in individuals of European ancestry, including KCNQ1, KCNH2, LITAF, and PLN.

CONCLUSIONS

We observed no difference in duration of cardiac repolarization with global genetic indices of African American ancestry. In addition, our genome-wide association study extends the association of polymorphisms at several loci associated with repolarization in individuals of European ancestry to include individuals of African ancestry.

Neuroanatomical assessment of biological maturity

Structural MRI allows unparalleled in vivo study of the anatomy of the developing human brain. For more than two decades, MRI research has revealed many new aspects of this multifaceted maturation process, significantly augmenting scientific knowledge gathered from postmortem studies. Postnatal brain development is notably protracted and involves considerable changes in cerebral cortical, subcortical, and cerebellar structures, as well as significant architectural changes in white matter fiber tracts (see [12]). Although much work has described isolated features of neuroanatomical development, it remains a critical challenge to characterize the multidimensional nature of brain anatomy, capturing different phases of development among individuals. Capitalizing on key advances in multisite, multimodal MRI, and using cross-validated nonlinear modeling, we demonstrate that developmental brain phase can be assessed with much greater precision than has been possible using other biological measures, accounting for more than 92% of the variance in age. Further, our composite metric of morphology, diffusivity, and signal intensity shows that the average difference in phase among children of the same age is only about 1 year, revealing for the first time a latent phenotype in the human brain for which maturation timing is tightly controlled.

Genetic variants and blood pressure in a population-based cohort: the Cardiovascular Risk in Young Finns study

Clinical relevance of a genetic predisposition to elevated blood pressure was quantified during the transition from childhood to adulthood in a population-based Finnish cohort (N=2357). Blood pressure was measured at baseline in 1980 (age 3-18 years) and in follow-ups in 1983, 1986, 2001, and 2007. Thirteen single nucleotide polymorphisms associated with blood pressure were genotyped, and 3 genetic risk scores associated with systolic and diastolic blood pressures and their combination were derived for all of the participants. Effects of the genetic risk score were 0.47 mm Hg for systolic and 0.53 mm Hg for diastolic blood pressures (both P<0.01). The combination genetic risk score was associated with diastolic blood pressure from age 9 years onward (β=0.68 mm Hg; P=0.015). Replications in 1194 participants of the Bogalusa Heart Study showed essentially similar results. The participants in the highest quintile of the combination genetic risk score had a 1.82-fold risk of hypertension in adulthood (P<0.0001) compared with the lowest quintile, independent of a family history of premature hypertension. These findings show that genetic variants are associated with preclinical blood pressure traits in childhood; individuals with several susceptibility alleles have, on average, a 0.5-mm Hg higher blood pressure, and this trajectory continues from childhood to adulthood.

Analysis of 94 candidate genes and 12 endophenotypes for schizophrenia from the Consortium on the Genetics of Schizophrenia

OBJECTIVE

The authors used a custom array of 1,536 single-nucleotide polymorphisms (SNPs) to interrogate 94 functionally relevant candidate genes for schizophrenia and identify associations with 12 heritable neurophysiological and neurocognitive endophenotypes in data collected by the Consortium on the Genetics of Schizophrenia.

METHOD

Variance-component association analyses of 534 genotyped subjects from 130 families were conducted by using Merlin software. A novel bootstrap total significance test was also developed to overcome the limitations of existing genomic multiple testing methods and robustly demonstrate significant associations in the context of complex family data and possible population stratification effects.

RESULTS

Associations with endophenotypes were observed for 46 genes of potential functional significance, with three SNPs at p<10(-4), 27 SNPs at p<10(-3), and 147 SNPs at p<0.01. The bootstrap analyses confirmed that the 47 SNP-endophenotype combinations with the strongest evidence of association significantly exceeded that expected by chance alone, with 93% of these findings expected to be true. Many of the genes interact on a molecular level, and eight genes (e.g., NRG1 and ERBB4) displayed evidence for pleiotropy, revealing associations with four or more endophenotypes. The results collectively support a strong role for genes related to glutamate signaling in mediating schizophrenia susceptibility.

CONCLUSIONS

This study supports use of relevant endophenotypes and the bootstrap total significance test for identifying genetic variation underlying the etiology of schizophrenia. In addition, the observation of extensive pleiotropy for some genes and singular associations for others suggests alternative, independent pathways mediating pathogenesis in the "group of schizophrenias."

Gene expression profiling of human whole blood samples with the Illumina WG-DASL assay

Background

Microarray-based gene expression analysis of peripheral whole blood is a common strategy in the development of clinically relevant biomarker panels for a variety of human diseases. However, the results of such an analysis are often plagued by decreased sensitivity and reliability due to the effects of relatively high levels of globin mRNA in whole blood. Globin reduction assays have been shown to overcome such effects, but they require large amounts of total RNA and may induce distinct gene expression profiles. The Illumina whole genome DASL assay can detect gene expression levels using partially degraded RNA samples and has the potential to detect rare transcripts present in highly heterogeneous whole blood samples without the need for globin reduction. We assessed the utility of the whole genome DASL assay in an analysis of peripheral whole blood gene expression profiles.

Results

We find that gene expression detection is significantly increased with the use of whole genome DASL compared to the standard IVT-based direct hybridization. Additionally, globin-probe negative whole genome DASL did not exhibit significant improvements over globin-probe positive whole genome DASL. Globin reduction further increases the detection sensitivity and reliability of both whole genome DASL and IVT-based direct hybridization with little effect on raw intensity correlations. Raw intensity correlations between total RNA and globin reduced RNA were 0.955 for IVT-based direct hybridization and 0.979 for whole genome DASL.

Conclusions

Overall, the detection sensitivity of the whole genome DASL assay is higher than the IVT-based direct hybridization assay, with or without globin reduction, and should be considered in conjunction with globin reduction methods for future blood-based gene expression studies.

Polymorphisms in melatonin synthesis pathways: possible influences on depression

Background

It has been reported that rs4446909, a single nucleotide polymorphism (SNP) in the promoter of acetylserotonin methyltransferase (ASMT), influences the expression of the ASMT enzyme. The common G allele is associated with lower ASMT activity, and therefore, diminishes conversion of N-acetylserotonin to melatonin. The G allele was associated with recurrent depressive disorder in a Polish group. ASMT might also affect bipolar relapse, given evidence that N-acetylserotonin might stimulate TRKB receptors, and TRKB may influence mood relapse in bipolar disorder. Additionally, arylalkylamine N-acetyltransferase (AANAT) polymorphisms have been reported associated with depression, perhaps through their influence upon N-acetylserotonin or melatonin synthesis.

Results

To replicate and further explore these ideas, rs4446909 was genotyped in four research groups, as part of a panel of 610 SNPs surveyed by an Illumina Golden Gate assay. In 768 cases with delayed sleep phase disorder or matched controls, rs4446909 was indeed associated with the depressive symptoms on a self-report scale (P = 0.01, R² = 0.007). However, there was no significant association of rs4446909 with self-reported depression in a sleep clinic patient group or with two groups of elderly men and women from multicenter studies, nor was the response to lithium treatment associated with rs4446909 in bipolar patients. No associations of two AANAT SNPs with depression were found.

Conclusions

The evidence did not support a strong influence of rs4446909 upon mood, but the partial replication may be consistent with a modest effect. It is possible that larger or younger subject groups with improved phenotype ascertainment might demonstrate more persuasive replication.

Thirty new loci for age at menarche identified by a meta-analysis of genome-wide association studies

To identify loci for age at menarche, we performed a meta-analysis of 32 genome-wide association studies in 87,802 women of European descent, with replication in up to 14,731 women. In addition to the known loci at LIN28B (P = 5.4 × 10⁻⁶⁰) and 9q31.2 (P = 2.2 × 10⁻³³), we identified 30 new menarche loci (all P < 5 × 10⁻⁸) and found suggestive evidence for a further 10 loci (P < 1.9 × 10⁻⁶). The new loci included four previously associated with body mass index (in or near FTO, SEC16B, TRA2B and TMEM18), three in or near other genes implicated in energy homeostasis (BSX, CRTC1 and MCHR2) and three in or near genes implicated in hormonal regulation (INHBA, PCSK2 and RXRG). Ingenuity and gene-set enrichment pathway analyses identified coenzyme A and fatty acid biosynthesis as biological processes related to menarche timing.

Longitudinal genome-wide association of cardiovascular disease risk factors in the Bogalusa heart study

Cardiovascular disease (CVD) is the leading cause of death worldwide. Recent genome-wide association (GWA) studies have pinpointed many loci associated with CVD risk factors in adults. It is unclear, however, if these loci predict trait levels at all ages, if they are associated with how a trait develops over time, or if they could be used to screen individuals who are pre-symptomatic to provide the opportunity for preventive measures before disease onset. We completed a genome-wide association study on participants in the longitudinal Bogalusa Heart Study (BHS) and have characterized the association between genetic factors and the development of CVD risk factors from childhood to adulthood. We report 7 genome-wide significant associations involving CVD risk factors, two of which have been previously reported. Top regions were tested for replication in the Young Finns Study (YF) and two associations strongly replicated: rs247616 in CETP with HDL levels (combined P = 9.7×10⁻²⁴), and rs445925 at APOE with LDL levels (combined P = 8.7×10⁻¹⁹). We show that SNPs previously identified in adult cross-sectional studies tend to show age-independent effects in the BHS with effect sizes consistent with previous reports. Previously identified variants were associated with adult trait levels above and beyond those seen in childhood; however, variants with time-dependent effects were also promising predictors. This is the first GWA study to evaluate the role of common genetic variants in the development of CVD risk factors in children as they advance through adulthood and highlights the utility of using longitudinal studies to identify genetic predictors of adult traits in children.

We have studied the association between genetic factors on a whole genome level and cardiovascular disease (CVD) risk factors in a population of individuals studied from childhood through adulthood. The longitudinal study design has enabled the investigation of genetic variation influencing trait values over time. We have identified DNA variants that are associated with CVD trait values consistently over time, and a second set of variants that are associated with CVD trait values in a time-dependent manner. We also show that variants previously identified in adult populations have consistent effects within our population and that these effects are usually similar across childhood through adulthood. The discovery of time-dependent variants that influence CVD trait values over time can potentially be used to screen young individuals who are pre-symptomatic and provide the opportunity for preventive measures decades before disease onset.

Genome-wide association identifies OBFC1 as a locus involved in human leukocyte telomere biology

Telomeres are engaged in a host of cellular functions, and their length is regulated by multiple genes. Telomere shortening, in the course of somatic cell replication, ultimately leads to replicative senescence. In humans, rare mutations in genes that regulate telomere length have been identified in monogenic diseases such as dyskeratosis congenita and idiopathic pulmonary fibrosis, which are associated with shortened leukocyte telomere length (LTL) and increased risk for aplastic anemia. Shortened LTL is observed in a host of aging-related complex genetic diseases and is associated with diminished survival in the elderly. We report results of a genome-wide association study of LTL in a consortium of four observational studies (n = 3,417 participants with LTL and genome-wide genotyping). SNPs in the regions of the oligonucleotide/oligosaccharide-binding folds containing one gene (OBFC1; rs4387287; P = 3.9 x 10(-9)) and chemokine (C-X-C motif) receptor 4 gene (CXCR4; rs4452212; P = 2.9 x 10(-8)) were associated with LTL at a genome-wide significance level (P < 5 x 10(-8)). We attempted replication of the top SNPs at these loci through de novo genotyping of 1,893 additional individuals and in silico lookup in another observational study (n = 2,876), and we confirmed the association findings for OBFC1 but not CXCR4. In addition, we confirmed the telomerase RNA component (TERC) as a gene associated with LTL (P = 1.1 x 10(-5)). The identification of OBFC1 through genome-wide association as a locus for interindividual variation in LTL in the general population advances the understanding of telomere biology in humans and may provide insights into aging-related disorders linked to altered LTL dynamics.

Accurate detection and genotyping of SNPs utilizing population sequencing data

Next-generation sequencing technologies have made it possible to sequence targeted regions of the human genome in hundreds of individuals. Deep sequencing represents a powerful approach for the discovery of the complete spectrum of DNA sequence variants in functionally important genomic intervals. Current methods for single nucleotide polymorphism (SNP) detection are designed to detect SNPs from single individual sequence data sets. Here, we describe a novel method SNIP-Seq (single nucleotide polymorphism identification from population sequence data) that leverages sequence data from a population of individuals to detect SNPs and assign genotypes to individuals. To evaluate our method, we utilized sequence data from a 200-kilobase (kb) region on chromosome 9p21 of the human genome. This region was sequenced in 48 individuals (five sequenced in duplicate) using the Illumina GA platform. Using this data set, we demonstrate that our method is highly accurate for detecting variants and can filter out false SNPs that are attributable to sequencing errors. The concordance of sequencing-based genotype assignments between duplicate samples was 98.8%. The 200-kb region was independently sequenced to a high depth of coverage using two sequence pools containing the 48 individuals. Many of the novel SNPs identified by SNIP-Seq from the individual sequencing were validated by the pooled sequencing data and were subsequently confirmed by Sanger sequencing. We estimate that SNIP-Seq achieves a low false-positive rate of approximately 2%, improving upon the higher false-positive rate for existing methods that do not utilize population sequence data. Collectively, these results suggest that analysis of population sequencing data is a powerful approach for the accurate detection of SNPs and the assignment of genotypes to individual samples.

Common vs. rare allele hypotheses for complex diseases

There has been growing debate over the nature of the genetic contribution to individual susceptibility to common complex diseases such as diabetes, osteoporosis, and cancer. The 'Common Disease, Common Variant (CDCV)' hypothesis argues that genetic variations with appreciable frequency in the population at large, but relatively low 'penetrance' (or the probability that a carrier of the relevant variants will express the disease), are the major contributors to genetic susceptibility to common diseases. The 'Common Disease, Rare Variant (CDRV)' hypothesis, on the contrary, argues that multiple rare DNA sequence variations, each with relatively high penetrance, are the major contributors to genetic susceptibility to common diseases. Both hypotheses have their place in current research efforts.

Human genetic variation and its contribution to complex traits

The last few years have seen extensive efforts to catalogue human genetic variation and correlate it with phenotypic differences. Most common SNPs have now been assessed in genome-wide studies for statistical associations with many complex traits, including many important common diseases. Although these studies have provided new biological insights, only a limited amount of the heritable component of any complex trait has been identified and it remains a challenge to elucidate the functional link between associated variants and phenotypic traits. Technological advances, such as the ability to detect rare and structural variants, and a clear understanding of the challenges in linking different types of variation with phenotype, will be essential for future progress.