Chromosome conformation markers of prostate cancer and lymphoma

Abstract

To provide processes for analyzing chromosomal regions and interactions related to prognosis of prostate cancer or DLBCL.SOLUTION: Provided is a process for identifying the prognosis or diagnosis of prostate cancer, by detecting chromosomal interaction to specifying the prognosis, or by detecting the presence or absence of the chromosomal interaction in an individual to specifying the diagnosis, where the presence or absence of the chromosomal interaction is detected by a method comprising the steps of: (i) cross-linking chromosome regions which have come together in a chromosome interaction; (ii) subjecting the cross-linked regions to cleavage; (iii) ligating the cross-linked and then cleaved DNA ends to form a ligated nucleic acid and (iv) detecting the presence or absence of a ligated nucleic acid corresponding to the chromosome interaction.SELECTED DRAWING: Figure 5

Description

[Technical field]

[0001] The present invention relates to disease processes. [Background technology]

[0002] The regulatory and causative aspects of the cancer disease process are complex and cannot be easily elucidated using available DNA and protein classification methods.

[0003] Diffuse large B-cell lymphoma (DLBCL) is a cancer of B cells, a type of white blood cell involved in antibody production. It is the most common type of non-Hodgkin's lymphoma among adults, with an incidence of 7-8 cases per 100,000 people per year in the United States and the United Kingdom. However, the outcome of the disease process is poorly understood.

[0004] Prostate cancer is caused by abnormal and uncontrolled growth of cells in the prostate gland. Although survival rates for prostate cancer have improved every decade, the disease is still considered largely incurable. According to the American Cancer Society, the one-year relative survival rate for all stages of prostate cancer combined is 20%, and the five-year relative survival rate is 7%. Summary of the Invention

[0005] The inventors identify subgroups of patients defined by chromosomal conformation signatures in prostate cancer, diffuse large B-cell lymphoma (DLBCL) and lymphoma.

[0006] According to the present invention, there is provided a process for detecting a chromosomal state representing a subgroup within a population, the process comprising the steps of determining whether a chromosomal interaction associated with the chromosomal state is present within a defined region of the genome; and - the chromosomal interactions are optionally identified by a method for determining which chromosomal interactions are associated with chromosomal states corresponding to subgroups of a population, the method comprising contacting a first set of nucleic acids from subgroups having different chromosomal states with a second set of index nucleic acids and allowing hybridization of complementary sequences, the nucleic acids in the first and second sets of nucleic acids representing ligated products comprising sequences from both chromosomal regions brought together in the chromosomal interaction, and a pattern of hybridization between the first and second sets of nucleic acids allows the determination of which chromosomal interactions are specific to the subgroup; and - Subgroups are associated with prognosis in prostate cancer, and chromosomal interactions are (i) in any of the regions or genes listed in Table 6; and / or (ii) corresponds to any of the chromosomal interactions represented by any of the probes shown in Table 6; and / or (iii) is present in a 4,000-base region that includes or is adjacent to (i) or (ii); or - Subgroups are associated with DLBCL prognosis and chromosomal interactions a) in any of the regions or genes listed in Table 5; and / or b) corresponds to any of the chromosomal interactions represented by any of the probes shown in Table 5; and / or c) is present in a 4,000 base region that includes or is adjacent to (a) or (b); or - Subgroups are associated with lymphoma prognosis, and chromosomal interactions are (iv) is in any of the regions or genes listed in Table 8; and / or (v) corresponds to any of the chromosomal interactions shown in Table 8; and / or (vi) is present in a 4,000-base region including (iv) or (v) or adjacent to (iv) or (v). [Brief description of the drawings]

[0007] [Figure 1] Principal component analysis (PCA) of prostate cancer studies is shown. [Diagram 2] VENN comparison of two PCA prognostic classifiers. [Diagram 3] PCA analysis of DLBCL is shown. [Figure 4] PCA for seven BTK markers (OBD RD051) in DLBCL is shown. [Diagram 5] 1 shows an example of how classification of chromosomal interactions can be performed. [Figure 6] Markers from a canine lymphoma study that can be used in the methods of the invention are shown. Figure 6 shows marker reduction. Seventy percent of the 38 samples were used as the training set (28) and were used for marker selection. The remaining 10 were used as the test set. Multiple training and test sets were used. Univariate analysis, Fisher's exact test (results in columns D and E), and multivariate analysis penalized logistic modeling (GLMNET, results in columns B and C). Markers 2-18 are lymphoma markers and 19-23 are controls. The top 11, all loops present in lymphoma, were selected for classification. [Figure 7] Dog markers relative to human genes. The table shows the top 11 dog markers mapped to the human genome (Hg38) with the closest mapping genomic region. The neighbor network is constructed using the 11 markers (dark), light colored nodes and linker proteins using the NCI database. [Figure 8] Showing dog markers versus human genes. As before, but with pathway enrichment for the network. Only the 11 dog mapping loci were used for enrichment, and ligated modes were omitted from enrichment. Light-colored nodes belong to KEGG CML pathways. [Figure 9]XGBoost11 Mark model for training set 1 and test set 1 is shown. [Figure 10] XGBoost11 Mark model for training set 2 and test set 2. [Figure 11] XGBoost11 mark model for training set 3 and test set 3 is shown. [Figure 12] 1 shows the logistic PCA of training set 1. [Figure 13] Logistic PCA of training set 1 and test set 1. The logistic PCA model was used to predict test set 1 (triangles). The darker triangles are lymphomas (labeled) from the test set and the lighter triangles are controls from the test set. The training lymphoma samples are the darker ones and the controls are the lighter ones. [Figure 14] The ROC and AUC for training set 1 and test set 1 are shown. [Figure 15] Showing patient PFS EpiSwitch™ calling and loop dynamics at NFKB1. 118 patients called ABC or GCB using the EpiSwitch™ 10 marker human model, PFS modeling using this calling and loop dynamics, no GCB with loop died, indicating that the human model works well for disease prognosis. [Figure 16] Shown is the PFS EpiSwitch™ call and loop dynamics of 118 patients with NFATC1. As before, but for NFATC1, again this shows that the prognostic human model using the marker as one of the 10 human markers is very good at classifying. [Figure 17] A three-step approach to identify, evaluate, and validate diagnostic and prognostic biomarkers for prostate cancer (PCa) is presented. [Figure 18]PCA for five markers applied to 78 samples comprising two groups: the first group is 49 known samples (24 PCa and 25 healthy controls (Cntrl)) combined with a second group of 29 samples including 24 PCa samples and 5 healthy Cntrl samples. [Figure 19] 1 shows a workflow for developing a classifier. [Figure 20] The gene clusters associated with the classifiers are shown. [Figure 21] Overlap of EpiSwitch DLBCL-CCS and Fluidigm subtype calling and ROC curves when applied to the Discovery cohort are shown. A. Subtype calling made by EpiSwitch DLBCL-CCS and Fluidigm assays for samples of known subtype. 60 of 60 samples were identified and called by both assays. B. Receiver operating curve (ROC) for DLBCL-CCS when applied to the discovery cohort. C. Kaplan-Meier analysis (by progression-free survival) of samples called ABC or GCB by DLBCL-CCS. Samples called ABC showed significantly poorer long-term survival than samples called GCB. [Figure 22] 1 shows assignment of DLBCL subtypes in type III samples by EpiSwitch and Fluidigm assays. [Figure 23] Comparison of baseline DLBCL subtype calling in type III samples using EpiSwitch and Fluidigm with long-term survival. Kaplan-Meier survival curves for 58 DLBCL patients classified by ABC, GCB, or unclassified by the Fluidigm assay (A) or EpiSwitch DLBCL-CCS (B). Fluidigm classified 15 samples as ABC, 22 as GCB, and 21 as UNC. EpiSwitch classified 34 as ABC and 24 as GCB. [Figure 24]Median survival times according to EpiSwitch and Fluidigm classification in the validation cohort are shown. [Diagram 25] 1 shows an initial assessment of possible DLBCL subtypes. [Figure 26] PCA of DLBCL patients with baseline ABC / GCB subtype calling by EpiSwitch in the discovery cohort. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0008] Aspects of the invention The present invention relates to the prognosis of prostate cancer, particularly whether prostate cancer is aggressive or indolent.This determination is carried out by classifying any of the relevant markers disclosed herein, for example in table 6, or a preferred combination of markers, or the markers of the defined specific region disclosed herein.Therefore, the present invention relates to a method for classifying patients with prostate cancer to identify whether prostate cancer is aggressive or indolent.

[0009] The present invention also relates to determining prognosis in DLBCL, in particular determining whether the prognosis is good or bad in terms of survival. This determination is carried out by classifying any of the relevant markers disclosed herein, for example in Table 5, or a preferred combination of markers, or markers in a defined specific region disclosed herein. Thus, the present invention relates to a method for classifying patients with DLBCL to identify whether the patient has a good or bad prognosis in terms of survival, for example, determining the expected incidence of disease and / or time to death.

[0010] Essentially, in the method of the present invention, subpopulations of prostate cancer or DLBCL are identified by classification of markers.Thus, the present invention relates to, for example, a panel of epigenetic markers associated with prognosis in these diseases.The present invention therefore allows for individualized treatment of patients that accurately reflects their needs.

[0011] The present invention also relates to determining the prognosis of lymphoma based on the classification of chromosomal interactions as defined by Tables 8 or 9.

[0012] Tables 5 to 7 preferably relate to determining the prognosis of humans, and Tables 8 and 9 preferably relate to determining the prognosis of dogs.

[0013] Any of the therapeutic agents, eg, drugs, mentioned herein can be administered to the individual based on the results of the method.

[0014] Marker sets are disclosed in the tables and figures. In one embodiment, at least 10 markers from any disclosed marker set find use in the present invention. In another embodiment, at least 20% of the markers from any disclosed marker set find use in the present invention.

[0015] The invention process The process of the present invention includes a classification system for detecting chromosomal interactions related to prognosis.This classification can be carried out using the EpiSwitch™ system referred to herein, which is based on the bridge region of chromosomes that are joined together in chromosomal interaction, and cleaves chromosomal DNA, and then ligates the nucleic acid present in the bridge entity to induce ligated nucleic acid that has sequences from both regions that form chromosomal interaction.The detection of this ligated nucleic acid allows the determination of the presence or absence of specific chromosomal interaction.

[0016] Chromosomal interactions can be identified using the methods described above, in which a population of first and second nucleic acids are used. These nucleic acids can also be generated using EpiSwitch™ technology.

[0017] Epigenetic interactions relevant to the present invention As used herein, the terms "epigenetic" and "chromosomal" interactions typically refer to interactions between distal regions of chromosomes, which interactions are dynamic and change, form or break depending on the state of the chromosomal regions.

[0018] In certain processes of the invention, chromosomal interactions are typically detected by first generating a ligated nucleic acid that contains sequences from both regions of the chromosome that are part of the interaction. In such processes, the regions can be crosslinked by any suitable means. In a preferred embodiment, the interaction is crosslinked using formaldehyde, but can also be crosslinked by any aldehyde, or D-biotinyl-e-aminocaproic acid-N-hydroxysuccinimide ester or digoxigenin-3-O-methylcarbonyl-e-aminocaproic acid-N-hydroxysuccinimide ester. Paraformaldehyde can crosslink DNA strands that are 4 angstroms apart. Preferably, the chromosomal interactions are on the same chromosome, optionally 2-10 angstroms apart.

[0019] Chromosomal interactions can reflect the state of a chromosomal region, for example, if it is transcribed or repressed in response to changes in physiological conditions. The chromosomal interactions specific to the subgroups defined herein have been found to be stable, thus providing a reliable means of measuring the differences between two subgroups.

[0020] Furthermore, chromosomal interactions specific for a trait (such as prognosis) usually occur earlier in a biological process compared to other epigenetic markers, such as, for example, changes in methylation or binding of histone proteins. Thus, the process of the present invention can detect early stages of a biological process. This allows for earlier intervention (e.g., treatment), which may be more effective as a result. Chromosomal interactions can also be used to evaluate changes in prognosis, since they reflect the current state of the individual. Furthermore, there is little change in relevant chromosomal interactions between individuals within the same subgroup. Since the detection of chromosomal interactions is highly informative with up to 50 different possible interactions per gene, the process of the present invention can collate 500,000 different interactions.

[0021] Preferred Marker Sets As used herein, the term "marker" or "biomarker" refers to a specific chromosomal interaction that can be detected (classified) in the present invention. Specific markers are disclosed herein, any of which can be used in the present invention. Additional sets of markers can be used, for example, in combinations or numbers disclosed herein. Specific markers disclosed in the table herein are preferred, as well as markers present in the genes and regions mentioned in the table herein. These can be classified by any suitable method, for example, PCR or probe-based methods disclosed herein, including qPCR methods. Markers are defined herein by location or by probe and / or primer sequence.

[0022] Location and causes of epigenetic interactions Epigenetic chromosomal interactions may overlap and include the regions of chromosomes that have been shown to code related genes or undescribed genes, but may equally exist in intergenic regions.It should be further noted that the inventors have found that the epigenetic interactions in all regions are equally important in determining the state of chromosomal locus.These interactions do not necessarily exist in the coding region of a particular gene located at the locus, but may exist in intergenic regions.

[0023] The chromosomal interactions detected by the present invention can be caused by changes to the underlying DNA sequence due to environmental factors, DNA methylation, non-coding antisense RNA transcripts, non-mutagenic carcinogens, histone modifications, chromatin remodeling, and specific local DNA interactions. The changes that lead to chromosomal interactions can be caused by changes to the underlying nucleic acid sequence, which do not directly affect gene products or the mode of gene expression themselves. Such changes can be, for example, intragenic and / or extragenic SNPs, gene fusions and / or intergenic DNA, microRNAs, and deletions of non-coding RNA. For example, the process described is also useful in non-coding situations, since it is known that approximately 20% of SNPs are present in non-coding regions. In one embodiment, the regions of chromosomes that come together to form an interaction are less than 5 kb, 3 kb, 1 kb, 500 base pairs, or 200 base pairs apart on the same chromosome.

[0024] The chromosomal interaction to be detected is preferably within any of the genes mentioned in Table 5. However, it may also be upstream or downstream of the gene, for example up to 50,000, up to 30,000, up to 20,000, up to 10,000, or up to 5000 bases upstream or downstream from the gene or coding sequence.

[0025] The chromosomal interaction to be detected is preferably within any of the genes mentioned in Table 6. However, it may also be upstream or downstream of the gene, for example up to 50,000, up to 30,000, up to 20,000, up to 10,000, or up to 5000 bases upstream or downstream from the gene or coding sequence.

[0026] The chromosomal interaction to be detected is preferably within any of the genes mentioned in Table 9. However, it may also be upstream or downstream of the gene, for example up to 50,000, up to 30,000, up to 20,000, up to 10,000, or up to 5000 bases upstream or downstream from the gene or coding sequence.

[0027] Subgroups, time points and individualized treatment The objective of the present invention is to determine the prognosis. This can be done at one or more defined time points, for example at least 1, 2, 5, 8 or 10 different time points. The period between the at least 1, 2, 5 or 8 time points can be at least 5, 10, 20, 50, 80 or 100 days.

[0028] As used herein, a "subgroup" preferably refers to a population subgroup (subgroup within a population), more preferably a subgroup within a population of a particular animal, such as a particular eukaryote, or mammal (e.g., human, non-human, non-human primate, or rodent, such as mouse or rat). Most preferably, a "subgroup" refers to a subgroup within a human population. The subgroup may be a subgroup of canines, such as dogs.

[0029] The present invention includes detecting and treating specific subgroups within a population. The inventors have discovered that chromosomal interactions are different between subsets (e.g., at least two subsets) within a given population. By identifying these differences, physicians can classify patients as part of one subset of the population, as described in the process. Thus, the present invention provides physicians with a process for personalizing medicines for patients based on epigenetic chromosomal interactions.

[0030] In one aspect, the present invention relates to a method for treating a chronic inflammatory bowel disease, comprising administering to an individual: are fast or slow "developers" and / or have aggressive or indolent disease, This relates to testing whether or not

[0031] The present invention can also determine the expected survival time of an individual.

[0032] Such testing may be used to select subsequent treatment methods for the patient, for example, type of drug and / or its dose and / or its frequency of administration.

[0033] Generation of ligated nucleic acids Certain aspects of the present invention utilize ligated nucleic acid, particularly ligated DNA. These contain sequences from both regions that come together in chromosomal interactions, thus providing information about the interaction. The EpiSwitch™ method described herein uses the generation of such ligated nucleic acid to detect chromosomal interactions.

[0034] Thus, the process of the invention can include generating ligated nucleic acid (e.g., DNA) by the following steps (including methods comprising these steps): (i) cross-linking epigenetic chromosomal interactions present at chromosomal loci, preferably in vitro; (ii) optionally isolating the cross-linked DNA from said chromosomal locus; (iii) subjecting the cross-linked DNA to cleavage, for example by restriction digestion, with an enzyme that cuts at least once, particularly an enzyme that cuts at least once within the chromosomal locus; (iv) ligating the cross-linked, cleaved DNA ends (particularly to form a DNA loop); and (v) Optionally, identifying the presence of the ligated DNA and / or the DNA loop, in particular using techniques such as PCR (polymerase chain reaction), to identify the presence of specific chromosomal interactions.

[0035] These steps can be performed to detect chromosomal interactions according to any of the embodiments referred to herein. These steps can also be performed to generate the first and / or second set of nucleic acids referred to herein.

[0036] PCR (polymerase chain reaction) can be used to detect or identify the ligated nucleic acid, for example, to identify the state of the locus, since the size of the PCR product generated can indicate the specific chromosomal interaction present. In a preferred embodiment, at least one, two, or three primers or primer pairs shown in Table 5 are used in the PCR reaction. In other embodiments, at least 1, 10, 20, 30, 50, or 80 primers or primer pairs, or primers or primer pairs shown in Table 6 are used in the PCR reaction. Those skilled in the art will recognize a number of restriction enzymes that can be used to cut DNA within a chromosomal locus of interest. It will be clear that the particular enzyme used will depend on the locus studied and the sequence of the DNA located therein. A non-limiting example of a restriction enzyme that can be used to cut DNA as described in the present invention is TaqI.

[0037] EpiSwitch™ Technology EpiSwitch™ technology also relates to the use of microarray EpiSwitch™ marker data in the detection of phenotype-specific epigenetic chromosome conformation signatures.Embodiments such as EpiSwitch™ that utilize ligated nucleic acids in the manner described herein have several advantages.They have low levels of stochastic noise, for example, because the nucleic acid sequence from a first set of nucleic acids of the present invention hybridizes or does not hybridize with a second set of nucleic acids.This provides a binary result that allows a relatively simple way to measure complex mechanisms at the epigenetic level.EpiSwitch™ technology also has a fast processing time and low cost.In one embodiment, the processing time is 3 to 6 hours.

[0038] Samples and sample processing The process of the invention is usually carried out on a sample. The sample can be obtained at a defined time point, for example at any time point defined herein. The sample usually contains DNA from an individual. It usually contains cells. In one embodiment, the sample is obtained by minimally invasive means, for example, it can be a blood sample. The DNA can be extracted and cut with standard restriction enzymes. This allows to predetermine which chromosome conformations are retained and detected by the EpiSwitch™ platform. Due to the synchronization of chromosomal interactions between tissues and blood, including horizontal transfer, blood samples can be used to detect chromosomal interactions in tissues, such as tissues associated with disease. For certain diseases, such as cancer, there is an advantage to using blood, since genetic noise due to mutations can affect the "signal" of chromosomal interactions in associated tissues.

[0039] Properties of the Nucleic Acids of the Invention The present invention relates to certain nucleic acids, such as the ligated nucleic acids described herein as used or produced in the process of the invention. These may be the same as the first and second nucleic acids referred to herein, or may have any of their properties. The nucleic acids of the present invention typically comprise two portions, each of which comprises a sequence from one of the two regions of a chromosome that come together in a chromosomal interaction. Typically, each portion is at least 8, 10, 15, 20, 30 or 40 nucleotides in length, for example 10-40 nucleotides in length. Preferred nucleic acids comprise a sequence from any of the genes referred to in any of the tables. Typically, preferred nucleic acids comprise the specific probe sequences referred to in Table 5, or comprise fragments and / or homologues of such sequences. Typically, preferred nucleic acids may comprise the specific probe sequences referred to in Table 6, or comprise fragments and / or homologues of such sequences.

[0040] Preferably, the nucleic acid is DNA. When a specific sequence is provided, it is understood that the present invention can use a complementary sequence as required in a particular embodiment. Preferably, the nucleic acid is DNA. When a specific sequence is provided, it is understood that the present invention can use a complementary sequence as required in a particular embodiment.

[0041] Primers shown in Table 5 may also be used in the present invention as described herein. In one embodiment, primers comprising any of the sequences shown in Table 5, or fragments and / or homologues of any of the sequences shown in Table 5 are used. Primers shown in Table 6 may also be used in the present invention as referred to herein. In one embodiment, primers comprising any of the sequences shown in Table 6, or fragments and / or homologues of any of the sequences shown in Table 6 are used. Primers shown in Table 8 may also be used in the present invention as described herein. In one embodiment, primers comprising any of the sequences shown in Table 8, or fragments and / or homologues of any of the sequences shown in Table 8 are used.

[0042] A second set of nucleic acids - "index" sequences The second set of nucleic acid sequences has the function of being a set of index sequences, and is essentially a set of nucleic acid sequences suitable for identifying specific sequences of a subgroup. They can represent "background" chromosomal interactions, which may or may not be selected in some way. They are generally a subset of all possible chromosomal interactions.

[0043] The second set of nucleic acids may be derived by any suitable process. They may be derived computationally or based on individual chromosomal interactions. They typically represent a larger population than the first set of nucleic acids. In one particular embodiment, the second set of nucleic acids represents all possible epigenetic chromosomal interactions in a specific set of genes. In another particular embodiment, the second set of nucleic acids represents a majority of all possible epigenetic chromosomal interactions present in the population described herein. In one particular embodiment, the second set of nucleic acids represents at least 50% or at least 80% of epigenetic chromosomal interactions in at least 20, 50, 100 or 500 genes, such as 20-100 or 50-500 genes.

[0044] The second set of nucleic acids typically represents at least 100 possible epigenetic chromosomal interactions that alter, regulate, or in some way mediate phenotypes in a population. The second set of nucleic acids may represent chromosomal interactions that affect disease states (typically associated with diagnosis or prognosis) in a species. The second set of nucleic acids typically includes sequences that represent epigenetic interactions that are both associated and not associated with prognostic subgroups.

[0045] In one particular embodiment, the second set of nucleic acids is at least partially derived from sequences of natural origin in the population, typically obtained by an in silico process. The nucleic acids may further comprise single or multiple mutations compared to the corresponding portion of the nucleic acid present in the nucleic acid of natural origin. Mutations include deletions, substitutions, and / or additions of one or more nucleotide base pairs. In one particular embodiment, the second set of nucleic acids may comprise sequences representing homologs and / or orthologs with at least 70% sequence identity to the corresponding portion of the nucleic acid present in the species of natural origin. In another particular embodiment, at least 80% sequence identity or at least 90% sequence identity to the corresponding portion of the nucleic acid present in the species of natural origin is provided.

[0046] A second set of properties of the nucleic acids In one particular embodiment, the second set of nucleic acids has at least 100 different nucleic acid sequences, preferably at least 1000, 2000 or 5000 different nucleic acid sequences, with up to 100,000, 1,000,000 or 10,000,000 different nucleic acid sequences. A typical number is between 100 and 1,000,000, such as between 1,000 and 100,000 different nucleic acid sequences. All or at least 90% or at least 50% of these correspond to different chromosomal interactions.

[0047] In one particular embodiment, the second set of nucleic acids represents chromosomal interactions at at least 20 different loci or genes, preferably at least 40 different loci or genes, more preferably at least 100, at least 500, at least 1000 or at least 5000 different loci or genes, such as 100 to 10,000 different loci or genes. The length of the nucleic acids of the second set is suitable to specifically hybridize to the nucleic acids of the first set according to Watson-Crick base pairs to allow identification of subgroup-specific chromosomal interactions. Typically, the second set of nucleic acids comprises two portions that correspond in sequence to two chromosomal regions that come together in a chromosomal interaction. The second set of nucleic acids typically comprises a nucleic acid sequence that is at least 10 bases (nucleotides) long, preferably 20 bases long, preferably also 30 bases long. In another embodiment, the nucleic acid sequence can be up to 500 base pairs in length, preferably up to 100 base pairs, preferably also up to 50 base pairs. In a preferred embodiment, the second set of nucleic acids comprises nucleic acid sequences between 17 base pairs and 25 base pairs. In one embodiment, at least 100%, 80% or 50% of the second set of nucleic acid sequences have a length as described above. Preferably, the different nucleic acids do not have overlapping sequences, e.g., at least 100%, 90%, 80% or 50% of the nucleic acids do not have the same sequence over at least 5 consecutive nucleotides.

[0048] The same set of second nucleic acids can be used with different sets of first nucleic acids representing subgroups for different traits, considering the second set of nucleic acids acting as an "index," i.e., the second set of nucleic acids can represent a "universal" collection of nucleic acids that can be used to identify chromosomal interactions associated with different traits.

[0049] A first set of nucleic acids The first set of nucleic acids is typically from a subgroup associated with prognosis. The first nucleic acid may have any of the characteristics and properties of the second set of nucleic acids mentioned herein. The first set of nucleic acids is usually derived from a sample from an individual that has undergone the treatment and processing described herein, particularly the crosslinking and cleavage steps of EpiSwitch™. Typically, the first set of nucleic acids represents all or at least 80% or 50% of the chromosomal interactions present in the sample taken from the individual.

[0050] Typically, the first set of nucleic acids represents a smaller population of chromosomal interactions across the loci or genes represented by the second set of nucleic acids compared to the chromosomal interactions represented by the second set of nucleic acids, i.e., the second set of nucleic acids represents a background or index set of interactions at a defined set of loci or genes.

[0051] Nucleic acid library Any type of nucleic acid population referred to herein may be in the form of a library containing at least 200, at least 500, at least 1000, at least 5000, or at least 10000 different nucleic acids of that type, such as a "first" or "second" nucleic acid. Such libraries may be in the form of arrays bound to them. The libraries may contain some or all of the probes or primer pairs shown in Tables 5 or 6. The libraries may contain all of the probe sequences from any of the tables disclosed herein.

[0052] Hybridization The present invention requires a means for allowing the fully or partially complementary nucleic acid sequences from the first set of nucleic acids and the second set of nucleic acids to hybridize.In one embodiment, all of the first set of nucleic acids are contacted with all of the second set of nucleic acids in a single assay, i.e., in a single hybridization step.However, any suitable assay can be used.

[0053] Labeled Nucleic Acids and Hybridization Patterns The nucleic acid referred to herein can preferably be labeled using an independent label, such as a fluorophore (fluorescent molecule) or a radioactive label, which aids in the detection of successful hybridization.Certain labels can be detected under UV light.The pattern of hybridization, for example on the array described herein, represents the difference in epigenetic chromosomal interactions between two subgroups, thus providing a process of comparing epigenetic chromosomal interactions and determining which epigenetic chromosomal interactions are specific to subgroups in the population of the present invention.

[0054] The term "pattern of hybridization" broadly covers the presence or absence of hybridization between a first set of nucleic acids and a second set of nucleic acids, i.e., which specific nucleic acids from the first set hybridize to which specific nucleic acids from the second set, and thus is not limited to any particular assay or technique, or the need to have a surface or array capable of detecting the "pattern."

[0055] Selection of subgroups with specific characteristics The present invention provides a process comprising detecting the presence or absence of a chromosomal interaction, typically 5-20 or 5-500 such interactions, preferably 20-300 or 50-100 interactions, to determine the presence or absence of a feature associated with the prognosis of an individual. Preferably, the chromosomal interaction is in any of the genes mentioned herein. In one embodiment, the chromosomal interactions being classified are those represented by the nucleic acids in Table 5. In another embodiment, the chromosomal interactions are those represented in Table 6. In a further embodiment, the chromosomal interactions being classified are those represented by the nucleic acids in Table 8. The column in the table entitled "Loops Detected" indicates which subgroups are detected by each probe. Detection can be the presence or absence of a chromosomal interaction within that subgroup, as indicated by "1" and "-1".

[0056] Individuals being tested Examples of species from which the individual to be tested originate are mentioned herein. In addition, the individual to be tested in the process of the present invention may be selected in any manner. The individual may be susceptible to any disease mentioned herein and / or in need of any treatment mentioned herein. The individual may be undergoing any treatment mentioned herein. In particular, the individual suffers from or is suspected of suffering from prostate cancer or DLBCL. The individual suffers from or is suspected of suffering from lymphoma.

[0057] Favorable gene regions, loci, genes and chromosomal interactions for prostate cancer For all aspects of the invention, preferred gene regions, loci, genes, and chromosomal interactions are referred to in tables, e.g., Table 6. Typically, in the process of the invention, chromosomal interactions are detected from at least 1, 2, 3, 4, or 5 of the relevant genes listed in Table 6. Preferably, the presence or absence of at least 1, 2, 3, 4, or 5 of the relevant specific chromosomal interactions represented by the probe sequences in Table 6 is detected. The chromosomal interactions may be upstream or downstream of any of the genes referred to herein, e.g., 50 kb upstream or 20 kb downstream from the coding sequence.

[0058] For all aspects of the invention, preferred gene regions, loci, genes, and chromosomal interactions are referred to in Table 25. Typically, in the process of the invention, chromosomal interactions are detected from at least 2, 4, 8, 10, 14, or all of the relevant genes listed in Table 25. Preferably, the presence or absence of at least 2, 4, 8, 10, 14, or all of the relevant specific chromosomal interactions shown in Table 25 are detected. The chromosomal interaction may be upstream or downstream of any of the genes referred to herein, for example, 50 kb upstream or 20 kb downstream from the coding sequence.

[0059] In one embodiment, a specific combination of markers disclosed herein and represented (identified) by the following combination of genes is classified: ETS1, MAP3K14, SLC22A3 and CASP2. This may be for determining a diagnosis. Preferably, at least two or three of these markers are classified.

[0060] In another embodiment, the specific marker combination disclosed herein is classified, which is represented (identified) by the following gene combinations: BMP6, ERG, MSR1, MUC1, ACAT1 and DAPK1. This may be for determining prognosis (high risk category 3 vs. low risk category 1 by nested PCR markers). Preferably, at least two or three of these markers are classified.

[0061] In a further embodiment, the specific marker combinations disclosed herein are classified as represented (identified) by the following gene combinations: HSD3B2, VEGFC, APAF1, MUC1, ACAT1 and DAPK1. This may be for determining prognosis (high risk category 3 vs. intermediate risk category 2). Preferably, at least two or three of these markers are classified.

[0062] Favorable gene regions, loci, genes, and chromosomal interactions for DLBCL Typically, at least 10, 20, 30, 50 or 80 chromosomal interactions are classified from any of the genes or regions disclosed in a table herein or part of a table disclosed herein. Preferably, at least 10, 20, 30, 50 or 80 chromosomal interactions are classified from any of the genes or regions disclosed in Table 5.

[0063] Preferably, at least 2, 3, 5, 8 of the markers in Table 7 are classified.

[0064] Preferably, the presence or absence of at least 10, 20, 30, 50 or 80 chromosomal interactions represented by the probe sequences in Table 5 is detected. The chromosomal interactions may be upstream or downstream of any of the genes mentioned herein, e.g., 50 kb upstream or 20 kb downstream from the coding sequence.

[0065] Preferably, at least one, two, five, eight or all of the first ten markers shown in Table 5 are classified. In one embodiment, at least one, two, three or six markers from Table 5 are classified, each of which corresponds to a different gene selected from STAT3, TNFRSF13B, ANXA11, MAP3K7, MEF2B and IFNAR1.

[0066] Favorable gene regions, loci, genes and chromosomal interactions for lymphoma Typically, at least 10, 20, 30 or 50 chromosomal interactions are classified from any of the genes or regions disclosed in a table herein or part of a table disclosed herein. Preferably, at least 10, 20, 30 or 50 chromosomal interactions are classified from any of the genes or regions disclosed in Table 8.

[0067] Preferably, at least 5, 10 or 15 of the markers in Table 9 are classified.

[0068] The chromosomal interaction may be upstream or downstream of any of the genes mentioned herein, for example 50 kb upstream or 20 kb downstream from the coding sequence.

[0069] In one embodiment, at least one of the first eleven markers shown in Figure 6 is classified. In another embodiment, at least one, two, three or six markers from Table 8 are classified, each of which corresponds to a different gene selected from the following: STAT3, TNFRSF13B, ANXA11, MAP3K7, MEF2B and IFNAR1.

[0070] Types of chromosomal interactions In one embodiment, the locus (including the gene and / or location where the chromosomal interaction is detected) may contain a CTCF binding site. This is any sequence that can bind to the transcriptional repressor CTCF. The sequence may consist of or include the sequence CCCTC, which may be present in one, two, or three copies at the locus. The CTCF binding site sequence may include the sequence CCGCGNGGNGGCAG (IUPAC notation). The CTCF binding site may be within at least 100, 500, 1000, or 4000 bases of the chromosomal interaction, or within any of the chromosomal regions shown in Tables 5 or 6. The CTCF binding site may be within at least 100, 500, 1000, or 4000 bases of the chromosomal interaction, or within any of the chromosomal regions shown in Tables 5 or 6.

[0071] In one embodiment, the chromosomal interaction detected is in any of the gene regions shown in Tables 5 or 6. If ligated nucleic acids are detected in the process, sequences shown in any of the probe sequences in Tables 5 or 6 can be detected.

[0072] Thus, typically sequences from both regions of the probe (i.e., from both sites of chromosomal interaction) can be detected. In preferred embodiments, a probe that contains or consists of the same or complementary sequence as any of the probes shown in the table is used in the process. In some embodiments, a probe that contains a sequence that is homologous to any of the probe sequences shown in the table is used.

[0073] Tables provided herein Tables 5 and 6 show probe (EpiSwitch™ marker) data and gene data representing chromosomal interactions associated with prognosis. The probe sequences indicate sequences that can be used to detect ligated products generated from both sites of the gene regions brought together in the chromosomal interaction, i.e., the probes contain sequences complementary to sequences in the ligated products. The first two sets of start-end positions indicate the probe positions and the second two sets of start-end positions indicate the relevant 4kb regions. The following information is provided in the tables of probe data: - HyperG_Stats: p-value for the probability of finding that number of significant EpiSwitch™ markers within a locus based on the hypergeometric enrichment parameters - Total Probes: The total number of EpiSwitch™ conformations tested at the locus - Number of significant probes: the number of EpiSwitch™ conformations found to be statistically significant at the locus - FDR HyperG: Hypergeometric p-value corrected for multiple trials (Fimmunoresposivenesse Discovery Rate) - Percent Significance: The percentage of significant EpiSwitch™ markers relative to the number of markers tested at the locus - logFC: base 2 logarithm of epigenetic ratio (FC) - AveExpr: Average log2 expression for the probes across all arrays and channels - T: Moderated t-statistic - p-value: raw p-value - adj.p-value: The adjusted p-value or q-value - The BB statistic (lod or B) is the log odds that the gene is differentially expressed. - FC - Non-logarithmic change - FC_1 - Non-log scale factor change around zero - LS-Binary value, which is related to the FC_1 value. For FC_1 values ​​less than -1.1, it is set to -1, for FC_1 values ​​greater than 1.1, it is set to 1. Values ​​in between these values ​​are 0.

[0074] Tables 5 and 6 show genes where relevant chromosomal interactions were found to occur. The p-values ​​in the locus tables are the same as HyperG Stats (p-value for the probability of finding that number of significant EpiSwitch™ within the locus based on the hypergeometric enrichment parameters). The LS column shows the presence or absence of relevant interactions (prognostic status) with that particular subgroup.

[0075] With respect to Table 5, DLBCL refers to the prognostic marker designated as 1 and healthy refers to the healthy control designated as -1.

[0076] The probe is designed to be 30 bp away from the Taq1 site. In the case of PCR, PCR primers are typically designed to detect the ligated product, but at a variable position from the Taq1 site.

[0077] Probe Location: Start1 - 30 bases upstream of the TaqI site on fragment 1 End1 - TaqI restriction site on fragment 1 TaqI restriction site on Start2-fragment 2 End2 - 30 bases downstream of the TaqI site on fragment 2

[0078] 4kb sequence location: Start1 - 4000 bases upstream of the TaqI site on fragment 1 End1 - TaqI restriction site on fragment 1 TaqI restriction site on Start2-fragment 2 End2 - 4000 bases downstream of the TaqI site on fragment 2

[0079] GLMNET values ​​related to the procedure of fitting the entire Lasso or Elastic Net regularization (lambda set to 0.5 (Elastic Net)).

[0080] In the tables herein, the aggressive subgroup of prostate cancer refers to class 3 patients with the following description: - PSA level greater than 20ng / ml, - A Gleason score between 8 and 10, and - T stage is T2c, T3, or T4.

[0081] In the tables herein, the indolent subgroup of prostate cancer refers to Class 1 patients with the following description: - PSA level less than 10 ng per ml, - Gleason score ≤ 6, and - T stage is between T1 and T2a.

[0082] Table 7 shows preferred markers for DLBCL. Tables 8 and 9 show preferred markers for lymphoma.

[0083] Tables 5-7 are preferred for classification of humans. Tables 8 and 9 are preferred for classification of canines, e.g. dogs.

[0084] Approaches taken to identify markers and panels of markers The invention described herein relates to chromosome conformation profiles that are closely linked to phenotypes and 3D architectures as a modality of regulation in itself. Biomarker discovery was based on annotation by pattern recognition and screening in representative cohorts of clinical samples that represent phenotypic differences. Significant portions of genomes were annotated and screened across coding and non-coding parts, and across large swings in non-coding 5'' and 3'' of known genes for the identification of statistically consistent conditional seeding chromosome conformations that anchor in non-coding sites within (introns) or outside of open reading frames.

[0085] The selection of the best markers is driven by statistical data and p-values ​​for the marker leads. Reference to a specific gene is used to facilitate positional reference, and the nearest gene is usually used as a reference. It is not possible to exclude the possibility that the cis-position from the gene and the associated nearby chromosome conformation may contribute a specific element of regulation to the expression of that particular gene. At the time of marker selection or validation, no expression parameters are required for the gene referenced as position coordinates in the name of the chromosome conformation. The selected and validated chromosome conformation in the signature is seeding a stratification entity by itself, regardless of the expression profile of the gene used in the reference. Further research may be performed on associated regulatory modalities, such as SNPs at anchor sites, changes in gene transcription profiles, changes at the level of H3K27ac, etc.

[0086] We address the problem of clinical phenotypic differences and their stratification from, for example, a regulatory network framework, based on fundamental biology and epigenetic control over the phenotype. Thus, to aid stratification, network changes can be captured, which is preferably done via a signature of several biomarkers, for example, by following a machine learning algorithm for marker reduction, including the evaluation of the optimal number of markers to stratify the test cohort with minimal noise. This usually ends up with 3-17 markers on a case-by-case basis. The selection of markers for the panel may be done by statistical performance of cross-validation (e.g., not by the functional relevance of neighboring genes used for reference names).

[0087] The panel of markers (with the names of the flanking genes) is the product of a clustered selection from a screen across a significant part of the genome in an unbiased manner that analyzes statistical seeding power across 14,000-60,000 annotated EpiSwitch sites across a significant part of the genome. This should not be perceived as a coordinated capture of chromosomal conformations on genes of known functional value for stratification problems. The total number of sites of chromosomal interactions is 1.2 million, so the potential number of combinations is 1.2 million to the power of 1.2 million. Nevertheless, the approach followed allows the identification of relevant chromosomal interactions.

[0088] The specific markers provided by this application are selected to be statistically (significantly) associated with the disease, as demonstrated by the p-values ​​in the associated table. Each marker can be seen as representing a biological epigenetic event as part of the deregulation of the network that appears in the associated disease. In practice, this means that these markers are prevalent in a group of patients compared to controls. By way of example, on average, an individual marker may typically be present in 80% of the patients tested and 10% of the controls tested.

[0089] The simple addition of all markers does not represent the network interrelationships between some of the deregulations. Here, a standard multivariate biomarker analysis GLMNET (R package) is introduced. The GLMNET package helps to identify the interdependencies between some of the markers, which reflect their joint role in achieving the deregulation that leads to the disease phenotype. Modeling and then testing the markers with the highest GLMNET score leads to the identification of not only the minimum number of markers that correctly identify the patient cohort, but also the minimum number that leads to the fewest false positives in the patient control group, due to the background statistical noise of the low prevalence in the control group. Typically, a group (combination) of selected markers (e.g., 3 to 10) provides an optimal balance between both sensitivity and specificity of detection, which arises in conjunction with a multivariate analysis from the individual characteristics of all statistically significant markers selected for the disease.

[0090] The table herein gives the reference name of the array probe (60mer) for array analysis overlapping the junction between the long-range interaction sites, the chromosome number, and the start and end of the two chromosomal fragments to be juxtaposed. The table also gives the standard array readout in competitive hybridization of disease against control samples (labeled with two different fluorescent colors) for each of the markers. As standard readout, indicated for each marker probe: - Mean expression signal - t-test for significant differences in fluorescent detection between control and disease samples - p-value of significance of marker readout - adjusted p-value (using Bonferroni correction for large datasets), B - background signal, FC - fold change for color detection in control samples -FC_1 - fold change for second color detection in case (disease or disease type) samples, LS (loop status) - typical fluorescent signal between the thresholds of the two colors in competitive hybridization, -1 means that the signal is typical for patient samples of the corresponding fluorescent color when tested against the probes of the CGH array. - Current locus - Total probe count - the number of different positional probes on the array that were tested across that locus - Number of significant probes - the number that were found to be significant in distinguishing case and control samples - Hypergeometric Stat is a statistic of enrichment of loci with significant probes for disease detection - FDR HyperG is the same statistic adjusted for large datasets by FDR (standard procedure) - The percentage of probes found to be significant at that locus - logFC is the logarithm of the fold change of the array readout for that probe. Focusing on loci with high enrichment of significant probes aids in the selection of top probes representing regulatory hubs with multiple inputs associated with, for example, disease providing markers with the best coverage of network deregulation.

[0091] Preferred Embodiments for Sample Preparation and Detection of Chromosomal Interactions Methods for preparing samples and detecting chromosome conformation are described herein. Optimized (non-conventional) versions of these methods can be used, for example, as described in this section.

[0092] Typically, a sample comprises at least 2×10 5 cells. A sample may comprise up to 5×10 5 cells. In one embodiment, a sample comprises between 2×10 5 and 5.5×10 5 cells.

[0093] Cross-linking of epigenetic chromosomal interactions present at chromosomal loci is described herein. This can be performed before cell lysis occurs. Cell lysis can be performed for 3-7 minutes, such as 4-6 minutes or about 5 minutes. In some embodiments, cell lysis is performed for at least 5 minutes and less than 10 minutes.

[0094] Digestion of DNA with restriction enzymes is described herein. Typically, DNA restriction is carried out at about 55° C. to about 70° C., such as about 65° C., for a period of about 10 to 30 minutes, such as about 20 minutes.

[0095] Preferably, frequent-cutter restriction enzymes are used, resulting in ligated DNA fragments having an average fragment size of up to 4000 base pairs. Optionally, restriction enzymes result in ligated DNA fragments having an average fragment size of about 200-300 base pairs, such as about 256 base pairs. In one embodiment, typical fragment sizes are 200-4000 base pairs, such as 400-2,000 or 500-1,000 base pairs.

[0096] In one embodiment of the EpiSwitch Method, no DNA precipitation step is performed between the DNA restriction digestion step and the DNA ligation step.

[0097] DNA ligation is described herein. Typically, DNA ligation is carried out for 5-30 minutes, such as about 10 minutes.

[0098] Proteins in the sample may be enzymatically digested, for example, using a proteinase, optionally proteinase K. Proteins may be enzymatically digested for about 30 minutes to 1 hour, for example, about 45 minutes. In one embodiment after protein digestion, for example proteinase K digestion, there is no crosslink reversal or phenol DNA extraction step.

[0099] In one embodiment, PCR detection can detect single copies of the ligated nucleic acid, preferably with a binary readout for the presence / absence of the ligated nucleic acid.

[0100] FIG. 5 shows a preferred method for detecting chromosomal interactions.

[0101] Processes and Uses of the Invention The process of the present invention can be described in various ways. It can be described as a method of making a ligated nucleic acid, which comprises the steps of (i) cross-linking in vitro chromosomal regions brought together by chromosomal interactions, (ii) exposing the cross-linked DNA to cleavage or restriction digestion cleavage, and (iii) ligating the cross-linked and cleaved DNA ends to form a ligated nucleic acid, the detection of which can be used to determine the chromosomal status at a locus, preferably - the locus may be any of the loci, regions, or genes mentioned in Table 5, and / or - the chromosomal interaction may be any of the chromosomal interactions corresponding to any of the probes mentioned herein or disclosed in Table 5, and / or - The ligated product may have or contain (i) a sequence that is identical or homologous to any of the probe sequences disclosed in Table 5 or (ii) a sequence that is complementary to (ii).

[0102] The process of the present invention can be described as a process for detecting chromosomal states representing different subgroups within a population, the process comprising determining whether chromosomal interactions are present or absent within defined epigenetically active regions of the genome, preferably comprising: - Subgroups are defined by the presence or absence of prognosis, and / or - the chromosomal status may be any locus, region, or gene mentioned in Table 5, and / or - the chromosomal interactions may be any of those mentioned in Table 5 or those corresponding to any of the probes disclosed in Table 5.

[0103] The process of the present invention can be described as a method of making a ligated nucleic acid, the method comprising the steps of: (i) cross-linking in vitro chromosomal regions brought together by chromosomal interactions; (ii) subjecting the cross-linked DNA to cleavage or restriction digestion cleavage; and (iii) ligating the cross-linked and cleaved DNA ends to form a ligated nucleic acid, the detection of which can be used to determine the chromosomal status at a locus, preferably - the locus may be any of the loci, regions, or genes mentioned in Table 6, and / or - the chromosomal interaction may be any of the chromosomal interactions corresponding to any of the probes mentioned herein or disclosed in Table 6, and / or - The ligated product may have or contain (i) a sequence that is identical or homologous to any of the probe sequences disclosed in Table 6 or (ii) a sequence that is complementary to (ii).

[0104] The process of the present invention can be described as a process for detecting chromosomal states representing different subgroups within a population, the process comprising determining whether chromosomal interactions are present or absent within defined epigenetically active regions of the genome, preferably comprising: - Subgroups are defined by the presence or absence of prognosis, and / or - the chromosomal status may be any locus, region, or gene mentioned in Table 6, and / or - the chromosomal interactions may be any of those mentioned in Table 6 or those corresponding to any of the probes disclosed in Table 6.

[0105] The invention includes the detection of chromosomal interactions at any locus, gene or region mentioned in Table 5. The invention includes the use of the nucleic acids and probes mentioned herein to detect chromosomal interactions, for example, the use of at least 1, 5, 10, 20 or 50 such nucleic acids or probes to detect chromosomal interactions. The nucleic acids or probes preferably detect chromosomal interactions at at least 1, 5, 10, 20 or 50 different loci or genes. The invention includes the detection of chromosomal interactions using any of the primers or primer pairs listed in Table 5, or using variants of these primers (sequences that include the primer sequences or that include fragments and / or homologues of the primer sequences) described herein.

[0106] The present invention includes the detection of chromosomal interactions at any locus, gene or region mentioned in Table 6. The present invention includes the use of the nucleic acids and probes mentioned herein to detect chromosomal interactions. The present invention includes the detection of chromosomal interactions using any of the primers or primer pairs listed in Table 6, or using variants of these primers described herein (sequences that include the primer sequences or that include fragments and / or homologues of the primer sequences).

[0107] When analyzing whether a chromosomal interaction occurs within a defined gene, region, or location, both portions of the chromosomes that come together in the interaction are within the defined gene, region, or location, or in some embodiments, only one portion of the chromosome is within the defined gene, region, or location.

[0108] Similarly, the chromosomal interactions of Tables 8 and 9 may be used in the processes and methods of the present invention.

[0109] Use of the methods of the present invention to identify new treatments Knowledge of chromosomal interactions can be used to identify new treatments for disease. The present invention provides methods and uses of the chromosomal interactions defined herein to identify or design new therapeutic agents, for example related to the treatment of prostate cancer or DLBCL.

[0110] homologue A homolog of a polynucleotide / nucleic acid (e.g., DNA) sequence is referred to herein. Such a homolog typically has at least 70% homology, preferably at least 80%, at least 85%, at least 90%, at least 95%, at least 97%, at least 98%, or at least 99% homology, for example, over a region of at least 10, 15, 20, 30, 100 or more consecutive nucleotides, or over a portion of a nucleic acid derived from a region of a chromosome involved in chromosomal interaction. Homology can be calculated based on nucleotide identity (sometimes referred to as "hard homology").

[0111] Thus, in certain embodiments, homologs of polynucleotide / nucleic acid (e.g., DNA) sequences are referred to herein by reference to percentage sequence identity. Typically, such homologs have at least 70% sequence identity, preferably at least 80%, at least 85%, at least 90%, at least 95%, at least 97%, at least 98%, or at least 99% sequence identity, for example, over a region of at least 10, 15, 20, 30, 100 or more contiguous nucleotides, or over a portion of the nucleic acid derived from a region of a chromosome involved in a chromosomal interaction.

[0112] For example, the UWGCG package provides the BESTFIT program, which can be used to calculate homology and / or % sequence identity (e.g., used with its default settings) (Devereux et al (1984) Nucleic Acids Research 12, p387-395). The PILEUP and BLAST algorithms can be used to calculate homology and / or % sequence identity and / or align sequences (such as to identify equivalent or corresponding sequences (typically with their default settings)), as described, for example, in Altschul SF (1993) J Mol Evol 36: 290-300; Altschul, S, F et al (1990) J Mol Biol 215: 403-10.

[0113] Software for performing BLAST analysis is publicly available through the National Center for Biotechnology Information. The algorithm involves first identifying high-scoring sequence pairs (HSPs) by identifying short words of length W in the matched sequences that match or meet some positive threshold score T when aligned with words of the same length in the database sequences. T is called the neighborhood word score threshold (Altschul et al, supra). These initial neighborhood word hits act as seeds to initiate searches to find HSPs that contain them. The word hits are extended in both directions along each sequence for as long as the cumulative alignment score can be increased. The extension of the word hits in each direction is stopped when: the cumulative alignment score falls by a quantity X from its maximum achieved value; the accumulation of one or more negatively scoring residue alignments causes the cumulative alignment score to fall below zero; or the end of either sequence is reached. The BLAST algorithm parameters W5 T and X determine the sensitivity and speed of the alignment. The BLAST program uses as defaults a wordlength (W) of 11, the BLOSUM62 scoring matrix (see Henikoff and Henikoff (1992) Proc. Natl. Acad. Sci. USA 89:10915-10919) of alignments (B) of 50, expectation (E) of 10, M=5, N=4, and a comparison of both strands.

[0114] The BLAST algorithm performs a statistical analysis of the similarity between two sequences (see, e.g., Karlin and Altschul (1993) Proc. Natl. Acad. Sci. USA 90: 5873-5787). One measure of similarity provided by the BLAST algorithm is the minimum sum probability (P(N)), which provides an indication of the probability that a match between two polynucleotide sequences would occur by chance. For example, a sequence is considered to be similar to another sequence if the minimum sum probability in the comparison of the first sequence to the second sequence is less than about 1, preferably less than about 0.1, more preferably less than about 0.01, and most preferably less than about 0.001.

[0115] Homologous sequences typically differ by 1, 2, 3, 4 or more bases, such as less than 10, 15, or 20 bases (which may be nucleotide substitutions, deletions, or insertions). These changes may be measured across any of the above-mentioned regions in conjunction with calculating homology and / or % sequence identity.

[0116] The homology of a "primer pair" can be calculated, for example, by considering the two sequences as a single sequence (as if the two sequences were joined together) for purposes of subsequent comparison to another primer pair (which is also considered as a single sequence).

[0117] array The second set of nucleic acids can be bound to an array, in one embodiment, there are at least 15,000, 45,000, 100,000 or 250,000 different second nucleic acids bound to the array, which preferably represent at least 300, 900, 2000 or 5000 loci. In one embodiment, one or more or all different populations of second nucleic acids are bound to more than one separate region of the array, effectively repeating on the array, allowing for error detection. The array can be based on the Agilent SurePrint G3 Custom CGH microarray platform. Detection of binding of the first nucleic acid to the array can be performed by a dual color system.

[0118] Therapeutic agents (e.g., selected based on classification of the individual or selected based on a test according to the invention) Therapeutic agents are referred to herein. The invention provides such agents for use in preventing or treating a condition in a particular individual, such as an individual identified by the processes of the invention. This may involve administering a therapeutically effective amount of the agent to an individual in need thereof. The invention provides for the use of the agent in the manufacture of a medicament for preventing or treating a disease in a particular individual.

[0119] The formulation of the drug depends on the nature of the drug. The drug is provided in the form of a pharmaceutical composition containing the drug and a pharma- ceutically acceptable carrier or diluent. Suitable carriers and diluents include isotonic saline, such as phosphate buffered saline. Typical oral dosage forms include tablets, capsules, liquid solutions, and liquid suspensions. The drug can be formulated for parenteral, intravenous, intramuscular, subcutaneous, transdermal, or oral administration.

[0120] The dosage of the drug can be determined according to various parameters, particularly according to the substance used; the age, weight, and disease of the individual to be treated; the route of administration; and the required regimen. A physician can determine the route of administration and dosage required for any particular drug. However, a suitable dosage may be, for example, 0.1 to 100 mg / kg of body weight, such as 1 to 40 mg / kg of body weight, taken one to three times a day.

[0121] The therapeutic agent may be any such agent disclosed herein or may target any "target" disclosed herein, including any protein or gene disclosed in any table herein, including Tables 5 or 6. It is understood that any agents disclosed in combination should also be considered to be disclosed for administration individually.

[0122] Prostate Cancer Treatment Prostate cancer treatment is recommended depending on the stage of disease progression. Radiation therapy, hormone therapy, and chemotherapy are the three options often used in prostate cancer treatment. A single treatment or a combination of treatments may be used.

[0123] chemotherapy Chemotherapy is often used to treat prostate cancer that has spread to other organs in the body (metastatic prostate cancer). Chemotherapy destroys cancer cells by interfering with the process of their growth. Chemotherapy does not cure prostate cancer, but it does keep it under control and relieve symptoms, therefore lessening its impact on daily life.

[0124] Radiation therapy This treatment can be used to cure localized and locally advanced prostate cancer. Radiation therapy can also be used to slow the progression and relieve symptoms of metastatic prostate cancer. Patients may receive hormone therapy before undergoing chemotherapy to increase the chances of successful treatment. Hormonal therapy may also be recommended after radiation therapy to reduce the chances of recurrence.

[0125] Hormone therapy Hormone therapy is often used in combination with radiation therapy. Hormone therapy should not usually be used alone to treat localized prostate cancer in men who are fit and willing to undergo surgery or radiation therapy. Hormone therapy can be used to slow the progression of advanced prostate cancer and relieve symptoms. Hormones control the growth of cells in the prostate. In particular, prostate cancer requires the hormone testosterone to grow. The goal of hormone therapy is to block the effects of testosterone by either stopping the production of testosterone or stopping the patient's body from using testosterone.

[0126] Other therapies that may be used to treat prostate cancer ·radical prostatectomy ·High-intensity focused ultrasound therapy ·Cryotherapy ·Brachytherapy ·Watchful waiting Transurethral resection of the prostate Treatment of advanced prostate cancer ·steroid

[0127] DLBCL therapy Four types of therapy may be used to treat DLBCL: - Chemotherapy - Radiation therapy - Monoclonal antibody therapy - Steroid therapy

[0128] Any of the above treatments may also be used to treat lymphoma.

[0129] Forms of matter referred to herein Any of the substances, such as nucleic acids or therapeutic agents, referred to herein may be in purified or isolated form. They may be in a form different from that found in nature, e.g., they may be present in combination with other substances that do not occur in nature. Nucleic acids (including portions of sequences defined herein) may have sequences different from those found in nature, e.g., have at least one, two, three, four or more nucleotide changes in the sequence as described in the section on homology. Nucleic acids may have heterologous sequences at the 5' or 3' ends. Nucleic acids may be chemically different from those found in nature, e.g., they may be modified in some way, but are still preferably capable of Watson-Crick base pairing. Where appropriate, nucleic acids are provided in double-stranded or single-stranded form. The invention provides all of the specific nucleic acid sequences referred to herein in single-stranded or double-stranded form, and thus includes the complementary strand to any sequence disclosed.

[0130] The present invention provides a kit for carrying out any process of the present invention, including the detection of chromosomal interactions related to prognosis.Such kits can include specific binding agents that can detect relevant chromosomal interactions, such as agents that can detect the ligated nucleic acid produced by the process of the present invention.Preferred agents present in the kit include probes that can hybridize to ligated nucleic acid or primer pairs that can amplify ligated nucleic acid in PCR reaction, for example as described herein.

[0131] The present invention provides a device capable of detecting relevant chromosomal interactions. The device preferably comprises any specific binding agent, probe or primer pair capable of detecting chromosomal interactions, such as any such agent, probe or primer pair described herein.

[0132] Detection Method In one embodiment, quantitative detection of ligated sequences associated with chromosomal interactions is carried out using a probe that is detectable upon activation during a PCR reaction, the ligated sequences comprising sequences from two chromosomal regions that come together in an epigenetic chromosomal interaction, the method comprising contacting the ligated sequences with a probe during a PCR reaction and detecting the degree of activation of the probe, the probe binding to the ligated site. The method typically uses a dual-labeled fluorescent hydrolysis probe to detect specific interactions in a MIQE-compliant manner.

[0133] Probe is generally labeled with a detectable label in inactive and active states, so that it is only detected when activated. The degree of activation is related to the degree of template (ligated product) present in PCR reaction. Detection can be carried out during all or part of PCR, for example, during at least 50% or 80% of the cycles of PCR.

[0134] The probe can include a fluorophore covalently attached to one end of the oligonucleotide and a quencher attached to the other end of the nucleotide, so that the fluorescence of the fluorophore is quenched by the quencher. In one embodiment, the fluorophore is attached to the 5' end of the oligonucleotide and the quencher is covalently attached to the 3' end of the oligonucleotide. Fluorophores that can be used in the method of the present invention include FAM, TET, JOE, Yakima Yellow, HEX, Cyanine 3, ATTO 550, TAMRA, ROX, Texas Red, Cyanine 3.5, LC 610, LC 640, ATTO 647N, Cyanine 5, Cyanine 5.5, and ATTO 680. Quenchers that can be used with the appropriate fluorophore include TAM, BHQ1, DAB, Eclipse, BHQ2, and BBQ650, and optionally the fluorophore is selected from HEX, Texas Red, and FAM. Preferred combinations of fluorophores and quenchers include FAM and BHQ1, and Texas Red and BHQ2.

[0135] Using probes in qPCR assays The hydrolysis probes of the present invention are typically optimized for temperature gradient with concentration-matched negative controls. Preferably, a single-step PCR reaction is optimized. More preferably, a standard curve is calculated. The advantage of using a specific probe that binds across the junction of the ligated sequences is that specificity for the ligated sequences can be achieved without using a nested PCR approach. The methods described herein allow accurate and precise quantification of low copy number targets. Prior to optimization of the temperature gradient, the ligated sequences of the target can be purified, e.g., gel purified. The ligated sequences of the target can be sequenced. Preferably, the PCR reaction is performed using about 10 ng, or 5-15 ng, or 10-20 ng, or 10-50 ng, or 10-200 ng of template DNA. The forward and reverse primers are designed such that one primer binds to one of the sequences of the chromosomal regions represented by the ligated DNA sequence and the other primer binds to the other chromosomal region represented by the ligated DNA sequence, e.g., by being complementary to the sequence.

[0136] Selection of ligated DNA targets The invention involves selecting primers and probes for use in the PCR methods defined herein, which includes selecting primers based on their ability to bind and amplify ligated sequences, and selecting the probe sequence based on the properties of the target sequence to which it binds, in particular the curvature of the target sequence.

[0137] Probes are typically designed / selected to bind to ligated sequences that are juxtaposed restriction fragments spanning restriction sites. In one embodiment of the invention, the predicted curvature of the potentially ligated sequences associated with a particular chromosomal interaction is calculated, for example, using certain algorithms referenced herein. Curvature can be expressed as degrees per helical turn, such as, for example, 10.5° per helical turn. Ligated sequences are selected for targeting if the ligated sequence has a curvature propensity peak score of at least 5° per helical turn, typically at least 10°, 15° or 20° per helical turn, for example, between 5° and 20° per helical turn. Preferably, the curvature propensity peak score per helical turn is calculated for at least 20, 50, 100, 200 or 400 bases, such as between 20 and 400 bases upstream and / or downstream of the ligated site. Thus, in one embodiment, the target sequence in the ligated product has any of these levels of curvature. Target sequences can also be selected based on the lowest thermodynamic structural free energy.

[0138] Specific Aspects In one embodiment, only intrachromosomal interactions are classified / detected and not extrachromosomal interactions (between different chromosomes).

[0139] In certain embodiments, a particular chromosomal interaction, e.g., any specific interaction referred to herein (e.g., as defined by any probe or primer pair referred to herein), is not classified. In some embodiments, a chromosomal interaction is not classified in any of the genes referred to herein.

[0140] The data provided herein show that the markers are "disseminating" and capable of distinguishing between cases and non-cases of the relevant disease state. Thus, when implementing the present invention, the skilled artisan will be able to determine which subgroup an individual is in by detecting interactions. In one embodiment, a threshold of detection of at least 70% of the tested markers in a form associated with the relevant disease state (either by absence or presence) can be used to determine whether an individual is in the relevant subgroup.

[0141] Screening Methods The present invention provides a method for determining which chromosomal interactions are associated with a chromosomal state corresponding to a prognostic subgroup of a population, the method comprising contacting a first set of nucleic acids from subgroups with different chromosomal states with a second set of index nucleic acids and allowing the complementary sequences to hybridize, the nucleic acids in the first and second sets of nucleic acids representing a ligated product comprising sequences from both chromosomal regions that are brought together in the chromosomal interaction, and the pattern of hybridization between the first set of nucleic acids and the second set of nucleic acids makes it possible to determine which chromosomal interactions are specific for a prognostic subgroup. The subgroups can be any of the specific subgroups defined herein, for example with reference to a particular disease or treatment.

[0142] Public bulletin The contents of all publications mentioned in this specification are incorporated herein by reference and may be used to further define features relevant to the present invention.

[0143] Specific Aspects The EpiSwitch™ platform technology detects epigenetic regulatory signatures of regulatory changes between normal and abnormal states at gene loci. The EpiSwitch™ platform identifies and monitors fundamental epigenetic levels of gene regulation associated with control of the higher order structure of human chromosomes, also known as chromosome conformation signatures. Chromosomal signatures are distinct key steps in the cascade of gene deregulation. These are higher order biomarkers with a unique set of advantages over biomarker platforms that utilize later epigenetic and gene expression biomarkers, such as DNA methylation and RNA profiling.

[0144] EpiSwitch™ Array Assays The custom EpiSwitch™ array screening platform has four densities of 15K, 45K, 100K, and 250K unique chromosome conformations, with each chimeric fragment repeated four times on the array, bringing the effective densities to 60K, 180K, 400K, and 1 million, respectively.

[0145] Custom-designed EpiSwitch™ arrays The 15K EpiSwitch™ arrays can screen the entire genome, including approximately 300 loci interrogated with the EpiSwitch™ Biomarker discovery technology. EpiSwitch™ arrays are built on the Agilent SurePrint G3 Custom CGH microarray platform, which offers four densities, 60K, 180K, 400K, and 1 million probes. The density per array is reduced to 15K, 45K, 100K, and 250K, as each EpiSwitch™ probe is displayed as a quadruplicate, thus allowing for statistical assessment of reproducibility. The average number of potential EpiSwitch™ markers interrogated per locus is 50, and the number of loci that can be investigated is 300, 900, 2000, and 5000.

[0146] EpiSwitch™ Custom Array Pipeline The EpiSwitch™ array is a dual-color system with one set of samples labeled in Cy5 and another set of samples to be compared / analyzed (controls) labeled in Cy3 after generation of the EpiSwitch™ library. The arrays are scanned using an Agilent SureScan Scanner and the resulting features are extracted using Agilent Feature Extraction software. The data is then processed using the EpiSwitch™ array processing script in R. The arrays are processed using the standard dual-color package in Bioconductor in R: Limma*. Array normalization is performed using the normalizedWithinArrays function in Limma* against the on-chip Agilent positive control and the EpiSwitch™ positive control. To analyze using Limma*, the data is filtered based on Agilent Flag calls, Agilent control probes are removed, and technical replicate probes are averaged. Probes are modeled based on the differences between the two scenarios being compared and corrected using the false discovery rate. Probes with coefficient of variation (CV) <=30% passing p-value of <=-1.1 or =>1.1 and p<=0.1 FDR are used for further screening. Further multiple factor analysis is performed using FactorMineR package in R to reduce the probe set.

[0147] *Note: LIMMA is a linear model and empirical Bayesian process for assessing differential expression in microarray experiments. Limma is an R package for analyzing gene expression data resulting from microarrays or RNA-Seq.

[0148] A pool of probes is initially selected based on the adjusted p-value, FC and CV < 30% (arbitrary cutoff point) parameters for final picking. Further analysis and the final list are made based only on the first two parameters (adjusted p-value; FC).

[0149] Statistics Pipeline The EpiSwitch™ screening arrays are processed using the EpiSwitch™ analysis package in R to select high value EpiSwitch™ markers for conversion to the EpiSwitch™ PCR platform.

[0150] Process 1 Probes are selected based on their corrected p-value (false discovery rate, FDR), which is the product of a modified linear regression model. Probes with p-value <=0.1 are selected and then further reduced by epigenetic ratio (ER), probes ER must be <=-1.1 or =>1.1 to be selected for further analysis. The final filter is the coefficient of variation (CV), probes must be <=0.3.

[0151] Process 2 The top 40 markers from the statistical list are selected based on their ER for selection as markers for PCR conversion. The top 20 markers with the highest negative ER load and the top 20 markers with the highest positive ER load form the list.

[0152] Process 3 The resulting markers, statistically significant probes from step 1, form the basis of an enrichment analysis using hypergeometric enrichment (HE), which allows marker reduction from the significant probe list and, together with the markers from step 2, forms a list of probes that are transferred to the EpiSwitch™ PCR platform.

[0153] Statistical probes are processed by HE to determine which genetic locations have statistically significant probe enrichment, indicating which genetic locations are hubs of epigenetic differences.

[0154] The most significant enrichment loci based on corrected p-values ​​are selected for creation of a probe list. Genetic positions with p-values ​​less than 0.3 or 0.2 are selected. Statistical probes that map to these genetic positions using markers from step 2 form high-value markers for EpiSwitch™ PCR conversion.

[0155] Array Design and Processing Array Design 1. Loci are processed using SII software (currently v3.2) as follows: a. Extract genomic sequence at these specific loci (50 kb upstream and 20 kb downstream gene sequence). b. Define the probability that sequences within this region are involved in CC. c. A specific RE is used to cut the sequence. d. Determine which restriction fragments tend to interact in particular orientations. e. Rank the likelihood that different CCs interact together. 2. Determine the array size and therefore the number of available probe positions (x). 3. Elicit x / 4 interactions. 4. For each interaction, define 30bp of sequence from part 1 to the restriction site and 30bp to the restriction site from part 2. Make sure those regions are not repeats, if they are, exclude them and write down the next interaction on the list. Join both 30bp to define the probe. 5. Create a list of x / 4 probes plus defined control probes, replicated four times, to be created on the array. 6. Upload the list of probes to the Agilent Sure design for custom CGH arrays website. 7. Design Agilent custom CGH arrays using the probe sets.

[0156] Array Processing 1. Process samples using EpiSwitch™ Standard Operating Procedures (SOPs) for template creation. 2. Clean up with ethanol precipitation by the array processing lab. 3. Process samples according to Agilent SureTag Complete DNA Labeling Kit - Agilent Oligonucleotide Array-based CGH for Genomic DNA Analysis Enzymatic labelling for Blood, Cells or Tissues. 4. Scan using an Agilent C scanner with Agilent feature extraction software.

[0157] EpiSwitch™ biomarker signatures demonstrate high robustness, sensitivity, and specificity in stratifying complex disease phenotypes. The technology takes advantage of the latest breakthrough in the science of epigenetics, the monitoring and evaluation of chromosome conformation signatures as a highly informative class of epigenetic biomarkers. Current research methodologies deployed in academic settings require 3-7 days of biochemical processing of cellular material to detect CCS. These procedures are limited in sensitivity and reproducibility and, furthermore, do not have the benefit of the targeted insights offered by the EpiSwitch™ analytical package at the design stage.

[0158] In silico marker identification for EpiSwitch™ arrays Genome-wide CCS sites are assessed by EpiSwitch™ arrays directly on clinical samples from the study cohort for identification of all relevant stratification lead biomarkers. The EpiSwitch™ array platform is used for marker identification due to its high throughput capacity and its ability to rapidly screen a large number of loci. The arrays used were Agilent custom CGH arrays, which allow for cross-checking of markers identified by the in silico software.

[0159] EpiSwitch(TM) PCR Potential markers identified by the EpiSwitch™ Array are then validated by either EpiSwitch™ PCR or a DNA sequencer (i.e., Roche 454, Nanopore MinION, etc.). Top PCR markers that are statistically significant and show the best reproducibility are selected for further reduction into a final EpiSwitch™ signature set and validated in an independent cohort of samples. EpiSwitch™ PCR can be performed by trained technicians following established standardized operating procedure protocols. Manufacturing of all protocols and reagents is performed under ISO 13485 and 9001 certification to ensure the quality of the research and the ability to transfer the protocols. The EpiSwitch™ PCR and EpiSwitch™ Array biomarker platforms are compatible with both whole blood and cell line analysis. The tests are sensitive enough to detect abnormalities at very low copy numbers using small amounts of blood.

[0160] Paragraph showing an embodiment of the present invention 1. A process for detecting a chromosomal state representing a subgroup within a population, the process comprising determining whether a chromosomal interaction associated with the chromosomal state is present within a defined region of the genome; and - the chromosomal interactions are optionally identified by a method of determining which chromosomal interactions are associated with chromosomal states corresponding to subgroups of a population, the method comprising contacting a first set of nucleic acids from subgroups with different chromosomal states with a second set of index nucleic acids and allowing hybridization of complementary sequences, the nucleic acids in the first and second sets of nucleic acids representing ligated products comprising sequences from both chromosomal regions that have come together in the chromosomal interaction, and a pattern of hybridization between the first set of nucleic acids and the second set of nucleic acids allows the determination of which chromosomal interactions are specific to the subgroup, and - Subgroups are associated with prognosis in prostate cancer, and chromosomal interactions are (i) in any of the regions or genes listed in Table 6; and / or (ii) corresponds to any of the chromosomal interactions represented by any of the probes shown in Table 6; and / or (iii) is present in a 4,000-base region including or adjacent to (i) or (ii); or - Subgroups are associated with DLBCL prognosis and chromosomal interactions a) in any of the regions or genes listed in Table 5; and / or b) corresponds to any of the chromosomal interactions represented by any of the probes shown in Table 5; and / or c) a process occurring in a 4,000 base region that includes or is adjacent to (a) or (b).

[0161] 2.- The prognosis of the above-mentioned prostate cancer is related to whether the cancer is aggressive or indolent; and / or - The above DLBCL prognosis is related to survival 2. The process described in paragraph 1.

[0162] 3. The specific combination of chromosomal interactions that are associated with prostate cancer in subgroups (i) containing all of the chromosomal interactions represented by the probes in Table 6, and / or (ii) comprising at least one, two, three, or four of the chromosomal interactions represented by the probes in Table 6; and / or (iii) are present together in at least one, two, three, or four of the regions or genes listed in Table 6; and / or (iv) at least one, two, three, or four of the classified chromosomal interactions are present in a 4,000 base region that includes or is adjacent to a chromosomal interaction represented by a probe in Table 6. 3. The process according to paragraph 1 or 2,

[0163] 4. Subgroups associated with DLBCL and specific combinations of chromosomal interactions (i) containing all of the chromosomal interactions represented by the probes in Table 5; and / or (ii) comprising at least 10, 20, 30, 50, or 80 of the chromosomal interactions represented by the probes in Table 5; and / or (iii) are present together in at least 10, 20, 30, or 50 of the regions or genes listed in Table 5; and / or (iv) at least 10, 20, 30, 50, or 80 of the classified chromosomal interactions are present in a 4,000 base region that includes or is adjacent to a chromosomal interaction represented by a probe in Table 5. 3. The process according to paragraph 1 or 2,

[0164] 5. Subgroups associated with DLBCL and specific combinations of chromosomal interactions (i) contains all of the chromosomal interactions shown in Table 7, and / or (ii) includes at least one, two, five, or eight of the chromosomal interactions set forth in Table 7. 3. The process according to paragraph 1 or 2,

[0165] 6. At least 10, 20, 30, 40 or 50 chromosomal interactions are classified, preferably at least 10 chromosomal interactions are classified 6. The process of any one of paragraphs 1 to 5.

[0166] 7. Chromosomal interactions - in a sample from an individual, and / or - by detecting the presence or absence of DNA loops at the sites of chromosomal interactions, and / or - by detecting the presence or absence of distal regions of chromosomes that are in conjunction with chromosomal conformation; and / or - sorting by detecting the presence of ligated nucleic acids, the ligated nucleic acids being generated during said sorting and comprising two regions each corresponding to regions of the chromosomes whose sequences come together in a chromosomal interaction, the detection of the ligated nucleic acids preferably comprising: (i) with a probe having at least 70% identity to any of the specific probe sequences referred to in Table 6, and / or (ii) with a primer pair having at least 70% identity to any of the primer pairs in Table 6, in the case of prostate cancer prognosis; or (ii) with a probe having at least 70% identity with any of the specific probe sequences mentioned in Table 5, and / or (b) with a primer pair having at least 70% identity with any of the primer pairs in Table 5, in the case of prognosis of DLBCL. 7. The process of any one of paragraphs 1 to 6.

[0167] 8.- The second set of nucleic acids is from a larger group of individuals than the first set of nucleic acids, and / or - the first set of nucleic acids is from at least eight individuals, and / or - the first set of nucleic acids is from at least four individuals from a first subgroup and at least four individuals from a second subgroup that preferably does not overlap with the first subgroup, and / or - a process is performed to select individuals for medical treatment 8. The process of any one of paragraphs 1 to 7.

[0168] 9.- The second set of nucleic acids represents an unselected group, and / or - a second set of nucleic acids is attached to the array at defined locations; and / or - the second set of nucleic acids represents chromosomal interactions in at least 100 different genes; and / or - the second set of nucleic acids comprises at least 1,000 different nucleic acids representing at least 1,000 different chromosomal interactions; and / or - the first set of nucleic acids and the second set of nucleic acids comprise at least 100 nucleic acids having a length between 10 and 100 nucleotide bases. 9. The process of any one of paragraphs 1 to 8.

[0169] 10. The first set of nucleic acids comprises: (i) bridging chromosomal regions that are brought together in chromosomal interactions; (ii) exposing said cross-linked region to cleavage, optionally by enzymatic restriction digestion; (iii) ligating the cross-linked cut DNA ends to form a first set of nucleic acids (particularly comprising ligated DNA). 10. The process of any one of paragraphs 1 to 9.

[0170] 11. The above-defined region of the genome (i) contains a single nucleotide polymorphism (SNP); and / or (ii) expresses microRNA (miRNA); and / or (iii) expresses a non-coding RNA (ncRNA); and / or (iv) expressing a nucleic acid sequence encoding at least 10 contiguous amino acid residues; and / or (v) represents a control element; and / or (vii) contains a CTCF binding site 11. The process of any one of paragraphs 1 to 10.

[0171] 12. The process of any one of paragraphs 1 to 11, carried out to determine whether prostate cancer is aggressive or indolent, comprising classifying at least five chromosomal interactions as defined in Table 6.

[0172] 13. A process according to any one of paragraphs 1 to 12, carried out for determining the prognosis of DLBLC, comprising classifying at least five chromosomal interactions as defined in Table 5.

[0173] 14.- Preferably, the above process is used to detect whether a candidate agent is capable of causing a change to a chromosomal state that is associated with different levels of prognosis, - the chromosomal interaction is represented by any of the probes of Table 6, and / or - the chromosomal interaction is present in any region or gene listed in Table 6, and optionally, - the chromosomal interaction is identified by a method for determining that the chromosomal interaction is associated with a chromosomal condition as defined in paragraph 1, and / or - carried out to identify or design therapeutic agents against prostate cancer, wherein changes in chromosomal interactions are monitored by (i) a probe having at least 70% identity with any of the probe sequences mentioned in Table 6, and / or (ii) a primer pair having at least 70% identity with any of the primer pairs in Table 6. 14. The process of any one of paragraphs 1 to 13.

[0174] 15.- Preferably, the above process is used to detect whether a candidate agent is capable of causing a change to a chromosomal state that is associated with different levels of prognosis, - the chromosomal interaction is represented by any of the probes of Table 5, and / or - the chromosomal interaction is present in any region or gene listed in Table 5, and optionally, - the chromosomal interaction is identified by a method for determining that the chromosomal interaction is associated with a chromosomal condition as defined in paragraph 1, and / or - a method for identifying or designing therapeutic agents against DLBCL, in which changes in chromosomal interactions are monitored by (i) a probe having at least 70% identity with any of the probe sequences mentioned in Table 5, and / or (ii) a primer pair having at least 70% identity with any of the primer pairs in Table 5; 14. The process of any one of paragraphs 1 to 13.

[0175] 16. A process according to paragraph 14 or 15, comprising a step of selecting a target based on detection of a chromosomal interaction, and preferably a step of screening modulators of the target to identify therapeutic agents for immunotherapy, said target optionally being a protein.

[0176] 17. The classification or detection includes specific detection of the ligated product by quantitative PCR (qPCR) using primers capable of amplifying the ligated product and a probe that binds to the ligated site during a PCR reaction, the probe comprising a sequence complementary to a sequence from each of the chromosomal regions that are brought together in the chromosomal interaction, preferably the probe comprises: an oligonucleotide that specifically binds to the ligated product; and / or a fluorophore covalently attached to the 5' end of the oligonucleotide, and / or comprising a quencher covalently attached to the 3' end of the oligonucleotide; and Optionally, the fluorophore is selected from HEX, Texas Red and FAM; and / or 17. The process according to any one of paragraphs 1 to 16, wherein the probe comprises a nucleic acid sequence of 10 to 40 nucleotide bases in length, preferably 20 to 30 nucleotide bases in length.

[0177] 18.- The results of the process are provided in a report, and / or - the results of the process are used to select a treatment schedule for the patient, and preferably to select a specific treatment for the individual; 18. The process of any one of paragraphs 1 to 17.

[0178] 19. A therapeutic agent for use in a method of treating prostate cancer or DLBCL in an individual by the process of any one of paragraphs 1-13 and 17, wherein the individual has been identified as in need of the therapeutic agent.

[0179] 2-1. A process for detecting a chromosomal state representing a subgroup within a population, comprising determining whether a chromosomal interaction associated with the chromosomal state is present within a defined region of the genome; - the chromosomal interactions are optionally identified by a method of determining which chromosomal interactions are associated with chromosomal states corresponding to subgroups of a population, the method comprising contacting a first set of nucleic acids from subgroups with different chromosomal states with a second set of index nucleic acids and allowing hybridization of complementary sequences, the nucleic acids in the first and second sets of nucleic acids representing ligated products comprising sequences from both chromosomal regions that have come together in the chromosomal interaction, and a pattern of hybridization between the first set of nucleic acids and the second set of nucleic acids allows determining which chromosomal interactions are specific to the subgroup, and - the subgroup is associated with prognosis of prostate cancer, and the chromosomal interaction is (i) in any of the regions or genes listed in Table 6; and / or (ii) corresponds to any of the chromosomal interactions represented by any of the probes shown in Table 6; and / or (iii) is present in a 4,000-base region that includes or is adjacent to (i) or (ii); or - the subgroup is associated with prognosis of DLBCL, and the chromosomal interaction is a) in any of the regions or genes listed in Table 5; and / or b) corresponds to any of the chromosomal interactions represented by any of the probes shown in Table 5; and / or c) is present in a 4,000 base region that includes or is adjacent to (a) or (b); or - the subgroup is associated with lymphoma prognosis, and the chromosomal interaction is (iv) is in any of the regions or genes listed in Table 8; and / or (v) corresponds to any of the chromosomal interactions shown in Table 8; and / or (vi) A process occurring in a 4,000 base region that includes (iv) or (v) or is adjacent to (iv) or (v). 2-2. The prognosis of the prostate cancer is related to whether the cancer is aggressive or indolent, and / or The process of paragraph 2-1, wherein the prognosis of DLBCL is related to survival. 2-3. The subgroup is associated with prostate cancer, and a specific combination of chromosomal interactions is (i) containing all of the chromosomal interactions represented by the probes in Table 6, and / or (ii) comprising at least one, two, three, or four of the chromosomal interactions represented by the probes in Table 6; and / or (iii) are present together in at least one, two, three, or four of the regions or genes listed in Table 6; and / or (iv) at least one, two, three, or four of the classified chromosomal interactions are present in a 4,000 base region that includes or is adjacent to a chromosomal interaction represented by a probe in Table 6. 2. The process according to paragraph 2-1 or 2-2, 2-4. The subgroup is associated with DLBCL, and a specific combination of chromosomal interactions is (i) containing all of the chromosomal interactions represented by the probes in Table 5; and / or (ii) comprising at least 10, 20, 30, 50, or 80 of the chromosomal interactions represented by the probes in Table 5; and / or (iii) are present together in at least 10, 20, 30, or 50 of the regions or genes listed in Table 5; and / or (iv) at least 10, 20, 30, 50, or 80 of the classified chromosomal interactions are present in a 4,000 base region that includes or is adjacent to a chromosomal interaction represented by a probe in Table 5. 2. The process according to paragraph 2-1 or 2-2, 2-5. The subgroup is associated with DLBCL, and a specific combination of chromosomal interactions is (i) contains all of the chromosomal interactions shown in Table 7, and / or (ii) comprising at least one, two, five or eight of the chromosomal interactions set forth in Table 7; 2. The process according to paragraph 2-1 or 2-2, 2-6. The subgroup is associated with lymphoma, and a specific combination of chromosomal interactions is (i) contains all of the chromosomal interactions shown in Table 8, and / or (ii) comprising at least 10, 20, 30 or 50 of the chromosomal interactions set forth in Table 8; and / or (iii) are present together in at least 10, 20, or 30 of the regions or genes listed in Table 8; and / or (iv) at least 10, 20, 30, or 50 of the chromosomal interactions typed are present in a 4,000 base region that includes or is adjacent to a chromosomal interaction set forth in Table 8; or preferably, the specific combination of chromosomal interactions is (a) contains all of the chromosomal interactions shown in Table 9, and / or (b) containing at least 5, 10 or 15 of the chromosomal interactions set forth in Table 9; and / or (c) are present together in at least 5, 10, or 15 of the regions or genes listed in Table 9; and / or (d) at least 5, 10, or 15 of the classified chromosomal interactions are present in a 4,000 base region that includes or is adjacent to a chromosomal interaction shown in Table 9. The process according to any one of paragraphs 2-1 to 2-5, 2-7. The process according to any of paragraphs 2-1 to 2-6, wherein at least 10, 20, 30, 40 or 50 chromosomal interactions are classified, preferably at least 10 chromosomal interactions are classified. 2-8.Chromosomal interactions - in a sample from an individual, and / or - by detecting the presence or absence of DNA loops at the sites of chromosomal interactions, and / or - by detecting the presence or absence of distal regions of chromosomes that are in conjunction with chromosomal conformation; and / or - by detecting the presence of ligated nucleic acids, said ligated nucleic acids being generated during said sorting and comprising two regions each corresponding to regions of chromosomes whose sequences come together in a chromosomal interaction, and said detection of the ligated nucleic acids preferably comprises: (i) with a probe having at least 70% identity to any of the specific probe sequences referred to in Table 6, and / or (ii) with a primer pair having at least 70% identity to any of the primer pairs in Table 6, in the case of prostate cancer prognosis; or (ii) with a probe having at least 70% identity to any of the specific probe sequences referred to in Table 5, and / or (ii) with a primer pair having at least 70% identity to any of the primer pairs in Table 5. 2-9. A chromosomal interaction is classified by detecting the presence of a ligated nucleic acid, the ligated nucleic acid being generated during the classification and comprising two regions whose sequences each correspond to a region of the chromosome that comes together in the chromosomal interaction, and the detection of the ligated nucleic acid in the case of lymphoma prognosis is - by a probe having at least 70% identity with any of the specific probe sequences mentioned in Table 5, and / or - by a primer pair having at least 70% identity with any primer pair mentioned in Table 5, and / or - the process of any of paragraphs 2-1 to 2-8, carried out with a primer pair having at least 70% identity to any primer pair of Table 8. 2-10.- The second set of nucleic acids is from a larger group of individuals than the first set of nucleic acids, and / or - said first set of nucleic acids is from at least eight individuals, and / or - the first set of nucleic acids is from at least four individuals from a first subgroup and at least four individuals from a second subgroup that preferably does not overlap with the first subgroup, and / or - The process according to any of paragraphs 2-1 to 2-9, wherein the process is carried out to select an individual for a medical procedure. 2-11.- The second set of nucleic acids represents an unselected group, and / or - the second set of nucleic acids is bound to the array at defined locations; and / or - the second set of nucleic acids represents chromosomal interactions in at least 100 different genes; and / or - the second set of nucleic acids comprises at least 1,000 different nucleic acids representing at least 1,000 different chromosomal interactions; and / or - the process according to any of paragraphs 2-1 to 2-10, wherein the first set of nucleic acids and the second set of nucleic acids comprise at least 100 nucleic acids having a length of between 10 and 100 nucleotide bases. 2-12. The first set of nucleic acids (i) bridging chromosomal regions that are brought together in chromosomal interactions; (ii) exposing the cross-linked region to cleavage, optionally by enzymatic restriction digestion; (iii) ligating the cross-linked cleaved DNA ends to form a first set of nucleic acids (particularly comprising ligated DNA). 2-13. The above-defined region of the genome (i) contains a single nucleotide polymorphism (SNP); and / or (ii) expresses microRNA (miRNA); and / or (iii) expresses a non-coding RNA (ncRNA); and / or (iv) expressing a nucleic acid sequence encoding at least 10 contiguous amino acid residues; and / or (v) represent a coordinating factor; and / or (vii) The process of any of paragraphs 2-1 to 2-12, comprising a CTCF binding site. 2-14. A process according to any one of paragraphs 2-1 to 2-13, carried out to determine whether a prostate cancer is aggressive or indolent, comprising classifying at least five chromosomal interactions as defined in Table 6. 2-15. A process according to any of paragraphs 2-1 to 2-14, carried out for determining the prognosis of DLBLC, comprising classifying at least five chromosomal interactions as defined in Table 5. 2-16.- Preferably, the process is used to detect whether a candidate agent is capable of causing a change to a chromosomal state that is associated with different levels of prognosis, - the chromosomal interaction is represented by any of the probes of Table 6, and / or - the chromosomal interaction is present in any region or gene listed in Table 6, and optionally, - the chromosomal interaction is identified by a method for determining that the chromosomal interaction is associated with a chromosomal condition as defined in claim 1, and / or - The process according to any of paragraphs 2-1 to 2-15, carried out for identifying or designing therapeutic agents for prostate cancer, wherein changes in chromosomal interactions are monitored by (i) a probe having at least 70% identity to any of the probe sequences mentioned in Table 6, and / or (ii) a primer pair having at least 70% identity to any of the primer pairs in Table 6. 2-17.- Preferably, the process is used to detect whether a candidate agent is capable of causing a change to a chromosomal state that is associated with different levels of prognosis, - the chromosomal interaction is represented by any of the probes of Table 5, and / or - the chromosomal interaction is present in any region or gene listed in Table 5, and optionally, - the chromosomal interaction is identified by a method for determining that the chromosomal interaction is associated with a chromosomal condition as defined in claim 1, and / or - the process according to any of paragraphs 2-1 to 2-15, carried out for identifying or designing therapeutic agents against DLBCL, wherein changes in chromosomal interactions are monitored by (i) a probe having at least 70% identity with any of the probe sequences mentioned in Table 5, and / or (ii) a primer pair having at least 70% identity with any of the primer pairs in Table 5. 2-18.- Preferably, the process is used to detect whether a candidate agent is capable of causing a change to a chromosomal state that is associated with different levels of prognosis, - the chromosomal interaction is represented by any of the probes of Tables 8 or 9, and / or - the chromosomal interaction is present in any region or gene listed in Table 8 or 9, and optionally, - the chromosomal interaction is identified by a method for determining that the chromosomal interaction is associated with a chromosomal condition as defined in claim 1, and / or - the process according to any of paragraphs 2-1 to 2-15, carried out for identifying or designing therapeutic agents against lymphoma, wherein changes in chromosomal interactions are monitored by (i) a probe having at least 70% identity with any of the probe sequences mentioned in Table 5, and / or (ii) a primer pair having at least 70% identity with any of the primer pairs of Tables 5 or 8. 2-19. A process according to any of paragraphs 2-16 to 2-18, comprising a step of selecting a target based on detection of a chromosomal interaction, and preferably a step of screening modulators of said target to identify a therapeutic agent for immunotherapy, said target optionally being a protein. 2-20. The process of any of paragraphs 2-1 to 2-19, wherein the prognosis is in humans or dogs. 2-21. The classification or detection comprises specific detection of the ligated product by quantitative PCR (qPCR) using primers capable of amplifying the ligated product and a probe that binds to the ligated site during a PCR reaction, the probe comprising a sequence complementary to a sequence from each of the chromosomal regions that are brought together in a chromosomal interaction, preferably the probe comprises: an oligonucleotide that specifically binds to the ligated product; and / or a fluorophore covalently attached to the 5' end of said oligonucleotide, and / or a quencher covalently attached to the 3' end of the oligonucleotide; and optionally, the fluorophore is selected from HEX, Texas Red and FAM; and / or 21. The process according to any of paragraphs 2-1 to 2-20, wherein the probe comprises a nucleic acid sequence of 10 to 40 nucleotide bases in length, preferably 20 to 30 nucleotide bases in length. 2-22.- The results of said process are provided in a report; and / or - The process according to any of paragraphs 2-1 to 2-21, wherein the results of said process are used to select a treatment schedule for the patient, and preferably to select a specific treatment for the individual. 2-23. A therapeutic agent for use in a method of treating prostate cancer, DLBCL or lymphoma in an individual by a process described in any of paragraphs 2-1 to 2-15, 2-20 and 2-21, wherein the individual has been identified as being in need of the therapeutic agent. 2-24.- The subgroup is associated with prostate cancer and at least one chromosomal interaction from table 25 is classified; and / or - said subgroup is related to prostate cancer and at least one of the following interaction combinations from table 25 is classified: (i) ETS1, MAP3K14, SLC22A3, and CASP2, or (ii) BMP6, ERG, MSR1, MUC1, ACAT1, and DAPK1, or (iii) HSD3B2, VEGFC, APAF1, MUC1, ACAT1, and DAPK1, and / or - said subgroup is associated with DLBCL, preferably classifying at least one of the first 10 markers shown in table 5 corresponding to one or more of the following genes: STAT3, TNFRSF13B, ANXA11, MAP3K7, MEF2B and IFNAR1, and / or - the subgroup is associated with lymphoma, preferably classifying at least one of the first eleven markers shown in Table 6 corresponding to one or more of the following genes: STAT3, TNFRSF13B, ANXA11, MAP3K7, MEF2B and IFNAR1, the process or therapeutic agent according to any of paragraphs 2-1 to 2-23. EXAMPLES

[0180] The invention is illustrated by the following examples. Example 1 Use of EpiSwitch™ (chromosome conformation signature) markers We have consistently observed highly popular EpiSwitch™ markers with high concordance and strong validation with primary and secondary diseased tissues. EpiSwitch™ biomarker signatures have demonstrated high robustness, sensitivity, and specificity in stratifying complex disease phenotypes.

[0181] EpiSwitch™ technology provides a highly effective screening tool, early detection, companion diagnostics, monitoring and prognostic analysis of major diseases associated with aberrant and responsive gene expression. The main advantages of the OBD approach are that it is non-invasive, rapid and relies on highly stable DNA-based targets as part of chromosomal signatures rather than unstable protein / RNA molecules.

[0182] CCS form a stable regulatory framework of epigenetic regulation and access to genetic information across the entire genome of a cell. Changes in CCS reflect early changes in the mode of regulation and gene expression long before their consequences spontaneously manifest as obvious abnormalities. A simple way of thinking about CCS is that they are topological arrangements in which different, distant regulatory parts of DNA are in close proximity to each other and affect each other's functions. These connections are not made randomly, but are highly regulated and well recognized as high-level regulatory mechanisms with significant biomarker stratification power.

[0183] Prognostic stratification of prostate cancer Markers were developed based on retrospective annotation of class I (low risk, indolent), class II (intermediate), and class III (high risk of aggressive). The markers show robust classification of patients versus healthy controls and differentiate between classes. Samples were from the UK.

[0184] Identification of EpiSwitch™ biomarkers capable of distinguishing blood from prostate cancer patients and healthy controls Initially, a custom EpiSwitch™ microarray study was used to identify and screen ~15,000 potential CCS across 425 loci to differentiate between 8 prostate cancer (PCa) and 8 control individuals. The top statistically significant markers were transferred to a nested PCR assay and screened in a larger sample cohort of 24 PCa and 25 healthy control samples. A classifier was developed using the top 5 CCS transferred from the microarray that classified PCa and control samples with a sensitivity and specificity of 100% (95% CI -86.2% to 100%) and 100% (95% CI -86.7% to 100%), respectively.

[0185] Figure 1 shows the principal components analysis of the top five markers for the 49 samples from the development sample cohort.

[0186] The diagnostic classifier was used to classify an additional blinded independent cohort (n=29) of 24 PCa and 5 healthy control samples with 83% accuracy. Further development of the EpiSwitch™ Prostate Cancer Assay was performed with an additional sample cohort (n=192) of 95 PCa and 97 controls. This was validated in a blinded sample cohort of 20 samples (10 PCa, 10 controls). The results of this validation are shown in Table 1.

[0187] [Table 1]

[0188] The latest project in the PCA program developed an alternative PCR format for PCa diagnosis utilizing hydrolysis probe-based real-time quantitative PCR (qPCR). The performance of the 6-marker model is shown in Table 2.

[0189] [Table 2]

[0190] overview Three independent, blinded validations of the EpiSwitch™ PCa diagnostic signature developed during the PCa diagnostic program using US and UK samples from various stages of disease achieve >80% sensitivity and specificity for diagnosing prostate cancer. The gold standard clinical assay for detecting PCa, the prostate-specific antigen (PSA) blood test, is itself dependent on a variety of other variables and typically has a sensitivity and specificity range of 32-68%. Furthermore, a parallel research track has resulted in the development of EpiSwitch™ assays that assess prostate cancer prognosis to aid in clinical management and treatment selection for individual patients diagnosed with PCA.

[0191] An additional custom EpiSwitch microarray study was performed to identify and screen ~15,000 potential CCS across 426 loci to distinguish between 8 aggressive prostate cancer patients (Class 3) and 8 indolent PCa patients (Class 1). A description of the PCa classes can be found in the attachment. The top statistically significant markers were transferred into nested PCR assays and screened in a larger sample cohort of 42 Class 1, 25 Class 2 and 19 Class 3 PCa samples.

[0192] The top six statistically significant markers were used to develop a prognostic classifier for classifying class 1 (low risk) and class 3 (high risk) PCa. The performance of the classifier on an independent sample cohort (n=27) of 42 class 1 and 25 class 3 samples is shown in Table 3.

[0193] [Table 3]

[0194] In an alternative analysis, six additional markers were found that stratified between class 2 and class 3 PCa. Two classifiers shared two markers, and each classifier also had four unique markers.

[0195] FIG. 2 shows a VENN comparison of the two PCA prognostic classifiers.

[0196] The performance of the PCa classifiers for class 2 vs. class 3 is shown in Table 4 .

[0197] [Table 4]

[0198] conclusion Diagnostic and prognostic biomarker development was achieved in multiple clinical sample cohorts. All marker screening and selection performed was based on systemic, blood-based epigenetic alterations monitored via chromosome conformation signatures in patients with different stages of prostate cancer (stages 1-3) versus healthy controls (diagnostic use), as well as in aggressive, high-risk category 3 versus indolent, low-risk category 1 prostate cancer (prognostic use), or intermediate-risk category 2 patients.

[0199] Results of stratification development for PCa versus healthy controls showed sensitivity and specificity of up to >80% in test cohorts and in series of blinded validation: stratification of PCa into high-risk category 3 versus low-risk category 1 showed sensitivity of up to 80% and specificity of up to 92% in cohorts of up to 67 samples, while stratification of high-risk category 3 versus intermediate-risk category 2 showed sensitivity of up to 84% and specificity of up to 88% in cohorts of up to 44 samples.

[0200] attached Low Risk - Category 1 Localized prostate cancer is classified as low risk if: A PSA level less than 10 ng per milliliter, and Gleason score ≤ 6, and The T stage is between T1 and T2a.

[0201] Medium Risk - Category 2 Localized prostate cancer is classified as intermediate risk if at least one of the following is present: PSA level is between 10 and 20 ng / ml. Gleason score of 7 The T stage is T2b.

[0202] High Risk - Category 3 Localized prostate cancer is classified as high risk if at least one of the following is present: A PSA level of 20 ng / ml or higher Gleason score between 8 and 10 T stage is T2c, T3, or T4.

[0203] If the cancer is at stage T3 or T4, this means it has broken through the prostate's outer fibrous covering (capsule), so it is classified as locally advanced prostate cancer.

[0204] Example 2. Identification of markers for DLBCL overview This is relevant for the identification of major groups of patients with poor and good prognosis for the selection of subsequent treatment (i.e. R-CHOP). Biomarkers have been developed based on retrospective overall survival. Usually, patients are classified by biopsy-based gene expression criteria such as Nanostring or Fluidigm according to disease subtypes such as ABC (poor prognosis) or GCB (good prognosis). However, not all patients could be classified as ABC or GCB (so-called type III, or unclassified patients). We identified biomarkers that provide a classification for the prognosis of survival at pre-treatment baseline, regardless of the standard classification of ABC or GCB.

[0205] Marker identification DLBCL shows clear differences in patient survival (poor vs. favorable prognosis) and is characterized by numerous molecular readouts into subtypes. The various subtypes are also treated differently in current clinical practice. This includes, for example, the combination of rituximab and CHOP for chemotherapy. Different approaches exist.

[0206] Currently implemented molecular readouts are based on gene expression profiling by arrays, performed on biological material obtained by direct biopsy. They include Nanostring and Fluidigm array-based tests for the extreme types ABC and GCB. ANC subtypes are usually associated with poor prognosis. Not all patients could be classified as ABC or GCB, and a large number of patients remain unclassified (or type III) with regard to the established gene expression profile and any association with poor prognosis of survival. We constructed a systemic biomarker that directly classifies patients with poor vs. good prognosis, regardless of transcriptional gene expression profiling by other modalities.

[0207] Step 1: Epigenetic profiles of groups of cell lines representing poor and good prognosis of survival for DLBCL were compared using Episwitch screening arrays, allowing for the identification of array-based markers and the design of nested PCR primers to be used against the same targets in a PCR format.

[0208] Step 2: We used the top 10 nested PCR-based markers read in baseline blood samples from 57–58 unclassified DLBCL patients with known retrospective survival annotation. Table 6 provides details on the markers, final signatures, and performance by the classifier model.

[0209] The study shows how the baseline call for these patients with poor / good prognosis compares with the clinical survival data. This is a Cox estimate of the hazard ratio, i.e., the baseline classification into poor prognosis has a certain value >1, indicating a higher probability of being in the poor survival group rather than the good group by clinical post annotation. The latter is of particular value and interest to the clinical team in the trial design.

[0210] Detailed commentary Diffuse large B-cell lymphoma (DLBCL) is the most common type of non-Hodgkin's lymphoma in adults. It can occur any time between adolescence and old age, affecting 7-8 per 100,000 people annually in the United States, but incidence increases with age. Gene expression profiling has revealed two major types of DLBCL: germinal center B cell-like (GCB) and activated B cell-like (ABC). GCB DLBCL arises from secondary lymphoid organs, e.g., lymph nodes, where naive B cells do not stop dividing after an infection is cleared. ABC DLBCL is thought to begin with a subset of B cells that are ready to leave the germinal center and become plasma cells, i.e., plasmablast B cells, but the reality is more complex, as different forms of DLBCL arise throughout the B cell life cycle.

[0211] The different subtypes have a variable prognosis, with a 5-year survival rate of 60% for GCB DLBCL, but only 35% for ABC DLBCL. Each of the subtypes is characterized by differences in gene expression. In GCB DLBCL, the transcriptional repressor BCL6 is often overexpressed, whereas in ABC DLBCL, the NF-κB pathway is often found constitutively activated. There is also a third type of DLBCL, called type III, which is currently less well understood, but which is thought to have a gene expression profile that falls between the two major types.

[0212] Current diagnostic methods include excision biopsy of affected lymph nodes followed by immunohistochemistry (IHC). Currently, the treatment procedure for DLBCL is the same regardless of subtype. Because the pathogenesis, treatment response, and outcomes of various subtypes vary significantly, there is still a need to develop robust, noninvasive assays to distinguish subtypes to aid in the development of differentiated treatment strategies. Although much research has been done to find predictive and prognostic biomarkers for DLBCL, there is no consensus on a single test that can be used to distinguish subtypes.

[0213] Identifying EpiSwitch™ biomarkers capable of distinguishing different subtypes of DLBCL in the blood of DLBCL patients Using the EpiSwitch™ array platform, we interrogated DLBCL cell lines and blood samples to identify biomarkers that were absent in healthy control patients, and then confirmed these biomarkers in a 70 patient cohort consisting of 30 ABC, 30 GCB, and 10 healthy control samples.

[0214] EpiSwitch™ Arrays EpiSwitch™ custom arrays allow for the screening of thousands of possible CCSs with probes designed using pattern recognition software. The various long-range chromosomal interactions captured by EpiSwitch™ technology reflect the epigenetic regulatory framework imposed on the loci of interest and correspond to the individual distinct inputs from signaling pathways that contribute to the co-regulation of these loci. Overall, the combination of distinct inputs regulates gene expression. The identification of aberrant or distinct chromosomal conformation signatures under specific physiological conditions provides critical evidence for the differential contribution to deregulation before all input signals are integrated into a gene expression profile.

[0215] Using data from several sources, we selected 98 loci for analysis in our proprietary software and tested 13,332 potential chromosome conformation probes. Looking at one locus does not equate to looking at one marker, as a specific locus may have one or more higher order epigenetic chromosome conformation markers or no higher order epigenetic chromosome conformation markers. After production, cell lines and blood samples from DLBCL patients and healthy controls were processed, labeled, and hybridized to the arrays using the EpiSwitch protocol.

[0216] Diagnostic Development Samples Sixteen cell lines were used, corresponding to different subtypes and with different levels of confidence in the subclassification. Cell lines with the most distinct ABC and GCB subtypes were used for the analysis. In addition, blood samples from four DLBCL patients and 11 healthy controls were used. After biomarker identification in part 1, 60 samples, consisting of 30 ABC and 30 GCB blood samples, were provided to OBD and fully characterized by Fluidigm testing, complemented by 10 healthy control samples provided by OBD.

[0217] result Array analysis Seventy-two chromosomal signature sites from the microarray were selected for screening based on two criteria: their ability to stratify between ABC and GCB cells (highABC_highGCB), and / or Low CV values ​​(median of 5 arrays analyzed, high ABC vs. high GCB, DLBCL1 vs. healthy controls, DLBCL2 vs. healthy controls, DLBCL3 vs. healthy controls and DLBCL4 vs. healthy controls)

[0218] Conversion of Arrays to the EpiSwitch™ PCR Platform After analysis of the sequences surrounding the probes of interest from the array, 69 sets of primers were designed to match the chromosomal signature sites, which were then tested on pooled DLBCL blood samples, 49 of which met the OBD criteria for PCR products to be used in the assay.

[0219] Each of these 49 potential markers was then tested on six DLBCL cell lines, three of which were ABC and three of which were GCB. The cell lines used were the ones with the most confidence, ABC or GCB, since the same classification was found using several different identification methods. This allowed us to select the markers that were most useful for distinguishing between ABC and GCB cell subtypes. 28 EpiSwitch™ markers were identified for use in the PCR platform, which were consistent with the EpiSwitch™ microarray results. In addition, the potential markers were also tested on four DLBCL patients and pooled healthy controls to identify those present in DLBCL patients but absent in healthy controls. 21 of the 28 EpiSwitch™ markers were absent in the healthy control samples, but present in the DLBCL samples, such that they could be used as markers for DLBCL and for subclassification.

[0220] Sample Test The 21 markers that were successfully translated to the EpiSwitch™ PCR platform were then tested across a cohort of 70 patient blood samples. Initially, each marker was tested on 6 new ABC samples and 6 new GCB samples, narrowing the 21 marker set to 10 markers that showed the greatest difference. These 10 markers were then tested on the remaining 24 ABC, 24 GCB, and 10 healthy control samples.

[0221] Each of the markers was then subjected to an analysis of its ability to distinguish subgroups, its collinearity with other markers, and also its ability to distinguish healthy from DLBCL. Six subsets of markers were identified that provided the greatest possible information, these being markers at the ANXA11 IFNAR, MAP3K7, MEF2B, NFATc1, and TNFRS13C loci. Figure 3 shows the ability of these markers to distinguish different groups of samples on a PCA plot. This six-marker panel was able to clearly distinguish healthy control patients from DLBCL patients, which is a key feature of a blood-based assay for DLBCL.

[0222] Figure 3 shows a PCA plot of 60 DLBCL and 10 healthy patients based on the six EpiSwitch™ marker binary data. Samples were characterized as ABC or GCB subtypes by Fluidigm data, and healthy controls are also shown.

[0223] Classification: Identification of ABC and GCB subtypes within a DLBCL patient cohort (60 samples) Classification was performed using a logistic regression classifier with 5-fold cross-validation and the following results were achieved: ABC subtype: 83.3% (95%CI -65.3%~94.3%) GCB subtype: 83.3% (95%CI -65.3% to 94.3%)

[0224] The resulting 6-marker logistic classifier model was further tested in 50 permutations of a dataset of 60 patients. Data were randomized each time and accuracy statistics were calculated using the ROC curve. The area under the curve (AUC) of 0.802 and p-value of 0.0000037 (H0 = "AUC equals 0.5") indicate that the model is accurate and efficiently implemented.

[0225] conclusion This study demonstrated the power of EpiSwitch™ technology to provide answers to challenging clinical questions, particularly the distinction between ABC and GCB subtypes of DLBCL. Using high-throughput array methods, and conversion to a simple, cost-effective PCR platform, over 13,000 potential CCSs were tested and refined into a 6-marker panel to differentiate DLBCL subtypes. This panel was able to distinguish DLBCL patients from healthy controls and correctly predict subtype 83.3% of the time. The study also had >80% concordance for class assignment between EpiSwitch™ (whole blood-based), LPS (cells, tissues of origin), and Fluidigm (cells, tissues of origin).

[0226] EpiSwitch™ technology detects changes in long-range gene-gene interactions, i.e. chromosomal conformation signatures, which result in altered regulation of the epigenetic state and expression mode of key genes involved in disease pathogenesis. The diagnostic procedure based on EpiSwitch™ technology is a simple and rapid technique that can be transferred to other laboratories. The test consists of several molecular biology reactions followed by detection by nested PCR. The test does not require complex procedures and can be performed in any laboratory that performs PCR-based assays.

[0227] Example 3 Further studies were performed in dogs. One objective was to investigate markers to aid in the early diagnosis of suspected lymphoma to inform veterinarians on the requirements to perform follow-up biopsies. In this study, the top 75 EpiSwitch microarray DLBCL markers (previously identified) were translated from the human genome build (Grch37) to the current canine genome. A total of 38 canine samples (consisting of 19 patients with possible lymphoma and 19 matched control samples) were screened using all 75 DLBCL markers. To perform this study, the following was performed: - Based on 75 human DLBCL markers (associated with specific genes), orthologues in the dog genome (CanFam3.1) were identified and loci were extracted from Biomart. - EpiSwitch™ software was implemented to identify potential interactions at these loci. -Added primer design software and other filters to reduce the list to 75 markers for study. The studies and results are shown in Figures 6-16 and Tables 8 and 9.

[0228] Example 4. Further studies on prostate cancer Current diagnostic blood tests for prostate cancer (PCa) are unreliable for early stage disease, resulting in a large number of unnecessary prostate biopsies in men with benign disease and false relief of negative biopsies in men with PCa. As the 5-year survival rate in low-risk groups is over 95% and most men would benefit from minimally invasive treatments, predicting risk for PCa is crucial to make informed decisions about treatment options. Three-dimensional genomic architecture and chromosomal structure undergo early changes during tumorigenesis in both tumor and circulating cells and can serve as biomarkers for disease.

[0229] In this prospective study, we performed chromosome conformation screening for 14,241 chromosomal loops at loci of 425 cancer-associated genes in whole blood of newly diagnosed, treatment-naïve PCa patients (n=140) and non-cancer controls (n=96).

[0230] Data show that peripheral blood mononuclear cells (PBMCs) from PCa patients acquired specific chromosomal conformational changes at the loci of the ETS1, MAP3K14, SLC22A3, and CASP2 genes. Blinded testing in an independent validation cohort yielded PCa detection with 80% sensitivity and 80% specificity. Further analysis between PCa risk groups yielded a prognostic validation set consisting of the genes BMP6, ERG, MSR1, MUC1, ACAT1, and DAPK1 for high-risk category 3 vs. low-risk category 1, and the genes HSD3B2, VEGFC, APAF1, MUC1, ACAT1, and DAPK1 for high-risk category 3 vs. intermediate-risk category 2, which had high similarity to the conformation in primary prostate tumors. These sets achieved a sensitivity of 80% and a specificity of 92% for stratifying high-risk category 3 versus low-risk category 1, and a sensitivity of 84% and a specificity of 88% for stratifying high-risk category 3 versus intermediate-risk category 2 disease.

[0231] Results demonstrate specific chromosomal conformations in the blood of PCa patients, which allow for the diagnosis and prognosis of PCa with high sensitivity and specificity. These conformations are shared between PBMCs and primary tumors. These epigenetic signatures may lead to the development of blood-based diagnostic and prognostic tests for PCa.

[0232] introduction In Western countries, prostate cancer (PCa) is now the most commonly diagnosed non-skin cancer in men and the second leading cause of cancer-related death. Many men aged 30 years have histological evidence of PCa, most of which is microscopic and probably will not show clinical symptoms. Prostate-specific antigen (PSA), invasive needle biopsy, Gleason score, and stage are used for diagnosis and prognosis. In a large multicenter trial of 2,299 patients, a 12-site biopsy scheme was superior to all other schemes, with an overall PCa detection rate of only 44.4%.

[0233] The only available blood test for PCa in widespread clinical use involves measuring circulating levels of PSA (sensitivity of 21% and specificity of 91%), although PSA levels can also be increased by prostate size, benign prostatic hyperplasia, and prostatitis. At the current cutoff limit of 4.0 ng / ml, only 20% of all PCa patients are detected. In early PCa, PSA testing is not specific enough to distinguish between early invasive cancer and occult non-lethal tumors, which may remain asymptomatic throughout a man's life. In advanced PCa, PSA kinetics are used as a clinical surrogate endpoint for outcome. However, while they provide a general prognosis, they lack individual specificity. A number of more specific blood tests are available for PCa detection, including the 4K blood test (AUC 0.8) and the PHI blood test (sensitivity of 90%, specificity of 17%). PSA levels, stage, and Gleason score are used to establish the severity of PCa and stratify patients into risk groups. To date, there are no prognostic blood tests that make it possible to distinguish between low-risk and high-risk PCa.

[0234] There are multiple genetic alterations associated with PCa, including mutations in the tumor suppressor genes p53 (up to 64% of tumors), p21 (up to 55%), p73 and MMAC1 / PTEN, but these mutations do not account for all observed effects on gene regulation. Epigenetic mechanisms involving dynamic and multi-layered chromosomal loop interactions are powerful regulators of gene expression. Chromosome Conformation Capture (3C) technology makes it possible to record these signatures. In this study, we used the EpiSwitch™ assay to screen, define and evaluate specific chromosomal conformations in the blood of PCa patients and to identify loci that may serve as diagnostic and prognostic markers.

[0235] method A total of 140 PCa patients and 96 controls were recruited in two cohorts. Cohort 1: Men with (n=105) or without (n=77) PCa diagnosis attending a urology clinic were prospectively recruited from October 2010 to September 2013. Cohort 2: Patient samples (19 controls and 35 PCa) were obtained from the United States. At the time of recruitment, a single blood sample (5 ml) was collected from the PCa patients into a BD Vacutainer® plastic EDTA tube using current practice of needle and blood collection methods. Blood samples were passively frozen and stored at -80°C until processed. Prostate tumor samples were obtained from previously recruited patients (n=5) who subsequently underwent radical prostatectomy. Patient clinical characteristics are shown in Table 17.

[0236] The primary endpoint of this study was to detect chromosomal conformational changes in PBMCs from PCa patients compared with controls. Therefore, all treatment-naïve PCa patients were eligible for this study, regardless of grade, stage, and PSA level. Patients who had received prior chemotherapy or had other cancers were excluded from the study. PCa diagnosis was established according to clinical routine, and patients were assigned to appropriate treatment. For the prognostic study (secondary endpoint), patients were stratified according to the relevant NCCN risk groups (Table 10). No follow-up studies were performed.

[0237] Based on the preliminary findings for melanoma, a priori power analysis was performed using the pwr.t.test function in the R package pwd. The test suggested that 15 patients per group would be sufficient to detect correlations between variables (β = 5% probability type II error, significance level; 95% power; 50% confidence interval and 40% standard deviation).

[0238] The EpiSwitch™ technology platform combines high-resolution 3C results with regression analysis and machine learning algorithms to advance disease classification. To select epigenetic biomarkers capable of diagnosing cancer, samples from patients with cancer compared to healthy (control) samples were screened for statistically significant differences in conditional and stable profiles of genomic architecture. The assay is performed on whole blood samples by first fixing the chromatin with formaldehyde to capture associations within the chromatin. The fixed chromatin is then digested into fragments with TaqI restriction enzyme and the DNA strands are joined to favor crosslinked fragments. The crosslinks are reversed and a polymerase chain reaction (PCR) is performed using primers previously established by the EpiSwitch™ software. EpiSwitch™ was used on blood samples in a three-step process to identify, evaluate, and validate statistically significant differences in chromosome conformation between PCa patients and healthy controls (Figure 17). In a first step, sequences from 425 manually curated PCa-associated genes (obtained from a public database (www.ensembl.org)) were used as templates for this computational probabilistic identification of regulatory signals involved in chromatin interactions (Table 18). A customized CGH Agilent microarray (8x60k) platform was designed to test technical and biological replicates for 14,241 potential chromosomal conformations across 425 loci. Eight PCa and eight control samples were competitively hybridized to the arrays, and differences in the presence or absence of each locus were defined by LIMMA linear modeling, followed by binary filtering and cluster analysis. This initially revealed 53 chromosomal interactions with the ability to optimally discriminate between PCa patients and controls (Figure 17).

[0239] In the second evaluation phase, 53 biomarkers selected from the array analysis were converted into EpiSwitch™ PCR-based detection probes and used for multiple rounds of biomarker evaluation. PCR primers were selected according to their ability to distinguish PCa from healthy controls (n=6 in each group). The identity of the PCR products generated using nested primers was confirmed by direct sequencing. Thus, the 53 selected biomarkers were reduced to 15 markers after initial statistical analysis, ultimately resulting in a 5-marker signature (Table 11). This set of selected chromosome conformation signature biomarkers was then tested in a known cohort (n=49). Furthermore, the 5-marker signature developed from the EpiSwitch™ PCR evaluation of the array marker leads was tested in an independent blinded validation cohort in which 29 samples were combined with the known 49 samples initially tested (total of 78 samples). Principal component analysis was also used to determine abundance levels and identify potential outliers (Figure 18).

[0240] In a final step, the set of 5 markers was tested in a blinded, independent (n=20) cohort of blood samples to further validate the chromosomal conformation signature used to inform PCa diagnosis. Results were analyzed using Bayesian logistic modeling, p-value null hypothesis (Pr(N|z|) analysis, Fisher's exact P-test, and Glmnet (Table 12). The sample cohort size in the 5 marker signature study was gradually increased to allow selection of the optimal markers to discriminate between PCa samples and healthy controls. The cohort size was expanded to 95 PCa and 96 healthy control samples. Analysis and presentation of the data was performed according to the CONSORT recommendations. All measurements were performed in a blinded manner. STARD criteria were used to validate the analytical procedures. A similar three-step approach was followed for the identification of prognostic markers (Table 13).

[0241] Sequence-specific oligonucleotides were designed around the selected sites to screen for potential markers by nested PCR using Primer3. All PCR-amplified samples were visualized by electrophoresis on a LabChip GX using the LabChip DNA 1K Version2 kit (Perkin Elmer, Beaconsfield, UK), and internal DNA markers were loaded onto the DNA chip using fluorescent dyes according to the manufacturer's protocol. Fluorescence was detected by laser and reading of the electropherograms was converted into simulated bands on the gel image using the instrument software. The threshold set for bands to be considered positive was 30 fluorescence units or higher.

[0242] Primary tumor samples were obtained from biopsies of selected patients (n=5). Ground tissue samples were incubated in 0.125% collagenase at 37°C with gentle agitation for 30 min. Resuspended cells (250ul) were then centrifuged at 800g for 5 min at room temperature in a fixed-arm centrifuge, the supernatant was removed, and the pellet was resuspended in phosphate-buffered saline (PBS). Primary tumors and matched blood samples were analyzed for the presence of a set of 6 markers for categories 3 vs. 1 and 3 vs. 2 over a range of assay sensitivities (dilution factor 1:2). A score of 1 was assigned if a matching PCR band of the correct size was detected, and a score of 0 if no band was detected (Table 14).

[0243] A stepwise diagnostic biomarker discovery process using EpiSwitch™ technology was applied as described in the methods. A customized CGH Agilent microarray (8×60k) platform was designed to test technical and biological replicates for 14,241 potential chromosomal conformations across 425 loci in 8 PCa and 8 control samples (Figure 17) (Figure 18). Presence or absence of each locus was defined by LIMMA linear modeling, followed by binary filtering and cluster analysis. In a second evaluation step, nested PCR was used for the 53 selected biomarkers, further reducing them to 15 markers and finally to a 5-marker signature (Figure 17). This distinct chromosome conformation disease classification signature for PCa consisted of chromosome interactions at five genomic loci: ETS proto-oncogene 1, transcription factor (ETS1), mitogen-activated protein kinase kinase kinase 14 (MAP3K14), solute carrier family 22 member 3 (SLC22A3), and caspase 2 (CASP2) (Table 11). The genomic locations of specific chromosome loops in the ETS1, MAP3K14, SLC22A3, and CASP2 genes in the chromosome conformation signature (Table 11) were mapped to their relative chromosomes. Two genomic sites corresponding to the junctions of each chromosome conformation signature locus for the ETS1, MAP3K14, SLC22A3, and CASP2 genes were mapped to chromosome 11 from 128,260,682 to 128,537,926, chromosome 17 from 43,303,603 to 43,432,282, chromosome 6 from 160,744,233 to 160,944,757, and chromosome 7 from 142,935,233 to 143,008,163. Circos plots of the ETS1, MAP3K14, SLC22A3, and CASP2 chromosome conformation signature markers showing chromosome loops were generated. Principal component analysis for the five markers was used to determine abundance levels and identify potential outliers. This analysis was applied to 78 samples including the two groups.The first group of 49 known samples (24 PCa and 25 healthy controls) was combined with a second group of 29 samples, including 24 PCa samples and 5 healthy control samples (Figure 18). A final training set was constructed using 95 PCa and 96 control samples, which were then tested in an independent blinded validation cohort of 20 samples (10 control and 10 PCa). The sensitivity and specificity for PCa detection using chromosomal interactions at the five genomic loci were 80% (CI 44.39%-97.48%) and 80% (CI 44.39%-97.48%), respectively (Table 12).

[0244] To select epigenetic biomarkers capable of stratifying PCa, samples from PCa patients classified into risk group categories 1-3 (low, medium, and high, respectively, Table 10) were screened for statistically significant differences in conditional and stable profiles of genomic architecture. EpiSwitch™ was used on blood samples in a three-step process to identify, evaluate, and validate statistically significant differences in chromosomal conformations between PCa patients at different stages of the disease (Figure 17). In the first step, the array used covered 425 loci and the test probes were for a total of 14,241 potential chromosomal conformations. Patients in high-risk PCa category 3 were compared with those in low-risk category 1 or medium-risk category 2. In total, 181 potential stratifying marker reads for PCR evaluation were identified using enrichment statistics (Table 19). The top 70 markers were then carried forward to the next stage of PCR detection to further evaluate the stratification of patient samples in high-risk category 3 vs. low-risk category 1, and finally a 6-marker set for high-risk category 3 vs. low-risk category 1 was established (Table 13). The best markers were identified using chi-square and then incorporated into the classifier for the test sets of category 1 (n=21) and category 3 (n=19). Independent cohorts of category 1 (n=21) and category 3 (n=6) that were not used for marker reduction were then used for the first round of blinded validation. Similarly, the 6-marker set was evaluated for high-risk category 3 vs. medium-risk category 2 for the test sets of category 3 and category 2, containing 25 and 19 samples, respectively. Independent cohorts of category 2 and category 3 (n=6 in each group) that were not used for marker reduction were then used for the first round of blinded validation.

[0245] In a final step, to further validate the chromosomal conformation signature used to inform PCa prognosis, the 6-marker set for high-risk category 3 versus low-risk category 1 was tested in a larger, more representative cohort. The original blinded cohort was expanded to 67 samples, including 40 samples used for marker reduction (Table 15). Similarly, the 6-marker set for high-risk category 3 versus intermediate-risk category 2 was tested in a larger, more representative cohort. The original blinded cohort was expanded to 43 samples (Table 16).

[0246] A six-marker set for category 3 versus category 1 was established. This set included bone morphogenetic protein 6 (BMP6), ETS transcription factor ERG (ERG), macrophage scavenger receptor 1 (MSR1), mucin 1 (MUC1), acetyl-CoA acetyltransferase 1 (ACAT1), and cell death-associated protein kinase 1 (DAPK1) genes (Table 13). Six biomarkers for high-risk category 3 versus intermediate-risk category 2 were identified, including hydroxy-delta-5-steroid dehydrogenase, 3 beta- and steroid delta-isomerase 2 (HSD3B2), vascular endothelial growth factor C (VEGFC), apoptotic peptidase activating factor 1 (APAF1), MUC1, ACAT1, and DAPK1. Notably, the last three biomarkers (MUC1, ACAT1, and DAPK1) were common between category 1 versus 3 and 3 versus 2 (Table 13). Stratification of PCa into high-risk category 3 versus low-risk category 1 using chromosomal interactions at six genomic loci demonstrated a sensitivity of 80% (CI 59.30%-93.17%) and a specificity of 92% (CI 80.52%-98.50%) in a blinded cohort of 67 samples (Table 15). Similarly, a six-marker set for high-risk category 3 versus intermediate-risk category 2 was tested in a larger, more representative cohort of 43 samples demonstrating a sensitivity of 84% (CI 63.92%-95.46%) and a specificity of 88% (CI 65.29%-98.62%) (Table 16).

[0247] Using five matched peripheral blood and primary tumor samples, epigenetic markers identified in peripheral circulation (Table 13) were compared with tumor tissue. Results showed that a number of deregulated markers detected in blood as part of the category 1 vs. 3 and category 2 vs. 3 stratification signatures could be detected in tumor tissue (Table 14). This demonstrates that chromosomal interactions that can be detected systemically could be detected under the same conditions at the primary site of tumor formation.

[0248] Timely diagnosis of prostate cancer is important to reduce mortality. European randomized studies of screening for PCa have shown that men who undergo regular PSA screening have a significantly reduced PCa mortality rate. However, comprehensive screening leads to overdiagnosis of clinically insignificant disease, and new minimally invasive tests that can distinguish low-risk from high-risk disease are urgently needed.

[0249] Epigenetic analysis approaches offer a potentially powerful means to address this need. The binary nature of the test (chromosomal loops are present or absent) and the enormous combinatorial power (>1010 combinations are possible with ~50,000 loops screened) may make it possible to create signatures that precisely fit clinically well-defined criteria; in PCa, discriminating between low-risk vs. high-risk disease, identifying small but aggressive tumors, and determining the most appropriate treatment options. Moreover, epigenetic changes are known to appear early in tumorigenesis, making them useful for both diagnosis and prognosis.

[0250] In this study, we identified and validated chromosome conformation as a distinctive biomarker of a noninvasive blood-based epigenetic signature for PCa. The data demonstrate the presence of stable chromatin loops at the loci of the genes ETS1, MAP3K14, SLC22A3, and CASP2, which are present exclusively in PCa patients (Table 11). Validation of these markers in an independent set of 20 blinded samples showed a sensitivity of 80% and a specificity of 80% (Table 12), which is remarkable for a PCa blood test. Interestingly, the expression of some of these genes has already been linked to cancer pathophysiology. ETS1 is a member of the ETS transcription factor family. Prostate tumors that overexpress ETS1 are associated with increased cell migration, invasion, and induction of epithelial-mesenchymal transition. MAP3K14, also known as nuclear factor kappa beta (NF-kβ)-inducing kinase (NIK), is a member of the MAP3K family (or MEKK). Physiologically, MAP3K14 / NIK can activate non-canonical NF-kβ signaling and induce canonical NF-kβ signaling, especially when MAP3K14 / NIK is overexpressed. A novel role for MAP3K14 / NIK to control mitochondrial dynamics to promote tumor cell invasion has been described. SLC22A3 (also known as organic cation transporter 3 (OCT3)) is a member of the SLC family of membrane transport proteins. SLC22A3 expression is associated with PCa progression. CASP2 is a member of the group of caspase activation and recruitment domains. Physiologically, CASP2 can function as an endogenous repressor of autophagy. Two of the identified genes (SLC22A3 and CASP2) were previously shown to be inversely correlated with cancer progression. Importantly, the presence of chromatin loops can have an uncertain impact on gene expression.

[0251] To screen for PCa prognostic markers, EpiSwitch™ custom arrays were performed to analyze competitive hybridization of DNA from peripheral blood of patients with low-risk PCa (classification 1) and high-risk PCa (classification 3). A set of 6 markers was identified for high-risk category 3 versus low-risk category 1, including BMP6, ERG, MSR1, MUC1, ACAT1, and DAPK1. Six biomarkers were identified for high-risk category 3 versus intermediate-risk category 2, including HSD3B2, VEGFC, APAF1, MUC1, ACAT1, and DAPK1. Three of these biomarkers (MUC1, ACAT1, and DAPK1) were shared between these sets. The data show a high concordance between chromosome conformations in primary tumors and blood of matched PCa patients at stages 1 and 3 (Table 14). The prognostic significance and diagnostic value of some of these genes have been previously suggested. BMP6 plays an important role in PCa bone metastasis. In addition to ETS1, ERG is another member of the ETS family of transcription factors. Overwhelming evidence suggests that ERG is involved in several processes related to PCa progression, including metastasis, epithelial-mesenchymal transition, epigenetic reprogramming, and inflammation. MSR1 may confer moderate risk for PCa. MUC1 is a membrane-bound glycoprotein that belongs to the mucin family. High expression of MUC1 in advanced PCa is associated with adverse clinicopathological tumor characteristics and poor outcome. ACAT1 expression is elevated in high-grade and advanced PCa and serves as an indicator of reduced biochemical recurrence-free survival. DAPK1 may function as a tumor suppressor or as an oncogenic molecule for different cells. HSD3B2 plays a key role in steroid hormone biosynthesis and is upregulated in a relevant fraction of PCa characterized by adverse tumor phenotype, increased androgen receptor signaling, and early biochemical recurrence. VEGFC is a member of the VEGF family, and its increased expression is associated with lymph node metastasis in PCa specimens. A comprehensive biochemical approach has described APAF1 as the core of the apoptosome.

[0252] Despite the identification of these loci, the mechanisms of cancer-associated epigenetic changes in PBMCs remain unidentified. However, interactions can be detected systemically and under the same conditions at the primary site of tumor formation (Table 14). Therefore, to be able to measure changes, chromatin conformation in PBMCs must be dictated by external factors, presumably those produced by cells of PCa tumors. It is known that a significant proportion of chromosome conformation is controlled by non-coding RNAs that control tumor-specific conformations. Tumor cells have been shown to secrete non-coding RNAs that can be invaginated by neighboring or circulating cells and alter their chromosome conformation, which in this case could be regulatory factors. Although RNA detection as a biomarker remains very challenging (low stability, background drift, ongoing foundation for statistical stratification analysis), chromosome conformation signatures offer well-recognized stable binary advantages for the use of biomarker targets, especially when tested in the nucleus, because circulating DNA present in plasma does not retain the topological structure of the 3D conformation present in intact cell nuclei. It is important to mention that looking at one locus does not equal looking at one marker, since there may be multiple chromosome conformations present, representing parallel pathways of epigenetic regulation across the locus of interest. One of the key challenges in the current clinical practice of PCa diagnosis is the time it takes to make a definitive diagnosis. So far, there is no single definitive test for PCa. High levels of PSA lead the patient on a long journey of uncertainty, where they undergo a magnetic resonance imaging scan followed by a biopsy, if necessary. Although biopsy is more reliable than PSA testing, it is still the main procedure where missing cancer lesions can be problematic. The five-set biomarker panel described here is based on a relatively inexpensive and well-established molecular biology technique (PCR). Samples are based on biofluids that are easy to collect and provide clinicians with a rapidly available clinical readout within a few hours.This in turn will result in significant time and cost savings and aid in valuable diagnostic decisions that fill a gap in current protocols for the definitive diagnosis of PCa.

[0253] Predicting the risk of PCa is crucial to make informed decisions regarding treatment options. The 5-year survival rate in the low-risk group is >95%, and most men would benefit from minimally invasive treatments. Currently, PCa risk stratification is based on a combined assessment of circulating PSA, tumor grade (from biopsy), and tumor stage (from imaging findings). The ability to derive similar information using a simple blood test would allow for significant reductions in costs and speed up the diagnostic process. Of particular importance in PCa treatment is the identification of a small number of tumors that initially present as low risk but subsequently progress to high risk. This subset would therefore benefit from faster and more radical interventions.

[0254] In conclusion, we have identified here a subset of chromosome conformations in patient PBMCs that are highly indicative of the presence and prognosis of PCa. These signatures offer great potential for the development of rapid diagnostic and prognostic blood tests for PCa, which significantly exceed the specificity of the currently used PSA test. Preferred markers and combinations include: - ETS1, MAP3K14, SLC22A3, and CASP2. This is diagnosed by the following nested PCR markers: - BMP6, ERG, MSR1, MUC1, ACAT1, and DAPK1. This is the prognostic signature (high-risk category 3 vs. low-risk category 1 by nested PCR markers): - HSD3B2, VEGFC, APAF1, MUC1, ACAT1, and DAPK1. This is prognostic (high risk category 3 vs. intermediate risk category 2)

[0255] Example 5. Further studies on DLBLC Diffuse large B-cell lymphoma (DLBCL) is a heterogeneous hematological cancer, but can be broadly classified into two major subtypes: germinal center B-cell-like (GCB) and activated B-cell-like (ABC). GCB and ABC subtypes have widely different clinical courses, with ABC having a much worse survival prognosis. It has been observed that patients with different subtypes also respond differently to therapeutic interventions, and in fact, some have argued that ABC and GCB can be considered completely separate diseases. Because of this variability in response to treatment, having an assay to determine DLBCL subtype has important implications in guiding clinical approaches to the use of existing therapies, as well as in the development of new drugs. The current gold standard assay for subtyping DLBCL uses gene expression profiling on formalin-fixed paraffin-embedded (FFPE) tissue to determine the "cell of origin" and therefore the disease subtype. However, this approach has several significant clinical limitations in that it 1) requires biopsy, 2) requires complex, costly and time-consuming analytical approaches, and 3) does not classify all DLBCL patients.

[0256] Here, we employed an epigenomic approach and developed a blood-based chromosome conformation signature (CCS) for identifying DLBCL subtypes. An iterative approach using clinical samples from 118 DLBCL patients was employed to define a panel of six markers (DLBCL-CCS) to subtype the disease. The performance of DLBCL-CCS was then compared to conventional gene expression profiling (GEX) from FFPE tissues.

[0257] DLBCL-CCS was accurate in classifying ABC and GCB in samples of known status, providing identical calls in 100% (60 / 60) of samples in the discovery cohort used to develop the classifier. Also, in the evaluation cohort, DLBCL-CCS was able to make a DLBCL subtype call in 100% (58 / 58) of samples with intermediate subtype (type III) as defined by GEX analysis. Most importantly, when these patients were followed longitudinally throughout the course of their disease, EpiSwitch™-related calls tracked better with known patterns of survival for ABC and GCB subtypes.

[0258] This study provides early indications that a simple, accurate, cost-effective and clinically adoptable blood-based diagnostic for identifying DLBCL subtypes is possible.

[0259] background Diffuse large B-cell lymphoma (DLBCL) is the most common type of hematological cancer and many studies using different methodologies have demonstrated that it is genetically and biologically heterogeneous. The two major DLBCL molecular subtypes are germinal center B cell-like (GCB) and activated B cell-like (ABC), although more detailed definitions of molecular subtypes have also been proposed. These two major subtypes have a high degree of clinical relevance, as it has been observed that they have dramatically different disease courses and that the ABC subtype has a much worse survival prognosis. Perhaps more importantly, as novel investigational agents for treating GCB and ABC (or non-GCB) subtypes are evaluated in the clinical setting, and due to the historical observation of poor overall response rates in unselected patients, there is a pressing need to identify patient subtypes prior to the initiation of treatment. Historically, DLBCL subtypes are determined by identifying the "cell of origin" (COO). The original COO classification was based on the observed similarity of DLBCL gene expression to activated peripheral blood B cells or normal germinal center B cells by hierarchical clustering analysis (3). This COO classification by whole genome expression profiling (GEP) divided DLBCL into activated B cell-like (ABC), germinal center B cell-like (GCB), and type III (unclassified) subtypes, with ABC-DLBCL characterized by poor prognosis and constitutive NF-kB activation. In their prospective study, Wright et al. identified 27 genes whose expression was most discriminative between ABC and GCB-DLBCL and developed the Linear Prediction Score (LPS) algorithm for COO classification. These original studies were based entirely on retrospective investigations of fresh frozen (FF) lymphoma tissues. The main challenge for the application of this COO classification in clinical practice was the establishment of a robust clinical assay suitable for routine formalin-fixed paraffin-embedded (FFPE) diagnostic biopsies.Several studies have also investigated the feasibility of COO classification of DLBCL using FFPE tissues by quantitative measurements of mRNA expression, including quantitative nuclease protection assays, GEP using the Affymetrix HG U133 Plus 2.0 platform or Illumina whole genome DASLassay, and NanoString Lymphoma Subtyping Test (LST) technology. Several immunohistochemistry (IHC)-based algorithms have also been investigated to summarize COO classification by GEP. Overall, these studies have demonstrated high reliability of COO classification of DLBCL using FFPE tissues and robust separation of overall survival between ABC and GCB subtypes, but suffer from reproducibility issues, especially lack of inter-assay agreement. Furthermore, any IHC-based measurement requires baseline tissue, which is not always available, and the current turnaround time from sample collection to assay readout is long, making it difficult to implement in clinical practice.

[0260] Among the approaches that have been used historically to subtype DLBCL, one method for COO assessment uses an assay that measures the expression of 27 genes from FFPE tissue by quantitative reverse transcription PCR (qRT-PCR) using the Fluidigm BioMark HD system. Although this methodology has several advantages over existing techniques, the approach still faces several major obstacles that limit its clinical application in that it 1) requires tissue biopsy and 2) relies on expensive, non-standard, and time-consuming testing procedures. Thus, the presence of a blood-based assay advances the field by providing a simple, reliable, and cost-effective method for DLBCL subtyping enriched for clinical applicability.

[0261] In this study, we used a novel blood-based assay to determine COO classification in DLBCL patients by focusing on the detection of alterations in genome architecture. As part of the epigenetic regulatory framework, genomic regions can modify their three-dimensional structure as a way to functionally control gene expression. The outcome of this regulatory mechanism is the formation of chromatin loops at distinct genomic loci. The presence or absence of these loops can be empirically measured using chromosome conformation capture (3C). Multiple genomic regions contribute to regulating epistasis through the formation of stable, conditional, long-range chromosomal interactions. The collective measurement of chromosome conformation at multiple genomic loci results in a chromosome conformation signature (CCS), or molecular barcode that reflects the genomic response to its external environment. To detect, screen, and monitor CCS, we utilized the EpiSwitch platform, an established, high-resolution, and high-throughput methodology for detecting CCS. Based on 3C, the EpiSwitch platform has been developed to assess chromatin structure changes at defined loci, as well as long-range non-coding cis- and trans-regulatory interactions. Among the advantages of using EpiSwitch for patient stratification are its binary nature, reproducibility, relatively low cost, rapid turnaround time (samples can be processed within 24 hours), requirement of only small amounts of blood (~50 mL), and compliance with FDA standards for PCR-based detection methods. Thus, chromosome conformation provides a stable binary readout of cell state and represents an emerging class of biomarker.

[0262] Here, we developed a blood-based diagnostic test for DLBCL C0O subtyping using an approach based on the assessment of chromosomal architectural alterations. We hypothesized that investigation of genomic architectural alterations in blood samples from DLBCL patients could offer an alternative to tissue-based COO classification approaches and provide a novel, non-invasive, and more clinically applicable methodology to guide clinical decision-making and study design.

[0263] In total, 118 DLBCL patients with known COO subtypes and 10 healthy controls (HC) were used in this study. Samples were a subset of samples collected in a phase III, randomized, placebo-controlled trial of rituximab and bevacizumab in advanced non-Hodgkin lymphoma. Briefly, adult patients aged 18 years or older with newly diagnosed CD20-positive DLBCL were randomized to R-CHOP or R-CHOP plus bevacizumab (RA-CHOP). Blood samples collected from 60 DLBCL patients were used as a development cohort to identify, evaluate, and refine CCS biomarker leads. All patients from this cohort were classified as high / strong GCB (30) or ABC (30) with high subtype-specific LPS (linear predictive score). The remaining 58 DLBCL samples had intermediate LPS and were determined by Fluidigm testing as ABC, GCB, or unclassified (Figure 25). These patient samples were not used in the discovery and development of CCS biomarkers, but were used at a later stage to evaluate the resulting classifiers. Fluidigm testing was performed using tissue obtained from lymph nodes (either as punch biopsies or removed during surgery), and EpiSwitch analysis was performed using matched peripheral whole blood collected from patients before receiving treatment.

[0264] In addition to patient samples, 12 cell lines (6 ABC and 6 GCB) were also used in the initial stage of biomarker screening to identify a set of chromosome conformations that could optimally distinguish between ABC and GCB disease subtypes (Table 20). Cell lines were obtained from the American Type Culture Collection (ATCC), the German Collection of Microbial Cell Cultures (DSMZ), and the Japan Health Science Foundation (JHSF).

[0265] RNA was isolated and purified from pretreatment FFPE biopsies. DLBCL subtype was determined by applying the algorithm of Wright et al. to expression data from a custom Fluidigm gene expression panel containing 27 genes of DLBCL subtype predictors. Validation of the COO assay by comparing Fluidigm qRT-PCR to Affymetrix data in a cohort of 15 non-study subjects revealed high correlation between qRT-PCR measurements from matched fresh frozen (FF) and FFPE samples across the 19 classification genes used. We also found high correlation between Affymetrix microarray and Fluidigm qRT-PCR measurements from the same FF samples. Classification gene weights calculated from qRT-PCR data from the Fluidigm COO assay were highly consistent with weights obtained from previous microarray data in an independent patient cohort. We observed high correlation (76% concordance) between LPS derived from Fluidigm assays, data in FFPE tumors, and LPS derived from Affymetrix microarray data in matched FF tissues in the technical registry cohort.

[0266] A pattern recognition algorithm was used to annotate the human genome for sites that may form long-range chromosomal conformations. The pattern recognition software works based on Bayesian modeling and provides a probabilistic score that a region is involved in long-range chromatin interactions. Sequences from 97 loci (Table 21) were processed through the pattern recognition software to generate a list of 13,322 chromosomal interactions that are most likely to be able to distinguish DLBCL subtypes. The initial screen involved an array-based comparison. 60-mer oligonucleotide probes were designed to match these potential interactions and uploaded to the Agilent SureDesign website as custom arrays. Each probe was present in quadruplicates on the EpiSwitch microarray. To subsequently evaluate potential CCS, nested PCR (EpiSwitch PCR) was performed using sequence-specific oligonucleotides designed using Primer3. Oligonucleotides were tested for specificity using oligonucleotide-specific BLAST.

[0267] The top 10 genomic loci identified as dysregulated in DLBCL were uploaded as protein lists into the Reactome Functional Interaction Network plugin in Cytoscape to generate a network of epigenetic dysregulation in DLBCL. The 10 loci were also uploaded into STRING (Search Tool for Retrieval of Interacting Genes / Protein DB) (https: / / string-db.org / ), a database containing over 9 million known and predicted protein-protein interactions. The main network (i.e., unconnected nodes were excluded) was generated by restricting to human interactions only. Top false discovery rate (FDR) corrected functional enrichments were identified by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. The top 10 genomic loci were also uploaded into the KEGG Pathway Database (https: / / www.genome.jp / kegg / pathway.html) to identify specific biological pathways that exhibit dysregulation in DLBCL.

[0268] Exact tests and Fisher's exact test (for categorical variables) were used to identify discriminatory markers. The level of statistical significance was set at p ≤ 0.05, and all tests were two-sided. A random forest classifier was used to evaluate the ability of EpiSwitch markers to identify DLBCL subtypes. Long-term survival analysis was performed by Kaplan-Meier analysis using the survival and survminer packages in R (38). Mean survival times were calculated using two-tailed t-tests.

[0269] A stepwise approach was taken to discover and validate a CCS biomarker panel capable of distinguishing DLBCL subtypes (Figure 19). As a first step in the discovery of EpiSwitch classifiers, 97 loci (Table 21) were selected and annotated for the predicted presence of chromosome conformation interaction sites and screened for their empirical presence using EpiSwitch CGH Agilent arrays. The annotated array design represented 13,322 chromosome interaction candidates, with an average of 99 distinct cis interactions tested at each locus (99 ± 64; mean ± SD). This discovery array was used to screen and identify a smaller pool of chromosome conformations capable of distinguishing the two major DLBCL subtypes. Samples used in this step were from GCB and ABC cell lines (Table 20), and whole blood from four classified DLBCL patients (2 GCB and 2 ABC) and four HCs. Cell lines were grouped into high ABC and GCB and low ABC and GCB based on gene expression analysis. Comparisons used on the array were: 1) DLBCL patients vs pooled HC, 2) pooled DLBCL samples vs pooled HC samples, 3) pooled high GCB cell lines vs pooled high ABC, and 4) pooled low ABC vs pooled low GCB cell lines.

[0270] From the array analysis, we identified 1,095 statistically significant chromosomal interactions that distinguished between ABC high and GCB cell lines and were present in blood samples from DLBCL patients but not in HCs. These were further reduced using a set of statistical filters to the top 293 interactions, of which 151 were associated with the ABC subtype and 143 with the GCB subtype. The top 72 interactions from either subtype (36 interactions for ABC and 36 interactions for GCB) were selected for further refinement using the EpiSwitch PCR platform on 60 classified DLBCL patient samples. For all 118 DLBCL samples, an initial subtype classification was assigned based on the Wright algorithm, which calculates a linear prediction score (LPS) from the expression of a panel of 27 genes. Sixty samples were classified as either ABC or GBC and were used to develop the EpiSwitch classifier (the "discovery cohort"), and 58 samples had intermediate LPS scores and were used to evaluate the performance of the EpiSwitch classifier (the "evaluation cohort") (Figure 19).

[0271] The 72 interactions identified in the initial screen were narrowed down to a smaller pool using both the DLBCL patient samples in the discovery process and a second cohort of 60 DLBCL typed (30 ABC and 30 GCB) patient samples along with 12 HCs (Figure 19). The DLBCL subtype call made by the EpiSwitch assay was confirmed using the Fluidigm platform. Fluidigm gene expression analysis was performed on tissue biopsy samples, while whole blood from the same patients was used for the EpiSwitch PCR assay. In the first step of purification, the 72 chromosomal interactions identified in the initial screen were confirmed by PCR to be specific for DLBCL and absent from the HC samples. This was first tested in six untyped DLBCL samples and two HCs, resulting in the identification of 21 interactions specific for DLBCL. EpiSwitch PCR was then used to test 24 blood samples from typed DLBCL patient samples (12 ABC and 12 GCB) and the Fisher exact test was used to identify chromosomal interactions specific for DLBCL. This resulted in a set of 10 discriminatory chromosome conformation interactions capable of accurately distinguishing between ABC and GCB subtypes, which were further evaluated on blood samples from an additional set of 36 DLBCL samples (18 ABC and 18 GCB) (Figure 19).

[0272] To test the accuracy, performance, and robustness of the 10-marker panel, we used direct testing for feature selection on 80% of the complete sample cohort (48 samples in total: 24 ABC and 24 GCB), and the remaining 20% ​​(12 samples, 6 ABC and 6 GCB) were used for subsequent testing of the final selected CCS markers. The data was split 10 times and direct testing was performed on each of the splits using 80% of the training set of each split. The markers were then ranked using a composite p-value for the 10 markers across the 10 splits. This analysis identified six chromosome conformations at the IFNAR1, MAP3K7, STAT3, TNFRSF13B, MEF2B, and ANXA11 loci. Collectively, these six interactions formed the DLBCL chromosome conformation signature (DLBCL-CCS) (Figure 20).

[0273] The six markers in DLBCL-CCS were used to create a random forest classifier model and applied to classify the test set (12 samples, 6 ABC and 6 GCB) for each of the data splits in the discovery cohort of known disease subtypes. Principal component analysis (PCA) allowed the DLBCL-CCS classifier to separate ABC and GCB patients from healthy controls (Figure 26). The composite predicted probability for DLBCL-CCS is shown in Table 22 along with the odds ratios for each marker and the model generated using logistic regression. The model provided predicted probability scores for ABC and GCB ranging from 0.186 to 0.81 (0=ABC, 1=GCB). Probability cutoff values ​​for correct classification were set at 0.30 or less for ABC and 0.70 or more for GCB. A score of 0.30 or less had a true positive rate (sensitivity) of 100% (95% confidence interval [95% CI] 88.4-100%), while a score of 0.70 or more had a true negative response rate (specificity) of 96.7% (95% CI 82.8-99.9%). Using the DLBCL-CCS classifier, 60 of 60 patients (100%) were correctly classified as ABC or GCB compared with Fluidigm calls for subtyping (Figure 21A, Table 22). The AUC under the receiver operating characteristic (ROC) curve for the DLBCL-CCS classifier in this sample cohort was 1 (Figure 21B). Finally, the DLBCL subtype calls made by DLBCL-CCS were compared to long-term survival curves of patients with known disease subtypes. Patients called as ABC had significantly worse survival than patients called as GBC (Figure 21C).

[0274] We next evaluated the performance of the DLBCL-CCS evaluation cohort of 58 DLBCL patients with more intermediate LPS values. We applied DLBCL-CCS to assign these patients to DLBCL subtypes and compared the readouts to those produced by Fluidigm. DLBCL-CCS made a subtyping call for all 58 samples, while the Fluidigm assay made a subtyping call for 37 samples, leaving 21 as "unclassified" (Figure 22). Of the 37 samples for which subtype calls for both assays were available, 15 samples (40%) were called similarly by both assays (8 ABC and 7 GCB) (Figure 22). We then evaluated the performance of DLBCL subtype calls made by DLBCL-CCS and Fluidigm by comparing the subtype calls made at diagnosis to the long-term survival curves of type III patients. As shown in the Kaplan-Meier survival curves in Figure 23, the ABC / GBC call made by DLBCL-CCS was able to separate the two populations based on known survival trends in DLBCL, with the ABC subtype having a poor prognosis. In contrast, the ABC and GCB populations defined by Fluidigm were the opposite of what has been observed clinically, with samples classified as ABC having a longer survival time than those classified as GCB. Although not statistically significant, the subtype call made by DLBCL-CCS was consistent with previous clinical observations of survival differences between subtypes by hazard ratio analysis. A significant difference in mean survival time was found between the two methods. The mean survival times of patients classified as ABC and GCB by Fluidigm were 651 and 626 days, respectively (p=0.854), whereas the mean survival times of patients classified as ABC and GCB by the DLBCL-CCS assay were 550 and 801 days (p=0.017) (Figure 24).

[0275] To explore the association between the loci observed to be epigenetically dysregulated in this study and biological mechanisms previously reported to be associated with DLBCL, a series of network and pathway analyses were performed using the top 10 dysregulated loci as input. We first explored how these loci are biologically related by constructing a Reactome Functional Interaction Network in Cytoscape, which revealed a network centered around NFKB1, STAT3, and NFATC1. A similar picture emerged when 10 loci were used to construct the network using STRING DB, with the most connected hubs centered around NFKB1, STAT3, MAP3K7, and CD40. The top enriched GO terms for biological process were “positive regulation of transcription, DNA template”, the top enriched GO terms for molecular function were “transcriptional activator activity, RNA polymerase II transcriptional control region sequence-specific binding”, and “Toll-like receptor signaling pathway” was the most enriched KEGG pathway (Table 22). Mapping the top 10 loci to the KEGG Toll-like receptor signaling pathway, we found specific cascades associated with the production of proinflammatory cytokines and costimulatory molecules via the NF-kB and interferon-mediated JAK-STAT signaling cascades.

[0276] As differences in disease progression are observed for the different DLBCL subtypes, there is an urgent clinical need for a simple and reliable test that can distinguish between ABC and GBC disease subtypes. Given the aggressiveness of the disease, DLBCL requires immediate treatment. The two major subtypes have different clinical management paradigms, and with several therapies in development targeting specific subtypes, a fast and accurate disease diagnosis is crucial if clinical management depends on knowledge of the disease subtype. The field of COO classification in DLBCL has expanded from IHC-based methodologies to DNA microarrays, parallel quantitative reverse transcription PCR (qRT-PCR), and digital gene expression. The currently favored method is based on the identification of COO by GEP on FFPE tissues and suffers from several technical and logistical limitations that limit its wide adoption in clinical practice. Furthermore, there are many factors that affect the performance and reliability of COO classification by GEP on FFPE tissues, and these include the nature / quality of the lymphoma specimen, the experimental method for data collection, data normalization and transformation, the type of classifier used, and the probability cutoff used for subtype assignment. Finally, going from sample collection to final readout using the Fluidigm approach is a complex and time-consuming process with many steps that can result in variability in performance. All of these factors impact the overall turnaround time of the assay, limiting the method's ability to be used clinically to select patients for late-stage testing of novel DLBCL therapeutics in addition to diagnosing the disease and informing its treatment with existing agents. Thus, there is a need for a simple, minimally invasive, and reliable assay to differentiate DLBCL subtypes.

[0277] Using a stepwise discovery approach, we identified a 6-marker epigenetic biomarker panel, DLBCL-CCS, that could accurately distinguish DLBCL subtypes. When compared to the subtype results derived from the gene expression signature, there was perfect agreement, which was expected as these were the samples used to develop the classifier. The agreement between the two assays was low (only >40%) when applied to samples with intermediate LPS. This is likely expected as the classification methods have been noted to lack overall agreement in calling different DLBCL subtypes, and type III samples are a heterogeneous population that likely reflects a more intermediate biology to begin with. However, when we evaluated the predictive classification ability of the EpiSwitch Assay in type III DLBCL patients followed longitudinally as their disease progressed, the baseline prediction of disease subtype using the EpiSwitch Assay was good at predicting the actual disease subtype based on the survival curves observed in patients with unclassified disease. The observation that the epigenetic readout based on controlled 3D genomics used here is more consistent with actual clinical outcomes than the gold standard molecular approach based on transcription represents a potential advance in the management of DLBCL. This is also consistent with the systems biology evaluation of controlled 3D genomics as a molecular modality closely related to phenotypic differences in oncological conditions. We note that DLBCL operates on a biological continuum, with significant heterogeneity in disease biology between subtypes. By design, we set up DLBCL-CCS to classify type III samples into either ABC or GCB subtypes. GEX analysis identified type III samples as having intermediate subtype biology, and thus may represent a more heterogeneous population of patients.However, the overall observation that DLBCL-CCS is a better predictor of disease subtype as measured by clinical progression than using a GEX-based approach, and the fact that the EpiSwitch assay was able to make subtype calls on all samples, provide a first indication that this approach could be applied in the clinical setting to inform prognostic outlook, guide treatment decisions, and provide predictions for response to novel therapeutics currently in development.

[0278] In network analysis, NF-kB and STAT3 signaling cascades emerged as putative mediators distinguishing DLBCL subtypes. The role of NF-kB signaling in DLBCL has been previously studied, and indeed, one of the discriminatory features of the ABC subtype is the constitutive expression of NF-kB target genes, a mechanism hypothesized for the poor prognosis of these patients. Moreover, mutations causing constitutive signaling activation have been predominantly observed in the ABC subtype in several NF-kB pathway genes, including TNFAIP3 and MYD88.

[0279] In addition to validating known mechanisms of DLBCL, our network analysis identified novel potential targets for therapeutic intervention in DLBCL. For example, ANXA11, a calcium-regulated phospholipid-binding protein, is involved in other oncological diseases such as colorectal, gastric, and ovarian cancers, and may be a novel therapeutic intervention point in DLBCL.

[0280] One of the main clinical advantages of the approach to DLBCL subtyping described here is the simplified testing method and workflow. Traditional gold-standard subtyping by GEP can be performed using a variety of commercial platforms, but all broadly follow (and require) the following four-step approach: 1) obtaining a tissue biopsy, 2) preparing FFPE tissue sections, 3) gene expression analysis, and 4) algorithmic classification of subtypes. Obtaining a fine-needle tissue biopsy of an enlarged peripheral lymph node requires an inpatient visit to a clinical site and undergoing an invasive medical procedure that requires anesthesia. Once obtained, the fresh biopsy must be prepared for paraffin embedding. This is a multi-step process, but broadly involves immersion in a liquid fixative (such as formalin) for a time sufficient to allow penetration throughout the specimen, sequential dehydration through an ethanol gradient, and subsequent clearing in xylene, a toxic chemical. Finally, the biospecimen must be infiltrated with paraffin wax and allowed to cool and solidify, and then cut into micrometer sections using a microtome and mounted on laboratory slides. The entire process of going from fresh tissue to FFPE sections on slides can take several days. Then, to perform gene expression analysis, intrinsically labile RNA is extracted from the slide-mounted tissue sections and prepared for hybridization to microarrays according to the array manufacturer's specifications, a process that can take more than a day. Following microarray hybridization, a digital readout of relative gene expression levels is obtained and fed into a classification algorithm to determine DLBCL subtype. Overall, the process of going from a patient with suspected DLBCL to a subtype readout can take up to a week or more and involves many different experimental steps using expensive technologies, each of which can introduce experimental variability along the way. The approach described here dramatically reduces the time and number of steps from biofluid collection to subtype readout. Patients with suspected DLBCL can present to an outpatient clinic for regular small (~1 mL) blood draws.Fresh frozen blood can then be sent to a central accredited reference laboratory for analysis for the presence or absence of the chromosome conformations identified in this study, a process that uses an additional small volume (~50 mL) of whole blood as input along with specific PCR primer sets and reaction conditions to detect chromosome conformations using concise, routine PCR measurements within 24 hours of receiving the sample. The approach to DLBCL subtyping described here offers the added advantage that the potential for further refinement exists using the proposed methodology. In this study, the final readout of DLBCL-CCS was performed using a set of nested PCR reactions to detect the chromosome conformations that constitute the classifier. This PCR-based output can be further refined to utilize quantitative PCR as the readout and operate under the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) guidelines, designed to increase the reproducibility and reliability of experiments across reference laboratories and testing locations. Finally, the approach described here is adaptable to the evolving understanding of the disease itself, such as its various physiologically heterogeneous forms.

[0281] In conclusion, here we developed a robust complementary method for non-invasive COO assignment from whole blood samples using the EpiSwitch CCS readout. We demonstrated the clinical validity of this classification approach in a large cohort of DLBCL patients. The EpiSwitch platform has several attractive features as a biomarker modality with clinical utility. CCS has very high biochemical stability, can be detected using very small volumes of blood (typically ~50 μl), and utilizes established laboratory methods and standard PCR readouts (including MIQE-compliant qPCR) for detection. Finally, the rapid turnaround time of the EpiSwitch assay (~8-16 hours) compares favorably to >48 hours for the Fluidigm platform.

[0282] Example 6 Further studies on canine DLBCL Here, we used the EpiSwitch™ platform technology to evaluate chromosome conformation signatures (CCSs) as biomarkers for the detection of canine diffuse large B-cell lymphoma (DLBCL). We investigated whether established systemic liquid cytology biomarkers previously characterized in human DLBCL patients by EpiSwitch™ could be translated to dogs with homologous immune diseases. Orthologous sequence translation of human-to-canine CCSs was first validated and validated in a cohort of control and lymphoma dogs.

[0283] Blood samples were obtained from dogs with DLBCL and apparently healthy dogs. All dogs diagnosed with DLBCL were part of the LICKing Lymphoma Clinical Trial. Blood samples were obtained from each dog before starting treatment and at day +5 after the experimental intervention but before starting doxorubicin chemotherapy. EpiSwitch technology was used to monitor systemic epigenetic biomarkers for CCS.

[0284] An 11-marker classifier was generated from a pool of 75 EpiSwitch CCSs identified in human DLBCL in whole blood from 28 dogs, 14 diagnosed with DLBCL and 14 controls without apparent disease. Validation of the developed diagnostic markers was performed in a second cohort of 10 dogs, 5 with DLBCL and 5 controls. The classifier provided stratification for DLBCL vs. non-DLBCL with 80% accuracy, 80% sensitivity, 80% specificity, 80% positive predictive value (PPV) and 80% negative predictive value (NPV) in the second cohort.

[0285] The established EpiSwitch classifier contains robust systemic binary markers of epigenetic deregulation with features normally attributed to genetic markers, and the binary status of these classifiers is statistically significant for diagnosis.

[0286] [Table 5] TIFF2025032205000007.tif69157 TIFF2025032205000008.tif165157 TIFF2025032205000009.tif161157 TIFF2025032205000010.tif72157 TIFF2025032205000011.tif254157 TIFF2025032205000012.tif235157 TIFF2025032205000013.tif90157 TIFF2025032205000014.tif135157 TIFF2025032205000015.tif237157 TIFF2025032205000016.tif237157 TIFF2025032205000017.tif231157 TIFF2025032205000018.tif88157 TIFF2025032205000019.tif136157 TIFF2025032205000020.tif177157 TIFF2025032205000021.tif48157 TIFF2025032205000022.tif238157 TIFF2025032205000023.tif43157 TIFF2025032205000024.tif177157 TIFF2025032205000025.tif119157 TIFF2025032205000026.tif96157 TIFF2025032205000027.tif147157 TIFF2025032205000028.tif77157 TIFF2025032205000029.tif180157 TIFF2025032205000030.tif54157 TIFF2025032205000031.tif198157 TIFF2025032205000032.tif254145 TIFF2025032205000033.tif236157 TIFF2025032205000034.tif237157 TIFF2025032205000035.tif251151 TIFF2025032205000036.tif251157 TIFF2025032205000037.tif251154 TIFF2025032205000038.tif237157

[0287]

Table 6

[0288]

Table 7

[0289]

Table 8

[0290]

Table 9

[0291]

Table 10

[0292]

Table 11

[0293]

Table 12

[0294]

Table 13

[0295]

Table 14

[0296]

Table 15

[0297]

Table 16

[0298]

Table 17

[0299]

Table 18

[0300]

Table 19

[0301]

Table 20

[0302]

Table 21

[0303]

Table 22

[0304]

Table 23

[0305]

Table 24

[0306]

Table 25

Claims

1. 1. A process for determining a prognosis or diagnosis of prostate cancer, comprising: all of the chromosomal interactions (a) to (f) for determining said prognosis, or All of the chromosomal interactions (g) to (k) for identifying the diagnosis the process of detecting the presence or absence in an individual, The presence or absence of the chromosomal interaction is (i) cross-linking chromosomal regions of said individuals that are joined together by chromosomal interactions; (ii) subjecting the crosslinked region to scission; (iii) ligating the cross-linked and cleaved DNA ends to form a ligated nucleic acid; (iv) detecting the presence or absence of ligated nucleic acids corresponding to each chromosomal interaction. is detected by a method comprising: In the prognosis, the individual is classified as having Class 1, low-risk, indolent prostate cancer, or Class 3, high-risk, aggressive prostate cancer, and the chromosomal interaction is (a) a chromosomal interaction formed on chromosome 6 of the BMP6 gene between the region of position numbers 7733465-7733496 and the region of position numbers 7806286-7806316; and (b) a chromosomal interaction formed on chromosome 11 of the ACAT1 gene between the region of position numbers 107960135 to 107960166 and the region of position numbers 108018337 to 108018367; and (c) a chromosomal interaction formed on chromosome 21 of the ERG gene between the region of position numbers 39895678-39895708 and the region of position numbers 39991875-39991905; and (d) a chromosomal interaction formed on chromosome 8 of the MSR1 gene between the region of position numbers 16203284 to 16203315 and the region of position numbers 16400368 to 16400398; and (e) a chromosomal interaction formed on chromosome 1 of the MUC1 gene between the region of position numbers 155149955-155149986 and the region of position numbers 155191807-155191837; and (f) a chromosomal interaction formed on chromosome 9 of the DAPK1 gene between the region of position numbers 90073586 to 90073617 and the region of position numbers 90140806 to 90140836; and The diagnosis identifies the presence or absence of prostate cancer and the chromosomal interaction comprises: (g) a chromosomal interaction formed on chromosome 11 of the ETS1 gene between the region of position numbers 128419843 to 128419873 and the region of position numbers 128481262 to 128481292; and (h) a chromosomal interaction formed on chromosome 6 of the SLC22A3 gene between the region of position numbers 160805748 to 160805778 and the region of position numbers 160839018 to 160839048; and (i) a chromosomal interaction formed on chromosome 6 of the SLC22A3 gene between the region of position numbers 160805748 to 160805778 and the region of position numbers 160884099 to 160884129; and (j) a chromosomal interaction formed on chromosome 17 of the MAP3K14 gene between the region of position numbers 43364251 to 43364282 and the region of position numbers 43409961 to 43409991; and (K) a chromosomal interaction formed on the seventh chromosome of the CASP2 gene between the region of position numbers 142947138 to 142947169 and the region of position numbers 142963973 to 142964003; process.

2. 2. The process of claim 1, wherein detecting the presence or absence of the chromosomal interaction comprises specific detection of the ligated nucleic acid by quantitative PCR (qPCR) using primers capable of amplifying the ligated nucleic acid and a probe that binds to the ligated site during a PCR reaction.

3. The process of claim 2, wherein a fluorophore is covalently attached to the 5' end of said probe.

4. The process of claim 2, wherein a quencher is covalently attached to the 3' end of said probe.

5. The process of claim 2 , wherein the probe comprises a nucleic acid sequence between 10 and 40 nucleotide bases in length.

Images