Coupling Seq-BSA and RNA-Seq Analyses Reveal the Molecular Pathway and Genes Associated with Heading Type in Chinese Cabbage

Dissection of quantitative trait nucleotides and candidate genes associated with agronomic and yield-related traits under drought stress in rapeseed varieties: integration of genome-wide association study and transcriptomic analysis

Introduction

An important strategy to combat yield loss challenge is the development of varieties with increased tolerance to drought to maintain production. Improvement of crop yield under drought stress is critical to global food security.

Methods

In this study, we performed multiomics analysis in a collection of 119 diverse rapeseed (Brassica napus L.) varieties to dissect the genetic control of agronomic traits in two watering regimes [well-watered (WW) and drought stress (DS)] for 3 years. In the DS treatment, irrigation continued till the 50% pod development stage, whereas in the WW condition, it was performed throughout the whole growing season.

Results

The results of the genome-wide association study (GWAS) using 52,157 single-nucleotide polymorphisms (SNPs) revealed 1,281 SNPs associated with traits. Six stable SNPs showed sequence variation for flowering time between the two irrigation conditions across years. Three novel SNPs on chromosome C04 for plant weight were located within drought tolerance-related gene ABCG16, and their pleiotropically effects on seed weight per plant and seed yield were characterized. We identified the C02 peak as a novel signal for flowering time, harboring 52.77% of the associated SNPs. The 288-kbps LD decay distance analysis revealed 2,232 candidate genes (CGs) associated with traits. The CGs BIG1-D, CAND1, DRG3, PUP10, and PUP21 were involved in phytohormone signaling and pollen development with significant effects on seed number, seed weight, and grain yield in drought conditions. By integrating GWAS and RNA-seq, 215 promising CGs were associated with developmental process, reproductive processes, cell wall organization, and response to stress. GWAS and differentially expressed genes (DEGs) of leaf and seed in the yield contrasting accessions identified BIG1-D, CAND1, and DRG3 genes for yield variation.

Discussion

The results of our study provide insights into the genetic control of drought tolerance and the improvement of marker-assisted selection (MAS) for breeding high-yield and drought-tolerant varieties.

Combined BSA-Seq and RNA-Seq Reveal Genes Associated with the Visual Stay-Green of Maize (Zea mays L.)

Maize has become one of the most widely grown grains in the world, and the stay-green mutant allows these plants to maintain their green leaves and photosynthetic potential for longer following anthesis than in non-mutated plants. As a result, stay-green plants have a higher production rate than non-stay-green varieties due to their prolonged grain-filling period. In this study, the candidate genes related to the visual stay-green at the maturation stage of maize were investigated. The F2 population was derived from the T01 (stay-green) and the Xin3 (non-stay-green) cross. Two bulked segregant analysis pools were constructed. According to the method of combining ED (Euclidean distance), Ridit (relative to an identified distribution unit), SmoothG, and SNP algorithms, a region containing 778 genes on chromosome 9 was recognized as the candidate region associated with the visual stay-green in maize. A total of eight modules were identified using WGCNA (weighted correlation network analysis), of which green, brown, pink, and salmon modules were significantly correlated with visual stay-green. BSA, combined with the annotation function, discovered 7 potential candidate genes, while WGCNA discovered 11 stay-green potential candidate genes. The candidate range was further reduced due through association analysis of BSA-seq and RNA-seq. We identified Zm00001eb378880, Zm00001eb383680, and Zm00001eb384100 to be the most likely candidate genes. Our results provide valuable insights into this new germplasm resource with reference to increasing the yield for maize.

Role of Brassica orphan gene BrLFM on leafy head formation in Chinese cabbage (Brassica rapa)

Brassica orphan gene BrFLM, identified by two allelic mutants, was involved in leafy head formation in Chinese cabbage. Leafy head formation is a unique agronomic trait of Chinese cabbage that determines its yield and quality. In our previous study, an EMS mutagenesis Chinese cabbage mutant library was constructed using the heading Chinese cabbage double haploid (DH) line FT as the wild-type. Here, we screened two extremely similar leafy head deficiency mutants lfm-1 and lfm-2 with geotropic growth leaves from the library to investigate the gene(s) related to leafy head formation. Reciprocal crossing results showed that these two mutants were allelic. We utilized lfm-1 to identify the mutant gene(s). Genetic analysis showed that the mutated trait was controlled by a single nuclear gene Brlfm. Mutmap analysis showed that Brlfm was located on chromosome A05, and BraA05g012440.3C or BraA05g021450.3C were the candidate gene. Kompetitive allele-specific PCR analysis eliminated BraA05g012440.3C from the candidates. Sanger sequencing identified an SNP from G to A at the 271st nucleotide on BraA05g021450.3C. The sequencing of lfm-2 detected another non-synonymous SNP (G to A) located at the 266st nucleotide on BraA05g021450.3C, which verified its function on leafy head formation. We blasted BraA05g021450.3C on database and found that it belongs to a Brassica orphan gene encoding an unknown 13.74 kDa protein, named BrLFM. Subcellular localization showed that BrLFM was located in the nucleus. These findings reveal that BrLFM is involved in leafy head formation in Chinese cabbage.

Auxin participates in regulating the leaf curl development of Wucai (Brassica campestris L.)

Wucai (Brassica campestris L. ssp. chinensis var. rosularis Tsen) belongs to the Brassica genus of the Cruciferae family, and its leaf curl is a typical feature that distinguishes Wucai from other nonheading cabbage subspecies. Our previous research found that plant hormones were involved in the development of the leaf curl in Wucai. However, the molecular mechanisms and the hormones regulating the formation of leaf curl in Wucai have not yet been reported. This study aimed to understand the molecular functions related to hormone metabolism during the formation of leaf curl in Wucai. A total of 386 differentially expressed genes (DEGs) were identified by transcriptome sequencing of two different morphological parts of the same leaf of Wucai germplasm W7-2, and 50 DEGs were found to be related to plant hormones, which were mainly involved in the auxin signal transduction pathway. Then, we measured the content of endogenous hormones in two different forms of the same leaf of Wucai germplasm W7-2. A total of 17 hormones with differential content were identified, including auxin, cytokinins, jasmonic acids, salicylic acids, and abscisic acid. And we found that treatment with auxin transport inhibitor N-1-naphthylphthalamic acid can affect the leaf curl phenotype of Wucai and pak choi (Brassica rapa L. subsp. Chinensis). These results indicated that plant hormones, especially auxin, are involved in developing the leaf curl of Wucai. Our findings provide a potentially valuable reference for future research on the development of leaf curls.

Identification and characterization analysis of candidate genes controlling mushroom leaf development in Chinese kale by BSA-seq

Mushroom leaves (MLs) are malformed leaves that develop from the leaf veins in some of Chinese kale genotypes. To study the genetic model and molecular mechanism of ML development in Chinese kale, the F2 segregation population was constructed by two inbred lines, genotype Boc52 with ML and genotype Boc55 with normal leaves (NL). In the present study, we have identified for the first time that the development of mushroom leaves may be affected by the change of adaxial-abaxial polarity of leaves. Examination of the phenotypes of F1 and F2 segregation populations suggested that ML development is controlled by two dominant major genes inherited independently. BSA-seq analysis showed that a major quantitative trait locus (QTL) qML4.1 that controls ML development is located within 7.4 Mb on chromosome kC4. The candidate region was further narrowed to 255 kb by linkage analysis combined with insertion/deletion (InDel) markers, and 37 genes were predicted in this region. According to the expression and annotation analysis, a B3 domain-containing transcription factor NGA1-like gene, BocNGA1, was identified as a key candidate gene for controlling ML development in Chinese kale. Fifteen single nucleotide polymorphisms (SNPs) were found in coding sequences and 21 SNPs and 3 InDels found in the promoter sequences of BocNGA1 from the genotype Boc52 with ML. The expression levels of BocNGA1 in ML genotypes are significantly lower than in the NL genotypes, which suggests that BocNGA1 may act as a negative regulator for ML genesis in Chinese kale. This study provides a new foundation for Chinese kale breeding and for the study of the molecular mechanism of plant leaf differentiation.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01364-6.

Targeted Mutagenesis of the Multicopy Myrosinase Gene Family in Allotetraploid Brassica juncea Reduces Pungency in Fresh Leaves across Environments

Recent breeding efforts in Brassica have focused on the development of new oilseed feedstock crop for biofuels (e.g., ethanol, biodiesel, bio-jet fuel), bio-industrial uses (e.g., bio-plastics, lubricants), specialty fatty acids (e.g., erucic acid), and producing low glucosinolates levels for oilseed and feed meal production for animal consumption. We identified a novel opportunity to enhance the availability of nutritious, fresh leafy greens for human consumption. Here, we demonstrated the efficacy of disarming the ‘mustard bomb’ reaction in reducing pungency upon the mastication of fresh tissue—a major source of unpleasant flavor and/or odor in leafy Brassica. Using gene-specific mutagenesis via CRISPR-Cas12a, we created knockouts of all functional copies of the type-I myrosinase multigene family in tetraploid Brassica juncea. Our greenhouse and field trials demonstrate, via sensory and biochemical analyses, a stable reduction in pungency in edited plants across multiple environments. Collectively, these efforts provide a compelling path toward boosting the human consumption of nutrient-dense, fresh, leafy green vegetables.

OCTOPUS regulates BIN2 to control leaf curvature in Chinese cabbage

Heading is one of the most important agronomic traits for Chinese cabbage crops. During the heading stage, leaf axial growth is an essential process. In the past, most genes predicted to be involved in the heading process have been based on leaf development studies in Arabidopsis. No genes that control leaf axial growth have been mapped and cloned via forward genetics in Chinese cabbage. In this study, we characterize the inward curling mutant ic1 in Brassica rapa ssp. pekinensis and identify a mutation in the OCTOPUS (BrOPS) gene by map-based cloning. OPS is involved in phloem differentiation in Arabidopsis, a functionalization of regulating leaf curvature that is differentiated in Chinese cabbage. In the presence of brassinosteroid (BR) at the early heading stage in ic1, the mutation of BrOPS fails to sequester brassinosteroid insensitive 2 (BrBIN2) from the nucleus, allowing BrBIN2 to phosphorylate and inactivate BrBES1, which in turn relieves the repression of BrAS1 and results in leaf inward curving. Taken together, the results of our findings indicate that BrOPS positively regulates BR signaling by antagonizing BrBIN2 to promote leaf epinastic growth at the early heading stage in Chinese cabbage.

BrAN contributes to leafy head formation by regulating leaf width in Chinese cabbage (Brassica rapa L. ssp. pekinensis)

Leafy head is an important agronomic trait that determines the yield and quality of Chinese cabbage. The molecular mechanism underlying heading in Chinese cabbage has been the focus of research, and wide leaves are a prerequisite for leafy head formation. In our study, two allelic leafy heading-deficient mutants (lhd1 and lhd2) with narrow leaf phenotypes were screened in an ethyl methanesulfonate mutagenized population from a heading Chinese cabbage double haploid line 'FT'. Genetic analysis revealed that the mutant trait was controlled by a recessive nuclear gene, which was found to be BraA10g000480.3C by MutMap and Kompetitive allele-specific PCR analyses. As BraA10g000480.3C was the ortholog of ANGUSTIFOLIA in Arabidopsis, which has been found to regulate leaf width by controlling cortical microtubule arrangement and pavement cell shape, we named it BrAN. BrAN in mutant lhd1 carried an SNP (G to A) on intron 2 that co-segregated with the mutant phenotype, and disrupted the exon-intron splice junction generating intron retention and a putative truncated protein. BrAN in mutant lhd2 carried an SNP (G to A) on exon 4 leading to a premature stop codon. The ectopic overexpression of BrAN restored normal leaf phenotype due to abnormal cortical microtubule arrangement and pavement cell shape in the Arabidopsis an-t1 mutant. However, transformation of Bran did not rescue the an-t1 phenotype. These results indicate that BrAN contributes to leafy head formation of Chinese cabbage.

Combined Analysis of BSA-Seq Based Mapping, RNA-Seq, and Metabolomic Unraveled Candidate Genes Associated with Panicle Grain Number in Rice (Oryza sativa L.)

Rice grain yield is a complex and highly variable quantitative trait consisting of several key components, including the grain weight, the effective panicles per unit area, and the grain number per panicle (GNPP). The GNPP is a significant contributor to grain yield controlled by multiple genes (QTL) and is crucial for improvement. Attempts have been made to find genes for this trait, which has always been a challenging and arduous task through conventional methods. We combined a BSA analysis, RNA profiling, and a metabolome analysis in the present study to identify new candidate genes involved in the GNPP. The F₂ population from crossing R4233 (high GNPP) and Ce679 (low GNPP) revealed a frequency distribution fitting two segregated genes. Three pools, including low, middle, and high GNPP, were constructed and a BSA analysis revealed six candidate regions spanning 5.38 Mb, containing 739 annotated genes. Further, a conjunctive analysis of BSA-Seq and RNA-Seq showed 31 differentially expressed genes (DEGs) in the candidate intervals. Subsequently, a metabolome analysis showed 1024 metabolites, with 71 significantly enriched, including 44 up and 27 downregulated in Ce679 vs. R4233. A KEGG enrichment analysis of these 31 DEGs and 71 differentially enriched metabolites (DEMs) showed two genes, Os12g0102100 and Os01g0580500, significantly enriched in the metabolic pathways’ biosynthesis of secondary metabolites, cysteine and methionine metabolism, and fatty acid biosynthesis. Os12g0102100, which encodes for the alcohol dehydrogenase superfamily and a zinc-containing protein, is a novel gene whose contribution to the GNPP is not yet elucidated. This gene coding for mitochondrial trans-2-enoyl-CoA reductase is involved in the biosynthesis of myristic acid, also known as tetradecanoic acid. The Os01g0580500 coding for the enzyme 1-aminoclopropane-1-carboxylate oxidase (OsACO7) is responsible for the final step of the ethylene biosynthesis pathway through the conversion of 1-aminocyclopropane-1-carboxylic acid (ACC) into ethylene. Unlike Os12g0102100, this gene was significantly upregulated in R4233, downregulated in Ce679, and significantly enriched in two of the three metabolite pathways. This result pointed out that these two genes are responsible for the difference in the GNPP in the two cultivars, which has never been identified. Further validation studies may disclose the physiological mechanisms through which they regulate the GNPP in rice.

The Adaxial/Abaxial Patterning of Auxin and Auxin Gene in Leaf Veins Functions in Leafy Head Formation of Chinese Cabbage

Leaf curling is an essential prerequisite for the formation of leafy heads in Chinese cabbage. However, the part or tissue that determines leaf curvature remains largely unclear. In this study, we first introduced the auxin-responsive marker DR5::GUS into the Chinese cabbage genome and visualized its expression during the farming season. We demonstrated that auxin response is adaxially/abaxially distributed in leaf veins. Together with the fact that leaf veins occupy considerable proportions of the Chinese cabbage leaf, we propose that leaf veins play a crucial supporting role as a framework for heading. Then, by combining analyses of QTL mapping and a time-course transcriptome from heading Chinese cabbage and non-heading pak choi during the farming season, we identified the auxin-related gene BrPIN5 as a strong candidate for leafy head formation. PIN5 displays an adaxial/abaxial expression pattern in leaf veins, similar to that of DR5::GUS, revealing an involvement of BrPIN5 in leafy head development. The association of BrPIN5 function with heading was further confirmed by its haplo-specificity to heading individuals in both a natural population and two segregating populations. We thus conclude that the adaxial/abaxial patterning of auxin and auxin genes in leaf veins functions in the formation of the leafy head in Chinese cabbage.

Detection of a major QTL and development of KASP markers for seed weight by combining QTL-seq, QTL-mapping and RNA-seq in peanut

Combining QTL-seq, QTL-mapping and RNA-seq identified a major QTL and candidate genes, which contributed to the development of KASP markers and understanding of molecular mechanisms associated with seed weight in peanut. Seed weight, as an important component of seed yield, is a significant target of peanut breeding. However, relatively little is known about the quantitative trait loci (QTLs) and candidate genes associated with seed weight in peanut. In this study, three major QTLs on chromosomes A05, B02, and B06 were determined by applying the QTL-seq approach in a recombinant inbred line (RIL) population. Based on conventional QTL-mapping, these three QTL regions were successfully narrowed down through newly developed single nucleotide polymorphism (SNP) and simple sequence repeat markers. Among these three QTL regions, qSWB06.3 exhibited stable expression, contributing mainly to phenotypic variance across environments. Furthermore, differentially expressed genes (DEGs) were identified at the three seed developmental stages between the two parents of the RIL population. It was found that the DEGs were widely distributed in the ubiquitin-proteasome pathway, the serine/threonine-protein pathway, signal transduction of hormones and transcription factors. Notably, DEGs at the early stage were mostly involved in regulating cell division, whereas DEGs at the middle and late stages were primarily involved in cell expansion during seed development. The expression patterns of candidate genes related to seed weight in qSWB06.3 were investigated using quantitative real-time PCR. In addition, the allelic diversity of qSWB06.3 was investigated in peanut germplasm accessions. The marker Ah011475 has higher efficiency for discriminating accessions with different seed weights, and it would be useful as a diagnostic marker in marker-assisted breeding. This study provided insights into the genetic and molecular mechanisms of seed weight in peanut.

A major yellow-seed QTL on chromosome A09 significantly increases the oil content and reduces the fiber content of seed in Brassica napus

A major yellow-seed QTL on chromosome A09 significantly increases the oil content and reduces the fiber content of seed in Brassica napus. The yellow-seed trait (YST) has always been a main breeding objective for rapeseed because yellow-seeded B. napus generally contains higher oil contents, fewer pigments and polyphenols and lower fiber content than black-seeded B. napus, although the mechanism controlling this correlation remains unclear. In this study, QTL mapping was implemented for YST based on a KN double haploid population derived from the hybridization of yellow-seeded B. napus N53-2 with a high oil content and black-seeded Ken-C8 with a relatively low oil content. Ten QTLs were identified, including four stable QTLs that could be detected in multiple environments. A major QTL, cqSC-A09, on chromosome A09 was identified by both QTL mapping and BSR-Seq technology, and explained more than 41% of the phenotypic variance. The major QTL cqSC-A09 for YST not only controls the seed color but also affects the oil and fiber contents in seeds. More importantly, the advantageous allele could increase the oil content and reduce the pigment and fiber content at the same time. This is the first QTL reported to control seed color, oil content and fiber content simultaneously with a large effect and has great application value for breeding high oil varieties with high seed quality. Important candidate genes, including BnaA09. JAZ1, BnaA09. GH3.3 and BnaA09. LOX3, were identified for cqSC-A09 by combining sequence variation annotation, expression differences and an interaction network, which lays a foundation for further cloning and breeding applications in the future.

Series-Spatial Transcriptome Profiling of Leafy Head Reveals the Key Transition Leaves for Head Formation in Chinese Cabbage

Chinese cabbage is an important leaf heading vegetable crop. At the heading stage, its leaves across inner to outer show significant morphological differentiation. However, the genetic control of this complex leaf morphological differentiation remains unclear. Here, we reported the transcriptome profiling of Chinese cabbage plant at the heading stage using 24 spatially dissected tissues representing different regions of the inner to outer leaves. Genome-wide transcriptome analysis clearly separated the inner leaf tissues from the outer leaf tissues. In particular, we identified the key transition leaf by the spatial expression analysis of key genes for leaf development and sugar metabolism. We observed that the key transition leaves were the first inwardly curved ones. Surprisingly, most of the heading candidate genes identified by domestication selection analysis obviously showed a corresponding expression transition, supporting that key transition leaves are related to leafy head formation. The key transition leaves were controlled by a complex signal network, including not only internal hormones and protein kinases but also external light and other stimuli. Our findings provide new insights and the rich resource to unravel the genetic control of heading traits.

Genetic analysis of the "head top shape" quality trait of Chinese cabbage and its association with rosette leaf variation

The agricultural and consumer quality of Chinese cabbage is determined by its shape. The shape is defined by the folding of the heading leaves, which defines the head top shape (HTS). The overlapping HTS, in which the heading leaves curve inward and overlap at the top, is the shape preferred by consumers. To understand the genetic regulation of HTS, we generated a large segregating F₂ population from a cross between pak choi and Chinese cabbage, with phenotypes ranging from nonheading to heading with either outward curving or inward curving overlapping heading leaves. HTS was correlated with plant height, outer/rosette leaf length, and petiole length. A high-density genetic map was constructed. Quantitative trait locus (QTL) analysis resulted in the identification of 22 QTLs for leafy head-related traits, which included five HTS QTLs. Bulked segregant analysis (BSA) was used to confirm HTS QTLs and identify candidate genes based on informative single-nucleotide polymorphisms. Interestingly, the HTS QTLs colocalized with QTLs for plant height, outer/rosette leaf, and petiole length, consistent with the observed phenotypic correlations. Combined QTL analysis and BSA laid a foundation for molecular marker-assisted breeding of Chinese cabbage HTS and directions for further research on the genetic regulation of this trait.

Physical mapping and InDel marker development for the restorer gene Rf2 in cytoplasmic male sterile CMS-D8 cotton

Background

Cytoplasmic male sterile (CMS) with cytoplasm from Gossypium Trilobum (D8) fails to produce functional pollen. It is useful for commercial hybrid cotton seed production. The restore line of CMS-D8 containing Rf₂ gene can restore the fertility of the corresponding sterile line. This study combined the whole genome resequencing bulked segregant analysis (BSA) with high-throughput SNP genotyping to accelerate the physical mapping of Rf₂ locus in CMS-D8 cotton.

Methods

The fertility of backcross population ((sterile line×restorer line)×maintainer line) comprising of 1623 individuals was investigated in the field. The fertile pool (100 plants with fertile phenotypes, F-pool) and the sterile pool (100 plants with sterile phenotypes, S-pool) were constructed for BSA resequencing. The selection of 24 single nucleotide polymorphisms (SNP) through high-throughput genotyping and the development insertion and deletion (InDel) markers were conducted to narrow down the candidate interval. The pentapeptide repeat (PPR) family genes and upregulated genes in restore line in the candidate interval were analysed by qRT-PCR.

Results

The fertility investigation results showed that fertile and sterile separation ratio was consistent with 1:1. BSA resequencing technology, high-throughput SNP genotyping, and InDel markers were used to identify Rf₂ locus on candidate interval of 1.48 Mb on chromosome D05. Furthermore, it was quantified in this experiment that InDel markers co-segregated with Rf₂ enhanced the selection of the restorer line. The qRT-PCR analysis revealed PPR family gene Gh_D05G3391 located in candidate interval had significantly lower expression than sterile and maintainer lines. In addition, utilization of anther RNA-Seq data of CMS-D8 identified that the expression level of Gh_D05G3374 encoding NB-ARC domain-containing disease resistance protein in restorer lines was significantly higher than that in sterile and maintainer lines.

Conclusions

This study not only enabled us to precisely locate the restore gene Rf₂ but also evaluated the utilization of InDel markers for marker assisted selection in the CMS-D8 Rf₂ cotton breeding line. The results of this study provide an important foundation for further studies on the mapping and cloning of restorer genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07342-y.

The mutation of ent-kaurene synthase, a key enzyme involved in gibberellin biosynthesis, confers a non-heading phenotype to Chinese cabbage (Brassica rapa L. ssp. pekinensis)

The presence of a leafy head is a vital agronomic trait that facilitates the evaluation of the yield and quality of Chinese cabbage. A non-heading mutant (nhm1) was identified in an ethyl methanesulfonate mutagenesis population of the heading Chinese cabbage double haploid line FT. Segregation analysis revealed that a single recessive gene, Brnhm1, controlled the mutant phenotype. Using MutMap, Kompetitive allele-specific PCR, and cloning analyses, we demonstrated that BraA07g042410.3C, which encodes an ent-kaurene synthase involved in the gibberellin biosynthesis pathway, is the nhm1 mutant candidate gene. A single-nucleotide mutation (C to T) in the fourth exon of BraA07g042410.3C caused an amino acid substitution from histidine to tyrosine. Compared to that of the wild-type FT, BraA07g042410.3C in the leaves of the nhm1 mutant had lower levels of expression. In addition, gibberellin contents were lower in the mutant than in the wild type, and the mutant plant phenotype could be restored to that of the wild type after exogenous GA₃ treatment. These results indicate that BraA07g042410.3C caused the non-heading mutation. This is the first study to demonstrate a relationship between gibberellin content in the leaves and leafy head formation in Chinese cabbage. These findings facilitate the understanding of the mechanisms underlying leafy head development in Chinese cabbage.

Fine Mapping to Identify the Functional Genetic Locus for Red Coloration in Pyropia yezoensis Thallus

Pyropia yezoensis, commonly known as "Nori" or "Laver" is an economically important marine crop. In natural or selected populations of P. yezoensis, coloration mutants are frequently observed. Various coloration mutants are excellent materials for genetic research and study photosynthesis. However, the candidate gene controlling the Pyropia coloration phenotype remains unclear to date. QTL-seq, in combination with kompetitive allele-specific PCR (KASP) and RNA-seq, can be generally applied to population genomics studies to rapidly identify genes that are responsible for phenotypes showing extremely opposite traits. Through cross experiments between the wild line RZ and red-mutant HT, offsprings with 1-4 sectors chimeric blade were generated. Statistical analyses revealed that the red thallus coloration phenotype is conferred by a single nuclear allele. Two-pair populations, consisting of 24 and 56 wild-type/red-type single-genotype sectors from F1 progeny, were used in QTL-seq to detect a genomic region in P. yezoensis harboring the red coloration locus. Based on a high-quality genome, we first identified the candidate region within a 3.30-Mb region at the end of chromosome 1. Linkage map-based QTL analysis was used to confirm the candidate region identified by QTL-seq. Then, four KASP markers developed in this region were used to narrow down the candidate region to a 1.42-Mb region. Finally, we conducted RNA-seq to focus on 13 differentially expressed genes and further predicted rcl-1, which contains one non-synonymous SNP [A/C] in the coding region that could be regulating red thallus coloration in P. yezoensis. Our results provide novel insights into the underlying mechanism controlling blade coloration, which is a desirable trait in algae.

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.

Using a transcriptome sequencing approach to explore candidate resistance genes against stemphylium blight in the wild lentil species Lens ervoides

BACKGROUND

Stemphylium blight (SB), caused by Stemphylium botryosum, is a devastating disease in lentil production. Although it is known that accessions of Lens ervoides possess superior SB resistance at much higher frequency than the cultivated lentil species, very little is known about the molecular basis regulating SB resistance in L. ervoides. Therefore, a comprehensive molecular study of SB resistance in L. ervoides was needed to exploit this wild resource available at genebanks for use by plant breeders in resistance breeding.

RESULTS

Microscopic and qPCR quantification of fungal growth revealed that 48, 96, and 144 h post-inoculation (hpi) were interesting time points for disease development in L. ervoides recombinant inbred lines (RILs) LR-66-637 (resistant to SB) and LR-66-577 (susceptible to SB). Results of transcriptome sequencing at 0, 48, 96 and 144 hpi showed that 8810 genes were disease-responsive genes after challenge by S. botryosum. Among them, 7526 genes displayed a similar expression trend in both RILs, and some of them were likely involved in non-host resistance. The remaining 1284 genes were differentially expressed genes (DEGs) between RILs. Of those, 712 DEGs upregulated in LR-66-637 were mostly enriched in 'carbohydrate metabolic process', 'cell wall organization or biogenesis', and 'polysaccharide metabolic process'. In contrast, there were another 572 DEGs that were upregulated in LR-66-577, and some of them were enriched in 'oxidation-reduction process', 'asparagine metabolic process' and 'asparagine biosynthetic process'. After comparing DEGs to genes identified in previously described quantitative trait loci (QTLs) for resistance to SB, nine genes were common and three of them showed differential gene expression between a resistant and a susceptible bulk consisting of five RILs each. Results showed that two genes encoding calcium-transporting ATPase and glutamate receptor3.2 were candidate resistance genes, whereas one gene with unknown function was a candidate susceptibility gene.

CONCLUSION

This study provides new insights into the mechanisms of resistance and susceptibility in L. ervoides RILs responding to S. botryosum infection. Furthermore, we identified candidate resistance or susceptibility genes which warrant further gene function analyses, and which could be valuable for resistance breeding, if their role in resistance or susceptibility can be confirmed.

Transcription Coactivator ANGUSTIFOLIA3 (AN3) Regulates Leafy Head Formation in Chinese Cabbage

Leafy head formation in Chinese cabbage (B. rapa ssp. pekinensis cv. Bre) results from leaf curvature, which is under the tight control of genes involved in the adaxial-abaxial patterning during leaf development. The transcriptional coactivator ANGUSTIFOLIA3 (AN3) binds to the SWI/SNF chromatin remodeling complexes formed around ATPases such as BRAHMA (BRM) in order to regulate transcription in various aspects of leaf development such as cell proliferation, leaf primordia expansion, and leaf adaxial/abaxial patterning in Arabidopsis. However, its regulatory function in Chinese cabbage remains poorly understood. Here, we analyzed the expression patterns of the Chinese cabbage AN3 gene (BrAN3) before and after leafy head formation, and produced BrAN3 gene silencing plants by using the turnip yellow mosaic virus (TYMV)-derived vector in order to explore its potential function in leafy head formation in Chinese cabbage. We found that BrAN3 had distinct expression patterns in the leaves of Chinese cabbage at the rosette and heading stages. We also found silencing of BrAN3 stimulated leafy head formation at the early stage. Transcriptome analysis indicated that silencing of BrAN3 modulated the hormone signaling pathways of auxin, ethylene, GA, JA, ABA, BR, CK, and SA in Chinese cabbage. Our study offers unique insights into the function of BrAN3 in leafy head formation in Chinese cabbage.

Characterization of Non-heading Mutation in Heading Chinese Cabbage (Brassica rapa L. ssp. pekinensis)

Heading is a key agronomic trait of Chinese cabbage. A non-heading mutant with flat growth of heading leaves (fg-1) was isolated from an EMS-induced mutant population of the heading Chinese cabbage inbred line A03. In fg-1 mutant plants, the heading leaves are flat similar to rosette leaves. The epidermal cells on the adaxial surface of these leaves are significantly smaller, while those on the abaxial surface are much larger than in A03 plants. The segregation of the heading phenotype in the F2 and BC1 population suggests that the mutant trait is controlled by a pair of recessive alleles. Phytohormone analysis at the early heading stage showed significant decreases in IAA, ABA, JA and SA, with increases in methyl IAA and trans-Zeatin levels, suggesting they may coordinate leaf adaxial-abaxial polarity, development and morphology in fg-1. RNA-sequencing analysis at the early heading stage showed a decrease in expression levels of several auxin transport (BrAUX1, BrLAXs, and BrPINs) and responsive genes. Transcript levels of important ABA responsive genes, including BrABF3, were up-regulated in mid-leaf sections suggesting that both auxin and ABA signaling pathways play important roles in regulating leaf heading. In addition, a significant reduction in BrIAMT1 transcripts in fg-1 might contribute to leaf epinastic growth. The expression profiles of 19 genes with known roles in leaf polarity were significantly different in fg-1 leaves compared to wild type, suggesting that these genes might also regulate leaf heading in Chinese cabbage. In conclusion, leaf heading in Chinese cabbage is controlled through a complex network of hormone signaling and abaxial-adaxial patterning pathways. These findings increase our understanding of the molecular basis of head formation in Chinese cabbage.

The Gene Structure and Expression Level Changes of the GH3 Gene Family in Brassica napus Relative to Its Diploid Ancestors

The GH3 gene family plays a vital role in the phytohormone-related growth and developmental processes. The effects of allopolyploidization on GH3 gene structures and expression levels have not been reported. In this study, a total of 38, 25, and 66 GH3 genes were identified in Brassica rapa (A_rA_r), Brassica oleracea (C_oC_o), and Brassica napus (A_nAC_nC_n), respectively. BnaGH3 genes were unevenly distributed on chromosomes with 39 on A_n and 27 on C_n, in which six BnaGH3 genes may appear as new genes. The whole genome triplication allowed the GH3 gene family to expand in diploid ancestors, and allopolyploidization made the GH3 gene family re-expand in B. napus. For most BnaGH3 genes, the exon-intron compositions were similar to diploid ancestors, while the cis-element distributions were obviously different from its ancestors. After allopolyploidization, the expression patterns of GH3 genes from ancestor species changed greatly in B. napus, and the orthologous gene pairs between A_n/A_r and C_n/C_o had diverged expression patterns across four tissues. Our study provides a comprehensive analysis of the GH3 gene family in B. napus, and these results could contribute to identifying genes with vital roles in phytohormone-related growth and developmental processes.

Molecular mapping of the Cf-10 gene by combining SNP/InDel-index and linkage analysis in tomato (Solanum lycopersicum)

BACKGROUND

Leaf mold, one of the major diseases of tomato caused by Cladosporium fulvum (C. fulvum), can dramatically reduce the yield and cause multimillion dollar losses annually worldwide. Mapping the resistance genes (R genes) of C. fulvum and devising MAS based strategies for breeding new cultivars is an effective approach to improve the resistance in tomato. Up to now, many C. fulvum genes or QTLs have been mapped using different genetic materials, but few studies focused on Cf-10 gene positioning.

RESULTS

In this study, we investigated the genetic rules for Cf-10 and used a novel combinatorial strategy to rapidly map the Cf-10 gene. Initially, the performance of F1, F2 and BC1F1 individuals after infection, demonstrated that the resistance against C. fulvum was controlled by a single dominant gene. Two pools of resistant and susceptible individuals from F2 population were investigated, using mapping by sequencing approach and Cf-10 was found to be localized to 3.35 Mb and 3.74 Mb on chromosome 1, employing SNP/InDel index methods, respectively. After accounting for overlapping regions, these two algorithms yielded a total length of 3.29 Mb, narrowing down the target region. We further developed five serviceable KASP markers for this region based on sequencing data and conducted local QTL mapping using individuals from the F2 population, except for mapping by sequencing as mentioned above. Finally Cf-10 gene was mapped spanning a region of 790 kb, where only one gene (Solyc01g007130.3) was annotated as probable receptor protein kinase TMK1 with a LRR motif, a common R gene characteristic. The RT-qPCR analysis further confirmed the localization and the relative expression of Solyc01g007130.3 in Ontario 792 and was found to be significantly higher than that in Moneymaker at 9 dpi and 12 dpi, respectively.

CONCLUSION

This study proposed a novel combinatorial strategy by combining SNP-index, InDel-index analyses and local QTL mapping using KASP genotyping approach to rapidly map genes responsible for specific traits and provided a robust base for cloning the Cf-10 gene. Furthermore, these analyses suggest that Solyc01g007130.3 is a potential candidate to be regarded as Cf-10 gene.

Mining of Brassica-Specific Genes (BSGs) and Their Induction in Different Developmental Stages and under Plasmodiophora brassicae Stress in Brassica rapa

Orphan genes, also called lineage-specific genes (LSGs), are important for responses to biotic and abiotic stresses, and are associated with lineage-specific structures and biological functions. To date, there have been no studies investigating gene number, gene features, or gene expression patterns of orphan genes in Brassica rapa. In this study, 1540 Brassica-specific genes (BSGs) and 1824 Cruciferae-specific genes (CSGs) were identified based on the genome of Brassica rapa. The genic features analysis indicated that BSGs and CSGs possessed a lower percentage of multi-exon genes, higher GC content, and shorter gene length than evolutionary-conserved genes (ECGs). In addition, five types of BSGs were obtained and 145 out of 529 real A subgenome-specific BSGs were verified by PCR in 51 species. In silico and semi-qPCR, gene expression analysis of BSGs suggested that BSGs are expressed in various tissue and can be induced by Plasmodiophora brassicae. Moreover, an A/C subgenome-specific BSG, BSGs1, was specifically expressed during the heading stage, indicating that the gene might be associated with leafy head formation. Our results provide valuable biological information for studying the molecular function of BSGs for Brassica-specific phenotypes and biotic stress in B. rapa.