Tumor characterization and stratification by integrated molecular profiles reveals essential pan-cancer features

Background

Identification of tumor heterogeneity and genomic similarities across different cancer types is essential to the design of effective stratified treatments and for the discovery of treatments that can be extended to different types of tumors. However, systematic investigations on comprehensive molecular profiles have not been fully explored to achieve this goal.

Results

Here, we performed a network-based integrative pan-cancer genomic analysis on >3000 samples from 12 cancer types to uncover novel stratifications among tumors. Our study not only revealed recurrently reported cross-cancer similarities, but also identified novel ones. The macro-scale stratification demonstrates strong clinical relevance and reveals consistent risk tendency among cancer types. The micro-scale stratification shows essential pan-cancer heterogeneity with subgroup-specific gene network characteristics and biological functions.

Conclusions

In summary, our comprehensive network-based pan-cancer stratification provides valuable information about inter- and intra- cancer stratification for patient clinical assessments and therapeutic strategies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1687-x) contains supplementary material, which is available to authorized users.     

Identification of murine poxvirus-specific CD8+ CTL epitopes with distinct functional profiles

Murine T cell epitopes against vaccinia virus (VV) have not been characterized to date in part due to the large and complex genome of VV. We have identified and characterized two CD8+ T cell epitopes on the A47L (modified VV Ankara strain (MVA)-029) and J6R (MVA-043) proteins of VV that are Db and Kb restricted, respectively. Following i.p. immunization with VV New York City Board of Health (NYCBH) strain, MVA-029 peptide-stimulated splenocytes secreted IFN-gamma from 7 days to 7 mo postimmunization, and virus-stimulated effectors were also able to lyse MVA-029-pulsed target cells at the same time points. In contrast, MVA-043 peptide-stimulated splenocytes secreted very low levels of IFN-gamma only at day 7 but maintained the ability to lyse target cells up to 2 mo postimmunization. Both MVA-029 and MVA-043 peptide-stimulated lymph node cells degranulated similarly as assessed by Ag-induced CD107 expression. T cell responses to whole-virus stimulation remained robust and steady during the acute and memory T cell response to VV. Identification of T cell epitopes on VV will enable further studies to increase our understanding of the role of CD8+ T cells in VV infection and assist in the design of new protective strategies.     

Populations of Populus angustifolia have evolved distinct metabolic profiles that influence their surrounding soil

Plant and Soil - Plant-microbial-soil interactions are key to understanding plant community succession, invasion success, patterns of biodiversity and aspects of ecosystem function. Yet root and...     

Human GPR17 missense variants identified in metabolic disease patients have distinct downstream signaling profiles

GPR17 is a G-protein-coupled receptor (GPCR) implicated in the regulation of glucose metabolism and energy homeostasis. Such evidence is primarily drawn from mouse knockout studies and suggests GPR17 as a potential novel therapeutic target for the treatment of metabolic diseases. However, links between human GPR17 genetic variants, downstream cellular signaling, and metabolic diseases have yet to be reported. Here, we analyzed GPR17 coding sequences from control and disease cohorts consisting of individuals with adverse clinical metabolic deficits including severe insulin resistance, hypercholesterolemia, and obesity. We identified 18 nonsynonymous GPR17 variants, including eight variants that were exclusive to the disease cohort. We characterized the protein expression levels, membrane localization, and downstream signaling profiles of nine GPR17 variants (F43L, V96M, V103M, D105N, A131T, G136S, R248Q, R301H, and G354V). These nine GPR17 variants had similar protein expression and subcellular localization as wild-type GPR17; however, they showed diverse downstream signaling profiles. GPR17-G136S lost the capacity for agonist-mediated cAMP, Ca²⁺, and β-arrestin signaling. GPR17-V96M retained cAMP inhibition similar to GPR17-WT, but showed impaired Ca²⁺ and β-arrestin signaling. GPR17-D105N displayed impaired cAMP and Ca²⁺ signaling, but unaffected agonist-stimulated β-arrestin recruitment. The identification and functional profiling of naturally occurring human GPR17 variants from individuals with metabolic diseases revealed receptor variants with diverse signaling profiles, including differential signaling perturbations that resulted in GPCR signaling bias. Our findings provide a framework for structure–function relationship studies of GPR17 signaling and metabolic disease.     

Human germline differentiation charts a new course

Understanding the molecular events of reproduction requires a system to differentiate human pluripotent stem cells to germline cells (gametes) in vitro. Such a system is not only critical to unlock the secrets of germline development; it may also allow screening for environmental agents that affect gametogenesis. Two recent papers, one in this issue of The EMBO Journal, have developed complementary approaches for generating human germline cells with unprecedented efficiency from pluripotent stem cells (Sugawa et al, 2015; Irie et al, 2015). This work illustrates the power and limitations of extrapolating molecular pathways for lineage differentiation from mice to humans and illuminates the importance of using human cell‐based models to study reproductive health.     

Human germline mutation in the factor IX gene

The molecular epidemiology of factor IX germline mutations in patients with hemophilia B has been studied in detail because it is an advantageous model for analyzing recent germline mutations in humans. It is estimated that mutations have been defined in the majority of nucleotides that are the target for mutation. The likelihood that a factor IX missense mutation will cause disease correlates with the degree of evolutionary conservation of the amino acid. Mutation rates per base-pair have been estimated after careful consideration and correction for biases, predicting about 76 de novo mutations per generation per individual resulting in 0.3 deleterious changes. The male-to-female sex ratio of mutation varies with the type of mutation. There is evidence for a maternal age effect and an excess of non-CpG G:C to A:T transitions. The factor IX mutation pattern is similar among geographically, racially and ethnically diverse human populations. The data support primarily endogenous mechanisms of germline mutation in the factor IX gene. Mutations at splice junctions are compatible with simple rules for predicting disease causing mutations.     

Integrative epigenetic and genetic pan-cancer somatic alteration portraits

Genetic and epigenetic alterations are required for carcinogenesis and the mutation burden across tumor types has been investigated. Here, we investigate epigenetic alterations with a novel measure of global DNA methylation dysregulation, the methylation dysregulation index (MDI), across 14 cancer types in The Cancer Genome Atlas (TCGA) database. DNA methylation data—obtained using Illumina HumanMethylation450 BeadChip—was accessed from TCGA. We calculated the MDI in 14 tumor types (n = 5,592 tumors), using adjacent normal tissues (n = 701) from each tumor site. Copy number alteration, and mutation burden were retrieved from cBioportal (n = 5,152). We tested the relation of subject MDI across tumors and with age, gender, tumor stage, estimated tumor purity, and copy number alterations for both overall MDI and genomic-context-specific MDI. We also investigated the top most dysregulated loci shared across tumor types. There was a broad range of extent in methylation dysregulation across tumor types (P < 2.2E-16). However, a consistent pattern of methylation dysregulation stratified by genomic context was observed across tumor types where the highest dysregulation occurred at non-CpG island regions. Considering other summary measures of somatic alteration, MDI was correlated with copy number alterations but not with mutation burden. Using the top dysregulated CpG sites in common across tumors, 4 classes of cancer types were observed, and the functional consequences of these alterations to gene expression were confirmed. This work identified the global DNA methylation dysregulation patterns across 14 cancer types showing a higher impact for the non-CpG island areas. The most dysregulated loci across cancer types identified common clusters across cancer types that may have implications for future treatment and prevention measures.     

Clinical characteristics and outcomes of hospitalized older patients with distinct risk profiles for functional decline: a prospective cohort study

Background

The aim of this research was to study the clinical characteristics and mortality and disability outcomes of patients who present distinct risk profiles for functional decline at admission.

Methods

Multicenter, prospective cohort study conducted between 2006 and 2009 in three hospitals in the Netherlands in consecutive patients of ≥65 years, acutely admitted and hospitalized for at least 48 hours. Nineteen geriatric conditions were assessed at hospital admission, and mortality and functional decline were assessed until twelve months after admission. Patients were divided into risk categories for functional decline (low, intermediate or high risk) according to the Identification of Seniors at Risk-Hospitalized Patients.

Results

A total of 639 patients were included, with a mean age of 78 years. Overall, 27%, 33% and 40% of the patients were at low, intermediate or high risk, respectively, for functional decline. Low-risk patients had fewer geriatric conditions (mean 2.2 [standard deviation [SD] 1.3]) compared with those at intermediate (mean 3.8 [SD 2.1]) or high risk (mean 5.1 [SD 1.8]) (p<0.001). Twelve months after admission, 39% of the low-risk group had an adverse outcome, compared with 50% in the intermediate risk group and 69% in the high risk group (p<0.001).

Conclusion

By using a simple risk assessment instrument at hospital admission, patients at low, intermediate or high risk for functional decline could be identified, with distinct clinical characteristics and outcomes. This approach should be tested in clinical practice and research and might help appropriately tailor patient care.     

Identification and functional characterization of the distinct plant pectin esterases PAE8 and PAE9 and their deletion mutants

Main conclusion

PAE8 and PAE9 have pectin acetylesterase activity and together remove one-third of the cell wall acetate associated with pectin formation in Arabidopsis leaves. In pae8 and pae9 mutants, substantial amounts of acetate accumulate in cell walls. In addition, the inflorescence stem height is decreased. Pectic polysaccharides constitute a significant part of the primary cell walls in dicotyledonous angiosperms. This diverse group of polysaccharides has been implicated in several physiological processes including cell-to-cell adhesion and pathogenesis. Several pectic polysaccharides contain acetyl-moieties directly affecting their physical properties such as gelling capacity, an important trait for the food industry. In order to gain further insight into the biological role of pectin acetylation, a reverse genetics approach was used to investigate the function of genes that are members of the Pectin AcetylEsterase gene family (PAE) in Arabidopsis. Mutations in two members of the PAE family (PAE8 and PAE9) lead to cell walls with an approximately 20 % increase in acetate content. High-molecular-weight fractions enriched in pectic rhamnogalacturonan I (RGI) extracted from the mutants had increased acetate content. In addition, the pae8 mutant displayed increased acetate content also in low-molecular-weight pectic fractions. The pae8/pae9-2 double mutant exhibited an additive effect by increasing wall acetate content by up to 37 %, suggesting that the two genes are not redundant and act on acetyl-substituents of different pectic domains. The pae8 and pae8/pae9-2 mutants exhibit reduced inflorescence growth underscoring the role of pectic acetylation in plant development. When heterologously expressed and purified, both gene products were shown to release acetate from the corresponding mutant pectic fractions in vitro. PAEs play a significant role in modulating the acetylation state of pectic polymers in the wall, highlighting the importance of apoplastic metabolism for the plant cell and plant growth.     

Cytochrome c oxidase-deficient patients have distinct subunit assembly profiles

Cytochrome c oxidase (COX) deficiency is the most common respiratory chain defect in childhood and is clinically heterogeneous. We report a study of six patients with COX deficiencies. Two of the patients had as yet undefined defects, three patients had Surf-1 mutations, and one patient had a 15-base pair deletion in the COX III subunit. We show that quantitative measurements of steady-state levels of subunits by monoclonal antibody reactivity, when used in combination with a discontinuous sucrose gradient methods, provide an improved diagnosis of COX deficiencies by distinguishing between kinetic, stability, and assembly defects. The two mutants of undefined etiology had a full complement of subunits with one stable and the other partially unstable to detergent solubilization. Both are likely to carry mutations in nuclear-encoded subunits of the complex. The three Surf-1 mutants and the COX III mutant each had reduced steady-state levels of subunits but variable associations of the residual subunits. This information, as well as aiding in diagnosis, helps in understanding the genotype-phenotype relationships of COX deficiencies and provides insight into the mechanism of assembly of the enzyme complex.     

Modern and prebiotic amino acids support distinct structural profiles in proteins

The earliest proteins had to rely on amino acids available on early Earth before the biosynthetic pathways for more complex amino acids evolved. In extant proteins, a significant fraction of the ‘late’ amino acids (such as Arg, Lys, His, Cys, Trp and Tyr) belong to essential catalytic and structure-stabilizing residues. How (or if) early proteins could sustain an early biosphere has been a major puzzle. Here, we analysed two combinatorial protein libraries representing proxies of the available sequence space at two different evolutionary stages. The first is composed of the entire alphabet of 20 amino acids while the second one consists of only 10 residues (ASDGLIPTEV) representing a consensus view of plausibly available amino acids through prebiotic chemistry. We show that compact conformations resistant to proteolysis are surprisingly similarly abundant in both libraries. In addition, the early alphabet proteins are inherently more soluble and refoldable, independent of the general Hsp70 chaperone activity. By contrast, chaperones significantly increase the otherwise poor solubility of the modern alphabet proteins suggesting their coevolution with the amino acid repertoire. Our work indicates that while both early and modern amino acids are predisposed to supporting protein structure, they do so with different biophysical properties and via different mechanisms.     

A GGCAGG motif in minisatellites affecting their germline instability

Mouse and human genomes contain hypervariable DNA regions consisting of tandem repeats of a short sequence referred to as minisatellites. This variation is thought to arise through processes such as unequal crossover or replication slippage. A mo-1 minisatellite probe comprising a 14-base pair repeat sequence reveals many polymorphic fragments even in DNA of BALB/c sublines. Oligonucleotide probes with single base substitution in the mo-1 have been synthesized and used for assessing sequence involved in generation of polymorphisms. The results indicate that the loci containing mo-1 homologues with mutation in the GGCAGG sequence are monomorphic despite the other mutants showing polymorphism. Reciprocally, locus-specific polymorphic clones, Pc-1 and Pc-2, have been isolated with hybridization to mo-1, and both are shown to contain repeated sequence comprising the GGCAGG sequence. They reveal high mutation rates of 8.8% and 3.3% per gamete, respectively. These results strongly suggest that the motif contributes to the germline instability of minisatellites.     

Human pex19p binds peroxisomal integral membrane proteins at regions distinct from their sorting sequences

The molecular machinery underlying peroxisomal membrane biogenesis is not well understood. The observation that cells deficient in the peroxins Pex3p, Pex16p, and Pex19p lack peroxisomal membrane structures suggests that these molecules are involved in the initial stages of peroxisomal membrane formation. Pex19p, a predominantly cytosolic protein that can be farnesylated, binds multiple peroxisomal integral membrane proteins, and it has been suggested that it functions as a soluble receptor for the targeting of peroxisomal membrane proteins (PMPs) to the peroxisome. An alternative view proposes that Pex19p functions as a chaperone at the peroxisomal membrane. Here, we show that the peroxisomal sorting determinants and the Pex19p-binding domains of a number of PMPs are distinct entities. In addition, we extend the list of peroxins with which human Pex19p interacts to include the PMP Pex16p and show that Pex19p's CaaX prenylation motif is an important determinant in the affinity of Pex19p for Pex10p, Pex11pbeta, Pex12p, and Pex13p.     

Ovaries and germline cysts and their evolution in Dermaptera (Insecta)

We studied the ovary structure and initial stages of oogenesis in 15 representatives of several dermapteran taxa, including the epizoic Arixeniina. In all examined species, the ovaries are meroistic–polytrophic. The ovaries of the basal taxa (‘Pygidicranidae’, ‘Diplatyidae’, and Labiduridae) are composed of elongated ovarioles, attached to short lateral oviducts. In these groups, ovarioles contain several (more than 30) ovarian follicles in a linear arrangement. In the Eudermaptera, the ovaries are composed of 1–6 (Spongiphoridae) or 20–40 (Forficulidae, Chelisochidae) short ovarioles (containing 2 ovarian follicles only) that open to strongly elongated lateral oviducts. In all investigated dermapterans, the ovarian follicles are composed of two germline cells only: an oocyte and a polyploid nurse cell that are covered by a simple follicular epithelium. Our studies indicate that despite a rather unique morphology of the ovarian follicles in the examined species, the processes leading to the formation of the oocyte and nurse cell units are significantly different in basal versus derived taxa.The ovaries of Arixenia esau are composed of 3 short ovarioles attached to a strongly dilated lateral oviduct, ‘the uterus’, containing developing embryos. Histological analysis suggests that the origin of the oocyte and nurse cell units in this species follows the pattern described in eudermapterans.The interpretation of our results in an evolutionary context supports the monophyly of the Dermaptera and Eudermaptera, and the inclusion of the Arixeniina and Hemimerina in the latter taxon.     

Suppression of distinct ovo phenotypes in the Drosophila female germline by maleless− and Sex-lethalM

Mutations in ovo result in several different phenotypes, which we show are due to the regulation of distinct developmental pathways. Two X (female) germ cells require ovo⁺ activity for viability, but 1X (male) germ cells do not. In our study, we observed suppression of the ovo germline-lethality phenotype in loss-of-function maleless (mle) females indicating that ovo⁺ and mle⁺ have opposing effects in female germ cells; or that they are hierarchically related. Gain-of-function Sex-lethal (Sxl) alleles and male specific lethal-2 alleles did not suppress the ovo germline death phenotype. Many of the surviving germ cells in females mutant for both ovo and mle showed ovarian tumors. In contrast to the germline viability phenotype, we did observe suppression of the tumor phenotype in females heterozygous for gain-of-function alleles of Sxl. Further, females mutant for some hypomorphic ovo alleles were rendered fertile by Sxl gain-of-function alleles. Thus, ovo⁺ is required for at least two distinct functions, one involving mle⁺, and one mediated by Sxl⁺ gene products. The existence of ovo⁺ functions independent of mle⁺ and Sxl⁺ is likely. Dev. Genet. 23:335–346, 1998. © 1998 Wiley-Liss, Inc.