Differences in gastric carcinoma microenvironment stratify according to EBV infection intensity: implications for possible immune adjuvant therapy.

Epstein-Barr virus (EBV) is associated with roughly 10% of gastric carcinomas worldwide (EBVaGC). Although previous investigations provide a strong link between EBV and gastric carcinomas, these studies were performed using selected EBV gene probes. Using a cohort of gastric carcinoma RNA-seq data sets from The Cancer Genome Atlas (TCGA), we performed a quantitative and global assessment of EBV gene expression in gastric carcinomas and assessed EBV associated cellular pathway alterations. EBV transcripts were detected in 17% of samples but these samples varied significantly in EBV coverage depth. In four samples with the highest EBV coverage (hiEBVaGC - high EBV associated gastric carcinoma), transcripts from the BamHI A region comprised the majority of EBV reads. Expression of LMP2, and to a lesser extent, LMP1 were also observed as was evidence of abortive lytic replication. Analysis of cellular gene expression indicated significant immune cell infiltration and a predominant IFNG response in samples expressing high levels of EBV transcripts relative to samples expressing low or no EBV transcripts. Despite the apparent immune cell infiltration, high levels of the cytotoxic T-cell (CTL) and natural killer (NK) cell inhibitor, IDO1, was observed in the hiEBVaGCs samples suggesting an active tolerance inducing pathway in this subgroup. These results were confirmed in a separate cohort of 21 Vietnamese gastric carcinoma samples using qRT-PCR and on tissue samples using in situ hybridization and immunohistochemistry. Lastly, a panel of tumor suppressors and candidate oncogenes were expressed at lower levels in hiEBVaGC versus EBV-low and EBV-negative gastric cancers suggesting the direct regulation of tumor pathways by EBV.

Modeling, simulation and analysis of methylation profiles from reduced representation bisulfite sequencing experiments.

The ENCODE project has funded the generation of a diverse collection of methylation profiles using reduced representation bisulfite sequencing (RRBS) technology, enabling the analysis of epigenetic variation on a genomic scale at single-site resolution. A standard application of RRBS experiments is in the location of differentially methylated regions (DMRs) between two sets of samples. Despite numerous publications reporting DMRs identified from RRBS datasets, there have been no formal analyses of the effects of experimental and biological factors on the performance of existing or newly developed analytical methods. These factors include variable read coverage, differing group sample sizes across genomic regions, uneven spacing between CpG dinucleotide sites, and correlation in methylation levels among sites in close proximity. To better understand the interplay among technical and biological variables in the analysis of RRBS methylation profiles, we have developed an algorithm for the generation of experimentally realistic RRBS datasets. Applying insights derived from our simulation studies, we present a novel procedure that can identify DMRs spanning as few as three CpG sites with both high sensitivity and specificity. Using RRBS data from muscle vs. non-muscle cell cultures as an example, we demonstrate that our method reveals many more DMRs that are likely to be of biological significance than previous methods.

TgMAPK1 is a Toxoplasma gondii MAP kinase that hijacks host MKK3 signals to regulate virulence and interferon-γ-mediated nitric oxide production.

The parasite Toxoplasma gondii controls tissue-specific nitric oxide (NO), thereby augmenting virulence and immunopathology through poorly-understood mechanisms. We now identify TgMAPK1, a Toxoplasma mitogen-activated protein kinase (MAPK), as a virulence factor regulating tissue-specific parasite burden by manipulating host interferon (IFN)-γ-mediated inducible nitric oxide synthase (iNOS). Toxoplasma with reduced TgMAPK1 expression (TgMAPK1(lo)) demonstrated that TgMAPK1 facilitates IFN-γ-driven p38 MAPK activation, reducing IFN-γ-generated NO in an MKK3-dependent manner, blunting IFN-γ-mediated parasite control. TgMAPK1(lo) infection in wild type mice produced ≥ten-fold lower parasite burden versus control parasites with normal TgMAPK1 expression (TgMAPK1(con)). Reduced parasite burdens persisted in IFN-γ KO mice, but equalized in normally iNOS-replete organs from iNOS KO mice. Parasite MAPKs are far less studied than other parasite kinases, but deserve additional attention as targets for immunotherapy and drug discovery.

Tumor necrosis factor receptor 1 functions as a tumor suppressor.

Tumor necrosis factor (TNF) is a key player in inflammatory bowel disease and has been variably associated with carcinogenesis, but details of the cross talk between inflammatory and tumorigenic pathways remain incompletely understood. It has been shown that, in C57BL/6 mice, signaling via TNF receptor 1 (TNFR1) is protective from injury and inflammation in experimental colitis. Therefore, we hypothesized that loss of TNFR1 signaling would confer increased risk of developing colitis-associated carcinoma. Using three models of murine tumorigenesis based on repeated bouts of inflammation or systemic tumor initiator, we sought to determine the roles of TNF and TNFR1 with regard to neoplastic transformation in the colon in wild-type (WT), TNFR1 knockout (R1KO), and TNF knockout (TNFKO) mice. We found R1KO animals to have more severe disease, as defined by weight loss, hematochezia, and histology. TNFKO mice demonstrated less weight loss but were consistently smaller, and rates and duration of hematochezia were comparable to WT mice. Histological inflammation scores were higher and neoplastic lesions occurred more frequently and earlier in R1KO mice. Apoptosis is not affected in R1KO mice although epithelial proliferation following injury is more ardent even before tumorigenesis is apparent. Lastly, there is earlier and more intense expression of activated ß-catenin in these mice, implying a connection between TNFR1 and Wnt signaling. Taken together, these findings show that in the context of colitis-associated carcinogenesis TNFR1 functions as a tumor suppressor, exerting this effect not via apoptosis but by modulating activation of ß-catenin and controlling epithelial proliferation.

In vitro high-capacity assay to quantify the clonal heterogeneity in trilineage potential of mesenchymal stem cells reveals a complex hierarchy of lineage commitment.

In regenerative medicine, bone marrow is a promising source of mesenchymal stem cells (MSCs) for a broad range of cellular therapies. This research addresses a basic prerequisite to realize the therapeutic potential of MSCs by developing a novel high-capacity assay to quantify the clonal heterogeneity in potency that is inherent to MSC preparations. The assay utilizes a 96-well format to (1) classify MSCs according to colony-forming efficiency as a measure of proliferation capacity and trilineage potential to exhibit adipo-, chondro-, and osteogenesis as a measure of multipotency and (2) preserve a frozen template of MSC clones of known potency for future use. The heterogeneity in trilineage potential of normal bone marrow MSCs is more complex than previously reported: all eight possible categories of trilineage potential were detected. In this study, the average colony-forming efficiency of MSC preparations was 55-62%, and tripotent MSCs accounted for nearly 50% of the colony-forming cells. The multiple phenotypes detected in this study infer a more convoluted hierarchy of lineage commitment than described in the literature. Greater cell amplification, colony-forming efficiency, and colony diameter for tri- versus unipotent clones suggest that MSC proliferation may be a function of potency. CD146 may be a marker of multipotency, with approximately 2-fold difference in mean fluorescence intensity between tri- and unipotent clones. The significance of these findings is discussed in the context of the efficacy of MSC therapies. The in vitro assay described herein will likely have numerous applications given the importance of heterogeneity to the therapeutic potential of MSCs.

Clonal analysis of the proliferation potential of human bone marrow mesenchymal stem cells as a function of potency.

Human mesenchymal stem cells (MSCs) from bone marrow are a heterogeneous ensemble of progenitors and lineage-committed cells, with a broad range of regenerative properties. Ex vivo expansion to produce sufficient quantities of MSCs is essential for most therapeutic applications. The present study resolves the relationship between proliferation potential of MSCs and their potency. Clonal analysis generated single-cell derived colonies of MSCs that were classified according to their trilineage potential to exhibit adipo- (A), chondro- (C), and osteogenesis (O) as a measure of potency. Multipotent OAC clones were highly proliferative with colony-forming efficiencies that ranged from 35% to 90%; whereas, O clones formed colonies with an efficiency of 5% or less (P < 0.01). Similar trends were evident during ex vivo expansion: for example, the median specific growth rate was 0.8 day(-1) (20 h doubling time) for cultures inoculated with OAC clones and was 5-fold less for inocula of O clones (P < 0.01). OA and OC clones had similar proliferation potentials. More than 75% of cells in subconfluent cultures inoculated with O clones stained positive for senescence-associated ß-galactosidase activity vs. less than 10% for OAC clones (P < 0.001). Apoptotic cells were in the minority for all potency groups. Preliminary data generated during clonal analysis suggest that osteogenic potential of MSCs to produce mineralized matrix is a function of potency, as well. These results are discussed in the context of the preparation of efficacious MSC therapies by ex vivo expansion.

Genetic requirements for the survival of tubercle bacilli in primates.

BACKGROUND: Tuberculosis (TB) leads to the death of 1.7 million people annually. The failure of the bacille Calmette-Guérin vaccine, synergy between AIDS and TB, and the emergence of drug resistance have worsened this situation. It is imperative to delineate the mechanisms employed by Mycobacterium tuberculosis to successfully infect and persist in mammalian lungs. METHODS: Nonhuman primates (NHPs) are arguably the best animal system to model critical aspects of human TB. We studied genes essential for growth and survival of M. tuberculosis in the lungs of NHPs experimentally exposed to aerosols of an M. tuberculosis transposon mutant library. RESULTS: Mutants in 108 M. tuberculosis genes (33.13% of all genes tested) were attenuated for in vivo growth. Comparable studies have reported the attenuation of only approximately 6% of mutants in mice. The M. tuberculosis mutants attenuated for in vivo survival in primates were involved in the transport of various biomolecules, including lipid virulence factors; biosynthesis of cell-wall arabinan and peptidoglycan; DNA repair; sterol metabolism; and mammalian cell entry. CONCLUSIONS: Our study highlights the various virulence mechanisms employed by M. tuberculosis to overcome the hostile environment encountered during infection of primates. Prophylactic approaches aimed against bacterial factors that respond to such in vivo stressors have the potential to prevent infection at an early stage, thus likely reducing the extent of transmission of M. tuberculosis.

Adult human mesenchymal stem cells enhance breast tumorigenesis and promote hormone independence.

Adult human mesenchymal stem cells (hMSCs) have been shown to home to sites of breast cancer and integrate into the tumor stroma. We demonstrate here the effect of hMSCs on primary breast tumor growth and the progression of these tumors to hormone independence. Co-injection of bone marrow-derived hMSCs enhances primary tumor growth of the estrogen receptor-positive, hormone-dependent breast carcinoma cell line MCF-7 in the presence or absence of estrogen in SCID/beige mice. We also show hormone-independent growth of MCF-7 cells when co-injected with hMSCs. These effects were found in conjunction with increased immunohistochemical staining of the progesterone receptor in the MCF-7/hMSC tumors as compared to MCF-7 control tumors. This increase in PgR expression indicates a link between MCF-7 cells and MSCs through ER-mediated signaling. Taken together, our data reveal the relationship between tumor microenvironment and tumor growth and the progression to hormone independence. This tumor stroma-cell interaction may provide a novel target for the treatment of estrogen receptor-positive, hormone-independent, and endocrine-resistant breast carcinoma.

Modeling dependence in methylation patterns with application to ovarian carcinomas.

Changes in cytosine methylation at CpG nucleotides are observed in many cancers and offer great potential for translational research. Diseases such as ovarian cancer that are especially challenging to diagnose and treat are of particular interest, and abnormal methylation in the tandem repeats Sat2 and NBL2 has been observed in a collection of ovarian carcinomas. In earlier analyses of double-stranded methylation patterns in 0.2 kb regions of Sat2 and NBL2, we detected clusters of identically methylated sites in close proximity. These clusters could not be explained by random variation, and our findings suggested a high degree of site-to-site dependence. However, previously developed stochastic models for methylation change have either treated CpG sites independently or employed a context dependent approach to adjust model parameters according to regional methylation levels. In this paper, we introduce a novel neighboring sites model as an alternative methodology for considering dependence in methylation patterns, and we compare the three models in their ability to generate simulated sequences statistically similar to our Sat2 and NBL2 carcinoma samples.

Modeling Methylation Patterns with Long Read Sequencing Data.

Variation in cytosine methylation at CpG dinucleotides is often observed in genomic regions, and analysis typically focuses on estimating the proportion of methylated sites observed in a given region and comparing these levels across samples to determine association with conditions of interest. While sites are tacitly treated as independent, when observed at the level of individual molecules methylation patterns exhibit strong evidence of local spatial dependence. We previously developed a neighboring sites model to account for correlation and clustering behavior observed in two tandem repeat regions in a collection of ovarian carcinomas. We now introduce extensions of the model that account for the effect of distance between sites as well as asymmetric correlation in de novo methylation and demethylation rates. We apply our models to published data from a whole genome bisulfite sequencing experiment using long reads, estimating model parameters for a selection of CpG-dense regions spanning between 21 and 67 sites. Our methods detect evidence of local spatial correlation as a function of site-to-site distance and demonstrate the added value of employing long read sequencing data in epigenetic research.

Serum cadmium levels in pancreatic cancer patients from the East Nile Delta region of Egypt.

UNLABELLED: The northeast Nile Delta region exhibits a high incidence of early-onset pancreatic cancer. It is well documented that this region has one of the highest levels of pollution in Egypt. Epidemiologic studies have suggested that cadmium, a prevalent pollutant in the northeast Nile Delta region, plays a role in the development of pancreatic cancer. OBJECTIVE: We aimed to assess serum cadmium levels as markers of exposure in pancreatic cancer patients and noncancer comparison subjects from the same region in Egypt. DESIGN AND PARTICIPANTS: We assessed serum cadmium levels of 31 newly diagnosed pancreatic cancer patients and 52 hospital comparison subjects from Mansoura, Egypt. EVALUATION/MEASUREMENTS: Serum cadmium levels were measured using a novel immunoassay procedure. RESULTS: We found a significant difference between the mean serum cadmium levels in patients versus comparison subjects (mean+/-SD, 11.1+/-7.7 ng/mL vs. 7.1+/-5.0 ng/mL, respectively; p=0.012) but not in age, sex, residence, occupation, or smoking status. The odds ratio (OR) for pancreatic cancer risk was significant for serum cadmium level [OR=1.12; 95% confidence interval (CI), 1.04-1.23; p=0.0089] and farming (OR=3.25; 95% CI, 1.03-11.64; p=0.0475) but not for age, sex, residence, or smoking status. CONCLUSIONS: The results from this pilot study suggest that pancreatic cancer in the East Nile Delta region is significantly associated with high levels of serum cadmium and farming. RELEVANCE TO CLINICAL PRACTICE/PUBLIC HEALTH: Future studies should further investigate the etiologic relationship between cadmium exposure and pancreatic carcinogenesis in cadmium-exposed populations.

A signal-to-noise analysis of phylogeny estimation by neighbor-joining: Insufficiency of polynomial length sequences.

Phylogeny reconstruction is the process of inferring evolutionary relationships from molecular sequences, and methods that are expected to accurately reconstruct trees from sequences of reasonable length are highly desirable. To formalize this concept, the property of fast-convergence has been introduced to describe phylogeny reconstruction methods that, with high probability, recover the true tree from sequences that grow polynomially in the number of taxa n. While provably fast-converging methods have been developed, the neighbor-joining (NJ) algorithm of Saitou and Nei remains one of the most popular methods used in practice. This algorithm is known to converge for sequences that are exponential in n, but no lower bound for its convergence rate has been established. To address this theoretical question, we analyze the performance of the NJ algorithm on a type of phylogeny known as a 'caterpillar tree'. We find that, for sequences of polynomial length in the number of taxa n, the variability of the NJ criterion is sufficiently high that the algorithm is likely to fail even in the first step of the phylogeny reconstruction process, regardless of the degree of polynomial considered. This result demonstrates that, for general n-taxa trees, the exponential bound cannot be improved.

RNA CoMPASS: a dual approach for pathogen and host transcriptome analysis of RNA-seq datasets.

High-throughput RNA sequencing (RNA-seq) has become an instrumental assay for the analysis of multiple aspects of an organism's transcriptome. Further, the analysis of a biological specimen's associated microbiome can also be performed using RNA-seq data and this application is gaining interest in the scientific community. There are many existing bioinformatics tools designed for analysis and visualization of transcriptome data. Despite the availability of an array of next generation sequencing (NGS) analysis tools, the analysis of RNA-seq data sets poses a challenge for many biomedical researchers who are not familiar with command-line tools. Here we present RNA CoMPASS, a comprehensive RNA-seq analysis pipeline for the simultaneous analysis of transcriptomes and metatranscriptomes from diverse biological specimens. RNA CoMPASS leverages existing tools and parallel computing technology to facilitate the analysis of even very large datasets. RNA CoMPASS has a web-based graphical user interface with intrinsic queuing to control a distributed computational pipeline. RNA CoMPASS was evaluated by analyzing RNA-seq data sets from 45 B-cell samples. Twenty-two of these samples were derived from lymphoblastoid cell lines (LCLs) generated by the infection of naïve B-cells with the Epstein Barr virus (EBV), while another 23 samples were derived from Burkitt's lymphomas (BL), some of which arose in part through infection with EBV. Appropriately, RNA CoMPASS identified EBV in all LCLs and in a fraction of the BLs. Cluster analysis of the human transcriptome component of the RNA CoMPASS output clearly separated the BLs (which have a germinal center-like phenotype) from the LCLs (which have a blast-like phenotype) with evidence of activated MYC signaling and lower interferon and NF-kB signaling in the BLs. Together, this analysis illustrates the utility of RNA CoMPASS in the simultaneous analysis of transcriptome and metatranscriptome data. RNA CoMPASS is freely available at http://rnacompass.sourceforge.net/.

Identification of biomarkers for tuberculosis susceptibility via integrated analysis of gene expression and longitudinal clinical data.

Tuberculosis (TB) is an infectious disease caused by the bacteria Mycobacterium tuberculosis (Mtb) that affects millions of people worldwide. The majority of individuals who are exposed to Mtb develop latent infections, in which an immunological response to Mtb antigens is present but there is no clinical evidence of disease. Because currently available tests cannot differentiate latent individuals who are at low risk from those who are highly susceptible to developing active disease, there is considerable interest in the identification of diagnostic biomarkers that can predict reactivation of latent TB. We present results from our analysis of a controlled longitudinal experiment in which a group of rhesus macaques were exposed to a low dose of Mtb to study their progression to latent infection or active disease. Subsets of the animals were then euthanized at scheduled time points, and granulomas taken from their lungs were assayed for gene expression using microarrays. The clinical profiles associated with the animals following Mtb exposure revealed considerable variability, and we developed models for the disease trajectory for each subject using a Bayesian hierarchical B-spline approach. Disease severity estimates were derived from these fitted curves and included as covariates in linear models to identify genes significantly associated with disease progression. Our results demonstrate that the incorporation of clinical data increases the value of information extracted from the expression profiles and contributes to the identification of predictive biomarkers for TB susceptibility.

Cell-surface expression of neuron-glial antigen 2 (NG2) and melanoma cell adhesion molecule (CD146) in heterogeneous cultures of marrow-derived mesenchymal stem cells.

Cellular heterogeneity of mesenchymal stem cells (MSCs) impedes their use in regenerative medicine. The objective of this research is to identify potential biomarkers for the enrichment of progenitors from heterogeneous MSC cultures. To this end, the present study examines variation in expression of neuron-glial antigen 2 (NG2) and melanoma cell adhesion molecule (CD146) on the surface of MSCs derived from human bone marrow in response to culture conditions and among cell populations. Multipotent cells isolated from heterogeneous MSC cultures exhibit a greater than three-fold increase in surface expression for NG2 and greater than two-fold increase for CD146 as compared with parental and lineage-committed MSCs. For both antigens, surface expression is downregulated by greater than or equal to six-fold when MSCs become confluent. During serial passage, maximum surface expression of NG2 and CD146 is associated with minimum doubling time. Upregulation of NG2 and CD146 during loss of adipogenic potential at early passage suggests some limits to their utility as potency markers. A potential relationship between proliferation and antigen expression was explored by sorting heterogeneous MSCs into rapidly and slowly dividing groups. Fluorescence-activated cell sorting revealed that rapidly dividing MSCs display lower scatter and 50% higher NG2 surface expression than slowly dividing cells, but CD146 expression is comparable in both groups. Heterogeneous MSCs were sorted based on scatter properties and surface expression of NG2 and CD146 into high (HI) and low (LO) groups. Sc(LO)NG2(HI) and Sc(LO)NG2(HI)CD146(HI) MSCs have the highest proliferative potential of the sorted groups, with colony-forming efficiencies that are 1.5-2.2 times the value for the parental controls. The Sc(LO) gate enriches for rapidly dividing cells. Addition of the NG2(HI) gate increases cell survival to 1.5 times the parental control. Further addition of the CD146(HI) gate does not significantly improve cell division or survival. The combination of low scatter and high NG2 surface expression is a promising selection criterion to enrich a proliferative phenotype from heterogeneous MSCs during ex vivo expansion, with potentially numerous applications.

Prognostic significance of cytoplasmic SOX9 in invasive ductal carcinoma and metastatic breast cancer.

SOX9, a high mobility group (HMG) box transcription factor, is required for development, differentiation and lineage commitment. It is known to exert its effects through nuclear translocation, such as cell cycle changes in response to retinoic acid treatment in breast cancer cells. However, it is not known whether SOX9 has prognostic significance in human breast cancer. Over-expression and cytoplasmic sequestration of nuclear proteins are implicated in tumor progression. To determine whether SOX9 has any prognostic significance in human breast cancer, its expression and subcellular localization were analyzed in more than 200 human breast carcinomas (BCs). SOX9 mRNA expression data for human BCs were computed from microarray studies available in public databases and correlated with known poor prognostic parameters of BCs. SOX9 protein expression and its correlation with Ki-67 staining in human BCs were assessed using immunohistochemistry. Higher SOX9 mRNA levels were significantly associated with estrogen receptor negative (P ≤ 0.001) and higher grade (P ≤ 0.01) human breast tumors. Patients with higher SOX9 mRNA level had significantly shorter overall survival (P ≤ 0.0001). SOX9 protein, which is normally nuclear, was instead localized in the cytoplasm of 25-30% invasive ductal carcinomas (IDCs) and lymph node metastases. Its cytoplasmic accumulation significantly correlated with enhanced proliferation in breast tumors (Kendall's tau = 0.337 with a P value < 0.0001). Cytoplasmic SOX9 can serve as a valuable prognostic marker for IDCs and metastatic breast cancer. Its significant correlation with breast tumor cell proliferation implies that SOX9 directly contributes to the poor clinical outcomes associated with invasive breast cancer.

Decompressive craniectomy for elevated intracranial pressure and its effect on the cumulative ischemic burden and therapeutic intensity levels after severe traumatic brain injury.

BACKGROUND: Increased intracranial pressure (ICP) can cause brain ischemia and compromised brain oxygen (PbtO2 < or = 20 mm Hg) after severe traumatic brain injury (TBI). OBJECTIVE: We examined whether decompressive craniectomy (DC) to treat elevated ICP reduces the cumulative ischemic burden (CIB) of the brain and therapeutic intensity level (TIL). METHODS: Ten severe TBI patients (mean age, 31.4 +/- 14.2 years) who had continuous PbtO2 monitoring before and after delayed DC were retrospectively identified. Patients were managed according to the guidelines for the management of severe TBI. The CIB was measured as the total time spent between a PbtO2 of 15 to 20, 10 to 15, and 0 to 10 mm Hg. The TIL was calculated every 12 hours. Mixed-effects models were used to estimate changes associated with DC. RESULTS: DC was performed on average 2.8 days after admission. DC was found to immediately reduce ICP (mean [SEM] decrease was 7.86 mm Hg [2.4 mm Hg]; P = .005). TIL, which was positively correlated with ICP (r = 0.46, P < or = .001), was reduced within 12 hours after surgery and continued to improve within the postsurgical monitoring period (P

A sharp error probability estimate for the reconstruction of phylogenetic quartets by the four-point method.

Despite the continued development of advanced algorithms for phylogeny reconstruction, the assessment of topological accuracy remains a challenging problem. New tools are needed to assist researchers in the prediction and evaluation of phylogenetic performance, particularly when short alignments are considered. We present a probabilistic analysis of quartet accuracy by the four-point method (FPM) for the Jukes-Cantor model for nucleotide substitution, developing a sharp error estimate as a function of the quartet edge lengths and the number of nucleotide positions available. Our multivariate product (MVP) estimate offers significant improvements over existing bounds and performs well even for short sequence lengths. Results demonstrate that the MVP may be employed in predicting the performance of several popular phylogeny reconstruction methods, and we illustrate a potential application of the MVP to larger datasets using a seven-taxa tree derived from a mitochondrial DNA distance matrix.

Toxoplasma gondii expresses two mitogen-activated protein kinase genes that represent distinct protozoan subfamilies.

All eukaryotes express mitogen-activated protein kinases (MAPKs) that govern diverse cellular processes including proliferation, differentiation, and survival. Even though these proteins are highly conserved throughout nature, MAPKs from closely related species often possess distinct signature sequences, making them well suited as drug discovery targets. Based on the central amino acid in the TXY dual phosphorylation loop, mammalian MAPKs are classified as extracellular signal-regulated kinases (ERKs), c-Jun amino-terminal kinases (JNKs), or p38 stress-response MAPKs. The presence of MAPKs in nonmetazoan eukaryotes suggests significant evolutionary conservation of these important signalling pathways. We recently cloned a novel stress-response MAPK gene (tgMAPK1) from Toxoplasma gondii, an obligate intracellular human parasite that can cause life-threatening infections in immunocompromised patients, and we now present data on a second T. gondii MAPK gene (tgMAPK2) that we cloned. We show that tgMAPK1 and tgMAPK2 are members of two distinct and previously unknown protozoan MAPK subfamilies that we have named pzMAPKl/pzMAPK3 and pzMAPK2. Our phylogenetic analysis of a collection of protozoan and metazoan MAPK genes in relation to ERK8-like genes demonstrates that an ERK8-like family, which includes the pzMAPK2 subfamily, is represented across a large variety of eukaryotic kingdoms and is evolutionarily very distant from other MAPK families.