A high-resolution linkage map of the Rfd1, a restorer-of-fertility locus for cytoplasmic male sterility in radish (Raphanus sativus L.) produced by a combination of bulked segregant analysis and RNA-Seq

An autopolyploid-suitable polyBSA-seq strategy for screening candidate genetic markers linked to leaf blight resistance in sugarcane

An autopolyploid-suitable polyBSA-seq strategy was developed for screening candidate genetic markers linked to leaf blight resistance in sugarcane. Due to the complex genome architecture, the quantitative trait loci mappings and linkage marker selections for agronomic traits of autopolyploid crops were mainly limited to the time-consuming and cost intensive construction of genetic maps. To map resistance-linked markers for sugarcane leaf blight (SLB) caused by Stagonospora tainanensis, the autopolyploid-suitable bulk-segregant analysis based on the sequencing (polyBSA-seq) strategy was successfully applied for the first time. Resistant- and susceptible-bulks (R- and S-bulks) constructed from the extreme-phenotypic sugarcane F₁ lines of YT93-159 × ROC22 were deep sequenced with 195.0 × for bulks and 74.4 × for parents. Informative single-dose variants (ISDVs) present as one copy in one parent and null in the other parent were detected based on the genome sequence of LA Purple, an autooctoploid Saccharum officinarum, to screen candidate linkage markers (CLMs). The proportion of the number of short reads harboring ISDVs in the total short reads covering a given genomic position was defined as ISDV index and the ISDVs with indices met the threshold set in this study (0.04-0.14) were selected as CLMs. In total, three resistance- and one susceptibility-related CLMs for SLB resistance were identified by the polyBSA-seq. Among them, two markers on chromosome 10 were less than 300 Kb apart. Furthermore, the RNA-seq was used to calculate the expression level of genes within 1.0 Mb from the aforementioned four CLMs, which demonstrated that twelve genes were differentially expressed between resistant and susceptible clones, including a receptor-like kinase and an ethylene-responsive transcription factor. This is the first reported polyBSA-seq in autopolyploid sugarcane, which specifically tailored for the fast selection of the CLMs and causal genes associated with important agronomic traits.

An insertion mutation located on putative enhancer regions of the MYB26-like gene induces inhibition of anther dehiscence resulting in novel genic male sterility in radish (Raphanus sativus L.)

A novel male-sterility trait was identified in a radish (Raphanus sativus L.) population. Although the size of male-sterile anthers was comparable to that of normal flowers, no pollen grain was observed during anther dehiscence. However, dissection of male-sterile anthers revealed an abundance of normal pollen grains. Analysis of segregating populations showed that a single recessive locus, designated RsMs1, conferred male sterility. Based on two radish draft genome sequences, molecular markers were developed to delimit the genomic region harboring the RsMs1. The region was narrowed down to approximately 24 kb after analyzing recombinants selected from 7511 individuals of a segregating population. Sequencing of the delimited region yielded six putative genes including four genes expressed in the floral tissue, and one gene with significant differential expression between male-fertile and male-sterile individuals of a segregating population. This differentially expressed gene was orthologous to the Arabidopsis MYB26 gene, which played a critical role in anther dehiscence. Excluding a synonymous single nucleotide polymorphism in exon3, no polymorphism involving coding and putative promoter regions was detected between alleles. A 955-bp insertion was identified 7.5 kb upstream of the recessive allele. Highly conserved motifs among four Brassicaceae species were identified around this insertion site, suggesting the presence of putative enhancer sequences. A functional marker was developed for genotyping of the RsMs1 based on the 955-bp insertion. A total of 120 PI accessions were analyzed using this marker, and 11 accessions were shown to carry the recessive rsms1 allele.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01254-9.

Genome-based breeding approaches in major vegetable crops

Vegetable crops are major nutrient sources for humanity and have been well-cultivated since thousands of years of domestication. With the rapid development of next-generation sequencing and high-throughput genotyping technologies, the reference genome of more than 20 vegetables have been well-assembled and published. Resequencing approaches on large-scale germplasm resources have clarified the domestication and improvement of vegetable crops by human selection; its application on genetic mapping and quantitative trait locus analysis has led to the discovery of key genes and molecular markers linked to important traits in vegetables. Moreover, genome-based breeding has been utilized in many vegetable crops, including Solanaceae, Cucurbitaceae, Cruciferae, and other families, thereby promoting molecular breeding at a single-nucleotide level. Thus, genome-wide SNP markers have been widely used, and high-throughput genotyping techniques have become one of the most essential methods in vegetable breeding. With the popularization of gene editing technology research on vegetable crops, breeding efficiency can be rapidly increased, especially by combining the genomic and variomic information of vegetable crops. This review outlines the present genome-based breeding approaches used for major vegetable crops to provide insights into next-generation molecular breeding for the increasing global population.

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.

Beyond balancing selection: frequent mitochondrial recombination contributes to high-female frequencies in gynodioecious Lobelia siphilitica (Campanulaceae)

Gynodioecy is a sexual system in which females and hermaphrodites co-occur. In most gynodioecious angiosperms, sex is determined by an interaction between mitochondrial male-sterility genes (CMS) that arise via recombination and nuclear restorer alleles that evolve to suppress them. In theory, gynodioecy occurs when multiple CMS types are maintained at equilibrium frequencies by balancing selection. However, some gynodioecious populations contain very high frequencies of females. High female frequencies are not expected under balancing selection, but could be explained by the repeated introduction of novel CMS types. To test for balancing selection and/or the repeated introduction of novel CMS, we characterised cytoplasmic haplotypes from 61 populations of Lobelia siphilitica that vary widely in female frequency. We confirmed that mitotype diversity and female frequency were positively correlated across populations, consistent with balancing selection. However, while low-female populations hosted mostly common mitotypes, high-female populations and female plants hosted mostly rare, recombinant mitotypes likely to carry novel CMS types. Our results suggest that balancing selection maintains established CMS types across this species, but extreme female frequencies result from frequent invasion by novel CMS types. We conclude that balancing selection alone cannot account for extreme population sex-ratio variation within a gynodioecious species.

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.

Genetic mapping reveals a marker for yellow skin in watermelon (Citrullus lanatus L.)

As a diverse species, watermelon [Citrullus lanatus (Thunb.) Matsum. &Nakai var. lanatus] has different kinds of fruit sizes, shapes, flesh colors and skin colors. Skin color is among the major objectives for breeding. Yellow skin is an important trait in watermelon, but the underlying genetic mechanism is unknown. In this study, we identified a locus for yellow skin through BSA-seq and GWAS. A segregation analysis in F₂ and BC₁ populations derived from a cross of two inbred lines ‘94E1’(yellow skin) and ‘Qingfeng’(green skin) suggested that skin color is a qualitative trait. BSA-seq mapping confirmed the locus in the F₂ population, which was detected on chromosome 4 by GWAS among 330 varieties. Several major markers, namely, 15 CAPS markers, 6 SSR markers and 2 SNP markers, were designed to delimit the region to 59.8 kb region on chromosome 4. Utilizing the two populations consisting of 10 yellow and 10 green skin watermelons, we found a tightly linked functional SNP marker for the yellow skin phenotype. The application of this marker as a selection tool in breeding programs will help to improve the breeder’s ability to make selections at early stages of growth, thus accelerating the breeding program.

Genetic mapping reveals a candidate gene (ClFS1) for fruit shape in watermelon (Citrullus lanatus L.)

A 159 bp deletion in ClFS1 gene encoding IQD protein is responsible for fruit shape in watermelon. Watermelon [Citrullus lanatus (Thunb.) Matsum. & Nakai] is known for its rich diversity in fruit size and shape. Fruit shape has been one of the major objectives of watermelon breeding. However, the candidate genes and the underlying genetic mechanism for such an important trait in watermelon are unknown. In this study, we identified a locus on chromosome 3 of watermelon genome controlling fruit shape. Segregation analysis in F₂ and BC₁ populations derived from a cross between two inbred lines "Duan125" (elongate fruit) and "Zhengzhouzigua" (spherical fruit) suggests that fruit shape of watermelon is controlled by a single locus and elongate fruit (OO) is incompletely dominant to spherical fruit (oo) with the heterozygote (Oo) being oval fruit. GWAS profiles among 315 accessions identified a major locus designated on watermelon chromosome 3, which was confirmed by BSA-seq mapping in the F₂ population. The candidate gene was mapped to a region 46 kb on chromosome 3. There were only four genes present in the corresponding region in the reference genome. Four candidate genes were sequenced in this region, revealing that the CDS of Cla011257 had a 159 bp deletion which resulted in the omission of 53 amino acids in elongate watermelon. An indel marker was derived from the 159 bp deletion to test the F₂ population and 105 watermelon accessions. The results showed that Cla011257 cosegregated with watermelon fruit shape. In addition, the Cla011257 expression was the highest at ovary formation stage. The predicted protein of the Cla011257 gene fitted in IQD protein family which was reported to have association with cell arrays and Ca²⁺-CaM signaling modules. Clear understanding of the genes facilitating the fruit shape along with marker association selection will be an effective way to develop new cultivars.

Genetic mapping of the LOBED LEAF 1 (ClLL1) gene to a 127.6-kb region in watermelon (Citrullus lanatus L.)

The lobed leaf character is a unique morphologic trait in crops, featuring many potential advantages for agricultural productivity. Although the majority of watermelon varieties feature lobed leaves, the genetic factors responsible for lobed leaf formation remain elusive. The F2:3 leaf shape segregating population offers the opportunity to study the underlying mechanism of lobed leaf formation in watermelon. Genetic analysis revealed that a single dominant allele (designated ClLL1) controlled the lobed leaf trait. A large-sized F3:4 population derived from F2:3 individuals was used to map ClLL1. A total of 5,966 reliable SNPs and indels were identified genome-wide via a combination of BSA and RNA-seq. Using the validated SNP and indel markers, the location of ClLL1 was narrowed down to a 127.6-kb region between markers W08314 and W07061, containing 23 putative ORFs. Expression analysis via qRT-PCR revealed differential expression patterns (fold-changes above 2-fold or below 0.5-fold) of three ORFs (ORF3, ORF11, and ORF18) between lobed and non-lobed leaf plants. Based on gene annotation and expression analysis, ORF18 (encoding an uncharacterized protein) and ORF22 (encoding a homeobox-leucine zipper-like protein) were considered as most likely candidate genes. Furthermore, sequence analysis revealed no polymorphisms in cDNA sequences of ORF18; however, two notable deletions were identified in ORF22. This study is the first report to map a leaf shape gene in watermelon and will facilitate cloning and functional characterization of ClLL1 in future studies.

Searching for an Accurate Marker-Based Prediction of an Individual Quantitative Trait in Molecular Plant Breeding

Molecular plant breeding with the aid of molecular markers has played an important role in modern plant breeding over the last two decades. Many marker-based predictions for quantitative traits have been made to enhance parental selection, but the trait prediction accuracy remains generally low, even with the aid of dense, genome-wide SNP markers. To search for more accurate trait-specific prediction with informative SNP markers, we conducted a literature review on the prediction issues in molecular plant breeding and on the applicability of an RNA-Seq technique for developing function-associated specific trait (FAST) SNP markers. To understand whether and how FAST SNP markers could enhance trait prediction, we also performed a theoretical reasoning on the effectiveness of these markers in a trait-specific prediction, and verified the reasoning through computer simulation. To the end, the search yielded an alternative to regular genomic selection with FAST SNP markers that could be explored to achieve more accurate trait-specific prediction. Continuous search for better alternatives is encouraged to enhance marker-based predictions for an individual quantitative trait in molecular plant breeding.

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.