INTEGRATIVE BIOMARKER IDENTIFICATION AND CLASSIFICATION USING HIGH THROUGHPUT ASSAYS

It is well accepted that tumorigenesis is a multi-step procedure involving aberrant functioning of genes regulating cell proliferation, differentiation, apoptosis, genome stability, angiogenesis and motility. To obtain a full understanding of tumorigenesis, it is necessary to collect information on all aspects of cell activity. Recent advances in high throughput technologies allow biologists to generate massive amounts of data, more than might have been imagined decades ago. These advances have made it possible to launch comprehensive projects such as (TCGA) and (ICGC) which systematically characterize the molecular fingerprints of cancer cells using gene expression, methylation, copy number, microRNA and SNP microarrays as well as next generation sequencing assays interrogating somatic mutation, insertion, deletion, translocation and structural rearrangements. Given the massive amount of data, a major challenge is to integrate information from multiple sources and formulate testable hypotheses. This thesis focuses on developing methodologies for integrative analyses of genomic assays profiled on the same set of samples. We have developed several novel methods for integrative biomarker identification and cancer classification. We introduce a regression-based approach to identify biomarkers predictive to therapy response or survival by integrating multiple assays including gene expression, methylation and copy number data through penalized regression. To identify key cancer-specific genes accounting for multiple mechanisms of regulation, we have developed the integIRTy software that provides robust and reliable inferences about gene alteration by automatically adjusting for sample heterogeneity as well as technical artifacts using Item Response Theory. To cope with the increasing need for accurate cancer diagnosis and individualized therapy, we have developed a robust and powerful algorithm called SIBER to systematically identify bimodally expressed genes using next generation RNAseq data. We have shown that prediction models built from these bimodal genes have the same accuracy as models built from all genes. Further, prediction models with dichotomized gene expression measurements based on their bimodal shapes still perform well. The effectiveness of outcome prediction using discretized signals paves the road for more accurate and interpretable cancer classification by integrating signals from multiple sources.

Un modelo de encaminamiento en las redes móviles ad hoc con QoS

En la última década ha aumentado en gran medida el interés por las redes móviles Ad Hoc. La naturaleza dinámica y sin infraestructura de estas redes, exige un nuevo conjunto de algoritmos y estrategias para proporcionar un servicio de comunicación fiable extremo a extremo. En el contexto de las redes móviles Ad Hoc, el encaminamiento surge como una de las áreas más interesantes para transmitir información desde una fuente hasta un destino, con la calidad de servicio de extremo a extremo. Debido a las restricciones inherentes a las redes móviles, los modelos de encaminamiento tradicionales sobre los que se fundamentan las redes fijas, no son aplicables a las redes móviles Ad Hoc. Como resultado, el encaminamiento en redes móviles Ad Hoc ha gozado de una gran atención durante los últimos años. Esto ha llevado al acrecentamiento de numerosos protocolos de encaminamiento, tratando de cubrir con cada uno de ellos las necesidades de los diferentes tipos de escenarios. En consecuencia, se hace imprescindible estudiar el comportamiento de estos protocolos bajo configuraciones de red variadas, con el fin de ofrecer un mejor encaminamiento respecto a los existentes. El presente trabajo de investigación muestra precisamente una solución de encaminamiento en las redes móviles Ad Hoc. Dicha solución se basa en el mejoramiento de un algoritmo de agrupamiento y la creación de un modelo de encaminamiento; es decir, un modelo que involucra la optimización de un protocolo de enrutamiento apoyado de un mecanismo de agrupación. El algoritmo mejorado, denominado GMWCA (Group Management Weighted Clustering Algorithm) y basado en el WCA (Weighted Clustering Algorithm), permite calcular el mejor número y tamaño de grupos en la red. Con esta mejora se evitan constantes reagrupaciones y que los jefes de clústeres tengan más tiempo de vida intra-clúster y por ende una estabilidad en la comunicación inter-clúster. En la tesis se detallan las ventajas de nuestro algoritmo en relación a otras propuestas bajo WCA. El protocolo de enrutamiento Ad Hoc propuesto, denominado QoS Group Cluster Based Routing Protocol (QoSG-CBRP), utiliza como estrategia el empleo de clúster y jerarquías apoyada en el algoritmo de agrupamiento. Cada clúster tiene un jefe de clúster (JC), quien administra la información de enrutamiento y la envía al destino cuando esta fuera de su área de cobertura. Para evitar que haya constantes reagrupamientos y llamados al algoritmo de agrupamiento se consideró agregarle un jefe de cluster de soporte (JCS), el que asume las funciones del JC, siempre y cuando este haya roto el enlace con los otros nodos comunes del clúster por razones de alejamiento o por desgaste de batería. Matemáticamente y a nivel de algoritmo se han demostrado las mejoras del modelo propuesto, el cual ha involucrado el mejoramiento a nivel de algoritmo de clustering y del protocolo de enrutamiento. El protocolo QoSG-CBRP, se ha implementado en la herramienta de simulación Network Simulator 2 (NS2), con la finalidad de ser comparado con el protocolo de enrutamiento jerárquico Cluster Based Routing Protocol (CBRP) y con un protocolo de enrutamiento Ad Hoc reactivo denominado Ad Hoc On Demand Distance Vector Routing (AODV). Estos protocolos fueron elegidos por ser los que mejor comportamiento presentaron dentro de sus categorías. Además de ofrecer un panorama general de los actuales protocolos de encaminamiento en redes Ad Hoc, este proyecto presenta un procedimiento integral para el análisis de capacidades de la propuesta del nuevo protocolo con respecto a otros, sobre redes que tienen un alto número de nodos. Estas prestaciones se miden en base al concepto de eficiencia de encaminamiento bajo parámetros de calidad de servicio (QoS), permitiendo establecer el camino más corto posible entre un nodo origen y un nodo destino. Con ese fin se han realizado simulaciones con diversos escenarios para responder a los objetivos de la tesis. La conclusiones derivadas del análisis de los resultados permiten evaluar cualitativamente las capacidades que presenta el protocolo dentro del modelo propuesto, al mismo tiempo que avizora un atractivo panorama en líneas futuras de investigación. ABSTRACT In the past decade, the interest in mobile Ad Hoc networks has greatly increased. The dynamic nature of these networks without infrastructure requires a new set of algorithms and strategies to provide a reliable end-to-end communication service. In the context of mobile Ad Hoc networks, routing emerges as one of the most interesting areas for transmitting information from a source to a destination, with the quality of service from end-to-end. Due to the constraints of mobile networks, traditional routing models that are based on fixed networks are not applicable to Ad Hoc mobile networks. As a result, the routing in mobile Ad Hoc networks has experienced great attention in recent years. This has led to the enhancement of many routing protocols, trying to cover with each one of them, the needs of different types of scenarios. Consequently, it is essential to study the behavior of these protocols under various network configurations, in order to provide a better routing scheme. Precisely, the present research shows a routing solution in mobile Ad Hoc networks. This solution is based on the improvement of a clustering algorithm, and the creation of a routing model, ie a model that involves optimizing a routing protocol with the support of a grouping mechanism. The improved algorithm called GMWCA (Group Management Weighted Clustering Algorithm) and based on the WCA (Weighted Clustering Algorithm), allows to calculate the best number and size of groups in the network. With this enhancement, constant regroupings are prevented and cluster heads are living longer intra-cluster lives and therefore stability in inter-cluster communication. The thesis details the advantages of our algorithm in relation to other proposals under WCA. The Ad Hoc routing protocol proposed, called QoS Group Cluster Based Routing Protocol (QoSG-CBRP), uses a cluster-employment strategy and hierarchies supported by the clustering algorithm. Each cluster has a cluster head (JC), who manages the routing information and sends it to the destination when is out of your coverage area. To avoid constant rearrangements and clustering algorithm calls, adding a support cluster head (JCS) was considered. The JCS assumes the role of the JC as long as JC has broken the link with the other nodes in the cluster for common restraining reasons or battery wear. Mathematically and at an algorithm level, the improvements of the proposed model have been showed, this has involved the improvement level clustering algorithm and the routing protocol. QoSG-CBRP protocol has been implemented in the simulation tool Network Simulator 2 (NS2), in order to be compared with the hierarchical routing protocol Cluster Based Routing Protocol (CBRP) and with the reactive routing protocol Ad Hoc On Demand Distance Vector Routing (AODV). These protocols were chosen because they showed the best individual performance in their categories. In addition to providing an overview of existing routing protocols in Ad Hoc networks, this project presents a comprehensive procedure for capacity analysis of the proposed new protocol with respect to others on networks that have a high number of nodes. These benefits are measured based on the concept of routing efficiency under the quality of service (QoS) parameters, thus allowing for the shortest possible path between a source node and a destination node. To meet the objectives of the thesis, simulations have been performed with different scenarios. The conclusions derived from the analysis of the results to assess qualitatively the protocol capabilities presented in the proposed model, while an attractive scenario for future research appears.

Genomic evolution during a 10,000-generation experiment with bacteria

Molecular methods are used widely to measure genetic diversity within populations and determine relationships among species. However, it is difficult to observe genomic evolution in action because these dynamics are too slow in most organisms. To overcome this limitation, we sampled genomes from populations of Escherichia coli evolving in the laboratory for 10,000 generations. We analyzed the genomes for restriction fragment length polymorphisms (RFLP) using seven insertion sequences (IS) as probes; most polymorphisms detected by this approach reflect rearrangements (including transpositions) rather than point mutations. The evolving genomes became increasingly different from their ancestor over time. Moreover, tremendous diversity accumulated within each population, such that almost every individual had a different genetic fingerprint after 10,000 generations. As has been often suggested, but not previously shown by experiment, the rates of phenotypic and genomic change were discordant, both across replicate populations and over time within a population. Certain pivotal mutations were shared by all descendants in a population, and these are candidates for beneficial mutations, which are rare and difficult to find. More generally, these data show that the genome is highly dynamic even over a time scale that is, from an evolutionary perspective, very brief.

Regulation of T cell receptor (TCR) β gene expression by CD3 complex signaling in immature thymocytes: Implications for TCRβ allelic exclusion

During αβ thymocyte development, clonotype-independent CD3 complexes are expressed at the cell surface before the pre-T cell receptor (TCR). Signaling through clonotype-independent CD3 complexes is required for expression of rearranged TCRβ genes. On expression of a TCRβ polypeptide chain, the pre-TCR is assembled, and TCRβ locus allelic exclusion is established. We investigated the putative contribution of clonotype-independent CD3 complex signaling to TCRβ locus allelic exclusion in mice single-deficient or double-deficient for CD3ζ/η and/or p56lck. These mice display defects in the expression of endogenous TCRβ genes in immature thymocytes, proportional to the severity of CD3 complex malfunction. Exclusion of endogenous TCRβ VDJ (variable, diversity, joining) rearrangements by a functional TCRβ transgene was severely compromised in the single-deficient and double-deficient mutant mice. In contrast to wild-type mice, most of the CD25+ double-negative (DN) thymocytes of the mutant mice failed to express the TCRβ transgene, suggesting defective expression of the TCRβ transgene similar to endogenous TCRβ genes. In the mutant mice, a proportion of CD25+ DN thymocytes that failed to express the transgene expressed endogenous TCRβ polypeptide chains. Many double-positive cells of the mutant mice coexpressed endogenous and transgenic TCRβ chains or more than one endogenous TCRβ chain. The data suggest that signaling through clonotype-independent CD3 complexes may contribute to allelic exclusion of the TCRβ locus by inducing the expression of rearranged TCRβ genes in CD25+ DN thymocytes.

A site-specific recombinase is required for competitive root colonization by Pseudomonas fluorescens WCS365

A colonization mutant of the efficient root-colonizing biocontrol strain Pseudomonas fluorescens WCS365 is described that is impaired in competitive root-tip colonization of gnotobiotically grown potato, radish, wheat, and tomato, indicating a broad host range mutation. The colonization of the mutant is also impaired when studied in potting soil, suggesting that the defective gene also plays a role under more natural conditions. A DNA fragment that is able to complement the mutation for colonization revealed a multicistronic transcription unit composed of at least six ORFs with similarity to lppL, lysA, dapF, orf235/233, xerC/sss, and the largely incomplete orf238. The transposon insertion in PCL1233 appeared to be present in the orf235/233 homologue, designated orf240. Introduction of a mutation in the xerC/sss homologue revealed that the xerC/sss gene homologue rather than orf240 is crucial for colonization. xerC in Escherichia coli and sss in Pseudomonas aeruginosa encode proteins that belong to the λ integrase family of site-specific recombinases, which play a role in phase variation caused by DNA rearrangements. The function of the xerC/sss homologue in colonization is discussed in terms of genetic rearrangements involved in the generation of different phenotypes, thereby allowing a bacterial population to occupy various habitats. Mutant PCL1233 is assumed to be locked in a phenotype that is not well suited to compete for colonization in the rhizosphere. Thus we show the importance of phase variation in microbe–plant interactions.

Self-antigen does not accelerate immature B cell apoptosis, but stimulates receptor editing as a consequence of developmental arrest

In pre-B lymphocytes, productive rearrangement of Ig light chain genes allows assembly of the B cell receptor (BCR), which selectively promotes further developmental maturation through poorly defined transmembrane signaling events. Using a novel in vitro system to study immune tolerance during development, we find that BCR reactivity to auto-antigen blocks this positive selection, preventing down-regulation of light chain gene recombination and promoting secondary light chain gene rearrangements that often alter BCR specificity, a process called receptor editing. Under these experimental conditions, self-antigen induces secondary light chain gene rearrangements in at least two-thirds of autoreactive immature B cells, but fails to accelerate cell death at this stage. These data suggest that in these cells the mechanism of immune tolerance is receptor selection rather than clonal selection.

Hodgkin and Reed–Sternberg cells in lymphocyte predominant Hodgkin disease represent clonal populations of germinal center-derived tumor B cells

Among the four subtypes of Hodgkin disease (HD), lymphocyte-predominant (LP) HD is now generally considered as a separate entity. The B cell nature of the typical Hodgkin and Reed–Sternberg (HRS) cells and their variants (L and H, lymphocytic and histiocytic cells) in LP HD has long been suspected, but the question of whether these cells represent a true tumor clone is unclear. We previously demonstrated clonal Ig gene rearrangements in one case of LP HD. In the present study, five cases of LP HD were analyzed by micromanipulation of single HRS cells from frozen tissue sections and DNA amplification of rearranged Ig heavy chain genes from those cells. Clonal V gene rearrangements harboring somatic mutations were detected in each case. In three cases ongoing somatic mutation was evident. This shows that HRS cells in LP HD are a clonal tumor population derived from germinal center B cells. The pattern of somatic mutation indicates that HRS cells in LP HD are selected for antibody expression. This, and the presence of ongoing mutation discriminates LP from classical HD.

Panhandle PCR strategy to amplify MLL genomic breakpoints in treatment-related leukemias

Panhandle PCR amplifies genomic DNA with known 5′ and unknown 3′ sequences from a template with an intrastrand loop schematically shaped like a pan with a handle. We used panhandle PCR to clone MLL genomic breakpoints in two pediatric treatment-related leukemias. The karyotype in a case of treatment-related acute lymphoblastic leukemia showed the t(4;11)(q21;q23). Panhandle PCR amplified the translocation breakpoint at position 2158 in intron 6 in the 5′ MLL breakpoint cluster region (bcr). The karyotype in a case of treatment-related acute myeloid leukemia was normal, but Southern blot analysis showed a single MLL gene rearrangement. Panhandle PCR amplified the breakpoint at position 1493 in MLL intron 6. Screening of somatic cell hybrid and radiation hybrid DNAs by PCR and reverse transcriptase-PCR analysis of the leukemic cells indicated that panhandle PCR identified a fusion of MLL intron 6 with a previously uncharacterized sequence in MLL intron 1, consistent with a partial duplication. In both cases, the breakpoints in the MLL bcr were in Alu repeats, and there were Alu repeats in proximity to the breakpoints in the partner DNAs, suggesting that Alu sequences were relevant to these rearrangements. This study shows that panhandle PCR is an effective method for cloning MLL genomic breakpoints in treatment-related leukemias. Analysis of additional pediatric cases will determine whether breakpoint distribution deviates from the predilection for 3′ distribution in the bcr that has been found in adult cases.

Highly conservative reciprocal translocations formed by apparent joining of exchanged DNA double-strand break ends

Chromosomal translocations induced by ionizing radiation and radiomimetic drugs are thought to arise by incorrect joining of DNA double-strand breaks. To dissect such misrepair events at a molecular level, large-scale, bleomycin-induced rearrangements in the aprt gene of Chinese hamster ovary D422 cells were mapped, the breakpoints were sequenced, and the original non-aprt parental sequences involved in each rearrangement were recovered from nonmutant cells. Of seven rearrangements characterized, six were reciprocal exchanges between aprt and unrelated sequences. Consistent with a mechanism involving joining of exchanged double-strand break ends, there was, in most cases, no homology between the two parental sequences, no overlap in sequences retained at the two newly formed junctions, and little or no loss of parental sequences (usually ≤2 bp) at the breakpoints. The breakpoints were strongly correlated (P < 0.0001) with expected sites of bleomycin-induced, double-strand breaks. Fluorescence in situ hybridization indicated that, in six of the mutants, the rearrangement was accompanied by a chromosomal translocation at the aprt locus, because upstream and downstream flanking sequences were detected on separate chromosomes. The results suggest that repair of free radical-mediated, double-strand breaks in confluence-arrested cells is effected by a conservative, homology-independent, end-joining pathway that does not involve single-strand intermediate and that misjoining of exchanged ends by this pathway can directly result in chromosomal translocations.

Induction of host signal transduction pathways by Helicobacter pylori

Adherence of Helicobacter pylori to cultured gastric epithelial cells is associated with several cellular events, including the tyrosine phosphorylation of a 145-kDa host protein; the reorganization of the host cell actin and associated cellular proteins, like vasodilator-stimulated phosphoprotein, adjacent to the attached bacterial cell; and the subsequent release of the cytokine, interleukin 8 (IL-8). H. pylori isolated from patients with ulcer disease and gastric cancer contain a DNA insertion, the cag pathogenicity island (PAI), that is not present in bacteria isolated from individuals with asymptomatic infection. Mutations in a number of PAI genes abolish tyrosine phosphorylation and IL-8 synthesis but not the cytoskeletal rearrangements. Kinase inhibition studies suggest there are two distinct pathways operative in stimulating IL-8 release from host cells and one of these H. pylori pathways is independent of the tyrosine phosphorylation step.

Developmental regulation of spine motility in the mammalian central nervous system

The function of dendritic spines, postsynaptic sites of excitatory input in the mammalian central nervous system (CNS), is still not well understood. Although changes in spine morphology may mediate synaptic plasticity, the extent of basal spine motility and its regulation and function remains controversial. We investigated spine motility in three principal neurons of the mouse CNS: cerebellar Purkinje cells, and cortical and hippocampal pyramidal neurons. Motility was assayed with time-lapse imaging by using two-photon microscopy of green fluorescent protein-labeled neurons in acute and cultured slices. In all three cell types, dendritic protrusions (filopodia and spines) were highly dynamic, exhibiting a diversity of morphological rearrangements over short (<1-min) time courses. The incidence of spine motility declined during postnatal maturation, but dynamic changes were still apparent in many spines in late-postnatal neurons. Although blockade or induction of neuronal activity did not affect spine motility, disruption of actin polymerization did. We hypothesize that this basal motility of dendritic protrusions is intrinsic to the neuron and underlies the heightened plasticity found in developing CNS.

Recombination-dependent deletion formation in mammalian cells deficient in the nucleotide excision repair gene ERCC1

Nucleotide excision repair proteins have been implicated in genetic recombination by experiments in Saccharomyces cerevisiae and Drosophila melanogaster, but their role, if any, in mammalian cells is undefined. To investigate the role of the nucleotide excision repair gene ERCC1, the hamster homologue to the S. cerevisiae RAD10 gene, we disabled the gene by targeted knockout. Partial tandem duplications of the adenine phosphoribosyltransferase (APRT) gene then were constructed at the endogenous APRT locus in ERCC1− and ERCC1+ cells. To detect the full spectrum of gene-altering events, we used a loss-of-function assay in which the parental APRT+ tandem duplication could give rise to APRT− cells by homologous recombination, gene rearrangement, or point mutation. Measurement of rates and analysis of individual APRT− products indicated that gene rearrangements (principally deletions) were increased at least 50-fold, whereas homologous recombination was affected little. The formation of deletions is not caused by a general effect of the ERCC1 deficiency on gene stability, because ERCC1− cell lines with a single wild-type copy of the APRT gene yielded no increase in deletions. Thus, deletion formation is dependent on the tandem duplication, and presumably the process of homologous recombination. Recombination-dependent deletion formation in ERCC1− cells is supported by a significant decrease in a particular class of crossover products that are thought to arise by repair of a heteroduplex intermediate in recombination. We suggest that the ERCC1 gene product in mammalian cells is involved in the processing of heteroduplex intermediates in recombination and that the misprocessed intermediates in ERCC1− cells are repaired by illegitimate recombination.

Chromosomal arm replacement generates a high level of intraspecific polymorphism in the terminal inverted repeats of the linear chromosomal DNA of Streptomyces ambofaciens

The chromosomal DNA of the bacteria Streptomyces ambofaciens DSM40697 is an 8-Mb linear molecule that ends in terminal inverted repeats (TIRs) of 210 kb. The sequences of the TIRs are highly variable between the different linear replicons of Streptomyces (plasmids or chromosomes). Two spontaneous mutant strains harboring TIRs of 480 and 850 kb were isolated. The TIR polymorphism seen is a result of the deletion of one chromosomal end and its replacement by 480 or 850 kb of sequence identical to the end of the undeleted chromosomal arm. Analysis of the wild-type sequences involved in these rearrangements revealed that a recombination event took place between the two copies of a duplicated DNA sequence. Each copy was mapped to one chromosomal arm, outside of the TIR, and encoded a putative alternative sigma factor. The two ORFs, designated hasR and hasL, were found to be 99% similar at the nucleotide level. The sequence of the chimeric regions generated by the recombination showed that the chromosomal structure of the mutant strains resulted from homologous recombination events between the two copies. We suggest that this mechanism of chromosomal arm replacement contributes to the rapid evolutionary diversification of the sequences of the TIR in Streptomyces.

Modification of EWS/WT1 functional properties by phosphorylation

In many human cancers, tumor-specific chromosomal rearrangements are known to create chimeric products with the ability to transform cells. The EWS/WT1 protein is such a fusion product, resulting from a t(11;22) chromosomal translocation in desmoplastic small round cell tumors, where 265 aa from the EWS amino terminus are fused to the DNA binding domain of the WT1 tumor suppressor gene. Herein, we find that EWS/WT1 is phosphorylated in vivo on serine and tyrosine residues and that this affects DNA binding and homodimerization. We also show that EWS/WT1 can interact with, and is a substrate for, modification on tyrosine residues by c-Abl. Tyrosine phosphorylation of EWS/WT1 by c-Abl negatively regulates its DNA binding properties. These results indicate that the biological activity of EWS/WT1 is closely linked to its phosphorylation status.

Genomic analysis of human and mouse TCL1 loci reveals a complex of tightly clustered genes

TCL1 and TCL1b genes on human chromosome 14q23.1 are activated in T cell leukemias by translocations and inversions at 14q32.1, juxtaposing them to regulatory elements of T cell receptor genes. In this report we present the cloning, mapping, and expression analysis of the human and murine TCL1/Tcl1 locus. In addition to TCL1 and TCL1b, the human locus contains two additional genes, TCL1-neighboring genes (TNG) 1 and 2, encoding proteins of 141 and 110 aa, respectively. Both genes show no homology to any known genes, but their expression profiles are very similar to those of TCL1 and TCL1b. TNG1 and TNG2 also are activated in T cell leukemias with rearrangements at 14q32.1. To aid in the development of a mouse model we also have characterized the murine Tcl1 locus and found five genes homologous to human TCL1b. Tcl1b1–Tcl1b5 proteins range from 117 to 123 aa and are 65–80% similar, but they show only a 30–40% similarity to human TCL1b. All five mouse Tcl1b and murine Tcl1 mRNAs are abundant in mouse oocytes and two-cell embryos but rare in various adult tissues and lymphoid cell lines. These data suggest a similar or complementary function of these proteins in early embryogenesis.