Insights from the first genome assembly of Onion (Allium cepa)

Identification of a candidate gene for the I locus determining the dominant white bulb color in onion (Allium cepa L.)

KEY MESSAGE

Through a map-based cloning approach, a gene coding for an R2R3-MYB transcription factor was identified as a causal gene for the I locus controlling the dominant white bulb color in onion. White bulb colors in onion (Allium cepa L.) are determined by either the C or I loci. The causal gene for the C locus was previously isolated, but the gene responsible for the I locus has not been identified yet. To identify candidate genes for the I locus, an approximately 7-Mb genomic DNA region harboring the I locus was obtained from onion and bunching onion (A. fistulosum) whole genome sequences using two tightly linked molecular markers. Within this interval, the AcMYB1 gene, known as a positive regulator of anthocyanin production, was identified. No polymorphic sequences were found between white and red AcMYB1 alleles in the 4,860-bp full-length genomic DNA sequences. However, a 4,838-bp LTR-retrotransposon was identified in the white allele, in the 79-bp upstream coding region from the stop codon. The insertion of this LTR-retrotransposon created a premature stop codon, resulting in the replacement of 26 amino acids with seven different residues. A molecular marker was developed based on the insertion of this LTR-retrotransposon to genotype the I locus. A perfect linkage between bulb color phenotypes and marker genotypes was observed among 5,303 individuals of segregating populations. The transcription of AcMYB1 appeared to be normal in both red and white onions, but the transcription of CHS-A, which encodes chalcone synthase and is involved in the first step of the anthocyanin biosynthesis pathway, was inactivated in the white onions. Taken together, an aberrant AcMYB1 protein produced from the mutant allele might be responsible for the dominant white bulb color in onions.

AlliumDB: a central portal for comparative and functional genomics in Allium

The genus Allium belongs to the botanical family Amaryllidaceae and includes economically important crops such as onion, garlic, bunching onion, and leek, used as vegetables, spices, and traditional medicines. The large sizes of Allium genomes hamper the genetic dissection of agronomically important traits and molecular breeding. With the growing accumulation of genomic, resequencing, transcriptome, and phenotypic data, the demand for an integrative Allium database is increasing. Here we present a user-friendly database, AlliumDB (https://allium.qau.edu.cn), as a functional genomics hub integrating public and in-house data. The database contains all currently available nuclear and organelle genomes for Allium species, with genes comprehensively annotated based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, orthology, gene families, protein families (Pfam), and non-coding RNA families (Rfam). Transcriptome and variation profiles are integrated into dynamic visualization tools. We took phenotypic photographs and generated trait records for hundreds of Allium germplasms collected worldwide, which are included in the database. We incorporated JBrowse for the visualization of gene structures, RNA sequencing data, and variation data. Analysis tools such as the basic local alignment search tool (BLAST), sequence fetch, enrichment, and motif analyses are available to explore potential gene functions. This database incorporates comprehensive Allium genotypic and phenotypic datasets. As the community assembles new genomes and generates resequencing data for Allium germplasms, the database will be improved and continuously updated with these multi-omics data and comparative genomic studies. We expect the AlliumDB database to become a key resource for the study of Allium crops.

A two-sequence motif-based method for the inventory of gene families in fragmented and poorly annotated genome sequences

Databases of genome sequences are growing exponentially, but, in some cases, assembly is incomplete and genes are poorly annotated. For evolutionary studies, it is important to identify all members of a given gene family in a genome. We developed a method for identifying most, if not all, members of a gene family from raw genomes in which assembly is of low quality, using the P-type ATPase superfamily as an example. The method is based on the translation of an entire genome in all six reading frames and the co-occurrence of two family-specific sequence motifs that are in close proximity to each other. To test the method's usability, we first used it to identify P-type ATPase members in the high-quality annotated genome of barley (Hordeum vulgare). Subsequently, after successfully identifying plasma membrane H+-ATPase family members (P3A ATPases) in various plant genomes of varying quality, we tested the hypothesis that the number of P3A ATPases correlates with the ability of the plant to tolerate saline conditions. In 19 genomes of glycophytes and halophytes, the total number of P3A ATPase genes was found to vary from 7 to 22, but no significant difference was found between the two groups. The method successfully identified P-type ATPase family members in raw genomes that are poorly assembled.

Comparative mitochondrial genome analysis reveals a candidate ORF for cytoplasmic male sterility in tropical onion

Cytoplasmic male sterility (CMS) has been widely exploited for hybrid seed production in onions (Allium cepa L.). In contrast to long-day onion cultivars, short-day onion has not yet been investigated for mitochondrial genome structure and DNA rearrangements associated with CMS activity. Here, we report the 3,16,321 bp complete circular mitochondrial genome of tropical onion CMS line (97A). Due to the substantial number of repetitive regions, the assembled mitochondrial genome of maintainer line (97B) remained linear with 15 scaffolds. Additionally, 13 and 20 chloroplast-derived fragments with a size ranging from 143 to 13,984 bp and 153-17,725 bp were identified in the 97A and 97B genomes, respectively. Genome annotation revealed 24 core protein-coding genes along with 24 and 28 tRNA genes in the mitochondrial genomes of 97A and 97B, respectively. Furthermore, comparative genome analysis of the 97A and 97B mitochondrial genomes showed that gene content was almost similar except for the chimeric ORF725 gene which is the extended form of the COX1 gene.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03850-2.

Three-dimensional chromatin architecture in plants - General features and novelties

Research on the three-dimensional (3D) structure of the genome and its distribution within the nuclear space has made a big leap in the last two decades. Work in the animal field has led to significant advances in our general understanding on eukaryotic genome organization. This did not only bring along insights into how the 3D genome interacts with the epigenetic landscape and the transcriptional machinery but also how 3D genome architecture is relevant for fundamental developmental processes, such as cell differentiation. In parallel, the 3D organization of plant genomes have been extensively studied, which resulted in both congruent and novel findings, contributing to a more complete view on how eukaryotic genomes are organized in multiple dimensions. Plant genomes are remarkably diverse in size, composition, and ploidy. Furthermore, as intrinsically sessile organisms without the possibility to relocate to more favorable environments, plants have evolved an elaborate epigenetic repertoire to rapidly respond to environmental challenges. The diversity in genome organization and the complex epigenetic programs make plants ideal study subjects to acquire a better understanding on universal features and inherent constraints of genome organization. Furthermore, considering a wide range of species allows us to study the evolutionary crosstalk between the various levels of genome architecture. In this article, we aim at summarizing important findings on 3D genome architecture obtained in various plant species. These findings cover many aspects of 3D genome organization on a wide range of levels, from gene loops to topologically associated domains and to global 3D chromosome configurations. We present an overview on plant 3D genome organizational features that resemble those in animals and highlight facets that have only been observed in plants to date.

Chromosome-level genomes of three key Allium crops and their trait evolution

Allium crop breeding remains severely hindered due to the lack of high-quality reference genomes. Here we report high-quality chromosome-level genome assemblies for three key Allium crops (Welsh onion, garlic and onion), which are 11.17 Gb, 15.52 Gb and 15.78 Gb in size with the highest recorded contig N50 of 507.27 Mb, 109.82 Mb and 81.66 Mb, respectively. Beyond revealing the genome evolutionary process of Allium species, our pathogen infection experiments and comparative metabolomic and genomic analyses showed that genes encoding enzymes involved in the metabolic pathway of Allium-specific flavor compounds may have evolved from an ancient uncharacterized plant defense system widely existing in many plant lineages but extensively boosted in alliums. Using in situ hybridization and spatial RNA sequencing, we obtained an overview of cell-type categorization and gene expression changes associated with spongy mesophyll cell expansion during onion bulb formation, thus indicating the functional roles of bulb formation genes.

CRISPR/Cas9-mediated editing of PHYTOENE DESATURASE gene in onion (Allium cepa L.)

Introduction

Clustered regularly interspaced short palindromic repeats (CRISPR)/ CRISPR-associated protein 9 (Cas9) is a precise genome editing tool used to introduce genetic modifications in a wide range of crop species. Thus far, there is no report of CRISPR/Cas9-mediated genome editing in onions (Allium cepa L.).

Methods

In the present study, we targeted two exons of the gene coding for Phytoene desaturase (AcPDS) in onion cv. Bhima Super. The sgRNA-carrying constructs were co-cultivated with 8-week-old embryogenic calli using an Agrobacterium-mediated transformation protocol and incubated on the media without hygromycin B selection.

Results and discussion

Out of the total 617 co-cultivated calli, 21 (3.4%) regenerated shoots exhibited three distinct phenotypes: albino, chimeric, and pale green; in comparison to the wild-type non-transformed regenerated shoots. Total chlorophyll content was drastically reduced in albino shoots and significantly decreased in chimeric shoots. Out of the six Cas9 gene PCR-confirmed regenerated shoots, two exhibited the albino phenotype due to insertions/deletions (InDels) and substitution-based mutations in and around the AcPDS target sites. Deep amplicon sequencing revealed a significantly variable InDel frequency between two sgRNAs, ranging from 1.2% to 63.4%, along with a 53.4% substitution frequency. The mutation of the AcPDS gene generated a visually detectable albino phenotype, thus confirming the successful editing of the AcPDS gene. This is the first time a CRISPR/Cas9-mediated genome editing protocol has been successfully established in onion, with the AcPDS gene serving as an example. This study will provide the necessary momentum for researchers to further basic and applied research on onions.

Development of SNP markers linked to purple blotch resistance for marker-assisted selection in onion (Allium cepa L.) breeding

Purple blotch (PB), caused by Alternaria porri (Ellis) Cifferi, is one of the most destructive diseases of onion worldwide. Rapid development and deployment of resistant onion varieties is the most effective approach to control this disease. A single dominant gene, ApR1 was previously linked to PB resistance in onion cultivar 'Arka Kalyan'. In this study, an advanced RIL population derived from a cross between the resistant (Arka Kalyan) and susceptible (Agrifound Rose) cultivar of onion was used to fine map the resistant locus with SNP markers. Twenty plants from the RIL population, ten each with disease resistance and susceptibility trait, were subjected to restriction site-associated DNA sequencing (RAD-Seq) and generated 7388 single nucleotide polymorphisms (SNPs). Correlation analysis between marker genotypes and PB disease phenotype on the 20 plants identified 27 SNPs as candidate markers linked to ApR1 gene for PB resistance. Six candidate SNPs were converted to Kompetitive Allele-Specific PCR (KASP) markers designated as ApRsnip5, ApRsnip8, ApRsnip14, ApRsnip21, ApRsnip23 and ApRsnip25. Marker-trait association based on disease phenotyping and KASP genotyping data on 153 RILs confirmed that all six KASP markers were tightly associated with ApR1 gene within the genetic distance of 1.3 CentiMorgan (cM). ApRsnip14 co-segregated with the ApR1 locus. Further, the six KASP markers were tested on 27 onion lines with different genetic backgrounds. ApRsnip14, ApRsnip21, ApRsnip5 and ApRsnip23 not only showed the correct resistance allele in 3 resistance genotypes, but also clustered together in the remaining 24 susceptible lines. Alternatively, ApRsnip8 and ApRsnip25 exhibited false positives in two onion lines which do not have the R-gene. Overall, our results suggest that ApRsnip14 and ApRsnip23 with their close linkage to ApR1 locus and greater applicability on breeding germplasm are recommended in marker-assisted selection for PB resistance in onion breeding program.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03562-7.

The Control of the Crossover Localization in Allium

Meiotic crossovers/chiasmata are not randomly distributed and strictly controlled. The mechanisms behind crossover (CO) patterning remain largely unknown. In Allium cepa, as in the vast majority of plants and animals, COs predominantly occur in the distal 2/3 of the chromosome arm, while in Allium fistulosum they are strictly localized in the proximal region. We investigated the factors that may contribute to the pattern of COs in A. cepa, A. fistulosum and their F1 diploid (2n = 2x = 8C + 8F) and F1 triploid (2n = 3x = 16F + 8C) hybrids. The genome structure of F1 hybrids was confirmed using genomic in situ hybridization (GISH). The analysis of bivalents in the pollen mother cells (PMCs) of the F1 triploid hybrid showed a significant shift in the localization of COs to the distal and interstitial regions. In F1 diploid hybrid, the COs localization was predominantly the same as that of the A. cepa parent. We found no differences in the assembly and disassembly of ASY1 and ZYP1 in PMCs between A. cepa and A. fistulosum, while F1 diploid hybrid showed a delay in chromosome pairing and a partial absence of synapsis in paired chromosomes. Immunolabeling of MLH1 (class I COs) and MUS81 (class II COs) proteins showed a significant difference in the class I/II CO ratio between A. fistulosum (50%:50%) and A. cepa (73%:27%). The MLH1:MUS81 ratio at the homeologous synapsis of F1 diploid hybrid (70%:30%) was the most similar to that of the A. cepa parent. F1 triploid hybrid at the A. fistulosum homologous synapsis showed a significant increase in MLH1:MUS81 ratio (60%:40%) compared to the A. fistulosum parent. The results suggest possible genetic control of CO localization. Other factors affecting the distribution of COs are discussed.

Turning Garlic into a Modern Crop: State of the Art and Perspectives

Garlic is cultivated worldwide for the value of its bulbs, but its cultivation is challenged by the infertility of commercial cultivars and the accumulation of pathogens over time, which occurs as a consequence of vegetative (clonal) propagation. In this review, we summarize the state of the art of garlic genetics and genomics, highlighting recent developments that will lead to its development as a modern crop, including the restoration of sexual reproduction in some garlic strains. The set of tools available to the breeder currently includes a chromosome-scale assembly of the garlic genome and multiple transcriptome assemblies that are furthering our understanding of the molecular processes underlying important traits like the infertility, the induction of flowering and bulbing, the organoleptic properties and resistance to various pathogens.

Vegetable biology and breeding in the genomics era

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

Construction of a high-resolution genetic map and identification of single nucleotide polymorphism markers relevant to flower stalk height in onion

INTRODUCTION

Onion (Allium cepa L., 2n=16) is an economically and nutritionally important vegetable crop worldwide. Construction of a high-resolution genetic map and map-based gene mining in onion have lagged behind other vegetable crops such as tomato and pepper.

METHODS

In this study, we constructed a high-resolution genetic map of onion using 321 F2 individuals from a cross between two double haploid lines DH-1×DH-17 and employing specific length amplified fragment (SLAF)-seq technology. The genetic map containing 10,584 polymorphic SLAFs with 21,250 single nucleotide polymorphism (SNP) markers and 8 linkage groups was developed for onion, which spanned 928.32 cM, with an average distance of 0.09 cM between adjacent markers.

RESULTS

Using this map, we carried out QTL mapping of Ms locus related to the male-fertile trait and reproduced previous mapping results, which proved that this map was of good quality. Then, four QTLs (located on LG2, LG5, and LG8) were detected for flower stalk height, explaining 26.60% of the phenotypic variance. Among them, we proposed that 20 SLAF markers (in three QTLs) of flower stalk height trait were effective favorable allelic variant markers associated with heterosis.

DISCUSSION

Overall, the genetic map was structured using SLAF-seq based on DH lines, and it is the highest-quality and highest-resolution linkage map of onion to date. It lays a foundation for the fine mapping and candidate gene identification of flower stalk height, and provides new insights into the developmental genetic mechanisms in onion breeding.

Targeted Metabolome Profiling of Indonesian Shallots and Japanese Long-Day/Short-Day Bulb Onions

In this study, targeted metabolome analysis was applied to identify the discriminative metabolites between Indonesian shallot landraces, Japanese long-day onion (LDO) varieties, and Japanese short-day onion (SDO) varieties. In total, 172 metabolite signal intensities were subjected to multivariate PLS-DA, VIP, and random forest modeling to gain further insight into genotype-specific metabolites. PLS-DA divides the examined genotypes into three different clusters, implying that shallot landraces exhibited a distinct metabolite profile compared with Japanese LDO and SDO varieties. The PLS-DA, VIP, and random forest results indicated that the shallot and LDO are richer in metabolite constituents in comparison with the SDO. Specifically, amino acids and organosulfur compounds were the key characteristic metabolites in shallot and LDO genotypes. The analysis of S-alk(en)yl-L-cysteine sulfoxide (ACSO) compounds showed higher accumulation in the shallot landraces relative to LDO and SDO varieties, which explains the stronger pungency and odor in shallots. In addition, the LDO showed higher ACSO content compared with the SDO, implying that long-day cultivation might enhance sulfur assimilation in the Japanese onion. The LDO 'Super Kitamomiji' and the shallots 'Probolinggo' and 'Thailand' showed higher ACSO content than other varieties, making it useful for Allium breeding to improve the flavor and stress tolerance of onions.

Chromosome-level genome assembly of bunching onion illuminates genome evolution and flavor formation in Allium crops

The Allium genus is cultivated globally as vegetables, condiments, or medicinal plants and is characterized by large genomes and strong pungency. However, the genome evolution and genomic basis underlying their unique flavor formation remain poorly understood. Herein, we report an 11.27-Gb chromosome-scale genome assembly for bunching onion (A. fistulosum). The uneven bursts of long-terminal repeats contribute to diversity in genome constituents, and dispersed duplication events largely account for gene expansion in Allium genomes. The extensive duplication and differentiation of alliinase and lachrymatory factor synthase manifest as important evolutionary events during flavor formation in Allium crops. Furthermore, differential selective preference for flavor-related genes likely lead to the variations in isoalliin content in bunching onions. Moreover, we reveal that China is the origin and domestication center for bunching onions. Our findings provide insights into Allium genome evolution, flavor formation and domestication history and enable future genome-assisted breeding of important traits in these crops.

De novo assembly and characterization of the draft genome of the cashew (Anacardium occidentale L.)

Cashew is the second most important tree nut crop in the global market. Cashew is a diploid and heterozygous species closely related to the mango and pistachio. Its improvement by conventional breeding is slow due to the long juvenile phase. Despite the economic importance, very little genomics/transcriptomics information is available for cashew. In this study, the Oxford nanopore reads and Illumina reads were used for de novo assembly of the cashew genome. The hybrid assembly yielded a 356.6 Mb genome corresponding to 85% of the estimated genome size (419 Mb). The BUSCO analysis showed 91.8% of genome completeness. Transcriptome mapping showed 92.75% transcripts aligned with the assembled genome. Gene predictions resulted in the identification of 31,263 genes coding for a total of 35,000 gene isoforms. About 46% (165 Mb) of the cashew genome comprised of repetitive sequences. Phylogenetic analyses of the cashew with nine species showed that it was closely related to Mangifera indica. Analysis of cashew genome revealed 3104 putative R-genes. The first draft assembly of the genome, transcriptome and R gene information generated in this study would be the foundation for understanding the molecular basis of economic traits and genomics-assisted breeding in cashew.

TVIR: a comprehensive vegetable information resource database for comparative and functional genomic studies

Vegetables are an indispensable part of the daily diet of humans. Therefore, it is vital to systematically study the genomic data of vegetables and build a platform for data sharing and analysis. In this study, a comprehensive platform for vegetables with a user-friendly Web interface-The Vegetable Information Resource (TVIR, http://tvir.bio2db.com)-was built based on the genomes of 59 vegetables. TVIR database contains numerous important functional genes, including 5215 auxin genes, 2437 anthocyanin genes, 15 002 flowering genes, 79 830 resistance genes, and 2639 glucosinolate genes of 59 vegetables. In addition, 2597 N6-methyladenosine (m6A) genes were identified, including 513 writers, 1058 erasers, and 1026 readers. A total of 2 101 501 specific clustered regularly interspaced short palindromic repeat (CRISPR) guide sequences and 17 377 miRNAs were detected and deposited in TVIR database. Information on gene synteny, duplication, and orthologs is also provided for 59 vegetable species. TVIR database contains 2 346 850 gene annotations by the Swiss-Prot, TrEMBL, Gene Ontology (GO), Pfam, and Non-redundant (Nr) databases. Synteny, Primer Design, Blast, and JBrowse tools are provided to facilitate users in conducting comparative genomic analyses. This is the first large-scale collection of vegetable genomic data and bioinformatic analysis. All genome and gene sequences, annotations, and bioinformatic results can be easily downloaded from TVIR. Furthermore, transcriptome data of 98 vegetables have been collected and collated, and can be searched by species, tissues, or different growth stages. TVIR is expected to become a key hub for vegetable research globally. The database will be updated with newly assembled vegetable genomes and comparative genomic studies in the future.

Genotyping-by-Sequencing Defines Genetic Structure within the “Acquaviva” Red Onion Landrace

Genetic structure and distinctive features of landraces, such as adaptability to local agro-ecosystems and specific qualitative profiles, can be substantially altered by the massive introduction of allochthonous germplasm. The landrace known as “Cipolla rossa di Acquaviva” (Acquaviva red onion, further referred to as ARO) is traditionally cultivated and propagated in a small area of the Apulia region (southern Italy). However, the recent rise of its market value and cultivation area is possibly causing genetic contamination with foreign propagating material. In this work, genotyping-by-sequencing (GBS) was used to characterize genetic variation of seven onion populations commercialized as ARO, as well as one population of the landrace “Montoro” (M), which is phenotypically similar, but originates from another cultivation area and displays different qualitative features. A panel of 5011 SNP markers was used to perform parametric and non-parametric genetic structure analyses, which supported the hypothesis of genetic contamination of germplasm commercialized as ARO with a gene pool including the M landrace. Four ARO populations formed a core genetic group, homogeneous and clearly distinct from the other ARO and M populations. Conversely, the remaining three ARO populations did not display significant differences with the M population. A set of private alleles for the ARO core genetic group was identified, indicating the possibility to trace the ARO landrace by means of a SNP-based molecular barcode. Overall, the results of this study provide a framework for further breeding activities and the traceability of the ARO landrace.

Integrating Genetic and Chromosome Maps of Allium cepa: From Markers Visualization to Genome Assembly Verification

The ability to directly look into genome sequences has opened great opportunities in plant breeding. Yet, the assembly of full-length chromosomes remains one of the most difficult problems in modern genomics. Genetic maps are commonly used in de novo genome assembly and are constructed on the basis of a statistical analysis of the number of recombinations. This may affect the accuracy of the ordering and orientation of scaffolds within the chromosome, especially in the region of recombination suppression. Moreover, it is impossible to assign contigs lacking DNA markers. Here, we report the use of Tyr-FISH to determine the position of the short DNA sequence of markers and non-mapped unique copy sequence on the physical chromosomes of a large-genome onion (Allium cepa L.). In order to minimize potential background masking of the target signal, we improved our earlier developed pipeline for probe design. A total of 23 markers were located on physical chromosomes 2 and 6. The order of markers was corrected by the integration of genetic, pseudochromosome maps and cytogenetic maps. Additionally, the position of the mlh1 gene, which was not on the genetic map, was defined on physical chromosome 2. Tyr-FISH mapping showed that the order of 23.1% (chromosome 2) and 27.3% (chromosome 6) of the tested genes differed between physical chromosomes and pseudochromosomes. The results can be used for the improvement of pseudochromosome 2 and 6 assembly. The present study aims to demonstrate the value of the in situ visualization of DNA sequences in chromosome-scaffold genome assembly.

Development of a genome-wide marker design workflow for onions and its application in target amplicon sequencing-based genotyping

Onions are one of the most widely cultivated vegetables worldwide; however, the development and utilization of molecular markers have been limited because of the large genome of this plant. We present a genome-wide marker design workflow for onions and its application in a high-throughput genotyping method based on target amplicon sequencing. The efficiency of the method was evaluated by genotyping of F₂ populations. In the marker design workflow, unigene and genomic sequence data sets were constructed, and polymorphisms between parental lines were detected through transcriptome sequence analysis. The positions of polymorphisms detected in the unigenes were mapped onto the genome sequence, and primer sets were designed. In total, 480 markers covering the whole genome were selected. By genotyping an F₂ population, 329 polymorphic sites were obtained from the estimated positions or the flanking sequences. However, missing or sparse marker regions were observed in the resulting genetic linkage map. We modified the markers to cover these regions by genotyping the other F₂ populations. The grouping and order of markers on the linkages were similar across the genetic maps. Our marker design workflow and target amplicon sequencing are useful for genome-wide genotyping of onions owing to their reliability, cost effectiveness, and flexibility.

Functional characterization and vacuolar localization of fructan exohydrolase derived from onion (Allium cepa)

Fructans such as inulin and levan accumulate in certain taxonomic groups of plants and are a reserve carbohydrate alternative to starch. Onion (Allium cepa L.) is a typical plant species that accumulates fructans, and it synthesizes inulin-type and inulin neoseries-type fructans in the bulb. Although genes for fructan biosynthesis in onion have been identified so far, no genes for fructan degradation had been found. In this study, phylogenetic analysis predicted that we isolated a putative vacuolar invertase gene (AcpVI1), but our functional analyses demonstrated that it encoded a fructan 1-exohydrolase (1-FEH) instead. Assessments of recombinant proteins and purified native protein showed that the protein had 1-FEH activity, hydrolyzing the β-(2,1)-fructosyl linkage in inulin-type fructans. Interestingly, AcpVI1 had an amino acid sequence close to those of vacuolar invertases and fructosyltransferases, unlike all other FEHs previously found in plants. We showed that AcpVI1 was localized in the vacuole, as are onion fructosyltransferases Ac1-SST and Ac6G-FFT. These results indicate that fructan-synthesizing and -degrading enzymes are both localized in the vacuole. In contrast to previously reported FEHs, our data suggest that onion 1-FEH evolved from a vacuolar invertase and not from a cell wall invertase. This demonstrates that classic phylogenetic analysis on its own is insufficient to discriminate between invertases and FEHs, highlighting the importance of functional markers in the nearby active site residues.

Cytotoxic activity of Nep1-like proteins on monocots

Necrosis- and ethylene-inducing peptide 1 (Nep1)-like proteins (NLPs) are found throughout several plant-associated microbial taxa and are typically considered to possess cytolytic activity exclusively on dicot plant species. However, cytolytic NLPs are also produced by pathogens of monocot plants such as the onion (Allium cepa) pathogen Botrytis squamosa. We determined the cytotoxic activity of B. squamosa BsNep1, as well as other previously characterized NLPs, on various monocot plant species and assessed the plant plasma membrane components required for NLP sensitivity. Leaf infiltration of NLPs showed that onion cultivars are differentially sensitive to NLPs, and analysis of their sphingolipid content revealed that the GIPC series A : series B ratio did not correlate to NLP sensitivity. A tri-hybrid population derived from a cross between onion and two wild relatives showed variation in NLP sensitivity within the population. We identified a quantitative trait locus (QTL) for NLP insensitivity that colocalized with a previously identified QTL for B. squamosa resistance and the segregating trait of NLP insensitivity correlated with the sphingolipid content. Our results demonstrate the cytotoxic activity of NLPs on several monocot plant species and legitimize their presence in monocot-specific plant pathogens.

Biosynthesis and Metabolism of Garlic Odor Compounds in Cultivated Chinese Chives (Allium tuberosum) and Wild Chinese Chives (Allium hookeri)

Chinese chives is a popular herb vegetable and medicine in Asian countries. Southwest China is one of the centers of origin, and the mountainous areas in this region are rich in wild germplasm. In this study, we collected four samples of germplasm from different altitudes: a land race of cultivated Chinese chives (Allium tuberosum), wide-leaf chives and extra-wide-leaf chives (Allium hookeri), and ovoid-leaf chives (Allium funckiaefolium). Leaf metabolites were detected and compared between A. tuberosum and A. hookeri. A total of 158 differentially accumulated metabolites (DAM) were identified by Gas Chromatography—Mass Spectrometry (GC-MS) and Liquid Chromatography—Mass Spectrometry (LC-MS), among which there was a wide range of garlic odor compounds, free amino acids, and sugars. A. hookeri contains a higher content of fructose, garlic odor compounds, and amino acids than A. tuberosum, which is supported by the higher expression level of biosynthetic genes revealed by transcriptome analysis. A. hookeri accumulates the same garlic odor compound precursors that A. tuberosum does (mainly methiin and alliin). We isolated full-length gene sequences of phytochelatin synthase (PCS), γ-glutamyltranspeptidases (GGT), flavin-containing monooxygenase (FMO), and alliinase (ALN). These sequences showed closer relations in phylogenetic analysis between A. hookeri and A. tuberosum (with sequence identities ranging from 86% to 90%) than with Allium cepa or Allium sativum (which had a lower sequence identity ranging from 76% to 88%). Among these assayed genes, ALN, the critical gene controlling the conversion of odorless precursors into odor compounds, was undetected in leaves, bulbs, and roots of A. tuberosum, which could account for its weaker garlic smell. Moreover, we identified a distinct FMO1 gene in extra-wide-leaf A. hookeri that is due to a CDS-deletion and frameshift mutation. These results above reveal the molecular and metabolomic basis of impressive strong odor in wild Chinese chives.

WRKY Gene Family Drives Dormancy Release in Onion Bulbs

Onion (Allium cepa L.) is an important bulb crop grown worldwide. Dormancy in bulbous plants is an important physiological state mainly regulated by a complex gene network that determines a stop of vegetative growth during unfavorable seasons. Limited knowledge on the molecular mechanisms that regulate dormancy in onion were available until now. Here, a comparison between uninfected and onion yellow dwarf virus (OYDV)-infected onion bulbs highlighted an altered dormancy in the virus-infected plants, causing several symptoms, such as leaf striping, growth reduction, early bulb sprouting and rooting, as well as a lower abscisic acid (ABA) level at the start of dormancy. Furthermore, by comparing three dormancy stages, almost five thousand four hundred (5390) differentially expressed genes (DEGs) were found in uninfected bulbs, while the number of DEGs was significantly reduced (1322) in OYDV-infected bulbs. Genes involved in cell wall modification, proteolysis, and hormone signaling, such as ABA, gibberellins (GAs), indole-3-acetic acid (IAA), and brassinosteroids (BRs), that have already been reported as key dormancy-related pathways, were the most enriched ones in the healthy plants. Interestingly, several transcription factors (TFs) were up-regulated in the uninfected bulbs, among them three genes belonging to the WRKY family, for the first time characterized in onion, were identified during dormancy release. The involvement of specific WRKY genes in breaking dormancy in onion was confirmed by GO enrichment and network analysis, highlighting a correlation between AcWRKY32 and genes driving plant development, cell wall modification, and division via gibberellin and auxin homeostasis, two key processes in dormancy release. Overall, we present, for the first time, a detailed molecular analysis of the dormancy process, a description of the WRKY-TF family in onion, providing a better understanding of the role played by AcWRKY32 in the bulb dormancy release. The TF co-expressed genes may represent targets for controlling the early sprouting in onion, laying the foundations for novel breeding programs to improve shelf life and reduce postharvest.

Using Transcriptome Analysis to Explore Gray Mold Resistance-Related Genes in Onion (Allium cepa L.)

Gray mold disease caused by Botrytis in onions (Allium cepa L.) during growth and storage negatively affects their yield and quality. Exploring the genes related to gray mold resistance in onion and their application to the breeding of resistant onion lines will support effective and ecological control methods of the disease. Here, the genetic relationship of 54 onion lines based on random amplified polymorphic DNA (RAPD) and in vitro-cultured onion lines infected with gray mold were used for screening resistance and susceptibility traits. Two genetically related onion lines were selected, one with a resistant and one with a susceptible phenotype. In vitro gray mold infection was repeated with these two lines, and leaf samples were collected for gene expression studies in time series. Transcript sequences obtained by RNA sequencing were subjected to DEG analysis, variant analysis, and KEGG mapping. Among the KEGG pathways, ‘α-linoleic acid metabolism’ was selected because the comparison of the time series expression pattern of Jasmonate resistant 1 (JAR1), Coronatine-insensitive protein 1 (COI 1), and transcription factor MYC2 (MYC2) genes between the resistant and susceptible lines revealed its significant relationship with gray-mold-resistant phenotypes. Expression pattern and SNP of the selected genes were verified by quantitative real-time PCR and high-resolution melting (HRM) analysis, respectively. The results of this study will be useful for the development of molecular marker and finally breeding of gray-mold-resistant onions.

Reference standards for flow cytometric estimation of absolute nuclear DNA content in plants

The estimation of nuclear DNA content has been by far the most popular application of flow cytometry in plants. Because flow cytometry measures relative fluorescence intensities of nuclei stained by a DNA fluorochrome, ploidy determination, and estimation of the nuclear DNA content in absolute units both require comparison to a reference standard of known DNA content. This implies that the quality of the results obtained depends on the standard selection and use. Internal standardization, when the nuclei of an unknown sample and the reference standard are isolated, stained, and measured simultaneously, is mandatory for precise measurements. As DNA peaks representing G₁/G₀ nuclei of the sample and standard appear on the same histogram of fluorescence intensity, the quotient of their position on the fluorescence intensity axis provides the quotient of DNA amounts. For the estimation of DNA amounts in absolute units, a number of well‐established standards are now available to cover the range of known plant genome sizes. Since there are different standards in use, the standard and the genome size assigned to it has always to be reported. When none of the established standards fits, the introduction of a new standard species is needed. For this purpose, the regression line approach or simultaneous analysis of the candidate standard with several established standards should be prioritized. Moreover, the newly selected standard organism has to fulfill a number of requirements: it should be easy to identify and maintain, taxonomically unambiguous, globally available, with known genome size stability, lacking problematic metabolites, suitable for isolation of sufficient amounts of nuclei, and enabling measurements with low coefficients of variation of DNA peaks, hence suitable for the preparation of high quality samples.