The HEK293 human cell lineage is widely used in cell biology and biotechnology. Here we use whole-genome resequencing of six 293 cell lines to study the dynamics of this aneuploid genome in response to the manipulations used to generate common 293 cell derivatives, such as transformation and stable clone generation (293T); suspension growth adaptation (293S); and cytotoxic lectin selection (293SG). Remarkably, we observe that copy number alteration detection could identify the genomic region that enabled cell survival under selective conditions (i.c. ricin selection). Furthermore, we present methods to detect human/vector genome breakpoints and a user-friendly visualization tool for the 293 genome data. We also establish that the genome structure composition is in steady state for most of these cell lines when standard cell culturing conditions are used. This resource enables novel and more informed studies with 293 cells, and we will distribute the sequenced cell lines to this effect.

Genome-wide association studies have reported 11 regions conferring risk of high-grade serous epithelial ovarian cancer (HGSOC). Expression quantitative trait locus (eQTL) analyses can identify candidate susceptibility genes at risk loci. Here we evaluate cis-eQTL associations at 47 regions associated with HGSOC risk (P≤10(-5)). For three cis-eQTL associations (P<1.4 × 10(-3), FDR<0.05) at 1p36 (CDC42), 1p34 (CDCA8) and 2q31 (HOXD9), we evaluate the functional role of each candidate by perturbing expression of each gene in HGSOC precursor cells. Overexpression of HOXD9 increases anchorage-independent growth, shortens population-doubling time and reduces contact inhibition. Chromosome conformation capture identifies an interaction between rs2857532 and the HOXD9 promoter, suggesting this SNP is a leading causal variant. Transcriptomic profiling after HOXD9 overexpression reveals enrichment of HGSOC risk variants within HOXD9 target genes (P=6 × 10(-10) for risk variants (P<10(-4)) within 10 kb of a HOXD9 target gene in ovarian cells), suggesting a broader role for this network in genetic susceptibility to HGSOC.

Common variants in 94 loci have been associated with breast cancer including 15 loci with genome-wide significant associations (P<5 × 10(-8)) with oestrogen receptor (ER)-negative breast cancer and BRCA1-associated breast cancer risk. In this study, to identify new ER-negative susceptibility loci, we performed a meta-analysis of 11 genome-wide association studies (GWAS) consisting of 4,939 ER-negative cases and 14,352 controls, combined with 7,333 ER-negative cases and 42,468 controls and 15,252 BRCA1 mutation carriers genotyped on the iCOGS array. We identify four previously unidentified loci including two loci at 13q22 near KLF5, a 2p23.2 locus near WDR43 and a 2q33 locus near PPIL3 that display genome-wide significant associations with ER-negative breast cancer. In addition, 19 known breast cancer risk loci have genome-wide significant associations and 40 had moderate associations (P<0.05) with ER-negative disease. Using functional and eQTL studies we implicate TRMT61B and WDR43 at 2p23.2 and PPIL3 at 2q33 in ER-negative breast cancer aetiology. All ER-negative loci combined account for ∼11% of familial relative risk for ER-negative disease and may contribute to improved ER-negative and BRCA1 breast cancer risk prediction.

A locus at 19p13 is associated with breast cancer (BC) and ovarian cancer (OC) risk. Here we analyse 438 SNPs in this region in 46,451 BC and 15,438 OC cases, 15,252 BRCA1 mutation carriers and 73,444 controls and identify 13 candidate causal SNPs associated with serous OC (P=9.2 × 10(-20)), ER-negative BC (P=1.1 × 10(-13)), BRCA1-associated BC (P=7.7 × 10(-16)) and triple negative BC (P-diff=2 × 10(-5)). Genotype-gene expression associations are identified for candidate target genes ANKLE1 (P=2 × 10(-3)) and ABHD8 (P<2 × 10(-3)). Chromosome conformation capture identifies interactions between four candidate SNPs and ABHD8, and luciferase assays indicate six risk alleles increased transactivation of the ADHD8 promoter. Targeted deletion of a region containing risk SNP rs56069439 in a putative enhancer induces ANKLE1 downregulation; and mRNA stability assays indicate functional effects for an ANKLE1 3'-UTR SNP. Altogether, these data suggest that multiple SNPs at 19p13 regulate ABHD8 and perhaps ANKLE1 expression, and indicate common mechanisms underlying breast and ovarian cancer risk.

Increased copy number alterations (CNAs) indicative of chromosomal instability (CIN) have been associated with poor cancer outcome. Here, we study CNAs as potential biomarkers of bevacizumab (BVZ) response in metastatic colorectal cancer (mCRC). We cluster 409 mCRCs in three subclusters characterized by different degrees of CIN. Tumors belonging to intermediate-to-high instability clusters have improved outcome following chemotherapy plus BVZ versus chemotherapy alone. In contrast, low instability tumors, which amongst others consist of POLE-mutated and microsatellite-instable tumors, derive no further benefit from BVZ. This is confirmed in 81 mCRC tumors from the phase 2 MoMa study involving BVZ. CNA clusters overlap with CRC consensus molecular subtypes (CMS); CMS2/4 xenografts correspond to intermediate-to-high instability clusters and respond to FOLFOX chemotherapy plus mouse avastin (B20), while CMS1/3 xenografts match with low instability clusters and fail to respond. Overall, we identify copy number load as a novel potential predictive biomarker of BVZ combination therapy.

Rejection remains the main cause of premature graft loss after kidney transplantation, despite the use of potent immunosuppression. This highlights the need to better understand the composition and the cell-to-cell interactions of the alloreactive inflammatory infiltrate. Here, we performed droplet-based single-cell RNA sequencing of 35,152 transcriptomes from 16 kidney transplant biopsies with varying phenotypes and severities of rejection and without rejection, and identified cell-type specific gene expression signatures for deconvolution of bulk tissue. A specific association was identified between recipient-derived FCGR3A+ monocytes, FCGR3A+ NK cells and the severity of intragraft inflammation. Activated FCGR3A+ monocytes overexpressed CD47 and LILR genes and increased paracrine signaling pathways promoting T cell infiltration. FCGR3A+ NK cells overexpressed FCRL3, suggesting that antibody-dependent cytotoxicity is a central mechanism of NK-cell mediated graft injury. Multiplexed immunofluorescence using 38 markers on 18 independent biopsy slides confirmed this role of FcγRIII+ NK and FcγRIII+ nonclassical monocytes in antibody-mediated rejection, with specificity to the glomerular area. These results highlight the central involvement of innate immune cells in the pathogenesis of allograft rejection and identify several potential therapeutic targets that might improve allograft longevity.

Identifying molecular alterations in normal tissue adjacent to cancer is important for understanding cancer aetiology and designing preventive measures. Here we analyse the DNA methylome of 569 breast tissue samples, including 50 from cancer-free women and 84 from matched normal cancer pairs. We use statistical algorithms for dissecting intra- and inter-sample cellular heterogeneity and demonstrate that normal tissue adjacent to breast cancer is characterized by tens to thousands of epigenetic alterations. We show that their genomic distribution is non-random, being strongly enriched for binding sites of transcription factors specifying chromatin architecture. We validate the field defects in an independent cohort and demonstrate that over 30% of the alterations exhibit increased enrichment within matched cancer samples. Breast cancers highly enriched for epigenetic field defects, exhibit adverse clinical outcome. Our data support a model where clonal epigenetic reprogramming towards reduced differentiation in normal tissue is an important step in breast carcinogenesis.

HNF1B is overexpressed in clear cell epithelial ovarian cancer, and we observed epigenetic silencing in serous epithelial ovarian cancer, leading us to hypothesize that variation in this gene differentially associates with epithelial ovarian cancer risk according to histological subtype. Here we comprehensively map variation in HNF1B with respect to epithelial ovarian cancer risk and analyse DNA methylation and expression profiles across histological subtypes. Different single-nucleotide polymorphisms associate with invasive serous (rs7405776 odds ratio (OR)=1.13, P=3.1 × 10(-10)) and clear cell (rs11651755 OR=0.77, P=1.6 × 10(-8)) epithelial ovarian cancer. Risk alleles for the serous subtype associate with higher HNF1B-promoter methylation in these tumours. Unmethylated, expressed HNF1B, primarily present in clear cell tumours, coincides with a CpG island methylator phenotype affecting numerous other promoters throughout the genome. Different variants in HNF1B associate with risk of serous and clear cell epithelial ovarian cancer; DNA methylation and expression patterns are also notably distinct between these subtypes. These findings underscore distinct mechanisms driving different epithelial ovarian cancer histological subtypes.

PTEN is a multifaceted tumor suppressor that is highly sensitive to alterations in expression or function. The PTEN C-tail domain, which is rich in phosphorylation sites, has been implicated in PTEN stability, localization, catalytic activity, and protein interactions, but its role in tumorigenesis remains unclear. To address this, we utilized several mouse strains with nonlethal C-tail mutations. Mice homozygous for a deletion that includes S370, S380, T382 and T383 contain low PTEN levels and hyperactive AKT but are not tumor prone. Analysis of mice containing nonphosphorylatable or phosphomimetic versions of S380, a residue hyperphosphorylated in human gastric cancers, reveal that PTEN stability and ability to inhibit PI3K-AKT depends on dynamic phosphorylation-dephosphorylation of this residue. While phosphomimetic S380 drives neoplastic growth in prostate by promoting nuclear accumulation of β-catenin, nonphosphorylatable S380 is not tumorigenic. These data suggest that C-tail hyperphosphorylation creates oncogenic PTEN and is a potential target for anti-cancer therapy.

Resistance to platinum-based chemotherapy represents a major clinical challenge for many tumors, including epithelial ovarian cancer. Patients often experience several response-relapse events, until tumors become resistant and life expectancy drops to 12-15 months. Despite improved knowledge of the molecular determinants of platinum resistance, the lack of clinical applicability limits exploitation of many potential targets, leaving patients with limited options. Serine biosynthesis has been linked to cancer growth and poor prognosis in various cancer types, however its role in platinum-resistant ovarian cancer is not known. Here, we show that a subgroup of resistant tumors decreases phosphoglycerate dehydrogenase (PHGDH) expression at relapse after platinum-based chemotherapy. Mechanistically, we observe that this phenomenon is accompanied by a specific oxidized nicotinamide adenine dinucleotide (NAD+) regenerating phenotype, which helps tumor cells in sustaining Poly (ADP-ribose) polymerase (PARP) activity under platinum treatment. Our findings reveal metabolic vulnerabilities with clinical implications for a subset of platinum resistant ovarian cancers.

Reduced BRCA1 expression causes homologous recombination (HR) repair defects in high-grade serous ovarian cancers (HGSOCs). Here, we demonstrate that BRCA1 is transcriptionally activated by a previously unknown function of ZC3H18. We show that ZC3H18 is a DNA-binding protein that interacts with an E2F site in the BRCA1 promoter where it facilitates recruitment of E2F4 to an adjacent E2F site to promote BRCA1 transcription. Consistent with ZC3H18 role in activating BRCA1 expression, ZC3H18 depletion induces BRCA1 promoter methylation, reduces BRCA1 expression, disrupts HR, and sensitizes cells to DNA crosslinkers and poly(ADP-ribose) polymerase inhibitors. Moreover, in patient-derived xenografts and primary HGSOC tumors, ZC3H18 and E2F4 mRNA levels are positively correlated with BRCA1 mRNA levels, further supporting ZC3H18 role in regulating BRCA1. Given that ZC3H18 lies within 16q24.2, a region with frequent copy number loss in HGSOC, these findings suggest that ZC3H18 copy number losses could contribute to HR defects in HGSOC.

Obesity is associated with an increased risk of developing breast cancer (BC) and worse prognosis in BC patients, yet its impact on BC biology remains understudied in humans. This study investigates how the biology of untreated primary BC differs according to patients' body mass index (BMI) using data from >2,000 patients. We identify several genomic alterations that are differentially prevalent in overweight or obese patients compared to lean patients. We report evidence supporting an ageing accelerating effect of obesity at the genetic level. We show that BMI-associated differences in bulk transcriptomic profile are subtle, while single cell profiling allows detection of more pronounced changes in different cell compartments. These analyses further reveal an elevated and unresolved inflammation of the BC tumor microenvironment associated with obesity, with distinct characteristics contingent on the estrogen receptor status. Collectively, our analyses imply that obesity is associated with an inflammaging-like phenotype. We conclude that patient adiposity may play a significant role in the heterogeneity of BC and should be considered for BC treatment tailoring.

GWAS have identified a breast cancer susceptibility locus on 2q35. Here we report the fine mapping of this locus using data from 101,943 subjects from 50 case-control studies. We genotype 276 SNPs using the 'iCOGS' genotyping array and impute genotypes for a further 1,284 using 1000 Genomes Project data. All but two, strongly correlated SNPs (rs4442975 G/T and rs6721996 G/A) are excluded as candidate causal variants at odds against >100:1. The best functional candidate, rs4442975, is associated with oestrogen receptor positive (ER+) disease with an odds ratio (OR) in Europeans of 0.85 (95% confidence interval=0.84-0.87; P=1.7 × 10(-43)) per t-allele. This SNP flanks a transcriptional enhancer that physically interacts with the promoter of IGFBP5 (encoding insulin-like growth factor-binding protein 5) and displays allele-specific gene expression, FOXA1 binding and chromatin looping. Evidence suggests that the g-allele confers increased breast cancer susceptibility through relative downregulation of IGFBP5, a gene with known roles in breast cell biology.

Mammographic density reflects the amount of stromal and epithelial tissues in relation to adipose tissue in the breast and is a strong risk factor for breast cancer. Here we report the results from meta-analysis of genome-wide association studies (GWAS) of three mammographic density phenotypes: dense area, non-dense area and percent density in up to 7,916 women in stage 1 and an additional 10,379 women in stage 2. We identify genome-wide significant (P<5 × 10(-8)) loci for dense area (AREG, ESR1, ZNF365, LSP1/TNNT3, IGF1, TMEM184B and SGSM3/MKL1), non-dense area (8p11.23) and percent density (PRDM6, 8p11.23 and TMEM184B). Four of these regions are known breast cancer susceptibility loci, and four additional regions were found to be associated with breast cancer (P<0.05) in a large meta-analysis. These results provide further evidence of a shared genetic basis between mammographic density and breast cancer and illustrate the power of studying intermediate quantitative phenotypes to identify putative disease-susceptibility loci.

The Ashkenazi Jewish (AJ) population is a genetic isolate close to European and Middle Eastern groups, with genetic diversity patterns conducive to disease mapping. Here we report high-depth sequencing of 128 complete genomes of AJ controls. Compared with European samples, our AJ panel has 47% more novel variants per genome and is eightfold more effective at filtering benign variants out of AJ clinical genomes. Our panel improves imputation accuracy for AJ SNP arrays by 28%, and covers at least one haplotype in ≈ 67% of any AJ genome with long, identical-by-descent segments. Reconstruction of recent AJ history from such segments confirms a recent bottleneck of merely ≈ 350 individuals. Modelling of ancient histories for AJ and European populations using their joint allele frequency spectrum determines AJ to be an even admixture of European and likely Middle Eastern origins. We date the split between the two ancestral populations to ≈ 12-25 Kyr, suggesting a predominantly Near Eastern source for the repopulation of Europe after the Last Glacial Maximum.

Quantifying the genetic correlation between cancers can provide important insights into the mechanisms driving cancer etiology. Using genome-wide association study summary statistics across six cancer types based on a total of 296,215 cases and 301,319 controls of European ancestry, here we estimate the pair-wise genetic correlations between breast, colorectal, head/neck, lung, ovary and prostate cancer, and between cancers and 38 other diseases. We observed statistically significant genetic correlations between lung and head/neck cancer (rg = 0.57, p = 4.6 × 10-8), breast and ovarian cancer (rg = 0.24, p = 7 × 10-5), breast and lung cancer (rg = 0.18, p =1.5 × 10-6) and breast and colorectal cancer (rg = 0.15, p = 1.1 × 10-4). We also found that multiple cancers are genetically correlated with non-cancer traits including smoking, psychiatric diseases and metabolic characteristics. Functional enrichment analysis revealed a significant excess contribution of conserved and regulatory regions to cancer heritability. Our comprehensive analysis of cross-cancer heritability suggests that solid tumors arising across tissues share in part a common germline genetic basis.

Understanding the pathology of COVID-19 is a global research priority. Early evidence suggests that the respiratory microbiome may be playing a role in disease progression, yet current studies report contradictory results. Here, we examine potential confounders in COVID-19 respiratory microbiome studies by analyzing the upper (n = 58) and lower (n = 35) respiratory tract microbiome in well-phenotyped COVID-19 patients and controls combining microbiome sequencing, viral load determination, and immunoprofiling. We find that time in the intensive care unit and type of oxygen support, as well as associated treatments such as antibiotic usage, explain the most variation within the upper respiratory tract microbiome, while SARS-CoV-2 viral load has a reduced impact. Specifically, mechanical ventilation is linked to altered community structure and significant shifts in oral taxa previously associated with COVID-19. Single-cell transcriptomics of the lower respiratory tract of COVID-19 patients identifies specific oral bacteria in physical association with proinflammatory immune cells, which show higher levels of inflammatory markers. Overall, our findings suggest confounders are driving contradictory results in current COVID-19 microbiome studies and careful attention needs to be paid to ICU stay and type of oxygen support, as bacteria favored in these conditions may contribute to the inflammatory phenotypes observed in severe COVID-19 patients.

Genome-wide association studies (GWAS) have identified more than 170 breast cancer susceptibility loci. Here we hypothesize that some risk-associated variants might act in non-breast tissues, specifically adipose tissue and immune cells from blood and spleen. Using expression quantitative trait loci (eQTL) reported in these tissues, we identify 26 previously unreported, likely target genes of overall breast cancer risk variants, and 17 for estrogen receptor (ER)-negative breast cancer, several with a known immune function. We determine the directional effect of gene expression on disease risk measured based on single and multiple eQTL. In addition, using a gene-based test of association that considers eQTL from multiple tissues, we identify seven (and four) regions with variants associated with overall (and ER-negative) breast cancer risk, which were not reported in previous GWAS. Further investigation of the function of the implicated genes in breast and immune cells may provide insights into the etiology of breast cancer.

Endometrial cancer is the most commonly diagnosed cancer of the female reproductive tract in developed countries. Through genome-wide association studies (GWAS), we have previously identified eight risk loci for endometrial cancer. Here, we present an expanded meta-analysis of 12,906 endometrial cancer cases and 108,979 controls (including new genotype data for 5624 cases) and identify nine novel genome-wide significant loci, including a locus on 12q24.12 previously identified by meta-GWAS of endometrial and colorectal cancer. At five loci, expression quantitative trait locus (eQTL) analyses identify candidate causal genes; risk alleles at two of these loci associate with decreased expression of genes, which encode negative regulators of oncogenic signal transduction proteins (SH2B3 (12q24.12) and NF1 (17q11.2)). In summary, this study has doubled the number of known endometrial cancer risk loci and revealed candidate causal genes for future study.

Epithelial ovarian cancer (EOC) has a heritable component that remains to be fully characterized. Most identified common susceptibility variants lie in non-protein-coding sequences. We hypothesized that variants in the 3' untranslated region at putative microRNA (miRNA)-binding sites represent functional targets that influence EOC susceptibility. Here, we evaluate the association between 767 miRNA-related single-nucleotide polymorphisms (miRSNPs) and EOC risk in 18,174 EOC cases and 26,134 controls from 43 studies genotyped through the Collaborative Oncological Gene-environment Study. We identify several miRSNPs associated with invasive serous EOC risk (odds ratio=1.12, P=10(-8)) mapping to an inversion polymorphism at 17q21.31. Additional genotyping of non-miRSNPs at 17q21.31 reveals stronger signals outside the inversion (P=10(-10)). Variation at 17q21.31 is associated with neurological diseases, and our collaboration is the first to report an association with EOC susceptibility. An integrated molecular analysis in this region provides evidence for ARHGAP27 and PLEKHM1 as candidate EOC susceptibility genes.