Whole genome profiling physical map and ancestral annotation of tobacco Hicks Broadleaf

Identification and editing of a hybrid lethality gene expands the range of interspecific hybridization potential in Nicotiana

KEY MESSAGE

Identification and inactivation of hybrid lethality genes can be used to expand the available gene pool for improvement of a cultivated crop species. Hybrid lethality is one genetic mechanism that contributes to reproductive isolation in plants and serves as a barrier to use of diverse germplasm for improvement of cultivated species. A classic example is the seedling lethality exhibited by progeny from the Nicotiana tabacum × N. africana interspecific cross. In order to increase the body of knowledge on mechanisms of hybrid lethality in plants, and to potentially develop tools to circumvent them, we utilized a transposon tagging strategy to identify a candidate gene involved in the control of this reaction. N. tabacum gene Nt6549g30 was identified to code for a class of coiled-coil nucleotide-binding site-leucine-rich repeat (CC-NBS-LRR) proteins, the largest class of plant defense proteins. Gene editing, along with other experiments, was used to verify that Nt6549g30 is the gene at the N. tabacum Hybrid Lethality 1 (NtHL1) locus controlling the hybrid lethality reaction in crosses with N. africana. Gene editing of Nt6549g30 was also used to reverse interspecific seedling lethality in crosses between N. tabacum and eight of nine additional tested species from section Suaveolentes. Results further implicate the role of disease resistance-like genes in the evolution of plant species and demonstrate the possibility of expanding the gene pool for a crop species through gene editing.

Characterization of Phn15.1, a Newly Identified Phytophthora nicotianae Resistance QTL in Nicotiana tabacum

Phytophthora nicotianae is an oomycete that causes black shank, one of the most economically important diseases affecting tobacco production worldwide. Identification and introgression of novel genetic variability affecting partial genetic resistance to this pathogen is important because of the increased durability of partial resistance over time as compared with genes conferring immunity. A previous mapping study identified a quantitative trait locus (QTL), hereafter designated as Phn15.1, with a major effect on P. nicotianae resistance in tobacco. In this research, we describe significantly improved resistance of nearly isogenic lines (NILs) of flue-cured tobacco carrying the introgressed Phn15.1 region derived from highly resistant cigar tobacco cultivar Beinhart 1000. The Phn15.1 region appeared to act in an additive or partially dominant manner to positively affect resistance. To more finely resolve the position of the gene or genes underlying the Phn15.1 effect, the QTL was mapped with an increased number of molecular markers (single-nucleotide polymorphisms) identified to reside within the region. Development and evaluation of subNILs containing varying amounts of Beinhart 1000-derived Phn15.1-associated genetic material permitted the localization of the QTL to a genetic interval of approximately 2.7 centimorgans. Importantly, we were able to disassociate the Beinhart 1000 Phn15.1 resistance alleles from a functional NtCPS2 allele(s) which contributes to the accumulation of a diterpene leaf surface exudate considered undesirable for flue-cured and burley tobacco. Information from this research should be of value for marker-assisted introgression of Beinhart 1000-derived partial black shank resistance into flue-cured and burley tobacco breeding programs.

Evolutionary and functional analyses of the 2-oxoglutarate-dependent dioxygenase genes involved in the flavonoid biosynthesis pathway in tobacco

MAIN CONCLUSION

This study illustrates the differences in the gene structure of 2-oxoglutarate-dependent oxygenase involved in flavonoid biosynthesis (2ODD-IFB), and their potential roles in regulating tobacco flavonoid biosynthesis and plant growth. Flavonol synthase (FLS), anthocyanidin synthase (ANS), and flavanone 3β-hydroxylase belong to the 2-oxoglutarate-dependent (2ODD) oxygenase family, and each performs crucial functions in the biosynthesis of flavonoids. We identified two NtFLS genes, two NtANS genes, and four NtF3H genes from Nicotiana tabacum genome, as well as their homologous genes in the N. sylvestris and N. tomentosiformis genomes. Our phylogenetic analysis indicated that these three types of genes split from each other before the divergence of gymnosperms and angiosperms. FLS evolved faster in the eudicot plants, whereas ANS evolved faster in the monocot plants. Gene structure analysis revealed two fragment insertions occurred at different times in the intron one position of tobacco FLS genes. Homologous protein modeling revealed distinct structures in the N terminus of the tobacco 2ODD oxygenases. We found that the expression patterns of genes encoding tobacco 2ODD oxygenases in flavonoids biosynthesis (2ODD-IFB) did not determine the accumulation patterns of flavonoids among various tobacco tissues, but strongly affected the concentration of flavonoids in the tissues, where they were biosynthesized. More carbon resource flowed to the flavonol biosynthesis when NtANS gene was silenced, otherwise more anthocyanidin accumulated when NtFLS gene was repressed. This study illustrates the 2ODD-IFB gene structure evolution, differences among their protein structures, and provides a foundation for regulating plant development and altering flavonoid content and/or composition through the manipulation of plant 2ODD-IFB genes.

Development of a Sequence-Based Reference Physical Map of Pea (Pisum sativum L.)

Whole genome profiling (WGP) is a sequence-based physical mapping technology and uses sequence tags generated by next generation sequencing for construction of bacterial artificial chromosome (BAC) contigs of complex genomes. The physical map provides a framework for assembly of genome sequence and information for localization of genes that are difficult to find through positional cloning. To address the challenges of accurate assembly of the pea genome (∼4.2 GB of which approximately 85% is repetitive sequences), we have adopted the WGP technology for assembly of a pea BAC library. Multi-dimensional pooling of 295,680 BAC clones and sequencing the ends of restriction fragments of pooled DNA generated 1,814 million high quality reads, of which 825 million were deconvolutable to 1.11 million unique WGP sequence tags. These WGP tags were used to assemble 220,013 BACs into contigs. Assembly of the BAC clones using the modified Fingerprinted Contigs (FPC) program has resulted in 13,040 contigs, consisting of 213,719 BACs, and 6,294 singleton BACs. The average contig size is 0.33 Mbp and the N₅₀ contig size is 0.62 Mbp. WGP^TM technology has proved to provide a robust physical map of the pea genome, which would have been difficult to assemble using traditional restriction digestion based methods. This sequence-based physical map will be useful to assemble the genome sequence of pea. Additionally, the 1.1 million WGP tags will support efficient assignment of sequence scaffolds to the BAC clones, and thus an efficient sequencing of BAC pools with targeted genome regions of interest.

Comparative genome-wide characterization leading to simple sequence repeat marker development for Nicotiana

BACKGROUND

Simple sequence repeats (SSRs) are tandem repeats of DNA that have been used to develop robust genetic markers. These molecular markers are powerful tools for basic and applied studies such as molecular breeding. In the model plants in Nicotiana genus e.g. N. benthamiana, a comprehensive assessment of SSR content has become possible now because several Nicotiana genomes have been sequenced. We conducted a genome-wide SSR characterization and marker development across seven Nicotiana genomes.

RESULTS

Here, we initially characterized 2,483,032 SSRs (repeat units of 1-10 bp) from seven genomic sequences of Nicotiana and developed SSR markers using the GMATA® software package. Of investigated repeat units, mono-, di- and tri-nucleotide SSRs account for 98% of all SSRs in Nicotiana. More complex SSR motifs, although rare, are highly variable between Nicotiana genomes. A total of 1,224,048 non-redundant Nicotiana (NIX) markers were developed, of which 99.98% are novel. An efficient and uniform genotyping protocol for NIX markers was developed and validated. We created a web-based database of NIX marker information including amplicon sizes of alleles in each genome for downloading and online analysis.

CONCLUSIONS

The present work constitutes the first deep characterization of SSRs in seven genomes of Nicotiana, and the development of NIX markers for these SSRs. Our online marker database and an efficient genotyping protocol facilitate the application of these markers. The NIX markers greatly expand Nicotiana marker resources, thus providing a useful tool for future research and breeding. We demonstrate a novel protocol for SSR marker development and utilization at the whole genome scale that can be applied to any lineage of organisms. The Tobacco Markers & Primers Database (TMPD) is available at http://biodb.sdau.edu.cn/tmpd/index.html.

pBACode: a random-barcode-based high-throughput approach for BAC paired-end sequencing and physical clone mapping

Applications that use Bacterial Artificial Chromosome (BAC) libraries often require paired-end sequences and knowledge of the physical location of each clone in plates. To facilitate obtaining this information in high-throughput, we generated pBACode vectors: a pool of BAC cloning vectors, each with a pair of random barcodes flanking its cloning site. In a pBACode BAC library, the BAC ends and their linked barcodes can be sequenced in bulk. Barcode pairs are determined by sequencing the empty pBACode vectors, which allows BAC ends to be paired according to their barcodes. For physical clone mapping, the barcodes are used as unique markers for their linked genomic sequence. After multi-dimensional pooling of BAC clones, the barcodes are sequenced and deconvoluted to locate each clone. We generated a pBACode library of 94,464 clones for the flounder Paralichthys olivaceus and obtained paired-end sequence from 95.4% of the clones. Incorporating BAC paired-ends into the genome preassembly improved its continuity by over 10-fold. Furthermore, we were able to use the barcodes to map the physical locations of each clone in just 50 pools, with up to 11 808 clones per pool. Our physical clone mapping located 90.2% of BAC clones, enabling targeted characterization of chromosomal rearrangements.

A reference genome for Nicotiana tabacum enables map-based cloning of homeologous loci implicated in nitrogen utilization efficiency

BACKGROUND

Tobacco (Nicotiana tabacum) is an important plant model system that has played a key role in the early development of molecular plant biology. The tobacco genome is large and its characterisation challenging because it is an allotetraploid, likely arising from hybridisation between diploid N. sylvestris and N. tomentosiformis ancestors. A draft assembly was recently published for N. tabacum, but because of the aforementioned genome complexities it was of limited utility due to a high level of fragmentation.

RESULTS

Here we report an improved tobacco genome assembly, which, aided by the application of optical mapping, achieves an N50 size of 2.17 Mb and enables anchoring of 64% of the genome to pseudomolecules; a significant increase from the previous value of 19%. We use this assembly to identify two homeologous genes that explain the differentiation of the burley tobacco market class, with potential for greater understanding of Nitrogen Utilization Efficiency and Nitrogen Use Efficiency in plants; an important trait for future sustainability of agricultural production.

CONCLUSIONS

Development of an improved genome assembly for N. tabacum enables what we believe to be the first successful map-based gene discovery for the species, and demonstrates the value of an improved assembly for future research in this model and commercially-important species.

MaGuS: a tool for quality assessment and scaffolding of genome assemblies with Whole Genome Profiling™ Data

Background

Scaffolding is an essential step in the genome assembly process. Current methods based on large fragment paired-end reads or long reads allow an increase in contiguity but often lack consistency in repetitive regions, resulting in fragmented assemblies. Here, we describe a novel tool to link assemblies to a genome map to aid complex genome reconstruction by detecting assembly errors and allowing scaffold ordering and anchoring.

Results

We present MaGuS (map-guided scaffolding), a modular tool that uses a draft genome assembly, a Whole Genome Profiling™ (WGP) map, and high-throughput paired-end sequencing data to estimate the quality and to enhance the contiguity of an assembly. We generated several assemblies of the Arabidopsis genome using different scaffolding programs and applied MaGuS to select the best assembly using quality metrics. Then, we used MaGuS to perform map-guided scaffolding to increase contiguity by creating new scaffold links in low-covered and highly repetitive regions where other commonly used scaffolding methods lack consistency.

Conclusions

MaGuS is a powerful reference-free evaluator of assembly quality and a WGP map-guided scaffolder that is freely available at https://github.com/institut-de-genomique/MaGuS. Its use can be extended to other high-throughput sequencing data (e.g., long-read data) and also to other map data (e.g., genetic maps) to improve the quality and the contiguity of large and complex genome assemblies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-0969-x) contains supplementary material, which is available to authorized users.

NtGNL1a ARF-GEF acts in endocytosis in tobacco cells

BACKGROUND

Processes of anterograde and retrograde membrane trafficking play an important role in cellular homeostasis and dynamic rearrangements of the plasma membrane (PM) in all eukaryotes. These processes depend on the activity of adenosine ribosylation factors (ARFs), a family of GTP-binding proteins and their guanine exchange factors (GEFs). However, knowledge on the function and specificity of individual ARF-GEFs for individual steps of membrane trafficking pathways is still limited in plants.

RESULTS

In this work, treatments with various trafficking inhibitors showed that the endocytosis of FM 4-64 is largely dynamin-dependent and relies on proteins containing endocytic tyrosine-based internalization motif and intact cytoskeleton. Interestingly, brefeldin A (BFA), reported previously as an inhibitor of anterograde membrane trafficking in plants, appeared to be the most potent inhibitor of endocytosis in tobacco. In concert with this finding, we demonstrate that the point mutation in the Sec7 domain of the GNOM-LIKE protein1a (NtGNL1a) confers intracellular trafficking pathway-specific BFA resistance. The internalization of FM 4-64 and trafficking of PIN-FORMED1 (PIN1) auxin efflux carrier in BY-2 tobacco cells were studied to reveal the function of the ARF-GEF NtGNL1a in these.

CONCLUSIONS

Altogether, our observations uncovered the role of NtGNL1a in endocytosis, including endocytosis of PM proteins (as PIN1 auxin efflux carrier). Moreover these data emphasize the need of careful evaluation of mode of action of non-native inhibitors in various species. In addition, they demonstrate the potential of tobacco BY-2 cells for selective mapping of ARF-GEF-regulated endomembrane trafficking pathways.

Engineering a Platform for Photosynthetic Pigment, Hormone and Cembrane-Related Diterpenoid Production in Nicotiana tabacum

Plants synthesize a large number of isoprenoids that are of nutritional, medicinal and industrial importance. 1-Deoxy-d-xylulose 5-phosphate reductoisomerase (DXR) catalyzes the first committed step for plastidial isoprenoid biosynthesis. Here, we identified two DXR isogenes, designated NtDXR1 and NtDXR2, from tetraploid common tobacco (Nicotiana tabacum L.). Southern blotting and genotyping analysis revealed that two NtDXR genes existed in the tetraploid tobacco genome; NtDXR1 and NtDXR2 were separately derived from N. tomentosiformis and N. sylvestris. Both NtDXRs were localized in chloroplasts. Expression patterns indicated that NtDXR1 and NtDXR2 had similar expression profiles. NtDXR genes were highly expressed in leaves with or without trichomes; expression was relatively reduced in flowers and stems, weak in leaf trichomes and marginal in roots and seeds. Overexpressing NtDXR1 under control of the 35S promoter resulted in longer primary roots and enhancement of various photosynthetic pigments and hormones in leaves. In contrast, there were no significant changes in cembrane-related diterpenoids synthesized in glandular trichomes. To elucidate further the function of DXR in the biosynthesis of diterpenoids, overexpression vectors for NtDXR1 under the control of a trichome-specific CYP promoter were transferred to tobacco plants. CYP:NtDXR1 tobacco exhibited larger glandular cells and increased cembrane-related diterpenoids in leaf glandular trichomes. Moreover, transcripts of eight MEP (2-C-methyl-d-erythritol 4-phosphate) pathway genes were significantly up-regulated in NtDXR1-overexpressing tobacco plants, indicating that overexpression of NtDXR could boost the expression of downstream genes in the MEP pathway. Our results suggested that overexpression of NtDXR1 could increase the levels of photosynthetic pigments, leaf surface exudates and hormones though the MEP pathway.

A comprehensive, genome-wide analysis of autophagy-related genes identified in tobacco suggests a central role of autophagy in plant response to various environmental cues

Autophagy is an evolutionarily conserved mechanism in both animals and plants, which has been shown to be involved in various essential developmental processes in plants. Nicotiana tabacum is considered to be an ideal model plant and has been widely used for the study of the roles of autophagy in the processes of plant development and in the response to various stresses. However, only a few autophagy-related genes (ATGs) have been identified in tobacco up to now. Here, we identified 30 ATGs belonging to 16 different groups in tobacco through a genome-wide survey. Comprehensive expression profile analysis reveals an abroad expression pattern of these ATGs, which could be detected in all tissues tested under normal growth conditions. Our series tests further reveal that majority of ATGs are sensitive and responsive to different stresses including nutrient starvation, plant hormones, heavy metal and other abiotic stresses, suggesting a central role of autophagy, likely as an effector, in plant response to various environmental cues. This work offers a detailed survey of all ATGs in tobacco and also suggests manifold functions of autophagy in both normal plant growth and plant response to environmental stresses.

Comprehensive analysis of cystatin family genes suggests their putative functions in sexual reproduction, embryogenesis, and seed formation

Cystatins are tightly bound and reversible inhibitors of cysteine proteases in C1A and C13 peptidase families, which have been identified in several species and shown to function in vegetative development and response to biotic/abiotic stresses in plants. Recent work revealed their critical role in regulating programmed cell death during embryogenesis in tobacco and suggested their more comprehensive roles in the process of sexual plant reproduction, although little is known about cystatin family genes in the processes. Here, 10 cystatin family genes in Nicotiana tabacum were identified using an expressed sequence tag (EST)-based gene clone strategy. Analysis of their biochemical properties showed that nine of them have the potency to inhibit the activities of both commercial cathepsin L-like proteases and extracted cysteine proteases from seeds, but with different K i values depending on the types of proteases and the developmental stages of the seed tested. This suggests that cystatin-dependent cathepsin L-like proteolytic pathways are probably important for early seed development. Comprehensive expression profile analysis revealed that cystatin family genes showed manifold variations in their transcription levels in different plant cell types, including the sperm, egg, and zygote, especially in the embryo and seed at different developmental stages. More interestingly, intracellular localization analysis of each cystatin revealed that most members of cystatin families are recognized as secretory proteins with signal peptides that direct them to the endoplasmic reticulum. These results suggest their widespread roles in cell fate determination and cell-cell communication in the process of sexual reproduction, especially in gamete and embryo development, as well as in seed formation.

Whole-genome profiling and shotgun sequencing delivers an anchored, gene-decorated, physical map assembly of bread wheat chromosome 6A

Bread wheat (Triticum aestivum L.) is the most important staple food crop for 35% of the world's population. International efforts are underway to facilitate an increase in wheat production, of which the International Wheat Genome Sequencing Consortium (IWGSC) plays an important role. As part of this effort, we have developed a sequence-based physical map of wheat chromosome 6A using whole-genome profiling (WGP™). The bacterial artificial chromosome (BAC) contig assembly tools fingerprinted contig (fpc) and linear topological contig (ltc) were used and their contig assemblies were compared. A detailed investigation of the contigs structure revealed that ltc created a highly robust assembly compared with those formed by fpc. The ltc assemblies contained 1217 contigs for the short arm and 1113 contigs for the long arm, with an L50 of 1 Mb. To facilitate in silico anchoring, WGP™ tags underlying BAC contigs were extended by wheat and wheat progenitor genome sequence information. Sequence data were used for in silico anchoring against genetic markers with known sequences, of which almost 79% of the physical map could be anchored. Moreover, the assigned sequence information led to the 'decoration' of the respective physical map with 3359 anchored genes. Thus, this robust and genetically anchored physical map will serve as a framework for the sequencing of wheat chromosome 6A, and is of immediate use for map-based isolation of agronomically important genes/quantitative trait loci located on this chromosome.

The tobacco genome sequence and its comparison with those of tomato and potato

The allotetraploid plant Nicotiana tabacum (common tobacco) is a major crop species and a model organism, for which only very fragmented genomic sequences are currently available. Here we report high-quality draft genomes for three main tobacco varieties. These genomes show both the low divergence of tobacco from its ancestors and microsynteny with other Solanaceae species. We identify over 90,000 gene models and determine the ancestral origin of tobacco mosaic virus and potyvirus disease resistance in tobacco. We anticipate that the draft genomes will strengthen the use of N. tabacum as a versatile model organism for functional genomics and biotechnology applications.

The widespread crucifer species Cardamine flexuosa is an allotetraploid with a conserved subgenomic structure

The origin of Cardamine flexuosa (Wavy Bittercress) has been a conundrum for more than six decades. Here we identify its parental species, analyse its genome structure in comparison to parental genomes and describe intergenomic structural variations in C. flexuosa. Genomic in situ hybridization (GISH) and comparative chromosome painting (CCP) uncovered the parental genomes and the chromosome composition of C. flexuosa and its presumed diploid progenitors. Cardamine flexuosa is an allotetraploid (2n = 4x = 32), originating from two diploid species, Cardamine amara and Cardamine hirsuta (2n = 2x = 16). The two parental species display almost perfectly conserved chromosomal collinearity for seven out of the eight chromosomes. A 13 Mb pericentric inversion distinguishes chromosome CA1 from CH1. A comparative cytomolecular map was established for C. flexuosa by CCP/GISH. Whereas conserved chromosome collinearity between the C. amara and C. hirsuta subgenomes might have promoted intergenomic rearrangements through homeologous recombination, only one reciprocal translocation between two homeologues has occurred since the origin of C. flexuosa. The genome of C. flexuosa demonstrates that allopolyploids can maintain remarkably stable subgenomes over 10(4) -10(5)  yr throughout a wide distribution range. By contrast, the rRNA genes underwent genome-specific elimination towards a diploid-like number of loci.

Reference genomes and transcriptomes of Nicotiana sylvestris and Nicotiana tomentosiformis

BACKGROUND

Nicotiana sylvestris and Nicotiana tomentosiformis are members of the Solanaceae family that includes tomato, potato, eggplant and pepper. These two Nicotiana species originate from South America and exhibit different alkaloid and diterpenoid production. N. sylvestris is cultivated largely as an ornamental plant and it has been used as a diploid model system for studies of terpenoid production, plastid engineering, and resistance to biotic and abiotic stress. N. sylvestris and N. tomentosiformis are considered to be modern descendants of the maternal and paternal donors that formed Nicotiana tabacum about 200,000 years ago through interspecific hybridization. Here we report the first genome-wide analysis of these two Nicotiana species.

RESULTS

Draft genomes of N. sylvestris and N. tomentosiformis were assembled to 82.9% and 71.6% of their expected size respectively, with N50 sizes of about 80 kb. The repeat content was 72-75%, with a higher proportion of retrotransposons and copia-like long terminal repeats in N. tomentosiformis. The transcriptome assemblies showed that 44,000-53,000 transcripts were expressed in the roots, leaves or flowers. The key genes involved in terpenoid metabolism, alkaloid metabolism and heavy metal transport showed differential expression in the leaves, roots and flowers of N. sylvestris and N. tomentosiformis.

CONCLUSIONS

The reference genomes of N. sylvestris and N. tomentosiformis represent a significant contribution to the SOL100 initiative because, as members of the Nicotiana genus of Solanaceae, they strengthen the value of the already existing resources by providing additional comparative information, thereby helping to improve our understanding of plant metabolism and evolution.