Comparative genomics of Giardia duodenalis sub-assemblage AI beaver (Be-2) and human (WB-C6) strains show remarkable homozygosity, sequence similarity, and conservation of VSP genes.

Giardia duodenalis, a major cause of waterborne infection, infects a wide range of mammalian hosts and is subdivided into eight genetically well-defined assemblages named A through H. However, fragmented genomes and a lack of comparative analysis within and between the assemblages render unclear the molecular mechanisms controlling host specificity and differential disease outcomes. To address this, we generated a near-complete de novo genome of AI assemblage using the Oxford Nanopore platform by sequencing the Be-2 genome. We generated 148,144 long-reads with quality scores of > 7. The final genome assembly consists of only nine contigs with an N50 of 3,045,186 bp. This assembly agrees closely with the assembly of another strain in the AI assemblage (WB-C6). However, a critical difference is that a region previously placed in the five-prime region of Chr5 belongs to Chr4 of Be-2. We find a high degree of conservation in the ploidy, homozygosity, and the presence of cysteine-rich variant-specific surface proteins (VSPs) within the AI assemblage. Our assembly provides a nearly complete genome of a member of the AI assemblage of G. duodenalis, aiding population genomic studies capable of elucidating Giardia transmission, host range, and pathogenicity.

Integrative analyses on the ciliates Colpoda illuminate the life history evolution of soil microorganisms.

Microorganisms play a central role in sustaining soil ecosystems and agriculture, and these functions are usually associated with their complex life history. Yet, the regulation and evolution of life history have remained enigmatic and poorly understood, especially in protozoa, the third most abundant group of organisms in the soil. Here, we explore the life history of a cosmopolitan species-Colpoda steinii. Our analysis has yielded a high-quality macronuclear genome for C. steinii, with size of 155 Mbp and 37,123 protein-coding genes, as well as mean intron length of ~93 bp, longer than most other studied ciliates. Notably, we identify two possible whole-genome duplication events in C. steinii, which may account for its genome being about twice the size of C. inflata's, another co-existing species. We further resolve the gene expression profiles in diverse life stages of C. steinii, which are also corroborated in C. inflata. During the resting cyst stage, genes associated with cell death and vacuole formation are upregulated, and translation-related genes are downregulated. While the translation-related genes are upregulated during the excystment of resting cysts. Reproductive cysts exhibit a significant reduction in cell adhesion. We also demonstrate that most genes expressed in specific life stages are under strong purifying selection. This study offers a deeper understanding of the life history evolution that underpins the extraordinary success and ecological functions of microorganisms in soil ecosystems.IMPORTANCEColpoda species, as a prominent group among the most widely distributed and abundant soil microorganisms, play a crucial role in sustaining soil ecosystems and promoting plant growth. This investigation reveals their exceptional macronuclear genomic features, including significantly large genome size, long introns, and numerous gene duplications. The gene expression profiles and the specific biological functions associated with the transitions between various life stages are also elucidated. The vast majority of genes linked to life stage transitions are subject to strong purifying selection, as inferred from multiple natural strains newly isolated and deeply sequenced. This substantiates the enduring and conservative nature of Colpoda's life history, which has persisted throughout the extensive evolutionary history of these highly successful protozoa in soil. These findings shed light on the evolutionary dynamics of microbial eukaryotes in the ever-fluctuating soil environments. This integrative research represents a significant advancement in understanding the life histories of these understudied single-celled eukaryotes.

Far-East Asian Toxoplasma isolates share ancestry with North and South/Central American recombinant lineages.

Toxoplasma gondii is a global protozoan pathogen. Clonal lineages predominate in Europe, North America, Africa, and China, whereas highly recombinant parasites are endemic in South/Central America. Far East Asian T. gondii isolates are not included in current global population genetic structure analyses at WGS resolution. Here we report a genome-wide population study that compared eight Japanese and two Chinese isolates against representative worldwide T. gondii genomes using POPSICLE, a novel population structure analyzing software. Also included were 7 genomes resurrected from non-viable isolates by target enrichment sequencing. Visualization of the genome structure by POPSICLE shows a mixture of Chinese haplogroup (HG) 13 haploblocks introgressed within the genomes of Japanese HG2 and North American HG12. Furthermore, two ancestral lineages were identified in the Japanese strains; one lineage shares a common ancestor with HG11 found in both Japanese strains and North American HG12. The other ancestral lineage, found in T. gondii isolates from a small island in Japan, is admixed with genetically diversified South/Central American strains. Taken together, this study suggests multiple ancestral links between Far East Asian and American T. gondii strains and provides insight into the transmission history of this cosmopolitan organism.

Genome-wide distribution of 5-hydroxymethyluracil and chromatin accessibility in the Breviolum minutum genome.

BACKGROUND: In dinoflagellates, a unique and extremely divergent genomic and nuclear organization has evolved. The highly unusual features of dinoflagellate nuclei and genomes include permanently condensed liquid crystalline chromosomes, primarily packaged by proteins other than histones, genes organized in very long unidirectional gene arrays, a general absence of transcriptional regulation, high abundance of the otherwise very rare DNA modification 5-hydroxymethyluracil (5-hmU), and many others. While most of these fascinating properties are originally identified in the 1970s and 1980s, they have not yet been investigated using modern genomic tools. RESULTS: In this work, we address some of the outstanding questions regarding dinoflagellate genome organization by mapping the genome-wide distribution of 5-hmU (using both immunoprecipitation-based and basepair-resolution chemical mapping approaches) and of chromatin accessibility in the genome of the Symbiodiniaceae dinoflagellate Breviolum minutum. We find that the 5-hmU modification is preferentially enriched over certain classes of repetitive elements, often coincides with the boundaries between gene arrays, and is generally correlated with decreased chromatin accessibility, the latter otherwise being largely uniform along the genome. We discuss the potential roles of 5-hmU in the functional organization of dinoflagellate genomes and its relationship to the transcriptional landscape of gene arrays. CONCLUSIONS: Our results provide the first window into the 5-hmU and chromatin accessibility landscapes in dinoflagellates.

Genomic insights into the cellular specialization of predation in raptorial protists.

BACKGROUND: Predation is a fundamental mechanism for organisms to acquire energy, and various species have evolved diverse tools to enhance their hunting abilities. Among protozoan predators, raptorial Haptorian ciliates are particularly fascinating as they possess offensive extrusomes known as toxicysts, which are rapidly discharged upon prey contact. However, our understanding of the genetic processes and specific toxins involved in toxicyst formation and discharge is still limited. RESULTS: In this study, we investigated the predation strategies and subcellular structures of seven Haptoria ciliate species and obtained their genome sequences using single-cell sequencing technology. Comparative genomic analysis revealed distinct gene duplications related to membrane transport proteins and hydrolytic enzymes in Haptoria, which play a crucial role in the production and discharge of toxicysts. Transcriptomic analysis further confirmed the abundant expression of genes related to membrane transporters and cellular toxins in Haptoria compared to Trichostomatia. Notably, polyketide synthases (PKS) and L-amino acid oxidases (LAAO) were identified as potentially toxin genes that underwent extensive duplication events in Haptoria. CONCLUSIONS: Our results shed light on the evolutionary and genomic adaptations of Haptorian ciliates for their predation strategies in evolution and provide insights into their toxic mechanisms.

Whole-genome assembly of a hybrid Trypanosoma cruzi strain assembled with Nanopore sequencing alone.

Trypanosoma cruzi is the causative agent of Chagas disease, which causes 10,000 deaths per year. Despite the high mortality associated with Chagas, relatively few parasite genomes have been assembled to date, with genome assemblies unavailable even for some commonly used laboratory strains. This is at least partially due to T. cruzi's highly complex and highly repetitive genome, which defies investigation using traditional short-read sequencing methods. In this study, we have generated a high-quality whole-genome assembly of the hybrid Tulahuen strain, a commercially available type VI strain, using long-read Nanopore sequencing without short-read scaffolding. The assembled genome contains 25% repeat regions, 17% variable multigene family members, and 27% transposable elements (TEs) and is of comparable quality with T. cruzi genome assemblies that utilized both long- and short-read data. Notably, we find that regions with TEs are significantly enriched for multicopy surface proteins, and that surface proteins are, on average, closer to TEs than to other coding regions. This finding suggests that mobile genetic elements such as transposons may drive recombination within surface protein gene families. This work demonstrates the feasibility of Nanopore sequencing to resolve complex regions of T. cruzi genomes, and with these resolved regions, provides support for a possible mechanism for genomic diversification.

Massive genome reduction predates the divergence of Symbiodiniaceae dinoflagellates.

Dinoflagellates in the family Symbiodiniaceae are taxonomically diverse, predominantly symbiotic lineages that are well-known for their association with corals. The ancestor of these taxa is believed to have been free-living. The establishment of symbiosis (i.e. symbiogenesis) is hypothesized to have occurred multiple times during Symbiodiniaceae evolution, but its impact on genome evolution of these taxa is largely unknown. Among Symbiodiniaceae, the genus Effrenium is a free-living lineage that is phylogenetically positioned between two robustly supported groups of genera within which symbiotic taxa have emerged. The apparent lack of symbiogenesis in Effrenium suggests that the ancestral features of Symbiodiniaceae may have been retained in this lineage. Here, we present de novo assembled genomes (1.2-1.9 Gbp in size) and transcriptome data from three isolates of Effrenium voratum and conduct a comparative analysis that includes 16 Symbiodiniaceae taxa and the other dinoflagellates. Surprisingly, we find that genome reduction, which is often associated with a symbiotic lifestyle, predates the origin of Symbiodiniaceae. The free-living lifestyle distinguishes Effrenium from symbiotic Symbiodiniaceae vis-à-vis their longer introns, more-extensive mRNA editing, fewer (~30%) lineage-specific gene sets, and lower (~10%) level of pseudogenization. These results demonstrate how genome reduction and the adaptation to distinct lifestyles intersect to drive diversification and genome evolution of Symbiodiniaceae.

Ancestral aneuploidy and stable chromosomal duplication resulting in differential genome structure and gene expression control in trypanosomatid parasites.

Aneuploidy is widely observed in both unicellular and multicellular eukaryotes, usually associated with adaptation to stress conditions. Chromosomal duplication stability is a tradeoff between the fitness cost of having unbalanced gene copies and the potential fitness gained from increased dosage of specific advantageous genes. Trypanosomatids, a family of protozoans that include species that cause neglected tropical diseases, are a relevant group to study aneuploidies. Their life cycle has several stressors that could select for different patterns of chromosomal duplications and/or losses, and their nearly universal use of polycistronic transcription increases their reliance on gene expansion/contraction, as well as post-transcriptional control as mechanisms for gene expression regulation. By evaluating the data from 866 isolates covering seven trypanosomatid genera, we have revealed that aneuploidy tolerance is an ancestral characteristic of trypanosomatids but has a reduced occurrence in a specific monophyletic clade that has undergone large genomic reorganization and chromosomal fusions. We have also identified an ancient chromosomal duplication that was maintained across these parasite's speciation, named collectively as the trypanosomatid ancestral supernumerary chromosome (TASC). TASC has most genes in the same coding strand, is expressed as a disomic chromosome (even having four copies), and has increased potential for functional variation, but it purges highly deleterious mutations more efficiently than other chromosomes. The evidence of stringent control over gene expression in this chromosome suggests that these parasites have adapted to mitigate the fitness cost associated with this ancient chromosomal duplication.

Comparative genomics analysis reveals sequence characteristics potentially related to host preference in Cryptosporidium xiaoi.

Cryptosporidium spp. are important diarrhea-associated pathogens in humans and livestock. Among the known species, Cryptosporidium xiaoi, which causes cryptosporidiosis in sheep and goats, was previously recognized as a genotype of the bovine-specific Cryptosporidium bovis based on their high sequence identity in the ssrRNA gene. However, the lack of genomic data has limited characterization of the genetic differences between the two closely related species. In this study, we sequenced the genomes of two C. xiaoi isolates and performed comparative genomic analysis to identify the sequence uniqueness of this ovine-adapted species compared with other Cryptosporidium spp. Our results showed that C. xiaoi is genetically related to C. bovis as shown by their 95.8% genomic identity and similar gene content. Consistent with this, both C. xiaoi and C. bovis appear to have fewer genes encoding mitochondrial metabolic enzymes and invasion-related protein families. However, they appear to possess several species-specific genes. Further analysis indicates that the sequence differences between these two Cryptosporidium spp. are mainly in 24 highly polymorphic genes, half of which are located in the subtelomeric regions. Some of these subtelomeric genes encode secretory proteins that have undergone positive selection. In addition, the genomes of two C. xiaoi isolates, identified as subtypes XXIIIf and XXIIIh, share 99.9% nucleotide sequence identity, with six highly divergent genes encoding putative secretory proteins. Therefore, these species-specific genes and sequence polymorphism in subtelomeric genes probably contribute to the different host preference of C. xiaoi and C. bovis.

PbARID-associated chromatin remodeling events are essential for gametocyte development in Plasmodium.

Gametocyte development of the Plasmodium parasite is a key step for transmission of the parasite. Male and female gametocytes are produced from a subpopulation of asexual blood-stage parasites, but the mechanisms that regulate the differentiation of sexual stages are still under investigation. In this study, we investigated the role of PbARID, a putative subunit of a SWI/SNF chromatin remodeling complex, in transcriptional regulation during the gametocyte development of P. berghei. PbARID expression starts in early gametocytes before the manifestation of male and female-specific features, and disruption of its gene results in the complete loss of gametocytes with detectable male features and the production of abnormal female gametocytes. ChIP-seq analysis of PbARID showed that it forms a complex with gSNF2, an ATPase subunit of the SWI/SNF chromatin remodeling complex, associating with the male cis-regulatory element, TGTCT. Further ChIP-seq of PbARID in gsnf2-knockout parasites revealed an association of PbARID with another cis-regulatory element, TGCACA. RIME and DNA-binding assays suggested that HDP1 is the transcription factor that recruits PbARID to the TGCACA motif. Our results indicated that PbARID could function in two chromatin remodeling events and paly essential roles in both male and female gametocyte development.

Genome assembly of a symbiotic balantidia (Balantidium ctenopharyngodoni) in fish hindgut.

Balantidium ctenopharyngodoni is identified as the sole ciliate species that exclusively resides within the hindgut of grass carp with high prevalence and intensity. In this study, the successful cultivation of B. ctenopharyngodoni enabled us to collect enough cells for genome sequencing. Consequently, we acquired a high-quality genome assembly spanning 68.66 Mb, encompassing a total of 22,334 nanochromosomes. Furthermore, we predicted 29,348 protein-coding genes, and 95.5% of them was supported by the RNA-seq data. The trend of GC content in the subtelomeric regions of single-gene chromosomes was similar to other ciliates containing nanochromosomes. A large number of genes encoding carbohydrate-binding modules with affinities for starch and peptidoglycans was identified. The identification of mitochondrion-related organelles (MROs) within genome indicates its well-suited adaptation to the anaerobic conditions in the hindgut environment. In summary, our results will offer resources for understanding the genetic basis and molecular adaptations of balantidia to hindgut of herbivorous fish.

Uncoupling programmed DNA cleavage and repair scrambles the Paramecium somatic genome.

In the ciliate Paramecium, precise excision of numerous internal eliminated sequences (IESs) from the somatic genome is essential at each sexual cycle. DNA double-strands breaks (DSBs) introduced by the PiggyMac endonuclease are repaired in a highly concerted manner by the non-homologous end joining (NHEJ) pathway, illustrated by complete inhibition of DNA cleavage when Ku70/80 proteins are missing. We show that expression of a DNA-binding-deficient Ku70 mutant (Ku70-6E) permits DNA cleavage but leads to the accumulation of unrepaired DSBs. We uncoupled DNA cleavage and repair by co-expressing wild-type and mutant Ku70. High-throughput sequencing of the developing macronucleus genome in these conditions identifies the presence of extremities healed by de novo telomere addition and numerous translocations between IES-flanking sequences. Coupling the two steps of IES excision ensures that both extremities are held together throughout the process, suggesting that DSB repair proteins are essential for assembly of a synaptic precleavage complex.

Conspicuous chloroplast with light harvesting-photosystem I/II megacomplex in marine Prorocentrum cordatum.

Marine photosynthetic (micro)organisms drive multiple biogeochemical cycles and display a large diversity. Among them, the bloom-forming, free-living dinoflagellate Prorocentrum cordatum CCMP 1329 (formerly P. minimum) stands out with its distinct cell biological features. Here, we obtained insights into the structural properties of the chloroplast and the photosynthetic machinery of P. cordatum using microscopic and proteogenomic approaches. High-resolution FIB/SEM analysis revealed a single large chloroplast (∼40% of total cell volume) with a continuous barrel-like structure, completely lining the inner face of the cell envelope and enclosing a single reticular mitochondrium, the Golgi apparatus, as well as diverse storage inclusions. Enriched thylakoid membrane fractions of P. cordatum were comparatively analyzed with those of the well-studied model-species Arabidopsis (Arabidopsis thaliana) using 2D BN DIGE. Strikingly, P. cordatum possessed a large photosystem-light harvesting megacomplex (>1.5 MDa), which is dominated by photosystems I and II (PSI, PSII), chloroplast complex I, and chlorophyll a-b binding light harvesting complex proteins. This finding parallels the absence of grana in its chloroplast and distinguishes from the predominant separation of PSI and PSII complexes in A. thaliana, indicating a different mode of flux balancing. Except for the core elements of the ATP synthase and the cytb6f-complex, the composition of the other complexes (PSI, PSII, and pigment-binding proteins, PBPs) of P. cordatum differed markedly from those of A. thaliana. Furthermore, a high number of PBPs was detected, accounting for a large share of the total proteomic data (∼65%) and potentially providing P. cordatum with flexible adaptation to changing light regimes.

Recent Advances in CRISPR/Cas9-Mediated Genome Editing in Leishmania Strains.

BACKGROUND: Genome manipulation of Leishmania species and the creation of modified strains are widely employed strategies for various purposes, including gene function studies, the development of live attenuated vaccines, and the engineering of host cells for protein production. OBJECTIVE: Despite the introduction of novel manipulation approaches like CRISPR/Cas9 technology with significant advancements in recent years, the development of a reliable protocol for efficiently and precisely altering the genes of Leishmania strains remains a challenging endeavor. Following the successful adaptation of the CRISPR/Cas9 system for higher eukaryotic cells, several research groups have endeavored to apply this system to manipulate the genome of Leishmania. RESULTS: Despite the substantial differences between Leishmania and higher eukaryotes, the CRISPR/Cas9 system has been effectively tested and applied in Leishmania.  CONCLUSION: This comprehensive review summarizes all the CRISPR/Cas9 systems that have been employed in Leishmania, providing details on their methods and the expression systems for Cas9 and gRNA. The review also explores the various applications of the CRISPR system in Leishmania, including the deletion of multicopy gene families, the development of the Leishmania vaccine, complete gene deletions, investigations into chromosomal translocations, protein tagging, gene replacement, large-scale gene knockout, genome editing through cytosine base replacement, and its innovative use in the detection of Leishmania. In addition, the review offers an up-to-date overview of all double-strand break repair mechanisms in Leishmania.

A nuclear orthologue of the dNTP triphosphohydrolase SAMHD1 controls dNTP homeostasis and genomic stability in Trypanosoma brucei.

Maintenance of dNTPs pools in Trypanosoma brucei is dependent on both biosynthetic and degradation pathways that together ensure correct cellular homeostasis throughout the cell cycle which is essential for the preservation of genomic stability. Both the salvage and de novo pathways participate in the provision of pyrimidine dNTPs while purine dNTPs are made available solely through salvage. In order to identify enzymes involved in degradation here we have characterized the role of a trypanosomal SAMHD1 orthologue denominated TbHD82. Our results show that TbHD82 is a nuclear enzyme in both procyclic and bloodstream forms of T. brucei. Knockout forms exhibit a hypermutator phenotype, cell cycle perturbations and an activation of the DNA repair response. Furthermore, dNTP quantification of TbHD82 null mutant cells revealed perturbations in nucleotide metabolism with a substantial accumulation of dATP, dCTP and dTTP. We propose that this HD domain-containing protein present in kinetoplastids plays an essential role acting as a sentinel of genomic fidelity by modulating the unnecessary and detrimental accumulation of dNTPs.

Selective whole-genome amplification reveals population genetics of Leishmania braziliensis directly from patient skin biopsies.

In Brazil, Leishmania braziliensis is the main causative agent of the neglected tropical disease, cutaneous leishmaniasis (CL). CL presents on a spectrum of disease severity with a high rate of treatment failure. Yet the parasite factors that contribute to disease presentation and treatment outcome are not well understood, in part because successfully isolating and culturing parasites from patient lesions remains a major technical challenge. Here we describe the development of selective whole genome amplification (SWGA) for Leishmania and show that this method enables culture-independent analysis of parasite genomes obtained directly from primary patient skin samples, allowing us to circumvent artifacts associated with adaptation to culture. We show that SWGA can be applied to multiple Leishmania species residing in different host species, suggesting that this method is broadly useful in both experimental infection models and clinical studies. SWGA carried out directly on skin biopsies collected from patients in Corte de Pedra, Bahia, Brazil, showed extensive genomic diversity. Finally, as a proof-of-concept, we demonstrated that SWGA data can be integrated with published whole genome data from cultured parasite isolates to identify variants unique to specific geographic regions in Brazil where treatment failure rates are known to be high. SWGA provides a relatively simple method to generate Leishmania genomes directly from patient samples, unlocking the potential to link parasite genetics with host clinical phenotypes.

Accessing the Variability of Multicopy Genes in Complex Genomes using Unassembled Next-Generation Sequencing Reads: The Case of Trypanosoma cruzi Multigene Families.

Repetitive elements cause assembly fragmentation in complex eukaryotic genomes, limiting the study of their variability. The genome of Trypanosoma cruzi, the parasite that causes Chagas disease, has a high repetitive content, including multigene families. Although many T. cruzi multigene families encode surface proteins that play pivotal roles in host-parasite interactions, their variability is currently underestimated, as their high repetitive content results in collapsed gene variants. To estimate sequence variability and copy number variation of multigene families, we developed a read-based approach that is independent of gene-specific read mapping and de novo assembly. This methodology was used to estimate the copy number and variability of MASP, TcMUC, and Trans-Sialidase (TS), the three largest T. cruzi multigene families, in 36 strains, including members of all six parasite discrete typing units (DTUs). We found that these three families present a specific pattern of variability and copy number among the distinct parasite DTUs. Inter-DTU hybrid strains presented a higher variability of these families, suggesting that maintaining a larger content of their members could be advantageous. In addition, in a chronic murine model and chronic Chagasic human patients, the immune response was focused on TS antigens, suggesting that targeting TS conserved sequences could be a potential avenue to improve diagnosis and vaccine design against Chagas disease. Finally, the proposed approach can be applied to study multicopy genes in any organism, opening new avenues to access sequence variability in complex genomes. IMPORTANCE Sequences that have several copies in a genome, such as multicopy-gene families, mobile elements, and microsatellites, are among the most challenging genomic segments to study. They are frequently underestimated in genome assemblies, hampering the correct assessment of these important players in genome evolution and adaptation. Here, we developed a new methodology to estimate variability and copy numbers of repetitive genomic regions and employed it to characterize the T. cruzi multigene families MASP, TcMUC, and transsialidase (TS), which are important virulence factors in this parasite. We showed that multigene families vary in sequence and content among the parasite's lineages, whereas hybrid strains have a higher sequence variability that could be advantageous to the parasite's survivability. By identifying conserved sequences within multigene families, we showed that the mammalian host immune response toward these multigene families is usually focused on the TS multigene family. These TS conserved and immunogenic peptides can be explored in future works as diagnostic targets or vaccine candidates for Chagas disease. Finally, this methodology can be easily applied to any organism of interest, which will aid in our understanding of complex genomic regions.

The in vivo RNA structurome of the malaria parasite Plasmodium falciparum, a protozoan with an A/U-rich transcriptome.

Plasmodium falciparum, a protozoan parasite and causative agent of human malaria, has one of the most A/T-biased genomes sequenced to date. This may give the genome and the transcriptome unusual structural features. Recent progress in sequencing techniques has made it possible to study the secondary structures of RNA molecules at the transcriptomic level. Thus, in this study we produced the in vivo RNA structurome of a protozoan parasite with a highly A/U-biased transcriptome. We showed that it is possible to probe the secondary structures of P. falciparum RNA molecules in vivo using two different chemical probes, and obtained structures for more than half of all transcripts in the transcriptome. These showed greater stability (lower free energy) than the same structures modelled in silico, and structural features appeared to influence translation efficiency and RNA decay. Finally, we compared the P. falciparum RNA structurome with the predicted RNA structurome of an A/U-balanced species, P. knowlesi, finding a bias towards lower overall transcript stability and more hairpins and multi-stem loops in P. falciparum. This unusual protozoan RNA structurome will provide a basis for similar studies in other protozoans and also in other unusual genomes.