The CORE of the matter: Identifying recurrent genomic regions to determine tumor phylogeny

June 24, 2013 by Stuart Mason Dambrot feature
CORE analysis of a simulated set of events. The red-hashed rectangles indicate the positions and event counts of the five recurrent regions used to simulate data. The dashed line gives the event count for 200 events assigned to recurrent regions and the solid line the joint event count for these 200 and additional 100 background events. The blue Pi shapes indicate the positions and CORE scores of the five significant (P < 0.05) cores. Copyright © PNAS, doi:10.1073/pnas.1306909110

(Medical Xpress)—Analyzing genome-wide data from organisms, tissues or cells generates lists of genomic intervals – continuous structures on a genomic sequence, such as a chromosome. An interval's meaning is dependent on its genomic context, and can identify a region of cancer-related DNA copy number variations (alterations of a genome's DNA that results in the cell having an abnormal number of copies of one or more sections of that DNA). Recently, scientists at Cold Spring Harbor Laboratory devised the computational biology method CORE (Cores Of Recurrent Events) that explains genome-wide data in terms of a small number of cores – that is, recurrent intervals – for studying tumor subpopulations and cancer type copy number aberrations.

Prof. Alexander Krasnitz discusses the research that he, Prof. Michael Wigler, graduate student Guoli Sun, and Dr. Peter Andrews conducted. "The main challenge in constructing our method was to find a suitable form of what we call explanation," Krasnitz tells Medical Xpress. The researchers' method utilizes a core to "explain" an observed interval event by assigning a measure of geometric association between the two. "Some of the forms we proposed and used identify recurrence of both broad and narrow genomic regions," Krasnitz adds, "allowing the method to work across a wide range of sizes of , from less than a single gene to entire ."

Krasnitz notes that CORE was implemented as a combinatorial optimization procedure that includes that had to be designed with great care. "In any study," he points out, "statistics is often the most treacherous part. Err in one direction, and you will miss important findings – but err in the opposite direction, and you will make false, irreproducible discoveries. It took us great effort to find the middle ground."

The first step in demonstrating CORE's ability to explain data with cores of widely varying lengths was to examine system performance with synthetic data. "Generating artificial data where the answers that the algorithm is supposed to provide are known in advance is standard practice in the field of data analysis," Krasnitz comments. "CORE passed this test to our satisfaction." The researchers then provided two demonstrations of CORE's utility with actual data. "By this we mean that on the one hand, CORE uncovers the genealogy of cells in a tumor, and on the other, can identify characteristic patterns of copy number variation in breast cancer."

Heat map of T10 incidence table, with rows corresponding to cells and columns to CORE cores. The amplification and deletion subtables are indicated by the blue and red horizontal bars. The order of cores in each subtable are left to right by the descending value of their CORE scores. Darker shades of gray correspond to higher values in the table. The order of the cells is clustered by the phylogenetic tree (on left) with horizontal distance related to distance. The tree yields a perfect separation of the (pseudo)diploid, hypodiploid, and the two aneuploid populations. The vertical color bar encodes the cell subpopulation label: yellow for aneuploid A, red for aneuploid B, blue for the hypodiploid, and green for diploid and pseudodiploid. Copyright © PNAS, doi:10.1073/pnas.1306909110

One of the ways the scientists applied their methodology was to determine tumor cell population phylogeny (that is, the history of cellular lineages as they change through time). To accomplish this, they applied CORE to a collection of DNA copy number profiles from single cells of a tumor. "These single-cell genomic profiles are now available thanks to highly innovative experimental work of Michael Wigler and his team," Krasnitz points out.

"Moreover," Krasnitz continues, "it's crucially important to determine how cells found in a tumor are related to each other genealogically." For example, he illustrates, they sometimes find that cells closely related to each other reside far apart in the tumor mass, which points to their ability to migrate rather than stay anchored as a healthy epithelial cell would. "However," he cautions, "our reconstructed genealogy must be accurate if we were to make such a claim."

Features that separate subpopulations are found among the cores. Indentifying such features also presented a statistical challenge. "What exactly do we mean," Krasnitz illustrates, "when we say that feature A is present in the genome of a cell? After all, cores are extended objects, not individual letters of the genome. What if a genomic region found in a cell is similar to the core but not quite identical to it? Can I still claim that the core is present in the genome of the cell? There are many different answers, depending on the degree of similarity required.' As a result, the researchers had to test many hypotheses, as is often the case when analyzing complex data. To alleviate the problem, he adds, they borrowed techniques from machine learning.

Another key difficulty was applying CORE to comparative genomic hybridization data from a large set of tumor samples in order to define regions of recurrent copy number aberration in breast cancer. "For each tumor analyzed," Krasnitz relates, "we had to deal with genomic data averaged over millions of tumor cells. This averaging weakens our ability to detect number changes present in some, but not all, the cells. We had to overcome this difficulty before we could apply CORE."

Heat map of T16 incidence table with rows corresponding to cells and columns to cores. The amplification and deletion subtables and row and column order and shading of the heat map are as described for Fig. 3. The vertical color bar encodes the cell subpopulation label: red for the primary aneuploid, black for the metastatic aneuploid, green for the primary diploid and pseudodiploid, and blue for the metastatic diploid. The phylogenetic tree for the cells is shown on the left, with a perfect separation of the primary and metastatic aneuploid populations. Copyright © PNAS, doi:10.1073/pnas.1306909110

In this case, Krasnitz explains, they were able to detect very narrow cores alongside cores that span entire chromosome arms. Moreover, some of the narrow cores point to known important breast cancer genes, such as ERBB2, a Herceptin® (trastuzumab) target. "This suggests that other narrow cores should be explored for presence of potential drug targets," notes Krasnitz.

Despite the array of significant challenges, Krasnitz sees one key innovative concept of the paper being of overriding value: explanation. "As one of our reviewers pointed out," he says, "maximizing explanation is an alternative to maximizing likelihood as an approach to modeling experimental data. We hope that this notion can take root in other areas, not only in genomics."

For example, Krasnitz points out that in network analysis, these ideas can be used to find an optimal set of hubs such that each node in the network is connected to at least one hub. Another example of a problem where CORE may be useful, he adds, is how to hire a fixed number of experts from a large pool of candidates so as to maximize their collective knowledge.

An interesting aspect of this research is that certain association measures are more favorable than others with regard to algorithmic complexity. "This is a good illustration of how interconnected modern science is," observes Krasnitz. "In the early stages of the work we were interested in one particular form of explanation. We first tried to solve the problem by brute force – that is, by exhaustively searching for the best among an enormous number of possible solutions – but soon ran out of computing power. We were surprised to find out that an analogous problem was very elegantly solved back in 1991 by two Israeli mathematicians working in the area of operations research, Refael Hassin and Arie Tamir, who demonstrated that in this particular case the complexity of the problem could be considerably reduced1. Unfortunately – as we later understood – this form of explanation, while advantageous computationally, is not ideal for problems in genomics."

Moving forward, the scientists are exploring other forms of explanation that are better able to take into account randomness present in genomic data. "First, we're planning to release CORE as public software to enable its use by other researchers," notes Krasnitz. "Internally, we have two ongoing studies – one using cores as markers for response to drug therapy in cancer, and another using CORE to improve molecular diagnostics of cancer patients."

Explore further: Abnormalities in new molecular pathway may increase breast cancer risk

More information: Target inference from collections of genomic intervals, PNAS June 18, 2013 vol. 110 no. 25 E2271-E2278, doi:10.1073/pnas.1306909110

Related
1Improved complexity bounds for location problems on the real line, Operations Research Letters 10 (1991), 395-402, doi:10.1016/0167-6377(91)90041-M (PDF)

Related Stories

Abnormalities in new molecular pathway may increase breast cancer risk

June 17, 2013
A new molecular pathway involving the gene ZNF365 has been identified and abnormalities in that pathway may predict worse outcomes for patients with breast cancer, according to data published in Cancer Discovery, a journal ...

Bacterial DNA may integrate into human genome more readily in tumor tissue

June 20, 2013
Bacterial DNA may integrate into the human genome more readily in tumors than in normal human tissue, according to a new study from the University of Maryland School of Medicine's Institute for Genome Sciences. Researchers ...

Scientists pair blood test and gene sequencing to detect cancer

November 28, 2012
Scientists at the Johns Hopkins Kimmel Cancer Center have combined the ability to detect cancer DNA in the blood with genome sequencing technology in a test that could be used to screen for cancers, monitor cancer patients ...

Mayo Clinic genomic analysis lends insight to prostate cancer

May 23, 2013
Mayo Clinic researchers have used next generation genomic analysis to determine that some of the more aggressive prostate cancer tumors have similar genetic origins, which may help in predicting cancer progression. The findings ...

Recommended for you

Scientists identify new way cells turn off genes

July 19, 2017
Cells have more than one trick up their sleeve for controlling certain genes that regulate fetal growth and development.

South Asian genomes could be boon for disease research, scientists say

July 18, 2017
The Indian subcontinent's massive population is nearing 1.5 billion according to recent accounts. But that population is far from monolithic; it's made up of nearly 5,000 well-defined sub-groups, making the region one of ...

Mutant yeast reveals details of the aberrant genomic machinery of children's high-grade gliomas

July 18, 2017
St. Jude Children's Research Hospital biologists have used engineered yeast cells to discover how a mutation that is frequently found in pediatric brain tumor high-grade glioma triggers a cascade of genomic malfunctions.

Late-breaking mutations may play an important role in autism

July 17, 2017
A study of nearly 6,000 families, combining three genetic sequencing technologies, finds that mutations that occur after conception play an important role in autism. A team led by investigators at Boston Children's Hospital ...

Newly identified genetic marker may help detect high-risk flu patients

July 17, 2017
Researchers have discovered an inherited genetic variation that may help identify patients at elevated risk for severe, potentially fatal influenza infections. The scientists have also linked the gene variant to a mechanism ...

Newly discovered gene variants link innate immunity and Alzheimer's disease

July 17, 2017
Three new gene variants, found in a genome wide association study of Alzheimer's disease (AD), point to the brain's immune cells in the onset of the disorder. These genes encode three proteins that are found in microglia, ...

0 comments

Please sign in to add a comment. Registration is free, and takes less than a minute. Read more

Click here to reset your password.
Sign in to get notified via email when new comments are made.