A multi-omic Nicotiana benthamiana resource for fundamental research and biotechnology.

Nicotiana benthamiana is an invaluable model plant and biotechnology platform with a ~3 Gb allotetraploid genome. To further improve its usefulness and versatility, we have produced high-quality chromosome-level genome assemblies, coupled with transcriptome, epigenome, microRNA and transposable element datasets, for the ubiquitously used LAB strain and a related wild accession, QLD. In addition, single nucleotide polymorphism maps have been produced for a further two laboratory strains and four wild accessions. Despite the loss of five chromosomes from the ancestral tetraploid, expansion of intergenic regions, widespread segmental allopolyploidy, advanced diploidization and evidence of recent bursts of Copia pseudovirus (Copia) mobility not seen in other Nicotiana genomes, the two subgenomes of N. benthamiana show large regions of synteny across the Solanaceae. LAB and QLD have many genetic, metabolic and phenotypic differences, including disparate RNA interference responses, but are highly interfertile and amenable to genome editing and both transient and stable transformation. The LAB/QLD combination has the potential to be as useful as the Columbia-0/Landsberg errecta partnership, utilized from the early pioneering days of Arabidopsis genomics to today.

A practical comparison of the next-generation sequencing platform and assemblers using yeast genome.

Assembling fragmented whole-genomic information from the sequencing data is an inevitable process for further genome-wide research. However, it is intricate to select the appropriate assembly pipeline for unknown species because of the species-specific genomic properties. Therefore, our study focused on relatively more static proclivities of sequencing platforms and assembly algorithms than the fickle genome sequences. A total of 212 draft and polished de novo assemblies were constructed under the different sequencing platforms and assembly algorithms with the repetitive yeast genome. Our comprehensive data indicated that sequencing reads from Oxford Nanopore with R7.3 flow cells generated more continuous assemblies than those derived from the PacBio Sequel, although the homopolymer-based assembly errors and chimeric contigs exist. In addition, the comparison between two second-generation sequencing platforms showed that Illumina NovaSeq 6000 provides more accurate and continuous assembly in the second-generation-sequencing-first pipeline, but MGI DNBSEQ-T7 provides a cheap and accurate read in the polishing process. Furthermore, our insight into the relationship among the computational time, read length, and coverage depth provided clues to the optimal pipelines of yeast assembly.

Transcriptome profiling of olive flounder responses under acute and chronic heat stress.

BACKGROUND: The olive flounder (Paralichthys olivaceus) is a saltwater fish, which is valuable to the economy. The olive flounder strives to adapt to environmental stressors through physiological, biochemical, and transcriptional responses. The rise in water temperature threatens the growth, development, reproduction, and survival of olive flounder. Each organ in the olive flounder can differentially respond to heat stress. OBJECTIVES: The purpose of this study is to investigate organ-specific transcriptional changes in olive flounder tissues during heat stress. METHODS: In this study, transcriptome dynamics of the gill, liver, and muscle of olive flounder to acute or chronic heat stress were investigated. RESULTS: Principal component analysis plotting revealed that the transcriptome of each organ is quite separated. K-means clustering, gene ontology, and Kyoto Encyclopedia of Genes and Genomes pathway analysis showed the differential transcriptome responses of each organ to heat stress. Heat stress commonly affects the pathways involved in the correct protein folding, DNA repair, and cell cycle. CONCLUSION: Our results may provide a valuable molecular basis of heat acclimation in fishes.

Twelve quick steps for genome assembly and annotation in the classroom.

Eukaryotic genome sequencing and de novo assembly, once the exclusive domain of well-funded international consortia, have become increasingly affordable, thus fitting the budgets of individual research groups. Third-generation long-read DNA sequencing technologies are increasingly used, providing extensive genomic toolkits that were once reserved for a few select model organisms. Generating high-quality genome assemblies and annotations for many aquatic species still presents significant challenges due to their large genome sizes, complexity, and high chromosome numbers. Indeed, selecting the most appropriate sequencing and software platforms and annotation pipelines for a new genome project can be daunting because tools often only work in limited contexts. In genomics, generating a high-quality genome assembly/annotation has become an indispensable tool for better understanding the biology of any species. Herein, we state 12 steps to help researchers get started in genome projects by presenting guidelines that are broadly applicable (to any species), sustainable over time, and cover all aspects of genome assembly and annotation projects from start to finish. We review some commonly used approaches, including practical methods to extract high-quality DNA and choices for the best sequencing platforms and library preparations. In addition, we discuss the range of potential bioinformatics pipelines, including structural and functional annotations (e.g., transposable elements and repetitive sequences). This paper also includes information on how to build a wide community for a genome project, the importance of data management, and how to make the data and results Findable, Accessible, Interoperable, and Reusable (FAIR) by submitting them to a public repository and sharing them with the research community.

Comparative Evaluation of Genome Assemblers from Long-Read Sequencing for Plants and Crops.

The availability of recent state-of-the-art long-read sequencing technologies has significantly increased the ease and speed of producing high-quality plant genome assemblies. A wide variety of genome-related software tools are now available and they are typically benchmarked using microbial or model eukaryotic genomes such as Arabidopsis and rice. However, many plant species have much larger and more complex genomes than these, and the choice of tools, parameters, and/or strategies that can be used is not always obvious. Thus, we have compared the metrics of assemblies generated by various pipelines to discuss how assembly quality can be affected by two different assembly strategies. First, we focused on optimizing read preprocessing and assembler variables using eight different de novo assemblers on five different Pacific Biosciences long-read datasets of diploid and tetraploid species. Then, we examined a single scaffolding tool (quickmerge) that has been employed for the postprocessing step. We then merged the outputs from multiple assemblies to produce a higher quality consensus assembly. Then, we benchmarked the assemblies for completeness and accuracy (assembly metrics and BUSCO), computer memory, and CPU times. Two lightweight assemblers, Miniasm/Minimap/Racon and WTDBG, were deemed good for novice users because they involved smaller required learning curves and light computational resources. However, two heavyweight tools, CANU and Flye, should be the first choice when the goal is to achieve accurate and complete assemblies. Our results will provide valuable guidance in future plant genome projects and beyond.

Tools and Strategies for Long-Read Sequencing and De Novo Assembly of Plant Genomes.

The commercial release of third-generation sequencing technologies (TGSTs), giving long and ultra-long sequencing reads, has stimulated the development of new tools for assembling highly contiguous genome sequences with unprecedented accuracy across complex repeat regions. We survey here a wide range of emerging sequencing platforms and analytical tools for de novo assembly, provide background information for each of their steps, and discuss the spectrum of available options. Our decision tree recommends workflows for the generation of a high-quality genome assembly when used in combination with the specific needs and resources of a project.

The Rise and Rise of Nicotiana benthamiana: A Plant for All Reasons.

A decade ago, the value of Nicotiana benthamiana as a tool for plant molecular biologists was beginning to be appreciated. Scientists were using it to study plant-microbe and protein-protein interactions, and it was the species of choice with which to activate plasmid-encoded viruses, screen for gene functions with virus-induced gene silencing (VIGS), and transiently express genes by leaf agroinfiltration. However, little information about the species' origin, diversity, genetics, and genomics was available, and biologists were asking the question of whether N. benthamiana is a second fiddle or virtuoso. In this review, we look at the increased knowledge about the species and its applications over the past decade. Although N. benthamiana may still be the sidekick to Arabidopsis, it shines ever more brightly with realized and yet-to-be-exploited potential.

In-Plant Protection against Helicoverpa armigera by Production of Long hpRNA in Chloroplasts.

Expressing double-stranded RNA (dsRNA) in transgenic plants to silence essential genes within herbivorous pests is referred to as trans-kingdom RNA interference (TK-RNAi) and has emerged as a promising strategy for crop protection. However, the dicing of dsRNA into siRNAs by the plant's intrinsic RNAi machinery may reduce this pesticidal activity. Therefore, genetic constructs, encoding ∼200 nt duplex-stemmed-hairpin (hp) RNAs, targeting the acetylcholinesterase gene of the cotton bollworm, Helicoverpa armigera, were integrated into either the nuclear or the chloroplast genome of Nicotiana benthamiana. Undiced, full-length hpRNAs accumulated in transplastomic lines of N. benthamiana and conferred strong protection against H. armigera herbivory while the hpRNAs of nuclear-transformed plants were processed into siRNAs and gave more modest anti-feeding activity. This suggests that there is little or no RNAi machinery or activity in the chloroplast, that hpRNAs produced within this organelle do not enter the cytoplasm, and that oral delivery of chloroplast-packaged intact hpRNA is a more effective means of delivering TK-RNAi than using nuclear encoded hpRNAs. This contrasts with a recently reported correlation between siRNA expression and effectiveness of TK-RNAi targeting the chitinase gene of H. armigera, but is consistent with reports of efficient TK-RNAi by dsRNA generated in chloroplasts by converging promoters flanking a pest gene sequence and from very small (21 nt-stem) hpRNAs resembling artificial miRNAs. Here we demonstrate that hpRNAs, constructed along the conventional design principles of plant RNAi constructs but integrated into the chloroplast genome, are stable and effective over multiple generations, and hold the promise of providing durable pest resistance in crops.

Optimizing Hybrid de Novo Transcriptome Assembly and Extending Genomic Resources for Giant Freshwater Prawns (Macrobrachium rosenbergii): The Identification of Genes and Markers Associated with Reproduction.

The giant freshwater prawn, Macrobrachium rosenbergii, a sexually dimorphic decapod crustacean is currently the world's most economically important cultured freshwater crustacean species. Despite its economic importance, there is currently a lack of genomic resources available for this species, and this has limited exploration of the molecular mechanisms that control the M. rosenbergii sex-differentiation system more widely in freshwater prawns. Here, we present the first hybrid transcriptome from M. rosenbergii applying RNA-Seq technologies directed at identifying genes that have potential functional roles in reproductive-related traits. A total of 13,733,210 combined raw reads (1720 Mbp) were obtained from Ion-Torrent PGM and 454 FLX. Bioinformatic analyses based on three state-of-the-art assemblers, the CLC Genomic Workbench, Trans-ABySS, and Trinity, that use single and multiple k-mer methods respectively, were used to analyse the data. The influence of multiple k-mers on assembly performance was assessed to gain insight into transcriptome assembly from short reads. After optimisation, de novo assembly resulted in 44,407 contigs with a mean length of 437 bp, and the assembled transcripts were further functionally annotated to detect single nucleotide polymorphisms and simple sequence repeat motifs. Gene expression analysis was also used to compare expression patterns from ovary and testis tissue libraries to identify genes with potential roles in reproduction and sex differentiation. The large transcript set assembled here represents the most comprehensive set of transcriptomic resources ever developed for reproduction traits in M. rosenbergii, and the large number of genetic markers predicted should constitute an invaluable resource for future genetic research studies on M. rosenbergii and can be applied more widely on other freshwater prawn species in the genus Macrobrachium.

Evaluation of potential candidate genes involved in salinity tolerance in striped catfish (Pangasianodon hypophthalmus) using an RNA-Seq approach.

Increasing salinity levels in freshwater and coastal environments caused by sea level rise linked to climate change is now recognized to be a major factor that can impact fish growth negatively, especially for freshwater teleost species. Striped catfish (Pangasianodon hypophthalmus) is an important freshwater teleost that is now widely farmed across the Mekong River Delta in Vietnam. Understanding the basis for tolerance and adaptation to raised environmental salinity conditions can assist the regional culture industry to mitigate predicted impacts of climate change across this region. Attempt of next generation sequencing using the ion proton platform results in more than 174 million raw reads from three tissue libraries (gill, kidney and intestine). Reads were filtered and de novo assembled using a variety of assemblers and then clustered together to generate a combined reference transcriptome. Downstream analysis resulted in a final reference transcriptome that contained 60,585 transcripts with an N50 of 683 bp. This resource was further annotated using a variety of bioinformatics databases, followed by differential gene expression analysis that resulted in 3062 transcripts that were differentially expressed in catfish samples raised under two experimental conditions (0 and 15 ppt). A number of transcripts with a potential role in salinity tolerance were then classified into six different functional gene categories based on their gene ontology assignments. These included; energy metabolism, ion transportation, detoxification, signal transduction, structural organization and detoxification. Finally, we combined the data on functional salinity tolerance genes into a hypothetical schematic model that attempted to describe potential relationships and interactions among target genes to explain the molecular pathways that control adaptive salinity responses in P. hypophthalmus. Our results indicate that P. hypophthalmus exhibit predictable plastic regulatory responses to elevated salinity by means of characteristic gene expression patterns, providing numerous candidate genes for future investigations.

Optimizing de novo transcriptome assembly and extending genomic resources for striped catfish (Pangasianodon hypophthalmus).

Striped catfish (Pangasianodon hypophthalmus) is a commercially important freshwater fish used in inland aquaculture in the Mekong Delta, Vietnam. The culture industry is facing a significant challenge however from saltwater intrusion into many low topographical coastal provinces across the Mekong Delta as a result of predicted climate change impacts. Developing genomic resources for this species can facilitate the production of improved culture lines that can withstand raised salinity conditions, and so we have applied high-throughput Ion Torrent sequencing of transcriptome libraries from six target osmoregulatory organs from striped catfish as a genomic resource for use in future selection strategies. We obtained 12,177,770 reads after trimming and processing with an average length of 97bp. De novo assemblies were generated using CLC Genomic Workbench, Trinity and Velvet/Oases with the best overall contig performance resulting from the CLC assembly. De novo assembly using CLC yielded 66,451 contigs with an average length of 478bp and N50 length of 506bp. A total of 37,969 contigs (57%) possessed significant similarity with proteins in the non-redundant database. Comparative analyses revealed that a significant number of contigs matched sequences reported in other teleost fishes, ranging in similarity from 45.2% with Atlantic cod to 52% with zebrafish. In addition, 28,879 simple sequence repeats (SSRs) and 55,721 single nucleotide polymorphisms (SNPs) were detected in the striped catfish transcriptome. The sequence collection generated in the current study represents the most comprehensive genomic resource for P. hypophthalmus available to date. Our results illustrate the utility of next-generation sequencing as an efficient tool for constructing a large genomic database for marker development in non-model species.

The extremophile Nicotiana benthamiana has traded viral defence for early vigour.

A single lineage of Nicotiana benthamiana is widely used as a model plant(1) and has been instrumental in making revolutionary discoveries about RNA interference (RNAi), viral defence and vaccine production. It is peerless in its susceptibility to viruses and its amenability in transiently expressing transgenes(2,3). These unparalleled characteristics have been associated both positively and negatively with a disruptive insertion in the RNA-dependent RNA polymerase 1 gene, Rdr1(4-6). For a plant so routinely used in research, the origin, diversity and evolution of the species, and the basis of its unusual abilities, have been relatively unexplored. Here, by comparison with wild accessions from across the spectrum of the species' natural distribution, we show that the laboratory strain of N. benthamiana is an extremophile originating from a population that has retained a mutation in Rdr1 for ∼0.8 Myr and thereby traded its defence capacity for early vigour and survival in the extreme habitat of central Australia. Reconstituting Rdr1 activity in this isolate provided protection. Silencing the functional allele in a wild strain rendered it hypersusceptible and was associated with a doubling of seed size and enhanced early growth rate. These findings open the way to a deeper understanding of the delicate balance between protection and vigour.

A transcriptomic analysis of striped catfish (Pangasianodon hypophthalmus) in response to salinity adaptation: De novo assembly, gene annotation and marker discovery.

The striped catfish (Pangasianodon hypophthalmus) culture industry in the Mekong Delta in Vietnam has developed rapidly over the past decade. The culture industry now however, faces some significant challenges, especially related to climate change impacts notably from predicted extensive saltwater intrusion into many low topographical coastal provinces across the Mekong Delta. This problem highlights a need for development of culture stocks that can tolerate more saline culture environments as a response to expansion of saline water-intruded land. While a traditional artificial selection program can potentially address this need, understanding the genomic basis of salinity tolerance can assist development of more productive culture lines. The current study applied a transcriptomic approach using Ion PGM technology to generate expressed sequence tag (EST) resources from the intestine and swim bladder from striped catfish reared at a salinity level of 9ppt which showed best growth performance. Total sequence data generated was 467.8Mbp, consisting of 4,116,424 reads with an average length of 112bp. De novo assembly was employed that generated 51,188 contigs, and allowed identification of 16,116 putative genes based on the GenBank non-redundant database. GO annotation, KEGG pathway mapping, and functional annotation of the EST sequences recovered with a wide diversity of biological functions and processes. In addition, more than 11,600 simple sequence repeats were also detected. This is the first comprehensive analysis of a striped catfish transcriptome, and provides a valuable genomic resource for future selective breeding programs and functional or evolutionary studies of genes that influence salinity tolerance in this important culture species.

A candidate gene association study for growth performance in an improved giant freshwater prawn (Macrobrachium rosenbergii ) culture line.

A candidate gene approach using type I single nucleotide polymorphism (SNP) markers can provide an effective method for detecting genes and gene regions that underlie phenotypic variation in adaptively significant traits. In the absence of available genomic data resources, transcriptomes were recently generated in Macrobrachium rosenbergii to identify candidate genes and markers potentially associated with growth. The characterisation of 47 candidate loci by ABI re-sequencing of four cultured and eight wild samples revealed 342 putative SNPs. Among these, 28 SNPs were selected in 23 growth-related candidate genes to genotype in 200 animals selected for improved growth performance in an experimental GFP culture line in Vietnam. The associations between SNP markers and individual growth performance were then examined. For additive and dominant effects, a total of three exonic SNPs in glycogen phosphorylase (additive), heat shock protein 90 (additive and dominant) and peroxidasin (additive), and a total of six intronic SNPs in ankyrin repeats-like protein (additive and dominant), rolling pebbles (dominant), transforming growth factor-ß induced precursor (dominant), and UTP-glucose-1-phosphate uridylyltransferase 2 (dominant) genes showed significant associations with the estimated breeding values in the experimental animals (P =0.001-0.031). Individually, they explained 2.6-4.8 % of the genetic variance (R²=0.026-0.048). This is the first large set of SNP markers reported for M. rosenbergii and will be useful for confirmation of associations in other samples or culture lines as well as having applications in marker-assisted selection in future breeding programs.

Genome Sequence of Marine Bacterium Idiomarina sp. Strain 28-8, Isolated from Korean Ark Shells.

Idiomarina sp. strain 28-8 is an aerobic, Gram-negative, flagellar bacterium isolated from the bodies of ark shells (Scapharca broughtonii) collected from underwater sediments in Gangjin Bay, South Korea. Here, we present the draft genome sequence of Idiomarina sp. 28-8 (2,971,606 bp, with a G+C content of 46.9%), containing 2,795 putative coding sequences.

Molecular and functional characterizations of a Kunitz-type serine protease inhibitor FcKuSPI of the shrimp Fenneropenaeus chinensis.

Serine proteinase inhibitors play important and diverse roles in biological processes such as coagulation, defense mechanisms, and immune responses. Here, we identified and characterized a Kunitz-type proteinase inhibitor, designated FcKuSPI, of the BPTI/Kunitz family of serine proteinase inhibitors from the hemocyte cDNA library of the shrimp Fenneropenaeus chinensis. The deduced amino acid sequence of FcKuSPI comprises 80 residues with a putative signal peptide of 15 amino acids. The predicted molecular weight of the mature peptide is 7.66 kDa and its predicted isoelectric point is 8.84. FcKuSPI includes a Kunitz domain containing six conserved cysteine residues that are predicted to form three disulfide bonds. FcKuSPI shares 44-53% homology with BPTI/Kunitz family members from other species. FcKuSPI mRNA was expressed highly in the hemocytes and moderately in muscle in healthy shrimp. Recombinant FcKuSPI protein demonstrated anti-protease activity against trypsin and anticoagulant activity against citrated human plasma in a dose-dependent manner in in vitro assays.

Analysis of genome survey sequences and SSR marker development for Siamese Mud Carp, Henicorhynchus siamensis, using 454 pyrosequencing.

Siamese mud carp (Henichorynchus siamensis) is a freshwater teleost of high economic importance in the Mekong River Basin. However, genetic data relevant for delineating wild stocks for management purposes currently are limited for this species. Here, we used 454 pyrosequencing to generate a partial genome survey sequence (GSS) dataset to develop simple sequence repeat (SSR) markers from H. siamensis genomic DNA. Data generated included a total of 65,954 sequence reads with average length of 264 nucleotides, of which 2.79% contain SSR motifs. Based on GSS-BLASTx results, 10.5% of contigs and 8.1% singletons possessed significant similarity (E value < 10(-5)) with the majority matching well to reported fish sequences. KEGG analysis identified several metabolic pathways that provide insights into specific potential roles and functions of sequences involved in molecular processes in H. siamensis. Top protein domains detected included reverse transcriptase and the top putative functional transcript identified was an ORF2-encoded protein. One thousand eight hundred and thirty seven sequences containing SSR motifs were identified, of which 422 qualified for primer design and eight polymorphic loci have been tested with average observed and expected heterozygosity estimated at 0.75 and 0.83, respectively. Regardless of their relative levels of polymorphism and heterozygosity, microsatellite loci developed here are suitable for further population genetic studies in H. siamensis and may also be applicable to other related taxa.

Transcriptomics of a giant freshwater prawn (Macrobrachium rosenbergii): de novo assembly, annotation and marker discovery.

BACKGROUND: Giant freshwater prawn (Macrobrachium rosenbergii or GFP), is the most economically important freshwater crustacean species. However, as little is known about its genome, 454 pyrosequencing of cDNA was undertaken to characterise its transcriptome and identify genes important for growth. METHODOLOGY AND PRINCIPAL FINDINGS: A collection of 787,731 sequence reads (244.37 Mb) obtained from 454 pyrosequencing analysis of cDNA prepared from muscle, ovary and testis tissues taken from 18 adult prawns was assembled into 123,534 expressed sequence tags (ESTs). Of these, 46% of the 8,411 contigs and 19% of 115,123 singletons possessed high similarity to sequences in the GenBank non-redundant database, with most significant (E value < 1e(-5)) contig (80%) and singleton (84%) matches occurring with crustacean and insect sequences. KEGG analysis of the contig open reading frames identified putative members of several biological pathways potentially important for growth. The top InterProScan domains detected included RNA recognition motifs, serine/threonine-protein kinase-like domains, actin-like families, and zinc finger domains. Transcripts derived from genes such as actin, myosin heavy and light chain, tropomyosin and troponin with fundamental roles in muscle development and construction were abundant. Amongst the contigs, 834 single nucleotide polymorphisms, 1198 indels and 658 simple sequence repeats motifs were also identified. CONCLUSIONS: The M. rosenbergii transcriptome data reported here should provide an invaluable resource for improving our understanding of this species' genome structure and biology. The data will also instruct future functional studies to manipulate or select for genes influencing growth that should find practical applications in aquaculture breeding programs.

Development of type I genetic markers from expressed sequence tags in highly polymorphic species.

Expressed sequence tag (EST) databases provide a primary source of nuclear DNA sequences for genetic marker development in non-model organisms. To date, the process has been relatively inefficient for several reasons: 1) priming site polymorphism in the template leads to inferior or erratic amplification; 2) introns in the target amplicon are too large and/or numerous to allow effective amplification under standard screening conditions; and 3) at least occasionally, a PCR primer straddles an exon-intron junction and is unable to bind to genomic DNA template. The first is only a minor issue for species or strains with low heterozygosity but becomes a significant problem for species with high genomic variation, such as marine organisms with extremely large effective population sizes. Problems arising from unanticipated introns are unavoidable but are most pronounced in intron-rich species, such as vertebrates and lophotrochozoans. We present an approach to marker development in the Pacific oyster Crassostrea gigas, a highly polymorphic and intron-rich species, which minimizes these problems, and should be applicable to other non-model species for which EST databases are available. Placement of PCR primers in the 3' end of coding sequence and 3' UTR improved PCR success rate from 51% to 97%. Almost all (37 of 39) markers developed for the Pacific oyster were polymorphic in a small test panel of wild and domesticated oysters.