The Madagascar palm genome provides new insights on the evolution of Apocynaceae specialized metabolism.

Specialized metabolites possess diverse interesting biological activities and some cardenolides- and monoterpene indole alkaloids- (MIAs) derived pharmaceuticals are currently used to treat human diseases such as cancers or hypertension. While these two families of biocompounds are produced by specific subfamilies of Apocynaceae, one member of this medicinal plant family, the succulent tree Pachypodium lamerei Drake (also known as Madagascar palm), does not produce such specialized metabolites. To explore the evolutionary paths that have led to the emergence and loss of cardenolide and MIA biosynthesis in Apocynaceae, we sequenced and assembled the P. lamerei genome by combining Oxford Nanopore Technologies long-reads and Illumina short-reads. Phylogenomics revealed that, among the Apocynaceae whose genomes have been sequenced, the Madagascar palm is so far the species closest to the common ancestor between MIA producers/non-MIA producers. Transposable elements, constituting 72.48% of the genome, emerge as potential key players in shaping genomic architecture and influencing specialized metabolic pathways. The absence of crucial MIA biosynthetic genes such as strictosidine synthase in P. lamerei and non-Rauvolfioideae species hints at a transposon-mediated mechanism behind gene loss. Phylogenetic analysis not only showcases the evolutionary divergence of specialized metabolite biosynthesis within Apocynaceae but also underscores the role of transposable elements in this intricate process. Moreover, we shed light on the low conservation of enzymes involved in the final stages of MIA biosynthesis in the distinct MIA-producing plant families, inferring independent gains of these specialized enzymes along the evolution of these medicinal plant clades. Overall, this study marks a leap forward in understanding the genomic dynamics underpinning the evolution of specialized metabolites biosynthesis in the Apocynaceae family, with transposons emerging as potential architects of genomics restructuring and gene loss.

The Rauvolfia tetraphylla genome suggests multiple distinct biosynthetic routes for yohimbane monoterpene indole alkaloids.

Monoterpene indole alkaloids (MIAs) are a structurally diverse family of specialized metabolites mainly produced in Gentianales to cope with environmental challenges. Due to their pharmacological properties, the biosynthetic modalities of several MIA types have been elucidated but not that of the yohimbanes. Here, we combine metabolomics, proteomics, transcriptomics and genome sequencing of Rauvolfia tetraphylla with machine learning to discover the unexpected multiple actors of this natural product synthesis. We identify a medium chain dehydrogenase/reductase (MDR) that produces a mixture of four diastereomers of yohimbanes including the well-known yohimbine and rauwolscine. In addition to this multifunctional yohimbane synthase (YOS), an MDR synthesizing mainly heteroyohimbanes and the short chain dehydrogenase vitrosamine synthase also display a yohimbane synthase side activity. Lastly, we establish that the combination of geissoschizine synthase with at least three other MDRs also produces a yohimbane mixture thus shedding light on the complex mechanisms evolved for the synthesis of these plant bioactives.

The Vinca minor genome highlights conserved evolutionary traits in monoterpene indole alkaloid synthesis.

Vinca minor, also known as the lesser periwinkle, is a well-known species from the Apocynaceae, native to central and southern Europe. This plant synthesizes monoterpene indole alkaloids, which are a class of specialized metabolites displaying a wide range of bioactive- and pharmacologically important properties. Within the almost 50 monoterpene indole alkaloids it produces, V. minor mainly accumulates vincamine, which is commercially used as a nootropic. Using a combination of Oxford Nanopore Technologies long read- and Illumina short-read sequencing, a 679,098 Mb V. minor genome was assembled into 296 scaffolds with an N50 scaffold length of 6 Mb, and encoding 29,624 genes. These genes were functionally annotated and used in a comparative genomic analysis to establish gene families and to investigate gene family expansion and contraction across the phylogenetic tree. Furthermore, homology-based monoterpene indole alkaloid gene predictions together with a metabolic analysis across 4 different V. minor tissue types guided the identification of candidate monoterpene indole alkaloid genes. These candidates were finally used to identify monoterpene indole alkaloid gene clusters, which combined with synteny analysis allowed for the discovery of a functionally validated vincadifformine-16-hydroxylase, reinforcing the potential of this dataset for monoterpene indole alkaloids gene discovery. It is expected that access to these resources will facilitate the elucidation of unknown monoterpene indole alkaloid biosynthetic routes with the potential of transferring these pathways to heterologous expression systems for large-scale monoterpene indole alkaloid production.

Genome Assembly of the Medicinal Plant Voacanga thouarsii.

The Apocynaceae tree Voacanga thouarsii, native to southern Africa and Madagascar, produces monoterpene indole alkaloids (MIA), which are specialized metabolites with a wide range of bioactive properties. Voacanga species mainly accumulates tabersonine in seeds making these species valuable medicinal plants currently used for industrial MIA production. Despite their importance, the MIA biosynthesis in Voacanga species remains poorly studied. Here, we report the first genome assembly and annotation of a Voacanga species. The combined assembly of Oxford Nanopore Technologies long-reads and Illumina short-reads resulted in 3,406 scaffolds with a total length of 1,354.26 Mb and an N50 of 3.04 Mb. A total of 33,300 protein-coding genes were predicted and functionally annotated. These genes were then used to establish gene families and to investigate gene family expansion and contraction across the phylogenetic tree. A transposable element (TE) analysis showed the highest proportion of TE in Voacanga thouarsii compared with all other MIA-producing plants. In a nutshell, this first reference genome of V. thouarsii will thus contribute to strengthen future comparative and evolutionary studies in MIA-producing plants leading to a better understanding of MIA pathway evolution. This will also allow the potential identification of new MIA biosynthetic genes for metabolic engineering purposes.

An updated version of the Madagascar periwinkle genome.

The Madagascar periwinkle, Catharanthus roseus, belongs to the Apocynaceae family. This medicinal plant, endemic to Madagascar, produces many important drugs including the monoterpene indole alkaloids (MIA) vincristine and vinblastine used to treat cancer worldwide. Here, we provide a new version of the C. roseus genome sequence obtained through the combination of Oxford Nanopore Technologies long-reads and Illumina short-reads. This more contiguous assembly consists of 173 scaffolds with a total length of 581.128 Mb and an N50 of 12.241 Mb. Using publicly available RNAseq data, 21,061 protein coding genes were predicted and functionally annotated. A total of 42.87% of the genome was annotated as transposable elements, most of them being long-terminal repeats. Together with the increasing access to MIA-producing plant genomes, this updated version should ease evolutionary studies leading to a better understanding of MIA biosynthetic pathway evolution.

Genome and transcriptome analysis of the beet armyworm Spodoptera exigua reveals targets for pest control.

The genus Spodoptera (Lepidoptera: Noctuidae) includes some of the most infamous insect pests of cultivated plants including Spodoptera frugiperda, Spodoptera litura, and Spodoptera exigua. To effectively develop targeted pest control strategies for diverse Spodoptera species, genomic resources are highly desired. To this aim, we provide the genome assembly and developmental transcriptome comprising all major life stages of S. exigua, the beet armyworm. Spodoptera exigua is a polyphagous herbivore that can feed on > 130 host plants, including several economically important crops. The 419 Mb beet armyworm genome was sequenced from a female S. exigua pupa. Using a hybrid genome sequencing approach (Nanopore long-read data and Illumina short read), a high-quality genome assembly was achieved (N50 = 1.1 Mb). An official gene set (18,477 transcripts) was generated by automatic annotation and by using transcriptomic RNA-seq datasets of 18 S. exigua samples as supporting evidence. In-depth analyses of developmental stage-specific expression combined with gene tree analyses of identified homologous genes across Lepidoptera genomes revealed four potential genes of interest (three of them Spodoptera-specific) upregulated during first- and third-instar larval stages for targeted pest-outbreak management. The beet armyworm genome sequence and developmental transcriptome covering all major developmental stages provide critical insights into the biology of this devastating polyphagous insect pest species worldwide. In addition, comparative genomic analyses across Lepidoptera significantly advance our knowledge to further control other invasive Spodoptera species and reveals potential lineage-specific target genes for pest control strategies.

Corrigendum: RNAseq Profiling of Leukocyte Populations in Zebrafish Larvae Reveals a cxcl11 Chemokine Gene as a Marker of Macrophage Polarization During Mycobacterial Infection.

[This corrects the article DOI: 10.3389/fimmu.2019.00832.].

RNAseq Profiling of Leukocyte Populations in Zebrafish Larvae Reveals a cxcl11 Chemokine Gene as a Marker of Macrophage Polarization During Mycobacterial Infection.

Macrophages are phagocytic cells from the innate immune system, which forms the first line of host defense against invading pathogens. These highly dynamic immune cells can adopt specific functional phenotypes, with the pro-inflammatory M1 and anti-inflammatory M2 polarization states as the two extremes. Recently, the process of macrophage polarization during inflammation has been visualized by real time imaging in larvae of the zebrafish. This model organism has also become widely used to study macrophage responses to microbial pathogens. To support the increasing use of zebrafish in macrophage biology, we set out to determine the complete transcriptome of zebrafish larval macrophages. We studied the specificity of the macrophage signature compared with other larval immune cells and the macrophage-specific expression changes upon infection. We made use of the well-established mpeg1, mpx, and lck fluorescent reporter lines to sort and sequence the transcriptome of larval macrophages, neutrophils, and lymphoid progenitor cells, respectively. Our results provide a complete dataset of genes expressed in these different immune cell types and highlight their similarities and differences. Major differences between the macrophage and neutrophil signatures were found within the families of proteinases. Furthermore, expression of genes involved in antigen presentation and processing was specifically detected in macrophages, while lymphoid progenitors showed expression of genes involved in macrophage activation. Comparison with datasets of in vitro polarized human macrophages revealed that zebrafish macrophages express a strongly homologous gene set, comprising both M1 and M2 markers. Furthermore, transcriptome analysis of low numbers of macrophages infected by the intracellular pathogen Mycobacterium marinum revealed that infected macrophages change their transcriptomic response by downregulation of M2-associated genes and overexpression of specific M1-associated genes. Among the infection-induced genes, a homolog of the human CXCL11 chemokine gene, cxcl11aa, stood out as the most strongly overexpressed M1 marker. Upregulation of cxcl11aa in Mycobacterium-infected macrophages was found to require the function of Myd88, a critical adaptor molecule in the Toll-like and interleukin 1 receptor pathways that are central to pathogen recognition and activation of the innate immune response. Altogether, our data provide a valuable data mining resource to support infection and inflammation research in the zebrafish model.

Author Correction: Rapid de novo assembly of the European eel genome from nanopore sequencing reads.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Deep transcriptome analysis of the heat shock response in an Atlantic sturgeon (Acipenser oxyrinchus) cell line.

Despite efforts to restore Atlantic sturgeon in European rivers, aquaculture techniques result in animals with high post-release mortality due to, among other reasons, their low tolerance to increasing water temperature. Marker genes to monitor heat stress are needed in order to identify heat-resistant fish. Therefore, an Atlantic sturgeon cell line was exposed to different heat shock protocols (30 °C and 35 °C) and differences in gene expression were investigated. In total 3020 contigs (∼1.5%) were differentially expressed. As the core of the upregulated contigs corresponded to heat shock proteins (HSP), the heat shock factor (HSF) and the HSP gene families were annotated in Atlantic sturgeon and mapped via Illumina RNA sequencing to identify heat-inducible family members. Up to 6 hsf and 76 hsp genes were identified in the Atlantic sturgeon transcriptome resources, 16 of which were significantly responsive to the applied heat shock. The previously studied hspa1 (hsp70) gene was only significantly upregulated at the highest heat shock (35 °C), while a set of 5 genes (hspc1, hsph3a, hspb1b, hspb11a, and hspb11b) was upregulated at all conditions. Although the hspc1 (hsp90a) gene was previously used as heat shock-marker in sturgeons, we found that hspb11a is the most heat-inducible gene, with up to 3296-fold higher expression in the treated cells, constituting the candidate gene markers for in vivo trials.

Rapid de novo assembly of the European eel genome from nanopore sequencing reads.

We have sequenced the genome of the endangered European eel using the MinION by Oxford Nanopore, and assembled these data using a novel algorithm specifically designed for large eukaryotic genomes. For this 860 Mbp genome, the entire computational process takes two days on a single CPU. The resulting genome assembly significantly improves on a previous draft based on short reads only, both in terms of contiguity (N50 1.2 Mbp) and structural quality. This combination of affordable nanopore sequencing and light weight assembly promises to make high-quality genomic resources accessible for many non-model plants and animals.

MinION Analysis and Reference Consortium: Phase 2 data release and analysis of R9.0 chemistry.

BACKGROUND: Long-read sequencing is rapidly evolving and reshaping the suite of opportunities for genomic analysis. For the MinION in particular, as both the platform and chemistry develop, the user community requires reference data to set performance expectations and maximally exploit third-generation sequencing. We performed an analysis of MinION data derived from whole genome sequencing of Escherichiacoli K-12 using the R9.0 chemistry, comparing the results with the older R7.3 chemistry. METHODS: We computed the error-rate estimates for insertions, deletions, and mismatches in MinION reads. RESULTS: Run-time characteristics of the flow cell and run scripts for R9.0 were similar to those observed for R7.3 chemistry, but with an 8-fold increase in bases per second (from 30 bps in R7.3 and SQK-MAP005 library preparation, to 250 bps in R9.0) processed by individual nanopores, and less drop-off in yield over time. The 2-dimensional ("2D") N50 read length was unchanged from the prior chemistry. Using the proportion of alignable reads as a measure of base-call accuracy, 99.9% of "pass" template reads from 1-dimensional ("1D")  experiments were mappable and ~97% from 2D experiments. The median identity of reads was ~89% for 1D and ~94% for 2D experiments. The total error rate (miscall + insertion + deletion ) decreased for 2D "pass" reads from 9.1% in R7.3 to 7.5% in R9.0 and for template "pass" reads from 26.7% in R7.3 to 14.5% in R9.0. CONCLUSIONS: These Phase 2 MinION experiments serve as a baseline by providing estimates for read quality, throughput, and mappability. The datasets further enable the development of bioinformatic tools tailored to the new R9.0 chemistry and the design of novel biological applications for this technology. ABBREVIATIONS: K: thousand, Kb: kilobase (one thousand base pairs), M: million, Mb: megabase (one million base pairs), Gb: gigabase (one billion base pairs).

De novo whole-genome assembly of a wild type yeast isolate using nanopore sequencing.

Background: The introduction of the MinION TM sequencing device by Oxford Nanopore Technologies may greatly accelerate whole genome sequencing. It has been shown that the nanopore sequence data, in combination with other sequencing technologies, is highly useful for accurate annotation of all genes in the genome. However, it also offers great potential for de novo assembly of complex genomes without using other technologies. In this manuscript we used nanopore sequencing as a tool to classify yeast strains. Methods: We compared various technical and software developments for the nanopore sequencing protocol, showing that the R9 chemistry is, as predicted, higher in quality than R7.3 chemistry. The R9 chemistry is an essential improvement for assembly of the extremely AT-rich mitochondrial genome. Results: In this study, we used this new technology to sequence and de novo assemble the genome of a recently isolated ethanologenic yeast strain, and compared the results with those obtained by classical Illumina short read sequencing. This strain was originally named Candida vartiovaarae ( Torulopsis vartiovaarae) based on ribosomal RNA sequencing. We show that the assembly using nanopore data is much more contiguous than the assembly using short read data. Conclusions: The mitochondrial and chromosomal genome sequences showed that our strain is clearly distinct from other yeast taxons and most closely related to published Cyberlindnera species. In conclusion, MinION-mediated long read sequencing can be used for high quality de novo assembly of new eukaryotic microbial genomes.

Changes in ovarian gene expression profiles and plasma hormone levels in maturing European eel (Anguilla anguilla); Biomarkers for broodstock selection.

Complete sexual maturation of European eels (Anguilla anguilla) in captivity can only be achieved via injections with gonadotropins. For female eels this procedure takes 4-6months and the response ranges from "unresponsive" to final maturation and ovulation. Reproductive success could be significantly increased via early selection of responders based on predictive markers and minimally invasive sampling methods. To get a better understanding of the genetic background of ovarian maturation of the European eel we performed a pilot deep-sequencing transcriptome analysis of ovarian tissue derived from a yellow eel, a prepubertal silver eel and a post-spawning matured eel. Two key players in steroidogenesis were strongly correlated with advanced sexual maturation, namely P450c17 and liver receptor homolog-1, suggesting that blood plasma steroids might qualify as minimally invasive markers for early detection of responders. Since the predictive value of plasma sex steroid levels for final maturation of the European eel had not yet been carefully examined, we performed an extensive artificial maturation trial. Farmed silver eels were treated with pituitary extracts and sampled at multiple time intervals. Expression of steroidogenesis-related genes in ovarian tissue of responding and non-responding eels after four weekly injections with pituitary extract was compared using a custom-built microarray and RNAseq. Increased expression of 17ß-hsd1 was strongly linked to sexual maturation. Blood plasma levels of sex steroids were measured using ELISAs. We show that a 2.5-fold increase in blood-plasma estradiol level after 4 weekly pituitary extract injections is a strong predictor of final sexual maturation of female European eel.

Testing tuberculosis drug efficacy in a zebrafish high-throughput translational medicine screen.

The translational value of zebrafish high-throughput screens can be improved when more knowledge is available on uptake characteristics of potential drugs. We investigated reference antibiotics and 15 preclinical compounds in a translational zebrafish-rodent screening system for tuberculosis. As a major advance, we have developed a new tool for testing drug uptake in the zebrafish model. This is important, because despite the many applications of assessing drug efficacy in zebrafish research, the current methods for measuring uptake using mass spectrometry do not take into account the possible adherence of drugs to the larval surface. Our approach combines nanoliter sampling from the yolk using a microneedle, followed by mass spectrometric analysis. To date, no single physicochemical property has been identified to accurately predict compound uptake; our method offers a great possibility to monitor how any novel compound behaves within the system. We have correlated the uptake data with high-throughput drug-screening data from Mycobacterium marinum-infected zebrafish larvae. As a result, we present an improved zebrafish larva drug-screening platform which offers new insights into drug efficacy and identifies potential false negatives and drugs that are effective in zebrafish and rodents. We demonstrate that this improved zebrafish drug-screening platform can complement conventional models of in vivo Mycobacterium tuberculosis-infected rodent assays. The detailed comparison of two vertebrate systems, fish and rodent, may give more predictive value for efficacy of drugs in humans.

MinION Analysis and Reference Consortium: Phase 1 data release and analysis.

The advent of a miniaturized DNA sequencing device with a high-throughput contextual sequencing capability embodies the next generation of large scale sequencing tools. The MinION™ Access Programme (MAP) was initiated by Oxford Nanopore Technologies™ in April 2014, giving public access to their USB-attached miniature sequencing device. The MinION Analysis and Reference Consortium (MARC) was formed by a subset of MAP participants, with the aim of evaluating and providing standard protocols and reference data to the community. Envisaged as a multi-phased project, this study provides the global community with the Phase 1 data from MARC, where the reproducibility of the performance of the MinION was evaluated at multiple sites. Five laboratories on two continents generated data using a control strain of Escherichia coli K-12, preparing and sequencing samples according to a revised ONT protocol. Here, we provide the details of the protocol used, along with a preliminary analysis of the characteristics of typical runs including the consistency, rate, volume and quality of data produced. Further analysis of the Phase 1 data presented here, and additional experiments in Phase 2 of E. coli from MARC are already underway to identify ways to improve and enhance MinION performance.

RNA sequencing of FACS-sorted immune cell populations from zebrafish infection models to identify cell specific responses to intracellular pathogens.

The zebrafish (Danio rerio) is increasingly used as a model for studying infectious diseases. This nonmammalian vertebrate host, which is transparent at the early life stages, is especially attractive for live imaging of interactions between pathogens and host cells. A number of useful fluorescent reporter lines have recently been developed and significant advances in RNA sequencing technology have been made, which now make it possible to apply the zebrafish model for investigating changes in transcriptional activity of specific immune cell types during the course of an infection process.Here we describe how to sequence RNA extracted from fluorescently labeled macrophages obtained by cell-sorting of 5-day-old zebrafish larvae of the transgenic Tg(mpeg1:Gal4-VP16);Tg(UAS-E1b:Kaede) line. This technique showed reproducible results and allowed to detect specific expression of macrophage markers in the mpeg1 positive cell population, whereas no markers specific for neutrophils or lymphoid cells were detected. This protocol has been also successfully extended to other immune cell types as well as cells infected by Mycobacterium marinum.

Advances in genomics of bony fish.

In this review, we present an overview of the recent advances of genomic technologies applied to studies of fish species belonging to the superclass of Osteichthyes (bony fish) with a major emphasis on the infraclass of Teleostei, also called teleosts. This superclass that represents more than 50% of all known vertebrate species has gained considerable attention from genome researchers in the last decade. We discuss many examples that demonstrate that this highly deserved attention is currently leading to new opportunities for answering important biological questions on gene function and evolutionary processes. In addition to giving an overview of the technologies that have been applied for studying various fish species we put the recent advances in genome research on the model species zebrafish and medaka in the context of its impact for studies of all fish of the superclass of Osteichthyes. We thereby want to illustrate how the combined value of research on model species together with a broad angle perspective on all bony fish species will have a huge impact on research in all fields of fundamental science and will speed up applications in many societally important areas such as the development of new medicines, toxicology test systems, environmental sensing systems and sustainable aquaculture strategies.

The king cobra genome reveals dynamic gene evolution and adaptation in the snake venom system.

Snakes are limbless predators, and many species use venom to help overpower relatively large, agile prey. Snake venoms are complex protein mixtures encoded by several multilocus gene families that function synergistically to cause incapacitation. To examine venom evolution, we sequenced and interrogated the genome of a venomous snake, the king cobra (Ophiophagus hannah), and compared it, together with our unique transcriptome, microRNA, and proteome datasets from this species, with data from other vertebrates. In contrast to the platypus, the only other venomous vertebrate with a sequenced genome, we find that snake toxin genes evolve through several distinct co-option mechanisms and exhibit surprisingly variable levels of gene duplication and directional selection that correlate with their functional importance in prey capture. The enigmatic accessory venom gland shows a very different pattern of toxin gene expression from the main venom gland and seems to have recruited toxin-like lectin genes repeatedly for new nontoxic functions. In addition, tissue-specific microRNA analyses suggested the co-option of core genetic regulatory components of the venom secretory system from a pancreatic origin. Although the king cobra is limbless, we recovered coding sequences for all Hox genes involved in amniote limb development, with the exception of Hoxd12. Our results provide a unique view of the origin and evolution of snake venom and reveal multiple genome-level adaptive responses to natural selection in this complex biological weapon system. More generally, they provide insight into mechanisms of protein evolution under strong selection.