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

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

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

Long-term Mortality in HIV-Positive Individuals Virally Suppressed for >3 Years With Incomplete CD4 Recovery

Virally suppressed HIV-positive individuals on combination antiretroviral therapy who do not achieve a CD4 count >200 cells/µL have substantially increased long-term mortality. The increased mortality was seen across different patient groups and for all causes of deat

The Architecture of Gene Regulatory Variation across Multiple Human Tissues: The MuTHER Study

While there have been studies exploring regulatory variation in one or more tissues, the complexity of tissue-specificity in multiple primary tissues is not yet well understood. We explore in depth the role of cis-regulatory variation in three human tissues: lymphoblastoid cell lines (LCL), skin, and fat. The samples (156 LCL, 160 skin, 166 fat) were derived simultaneously from a subset of well-phenotyped healthy female twins of the MuTHER resource. We discover an abundance of cis-eQTLs in each tissue similar to previous estimates (858 or 4.7% of genes). In addition, we apply factor analysis (FA) to remove effects of latent variables, thus more than doubling the number of our discoveries (1,822 eQTL genes). The unique study design (Matched Co-Twin Analysis—MCTA) permits immediate replication of eQTLs using co-twins (93%–98%) and validation of the considerable gain in eQTL discovery after FA correction. We highlight the challenges of comparing eQTLs between tissues. After verifying previous significance threshold-based estimates of tissue-specificity, we show their limitations given their dependency on statistical power. We propose that continuous estimates of the proportion of tissue-shared signals and direct comparison of the magnitude of effect on the fold change in expression are essential properties that jointly provide a biologically realistic view of tissue-specificity. Under this framework we demonstrate that 30% of eQTLs are shared among the three tissues studied, while another 29% appear exclusively tissue-specific. However, even among the shared eQTLs, a substantial proportion (10%–20%) have significant differences in the magnitude of fold change between genotypic classes across tissues. Our results underline the need to account for the complexity of eQTL tissue-specificity in an effort to assess consequences of such variants for complex traits