MinION Long-read Sequencing Research Articles

Chromosome-level subgenome-aware de novo assembly provides insight into Saccharomyces bayanus genome divergence after hybridization.

Interspecies hybridization is prevalent in various eukaryotic lineages and plays important roles in phenotypic diversification, adaptation, and speciation. To better understand the changes that occurred in the different subgenomes of a hybrid species and how they facilitate adaptation, we have completed chromosome-level de novo assemblies of all chromosomes for a recently formed hybrid yeast, Saccharomyces bayanus strain CBS380, using Oxford Nanopore Technologies' MinION long-read sequencing. We characterize the S. bayanus genome and compare it with its parent species, Saccharomyces uvarum and Saccharomyces eubayanus, and other S. bayanus genomes to better understand genome evolution after a relatively recent hybridization event. We observe multiple recombination events between the subgenomes in each chromosome, followed by loss of heterozygosity (LOH) in nine chromosome pairs. In addition to maintaining nearly all gene content and synteny from its parental genomes, S. bayanus has acquired many genes from other yeast species, primarily through the introgression of Saccharomyces cerevisiae, such as those involved in the maltose metabolism. Finally, the patterns of recombination and LOH suggest an allotetraploid origin of S. bayanus The gene acquisition and rapid LOH in the hybrid genome probably facilitated its adaptation to maltose brewing environments and mitigated the maladaptive effect of hybridization. This paper describes the first in-depth study using long-read sequencing technology of an S. bayanus hybrid genome, which may serve as an excellent reference for future studies of this important yeast and other yeast strains.

New Variants in the Chloroplast Genome Sequence of Two Colombian Individuals of the Cedar Timber Species (Cedrela odorata L.), Using Long-Read Oxford Nanopore Technology

The plant species Cedrela odorata has been largely exploited in the timber industry due to the high demand for its wood. Therefore, C. odorata has been considered a vulnerable species since 1994, under the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES). C. odorata is a key timber species included in the management and conservation plans for the Amazon and Central American rainforests. These plans include the development of genetic and genomic resources to study local populations of the species in Colombia. In this study, two novel chloroplast (cp) genomes were assembled and annotated using the MinION long-read sequencing technology. The new cp genomes were screened for sequence variants (SVs), and a total of 16 SNPs were identified, presumably unique to populations from the Amazon region in Colombia. Comparative genomics with other reported cp genomes from different populations of C. odorata support the hypothesis of intraspecific diversity associated with their place of origin. These cp genome sequences of C. odorata from Colombian individuals represent valuable genomic resources for the species, suitable for identifying novel DNA fingerprinting and barcoding applications.

Longitudinal Analysis of Urban Stormwater Microbiome and Resistome from Watersheds with and without Green Infrastructure using Long-Read Sequencing

Longitudinal Analysis of Urban Stormwater Microbiome and Resistome from Watersheds with and without Green Infrastructure using Long-Read Sequencing

Distribution and expression of the aac(6')-Im (aacA16) aminoglycoside resistance gene.

The aac(6')-Im (aacA16) amikacin, netilmicin and tobramycin resistance gene cassette had been circulating globally undetected for many years in a sublineage of Acinetobacter baumannii global clone 2. To identify sources for the aac(6')-Im fragment found in A. baumannii. MinION long-read sequencing and Unicycler hybrid assemblies were used to determine the genetic context of the aac(6')-Im gene. Quantitative reverse transcriptase PCR was used to measure expression. Among >60 000 non-Acinetobacter draft genomes in the MRSN collection, the aac(6')-Im gene was detected in Pseudomonas putida and Enterobacter hormaechei isolates recovered from patients in Thailand between 2016 and 2019. Genomes of multiply resistant P. putida MRSN365855 and E. hormaechei MRSN791417 were completed. The class 1 integron containing the aac(6')-Im cassette was in the chromosome in MRSN365855, and in an HI2 plasmid in MRSN791417. However, MRSN791417 was amikacin susceptible and the gene was not expressed due to loss of the Pc promoter of the integron. Further examples of aac(6')-Im in plasmids from or the chromosome of various Gram-negative species were found in the GenBank nucleotide database. The aac(6')-Im context in integrons in pMRSN791417-8 and a Klebsiella plasmid pAMR200031 shared similarities with the aac(6')-Im region of AbGRI2-Im islands in A. baumannii. In other cases, the cassette array including the aac(6')-Im cassette was different. The aac(6')-Im gene is widespread, being found so far in several different species and in several different gene cassette arrays. The lack of amikacin resistance in E. hormaechei highlights the importance of correlating resistance gene content and antibiotic resistance phenotype.

Household environment and animal fecal contamination are critical modifiers of the gut microbiome and resistome in young children from rural Nicaragua

BackgroundEarly life plays a vital role in the development of the gut microbiome and subsequent health. While many factors that shape the gut microbiome have been described, including delivery mode, breastfeeding, and antibiotic use, the role of household environments is still unclear. Furthermore, the development of the gut antimicrobial resistome and its role in health and disease is not well characterized, particularly in settings with water insecurity and less sanitation infrastructure.ResultsThis study investigated the gut microbiome and resistome of infants and young children (ages 4 days-6 years) in rural Nicaragua using Oxford Nanopore Technology’s MinION long-read sequencing. Differences in gut microbiome diversity and antibiotic resistance gene (ARG) abundance were examined for associations with host factors (age, sex, height for age z-score, weight for height z-score, delivery mode, breastfeeding habits) and household environmental factors (animals inside the home, coliforms in drinking water, enteric pathogens in household floors, fecal microbial source tracking markers in household floors). We identified anticipated associations of higher gut microbiome diversity with participant age and vaginal delivery. However, novel to this study were the significant, positive associations between ruminant and dog fecal contamination of household floors and gut microbiome diversity. We also identified greater abundance of potential pathogens in the gut microbiomes of participants with higher fecal contamination on their household floors. Path analysis revealed that water quality and household floor contamination independently and significantly influenced gut microbiome diversity when controlling for age. These gut microbiome contained diverse resistome, dominated by multidrug, tetracycline, macrolide/lincosamide/streptogramin, and beta-lactam resistance. We found that the abundance of ARGs in the gut decreased with age. The bacterial hosts of ARGs were mainly from the family Enterobacteriaceae, particularly Escherichia coli.ConclusionsThis study identified the role of household environmental contamination in the developing gut microbiome and resistome of young children and infants with a One Health perspective. We found significant relationships between host age, gut microbiome diversity, and the resistome. Understanding the impact of the household environment on the development of the resistome and microbiome in early life is essential to optimize the relationship between environmental exposure and human health.92Jf2w8Y_xc1zDLBRpWCqWVideo Graphical

Division of Blastocystis ST10 into three new subtypes: ST42-ST44.

The Blastocystis subtype ST10 has been recognized to contain a great deal of diversity at the sequence level, potentially indicating the presence of multiple new STs within the clade. However, the data needed to validate these new STs were not available. To help resolve this diversity, full-length small subunit (SSU) rRNA gene reference sequences were generated using Oxford Nanopore MinION long-read sequencing from 21 samples representing multiple domestic and wild hosts and geographic regions and covering the sequence diversity previously described using fragments of the SSU rRNA gene. Phylogenetic and pairwise distance analyses were used to compare full-length sequences of the SSU rRNA gene generated in this study with all other valid STs of Blastocystis. We present data supporting the division of ST10/ST23 cluster into five subtypes, ST10, ST23, and three new subtypes with the proposed ST designations of ST42, ST43, and ST44. As the host range of Blastocystis continues to expand with new subtypes and new hosts being frequently identified, the reference sequences provided in this study will assist in accurate sequence classification and help to clarify the epidemiology of this common intestinal microeukaryote.

Identification and validation of novel Blastocystis subtype ST41 in a Colombian patient undergoing colorectal cancer screening.

Blastocystis sp. is among the most frequent intestinal protists identified in humans globally. However, characterization of Blastocystis subtype diversity in humans is ongoing. We report here the identification of novel Blastocystis subtype ST41 in a Colombian patient undergoing colorectal cancer screening involving colonoscopy and fecal testing (microscopy, culture, PCR). The full-length ssu rRNA gene sequence of the protist was generated using MinION long-read sequencing technology. The validity of the novel subtype was confirmed via phylogenetic and pairwise distance analyses of the full-length ST41 sequence and all other valid subtypes. The study provides reference material essential for conducting subsequent experimental studies.

A preliminary study of the use of MinION sequencing to specifically detect Shiga toxin-producing Escherichia coli in culture swipes containing multiple serovars of this species

An important challenge relating to clinical diagnostics of the foodborne pathogen Shiga toxin-producing E. coli (STEC), is that PCR-detection of the shiga-toxin gene (stx) in DNA from stool samples can be accompanied by a failure to identify an STEC isolate in pure culture on agar. In this study, we have explored the use of MinION long-read sequencing of DNA from bacterial culture swipes to detect the presence of STEC, and bioinformatic tools to characterize the STEC virulence factors. The online workflow “What’s in my pot” (WIMP) in the Epi2me cloud service, rapidly identified STEC also when it was present in culture swipes together with multiple other E. coli serovars, given sufficient abundance. These preliminary results provide useful information about the sensitivity of the method, which has potential to be used in clinical diagnostic of STEC, particularly in cases where a pure culture of the STEC isolate is not obtained due to the ‘STEC lost Shiga toxin’ phenomenon.

Structural genome variants of Pseudomonas aeruginosa clone C and PA14 strains.

Plasticity of Pseudomonas aeruginosa chromosomes is mainly driven by an extended accessory genome that is shaped by insertion and deletion events. Further modification of the genome composition can be induced by chromosomal inversion events which lead to relocation of genes in the affected genomic DNA segments, modify the otherwise highly conserved core genome synteny and could even alter the location of the replication terminus. Although the genome of the first sequenced strain, PAO1, displayed such a large genomic inversion, knowledge on such recombination events in the P. aeruginosa population is limited. Several large inversions had been discovered in the late 1990s in cystic fibrosis isolates of the major clonal lineage C by physical genome mapping, and subsequent work on these examples led to the characterization of the DNA at the recombination breakpoints and a presumed recombination mechanism. Since then, the topic was barely addressed in spite of the compilation of thousands of P. aeruginosa genome sequences that are deposited in databases. Due to the use of second-generation sequencing, genome contig assembly had usually followed synteny blueprints provided by the existing reference genome sequences. Inversion detection was not feasible by these approaches, as the respective read lengths did not allow reliable resolution of sequence repeats that are typically found at the borders of inverted segments. In this study, we applied PacBio and MinION long-read sequencing to isolates of the mentioned clone C collection. Confirmation of inversions predicted from the physical mapping data demonstrated that unbiased sequence assembly of such read datasets allows the detection of genomic inversions and the resolution of the recombination breakpoint regions. Additional long-read sequencing of representatives of the other major clonal lineage, PA14, revealed large inversions in several isolates, from cystic fibrosis origin as well as from other sources. These findings indicated that inversion events are not restricted to strains from chronic infection background, but could be widespread in the P. aeruginosa population and contribute to genome plasticity. Moreover, the monitored examples emphasized the role of small mobile DNA units, such as IS elements or transposons, and accessory DNA elements in the inversion-related recombination processes.

Functional genome annotation and transcriptome analysis of Pseudozyma hubeiensis BOT-O, an oleaginous yeast that utilizes glucose and xylose at equal rates

Pseudozyma hubeiensis is a basidiomycete yeast that has the highly desirable traits for lignocellulose valorisation of being equally efficient at utilization of glucose and xylose, and capable of their co-utilization. The species has previously mainly been studied for its capacity to produce secreted biosurfactants in the form of mannosylerythritol lipids, but it is also an oleaginous species capable of accumulating high levels of triacylglycerol storage lipids during nutrient starvation. In this study, we aimed to further characterize the oleaginous nature of P. hubeiensis by evaluating metabolism and gene expression responses during storage lipid formation conditions with glucose or xylose as a carbon source.The genome of the recently isolated P. hubeiensis BOT-O strain was sequenced using MinION long-read sequencing and resulted in the most contiguous P. hubeiensis assembly to date with 18.95 Mb in 31 contigs. Using transcriptome data as experimental support, we generated the first mRNA-supported P. hubeiensis genome annotation and identified 6540 genes. 80% of the predicted genes were assigned functional annotations based on protein homology to other yeasts. Based on the annotation, key metabolic pathways in BOT-O were reconstructed, including pathways for storage lipids, mannosylerythritol lipids and xylose assimilation. BOT-O was confirmed to consume glucose and xylose at equal rates, but during mixed glucose-xylose cultivation glucose was found to be taken up faster. Differential expression analysis revealed that only a total of 122 genes were significantly differentially expressed at a cut-off of |log2 fold change| ≥ 2 when comparing cultivation on xylose with glucose, during exponential growth and during nitrogen-starvation. Of these 122 genes, a core-set of 24 genes was identified that were differentially expressed at all time points. Nitrogen-starvation resulted in a larger transcriptional effect, with a total of 1179 genes with significant expression changes at the designated fold change cut-off compared with exponential growth on either glucose or xylose.

Evaluation of high molecular weight DNA extraction methods for long-read sequencing of Shiga toxin-producing Escherichia coli.

Next generation sequencing has become essential for pathogen characterization and typing. The most popular second generation sequencing technique produces data of high quality with very low error rates and high depths. One major drawback of this technique is the short reads. Indeed, short-read sequencing data of Shiga toxin-producing Escherichia coli (STEC) are difficult to assemble because of the presence of numerous mobile genetic elements (MGEs), which contain repeated elements. The resulting draft assemblies are often highly fragmented, which results in a loss of information, especially concerning MGEs or large structural variations. The use of long-read sequencing can circumvent these problems and produce complete or nearly complete genomes. The ONT MinION, for its small size and minimal investment requirements, is particularly popular. The ultra-long reads generated with the MinION can easily span prophages and repeat regions. In order to take full advantage of this technology it requires High Molecular Weight (HMW) DNA of high quality in high quantity. In this study, we have tested three different extraction methods: bead-based, solid-phase and salting-out, and evaluated their impact on STEC DNA yield, quality and integrity as well as performance in MinION long-read sequencing. Both the bead-based and salting-out methods allowed the recovery of large quantities of HMW STEC DNA suitable for MinION library preparation. The DNA extracted using the salting-out method consistently produced longer reads in the subsequent MinION runs, compared with the bead-based methods. While both methods performed similarly in subsequent STEC genome assembly, DNA extraction based on salting-out appeared to be the overall best method to produce high quantity of pure HMW STEC DNA for MinION sequencing.

Staphylococcal cassette chromosome mec amplification as a mechanism for ceftobiprole resistance in clinical methicillin-resistant Staphylococcus aureus isolates

Staphylococcal cassette chromosome mec amplification as a mechanism for ceftobiprole resistance in clinical methicillin-resistant Staphylococcus aureus isolates

Genome Resources for the Ex-Type of Phytophthora citricola, and Well-Authenticated Isolates of P. hibernalis, P. nicotianae, and P. syringae.

Phytophthora is one of the most important genera of plant pathogens, with many members causing high economic losses worldwide. To build robust molecular identification systems, it is very important to have information from well-authenticated specimens and, in preference, the ex-type specimens. The reference genomes of well-authenticated specimens form a critical foundation for genetics, biological research, and diagnostic applications. In this study, we describe four draft Phytophthora genome resources for the ex-type of Phytophthora citricola BL34 (P0716 WPC) (118 contigs for 50 Mb), and well-authenticated specimens of P. syringae BL57G (P10330 WPC) (591 contigs for 75 Mb), P. hibernalis BL41G (P3822 WPC) (404 contigs for 84 Mb), and P. nicotianae BL162 (P6303 WPC) (3,984 contigs for 108 Mb) generated with MinION long-read high-throughput sequencing technology (Oxford Nanopore Technologies). Using the quality reads, we assembled high-coverage genomes of P. citricola with 291× coverage and 16,662 annotated genes; P. nicotianae with 205× coverage and 29,271 annotated genes; P. syringae with 76× coverage and 23,331 annotated genes, and P. hibernalis with 42× coverage and 21,762 annotated genes. With the availability of genome sequences and their annotations, we predict that these draft genomes will be accommodating for various basic and applied research, including diagnostics to protect global agriculture.

Integrative profiling of Epstein\u2013Barr virus transcriptome using a multiplatform approach

BackgroundEpstein–Barr virus (EBV) is an important human pathogenic gammaherpesvirus with carcinogenic potential. The EBV transcriptome has previously been analyzed using both Illumina-based short read-sequencing and Pacific Biosciences RS II-based long-read sequencing technologies. Since the various sequencing methods have distinct strengths and limitations, the use of multiplatform approaches have proven to be valuable. The aim of this study is to provide a more complete picture on the transcriptomic architecture of EBV.MethodsIn this work, we apply the Oxford Nanopore Technologies MinION (long-read sequencing) platform for the generation of novel transcriptomic data, and integrate these with other’s data generated by another LRS approach, Pacific BioSciences RSII sequencing and Illumina CAGE-Seq and Poly(A)-Seq approaches. Both amplified and non-amplified cDNA sequencings were applied for the generation of sequencing reads, including both oligo-d(T) and random oligonucleotide-primed reverse transcription. EBV transcripts are identified and annotated using the LoRTIA software suite developed in our laboratory.ResultsThis study detected novel genes embedded into longer host genes containing 5′-truncated in-frame open reading frames, which potentially encode N-terminally truncated proteins. We also detected a number of novel non-coding RNAs and transcript length isoforms encoded by the same genes but differing in their start and/or end sites. This study also reports the discovery of novel splice isoforms, many of which may represent altered coding potential, and of novel replication-origin-associated transcripts. Additionally, novel mono- and multigenic transcripts were identified. An intricate meshwork of transcriptional overlaps was revealed.ConclusionsAn integrative approach applying multi-technique sequencing technologies is suitable for reliable identification of complex transcriptomes because each techniques has different advantages and limitations, and the they can be used for the validation of the results obtained by a particular approach.

Accelerating Complete Phytoplasma Genome Assembly by Immunoprecipitation-Based Enrichment and MinION-Based DNA Sequencing for Comparative Analyses.

Phytoplasmas are uncultivated plant-pathogenic bacteria with agricultural importance. Those belonging to the 16SrII group, represented by ‘Candidatus P. aurantifolia’, have a wide range of plant hosts and cause significant yield losses in valuable crops, such as pear, sweet potato, peanut, and soybean. In this study, a method that combines immunoprecipitation-based enrichment and MinION long-read DNA sequencing was developed to solve the challenge of phytoplasma genome studies. This approach produced long reads with high mapping rates and high genomic coverage that can be combined with Illumina reads to produce complete genome assemblies with high accuracy. We applied this method to strain NCHU2014 and determined its complete genome sequence, which consists of one circular chromosome with 635,584 bp and one plasmid with 4,224 bp. Although ‘Ca. P. aurantifolia’ NCHU2014 has a small chromosome with only 471 protein-coding genes, it contains 33 transporter genes and 27 putative effector genes, which may contribute to obtaining nutrients from hosts and manipulating host developments for their survival and multiplication. Two effectors, the homologs of SAP11 and SAP54/PHYL1 identified in ‘Ca. P. aurantifolia’ NCHU2014, have the biochemical activities in destabilizing host transcription factors, which can explain the disease symptoms observed in infected plants. Taken together, this study provides the first complete genome available for the 16SrII phytoplasmas and contributes to the understanding of phytoplasma pathogenicity.

Cunning plasmid fusion mediates antibiotic resistance genes represented by ESBLs encoding genes transfer in foodborne Salmonella

Cunning plasmid fusion mediates antibiotic resistance genes represented by ESBLs encoding genes transfer in foodborne Salmonella

Molecular Characterization of Multidrug-Resistant Yersinia enterocolitica From Foodborne Outbreaks in Sweden.

The foodborne pathogen Yersinia enterocolitica causes gastrointestinal infections worldwide. In the spring of 2019, the Swedish Public Health Agency and Statens Serum Institut in Denmark independently identified an outbreak caused by Yersinia enterocolitica 4/O:3 that after sequence comparison turned out to be a cross-border outbreak. A trace-back investigation suggested shipments of fresh prewashed spinach from Italy as a common source for the outbreak. Here, we determined the genome sequences of five Y. enterocolitica clinical isolates during the Swedish outbreak using a combination of Illumina HiSeq short-read and Nanopore Technologies’ MinION long-read whole-genome sequencing. WGS results showed that all clinical strains have a fully assembled chromosome of approximately 4.6 Mbp in size and a 72-kbp virulence plasmid; one of the strains was carrying an additional 5.7-kbp plasmid, pYE-tet. All strains showed a high pathogen probability score (87.5%) with associated genes for virulence, all of which are closely related to an earlier clinical strain Y11 from Germany. In addition, we identified a chromosomally encoded multidrug-resistance cassette carrying resistance genes against chloramphenicol (catA1), streptomycin (aadA1), sulfonamides (sul1), and a mercury resistance module. This chromosomally encoded Tn2670 transposon has previously been reported associated with IncFII plasmids in Enterobacteriaceae: a Shigella flexneri clinical isolate from Japan in 1950s, a Klebsiella pneumoniae outbreak from Australia in 1997, and Salmonella enterica serovar Typhimurium. Interestingly, we identified an additional 5.7-kbp plasmid with tetB (encoding an ABC transporter), Rep, and its own ORI and ORIt sites, sharing high homology with small tetB-Rep plasmids from Pasteurellaceae. This is the first time that Tn2670 and Pasteurellaceae plasmids have been reported in Y. enterocolitica. Taken together, our study showed that the Swedish Y. enterocolitica outbreak strains acquired multi-antibiotic and metal-resistance genes through horizontal gene transfer, suggesting a potential reservoir of intraspecies dissemination of multidrug-resistance genes among foodborne pathogens. This study also highlights the concern of food-chain contamination of prewashed vegetables as a perpetual hazard against public health.

Mind the Gap: New Full-Length Sequences of Blastocystis Subtypes Generated via Oxford Nanopore Minion Sequencing Allow for Comparisons between Full-Length and Partial Sequences of the Small Subunit of the Ribosomal RNA Gene.

Blastocystis is a common food- and water-borne intestinal protist parasite of humans and many other animals. Blastocystis comprises multiple subtypes (STs) based on variability within the small subunit ribosomal (SSU rRNA) RNA gene. Though full-length reference sequences of the SSU rRNA gene are a current requirement to name a novel Blastocystis subtype, full-length reference sequences are not currently available for all subtypes. In the present study, Oxford Nanopore MinION long-read sequencing was employed to generate full-length SSU rRNA sequences for seven new Blastocystis subtypes for which no full-length references currently exist: ST21, ST23, ST24, ST25, ST26, ST27, and ST28. Phylogenetic analyses and pairwise distance matrixes were used to compare full-length and partial sequences of the two regions that are most commonly used for subtyping. Analyses included Blastocystis nucleotide sequences obtained in this study (ST21 and ST23–ST28) and existing subtypes for which full-length reference sequences were available (ST1–ST17 and ST29). The relationships and sequence variance between new and existing subtypes observed in analyses of different portions of the SSU rRNA gene are discussed. The full-length SSU rRNA reference sequences generated in this study provide essential new data to study and understand the relationships between the genetic complexity of Blastocystis and its host specificity, pathogenicity, and epidemiology.

CaBagE: A Cas9-based Background Elimination strategy for targeted, long-read DNA sequencing.

A substantial fraction of the human genome is difficult to interrogate with short-read DNA sequencing technologies due to paralogy, complex haplotype structures, or tandem repeats. Long-read sequencing technologies, such as Oxford Nanopore's MinION, enable direct measurement of complex loci without introducing many of the biases inherent to short-read methods, though they suffer from relatively lower throughput. This limitation has motivated recent efforts to develop amplification-free strategies to target and enrich loci of interest for subsequent sequencing with long reads. Here, we present CaBagE, a method for target enrichment that is efficient and useful for sequencing large, structurally complex targets. The CaBagE method leverages the stable binding of Cas9 to its DNA target to protect desired fragments from digestion with exonuclease. Enriched DNA fragments are then sequenced with Oxford Nanopore's MinION long-read sequencing technology. Enrichment with CaBagE resulted in a median of 116X coverage (range 39-416) of target loci when tested on five genomic targets ranging from 4-20kb in length using healthy donor DNA. Four cancer gene targets were enriched in a single reaction and multiplexed on a single MinION flow cell. We further demonstrate the utility of CaBagE in two ALS patients with C9orf72 short tandem repeat expansions to produce genotype estimates commensurate with genotypes derived from repeat-primed PCR for each individual. With CaBagE there is a physical enrichment of on-target DNA in a given sample prior to sequencing. This feature allows adaptability across sequencing platforms and potential use as an enrichment strategy for applications beyond sequencing. CaBagE is a rapid enrichment method that can illuminate regions of the 'hidden genome' underlying human disease.

Rapid protein evolution, organellar reductions, and invasive intronic elements in the marine aerobic parasite dinoflagellate Amoebophrya spp

BackgroundDinoflagellates are aquatic protists particularly widespread in the oceans worldwide. Some are responsible for toxic blooms while others live in symbiotic relationships, either as mutualistic symbionts in corals or as parasites infecting other protists and animals. Dinoflagellates harbor atypically large genomes (~ 3 to 250 Gb), with gene organization and gene expression patterns very different from closely related apicomplexan parasites. Here we sequenced and analyzed the genomes of two early-diverging and co-occurring parasitic dinoflagellate Amoebophrya strains, to shed light on the emergence of such atypical genomic features, dinoflagellate evolution, and host specialization.ResultsWe sequenced, assembled, and annotated high-quality genomes for two Amoebophrya strains (A25 and A120), using a combination of Illumina paired-end short-read and Oxford Nanopore Technology (ONT) MinION long-read sequencing approaches. We found a small number of transposable elements, along with short introns and intergenic regions, and a limited number of gene families, together contribute to the compactness of the Amoebophrya genomes, a feature potentially linked with parasitism. While the majority of Amoebophrya proteins (63.7% of A25 and 59.3% of A120) had no functional assignment, we found many orthologs shared with Dinophyceae. Our analyses revealed a strong tendency for genes encoded by unidirectional clusters and high levels of synteny conservation between the two genomes despite low interspecific protein sequence similarity, suggesting rapid protein evolution. Most strikingly, we identified a large portion of non-canonical introns, including repeated introns, displaying a broad variability of associated splicing motifs never observed among eukaryotes. Those introner elements appear to have the capacity to spread over their respective genomes in a manner similar to transposable elements. Finally, we confirmed the reduction of organelles observed in Amoebophrya spp., i.e., loss of the plastid, potential loss of a mitochondrial genome and functions.ConclusionThese results expand the range of atypical genome features found in basal dinoflagellates and raise questions regarding speciation and the evolutionary mechanisms at play while parastitism was selected for in this particular unicellular lineage.