Identification of Genes Potentially Associated with the Fertility Instability of S-Type Cytoplasmic Male Sterility in Maize via Bulked Segregant RNA-Seq

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.

Molecular genetics of cytoplasmic male sterility and restorer-of-fertility for the fine tuning of pollen production in crops

Cytoplasmic male sterility (CMS) is an increasingly important issue within the context of hybrid seed production. Its genetic framework is simple: S-cytoplasm for male sterility induction and dominant allele of the restorer-of-fertility gene (Rf) for suppression of S. However, breeders sometimes encounter a phenotype of CMS plants too complex to be explained via this simple model. The molecular basis of CMS provides clue to the mechanisms that underlie the expression of CMS. Mitochondria have been associated with S, and several unique ORFs to S-mitochondria are thought to be responsible for the induction of male sterility in various crops. Their functions are still the subject of debate, but they have been hypothesized to emit elements that trigger sterility. Rf suppresses the action of S by various mechanisms. Some Rfs, including those that encode the pentatricopeptide repeat (PPR) protein and other proteins, are now considered members of unique gene families that are specific to certain lineages. Additionally, they are thought to be complex loci in which several genes in a haplotype simultaneously counteract an S-cytoplasm and differences in the suite of genes in a haplotype can lead to multiple allelism including strong and weak Rf at phenotypic level. The stability of CMS is influenced by factors such as the environment, cytoplasm, and genetic background; the interaction of these factors is also important. In contrast, unstable CMS becomes inducible CMS if its expression can be controlled. CMS becomes environmentally sensitive in a genotype-dependent manner, suggesting the feasibility of controlling the expression of CMS.

Mechanism and Utilization of Ogura Cytoplasmic Male Sterility in Cruciferae Crops

Hybrid production using lines with cytoplasmic male sterility (CMS) has become an important way to utilize heterosis in vegetables. Ogura CMS, with the advantages of complete pollen abortion, ease of transfer and a progeny sterility rate reaching 100%, is widely used in cruciferous crop breeding. The mapping, cloning, mechanism and application of Ogura CMS and fertility restorer genes in Brassica napus, Brassica rapa, Brassica oleracea and other cruciferous crops are reviewed herein, and the existing problems and future research directions in the application of Ogura CMS are discussed.

PNGSeqR: An R Package for Rapid Candidate Gene Selection through Pooled Next-Generation Sequencing

Although bulked segregant analysis (BSA) has been used extensively in genetic mapping, user-friendly tools which can integrate current algorithms for researchers with no background in bioinformatics are scarce. To address this issue, we developed an R package, PNGSeqR, which takes single-nucleotide polymorphism (SNP) markers from next-generation sequencing (NGS) data in variant call format (VCF) as the input file, provides four BSA algorithms to indicate the magnitude of genome-wide signals, and rapidly defines the candidate region through the permutation test and fractile quantile. Users can choose the analysis methods according to their data and experimental design. In addition, it also supports differential expression gene analysis (DEG) and gene ontology analysis (GO) to prioritize the target gene. Once the analysis is completed, the plots can conveniently be exported.

Fine mapping and analysis of candidate genes for qFT7.1, a major quantitative trait locus controlling flowering time in Brassica rapa L

qFT7.1, a major QTL for flowering time in Brassica rapa was fine-mapped to chromosome A07 in a 56.4-kb interval, in which the most likely candidate gene is BraA07g018240.3C. In Brassica rapa, flowering time (FT) is an important agronomic trait that affects the yield, quality, and adaption. FT is a complicated trait that is regulated by many genes and is affected greatly by the environment. In this study, a chromosome segment substitution line (CSSL), CSSL16, was selected that showed later flowering than the recurrent parent, a rapid-cycling inbred line of B. rapa (RcBr). Using Bulked Segregant RNA sequencing, we identified a late flowering quantitative trait locus (QTL), designated as qFT7.1, on chromosome A07, based on a secondary-F₂ population derived from the cross between CSSL16 and RcBr. qFT7.1 was further validated by conventional QTL mapping. This QTL explained 39.9% (logarithm of odds = 32.2) of the phenotypic variations and was fine mapped to a 56.4-kb interval using recombinant analysis. Expression analysis suggested that BraA07g018240.3C, which is homologous to ATC (encoding Arabidopsis thaliana CENTRORADIALIS homologue), a gene for delayed flowering in Arabidopsis, as the most promising candidate gene. Sequence analysis demonstrated that two synonymous mutations existed in the coding region and numerous bases replacements existed in promoter region between BraA07g018240.3C from CSSL16 and RcBr. The results will increase our knowledge related to the molecular mechanism of late flowering in B. rapa and lays a solid foundation for the breeding of late bolting B. rapa.

Conjunctive Analyses of Bulk Segregant Analysis Sequencing and Bulk Segregant RNA Sequencing to Identify Candidate Genes Controlling Spikelet Sterility of Foxtail Millet

Foxtail millet has gradually become a model gramineous C4 crop owing to its short growth period and small genome. Research on the development of its spikelets is not only directly related to the yield and economic value of foxtail millet but also can provide a reference for studying the fertility of other C4 crops. In this study, a hybrid population containing 200 offspring was constructed from the Xinong8852 and An15 parental lines, and two extreme trait populations were constructed from the F2 generation for analysis of the spikelet sterility. The F2 population conformed to a 3:1 Mendelian segregation ratio, and it was thus concluded that this trait is likely controlled by a single recessive gene. Bulk segregant analysis sequencing (BSA-Seq) was used to determine the candidate regions and candidate genes related to the development of foxtail millet spikelets. Additionally, the functional analysis of differentially expressed genes in populations with different traits was conducted by bulk segregant RNA sequencing (BSR-Seq). Finally, conjunctive analysis of BSA-Seq and BSR-Seq results, combined with biological information analysis, revealed six genes on chromosome VII that were ultimately identified as candidate genes controlling foxtail millet spikelet development. This study provides a new reference for research on foxtail millet sterility and lays a solid foundation for the examination of fertility in other gramineous crops.

Allelically and Differentially Expressed Genes After Infection of Edwardsiella ictaluri in Channel Catfish as Determined by Bulk Segregant RNA-Seq

Identification of genetic markers associated with resistance against enteric septicemia of catfish (ESC) is of great interest for genetic enhancement programs of catfish. In the present study, bulk segregant RNA-Seq analysis was applied to determine differentially expressed genes and alleles after ESC infection. Here we report three genomic regions on LG1, LG12, and LG26, containing significant single-nucleotide polymorphisms (SNPs). These genomic regions aligned well with quantitative trait loci (QTL) previously identified. Within the QTL regions, eleven genes were found to be differentially regulated between phenotypic bulks. Importantly, the QTL on linkage group 1 (LG1) were found to be expressed in the liver, whereas the QTL on LG12 and LG26 were expressed in the intestine, suggesting multiple mechanisms of ESC resistance. It is apparent that apolipoproteins may be important for ESC resistance as the QTL on LG1 included the 14-kDa apolipoprotein genes that are both allelically expressed and differentially expressed between the resistant and susceptible bulks. Traf2 and NCK-interacting protein kinase (TNIK) were found in the QTL on LG12, and it was downregulated in resistant fish, suggesting the importance of NCK downregulation in ESC resistance, as previously reported. In addition, we observed divergent gene expression patterns between the liver and intestine after infection. Immune/inflammatory-related processes were overrepresented from liver DEGs, while those DEGs identified from intestine were enriched for proteolysis and wounding processes. Taken together, the BSR-Seq analysis presented here advanced the knowledge of ESC resistance, providing information of not only positions of QTL but also genes and their differential expression between resistant and susceptible fish, making it one step closer to the identification of the causal genes for ESC resistance.

Bulk segregation analysis in the NGS era: a review of its teenage years

Bulk segregation analysis (BSA) utilizes a strategy of pooling individuals with extreme phenotypes to conduct economical and rapidly linked marker screening or quantitative trait locus (QTL) mapping. With the development of next-generation sequencing (NGS) technology in the past 10 years, BSA methods and technical systems have been gradually developed and improved. At the same time, the ever-decreasing costs of sequencing accelerate NGS-based BSA application in different species, including eukaryotic yeast, grain crops, economic crops, horticultural crops, trees, aquatic animals, and insects. This paper provides a landscape of BSA methods and reviews the BSA development process in the past decade, including the sequencing method for BSA, different populations, different mapping algorithms, associated region threshold determination, and factors affecting BSA mapping. Finally, we summarize related strategies in QTL fine mapping combining BSA.

Development of genic KASP SNP markers from RNA-Seq data for map-based cloning and marker-assisted selection in maize

BACKGROUND

Maize is one of the most important field crops in the world. Most of the key agronomic traits, including yield traits and plant architecture traits, are quantitative. Fine mapping of genes/ quantitative trait loci (QTL) influencing a key trait is essential for marker-assisted selection (MAS) in maize breeding. However, the SNP markers with high density and high polymorphism are lacking, especially kompetitive allele specific PCR (KASP) SNP markers that can be used for automatic genotyping. To date, a large volume of sequencing data has been produced by the next generation sequencing technology, which provides a good pool of SNP loci for development of SNP markers. In this study, we carried out a multi-step screening method to identify kompetitive allele specific PCR (KASP) SNP markers based on the RNA-Seq data sets of 368 maize inbred lines.

RESULTS

A total of 2,948,985 SNPs were identified in the high-throughput RNA-Seq data sets with the average density of 1.4 SNP/kb. Of these, 71,311 KASP SNP markers (the average density of 34 KASP SNP/Mb) were developed based on the strict criteria: unique genomic region, bi-allelic, polymorphism information content (PIC) value ≥0.4, and conserved primer sequences, and were mapped on 16,161 genes. These 16,161 genes were annotated to 52 gene ontology (GO) terms, including most of primary and secondary metabolic pathways. Subsequently, the 50 KASP SNP markers with the PIC values ranging from 0.14 to 0.5 in 368 RNA-Seq data sets and with polymorphism between the maize inbred lines 1212 and B73 in in silico analysis were selected to experimentally validate the accuracy and polymorphism of SNPs, resulted in 46 SNPs (92.00%) showed polymorphism between the maize inbred lines 1212 and B73. Moreover, these 46 polymorphic SNPs were utilized to genotype the other 20 maize inbred lines, with all 46 SNPs showing polymorphism in the 20 maize inbred lines, and the PIC value of each SNP was 0.11 to 0.50 with an average of 0.35. The results suggested that the KASP SNP markers developed in this study were accurate and polymorphic.

CONCLUSIONS

These high-density polymorphic KASP SNP markers will be a valuable resource for map-based cloning of QTL/genes and marker-assisted selection in maize. Furthermore, the method used to develop SNP markers in maize can also be applied in other species.

Identification of fertility-related genes for maize CMS-S via Bulked Segregant RNA-Seq

Cytoplasmic male sterility (CMS) is extensively used in maize hybrid production, and identification of genes related to fertility restoration for CMS is important for hybrid breeding. The fertility restoration of S type CMS is governed by several loci with major and minor effects, while the mechanism of fertility restoration for CMS-S is still unknown. In this study, BSR-Seq was conducted with two backcrossing populations with the fertility restoration genes, Rf3 and Rf10, respectively. Genetic mapping via BSR-Seq verified the positions of the two loci. A total of 353 and 176 differentially expressed genes (DEGs) between the male fertility and male sterile pools were identified in the populations with Rf3 and Rf10, respectively. In total, 265 DEGs were co-expressed in the two populations, which were up-regulated in the fertile plants, and they might be related to male fertility involving in anther or pollen development. Moreover, 35 and seven DEGs were specifically up-regulated in the fertile plants of the population with Rf3 and Rf10, respectively. Function analysis of these DEGs revealed that jasmonic acid (JA) signal pathway might be involved in the Rf3 mediated fertility restoration for CMS-S, while the small ubiquitin-related modifier system could play a role in the fertility restoration of Rf10.

Role of BrSDG8 on bolting in Chinese cabbage (Brassica rapa)

Mapping and resequencing of two allelic early bolting mutants ebm5-1 and ebm5-2 revealed that the BrSDG8 gene is related to bolting in Chinese cabbage (Brassica rapa ssp. pekinensis). Bolting influences the leafy head formation and seed yield of Chinese cabbage therefore being an important agronomic trait. Herein, two allelic early bolting mutants, ebm5-1 and ebm5-2, stably inherited in Chinese cabbage were obtained from wild-type 'FT' seeds by ethyl methane sulfonate mutagenesis. Both mutants flowered significantly earlier than 'FT,' and genetic analysis revealed that the early bolting of the two mutants was controlled by one recessive nuclear gene. With BSR-seq, the mutations originating lines ebm5-1 and ebm5-2 were located to the same region in chromosome A07. Using the 1741 F₂ individuals with the ebm5-1 phenotype as the mapping population, this region was narrowed to 56.24 kb between markers InDel18 and InDel45. A single-nucleotide polymorphism (SNP) was aligned to the BraA07g040740.3C (BrSDG8) region by whole-genome resequencing of ebm5-1 mutant and 'FT.' BrSDG8 is a homolog of Arabidopsis thaliana SDG8 encoding a histone methyltransferase affecting H3K4 trimethylation in FLOWERING LOCUS C chromatin. Comparative sequencing established that the SNP occurred on BrSDG8 17th exon in ebm5-1. Genotype analysis showed full co-segregation of the early bolting phenotype with this SNP. Cloning of allelic mutant ebm5-2 indicated that it harbors a deletion mutation on the 12th exon of BrSDG8. Quantitative real-time PCR analysis indicated that BrSDG8 expression level was observably lower in mutant ebm5-1 than in 'FT.' Overall, the present results provide strong evidence that BrSDG8 mutation leads to early bolting in Chinese cabbage, thereby providing a basis to understand the molecular mechanisms underlying this phenotype.

Fine mapping of a leaf flattening gene Bralcm through BSR-Seq in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

Leaf flattening influences plant photosynthesis, thereby affecting product yield and quality. Here, we obtained a stably inherited leaf crinkled mutant (lcm), derived from the Chinese cabbage doubled haploid (DH) ‘FT’ line using EMS mutagenesis combined with isolated microspore culture. The crinkled phenotype was controlled by a single recessive nuclear gene, namely Bralcm, which was preliminarily mapped to chromosome A01 by bulked segregant analysis RNA-seq, and further between markers SSRS-1 and IndelD-20 using 1,575 recessive homozygous individuals in F₂ population by a map-based cloning method. The target region physical distance was 126.69 kb, containing 23 genes; the marker SSRMG-4 co-segregated with the crinkled trait. Further, we found SSRMG-4 to be located on BraA01g007510.3C, a homolog of AHA2, which encodes H⁺-ATPase2, an essential enzyme in plant growth and development. Sequence analysis indicated a C to T transition in exon 7 of BraA01g007510.3C, resulting in a Thr (ACT) to Ile (ATT) amino acid change. Genotyping revealed that the leaf crinkled phenotype fully co-segregated with this SNP within the recombinants. qRT-PCR demonstrated that BraA01g007510.3C expression in lcm mutant leaves was dramatically higher than that in wild-type ‘FT’. Thus, BraA01g007510.3C is a strong candidate gene for Bralcm, and AHA2 is possibly associated with leaf flattening in Chinese cabbage.

Analysis of differentially expressed genes and pathways associated with male sterility lines in watermelon via bulked segregant RNA-seq

Genic male sterility (GMS) is a common and important trait, which is widely used for the production of hybrid seeds. However, the molecular mechanism of GMS in watermelon remains poorly understood. In this study, we comparatively analyzed the transcriptome profiles of sterile and fertile floral buds using the bulked segregant analysis (BSA) and transcriptome sequencing (RNA-seq). A total of 2507 differentially expressed genes (DEGs) including 593 up-regulated and 1914 down-regulated, were identified to be related to male sterility in watermelon line Se18. Gene ontology (GO) analysis showed that 57 GO terms were significantly enriched, while Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis revealed plant hormone signal transduction, glycolysis/gluconeogenesis, starch and sucrose metabolism, plant-pathogen interaction, phenylpropanoid biosynthesis pathways were obviously enriched. Furthermore, the efficiency of the RNA-seq analysis was validated by quantitative real-time PCR (qRT-PCR). Among the DEGs, some valuable candidate genes involved in pollen development were identified, such as gene Cla000029, a bHLH transcription factor and homologous to MS1 in Arabidopsis. Moreover, other DEGs including MYB gene Cla012590 (MYB26), Cla017100 (MYB21), etc., also provide useful information for further understanding the function of key genes involved in pollen development. This study provides new insights into the global network of male sterility in watermelon.

Identification of two recessive etiolation genes (py1, py2) in pakchoi (Brassica rapa L. ssp. chinensis)

BACKGROUND

Leaf color is a major agronomic trait, which has a strong influence on crop yields. Isolating leaf color mutants can represent valuable materials for research in chlorophyll (Chl) biosynthesis and metabolism regulation.

RESULTS

In this study, we identified a stably inherited yellow leaf mutant derived from 'Huaguan' pakchoi variety via isolated microspore culture and designated as pylm. This mutant displayed yellow leaves after germination. Its etiolated phenotype was nonlethal and stable during the whole growth period. Its growth was weak and its hypocotyls were markedly elongated. Genetic analysis revealed that two recessive nuclear genes, named py1 and py2, are responsible for the etiolation phenotype. Bulked segregant RNA sequencing (BSR-Seq) showed that py1 and py2 were mapped on chromosomes A09 and A07, respectively. The genes were single Mendelian factors in F_3:4 populations based on a 3:1 phenotypic segregation ratio. The py1 was localized to a 258.3-kb interval on a 34-gene genome. The differentially expressed gene BraA09004189 was detected in the py1 mapping region and regulated heme catabolism. One single-nucleotide polymorphism (SNP) of BraA09004189 occurred in pylm. A candidate gene-specific SNP marker in 1520 F_3:4 yellow-colored individuals co-segregated with py1. For py2, 1860 recessive homozygous F_3:4 individuals were investigated and localized py2 to a 4.4-kb interval. Of the five genes in this region, BraA07001774 was predicted as a candidate for py2. It encoded an embryo defective 1187 and a phosphotransferase related to chlorophyll deficiency and hypocotyl elongation. One SNP of BraA07001774 occurred in pylm. It caused a single amino acid mutation from Asp to Asn. According to quantitative real-time polymerase chain reaction (qRT-PCR), BraA07001774 was downregulated in pylm.

CONCLUSIONS

Our study identified a Chl deficiency mutant pylm in pakchoi. Two recessive nuclear genes named py1 and py2 had a significant effect on etiolation. Candidate genes regulating etiolation were identified as BraA09004189 and BraA07001774, respectively. These findings will elucidate chlorophyll metabolism and the molecular mechanisms of the gene interactions controlling pakchoi etiolation.

Construction of restorer lines and molecular mapping for restorer gene of hau cytoplasmic male sterility in Brassica napus

Successfully constructing restorer lines for the hau CMS line and molecular mapping of Rfh to a 94 kb candidate region on chromosome A03 in Brassica napus. Cytoplasmic male sterility is a general phenomenon in almost 200 species, and the interaction between chimeric genes in mitochondria and restorer genes in nucleus may be responsible for restoration of male fertility. Orf288 has been identified as a CMS-associated gene in the hau CMS line of Brassica napus and Brassica juncea; however, the restorer lines/genes have not been found yet. We therefore have successfully constructed two restorer lines in B. napus by extensive testcrossing and have mapped a major restorer gene Rfh to a physical distance of 94 kb on chromosome A03 by whole-genome resequencing and molecular markers. We found that the restorer line is indeed restored to male fertility at histological level. Comparative genomics and collinearity analysis between close relatives revealed that rearrangements and recombination may have happened and thus caused the production of Rfh or components of the restoration of fertility complex. Meanwhile, nuclear backgrounds with multiple loci and temperature were related to the variation and instability of restoration of fertility in three different populations. Our study provides new sights into the coevolution between restorer genes and CMS-associated genes as well as the cultivation of superior hybrids via molecular breeding.

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.

RNA Sequencing-Based Bulked Segregant Analysis Facilitates Efficient D-genome Marker Development for a Specific Chromosomal Region of Synthetic Hexaploid Wheat

Common wheat originated from interspecific hybridization between cultivated tetraploid wheat and its wild diploid relative Aegilops tauschii followed by amphidiploidization. This evolutionary process can be reproduced artificially, resulting in synthetic hexaploid wheat lines. Here we performed RNA sequencing (RNA-seq)-based bulked segregant analysis (BSA) using a bi-parental mapping population of two synthetic hexaploid wheat lines that shared identical A and B genomes but included with D-genomes of distinct origins. This analysis permitted identification of D-genome-specific polymorphisms around the Net2 gene, a causative locus to hybrid necrosis. The resulting single nucleotide polymorphisms (SNPs) were classified into homoeologous polymorphisms and D-genome allelic variations, based on the RNA-seq results of a parental tetraploid and two Ae. tauschii accessions. The difference in allele frequency at the D-genome-specific SNP sites between the contrasting bulks (ΔSNP-index) was higher on the target chromosome than on the other chromosomes. Several SNPs with the highest ΔSNP-indices were converted into molecular markers and assigned to the Net2 chromosomal region. These results indicated that RNA-seq-based BSA can be applied efficiently to a synthetic hexaploid wheat population to permit molecular marker development in a specific chromosomal region of the D genome.

Identification of the dwarf gene GmDW1 in soybean (Glycine max L.) by combining mapping-by-sequencing and linkage analysis

KEY MESSAGE

GmDW1 encodes an ent-kaurene synthase (KS) acting at the early step of the biosynthesis pathway for gibberellins (GAs) and regulates the development of plant height in soybean. Plant height is an important component of plant architecture, and significantly affects crop breeding practices and yield. Here, we report the characterization of an EMS-induced dwarf mutant (dw) of the soybean cultivar Zhongpin 661 (ZDD23893). The dw mutant displayed reduced plant height and shortened internodes, both of which were mainly attributed to the longitudinally decreased cell length. The bioactive GA1 (gibberellin A1) and GA4 (gibberellin A4) were not detectable in the stem of dw, and the dwarf phenotype could be rescued by treatment with exogenous GA3. Genetic analysis showed that the dwarf trait of dw was controlled by a recessive nuclear gene. By combining linkage analysis and mapping-by-sequencing, we mapped the GmDW1 gene to an approximately 460-kb region on chromosome (Chr.) 8, containing 36 annotated genes in the reference Willliams 82 genome. Of these genes, we identified two nonsynonymous single nucleotide polymorphisms (SNPs) that are present in the encoding regions of Gmdw1 and Glyma.08G165100 in dw, respectively. However, only the SNP mutation (T>A) at nucleotide 1224 in Gmdw1 cosegregated with the dwarf phenotype. GmDW1 encodes an ent-kaurene synthase, and was expressed in various tissues including root, stem, and leaf. Further phenotypic analysis of the allelic variations in soybean accessions strongly indicated that GmDW1 is responsible for the dwarf phenotype in dw. Our results provide important information for improving our understanding of the genetics of soybean plant height and crop breeding.

Morphological, transcriptomics and biochemical characterization of new dwarf mutant of Brassica napus

Plant height is a key trait of plant architecture, and is responsible for both yield and lodging resistance in Brassica napus. A dwarf mutant line (bnaC.dwf) was obtained by chemical mutagenesis of an inbred line T6. However, the molecular mechanisms and changed biological processes of the dwarf mutant remain to be determined. In this study, a comparative transcriptome analysis between bnaC.dwf and T6 plants was performed to identify genome-wide differentially expressed genes (DEGs) and possible biological processes that may explain the phenotype variations in bnaC.dwf. As a result of this analysis, 60,134,746-60,301,384 clean reads were aligned to 60,074 genes in the B. napus genome, and accounted for 60.03% of the annotated genes. In total, 819 differentially expressed genes were used for GO (Gene Ontology) term and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analyses with a FDR (false discovery rate) criterion of <0.001, |log₂Ratio| ≥ 1. We focused on plant hormone signal transduction pathways, plant-pathogen interaction pathway, protein phosphorylation and degradation pathways and sugar metabolism pathways. Taken together, the decrease in local auxin (IAA) levels, the variation in BnTCH4, BnKAN1, BnERF109, COI1-JAZ9-MYC2, auxin response genes (BnGH3.11, BnSAUR78, and AUX/IAA19), and ABA (abscisic acid) signaling genes (BnADP5, BnSnRK2.1, BnABF3.1) partially accounted for variations of cell proliferation in internodes, shoot and root apical meristem maintenance, abiotic and biotic stress resistance, and pre-harvest sprouting. As a comprehensive consequence of the cross-talk between plant hormones, sugar metabolism, plant-pathogen interactions and protein metabolism, bnaC.dwf presents distinct phenotypes from T6. These results will be helpful for shedding light on molecular mechanisms in the dwarf mutant, and give insight into further molecular breeding of semi-dwarf B. napus.

Advantages of Amplifluor-like SNP markers over KASP in plant genotyping

BACKGROUND

KASP (KBioscience Competitive Allele Specific PCR) and Amplifluor (Amplification with fluorescence) SNP markers are two prominent technologies based upon a shared identical Allele-specific PCR platform.

METHODS

Amplifluor-like SNP and KASP analysis was carried out using published and own design of Universal probes (UPs) and Gene-specific primers (GSPs).

RESULTS

Advantages of the Amplifluor-like system over KASP include the significantly lower costs and much greater flexibility in the adjustment and development of 'self-designed' dual fluorescently-labelled UPs and regular GSPs. The presented results include optimisation of 'tail' length in UPs and GSPs, protocol adjustment, and the use of various fluorophores in different qPCR instruments. The compatibility of the KASP Master-mix in both original and Amplifluor-like systems has been demonstrated in the presented results, proving their similar principles. Results of SNP scoring with rare alleles in addition to more frequently occurring alleles are shown.

CONCLUSIONS

The Amplifluor-like system produces SNP genotyping results with a level of sensitivity and accuracy comparable to KASP but at a significantly cheaper cost and with much greater flexibility for UPs with self-designed GSPs.

Mapping and Preliminary Analysis of ABORTED MICROSPORES (AMS) as the Candidate Gene Underlying the Male Sterility (MS-5) Mutant in Melon (Cucumis melo L.)

Melon is an important agricultural and economic vegetable crop worldwide. The genetic male sterility mutant (ms-5) has a recessive nuclear gene that controls the male sterility germplasm. Male sterility could reduce the cost of F₁ seed production in melon, but heterozygous fertile plants should be removed before pollination. In this study, bulked segregant analysis combined with specific length amplified fragment sequencing was applied to map the single nuclear male sterility recessive gene. A 30-kb candidate region on chromosome 9 located on scaffold 000048 and spanning 2,522,791 to 2,555,104 bp was identified and further confirmed by cleavage amplified polymorphic sequence markers based on parental line resequencing data and classical mapping of 252 F₂ individuals. Gene prediction indicated that six annotated genes are present in the 30-kb candidate region. Quantitative RT-PCR revealed significant differences in the expression level of the LOC103498166 ABORTED MICROSPORES (AMS) gene in male-sterile lines (ms-5) and male-fertile (HM1-1) lines during the 2-mm (tetrad) and 5-mm (the first pollen mitosis) periods, and negative regulation of the AMS candidate gene transcription factor was also detected. Sequencing and cluster analysis of the AMS transcription factor revealed five single-nucleotide polymorphisms between the parental lines. The data presented herein suggest that the AMS transcription factor is a possible candidate gene for single nuclear male sterility in melon. The results of this study will help breeders to identify male-sterile and -fertile plants at seeding as marker-assisted selection methods, which would reduce the cost of seed production and improve the use of male-sterile lines in melon.