De Novo Assembly of a New Solanum pennellii Accession Using Nanopore Sequencing

Third generation sequencing transforming plant genome research: Current trends and challenges

In recent years, third-generation sequencing (TGS) technologies have transformed genomics and transcriptomics research, providing novel opportunities for significant discoveries. The long-read sequencing platforms, with their unique advantages over next-generation sequencing (NGS), including a definitive protocol, reduced operational time, and real-time sequencing, possess the potential to transform plant genomics. TGS optimizes and enhances the efficiency of data analysis by removing the necessity for time-consuming assembly tools. The current review examines the development and application of bioinformatics tools for data analysis and annotation, driven by the rapid advancement of TGS platforms like Oxford Nanopore Technologies and Pacific Biosciences. Transcriptome analysis utilizing TGS has been extensively employed to elucidate complex plant transcriptomes and genomes, particularly those characterized by high frequencies of duplicated genomes and repetitive sequences. As a result, current methodologies that allow for generating transcriptomes and comprehensive whole-genome sequences of complex plant genomes employing tailored hybrid sequencing techniques that integrate NGS and TGS technologies have been emphasized herein. This paper, thus, articulates a vision for a future in which TGS effectively addresses the challenges faced in plant research, offering a comprehensive understanding of its advantages, applications, limitations, and promising prospects.

Nanopore sequencing: flourishing in its teenage years

Over the past decade, nanopore sequencing has experienced significant advancements and changes, transitioning from an initially emerging technology to a significant instrument in the field of genomic sequencing. However, as advancements in next-generation sequencing technology persist, nanopore sequencing also improves. This paper reviews the developments, applications, and outlook on nanopore sequencing technology. Currently, nanopore sequencing supports both DNA and RNA sequencing, making it widely applicable in areas such as telomere-to-telomere (T2T) genome assembly, direct RNA sequencing (DRS), and metagenomics. The openness and versatility of nanopore sequencing have established it as a preferred option for an increasing number of research teams, signaling a transformative influence on life science research. As the nanopore sequencing technology advances, it provides a faster, more cost-effective approach with extended read lengths, demonstrating the significant potential for complex genome assembly, pathogen detection, environmental monitoring, and human disease research, offering a fresh perspective in sequencing technologies.

Triumphs of genomic-assisted breeding in crop improvement

Conventional breeding approaches have played a significant role in meeting the food demand remarkably well until now. However, the increasing population, yield plateaus in certain crops, and limited recombination necessitate using genomic resources for genomics-assisted crop improvement programs. As a result of advancements in the next-generation sequence technology, GABs have developed dramatically to characterize allelic variants and facilitate their rapid and efficient incorporation in crop improvement programs. Genomics-assisted breeding (GAB) has played an important role in harnessing the potential of modern genomic tools, exploiting allelic variation from genetic resources and developing cultivars over the past decade. The availability of pangenomes for major crops has been a significant development, albeit with varying degrees of completeness. Even though adopting these technologies is essentially determined on economic grounds and cost-effective assays, which create a wealth of information that can be successfully used to exploit the latent potential of crops. GAB has been instrumental in harnessing the potential of modern genomic resources and exploiting allelic variation for genetic enhancement and cultivar development. GAB strategies will be indispensable for designing future crops and are expected to play a crucial role in breeding climate-smart crop cultivars with higher nutritional value.

Solanum pennellii (LA5240) backcross inbred lines (BILs) for high resolution mapping in tomato

Wild species are an invaluable source of new traits for crop improvement. Over the years, the tomato community bred cultivated lines that carry introgressions from different species of the tomato tribe to facilitate trait discovery and mapping. The next phase in such projects is to find the genes that drive the identified phenotypes. This can be achieved by genotyping a few thousand individuals resulting in fine mapping that can potentially identify the causative gene. To couple trait discovery and fine mapping, we are presenting large, recombination-rich, Backcross Inbred Line (BIL) populations involving an unexplored accession of the wild, green-fruited species Solanum pennellii (LA5240; the 'Lost' Accession) with two modern tomato inbreds: LEA, determinate, and TOP, indeterminate. The LEA and TOP BILs are in BC2F6-8 generation and include 1400 and 500 lines, respectively. The BILs were genotyped with 5000 SPET markers, showing that in the euchromatic regions there was one recombinant every 17-18 Kb while in the heterochromatin a recombinant every 600-700 Kb (TOP and LEA respectively). To gain perspective on the topography of recombination we compared five independent members of the Self-pruning gene family with their respective neighboring genes; based on PCR markers, in all cases we found recombinants. Further mapping analysis of two known morphological mutations that segregated in the BILs (self-pruning and hairless) showed that the maximal delimited intervals were 73 Kb and 210 Kb, respectively, and included the known causative genes. The 'Lost'_BILs provide a solid framework to study traits derived from a drought-tolerant wild tomato.

Deep learning the cis-regulatory code for gene expression in selected model plants

Elucidating the relationship between non-coding regulatory element sequences and gene expression is crucial for understanding gene regulation and genetic variation. We explored this link with the training of interpretable deep learning models predicting gene expression profiles from gene flanking regions of the plant species Arabidopsis thaliana, Solanum lycopersicum, Sorghum bicolor, and Zea mays. With over 80% accuracy, our models enabled predictive feature selection, highlighting e.g. the significant role of UTR regions in determining gene expression levels. The models demonstrated remarkable cross-species performance, effectively identifying both conserved and species-specific regulatory sequence features and their predictive power for gene expression. We illustrated the application of our approach by revealing causal links between genetic variation and gene expression changes across fourteen tomato genomes. Lastly, our models efficiently predicted genotype-specific expression of key functional gene groups, exemplified by underscoring known phenotypic and metabolic differences between Solanum lycopersicum and its wild, drought-resistant relative, Solanum pennellii.

The genome of Haberlea rhodopensis provides insights into the mechanisms for tolerance to multiple extreme environments

Haberlea rhodopensis, a resurrection species, is the only plant known to be able to survive multiple extreme environments, including desiccation, freezing temperatures, and long-term darkness. However, the molecular mechanisms underlying tolerance to these stresses are poorly studied. Here, we present a high-quality genome of Haberlea and found that ~ 23.55% of the 44,306 genes are orphan. Comparative genomics analysis identified 89 significantly expanded gene families, of which 25 were specific to Haberlea. Moreover, we demonstrated that Haberlea preserves its resurrection potential even in prolonged complete darkness. Transcriptome profiling of plants subjected to desiccation, darkness, and low temperatures revealed both common and specific footprints of these stresses, and their combinations. For example, PROTEIN PHOSPHATASE 2C (PP2C) genes were substantially induced in all stress combinations, while PHYTOCHROME INTERACTING FACTOR 1 (PIF1) and GROWTH RESPONSE FACTOR 4 (GRF4) were induced only in darkness. Additionally, 733 genes with unknown functions and three genes encoding transcription factors specific to Haberlea were specifically induced/repressed upon combination of stresses, rendering them attractive targets for future functional studies. The study provides a comprehensive understanding of the genomic architecture and reports details of the mechanisms of multi-stress tolerance of this resurrection species that will aid in developing strategies that allow crops to survive extreme and multiple abiotic stresses.

Quantitative detection of the Ralstonia solanacearum species complex in soil by qPCR combined with a recombinant internal control strain

IMPORTANCE

DNA-based detection and quantification of soil-borne pathogens, such as the Ralstonia solanacearum species complex (RSSC), plays a vital role in risk assessment, but meanwhile, precise quantification is difficult due to the poor purity and yield of the soil DNA retrieved. The internal sample process control (ISPC) strain RsPC we developed solved this problem and significantly improved the accuracy of quantification of RSSC in different soils. ISPC-based quantitative PCR detection is a method especially suitable for the quantitative detection of microbes in complex matrices (such as soil and sludge) containing various PCR inhibitors and for those not easy to lyse (like Gram-positive bacteria, fungi, and thick-wall cells like resting spores). In addition, the use of ISPC strains removes additional workload on the preparation of high-quality template DNA and facilitates the development of high-throughput quantitative detection techniques for soil microbes.

Two long read-based genome assembly and annotation of polyploidy woody plants, Hibiscus syriacus L. using PacBio and Nanopore platforms

Improvements in long read DNA sequencing and related techniques facilitated the generation of complex eukaryotic genomes. Despite these advances, the quality of constructed plant reference genomes remains relatively poor due to the large size of genomes, high content of repetitive sequences, and wide variety of ploidy. Here, we developed the de novo sequencing and assembly of high polyploid plant genome, Hibiscus syriacus, a flowering plant species of the Malvaceae family, using the Oxford Nanopore Technologies and Pacific Biosciences Sequel sequencing platforms. We investigated an efficient combination of high-quality and high-molecular-weight DNA isolation procedure and suitable assembler to achieve optimal results using long read sequencing data. We found that abundant ultra-long reads allow for large and complex polyploid plant genome assemblies with great recovery of repetitive sequences and error correction even at relatively low depth Nanopore sequencing data and polishing compared to previous studies. Collectively, our combination provides cost effective methods to improve genome continuity and quality compared to the previously reported reference genome by accessing highly repetitive regions. The application of this combination may enable genetic research and breeding of polyploid crops, thus leading to improvements in crop production.

Identification of errors in draft genome assemblies at single-nucleotide resolution for quality assessment and improvement

Assembly of a high-quality genome is important for downstream comparative and functional genomic studies. However, most tools for genome assembly assessment only give qualitative reports, which do not pinpoint assembly errors at specific regions. Here, we develop a new reference-free tool, Clipping information for Revealing Assembly Quality (CRAQ), which maps raw reads back to assembled sequences to identify regional and structural assembly errors based on effective clipped alignment information. Error counts are transformed into corresponding assembly evaluation indexes to reflect the assembly quality at single-nucleotide resolution. Notably, CRAQ distinguishes assembly errors from heterozygous sites or structural differences between haplotypes. This tool can clearly indicate low-quality regions and potential structural error breakpoints; thus, it can identify misjoined regions that should be split for further scaffold building and improvement of the assembly. We have benchmarked CRAQ on multiple genomes assembled using different strategies, and demonstrated the misjoin correction for improving the constructed pseudomolecules.

Characterization of the horse chestnut genome reveals the evolution of aescin and aesculin biosynthesis

Horse chestnut (Aesculus chinensis) is an important medicinal tree that contains various bioactive compounds, such as aescin, barrigenol-type triterpenoid saponins (BAT), and aesculin, a glycosylated coumarin. Herein, we report a 470.02 Mb genome assembly and characterize an Aesculus-specific whole-genome duplication event, which leads to the formation and duplication of two triterpenoid biosynthesis-related gene clusters (BGCs). We also show that AcOCS6, AcCYP716A278, AcCYP716A275, and AcCSL1 genes within these two BGCs along with a seed-specific expressed AcBAHD6 are responsible for the formation of aescin. Furthermore, we identify seven Aesculus-originated coumarin glycoside biosynthetic genes and achieve the de novo synthesis of aesculin in E. coli. Collinearity analysis shows that the collinear BGC segments can be traced back to early-diverging angiosperms, and the essential gene-encoding enzymes necessary for BAT biosynthesis are recruited before the splitting of Aesculus, Acer, and Xanthoceras. These findings provide insight on the evolution of gene clusters associated with medicinal tree metabolites.

Regulation of plant epigenetic memory in response to cold and heat stress: towards climate resilient agriculture

Plants have evolved to adapt and grow in hot and cold climatic conditions. Some also adapt to daily and seasonal temperature changes. Epigenetic modifications play an important role in regulating plant tolerance under such conditions. DNA methylation and post-translational modifications of histone proteins influence gene expression during plant developmental stages and under stress conditions, including cold and heat stress. While short-term modifications are common, some modifications may persist and result in stress memory that can be inherited by subsequent generations. Understanding the mechanisms of epigenomes responding to stress and the factors that trigger stress memory is crucial for developing climate-resilient agriculture, but such an integrated view is currently limited. This review focuses on the plant epigenetic stress memory during cold and heat stress. It also discusses the potential of machine learning to modify stress memory through epigenetics to develop climate-resilient crops.

Estimated Nucleotide Reconstruction Quality Symbols of Basecalling Tools for Oxford Nanopore Sequencing

Currently, one of the fastest-growing DNA sequencing technologies is nanopore sequencing. One of the key stages involved in processing sequencer data is the basecalling process, where the input sequence of currents measured on the nanopores of the sequencer reproduces the DNA sequences, called DNA reads. Many of the applications dedicated to basecalling, together with the DNA sequence, provide the estimated quality of the reconstruction of a given nucleotide (quality symbols are contained on every fourth line of the FASTQ file; each nucleotide in the FASTQ file corresponds to exactly one estimated nucleotide reconstruction quality symbol). Herein, we compare the estimated nucleotide reconstruction quality symbols (signs from every fourth line of the FASTQ file) reported by other basecallers. The conducted experiments consisted of basecalling the same raw datasets from the nanopore device by other basecallers and comparing the provided quality symbols, denoting the estimated quality of the nucleotide reconstruction. The results show that the estimated quality reported by different basecallers may vary, depending on the tool used, particularly in terms of range and distribution. Moreover, we mapped basecalled DNA reads to reference genomes and calculated matched and mismatched rates for groups of nucleotides with the same quality symbol. Finally, the presented paper shows that the estimated nucleotide reconstruction quality reported in the basecalling process is not used in any investigated tool for processing nanopore DNA reads.

Comparative Functional Genome Analysis Reveals the Habitat Adaptation and Biocontrol Characteristics of Plant Growth-Promoting Bacteria in NCBI Databases

Plant growth-promoting bacteria (PGPB) are a group of beneficial microorganisms that include 60 bacterial genera, such as Bacillus, Pseudomonas, and Burkholderia, which widely colonize plant leaves and soil, promote plant growth, and/or inhibit pathogen infection. However, the genetic factors underpinning adaptation of PGPB to plant leaves and soil remain poorly understood. In this study, we performed a comparative functional genome analysis approach to investigate the functional genes of 195 leaf-associated (LA) and 283 soil-associated (SA) PGPB strains and their roles in adapting to their environment, using 95 strains from other-associated (OA) environmental habitats with growth-promoting or antimicrobial functions as negative controls. Comparison analysis of the enrichment of nonredundant (NR) protein sequence databases showed that cytochrome P450, DNA repair, and motor chemotaxis genes were significantly enriched in LA PGPB strains related to environmental adaptation, while cell wall-degrading enzymes, TetR transcriptional regulatory factors, and sporulation-related genes were highly enriched in SA PGPB strains. Additionally, analysis of carbohydrate-active enzymes demonstrated that glycosyltransferases (GTs) and glycoside hydrolases (GHs) were abundant families in all PGPB strains, which is in favor of plant growth, and enriched in SA PGPB strains. Except for most Bacillus strains, SA PGPB genomes contained significantly more secondary metabolism clusters than LA PGPB. Most LA PGPB contained hormone biosynthesis genes, which may contribute to plant growth promotion, while SA PGPB harbored numerous carbohydrate and antibiotic metabolism genes. In summary, this study further deepens our understanding of the habitat adaptation and biocontrol characteristics of LA and SA PGPB strains. IMPORTANCE Plant growth-promoting bacteria (PGPB) are essential for the effectiveness of biocontrol agents in plant phyllosphere and rhizosphere. However, little is known about the ecological adaptation of PGPB to different habitats. In this study, comparative functional genome analysis of leaf-associated (LA), soil-associated (SA), and other-associated (OA) PGPB strains was performed. We found that genes related to the metabolism of hormones were enriched in LA PGPB. Carbohydrate and antibiotic metabolism genes were enriched in SA PGPB, which likely facilitated their adaptation to the plant growth environment. Our findings provide genetic insights on LA and SA PGPB strains' ecological adaptation and biocontrol characteristics.

Insights into the Transcriptomics of Crop Wild Relatives to Unravel the Salinity Stress Adaptive Mechanisms

The narrow genomic diversity of modern cultivars is a major bottleneck for enhancing the crop's salinity stress tolerance. The close relatives of modern cultivated plants, crop wild relatives (CWRs), can be a promising and sustainable resource to broaden the diversity of crops. Advances in transcriptomic technologies have revealed the untapped genetic diversity of CWRs that represents a practical gene pool for improving the plant's adaptability to salt stress. Thus, the present study emphasizes the transcriptomics of CWRs for salinity stress tolerance. In this review, the impacts of salt stress on the plant's physiological processes and development are overviewed, and the transcription factors (TFs) regulation of salinity stress tolerance is investigated. In addition to the molecular regulation, a brief discussion on the phytomorphological adaptation of plants under saline environments is provided. The study further highlights the availability and use of transcriptomic resources of CWR and their contribution to pangenome construction. Moreover, the utilization of CWRs' genetic resources in the molecular breeding of crops for salinity stress tolerance is explored. Several studies have shown that cytoplasmic components such as calcium and kinases, and ion transporter genes such as Salt Overly Sensitive 1 (SOS1) and High-affinity Potassium Transporters (HKTs) are involved in the signaling of salt stress, and in mediating the distribution of excess Na⁺ ions within the plant cells. Recent comparative analyses of transcriptomic profiling through RNA sequencing (RNA-Seq) between the crops and their wild relatives have unraveled several TFs, stress-responsive genes, and regulatory proteins for generating salinity stress tolerance. This review specifies that the use of CWRs transcriptomics in combination with modern breeding experimental approaches such as genomic editing, de novo domestication, and speed breeding can accelerate the CWRs utilization in the breeding programs for enhancing the crop's adaptability to saline conditions. The transcriptomic approaches optimize the crop genomes with the accumulation of favorable alleles that will be indispensable for designing salt-resilient crops.

6mA DNA Methylation on Genes in Plants Is Associated with Gene Complexity, Expression and Duplication

N6-methyladenine (6mA) DNA methylation has emerged as an important epigenetic modification in eukaryotes. Nevertheless, the evolution of the 6mA methylation of homologous genes after species and after gene duplications remains unclear in plants. To understand the evolution of 6mA methylation, we detected the genome-wide 6mA methylation patterns of four lotus plants (Nelumbo nucifera) from different geographic origins by nanopore sequencing and compared them to patterns in Arabidopsis and rice. Within lotus, the genomic distributions of 6mA sites are different from the widely studied 5mC methylation sites. Consistently, in lotus, Arabidopsis and rice, 6mA sites are enriched around transcriptional start sites, positively correlated with gene expression levels, and preferentially retained in highly and broadly expressed orthologs with longer gene lengths and more exons. Among different duplicate genes, 6mA methylation is significantly more enriched and conserved in whole-genome duplicates than in local duplicates. Overall, our study reveals the convergent patterns of 6mA methylation evolution based on both lineage and duplicate gene divergence, which underpin their potential role in gene regulatory evolution in plants.

Super-pangenome analyses highlight genomic diversity and structural variation across wild and cultivated tomato species

Effective utilization of wild relatives is key to overcoming challenges in genetic improvement of cultivated tomato, which has a narrow genetic basis; however, current efforts to decipher high-quality genomes for tomato wild species are insufficient. Here, we report chromosome-scale tomato genomes from nine wild species and two cultivated accessions, representative of Solanum section Lycopersicon, the tomato clade. Together with two previously released genomes, we elucidate the phylogeny of Lycopersicon and construct a section-wide gene repertoire. We reveal the landscape of structural variants and provide entry to the genomic diversity among tomato wild relatives, enabling the discovery of a wild tomato gene with the potential to increase yields of modern cultivated tomatoes. Construction of a graph-based genome enables structural-variant-based genome-wide association studies, identifying numerous signals associated with tomato flavor-related traits and fruit metabolites. The tomato super-pangenome resources will expedite biological studies and breeding of this globally important crop.

A high-quality chromosome-level Eutrema salsugineum genome, an extremophile plant model

BACKGROUND

Eutrema salsugineum (2n = 14), a halophyte in the family Brassicaceae, is an attractive model to study abiotic stress tolerance in plants. Two versions of E. salsugineum genomes that previously reported were based on relatively short reads; thus, the repetitive regions were difficult to characterize.

RESULTS

We report the sequencing and assembly of the E. salsugineum (Shandong accession) genome using long-read sequencing and chromosome conformation capture data. We generated Oxford Nanopore long reads at high depth (> 60X) of genome coverage with additional short reads for error correction. The new assembly has a total size of 295.5 Mb with 52.8% repetitive sequences, and the karyotype of E. salsugineum is consistent with the ancestral translocation Proto-Calepineae Karyotype structure in both order and orientation. Compared with previous assemblies, this assembly has higher contiguity, especially in the centromere region. Based on this new assembly, we predicted 25,399 protein-coding genes and identified the positively selected genes associated with salt and drought stress responses.

CONCLUSION

The new genome assembly will provide a valuable resource for future genomic studies and facilitate comparative genomic analysis with other plants.

Revealing evolution of tropane alkaloid biosynthesis by analyzing two genomes in the Solanaceae family

Tropane alkaloids (TAs) are widely distributed in the Solanaceae, while some important medicinal tropane alkaloids (mTAs), such as hyoscyamine and scopolamine, are restricted to certain species/tribes in this family. Little is known about the genomic basis and evolution of TAs biosynthesis and specialization in the Solanaceae. Here, we present chromosome-level genomes of two representative mTAs-producing species: Atropa belladonna and Datura stramonium. Our results reveal that the two species employ a conserved biosynthetic pathway to produce mTAs despite being distantly related within the nightshade family. A conserved gene cluster combined with gene duplication underlies the wide distribution of TAs in this family. We also provide evidence that branching genes leading to mTAs likely have evolved in early ancestral Solanaceae species but have been lost in most of the lineages, with A. belladonna and D. stramonium being exceptions. Furthermore, we identify a cytochrome P450 that modifies hyoscyamine into norhyoscyamine. Our results provide a genomic basis for evolutionary insights into the biosynthesis of TAs in the Solanaceae and will be useful for biotechnological production of mTAs via synthetic biology approaches.

Shoot-root signal circuit: Phytoremediation of heavy metal contaminated soil

High concentrations of heavy metals in the environment will cause serious harm to ecosystems and human health. It is urgent to develop effective methods to control soil heavy metal pollution. Phytoremediation has advantages and potential for soil heavy metal pollution control. However, the current hyperaccumulators have the disadvantages of poor environmental adaptability, single enrichment species and small biomass. Based on the concept of modularity, synthetic biology makes it possible to design a wide range of organisms. In this paper, a comprehensive strategy of "microbial biosensor detection - phytoremediation - heavy metal recovery" for soil heavy metal pollution control was proposed, and the required steps were modified by using synthetic biology methods. This paper summarizes the new experimental methods that promote the discovery of synthetic biological elements and the construction of circuits, and combs the methods of producing transgenic plants to facilitate the transformation of constructed synthetic biological vectors. Finally, the problems that should be paid more attention to in the remediation of soil heavy metal pollution based on synthetic biology were discussed.

Vegetable biology and breeding in the genomics era

Vegetable crops provide a rich source of essential nutrients for humanity and represent critical economic values to global rural societies. However, genetic studies of vegetable crops have lagged behind major food crops, such as rice, wheat and maize, thereby limiting the application of molecular breeding. In the past decades, genome sequencing technologies have been increasingly applied in genetic studies and breeding of vegetables. In this review, we recapitulate recent progress on reference genome construction, population genomics and the exploitation of multi-omics datasets in vegetable crops. These advances have enabled an in-depth understanding of their domestication and evolution, and facilitated the genetic dissection of numerous agronomic traits, which jointly expedites the exploitation of state-of-the-art biotechnologies in vegetable breeding. We further provide perspectives of further directions for vegetable genomics and indicate how the ever-increasing omics data could accelerate genetic, biological studies and breeding in vegetable crops.

The LeEIX Locus Determines Pathogen Resistance in Tomato

The mechanisms underlying the ability of plants to differentiate between pathogens and commensals in their environment are currently unresolved. It has been suggested that spatiotemporal regulation of pattern-recognition receptor (PRR) content could be one of the components providing plants with the ability to distinguish between pathogens and nonpathogenic microbes. The LeEIX PRRs recognize xylanases derived from beneficial or commensal plant colonizers of Trichoderma species, including the xylanase known as EIX. Here, we investigated possible general roles of PRRs from the LeEIX locus in immunity and pathogen resistance in tomato. Mutating the inhibitory PRR LeEIX1, or overexpressing the activating PRR LeEIX2, resulted in resistance to a wide range of pathogens and increased basal and elicited immunity. LeEIX1 knockout caused increases in the expression level of several tested PRRs, including FLS2, as well as bacterial pathogen resistance coupled with an increase in flg22-mediated immunity. The wild tomato relative Solanum pennellii contains inactive LeEIX PRR variants. S. pennellii does not respond to elicitation with the LeEIX PRR ligand EIX. Given that EIX is derived from a mostly nonpathogenic microbe, the connection of its PRRs to disease resistance has not previously been investigated directly. Here, we observed that compared with S. lycopersicum cultivar M82, S. pennellii was more sensitive to several fungal and bacterial pathogens. Our results suggest that the LeEIX locus might determine resistance to fungal necrotrophs, whereas the resistance to biotrophs is effected in combination with a gene/quantitative trait locus not within the LeEIX locus.

Going broad and deep: sequencing-driven insights into plant physiology, evolution, and crop domestication

Deep sequencing is a term that has become embedded in the plant genomic literature in recent years and with good reason. A torrent of (largely) high-quality genomic and transcriptomic data has been collected and most of this has been publicly released. Indeed, almost 1000 plant genomes have been reported (www.plabipd.de) and the 2000 Plant Transcriptomes Project has long been completed. The EarthBioGenome project will dwarf even these milestones. That said, massive progress in understanding plant physiology, evolution, and crop domestication has been made by sequencing broadly (across a species) as well as deeply (within a single individual). We will outline the current state of the art in genome and transcriptome sequencing before we briefly review the most visible of these broad approaches, namely genome-wide association and transcriptome-wide association studies, as well as the compilation of pangenomes. This will include both (i) the most commonly used methods reliant on single nucleotide polymorphisms and short InDels and (ii) more recent examples which consider structural variants. We will subsequently present case studies exemplifying how their application has brought insight into either plant physiology or evolution and crop domestication. Finally, we will provide conclusions and an outlook as to the perspective for the extension of such approaches to different species, tissues, and biological processes.

Advances in sequencing technologies for amyotrophic lateral sclerosis research

Amyotrophic lateral sclerosis (ALS) is caused by upper and lower motor neuron loss and has a fairly rapid disease progression, leading to fatality in an average of 2-5 years after symptom onset. Numerous genes have been implicated in this disease; however, many cases remain unexplained. Several technologies are being used to identify regions of interest and investigate candidate genes. Initial approaches to detect ALS genes include, among others, linkage analysis, Sanger sequencing, and genome-wide association studies. More recently, next-generation sequencing methods, such as whole-exome and whole-genome sequencing, have been introduced. While those methods have been particularly useful in discovering new ALS-linked genes, methodological advances are becoming increasingly important, especially given the complex genetics of ALS. Novel sequencing technologies, like long-read sequencing, are beginning to be used to uncover the contribution of repeat expansions and other types of structural variation, which may help explain missing heritability in ALS. In this review, we discuss how popular and/or upcoming methods are being used to discover ALS genes, highlighting emerging long-read sequencing platforms and their role in aiding our understanding of this challenging disease.

Dynamic Pooling Improves Nanopore Base Calling Accuracy

In nanopore sequencing, electrical signal is measured as DNA molecules pass through the sequencing pores. Translating these signals into DNA bases (base calling) is a highly non-trivial task, and its quality has a large impact on the sequencing accuracy. The most successful nanopore base callers to date use convolutional neural networks (CNN) to accomplish the task. Convolutional layers in CNNs are typically composed of filters with constant window size, performing best in analysis of signals with uniform speed. However, the speed of nanopore sequencing varies greatly both within reads and between sequencing runs. Here, we present dynamic pooling, a novel neural network component, which addresses this problem by adaptively adjusting the pooling ratio. To demonstrate the usefulness of dynamic pooling, we developed two base callers: Heron and Osprey. Heron improves the accuracy beyond the experimental high-accuracy base caller Bonito developed by Oxford Nanopore. Osprey is a fast base caller that can compete in accuracy with Guppy high-accuracy mode, but does not require GPU acceleration and achieves a near real-time speed on common desktop CPUs. Availability: https://github.com/fmfi-compbio/osprey, https://github.com/fmfi-compbio/heron.

The Chromosome-level genome of Aesculus wilsonii provides new insights into terpenoid biosynthesis and Aesculus evolution

Aesculus L. (buckeye and horse chestnut) are woody plant species with important horticultural and medicinal values. Aesculus seeds are widely used as biomedicine and cosmetic ingredients due to their saponins. We report a chromosomal-scale genome of Aesculus wilsonii. Sequences amounting to a total of 579.01 Mb were assembled into 20 chromosomes. More than half of the genome (54.46%) were annotated as repetitive sequences, and 46,914 protein-coding genes were predicted. In addition to the widespread gamma event with core eudicots, a unique whole-genome duplication (WGD) event (17.69 Mya) occurred in Aesculus after buckeye differentiated from longan. Due to WGD events and tandem duplications, the related synthetic genes of triterpene saponins unique to Aesculus increased significantly. Combined with transcriptome characterization, the study preliminarily resolved the biosynthetic pathway of triterpenoid saponins like aescin in A. wilsonii genome. Analyses of the resequencing of 104 buckeye accessions revealed clear relationship between the geographic distribution and genetic differentiation of buckeye trees in China. We found that the buckeye species found in southern Shaanxi is A. wilsonii rather than A. chinensis. Population dynamics analysis further suggests that the population size and evolution of existing buckeye species have been influenced by climate fluctuations during the Pleistocene and recent domestication events. The genome of A. wilsonii and population genomics of Aesculus provide a resource for future research on Hippocastanaceae. These findings will contribute to the utilization and diversity protection of Aesculus.

Genome sequence for the blue-flowered Andean shrub Iochroma cyaneum reveals extensive discordance across the berry clade of Solanaceae

The tomato (Solanum lycopersicum L.) family, Solanaceae, is a model clade for a wide range of applied and basic research questions. Currently, reference-quality genomes are available for over 30 species from seven genera, and these include numerous crops as well as wild species [e.g., Jaltomata sinuosa (Miers) Mione and Nicotiana attenuata Torr. ex S. Watson]. Here we present the genome of the showy-flowered Andean shrub Iochroma cyaneum (Lindl.) M. L. Green, a woody lineage from the tomatillo (Physalis philadelphica Lam.) subfamily Physalideae. The assembled size of the genome (2.7 Gb) is more similar in size to pepper (Capsicum annuum L.) (2.6 Gb) than to other sequenced diploid members of the berry clade of Solanaceae [e.g., potato (Solanum tuberosum L.), tomato, and Jaltomata]. Our assembly recovers 92% of the conserved orthologous set, suggesting a nearly complete genome for this species. Most of the genomic content is repetitive (69%), with Gypsy elements alone accounting for 52% of the genome. Despite the large amount of repetitive content, most of the 12 I. cyaneum chromosomes are highly syntenic with tomato. Bayesian concordance analysis provides strong support for the berry clade, including I. cyaneum, but reveals extensive discordance along the backbone, with placement of chili pepper and Jaltomata being highly variable across gene trees. The I. cyaneum genome contributes to a growing wealth of genomic resources in Solanaceae and underscores the need for expanded sampling of diverse berry genomes to dissect major morphological transitions.

A chromosome-level genome assembly provides insights into the environmental adaptability and outbreaks of Chlorops oryzae

Chlorops oryzae is a pest of rice that has caused severe damage to crops in major rice-growing areas in recent years. We generated a 447.60 Mb high-quality chromosome-level genome with contig and scaffold N50 values of 1.17 Mb and 117.57 Mb, respectively. Hi-C analysis anchored 93.22% scaffolds to 4 chromosomes. The relatively high expression level of Heat Shock Proteins (HSPs) and antioxidant genes in response to thermal stress suggests these genes may play a role in the environmental adaptability of C. oryzae. The identification of multiple pathways that regulate reproductive development (juvenile hormone, 20-hydroxyecdsone, and insulin signaling pathways) provides evidence that these pathways also play an important role in vitellogenesis and thus insect population maintenance. These findings identify possible reasons for the increased frequency of outbreaks of C. oryzae in recent years. Our chromosome-level genome assembly may provide a basis for further genetic studies of C. oryzae, and promote the development of novel, sustainable strategies to control this pest.

A chromosome-level genome assembly for the rice pest, Chlorops oryzae, pinpoints molecular pathways that might contribute toward increased outbreaks for this important crop pest.

The chromosome-level genome provides insight into the molecular mechanism underlying the tortuous-branch phenotype of Prunus mume

Plant with naturally twisted branches is referred to as a tortuous-branch plant, which have extremely high ornamental value due to their zigzag shape and the natural twisting of their branches. Prunus mume is an important woody ornamental plant. However, the molecular mechanism underlying this unique trait in Prunus genus is unknown. Here, we present a chromosome-level genome assembly of the cultivated P. mume var. tortuosa created using Oxford Nanopore combined with Hi-C scaffolding, which resulted in a 237.8 Mb genome assembly being anchored onto eight pseudochromosomes. Molecular dating indicated that P. mume is the most recently differentiated species in Prunus. Genes associated with cell division, development and plant hormones play essential roles in the formation of tortuous branch trait. A putative regulatory pathway for the tortuous branch trait was constructed based on gene expression levels. Furthermore, after transferring candidate PmCYCD genes into Arabidopsis thaliana, we found that seedlings overexpressing these genes exhibited curled rosette leaves. Our results provide insights into the evolutionary history of recently differentiated species in Prunus genus, the molecular basis of stem morphology, and the molecular mechanism underlying the tortuous branch trait and highlight the utility of multi-omics in deciphering the properties of P. mume plant architecture.

A chromosome-level genome assembly of Paracymoriza distinctalis (Lepidoptera: Crambidae: Acentropinae)

Paracymoriza distinctalis is a semiaquatic lepidopteran insect, which is of great value for studying the differentiation of the Pyraloidea superfamily. However, the understanding of heredity, evolution, and functional genomics of P. distinctalis are limited by few genome-wide resources. Here, we applied PacBio sequencing and the chromosome capture technique to assemble the first P. distinctalis genome from a single female individual. The genome size is 1.2 Gb with 32 chromosomes and the N50 is 38.91 Mb. Approximately 576.37 Mb, accounting for 48.93% of the genome, was identified as repeats. The genome comprises 39,003 protein-coding genes, 66.56% of which were functionally annotated. Comparative genomics analysis suggested that the common ancestor of P. distinctalis and Chilo suppressalis lived ~83.5 million years ago. This chromosome-level genome assembly work is not only conducive to the understanding of P. distinctalis, but also may promote the study of the genomes of other lepidopteran species.

Advances and applications of single-cell omics technologies in plant research

Single-cell sequencing approaches reveal the intracellular dynamics of individual cells and answer biological questions with high-dimensional catalogs of millions of cells, including genomics, transcriptomics, chromatin accessibility, epigenomics, and proteomics data across species. These emerging yet thriving technologies have been fully embraced by the field of plant biology, with a constantly expanding portfolio of applications. Here, we introduce the current technical advances used for single-cell omics, especially single-cell genome and transcriptome sequencing. Firstly, we overview methods for protoplast and nucleus isolation and genome and transcriptome amplification. Subsequently, we use well-executed benchmarking studies to highlight advances made through the application of single-cell omics techniques. Looking forward, we offer a glimpse of additional hurdles and future opportunities that will introduce broad adoption of single-cell sequencing with revolutionary perspectives in plant biology.

A Solanum lycopersicoides reference genome facilitates insights into tomato specialized metabolism and immunity

Wild relatives of tomato are a valuable source of natural variation in tomato breeding, as many can be hybridized to the cultivated species (Solanum lycopersicum). Several, including Solanum lycopersicoides, have been crossed to S. lycopersicum for the development of ordered introgression lines (ILs), facilitating breeding for desirable traits. Despite the utility of these wild relatives and their associated ILs, few finished genome sequences have been produced to aid genetic and genomic studies. Here we report a chromosome-scale genome assembly for S. lycopersicoides LA2951, which contains 37 938 predicted protein-coding genes. With the aid of this genome assembly, we have precisely delimited the boundaries of the S. lycopersicoides introgressions in a set of S. lycopersicum cv. VF36 × LA2951 ILs. We demonstrate the usefulness of the LA2951 genome by identifying several quantitative trait loci for phenolics and carotenoids, including underlying candidate genes, and by investigating the genome organization and immunity-associated function of the clustered Pto gene family. In addition, syntenic analysis of R2R3MYB genes sheds light on the identity of the Aubergine locus underlying anthocyanin production. The genome sequence and IL map provide valuable resources for studying fruit nutrient/quality traits, pathogen resistance, and environmental stress tolerance. We present a new genome resource for the wild species S. lycopersicoides, which we use to shed light on the Aubergine locus responsible for anthocyanin production. We also provide IL boundary mappings, which facilitated identifying novel carotenoid quantitative trait loci of which one was likely driven by an uncharacterized lycopene β-cyclase whose function we demonstrate.

Chromosome-level pepino genome provides insights into genome evolution and anthocyanin biosynthesis in Solanaceae

Pepino (Solanum muricatum, 2n = 2x = 24), a member of the Solanaceae family, is an important globally grown fruit. Herein, we report high-quality, chromosome-level pepino genomes. The 91.67% genome sequence is anchored to 12 chromosomes, with a total length of 1.20 Gb and scaffold N50 of 87.03 Mb. More than half the genome comprises repetitive sequences. In addition to the shared ancient whole-genome triplication (WGT) event in eudicots, an additional new WGT event was present in the pepino. Our findings suggest that pepinos experienced chromosome rearrangements, fusions, and gene loss after a WGT event. The large number of gene removals indicated the instability of Solanaceae genomes, providing opportunities for species divergence and natural selection. The paucity of disease-resistance genes (NBS) in pepino and eggplant has been explained by extensive loss and limited generation of genes after WGT events in Solanaceae. The outbreak of NBS genes was not synchronized in Solanaceae species, which occurred before the Solanaceae WGT event in pepino, tomato, and tobacco, whereas it was almost synchronized with WGT events in the other four Solanaceae species. Transcriptome and comparative genomic analyses revealed several key genes involved in anthocyanin biosynthesis. Although an extra WGT event occurred in Solanaceae, CHS genes related to anthocyanin biosynthesis in grapes were still significantly expanded compared with those in Solanaceae species. Proximal and tandem duplications contributed to the expansion of CHS genes. In conclusion, the pepino genome and annotation facilitate further research into important gene functions and comparative genomic analysis in Solanaceae.

A chromosome scale tomato genome built from complementary PacBio and Nanopore sequences alone reveals extensive linkage drag during breeding

The assembly and scaffolding of plant crop genomes facilitate the characterization of genetically diverse cultivated and wild germplasm. The cultivated tomato (Solanum lycopersicum) has been improved through the introgression of genetic material from related wild species, including resistance to pandemic strains of tobacco mosaic virus (TMV) from Solanum peruvianum. Here we applied PacBio HiFi and ONT Nanopore sequencing to develop independent, highly contiguous and complementary assemblies of an inbred TMV-resistant tomato variety. We show specific examples of how HiFi and ONT datasets can complement one another to improve assembly contiguity. We merged the HiFi and ONT assemblies to generate a long-read-only assembly where all 12 chromosomes were represented as 12 contiguous sequences (N50 = 68.5 Mbp). This chromosome scale assembly did not require scaffolding using an orthogonal data type. The merged assembly was validated by chromosome conformation capture data and is highly consistent with previous tomato genome assemblies that made use of genetic maps and Hi-C for scaffolding. Our long-read-only assembly reveals that a complex series of structural variants linked to the TMV resistance gene likely contributed to linkage drag of a 64.1-Mbp region of the S. peruvianum genome during tomato breeding. Through marker studies and ONT-based comprehensive haplotyping we show that this minimal introgression region is present in six cultivated tomato hybrid varieties developed in three commercial breeding programs. Our results suggest that complementary long read technologies can facilitate the rapid generation of near-complete genome sequences.

Designing Uncorrelated Address Constrain for DNA Storage by DMVO Algorithm

At present, huge amounts of data are being produced every second, a situation that will gradually overwhelm current storage technology. DNA is a storage medium that features high storage density and long-term stability and is now considered to be a feasible storage solution. Errors are easily made during the sequencing and synthesis of DNA, however. In order to reduce the error rate, novel uncorrelated address constrain are reported, and a Damping Multi-Verse Optimizer (DMVO)algorithm is proposed to construct a set of DNA coding, which is used as the non-payload. The DMVO algorithm exchanges objects through black/white holes in order to achieve a stable state and adds damping factors as disturbances. Compared with previous work, the coding set obtained by the DMVO algorithm is larger in size and of higher quality. The results of this study reveal that the size of the DNA storage coding set obtained by the DMVO algorithm increased by 4-16 percent, and the variance of the melting temperature decreased by 3-18 percent.

A High-Quality, Chromosome-Level Genome Provides Insights Into Determinate Flowering Time and Color of Cotton Rose (Hibiscus mutabilis)

Hibiscus mutabilis (cotton rose) is a deciduous shrub or small tree of the Malvaceae family. Here, we report a chromosome-scale assembly of the H. mutabilis genome based on a combination of single-molecule sequencing and Hi-C technology. We obtained an optimized assembly of 2.68 Gb with a scaffold N50 length of 54.7 Mb. An integrated strategy of homology-based, de novo, and transcriptome-based gene predictions identified 118,222 protein-coding genes. Repetitive DNA sequences made up 58.55% of the genome, and LTR retrotransposons were the most common repetitive sequence type, accounting for 53.15% of the genome. Through the use of Hi-C data, we constructed a chromosome-scale assembly in which Nanopore scaffolds were assembled into 46 pseudomolecule sequences. We identified important genes involved in anthocyanin biosynthesis and documented copy number variation in floral regulators. Phylogenetic analysis indicated that H. mutabilis was closely related to H. syriacus, from which it diverged approximately 15.3 million years ago. The availability of cotton rose genome data increases our understanding of the species' genetic evolution and will support further biological research and breeding in cotton rose, as well as other Malvaceae species.

Hybrid Assembly Improves Genome Quality and Completeness of Trametes villosa CCMB561 and Reveals a Huge Potential for Lignocellulose Breakdown

Trametes villosa is a wood-decaying fungus with great potential to be used in the bioconversion of agro-industrial residues and to obtain high-value-added products, such as biofuels. Nonetheless, the lack of high-quality genomic data hampers studies investigating genetic mechanisms and metabolic pathways in T. villosa, hindering its application in industry. Herein, applying a hybrid assembly pipeline using short reads (Illumina HiSeq) and long reads (Oxford Nanopore MinION), we obtained a high-quality genome for the T. villosa CCMB561 and investigated its genetic potential for lignocellulose breakdown. The new genome possesses 143 contigs, N50 of 1,009,271 bp, a total length of 46,748,415 bp, 14,540 protein-coding genes, 22 secondary metabolite gene clusters, and 426 genes encoding Carbohydrate-Active enzymes. Our CAZome annotation and comparative genomic analyses of nine Trametes spp. genomes revealed T. villosa CCMB561 as the species with the highest number of genes encoding lignin-modifying enzymes and a wide array of genes encoding proteins for the breakdown of cellulose, hemicellulose, and pectin. These results bring to light the potential of this isolate to be applied in the bioconversion of lignocellulose and will support future studies on the expression, regulation, and evolution of genes, proteins, and metabolic pathways regarding the bioconversion of lignocellulosic residues.

Assembling Plant Genomes with Long-Read Sequencing

Continuous improvements in long-read sequencing allow us to tackle increasingly big and complex genomes. Here we present the principles of long-read genome assembly, taking Solanum pennellii nanopore sequencing as an example.

High quality genome assembly of the amitochondriate eukaryote Monocercomonoides exilis

Monocercomonoides exilis is considered the first known eukaryote to completely lack mitochondria. This conclusion is based primarily on a genomic and transcriptomic study which failed to identify any mitochondrial hallmark proteins. However, the available genome assembly has limited contiguity and around 1.5 % of the genome sequence is represented by unknown bases. To improve the contiguity, we re-sequenced the genome and transcriptome of M. exilis using Oxford Nanopore Technology (ONT). The resulting draft genome is assembled in 101 contigs with an N50 value of 1.38 Mbp, almost 20 times higher than the previously published assembly. Using a newly generated ONT transcriptome, we further improve the gene prediction and add high quality untranslated region (UTR) annotations, in which we identify two putative polyadenylation signals present in the 3′UTR regions and characterise the Kozak sequence in the 5′UTR regions. All these improvements are reflected by higher BUSCO genome completeness values. Regardless of an overall more complete genome assembly without missing bases and a better gene prediction, we still failed to identify any mitochondrial hallmark genes, thus further supporting the hypothesis on the absence of mitochondrion.

Genomic prediction using low-coverage portable Nanopore sequencing

Most traits in livestock, crops and humans are polygenic, that is, a large number of loci contribute to genetic variation. Effects at these loci lie along a continuum ranging from common low-effect to rare high-effect variants that cumulatively contribute to the overall phenotype. Statistical methods to calculate the effect of these loci have been developed and can be used to predict phenotypes in new individuals. In agriculture, these methods are used to select superior individuals using genomic breeding values; in humans these methods are used to quantitatively measure an individual’s disease risk, termed polygenic risk scores. Both fields typically use SNP array genotypes for the analysis. Recently, genotyping-by-sequencing has become popular, due to lower cost and greater genome coverage (including structural variants). Oxford Nanopore Technologies’ (ONT) portable sequencers have the potential to combine the benefits genotyping-by-sequencing with portability and decreased turn-around time. This introduces the potential for in-house clinical genetic disease risk screening in humans or calculating genomic breeding values on-farm in agriculture. Here we demonstrate the potential of the later by calculating genomic breeding values for four traits in cattle using low-coverage ONT sequence data and comparing these breeding values to breeding values calculated from SNP arrays. At sequencing coverages between 2X and 4X the correlation between ONT breeding values and SNP array-based breeding values was > 0.92 when imputation was used and > 0.88 when no imputation was used. With an average sequencing coverage of 0.5x the correlation between the two methods was between 0.85 and 0.92 using imputation, depending on the trait. This suggests that ONT sequencing has potential for in clinic or on-farm genomic prediction, however, further work to validate these findings in a larger population still remains.

Molecular Methods for Pathogenic Bacteria Detection and Recent Advances in Wastewater Analysis

With increasing concerns about public health and the development of molecular techniques, new detection tools and the combination of existing approaches have increased the abilities of pathogenic bacteria monitoring by exploring new biomarkers, increasing the sensitivity and accuracy of detection, quantification, and analyzing various genes such as functional genes and antimicrobial resistance genes (ARG). Molecular methods are gradually emerging as the most popular detection approach for pathogens, in addition to the conventional culture-based plate enumeration methods. The analysis of pathogens in wastewater and the back-estimation of infections in the community, also known as wastewater-based epidemiology (WBE), is an emerging methodology and has a great potential to supplement current surveillance systems for the monitoring of infectious diseases and the early warning of outbreaks. However, as a complex matrix, wastewater largely challenges the analytical performance of molecular methods. This review synthesized the literature of typical pathogenic bacteria in wastewater, types of biomarkers, molecular methods for bacterial analysis, and their recent advances in wastewater analysis. The advantages and limitation of these molecular methods were evaluated, and their prospects in WBE were discussed to provide insight for future development.

MinION sequencing from sea ice cryoconites leads to de novo genome reconstruction from metagenomes

Genome reconstruction from metagenomes enables detailed study of individual community members, their metabolisms, and their survival strategies. Obtaining high quality metagenome-assembled genomes (MAGs) is particularly valuable in extreme environments like sea ice cryoconites, where the native consortia are recalcitrant to culture and strong astrobiology analogues. We evaluated three separate approaches for MAG generation from Allen Bay, Nunavut sea ice cryoconites-HiSeq-only, MinION-only, and hybrid (HiSeq + MinION)-where field MinION sequencing yielded a reliable metagenome. The hybrid assembly produced longer contigs, more coding sequences, and more total MAGs, revealing a microbial community dominated by Bacteroidetes. The hybrid MAGs also had the highest completeness, lowest contamination, and highest N50. A putatively novel species of Octadecabacter is among the hybrid MAGs produced, containing the genus's only known instances of genomic potential for nitrate reduction, denitrification, sulfate reduction, and fermentation. This study shows that the inclusion of MinION reads in traditional short read datasets leads to higher quality metagenomes and MAGs for more accurate descriptions of novel microorganisms in this extreme, transient habitat and has produced the first hybrid MAGs from an extreme environment.

Reap the crop wild relatives for breeding future crops

Crop wild relatives (CWRs) have provided breeders with several 'game-changing' traits or genes that have boosted crop resilience and global agricultural production. Advances in breeding and genomics have accelerated the identification of valuable CWRs for use in crop improvement. The enhanced genetic diversity of breeding pools carrying optimum combinations of favorable alleles for targeted crop-growing regions is crucial to sustain genetic gain. In parallel, growing sequence information on wild genomes in combination with precise gene-editing tools provide a fast-track route to transform CWRs into ideal future crops. Data-informed germplasm collection and management strategies together with adequate policy support will be equally important to improve access to CWRs and their sustainable use to meet food and nutrition security targets.

Draft genome assemblies using sequencing reads from Oxford Nanopore Technology and Illumina platforms for four species of North American Fundulus killifish

BACKGROUND

Whole-genome sequencing data from wild-caught individuals of closely related North American killifish species (Fundulus xenicus, Fundulus catenatus, Fundulus nottii, and Fundulus olivaceus) were obtained using long-read Oxford Nanopore Technology (ONT) PromethION and short-read Illumina platforms.

FINDINGS

Draft de novo reference genome assemblies were generated using a combination of long and short sequencing reads. For each species, the PromethION platform was used to generate 30-45× sequence coverage, and the Illumina platform was used to generate 50-160× sequence coverage. Illumina-only assemblies were fragmented with high numbers of contigs, while ONT-only assemblies were error prone with low BUSCO scores. The highest N50 values, ranging from 0.4 to 2.7 Mb, were from assemblies generated using a combination of short- and long-read data. BUSCO scores were consistently >90% complete using the Eukaryota database.

CONCLUSIONS

High-quality genomes can be obtained from a combination of using short-read Illumina data to polish assemblies generated with long-read ONT data. Draft assemblies and raw sequencing data are available for public use. We encourage use and reuse of these data for assembly benchmarking and other analyses.

Tomato leaves under stress: a comparison of stress response to mild abiotic stress between a cultivated and a wild tomato species

Tomato is one of the most produced crop plants on earth and growing in the fields and greenhouses all over the world. Breeding with known traits of wild species can enhance stress tolerance of cultivated crops. In this study, we investigated responses of the transcriptome as well as primary and secondary metabolites in leaves of a cultivated and a wild tomato to several abiotic stresses such as nitrogen deficiency, chilling or warmer temperatures, elevated light intensities and combinations thereof. The wild species responded different to varied temperature conditions compared to the cultivated tomato. Nitrogen deficiency caused the strongest responses and induced in particular the secondary metabolism in both species but to much higher extent in the cultivated tomato. Our study supports the potential of a targeted induction of valuable secondary metabolites in green residues of horticultural production, that will otherwise only be composted after fruit harvest. In particular, the cultivated tomato showed a strong induction in the group of mono caffeoylquinic acids in response to nitrogen deficiency. In addition, the observed differences in stress responses between cultivated and wild tomato can lead to new breeding targets for better stress tolerance.

Telomere-to-telomere gapless chromosomes of banana using nanopore sequencing

Long-read technologies hold the promise to obtain more complete genome assemblies and to make them easier. Coupled with long-range technologies, they can reveal the architecture of complex regions, like centromeres or rDNA clusters. These technologies also make it possible to know the complete organization of chromosomes, which remained complicated before even when using genetic maps. However, generating a gapless and telomere-to-telomere assembly is still not trivial, and requires a combination of several technologies and the choice of suitable software. Here, we report a chromosome-scale assembly of a banana genome (Musa acuminata) generated using Oxford Nanopore long-reads. We generated a genome coverage of 177X from a single PromethION flowcell with near 17X with reads longer than 75 kbp. From the 11 chromosomes, 5 were entirely reconstructed in a single contig from telomere to telomere, revealing for the first time the content of complex regions like centromeres or clusters of paralogous genes.

De Novo Genome Assembly Of Solanum Sitiens Reveals Structural Variation Associated With Drought And Salinity Tolerance

MOTIVATION

Solanum sitiens is a self-incompatible wild relative of tomato, characterised by salt and drought resistance traits, with the potential to contribute through breeding programmes to crop improvement in cultivated tomato. This species has a distinct morphology, classification and ecotype compared to other stress resistant wild tomato relatives such as S. pennellii and S. chilense. Therefore, the availability of a reference genome for S. sitiens will facilitate the genetic and molecular understanding of salt and drought resistance.

RESULTS

A high-quality de novo genome and transcriptome assembly for S. sitiens (Accession LA1974) has been developed. A hybrid assembly strategy was followed using Illumina short reads (∼159X coverage) and PacBio long reads (∼44X coverage), generating a total of ∼262 Gbp of DNA sequence. A reference genome of 1,245 Mbp, arranged in 1,483 scaffolds with a N50 of 1.826 Mbp was generated. Genome completeness was estimated at 95% using the Benchmarking Universal Single-Copy Orthologs (BUSCO) and the K-mer Analysis Tool (KAT). In addition, ∼63 Gbp of RNA-Seq were generated to support the prediction of 31,164 genes from the assembly, and to perform a de novo transcriptome. Lastly, we identified three large inversions compared to S. lycopersicum, containing several drought resistance related genes, such as beta-amylase 1 and YUCCA7.

AVAILABILITY

S. sitiens (LA1974) raw sequencing, transcriptome and genome assembly have been deposited at the NCBI's Sequence Read Archive, under the BioProject number "PRJNA633104".All the commands and scripts necessary to generate the assembly are available at the following github repository: https://github.com/MCorentin/Solanum_sitiens_assembly.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genomic basis of high-altitude adaptation in Tibetan Prunus fruit trees

The Great Himalayan Mountains and their foothills are believed to be the place of origin and development of many plant species. The genetic basis of adaptation to high plateaus is a fascinating topic that is poorly understood at the population level. We comprehensively collected and sequenced 377 accessions of Prunus germplasm along altitude gradients ranging from 2,067 to 4,492 m in the Himalayas. We de novo assembled three high-quality genomes of Tibetan Prunus species. A comparative analysis of Prunus genomes indicated a remarkable expansion of the SINE retrotransposons occurred in the genomes of Tibetan species. We observed genetic differentiation between Tibetan peaches from high and low altitudes and that genes associated with light stress signaling, especially UV stress signaling, were enriched in the differentiated regions. By profiling the metabolomes of Tibetan peach fruit, we determined 379 metabolites had significant genetic correlations with altitudes and that in particular phenylpropanoids were positively correlated with altitudes. We identified 62 Tibetan peach-specific SINEs that colocalized with metabolites differentially accumualted in Tibetan relative to cultivated peach. We demonstrated that two SINEs were inserted in a locus controlling the accumulation of 3-O-feruloyl quinic acid. SINE1 was specific to Tibetan peach. SINE2 was predominant in high altitudes and associated with the accumulation of 3-O-feruloyl quinic acid. These genomic and metabolic data for Prunus populations native to the Himalayan region indicate that the expansion of SINE retrotransposons helped Tibetan Prunus species adapt to the harsh environment of the Himalayan plateau by promoting the accumulation of beneficial metabolites.

A chromosome-level reference genome of red swamp crayfish Procambarus clarkii provides insights into the gene families regarding growth or development in crustaceans

Red swamp crayfish Procambarus clarkii is an ecologically and economically important crustacean species. Here, based on a de novo assembly strategy combining PacBio with Hi-C sequencing, we presented a high quality chromosome-level P. clarkii genome. The assembled genome is 2.75 Gb in size with a contig N50 of 216.75 kb. Transposable elements (TEs) make up the largest fraction of the genome (~79.61%), and LINEs comprise the majority of the TEs. Frequent molting and rapid growth of the red swamp crayfish may be explained by the expansion of multiple gene families regarding growth or development. Phylogenetic analysis revealed that P. clarkii diverged from Portunus trituberculatus at 278-407 million years ago (Mya). PSMC analysis identified multiple bottleneck events of the P. clarkii population between 2 kaBP to 14 kaBP. The obtained P. clarkii genome should not only facilitate us understanding the development and evolution of the crayfish species, but also contribute to the genetic improvement in future breeding selections.