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.

A novel study of copy number variations in Hirschsprung disease using the multiple ligation-dependent probe amplification (MLPA) technique.

BACKGROUND Hirschsprung disease (HSCR) is a congenital malformation of the hindgut produced by a disruption in neural crest cell migration during embryonic development. HSCR has a complex genetic etiology and mutations in several genes, mainly the RET proto-oncogene, have been related to the disease. There is a clear predominance of missense/nonsense mutations in these genes whereas copy number variations (CNVs) have been seldom described, probably due to the limitations of conventional techniques usually employed for mutational analysis. METHODS In this study we have aimed to analyze the presence of CNVs in some HSCR genes (RET, EDN3, GDNF and ZFHX1B) using the Multiple Ligation-dependent Probe Amplification (MLPA) approach. RESULTS Two alterations in the MLPA profiles of RET and EDN3 were detected, but a detailed inspection showed that the decrease in the corresponding dosages were due to point mutations affecting the hybridization probes regions. CONCLUSION Our results indicate that CNVs of the gene coding regions analyzed here are not a common molecular cause of Hirschsprung disease. However, further studies are required to determine the presence of CNVs affecting non-coding regulatory regions, as well as other candidate genes.

Extracellular tumor-related mRNA in plasma of lymphoma patients and survival implications.

BACKGROUND We studied anomalous extracellular mRNAs in plasma from patients with diffuse large B-cell lymphoma (DLBCL) and their survival implications. mRNAs studied have been reported in the literature as markers of poor (BCL2, CCND2, MYC) and favorable outcome (LMO2, BCL6, FN1) in tumors. These markers were also analyzed in lymphoma tissues to test possible associations with their presence in plasma. METHODOLOGY/PRINCIPAL FINDINGS mRNA from 42 plasma samples and 12 tumors from patients with DLBCL was analyzed by real-time PCR. Samples post-treatment were studied. The immunohistochemistry of BCL2 and BCL6 was defined. Presence of circulating tumor cells was determined by analyzing the clonality of the immunoglobulin heavy-chain genes by PCR. In DLBCL, MYC mRNA was associated with short overall survival. mRNA targets with unfavorable outcome in tumors were associated with characteristics indicative of poor prognosis, with partial treatment response and with short progression-free survival in patients with complete response. In patients with low IPI score, unfavorable mRNA targets were related to shorter overall survival, partial response, high LDH levels and death. mRNA disappeared in post-treatment samples of patients with complete response, and persisted in those with partial response or death. No associations were found between circulating tumor cells and plasma mRNA. Absence of BCL6 protein in tumors was associated with presence of unfavorable plasma mRNA. CONCLUSIONS/SIGNIFICANCE Through a non-invasive procedure, tumor-derived mRNAs can be obtained in plasma. mRNA detected in plasma did not proceed from circulating tumor cells. In our study, unfavorable targets in plasma were associated with poor prognosis in B-cell lymphomas, mainly MYC mRNA. Moreover, the unfavorable targets in plasma could help us to classify patients with poor outcome within the good prognosis group according to IPI.

The effect of homozygous deletion of the BBOX1 and Fibin genes on carnitine level and acyl carnitine profile.

BACKGROUND: Carnitine is a key molecule in energy metabolism that helps transport activated fatty acids into the mitochondria. Its homeostasis is achieved through oral intake, renal reabsorption and de novo biosynthesis. Unlike dietary intake and renal reabsorption, the importance of de novo biosynthesis pathway in carnitine homeostasis remains unclear, due to lack of animal models and description of a single patient defective in this pathway. CASE PRESENTATION: We identified by array comparative genomic hybridization a 42 months-old girl homozygote for a 221 Kb interstitial deletions at 11p14.2, that overlaps the genes encoding Fibin and butyrobetaine-gamma 2-oxoglutarate dioxygenase 1 (BBOX1), an enzyme essential for the biosynthesis of carnitine de novo. She presented microcephaly, speech delay, growth retardation and minor facial anomalies. The levels of almost all evaluated metabolites were normal. Her serum level of free carnitine was at the lower limit of the reference range, while her acylcarnitine to free carnitine ratio was normal. CONCLUSIONS: We present an individual with a completely defective carnitine de novo biosynthesis. This condition results in mildly decreased free carnitine level, but not in clinical manifestations characteristic of carnitine deficiency disorders, suggesting that dietary carnitine intake and renal reabsorption are sufficient to carnitine homeostasis. Our results also demonstrate that haploinsufficiency of BBOX1 and/or Fibin is not associated with Primrose syndrome as previously suggested.

Incidence of liver damage of uncertain origin in HIV patients not co-infected with HCV/HBV.

BACKGROUND AND AIMS Several studies have reported that a significant number of HIV patients not co-infected with HCV/HBV develop liver damage of uncertain origin (LDUO). The objective of our study was to evaluate the incidence of and risk factors for the development of LDUO in HIV infected patients not co-infected with HCV/HBV. METHODS Prospective longitudinal study that included HIV-infected patients free of previous liver damage and viral hepatitis B or C co-infections. Patients were followed up at 6-monthly intervals. Liver stiffness was measured at each visit. Abnormal liver stiffness (ALS) was defined as a liver stiffness value greater than 7.2 kPa at two consecutive measurements. For patients who developed ALS, a protocol was followed to diagnose the cause of liver damage. Those patients who could not be diagnosed with any specific cause of liver disease were diagnosed as LDUO and liver biopsy was proposed. RESULTS 210 patients matched the inclusion criteria and were included. 198 patients completed the study. After a median (Q1-Q3) follow-up of 18 (IQR 12-26) months, 21 patients (10.6%) developed ALS. Of these, fifteen patients were diagnosed as LDUO. The incidence of LDUO was 7.64 cases/100 patient-years. Histological studies were performed on ten (66.6%) patients and all showed liver steatosis. A higher HOMA-IR value and body mass index were independently associated with the development of LDUO. CONCLUSION We found a high incidence of LDUO in HIV-infected patients associated with metabolic risk factors. The leading cause of LDUO in our study was non-alcoholic fatty liver disease.

Gene expression profiling in breast cancer: understanding the molecular basis of histologic grade to improve prognosis.

BACKGROUND: Histologic grade in breast cancer provides clinically important prognostic information. However, 30%-60% of tumors are classified as histologic grade 2. This grade is associated with an intermediate risk of recurrence and is thus not informative for clinical decision making. We examined whether histologic grade was associated with gene expression profiles of breast cancers and whether such profiles could be used to improve histologic grading. METHODS: We analyzed microarray data from 189 invasive breast carcinomas and from three published gene expression datasets from breast carcinomas. We identified differentially expressed genes in a training set of 64 estrogen receptor (ER)-positive tumor samples by comparing expression profiles between histologic grade 3 tumors and histologic grade 1 tumors and used the expression of these genes to define the gene expression grade index. Data from 597 independent tumors were used to evaluate the association between relapse-free survival and the gene expression grade index in a Kaplan-Meier analysis. All statistical tests were two-sided. RESULTS: We identified 97 genes in our training set that were associated with histologic grade; most of these genes were involved in cell cycle regulation and proliferation. In validation datasets, the gene expression grade index was strongly associated with histologic grade 1 and 3 status; however, among histologic grade 2 tumors, the index spanned the values for histologic grade 1-3 tumors. Among patients with histologic grade 2 tumors, a high gene expression grade index was associated with a higher risk of recurrence than a low gene expression grade index (hazard ratio = 3.61, 95% confidence interval = 2.25 to 5.78; P < .001, log-rank test). CONCLUSIONS: Gene expression grade index appeared to reclassify patients with histologic grade 2 tumors into two groups with high versus low risks of recurrence. This approach may improve the accuracy of tumor grading and thus its prognostic value.

HLA and non-HLA genes in Behçet's disease: a multicentric study in the Spanish population.

INTRODUCTION According to genome wide association (GWA) studies as well as candidate gene approaches, Behçet's disease (BD) is associated with human leukocyte antigen (HLA)-A and HLA-B gene regions. The HLA-B51 has been consistently associated with the disease, but the role of other HLA class I molecules remains controversial. Recently, variants in non-HLA genes have also been associated with BD. The aims of this study were to further investigate the influence of the HLA region in BD and to explore the relationship with non-HLA genes recently described to be associated in other populations. METHODS This study included 304 BD patients and 313 ethnically matched controls. HLA-A and HLA-B low resolution typing was carried out by PCR-SSOP Luminex. Eleven tag single nucleotide polymorphisms (SNPs) located outside of the HLA-region, previously described associated with the disease in GWA studies and having a minor allele frequency in Caucasians greater than 0.15 were genotyped using TaqMan assays. Phenotypic and genotypic frequencies were estimated by direct counting and distributions were compared using the χ(2) test. RESULTS In addition to HLA-B*51, HLA-B*57 was found as a risk factor in BD, whereas, B*35 was found to be protective. Other HLA-A and B specificities were suggestive of association with the disease as risk (A*02 and A*24) or protective factors (A*03 and B*58). Regarding the non-HLA genes, the three SNPs located in IL23R and one of the SNPs in IL10 were found to be significantly associated with susceptibility to BD in our population. CONCLUSION Different HLA specificities are associated with Behçet's disease in addition to B*51. Other non-HLA genes, such as IL23R and IL-10, play a role in the susceptibility to the disease.

TRAF1/C5 but not PTPRC variants are potential predictors of rheumatoid arthritis response to anti-tumor necrosis factor therapy.

BACKGROUND The aim of our work was to replicate, in a Southern European population, the association reported in Northern populations between PTPRC locus and response to anti-tumor necrosis factor (anti-TNF) treatment in rheumatoid arthritis (RA). We also looked at associations between five RA risk alleles and treatment response. METHODS We evaluated associations between anti-TNF treatment responses assessed by DAS28 change and by EULAR response at six months in 383 Portuguese patients. Univariate and multivariate linear and logistic regression analyses were performed. In a second step to confirm our findings, we pooled our population with 265 Spanish patients. RESULTS No association was found between PTPRC rs10919563 allele and anti-TNF treatment response, neither in Portuguese modeling for several clinical variables nor in the overall population combining Portuguese and Spanish patients. The minor allele for RA susceptibility, rs3761847 SNP in TRAF1/C5 region, was associated with a poor response in linear and logistic univariate and multivariate regression analyses. No association was observed with the other allellic variants. Results were confirmed in the pooled analysis. CONCLUSION This study did not replicate the association between PTPRC and the response to anti-TNF treatment in our Southern European population. We found that TRAF1/C5 risk RA variants potentially influence anti-TNF treatment response.

Evolution of the correlation between expression divergence and protein divergence in mammals.

Divergence of protein sequences and gene expression patterns are two fundamental mechanisms that generate organismal diversity. Here, we have used genome and transcriptome data from eight mammals and one bird to study the positive correlation of these two processes throughout mammalian evolution. We demonstrate that the correlation is stable over time and most pronounced in neural tissues, which indicates that it is the result of strong negative selection. The correlation is not driven by genes with specific functions and may instead best be viewed as an evolutionary default state, which can nevertheless be evaded by certain gene types. In particular, genes with developmental and neural functions are skewed toward changes in gene expression, consistent with selection against pleiotropic effects associated with changes in protein sequences. Surprisingly, we find that the correlation between expression divergence and protein divergence is not explained by between-gene variation in expression level, tissue specificity, protein connectivity, or other investigated gene characteristics, suggesting that it arises independently of these gene traits. The selective constraints on protein sequences and gene expression patterns also fluctuate in a coordinate manner across phylogenetic branches: We find that gene-specific changes in the rate of protein evolution in a specific mammalian lineage tend to be accompanied by similar changes in the rate of expression evolution. Taken together, our findings highlight many new aspects of the correlation between protein divergence and expression divergence, and attest to its role as a fundamental property of mammalian genome evolution.

Fas ligand elicits a caspase-independent proinflammatory response in human keratinocytes: implications for dermatitis.

Fas ligand (FasL) causes apoptosis of epidermal keratinocytes and triggers the appearance of spongiosis in eczematous dermatitis. We demonstrate here that FasL also aggravates inflammation by triggering the expression of proinflammatory cytokines, chemokines, and adhesion molecules in keratinocytes. In HaCaT cells and in reconstructed human epidermis (RHE), FasL triggered a NF-kappaB-dependent mRNA accumulation of inflammatory cytokines (tumor necrosis factor-alpha, IL-6, and IL-1beta), chemokines (CCL2/MCP-1, CXCL1/GROalpha, CXCL3/GROgamma, and CXCL8/IL-8), and the adhesion molecule ICAM-1. Oligomerization of Fas was required both for apoptosis and for gene expression. Inhibition of caspase activity abolished FasL-dependent apoptosis; however, it failed to suppress the expression of FasL-induced genes. Additionally, in the presence of caspase inhibitors, but not in their absence, FasL triggered the accumulation of CCL5/RANTES (regulated on activation normal T cell expressed and secreted) mRNA. Our findings identify a novel proinflammatory role of FasL in keratinocytes that is independent of caspase activity and is separable from apoptosis. Thus, in addition to causing spongiosis, FasL may play a direct role in triggering and/or sustaining inflammation in eczemas.

Infrequent expression of the MAGE gene family in uveal melanomas.

It has previously been reported that MAGE-1, -2, -3 and -4 genes are expressed in human cancers including cutaneous melanoma. MAGE-1 and MAGE-3 represent targets for specific immunotherapy because they encode peptide antigens which are recognised by cytotoxic T lymphocytes (CTL) when presented by HLA class I molecules, and pilot clinical trials with these peptides are currently in progress. It is likely that other members of the MAGE gene family may also encode antigens recognised by CTL. Uveal melanomas, like cutaneous melanomas, arise from melanocytes that are derived from the neural crest. To determine if uveal melanoma patients would be suitable for MAGE-peptide immunotherapy, the expression of MAGE-1, -2, -3 and -4 genes was assessed by reverse transcription followed by polymerase chain reaction (RT-PCR) amplification and ethidium bromide staining. Expression of MAGE genes was not detected in any of 27 primary tumours. Either MAGE-1 or MAGE-4 was expressed in only 2 of 26 metastatic samples, but expression of MAGE-2 or -3 was not detected. Our data suggest that, unlike cutaneous melanomas, uveal melanomas may not be suitable candidates for MAGE-peptide immunotherapy.

Network-guided analysis of genes with altered somatic copy number and gene expression reveals pathways commonly perturbed in metastatic melanoma.

Cancer genomes frequently contain somatic copy number alterations (SCNA) that can significantly perturb the expression level of affected genes and thus disrupt pathways controlling normal growth. In melanoma, many studies have focussed on the copy number and gene expression levels of the BRAF, PTEN and MITF genes, but little has been done to identify new genes using these parameters at the genome-wide scale. Using karyotyping, SNP and CGH arrays, and RNA-seq, we have identified SCNA affecting gene expression ('SCNA-genes') in seven human metastatic melanoma cell lines. We showed that the combination of these techniques is useful to identify candidate genes potentially involved in tumorigenesis. Since few of these alterations were recurrent across our samples, we used a protein network-guided approach to determine whether any pathways were enriched in SCNA-genes in one or more samples. From this unbiased genome-wide analysis, we identified 28 significantly enriched pathway modules. Comparison with two large, independent melanoma SCNA datasets showed less than 10% overlap at the individual gene level, but network-guided analysis revealed 66% shared pathways, including all but three of the pathways identified in our data. Frequently altered pathways included WNT, cadherin signalling, angiogenesis and melanogenesis. Additionally, our results emphasize the potential of the EPHA3 and FRS2 gene products, involved in angiogenesis and migration, as possible therapeutic targets in melanoma. Our study demonstrates the utility of network-guided approaches, for both large and small datasets, to identify pathways recurrently perturbed in cancer.

An assessment of gene prediction accuracy in large DNA sequences

One of the first useful products from the human genome will be a set of predicted genes. Besides its intrinsic scientific interest, the accuracy and completeness of this data set is of considerable importance for human health and medicine. Though progress has been made on computational gene identification in terms of both methods and accuracy evaluation measures, most of the sequence sets in which the programs are tested are short genomic sequences, and there is concern that these accuracy measures may not extrapolate well to larger, more challenging data sets. Given the absence of experimentally verified large genomic data sets, we constructed a semiartificial test set comprising a number of short single-gene genomic sequences with randomly generated intergenic regions. This test set, which should still present an easier problem than real human genomic sequence, mimics the approximately 200kb long BACs being sequenced. In our experiments with these longer genomic sequences, the accuracy of GENSCAN, one of the most accurate ab initio gene prediction programs, dropped significantly, although its sensitivity remained high. Conversely, the accuracy of similarity-based programs, such as GENEWISE, PROCRUSTES, and BLASTX was not affected significantly by the presence of random intergenic sequence, but depended on the strength of the similarity to the protein homolog. As expected, the accuracy dropped if the models were built using more distant homologs, and we were able to quantitatively estimate this decline. However, the specificities of these techniques are still rather good even when the similarity is weak, which is a desirable characteristic for driving expensive follow-up experiments. Our experiments suggest that though gene prediction will improve with every new protein that is discovered and through improvements in the current set of tools, we still have a long way to go before we can decipher the precise exonic structure of every gene in the human genome using purely computational methodology.

Comparative gene prediction in human and mouse

The completion of the sequencing of the mouse genome promises to help predict human genes with greater accuracy. While current ab initio gene prediction programs are remarkably sensitive (i.e., they predict at least a fragment of most genes), their specificity is often low, predicting a large number of false-positive genes in the human genome. Sequence conservation at the protein level with the mouse genome can help eliminate some of those false positives. Here we describe SGP2, a gene prediction program that combines ab initio gene prediction with TBLASTX searches between two genome sequences to provide both sensitive and specific gene predictions. The accuracy of SGP2 when used to predict genes by comparing the human and mouse genomes is assessed on a number of data sets, including single-gene data sets, the highly curated human chromosome 22 predictions, and entire genome predictions from ENSEMBL. Results indicate that SGP2 outperforms purely ab initio gene prediction methods. Results also indicate that SGP2 works about as well with 3x shotgun data as it does with fully assembled genomes. SGP2 provides a high enough specificity that its predictions can be experimentally verified at a reasonable cost. SGP2 was used to generate a complete set of gene predictions on both the human and mouse by comparing the genomes of these two species. Our results suggest that another few thousand human and mouse genes currently not in ENSEMBL are worth verifying experimentally.

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.