Claudio Lottaz

Substantial biases in ultra-short read data sets from high-throughput DNA sequencing

Novel sequencing technologies permit the rapid production of large sequence data sets. These technologies are likely to revolutionize genetics and biomedical research, but a thorough characterization of the ultra-short read output is necessary. We generated and analyzed two Illumina 1G ultra-short read data sets, i.e. 2.8 million 27mer reads from a Beta vulgaris genomic clone and 12.3 million 36mers from the Helicobacter acinonychis genome. We found that error rates range from 0.3% at the beginning of reads to 3.8% at the end of reads. Wrong base calls are frequently preceded by base G. Base substitution error frequencies vary by 10- to 11-fold, with A > C transversion being among the most frequent and C > G transversions among the least frequent substitution errors. Insertions and deletions of single bases occur at very low rates. When simulating re-sequencing we found a 20-fold sequencing coverage to be sufficient to compensate errors by correct reads. The read coverage of the sequenced regions is biased; the highest read density was found in intervals with elevated GC content. High Solexa quality scores are over-optimistic and low scores underestimate the data quality. Our results show different types of biases and ways to detect them. Such biases have implications on the use and interpretation of Solexa data, for de novo sequencing, re-sequencing, the identification of single nucleotide polymorphisms and DNA methylation sites, as well as for transcriptome analysis.

T cells become licensed in the lung to enter the central nervous system

The blood–brain barrier (BBB) and the environment of the central nervous system (CNS) guard the nervous tissue from peripheral immune cells. In the autoimmune disease multiple sclerosis, myelin-reactive T-cell blasts are thought to transgress the BBB and create a pro-inflammatory environment in the CNS, thereby making possible a second autoimmune attack that starts from the leptomeningeal vessels and progresses into the parenchyma. Using a Lewis rat model of experimental autoimmune encephalomyelitis, we show here that contrary to the expectations of this concept, T-cell blasts do not efficiently enter the CNS and are not required to prepare the BBB for immune-cell recruitment. Instead, intravenously transferred T-cell blasts gain the capacity to enter the CNS after residing transiently within the lung tissues. Inside the lung tissues, they move along and within the airways to bronchus-associated lymphoid tissues and lung-draining mediastinal lymph nodes before they enter the blood circulation from where they reach the CNS. Effector T cells transferred directly into the airways showed a similar migratory pattern and retained their full pathogenicity. On their way the T cells fundamentally reprogrammed their gene-expression profile, characterized by downregulation of their activation program and upregulation of cellular locomotion molecules together with chemokine and adhesion receptors. The adhesion receptors include ninjurin 1, which participates in T-cell intravascular crawling on cerebral blood vessels. We detected that the lung constitutes a niche not only for activated T cells but also for resting myelin-reactive memory T cells. After local stimulation in the lung, these cells strongly proliferate and, after assuming migratory properties, enter the CNS and induce paralytic disease. The lung could therefore contribute to the activation of potentially autoaggressive T cells and their transition to a migratory mode as a prerequisite to entering their target tissues and inducing autoimmune disease.

Commonly altered genomic regions in acute myeloid leukemia are enriched for somatic mutations involved in chromatin remodeling and splicing

Acute myeloid leukemia (AML) is characterized by molecular heterogeneity. As commonly altered genomic regions point to candidate genes involved in leukemogenesis, we used microarray-based comparative genomic hybridization and single nucleotide polymorphism profiling data of 391 AML cases to further narrow down genomic regions of interest. Targeted resequencing of 1000 genes located in the critical regions was performed in a representative cohort of 50 AML samples comprising all major cytogenetic subgroups. We identified 120 missense/nonsense mutations as well as 60 insertions/deletions affecting 73 different genes (∼ 3.6 tumor-specific aberrations/AML). While most of the newly identified alterations were nonrecurrent, we observed an enrichment of mutations affecting genes involved in epigenetic regulation including known candidates like TET2, TET1, DNMT3A, and DNMT1, as well as mutations in the histone methyltransferases NSD1, EZH2, and MLL3. Furthermore, we found mutations in the splicing factor SFPQ and in the nonclassic regulators of mRNA processing CTCF and RAD21. These splicing-related mutations affected 10% of AML patients in a mutually exclusive manner. In conclusion, we could identify a large number of alterations in genes involved in aberrant splicing and epigenetic regulation in genomic regions commonly altered in AML, highlighting their important role in the molecular pathogenesis of AML.

Activation of the HIF pathway in childhood ALL, prognostic implications of VEGF.

Hypoxia-inducible factor 1 (HIF-1) controls angiogenesis and glycolysis, two leading characteristics of solid tumor invasion, metastasis, and lethality. Increased angiogenesis is also found in the bone marrow (BM) of leukemias. Less is known in leukemia about the role of HIF-1 and vascular endothelial growth factor (VEGF), the most important proangiogenic target gene of HIF-1. We show by immunohistochemistry that the oxygen-regulated component of HIF-1 (HIF-1α) is overexpressed in clusters of leukemic cells in BM specimens of childhood acute lymphoblastic leukemia (ALL) and absent in biopsies of normal BM. Half the HIF-1α-positive ALL biopsies exhibited VEGF coexpression. Among 96 children with relapsed ALL, diagnostic BM aspirates with high VEGF mRNA levels were associated with a significantly lower probability of event-free survival at 3 years (0.31±0.08 vs 0.65±0.07, P=0.003). Those with poor molecular response to therapy (evaluated by MRD assessment) had 2.2-fold higher VEGF levels than those responding well to chemotherapy (P=0.005). In conclusion, the data demonstrate activation of the HIF pathway in the BM of ALL patients and indicate that the expression of HIF target genes, such as VEGF, play an important role in leukemia progression, therapy response, and outcome.

Integrative normalization and comparative analysis for metabolic fingerprinting by comprehensive two-dimensional gas chromatography-time-of-flight mass spectrometry.

Comprehensive two-dimensional gas chromatography-time-of-flight mass spectrometry (GC x GC-TOF-MS) was applied to the comparative metabolic fingerprinting of a wild-type versus a double mutant strain of Escherichia coli lacking the transhydrogenases UdhA and PntAB. Using peak lists generated with the Leco ChromaTOF software as input, we developed retention time correction and data alignment tools (INCA). The accuracy of peak alignment and detection of 1.1- to 4-fold changes in metabolite concentration was validated by a spike-in experiment with 20 standard compounds. A list of 48 significant features that differentiated the two E. coli strains was obtained with an estimated false discovery rate (FDR) of <0.05. A total of 27 metabolites, mainly from the citrate cycle, were identified. That the signal intensity of the m/z 73 trace of the trimethylsilyl (TMS) group reflected true differences in metabolite abundance was confirmed by quantification of pyruvate, fumarate, malate, succinate, alpha-ketoglutarate, citrate, cis-aconitate, myo-inositol, and glucose-6-phosphate using compound specific fragment ions and stable isotope labeled standards. Relative standard deviations for metabolite extraction and GC x GC-TOF-MS analysis of those analytes ranged from 13.2 to 26.3% for the universal m/z 73 trace and 7.4 to 24.5% for the analyte specific fragment ion trace.

Prognostic value of genetic alterations in children with first bone marrow relapse of childhood B-cell precursor acute lymphoblastic leukemia.

Despite risk-adapted treatment, survival of children with relapse of acute lymphoblastic leukemia (ALL) remains poor compared with that of patients with initial diagnosis of ALL. Leukemia-associated genetic alterations may provide novel prognostic factors to refine present relapse treatment strategies. Therefore, we investigated the clinical relevance of 13 recurrent genetic alterations in 204 children treated uniformly for relapsed B-cell precursor ALL according to the ALL-REZ BFM 2002 protocol. The most common alterations were deletions of CDKN2A/2B, IKZF1, PAX5, ETV6, fusion of ETV6-RUNX1 and deletions and/or mutations of TP53. Multivariate analysis identified IKZF1 deletion and TP53 alteration as independent predictors of inferior outcome (P=0.002 and P=0.001). Next, we investigated how both alterations can improve the established risk stratification in relapsed ALL. Intermediate-risk relapse patients with low minimal residual disease are currently considered to have a good prognosis. In this group, deletion of IKZF1 and alteration of TP53 identify patients with significantly inferior outcome (P<0.001). In high-risk relapse patients, deletion of IKZF1 is strongly predictive of a second relapse after stem cell transplantation (P<0.001). We conclude that IKZF1 and TP53 represent relevant prognostic factors that should be considered in future risk assessment of children with relapsed ALL to indicate treatment intensification or intervention.

Expression of Late Cell Cycle Genes and an Increased Proliferative Capacity Characterize Very Early Relapse of Childhood Acute Lymphoblastic Leukemia

Purpose: In childhood acute lymphoblastic leukemia (ALL), ∼25% of patients suffer from relapse. In recurrent disease, despite intensified therapy, overall cure rates of 40% remain unsatisfactory and survival rates are particularly poor in certain subgroups. The probability of long-term survival after relapse is predicted from well-established prognostic factors (i.e., time and site of relapse, immunophenotype, and minimal residual disease). However, the underlying biological determinants of these prognostic factors remain poorly understood. Experimental Design: Aiming at identifying molecular pathways associated with these clinically well-defined prognostic factors, we did gene expression profiling on 60 prospectively collected samples of first relapse patients enrolled on the relapse trial ALL-REZ BFM 2002 of the Berlin-Frankfurt-Münster study group. Results: We show here that patients with very early relapse of ALL are characterized by a distinctive gene expression pattern. We identified a set of 83 genes differentially expressed in very early relapsed ALL compared with late relapsed disease. The vast majority of genes were up-regulated and many were late cell cycle genes with a function in mitosis. In addition, samples from patients with very early relapse showed a significant increase in the percentage of S and G2-M phase cells and this correlated well with the expression level of cell cycle genes. Conclusions: Very early relapse of ALL is characterized by an increased proliferative capacity of leukemic blasts and up-regulated mitotic genes. The latter suggests that novel drugs, targeting late cell cycle proteins, might be beneficial for these patients that typically face a dismal prognosis.

LMX1B is Essential for the Maintenance of Differentiated Podocytes in Adult Kidneys

Mutations of the LMX1B gene cause nail-patella syndrome, a rare autosomal-dominant disorder affecting the development of the limbs, eyes, brain, and kidneys. The characterization of conventional Lmx1b knockout mice has shown that LMX1B regulates the development of podocyte foot processes and slit diaphragms, but studies using podocyte-specific Lmx1b knockout mice have yielded conflicting results regarding the importance of LMX1B for maintaining podocyte structures. In order to address this question, we generated inducible podocyte-specific Lmx1b knockout mice. One week of Lmx1b inactivation in adult mice resulted in proteinuria with only minimal foot process effacement. Notably, expression levels of slit diaphragm and basement membrane proteins remained stable at this time point, and basement membrane charge properties also did not change, suggesting that alternative mechanisms mediate the development of proteinuria in these mice. Cell biological and biophysical experiments with primary podocytes isolated after 1 week of Lmx1b inactivation indicated dysregulation of actin cytoskeleton organization, and time-resolved DNA microarray analysis identified the genes encoding actin cytoskeleton-associated proteins, including Abra and Arl4c, as putative LMX1B targets. Chromatin immunoprecipitation experiments in conditionally immortalized human podocytes and gel shift assays showed that LMX1B recognizes AT-rich binding sites (FLAT elements) in the promoter regions of ABRA and ARL4C, and knockdown experiments in zebrafish support a model in which LMX1B and ABRA act in a common pathway during pronephros development. Our report establishes the importance of LMX1B in fully differentiated podocytes and argues that LMX1B is essential for the maintenance of an appropriately structured actin cytoskeleton in podocytes.

Lactate-Modulated Induction of THBS-1 Activates Transforming Growth Factor (TGF)-beta2 and Migration of Glioma Cells In Vitro

Background An important phenomenon observed in glioma metabolism is increased aerobic glycolysis in tumor cells, which is generally referred to as the Warburg effect. Transforming growth factor (TGF)-beta2, which we previously showed to be induced by lactic acid, is a key pathophysiological factor in glioblastoma, leading to increased invasion and severe local immunosuppression after proteolytic cleavage from its latency associated peptide. In this study we tested the hypothesis, that lactate regulates TGF-beta2 expression and glioma cell migration via induction of Thrombospondin-1 (THBS-1), a TGF-beta activating protein. Methods Lactate levels were reduced by knockdown of LDH-A using specific small interfering RNA (siRNA) and competitive inhibition of LDH-A by sodium oxamate. Knockdown of THBS-1 was performed using specific siRNA. Western Blot, qRT-PCR, and ELISA were used to investigate expression levels of LDH-A, LDH-B, TGF-beta2 and THBS-1. Migration of cells was examined by Spheroid, Scratch and Boyden Chamber assays. Results Knockdown of LDH-A with subsequent decrease of lactate concentration leads to reduced levels of THBS-1 and TGF-beta2 in glioma cells. Lactate addition increases THBS-1 protein, leading to increased activation of TGF-beta2. Inhibition of THBS-1 reduces TGF-beta2 protein and migration of glioma cells. Addition of synthetic THBS-1 can rescue reduced TGF-beta2 protein levels and glioma cell migration in siLDH-A treated cells. Conclusion We define a regulatory cascade between lactate, THBS-1 and TGF-beta2, leading to enhanced migration of glioma cells. Our results demonstrate a specific interaction between tumor metabolism and migration and provide a better understanding of the mechanisms underlying glioma cell invasion.

ReseqChip: automated integration of multiple local context probe data from the MitoChip array in mitochondrial DNA sequence assembly

BackgroundThe Affymetrix MitoChip v2.0 is an oligonucleotide tiling array for the resequencing of the human mitochondrial (mt) genome. For each of 16,569 nucleotide positions of the mt genome it holds two sets of four 25-mer probes each that match the heavy and the light strand of a reference mt genome and vary only at their central position to interrogate all four possible alleles. In addition, the MitoChip v2.0 carries alternative local context probes to account for known mtDNA variants. These probes have been neglected in most studies due to the lack of software for their automated analysis.ResultsWe provide ReseqChip, a free software that automates the process of resequencing mtDNA using multiple local context probes on the MitoChip v2.0. ReseqChip significantly improves base call rate and sequence accuracy. ReseqChip is available at http://code.open-bio.org/svnweb/index.cgi/bioperl/browse/bioperl-live/trunk/Bio/Microarray/Tools/.ConclusionsReseqChip allows for the automated consolidation of base calls from alternative local mt genome context probes. It thereby improves the accuracy of resequencing, while reducing the number of non-called bases.