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

Repeat turnover meets stable chromosomes: repetitive DNA sequences mark speciation and gene pool boundaries in sugar beet and wild beets

Sugar beet and its wild relatives share a base chromosome number of nine and similar chromosome morphologies. Yet, interspecific breeding is impeded by chromosome and sequence divergence that is still not fully understood. Since repetitive DNAs are among the fastest evolving parts of the genome, we investigated, if repeatome innovations and losses are linked to chromosomal differentiation and speciation. We traced genome and chromosome-wide evolution across 13 beet species comprising all sections of the genera Beta and Patellifolia. For this, we combined short and long read sequencing, flow cytometry, and cytogenetics to build a comprehensive framework that spans the complete scale from DNA to chromosome to genome. Genome sizes and repeat profiles reflect the separation into three gene pools with contrasting evolutionary patterns. Among all repeats, satellite DNAs harbor most genomic variability, leading to fundamentally different centromere architectures, ranging from chromosomal uniformity in Beta and Patellifolia to the formation of patchwork chromosomes in Corollinae/Nanae. We show that repetitive DNAs are causal for the genome expansions and contractions across the beet genera, providing insights into the genomic underpinnings of beet speciation. Satellite DNAs in particular vary considerably between beet genomes, leading to the evolution of distinct chromosomal setups in the three gene pools, likely contributing to the barriers in beet breeding. Thus, with their isokaryotypic chromosome sets, beet genomes present an ideal system for studying the link between repeats, genomic variability, and chromosomal differentiation and provide a theoretical fundament for understanding barriers in any crop breeding effort.

Assessment of co-infection with BNYVV and BSCTV on resistance against Rhizomania disease in transgenic sugar beet plants

Sugar beet is an economically important crop and one of the major sources of sucrose around the world. Beet necrotic yellow vein virus (BNYVV) and Beet severe curly top virus (BSCTV) are two widespread viruses in sugar beet that cause severe damage to its performance. Previously, we have successfully produced resistance to BNYVV based on RNA silencing in sugar beet by introducing constructs carrying the viral coat-protein-encoding DNA sequence, CP21, in sense and anti-sense orientations. Yet, the RNA silencing-mediated resistance to a specific virus could be affected by other ones as a part of synergistic interactions. In this study, we assayed the specificity of the induced resistance against BNYVV in two sets of transgenic events, S3 and S6 carrying 5'-UTR with or without CP21-coding sequences, respectively. These events were subjected to viral challenges with either BNYVV, an Iranian isolate of BSCTV (BSCTV-Ir) or both. All the plants inoculated with just BSCTV-Ir displayed curly-leaf symptoms. However, partial resistance was evident in S3 events as shown by mild symptoms and reduced PCR amplification of the BSCTV-Ir coat protein encoding sequence. Based on the presented data, resistance to BNYVV was stable in almost all the transgenic plants co-infected with BSCTV-Ir, except for one event, S3-229. In general, it seems that the co-infection does not affect the resistance to BNYVV in transgenic plants. These findings demonstrated that the introduced RNA silencing-mediated resistance against BNYVV in transgenic sugar beets is specific and is not suppressed after co-infection with a heterologous virus.

The m6 A reader MhYTP2 negatively modulates apple Glomerella leaf spot resistance by binding to and degrading MdRGA2L mRNA

Glomerella leaf spot (GLS), caused by the fungal pathogen Colletotrichum fructicola, significantly threatens apple production. Some resistances to plant disease are mediated by the accumulation of nucleotide-binding site and leucine-rich repeat (NBS-LRR) proteins that are encoded by a major class of plant disease resistance genes (R genes). However, the R genes that confer resistance to GLS in apple remain largely unclear. Malus hupehensis YT521-B homology domain-containing protein 2 (MhYTP2) was identified as an N6 -methyladenosine RNA methylation (m6 A) modified RNA reader in our previous study. However, whether MhYTP2 binds to mRNAs without m6 A RNA modifications remains unknown. In this study, we discovered that MhYTP2 exerts both m6 A-dependent and -independent functions by analysing previously obtained RNA immunoprecipitation sequencing results. The overexpression of MhYTP2 significantly reduced the resistance of apple to GLS and down-regulated the transcript levels of some R genes whose transcripts do not contain m6 A modifications. Further analysis indicated that MhYTP2 binds to and reduces the stability of MdRGA2L mRNA. MdRGA2L positively regulates resistance to GLS by activating salicylic acid signalling. Our findings revealed that MhYTP2 plays an essential role in the regulation of resistance to GLS and identified a promising R gene, MdRGA2L, for use in developing apple cultivars with GLS resistance.

Nuclear DNA segments homologous to mitochondrial DNA are obstacles for detecting heteroplasmy in sugar beet (Beta vulgaris L.)

Heteroplasmy, the coexistence of multiple mitochondrial DNA (mtDNA) sequences in a cell, is well documented in plants. Next-generation sequencing technology (NGS) has made it feasible to sequence entire genomes. Thus, NGS has the potential to detect heteroplasmy; however, the methods and pitfalls in heteroplasmy detection have not been fully investigated and identified. One obstacle for heteroplasmy detection is the sequence homology between mitochondrial-, plastid-, and nuclear DNA, of which the influence of nuclear DNA segments homologous to mtDNA (numt) need to be minimized. To detect heteroplasmy, we first excluded nuclear DNA sequences of sugar beet (Beta vulgaris) line EL10 from the sugar beet mtDNA sequence. NGS reads were obtained from single plants of sugar beet lines NK-195BRmm-O and NK-291BRmm-O and mapped to the unexcluded mtDNA regions. More than 1000 sites exhibited intra-individual polymorphism as detected by genome browsing analysis. We focused on a 309-bp region where 12 intra-individual polymorphic sites were closely linked to each other. Although the existence of DNA molecules having variant alleles at the 12 sites was confirmed by PCR amplification from NK-195BRmm-O and NK-291BRmm-O, these variants were not always called by six variant-calling programs, suggesting that these programs are inappropriate for intra-individual polymorphism detection. When we changed the nuclear DNA reference, a numt absent from EL10 was found to include the 309-bp region. Genetic segregation of an F2 population from NK-195BRmm-O x NK-291BRmm-O supported the numt origin of the variant alleles. Using four references, we found that numt detection exhibited reference dependency, and extreme polymorphism of numts exists among sugar beet lines. One of the identified numts absent from EL10 is also associated with another intra-individual polymorphic site in NK-195mm-O. Our data suggest that polymorphism among numts is unexpectedly high within sugar beets, leading to confusion about the true degree of heteroplasmy.

Engineering Resistance against Sclerotinia sclerotiorum Using a Truncated NLR (TNx) and a Defense-Priming Gene

The association of both cell-surface PRRs (Pattern Recognition Receptors) and intracellular receptor NLRs (Nucleotide-Binding Leucine-Rich Repeat) in engineered plants have the potential to activate strong defenses against a broad range of pathogens. Here, we describe the identification, characterization, and in planta functional analysis of a novel truncated NLR (TNx) gene from the wild species Arachis stenosperma (AsTIR19), with a protein structure lacking the C-terminal LRR (Leucine Rich Repeat) domain involved in pathogen perception. Overexpression of AsTIR19 in tobacco plants led to a significant reduction in infection caused by Sclerotinia sclerotiorum, with a further reduction in pyramid lines containing an expansin-like B gene (AdEXLB8) potentially involved in defense priming. Transcription analysis of tobacco transgenic lines revealed induction of hormone defense pathways (SA; JA-ET) and PRs (Pathogenesis-Related proteins) production. The strong upregulation of the respiratory burst oxidase homolog D (RbohD) gene in the pyramid lines suggests its central role in mediating immune responses in plants co-expressing the two transgenes, with reactive oxygen species (ROS) production enhanced by AdEXLB8 cues leading to stronger defense response. Here, we demonstrate that the association of potential priming elicitors and truncated NLRs can produce a synergistic effect on fungal resistance, constituting a promising strategy for improved, non-specific resistance to plant pathogens.

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.

A TIR-NBS-LRR Gene MdTNL1 Regulates Resistance to Glomerella Leaf Spot in Apple

Glomerella leaf spot (GLS), caused by the fungus Colletotrichum fructicola, is one of the most devastating apple diseases. Our previous study reported that the GLS resistance locus was defined on the chromosome 15 region. Here, we further found a single-nucleotide polymorphism (SNP) site (SNP_7309212) in the GLS resistance that was able to distinguish resistant cultivars (lines) from susceptible ones. On the basis of the SNP site, we cloned a TNL gene from the GLS resistant locus and named it MdTNL1 (NCBI Accession Number: ON402514). This gene contains a toll/interleukin-1 receptor transmembrane domain (TIR), nucleotide-binding sites (NBS), and leucine-rich repeat (LRR) domain. Subcellular location indicated that MdTNL1 was expressed in the nucleus and cell membrane. Ectopic overexpression of MdTNL1 in Nicotiana benthamiana caused cell death. We further demonstrated allelic polymorphisms in MdTNL1. It is noteworthy that NBS and LRR domains of the MdTNL1 protein serve as the repository for generating allelic diversity. Quantitative real-time PCR (qRT-PCR) assay revealed that MdTNL1 was highly expressed in resistant apple cultivar ‘Fuji’ after inoculation with C. fructicola, whereas susceptible cultivar ‘Golden Delicious’ exhibited low expression after inoculation. Over-expression of MdTNL1-1 in susceptible apple fruits and leaves improved disease resistance, while in ‘Orin’ calli, silencing the MdTNL1-1 gene conversely decreased GLS resistance. In conclusion, we identified a GLS associated with SNP_7309212 and demonstrated that a TIR-NBS-LRR gene MdTNL1-1 positively regulates GLS resistance in apple.

Prediction of NB-LRR resistance genes based on full-length sequence homology

The activation of plant immunity is mediated by resistance (R)-gene receptors, also known as nucleotide-binding leucine-rich repeat (NB-LRR) genes, which in turn trigger the authentic defense response. R-gene identification is a crucial goal for both classic and modern plant breeding strategies for disease resistance. The conventional method identifies NB-LRR genes using a protein motif/domain-based search (PDS) within an automatically predicted gene set of the respective genome assembly. PDS proved to be imprecise since repeat masking prior to automatic genome annotation unwittingly prevented comprehensive NB-LRR gene detection. Furthermore, R-genes have diversified in a species-specific manner, so that NB-LRR gene identification cannot be universally standardized. Here, we present the full-length Homology-based R-gene Prediction (HRP) method for the comprehensive identification and annotation of a genome's R-gene repertoire. Our method has substantially addressed the complex genomic organization of tomato (Solanum lycopersicum) NB-LRR gene loci, proving to be more performant than the well-established RenSeq approach. HRP efficiency was also tested on three differently assembled and annotated Beta sp. genomes. Indeed, HRP identified up to 45% more full-length NB-LRR genes compared to previous approaches. HRP also turned out to be a more refined strategy for R-gene allele mining, testified by the identification of hitherto undiscovered Fom-2 homologs in five Cucurbita sp. genomes. In summary, our high-performance method for full-length NB-LRR gene discovery will propel the identification of novel R-genes towards development of improved cultivars.

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.

ECCsplorer: a pipeline to detect extrachromosomal circular DNA (eccDNA) from next-generation sequencing data

Background

Extrachromosomal circular DNAs (eccDNAs) are ring-like DNA structures physically separated from the chromosomes with 100 bp to several megabasepairs in size. Apart from carrying tandemly repeated DNA, eccDNAs may also harbor extra copies of genes or recently activated transposable elements. As eccDNAs occur in all eukaryotes investigated so far and likely play roles in stress, cancer, and aging, they have been prime targets in recent research—with their investigation limited by the scarcity of computational tools.

Results

Here, we present the ECCsplorer, a bioinformatics pipeline to detect eccDNAs in any kind of organism or tissue using next-generation sequencing techniques. Following Illumina-sequencing of amplified circular DNA (circSeq), the ECCsplorer enables an easy and automated discovery of eccDNA candidates. The data analysis encompasses two major procedures: first, read mapping to the reference genome allows the detection of informative read distributions including high coverage, discordant mapping, and split reads. Second, reference-free comparison of read clusters from amplified eccDNA against control sample data reveals specifically enriched DNA circles. Both software parts can be run separately or jointly, depending on the individual aim or data availability. To illustrate the wide applicability of our approach, we analyzed semi-artificial and published circSeq data from the model organisms Homo sapiens and Arabidopsis thaliana, and generated circSeq reads from the non-model crop plant Beta vulgaris. We clearly identified eccDNA candidates from all datasets, with and without reference genomes. The ECCsplorer pipeline specifically detected mitochondrial mini-circles and retrotransposon activation, showcasing the ECCsplorer’s sensitivity and specificity.

Conclusion

The ECCsplorer (available online at https://github.com/crimBubble/ECCsplorer) is a bioinformatics pipeline to detect eccDNAs in any kind of organism or tissue using next-generation sequencing data. The derived eccDNA targets are valuable for a wide range of downstream investigations—from analysis of cancer-related eccDNAs over organelle genomics to identification of active transposable elements.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04545-2.

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.

Development of an SNP Assay for Marker-Assisted Selection of Soil-Borne Rhizoctonia solani AG-2-2-IIIB Resistance in Sugar Beet

Simple Summary

Sustainable breeding of sugar beet against Rhizoctonia solani relies on the continuous identification of resistance genes to allow their integration into new and modern cultivars. Better control of the disease may thus be achieved by a combination of tolerant or resistant cultivars selected based on molecular markers such as SNPs. The utility of one such marker, RsBv1 (Chromosome 6, 9,000,093 bp, C/T), located in an ADP-ribosylation factor and associated with Rhizoctonia resistance resulting from validation of three geographically diverse plant materials is reported.

Abstract

Rhizoctonia solani, causing Rhizoctonia crown and root rot, is a major risk to sugar beet (Beta vulgaris L.) cultivation. The development of resistant varieties accelerated by marker-assisted selection is a priority of breeding programs. We report the identification of a single-nucleotide polymorphism (SNP) marker linked to Rhizoctonia resistance using restriction site-associated DNA (RAD) sequencing of two geographically discrete sets of plant materials with different degrees of resistance/susceptibility to enable a wider selection of superior genotypes. The variant calling pipeline utilized SAMtools for variant calling and the resulting raw SNPs from RAD sequencing (15,988 and 22,439 SNPs) were able to explain 13.40% and 25.45% of the phenotypic variation in the two sets of material from different sources of origin, respectively. An association analysis was carried out independently on both the datasets and mutually occurring significant SNPs were filtered depending on their contribution to the phenotype using principal component analysis (PCA) biplots. To provide a ready-to-use marker for the breeding community, a systematic molecular validation of significant SNPs distributed across the genome was undertaken to combine high-resolution melting, Sanger sequencing, and rhAmp SNP genotyping. We report that RsBv1 located on Chromosome 6 (9,000,093 bp) is significantly associated with Rhizoctonia resistance (p < 0.01) and able to explain 10% of the phenotypic disease variance. The related SNP assay is thus ready for marker-assisted selection in sugar beet breeding for Rhizoctonia resistance.

Beta vulgaris ssp. vulgaris chromosome 8 shows significant association with geosmin concentration in table beet

Geosmin, a degraded sesquiterpene molecule with earthy and musty odor, imbues table beet with its characteristic aroma. Geosmin is heritable and endogenously produced in table beet; its earthy aroma is sought by some consumers but deters others. Geosmin biosynthesis is catalyzed by a bifunctional geosmin synthase enzyme in diverse bacteria and fungi, but a mechanism for geosmin biosynthesis in plants has not been reported. This work employed association analysis and selective genotyping of a segregating F2:3 mapping population to seek QTL associated with geosmin concentration in table beet. GBS reads were aligned to sugar beet reference genome EL10.2, and association analysis revealed two QTL for geosmin concentration on Beta vulgaris ssp. vulgaris chromosome 8. QTL at EL10.2 positions 28,017,624 and 38,488,687 each show effect size 8.7 μg·kg-1 geosmin and explain 8.5% and 6.4% of total variation in geosmin concentration, respectively. Resolution was low due to large recombination bin size and imperfect alignment between the reference genome and mapping population, but population size and selection proportion were sufficient to detect moderate to large effect QTL. This study, the first molecular genetic mapping experiment in table beet, succeeded in finding QTL for geosmin concentration in table beet, and it provides the basis for fine mapping or candidate gene investigation of functional loci for this distinctive sensory trait.

Broken, silent, and in hiding: tamed endogenous pararetroviruses escape elimination from the genome of sugar beet (Beta vulgaris)

BACKGROUND AND AIMS

Endogenous pararetroviruses (EPRVs) are widespread components of plant genomes that originated from episomal DNA viruses of the Caulimoviridae family. Due to fragmentation and rearrangements, most EPRVs have lost their ability to replicate through reverse transcription and to initiate viral infection. Similar to the closely related retrotransposons, extant EPRVs were retained and often amplified in plant genomes for several million years. Here, we characterize the complete genomic EPRV fraction of the crop sugar beet (Beta vulgaris, Amaranthaceae) to understand how they shaped the beet genome and to suggest explanations for their absent virulence.

METHODS

Using next- and third-generation sequencing data and genome assembly, we reconstructed full-length in silico representatives for the three host-specific EPRVs (beetEPRVs) in the B. vulgaris genome. Focusing on the endogenous caulimovirid beetEPRV3, we investigated its chromosomal localization, abundance and distribution by fluorescent in situ and Southern hybridization.

KEY RESULTS

Full-length beetEPRVs range between 7.5 and 10.7 kb in size, are heterogeneous in structure and sequence, and occupy about 0.3 % of the beet genome. Although all three beetEPRVs were assigned to the florendoviruses, they showed variably arranged protein-coding domains, different fragmentation, and preferences for diverse sequence contexts. We observed small RNAs that specifically target the individual beetEPRVs, indicating stringent epigenetic suppression. BeetEPRV3 sequences occur along all sugar beet chromosomes, preferentially in the vicinity of each other and are associated with heterochromatic, centromeric and intercalary satellite DNAs. BeetEPRV3 members also exist in genomes of related wild species, indicating an initial beetEPRV3 integration 13.4-7.2 million years ago.

CONCLUSIONS

Our study in beet illustrates the variability of EPRV structure and sequence in a single host genome. Evidence of sequence fragmentation and epigenetic silencing implies possible plant strategies to cope with long-term persistence of EPRVs, including amplification, fixation in the heterochromatin, and containment of EPRV virulence.

A spinach genome assembly with remarkable completeness, and its use for rapid identification of candidate genes for agronomic traits

Spinach (Spinacia oleracea) is grown as a nutritious leafy vegetable worldwide. To accelerate spinach breeding efficiency, a high-quality reference genome sequence with great completeness and continuity is needed as a basic infrastructure. Here, we used long-read and linked-read technologies to construct a de novo spinach genome assembly, designated SOL_r1.1, which was comprised of 287 scaffolds (total size: 935.7 Mb; N50 = 11.3 Mb) with a low proportion of undetermined nucleotides (Ns = 0.34%) and with high gene completeness (BUSCO complete 96.9%). A genome-wide survey of resistance gene analogues identified 695 genes encoding nucleotide-binding site domains, receptor-like protein kinases, receptor-like proteins and transmembrane-coiled coil domains. Based on a high-density double-digest restriction-site associated DNA sequencing-based linkage map, the genome assembly was anchored to six pseudomolecules representing ∼73.5% of the whole genome assembly. In addition, we used SOL_r1.1 to identify quantitative trait loci for bolting timing and fruit/seed shape, which harbour biologically plausible candidate genes, such as homologues of the FLOWERING LOCUS T and EPIDERMAL PATTERNING FACTOR-LIKE genes. The new genome assembly, SOL_r1.1, will serve as a useful resource for identifying loci associated with important agronomic traits and for developing molecular markers for spinach breeding/selection programs.

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.

Assembly and characterization of the genome of chard (Beta vulgaris ssp. vulgaris var. cicla)

Chard (Beta vulgaris ssp. vulgaris var. cicla) is a member of one of four different cultigroups of beets. While the genome of sugar beet, the most prominent beet crop, has been studied extensively, molecular data on other beet cultivars is scant. Here, we present a genome assembly of chard, a vegetable crop grown for its fleshy leaves. We report a de novo genome assembly of 604 Mbp, slightly larger than sugar beet assemblies presented so far. About 57 % of the assembly was annotated as repetitive sequence, of which LTR retrotransposons were the most abundant. Based on the presence of conserved genes, the chard assembly was estimated to be at least 96 % complete regarding its gene space. We predicted 34,521 genes of which 27,582 genes were supported by evidence from transcriptomic sequencing reads, and 5503 of the evidence-supported genes had multiple isoforms. We compared the chard gene set with gene sets from sugar beet and two wild beets (i.e. Beta vulgaris ssp. maritima and Beta patula) to find orthology relationships and identified genome-wide syntenic regions between chard and sugar beet. Lastly, we determined genomic variants that distinguish sugar beet and chard. Assessing the variation distribution along the chard chromosomes, we found extensive haplotype sharing between the two cultivars. In summary, our work provides a foundation for the molecular analysis of Beta vulgaris cultigroups as a basis for chard genomics and to unravel the domestication history of beet crops.

A high-quality Brassica napus genome reveals expansion of transposable elements, subgenome evolution and disease resistance

Rapeseed (Brassica napus L.) is a recent allotetraploid crop, which is well known for its high oil production. Here, we report a high-quality genome assembly of a typical semi-winter rapeseed cultivar, 'Zhongshuang11' (hereafter 'ZS11'), using a combination of single-molecule sequencing and chromosome conformation capture (Hi-C) techniques. Most of the high-confidence sequences (93.1%) were anchored to the individual chromosomes with a total of 19 centromeres identified, matching the exact chromosome count of B. napus. The repeat sequences in the A and C subgenomes in B. napus expanded significantly from 500 000 years ago, especially over the last 100 000 years. These young and recently amplified LTR-RTs showed dispersed chromosomal distribution but significantly preferentially clustered into centromeric regions. We exhaustively annotated the nucleotide-binding leucine-rich repeat (NLR) gene repertoire, yielding a total of 597 NLR genes in B. napus genome and 17.4% of which are paired (head-to-head arrangement). Based on the resequencing data of 991 B. napus accessions, we have identified 18 759 245 single nucleotide polymorphisms (SNPs) and detected a large number of genomic regions under selective sweep among the three major ecotype groups (winter, semi-winter and spring) in B. napus. We found 49 NLR genes and five NLR gene pairs colocated in selective sweep regions with different ecotypes, suggesting a rapid diversification of NLR genes during the domestication of B. napus. The high quality of our B. napus 'ZS11' genome assembly could serve as an important resource for the study of rapeseed genomics and reveal the genetic variations associated with important agronomic traits.

Major latex protein-like encoding genes contribute to Rhizoctonia solani defense responses in sugar beet

Sugar beets are attacked by several pathogens that cause root damages. Rhizoctonia (Greek for “root killer”) is one of them. Rhizoctonia root rot has become an increasing problem for sugar beet production and to decrease yield losses agronomical measures are adopted. Here, two partially resistant and two susceptible sugar beet genotypes were used for transcriptome analysis to discover new defense genes to this fungal disease, information to be implemented in molecular resistance breeding. Among 217 transcripts with increased expression at 2 days post-infection (dpi), three resistance-like genes were found. These genes were not significantly elevated at 5 dpi, a time point when increased expression of three Bet v I/Major latex protein (MLP) homologous genes BvMLP1, BvMLP2 and BvML3 was observed in the partially resistant genotypes. Quantitative RT-PCR analysis on diseased sugar beet seedlings validated the activity of BvMLP1 and BvMLP3 observed in the transcriptome during challenge by R. solani. The three BvMLP genes were cloned and overexpressed in Arabidopsis thaliana to further dissect their individual contribution. Transgenic plants were also compared to T-DNA mutants of orthologous MLP genes. Plants overexpressing BvMLP1 and BvMLP3 showed significantly less infection whereas additive effects were seen on Atmlp1/Atmlp3 double mutants. The data suggest that BvMLP1 and BvMLP3 may contribute to the reduction of the Rhizoctonia root rot disease in sugar beet. Impact on the defense reaction from other differential expressed genes observed in the study is discussed.

Cercospora beticola: The intoxicating lifestyle of the leaf spot pathogen of sugar beet

Cercospora leaf spot, caused by the fungal pathogen Cercospora beticola, is the most destructive foliar disease of sugar beet worldwide. This review discusses C. beticola genetics, genomics, and biology and summarizes our current understanding of the molecular interactions that occur between C. beticola and its sugar beet host. We highlight the known virulence arsenal of C. beticola as well as its ability to overcome currently used disease management strategies. Finally, we discuss future prospects for the study and management of C. beticola infections in the context of newly employed molecular tools to uncover additional information regarding the biology of this pathogen.

TAXONOMY

Cercospora beticola Sacc.; Kingdom Fungi, Phylum Ascomycota, Class Dothideomycetes, Order Capnodiales, Family Mycosphaerellaceae, Genus Cercospora.

HOST RANGE

Well-known pathogen of sugar beet (Beta vulgaris subsp. vulgaris) and most species of the Beta genus. Reported as pathogenic on other members of the Chenopodiaceae (e.g., lamb's quarters, spinach) as well as members of the Acanthaceae (e.g., bear's breeches), Apiaceae (e.g., Apium), Asteraceae (e.g., chrysanthemum, lettuce, safflower), Brassicaceae (e.g., wild mustard), Malvaceae (e.g., Malva), Plumbaginaceae (e.g., Limonium), and Polygonaceae (e.g., broad-leaved dock) families.

DISEASE SYMPTOMS

Leaves infected with C. beticola exhibit circular lesions that are coloured tan to grey in the centre and are often delimited by tan-brown to reddish-purple rings. As disease progresses, spots can coalesce to form larger necrotic areas, causing severely infected leaves to wither and die. At the centre of these spots are black spore-bearing structures (pseudostromata). Older leaves often show symptoms first and younger leaves become infected as the disease progresses.

MANAGEMENT

Application of a mixture of fungicides with different modes of action is currently performed although elevated resistance has been documented in most employed fungicide classes. Breeding for high-yielding cultivars with improved host resistance is an ongoing effort and prudent cultural practices, such as crop rotation, weed host management, and cultivation to reduce infested residue levels, are widely used to manage disease. USEFUL WEBSITE: https://www.ncbi.nlm.nih.gov/genome/11237?genome_assembly_id=352037.

A Chromosome-Scale Assembly of the Garden Orach (Atriplex hortensis L.) Genome Using Oxford Nanopore Sequencing

Atriplex hortensis (2n = 2x = 18, 1C genome size ∼1.1 gigabases), also known as garden orach and mountain-spinach, is a highly nutritious, broadleaf annual of the Amaranthaceae-Chenopodiaceae alliance (Chenopodiaceae sensu stricto, subfam. Chenopodioideae) that has spread in cultivation from its native primary domestication area in Eurasia to other temperate and subtropical regions worldwide. Atriplex L. is a highly complex but, as understood now, a monophyletic group of mainly halophytic and/or xerophytic plants, of which A. hortensis has been a vegetable of minor importance in some areas of Eurasia (from Central Asia to the Mediterranean) at least since antiquity. Nonetheless, it is a crop with tremendous nutritional potential due primarily to its exceptional leaf and seed protein quantities (approaching 30%) and quality (high levels of lysine). Although there is some literature describing the taxonomy and production of A. hortensis, there is a general lack of genetic and genomic data that would otherwise help elucidate the genetic variation, phylogenetic positioning, and future potential of the species. Here, we report the assembly of the first high-quality, chromosome-scale reference genome for A. hortensis cv. "Golden." Long-read data from Oxford Nanopore's MinION DNA sequencer was assembled with the program Canu and polished with Illumina short reads. Contigs were scaffolded to chromosome scale using chromatin-proximity maps (Hi-C) yielding a final assembly containing 1,325 scaffolds with a N50 of 98.9 Mb - with 94.7% of the assembly represented in the nine largest, chromosome-scale scaffolds. Sixty-six percent of the genome was classified as highly repetitive DNA, with the most common repetitive elements being Gypsy-(32%) and Copia-like (11%) long-terminal repeats. The annotation was completed using MAKER which identified 37,083 gene models and 2,555 tRNA genes. Completeness of the genome, assessed using the Benchmarking Universal Single Copy Orthologs (BUSCO) metric, identified 97.5% of the conserved orthologs as complete, with only 2.2% being duplicated, reflecting the diploid nature of A. hortensis. A resequencing panel of 21 wild, unimproved and cultivated A. hortensis accessions revealed three distinct populations with little variation within subpopulations. These resources provide vital information to better understand A. hortensis and facilitate future study.

Integrated annotations and analyses of small RNA-producing loci from 47 diverse plants

Plant endogenous small RNAs (sRNAs) are important regulators of gene expression. There are two broad categories of plant sRNAs: microRNAs (miRNAs) and endogenous short interfering RNAs (siRNAs). MicroRNA loci are relatively well-annotated but compose only a small minority of the total sRNA pool; siRNA locus annotations have lagged far behind. Here, we used a large data set of published and newly generated sRNA sequencing data (1333 sRNA-seq libraries containing more than 20 billion reads) and a uniform bioinformatic pipeline to produce comprehensive sRNA locus annotations of 47 diverse plants, yielding more than 2.7 million sRNA loci. The two most numerous classes of siRNA loci produced mainly 24- and 21-nucleotide (nt) siRNAs, respectively. Most often, 24-nt-dominated siRNA loci occurred in intergenic regions, especially at the 5′-flanking regions of protein-coding genes. In contrast, 21-nt-dominated siRNA loci were most often derived from double-stranded RNA precursors copied from spliced mRNAs. Genic 21-nt-dominated loci were especially common from disease resistance genes, including from a large number of monocots. Individual siRNA sequences of all types showed very little conservation across species, whereas mature miRNAs were more likely to be conserved. We developed a web server where our data and several search and analysis tools are freely accessible.

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.

How did a duplicated gene copy evolve into a restorer-of-fertility gene in a plant? The case of Oma1

Restorer-of-fertility (Rf) is a suppressor of cytoplasmic male sterility (CMS), a mitochondrion-encoded trait that has been reported in many plant species. The occurrence of CMS is considered to be independent in each lineage; hence, the question of how Rf evolved was raised. Sugar beet Rf resembles Oma1, a gene for quality control of the mitochondrial inner membrane. Oma1 homologues comprise a small gene family in the sugar beet genome, unlike Arabidopsis and other eukaryotes. The sugar beet sequence that best matched Arabidopsis atOma1 was named bvOma1; sugar beet Rf (RF1-Oma1) was another member. During anther development, atOma1 mRNA was detected from the tetrad to the microspore stages, whereas bvOma1 mRNA was detected at the microspore stage and RF1-Oma1 mRNA was detected during the meiosis and tetrad stages. A transgenic study revealed that, whereas RF1-Oma1 can bind to a CMS-specific protein and alter the higher-order structure of the CMS-specific protein complex, neither bvOma1 nor atOma1 show such activity. We favour the hypothesis that an ancestral Oma1 gene duplicated to form a small gene family, and that one of the copies evolved and acquired a novel expression pattern and protein function as an Rf, i.e. RF1-Oma1 evolved via neofunctionalization.

The genome of Chenopodium pallidicaule: An emerging Andean super grain

PREMISE

Cañahua is a semi-domesticated crop grown in high-altitude regions of the Andes. It is an A-genome diploid (2n = 2x = 18) relative of the allotetraploid (AABB) Chenopodium quinoa and shares many of its nutritional benefits. Cañahua seed contains a complete protein, a low glycemic index, and offers a wide variety of nutritionally important vitamins and minerals.

METHODS

The reference assembly was developed using a combination of short- and long-read sequencing techniques, including multiple rounds of Hi-C-based proximity-guided assembly.

RESULTS

The final assembly of the ~363-Mbp genome consists of 4633 scaffolds, with 96.6% of the assembly contained in nine scaffolds representing the nine haploid chromosomes of the species. Repetitive element analysis classified 52.3% of the assembly as repetitive, with the most common repeat identified as long terminal repeat retrotransposons. MAKER annotation of the final assembly yielded 22,832 putative gene models.

DISCUSSION

When compared with quinoa, strong patterns of synteny support the hypothesis that cañahua is a close A-genome diploid relative, and thus potentially a simplified model diploid species for genetic analysis and improvement of quinoa. Resequencing and phylogenetic analysis of a diversity panel of cañahua accessions suggests that coordinated efforts are needed to enhance genetic diversity conservation within ex situ germplasm collections.

Disease Resistance Genetics and Genomics in Octoploid Strawberry

Octoploid strawberry (Fragaria ×ananassa) is a valuable specialty crop, but profitable production and availability are threatened by many pathogens. Efforts to identify and introgress useful disease resistance genes (R-genes) in breeding programs are complicated by strawberry’s complex octoploid genome. Recently-developed resources in strawberry, including a complete octoploid reference genome and high-resolution octoploid genotyping, enable new analyses in strawberry disease resistance genetics. This study characterizes the complete R-gene collection in the genomes of commercial octoploid strawberry and two diploid ancestral relatives, and introduces several new technological and data resources for strawberry disease resistance research. These include octoploid R-gene transcription profiling, dN/dS analysis, expression quantitative trait loci (eQTL) analysis and RenSeq analysis in cultivars. Octoploid fruit eQTL were identified for 76 putative R-genes. R-genes from the ancestral diploids Fragaria vesca and Fragaria iinumae were compared, revealing differential inheritance and retention of various octoploid R-gene subtypes. The mode and magnitude of natural selection of individual F. ×ananassa R-genes was also determined via dN/dS analysis. R-gene sequencing using enriched libraries (RenSeq) has been used recently for R-gene discovery in many crops, however this technique somewhat relies upon a priori knowledge of desired sequences. An octoploid strawberry capture-probe panel, derived from the results of this study, is validated in a RenSeq experiment and is presented for community use. These results give unprecedented insight into crop disease resistance genetics, and represent an advance toward exploiting variation for strawberry cultivar improvement.

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.

Genome-wide identification, characterization, and expression analysis of nucleotide-binding leucine-rich repeats gene family under environmental stresses in tea (Camellia sinensis)

Plants often use nucleotide-binding leucine-rich repeats (NLRs) to recognize specific virulence proteins and activate the hypersensitive response thereby defending against invaders. However, data on NLRs and the resistance mechanism of NLR protein mediation in tea plant are extremely limited. In this study, 400 and 303 CsNLRs were identified from the genomes of C. sinensis var. sinensis (CSS) and C. sinensis var. assamica (CSA), respectively. Phylogenetic analysis revealed that the numbers in CNL groups are predominant in both CSS and CSA. RNA-Seq revealed that the expression of CsNLRs is induced by Colletotrichum fructicola, cold, drought, salt stress and exogenous methyl jasmonate. The 21 CsCNLs that are highly expressed in tea plant under biotic and abiotic stresses as well as during bud dormancy and in different tissues are identified. Gene structure analysis revealed several cis-regulatory elements associated with phytohormones and light responsiveness in the promoter regions of these 21 CsCNLs.

Differential regulation of TNL-mediated immune signaling by redundant helper CNLs

Higher plants utilize nucleotide-binding leucine-rich repeat domain proteins (NLRs) as intracellular immune receptors to recognize pathogen-derived effectors and trigger a robust defense. The Activated Disease Resistance 1 (ADR1) family of coiled-coil NLRs (CNLs) have evolved as helper NLRs that function downstream of many TIR-type sensor NLRs (TNLs). Close homologs of ADR1s form the N REQUIREMENT GENE 1 (NRG1) family in Arabidopsis, the function of which is unclear. Through CRISPR/Cas9 gene editing methods, we discovered that the tandemly repeated NRG1A and NRG1B are functionally redundant and operate downstream of TNLs with differential strengths. Interestingly, ADR1s and NRG1s function in two distinct parallel pathways contributing to TNL-specific immunity. Synergistic effects on basal and TNL-mediated defense were detected among ADR1s and NRG1s. An intact P-loop of NRG1s is not required for mediating signals from sensor TNLs, whereas auto-active NRG1A exhibits autoimmunity. Importantly, NRG1s localize to the cytosol and endomembrane network regardless of the presence of effectors, suggesting a cytosolic activation mechanism. Taken together, different sensor TNLs differentially use two groups of helper NLRs, ADR1s and NRG1s, to transduce downstream defense signals.