Post-genome-wide association study challenges for lipid traits: describing age as a modifier of gene-lipid associations in the Population Architecture using Genomics and Epidemiology (PAGE) study

Numerous common genetic variants that influence plasma high-density lipoprotein cholesterol, low-density lipoprotein cholesterol (LDL-C), and triglyceride distributions have been identified via genome-wide association studies (GWAS). However, whether or not these associations are age-dependent has largely been overlooked. We conducted an association study and meta-analysis in more than 22,000 European Americans between 49 previously identified GWAS variants and the three lipid traits, stratified by age (males: <50 or ≥50 years of age; females: pre- or postmenopausal). For each variant, a test of heterogeneity was performed between the two age strata and significant Phet values were used as evidence of age-specific genetic effects. We identified seven associations in females and eight in males that displayed suggestive heterogeneity by age (Phet < 0.05). The association between rs174547 (FADS1) and LDL-C in males displayed the most evidence for heterogeneity between age groups (Phet = 1.74E-03, I(2) = 89.8), with a significant association in older males (P = 1.39E-06) but not younger males (P = 0.99). However, none of the suggestive modifying effects survived adjustment for multiple testing, highlighting the challenges of identifying modifiers of modest SNP-trait associations despite large sample sizes.

The influence of obesity-related single nucleotide polymorphisms on BMI across the life course: the PAGE study

Evidence is limited as to whether heritable risk of obesity varies throughout adulthood. Among >34,000 European Americans, aged 18-100 years, from multiple U.S. studies in the Population Architecture using Genomics and Epidemiology (PAGE) Consortium, we examined evidence for heterogeneity in the associations of five established obesity risk variants (near FTO, GNPDA2, MTCH2, TMEM18, and NEGR1) with BMI across four distinct epochs of adulthood: 1) young adulthood (ages 18-25 years), adulthood (ages 26-49 years), middle-age adulthood (ages 50-69 years), and older adulthood (ages ≥70 years); or 2) by menopausal status in women and stratification by age 50 years in men. Summary-effect estimates from each meta-analysis were compared for heterogeneity across the life epochs. We found heterogeneity in the association of the FTO (rs8050136) variant with BMI across the four adulthood epochs (P = 0.0006), with larger effects in young adults relative to older adults (β [SE] = 1.17 [0.45] vs. 0.09 [0.09] kg/m², respectively, per A allele) and smaller intermediate effects. We found no evidence for heterogeneity in the association of GNPDA2, MTCH2, TMEM18, and NEGR1 with BMI across adulthood. Genetic predisposition to obesity may have greater effects on body weight in young compared with older adulthood for FTO, suggesting changes by age, generation, or secular trends. Future research should compare and contrast our findings with results using longitudinal data.

Investigation of gene-by-sex interactions for lipid traits in diverse populations from the population architecture using genomics and epidemiology study

BACKGROUND

High-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglyceride (TG) levels are influenced by both genes and the environment. Genome-wide association studies (GWAS) have identified ~100 common genetic variants associated with HDL-C, LDL-C, and/or TG levels, mostly in populations of European descent, but little is known about the modifiers of these associations. Here, we investigated whether GWAS-identified SNPs for lipid traits exhibited heterogeneity by sex in the Population Architecture using Genomics and Epidemiology (PAGE) study.

RESULTS

A sex-stratified meta-analysis was performed for 49 GWAS-identified SNPs for fasting HDL-C, LDL-C, and ln(TG) levels among adults self-identified as European American (25,013). Heterogeneity by sex was established when phet < 0.001. There was evidence for heterogeneity by sex for two SNPs for ln(TG) in the APOA1/C3/A4/A5/BUD13 gene cluster: rs28927680 (p(het) = 7.4 x 10(-7)) and rs3135506 (p(het) = 4.3 x 10(-4)one SNP in PLTP for HDL levels (rs7679; p(het) = 9.9 x 10(-4)), and one in HMGCR for LDL levels (rs12654264; p(het) = 3.1 x 10(-5)). We replicated heterogeneity by sex in five of seventeen loci previously reported by genome-wide studies (binomial p = 0.0009). We also present results for other racial/ethnic groups in the supplementary materials, to provide a resource for future meta-analyses.

CONCLUSIONS

We provide further evidence for sex-specific effects of SNPs in the APOA1/C3/A4/A5/BUD13 gene cluster, PLTP, and HMGCR on fasting triglyceride levels in European Americans from the PAGE study. Our findings emphasize the need for considering context-specific effects when interpreting genetic associations emerging from GWAS, and also highlight the difficulties in replicating interaction effects across studies and across racial/ethnic groups.

Genetic risk factors for BMI and obesity in an ethnically diverse population: results from the population architecture using genomics and epidemiology (PAGE) study

OBJECTIVE

Several genome-wide association studies (GWAS) have demonstrated that common genetic variants contribute to obesity. However, studies of this complex trait have focused on ancestrally European populations, despite the high prevalence of obesity in some minority groups.

DESIGN AND METHODS

As part of the "Population Architecture using Genomics and Epidemiology (PAGE)" Consortium, we investigated the association between 13 GWAS-identified single-nucleotide polymorphisms (SNPs) and BMI and obesity in 69,775 subjects, including 6,149 American Indians, 15,415 African-Americans, 2,438 East Asians, 7,346 Hispanics, 604 Pacific Islanders, and 37,823 European Americans. For the BMI-increasing allele of each SNP, we calculated β coefficients using linear regression (for BMI) and risk estimates using logistic regression (for obesity defined as BMI ≥ 30) followed by fixed-effects meta-analysis to combine results across PAGE sites. Analyses stratified by racial/ethnic group assumed an additive genetic model and were adjusted for age, sex, and current smoking. We defined "replicating SNPs" (in European Americans) and "generalizing SNPs" (in other racial/ethnic groups) as those associated with an allele frequency-specific increase in BMI.

RESULTS

By this definition, we replicated 9/13 SNP associations (5 out of 8 loci) in European Americans. We also generalized 8/13 SNP associations (5/8 loci) in East Asians, 7/13 (5/8 loci) in African Americans, 6/13 (4/8 loci) in Hispanics, 5/8 in Pacific Islanders (5/8 loci), and 5/9 (4/8 loci) in American Indians.

CONCLUSION

Linkage disequilibrium patterns suggest that tagSNPs selected for European Americans may not adequately tag causal variants in other ancestry groups. Accordingly, fine-mapping in large samples is needed to comprehensively explore these loci in diverse populations.

Trans-ethnic fine-mapping of lipid loci identifies population-specific signals and allelic heterogeneity that increases the trait variance explained

Genome-wide association studies (GWAS) have identified ∼100 loci associated with blood lipid levels, but much of the trait heritability remains unexplained, and at most loci the identities of the trait-influencing variants remain unknown. We conducted a trans-ethnic fine-mapping study at 18, 22, and 18 GWAS loci on the Metabochip for their association with triglycerides (TG), high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C), respectively, in individuals of African American (n = 6,832), East Asian (n = 9,449), and European (n = 10,829) ancestry. We aimed to identify the variants with strongest association at each locus, identify additional and population-specific signals, refine association signals, and assess the relative significance of previously described functional variants. Among the 58 loci, 33 exhibited evidence of association at P<1×10⁻⁴ in at least one ancestry group. Sequential conditional analyses revealed that ten, nine, and four loci in African Americans, Europeans, and East Asians, respectively, exhibited two or more signals. At these loci, accounting for all signals led to a 1.3- to 1.8-fold increase in the explained phenotypic variance compared to the strongest signals. Distinct signals across ancestry groups were identified at PCSK9 and APOA5. Trans-ethnic analyses narrowed the signals to smaller sets of variants at GCKR, PPP1R3B, ABO, LCAT, and ABCA1. Of 27 variants reported previously to have functional effects, 74% exhibited the strongest association at the respective signal. In conclusion, trans-ethnic high-density genotyping and analysis confirm the presence of allelic heterogeneity, allow the identification of population-specific variants, and limit the number of candidate SNPs for functional studies.

Lipid traits are heritable, but many of the DNA variants that influence lipid levels remain unknown. In a genomic region, more than one variant may affect gene expression or function, and the frequencies of these variants can differ across populations. Genotyping densely spaced variants in individuals with different ancestries may increase the chance of identifying variants that affect gene expression or function. We analyzed high-density genotyped variants for association with TG, HDL-C, and LDL-C in African Americans, East Asians, and Europeans. At several genomic regions, we provide evidence that two or more variants can influence lipid traits; across loci, these additional signals increase the proportion of trait variation that can be explained by genes. At some association signals shared across populations, combining data from individuals of different ancestries narrowed the set of likely functional variants. At PCSK9 and APOA5, the data suggest that different variants influence trait levels in different populations. Variants previously reported to alter gene expression or function frequently exhibited the strongest association at those signals. The multiple signals and population-specific characteristics of the loci described here may be shared by genetic loci for other complex traits.

Replication of genetic loci for ages at menarche and menopause in the multi-ethnic Population Architecture using Genomics and Epidemiology (PAGE) study

STUDY QUESTION

Do genetic associations identified in genome-wide association studies (GWAS) of age at menarche (AM) and age at natural menopause (ANM) replicate in women of diverse race/ancestry from the Population Architecture using Genomics and Epidemiology (PAGE) Study?

SUMMARY ANSWER

We replicated GWAS reproductive trait single nucleotide polymorphisms (SNPs) in our European descent population and found that many SNPs were also associated with AM and ANM in populations of diverse ancestry.

WHAT IS KNOWN ALREADY

Menarche and menopause mark the reproductive lifespan in women and are important risk factors for chronic diseases including obesity, cardiovascular disease and cancer. Both events are believed to be influenced by environmental and genetic factors, and vary in populations differing by genetic ancestry and geography. Most genetic variants associated with these traits have been identified in GWAS of European-descent populations.

STUDY DESIGN, SIZE, DURATION

A total of 42 251 women of diverse ancestry from PAGE were included in cross-sectional analyses of AM and ANM.

MATERIALS, SETTING, METHODS

SNPs previously associated with ANM (n = 5 SNPs) and AM (n = 3 SNPs) in GWAS were genotyped in American Indians, African Americans, Asians, European Americans, Hispanics and Native Hawaiians. To test SNP associations with ANM or AM, we used linear regression models stratified by race/ethnicity and PAGE sub-study. Results were then combined in race-specific fixed effect meta-analyses for each outcome. For replication and generalization analyses, significance was defined at P < 0.01 for ANM analyses and P < 0.017 for AM analyses.

MAIN RESULTS AND THE ROLE OF CHANCE

We replicated findings for AM SNPs in the LIN28B locus and an intergenic region on 9q31 in European Americans. The LIN28B SNPs (rs314277 and rs314280) were also significantly associated with AM in Asians, but not in other race/ethnicity groups. Linkage disequilibrium (LD) patterns at this locus varied widely among the ancestral groups. With the exception of an intergenic SNP at 13q34, all ANM SNPs replicated in European Americans. Three were significantly associated with ANM in other race/ethnicity populations: rs2153157 (6p24.2/SYCP2L), rs365132 (5q35/UIMC1) and rs16991615 (20p12.3/MCM8). While rs1172822 (19q13/BRSK1) was not significant in the populations of non-European descent, effect sizes showed similar trends.

LIMITATIONS, REASONS FOR CAUTION

Lack of association for the GWAS SNPs in the non-European American groups may be due to differences in locus LD patterns between these groups and the European-descent populations included in the GWAS discovery studies; and in some cases, lower power may also contribute to non-significant findings.

WIDER IMPLICATIONS OF THE FINDINGS

The discovery of genetic variants associated with the reproductive traits provides an important opportunity to elucidate the biological mechanisms involved with normal variation and disorders of menarche and menopause. In this study we replicated most, but not all reported SNPs in European descent populations and examined the epidemiologic architecture of these early reported variants, describing their generalizability and effect size across differing ancestral populations. Such data will be increasingly important for prioritizing GWAS SNPs for follow-up in fine-mapping and resequencing studies, as well as in translational research.

Pain is associated with short leukocyte telomere length in women with fibromyalgia

Telomere length, considered a measure of biological aging, is linked to morbidity and mortality. Psychosocial factors associated with shortened telomeres are also common in chronic pain; yet, little is known about telomere length in pain populations. Leukocyte telomere length was evaluated in 66 women with fibromyalgia and 22 healthy female controls. Participants completed questionnaires and a subgroup of fibromyalgia patients underwent quantitative sensory testing (QST; n = 12) and neuroimaging (n = 12). Telomere length was measured using the quantitative polymerase chain reaction method. Although patients had shorter telomere length than controls, the difference was not statistically significant. However, higher levels of pain within fibromyalgia were associated with shorter telomere length (P = .039). When pain and depression were combined, patients categorized as high-pain/high-depression had an age-adjusted telomere length 265 base pairs shorter than those with low-pain/low-depression (P = .043), a difference consistent with approximately 6 years of chronological aging. In the subset tested, telomere length was also related to pain threshold and pain sensitivity, as well as gray matter volume, such that patients with shorter telomeres were more sensitive to evoked pain and had less gray matter in brain regions associated with pain processing (eg, primary somatosensory cortex). These preliminary data support a relationship between pain and telomere length.

PERSPECTIVE

Our findings support a link between premature cellular aging and chronic pain. These preliminary data imply that chronic pain is a more serious condition than has typically been recognized in terms of bodily aging.

Genetic variation and reproductive timing: African American women from the Population Architecture using Genomics and Epidemiology (PAGE) Study

Age at menarche (AM) and age at natural menopause (ANM) define the boundaries of the reproductive lifespan in women. Their timing is associated with various diseases, including cancer and cardiovascular disease. Genome-wide association studies have identified several genetic variants associated with either AM or ANM in populations of largely European or Asian descent women. The extent to which these associations generalize to diverse populations remains unknown. Therefore, we sought to replicate previously reported AM and ANM findings and to identify novel AM and ANM variants using the Metabochip (n = 161,098 SNPs) in 4,159 and 1,860 African American women, respectively, in the Women's Health Initiative (WHI) and Atherosclerosis Risk in Communities (ARIC) studies, as part of the Population Architecture using Genomics and Epidemiology (PAGE) Study. We replicated or generalized one previously identified variant for AM, rs1361108/CENPW, and two variants for ANM, rs897798/BRSK1 and rs769450/APOE, to our African American cohort. Overall, generalization of the majority of previously-identified variants for AM and ANM, including LIN28B and MCM8, was not observed in this African American sample. We identified three novel loci associated with ANM that reached significance after multiple testing correction (LDLR rs189596789, p = 5×10⁻⁰⁸; KCNQ1 rs79972789, p = 1.9×10⁻⁰⁷; COL4A3BP rs181686584, p = 2.9×10⁻⁰⁷). Our most significant AM association was upstream of RSF1, a gene implicated in ovarian and breast cancers (rs11604207, p = 1.6×10⁻⁰⁶). While most associations were identified in either AM or ANM, we did identify genes suggestively associated with both: PHACTR1 and ARHGAP42. The lack of generalization coupled with the potentially novel associations identified here emphasize the need for additional genetic discovery efforts for AM and ANM in diverse populations.

Genetic association of GABA-A receptor alpha-2 and mu opioid receptor with cocaine cue-reactivity: evidence for inhibitory synaptic neurotransmission involvement in cocaine dependence

BACKGROUND

This pilot feasibility study examined the role of genetics in laboratory-induced cocaine craving.

METHODS

Thirty-four African American, cocaine-depend- ent male subjects underwent a baseline assessment, cue-exposure session, and genetic analysis. Subjects were classified as either cue-reactive or nonreactive.

RESULTS

Among single nucleotide polymorphism markers in 13 candidate genes examined for association with cocaine cue-reactivity, two were statistically significant: GABRA2 (coding for GABA-A receptor alpha-2 subunit; rs11503014, nominal p= .001) and OPRM1 (coding for mu opioid receptor; rs2236256, nominal p= .03).

CONCLUSIONS

These pilot results suggest that cocaine craving shows variability among cocaine-dependent subjects, and that GABRA2 and OPRM1 polymorphisms have differential influences on cocaine cue-reactivity, warranting studies in future research.

Fine-mapping and initial characterization of QT interval loci in African Americans

The QT interval (QT) is heritable and its prolongation is a risk factor for ventricular tachyarrhythmias and sudden death. Most genetic studies of QT have examined European ancestral populations; however, the increased genetic diversity in African Americans provides opportunities to narrow association signals and identify population-specific variants. We therefore evaluated 6,670 SNPs spanning eleven previously identified QT loci in 8,644 African American participants from two Population Architecture using Genomics and Epidemiology (PAGE) studies: the Atherosclerosis Risk in Communities study and Women's Health Initiative Clinical Trial. Of the fifteen known independent QT variants at the eleven previously identified loci, six were significantly associated with QT in African American populations (P≤1.20×10(-4)): ATP1B1, PLN1, KCNQ1, NDRG4, and two NOS1AP independent signals. We also identified three population-specific signals significantly associated with QT in African Americans (P≤1.37×10(-5)): one at NOS1AP and two at ATP1B1. Linkage disequilibrium (LD) patterns in African Americans assisted in narrowing the region likely to contain the functional variants for several loci. For example, African American LD patterns showed that 0 SNPs were in LD with NOS1AP signal rs12143842, compared with European LD patterns that indicated 87 SNPs, which spanned 114.2 Kb, were in LD with rs12143842. Finally, bioinformatic-based characterization of the nine African American signals pointed to functional candidates located exclusively within non-coding regions, including predicted binding sites for transcription factors such as TBX5, which has been implicated in cardiac structure and conductance. In this detailed evaluation of QT loci, we identified several African Americans SNPs that better define the association with QT and successfully narrowed intervals surrounding established loci. These results demonstrate that the same loci influence variation in QT across multiple populations, that novel signals exist in African Americans, and that the SNPs identified as strong candidates for functional evaluation implicate gene regulatory dysfunction in QT prolongation.

Consistent directions of effect for established type 2 diabetes risk variants across populations: the population architecture using Genomics and Epidemiology (PAGE) Consortium

Common genetic risk variants for type 2 diabetes (T2D) have primarily been identified in populations of European and Asian ancestry. We tested whether the direction of association with 20 T2D risk variants generalizes across six major racial/ethnic groups in the U.S. as part of the Population Architecture using Genomics and Epidemiology Consortium (16,235 diabetes case and 46,122 control subjects of European American, African American, Hispanic, East Asian, American Indian, and Native Hawaiian ancestry). The percentage of positive (odds ratio [OR] >1 for putative risk allele) associations ranged from 69% in American Indians to 100% in European Americans. Of the nine variants where we observed significant heterogeneity of effect by racial/ethnic group (P(heterogeneity) < 0.05), eight were positively associated with risk (OR >1) in at least five groups. The marked directional consistency of association observed for most genetic variants across populations implies a shared functional common variant in each region. Fine-mapping of all loci will be required to reveal markers of risk that are important within and across populations.

Evaluation of the metabochip genotyping array in African Americans and implications for fine mapping of GWAS-identified loci: the PAGE study

The Metabochip is a custom genotyping array designed for replication and fine mapping of metabolic, cardiovascular, and anthropometric trait loci and includes low frequency variation content identified from the 1000 Genomes Project. It has 196,725 SNPs concentrated in 257 genomic regions. We evaluated the Metabochip in 5,863 African Americans; 89% of all SNPs passed rigorous quality control with a call rate of 99.9%. Two examples illustrate the value of fine mapping with the Metabochip in African-ancestry populations. At CELSR2/PSRC1/SORT1, we found the strongest associated SNP for LDL-C to be rs12740374 (p = 3.5 × 10(-11)), a SNP indistinguishable from multiple SNPs in European ancestry samples due to high correlation. Its distinct signal supports functional studies elsewhere suggesting a causal role in LDL-C. At CETP we found rs17231520, with risk allele frequency 0.07 in African Americans, to be associated with HDL-C (p = 7.2 × 10(-36)). This variant is very rare in Europeans and not tagged in common GWAS arrays, but was identified as associated with HDL-C in African Americans in a single-gene study. Our results, one narrowing the risk interval and the other revealing an associated variant not found in Europeans, demonstrate the advantages of high-density genotyping of common and rare variation for fine mapping of trait loci in African American samples.

Abstract 1662: Pleiotropic effects on lung cancer of genetic susceptibility variants identified for other malignancies: The Population Architecture using Genomics and Epidemiology Study

Genome-wide association studies of cancer have succeeded in identifying over 100 susceptibility loci, including some, such as those found at chromosome 8q24 and CPTM1L-TERT that influence the risk of several malignancies. As part of the Population Architecture using Genomics and Epidemiology (PAGE) Study, we evaluated whether established risk variants for 18 different cancer sites, excluding the lung, impact the risk of lung cancer. A total of 152 risk variants for these 18 malignancies were selected from genome-wide association studies of cancer published through May, 2009. These SNPs were genotyped in lung cancer case-control studies of the Atherosclerosis Risk in Communities (ARIC; cases/controls=297/3088), Epidemiologic Architecture for Genes Linked to Environment (EAGLE; cases/controls=432/1567), Multiethnic Cohort (MEC; cases/controls=461/9029), and the Women's Health Initiative (WHI; cases/controls=1600/5848). The total study population was comprised of 2,790 lung cancer cases and 19,532 controls of African, American Indian, Asian, European, Latino, and Pacific Islander ancestry. For each study site, unconditional logistic regression was performed to evaluate the association between SNPs and lung cancer risk, adjusting for age, sex, race/ethnicity, and smoking status. A fixed-effect meta-analysis was performed and heterogeneity across study sites was tested. Seventeen of the 152 SNPs were associated with lung cancer risk (nominal P&lt;0.05): five breast cancer SNPs, one glioma SNP, three acute lymphocytic leukemia SNPs, one melanoma SNP, one non-Hodgkin lymphoma SNP, and six prostate cancer SNPs. By chance, only eight associations would be expected for the 152 tests (152 x 0.05=8) performed, suggesting robust significant associations among the 17 associated SNPs. The top associations (P&lt;0.01) were observed for: breast cancer SNP (rs3803662 at TOX3: OR=1.13; P=1.1x10−3), prostate cancer SNP (rs7837688 at 8q24: OR=1.19; P=5.2x10−3), and melanoma SNP (rs910873 at PIGU: OR=0.79; P=5.9x10−3). No evidence of heterogeneous effects across study sites was observed. Other analyses will test for independent genetic effects among the associated loci and examine stratified effects by sex, cell-type, and smoking status. In conclusion, findings from our meta-analysis suggest that lung cancer may share genetic risk factors and common etiologic pathways with several cancer sites.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 1662. doi:1538-7445.AM2012-1662

Genotype imputation of Metabochip SNPs using a study-specific reference panel of ~4,000 haplotypes in African Americans from the Women's Health Initiative

Genetic imputation has become standard practice in modern genetic studies. However, several important issues have not been adequately addressed including the utility of study-specific reference, performance in admixed populations, and quality for less common (minor allele frequency [MAF] 0.005-0.05) and rare (MAF < 0.005) variants. These issues only recently became addressable with genome-wide association studies (GWAS) follow-up studies using dense genotyping or sequencing in large samples of non-European individuals. In this work, we constructed a study-specific reference panel of 3,924 haplotypes using African Americans in the Women's Health Initiative (WHI) genotyped on both the Metabochip and the Affymetrix 6.0 GWAS platform. We used this reference panel to impute into 6,459 WHI SNP Health Association Resource (SHARe) study subjects with only GWAS genotypes. Our analysis confirmed the imputation quality metric Rsq (estimated r(2) , specific to each SNP) as an effective post-imputation filter. We recommend different Rsq thresholds for different MAF categories such that the average (across SNPs) Rsq is above the desired dosage r(2) (squared Pearson correlation between imputed and experimental genotypes). With a desired dosage r(2) of 80%, 99.9% (97.5%, 83.6%, 52.0%, 20.5%) of SNPs with MAF > 0.05 (0.03-0.05, 0.01-0.03, 0.005-0.01, and 0.001-0.005) passed the post-imputation filter. The average dosage r(2) for these SNPs is 94.7%, 92.1%, 89.0%, 83.1%, and 79.7%, respectively. These results suggest that for African Americans imputation of Metabochip SNPs from GWAS data, including low frequency SNPs with MAF 0.005-0.05, is feasible and worthwhile for power increase in downstream association analysis provided a sizable reference panel is available.

TDT-HET: a new transmission disequilibrium test that incorporates locus heterogeneity into the analysis of family-based association data

BACKGROUND

Locus heterogeneity is one of the most documented phenomena in genetics. To date, relatively little work had been done on the development of methods to address locus heterogeneity in genetic association analysis. Motivated by Zhou and Pan's work, we present a mixture model of linked and unlinked trios and develop a statistical method to estimate the probability that a heterozygous parent transmits the disease allele at a di-allelic locus, and the probability that any trio is in the linked group. The purpose here is the development of a test that extends the classic transmission disequilibrium test (TDT) to one that accounts for locus heterogeneity.

RESULTS

Our simulations suggest that, for sufficiently large sample size (1000 trios) our method has good power to detect association even the proportion of unlinked trios is high (75%). While the median difference (TDT-HET empirical power - TDT empirical power) is approximately 0 for all MOI, there are parameter settings for which the power difference can be substantial. Our multi-locus simulations suggest that our method has good power to detect association as long as the markers are reasonably well-correlated and the genotype relative risk are larger. Results of both single-locus and multi-locus simulations suggest our method maintains the correct type I error rate.Finally, the TDT-HET statistic shows highly significant p-values for most of the idiopathic scoliosis candidate loci, and for some loci, the estimated proportion of unlinked trios approaches or exceeds 50%, suggesting the presence of locus heterogeneity.

CONCLUSIONS

We have developed an extension of the TDT statistic (TDT-HET) that allows for locus heterogeneity among coded trios. Benefits of our method include: estimates of parameters in the presence of heterogeneity, and reasonable power even when the proportion of linked trios is small. Also, we have extended multi-locus methods to TDT-HET and have demonstrated that the empirical power may be high to detect linkage. Last, given that we obtain PPBs, we conjecture that the TDT-HET may be a useful method for correctly identifying linked trios. We anticipate that researchers will find this property increasingly useful as they apply next-generation sequencing data in family based studies.

Development and validation of a short form of the valued life activities disability questionnaire for rheumatoid arthritis

OBJECTIVE

To develop and validate a shortened version of the Valued Life Activities disability and accommodations scale (VLA) for individuals with rheumatoid arthritis (RA).

METHODS

To shorten the existing VLA measure, item response theory analyses were conducted using data from 449 patients with RA. Next, the resulting 14-item shortened version of the VLA scale (S-VLA) was evaluated by structured interviews among 20 RA patients. Lastly, the S-VLA was administered to 150 RA patients along with other measures, including the Health Assessment Questionnaire (HAQ) and Short Form 36 (SF-36). A random sample of 50 patients completed the S-VLA 2 weeks later to assess reliability. Item statistics were calculated to evaluate correlations between individual items and the S-VLA total score. Correlations between the S-VLA and other measures were used to evaluate validity.

RESULTS

Test-retest reliability was 0.91, while Cronbach's alpha for the S-VLA was 0.95. None of the 14 items was associated with improved alpha coefficients when omitted. All of the items were strongly correlated with the S-VLA total score. S-VLA scores were highly positively correlated with the HAQ (r = 0.81, P ≤ 0.001), patient-reported disease activity (r = 0.71, P ≤ 0.001), satisfaction with abilities (r = 0.82, P ≤ 0.001), and number of days with activity limitations (r = 0.65, P ≤ 0.001). In addition, as hypothesized, the S-VLA was inversely correlated with the SF-36 physical component summary score (r = -0.78, P ≤ 0.001) and the physical functioning (r = -0.80, P ≤ 0.001), role physical (r = -0.67, P ≤ 0.001), and social functioning (r = -0.72, P ≤ 0.001) subscales.

CONCLUSION

The S-VLA is a short, valid, and reliable instrument that may prove useful for monitoring disability among individuals with RA.

The Next PAGE in understanding complex traits: design for the analysis of Population Architecture Using Genetics and Epidemiology (PAGE) Study

Genetic studies have identified thousands of variants associated with complex traits. However, most association studies are limited to populations of European descent and a single phenotype. The Population Architecture using Genomics and Epidemiology (PAGE) Study was initiated in 2008 by the National Human Genome Research Institute to investigate the epidemiologic architecture of well-replicated genetic variants associated with complex diseases in several large, ethnically diverse population-based studies. Combining DNA samples and hundreds of phenotypes from multiple cohorts, PAGE is well-suited to address generalization of associations and variability of effects in diverse populations; identify genetic and environmental modifiers; evaluate disease subtypes, intermediate phenotypes, and biomarkers; and investigate associations with novel phenotypes. PAGE investigators harmonize phenotypes across studies where possible and perform coordinated cohort-specific analyses and meta-analyses. PAGE researchers are genotyping thousands of genetic variants in up to 121,000 DNA samples from African-American, white, Hispanic/Latino, Asian/Pacific Islander, and American Indian participants. Initial analyses will focus on single nucleotide polymorphisms (SNPs) associated with obesity, lipids, cardiovascular disease, type 2 diabetes, inflammation, various cancers, and related biomarkers. PAGE SNPs are also assessed for pleiotropy using the “phenome-wide association study” approach, testing each SNP for associations with hundreds of phenotypes. PAGE data will be deposited into the National Center for Biotechnology Information's Database of Genotypes and Phenotypes and made available via a custom browser.

Genetic determinants of lipid traits in diverse populations from the population architecture using genomics and epidemiology (PAGE) study

For the past five years, genome-wide association studies (GWAS) have identified hundreds of common variants associated with human diseases and traits, including high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglyceride (TG) levels. Approximately 95 loci associated with lipid levels have been identified primarily among populations of European ancestry. The Population Architecture using Genomics and Epidemiology (PAGE) study was established in 2008 to characterize GWAS–identified variants in diverse population-based studies. We genotyped 49 GWAS–identified SNPs associated with one or more lipid traits in at least two PAGE studies and across six racial/ethnic groups. We performed a meta-analysis testing for SNP associations with fasting HDL-C, LDL-C, and ln(TG) levels in self-identified European American (∼20,000), African American (∼9,000), American Indian (∼6,000), Mexican American/Hispanic (∼2,500), Japanese/East Asian (∼690), and Pacific Islander/Native Hawaiian (∼175) adults, regardless of lipid-lowering medication use. We replicated 55 of 60 (92%) SNP associations tested in European Americans at p<0.05. Despite sufficient power, we were unable to replicate ABCA1 rs4149268 and rs1883025, CETP rs1864163, and TTC39B rs471364 previously associated with HDL-C and MAFB rs6102059 previously associated with LDL-C. Based on significance (p<0.05) and consistent direction of effect, a majority of replicated genotype-phentoype associations for HDL-C, LDL-C, and ln(TG) in European Americans generalized to African Americans (48%, 61%, and 57%), American Indians (45%, 64%, and 77%), and Mexican Americans/Hispanics (57%, 56%, and 86%). Overall, 16 associations generalized across all three populations. For the associations that did not generalize, differences in effect sizes, allele frequencies, and linkage disequilibrium offer clues to the next generation of association studies for these traits.

Low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and triglyceride (TG) levels are well known independent risk factors for cardiovascular disease. Lipid-associated genetic variants are being discovered in genome-wide association studies (GWAS) in samples of European descent, but an insufficient amount of data exist in other populations. Therefore, there is a strong need to characterize the effect of these GWAS–identified variants in more diverse cohorts. In this study, we selected over forty genetic loci previously associated with lipid levels and tested for replication in a large European American cohort. We also investigated if the effect of these variants generalizes to non-European descent populations, including African Americans, American Indians, and Mexican Americans/Hispanics. A majority of these GWAS–identified associations replicated in our European American cohort. However, the ability of associations to generalize across other racial/ethnic populations varied greatly, indicating that some of these GWAS–identified variants may not be functional and are more likely to be in linkage disequilibrium with the functional variant(s).

The use of phenome-wide association studies (PheWAS) for exploration of novel genotype-phenotype relationships and pleiotropy discovery

The field of phenomics has been investigating network structure among large arrays of phenotypes, and genome-wide association studies (GWAS) have been used to investigate the relationship between genetic variation and single diseases/outcomes. A novel approach has emerged combining both the exploration of phenotypic structure and genotypic variation, known as the phenome-wide association study (PheWAS). The Population Architecture using Genomics and Epidemiology (PAGE) network is a National Human Genome Research Institute (NHGRI)-supported collaboration of four groups accessing eight extensively characterized epidemiologic studies. The primary focus of PAGE is deep characterization of well-replicated GWAS variants and their relationships to various phenotypes and traits in diverse epidemiologic studies that include European Americans, African Americans, Mexican Americans/Hispanics, Asians/Pacific Islanders, and Native Americans. The rich phenotypic resources of PAGE studies provide a unique opportunity for PheWAS as each genotyped variant can be tested for an association with the wide array of phenotypic measurements available within the studies of PAGE, including prevalent and incident status for multiple common clinical conditions and risk factors, as well as clinical parameters and intermediate biomarkers. The results of PheWAS can be used to discover novel relationships between SNPs, phenotypes, and networks of interrelated phenotypes; identify pleiotropy; provide novel mechanistic insights; and foster hypothesis generation. The PAGE network has developed infrastructure to support and perform PheWAS in a high-throughput manner. As implementing the PheWAS approach has presented several challenges, the infrastructure and methodology, as well as insights gained in this project, are presented herein to benefit the larger scientific community.

Enhanced genetic maps from family-based disease studies: population-specific comparisons

Background

Accurate genetic maps are required for successful and efficient linkage mapping of disease genes. However, most available genome-wide genetic maps were built using only small collections of pedigrees, and therefore have large sampling errors. A large set of genetic studies genotyped by the NHLBI Mammalian Genotyping Service (MGS) provide appropriate data for generating more accurate maps.

Results

We collected a large sample of uncleaned genotype data for 461 markers generated by the MGS using the Weber screening sets 9 and 10. This collection includes genotypes for over 4,400 pedigrees containing over 17,000 genotyped individuals from different populations. We identified and cleaned numerous relationship and genotyping errors, as well as verified the marker orders. We used this dataset to test for population-specific genetic maps, and to re-estimate the genetic map distances with greater precision; standard errors for all intervals are provided. The map-interval sizes from the European (or European descent), Chinese, and Hispanic samples are in quite good agreement with each other. We found one map interval on chromosome 8p with a statistically significant size difference between the European and Chinese samples, and several map intervals with significant size differences between the African American and Chinese samples. When comparing Palauan with European samples, a statistically significant difference was detected at the telomeric region of chromosome 11p. Several significant differences were also identified between populations in chromosomal and genome lengths.

Conclusions

Our new population-specific screening set maps can be used to improve the accuracy of disease-mapping studies. As a result of the large sample size, the average length of the 95% confidence interval (CI) for a 10 cM map interval is only 2.4 cM, which is considerably smaller than on previously published maps.

Altered microRNA expression profiles in postmortem brain samples from individuals with schizophrenia and bipolar disorder

BACKGROUND

MicroRNAs (miRNAs) are potent regulators of gene expression with proposed roles in brain development and function. We hypothesized that miRNA expression profiles are altered in individuals with severe psychiatric disorders.

METHODS

With real-time quantitative polymerase chain reaction, we compared the expression of 435 miRNAs and 18 small nucleolar RNAs in postmortem brain tissue samples from individuals with schizophrenia, individuals with bipolar disorder, and psychiatrically healthy control subjects (n = 35 each group). Detailed demographic data, sample selection and storage conditions, and drug and substance exposure histories were available for all subjects. Bayesian model averaging was used to simultaneously assess the impact of these covariates as well as the psychiatric phenotype on miRNA expression profiles.

RESULTS

Of the variables considered, sample storage time, brain pH, alcohol at time of death, and postmortem interval were found to affect the greatest proportion of miRNAs. Of miRNAs analyzed, 19% exhibited positive evidence of altered expression due to a diagnosis of schizophrenia or bipolar disorder. Both conditions were associated with reduced miRNA expression levels, with a much more pronounced effect observed for bipolar disorder.

CONCLUSIONS

This study suggests that modest underexpression of several miRNAs might be involved in the complex pathogenesis of major psychosis.

Genetic variant of glutathione peroxidase 1 in autism

Genetic factors can contribute to autistic disorder (AD). Abnormal genes of oxidative stress pathways and increased oxidative stress have been reported in autism spectrum disorders. Polymorphisms of genes involved in glutathione metabolism, e.g. GSTP1 and GSTM1 are reportedly associated with autistic disorder. We investigated a GCG repeat polymorphism of a human glutathione peroxidase (GPX1) polyalanine repeat (ALA5, ALA6 and ALA7) in 103 trios of AD (probands and parents) using the transmission disequilibrium test. Significant transmission disequilibrium (p=0.044) was found in the overall transmission of the three alleles. The ALA6 allele was under transmitted (p=0.017). These results suggest that possessing this ALA6 allele may be protective for AD. Future study of interaction of the GPX1 GCG repeat and other gene polymorphisms such as the MnSOD ALA16 or the GPX1 Pro198Leu polymorphism in this cohort of AD families may shed light in whether the combination of the ALA6 allele with another polymorphism of antioxidant allele contributes to the increased oxidative stress in autism.