BSTA: a targeted approach combines bulked segregant analysis with next- generation sequencing and de novo transcriptome assembly for SNP discovery in sunflower

Identification of salt stress-tolerant candidate genes in the BC2F2 population at the seedling stages of G. hirsutum and G. darwinii using NGS-based bulked segregant analysis

Salinity is a major threat to the yield and productivity of cotton seedlings. In the present study, we developed a BC₂F₂ population of cotton plants from Gossypium darwinii (5-7) and Gossypium hirsutum (CCRI 12-4) salt-susceptible parents to identify salt-resistant candidate genes. The Illumina HiSeq™ strategy was used with bulked segregant analysis. Salt-resistant and salt-susceptible DNA bulks were pooled by using 30 plants from a BC₂F₂ population. Next-generation sequencing (NGS) technology was used for the sequencing of parents and both bulks. Four significant genomic regions were identified: the first genomic region was located on chromosome 18 (1.86 Mb), the second and third genomic regions were on chromosome 25 (1.06 Mb and 1.94 Mb, respectively), and the fourth was on chromosome 8 (1.41 Mb). The reads of bulk1 and bulk2 were aligned to the G. darwinii and G. hirsutum genomes, respectively, leading to the identification of 20,664,007 single-nucleotide polymorphisms (SNPs) and insertions/deletions (indels). After the screening, 6,573 polymorphic markers were obtained after filtration of the candidate regions. The SNP indices in resistant and susceptible bulks and Δ(SNP-index) values of resistant and susceptible bulks were measured. Based on the higher Δ(SNP-index) value, six effective polymorphic SNPs were selected in a different chromosome. Six effective SNPs were linked to five candidate genes in four genomic regions. Further validation of these five candidate genes was carried out using reverse transcription-quantitative polymerase chain reaction (RT-qPCR), resulting in an expression profile that showed two highly upregulated genes in the salt-tolerant species G. darwinii, i.e., Gohir.D05G367800 and Gohir.D12G239100; however, the opposite was shown in G. hirsutum, for which all genes, except one, showed partial expression. The results indicated that Gohir.D05G367800 and Gohir.D12G239100 may be salt-tolerant genes. We are confident that this study could be helpful for the cloning, transformation, and development of salt-resistant cotton varieties.

BSR and Full-Length Transcriptome Approaches Identified Candidate Genes for High Seed Ratio in Camellia vietnamensis

(1) Background: C. vietnamensis is very suitable for growth in the low hilly areas of southern subtropical regions. Under appropriate conditions, the oil yield of C. vietnamensis can reach 1125 kg/ha (the existing varieties can reach 750 kg/ha). Moreover, the fruit of C. vietnamensis is large and the pericarp is thick (>5 cm). Therefore, a high seed ratio has become the main target economic trait for the breeding of C. vietnamensis. (2) Methods: A half-sibling population of C. vietnamensis plants with a combination of high and low seed ratios was constructed by crossing a C. vietnamensis female parent. Bulked segregant RNA analysis and full-length transcriptome sequencing were performed to determine the molecular mechanisms underlying a high seed ratio. (3) Results: Seed ratio is a complex quantitative trait with a normal distribution, which is significantly associated with four other traits of fruit (seed weight, seed number, fruit diameter, and pericarp thickness). Two candidate regions related to high seed ratio (HSR) were predicted. One spanned 140.8−148.4 Mb of chromosome 2 and was associated with 97 seed-yield-related candidate genes ranging in length from 278 to 16,628 bp. The other spanned 35.3−37.3 Mb on chromosome 15 and was associated with 38 genes ranging in length from 221 to 16,928 bp. Using the full-length transcript as a template, a total of 115 candidate transcripts were obtained, and 78 transcripts were predicted to be functionally annotated. The DEGs from two set pairs of cDNA sequencing bulks were enriched to cytochrome p450 CYP76F14 (KOG0156; GO:0055114, HSR4, HSR7), the gibberellin phytohormone pathway (GO:0016787, HSR5), the calcium signaling pathway (GO:0005509, HSR6), the polyubiquitin-PPAR signaling pathway (GO:0005515, HSR2, HSR3), and several main transcription factors (bZIP transcription factor, HSR1) in C. vietnamensis.

Conjunctive Analyses of BSA-Seq and BSR-Seq Unveil the Msβ-GAL and MsJMT as Key Candidate Genes for Cytoplasmic Male Sterility in Alfalfa (Medicago sativa L.)

Knowing the molecular mechanism of male sterility in alfalfa is important to utilize the heterosis more effectively. However, the molecular mechanisms of male sterility in alfalfa are still unclear. In this study, the bulked segregant analysis (BSA) and bulked segregant RNA-seq (BSR) were performed with F2 separation progeny to study the molecular mechanism of male sterility in alfalfa. The BSA-seq analysis was located in a candidate region on chromosome 5 containing 626 candidate genes which were associated with male sterility in alfalfa, while the BSR-seq analysis filtered seven candidate DEGs related to male sterility, and these candidate genes including EF-Tu, β-GAL, CESA, PHGDH, and JMT. The conjunctive analyses of BSR and BSA methods revealed that the genes of Msβ-GAL and MsJMT are the common detected candidate genes involved in male sterility in alfalfa. Our research provides a theory basis for further study of the molecular mechanism of male sterility in alfalfa and significant information for the genetic breeding of Medicago sativa.

Identification of genes associated with the biosynthesis of unsaturated fatty acid and oil accumulation in herbaceous peony 'Hangshao' (Paeonia lactiflora 'Hangshao') seeds based on transcriptome analysis

Background

Paeonia lactiflora ‘Hangshao’ is widely cultivated in China as a traditional Chinese medicine ‘Radix Paeoniae Alba’. Due to the abundant unsaturated fatty acids in its seed, it can also be regarded as a new oilseed plant. However, the process of the biosynthesis of unsaturated fatty acids in it has remained unknown. Therefore, transcriptome analysis is helpful to better understand the underlying molecular mechanisms.

Results

Five main fatty acids were detected, including stearic acid, palmitic acid, oleic acid, linoleic acid and α-linolenic acid, and their absolute contents first increased and then decreased during seed development. A total of 150,156 unigenes were obtained by transcriptome sequencing. There were 15,005 unigenes annotated in the seven functional databases, including NR, NT, GO, KOG, KEGG, Swiss-Prot and InterPro. Based on the KEGG database, 1766 unigenes were annotated in the lipid metabolism. There were 4635, 12,304, and 18,291 DEGs in Group I (60 vs 30 DAF), Group II (90 vs 60 DAF) and Group III (90 vs 30 DAF), respectively. A total of 1480 DEGs were detected in the intersection of the three groups. In 14 KEGG pathways of lipid metabolism, 503 DEGs were found, belonging to 111 enzymes. We screened out 123 DEGs involved in fatty acid biosynthesis (39 DEGs), fatty acid elongation (33 DEGs), biosynthesis of unsaturated fatty acid (24 DEGs), TAG assembly (17 DEGs) and lipid storage (10 DEGs). Furthermore, qRT-PCR was used to analyze the expression patterns of 16 genes, including BBCP, BC, MCAT, KASIII, KASII, FATA, FATB, KCR, SAD, FAD2, FAD3, FAD7, GPAT, DGAT, OLE and CLO, most of which showed the highest expression at 45 DAF, except for DGAT, OLE and CLO, which showed the highest expression at 75 DAF.

Conclusions

We predicted that MCAT, KASIII, FATA, SAD, FAD2, FAD3, DGAT and OLE were the key genes in the unsaturated fatty acid biosynthesis and oil accumulation in herbaceous peony seed. This study provides the first comprehensive genomic resources characterizing herbaceous peony seed gene expression at the transcriptional level. These data lay the foundation for elucidating the molecular mechanisms of fatty acid biosynthesis and oil accumulation for herbaceous peony.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07339-7.

QTL-Seq identifies quantitative trait loci of relative electrical conductivity associated with heat tolerance in bottle gourd (Lagenaria siceraria)

Heat is a major abiotic stress that seriously affects watermelon (Citrullus lanatus) production. However, its effects may be mitigated through grafting watermelon to heat tolerant bottle gourd (Lagenaria siceraria) rootstocks. Understanding the genetic basis of heat tolerance and development of reliable DNA markers to indirectly select for the trait are necessary in breeding for new varieties with heat tolerance. The objectives of this study were to investigate the inheritance of heat tolerance and identify molecular markers associated with heat tolerance in bottle gourd. A segregating F₂ population was developed from a cross between two heat tolerant and sensitive inbred lines. The population was phenotyped for relative electrical conductivity (REC) upon high temperature treatment which was used as an indicator for heat tolerance. QTL-seq was performed to identify regions associated with heat tolerance. We found that REC-based heat tolerance in this population exhibited recessive inheritance. Seven heat-tolerant quantitative trait loci (qHT1.1, qHT2.1, qHT2.2, qHT5.1, qHT6.1, qHT7.1, and qHT8.1) were identified with qHT2.1 being a promising major-effect QTL. In the qHT2.1 region, we identified three non-synonymous SNPs that were potentially associated with heat tolerance. These SNPs were located in the genes that may play roles in pollen sterility, intracellular transport, and signal recognition. Association of the three SNPs with heat tolerance was verified in segregating F₂ populations, which could be candidate markers for marker assisted selection for heat tolerance in bottle gourd. The qHT2.1 region is an important finding that may be used for fine mapping and discovery of novel genes associated with heat tolerance in bottle gourd.

Conjunctive Analyses of BSA-Seq and BSR-Seq to Reveal the Molecular Pathway of Leafy Head Formation in Chinese Cabbage

As the storage organ of Chinese cabbage, the leafy head was harvested as a commercial product due to its edible value. In this study, the bulked segregant analysis (BSA) and bulked segregant RNA-Seq (BSR) were performed with F2 separation progeny to study the molecular mechanism of leafy head formation in Chinese cabbage. BSA-Seq analysis located four candidate regions containing 40 candidate genes, while BSR-Seq analysis revealed eight candidate regions containing 607 candidate genes. The conjunctive analyses of these two methods identified that Casein kinase gene BrCKL8 (Bra035974) is the common candidate gene related with leafy head formation in Chinese cabbage, and it showed high expression levels at the three segments of heading type plant leaves. The differentially expressed genes (DEGs) between two set pairs of cDNA sequencing bulks were divided into two categories: one category was related with five hormone pathways (Auxin, Ethylene, Abscisic acid, Jasmonic acid and Gibberellin), the other category was composed of genes that associate with the calcium signaling pathway. Moreover, a series of upregulated transcriptional factors (TFs) were also identified by the association analysis of BSR-Seq analysis. The leafy head development was regulated by various biological processes and effected by diverse external environment factors, so our research will contribute to the breeding of perfect leaf-heading types of Chinese cabbage.

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.

Mapping of the male sterile mutant gene ftms in Brassica rapa L. ssp. pekinensis via BSR-Seq combined with whole-genome resequencing

A male sterile mutant was created by ⁶⁰Co γ-rays of microspores isolated from Chinese cabbage DH line 'FT'. A candidate gene for the male sterile trait was identified as Bra010198. Male sterility is used for hybrid seed production in Chinese cabbage. In this study, we derived a male sterile mutant (ftms) from Chinese cabbage DH line 'FT' by irradiating microspores with ⁶⁰Co γ-rays and realized the rapid trait transformation from male fertility to sterility for creating valuable breeding materials. Genetic analysis indicated that the male sterile trait is controlled by a single recessive nuclear gene, ftms. Microspore development in mutant ftms was aborted at the tetrad stage and associated with severely retarded degeneration and vacuolation of tapetum. Using BSR-seq analysis, the candidate region for ftms was mapped on chromosome A05. A large F₂ population was created, and the region was narrowed to approximately 1.7-Mb between markers Indel20 and Indel14 via linkage analysis. The recombination frequency was extremely suppressed because the region was located on the chromosome A05 centromere. Whole-genome resequencing of mutant ftms and wild-type 'FT' aligned only one nonsynonymous SNP to Bra010198; this gene is a homolog of Arabidopsis KNS4/UPEX1, which encodes a putative β-(1,3)-galactosyltransferase that controls pollen exine development. Comparative sequencing verified the SNP position on the fifth exon of Bra010198 in mutant ftms. Further genotyping revealed that the male sterile phenotype was fully co-segregated with this SNP. Quantitative real-time PCR indicated that Bra0101918 specifically expressed in stamen. The data presented herein suggested that Bra010198 is a strong candidate gene for ftms. Hence, we developed a male sterile line for potential application in breeding and expanded the knowledge about the molecular mechanism underlying male sterility in Chinese cabbage.

Transcriptome Profiling and Characterization of Drought-Tolerant Potato Plant (Solanum tuberosum L.)

Potato (Solanum tuberosum L.) is the third most important food crop, and breeding drought-tolerant varieties is vital research goal. However, detailed molecular mechanisms in response to drought stress in potatoes are not well known. In this study, we developed EMS-mutagenized potatoes that showed significant tolerance to drought stress compared to the wild-type (WT) 'Desiree' cultivar. In addition, changes to transcripts as a result of drought stress in WT and drought-tolerant (DR) plants were investigated by de novo assembly using the Illumina platform. One-week-old WT and DR plants were treated with -1.8 Mpa polyethylene glycol-8000, and total RNA was prepared from plants harvested at 0, 6, 12, 24, and 48 h for subsequent RNA sequencing. In total, 61,100 transcripts and 5,118 differentially expressed genes (DEGs) displaying up- or down-regulation were identified in pairwise comparisons of WT and DR plants following drought conditions. Transcriptome profiling showed the number of DEGs with up-regulation and down-regulation at 909, 977, 1181, 1225 and 826 between WT and DR plants at 0, 6, 12, 24, and 48 h, respectively. Results of KEGG enrichment showed that the drought tolerance mechanism of the DR plant can mainly be explained by two aspects, the 'photosynthetic-antenna protein' and 'protein processing of the endoplasmic reticulum'. We also divided eight expression patterns in four pairwise comparisons of DR plants (DR0 vs DR6, DR12, DR24, DR48) under PEG treatment. Our comprehensive transcriptome data will further enhance our understanding of the mechanisms regulating drought tolerance in tetraploid potato cultivars.

Identification of a cold-tolerant locus in rice (Oryza sativa L.) using bulked segregant analysis with a next-generation sequencing strategy

BACKGROUND

Cold stress can cause serious abiotic damage that limits the growth, development and yield of rice. Cold tolerance during the booting stage of rice is a key factor that can guarantee a high and stable yield under cold stress. The cold tolerance of rice is controlled by quantitative trait loci (QTLs). Based on the complex genetic basis of cold tolerance in rice, additional efforts are needed to detect reliable QTLs and identify candidate genes. In this study, recombinant inbred lines (RILs) derived from a cross between a cold sensitive variety, Dongnong422, and strongly cold-tolerant variety, Kongyu131, were used to screen for cold-tolerant loci at the booting stage of rice.

RESULTS

A novel major QTL, qPSST6, controlling the percent seed set under cold water treatment (PSST) under the field conditions of 17 °C cold water irrigation was located on the 28.4 cM interval on chromosome 6. Using the combination of bulked-segregant analysis (BSA) and next-generation sequencing (NGS) technology (Seq-BSA), a 1.81 Mb region that contains 269 predicted genes on chromosome 6 was identified as the candidate region of qPSST6. Two genes, LOC_Os06g39740 and LOC_Os06g39750, were annotated as "response to cold" by gene ontology (GO) analysis. qRT-PCR analysis revealed that LOC_Os06g39750 was strongly induced by cold stress. Haplotype analysis also demonstrate a key role of LOC_Os06g39750 in regulating the PSST of rice, suggesting that it was the candidate gene of qPSST6.

CONCLUSIONS

The information obtained in this study is useful for gene cloning of qPSST6 and for breeding cold-tolerant varieties of rice using marker assisted selection (MAS).

Combining QTL-seq and linkage mapping to fine map a wild soybean allele characteristic of greater plant height

BACKGROUND

Plant height (PH) is an important agronomic trait and is closely related to yield in soybean [Glycine max (L.) Merr.]. Previous studies have identified many QTLs for PH. Due to the complex genetic background of PH in soybean, there are few reports on its fine mapping.

RESULTS

In this study, we used a mapping population derived from a cross between a chromosome segment substitution line CSSL3228 (donor N24852 (G. Soja), a receptor NN1138-2 (G. max)) and NN1138-2 to fine map a wild soybean allele of greater PH by QTL-seq and linkage mapping. We identified a QTL for PH in a 1.73 Mb region on soybean chromosome 13 through QTL-seq, which was confirmed by SSR marker-based classical QTL mapping in the mapping population. The linkage analysis showed that the QTLs of PH were located between the SSR markers BARCSOYSSR_13_1417 and BARCSOYSSR_13_1421 on chromosome 13, and the physical distance was 69.3 kb. RT-PCR and sequence analysis of possible candidate genes showed that Glyma.13 g249400 revealed significantly higher expression in higher PH genotypes, and the gene existed 6 differences in the amino acids encoding between the two parents.

CONCLUSIONS

Data presented here provide support for Glyma.13 g249400 as a possible candidate genes for higher PH in wild soybean line N24852.

Screening of Genes Related to Early and Late Flowering in Tree Peony Based on Bulked Segregant RNA Sequencing and Verification by Quantitative Real-Time PCR

Tree peony (Paeonia suffruticosa Andrews) is a perennial woody shrub bearing large and colorful flowers. However, the flowering period is short and relatively uniform, which to an important extent hinders the cultivation and exploitation of ornamental peonies. In this study, the segregation of an F₁ population derived from P. ostti ‘Feng Dan’ (an early-flowering cultivar) × P. suffruticosa ‘Xin Riyuejin’ (a late-flowering cultivar) was used to screen and analyze candidate genes associated with flowering period of the two parents. Extreme early- and late-flowering genotypes of the F₁ population at full-bloom stage were sampled to establish an early-flowering mixed pool (T03), a late-flowering mixed pool (T04), a late-flowering male pool (T01), and an early-flowering female pool (T02), using the Sequencing By Synthesis (SBS) technology on the Illumina HiSeq TM2500 platform. A total of 56.51 Gb of clean reads data, comprising at least 87.62% of Quality30 (Q30), was generated, which was then combined into 173,960 transcripts (N50 = 1781) and 78,645 (N50 = 1282) unigenes, with a mean length of 1106.76 and 732.27 base pairs (bp), respectively. Altogether, 58,084 genes were annotated by comparison with public databases, based on an E-value parameter of less than 10⁻⁵ and 10⁻¹⁰ for BLAST and HMMER, respectively. In total, 291 unigene sequences were finally screened out by BSR-seq (bulked segregant RNA-seq) association analysis. To validate these unigenes, we finally confirmed seven unigenes that were related to early and late flowering, which were then verified by quantitative real-time PCR (qRT-PCR). This is the first reported study to screen genes associated with early and late flowering of tree peony by the BSA (bulked sample analysis) method of transcriptome sequencing, and to construct a high-quality transcriptome database. A set of candidate functional genes related to flowering time was successfully obtained, providing an important genetic resource for further studies of flowering in peony and the mechanism of regulation of flowering time in tree peony.

Identification of a locus conferring dominant resistance to maize rough dwarf disease in maize

Maize rough dwarf disease (MRDD) is a severe viral disease of maize that occurs worldwide, particularly in the summer maize-growing areas in China, resulting in yield losses and quality deterioration in susceptible maize varieties. An effective solution to control MRDD is to use resistance genes to improve the behavior of susceptible genotypes. Here, we employed maize F₂ populations derived from a cross between susceptible line S221 and resistant line K36 for the deep sequencing of the two DNA pools containing extremely resistant and susceptible F₂ individuals, and used traditional linkage analysis to locate the resistance-related genomic region. The results showed that MRDD resistance in K36 was controlled by a single dominant locus, and an associated region was identified within the genomic interval of 68,396,487 bp and 69,523,478 bp on chromosome 6. Two simple sequence repeat (SSR) markers 6F29R29 and 6F34R34 were tightly linked to the MRDD resistance locus. The findings of the present study improve our understanding of the inheritance patterns of MRDD resistance, and should inform MRDD-resistant maize breeding programs.

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.

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.

Rapid identification of QTLs underlying resistance to Cucumber mosaic virus in pepper (Capsicum frutescens)

Next-generation sequencing enabled a fast discovery of QTLs controlling CMV resistant in pepper. The gene CA02g19570 as a possible candidate gene of qCmr2.1 was identified for resistance to CMV in pepper. Cucumber mosaic virus (CMV) is one of the most important viruses infecting pepper, but the genetic basis of CMV resistance in pepper is elusive. In this study, we identified a candidate gene for CMV resistance QTL, qCmr2.1 through SLAF-seq. Segregation analysis in F₂, BC₁ and F_2:3 populations derived from a cross between two inbred lines 'PBC688' (CMV-resistant) and 'G29' (CMV-susceptible) suggested quantitative inheritance of resistance to CMV in pepper. Genome-wide comparison of SNP profiles between the CMV-resistant and CMV-susceptible bulks constructed from an F₂ population identified two QTLs, designated as qCmr2.1 on chromosome 2 and qCmr11.1 on chromosome 11 for resistance to CMV in PBC688, which were confirmed by InDel marker-based classical QTL mapping in the F₂ population. As a major QTL, joint SLAF-seq and traditional QTL analysis delimited qCmr2.1 to a 330 kb genomic region. Two pepper genes, CA02g19570 and CA02g19600, were identified in this region, which are homologous with the genes LOC104113703, LOC104248995, LOC102603934 and LOC101248357, which were predicted to encode N-like protein associated with TMV-resistant in Solanum crops. Quantitative RT-PCR revealed higher expression levels of CA02g19570 in CMV resistance genotypes. The CA02g19600 did not exhibit obvious regularity in expression patterns. Higher relative expression levels of CA02g19570 in PBC688 and F₁ were compared with those in G29 during days after inoculation. These results provide support for CA02g19570 as a possible candidate gene of qCmr2.1 for resistance to CMV in pepper.

Bulked sample analysis in genetics, genomics and crop improvement

Biological assay has been based on analysis of all individuals collected from sample populations. Bulked sample analysis (BSA), which works with selected and pooled individuals, has been extensively used in gene mapping through bulked segregant analysis with biparental populations, mapping by sequencing with major gene mutants and pooled genomewide association study using extreme variants. Compared to conventional entire population analysis, BSA significantly reduces the scale and cost by simplifying the procedure. The bulks can be built by selection of extremes or representative samples from any populations and all types of segregants and variants that represent wide ranges of phenotypic variation for the target trait. Methods and procedures for sampling, bulking and multiplexing are described. The samples can be analysed using individual markers, microarrays and high-throughput sequencing at all levels of DNA, RNA and protein. The power of BSA is affected by population size, selection of extreme individuals, sequencing strategies, genetic architecture of the trait and marker density. BSA will facilitate plant breeding through development of diagnostic and constitutive markers, agronomic genomics, marker-assisted selection and selective phenotyping. Applications of BSA in genetics, genomics and crop improvement are discussed with their future perspectives.

Genetic mapping and molecular marker development for Pi65(t), a novel broad-spectrum resistance gene to rice blast using next-generation sequencing

A novel R gene was mapped to a locus on chromosome 11 from 30.42 to 30.85 Mb, which was proven to be efficient in the improvement of rice blast resistance. Rice blast is a devastating fungal disease worldwide. The use of blast resistance (R) genes is the most important approach to control the disease in rice breeding. In the present study, we finely mapped a novel resistance gene Pi65(t), conferring a broad-spectrum resistance to the fungus Magnaporthe oryzae, using bulked segregant analysis in combination with next-generation sequencing technology. Segregation in a doubled haploid (DH) population and a BC1F2 population suggested that resistance to blast in Gangyu129 was likely conferred by a single dominant gene, designated Pi65(t); it was located on chromosome 11 from 30.20 to 31.20 Mb using next-generation sequencing. After screening recombinants with newly developed molecular markers, the region was narrowed down to 0.43 Mb, flanked by SNP-2 and SNP-8 at the physical location from 30.42 to 30.85 Mb based on the Nipponbare reference database in build 5. Using the software QTL IciMapping, Pi65(t) was further mapped to a locus between InDel-1 and SNP-4 with genetic distances of 0.11 and 0.98 cM, respectively. Within this region, 4 predicted R genes were found with nucleotide binding site and leucine-rich repeat (NBS-LRR) domains. We developed molecular markers to genotype 305 DH lines and found that InDel-1 was closely linked with Pi65(t). Using InDel-1, a new rice variety Chuangxin1 containing Pi65(t) was developed, and it is highly resistant to rice blast and produces a high yield in Liaoning province of China. This indicated that Pi65(t) could play a key role in the improvement of rice blast resistance.

Mining candidate genes associated with powdery mildew resistance in cucumber via super-BSA by specific length amplified fragment (SLAF) sequencing

BACKGROUND

Powdery mildew (PM) is the most common fungal disease of cucumber and other cucurbit crops, while breeding the PM-resistant materials is the effective way to defense this disease, and the recent development of modern genetics and genomics make us aware of that studying the resistance genes is the essential way to breed the PM high-resistance plant. With the ever increasing throughput of next-generation sequencing (NGS), the development of specific length amplified fragment sequencing (SLAF-seq) as a high-resolution strategy for large-scale de novo SNP discovery is gradually applied for functional gene mining. Here we combined the bulked segregant analysis (BSA) with SLAF-seq to identify candidate genes associated with PM resistance in cucumber.

METHODS

A segregating population comprising 251 F2 individuals was developed using H136 (female parent) as susceptible parent and BK2 (male parent) as resistance donor. After PMR test, total genomic DNA was prepared from each plant. Systemic genomic analysis of the GC content, repeat sequence, etc. was carried out by prediction software SLAF_Predict to establish condition to ensure the uniformity and density of the molecular markers. After samples were gel purified, SLAFs were generated at Biomarker Technologies Corporation in Beijing. Based on SLAF tags and the PMR test result, the hot region were annotated.

RESULTS

A total of 73,100 high-quality SLAF tags with an average depth of 99.11× were sequenced. Among these, 5,355 polymorphic tags were identified with a polymorphism rate of 7.34 %, including 7.09 % SNPs and other polymorphism types. Finally, 140 associated SLAFs were identified, and two main Hot Regions were detected on chromosome 1 and 6, which contained five genes invovled in defense response, toxin metabolism, cell stress response, and injury response in cucumber.

CONCLUSIONS

Associated markers identified by super-BSA in this study, could not only speed up the study of the PMR genes, but also provide a feasible solution for breeding the marker-assisted PMR cucumber. Moreover, this study could also be extended to any other species with reference genome.

Comparative quantitative trait loci for silique length and seed weight in Brassica napus

Silique length (SL) and seed weight (SW) are important yield-associated traits in rapeseed (Brassica napus). Although many quantitative trait loci (QTL) for SL and SW have been identified in B. napus, comparative analysis for those QTL is seldom performed. In the present study, 20 and 21 QTL for SL and SW were identified in doubled haploid (DH) and DH-derived reconstructed F2 populations in rapeseed, explaining 55.1-74.3% and 24.4-62.9% of the phenotypic variation across three years, respectively. Of which, 17 QTL with partially or completely overlapped confidence interval on chromosome A09, were homologous with two overlapped QTL on chromosome C08 by aligning QTL confidence intervals with the reference genomes of Brassica crops. By high density selective genotyping of DH lines with extreme phenotypes, using a Brassica single-nucleotide polymorphism (SNP) array, the QTL on chromosome A09 was narrowed, and aligned into 1.14-Mb region from 30.84 to 31.98 Mb on chromosome R09 of B. rapa and 1.05-Mb region from 27.21 to 28.26 Mb on chromosome A09 of B. napus. The alignment of QTL with Brassica reference genomes revealed homologous QTL on A09 and C08 for SL. The narrowed QTL region provides clues for gene cloning and breeding cultivars by marker-assisted selection.

Transcriptome Profiling of the Potato (Solanum tuberosum L.) Plant under Drought Stress and Water-Stimulus Conditions

Drought stress can seriously affect tuberization, yield and quality of potato plant. However, the precise molecular mechanisms governing potato stolon’s response to drought stress and water supply are not very well understood. In this work, a potato (Solanum tuberosum L.) variant, Ningshu 4, was subjected to severe drought stress treatment (DT) and re-watering treatment (RWT) at tuber bulking stage. Strand-specific cDNA libraries of stolon materials were constructed for paired-end transcriptome sequencing analyses and differentially expressed gene (DEG) examination. In comparison to untreated-control (CT) plants, 3189 and 1797 DEGs were identified in DT and RWT plants and 4154 solely expressed DEGs were screened out from these two comparison groups. Interestingly, 263 genes showed opposite expression patterns in DT and RWT plants. Among them, genes homologous to Protein Phosphatase 2C (PP2C), Aspartic protease in guard cell 1 (ASPG1), auxin-responsive protein, Arabidopsis pseudo response regualtor 2 (APRR2), GA stimulated transcripts in Arabidopsis 6 (GASA6), Calmodulin-like protein 19 (CML19), abscisic acid 8'-hydroxylases and calcium-transporting ATPase, et al. were related with drought-stress and water stimulus response. Sixteen DEGs involved in starch synthesis, accumulation and tuber formation exhibited significantly different expression upon re-watering. In addition, 1630, 1527 and 1596 transcription factor encoding genes were detected in CT, DT and RWT. DEGs of ERF, bHLH, MYB, NAC, WRKY, C2H2, bZIP and HD-ZIP families accounted for 50% in three comparison groups, respectively. Furthermore, characteristics of 565 gene ontology (GO) and 108 Kyoto Encyclopedia of Genes and Genomes pathways (KEGG) were analyzed with the 4154 DEGs. All these results suggest that the drought- and water-stimulus response could be implemented by the regulated expression of metabolic pathway DEGs, and these genes were involved in the endogenous hormone biosynthesis and signal transduction pathways. Our data provide more direct information for future study on the interaction between key genes involved in various metabolic pathways under drought stress in potato.

Combining Next Generation Sequencing with Bulked Segregant Analysis to Fine Map a Stem Moisture Locus in Sorghum (Sorghum bicolor L. Moench)

Sorghum is one of the most promising bioenergy crops. Stem juice yield, together with stem sugar concentration, determines sugar yield in sweet sorghum. Bulked segregant analysis (BSA) is a gene mapping technique for identifying genomic regions containing genetic loci affecting a trait of interest that when combined with deep sequencing could effectively accelerate the gene mapping process. In this study, a dry stem sorghum landrace was characterized and the stem water controlling locus, qSW6, was fine mapped using QTL analysis and the combined BSA and deep sequencing technologies. Results showed that: (i) In sorghum variety Jiliang 2, stem water content was around 80% before flowering stage. It dropped to 75% during grain filling with little difference between different internodes. In landrace G21, stem water content keeps dropping after the flag leaf stage. The drop from 71% at flowering time progressed to 60% at grain filling time. Large differences exist between different internodes with the lowest (51%) at the 7^th and 8^th internodes at dough stage. (ii) A quantitative trait locus (QTL) controlling stem water content mapped on chromosome 6 between SSR markers Ch6-2 and gpsb069 explained about 34.7-56.9% of the phenotypic variation for the 5^th to 10^th internodes, respectively. (iii) BSA and deep sequencing analysis narrowed the associated region to 339 kb containing 38 putative genes. The results could help reveal molecular mechanisms underlying juice yield of sorghum and thus to improve total sugar yield.

A genome scan for quantitative trait loci associated with Vibrio anguillarum infection resistance in Japanese flounder (Paralichthys olivaceus) by bulked segregant analysis

A recent genetic linkage map was employed to detect quantitative trait loci (QTLs) associated with Vibrio anguillarum resistance in Japanese flounder. An F1 family established and challenged with V. anguillarum in 2009 was used for QTL mapping. Of the 221 simple sequence repeat (SSR) markers used to detect polymorphisms in the parents of F1, 170 were confirmed to be polymorphic. The average distance between the markers was 10.6 cM. Equal amounts of genomic DNA from 15 fry that died early and from 15 survivors were pooled separately to constitute susceptible bulk and resistance bulk DNA. Bulked segregant analysis and QTL mapping were combined to detect candidate SSR markers and regions associated with the disease. A genome scan identified four polymorphic SSR markers, two of which were significantly different between susceptible and resistance bulk (P=0.008). These two markers were located in linkage group (LG) 7; therefore, all the SSR markers in LG7 were genotyped in all the challenged fry by single marker analysis. Using two different models, 11-17 SSR markers were detected with different levels of significance. To confirm the associations of these markers with the disease, composite interval mapping was employed to genotype all the challenged individuals. One and three QTLs, which explained more than 60 % of the phenotypic variance, were detected by the two models. Two of the QTLs were located at 48.6 cM. The common QTL may therefore be a major candidate region for disease resistance against V. anguillarum infection.

Using next-generation sequencing to isolate mutant genes from forward genetic screens

The long-lasting success of forward genetic screens relies on the simple molecular basis of the characterized phenotypes, which are typically caused by mutations in single genes. Mapping the location of causal mutations using genetic crosses has traditionally been a complex, multistep procedure, but next-generation sequencing now allows the rapid identification of causal mutations at single-nucleotide resolution even in complex genetic backgrounds. Recent advances of this mapping-by-sequencing approach include methods that are independent of reference genome sequences, genetic crosses and any kind of linkage information, which make forward genetics amenable for species that have not been considered for forward genetic screens so far.

Early embryonic androgen exposure induces transgenerational epigenetic and metabolic changes

Androgen excess is a central feature of polycystic ovary syndrome (PCOS), which affects 6% to 10% of young women. Mammals exposed to elevated androgens in utero develop PCOS-like phenotypes in adulthood, suggesting fetal origins of PCOS. We hypothesize that excess androgen exposure during early embryonic development may disturb the epigenome and disrupt metabolism in exposed and unexposed subsequent generations. Zebrafish were used to study the underlying mechanism of fetal origins. Embryos were exposed to androgens (testosterone and dihydrotestosterone) early at 26 to 56 hours post fertilization or late at 21 to 28 days post fertilization. Exposed zebrafish (F0) were grown to adults and crossed to generate unexposed offspring (F1). For both generations, global DNA methylation levels were examined in ovaries using a luminometric methylation assay, and fasting and postprandial blood glucose levels were measured. We found that early but not late androgen exposure induced changes in global methylation and glucose homeostasis in both generations. In general, F0 adult zebrafish exhibited altered global methylation levels in the ovary; F1 zebrafish had global hypomethylation. Fasting blood glucose levels were decreased in F0 but increased in F1; postprandial glucose levels were elevated in both F0 and F1. This androgenized zebrafish study suggests that transient excess androgen exposure during early development can result in transgenerational alterations in the ovarian epigenome and glucose homeostasis. Current data cannot establish a causal relationship between epigenetic changes and altered glucose homeostasis. Whether transgenerational epigenetic alteration induced by prenatal androgen exposure plays a role in the development of PCOS in humans deserves study.

QTL-seq identifies an early flowering QTL located near Flowering Locus T in cucumber

Next-generation sequencing enabled a fast discovery of a major QTL controlling early flowering in cucumber, corresponding to the FT gene conditioning flowering time in Arabidopsis. Next-generation sequencing technologies are making it faster and more efficient to establish the association of agronomic traits with molecular markers or candidate genes, which is the requirement for marker-assisted selection in molecular breeding. Early flowering is an important agronomic trait in cucumber (Cucumis sativus L.), but the underlying genetic mechanism is unknown. In this study, we identified a candidate gene for early flowering QTL, Ef1.1 through QTL-seq. Segregation analysis in F2 and BC1 populations derived from a cross between two inbred lines "Muromskij" (early flowering) and "9930" (late flowering) suggested quantitative nature of flowering time in cucumber. Genome-wide comparison of SNP profiles between the early and late-flowering bulks constructed from F2 plants identified a major QTL, designated Ef1.1 on cucumber chromosome 1 for early flowering in Muromskij, which was confirmed by microsatellite marker-based classical QTL mapping in the F2 population. Joint QTL-seq and traditional QTL analysis delimited Ef1.1 to an 890 kb genomic region. A cucumber gene, Csa1G651710, was identified in this region, which is a homolog of the FLOWERING LOCUS T (FT), the main flowering switch gene in Arabidopsis. Quantitative RT-PCR study of the expression level of Csa1G651710 revealed significantly higher expression in early flowering genotypes. Data presented here provide support for Csa1G651710 as a possible candidate gene for early flowering in the cucumber line Muromskij.