Deciphering oxygen distribution and hypoxia profiles in the tumor microenvironment: a data-driven mechanistic modeling approach.

Objective. The distribution of hypoxia within tissues plays a critical role in tumor diagnosis and prognosis. Recognizing the significance of tumor oxygenation and hypoxia gradients, we introduce mathematical frameworks grounded in mechanistic modeling approaches for their quantitative assessment within a tumor microenvironment. By utilizing known blood vasculature, we aim to predict hypoxia levels across different tumor types.Approach. Our approach offers a computational method to measure and predict hypoxia using known blood vasculature. By formulating a reaction-diffusion model for oxygen distribution, we derive the corresponding hypoxia profile.Main results. The framework successfully replicates observed inter- and intra-tumor heterogeneity in experimentally obtained hypoxia profiles across various tumor types (breast, ovarian, pancreatic). Additionally, we propose a data-driven method to deduce partial differential equation models with spatially dependent parameters, which allows us to comprehend the variability of hypoxia profiles within tissues. The versatility of our framework lies in capturing diverse and dynamic behaviors of tumor oxygenation, as well as categorizing states of vascularization based on the dynamics of oxygen molecules, as identified by the model parameters.Significance. The proposed data-informed mechanistic method quantitatively assesses hypoxia in the tumor microenvironment by integrating diverse histopathological data and making predictions across different types of data. The framework provides valuable insights from both modeling and biological perspectives, advancing our comprehension of spatio-temporal dynamics of tumor oxygenation.

A refined molecular taxonomy of breast cancer.

The current histoclinical breast cancer classification is simple but imprecise. Several molecular classifications of breast cancers based on expression profiling have been proposed as alternatives. However, their reliability and clinical utility have been repeatedly questioned, notably because most of them were derived from relatively small initial patient populations. We analyzed the transcriptomes of 537 breast tumors using three unsupervised classification methods. A core subset of 355 tumors was assigned to six clusters by all three methods. These six subgroups overlapped with previously defined molecular classes of breast cancer, but also showed important differences, notably the absence of an ERBB2 subgroup and the division of the large luminal ER+ group into four subgroups, two of them being highly proliferative. Of the six subgroups, four were ER+/PR+/AR+, one was ER-/PR-/AR+ and one was triple negative (AR-/ER-/PR-). ERBB2-amplified tumors were split between the ER-/PR-/AR+ subgroup and the highly proliferative ER+ LumC subgroup. Importantly, each of these six molecular subgroups showed specific copy-number alterations. Gene expression changes were correlated to specific signaling pathways. Each of these six subgroups showed very significant differences in tumor grade, metastatic sites, relapse-free survival or response to chemotherapy. All these findings were validated on large external datasets including more than 3000 tumors. Our data thus indicate that these six molecular subgroups represent well-defined clinico-biological entities of breast cancer. Their identification should facilitate the detection of novel prognostic factors or therapeutical targets in breast cancer.

Lobular and ductal carcinomas of the breast have distinct genomic and expression profiles.

Invasive ductal carcinomas (IDCs) and invasive lobular carcinomas (ILCs) are the two major pathological types of breast cancer. Epidemiological and histoclinical data suggest biological differences, but little is known about the molecular alterations involved in ILCs. We undertook a comparative large-scale study by both array-compared genomic hybridization and cDNA microarray of a set of 50 breast tumors (21 classic ILCs and 29 IDCs) selected on homogeneous histoclinical criteria. Results were validated on independent tumor sets, as well as by quantitative RT-PCR. ILCs and IDCs presented differences at both the genomic and expression levels with ILCs being less rearranged and heterogeneous than IDCs. Supervised analysis defined a 75-BACs signature discriminating accurately ILCs from IDCs. Expression profiles identified two subgroups of ILCs: typical ILCs ( approximately 50%), which were homogeneous and displayed a normal-like molecular pattern, and atypical ILCs, more heterogeneous with features intermediate between ILCs and IDCs. Supervised analysis identified a 75-gene expression signature that discriminated ILCs from IDCs, with many genes involved in cell adhesion, motility, apoptosis, protein folding, extracellular matrix and protein phosphorylation. Although ILCs and IDCs share common alterations, our data show that ILCs and IDCs could be distinguished on the basis of their genomic and expression profiles suggesting that they evolve along distinct genetic pathways.

Genetic profiling of chromosome 1 in breast cancer: mapping of regions of gains and losses and identification of candidate genes on 1q.

Chromosome 1 is involved in quantitative anomalies in 50-60% of breast tumours. However, the structure of these anomalies and the identity of the affected genes remain to be determined. To characterise these anomalies and define their consequences on gene expression, we undertook a study combining array-CGH analysis and expression profiling using specialised arrays. Array-CGH data showed that 1p was predominantly involved in losses and 1q almost exclusively in gains. Noticeably, high magnitude amplification was infrequent. In an attempt to fine map regions of copy number changes, we defined 19 shortest regions of overlap (SROs) for gains (one at 1p and 18 at 1q) and of 20 SROs for losses (all at 1p). These SROs, whose sizes ranged from 170 kb to 3.2 Mb, represented the smallest genomic intervals possible based on the resolution of our array. The elevated incidence of gains at 1q, added to the well-established concordance between DNA copy increase and augmented RNA expression, made us focus on gene expression changes at this chromosomal arm. To identify candidate oncogenes, we studied the RNA expression profiles of 307 genes located at 1q using a home-made built cDNA array. We identified 30 candidate genes showing significant overexpression correlated to copy number increase. In order to substantiate their involvement, RNA expression levels of these candidate genes were measured by quantitative (Q)-RT-PCR in a panel of 25 breast cancer cell lines previously typed by array-CGH. Q-PCR showed that 11 genes were significantly overexpressed in the presence of a genomic gain in these cell lines, and 20 overexpressed when compared to normal breast.

Chordin is underexpressed in ovarian tumors and reduces tumor cell motility.

Ovarian cancers mostly derive from the monolayer epithelium that covers the ovary. There are currently very few molecular clues to the etiology of this cancer. Bone morphogenetic proteins (BMPs) are required for follicular development and female fertility and are expressed in the ovarian surface epithelium (OSE). We previously reported the expression of human chordin (CHRD), a BMP extracellular regulator, in the ovary. Here we show that CHRD is underexpressed in epithelium ovary cancer and epithelial cancer cell lines as compared with normal tissues and OSE, respectively. Besides, we detected BMP expression in all ovarian cell lines analyzed. To determine the functional relevance of the absence of CHRD mRNA in tumors and cancer cell lines, we studied the effects of CHRD on two cancer cell lines, BG1 and PEO14. Migratory and invasive properties were greatly reduced, whereas cell adhesion to the support was enhanced. In addition, we detected chordin (Chrd) expression in OSE of rat ovaries in a pattern similar to that of BMP4. Altogether, these results suggest that CHRD could participate in regulating BMP activity in normal OSE physiology, and that its mis-expression in OSE may facilitate cancer incidence and/or progression.

The human chordin gene encodes several differentially expressed spliced variants with distinct BMP opposing activities.

During early embryogenesis of both vertebrates and invertebrates, antagonism between bone morphogenetic proteins (BMPs) and several unrelated secreted factors including Chordin (Chd) is a general mechanism by which the dorso-ventral axis is established. High affinity binding of Chd sequesters the BMP ligands in the extracellular space, preventing interactions with their membrane receptors. Another level of regulation consists in processing of vertebrate Chd or its Drosophila counterpart Sog by astacine metalloproteases like Xolloid-BMP-1/Tolloid, respectively, which releases an active BMP. Recently, it was shown that cleavage of Sog by Tolloid could generate novel BMP inhibitory activity and that sog is also capable of stimulation of BMP activity in a tolloid-dependant way. Activity and/or cleavage of Chd/Sog are influenced by other secreted factors like twisted gastrulation. In this study, we have cloned cDNAs of the human chordin gene (CHRD) and characterized alternative splice variants that code for C-truncated forms of the protein. We have found that CHRD is expressed in fetal as well as in adult tissues with relatively high levels in liver, cerebellum and female genital tract, suggesting functions in late embryogenesis and adult physiology. We also show that spliced variants are present with specific patterns in various tissues. When tested in an axis-duplication assay in Xenopus, we find that these variants can antagonize BMP activity. Altogether, these results suggest that, in addition to processing by metalloproteases, alternative splicing (AS) is another mechanism by which sub-products of CHRD can be generated to influence BMP activity in different developmental and physiological situations.

[Genomic profiling: from molecular cytogenetics to DNA arrays]. / Le typage génomique: de la cytogénétique moléculaire aux puces á ADN.

Genetic instability results, in a large majority of solid tumors, in deep chromosomal rearrangements. However, because chromosomal instability produces highly complex caryotypes, rarely showing stereotypic aberrations, it has not been possible to characterize solid cancers according to specific patterns of chromosomal rearrangements. This contrasts with the situation in hematological malignancies, where cytogenetics has allowed to lay out the basis of a renewed classification. New insights have been brought by the development of comparative genomic hybridization (CGH). This molecular cytogenetics approach was originally devised to detect regions in the genome of tumor cells undergoing quantitative changes, i.e. gains or losses of copy numbers. The large body of studies based on CGH show that solid tumors undergo frequent gains and losses and that every chromosomes show at least one region of anomaly. Furthermore, different tumor types present distinct CGH patterns of gains and losses. These observations favor the idea that it may be possible to type human solid cancers according to their patterns of genomic aberrations. However, despite the fact that a number of CGH based studies present data suggesting that different tumor types or cancers at different stages of evolution show distinct patterns of gains and losses, it has proven difficult to be conclusive. This can be mainly attributed to the lack of spatial resolution of CGH. Indeed, CGH uses metaphase chromosomes as hybridization targets and therefore its resolution is at the level of chromosomal banding. The recent adaptation of DNA array technology to CGH will allow to pass this limitation. In DNA array based CGH (array-CGH) metaphase chromosomes have been replaced by spots of cloned DNA. These DNA clones may either be genomic (BACs, YACs or cosmids) or coding (cDNAs). The resolution of array-CGH is therefore determined by the size of the cloned DNA insert (100 Kb for BACs, 1-2 kb for cDNAs). Data corresponding to each of these clones is or will be in a near future linked to DNA sequence data. Hence, in a near future, array-CGH will allow to increase the resolution from a cytogenetic level to a molecular level. Finally, because array technology is highly adaptable to automation, going from classical CGH to array-CGH will produce a quantum leap in throughput.

Molecular cytogenetic analysis of breast cancer cell lines.

The extensive chromosome rearrangements of breast carcinomas must contribute to tumour development, but have been largely intractable to classical cytogenetic banding. We report here the analysis by 24-colour karyotyping and comparative genomic hybridization (CGH) of 19 breast carcinoma cell lines and one normal breast epithelial cell line, which provide model examples of karyotype patterns and translocations present in breast carcinomas. The CGH was compared with CGH of 106 primary breast cancers. The lines varied from perfectly diploid to highly aneuploid. Translocations were very varied and over 98% were unbalanced. The most frequent in the carcinomas were 8;11 in five lines; and 8;17, 1;4 and 1;10 in four lines. The most frequently involved chromosome was 8. Several lines showed complex multiply-translocated chromosomes. The very aneuploid karyotypes appeared to fall into two groups that evolved by different routes: one that steadily lost chromosomes and at one point doubled their entire karyotype; and another that steadily gained chromosomes, together with abnormalities. All karyotypes fell within the range seen in fresh material and CGH confirmed that the lines were broadly representative of fresh tumours. The karyotypes provide a resource for the cataloguing and analysis of translocations in these tumours, accessible at http://www.path.cam.ac.uk/ approximately pawefish.

Human E2F5 gene is oncogenic in primary rodent cells and is amplified in human breast tumors.

E2F transcription factors (E2F1 to 6) are central players in the control of animal cell proliferation as regulators of genes involved in cell cycle progression and in transformation. In this report, we have investigated the potential involvement of the E2F5 gene in tumorigenesis. We show that E2F5 can promote the formation of morphologically transformed foci in primary baby rat kidney cells (BRK) when it is overexpressed in the presence of its heterodimeric partner DP1 and activated RAS. This suggests that E2F5 behaves like a MYC-type cooperating oncogene in functional assays, prompting us to monitor potential amplifications of the E2F5 gene in primary human tumors. We mapped the human E2F5 gene to 8q21.1-21.3 equidistant from the MOS (8q12) and MYC (8q24) oncogenes. Since the long arm of chromosome 8 is frequently the site of increased gene copy number (ICN) in breast cancer, we screened 442 breast tumor DNAs for gains of E2F5, MOS, and MYC genes. The three genes showed ICN, albeit at variable incidence and levels of amplification, with the ICN of E2F5 occurring concomitantly with those of MOS and/or MYC in almost half of the cases. Moreover, a marked increase of the 2. 5-kb E2F5 transcript was also detected in some tumors and tumor cell lines. In conclusion, the evidence that sustained unregulated expression of E2F5 can cooperate with other oncogenes to promote cell transformation in functional assays, together with the detection of chromosomal amplifications and overexpressions of the E2F5 gene in breast tumors, provides the first indications that E2F5 deregulation may have a role in human tumor development.

17q21-q25 aberrations in breast cancer: combined allelotyping and CGH analysis reveals 5 regions of allelic imbalance among which two correspond to DNA amplification.

Chromosome 17q is frequently rearranged in breast cancer. Allelotyping studies have proposed the existence of at least four regions of allelic imbalance (AI). Here we present a study combining allelotyping using 19 CA repeat markers mapping in the 17q21-25 region and molecular cytogenetics (CGH and FISH). Allelotyping was undertaken on 178 pairs of cognate tumor and normal DNA in order to determine the number of regions of AI and define the shortest overlaps. AI ranged from 34-54% of the informative cases according to the marker and, overall, 66% of the tumors presented AI at one of the markers tested. Analysis of the patterns of imbalances revealed at least five common regions of imbalance respectively defined by markers: D17S855, which is intragenic of BRCA1 (SRO 1), D17S1607 (SRO 2), D17S1855 (SRO 3), between D17S789 and D17S785 (SRO 4) and D17S784 (SRO 5). In order to characterize the nature of the genetic events revealed by allelotyping we performed CGH analysis on a subset of 43 tumors presenting variable patterns of imbalance. CGH showed that AI at 17q could represent four different types of genetic events: loss of chromosome 17, gain of 17q, gain of 17q22-q24, loss of 17q11-q21 and/or 17q25-qter. Some of these anomalies could occur concomitantly within the same tumor. Since 35% of the tumors analysed by CGH presented gains, these data indicated that AI at 17q were not solely indicative of losses of genetic material and could also represent DNA amplification. Gains were most commonly observed in the 17q23-q24 regions. This suggested that AI in SRO 2 and SRO 3 corresponded to DNA amplification. To assess this, we isolated BAC clones by PCR screening for markers D17S1607 and D17S1855 and used these in FISH experiments on six breast tumor cell lines and 14 breast cancer specimens. FISH results showed that both D17S1607 and D17S1855 were frequently involved in DNA amplification (8-30 copies). Altogether, our data show that allelotyping can be efficiently used in amplicon mapping. Clinico-pathological correlations indicated that imbalance at 17q preferentially occurred in high grade, PR- and ERBB2 amplified tumors.

[Application of the PRINS technique for chromosome examination in fetal cells present in maternal blood]. / Application de la technique PRINS à l'examen chromosomique des cellules foetales présentes dans le sang maternal.

Both detection and chromosomal analysis of fetal cells present in the maternal circulation can be performed using Histopaque double density gradient centrifugation followed by primed in situ (PRINS) labeling technique. This approach has been tested on blood samples from 15 pregnant women and 6 control donors with primers specific for chromosomes 9, X and Y. The cell separation technique allows recovery of both mononuclear cells and polynuclear cells with a 97% efficiency. PRINS labeling was successful in 100% cells from control blood samples. Among patient samples, 2 "false-negative" results were observed. These preliminary results suggest that the present protocol might be efficient for non-invasive prenatal chromosome analysis.

Fetal cells in maternal blood: the use of primed in situ (PRINS) labelling technique for fetal cell detection and sex assessment.

Prenatal diagnosis is presently performed following invasive procedures with variable risks of fetal loss; non-invasive procedures using fetal cells in maternal blood would be welcome for the early detection of fetal sex or aneuploidy. We describe a simple and rapid protocol to detect fetal cells and thus to assess fetal sex. In a first step, nucleated blood cells were separated into mononuclear and polynuclear cells using a double density gradient centrifugation. In a second step primed in situ (PRINS) labelling technique was performed to label Y-chromosomes. 15 samples were studied and correct gender assignment was made in 13/15. The number of labelled nuclei was higher in polynuclear cell phases than in mononuclear cell phases. Moreover, the polylobular aspect of labelled nuclei from polynuclear cell phases strongly suggested that they could belong to fetal polynuclear cells. The PRINS technique combines some advantages of FISH, such as visual assessment of in situ chromosome labelling and the powerful specificity and sensitivity of PCR. In association with a simple enrichment procedure it constitutes a rapid protocol for fetal cell detection, non-invasive early prenatal sex assessment, and could further be applied to detect the main viable aneuploidies.