The development and progression of cancer, a collection of diseases with complex genetic architectures, is facilitated by the interplay of multiple etiological factors. This complexity challenges the traditional single-platform study design and calls for an integrated approach to data analysis. However, integration of heterogeneous measurements of biological variation is a non-trivial exercise due to the diversity of the human genome and the variety of output data formats and genome coverage obtained from the commonly used molecular platforms. This review article will provide an introduction to integration strategies used for analyzing genetic risk factors for cancer. We critically examine the ability of these strategies to handle the complexity of the human genome and also accommodate information about the biological and functional interactions between the elements that have been measured-making the assessment of disease risk against a composite genomic factor possible. The focus of this review is to provide an overview and introduction to the main strategies and to discuss where there is a need for further development.

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.

Genetic variations, such as single nucleotide polymorphisms (SNPs) in microRNAs (miRNA) or in the miRNA binding sites may affect the miRNA dependent gene expression regulation, which has been implicated in various cancers, including breast cancer, and may alter individual susceptibility to cancer. We investigated associations between miRNA related SNPs and breast cancer risk. First we evaluated 2,196 SNPs in a case-control study combining nine genome wide association studies (GWAS). Second, we further investigated 42 SNPs with suggestive evidence for association using 41,785 cases and 41,880 controls from 41 studies included in the Breast Cancer Association Consortium (BCAC). Combining the GWAS and BCAC data within a meta-analysis, we estimated main effects on breast cancer risk as well as risks for estrogen receptor (ER) and age defined subgroups. Five miRNA binding site SNPs associated significantly with breast cancer risk: rs1045494 (odds ratio (OR) 0.92; 95% confidence interval (CI): 0.88-0.96), rs1052532 (OR 0.97; 95% CI: 0.95-0.99), rs10719 (OR 0.97; 95% CI: 0.94-0.99), rs4687554 (OR 0.97; 95% CI: 0.95-0.99, and rs3134615 (OR 1.03; 95% CI: 1.01-1.05) located in the 3' UTR of CASP8, HDDC3, DROSHA, MUSTN1, and MYCL1, respectively. DROSHA belongs to miRNA machinery genes and has a central role in initial miRNA processing. The remaining genes are involved in different molecular functions, including apoptosis and gene expression regulation. Further studies are warranted to elucidate whether the miRNA binding site SNPs are the causative variants for the observed risk effects.

Stratification of women according to their risk of breast cancer based on polygenic risk scores (PRSs) could improve screening and prevention strategies. Our aim was to develop PRSs, optimized for prediction of estrogen receptor (ER)-specific disease, from the largest available genome-wide association dataset and to empirically validate the PRSs in prospective studies. The development dataset comprised 94,075 case subjects and 75,017 control subjects of European ancestry from 69 studies, divided into training and validation sets. Samples were genotyped using genome-wide arrays, and single-nucleotide polymorphisms (SNPs) were selected by stepwise regression or lasso penalized regression. The best performing PRSs were validated in an independent test set comprising 11,428 case subjects and 18,323 control subjects from 10 prospective studies and 190,040 women from UK Biobank (3,215 incident breast cancers). For the best PRSs (313 SNPs), the odds ratio for overall disease per 1 standard deviation in ten prospective studies was 1.61 (95%CI: 1.57-1.65) with area under receiver-operator curve (AUC) = 0.630 (95%CI: 0.628-0.651). The lifetime risk of overall breast cancer in the top centile of the PRSs was 32.6%. Compared with women in the middle quintile, those in the highest 1% of risk had 4.37- and 2.78-fold risks, and those in the lowest 1% of risk had 0.16- and 0.27-fold risks, of developing ER-positive and ER-negative disease, respectively. Goodness-of-fit tests indicated that this PRS was well calibrated and predicts disease risk accurately in the tails of the distribution. This PRS is a powerful and reliable predictor of breast cancer risk that may improve breast cancer prevention programs.

Many common diseases have a complex genetic basis in which large numbers of genetic variations combine with environmental factors to determine risk. However, quantifying such polygenic effects has been challenging. In order to address these difficulties we developed a global measure of the information content of an individual's genome relative to a reference population, which may be used to assess differences in global genome structure between cases and appropriate controls. Informally this measure, which we call relative genome information (RGI), quantifies the relative "disorder" of an individual's genome. In order to test its ability to predict disease risk we used RGI to compare single-nucleotide polymorphism genotypes from two independent samples of women with early-onset breast cancer with three independent sets of controls. We found that RGI was significantly elevated in both sets of breast cancer cases in comparison with all three sets of controls, with disease risk rising sharply with RGI. Furthermore, these differences are not due to associations with common variants at a small number of disease-associated loci, but rather are due to the combined associations of thousands of markers distributed throughout the genome. Our results indicate that the information content of an individual's genome may be used to measure the risk of a complex disease, and suggest that early-onset breast cancer has a strongly polygenic component.

Limited evidence exists on the link between inflammation and epigenetic aging. We aimed to (a) assess the cross-sectional and prospective associations of 22 inflammation-related plasma markers and a signature of inflammaging with epigenetic aging and (b) determine whether epigenetic aging and inflammaging are independently associated with mortality. Blood samples from 940 participants in the Melbourne Collaborative Cohort Study collected at baseline (1990-1994) and follow-up (2003-2007) were assayed for DNA methylation and 22 inflammation-related markers, including well-established markers (eg, interleukins and C-reactive protein) and metabolites of the tryptophan-kynurenine pathway. Four measures of epigenetic aging (PhenoAge, GrimAge, DunedinPoAm, and Zhang) and a signature of inflammaging were considered, adjusted for age, and transformed to Z scores. Associations were assessed using linear regression, and mortality hazard ratios (HR) and 95% confidence intervals (95% CI) were estimated using Cox regression. Cross-sectionally, most inflammation-related markers were associated with epigenetic aging measures, although with generally modest effect sizes (regression coefficients per SD ≤ 0.26) and explaining altogether between 1% and 11% of their variation. Prospectively, baseline inflammation-related markers were not, or only weakly, associated with epigenetic aging after 11 years of follow-up. Epigenetic aging and inflammaging were strongly and independently associated with mortality, for example, inflammaging: HR = 1.41, 95% CI = 1.27-1.56, p = 2 × 10-10, which was only slightly attenuated after adjustment for 4 epigenetic aging measures: HR = 1.35, 95% CI = 1.22-1.51, p = 7 × 10-9). Although cross-sectionally associated with epigenetic aging, inflammation-related markers accounted for a modest proportion of its variation. Inflammaging and epigenetic aging are essentially nonoverlapping markers of biological aging and may be used jointly to predict mortality.

We examined the associations between germline variants and breast cancer mortality using a large meta-analysis of women of European ancestry.

Genome-wide association studies (GWAS) and large-scale replication studies have identified common variants in 79 loci associated with breast cancer, explaining ∼14% of the familial risk of the disease. To identify new susceptibility loci, we performed a meta-analysis of 11 GWAS, comprising 15,748 breast cancer cases and 18,084 controls together with 46,785 cases and 42,892 controls from 41 studies genotyped on a 211,155-marker custom array (iCOGS). Analyses were restricted to women of European ancestry. We generated genotypes for more than 11 million SNPs by imputation using the 1000 Genomes Project reference panel, and we identified 15 new loci associated with breast cancer at P < 5 × 10(-8). Combining association analysis with ChIP-seq chromatin binding data in mammary cell lines and ChIA-PET chromatin interaction data from ENCODE, we identified likely target genes in two regions: SETBP1 at 18q12.3 and RNF115 and PDZK1 at 1q21.1. One association appears to be driven by an amino acid substitution encoded in EXO1.

Breast cancer is the most common cancer among women. To date, 22 common breast cancer susceptibility loci have been identified accounting for ∼8% of the heritability of the disease. We attempted to replicate 72 promising associations from two independent genome-wide association studies (GWAS) in ∼70,000 cases and ∼68,000 controls from 41 case-control studies and 9 breast cancer GWAS. We identified three new breast cancer risk loci at 12p11 (rs10771399; P = 2.7 × 10(-35)), 12q24 (rs1292011; P = 4.3 × 10(-19)) and 21q21 (rs2823093; P = 1.1 × 10(-12)). rs10771399 was associated with similar relative risks for both estrogen receptor (ER)-negative and ER-positive breast cancer, whereas the other two loci were associated only with ER-positive disease. Two of the loci lie in regions that contain strong plausible candidate genes: PTHLH (12p11) has a crucial role in mammary gland development and the establishment of bone metastasis in breast cancer, and NRIP1 (21q21) encodes an ER cofactor and has a role in the regulation of breast cancer cell growth.

Lipid dyshomeostasis is associated with the most common form of dementia, Alzheimer's disease (AD). Substantial progress has been made in identifying positron emission tomography and cerebrospinal fluid biomarkers for AD, but they have limited use as front-line diagnostic tools. Extracellular vesicles (EVs) are released by all cells and contain a subset of their parental cell composition, including lipids. EVs are released from the brain into the periphery, providing a potential source of tissue and disease specific lipid biomarkers. However, the EV lipidome of the central nervous system is currently unknown and the potential of brain-derived EVs (BDEVs) to inform on lipid dyshomeostasis in AD remains unclear. The aim of this study was to reveal the lipid composition of BDEVs in human frontal cortex, and to determine whether BDEVs have an altered lipid profile in AD. Using semi-quantitative mass spectrometry, we describe the BDEV lipidome, covering four lipid categories, 17 lipid classes and 692 lipid molecules. BDEVs were enriched in glycerophosphoserine (PS) lipids, a characteristic of small EVs. Here we further report that BDEVs are enriched in ether-containing PS lipids, a finding that further establishes ether lipids as a feature of EVs. BDEVs in the AD frontal cortex offered improved detection of dysregulated lipids in AD over global lipid profiling of this brain region.  AD BDEVs had significantly altered glycerophospholipid and sphingolipid levels, specifically increased plasmalogen glycerophosphoethanolamine and decreased polyunsaturated fatty acyl containing lipids, and altered amide-linked acyl chain content in sphingomyelin and ceramide lipids relative to CTL. The most prominent alteration was a two-fold decrease in lipid species containing anti-inflammatory/pro-resolving docosahexaenoic acid. The in-depth lipidome analysis provided in this study highlights the advantage of EVs over more complex tissues for improved detection of dysregulated lipids that may serve as potential biomarkers in the periphery.

Previous transcriptome-wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome-wide association studies (GWAS), but analyses of breast cancer subtype-specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta-analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER-). We further compared associations with ER+ and ER- subtypes, using a case-only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER- breast cancer. We further identified 30 TWAS-significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast-cancer gene in three of six regions harboring multiple TWAS-significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER- breast cancer.

Glioma accounts for approximately 80% of malignant adult brain cancer and its most common subtype, glioblastoma, has one of the lowest 5-year cancer survivals. Fifty risk-associated variants within 34 glioma genetic risk regions have been found by genome-wide association studies (GWAS) with a sex difference reported for 8q24.21 region. We conducted an Australian GWAS by glioma subtype and sex.

In Australia, there is a high burden of acute rheumatic fever (ARF) among Aboriginal and Torres Strait Islander peoples. Clinical diagnostic criteria can result in a diagnosis of 'definite', 'probable' or 'possible' ARF and outcomes range from recovery to severe rheumatic heart disease (RHD). We compared outcomes by ARF diagnosis, where the main outcome was defined as disease progression from: possible to probable ARF, definite ARF or RHD; probable to definite ARF or RHD; or definite ARF to definite ARF recurrence or RHD. Data were extracted from the Northern Territory RHD register for Indigenous Australians with an initial diagnosis of ARF during the 5.5-year study period (01/01/2013-30/06/2019). Descriptive statistics were used to describe cohort characteristics, probability of survival, and cumulative incidence risk of disease progression. Cox proportional hazards regression was used to determine whether time to disease progression differed according to ARF diagnosis. Sub-analyses on RHD outcome, clinical manifestations, and antibiotic adherence were also performed. In total there were 913 cases with an initial ARF diagnosis. Of these, 92 (13%) experienced disease progression. The probability of disease progression significantly differed between ARF diagnoses (p = 0.0043; log rank test). Cumulative incidence risk of disease progression at 5.5 years was 33.6% (95% CI 23.6-46.2) for definite, 13.5% (95% CI 8.8-20.6) for probable and 11.4% (95% CI 6.0-21.3) for possible ARF. Disease progression was 2.19 times more likely in those with definite ARF than those with possible ARF (p = 0.026). Progression to RHD was reported in 52/732 (7%) of ARF cases with normal baseline echocardiography. There was a significantly higher risk of progression from no RHD to RHD if the initial diagnosis was definite compared to possible ARF (p<0.001). These data provide a useful way to stratify risk and guide prognosis for people diagnosed with ARF and can help inform practice.

Genome-wide association studies (GWAS) are a common approach for systematic discovery of single nucleotide polymorphisms (SNPs) which are associated with a given disease. Univariate analysis approaches commonly employed may miss important SNP associations that only appear through multivariate analysis in complex diseases. However, multivariate SNP analysis is currently limited by its inherent computational complexity. In this work, we present a computational framework that harnesses supercomputers. Based on our results, we estimate a three-way interaction analysis on 1.1 million SNP GWAS data requiring over 5.8 years on the full "Avoca" IBM Blue Gene/Q installation at the Victorian Life Sciences Computation Initiative. This is hundreds of times faster than estimates for other CPU based methods and four times faster than runtimes estimated for GPU methods, indicating how the improvement in the level of hardware applied to interaction analysis may alter the types of analysis that can be performed. Furthermore, the same analysis would take under 3 months on the currently largest IBM Blue Gene/Q supercomputer "Sequoia" at the Lawrence Livermore National Laboratory assuming linear scaling is maintained as our results suggest. Given that the implementation used in this study can be further optimised, this runtime means it is becoming feasible to carry out exhaustive analysis of higher order interaction studies on large modern GWAS.

Various common genetic susceptibility loci have been identified for breast cancer; however, it is unclear how they combine with lifestyle/environmental risk factors to influence risk. We undertook an international collaborative study to assess gene-environment interaction for risk of breast cancer. Data from 24 studies of the Breast Cancer Association Consortium were pooled. Using up to 34,793 invasive breast cancers and 41,099 controls, we examined whether the relative risks associated with 23 single nucleotide polymorphisms were modified by 10 established environmental risk factors (age at menarche, parity, breastfeeding, body mass index, height, oral contraceptive use, menopausal hormone therapy use, alcohol consumption, cigarette smoking, physical activity) in women of European ancestry. We used logistic regression models stratified by study and adjusted for age and performed likelihood ratio tests to assess gene-environment interactions. All statistical tests were two-sided. We replicated previously reported potential interactions between LSP1-rs3817198 and parity (Pinteraction = 2.4 × 10(-6)) and between CASP8-rs17468277 and alcohol consumption (Pinteraction = 3.1 × 10(-4)). Overall, the per-allele odds ratio (95% confidence interval) for LSP1-rs3817198 was 1.08 (1.01-1.16) in nulliparous women and ranged from 1.03 (0.96-1.10) in parous women with one birth to 1.26 (1.16-1.37) in women with at least four births. For CASP8-rs17468277, the per-allele OR was 0.91 (0.85-0.98) in those with an alcohol intake of <20 g/day and 1.45 (1.14-1.85) in those who drank ≥ 20 g/day. Additionally, interaction was found between 1p11.2-rs11249433 and ever being parous (Pinteraction = 5.3 × 10(-5)), with a per-allele OR of 1.14 (1.11-1.17) in parous women and 0.98 (0.92-1.05) in nulliparous women. These data provide first strong evidence that the risk of breast cancer associated with some common genetic variants may vary with environmental risk factors.

Candidate variant association studies have been largely unsuccessful in identifying common breast cancer susceptibility variants, although most studies have been underpowered to detect associations of a realistic magnitude. We assessed 41 common non-synonymous single-nucleotide polymorphisms (nsSNPs) for which evidence of association with breast cancer risk had been previously reported. Case-control data were combined from 38 studies of white European women (46 450 cases and 42 600 controls) and analyzed using unconditional logistic regression. Strong evidence of association was observed for three nsSNPs: ATXN7-K264R at 3p21 [rs1053338, per allele OR = 1.07, 95% confidence interval (CI) = 1.04-1.10, P = 2.9 × 10(-6)], AKAP9-M463I at 7q21 (rs6964587, OR = 1.05, 95% CI = 1.03-1.07, P = 1.7 × 10(-6)) and NEK10-L513S at 3p24 (rs10510592, OR = 1.10, 95% CI = 1.07-1.12, P = 5.1 × 10(-17)). The first two associations reached genome-wide statistical significance in a combined analysis of available data, including independent data from nine genome-wide association studies (GWASs): for ATXN7-K264R, OR = 1.07 (95% CI = 1.05-1.10, P = 1.0 × 10(-8)); for AKAP9-M463I, OR = 1.05 (95% CI = 1.04-1.07, P = 2.0 × 10(-10)). Further analysis of other common variants in these two regions suggested that intronic SNPs nearby are more strongly associated with disease risk. We have thus identified a novel susceptibility locus at 3p21, and confirmed previous suggestive evidence that rs6964587 at 7q21 is associated with risk. The third locus, rs10510592, is located in an established breast cancer susceptibility region; the association was substantially attenuated after adjustment for the known GWAS hit. Thus, each of the associated nsSNPs is likely to be a marker for another, non-coding, variant causally related to breast cancer risk. Further fine-mapping and functional studies are required to identify the underlying risk-modifying variants and the genes through which they act.

Breast cancer risk is influenced by rare coding variants in susceptibility genes, such as BRCA1, and many common, mostly non-coding variants. However, much of the genetic contribution to breast cancer risk remains unknown. Here we report the results of a genome-wide association study of breast cancer in 122,977 cases and 105,974 controls of European ancestry and 14,068 cases and 13,104 controls of East Asian ancestry. We identified 65 new loci that are associated with overall breast cancer risk at P < 5 × 10-8. The majority of credible risk single-nucleotide polymorphisms in these loci fall in distal regulatory elements, and by integrating in silico data to predict target genes in breast cells at each locus, we demonstrate a strong overlap between candidate target genes and somatic driver genes in breast tumours. We also find that heritability of breast cancer due to all single-nucleotide polymorphisms in regulatory features was 2-5-fold enriched relative to the genome-wide average, with strong enrichment for particular transcription factor binding sites. These results provide further insight into genetic susceptibility to breast cancer and will improve the use of genetic risk scores for individualized screening and prevention.

Genome-wide association studies have identified breast cancer risk variants in over 150 genomic regions, but the mechanisms underlying risk remain largely unknown. These regions were explored by combining association analysis with in silico genomic feature annotations. We defined 205 independent risk-associated signals with the set of credible causal variants in each one. In parallel, we used a Bayesian approach (PAINTOR) that combines genetic association, linkage disequilibrium and enriched genomic features to determine variants with high posterior probabilities of being causal. Potentially causal variants were significantly over-represented in active gene regulatory regions and transcription factor binding sites. We applied our INQUSIT pipeline for prioritizing genes as targets of those potentially causal variants, using gene expression (expression quantitative trait loci), chromatin interaction and functional annotations. Known cancer drivers, transcription factors and genes in the developmental, apoptosis, immune system and DNA integrity checkpoint gene ontology pathways were over-represented among the highest-confidence target genes.

Clustering of breast and colorectal cancer has been observed within some families and cannot be explained by chance or known high-risk mutations in major susceptibility genes. Potential shared genetic susceptibility between breast and colorectal cancer, not explained by high-penetrance genes, has been postulated. We hypothesized that yet undiscovered genetic variants predispose to a breast-colorectal cancer phenotype.

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.