Christiane Bouchier

Sequencing of Aspergillus nidulans and comparative analysis with A. fumigatus and A. oryzae

The aspergilli comprise a diverse group of filamentous fungi spanning over 200 million years of evolution. Here we report the genome sequence of the model organism Aspergillus nidulans, and a comparative study with Aspergillus fumigatus, a serious human pathogen, and Aspergillus oryzae, used in the production of sake, miso and soy sauce. Our analysis of genome structure provided a quantitative evaluation of forces driving long-term eukaryotic genome evolution. It also led to an experimentally validated model of mating-type locus evolution, suggesting the potential for sexual reproduction in A. fumigatus and A. oryzae. Our analysis of sequence conservation revealed over 5,000 non-coding regions actively conserved across all three species. Within these regions, we identified potential functional elements including a previously uncharacterized TPP riboswitch and motifs suggesting regulation in filamentous fungi by Puf family genes. We further obtained comparative and experimental evidence indicating widespread translational regulation by upstream open reading frames. These results enhance our understanding of these widely studied fungi as well as provide new insight into eukaryotic genome evolution and gene regulation.

Genome evolution in yeasts

Identifying the mechanisms of eukaryotic genome evolution by comparative genomics is often complicated by the multiplicity of events that have taken place throughout the history of individual lineages, leaving only distorted and superimposed traces in the genome of each living organism. The hemiascomycete yeasts, with their compact genomes, similar lifestyle and distinct sexual and physiological properties, provide a unique opportunity to explore such mechanisms. We present here the complete, assembled genome sequences of four yeast species, selected to represent a broad evolutionary range within a single eukaryotic phylum, that after analysis proved to be molecularly as diverse as the entire phylum of chordates. A total of approximately 24,200 novel genes were identified, the translation products of which were classified together with Saccharomyces cerevisiae proteins into about 4,700 families, forming the basis for interspecific comparisons. Analysis of chromosome maps and genome redundancies reveal that the different yeast lineages have evolved through a marked interplay between several distinct molecular mechanisms, including tandem gene repeat formation, segmental duplication, a massive genome duplication and extensive gene loss.

Organised Genome Dynamics in the Escherichia coli Species Results in Highly Diverse Adaptive Paths

The Escherichia coli species represents one of the best-studied model organisms, but also encompasses a variety of commensal and pathogenic strains that diversify by high rates of genetic change. We uniformly (re-) annotated the genomes of 20 commensal and pathogenic E. coli strains and one strain of E. fergusonii (the closest E. coli related species), including seven that we sequenced to completion. Within the ∼18,000 families of orthologous genes, we found ∼2,000 common to all strains. Although recombination rates are much higher than mutation rates, we show, both theoretically and using phylogenetic inference, that this does not obscure the phylogenetic signal, which places the B2 phylogenetic group and one group D strain at the basal position. Based on this phylogeny, we inferred past evolutionary events of gain and loss of genes, identifying functional classes under opposite selection pressures. We found an important adaptive role for metabolism diversification within group B2 and Shigella strains, but identified few or no extraintestinal virulence-specific genes, which could render difficult the development of a vaccine against extraintestinal infections. Genome flux in E. coli is confined to a small number of conserved positions in the chromosome, which most often are not associated with integrases or tRNA genes. Core genes flanking some of these regions show higher rates of recombination, suggesting that a gene, once acquired by a strain, spreads within the species by homologous recombination at the flanking genes. Finally, the genomes long-scale structure of recombination indicates lower recombination rates, but not higher mutation rates, at the terminus of replication. The ensuing effect of background selection and biased gene conversion may thus explain why this region is A+T-rich and shows high sequence divergence but low sequence polymorphism. Overall, despite a very high gene flow, genes co-exist in an organised genome.

A molecular marker of artemisinin-resistant Plasmodium falciparum malaria

Plasmodium falciparum resistance to artemisinin derivatives in southeast Asia threatens malaria control and elimination activities worldwide. To monitor the spread of artemisinin resistance, a molecular marker is urgently needed. Here, using whole-genome sequencing of an artemisinin-resistant parasite line from Africa and clinical parasite isolates from Cambodia, we associate mutations in the PF3D7_1343700 kelch propeller domain (‘K13-propeller’) with artemisinin resistance in vitro and in vivo. Mutant K13-propeller alleles cluster in Cambodian provinces where resistance is prevalent, and the increasing frequency of a dominant mutant K13-propeller allele correlates with the recent spread of resistance in western Cambodia. Strong correlations between the presence of a mutant allele, in vitro parasite survival rates and in vivo parasite clearance rates indicate that K13-propeller mutations are important determinants of artemisinin resistance. K13-propeller polymorphism constitutes a useful molecular marker for large-scale surveillance efforts to contain artemisinin resistance in the Greater Mekong Subregion and prevent its global spread.

Resistance of Plasmodium falciparum field isolates to in-vitro artemether and point mutations of the SERCA-type PfATPase6

Artemisinin derivatives are an essential component of treatment against multidrug-resistant Plasmodium falciparum malaria. We aimed to investigate in-vitro resistance to artemisinin derivatives in field isolates. In-vitro susceptibility of 530 P falciparum isolates from three countries (Cambodia, French Guiana, and Senegal) with different artemisinin use was assessed with an isotopic microtest. Artemether IC50 up to 117 and 45 nmol/L was seen in French Guiana and Senegal, respectively. DNA sequencing in a subsample of 60 isolates lends support to SERCA-PfATPase6 as the target for artemisinins. The S769N PfATPase6 mutation, noted exclusively in French Guiana, was associated with raised (>30 nmol/L) artemether IC50s (p<0.0001, Mann-Whitney). All resistant isolates came from areas with uncontrolled use of artemisinin derivatives. This rise in resistance indicates the need for increased vigilance and a coordinated and rapid deployment of drug combinations.

Evolutionary dynamics of human Toll-like receptors and their different contributions to host defense.

Infectious diseases have been paramount among the threats to health and survival throughout human evolutionary history. Natural selection is therefore expected to act strongly on host defense genes, particularly on innate immunity genes whose products mediate the direct interaction between the host and the microbial environment. In insects and mammals, the Toll-like receptors (TLRs) appear to play a major role in initiating innate immune responses against microbes. In humans, however, it has been speculated that the set of TLRs could be redundant for protective immunity. We investigated how natural selection has acted upon human TLRs, as an approach to assess their level of biological redundancy. We sequenced the ten human TLRs in a panel of 158 individuals from various populations worldwide and found that the intracellular TLRs—activated by nucleic acids and particularly specialized in viral recognition—have evolved under strong purifying selection, indicating their essential non-redundant role in host survival. Conversely, the selective constraints on the TLRs expressed on the cell surface—activated by compounds other than nucleic acids—have been much more relaxed, with higher rates of damaging nonsynonymous and stop mutations tolerated, suggesting their higher redundancy. Finally, we tested whether TLRs have experienced spatially-varying selection in human populations and found that the region encompassing TLR10-TLR1-TLR6 has been the target of recent positive selection among non-Africans. Our findings indicate that the different TLRs differ in their immunological redundancy, reflecting their distinct contributions to host defense. The insights gained in this study foster new hypotheses to be tested in clinical and epidemiological genetics of infectious disease.

Genome Sequence of the Saprophyte Leptospira biflexa Provides Insights into the Evolution of Leptospira and the Pathogenesis of Leptospirosis

Leptospira biflexa is a free-living saprophytic spirochete present in aquatic environments. We determined the genome sequence of L. biflexa, making it the first saprophytic Leptospira to be sequenced. The L. biflexa genome has 3,590 protein-coding genes distributed across three circular replicons: the major 3,604 chromosome, a smaller 278-kb replicon that also carries essential genes, and a third 74-kb replicon. Comparative sequence analysis provides evidence that L. biflexa is an excellent model for the study of Leptospira evolution; we conclude that 2052 genes (61%) represent a progenitor genome that existed before divergence of pathogenic and saprophytic Leptospira species. Comparisons of the L. biflexa genome with two pathogenic Leptospira species reveal several major findings. Nearly one-third of the L. biflexa genes are absent in pathogenic Leptospira. We suggest that once incorporated into the L. biflexa genome, laterally transferred DNA undergoes minimal rearrangement due to physical restrictions imposed by high gene density and limited presence of transposable elements. In contrast, the genomes of pathogenic Leptospira species undergo frequent rearrangements, often involving recombination between insertion sequences. Identification of genes common to the two pathogenic species, L. borgpetersenii and L. interrogans, but absent in L. biflexa, is consistent with a role for these genes in pathogenesis. Differences in environmental sensing capacities of L. biflexa, L. borgpetersenii, and L. interrogans suggest a model which postulates that loss of signal transduction functions in L. borgpetersenii has impaired its survival outside a mammalian host, whereas L. interrogans has retained environmental sensory functions that facilitate disease transmission through water.

Plasmodium vivax clinical malaria is commonly observed in Duffy-negative Malagasy people

Malaria therapy, experimental, and epidemiological studies have shown that erythrocyte Duffy blood group-negative people, largely of African ancestry, are resistant to erythrocyte Plasmodium vivax infection. These findings established a paradigm that the Duffy antigen is required for P. vivax erythrocyte invasion. P. vivax is endemic in Madagascar, where admixture of Duffy-negative and Duffy-positive populations of diverse ethnic backgrounds has occurred over 2 millennia. There, we investigated susceptibility to P. vivax blood-stage infection and disease in association with Duffy blood group polymorphism. Duffy blood group genotyping identified 72% Duffy-negative individuals (FY*BES/*BES) in community surveys conducted at eight sentinel sites. Flow cytometry and adsorption–elution results confirmed the absence of Duffy antigen expression on Duffy-negative erythrocytes. P. vivax PCR positivity was observed in 8.8% (42/476) of asymptomatic Duffy-negative people. Clinical vivax malaria was identified in Duffy-negative subjects with nine P. vivax monoinfections and eight mixed Plasmodium species infections that included P. vivax (4.9 and 4.4% of 183 participants, respectively). Microscopy examination of blood smears confirmed blood-stage development of P. vivax, including gametocytes. Genotyping of polymorphic surface and microsatellite markers suggested that multiple P. vivax strains were infecting Duffy-negative people. In Madagascar, P. vivax has broken through its dependence on the Duffy antigen for establishing human blood-stage infection and disease. Further studies are necessary to identify the parasite and host molecules that enable this Duffy-independent P. vivax invasion of human erythrocytes.

Genomic analysis of smooth tubercle bacilli provides insights into ancestry and pathoadaptation of Mycobacterium tuberculosis

Global spread and limited genetic variation are hallmarks of M. tuberculosis, the agent of human tuberculosis. In contrast, Mycobacterium canettii and related tubercle bacilli that also cause human tuberculosis and exhibit unusual smooth colony morphology are restricted to East Africa. Here, we sequenced and analyzed the whole genomes of five representative strains of smooth tubercle bacilli (STB) using Sanger (4–5× coverage), 454/Roche (13–18× coverage) and/or Illumina DNA sequencing (45–105× coverage). We show that STB isolates are highly recombinogenic and evolutionarily early branching, with larger genome sizes, higher rates of genetic variation, fewer molecular scars and distinct CRISPR-Cas systems relative to M. tuberculosis. Despite the differences, all tuberculosis-causing mycobacteria share a highly conserved core genome. Mouse infection experiments showed that STB strains are less persistent and virulent than M. tuberculosis. We conclude that M. tuberculosis emerged from an ancestral STB-like pool of mycobacteria by gain of persistence and virulence mechanisms, and we provide insights into the molecular events involved.

Pfmdr1 copy number and arteminisin derivatives combination therapy failure in falciparum malaria in Cambodia

BackgroundThe combination of artesunate and mefloquine was introduced as the national first-line treatment for Plasmodium falciparum malaria in Cambodia in 2000. However, recent clinical trials performed at the Thai-Cambodian border have pointed to the declining efficacy of both artesunate-mefloquine and artemether-lumefantrine. Since pfmdr1 modulates susceptibility to mefloquine and artemisinin derivatives, the aim of this study was to assess the link between pfmdr1 copy number, in vitro susceptibility to individual drugs and treatment failure to combination therapy.MethodsBlood samples were collected from P. falciparum-infected patients enrolled in two in vivo efficacy studies in north-western Cambodia: 135 patients were treated with artemether-lumefantrine (AL group) in Sampovloun in 2002 and 2003, and 140 patients with artesunate-mefloquine (AM group) in Sampovloun and Veal Veng in 2003 and 2004. At enrollment, the in vitro IC50 was tested and the strains were genotyped for pfmdr1 copy number by real-time PCR.ResultsThe pfmdr1 copy number was analysed for 115 isolates in the AM group, and for 109 isolates in the AL group. Parasites with increased pfmdr1 copy number had significantly reduced in vitro susceptibility to mefloquine, lumefantrine and artesunate. There was no association between pfmdr1 polymorphisms and in vitro susceptibilities. In the patients treated with AM, the mean pfmdr1 copy number was lower in subjects with adequate clinical and parasitological response compared to those who experienced late treatment failure (n = 112, p < 0.001). This was not observed in the patients treated with AL (n = 96, p = 0.364). The presence of three or more copies of pfmdr1 were associated with recrudescence in artesunate-mefloquine treated patients (hazard ratio (HR) = 7.80 [95%CI: 2.09–29.10], N = 115), p = 0.002) but not with recrudescence in artemether-lumefantrine treated patients (HR = 1.03 [95%CI: 0.24–4.44], N = 109, p = 0.969).ConclusionThis study shows that pfmdr1 copy number is a molecular marker of AM treatment failure in falciparum malaria on the Thai-Cambodian border. However, while it is associated with increased IC50 for lumefantrine, pfmdr1 copy number is not associated with AL treatment failure in the area, suggesting involvement of other molecular mechanisms in AL treatment failures in Cambodia.