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.

Frequency of Fabry disease in male and female haemodialysis patients in Spain

BACKGROUND: Fabry disease (FD), an X-linked lysosomal storage disorder, is caused by a reduced activity of the lysosomal enzyme alpha-galactosidase A. The disorder ultimately leads to organ damage (including renal failure) in males and females. However, heterozygous females usually present a milder phenotype with a later onset and a slower progression. METHODS: A combined enzymatic and genetic strategy was used, measuring the activity of alpha-galactosidase A and genotyping the alpha-galactosidase A gene (GLA) in dried blood samples (DBS) of 911 patients undergoing haemodialysis in centers across Spain. RESULTS: GLA alterations were found in seven unrelated patients (4 males and 3 females). Two novel mutations (p.Gly346AlafsX347 and p.Val199GlyfsX203) were identified as well as a previously described mutation, R118C. The R118C mutation was present in 60% of unrelated patients with GLA causal mutations. The D313Y alteration, considered by some authors as a pseudo-deficiency allele, was also found in two out of seven patients. CONCLUSIONS: Excluding the controversial D313Y alteration, FD presents a frequency of one in 182 individuals (0.55%) within this population of males and females undergoing haemodialysis. Moreover, our findings suggest that a number of patients with unexplained and atypical symptoms of renal disease may have FD. Screening programmes for FD in populations of individuals presenting severe kidney dysfunction, cardiac alterations or cerebrovascular disease may lead to the diagnosis of FD in those patients, the study of their families and eventually the implementation of a specific therapy.

Mutant prevention concentrations of fluoroquinolones for Enterobacteriaceae expressing the plasmid-carried quinolone resistance determinant qnrA1.

The influence of qnrA1 on the development of quinolone resistance in Enterobacteriaceae was evaluated by using the mutant prevention concentration parameter. The expression of qnrA1 considerably increased the mutant prevention concentration compared to strains without this gene. In the presence of qnrA1, mutations in gyrA and parC genes were easily selected to produce high levels of quinolone resistance.

In Vivo Selection of Fluoroquinolone-Resistant Escherichia coli Isolates Expressing Plasmid-Mediated Quinolone Resistance and Expanded-Spectrum β-Lactamase

A ciprofloxacin-resistant Escherichia coli isolate, isolate 1B, was obtained from a urinary specimen of a Canadian patient treated with norfloxacin for infection due to a ciprofloxacin-susceptible isolate, isolate 1A. Both isolates harbored a plasmid-encoded sul1-type integron with qnrA1 and blaVEB-1 genes. Isolate 1B had amino acid substitutions in gyrase and topoisomerase.

Plasmid-mediated quinolone resistance in Australia

The aim of this study was to search for plasmid-encoded quinolone resistance determinants QnrA and QnrS in fluoroquinolone-resistant and extended-spectrum beta-lactamase (ESBL)-producing enterobacterial isolates recovered in Sydney, Australia, in 2002. Twenty-three fluoroquinolone-resistant, of which 16 were also ESBL-positive, enterobacterial and nonrelated isolates were studied. PCR with primers specific for qnrA and qnrS genes and primers specific for a series of ESBL genes were used. A qnrA gene was identified in two ESBL-positive isolates, whereas no qnrS-positive strain was found. The QnrA1 determinant was identified in an Enterobacter cloacae isolate and in a carbapenem-resistant Klebsiella pneumoniae isolate, both of which expressed the same ESBL SHV- 12. Whereas no plasmid was identified in the E. cloacae isolate, K. pneumoniae K149 possessed two conjugative plasmids, one that harbored the qnrA and bla (SHV)-12 genes whereas the other expressed the carbapenemase gene bla (IMP-4). The qnrA gene, was located in both cases downstream of the orf513 recombinase gene and upstream of the qnrA1 gene, a structure identical to that found in sul1-type integron In36 and qnrA-positive strains from Shanghai, China. However, the gene cassettes of the sul1-type integrons were different. This study identified the first plasmid-mediated quinolone resistance determinant in Enterobacteriaceae in Australia.

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.

Multiple sclerosis risk variant HLA-DRB1*1501 associates with high expression of DRB1 gene in different human populations.

The human leukocyte antigen (HLA) DRB1*1501 has been consistently associated with multiple sclerosis (MS) in nearly all populations tested. This points to a specific antigen presentation as the pathogenic mechanism though this does not fully explain the disease association. The identification of expression quantitative trait loci (eQTL) for genes in the HLA locus poses the question of the role of gene expression in MS susceptibility. We analyzed the eQTLs in the HLA region with respect to MS-associated HLA-variants obtained from genome-wide association studies (GWAS). We found that the Tag of DRB1*1501, rs3135388 A allele, correlated with high expression of DRB1, DRB5 and DQB1 genes in a Caucasian population. In quantitative terms, the MS-risk AA genotype carriers of rs3135388 were associated with 15.7-, 5.2- and 8.3-fold higher expression of DQB1, DRB5 and DRB1, respectively, than the non-risk GG carriers. The haplotype analysis of expression-associated variants in a Spanish MS cohort revealed that high expression of DRB1 and DQB1 alone did not contribute to the disease. However, in Caucasian, Asian and African American populations, the DRB1*1501 allele was always highly expressed. In other immune related diseases such as type 1 diabetes, inflammatory bowel disease, ulcerative colitis, asthma and IgA deficiency, the best GWAS-associated HLA SNPs were also eQTLs for different HLA Class II genes. Our data suggest that the DR/DQ expression levels, together with specific structural properties of alleles, seem to be the causal effect in MS and in other immunopathologies rather than specific antigen presentation alone.

Hypermethylated 14-3-3-sigma and ESR1 gene promoters in serum as candidate biomarkers for the diagnosis and treatment efficacy of breast cancer metastasis

Background: Numerous hypermethylated genes have been reported in breast cancer, and the silencing of these genes plays an important role in carcinogenesis, tumor progression and diagnosis. These hypermethylated promoters are very rarely found in normal breast. It has been suggested that aberrant hypermethylation may be useful as a biomarker, with implications for breast cancer etiology, diagnosis, and management. The relationship between primary neoplasm and metastasis remains largely unknown. There has been no comprehensive comparative study on the clinical usefulness of tumor-associated methylated DNA biomarkers in primary breast carcinoma and metastatic breast carcinoma. The objective of the present study was to investigate the association between clinical extension of breast cancer and methylation status of Estrogen Receptor1 (ESR1) and Stratifin (14-3-3-σ) gene promoters in disease-free and metastatic breast cancer patients. Methods: We studied two cohorts of patients: 77 patients treated for breast cancer with no signs of disease, and 34 patients with metastatic breast cancer. DNA was obtained from serum samples, and promoter methylation status was determined by using DNA bisulfite modification and quantitative methylation-specific PCR. Results: Serum levels of methylated gene promoter 14-3-3-σ significantly differed between Control and Metastatic Breast Cancer groups (P < 0.001), and between Disease-Free and Metastatic Breast Cancer groups (P < 0.001). The ratio of the 14-3-3-σ level before the first chemotherapy cycle to the level just before administration of the second chemotherapy cycle was defined as the Biomarker Response Ratio [BRR]. We calculated BRR values for the "continuous decline" and "rise-and-fall" groups. Subsequent ROC analysis showed a sensitivity of 75% (95% CI: 47.6 - 86.7) and a specificity of 66.7% (95% CI: 41.0 - 86.7) to discriminate between the groups for a cut-off level of BRR = 2.39. The area under the ROC curve (Z = 0.804 ± 0.074) indicates that this test is a good approach to post-treatment prognosis. Conclusions: The relationship of 14-3-3-σ with breast cancer metastasis and progression found in this study suggests a possible application of 14-3-3-σ as a biomarker to screen for metastasis and to follow up patients treated for metastatic breast cancer, monitoring their disease status and treatment response.

CD40: Novel Association with Crohn's Disease and Replication in Multiple Sclerosis Susceptibility

Background: A functional polymorphism located at 21 from the start codon of the CD40 gene, rs1883832, was previously reported to disrupt a Kozak sequence essential for translation. It has been consistently associated with Graves’ disease risk in populations of different ethnicity and genetic proxies of this variant evaluated in genome-wide association studies have shown evidence of an effect in rheumatoid arthritis and multiple sclerosis (MS) susceptibility. However, the protective allele associated with Graves’ disease or rheumatoid arthritis has shown a risk role in MS, an effect that we aimed to replicate in the present work. We hypothesized that this functional polymorphism might also show an association with other complex autoimmune condition such as inflammatory bowel disease, given the CD40 overexpression previously observed in Crohn’s disease (CD) lesions. Methodology: Genotyping of rs1883832C.T was performed in 1564 MS, 1102 CD and 969 ulcerative colitis (UC) Spanish patients and in 2948 ethnically matched controls by TaqMan chemistry. Principal Findings: The observed effect of the minor allele rs1883832T was replicated in our independent Spanish MS cohort [p= 0.025; OR (95% CI)= 1.12 (1.01–1.23)]. The frequency of the minor allele was also significantly higher in CD patients than in controls [p= 0.002; OR (95% CI)= 1.19 (1.06–1.33)]. This increased predisposition was not detected in UC patients [p= 0.5; OR (95% CI)= 1.04 (0.93–1.17)]. Conclusion: The impact of CD40 rs1883832 on MS and CD risk points to a common signaling shared by these autoimmune conditions

Distinct mechanisms of loss of IFN-gamma mediated HLA class I inducibility in two melanoma cell lines

BACKGROUND The inability of cancer cells to present antigen on the cell surface via MHC class I molecules is one of the mechanisms by which tumor cells evade anti-tumor immunity. Alterations of Jak-STAT components of interferon (IFN)-mediated signaling can contribute to the mechanism of cell resistance to IFN, leading to lack of MHC class I inducibility. Hence, the identification of IFN-gamma-resistant tumors may have prognostic and/or therapeutic relevance. In the present study, we investigated a mechanism of MHC class I inducibility in response to IFN-gamma treatment in human melanoma cell lines. METHODS Basal and IFN-induced expression of HLA class I antigens was analyzed by means of indirect immunofluorescence flow cytometry, Western Blot, RT-PCR, and quantitative real-time RT-PCR (TaqMan(R) Gene Expression Assays). In demethylation studies cells were cultured with 5-aza-2'-deoxycytidine. Electrophoretic Mobility Shift Assay (EMSA) was used to assay whether IRF-1 promoter binding activity is induced in IFN-gamma-treated cells. RESULTS Altered IFN-gamma mediated HLA-class I induction was observed in two melanoma cells lines (ESTDAB-004 and ESTDAB-159) out of 57 studied, while treatment of these two cell lines with IFN-alpha led to normal induction of HLA class I antigen expression. Examination of STAT-1 in ESTDAB-004 after IFN-gamma treatment demonstrated that the STAT-1 protein was expressed but not phosphorylated. Interestingly, IFN-alpha treatment induced normal STAT-1 phosphorylation and HLA class I expression. In contrast, the absence of response to IFN-gamma in ESTDAB-159 was found to be associated with alterations in downstream components of the IFN-gamma signaling pathway. CONCLUSION We observed two distinct mechanisms of loss of IFN-gamma inducibility of HLA class I antigens in two melanoma cell lines. Our findings suggest that loss of HLA class I induction in ESTDAB-004 cells results from a defect in the earliest steps of the IFN-gamma signaling pathway due to absence of STAT-1 tyrosine-phosphorylation, while absence of IFN-gamma-mediated HLA class I expression in ESTDAB-159 cells is due to epigenetic blocking of IFN-regulatory factor 1 (IRF-1) transactivation.

Development and Preliminary Evaluation of a Rapid Oligochromatographic Assay for Specific Detection of New Human Influenza A H1N1 Virus

A new oligochromatographic assay, Speed-Oligo Novel Influenza A H1N1, was designed and optimized for the specific detection of the 2009 influenza A H1N1 virus. The assay is based on a PCR method coupled to detection of PCR products by means of a dipstick device. The target sequence is a 103-bp fragment within the hemagglutinin gene. The analytical sensitivity of the new assay was measured with serial dilutions of a plasmid that contained the target sequence, and we determined that down to one copy per reaction of the plasmid was reliably detected. Diagnostic performance was assessed with 103 RNAs from suspected cases (40 positive and 63 negative results) previously analyzed with a reference real-time PCR technique. All positive cases were confirmed, and no false-positive results were detected with the new assay. No cross-reactions were observed when other viral strains or clinical samples with other respiratory viruses were tested. According to these results, this new assay has 100% sensitivity and specificity. The turnaround time for the whole procedure was 140 min. The assay may be especially useful for the specific detection of 2009 H1N1 virus in laboratories not equipped with real-time PCR instruments

Systematic survey of clonal complexity in tuberculosis at a populational level and detailed characterization of the isolates involved.

Clonally complex infections by Mycobacterium tuberculosis are progressively more accepted. Studies of their dimension in epidemiological scenarios where the infective pressure is not high are scarce. Our study systematically searched for clonally complex infections (mixed infections by more than one strain and simultaneous presence of clonal variants) by applying mycobacterial interspersed repetitive-unit (MIRU)-variable-number tandem-repeat (VNTR) analysis to M. tuberculosis isolates from two population-based samples of respiratory (703 cases) and respiratory-extrapulmonary (R+E) tuberculosis (TB) cases (71 cases) in a context of moderate TB incidence. Clonally complex infections were found in 11 (1.6%) of the respiratory TB cases and in 10 (14.1%) of those with R+E TB. Among the 21 cases with clonally complex TB, 9 were infected by 2 independent strains and the remaining 12 showed the simultaneous presence of 2 to 3 clonal variants. For the 10 R+E TB cases with clonally complex infections, compartmentalization (different compositions of strains/clonal variants in independent infected sites) was found in 9 of them. All the strains/clonal variants were also genotyped by IS6110-based restriction fragment length polymorphism analysis, which split two MIRU-defined clonal variants, although in general, it showed a lower discriminatory power to identify the clonal heterogeneity revealed by MIRU-VNTR analysis. The comparative analysis of IS6110 insertion sites between coinfecting clonal variants showed differences in the genes coding for a cutinase, a PPE family protein, and two conserved hypothetical proteins. Diagnostic delay, existence of previous TB, risk for overexposure, and clustered/orphan status of the involved strains were analyzed to propose possible explanations for the cases with clonally complex infections. Our study characterizes in detail all the clonally complex infections by M. tuberculosis found in a systematic survey and contributes to the characterization that these phenomena can be found to an extent higher than expected, even in an unselected population-based sample lacking high infective pressure.

Prospective universal application of mycobacterial interspersed repetitive-unit-variable-number tandem-repeat genotyping to characterize Mycobacterium tuberculosis isolates for fast identification of clustered and orphan cases.

The use of molecular tools for genotyping Mycobacterium tuberculosis isolates in epidemiological surveys in order to identify clustered and orphan strains requires faster response times than those offered by the reference method, IS6110 restriction fragment length polymorphism (RFLP) genotyping. A method based on PCR, the mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) genotyping technique, is an option for fast fingerprinting of M. tuberculosis, although precise evaluations of correlation between MIRU-VNTR and RFLP findings in population-based studies in different contexts are required before the methods are switched. In this study, we evaluated MIRU-VNTR genotyping (with a set of 15 loci [MIRU-15]) in parallel to RFLP genotyping in a 39-month universal population-based study in a challenging setting with a high proportion of immigrants. For 81.9% (281/343) of the M. tuberculosis isolates, both RFLP and MIRU-VNTR types were obtained. The percentages of clustered cases were 39.9% (112/281) and 43.1% (121/281) for RFLP and MIRU-15 analyses, and the numbers of clusters identified were 42 and 45, respectively. For 85.4% of the cases, the RFLP and MIRU-15 results were concordant, identifying the same cases as clustered and orphan (kappa, 0.7). However, for the remaining 14.6% of the cases, discrepancies were observed: 16 of the cases clustered by RFLP analysis were identified as orphan by MIRU-15 analysis, and 25 cases identified as orphan by RFLP analysis were clustered by MIRU-15 analysis. When discrepant cases showing subtle genotypic differences were tolerated, the discrepancies fell from 14.6% to 8.6%. Epidemiological links were found for 83.8% of the cases clustered by both RFLP and MIRU-15 analyses, whereas for the cases clustered by RFLP or MIRU-VNTR analysis alone, links were identified for only 30.8% or 38.9% of the cases, respectively. The latter group of cases mainly comprised isolates that could also have been clustered, if subtle genotypic differences had been tolerated. MIRU-15 genotyping seems to be a good alternative to RFLP genotyping for real-time interventional schemes. The correlation between MIRU-15 and IS6110 RFLP findings was reasonable, although some uncertainties as to the assignation of clusters by MIRU-15 analysis were identified.

Characterization of microevolution events in Mycobacterium tuberculosis strains involved in recent transmission clusters.

Under certain circumstances, it is possible to identify clonal variants of Mycobacterium tuberculosis infecting a single patient, probably as a result of subtle genetic rearrangements in part of the bacillary population. We systematically searched for these microevolution events in a different context, namely, recent transmission chains. We studied the clustered cases identified using a population-based universal molecular epidemiology strategy over a 5-year period. Clonal variants of the reference strain defining the cluster were found in 9 (12%) of the 74 clusters identified after the genotyping of 612 M. tuberculosis isolates by IS6110 restriction fragment length polymorphism analysis and mycobacterial interspersed repetitive units-variable-number tandem repeat typing. Clusters with microevolution events were epidemiologically supported and involved 4 to 9 cases diagnosed over a 1- to 5-year period. The IS6110 insertion sites from 16 representative isolates of reference and microevolved variants were mapped by ligation-mediated PCR in order to characterize the genetic background involved in microevolution. Both intragenic and intergenic IS6110 locations resulted from these microevolution events. Among those cases of IS6110 locations in intergenic regions which could have an effect on the regulation of adjacent genes, we identified the overexpression of cytochrome P450 in one microevolved variant using quantitative real-time reverse transcription-PCR. Our results help to define the frequency with which microevolution can be expected in M. tuberculosis transmission chains. They provide a snapshot of the genetic background of these subtle rearrangements and identify an event in which IS6110-mediated microevolution in an isogenic background has functional consequences.

Genetic analysis of caveolin-1 and eNOS genes in colorectal cancer.

Caveolae are involved in physical compartmentalization between different groups of signaling events. Its main component, CAV1, modulates different pathways in cellular physiology. The emerging evidence pointing to the role of CAV1 in cancer led us to study whether different alleles of this gene are associated with colorectal cancer (CRC). Since one of the most characterized enzymes regulated by CAV1 is eNOS, we decided to include both genes in this study. We analyzed five SNPs in 360 unrelated CRC patients and 550 controls from the general population. Two of these SNPs were located within eNOS and three within the CAV1 gene. Although haplotype distribution was not associated with CRC, haplotype TiA (CAV1) was associated with familiar forms of CRC (p<0.05). This was especially evident in CRC antecedents and nuclear forms of CRC. If both CG (eNOS) and TiA (CAV1) haplotypes were taken together, this association increased in significance. Thus, we propose that CAV1, either alone or together with eNOS alleles, might modify CRC heritability.