Fast-evolving noncoding sequences in the human genome

BACKGROUND: Gene regulation is considered one of the driving forces of evolution. Although protein-coding DNA sequences and RNA genes have been subject to recent evolutionary events in the human lineage, it has been hypothesized that the large phenotypic divergence between humans and chimpanzees has been driven mainly by changes in gene regulation rather than altered protein-coding gene sequences. Comparative analysis of vertebrate genomes has revealed an abundance of evolutionarily conserved but noncoding sequences. These conserved noncoding (CNC) sequences may well harbor critical regulatory variants that have driven recent human evolution. RESULTS: Here we identify 1,356 CNC sequences that appear to have undergone dramatic human-specific changes in selective pressures, at least 15% of which have substitution rates significantly above that expected under neutrality. The 1,356 'accelerated CNC' (ANC) sequences are enriched in recent segmental duplications, suggesting a recent change in selective constraint following duplication. In addition, single nucleotide polymorphisms within ANC sequences have a significant excess of high frequency derived alleles and high F(ST)values relative to controls, indicating that acceleration and positive selection are recent in human populations. Finally, a significant number of single nucleotide polymorphisms within ANC sequences are associated with changes in gene expression. The probability of variation in an ANC sequence being associated with a gene expression phenotype is fivefold higher than variation in a control CNC sequence. CONCLUSION: Our analysis suggests that ANC sequences have until very recently played a role in human evolution, potentially through lineage-specific changes in gene regulation

Competing risks of death in women treated with adjuvant aromatase inhibitors for early breast cancer on NCIC CTG MA.27

Baseline patient and tumor characteristics differentially affected type of death in the MA.17 placebo-controlled letrozole trial where cardiovascular death was not separately identified. The MA.27 trial allowed competing risks analysis of breast cancer (BC), cardiovascular, and other type (OT) of death. MA.27 was a phase III adjuvant breast cancer trial of exemestane versus anastrozole. Effects of baseline patient and tumor characteristics were tested for whether factors were associated with (1) all cause mortality and (2) cause-specific mortality. We also fit step-wise forward cause-specific-adjusted models. 7576 women (median age 64 years; 5417 (72 %) < 70 years and 2159 (28 %) ≥ 70 years) were enrolled and followed for median 4.1 years. The 432 deaths comprised 187 (43 %) BC, 66 (15 %) cardiovascular, and 179 (41 %) OT. Five baseline factors were differentially associated with type of death. Older patients had greater BC (p = 0.03), cardiovascular (p < 0.001), and other types (p < 0.001) of mortality. Patients with pre-existing cardiovascular history had worse cardiovascular mortality (p < 0.001); those with worse ECOG performance status had worse OT mortality (p < 0.001). Patients with T1 tumors (p < 0.001) and progesterone receptor positive had less BC mortality (p < 0.001). Fewer BC deaths occurred with node-negative disease (p < 0.001), estrogen receptor-positive tumors (p = 0.001), and without adjuvant chemotherapy (p = 0.005); worse cardiovascular mortality (p = 0.01), with trastuzumab; worse OT mortality, for non-whites (p = 0.03) and without adjuvant radiotherapy (p = 0.003). Overall, 57 % of deaths in MA.27 AI-treated patients were non-breast cancer related. Baseline patient and tumor characteristics differentially affected type of death with women 70 or older experiencing more non-breast cancer death

Modifier Effects between Regulatory and Protein-Coding Variation

Genome-wide associations have shown a lot of promise in dissecting the genetics of complex traits in humans with single variants, yet a large fraction of the genetic effects is still unaccounted for. Analyzing genetic interactions between variants (epistasis) is one of the potential ways forward. We investigated the abundance and functional impact of a specific type of epistasis, namely the interaction between regulatory and protein-coding variants. Using genotype and gene expression data from the 210 unrelated individuals of the original four HapMap populations, we have explored the combined effects of regulatory and protein-coding single nucleotide polymorphisms (SNPs). We predict that about 18% (1,502 out of 8,233 nsSNPs) of protein-coding variants are differentially expressed among individuals and demonstrate that regulatory variants can modify the functional effect of a coding variant in cis. Furthermore, we show that such interactions in cis can affect the expression of downstream targets of the gene containing the protein-coding SNP. In this way, a cis interaction between regulatory and protein-coding variants has a trans impact on gene expression. Given the abundance of both types of variants in human populations, we propose that joint consideration of regulatory and protein-coding variants may reveal additional genetic effects underlying complex traits and disease and may shed light on causes of differential penetrance of known disease variants

Effects of adjuvant exemestane versus anastrozole on bone mineral density for women with early breast cancer (MA.27B): a companion analysis of a randomised controlled trial

Treatment of breast cancer with aromatase inhibitors is associated with damage to bones. NCIC CTG MA.27 was an open-label, phase 3, randomised controlled trial in which women with breast cancer were assigned to one of two adjuvant oral aromatase inhibitors—exemestane or anastrozole. We postulated that exemestane—a mildly androgenic steroid—might have a less detrimental effect on bone than non-steroidal anastrozole. In this companion study to MA.27, we compared changes in bone mineral density (BMD) in the lumbar spine and total hip between patients treated with exemestane and patients treated with anastrozole

Prediction of Late Disease Recurrence and Extended Adjuvant Letrozole Benefit by the HOXB13/IL17BR Biomarker

BackgroundBiomarkers to optimize extended adjuvant endocrine therapy for women with estrogen receptor (ER)–positive breast cancer are limited. The HOXB13/IL17BR (H/I) biomarker predicts recurrence risk in ER-positive, lymph node–negative breast cancer patients. H/I was evaluated in MA.17 trial for prognostic performance for late recurrence and treatment benefit from extended adjuvant letrozole.MethodsA prospective–retrospective, nested case-control design of 83 recurrences matched to 166 nonrecurrences from letrozole- and placebo-treated patients within MA.17 was conducted. Expression of H/I within primary tumors was determined by reverse-transcription polymerase chain reaction with a prespecified cutpoint. The predictive ability of H/I for ascertaining benefit from letrozole was determined using multivariable conditional logistic regression including standard clinicopathological factors as covariates. All statistical tests were two-sided.ResultsHigh H/I was statistically significantly associated with a decrease in late recurrence in patients receiving extended letrozole therapy (odds ratio [OR] = 0.35; 95% confidence interval [CI] = 0.16 to 0.75; P = .007). In an adjusted model with standard clinicopathological factors, high H/I remained statistically significantly associated with patient benefit from letrozole (OR = 0.33; 95% CI = 0.15 to 0.73; P = .006). Reduction in the absolute risk of recurrence at 5 years was 16.5% for patients with high H/I (P = .007). The interaction between H/I and letrozole treatment was statistically significant (P = .03).ConclusionsIn the absence of extended letrozole therapy, high H/I identifies a subgroup of ER-positive patients disease-free after 5 years of tamoxifen who are at risk for late recurrence. When extended endocrine therapy with letrozole is prescribed, high H/I predicts benefit from therapy and a decreased probability of late disease recurrence

Patterns of Cis Regulatory Variation in Diverse Human Populations

The genetic basis of gene expression variation has long been studied with the aim to understand the landscape of regulatory variants, but also more recently to assist in the interpretation and elucidation of disease signals. To date, many studies have looked in specific tissues and population-based samples, but there has been limited assessment of the degree of inter-population variability in regulatory variation. We analyzed genome-wide gene expression in lymphoblastoid cell lines from a total of 726 individuals from 8 global populations from the HapMap3 project and correlated gene expression levels with HapMap3 SNPs located in cis to the genes. We describe the influence of ancestry on gene expression levels within and between these diverse human populations and uncover a non-negligible impact on global patterns of gene expression. We further dissect the specific functional pathways differentiated between populations. We also identify 5,691 expression quantitative trait loci (eQTLs) after controlling for both non-genetic factors and population admixture and observe that half of the cis-eQTLs are replicated in one or more of the populations. We highlight patterns of eQTL-sharing between populations, which are partially determined by population genetic relatedness, and discover significant sharing of eQTL effects between Asians, European-admixed, and African subpopulations. Specifically, we observe that both the effect size and the direction of effect for eQTLs are highly conserved across populations. We observe an increasing proximity of eQTLs toward the transcription start site as sharing of eQTLs among populations increases, highlighting that variants close to TSS have stronger effects and therefore are more likely to be detected across a wider panel of populations. Together these results offer a unique picture and resource of the degree of differentiation among human populations in functional regulatory variation and provide an estimate for the transferability of complex trait variants across populations

Characterization of functional methylomes by next-generation capture sequencing identifies novel disease-associated variants.

Most genome-wide methylation studies (EWAS) of multifactorial disease traits use targeted arrays or enrichment methodologies preferentially covering CpG-dense regions, to characterize sufficiently large samples. To overcome this limitation, we present here a new customizable, cost-effective approach, methylC-capture sequencing (MCC-Seq), for sequencing functional methylomes, while simultaneously providing genetic variation information. To illustrate MCC-Seq, we use whole-genome bisulfite sequencing on adipose tissue (AT) samples and public databases to design AT-specific panels. We establish its efficiency for high-density interrogation of methylome variability by systematic comparisons with other approaches and demonstrate its applicability by identifying novel methylation variation within enhancers strongly correlated to plasma triglyceride and HDL-cholesterol, including at CD36. Our more comprehensive AT panel assesses tissue methylation and genotypes in parallel at ∼4 and ∼3 M sites, respectively. Our study demonstrates that MCC-Seq provides comparable accuracy to alternative approaches but enables more efficient cataloguing of functional and disease-relevant epigenetic and genetic variants for large-scale EWAS.This work was supported by a Canadian Institute of Health Research (CIHR) team grant awarded to E.G., A.T., M.C.V. and M.L. (TEC-128093) and the CIHR funded Epigeneome Mapping Centre at McGill University (EP1-120608) awarded to T.P. and M.L., and the Swedish Research Council, Knut and Alice Wallenberg Foundation and the Torsten Söderberg Foundation awarded to L.R. F.A. holds studentship from The Research Institute of the McGill University Health Center (MUHC). F.G. is a recipient of a research fellowship award from the Heart and Stroke Foundation of Canada. A.T. is the director of a Research Chair in Bariatric and Metabolic Surgery. M.C.V. is the recipient of the Canada Research Chair in Genomics Applied to Nutrition and Health (Tier 1). E.G. and T.P. are recipients of a Canada Research Chair Tier 2 award. The MuTHER Study was funded by a programme grant from the Wellcome Trust (081917/Z/07/Z) and core funding for the Wellcome Trust Centre for Human Genetics (090532). TwinsUK was funded by the Wellcome Trust; European Community's Seventh Framework Programme (FP7/2007-2013). The study also receives support from the National Institute for Health Research (NIHR)-funded BioResource, Clinical Research Facility and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. T.D.S. is a holder of an ERC Advanced Principal Investigator award. SNP genotyping was performed by The Wellcome Trust Sanger Institute and National Eye Institute via NIH/CIDR. Finally, we thank the NIH Roadmap Epigenomics Consortium and the Mapping Centers (http://nihroadmap.nih.gov/epigenomics/) for the production of publicly available reference epigenomes. Specifically, we thank the mapping centre at MGH/BROAD for generation of human adipose reference epigenomes used in this study.This is the final version. It was first published by NPG at http://www.nature.com/ncomms/2015/150529/ncomms8211/full/ncomms8211.html#abstrac

New genetic loci link adipose and insulin biology to body fat distribution.

Body fat distribution is a heritable trait and a well-established predictor of adverse metabolic outcomes, independent of overall adiposity. To increase our understanding of the genetic basis of body fat distribution and its molecular links to cardiometabolic traits, here we conduct genome-wide association meta-analyses of traits related to waist and hip circumferences in up to 224,459 individuals. We identify 49 loci (33 new) associated with waist-to-hip ratio adjusted for body mass index (BMI), and an additional 19 loci newly associated with related waist and hip circumference measures (P < 5 × 10(-8)). In total, 20 of the 49 waist-to-hip ratio adjusted for BMI loci show significant sexual dimorphism, 19 of which display a stronger effect in women. The identified loci were enriched for genes expressed in adipose tissue and for putative regulatory elements in adipocytes. Pathway analyses implicated adipogenesis, angiogenesis, transcriptional regulation and insulin resistance as processes affecting fat distribution, providing insight into potential pathophysiological mechanisms

Genome-wide association identifies nine common variants associated with fasting proinsulin levels and provides new insights into the pathophysiology of type 2 diabetes.

OBJECTIVE: Proinsulin is a precursor of mature insulin and C-peptide. Higher circulating proinsulin levels are associated with impaired β-cell function, raised glucose levels, insulin resistance, and type 2 diabetes (T2D). Studies of the insulin processing pathway could provide new insights about T2D pathophysiology. RESEARCH DESIGN AND METHODS: We have conducted a meta-analysis of genome-wide association tests of ∼2.5 million genotyped or imputed single nucleotide polymorphisms (SNPs) and fasting proinsulin levels in 10,701 nondiabetic adults of European ancestry, with follow-up of 23 loci in up to 16,378 individuals, using additive genetic models adjusted for age, sex, fasting insulin, and study-specific covariates. RESULTS: Nine SNPs at eight loci were associated with proinsulin levels (P < 5 × 10(-8)). Two loci (LARP6 and SGSM2) have not been previously related to metabolic traits, one (MADD) has been associated with fasting glucose, one (PCSK1) has been implicated in obesity, and four (TCF7L2, SLC30A8, VPS13C/C2CD4A/B, and ARAP1, formerly CENTD2) increase T2D risk. The proinsulin-raising allele of ARAP1 was associated with a lower fasting glucose (P = 1.7 × 10(-4)), improved β-cell function (P = 1.1 × 10(-5)), and lower risk of T2D (odds ratio 0.88; P = 7.8 × 10(-6)). Notably, PCSK1 encodes the protein prohormone convertase 1/3, the first enzyme in the insulin processing pathway. A genotype score composed of the nine proinsulin-raising alleles was not associated with coronary disease in two large case-control datasets. CONCLUSIONS: We have identified nine genetic variants associated with fasting proinsulin. Our findings illuminate the biology underlying glucose homeostasis and T2D development in humans and argue against a direct role of proinsulin in coronary artery disease pathogenesis