Complete assembly, annotation of virulence genes and CRISPR editing of the genome of Leishmania amazonensis PH8 strain

We report the sequencing and assembly of the PH8 strain of Leishmania amazonensis one of the etiological agents of leishmaniasis. After combining data from long Pacbio reads, short Illumina reads and synteny with the Leishmania mexicana genome, the sequence of 34 chromosomes with 8317 annotated genes was generated. Multigene families encoding three virulence factors, A2, amastins and the GP63 metalloproteases, were identified and compared to their annotation in other Leishmania species. As they have been recently recognized as virulence factors essential for disease establishment and progression of the infection, we also identified 14 genes encoding proteins involved in parasite iron and heme metabolism and compared to genes from other Trypanosomatids. To follow these studies with a genetic approach to address the role of virulence factors, we tested two CRISPR-Cas9 protocols to generate L. amazonensis knockout cell lines, using the Miltefosine transporter gene as a proof of concept.

Holosteans contextualize the role of the teleost genome duplication in promoting the rise of evolutionary novelties in the ray-finned fish innate immune system

Over 99% of ray-finned fishes (Actinopterygii) are teleosts, a clade that comprises half of all living vertebrate species that have diversified across virtually all fresh and saltwater ecosystems. This ecological breadth raises the question of how the immunogenetic diversity required to persist under heterogeneous pathogen pressures evolved. The teleost genome duplication (TGD) has been hypothesized as the evolutionary event that provided the substrate for rapid genomic evolution and innovation. However, studies of putative teleost-specific innate immune receptors have been largely limited to comparisons either among teleosts or between teleosts and distantly related vertebrate clades such as tetrapods. Here we describe and characterize the receptor diversity of two clustered innate immune gene families in the teleost sister lineage: Holostei (bowfin and gars). Using genomic and transcriptomic data, we provide a detailed investigation of the phylogenetic history and conserved synteny of gene clusters encoding diverse immunoglobulin domain-containing proteins (DICPs) and novel immune-type receptors (NITRs). These data demonstrate an ancient linkage of DICPs to the major histocompatibility complex (MHC) and reveal an evolutionary origin of NITR variable-joining (VJ) exons that predate the TGD by at least 50 million years. Further characterizing the receptor diversity of Holostean DICPs and NITRs illuminates a sequence diversity that rivals the diversity of these innate immune receptor families in many teleosts. Taken together, our findings provide important historical context for the evolution of these gene families that challenge prevailing expectations concerning the consequences of the TGD during actinopterygiian evolution.

The thermal proteome stability profile of Trypanosoma cruzi in epimastigote and trypomastigote life stages

Trypanosoma cruzi is a flagellate protozoa being the etiological agent of Chagas disease, a neglected tropical disease, which still poses a public health problem worldwide. The intricate molecular changes during T. cruzi-host interaction have been explored using different largescale omics techniques. However, protein stability is largely unknown. Thermal proteome profiling (TPP) methodology has the potential to characterize proteome-wide stability highlighting key proteins during T. cruzi infection and life stage transition from the invertebrate to the mammalian host. In the present work, T. cruzi epimastigotes and trypomastigotes cell lysates were subjected to TPP workflow and analyzed by quantitative large-scale mass spectrometry-based proteomics to fit a melting profile for each protein. A total of 2884 proteins were identified and associated to 1741 melting curves being 1370 in trypomastigotes (Tm_AVG 53.53 °C) and 1279 in epimastigotes (Tm_AVG 50.89 °C). A total of 453 proteins were identified with statistically different melting profiles between the two life stages. Proteins associated to pathogenesis and intracellular transport had regulated melting temperatures. Membrane and glycosylated proteins had a higher average Tm in trypomastigotes compared to epimastigotes. This study represents the first large-scale comparison of parasite protein stability between life stages. SIGNIFICANCE: Trypanosoma cruzi, a unicellular flagellate parasite, is the etiological agent of Chagas disease, endemic in South America and affecting more that 7 million people worldwide. There is an intense research to identify novel chemotherapeutic and diagnostic targets of Chagas disease. Proteomic approaches have helped in elucidating the quantitative proteome and PTMs changes of T. cruzi during life cycle transition and upon different biotic and abiotic stimuli. However, a comprehensive knowledge of the protein-protein interaction and protein conformation is still missing. In order to fill this gap, this manuscript elucidates the T. cruzi Y strain proteome-wide thermal stability map in the epimastigote and trypomastigote life stages. Comparison between life stages showed a higher average melting temperature stability for trypomastigotes than epimastigotes indicating a host temperature adaptation. Both presented a selective thermal stability shift for cellular compartments, molecular functions and biological processes based on the T. cruzi life stage. Membrane and glycosylated proteins presented a higher thermal stability in trypomastigotes when compared to the epimastigotes.

Plasmodium vinckei genomes provide insights into the pan-genome and evolution of rodent malaria parasites

Background

Rodent malaria parasites (RMPs) serve as tractable tools to study malaria parasite biology and host-parasite-vector interactions. Among the four RMPs originally collected from wild thicket rats in sub-Saharan Central Africa and adapted to laboratory mice, Plasmodium vinckei is the most geographically widespread with isolates collected from five separate locations. However, there is a lack of extensive phenotype and genotype data associated with this species, thus hindering its use in experimental studies.

Results

We have generated a comprehensive genetic resource for P. vinckei comprising of five reference-quality genomes, one for each of its subspecies, blood-stage RNA sequencing data for five P. vinckei isolates, and genotypes and growth phenotypes for ten isolates. Additionally, we sequenced seven isolates of the RMP species Plasmodium chabaudi and Plasmodium yoelii, thus extending genotypic information for four additional subspecies enabling a re-evaluation of the genotypic diversity and evolutionary history of RMPs.

The five subspecies of P. vinckei have diverged widely from their common ancestor and have undergone large-scale genome rearrangements. Comparing P. vinckei genotypes reveals region-specific selection pressures particularly on genes involved in mosquito transmission. Using phylogenetic analyses, we show that RMP multigene families have evolved differently across the vinckei and berghei groups of RMPs and that family-specific expansions in P. chabaudi and P. vinckei occurred in the common vinckei group ancestor prior to speciation. The erythrocyte membrane antigen 1 and fam-c families in particular show considerable expansions among the lowland forest-dwelling P. vinckei parasites. The subspecies from the highland forests of Katanga, P. v. vinckei, has a uniquely smaller genome, a reduced multigene family repertoire and is also amenable to transfection making it an ideal parasite for reverse genetics. We also show that P. vinckei parasites are amenable to genetic crosses.

Conclusions

Plasmodium vinckei isolates display a large degree of phenotypic and genotypic diversity and could serve as a resource to study parasite virulence and immunogenicity. Inclusion of P. vinckei genomes provide new insights into the evolution of RMPs and their multigene families. Amenability to genetic crossing and transfection make them also suitable for classical and functional genetics to study Plasmodium biology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-00995-5.

Evolutionary history of the human multigene families reveals widespread gene duplications throughout the history of animals

BACKGROUND

The hypothesis that vertebrates have experienced two ancient, whole genome duplications (WGDs) is of central interest to evolutionary biology and has been implicated in evolution of developmental complexity. Three-way and Four-way paralogy regions in human and other vertebrate genomes are considered as vital evidence to support this hypothesis. Alternatively, it has been proposed that such paralogy regions are created by small-scale duplications that occurred at different intervals over the evolution of life.

RESULTS

To address this debate, the present study investigates the evolutionary history of multigene families with at least three-fold representation on human chromosomes 1, 2, 8 and 20. Phylogenetic analysis and the tree topology comparisons classified the members of 36 multigene families into four distinct co-duplicated groups. Gene families falling within the same co-duplicated group might have duplicated together, whereas genes belong to different co-duplicated groups might have distinct evolutionary origins.

CONCLUSION

Taken together with previous investigations, the current study yielded no proof in favor of WGDs hypothesis. Rather, it appears that the vertebrate genome evolved as a result of small-scale duplication events, that cover the entire span of the animals' history.

Chromosome Mapping of H1 and H4 Histones in Parodontidae (Actinopterygii: Characiformes): Dispersed and/or Co-Opted Transposable Elements?

The karyotypes of the family Parodontidae consist of 2n = 54 chromosomes. The main chromosomal evolutionary changes of its species are attributed to chromosome rearrangements in repetitive DNA regions in their genomes. Physical mapping of the H1 and H4 histones was performed in 7 Parodontidae species to analyze the chromosome rearrangements involved in karyotype diversification in the group. In parallel, the observation of a partial sequence of an endogenous retrovirus (ERV) retrotransposon in the H1 histone sequence was evaluated to verify molecular co-option of the transposable elements (TEs) and to assess paralogous sequence dispersion in the karyotypes. Six of the studied species had an interstitial histone gene cluster in the short arm of the autosomal pair 13. Besides this interstitial cluster, in Apareiodon davisi, a probable further site was detected in the terminal region of the long arm in the same chromosome pair. The H1/H4 clusters in Parodon cf. pongoensis were located in the smallest chromosomes (pair 20). In addition, scattered H1 signals were observed on the chromosomes in all species. The H1 sequence showed an ERV in the open reading frame (ORF), and the scattered H1 signals on the chromosomes were attributed to the ERV's location. The H4 sequence had no similarity to the TEs and displayed no dispersed signals. Furthermore, the degeneration of the inner ERV in the H1 sequence (which overlapped a stretch of the H1 ORF) was discussed regarding the likelihood of molecular co-option of this retroelement in histone gene function in Parodontidae.

Guide RNA selection for CRISPR-Cas9 transfections in Plasmodium falciparum

CRISPR-Cas9 mediated genome editing is addressing key limitations in the transfection of malaria parasites. While this method has already simplified the needed molecular cloning and reduced the time required to generate mutants in the human pathogen Plasmodium falciparum, optimal selection of required guide RNAs and guidelines for successful transfections have not been well characterised, leading workers to use time-consuming trial and error approaches. We used a genome-wide computational approach to create a comprehensive and publicly accessible database of possible guide RNA sequences in the P. falciparum genome. For each guide, we report on-target efficiency and specificity scores as well as information about the genomic site relevant to optimal design of CRISPR-Cas9 transfections to modify, disrupt, or conditionally knockdown any gene. As many antimalarial drug and vaccine targets are encoded by multigene families, we also developed a new paralog specificity score that should facilitate modification of either a single family member of interest or multiple paralogs that serve overlapping roles. Finally, we tabulated features of successful transfections in our laboratory, providing broadly useful guidelines for parasite transfections. Molecular studies aimed at understanding parasite biology or characterising drug and vaccine targets in P. falciparum should be facilitated by this comprehensive database.

Retracing the path of planar cell polarity

BACKGROUND

The Planar Cell Polarity pathway (PCP) has been described as the main feature involved in patterning cell orientation in bilaterian tissues. Recently, a similar phenomenon was revealed in cnidarians, in which the inhibition of this pathway results in the absence of cilia orientation in larvae, consequently proving the functional conservation of PCP signaling between Cnidaria and Bilateria. Nevertheless, despite the growing accumulation of databases concerning basal lineages of metazoans, very few information concerning the existence of PCP components have been gathered outside of Bilateria and Cnidaria. Thus, the origin of this module or its prevalence in early emerging metazoans has yet to be elucidated.

RESULTS

The present study addresses this question by investigating the genomes and transcriptomes from all poriferan lineages in addition to Trichoplax (Placozoa) and Mnemiopsis (Ctenophora) genomes for the presence of the core components of this pathway. Our results confirm that several PCP components are metazoan innovations. In addition, we show that all members of the PCP pathway, including a bona fide Strabismus ortholog (Van gogh), are retrieved only in one sponge lineage (Homoscleromorpha) out of four. This highly suggests that the full PCP pathway dates back at least to the emergence of homoscleromorph sponges. Consequently, several secondary gene losses would have occurred in the three other poriferan lineages including Amphimedon queenslandica (Demospongiae). Several proteins were not retrieved either in placozoans or ctenophores leading us to discuss the difficulties to predict orthologous proteins in basally branching animals. Finally, we reveal how the study of multigene families may be helpful to unravel the relationships at the base of the metazoan tree.

CONCLUSION

The PCP pathway antedates the radiation of Porifera and may have arisen in the last common ancestor of animals. Oscarella species now appear as key organisms to understand the ancestral function of PCP signaling and its potential links with Wnt pathways.

High similarity of U2 snDNA sequence between A and B chromosomes in the grasshopper Abracris flavolineata

B chromosomes are frequently enriched for a wide variety of repetitive DNAs. Among grasshoppers in the species Abracris flavolineata (Ommatolampidinae) the B chromosomes are submetacentric, C-negative and harbor repetitive DNAs such as, U2 snDNA, C 0 t-1 DNA, two Mariner-like elements and some microsatellites. Here, we provide evidence showing the intragenome similarity between the B chromosome and the A complement in A. flavolineata, combining analysis of microdissection and chromosome painting and B chromosome-specific amplification through polymerase chain reaction (PCR) of U2 snDNA. Chromosome painting revealed signals spread through the C-negative regions, including the A and B chromosomes. Moreover, significant clustered signals forming bands were observed in some A chromosomes, and for the B chromosome, significant signals were located on both arms, which could be caused by accumulation of repetitive DNA sequences. The C-positive regions did not reveal any signals. Sequence comparison of U2 snDNA between that obtained from a genome without the B chromosome and that from µB-DNA revealed high similarity with the occurrence of four shared haplotypes, one of them (i.e., Hap1) being highly prevalent and putatively ancestral. The highest divergence from Hap1 was observed for Hap3, which was caused by only six mutational steps. These data support an intraspecific origin of the B chromosome in A. flavolineata that is highly similar with the A complement, and the low U2 snDNA sequence diversity observed in the B chromosome could be related to its recent origin, besides intrachromosomal concerted evolution for U2 snDNA repeats in the B chromosome.

The MHC class I genes of zebrafish

Major histocompatibility complex (MHC) molecules play a central role in the immune response and in the recognition of non-self. Found in all jawed vertebrate species, including zebrafish and other teleosts, MHC genes are considered the most polymorphic of all genes. In this review we focus on the multi-faceted diversity of zebrafish MHC class I genes, which are classified into three sequence lineages: U, Z, and L. We examine the polygenic, polymorphic, and haplotypic diversity of the zebrafish MHC class I genes, discussing known and postulated functional differences between the different class I lineages. In addition, we provide the first comprehensive nomenclature for the L lineage genes in zebrafish, encompassing at least 15 genes, and characterize their sequence properties. Finally, we discuss how recent findings have shed new light on the remarkably diverse MHC loci of this species.