A high-density SNP Map of sunflower derived from RAD-sequencing facilitating fine-mapping of the rust resistance gene R12

Construction of a high-density genetic map based on large-scale marker development in Coix lacryma-jobi L. using specific-locus amplified fragment sequencing (slaf-seq)

Coix lacryma-jobi L. is one of the most economically and medicinally important corns. This study constructed a high-density genetic linkage map of C. lacryma-jobi based on a cross between the parents 'Qianyi No. 2' × 'Wenyi No. 2' and their F2 progeny through high-throughput sequencing and the construction of a specific-locus amplified fragment (SLAF) library. After pre-processing, 325.49 GB of raw data containing 1628 M reads were obtained. A total of 22,944 high-quality SLAFs were identified, among which 3952 SLAFs and 3646 polymorphic markers met the requirements for the construction of a genetic linkage map. The integrated map contained 3605 high-quality SLAFs, which were grouped into ten genetic linkage groups. The total length of the map was 1620.39 cM, with an average distance of 0.45 cM and an average of 360.5 markers per linkage group. This report presents the first high-density genetic map of C. lacryma-jobi. This map was constructed using an F2 population and SLAF-seq approach, which allows the development of a large number of polymorphic markers in a short period. These results provide a platform for precise gene/quantitative trait locus (QTL) mapping, map-based gene separation, and molecular breeding in C. lacryma-jobi. They also help identify a target gene for tracking, splitting quantitative traits, and estimating the phenotypic effects of each QTL for QTL mapping. They are of great significance for improving the efficiency of discovering and utilizing excellent gene resources of C. lacryma-jobi.

Whole-genome sequencing enables molecular dissection and candidate gene identification of the rust resistance gene R12 in sunflower (Helianthus annuus L.)

KEY MESSAGE

We finely mapped the rust resistance gene R₁₂ to a 0.1248-cM region, identified a potential R₁₂ candidate gene in the XRQ reference genome, and developed three diagnostic SNP markers for R₁₂. Rust is a devastating disease in sunflower that is damaging to the sunflower production globally. Identification and utilization of host-plant resistance are proven to be preferable means for disease control. The rust resistance gene R₁₂ with broad-spectrum specificity to rust was previously localized to a 2.4 Mb region on sunflower chromosome 11. To understand the molecular mechanism of resistance, we conducted whole-genome sequencing of RHA 464 (R₁₂ donor line) and reference genome-based fine mapping of the gene R₁₂. Overall, the 213 markers including 186 SNPs and 27 SSRs' were identified from RHA 464 sequences and used to survey polymorphisms between the parents HA 89 and RHA 464. Saturation mapping identified 26 new markers positioned in the R₁₂ region, and fine mapping with a large population of 2004 individuals positioned R₁₂ at a genetic distance of 0.1248 cM flanked by SNP markers C11_150451336 and S11_189205190. One gene, HanXRQChr11g0348661, with a defense-related NB-ARC-LRR domain, was identified in the XRQr1.0 genome assembly in the R₁₂ region; it is predicted to be a potential R₁₂ candidate gene. Comparative analysis clearly distinguished R₁₂ from the rust R₁₄ gene located in the vicinity of the R₁₂ gene on chromosome 11. Three diagnostic SNP markers, C11_147181749, C11_147312085, and C11_149085167, specific for R₁₂ were developed in the current study, facilitating more accurate and efficient selection in sunflower rust resistance breeding. The current study provides a new genetic resource and starting point for cloning R₁₂ in the future.

Introgression and targeting of the Pl37 and Pl38 genes for downy mildew resistance from wild Helianthus annuus and H. praecox into cultivated sunflower (Helianthus annuus L.)

Two new downy mildew resistance genes, Pl₃₇ and Pl₃₈, were introgressed from wild sunflower species into cultivated sunflower and mapped to sunflower chromosomes 4 and 2, respectively Downy mildew (DM), caused by the oomycete pathogen Plasmopara halstedii (Farl.) Berl. & de Toni, is known as the most prevalent disease occurring in global sunflower production areas, especially in North America and Europe. In this study, we report the introgression and molecular mapping of two new DM resistance genes from wild sunflower species, Helianthus annuus and H. praecox, into cultivated sunflower. Two mapping populations were developed from the crosses of HA 89/H. annuus PI 435417 (Pop1) and CMS HA 89/H. praecox PRA-417 (Pop2). The phenotypic evaluation of DM resistance/susceptibility was conducted in the BC₁F₂-derived BC₁F₃ populations using P. halstedii race 734. The BC₁F₂ segregating Pop1 was genotyped using an Optimal GBS AgriSeq™ Panel consisting of 768 mapped SNP markers, while the BC₁F₂ segregating Pop2 was genotyped using a genotyping-by-sequencing approach. Linkage analysis and subsequent saturation mapping placed the DM resistance gene, designated Pl₃₇, derived from H. annuus PI 435417 in a 1.6 cM genetic interval on sunflower chromosome 4. Pl₃₇ co-segregated with SNP markers SPB0003 and C4_5738736. Similarly, linkage analysis and subsequent saturation mapping placed the DM resistance gene, designated Pl₃₈, derived from H. praecox PRA-417 in a 0.8 cM genetic interval on sunflower chromosome 2. Pl₃₈ co-segregated with seven SNP markers. Multi-pathotype tests revealed that lines with Pl₃₇ or Pl₃₈ are immune to the most prevalent and virulent P. halstedii races tested. Two germplasm lines, HA-DM15 with Pl₃₇ and HA-DM16 with Pl₃₈, were developed for use in sunflower DM-resistance breeding.

Genetic mapping of a pollinator preference trait: Nectar volume in sunflower (Helianthus annuus L.)

Although high pollinator visitation is crucial to ensure the yields of pollinator-dependent crops, the quantitative trait loci (QTL) controlling nectar volume in sunflower (Helianthus annuus L.), a pollinator preference trait, have yet to be identified. To address this, a recombinant inbred line mapping population, derived from lines with contrasting nectar volume, was used to identify loci responsible for the phenotype. As a result, linkage mapping and QTL analysis discovered major loci on chromosomes 2 and 16 that are associated with variation in nectar volume in sunflower. Increased nectar volume is also associated with increased sugars and total energy available per floret. The regions on chromosomes 2 and 16 associated with the nectar phenotype exhibit indications of chromosome structural variation, such that the phenotype is associated with rearrangements affecting regions containing hundreds of genes. Candidate genes underlying QTL on chromosomes 9 and 16 are homologous to genes with nectary function in Arabidopsis. These results have implications for sunflower breeding, to enhance pollination efficiency in sunflower, as well as current and future studies on sunflower evolution.

A Quantitative Genetic Study of Sclerotinia Head Rot Resistance Introgressed from the Wild Perennial Helianthus maximiliani into Cultivated Sunflower (Helianthus annuus L.)

Sclerotinia head rot (HR), caused by Sclerotinia sclerotiorum, is an economically important disease of sunflower with known detrimental effects on yield and quality in humid climates worldwide. The objective of this study was to gain insight into the genetic architecture of HR resistance from a sunflower line HR21 harboring HR resistance introgressed from the wild perennial Helianthus maximiliani. An F2 population derived from the cross of HA 234 (susceptible-line)/HR21 (resistant-line) was evaluated for HR resistance at two locations during 2019−2020. Highly significant genetic variations (p < 0.001) were observed for HR disease incidence (DI) and disease severity (DS) in both individual and combined analyses. Broad sense heritability (H2) estimates across environments for DI and DS were 0.51 and 0.62, respectively. A high-density genetic map of 1420.287 cM was constructed with 6315 SNP/InDel markers developed using genotype-by-sequencing technology. A total of 16 genomic regions on eight sunflower chromosomes, 1, 2, 10, 12, 13, 14, 16 and 17 were associated with HR resistance, each explaining between 3.97 to 16.67% of the phenotypic variance for HR resistance. Eleven of these QTL had resistance alleles from the HR21 parent. Molecular markers flanking the QTL will facilitate marker-assisted selection breeding for HR resistance in sunflower.

Genomic Insights Into Sclerotinia Basal Stalk Rot Resistance Introgressed From Wild Helianthus praecox Into Cultivated Sunflower (Helianthus annuus L.)

Crop wild relatives of the cultivated sunflower (Helianthus annuus L.) are a valuable resource for its sustainable production. Helianthus praecox ssp. runyonii is a wild sunflower known for its resistance against diseases caused by the fungus, Sclerotinia sclerotiorum (Lib.) de Bary, which infects over 400 broadleaf hosts including many important food crops. The objective of this research was to dissect the Sclerotinia basal stalk rot (BSR) resistance introgressed from H. praecox ssp. runyonii into cultivated sunflower. An advanced backcross quantitative trait loci (AB-QTL) mapping population was developed from the cross of a H. praecox accession with cultivated sunflower lines. The AB-QTL population was evaluated for BSR resistance in the field during the summers of 2017-2018 and in the greenhouse in the spring of 2018. Highly significant genetic variations (p < 0.001) were observed for the BSR disease in the field and greenhouse with a moderately high broad-sense heritability (H ²) ranging from 0.66 to 0.73. Genotyping-by-sequencing approach was used to genotype the parents and the progeny lines of the AB-QTL population. A genetic linkage map spanning 1,802.95 cM was constructed using 1,755 single nucleotide polymorphism (SNP) markers mapped on 17 sunflower chromosomes. A total of 19 BSR resistance QTL were detected on nine sunflower chromosomes, each explaining 2.21%-16.99% of the phenotypic variance for resistance in the AB-QTL population. Sixteen of the 19 QTL had alleles conferring increased BSR resistance derived from the H. praecox parent. SNP markers flanking the identified QTL will facilitate marker-assisted breeding to combat the disease in sunflower.

Genetic mapping of loci involved in oil tocopherol composition control in Russian sunflower (Helianthus annuus L.) lines

Tocopherols are antioxidants that preserve oil lipids against oxidation and serve as a natural source of vitamin E in the human diet. Compared with other major oilseeds like rapeseed and soybean, sunflower (Helianthus annuus L.) exhibits low phenotypic diversity of tocopherol composition, both in wild and cultivated accessions from germplasm collections. Two major mutations that alter tocopherol composition were identified in genetic collections, and several studies suggested additional loci controlling tocopherol composition, with their expression possibly depending on the genetic background. In the present study, we performed QTL mapping of tocopherol composition in two independent F2 crosses between lines with contrasting tocopherol composition from the Pustovoit All-Russia Research Institute of Oil Crops (VNIIMK) collection. We used genotyping-bysequencing (GBS) to construct single nucleotide polymorphism-based genetic maps, and performed QTL mapping using quantitative and qualitative encoding for phenotypic traits. Our results support the notion that the tocopherol composition in the assessed crosses is controlled by two loci. We additionally selected and validated two single nucleotide polymorphism markers for each cross which could be used for marker-assisted selection.

Inflorescence Transcriptome Sequencing and Development of New EST-SSR Markers in Common Buckwheat (Fagopyrum esculentum)

Common buckwheat (Fagopyrum esculentum M.) is known for its adaptability, good nutrition, and medicinal and health care value. However, genetic studies of buckwheat have been hindered by limited genomic resources and genetic markers. In this study, Illumina HiSeq 4000 high-throughput sequencing technology was used to sequence the transcriptome of green-flower common buckwheat (Gr) with coarse pedicels and white-flower Ukrainian daliqiao (UD) with fine pedicels. A total of 118,448 unigenes were obtained, with an average length of 1248 bp and an N50 of 1850 bp. A total of 39,432 differentially expressed genes (DEGs) were identified, and the DEGs of the porphyrins and chlorophyll metabolic pathway had significantly upregulated expression in Gr. Then, a total of 17,579 sequences containing SSR loci were detected, and 20,756 EST-SSR loci were found. The distribution frequency of EST-SSR in the transcriptome was 17.52%, and the average distribution density was 8.21 kb. A total of 224 pairs of primers were randomly selected for synthesis; 35 varieties of common buckwheat and 13 varieties of Tartary buckwheat were verified through these primers. The clustering results well verified the previous conclusion that common buckwheat and Tartary buckwheat had a distant genetic relationship. The EST-SSR markers identified and developed in this study will be helpful to enrich the transcriptome information and marker-assisted selection breeding of buckwheat.

Characterization and mapping of a downy mildew resistance gene, Pl36, in sunflower (Helianthus annuus L.)

Downy mildew (DM) is one of the most serious diseases in sunflower-growing regions worldwide, often significantly reducing sunflower yields. The causal agent of sunflower DM, the oomycete pathogen Plasmopara halstedii, is highly virulent and aggressive. Studying regional disease spread and virulence evolution in the DM pathogen population is important for the development of new sunflower inbred lines with resistance to the existing DM pathogen. The sunflower line 803-1, as one of nine international differential hosts, has been used in the identification of P. halstedii virulent pathotypes in sunflower since 2000. The DM resistance gene in 803-1 was temporally designated Pl5 + based on allelic analysis but has not been molecularly characterized. In the present study, bulked segregant analysis and genetic mapping confirmed the presence of the Pl gene within a large gene cluster on sunflower chromosome 13 in 803-1, as previously reported. Subsequent saturation mapping in the gene target region with single nucleotide polymorphism (SNP) markers placed this gene at an interval of 3.4 Mb in the XRQ reference genome assembly, a location different from that of Pl5. Therefore, the Pl gene in 803-1 was re-designated Pl36 because it is not allelic with Pl5. Four SNP markers co-segregated with Pl36, and SNP SFW05743 was 1.1 cM proximal to Pl36. The relationship of eight Pl genes in the cluster is discussed based on their origin, map position, and specificity of resistance/susceptibility to DM infection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01280-1.

Construction of a High-Density Genetic Linkage Map for the Mapping of QTL Associated with Growth-Related Traits in Sea Cucumber (Apostichopus japonicus)

Simple Summary

Slow growth and germplasm degradation have restricted the sustainable commercial development of the sea cucumber industry. To analyze the genetic mechanism of growth traits of sea cucumbers, we constructed a high-density genetic linkage map based on single nucleotide polymorphism (SNP) molecular markers and performed a quantitative trait loci (QTL) mapping analysis. We annotated a critical candidate gene related to growth traits and explored mRNA expression levels. The results showed that the gene was significantly highly expressed during the larval developmental stages. These results can be used to genetically improve the growth traits of sea cucumbers.

Abstract

Genetic linkage maps have become an indispensable tool for genetics and genomics research. Sea cucumber (Apostichopus japonicus), which is an economically important mariculture species in Asia, is an edible echinoderm with medicinal properties. In this study, the first SNP-based high-density genetic linkage map was constructed by sequencing 132 A. japonicus individuals (2 parents and 130 offspring) according to a genotyping-by-sequencing (GBS) method. The consensus map was 3181.54 cM long, with an average genetic distance of 0.52 cM. A total of 6144 SNPs were assigned to 22 linkage groups (LGs). A Pearson analysis and QTL mapping revealed the correlations among body weight, body length, and papillae number. An important growth-related candidate gene, protein still life, isoforms C/SIF type 2 (sif), was identified in LG18. The gene was significantly highly expressed during the larval developmental stages. Its encoded protein reportedly functions as a guanine nucleotide exchange factor. These results would facilitate the genetic analysis of growth traits and provide valuable genomic resources for the selection and breeding of new varieties of sea cucumbers with excellent production traits.

RAD-R scripts: R pipeline for RAD-seq from FASTQ files to linkage maps construction and run R/QTL, operating only at copying and pasting scripts into R console

Coupled with the reduction in sequencing costs, the number of RAD-seq analysis have been surging, generating vast genetic knowledge in relation with many crops. Specialized platforms might be intimidating to non-expert users and difficult to implement on each computer despite the growing interest in the usage of the dataset obtained by high-throughput sequencing. Therefore, RAD-R scripts were developed on Windows10 for RAD-seq analysis, allowing users who are not familiar with bioinformatics to easily analyze big sequence data. These RAD-R scripts that run a flow from raw sequence reads of F₂ population for the self-fertilization plants to the linkage map construction as well as the QTL analysis can be also useful to many users with limited experience due to the simplicity of copying Excel cells into the R console. During the comparison of linkage maps constructed by RAD-R scripts and Stacks, RAD-R scripts were shown to construct the linkage map with less missing genotype data and a shorter total genetic distance. QTL analysis results can be easily obtained by selecting the reliable genotype data that is visually inferred to be appropriate for error correction from the genotype data files created by RAD-R scripts.

Discovery and mapping of two new rust resistance genes, R17 and R18, in sunflower using genotyping by sequencing

Discovery of two rust resistance genes, R₁₇ and R₁₈, from the sunflower lines introduced from South Africa and genetic mapping of them to sunflower chromosome 13. Rust, caused by the fungus Puccinia helianthi Schw., is one of the most serious diseases of sunflower in the world. The rapid changes that occur in the virulence characteristics of pathogen populations present a continuous threat to the effectiveness of existing rust-resistant hybrids. Thus, there is a continued need for the characterization of genetically diverse sources of rust resistance. In this study, we report to identify two new rust resistance genes, R₁₇ and R₁₈, from the sunflower lines, KP193 and KP199, introduced from South Africa. The inheritance of rust resistance was investigated in both lines using two mapping populations developed by crossing the resistant plants selected from KP193 and KP199 with a common susceptible parent HA 89. The F₂ populations were first genotyped using genotyping by sequencing for mapping of the rust genes and further saturated with markers in the target region. Molecular mapping positioned the two genes at the lower end of sunflower chromosome 13 within a large gene cluster. Two co-segregating SNP markers, SFW01497 and SFW08875, were distal to R₁₇ at a 1.9 cM genetic distance, and a cluster of five co-segregating SNPs was proximal to R₁₇ at 0.7 cM. R₁₈ co-segregated with the SNP marker SFW04317 and was proximal to two cosegregating SNPs, SFW01497 and SFW05453, at 1.9 cM. These maps provide markers for stacking R₁₇ or R₁₈ with other broadly effective rust resistance genes to extend the durability of rust resistance. The relationship of the six rust resistance genes in the cluster was discussed.

Unraveling the Sclerotinia Basal Stalk Rot Resistance Derived From Wild Helianthus argophyllus Using a High-Density Single Nucleotide Polymorphism Linkage Map

Basal stalk rot (BSR), caused by the fungus Sclerotinia sclerotiorum, is a serious disease of sunflower (Helianthus annuus L.) in the humid temperate growing areas of the world. BSR resistance is quantitative and conditioned by multiple genes. Our objective was to dissect the BSR resistance introduced from the wild annual species Helianthus argophyllus using a quantitative trait loci (QTL) mapping approach. An advanced backcross population (AB-QTL) with 134 lines derived from the cross of HA 89 with a H. argophyllus Torr. and Gray accession, PI 494573, was evaluated for BSR resistance in three field and one greenhouse growing seasons of 2017-2019. Highly significant genetic variations (p < 0.001) were observed for BSR disease incidence (DI) in all field screening tests and disease rating and area under the disease progress curve in the greenhouse. The AB-QTL population and its parental lines were genotyped using the genotyping-by-sequencing method. A genetic linkage map spanning 2,045.14 cM was constructed using 3,110 SNP markers mapped on 17 sunflower chromosomes. A total of 21 QTL associated with BSR resistance were detected on 11 chromosomes, each explaining a phenotypic variation ranging from 4.5 to 22.6%. Of the 21 QTL, eight were detected for BSR DI measured in the field, seven were detected for traits measured in the greenhouse, and six were detected from both field and greenhouse tests. Thirteen of the 21 QTL had favorable alleles from the H. argophyllus parent conferring increased BSR resistance.

A Unique Cytoplasmic-Nuclear Interaction in Sunflower (Helianthus annuus L.) Causing Reduced-Vigor Plants and the Genetics of Vigor Restoration

Wild Helianthus species are an important genetic resource for sunflower improvement, but sometimes there are adverse interactions between the wild and cultivated sunflowers. This study reports the inheritance of reduced vigor and its restoration resulting from an interaction of perennial Helianthus cytoplasms with nuclear genes of cultivated sunflower lines. The large number of vigor restoration (V) genes identified in cultivated lines are all located at the same locus, designated V₁ , suggesting a common origin of these genes. Additional V genes derived from the wild perennial species H. giganteus L. and H. hirsutus Raf. are located at a different locus than V₁ , designated V₂ . A major difference between the wild annual Helianthus cytoplasms and perennial cytoplasms is the lack of the vigor-reducing cytoplasms, but surprisingly V genes were observed in wild annual H. annuus L. and H. petiolaris Nutt. which were at the same locus as V₁ . A common vigor-reducing cytoplasmic effect of the perennial Helianthus species and the existence of a common vigor restoration V gene in most perennial Helianthus species could be explained as a result of vigor selection during Helianthus speciation. V₁ was mapped on linkage group (LG) 7 of the sunflower genome, using an F₂ population derived from MOL-RV/HA 821. V₁ co-segregated with an InDel marker ZVG31, with three single-nucleotide polymorphism (SNP) markers, SFW01024, SFW07230, and SFW00604, located above it on the map at a genetic distance of 0.8 cM, and another SNP marker, SFW08671, below it at a distance of 0.4 cM. The physical distance between the two closest flanking SNP markers corresponds to 0.56 and 1.37 Mb on the HA 412-HO and XRQ assemblies, respectively. The tightly linked markers will help select normal vigor progenies when using perennial Helianthus cytoplasms in a breeding program, which will also provide a basis for studying the mechanism of the cytonuclear interaction, and the speciation of annual and perennial Helianthus species.

Inheritance and molecular mapping of powdery mildew (Golovinomyces orontii) resistance gene(s) in sunflower (Helianthus annuus L.)

Sources of resistance to powdery mildew incited by Golovinomyces orontii have been identified in wild sunflowers and few exotic lines. The present investigation has been undertaken to study the inheritance of powdery mildew resistance and to map the quantitative trait loci (QTLs) governing resistance to powdery mildew in a multiple disease resistance line, TX16R (PI 642072). The inheritance was observed as a continuous distribution in a set of 264 F₂ population and 93 recombinant inbred lines (RILs) of a cross between a highly susceptible accession PS 2023 and TX16R. Screening of the two population sets was done with 484 sunflower-specific SSR primers of which 175 primers showed polymorphism between the parents. Based on the phenotyping and genotyping data, the linkage map was constructed with 93 RILs. The map spanned 1200 cM and included 64 markers distributed along the 17 sunflower chromosomes in the haploid set. Quantitative trait loci (QTL) analysis identified three genomic regions for resistance to powdery mildew, two of which mapped on chromosome 10 and one on chromosome 5. This is the first report on mapping of powdery mildew resistance in sunflower and paves the way in fine mapping and introgression of resistance for powdery mildew in sunflower through marker-assisted breeding.

Genetic Dissection of Phomopsis Stem Canker Resistance in Cultivated Sunflower Using High Density SNP Linkage Map

Phomopsis stem canker (PSC) caused by Diaporthe helianthi is increasingly becoming a global threat for sunflower production. In this study, the genetic basis of PSC resistance was investigated in a recombinant inbred line (RIL) population developed from a cross between HA 89 (susceptible) and HA-R3 (resistant). The RIL population was evaluated for PSC disease incidence (DI) in seven screening trials at multiple locations during 2016–2018. The distribution of PSC DI in the RIL population was continuous, confirming a polygenic inheritance of the trait. A moderately high broad-sense heritability (H², 0.76) was estimated for the trait across environments. In the combined analysis, both the genotype and the genotype × environment interactions were highly significant. A linkage map spanning 1505.33 cM was constructed using genotyping-by-sequencing derived markers. Marker–trait association analysis identified a total of 15 quantitative trait loci (QTL) associated with PSC resistance on 11 sunflower chromosomes, each explaining between 5.24 and 17.39% of the phenotypic variation. PSC resistance QTL were detected in two genomic regions each on chromosomes 3, 5, 13, and 17, while one QTL each was detected in the remaining seven chromosomes. Tightly linked single nucleotide polymorphism (SNP) markers flanking the PSC resistance QTL will facilitate marker-assisted selection in PSC resistance sunflower breeding.

Marker-Assisted Gene Pyramiding and the Reliability of Using SNP Markers Located in the Recombination Suppressed Regions of Sunflower (Helianthus annuus L.)

Rust caused by the fungus Puccinia helianthi and downy mildew (DM) caused by the obligate pathogen Plasmopara halstedii are two of the most globally important sunflower diseases. Resistance to rust and DM is controlled by race-specific single dominant genes. The present study aimed at pyramiding rust resistance genes combined with a DM resistance gene, using molecular markers. Four rust resistant lines, HA-R3 (carrying the R4 gene), HA-R2 (R5), HA-R8 (R15), and RHA 397 (R13b), were each crossed with a common line, RHA 464, carrying a rust gene R12 and a DM gene PlArg. An additional cross was made between HA-R8 and RHA 397. Co-dominant simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers linked to the target genes were used to discriminate between homozygotes and heterozygotes in F2 populations. Five pyramids with different combinations of rust resistance genes were selected in the homozygous condition through marker-assisted selection, and three of them were combined with a DM resistance gene PlArg: R4/R12/PlArg, R5/R12/PlArg, R13b/R12/PlArg, R15/R12, and R13b/R15. The pyramiding lines with the stacking of two rust and one DM genes were resistant to all known races of North American sunflower rust and all known races of the pathogen causing DM, potentially providing multiple and durable resistance to both rust and DM. A cluster of 12 SNP markers spanning a region of 34.5 Mb on chromosome 1, which co-segregate with PlArg, were tested in four populations. Use of those markers, located in a recombination suppressed region in marker selection, is discussed.

Construction of the first high-density genetic linkage map and identification of seed yield-related QTLs and candidate genes in Elymus sibiricus, an important forage grass in Qinghai-Tibet Plateau

Background

Elymus sibiricus is an ecologically and economically important perennial, self-pollinated, and allotetraploid (StStHH) grass, widely used for forage production and animal husbandry in Western and Northern China. However, it has low seed yield mainly caused by seed shattering, which makes seed production difficult for this species. The goals of this study were to construct the high-density genetic linkage map, and to identify QTLs and candidate genes for seed-yield related traits.

Results

An F₂ mapping population of 200 individuals was developed from a cross between single genotype from “Y1005” and “ZhN06”. Specific-locus amplified fragment sequencing (SLAF-seq) was applied to construct the first genetic linkage map. The final genetic map included 1971 markers on the 14 linkage groups (LGs) and was 1866.35 cM in total. The length of each linkage group varied from 87.67 cM (LG7) to 183.45 cM (LG1), with an average distance of 1.66 cM between adjacent markers. The marker sequences of E. sibiricus were compared to two grass genomes and showed 1556 (79%) markers mapped to wheat, 1380 (70%) to barley. Phenotypic data of eight seed-related traits (2016–2018) were used for QTL identification. A total of 29 QTLs were detected for eight seed-related traits on 14 linkage groups, of which 16 QTLs could be consistently detected for two or three years. A total of 6 QTLs were associated with seed shattering. Based on annotation with wheat and barley genome and transcriptome data of abscission zone in E. sibiricus, we identified 30 candidate genes for seed shattering, of which 15, 7, 6 and 2 genes were involved in plant hormone signal transcription, transcription factor, hydrolase activity and lignin biosynthetic pathway, respectively.

Conclusion

This study constructed the first high-density genetic linkage map and identified QTLs and candidate genes for seed-related traits in E. sibiricus. Results of this study will not only serve as genome-wide resources for gene/QTL fine mapping, but also provide a genetic framework for anchoring sequence scaffolds on chromosomes in future genome sequence assembly of E. sibiricus.

Diversification of the downy mildew resistance gene pool by introgression of a new gene, Pl35, from wild Helianthus argophyllus into oilseed and confection sunflowers (Helianthus annuus L.)

We have mapped a new downy mildew resistance gene, Pl₃₅, derived from wild Helianthus argophyllus to sunflower linkage group 1. New germplasms incorporating the Pl₃₅ gene were developed for both oilseed and confection sunflower Sunflower downy mildew (DM), caused by the oomycete pathogen Plasmopara halstedii, is an economically important and widespread sunflower disease worldwide. Non-race-specific resistance is not available in sunflower, and breeding for DM resistance relies on race-specific resistance to control this disease. The discovery of the novel DM resistance genes is a long-term task due to the highly virulent and aggressive nature of the P. halstedii pathogen, which reduces the effectiveness of resistance genes. The objectives of this study were to: (1) transfer DM resistance from a wild sunflower species Helianthus argophyllus (PI 494576) into cultivated sunflowers; (2) map the resistance gene; and (3) develop diagnostic single-nucleotide polymorphism (SNP) markers for efficient targeting of the gene in breeding programs. The H. argophyllus accession PI 494576 previously identified with resistance to the most virulent P. halstedii race 777 was crossed with oilseed and confection sunflower in 2012. Molecular mapping using the BC₂F₂ and BC₂F₃ populations derived from the cross CONFSCLB1/PI 494576 located a new resistance gene Pl₃₅ on linkage group 1 of the sunflower genome. The new gene Pl₃₅ was successfully transferred from PI 494576 into cultivated sunflowers. SNP markers flanking Pl₃₅ were surveyed in a validation panel of 548 diversified sunflower lines collected globally. Eleven SNP markers were found to be diagnostic for Pl₃₅ SNP alleles, with four co-segregating with Pl₃₅. The developed oilseed and confection germplasms with diagnostic SNP markers for Pl₃₅ will be very useful resources for breeding of DM resistance in sunflower.

Linkage Mapping and Genome-Wide Association Studies of the Rf Gene Cluster in Sunflower (Helianthus annuus L.) and Their Distribution in World Sunflower Collections

Commercial hybrid seed production in sunflower currently relies on a single cytoplasmic male sterility (CMS) source, PET1 and the major fertility restoration gene, Rf1, leaving the crop highly vulnerable to issues with genetic bottlenecks. Therefore, having multiple CMS/Rf systems is important for sustainable sunflower production. Here, we report the identification of a new fertility restoration gene, Rf7, which is tightly linked to a new downy mildew (DM) resistance gene, Pl₃₄ , in the USDA sunflower inbred line, RHA 428. The Rf7 gene was genetically mapped to an interval of 0.6 cM on the lower end of linkage group (LG) 13, while Pl₃₄ was mapped 2.1 cM proximal to the Rf7. Both the genes are located in a cluster of Rf and Pl genes. To gain further insights into the distribution of Rf genes in the sunflower breeding lines, we used a genome-wide association study (GWAS) approach to identify markers associated with the fertility restoration trait in a panel of 333 sunflower lines genotyped with 8,723 single nucleotide polymorphism (SNP) markers. Twenty-four SNP markers on the lower end of LG13 spanning a genomic region of 2.47 cM were significantly associated with the trait. The significant markers were surveyed in a world collection panel of 548 sunflower lines and validated to be associated with the Rf1 gene. The SNP haplotypes for the Rf1 gene are different from Rf5 and the Rf7gene located in the Rf gene cluster on LG13. The SNP and SSR markers tightly flanking the Rf7 gene and the Pl₃₄ gene would benefit the sunflower breeders in facilitating marker assisted selection (MAS) of Rf and Pl genes.

Linkage Mapping and Genome-Wide Association Studies of the Rf Gene Cluster in Sunflower (Helianthus annuus L.) and Their Distribution in World Sunflower Collections

Commercial hybrid seed production in sunflower currently relies on a single cytoplasmic male sterility (CMS) source, PET1 and the major fertility restoration gene, Rf1, leaving the crop highly vulnerable to issues with genetic bottlenecks. Therefore, having multiple CMS/Rf systems is important for sustainable sunflower production. Here, we report the identification of a new fertility restoration gene, Rf7, which is tightly linked to a new downy mildew (DM) resistance gene, Pl₃₄ , in the USDA sunflower inbred line, RHA 428. The Rf7 gene was genetically mapped to an interval of 0.6 cM on the lower end of linkage group (LG) 13, while Pl₃₄ was mapped 2.1 cM proximal to the Rf7. Both the genes are located in a cluster of Rf and Pl genes. To gain further insights into the distribution of Rf genes in the sunflower breeding lines, we used a genome-wide association study (GWAS) approach to identify markers associated with the fertility restoration trait in a panel of 333 sunflower lines genotyped with 8,723 single nucleotide polymorphism (SNP) markers. Twenty-four SNP markers on the lower end of LG13 spanning a genomic region of 2.47 cM were significantly associated with the trait. The significant markers were surveyed in a world collection panel of 548 sunflower lines and validated to be associated with the Rf1 gene. The SNP haplotypes for the Rf1 gene are different from Rf5 and the Rf7gene located in the Rf gene cluster on LG13. The SNP and SSR markers tightly flanking the Rf7 gene and the Pl₃₄ gene would benefit the sunflower breeders in facilitating marker assisted selection (MAS) of Rf and Pl genes.

Construction of a high-density genetic linkage map and QTL mapping of oleic acid content and three agronomic traits in sunflower (Helianthus annuus L.) using specific-locus amplified fragment sequencing (SLAF-seq)

High-density genetic linkage maps are particularly important for quantitative trait loci (QTL) mapping, genome assembly, and marker-assisted selection (MAS) in plants. In this study, a high-density genetic linkage map of sunflower (Helianthus annuus L.) was constructed using an F₂ population generated from a cross between Helianthus annuus L. '86-1' and 'L-1-OL-1' via specific-locus amplified fragment sequencing (SLAF-seq). After sequence preprocessing, 530.50 M reads (105.60 Gb) were obtained that contained a total of 343,197 SLAFs, of which 39,589 were polymorphic. Of the polymorphic SLAFs, 6,136 were organized into a linkage map consisting of 17 linkage groups (LGs) spanning 2,221.86 cM, with an average genetic distance of 0.36 cM between SLAFs. Based on this high-density genetic map, QTL analysis was performed that focused on four sunflower phenotypic traits: oleic acid content (OAC), plant height (PH), head diameter (HD), and stem diameter (SD). Subsequently, for these four traits eight QTLs were detected that will likely be useful for increasing our understanding of genetic factors underlying these traits and for use in marker-assisted selection (MAS) for future sunflower breeding.

High-throughput genotyping-by-sequencing facilitates molecular tagging of a novel rust resistance gene, R 15 , in sunflower (Helianthus annuus L.)

A novel rust resistance gene, R ₁₅ , derived from the cultivated sunflower HA-R8 was assigned to linkage group 8 of the sunflower genome using a genotyping-by-sequencing approach. SNP markers closely linked to R ₁₅ were identified, facilitating marker-assisted selection of resistance genes. The rust virulence gene is co-evolving with the resistance gene in sunflower, leading to the emergence of new physiologic pathotypes. This presents a continuous threat to the sunflower crop necessitating the development of resistant sunflower hybrids providing a more efficient, durable, and environmentally friendly host plant resistance. The inbred line HA-R8 carries a gene conferring resistance to all known races of the rust pathogen in North America and can be used as a broad-spectrum resistance resource. Based on phenotypic assessments of 140 F₂ individuals derived from a cross of HA 89 with HA-R8, rust resistance in the population was found to be conferred by a single dominant gene (R ₁₅ ) originating from HA-R8. Genotypic analysis with the currently available SSR markers failed to find any association between rust resistance and any markers. Therefore, we used genotyping-by-sequencing (GBS) analysis to achieve better genomic coverage. The GBS data showed that R ₁₅ was located at the top end of linkage group (LG) 8. Saturation with 71 previously mapped SNP markers selected within this region further showed that it was located in a resistance gene cluster on LG8, and mapped to a 1.0-cM region between three co-segregating SNP makers SFW01920, SFW00128, and SFW05824 as well as the NSA_008457 SNP marker. These closely linked markers will facilitate marker-assisted selection and breeding in sunflower.

Construction of a SNP-Based High-Density Genetic Map Using Genotyping by Sequencing (GBS) and QTL Analysis of Nut Traits in Chinese Chestnut (Castanea mollissima Blume)

Chinese chestnut is a wildly distributed nut species with importantly economic value. The nut size and ripening period are mainly desired breeding objectives in Chinese chestnut. However, high-density linkage maps and quantitative trait loci (QTL) analyses related to nut traits are less than satisfactory, which hinders progress in the breeding of Chinese chestnut. Here, a single nucleotide polymorphism (SNP)-based high-density linkage map was constructed through genotyping-by-sequencing (GBS) of an F1 cross between the two widely grown Chinese chestnut cultivars 'Yanshanzaofeng' and 'Guanting No. 10'. The genetic linkage map consists of 2,620 SNP markers with a total length of 1078.06 cM in 12 linkage groups (LGs) and an average marker distance of 0.41 cM. 17 QTLs were identified for five nut traits, specifically single-nut weight (SNW), nut width (NW), nut thickness (NT), nut height (NH), and ripening period (RP), based on phenotypic data from two successive years. Of the 17 QTLs, two major QTLs, i.e., qNT-I-1 and qRP-B-1 related to the NT and RP traits, respectively, were exploited. Moreover, the data revealed one pleiotropic QTL at 23.97 cM on LG I, which might simultaneously control SNW, NT, and NW. This study provides useful benchmark information concerning high-density genetic mapping and QTLs identification related to nut size and ripening period, and will accelerate genetic improvements for nuts in the marker-assisted selection (MAS) breeding of Chinese chestnut.

Genetic and Genomic Tools to Asssist Sugar Beet Improvement: The Value of the Crop Wild Relatives

Sugar beet (Beta vulgaris L. ssp. vulgaris) is one of the most important European crops for both food and sugar production. Crop improvement has been developed to enhance productivity, sugar content or other breeder's desirable traits. The introgression of traits from Crop Wild Relatives (CWR) has been done essentially for lessening biotic stresses constraints, namely using Beta and Patellifolia species which exhibit disease resistance characteristics. Several studies have addressed crop-to-wild gene flow, yet, for breeding programs genetic variability associated with agronomically important traits remains unexplored regarding abiotic factors. To accomplish such association from phenotype-to-genotype, screening for wild relatives occurring in habitats where selective pressures are in play (i.e., populations in salt marshes for salinity tolerance; populations subjected to pathogen attacks and likely evolved resistance to pathogens) are the most appropriate streamline to identify causal genetic information. By selecting sugar beet CWR species based on genomic tools, rather than random variations, is a promising but still seldom explored route toward the development of improved crops. In this perspective, a viable streamline for sugar beet improvement is proposed through the use of different genomic tools by recurring to sugar beet CWRs and focusing on agronomic traits associated with abiotic stress tolerance. Overall, identification of genomic and epigenomic landscapes associated to adaptive ecotypes, along with the cytogenetic and habitat characterization of sugar beet CWR, will enable to identify potential hotspots for agrobiodiversity of sugar beet crop improvement toward abiotic stress tolerance.

Sunflower Hybrid Breeding: From Markers to Genomic Selection

In sunflower, molecular markers for simple traits as, e.g., fertility restoration, high oleic acid content, herbicide tolerance or resistances to Plasmopara halstedii, Puccinia helianthi, or Orobanche cumana have been successfully used in marker-assisted breeding programs for years. However, agronomically important complex quantitative traits like yield, heterosis, drought tolerance, oil content or selection for disease resistance, e.g., against Sclerotinia sclerotiorum have been challenging and will require genome-wide approaches. Plant genetic resources for sunflower are being collected and conserved worldwide that represent valuable resources to study complex traits. Sunflower association panels provide the basis for genome-wide association studies, overcoming disadvantages of biparental populations. Advances in technologies and the availability of the sunflower genome sequence made novel approaches on the whole genome level possible. Genotype-by-sequencing, and whole genome sequencing based on next generation sequencing technologies facilitated the production of large amounts of SNP markers for high density maps as well as SNP arrays and allowed genome-wide association studies and genomic selection in sunflower. Genome wide or candidate gene based association studies have been performed for traits like branching, flowering time, resistance to Sclerotinia head and stalk rot. First steps in genomic selection with regard to hybrid performance and hybrid oil content have shown that genomic selection can successfully address complex quantitative traits in sunflower and will help to speed up sunflower breeding programs in the future. To make sunflower more competitive toward other oil crops higher levels of resistance against pathogens and better yield performance are required. In addition, optimizing plant architecture toward a more complex growth type for higher plant densities has the potential to considerably increase yields per hectare. Integrative approaches combining omic technologies (genomics, transcriptomics, proteomics, metabolomics and phenomics) using bioinformatic tools will facilitate the identification of target genes and markers for complex traits and will give a better insight into the mechanisms behind the traits.

Genetic Architecture of Capitate Glandular Trichome Density in Florets of Domesticated Sunflower (Helianthus annuus L.)

Capitate glandular trichomes (CGT), one type of glandular trichomes, are most common in Asteraceae species. CGT can produce various secondary metabolites such as sesquiterpene lactones (STLs) and provide durable resistance to insect pests. In sunflower, CGT-based host resistance is effective to combat the specialist pest, sunflower moth. However, the genetic basis of CGT density is not well understood in sunflower. In this study, we identified two major QTL controlling CGT density in sunflower florets by using a F₄ mapping population derived from the cross HA 300 × RHA 464 with a genetic linkage map constructed from genotyping-by-sequencing data and composed of 2121 SNP markers. One major QTL is located on chromosome 5, which explained 11.61% of the observed phenotypic variation, and the second QTL is located on chromosome 6, which explained 14.06% of the observed phenotypic variation. The QTL effects and the association between CGT density and QTL support interval were confirmed in a validation population which included 39 sunflower inbred lines with diverse genetic backgrounds. We also identified two strong candidate genes in the QTL support intervals, and the functions of their orthologs in other plant species suggested their potential roles in regulating capitate glandular trichome density in sunflower. Our results provide valuable information to sunflower breeding community for developing host resistance to sunflower insect pests.

Revolution in Genotyping Platforms for Crop Improvement

In the past decade, the application of high-throughput sequencing to crop genotyping has given rise to novel platforms capable of genotyping tens of thousands of genome-wide DNA markers. Coupled with the decreasing costs of sequencing, this rapid increase in markers allows accelerated and highly accurate genotyping of entire crop populations and diversity sets using single nucleotide polymorphisms (SNPs). These revolutionary advances accelerate crop improvement by facilitating a more precise connection of phenotype to genotype through association studies, linkage mapping and diversity analysis. The platforms driving the advances in genotyping are array technologies and genotyping by sequencing (GBS) methods, which include both low-coverage whole genome resequencing (skim sequencing) and reduced representation sequencing (RRS) approaches. Here, we outline and compare these genotyping platforms and provide a perspective on the promising future of crop genotyping. While SNP arrays provide high quality, simple handling, and unchallenging analysis, the lower cost of RRS and the greater data volume produced by skim sequencing suggest that use of GBS will become more prevalent in crop genomics as sequencing costs decrease and data analysis becomes more streamlined. Graphical Abstract.

MiSNPDb: a web-based genomic resources of tropical ecology fruit mango (Mangifera indica L.) for phylogeography and varietal differentiation

Mango is one of the most important fruits of tropical ecological region of the world, well known for its nutritive value, aroma and taste. Its world production is >45MT worth >200 billion US dollars. Genomic resources are required for improvement in productivity and management of mango germplasm. There is no web-based genomic resources available for mango. Hence rapid and cost-effective high throughput putative marker discovery is required to develop such resources. RAD-based marker discovery can cater this urgent need till whole genome sequence of mango becomes available. Using a panel of 84 mango varieties, a total of 28.6 Gb data was generated by ddRAD-Seq approach on Illumina HiSeq 2000 platform. A total of 1.25 million SNPs were discovered. Phylogenetic tree using 749 common SNPs across these varieties revealed three major lineages which was compared with geographical locations. A web genomic resources MiSNPDb, available at http://webtom.cabgrid.res.in/mangosnps/ is based on 3-tier architecture, developed using PHP, MySQL and Javascript. This web genomic resources can be of immense use in the development of high density linkage map, QTL discovery, varietal differentiation, traceability, genome finishing and SNP chip development for future GWAS in genomic selection program. We report here world’s first web-based genomic resources for genetic improvement and germplasm management of mango.

Identification of massive molecular markers in Echinochloa phyllopogon using a restriction-site associated DNA approach

Echinochloa phyllopogon proliferation seriously threatens rice production worldwide. We combined a restriction-site associated DNA (RAD) approach with Illumina DNA sequencing for rapid and mass discovery of simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers for E. phyllopogon. RAD tags were generated from the genomic DNA of two E. phyllopogon plants, and sequenced to produce 5197.7 Mb and 5242.9 Mb high quality sequences, respectively. The GC content of E. phyllopogon was 45.8%, which is high for monocots. In total, 4710 putative SSRs were identified in 4132 contigs, which permitted the design of PCR primers for E. phyllopogon. Most repeat motifs among the SSRs identified were dinucleotide (>82%), and most of these SSRs were four motif-repeats (>75%). The most frequent motif was AT, accounting for 36.3%-37.2%, followed by AG and AC. In total, 78 putative polymorphic SSR loci were found. A total of 49,179 SNPs were discovered between the two samples of E. phyllopogon, 67.1% of which were transversions and 32.9% were transitions. We used eight SSRs to study the genetic diversity of four E. phyllopogon populations collected from rice fields in China and all eight loci tested were polymorphic.

Effect of Co-segregating Markers on High-Density Genetic Maps and Prediction of Map Expansion Using Machine Learning Algorithms

Advances in sequencing and genotyping methods have enable cost-effective production of high throughput single nucleotide polymorphism (SNP) markers, making them the choice for linkage mapping. As a result, many laboratories have developed high-throughput SNP assays and built high-density genetic maps. However, the number of markers may, by orders of magnitude, exceed the resolution of recombination for a given population size so that only a minority of markers can accurately be ordered. Another issue attached to the so-called 'large p, small n' problem is that high-density genetic maps inevitably result in many markers clustering at the same position (co-segregating markers). While there are a number of related papers, none have addressed the impact of co-segregating markers on genetic maps. In the present study, we investigated the effects of co-segregating markers on high-density genetic map length and marker order using empirical data from two populations of wheat, Mohawk × Cocorit (durum wheat) and Norstar × Cappelle Desprez (bread wheat). The maps of both populations consisted of 85% co-segregating markers. Our study clearly showed that excess of co-segregating markers can lead to map expansion, but has little effect on markers order. To estimate the inflation factor (IF), we generated a total of 24,473 linkage maps (8,203 maps for Mohawk × Cocorit and 16,270 maps for Norstar × Cappelle Desprez). Using seven machine learning algorithms, we were able to predict with an accuracy of 0.7 the map expansion due to the proportion of co-segregating markers. For example in Mohawk × Cocorit, with 10 and 80% co-segregating markers the length of the map inflated by 4.5 and 16.6%, respectively. Similarly, the map of Norstar × Cappelle Desprez expanded by 3.8 and 11.7% with 10 and 80% co-segregating markers. With the increasing number of markers on SNP-chips, the proportion of co-segregating markers in high-density maps will continue to increase making map expansion unavoidable. Therefore, we suggest developers improve linkage mapping algorithms for efficient analysis of high-throughput data. This study outlines a practical strategy to estimate the IF due to the proportion of co-segregating markers and outlines a method to scale the length of the map accordingly.

Genotyping-by-sequencing targeting of a novel downy mildew resistance gene Pl 20 from wild Helianthus argophyllus for sunflower (Helianthus annuus L.)

Genotyping-by-sequencing revealed a new downy mildew resistance gene, Pl ₂₀ , from wild Helianthus argophyllus located on linkage group 8 of the sunflower genome and closely linked to SNP markers that facilitate the marker-assisted selection of resistance genes. Downy mildew (DM), caused by Plasmopara halstedii, is one of the most devastating and yield-limiting diseases of sunflower. Downy mildew resistance identified in wild Helianthus argophyllus accession PI 494578 was determined to be effective against the predominant and virulent races of P. halstedii occurring in the United States. The evaluation of 114 BC₁F_2:3 families derived from the cross between HA 89 and PI 494578 against P. halstedii race 734 revealed that single dominant gene controls downy mildew resistance in the population. Genotyping-by-sequencing analysis conducted in the BC₁F₂ population indicated that the DM resistance gene derived from wild H. argophyllus PI 494578 is located on the upper end of the linkage group (LG) 8 of the sunflower genome, as was determined single nucleotide polymorphism (SNP) markers associated with DM resistance. Analysis of 11 additional SNP markers previously mapped to this region revealed that the resistance gene, named Pl ₂₀ , co-segregated with four markers, SFW02745, SFW09076, S8_11272025, and S8_11272046, and is flanked by SFW04358 and S8_100385559 at an interval of 1.8 cM. The newly discovered P. halstedii resistance gene has been introgressed from wild species into cultivated sunflower to provide a novel gene with DM resistance. The homozygous resistant individuals were selected from BC₂F₂ progenies with the use of markers linked to the Pl ₂₀ gene, and these lines should benefit the sunflower community for Helianthus improvement.

Linkage Mapping and Comparative Genomics of Red Drum (Sciaenops ocellatus) Using Next-Generation Sequencing

Developments in next-generation sequencing allow genotyping of thousands of genetic markers across hundreds of individuals in a cost-effective manner. Because of this, it is now possible to rapidly produce dense genetic linkage maps for nonmodel species. Here, we report a dense genetic linkage map for red drum, a marine fish species of considerable economic importance in the southeastern United States and elsewhere. We used a prior microsatellite-based linkage map as a framework and incorporated 1794 haplotyped contigs derived from high-throughput, reduced representation DNA sequencing to produce a linkage map containing 1794 haplotyped restriction-site associated DNA (RAD) contigs, 437 anonymous microsatellites, and 44 expressed sequence-tag-linked microsatellites (EST-SSRs). A total of 274 candidate genes, identified from transcripts from a preliminary hydrocarbon exposure study, were localized to specific chromosomes, using a shared synteny approach. The linkage map will be a useful resource for red drum commercial and restoration aquaculture, and for better understanding and managing populations of red drum in the wild.

Development and dissection of diagnostic SNP markers for the downy mildew resistance genes Pl Arg and Pl 8 and maker-assisted gene pyramiding in sunflower (Helianthus annuus L.)

Diagnostic DNA markers are an invaluable resource in breeding programs for successful introgression and pyramiding of disease resistance genes. Resistance to downy mildew (DM) disease in sunflower is mediated by Pl genes which are known to be effective against the causal fungus, Plasmopara halstedii. Two DM resistance genes, Pl _Arg and Pl ₈ , are highly effective against P. halstedii races in the USA, and have been previously mapped to the sunflower linkage groups (LGs) 1 and 13, respectively, using simple sequence repeat (SSR) markers. In this study, we developed high-density single nucleotide polymorphism (SNP) maps encompassing the Pl _arg and Pl ₈ genes and identified diagnostic SNP markers closely linked to these genes. The specificity of the diagnostic markers was validated in a highly diverse panel of 548 sunflower lines. Dissection of a large marker cluster co-segregated with Pl _Arg revealed that the closest SNP markers NSA_007595 and NSA_001835 delimited Pl _Arg to an interval of 2.83 Mb on the LG1 physical map. The SNP markers SFW01497 and SFW06597 delimited Pl ₈ to an interval of 2.85 Mb on the LG13 physical map. We also developed sunflower lines with homozygous, three gene pyramids carrying Pl _Arg , Pl ₈ , and the sunflower rust resistance gene R ₁₂ using the linked SNP markers from a segregating F₂ population of RHA 340 (carrying Pl ₈ )/RHA 464 (carrying Pl _Arg and R ₁₂ ). The high-throughput diagnostic SNP markers developed in this study will facilitate marker-assisted selection breeding, and the pyramided sunflower lines will provide durable resistance to downy mildew and rust diseases.

Discovery and introgression of the wild sunflower-derived novel downy mildew resistance gene Pl 19 in confection sunflower (Helianthus annuus L.)

A new downy mildew resistance gene, Pl ₁₉ , was identified from wild Helianthus annuus accession PI 435414, introduced to confection sunflower, and genetically mapped to linkage group 4 of the sunflower genome. Wild Helianthus annuus accession PI 435414 exhibited resistance to downy mildew, which is one of the most destructive diseases to sunflower production globally. Evaluation of the 140 BC₁F_2:3 families derived from the cross of CMS CONFSCLB1 and PI 435414 against Plasmopara halstedii race 734 revealed that a single dominant gene controls downy mildew resistance in the population. Bulked segregant analysis conducted in the BC₁F₂ population with 860 simple sequence repeat (SSR) markers indicated that the resistance derived from wild H. annuus was associated with SSR markers located on linkage group (LG) 4 of the sunflower genome. To map and tag this resistance locus, designated Pl ₁₉ , 140 BC₁F₂ individuals were used to construct a linkage map of the gene region. Two SSR markers, ORS963 and HT298, were linked to Pl ₁₉ within a distance of 4.7 cM. After screening 27 additional single nucleotide polymorphism (SNP) markers previously mapped to this region, two flanking SNP markers, NSA_003564 and NSA_006089, were identified as surrounding the Pl ₁₉ gene at a distance of 0.6 cM from each side. Genetic analysis indicated that Pl ₁₉ is different from Pl ₁₇ , which had previously been mapped to LG4, but is closely linked to Pl ₁₇ . This new gene is highly effective against the most predominant and virulent races of P. halstedii currently identified in North America and is the first downy mildew resistance gene that has been transferred to confection sunflower. The selected resistant germplasm derived from homozygous BC₂F₃ progeny provides a novel gene for use in confection sunflower breeding programs.

SNP Discovery and QTL Mapping of Sclerotinia Basal Stalk Rot Resistance in Sunflower using Genotyping-by-Sequencing

Basal stalk rot (BSR), caused by the ascomycete fungus (Lib.) de Bary, is a serious disease of sunflower ( L.) in the cool and humid production areas of the world. Quantitative trait loci (QTL) for BSR resistance were identified in a sunflower recombinant inbred line (RIL) population derived from the cross HA 441 × RHA 439. A genotyping-by-sequencing (GBS) approach was adapted to discover single nucleotide polymorphism (SNP) markers. A genetic linkage map was developed comprised of 1053 SNP markers on 17 linkage groups (LGs) spanning 1401.36 cM. The RILs were tested in five environments (locations and years) for resistance to BSR. Quantitative trait loci were identified in each environment separately and also with integrated data across environments. A total of six QTL were identified in all five environments: one of each on LGs 4, 9, 10, 11, 16, and 17. The most significant QTL, and , were identified at multiple environments on LGs 10 and 17, explaining 31.6 and 20.2% of the observed phenotypic variance, respectively. The remaining four QTL, , , , and , were detected in only one environment on LGs 4, 9, 11, and 16, respectively. Each of these QTL explains between 6.4 and 10.5% of the observed phenotypic variation in the RIL population. Alleles conferring increased resistance were contributed by both parents. The potential of the and in marker-assisted selection (MAS) breeding are discussed.

Genetics and mapping of a novel downy mildew resistance gene, Pl(18), introgressed from wild Helianthus argophyllus into cultivated sunflower (Helianthus annuus L.)

A novel downy mildew resistance gene, Pl(18), was introgressed from wild Helianthus argophyllus into cultivated sunflower and genetically mapped to linkage group 2 of the sunflower genome. The new germplasm, HA-DM1, carrying Pl(18) has been released to the public. Sunflower downy mildew (DM) is considered to be the most destructive foliar disease that has spread to every major sunflower-growing country of the world, except Australia. A new dominant downy mildew resistance gene (Pl 18) transferred from wild Helianthus argophyllus (PI 494573) into cultivated sunflower was mapped to linkage group (LG) 2 of the sunflower genome using bulked segregant analysis with 869 simple sequence repeat (SSR) markers. Phenotyping 142 BC1F2:3 families derived from the cross of HA 89 and H. argophyllus confirmed the single gene inheritance of resistance. Since no other Pl gene has been mapped to LG2, this gene was novel and designated as Pl (18). SSR markers CRT214 and ORS203 flanked Pl(18) at a genetic distance of 1.1 and 0.4 cM, respectively. Forty-six single nucleotide polymorphism (SNP) markers that cover the Pl(18) region were surveyed for saturation mapping of the region. Six co-segregating SNP markers were 1.2 cM distal to Pl(18), and another four co-segregating SNP markers were 0.9 cM proximal to Pl(18). The new BC2F4-derived germplasm, HA-DM1, carrying Pl(18) has been released to the public. This new line is highly resistant to all Plasmopara halstedii races identified in the USA providing breeders with an effective new source of resistance against downy mildew in sunflower. The molecular markers that were developed will be especially useful in marker-assisted selection and pyramiding of Pl resistance genes because of their close proximity to the gene and the availability of high-throughput SNP detection assays.

SNP discovery in common bean by restriction-associated DNA (RAD) sequencing for genetic diversity and population structure analysis

Researchers have made great advances into the development and application of genomic approaches for common beans, creating opportunities to driving more real and applicable strategies for sustainable management of the genetic resource towards plant breeding. This work provides useful polymorphic single-nucleotide polymorphisms (SNPs) for high-throughput common bean genotyping developed by RAD (restriction site-associated DNA) sequencing. The RAD tags were generated from DNA pooled from 12 common bean genotypes, including breeding lines of different gene pools and market classes. The aligned sequences identified 23,748 putative RAD-SNPs, of which 3357 were adequate for genotyping; 1032 RAD-SNPs with the highest ADT (assay design tool) score are presented in this article. The RAD-SNPs were structurally annotated in different coding (47.00 %) and non-coding (53.00 %) sequence components of genes. A subset of 384 RAD-SNPs with broad genome distribution was used to genotype a diverse panel of 95 common bean germplasms and revealed a successful amplification rate of 96.6 %, showing 73 % of polymorphic SNPs within the Andean group and 83 % in the Mesoamerican group. A slightly increased He (0.161, n = 21) value was estimated for the Andean gene pool, compared to the Mesoamerican group (0.156, n = 74). For the linkage disequilibrium (LD) analysis, from a group of 580 SNPs (289 RAD-SNPs and 291 BARC-SNPs) genotyped for the same set of genotypes, 70.2 % were in LD, decreasing to 0.10 %in the Andean group and 0.77 % in the Mesoamerican group. Haplotype patterns spanning 310 Mb of the genome (60 %) were characterized in samples from different origins. However, the haplotype frameworks were under-represented for the Andean (7.85 %) and Mesoamerican (5.55 %) gene pools separately. In conclusion, RAD sequencing allowed the discovery of hundreds of useful SNPs for broad genetic analysis of common bean germplasm. From now, this approach provides an excellent panel of molecular tools for whole genome analysis, allowing integrating and better exploring the common bean breeding practices.

Diversity analysis and genomic prediction of Sclerotinia resistance in sunflower using a new 25 K SNP genotyping array

KEY MESSAGE

We have developed a SNP array for sunflower containing more than 25 K markers, representing single loci mostly in or near transcribed regions of the genome. The array was successfully applied to genotype a diversity panel of lines, hybrids, and mapping populations and represented well the genetic diversity of cultivated sunflower. Results of PCoA and population substructure analysis underlined the complexity of the genetic composition of current elite breeding material. The performance of this genotyping platform for genome-based prediction of phenotypes and detection of QTL with improved resolution could be demonstrated based on the re-evaluation of a population segregating for resistance to Sclerotinia midstalk rot. Given our results, the newly developed 25 K SNP array is expected to be of great utility for the most important applications in genome-based sunflower breeding and research.

ABSTRACT

Genotyping with a large number of molecular markers is a prerequisite to conduct genome-based genetic analyses with high precision. Here, we report the design and performance of a 25 K SNP genotyping array for sunflower (Helianthus annuus L.). SNPs were discovered based on variant calling in de novo assembled, UniGene-based contigs of sunflower derived from whole genome sequencing and amplicon sequences originating from four and 48 inbred lines, respectively. After inclusion of publically available transcriptome-derived SNPs, in silico design of the Illumina(®) Infinium iSelect HD BeadChip yielded successful assays for 22,299 predominantly haplotype-specific SNPs. The array was validated in a sunflower diversity panel including inbred lines, open-pollinated varieties, introgression lines, landraces, recombinant inbred lines, and F2 populations. Validation provided 20,502 high-quality bi-allelic SNPs with stable cluster performance whereby each SNP marker represents a single locus mostly in or near transcribed regions of the sunflower genome. Analyses of population structure and quantitative resistance to Sclerotinia midstalk rot demonstrate that this array represents a significant improvement over currently available genomic tools for genetic diversity analyses, genome-wide marker-trait association studies, and genetic mapping in sunflower.

Large-Scale SNP Discovery and Genotyping for Constructing a High-Density Genetic Map of Tea Plant Using Specific-Locus Amplified Fragment Sequencing (SLAF-seq)

Genetic maps are important tools in plant genomics and breeding. The present study reports the large-scale discovery of single nucleotide polymorphisms (SNPs) for genetic map construction in tea plant. We developed a total of 6,042 valid SNP markers using specific-locus amplified fragment sequencing (SLAF-seq), and subsequently mapped them into the previous framework map. The final map contained 6,448 molecular markers, distributing on fifteen linkage groups corresponding to the number of tea plant chromosomes. The total map length was 3,965 cM, with an average inter-locus distance of 1.0 cM. This map is the first SNP-based reference map of tea plant, as well as the most saturated one developed to date. The SNP markers and map resources generated in this study provide a wealth of genetic information that can serve as a foundation for downstream genetic analyses, such as the fine mapping of quantitative trait loci (QTL), map-based cloning, marker-assisted selection, and anchoring of scaffolds to facilitate the process of whole genome sequencing projects for tea plant.

The First High-Density Genetic Map Construction in Tree Peony (Paeonia Sect. Moutan) using Genotyping by Specific-Locus Amplified Fragment Sequencing

Genetic linkage maps, permitting the elucidation of genome structure, are one of most powerful genomic tools to accelerate marker-assisted breeding. However, due to a lack of sufficient user-friendly molecular markers, no genetic linkage map has been developed for tree peonies (Paeonia Sect. Moutan), a group of important horticultural plants worldwide. Specific-locus amplified fragment sequencing (SLAF-seq) is a recent molecular marker development technology that enable the large-scale discovery and genotyping of sequence-based marker in genome-wide. In this study, we performed SLAF sequencing of an F1 population, derived from the cross P. ostti ‘FenDanBai’ × P. × suffruticosa ‘HongQiao’, to identify sufficient high-quality markers for the construction of high-density genetic linkage map in tree peonies. After SLAF sequencing, a total of 78 Gb sequencing data and 285,403,225 pair-end reads were generated. We detected 309,198 high-quality SLAFs from these data, of which 85,124 (27.5%) were polymorphic. Subsequently, 3518 of the polymorphic markers, which were successfully encoded in to Mendelian segregation types, and were in conformity with the criteria of high-quality markers, were defined as effective markers and used for genetic linkage mapping. Finally, we constructed an integrated genetic map, which comprised 1189 markers on the five linkage groups, and spanned 920.699 centiMorgans (cM) with an average inter-marker distance of 0.774 cM. There were 1115 ‘SNP-only’ markers, 18 ‘InDel-only’ markers, and 56 ‘SNP&InDel’ markers on the map. Among these markers, 450 (37.85%) showed significant segregation distortion (P < 0.05). In conclusion, this investigation reported the first large-scale marker development and high-density linkage map construction for tree peony. The results of this study will serve as a solid foundation not only for marker-assisted breeding, but also for genome sequence assembly for tree peony.

Molecular Mapping of Restriction-Site Associated DNA Markers In Allotetraploid Upland Cotton

Upland cotton (Gossypium hirsutum L., 2n = 52, AADD) is an allotetraploid, therefore the discovery of single nucleotide polymorphism (SNP) markers is difficult. The recent emergence of genome complexity reduction technologies based on the next-generation sequencing (NGS) platform has greatly expedited SNP discovery in crops with highly repetitive and complex genomes. Here we applied restriction-site associated DNA (RAD) sequencing technology for de novo SNP discovery in allotetraploid cotton. We identified 21,109 SNPs between the two parents and used these for genotyping of 161 recombinant inbred lines (RILs). Finally, a high dense linkage map comprising 4,153 loci over 3500-cM was developed based on the previous result. Using this map quantitative trait locus (QTLs) conferring fiber strength and Verticillium Wilt (VW) resistance were mapped to a more accurate region in comparison to the 1576-cM interval determined using the simple sequence repeat (SSR) genetic map. This suggests that the newly constructed map has more power and resolution than the previous SSR map. It will pave the way for the rapid identification of the marker-assisted selection in cotton breeding and cloning of QTL of interest traits.

Pl(17) is a novel gene independent of known downy mildew resistance genes in the cultivated sunflower (Helianthus annuus L.)

Pl 17, a novel downy mildew resistance gene independent of known downy mildew resistance genes in sunflowers, was genetically mapped to linkage group 4 of the sunflower genome. Downy mildew (DM), caused by Plasmopara halstedii (Farl.). Berl. et de Toni, is one of the serious sunflower diseases in the world due to its high virulence and the variability of the pathogen. DM resistance in the USDA inbred line, HA 458, has been shown to be effective against all virulent races of P. halstedii currently identified in the USA. To determine the chromosomal location of this resistance, 186 F 2:3 families derived from a cross of HA 458 with HA 234 were phenotyped for their resistance to race 734 of P. halstedii. The segregation ratio of the population supported that the resistance was controlled by a single dominant gene, Pl 17. Simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) primers were used to identify molecular markers linked to Pl 17. Bulked segregant analysis using 849 SSR markers located Pl 17 to linkage group (LG) 4, which is the first DM gene discovered in this linkage group. An F2 population of 186 individuals was screened with polymorphic SSR and SNP primers from LG4. Two flanking markers, SNP SFW04052 and SSR ORS963, delineated Pl 17 in an interval of 3.0 cM. The markers linked to Pl 17 were validated in a BC3 population. A search for the physical location of flanking markers in sunflower genome sequences revealed that the Pl 17 region had a recombination frequency of 0.59 Mb/cM, which was a fourfold higher recombination rate relative to the genomic average. This region can be considered amenable to molecular manipulation for further map-based cloning of Pl 17.

Relocation of a rust resistance gene R 2 and its marker-assisted gene pyramiding in confection sunflower (Helianthus annuus L.)

The rust resistance gene R 2 was reassigned to linkage group 14 of the sunflower genome. DNA markers linked to R 2 were identified and used for marker-assisted gene pyramiding in a confection type genetic background. Due to the frequent evolution of new pathogen races, sunflower rust is a recurring threat to sunflower production worldwide. The inbred line Morden Cross 29 (MC29) carries the rust resistance gene, R 2 , conferring resistance to numerous races of rust fungus in the US, Canada, and Australia, and can be used as a broad-spectrum resistance resource. Based on phenotypic assessments and SSR marker analyses on the 117 F2 individuals derived from a cross of HA 89 with MC29 (USDA), R 2 was mapped to linkage group (LG) 14 of the sunflower, and not to the previously reported location on LG9. The closest SSR marker HT567 was located at 4.3 cM distal to R 2 . Furthermore, 36 selected SNP markers from LG14 were used to saturate the R 2 region. Two SNP markers, NSA_002316 and SFW01272, flanked R 2 at a genetic distance of 2.8 and 1.8 cM, respectively. Of the three closely linked markers, SFW00211 amplified an allele specific for the presence of R 2 in a marker validation set of 46 breeding lines, and SFW01272 was also shown to be diagnostic for R 2 . These newly developed markers, together with the previously identified markers linked to the gene R 13a , were used to screen 524 F2 individuals from a cross of a confection R 2 line and HA-R6 carrying R 13a . Eleven homozygous double-resistant F2 plants with the gene combination of R 2 and R 13a were obtained. This double-resistant line will be extremely useful in confection sunflower, where few rust R genes are available, risking evolution of new virulence phenotypes and further disease epidemics.

De novo assembly of transcriptome sequencing in Caragana korshinskii Kom. and characterization of EST-SSR markers

Caragana korshinskii Kom. is widely distributed in various habitats, including gravel desert, clay desert, fixed and semi-fixed sand, and saline land in the Asian and African deserts. To date, no previous genomic information or EST-SSR marker has been reported in Caragana Fabr. genus. In this study, more than two billion bases of high-quality sequence of C. korshinskii were generated by using illumina sequencing technology and demonstrated the de novo assembly and annotation of genes without prior genome information. These reads were assembled into 86,265 unigenes (mean length = 709 bp). The similarity search indicated that 33,955 and 21,978 unigenes showed significant similarities to known proteins from NCBI non-redundant and Swissprot protein databases, respectively. Among these annotated unigenes, 26,232 a unigenes were separately assigned to Gene Ontology (GO) database. When 22,756 unigenes searched against the Kyoto Encyclopedia of Genes and Genomes Pathway (KEGG) database, 5,598 unigenes were assigned to 5 main categories including 32 KEGG pathways. Among the main KEGG categories, metabolism was the biggest category (2,862, 43.7%), suggesting the active metabolic processes in the desert tree. In addition, a total of 19,150 EST-SSRs were identified from 15,484 unigenes, and the characterizations of EST-SSRs were further compared with other four species in Fabraceae. 126 potential marker sites were randomly selected to validate the assembly quality and develop EST-SSR markers. Among the 9 germplasms in Caranaga Fabr. genus, PCR success rate were 93.7% and the phylogenic tree was constructed based on the genotypic data. This research generated a substantial fraction of transcriptome sequences, which were very useful resources for gene annotation and discovery, molecular markers development, genome assembly and annotation. The EST-SSR markers identified and developed in this study will facilitate marker-assisted selection breeding.

De novo assembly of the desert tree Haloxylon ammodendron (C. A. Mey.) based on RNA-Seq data provides insight into drought response, gene discovery and marker identification

Background

Haloxylon ammodendron (C. A. Mey.) is widely distributed across a range of habitats, including gravel desert, clay desert, fixed and semi-fixed sand, and saline land in Asian and African deserts. To date, no genomic information or expressed sequence tag-simple sequence repeat (EST-SSR) marker has been reported for H. ammodendron plants.

Results

Using Illumina sequencing technology, we generated over two billion bases of high-quality sequence data on H. ammodendron and conducted de novo assembly and annotation of genes without prior genome information. These reads were assembled into 79,918 unigenes (mean length = 728 bp). Based on similarity searches comparing these unigenes with known proteins in the non-redundant (nr) protein database, 25,619 unigenes were functionally annotated with a cut-off E-value of 10^-5. In addition, DGE reads were mapped to the assembled transcriptome for gene expression analysis under drought stress. In total, 1,060 differentially expressed genes were identified. Among these genes, 356 genes were upregulated after drought treatment, and 704 genes were downregulated. We used the KEGG database to annotate these drought-induced genes; 207 unigenes were identified in the KEGG pathway annotation, and approximately 12.1% of the unigenes with known function fell into categories related to fatty acid metabolism, starch and sucrose metabolism, and nitrogen metabolism, suggesting that these pathways or processes may be involved in the drought response. Together, a total of 35 drought-inducible transcription factors were identified, including WRKY, MYB and bZIP family members.

Conclusions

Our study is the first to provide a transcriptome sequence resource for H. ammodendron plants and to determine its digital gene expression profile under drought conditions using the assembled transcriptome data for reference. These data provide a valuable resource for genetic and genomic studies of desert plants under abiotic conditions.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1111) contains supplementary material, which is available to authorized users.