A genome-wide SNP-based genetic map and QTL mapping for agronomic traits in Chinese cabbage

QTL mapping and candidate gene analysis reveal two major loci regulating green leaf color in non-heading Chinese cabbage

KEY MESSAGE

Two major-effect QTL GlcA07.1 and GlcA09.1 for green leaf color were fine mapped into 170.25 kb and 191.41 kb intervals on chromosomes A07 and A09, respectively, and were validated by transcriptome analysis. Non-heading Chinese cabbage (NHCC) is a leafy vegetable with a wide range of green colors. Understanding the genetic mechanism behind broad spectrum of green may facilitate the breeding of high-quality NHCC. Here, we used F2 and F7:8 recombination inbred line (RIL) population from a cross between Wutacai (dark-green) and Erqing (lime-green) to undertake the genetic analysis and quantitative trait locus (QTL) mapping in NHCC. The genetic investigation of the F2 population revealed that the variation of green leaf color was controlled by two recessive genes. Six pigments associated with green leaf color, including total chlorophyll, chlorophyll a, chlorophyll b, total carotenoids, lutein, and carotene were quantified and applied for QTL mapping in the RIL population. A total of 7 QTL were detected across the whole genome. Among them, two major-effect QTL were mapped on chromosomes A07 (GlcA07.1) and A09 (GlcA09.1) corresponding to two QTL identified in the F2 population. The QTL GlcA07.1 and GlcA09.1 were further fine mapped into 170.25 kb and 191.41 kb genomic regions, respectively. By comparing gene expression level and gene annotation, BraC07g023810 and BraC07g023970 were proposed as the best candidates for GlcA07.1, while BraC09g052220 and BraC09g052270 were suggested for GlcA09.1. Two InDel molecular markers (GlcA07.1-BcGUN4 and GlcA09.1-BcSG1) associated with BraC07gA023810 and BraC09g052220 were developed and could effectively identify leaf color in natural NHCC accessions, suggesting their potential for marker-assisted leaf color selection in NHCC breeding.

Genetic and Transcriptome Analysis of Leaf Trichome Development in Chinese Cabbage (Brassica rapa L. subsp. pekinensis) and Molecular Marker Development

Chinese cabbage (Brassica rapa L. subsp. pekinensis) is one of the vegetables with the largest cultivated area in China and has been a great addition to the daily diet of Chinese people. A genetic map has been constructed in our previous study using the F₂ population of two inbred lines of Chinese cabbage, namely "G291" (a hairy line) and "ZHB" (a hairless line), based on which a candidate gene related to trichome traits was identified on chromosome A06 with a phenotypic variance of 47%. A molecular marker was found to co-segregate with the trichome traits of the F₂ population, which is in the 5'-flanking region of BrGL1, and a corresponding patent has been granted (NO. CN 108545775 B). Transcriptome analysis was carried out on the cotyledon, the first true leaf and the leaf closest to each inflorescence of F₂ individuals of "G291 × ZHB" with or without trichomes, respectively. Ten pathways, including 189 DEGs, were identified to be involved in the development of trichomes in Chinese cabbage, which may be specifically related to the development of leaf trichomes. Most of the pathways were related to the biosynthesis of the secondary metabolites, which may help plants to adapt to the ever-changing external environment. DEGs also enriched the "plant-pathogen interaction" pathway, which is consistent with the conclusion that trichomes are related to the disease resistance of plants. Our study provides a basis for future research on the occurrence and development of trichomes in Chinese cabbage.

RAD-Seq-Based High-Density Linkage Maps Construction and Quantitative Trait Loci Mapping of Flowering Time Trait in Alfalfa (Medicago sativa L.)

Alfalfa (Medicago sativa L.) is a perennial forage crop known as the "Queen of Forages." To dissect the genetic mechanism of flowering time (FT) in alfalfa, high-density linkage maps were constructed for both parents of an F1 mapping population derived from a cross between Cangzhou (P1) and ZhongmuNO.1 (P2), consisting of 150 progenies. The FT showed a transgressive segregation pattern in the mapping population. A total of 13,773 single-nucleotide polymorphism markers was obtained by using restriction-site associated DNA sequencing and distributed on 64 linkage groups, with a total length of 3,780.49 and 4,113.45 cM and an average marker interval of 0.58 and 0.59 cM for P1 and P2 parent, respectively. Quantitative trait loci (QTL) analyses were performed using the least square means of each year as well as the best linear unbiased prediction values across 4 years. Sixteen QTLs for FT were detected for P1 and 22 QTLs for P2, accounting for 1.40-16.04% of FT variation. RNA-Seq analysis at three flowering stages identified 5,039, 7,058, and 7,996 genes that were differentially expressed between two parents, respectively. Based on QTL mapping, DEGs analysis, and functional annotation, seven candidate genes associated with flowering time were finally detected. This study discovered QTLs and candidate genes for alfalfa FT, making it a useful resource for breeding studies on this essential crop.

Construction of an Intragenic SSR-Based Linkage Map and QTL Mapping for Agronomic Traits in Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

Chinese cabbage (Brassica rapa L. ssp. pekinensis) is one of the most widely cultivated and economically important vegetables in China. Constructing an effective genetic linkage map and mapping quantitative trait loci (QTLs) related to yield and leafy head morphology is of great importance for molecular breeding of Chinese cabbage. Using two diverse Chinese cabbage inbred lines, ZHB and G291, as parents, an F2 segregating population consisting of 240 individuals was prepared for genetic map construction and phenotype investigation in this study. The two parents are significantly different in both shape and size. Sixteen important agronomic traits of F2 individuals were investigated. A genetic map of 105 intragenic simple sequence repeat (SSR) markers distributed across 10 linkage groups (LGs) was constructed, which was 2034.1 cM in length and had an average inter-locus distance of 21.75 cM. We identified 48 QTLs for the tested important agronomic traits on the studied LGs, with LOD scores of 2.51–12.49, which explained the phenotypic variance of 3.41–26.66%. The QTLs identified in this study will facilitate further genetic analysis and marker-assisted genetic improvement of Chinese cabbage.

An extended KASP-SNP resource for molecular breeding in Chinese cabbage(Brassica rapa L. ssp. pekinensis)

Kompetitive allele-specific PCR (KASP) is a cost-effective single-step SNP genotyping technology, With an objective to enhance the marker repertoire and develop high efficient KASP-SNP markers in Chinese cabbage, we re-sequenced four Chinese cabbage doubled haploid (DH) lines, Y177-47, Y635-10, Y510-1 and Y510-9, and generated a total of more than 38.5 billion clean base pairs. A total of 827,720 SNP loci were identified with an estimated density of 3,217 SNPs/Mb. Further, a total of 387,354 SNPs with at least 30 bp to the next most adjacent SNPs on either side were selected as resource for KASP markers. From this resource, 258 (96.27%) of 268 SNP loci were successfully transformed into KASP-SNP markers using a Roche LightCycler 480-II instrument. Among these markers, 221 (85.66%) were co-dominant markers, 220 (85.27%) were non-synonymous SNPs, and 257 (99.6%) were newly developed markers. In addition, 53 markers were applied for genotyping of 34 Brassica rapa accessions. Cluster analysis separated these 34 accessions into three clusters based on heading types. The millions of SNP loci, a large set of resource for KASP markers, as well as the newly developed KASP markers in this study may facilitate further genetic and molecular breeding studies in Brassica rapa.

Insights on SNP types, detection methods and their utilization in Brassica species: Recent progress and future perspectives

The genus Brassica, family Brassicaceae (Cruciferae), comprises many important species of oil crops, vegetables and medicinal plants including B. rapa, B. oleracea, B. nigra, B. napus, B. juncea, B. carinata. Genomic researches in Brassica species is constrained by polyploidization, mainly due to its complicated genomic structure. However, rapid development of methods for detecting single nucleotide polymorphisms (SNP), such as next generation sequencing and SNP microarray, has accelerated release of reference Brassica species genomes as well as discovery of large numbers and genome-wide SNPs, thus intensifying forward genetics in this genus. In this review, we summarize biological characteristics, classification and various methods for detecting SNPs, focusing on high-throughput techniques. Moreover, we describe the pivotal roles of SNPs in genetic diversity, linkage map construction and QTL mapping, comparative genomics, linkage disequilibrium and genome-wide association studies. These insights are expected to deepen our understanding and guide further advancements in Brassica species research.

Genetic Analysis and Fine Mapping of a Spontaneously Mutated Male Sterility Gene in Brassica rapa ssp. chinensis

Male sterility has been widely used in hybrid seed production in Brassica, but not in B. rapa ssp. chinensis, and genetic models of male sterility for this subspecies are unclear. We discovered a spontaneous mutant in B. rapa ssp. chinensis A series of progeny tests indicated that male sterility in B. rapa ssp. chinensis follows a three-allele model with BrMsa , BrMsb , and BrMsc The male sterility locus has been mapped to chromosome A07 in BC1 and F2 populations through genotyping by sequencing. Fine mapping in a total of 1,590 F2 plants narrowed the male sterility gene BrMs to a 400 kb region, with two SNP markers only 0.3 cM from the gene. Comparative gene mapping shows that the Ms gene in B. rapa ssp. pekinensis is different from the BrMs gene of B. rapa ssp. chinensis, despite that both genes are located on chromosome A07. Interestingly, the DNA sequence orthologous to a male sterile gene in Brassica napus, BnRf, is within 400 kb of the BrMs locus. The BnRf orthologs of B. rapa ssp. chinensis were sequenced, and one KASP marker (BrMs_indel) was developed for genotyping based on a 14 bp indel at intron 4. Cosegregation of male sterility and BrMs_indel genotypes in the F2 population indicated that BnRf from B. napus and BrMs from B. rapa are likely to be orthologs. The BrMs_indel marker developed in this study will be useful in marker-assisted selection for the male sterility trait.

A high-density genetic map developed by specific-locus amplified fragment (SLAF) sequencing and identification of a locus controlling anthocyanin pigmentation in stalk of Zicaitai (Brassica rapa L. ssp. chinensis var. purpurea)

BACKGROUND

Caixin and Zicaitai (Brassica rapa) belong to Southern and Central China respectively. Zicaitai contains high amount of anthocyanin in leaf and stalk resulting to the purple color. Stalk is the major edible part and stalk color is an economically important trait for the two vegetables. The aim of this study is to construct a high density genetic map using the specific length amplified fragment sequencing (SLAF-seq) technique to explore genetic basis for anthocyanin pigmentation traits via quantitative trait loci (QTL) mapping.

RESULTS

We constructed a high generation linkage map with a mapping panel of F2 populations derived from 150 individuals of parental lines "Xianghongtai 01" and "Yinong 50D" with purple and green stalk respectively. The map was constructed containing 4253 loci, representing 10,940 single nucleotide polymorphism (SNP) markers spanning 1030.04 centiMorgans (cM) over 10 linkage groups (LGs), with an average distance between markers of 0.27 cM. Quantitative trait loci (QTL) analysis revealed that a major locus on chromosome 7 and 4 minor QTLs explaining 2.69-61.21% of phenotypic variation (PVE) were strongly responsible for variation in stalk color trait. Bioinformatics analysis of the major locus identified 62 protein-coding genes. Among the major locus, there were no biosynthetic genes related to anthocyanin. However, there were several transcription factors like helix-loop-helix (bHLH) bHLH, MYB in the locus. Seven predicted candidate genes were selected for the transcription level analysis. Only bHLH49 transcription factor, was significantly higher expressed in both stalks and young leaves of Xianghongtai01 than Yinong50D. An insertion and deletion (InDel) marker developed from deletion/insertion in the promoter region of bHLH49 showed significant correlation with the stalk color trait in the F2 population.

CONCLUSION

Using the constructed high-qualified linkage map, this study successfully identified QTLs for stalk color trait. The identified valuable markers and candidate genes for anthocyanin accumulation in stalk will provide useful information for molecular regulation of anthocyanin biosynthesis. Overall our findings will lay a foundation for functional gene cloning, marker-assisted selection (MAS) and molecular breeding of important economic traits in B. rapa.

Transcriptomic analysis of contrasting inbred lines and F2 segregant of Chinese cabbage provides valuable information on leaf morphology

BACKGROUND

Leaf morphology influences plant growth and productivity and is controlled by genetic and environmental cues. The various morphotypes of Brassica rapa provide an excellent resource for genetic and molecular studies of morphological traits.

OBJECTIVE

This study aimed to identify genes regulating leaf morphology using segregating B. rapa p F₂ population.

METHODS

Phenotyping and transcriptomic analyses were performed on an F₂ population derived from a cross between Rapid cycling B. rapa (RCBr) and B. rapa ssp. penkinensis, inbred line Kenshin. Analyses focused on four target traits: lamina (leaf) length (LL), lamina width (LW), petiole length (PL), and leaf margin (LM).

RESULTS

All four traits were controlled by multiple QTLs, and expression of 466 and 602 genes showed positive and negative correlation with leaf phenotypes, respectively. From this microarray analysis, large numbers of genes were putatively identified as leaf morphology-related genes. The Gene Ontology (GO) category containing the highest number of differentially expressed genes (DEGs) was "phytohormones". The sets of genes enriched in the four leaf phenotypes did not overlap, indicating that each phenotype was regulated by a different set of genes. The expression of BrAS2, BrAN3, BrCYCB1;2, BrCYCB2;1,4, BrCYCB3;1, CrCYCBD3;2, BrULT1, and BrANT seemed to be related to leaf size traits (LL and LW), whereas BrCUC1, BrCUC2, and BrCUC3 expression for LM trait.

CONCLUSION

An analysis integrating the results of the current study with previously published data revealed that Kenshin alleles largely determined LL and LW but LM resulted from RCBr alleles. Genes identified in this study could be used to develop molecular markers for use in Brassica breeding projects and for the dissection of gene function.

Meta-analysis of chromosome-scale crossover rate variation in eukaryotes and its significance to evolutionary genomics

Understanding the distribution of crossovers along chromosomes is crucial to evolutionary genomics because the crossover rate determines how strongly a genome region is influenced by natural selection on linked sites. Nevertheless, generalities in the chromosome-scale distribution of crossovers have not been investigated formally. We fill this gap by synthesizing joint information on genetic and physical maps across 62 animal, plant and fungal species. Our quantitative analysis reveals a strong and taxonomically widespread reduction of the crossover rate in the centre of chromosomes relative to their peripheries. We demonstrate that this pattern is poorly explained by the position of the centromere, but find that the magnitude of the relative reduction in the crossover rate in chromosome centres increases with chromosome length. That is, long chromosomes often display a dramatically low crossover rate in their centre, whereas short chromosomes exhibit a relatively homogeneous crossover rate. This observation is compatible with a model in which crossover is initiated from the chromosome tips, an idea with preliminary support from mechanistic investigations of meiotic recombination. Consequently, we show that organisms achieve a higher genome-wide crossover rate by evolving smaller chromosomes. Summarizing theory and providing empirical examples, we finally highlight that taxonomically widespread and systematic heterogeneity in crossover rate along chromosomes generates predictable broad-scale trends in genetic diversity and population differentiation by modifying the impact of natural selection among regions within a genome. We conclude by emphasizing that chromosome-scale heterogeneity in crossover rate should urgently be incorporated into analytical tools in evolutionary genomics, and in the interpretation of resulting patterns.