Molecular markers allow sibling species identification in red wood ants (Formica rufa group)

Red wood ants (Formica rufa group) constitute a group of species that are considered to be among the most promising bioindicators in forest ecosystems. However, because of their morphological similarity and intraspecific variability, morphological species identification can be difficult. Considerable expertise is necessary to discriminate between the sibling species F. lugubris and F. paralugubris, two species that often live in sympatry in the same Alpine forests. New taxonomic tools providing rapid and reliable species identification are needed. We present a simple and reliable molecular technique based on mtDNA (COI gene) and a restriction enzyme for discriminating between F. lugubris and F. paralugubris. We confirm the validity of this method with a Bayesian analysis based on microsatellites. This new molecular tool represents a clear breakthrough for discriminating between F. lugubris and F. paralugubris and is likely to be helpful in large-scale biomonitoring.

A novel method for automated classification of epileptiform activity in the human electroencephalogram-based on independent component analysis.

Diagnosis of several neurological disorders is based on the detection of typical pathological patterns in the electroencephalogram (EEG). This is a time-consuming task requiring significant training and experience. Automatic detection of these EEG patterns would greatly assist in quantitative analysis and interpretation. We present a method, which allows automatic detection of epileptiform events and discrimination of them from eye blinks, and is based on features derived using a novel application of independent component analysis. The algorithm was trained and cross validated using seven EEGs with epileptiform activity. For epileptiform events with compensation for eyeblinks, the sensitivity was 65 +/- 22% at a specificity of 86 +/- 7% (mean +/- SD). With feature extraction by PCA or classification of raw data, specificity reduced to 76 and 74%, respectively, for the same sensitivity. On exactly the same data, the commercially available software Reveal had a maximum sensitivity of 30% and concurrent specificity of 77%. Our algorithm performed well at detecting epileptiform events in this preliminary test and offers a flexible tool that is intended to be generalized to the simultaneous classification of many waveforms in the EEG.

Subject-based steroid profiling and the determination of novel biomarkers for DHT and DHEA misuse in sports.

Doping with natural steroids can be detected by evaluating the urinary concentrations and ratios of several endogenous steroids. Since these biomarkers of steroid doping are known to present large inter-individual variations, monitoring of individual steroid profiles over time allows switching from population-based towards subject-based reference ranges for improved detection. In an Athlete Biological Passport (ABP), biomarkers data are collated throughout the athlete's sporting career and individual thresholds defined adaptively. For now, this approach has been validated on a limited number of markers of steroid doping, such as the testosterone (T) over epitestosterone (E) ratio to detect T misuse in athletes. Additional markers are required for other endogenous steroids like dihydrotestosterone (DHT) and dehydroepiandrosterone (DHEA). By combining comprehensive steroid profiles composed of 24 steroid concentrations with Bayesian inference techniques for longitudinal profiling, a selection was made for the detection of DHT and DHEA misuse. The biomarkers found were rated according to relative response, parameter stability, discriminative power, and maximal detection time. This analysis revealed DHT/E, DHT/5β-androstane-3α,17β-diol and 5α-androstane-3α,17β-diol/5β-androstane-3α,17β-diol as best biomarkers for DHT administration and DHEA/E, 16α-hydroxydehydroepiandrosterone/E, 7β-hydroxydehydroepiandrosterone/E and 5β-androstane-3α,17β-diol/5α-androstane-3α,17β-diol for DHEA. The selected biomarkers were found suitable for individual referencing. A drastic overall increase in sensitivity was obtained. The use of multiple markers as formalized in an Athlete Steroidal Passport (ASP) can provide firm evidence of doping with endogenous steroids.

Biochemical markers of fatal hypothermia.

The aim of this study was to investigate the usefulness of postmortem biochemical investigations in the diagnosis of fatal hypothermia. 10 cases of fatal hypothermia and 30 control cases were selected. A series of biochemical parameters, such as glucose, acetone, 3-beta-hydroxybutyrate, isopropyl alcohol, free fatty acids, adrenaline, growth hormone, adrenocorticotropic hormone, thyroid-stimulating hormone, cortisol, calcium, magnesium, C-reactive protein, procalcitonin as well as markers of renal and cardiac functions were measured in blood, postmortem serum from femoral blood, urine, vitreous and pericardial fluid. The results suggested that deaths due to hypothermia, especially in free-ethanol cases, are characterized by increased ketone levels in blood and other biological fluids, increased adrenaline concentrations in urine, increased cortisol levels in postmortem serum from femoral blood and increased free cortisol values in urine. Increased or decreased levels of other biological parameters are either the result of terminal metabolic changes or the expression of preexisting diseases and may provide information to elucidate the death process on a case-by-case basis.

Functional analysis: evaluation of response intensities--tailoring ANOVA for lists of expression subsets.

Background: Microarray data is frequently used to characterize the expression profile of a whole genome and to compare the characteristics of that genome under several conditions. Geneset analysis methods have been described previously to analyze the expression values of several genes related by known biological criteria (metabolic pathway, pathology signature, co-regulation by a common factor, etc.) at the same time and the cost of these methods allows for the use of more values to help discover the underlying biological mechanisms. Results: As several methods assume different null hypotheses, we propose to reformulate the main question that biologists seek to answer. To determine which genesets are associated with expression values that differ between two experiments, we focused on three ad hoc criteria: expression levels, the direction of individual gene expression changes (up or down regulation), and correlations between genes. We introduce the FAERI methodology, tailored from a two-way ANOVA to examine these criteria. The significance of the results was evaluated according to the self-contained null hypothesis, using label sampling or by inferring the null distribution from normally distributed random data. Evaluations performed on simulated data revealed that FAERI outperforms currently available methods for each type of set tested. We then applied the FAERI method to analyze three real-world datasets on hypoxia response. FAERI was able to detect more genesets than other methodologies, and the genesets selected were coherent with current knowledge of cellular response to hypoxia. Moreover, the genesets selected by FAERI were confirmed when the analysis was repeated on two additional related datasets. Conclusions: The expression values of genesets are associated with several biological effects. The underlying mathematical structure of the genesets allows for analysis of data from several genes at the same time. Focusing on expression levels, the direction of the expression changes, and correlations, we showed that two-step data reduction allowed us to significantly improve the performance of geneset analysis using a modified two-way ANOVA procedure, and to detect genesets that current methods fail to detect.

Prospective analysis of KRAS, BRAF, and PIK3CA mutations and EGFR copy number in patients (pts) with locally advanced rectal cancer: A translational substudy of a clinical trial (SAKK 41/07) evaluating the effect of neoadjuvant chemoradiation (CRT) with or without panitumumab

Background: Prognostic and predictive markers are of great importance for future study designs and essential for the interpretation of clinical trials incorporating an EGFR-inhibitor. The current study prospectively assessed and validated KRAS, BRAF and PIK3CA mutations in rectal cancer patients screened for the trial SAKK41/07 of concomitant preoperative radio-chemotherapy with or without panitumumab.Methods: Macrodissection was performed on pretreatment formalin fixed paraffin embedded biopsy tissue sections to arrive at a minimum of 50% of tumor cells. DNA was extracted with the Maxwell 16 FFPE Tissue LEV DNA purification kit. After PCR amplification, mutations were identified by pyrosequencing. We prospectively analysed pretreatment biopsy material from 149 rectal cancer pts biopsies for KRAS (exon 2 codon 12 [2-12] and 13 [2-13], exon 3 codon 59 [3-59]) and 61 [3-61], exon 4 codon 117 [4-117] and 146 [4-146]). Sixty-eight pts (KRASwt exon 2, 3 only) were further analysed for BRAF (exon 15 codon 600) and PIK3CA (exon 9 codon 542, 545 and 546, exon 20 codon 1043 [20-1043] and 1047 [20-1047]) mutations, and EGFR copy number by qPCR. For the calculation of the EGFR copy number, we used KRAS copy number as internal reference standard. The calculation was done on the basis of the two standard curves relative quantification method.Results: In 149 screened pts with rectal cancer, the prevalence of KRAS mutations was 36%. Among the 68 pts enrolled in SAKK 41/07 based on initially presumed KRASwt status (exon 2/codons 12+13), 18 pts (26%) had a total of 23 mutations in the RAS/PIK3CA-pathways upon validation analysis. Twelve pts had a KRAS mutation, 7 pts had a PIK3CA mutation, 3 pts had a NRAS mutation, 1 patient a BRAF mutation. Surprisingly, five of these pts had double- mutations, including 4 pts with KRAS plus PIK3CA mutations, and 1 pt with NRAS plus PIK3CA mutations. The median normalized EGFR copy number was 1. Neither mutations of KRAS, BRAF, and PIK3CA, nor EGFR copy number were statistically associated with the primary study endpoint pCR (pathological complete regression).Conclusions: The prevalence of KRAS mutations in rectal and in colon cancer appears to be similar. BRAF mutations are rare; PIK3CA mutations are more common (10%). EGFR copy number is not increased in rectal cancer. A considerable number or KRAS exon 2 wt tumors harbored KRAS exon 3+4 mutations. Further study is needed to determine if KRAS testing should include exons 2-4.

A 2D/3D image analysis system to track fluorescently labeled structures in rod-shaped cells: application to measure spindle pole asymmetry during mitosis.

BACKGROUND: The yeast Schizosaccharomyces pombe is frequently used as a model for studying the cell cycle. The cells are rod-shaped and divide by medial fission. The process of cell division, or cytokinesis, is controlled by a network of signaling proteins called the Septation Initiation Network (SIN); SIN proteins associate with the SPBs during nuclear division (mitosis). Some SIN proteins associate with both SPBs early in mitosis, and then display strongly asymmetric signal intensity at the SPBs in late mitosis, just before cytokinesis. This asymmetry is thought to be important for correct regulation of SIN signaling, and coordination of cytokinesis and mitosis. In order to study the dynamics of organelles or large protein complexes such as the spindle pole body (SPB), which have been labeled with a fluorescent protein tag in living cells, a number of the image analysis problems must be solved; the cell outline must be detected automatically, and the position and signal intensity associated with the structures of interest within the cell must be determined. RESULTS: We present a new 2D and 3D image analysis system that permits versatile and robust analysis of motile, fluorescently labeled structures in rod-shaped cells. We have designed an image analysis system that we have implemented as a user-friendly software package allowing the fast and robust image-analysis of large numbers of rod-shaped cells. We have developed new robust algorithms, which we combined with existing methodologies to facilitate fast and accurate analysis. Our software permits the detection and segmentation of rod-shaped cells in either static or dynamic (i.e. time lapse) multi-channel images. It enables tracking of two structures (for example SPBs) in two different image channels. For 2D or 3D static images, the locations of the structures are identified, and then intensity values are extracted together with several quantitative parameters, such as length, width, cell orientation, background fluorescence and the distance between the structures of interest. Furthermore, two kinds of kymographs of the tracked structures can be established, one representing the migration with respect to their relative position, the other representing their individual trajectories inside the cell. This software package, called "RodCellJ", allowed us to analyze a large number of S. pombe cells to understand the rules that govern SIN protein asymmetry. CONCLUSIONS: "RodCell" is freely available to the community as a package of several ImageJ plugins to simultaneously analyze the behavior of a large number of rod-shaped cells in an extensive manner. The integration of different image-processing techniques in a single package, as well as the development of novel algorithms does not only allow to speed up the analysis with respect to the usage of existing tools, but also accounts for higher accuracy. Its utility was demonstrated on both 2D and 3D static and dynamic images to study the septation initiation network of the yeast Schizosaccharomyces pombe. More generally, it can be used in any kind of biological context where fluorescent-protein labeled structures need to be analyzed in rod-shaped cells. AVAILABILITY: RodCellJ is freely available under http://bigwww.epfl.ch/algorithms.html, (after acceptance of the publication).

Functional analysis of different haematopoietic stem cell subsets

RÉSUMÉ : Elucider les bases moléculaires et cellulaires du fonctionnement des cellules souches s'avère crucial dans la compréhension de l'organisation cellulaire au sein des tissus et des organes ainsi que pour le développement de nouvelles stratégies thérapeutiques en médecine régénérative et en oncologie. Les cellules souches adultes les mieux connues sont celles responsables de l'hématopoïèse, les cellules souches hématopoïétiques (CSH). Durant ces dernières années, la recherche a porté une attention particulière à l'isolation prospective de CSH dérivées de la moelle osseuse de souris en utilisant des marqueurs de surface cellulaire ainsi que des propriétés fonctionnelles alléguées. Par la suite, la capacité fonctionnelle des CSH a été vérifiée classiquement par leur transplantation intraveineuse dans des souris réceptrices conditionnées et par l'analyse de leur aptitude à reconstituer le système hématopoïétique à long terme. Des études récentes suggèrent que la transplantation des cellules directement dans la moelle osseuse pourrait non seulement aboutir à une prise de greffe plus rapide et plus efficace, mais pourrait même aider à l'identification de cellules qui ont certes des propriétés intrinsèques de CSH, mais qui n'ont pas la capacité de trouver leur niche au sein de la moelle osseuse et ont donc échoué dans les analyses classiques de reconstitution. Dans cette étude, nous comparons à deux niveaux la fonction de différents sous-groupes de cellules souches de la moelle osseuse, définis par leur phénotype de surface cellulaire. Premièrement, nous étudions leur capacité à reconstituer des souris létalement irradiées après injection soit intraveineuse soit intrafémorale. Deuxièmement, par analyse cytométrique de flux à 8 couleurs, nous comparons leur activité relative de « side population » (SP) par exclusion du colorant fluorescent Hoechst 33342. Nos résultats préliminaires renforcent en effet l'idée que la transplantation intrafémorale aboutit à une greffe plus rapide et plus efficace. Par contre, en utilisant cette approche, nous n'arrivons pas à identifier des cellules capables de prendre greffe spécifiquement quand elles sont injectées en intrafémorale. Finalement, bien qu'une confirmation in vivo soit encore nécessaire, nous suggérons sur la base de nos analyses cytométriques de flux, que les cellules SP Sca1t~és éie~~ CD48t~és bas sont très enrichies en CSH. Ceci permettrait l'isolation ex vivo de CSH de la moelle osseuse de souris par une stratégie à la fois nouvelle et simple. SUMMARY : Elucidating the molecular and cellular bases of stem cell function is crucial for the understanding of cellular organisation within tissues and organs as well as for the development of new therapeutic strategies in regenerative medicine and oncology. The best-known adult stem cells are those responsible for haematopoiesis, the haematopoietic stem cells (HSCs). In recent years, much effort has been put into the prospective isolation of mouse bone marrow (BM)-derived HSCs using cell-surface markers and alleged functional properties. Upon isolation, the functional capacity of putative HSCs has been classically assessed by intravenous transplantation into conditioned recipient mice and analysis of their ability to reconstitute the haematopoietic system at long-term. It has recently been suggested that transplanting the cells directly into the BM might not only result in more rapid and more effective engraftment, but even help to identify cells that have intrinsic HSC properties but lack the ability to home to their BM niche and have thus failed to succeed in classical reconstitution assays. In this study, we compare the function of different BM cell subsets, as defined by their cell surface phenotype, on two levels. Firstly, we assess their ability to reconstitute lethally irradiated mice, when injected either intravenously or intrafemorally. Secondly, using 8-colour flow cytometric analysis, we compare their relative side population (SP) activity by exclusion of the fluorescent dye Hoechst 33342. Our preliminary results indeed reinforce the idea that intrafemoral transplantation results in faster and more effective engraftment, however, using this approach, we are unable to identify cells that are capable of engrafting specifically when injected intrafemorally. Finally, although in vivo confirmation is still required, we propose, based on the results of our flow cytometric analyses, that SP Scat Very h'9h CD48Very'°W cells should be highly enriched for HSCs. This would allow for a simple new strategy for the isolation of mouse BM HSCs ex vivo.

FGF21 signalling pathway and metabolic traits - genetic association analysis.

Fibroblast growth factor 21 (FGF21) is a novel master regulator of metabolic profile. The biological actions of FGF21 are elicited upon its klotho beta (KLB)-facilitated binding to FGF receptor 1 (FGFR1), FGFR2 and FGFR3. We hypothesised that common polymorphisms in the FGF21 signalling pathway may be associated with metabolic risk. At the screening stage, we examined associations between 63 common single-nucleotide polymorphisms (SNPs) in five genes of this pathway (FGF21, KLB, FGFR1, FGFR2, FGFR3) and four metabolic phenotypes (LDL cholesterol - LDL-C, HDL-cholesterol - HDL-C, triglycerides and body mass index) in 629 individuals from Silesian Hypertension Study (SHS). Replication analyses were performed in 5478 unrelated individuals of the Swiss CoLaus cohort (imputed genotypes) and in 3030 directly genotyped individuals of the German Myocardial Infarction Family Study (GerMIFS). Of 54 SNPs that met quality control criteria after genotyping in SHS, 4 (rs4733946 and rs7012413 in FGFR1; rs2071616 in FGFR2 and rs7670903 in KLB) showed suggestive association with LDL-C (P=0.0006, P=0.0013, P=0.0055, P=0.011, respectively) and 1 (rs2608819 in KLB) was associated with body mass index (P=0.011); all with false discovery rate q<0.5. Of these, only one FGFR2 polymorphism (rs2071616) showed replicated association with LDL-C in both CoLaus (P=0.009) and men from GerMIFS (P=0.017). The direction of allelic effect of rs2071616 upon LDL-C was consistent in all examined populations. These data show that common genetic variations in FGFR2 may be associated with LDL-C in subjects of white European ancestry.

Sequencing and characterizing the genome of Estrella lausannensis as an undergraduate project: training students and biological insights.

With the widespread availability of high-throughput sequencing technologies, sequencing projects have become pervasive in the molecular life sciences. The huge bulk of data generated daily must be analyzed further by biologists with skills in bioinformatics and by "embedded bioinformaticians," i.e., bioinformaticians integrated in wet lab research groups. Thus, students interested in molecular life sciences must be trained in the main steps of genomics: sequencing, assembly, annotation and analysis. To reach that goal, a practical course has been set up for master students at the University of Lausanne: the "Sequence a genome" class. At the beginning of the academic year, a few bacterial species whose genome is unknown are provided to the students, who sequence and assemble the genome(s) and perform manual annotation. Here, we report the progress of the first class from September 2010 to June 2011 and the results obtained by seven master students who specifically assembled and annotated the genome of Estrella lausannensis, an obligate intracellular bacterium related to Chlamydia. The draft genome of Estrella is composed of 29 scaffolds encompassing 2,819,825 bp that encode for 2233 putative proteins. Estrella also possesses a 9136 bp plasmid that encodes for 14 genes, among which we found an integrase and a toxin/antitoxin module. Like all other members of the Chlamydiales order, Estrella possesses a highly conserved type III secretion system, considered as a key virulence factor. The annotation of the Estrella genome also allowed the characterization of the metabolic abilities of this strictly intracellular bacterium. Altogether, the students provided the scientific community with the Estrella genome sequence and a preliminary understanding of the biology of this recently-discovered bacterial genus, while learning to use cutting-edge technologies for sequencing and to perform bioinformatics analyses.

Clinical utility of g-cimp and idh1 status as dual prognostic markers in glioblastoma

The glioma CpG island methylator phenotype (G-CIMP) has been shown to be highly correlated with prognosis andwas noted to be highly concordant with IDH1mutation in malignant glioma in the limited number of samples analyzed. To better understand the relationship of G-CIMP with IDH1 mutation status and patient outcome, we examined G-CIMP status in detail in a larger retrospective series of glioblastomas as well as tumor samples from the RTOG 0525 clinical trial. Sampleswere tested for 6 CIMPmarkers andwere correlated with patient outcomes. In the retrospective tumor set (n ¼ 301),we found 3 distinct survival groups based on the number of CIMP markers: 0-1 (CIMP-negative), 2-4 (CIMP-intermediate), and 5 or greater (CIMP-positive) with median survivals 13.8, 20.1, and 90.6 months, respectively. This finding was validated in the RTOG 0525 samples (median survivals 15.0, 20.3, and 37.0 months). Among 787 cases with both IDH and CIMP data, 617 were CIMP-negative, 136 were CIMP-intermediate, and 34 were CIMP-positive. Seven hundred forty-four were wild type for IDH1 mutation, and 43 were mutant. CIMP and IDH status were positively correlated but outliers were found. Among the 610 CIMP-negative tumors, there were 7 IDH-mutant tumors, which showed no difference in outcome. Similarly, among the 34 CIMP-positive tumors, there were 21 IDH-mutant cases, which also showed no difference in outcome. However, among the CIMP-intermediate cases, there were 15 IDH-mutant cases with significantly (p ¼ 0.0003) improved outcome (medians not reached vs. 18.5 months, 2 year survival 87% vs. 32%). Multivariate analysis showed that both IDH1 mutation status and CIMP status were independent predictors of outcome. These findings suggest the clinical utility of refining the CIMP status into negative, intermediate, and positive groups and the finding that both IDH1 and CIMPstatus are important molecular markers in GBM.

S-system parameter estimation for noisy metabolic profiles using newton-flow analysis.

Biochemical systems are commonly modelled by systems of ordinary differential equations (ODEs). A particular class of such models called S-systems have recently gained popularity in biochemical system modelling. The parameters of an S-system are usually estimated from time-course profiles. However, finding these estimates is a difficult computational problem. Moreover, although several methods have been recently proposed to solve this problem for ideal profiles, relatively little progress has been reported for noisy profiles. We describe a special feature of a Newton-flow optimisation problem associated with S-system parameter estimation. This enables us to significantly reduce the search space, and also lends itself to parameter estimation for noisy data. We illustrate the applicability of our method by applying it to noisy time-course data synthetically produced from previously published 4- and 30-dimensional S-systems. In addition, we propose an extension of our method that allows the detection of network topologies for small S-systems. We introduce a new method for estimating S-system parameters from time-course profiles. We show that the performance of this method compares favorably with competing methods for ideal profiles, and that it also allows the determination of parameters for noisy profiles.

Proteomic analysis of peach fruit mesocarp softening and chilling injury using difference gel electrophoresis (DIGE)

Background: Peach fruit undergoes a rapid softening process that involves a number of metabolic changes. Storing fruit at low temperatures has been widely used to extend its postharvest life. However, this leads to undesired changes, such as mealiness and browning, which affect the quality of the fruit. In this study, a 2-D DIGE approach was designed to screen for differentially accumulated proteins in peach fruit during normal softening as well as under conditions that led to fruit chilling injury. Results:The analysis allowed us to identify 43 spots -representing about 18% of the total number analyzed- that show statistically significant changes. Thirty-nine of the proteins could be identified by mass spectrometry. Some of the proteins that changed during postharvest had been related to peach fruit ripening and cold stress in the past. However, we identified other proteins that had not been linked to these processes. A graphical display of the relationship between the differentially accumulated proteins was obtained using pairwise average-linkage cluster analysis and principal component analysis. Proteins such as endopolygalacturonase, catalase, NADP-dependent isocitrate dehydrogenase, pectin methylesterase and dehydrins were found to be very important for distinguishing between healthy and chill injured fruit. A categorization of the differentially accumulated proteins was performed using Gene Ontology annotation. The results showed that the 'response to stress', 'cellular homeostasis', 'metabolism of carbohydrates' and 'amino acid metabolism' biological processes were affected the most during the postharvest. Conclusions: Using a comparative proteomic approach with 2-D DIGE allowed us to identify proteins that showed stage-specific changes in their accumulation pattern. Several proteins that are related to response to stress, cellular homeostasis, cellular component organization and carbohydrate metabolism were detected as being differentially accumulated. Finally, a significant proportion of the proteins identified had not been associated with softening, cold storage or chilling injury-altered fruit before; thus, comparative proteomics has proven to be a valuable tool for understanding fruit softening and postharvest.

Contrasting the form and use of reformulation markers

This paper deals with the form and use of reformulation markers in research papers written in English, Spanish and Catalan. Considering the form and frequency of themarkers, English papers tends to prefer simple fixed markers and includes less reformulators than Spanish and Catalan. On the contrary, formal Catalan and Spanish papers include more markers, some of which are complex and allow for some structural variability. As for use, reformulation markers establish dynamic relationships between portions of discourse which can be identified in our corpus with expansion, reduction, and permutation. The analysis of the corpus shows that English authors usually reformulate to add more information to the concept (expansion), whereas Catalan and Spanish authors reduce the contents or the implicatures of the previous formulation more frequently than English.

Short term impact of Tribulus terrestris intake on doping control analysis of endogenous steroids.

Tribulus terrestris is a nutritional supplement highly debated regarding its physiological and actual effects on the organism. The main claimed effect is an increase of testosterone anabolic and androgenic action through the activation of endogenous testosterone production. Even if this biological pathway is not entirely proven, T. terrestris is regularly used by athletes. Recently, the analysis of two female urine samples by GC/C/IRMS (gas chromatography/combustion/isotope-ratio-mass-spectrometry) conclusively revealed the administration of exogenous testosterone or its precursors, even if the testosterone glucuronide/epitestosterone glucuronide (T/E) ratio and steroid marker concentrations were below the cut-off values defined by World Anti-Doping Agency (WADA). To argue against this adverse analytical finding, the athletes recognized having used T. terrestris in their diet. In order to test this hypothesis, two female volunteers ingested 500 mg of T. terrestris, three times a day and for two consecutive days. All spot urines were collected during 48 h after the first intake. The (13)C/(12)C ratio of ketosteroids was determined by GC/C/IRMS, the T/E ratio and DHEA concentrations were measured by GC/MS and LH concentrations by radioimmunoassay. None of these parameters revealed a significant variation or increased above the WADA cut-off limits. Hence, the short-term treatment with T. terrestris showed no impact on the endogenous testosterone metabolism of the two subjects.