Integrative structural annotation of de novo RNA-Seq provides an accurate reference gene set of the enormous genome of the onion (Allium cepa L.)

Purification and Characterization of α-Mannosidase from Onion, Allium cepa

α-Mannosidase (ALMAN) extracted from onion (Allium cepa) was purified by column chromatography such as hydrophobic and gel filtration. ALMAN is an acidic α-mannosidase that exhibits maximum activity against pNP-α-Man at pH 4.0-5.0 at 50°C. Amino acid sequence analysis of ALMAN was consistent with α-mannosidase deduced from Allium cepa transcriptome analysis. The gene alman was amplified by PCR using mRNA extracted from onions, and a full-length gene of 3,054 bp encoding a protein of 1,018 amino acid residues was revealed. ALMAN is classified as Glycoside Hydrolase Family (GH) 38 and showed homology with other plant-derived α-mannosidases such as tomato and hot pepper.

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.

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.

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.

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.

Genotyping-by-Sequencing Derived Genetic Linkage Map and Quantitative Trait Loci for Sugar Content in Onion (Allium cepa L.)

Onion (2n = 2x = 16) has been a nutritional, medicinal and economically valuable vegetable crop all over the world since ancient times. To accelerate the molecular breeding in onion, genetic linkage maps are prerequisite. However, construction of genetic linkage maps of onion remains relatively rudimentary due to a large genome (about 16.3 Gbp) as well as biennial life cycle, cross-pollinated nature, and high inbreeding depression. In this study, we constructed single nucleotide polymorphism (SNP)-based genetic linkage map of onion in an F₂ segregating population derived from a cross between the doubled haploid line '16P118' and inbred line 'Sweet Green' through genotyping by sequencing (GBS). A total of 207.3 Gbp of raw sequences were generated using an Illumina HiSeq X system, and 24,341 SNPs were identified with the criteria based on three minimum depths, lower than 30% missing rate, and more than 5% minor allele frequency. As a result, an onion genetic linkage map consisting of 216 GBS-based SNPs were constructed comprising eight linkage groups spanning a genetic length of 827.0 cM. Furthermore, we identified the quantitative trait loci (QTLs) for the sucrose, glucose, fructose, and total sugar content across the onion genome. We identified a total of four QTLs associated with sucrose (qSC4.1), glucose (qGC5.1), fructose (qFC5.1), and total sugar content (qTSC5.1) explaining the phenotypic variation (R²%) ranging from 6.07-11.47%. This map and QTL information will contribute to develop the molecular markers to breed the cultivars with high sugar content in onion.

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

Onion is an important vegetable crop with an estimated genome size of 16 Gb. We describe the de novo assembly and ab initio annotation of the genome of a doubled haploid onion line DHCU066619, which resulted in a final assembly of 14.9 Gb with an N50 of 464 Kb. Of this, 2.4 Gb was ordered into eight pseudomolecules using four genetic linkage maps. The remainder of the genome is available in 89.6 K scaffolds. Only 72.4% of the genome could be identified as repetitive sequences and consist, to a large extent, of (retro) transposons. In addition, an estimated 20% of the putative (retro) transposons had accumulated a large number of mutations, hampering their identification, but facilitating their assembly. These elements are probably already quite old. The ab initio gene prediction indicated 540,925 putative gene models, which is far more than expected, possibly due to the presence of pseudogenes. Of these models, 47,066 showed RNASeq support. No gene rich regions were found, genes are uniformly distributed over the genome. Analysis of synteny with Allium sativum (garlic) showed collinearity but also major rearrangements between both species. This assembly is the first high-quality genome sequence available for the study of onion and will be a valuable resource for further research.

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

Onion is an important vegetable crop with an estimated genome size of 16 Gb. We describe the de novo assembly and ab initio annotation of the genome of a doubled haploid onion line DHCU066619, which resulted in a final assembly of 14.9 Gb with an N50 of 464 Kb. Of this, 2.4 Gb was ordered into eight pseudomolecules using four genetic linkage maps. The remainder of the genome is available in 89.6 K scaffolds. Only 72.4% of the genome could be identified as repetitive sequences and consist, to a large extent, of (retro) transposons. In addition, an estimated 20% of the putative (retro) transposons had accumulated a large number of mutations, hampering their identification, but facilitating their assembly. These elements are probably already quite old. The ab initio gene prediction indicated 540,925 putative gene models, which is far more than expected, possibly due to the presence of pseudogenes. Of these models, 47,066 showed RNASeq support. No gene rich regions were found, genes are uniformly distributed over the genome. Analysis of synteny with Allium sativum (garlic) showed collinearity but also major rearrangements between both species. This assembly is the first high-quality genome sequence available for the study of onion and will be a valuable resource for further research.

Determination of reliable reference genes for gene expression studies in Chinese chive (Allium tuberosum) based on the transcriptome profiling

Chinese chive (Allium tuberosum) is widely cultivated around the world for its unique flavor, nutrient, and medicinal values, yet its molecular mechanism on flavor formation and other metabolic pathways remains intangible. The elucidation of these complex processes begins with investigating the expression of the genes of interest, however the appropriate reference genes (RGs) for normalizing the gene expression are still unavailable in A. tuberosum. To fill this lacuna, transcriptome-wide screening was undertaken to identify the most stable genes according to the analysis of their FPKM values. The expression stability of the RGs was further evaluated using geNorm, NormFinder, BestKeeper, and RefFinder algorithms. The comprehensive analysis showed that GLY1 and SKP1, instead of two traditionally used RGs (eIF1α and ACT2), were the most stable genes across diverse A. tuberosum tissues, indicating the necessity to carefully validate the stability of RGs prior to their use for normalizations. As indicated by geNorm, the normalizations with at least two RGs could give more accurate results. qRT-PCR experiments were conducted with randomly selected genes, demonstrating that normalization with a combination of GLY1 and SKP1 resulted in reliable normalization results. Our finding represents the first attempt toward establishing a standardized qRT-PCR analysis in this economically important vegetable.

Systematic Identification and Validation of Housekeeping and Tissue-Specific Genes in Allotetraploid Chenopodium quinoa

Quinoa is a gluten-free food crop that contains all the essential amino acids and vitamins. The selection of proper housekeeping and tissue-specific genes is the crucial prerequisite for gene expression analysis using the common approach, real-time quantitative PCR (RT-qPCR). In this study, we identified 40 novel candidate housekeeping genes by the minimum transcript per million (TPM), coefficient of variation (CV) and maximum fold change (MFC) methods and 19 candidate tissue-specific genes by the co-expression network method based on an RNA-seq dataset that included 53 stem, leaf, flower and seed samples, as well as additional shoot and root samples under different stresses. The expression stability of 12 housekeeping and tissue-specific genes, as well as that of another two traditionally used housekeeping genes, was further evaluated using qPCR and ranked using NormFinder, BestKeeper and the comparative delta-Ct method. The results demonstrated that MIF, RGGA, VATE and UBA2B were ranked as the top four most stable candidate housekeeping genes. qPCR analysis also revealed three leaf-specific genes and five root-specific genes, but no stem-specific gene was identified. Gene Ontology (GO) enrichment analysis identified that housekeeping genes were mainly enriched in the small molecule metabolic process, organonitrogen compound metabolic process, NAD binding and ligase activity. In addition, tissue-specific genes are closely associated with the major functions of a specific tissue. Specifically, GO terms “photosynthesis” and “thylakoid” were most significantly overrepresented in candidate leaf-specific genes. The novel housekeeping and tissue-specific genes in our study will enable better normalization and quantification of transcript levels in quinoa.

Comparative analysis of de novo genomes reveals dynamic intra-species divergence of NLRs in pepper

BACKGROUND

Peppers (Capsicum annuum L.) containing distinct capsaicinoids are the most widely cultivated spices in the world. However, extreme genomic diversity among species represents an obstacle to breeding pepper.

RESULTS

Here, we report de novo genome assemblies of Capsicum annuum 'Early Calwonder (non-pungent, ECW)' and 'Small Fruit (pungent, SF)' along with their annotations. In total, we assembled 2.9 Gb of ECW and SF genome sequences, representing over 91% of the estimated genome sizes. Structural and functional annotation of the two pepper genomes generated about 35,000 protein-coding genes each, of which 93% were assigned putative functions. Comparison between newly and publicly available pepper gene annotations revealed both shared and specific gene content. In addition, a comprehensive analysis of nucleotide-binding and leucine-rich repeat (NLR) genes through whole-genome alignment identified five significant regions of NLR copy number variation (CNV). Detailed comparisons of those regions revealed that these CNVs were generated by intra-specific genomic variations that accelerated diversification of NLRs among peppers.

CONCLUSIONS

Our analyses unveil an evolutionary mechanism responsible for generating CNVs of NLRs among pepper accessions, and provide novel genomic resources for functional genomics and molecular breeding of disease resistance in Capsicum species.

Comprehensive Transcriptome and Metabolic Profiling of Petal Color Development in Lycoris sprengeri

Lycoris sprengeri (L. sprengeri) is an important ornamental bulbous plant, and its numerous varieties in different color forms are widely planted. Multiple color types of petals in L. sprengeri provide us with possibilities to delineate the complicated metabolic networks underlying the biochemical traits behind color formation in this plant species, especially petal color. In this study, we sequenced and annotated a reference transcriptome of pink and white petals of L. sprengeri and analyzed the metabolic role of anthocyanin biosynthesis in regulating color pigment metabolism. Briefly, white and pink petal samples were sequenced with an Illumina platform, to obtain the reads that could be assembled into 100,778 unique sequences. Sequences expressed differentially between white vs. pink petals were further annotated with the terms of Gene Ontology (GO), Clusters of Orthologous Groups (COG), Kyoto Encyclopedia of Genes and Genomes (KEGG), and eggNOG. Gene expression analyses revealed the repression of anthocyanin and steroid biosynthesis enzymes and R2R3 MYB transcription factor (TF) genes in white petals compared to pink petals. Furthermore, the targeted metabolic profiling of anthocyanins revealed that color-related delphinidin (Del) and cyanidin (Cy) pigments are lower in white petals, which correlate well with the reduced gene expression levels of anthocyanin biosynthesis genes. Taken together, it is hypothesized that anthocyanin biosynthesis, steroid biosynthesis, and R2R3 MYB TFs may play vital regulatory roles in petal color development in L. sprengeri. This work provides a valuable genomic resource for flower breeding and metabolic engineering in horticulture and markers for studying the flower trait evolution of L. sprengeri.

Omics approaches in Allium research: Progress and way ahead

BACKGROUND

The genus Allium (Family: Amaryllidaceae) is an economically important group of crops cultivated worldwide for their use as a vegetable and spices. Alliums are also well known for their nutraceutical properties. Among alliums, onion, garlic, leek, and chives cultivated worldwide. Despite their substantial economic and medicinal importance, the genome sequence of any of the Allium is not available, probably due to their large genome sizes. Recently evolved omics technologies are highly efficient and robust in elucidating molecular mechanisms of several complex life processes in plants. Omics technologies, such as genomics, transcriptomics, proteomics, metabolomics, metagenomics, etc. have the potential to open new avenues in research and improvement of allium crops where genome sequence information is limited. A significant amount of data has been generated using these technologies for various Allium species; it will help in understanding the key traits in Allium crops such as flowering, bulb development, flavonoid biosynthesis, male sterility and stress tolerance at molecular and metabolite level. This information will ultimately assist us in speeding up the breeding in Allium crops.

METHOD

In the present review, major omics approaches, and their progress, as well as potential applications in Allium crops, could be discussed in detail.

RESULTS

Here, we have discussed the recent progress made in Allium research using omics technologies such as genomics, transcriptomics, micro RNAs, proteomics, metabolomics, and metagenomics. These omics interventions have been used in alliums for marker discovery, the study of the biotic and abiotic stress response, male sterility, organ development, flavonoid and bulb color, micro RNA discovery, and microbiome associated with Allium crops. Further, we also emphasized the integrated use of these omics platforms for a better understanding of the complex molecular mechanisms to speed up the breeding programs for better cultivars.

CONCLUSION

All the information and literature provided in the present review throws light on the progress and potential of omics platforms in the research of Allium crops. We also mentioned a few research areas in Allium crops that need to be explored using omics technologies to get more insight. Overall, alliums are an under-studied group of plants, and thus, there is tremendous scope and need for research in Allium species.

TGFam-Finder: a novel solution for target-gene family annotation in plants

Whole-genome annotation error that omits essential protein-coding genes hinders further research. We developed Target Gene Family Finder (TGFam-Finder), an alternative tool for the structural annotation of protein-coding genes containing target domain(s) of interest in plant genomes. TGFam-Finder took considerably reduced annotation run-time and improved accuracy compared to conventional annotation tools. Large-scale re-annotation of 50 plant genomes identified an average of 150, 166 and 86 additional far-red-impaired response 1, nucleotide-binding and leucine-rich-repeat, and cytochrome P450 genes, respectively, that were missed in previous annotations. We detected significantly higher number of translated genes in the new annotations using mass spectrometry data from seven plant species compared to previous annotations. TGFam-Finder along with the new gene models can provide an optimized platform for comprehensive functional, comparative, and evolutionary studies in plants.

Construction of an Onion (Allium cepa L.) Genetic Linkage Map using Genotyping-by-Sequencing Analysis with a Reference Gene Set and Identification of QTLs Controlling Anthocyanin Synthesis and Content

Anthocyanins, the pigmented flavonoids responsible for red and blue colors in horticultural products, promote human health by preventing cancers and lowering the risk of cardiovascular disease. Red onions contain several cyanidin- and peonidin-based anthocyanins. In this study, we constructed a single-nucleotide polymorphism (SNP)-based genetic linkage map in an F₂ segregating population derived from a cross between the inbred line ‘SP3B’ (yellow bulb) and the doubled haploid line ‘H6′ (red bulb) to identify quantitative trait loci (QTLs) for total anthocyanin content of onion bulbs using a genotyping-by-sequencing (GBS) analysis based on a reference gene set. A total of 101.9 Gbp of raw sequences were generated using an Illumina HiSeq 2500 system and a total of 1625 SNP loci were identified with the criteria of three minimum depths, lower than 30% missing rate, and more than 5% minor allele frequency. As a result, an onion genetic linkage map consisting of 319 GBS-based SNP loci and 34 high-resolution melting (HRM) markers was constructed with eight linkage groups and a total genetic distance of 881.4 cM. In addition, the linkage groups were assigned to corresponding chromosomes by comparison with the reference genetic map OH1×5225 through marker development based on common transcripts. The analysis revealed one major QTL, qAS7.1, for anthocyanin synthesis and two significant QTLs, qAC4.1 and qAC4.2, for anthocyanin content. The QTL qAS7.1, located on chromosome 7 with a phenotypic variation of 87.61%, may be a dihydroflavonol 4-reductase (DFR) gene that determines whether the bulb color is red or yellow. The QTLs qAC4.1 and qAC4.2 are separately positioned on chromosome 4 with R² values of 19.43% and 26.28%, respectively. This map and QTL information will contribute to marker development and breeding for high anthocyanin content in bulb onion.

Functional Allium fistulosum Centromeres Comprise Arrays of a Long Satellite Repeat, Insertions of Retrotransposons and Chloroplast DNA

The centromere is a unique part of the chromosome combining a conserved function with an extreme variability in its DNA sequence. Most of our knowledge about the functional centromere organization is obtained from species with small and medium genome/chromosome sizes while the progress in plants with big genomes and large chromosomes is lagging behind. Here, we studied the genomic organization of the functional centromere in Allium fistulosum and A. cepa, both species with a large genome (13 Gb and 16 Gb/1C, 2n = 2x = 16) and large-sized chromosomes. Using low-depth DNA sequencing for these two species and previously obtained CENH3 immunoprecipitation data we identified two long (1.2 Kb) and high-copy repeats, AfCen1K and AcCen1K. FISH experiments showed that AfCen1K is located in all centromeres of A. fistulosum chromosomes while no AcCen1K FISH signals were identified on A. cepa chromosomes. Our molecular cytogenetic and bioinformatics survey demonstrated that these repeats are partially similar but differ in chromosomal location, sequence structure and genomic organization. In addition, we could conclude that the repeats are transcribed and their RNAs are not polyadenylated. We also observed that these repeats are associated with insertions of retrotransposons and plastidic DNA and the landscape of A. cepa and A. fistulosum centromeric regions possess insertions of plastidic DNA. Finally, we carried out detailed comparative satellitome analysis of A. cepa and A. fistulosum genomes and identified a new chromosome- and A. cepa-specific tandem repeat, TR2CL137, located in the centromeric region. Our results shed light on the Allium centromere organization and provide unique data for future application in Allium genome annotation.

Transposition of a non-autonomous DNA transposon in the gene coding for a bHLH transcription factor results in a white bulb color of onions (Allium cepa L.)

A DNA transposon was found in the gene encoding a bHLH transcription factor. Genotypes of the marker tagging this DNA transposon perfectly co-segregated with color phenotypes in large F_2:3 populations A combined approach of bulked segregant analysis and RNA-Seq was used to isolate causal gene for C locus controlling white bulb color in onions (Allium cepa L.). A total of 114 contigs containing homozygous single nucleotide polymorphisms (SNPs) between white and yellow bulked RNAs were identified. Four of them showed high homologies with loci clustered in the middle of chromosome 5. SNPs in 34 contigs were confirmed by sequencing of PCR products. One of these contigs showed perfect linkage to the C locus in F_2:3 populations consisting of 2491 individuals. However, genotypes of molecular marker tagging this contig were inconsistent with color phenotypes of diverse breeding lines. A total of 146 contigs showed differential expression between yellow and white bulks. Among them, transcription levels of B2 gene encoding a bHLH transcription factor were significantly reduced in white RNA bulk and F_2:3 individuals, although there was no SNP in the coding region. Phylogenetic analysis showed that onion B2 was orthologous to bHLH-coding genes regulating anthocyanin biosynthesis pathway in other plant species. Promoter regions of B2 gene were obtained by genome walking and a 577-bp non-autonomous DNA transposon designated as AcWHITE was found in the white allele. Molecular marker tagging AcWHITE showed perfect linkage with the C locus. Marker genotypes of the white allele were detected in some white accessions. However, none of tested red or yellow onions contained AcWHITE insertion, implying that B2 gene was likely to be a casual gene for the C locus.

Transcriptome Analysis and Metabolic Profiling of Lycoris Radiata

Lycoris radiata belongs to the Amaryllidaceae family and is a bulbous plant native to South Korea, China, and Japan. Galantamine, a representative alkaloid of Amaryllidaceae plants, including L. radiata, exhibits selective and dominant acetylcholinesterase inhibition. In spite of the economic and officinal importance of L. radiata, the molecular biological and biochemical information on L. radiata is relatively deficient. Therefore, this study provides functional information of L. radiata, describe galantamine biosynthesis in the various organs, and provide transcriptomic and metabolic datasets to support elucidation of galantamine biosynthesis pathway in future studies. The results of studies conducted in duplicate revealed the presence of a total of 325,609 and 404,019 unigenes, acquired from 9,913,869,968 and 10,162,653,038 raw reads, respectively, after trimming the raw reads using CutAdapt, assembly using Trinity package, and clustering using CD-Hit-EST. All of the assembled unigenes were aligned to the public databases, including National Center for Biotechnology Information (NCBI) non-redundant protein (NR) and nucleotide (Nt) database, SWISS-PROT (UniProt) protein sequence data bank, The Arabidopsis Information Resource (TAIR), the Swiss-Prot protein database, Gene Ontology (GO), and Clusters of Orthologous Groups (COG) database to predict potential genes and provide their functional information. Based on our transcriptome data and published literatures, eight full-length cDNA clones encoding LrPAL2, LrPAL3, LrC4H2, LrC3H, LrTYDC2, LrNNR, LrN4OMT, and LrCYP96T genes, involved in galantamine biosynthesis, were identified in L. radiata. In order to investigate galantamine biosynthesis in different plant parts of L. radiata grown in a growth chamber, gene expression levels were measured through quantitative real-time polymerase chain reaction (qRT-PCR) analysis using these identified genes and galantamine levels were quantified by high-performance liquid chromatography (HPLC) analysis. The qRT-PCR data revealed high expression levels of LrNNR, LrN4OMT, and LrCYP96T in the bulbs, and, as expected, we observed higher amounts of galantamine in the bulbs than in the root and leaves. Additionally, a total of 40 hydrophilic metabolites were detected in the different organs using gas-chromatography coupled with time-of-flight mass spectrometry. In particular, a strong positive correlation between galantamine and sucrose, which provides energy for the secondary metabolite biosynthesis, was observed.

Comparative Tyramide-FISH mapping of the genes controlling flavor and bulb color in Allium species revealed an altered gene order

Evolutionarily related species often share a common order of genes along homeologous chromosomes. Here we report the collinearity disruption of genes located on homeologous chromosome 4 in Allium species. Ultra-sensitive fluorescence in situ hybridization with tyramide signal amplification (tyr-FISH) allowed the visualization of the alliinase multigene family, chalcon synthase gene and EST markers on Allium cepa and Allium fistulosum chromosomes. In A. cepa, bulb alliinase, root alliinase (ALL1) and chalcon synthase (CHS-B) genes were located in the long arm but EST markers (API18 and ACM082) were located in the short arm. In A. fistulosum, all the visualized genes and markers were located in the short arm. Moreover, root alliinase genes (ALL1 and AOB249) showed contrast patterns in number of loci. We suppose that the altered order of the genes/markers is the result of a large pericentric inversion. To get insight into the evolution of the chromosome rearrangement, we mapped the bulb alliinase gene in phylogenetically close and distant species. In the taxonomic clade including A. fistulosum, A. altaicum, A. oschaninii and A. pskemense and in phylogenetically distant species A. roylei and A. nutans, the bulb alliinase gene was located on the short arm of chromosome 4 while, in A. cepa and A. schoenoprasum, the bulb alliinase gene was located on the long arm of chromosome 4. These results have encouraging implications for the further tracing of inverted regions in meiosis of interspecific hybrids and studding chromosome evolution. Also, this finding may have a practical benefit as closely related species are actively used for improving onion crop stock.

Interplay between Jasmonic Acid, Phosphate Signaling and the Regulation of Glycerolipid Homeostasis in Arabidopsis

Jasmonic acid (JA) biosynthesis and signaling are activated in Arabidopsis cultivated in phosphate (Pi) deprived conditions. This activation occurs mainly in photosynthetic tissues and is less important in roots. In leaves, the enhanced biosynthesis of JA coincides with membrane glycerolipid remodeling triggered by the lack of Pi. We addressed the possible role of JA on the dynamics and magnitude of glycerolipid remodeling in response to Pi deprivation and resupply. Based on combined analyses of gene expression, JA biosynthesis and glycerolipid remodeling in wild-type Arabidopsis and in the coi1-16 mutant, JA signaling seems important in the determination of the basal levels of phosphatidylcholine, phosphatidic acid (PA), monogalactosyldiacylglycerol (MGDG) and digalactosyldiacylglycerol. JA impact on MGDG steady state level and fluctuations seem contradictory. In the coi1-16 mutant, the steady state level of MGDG is higher, possibly due to a higher level of PA in the mutant, activating MGD1, and to an increased expression of MGD3. These results support a possible impact of JA in limiting the overall content of this lipid. Concerning lipid variations, upon Pi deprivation, JA seems rather associated with a specific MGDG increase. Following Pi resupply, whereas the expression of glycerolipid remodeling genes returns to basal level, JA biosynthesis and signaling genes are still upregulated, likely due to a JA-induced positive feedback remaining active. Distinct impacts on enzymes synthesizing MGDG, that is, downregulating MGD3, possibly activating MGD1 expression and limiting the activation of MGD1 via PA, might allow JA playing a role in a sophisticated fine tuning of galactolipid variations.

The Power of Genomic in situ Hybridization (GISH) in Interspecific Breeding of Bulb Onion (Allium cepa L.) Resistant to Downy Mildew (Peronospora destructor [Berk.] Casp.)

We exploited the advantages of genomic in situ hybridization (GISH) to monitor the introgression process at the chromosome level using a simple and robust molecular marker in the interspecific breeding of bulb onion (Allium cepa L.) that is resistant to downy mildew. Downy mildew (Peronospora destructor [Berk.] Casp.) is the most destructive fungal disease for bulb onions. With the application of genomic in situ hybridization (GISH) and previously developed DMR1 marker, homozygous introgression lines that are resistant to downy mildew were successfully produced in a rather short breeding time. Considering that the bulb onion is a biennial plant, it took seven years from the F₁ hybrid production to the creation of S₂BC2 homozygous lines that are resistant to downy mildew. Using GISH, it was shown that three progeny plants of S₂BC₂ possessed an A. roylei homozygous fragment in the distal region of the long arm of chromosomes 3 in an A. cepa genetic background. Previously, it was hypothesized that a lethal gene(s) was linked to the downy mildew resistance gene. With the molecular cytogenetic approach, we physically mapped more precisely the lethal gene(s) using the homozygous introgression lines that differed in the size of the A. roylei fragments on chromosome 3.

SNP discovery of Korean short day onion inbred lines using double digest restriction site-associated DNA sequencing

Onion (Allium cepa L.) is an economically important vegetable crop around the world. Genetic and genomic research into various onion accessions will provide insights into the onion genome to enhance breeding strategies and improve crops. However, the onion’s large genome size means that studies of molecular markers are limited in onion. This study aimed to discover high quality single nucleotide polymorphisms (SNPs) from 192 onion inbred lines relating to short-day cultivation in Korea. Paired-end (PE) double digested restriction site-associated DNA sequencing (ddRAD-seq) was used to discover SNPs in onion. A total of 538,973,706 reads (25.9 GB), with an average of 2,658,491 high-quality reads, were generated using ddRAD-seq. With stringent filtering, 1904 SNPs were discovered based on onion reference scaffolds. Further, population structure and genetic relationship studies suggested that two well-differentiated sub-populations exist in onion lines. SNP-associated flanking sequences were also compared with a public non-redundant database for gene ontology and pathway analysis. To our knowledge, this is the first report to identify high-quality SNPs in onion based on reference sequences using the ddRAD-seq platform. The SNP markers identified will be useful for breeders and the research community to deepen their understanding, enhance breeding programs, and support the management of onion genomic resources.

New reference genome sequences of hot pepper reveal the massive evolution of plant disease-resistance genes by retroduplication

Background

Transposable elements are major evolutionary forces which can cause new genome structure and species diversification. The role of transposable elements in the expansion of nucleotide-binding and leucine-rich-repeat proteins (NLRs), the major disease-resistance gene families, has been unexplored in plants.

Results

We report two high-quality de novo genomes (Capsicum baccatum and C. chinense) and an improved reference genome (C. annuum) for peppers. Dynamic genome rearrangements involving translocations among chromosomes 3, 5, and 9 were detected in comparison between C. baccatum and the two other peppers. The amplification of athila LTR-retrotransposons, members of the gypsy superfamily, led to genome expansion in C. baccatum. In-depth genome-wide comparison of genes and repeats unveiled that the copy numbers of NLRs were greatly increased by LTR-retrotransposon-mediated retroduplication. Moreover, retroduplicated NLRs are abundant across the angiosperms and, in most cases, are lineage-specific.

Conclusions

Our study reveals that retroduplication has played key roles for the massive emergence of NLR genes including functional disease-resistance genes in pepper plants.

Electronic supplementary material

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

New reference genome sequences of hot pepper reveal the massive evolution of plant disease-resistance genes by retroduplication

BACKGROUND

Transposable elements are major evolutionary forces which can cause new genome structure and species diversification. The role of transposable elements in the expansion of nucleotide-binding and leucine-rich-repeat proteins (NLRs), the major disease-resistance gene families, has been unexplored in plants.

RESULTS

We report two high-quality de novo genomes (Capsicum baccatum and C. chinense) and an improved reference genome (C. annuum) for peppers. Dynamic genome rearrangements involving translocations among chromosomes 3, 5, and 9 were detected in comparison between C. baccatum and the two other peppers. The amplification of athila LTR-retrotransposons, members of the gypsy superfamily, led to genome expansion in C. baccatum. In-depth genome-wide comparison of genes and repeats unveiled that the copy numbers of NLRs were greatly increased by LTR-retrotransposon-mediated retroduplication. Moreover, retroduplicated NLRs are abundant across the angiosperms and, in most cases, are lineage-specific.

CONCLUSIONS

Our study reveals that retroduplication has played key roles for the massive emergence of NLR genes including functional disease-resistance genes in pepper plants.

Extensive transcriptomic and epigenomic remodelling occurs during Arabidopsis thaliana germination

BACKGROUND

Seed germination involves progression from complete metabolic dormancy to a highly active, growing seedling. Many factors regulate germination and these interact extensively, forming a complex network of inputs that control the seed-to-seedling transition. Our understanding of the direct regulation of gene expression and the dynamic changes in the epigenome and small RNAs during germination is limited. The interactions between genome, transcriptome and epigenome must be revealed in order to identify the regulatory mechanisms that control seed germination.

RESULTS

We present an integrated analysis of high-resolution RNA sequencing, small RNA sequencing and MethylC sequencing over ten developmental time points in Arabidopsis thaliana seeds, finding extensive transcriptomic and epigenomic transformations associated with seed germination. We identify previously unannotated loci from which messenger RNAs are expressed transiently during germination and find widespread alternative splicing and divergent isoform abundance of genes involved in RNA processing and splicing. We generate the first dynamic transcription factor network model of germination, identifying known and novel regulatory factors. Expression of both microRNA and short interfering RNA loci changes significantly during germination, particularly between the seed and the post-germinative seedling. These are associated with changes in gene expression and large-scale demethylation observed towards the end of germination, as the epigenome transitions from an embryo-like to a vegetative seedling state.

CONCLUSIONS

This study reveals the complex dynamics and interactions of the transcriptome and epigenome during seed germination, including the extensive remodelling of the seed DNA methylome from an embryo-like to vegetative-like state during the seed-to-seedling transition. Data are available for exploration in a user-friendly browser at https://jbrowse.latrobe.edu.au/germination_epigenome .

Selection and validation of reference genes for quantitative real-time PCR analysis under different experimental conditions in the leafminer Liriomyza trifolii (Diptera: Agromyzidae)

Liriomyza trifolii is a highly-invasive leafmining insect that causes significant damage to vegetables and horticultural crops worldwide. Relatively few studies have quantified gene expression in L. trifolii using real-time quantitative PCR (RT-qPCR), which is a reliable and sensitive technique for measuring gene expression. RT-qPCR requires the selection of reference genes to normalize gene expression data and control for internal differences between samples. In this study, nine housekeeping genes from L. trifolii were selected for their suitability in normalizing gene expression using geNorm, Normfinder, BestKeeper, the ΔCt method and RefFinder. HSP21.7, which encodes heat shock protein 21.7, was used as a target gene to validate the expression of candidate reference genes. Results indicated that ACTIN and 18S were optimal for developmental stage and low temperature, TUB and 18S showed the most stable expression for sex, and GAPDH and ACTIN were the best reference genes for monitoring gene expression at high temperature. Selection and validation of appropriate reference genes are critical steps in normalizing gene expression levels, which improve the accuracy and quality of expression data. Results of this study provide vital information on reference genes and is valuable in developing a standardized RT-qPCR protocol for functional genomics research in L. trifolii.

Transcriptome analysis of the oriental melon (Cucumis melo L. var. makuwa) during fruit development

Background

The oriental melon (Cucumis melo L. var. makuwa) is one of the most important cultivated cucurbits grown widely in Korea, Japan, and northern China. It is cultivated because its fruit has a sweet aromatic flavor and is rich in soluble sugars, organic acids, minerals, and vitamins. In order to elucidate the genetic and molecular basis of the developmental changes that determine size, color, and sugar contents of the fruit, we performed de novo transcriptome sequencing to analyze the genes expressed during fruit development.

Results

We identified a total of 47,666 of representative loci from 100,875 transcripts and functionally annotated 33,963 of the loci based on orthologs in Arabidopsis thaliana. Among those loci, we identified 5,173 differentially expressed genes, which were classified into 14 clusters base on the modulation of their expression patterns. The expression patterns suggested that the differentially expressed genes were related to fruit development and maturation through diverse metabolic pathways. Analyses based on gene set enrichment and the pathways described in the Kyoto Encyclopedia of Genes and Genomes suggested that the expression of genes involved in starch and sucrose metabolism and carotenoid biosynthesis were regulated dynamically during fruit development and subsequent maturation.

Conclusion

Our results provide the gene expression patterns related to different stages of fruit development and maturation in the oriental melon. The expression patterns give clues about important regulatory mechanisms, especially those involving starch, sugar, and carotenoid biosynthesis, in the development of the oriental melon fruit.

Doubled Haploid 'CUDH2107' as a Reference for Bulb Onion (Allium cepa L.) Research: Development of a Transcriptome Catalogue and Identification of Transcripts Associated with Male Fertility

Researchers working on model plants have derived great benefit from developing genomic and genetic resources using ‘reference’ genotypes. Onion has a large and highly heterozygous genome making the sharing of germplasm and analysis of sequencing data complicated. To simplify the discovery and analysis of genes underlying important onion traits, we are promoting the use of the homozygous double haploid line ‘CUDH2107’ by the onion research community. In the present investigation, we performed transcriptome sequencing on vegetative and reproductive tissues of CUDH2107 to develop a multi-organ reference transcriptome catalogue. A total of 396 million 100 base pair paired reads was assembled using the Trinity pipeline, resulting in 271,665 transcript contigs. This dataset was analysed for gene ontology and transcripts were classified on the basis of putative biological processes, molecular function and cellular localization. Significant differences were observed in transcript expression profiles between different tissues. To demonstrate the utility of our CUDH2107 transcriptome catalogue for understanding the genetic and molecular basis of various traits, we identified orthologues of rice genes involved in male fertility and flower development. These genes provide an excellent starting point for studying the molecular regulation, and the engineering of reproductive traits.

Completion of the mitochondrial genome sequence of onion (Allium cepa L.) containing the CMS-S male-sterile cytoplasm and identification of an independent event of the ccmF N gene split

Cytoplasmic male-sterility (CMS) conferred by the CMS-S cytoplasm has been most commonly used for onion (Allium cepa L.) F₁ hybrid seed production. We first report the complete mitochondrial genome sequence containing CMS-S cytoplasm in this study. Initially, seven contigs were de novo assembled from 150-bp paired-end raw reads produced from the total genomic DNA using the Illumina NextSeq500 platform. These contigs were connected into a single circular genome consisting of 316,363 bp (GenBank accession: KU318712) by PCR amplification. Although all 24 core protein-coding genes were present, no ribosomal protein-coding genes, except rps12, were identified in the onion mitochondrial genome. Unusual trans-splicing of the cox2 gene was verified, and the cox1 gene was identified as part of the chimeric orf725 gene, which is a candidate gene responsible for inducing CMS. In addition to orf725, two small chimeric genes were identified, but no transcripts were detected for these two open reading frames. Thirteen chloroplast-derived sequences, with sizes of 126-13,986 bp, were identified in the intergenic regions. Almost 10 % of the onion mitochondrial genome was composed of repeat sequences. The vast majority of repeats were short repeats of <100 base pairs. Interestingly, the gene encoding ccmF_N was split into two genes. The ccmF _N gene split is first identified outside the Brassicaceae family. The breakpoint in the onion ccmF _N gene was different from that of other Brassicaceae species. This split of the ccmF _N gene was also present in 30 other Allium species. The complete onion mitochondrial genome sequence reported in this study would be fundamental information for elucidation of onion CMS evolution.

An RNA-Seq-based reference transcriptome for Citrus

Previous RNA-Seq studies in citrus have been focused on physiological processes relevant to fruit quality and productivity of the major species, especially sweet orange. Less attention has been paid to vegetative or reproductive tissues, while most Citrus species have never been analysed. In this work, we characterized the transcriptome of vegetative and reproductive tissues from 12 Citrus species from all main phylogenetic groups. Our aims were to acquire a complete view of the citrus transcriptome landscape, to improve previous functional annotations and to obtain genetic markers associated with genes of agronomic interest. 28 samples were used for RNA-Seq analysis, obtained from 12 Citrus species: C. medica, C. aurantifolia, C. limon, C. bergamia, C. clementina, C. deliciosa, C. reshni, C. maxima, C. paradisi, C. aurantium, C. sinensis and Poncirus trifoliata. Four different organs were analysed: root, phloem, leaf and flower. A total of 3421 million Illumina reads were produced and mapped against the reference C. clementina genome sequence. Transcript discovery pipeline revealed 3326 new genes, the number of genes with alternative splicing was increased to 19,739, and a total of 73,797 transcripts were identified. Differential expression studies between the four tissues showed that gene expression is overall related to the physiological function of the specific organs above any other variable. Variants discovery analysis revealed the presence of indels and SNPs in genes associated with fruit quality and productivity. Pivotal pathways in citrus such as those of flavonoids, flavonols, ethylene and auxin were also analysed in detail.

Major Quantitative Trait Loci and Putative Candidate Genes for Powdery Mildew Resistance and Fruit-Related Traits Revealed by an Intraspecific Genetic Map for Watermelon (Citrullus lanatus var. lanatus)

An intraspecific genetic map for watermelon was constructed using an F₂ population derived from ‘Arka Manik’ × ‘TS34’ and transcript sequence variants and quantitative trait loci (QTL) for resistance to powdery mildew (PMR), seed size (SS), and fruit shape (FS) were analyzed. The map consists of 14 linkage groups (LGs) defined by 174 cleaved amplified polymorphic sequences (CAPS), 2 derived-cleaved amplified polymorphic sequence markers, 20 sequence-characterized amplified regions, and 8 expressed sequence tag-simple sequence repeat markers spanning 1,404.3 cM, with a mean marker interval of 6.9 cM and an average of 14.6 markers per LG. Genetic inheritance and QTL analyses indicated that each of the PMR, SS, and FS traits is controlled by an incompletely dominant effect of major QTLs designated as pmr2.1, ss2.1, and fsi3.1, respectively. The pmr2.1, detected on chromosome 2 (Chr02), explained 80.0% of the phenotypic variation (LOD = 30.76). This QTL was flanked by two CAPS markers, wsb2-24 (4.00 cM) and wsb2-39 (13.97 cM). The ss2.1, located close to pmr2.1 and CAPS marker wsb2-13 (1.00 cM) on Chr02, explained 92.3% of the phenotypic variation (LOD = 68.78). The fsi3.1, detected on Chr03, explained 79.7% of the phenotypic variation (LOD = 31.37) and was flanked by two CAPS, wsb3-24 (1.91 cM) and wsb3-9 (7.00 cM). Candidate gene-based CAPS markers were developed from the disease resistance and fruit shape gene homologs located on Chr.02 and Chr03 and were mapped on the intraspecific map. Colocalization of these markers with the major QTLs indicated that watermelon orthologs of a nucleotide-binding site-leucine-rich repeat class gene containing an RPW8 domain and a member of SUN containing the IQ67 domain are candidate genes for pmr2.1 and fsi3.1, respectively. The results presented herein provide useful information for marker-assisted breeding and gene cloning for PMR and fruit-related traits.

Identification of candidate genes associated with fertility restoration of cytoplasmic male-sterility in onion (Allium cepa L.) using a combination of bulked segregant analysis and RNA-seq

A combination of BSA and RNA-seq was performed to identify candidates for the restorer-of-fertility gene in onion. The AcPMS1 involved in DNA mismatch repair was identified as the best candidate. To identify candidate genes of the restorer-of-fertility gene (Ms) responsible for fertility restoration of onion cytoplasmic male-sterility, a combined approach of bulked segregant analysis and RNA-seq was employed. From 32,674 de novo assembled contigs, 430 perfectly homozygous SNPs between male-fertile (MF) and male-sterile (MS) bulks were identified in 141 contigs. After verifying the homozygosity of the SNPs by PCR amplification and sequencing, the SNPs on 139 of the contigs were genotypes for the two recombinants which contained crossover events between the Ms locus and two tightly linked molecular markers. As a result, 30 contigs showing perfect linkage with the Ms locus in the large-sized segregating population were identified. Among them, 14 showed perfect linkage disequilibrium (LD) with the Ms locus, as determined by genotyping 251 domestic breeding lines. Furthermore, molecular markers tagging the 14 contigs also showed almost perfect LD with each other in 124 exotic accessions introduced from 21 countries, except for one accession which contained a crossover event by which the 14 markers were divided into two groups. After sequencing of the full-length cDNA of the 14 contigs showing perfect LD, the deduced amino acids sequences of the MF and MS alleles were compared. Four genes were shown to harbor putative critical amino acid changes in the known domains. Among them, the gene encoding PMS1, involved in the DNA mismatch repair pathway, was assumed to be the best candidate gene responsible for fertility restoration of male-sterility in onion.

De Novo Assembly and Annotation of the Chinese Chive (Allium tuberosum Rottler ex Spr.) Transcriptome Using the Illumina Platform

Chinese chive (A. tuberosum Rottler ex Spr.) is one of the most widely cultivated Allium species in China. However, minimal transcriptomic and genomic data are available to reveal its evolution and genetic diversity. In this study, de novo transcriptome sequencing was performed to produce large transcript sequences using an Illumina HiSeq 2000 instrument. We produced 51,968,882 high-quality clean reads and assembled them into 150,154 contigs. A total of 60,031 unigenes with an average length of 631 bp were identified. Of these, 36,523 unigenes were homologous to existing database sequences, 35,648 unigenes were annotated in the NCBI non-redundant (Nr) sequence database, and 23,509 unigenes were annotated in the Swiss-Prot database. A total of 26,798 unigenes were assigned to 57 Gene Ontology (GO) terms, and 13,378 unigenes were assigned to Cluster of Orthologous Group categories. Using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database, we mapped 21,361 unigenes onto 128 pathways. Furthermore, 2,125 sequences containing simple sequence repeats (SSRs) were identified. This new dataset provides the most comprehensive resource currently available for gene expression, gene discovery, and future genomic research on Chinese chive. The sequence resources developed in this study can be used to develop molecular markers that will facilitate further genetic research on Chinese chive and related species.