Thermal niches are more conserved at cold than warm limits in arctic-alpine plant species

Aim Understanding the stability of realised niches is crucial for predicting the responses of species to climate change. One approach is to evaluate the niche differences of populations of the same species that occupy regions that are geographically disconnected. Here, we assess niche conservatism along thermal gradients for 26 plant species with a disjunct distribution between the Alps and the Arctic. Location European Alps and Norwegian Finnmark. Methods We collected a comprehensive dataset of 26 arctic-alpine plant occurrences in two regions. We assessed niche conservatism through a multi-species comparison and analysed species rankings at cold and warm thermal limits along two distinct gradients corresponding to (1) air temperatures at 2 meters above ground level and (2) elevation distances to the treeline (TLD) for the two regions. We assessed whether observed relationships were close to those predicted under thermal limit conservatism. Results We found a weak similarity in species ranking at the warm thermal limits. The range of warm thermal limits for the 26 species was much larger in the Alps than in Finnmark. We found a stronger similarity in species ranking and correspondence at the cold thermal limit along the gradients of 2-m temperature and TLD. Yet, along the 2-m temperature gradient, the cold thermal limits of species in the Alps were lower on average than those in Finnmark. Main conclusion We found low conservatism of the warm thermal limits but a stronger conservatism of the cold thermal limits. We suggest that biotic interactions at the warm thermal limit likely modulate species responses more strongly than at the cold limit. The differing biotic context between the two regions is likely responsible for the observed differences in realised niches.

Species delimitation: A case study in a problematic ant taxon

Species delimitation has been invigorated as a discipline in systematics by an influx of new character sets, analytical methods, and conceptual advances. We use genetic data from 68 markers, combined with distributional, bioclimatic, and coloration information, to hypothesize boundaries of evolutionarily independent lineages (species) within the widespread and highly variable nominal fire ant species Solenopsis saevissima, a member of a species group containing invasive pests as well as species that are models for ecological and evolutionary research. Our integrated approach uses diverse methods of analysis to sequentially test whether populations meet specific operational criteria (contingent properties) for candidacy as morphologically cryptic species, including genetic clustering, monophyly, reproductive isolation, and occupation of distinctive niche space. We hypothesize that nominal S. saevissima comprises at least 4-6 previously unrecognized species, including several pairs whose parapatric distributions implicate the development of intrinsic premating or postmating barriers to gene flow. Our genetic data further suggest that regional genetic differentiation in S. saevissima has been influenced by hybridization with other nominal species occurring in sympatry or parapatry, including the quite distantly related Solenopsis geminata. The results of this study illustrate the importance of employing different classes of genetic data (coding and noncoding regions and nuclear and mitochondrial DNA [mtDNA] markers), different methods of genetic data analysis (tree-based and non-tree based methods), and different sources of data (genetic, morphological, and ecological data) to explicitly test various operational criteria for species boundaries in clades of recently diverged lineages, while warning against over reliance on any single data type (e.g., mtDNA sequence variation) when drawing inferences.

A molecular framework for light and gibberellin control of cell elongation.

Cell elongation during seedling development is antagonistically regulated by light and gibberellins (GAs). Light induces photomorphogenesis, leading to inhibition of hypocotyl growth, whereas GAs promote etiolated growth, characterized by increased hypocotyl elongation. The mechanism underlying this antagonistic interaction remains unclear. Here we report on the central role of the Arabidopsis thaliana nuclear transcription factor PIF4 (encoded by PHYTOCHROME INTERACTING FACTOR 4) in the positive control of genes mediating cell elongation and show that this factor is negatively regulated by the light photoreceptor phyB (ref. 4) and by DELLA proteins that have a key repressor function in GA signalling. Our results demonstrate that PIF4 is destabilized by phyB in the light and that DELLAs block PIF4 transcriptional activity by binding the DNA-recognition domain of this factor. We show that GAs abrogate such repression by promoting DELLA destabilization, and therefore cause a concomitant accumulation of free PIF4 in the nucleus. Consistent with this model, intermediate hypocotyl lengths were observed in transgenic plants over-accumulating both DELLAs and PIF4. Destabilization of this factor by phyB, together with its inactivation by DELLAs, constitutes a protein interaction framework that explains how plants integrate both light and GA signals to optimize growth and development in response to changing environments.

Warming experiments underpredict plant phenological responses to climate change.

Warming experiments are increasingly relied on to estimate plant responses to global climate change. For experiments to provide meaningful predictions of future responses, they should reflect the empirical record of responses to temperature variability and recent warming, including advances in the timing of flowering and leafing. We compared phenology (the timing of recurring life history events) in observational studies and warming experiments spanning four continents and 1,634 plant species using a common measure of temperature sensitivity (change in days per degree Celsius). We show that warming experiments underpredict advances in the timing of flowering and leafing by 8.5-fold and 4.0-fold, respectively, compared with long-term observations. For species that were common to both study types, the experimental results did not match the observational data in sign or magnitude. The observational data also showed that species that flower earliest in the spring have the highest temperature sensitivities, but this trend was not reflected in the experimental data. These significant mismatches seem to be unrelated to the study length or to the degree of manipulated warming in experiments. The discrepancy between experiments and observations, however, could arise from complex interactions among multiple drivers in the observational data, or it could arise from remediable artefacts in the experiments that result in lower irradiance and drier soils, thus dampening the phenological responses to manipulated warming. Our results introduce uncertainty into ecosystem models that are informed solely by experiments and suggest that responses to climate change that are predicted using such models should be re-evaluated.

Protein-protein interaction investigated by steered molecular dynamics: the TCR-pMHC complex.

We present a novel steered molecular dynamics scheme to induce the dissociation of large protein-protein complexes. We apply this scheme to study the interaction of a T cell receptor (TCR) with a major histocompatibility complex (MHC) presenting a peptide (p). Two TCR-pMHC complexes are considered, which only differ by the mutation of a single amino acid on the peptide; one is a strong agonist that produces T cell activation in vivo, while the other is an antagonist. We investigate the interaction mechanism from a large number of unbinding trajectories by analyzing van der Waals and electrostatic interactions and by computing energy changes in proteins and solvent. In addition, dissociation potentials of mean force are calculated with the Jarzynski identity, using an averaging method developed for our steering scheme. We analyze the convergence of the Jarzynski exponential average, which is hampered by the large amount of dissipative work involved and the complexity of the system. The resulting dissociation free energies largely underestimate experimental values, but the simulations are able to clearly differentiate between wild-type and mutated TCR-pMHC and give insights into the dissociation mechanism.

Structural and functional interaction sites between Na,K-ATPase and FXYD proteins.

Several members of the FXYD protein family are tissue-specific regulators of Na,K-ATPase that produce distinct effects on its apparent K(+) and Na(+) affinity. Little is known about the interaction sites between the Na,K-ATPase alpha subunit and FXYD proteins that mediate the efficient association and/or the functional effects of FXYD proteins. In this study, we have analyzed the role of the transmembrane segment TM9 of the Na,K-ATPase alpha subunit in the structural and functional interaction with FXYD2, FXYD4, and FXYD7. Mutational analysis combined with expression in Xenopus oocytes reveals that Phe(956), Glu(960), Leu(964), and Phe(967) in TM9 of the Na,K-ATPase alpha subunit represent one face interacting with the three FXYD proteins. Leu(964) and Phe(967) contribute to the efficient association of FXYD proteins with the Na,K-ATPase alpha subunit, whereas Phe(956) and Glu(960) are essential for the transmission of the functional effect of FXYD proteins on the apparent K(+) affinity of Na,K-ATPase. The relative contribution of Phe(956) and Glu(960) to the K(+) effect differs for different FXYD proteins, probably reflecting the intrinsic differences of FXYD proteins on the apparent K(+) affinity of Na,K-ATPase. In contrast to the effect on the apparent K(+) affinity, Phe(956) and Glu(960) are not involved in the effect of FXYD2 and FXYD4 on the apparent Na(+) affinity of Na,K-ATPase. The mutational analysis is in good agreement with a docking model of the Na,K-ATPase/FXYD7 complex, which also predicts the importance of Phe(956), Glu(960), Leu(964), and Phe(967) in subunit interaction. In conclusion, by using mutational analysis and modeling, we show that TM9 of the Na,K-ATPase alpha subunit exposes one face of the helix that interacts with FXYD proteins and contributes to the stable interaction with FXYD proteins, as well as mediating the effect of FXYD proteins on the apparent K(+) affinity of Na,K-ATPase.

Phylogenetic conservatism in plant phenology

Phenological events - defined points in the life cycle of a plant or animal - have been regarded as highly plastic traits, reflecting flexible responses to various environmental cues. The ability of a species to track, via shifts in phenological events, the abiotic environment through time might dictate its vulnerability to future climate change. Understanding the predictors and drivers of phenological change is therefore critical. Here, we evaluated evidence for phylogenetic conservatism - the tendency for closely related species to share similar ecological and biological attributes - in phenological traits across flowering plants. We aggregated published and unpublished data on timing of first flower and first leaf, encompassing 4000 species at 23 sites across the Northern Hemisphere. We reconstructed the phylogeny for the set of included species, first, using the software program Phylomatic, and second, from DNA data. We then quantified phylogenetic conservatism in plant phenology within and across sites. We show that more closely related species tend to flower and leaf at similar times. By contrasting mean flowering times within and across sites, however, we illustrate that it is not the time of year that is conserved, but rather the phenological responses to a common set of abiotic cues. Our findings suggest that species cannot be treated as statistically independent when modelling phenological responses.Synthesis. Closely related species tend to resemble each other in the timing of their life-history events, a likely product of evolutionarily conserved responses to environmental cues. The search for the underlying drivers of phenology must therefore account for species' shared evolutionary histories.

Three-dimensional models of 14α-sterol demethylase (Cyp51A) from Aspergillus lentulus and Aspergillus fumigatus: an insight into differences in voriconazole interaction.

Aspergillus lentulus, an Aspergillus fumigatus sibling species, is increasingly reported in corticosteroid-treated patients. Its clinical significance is unknown, but the fact that A. lentulus shows reduced antifungal susceptibility, mainly to voriconazole, is of serious concern. Heterologous expression of cyp51A from A. fumigatus and A. lentulus was performed in Saccharomyces cerevisiae to assess differences in the interaction of Cyp51A with the azole drugs. The absence of endogenous ERG11 was efficiently complemented in S. cerevisiae by the expression of either Aspergillus cyp51A allele. There was a marked difference between azole minimum inhibitory concentration (MIC) values of the clones expressing each Aspergillus spp. cyp51A. Saccharomyces cerevisiae clones expressing A. lentulus alleles showed higher MICs to all of the azoles tested, supporting the hypothesis that the intrinsic azole resistance of A. lentulus could be associated with Cyp51A. Homology models of A. fumigatus and A. lentulus Cyp51A protein based on the crystal structure of Cyp51p from Mycobacterium tuberculosis in complex with fluconazole were almost identical owing to their mutual high sequence identity. Molecular dynamics (MD) was applied to both three-dimensional protein models to refine the homology modelling and to explore possible differences in the Cyp51A-voriconazole interaction. After 20ns of MD modelling, some critical differences were observed in the putative closed form adopted by the protein upon voriconazole binding. A closer study of the A. fumigatus and A. lentulus voriconazole putative binding site in Cyp51A suggested that some major differences in the protein's BC loop could differentially affect the lock-up of voriconazole, which in turn could correlate with their different azole susceptibility profiles.

LILRB2 interaction with HLA class I correlates with control of HIV-1 infection.

Natural progression of HIV-1 infection depends on genetic variation in the human major histocompatibility complex (MHC) class I locus, and the CD8+ T cell response is thought to be a primary mechanism of this effect. However, polymorphism within the MHC may also alter innate immune activity against human immunodeficiency virus type 1 (HIV-1) by changing interactions of human leukocyte antigen (HLA) class I molecules with leukocyte immunoglobulin-like receptors (LILR), a group of immunoregulatory receptors mainly expressed on myelomonocytic cells including dendritic cells (DCs). We used previously characterized HLA allotype-specific binding capacities of LILRB1 and LILRB2 as well as data from a large cohort of HIV-1-infected individuals (N = 5126) to test whether LILR-HLA class I interactions influence viral load in HIV-1 infection. Our analyses in persons of European descent, the largest ethnic group examined, show that the effect of HLA-B alleles on HIV-1 control correlates with the binding strength between corresponding HLA-B allotypes and LILRB2 (p = 10(-2)). Moreover, overall binding strength of LILRB2 to classical HLA class I allotypes, defined by the HLA-A/B/C genotypes in each patient, positively associates with viral replication in the absence of therapy in patients of both European (p = 10(-11)-10(-9)) and African (p = 10(-5)-10(-3)) descent. This effect appears to be driven by variations in LILRB2 binding affinities to HLA-B and is independent of individual class I allelic effects that are not related to the LILRB2 function. Correspondingly, in vitro experiments suggest that strong LILRB2-HLA binding negatively affects antigen-presenting properties of DCs. Thus, we propose an impact of LILRB2 on HIV-1 disease outcomes through altered regulation of DCs by LILRB2-HLA engagement.

The degradation of HFR1, a putative bHLH class transcription factor involved in light signaling, is regulated by phosphorylation and requires COP1.

All developmental transitions throughout the life cycle of a plant are influenced by light. In Arabidopsis, multiple photoreceptors including the UV-A/blue-sensing cryptochromes (cry1-2) and the red/far-red responsive phytochromes (phyA-E) monitor the ambient light conditions. Light-regulated protein stability is a major control point of photomorphogenesis. The ubiquitin E3 ligase COP1 (constitutively photomorphogenic 1) regulates the stability of several light-signaling components. HFR1 (long hypocotyl in far-red light) is a putative transcription factor with a bHLH domain acting downstream of both phyA and the cryptochromes. HFR1 is closely related to PIF1, PIF3, and PIF4 (phytochrome interacting factor 1, 3 and 4), but in contrast to the latter three, there is no evidence for a direct interaction between HFR1 and the phytochromes. Here, we show that the protein abundance of HFR1 is tightly controlled by light. HFR1 is an unstable phosphoprotein, particularly in the dark. The proteasome and COP1 are required in vivo to degrade phosphorylated HFR1. In addition, HFR1 can interact with COP1, consistent with the idea of COP1 directly mediating HFR1 degradation. We identify a domain, conserved among several bHLH class proteins involved in light signaling , as a determinant of HFR1 stability. Our physiological experiments indicate that the control of HFR1 protein abundance is important for a normal de-etiolation response.

Genomics of the divergence continuum in an African plant biodiversity hotspot, I: drivers of population divergence in Restio capensis (Restionaceae).

Understanding the drivers of population divergence, speciation and species persistence is of great interest to molecular ecology, especially for species-rich radiations inhabiting the world's biodiversity hotspots. The toolbox of population genomics holds great promise for addressing these key issues, especially if genomic data are analysed within a spatially and ecologically explicit context. We have studied the earliest stages of the divergence continuum in the Restionaceae, a species-rich and ecologically important plant family of the Cape Floristic Region (CFR) of South Africa, using the widespread CFR endemic Restio capensis (L.) H.P. Linder & C.R. Hardy as an example. We studied diverging populations of this morphotaxon for plastid DNA sequences and >14 400 nuclear DNA polymorphisms from Restriction site Associated DNA (RAD) sequencing and analysed the results jointly with spatial, climatic and phytogeographic data, using a Bayesian generalized linear mixed modelling (GLMM) approach. The results indicate that population divergence across the extreme environmental mosaic of the CFR is mostly driven by isolation by environment (IBE) rather than isolation by distance (IBD) for both neutral and non-neutral markers, consistent with genome hitchhiking or coupling effects during early stages of divergence. Mixed modelling of plastid DNA and single divergent outlier loci from a Bayesian genome scan confirmed the predominant role of climate and pointed to additional drivers of divergence, such as drift and ecological agents of selection captured by phytogeographic zones. Our study demonstrates the usefulness of population genomics for disentangling the effects of IBD and IBE along the divergence continuum often found in species radiations across heterogeneous ecological landscapes.

Mechanisms of reproductive isolation between an ant species of hybrid origin and one of its parents.

The establishment of new species by hybridization is difficult because it requires the development of reproductive isolation (RI) in sympatry to escape the homogenizing effects of gene flow from the parental species. Here we investigated the role of two pre- and two postzygotic mechanisms of RI in a system comprising two interdependent Pogonomyrmex harvester ant lineages (the H1 and H2 lineages) of hybrid origin and one of their parental species (P. rugosus). Similar to most other ants, P. rugosus is characterized by an environmental system of caste determination with female brood developing either into queens or workers depending on nongenetic factors. By contrast, there is a strong genetic component to caste determination in the H1 and H2 lineages because the developmental fate of female brood depends on the genetic origin of the parents, with interlineage eggs developing into workers and intralineage eggs developing into queens. The study of a mixed mating aggregation revealed strong differences in mating flight timing between P. rugosus and the two lineages as a first mechanism of RI. A second important prezygotic mechanism was assortative mating. Laboratory experiments also provided support for one of the two investigated mechanisms of postzygotic isolation. The majority of offspring produced from the few matings between P. rugosus and the lineages aborted at the egg stage. This hybrid inviability was under maternal influence, with hybrids produced by P. rugosus queens being always inviable whereas a small proportion of H2 lineage queens produced large numbers of adult hybrid offspring. Finally, we found no evidence that genetic caste determination acted as a second postzygotic mechanism reducing gene flow between P. rugosus and the H lineages. The few viable P. rugosus-H hybrids were not preferentially shunted into functionally sterile workers but developed into both workers and queens. Overall, these results reveal that the nearly complete (99.5%) RI between P. rugosus and the two hybrid lineages stems from the combination of two typical prezygotic mechanisms (mating time divergence and assortative mating) and one postzygotic mechanism (hybrid inviability).

Systems biology of plants : from intraspecific genome variation to computational morphodynamics

Recently, the introduction of second generation sequencing and further advance-ments in confocal microscopy have enabled system-level studies for the functional characterization of genes. The degree of complexity intrinsic to these approaches needs the development of bioinformatics methodologies and computational models for extracting meaningful biological knowledge from the enormous amount of experi¬mental data which is continuously generated. This PhD thesis presents several novel bioinformatics methods and computational models to address specific biological questions in Plant Biology by using the plant Arabidopsis thaliana as a model system. First, a spatio-temporal qualitative analysis of quantitative transcript and protein profiles is applied to show the role of the BREVIS RADIX (BRX) protein in the auxin- cytokinin crosstalk for root meristem growth. Core of this PhD work is the functional characterization of the interplay between the BRX protein and the plant hormone auxin in the root meristem by using a computational model based on experimental evidence. Hyphotesis generated by the modelled to the discovery of a differential endocytosis pattern in the root meristem that splits the auxin transcriptional response via the plasma membrane to nucleus partitioning of BRX. This positional information system creates an auxin transcriptional pattern that deviates from the canonical auxin response and is necessary to sustain the expression of a subset of BRX-dependent auxin-responsive genes to drive root meristem growth. In the second part of this PhD thesis, we characterized the genome-wide impact of large scale deletions on four divergent Arabidopsis natural strains, through the integration of Ultra-High Throughput Sequencing data with data from genomic hybridizations on tiling arrays. Analysis of the identified deletions revealed a considerable portion of protein coding genes affected and supported a history of genomic rearrangements shaped by evolution. In the last part of the thesis, we showed that VIP3 gene in Arabidopsis has an evo-lutionary conserved role in the 3' to 5' mRNA degradation machinery, by applying a novel approach for the analysis of mRNA-Seq data from random-primed mRNA. Altogether, this PhD research contains major advancements in the study of natural genomic variation in plants and in the application of computational morphodynamics models for the functional characterization of biological pathways essential for the plant. - Récemment, l'introduction du séquençage de seconde génération et les avancées dans la microscopie confocale ont permis des études à l'échelle des différents systèmes cellulaires pour la caractérisation fonctionnelle de gènes. Le degrés de complexité intrinsèque à ces approches ont requis le développement de méthodologies bioinformatiques et de modèles mathématiques afin d'extraire de la masse de données expérimentale générée, des information biologiques significatives. Ce doctorat présente à la fois des méthodes bioinformatiques originales et des modèles mathématiques pour répondre à certaines questions spécifiques de Biologie Végétale en utilisant la plante Arabidopsis thaliana comme modèle. Premièrement, une analyse qualitative spatio-temporelle de profiles quantitatifs de transcripts et de protéines est utilisée pour montrer le rôle de la protéine BREVIS RADIX (BRX) dans le dialogue entre l'auxine et les cytokinines, des phytohormones, dans la croissance du méristème racinaire. Le noyau de ce travail de thèse est la caractérisation fonctionnelle de l'interaction entre la protéine BRX et la phytohormone auxine dans le méristème de la racine en utilisant des modèles informatiques basés sur des preuves expérimentales. Les hypothèses produites par le modèle ont mené à la découverte d'un schéma différentiel d'endocytose dans le méristème racinaire qui divise la réponse transcriptionnelle à l'auxine par le partitionnement de BRX de la membrane plasmique au noyau de la cellule. Cette information positionnelle crée une réponse transcriptionnelle à l'auxine qui dévie de la réponse canonique à l'auxine et est nécessaire pour soutenir l'expression d'un sous ensemble de gènes répondant à l'auxine et dépendant de BRX pour conduire la croissance du méristème. Dans la seconde partie de cette thèse de doctorat, nous avons caractérisé l'impact sur l'ensemble du génome des délétions à grande échelle sur quatre souches divergentes naturelles d'Arabidopsis, à travers l'intégration du séquençage à ultra-haut-débit avec l'hybridation génomique sur puces ADN. L'analyse des délétions identifiées a révélé qu'une proportion considérable de gènes codant était affectée, supportant l'idée d'un historique de réarrangement génomique modelé durant l'évolution. Dans la dernière partie de cette thèse, nous avons montré que le gène VÏP3 dans Arabidopsis a conservé un rôle évolutif dans la machinerie de dégradation des ARNm dans le sens 3' à 5', en appliquant une nouvelle approche pour l'analyse des données de séquençage d'ARNm issue de transcripts amplifiés aléatoirement. Dans son ensemble, cette recherche de doctorat contient des avancées majeures dans l'étude des variations génomiques naturelles des plantes et dans l'application de modèles morphodynamiques informatiques pour la caractérisation de réseaux biologiques essentiels à la plante. - Le développement des plantes est écrit dans leurs codes génétiques. Pour comprendre comment les plantes sont capables de s'adapter aux changements environnementaux, il est essentiel d'étudier comment leurs gènes gouvernent leur formation. Plus nous essayons de comprendre le fonctionnement d'une plante, plus nous réalisons la complexité des mécanismes biologiques, à tel point que l'utilisation d'outils et de modèles mathématiques devient indispensable. Dans ce travail, avec l'utilisation de la plante modèle Arabidopsis thalicinci nous avons résolu des problèmes biologiques spécifiques à travers le développement et l'application de méthodes informatiques concrètes. Dans un premier temps, nous avons investigué comment le gène BREVIS RADIX (BRX) régule le développement de la racine en contrôlant la réponse à deux hormones : l'auxine et la cytokinine. Nous avons employé une analyse statistique sur des mesures quantitatives de transcripts et de produits de gènes afin de démontrer que BRX joue un rôle antagonisant dans le dialogue entre ces deux hormones. Lorsque ce-dialogue moléculaire est perturbé, la racine primaire voit sa longueur dramatiquement réduite. Pour comprendre comment BRX répond à l'auxine, nous avons développé un modèle informatique basé sur des résultats expérimentaux. Les simulations successives ont mené à la découverte d'un signal positionnel qui contrôle la réponse de la racine à l'auxine par la régulation du mouvement intracellulaire de BRX. Dans la seconde partie de cette thèse, nous avons analysé le génome entier de quatre souches naturelles d'Arabidopsis et nous avons trouvé qu'une grande partie de leurs gènes étaient manquant par rapport à la souche de référence. Ce résultat indique que l'historique des modifications génomiques conduites par l'évolution détermine une disponibilité différentielle des gènes fonctionnels dans ces plantes. Dans la dernière partie de ce travail, nous avons analysé les données du transcriptome de la plante où le gène VIP3 était non fonctionnel. Ceci nous a permis de découvrir le rôle double de VIP3 dans la régulation de l'initiation de la transcription et dans la dégradation des transcripts. Ce rôle double n'avait jusqu'alors été démontrée que chez l'homme. Ce travail de doctorat supporte le développement et l'application de méthodologies informatiques comme outils inestimables pour résoudre la complexité des problèmes biologiques dans la recherche végétale. L'intégration de la biologie végétale et l'informatique est devenue de plus en plus importante pour l'avancée de nos connaissances sur le fonctionnement et le développement des plantes.

Pharmacokinetic interaction between methotrexate and chloral hydrate.

We report the case of a drug interaction between methotrexate (MTX) and chloral hydrate (CH) observed in a child treated for acute leukemia. Significantly slower MTX clearance and increased MTX exposure occurred on the first three courses of a high-dose chemotherapy when co-administered with CH despite normal renal function, adequate hydration, and alkalinization. Mean MTX area under the curve associated with CH administration was 1,134 µmol hours/L, compared to 608 µmol hours/L after discontinuation of CH. This interaction possibly resulted from a competition between anionic CH metabolites and MTX for renal tubular excretion.