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.

De novo whole-genome assembly of Chrysanthemum makinoi, a key wild chrysanthemum.

Chrysanthemum is among the top 10 cut, potted, and perennial garden flowers in the world. Despite this, to date, only the genomes of two wild diploid chrysanthemums have been sequenced and assembled. Here, we present the most complete and contiguous chrysanthemum de novo assembly published so far, as well as a corresponding ab initio annotation. The cultivated hexaploid varieties are thought to originate from a hybrid of wild chrysanthemums, among which the diploid Chrysanthemum makinoi has been mentioned. Using a combination of Oxford Nanopore long reads, Pacific Biosciences long reads, Illumina short reads, Dovetail sequences, and a genetic map, we assembled 3.1 Gb of its sequence into nine pseudochromosomes, with an N50 of 330 Mb and a BUSCO complete score of 92.1%. Our ab initio annotation pipeline predicted 95,074 genes and marked 80.0% of the genome as repetitive. This genome assembly of C. makinoi provides an important step forward in understanding the chrysanthemum genome, evolution, and history.

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.

Fine mapping of a thrips resistance QTL in Capsicum and the role of diterpene glycosides in the underlying mechanism.

KEY MESSAGE: A major thrips resistance QTL in Capsicum was fine-mapped to a region of 0.4 Mbp, and a multidisciplinary approach has been used to study putative underlying mechanisms. Resistance to thrips is an important trait for pepper growers. These insects can cause extensive damage to fruits, flowers and leaves on field and greenhouse grown plants worldwide. Two independent studies in Capsicum identified diterpene glycosides as metabolites that are correlated with thrips resistance. In this study, we fine-mapped a previously defined thrips resistance QTL on chromosome 6, to a region of 0.4 Mbp harbouring 15 genes. Two of these 15 candidate genes showed differences in gene expression upon thrips induction, when comparing plants carrying the resistance allele in homozygous state to plants with the susceptibility allele in homozygous state for the QTL region. Three genes, including the two genes that showed difference in gene expression, contained a SNP that was predicted to lead to changes in protein structure. Therefore, these three genes, i.e. an acid phosphatase 1 (APS1), an organic cation/carnitine transporter 7 (OCT7) and an uncharacterized locus LOC107874801, are the most likely candidates for playing a role in thrips resistance and are a first step in elucidating the genetic basis of thrips resistance in Capsicum. In addition, we show that the diterpene glycoside profiles did not differ between plants with the resistance and susceptibility allele for the chromosome 6 QTL, suggesting that these compounds do not play a role in the resistance conferred by the genes located in the major thrips resistance QTL studied.

Exploration of a Resequenced Tomato Core Collection for Phenotypic and Genotypic Variation in Plant Growth and Fruit Quality Traits.

A tomato core collection consisting of 122 gene bank accessions, including landraces, old cultivars, and wild relatives, was explored for variation in several plant growth, yield and fruit quality traits. The resequenced accessions were also genotyped with respect to a number of mutations or variations in key genes known to underlie these traits. The yield-related traits fruit number and fruit weight were much higher in cultivated varieties when compared to wild accessions, while, in wild tomato accessions, Brix was higher than in cultivated varieties. Known mutations in fruit size and shape genes could well explain the fruit size variation, and fruit colour variation could be well explained by known mutations in key genes of the carotenoid and flavonoid pathway. The presence and phenotype of several plant architecture affecting mutations, such as self-pruning (sp), compound inflorescence (s), jointless-2 (j-2), and potato leaf (c) were also confirmed. This study provides valuable phenotypic information on important plant growth- and quality-related traits in this collection. The allelic distribution of known genes that underlie these traits provides insight into the role and importance of these genes in tomato domestication and breeding. This resource can be used to support (precision) breeding strategies for tomato crop improvement.

Solyntus, the New Highly Contiguous Reference Genome for Potato (Solanum tuberosum).

With the rapid expansion of the application of genomics and sequencing in plant breeding, there is a constant drive for better reference genomes. In potato (Solanum tuberosum), the third largest food crop in the world, the related species S. phureja, designated "DM", has been used as the most popular reference genome for the last 10 years. Here, we introduce the de novo sequenced genome of Solyntus as the next standard reference in potato genome studies. A true Solanum tuberosum made up of 116 contigs that is also highly homozygous, diploid, vigorous and self-compatible, Solyntus provides a more direct and contiguous reference then ever before available. It was constructed by sequencing with state-of-the-art long and short read technology and assembled with Canu. The 116 contigs were assembled into scaffolds to form each pseudochromosome, with three contigs to 17 contigs per chromosome. This assembly contains 93.7% of the single-copy gene orthologs from the Solanaceae set and has an N50 of 63.7 Mbp. The genome and related files can be found at https://www.plantbreeding.wur.nl/Solyntus/ With the release of this research line and its draft genome we anticipate many exciting developments in (diploid) potato research.

The ability to manipulate ROS metabolism in pepper may affect aphid virulence.

Myzus persicae has severe economic impact on pepper (Capsicum) cultivation. Previously, we identified two populations of M. persicae, NL and SW, that were avirulent and virulent, respectively on C. baccatum accession PB2013071. The transcriptomics approach used in the current study, which is the first study to explore the pepper-aphid interaction at the whole genome gene expression level, revealed genes whose expression is differentially regulated in pepper accession PB2013071 upon infestation with these M. persicae populations. The NL population induced ROS production genes, while the SW population induced ROS scavenging genes and repressed ROS production genes. We also found that the SW population can induce the removal of ROS which accumulated in response to preinfestion with the NL population, and that preinfestation with the SW population significantly improved the performance of the NL population. This paper supports the hypothesis that M. persicae can overcome the resistance in accession PB2013071 probably because of its ability to manipulate plant defense response especially the ROS metabolism and such ability may benefit avirulent conspecific aphids.

Aphid resistance in Capsicum maps to a locus containing LRR-RLK gene analogues.

KEY MESSAGE: A QTL for aphid resistance on pepper chromosome 2 was identified and validated. This QTL affects aphid survival and reproduction, and was fine mapped to a locus containing LRR-RLK analogues. Myzus persicae is one of the most threatening insect pests that adversely affects pepper (Capsicum) cultivation. Resistance to aphids was previously identified in Capsicum baccatum. This study aimed at elucidating the genetics of aphid resistance in C. baccatum. A QTL analysis was carried out for M. persicae resistance in an F2 population derived from an intraspecific cross between a highly resistant plant and a susceptible plant. Survival and reproduction were used as resistance parameters. Interval mapping detected two QTLs affecting aphid survival (Rmpas-1) and reproduction (Rmprp-1), respectively, both localized in the same area and sharing the same top marker on chromosome 2. Use of this marker as co-factor in multiple-QTL mapping analysis revealed a second, minor QTL (Rmprp-2) only affecting aphid reproduction, on chromosome 4. Fine mapping confirmed the effects of Rmpas-1 and Rmprp-1 and narrowed the major QTL Rmprp-1 down to a genomic region of 96 kb which is predicted to encode four analogues of resistance genes of the receptor-like kinase family containing a leucine-rich repeat domain (LRR-RLKs). This work provides not only initial information for breeding aphid-resistant pepper varieties, but also forms the basis for future molecular analysis of gene(s) involved in aphid resistance.

Patterns of Transmission Ratio Distortion in Interspecific Lettuce Hybrids Reveal a Sex-Independent Gametophytic Barrier.

Interspecific crosses can result in progeny with reduced vitality or fertility due to genetic incompatibilities between species, a phenomenon known as hybrid incompatibility (HI). HI is often caused by a bias against deleterious allele combinations, which results in transmission ratio distortion (TRD). Here, we determined the genome-wide distribution of HI between wild lettuce, Lactuca saligna, and cultivated lettuce, L. sativa, in a set of backcross inbred lines (BILs) with single introgression segments from L. saligna introgressed into a L. sativa genetic background. Almost all BILs contained an introgression segment in a homozygous state except a few BILs, for which we were able to obtain only a single heterozygous introgression. Their inbred progenies displayed severe TRD with a bias toward the L. sativa allele and complete nontransmission of the homozygous L. saligna introgression, i.e., absolute HI. These HI might be caused by deleterious heterospecific allele combinations at two loci. We used an multilocus segregating interspecific F2 population to identify candidate conspecific loci that can nullify the HI in BILs. Segregation analysis of developed double-introgression progenies showed nullification of three HI and proved that these HI are explained by nuclear pairwise incompatibilities. One of these digenic HI showed 29% reduced seed set and its pattern of TRD pointed to a sex-independent gametophytic barrier. Namely, this HI was caused by complete nontransmission of one heterospecific allele combination at the haploid stage, surprisingly in both male and female gametophytes. Our study shows that two-locus incompatibility systems contribute to reproductive barriers among Lactuca species.

Combining QTL mapping with transcriptome and metabolome profiling reveals a possible role for ABA signaling in resistance against the cabbage whitefly in cabbage.

Whiteflies are among the world's most significant agricultural pests and chemical insecticides are extensively used to reduce crop damage to acceptable levels. However, nearly all insecticides pose a threat to the environment and alternative control methods, such as breeding of crop varieties that are inherently insect-resistant, are needed. Previously, a strong source of plant-age dependent resistance to the cabbage whitefly (Aleyrodes proletella) has been identified in the modern white cabbage (Brassica oleracea var. capitata) variety Rivera. However, nothing is known about the molecular mechanisms or the genes involved in this resistance. In the present study, a multidisciplinary approach combining transcriptome and metabolome profiling with genetic mapping was used to identify the molecular players of whitefly resistance in cabbage. Transcriptome profiles of young (susceptible) and older (resistant) Rivera plants were analyzed using RNA sequencing. While many genes involved in general processes were differentially expressed between both ages, several defense-related processes were overrepresented in the transcriptome profile of older plants. Hormone measurements revealed that jasmonic acid (JA) levels decreased upon whitefly infestation at both plant ages. Interestingly, abscisic acid (ABA) levels showed contrasting effects in response to whitefly infestation: ABA levels were reduced in young plants but induced in older plants upon whitefly feeding. Auxin levels were significantly lower in older plants compared with young plants, independent of whitefly presence, while glucosinolate levels were higher. Additionally, whitefly performance was monitored in an F2 population derived from a cross between Rivera and the susceptible white cabbage variety Christmas Drumhead. Significant QTL intervals were mapped on chromosome 2 and 9 for oviposition rate and whitefly adult survival, respectively. Several genes that were higher expressed in older plants and located in the identified QTL intervals were orthologous to Arabidopsis genes that have been related to ABA signaling, suggesting a role for ABA in the regulation of resistance towards whiteflies. Our results show that combining different omics approaches is a useful strategy to identify candidate genes underlying insect resistance.

Bidirectional backcrosses between wild and cultivated lettuce identify loci involved in nonhost resistance to downy mildew.

KEY MESSAGE: The nonhost resistance of wild lettuce to lettuce downy mildew seems explained by four components of a putative set of epistatic genes. The commonplace observation that plants are immune to most potential pathogens is known as nonhost resistance (NHR). The genetic basis of NHR is poorly understood. Inheritance studies of NHR require crosses of nonhost species with a host, but these crosses are usually unsuccessful. The plant-pathosystem of lettuce and downy mildew, Bremia lactucae, provides a rare opportunity to study the inheritance of NHR, because the nonhost wild lettuce species Lactuca saligna is sufficiently cross-compatible with the cultivated host Lactuca sativa. Our previous studies on NHR in one L. saligna accession led to the hypothesis that multi-locus epistatic interactions might explain NHR. Here, we studied NHR at the species level in nine accessions. Besides the commonly used approach of studying a target trait from a wild donor species in a cultivar genetic background, we also explored the opposite, complementary approach of cultivar introgression in a wild species background. This bidirectional approach encompassed (1) nonhost into host introgression: identification of L. saligna derived chromosome regions that were overrepresented in highly resistant BC1 plants (F1 × L. sativa), (2) host into nonhost introgression: identification of L. sativa derived chromosome regions that were overrepresented in BC1 inbred lines (F1 × L. saligna) with relatively high infection levels. We demonstrated that NHR is based on resistance factors from L. saligna and the genetic dose for NHR differs between accessions. NHR seemed explained by combinations of epistatic genes on three or four chromosome segments, of which one chromosome segment was validated by the host into nonhost approach.

SNP-markers in Allium species to facilitate introgression breeding in onion.

BACKGROUND: Within onion, Allium cepa L., the availability of disease resistance is limited. The identification of sources of resistance in related species, such as Allium roylei and Allium fistulosum, was a first step towards the improvement of onion cultivars by breeding. SNP markers linked to resistance and polymorphic between these related species and onion cultivars are a valuable tool to efficiently introgress disease resistance genes. In this paper we describe the identification and validation of SNP markers valuable for onion breeding. RESULTS: Transcriptome sequencing resulted in 192 million RNA seq reads from the interspecific F1 hybrid between A. roylei and A. fistulosum (RF) and nine onion cultivars. After assembly, reliable SNPs were discovered in about 36 % of the contigs. For genotyping of the interspecific three-way cross population, derived from a cross between an onion cultivar and the RF (CCxRF), 1100 SNPs that are polymorphic in RF and monomorphic in the onion cultivars (RF SNPs) were selected for the development of KASP assays. A molecular linkage map based on 667 RF-SNP markers was constructed for CCxRF. In addition, KASP assays were developed for 1600 onion-SNPs (SNPs polymorphic among onion cultivars). A second linkage map was constructed for an F2 of onion x A. roylei (F2(CxR)) that consisted of 182 onion-SNPs and 119 RF-SNPs, and 76 previously mapped markers. Markers co-segregating in both the F2(CxR) and the CCxRF population were used to assign the linkage groups of RF to onion chromosomes. To validate usefulness of these SNP markers, QTL mapping was applied in the CCxRF population that segregates for resistance to Botrytis squamosa and resulted in a QTL for resistance on chromosome 6 of A. roylei. CONCLUSIONS: Our research has more than doubled the publicly available marker sequences of expressed onion genes and two onion-related species. It resulted in a detailed genetic map for the interspecific CCxRF population. This is the first paper that reports the detection of a QTL for resistance to B. squamosa in A. roylei.

QualitySNPng: a user-friendly SNP detection and visualization tool.

QualitySNPng is a new software tool for the detection and interactive visualization of single-nucleotide polymorphisms (SNPs). It uses a haplotype-based strategy to identify reliable SNPs; it is optimized for the analysis of current RNA-seq data; but it can also be used on genomic DNA sequences derived from next-generation sequencing experiments. QualitySNPng does not require a sequenced reference genome and delivers reliable SNPs for di- as well as polyploid species. The tool features a user-friendly interface, multiple filtering options to handle typical sequencing errors, support for SAM and ACE files and interactive visualization. QualitySNPng produces high-quality SNP information that can be used directly in genotyping by sequencing approaches for application in QTL and genome-wide association mapping as well as to populate SNP arrays. The software can be used as a stand-alone application with a graphical user interface or as part of a pipeline system like Galaxy. Versions for Windows, Mac OS X and Linux, as well as the source code, are available from http://www.bioinformatics.nl/QualitySNPng.

Generation and analysis of expressed sequence tags in the extreme large genomes Lilium and Tulipa.

BACKGROUND: Bulbous flowers such as lily and tulip (Liliaceae family) are monocot perennial herbs that are economically very important ornamental plants worldwide. However, there are hardly any genetic studies performed and genomic resources are lacking. To build genomic resources and develop tools to speed up the breeding in both crops, next generation sequencing was implemented. We sequenced and assembled transcriptomes of four lily and five tulip genotypes using 454 pyro-sequencing technology. RESULTS: Successfully, we developed the first set of 81,791 contigs with an average length of 514 bp for tulip, and enriched the very limited number of 3,329 available ESTs (Expressed Sequence Tags) for lily with 52,172 contigs with an average length of 555 bp. The contigs together with singletons covered on average 37% of lily and 39% of tulip estimated transcriptome. Mining lily and tulip sequence data for SSRs (Simple Sequence Repeats) showed that di-nucleotide repeats were twice more abundant in UTRs (UnTranslated Regions) compared to coding regions, while tri-nucleotide repeats were equally spread over coding and UTR regions. Two sets of single nucleotide polymorphism (SNP) markers suitable for high throughput genotyping were developed. In the first set, no SNPs flanking the target SNP (50 bp on either side) were allowed. In the second set, one SNP in the flanking regions was allowed, which resulted in a 2 to 3 fold increase in SNP marker numbers compared with the first set. Orthologous groups between the two flower bulbs: lily and tulip (12,017 groups) and among the three monocot species: lily, tulip, and rice (6,900 groups) were determined using OrthoMCL. Orthologous groups were screened for common SNP markers and EST-SSRs to study synteny between lily and tulip, which resulted in 113 common SNP markers and 292 common EST-SSR. Lily and tulip contigs generated were annotated and described according to Gene Ontology terminology. CONCLUSIONS: Two transcriptome sets were built that are valuable resources for marker development, comparative genomic studies and candidate gene approaches. Next generation sequencing of leaf transcriptome is very effective; however, deeper sequencing and using more tissues and stages is advisable for extended comparative studies.

Diversity arrays technology (DArT) markers in apple for genetic linkage maps.

Diversity Arrays Technology (DArT) provides a high-throughput whole-genome genotyping platform for the detection and scoring of hundreds of polymorphic loci without any need for prior sequence information. The work presented here details the development and performance of a DArT genotyping array for apple. This is the first paper on DArT in horticultural trees. Genetic mapping of DArT markers in two mapping populations and their integration with other marker types showed that DArT is a powerful high-throughput method for obtaining accurate and reproducible marker data, despite the low cost per data point. This method appears to be suitable for aligning the genetic maps of different segregating populations. The standard complexity reduction method, based on the methylation-sensitive PstI restriction enzyme, resulted in a high frequency of markers, although there was 52-54% redundancy due to the repeated sampling of highly similar sequences. Sequencing of the marker clones showed that they are significantly enriched for low-copy, genic regions. The genome coverage using the standard method was 55-76%. For improved genome coverage, an alternative complexity reduction method was examined, which resulted in less redundancy and additional segregating markers. The DArT markers proved to be of high quality and were very suitable for genetic mapping at low cost for the apple, providing moderate genome coverage. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11032-011-9579-5) contains supplementary material, which is available to authorized users.

Identification and functional characterization of the Lactococcus lactis rfb operon, required for dTDP-rhamnose Biosynthesis.

dTDP-rhamnose is an important precursor of cell wall polysaccharides and rhamnose-containing exopolysaccharides (EPS) in Lactococcus lactis. We cloned the rfbACBD operon from L. lactis MG1363, which comprises four genes involved in dTDP-rhamnose biosynthesis. When expressed in Escherichia coli, the lactococcal rfbACBD genes could sustain heterologous production of the Shigella flexneri O antigen, providing evidence of their functionality. Overproduction of the RfbAC proteins in L. lactis resulted in doubled dTDP-rhamnose levels, indicating that the endogenous RfbAC activities control the intracellular dTDP-rhamnose biosynthesis rate. However, RfbAC overproduction did not affect rhamnose-containing B40-EPS production levels. A nisin-controlled conditional RfbBD mutant was unable to grow in media lacking the inducer nisin, indicating that the rfb genes have an essential role in L. lactis. Limitation of RfbBD activities resulted in the production of altered EPS. The monomeric sugar of the altered EPS consisted of glucose, galactose, and rhamnose at a molar ratio of 1:0.3:0.2, which is clearly different from the ratio in the native sugar. Biophysical analysis revealed a fourfold-greater molecular mass and a twofold-smaller radius of gyration for the altered EPS, indicating that these EPS are more flexible polymers with changed viscosifying properties. This is the first indication that enzyme activity at the level of central carbohydrate metabolism affects EPS composition.