Genome evolution in intracellular parasites: Microsporidia and Apicomplexa.

Microsporidia and Apicomplexa are eukaryotic, single-celled, intracellular parasites with huge public health and economic importance. Typically, these parasites are studied separately, emphasizing their uniqueness and diversity. In this review, we explore the huge amount of genomic data that has recently become available for the two groups. We compare and contrast their genome evolution and discuss how their transitions to intracellular life may have shaped it. In particular, we explore genome reduction and compaction, genome expansion and ploidy, gene shuffling and rearrangements, and the evolution of centromeres and telomeres.

The genome sequence of the Tipped Oak Case-bearer, Coleophora flavipennella (Duponchel 1843).

We present a genome assembly from an individual female Coleophora flavipennella (the Tipped Oak Case-bearer; Arthropoda; Insecta; Lepidoptera; Coleophoridae). The genome sequence is 989.3 megabases in span. Most of the assembly is scaffolded into 57 chromosomal pseudomolecules, including the W and Z sex chromosomes. The mitochondrial genome has also been assembled and is 15.77 kilobases in length.

MarkerScan: Separation and assembly of cobionts sequenced alongside target species in biodiversity genomics projects.

Contamination of public databases by mislabelled sequences has been highlighted for many years and the avalanche of novel sequencing data now being deposited has the potential to make databases difficult to use effectively. It is therefore crucial that sequencing projects and database curators perform pre-submission checks to remove obvious contamination and avoid propagating erroneous taxonomic relationships. However, it is important also to recognise that biological contamination of a target sample with unexpected species' DNA can also lead to the discovery of fascinating biological phenomena through the identification of environmental organisms or endosymbionts. Here, we present a novel, integrated method for detection and generation of high-quality genomes of all non-target genomes co-sequenced in eukaryotic genome sequencing projects. After performing taxonomic profiling of an assembly from the raw data, and leveraging the identity of small rRNA sequences discovered therein as markers, a targeted classification approach retrieves and assembles high-quality genomes. The genomes of these cobionts are then not only removed from the target species' genome but also available for further interrogation. Source code is available from https://github.com/CobiontID/MarkerScan. MarkerScan is written in Python and is deployed as a Docker container.

This article addresses a common issue in genetic research: the accidental mixing of genetic information from different species in public databases, often due to mislabelling or contamination. Interestingly, this 'contamination' can sometimes lead to exciting discoveries, like identifying DNA from unexpected species in a sample, revealing insights about organisms that live in the environment of the target organism. In our study, we developed a tool called MarkerScan for identifying these additional species found alongside the target species in eukaryotic genome sequencing projects. The method includes a way to sequence the whole genomes of the additional species. Our method involves sorting through the genetic data to identify certain small RNA sequences, which we then use as markers. These markers help to classify and assemble high-quality genomes from these additional species. This not only cleans up the main target species' genome data but also provides new, valuable genomes for further exploration.

Polyploidy is widespread in Microsporidia.

Microsporidia are obligate intracellular eukaryotic parasites with an extremely broad host range. They have both economic and public health importance. Ploidy in microsporidia is variable, with a few species formally identified as diploid and one as polyploid. Given the increase in the number of studies sequencing microsporidian genomes, it is now possible to assess ploidy levels across all currently explored microsporidian diversity. We estimate ploidy for all microsporidian data sets available on the Sequence Read Archive using k-mer-based analyses, indicating that polyploidy is widespread in Microsporidia and that ploidy change is dynamic in the group. Using genome-wide heterozygosity estimates, we also show that polyploid microsporidian genomes are relatively homozygous, and we discuss the implications of these findings on the timing of polyploidization events and their origin.IMPORTANCEMicrosporidia are single-celled intracellular parasites, distantly related to fungi, that can infect a broad range of hosts, from humans all the way to protozoans. Exploiting the wealth of microsporidian genomic data available, we use k-mer-based analyses to assess ploidy status across the group. Understanding a genome's ploidy is crucial in order to assemble it effectively and may also be relevant for better understanding a parasite's behavior and life cycle. We show that tetraploidy is present in at least six species in Microsporidia and that these polyploidization events are likely to have occurred independently. We discuss why these findings may be paradoxical, given that Microsporidia, like other intracellular parasites, have extremely small, reduced genomes.

The genome sequence of the common earthworm, Lumbricus terrestris (Linnaeus, 1758).

We present a genome assembly from an individual Lumbricus terrestris (the common earthworm; Annelida; Clitellata; Haplotaxida; Lumbricidae). The genome sequence is 1,056.5 megabases in span. Most of the assembly is scaffolded into 18 chromosomal pseudomolecules. The mitochondrial genome has also been assembled and is 15.93 kilobases in length.

Standards recommendations for the Earth BioGenome Project.

A global international initiative, such as the Earth BioGenome Project (EBP), requires both agreement and coordination on standards to ensure that the collective effort generates rapid progress toward its goals. To this end, the EBP initiated five technical standards committees comprising volunteer members from the global genomics scientific community: Sample Collection and Processing, Sequencing and Assembly, Annotation, Analysis, and IT and Informatics. The current versions of the resulting standards documents are available on the EBP website, with the recognition that opportunities, technologies, and challenges may improve or change in the future, requiring flexibility for the EBP to meet its goals. Here, we describe some highlights from the proposed standards, and areas where additional challenges will need to be met.

A dense linkage map for a large repetitive genome: discovery of the sex-determining region in hybridizing fire-bellied toads (Bombina bombina and Bombina variegata).

Genomic analysis of hybrid zones offers unique insights into emerging reproductive isolation and the dynamics of introgression. Because hybrid genomes consist of blocks inherited from one or the other parental taxon, linkage information is essential. In most cases, the spectrum of local ancestry tracts can be efficiently uncovered from dense linkage maps. Here, we report the development of such a map for the hybridizing toads, Bombina bombina and Bombina variegata (Anura: Bombinatoridae). Faced with the challenge of a large (7-10 Gb), repetitive genome, we set out to identify a large number of Mendelian markers in the nonrepetitive portion of the genome that report B. bombina vs B. variegata ancestry with appropriately quantified statistical support. Bait sequences for targeted enrichment were selected from a draft genome assembly, after filtering highly repetitive sequences. We developed a novel approach to infer the most likely diplotype per sample and locus from the raw read mapping data, which is robust to over-merging and obviates arbitrary filtering thresholds. Validation of the resulting map with 4755 markers underscored the large-scale synteny between Bombina and Xenopus tropicalis. By assessing the sex of late-stage F2 tadpoles from histological sections, we identified the sex-determining region in the Bombina genome to 7 cM on LG5, which is homologous to X. tropicalis chromosome 5, and inferred male heterogamety. Interestingly, chromosome 5 has been repeatedly recruited as a sex chromosome in anurans with XY sex determination.

Author Correction: A transcriptomic analysis of the phylum Nematoda.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Genome sequence of the root-knot nematode Meloidogyne luci.

Root-knot nematodes from the genus Meloidogyne are polyphagous plant endoparasites and agricultural pests of global importance. Here, we report the high-quality genome sequence of Meloidogyne luci population SI-Smartno V13. The resulting genome assembly of M. luci SI-Smartno V13 consists of 327 contigs, with an N50 contig length of 1,711,905 bp and a total assembly length of 209.16 Mb.Root-knot nematodes from the genus Meloidogyne are polyphagous plant endoparasites and agricultural pests of global importance. Here, we report the high-quality genome sequence of Meloidogyne luci population SI-Smartno V13. The resulting genome assembly of M. luci SI-Smartno V13 consists of 327 contigs, with an N50 contig length of 1,711,905 bp and a total assembly length of 209.16 Mb.

Chromosome-Wide Evolution and Sex Determination in the Three-Sexed Nematode Auanema rhodensis.

Trioecy, a mating system in which males, females and hermaphrodites co-exist, is a useful system to investigate the origin and maintenance of alternative mating strategies. In the trioecious nematode Auanema rhodensis, males have one X chromosome (XO), whereas females and hermaphrodites have two (XX). The female vs. hermaphrodite sex determination mechanisms have remained elusive. In this study, RNA-seq analyses show a 20% difference between the L2 hermaphrodite and female gene expression profiles. RNAi experiments targeting the DM (doublesex/mab-3) domain transcription factor dmd-10/11 suggest that the hermaphrodite sexual fate requires the upregulation of this gene. The genetic linkage map (GLM) shows that there is chromosome-wide heterozygosity for the X chromosome in F2 hermaphrodite-derived lines originated from crosses between two parental inbred strains. These results confirm the lack of recombination of the X chromosome in hermaphrodites, as previously reported. We also describe conserved chromosome elements (Nigon elements), which have been mostly maintained throughout the evolution of Rhabditina nematodes. The seven-chromosome karyotype of A. rhodensis, instead of the typical six found in other rhabditine species, derives from fusion/rearrangements events involving three Nigon elements. The A. rhodensis X chromosome is the smallest and most polymorphic with the least proportion of conserved genes. This may reflect its atypical mode of father-to-son transmission and its lack of recombination in hermaphrodites and males. In conclusion, this study provides a framework for studying the evolution of chromosomes in rhabditine nematodes, as well as possible mechanisms for the sex determination in a three-sexed species.

High genetic diversity in the Dirofilaria repens species complex revealed by mitochondrial genomes of feline microfilaria samples from Narathiwat, Thailand.

Dirofilaria repens is a zoonotic, mosquito-borne filaria infecting carnivores, particularly dogs. It is expanding its range in Europe but epidemiological information is sparse for other Eurasian regions. In Hong Kong and India, the closely related species Candidatus Dirofilaria hongkongensis was proposed. Previous analysis of 2.5 kb partial mitochondrial genome sequences containing the particularly variable non-coding control region revealed low diversity in European D. repens while Asian nematodes showed high diversity. Sequences derived from feline blood samples from Narathiwat (Thailand) led to the proposal of a third potential species, Dirofilaria sp. "Thailand II". To avoid bias from rapidly evolving non-coding regions, this study aimed to compare Dirofilaria sp. "Thailand II" with D. repens and C. D. hongkongensis based on complete mitochondrial genomes. Using PCRs and Sanger sequencing, three complete mitochondrial genomes (13,651 bp) were assembled from DNA obtained from different feline blood samples. Mitochondrial genome organization was identical to other onchocercids with eleven protein-coding, two rRNA and 22 tRNA genes and no atp-8 gene. All genes were on the same strand showing an extremely high thymidine content (56.7%). Maximum-likelihood phylogenetic analysis using protein and rRNA sequences confirmed closer relationship of Dirofilaria sp. "Thailand II" to C. D. hongkongensis than to D. repens. All distances between these three putative species were considerably larger than the distance between the valid sibling species Onchocerca volvulus and Onchocerca ochengi. Sequencing of a 2.5 kb fragment containing the control region from microfilarial DNA from additional feline blood samples from Narathiwat 3-4 years later revealed that these also fell into the C. D. hongkongensis clade but were remarkably different from C. D. hongkongensis and Dirofilaria sp. "Thailand II". Since D. repens-like filaria are absent from dogs in Narathiwat, further field studies are required to confirm if these genotypes represent locally circulating cat-specific Dirofilaria genotypes or species.

Sex- and Gamete-Specific Patterns of X Chromosome Segregation in a Trioecious Nematode.

Three key steps in meiosis allow diploid organisms to produce haploid gametes: (1) homologous chromosomes (homologs) pair and undergo crossovers; (2) homologs segregate to opposite poles; and (3) sister chromatids segregate to opposite poles. The XX/XO sex determination system found in many nematodes [1] facilitates the study of meiosis because variation is easily recognized [2-4]. Here we show that meiotic segregation of X chromosomes in the trioecious nematode Auanema rhodensis [5] varies according to sex (hermaphrodite, female, or male) and type of gametogenesis (oogenesis or spermatogenesis). In this species, XO males exclusively produce X-bearing sperm [6, 7]. The unpaired X precociously separates into sister chromatids, which co-segregate with the autosome set to generate a functional haplo-X sperm. The other set of autosomes is discarded into a residual body. Here we explore the X chromosome behavior in female and hermaphrodite meioses. Whereas X chromosomes segregate following the canonical pattern during XX female oogenesis to yield haplo-X oocytes, during XX hermaphrodite oogenesis they segregate to the first polar body to yield nullo-X oocytes. Thus, crosses between XX hermaphrodites and males yield exclusively male progeny. During hermaphrodite spermatogenesis, the sister chromatids of the X chromosomes separate during meiosis I, and homologous X chromatids segregate to the functional sperm to create diplo-X sperm. Given these intra-species, intra-individual, and intra-gametogenesis variations in the meiotic program, A. rhodensis is an ideal model for studying the plasticity of meiosis and how it can be modulated.

Genome comparisons indicate recent transfer of wRi-like Wolbachia between sister species Drosophila suzukii and D. subpulchrella.

Wolbachia endosymbionts may be acquired by horizontal transfer, by introgression through hybridization between closely related species, or by cladogenic retention during speciation. All three modes of acquisition have been demonstrated, but their relative frequency is largely unknown. Drosophila suzukii and its sister species D. subpulchrella harbor Wolbachia, denoted wSuz and wSpc, very closely related to wRi, identified in California populations of D. simulans. However, these variants differ in their induced phenotypes: wRi causes significant cytoplasmic incompatibility (CI) in D. simulans, but CI has not been detected in D. suzukii or D. subpulchrella. Our draft genomes of wSuz and wSpc contain full-length copies of 703 of the 734 single-copy genes found in wRi. Over these coding sequences, wSuz and wSpc differ by only 0.004% (i.e., 28 of 704,883 bp); they are sisters relative to wRi, from which each differs by 0.014%-0.015%. Using published data from D. melanogaster, Nasonia wasps and Nomada bees to calibrate relative rates of Wolbachia versus host nuclear divergence, we conclude that wSuz and wSpc are too similar-by at least a factor of 100-to be plausible candidates for cladogenic transmission. These three wRi-like Wolbachia, which differ in CI phenotype in their native hosts, have different numbers of orthologs of genes postulated to contribute to CI; and the CI loci differ at several nucleotides that may account for the CI difference. We discuss the general problem of distinguishing alternative modes of Wolbachia acquisition, focusing on the difficulties posed by limited knowledge of variation in absolute and relative rates of molecular evolution for host nuclear genomes, mitochondria, and Wolbachia.

Comparative Genomics of Apomictic Root-Knot Nematodes: Hybridization, Ploidy, and Dynamic Genome Change.

The root-knot nematodes (genus Meloidogyne) are important plant parasites causing substantial agricultural losses. The Meloidogyne incognita group (MIG) of species, most of which are obligatory apomicts (mitotic parthenogens), are extremely polyphagous and important problems for global agriculture. While understanding the genomic basis for their variable success on different crops could benefit future agriculture, analyses of their genomes are challenging due to complex evolutionary histories that may incorporate hybridization, ploidy changes, and chromosomal fragmentation. Here, we sequence 19 genomes, representing five species of key root-knot nematodes collected from different geographic origins. We show that a hybrid origin that predated speciation within the MIG has resulted in each species possessing two divergent genomic copies. Additionally, the apomictic MIG species are hypotriploids, with a proportion of one genome present in a second copy. The hypotriploid proportion varies among species. The evolutionary history of the MIG genomes is revealed to be very dynamic, with noncrossover recombination both homogenizing the genomic copies, and acting as a mechanism for generating divergence between species. Interestingly, the automictic MIG species M. floridensis differs from the apomict species in that it has become homozygous throughout much of its genome.

KinFin: Software for Taxon-Aware Analysis of Clustered Protein Sequences.

The field of comparative genomics is concerned with the study of similarities and differences between the information encoded in the genomes of organisms. A common approach is to define gene families by clustering protein sequences based on sequence similarity, and analyze protein cluster presence and absence in different species groups as a guide to biology. Due to the high dimensionality of these data, downstream analysis of protein clusters inferred from large numbers of species, or species with many genes, is nontrivial, and few solutions exist for transparent, reproducible, and customizable analyses. We present KinFin, a streamlined software solution capable of integrating data from common file formats and delivering aggregative annotation of protein clusters. KinFin delivers analyses based on systematic taxonomy of the species analyzed, or on user-defined, groupings of taxa, for example, sets based on attributes such as life history traits, organismal phenotypes, or competing phylogenetic hypotheses. Results are reported through graphical and detailed text output files. We illustrate the utility of the KinFin pipeline by addressing questions regarding the biology of filarial nematodes, which include parasites of veterinary and medical importance. We resolve the phylogenetic relationships between the species and explore functional annotation of proteins in clusters in key lineages and between custom taxon sets, identifying gene families of interest. KinFin can easily be integrated into existing comparative genomic workflows, and promotes transparent and reproducible analysis of clustered protein data.

Genetic Drift, Not Life History or RNAi, Determine Long-Term Evolution of Transposable Elements.

Transposable elements (TEs) are a major source of genome variation across the branches of life. Although TEs may play an adaptive role in their host's genome, they are more often deleterious, and purifying selection is an important factor controlling their genomic loads. In contrast, life history, mating system, GC content, and RNAi pathways have been suggested to account for the disparity of TE loads in different species. Previous studies of fungal, plant, and animal genomes have reported conflicting results regarding the direction in which these genomic features drive TE evolution. Many of these studies have had limited power, however, because they studied taxonomically narrow systems, comparing only a limited number of phylogenetically independent contrasts, and did not address long-term effects on TE evolution. Here, we test the long-term determinants of TE evolution by comparing 42 nematode genomes spanning over 500 million years of diversification. This analysis includes numerous transitions between life history states, and RNAi pathways, and evaluates if these forces are sufficiently persistent to affect the long-term evolution of TE loads in eukaryotic genomes. Although we demonstrate statistical power to detect selection, we find no evidence that variation in these factors influence genomic TE loads across extended periods of time. In contrast, the effects of genetic drift appear to persist and control TE variation among species. We suggest that variation in the tested factors are largely inconsequential to the large differences in TE content observed between genomes, and only by these large-scale comparisons can we distinguish long-term and persistent effects from transient or random changes.

Para-allopatry in hybridizing fire-bellied toads (Bombina bombina and B. variegata): Inference from transcriptome-wide coalescence analyses.

Ancient origins, profound ecological divergence, and extensive hybridization make the fire-bellied toads Bombina bombina and B. variegata (Anura: Bombinatoridae) an intriguing test case of ecological speciation. Previous modeling has proposed that the narrow Bombina hybrid zones represent strong barriers to neutral introgression. We test this prediction by inferring the rate of gene exchange between pure populations on either side of the intensively studied Kraków transect. We developed a method to extract high confidence sets of orthologous genes from de novo transcriptome assemblies, fitted a range of divergence models to these data and assessed their relative support with analytic likelihood calculations. There was clear evidence for postdivergence gene flow, but, as expected, no perceptible signal of recent introgression via the nearby hybrid zone. The analysis of two additional Bombina taxa (B. v. scabra and B. orientalis) validated our parameter estimates against a larger set of prior expectations. Despite substantial cumulative introgression over millions of years, adaptive divergence of the hybridizing taxa is essentially unaffected by their lack of reproductive isolation. Extended distribution ranges also buffer them against small-scale environmental perturbations that have been shown to reverse the speciation process in other, more recent ecotypes.

Stage-specific Proteomes from Onchocerca ochengi, Sister Species of the Human River Blindness Parasite, Uncover Adaptations to a Nodular Lifestyle.

Despite 40 years of control efforts, onchocerciasis (river blindness) remains one of the most important neglected tropical diseases, with 17 million people affected. The etiological agent, Onchocerca volvulus, is a filarial nematode with a complex lifecycle involving several distinct stages in the definitive host and blackfly vector. The challenges of obtaining sufficient material have prevented high-throughput studies and the development of novel strategies for disease control and diagnosis. Here, we utilize the closest relative of O. volvulus, the bovine parasite Onchocerca ochengi, to compare stage-specific proteomes and host-parasite interactions within the secretome. We identified a total of 4260 unique O. ochengi proteins from adult males and females, infective larvae, intrauterine microfilariae, and fluid from intradermal nodules. In addition, 135 proteins were detected from the obligate Wolbachia symbiont. Observed protein families that were enriched in all whole body extracts relative to the complete search database included immunoglobulin-domain proteins, whereas redox and detoxification enzymes and proteins involved in intracellular transport displayed stage-specific overrepresentation. Unexpectedly, the larval stages exhibited enrichment for several mitochondrial-related protein families, including members of peptidase family M16 and proteins which mediate mitochondrial fission and fusion. Quantification of proteins across the lifecycle using the Hi-3 approach supported these qualitative analyses. In nodule fluid, we identified 94 O. ochengi secreted proteins, including homologs of transforming growth factor-ß and a second member of a novel 6-ShK toxin domain family, which was originally described from a model filarial nematode (Litomosoides sigmodontis). Strikingly, the 498 bovine proteins identified in nodule fluid were strongly dominated by antimicrobial proteins, especially cathelicidins. This first high-throughput analysis of an Onchocerca spp. proteome across the lifecycle highlights its profound complexity and emphasizes the extremely close relationship between O. ochengi and O. volvulus The insights presented here provide new candidates for vaccine development, drug targeting and diagnostic biomarkers.