The degree of segmental aneuploidy measured by total copy number abnormalities predicts survival and recurrence in superficial gastroesophageal adenocarcinoma.

BACKGROUND: Prognostic biomarkers are needed for superficial gastroesophageal adenocarcinoma (EAC) to predict clinical outcomes and select therapy. Although recurrent mutations have been characterized in EAC, little is known about their clinical and prognostic significance. Aneuploidy is predictive of clinical outcome in many malignancies but has not been evaluated in superficial EAC. METHODS: We quantified copy number changes in 41 superficial EAC using Affymetrix SNP 6.0 arrays. We identified recurrent chromosomal gains and losses and calculated the total copy number abnormality (CNA) count for each tumor as a measure of aneuploidy. We correlated CNA count with overall survival and time to first recurrence in univariate and multivariate analyses. RESULTS: Recurrent segmental gains and losses involved multiple genes, including: HER2, EGFR, MET, CDK6, KRAS (recurrent gains); and FHIT, WWOX, CDKN2A/B, SMAD4, RUNX1 (recurrent losses). There was a 40-fold variation in CNA count across all cases. Tumors with the lowest and highest quartile CNA count had significantly better overall survival (p = 0.032) and time to first recurrence (p = 0.010) compared to those with intermediate CNA counts. These associations persisted when controlling for other prognostic variables. SIGNIFICANCE: SNP arrays facilitate the assessment of recurrent chromosomal gain and loss and allow high resolution, quantitative assessment of segmental aneuploidy (total CNA count). The non-monotonic association of segmental aneuploidy with survival has been described in other tumors. The degree of aneuploidy is a promising prognostic biomarker in a potentially curable form of EAC.

Serum protein profiles predict coronary artery disease in symptomatic patients referred for coronary angiography.

BACKGROUND: More than a million diagnostic cardiac catheterizations are performed annually in the US for evaluation of coronary artery anatomy and the presence of atherosclerosis. Nearly half of these patients have no significant coronary lesions or do not require mechanical or surgical revascularization. Consequently, the ability to rule out clinically significant coronary artery disease (CAD) using low cost, low risk tests of serum biomarkers in even a small percentage of patients with normal coronary arteries could be highly beneficial. METHODS: Serum from 359 symptomatic subjects referred for catheterization was interrogated for proteins involved in atherogenesis, atherosclerosis, and plaque vulnerability. Coronary angiography classified 150 patients without flow-limiting CAD who did not require percutaneous intervention (PCI) while 209 required coronary revascularization (stents, angioplasty, or coronary artery bypass graft surgery). Continuous variables were compared across the two patient groups for each analyte including calculation of false discovery rate (FDR ≤ 1%) and Q value (P value for statistical significance adjusted to ≤ 0.01). RESULTS: Significant differences were detected in circulating proteins from patients requiring revascularization including increased apolipoprotein B100 (APO-B100), C-reactive protein (CRP), fibrinogen, vascular cell adhesion molecule 1 (VCAM-1), myeloperoxidase (MPO), resistin, osteopontin, interleukin (IL)-1ß, IL-6, IL-10 and N-terminal fragment protein precursor brain natriuretic peptide (NT-pBNP) and decreased apolipoprotein A1 (APO-A1). Biomarker classification signatures comprising up to 5 analytes were identified using a tunable scoring function trained against 239 samples and validated with 120 additional samples. A total of 14 overlapping signatures classified patients without significant coronary disease (38% to 59% specificity) while maintaining 95% sensitivity for patients requiring revascularization. Osteopontin (14 times) and resistin (10 times) were most frequently represented among these diagnostic signatures. The most efficacious protein signature in validation studies comprised osteopontin (OPN), resistin, matrix metalloproteinase 7 (MMP7) and interferon γ (IFNγ) as a four-marker panel while the addition of either CRP or adiponectin (ACRP-30) yielded comparable results in five protein signatures. CONCLUSIONS: Proteins in the serum of CAD patients predominantly reflected (1) a positive acute phase, inflammatory response and (2) alterations in lipid metabolism, transport, peroxidation and accumulation. There were surprisingly few indicators of growth factor activation or extracellular matrix remodeling in the serum of CAD patients except for elevated OPN. These data suggest that many symptomatic patients without significant CAD could be identified by a targeted multiplex serum protein test without cardiac catheterization thereby eliminating exposure to ionizing radiation and decreasing the economic burden of angiographic testing for these patients.

MicroRNA-137 promoter methylation is associated with poorer overall survival in patients with squamous cell carcinoma of the head and neck.

BACKGROUND: The overall 5-year survival rate of approximately 60% for head and neck cancer patients has remained essentially unchanged over the past 30 years. MicroRNA-137 (miR-137) plays an essential role in cell-cycle control at the G1/S-phase checkpoint. However, the aberrant miR-137 promoter methylation observed in squamous cell carcinoma of the head and neck (SCCHN) suggests a tumor-specific molecular defect that may contribute to disease progression. METHODS: The goal of this study was to assess, in formalin-fixed, paraffin-embedded tumor tissue, the association between miR-137 promoter methylation and survival (both overall and disease free) and with prognostic factors including stage, tumor size, lymph node positivity, tumor grade, and surgical tumor margin positivity. RESULTS: The promoter methylation status of miR-137 was ascertained by methylation-specific polymerase chain reaction and detected in 11 of 67 SCCHN patients (16.4%), with no significant differences according to site (oral cavity, pharynx, larynx). Methylation of the miR-137 promoter was significantly associated with overall survival (hazard ratio, 3.68; 95% confidence interval, 1.01-13.38) but not with disease-free survival or any of the prognostic factors evaluated. CONCLUSIONS: The results of this study indicate that miR-137 is methylated in tumor tissue from pharyngeal and laryngeal squamous cancers, in addition to oral squamous cell carcinoma, and that miR-137 promoter methylation has potential utility as a prognostic marker for SCCHN.

A novel SNP analysis method to detect copy number alterations with an unbiased reference signal directly from tumor samples.

BACKGROUND: Genomic instability in cancer leads to abnormal genome copy number alterations (CNA) as a mechanism underlying tumorigenesis. Using microarrays and other technologies, tumor CNA are detected by comparing tumor sample CN to normal reference sample CN. While advances in microarray technology have improved detection of copy number alterations, the increase in the number of measured signals, noise from array probes, variations in signal-to-noise ratio across batches and disparity across laboratories leads to significant limitations for the accurate identification of CNA regions when comparing tumor and normal samples. METHODS: To address these limitations, we designed a novel "Virtual Normal" algorithm (VN), which allowed for construction of an unbiased reference signal directly from test samples within an experiment using any publicly available normal reference set as a baseline thus eliminating the need for an in-lab normal reference set. RESULTS: The algorithm was tested using an optimal, paired tumor/normal data set as well as previously uncharacterized pediatric malignant gliomas for which a normal reference set was not available. Using Affymetrix 250K Sty microarrays, we demonstrated improved signal-to-noise ratio and detected significant copy number alterations using the VN algorithm that were validated by independent PCR analysis of the target CNA regions. CONCLUSIONS: We developed and validated an algorithm to provide a virtual normal reference signal directly from tumor samples and minimize noise in the derivation of the raw CN signal. The algorithm reduces the variability of assays performed across different reagent and array batches, methods of sample preservation, multiple personnel, and among different laboratories. This approach may be valuable when matched normal samples are unavailable or the paired normal specimens have been subjected to variations in methods of preservation.

IDH1 mutations are common in malignant gliomas arising in adolescents: a report from the Children's Oncology Group.

PURPOSE: Recent studies have demonstrated a high frequency of IDH mutations in adult "secondary" malignant gliomas arising from preexisting lower grade lesions, often in young adults, but not in "primary" gliomas. Because pediatric malignant gliomas share some molecular features with adult secondary gliomas, we questioned whether a subset of these tumors also exhibited IDH mutations. EXPERIMENTAL DESIGN: We examined the frequency of IDH mutations, using real-time polymerase chain reaction and sequencing analysis, in a cohort of 43 pediatric primary malignant gliomas treated on the Children's Oncology Group ACNS0423 study. The relationship between IDH mutations and other molecular and clinical factors, and outcome, was evaluated. RESULTS: IDH1 mutations were observed in 7 of 43 (16.3%) tumors; no IDH2 mutations were observed. A striking age association was apparent in that mutations were noted in 7 of 20 tumors (35%) from children ≥14 years, but in 0 of 23 (0%) younger children (p = 0.0024). No association was observed with clinical factors other than age. One-year event-free survival was 86 ± 15% in the IDH-mutated group versus 64 ± 8% in the non-mutated group (p = 0.03, one-sided logrank test). One-year overall survival was 100% in patients with mutations versus 81 ± 6.7% in those without mutations (p = 0.035, one-sided logrank test). CONCLUSIONS: IDH1 mutations are common in malignant gliomas in older children, suggesting that a subset of these lesions may be biologically similar to malignant gliomas arising in younger adults and may be associated with a more favorable prognosis.

Reproducibility and performance of virtual karyotyping with SNP microarrays for the detection of chromosomal imbalances in formalin-fixed paraffin-embedded tissues.

BACKGROUND: Chromosomal imbalances are commonly seen in cancer and inherited genetic diseases. These imbalances may assist in the diagnosis, prognosis, and/or therapeutic management of certain neoplasms. Several methods for detecting chromosomal imbalances, such as, fluorescent in situ hybridization, array comparative genomic hybridization, and single nucleotide polymorphism (SNP) arrays have proven useful in formalin-fixed paraffin-embedded (FFPE) tissues. Here, we report the performance and reproducibility of virtual karyotyping of FFPE tissues with Affymetrix SNP arrays. METHODS: Virtual karyotypes from 442 FFPE tumor samples were generated using the Affymetrix GeneChip Mapping 10K Xba 2.0 and/or 250K Nsp SNP mapping arrays. Samples ranged from a few weeks to 17 years in archival storage. Virtual karyotypes were assessed for copy number changes, loss of heterozygosity, and acquired uniparental disomy. RESULTS: Overall, 75.3% of samples produced interpretable virtual karyotypes with the 10K arrays and 76.7% in the 250K arrays. Parameters for the selection of samples for hybridization were determined, which increased the success rate in both platforms to 81.3 and 92.6%, respectively. FFPE virtual karyotypes generated with both 10K Xba 2.0 and 250K Nsp arrays showed 100% concordance in intralaboratory and interlaboratory reproducibility studies. Samples older than 7 years showed decreased performance. CONCLUSIONS: SNP arrays are a reliable, reproducible, and robust platform for the virtual karyotyping of FFPE tumor tissues with performance characteristics adequate for clinical application. Parameters that most significantly affected sample performance were sample age and storage conditions.

Mismatch repair deficiency is an uncommon mechanism of alkylator resistance in pediatric malignant gliomas: a report from the Children's Oncology Group.

BACKGROUND: Alkylating agents are commonly used in the treatment of childhood malignant gliomas. Overexpression of O(6)-methylguanine-DNA methyltransferase (MGMT) constitutes an important mechanism for resistance to such agents, and MGMT status has been associated with outcome in several recent trials. Deficiency in mismatch repair (MMR) function has been implicated in preclinical studies as an additional potential mechanism of resistance to methylating agents, such as temozolomide, independent of tumor MGMT status. However, the frequency of this abnormality as a clinical resistance mechanism in childhood malignant gliomas has not been well characterized. METHODS: To address this issue, we examined the frequency of microsatellite instability (MSI), a marker of defective MMR, in a series of 68 tumors, derived from newly diagnosed patients treated on the Children's Cancer Group 945 study, and the Children's Oncology Group ACNS0126 and 0423 studies. MSI was assessed using a panel of six microsatellite markers, including BAT-25, BAT-26, CAT-25, D2S123, D5S346, and D17S250. MGMT immunoreactivity was assessed in parallel to allow comparison of the relative incidence of MGMT overexpression and MSI. RESULTS: Only three tumors had high-level MSI involving three or more markers; the remainder had no MSI at any of the loci examined. These children did not have unusual features in terms of their outcome. In contrast to the infrequency of MSI, 25 tumors (37%) exhibited MGMT overexpression as assessed by immunohistochemistry. None of the tumors with MSI exhibited overexpression of MGMT. CONCLUSION: MMR deficiency is an infrequent contributor to initial alkylator resistance in children with malignant gliomas.

Identification of tissue of origin in carcinoma of unknown primary with a microarray-based gene expression test.

BACKGROUND: Carcinomas of unknown primary (CUP) represent approximately 3%-5% of malignant neoplasms. Identifying the tissue of origin (TOO) in these tumors allows for more specific treatment and improves outcomes. However, primary classification remains a challenge in many cases. We evaluated the ability of a microarray-based gene expression test to identify the TOO in tumor specimens from 21 patients with a diagnosis of CUP. METHODS: The Pathwork TOO Test was used to measure gene expression patterns for 1550 genes; these were compared for similarity to patterns from 15 known tissue types. RESULTS: The TOO Test yielded a clear single positive call for the primary site in 16 of 21 (76%) specimens and was indeterminate in 5 (24%). The positive results were consistent with clinicopathologic suggestions in 10 of the 16 cases (62%). In the remaining six cases the positive results were considered plausible based on clinical information. Positive calls included colorectal (5), breast (4), ovarian (3), lung (2), and pancreas (2). The TOO Test ruled out an average of 11 primary tissues in each CUP specimen. CONCLUSION: The Pathwork TOO Test reduced diagnostic uncertainty in all CUP cases and could be a valuable addition or alternative to current diagnostic methods for classifying uncertain primary cancers.

Multicenter validation of a 1,550-gene expression profile for identification of tumor tissue of origin.

PURPOSE: Malignancies found in unexpected locations or with poorly differentiated morphologies can pose a significant challenge for tissue of origin determination. Current histologic and imaging techniques fail to yield definitive identification of the tissue of origin in a significant number of cases. The aim of this study was to validate a predefined 1,550-gene expression profile for this purpose. METHODS: Four institutions processed 547 frozen specimens representing 15 tissues of origin using oligonucleotide microarrays. Half of the specimens were metastatic tumors, with the remainder being poorly differentiated and undifferentiated primary cancers chosen to resemble those that present as a clinical challenge. RESULTS: In this blinded multicenter validation study the 1,550-gene expression profile was highly informative in tissue determination. The study found overall sensitivity (positive percent agreement with reference diagnosis) of 87.8% (95% CI, 84.7% to 90.4%) and overall specificity (negative percent agreement with reference diagnosis) of 99.4% (95% CI, 98.3% to 99.9%). Performance within the subgroup of metastatic tumors (n = 258) was found to be slightly lower than that of the poorly differentiated and undifferentiated primary tumor subgroup, 84.5% and 90.7%, respectively (P = .04). Differences between individual laboratories were not statistically significant. CONCLUSION: This study represents the first adequately sized, multicenter validation of a gene-expression profile for tissue of origin determination restricted to poorly differentiated and undifferentiated primary cancers and metastatic tumors. These results indicate that this profile should be a valuable addition or alternative to currently available diagnostic methods for the evaluation of uncertain primary cancers.

Virtual karyotyping with SNP microarrays reduces uncertainty in the diagnosis of renal epithelial tumors.

BACKGROUND: Renal epithelial tumors are morphologically, biologically, and clinically heterogeneous. Different morphologic subtypes require specific management due to markedly different prognosis and response to therapy. Each common subtype has characteristic chromosomal gains and losses, including some with prognostic value. However, copy number information has not been readily accessible for clinical purposes and thus has not been routinely used in the diagnostic evaluation of these tumors. This information can be useful for classification of tumors with complex or challenging morphology. 'Virtual karyotypes' generated using SNP arrays can readily detect characteristic chromosomal lesions in paraffin embedded renal tumors and can be used to correctly categorize the common subtypes with performance characteristics that are amenable for routine clinical use. METHODS: To investigate the use of virtual karyotypes for diagnostically challenging renal epithelial tumors, we evaluated 25 archived renal neoplasms where sub-classification could not be definitively rendered based on morphology and other ancillary studies. We generated virtual karyotypes with the Affymetrix 10 K 2.0 mapping array platform and identified the presence of genomic lesions across all 22 autosomes. RESULTS: In 91% of challenging cases the virtual karyotype unambiguously detected the presence or absence of chromosomal aberrations characteristic of one of the common subtypes of renal epithelial tumors, while immunohistochemistry and fluorescent in situ hybridization had no or limited utility in the diagnosis of these tumors. CONCLUSION: These results show that virtual karyotypes generated by SNP arrays can be used as a practical ancillary study for the classification of renal epithelial tumors with complex or ambiguous morphology.

Whole genome SNP arrays as a potential diagnostic tool for the detection of characteristic chromosomal aberrations in renal epithelial tumors.

Renal tumors with complex or unusual morphology require extensive workup for accurate classification. Chromosomal aberrations that define subtypes of renal epithelial neoplasms have been reported. We explored if whole-genome chromosome copy number and loss-of-heterozygosity analysis with single nucleotide polymorphism (SNP) arrays can be used to identify these aberrations and classify renal epithelial tumors. We analyzed 20 paraffin-embedded tissues representing clear cell, papillary renal and chromophobe renal cell carcinoma, as well as oncocytoma with Affymetrix GeneChip 10K 2.0 Mapping arrays. SNP array results were in concordance with known genetic aberrations for each renal tumor subtype. Additional chromosomal aberrations were detected in all renal cell tumor types. The unique patterns allowed 19 out of 20 tumors to be readily categorized by their chromosomal copy number aberrations. One papillary renal cell carcinoma type 2 did not show the characteristic 7/17 trisomies. Clustering using the median copy number of each chromosomal arm correlated with histological class when using a restricted set of chromosomes. In addition, three morphologically challenging tumors were analyzed to explore the potential clinical utility of this method. In these cases, the SNP array-based copy number evaluation yielded information with potential clinical value. These results show that SNP arrays can detect characteristic chromosomal aberrations in paraffin-embedded renal tumors, and thus offer a high-resolution, genome-wide method that can be used as an ancillary study for classification and potentially for prognostic stratification of these tumors.

Optimization of the Affymetrix GeneChip Mapping 10K 2.0 Assay for routine clinical use on formalin-fixed paraffin-embedded tissues.

The use of chromosomal copy number changes as markers for tumor behavior or as prognostic markers for patient outcome has been suggested. However, current clinically used technologies cannot perform genome-wide assessment of chromosome copy number and analysis of loss of heterozygosity in the same assay for paraffin-embedded tissue. We have optimized the Affymetrix GeneChip Mapping Assay for the 10K 2.0 array for use with formalin-fixed, paraffin-embedded (FFPE) tissues. This technology uses single nucleotide polymorphism (SNP) arrays to assess the changes in chromosomal copy number and loss of heterozygosity. DNA from 3 paired tumor/normal samples of adrenal tumors and 4 samples of renal tumors were processed with modifications to the manufacturer's protocol. Modifications at different steps were evaluated for their effects on SNP signal-detection and call rates. Frozen samples showed 99.6%+/-0.3% signal-detection rates and 94.7%+/-3.0% SNP call rates. FFPE samples labeled with the original protocol failed to produce enough polymerase chain reaction products for hybridization, whereas the same samples processed with the optimized protocol had signal-detection rates of 97.4%+/-0.018% and SNP call rates of 90.9%+/-0.034%. The average SNP call concordance between fresh and matching FFPE samples was 96%. Chromosomal aberrations detected in the frozen tumors were also detected in the FFPE tissues. Our optimized protocol significantly improves the performance of the FFPE samples in the Affymetrix GeneMapping Assay with the 10K 2.0 SNP array. This optimized protocol opens up the potential for the GeneChip Mapping assay to be used in the development of clinical test assays.

Interlaboratory performance of a microarray-based gene expression test to determine tissue of origin in poorly differentiated and undifferentiated cancers.

Clinical workup of metastatic malignancies of unknown origin is often arduous and expensive and is reported to be unsuccessful in 30 to 60% of cases. Accurate classification of uncertain primary cancers may improve with microarray-based gene expression testing. We evaluated the analytical performance characteristics of the Pathwork tissue of origin test, which uses expression signals from 1668 probe sets in a gene expression microarray, to quantify the similarity of tumor specimens to 15 known tissues of origin. Sixty archived tissue specimens from poorly and undifferentiated tumors (metastatic and primary) were analyzed at four laboratories representing a wide range of preanalytical conditions (eg, personnel, reagents, instrumentation, and protocols). Cross-laboratory comparisons showed highly reproducible results between laboratories, with correlation coefficients between 0.95 to 0.97 for measurements of similarity scores, and an average 93.8% overall concordance between laboratories in terms of final tissue calls. Bland-Altman plots (mean coefficients of reproducibility of 32.48+/-3.97) and kappa statistics (kappa >0.86) also indicated a high level of agreement between laboratories. We conclude that the Pathwork tissue of origin test is a robust assay that produces consistent results in diverse laboratory conditions reflecting the preanalytical variations found in the everyday clinical practice of molecular diagnostics laboratories.

Gene expression profiles of prostate cancer reveal involvement of multiple molecular pathways in the metastatic process.

BACKGROUND: Prostate cancer is characterized by heterogeneity in the clinical course that often does not correlate with morphologic features of the tumor. Metastasis reflects the most adverse outcome of prostate cancer, and to date there are no reliable morphologic features or serum biomarkers that can reliably predict which patients are at higher risk of developing metastatic disease. Understanding the differences in the biology of metastatic and organ confined primary tumors is essential for developing new prognostic markers and therapeutic targets. METHODS: Using Affymetrix oligonucleotide arrays, we analyzed gene expression profiles of 24 androgen-ablation resistant metastatic samples obtained from 4 patients and a previously published dataset of 64 primary prostate tumor samples. Differential gene expression was analyzed after removing potentially uninformative stromal genes, addressing the differences in cellular content between primary and metastatic tumors. RESULTS: The metastatic samples are highly heterogenous in expression; however, differential expression analysis shows that 415 genes are upregulated and 364 genes are downregulated at least 2 fold in every patient with metastasis. The expression profile of metastatic samples reveals changes in expression of a unique set of genes representing both the androgen ablation related pathways and other metastasis related gene networks such as cell adhesion, bone remodelling and cell cycle. The differentially expressed genes include metabolic enzymes, transcription factors such as Forkhead Box M1 (FoxM1) and cell adhesion molecules such as Osteopontin (SPP1). CONCLUSION: We hypothesize that these genes have a role in the biology of metastatic disease and that they represent potential therapeutic targets for prostate cancer.

In vitro transcription amplification and labeling methods contribute to the variability of gene expression profiling with DNA microarrays.

The effect of different amplification and labeling methods on DNA microarray expression results has not been previously delineated. To analyze the variation associated with widely accepted T7-based RNA amplificationand labeling methods, aliquots of the Stratagene Human Universal Reference RNA were labeled using three eukaryotic target preparation methods followed by uniform replicate array hybridization (Affymetrix U95Av2). Method-dependent variability was observed in the yield and size distribution of labeled products, as well as in the gene expression results. A significant increase in short transcripts, when compared to unamplified mRNA, was observed in methods with long in vitro transcription reactions. Intramethod reproducibility showed correlation coefficients >0.99, whereas intermethod comparisons showed coefficients ranging from 0.94 to 0.98 and a nearly twofold increase in coefficient of variation. Fold amplification for each method positively correlated with the number of genes present. Our experiments uncovered two factors that introduced significant bias in gene expression data: the number of labeled nucleotides, which introduces sequence-dependent bias, and the length of the in vitro transcription reaction, which introduces transcript size-dependent bias. This study provides evidence that variability in expression data may be caused, in part, by differences in amplification and labeling protocols.

The colonization of land by animals: molecular phylogeny and divergence times among arthropods.

BACKGROUND: The earliest fossil evidence of terrestrial animal activity is from the Ordovician, approximately 450 million years ago (Ma). However, there are earlier animal fossils, and most molecular clocks suggest a deep origin of animal phyla in the Precambrian, leaving open the possibility that animals colonized land much earlier than the Ordovician. To further investigate the time of colonization of land by animals, we sequenced two nuclear genes, glyceraldehyde-3-phosphate dehydrogenase and enolase, in representative arthropods and conducted phylogenetic and molecular clock analyses of those and other available DNA and protein sequence data. To assess the robustness of animal molecular clocks, we estimated the deuterostome-arthropod divergence using the arthropod fossil record for calibration and tunicate instead of vertebrate sequences to represent Deuterostomia. Nine nuclear and 15 mitochondrial genes were used in phylogenetic analyses and 61 genes were used in molecular clock analyses. RESULTS: Significant support was found for the unconventional pairing of myriapods (millipedes and centipedes) with chelicerates (spiders, scorpions, horseshoe crabs, etc.) using nuclear and mitochondrial genes. Our estimated time for the divergence of millipedes (Diplopoda) and centipedes (Chilopoda) was 442 +/- 50 Ma, and the divergence of insects and crustaceans was estimated as 666 +/- 58 Ma. Our results also agree with previous studies suggesting a deep divergence (approximately 1100 - 900 Ma) for arthropods and deuterostomes, considerably predating the Cambrian Explosion seen in the animal fossil record. CONCLUSIONS: The consistent support for a close relationship between myriapods and chelicerates, using mitochondrial and nuclear genes and different methods of analysis, suggests that this unexpected result is not an artefact of analysis. We propose the name Myriochelata for this group of animals, which includes many that immobilize prey with venom. Our molecular clock analyses using arthropod fossil calibrations support earlier studies using vertebrate calibrations in finding that deuterostomes and arthropods diverged hundreds of millions of years before the Cambrian explosion. However, our molecular time estimate for the divergence of millipedes and centipedes is close to the divergence time inferred from fossils. This suggests that arthropods may have adapted to the terrestrial environment relatively late in their evolutionary history.