Proteogenomic characterization of primary colorectal cancer and metastatic progression identifies proteome-based subtypes and signatures.

Metastatic progression of colorectal adenocarcinoma (CRC) remains poorly understood and poses significant challenges for treatment. To overcome these challenges, we performed multiomics analyses of primary CRC and liver metastases. Genomic alterations, such as structural variants or copy number alterations, were enriched in oncogenes and tumor suppressor genes and increased in metastases. Unsupervised mass spectrometry-based proteomics of 135 primary and 123 metastatic CRCs uncovered distinct proteomic subtypes, three each for primary and metastatic CRCs, respectively. Integrated analyses revealed that hypoxia, stemness, and immune signatures characterize these 6 subtypes. Hypoxic CRC harbors high epithelial-to-mesenchymal transition features and metabolic adaptation. CRC with a stemness signature shows high oncogenic pathway activation and alternative telomere lengthening (ALT) phenotype, especially in metastatic lesions. Tumor microenvironment analysis shows immune evasion via modulation of major histocompatibility complex (MHC) class I/II and antigen processing pathways. This study characterizes both primary and metastatic CRCs and provides a large proteogenomics dataset of metastatic progression.

Quantum processor-inspired machine learning in the biomedical sciences.

Recent advances in high-throughput genomic technologies coupled with exponential increases in computer processing and memory have allowed us to interrogate the complex molecular underpinnings of human disease from a genome-wide perspective. While the deluge of genomic information is expected to increase, a bottleneck in conventional high-performance computing is rapidly approaching. Inspired by recent advances in physical quantum processors, we evaluated several unconventional machine-learning (ML) strategies on actual human tumor data, namely "Ising-type" methods, whose objective function is formulated identical to simulated annealing and quantum annealing. We show the efficacy of multiple Ising-type ML algorithms for classification of multi-omics human cancer data from The Cancer Genome Atlas, comparing these classifiers to a variety of standard ML methods. Our results indicate that Ising-type ML offers superior classification performance with smaller training datasets, thus providing compelling empirical evidence for the potential future application of unconventional computing approaches in the biomedical sciences.

Chronic mTOR activation induces a degradative smooth muscle cell phenotype.

Smooth muscle cell (SMC) proliferation has been thought to limit the progression of thoracic aortic aneurysm and dissection (TAAD) because loss of medial cells associates with advanced disease. We investigated effects of SMC proliferation in the aortic media by conditional disruption of Tsc1, which hyperactivates mTOR complex 1. Consequent SMC hyperplasia led to progressive medial degeneration and TAAD. In addition to diminished contractile and synthetic functions, fate-mapped SMCs displayed increased proteolysis, endocytosis, phagocytosis, and lysosomal clearance of extracellular matrix and apoptotic cells. SMCs acquired a limited repertoire of macrophage markers and functions via biogenesis of degradative organelles through an mTOR/ß-catenin/MITF-dependent pathway, but were distinguishable from conventional macrophages by an absence of hematopoietic lineage markers and certain immune effectors even in the context of hyperlipidemia. Similar mTOR activation and induction of a degradative SMC phenotype in a model of mild TAAD due to Fbn1 mutation greatly worsened disease with near-uniform lethality. The finding of increased lysosomal markers in medial SMCs from clinical TAAD specimens with hyperplasia and matrix degradation further supports the concept that proliferation of degradative SMCs within the media causes aortic disease, thus identifying mTOR-dependent phenotypic modulation as a therapeutic target for combating TAAD.

A study based on whole-genome sequencing yields a rare variant at 8q24 associated with prostate cancer.

In Western countries, prostate cancer is the most prevalent cancer of men and one of the leading causes of cancer-related death in men. Several genome-wide association studies have yielded numerous common variants conferring risk of prostate cancer. Here, we analyzed 32.5 million variants discovered by whole-genome sequencing 1,795 Icelanders. We identified a new low-frequency variant at 8q24 associated with prostate cancer in European populations, rs188140481[A] (odds ratio (OR) = 2.90; P(combined) = 6.2 × 10(-34)), with an average risk allele frequency in controls of 0.54%. This variant is only very weakly correlated (r(2) ≤ 0.06) with previously reported risk variants at 8q24, and its association remains significant after adjustment for all known risk-associated variants. Carriers of rs188140481[A] were diagnosed with prostate cancer 1.26 years younger than non-carriers (P = 0.0059). We also report results for a previously described HOXB13 variant (rs138213197[T]), confirming it as a prostate cancer risk variant in populations from across Europe.

Microsatellite markers for linkage and association studies.

During the 1990s and the first several years of this century, microsatellites or short tandem repeats were the workhorse genetic markers for hypothesis-independent studies in human genetics, facilitating genome-wide linkage studies and allelic imbalance studies. However, the rise of higher throughput and cost-effective single-nucleotide polymorphism (SNP) platforms led to the era of the SNP for genome scans. Nevertheless, it is important to note that microsatellites remain highly informative and useful measures of genomic variation for linkage and association studies. Their continued advantage in complementing SNPs lies in their greater allelic diversity than biallelic SNPs as well as in their population history, in which single-step expansion or contraction of the tandem repeat on the background of ancestral SNP haplotypes can break up common haplotypes, leading to greater haplotype diversity within the linkage disequilibrium block of interest. In fact, microsatellites have starred in association studies leading to widely replicated discoveries of type 2 diabetes (TCF7L2) and prostate cancer genes (the 8q21 region). At the end of the day, it will be important to catalog all variation, including SNPs, microsatellites, copy number variations, and polymorphic inversions in human genetic studies. This article describes the utilities of microsatellites and experimental approaches in their use.

Discovery of common variants associated with low TSH levels and thyroid cancer risk.

To search for sequence variants conferring risk of nonmedullary thyroid cancer, we focused our analysis on 22 SNPs with a P < 5 × 10(-8) in a genome-wide association study on levels of thyroid stimulating hormone (TSH) in 27,758 Icelanders. Of those, rs965513 has previously been shown to associate with thyroid cancer. The remaining 21 SNPs were genotyped in 561 Icelandic individuals with thyroid cancer (cases) and up to 40,013 controls. Variants suggestively associated with thyroid cancer (P < 0.05) were genotyped in an additional 595 non-Icelandic cases and 2,604 controls. After combining the results, three variants were shown to associate with thyroid cancer: rs966423 on 2q35 (OR = 1.34; P(combined) = 1.3 × 10(-9)), rs2439302 on 8p12 (OR = 1.36; P(combined) = 2.0 × 10(-9)) and rs116909374 on 14q13.3 (OR = 2.09; P(combined) = 4.6 × 10(-11)), a region previously reported to contain an uncorrelated variant conferring risk of thyroid cancer. A strong association (P = 9.1 × 10(-91)) was observed between rs2439302 on 8p12 and expression of NRG1, which encodes the signaling protein neuregulin 1, in blood.

Mutations in BRIP1 confer high risk of ovarian cancer.

Ovarian cancer causes more deaths than any other gynecologic malignancy in developed countries. Sixteen million sequence variants, identified through whole-genome sequencing of 457 Icelanders, were imputed to 41,675 Icelanders genotyped using SNP chips, as well as to their relatives. Sequence variants were tested for association with ovarian cancer (N of affected individuals = 656). We discovered a rare (0.41% allelic frequency) frameshift mutation, c.2040_2041insTT, in the BRIP1 (FANCJ) gene that confers an increase in ovarian cancer risk (odds ratio (OR) = 8.13, P = 2.8 × 10(-14)). The mutation was also associated with increased risk of cancer in general and reduced lifespan by 3.6 years. In a Spanish population, another frameshift mutation in BRIP1, c.1702_1703del, was seen in 2 out of 144 subjects with ovarian cancer and 1 out of 1,780 control subjects (P = 0.016). This allele was also associated with breast cancer (seen in 6/927 cases; P = 0.0079). Ovarian tumors from heterozygous carriers of the Icelandic mutation show loss of the wild-type allele, indicating that BRIP1 behaves like a classical tumor suppressor gene in ovarian cancer.

Genetic correction of PSA values using sequence variants associated with PSA levels.

Measuring serum levels of the prostate-specific antigen (PSA) is the most common screening method for prostate cancer. However, PSA levels are affected by a number of factors apart from neoplasia. Notably, around 40% of the variability of PSA levels in the general population is accounted for by inherited factors, suggesting that it may be possible to improve both sensitivity and specificity by adjusting test results for genetic effects. To search for sequence variants that associate with PSA levels, we performed a genome-wide association study and follow-up analysis using PSA information from 15,757 Icelandic and 454 British men not diagnosed with prostate cancer. Overall, we detected a genome-wide significant association between PSA levels and single-nucleotide polymorphisms (SNPs) at six loci: 5p15.33 (rs2736098), 10q11 (rs10993994), 10q26 (rs10788160), 12q24 (rs11067228), 17q12 (rs4430796), and 19q13.33 [rs17632542 (KLK3: I179T)], each with P(combined) <3 × 10(-10). Among 3834 men who underwent a biopsy of the prostate, the 10q26, 12q24, and 19q13.33 alleles that associate with high PSA levels are associated with higher probability of a negative biopsy (odds ratio between 1.15 and 1.27). Assessment of association between the six loci and prostate cancer risk in 5325 cases and 41,417 controls from Iceland, the Netherlands, Spain, Romania, and the United States showed that the SNPs at 10q26 and 12q24 were exclusively associated with PSA levels, whereas the other four loci also were associated with prostate cancer risk. We propose that a personalized PSA cutoff value, based on genotype, should be used when deciding to perform a prostate biopsy.

Ancestry-shift refinement mapping of the C6orf97-ESR1 breast cancer susceptibility locus.

We used an approach that we term ancestry-shift refinement mapping to investigate an association, originally discovered in a GWAS of a Chinese population, between rs2046210[T] and breast cancer susceptibility. The locus is on 6q25.1 in proximity to the C6orf97 and estrogen receptor alpha (ESR1) genes. We identified a panel of SNPs that are correlated with rs2046210 in Chinese, but not necessarily so in other ancestral populations, and genotyped them in breast cancer case:control samples of Asian, European, and African origin, a total of 10,176 cases and 13,286 controls. We found that rs2046210[T] does not confer substantial risk of breast cancer in Europeans and Africans (OR = 1.04, P = 0.099, and OR = 0.98, P = 0.77, respectively). Rather, in those ancestries, an association signal arises from a group of less common SNPs typified by rs9397435. The rs9397435[G] allele was found to confer risk of breast cancer in European (OR = 1.15, P = 1.2 x 10(-3)), African (OR = 1.35, P = 0.014), and Asian (OR = 1.23, P = 2.9 x 10(-4)) population samples. Combined over all ancestries, the OR was 1.19 (P = 3.9 x 10(-7)), was without significant heterogeneity between ancestries (P(het) = 0.36) and the SNP fully accounted for the association signal in each ancestry. Haplotypes bearing rs9397435[G] are well tagged by rs2046210[T] only in Asians. The rs9397435[G] allele showed associations with both estrogen receptor positive and estrogen receptor negative breast cancer. Using early-draft data from the 1,000 Genomes project, we found that the risk allele of a novel SNP (rs77275268), which is closely correlated with rs9397435, disrupts a partially methylated CpG sequence within a known CTCF binding site. These studies demonstrate that shifting the analysis among ancestral populations can provide valuable resolution in association mapping.

Genome-wide association and replication studies identify four variants associated with prostate cancer susceptibility.

We report a prostate cancer genome-wide association follow-on study. We discovered four variants associated with susceptibility to prostate cancer in several European populations: rs10934853[A] (OR = 1.12, P = 2.9 x 10(-10)) on 3q21.3; two moderately correlated (r2 = 0.07) variants, rs16902094[G] (OR = 1.21, P = 6.2 x 10(-15)) and rs445114[T] (OR = 1.14, P = 4.7 x 10(-10)), on 8q24.21; and rs8102476[C] (OR = 1.12, P = 1.6 x 10(-11)) on 19q13.2. We also refined a previous association signal on 11q13 with the SNP rs11228565[A] (OR = 1.23, P = 6.7 x 10(-12)). In a multivariate analysis using 22 prostate cancer risk variants typed in the Icelandic population, we estimated that carriers in the top 1.3% of the risk distribution are at a 2.5 times greater risk of developing the disease than members of the general population.

Common variants on 9q22.33 and 14q13.3 predispose to thyroid cancer in European populations.

In order to search for sequence variants conferring risk of thyroid cancer we conducted a genome-wide association study in 192 and 37,196 Icelandic cases and controls, respectively, followed by a replication study in individuals of European descent. Here we show that two common variants, located on 9q22.33 and 14q13.3, are associated with the disease. Overall, the strongest association signals were observed for rs965513 on 9q22.33 (OR = 1.75; P = 1.7 x 10(-27)) and rs944289 on 14q13.3 (OR = 1.37; P = 2.0 x 10(-9)). The gene nearest to the 9q22.33 locus is FOXE1 (TTF2) and NKX2-1 (TTF1) is among the genes located at the 14q13.3 locus. Both variants contribute to an increased risk of both papillary and follicular thyroid cancer. Approximately 3.7% of individuals are homozygous for both variants, and their estimated risk of thyroid cancer is 5.7-fold greater than that of noncarriers. In a study on a large sample set from the general population, both risk alleles are associated with low concentrations of thyroid stimulating hormone (TSH), and the 9q22.33 allele is associated with low concentration of thyroxin (T(4)) and high concentration of triiodothyronine (T(3)).

ASIP and TYR pigmentation variants associate with cutaneous melanoma and basal cell carcinoma.

Fair color increases risk of cutaneous melanoma (CM) and basal cell carcinoma (BCC). Recent genome-wide association studies have identified variants affecting hair, eye and skin pigmentation in Europeans. Here, we assess the effect of these variants on risk of CM and BCC in European populations comprising 2,121 individuals with CM, 2,163 individuals with BCC and over 40,000 controls. A haplotype near ASIP, known to affect a similar spectrum of pigmentation traits as MC1R variants, conferred significant risk of CM (odds ratio (OR) = 1.45, P = 1.2 x 10(-9)) and BCC (OR = 1.33, P = 1.2 x 10(-6)). The variant in TYR encoding the R402Q amino acid substitution, previously shown to affect eye color and tanning response, conferred risk of CM (OR = 1.21, P = 2.8 x 10(-7)) and BCC (OR = 1.14, P = 6.1 x 10(-4)). An eye color variant in TYRP1 was associated with risk of CM (OR = 1.15, P = 4.6 x 10(-4)). The association of all three variants is robust with respect to adjustment for the effect of pigmentation.

Multiple genetic loci for bone mineral density and fractures.

BACKGROUND: Bone mineral density influences the risk of osteoporosis later in life and is useful in the evaluation of the risk of fracture. We aimed to identify sequence variants associated with bone mineral density and fracture. METHODS: We performed a quantitative trait analysis of data from 5861 Icelandic subjects (the discovery set), testing for an association between 301,019 single-nucleotide polymorphisms (SNPs) and bone mineral density of the hip and lumbar spine. We then tested for an association between 74 SNPs (most of which were implicated in the discovery set) at 32 loci in replication sets of Icelandic, Danish, and Australian subjects (4165, 2269, and 1491 subjects, respectively). RESULTS: Sequence variants in five genomic regions were significantly associated with bone mineral density in the discovery set and were confirmed in the replication sets (combined P values, 1.2x10(-7) to 2.0x10(-21)). Three regions are close to or within genes previously shown to be important to the biologic characteristics of bone: the receptor activator of nuclear factor-kappaB ligand gene (RANKL) (chromosomal location, 13q14), the osteoprotegerin gene (OPG) (8q24), and the estrogen receptor 1 gene (ESR1) (6q25). The two other regions are close to the zinc finger and BTB domain containing 40 gene (ZBTB40) (1p36) and the major histocompatibility complex region (6p21). The 1p36, 8q24, and 6p21 loci were also associated with osteoporotic fractures, as were loci at 18q21, close to the receptor activator of the nuclear factor-kappaB gene (RANK), and loci at 2p16 and 11p11. CONCLUSIONS: We have discovered common sequence variants that are consistently associated with bone mineral density and with low-trauma fractures in three populations of European descent. Although these variants alone are not clinically useful in the prediction of risk to the individual person, they provide insight into the biochemical pathways underlying osteoporosis.

Common variants on chromosomes 2q35 and 16q12 confer susceptibility to estrogen receptor-positive breast cancer.

Familial clustering studies indicate that breast cancer risk has a substantial genetic component. To identify new breast cancer risk variants, we genotyped approximately 300,000 SNPs in 1,600 Icelandic individuals with breast cancer and 11,563 controls using the Illumina Hap300 platform. We then tested selected SNPs in five replication sample sets. Overall, we studied 4,554 affected individuals and 17,577 controls. Two SNPs consistently associated with breast cancer: approximately 25% of individuals of European descent are homozygous for allele A of rs13387042 on chromosome 2q35 and have an estimated 1.44-fold greater risk than noncarriers, and for allele T of rs3803662 on 16q12, about 7% are homozygous and have a 1.64-fold greater risk. Risk from both alleles was confined to estrogen receptor-positive tumors. At present, no genes have been identified in the linkage disequilibrium block containing rs13387042. rs3803662 is near the 5' end of TNRC9 , a high mobility group chromatin-associated protein whose expression is implicated in breast cancer metastasis to bone.

A common variant on chromosome 9p21 affects the risk of myocardial infarction.

The global endemic of cardiovascular diseases calls for improved risk assessment and treatment. Here, we describe an association between myocardial infarction (MI) and a common sequence variant on chromosome 9p21. This study included a total of 4587 cases and 12,767 controls. The identified variant, adjacent to the tumor suppressor genes CDKN2A and CDKN2B, was associated with the disease with high significance. Approximately 21% of individuals in the population are homozygous for this variant, and their estimated risk of suffering myocardial infarction is 1.64 times as great as that of noncarriers. The corresponding risk is 2.02 times as great for early-onset cases. The population attributable risk is 21% for MI in general and 31% for early-onset cases.

Genome-wide association study identifies a second prostate cancer susceptibility variant at 8q24.

Prostate cancer is the most prevalent noncutaneous cancer in males in developed regions, with African American men having among the highest worldwide incidence and mortality rates. Here we report a second genetic variant in the 8q24 region that, in conjunction with another variant we recently discovered, accounts for about 11%-13% of prostate cancer cases in individuals of European descent and 31% of cases in African Americans. We made the current discovery through a genome-wide association scan of 1,453 affected Icelandic individuals and 3,064 controls using the Illumina HumanHap300 BeadChip followed by four replication studies. A key step in the discovery was the construction of a 14-SNP haplotype that efficiently tags a relatively uncommon (2%-4%) susceptibility variant in individuals of European descent that happens to be very common (approximately 42%) in African Americans. The newly identified variant shows a stronger association with affected individuals who have an earlier age at diagnosis.

The BARD1 Cys557Ser variant and breast cancer risk in Iceland.

BACKGROUND: Most, if not all, of the cellular functions of the BRCA1 protein are mediated through heterodimeric complexes composed of BRCA1 and a related protein, BARD1. Some breast-cancer-associated BRCA1 missense mutations disrupt the function of the BRCA1/BARD1 complex. It is therefore pertinent to determine whether variants of BARD1 confer susceptibility to breast cancer. Recently, a missense BARD1 variant, Cys557Ser, was reported to be at increased frequencies in breast cancer families. We investigated the role of the BARD1 Cys557Ser variant in a population-based cohort of 1,090 Icelandic patients with invasive breast cancer and 703 controls. We then used a computerized genealogy of the Icelandic population to study the relationships between the Cys557Ser variant and familial clustering of breast cancer. METHODS AND FINDINGS: The Cys557Ser allele was present at a frequency of 0.028 in patients with invasive breast cancer and 0.016 in controls (odds ratio [OR] = 1.82, 95% confidence interval [CI] 1.11-3.01, p = 0.014). The alleleic frequency was 0.037 in a high-predisposition group of cases defined by having a family history of breast cancer, early onset of breast cancer, or multiple primary breast cancers (OR = 2.41, 95% CI 1.22-4.75, p = 0.015). Carriers of the common Icelandic BRCA2 999del5 mutation were found to have their risk of breast cancer further increased if they also carried the BARD1 variant: the frequency of the BARD1 variant allele was 0.047 (OR = 3.11, 95% CI 1.16-8.40, p = 0.046) in 999del5 carriers with breast cancer. This suggests that the lifetime probability of a BARD1 Cys557Ser/BRCA2 999del5 double carrier developing breast cancer could approach certainty. Cys557Ser carriers, with or without the BRCA2 mutation, had an increased risk of subsequent primary breast tumors after the first breast cancer diagnosis compared to non-carriers. Lobular and medullary breast carcinomas were overrepresented amongst Cys557Ser carriers. We found that an excess of ancestors of contemporary carriers lived in a single county in the southeast of Iceland and that all carriers shared a SNP haplotype, which is suggestive of a founder event. Cys557Ser was found on the same SNP haplotype background in the HapMap Project CEPH sample of Utah residents. CONCLUSIONS: Our findings suggest that BARD1 Cys557Ser is an ancient variant that confers risk of single and multiple primary breast cancers, and this risk extends to carriers of the BRCA2 999del5 mutation.

A common variant associated with prostate cancer in European and African populations.

With the increasing incidence of prostate cancer, identifying common genetic variants that confer risk of the disease is important. Here we report such a variant on chromosome 8q24, a region initially identified through a study of Icelandic families. Allele -8 of the microsatellite DG8S737 was associated with prostate cancer in three case-control series of European ancestry from Iceland, Sweden and the US. The estimated odds ratio (OR) of the allele is 1.62 (P = 2.7 x 10(-11)). About 19% of affected men and 13% of the general population carry at least one copy, yielding a population attributable risk (PAR) of approximately 8%. The association was also replicated in an African American case-control group with a similar OR, in which 41% of affected individuals and 30% of the population are carriers. This leads to a greater estimated PAR (16%) that may contribute to higher incidence of prostate cancer in African American men than in men of European ancestry.

A genome-wide study of allelic imbalance in human testicular germ cell tumors using microsatellite markers.

Testicular germ cell tumors (TGCT) arise by multistep carcinogenesis pathways involving selective losses and gains of chromosome material. To locate cancer genes underlying this selection, we performed a genome-wide study of allelic imbalance (AI) in 32 tumors, using 710 microsatellite markers. The highest prevalence of AI was found at 12p, in line with previous studies finding consistent gain of the region in TGCTs. High frequency of AI was also observed at chromosome arms 4p, 9q, 10p, 11q, 11p, 13q, 16q, 18p, and 22q. Within 39 candidate regions identified by mapping of smallest regions of overlap (SROs), the highest frequency of AI was at 12p11.21 approximately p11.22 (62%), 12p12.1 approximately p13.1 (53%), 12p13.1 approximately p13.2 (53%), 11q14.1 approximately q14.2 (53%), 11p13 approximately p14.3 (47%), 9q21.13 approximately q21.32 (47%), and 4p15.1 approximately p15.2 (44%). Two genes known to be involved in cancer reside in these regions, ETV6 at 12p13.2 (TEL oncogene) and WT1 at 11p13. We also found a significant association (P = 0.02) between AI at 10q21.1 approximately q22.2 and higher clinical stage. This study contributes to the ongoing search for genes involved in transformation of germ cells and provides a useful reference point to previous studies using cytogenetic techniques to map chromosome changes in TGCTs.

Profiling of genes expressed in peripheral blood mononuclear cells predicts glucocorticoid sensitivity in asthma patients.

Gene expression profiles were examined in freshly isolated peripheral blood mononuclear cells (PBMC) from two independent cohorts (training and test sets) of glucocorticoid (GC)-sensitive (n = 64) and GC-resistant (n = 42) asthma patients in search of genes that accurately predict responders and nonresponders to inhaled corticosteroids. A total of 11,812 genes were examined with high-density oligonucleotide microarrays in both resting PBMC (106 patients) and cells treated in vitro with IL-1beta and TNF-alpha combined (88 patients), with or without GC. A total of 5,011 genes were expressed at significant levels in the PBMC, and 1,334 of those were notably up-regulated or down-regulated by IL-1beta/TNF-alpha treatment. The expression changes of 923 genes were significantly reversed in GC responders in the presence of GC. The expression pattern of 15 of these 923 genes that most accurately separated GC responders (n = 26) from the nonresponders (n = 18) in the training set, based on the weighted voting algorithm, predicted the independent test set of equal size with 84% accuracy. The expression accuracy of these genes was confirmed by real-time-quantitative PCR, wherein 11 of the 15 genes predicted GC sensitivity at baseline with 84% accuracy, with one gene predicting at 81% in an independent cohort of 79 patients. We conclude that we have uncovered gene expression profiles in PBMC that predict clinical response to inhaled GC therapy with meaningful accuracy. Upon validation in an independent study, these results support the development of a diagnostic test to guide GC therapy in asthma patients.