Arabidopsis REM16 acts as a B3 domain transcription factor to promote flowering time via directly binding to the promoters of SOC1 and FT

Large-scale transcriptomic and genomic analyses reveal a novel functional gene SERPINB6 for chicken carcass traits

BACKGROUND

Carcass traits are crucial indicators of meat production efficiency. However, the molecular regulatory mechanisms associated with these traits remain unclear.

RESULTS

In this study, we conducted comprehensive transcriptomic and genomic analyses on 399 Tiannong partridge chickens to identify key genes and variants associated with carcass traits and to elucidate the underlying regulatory mechanisms. Based on association analyses with the elastic net (EN) model, we identified 12 candidate genes (AMY1A, AP3B2, CEBPG, EEF2, EIF4EBP1, FGFR1, FOXD3, GOLM1, LOC107052698, PABPC1, SERPINB6 and TBC1D16) for 4 carcass-related traits, namely live weight, dressed weight, eviscerated weight, and breast muscle weight. SERPINB6 was identified as the only overlapping gene by 3 analyses, EN model analysis, weighted gene co-expression network analysis and differential expression analysis. Cell-level experiments confirmed that SERPINB6 promotes the proliferation of chicken DF1 cells and primary myoblasts. Further expression genome-wide association study and association analysis indicated that rs317934171 is the critical site that enhances SERPINB6 expression. Furthermore, a dual-luciferase reporter assay proved that gga-miR-1615 targets the 3'UTR of SERPINB6.

CONCLUSIONS

Collectively, our findings reveal that SERPINB6 serves as a novel gene for chicken carcass traits by promoting fibroblast and myoblast proliferation. Additionally, the downstream variant rs317934171 regulates SERPINB6 expression. These results identify a new target gene and molecular marker for the molecular mechanisms of chicken carcass traits.

The Arabidopsis RTH plays an important role in regulation of iron (Fe) absorption and transport

KEY MESSAGE

RTH may activate Fe assimilation related genes to promote Fe absorption, transport and accumulation in Arabidopsis. Iron (Fe) is an important nutrient element. The Fe absorption and transport in plants are well investigated over the past decade. Our previous work indicated that RTE1-HOMOLOG (RTH), the homologous gene of reversion-to-ethylene sensitivity 1 (RTE1), plays a role in ethylene signaling pathway. However, its function in Fe absorption and transport is largely unknown. In the present study, we found that RTH was expressed in absorptive tissue and conducting tissue, including root hairs, root vascular bundle, and leaf veins. Under high Fe concentration, the seedling growth of rth-1 mutant was better, while the RTH overexpression lines were retarded compared to the wild type (Col-0). When treated with EDTA-Fe3+ (400 μM), the chlorophyll content and ion leakage rate were higher and lower in rth-1 than those of Col-0, respectively. By contrast, the chlorophyll contents and ion leakage rates of RTH overexpression lines were decreased and hastened compared with Col-0, respectively. Fe measurement indicated that the Fe contents of rth-1 were lower than those of Col-0, whereas those of RTH overexpression lines were comparably higher. Gene expression analysis revealed that Fe absorption and transport genes AHA2, IRT1, FIT, FPN1, and YSL1 decreased in rth-1 but increased in RTH overexpression lines compared with Col-0. Additionally, Y2H (yeast two-hybrid) and BiFC (bimolecular fluorescence complementation) assays showed that RTH can physically interact with hemoglobin 1 (HB1) and HB2. All these findings suggest that RTH may play an important role in regulation of Fe absorption, transport, and accumulation in Arabidopsis.

Comparative Genomic Analysis of PEBP Genes in Cucurbits Explores the Interactors of Cucumber CsPEBPs Related to Flowering Time

The family of phosphatidylethanolamine-binding proteins (PEBPs) participates in various plant biological processes, mainly flowering regulation and seed germination. In cucurbit crops, several PEBP genes have been recognized to be responsible for flowering time. However, the investigation of PEBP family members across the genomes of cucurbit species has not been reported, and their conservation and divergence in structure and function remain largely unclear. Herein, PEBP genes were identified from seven cucurbit crops and were used to perform a comparative genomics analysis. The cucurbit PEBP proteins could be classified into MFT, FT, TFL, and PEBP clades, and further, the TFL clade was divided into BFT-like, CEN-like, and TFL1-like subclades. The MFT-like, FT-like, and TFL-like proteins were clearly distinguished by a critical amino acid residue at the 85th position of the Arabidopsis FT protein. In gene expression analysis, CsaPEBP1 was highly expressed in flowers, and its expression levels in females and males were 70.5 and 89.2 times higher, respectively, than those in leaves. CsaPEBP5, CsaPEBP6, and CsaPEBP7 were specifically expressed in male flowers, with expression levels 58.1, 17.3, and 15.7 times higher, respectively, than those of leaves. At least five CsaPEBP genes exhibited the highest expression during the later stages of corolla opening. Through clustering of time-series-based RNA-seq data, several potential transcription factors (TFs) interacting with four CsaPEBPs were identified during cucumber corolla opening. Because of the tandem repeats of binding sites in promoters, NF-YB (Csa4G037610) and GATA (Csa7G64580) TFs appeared to be better able to regulate the CsaPEBP2 and CsaPEBP5 genes, respectively. This study would provide helpful information for further investigating the roles of PEBP genes and their interacting TFs in growth and development processes, such as flowering time regulation in cucurbit crops.

Plant-TFClass: a structural classification for plant transcription factors

Transcription factors (TFs) bind DNA at specific sequences to regulate gene expression. This universal process is achieved via their DNA-binding domain (DBD). In mammals, the vast diversity of DBD structural conformations and the way in which they contact DNA has been used to organize TFs in the TFClass hierarchical classification. However, the numerous DBD types present in plants but absent from mammalian genomes were missing from this classification. We reviewed DBD 3D structures and models available for plant TFs to classify most of the 56 recognized plant TF types within the TFClass framework. This extended classification adds eight new classes and 37 new families corresponding to DBD structures absent in mammals. Plant-TFClass provides a unique resource for TF comparison across families and organisms.

Comparative and population genomics of buckwheat species reveal key determinants of flavor and fertility

Common buckwheat (Fagopyrum esculentum) is an ancient crop with a world-wide distribution. Due to its excellent nutritional quality and high economic and ecological value, common buckwheat is becoming increasingly important throughout the world. The availability of a high-quality reference genome sequence and population genomic data will accelerate the breeding of common buckwheat, but the high heterozygosity due to the outcrossing nature has greatly hindered the genome assembly. Here we report the assembly of a chromosome-scale high-quality reference genome of F. esculentum var. homotropicum, a homozygous self-pollinating variant of common buckwheat. Comparative genomics revealed that two cultivated buckwheat species, common buckwheat (F. esculentum) and Tartary buckwheat (F. tataricum), underwent metabolomic divergence and ecotype differentiation. The expansion of several gene families in common buckwheat, including FhFAR genes, is associated with its wider distribution than Tartary buckwheat. Copy number variation of genes involved in the metabolism of flavonoids is associated with the difference of rutin content between common and Tartary buckwheat. Furthermore, we present a comprehensive atlas of genomic variation based on whole-genome resequencing of 572 accessions of common buckwheat. Population and evolutionary genomics reveal genetic variation associated with environmental adaptability and floral development between Chinese and non-Chinese cultivated groups. Genome-wide association analyses of multi-year agronomic traits with the content of flavonoids revealed that Fh05G014970 is a potential major regulator of flowering period, a key agronomic trait controlling the yield of outcrossing crops, and that Fh06G015130 is a crucial gene underlying flavor-associated flavonoids. Intriguingly, we found that the gene translocation and sequence variation of FhS-ELF3 contribute to the homomorphic self-compatibility of common buckwheat. Collectively, our results elucidate the genetic basis of speciation, ecological adaptation, fertility, and unique flavor of common buckwheat, and provide new resources for future genomics-assisted breeding of this economically important crop.

Epigenetic silencing of callose synthase by VIL1 promotes bud-growth transition in lily bulbs

In plants, restoring intercellular communication is required for cell activity in buds during the growth transition from slow to fast growth after dormancy release. However, the epigenetic regulation of this phenomenon is far from understood. Here we demonstrate that lily VERNALIZATION INSENSITIVE 3-LIKE 1 (LoVIL1) confers growth transition by mediating plasmodesmata opening via epigenetic repression of CALLOSE SYNTHASE 3 (LoCALS3). Moreover, we found that a novel transcription factor, NUCLEAR FACTOR Y, SUBUNIT A7 (LoNFYA7), is capable of recruiting the LoVIL1-Polycomb Repressive Complex 2 (PRC2) and enhancing H3K27me3 at the LoCALS3 locus by recognizing the CCAAT cis-element (Cce) of its promoter. The LoNFYA7-LoVIL1 module serves as a key player in orchestrating the phase transition from slow to fast growth in lily bulbs. These studies also indicate that LoVIL1 is a suitable marker for the bud-growth-transition trait following dormancy release in lily cultivars.

Comparatively Evolution and Expression Analysis of GRF Transcription Factor Genes in Seven Plant Species

Growth regulatory factors (GRF) are plant-specific transcription factors that play pivotal roles in growth and various abiotic stresses regulation. However, adaptive evolution of GRF gene family in land plants are still being elucidated. Here, we performed the evolutionary and expression analysis of GRF gene family from seven representative species. Extensive phylogenetic analyses and gene structure analysis revealed that the number of genes, QLQ domain and WRC domain identified in higher plants was significantly greater than those identified in lower plants. Besides, dispersed duplication and WGD/segmental duplication effectively promoted expansion of the GRF gene family. The expression patterns of GRF gene family and target genes were found in multiple floral organs and abundant in actively growing tissues. They were also found to be particularly expressed in response to various abiotic stresses, with stress-related elements in promoters, implying potential roles in floral development and abiotic stress. Our analysis in GRF gene family interaction network indicated the similar results that GRFs resist to abiotic stresses with the cooperation of other transcription factors like GIFs. This study provides insights into evolution in the GRF gene family, together with expression patterns valuable for future functional researches of plant abiotic stress biology.

TM3 and STM3 Promote Flowering Together with FUL2 and MBP20, but Act Antagonistically in Inflorescence Branching in Tomato

The moment at which a plant transitions to reproductive development is paramount to its life cycle and is strictly controlled by many genes. The transcription factor SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1) plays a central role in this process in Arabidopsis. However, the role of SOC1 in tomato (Solanum lycopersicum) has been sparsely studied. Here, we investigated the function of four tomato SOC1 homologs in the floral transition and inflorescence development. We thoroughly characterized the SOC1-like clade throughout the Solanaceae and selected four tomato homologs that are dynamically expressed upon the floral transition. We show that of these homologs, TOMATO MADS 3 (TM3) and SISTER OF TM3 (STM3) promote the primary and sympodial transition to flowering, while MADS-BOX PROTEIN 23 (MBP23) and MBP18 hardly contribute to flowering initiation in the indeterminate cultivar Moneyberg. Protein-protein interaction assays and whole-transcriptome analysis during reproductive meristem development revealed that TM3 and STM3 interact and share many targets with FRUITFULL (FUL) homologs, including cytokinin regulators. Furthermore, we observed that mutating TM3/STM3 affects inflorescence development, but counteracts the inflorescence-branching phenotype of ful2 mbp20. Collectively, this indicates that TM3/STM3 promote the floral transition together with FUL2/MBP20, while these transcription factors have opposite functions in inflorescence development.

Assessing the role of REM13, REM34 and REM46 during the transition to the reproductive phase in Arabidopsis thaliana

REM (reproductive meristem) transcription factors have been proposed as regulators of plant reproductive development mainly based on their specific expression patterns in reproductive structures, but their roles are still largely unknown probably because of their redundancy. We selected three REM genes (REM13, REM34 and REM46) for functional analysis, based on their genome position and/or co-expression data.Our results suggest that these genes have a role in flowering time regulation and may modulate cell cycle progression. In addition, protein interaction experiments revealed that REM34 and REM46 interact with each other, suggesting that they might work cooperatively to regulate cell division during inflorescence meristem commitment.Previous attempts of using co-expression data as a guide for functional analysis of REMs were limited by the transcriptomic data available at the time. Our results uncover previously unknown functions of three members of the REM family of Arabidopsis thaliana and open the door to more comprehensive studies of the REM family, where the combination of co-expression analysis followed by functional studies might contribute to uncovering the biological roles of these proteins and the relationship among them.

The nucleoporin NUP160 and NUP96 regulate nucleocytoplasmic export of mRNAs and participate in ethylene signaling and response in Arabidopsis

Arabidopsis nucleoporin involved in the regulation of ethylene signaling via controlling of nucleocytoplasmic transport of mRNAs. The two-way transport of mRNAs between the nucleus and cytoplasm are controlled by the nuclear pore complex (NPC). In higher plants, the NPC contains at least 30 nucleoporins. The Arabidopsis nucleoporins are involved in various biological processes such as pathogen interaction, nodulation, cold response, flowering, and hormone signaling. However, little is known about the regulatory functions of the nucleoporin NUP160 and NUP96 in ethylene signaling pathway. In the present study, we provided data showing that the Arabidopsis nucleoporin NUP160 and NUP96 participate in ethylene signaling-related mRNAs nucleocytoplasmic transport. The Arabidopsis nucleoporin mutants (nup160, nup96-1, nup96-2) exhibited enhanced ethylene sensitivity. Nuclear qRT-PCR analysis and poly(A)-mRNA in situ hybridization showed that the nucleoporin mutants affected the nucleocytoplasmic transport of all the examined mRNAs, including the ethylene signaling-related mRNAs such as ETR2, ERS1, ERS2, EIN4, CTR1, EIN2, and EIN3. Transcriptome analysis of the nucleoporin mutants provided clues suggesting that the nucleoporin NUP160 and NUP96 may participate in ethylene signaling via various molecular mechanisms. These observations significantly advance our understanding of the regulatory mechanisms of nucleoporin proteins in ethylene signaling and ethylene response.

Genomic targets of HOP2 are enriched for features found at recombination hotspots

HOP2 is a conserved protein that plays a positive role in homologous chromosome pairing and a separable role in preventing illegitimate connections between nonhomologous chromosome regions during meiosis. We employed ChIP-seq to discover that Arabidopsis HOP2 binds along the length of all chromosomes, except for centromeric and nucleolar organizer regions, and no binding sites were detected in the organelle genomes. A large number of reads were assigned to the HOP2 locus itself, yet TAIL-PCR and SNP analysis of the aligned sequences indicate that many of these reads originate from the transforming T-DNA, supporting the role of HOP2 in preventing nonhomologous exchanges. The 292 ChIP-seq peaks are largely found in promoter regions and downstream from genes, paralleling the distribution of recombination hotspots, and motif analysis revealed that there are several conserved sequences that are also enriched at crossover sites. We conducted coimmunoprecipitation of HOP2 followed by LC-MS/MS and found enrichment for several proteins, including some histone variants and modifications that are also known to be associated with recombination hotspots. We propose that HOP2 may be directed to chromatin motifs near double strand breaks, where homology checks are proposed to occur.

Anatomical observation and transcriptome analysis of buds reveal the association between the AP2 gene family and reproductive induction in hybrid larch (Larix kaempferi × Larix olgensis)

Hybrid larch is an excellent afforestation species in northern China. The instability of seed yield is an urgent problem to be solved. The biological characteristics related to seed setting in larch are different from those in angiosperms and other gymnosperms. Studying the developmental mechanism of the larch sporophyll can deepen our understanding of conifer reproductive development and help to ensure an adequate supply of seeds in the seed orchard. The results showed that the formation of microstrobilus primordia in hybrid larch could be observed in anatomical sections collected in the middle of July. The contents of endogenous gibberellin 3 (GA3) and abscisic acid (ABA) were higher and the contents of GA4, GA7, jasmonic acid and salicylic acid were lower in multiseeded larch. Transcriptome analysis showed that transcription factors were significantly enriched in the AP2 family. There were 23 differentially expressed genes in the buds of the multiseeded and less-seeded types, and the expression of most of these genes was higher in the buds than in the needles. We conclude that mid-July is the early stage of reproductive organ development in hybrid larch and is suitable for the study of reproductive development. GA3 and ABA may be helpful for improving seed setting in larch, and 23 AP2/EREBP family genes are involved in the regulation of reproductive development in larch.

Integrated analysis of miRNAome transcriptome and degradome reveals miRNA-target modules governing floral florescence development and senescence across early- and late-flowering genotypes in tree peony

As a candidate national flower of China, tree peony has extremely high ornamental, medicinal and oil value. However, the short florescence and rarity of early-flowering and late-flowering varieties restrict further improvement of the economic value of tree peony. Specific miRNAs and their target genes engaged in tree peony floral florescence, development and senescence remain unknown. This report presents the integrated analysis of the miRNAome, transcriptome and degradome of tree peony petals collected from blooming, initial flowering, full blooming and decay stages in early-flowering variety Paeonia ostii 'Fengdan', an early-flowering mutant line of Paeonia ostii 'Fengdan' and late-flowering variety Paeonia suffruticosa 'Lianhe'. Transcriptome analysis revealed a transcript ('psu.G.00014095') which was annotated as a xyloglucan endotransglycosylase/hydrolase precursor XTH-25 and found to be differentially expressed across flower developmental stages in Paeonia ostii 'Fengdan' and Paeonia suffruticosa 'Lianhe'. The miRNA-mRNA modules were presented significant enrichment in various pathways such as plant hormone signal transduction, indole alkaloid biosynthesis, arachidonic acid metabolism, folate biosynthesis, fatty acid elongation, and the MAPK signaling pathway. Multiple miRNA-mRNA-TF modules demonstrated the potential functions of MYB-related, bHLH, Trihelix, NAC, GRAS and HD-ZIP TF families in floral florescence, development, and senescence of tree peony. Comparative spatio-temporal expression investigation of eight floral-favored miRNA-target modules suggested that transcript 'psu.T.00024044' and microRNA mtr-miR166g-5p are involved in the floral florescence, development and senescence associated agronomic traits of tree peony. The results might accelerate the understanding of the potential regulation mechanism in regards to floral florescence, development and abscission, and supply guidance for tree peony breeding of varieties with later and longer florescence characteristics.

Identification of Brachypodium distachyon B3 genes reveals that BdB3-54 regulates primary root growth

B3 is a class of plant-specific transcription factors with important roles in plant development and growth. Here, we identified 69 B3 transcription factors in Brachypodium distachyon that were unevenly distributed across all five chromosomes. The ARF, REM, LAV, and RAV subfamilies were grouped based on sequence characteristics and phylogenetic relationships. The phylogenetically related members in the B3 family shared conserved domains and gene structures. Expression profiles showed that B3 genes were widely expressed in different tissues and varied in response to different abiotic stresses. BdB3-54 protein from the REM subfamily was located in the nucleus by subcellular localization and processed transcriptional activation activity. Overexpression of BdB3-54 in Arabidopsis increased primary root length. Our study provides a basis for further research on the functions of BdB3 genes.

Genome-wide identification and characterization profile of phosphatidy ethanolamine-binding protein family genes in carrot

Members of the family of Phosphatidy Ethanolamine-Binding Protein (PEBP) have been shown to be key regulators of the transition of plants from vegetative to reproductive phases. Here, a total of 12 PEBP proteins were identified in the carrot (Daucus carota L.) genome and classified into FT-like (4), TFL1-like (6), and MFT-like 2) subfamilies, that had different lengths (110-267 aa) and were distributed unevenly across seven chromosomes. Moreover, 13 and 31 PEBP proteins were identified in other two Apiaceae species, celery (Apium graveolens L.) and coriander (Coriandrum sativum L.). The phylogenetic and evolutionary results of these PEBP family proteins were obtained based on the protein sequences. In the three Apiaceae species, purifying selection was the main evolutionary force, and WGD, segmental duplication, and dispersed duplication have played key roles in the PEBP family expansion. The expression analysis showed that carrot PEBP genes exhibited relatively broad expression patterns across various tissues. In the period of bolting to flowering, the carrot FT-like subfamily genes were upregulated as positive regulators, and TFL1-like subfamily genes remained at lower expression levels as inhibitors. More interestingly, the members of carrot FT-like genes had different temporal-spatial expression characteristics, suggesting that they have different regulatory functions in the carrot reproductive phase. In summary, this study contributes to our understanding of the PEBP family proteins and provides a foundation for exploring the mechanism of carrot bolting and flowering for the breeding of cultivars with bolting resistance.

JAZ1-3 and MYC2-1 Synergistically Regulate the Transformation from Completely Mixed Flower Buds to Female Flower Buds in Castanea mollisima

Chestnut (Castanea mollisima) is an important woody food crop, but its yield has been low in cultivation, mainly due to the problems of fewer female flowers and more male flowers. Therefore, regulating the transition of chestnut flowers and effectively balancing the proportion of male and female to improve the yield are key factor to be solved in production. In this study, the chestnut floral buds in pre- and post-winter were used as materials. The data of metabolites, hormones, and gene expression during flower bud differentiation of chestnut were analyzed by transcriptomics and metabolomics to preliminarily reveal the possible reason of male and female flower bud transformation in pre- and post-winter. The analysis of Differentially Expressed Genes (DEGs) showed that there were 6323 DEGs in the Complete mixed flower bud (CMF) group in pre- and post-winter, of which 3448 genes were up-regulated and 2875 genes were down-regulated. There were 8037 DEGs in the Incomplete mixed flower bud (IMF) in pre- and post-winter, of which 4546 genes were up-regulated and 3491 genes were down-regulated. A total of 726 genes from the two flower buds were enriched into 251 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in post winter, of which plant hormone signal transduction accounted for 4.13%. The analysis results of differential metabolites showed that the differential metabolites of the two flower buds were mainly concentrated in the secondary metabolic synthesis pathway. The difference of hormone content showed that the content of Gibberellin 9 (GA9) and GA19 in CMF was higher than that in IMF in pre-winter, but the opposite in post-winter. Methyl jasmonate (MeJA) content was only very high in CMF in pre-winter, while Jasmonoyl-(l)-Isoleucine (JA-ILE) showed high content in CMF in post-winter. In post-winter, higher concentration of JA-ILE was positively correlated with the expression of Flowering Locus T (CmFT), and CmFT gene was significantly positively correlated with the expression levels of MYC2-1, MYC2-2 and LFY 3 (LEAFY 3). The higher concentration of JA-ILE was negatively correlated with the transcription level of JAZ1-3. In vitro experiments further verified that Jasmonate-Zim 1–3 (JAZ 1–3) combined with MYC2-1 inhibited the transcription of CmFT gene, while MYC2-1 alone promoted the expression of FT. The results suggested that a higher concentration of GA is conducive to breaking the dormancy of flower buds and promoting the development of male flower buds, while a lower concentration of GA and a higher concentration of JA-ILE are conducive to the differentiation and formation of female flower buds in post-winter, in which JAZ1-3 and MYC2-1 play a key role in the differentiation of female flower buds of chestnut.

Arabidopsis CPR5 plays a role in regulating nucleocytoplasmic transport of mRNAs in ethylene signaling pathway

Arabidopsis CPR5 is involved in regulation of ethylene signaling via two different ways: interacting with the ETR1 N-terminal domains, and controlling nucleocytoplasmic transport of ethylene-related mRNAs. The ETR1 receptor plays a predominant role in ethylene signaling in Arabidopsis thaliana. Previous studies showed that both RTE1 and CPR5 can directly bind to the ETR1 receptor and regulate ethylene signaling. RTE1 was suggested to promote the ETR1 receptor signaling by influencing its conformation, but little is known about the regulatory mechanism of CPR5 in ethylene signaling. In this study, we presented the data showing that both RTE1 and CPR5 bound to the N-terminal domains of ETR1, and regulated ethylene signaling via the ethylene receptor. On the other hand, the research provided evidence indicating that CPR5 could act as a nucleoporin to regulate the ethylene-related mRNAs export out of the nucleus, while RTE1 or its homolog (RTH) had no effect on the nucleocytoplasmic transport of mRNAs. Nuclear qRT-PCR analysis and poly(A)-mRNA in situ hybridization showed that defect of CPR5 restricted nucleocytoplasmic transport of mRNAs. These results advance our understanding of the regulatory mechanism of CPR5 in ethylene signaling.

The spinach YY genome reveals sex chromosome evolution, domestication, and introgression history of the species

Background

Spinach (Spinacia oleracea L.) is a dioecious species with an XY sex chromosome system, but its Y chromosome has not been fully characterized. Our knowledge about the history of its domestication and improvement remains limited.

Results

A high-quality YY genome of spinach is assembled into 952 Mb in six pseudo-chromosomes. By a combination of genetic mapping, Genome-Wide Association Studies, and genomic analysis, we characterize a 17.42-Mb sex determination region (SDR) on chromosome 1. The sex chromosomes of spinach evolved when an insertion containing sex determination genes occurred, followed by a large genomic inversion about 1.98 Mya. A subsequent burst of SDR-specific repeats (0.1–0.15 Mya) explains the large size of this SDR. We identify a Y-specific gene, NRT1/PTR 6.4 which resides in this insertion, as a strong candidate for the sex determination or differentiation factor. Resequencing of 112 spinach genomes reveals a severe domestication bottleneck approximately 10.87 Kya, which dates the domestication of spinach 7000 years earlier than the archeological record. We demonstrate that a strong selection signal associated with internode elongation and leaf area expansion is associated with domestication of edibility traits in spinach. We find that several strong genomic introgressions from the wild species Spinacia turkestanica and Spinacia tetrandra harbor desirable alleles of genes related to downy mildew resistance, frost resistance, leaf morphology, and flowering-time shift, which likely contribute to spinach improvement.

Conclusions

Analysis of the YY genome uncovers evolutionary forces shaping nascent sex chromosome evolution in spinach. Our findings provide novel insights about the domestication and improvement of spinach.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02633-x.

Identification of Flowering Regulatory Networks and Hub Genes Expressed in the Leaves of Elymus sibiricus L. Using Comparative Transcriptome Analysis

Flowering is a significant stage from vegetative growth to reproductive growth in higher plants, which impacts the biomass and seed yield. To reveal the flowering time variations and identify the flowering regulatory networks and hub genes in Elymus sibiricus, we measured the booting, heading, and flowering times of 66 E. sibiricus accessions. The booting, heading, and flowering times varied from 136 to 188, 142 to 194, and 148 to 201 days, respectively. The difference in flowering time between the earliest- and the last-flowering accessions was 53 days. Furthermore, transcriptome analyses were performed at the three developmental stages of six accessions with contrasting flowering times. A total of 3,526 differentially expressed genes (DEGs) were predicted and 72 candidate genes were identified, including transcription factors, known flowering genes, and plant hormone-related genes. Among them, four candidate genes (LATE, GA2OX6, FAR3, and MFT1) were significantly upregulated in late-flowering accessions. LIMYB, PEX19, GWD3, BOR7, PMEI28, LRR, and AIRP2 were identified as hub genes in the turquoise and blue modules which were related to the development time of flowering by weighted gene co-expression network analysis (WGCNA). A single-nucleotide polymorphism (SNP) of LIMYB found by multiple sequence alignment may cause late flowering. The expression pattern of flowering candidate genes was verified in eight flowering promoters (CRY, COL, FPF1, Hd3, GID1, FLK, VIN3, and FPA) and four flowering suppressors (CCA1, ELF3, Ghd7, and COL4) under drought and salt stress by qRT-PCR. The results suggested that drought and salt stress activated the flowering regulation pathways to some extent. The findings of the present study lay a foundation for the functional verification of flowering genes and breeding of new varieties of early- and late-flowering E. sibiricus.

GhADF6-mediated actin reorganization is associated with defence against Verticillium dahliae infection in cotton

Several studies have revealed that actin depolymerizing factors (ADFs) participate in plant defence responses; however, the functional mechanisms appear intricate and need further exploration. In this study, we identified an ADF6 gene in upland cotton (designated as GhADF6) that is evidently involved in cotton's response to the fungal pathogen Verticillium dahliae. GhADF6 binds to actin filaments and possesses actin severing and depolymerizing activities in vitro and in vivo. When cotton root (the site of the fungus invasion) was inoculated with the pathogen, the expression of GhADF6 was markedly down-regulated in the epidermal cells. By virus-induced gene silencing analysis, the down-regulation of GhADF6 expression rendered the cotton plants tolerant to V. dahliae infection. Accordingly, the abundance of actin filaments and bundles in the root cells was significantly higher than that in the control plant, which phenocopied that of the V. dahliae-challenged wild-type cotton plant. Altogether, our results provide evidence that an increase in filament actin (F-actin) abundance as well as dynamic actin remodelling are required for plant defence against the invading pathogen, which are likely to be fulfilled by the coordinated expressional regulation of the actin-binding proteins, including ADF.

Genome-wide analysis reveals widespread roles for RcREM genes in floral organ development in Rosa chinensis

Members of the REM (Reproductive Meristem) gene family are expressed primarily in reproductive meristems and floral organs. However, their evolution and their functional profiles in flower development remain poorly understood. Here, we performed genome-wide identification and evolutionary analysis of the REM gene family in Rosaceae. This family has been greatly expanded in rose (Rosa chinensis) compared to other species, primarily through tandem duplication. Expression analysis revealed that most RcREM genes are specifically expressed in reproductive organs and that their specific expression patterns are dramatically altered in rose plants with mutations affecting floral organs. Protein-protein interaction analysis indicated that RcREM14 interact with RcAP1 (one of the homology of A class genes in ABCDE model), highlighting the roles of RcREM genes in floral organ identity. Finally, co-expression network analysis indicated that RcREM genes are co-expressed with a high proportion of key genes that regulate flowering time, floral organ development, and cell proliferation and expansion in R. chinensis.

Transcriptome profiling of developing leaf and shoot apices to reveal the molecular mechanism and co-expression genes responsible for the wheat heading date

Background

Wheat is one of the most widely planted crops worldwide. The heading date is important for wheat environmental adaptability, as it not only controls flowering time but also determines the yield component in terms of grain number per spike.

Results

In this research, homozygous genotypes with early and late heading dates derived from backcrossed progeny were selected to conduct RNA-Seq analysis at the double ridge stage (W2.0) and androgynous primordium differentiation stage (W3.5) of the leaf and apical meristem, respectively. In total, 18,352 differentially expressed genes (DEGs) were identified, many of which are strongly associated with wheat heading date genes. Gene Ontology (GO) enrichment analysis revealed that carbohydrate metabolism, trehalose metabolic process, photosynthesis, and light reaction are closely related to the flowering time regulation pathway. Based on MapMan metabolic analysis, the DEGs are mainly involved in the light reaction, hormone signaling, lipid metabolism, secondary metabolism, and nucleotide synthesis. In addition, 1,225 DEGs were annotated to 45 transcription factor gene families, including LFY, SBP, and MADS-box transcription factors closely related to flowering time. Weighted gene co-expression network analysis (WGCNA) showed that 16, 336, 446, and 124 DEGs have biological connections with Vrn1-5 A, Vrn3-7B, Ppd-1D, and WSOC1, respectively. Furthermore, TraesCS2D02G181400 encodes a MADS-MIKC transcription factor and is co-expressed with Vrn1-5 A, which indicates that this gene may be related to flowering time.

Conclusions

RNA-Seq analysis provided transcriptome data for the wheat heading date at key flower development stages of double ridge (W2.0) and androgynous primordium differentiation (W3.5). Based on the DEGs identified, co-expression networks of key flowering time genes in Vrn1-5 A, Vrn3-7B, WSOC1, and Ppd-1D were established. Moreover, we discovered a potential candidate flowering time gene, TraesCS2D02G181400. Taken together, these results serve as a foundation for further study on the regulatory mechanism of the wheat heading date.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07797-7.

Development of EST-SSR Markers Linked to Flowering Candidate Genes in Elymus sibiricus L. Based on RNA Sequencing

Elymus sibiricus L. is an important cold-season grass with excellent cold and drought tolerance, good palatability, and nutrition. Flowering time is a key trait that affects forage and seed yield. Development of EST-SSR (expressed sequence tag simple sequence repeat) markers based on flowering genes contributes to the improvement of flowering traits. In the study, we detected 155 candidate genes related to flowering traits from 10,591 unigenes via transcriptome sequencing in early- and late-flowering genotypes. These candidate genes were mainly involved in the photoperiodic pathway, vernalization pathway, central integrator, and gibberellin pathway. A total of 125 candidate gene-based EST-SSRs were developed. Further, 15 polymorphic EST-SSRs closely associated to 13 candidate genes were used for genetic diversity and population structure analysis among 20 E. sibiricus accessions, including two contrasting panels (early-flowering and late-flowering). Among them, primer 28366, designed from heading date 3a (HD3a), effectively distinguished early- and late-flowering genotypes using a specifically amplified band of 175 bp. The polymorphic information content (PIC) value ranged from 0.12 to 0.48, with an average of 0.25. The unweighted pair group method analysis (UPGMA) cluster and structure analysis showed that the 20 E. sibiricus genotypes with similar flowering times tended to group together. These newly developed EST-SSR markers have the potential to be used for molecular markers assisted selection and germplasm evaluation of flowering traits in E. sibiricus.