Combined Omics Approaches Reveal Distinct Mechanisms of Resistance and/or Susceptibility in Sugar Beet Double Haploid Genotypes at Early Stages of Beet Curly Top Virus Infection.

Sugar beet is susceptible to Beet curly top virus (BCTV), which significantly reduces yield and sugar production in the semi-arid growing regions worldwide. Sources of genetic resistance to BCTV is limited and control depends upon insecticide seed treatments with neonicotinoids. Through double haploid production and genetic selection, BCTV resistant breeding lines have been developed. Using BCTV resistant (R) [KDH13; Line 13 and KDH4-9; Line 4] and susceptible (S) [KDH19-17; Line 19] lines, beet leafhopper mediated natural infection, mRNA/sRNA sequencing, and metabolite analyses, potential mechanisms of resistance against the virus and vector were identified. At early infection stages (2- and 6-days post inoculation), examples of differentially expressed genes highly up-regulated in the 'R' lines (vs. 'S') included EL10Ac5g10437 (inhibitor of trypsin and hageman factor), EL10Ac6g14635 (jasmonate-induced protein), EL10Ac3g06016 (ribosome related), EL10Ac2g02812 (probable prolyl 4-hydroxylase 10), etc. Pathway enrichment analysis showed differentially expressed genes were predominantly involved with peroxisome, amino acids metabolism, fatty acid degradation, amino/nucleotide sugar metabolism, etc. Metabolite analysis revealed significantly higher amounts of specific isoflavonoid O-glycosides, flavonoid 8-C glycosides, triterpenoid, and iridoid-O-glycosides in the leaves of the 'R' lines (vs. 'S'). These data suggest that a combination of transcriptional regulation and production of putative antiviral metabolites might contribute to BCTV resistance. In addition, genome divergence among BCTV strains differentially affects the production of small non-coding RNAs (sncRNAs) and small peptides which may potentially affect pathogenicity and disease symptom development.

Leaf Bacteriome in Sugar Beet Shows Differential Response against Beet curly top virus during Resistant and Susceptible Interactions.

Beet curly top virus (BCTV) significantly reduces sugar beet yield in semi-arid production areas. Genetic resistance to BCTV is limited; therefore, identification of additional resistance-associated factors is highly desired. Using 16S rRNA sequencing and BCTV resistant (R) genotypes (KDH13, KDH4-9) along with a susceptible (S) genotype (KDH19-17), we investigated leaf bacteriome changes during BCTV post inoculation (pi). At day 6 (~6-week-old plants), Cyanobacteria were predominant (~90%); whereas, at week 4 (~10-week-old plants) Firmicutes (11-66%), Bacteroidetes (17-26%), and Verrucomicrobia (12-29%) were predominant phyla and genotype dependent. Both Bacteroidetes and Verrucomicrobia, increased post infection only in the R lines. The bacterial genera Brevibacillus increased at 6 dpi, and Akkermansia and Bacteroides at 4 wkpi in the R lines. Linear discriminant analysis effect size (LEfSe) identified potential biomarkers in the R vs. S lines. Functional profiling revealed bacterial enrichment associated with the TCA cycle, polyisoprenoid, and L-methionine biosynthesis pathways only in KDH4-9 at 6 dpi. At 4 wkpi, bacteria associated with tryptophan and palmitate biosynthesis in the R lines, and uridine monophosphate, phosphatidyl glycerol, and phospholipid biosynthesis in the S line, were enriched. Future characterization of bacterial genera with antiviral properties will help establish their use as biocontrol agents/biomarkers against BCTV.

Cell-Wall-Degrading Enzymes-Related Genes Originating from Rhizoctonia solani Increase Sugar Beet Root Damage in the Presence of Leuconostoc mesenteroides.

Sugar beet crown and root rot caused by Rhizoctonia solani is a major yield constraint. Root rot is highly increased when R. solani and Leuconostoc mesenteroides co-infect roots. We hypothesized that the absence of plant cell-wall-degrading enzymes in L. mesenteroides and their supply by R. solani during close contact, causes increased damage. In planta root inoculation with or without cell-wall-degrading enzymes showed greater rot when L. mesenteroides was combined with cellulase (22 mm rot), polygalacturonase (47 mm), and pectin lyase (57 mm) versus these enzymes (0-26 mm), R. solani (20 mm), and L. mesenteroides (13 mm) individually. Carbohydrate analysis revealed increased simpler carbohydrates (namely glucose + galactose, and fructose) in the infected roots versus mock control, possibly due to the degradation of complex cell wall carbohydrates. Expression of R. solani cellulase, polygalacturonase, and pectin lyase genes during root infection corroborated well with the enzyme data. Global mRNAseq analysis identified candidate genes and highly co-expressed gene modules in all three organisms that might be critical in host plant defense and pathogenesis. Targeting R. solani cell-wall-degrading enzymes in the future could be an effective strategy to mitigate root damage during its interaction with L. mesenteroides.

Genetic diversity among cultivated beets (Beta vulgaris) assessed via population-based whole genome sequences.

BACKGROUND: Diversification on the basis of utilization is a hallmark of Beta vulgaris (beet), as well as other crop species. Often, crop improvement and management activities are segregated by crop type, thus preserving unique genome diversity and organization. Full interfertility is typically retained in crosses between these groups and more traits may be accessible if the genetic basis of crop type lineage were known, along with available genetic markers to effect efficient transfer (e.g., via backcrossing). Beta vulgaris L. (2n =18) is a species complex composed of diverged lineages (e.g., crop types), including the familiar table, leaf (chard), fodder, and sugar beet crop types. Using population genetic and statistical methods with whole genome sequence data from pooled samples of 23 beet cultivars and breeding lines, relationships were determined between accessions based on identity-by-state metrics and shared genetic variation among lineages. RESULTS: Distribution of genetic variation within and between crop types showed extensive shared (e.g. non-unique) genetic variation. Lineage specific variation (e.g. apomorphy) within crop types supported a shared demographic history within each crop type, while principal components analysis revealed strong crop type differentiation. Relative contributions of specific chromosomes to genome wide differentiation were ascertained, with each chromosome revealing a different pattern of differentiation with respect to crop type. Inferred population size history for each crop type helped integrate selection history for each lineage, and highlighted potential genetic bottlenecks in the development of cultivated beet lineages. CONCLUSIONS: A complex evolutionary history of cultigroups in Beta vulgaris was demonstrated, involving lineage divergence as a result of selection and reproductive isolation. Clear delineation of crop types was obfuscated by historical gene flow and common ancestry (e.g. admixture and introgression, and sorting of ancestral polymorphism) which served to share genome variation between crop types and, likely, important phenotypic characters. Table beet was well differentiated as a crop type, and shared more genetic variation within than among crop types. The sugar beet group was not quite as well differentiated as the table beet group. Fodder and chard groups were intermediate between table and sugar groups, perhaps the result of less intensive selection for end use.

Nucleotide-binding resistance gene signatures in sugar beet, insights from a new reference genome.

Nucleotide-binding (NB-ARC), leucine-rich-repeat genes (NLRs) account for 60.8% of resistance (R) genes molecularly characterized from plants. NLRs exist as large gene families prone to tandem duplication and transposition, with high sequence diversity among crops and their wild relatives. This diversity can be a source of new disease resistance, but difficulty in distinguishing specific sequences from homologous gene family members hinders characterization of resistance for improving crop varieties. Current genome sequencing and assembly technologies, especially those using long-read sequencing, are improving resolution of repeat-rich genomic regions and clarifying locations of duplicated genes, such as NLRs. Using the conserved NB-ARC domain as a model, 231 tentative NB-ARC loci were identified in a highly contiguous genome assembly of sugar beet, revealing diverged and truncated NB-ARC signatures as well as full-length sequences. The NB-ARC-associated proteins contained NLR resistance gene domains, including TIR, CC and LRR, as well as other integrated domains. Phylogenetic relationships of partial and complete domains were determined, and patterns of physical clustering in the genome were evaluated. Comparison of sugar beet NB-ARC domains to validated R-genes from monocots and eudicots suggested extensive Beta vulgaris-specific subfamily expansions. The NLR landscape in the rhizomania resistance conferring Rz region of Chromosome 3 was characterized, identifying 26 NLR-like sequences spanning 20 MB. This work presents the first detailed view of NLR family composition in a member of the Caryophyllales, builds a foundation for additional disease resistance work in B. vulgaris, and demonstrates an additional nucleic-acid-based method for NLR prediction in non-model plant species.

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.

Select and Sequence of a Segregating Sugar Beet Population Provides Genomic Perspective of Host Resistance to Seedling Rhizoctonia solani Infection.

Understanding the genetic basis of polygenic traits is a major challenge in agricultural species, especially in non-model systems. Select and sequence (SnS) experiments carried out within existing breeding programs provide a means to simultaneously identify the genomic background of a trait while improving the mean phenotype for a population. Using pooled whole genome sequencing (WGS) of selected and unselected bulks derived from a synthetic outcrossing sugar beet population EL57 (PI 663212), which segregates for seedling rhizoctonia resistance, we identified a putative genomic background involved in conditioning a resistance phenotype. Population genomic parameters were estimated to measure fixation (He), genome divergence (F ST ), and allele frequency changes between bulks (DeltaAF). We report on the genome wide patterns of variation resulting from selection and highlight specific genomic features associated with resistance. Expected heterozygosity (He) showed an increased level of fixation in the resistant bulk, indicating a greater selection pressure was applied. In total, 1,311 biallelic loci were detected as significant FST outliers (p < 0.01) in comparisons between the resistant and susceptible bulks. These loci were detected in 206 regions along the chromosomes and contained 275 genes. We estimated changes in allele frequency between bulks resulting from selection for resistance by leveraging the allele frequencies of an unselected bulk. DeltaAF was a more stringent test of selection and recovered 186 significant loci, representing 32 genes, all of which were also detected using FST. Estimates of population genetic parameters and statistical significance were visualized with respect to the EL10.2 physical map and produced a candidate gene list that was enriched for function in cell wall metabolism and plant disease resistance, including pathogen perception, signal transduction, and pathogen response. Specific variation associated with these genes was also reported and represents genetic markers for validation and prediction of resistance to Rhizoctonia. Select and sequence experiments offer a means to characterize the genetic base of sugar beet, inform selection within breeding programs, and prioritize candidate variation for functional studies.

Regulatory Roles of Small Non-coding RNAs in Sugar Beet Resistance Against Beet curly top virus.

Beet curly top virus (BCTV) mediated yield loss in sugar beets is a major problem worldwide. The circular single-stranded DNA virus is transmitted by the beet leafhopper. Genetic sources of BCTV resistance in sugar beet are limited and commercial cultivars rely on chemical treatments versus durable genetic resistance. Phenotypic selection and double haploid production have resulted in sugar beet germplasm (KDH13; 13 and KDH4-9; 4) that are highly resistant to BCTV. The molecular mechanism of resistance to the virus is unknown, especially the role of small non-coding RNAs (sncRNAs) during early plant-viral interaction. Using the resistant lines along with a susceptible line (KDH19-17; 19), we demonstrate the role of sugar beet microRNAs (miRNAs) in BCTV resistance during early infection stages when symptoms are not yet visible. The differentially expressed miRNAs altered the expression of their corresponding target genes such as pyruvate dehydrogenase (EL10Ac1g02046), carboxylesterase (EL10Ac1g01087), serine/threonine protein phosphatase (EL10Ac1g01374), and leucine-rich repeats (LRR) receptor-like (EL10Ac7g17778), that were highly expressed in the resistant lines versus susceptible lines. Pathway enrichment analysis of the miRNA target genes showed an enrichment of genes involved in glycolysis/gluconeogenesis, galactose metabolism, starch, and sucrose metabolism to name a few. Carbohydrate analysis revealed altered glucose, galactose, fructose, and sucrose concentrations in the infected leaves of resistant versus susceptible lines. We also demonstrate differential regulation of BCTV derived sncRNAs in the resistant versus susceptible lines that target sugar beet genes such as LRR (EL10Ac1g01206), 7-deoxyloganetic acid glucosyltransferase (EL10Ac5g12605), and transmembrane emp24 domain containing (EL10Ac6g14074) and altered their expression. In response to viral infection, we found that plant derived miRNAs targeted BCTV capsid protein/replication related genes and showed differences in expression among resistant and susceptible lines. The data presented here demonstrate the contribution of miRNA mediated regulation of metabolic pathways and cross-kingdom RNA interference (RNAi) in sugar beet BCTV resistance.