ploidyNGS: visually exploring ploidy with Next Generation Sequencing data

Genetic characteristics of the diploid offsprings in potato Cooperation 88 induced by diploid donor IVP101

Diploid lines (2n = 2x = 24) derived from tetraploid potato cultivars have been utilized to hybridize with wild diploid potato species, yielding fertile offsprings. Utilizing the pollen of Solanum tuberosum Group Phureja, such as IVP101, IVP35 and IVP48, as an inducer for wide hybridization with tetraploid cultivars represents a common method for producing diploids. In this study, we created a distant hybridization induced population of tetraploid potato cultivar Cooperation 88 (C88) and IVP101, and screened all diploids using flow cytometry and ploidyNGS. We investigated the genetic composition of chloroplast and nuclear genomes in 43 diploid offsprings. We found that all diploid offsprings share the same chloroplast genomic sequence as C88 and no evidence of paternal chloroplast inheritance was found. Used SNP data to calculate the theoretical introgression index of IVP101 with diploid offsprings. The results showed that the inducer's nuclear genome was involved in the nuclear genome of the diploid offsprings with purple stem trait, indicating that the inducer nuclear genome was not completely eliminated in the nuclear genome during distant hybridization. Furthermore, we conducted a comparative analysis of the chloroplast genomes of the Solanum genus. The results indicated that (1) the chloroplast genome sizes of the 14 Solanum species ranged from 154,289 bp to 155,614 bp, with a total number of genes ranging 128-141, and with ycf1 and rps19 pseudogenes appearing at the IRB/SSC and IRA/LSC boundaries, respectively; (2) eight divergent hotspots distributed in the LSC and SSC regions of the Solanum chloroplast genomes were identified; (3) positive selection was detected in the clpP, rbcL, rps15, and rps4 genes, likely contributing to the adaptation of Solanum species to different habitats. These results reveal the variation and evolutionary characteristics of chloroplast genomes in Solanum plants.

Comparative transcriptome analysis of Cyperus esculentus and C. rotundus with contrasting oil contents in tubers defines genes and regulatory networks involved in oil accumulation

Plant vegetative organs present great potential for lipid storage, with tubers of Cyperus esculentus as a unique example. To investigate the genome and transcriptomic features of C. esculentus and related species, we sequenced and assembled the C. esculentus genome at the contig level. Through a comparative study of high-quality transcriptomes across 36 tissues from high-oil and intermediate-oil C. esculentus and low-oil Cyperus rotundus, we identified potential genes and regulatory networks related to tuber oil accumulation. First, we identified tuber-specific genes in two C. esculentus cultivars. Second, genes involved in fatty acid (FA) biosynthesis, triacylglycerol synthesis, and TAG packaging presented increased activity in the later stages of tuber development. Notably, tubers with high oil contents presented higher levels of these genes than those with intermediate oil contents did, whereas tubers with low oil contents presented minimal gene expression. Notably, a large fragment of the FA biosynthesis rate-limiting enzyme-encoding gene BCCP1 was missing from the C. rotundus transcript, which might be responsible for blocking FA biosynthesis in its tubers. WGCNA pinpointed a gene module linked to tuber oil accumulation, with a coexpression network involving the transcription factors WRI1, MYB4, and bHLH68. The ethylene-related genes in this module suggest a role for ethylene signaling in oil accumulation, which is supported by the finding that ethylene (ETH) treatment increases the oil content in C. esculentus tubers. This study identified potential genes and networks associated with tuber oil accumulation in C. esculentus, highlighting the role of specific genes, transcription factors, and ethylene signaling in this process.

Towards Digital Quantification of Ploidy from Pan-Cancer Digital Pathology Slides using Deep Learning

Abnormal DNA ploidy, found in numerous cancers, is increasingly being recognized as a contributor in driving chromosomal instability, genome evolution, and the heterogeneity that fuels cancer cell progression. Furthermore, it has been linked with poor prognosis of cancer patients. While next-generation sequencing can be used to approximate tumor ploidy, it has a high error rate for near-euploid states, a high cost and is time consuming, motivating alternative rapid quantification methods. We introduce PloiViT, a transformer-based model for tumor ploidy quantification that outperforms traditional machine learning models, enabling rapid and cost-effective quantification directly from pathology slides. We trained PloiViT on a dataset of fifteen cancer types from The Cancer Genome Atlas and validated its performance in multiple independent cohorts. Additionally, we explored the impact of self-supervised feature extraction on performance. PloiViT, using self-supervised features, achieved the lowest prediction error in multiple independent cohorts, exhibiting better generalization capabilities. Our findings demonstrate that PloiViT predicts higher ploidy values in aggressive cancer groups and patients with specific mutations, validating PloiViT potential as complementary for ploidy assessment to next-generation sequencing data. To further promote its use, we release our models as a user-friendly inference application and a Python package for easy adoption and use.

Chromosome-level genome assemblies reveal genome evolution of an invasive plant Phragmites australis

Biological invasions pose a significant threat to ecosystems, disrupting local biodiversity and ecosystem functions. The genomic underpinnings of invasiveness, however, are still largely unknown, making it difficult to predict and manage invasive species effectively. The common reed (Phragmites australis) is a dominant grass species in wetland ecosystems and has become particularly invasive when transferred from Europe to North America. Here, we present a high-quality gap-free, telomere-to-telomere genome assembly of Phragmites australis consisting of 24 pseudochromosomes and a B chromosome. Fully phased subgenomes demonstrated considerable subgenome dominance and revealed the divergence of diploid progenitors approximately 30.9 million years ago. Comparative genomics using chromosome-level scaffolds for three other lineages and a previously published draft genome assembly of an invasive lineage revealed that gene family expansions in the form of tandem duplications may have contributed to the invasiveness of the lineage. This study sheds light on the genome evolution of Arundinoideae grasses and suggests that genetic drivers, such as gene family expansions and tandem duplications, may underly the processes of biological invasion in plants. These findings provide a crucial step toward understanding and managing the genetic basis of invasiveness in plant species.

Variant calling in polyploids for population and quantitative genetics

Advancements in genome assembly and sequencing technology have made whole genome sequence (WGS) data and reference genomes accessible to study polyploid species. Compared to popular reduced-representation sequencing approaches, the genome-wide coverage and greater marker density provided by WGS data can greatly improve our understanding of polyploid species and polyploid biology. However, biological features that make polyploid species interesting also pose challenges in read mapping, variant identification, and genotype estimation. Accounting for characteristics in variant calling like allelic dosage uncertainty, homology between subgenomes, and variance in chromosome inheritance mode can reduce errors. Here, I discuss the challenges of variant calling in polyploid WGS data and discuss where potential solutions can be integrated into a standard variant calling pipeline.

Haematococcus lacustris genome assembly and annotation reveal diploid genetic traits and stress-induced gene expression patterns

The green alga Haematococcus lacustris (formerly Haematococcus pluvialis) is a primary source of astaxanthin, a ketocarotenoid with high antioxidant activity and several industrial applications. Here, the Haematococcus lacustris highly repetitive genome was reconstructed by exploiting next-generation sequencing integrated with Hi-C scaffolding, obtaining a 151 Mb genome assembly in 32 scaffolds at a near-chromosome level with high continuity. Surprisingly, the distribution of the single-nucleotide-polymorphisms identified demonstrates a diploid configuration for the Haematococcus genome, further validated by Sanger sequencing of heterozygous regions. Functional annotation and RNA-seq data enabled the identification of 13,946 nuclear genes, with >5000 genes not previously identified in this species, providing insights into the molecular basis for metabolic rear-rangement in stressing conditions such as high light and/or nitrogen starvation, where astaxanthin biosynthesis is triggered. These data constitute a rich genetic resource for biotechnological manipulation of Haematococcus lacustris highlighting potential targets to improve astaxanthin and carotenoid productivity.

Genomic characterisation and ecological distribution of Mantoniella tinhauana: a novel Mamiellophycean green alga from the Western Pacific

Mamiellophyceae are dominant marine algae in much of the ocean, the most prevalent genera belonging to the order Mamiellales: Micromonas, Ostreococcus and Bathycoccus, whose genetics and global distributions have been extensively studied. Conversely, the genus Mantoniella, despite its potential ecological importance, remains relatively under-characterised. In this study, we isolated and characterised a novel species of Mamiellophyceae, Mantoniella tinhauana, from subtropical coastal waters in the South China Sea. Morphologically, it resembles other Mantoniella species; however, a comparative analysis of the 18S and ITS2 marker genes revealed its genetic distinctiveness. Furthermore, we sequenced and assembled the first genome of Mantoniella tinhauana, uncovering significant differences from previously studied Mamiellophyceae species. Notably, the genome lacked any detectable outlier chromosomes and exhibited numerous unique orthogroups. We explored gene groups associated with meiosis, scale and flagella formation, shedding light on species divergence, yet further investigation is warranted. To elucidate the biogeography of Mantoniella tinhauana, we conducted a comprehensive analysis using global metagenomic read mapping to the newly sequenced genome. Our findings indicate this species exhibits a cosmopolitan distribution with a low-level prevalence worldwide. Understanding the intricate dynamics between Mamiellophyceae and the environment is crucial for comprehending their impact on the ocean ecosystem and accurately predicting their response to forthcoming environmental changes.

Target capture and genome skimming for plant diversity studies

Recent technological advances in long-read high-throughput sequencing and assembly methods have facilitated the generation of annotated chromosome-scale whole-genome sequence data for evolutionary studies; however, generating such data can still be difficult for many plant species. For example, obtaining high-molecular-weight DNA is typically impossible for samples in historical herbarium collections, which often have degraded DNA. The need to fast-freeze newly collected living samples to conserve high-quality DNA can be complicated when plants are only found in remote areas. Therefore, short-read reduced-genome representations, such as target capture and genome skimming, remain important for evolutionary studies. Here, we review the pros and cons of each technique for non-model plant taxa. We provide guidance related to logistics, budget, the genomic resources previously available for the target clade, and the nature of the study. Furthermore, we assess the available bioinformatic analyses, detailing best practices and pitfalls, and suggest pathways to combine newly generated data with legacy data. Finally, we explore the possible downstream analyses allowed by the type of data generated using each technique. We provide a practical guide to help researchers make the best-informed choice regarding reduced genome representation for evolutionary studies of non-model plants in cases where whole-genome sequencing remains impractical.

Genomic evidence for homoploid hybrid speciation in a marine mammal apex predator

Hybridization is widespread and constitutes an important source of genetic variability and evolution. In animals, its role in generating novel and independent lineages (hybrid speciation) has been strongly debated, with only a few cases supported by genomic data. The South American fur seal (SAfs) Arctocephalus australis is a marine apex predator of Pacific and Atlantic waters, with a disjunct set of populations in Peru and Northern Chile [Peruvian fur seal (Pfs)] with controversial taxonomic status. We demonstrate, using complete genome and reduced representation sequencing, that the Pfs is a genetically distinct species with an admixed genome that originated from hybridization between the SAfs and the Galapagos fur seal (Arctocephalus galapagoensis) ~400,000 years ago. Our results strongly support the origin of Pfs by homoploid hybrid speciation over alternative introgression scenarios. This study highlights the role of hybridization in promoting species-level biodiversity in large vertebrates.

Population genetic analysis of the microsporidium Ordospora colligata reveals the role of natural selection and phylogeography on its extremely compact and reduced genome

The determinants of variation in a species' genome-wide nucleotide diversity include historical, environmental, and stochastic aspects. This diversity can inform us about the species' past and present evolutionary dynamics. In parasites, the mode of transmission and the interactions with the host might supersede the effects of these aspects in shaping parasite genomic diversity. We used genomic samples from 10 populations of the microsporidian parasite Ordospora colligata to investigate present genomic diversity and how it was shaped by evolutionary processes, specifically, the role of phylogeography, co-phylogeography (with the host), natural selection, and transmission mode. Although very closely related microsporidia cause diseases in humans, O. colligata is specific to the freshwater crustacean Daphnia magna and has one of the smallest known eukaryotic genomes. We found an overlapping phylogeography between O. colligata and its host highlighting the long-term, intimate relationship between them. The observed geographic distribution reflects previous findings that O. colligata exhibits adaptations to colder habitats, which differentiates it from other microsporidian gut parasites of D. magna predominantly found in warmer areas. The co-phylogeography allowed us to calibrate the O. colligata phylogeny and thus estimate its mutation rate. We identified several genetic regions under potential selection. Our whole-genome study provides insights into the evolution of one of the most reduced eukaryotic genomes and shows how different processes shape genomic diversity of an obligate parasite.

ploidyfrost: Reference-free estimation of ploidy level from whole genome sequencing data based on de Bruijn graphs

Polyploidy is ubiquitous and its consequences are complex and variable. A change of ploidy level generally influences genetic diversity and results in morphological, physiological and ecological differences between cells or organisms with different ploidy levels. To avoid cumbersome experiments and take advantage of the less biased information provided by the vast amounts of genome sequencing data, computational tools for ploidy estimation are urgently needed. Until now, although a few such tools have been developed, many aspects of this estimation, such as the requirement of a reference genome, the lack of informative results and objective inferences, and the influence of false positives from errors and repeats, need further improvement. We have developed ploidyfrost, a de Bruijn graph-based method, to estimate ploidy levels from whole genome sequencing data sets without a reference genome. ploidyfrost provides a visual representation of allele frequency distribution generated using the ggplot2 package as well as quantitative results using the Gaussian mixture model. In addition, it takes advantage of colouring information encoded in coloured de Bruijn graphs to analyse multiple samples simultaneously and to flexibly filter putative false positives. We evaluated the performance of ploidyfrost by analysing highly heterozygous or repetitive samples of Cyclocarya paliurus and a complex allooctoploid sample of Fragaria × ananassa. Moreover, we demonstrated that the accuracy of analysis results can be improved by constraining a threshold such as Cramér's V coefficient on variant features, which may significantly reduce the side effects of sequencing errors and annoying repeats on the graphical structure constructed.

Systematic multi-omics cell line profiling uncovers principles of Ewing sarcoma fusion oncogene-mediated gene regulation

Ewing sarcoma (EwS) is characterized by EWSR1-ETS fusion transcription factors converting polymorphic GGAA microsatellites (mSats) into potent neo-enhancers. Although the paucity of additional mutations makes EwS a genuine model to study principles of cooperation between dominant fusion oncogenes and neo-enhancers, this is impeded by the limited number of well-characterized models. Here we present the Ewing Sarcoma Cell Line Atlas (ESCLA), comprising whole-genome, DNA methylation, transcriptome, proteome, and chromatin immunoprecipitation sequencing (ChIP-seq) data of 18 cell lines with inducible EWSR1-ETS knockdown. The ESCLA shows hundreds of EWSR1-ETS-targets, the nature of EWSR1-ETS-preferred GGAA mSats, and putative indirect modes of EWSR1-ETS-mediated gene regulation, converging in the duality of a specific but plastic EwS signature. We identify heterogeneously regulated EWSR1-ETS-targets as potential prognostic EwS biomarkers. Our freely available ESCLA (http://r2platform.com/escla/) is a rich resource for EwS research and highlights the power of comprehensive datasets to unravel principles of heterogeneous gene regulation by chimeric transcription factors.

Genotyping-by-sequencing reveals range expansion of Adonis vernalis (Ranunculaceae) from Southeastern Europe into the zonal Euro-Siberian steppe

The Euro-Siberian steppe flora consists of warm- and cold-adapted species, which may have responded differently to Pleistocene glacials and interglacials. Genotyping-by-sequencing individuals from across the distribution range of the pheasant's eye (Adonis vernalis), we aimed to gain insight into steppe florogenesis based on the species' evolutionary history. Although the primary area of origin of the species group comprising A. vernalis, A. villosa and A. volgensis is in Asia, our results indicate that recent populations of A. vernalis are not of Asian origin but evolved in the southern part of Europe during the Pleistocene, with Spanish populations clearly genetically distinct from the Southeastern European populations. We inferred that A. vernalis migrated eastwards from the sub-Mediterranean forest-steppes of Southeastern Europe into the continental forest-steppe zone. Eastern European populations had the highest private allelic richness, indicating long-term large population sizes in this region. As a thermophilic species, A. vernalis seems unlikely to have survived in the cold deserts of the Last Glacial Maximum in Western Siberia, so this region was likely (re)colonized postglacially. Overall, our results reinforce the importance of identifying the area of origin and the corresponding ecological requirements of steppe plants in order to understand the composition of today's steppe flora.

Genome sequence data of the antagonistic soil-borne yeast Cyberlindnera sargentensis (SHA 17.2)

Cyberlindnera sargentensis strain SHA 17.2, isolated from a Swiss soil sample, exhibited strong antagonistic activity against several plant pathogenic fungi in vitro and was highly competitive against other yeasts in soil. As a basis for identifying the mechanisms underlying its strong antagonistic activity, we have sequenced the genome of C. sargentensis (SHA 17.2) by long- and short read sequencing, de novo assembled them into seven contigs/chromosomes and a mitogenome (total genome size 11.4 Mbp), and annotated 5455 genes. This high-quality genome is the reference for transcriptome and proteome analyses aiming at elucidating the mode of action of C. sargentensis against fungal plant pathogens. It will thus serve as a resource for identifying potential biocontrol genes and performing comparative genomics analyses of yeast genomes.

Andean uplift, drainage basin formation, and the evolution of plants living in fast-flowing aquatic ecosystems in northern South America

Northern South America is a geologically dynamic and species-rich region. Fossil and stratigraphic data show that mountain uplift in the tropical Andes reconfigured river drainages. These landscape changes shaped the evolution of the flora in the region, yet the impacts on aquatic taxa have been overlooked. We explore the role of landscape change on the evolution of plants living strictly in rivers across drainage basins in northern South America by conducting population structure, phylogenetic inference, and divergence-dating analyses for two species in the genus Marathrum (Podostemaceae). Mountain uplift and drainage basin formation isolated populations of M. utile and M. foeniculaceum in northern South America and created barriers to gene flow across river drainages. Sympatric species hybridize and the hybrids show the phenotype of one parental line. We propose that the pattern of divergence of populations reflects the formation of river drainages, which was not complete until < 4.1 million yr ago (Ma). Our study provides a clear picture of the role of landscape change on the evolution of plants living strictly in rivers in northern South America. By shifting the focus to aquatic taxa, we provide a novel perspective on the processes shaping the evolution of the Neotropical flora.

Neither connectivity nor genetic diversity matter in the conservation of a rare fern and a moss on insular erratic boulders

Erratic boulders provide habitat for rock-dwelling species and contribute to the biodiversity of landscapes. In the calcareous Swiss lowlands, siliceous erratic boulders are exclusive habitat islands for the regionally critically endangered fern Asplenium septentrionale, about 20 bryophyte species and numerous lichens. Focusing on island biogeographical processes, we analysed the conservation genomics of A. septentrionale and the moss Hedwigia ciliata on insular erratic boulders in the Swiss lowlands and the adjacent “mainland” in siliceous mountains. We genotyped both species using double digest restriction associated DNA sequencing (ddRAD). For the tetraploid A. septentrionale, abundant identical multilocus genotypes within populations suggested prevalent intragametophytic selfing, and six out of eight boulder populations consisting of a single multilocus genotype each indicated single spore founder events. The genetic structure of A. septentrionale mainland populations coincided with Pleistocene glacial refugia. Four genetic lineages of H. ciliata were identified, and populations consisting of a single multilocus genotype were less common than in A. septentrionale. For both taxa, multilocus genotype diversity on boulders was lower than in mainland populations. The absence of common genetic groups among boulder populations, and the absence of isolation by distance patterns, suggested colonisation of boulders through independent long-distance dispersal events. Successful boulder colonisation of A. septentrionale seems to be rare, while colonisation by H. ciliata appears to be more frequent. We conclude that pivotal principles of conservation biology, such as connectivity and genetic diversity, are of less importance for the studied cryptogams on insular erratic boulders because of long-distance dispersal, intragametophytic selfing and polyploidy.

Analysis of Aneuploidy Spectrum From Whole-Genome Sequencing Provides Rapid Assessment of Clonal Variation Within Established Cancer Cell Lines

BACKGROUND

The revolution in next-generation sequencing (NGS) technology has allowed easy access and sharing of high-throughput sequencing datasets of cancer cell lines and their integrative analyses. However, long-term passaging and culture conditions introduce high levels of genomic and phenotypic diversity in established cell lines resulting in strain differences. Thus, clonal variation in cultured cell lines with respect to the reference standard is a major barrier in systems biology data analyses. Therefore, there is a pressing need for a fast and entry-level assessment of clonal variations within cell lines using their high-throughput sequencing data.

RESULTS

We developed a Python-based software, AStra, for de novo estimation of the genome-wide segmental aneuploidy to measure and visually interpret strain-level similarities or differences of cancer cell lines from whole-genome sequencing (WGS). We demonstrated that aneuploidy spectrum can capture the genetic variations in 27 strains of MCF7 breast cancer cell line collected from different laboratories. Performance evaluation of AStra using several cancer sequencing datasets revealed that cancer cell lines exhibit distinct aneuploidy spectra which reflect their previously-reported karyotypic observations. Similarly, AStra successfully identified large-scale DNA copy number variations (CNVs) artificially introduced in simulated WGS datasets.

CONCLUSIONS

AStra provides an analytical and visualization platform for rapid and easy comparison between different strains or between cell lines based on their aneuploidy spectra solely using the raw BAM files representing mapped reads. We recommend AStra for rapid first-pass quality assessment of cancer cell lines before integrating scientific datasets that employ deep sequencing. AStra is an open-source software and is available at https://github.com/AISKhalil/AStra.

Genomic analysis finds no evidence of canonical eukaryotic DNA processing complexes in a free-living protist

Cells replicate and segregate their DNA with precision. Previous studies showed that these regulated cell-cycle processes were present in the last eukaryotic common ancestor and that their core molecular parts are conserved across eukaryotes. However, some metamonad parasites have secondarily lost components of the DNA processing and segregation apparatuses. To clarify the evolutionary history of these systems in these unusual eukaryotes, we generated a genome assembly for the free-living metamonad Carpediemonas membranifera and carried out a comparative genomics analysis. Here, we show that parasitic and free-living metamonads harbor an incomplete set of proteins for processing and segregating DNA. Unexpectedly, Carpediemonas species are further streamlined, lacking the origin recognition complex, Cdc6 and most structural kinetochore subunits. Carpediemonas species are thus the first known eukaryotes that appear to lack this suite of conserved complexes, suggesting that they likely rely on yet-to-be-discovered or alternative mechanisms to carry out these fundamental processes.

Unveiling the Patterns of Reticulated Evolutionary Processes with Phylogenomics: Hybridization and Polyploidy in the genus Rosa

Reticulation, caused by hybridization and allopolyploidization, is considered an important and frequent phenomenon in the evolution of numerous plant lineages. Although both processes represent important driving forces of evolution, they are mostly ignored in phylogenetic studies involving a large number of species. Indeed only a scattering of methods exists to recover a comprehensive reticulated evolutionary history for a broad taxon sampling. Among these methods, comparisons of topologies obtained from plastid markers with those from a few nuclear sequences are favored, even though they restrict in-depth studies of hybridization and polyploidization. The genus Rosa encompasses c. 150 species widely distributed throughout the northern hemisphere and represents a challenging taxonomic group in which hybridization and polyploidization are prominent. Our main objective was to develop a general framework that would take patterns of reticulation into account in the study of the phylogenetic relationships among Rosa species. Using amplicon sequencing we targeted allele variation in the nuclear genome as well as haploid sequences in the chloroplast genome. We successfully recovered robust plastid and nuclear phylogenies and performed in-depth tests for several scenarios of hybridization using a maximum pseudo-likelihood approach on taxon subsets. Our diploid-first approach followed by hybrid and polyploid grafting resolved most of the evolutionary relationships among Rosa subgenera, sections, and selected species. Based on these results, we provide new directions for a future revision of the infrageneric classification in Rosa. The stepwise strategy proposed here can be used to reconstruct the phylogenetic relationships of other challenging taxonomic groups with large numbers of hybrid and polyploid taxa.

Antarctic Thraustochytrids as Sources of Carotenoids and High-Value Fatty Acids

Eicosapentaenoic acid (EPA), docosahexaenoic acid (DHA), and carotenoids are needed as human dietary supplements and are essential components in commercial feeds for the production of aquacultured seafood. Microorganisms such as thraustochytrids are potential natural sources of these compounds. This research reports on the lipid and carotenoid production capacity of thraustochytrids that were isolated from coastal waters of Antarctica. Of the 22 isolates, 21 produced lipids containing EPA+DHA, and the amount of these fatty acids exceeded 20% of the total fatty acids in 12 isolates. Ten isolates were shown to produce carotenoids (27.4-63.9 μg/g dry biomass). The isolate RT2316-16, identified as Thraustochytrium sp., was the best producer of biomass (7.2 g/L in five days) rich in carotenoids (63.9 μg/g) and, therefore, became the focus of this investigation. The main carotenoids in RT2316-16 were β-carotene and canthaxanthin. The content of EPA+DHA in the total lipids (34 ± 3% w/w in dry biomass) depended on the stage of growth of RT2316-16. Lipid and carotenoid content of the biomass and its concentration could be enhanced by modifying the composition of the culture medium. The estimated genome size of RT2316-16 was 44 Mb. Of the 5656 genes predicted from the genome, 4559 were annotated. These included genes of most of the enzymes in the elongation and desaturation pathway of synthesis of ω-3 polyunsaturated fatty acids. Carotenoid precursors in RT2316-16 were synthesized through the mevalonate pathway. A β-carotene synthase gene, with a different domain organization compared to the gene in other thraustochytrids, explained the carotenoid profile of RT2316-16.

Genome, transcriptome and secretome analyses of the antagonistic, yeast-like fungus Aureobasidium pullulans to identify potential biocontrol genes

Aureobasidium pullulans is an extremotolerant, cosmopolitan yeast-like fungus that successfully colonises vastly different ecological niches. The species is widely used in biotechnology and successfully applied as a commercial biocontrol agent against postharvest diseases and fireblight. However, the exact mechanisms that are responsible for its antagonistic activity against diverse plant pathogens are not known at the molecular level. Thus, it is difficult to optimise and improve the biocontrol applications of this species. As a foundation for elucidating biocontrol mechanisms, we have de novo assembled a high-quality reference genome of a strongly antagonistic A. pullulans strain, performed dual RNA-seq experiments, and analysed proteins secreted during the interaction with the plant pathogen Fusarium oxysporum. Based on the genome annotation, potential biocontrol genes were predicted to encode secreted hydrolases or to be part of secondary metabolite clusters (e.g., NRPS-like, NRPS, T1PKS, terpene, and β-lactone clusters). Transcriptome and secretome analyses defined a subset of 79 A. pullulans genes (among the 10,925 annotated genes) that were transcriptionally upregulated or exclusively detected at the protein level during the competition with F. oxysporum. These potential biocontrol genes comprised predicted secreted hydrolases such as glycosylases, esterases, and proteases, as well as genes encoding enzymes, which are predicted to be involved in the synthesis of secondary metabolites. This study highlights the value of a sequential approach starting with genome mining and consecutive transcriptome and secretome analyses in order to identify a limited number of potential target genes for detailed, functional analyses.

Re-examination of two diatom reference genomes using long-read sequencing

Background

The marine diatoms Thalassiosira pseudonana and Phaeodactylum tricornutum are valuable model organisms for exploring the evolution, diversity and ecology of this important algal group. Their reference genomes, published in 2004 and 2008, respectively, were the product of traditional Sanger sequencing. In the case of T. pseudonana, optical restriction site mapping was employed to further clarify and contextualize chromosome-level scaffolds. While both genomes are considered highly accurate and reasonably contiguous, they still contain many unresolved regions and unordered/unlinked scaffolds.

Results

We have used Oxford Nanopore Technologies long-read sequencing to update and validate the quality and contiguity of the T. pseudonana and P. tricornutum genomes. Fine-scale assessment of our long-read derived genome assemblies allowed us to resolve previously uncertain genomic regions, further characterize complex structural variation, and re-evaluate the repetitive DNA content of both genomes. We also identified 1862 previously undescribed genes in T. pseudonana. In P. tricornutum, we used transposable element detection software to identify 33 novel copia-type LTR-RT insertions, indicating ongoing activity and rapid expansion of this superfamily as the organism continues to be maintained in culture. Finally, Bionano optical mapping of P. tricornutum chromosomes was combined with long-read sequence data to explore the potential of long-read sequencing and optical mapping for resolving haplotypes.

Conclusion

Despite its potential to yield highly contiguous scaffolds, long-read sequencing is not a panacea. Even for relatively small nuclear genomes such as those investigated herein, repetitive DNA sequences cause problems for current genome assembly algorithms. Determining whether a long-read derived genomic assembly is ‘better’ than one produced using traditional sequence data is not straightforward. Our revised reference genomes for P. tricornutum and T. pseudonana nevertheless provide additional insight into the structure and evolution of both genomes, thereby providing a more robust foundation for future diatom research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07666-3.

The phylogeographic history of Krascheninnikovia reflects the development of dry steppes and semi-deserts in Eurasia

Constituting one of Earth's major biomes, steppes are characterised by naturally treeless extra-tropical vegetation. The formation of the Eurasian steppe belt, the largest steppe region in the world, began in Central Asia during the Neogene. In the glacial stages of the Pleistocene, steppe displaced forest vegetation, which in turn recolonised the area during the warmer interglacial periods, thus affecting the distribution of plants adapted to these habitats. Krascheninnikovia ceratoides (Chenopodiaceae) is a plant characteristic of dry steppe and semi-desert formations. Earlier studies showed that the ancestor of this autochthonous steppe element originated in Central Asia during the Miocene/Pliocene, i.e., in the same region and at the same time as the first appearance of steppe vegetation. However, as the extant lineages of Krascheninnikovia ceratoides diversified only 2.2 ± 0.9 Mya, it may represent a modern element of current dry steppe and semi-desert formations, rather than a component of the first steppe precursors of the Miocene. As such, it may have capitalised on the climatic conditions of the cold stages of the Quaternary to expand its range and colonise suitable habitats outside of its area of origin. To test this hypothesis, phylogeographic methods were applied to high-resolution genotyping-by-sequencing data. Our results indicate that Krascheninnikovia originated in western Central Asia and the Russian Altai, then spread to Europe in the West, and reached North America in the East. The populations of eastern Central Asia and North America belong to the same clade and are genetically clearly distinct from the Euro-Siberian populations. Among the populations west of the Altai Mountains, the European populations are genetically distinct from all others, which could be the result of the separation of populations east and west of the Urals caused by the Pleistocene transgressions of the Caspian Sea.

Estimating the ploidy of Gracilariopsis lemaneiformis at both the cellular and genomic level1

The determination of the ploidy level of an organism is a prerequisite for studies of evolution, cellular function, and genomic construction. Identification of the ploidy of the economically important red alga Gracilariopsis lemaneiformis has been hindered by its small genome and large number of chromosomes. Therefore, in the current study, PloidyNGS, a tool that calculates the number of reads supporting different alleles at each position along the genome sequence, and fluorescence in situ hybridization coupled with tyramide signal amplification (TSA-FISH) were used to clarify the ploidy of G. lemaneiformis. In addition, flow cytometry (FCM) was used to estimate the ploidy of different somatic cells. The PloidyNGS results showed that most of the alleles in the gametophyte were monomorphic, whereas the TSA-FISH results showed that one hybridization signal was observed in gametophytic nuclei and two in tetrasporophytic nuclei when the nuclei were hybridized by single copy gene probes. These results confirmed that G. lemaneiformis is a haploid in the gametophytic generation and diploid in the sporophytic generation. Moreover, the FCM result suggested that G. lemaneiformis was not an endopolyploid. Based on previous studies, we hypothesize that the nuclear number is important for the cellular differentiation and function of this species. We also suggest that G. lemaneiformis evolved from a paleopolyploid, the genome of which has been diploidized, and that traces of genomic doubling are no longer apparent. Thus, this study provides important evidence for further studies on the evolution and genomes of red algae.

Comparative genomics confirms a rare melioidosis human-to-human transmission event and reveals incorrect phylogenomic reconstruction due to polyclonality

Human-to-human transmission of the melioidosis bacterium, Burkholderia pseudomallei, is exceedingly rare, with only a handful of suspected cases documented to date. Here, we used whole-genome sequencing (WGS) to characterize one such unusual B. pseudomallei transmission event, which occurred between a breastfeeding mother with mastitis and her child. Two strains corresponding to multilocus sequence types (STs)-259 and -261 were identified in the mother's sputum from both the primary culture sweep and in purified colonies, confirming an unusual polyclonal infection in this patient. In contrast, primary culture sweeps of the mother's breast milk and the child's cerebrospinal fluid and blood samples contained only ST-259, indicating monoclonal transmission to the child. Analysis of purified ST-259 isolates showed no genetic variation between mother and baby isolates, providing the strongest possible evidence of B. pseudomallei human-to-human transmission, probably via breastfeeding. Next, phylogenomic analysis of all isolates, including the mother's mixed ST-259/ST-261 sputum sample, was performed to investigate the effects of mixtures on phylogenetic inference. Inclusion of this mixture caused a dramatic reduction in the number of informative SNPs, resulting in branch collapse of ST-259 and ST-261 isolates, and several instances of incorrect topology in a global B. pseudomallei phylogeny, resulting in phylogenetic incongruence. Although phylogenomics can provide clues about the presence of mixtures within WGS datasets, our results demonstrate that this methodology can lead to phylogenetic misinterpretation if mixed genomes are not correctly identified and omitted. Using current bioinformatic tools, we demonstrate a robust method for bacterial mixture identification and strain parsing that avoids these pitfalls.

A Target Capture-Based Method to Estimate Ploidy From Herbarium Specimens

Whole genome duplication (WGD) events are common in many plant lineages, but the ploidy status and possible occurrence of intraspecific ploidy variation are unknown for most species. Standard methods for ploidy determination are chromosome counting and flow cytometry approaches. While flow cytometry approaches typically use fresh tissue, an increasing number of studies have shown that recently dried specimens can be used to yield ploidy data. Recent studies have started to explore whether high-throughput sequencing (HTS) data can be used to assess ploidy levels by analyzing allelic frequencies from single copy nuclear genes. Here, we compare different approaches using a range of yam (Dioscorea) tissues of varying ages, drying methods and quality, including herbarium tissue. Our aims were to: (1) explore the limits of flow cytometry in estimating ploidy level from dried samples, including herbarium vouchers collected between 1831 and 2011, and (2) optimize a HTS-based method to estimate ploidy by considering allelic frequencies from nuclear genes obtained using a target-capture method. We show that, although flow cytometry can be used to estimate ploidy levels from herbarium specimens collected up to fifteen years ago, success rate is low (5.9%). We validated our HTS-based estimates of ploidy using 260 genes by benchmarking with dried samples of species of known ploidy (Dioscorea alata, D. communis, and D. sylvatica). Subsequently, we successfully applied the method to the 85 herbarium samples analyzed with flow cytometry, and successfully provided results for 91.7% of them, comprising species across the phylogenetic tree of Dioscorea. We also explored the limits of using this HTS-based approach for identifying high ploidy levels in herbarium material and the effects of heterozygosity and sequence coverage. Overall, we demonstrated that ploidy diversity within and between species may be ascertained from historical collections, allowing the determination of polyploidization events from samples collected up to two centuries ago. This approach has the potential to provide insights into the drivers and dynamics of ploidy level changes during plant evolution and crop domestication.

Genome-wide sequencing and metabolic annotation of Pythium irregulare CBS 494.86: understanding Eicosapentaenoic acid production

BACKGROUND

Pythium irregulare is an oleaginous Oomycete able to accumulate large amounts of lipids, including Eicosapentaenoic acid (EPA). EPA is an important and expensive dietary supplement with a promising and very competitive market, which is dependent on fish-oil extraction. This has prompted several research groups to study biotechnological routes to obtain specific fatty acids rather than a mixture of various lipids. Moreover, microorganisms can use low cost carbon sources for lipid production, thus reducing production costs. Previous studies have highlighted the production of EPA by P. irregulare, exploiting diverse low cost carbon sources that are produced in large amounts, such as vinasse, glycerol, and food wastewater. However, there is still a lack of knowledge about its biosynthetic pathways, because no functional annotation of any Pythium sp. exists yet. The goal of this work was to identify key genes and pathways related to EPA biosynthesis, in P. irregulare CBS 494.86, by sequencing and performing an unprecedented annotation of its genome, considering the possibility of using wastewater as a carbon source.

RESULTS

Genome sequencing provided 17,727 candidate genes, with 3809 of them associated with enzyme code and 945 with membrane transporter proteins. The functional annotation was compared with curated information of oleaginous organisms, understanding amino acids and fatty acids production, and consumption of carbon and nitrogen sources, present in the wastewater. The main features include the presence of genes related to the consumption of several sugars and candidate genes of unsaturated fatty acids production.

CONCLUSIONS

The whole metabolic genome presented, which is an unprecedented reconstruction of P. irregulare CBS 494.86, shows its potential to produce value-added products, in special EPA, for food and pharmaceutical industries, moreover it infers metabolic capabilities of the microorganism by incorporating information obtained from literature and genomic data, supplying information of great importance to future work.

Extensive loss of cell-cycle and DNA repair genes in an ancient lineage of bipolar budding yeasts

Cell-cycle checkpoints and DNA repair processes protect organisms from potentially lethal mutational damage. Compared to other budding yeasts in the subphylum Saccharomycotina, we noticed that a lineage in the genus Hanseniaspora exhibited very high evolutionary rates, low Guanine-Cytosine (GC) content, small genome sizes, and lower gene numbers. To better understand Hanseniaspora evolution, we analyzed 25 genomes, including 11 newly sequenced, representing 18/21 known species in the genus. Our phylogenomic analyses identify two Hanseniaspora lineages, a faster-evolving lineage (FEL), which began diversifying approximately 87 million years ago (mya), and a slower-evolving lineage (SEL), which began diversifying approximately 54 mya. Remarkably, both lineages lost genes associated with the cell cycle and genome integrity, but these losses were greater in the FEL. E.g., all species lost the cell-cycle regulator WHIskey 5 (WHI5), and the FEL lost components of the spindle checkpoint pathway (e.g., Mitotic Arrest-Deficient 1 [MAD1], Mitotic Arrest-Deficient 2 [MAD2]) and DNA-damage-checkpoint pathway (e.g., Mitosis Entry Checkpoint 3 [MEC3], RADiation sensitive 9 [RAD9]). Similarly, both lineages lost genes involved in DNA repair pathways, including the DNA glycosylase gene 3-MethylAdenine DNA Glycosylase 1 (MAG1), which is part of the base-excision repair pathway, and the DNA photolyase gene PHotoreactivation Repair deficient 1 (PHR1), which is involved in pyrimidine dimer repair. Strikingly, the FEL lost 33 additional genes, including polymerases (i.e., POLymerase 4 [POL4] and POL32) and telomere-associated genes (e.g., Repressor/activator site binding protein-Interacting Factor 1 [RIF1], Replication Factor A 3 [RFA3], Cell Division Cycle 13 [CDC13], Pbp1p Binding Protein [PBP2]). Echoing these losses, molecular evolutionary analyses reveal that, compared to the SEL, the FEL stem lineage underwent a burst of accelerated evolution, which resulted in greater mutational loads, homopolymer instabilities, and higher fractions of mutations associated with the common endogenously damaged base, 8-oxoguanine. We conclude that Hanseniaspora is an ancient lineage that has diversified and thrived, despite lacking many otherwise highly conserved cell-cycle and genome integrity genes and pathways, and may represent a novel, to our knowledge, system for studying cellular life without them.

Ploidy Determination in the Pathogenic Fungus Sporothrix spp

The pathogenic clade of the Sporothrix genus comprises the etiological agents of sporotrichosis, a worldwide emergent disease. Despite the growing understanding of their successful pathogen traits, there is little information on genome sizes and ploidy within the genus. Therefore, in this work, we evaluated the ploidy of four species of the Sporothrix genus, specifically Sporothrix brasiliensis, Sporothrix schenckii, Sporothrix globosa, and Sporothrix pallida. Through cell cycle analysis of the yeast-phase cells, we showed that the DNA content of G0/G1 cells was similar to the genome size determined by whole genome sequencing. Moreover, ploidy of S. schenckii, S. brasiliensis, and S. pallida that was determined by allele composition using next-generation sequencing (NGS) data is consistent with monomorphic positions at each allele. These data show that the analyzed strains of Sporothrix are haploid, or at least aneuploid, thereby laying the foundation for the development of a molecular toolbox for Sporothrix spp.

Draft Genome Sequence for the Tree Pathogen Phytophthora plurivora

Species from the genus Phytophthora are well represented among organisms causing serious diseases on trees. Phytophthora plurivora has been implicated in long-term decline of woodland trees across Europe. Here we present a draft genome sequence of P. plurivora, originally isolated from diseased European beech (Fagus sylvatica) in Malmö, Sweden. When compared with other sequenced Phytophthora species, the P. plurivora genome assembly is relatively compact, spanning 41 Mb. This is organized in 1,919 contigs and 1,898 scaffolds, encompassing 11,741 predicted genes, and has a repeat content of approximately 15%. Comparison of allele frequencies revealed evidence for tetraploidy in the sequenced isolate. As in other sequenced Phytophthora species, P. plurivora possesses genes for pathogenicity-associated RXLR and Crinkle and Necrosis effectors, predominantly located in gene-sparse genomic regions. Comparison of the P. plurivora RXLR effectors with orthologs in other sequenced species in the same clade (Phytophthora multivora and Phytophthora capsici) revealed that the orthologs were likely to be under neutral or purifying selection.

nQuire: a statistical framework for ploidy estimation using next generation sequencing

Background

Intraspecific variation in ploidy occurs in a wide range of species including pathogenic and nonpathogenic eukaryotes such as yeasts and oomycetes. Ploidy can be inferred indirectly - without measuring DNA content - from experiments using next-generation sequencing (NGS). We present nQuire, a statistical framework that distinguishes between diploids, triploids and tetraploids using NGS. The command-line tool models the distribution of base frequencies at variable sites using a Gaussian Mixture Model, and uses maximum likelihood to select the most plausible ploidy model. nQuire handles large genomes at high coverage efficiently and uses standard input file formats.

Results

We demonstrate the utility of nQuire analyzing individual samples of the pathogenic oomycete Phytophthora infestans and the Baker’s yeast Saccharomyces cerevisiae. Using these organisms we show the dependence between reliability of the ploidy assignment and sequencing depth. Additionally, we employ normalized maximized log- likelihoods generated by nQuire to ascertain ploidy level in a population of samples with ploidy heterogeneity. Using these normalized values we cluster samples in three dimensions using multivariate Gaussian mixtures. The cluster assignments retrieved from a S. cerevisiae population recovered the true ploidy level in over 96% of samples. Finally, we show that nQuire can be used regionally to identify chromosomal aneuploidies.

Conclusions

nQuire provides a statistical framework to study organisms with intraspecific variation in ploidy. nQuire is likely to be useful in epidemiological studies of pathogens, artificial selection experiments, and for historical or ancient samples where intact nuclei are not preserved. It is implemented as a stand-alone Linux command line tool in the C programming language and is available at https://github.com/clwgg/nQuireunder the MIT license.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2128-z) contains supplementary material, which is available to authorized users.

nQuire: a statistical framework for ploidy estimation using next generation sequencing

BACKGROUND

Intraspecific variation in ploidy occurs in a wide range of species including pathogenic and nonpathogenic eukaryotes such as yeasts and oomycetes. Ploidy can be inferred indirectly - without measuring DNA content - from experiments using next-generation sequencing (NGS). We present nQuire, a statistical framework that distinguishes between diploids, triploids and tetraploids using NGS. The command-line tool models the distribution of base frequencies at variable sites using a Gaussian Mixture Model, and uses maximum likelihood to select the most plausible ploidy model. nQuire handles large genomes at high coverage efficiently and uses standard input file formats.

RESULTS

We demonstrate the utility of nQuire analyzing individual samples of the pathogenic oomycete Phytophthora infestans and the Baker's yeast Saccharomyces cerevisiae. Using these organisms we show the dependence between reliability of the ploidy assignment and sequencing depth. Additionally, we employ normalized maximized log- likelihoods generated by nQuire to ascertain ploidy level in a population of samples with ploidy heterogeneity. Using these normalized values we cluster samples in three dimensions using multivariate Gaussian mixtures. The cluster assignments retrieved from a S. cerevisiae population recovered the true ploidy level in over 96% of samples. Finally, we show that nQuire can be used regionally to identify chromosomal aneuploidies.

CONCLUSIONS

nQuire provides a statistical framework to study organisms with intraspecific variation in ploidy. nQuire is likely to be useful in epidemiological studies of pathogens, artificial selection experiments, and for historical or ancient samples where intact nuclei are not preserved. It is implemented as a stand-alone Linux command line tool in the C programming language and is available at https://github.com/clwgg/nQuire under the MIT license.

Current Strategies of Polyploid Plant Genome Sequence Assembly

Polyploidy or duplication of an entire genome occurs in the majority of angiosperms. The understanding of polyploid genomes is important for the improvement of those crops, which humans rely on for sustenance and basic nutrition. As climate change continues to pose a potential threat to agricultural production, there will increasingly be a demand for plant cultivars that can resist biotic and abiotic stresses and also provide needed and improved nutrition. In the past decade, Next Generation Sequencing (NGS) has fundamentally changed the genomics landscape by providing tools for the exploration of polyploid genomes. Here, we review the challenges of the assembly of polyploid plant genomes, and also present recent advances in genomic resources and functional tools in molecular genetics and breeding. As genomes of diploid and less heterozygous progenitor species are increasingly available, we discuss the lack of complexity of these currently available reference genomes as they relate to polyploid crops. Finally, we review recent approaches of haplotyping by phasing and the impact of third generation technologies on polyploid plant genome assembly.

Draft Genome Sequence of Saccharomyces cerevisiae Barra Grande (BG-1), a Brazilian Industrial Bioethanol-Producing Strain

Here, we present the draft genome sequence of Saccharomyces cerevisiae BG-1, a Brazilian industrial strain widely used for bioethanol production from sugarcane. The 11.7-Mb genome sequence consists of 216 scaffolds and harbors 5,607 predicted protein-coding genes.