Palaeohexaploid ancestry for Caryophyllales inferred from extensive gene-based physical and genetic mapping of the sugar beet genome (Beta vulgaris)

The link between ancient whole-genome duplications and cold adaptations in the Caryophyllaceae

PREMISE

The Caryophyllaceae (the carnation family) have undergone multiple transitions into colder climates and convergence on cushion plant adaptation, indicating that they may provide a natural system for cold adaptation research. Previous research has suggested that putative ancient whole-genome duplications (WGDs) are correlated with niche shifts into colder climates across the Caryophyllales. Here, we explored the genomic changes potentially involved in one of these discovered shifts in the Caryophyllaceae.

METHODS

We constructed a data set combining 26 newly generated transcriptomes with 45 published transcriptomes, including 11 cushion plant species across seven genera. With this data set, we inferred a dated phylogeny for the Caryophyllaceae and mapped ancient WGDs and gene duplications onto the phylogeny. We also examined functional groups enriched for gene duplications related to the climatic shift.

RESULTS

The ASTRAL topology was mostly congruent with the current consensus of relationships within the family. We inferred 15 putative ancient WGDs in the family, including eight that have not been previously published. The oldest ancient WGD (ca. 64.4-56.7 million years ago), WGD1, was found to be associated with a shift into colder climates by previous research. Gene regions associated with ubiquitination were overrepresented in gene duplications retained after WGD1 and those convergently retained by cushion plants in Colobanthus and Eremogone, along with other functional annotations.

CONCLUSIONS

Gene family expansions induced by ancient WGDs may have contributed to the shifts to cold climatic niches in the Caryophyllaceae. Transcriptomic data are crucial resources that help unravel heterogeneity in deep-time evolutionary patterns in plants.

Genome-wide identification, phylogenetic classification of histone acetyltransferase genes, and their expression analysis in sugar beet (Beta vulgaris L.) under salt stress

MAIN CONCLUSION

This study identified seven histone acetyltransferase-encoding genes (HATs) from Beta vulgaris L. (sugar beet) genome through bioinformatics tools and analyzed their expression profiles under salt stress. Sugar beet HATs are phylogenetically divided into four families: GNAT, MYST, CBP, and TAFII250. The BvHAT genes were differentially transcribed in leaves, stems, and roots of B. vulgaris salt-resistant (Casino) and -sensitive (Bravo) cultivars under salt stress. Histone acetylation is regulated by histone acetyltransferases (HATs), which catalyze ɛ-amino bond formation between lysine residues and acetyl groups with a cofactor, acetyl-CoA. Even though the HATs are known to participate in stress response and development in model plants, little is known about the functions of HATs in crops. In sugar beet (Beta vulgaris L.), they have not yet been identified and characterized. Here, an in silico analysis of the HAT gene family in sugar beet was performed, and their expression patterns in leaves, stems, and roots of B. vulgaris were analyzed under salt stress. Salt-resistant (Casino) and -sensitive (Bravo) beet cultivars were used for gene expression assays. Seven HATs were identified from sugar beet genome, and named BvHAG1, BvHAG2, BvHAG3, BvHAG4, BvHAC1, BvHAC2, and BvHAF1. The HAT proteins were divided into 4 groups including MYST, GNAT (GCN5, HAT1, ELP3), CBP and TAFII250. Analysis of cis-acting elements indicated that the BvHAT genes might be involved in hormonal regulation, light response, plant development, and abiotic stress response. The BvHAT genes were differentially expressed in leaves, stems, and roots under control and 300 mM NaCl. In roots of B. vulgaris cv. Bravo, the BvHAG1, BvHAG2, BvHAG4, BvHAF1, and BvHAC1 genes were dramatically expressed after 7 and 14 days of salt stress. Interestingly, the BvHAC2 gene was not expressed under both control and stress conditions. However, the expression of BvHAG2, BvHAG3, BvHAG4, BvHAC1, BvHAC2 genes showed a significant increase in response to salt stress in the roots of cv. Casino. This study provides new insights into the potential roles of histone acetyltransferases in sugar beet.

Association analysis of agronomic traits and construction of genetic networks by resequencing of 306 sugar beet (Beta vulgaris L.) lines

Due to the relatively brief domestication history of sugar beet (Beta vulgaris ssp. vulgaris), our understanding of the genomic diversity and functional genes in its cultivars is limited, resulting in slow breeding progress. To address this issue, a total of 306 germplasm materials of major cultivars and breeding lines from China, the USA, and Europe were selected for genome resequencing. We investigated population structure and genetic diversity and performed selective scanning of genomic regions, identifying six novel genes associated with important agronomic traits: the candidate genes DFAX2 and P5CS for skin roughness; the candidate genes FRO5, GL24, and PPR91 for root yield and sugar yield, and the pleiotropic candidate gene POLX for flourishing growth vigour, plant height, crown size, flesh coarseness, and sugar yield. In addition, we constructed a protein-protein interaction network map and a phenotype-gene network map, which provide valuable information for identifying and characterizing functional genes affecting agronomic traits in sugar beet. Overall, our study sheds light on the future improvement of sugar beet agronomic traits at the molecular level.

Genomic and transcriptomic-based analysis of agronomic traits in sugar beet (Beta vulgaris L.) pure line IMA1

Sugar beet (Beta vulgaris L.) is an important sugar-producing and energy crop worldwide. The sugar beet pure line IMA1 independently bred by Chinese scientists is a standard diploid parent material that is widely used in hybrid-breeding programs. In this study, a high-quality, chromosome-level genome assembly for IMA1was conducted, and 99.1% of genome sequences were assigned to nine chromosomes. A total of 35,003 protein-coding genes were annotated, with 91.56% functionally annotated by public databases. Compared with previously released sugar beet assemblies, the new genome was larger with at least 1.6 times larger N50 size, thereby substantially improving the completeness and continuity of the sugar beet genome. A Genome-Wide Association Studies analysis identified 10 disease-resistance genes associated with three important beet diseases and five genes associated with sugar yield per hectare, which could be key targets to improve sugar productivity. Nine highly expressed genes associated with pollen fertility of sugar beet were also identified. The results of this study provide valuable information to identify and dissect functional genes affecting sugar beet agronomic traits, which can increase sugar beet production and help screen for excellent sugar beet breeding materials. In addition, information is provided that can precisely incorporate biotechnology tools into breeding efforts.

A contiguous de novo genome assembly of sugar beet EL10 (Beta vulgaris L.)

A contiguous assembly of the inbred 'EL10' sugar beet (Beta vulgaris ssp. vulgaris) genome was constructed using PacBio long-read sequencing, BioNano optical mapping, Hi-C scaffolding, and Illumina short-read error correction. The EL10.1 assembly was 540 Mb, of which 96.2% was contained in nine chromosome-sized pseudomolecules with lengths from 52 to 65 Mb, and 31 contigs with a median size of 282 kb that remained unassembled. Gene annotation incorporating RNA-seq data and curated sequences via the MAKER annotation pipeline generated 24,255 gene models. Results indicated that the EL10.1 genome assembly is a contiguous genome assembly highly congruent with the published sugar beet reference genome. Gross duplicate gene analyses of EL10.1 revealed little large-scale intra-genome duplication. Reduced gene copy number for well-annotated gene families relative to other core eudicots was observed, especially for transcription factors. Variation in genome size in B. vulgaris was investigated by flow cytometry among 50 individuals producing estimates from 633 to 875 Mb/1C. Read-depth mapping with short-read whole-genome sequences from other sugar beet germplasm suggested that relatively few regions of the sugar beet genome appeared associated with high-copy number variation.

Genomic distances reveal relationships of wild and cultivated beets

Cultivated beets (Beta vulgaris ssp. vulgaris), including sugar beet, rank among the most important crops. The wild ancestor of beet crops is the sea beet Beta vulgaris ssp. maritima. Species and subspecies of wild beets are readily crossable with cultivated beets and are thus available for crop improvement. To study genomic relationships in the genus Beta, we sequence and analyse 606 beet genomes, encompassing sugar beet, sea beet, B. v. adanensis, B. macrocarpa, and B. patula. We observe two genetically distinct groups of sea beets, one from the Atlantic coast and the other from the Mediterranean area. Genomic comparisons based on k-mers identify sea beets from Greece as the closest wild relatives of sugar beet, suggesting that domestication of the ancestors of sugar beet may be traced to this area. Our work provides comprehensive insight into the phylogeny of wild and cultivated beets and establishes a framework for classification of further accessions of unknown (sub-)species assignment.

An anchored chromosome-scale genome assembly of spinach improves annotation and reveals extensive gene rearrangements in euasterids

Spinach (Spinacia oleracea L.) is a member of the Caryophyllales family, a basal eudicot asterid that consists of sugar beet (Beta vulgaris L. subsp. vulgaris), quinoa (Chenopodium quinoa Willd.), and amaranth (Amaranthus hypochondriacus L.). With the introduction of baby leaf types, spinach has become a staple food in many homes. Production issues focus on yield, nitrogen-use efficiency and resistance to downy mildew (Peronospora effusa). Although genomes are available for the above species, a chromosome-level assembly exists only for quinoa, allowing for proper annotation and structural analyses to enhance crop improvement. We independently assembled and annotated genomes of the cultivar Viroflay using short-read strategy (Illumina) and long-read strategies (Pacific Biosciences) to develop a chromosome-level, genetically anchored assembly for spinach. Scaffold N50 for the Illumina assembly was 389 kb, whereas that for Pacific BioSciences was 4.43 Mb, representing 911 Mb (93% of the genome) in 221 scaffolds, 80% of which are anchored and oriented on a sequence-based genetic map, also described within this work. The two assemblies were 99.5% collinear. Independent annotation of the two assemblies with the same comprehensive transcriptome dataset show that the quality of the assembly directly affects the annotation with significantly more genes predicted (26,862 vs. 34,877) in the long-read assembly. Analysis of resistance genes confirms a bias in resistant gene motifs more typical of monocots. Evolutionary analysis indicates that Spinacia is a paleohexaploid with a whole-genome triplication followed by extensive gene rearrangements identified in this work. Diversity analysis of 75 lines indicate that variation in genes is ample for hypothesis-driven, genomic-assisted breeding enabled by this work.

The effect of heat stress on sugar beet recombination

Meiotic recombination plays a crucial role in plant breeding through the creation of new allelic combinations. Therefore, lack of recombination in some genomic regions constitutes a constraint for breeding programmes. In sugar beet, one of the major crops in Europe, recombination occurs mainly in the distal portions of the chromosomes, and so the development of simple approaches to change this pattern is of considerable interest for future breeding and genetics. In the present study, the effect of heat stress on recombination in sugar beet was studied by treating F₁ plants at 28 °C/25 °C (day/night) and genotyping the progeny. F₁ plants were reciprocally backcrossed allowing the study of male and female meiosis separately. Genotypic data indicated an overall increase in crossover frequency of approximately one extra crossover per meiosis, with an associated increase in pericentromeric recombination under heat treatment. Our data indicate that the changes were mainly induced by alterations in female meiosis only, showing that heterochiasmy in sugar beet is reduced under heat stress. Overall, despite the associated decrease in fertility, these data support the potential use of heat stress to foster recombination in sugar beet breeding programmes.

Genome-wide sequence identification and expression analysis of ARF family in sugar beet (Beta vulgaris L.) under salinity stresses

Auxin response factor (ARF) proteins respond to biological and abiotic stresses and play important roles in regulating plant growth and development. In this study, based on the genome-wide database of sugar beet, 16 BvARF proteins were identified. A detailed investigation into the BvARF family is performed, including analysis of the conserved domains, chromosomal locations, phylogeny, exon-intron structure, conserved motifs, subcellular localization, gene ontology (GO) annotations and expression profiles of BvARF under salt-tolerant condition. The majority of BvARF proteins contain B3 domain, AUX_RESP domain and AUX/IAA domain and a few lacked of AUX/IAA domain. Phylogenetic analysis suggests that the 16 BvARF proteins are clustered into six groups. Expression profile analysis shows that most of these BvARF genes in sugar beet under salinity stress were up-regulated or down-regulated to varying degrees and nine of the BvARF genes changed significantly. They were thought to have a significant response to salinity stress. The current study provides basic information for the BvARF genes and will pave the way for further studies on the roles of BvARF genes in regulating sugar beet's growth, development and responses to salinity stress.

A Partially Phase-Separated Genome Sequence Assembly of the Vitis Rootstock 'Börner' (Vitis riparia × Vitis cinerea) and Its Exploitation for Marker Development and Targeted Mapping

Grapevine breeding has become highly relevant due to upcoming challenges like climate change, a decrease in the number of available fungicides, increasing public concern about plant protection, and the demand for a sustainable production. Downy mildew caused by Plasmopara viticola is one of the most devastating diseases worldwide of cultivated Vitis vinifera. In modern breeding programs, therefore, genetic marker technologies and genomic data are used to develop new cultivars with defined and stacked resistance loci. Potential sources of resistance are wild species of American or Asian origin. The interspecific hybrid of Vitis riparia Gm 183 x Vitis cinerea Arnold, available as the rootstock cultivar 'Börner,' carries several relevant resistance loci. We applied next-generation sequencing to enable the reliable identification of simple sequence repeats (SSR), and we also generated a draft genome sequence assembly of 'Börner' to access genome-wide sequence variations in a comprehensive and highly reliable way. These data were used to cover the 'Börner' genome with genetic marker positions. A subset of these marker positions was used for targeted mapping of the P. viticola resistance locus, Rpv14, to validate the marker position list. Based on the reference genome sequence PN40024, the position of this resistance locus can be narrowed down to less than 0.5 Mbp on chromosome 5.

Genomes of the wild beets Beta patula and Beta vulgaris ssp. maritima

We present draft genome assemblies of Beta patula, a critically endangered wild beet endemic to the Madeira archipelago, and of the closely related Beta vulgaris ssp. maritima (sea beet). Evidence-based reference gene sets for B. patula and sea beet were generated, consisting of 25 127 and 27 662 genes, respectively. The genomes and gene sets of the two wild beets were compared with their cultivated sister taxon B. vulgaris ssp. vulgaris (sugar beet). Large syntenic regions were identified, and a display tool for automatic genome-wide synteny image generation was developed. Phylogenetic analysis based on 9861 genes showing 1:1:1 orthology supported the close relationship of B. patula to sea beet and sugar beet. A comparative analysis of the Rz2 locus, responsible for rhizomania resistance, suggested that the sequenced B. patula accession was rhizomania susceptible. Reference karyotypes for the two wild beets were established, and genomic rearrangements were detected. We consider our data as highly valuable and comprehensive resources for wild beet studies, B. patula conservation management, and sugar beet breeding research.

Meta-analysis of chromosome-scale crossover rate variation in eukaryotes and its significance to evolutionary genomics

Understanding the distribution of crossovers along chromosomes is crucial to evolutionary genomics because the crossover rate determines how strongly a genome region is influenced by natural selection on linked sites. Nevertheless, generalities in the chromosome-scale distribution of crossovers have not been investigated formally. We fill this gap by synthesizing joint information on genetic and physical maps across 62 animal, plant and fungal species. Our quantitative analysis reveals a strong and taxonomically widespread reduction of the crossover rate in the centre of chromosomes relative to their peripheries. We demonstrate that this pattern is poorly explained by the position of the centromere, but find that the magnitude of the relative reduction in the crossover rate in chromosome centres increases with chromosome length. That is, long chromosomes often display a dramatically low crossover rate in their centre, whereas short chromosomes exhibit a relatively homogeneous crossover rate. This observation is compatible with a model in which crossover is initiated from the chromosome tips, an idea with preliminary support from mechanistic investigations of meiotic recombination. Consequently, we show that organisms achieve a higher genome-wide crossover rate by evolving smaller chromosomes. Summarizing theory and providing empirical examples, we finally highlight that taxonomically widespread and systematic heterogeneity in crossover rate along chromosomes generates predictable broad-scale trends in genetic diversity and population differentiation by modifying the impact of natural selection among regions within a genome. We conclude by emphasizing that chromosome-scale heterogeneity in crossover rate should urgently be incorporated into analytical tools in evolutionary genomics, and in the interpretation of resulting patterns.

Disparity, diversity, and duplications in the Caryophyllales

The role played by whole genome duplication (WGD) in plant evolution is actively debated. WGDs have been associated with advantages such as superior colonization, various adaptations, and increased effective population size. However, the lack of a comprehensive mapping of WGDs within a major plant clade has led to uncertainty regarding the potential association of WGDs and higher diversification rates. Using seven chloroplast and nuclear ribosomal genes, we constructed a phylogeny of 5036 species of Caryophyllales, representing nearly half of the extant species. We phylogenetically mapped putative WGDs as identified from analyses on transcriptomic and genomic data and analyzed these in conjunction with shifts in climatic occupancy and lineage diversification rate. Thirteen putative WGDs and 27 diversification shifts could be mapped onto the phylogeny. Of these, four WGDs were concurrent with diversification shifts, with other diversification shifts occurring at more recent nodes than WGDs. Five WGDs were associated with shifts to colder climatic occupancy. While we find that many diversification shifts occur after WGDs, it is difficult to consider diversification and duplication to be tightly correlated. Our findings suggest that duplications may often occur along with shifts in either diversification rate, climatic occupancy, or rate of evolution.

Improved transcriptome sampling pinpoints 26 ancient and more recent polyploidy events in Caryophyllales, including two allopolyploidy events

Studies of the macroevolutionary legacy of polyploidy are limited by an incomplete sampling of these events across the tree of life. To better locate and understand these events, we need comprehensive taxonomic sampling as well as homology inference methods that accurately reconstruct the frequency and location of gene duplications. We assembled a data set of transcriptomes and genomes from 168 species in Caryophyllales, of which 43 transcriptomes were newly generated for this study, representing one of the most densely sampled genomic-scale data sets available. We carried out phylogenomic analyses using a modified phylome strategy to reconstruct the species tree. We mapped the phylogenetic distribution of polyploidy events by both tree-based and distance-based methods, and explicitly tested scenarios for allopolyploidy. We identified 26 ancient and more recent polyploidy events distributed throughout Caryophyllales. Two of these events were inferred to be allopolyploidy. Through dense phylogenomic sampling, we show the propensity of polyploidy throughout the evolutionary history of Caryophyllales. We also provide a framework for utilizing transcriptome data to detect allopolyploidy, which is important as it may have different macroevolutionary implications compared with autopolyploidy.

A high-quality genome assembly of quinoa provides insights into the molecular basis of salt bladder-based salinity tolerance and the exceptional nutritional value

Chenopodium quinoa is a halophytic pseudocereal crop that is being cultivated in an ever-growing number of countries. Because quinoa is highly resistant to multiple abiotic stresses and its seed has a better nutritional value than any other major cereals, it is regarded as a future crop to ensure global food security. We generated a high-quality genome draft using an inbred line of the quinoa cultivar Real. The quinoa genome experienced one recent genome duplication about 4.3 million years ago, likely reflecting the genome fusion of two Chenopodium parents, in addition to the γ paleohexaploidization reported for most eudicots. The genome is highly repetitive (64.5% repeat content) and contains 54 438 protein-coding genes and 192 microRNA genes, with more than 99.3% having orthologous genes from glycophylic species. Stress tolerance in quinoa is associated with the expansion of genes involved in ion and nutrient transport, ABA homeostasis and signaling, and enhanced basal-level ABA responses. Epidermal salt bladder cells exhibit similar characteristics as trichomes, with a significantly higher expression of genes related to energy import and ABA biosynthesis compared with the leaf lamina. The quinoa genome sequence provides insights into its exceptional nutritional value and the evolution of halophytes, enabling the identification of genes involved in salinity tolerance, and providing the basis for molecular breeding in quinoa.

Single-molecule sequencing and Hi-C-based proximity-guided assembly of amaranth (Amaranthus hypochondriacus) chromosomes provide insights into genome evolution

Background

Amaranth (Amaranthus hypochondriacus) was a food staple among the ancient civilizations of Central and South America that has recently received increased attention due to the high nutritional value of the seeds, with the potential to help alleviate malnutrition and food security concerns, particularly in arid and semiarid regions of the developing world. Here, we present a reference-quality assembly of the amaranth genome which will assist the agronomic development of the species.

Results

Utilizing single-molecule, real-time sequencing (Pacific Biosciences) and chromatin interaction mapping (Hi-C) to close assembly gaps and scaffold contigs, respectively, we improved our previously reported Illumina-based assembly to produce a chromosome-scale assembly with a scaffold N50 of 24.4 Mb. The 16 largest scaffolds contain 98% of the assembly and likely represent the haploid chromosomes (n = 16). To demonstrate the accuracy and utility of this approach, we produced physical and genetic maps and identified candidate genes for the betalain pigmentation pathway. The chromosome-scale assembly facilitated a genome-wide syntenic comparison of amaranth with other Amaranthaceae species, revealing chromosome loss and fusion events in amaranth that explain the reduction from the ancestral haploid chromosome number (n = 18) for a tetraploid member of the Amaranthaceae.

Conclusions

The assembly method reported here minimizes cost by relying primarily on short-read technology and is one of the first reported uses of in vivo Hi-C for assembly of a plant genome. Our analyses implicate chromosome loss and fusion as major evolutionary events in the 2n = 32 amaranths and clearly establish the homoeologous relationship among most of the subgenome chromosomes, which will facilitate future investigations of intragenomic changes that occurred post polyploidization.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0412-4) contains supplementary material, which is available to authorized users.

Functional characterisation and cell specificity of BvSUT1, the transporter that loads sucrose into the phloem of sugar beet (Beta vulgaris L.) source leaves

Sugar beet (Beta vulgaris L.) is one of the most important sugar-producing plants worldwide and provides about one third of the sugar consumed by humans. Here we report on molecular characterisation of the BvSUT1 gene and on the functional characterisation of the encoded transporter. In contrast to the recently identified tonoplast-localised sucrose transporter BvTST2.1 from sugar beet taproots, which evolved within the monosaccharide transporter (MST) superfamily, BvSUT1 represents a classical sucrose transporter and is a typical member of the disaccharide transporter (DST) superfamily. Transgenic Arabidopsis plants expressing the β-GLUCURONIDASE (GUS) reporter gene under control of the BvSUT1-promoter showed GUS histochemical staining of their phloem; an anti-BvSUT1-antiserum identified the BvSUT1 transporter specifically in phloem companion cells. After expression of BvSUT1 cDNA in bakers' yeasts (Saccharomyces cerevisiae) uptake characteristics of the BvSUT1 protein were studied. Moreover, the sugar beet transporter was characterised as a proton-coupled sucrose symporter in Xenopus laevis oocytes. Our findings indicate that BvSUT1 is the sucrose transporter that is responsible for loading of sucrose into the phloem of sugar beet source leaves delivering sucrose to the storage tissue in sugar beet taproot sinks.

Comparison of Spinach Sex Chromosomes with Sugar Beet Autosomes Reveals Extensive Synteny and Low Recombination at the Male-Determining Locus

Spinach (Spinacia oleracea, 2n = 12) and sugar beet (Beta vulgaris, 2n = 18) are important crop members of the family Chenopodiaceae ss Sugar beet has a basic chromosome number of 9 and a cosexual breeding system, as do most members of the Chenopodiaceae ss. family. By contrast, spinach has a basic chromosome number of 6 and, although certain cultivars and genotypes produce monoecious plants, is considered to be a dioecious species. The loci determining male and monoecious sexual expression were mapped to different loci on the spinach sex chromosomes. In this study, a linkage map with 46 mapped protein-coding sequences was constructed for the spinach sex chromosomes. Comparison of the linkage map with a reference genome sequence of sugar beet revealed that the spinach sex chromosomes exhibited extensive synteny with sugar beet chromosomes 4 and 9. Tightly linked protein-coding genes linked to the male-determining locus in spinach corresponded to genes located in or around the putative pericentromeric and centromeric regions of sugar beet chromosomes 4 and 9, supporting the observation that recombination rates were low in the vicinity of the male-determining locus. The locus for monoecism was confined to a chromosomal segment corresponding to a region of approximately 1.7Mb on sugar beet chromosome 9, which may facilitate future positional cloning of the locus.

A major QTL corresponding to the Rk locus for resistance to root-knot nematodes in cowpea (Vigna unguiculata L. Walp.)

Genome resolution of a major QTL associated with the Rk locus in cowpea for resistance to root-knot nematodes has significance for plant breeding programs and R gene characterization. Cowpea (Vigna unguiculata L. Walp.) is a susceptible host of root-knot nematodes (Meloidogyne spp.) (RKN), major plant-parasitic pests in global agriculture. To date, breeding for host resistance in cowpea has relied on phenotypic selection which requires time-consuming and expensive controlled infection assays. To facilitate marker-based selection, we aimed to identify and map quantitative trait loci (QTL) conferring the resistance trait. One recombinant inbred line (RIL) and two F2:3 populations, each derived from a cross between a susceptible and a resistant parent, were genotyped with genome-wide single nucleotide polymorphism (SNP) markers. The populations were screened in the field for root-galling symptoms and/or under growth-chamber conditions for nematode reproduction levels using M. incognita and M. javanica biotypes. One major QTL was mapped consistently on linkage group VuLG11 of each population. By genotyping additional cowpea lines and near-isogenic lines derived from conventional backcrossing, we confirmed that the detected QTL co-localized with the genome region associated with the Rk locus for RKN resistance that has been used in conventional breeding for many decades. This chromosomal location defined with flanking markers will be a valuable target in marker-assisted breeding and for positional cloning of genes controlling RKN resistance.

SBMDb: first whole genome putative microsatellite DNA marker database of sugarbeet for bioenergy and industrial applications

DNA marker plays important role as valuable tools to increase crop productivity by finding plausible answers to genetic variations and linking the Quantitative Trait Loci (QTL) of beneficial trait. Prior approaches in development of Short Tandem Repeats (STR) markers were time consuming and inefficient. Recent methods invoking the development of STR markers using whole genomic or transcriptomics data has gained wide importance with immense potential in developing breeding and cultivator improvement approaches. Availability of whole genome sequences and in silico approaches has revolutionized bulk marker discovery. We report world's first sugarbeet whole genome marker discovery having 145 K markers along with 5 K functional domain markers unified in common platform using MySQL, Apache and PHP in SBMDb. Embedded markers and corresponding location information can be selected for desired chromosome, location/interval and primers can be generated using Primer3 core, integrated at backend. Our analyses revealed abundance of 'mono' repeat (76.82%) over 'di' repeats (13.68%). Highest density (671.05 markers/Mb) was found in chromosome 1 and lowest density (341.27 markers/Mb) in chromosome 6. Current investigation of sugarbeet genome marker density has direct implications in increasing mapping marker density. This will enable present linkage map having marker distance of ∼2 cM, i.e. from 200 to 2.6 Kb, thus facilitating QTL/gene mapping. We also report e-PCR-based detection of 2027 polymorphic markers in panel of five genotypes. These markers can be used for DUS test of variety identification and MAS/GAS in variety improvement program. The present database presents wide source of potential markers for developing and implementing new approaches for molecular breeding required to accelerate industrious use of this crop, especially for sugar, health care products, medicines and color dye. Identified markers will also help in improvement of bioenergy trait of bioethanol and biogas production along with reaping advantage of crop efficiency in terms of low water and carbon footprint especially in era of climate change. Database URL: http://webapp.cabgrid.res.in/sbmdb/.

Exploiting single-molecule transcript sequencing for eukaryotic gene prediction

We develop a method to predict and validate gene models using PacBio single-molecule, real-time (SMRT) cDNA reads. Ninety-eight percent of full-insert SMRT reads span complete open reading frames. Gene model validation using SMRT reads is developed as automated process. Optimized training and prediction settings and mRNA-seq noise reduction of assisting Illumina reads results in increased gene prediction sensitivity and precision. Additionally, we present an improved gene set for sugar beet (Beta vulgaris) and the first genome-wide gene set for spinach (Spinacia oleracea). The workflow and guidelines are a valuable resource to obtain comprehensive gene sets for newly sequenced genomes of non-model eukaryotes.

Lineage-specific gene radiations underlie the evolution of novel betalain pigmentation in Caryophyllales

Betalain pigments are unique to the Caryophyllales and structurally and biosynthetically distinct from anthocyanins. Two key enzymes within the betalain synthesis pathway have been identified: 4,5-dioxygenase (DODA) that catalyzes the formation of betalamic acid and CYP76AD1, a cytochrome P450 gene that catalyzes the formation of cyclo-DOPA. We performed phylogenetic analyses to reveal the evolutionary history of the DODA and CYP76AD1 lineages and in the context of an ancestral reconstruction of pigment states we explored the evolution of these genes in relation to the complex evolution of pigments in Caryophylalles. Duplications within the CYP76AD1 and DODA lineages arose just before the origin of betalain pigmentation in the core Caryophyllales. The duplications gave rise to DODA-α and CYP76AD1-α isoforms that appear specific to betalain synthesis. Both betalain-specific isoforms were then lost or downregulated in the anthocyanic Molluginaceae and Caryophyllaceae. Our findings suggest a single origin of the betalain synthesis pathway, with neofunctionalization following gene duplications in the CYP76AD1 and DODA lineages. Loss of DODA-α and CYP76AD1-α in anthocyanic taxa suggests that betalain pigmentation has been lost twice in Caryophyllales, and exclusion of betalain pigments from anthocyanic taxa is mediated through gene loss or downregulation. [Correction added after online publication 13 May 2015: in the last two paragraphs of the Summary the gene name CYP761A was changed to CYP76AD1.].

Transcriptome Analysis of Beta macrocarpa and Identification of Differentially Expressed Transcripts in Response to Beet Necrotic Yellow Vein Virus Infection

BACKGROUND

Rhizomania is one of the most devastating diseases of sugar beet. It is caused by Beet necrotic yellow vein virus (BNYVV) transmitted by the obligate root-infecting parasite Polymyxa betae. Beta macrocarpa, a wild beet species widely used as a systemic host in the laboratory, can be rub-inoculated with BNYVV to avoid variation associated with the presence of the vector P. betae. To better understand disease and resistance between beets and BNYVV, we characterized the transcriptome of B. macrocarpa and analyzed global gene expression of B. macrocarpa in response to BNYVV infection using the Illumina sequencing platform.

RESULTS

The overall de novo assembly of cDNA sequence data generated 75,917 unigenes, with an average length of 1054 bp. Based on a BLASTX search (E-value ≤ 10-5) against the non-redundant (NR, NCBI) protein, Swiss-Prot, the Gene Ontology (GO), Clusters of Orthologous Groups of proteins (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, there were 39,372 unigenes annotated. In addition, 4,834 simple sequence repeats (SSRs) were also predicted, which could serve as a foundation for various applications in beet breeding. Furthermore, comparative analysis of the two transcriptomes revealed that 261 genes were differentially expressed in infected compared to control plants, including 128 up- and 133 down-regulated genes. GO analysis showed that the changes in the differently expressed genes were mainly enrichment in response to biotic stimulus and primary metabolic process.

CONCLUSION

Our results not only provide a rich genomic resource for beets, but also benefit research into the molecular mechanisms of beet- BNYV Vinteraction.

Dissecting Molecular Evolution in the Highly Diverse Plant Clade Caryophyllales Using Transcriptome Sequencing

Many phylogenomic studies based on transcriptomes have been limited to “single-copy” genes due to methodological challenges in homology and orthology inferences. Only a relatively small number of studies have explored analyses beyond reconstructing species relationships. We sampled 69 transcriptomes in the hyperdiverse plant clade Caryophyllales and 27 outgroups from annotated genomes across eudicots. Using a combined similarity- and phylogenetic tree-based approach, we recovered 10,960 homolog groups, where each was represented by at least eight ingroup taxa. By decomposing these homolog trees, and taking gene duplications into account, we obtained 17,273 ortholog groups, where each was represented by at least ten ingroup taxa. We reconstructed the species phylogeny using a 1,122-gene data set with a gene occupancy of 92.1%. From the homolog trees, we found that both synonymous and nonsynonymous substitution rates in herbaceous lineages are up to three times as fast as in their woody relatives. This is the first time such a pattern has been shown across thousands of nuclear genes with dense taxon sampling. We also pinpointed regions of the Caryophyllales tree that were characterized by relatively high frequencies of gene duplication, including three previously unrecognized whole-genome duplications. By further combining information from homolog tree topology and synonymous distance between paralog pairs, phylogenetic locations for 13 putative genome duplication events were identified. Genes that experienced the greatest gene family expansion were concentrated among those involved in signal transduction and oxidoreduction, including a cytochrome P450 gene that encodes a key enzyme in the betalain synthesis pathway. Our approach demonstrates a new approach for functional phylogenomic analysis in nonmodel species that is based on homolog groups in addition to inferred ortholog groups.

Reliable in silico identification of sequence polymorphisms and their application for extending the genetic map of sugar beet (Beta vulgaris)

Molecular markers are a highly valuable tool for creating genetic maps. Like in many other crops, sugar beet (Beta vulgaris L.) breeding is increasingly supported by the application of such genetic markers. Single nucleotide polymorphism (SNP) based markers have a high potential for automated analysis and high-throughput genotyping. We developed a bioinformatics workflow that uses Sanger and 2nd-generation sequence data for detection, evaluation and verification of new transcript-associated SNPs from sugar beet. RNAseq data from one parent of an established mapping population were produced by 454-FLX sequencing and compared to Sanger ESTs derived from the other parent. The workflow established for SNP detection considers the quality values of both types of reads, provides polymorphic alignments as well as selection criteria for reliable SNP detection and allows painless generation of new genetic markers within genes. We obtained a total of 14,323 genic SNPs and InDels. According to empirically optimised settings for the quality parameters, we classified these SNPs into four usability categories. Validation of a subset of the in silico detected SNPs by genotyping the mapping population indicated a high success rate of the SNP detection. Finally, a total of 307 new markers were integrated with existing data into a new genetic map of sugar beet which offers improved resolution and the integration of terminal markers.

Genome-wide identification and characterisation of R2R3-MYB genes in sugar beet (Beta vulgaris)

BACKGROUND

The R2R3-MYB genes comprise one of the largest transcription factor gene families in plants, playing regulatory roles in plant-specific developmental processes, metabolite accumulation and defense responses. Although genome-wide analysis of this gene family has been carried out in some species, the R2R3-MYB genes in Beta vulgaris ssp. vulgaris (sugar beet) as the first sequenced member of the order Caryophyllales, have not been analysed heretofore.

RESULTS

We present a comprehensive, genome-wide analysis of the MYB genes from Beta vulgaris ssp. vulgaris (sugar beet) which is the first species of the order Caryophyllales with a sequenced genome. A total of 70 R2R3-MYB genes as well as genes encoding three other classes of MYB proteins containing multiple MYB repeats were identified and characterised with respect to structure and chromosomal organisation. Also, organ specific expression patterns were determined from RNA-seq data. The R2R3-MYB genes were functionally categorised which led to the identification of a sugar beet-specific clade with an atypical amino acid composition in the R3 domain, putatively encoding betalain regulators. The functional classification was verified by experimental confirmation of the prediction that the R2R3-MYB gene Bv_iogq encodes a flavonol regulator.

CONCLUSIONS

This study provides the first step towards cloning and functional dissection of the role of MYB transcription factor genes in the nutritionally and evolutionarily interesting species B. vulgaris. In addition, it describes the flavonol regulator BvMYB12, being the first sugar beet R2R3-MYB with an experimentally proven function.

The B2 flowering time locus of beet encodes a zinc finger transcription factor

Sugar beet (Beta vulgaris) is a biennial root crop that grows vegetatively in the first year and starts shoot elongation (bolting) and flowering after exposure to cold temperatures over winter. Early bolting before winter is controlled by the dominant allele of the B locus. Recently, the BOLTING time control 1 (BTC1) gene has been cloned from this locus. BTC1 promotes early bolting through repression of the downstream bolting repressor B. vulgaris flowering locus T1 (BvFT1) and activation of the downstream floral activator BvFT2. We have identified a new bolting locus B2 acting epistatically to B. B2 houses a transcription factor which is diurnally regulated and acts like BTC1 upstream of BvFT1 and BvFT2. It was termed BvBBX19 according to its closest homolog from Arabidopsis thaliana. The encoded protein has two conserved domains with homology to zinc finger B-boxes. Ethyl methanesulfonate-induced mutations within the second B-box caused up-regulation of BvFT1 and complete down-regulation of BvFT2. In Arabidopsis, the expression of FT is promoted by the B-box containing protein CONSTANS (CO). We performed a phylogenetic analysis with B-box genes from beet and A. thaliana but only BvCOL1 clustered with CO. However, BvCOL1 had been excluded as a CO ortholog by previous studies. Therefore, a new model for flowering induction in beet is proposed in which BTC1 and BvBBX19 complement each other and thus acquire a CO function to regulate their downstream targets BvFT1 and BvFT2.

A sequence-ready physical map of barley anchored genetically by two million single-nucleotide polymorphisms

Barley (Hordeum vulgare) is an important cereal crop and a model species for Triticeae genomics. To lay the foundation for hierarchical map-based sequencing, a genome-wide physical map of its large and complex 5.1 billion-bp genome was constructed by high-information content fingerprinting of almost 600,000 bacterial artificial chromosomes representing 14-fold haploid genome coverage. The resultant physical map comprises 9,265 contigs with a cumulative size of 4.9 Gb representing 96% of the physical length of the barley genome. The reliability of the map was verified through extensive genetic marker information and the analysis of topological networks of clone overlaps. A minimum tiling path of 66,772 minimally overlapping clones was defined that will serve as a template for hierarchical clone-by-clone map-based shotgun sequencing. We integrated whole-genome shotgun sequence data from the individuals of two mapping populations with published bacterial artificial chromosome survey sequence information to genetically anchor the physical map. This novel approach in combination with the comprehensive whole-genome shotgun sequence data sets allowed us to independently validate and improve a previously reported physical and genetic framework. The resources developed in this study will underpin fine-mapping and cloning of agronomically important genes and the assembly of a draft genome sequence.

The genome of the recently domesticated crop plant sugar beet (Beta vulgaris)

Sugar beet (Beta vulgaris ssp. vulgaris) is an important crop of temperate climates which provides nearly 30% of the world's annual sugar production and is a source for bioethanol and animal feed. The species belongs to the order of Caryophylalles, is diploid with 2n = 18 chromosomes, has an estimated genome size of 714-758 megabases and shares an ancient genome triplication with other eudicot plants. Leafy beets have been cultivated since Roman times, but sugar beet is one of the most recently domesticated crops. It arose in the late eighteenth century when lines accumulating sugar in the storage root were selected from crosses made with chard and fodder beet. Here we present a reference genome sequence for sugar beet as the first non-rosid, non-asterid eudicot genome, advancing comparative genomics and phylogenetic reconstructions. The genome sequence comprises 567 megabases, of which 85% could be assigned to chromosomes. The assembly covers a large proportion of the repetitive sequence content that was estimated to be 63%. We predicted 27,421 protein-coding genes supported by transcript data and annotated them on the basis of sequence homology. Phylogenetic analyses provided evidence for the separation of Caryophyllales before the split of asterids and rosids, and revealed lineage-specific gene family expansions and losses. We sequenced spinach (Spinacia oleracea), another Caryophyllales species, and validated features that separate this clade from rosids and asterids. Intraspecific genomic variation was analysed based on the genome sequences of sea beet (Beta vulgaris ssp. maritima; progenitor of all beet crops) and four additional sugar beet accessions. We identified seven million variant positions in the reference genome, and also large regions of low variability, indicating artificial selection. The sugar beet genome sequence enables the identification of genes affecting agronomically relevant traits, supports molecular breeding and maximizes the plant's potential in energy biotechnology.

A sugar beet (Beta vulgaris L.) reference FISH karyotype for chromosome and chromosome-arm identification, integration of genetic linkage groups and analysis of major repeat family distribution

We developed a reference karyotype for B. vulgaris which is applicable to all beet cultivars and provides a consistent numbering of chromosomes and genetic linkage groups. Linkage groups of sugar beet were assigned to physical chromosome arms by FISH (fluorescent in situ hybridization) using a set of 18 genetically anchored BAC (bacterial artificial chromosome) markers. Genetic maps of sugar beet were correlated to chromosome arms, and North-South orientation of linkage groups was established. The FISH karyotype provides a technical platform for genome studies and can be applied for numbering and identification of chromosomes in related wild beet species. The discrimination of all nine chromosomes by BAC probes enabled the study of chromosome-specific distribution of the major repetitive components of sugar beet genome comprising pericentromeric, intercalary and subtelomeric satellites and 18S-5.8S-25S and 5S rRNA gene arrays. We developed a multicolor FISH procedure allowing the identification of all nine sugar beet chromosome pairs in a single hybridization using a pool of satellite DNA probes. Fiber-FISH was applied to analyse five chromosome arms in which the furthermost genetic marker of the linkage group was mapped adjacently to terminal repetitive sequences on pachytene chromosomes. Only on two arms telomere arrays and the markers are physically linked, hence these linkage groups can be considered as terminally closed making the further identification of distal informative markers difficult. The results support genetic mapping by marker localization, the anchoring of contigs and scaffolds for the annotation of the sugar beet genome sequence and the analysis of the chromosomal distribution patterns of major families of repetitive DNA.

Evolutionary reshuffling in the Errantivirus lineage Elbe within the Beta vulgaris genome

LTR retrotransposons and retroviruses are closely related. Although a viral envelope gene is found in some LTR retrotransposons and all retroviruses, only the latter show infectivity. The identification of Ty3-gypsy-like retrotransposons possessing putative envelope-like open reading frames blurred the taxonomical borders and led to the establishment of the Errantivirus, Metavirus and Chromovirus genera within the Metaviridae. Only a few plant Errantiviruses have been described, and their evolutionary history is not well understood. In this study, we investigated 27 retroelements of four abundant Elbe retrotransposon families belonging to the Errantiviruses in Beta vulgaris (sugar beet). Retroelements of the Elbe lineage integrated between 0.02 and 5.59 million years ago, and show family-specific variations in autonomy and degree of rearrangements: while Elbe3 members are highly fragmented, often truncated and present in a high number of solo LTRs, Elbe2 members are mainly autonomous. We observed extensive reshuffling of structural motifs across families, leading to the formation of new retrotransposon families. Elbe retrotransposons harbor a typical envelope-like gene, often encoding transmembrane domains. During the course of Elbe evolution, the additional open reading frames have been strongly modified or independently acquired. Taken together, the Elbe lineage serves as retrotransposon model reflecting the various stages in Errantivirus evolution, and allows a detailed analysis of retrotransposon family formation.

Unusual and typical features of a novel restorer-of-fertility gene of sugar beet (Beta vulgaris L.)

Male gametogenesis in plants can be impaired by an incompatibility between nuclear and mitochondrial genomes, termed cytoplasmic male sterility (CMS). A sterilizing factor resides in mitochondria, whereas a nuclear factor, Restorer-of-fertility (Rf), restores male fertility. Although a majority of plant Rf genes are thought to encode a family of RNA-binding proteins called pentatrico-peptide repeat (PPR) proteins, we isolated a novel type of Rf from sugar beet. Two BACs and one cosmid clone that constituted a 383-kbp contig covering the sugar beet Rf1 locus were sequenced. Of 41 genes borne by the contig, quadruplicated genes were found to be associated with specific transcripts in Rf1 flower buds. The quadruplicated genes encoded a protein resembling OMA1, a protein known from yeast and mammals to be involved in mitochondrial protein quality control. Construction of transgenic plants revealed that one of the four genes (bvORF20) was capable of restoring partial pollen fertility to CMS sugar beet; the level of restoration was comparable to that evaluated by a crossing experiment. However, the other genes lacked such a capability. A GFP-fusion experiment showed that bvORF20 encoded a mitochondrial protein. The corresponding gene was cloned from rf1rf1 sugar beet and sequenced, and a solitary gene that was similar but not identical to bvORF20 was found. Genetic features exhibited by sugar beet Rf1, such as gene clustering and copy-number variation between Rf1 and rf, were reminiscent of PPR-type Rf, suggesting that a common evolutionary mechanism(s) operates on plant Rfs irrespective of the translation product.

The role of a pseudo-response regulator gene in life cycle adaptation and domestication of beet

Life cycle adaptation to latitudinal and seasonal variation in photoperiod and temperature is a major determinant of evolutionary success in flowering plants. Whereas the life cycle of the dicotyledonous model species Arabidopsis thaliana is controlled by two epistatic genes, FLOWERING LOCUS C and FRIGIDA, three unrelated loci (VERNALIZATION) determine the spring and winter habits of monocotyledonous plants such as temperate cereals. In the core eudicot species Beta vulgaris, whose lineage diverged from that leading to Arabidopsis shortly after the monocot-dicot split 140 million years ago, the bolting locus B is a master switch distinguishing annuals from biennials. Here, we isolated B and show that the pseudo-response regulator gene BOLTING TIME CONTROL 1 (BvBTC1), through regulation of the FLOWERING LOCUS T genes, is absolutely necessary for flowering and mediates the response to both long days and vernalization. Our results suggest that domestication of beets involved the selection of a rare partial loss-of-function BvBTC1 allele that imparts reduced sensitivity to photoperiod that is restored by vernalization, thus conferring bienniality, and illustrate how evolutionary plasticity at a key regulatory point can enable new life cycle strategies.