Spatiotemporal profiles of gene activity in stamen delineate nucleo-cytoplasmic interaction in a male-sterile somatic cybrid citrus.

Cytoplasmic male sterility (CMS) has long been used to produce seedless fruits in perennial woody crops like citrus. A male-sterile somatic cybrid citrus (G1 + HBP) was generated by protoplast fusion between a CMS callus parent 'Guoqing No. 1' Satsuma mandarin (Citrus unshiu, G1) and a fertile mesophyll parent Hirado Buntan pummelo (Citrus grandis, HBP). To uncover the male-sterile mechanism of G1 + HBP, we compared the transcriptome profiles of stamen organ and cell types at five stages between G1 + HBP and HBP, including the initial stamen primordia, enlarged stamen primordia, pollen mother cells, tetrads, and microspores captured by laser microdissection. The stamen organ and cell types showed distinct gene expression profiles. A majority of genes involved in stamen development were differentially expressed, especially CgAP3.2, which was downregulated in enlarged stamen primordia and upregulated in tetrads of G1 + HBP compared with HBP. Jasmonic acid- and auxin-related biological processes were enriched among the differentially expressed genes of stamen primordia, and the content of jasmonic acid biosynthesis metabolites was higher in flower buds and anthers of G1 + HBP. In contrast, the content of auxin biosynthesis metabolites was lower in G1 + HBP. The mitochondrial tricarboxylic acid cycle and oxidative phosphorylation processes were enriched among the differentially expressed genes in stamen primordia, meiocytes, and microspores, indicating the dysfunction of mitochondria in stamen organ and cell types of G1 + HBP. Taken together, the results indicate that malfunction of mitochondria-nuclear interaction might cause disorder in stamen development, and thus lead to male sterility in the citrus cybrid.

miR171-targeted SCARECROW-LIKE genes CsSCL2 and CsSCL3 regulate somatic embryogenesis in citrus.

Somatic embryogenesis (SE) is a key regeneration pathway in various biotechnology approaches to crop improvement, especially for economically important perennial woody crops like citrus. However, maintenance of SE capability has long been a challenge and becomes a bottleneck in biotechnology-facilitated plant improvement. In the embryogenic callus (EC) of citrus, we identified 2 csi-miR171c-targeted SCARECROW-LIKE genes CsSCL2 and CsSCL3 (CsSCL2/3), which exert positive feedback regulation on csi-miR171c expression. Suppression of CsSCL2 expression by RNA interference (RNAi) enhanced SE in citrus callus. A thioredoxin superfamily protein CsClot was identified as an interactive protein of CsSCL2/3. Overexpression of CsClot disturbed reactive oxygen species (ROS) homeostasis in EC and enhanced SE. Chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA-Seq identified 660 genes directly suppressed by CsSCL2 that were enriched in biological processes including development-related processes, auxin signaling pathway, and cell wall organization. CsSCL2/3 bound to the promoters of regeneration-related genes, such as WUSCHEL-RELATED HOMEOBOX 2 (CsWOX2), CsWOX13, and Lateral Organ Boundaries Domain 40 (LBD40), and repressed their expression. Overall, CsSCL2/3 modulate ROS homeostasis through the interactive protein CsClot and directly suppress the expression of regeneration-related genes, thus regulating SE in citrus. We uncovered a regulatory pathway of miR171c-targeted CsSCL2/3 in SE, which shed light on the mechanism of SE and regeneration capability maintenance in citrus.

MADS-box protein PpDAM6 regulates chilling requirement-mediated dormancy and bud break in peach.

Bud dormancy is crucial for winter survival and is characterized by the inability of the bud meristem to respond to growth-promotive signals before the chilling requirement (CR) is met. However, our understanding of the genetic mechanism regulating CR and bud dormancy remains limited. This study identified PpDAM6 (DORMANCY-ASSOCIATED MADS-box) as a key gene for CR using a genome-wide association study analysis based on structural variations in 345 peach (Prunus persica (L.) Batsch) accessions. The function of PpDAM6 in CR regulation was demonstrated by transiently silencing the gene in peach buds and stably overexpressing the gene in transgenic apple (Malus × domestica) plants. The results showed an evolutionarily conserved function of PpDAM6 in regulating bud dormancy release, followed by vegetative growth and flowering, in peach and apple. The 30-bp deletion in the PpDAM6 promoter was substantially associated with reducing PpDAM6 expression in low-CR accessions. A PCR marker based on the 30-bp indel was developed to distinguish peach plants with non-low and low CR. Modification of the H3K27me3 marker at the PpDAM6 locus showed no apparent change across the dormancy process in low- and non-low- CR cultivars. Additionally, H3K27me3 modification occurred earlier in low-CR cultivars on a genome-wide scale. PpDAM6 could mediate cell-cell communication by inducing the expression of the downstream genes PpNCED1 (9-cis-epoxycarotenoid dioxygenase 1), encoding a key enzyme for ABA biosynthesis, and CALS (CALLOSE SYNTHASE), encoding callose synthase. We shed light on a gene regulatory network formed by PpDAM6-containing complexes that mediate CR underlying dormancy and bud break in peach. A better understanding of the genetic basis for natural variations of CR can help breeders develop cultivars with different CR for growing in different geographical regions.

Genome-wide identification of the C2H2-Zinc finger gene family and functional validation of CsZFP7 in citrus nucellar embryogenesis.

KEY MESSAGE: Genome-wide identification of C2H2-ZF gene family in the poly- and mono-embryonic citrus species and validation of the positive role of CsZFP7 in sporophytic apomixis. The C2H2 zinc finger (C2H2-ZF) gene family is involved in plant vegetative and reproductive development. Although a large number of C2H2 zinc-finger proteins (C2H2-ZFPs) have been well characterized in some horticultural plants, little is known about the C2H2-ZFPs and their function in citrus. In this work, we performed a genome-wide sequence analysis and identified 97 and 101 putative C2H2-ZF gene family members in the genomes of sweet orange (C. sinensis, poly-embryonic) and pummelo (C. grandis, mono-embryonic), respectively. Phylogenetic analysis categorized citrus C2H2-ZF gene family into four clades, and their possible functions were inferred. According to the numerous regulatory elements on promoter, citrus C2H2-ZFPs can be divided into five different regulatory function types that indicate functional differentiation. RNA-seq data revealed 20 differentially expressed C2H2-ZF genes between poly-embryonic and mono-embryonic ovules at two stages of citrus nucellar embryogenesis, among them CsZFP52 specifically expressed in mono-embryonic pummelo ovules, while CsZFP7, 37, 44, 45, 67 and 68 specifically expressed in poly-embryonic sweet orange ovules. RT-qPCR further validated that CsZFP7 specifically expressed at higher levels in poly-embryonic ovules, and down-regulation of CsZFP7 in the poly-embryonic mini citrus (Fortunella hindsii) increased rate of mono-embryonic seeds compared with the wild type, indicating the regulatory potential of CsZFP7 in nucellar embryogenesis of citrus. This work provided a comprehensive analysis of C2H2-ZF gene family in citrus, including genome organization and gene structure, phylogenetic relationships, gene duplications, possible cis-elements on promoter regions and expression profiles, especially in the poly- and mono-embryogenic ovules, and suggested that CsZFP7 is involved in nucellar embryogenesis.

Acetylome reprograming participates in the establishment of fruit metabolism during polyploidization in citrus.

Polyploidization leads to novel phenotypes and is a major force in evolution. However, the relationship between the evolution of new traits and variations in the post-translational modifications (PTM) of proteins during polyploidization has not been studied. Acetylation of lysine residues is a common protein PTM that plays a critical regulatory role in central metabolism. To test whether changes in metabolism in citrus fruit is associated with the reprogramming of lysine acetylation (Kac) in non-histone proteins during allotetraploidization, we performed a global acetylome analysis of fruits from a synthetic allotetraploid citrus and its diploid parents. A total of 4,175 Kac sites were identified on 1,640 proteins involved in a wide range of fruit traits. In the allotetraploid, parental dominance (i.e. resemblance to one of the two parents) in specific fruit traits, such as fruit acidity and flavonol metabolism, was highly associated with parental Kac level dominance in pertinent enzymes. This association is due to Kac-mediated regulation of enzyme activity. Moreover, protein Kac probably contributes to the discordance between the transcriptomic and proteomic variations during allotetraploidization. The acetylome reprogramming can be partially explained by the expression pattern of several lysine deacetylases (KDACs). Overexpression of silent information regulator 2 (CgSRT2) and histone deacetylase 8 (CgHDA8) diverted metabolic flux from primary metabolism to secondary metabolism and partially restored a metabolic status to the allotetraploid, which expressed attenuated levels of CgSRT2 and CgHDA8. Additionally, KDAC inhibitor treatment greatly altered metabolism in citrus fruit. Collectively, these findings reveal the important role of acetylome reprogramming in trait evolution during polyploidization.

miR156 regulates somatic embryogenesis by modulating starch accumulation in citrus.

Somatic embryogenesis (SE) is a major regeneration approach for in vitro cultured tissues of plants, including citrus. However, SE capability is difficult to maintain, and recalcitrance to SE has become a major obstacle to plant biotechnology. We previously reported that miR156-SPL modules regulate SE in citrus callus. However, the downstream regulatory pathway of the miR156-SPL module in SE remains unclear. In this study, we found that transcription factors CsAGL15 and CsFUS3 bind to the CsMIR156A promoter and activate its expression. Suppression of csi-miR156a function leads to up-regulation of four target genes, SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (CsSPL) genes, and reduction of SE efficiency. In the short tandem target mimic (STTM)-miR156a overexpression callus (MIM156), the number of amyloplasts and starch content were significantly reduced, and genes involved in starch synthesis and transport were down-regulated. csi-miR172d was down-regulated, whereas the target genes, CsTOE1.1 and CsTOE1.2, which inhibit the expression of starch biosynthesis genes, were up-regulated. In our working model, CsAGL15 and CsFUS3 activate csi-miR156a, which represses CsSPLs and further regulates csi-miR172d and CsTOEs, thus altering starch accumulation in callus cells and regulating SE in citrus. This study elucidates the pathway of miR156-SPLs and miR172-TOEs-mediated regulation of SE, and provides new insights into enhancing SE capability in citrus.

miR171 modulates induction of somatic embryogenesis in citrus callus.

KEY MESSAGE: Overexpression of miR171 restored SE competence in the recalcitrant citrus callus, and inhibition of miR171 function weakened SE competence in the strongly embryogenic citrus callus. Somatic embryogenesis (SE) is an important way of in vitro regeneration for plants. For perennial woody crops such as citrus, embryogenic callus is usually induced from unfertilized aborted ovules and widely used in biotechnology aided breeding. However, SE capacity always declines in callus during subculture, which makes regeneration difficult and hinders the application of biotechnology. We previously found that miR171 may be a regulator of SE in citrus, based on the abundant expression of csi-miR171c in the embryogenic callus and during SE of citrus. Here, we report that miR171 promotes SE and is required for SE in citrus. Overexpression of miR171 restored SE competence in the recalcitrant callus of 'Guoqing No.1' Satsuma mandarin (G1), whereas inhibition of miR171 function by Short Tandem Target Mimic (STTM) weakened SE competence in the strongly embryogenic callus of 'Valencia' sweet orange (V). The comparative transcriptomic analysis in miR171 overexpressed callus line (OE) and the wild type callus (WT) indicated that overexpression of miR171 decreased the expression level of its SCARECROW-LIKE (CsSCL) targets, and activated stress response related biological processes and metabolic processes that are required for cell differentiation. However, CsSCLs were up-regulated in the OE callus during SE induction process, which activated the cell division and developmental processes that are required for embryogenesis progress. Our results validate the function of miR171 in regulation of SE and reveal the biological responses provoked by miR171 in citrus that may promote SE.

Transcriptomes and DNA methylomes in apomictic cells delineate nucellar embryogenesis initiation in citrus.

Citrus nucellar poly-embryony (NPE) is a mode of sporophytic apomixis that asexual embryos formed in the seed through adventitious embryogenesis from the somatic nucellar cells. NPE allows clonal propagation of rootstocks, but it impedes citrus cross breeding. To understand the cellular processes involved in NPE initiation, we profiled the transcriptomes and DNA methylomes in laser microdissection captured citrus apomictic cells. In apomictic cells, ribosome biogenesis and protein degradation were activated, whereas auxin polar transport was repressed. Reactive oxygen species (ROS) accumulated in the poly-embryonic ovules, and response to oxidative stress was provoked. The global DNA methylation level, especially that of CHH context, was decreased, whereas the methylation level of the NPE-controlling key gene CitRWP was increased. A C2H2 domain-containing transcription factor gene and CitRWP co-expressed specifically in apomictic cells may coordinate to initiate NPE. The activated embryogenic development and callose deposition processes indicated embryogenic fate of nucellar embryo initial (NEI) cells. In our working model for citrus NPE initiation, DNA hyper-methylation may activate transcription of CitRWP, which increases C2H2 expression and ROS accumulation, triggers epigenetic regulation and regulates cell fate transition and NEI cell identity in the apomictic cells.

GRAS transcription factor LOSS OF AXILLARY MERISTEMS is essential for stamen and runner formation in wild strawberry.

Axillary bud development is a major factor that impacts plant architecture. A runner is an elongated shoot that develops from axillary bud and is frequently used for clonal propagation of strawberry. However, the genetic control underlying runner production is largely unknown. Here, we identified and characterized loss of axillary meristems (lam), an ethyl methanesulfonate-induced mutant of the diploid woodland strawberry (Fragaria vesca) that lacked stamens in flowers and had reduced numbers of branch crowns and runners. The reduced branch crown and runner phenotypes were caused by a failure of axillary meristem initiation. The causative mutation of lam was located in FvH4_3g41310, which encodes a GRAS transcription factor, and was validated by a complementation test. lamCR mutants generated by CRISPR/Cas9 produced flowers without stamens and had fewer runners than the wild-type. LAM was broadly expressed in meristematic tissues. Gibberellic acid (GA) application induced runner outgrowth from the remaining buds in lam, but failed to do so at the empty axils of lam. In contrast, treatment with the GA biosynthesis inhibitor paclobutrazol converted the runners into branch crowns. Moreover, genetic studies indicated that lam is epistatic to suppressor of runnerless (srl), a mutant of FveRGA1 in the GA pathway, during runner formation. Our results demonstrate that LAM is required for stamen and runner formation and acts sequentially with GA from bud initiation to runner outgrowth, providing insights into the molecular regulation of these economically important organs in strawberry.

Identification of a delayed leaf greening gene from a mutation of pummelo.

Delayed greening of young leaves is an unusual phenomenon of plants in nature. Citrus are mostly evergreen tree species. Here, a natural mutant of "Guanxi" pummelo (Citrus maxima), which shows yellow leaves at the young stage, was characterized to identify the genes underlying the trait of delayed leaf greening in plants. A segregating population with this mutant as the seed parent and a normal genotype as the pollen parent was generated. Two DNA pools respectively from the leaves of segregating seedlings with extreme phenotypes of normal leaf greening and delayed leaf greening were collected for sequencing. Bulked segregant analysis (BSA) and InDel marker analysis demonstrated that the delayed leaf greening trait is governed by a 0.3 Mb candidate region on chromosome 6. Gene expression analysis further identified a key candidate gene (Citrus Delayed Greening gene 1, CDG1) in the 0.3 Mb region, which showed significantly differential expression between the genotypes with delayed and normal leaf greening phenotypes. There was a 67 bp InDel region difference in the CDG1 promoter and the InDel region contains a TATA-box element. Confocal laser-scanning microscopy revealed that the CDG1-GFP fusion protein signals were co-localized with the chloroplast signals in the protoplasts. Overexpression of CDG1 in tobacco and Arabidopsis led to the phenotype of delayed leaf greening. These results suggest that the CDG1 gene is involved in controlling the delayed leaf greening phenotype with important functions in chloroplast development.

Localization and characterization of Citrus centromeres by combining half-tetrad analysis and CenH3-associated sequence profiling.

KEY MESSAGE: The physical locations of citrus centromere are revealed by combining genetic and immunological assays for the first time and nine citrus centromere-specific markers for cytogenetics are mined. Centromere localization is challenging, because highly redundant repetitive sequences in centromeric regions make sequence assembly difficult. Although several citrus genomes have been released, the centromeric regions and their characteristics remain to be elucidated. Here, we mapped citrus centromeres through half-tetrad analysis (HTA) that included the genotyping of 54 tetraploid hybrids derived from 2n megagametophytes of Nadorcott tangor with 212 single nucleotide polymorphism (SNP) markers. The sizes of centromeric regions, which estimated based on the heterozygosity restitution rate pattern along the chromosomes, ranged from 1.12 to 18.19 Mb. We also profiled the binding sequences with the centromere-specific histone variant CenH3 by chromatin immunoprecipitation sequencing (ChIP-seq). Based on the positions of the top ten CenH3-enriched contigs, the sizes of centromeric regions were estimated to range from 0.01 to 7.60 Mb and were either adjacent to or included in the centromeric regions identified by HTA. We used DNA probes from two repeats selected from the centromeric regions and seven CenH3-binding centromeric repeats to verify centromeric locations by fluorescence in situ hybridization (FISH). Centromere localization in citrus will contribute to the mining of centromeric/pericentromeric markers, thus to facilitate the rapid identification of mechanisms underlying 2n gamete formation and serve the polyploidy breeding.

The miR399-CsUBC24 Module Regulates Reproductive Development and Male Fertility in Citrus.

MicroRNA399 (miR399) regulates phosphate homeostasis in plants by down-regulating the expression of PHOSPHATE2 (PHO2, or UBC24 encoding the ubiquitin-conjugating E2 enzyme). We previously identified CsmiR399a.1 in a small RNA sequencing screen of a male-sterile somatic cytoplasmic hybrid (or cybrid) of pummelo (Citrus grandis). Here, we report that miR399 affects reproductive development and male fertility in citrus. Down-regulation of CsmiR399a.1 using a short tandem target mimic (STTM) led to abnormal floral development, inhibition of anther dehiscence, and decreased pollen fertility. When grown in inorganic phosphate (Pi)-sufficient conditions, CsmiR399a.1-STTM plants had lower total phosphorus content in their leaves than the wild type and showed typical symptoms of Pi deficiency. In CsmiR399a.1-STTM plants, the expression of genes involved in starch metabolism and Pi homeostasis was significantly different than in the wild type. Thus, we conclude that miR399-STTM mimicked Pi deficiency, disturbed starch metabolism, and was responsible for pollen grain collapse in the transgenic lines. We identified CsUBC24, a citrus homolog of Arabidopsis (Arabidopsis thaliana) AtUBC24 (PHO2), as a target of CsmiR399a.1 that physically interacts with the floral development regulators SEPALLATA family (CsSEP1.1, CsSEP1.2, and CsSEP3) and the anther dehiscence regulator INDUCER OF CBF EXPRESSION1 (CsICE1). We hypothesize that CsUBC24 downregulates the CsSEPs, which disrupts the floral meristem identity regulatory network and leads to developmental abnormalities in flowers. By interacting with CsICE1, CsUBC24 disturbs stomate function on the anther surface, which inhibits anther dehiscence. These findings indicate that a miR399-based mechanism influences both reproductive development and male fertility in citrus.

Genome-wide analysis of the citrus B3 superfamily and their association with somatic embryogenesis.

BACKGROUND: In citrus, genetic improvement via biotechnology is hindered by the obstacle of in vitro regeneration via somatic embryogenesis (SE). Although a few B3 transcription factors are reported to regulate embryogenesis, little is known about the B3 superfamily in citrus, and which members might be involved in SE. RESULTS: Genome-wide sequence analysis identified 72 (CsB3) and 69 (CgB3) putative B3 superfamily members in the genomes of sweet orange (Citrus sinensis, polyembryonic) and pummelo (C. grandis, monoembryonic), respectively. Genome duplication analysis indicated that segmental and tandem duplication events contributed to the expansion of the B3 superfamily in citrus, and that the B3 superfamily evolved under the effect of purifying selection. Phylogenetic relationships were well supported by conserved gene structure and motifs outside the B3 domain, which allowed possible functions to be inferred by comparison with homologous genes from Arabidopsis. Expression analysis identified 23 B3 superfamily members that were expressed during SE in citrus and 17 that may play functional roles at late SE stages. Eight B3 genes were identified that were specific to the genome of polyembryonic sweet orange compared to monoembryonic pummelo. Of these eight B3 genes, CsARF19 was found to be specifically expressed at higher levels in embryogenic callus (EC), implying its possible involvement in EC initiation. CONCLUSIONS: This study provides a genome-wide analysis of the citrus B3 superfamily, including its genome organization, evolutionary features and expression profiles, and identifies specific family members that may be associated with SE.

Polyploidy remodels fruit metabolism by modifying carbon source utilization and metabolic flux in Ponkan mandarin (Citrus reticulata Blanco).

The phenotypic variations that follow polyploidization are expected to improve agricultural productivity and efficiency [1]. However, the effect of polyploidization on plant metabolism has rarely been studied. This study evaluated the metabolic alterations that followed autotetraploidization in the fruit of Ponkan mandarin (C. reticulata Blanco) for three consecutive years and explored the underlying changes to the transcriptome. The autotetraploid (4x) Ponkan fruit had higher levels of total acids, ascorbic acid and total phenolic compounds than the diploid (2x). The primary metabolites especially the organic acids tended to accumulate at higher levels in the 4x fruit. Conversely, two major groups of secondary metabolites (i.e. flavonoids and carotenoids) tended to accumulate at lower levels. The expression levels of citric acid biosynthesis-related genes were unaltered in 4x fruit compared to the 2x fruit. Additionally, genes associated with the transport and utilization of citric acid were significantly down-regulated during ripening, which might induce increases in the levels of citric acid in the 4x fruit. Lower levels of flavonoids and carotenoids in the 4x fruit are potentially associated with decreases in the transport and utilization of citric acid, which is an important metabolite. Citric acid contributes to respiration by serving as an intermediated in the tricarboxylic acid cycle (TCA) and also provides carbon for the production of secondary metabolites. This study demonstrates that polyploidization can influence metabolism in plants.

Overexpression of the CsFUS3 gene encoding a B3 transcription factor promotes somatic embryogenesis in Citrus.

In citrus, genetic improvement via biotechnology is challenging due to insufficient understanding of molecular barriers that prevent regeneration by somatic embryogenesis (SE). Our previous study indicated that LEC genes were involved in SE in citrus, but their regulatory roles remain to be elucidated. Here, we cloned one of the LEC genes, CsFUS3, and show that it is preferentially expressed during SE and in the embryogenic callus (EC) derived from citrus varieties with strong embryogenic competence. The overexpression of CsFUS3 in recalcitrant citrus callus restored embryogenic competence. Complementation of the loss-of-function Arabidopsis fus3 mutant with the CsFUS3 gene restored normal late embryogenesis, which is consistent with the CsFUS3 and AtFUS3 proteins contributing to the same regulatory network in Arabidopsis. Transcriptome profiling revealed that the expression of particular TFs that promote SE was up-regulated in the citrus overexpression (OE) line. The 104 differentially expressed genes associated with hormone biosynthesis, catabolism, and signaling are particularly noteworthy. The dynamic change in the ratio of ABA to GA during SE in wild-type callus mirrored the expression pattern of CsFUS3. In contrast, in the OE line, the ratio of ABA to GA was higher and the capacity for SE was greater when the OE line was separately treated with ABA and GA biosynthesis inhibitors. Taken together, our results demonstrate that the overexpression of CsFUS3 appears to establish a cellular environment favorable to SE, at least in part by promoting a high ABA to GA ratio and by regulating the expression of TFs that promote SE.

miR156-SPL modules regulate induction of somatic embryogenesis in citrus callus.

miR156 is a highly conserved plant miRNA and has been extensively studied because of its versatile roles in plant development. Here, we report a novel role of miR156 in regulating somatic embryogenesis (SE) in citrus, one of the most widely cultivated fruit crops in the world. SE is an important means of in vitro regeneration, but over the course of long-term sub-culturing there is always a decline in the SE potential of the preserved citrus embryogenic callus, and this represents a key obstacle for citrus biotechnology. In this study, the SE competence of citrus callus of wild kumquat (Fortunella hindsii) was significantly enhanced by either overexpression of csi-miR156a or by individual knock-down of the two target genes, CsSPL3 and CsSPL14, indicating that the effect of miR156-SPL modules was established during the initial phases of SE induction. Biological processes that might promote SE in response to miR156 overexpression were explored using RNA-seq, and mainly included hormone signaling pathways, stress responses, DNA methylation, and the cell cycle. CsAKIN10 was identified as interacting protein of CsSPL14. Our results provide insights into the regulatory pathway through which miR156-SPL modules enhance the SE potential of citrus callus, and provide a theoretical basis for improvement of plant SE competence.

Genome-scale mRNA and small RNA transcriptomic insights into initiation of citrus apomixis.

Nucellar embryony (NE) is an adventitious form of apomixis common in citrus, wherein asexual embryos initiate directly from nucellar cells surrounding the embryo sac. NE enables the fixation of desirable agronomic traits and the production of clonal offspring of virus-free rootstock, but impedes progress in hybrid breeding. In spite of the great importance of NE in citrus breeding and commercial production, little is understood about the underlying molecular mechanisms. In this study, the stages of nucellar embryo initiation (NEI) were determined for two polyembryonic citrus cultivars via histological observation. To explore the genes and regulatory pathways involved in NEI, we performed mRNA-seq and sRNA-seq analyses of ovules immediately prior to and at stages during NEI in the two pairs of cultivars. A total of 305 differentially expressed genes (DEGs) were identified between the poly- and monoembryonic ovules. Gene ontology (GO) analysis revealed that several processes are significantly enriched based on DEGs. In particular, response to stress, and especially response to oxidative stress, was over-represented in polyembryonic ovules. Nearly 150 miRNAs, comprising ~90 conserved and ~60 novel miRNAs, were identified in the ovules of either cultivar pair. Only two differentially expressed miRNAs (DEMs) were identified, of which the novel miRN23-5p was repressed whereas the targets accumulated in the polyembryonic ovules. This integrated study on the transcriptional and post-transcriptional regulatory profiles between poly- and monoembryonic citrus ovules provides new insights into the mechanism of NE, which should contribute to revealing the regulatory mechanisms of plant apomixis.

High-throughput sequencing and degradome analysis reveal altered expression of miRNAs and their targets in a male-sterile cybrid pummelo (Citrus grandis).

BACKGROUND: G1 + HBP is a male sterile cybrid line with nuclear genome from Hirado Buntan pummelo (C. grandis Osbeck) (HBP) and mitochondrial genome from "Guoqing No.1" (G1, Satsuma mandarin), which provides a good opportunity to study male sterility and nuclear-cytoplasmic cross talk in citrus. High-throughput sRNA and degradome sequencing were applied to identify miRNAs and their targets in G1 + HBP and its fertile type HBP during reproductive development. RESULTS: A total of 184 known miRNAs, 22 novel miRNAs and 86 target genes were identified. Some of the targets are transcription factors involved in floral development, such as auxin response factors (ARFs), SQUAMOSA promoter binding protein box (SBP-box), MYB, basic region-leucine zipper (bZIP), APETALA2 (AP2) and transport inhibitor response 1 (TIR1). Eight target genes were confirmed to be sliced by corresponding miRNAs using 5' RACE technology. Based on the sequencing abundance, 42 differentially expressed miRNAs between sterile line G1 + HBP and fertile line HBP were identified. Differential expression of miRNAs and their target genes between two lines was validated by quantitative RT-PCR, and reciprocal expression patterns between some miRNAs and their targets were demonstrated. The regulatory mechanism of miR167a was investigated by yeast one-hybrid and dual-luciferase assays that one dehydrate responsive element binding (DREB) transcription factor binds to miR167a promoter and transcriptionally repress miR167 expression. CONCLUSION: Our study reveals the altered expression of miRNAs and their target genes in a male sterile line of pummelo and highlights that miRNA regulatory network may be involved in floral bud development and cytoplasmic male sterility in citrus.

Fragile Sites of 'Valencia' Sweet Orange (Citrus sinensis) Chromosomes Are Related with Active 45s rDNA.

Citrus sinensis chromosomes present a morphological differentiation of bands after staining by the fluorochromes CMA and DAPI, but there is still little information on its chromosomal characteristics. In this study, the chromosomes in 'Valencia' C. sinensis were analyzed by fluorescence in situ hybridization (FISH) using telomere DNA and the 45S rDNA gene as probes combining CMA/DAPI staining, which showed that there were two fragile sites in sweet orange chromosomes co-localizing at distended 45S rDNA regions, one proximally locating on B-type chromosome and the other subterminally locating on D-type chromosome. While the chromosomal CMA banding and 45S rDNA FISH mapping in the doubled haploid line of 'Valencia' C. sinensis indicated six 45S rDNA regions, four were identified as fragile sites as doubled comparing its parental line, which confirmed the cytological heterozygosity and chromosomal heteromorphisms in sweet orange. Furthermore, Ag-NOR identified two distended 45S rDNA regions to be active nucleolar organizing regions (NORs) in diploid 'Valencia' C. sinensis. The occurrence of quadrivalent in meiosis of pollen mother cells (PMCs) in 'Valencia' sweet orange further confirmed it was a chromosomal reciprocal translocation line. We speculated this chromosome translocation was probably related to fragile sites. Our data provide insights into the chromosomal characteristics of the fragile sites in 'Valencia' sweet orange and are expected to facilitate the further investigation of the possible functions of fragile sites.