Replication of an association of variation in the FOXO3A gene with human longevity using both case–control and longitudinal data

Genetic variation in FOXO3A has previously been associated with human longevity. Studies published so far have been case–control studies and hence vulnerable to bias introduced by cohort effects. In this study we extended the previous findings in the cohorts of oldest old Danes (the Danish 1905 cohort, N = 1089) and middle‐aged Danes (N = 736), applying a longitudinal study design as well as the case–control study design. Fifteen SNPs were chosen in order to cover the known common variation in FOXO3A. Comparing SNP frequencies in the oldest old with middle‐aged individuals, we found association (after correction for multiple testing) of eight SNPs; 4 (rs13217795, rs2764264, rs479744, and rs9400239) previously reported to be associated with longevity and four novel SNPs (rs12206094, rs13220810, rs7762395, and rs9486902 (corrected P‐values 0.001–0.044). Moreover, we found association of the haplotypes TAC and CAC of rs9486902, rs10499051, and rs12206094 (corrected P‐values: 0.01–0.03) with longevity. Finally, we here present data applying a longitudinal study design; when using follow‐up survival data on the oldest old in a longitudinal analysis, we found no SNPs to remain significant after the correction for multiple testing (Bonferroni correction). Hence, our results support and extent the proposed role of FOXO3A as a candidate longevity gene for survival from younger ages to old age, yet not during old age.     

Evaluation of reduced subsets of single nucleotide polymorphisms for the prediction of age at puberty in sows

Genomic information could be used efficiently to improve traits that are expensive to measure, sex limited or expressed late in life. This study analyzed the phenotypic variation explained by major SNPs and windows for age at puberty in gilts, an indicator of reproductive longevity. A genome‐wide association study using 56 424 SNPs explained 25.2% of the phenotypic variation in age at puberty in a training set (n = 820). All SNPs from the top 10% of 1‐Mb windows explained 33.5% of the phenotypic variance compared to 47.1% explained by the most informative markers (n = 261). In an evaluation population, consisting of subsequent batches (n = 412), the predictive ability of all SNPs from the major 1‐Mb windows was higher compared to the variance captured by the most informative SNP from each of these windows. The phenotypic variance explained in the evaluation population varied from 12.3% to 36.8% when all SNPs from major windows were used compared to 6.5–23.7% explained by most informative SNPs. The correlation between phenotype and genomic prediction values based on SNP effects estimated in the training population was marginal compared to their effects retrained in the evaluation population for all (0.46–0.81) or most informative SNPs (0.30–0.65) from major windows. An increase in genetic gain of 20.5% could be obtained if genomic selection included both sexes compared to females alone. The pleiotropic role of major genes such as AVPR1A could be exploited in selection of both age at puberty and reproductive longevity.     

The IGF1 pathway genes and their association with age of puberty in cattle

Insulin‐like growth factor I (somatomedin C) (IGF1) influences gonadotrophin‐releasing hormone (GnRH) neurons during puberty, and GnRH release guides pubertal development. Therefore, genes of the IGF1 pathway are biological candidates for the identification of single‐nucleotide polymorphisms (SNPs) affecting age of puberty. In a genome‐wide association study, genotyped heifers were Tropical Composite (TCOMP, n = 866) or Brahman (BRAH, n = 843), with observation of age at first corpus luteum defining puberty. We examined SNPs in or near genes of the IGF1 pathway and report seven genes associated with age at puberty in cattle: IGF1R, IGFBP2, IGFBP4, PERK (HUGO symbol EIF2AK3), PIK3R1, GSK3B and IRS1. SNPs in the IGF1 receptor (IGF1R) showed the most promising associations: two SNPs were associated with puberty in TCOMP (P < 0.05) and one in BRAH (P = 0.00009). This last SNP explained 2% of the genetic variation (R² = 2.04%) for age of puberty in BRAH. Hence, IGF1R was examined further. Additional SNPs were genotyped, and haplotypes were analysed. To test more SNPs in this gene, four new SNPs from dbSNP were selected and genotyped. Single SNP and haploytpe analysis revealed associations with age of puberty in both breeds. There were two haplotypes of 12 IGF1R SNPs associated with puberty in BRAH (P < 0.05) and one in TCOMP (P < 0.05). One haplotype of two SNPs was associated (P < 0.01) with puberty in BRAH, but not in TCOMP. In conclusion, the IGF1 pathway appeared more relevant for age of puberty in Brahman cattle, and IGF1R showed higher significance when compared with other genes from the pathway.     

The FoxO3 gene and cause‐specific mortality

The G allele of the FOXO3 single nucleotide polymorphism (SNP) rs2802292 exhibits a consistently replicated genetic association with longevity in multiple populations worldwide. The aims of this study were to quantify the mortality risk for the longevity‐associated genotype and to discover the particular cause(s) of death associated with this allele in older Americans of diverse ancestry. It involved a 17‐year prospective cohort study of 3584 older American men of Japanese ancestry from the Honolulu Heart Program cohort, followed by a 17‐year prospective replication study of 1595 white and 1056 black elderly individuals from the Health Aging and Body Composition cohort. The relation between FOXO3 genotype and cause‐specific mortality was ascertained for major causes of death including coronary heart disease (CHD), cancer, and stroke. Age‐adjusted and multivariable Cox proportional hazards models were used to compute hazard ratios (HRs) for all‐cause and cause‐specific mortality. We found G allele carriers had a combined (Japanese, white, and black populations) risk reduction of 10% for total (all‐cause) mortality (HR = 0.90; 95% CI, 0.84–0.95; P = 0.001). This effect size was consistent across populations and mostly contributed by 26% lower risk for CHD death (HR = 0.74; 95% CI, 0.64–0.86; P = 0.00004). No other causes of death made a significant contribution to the survival advantage for G allele carriers. In conclusion, at older age, there is a large risk reduction in mortality for G allele carriers, mostly due to lower CHD mortality. The findings support further research on FOXO3 and FoxO3 protein as potential targets for therapeutic intervention in aging‐related diseases, particularly cardiovascular disease.     

The genetic component of human longevity: New insights from the analysis of pathway‐based SNP‐SNP interactions

In human longevity studies, single nucleotide polymorphism (SNP) analysis identified a large number of genetic variants with small effects, yet not easily replicable in different populations. New insights may come from the combined analysis of different SNPs, especially when grouped by metabolic pathway. We applied this approach to study the joint effect on longevity of SNPs belonging to three candidate pathways, the insulin/insulin‐like growth factor signalling (IIS), DNA repair and pro/antioxidant. We analysed data from 1,058 tagging SNPs in 140 genes, collected in 1825 subjects (1,089 unrelated nonagenarians from the Danish 1905 Birth Cohort Study and 736 Danish controls aged 46–55 years) for evaluating synergic interactions by SNPsyn. Synergies were further tested by the multidimensional reduction (MDR) approach, both intra‐ and interpathways. The best combinations (FDR<0.0001) resulted those encompassing IGF1R‐rs12437963 and PTPN1‐rs6067484, TP53‐rs2078486 and ERCC2‐rs50871, TXNRD1‐rs17202060 and TP53‐rs2078486, the latter two supporting a central role of TP53 in mediating the concerted activation of the DNA repair and pro‐antioxidant pathways in human longevity. Results were consistently replicated with both approaches, as well as a significant effect on longevity was found for the GHSR gene, which also interacts with partners belonging to both IIS and DNA repair pathways (PAPPA, PTPN1, PARK7, MRE11A). The combination GHSR‐MREA11, positively associated with longevity by MDR, was further found influencing longitudinal survival in nonagenarian females (p = .026). Results here presented highlight the validity of SNP‐SNP interactions analyses for investigating the genetics of human longevity, confirming previously identified markers but also pointing to novel genes as central nodes of additional networks involved in human longevity.     

Evaluation of polymorphisms in microRNA‐binding sites and pancreatic cancer risk in Chinese population

As promising biomarkers and therapy targets, microRNAs (miRNAs) are involved in various physiological and tumorigenic processes. Genetic variants in miRNA‐binding sites can lead to dysfunction of miRNAs and contribute to disease. However, systematic investigation of the miRNA‐related single nucleotide polymorphisms (SNPs) for pancreatic cancer (PC) risk remains elusive. We performed integrative bioinformatics analyses to select 31 SNPs located in miRNA‐target binding sites using the miRNASNP v2.0, a solid database providing miRNA‐related SNPs for genetic research, and investigated their associations with risk of PC in two large case‐control studies totally including 1847 cases and 5713 controls. We observed that the SNP rs3802266 is significantly associated with increased risk of PC (odds ratio (OR) = 1.21, 95% confidence intervals (CI) = 1.11‐1.31, P = 1.29E‐05). Following luciferase reporter gene assays show that rs3802266‐G creates a stronger binding site for miR‐181a‐2‐3p in 3′ untranslated region (3′UTR) of the gene ZHX2. Expression quantitative trait loci (eQTL) analysis suggests that ZHX2 expression is lower in individuals carrying rs3802266‐G with increased PC risk. In conclusion, our findings highlight the involvement of miRNA‐binding SNPs in PC susceptibility and provide new clues for PC carcinogenesis.     

Analysis of genetic variation contributing to measured speed in Thoroughbreds identifies genomic regions involved in the transcriptional response to exercise

Despite strong selection for athletic traits in Thoroughbred horses, there is marked variation in speed and aptitude for racing performance within the breed. Using global positioning system monitoring during exercise training, we measured speed variables and temporal changes in speed with age to derive phenotypes for GWAS. The aim of the study was to test the hypothesis that genetic variation contributes to variation in end‐point physiological traits, in this case galloping speed measured during field exercise tests. Standardisation of field‐measured phenotypes was attempted by assessing horses exercised on the same gallop track and managed under similar conditions by a single trainer. PCA of six key speed indices captured 73.9% of the variation with principal component 1 (PC1). Verifying the utility of the phenotype, we observed that PC1 (median) in 2‐year‐old horses was significantly different among elite, non‐elite and unraced horses (P < 0.001) and the temporal change with age in PC1 varied among horses with different myostatin (MSTN) g.66493737C>T SNP genotypes. A GWAS for PC1 in 2‐year‐old horses (n = 122) identified four SNPs reaching the suggestive threshold for association (P < 4.80 × 10^?5), defining a 1.09 Mb candidate region on ECA8 containing the myosin XVIIIB (MYO18B) gene. In a GWAS for temporal change in PC1 with age (n = 168), five SNPs reached the suggestive threshold for association and defined candidate regions on ECA2 and ECA11. Both regions contained genes that are significantly differentially expressed in equine skeletal muscle in response to acute exercise and training stimuli, including MYO18A. As MYO18A plays a regulatory role in the skeletal muscle response to exercise, the identified genomic variation proximal to the myosin family genes may be important for the regulation of the response to exercise and training.     

SNPs in SNCA,MCCC1, DLG2, GBF1 and MBNL2 are associated with Parkinson's disease in southern Chinese population

Numerous single nucleotide polymorphisms (SNPs), which have been identified as susceptibility factors for Parkinson's disease (PD) as per genome‐wide association studies, have not been fully characterized for PD patients in China. This study aimed to replicate the relationship between 12 novel SNPs of 12 genes and PD risk in southern Chinese population. Twelve SNPs of 12 genes were detected in 231 PD patients and 249 controls, using the SNaPshot technique. Meta‐analysis was used to assess heterogeneity of effect sizes between this study and published data. The impact of SNPs on gene expression was investigated by analysing the SNP‐gene association in the expression quantitative trait loci (eQTL) data sets. rs8180209 of SNCA (allele model: P = .047, OR = 0.77; additive model: P = .047, OR = 0.77), rs2270968 of MCCC1 (dominant model: P = .024, OR = 1.52), rs7479949 of DLG2 (recessive model; P = .019, OR = 1.52), rs10748818 of GBF1 (additive model: P < .001, OR = 0.37), and rs4771268 of MBNL2 (recessive model: P = .003, OR = 0.48) were replicated to be significantly associated with the increased risk of PD. Noteworthy, a meta‐analysis of previous studies suggested rs8180209, rs2270968, rs7479949 and rs4771268 were in line with those of our cohort. Our study replicated five novel functional SNPs in SNCA, MCCC1, DLG2, GBF1 and MBNL2 could be associated with increased risk of PD in southern Chinese population.     

Common genetic variants of the β2‐adrenergic receptor affect its translational efficiency and are associated with human longevity

β‐adrenoceptors are the common pharmacological targets for the treatment of cardiovascular diseases and asthma. Genetic modifications of β‐adrenergic system in engineered mice affect their lifespan. Here, we tested whether genes encoding for key components of the β‐adrenergic signaling pathway are associated with human longevity. We performed a 10‐year follow‐up study of the Chinese longitudinal healthy longevity survey. The Han Chinese population in this study consisted of 963 long‐lived and 1028 geography‐matched young individuals. Sixteen SNPs from ADRB1, ADRB2, ADCY5, ADCY6, and MAPK1 were selected and genotyped. Two SNPs, rs1042718 (C/A) and rs1042719 (G/C), of ADRB2 in linkage disequilibrium (D' = 1.0; r² = 0.67) were found to be associated with enhanced longevity in men in two geographically isolated populations. Bonferroni‐corrected P‐values in a combined analysis were 0.00053–0.010. Men with haplotype A‐C showed an increased probability to become centenarians (the frequency of A‐C in long‐lived and young individuals are 0.332 and 0.250, respectively, OR = 1.49, CI 95% = 1.17–1.88, P = 0.0007), in contrast to those with haplotype C‐G (the frequency of C‐G in long‐lived and young individuals are 0.523 and 0.635, respectively, OR = 0.63, CI 95% = 0.51–0.78, P = 0.000018). The permuted P‐values were 0.00005 and 0.0009, respectively. ADRB2 encodes the β2‐adrenergic receptor; the haplotype A‐C markedly reduced its translational efficiency compared with C‐G (P = 0.002) in transfected HEK293 cells. Thus, our data indicate that enhanced production of β2‐adrenergic receptors caused by genetic variants is inversely associated with human lifespan.     

Genome‐wide association study discovered candidate genes of Verticillium wilt resistance in upland cotton (Gossypium hirsutum L.)

Verticillium wilt (VW), caused by infection by Verticillium dahliae, is considered one of the most yield‐limiting diseases in cotton. To examine the genetic architecture of cotton VW resistance, we performed a genome‐wide association study (GWAS) using a panel of 299 accessions and 85 630 single nucleotide polymorphisms (SNPs) detected using the specific‐locus amplified fragment sequencing (SLAF‐seq) approach. Trait–SNP association analysis detected a total of 17 significant SNPs at P < 1.17 × 10^–5 (P = 1/85 630, –log₁₀P = 4.93); the peaks of SNPs associated with VW resistance on A10 were continuous and common in three environments (RDIG2015, RDIF2015 and RDIF2016). Haplotype block structure analysis predicted 22 candidate genes for VW resistance based on A10_99672586 with a minimum P‐value (–log₁₀P = 6.21). One of these genes (CG02) was near the significant SNP A10_99672586 (0.26 Mb), located in a 372‐kb haplotype block, and its Arabidopsis AT3G25510 homologues contain TIR‐NBS‐LRR domains that may be involved in disease resistance response. Real‐time quantitative PCR and virus‐induced gene silencing (VIGS) analysis showed that CG02 was specific to up‐regulation in the resistant (R) genotype Zhongzhimian2 (ZZM2) and that silenced plants were more susceptible to V. dahliae. These results indicate that CG02 is likely the candidate gene for resistance against V. dahliae in cotton. The identified locus or gene may serve as a promising target for genetic engineering and selection for improving resistance to VW in cotton.     

Genome‐wide association study of body morphological traits in Sudanese goats

Long‐term selection of goats for a certain production system and/or different environmental conditions will be reflected in the body morphology of the animals under selection. To investigate the variation contributing to different morphological traits and to identify genomic regions that are associated with body morphological traits in Sudanese goats, we genotyped 96 females belonging to four Sudanese goat breeds with the SNP52 BeadChip. After quality control of the data, the genome‐wide association study was performed using 95 goats and 24 027 informative single nucleotide polymorphisms (SNPs). Bicoastal diameter was significantly associated (LOD = 6.32) with snp10185‐scaffold1365‐620922 on chromosome 2. The minor allele has an additive effect, increasing the bicoastal diameter by 2.6 cm. A second significant association was found between body length and snp56482‐scaffold89‐467312 on chromosome 3 (LOD = 5.65). The minor allele is associated with increased body length. Additionally, five regions were suggestive for cannon bone, head width, rump length and withers height (LOD > 5). Only one gene (CNTNAP5) is located within the 1‐Mb region surrounding the significant SNP for bicoastal diameter on chromosome 2. The body length QTL on chromosome 3 harbors 49 genes. Further research is required to validate the observed associations and to prioritize candidate genes.     

A genome‐wide association analysis for carcass traits in a commercial Duroc pig population

We performed a genome‐wide association study to map the genetic determinants of carcass traits in 350 Duroc pigs typed with the Porcine SNP60 BeadChip. Association analyses were carried out using the gemma software. The proportion of phenotypic variance explained by the SNPs ranged between negligible to moderate (= 0.01–0.30) depending on the trait under consideration. At the genome‐wide level, we detected one significant association between backfat thickness between the 3^rd and 4^th ribs and six SNPs mapping to SSC12 (37–40 Mb). We also identified several chromosome‐wide significant associations for ham weight (SSC11: 51–53 Mb, three SNPs; 67–68 Mb, two SNPs), carcass weight (SSC11: 66–68 Mb, two SNPs), backfat thickness between the 3^rd and 4^th ribs (SSC12: 21 Mb, one SNP; 33–40 Mb, 17 SNPs; 51–58 Mb, two SNPs), backfat thickness in the last rib (SSC12: 37 Mb, one SNP; 40–41 Mb, nine SNPs) and lean meat content (SSC13: 34 Mb, three SNPs and SSC16: 45.1 Mb, one SNP; 62–63 Mb, 10 SNPs; 71–75 Mb, nine SNPs). The ham weight trait‐associated region on SSC11 contains two genes (UCHL3 and LMO7) related to muscle development. In addition, the ACACA gene, which encodes an enzyme for the catalysis of fatty acid synthesis, maps to the SSC12 (37–41 Mb) region harbouring trait‐associated regions for backfat thickness traits. Sequencing of these candidate genes may help to uncover the causal mutations responsible for the associations found in the present study.     

Copy number variations and polymorphisms in HSP90AB1 and risk of systemic lupus erythematosus and efficacy of glucocorticoids

The aim of our study was to assess the associations of HSP90AB1 copy number variations (CNVs) with systemic lupus erythematosus (SLE) risk and glucocorticoids (GCs) efficacy, as well as the relationship between HSP90AB1 single‐nucleotide polymorphisms (SNPs) and GCs efficacy. HSP90AB1 CNVs and SLE risk were analysed in 519 patients and 538 controls. Patients treated with GCs were followed up for 12 weeks and were divided into sensitive and insensitive groups to investigate the effects of CNVs (419 patients) and SNPs (457 patients) on the efficacy of GCs. Health‐related quality of life (HRQoL) was also measured by SF‐36 at baseline and week 12 to explore the relationship between CNVs/SNPs and HRQoL improvements in Chinese SLE patients. Our results indicated a statistically significant association between HSP90AB1 CNVs and SLE (P_BH = 0.039), and this association was more pronounced in the female subgroup (P_BH = 0.039). However, we did not detect association of HSP90AB1 CNVs/SNPs with efficacy of GCs. But we found a marginal association between SNP rs13296 and improvement in Role‐emotional, while this association was not strong enough to survive in the multiple testing corrections. Collectively, our findings suggest that the copy number of HSP90AB1 is associated with SLE susceptibility. But copy number and polymorphisms of HSP90AB1 may not be associated with efficacy of GCs.     

A large-scale candidate gene association study of age at menarche and age at natural menopause

Recent genome-wide association (GWA) studies have identified several novel genetic loci associated with age at menarche and age at natural menopause. However, the stringent significance threshold used in GWA studies potentially led to false negatives and true associations may have been overlooked. Incorporating biologically relevant information, we examined whether common genetic polymorphisms in candidate genes of nine groups of biologically plausible pathways and related phenotypes are associated with age at menarche and age at natural menopause. A total of 18,862 genotyped and imputed single nucleotide polymorphisms (SNPs) in 278 genes were assessed for their associations with these two traits among a total of 24,341 women from the Nurses’ Health Study (NHS, N = 2,287) and the Women’s Genome Health Study (WGHS, N = 22,054). Linear regression was used to assess the marginal association of each SNP with each phenotype. We adjusted for multiple testing within each gene to identify statistically significant SNP associations at the gene level. To evaluate the overall evidence for an excess of statistically significant gene associations over the proportion expected by chance, we applied a one-sample test of proportion to each group of candidate genes. The steroid-hormone metabolism and biosynthesis pathway was found significantly associated with both age at menarche and age at natural menopause (P = 0.040 and 0.011, respectively). In addition, the group of genes associated with precocious or delayed puberty was found significantly associated with age at menarche (P = 0.013), and the group of genes involved in premature ovarian failure with age at menopause (P = 0.025).