Cytochrome P450 CYP78A9 is involved in Arabidopsis reproductive development

Identification of anther thermotolerance genes by the integration of linkage and association analysis in maize

Maize is one of the world's most important staple crops, yet its production is increasingly threatened by the rising frequency of high-temperature stress (HTS). To investigate the genetic basis of anther thermotolerance under field conditions, we performed linkage and association analysis to identify HTS response quantitative trait loci (QTL) using three recombinant inbred line (RIL) populations and an association panel containing 375 diverse maize inbred lines. These analyses resulted in the identification of 16 co-located large QTL intervals. Among the 37 candidate genes identified in these QTL intervals, five have rice or Arabidopsis homologs known to influence pollen and filament development. Notably, one of the candidate genes, ZmDUP707, has been subject to selection pressure during breeding. Its expression is suppressed by HTS, leading to pollen abortion and barren seeds. We also identified several additional candidate genes potentially underly QTL previously reported by other researchers. Taken together, our results provide a pool of valuable candidate genes that could be employed by future breeding programs aiming at enhancing maize HTS tolerance.

Parthenocarpy, a pollination-independent fruit set mechanism to ensure yield stability

Fruit development is essential for flowering plants' reproduction and a significant food source. Climate change threatens fruit yields due to its impact on pollination and fertilization processes, especially vulnerable to extreme temperatures, insufficient light, and pollinator decline. Parthenocarpy, the development of fruit without fertilization, offers a solution, ensuring yield stability in adverse conditions and enhancing fruit quality. Parthenocarpic fruits not only secure agricultural production but also exhibit improved texture, appearance, and shelf life, making them desirable for food processing and other applications. Recent research unveils the molecular mechanisms behind parthenocarpy, implicating transcription factors (TFs), noncoding RNAs, and phytohormones such as auxin, gibberellin (GA), and cytokinin (CK). Here we review recent findings, construct regulatory models, and identify areas for further research.

Pleiotropic phenotypic effects of the TaCYP78A family on multiple yield-related traits in wheat

Increasing crop yield depends on selecting and utilizing pleiotropic genes/alleles to improve multiple yield-related traits (YRTs) during crop breeding. However, synergistic improvement of YRTs is challenging due to the trade-offs between YRTs in breeding practices. Here, the favourable haplotypes of the TaCYP78A family are identified by analysing allelic variations in 1571 wheat accessions worldwide, demonstrating the selection and utilization of pleiotropic genes to improve yield and related traits during wheat breeding. The TaCYP78A family members, including TaCYP78A3, TaCYP78A5, TaCYP78A16, and TaCYP78A17, are organ size regulators expressed in multiple organs, and their allelic variations associated with various YRTs. However, due to the trade-offs between YRTs, knockdown or overexpression of TaCYP78A family members does not directly increase yield. Favourable haplotypes of the TaCYP78A family, namely A3/5/16/17Ap-Hap II, optimize the expression levels of TaCYP78A3/5/16/17-A across different wheat organs to overcome trade-offs and improve multiple YRTs. Different favourable haplotypes have both complementary and specific functions in improving YRTs, and their aggregation in cultivars under strong artificial selection greatly increase yield, even under various planting environments and densities. These findings provide new support and valuable genetic resources for molecular breeding of wheat and other crops in the era of Breeding 4.0.

Population genomics and epigenomics of Spirodela polyrhiza provide insights into the evolution of facultative asexuality

Many plants are facultatively asexual, balancing short-term benefits with long-term costs of asexuality. During range expansion, natural selection likely influences the genetic controls of asexuality in these organisms. However, evidence of natural selection driving asexuality is limited, and the evolutionary consequences of asexuality on the genomic and epigenomic diversity remain controversial. We analyzed population genomes and epigenomes of Spirodela polyrhiza, (L.) Schleid., a facultatively asexual plant that flowers rarely, revealing remarkably low genomic diversity and DNA methylation levels. Within species, demographic history and the frequency of asexual reproduction jointly determined intra-specific variations of genomic diversity and DNA methylation levels. Genome-wide scans revealed that genes associated with stress adaptations, flowering and embryogenesis were under positive selection. These data are consistent with the hypothesize that natural selection can shape the evolution of asexuality during habitat expansions, which alters genomic and epigenomic diversity levels.

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.

Comparative transcriptome analysis reveals candidate genes for cold stress response and early flowering in pineapple

Pineapple originates from tropical regions in South America and is therefore significantly impacted by cold stress. Periodic cold events in the equatorial regions where pineapple is grown may induce early flowering, also known as precocious flowering, resulting in monetary losses due to small fruit size and the need to make multiple passes for harvesting a single field. Currently, pineapple is one of the most important tropical fruits in the world in terms of consumption, and production losses caused by weather can have major impacts on worldwide exportation potential and economics. To further our understanding of and identify mechanisms for low-temperature tolerance in pineapple, and to identify the relationship between low-temperature stress and flowering time, we report here a transcriptomic analysis of two pineapple genotypes in response to low-temperature stress. Using meristem tissue collected from precocious flowering-susceptible MD2 and precocious flowering-tolerant Dole-17, we performed pairwise comparisons and weighted gene co-expression network analysis (WGCNA) to identify cold stress, genotype, and floral organ development-specific modules. Dole-17 had a greater increase in expression of genes that confer cold tolerance. The results suggested that low temperature stress in Dole-17 plants induces transcriptional changes to adapt and maintain homeostasis. Comparative transcriptomic analysis revealed differences in cuticular wax biosynthesis, carbohydrate accumulation, and vernalization-related gene expression between genotypes. Cold stress induced changes in ethylene and abscisic acid-mediated pathways differentially between genotypes, suggesting that MD2 may be more susceptible to hormone-mediated early flowering. The differentially expressed genes and module hub genes identified in this study are potential candidates for engineering cold tolerance in pineapple to develop new varieties capable of maintaining normal reproduction cycles under cold stress. In addition, a total of 461 core genes involved in the development of reproductive tissues in pineapple were also identified in this study. This research provides an important genomic resource for understanding molecular networks underlying cold stress response and how cold stress affects flowering time in pineapple.

Sweet cherry AP2/ERF transcription factor, PavRAV2, negatively modulates fruit size by directly repressing PavKLUH expression

For sweet cherry, fruit size is one of the main targets in breeding programs owing to the high market value of larger fruits. KLUH/CYP78A5 is an important regulator of seed/fruit size in several plant species, but its molecular mechanism is largely unknown. In this study, we characterized the function of PavKLUH in the regulation of sweet cherry fruit size. The ectopic overexpression of PavKLUH in Arabidopsis increased the size of its siliques and seeds, whereas virus-induced gene silencing of PavKLUH in sweet cherry significantly decreased fruit size by restricting mesocarp cell expansion. We screened out an AP2/ERF transcription factor containing a B3-like domain, designated as PavRAV2, which was able to physically interact with PavKLUH promoter in a yeast one-hybrid (Y1H) system. In Y1H assays, electrophoretic mobility shift assays, and dual-luciferase reporter analyses, PavRAV2 directly bound to the promoter of PavKLUH in vitro and in vivo, and suppressed PavKLUH expression. Silencing of PavRAV2 resulted in enlarged fruit as a result of enhanced mesocarp cell expansion. Together, our results provide new insights into signaling pathways related to fruit size, and outline a possible mechanism for how the RAV transcription factor directly regulates CYP78A family members to influence fruit size and development.

Cytochrome P450 Gene Families: Role in Plant Secondary Metabolites Production and Plant Defense

Cytochrome P450s (CYPs) are the most prominent family of enzymes involved in NADPH- and O2-dependent hydroxylation processes throughout all spheres of life. CYPs are crucial for the detoxification of xenobiotics in plants, insects, and other organisms. In addition to performing this function, CYPs serve as flexible catalysts and are essential for producing secondary metabolites, antioxidants, and phytohormones in higher plants. Numerous biotic and abiotic stresses frequently affect the growth and development of plants. They cause a dramatic decrease in crop yield and a deterioration in crop quality. Plants protect themselves against these stresses through different mechanisms, which are accomplished by the active participation of CYPs in several biosynthetic and detoxifying pathways. There are immense potentialities for using CYPs as a candidate for developing agricultural crop species resistant to biotic and abiotic stressors. This review provides an overview of the plant CYP families and their functions to plant secondary metabolite production and defense against different biotic and abiotic stresses.

Genomic survey of high-throughput RNA-Seq data implicates involvement of long intergenic non-coding RNAs (lincRNAs) in cytoplasmic male-sterility and fertility restoration in pigeon pea

BACKGROUND

Long-intergenic non-coding RNAs (lincRNAs) originate from intergenic regions and have no coding potential. LincRNAs have emerged as key players in the regulation of various biological processes in plant development. Cytoplasmic male-sterility (CMS) in association with restorer-of-fertility (Rf) systems makes it a highly reliable tool for exploring heterosis for producing commercial hybrid seeds. To date, there have been no reports of lincRNAs during pollen development in CMS and fertility restorer lines in pigeon pea.

OBJECTIVE

Identification of lincRNAs in the floral buds of cytoplasmic male-sterile (AKCMS11) and fertility restorer (AKPR303) pigeon pea lines.

METHODS

We employed a computational approach to identify lincRNAs in the floral buds of cytoplasmic male-sterile (AKCMS11) and fertility restorer (AKPR303) pigeon pea lines using RNA-Seq data.

RESULTS

We predicted a total of 2145 potential lincRNAs of which 966 were observed to be differentially expressed between the sterile and fertile pollen. We identified, 927 cis-regulated and 383 trans-regulated target genes of the lincRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of the target genes revealed that these genes were specifically enriched in pathways like pollen and pollen tube development, oxidative phosphorylation, etc. We detected 23 lincRNAs that were co-expressed with 17 pollen-related genes with known functions. Fifty-nine lincRNAs were predicted to be endogenous target mimics (eTMs) for 25 miRNAs, and found to be associated with pollen development. The, lincRNA regulatory networks revealed that different lincRNA-miRNA-mRNA networks might be associated with CMS and fertility restoration.

CONCLUSION

Thus, this study provides valuable information by highlighting the functions of lincRNAs as regulators during pollen development in pigeon pea and utilization in hybrid seed production.

Whole-Genome Comparison Reveals Structural Variations behind Heading Leaf Trait in Brassica oleracea

Brassica oleracea displays remarkable morphological variations. It intrigued researchers to study the underlying cause of the enormous diversification of this organism. However, genomic variations in complex heading traits are less known in B. oleracea. Herein, we performed a comparative population genomics analysis to explore structural variations (SVs) responsible for heading trait formation in B. oleracea. Synteny analysis showed that chromosomes C1 and C2 of B. oleracea (CC) shared strong collinearity with A01 and A02 of B. rapa (AA), respectively. Two historical events, whole genome triplication (WGT) of Brassica species and differentiation time between AA and CC genomes, were observed clearly by phylogenetic and Ks analysis. By comparing heading and non-heading populations of B. oleracea genomes, we found extensive SVs during the diversification of the B. oleracea genome. We identified 1205 SVs that have an impact on 545 genes and might be associated with the heading trait of cabbage. Overlapping the genes affected by SVs and the differentially expressed genes identified by RNA-seq analysis, we identified six vital candidate genes that may be related to heading trait formation in cabbage. Further, qRT-PCR experiments also verified that six genes were differentially expressed between heading leaves and non-heading leaves, respectively. Collectively, we used available genomes to conduct a comparison population genome analysis and identify candidate genes for the heading trait of cabbage, which provides insight into the underlying reason for heading trait formation in B. oleracea.

Construction of a high density genetic linkage map to define the locus conferring seedlessness from Mukaku Kishu mandarin

Mukaku Kishu ('MK'), a small sized mandarin, is an important source of seedlessness in citrus breeding. Identification and mapping the gene(s) governing 'MK' seedlessness will expedite seedless cultivar development. In this study, two 'MK'-derived mapping populations- LB8-9 Sugar Belle^® ('SB') × 'MK' (N=97) and Daisy ('D') × 'MK' (N=68) were genotyped using an Axiom_Citrus56 Array encompassing 58,433 SNP probe sets, and population specific male and female parent linkage maps were constructed. The parental maps of each population were integrated to produce sub-composite maps, which were further merged to develop a consensus linkage map. All the parental maps (except 'MK_D') had nine major linkage groups, and contained 930 ('SB'), 810 ('MK_SB'), 776 ('D') and 707 ('MK_D') SNPs. The linkage maps displayed 96.9 ('MK_D') to 98.5% ('SB') chromosomal synteny with the reference Clementine genome. The consensus map was comprised of 2588 markers including a phenotypic seedless (Fs)-locus and spanned a genetic distance of 1406.84 cM, with an average marker distance of 0.54 cM, which is substantially lower than the reference Clementine map. For the phenotypic Fs-locus, the distribution of seedy and seedless progenies in both 'SB' × 'MK' (55:42, χ² = 1.74) and 'D' × 'MK' populations (33:35, χ² = 0.06) followed a test cross pattern. The Fs-locus mapped on chromosome 5 with SNP marker 'AX-160417325' at 7.4 cM in 'MK_SB' map and between two SNP markers 'AX-160536283' and 'AX-160906995' at a distance of 2.4 and 4.9 cM, respectively in 'MK_D' map. The SNPs 'AX-160417325' and 'AX-160536283' correctly predicted seedlessness of 25-91.9% progenies in this study. Based on the alignment of flanking SNP markers to the Clementine reference genome, the candidate gene for seedlessness hovered in a ~ 6.0 Mb region between 3.97 Mb (AX-160906995) to 10.00 Mb (AX-160536283). This region has 131 genes of which 13 genes (belonging to seven gene families) reportedly express in seed coat or developing embryo. The findings of the study will prove helpful in directing future research for fine mapping this region and eventually underpinning the exact causative gene governing seedlessness in 'MK'.

BSR and Full-Length Transcriptome Approaches Identified Candidate Genes for High Seed Ratio in Camellia vietnamensis

(1) Background: C. vietnamensis is very suitable for growth in the low hilly areas of southern subtropical regions. Under appropriate conditions, the oil yield of C. vietnamensis can reach 1125 kg/ha (the existing varieties can reach 750 kg/ha). Moreover, the fruit of C. vietnamensis is large and the pericarp is thick (>5 cm). Therefore, a high seed ratio has become the main target economic trait for the breeding of C. vietnamensis. (2) Methods: A half-sibling population of C. vietnamensis plants with a combination of high and low seed ratios was constructed by crossing a C. vietnamensis female parent. Bulked segregant RNA analysis and full-length transcriptome sequencing were performed to determine the molecular mechanisms underlying a high seed ratio. (3) Results: Seed ratio is a complex quantitative trait with a normal distribution, which is significantly associated with four other traits of fruit (seed weight, seed number, fruit diameter, and pericarp thickness). Two candidate regions related to high seed ratio (HSR) were predicted. One spanned 140.8−148.4 Mb of chromosome 2 and was associated with 97 seed-yield-related candidate genes ranging in length from 278 to 16,628 bp. The other spanned 35.3−37.3 Mb on chromosome 15 and was associated with 38 genes ranging in length from 221 to 16,928 bp. Using the full-length transcript as a template, a total of 115 candidate transcripts were obtained, and 78 transcripts were predicted to be functionally annotated. The DEGs from two set pairs of cDNA sequencing bulks were enriched to cytochrome p450 CYP76F14 (KOG0156; GO:0055114, HSR4, HSR7), the gibberellin phytohormone pathway (GO:0016787, HSR5), the calcium signaling pathway (GO:0005509, HSR6), the polyubiquitin-PPAR signaling pathway (GO:0005515, HSR2, HSR3), and several main transcription factors (bZIP transcription factor, HSR1) in C. vietnamensis.

QTL Detection for Rice Grain Length and Fine Mapping of a Novel Locus qGL6.1

BACKGROUND

Grain length (GL) that is directly associated with appearance quality is a key target of selection in rice breeding. Although abundant quantitative trait locus (QTL) associated with GL have been identified, it was still relatively weak to identify QTL for GL from japonica genetic background, as the shortage of japonica germplasms with long grains. We performed QTLs analysis for GL using a recombinant inbred lines (RILs) population derived from the cross between japonica variety GY8 (short grains) and LX1 (long grains) in four environments.

RESULTS

A total of 197 RILs were genotyped with 285 polymorphic SNP markers. Three QTLs qGL5.3, qGL6.1 and qGL11 were detected to control GL by individual environmental analyses and multi-environment joint analysis. Of these, a major-effect and stable QTL qGL6.1 was identified to be a novel QTL, and its LX1 allele had a positive effect on GL. For fine-mapping qGL6.1, a BC₁F₂ population consisting of 2,487 individuals was developed from a backcross between GY8 and R176, one line with long grain. Eight key informative recombinants were identified by nine kompetitive allele specific PCR (KASP) markers. By analyzing key recombinants, the qGL6.1 locus was narrowed down to a 40.41 kb genomic interval on chromosome 6. One candidate gene LOC_Os06g43304.1 encoding cytochrome P450 (CYP71D55) was finally selected based on the difference in the transcriptional expression and variations in its upstream and downstream region.

CONCLUSIONS

Three QTLs qGL5.3, qGL6.1 and qGL11 were identified to control grain length in rice. One novel QTL qGL6.1 was fine mapped within 40.41 kb region, and LOC_Os06g43304.1 encoding cytochrome P450 (CYP71D55) may be its candidate gene. We propose that the further cloning of the qGL6.1 will facilitate improving appearance quality in japonica varieties.

Overexpression of Rice BSR2 Confers Disease Resistance and Induces Enlarged Flowers in Torenia fournieri Lind

Plant pathogens evade basal defense systems and attack different organs and tissues of plants. Genetic engineering of plants with genes that confer resistance against pathogens is very effective in pathogen control. Conventional breeding for disease resistance in ornamental crops is difficult and lagging relative to that in non-ornamental crops due to an inadequate number of disease-resistant genes. Therefore, genetic engineering of these plants with defense-conferring genes is a practical approach. We used rice BSR2 encoding CYP78A15 for developing transgenic Torenia fournieri Lind. lines. The overexpression of BSR2 conferred resistance against two devastating fungal pathogens, Rhizoctonia solani and Botrytis cinerea. In addition, BSR2 overexpression resulted in enlarged flowers with enlarged floral organs. Histological observation of the petal cells suggested that the enlargement in the floral organs could be due to the elongation and expansion of the cells. Therefore, the overexpression of BSR2 confers broad-spectrum disease resistance and induces the production of enlarged flowers simultaneously. Therefore, this could be an effective strategy for developing ornamental crops that are disease-resistant and economically more valuable.

TaKLU Plays as a Time Regulator of Leaf Growth via Auxin Signaling

The growth of leaves is subject to strict time regulation. Several genes influencing leaf growth have been identified, but little is known about how genes regulate the orderly initiation and growth of leaves. Here, we demonstrate that TaKLU/TaCYP78A5 contributes to a time regulation mechanism in leaves from initiation to expansion. TaKLU encodes the cytochrome P450 CYP78A5, and its homolog AtKLU has been described whose deletion is detrimental to organ growth. Our results show that TaKLU overexpression increases leaf size and biomass by altering the time of leaf initiation and expansion. TaKLU-overexpressing plants have larger leaves with more cells. Further dynamic observations indicate that enlarged wheat leaves have experienced a longer expansion time. Different from AtKLU inactivation increases leaf number and initiation rates, TaKLU overexpression only smooths the fluctuations of leaf initiation rates by adjusting the initiation time of local leaves, without affecting the overall leaf number and initiation rates. In addition, complementary analyses suggest TaKLU is functionally conserved with AtKLU in controlling the leaf initiation and size and may involve auxin accumulation. Our results provide a new insight into the time regulation mechanisms of leaf growth in wheat.

Genome- and transcriptome-wide association studies reveal the genetic basis and the breeding history of seed glucosinolate content in Brassica napus

A high content of seed glucosinolates and their degradation products imposes anti-nutritional effects on livestock; therefore, persistent efforts are made to reduce the seed GSL content to increase the commercial value of rapeseed meal. Here, we dissected the genetic structure of SGC by genome-wide association studies (GWAS) combined with transcriptome-wide association studies (TWAS). Fifteen reliable quantitative trait loci (QTLs) were identified to be associated with the reduced SGC in modern B. napus cultivars by GWAS. Analysis of the selection strength and haplotypes at these QTLs revealed that low SGC was predominantly generated by the co-selection of qGSL.A02.2, qGSL.C02.1, qGSL.A09.2, and qGSL.C09.1. Integration of the results from TWAS, comprehensive bioinformatics, and POCKET algorithm analyses indicated that BnaC02.GTR2 (BnaC02g42260D) is a candidate gene underlying qGSL.C02.1. Using CRISPR/Cas9-derived Bna.gtr2s knockout mutants, we experimentally verified that both BnaC02.GTR2 and its three paralogs positively regulate seed GSL accumulation but negatively regulated vegetative tissue GSL contents. In addition, we observed smaller seeds with higher seed oil content in these Bna.gtr2 mutants. Furthermore, both RNA-seq and correlation analyses suggested that Bna.GTR2s might play a comprehensive role in seed development, such as amino acid accumulation, GSL synthesis, sugar assimilation, and oil accumulation. This study unravels the breeding selection history of low-SGC improvement and provides new insights into the molecular function of Bna.GTR2s in both seed GSL accumulation and seed development in B. napus.

Specialized metabolism and development: An unexpected friendship

Plants produce a myriad of metabolites. Some of them have been regarded for a long time as secondary or specialized metabolites and are considered to have functions mostly in defense and the adaptation of plants to their environment. However, in the last years, new research has shown that these metabolites can also have roles in the regulation of plant growth and development, some acting as signals, through the interaction with hormonal pathways, and some independently of them. These reports provide a glimpse of the functional possibilities that specialized metabolites present in the modulation of plant development and encourage more research in this direction.

Strigolactones (SLs) modulate the plastochron by regulating KLUH (KLU) transcript abundance in Arabidopsis

The timing of leaf emergence at the shoot apical meristem, or plastochron, is highly regulated in plants. Among the genes known to regulate the plastochron in Arabidopsis (Arabidopsis thaliana), KLUH (KLU), orthologous to the rice (Oryza sativa) PLASTOCHRON1, encodes the cytochrome P450 CYP78A5, and is thought to act through generation of a still unknown mobile signal. As klu mutants display not only a short plastochron but also a branching phenotype reminiscent of strigolactone (SL) mutants, we investigated whether KLU/CYP78A5 is involved in SL biosynthesis. We combined a genetic approach, a parasitic plant seed germination bioassay to test klu root exudates, and analysis of transcript abundances of SL-biosynthesis genes in the Arabidopsis klu mutants. We demonstrate that KLU is not involved in the SL-biosynthesis pathway. Moreover, this work allowed us to uncover a new role for SL during Arabidopsis development in modulating plastochron via a KLU-dependent pathway. Globally our data reveal that KLU is required for plastochron-specific SL responses, a first indication of crosstalk between SL and the KLU-derived signal.

QTL Analysis of Five Silique-Related Traits in Brassica napus L. Across Multiple Environments

As an important physiological and reproductive organ, the silique is a determining factor of seed yield and a breeding target trait in rapeseed (Brassica napus L.). Genetic studies of silique-related traits are helpful for rapeseed marker-assisted high-yield breeding. In this study, a recombinant inbred population containing 189 lines was used to perform a quantitative trait loci (QTLs) analysis for five silique-related traits in seven different environments. As a result, 120 consensus QTLs related to five silique-related traits were identified, including 23 for silique length, 25 for silique breadth, 29 for silique thickness, 22 for seed number per silique and 21 for silique volume, which covered all the chromosomes, except C5. Among them, 13 consensus QTLs, one, five, two, four and one for silique length, silique breadth, silique thickness, seed number per silique and silique volume, respectively, were repeatedly detected in multiple environments and explained 4.38-13.0% of the phenotypic variation. On the basis of the functional annotations of Arabidopsis homologous genes and previously reported silique-related genes, 12 potential candidate genes underlying these 13 QTLs were screened and found to be stable in multiple environments by analyzing the re-sequencing results of the two parental lines. These findings provide new insights into the gene networks affecting silique-related traits at the QTL level in rapeseed.

Resequencing of 388 cassava accessions identifies valuable loci and selection for variation in heterozygosity

Background

Heterozygous genomes are widespread in outcrossing and clonally propagated crops. However, the variation in heterozygosity underlying key agronomic traits and crop domestication remains largely unknown. Cassava is a staple crop in Africa and other tropical regions and has a highly heterozygous genome.

Results

We describe a genomic variation map from 388 resequenced genomes of cassava cultivars and wild accessions. We identify 52 loci for 23 agronomic traits through a genome-wide association study. Eighteen allelic variations in heterozygosity for nine candidate genes are significantly associated with seven key agronomic traits. We detect 81 selective sweeps with decreasing heterozygosity and nucleotide diversity, harboring 548 genes, which are enriched in multiple biological processes including growth, development, hormone metabolisms and responses, and immune-related processes. Artificial selection for decreased heterozygosity has contributed to the domestication of the large starchy storage root of cassava. Selection for homozygous GG allele in MeTIR1 during domestication contributes to increased starch content. Selection of homozygous AA allele in MeAHL17 is associated with increased storage root weight and cassava bacterial blight (CBB) susceptibility. We have verified the positive roles of MeTIR1 in increasing starch content and MeAHL17 in resistance to CBB by transient overexpression and silencing analysis. The allelic combinations in MeTIR1 and MeAHL17 may result in high starch content and resistance to CBB.

Conclusions

This study provides insights into allelic variation in heterozygosity associated with key agronomic traits and cassava domestication. It also offers valuable resources for the improvement of cassava and other highly heterozygous crops.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02524-7.

Genome-Wide Identification and Analyses of Drought/Salt-Responsive Cytochrome P450 Genes in Medicago truncatula

Cytochrome P450 monooxygenases (P450s) catalyze a great number of biochemical reactions and play vital roles in plant growth, development and secondary metabolism. As yet, the genome-scale investigation on P450s is still lacking in the model legume Medicago truncatula. In particular, whether and how many MtP450s are involved in drought and salt stresses for Medicago growth, development and yield remain unclear. In this study, a total of 346 MtP450 genes were identified and classified into 10 clans containing 48 families. Among them, sixty-one MtP450 genes pairs are tandem duplication events and 10 MtP450 genes are segmental duplication events. MtP450 genes within one family exhibit high conservation and specificity in intron–exon structure. Meanwhile, many Mt450 genes displayed tissue-specific expression pattern in various tissues. Specifically, the expression pattern of 204 Mt450 genes under drought/NaCl treatments were analyzed by using the weighted correlation network analysis (WGCNA). Among them, eight genes (CYP72A59v1, CYP74B4, CYP71AU56, CYP81E9, CYP71A31, CYP704G6, CYP76Y14, and CYP78A126), and six genes (CYP83D3, CYP76F70, CYP72A66, CYP76E1, CYP74C12, and CYP94A52) were found to be hub genes under drought/NaCl treatments, respectively. The expression levels of these selected hub genes could be induced, respectively, by drought/NaCl treatments, as validated by qPCR analyses, and most of these genes are involved in the secondary metabolism and fatty acid pathways. The genome-wide identification and co-expression analyses of M. truncatulaP450 superfamily genes established a gene atlas for a deep and systematic investigation of P450 genes in M. truncatula, and the selected drought-/salt-responsive genes could be utilized for further functional characterization and molecular breeding for resistance in legume crops.

Downstream of GA4, PbCYP78A6 participates in regulating cell cycle-related genes and parthenogenesis in pear (Pyrus bretshneideri Rehd.)

BACKGROUND

Parthenocarpy results in traits attractive to both consumers and breeders, and it overcomes the obstacle of self-incompatibility in the fruit set of horticultural crops, including pear (Pyrus bretshneider). However, there is limited knowledge regarding the genetic and molecular mechanisms that regulate parthenogenesis.

RESULTS

Here, in a transcriptional comparison between pollination-dependent fruit and GA₄-induced parthenocarpy, PbCYP78A6 was identified and proposed as a candidate gene involved in parthenocarpy. PbCYP78A6 is similar to Arabidopsis thaliana CYP78A6 and highly expressed in pear hypanthia. The increased PbCYP78A6 expression, as assessed by RT-qPCR, was induced by pollination and GA₄ exposure. The ectopic overexpression of PbCYP78A6 contributed to parthenocarpic fruit production in tomato. The PbCYP78A6 expression coincided with fertilized and parthenocarpic fruitlets development and the expression of fruit development-related genes as assessed by cytological observations and RT-qPCR, respectively. PbCYP78A6 RNA interference and overexpression in pear calli revealed that the gene is an upstream regulator of specific fruit development-related genes in pear.

CONCLUSIONS

Our findings indicate that PbCYP78A6 plays a critical role in fruit formation and provide insights into controlling parthenocarpy.

A global survey of the gene network and key genes for oil accumulation in cultivated tetraploid cottons

To enrich our knowledge about gene network of fatty acid biosynthesis in cottonseed, we conducted comparative transcriptome to reveal the differences in gene expression between Gossypium hirsutum and Gossypium barbadense during cottonseed development. The prolonged expression period and increased expression abundance of oil-related genes are the main reasons for producing high seed oil content (SOC) in G. barbadense, which manifested as the bias of homeologous gene expression in Dt-subgenome after 25 day postanthesis (DPA). The dynamic expression profile showed that SAD6 and FATA are more important for oil biosynthesis in G. barbadense than that in G. hirsutum. Three key transcription factors, WRI1, NF-YB6 and DPBF2, showed their elite roles in regulating seed oil in cotton. We observed that sequence variations in the promoter region of BCCP2 genes might contribute to its divergence in expression level between the two species. Based on the quantitative trait loci (QTL) information of the seed oil content and utilizing additional G. barbadense introgression lines (ILs), we propose 21 candidate genes on the basis of their differential expression level, of which the GbSWEET and the GbACBP6 showed the potential functional to improve the oil content. Taken together, studying the different expression of oil-related genes and their genetic regulation mechanisms between G. hirsutum and G. barbadense provide new insights to understanding the mechanism of fatty acid biosynthesis network and fatty acid genetic improving breeding in cotton.

Highly pleiotropic functions of CYP78As and AMP1 are regulated in non-cell-autonomous/organ-specific manners

The cytochrome P450 CYP78A5/KLUH in Arabidopsis thaliana is predicted to be involved in the synthesis of a mobile signal molecule that has a pleiotropic function that is distinct from classical phytohormones. CYP78A5 has five close relatives in Arabidopsis. We first investigated their functions, focusing on the plastochron, leaf size, and leaf senescence. Our analyses revealed that CYP78A5 and CYP78A7 are involved in the plastochron and leaf size, and CYP78A6 and CYP78A9 are involved in leaf senescence. Complementation analyses using heterologous promoters and expression analyses suggested that CYP78A isoforms have a common biochemical function and are functionally differentiated via organ-specific expression. The altered meristem program1 (amp1) carboxypeptidase mutant shows a phenotype very similar to that of the cyp78a5 mutant. Complementation analyses using boundary and organizing center-specific promoters suggested that both CYP78A5 and AMP1 act in a non-cell-autonomous manner. Analyses of multiple cyp78a mutants and crosses between cyp78a and amp1 mutants revealed that AMP1/LIKE AMP1 (LAMP1) and CYP78A isoforms regulate plastochron length and leaf senescence in the same genetic pathway, whereas leaf size is independently regulated. Furthermore, we detected feedback regulation between CYP78A6/CYP78A9 and AMP1 at the gene expression level. These observations raise the possibility that AMP1 and CYP78A isoforms are involved in the synthesis of the same mobile signal molecule, and suggest that AMP1 and CYP78A signaling pathways have a very close, albeit complex, functional relationship.

Shy Girl, a kiwifruit suppressor of feminization, restricts gynoecium development via regulation of cytokinin metabolism and signalling

Kiwifruit (Actinidia chinensis) is a dioecious, long-living woody perennial vine. Reduced generation time and induction of hermaphroditism can accelerate crop improvement and facilitate alternative farming for better food security in the face of climate change. Previous studies identified that CENTRORADIALIS genes CEN and CEN4 act to repress flowering, whilst the male-specific Shy Girl (SyGl) gene with homology to type-C cytokinin response regulators could repress gynoecium development in model plants. Here we use CRISPR/Cas9 to mutagenize CEN, CEN4 and SyGl in the male kiwifruit A. chinensis 'Bruce'. Biallelic mutations of CEN and CEN4 generated rapid-flowering male plants, and simultaneous targeting of CEN4 and SyGl gave rise to rapid-flowering hermaphrodites with restored gynoecial function and viable pollen, providing functional evidence for the role of SyGl in suppression of feminization. Analysis of ovary tissues identified genes that contribute to carpel development and revealed that SyGl affected both cytokinin profiles and the expression of genes involved in cytokinin metabolism and signalling. The plant lines generated by CEN4/SyGl knockout could self-pollinate and produce fast-flowering offspring. These results establish that SyGI acts as the suppressor of feminization in kiwifruit and demonstrate the potential for accelerated breeding in an outcrossing horticultural woody perennial.

Tomato agamous-like6 parthenocarpy is facilitated by ovule integument reprogramming involving the growth regulator KLUH

Fruit set is established during and soon after fertilization of the ovules inside the quiescent ovary, but the signaling pathways involved remain obscure. The tomato (Solanum lycopersicum) CRISPR loss-of-function mutant of the transcription factor gene AGAMOUS-like6 (SlAGL6; slagl6CR-sg1) is capable of fertilization-independent setting of normal, yet seedless (parthenocarpic), fruit. To gain insight into the mechanism of fleshy fruit set, in this study, we investigated how slagl6CR-sg1 uncouples fruit set from fertilization. We found that mutant ovules were enlarged due to integument over-proliferation and failed to differentiate an endothelium, the integument's innermost layer, upon maturation. A causal relationship between slagl6 loss-of-function and these abnormal phenotypes is inferred from the observation that SlAGL6 is predominantly expressed in the immature ovule integument, and upon ovule maturation, its expression shifts to the endothelium. The transcriptome of unfertilized mutant ovules profoundly differs from that of wild-type and exhibits substantial overlap with the transcriptomes of fertilized ovules sporophytic tissues. One prominent upregulated gene was the fertilization-induced cytochrome P450 cell proliferation regulator SlKLUH. Indeed, ectopic overexpression of SlKLUH stimulated both integument growth in unfertilized ovules and parthenocarpy, suggesting that its suppression by SlAGL6 is paramount for preventing fertilization-independent fruit set. Taken together, our study informs on the transcriptional programs that are regulated by SlAGL6 and demonstrates that it acts from within the ovule integument to inhibit ovary growth beyond anthesis. That by suppressing components of the fertilization-induced ovule reprogramming underlying fruit set.

Somatic variants for seed and fruit set in grapevine

BACKGROUND

Grapevine reproductive development has direct implications on yield. It also impacts on berry and wine quality by affecting traits like seedlessness, berry and bunch size, cluster compactness and berry skin to pulp ratio. Seasonal fluctuations in yield, fruit composition and wine attributes, which are largely driven by climatic factors, are major challenges for worldwide table grape and wine industry. Accordingly, a better understanding of reproductive processes such as gamete development, fertilization, seed and fruit set is of paramount relevance for managing yield and quality. With the aim of providing new insights into this field, we searched for clones with contrasting seed content in two germplasm collections.

RESULTS

We identified eight variant pairs that seemingly differ only in seed-related characteristics while showing identical genotype when tested with the GrapeReSeq_Illumina_20K_SNP_chip and several microsatellites. We performed multi-year observations on seed and fruit set deriving from different pollination treatments, with special emphasis on the pair composed by Sangiovese and its seedless variant locally named Corinto Nero. The pollen of Corinto Nero failed to germinate in vitro and gave poor berry set when used to pollinate other varieties. Most berries from both open- and cross-pollinated Corinto Nero inflorescences did not contain seeds. The genetic analysis of seedlings derived from occasional Corinto Nero normal seeds revealed that the few Corinto Nero functional gametes are mostly unreduced. Moreover, three genotypes, including Sangiovese and Corinto Nero, were unexpectedly found to develop fruits without pollen contribution and occasionally showed normal-like seeds. Five missense single nucleotide polymorphisms were identified between Corinto Nero and Sangiovese from transcriptomic data.

CONCLUSIONS

Our observations allowed us to attribute a seedlessness type to some variants for which it was not documented in the literature. Interestingly, the VvAGL11 mutation responsible for Sultanina stenospermocarpy was also discovered in a seedless mutant of Gouais Blanc. We suggest that Corinto Nero parthenocarpy is driven by pollen and/or embryo sac defects, and both events likely arise from meiotic anomalies. The single nucleotide polymorphisms identified between Sangiovese and Corinto Nero are suitable for testing as traceability markers for propagated material and as functional candidates for the seedless phenotype.

Targeted mutagenesis of EOD3 gene in Brassica napus L. regulates seed production

Seed size and number are central to the evolutionary fitness of plants and are also crucial for seed production of crops. However, the molecular mechanisms of seed production control are poorly understood in Brassica crops. Here, we report the gene cloning, expression analysis, and functional characterization of the EOD3/CYP78A6 gene in rapeseed. BnaEOD3 has four copies located in two subgenomes, which exhibited a steady higher expression during seed development with differential expression among copies. The targeted mutations of BnaEOD3 gene were efficiently generated by stable transformation of the CRISPR/Cas9 (clustered regularly interspaced short palindromic repeat) vector. These mutations were stably transmitted to T1 and T2 generations and a large collection of homozygous mutants with combined loss-of-function alleles across four BnaEOD3 copies were created for phenotyping. All mutant T1 lines had shorter siliques, smaller seeds, and an increased number of seeds per silique, in which the quadrable mutants showed the most significant changes in these traits. Consequently, the seed weight per plant in the quadrable mutants increased by 13.9% on average compared with that of wild type, indicating that these BnaEOD3 copies have redundant functions in seed development in rapeseed. The phenotypes of the different allelic combinations of BnaEOD3 copies also revealed gene functional differentiation among the two subgenomes. Cytological observations indicated that the BnaEOD3 could act maternally to promote cotyledon cell expansion and proliferation to regulate seed growth in rapeseed. Collectively, our findings reveal the quantitative involvement of the different BnaEOD3 copies function in seed development, but also provided valuable resources for rapeseed breeding programs.

Differential expression of SlKLUH controlling fruit and seed weight is associated with changes in lipid metabolism and photosynthesis-related genes

The sizes of plant organs such as fruit and seed are crucial yield components. Tomato KLUH underlies the locus fw3.2, an important regulator of fruit and seed weight. However, the mechanism by which the expression levels of KLUH affect organ size is poorly understood. We found that higher expression of SlKLUH increased cell proliferation in the pericarp within 5 d post-anthesis in tomato near-isogenic lines. Differential gene expression analyses showed that lower expression of SlKLUH was associated with increased expression of genes involved in lipid metabolism. Lipidomic analysis revealed that repression of SlKLUH mainly increased the contents of certain non-phosphorus glycerolipids and phospholipids and decreased the contents of four unknown lipids. Co-expression network analyses revealed that lipid metabolism was possibly associated with but not directly controlled by SlKLUH, and that this gene instead controls photosynthesis-related processes. In addition, many transcription factors putatively involved in the KLUH pathway were identified. Collectively, we show that SlKLUH regulates fruit and seed weight which is associated with altered lipid metabolism. The results expand our understanding of fruit and seed weight regulation and offer a valuable resource for functional studies of candidate genes putatively involved in regulation of organ size in tomato and other crops.

Enhanced resistance to fungal and bacterial diseases in tomato and Arabidopsis expressing BSR2 from rice

The overexpression of rice BSR2 would offer a simple and effective strategy to protect plants from multiple devastating diseases in tomato and Arabidopsis. Many devastating plant diseases are caused by pathogens possessing a wide host range. Fungal Botrytis cinerea and Rhizoctonia solani, as well as bacterial Pseudomonas syringae and Ralstonia pseudosolanacearum are four such pathogens that infect hundreds of plant species, including agronomically important crops, and cause serious diseases, leading to severe economic losses. However, reports of genes that can confer resistance to broad host-range pathogens via traditional breeding methods are currently limited. We previously reported that Arabidopsis plants overexpressing rice BROAD-SPECTRUM RESISTANCE2 (BSR2/CYP78A15) showed tolerance not only to bacterial P. syringae pv. tomato DC3000 but also to fungal Colletotrichum higginsianum and R. solani. Rice plants overexpressing BSR2 displayed tolerance to two R. solani anastomosis groups. In the present study, first, BSR2-overexpressing (OX) Arabidopsis plants were shown to be additionally tolerant to B. cinerea, R. solani, and R. pseudosolanacearum. Next, tomato 'Micro-Tom' was used as a model to determine whether such tolerance by BSR2 can be introduced into dicot crops to prevent infection from pathogens possessing wide host range. BSR2-OX tomato displayed broad-spectrum disease tolerance to fungal B. cinerea and R. solani, as well as to bacterial P. syringae and R. pseudosolanacearum. Additionally, undesirable traits such as morphological changes were not detected. Thus, BSR2 overexpression can offer a simple and effective strategy to protect crops from multiple destructive diseases.

Genome-wide analysis of PmTCP4 transcription factor binding sites by ChIP-Seq during pistil abortion in Japanese apricot

The TCP4 transcription factor plays an important role in plant growth and development, especially in flower development. PmTCP4 is involved in the process of pistil abortion in Japanese apricot, but its molecular mechanism, particularly the DNA binding sites and co-regulatory genes, are quite unknown. Therefore, to identify the genome-wide binding sites of PmTCP4 transcription factors and their co-regulatory genes, chromatin immunoprecipitation sequencing (ChIP-Seq) was carried out. ChIP-Seq data produced the maximum enriched peaks in two Japanese apricot cultivars 'Daqiandi' (DQD) and 'Longyan' (LY), which showed that the majority of DNA-protein interactions are relevant and have a significant function in binding sites. Moreover, 720 and 251 peak-associated genes regulated by PmTCP4 were identified in DQD and LY, respectively, and most of them were involved in the flower and pistil development process. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that photosynthesis and oxidative phosphorylation were the most enriched pathways in both cultivars and all identified genes related to these pathways were down-regulated. This study will provide a reference for a better understanding of the PmTCP4 regulatory mechanism during pistil abortion in Japanese apricot.

VviERF6Ls: an expanded clade in Vitis responds transcriptionally to abiotic and biotic stresses and berry development

Background

VviERF6Ls are an uncharacterized gene clade in Vitis with only distant Arabidopsis orthologs. Preliminary data indicated these transcription factors may play a role in berry development and extreme abiotic stress responses. To better understand this highly duplicated, conserved clade, additional members of the clade were identified in four Vitis genotypes. A meta-data analysis was performed on publicly available microarray and RNA-Seq data (confirmed and expanded with RT-qPCR), and Vitis VviERF6L1 overexpression lines were established and characterized with phenotyping and RNA-Seq.

Results

A total of 18 PN40024 VviERF6Ls were identified; additional VviERF6Ls were identified in Cabernet Sauvignon, Chardonnay, and Carménère. The amino acid sequences of VviERF6Ls were found to be highly conserved. VviERF6L transcripts were detected in numerous plant organs and were differentially expressed in response to numerous abiotic stresses including water deficit, salinity, and cold as well as biotic stresses such as red blotch virus, N. parvum, and E. necator. VviERF6Ls were differentially expressed across stages of berry development, peaking in the pre-veraison/veraison stage and retaining conserved expression patterns across different vineyards, years, and Vitis cultivars. Co-expression network analysis identified a scarecrow-like transcription factor and a calmodulin-like gene with highly similar expression profiles to the VviERF6L clade. Overexpression of VviERF6L1 in a Seyval Blanc background did not result in detectable morphological phenotypes. Genes differentially expressed in response to VviERF6L1 overexpression were associated with abiotic and biotic stress responses.

Conclusions

VviERF6Ls represent a large and distinct clade of ERF transcription factors in grapevine. The high conservation of protein sequence between these 18 transcription factors may indicate these genes originate from a duplication event in Vitis. Despite high sequence similarity and similar expression patterns, VviERF6Ls demonstrate unique levels of expression supported by similar but heterogeneous promoter sequences. VviERF6L gene expression differed between Vitis species, cultivars and organs including roots, leaves and berries. These genes respond to berry development and abiotic and biotic stresses. VviERF6L1 overexpression in Vitis vinifera results in differential expression of genes related to phytohormone and immune system signaling. Further investigation of this interesting gene family is warranted.

Generation of High Yielding and Fragrant Rice (Oryza sativa L.) Lines by CRISPR/Cas9 Targeted Mutagenesis of Three Homoeologs of Cytochrome P450 Gene Family and OsBADH2 and Transcriptome and Proteome Profiling of Revealed Changes Triggered by Mutations

The significant increase in grain yield and quality are often antagonistic but a constant demand for breeders and consumers. Some genes related to cytochrome P450 family are known for rice organ growth but their role in controlling grain yield is still unknown. Here, we generated new rice mutants with high yield and improved aroma by simultaneously editing three cytochrome P450 homoeologs (Os03g0603100, Os03g0568400, and GL3.2) and OsBADH2 with the CRISPR/Cas9 system, and RNA-sequencing and proteomic analysis were performed to unveil the subsequent changes. High mutation efficiency was achieved in both target sites of each gene and the mutations were predominantly only deletions, while insertions were rare, and no mutations were detected in the five most likely off-target sites against each sgRNA. Mutants exhibited increased grain size, 2-acetyl-1-pyrroline (2AP) content, and grain cell numbers while there was no change in other agronomic traits. Transgene-DNA-free mutant lines appeared with a frequency of 44.44% and homozygous mutations were stably transmitted, and bi-allelic and heterozygous mutations followed Mendelian inheritance, while the inheritance of chimeric mutations was unpredictable. Deep RNA sequencing and proteomic results revealed the regulation of genes and proteins related to cytochrome P450 family, grain size and development, and cell cycle. The KEGG and hub-gene and protein network analysis showed that the gene and proteins related to ribosomal and photosynthesis pathways were mainly enriched, respectively. Our findings provide a broad and detailed basis to understand the role of CRISPR/Cas9 in rice yield and quality improvement.

Genetic and signalling pathways of dry fruit size: targets for genome editing-based crop improvement

Fruit is seed-bearing structures specific to angiosperm that form from the gynoecium after flowering. Fruit size is an important fitness character for plant evolution and an agronomical trait for crop domestication/improvement. Despite the functional and economic importance of fruit size, the underlying genes and mechanisms are poorly understood, especially for dry fruit types. Improving our understanding of the genomic basis for fruit size opens the potential to apply gene-editing technology such as CRISPR/Cas to modulate fruit size in a range of species. This review examines the genes involved in the regulation of fruit size and identifies their genetic/signalling pathways, including the phytohormones, transcription and elongation factors, ubiquitin-proteasome and microRNA pathways, G-protein and receptor kinases signalling, arabinogalactan and RNA-binding proteins. Interestingly, different plant taxa have conserved functions for various fruit size regulators, suggesting that common genome edits across species may have similar outcomes. Many fruit size regulators identified to date are pleiotropic and affect other organs such as seeds, flowers and leaves, indicating a coordinated regulation. The relationships between fruit size and fruit number/seed number per fruit/seed size, as well as future research questions, are also discussed.

Comparative functional genomics analysis of cytochrome P450 gene superfamily in wheat and maize

BACKGROUND

The cytochrome P450s (CYP450s) as the largest enzyme family of plant metabolism participate in various physiological processes, whereas no study has demonstrated interest in comprehensive comparison of the genes in wheat and maize. Genome-wide survey, characterization and comparison of wheat and maize CYP450 gene superfamily are useful for genetic manipulation of the Gramineae crops.

RESULTS

In total, 1285 and 263 full-length CYP450s were identified in wheat and maize, respectively. According to standard nomenclature, wheat CYP450s (TaCYP450s) were categorized into 45 families, while maize CYP450s (ZmCYP450s) into 43 families. A comprehensive analysis of wheat and maize CYP450s, involved in functional domains, conserved motifs, phylogeny, gene structures, chromosome locations and duplicated events was performed. The result showed that each family/subfamily in both species exhibited characteristic features, suggesting their phylogenetic relationship and the potential divergence in their functions. Functional divergence analysis at the amino acid level of representative clans CYP51, CYP74 and CYP97 in wheat, maize and rice identified some critical amino acid sites that are responsible for functional divergence of a gene family. Expression profiles of Ta-, ZmCYP450s were investigated using RNA-seq data, which contribute to infer the potential functions of the genes during development and stress responses. We found in both species CYP450s had preferential expression in specific tissues, and many tissue-specific genes were identified. Under water-deficit condition, 82 and 39 significantly differentially expressed CYP450s were respectively detected in wheat and maize. These genes may have some roles in protecting plants against drought damage. Thereinto, fourteen CYP450s were selected to validate their expression level through qRT-PCR. To further elucidating molecular mechanisms of CYP450 action, gene co-expression network was constructed. In total, 477 TaCYP450s were distributed in 22 co-expression modules, and some co-expressed genes that likely take part in the same biochemical pathway were identified. For instance, the expression of TaCYP74A98_4D was highly correlated with TaLOX9, TaLOX36, TaLOX39, TaLOX44 and TaOPR8, and all of them may be involved in jasmonate (JA) biosynthesis. TaCYP73A201_3A showed coexpression with TaPAL1.25, TaCCoAOMT1.2, TaCOMT.1, TaCCR1.6 and TaLAC5, which probably act in the wheat stem and/or root lignin synthesis pathway.

CONCLUSION

Our study first established systematic information about evolutionary relationship, expression pattern and function characterization of CYP450s in wheat and maize.

Dominant inhibition of awn development by a putative zinc-finger transcriptional repressor expressed at the B1 locus in wheat

Awns, bristle-like structures extending from grass lemmas, provide protection against predators, contribute to photosynthesis and aid in grain dispersal. In wheat, selection of awns with minimal extension, termed awnletted, has occurred during domestication by way of loci that dominantly inhibit awn development, such as Tipped1 (B1), Tipped2 (B2), and Hooded (Hd). Here we identify and characterize the B1 gene. B1 was identified using bulked segregant RNA-sequencing of an F₂ durum wheat population and through deletion mapping of awned bread wheat mutants. Functional characterization was accomplished by gene overexpression while haplotype analyses assessed B1 polymorphisms and genetic variation. Located on chromosome 5A, B1 is a C2H2 zinc finger encoding gene with ethylene-responsive element binding factor-associated amphiphilic repression (EAR) motifs. Constitutive overexpression of B1 in awned wheat produced an awnletted phenotype with pleiotropic effects on plant height and fertility. Transcriptome analysis of B1 overexpression plants suggests a role as transcriptional repressor, putatively targeting pathways involved in cell proliferation. Haplotype analysis revealed a conserved B1 coding region with proximal polymorphisms and supported the contention that B1 is mainly responsible for awnletted wheats globally. B1, predominantly responsible for awn inhibition in wheat, encodes a C2H2 zinc finger protein with EAR motifs which putatively functions as a transcriptional repressor.

Growth dynamics of the Arabidopsis fruit is mediated by cell expansion

Fruit have evolved a sophisticated tissue and cellular architecture to secure plant reproductive success. Postfertilization growth is perhaps the most dramatic event during fruit morphogenesis. Several studies have proposed that fertilized ovules and developing seeds initiate signaling cascades to coordinate and promote the growth of the accompanying fruit tissues. This dynamic process allows the fruit to conspicuously increase its size and acquire its final shape and means for seed dispersal. All these features are key for plant survival and crop yield. Despite its importance, we lack a high-resolution spatiotemporal map of how postfertilization fruit growth proceeds at the cellular level. In this study, we have combined live imaging, mutant backgrounds in which fertilization can be controlled, and computational modeling to monitor and predict postfertilization fruit growth in Arabidopsis We have uncovered that, unlike leaves, sepals, or roots, fruit do not exhibit a spatial separation of cell division and expansion domains; instead, there is a separation into temporal stages with fertilization as the trigger for transitioning to cell expansion, which drives postfertilization fruit growth. We quantified the coordination between fertilization and fruit growth by imaging no transmitting tract (ntt) mutants, in which fertilization fails in the bottom half of the fruit. By combining our experimental data with computational modeling, we delineated the mobility properties of the seed-derived signaling cascades promoting growth in the fruit. Our study provides the basis for generating a comprehensive understanding of the molecular and cellular mechanisms governing fruit growth and shape.

Fine mapping of a silique length- and seed weight-related gene in Brassica napus

Using microarray analysis combined with map-based cloning, a major locus positively regulating SL and SW was mapped to a 98.47 kb interval on A09 in rapeseed. In rapeseed, seed yield is closely associated with silique-related traits such as silique length (SL) and seed weight (SW). Previously identified quantitative trait loci (QTLs) revealed that SL and SW are complex traits and many QTLs overlap. However, the genetic characterization of the association between SL and SW is poorly understood. In the present study, a BC₃F₃ near isogenic line developed from a short silique plant and the long silique cultivar 'ZS11' was analyzed to identify the locus related to SL. Map-based cloning indicated that a major locus acting as a single Mendelian factor was mapped to a 98.47 kb region on chromosome A09. BLAST analysis and DNA sequencing showed SNP variations and a fragment replacement in the upstream region of the candidate gene BnaA09g55530D may alter gene expression and influence SL. The results showed that this SL locus may also positively affect SW as well as in the 186 rapeseed accessions identified by the associated markers. Therefore, selecting plants with appropriate SL and developing functional markers for the associated gene could play important roles in the molecular breeding of high-yield rapeseed varieties.

Natural Variation in Ovule Morphology Is Influenced by Multiple Tissues and Impacts Downstream Grain Development in Barley (Hordeum vulgare L.)

The ovule plays a critical role in cereal yield as it is the site of fertilization and the progenitor of the grain. The ovule primordium is generally comprised of three domains, the funiculus, chalaza, and nucellus, which give rise to distinct tissues including the integuments, nucellar projection, and embryo sac. The size and arrangement of these domains varies significantly between model eudicots, such as Arabidopsis thaliana, and agriculturally important monocotyledonous cereal species, such as Hordeum vulgare (barley). However, the amount of variation in ovule development among genotypes of a single species, and its functional significance, remains unclear. To address this, wholemount clearing was used to examine the details of ovule development in barley. Nine sporophytic and gametophytic features were examined at ovule maturity in a panel of 150 European two-row spring barley genotypes, and compared with grain traits from the preceding and same generation. Correlations were identified between ovule traits and features of grain they produced, which in general highlighted a negative correlation between nucellus area, ovule area, and grain weight. We speculate that the amount of ovule tissue, particularly the size of the nucellus, may affect the timing of maternal resource allocation to the fertilized embryo sac, thereby influencing subsequent grain development.

A CACTA-like transposable element in the upstream region of BnaA9.CYP78A9 acts as an enhancer to increase silique length and seed weight in rapeseed

Rapeseed (Brassica napus L.) is a model plant for polyploid crop research and the second-leading source of vegetable oil worldwide. Silique length (SL) and seed weight are two important yield-influencing traits in rapeseed. Using map-based cloning, we isolated qSLWA9, which encodes a P450 monooxygenase (BnaA9.CYP78A9) and functions as a positive regulator of SL. The expression level of BnaA9.CYP78A9 in silique valves of the long-silique variety is much higher than that in the regular-silique variety, which results in elongated cells and a prolonged phase of silique elongation. Plants of the long-silique variety and transgenic plants with high expression of BnaA9.CYP78A9 had a higher concentration of auxin in the developing silique; this induced a number of auxin-related genes but no genes in well-known auxin biosynthesis pathways, suggesting that BnaA9.CYP78A9 may influence auxin concentration by affecting auxin metabolism or an unknown auxin biosynthesis pathway. A 3.7-kb CACTA-like transposable element (TE) inserted in the 3.9-kb upstream regulatory sequence of BnaA9.CYP78A9 elevates the expression level, suggesting that the CACTA-like TE acts as an enhancer to stimulate high gene expression and silique elongation. Marker and sequence analysis revealed that the TE in B. napus had recently been introgressed from Brassica rapa by interspecific hybridization. The insertion of the TE is consistently associated with long siliques and large seeds in both B. napus and B. rapa collections. However, the frequency of the CACTA-like TE in rapeseed varieties is still very low, suggesting that this allele has not been widely used in rapeseed breeding programs and would be invaluable for yield improvement in rapeseed breeding.

Maternal control of seed weight in rapeseed (Brassica napus L.): the causal link between the size of pod (mother, source) and seed (offspring, sink)

Seed size/weight is one of the key traits related to plant domestication and crop improvement. In rapeseed (Brassica napus L.) germplasm, seed weight shows extensive variation, but its regulatory mechanism is poorly understood. To identify the key mechanism of seed weight regulation, a systematic comparative study was performed. Genetic, morphological and cytological evidence showed that seed weight was controlled by maternal genotype, through the regulation of seed size mainly via cell number. The physiological evidence indicated that differences in the pod length might result in differences in pod wall photosynthetic area, carbohydrates and the final seed weight. We also identified two pleiotropic major quantitative trait loci that acted indirectly on seed weight via their effects on pod length. RNA-seq results showed that genes related to pod development and hormones were significantly differentially expressed in the pod wall; genes related to development, cell division, nutrient reservoir and ribosomal proteins were all up-regulated in the seeds of the large-seed pool. Finally, we proposed a potential seed weight regulatory mechanism that is specific to rapeseed and novel in plants. The results demonstrate a causal link between the size of the pod (mother, source) and the seed (offspring, sink) in rapeseed, which provides novel insight into the maternal control of seed weight and will open a new research field in plants.

The rice CYP78A gene BSR2 confers resistance to Rhizoctonia solani and affects seed size and growth in Arabidopsis and rice

The fungal pathogen Rhizoctonia solani causes devastating diseases in hundreds of plant species. Among these, R. solani causes sheath blight, one of the three major diseases in rice. To date, few genes have been reported that confer resistance to R. solani. Here, rice-FOX Arabidopsis lines identified as having resistance to a bacterial pathogen, Pseudomonas syringae pv. tomato DC3000, and a fungal pathogen, Colletotrichum higginsianum were screened for disease resistance to R. solani. BROAD-SPECTRUM RESISTANCE2 (BSR2), a gene encoding an uncharacterized cytochrome P450 protein belonging to the CYP78A family, conferred resistance to R. solani in Arabidopsis. When overexpressed in rice, BSR2 also conferred resistance to two R. solani anastomosis groups. Both Arabidopsis and rice plants overexpressing BSR2 had slower growth and produced longer seeds than wild-type control plants. In contrast, BSR2-knockdown rice plants were more susceptible to R. solani and displayed faster growth and shorter seeds in comparison with the control. These results indicate that BSR2 is associated with disease resistance, growth rate and seed size in rice and suggest that its function is evolutionarily conserved in both monocot rice and dicot Arabidopsis.

New roles of NO TRANSMITTING TRACT and SEEDSTICK during medial domain development in Arabidopsis fruits

The gynoecium, the female reproductive part of the flower, is key for plant sexual reproduction. During its development, inner tissues such as the septum and the transmitting tract tissue, important for pollen germination and guidance, are formed. In Arabidopsis, several transcription factors are known to be involved in the development of these tissues. One of them is NO TRANSMITTING TRACT (NTT), essential for transmitting tract formation. We found that the NTT protein can interact with several gynoecium-related transcription factors, including several MADS-box proteins, such as SEEDSTICK (STK), known to specify ovule identity. Evidence suggests that NTT and STK control enzyme and transporter-encoding genes involved in cell wall polysaccharide and lipid distribution in gynoecial medial domain cells. The results indicate that the simultaneous loss of NTT and STK activity affects polysaccharide and lipid deposition and septum fusion, and delays entry of septum cells to their normal degradation program. Furthermore, we identified KAWAK, a direct target of NTT and STK, which is required for the correct formation of fruits in Arabidopsis These findings position NTT and STK as important factors in determining reproductive competence.

Transcriptomic response of Arabidopsis thaliana roots to naproxen and praziquantel

Exposition to pharmaceutical compounds released to the environment is considered as a potential risk for various organisms. We exposed Arabidopsis thaliana plants to naproxen (NAP) and praziquantel (PZQ) in 5 µM concentration for 2 days and recorded transcriptomic response in their roots with the aim to estimate ecotoxicity and to identify gene candidates potentially involved in metabolism of both compounds. Nonsteroidal anti-inflammatory drug NAP up-regulated 105 and down-regulated 29 genes (p-value ≤ 0.1, fold change ≥ 2), while anthelmintic PZQ up-regulated 389 and down-regulated 353 genes with more rigorous p-value ≤ 0.001 (fold change ≥ 2). High number of up-regulated genes coding for heat shock proteins and other genes involved in response to biotic and abiotic stresses as well as down-regulation of genes involved in processes such as cell proliferation, transcription and water transport indicates serious negative effect of PZQ. NAP up-regulated mostly genes involved in various biological processes and signal transduction and down-regulated mainly genes involved in signal transduction and electron transport or energy pathways. Further, two cytochrome P450s (demethylation) and one methyltransferase (methylation of carboxyl group) were identified as candidates for phase I and several glutathione- and glycosyltransferases (conjugation) for phase II of NAP metabolism. Cytochrome P450s, glutathione and glycosyltransferases seem to play role also in metabolism of PZQ. Up-regulation of several ABC and MATE transporters by NAP and PZQ indicated their role in transport of both compounds.

Analysis of ethylene-induced gene regulation during carposporogenesis in the red seaweed Grateloupia imbricata (Rhodophyta)

Ethylene favors carposporogenesis in the red seaweed Grateloupia imbricata. Analyses of cystocarp development in vitro in thalli treated with ethylene suggest an interconnection between polyamine and ethylene biosynthesis pathways. Yet, little is known about molecular mechanisms underlying carposporogenesis. Here, we used droplet digital PCR to analyze genes encoding enzymes related to polyamine (Spermidine [Spd] synthase) and ethylene (ACC synthase) synthesis; a pivotal compound of both pathways (S-adenosyl methionine synthase, SAMS); the gene that encodes amine oxidase, which is involved in polyamine degradation, and a candidate gene involved in seaweed reproduction (ornithine decarboxylase, ODC). In addition, we analyzed genes encoding proteins related to stress and reactive oxygen species, ascorbate peroxidase (APX), cytochrome P450 and WD 40. We characterized gene expression in fertilized and fertile thalli from G. imbricata that were exposed to ethylene for 15 min at two time points after treatment (1 and 7 d). The differential gene expression of SAMS, Spd synthase, ACC synthase, and cytochrome P450 was related to disclosure and development of cystocarps in fertilized thalli that transitioned from having no visible cystocarps at 1 d to developing cystocarps at 7 d. Likewise, cytochrome P450 was associated with cystocarp disclosure and maturation. In addition, amine oxidase and APX were involved in fine-tuning polyamine and reactive oxygen species during carposporogenesis, respectively, whereas WD 40 did so in relation to ethylene signaling. Expression of the candidate gene ODC was increased when cystocarps were not visible (fertilized thalli, 1d), as previously described. This analysis suggests developmental stage-specific roles for these genes during carposporogenesis.

Arabidopsis thaliana ambient temperature responsive lncRNAs

BACKGROUND

Long non-coding RNAs (lncRNAs) have emerged as new class of regulatory molecules in animals where they regulate gene expression at transcriptional and post-transcriptional level. Recent studies also identified lncRNAs in plant genomes, revealing a new level of transcriptional complexity in plants. Thousands of lncRNAs have been predicted in the Arabidopsis thaliana genome, but only a few have been studied in depth.

RESULTS

Here we report the identification of Arabidopsis lncRNAs that are expressed during the vegetative stage of development in either the shoot apical meristem or in leaves. We found that hundreds of lncRNAs are expressed in these tissues, of which 50 show differential expression upon an increase in ambient temperature. One of these lncRNAs, FLINC, is down-regulated at higher ambient temperature and affects ambient temperature-mediated flowering in Arabidopsis.

CONCLUSION

A number of ambient temperature responsive lncRNAs were identified with potential roles in the regulation of temperature-dependent developmental changes, such as the transition from the vegetative to the reproductive (flowering) phase. The challenge for the future is to characterize the biological function and molecular mode of action of the large number of ambient temperature-regulated lncRNAs that have been identified in this study.

Comparative transcriptome analysis reveals the regulatory networks of cytokinin in promoting the floral feminization in the oil plant Sapium sebiferum

BACKGROUND

Sapium sebiferum, whose seeds contain high level of fatty acids, has been considered as one of the most important oil plants. However, the high male to female flower ratio limited the seed yield improvement and its industrial potentials. Thus, the study of the sex determination in S. sebiferum is of significant importance in increasing the seed yield.

RESULTS

In this study, we demonstrated that in S. sebiferum, cytokinin (CK) had strong feminization effects on the floral development. Exogenous application with 6-benzylaminopurine (6-BA) or thidiazuron (TDZ) significantly induced the development of female flowers and increased the fruit number. Interestingly, the feminization effects of cytokinin were also detected on the androecious genotype of S. sebiferum which only produce male flowers. To further investigate the mechanism underlying the role of cytokinin in the flower development and sex differentiation, we performed the comparative transcriptome analysis of the floral buds of the androecious plants subjected to 6-BA. The results showed that there were separately 129, 352 and 642 genes differentially expressed at 6 h, 12 h and 24 h after 6-BA treatment. Functional analysis of the differentially expressed genes (DEGs) showed that many genes are related to the hormonal biosynthesis and signaling, nutrients translocation and cell cycle. Moreover, there were twenty one flowering-related genes identified to be differentially regulated by 6-BA treatment. Specifically, the gynoecium development-related genes SPATULA (SPT), KANADI 2 (KAN2), JAGGED (JAG) and Cytochrome P450 78A9 (CYP79A9) were significantly up-regulated, whereas the expression of PISTILLATA (PI), TATA Box Associated Factor II 59 (TAFII59) and MYB Domain Protein 108 (MYB108) that were important for male organ development was down-regulated in response to 6-BA treatment, demonstrating that cytokinin could directly target the floral organ identity genes to regulate the flower sex.

CONCLUSIONS

Our work demonstrated that cytokinin is a potential regulator in female flower development in S. sebiferum. The transcriptome analysis of the floral sex transition from androecious to monoecious in response to cytokinin treatment on the androecious S. sebiferum provided valuable information related to the mechanism of sex determination in the perennial woody plants.

Genome-Wide Analysis of DNA Methylation During Ovule Development of Female-Sterile Rice fsv1

The regulation of female fertility is an important field of rice sexual reproduction research. DNA methylation is an essential epigenetic modification that dynamically regulates gene expression during development processes. However, few reports have described the methylation profiles of female-sterile rice during ovule development. In this study, ovules were continuously acquired from the beginning of megaspore mother cell meiosis until the mature female gametophyte formation period, and global DNA methylation patterns were compared in the ovules of a high-frequency female-sterile line (fsv1) and a wild-type rice line (Gui99) using whole-genome bisulfite sequencing (WGBS). Profiling of the global DNA methylation revealed hypo-methylation, and 3471 significantly differentially methylated regions (DMRs) were observed in fsv1 ovules compared with Gui99. Based on functional annotation and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis of differentially methylated genes (DMGs), we observed more DMGs enriched in cellular component, reproduction regulation, metabolic pathway, and other pathways. In particular, many ovule development genes and plant hormone-related genes showed significantly different methylation patterns in the two rice lines, and these differences may provide important clues for revealing the mechanism of female gametophyte abortion.

PaCYP78A9, a Cytochrome P450, Regulates Fruit Size in Sweet Cherry (Prunus avium L.)

Sweet cherry (Prunus avium L.) is an important fruit crop in which fruit size is strongly associated with commercial value; few genes associated with fruit size have, however, been identified in sweet cherry. Members of the CYP78A subfamily, a group of important cytochrome P450s, have been found to be involved in controlling seed size and development in Arabidopsis thaliana, rice, soybean, and tomato. However, the influence of CYP78A members in controlling organ size and the underlying molecular mechanisms in sweet cherry and other fruit trees remains unclear. Here, we characterized a P. avium CYP78A gene PaCYP78A9 that is thought to be involved in the regulation of fruit size and organ development using overexpression and silencing approaches. PaCYP78A9 was significantly expressed in the flowers and fruit of sweet cherry. RNAi silencing of PaCYP78A9 produced small cherry fruits and PaCYP78A9 was found to affect fruit size by mediating mesocarp cell proliferation and expansion during fruit growth and development. Overexpression of PaCYP78A9 in Arabidopsis resulted in increased silique and seed size and PaCYP78A9 was found to be highly expressed in the inflorescences and siliques of transgenic plants. Genes related to cell cycling and proliferation were downregulated in fruit from sweet cherry TRV::PaCYP78A9-silencing lines, suggesting that PaCYP78A9 is likely to be an important upstream regulator of cell cycle processes. Together, our findings indicate that PaCYP78A9 plays an essential role in the regulation of cherry fruit size and provide insights into the molecular basis of the mechanisms regulating traits such as fruit size in P. avium.

Altered expression of the bZIP transcription factor DRINK ME affects growth and reproductive development in Arabidopsis thaliana

Here we describe an uncharacterized gene that negatively influences Arabidopsis growth and reproductive development. DRINK ME (DKM; bZIP30) is a member of the bZIP transcription factor family, and is expressed in meristematic tissues such as the inflorescence meristem (IM), floral meristem (FM), and carpel margin meristem (CMM). Altered DKM expression affects meristematic tissues and reproductive organ development, including the gynoecium, which is the female reproductive structure and is determinant for fertility and sexual reproduction. A microarray analysis indicates that DKM overexpression affects the expression of cell cycle, cell wall, organ initiation, cell elongation, hormone homeostasis, and meristem activity genes. Furthermore, DKM can interact in yeast and in planta with proteins involved in shoot apical meristem maintenance such as WUSCHEL, KNAT1/BP, KNAT2 and JAIBA, and with proteins involved in medial tissue development in the gynoecium such as HECATE, BELL1 and NGATHA1. Taken together, our results highlight the relevance of DKM as a negative modulator of Arabidopsis growth and reproductive development.