Biochemical characterization of PICK1, a PKC binding protein

Protein kinase C (PKC) is a family of serine-threonine kinases that are activated by a wide variety of hormones, neurotransmitters and growth factors. A single cell type contains multiple isoforms that are translocated to distinct and different subcellular sites upon mitogenic stimulus. Many different cellular responses are attributed to PKC activity though relatively few substrates or binding proteins have been definitively characterized. We used the hinge and catalytic domain of PKC$\alpha$ (PKC7) in a yeast two-hybrid screen to clone proteins that interact with C-kinase (PICKs). One protein which we have termed PICK1 may be involved in PKC$\alpha$-specific function at the level of the nuclear membrane after activation. Binding of PICK1 to PKC$\alpha$ has been shown to be isoform specific as it does not bind to PKC$\beta$II or PKC$\alpha$ in the yeast two-hybrid system. PICK1 mRNA expression level is highest in testis and brain with lower levels of expression in skeletal muscle, heart, kidney, lung and liver. PICK1 protein contains five PKC consensus phosphorylation sites and serves as an in vitro substrate for PKC. The PICK1 protein also contains a P-Loop motif that has been shown to bind ATP or GTP in the Ras family of oncoproteins as well as the G-Protein family. Proteins which bind ATP or GTP using this motif all have some sort of catalytic function although none has been identified for PICK1 as yet. PICK1 contains a DHR/GLGF motif at the N-terminus of the protein. The DHR/GLGF motif is contained in a number of recently described proteins and has been shown to mediate protein-protein interactions at the level of membranes and cytoskeleton. When both PKC$\alpha$ and PICK1 are co-expressed in Cos1 cells the two proteins co-localize to the perinucleus in immunoflouresence studies and co-immunoprecipitate. The binding site for PKC7 has been localized to amino acids 1-358 on PICK1 which contains the DHR/GLGF motif. Binding of PICK1 to PKC$\alpha$ requires the hinge and C-terminal domains of PKC$\alpha$. In vitro, PICK1 binds to PKC$\alpha$ and inhibits its activity as assayed by myelin basic protein phosphorylation. PICK1 also binds to TIS21, a primary response gene that is expressed in response to phorbol ester and growth factor treatment. The Caenorhabditis elegans homologue of PICK1 has been cloned and sequenced revealing a high degree of conservation in the DHR/GLGF motif. A more C-terminal region also shows a high degree of conservation, and the C. elegans PICK1 homologue binds to PKC7 suggesting a conservation of function. Taken together these results suggest that PICK1 may be involved in a PKC$\alpha$-specific function at the level of the nuclear membrane. ^

Functional studies of the homeobox gene {\it Goosecoid\/} during mouse embryogenesis

A fundamental question in developmental biology is to understand the mechanisms that govern the development of an adult individual from a single cell. Goosecoid (Gsc) is an evolutionarily conserved homeobox gene that has been cloned in vertebrates and in Drosophila. In mice, Gsc is first expressed during gastrulation stages where it marks anterior structures of the embryo, this pattern of expression is conserved among vertebrates. Later, expression is observed during organogenesis of the head, limbs and the trunk. The conserved pattern of expression of Gsc during gastrulation and gain of function experiments in Xenopus suggested a function for Gsc in the development of anterior structures in vertebrates. Also, its expression pattern in mouse suggested a role in morphogenesis of the head, limbs and trunk. To determine the functional requirement of Gsc in mice a loss of function mutation was generated by homologous recombination in embryonic stem cells and mice mutant for Gsc were generated.^ Gsc-null mice survived to birth but died hours after delivery. Phenotypic analysis revealed craniofacial and rib cage abnormalities that correlated with the second phase of Gsc expression in the head and trunk but no anomalies were found that correlated with its pattern of expression during gastrulation or limb development.^ To determine the mode of action of Gsc during craniofacial development aggregation chimeras were generated between Gsc-null and wild-type embryos. Chimeras were generated by the aggregation of cleavage stage embryos, taking advantage of two different Gsc-null alleles generated during gene targeting. Chimeras demonstrated a cell-autonomous function for Gsc during craniofacial development and a requirement for Gsc function in cartilage and mesenchymal tissues.^ Thus, during embryogenesis in mice, Gsc is not an essential component of gastrulation as had been suggested in previous experiments. Gsc is required for craniofacial development where it acts cell autonomously in cartilage and mesenchymal tissues. Gsc is also required for proper development of the rib cage but it is dispensable for limb development in mice. ^

Quantifying the turnover of transcriptional subclasses of HIV-1-infected cells

HIV-1-infected cells in peripheral blood can be grouped into different transcriptional subclasses. Quantifying the turnover of these cellular subclasses can provide important insights into the viral life cycle and the generation and maintenance of latently infected cells. We used previously published data from five patients chronically infected with HIV-1 that initiated combination antiretroviral therapy (cART). Patient-matched PCR for unspliced and multiply spliced viral RNAs combined with limiting dilution analysis provided measurements of transcriptional profiles at the single cell level. Furthermore, measurement of intracellular transcripts and extracellular virion-enclosed HIV-1 RNA allowed us to distinguish productive from non-productive cells. We developed a mathematical model describing the dynamics of plasma virus and the transcriptional subclasses of HIV-1-infected cells. Fitting the model to the data allowed us to better understand the phenotype of different transcriptional subclasses and their contribution to the overall turnover of HIV-1 before and during cART. The average number of virus-producing cells in peripheral blood is small during chronic infection. We find that a substantial fraction of cells can become defectively infected. Assuming that the infection is homogenous throughout the body, we estimate an average in vivo viral burst size on the order of 104 virions per cell. Our study provides novel quantitative insights into the turnover and development of different subclasses of HIV-1-infected cells, and indicates that cells containing solely unspliced viral RNA are a good marker for viral latency. The model illustrates how the pool of latently infected cells becomes rapidly established during the first months of acute infection and continues to increase slowly during the first years of chronic infection. Having a detailed understanding of this process will be useful for the evaluation of viral eradication strategies that aim to deplete the latent reservoir of HIV-1.

Aggravation of cardiac myofibroblast arrhythmogeneicity by mechanical stress

Aims Myofibroblasts (MFBs) as appearing in the myocardium during fibrotic remodelling induce slow conduction following heterocellular gap junctional coupling with cardiomyocytes (CMCs) in bioengineered tissue preparations kept under isometric conditions. In this study, we investigated the hypothesis that strain as developed during diastolic filling of the heart chambers may modulate MFB-dependent slow conduction. Methods and results Effects of defined levels of strain on single-cell electrophysiology (patch clamp) and impulse conduction in patterned growth cell strands (optical mapping) were investigated in neonatal rat ventricular cell cultures (Wistar) grown on flexible substrates. While 10.5% strain only minimally affected conduction times in control CMC strands (+3.2%, n.s.), it caused a significant slowing of conduction in the fibrosis model consisting of CMC strands coated with MFBs (conduction times +26.3%). Increased sensitivity to strain of the fibrosis model was due to activation of mechanosensitive channels (MSCs) in both CMCs and MFBs that aggravated the MFB-dependent baseline depolarization of CMCs. As found in non-strained preparations, baseline depolarization of CMCs was partly due to the presence of constitutively active MSCs in coupled MFBs. Constitutive activity of MSCs was not dependent on the contractile state of MFBs, because neither stimulation (thrombin) nor suppression (blebbistatin) thereof significantly affected conduction velocities in the non-strained fibrosis model. Conclusions The findings demonstrate that both constitutive and strain-induced activity of MSCs in MFBs significantly enhance their depolarizing effect on electrotonically coupled CMCs. Ensuing aggravation of slow conduction may contribute to the precipitation of strain-related arrhythmias in fibrotically remodelled hearts.

Toxicity of aged gasoline exhaust particles to normal and diseased airway epithelia

Particulate matter (PM) pollution is a leading cause of premature death, particularly in those with pre-existing lung disease. A causative link between particle properties and adverse health effects remains unestablished mainly due to complex and variable physico-chemical PM parameters. Controlled laboratory experiments are required. Generating atmospherically realistic Aerosols and performing cell-exposure studies at relevant particle-doses are challenging. Here we examine gasoline-exhaust particle toxicity from a Euro-5 passenger car in a uniquely realistic exposure scenario, combining a smog chamber simulating atmospheric ageing, an aerosol enrichment System varying particle number concentration independent of particle chemistry, and an aerosol Deposition chamber physiologically delivering particles on air-liquid interface (ALI) cultures reproducing normal and susceptible health status. Gasoline-exhaust is an important PM source with largely unknown health effects. We investigated acute responses of fully-differentiated normal, distressed (antibiotics treated) normal, and cystic fibrosis human bronchial epithelia (HBE), and a proliferating, single-cell type bronchial epithelial cell-line (BEAS-2B). We show that a single, short-term exposure to realistic doses of atmospherically-aged gasoline-exhaust particles impairs epithelial key-defence mechanisms, rendering it more vulnerable to subsequent hazards. We establish dose-response curves at realistic particle-concentration levels. Significant differences between cell models suggest the use of fully differentiated HBE is most appropriate in future toxicity studies.

Expression of the hyaluronan-mediated motility receptor RHAMM in tumor budding cells identifies aggressive colorectal cancers.

Expression of the hyaluronan-mediated motility receptor (RHAMM, CD168) predicts adverse clinicopathological features and decreased survival for colorectal cancer (CRC) patients. Using full tissue sections, we investigated the expression of RHAMM in tumor budding cells of 103 primary CRCs to characterize the biological processes driving single-cell invasion and early metastatic dissemination. RHAMM expression in tumor buds was analyzed with clinicopathological data, molecular features and survival. Tumor budding cells at the invasive front of CRC expressed RHAMM in 68% of cases. Detection of RHAMM-positive tumor budding cells was significantly associated with poor survival outcome (P = .0312), independent of TNM stage and adjuvant therapy in multivariate analysis (P = .0201). RHAMM-positive tumor buds were associated with frequent lymphatic invasion (P = .0007), higher tumor grade (P = .0296), and nodal metastasis (P = .0364). Importantly, the prognostic impact of RHAMM expression in tumor buds was maintained independently of the number of tumor buds found in an individual case (P = .0246). No impact of KRAS/BRAF mutation, mismatch repair deficiency and CpG island methylation was observed. RHAMM expression identifies an aggressive subpopulation of tumor budding cells and is an independent adverse prognostic factor for CRC patients. These data support ongoing efforts to develop RHAMM as a target for precision therapy.

Identification and characterisation of ncRNAs in dormant bacteria

Recent years have led to increasing interest and appreciation of the possible importance of single cell heterogeneity in various biological processes. One of the examples of phenotypic heterogeneity in bacterial populations is antibiotic tolerant persister cells. Such an antibiotic tolerance phenotype is of considerable clinical relevance since dormant bacteria can re-establish infections rapidly after the antibiotic treatment has been terminated. Up to now mechanisms for establishing the persistence phenomenon in bacteria have remained largely enigmatic. Persisters are cells considered to be in a dormant state with down regulated gene expression. Only recently small regulatory RNAs (sRNAs) have been appreciated as important regulators of gene expression in response to environmental stimuli and several theoretical studies have suggested a possible involvement of sRNAs in the mechanisms of regulated heterogeneity in bacteria. We have experimentally addressed this potential link between sRNAs and persistence/dormancy in E. coli as an example of heterogeneity. Beside classical sRNAs we are focusing also on sRNAs directly associating with and possibly regulating the ribosome, the central enzyme of gene expression. The persister and dormant cell specific sRNA profile is studied by the comparative analysis of sRNA profile changes of the whole bacterial population after antibiotic killing. From RNA-Seq data ~ 25 000 potentially stable RNA fragments were identified and initial analysis predicted ~300 of them to be dormant/persister cell specific. After further evaluation the most prominent dormant/persister cell specific sRNAs are functionally characterized and their potential role in the persistence/dormancy will be evaluated by applying genetic, molecular and biochemical tools. The potential results of this project will provide a better understanding on the molecular mechanism of bacterial persistence/dormancy and on the role of ribosome-bound sRNA molecules in fine-tuning gene expression.

Identification and characterisation of ncRNAs in dormant bacteria

Recent years have led to increasing interest and appreciation of the possible importance of single cell heterogeneity in various biological processes. One of the examples of phenotypic heterogeneity in bacterial populations is antibiotic tolerant persister cells. Such an antibiotic tolerance phenotype is of considerable clinical relevance since dormant bacteria can re-establish infections rapidly after the antibiotic treatment has been terminated. Up to now mechanisms for establishing the persistence phenomenon in bacteria have remained largely enigmatic. Persisters are cells considered to be in a dormant state with down regulated gene expression. Only recently small regulatory RNAs (sRNAs) have been appreciated as important regulators of gene expression in response to environmental stimuli and several theoretical studies have suggested a possible involvement of sRNAs in the mechanisms of regulated heterogeneity in bacteria. We have experimentally addressed this potential link between sRNAs and persistence/dormancy in E. coli as an example of heterogeneity. Beside classical sRNAs we are focusing also on sRNAs directly associating with and possibly regulating the ribosome, the central enzyme of gene expression. The persister and dormant cell specific sRNA profile is studied by the comparative analysis of sRNA profile changes of the whole bacterial population after antibiotic killing. From RNA-Seq data ~ 25 000 potentially stable RNA fragments were identified and initial analysis predicted ~300 of them to be dormant/persister cell specific. After further evaluation the most prominent dormant/persister cell specific sRNAs are functionally characterized and their potential role in the persistence/dormancy will be evaluated by applying genetic, molecular and biochemical tools. The potential results of this project will provide a better understanding on the molecular mechanism of bacterial persistence/dormancy and on the role of ribosome-bound sRNA molecules in fine-tuning gene expression.

Identification and characterisation of ncRNAs in dormant bacteria

Recent years have led to increasing interest and appreciation of the possible importance of single cell heterogeneity in various biological processes. One of the examples of phenotypic heterogeneity in bacterial populations is antibiotic tolerant persister cells. Such an antibiotic tolerance phenotype is of considerable clinical relevance since dormant bacteria can re-establish infections rapidly after the antibiotic treatment has been terminated. Up to now mechanisms for establishing the persistence phenomenon in bacteria have remained largely enigmatic. Persisters are cells considered to be in a dormant state with down regulated gene expression. Only recently small regulatory RNAs (sRNAs) have been appreciated as important regulators of gene expression in response to environmental stimuli and several theoretical studies have suggested a possible involvement of sRNAs in the mechanisms of regulated heterogeneity in bacteria. We have experimentally addressed this potential link between sRNAs and persistence/dormancy in E. coli as an example of heterogeneity. Beside classical sRNAs we are focusing also on sRNAs directly associating with and possibly regulating the ribosome, the central enzyme of gene expression. The persister and dormant cell specific sRNA profile is studied by the comparative analysis of sRNA profile changes of the whole bacterial population after antibiotic killing. From RNA-Seq data ~ 25 000 potentially stable RNA fragments were identified and initial analysis predicted ~300 of them to be dormant/persister cell specific. After further evaluation the most prominent dormant/persister cell specific sRNAs are functionally characterized and their potential role in the persistence/dormancy will be evaluated by applying genetic, molecular and biochemical tools. The potential results of this project will provide a better understanding on the molecular mechanism of bacterial persistence/dormancy and on the role of ribosome-bound sRNA molecules in fine-tuning gene expression.

Tissue Crowding Induces Caspase-Dependent Competition for Space.

Regulation of tissue size requires fine tuning at the single-cell level of proliferation rate, cell volume, and cell death. Whereas the adjustment of proliferation and growth has been widely studied [1, 2, 3, 4 and 5], the contribution of cell death and its adjustment to tissue-scale parameters have been so far much less explored. Recently, it was shown that epithelial cells could be eliminated by live-cell delamination in response to an increase of cell density [6]. Cell delamination was supposed to occur independently of caspase activation and was suggested to be based on a gradual and spontaneous disappearance of junctions in the delaminating cells [6]. Studying the elimination of cells in the midline region of the Drosophila pupal notum, we found that, contrary to what was suggested before, Caspase 3 activation precedes and is required for cell delamination. Yet, using particle image velocimetry, genetics, and laser-induced perturbations, we confirmed [ 6] that local tissue crowding is necessary and sufficient to drive cell elimination and that cell elimination is independent of known fitness-dependent competition pathways [ 7, 8 and 9]. Accordingly, activation of the oncogene Ras in clones was sufficient to compress the neighboring tissue and eliminate cells up to several cell diameters away from the clones. Mechanical stress has been previously proposed to contribute to cell competition [ 10 and 11]. These results provide the first experimental evidences that crowding-induced death could be an alternative mode of super-competition, namely mechanical super-competition, independent of known fitness markers [ 7, 8 and 9], that could promote tumor growth.

XYLITOL METABOLISM IN RAT LIVER AND HEPATOCELLULAR CARCINOMAS

The significance of nutritional factors in cancer research has been strongly emphasized. Such research is concerned not only with epidemiological effects relative to dietary factors on the causation of cancer, but with nutritional effects as an energy source on the prevention of cancer. Many studies speculate that the energy flow between tumor and host can be regulated by dietary intake. However, little knowledge on the comparison of the specific nutritional and energy requirements of different cells and tissues is available. Most popular and essential energy sources for the body are the carbohydrates. Among them, xylitol is known as efficient an energy source as glucose. In carbohydrate metabolism, glycolysis is one of the major energy producing pathways. However, recently the existence of an alternate catabolic pathway in mammals for carbohydrate besides glycolysis, i.e. bypass through triosephosphates to lactate via methylglyoxal has been suggested. This bypass was implicated to regulate glycolysis and also be responsible for the fluctuation in the levels of a regulator of cell growth. Methylglyoxal itself is known as a cancerostatic agent. The alterations of biochemical parameters in xylitol metabolism in animals indicated that xylitol may be metabolized through a methylglyoxal pathway.^ To elucidate the biological effect of xylitol as an energy source and the biological effect of its metabolites as a cancerostatis agent, the mode and extent of metabolism must be understood in tumor-bearing animals. Differential utilization of xylitol and glucose, if any, between tumor and host in such animals may exert tissue selective effects on both in terms of methylglyoxal formation and energy provision. The aim of this work was to assess the extent to which the differential utilization of xylitol might be used to generate different metabolic pathways in tumor and host, and to consider a role of nutrition in cancer.^ The results disclose that the existence of a pathway for biological methylglyoxal formation in normal rat liver has been confirmed in single cell suspension; the metabolic significance of the methylglyoxal pathway in the metabolism of glucose and xylitol has been evaluated quantitatively in normal rat liver and the differential metabolism of glucose and xylitol through overall catabolic pathways of carbohydrates has been studied in normal hepatic cells, AS-30D hepatoma and other several hepatoma lines. ^

CHOLINERGIC MODULATION OF THALAMIC INPUT INTO THE CINGULATE CORTEX

This research demonstrates cholinergic modulation of thalamic input into the limbic cortex. A projection from the mediodorsal thalamus (MD) to the anterior cingulate cortex was defined anatomically and physiologically. Injections of horse-radish peroxidase into the anterior cingulate cortex labels neurons in the lateral, parvocellular, region of MD. Electrical Stimulation of this area produces a complex field potential in the anterior cingulate cortex which was further characterized by current density analysis and single cell recordings.^ The monsynaptic component of the response was identified as a large negative field which is maximal in layer IV of the anterior cingulate cortex. This response shows remarkable tetanic potentiation of frequencies near 7 Hz. During a train of 50 or more stimuli, the response would grow quickly and remain at a fairly stable potentiated level throughout the train.^ Cholinergic modulation of this thalamic response was demonstrated by iontophoretic application of the cholinergic agonist carbachol decreased the effectiveness of the thalamic imput by rapidly attenuation the response during a train of stimuli. The effect was apparently mediated by muscarinic receptors since the effect of carbachol was blocked by atropine but not by hexamethonium.^ To determine the source of the cingulate cortex cholinergic innervation, lesions were made in the anterior and medial thalamus and in the nucleus of the diagonal band of Broca. The effects of these lesions on choline acetyltranferase activity in the cingulate cortex were determined by a micro-radio-enzymatical assay. Only the lesions of the nucleus of the diagonal band significantly decreased the choline acetyltransferase activity in the cingulate cortex regions. Therefore, the diagonal band appears to be a major source of sensory cholinergic innervation and may be involved in gating of sensory information from the thalamus into the limbic cortex. Attempts to modulate the cingulate response to MD stimulation with electrical stimulation of the diagonal band, however were not successful.^

Cellular mechanisms of operant and classical conditioning

The ability to associate a predictive stimulus with a subsequent salient event (i.e., classical conditioning) and the ability to associate an expressed behavior with the consequences (i.e., operant conditioning) allow for a predictive understanding of a changing environment. Although they are operationally distinct, there has been considerable debate whether at some fundamental level classical and operant conditioning are mechanistically distinct or similar. Feeding behavior of Aplysia (i.e., biting) was chosen as the model system and was successfully conditioned with appetitive forms of both operant and classical conditioning. The neuronal circuitry responsible for feeding is well understood and is suitable for cellular analyses, thus providing for a mechanistic comparison between these two forms of associative learning. ^ Neuron B51 is part of the feeding circuitry of Aplysia and is critical for the expression of ingestive behaviors. B51 also is a locus of plasticity following both operant and classical conditioning. Both in vivo and in vitro operant conditioning increased the input resistance and the excitability of B51. No pairing-specific changes in the input resistance were observed following both in vivo and in vitro classical conditioning. However, classical conditioning decreased the excitability of B51. Thus, both operant and classical conditioning modified the threshold level for activation of neuron B51, but in opposite directions, revealing key differences in the cellular mechanisms underlying these two forms of associative learning. ^ Next, the cellular mechanisms underlying operant conditioning were investigated in more detail using a single-cell analogue. The single-cell analogue successfully recapitulated the previous in vivo and in vitro operant conditioning results by increasing the input resistance and the excitability of B51. Both PKA and PKC were necessary for operant conditioning. Dopamine appears to be the transmitter mediating the reinforcement signal in this form of conditioning. A D1 dopamine receptor antibody revealed that the D1receptor localizes to the axon hillock, which is also the region that gives the strongest response when iontophoresing dopamine. ^ The studies presented herein, thus, provide for a greater understanding of the mechanisms underlying both of these forms of associative learning and demonstrate that they likely operate through distinct cellular mechanisms. ^

IDENTIFICATION AND CHARACTERIZATION OF DISTINCT POPULATIONS OF CLONOGENIC BONE MARROW STROMAL CELLS CAPABLE OF TRANSFERRING THE HEMATOPOIETIC MICROENVIRONMENT IN VIVO AND SUPPORTING LT-HSCs IN VITRO

Bone marrow (BM) stromal cells are ascribed two key functions, 1) stem cells for non-hematopoietic tissues (MSC) and 2) as components of the hematopoietic stem cell niche. Current approaches studying the stromal cell system in the mouse are complicated by the low yield of clonogenic progenitors (CFU-F). Given the perivascular location of MSC in BM, we developed an alternative methodology to isolate MSC from mBM. An intact ‘plug’ of bone marrow is expelled from bones and enzymatically disaggregated to yield a single cell suspension. The recovery of CFU-F (1917.95+199) reproducibly exceeds that obtained using the standard BM flushing technique (14.32+1.9) by at least 2 orders of magnitude (P<0.001; N = 8) with an accompanying 196-fold enrichment of CFU-F frequency. Purified BM stromal and vascular endothelial cell populations are readily obtained by FACS. A detailed immunophenotypic analysis of lineage depleted BM identified PDGFRαβPOS stromal cell subpopulations distinguished by their expression of CD105. Both subpopulations retained their original phenotype of CD105 expression in culture and demonstrate MSC properties of multi-lineage differentiation and the ability to transfer the hematopoietic microenvironment in vivo. To determine the capacity of either subpopulation to support long-term multi-lineage reconstituting HSCs, we fractionated BM stromal cells into either the LinNEGPDGFRαβPOSCD105POS and LINNEGPDGFRαβPOSCD105LOW/- populations and tested their capacity to support LT-HSC by co-culturing each population with either 1 or 10 HSCs for 10 days. Following the 10 day co-culture period, both populations supported transplantable HSCs from 10 cells/well co-cultures demonstrating high levels of donor repopulation with an average of 65+23.6% chimerism from CD105POS co-cultures and 49.3+19.5% chimerism from the CD105NEG co-cultures. However, we observed a significant difference when mice were transplanted with the progeny of a single co-cultured HSC. In these experiments, CD105POS co-cultures (100%) demonstrated long-term multi- lineage reconstitution, while only 4 of 8 mice (50%) from CD105NEG -single HSC co-cultures demonstrated long-term reconstitution, suggesting a more limited expansion of functional stem cells. Taken together, these results demonstrate that the PDGFRαβCD105POS stromal cell subpopulation is distinguished by a unique capacity to support the expansion of long-term reconstituting HSCs in vitro.

DEVELOPMENT OF NOVEL METHODS TO MINIMIZE THE IMPACT OF SEQUENCING ERRORS IN THE NEXT-GENERATION SEQUENCING DATA ANALYSIS

Next-generation sequencing (NGS) technology has become a prominent tool in biological and biomedical research. However, NGS data analysis, such as de novo assembly, mapping and variants detection is far from maturity, and the high sequencing error-rate is one of the major problems. . To minimize the impact of sequencing errors, we developed a highly robust and efficient method, MTM, to correct the errors in NGS reads. We demonstrated the effectiveness of MTM on both single-cell data with highly non-uniform coverage and normal data with uniformly high coverage, reflecting that MTM’s performance does not rely on the coverage of the sequencing reads. MTM was also compared with Hammer and Quake, the best methods for correcting non-uniform and uniform data respectively. For non-uniform data, MTM outperformed both Hammer and Quake. For uniform data, MTM showed better performance than Quake and comparable results to Hammer. By making better error correction with MTM, the quality of downstream analysis, such as mapping and SNP detection, was improved. SNP calling is a major application of NGS technologies. However, the existence of sequencing errors complicates this process, especially for the low coverage (