Evaluation of a PCR-RFLP assay for the identification of dermatophytes in situ

Dermatophytes are the main cause of superficial mycoses. These fungi have the capacity to invade keratinized tissue of humans or animals to produce infections that are generally restricted to the corneocytes of the skin, hair, and nails. Nevertheless, it is common to obtain negative results from fungal cultures of dermatological specimens where direct mycological examination showed fungal elements (30-40%). However, correct identification of the isolated dermatophytes from Tinea is important to choose the appropriate treatment. Therefore, we aim to develop a rapid polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assay based on 28S rDNA that is able to identify dermatophytes species in positive dermatological samples. PCR-RFLP identification of dermatophytes in skin or hair allowed validation of the results obtained in culture. It was also possible to identify the infectious dermatophytes when direct hair/ skin mycological examination showed fungal elements, but negative results were obtained from fungal culture. As a conclusion, PCR methods may provide significant benefits in the rapid diagnosis of Tinea. First, there is an increase in sensitivity of dermatophytes identification when enough material is available. Secondly, identification of the infecting agent can be obtained in 24 hours with PCR-RFLP or sequencing, whereas results from fungal cultures can take 2-3 weeks.

Onychomycosis insensitive to systemic terbinafine and azole treatments reveals non-dermatophyte moulds as infectious agents.

BACKGROUND: Dermatophytes are the main cause of onychomycoses, but various non-dermatophyte filamentous fungi are often isolated from abnormal nails. OBJECTIVE: Our aim was the in situ identification of the fungal infectious agent in 8 cases of onychomycoses which could not be cured after systemic terbinafine and itraconazole treatment. METHODS: Fungal DNA was extracted from nail samples, and infectious fungi were identified by restriction fragment length polymorphism (RFLP) of amplified fungal ribosomal DNA using a previously described PCR/RFLP assay. RESULTS: PCR/RFLP identification of fungi in nails allows the identification of the infectious agent: Fusarium sp., Acremonium sp. and Aspergillus sp. were found as a sole infectious agent in 5, 2 and 1 cases, respectively. CONCLUSIONS: Fusarium spp. and other non-dermatophyte filamentous fungi are especially difficult to cure in onychomycoses utilising standard treatment with terbinafine and itraconazole. PCR fungal identification helps demonstrate the presence of moulds in order to prescribe alternative antifungal treatments.

Decreased PM10 exposure attenuates age-related lung function decline: genetic variants in p53, p21, and CCND1 modify this effect.

BACKGROUND: Decreasing exposure to airborne particulates was previously associated with reduced age-related decline in lung function. However, whether the benefit from improved air quality depends on genetic background is not known. Recent evidence points to the involvement of the genes p53 and p21 and of the cell cycle control gene cyclin D1 (CCND1) in the response of bronchial cells to air pollution. OBJECTIVE: We determined in 4,326 participants of the Swiss Cohort Study on Air Pollution and Lung and Heart Diseases in Adults (SAPALDIA) whether four single-nucleotide polymorphisms in three genes [CCND1 (rs9344 [P242P], rs667515), p53 (rs1042522 [R72P]), and p21 (rs1801270 [S31R])] modified the previously observed attenuation of the decline in the forced expiratory flow between 25% and 75% of the forced vital capacity (FEF(25-75)) associated with improved air quality. METHODS: Subjects of the prospective population-based SAPALDIA cohort were assessed in 1991 and 2002 by spirometry, questionnaires, and biological sample collection for genotyping. We assigned spatially resolved concentrations of particulate matter with aerodynamic diameter < or = 10 microm (PM(10)) to each participant's residential history 12 months before the baseline and follow-up assessments. RESULTS: The effect of diminishing PM(10) exposure on FEF(25-75) decline appeared to be modified by p53 R72P, CCND1 P242P, and CCND1 rs667515. For example, a 10-microg/m(3) decline in average PM(10) exposure over an 11-year period attenuated the average annual decline in FEF(25-75) by 21.33 mL/year (95% confidence interval, 10.57-32.08) among participants homozygous for the CCND1 (P242P) GG genotype, by 13.72 mL/year (5.38-22.06) among GA genotypes, and by 6.00 mL/year (-4.54 to 16.54) among AA genotypes. CONCLUSIONS: Our results suggest that cell cycle control genes may modify the degree to which improved air quality may benefit respiratory function in adults.

Morphological, ecological and genetic aspects associated with endemism in the Fly Orchid group.

The European genus Ophrys (Orchidaceae) is famous for its insect-like floral morphology, an adaptation for a pseudocopulatory pollination strategy involving Hymenoptera males. A large number of endemic Ophrys species have recently been described, especially within the Mediterranean Basin, which is one of the major species diversity hotspots. Subtle morphological variation and specific pollinator dependence are the two main perceptible criteria for describing numerous endemic taxa. However, the degree to which endemics differ genetically remains a challenging question. Additionally, knowledge regarding the factors underlying the emergence of such endemic entities is limited. To achieve new insights regarding speciation processes in Ophrys, we have investigated species boundaries in the Fly Orchid group (Ophrys insectifera sensu lato) by examining morphological, ecological and genetic evidence. Classically, authors have recognized one widespread taxon (O. insectifera) and two endemics (O. aymoninii from France and O. subinsectifera from Spain). Our research has identified clear morphological and ecological factors segregating among these taxa; however, genetic differences were more ambiguous. Insights from cpDNA sequencing and amplified fragment length polymorphisms genotyping indicated a recent diversification in the three extant Fly Orchid species, which may have been further obscured by active migration and admixture across the European continent. Our genetic results still indicate weak but noticeable phylogeographic clustering that partially correlates with the described species. Particularly, we report several isolated haplotypes and genetic clusters in central and southeastern Europe. With regard to the morphological, ecological and genetic aspects, we discuss the endemism status within the Fly Orchid group from evolutionary, taxonomical and conservation perspectives.

Inference about the number of contributors to a DNA mixture: Comparative analyses of a Bayesian network approach and the maximum allele count method.

In the forensic examination of DNA mixtures, the question of how to set the total number of contributors (N) presents a topic of ongoing interest. Part of the discussion gravitates around issues of bias, in particular when assessments of the number of contributors are not made prior to considering the genotypic configuration of potential donors. Further complication may stem from the observation that, in some cases, there may be numbers of contributors that are incompatible with the set of alleles seen in the profile of a mixed crime stain, given the genotype of a potential contributor. In such situations, procedures that take a single and fixed number contributors as their output can lead to inferential impasses. Assessing the number of contributors within a probabilistic framework can help avoiding such complication. Using elements of decision theory, this paper analyses two strategies for inference on the number of contributors. One procedure is deterministic and focuses on the minimum number of contributors required to 'explain' an observed set of alleles. The other procedure is probabilistic using Bayes' theorem and provides a probability distribution for a set of numbers of contributors, based on the set of observed alleles as well as their respective rates of occurrence. The discussion concentrates on mixed stains of varying quality (i.e., different numbers of loci for which genotyping information is available). A so-called qualitative interpretation is pursued since quantitative information such as peak area and height data are not taken into account. The competing procedures are compared using a standard scoring rule that penalizes the degree of divergence between a given agreed value for N, that is the number of contributors, and the actual value taken by N. Using only modest assumptions and a discussion with reference to a casework example, this paper reports on analyses using simulation techniques and graphical models (i.e., Bayesian networks) to point out that setting the number of contributors to a mixed crime stain in probabilistic terms is, for the conditions assumed in this study, preferable to a decision policy that uses categoric assumptions about N.

Seasonal fluctuations of bacterial community diversity in agricultural soil and experimental validation by laboratory disturbance experiments.

Natural fluctuations in soil microbial communities are poorly documented because of the inherent difficulty to perform a simultaneous analysis of the relative abundances of multiple populations over a long time period. Yet, it is important to understand the magnitudes of community composition variability as a function of natural influences (e.g., temperature, plant growth, or rainfall) because this forms the reference or baseline against which external disturbances (e.g., anthropogenic emissions) can be judged. Second, definition of baseline fluctuations in complex microbial communities may help to understand at which point the systems become unbalanced and cannot return to their original composition. In this paper, we examined the seasonal fluctuations in the bacterial community of an agricultural soil used for regular plant crop production by using terminal restriction fragment length polymorphism profiling (T-RFLP) of the amplified 16S ribosomal ribonucleic acid (rRNA) gene diversity. Cluster and statistical analysis of T-RFLP data showed that soil bacterial communities fluctuated very little during the seasons (similarity indices between 0.835 and 0.997) with insignificant variations in 16S rRNA gene richness and diversity indices. Despite overall insignificant fluctuations, between 8 and 30% of all terminal restriction fragments changed their relative intensity in a significant manner among consecutive time samples. To determine the magnitude of community variations induced by external factors, soil samples were subjected to either inoculation with a pure bacterial culture, addition of the herbicide mecoprop, or addition of nutrients. All treatments resulted in statistically measurable changes of T-RFLP profiles of the communities. Addition of nutrients or bacteria plus mecoprop resulted in bacteria composition, which did not return to the original profile within 14 days. We propose that at less than 70% similarity in T-RFLP, the bacterial communities risk to drift apart to inherently different states.

Genetic diversity and host plant preferences revealed by simple sequence repeat and mitochondrial markers in a population of the arbuscular mycorrhizal fungus Glomus intraradices.

Arbuscular mycorrhizal fungi (AMF) are important symbionts of plants that improve plant nutrient acquisition and promote plant diversity. Although within-species genetic differences among AMF have been shown to differentially affect plant growth, very little is actually known about the degree of genetic diversity in AMF populations. This is largely because of difficulties in isolation and cultivation of the fungi in a clean system allowing reliable genotyping to be performed. A population of the arbuscular mycorrhizal fungus Glomus intraradices growing in an in vitro cultivation system was studied using newly developed simple sequence repeat (SSR), nuclear gene intron and mitochondrial ribosomal gene intron markers. The markers revealed a strong differentiation at the nuclear and mitochondrial level among isolates. Genotypes were nonrandomly distributed among four plots showing genetic subdivisions in the field. Meanwhile, identical genotypes were found in geographically distant locations. AMF genotypes showed significant preferences to different host plant species (Glycine max, Helianthus annuus and Allium porrum) used before the fungal in vitro culture establishment. Host plants in a field could provide a heterogeneous environment favouring certain genotypes. Such preferences may partly explain within-population patterns of genetic diversity.

First insights into the transcriptome and development of new genomic tools of a widespread circum-Mediterranean tree species, Pinus halepensis Mill.

Aleppo pine (Pinus halepensis Mill.) is a relevant conifer species for studying adaptive responses to drought and fire regimes in the Mediterranean region. In this study, we performed Illumina next-generation sequencing of two phenotypically divergent Aleppo pine accessions with the aims of (i) characterizing the transcriptome through Illumina RNA-Seq on trees phenotypically divergent for adaptive traits linked to fire adaptation and drought, (ii) performing a functional annotation of the assembled transcriptome, (iii) identifying genes with accelerated evolutionary rates, (iv) studying the expression levels of the annotated genes and (v) developing gene-based markers for population genomic and association genetic studies. The assembled transcriptome consisted of 48,629 contigs and covered about 54.6 Mbp. The comparison of Aleppo pine transcripts to Picea sitchensis protein-coding sequences resulted in the detection of 34,014 SNPs across species, with a Ka /Ks average value of 0.216, suggesting that the majority of the assembled genes are under negative selection. Several genes were differentially expressed across the two pine accessions with contrasted phenotypes, including a glutathione-s-transferase, a cellulose synthase and a cobra-like protein. A large number of new markers (3334 amplifiable SSRs and 28,236 SNPs) have been identified which should facilitate future population genomics and association genetics in this species. A 384-SNP Oligo Pool Assay for genotyping with the Illumina VeraCode technology has been designed which showed an high overall SNP conversion rate (76.6%). Our results showed that Illumina next-generation sequencing is a valuable technology to obtain an extensive overview on whole transcriptomes of nonmodel species with large genomes.

Discovery and refinement of loci associated with lipid levels.

Levels of low-density lipoprotein (LDL) cholesterol, high-density lipoprotein (HDL) cholesterol, triglycerides and total cholesterol are heritable, modifiable risk factors for coronary artery disease. To identify new loci and refine known loci influencing these lipids, we examined 188,577 individuals using genome-wide and custom genotyping arrays. We identify and annotate 157 loci associated with lipid levels at P < 5 × 10(-8), including 62 loci not previously associated with lipid levels in humans. Using dense genotyping in individuals of European, East Asian, South Asian and African ancestry, we narrow association signals in 12 loci. We find that loci associated with blood lipid levels are often associated with cardiovascular and metabolic traits, including coronary artery disease, type 2 diabetes, blood pressure, waist-hip ratio and body mass index. Our results demonstrate the value of using genetic data from individuals of diverse ancestry and provide insights into the biological mechanisms regulating blood lipids to guide future genetic, biological and therapeutic research.

Gene duplication, population genomics, and species-level differentiation within a tropical mountain shrub.

Gene duplication leads to paralogy, which complicates the de novo assembly of genotyping-by-sequencing (GBS) data. The issue of paralogous genes is exacerbated in plants, because they are particularly prone to gene duplication events. Paralogs are normally filtered from GBS data before undertaking population genomics or phylogenetic analyses. However, gene duplication plays an important role in the functional diversification of genes and it can also lead to the formation of postzygotic barriers. Using populations and closely related species of a tropical mountain shrub, we examine 1) the genomic differentiation produced by putative orthologs, and 2) the distribution of recent gene duplication among lineages and geography. We find high differentiation among populations from isolated mountain peaks and species-level differentiation within what is morphologically described as a single species. The inferred distribution of paralogs among populations is congruent with taxonomy and shows that GBS could be used to examine recent gene duplication as a source of genomic differentiation of nonmodel species.

Genotyping faeces reveals facultative kin association on capercaillie's leks

The role that kin selection might play in the evolution of lekking in birds remains controversial. Recent molecular data suggest that males displaying on leks are related. Here we investigated the genetic structure and pattern of relatedness on leks of a declining population of capercaillie (Tetrao urogallus) using microsatellite genetic markers. Since the species is highly sensitive to disturbance, we adopted a non-invasive method by using faecal samples collected in the field. Based on a dataset of 50 males distributed in 6 sub-populations, we found significant genetic structuring among sub-populations, and a significant pattern of isolation by distance among leks. Estimates of relatedness showed that males displaying on the same lek were related, even when controlling for the effects of genetical differentiation among sub-populations. In addition, the frequency distribution of relatedness values indicated that leks contain a mixture of close kin and unrelated individuals (34 and 66%, respectively). This pattern is consistent with the hypothesis that leks often contain kin associations, which might be due to very restricted dispersal of some of the males or to joint dispersal of kin. The results are discussed with respect to their implication for the conservation of endangered populations.

Highly parallel genetic approaches in Mendelian and complex diseases

The recent advance in high-throughput sequencing and genotyping protocols allows rapid investigation of Mendelian and complex diseases on a scale not previously been possible. In my thesis research I took advantage of these modern techniques to study retinitis pigmentosa (RP), a rare inherited disease characterized by progressive loss of photoreceptors and leading to blindness; and hypertension, a common condition affecting 30% of the adult population. Firstly, I compared the performance of different next generation sequencing (NGS) platforms in the sequencing of the RP-linked gene PRPF31. The gene contained a mutation in an intronic repetitive element, which presented difficulties for both classic sequencing methods and NGS. We showed that all NGS platforms are powerful tools to identify rare and common DNA variants, also in case of more complex sequences. Moreover, we evaluated the features of different NGS platforms that are important in re-sequencing projects. The main focus of my thesis was then to investigate the involvement of pre-mRNA splicing factors in autosomal dominant RP (adRP). I screened 5 candidate genes in a large cohort of patients by using long-range PCR as enrichment step, followed by NGS. We tested two different approaches: in one, all target PCRs from all patients were pooled and sequenced as a single DNA library; in the other, PCRs from each patient were separated within the pool by DNA barcodes. The first solution was more cost-effective, while the second one allowed obtaining faster and more accurate results, but overall they both proved to be effective strategies for gene screenings in many samples. We could in fact identify novel missense mutations in the SNRNP200 gene, encoding an essential RNA helicase for splicing catalysis. Interestingly, one of these mutations showed incomplete penetrance in one family with adRP. Thus, we started to study the possible molecular causes underlying phenotypic differences between asymptomatic and affected members of this family. For the study of hypertension, I joined a European consortium to perform genome-wide association studies (GWAS). Thanks to the use of very informative genotyping arrays and of phenotipically well-characterized cohorts, we could identify a novel susceptibility locus for hypertension in the promoter region of the endothelial nitric oxide synthase gene (NOS3). Moreover, we have proven the direct causality of the associated SNP using three different methods: 1) targeted resequencing, 2) luciferase assay, and 3) population study. - Le récent progrès dans le Séquençage à haut Débit et les protocoles de génotypage a permis une plus vaste et rapide étude des maladies mendéliennes et multifactorielles à une échelle encore jamais atteinte. Durant ma thèse de recherche, j'ai utilisé ces nouvelles techniques de séquençage afin d'étudier la retinite pigmentale (RP), une maladie héréditaire rare caractérisée par une perte progressive des photorécepteurs de l'oeil qui entraine la cécité; et l'hypertension, une maladie commune touchant 30% de la population adulte. Tout d'abord, j'ai effectué une comparaison des performances de différentes plateformes de séquençage NGS (Next Generation Sequencing) lors du séquençage de PRPF31, un gène lié à RP. Ce gène contenait une mutation dans un élément répétable intronique, qui présentait des difficultés de séquençage avec la méthode classique et les NGS. Nous avons montré que les plateformes de NGS analysées sont des outils très puissants pour identifier des variations de l'ADN rares ou communes et aussi dans le cas de séquences complexes. De plus, nous avons exploré les caractéristiques des différentes plateformes NGS qui sont importantes dans les projets de re-séquençage. L'objectif principal de ma thèse a été ensuite d'examiner l'effet des facteurs d'épissage de pre-ARNm dans une forme autosomale dominante de RP (adRP). Un screening de 5 gènes candidats issus d'une large cohorte de patients a été effectué en utilisant la long-range PCR comme étape d'enrichissement, suivie par séquençage avec NGS. Nous avons testé deux approches différentes : dans la première, toutes les cibles PCRs de tous les patients ont été regroupées et séquencées comme une bibliothèque d'ADN unique; dans la seconde, les PCRs de chaque patient ont été séparées par code barres d'ADN. La première solution a été la plus économique, tandis que la seconde a permis d'obtenir des résultats plus rapides et précis. Dans l'ensemble, ces deux stratégies se sont démontrées efficaces pour le screening de gènes issus de divers échantillons. Nous avons pu identifier des nouvelles mutations faux-sens dans le gène SNRNP200, une hélicase ayant une fonction essentielle dans l'épissage. Il est intéressant de noter qu'une des ces mutations montre une pénétrance incomplète dans une famille atteinte d'adRP. Ainsi, nous avons commencé une étude sur les causes moléculaires entrainant des différences phénotypiques entre membres affectés et asymptomatiques de cette famille. Lors de l'étude de l'hypertension, j'ai rejoint un consortium européen pour réaliser une étude d'association Pangénomique ou genome-wide association study Grâce à l'utilisation de tableaux de génotypage très informatifs et de cohortes extrêmement bien caractérisées au niveau phénotypique, un nouveau locus lié à l'hypertension a été identifié dans la région promotrice du gène endothélial nitric oxide sinthase (NOS3). Par ailleurs, nous avons prouvé la cause directe du SNP associé au moyen de trois méthodes différentes: i) en reséquençant la cible avec NGS, ii) avec des essais à la luciférase et iii) une étude de population.

HLA Class I alleles associated with HCV polymorphisms and response to antiviral therapy in HCV genotype 1-infected patients

Aims: The adaptive immune response against hepatitis C virus (HCV) is significantly shaped by the host's composition of HLA alleles. Thus, the HLA phenotype is a critical determinant of viral evolution during adaptive immune pressure. Potential associations of HLA class I alleles with polymorphisms of HCV immune escape variants are largely unknown. Methods: Direct sequence analysis of the genes encoding the HCV proteins E2, NS3 and NS5B in a cohort of 159 patients with chronic HCV genotype 1 infection who were treated with pegylated interferon-alfa 2b and ribavirin in a prospective controlled trial for 48 weeks was exhibited. HLA class I genotyping was performed by strand-specific reverse hybridization with the INNO-LiPA line probe assays for HLA-A and HLA-B and by strand-specific PCR-SSP. We analyzed each amino acid position of HCV proteins using an extension of Fisher's exact test for associations with HLA alleles. In addition, associations of specific HLA alleles with inflammatory activity, liver fibrosis, HCV RNA viral load and virologic treatment outcome were investigated. Results: Separate analyses of HCV subtype 1a and 1b isolates revealed substantially different patterns of HLA-restricted polymorphisms between subtypes. Only one polymorphism within NS5B (V2758x) was significantly associated with HLA B*15 in HCV genotype 1b infected patients (adjusted p=0,048). However, a number of HLA class I-restricted polymorphisms within novel putative HCV CD8+ T cell epitopes (genotype 1a: HLA-A*11 GTRTIASPK1086-1094 [NS3], HLA-B*07 WPAPQGARSL1111-1120 [NS3]; genotype 1b: HLA-A*24 HYAPRPCGI488-496 [E2], HLA-B*44 GENETDVLL530-538 [E2], HLA-B*15 RVFTEAMTRY2757-2766 [NS5B]) were observed with high predicted epitope binding scores assessed by the web-based software SYFPEITHI (>21). Most of the identified putative epitopes were overlapping with already otherwise published epitopes, indicating a high immunogenicity of the accordant HCV protein region. In addition, certain HLA class I alleles were associated with inflammatory activity, stage of liver fibrosis, and sustained virologic response to antiviral therapy. Conclusions: HLA class I restricted HCV sequence polymorphisms are rare. HCV polymorphisms identified within putative HCV CD8+ T cell epitopes in the present study differ in their genomic distribution between genotype 1a and 1b isolates, implying divergent adaptation to the host's immune pressure on the HCV subtype level.

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.

Effects of KRAS, BRAF, NRAS, and PIK3CA mutations on the efficacy of cetuximab plus chemotherapy in chemotherapy-refractory metastatic colorectal cancer: a retrospective consortium analysis.

Background Following the discovery that mutant KRAS is associated with resistance to anti-epidermal growth factor receptor (EGFR) antibodies, the tumours of patients with metastatic colorectal cancer are now profiled for seven KRAS mutations before receiving cetuximab or panitumumab. However, most patients with KRAS wild-type tumours still do not respond. We studied the effect of other downstream mutations on the efficacy of cetuximab in, to our knowledge, the largest cohort to date of patients with chemotherapy-refractory metastatic colorectal cancer treated with cetuximab plus chemotherapy in the pre-KRAS selection era. Methods 1022 tumour DNA samples (73 from fresh-frozen and 949 from formalin-fixed, paraffin-embedded tissue) from patients treated with cetuximab between 2001 and 2008 were gathered from 11 centres in seven European countries. 773 primary tumour samples had sufficient quality DNA and were included in mutation frequency analyses; mass spectrometry genotyping of tumour samples for KRAS, BRAF, NRAS, and PIK3CA was done centrally. We analysed objective response, progression-free survival (PFS), and overall survival in molecularly defined subgroups of the 649 chemotherapy-refractory patients treated with cetuximab plus chemotherapy. Findings 40.0% (299/747) of the tumours harboured a KRAS mutation, 14.5% (108/743) harboured a PIK3CA mutation (of which 68.5% [74/108] were located in exon 9 and 20.4% [22/108] in exon 20), 4.7% (36/761) harboured a BRAF mutation, and 2.6% (17/644) harboured an NRAS mutation. KRAS mutants did not derive benefit compared with wild types, with a response rate of 6.7% (17/253) versus 35.8% (126/352; odds ratio [OR] 0.13, 95% CI 0.07-0.22; p<0.0001), a median PFS of 12. weeks versus 24 weeks (hazard ratio [HR] 1 98, 1.66-2.36; p<0.0001), and a median overall survival of 32 weeks versus 50 weeks (1.75, 1.47-2.09; p<0.0001). In KRAS wild types, carriers of BRAF and NRAS mutations had a significantly lower response rate than did BRAF and NRAS wild types, with a response rate of 8.3% (2/24) in carriers of BRAF mutations versus 38.0% in BRAF wild types (124/326; OR 0.15, 95% CI 0.02-0.51; p=0.0012); and 7.7% (1/13) in carriers of NRAS mutations versus 38.1% in NRAS wild types (110/289; OR 0.14, 0.007-0.70; p=0.013). PIK3CA exon 9 mutations had no effect, whereas exon 20 mutations were associated with a worse outcome compared with wild types, with a response rate of 0.0% (0/9) versus 36.8% (121/329; OR 0.00,0.00-0.89; p=0.029), a median PFS of 11.5 weeks versus 24 weeks (HR 2.52, 1.33-4.78; p=0.013), and a median overall survival of 34 weeks versus 51 weeks (3.29, 1.60-6.74; p=0.0057). Multivariate analysis and conditional inference trees confirmed that, if KRAS is not mutated, assessing BRAF, NRAS, and PIK3CA exon 20 mutations (in that order) gives additional information about outcome. Objective response rates in our series were 24.4% in the unselected population, 36.3% in the KRAS wild-type selected population, and 41.2% in the KRAS, BRAF, NRAS, and PIK3CA exon 20 wild-type population. Interpretation While confirming the negative effect of KRAS mutations on outcome after cetuximab, we show that BRAF, NRAS, and PIK3CA,exon 20 mutations are significantly associated with a low response rate. Objective response rates could be improved by additional genotyping of BRAF, NRAS, and PIK3CA exon 20 mutations in a KRAS wild-type population.