Mathieu Lajoie

The Artiodactyl APOBEC3 Innate Immune Repertoire Shows Evidence for a Multi-Functional Domain Organization that Existed in the Ancestor of Placental Mammals

BackgroundAPOBEC3 (A3) proteins deaminate DNA cytosines and block the replication of retroviruses and retrotransposons. Each A3 gene encodes a protein with one or two conserved zinc-coordinating motifs (Z1, Z2 or Z3). The presence of one A3 gene in mice (Z2–Z3) and seven in humans, A3A-H (Z1a, Z2a-Z1b, Z2b, Z2c-Z2d, Z2e-Z2f, Z2g-Z1c, Z3), suggests extraordinary evolutionary flexibility. To gain insights into the mechanism and timing of A3 gene expansion and into the functional modularity of these genes, we analyzed the genomic sequences, expressed cDNAs and activities of the full A3 repertoire of three artiodactyl lineages: sheep, cattle and pigs.ResultsSheep and cattle have three A3 genes, A3Z1, A3Z2 and A3Z3, whereas pigs only have two, A3Z2 and A3Z3. A comparison between domestic and wild pigs indicated that A3Z1 was deleted in the pig lineage. In all three species, read-through transcription and alternative splicing also produced a catalytically active double domain A3Z2-Z3 protein that had a distinct cytoplasmic localization. Thus, the three A3 genes of sheep and cattle encode four conserved and active proteins. These data, together with phylogenetic analyses, indicated that a similar, functionally modular A3 repertoire existed in the common ancestor of artiodactyls and primates (i.e., the ancestor of placental mammals). This mammalian ancestor therefore possessed the minimal A3 gene set, Z1-Z2-Z3, required to evolve through a remarkable series of eight recombination events into the present day eleven Z domain human repertoire.ConclusionThe dynamic recombination-filled history of the mammalian A3 genes is consistent with the modular nature of the locus and a model in which most of these events (especially the expansions) were selected by ancient pathogenic retrovirus infections.

Duplication and inversion history of a tandemly repeated genes family

Given a phylogenetic tree for a family of tandemly repeated genes and their signed order on the chromosome, we aim to find the minimum number of inversions compatible with an evolutionary history of this family. This is the first attempt to account for inversions in an evolutionary model of tandemly repeated genes. We present a branch-and-bound algorithm that finds the exact solution, and a polynomial-time heuristic based on the breakpoint distance. We show, on simulated data, that those algorithms can be used to improve phylogenetic inference of tandemly repeated gene families. An application on a published phylogeny of KRAB zinc finger genes is presented.

Inferring ancestral gene orders for a family of tandemly arrayed genes.

Tandemly arrayed genes (TAG) constitute a large fraction of most genomes and play important biological roles. They evolve through unequal recombination, which places duplicated genes next to the original ones (tandem duplications). Many algorithms have been proposed to infer a tandem duplication history for a TAG cluster. However, the presence of different transcriptional orientations in many clusters highlights the fact that processes such as inversions also contribute to their evolution. Moreover, existing algorithms are restricted to the study of TAGs evolution in a single species (only paralogous genes are considered). To circumvent these limitations, we consider an evolutionary model for TAGs involving duplication, gene loss, inversion, and speciation events. A general framework to infer ancestral gene orders that minimize the number of inversions in the whole evolutionary history is presented. At the methodological level, this paper integrates three approaches to genome evolution: the duplication tree reconstruction, the gene tree/species tree reconciliation theory, and the concept of inversion median used in order-based phylogeny reconstruction. An application on a cluster of olfactory receptor genes in four mammals is presented.

Inferring the Evolutionary History of Gene Clusters from Phylogenetic and Gene Order Data

Gene duplication is frequent within gene clusters and plays a fundamental role in evolution by providing a source of new genetic material upon which natural selection can act. Although classical phylogenetic inference methods provide some insight into the evolutionary history of a gene cluster, they are not sufficient alone to differentiate single- from multiple gene duplication events and to answer other questions regarding the nature and size of evolutionary events. In this paper, we present an algorithm allowing to infer a set of optimal evolutionary histories for a gene cluster in a single species, according to a general cost model involving variable length duplications (in tandem or inverted), deletions, and inversions. We applied our algorithm to the human olfactory receptor and protocadherin gene clusters, showing that the duplication size distribution differs significantly between the two gene families. The algorithm is available through a web interface at http://www-lbit.iro.umontreal.ca/DILTAG/.

Promiscuous DNA-binding of a mutant zinc finger protein corrupts the transcriptome and diminishes cell viability

Abstract The rules of engagement between zinc finger transcription factors and DNA have been partly defined by in vitro DNA-binding and structural studies, but less is known about how these rules apply in vivo. Here, we demonstrate how a missense mutation in the second zinc finger of Krüppel-like factor-1 (KLF1) leads to degenerate DNA-binding specificity in vivo, resulting in ectopic transcription and anemia in the Nan mouse model. We employed ChIP-seq and 4sU-RNA-seq to identify aberrant DNA-binding events genome wide and ectopic transcriptional consequences of this binding. We confirmed novel sequence specificity of the mutant recombinant zinc finger domain by performing biophysical measurements of in vitro DNA-binding affinity. Together, these results shed new light on the mechanisms by which missense mutations in DNA-binding domains of transcription factors can lead to autosomal dominant diseases.

Evolution of tandemly repeated sequences through duplication and inversion

Given a phylogenetic tree T for a family of tandemly repeated genes and their signed order O on the chromosome, we aim to find the minimum number of inversions compatible with an evolutionary history of this family. This is the first attempt to account for inversions in an evolutionary model of tandemly repeated genes. We present a time-efficient branch-and-bound algorithm and show, using simulated data, that it can be used to detect “wrong” phylogenies among a set of putative ones for a given gene family. An application on a published phylogeny of KRAB zinc finger genes is presented.

A childhood acute lymphoblastic leukemia-specific lncRNA implicated in prednisolone resistance, cell proliferation, and migration

Childhood acute lymphoblastic leukemia (cALL) is the most common pediatric cancer and, despite an 85% cure rate, still represents a major cause of disease-related death in children. Recent studies have implicated long non-coding RNAs (lncRNAs) in cALL etiology, progression, and treatment response. However, barring some exceptions little is known about the functional impact of lncRNAs on cancer biology, which limits their potential as potential therapeutic targets. We wanted to investigate the functional role of lncRNAs identified as specifically overexpressed in pre-B cALL by whole-transcriptome sequencing. Here we report five lncRNAs specifically upregulated in pre-B cALL that had significant impacts on cancer hallmark traits such as cell proliferation, migration, apoptosis, and treatment response. In particular, silencing of the RP11-137H2.4 lncRNA effectively restored normal glucocorticoid (GC) response in a GC-resistant pre-B cALL cell line and specifically modulated expression of members of both the NRAS/BRAF/NF-?B MAPK cascade and cell cycle pathways. Since GC form the cornerstone of cALL chemotherapy and resistance in cALL confers a dismal prognosis, characterizing RP11-137H2.4sexact role and function in this process will be critical to the development of new therapeutic approaches to overcome GC resistance in children treated for cALL.

An overlapping set of genes is regulated by both NFIB and the glucocorticoid receptor during lung maturation

BackgroundLung maturation is a late fetal developmental event in both mice and humans. Because of this, lung immaturity is a serious problem in premature infants. Disruption of genes for either the glucocorticoid receptor (Nr3c1) or the NFIB transcription factors results in perinatal lethality due to lung immaturity. In both knockouts, the phenotype includes excess cell proliferation, failure of saccularization and reduced expression of markers of epithelial differentiation. This similarity suggests that the two genes may co-regulate a specific set of genes essential for lung maturation.ResultsWe analyzed the roles of these two transcription factors in regulating transcription using ChIP-seq data for NFIB, and RNA expression data and motif analysis for both. Our new ChIP-seq data for NFIB in lung at E16.5 shows that NFIB binds to a NFI motif. This motif is over-represented in the promoters of genes that are under-expressed in Nfib-KO mice at E18.5, suggesting an activator role for NFIB. Using available microarray data from Nr3c1-KO mice, we further identified 52 genes that are under-expressed in both Nfib and Nr3c1 knockouts, an overlap which is 13.1 times larger than what would be expected by chance. Finally, we looked for enrichment of 738 recently published transcription factor motifs in the promoters of these putative target genes and found that the NFIB and glucocorticoid receptor motifs were among the most enriched, suggesting that a subset of these genes may be directly activated by Nfib and Nr3c1.ConclusionsOur data provide the first evidence for Nfib and Nr3c1 co-regulating genes related to lung maturation. They also establish that the in vivo DNA-binding specificity of NFIB is the same as previously seen in vitro, and highly similar to that of the other NFI-family members NFIA, NFIC and NFIX.

Recovering haplotype structure through recombination and gene conversion

MOTIVATION Understanding haplotype evolution subject to mutation, recombination and gene conversion is fundamental to understand genetic specificities of human populations and hereditary bases of complex disorders. The goal of this project is to develop new algorithmic tools assisting the reconstruction of historical relationships between haplotypes and the inference of haplotypes from genotypes. RESULTS We present two new algorithms. The first one finds an optimal pathway of mutations, recombinations and gene conversions leading to a given haplotype of size m from a population of h haplotypes. It runs in time O(mhs(2)), where s is the maximum number of contiguous sites that can be exchanged in a single gene conversion. The second one finds an optimal pathway of mutations and recombinations leading to a given genotype, and runs in time O(mh(2)). Both algorithms are based on a penalty score model and use a dynamic programming approach. We apply the second one to the problem of inferring haplotypes from genotypes, and show how it can be used as an independent tool, or to improve the performance of existing methods. AVAILABILITY The algorithms have been implemented in JAVA and are available on request.

CLIC5: a novel ETV6 target gene in childhood acute lymphoblastic leukemia

The most common rearrangement in childhood precursor B-cell acute lymphoblastic leukemia is the t(12;21)(p13;q22) translocation resulting in the ETV6-AML1 fusion gene. A frequent concomitant event is the loss of the residual ETV6 allele suggesting a critical role for the ETV6 transcriptional repressor in the etiology of this cancer. However, the precise mechanism through which loss of functional ETV6 contributes to disease pathogenesis is still unclear. To investigate the impact of ETV6 loss on the transcriptional network and to identify new transcriptional targets of ETV6, we used whole transcriptome analysis of both pre-B leukemic cell lines and patients combined with chromatin immunoprecipitation. Using this integrative approach, we identified 4 novel direct ETV6 target genes: CLIC5, BIRC7, ANGPTL2 and WBP1L. To further evaluate the role of chloride intracellular channel protein CLIC5 in leukemogenesis, we generated cell lines overexpressing CLIC5 and demonstrated an increased resistance to hydrogen peroxide-induced apoptosis. We further described the implications of CLIC5’s ion channel activity in lysosomal-mediated cell death, possibly by modulating the function of the transferrin receptor with which it colocalizes intracellularly. For the first time, we showed that loss of ETV6 leads to significant overexpression of CLIC5, which in turn leads to decreased lysosome-mediated apoptosis. Our data suggest that heightened CLIC5 activity could promote a permissive environment for oxidative stress-induced DNA damage accumulation, and thereby contribute to leukemogenesis.