Tall cell papillary thyroid carcinoma: new diagnostic criteria and mutations in BRAF and TERT

The tall cell (TC) variant of papillary thyroid carcinoma (PTC) has an unfavorable prognosis. The diagnostic criteria remain inconsistent, and the role of a minor TC component is unclear. Molecular diagnostic markers are not available; however, there are two potential candidates: BRAF V600E and telomerase reverse transcriptase (TERT) promoter mutations. Using a novel approach, we enriched a collective with PTCs that harbored an adverse outcome, which overcame the limited statistical power of most studies. This enabled us to review 125 PTC patients, 57 of which had an adverse outcome. The proportion of TCs that constituted a poor prognosis was assessed. All of the tumors underwent sequencing for TERT promoter and BRAF V600E mutational status and were stained with an antibody to detect the BRAF V600E mutation. A 10% cutoff for TCs was significantly associated with advanced tumor stage and lymph node metastasis. Multivariate analysis showed that TCs above 10% were the only significant factor for overall, tumor-specific, and relapse-free survival. Seven percent of the cases had a TERT promoter mutation, whereas 61% demonstrated a BRAF mutation. The presence of TC was significantly associated with TERT promoter and BRAF mutations. TERT predicted highly significant tumor relapse (P<0.001). PTCs comprised of at least 10% TCs are associated with an adverse clinical outcome and should be reported accordingly. BRAF did not influence patient outcome. Nevertheless, a positive status should encourage the search for TCs. TERT promoter mutations are a strong predictor of tumor relapse, but their role as a surrogate marker for TCs is limited.

Towards person-centered neuroimaging markers for resilience and vulnerability in Bipolar Disorder

Improved clinical care for Bipolar Disorder (BD) relies on the identification of diagnostic markers that can reliably detect disease-related signals in clinically heterogeneous populations. At the very least, diagnostic markers should be able to differentiate patients with BD from healthy individuals and from individuals at familial risk for BD who either remain well or develop other psychopathology, most commonly Major Depressive Disorder (MDD). These issues are particularly pertinent to the development of translational applications of neuroimaging as they represent challenges for which clinical observation alone is insufficient. We therefore applied pattern classification to task-based functional magnetic resonance imaging (fMRI) data of the n-back working memory task, to test their predictive value in differentiating patients with BD (n=30) from healthy individuals (n=30) and from patients' relatives who were either diagnosed with MDD (n=30) or were free of any personal lifetime history of psychopathology (n=30). Diagnostic stability in these groups was confirmed with 4-year prospective follow-up. Task-based activation patterns from the fMRI data were analyzed with Gaussian Process Classifiers (GPC), a machine learning approach to detecting multivariate patterns in neuroimaging datasets. Consistent significant classification results were only obtained using data from the 3-back versus 0-back contrast. Using contrast, patients with BD were correctly classified compared to unrelated healthy individuals with an accuracy of 83.5%, sensitivity of 84.6% and specificity of 92.3%. Classification accuracy, sensitivity and specificity when comparing patients with BD to their relatives with MDD, were respectively 73.1%, 53.9% and 94.5%. Classification accuracy, sensitivity and specificity when comparing patients with BD to their healthy relatives were respectively 81.8%, 72.7% and 90.9%. We show that significant individual classification can be achieved using whole brain pattern analysis of task-based working memory fMRI data. The high accuracy and specificity achieved by all three classifiers suggest that multivariate pattern recognition analyses can aid clinicians in the clinical care of BD in situations of true clinical uncertainty regarding the diagnosis and prognosis.

Right Orbitofrontal Corticolimbic and Left Corticocortical White Matter Connectivity Differentiate Bipolar and Unipolar Depression

Objectives: The absence of pathophysiologically relevant diagnostic markers of bipolar disorder (BD) leads to its frequent misdiagnosis as unipolar depression (UD). We aimed to determine whether whole brain white matter connectivity differentiated BD from UD depression. Methods: We employed a three-way analysis of covariance, covarying for age, to examine whole brain fractional anisotropy (FA), and corresponding longitudinal and radial diffusivity, in currently depressed adults: 15 with BD-type I (mean age 36.3 years, SD 12.0 years), 16 with recurrent UD (mean age 32.3 years, SD 10.0 years), and 24 healthy control adults (HC) (mean age 29.5 years, SD 9.43 years). Depressed groups did not differ in depression severity, age of illness onset, and illness duration. Results: There was a main effect of group in left superior and inferior longitudinal fasciculi (SLF and ILF) (all F >= 9.8; p <= .05, corrected). Whole brain post hoc analyses (all t >= 4.2; p <= .05, corrected) revealed decreased FA in left SLF in BD, versus UD adults in inferior temporal cortex and, versus HC, in primary sensory cortex (associated with increased radial and decreased longitudinal diffusivity, respectively); and decreased FA in left ILF in UD adults versus HC. A main effect of group in right uncinate fasciculus (in orbitofrontal cortex) just failed to meet significance in all participants but was present in women. Post hoc analyses revealed decreased right uncinate fasciculus FA in all and in women, BD versus HC. Conclusions: White matter FA in left occipitotemporal and primary sensory regions supporting visuospatial and sensory processing differentiates BD from UD depression. Abnormally reduced FA in right fronto-temporal regions supporting mood regulation, might underlie. predisposition to depression in BD. These measures might help differentiate pathophysiologic processes of BD versus UD depression.

Carapace dentition patterns, morphometrics and allozyme differentiation amongst two toothed freshwater crab species (Potamonautes warreni and P-unispinus) (Decapoda : Brachyura : Potamonautidae) from river systems in South Africa

The taxonomic relationship between two toothed South African river crabs, Potamonautes warreni and P. unispinus, is unclear. The problem stems from the widespread variation in carapace dentition patterns amongst P. warreni individuals over its biogeographic range, where single toothed individuals may appear similar in carapace morphology to P. unispinus. Ten populations of P. warreni and 18 populations of P. unispinus were collected and the morphometric and genetic differentiation between the two taxa quantified. Patterns of morphometric and genetic variation were examined using multivariate statistics and protein gel electrophoresis, respectively. Principal component analyses of carapace characters showed that the two species are morphologically indistinguishable. However, discriminate functions analyses and additional statistical results corroborate the morphological distinction between the two taxa. Allozyme electrophoresis of 17 protein coding loci, indicated a close genetic similarity between the two species (I = 0.92). A fixed allelic difference at one locus (LT-2) and extensive genetic variability at another locus (PGM-1) indicate that two gene pools are present and that the two taxa are genetically isolated. Intraspecific genetic I values for both species were > 0.97 and indicated no apparent genetic structuring on a micro or macro-geographic scale. The variation in carapace dentition among P. warreni populations possesses no genetic basis and may possibly toe the product of ecogenesis. The value of dentition patterns in the systematics of river crabs is discussed. Dentition patterns among river crab species appear to be conserved and reliable as species specific diagnostic markers, but should ideally be used in combination with other morphological data sets and genetic evidence.

Combining the analyses of introgressive hybridisation and linkage mapping to investigate the genetic architecture of population divergence in the lake whitefish (Coregonus clupeaformis, Mitchill)

Adaptation and reproductive isolation, the engines of biological diversity, are still elusive when discussing the genetic bases of speciation. Namely, the number of genes and magnitude of selection acting positively or negatively on genomic traits implicated in speciation is contentious. Here, we describe the first steps of an ongoing research program aimed at understanding the genetic bases of population divergence and reproductive isolation in the lake whitefish (Coregonus clupeaformis). A preliminary linkage map originating from a hybrid cross between dwarf and normal ecotypes is presented, whereby some of the segregating AFLP markers were found to be conserved among natural populations. Maximum-likelihood was used to estimate hybrid indices from non-diagnostic markers at 998 AFLP loci. This allowed identification of the most likely candidate loci that have been under the influence of selection during the natural hybridisation of whitefish originating from different glacial races. As some of these loci could be identified on the linkage map, the possibility that selection of traits in natural populations may eventually be correlated to specific chromosomal regions was demonstrated. The future prospects and potential of these approaches to elucidate the genetic bases of adaptation and reproductive isolation among sympatric ecotypes of lake whitefish is discussed.

Coactivated platelet-derived growth factor receptor {alpha} and epidermal growth factor receptor are potential therapeutic targets in intimal sarcoma.

Intimal sarcoma (IS) is a rare, malignant, and aggressive tumor that shows a relentless course with a concomitant low survival rate and for which no effective treatment is available. In this study, 21 cases of large arterial blood vessel IS were analyzed by immunohistochemistry and fluorescence in situ hybridization and selectively by karyotyping, array comparative genomic hybridization, sequencing, phospho-kinase antibody arrays, and Western immunoblotting in search for novel diagnostic markers and potential molecular therapeutic targets. Ex vivo immunoassays were applied to test the sensitivity of IS primary tumor cells to the receptor tyrosine kinase (RTK) inhibitors imatinib and dasatinib. We showed that amplification of platelet-derived growth factor receptor α (PDGFRA) is a common finding in IS, which should be considered as a molecular hallmark of this entity. This amplification is consistently associated with PDGFRA activation. Furthermore, the tumors reveal persistent activation of the epidermal growth factor receptor (EGFR), concurrent to PDGFRA activation. Activated PDGFRA and EGFR frequently coexist with amplification and overexpression of the MDM2 oncogene. Ex vivo immunoassays on primary IS cells from one case showed the potency of dasatinib to inhibit PDGFRA and downstream signaling pathways. Our findings provide a rationale for investigating therapies that target PDGFRA, EGFR, or MDM2 in IS. Given the clonal heterogeneity of this tumor type and the potential cross-talk between the PDGFRA and EGFR signaling pathways, targeting multiple RTKs and aberrant downstream effectors might be required to improve the therapeutic outcome for patients with this disease.

Pancreatic stone protein as a novel marker for neonatal sepsis.

PURPOSE: Early-onset sepsis (EOS) is one of the main causes for the admission of newborns to the neonatal intensive care unit. However, traditional infection markers are poor diagnostic markers of EOS. Pancreatic stone protein (PSP) is a promising sepsis marker in adults. The aim of this study was to investigate whether determining PSP improves the diagnosis of EOS in comparison with other infection markers. METHODS: This was a prospective multicentre study involving 137 infants with a gestational age of >34 weeks who were admitted with suspected EOS. PSP, procalcitonin (PCT), soluble human triggering receptor expressed on myeloid cells-1 (sTREM-1), macrophage migration inhibitory factor (MIF) and C-reactive protein (CRP) were measured at admission. Receiver-operating characteristic (ROC) curve analysis was performed. RESULTS: The level of PSP in infected infants was significantly higher than that in uninfected ones (median 11.3 vs. 7.5 ng/ml, respectively; p = 0.001). The ROC area under the curve was 0.69 [95 % confidence interval (CI) 0.59-0.80; p < 0.001] for PSP, 0.77 (95 % CI 0.66-0.87; p < 0.001) for PCT, 0.66 (95 % CI 0.55-0.77; p = 0.006) for CRP, 0.62 (0.51-0.73; p = 0.055) for sTREM-1 and 0.54 (0.41-0.67; p = 0.54) for MIF. PSP independently of PCT predicted EOS (p < 0.001), and the use of both markers concomitantly significantly increased the ability to diagnose EOS. A bioscore combining PSP (>9 ng/ml) and PCT (>2 ng/ml) was the best predictor of EOS (0.83; 95 % CI 0.74-0.93; p < 0.001) and resulted in a negative predictive value of 100 % and a positive predictive value of 71 %. CONCLUSIONS: In this prospective study, the diagnostic performance of PSP and PCT was superior to that of traditional markers and a combination bioscore improved the diagnosis of sepsis. Our findings suggest that PSP is a valuable biomarker in combination with PCT in EOS.

The value of serum procalcitonin level for differentiation of infectious from noninfectious causes of fever after orthopaedic surgery.

BACKGROUND: Early diagnosis of postoperative orthopaedic infections is important in order to rapidly initiate adequate antimicrobial therapy. There are currently no reliable diagnostic markers to differentiate infectious from noninfectious causes of postoperative fever. We investigated the value of the serum procalcitonin level in febrile patients after orthopaedic surgery. METHODS: We prospectively evaluated 103 consecutive patients with new onset of fever within ten days after orthopaedic surgery. Fever episodes were classified by two independent investigators who were blinded to procalcitonin results as infectious or noninfectious origin. White blood-cell count, C-reactive protein level, and procalcitonin level were assessed on days 0, 1, and 3 of the postoperative fever. RESULTS: Infection was diagnosed in forty-five (44%) of 103 patients and involved the respiratory tract (eighteen patients), urinary tract (eighteen), joints (four), surgical site (two), bloodstream (two), and soft tissues (one). Unlike C-reactive protein levels and white blood-cell counts, procalcitonin values were significantly higher in patients with infection compared with patients without infection on the day of fever onset (p = 0.04), day 1 (p = 0.07), and day 3 (p = 0.003). Receiver-operating characteristics demonstrated that procalcitonin had the highest diagnostic accuracy, with a value of 0.62, 0.62, and 0.71 on days 0, 1, and 3, respectively. In a multivariate logistic regression analysis, procalcitonin was a significant predictor for postoperative infection on days 0, 1, and 3 of fever with an odds ratio of 2.3 (95% confidence interval, 1.1 to 4.4), 2.3 (95% confidence interval, 1.1 to 5.2), and 3.3 (95% confidence interval, 1.2 to 9.0), respectively. CONCLUSIONS: Serum procalcitonin is a helpful diagnostic marker supporting clinical and microbiological findings for more reliable differentiation of infectious from noninfectious causes of fever after orthopaedic surgery.

A gene expression signature that distinguishes desmoid tumours from nodular fasciitis.

Nodular fasciitis (NF) is a rapidly growing cellular mass composed of fibroblasts/myofibroblasts, usually localized in subcutaneous tissues, that typically undergoes fibrosis and almost never recurs. Desmoid tumours (DTs) are rare forms of fibroblastic/myofibroblastic growth that arise in deep soft tissues, display a propensity for local infiltration and recurrence, but fail to metastasize. Given that both entities are primarily fibroblastic/myofibroblastic lesions with overlapping histological features, their gene expression profiles were compared to identify differentially expressed genes that may provide not only potential diagnostic markers, but also clues as to the pathogenesis of each disorder. Differentially expressed transcripts (89 clones displaying increased expression in DTs and 246 clones displaying increased expression in NF) included genes encoding several receptor and non-receptor tyrosine kinases (EPHB3, PTPRF, GNAZ, SYK, LYN, EPHA4, BIRC3), transcription factors (TWIST1, PITX2, EYA2, OAS1, MITF, TCF20), and members of the Wnt signalling pathway (AXIN2, WISP1, SFRP). Remarkably, almost one-quarter of the differentially expressed genes encode proteins associated with inflammation and tissue remodelling, including members of the interferon (IFN), tumour necrosis factor (TNF), and transforming growth factor beta (TGF-beta) signalling pathways as well as metalloproteinases (MMP1, 9, 13, 23), urokinase plasminogen activator (PLAU), and cathepsins. The observations provide the first comparative molecular characterization of desmoid tumours and nodular fasciitis and suggest that selected tyrosine kinases, transcription factors, and members of the Wnt, TGF-beta, IFN, and TNF signalling pathways may be implicated in influencing and distinguishing their fate.

Pathogenesis of hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is one of the most common malignant tumours worldwide. The major aetiologies and risk factors for the development of HCC are well defined and some of the multiple steps involved in hepatocarcinogenesis have been elucidated in recent years. However, no clear picture of how and in what sequence these factors interact at the molecular level has emerged yet. Malignant transformation of hepatocytes may occur as a consequence of various aetiologies, such as chronic viral hepatitis, alcohol, and metabolic disorders, in the context of increased cellular turnover induced by chronic liver injury, regeneration and cirrhosis. Activation of cellular oncogenes, inactivation of tumour suppressor genes, genomic instability, including DNA mismatch repair defects and impaired chromosomal segregation, overexpression of growth and angiogenic factors, and telomerase activation may contribute to the development of HCC. Overall, HCCs are genetically very heterogeneous tumours. New technologies, including gene expression profiling and proteomic analyses, should allow us to further elucidate the molecular events underlying HCC development and identify novel diagnostic markers as well as therapeutic targets.

Genetic variability of the Brazilian hair sheep breeds

The objectives of this work were to investigate the genetic structure of the Brazilian hair sheep breeds and to determine the origin of the Santa Inês breed. Molecular similarity was determined using Randomly Amplified Polymorphic DNA - Polymerase Chain Reaction markers in 238 individuals from five naturalized sheep breeds: Santa Inês (48 animals), Rabo Largo (48), Somali (48), Morada Nova (48) and Bergamasca (46), collected in Goiás, Sergipe, Bahia, and Ceará States as well as in the Federal District. Fifty-four loci were selected from 19 primers, after a pilot test using 140 primers. Qualitative analyses indicate diagnostic markers for all breeds. All breeds were significantly different from each other. Interbreed differences were explained by 14.92% of the total variation. Santa Inês clustered with Bergamasca (97% bootstrap) and with Rabo Largo, composing the third member of the group (81% bootstrap) while Morada Nova and Somali breeds clustered separately. Each breed should be considered as a separate management and conservation unit, and special care should be taken with Rabo Largo, Morada Nova and Somali breeds, represented by small herds in Brazil.

Urinary pyridinoline cross-links as biomarkers of osteogenesis imperfecta.

Osteogenesis imperfecta (OI) is a group of genetic heterogeneous connective tissue disorders characterized by increased bone fragility and susceptibility to fractures. Laboratory diagnosis relies on time-consuming and cost-intensive biochemical and molecular genetics analyses. Therefore, it is desirable to identify and establish new diagnostic markers for OI that are reliable, cost-effective and easily accessible. In our study we have identified the ratio of the urinary pyridinoline cross-links lysyl-pyridinoline and hydroxylysyl-pyridinoline as a promising, time- and cost-effective biomarker for osteogenesis imperfecta, that could be used furthermore to investigate cases of suspected non-accidental injury in infants.

Differential gene expression profiles of hepatocellular carcinomas associated or not with viral infection

Chronic hepatitis B (HBV) and C (HCV) virus infections are the most important factors associated with hepatocellular carcinoma (HCC), but tumor prognosis remains poor due to the lack of diagnostic biomarkers. In order to identify novel diagnostic markers and therapeutic targets, the gene expression profile associated with viral and non-viral HCC was assessed in 9 tumor samples by oligo-microarrays. The differentially expressed genes were examined using a z-score and KEGG pathway for the search of ontological biological processes. We selected a non-redundant set of 15 genes with the lowest P value for clustering samples into three groups using the non-supervised algorithm k-means. Fisher’s linear discriminant analysis was then applied in an exhaustive search of trios of genes that could be used to build classifiers for class distinction. Different transcriptional levels of genes were identified in HCC of different etiologies and from different HCC samples. When comparing HBV-HCC vs HCV-HCC, HBV-HCC/HCV-HCC vs non-viral (NV)-HCC, HBC-HCC vs NV-HCC, and HCV-HCC vs NV-HCC of the 58 non-redundant differentially expressed genes, only 6 genes (IKBKβ, CREBBP, WNT10B, PRDX6, ITGAV, and IFNAR1) were found to be associated with hepatic carcinogenesis. By combining trios, classifiers could be generated, which correctly classified 100% of the samples. This expression profiling may provide a useful tool for research into the pathophysiology of HCC. A detailed understanding of how these distinct genes are involved in molecular pathways is of fundamental importance to the development of effective HCC chemoprevention and treatment.

Développement et fonctionnement des mécanismes de résonance motrice chez l'humain

La découverte dans le cerveau du singe macaque de cellules visuo-motrices qui répondent de façon identique à la production et la perception d’actes moteurs soutient l’idée que ces cellules, connues sous le nom de neurones-miroirs, encoderaient la représentation d’actes moteurs. Ces neurones, et le système qu’ils forment, constitueraient un système de compréhension moteur; par delà la simple représentation motrice, il est également possible que ce système participe à des processus de haut niveau en lien avec la cognition sociale. Chez l’humain adulte, des études d’imagerie récentes montrent d’importants chevauchements entre les patrons d’activité liés à l’exécution d’actes moteurs et ceux associés à la perception d’actions. Cependant, malgré le nombre important d’études sur ce système de résonance motrice, étonnamment peu se sont penchées sur les aspects développementaux de ce mécanisme, de même que sa relation avec certaines habiletés sociales dans la population neurotypique. De plus, malgré l’utilisation répandue de certaines techniques neurophysiologiques pour quantifier l’activité de ce système, notamment l’électroencéphalographie et la stimulation magnétique transcrânienne, on ignore en grande partie la spécificité et la convergence de ces mesures dans l’étude des processus de résonance motrice. Les études rassemblées ici visent à combler ces lacunes, c'est-à-dire (1) définir l’existence et les propriétés fonctionnelles du système de résonance motrice chez l’enfant humain, (2) établir le lien entre ce système et certaines habiletés sociales spécifiques et (3) déterminer la validité des outils d’investigation couramment utilisés pour mesurer son activité. Dans l’article 1, l’électroencéphalographie quantitative est utilisée afin de mesurer l’activité des régions sensorimotrices chez un groupe d’enfants d’âge scolaire durant la perception d’actions de la main. On y démontre une modulation de l’activité du rythme mu aux sites centraux non seulement lors de l’exécution de tâches motrices, mais également lors de l’observation passive d’actions. Ces résultats soutiennent l’hypothèse de l’existence d’un système de résonance motrice sensible aux représentations visuelles d’actes moteurs dans le cerveau immature. L’article 2 constitue une étude de cas réalisée chez une jeune fille de 12 ans opérée pour épilepsie réfractaire aux médicaments. L’électroencéphalographie intracrânienne est utilisée afin d’évaluer le recrutement du cortex moteur lors de la perception de sons d’actions. On y montre une modulation de l’activité du cortex moteur, visible dans deux périodes distinctes, qui se reflètent par une diminution de la puissance spectrale des fréquences beta et alpha. Ces résultats soutiennent l’hypothèse de l’existence d’un système de résonance motrice sensible aux représentations auditives d’actions chez l’enfant. L’article 3 constitue une recension des écrits portant sur les données comportementales et neurophysiologiques qui suggèrent la présence d’un système de compréhension d’action fonctionnel dès la naissance. On y propose un modèle théorique où les comportements d’imitation néonataux sont vus comme la résultante de mécanismes d’appariement moteurs non inhibés. Afin de mesurer adéquatement la présence de traits empathiques et autistique dans le but de les mettre en relation avec l’activité du système de résonance motrice, l’article 4 consiste en une validation de versions françaises des échelles Empathy Quotient (Baron-Cohen & Wheelwright, 2004) et Autism Spectrum Quotient (Baron-Cohen et al., 2001) qui seront utilisées dans l’article 5. Les versions traduites de ces échelles ont été administrées à 100 individus sains et 23 personnes avec un trouble du spectre autistique. Les résultats répliquent fidèlement ceux obtenus avec les questionnaires en version anglaise, ce qui suggère la validité des versions françaises. Dans l’article 5, on utilise la stimulation magnétique transcrânienne afin d’investiguer le décours temporel de l’activité du cortex moteur durant la perception d’action et le lien de cette activité avec la présence de traits autistiques et empathiques chez des individus normaux. On y montre que le cortex moteur est rapidement activé suivant la perception d’un mouvement moteur, et que cette activité est corrélée avec les mesures sociocognitives utilisées. Ces résultats suggèrent l’existence d’un système d’appariement moteur rapide dans le cerveau humain dont l’activité est associée aux aptitudes sociales. L’article 6 porte sur la spécificité des outils d’investigation neurophysiologique utilisés dans les études précédentes : la stimulation magnétique transcrânienne et l’électroencéphalographie quantitative. En utilisant ces deux techniques simultanément lors d’observation, d’imagination et d’exécution d’actions, on montre qu’elles évaluent possiblement des processus distincts au sein du système de résonance motrice. En résumé, cette thèse vise à documenter l’existence d’un système de résonance motrice chez l’enfant, d’établir le lien entre son fonctionnement et certaines aptitudes sociales et d’évaluer la validité et la spécificité des outils utilisés pour mesurer l’activité au sein de ce système. Bien que des recherches subséquentes s’avèrent nécessaires afin de compléter le travail entamé ici, les études présentées constituent une avancée significative dans la compréhension du développement et du fonctionnement du système de résonance motrice, et pourraient éventuellement contribuer à l’élaboration d’outils diagnostiques et/ou de thérapeutiques chez des populations où des anomalies de ce système ont été répertoriées.

Chemical characterisation of wild populations of Thymus from different climatic regions in southeast Spain

Thymus is taxonomically a very complex genus with a high frequency of hybridisation and introgression among sympatric species. The variation in accumulation of leaf-surface flavonoids was investigated in 71 wild populations of Thymus front different putative hybrid swarm areas in Andalucia, Spain. Twenty-two flavones, five flavanones, two dihydroflavonols, a flavonol and two unknowns were detected by HPLC-DAD combined with LC-APCI-MS analysis. The majority of compounds were flavones with a lutelin-type substitution of the B-ring, in contrast to previous reports on Macedonian taxa, which predominantly accumulate flavones with apigenin-type substitution of the B-ring. Anatomical and morphometric studies, supported by cluster analysis, identified pure Thymus hyemalis and Thymus baeticus populations, and a large number of putative hybrids. Flavonoid variation was closely related to morphological variation in all populations and is suspected to be a result of genetic polymorphism. Principal component analysis identified the presence of species-specific and geographically linked chemotypes and putative hybrids with mixed morphological and chemical characteristics. Qualitative and quantitative flavonoid accumulation appears to be genetically regulated, while external factors play a secondary role. Flavonoid profiles can thus provide diagnostic markers for the taxonomy of Thymus and are also useful in detecting hybridising taxa. (C) 2007 Elsevier Ltd. All rights reserved.