Phased chromosome-scale genome assembly of an asexual, allopolyploid root-knot nematode reveals complex subgenomic structure.

We present the chromosome-scale genome assembly of the allopolyploid root-knot nematode Meloidogyne javanica. We show that the M. javanica genome is predominantly allotetraploid, comprising two subgenomes, A and B, that most likely originated from hybridisation of two ancestral parental species. The assembly was annotated using full-length non-chimeric transcripts, comparison to reference databases, and ab initio prediction techniques, and the subgenomes were phased using ancestral k-mer spectral analysis. Subgenome B appears to show fission of chromosomal contigs, and while there is substantial synteny between subgenomes, we also identified regions lacking synteny that may have diverged in the ancestral genomes prior to or following hybridisation. This annotated and phased genome assembly forms a significant resource for understanding the origins and genetics of these globally important plant pathogens.

Identification of individual root-knot nematodes using low coverage long-read sequencing.

Root-knot nematodes (RKN; genus Meloidogyne) are polyphagous plant pathogens of great economic importance to agriculturalists globally. These species are small, diverse, and can be challenging for accurate taxonomic identification. Many of the most important crop pests confound analysis with simple genetic marker loci as they are polyploids of likely hybrid origin. Here we take a low-coverage, long-read genome sequencing approach to characterisation of individual root-knot nematodes. We demonstrate library preparation for Oxford Nanopore Technologies Flongle sequencing of low input DNA from individual juveniles and immature females, multiplexing up to twelve samples per flow cell. Taxonomic identification with Kraken 2 (a k-mer-based taxonomic assignment tool) is shown to reliably identify individual nematodes to species level, even within the very closely related Meloidogyne incognita group. Our approach forms a robust, low-cost, and scalable method for accurate RKN species diagnostics.

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.

The Genomic Substrate for Adaptive Radiation: Copy Number Variation across 12 Tribes of African Cichlid Species.

The initial sequencing of five cichlid genomes revealed an accumulation of genetic variation, including extensive copy number variation in cichlid lineages particularly those that have undergone dramatic evolutionary radiation. Gene duplication has the potential to generate substantial molecular substrate for the origin of evolutionary novelty. We use array-based comparative heterologous genomic hybridization to identify copy number variation events (CNVEs) for 168 samples representing 53 cichlid species including the 5 species for which full genome sequence is available. We identify an average of 50-100 CNVEs per individual. For those species represented by multiple samples, we identify 150-200 total CNVEs suggesting a substantial amount of intraspecific variation. For these species, only ∼10% of the detected CNVEs are fixed. Hierarchical clustering of species according to CNVE data recapitulates phylogenetic relationships fairly well at both the tribe and radiation level. Although CNVEs are detected on all linkage groups, they tend to cluster in "hotspots" and are likely to contain and be flanked by transposable elements. Furthermore, we show that CNVEs impact functional categories of genes with potential roles in adaptive phenotypes that could reasonably promote divergence and speciation in the cichlid clade. These data contribute to a more complete understanding of the molecular basis for adaptive natural selection, speciation, and evolutionary radiation.

Detecting host-parasitoid interactions in an invasive Lepidopteran using nested tagging DNA metabarcoding.

Determining the host-parasitoid interactions and parasitism rates for invasive species entering novel environments is an important first step in assessing potential routes for biocontrol and integrated pest management. Conventional insect rearing techniques followed by taxonomic identification are widely used to obtain such data, but this can be time-consuming and prone to biases. Here, we present a next-generation sequencing approach for use in ecological studies which allows for individual-level metadata tracking of large numbers of invertebrate samples through the use of hierarchically organised molecular identification tags. We demonstrate its utility using a sample data set examining both species identity and levels of parasitism in late larval stages of the oak processionary moth (Thaumetopoea processionea-Linn. 1758), an invasive species recently established in the United Kingdom. Overall, we find that there are two main species exploiting the late larval stages of oak processionary moth in the United Kingdom with the main parasitoid (Carcelia iliaca-Ratzeburg, 1840) parasitising 45.7% of caterpillars, while a rare secondary parasitoid (Compsilura concinnata-Meigen, 1824) was also detected in 0.4% of caterpillars. Using this approach on all life stages of the oak processionary moth may demonstrate additional parasitoid diversity. We discuss the wider potential of nested tagging DNA metabarcoding for constructing large, highly resolved species interaction networks.

Climate warming alters the structure of farmland tritrophic ecological networks and reduces crop yield.

It is unclear how sustained increases in temperature and changes in precipitation, as a result of climate change, will affect crops and their interactions with agricultural weeds, insect pests and predators, due to the difficulties in quantifying changes in such complex relationships. We simulated the combined effects of increasing temperature (by an average of 1.4°C over a growing season) and applying additional rainwater (10% of the monthly mean added weekly, 40% total) using a replicated, randomized block experiment within a wheat crop. We examined how this affected the structure of 24 quantitative replicate plant-aphid-parasitoid networks constructed using DNA-based methods. Simulated climate warming affected species richness, significantly altered consumer-resource asymmetries and reduced network complexity. Increased temperature induced an aphid outbreak, but the parasitism rates of aphids by parasitoid wasps remained unchanged. It also drove changes in the crop, altering in particular the phenology of the wheat as well as its quality (i.e., fewer, lighter seeds). We discuss the importance of considering the wider impacts of climate change on interacting species across trophic levels in agroecosystems.

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.

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.

ReproPhylo: An Environment for Reproducible Phylogenomics.

The reproducibility of experiments is key to the scientific process, and particularly necessary for accurate reporting of analyses in data-rich fields such as phylogenomics. We present ReproPhylo, a phylogenomic analysis environment developed to ensure experimental reproducibility, to facilitate the handling of large-scale data, and to assist methodological experimentation. Reproducibility, and instantaneous repeatability, is built in to the ReproPhylo system and does not require user intervention or configuration because it stores the experimental workflow as a single, serialized Python object containing explicit provenance and environment information. This 'single file' approach ensures the persistence of provenance across iterations of the analysis, with changes automatically managed by the version control program Git. This file, along with a Git repository, are the primary reproducibility outputs of the program. In addition, ReproPhylo produces an extensive human-readable report and generates a comprehensive experimental archive file, both of which are suitable for submission with publications. The system facilitates thorough experimental exploration of both parameters and data. ReproPhylo is a platform independent CC0 Python module and is easily installed as a Docker image or a WinPython self-sufficient package, with a Jupyter Notebook GUI, or as a slimmer version in a Galaxy distribution.

Determining plant-leaf miner-parasitoid interactions: a DNA barcoding approach.

A major challenge in network ecology is to describe the full-range of species interactions in a community to create highly-resolved food-webs. We developed a molecular approach based on DNA full barcoding and mini-barcoding to describe difficult to observe plant-leaf miner-parasitoid interactions, consisting of animals commonly regarded as agricultural pests and their natural enemies. We tested the ability of universal primers to amplify the remaining DNA inside leaf miner mines after the emergence of the insect. We compared the results of a) morphological identification of adult specimens; b) identification based on the shape of the mines; c) the COI Mini-barcode (130 bp) and d) the COI full barcode (658 bp) fragments to accurately identify the leaf-miner species. We used the molecular approach to build and analyse a tri-partite ecological network of plant-leaf miner-parasitoid interactions. We were able to detect the DNA of leaf-mining insects within their feeding mines on a range of host plants using mini-barcoding primers: 6% for the leaves collected empty and 33% success after we observed the emergence of the leaf miner. We suggest that the low amplification success of leaf mines collected empty was mainly due to the time since the adult emerged and discuss methodological improvements. Nevertheless our approach provided new species-interaction data for the ecological network. We found that the 130 bp fragment is variable enough to identify all the species included in this study. Both COI fragments reveal that some leaf miner species could be composed of cryptic species. The network built using the molecular approach was more accurate in describing tri-partite interactions compared with traditional approaches based on morphological criteria.

The evolution of tyrosine-recombinase elements in Nematoda.

Transposable elements can be categorised into DNA and RNA elements based on their mechanism of transposition. Tyrosine recombinase elements (YREs) are relatively rare and poorly understood, despite sharing characteristics with both DNA and RNA elements. Previously, the Nematoda have been reported to have a substantially different diversity of YREs compared to other animal phyla: the Dirs1-like YRE retrotransposon was encountered in most animal phyla but not in Nematoda, and a unique Pat1-like YRE retrotransposon has only been recorded from Nematoda. We explored the diversity of YREs in Nematoda by sampling broadly across the phylum and including 34 genomes representing the three classes within Nematoda. We developed a method to isolate and classify YREs based on both feature organization and phylogenetic relationships in an open and reproducible workflow. We also ensured that our phylogenetic approach to YRE classification identified truncated and degenerate elements, informatively increasing the number of elements sampled. We identified Dirs1-like elements (thought to be absent from Nematoda) in the nematode classes Enoplia and Dorylaimia indicating that nematode model species do not adequately represent the diversity of transposable elements in the phylum. Nematode Pat1-like elements were found to be a derived form of another Pat1-like element that is present more widely in animals. Several sequence features used widely for the classification of YREs were found to be homoplasious, highlighting the need for a phylogenetically-based classification scheme. Nematode model species do not represent the diversity of transposable elements in the phylum.

The complex hybrid origins of the root knot nematodes revealed through comparative genomics.

Root knot nematodes (RKN) can infect most of the world's agricultural crop species and are among the most important of all plant pathogens. As yet however we have little understanding of their origins or the genomic basis of their extreme polyphagy. The most damaging pathogens reproduce by obligatory mitotic parthenogenesis and it has been suggested that these species originated from interspecific hybridizations between unknown parental taxa. We have sequenced the genome of the diploid meiotic parthenogen Meloidogyne floridensis, and use a comparative genomic approach to test the hypothesis that this species was involved in the hybrid origin of the tropical mitotic parthenogen Meloidogyne incognita. Phylogenomic analysis of gene families from M. floridensis, M. incognita and an outgroup species Meloidogyne hapla was carried out to trace the evolutionary history of these species' genomes, and we demonstrate that M. floridensis was one of the parental species in the hybrid origins of M. incognita. Analysis of the M. floridensis genome itself revealed many gene loci present in divergent copies, as they are in M. incognita, indicating that it too had a hybrid origin. The triploid M. incognita is shown to be a complex double-hybrid between M. floridensis and a third, unidentified, parent. The agriculturally important RKN have very complex origins involving the mixing of several parental genomes by hybridization and their extreme polyphagy and success in agricultural environments may be related to this hybridization, producing transgressive variation on which natural selection can act. It is now clear that studying RKN variation via individual marker loci may fail due to the species' convoluted origins, and multi-species population genomics is essential to understand the hybrid diversity and adaptive variation of this important species complex. This comparative genomic analysis provides a compelling example of the importance and complexity of hybridization in generating animal species diversity more generally.

Gene duplication in an African cichlid adaptive radiation.

BACKGROUND: Gene duplication is a source of evolutionary innovation and can contribute to the divergence of lineages; however, the relative importance of this process remains to be determined. The explosive divergence of the African cichlid adaptive radiations provides both a model for studying the general role of gene duplication in the divergence of lineages and also an exciting foray into the identification of genomic features that underlie the dramatic phenotypic and ecological diversification in this particular lineage. We present the first genome-wide study of gene duplication in African cichlid fishes, identifying gene duplicates in three species belonging to the Lake Malawi adaptive radiation (Metriaclima estherae, Protomelas similis, Rhamphochromis "chilingali") and one closely related species from a non-radiated riverine lineage (Astatotilapia tweddlei). RESULTS: Using Astatotilapia burtoni as reference, microarray comparative genomic hybridization analysis of 5689 genes reveals 134 duplicated genes among the four cichlid species tested. Between 51 and 55 genes were identified as duplicated in each of the three species from the Lake Malawi radiation, representing a 38%-49% increase in number of duplicated genes relative to the non-radiated lineage (37 genes). Duplicated genes include several that are involved in immune response, ATP metabolism and detoxification. CONCLUSIONS: These results contribute to our understanding of the abundance and type of gene duplicates present in cichlid fish lineages. The duplicated genes identified in this study provide candidates for the analysis of functional relevance with regard to phenotype and divergence. Comparative sequence analysis of gene duplicates can address the role of positive selection and adaptive evolution by gene duplication, while further study across the phylogenetic range of cichlid radiations (and more generally in other adaptive radiations) will determine whether the patterns of gene duplication seen in this study consistently accompany rapid radiation.

Divergence of scent pheromones in allopatric populations of Acanthodactylus boskianus (Squamata: Lacertidae).

The evolutionary phenomena associated with divergence in chemical signals between populations of the same species help to understand the process of speciation. Animals detect and react to semiochemicals and pheromones used in communication. Comparison between populations of the same species that are geographically isolated from one another allows us to determine the genetic or environmental factors responsible for chemical differentiation. Acanthodactylus boskianus from the east and west of Egypt were used as an example to compare the geographical diversity in chemical fingerprints of this species' femoral gland secretions and its phylogeography. Chemical analysis via GC-MS showed that the two geographically distinct populations' odor fingerprints are quantitatively different despite sharing the same components of the secretions. Phylogenetic analysis showed that the eastern and western Egyptian populations are genetically distinct and that chemical divergence of these lizards' odor profiles may be an example of signal evolution.

Comparative analysis of teleost genome sequences reveals an ancient intron size expansion in the zebrafish lineage.

We have developed a bioinformatics pipeline for the comparative evolutionary analysis of Ensembl genomes and have used it to analyze the introns of the five available teleost fish genomes. We show our pipeline to be a powerful tool for revealing variation between genomes that may otherwise be overlooked with simple summary statistics. We identify that the zebrafish, Danio rerio, has an unusual distribution of intron sizes, with a greater number of larger introns in general and a notable peak in the frequency of introns of approximately 500 to 2,000 bp compared with the monotonically decreasing frequency distributions of the other fish. We determine that 47% of D. rerio introns are composed of repetitive sequences, although the remainder, over 331 Mb, is not. Because repetitive elements may be the origin of the majority of all noncoding DNA, it is likely that the remaining D. rerio intronic sequence has an ancient repetitive origin and has since accumulated so many mutations that it can no longer be recognized as such. To study such an ancient expansion of repeats in the Danio, lineage will require further comparative analysis of fish genomes incorporating a broader distribution of teleost lineages.

Low endemism, continued deep-shallow interchanges, and evidence for cosmopolitan distributions in free-living marine nematodes (order Enoplida).

BACKGROUND: Nematodes represent the most abundant benthic metazoa in one of the largest habitats on earth, the deep sea. Characterizing major patterns of biodiversity within this dominant group is a critical step towards understanding evolutionary patterns across this vast ecosystem. The present study has aimed to place deep-sea nematode species into a phylogenetic framework, investigate relationships between shallow water and deep-sea taxa, and elucidate phylogeographic patterns amongst the deep-sea fauna. RESULTS: Molecular data (18 S and 28 S rRNA) confirms a high diversity amongst deep-sea Enoplids. There is no evidence for endemic deep-sea lineages in Maximum Likelihood or Bayesian phylogenies, and Enoplids do not cluster according to depth or geographic location. Tree topologies suggest frequent interchanges between deep-sea and shallow water habitats, as well as a mixture of early radiations and more recently derived lineages amongst deep-sea taxa. This study also provides convincing evidence of cosmopolitan marine species, recovering a subset of Oncholaimid nematodes with identical gene sequences (18 S, 28 S and cox1) at trans-Atlantic sample sites. CONCLUSIONS: The complex clade structures recovered within the Enoplida support a high global species richness for marine nematodes, with phylogeographic patterns suggesting the existence of closely related, globally distributed species complexes in the deep sea. True cosmopolitan species may additionally exist within this group, potentially driven by specific life history traits of Enoplids. Although this investigation aimed to intensively sample nematodes from the order Enoplida, specimens were only identified down to genus (at best) and our sampling regime focused on an infinitesimal small fraction of the deep-sea floor. Future nematode studies should incorporate an extended sample set covering a wide depth range (shelf, bathyal, and abyssal sites), utilize additional genetic loci (e.g. mtDNA) that are informative at the species level, and apply high-throughput sequencing methods to fully assay community diversity. Finally, further molecular studies are needed to determine whether phylogeographic patterns observed in Enoplids are common across other ubiquitous marine groups (e.g. Chromadorida, Monhysterida).

Moving towards a complete molecular framework of the Nematoda: a focus on the Enoplida and early-branching clades.

BACKGROUND: The subclass Enoplia (Phylum Nematoda) is purported to be the earliest branching clade amongst all nematode taxa, yet the deep phylogeny of this important lineage remains elusive. Free-living marine species within the order Enoplida play prominent roles in marine ecosystems, but previous molecular phylogenies have provided only the briefest evolutionary insights; this study aimed to firmly resolve internal relationships within the hyper-diverse but poorly understood Enoplida. In addition, we revisited the molecular framework of the Nematoda using a rigorous phylogenetic approach in order to investigate patterns of early splits amongst the oldest lineages (Dorylaimia and Enoplia). RESULTS: Morphological identifications, nuclear gene sequences (18S and 28S rRNA), and mitochondrial gene sequences (cox1) were obtained from marine Enoplid specimens representing 37 genera. The 18S gene was used to resolve deep splits within the Enoplia and evaluate the branching order of major clades in the nematode tree; multiple phylogenetic methods and rigorous empirical tests were carried out to assess tree topologies under different parameters and combinations of taxa. Significantly increased taxon sampling within the Enoplida resulted in a well-supported, robust phylogenetic topology of this group, although the placement of certain clades was not fully resolved. Our analysis could not unequivocally confirm the earliest splits in the nematode tree, and outgroup choice significantly affected the observed branching order of the Dorylaimia and Enoplia. Both 28S and cox1 were too variable to infer deep phylogeny, but provided additional insight at lower taxonomic levels. CONCLUSIONS: Analysis of internal relationships reveals that the Enoplia is split into two main clades, with groups consisting of terrestrial (Triplonchida) and primarily marine fauna (Enoplida). Five independent lineages were recovered within the Enoplida, containing a mixture of marine and terrestrial species; clade structure suggests that habitat transitions have occurred at least four times within this group. Unfortunately, we were unable to obtain a consistent or well-supported topology amongst early-branching nematode lineages. It appears unlikely that single-gene phylogenies using the conserved 18S gene will be useful for confirming the branching order at the base of the nematode tree-future efforts will require multi-gene analyses or phylogenomic methods.

Genetic tests of ancient asexuality in root knot nematodes reveal recent hybrid origins.

BACKGROUND: The existence of "ancient asexuals", taxa that have persisted for long periods of evolutionary history without sexual recombination, is both controversial and important for our understanding of the evolution and maintenance of sexual reproduction. A lack of sex has consequences not only for the ecology of the asexual organism but also for its genome. Several genetic signatures are predicted from long-term asexual (apomictic) reproduction including (i) large "allelic" sequence divergence (ii) lack of phylogenetic clustering of "alleles" within morphological species and (iii) decay and loss of genes specific to meiosis and sexual reproduction. These genetic signatures can be hard to assess since it is difficult to demonstrate the allelic nature of very divergent sequences, divergence levels may be complicated by processes such as inter-specific hybridization, and genes may have secondary roles unrelated to sexual reproduction. Apomictic species of Meloidogyne root knot nematodes have been suggested previously to be ancient asexuals. Their relatives reproduce by meiotic parthenogenesis or facultative sexuality, which in combination with the abundance of nematode genomic sequence data, makes them a powerful system in which to study the consequences of reproductive mode on genomic divergence. RESULTS: Here, sequences from nuclear protein-coding genes are used to demonstrate that the first two predictions of ancient asexuality are found within the apomictic root knot nematodes. Alleles are more divergent in the apomictic taxa than in those species exhibiting recombination and do not group phylogenetically according to recognized species. In contrast some nuclear alleles, and mtDNA, are almost identical across species. Sequencing of Major Sperm Protein, a gamete-specific gene, from both meiotic and ameiotic species reveals no increase in evolutionary rate nor change in substitution pattern in the apomictic taxa, indicating that the locus has been maintained by selection. CONCLUSION: The data strongly suggests the tropical root knot nematode apomicts have a recent origin and are not anciently asexual. The results support that interspecific hybridization has been involved in the origin of this asexual group and has played a role in shaping the patterns of genetic diversity observed. This study suggests that genetic signatures of ancient asexuality must be taken with caution due to the confounding effect of interspecific hybridization, which has long been implicated in the origins of apomictic species.