Clinical pharmacogenomic testing of KRAS, BRAF and EGFR mutations by high resolution melting analysis and ultra-deep pyrosequencing

BACKGROUND: Epidermal growth factor receptor (EGFR) and its downstream factors KRAS and BRAF are mutated in several types of cancer, affecting the clinical response to EGFR inhibitors. Mutations in the EGFR kinase domain predict sensitivity to the tyrosine kinase inhibitors gefitinib and erlotinib in lung adenocarcinoma, while activating point mutations in KRAS and BRAF confer resistance to the anti-EGFR monoclonal antibody cetuximab in colorectal cancer. The development of new generation methods for systematic mutation screening of these genes will allow more appropriate therapeutic choices. METHODS: We describe a high resolution melting (HRM) assay for mutation detection in EGFR exons 19-21, KRAS codon 12/13 and BRAF V600 using formalin-fixed paraffin-embedded samples. Somatic variation of KRAS exon 2 was also analysed by massively parallel pyrosequencing of amplicons with the GS Junior 454 platform. RESULTS: We tested 120 routine diagnostic specimens from patients with colorectal or lung cancer. Mutations in KRAS, BRAF and EGFR were observed in 41.9%, 13.0% and 11.1% of the overall samples, respectively, being mutually exclusive. For KRAS, six types of substitutions were detected (17 G12D, 9 G13D, 7 G12C, 2 G12A, 2 G12V, 2 G12S), while V600E accounted for all the BRAF activating mutations. Regarding EGFR, two cases showed exon 19 deletions (delE746-A750 and delE746-T751insA) and another two substitutions in exon 21 (one showed L858R with the resistance mutation T590M in exon 20, and the other had P848L mutation). Consistent with earlier reports, our results show that KRAS and BRAF mutation frequencies in colorectal cancer were 44.3% and 13.0%, respectively, while EGFR mutations were detected in 11.1% of the lung cancer specimens. Ultra-deep amplicon pyrosequencing successfully validated the HRM results and allowed detection and quantitation of KRAS somatic mutations. CONCLUSIONS: HRM is a rapid and sensitive method for moderate-throughput cost-effective screening of oncogene mutations in clinical samples. Rather than Sanger sequence validation, next-generation sequencing technology results in more accurate quantitative results in somatic variation and can be achieved at a higher throughput scale.

Comprehensive analysis of RET common and rare variant patientsin a series of Spanish Hirschsprung patients confirms a synergistic effect of both kinds of events

Background. RET is the major gene associated to Hirschsprung disease (HSCR) with differential contributions of its rare and common, coding and noncoding mutations to the multifactorial nature of this pathology. In the present study, we have performed a comprehensive study of our HSCR series evaluating the involvement of both RET rare variants (RVs) and common variants (CVs) in the context of the disease. Methods. RET mutational screening was performed by dHPLC and direct sequencing for the identification of RVs. In addition Taqman technology was applied for the genotyping of 3 RET CVs previously associated to HSCR, including a variant lying in an enhancer domain within RET intron 1 (rs2435357). Statistical analyses were performed using the SPSS v.17.0 to analyze the distribution of the variants. Results. Our results confirm the strongest association to HSCR for the "enhancer" variant, and demonstrate a significantly higher impact of it in male versus female patients. Integration of the RET RVs and CVs analysis showed that in 91.66% of cases with both kinds of mutational events, the enhancer allele is in trans with the allele bearing the RET RV. Conclusions. A gender effect exists on both the transmission and distribution of rare coding and common HSCR causing mutations. In addition, these RET CVs and RVs seem to act in a synergistic way leading to HSCR phenotype.

Mutation analysis of genes that control the G1/S cell cycle in melanoma: TP53, CDKN1A, CDKN2A, and CDKN2B.

BACKGROUND The role of genes involved in the control of progression from the G1 to the S phase of the cell cycle in melanoma tumors in not fully known. The aim of our study was to analyse mutations in TP53, CDKN1A, CDKN2A, and CDKN2B genes in melanoma tumors and melanoma cell lines METHODS We analysed 39 primary and metastatic melanomas and 9 melanoma cell lines by single-stranded conformational polymorphism (SSCP). RESULTS The single-stranded technique showed heterozygous defects in the TP53 gene in 8 of 39 (20.5%) melanoma tumors: three new single point mutations in intronic sequences (introns 1 and 2) and exon 10, and three new single nucleotide polymorphisms located in introns 1 and 2 (C to T transition at position 11701 in intron 1; C insertion at position 11818 in intron 2; and C insertion at position 11875 in intron 2). One melanoma tumor exhibited two heterozygous alterations in the CDKN2A exon 1 one of which was novel (stop codon, and missense mutation). No defects were found in the remaining genes. CONCLUSION These results suggest that these genes are involved in melanoma tumorigenesis, although they may be not the major targets. Other suppressor genes that may be informative of the mechanism of tumorigenesis in skin melanomas should be studied.

Comprehensive analysis of NRG1 common and rare variants in Hirschsprung patients.

Hirschsprung disease (HSCR, OMIM 142623) is a developmental disorder characterized by the absence of ganglion cells along variable lengths of the distal gastrointestinal tract, which results in tonic contraction of the aganglionic gut segment and functional intestinal obstruction. The RET proto-oncogene is the major gene for HSCR with differential contributions of its rare and common, coding and noncoding mutations to the multifactorial nature of this pathology. Many other genes have been described to be associated with the pathology, as NRG1 gene (8p12), encoding neuregulin 1, which is implicated in the development of the enteric nervous system (ENS), and seems to contribute by both common and rare variants. Here we present the results of a comprehensive analysis of the NRG1 gene in the context of the disease in a series of 207 Spanish HSCR patients, by both mutational screening of its coding sequence and evaluation of 3 common tag SNPs as low penetrance susceptibility factors, finding some potentially damaging variants which we have functionally characterized. All of them were found to be associated with a significant reduction of the normal NRG1 protein levels. The fact that those mutations analyzed alter NRG1 protein would suggest that they would be related with HSCR disease not only in Chinese but also in a Caucasian population, which reinforces the implication of NRG1 gene in this pathology.

Analysis of genes encoding penicillin-binding proteins in clinical isolates of Acinetobacter baumannii.

There is limited information on the role of penicillin-binding proteins (PBPs) in the resistance of Acinetobacter baumannii to β-lactams. This study presents an analysis of the allelic variations of PBP genes in A. baumannii isolates. Twenty-six A. baumannii clinical isolates (susceptible or resistant to carbapenems) from three teaching hospitals in Spain were included. The antimicrobial susceptibility profile, clonal pattern, and genomic species identification were also evaluated. Based on the six complete genomes of A. baumannii, the PBP genes were identified, and primers were designed for each gene. The nucleotide sequences of the genes identified that encode PBPs and the corresponding amino acid sequences were compared with those of ATCC 17978. Seven PBP genes and one monofunctional transglycosylase (MGT) gene were identified in the six genomes, encoding (i) four high-molecular-mass proteins (two of class A, PBP1a [ponA] and PBP1b [mrcB], and two of class B, PBP2 [pbpA or mrdA] and PBP3 [ftsI]), (ii) three low-molecular-mass proteins (two of type 5, PBP5/6 [dacC] and PBP6b [dacD], and one of type 7 (PBP7/8 [pbpG]), and (iii) a monofunctional enzyme (MtgA [mtgA]). Hot spot mutation regions were observed, although most of the allelic changes found translated into silent mutations. The amino acid consensus sequences corresponding to the PBP genes in the genomes and the clinical isolates were highly conserved. The changes found in amino acid sequences were associated with concrete clonal patterns but were not directly related to susceptibility or resistance to β-lactams. An insertion sequence disrupting the gene encoding PBP6b was identified in an endemic carbapenem-resistant clone in one of the participant hospitals.

Added Value of Deep Sequencing Relative to Population Sequencing in Heavily Pre-Treated HIV-1-Infected Subjects.

OBJECTIVE: To explore the potential of deep HIV-1 sequencing for adding clinically relevant information relative to viral population sequencing in heavily pre-treated HIV-1-infected subjects. METHODS: In a proof-of-concept study, deep sequencing was compared to population sequencing in HIV-1-infected individuals with previous triple-class virological failure who also developed virologic failure to deep salvage therapy including, at least, darunavir, tipranavir, etravirine or raltegravir. Viral susceptibility was inferred before salvage therapy initiation and at virological failure using deep and population sequencing genotypes interpreted with the HIVdb, Rega and ANRS algorithms. The threshold level for mutant detection with deep sequencing was 1%. RESULTS: 7 subjects with previous exposure to a median of 15 antiretrovirals during a median of 13 years were included. Deep salvage therapy included darunavir, tipranavir, etravirine or raltegravir in 4, 2, 2 and 5 subjects, respectively. Self-reported treatment adherence was adequate in 4 and partial in 2; one individual underwent treatment interruption during follow-up. Deep sequencing detected all mutations found by population sequencing and identified additional resistance mutations in all but one individual, predominantly after virological failure to deep salvage therapy. Additional genotypic information led to consistent decreases in predicted susceptibility to etravirine, efavirenz, nucleoside reverse transcriptase inhibitors and indinavir in 2, 1, 2 and 1 subject, respectively. Deep sequencing data did not consistently modify the susceptibility predictions achieved with population sequencing for darunavir, tipranavir or raltegravir. CONCLUSIONS: In this subset of heavily pre-treated individuals, deep sequencing improved the assessment of genotypic resistance to etravirine, but did not consistently provide additional information on darunavir, tipranavir or raltegravir susceptibility. These data may inform the design of future studies addressing the clinical value of minority drug-resistant variants in treatment-experienced subjects.

Prevalence and resistance mutations of non-B HIV-1 subtypes among immigrants in Southern Spain along the decade 2000-2010.

Most of the non-B HIV-1 subtypes are predominant in Sub-Saharan Africa and India although they have been found worldwide. In the last decade, immigration from these areas has increased considerably in Spain. The objective of this study was to evaluate the prevalence of non-B subtypes circulating in a cohort of HIV-1-infected immigrants in Seville, Southern Spain and to identify drug resistance-associated mutations. METHODS: Complete protease and first 220 codons of the reverse transcriptase coding regions were amplified and sequenced by population sequencing. HIV-1 subtypes were determined using Stanford University Drug Resistance Database, and phylogenetic analysis was performed comparing multiple reported sequences. Drug resistance mutations were defined according to the International AIDS Society-USA. RESULTS: From 2000 to 2010 a total of 1,089 newly diagnosed HIV-1-infected patients were enrolled in our cohort. Of these, 121 were immigrants, of which 98 had ethical approval and informed consent to include in our study. Twenty-nine immigrants (29/98, 29.6%) were infected with non-B subtypes, of which 15/29 (51.7%) were CRF02-AG, mostly from Sub-Saharan Africa, and 2/29 (6.9%) were CRF01-AE from Eastern Europe. A, C, F, J and G subtypes from Eastern Europe, Central-South America and Sub-Saharan Africa were also present. Some others harboured recombinant forms CRF02-AG/CRF01-AE, CRF2-AG/G and F/B, B/C, and K/G, in PR and RT-coding regions. Patients infected with non-B subtypes showed a high frequency of minor protease inhibitor resistance mutations, M36I, L63P, and K20R/I. Only one patient, CRF02_AG, showed major resistance mutation L90M. Major RT inhibitor resistance mutations K70R and A98G were present in one patient with subtype G, L100I in one patient with CRF01_AE, and K103N in another patient with CRF01_AE. Three patients had other mutations such as V118I, E138A and V90I. CONCLUSIONS: The circulation of non-B subtypes has significantly increased in Southern Spain during the last decade, with 29.6% prevalence, in association with demographic changes among immigrants. This could be an issue in the treatment and management of these patients. Resistance mutations have been detected in these patients with a prevalence of 7% among treatment-naïve patients compared with the 21% detected among patients under HAART or during treatment interruption.

Functional Variants in NOS1 and NOS2A Are Not Associated with Progressive Hearing Loss in Ménière's Disease in a European Caucasian Population

Hearing loss in Meniere's disease (MD) is associated with loss of spiral ganglion neurons and hair cells. In a guinea pig model of endolymphatic hydrops, nitric oxide synthases (NOS) and oxidative stress mediate loss of spiral ganglion neurons. To test the hypothesis that functional variants of NOS1 and NOS2A are associated with MD, wed genotyped three functional variants of NOS1 (rs41279104,rs2682826, and a cytosine-adenosine microsatellite repeat in exon 1f) and the CCTTT repeat in the promoter of NOS2A gene (rs3833912) in two independent MD sets(273 patients in total) and 550 controls. A third cohort of American patients was genotyped as replication cohort for the CCTTT repeat. Neither allele nor genotype frequencies of rs41279104 and rs2682826 were associated with MD, although longer alleles of the cytosine-adenosine microsatellite repeat were marginally significant (corrected p = 0.05) in the Mediterranean cohort but not in a second Galicia cohort. Shorter numbers of the CCTTT repeat in NOS2A were significantly more frequent in Galicia controls (OR = 0.37 [CI, 0.18-0.76], corrected p =0.04), but this finding could not be replicated in Mediterranean or American case-control populations. Meta-analysis did not support an association between CCTTT repeats and risk for MD. Severe hearing loss (>75 dB) was also not associated with any functional variants studied. Functional variants of NOS1 and and NOS2A do not confer susceptibility for MD.

Toscana virus in Spain

Toscana virus (TOSV, Phlebovirus, family Bunyaviridae) infection is one of the most prevalent arboviruses in Spain. Within the objectives of a multidisciplinary network, a study on the epidemiology of TOSV was conducted in Granada, in southern Spain. The overall seroprevalence rate was 24.9%, significantly increasing with age. TOSV was detected in 3 of 103 sandfly pools by viral culture or reverse transcription-polymerase chain reaction from a region of the L gene. Nucleotide sequence homology was 99%-100% in TOSV from vectors and patients and 80%-81% compared to the Italian strain ISS Phl.3. Sequencing of the N gene of TOSV isolates from patients and vectors indicated 87%-88% and 100% homology at the nucleotide and amino acid levels, respectively, compared to the Italian strain. These findings demonstrate the circulation of at least 2 different lineages of TOSV in the Mediterranean basin, the Italian lineage and the Spanish lineage.

Lack of evidence for KRAS oncogenic mutations in triple-negative breast cancer.

BACKGROUND Mutational analysis of the KRAS gene has recently been established as a complementary in vitro diagnostic tool for the identification of patients with colorectal cancer who will not benefit from anti-epidermal growth factor receptor (EGFR) therapies. Assessment of the mutation status of KRAS might also be of potential relevance in other EGFR-overexpressing tumors, such as those occurring in breast cancer. Although KRAS is mutated in only a minor fraction of breast tumors (5%), about 60% of the basal-like subtype express EGFR and, therefore could be targeted by EGFR inhibitors. We aimed to study the mutation frequency of KRAS in that subtype of breast tumors to provide a molecular basis for the evaluation of anti-EGFR therapies. METHODS Total, genomic DNA was obtained from a group of 35 formalin-fixed paraffin-embedded, triple-negative breast tumor samples. Among these, 77.1% (27/35) were defined as basal-like by immunostaining specific for the established surrogate markers cytokeratin (CK) 5/6 and/or EGFR. KRAS mutational status was determined in the purified DNA samples by Real Time (RT)-PCR using primers specific for the detection of wild-type KRAS or the following seven oncogenic somatic mutations: Gly12Ala, Gly12Asp, Gly12Arg, Gly12Cys, Gly12Ser, Gly12Val and Gly13Asp. RESULTS We found no evidence of KRAS oncogenic mutations in all analyzed tumors. CONCLUSIONS This study indicates that KRAS mutations are very infrequent in triple-negative breast tumors and that EGFR inhibitors may be of potential benefit in the treatment of basal-like breast tumors, which overexpress EGFR in about 60% of all cases.

Mitochondrial haplogroups and polymorphisms reveal no association with sporadic prostate cancer in a southern European population.

BACKGROUND It is known that mitochondria play an important role in certain cancers (prostate, renal, breast, or colorectal) and coronary disease. These organelles play an essential role in apoptosis and the production of reactive oxygen species; in addition, mtDNA also reveals the history of populations and ancient human migration. All these events and variations in the mitochondrial genome are thought to cause some cancers, including prostate cancer, and also help us to group individuals into common origin groups. The aim of the present study is to analyze the different haplogroups and variations in the sequence in the mitochondrial genome of a southern European population consisting of subjects affected (n = 239) and non-affected (n = 150) by sporadic prostate cancer. METHODOLOGY AND PRINCIPAL FINDINGS Using primer extension analysis and DNA sequencing, we identified the nine major European haplogroups and CR polymorphisms. The frequencies of the haplogroups did not differ between patients and control cohorts, whereas the CR polymorphism T16356C was significantly higher in patients with PC compared to the controls (p = 0.029). PSA, staging, and Gleason score were associated with none of the nine major European haplogroups. The CR polymorphisms G16129A (p = 0.007) and T16224C (p = 0.022) were significantly associated with Gleason score, whereas T16311C (p = 0.046) was linked with T-stage. CONCLUSIONS AND SIGNIFICANCE Our results do not suggest that mtDNA haplogroups could be involved in sporadic prostate cancer etiology and pathogenesis as previous studies performed in middle Europe population. Although some significant associations have been obtained in studying CR polymorphisms, further studies should be performed to validate these results.

Unusual manifestation of toscana virus infection, Spain.

This study was supported in part by project 05/305, Junta de Andalucía, Spain

Introduction of an automated system for the diagnosis and quantification of hepatitis B and hepatitis C viruses.

Hepatitis B virus (HBV) and Hepatitis C virus (HCV) infections pose major public health problems because of their prevalence worldwide. Consequently, screening for these infections is an important part of routine laboratory activity. Serological and molecular markers are key elements in diagnosis, prognosis and treatment monitoring for HBV and HCV infections. Today, automated chemiluminescence immunoassay (CLIA) analyzers are widely used for virological diagnosis, particularly in high-volume clinical laboratories. Molecular biology techniques are routinely used to detect and quantify viral genomes as well as to analyze their sequence; in order to determine their genotype and detect resistance to antiviral drugs. Real-time PCR, which provides high sensitivity and a broad dynamic range, has gradually replaced other signal and target amplification technologies for the quantification and detection of nucleic acid. The next-generation DNA sequencing techniques are still restricted to research laboratories.The serological and molecular marker methods available for HBV and HCV are discussed in this article, along with their utility and limitations for use in Chronic Hepatitis B (CHB) diagnosis and monitoring.

Mutational spectrum of the SPG4 (SPAST) and SPG3A (ATL1) genes in Spanish patients with hereditary spastic paraplegia.

BACKGROUND Hereditary Spastic Paraplegias (HSP) are characterized by progressive spasticity and weakness of the lower limbs. At least 45 loci have been identified in families with autosomal dominant (AD), autosomal recessive (AR), or X-linked hereditary patterns. Mutations in the SPAST (SPG4) and ATL1 (SPG3A) genes would account for about 50% of the ADHSP cases. METHODS We defined the SPAST and ATL1 mutational spectrum in a total of 370 unrelated HSP index cases from Spain (83% with a pure phenotype). RESULTS We found 50 SPAST mutations (including two large deletions) in 54 patients and 7 ATL1 mutations in 11 patients. A total of 33 of the SPAST and 3 of the ATL1 were new mutations. A total of 141 (31%) were familial cases, and we found a higher frequency of mutation carriers among these compared to apparently sporadic cases (38% vs. 5%). Five of the SPAST mutations were predicted to affect the pre-mRNA splicing, and in 4 of them we demonstrated this effect at the cDNA level. In addition to large deletions, splicing, frameshifting, and missense mutations, we also found a nucleotide change in the stop codon that would result in a larger ORF. CONCLUSIONS In a large cohort of Spanish patients with spastic paraplegia, SPAST and ATL1 mutations were found in 15% of the cases. These mutations were more frequent in familial cases (compared to sporadic), and were associated with heterogeneous clinical manifestations.

Functional Characterization of a Novel Frameshift Mutation in the C-terminus of the Nav1.5 Channel Underlying a Brugada Syndrome with Variable Expression in a Spanish Family.

INTRODUCTION We functionally analyzed a frameshift mutation in the SCN5A gene encoding cardiac Na(+) channels (Nav1.5) found in a proband with repeated episodes of ventricular fibrillation who presented bradycardia and paroxysmal atrial fibrillation. Seven relatives also carry the mutation and showed a Brugada syndrome with an incomplete and variable expression. The mutation (p.D1816VfsX7) resulted in a severe truncation (201 residues) of the Nav1.5 C-terminus. METHODS AND RESULTS Wild-type (WT) and mutated Nav1.5 channels together with hNavβ1 were expressed in CHO cells and currents were recorded at room temperature using the whole-cell patch-clamp. Expression of p.D1816VfsX7 alone resulted in a marked reduction (≈90%) in peak Na(+) current density compared with WT channels. Peak current density generated by p.D1816VfsX7+WT was ≈50% of that generated by WT channels. p.D1816VfsX7 positively shifted activation and inactivation curves, leading to a significant reduction of the window current. The mutation accelerated current activation and reactivation kinetics and increased the fraction of channels developing slow inactivation with prolonged depolarizations. However, late INa was not modified by the mutation. p.D1816VfsX7 produced a marked reduction of channel trafficking toward the membrane that was not restored by decreasing incubation temperature during cell culture or by incubation with 300 μM mexiletine and 5 mM 4-phenylbutirate. CONCLUSION Despite a severe truncation of the C-terminus, the resulting mutated channels generate currents, albeit with reduced amplitude and altered biophysical properties, confirming the key role of the C-terminal domain in the expression and function of the cardiac Na(+) channel.