Epistasis for Growth Rate and Total Metabolic Flux in Yeast

Studies of interactions between gene deletions repeatedly show that the effect of epistasis on the growth of yeast cells is roughly null or barely positive. These observations relate generally to the pace of growth, its costs in terms of required metabolites and energy are unknown. We measured the maximum rate at which yeast cultures grow and amounts of glucose they consume per synthesized biomass for strains with none, single, or double gene deletions. Because all strains were maintained under a fermentative mode of growth and thus shared a common pattern of metabolic processes, we used the rate of glucose uptake as a proxy for the total flux of metabolites and energy. In the tested sample, the double deletions showed null or slightly positive epistasis both for the mean growth and mean flux. This concordance is explained by the fact that average efficiency of converting glucose into biomass was nearly constant, that is, it did not change with the strength of growth effect. Individual changes in the efficiency caused by gene deletions did have a genetic basis as they were consistent over several environments and transmitted between single and double deletion strains indicating that the efficiency of growth, although independent of its rate, was appreciably heritable. Together, our results suggest that data on the rate of growth can be used as a proxy for the rate of total metabolism when the goal is to find strong individual interactions or estimate the mean epistatic effect. However, it may be necessary to assay both growth and flux in order to detect smaller individual effects of epistasis

Treponema denticola chymotrypsin-like proteinase may contribute to orodigestive carcinogenesis through immunomodulation

Background: Periodontal pathogens have been linked to oral and gastrointestinal (orodigestive) carcinogenesis. However, the exact mechanisms remain unknown. Treponema denticola (Td) is associated with severe periodontitis, a chronic inflammatory disease leading to tooth loss. The anaerobic spirochete Td is an invasive bacteria due to its major virulence factor chymotrypsin-like proteinase. Here we aimed to investigate the presence of Td chymotrypsin-like proteinase (Td-CTLP) in major orodigestive tumours and to elucidate potential mechanisms for Td to contribute to carcinogenesis. Methods: The presence of Td-CTLP within orodigestive tumour tissues was examined using immunohistochemistry. Oral, tonsillar, and oesophageal squamous cell carcinomas, alongside gastric, pancreatic, and colon adenocarcinomas were stained with a Td-CTLP-specific antibody. Gingival tissue from periodontitis patients served as positive controls. SDS-PAGE and immunoblot were used to analyse the immumodulatory activity of Td-CTLP in vitro. Results: Td-CTLP was present in majority of orodigestive tumour samples. Td-CTLP was found to convert pro MMP-8 and -9 into their active forms. In addition, Td-CTLP was able to degrade the proteinase inhibitors TIMP-1, TIMP-2, and alpha-1-antichymotrypsin, as well as complement C1q. Conclusions: Because of its presence within tumours and regulatory activity on proteins critical for the regulation of tumour microenvironment and inflammation, the Td-CTLP may contribute to orodigestive carcinogenesis.Peer reviewe

Inheritance of deleterious mutations at both BRCA1 and BRCA2 in an international sample of 32,295 women

Background: Most $\textit{BRCA1}$ or $\textit{BRCA2}$ mutation carriers have inherited a single (heterozygous) mutation. Transheterozygotes (TH) who have inherited deleterious mutations in both $\textit{BRCA1}$ and $\textit{BRCA2}$ are rare, and the consequences of transheterozygosity are poorly understood. Methods: From 32,295 female $\textit{BRCA1/2}$ mutation carriers, we identified 93 TH (0.3 %). "Cases" were defined as TH, and "controls" were single mutations at $\textit{BRCA1}$ (SH1) or $\textit{BRCA2}$ (SH2). Matched SH1 "controls" carried a BRCA1 mutation found in the TH "case". Matched SH2 "controls" carried a BRCA2 mutation found in the TH "case". After matching the TH carriers with SH1 or SH2, 91 TH were matched to 9316 SH1, and 89 TH were matched to 3370 SH2. Results: The majority of TH (45.2 %) involved the three common Jewish mutations. TH were more likely than SH1 and SH2 women to have been ever diagnosed with breast cancer (BC; $p$ = 0.002). TH were more likely to be diagnosed with ovarian cancer (OC) than SH2 ($p$ = 0.017), but not SH1. Age at BC diagnosis was the same in TH vs. SH1 ($p$ = 0.231), but was on average 4.5 years younger in TH than in SH2 ($p$ < 0.001). BC in TH was more likely to be estrogen receptor (ER) positive ($p$ = 0.010) or progesterone receptor (PR) positive ($p$ = 0.013) than in SH1, but less likely to be ER positive ($p$ < 0.001) or PR positive ($p$ = 0.012) than SH2. Among 15 tumors from TH patients, there was no clear pattern of loss of heterozygosity (LOH) for $\textit{BRCA1}$ or $\textit{BRCA2}$ in either BC or OC. Conclusions: Our observations suggest that clinical TH phenotypes resemble SH1. However, TH breast tumor marker characteristics are phenotypically intermediate to SH1 and SH2.ACA and the CIMBA data management are funded by Cancer Research UK (C12292/A20861 and C12292/A11174). TRR was supported by R01-CA083855, R01-CA102776, and P50-CA083638. KLN, TMF, and SMD are supported by the Basser Research Center at the University of Pennsylvania. BP is supported by R01-CA112520. Cancer Research UK provided financial support for this work. ACA is a Senior Cancer Research UK Cancer Research Fellow. DFE is Cancer Research UK Principal Research Fellow. Tumor analysis was funded by STOP CANCER (to SJR). Study-specific acknowledgements are as provided in the manuscript

An intergenic risk locus containing an enhancer deletion in 2q35 modulates breast cancer risk by deregulating IGFBP5 expression.

Breast cancer is the most diagnosed malignancy and the second leading cause of cancer mortality in females. Previous association studies have identified variants on 2q35 associated with the risk of breast cancer. To identify functional susceptibility loci for breast cancer, we interrogated the 2q35 gene desert for chromatin architecture and functional variation correlated with gene expression. We report a novel intergenic breast cancer risk locus containing an enhancer copy number variation (enCNV; deletion) located approximately 400Kb upstream to IGFBP5, which overlaps an intergenic ERα-bound enhancer that loops to the IGFBP5 promoter. The enCNV is correlated with modified ERα binding and monoallelic-repression of IGFBP5 following estrogen treatment. We investigated the association of enCNV genotype with breast cancer in 1,182 cases and 1,362 controls, and replicate our findings in an independent set of 62,533 cases and 60,966 controls from 41 case control studies and 11 GWAS. We report a dose-dependent inverse association of 2q35 enCNV genotype (percopy OR=0.68 95%CI 0.55-0.83, P=0.0002; replication OR=0.77 95%CI 0.73-0.82, P=2.1x10(-19)) and identify 13 additional linked variants (r(2)>0.8) in the 20Kb linkage block containing the enCNV (P=3.2x10(-15) - 5.6x10(-17)). These associations were independent of previously reported 2q35 variants, rs13387042/rs4442975 and rs16857609, and were stronger for ER-positive than ER-negative disease. Together, these results suggest that 2q35 breast cancer risk loci may be mediating their effect through IGFBP5

An intergenic risk locus containing an enhancer deletion in 2q35 modulates breast cancer risk by deregulating IGFBP5 expression.

Breast cancer is the most diagnosed malignancy and the second leading cause of cancer mortality in females. Previous association studies have identified variants on 2q35 associated with the risk of breast cancer. To identify functional susceptibility loci for breast cancer, we interrogated the 2q35 gene desert for chromatin architecture and functional variation correlated with gene expression. We report a novel intergenic breast cancer risk locus containing an enhancer copy number variation (enCNV; deletion) located approximately 400Kb upstream to IGFBP5, which overlaps an intergenic ERα-bound enhancer that loops to the IGFBP5 promoter. The enCNV is correlated with modified ERα binding and monoallelic-repression of IGFBP5 following estrogen treatment. We investigated the association of enCNV genotype with breast cancer in 1,182 cases and 1,362 controls, and replicate our findings in an independent set of 62,533 cases and 60,966 controls from 41 case control studies and 11 GWAS. We report a dose-dependent inverse association of 2q35 enCNV genotype (percopy OR=0.68 95%CI 0.55-0.83, P=0.0002; replication OR=0.77 95%CI 0.73-0.82, P=2.1x10(-19)) and identify 13 additional linked variants (r(2)>0.8) in the 20Kb linkage block containing the enCNV (P=3.2x10(-15) - 5.6x10(-17)). These associations were independent of previously reported 2q35 variants, rs13387042/rs4442975 and rs16857609, and were stronger for ER-positive than ER-negative disease. Together, these results suggest that 2q35 breast cancer risk loci may be mediating their effect through IGFBP5

A transcriptome-wide association study among 97,898 women to identify candidate susceptibility genes for epithelial ovarian cancer risk

Large-scale genome-wide association studies (GWAS) have identified approximately 35 loci associated with epithelial ovarian cancer (EOC) risk. The majority of GWAS-identified disease susceptibility variants are located in non-coding regions, and causal genes underlying these associations remain largely unknown. Here we performed a transcriptome-wide association study to search for novel genetic loci and plausible causal genes at known GWAS loci. We used RNA sequencing data (68 normal ovarian-tissue samples from 68 individuals and 6,124 cross-tissue samples from 369 individuals) and high-density genotyping data from European descendants of the Genotype-Tissue Expression (GTEx V6) project to build ovarian and cross-tissue models of genetically regulated expression using elastic net methods. We evaluated 17,121 genes for their cis-predicted gene expression in relation to EOC risk using summary statistics data from GWAS of 97,898 women, including 29,396 EOC cases. With a Bonferroni-corrected significance level of P<2.2×10-6, we identified 35 genes including FZD4 at 11q14.2 (Z=5.08, P=3.83×10-7, the cross-tissue model; 1 Mb away from any GWAS-identified EOC risk variant), a potential novel locus for EOC risk. All other 34 significantly-associated genes were located within 1 Mb of known GWAS-identified loci, including 23 genes at 6 loci not previously linked to EOC risk. Upon conditioning on nearby known EOC GWAS-identified variants, the associations for 31 genes disappeared and 3 genes remained (P<1.47 x 10-3). These data identify one novel locus (FZD4) and 34 genes at 13 known EOC risk loci associated with EOC risk, providing new insights into EOC carcinogenesis

Association of breast cancer risk in BRCA1 and BRCA2 mutation carriers with genetic variants showing differential allelic expression:Identification of a modifier of breast cancer risk at locus 11q22.3

Cis-acting regulatory SNPs resulting in differential allelic expression (DAE) may, in part, explain the underlying phenotypic variation associated with many complex diseases. To investigate whether common variants associated with DAE were involved in breast cancer susceptibility among BRCA1 and BRCA2 mutation carriers, a list of 175 genes was developed based of their involvement in cancer-related pathways.Using data from a genome-wide map of SNPs associated with allelic expression, we assessed the association of similar to 320 SNPs located in the vicinity of these genes with breast and ovarian cancer risks in 15,252 BRCA1 and 8211 BRCA2 mutation carriers ascertained from 54 studies participating in the Consortium of Investigators of Modifiers of BRCA1/2.We identified a region on 11q22.3 that is significantly associated with breast cancer risk in BRCA1 mutation carriers (most significant SNP rs228595 p = 7 x 10(-6)). This association was absent in BRCA2 carriers (p = 0.57). The 11q22.3 region notably encompasses genes such as ACAT1, NPAT, and ATM. Expression quantitative trait loci associations were observed in both normal breast and tumors across this region, namely for ACAT1, ATM, and other genes. In silico analysis revealed some overlap between top risk-associated SNPs and relevant biological features in mammary cell data, which suggests potential functional significance.We identified 11q22.3 as a new modifier locus in BRCA1 carriers. Replication in larger studies using estrogen receptor (ER)-negative or triple-negative (i.e., ER-, progesterone receptor-, and HER2-negative) cases could therefore be helpful to confirm the association of this locus with breast cancer risk.</p

Evaluation of copy-number variants as modifiers of breast and ovarian cancer risk for BRCA1 pathogenic variant carriers.

Genome-wide studies of patients carrying pathogenic variants (mutations) in BRCA1 or BRCA2 have reported strong associations between single-nucleotide polymorphisms (SNPs) and cancer risk. To conduct the first genome-wide association analysis of copy-number variants (CNVs) with breast or ovarian cancer risk in a cohort of 2500 BRCA1 pathogenic variant carriers, CNV discovery was performed using multiple calling algorithms and Illumina 610k SNP array data from a previously published genome-wide association study. Our analysis, which focused on functionally disruptive genomic deletions overlapping gene regions, identified a number of loci associated with risk of breast or ovarian cancer for BRCA1 pathogenic variant carriers. Despite only including putative deletions called by at least two or more algorithms, detection of selected CNVs by ancillary molecular technologies only confirmed 40% of predicted common (>1% allele frequency) variants. These include four loci that were associated (unadjusted P<0.05) with breast cancer (GTF2H2, ZNF385B, NAALADL2 and PSG5), and two loci associated with ovarian cancer (CYP2A7 and OR2A1). An interesting finding from this study was an association of a validated CNV deletion at the CYP2A7 locus (19q13.2) with decreased ovarian cancer risk (relative risk=0.50, P=0.007). Genomic analysis found this deletion coincides with a region displaying strong regulatory potential in ovarian tissue, but not in breast epithelial cells. This study highlighted the need to verify CNVs in vitro, but also provides evidence that experimentally validated CNVs (with plausible biological consequences) can modify risk of breast or ovarian cancer in BRCA1 pathogenic variant carriers.European Journal of Human Genetics advance online publication, 1 February 2017; doi:10.1038/ejhg.2016.203.Multiple funders listed in the article

Breast cancer risk variants at 6q25 display different phenotype associations and regulate ESR1, RMND1 and CCDC170.

We analyzed 3,872 common genetic variants across the ESR1 locus (encoding estrogen receptor α) in 118,816 subjects from three international consortia. We found evidence for at least five independent causal variants, each associated with different phenotype sets, including estrogen receptor (ER(+) or ER(-)) and human ERBB2 (HER2(+) or HER2(-)) tumor subtypes, mammographic density and tumor grade. The best candidate causal variants for ER(-) tumors lie in four separate enhancer elements, and their risk alleles reduce expression of ESR1, RMND1 and CCDC170, whereas the risk alleles of the strongest candidates for the remaining independent causal variant disrupt a silencer element and putatively increase ESR1 and RMND1 expression.This is the author accepted manuscript. The final version is available from Nature Publishing Group via http://dx.doi.org/10.1038/ng.352