Pulmonary function measures reflect respiratory health and predict mortality, and are used in the diagnosis of chronic obstructive pulmonary disease (COPD). We tested genome-wide association with the forced expiratory volume in 1 second (FEV1) and the ratio of FEV1 to forced vital capacity (FVC) in 48,201 individuals of European ancestry, with follow-up of top associations in up to an additional 46,411 individuals. We identified new regions showing association (combined P<5×10−8) with pulmonary function, in or near MFAP2, TGFB2, HDAC4, RARB, MECOM (EVI1), SPATA9, ARMC2, NCR3, ZKSCAN3, CDC123, C10orf11, LRP1, CCDC38, MMP15, CFDP1, and KCNE2. Identification of these 16 new loci may provide insight into the molecular mechanisms regulating pulmonary function and into molecular targets for future therapy to alleviate reduced lung function.

Rationale: Genomic loci are associated with FEV1 or the ratio of FEV1 to FVC in population samples, but their association with chronic obstructive pulmonary disease (COPD) has not yet been proven, nor have their combined effects on lung function and COPD been studied.

Objectives: To test association with COPD of variants at five loci (TNS1, GSTCD, HTR4, AGER, and THSD4) and to evaluate joint effects on lung function and COPD of these single-nucleotide polymorphisms (SNPs), and variants at the previously reported locus near HHIP.

Methods: By sampling from 12 population-based studies (n = 31,422), we obtained genotype data on 3,284 COPD case subjects and 17,538 control subjects for sentinel SNPs in TNS1, GSTCD, HTR4, AGER, and THSD4. In 24,648 individuals (including 2,890 COPD case subjects and 13,862 control subjects), we additionally obtained genotypes for rs12504628 near HHIP. Each allele associated with lung function decline at these six SNPs contributed to a risk score. We studied the association of the risk score to lung function and COPD.

Measurements and Main Results: Association with COPD was significant for three loci (TNS1, GSTCD, and HTR4) and the previously reported HHIP locus, and suggestive and directionally consistent for AGER and TSHD4. Compared with the baseline group (7 risk alleles), carrying 10–12 risk alleles was associated with a reduction in FEV1 (β = –72.21 ml, P = 3.90 × 10−4) and FEV1/FVC (β = –1.53%, P = 6.35 × 10−6), and with COPD (odds ratio = 1.63, P = 1.46 × 10−5).

Conclusions: Variants in TNS1, GSTCD, and HTR4 are associated with COPD. Our highest risk score category was associated with a 1.6-fold higher COPD risk than the population average score.

Genome-wide association studies have identified 32 loci associated with body mass index (BMI), a measure that does not allow distinguishing lean from fat mass. To identify adiposity loci, we meta-analyzed associations between ~2.5 million SNPs and body fat percentage from 36,626 individuals, and followed up the 14 most significant (P<10−6) independent loci in 39,576 individuals. We confirmed the previously established adiposity locus in FTO (P=3×10−26), and identified two new loci associated with body fat percentage, one near IRS1 (P=4×10−11) and one near SPRY2 (P=3×10−8). Both loci harbour genes with a potential link to adipocyte physiology, of which the locus near IRS1 shows an intriguing association pattern. The body-fat-decreasing allele associates with decreased IRS1 expression and with an impaired metabolic profile, including decreased subcutaneous-to-visceral fat ratio, increased insulin resistance, dyslipidemia, risk of diabetes and coronary artery disease, and decreased adiponectin levels. Our findings provide new insights into adiposity and insulin resistance.

Evidence for a major role of genetic factors in the determination of body mass index (BMI) comes from studies of related individuals. Despite consistent evidence for a heritable component of BMI, estimates of BMI heritability vary widely between studies and the reasons for this remain unclear. While some variation is natural due to differences between populations and settings, study design factors may also explain some of the heterogeneity. We performed a systematic review that identified 88 independent estimates of BMI heritability from twin studies (total 140,525 twins) and 27 estimates from family studies (42,968 family members). BMI heritability estimates from twin studies ranged from 0.47 to 0.90 (5th/50th/95th centiles: 0.58/0.75/0.87) and were generally higher than those from family studies (range: 0.24–0.81; 5th/50th/95th centiles: 0.25/0.46/0.68). Meta-regression of the results from twin studies showed that BMI heritability estimates were 0.07 (P = 0.001) higher in children than in adults; estimates increased with mean age among childhood studies (+0.012/year, P = 0.002), but decreased with mean age in adult studies (−0.002/year, P = 0.002). Heritability estimates derived from AE twin models (which assume no contribution of shared environment) were 0.12 higher than those from ACE models (P < 0.001), whilst lower estimates were associated with self reported versus DNA-based determination of zygosity (−0.04, P = 0.02), and with self reported versus measured BMI (−0.05, P = 0.03). Although the observed differences in heritability according to aspects of study design are relatively small, together, the above factors explained 47% of the heterogeneity in estimates of BMI heritability from twin studies. In summary, while some variation in BMI heritability is expected due to population-level differences, study design factors explained nearly half the heterogeneity reported in twin studies. The genetic contribution to BMI appears to vary with age and may have a greater influence during childhood than adult life.

OBJECTIVE

Large-scale genome-wide association (GWA) studies have thus far identified 16 loci incontrovertibly associated with obesity-related traits in adults. We examined associations of variants in these loci with anthropometric traits in children and adolescents.

RESEARCH DESIGN AND METHODS

Seventeen variants representing 16 obesity susceptibility loci were genotyped in 1,252 children (mean ± SD age 9.7 ± 0.4 years) and 790 adolescents (15.5 ± 0.5 years) from the European Youth Heart Study (EYHS). We tested for association of individual variants and a genetic predisposition score (GPS-17), calculated by summing the number of effect alleles, with anthropometric traits. For 13 variants, summary statistics for associations with BMI were meta-analyzed with previously reported data (Ntotal = 13,071 children and adolescents).

RESULTS

In EYHS, 15 variants showed associations or trends with anthropometric traits that were directionally consistent with earlier reports in adults. The meta-analysis showed directionally consistent associations with BMI for all 13 variants, of which 9 were significant (0.033–0.098 SD/allele; P < 0.05). The near-TMEM18 variant had the strongest effect (0.098 SD/allele P = 8.5 × 10−11). Effect sizes for BMI tended to be more pronounced in children and adolescents than reported earlier in adults for variants in or near SEC16B, TMEM18, and KCTD15, (0.028–0.035 SD/allele higher) and less pronounced for rs925946 in BDNF (0.028 SD/allele lower). Each additional effect allele in the GPS-17 was associated with an increase of 0.034 SD in BMI (P = 3.6 × 10−5), 0.039 SD, in sum of skinfolds (P = 1.7 × 10−7), and 0.022 SD in waist circumference (P = 1.7 × 10−4), which is comparable with reported results in adults (0.039 SD/allele for BMI and 0.033 SD/allele for waist circumference).

CONCLUSIONS

Most obesity susceptibility loci identified by GWA studies in adults are already associated with anthropometric traits in children/adolescents. Whereas the association of some variants may differ with age, the cumulative effect size is similar.

The timing of puberty is highly variable1. We carried out a genome-wide association study for age at menarche in 4,714 women and report an association in LIN28B on chromosome 6 (rs314276, minor allele frequency (MAF) = 0.33, P = 1.5 × 10−8). In independent replication studies in 16,373 women, each major allele was associated with 0.12 years earlier menarche (95% CI = 0.08–0.16; P = 2.8 × 10−10; combined P = 3.6 × 10−16). This allele was also associated with earlier breast development in girls (P = 0.001; N = 4,271); earlier voice breaking (P = 0.006, N = 1,026) and more advanced pubic hair development in boys (P = 0.01; N = 4,588); a faster tempo of height growth in girls (P = 0.00008; N = 4,271) and boys (P = 0.03; N = 4,588); and shorter adult height in women (P = 3.6 × 10−7; N = 17,274) and men (P = 0.006; N = 9,840) in keeping with earlier growth cessation. These studies identify variation in LIN28B, a potent and specific regulator of microRNA processing2, as the first genetic determinant regulating the timing of human pubertal growth and development.

Recently, an indirect genetic association approach that compares genotype frequencies in offspring of long-lived subjects and offspring from random families has been introduced to study gene-longevity associations. Although the indirect genetic association has certain advantages over the direct association approach that compares genotype frequency between centenarians and young controls, the power has been of concern. This paper reports a power study performed on the indirect approach using computer simulation. We perform our simulation study by introducing the current Danish population life table and the proportional hazard model for generating individual lifespan. Family genotype data is generated using a genetic linkage program for given SNP allele frequency. Power is estimated by setting the type I error rate at 0.05 and by calculating the Armitage’s chi-squared test statistic for 200 replicate samples for each setting of the specified allele risk and frequency parameters under different modes of inheritance and for different sample sizes. The indirect genetic association analysis is a valid approach for studying gene-longevity association, but the sample size requirement is about 3–4 time larger than the direct approach. It also has low power in detecting non-additive effect genes. Indirect genetic association using offspring from families with both parents as nonagenarians is nearly as powerful as using offspring from families with one centenarian parent. In conclusion, the indirect design can be a good choice for studying longevity in comparison with other alternatives, when relatively large sample size is available.

Elevated blood pressure is a common, heritable cause of cardiovascular disease worldwide. To date, identification of common genetic variants influencing blood pressure has proven challenging. We tested 2.5m genotyped and imputed SNPs for association with systolic and diastolic blood pressure in 34,433 subjects of European ancestry from the Global BPgen consortium and followed up findings with direct genotyping (N≤71,225 European ancestry, N=12,889 Indian Asian ancestry) and in silico comparison (CHARGE consortium, N=29,136). We identified association between systolic or diastolic blood pressure and common variants in 8 regions near the CYP17A1 (P=7×10−24), CYP1A2 (P=1×10−23), FGF5 (P=1×10−21), SH2B3 (P=3×10−18), MTHFR (P=2×10−13), c10orf107 (P=1×10−9), ZNF652 (P=5×10−9) and PLCD3 (P=1×10−8) genes. All variants associated with continuous blood pressure were associated with dichotomous hypertension. These associations between common variants and blood pressure and hypertension offer mechanistic insights into the regulation of blood pressure and may point to novel targets for interventions to prevent cardiovascular disease.

Shengxu Li and colleagues use data from a large prospective observational cohort to examine the extent to which a genetic predisposition toward obesity may be modified by living a physically active lifestyle.

Background

We have previously shown that multiple genetic loci identified by genome-wide association studies (GWAS) increase the susceptibility to obesity in a cumulative manner. It is, however, not known whether and to what extent this genetic susceptibility may be attenuated by a physically active lifestyle. We aimed to assess the influence of a physically active lifestyle on the genetic predisposition to obesity in a large population-based study.

Methods and Findings

We genotyped 12 SNPs in obesity-susceptibility loci in a population-based sample of 20,430 individuals (aged 39–79 y) from the European Prospective Investigation of Cancer (EPIC)-Norfolk cohort with an average follow-up period of 3.6 y. A genetic predisposition score was calculated for each individual by adding the body mass index (BMI)-increasing alleles across the 12 SNPs. Physical activity was assessed using a self-administered questionnaire. Linear and logistic regression models were used to examine main effects of the genetic predisposition score and its interaction with physical activity on BMI/obesity risk and BMI change over time, assuming an additive effect for each additional BMI-increasing allele carried. Each additional BMI-increasing allele was associated with 0.154 (standard error [SE] 0.012) kg/m2 (p = 6.73×10−37) increase in BMI (equivalent to 445 g in body weight for a person 1.70 m tall). This association was significantly (pinteraction = 0.005) more pronounced in inactive people (0.205 [SE 0.024] kg/m2 [p = 3.62×10−18; 592 g in weight]) than in active people (0.131 [SE 0.014] kg/m2 [p = 7.97×10−21; 379 g in weight]). Similarly, each additional BMI-increasing allele increased the risk of obesity 1.116-fold (95% confidence interval [CI] 1.093–1.139, p = 3.37×10−26) in the whole population, but significantly (pinteraction = 0.015) more in inactive individuals (odds ratio [OR] = 1.158 [95% CI 1.118–1.199; p = 1.93×10−16]) than in active individuals (OR = 1.095 (95% CI 1.068–1.123; p = 1.15×10−12]). Consistent with the cross-sectional observations, physical activity modified the association between the genetic predisposition score and change in BMI during follow-up (pinteraction = 0.028).

Conclusions

Our study shows that living a physically active lifestyle is associated with a 40% reduction in the genetic predisposition to common obesity, as estimated by the number of risk alleles carried for any of the 12 recently GWAS-identified loci.

Please see later in the article for the Editors' Summary

Editors' Summary

Background

In the past few decades, the global incidence of obesity—defined as a body mass index (BMI, a simple index of weight-for-height that uses the weight in kilograms divided by the square of the height in meters) of 30 and over, has increased so much that this growing public health concern is now commonly referred to as the “obesity epidemic.” Once considered prevalent only in high-income countries, obesity is an increasing health problem in low- and middle-income countries, particularly in urban settings. In 2005, at least 400 million adults world-wide were obese, and the projected figure for 2015 is a substantial increase of 300 million to around 700 million. Childhood obesity is also a growing concern. Contributing factors to the obesity epidemic are a shift in diet to an increased intake of energy-dense foods that are high in fat and sugars and a trend towards decreased physical activity due to increasingly sedentary lifestyles.

However, genetics are also thought to play a critical role as genetically predisposed individuals may be more prone to obesity if they live in an environment that has abundant access to energy-dense food and labor-saving devices.

Why Was This Study Done?

Although recent genetic studies (genome-wide association studies) have identified 12 alleles (a DNA variant that is located at a specific position on a specific chromosome) associated with increased BMI, there has been no convincing evidence of the interaction between genetics and lifestyle. In this study the researchers examined the possibility of such an interaction by assessing whether individuals with a genetic predisposition to increased obesity risk could modify this risk by increasing their daily physical activity.

What Did the Researchers Do and Find?

The researchers used a population-based cohort study of 25,631 people living in Norwich, UK (The EPIC-Norfolk study) and identified individuals who were 39 to 79 years old during a health check between 1993 and 1997. The researchers invited these people to a second health examination. In total, 20,430 individuals had baseline data available, of which 11,936 had BMI data at the second health check. The researchers used genotyping methods and then calculated a genetic predisposition score for each individual and their occupational and leisure-time physical activities were assessed by using a validated self-administered questionnaire. Then, the researchers used modeling techniques to examine the main effects of the genetic predisposition score and its interaction with physical activity on BMI/obesity risk and BMI change over time. The researchers found that each additional BMI-increasing allele was associated with an increase in BMI equivalent to 445 g in body weight for a person 1.70 m tall and that the size of this effect was greater in inactive people than in active people. In individuals who have a physically active lifestyle, this increase was only 379 g/allele, or 36% lower than in physically inactive individuals in whom the increase was 592 g/allele. Furthermore, in the total sample each additional obesity-susceptibility allele increased the odds of obesity by 1.116-fold. However, the increased odds per allele for obesity risk were 40% lower in physically active individuals (1.095 odds/allele) compared to physically inactive individuals (1.158 odds/allele).

What Do These Findings Mean?

The findings of this study indicate that the genetic predisposition to obesity can be reduced by approximately 40% by having a physically active lifestyle. The findings of this study suggest that, while the whole population benefits from increased physical activity levels, individuals who are genetically predisposed to obesity would benefit more than genetically protected individuals. Furthermore, these findings challenge the deterministic view of the genetic predisposition to obesity that is often held by the public, as they show that even the most genetically predisposed individuals will benefit from adopting a healthy lifestyle. The results are limited by participants self-reporting their physical activity levels, which is less accurate than objective measures of physical activity.

Additional Information

Please access these Web sites via the online version of this summary at http://dx.doi.org/10.1371/journal.pmed.1000332.

This study relies on the results of previous genome-wide association studies The National Human Genome Research Institute provides an easy-to-follow guide to understanding such studies

The International Association for the Study of Obesity aims to improve global health by promoting the understanding of obesity and weight-related diseases through scientific research and dialogue

The International Obesity Taskforce is the research-led think tank and advocacy arm of the International Association for the Study of Obesity

The Global Alliance for the Prevention of Obesity and Related Chronic Disease is a global action program that addresses the issues surrounding the prevention of obesity

The National Institutes of Health has its own obesity task force, which includes 26 institutes

Pulmonary function measures are heritable traits that predict morbidity and mortality and define chronic obstructive pulmonary disease (COPD). We tested genome-wide association with forced expiratory volume in 1 s (FEV1) and the ratio of FEV1 to forced vital capacity (FVC) in the SpiroMeta consortium (n = 20,288 individuals of European ancestry). We conducted a meta-analysis of top signals with data from direct genotyping (n ≤ 32,184 additional individuals) and in silico summary association data from the CHARGE Consortium (n = 21,209) and the Health 2000 survey (n ≤ 883). We confirmed the reported locus at 4q31 and identified associations with FEV1 or FEV1/FVC and common variants at five additional loci: 2q35 in TNS1 (P = 1.11 × 10−12), 4q24 in GSTCD (2.18 × 10−23), 5q33 in HTR4 (P = 4.29 × 10−9), 6p21 in AGER (P = 3.07 × 10−15) and 15q23 in THSD4 (P = 7.24 × 10−15). mRNA analyses showed expression of TNS1, GSTCD, AGER, HTR4 and THSD4 in human lung tissue. These associations offer mechanistic insight into pulmonary function regulation and indicate potential targets for interventions to alleviate respiratory disease.

Mutation within the Leucine-rich repeat kinase 2 (LRRK2) gene has been identified as a cause of autosomal dominant Parkinson’s disease (PD). The purpose of this study was to determine the frequency of G2019S mutation and whether the differences in the allele and genotype distribution of six SNPs within LRRK2 gene are associated with PD in an American non-Hispanic white population. The sample included 350 sporadic PD (SPD), 225 familial PD (FPD) patients and 186 controls of the same race and ethnicity. The frequency of LRRK2 G2019S mutation in our total sample of PD (FPD and SPD) was 1.56%. The frequency of this mutation was 3.5% in the FPD and 0.3% in the SPD groups, respectively. Allele and genotype frequencies of six SNPs were compared between PD and control samples. In addition, PD groups were categorized by sporadic PD (no family history), familial PD (first degree relative with PD) and age of onset (AON, ≤50 or ≥51 years). The haplotypes of the six SNPs were also constructed for association analysis. After correction for multiple comparisons, there was no association between any SNPs (allele or genotype) and PD groups. One of the haplotypes was modestly associated with the combined PD (SPD & FPD) sample. There was also no association with age at onset of PD. Our study suggests that the LRRK2 gene may be a risk factor or the cause for a very small fraction of PD in American white population.

Obesity has become an epidemic in many countries and is one of the major risk conditions for disease including type 2 diabetes, coronary heart disease, stroke, dyslipidemia, and hypertension. Recent genome-wide association studies have identified two genes (FTO and near MC4R) that were unequivocally associated with body mass index (BMI) and obesity. For the Genetic Analysis Workshop 16, data from the Framingham Heart Study were made available, including longitudinal anthropometric and metabolic traits for 7130 Caucasian individuals over three generations, each with follow-up data at up to four time points. We explored the associations between four single-nucleotide polymorphisms (SNPs) on FTO (rs1121980, rs9939609) or near MC4R (rs17782313, rs17700633) with weight and BMI under an additive model. We applied multilevel linear mixed model for continuous outcomes, using the Affymetrix 500k genome-wide genotype data for the four SNPs. The results of the multilevel modeling in the entire sample indicated that the minor alleles of the four SNPs were associated with higher weight and higher BMI. The most significant associations were between rs9939609 and weight (p = 4.7 × 10-6) and BMI (p = 8.9 × 10-8). The results also showed that, for SNPs at FTO, the homozygotes for the minor allele had the most pronounced increase in weight and BMI, while the common allele homozygotes gained less weight and BMI during the follow-up period. Linkage disequilibrium (LD) between the two FTO SNPs was strong (D' = 0.997, r2 = 0.875) but their haplotype was not significantly associated with either weight or BMI. The two SNPs near MC4R were in weak LD (D' = 0.487, r2 = 0.166).

The efficiency of the popular case-control design in gene-longevity association studies needs to be verified because, different from a binary trait, longevity represents only the extreme end of the continuous life span distribution without a clear cutoff for defining the phenotype. In this paper, the authors use the current Danish life tables to simulate individual life span by using a variety of scenarios and assess the empirical power for different sample sizes when cases are defined as centenarians or as nonagenarians. Results show that, although using small samples of centenarians (several hundred) provides power to detect only common alleles with large effects (a >20% reduction in hazard rate), large samples of centenarians (>1,000) achieve power to capture genes responsible for minor effects (5%–10% hazard reduction depending on the mode of inheritance). Although the method provides good power for rare alleles with multiplicative or dominant effects, it performs poorly for rare recessive alleles. Power is drastically reduced when nonagenarians are considered cases, with a more than 5-fold difference in the size of the case sample required to achieve comparable power as that found with centenarians.

To identify common variants influencing body mass index (BMI), we analyzed genome-wide association data from 16,876 individuals of European descent. After previously reported variants in FTO, the strongest association signal (rs17782313, P = 2.9 × 10−6) mapped 188 kb downstream of MC4R (melanocortin-4 receptor), mutations of which are the leading cause of monogenic severe childhood-onset obesity. We confirmed the BMI association in 60,352 adults (per-allele effect = 0.05 Z-score units; P = 2.8 × 10−15) and 5,988 children aged 7–11 (0.13 Z-score units; P = 1.5 × 10−8). In case-control analyses (n = 10,583), the odds for severe childhood obesity reached 1.30 (P = 8.0 × 10−11). Furthermore, we observed overtransmission of the risk allele to obese offspring in 660 families (P (pedigree disequilibrium test average; PDT-avg) = 2.4 × 10−4). The SNP location and patterns of phenotypic associations are consistent with effects mediated through altered MC4R function. Our findings establish that common variants near MC4R influence fat mass, weight and obesity risk at the population level and reinforce the need for large-scale data integration to identify variants influencing continuous biomedical traits.

Summary

Background

LDL cholesterol has a causal role in the development of cardiovascular disease. Improved understanding of the biological mechanisms that underlie the metabolism and regulation of LDL cholesterol might help to identify novel therapeutic targets. We therefore did a genome-wide association study of LDL-cholesterol concentrations.

Methods

We used genome-wide association data from up to 11 685 participants with measures of circulating LDL-cholesterol concentrations across five studies, including data for 293 461 autosomal single nucleotide polymorphisms (SNPs) with a minor allele frequency of 5% or more that passed our quality control criteria. We also used data from a second genome-wide array in up to 4337 participants from three of these five studies, with data for 290 140 SNPs. We did replication studies in two independent populations consisting of up to 4979 participants. Statistical approaches, including meta-analysis and linkage disequilibrium plots, were used to refine association signals; we analysed pooled data from all seven populations to determine the effect of each SNP on variations in circulating LDL-cholesterol concentrations.

Findings

In our initial scan, we found two SNPs (rs599839 [p=1·7×10−15] and rs4970834 [p=3·0×10−11]) that showed genome-wide statistical association with LDL cholesterol at chromosomal locus 1p13.3. The second genome screen found a third statistically associated SNP at the same locus (rs646776 [p=4·3×10−9]). Meta-analysis of data from all studies showed an association of SNPs rs599839 (combined p=1·2×10−33) and rs646776 (p=4·8×10−20) with LDL-cholesterol concentrations. SNPs rs599839 and rs646776 both explained around 1% of the variation in circulating LDL-cholesterol concentrations and were associated with about 15% of an SD change in LDL cholesterol per allele, assuming an SD of 1 mmol/L.

Interpretation

We found evidence for a novel locus for LDL cholesterol on chromosome 1p13.3. These results potentially provide insight into the biological mechanisms that underlie the regulation of LDL cholesterol and might help in the discovery of novel therapeutic targets for cardiovascular disease.

The Genetic Analysis Workshop 15 (GAW15) Problem 1 contained baseline expression levels of 8793 genes in immortalized B cells from 194 individuals in 14 Centre d'Etude du Polymorphisme Humain (CEPH) Utah pedigrees. Previous analysis of the data showed linkage and association and evidence of substantial individual variations. In particular, correlation was examined on expression levels of 31 genes and 25 target genes corresponding to two master regulatory regions. In this analysis, we apply Bayesian network analysis to gain further insight into these findings. We identify strong dependences and therefore provide additional insight into the underlying relationships between the genes involved. More generally, the approach is expected to be applicable for integrated analysis of genes on biological pathways.

Background

The etiology of multifactorial human diseases involves complex interactions between numerous environmental factors and alleles of many genes. Efficient statistical tools are demanded in identifying the genetic and environmental variants that affect the risk of disease development. This paper introduces a retrospective polytomous logistic regression model to measure both the main and interaction effects in genetic association studies of human discrete and continuous complex traits. In this model, combinations of genotypes at two interacting loci or of environmental exposure and genotypes at one locus are treated as nominal outcomes of which the proportions are modeled as a function of the disease trait assigning both main and interaction effects and with no assumption of normality in the trait distribution. Performance of our method in detecting interaction effect is compared with that of the case-only model.

Results

Results from our simulation study indicate that our retrospective model exhibits high power in capturing even relatively small effect with reasonable sample sizes. Application of our method to data from an association study on the catalase -262C/T promoter polymorphism and aging phenotypes detected significant main and interaction effects for age-group and allele T on individual's cognitive functioning and produced consistent results in estimating the interaction effect as compared with the popular case-only model.

Conclusion

The retrospective polytomous logistic regression model can be used as a convenient tool for assessing both main and interaction effects in genetic association studies of human multifactorial diseases involving genetic and non-genetic factors as well as categorical or continuous traits.

We examine the efficiency of a number of schemes to select cases from nuclear families for case-control association analysis using the Genetic Analysis Workshop 14 simulated dataset. We show that with this simulated dataset comparing all affected siblings with unrelated controls is considerably more powerful than all of the other approaches considered. We find that the test statistic is increased by almost 3-fold compared to the next best sampling schemes of selecting all affected sibs only from families with affected parents (AFaff), one affected sib with most evidence of allele-sharing from each family (SF), and all affected sibs from families with evidence for linkage (AFL). We consider accounting for biological relatedness of samples in the association analysis to maintain the correct type I error. We also discuss the relative efficiencies of increasing the ratio of unrelated cases to controls, methods to confirm associations and issues to consider when applying our conclusions to other complex disease datasets.

The presence of disease is commonly used in genetic studies; however, the time to onset often provides additional information. To apply the popular Cox model for such data, it is desirable to consider the familial correlation, which involves kinship or identity by descent (IBD) information between family members. Recently, such a framework has been developed and implemented in a UNIX-based S-PLUS package called kinship, extending the Cox model with mixed effects and familial relationship. The model is of great potential in joint analysis of family data with genetic and environmental factors. We apply this framework to data from the Collaborative Study on the Genetics of Alcoholism data as part of Genetic Analysis Workshop 14. We use the S-PLUS package, ported into the R environment , for the analysis of microsatellite data on chromosomes 4 and 7. In these analyses, IBD information at those markers is used in addition to the basic Cox model with mixed effects, which provides estimates of the relative contribution of specific genetic markers. D4S1645 had the largest variance and contribution to the log-likelihood on chromosome 4, but the significance of this finding requires further investigation.

Genetic data are now widely available. There is, however, an apparent lack of concerted effort to produce software systems for statistical analysis of genetic data compared with other fields of statistics. It is often a tremendous task for end-users to tailor them for particular data, especially when genetic data are analysed in conjunction with a large number of covariates. Here, R http://www.r-project.org, a free, flexible and platform-independent environment for statistical modelling and graphics is explored as an integrated system for genetic data analysis. An overview of some packages currently available for analysis of genetic data is given. This is followed by examples of package development and practical applications. With clear advantages in data management, graphics, statistical analysis, programming, internet capability and use of available codes, it is a feasible, although challenging, task to develop it into an integrated platform for genetic analysis; this will require the joint efforts of many researchers.

Genome-wide association studies have identified numerous genetic loci for spirometic measures of pulmonary function, forced expiratory volume in one second (FEV1), and its ratio to forced vital capacity (FEV1/FVC). Given that cigarette smoking adversely affects pulmonary function, we conducted genome-wide joint meta-analyses (JMA) of single nucleotide polymorphism (SNP) and SNP-by-smoking (ever-smoking or pack-years) associations on FEV1 and FEV1/FVC across 19 studies (total N = 50,047). We identified three novel loci not previously associated with pulmonary function. SNPs in or near DNER (smallest PJMA = 5.00×10−11), HLA-DQB1 and HLA-DQA2 (smallest PJMA = 4.35×10−9), and KCNJ2 and SOX9 (smallest PJMA = 1.28×10−8) were associated with FEV1/FVC or FEV1 in meta-analysis models including SNP main effects, smoking main effects, and SNP-by-smoking (ever-smoking or pack-years) interaction. The HLA region has been widely implicated for autoimmune and lung phenotypes, unlike the other novel loci, which have not been widely implicated. We evaluated DNER, KCNJ2, and SOX9 and found them to be expressed in human lung tissue. DNER and SOX9 further showed evidence of differential expression in human airway epithelium in smokers compared to non-smokers. Our findings demonstrated that joint testing of SNP and SNP-by-environment interaction identified novel loci associated with complex traits that are missed when considering only the genetic main effects.

Author Summary

Measures of pulmonary function provide important clinical tools for evaluating lung disease and its progression. Genome-wide association studies have identified numerous genetic risk factors for pulmonary function but have not considered interaction with cigarette smoking, which has consistently been shown to adversely impact pulmonary function. In over 50,000 study participants of European descent, we applied a recently developed joint meta-analysis method to simultaneously test associations of gene and gene-by-smoking interactions in relation to two major clinical measures of pulmonary function. Using this joint method to incorporate genetic main effects plus gene-by-smoking interaction, we identified three novel gene regions not previously related to pulmonary function: (1) DNER, (2) HLA-DQB1 and HLA-DQA2, and (3) KCNJ2 and SOX9. Expression analyses in human lung tissue from ours or prior studies indicate that these regions contain genes that are plausibly involved in pulmonary function. This work highlights the utility of employing novel methods for incorporating environmental interaction in genome-wide association studies to identify novel genetic regions.

Rationale

Asthma has substantial morbidity and mortality and a strong genetic component, but identification of genetic risk factors is limited by availability of suitable studies.

Objectives

To test if population-based cohorts with self-reported physician-diagnosed asthma and genome-wide association (GWA) data could be used to validate known associations with asthma and identify novel associations.

Methods

The APCAT (Analysis in Population-based Cohorts of Asthma Traits) consortium consists of 1,716 individuals with asthma and 16,888 healthy controls from six European-descent population-based cohorts. We examined associations in APCAT of thirteen variants previously reported as genome-wide significant (P<5x10−8) and three variants reported as suggestive (P<5×10−7). We also searched for novel associations in APCAT (Stage 1) and followed-up the most promising variants in 4,035 asthmatics and 11,251 healthy controls (Stage 2). Finally, we conducted the first genome-wide screen for interactions with smoking or hay fever.

Main Results

We observed association in the same direction for all thirteen previously reported variants and nominally replicated ten of them. One variant that was previously suggestive, rs11071559 in RORA, now reaches genome-wide significance when combined with our data (P = 2.4×10−9). We also identified two genome-wide significant associations: rs13408661 near IL1RL1/IL18R1 (PStage1+Stage2 = 1.1x10−9), which is correlated with a variant recently shown to be associated with asthma (rs3771180), and rs9268516 in the HLA region (PStage1+Stage2 = 1.1x10−8), which appears to be independent of previously reported associations in this locus. Finally, we found no strong evidence for gene-environment interactions with smoking or hay fever status.

Conclusions

Population-based cohorts with simple asthma phenotypes represent a valuable and largely untapped resource for genetic studies of asthma.

Numerous genetic loci influence systolic blood pressure (SBP) and diastolic blood pressure (DBP) in Europeans 1-3. We now report genome-wide association studies of pulse pressure (PP) and mean arterial pressure (MAP). In discovery (N=74,064) and follow-up studies (N=48,607), we identified at genome-wide significance (P= 2.7×10-8 to P=2.3×10-13) four novel PP loci (at 4q12 near CHIC2/PDGFRAI, 7q22.3 near PIK3CG, 8q24.12 in NOV, 11q24.3 near ADAMTS-8), two novel MAP loci (3p21.31 in MAP4, 10q25.3 near ADRB1) and one locus associated with both traits (2q24.3 near FIGN) which has recently been associated with SBP in east Asians. For three of the novel PP signals, the estimated effect for SBP was opposite to that for DBP, in contrast to the majority of common SBP- and DBP-associated variants which show concordant effects on both traits. These findings indicate novel genetic mechanisms underlying blood pressure variation, including pathways that may differentially influence SBP and DBP.

Genetic association studies are now routinely used to identify single nucleotide polymorphisms (SNPs) linked with human diseases or traits through single SNP-single trait tests. Here we introduced partial least squares path modeling (PLSPM) for association between single or multiple SNPs and a latent trait that can involve single or multiple correlated measurement(s). Furthermore, the framework naturally provides estimators of polygenic effect by appropriately weighting trait-attributing alleles. We conducted computer simulations to assess the performance via multiple SNPs and human obesity-related traits as measured by body mass index (BMI), waist and hip circumferences. Our results showed that the associate statistics had type I error rates close to nominal level and were powerful for a range of effect and sample sizes. When applied to 12 candidate regions in data (N = 2,417) from the European Prospective Investigation of Cancer (EPIC)-Norfolk study, a region in FTO was found to have stronger association (rs7204609∼rs9939881 at the first intron P = 4.29×10−7) than single SNP analysis (all with P>10−4) and a latent quantitative phenotype was obtained using a subset sample of EPIC-Norfolk (N = 12,559). We believe our method is appropriate for assessment of regional association and polygenic effect on a single or multiple traits.

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−60) and 9q31.2 (P=2.2×10−33), we identified 30 novel menarche loci (all P<5×10−8) and found suggestive evidence for a further 10 loci (P<1.9×10−6). New loci included four previously associated with BMI (in/near FTO, SEC16B, TRA2B and TMEM18), three in/near other genes implicated in energy homeostasis (BSX, CRTC1, and MCHR2), and three in/near genes implicated in hormonal regulation (INHBA, PCSK2 and RXRG). Ingenuity and MAGENTA pathway analyses identified coenzyme A and fatty acid biosynthesis as biological processes related to menarche timing.