Construction of a high-density genetic map and the X/Y sex-determining gene mapping in spinach based on large-scale markers developed by specific-locus amplified fragment sequencing (SLAF-seq)

Fine mapping and analysis of a candidate gene controlling Phytophthora blight resistance in cucumber

Cucumber blight is a destructive disease. The best way to control this disease is resistance breeding. This study focuses on disease resistance gene mapping and molecular marker development. We used the resistant cucumber, JSH, and susceptible cucumber, B80, as parents to construct F2 populations. Bulked segregant analysis (BSA) combined with specific length amplified fragment sequencing (SLAF-seq) were used, from which we developed cleaved amplified polymorphic sequence (CAPs) markers to map the resistance gene. Resistance in F2 individuals showed a segregation ratio of resistance:susceptibility close to 3:1. The gene in JSH resistant cucumber was mapped to an interval of 9.25 kb, and sequencing results for the three genes in the mapped region revealed three mutations at base sites 225, 302, and 591 in the coding region of Csa5G139130 between JSH and B80, but no mutations in coding regions of Csa5G139140 and Csa5G139150. The mutations caused changes in amino acids 75 and 101 of the protein encoded by Csa5G139130, suggesting that Csa5G139130 is the most likely resistance candidate gene. We developed a molecular marker, CAPs-4, as a closely linked marker for the cucumber blight resistance gene. This is the first report on mapping of a cucumber blight resistance gene and will provideg a useful marker for molecular breeding of cucumber resistance to Phytophthora blight.

Genome-Wide Identification and Characterization of MYB Gene Family and Analysis of Its Sex-Biased Expression Pattern in Spinacia oleracea L

The members of the myeloblastosis (MYB) family of transcription factors (TFs) participate in a variety of biological regulatory processes in plants, such as circadian rhythm, metabolism, and flower development. However, the characterization of MYB genes across the genomes of spinach Spinacia oleracea L. has not been reported. Here, we identified 140 MYB genes in spinach and described their characteristics using bioinformatics approaches. Among the MYB genes, 54 were 1R-MYB, 80 were 2R-MYB, 5 were 3R-MYB, and 1 was 4R-MYB. Almost all MYB genes were located in the 0-30 Mb region of autosomes; however, the 20 MYB genes were enriched at both ends of the sex chromosome (chromosome 4). Based on phylogeny, conserved motifs, and the structure of genes, 2R-MYB exhibited higher conservation relative to 1R-MYB genes. Tandem duplication and collinearity of spinach MYB genes drive their evolution, enabling the functional diversification of spinach genes. Subcellular localization prediction indicated that spinach MYB genes were mainly located in the nucleus. Cis-acting element analysis confirmed that MYB genes were involved in various processes of spinach growth and development, such as circadian rhythm, cell differentiation, and reproduction through hormone synthesis. Furthermore, through the transcriptome data analysis of male and female flower organs at five different periods, ten candidate genes showed biased expression in spinach males, suggesting that these genes might be related to the development of spinach anthers. Collectively, this study provides useful information for further investigating the function of MYB TFs and novel insights into the regulation of sex determination in spinach.

Evolution of the spinach sex-linked region within a rarely recombining pericentromeric region

Sex chromosomes have evolved independently in many different plant lineages. Here, we describe reference genomes for spinach (Spinacia oleracea) X and Y haplotypes by sequencing homozygous XX females and YY males. The long arm of 185-Mb chromosome 4 carries a 13-Mb X-linked region (XLR) and 24.1-Mb Y-linked region (YLR), of which 10 Mb is Y specific. We describe evidence that this reflects insertions of autosomal sequences creating a "Y duplication region" or "YDR" whose presence probably directly reduces genetic recombination in the immediately flanking regions, although both the X and Y sex-linked regions are within a large pericentromeric region of chromosome 4 that recombines rarely in meiosis of both sexes. Sequence divergence estimates using synonymous sites indicate that YDR genes started diverging from their likely autosomal progenitors about 3 MYA, around the time when the flanking YLR stopped recombining with the XLR. These flanking regions have a higher density of repetitive sequences in the YY than the XX assembly and include slightly more pseudogenes compared with the XLR, and the YLR has lost about 11% of the ancestral genes, suggesting some degeneration. Insertion of a male-determining factor would have caused Y linkage across the entire pericentromeric region, creating physically small, highly recombining, terminal pseudoautosomal regions. These findings provide a broader understanding of the origin of sex chromosomes in spinach.

Identification of Sex Differentiation-Related microRNAs in Spinach Female and Male Flower

Sex determination and differentiation is an important biological process for unisexual flower development. Spinach is a model plant to study the mechanism of sex determination and differentiation of dioecious plant. Till now, little is known about spinach sex determination and differentiation mechanism. MicroRNAs are key factors in flower development. Herein, small RNA sequencing was performed to explore the roles of microRNAs in spinach sex determination and differentiation. As a result, 92 known and 3402 novel microRNAs were identified in 18 spinach female and male flower samples. 74 differentially expressed microRNAs were identified between female and male flowers, including 20 female-biased and 48 male-biased expression microRNAs. Target prediction identified 22 sex-biased microRNA-target pairs, which may be involved in spinach sex determination or differentiation. Among the differentially expressed microRNAs between FNS and M03, 55 microRNAs were found to reside in sex chromosome; one of them, sol-miR2550n, was functionally studied via genetic transformation. Silencing of sol-miR2550n resulted in abnormal anther while overexpression of sol-miR2550n induced early flowering, indicating sol-miR2550n was a male-promoting factor and validating the reliability of our small RNA sequencing data. Conclusively, this work can supply valuable information for exploring spinach sex determination and differentiation and provide a new insight in studying unisexual flower development.

Quantitative Trait Locus Mapping and Identification of Candidate Genes Controlling Bolting in Spinach (Spinacia oleracea L.)

Spinach is a typical light-sensitive plant. Long days can induce early bolting, thereby influencing the regional adaptation, quality, and vegetative yield of spinach. However, the genes and genetic mechanisms underlying this trait in spinach remain unclear. In this study, a major quantitative trait locus (QTL) qBT1.1, was mapped on chromosome 1 using a BC₁ population (BC_1a) derived from 12S3 (late-bolting recurrent lines) and 12S4 (early bolting lines) with specific-locus amplified fragment (SLAF) markers and Kompetitive Allele Specific PCR (KASP) markers. The qBT1.1 locus was further confirmed and narrowed down to 0.56 Mb by using a large BC₁ (BC_1b) population and an F₂ population using the above KASP markers and the other 20 KASP markers. Within this region, two putative genes, namely, SpFLC and SpCOL14, were of interest due to their relationship with flower regulatory pathways. For SpCOL14, we found multiple variations in the promoter, and the expression pattern was consistent with bolting stages. SpCOL14 was therefore assumed to the best candidate gene for bolting. Overall, our results provide a basis for understanding the molecular mechanisms of bolting in spinach and contribute to the breeding of diverse spinach germplasms for adaptation to different regions.

Identification of structural variation and polymorphisms of a sex co-segregating scaffold in spinach

Spinach is a common vegetable, and dioecy is maintained by a pair of XY sex chromosomes. Due to limited genomic resources and its highly repetitive genome, limited studies were conducted to investigate the genomic landscape of the region near sex-determining loci. In this study, we screened the structure variations (SVs) between Y-linked contigs and a 1.78-Mb X scaffold (Super_scaffold 66), which enabled the development of 12 sex co-segregating DNA markers. These markers were tested in one F₁ mapping population and 40 spinach accessions, which comprised 692 individual plants with the strong sex linkage pattern. In addition, we found that Super_scaffold 66 was highly repetitive along with the enriched LTR-RTs insertions and decreased microsatellite distribution compared with the rest genome, which matches extremely low gene density featured by only nine annotated genes. Synteny analysis between Y contigs and Superscaffold_66 revealed a 340-Kb accumulative Y contig (non-continuous) and a 500-Kb X counterpart along with SVs and wide-spread tandem duplications. Among the nine genes, one ABC transporter gene revealed noticeable SVs between Y contig and X counterpart, as an approximate 5-Kb recent Gypsy LTR-RT insertion in the Y-linked allele, but not the X allele. The gene paucity, SVs, and sex-linked polymorphisms attributed to the recombination suppression. We proposed that Super_scaffold 66 is part of the non-recombining region containing the sex determination genes. The spread of 12 sex co-segregating markers from this 1.78 Mb genomic region indicated the existence and expansion of sex determination region during progression of the Y chromosome.

High resolution mapping and candidate gene identification of downy mildew race 16 resistance in spinach

Background

Downy mildew, the most devastating disease of spinach (Spinacia oleracea L.), is caused by the oomycete Peronospora effusa [=P. farinosa f. sp. spinaciae]. The P. effusa shows race specificities to the resistant host and comprises 19 reported races and many novel isolates. Sixteen new P. effusa races were identified during the past three decades, and the new pathogen races are continually overcoming the genetic resistances used in commercial cultivars. A spinach breeding population derived from the cross between cultivars Whale and Lazio was inoculated with P. effusa race 16 in an environment-controlled facility; disease response was recorded and genotyped using genotyping by sequencing (GBS). The main objective of this study was to identify resistance-associated single nucleotide polymorphism (SNP) markers from the cultivar Whale against the P. effusa race 16.

Results

Association analysis conducted using GBS markers identified six significant SNPs (S3_658,306, S3_692697, S3_1050601, S3_1227787, S3_1227802, S3_1231197). The downy mildew resistance locus from cultivar Whale was mapped to a 0.57 Mb region on chromosome 3, including four disease resistance candidate genes (Spo12736, Spo12784, Spo12908, and Spo12821) within 2.69–11.28 Kb of the peak SNP.

Conclusions

Genomewide association analysis approach was used to map the P. effusa race 16 resistance loci and identify associated SNP markers and the candidate genes. The results from this study could be valuable in understanding the genetic basis of downy mildew resistance, and the SNP marker will be useful in spinach breeding to select resistant lines.

Quantitative trait loci (QTL) analysis of leaf related traits in spinach (Spinacia oleracea L.)

BACKGROUND

Spinach (Spinacia oleracea L.) is an important leafy vegetable crop, and leaf-related traits including leaf length, leaf width, and petiole length, are important commercial traits. However, the underlying genes remain unclear. The objective of the study was to conduct QTL mapping of leaf-related traits in spinach.

RESULTS

A BC₁ population was used to construct the linkage map and for QTL mapping of leaf length, leaf width, petiole length, and the ratio of leaf length to width in 2015 and 2019. Two genetic linkage maps were constructed by specific locus amplified fragment sequencing (SLAF-seq), and kompetitive allele specific PCR (KASP) technology, respectively using BC₁ population in 2015. Based on the results of 2015, the specific linkage groups (LG) detected QTLs were generated using BC₁ population in 2019. A total of 13 QTLs were detected for leaf-related traits, only five QTLs being repeatedly detected in multiple years or linkage maps. Interestingly, the major QTLs of leaf length, petiole length, and the ratio of leaf length to width were highly associated with the same SNP markers (KM3102838, KM1360385 and KM2191098). A major QTL of leaf width was mapped on chromosome 1 from 41.470-42.045 Mb. And 44 genes were identified within the region. Based on the GO analysis, these genes were significantly enriched on ribonuclease, lyase activity, phosphodiester bond hydrolysis process, and cell wall component, thus it might change cell size to determine leaves shape.

CONCLUSIONS

Five QTLs for leaf-related traits were repeatedly detected at least two years or linkage maps. The major QTLs of leaf length, petiole length, and the ratio of leaf length to width were mapped on the same loci. And three genes (Spo10792, Spo21018, and Spo21019) were identified as important candidate genes for leaf width.

Fine mapping and molecular marker development of the Fs gene controlling fruit spines in spinach (Spinacia oleracea L.)

The Fs gene, which controls spinach fruit spines, was fine mapped to a 0.27 Mb interval encompassing four genes on chromosome 3. There are two types of fruit of spinach (Spinacia oleracea L.), spiny and spineless, which are visually distinguishable by the spines of fruit coat. In spinach breeding, the fruit characteristic is an important agronomic trait that have impacts on "seed" treatment and mechanized sowing. However, the gene(s) controlling the fruit spiny trait have not been characterized and the genetic mechanism of this trait remained unclear. The objectives of the study were to fine map the gene controlling fruit spines and develop molecular markers for marker-assisted selection purpose. Genetic analysis of the spiny trait in segregating populations indicated that fruit spines were controlled by a single dominant gene, designated as Fs. Using a super-BSA method and recombinants analysis in a BC₁ population, Fs was mapped to a 1.9-Mb interval on chromosome 3. The Fs gene was further mapped to a 0.27-Mb interval using a recombinant inbred line (RIL) population with 120 lines. From this 0.27 Mb region, four candidate genes were identified in the reference genome. The structure and expression of the four genes were compared between the spiny and spineless parents. A co-dominant marker YC-15 was found to be co-segregating with the fruit spines trait, which produced a 129-bp fragment specific to spiny trait and a 108-bp fragment for spineless fruit. This marker can predict spiny trait with a 94.8% accuracy rate when tested with 100 diverse germplasm, suggesting that this marker would be valuable for marker-assisted selection in spinach breeding.

Transcriptome profiling of differentially expressed genes of male and female inflorescences in spinach (Spinacia oleracea L.)

Spinach (Spinacia oleracea L.) is commonly considered a dioecious plant with heterogametic (XY) and homogametic (XX) sex chromosomes. The characteristic is also utilized for the production of spinach hybrid seeds. However, the molecular mechanisms of sex determination in spinach are still unclear because of a lack of genomic and transcriptomic information. In this study, RNA sequencing (RNA-seq) was performed in male and female inflorescences to provide insight into the molecular basis of sex determination in spinach. Comparative transcriptome analyses showed that 2278 differentially expressed genes (DEGs) were identified between male and female inflorescences. A high correlation between the RNA-Seq and qRT-PCR validation for DEGs was observed. Among these, 182 DEGs were annotated to transcription factors including the MYB family protein, bHLH family, and MADS family, suggesting these factors might play a vital role in sex determination. Moreover, 26 DEGs related to flower development, including nine ABCE class genes, were detected. Expression analyses of hormone pathways showed that brassinosteroids may be key hormones related to sex determination in spinach. Overall, this study provides a large amount of DEGs related to sexual expression and lays a foundation for unraveling the regulatory mechanism of sex determination in spinach.

Construction of high-density genetic maps defined sex determination region of the Y chromosome in spinach

Spinach (Spinacia olracea L.) is a dioecious leafy vegetable with a highly repetitive genome of around 990 Mb, which is challenging for de-novo genome assembly. In our study, a segregating F₁ (double pseudo-testcross) population from 'Viroflay' × 'Cornell-NO. 9' was used for genetic mapping by resequencing genotyping. In the paternal 'Cornell-NO. 9' map, 212,414 SNPs were mapped, and the total linkage distance was 476.83 cM; the maternal 'Viroflay' map included 29,282 SNPs with 401.28 cM total genetic distance. Both paternal and maternal maps have the expected number of six linkage groups (LGs). A non-recombining region with 5678 SNPs (39 bin markers) co-segregates with sex type which located at 45.2 cM of LG1 in the 'Cornell-NO. 9' map while indicates the sex determination region (SDR). Integration of two maps into a consensus map guided us to anchor additional 1242 contigs to six pseudomolecules from the published reference genome, which improved additional 233 Mb (23.4%) assembly based on spinach estimated genome size. Particularly, the X counterpart of SDR in our assembly is estimated around 18.4 Mb which locates at the largest chromosome, as consensus with sex-biased FISH signals from previous cytogenetics studies. The region is featured by reduced gene density, higher percentage of repetitive sequences, and no recombination. Our linkage maps provide the resource for improving spinach genome de-novo assembly and identification of sex-determining genes in spinach.

Identification of a quantitative trait loci (QTL) associated with ammonia tolerance in the Pacific white shrimp (Litopenaeus vannamei)

Background

Ammonia is one of the most common toxicological environment factors affecting shrimp health. Although ammonia tolerance in shrimp is closely related to successful industrial production, few genetic studies of this trait are available.

Results

In this study, we constructed a high-density genetic map of the Pacific white shrimp (Litopenaeus vannamei) using specific length amplified fragment sequencing (SLAF-seq). The constructed genetic map contained 17,338 polymorphic markers spanning 44 linkage groups, with a total distance of 6360.12 centimorgans (cM) and an average distance of 0.37 cM. Using this genetic map, we identified a quantitative trait locus (QTL) that explained 7.41–8.46% of the phenotypic variance in L. vannamei survival time under acute ammonia stress. We then sequenced the transcriptomes of the most ammonia-tolerant and the most ammonia-sensitive individuals from each of four genetically distinct L. vannamei families. We found that 7546 genes were differentially expressed between the ammonia-tolerant and ammonia-sensitive individuals. Using QTL analysis and the transcriptomes, we identified one candidate gene (annotated as an ATP synthase g subunit) associated with ammonia tolerance.

Conclusions

In this study, we constructed a high-density genetic map of L. vannamei and identified a QTL for ammonia tolerance. By combining QTL and transcriptome analyses, we identified a candidate gene associated with ammonia tolerance. Our work provides the basis for future genetic studies focused on molecular marker-assisted selective breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07254-x.

Comparative transcriptome analysis of male and female flowers in Spinacia oleracea L

Background

Dioecious spinach (Spinacia oleracea L.), a commercial and nutritional vegetable crop, serves as a model for studying the mechanisms of sex determination and differentiation in plants. However, this mechanism is still unclear. Herein, based on PacBio Iso-seq and Illumina RNA-seq data, comparative transcriptome analysis of male and female flowers were performed to explore the sex differentiation mechanism in spinach.

Results

Compared with published genome of spinach, 10,800 transcripts were newly annotated; alternative splicing, alternative polyadenylation and lncRNA were analyzed for the first time, increasing the diversity of spinach transcriptome. A total of 2965 differentially expressed genes were identified between female and male flowers at three early development stages. The differential expression of RNA splicing-related genes, polyadenylation-related genes and lncRNAs suggested the involvement of alternative splicing, alternative polyadenylation and lncRNA in sex differentiation. Moreover, 1946 male-biased genes and 961 female-biased genes were found and several candidate genes related to gender development were identified, providing new clues to reveal the mechanism of sex differentiation. In addition, weighted gene co-expression network analysis showed that auxin and gibberellin were the common crucial factors in regulating female or male flower development; however, the closely co-expressed genes of these two factors were different between male and female flower, which may result in spinach sex differentiation.

Conclusions

In this study, 10,800 transcripts were newly annotated, and the alternative splicing, alternative polyadenylation and long-noncoding RNA were comprehensively analyzed for the first time in spinach, providing valuable information for functional genome study. Moreover, candidate genes related to gender development were identified, shedding new insight on studying the mechanism of sex determination and differentiation in plant.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07277-4.

Identification of candidate genes for key fibre-related QTLs and derivation of favourable alleles in Gossypium hirsutum recombinant inbred lines with G. barbadense introgressions

Fine mapping QTLs and identifying candidate genes for cotton fibre-quality and yield traits would be beneficial to cotton breeding. Here, we constructed a high-density genetic map by specific-locus amplified fragment sequencing (SLAF-seq) to identify QTLs associated with fibre-quality and yield traits using 239 recombinant inbred lines (RILs), which was developed from LMY22 (a high-yield Gossypium hirsutumL. cultivar) × LY343 (a superior fibre-quality germplasm with G. barbadenseL. introgressions). The genetic map spanned 3426.57 cM, including 3556 SLAF-based SNPs and 199 SSR marker loci. A total of 104 QTLs, including 67 QTLs for fibre quality and 37 QTLs for yield traits, were identified with phenotypic data collected from 7 environments. Among these, 66 QTLs were co-located in 19 QTL clusters on 12 chromosomes, and 24 QTLs were detected in three or more environments and determined to be stable. We also investigated the genomic components of LY343 and their contributions to fibre-related traits by deep sequencing the whole genome of LY343, and we found that genomic components from G. hirsutum races (which entered LY343 via its G. barbadense parent) contributed more favourable alleles than those from G. barbadense. We further identified six putative candidate genes for stable QTLs, including Gh_A03G1147 (GhPEL6), Gh_D07G1598 (GhCSLC6) and Gh_D13G1921 (GhTBL5) for fibre-length QTLs and Gh_D03G0919 (GhCOBL4), Gh_D09G1659 (GhMYB4) and Gh_D09G1690 (GhMYB85) for lint-percentage QTLs. Our results provide comprehensive insight into the genetic basis of the formation of fibre-related traits and would be helpful for cloning fibre-development-related genes as well as for marker-assisted genetic improvement in cotton.

Young sex chromosomes in plants and animals

A major reason for studying plant sex chromosomes is that they may often be 'young' systems. There is considerable evidence for the independent evolution of separate sexes within plant families or genera, in some cases showing that the maximum possible time during which their sex-determining genes have existed must be much shorter than those of several animal taxa. Consequently, their sex-linked regions could either have evolved soon after genetic sex determination arose or considerably later. Plants, therefore, include species with both young and old systems. I review several questions about the evolution of sex-determining systems and sex chromosomes that require studies of young systems, including: the kinds of mutations involved in the transition to unisexual reproduction from hermaphroditism or monoecy (a form of functional hermaphroditism); the times when they arose; and the extent to which the properties of sex-linked regions of genomes reflect responses to new selective situations created by the presence of a sex-determining locus. I also evaluate which questions are best studied in plants, vs other suitable candidate organisms. Studies of young plant systems can help understand general evolutionary processes that are shared with the sex chromosomes of other organisms.

Construction, characteristics and high throughput molecular screening methodologies in some special breeding populations: a horticultural perspective

Advanced marker technologies are widely used for evaluation of genetic diversity in cultivated crops, wild ancestors, landraces or any special plant genotypes. Developing agricultural cultivars requires the following steps: (i) determining desired characteristics to be improved, (ii) screening genetic resources to help find a superior cultivar, (iii) intercrossing selected individuals, (iv) generating genetically hybrid populations and screening them for agro-morphological or molecular traits, (v) evaluating the superior cultivar candidates, (vi) testing field performance at different locations, and (vii) certifying. In the cultivar development process valuable genes can be identified by creating special biparental or multiparental populations and analysing their association using suitable markers in given populations. These special populations and advanced marker technologies give us a deeper knowledge about the inherited agronomic characteristics. Unaffected by the changing environmental conditions, these provide a higher understanding of genome dynamics in plants. The last decade witnessed new applications for advanced molecular techniques in the area of breeding,with low costs per sample. These, especially, include next-generation sequencing technologies like reduced representation genome sequencing (genotyping by sequencing, restriction site-associated DNA). These enabled researchers to develop new markers, such as simple sequence repeat and single- nucleotide polymorphism, for expanding the qualitative and quantitative information onpopulation dynamics. Thus, the knowledge acquired from novel technologies is a valuable asset for the breeding process and to better understand the population dynamics, their properties, and analysis methods.

QTL mapping pod dehiscence resistance in soybean (Glycine max L. Merr.) using specific-locus amplified fragment sequencing

KEY MESSAGE

We constructed a high-density genetic linkage map comprising 4,593 SLAF markers using specific-locus amplified fragment sequencing and identified six quantitative trait loci for pod dehiscence resistance in soybean. Pod dehiscence is necessary for propagation in wild soybean (Glycine soja). It is a major component causing yield losses in cultivated soybean, however, and thus, cultivated soybean varieties have been artificially selected for resistance to pod dehiscence. Detecting quantitative trait loci (QTLs) related to pod dehiscence is required for molecular marker-assisted selection for breeding new varieties with pod dehiscence resistance. In this study, we constructed a high-density genetic linkage map using 260 recombinant inbred lines derived from the cultivars of Heihe 43 (pod-indehiscent) (ZDD24325) and Heihe 18 (pod-dehiscent) (ZDD23620). The map contained 4953 SLAF markers spanning 1478.86 cM on 20 linkage groups with an average distance between adjacent markers of 0.53 cM. In total, six novel QTLs related to pod dehiscence were mapped using inclusive composite interval mapping, explaining 7.22-24.44% of the phenotypic variance across 3 years, including three stable QTLs (qPD01, qPD05-1 and qPD08-1), that had been validated by developing CAPS/dCAPS markers. Based on the SNP/Indel and significant differential expression analyses of two parents, seven genes were selected as candidate genes for future study. The high-density map, three stable QTLs and their molecular markers will be helpful for map-based cloning of pod dehiscence resistance genes and marker-assisted selection of pod dehiscence resistance in soybean breeding.

High-density genetic map construction and QTL mapping of first flower node in pepper (Capsicum annuum L.)

BACKGROUND

First flower node (FFN) is an important trait for evaluating fruit earliness in pepper (Capsicum annuum L.). The trait is controlled by quantitative trait loci (QTL); however, studies have been limited on QTL mapping and genes contributing to the trait.

RESULTS

In this study, we developed a high density genetic map using specific-locus amplified fragment sequencing (SLAF-seq), a high-throughput strategy for de novo single nucleotide polymorphism discovery, based on 146 recombinant inbred lines (RILs) derived from an intraspecific cross between PM702 and FS871. The map contained 9328 SLAF markers on 12 linkage groups (LGs), and spanned a total genetic distance of 2009.69 centimorgan (cM) with an average distance of 0.22 cM. The sequencing depth for the map was 72.39-fold in the male parent, 57.04-fold in the female parent, and 15.65-fold in offspring. Using the genetic map, two major QTLs, named Ffn2.1 and Ffn2.2, identified on LG02 were strongly associated with FFN, with a phenotypic variance explanation of 28.62 and 19.56%, respectively. On the basis of the current annotation of C. annuum cv. Criollo de Morelos (CM334), 59 candidate genes were found within the Ffn2.1 and Ffn2.2 region, but only 3 of 59 genes were differentially expressed according to the RNA-seq results. Eventually we identified one gene associated with the FFN based on the function through GO, KEGG, and Swiss-prot analysis.

CONCLUSIONS

Our research showed that the construction of high-density genetic map using SLAF-seq is a valuable tool for fine QTL mapping. The map we constructed is by far the most saturated complete genetic map of pepper, and using it we conducted fine QTL mapping for the important trait, FFN. QTLs and candidate genes obtained in this study lay a good foundation for the further research on FFN-related genes and other genetic applications in pepper.

Construction of a High-Density Genetic Map and Analysis of Seed-Related Traits Using Specific Length Amplified Fragment Sequencing for Cucurbita maxima

Seed traits are agronomically important for Cucurbita breeding, but the genes controlling seed size, seed weight and seed number have not been mapped in Cucurbita maxima (C. maxima). In this study, 100 F₂ individual derived from two parental lines, "2013-12" and "9-6", were applied to construct a 3,376.87-cM genetic map containing 20 linkage groups (LGs) with an average genetic distance of 0.47 cM using a total of 8,406 specific length amplified fragment (SLAF) markers in C. maxima. Ten quantitative trait loci (QTLs) of seed width (SW), seed length (SL) and hundred-seed weight (HSW) were identified using the composite interval mapping (CIM) method. The QTLs affecting SW, SL and HSW explained a maximum of 38.6%, 28.9% and 17.2% of the phenotypic variation and were detected in LG6, LG6 and LG17, respectively. To validate these results, an additional 150 F₂ individuals were used for QTL mapping of SW and SL with cleaved amplified polymorphic sequence (CAPS) markers. We found that two major QTLs, SL6-1 and SW6-1, could be detected in both SLAF-seq and CAPS markers in an overlapped region. Based on gene annotation and non-synonymous single-nucleotide polymorphisms (SNPs) in the major SWand SL-associated regions, we found that two genes encoding a VQ motif and an E3 ubiquitin-protein ligase may be candidate genes influencing SL, while an F-box and leucinerich repeat (LRR) domain-containing protein is the potential regulator for SW in C. maxima. This study provides the first high-density linkage map of C. maxima using SNPs developed by SLAF-seq technology, which is a powerful tool for associated mapping of important agronomic traits, map-based gene cloning and marker-assisted selection (MAS)-based breeding in C. maxima.

Fine mapping and candidate gene screening of the downy mildew resistance gene RPF1 in Spinach

A SLAF-BSA approach was used to locate the RPF1 locus. The three most likely candidate genes were identified which provide a basic for cloning the resistance gene at the RPF1 locus. Spinach downy mildew is a globally devastating oomycete disease. The use of downy mildew resistance genes constitutes the most effective approach for disease management. Hence, the objective of the present study was to fine map the first-reported resistance locus RPF1. The resistance allele at this resistance locus was effective against races 1-7, 9, 11, 13, and 15 of Peronospora farinosa f. sp. spinaciae (P. effusa). The approach fine mapped RPF1 using specific-locus amplified fragment sequencing (SLAF-Seq) technology combined with bulked segregant analysis. A 1.72 Mb region localized on chromosome 3 was found to contain RPF1 based on association analysis. After screening recombinants with the SLAF markers within the region, the region was narrowed down to 0.89 Mb. Within this region, 14 R genes were identified based on the annotation information. To identify the genes involved in resistance, resequencing of two resistant inbred lines (12S2 and 12S3) and three susceptible inbred lines (12S1, 12S4, and 10S2) was performed. The three most likely candidate genes were identified via amino acid sequence analysis and conserved domain analysis between resistant and susceptible inbred lines. These included Spo12729, encoding a receptor-like protein, and Spo12784 and Spo12903, encoding a nucleotide-binding site and leucine-rich repeat domains. Additionally, based on the sequence variation in the three genes between the resistant and susceptible lines, molecular markers were developed for marker-assisted selection. The results could be valuable in cloning the RPF1 alleles and improving our understanding of the interaction between the host and pathogen.