Yves Lemoine

Transcriptomic and genetic studies identify IL-33 as a candidate gene for Alzheimer's disease

The only recognized genetic determinant of the common forms of Alzheimers disease (AD) is the ɛ4 allele of the apolipoprotein E gene (APOE). To identify new candidate genes, we recently performed transcriptomic analysis of 2741 genes in chromosomal regions of interest using brain tissue of AD cases and controls. From 82 differentially expressed genes, 1156 polymorphisms were genotyped in two independent discovery subsamples (n=945). Seventeen genes exhibited at least one polymorphism associated with AD risk, and following correction for multiple testing, we retained the interleukin (IL)-33 gene. We first confirmed that the IL-33 expression was decreased in the brain of AD cases compared with that of controls. Further genetic analysis led us to select three polymorphisms within this gene, which we analyzed in three independent case–control studies. These polymorphisms and a resulting protective haplotype were systematically associated with AD risk in non-APOE ɛ4 carriers. Using a large prospective study, these associations were also detected when analyzing the prevalent and incident AD cases together or the incident AD cases alone. These polymorphisms were also associated with less cerebral amyloid angiopathy (CAA) in the brain of non-APOE ɛ4 AD cases. Immunohistochemistry experiments finally indicated that the IL-33 expression was consistently restricted to vascular capillaries in the brain. Moreover, IL-33 overexpression in cellular models led to a specific decrease in secretion of the Aβ40 peptides, the main CAA component. In conclusion, our data suggest that genetic variants in IL-33 gene may be associated with a decrease in AD risk potentially in modulating CAA formation.

Characterization of the biotin biosynthesis pathway in Saccharomyces cerevisiae and evidence for a cluster containing BIO5, a novel gene involved in vitamer uptake.

An engineered mutant of Saccharomyces cerevisiae affected in biotin biosynthesis has been isolated. This mutant allowed the characterization of a bio cluster (BIO3-4-5). We demonstrate that BIO3 (YNR058w) and BIO4 (YNR057c) encode, respectively, a 7, 8-diaminopelargonic acid aminotransferase and a dethiobiotin synthase, involved in the biotin biosynthesis pathway. A novel gene, BIO5 (YNR056c), is present immediately downstream from BIO4. This gene encodes Bio5p, a protein with 11 putative transmembrane regions. Uptake experiments performed with labeled 7-keto 8-aminopelargonic acid indicate that Bio5p is responsible for transport into the cell of 7-keto 8-aminopelargonic acid.

Comparison of mapping algorithms used in high-throughput sequencing: application to Ion Torrent data

BackgroundThe rapid evolution in high-throughput sequencing (HTS) technologies has opened up new perspectives in several research fields and led to the production of large volumes of sequence data. A fundamental step in HTS data analysis is the mapping of reads onto reference sequences. Choosing a suitable mapper for a given technology and a given application is a subtle task because of the difficulty of evaluating mapping algorithms.ResultsIn this paper, we present a benchmark procedure to compare mapping algorithms used in HTS using both real and simulated datasets and considering four evaluation criteria: computational resource and time requirements, robustness of mapping, ability to report positions for reads in repetitive regions, and ability to retrieve true genetic variation positions. To measure robustness, we introduced a new definition for a correctly mapped read taking into account not only the expected start position of the read but also the end position and the number of indels and substitutions. We developed CuReSim, a new read simulator, that is able to generate customized benchmark data for any kind of HTS technology by adjusting parameters to the error types. CuReSim and CuReSimEval, a tool to evaluate the mapping quality of the CuReSim simulated reads, are freely available. We applied our benchmark procedure to evaluate 14 mappers in the context of whole genome sequencing of small genomes with Ion Torrent data for which such a comparison has not yet been established.ConclusionsA benchmark procedure to compare HTS data mappers is introduced with a new definition for the mapping correctness as well as tools to generate simulated reads and evaluate mapping quality. The application of this procedure to Ion Torrent data from the whole genome sequencing of small genomes has allowed us to validate our benchmark procedure and demonstrate that it is helpful for selecting a mapper based on the intended application, questions to be addressed, and the technology used. This benchmark procedure can be used to evaluate existing or in-development mappers as well as to optimize parameters of a chosen mapper for any application and any sequencing platform.

Transcriptomic analysis in the leech Theromyzon tessulatum: involvement of cystatin B in innate immunity.

At the present time, there is little information on mechanisms of innate immunity in invertebrate groups other than insects, especially annelids. In the present study, we have performed a transcriptomic study of the immune response in the leech Theromyzon tessulatum after bacterial challenge, by a combination of differential display RT (reverse transcriptase)-PCR and cDNA microarrays. The results show relevant modulations concerning several known and unknown genes. Indeed, threonine deaminase, malate dehydrogenase, cystatin B, polyadenylate-binding protein and alpha-tubulin-like genes are up-regulated after immunostimulation. We focused on cystatin B (stefin B), which is an inhibitor of cysteine proteinases involved in the vertebrate immune response. We have cloned the full-length cDNA and named the T. tessulatum gene as Tt-cysb. Main structural features of cystatins were identified in the derived amino acid sequence of Tt-cysb cDNA; namely, a glycine residue in the N-terminus and a consensus sequence of Gln-Xaa-Val-Xaa-Gly (QXVXG) corresponding to the catalytic site. Moreover, Tt-cysb is the first cystatin B gene characterized in invertebrates. We have determined by in situ hybridization and immunocytochemistry that Tt-cysb is only expressed in large coelomic cells. In addition, this analysis confirmed that Tt-cysb is up-regulated after bacterial challenge, and that increased expression occurs only in coelomic cells. These data demonstrate that the innate immune response in the leech involves a cysteine proteinase inhibitor that is not found in ecdysozoan models, such as Drosophila melanogaster or Caenorhabditis elegans, and so underlines the great need for information about innate immunity mechanisms in different invertebrate groups.

Genetic stability of the live attenuated Bordetella pertussis vaccine candidate BPZE1

Despite the extensive use of efficacious pertussis vaccines, Bordetella pertussis infections are still among the main causes for childhood morbidity and mortality. Severe pertussis occurs mostly in very young children, often too young to be sufficiently protected by current vaccines, which require several administrations in regimens that vary between countries. Since natural infection with B. pertussis is able to induce protection, we have developed the live attenuated B. pertussis vaccine strain BPZE1 that protects mice upon a single intranasal administration. This strain was obtained by genetically inactivating pertussis toxin via two point mutations in the ptx gene, by deleting dnt encoding dermonecrotic toxin, and by replacing the B. pertussis ampG gene by Escherichia coli ampG, resulting in the removal of tracheal cytotoxin. Here, we assessed the genetic stability of BPZE1 after 20 and 27 weeks of continuous passaging in vitro and in vivo, respectively. BPZE1 was passaged 20 times in vitro and 9 times in vivo in Balb/C mice. After these passages, 8 hemolytic colonies were analyzed by PCR for the absence of dnt and B. pertussis ampG and the presence of E. coli ampG, by DNA sequencing for the presence of the two ptx point mutations and by DNA microarrays for the global genomic stability. In addition, the protective capacity of BPZE1 was evaluated after the passages. No genetic or protective difference was detected between the passaged bacteria and non-passaged BPZE1, indicating that stability of the vaccine strain is not a concern for BPZE1 to be considered as an attenuated live vaccine against whooping cough.

Evidence for induction of the ornithine transcarbamylase expression in Alzheimer's disease

To more rapidly identify candidate genes located within chromosomal regions of interest defined by genome scan studies in Alzheimers disease (AD), we have developed a customized microarray containing all the ORFs (n=2741) located within nine of these regions. Levels of gene expression were assessed in total RNA from brain tissue of 12 controls and 12 AD patients. Of all genes showing differential expression, we focused on the ornithine transcarbamylase (OTC) gene on Xp21.1., a key enzyme of the urea cycle which we found to be expressed in AD brains but not in controls, as confirmed by RT–PCR. We also detected mRNA expression of all the other urea cycle enzymes in AD brains. Immunochemistry experiments revealed that the OTC expression was strictly restricted to vascular endothelial cells in brain. Furthermore, OTC activity was 880% increased in the CSF of probable AD cases compared with controls. We analysed the association of the OTC −389 G/A and −241 A/G promoter polymorphisms with the risk of developing AD. We observed that rare haplotypes may be associated with the risk of AD through a possible modulation of the methylation of the OTC promoter. In conclusion, our results suggest the involvement of a new pathway in AD brains involving the urea cycle.

Genomic island excisions in Bordetella petrii

BackgroundAmong the members of the genus Bordetella B. petrii is unique, since it is the only species isolated from the environment, while the pathogenic Bordetellae are obligately associated with host organisms. Another feature distinguishing B. petrii from the other sequenced Bordetellae is the presence of a large number of mobile genetic elements including several large genomic regions with typical characteristics of genomic islands collectively known as integrative and conjugative elements (ICEs). These elements mainly encode accessory metabolic factors enabling this bacterium to grow on a large repertoire of aromatic compounds.ResultsDuring in vitro culture of Bordetella petrii colony variants appear frequently. We show that this variability can be attributed to the presence of a large number of metastable mobile genetic elements on its chromosome. In fact, the genome sequence of B. petrii revealed the presence of at least seven large genomic islands mostly encoding accessory metabolic functions involved in the degradation of aromatic compounds and detoxification of heavy metals. Four of these islands (termed GI1 to GI3 and GI6) are highly related to ICEclc of Pseudomonas knackmussii sp. strain B13. Here we present first data about the molecular characterization of these islands. We defined the exact borders of each island and we show that during standard culture of the bacteria these islands get excised from the chromosome. For all but one of these islands (GI5) we could detect circular intermediates. For the clc-like elements GI1 to GI3 of B. petrii we provide evidence that tandem insertion of these islands which all encode highly related integrases and attachment sites may also lead to incorporation of genomic DNA which originally was not part of the island and to the formation of huge composite islands. By integration of a tetracycline resistance cassette into GI3 we found this island to be rather unstable and to be lost from the bacterial population within about 100 consecutive generations. Furthermore, we show that GI3 is self transmissible and by conjugation can be transferred to B. bronchiseptica thus proving it to be an active integrative and conjugative elementConclusionThe results show that phenotypic variation of B. petrii is correlated with the presence of genomic islands. Tandem integration of related islands may contribute to island evolution by the acquisition of genes originally belonging to the bacterial core genome. In conclusion, B. petrii appears to be the first member of the genus in which horizontal gene transfer events have massively shaped its genome structure.

Bordetella pertussis: from functional genomics to intranasal vaccination

Whooping cough still represents a major health problem, despite the use of effective vaccines for several decades. Being classically a typical childhood disease, whooping cough in young adults is now more common than it used to be, suggesting that protection after vaccination wanes during adolescence. As an alternative to the current vaccines, we wish to develop live attenuated vaccines to be delivered by the nasal route, such as to mimic the natural route of infection and to induce long lasting immunity. Bordetella pertussis, the etiological agent of whooping cough, produces a number of virulence factors, including toxins. Its recently determined genome sequence makes it now possible to apply functional genomics, such as transcriptomics and systematic knock-out mutagenesis. The expression of most known B. pertussis virulence genes is controlled by the two-component system BvgA/S. DNA microarray analyses have led to the identification of novel genes in the BvgA/S regulon, some of which are activated by BvgA/S and others are repressed by BvgA/S. In addition, some genes appear to be differentially modulated by nicotinic acid and MgSO4, both known to modulate the expression of BvgA/S-regulated genes. Among others, the functional genomics approach has uncovered two strongly BvgA/S-activated genes, named hotA and hotB (for homolog of toxin), the products of which show high sequence similarities to pertussis toxin subunits. The identification of the full array of virulence factors, as well as an integrated understanding of the bacterial physiology should allow us to design attenuated B. pertussis strains useful for intranasal vaccination. A first generation of attenuated strains has already shown full protection in mice after a single intranasal administration. Such strains may also serve as vaccine carriers for heterologous antigens, in order to vaccinate against several different pathogens simultaneously.

Assessment of Common and Emerging Bioinformatics Pipelines for Targeted Metagenomics

Targeted metagenomics, also known as metagenetics, is a high-throughput sequencing application focusing on a nucleotide target in a microbiome to describe its taxonomic content. A wide range of bioinformatics pipelines are available to analyze sequencing outputs, and the choice of an appropriate tool is crucial and not trivial. No standard evaluation method exists for estimating the accuracy of a pipeline for targeted metagenomics analyses. This article proposes an evaluation protocol containing real and simulated targeted metagenomics datasets, and adequate metrics allowing us to study the impact of different variables on the biological interpretation of results. This protocol was used to compare six different bioinformatics pipelines in the basic user context: Three common ones (mothur, QIIME and BMP) based on a clustering-first approach and three emerging ones (Kraken, CLARK and One Codex) using an assignment-first approach. This study surprisingly reveals that the effect of sequencing errors has a bigger impact on the results that choosing different amplified regions. Moreover, increasing sequencing throughput increases richness overestimation, even more so for microbiota of high complexity. Finally, the choice of the reference database has a bigger impact on richness estimation for clustering-first pipelines, and on correct taxa identification for assignment-first pipelines. Using emerging assignment-first pipelines is a valid approach for targeted metagenomics analyses, with a quality of results comparable to popular clustering-first pipelines, even with an error-prone sequencing technology like Ion Torrent. However, those pipelines are highly sensitive to the quality of databases and their annotations, which makes clustering-first pipelines still the only reliable approach for studying microbiomes that are not well described.

Identification of novel genes potentially involved in somatic embryogenesis in chicory ( Cichorium intybus L.)

BackgroundIn our laboratory we use cultured chicory (Cichorium intybus) explants as a model to investigate cell reactivation and somatic embryogenesis and have produced 2 chicory genotypes (K59, C15) sharing a similar genetic background. K59 is a responsive genotype (embryogenic) capable of undergoing complete cell reactivation i.e. cell de- and re-differentiation leading to somatic embryogenesis (SE), whereas C15 is a non-responsive genotype (non-embryogenic) and is unable to undergo SE. Previous studies [1] showed that the use of the β-D-glucosyl Yariv reagent (β-GlcY) that specifically binds arabinogalactan-proteins (AGPs) blocked somatic embryo production in chicory root explants. This observation indicates that β-GlcY is a useful tool for investigating somatic embryogenesis (SE) in chicory. In addition, a putative AGP (DT212818) encoding gene was previously found to be significantly up-regulated in the embryogenic K59 chicory genotype as compared to the non-embryogenic C15 genotype suggesting that this AGP could be involved in chicory re-differentiation [2]. In order to improve our understanding of the molecular and cellular regulation underlying SE in chicory, we undertook a detailed cytological study of cell reactivation events in K59 and C15 genotypes, and used microarray profiling to compare gene expression in these 2 genotypes. In addition we also used β-GlcY to block SE in order to identify genes potentially involved in this process.ResultsMicroscopy confirmed that only the K59, but not the C15 genotype underwent complete cell reactivation leading to SE formation. β-GlcY-treatment of explants blocked in vitro SE induction, but not cell reactivation, and induced cell wall modifications. Microarray analyses revealed that 78 genes were differentially expressed between induced K59 and C15 genotypes. The expression profiles of 19 genes were modified by β-GlcY-treatment. Eight genes were both differentially expressed between K59 and C15 genotypes during SE induction and transcriptionally affected by β-GlcY-treatment: AGP (DT212818), 26 S proteasome AAA ATPase subunit 6 (RPT6), remorin (REM), metallothionein-1 (MT1), two non-specific lipid transfer proteins genes (SDI-9 and DEA1), 3-hydroxy-3-methylglutaryl-CoA reductase (HMG-CoA reductase), and snakin 2 (SN2). These results suggest that the 8 genes, including the previously-identified AGP gene (DT212818), could be involved in cell fate determination events leading to SE commitment in chicory.ConclusionThe use of two different chicory genotypes differing in their responsiveness to SE induction, together with β-GlcY-treatment represented an efficient tool to discriminate cell reactivation from the SE morphogenetic pathway. Such an approach, together with microarray analyses, permitted us to identify several putative key genes related to the SE morphogenetic pathway in chicory.