Single-cell transcriptomic analysis of flowering regulation and vernalization in Chinese cabbage shoot apex.

In Chinese cabbage development the interplay between shoot apex activity and vernalization is pivotal for flowering timing. The intricate relationship between various cell types in the shoot apex meristem and their roles in regulating flowering gene expression in Chinese cabbage is not yet fully understood. A thorough analysis of single-cell types in the Chinese cabbage shoot apex and their influence on flowering genes and vernalization is essential for deeper insight. Our study first established a single-cell transcriptomic atlas of Chinese cabbage after 25 days of non-vernalization. Analyzing 19 602 single cells, we differentiated them into 15 distinct cell clusters using established marker genes. We found that key genes in shoot apex development and flowering were primarily present in shoot meristematic cells (SMCs), companion cells (CCs), and mesophyll cells (MCs). MADS-box protein FLOWERING LOCUS C 2 (BrFLC2), a gene suppressing flowering, was observed in CCs, mirroring patterns found in Arabidopsis. By mapping developmental trajectories of SMCs, CCs, and MCs, we elucidated the evolutionary pathways of crucial genes in shoot apex development and flowering. The creation of a single-cell transcriptional atlas of the Chinese cabbage shoot apex under vernalization revealed distinct alterations in the expression of known flowering genes, such as VERNALIZATION INSENSITIVE 3 (VIN3), VERNALIZATION 1 (VRN1), VERNALIZATION 2 (VRN2), BrFLC, and FLOWERING LOCUS T (FT), which varied by cell type. Our study underscores the transformative impact of single-cell RNA sequencing (scRNA-seq) for unraveling the complex differentiation and vernalization processes in the Chinese cabbage shoot apex. These insights are pivotal for enhancing breeding strategies and cultivation management of this vital vegetable.

Genome-wide identification of Brassinosteroid insensitive 1-associated receptor kinase 1 genes and expression analysis in response to pathogen infection in cucumber (Cucumis sativus L.).

BACKGROUND: BAK1 (Brassinosteroid insensitive 1-associated receptor kinase 1) plays an important role in disease resistance in plants. However, the function of BAK1 family in cucumber and the decisive genes for disease-resistance remain elusive. RESULTS: Here, we identified 27 CsBAK1s in cucumber, and classified them into five subgroups based on phylogenetic analysis and gene structure. CsBAK1s in the same subgroup shared the similar motifs, but different gene structures. Cis-elements analysis revealed that CsBAK1s might respond to various stress and growth regulation. Three segmentally duplicated pairwise genes were identified in cucumber. In addition, Ka/Ks analysis indicated that CsBAK1s were under positive selection during evolution. Tissue expression profile showed that most CsBAK1s in Subgroup II and IV showed constitutive expression, members in other subgroups showed tissue-specific expression. To further explore whether CsBAK1s were involved in the resistance to pathogens, the expression patterns of CsBAK1s to five pathogens (gummy stem blight, powdery mildew, downy mildew, grey mildew, and fusarium wilt) reveled that different CsBAK1s had specific roles in different pathogen infections. The expression of CsBAK1-14 was induced/repressed significantly by five pathogens, CsBAK1-14 might play an important role in disease resistance in cucumber. CONCLUSIONS: 27 BAK1 genes were identified in cucumber from a full perspective, which have important functions in pathogen infection. Our study provided a theoretical basis to further clarify the function of BAK1s to disease resistance in cucumber.

Genome-wide association analysis reveal candidate genes and haplotypes related to root weight in cucumber (Cucumis sativus L.).

Background: The plant root system is critical for the absorption of water and nutrients, and have a direct influence on growth and yield. In cucumber, a globally consumed crop, the molecular mechanism of root development remains unclear, and this has implications for developing stress tolerant varieties. This study sought to determine the genetic patterns and related genes of cucumber root weight. A core cucumber germplasms population was used to do the GWAS analysis in three environments. Results: Here, we investigated four root-weight related traits including root fresh weight (RFW), root dry weight (RDW), ratio of root dry weight to root fresh weight (RDFW) and the comprehensive evaluation index, D-value of root weight (DRW) deduced based on the above three traits for the core germplasm of the cucumber global repository. According to the D-value, we identified 21 and 16 accessions with light and heavy-root, respectively. We also found that the East Asian ecotype accessions had significantly heavier root than other three ecotypes. The genome-wide association study (GWAS) for these four traits reveals that 4 of 10 significant loci (gDRW3.1, gDRW3.2, gDRW4.1 and gDRW5.1) were repeatedly detected for at least two traits. Further haplotype and expression analysis for protein-coding genes positioned within these 4 loci between light and heavy-root accessions predicted five candidate genes (i.e., Csa3G132020 and Csa3G132520 both encoding F-box protein PP2-B1 for gDRW3.1, Csa3G629240 encoding a B-cell receptor-associated protein for gDRW3.2, Csa4G499330 encodes a GTP binding protein for gDRW4.1, and Csa5G286040 encodes a proteinase inhibitor for gDRW5.1). Conclusions: We conducted a systematic analysis of the root genetic basis and characteristics of cucumber core germplasms population. We detected four novel loci, which regulate the root weight in cucumber. Our study provides valuable candidate genes and haplotypes for the improvement of root system in cucumber breeding.

CsMLO8/11 are required for full susceptibility of cucumber stem to powdery mildew and interact with CsCRK2 and CsRbohD.

Powdery mildew (PM) is one of the most destructive diseases that threaten cucumber production globally. Efficient breeding of novel PM-resistant cultivars will require a robust understanding of the molecular mechanisms of cucumber resistance against PM. Using a genome-wide association study, we detected a locus significantly correlated with PM resistance in cucumber stem, pm-s5.1. A 1449-bp insertion in the CsMLO8 coding region at the pm-s5.1 locus resulted in enhanced stem PM resistance. Knockout mutants of CsMLO8 and CsMLO11 generated by CRISPR/Cas9 both showed improved PM resistance in the stem, hypocotyl, and leaves, and the double mutant mlo8mlo11 displayed even stronger resistance. We found that reactive oxygen species (ROS) accumulation was higher in the stem of these mutants. Protein interaction assays suggested that CsMLO8 and CsMLO11 could physically interact with CsRbohD and CsCRK2, respectively. Further, we showed that CsMLO8 and CsCRK2 competitively interact with the C-terminus of CsRbohD to affect CsCRK2-CsRbohD module-mediated ROS production during PM defense. These findings provide new insights into the understanding of CsMLO proteins during PM defense responses.

Natural variation in STAYGREEN contributes to low-temperature tolerance in cucumber.

Low-temperature (LT) stress threatens cucumber production globally; however, the molecular mechanisms underlying LT tolerance in cucumber remain largely unknown. Here, using a genome-wide association study (GWAS), we found a naturally occurring single nucleotide polymorphism (SNP) in the STAYGREEN (CsSGR) coding region at the gLTT5.1 locus associated with LT tolerance. Knockout mutants of CsSGR generated by clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated nuclease 9 exhibit enhanced LT tolerance, in particularly, increased chlorophyll (Chl) content and reduced reactive oxygen species (ROS) accumulation in response to LT. Moreover, the C-repeat Binding Factor 1 (CsCBF1) transcription factor can directly activate the expression of CsSGR. We demonstrate that the LT-sensitive haplotype CsSGRHapA , but not the LT-tolerant haplotype CsSGRHapG could interact with NON-YELLOW COLORING 1 (CsNYC1) to mediate Chl degradation. Geographic distribution of the CsSGR haplotypes indicated that the CsSGRHapG was selected in cucumber accessions from high latitudes, potentially contributing to LT tolerance during cucumber cold-adaptation in these regions. CsSGR mutants also showed enhanced tolerance to salinity, water deficit, and Pseudoperonospora cubensis, thus CsSGR is an elite target gene for breeding cucumber varieties with broad-spectrum stress tolerance. Collectively, our findings provide new insights into LT tolerance and will ultimately facilitate cucumber molecular breeding.

Investigation of the Role of a Zinc Uptake Regulator (Zur) in the Virulence of Pectobacterium odoriferum.

Pectobacterium spp. infect many horticultural crops worldwide and lead to serious crop losses. Zinc-uptake-regulator (Zur) proteins are present widely in prokaryotes and play an important role in pathogenicity. To uncover the role of Zur in P. odoriferum, we constructed mutant (ΔZur) and overexpression [Po (Zur)] strains of a Zur, and a virulence assay showed that the Po (Zur) was of significantly lower virulence, while the ΔZur displayed significantly increased virulence on Chinese cabbage compared to their respective control strains, wild-type P. odoriferum (Po WT) and P. odoriferum harboring an empty vector (Po (EV)) (p < 0.05). The growth curves of the ΔZur and Po (Zur) showed no obvious differences from those of the control strains. Comparative transcriptome analysis showed that Zur overexpression in P. odoriferum induced differentially expressed genes (DEGs) related to flagellum and cell motility, while mutating Zur resulted in DEGs mainly corresponding to divalent-metal-ion transport and membrane transport. Phenotypic experiments on the Po (Zur) showed that flagellum numbers and cell motility were reduced in comparison with the control, while those of the ΔZur did not change. Collectively, these results show that the Zur negatively regulates the virulence of P. odoriferum and might function via a dual mechanism dependent on dose.

GWAS reveals novel loci and identifies a pentatricopeptide repeat-containing protein (CsPPR) that improves low temperature germination in cucumber.

Low temperatures (LTs) negatively affect the percentage and rate of cucumber (Cucumis sativus L.) seed germination, which has deleterious effects on yield. Here, a genome-wide association study (GWAS) was used to identify the genetic loci underlying low temperature germination (LTG) in 151 cucumber accessions that represented seven diverse ecotypes. Over two years, phenotypic data for LTG i.e., relative germination rate (RGR), relative germination energy (RGE), relative germination index (RGI) and relative radical length (RRL), were collected in two environments, and 17 of the 151 accessions were found to be highly cold tolerant using cluster analysis. A total of 1,522,847 significantly associated single-nucleotide polymorphism (SNP) were identified, and seven loci associated with LTG, on four chromosomes, were detected: gLTG1.1, gLTG1.2, gLTG1.3, gLTG4.1, gLTG5.1, gLTG5.2, and gLTG6.1 after resequencing of the accessions. Of the seven loci, three, i.e., gLTG1.2, gLTG4.1, and gLTG5.2, showed strong signals that were consistent over two years using the four germination indices, and are thus strong and stable for LTG. Eight candidate genes associated with abiotic stress were identified, and three of them were potentially causal to LTG: CsaV3_1G044080 (a pentatricopeptide repeat-containing protein) for gLTG1.2, CsaV3_4G013480 (a RING-type E3 ubiquitin transferase) for gLTG4.1, and CsaV3_5G029350 (a serine/threonine-protein kinase) for gLTG5.2. The function for CsPPR (CsaV3_1G044080) in regulating LTG was confirmed, as Arabidopsis lines ectopically expressing CsPPR showed higher germination and survival rates at 4°C compared to the wild-type, which preliminarily illustrates that CsPPR positively regulates cucumber cold tolerance at the germination stage. This study will provide insights into cucumber LT-tolerance mechanisms and further promote cucumber breeding development.

The roles of microRNAs in horticultural plant disease resistance.

The development of the horticultural industry is largely limited by disease and excessive pesticide application. MicroRNAs constitute a major portion of the transcriptomes of eukaryotes. Various microRNAs have been recognized as important regulators of the expression of genes involved in essential biological processes throughout the whole life cycle of plants. Recently, small RNA sequencing has been applied to study gene regulation in horticultural plants. In this review, we summarize the current understanding of the biogenesis and contributions of microRNAs in horticultural plant disease resistance. These microRNAs may potentially be used as genetic resources for improving disease resistance and for molecular breeding. The challenges in understanding horticultural plant microRNA biology and the possibilities to make better use of these horticultural plant gene resources in the future are discussed in this review.

Genomic analyses of rice bean landraces reveal adaptation and yield related loci to accelerate breeding.

Rice bean (Vigna umbellata) is an underexploited domesticated legume crop consumed for dietary protein in Asia, yet little is known about the genetic diversity of this species. Here, we present a high-quality reference genome for a rice bean landrace (FF25) built using PacBio long-read data and a Hi-C chromatin interaction map, and assess the phylogenetic position and speciation time of rice bean within the Vigna genus. We sequence 440 landraces (two core collections), and GWAS based on data for growth sites at three widely divergent latitudes reveal loci associated with flowering and yield. Loci harboring orthologs of FUL (FRUITFULL), FT (FLOWERING LOCUS T), and PRR3 (PSEUDO-RESPONSE REGULATOR 3) contribute to the adaptation of rice bean from its low latitude center of origin towards higher latitudes, and the landraces which pyramid early-flowering alleles for these loci display maximally short flowering times. We also demonstrate that copy-number-variation for VumCYP78A6 can regulate seed-yield traits. Intriguingly, 32 landraces collected from a mountainous region in South-Central China harbor a recently acquired InDel in TFL1 (TERMINAL FLOWER1) affecting stem determinacy; these materials also have exceptionally high values for multiple human-desired traits and could therefore substantially advance breeding efforts to improve rice bean.

InDels Identification and Association Analysis with Spike and Awn Length in Chinese Wheat Mini-Core Collection.

Diversity surveys of germplasm are important for gaining insight into the genomic basis for crop improvement; especially InDels, which are poorly understood in hexaploid common wheat. Here, we describe a map of 89,923 InDels from exome sequencing of 262 accessions of a Chinese wheat mini-core collection. Population structure analysis, principal component analysis and selective sweep analysis between landraces and cultivars were performed. Further genome-wide association study (GWAS) identified five QTL (Quantitative Trait Loci) that were associated with spike length, two of them, on chromosomes 2B and 6A, were detected in 10 phenotypic data sets. Assisted with RNA-seq data, we identified 14 and 21 genes, respectively that expressed in spike and rachis within the two QTL regions that can be further investigated for candidate genes discovery. Moreover, InDels were found to be associated with awn length on chromosomes 5A, 6B and 4A, which overlapped with previously reported genetic loci B1 (Tipped 1), B2 (Tipped 2) and Hd (Hooded). One of the genes TaAGL6 that was previously shown to affect floral organ development was found at the B2 locus to affect awn length development. Our study shows that trait-associated InDels may contribute to wheat improvement and may be valuable molecular markers for future wheat breeding.

Transcriptome Analysis of Developing Wheat Grains at Rapid Expanding Phase Reveals Dynamic Gene Expression Patterns.

Grain development, as a vital process in the crop's life cycle, is crucial for determining crop quality and yield. The wheat grain expanding phase is the early process involving the rapid morphological changes and initiation of grain filling. However, little is known about the molecular basis of grain development at this stage. Here, we provide a time-series transcriptome profile of developing wheat grain at 0, 2, 4, 6, 8, and 10 days after pollination of the wheat landrace Chinese Spring. A total of 26,892 differentially expressed genes, including 1468 transcription factors, were found between adjacent time points. Co-expression cluster analysis and Gene Ontology enrichment revealed dynamic expressions of cell division and starch biosynthesis related structural genes and transcription factors. Moreover, diverse, differential and drastically varied expression trends of the key genes related to hormone metabolism were identified. Furthermore, ~30% of triads showed unbalanced expression patterns enriching for genes in multiple pivotal metabolic pathways. Hormone metabolism related genes, such as YUC10 (YUCCA flavin-containing monooxygenase 10), AOS2 (allene oxide synthase 2), CYP90D2 (cytochrome P450 90D2), and CKX1 (cytokinin dehydrogenase 1), were dominantly contributed by A or D homoeologs of the triads. Our study provided a systematic picture of transcriptional regulation of wheat grains at the early grain expanding phase which should deepen our understanding of wheat grain development and help in wheat yield improvement.

The wheat AGL6-like MADS-box gene is a master regulator for floral organ identity and a target for spikelet meristem development manipulation.

The AGAMOUS-LIKE6 (AGL6)-like genes are ancient MADS-box genes and are functionally studied in a few model plants. The knowledge of these genes in wheat remains limited. Here, by studying a 'double homoeolog mutant' of the AGL6 gene in tetraploid wheat, we showed that AGL6 was required for the development of all four whorls of floral organs with dosage-dependent effect on floret fertility. Yeast two-hybrid analyses detected interactions of AGL6 with all classes of MADS-box proteins in the ABCDE model for floral organ development. AGL6 was found to interact with several additional proteins, including the G protein ß and γ (DEP1) subunits. Analysis of the DEP1-B mutant showed a significant reduction in spikelet number per spike in tetraploid wheat, while overexpression of AGL6 in common wheat increased the spikelet number per spike and hence the grain number per spike. RNA-seq analysis identified the regulation of several meristem activity genes by AGL6, such as FUL2 and TaMADS55. Our work therefore extensively updated the wheat ABCDE model and proposed an alternative approach to improve wheat grain yield by manipulating the AGL6 gene.

Population-scale peach genome analyses unravel selection patterns and biochemical basis underlying fruit flavor.

A narrow genetic basis in modern cultivars and strong linkage disequilibrium in peach (Prunus persica) has restricted resolution power for association studies in this model fruit species, thereby limiting our understanding of economically important quality traits including fruit flavor. Here, we present a high-quality genome assembly for a Chinese landrace, Longhua Shui Mi (LHSM), a representative of the Chinese Cling peaches that have been central in global peach genetic improvement. We also map the resequencing data for 564 peach accessions to this LHSM assembly at an average depth of 26.34× per accession. Population genomic analyses reveal a fascinating history of convergent selection for sweetness yet divergent selection for acidity in eastern vs. western modern cultivars. Molecular-genetics and biochemical analyses establish that PpALMT1 (aluminum-activated malate transporter 1) contributes to their difference of malate content and that increases fructose content accounts for the increased sweetness of modern peach fruits, as regulated by PpERDL16 (early response to dehydration 6-like 16). Our study illustrates the strong utility of the genomics resources for both basic and applied efforts to understand and exploit the genetic basis of fruit quality in peach.

Genome structure variation analyses of peach reveal population dynamics and a 1.67 Mb causal inversion for fruit shape.

BACKGROUND: Structural variations (SVs), a major resource of genomic variation, can have profound consequences on phenotypic variation, yet the impacts of SVs remain largely unexplored in crops. RESULTS: Here, we generate a high-quality de novo genome assembly for a flat-fruit peach cultivar and produce a comprehensive SV map for peach, as a high proportion of genomic sequence is occupied by heterozygous SVs in the peach genome. We conduct population-level analyses that indicate SVs have undergone strong purifying selection during peach domestication, and find evidence of positive selection, with a significant preference for upstream and intronic regions during later peach improvement. We perform a SV-based GWAS that identifies a large 1.67-Mb heterozygous inversion that segregates perfectly with flat-fruit shape. Mechanistically, this derived allele alters the expression of the PpOFP2 gene positioned near the proximal breakpoint of the inversion, and we confirm in transgenic tomatoes that PpOFP2 is causal for flat-fruit shape. CONCLUSIONS: Thus, beyond introducing new genomics resources for peach research, our study illustrates how focusing on SV data can drive basic functional discoveries in plant science.

The Battle to Sequence the Bread Wheat Genome: A Tale of the Three Kingdoms.

In the year 2018, the world witnessed the finale of the race to sequence the genome of the world's most widely grown crop, the common wheat. Wheat has been known to bear a notoriously large and complicated genome of a polyploidy nature. A decade competition to sequence the wheat genome initiated with a single consortium of multiple countries, taking a conventional strategy similar to that for sequencing Arabidopsis and rice, became ferocious over time as both sequencing technologies and genome assembling methodologies advanced. At different stages, multiple versions of genome sequences of the same variety (e.g., Chinese Spring) were produced by several groups with their special strategies. Finally, 16 years after the rice genome was finished and 9 years after that of maize, the wheat research community now possesses its own reference genome. Armed with these genomics tools, wheat will reestablish itself as a model for polyploid plants in studying the mechanisms of polyploidy evolution, domestication, genetic and epigenetic regulation of homoeolog expression, as well as defining its genetic diversity and breeding on the genome level. The enhanced resolution of the wheat genome should also help accelerate development of wheat cultivars that are more tolerant to biotic and/or abiotic stresses with better quality and higher yield.

A transcriptomic view of the ability of nascent hexaploid wheat to tolerate aneuploidy.

BACKGROUND: In contrast to most animal species, polyploid plant species are quite tolerant of aneuploidy. Here, the global transcriptome of four aneuploid derivatives of a synthetic hexaploid wheat line was acquired, with the goal of characterizing the relationship between gene copy number and transcript abundance. RESULTS: For most of the genes mapped to the chromosome involved in aneuploidy, the abundance of transcripts reflected the gene copy number. Aneuploidy had a greater effect on the strength of transcription of genes mapped to the chromosome present in a noneuploid dose than on that of genes mapped elsewhere in the genome. Overall, changing the copy number of one member of a homeologous set had little effect on the abundance of transcripts generated from the set of homeologs as a whole, consistent with the tolerance of aneuploidy exhibited by allopolyploids, whether in the form of a chromosomal deficit (monosomy) or chromosomal excess (trisomy). CONCLUSIONS: Our findings shed new light on the genetic regulation of homeoallele transcription and contribute to a deeper understanding of allopolyploid genome evolution, with implications for the breeding of polyploid crops.

A high-throughput BAC end analysis protocol (BAC-anchor) for profiling genome assembly and physical mapping.

Traditional approaches for sequencing insertion ends of bacterial artificial chromosome (BAC) libraries are laborious and expensive, which are currently some of the bottlenecks limiting a better understanding of the genomic features of auto- or allopolyploid species. Here, we developed a highly efficient and low-cost BAC end analysis protocol, named BAC-anchor, to identify paired-end reads containing large internal gaps. Our approach mainly focused on the identification of high-throughput sequencing reads carrying restriction enzyme cutting sites and searching for large internal gaps based on the mapping locations of both ends of the reads. We sequenced and analysed eight libraries containing over 3 200 000 BAC end clones derived from the BAC library of the tetraploid potato cultivar C88 digested with two restriction enzymes, Cla I and Mlu I. About 25% of the BAC end reads carrying cutting sites generated a 60-100 kb internal gap in the potato DM reference genome, which was consistent with the mapping results of Sanger sequencing of the BAC end clones and indicated large differences between autotetraploid and haploid genotypes in potato. A total of 5341 Cla I- and 165 Mlu I-derived unique reads were distributed on different chromosomes of the DM reference genome and could be used to establish a physical map of target regions and assemble the C88 genome. The reads that matched different chromosomes are especially significant for the further assembly of complex polyploid genomes. Our study provides an example of analysing high-coverage BAC end libraries with low sequencing cost and is a resource for further genome sequencing studies.

Genetic mapping of folate QTLs using a segregated population in maize.

As essential B vitamin for humans, folates accumulation in edible parts of crops, such as maize kernels, is of great importance for human health. But its breeding is always limited by the prohibitive cost of folate profiling. The molecular breeding is a more executable and efficient way for folate fortification, but is limited by the molecular knowledge of folate regulation. Here we report the genetic mapping of folate quantitative trait loci (QTLs) using a segregated population crossed by two maize lines, one high in folate (GEMS31) and the other low in folate (DAN3130). Two folate QTLs on chromosome 5 were obtained by the combination of F2 whole-exome sequencing and F3 kernel-folate profiling. These two QTLs had been confirmed by bulk segregant analysis using F6 pooled DNA and F7 kernel-folate profiling, and were overlapped with QTLs identified by another segregated population. These two QTLs contributed 41.6% of phenotypic variation of 5-formyltetrahydrofolate, the most abundant storage form among folate derivatives in dry maize grains, in the GEMS31×DAN3130 population. Their fine mapping and functional analysis will reveal details of folate metabolism, and provide a basis for marker-assisted breeding aimed at the enrichment of folates in maize kernels.