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

The Roles of Gibberellins in Regulating Leaf Development

Plant growth and development are correlated with many aspects, including phytohormones, which have specific functions. However, the mechanism underlying the process has not been well elucidated. Gibberellins (GAs) play fundamental roles in almost every aspect of plant growth and development, including cell elongation, leaf expansion, leaf senescence, seed germination, and leafy head formation. The central genes involved in GA biosynthesis include GA20 oxidase genes (GA20oxs), GA3oxs, and GA2oxs, which correlate with bioactive GAs. The GA content and GA biosynthesis genes are affected by light, carbon availability, stresses, phytohormone crosstalk, and transcription factors (TFs) as well. However, GA is the main hormone associated with BR, ABA, SA, JA, cytokinin, and auxin, regulating a wide range of growth and developmental processes. DELLA proteins act as plant growth suppressors by inhibiting the elongation and proliferation of cells. GAs induce DELLA repressor protein degradation during the GA biosynthesis process to control several critical developmental processes by interacting with F-box, PIFS, ROS, SCLl3, and other proteins. Bioactive GA levels are inversely related to DELLA proteins, and a lack of DELLA function consequently activates GA responses. In this review, we summarized the diverse roles of GAs in plant development stages, with a focus on GA biosynthesis and signal transduction, to develop new insight and an understanding of the mechanisms underlying plant development.

Screening of Candidate Genes Associated with Brown Stripe Resistance in Sugarcane via BSR-seq Analysis

Sugarcane brown stripe (SBS), caused by the fungal pathogen Helminthosporium stenospilum, is one of the most serious threats to sugarcane production. However, its outbreaks and epidemics require suitable climatic conditions, resulting in the inefficient improvement of the SBS resistance by phenotype selection. The sugarcane F₁ population of SBS-resistant YT93-159 × SBS-susceptible ROC22 was used for constructing the bulks. Bulked segregant RNA-seq (BSR-seq) was then performed on the parents YT93-159 (T01) and ROC22 (T02), and the opposite bulks of 30 SBS-susceptible individuals mixed bulk (T03) and 30 SBS-resistant individuals mixed bulk (T04) collected from 287 F₁ individuals. A total of 170.00 Gb of clean data containing 297,921 SNPs and 70,426 genes were obtained. Differentially expressed genes (DEGs) analysis suggested that 7787 and 5911 DEGs were identified in the parents (T01 vs. T02) and two mixed bulks (T03 vs. T04), respectively. In addition, 25,363 high-quality and credible SNPs were obtained using the genome analysis toolkit GATK for SNP calling. Subsequently, six candidate regions with a total length of 8.72 Mb, which were located in the chromosomes 4B and 7C of sugarcane wild species Saccharum spontaneum, were identified, and 279 genes associated with SBS-resistance were annotated by ED algorithm and ΔSNP-index. Furthermore, the expression profiles of candidate genes were verified by quantitative real-time PCR (qRT-PCR) analysis, and the results showed that eight genes (LRR-RLK, DHAR1, WRKY7, RLK1, BLH4, AK3, CRK34, and NDA2) and seven genes (WRKY31, CIPK2, CKA1, CDPK6, PFK4, CBL2, and PR2) of the 20 tested genes were significantly up-regulated in YT93-159 and ROC22, respectively. Finally, a potential molecular mechanism of sugarcane response to H. stenospilum infection is illustrate that the activations of ROS signaling, MAPK cascade signaling, Ca²⁺ signaling, ABA signaling, and the ASA-GSH cycle jointly promote the SBS resistance in sugarcane. This study provides abundant gene resources for the SBS resistance breeding in sugarcane.

Phenotypic analysis of Longya-10 × pale flax hybrid progeny and identification of candidate genes regulating prostrate/erect growth in flax plants

Flax is a dual-purpose crop that is important for oil and fiber production. The growth habit is one of the crucial targets of selection during flax domestication. Wild hybridization between cultivated flax and wild flax can produce superior germplasms for flax breeding and facilitate the study of the genetic mechanism underlying agronomically important traits. In this study, we used pale flax, Linum grandiflorum, and L. perenne to pollinate Longya-10. Only pale flax interspecific hybrids were obtained, and the trait analysis of the F₁ and F₂ generations showed that the traits analyzed in this study exhibited disparate genetic characteristics. In the F₁ generation, only one trait, i.e., the number of capsules per plant (140) showed significant heterosis, while the characteristics of other traits were closely associated with those of the parents or a decline in hybrid phenotypes. The traits of the F₂ generation were widely separated, and the variation coefficient ranged from 9.96% to 146.15%. The quantitative trait locus underlying growth habit was preliminarily found to be situated on chromosome 2 through Bulked-segregant analysis sequencing. Then linkage mapping analysis was performed to fine-map GH2.1 to a 23.5-kb interval containing 4 genes. Among them, L.us.o.m.scaffold22.109 and L.us.o.m.scaffold22.112 contained nonsynonymous SNPs with Δindex=1. Combined with the qRT-PCR results, the two genes might be possible candidate genes for GH2.1. This study will contribute to the development of important germplasms for flax breeding, which would facilitate the elucidation of the genetic mechanisms regulating the growth habit and development of an ideal architecture for the flax plant.

Conuping BSA-Seq and RNA-Seq Reveal the Molecular Pathway and Genes Associated with the Plant Height of Foxtail Millet (Setaria italica)

Foxtail millet (Setaria italica) plays an important role in C4 crop research and agricultural development in arid areas due to its short growth period, drought tolerance, and barren tolerance. Exploration of the dwarfing mechanism and the dwarf genes of foxtail millet can provide a reference for dwarf breeding and dwarf research of other C4 crops. In this study, genetic analysis was performed using phenotypic data, candidate genes were screened by bulk segregant analysis sequencing (BSA-Seq); differentially expressed genes and metabolic pathways in different strains of high samples were analyzed by RNA sequencing (RNA-Seq). The association analysis of BSA-Seq and RNA-Seq further narrowed the candidate range. As a result, a total of three quantitative trait loci (QTLs) and nine candidate genes related to plant height were obtained on chromosomes I and IX. Based on the functional prediction of the candidate genes, we propose a hypothetical mechanism for the formation of millet dwarfing, in which, metabolism and MAPK signaling play important roles in the formation of foxtail millet plant height.

Epidermal Patterning Factor 2-like (McEPFL2): A Putative Candidate for the Continuous Ridge (cr) Fruit Skin Locus in Bitter Gourd (Momordica charantia L.)

Bitter gourd (Momordica charantia L.) is an economically important vegetable and medicinal crop in many Asian countries. Limited work has been conducted in understanding the genetic basis of horticulturally important traits in bitter gourd. Bitter gourd is consumed primarily for its young, immature fruit, and fruit appearance plays an important role in market acceptability. One such trait is the ridges on the fruit skin. In the present study, molecular mapping of a locus underlying fruit ridge continuity was conducted. Genetic analysis in segregating populations, derived from the crosses between two inbred lines Y1 with continuous ridges (CR) and Z-1-4 with discontinuous ridges (DCR), suggested that CR was controlled by a single recessive gene (cr). High-throughput genome sequencing of CR and DCR bulks combined with high-resolution genetic mapping in an F₂ population delimited cr into a 108 kb region with 16 predicted genes. Sequence variation analysis and expression profiling supported the epidermal patterning factor 2-like (McEPFL2) gene as the best candidate of the cr locus. A 1 bp deletion in the first exon of McEPFL2 in Y1 which would result in a truncated McEPFL2 protein may be the causal polymorphism for the phenotypic difference between Y1 and Z-1-4. The association of this 1 bp deletion with CR was further supported by gDNA sequencing of McEPFL2 among 31 bitter gourd accessions. This work provides a foundation for understanding the genetic and molecular control of fruit epidermal pattering and development, which also facilitates marker-assisted selection in bitter melon breeding.

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.

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.

Frequent gain- and loss-of-function mutations of the BjMYB113 gene accounted for leaf color variation in Brassica juncea

BACKGROUND

Mustard (Brassica juncea) is an important economic vegetable, and some cultivars have purple leaves and accumulate more anthocyanins than the green. The genetic and evolution of purple trait in mustard has not been well studied.

RESULT

In this study, free-hand sections and metabolomics showed that the purple leaves of mustard accumulated more anthocyanins than green ones. The gene controlling purple leaves in mustard, Mustard Purple Leaves (MPL), was genetically mapped and a MYB113-like homolog was identified as the candidate gene. We identified three alleles of the MYB113-like gene, BjMYB113a from a purple cultivar, BjMYB113b and BjMYB113c from green cultivars. A total of 45 single nucleotide polymorphisms (SNPs) and 8 InDels were found between the promoter sequences of the purple allele BjMYB113a and the green allele BjMYB113b. On the other hand, the only sequence variation between the purple allele BjMYB113a and the green allele BjMYB113c is an insertion of 1,033-bp fragment in the 3'region of BjMYB113c. Transgenic assay and promoter activity studies showed that the polymorphism in the promoter region was responsible for the up-regulation of the purple allele BjMYB113a and high accumulation of anthocyanin in the purple cultivar. The up-regulation of BjMYB113a increased the expression of genes in the anthocyanin biosynthesis pathway including BjCHS, BjF3H, BjF3'H, BjDFR, BjANS and BjUGFT, and consequently led to high accumulation of anthocyanin. However, the up-regulation of BjMYB113 was compromised by the insertion of 1,033-bp in 3'region of the allele BjMYB113c.

CONCLUSIONS

Our results contribute to a better understanding of the genetics and evolution of the BjMYB113 gene controlling purple leaves and provide useful information for further breeding programs of mustard.

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.

Mapping of Genetic Locus for Leaf Trichome Formation in Chinese Cabbage Based on Bulked Segregant Analysis

Chinese cabbage is a leafy vegetable, and its leaves are the main edible organs. The formation of trichomes on the leaves can significantly affect its taste, so studying this phenomenon is of great significance for improving the quality of Chinese cabbage. In this study, two varieties of Chinese cabbage, W30 with trichome leaves and 082 with glabrous leaves, were crossed to generate F₁ and F₁ plants, which were self-fertilized to develop segregating populations with trichome or glabrous morphotypes. The two bulks of the different segregating populations were used to conduct bulked segregant analysis (BSA). A total of 293.4 M clean reads were generated from the samples, and plants from the trichome leaves (AL) bulk and glabrous leaves (GL) bulk were identified. Between the two DNA pools generated from the trichome and glabrous plants, 55,048 SNPs and 272 indels were generated. In this study, three regions (on chromosomes 6, 10 and scaffold000100) were identified, and the annotation revealed three candidate genes that may participate in the formation of leaf trichomes. These findings suggest that the three genes-Bra025087 encoding a cyclin family protein, Bra035000 encoding an ATP-binding protein/kinase/protein kinase/protein serine/threonine kinase and Bra033370 encoding a WD-40 repeat family protein-influence the formation of trichomes by participating in trichome morphogenesis (GO: 0010090). These results demonstrate that BSA can be used to map genes associated with traits and provide new insights into the molecular mechanism of leafy trichome formation in Chinese cabbage.

Comprehensive Transcriptome Analyses Reveal Candidate Genes for Variation in Seed Size/Weight During Peanut (Arachis hypogaea L.) Domestication

Seed size/weight, a key domestication trait, is also an important selection target during peanut breeding. However, the mechanisms that regulate peanut seed development are unknown. We re-sequenced 12 RNA samples from developing seeds of two cultivated peanut accessions (Lines 8106 and 8107) and wild Arachis monticola at 15, 30, 45, and 60 days past flowering (DPF). Transcriptome analyses showed that ∼36,000 gene loci were expressed in each of the 12 RNA samples, with nearly half exhibiting moderate (2 ≤ FPKM < 10) expression levels. Of these genes, 12.2% (4,523) were specifically expressed during seed development, mainly at 15 DPF. Also, ∼12,000 genes showed significant differential expression at 30, 45, and/or 60 DPF within each of the three peanut accessions, accounting for 31.8-34.1% of the total expressed genes. Using a method that combined comprehensive transcriptome analysis and previously mapped QTLs, we identified several candidate genes that encode transcription factor TGA7, topless-related protein 2, IAA-amino acid hydrolase ILR1-like 5, and putative pentatricopeptide repeat-containing (PPR) protein. Based on sequence variations identified in these genes, SNP markers were developed and used to genotype both 30 peanut landraces and a genetic segregated population, implying that EVM0025654 encoding a PPR protein may be associated with the increased seed size/weight of the cultivated accessions in comparison with the allotetraploid wild peanut. Our results provide additional knowledge for the identification and functional research into candidate genes responsible for the seed size/weight phenotype in peanut.