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

The role of microRNAs in responses to drought and heat stress in peanut (Arachis hypogaea)

MicroRNAs (miRNAs) are 21-24 nt small RNAs (sRNAs) that negatively regulate protein-coding genes and/or trigger phased small-interfering RNA (phasiRNA) production. Two thousand nine hundred miRNA families, of which ∼40 are deeply conserved, have been identified in ∼80 different plant species genomes. miRNA functions in response to abiotic stresses is less understood than their roles in development. Only seven peanut MIRNA families are documented in miRBase, yet a reference genome assembly is now published and over 480 plant-like MIRNA loci were predicted in the diploid peanut progenitor Arachis duranensis genome. We explored by computational analysis of a leaf sRNA library and publicly available sRNA, degradome, and transcriptome datasets the miRNA and phasiRNA space associated with drought and heat stresses in peanut. We characterized 33 novel candidate and 33 ancient conserved families of MIRNAs and present degradome evidence for their cleavage activities on mRNA targets, including several noncanonical targets and novel phasiRNA-producing noncoding and mRNA loci with validated novel targets such as miR1509 targeting serine/threonine-protein phosphatase7 and miRc20 and ahy-miR3514 targeting penta-tricopeptide repeats (PPRs), in contradistinction to other claims of miR1509/173/7122 superfamily miRNAs indirectly targeting PPRs via TAS-like noncoding RNA loci. We characterized the inverse correlations of significantly differentially expressed drought- and heat-regulated miRNAs, assayed by sRNA blots or transcriptome datasets, with target mRNA expressions in the same datasets. Meta-analysis of an expression atlas and over representation of miRNA target genes in co-expression networks suggest that miRNAs have functions in unique aspects of peanut gynophore development. Genome-wide MIRNA annotation of the published allopolyploid peanut genome can facilitate molecular breeding of value-added traits.

Apomixis: oh, what a tangled web we have!

Apomixis is a complex evolutionary trait with many possible origins. Here we discuss various clues and causes, ultimately proposing a model harmonizing the three working hypotheses on the topic. Asexual reproduction through seeds, i.e., apomixis, is the holy grail of plant biology. Its implementation in modern breeding could be a game-changer for agriculture. It has the potential to generate clonal crops and maintain valuable complex genotypes and their associated heterotic traits without inbreeding depression. The genetic basis and origins of apomixis are still unclear. There are three central hypothesis for the development of apomixis that could be: i) a deviation from the sexual developmental program caused by an asynchronous development, ii) environmentally triggered through epigenetic regulations (a polyphenism of sex), iii) relying on one or more genes/alleles. Because of the ever-increasing complexity of the topic, the path toward a detailed understanding of the mechanisms underlying apomixis remains unclear. Here, we discuss the most recent advances in the evolution perspective of this multifaceted trait. We incorporated our understanding of the effect of endogenous effectors, such as small RNAs, epigenetic regulation, hormonal pathways, protein turnover, and cell wall modification in response to an upside stress. This can be either endogenous (hybridization or polyploidization) or exogenous environmental stress, mainly due to oxidative stress and the corresponding ROS (Reacting Oxygen Species) effectors. Finally, we graphically represented this tangled web.

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.

Data mining of transcriptional biomarkers at different cotton fiber developmental stages

Advancement of the gene expression study provides comprehensive information on pivotal genes at different cotton fiber development stages. For the betterment of cotton fiber yield and their quality, genetic improvement is a major target point for the cotton community. Therefore, various studies were carried out to understand the transcriptional machinery of fiber leading to the detailed integrative as well as innovative study. Through data mining and statistical approaches, we identified and validated the transcriptional biomarkers for staged specific differentiation of fiber. With the unique mapping read matrix of ~ 200 cotton transcriptome data and sequential statistical analysis, we identified several important genes that have a deciding and specific role in fiber cell commitment, initiation and elongation, or secondary cell wall synthesis stage. Based on the importance score and validation analysis, IQ domain 26, Aquaporin, Gibberellin regulated protein, methionine gamma lyase, alpha/beta hydrolases, and HAD-like superfamily have shown the specific and determining role for fiber developmental stages. These genes are represented as transcriptional biomarkers that provide a base for molecular characterization for cotton fiber development which will ultimately determine the high yield.

Analysis of circRNAs and circRNA-associated competing endogenous RNA networks in β-thalassemia

The involvement of circRNAs in β-thalassemia and their actions on fetal hemoglobin (HbF) is unclear. Here, the circRNAs in β-thalassemia carriers with high HbF levels were comprehensively analyzed and compared with those of healthy individuals. Differential expression of 2183 circRNAs was observed and their correlations with hematological parameters were investigated. Down-regulated hsa-circRNA-100466 had a strong negative correlation with HbF and HbA₂. Bioinformatics was employed to construct a hsa-circRNA-100466‑associated competing endogenous RNA (ceRNA) network to identify hub genes and associated miRNAs. The hsa-circRNA-100466▁miR-19b-3p▁SOX6 pathway was identified using both present and previously published data. The ceRNA network was verified by qRT-PCR analysis of β-thalassemia samples, RNA immunoprecipitation of K562 cell lysates, and dual-luciferase reporter analysis. qRT-PCR confirmed that hsa-circRNA-100466 and SOX6 were significantly down-regulated, while miR-19b-3p was up-regulated. Hsa-circRNA-100466, miR-19b-3p, and SOX6 were co-immunoprecipitated by anti-argonaute antibodies, indicating involvement with HbF induction. A further dual-luciferase reporter assay verified that miR-19b-3p interacted directly with hsa-circRNA-100466 and SOX6. Furthermore, spearman correlation coefficients revealed their significant correlations with HbF. In conclusion, a novel hsa-circRNA-100466▁miR-19b-3p▁SOX6 pathway was identified, providing insight into HbF induction and suggesting targets β-thalassemia treatment.

The emerging role of small RNAs in ovule development, a kind of magic

In plants, small RNAs have been recognized as key genetic and epigenetic regulators of development. Small RNAs are usually 20 to 30 nucleotides in length and they control, in a sequence specific manner, the transcriptional or post-transcriptional expression of genes. In this review, we present a comprehensive overview of the most recent findings about the function of small RNAs in ovule development, including megasporogenesis and megagametogenesis, both in sexual and apomictic plants. We discuss recent studies on the role of miRNAs, siRNAs and trans-acting RNAs (ta-siRNAs) in early female germline differentiation. The mechanistic complexity and unique regulatory features are reviewed, and possible directions for future research are provided.

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.

MicroRNAs play important roles in regulating the rapid growth of the Phyllostachys edulis culm internode

Moso bamboo (Phyllostachys edulis) is a fast-growing species with uneven growth and lignification from lower to upper segments within one internode. MicroRNAs (miRNAs) play a vital role in post-transcriptional regulation in plants. However, how miRNAs regulate fast growth in bamboo internodes is poorly understood. In this study, one moso bamboo internode was divided during early rapid growth into four segments called F4 (bottom) to F1 (upper) and these were then analysed for transcriptomes, miRNAs and degradomes. The F4 segment had a higher number of actively dividing cells as well as a higher content of auxin (IAA), cytokinin (CK) and gibberellin (GA) compared with the F1 segment. RNA-seq analysis showed DNA replication and cell division-associated genes highly expressed in F4 rather than in F1. In total, 63 miRNAs (DEMs) were identified as differentially expressed between F4 and F1. The degradome and the transcriptome indicated that many downstream transcription factors and hormonal responses genes were modulated by DEMs. Several miR-target interactions were further validated by tobacco co-infiltration. Our findings give new insights into miRNA-mediated regulatory pathways in bamboo, and will contribute to a comprehensive understanding of the molecular mechanisms governing rapid growth.

Chromosome-level genome assembly of Sichuan pepper provides insights into apomixis, drought tolerance, and alkaloid biosynthesis

Sichuan pepper is a commonly used spice in Asian cuisine. Sanshools and wgx-50/gx-50 isolated from it have been shown to possess a wide spectrum of medicinal properties. Here we generated a chromosome-level genome assembly of one Sichuan pepper species Zanthoxylum armatum characterized by drought tolerance and apomixis. Analyses of functionally related genes suggested that increased gene copy number and expression level of drought-tolerant genes might play an important role in improving drought tolerance of Z. armatum. Moreover, a gene encoding an RWP-RK domain-containing protein was shown to contribute to apomixis in Z. armatum, which was further characterized by overexpression in Arabidopsis thaliana. Furthermore, based on gene homology searching and co-expression patterns of metabolite concentration and gene expressions, we identified a number of candidate genes involved in the biosynthesis of sanshools and wgx-50/gx-50. Taken together, our results yield valuable insights for understanding the evolution of apomixis, drought tolerance, and alkaloid biosynthesis in Z. armatum.

Small RNA sequencing provides candidate miRNA-target pairs for revealing the mechanism of apomixis in Zanthoxylum bungeanum

BACKGROUND

Apomixis is a form of asexual reproduction that produces offspring without the need for combining male and female gametes, and the offspring have the same genetic makeup as the mother. Therefore, apomixis technology has great application potential in plant breeding. To identify the apomixis types and critical period, embryonic development at different flower development stages of Zanthoxylum bungeanum was observed by cytology.

RESULTS

The results show that the S3 stage is the critical period of apomixis, during which the nucellar cells develop into an adventitious primordial embryo. Cytological observations showed that the type of apomixis in Z. bungeanum is sporophytic apomixis. Furthermore, miRNA sequencing, miRNA-target gene interaction, dual luciferase reporter assay, and RT-qPCR verification were used to reveal the dynamic regulation of miRNA-target pairs in Z. bungeanum apomixis. The miRNA sequencing identified 96 mature miRNAs, of which 40 were known and 56 were novel. Additionally, 29 differentially expressed miRNAs were screened according to the miRNAs expression levels at the different developmental stages. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analyses showed that the target genes of the differentially expressed miRNAs were mainly enriched in plant hormone signal transduction, RNA biosynthetic process, and response to hormone pathways.

CONCLUSIONS

During the critical period of apomictic embryonic development, miR172c significantly reduces the expression levels of TOE3 and APETALA 2 (AP2) genes, thereby upregulating the expression of the AGAMOUS gene. A molecular regulation model of miRNA-target pairs was constructed based on their interactions and expression patterns to further understand the role of miRNA-target pairs in apomixis. Our data suggest that miR172c may regulates AGAMOUS expression by inhibiting TOE3 in the critical period of apomixis.

A study of the heterochronic sense/antisense RNA representation in florets of sexual and apomictic Paspalum notatum

Background

Apomixis, an asexual mode of plant reproduction, is a genetically heritable trait evolutionarily related to sexuality, which enables the fixation of heterozygous genetic combinations through the development of maternal seeds. Recently, reference floral transcriptomes were generated from sexual and apomictic biotypes of Paspalum notatum, one of the most well-known plant models for the study of apomixis. However, the transcriptome dynamics, the occurrence of apomixis vs. sexual expression heterochronicity across consecutive developmental steps and the orientation of transcription (sense/antisense) remain unexplored.

Results

We produced 24 Illumina TruSeq®/ Hiseq 1500 sense/antisense floral transcriptome libraries covering four developmental stages (premeiosis, meiosis, postmeiosis, and anthesis) in biological triplicates, from an obligate apomictic and a full sexual genotype. De novo assemblies with Trinity yielded 103,699 and 100,114 transcripts for the apomictic and sexual samples respectively. A global comparative analysis involving reads from all developmental stages revealed 19,352 differentially expressed sense transcripts, of which 13,205 (68%) and 6147 (32%) were up- and down-regulated in apomictic samples with respect to the sexual ones. Interestingly, 100 differentially expressed antisense transcripts were detected, 55 (55%) of them up- and 45 (45%) down-regulated in apomictic libraries. A stage-by-stage comparative analysis showed a higher number of differentially expressed candidates due to heterochronicity discrimination: the highest number of differential sense transcripts was detected at premeiosis (23,651), followed by meiosis (22,830), postmeiosis (19,100), and anthesis (17,962), while the highest number of differential antisense transcripts were detected at anthesis (495), followed by postmeiosis (164), meiosis (120) and premeiosis (115). Members of the AP2, ARF, MYB and WRKY transcription factor families, as well as the auxin, jasmonate and cytokinin plant hormone families appeared broadly deregulated. Moreover, the chronological expression profile of several well-characterized apomixis controllers was examined in detail.

Conclusions

This work provides a quantitative sense/antisense gene expression catalogue covering several subsequent reproductive developmental stages from premeiosis to anthesis for apomictic and sexual P. notatum, with potential to reveal heterochronic expression between reproductive types and discover sense/antisense mediated regulation. We detected a contrasting transcriptional and hormonal control in apomixis and sexuality as well as specific sense/antisense modulation occurring at the onset of parthenogenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07450-3.

Small Non-Coding RNAs at the Crossroads of Regulatory Pathways Controlling Somatic Embryogenesis in Seed Plants

Small non-coding RNAs (sncRNAs) are molecules with important regulatory functions during development and environmental responses across all groups of terrestrial plants. In seed plants, the development of a mature embryo from the zygote follows a synchronized cell division sequence, and growth and differentiation events regulated by highly regulated gene expression. However, given the distinct features of the initial stages of embryogenesis in gymnosperms and angiosperms, it is relevant to investigate to what extent such differences emerge from differential regulation mediated by sncRNAs. Within these, the microRNAs (miRNAs) are the best characterized class, and while many miRNAs are conserved and significantly represented across angiosperms and other seed plants during embryogenesis, some miRNA families are specific to some plant lineages. Being a model to study zygotic embryogenesis and a relevant biotechnological tool, we systematized the current knowledge on the presence and characterization of miRNAs in somatic embryogenesis (SE) of seed plants, pinpointing the miRNAs that have been reported to be associated with SE in angiosperm and gymnosperm species. We start by conducting an overview of sncRNA expression profiles in the embryonic tissues of seed plants. We then highlight the miRNAs described as being involved in the different stages of the SE process, from its induction to the full maturation of the somatic embryos, adding references to zygotic embryogenesis when relevant, as a contribution towards a better understanding of miRNA-mediated regulation of SE.

Regulation of nucellar embryony, a mode of sporophytic apomixis in Citrus resembling somatic embryogenesis

Apomixis, is an asexual mode of seed formation resulting in genetically identical or clonal seed with a maternal genotype. Apomixis has not been reported in seed crops where its flexible application in plant breeding could accelerate delivery of new varieties. By contrast, a sporophytic form of apomixis termed nucellar or adventitious embryony is common in the Rutaceae containing Citrus crop species. Here, multiple embryos develop from the maternal, somatic, nucellar cells of the ovule. They are incorporated into the enlarging embryo sac containing the sexually derived zygotic embryo and endosperm, which are products of double fertilization. Recent research has provided insights to the molecular basis for nucellar embryony. Here, we review the current understanding of the initiation, genetic basis and evolution of nucellar embryony in Citrus, and discuss prospects for future study and breeding applications of Citrus sporophytic apomixis.

ZbAGL11, a class D MADS-box transcription factor of Zanthoxylum bungeanum, is involved in sporophytic apomixis

Apomixis is a reproductive model that bypasses sexual reproduction, so it does not require the combination of paternal and maternal gametes but instead results in the production of offspring directly from maternal tissues. This reproductive mode results in the same genetic material in the mother and the offspring and has significant applications in agricultural breeding. Molecular and cytological methods were used to identify the reproductive type of Zanthoxylum bungeanum (ZB). Fluorescence detection of the amplified products of 12 pairs of polymorphic SSR primers showed consistent fluorescence signals for mother and offspring, indicating that no trait separation occurred during reproduction. In addition, the cytological observation results showed differentiation of ZB embryos (2n) from nucellar cells (2n) to form indefinite embryonic primordia and then form adventitious embryos (2n), indicating that the apomictic type of ZB is sporophytic apomixis. The MADS-box transcription factor ZbAGL11 was highly expressed during the critical period of nucellar embryo development in ZB. Unpollinated ZbAGL11-OE Arabidopsis produced fertile offspring and exhibited an apomictic phenotype. The overexpression of ZbAGL11 increased the callus induction rate of ZB tissue. In addition, the results of the yeast two-hybrid experiment showed that ZbAGL11 could interact with the ZbCYP450 and ZbCAD11 proteins. Our results demonstrate that ZbAGL11 can cause developmental disorders of Arabidopsis flower organs and result in apomixis-like phenotypes.

Genome-wide sRNA and mRNA transcriptomic profiling insights into dynamic regulation of taproot thickening in radish (Raphanus sativus L.)

BACKGROUND

Taproot is the main edible organ and ultimately determines radish yield and quality. However, the precise molecular mechanism underlying taproot thickening awaits further investigation in radish. Here, RNA-seq was performed to identify critical genes involved in radish taproot thickening from three advanced inbred lines with different root size.

RESULTS

A total of 2606 differentially expressed genes (DEGs) were shared between 'NAU-DY' (large acicular) and 'NAU-YB' (medium obovate), which were significantly enriched in 'phenylpropanoid biosynthesis', 'glucosinolate biosynthesis', and 'starch and sucrose metabolism' pathway. Meanwhile, a total of 16 differentially expressed miRNAs (DEMs) were shared between 'NAU-DY' and 'NAU-YH' (small circular), whereas 12 miRNAs exhibited specific differential expression in 'NAU-DY'. Association analysis indicated that miR393a-bHLH77, miR167c-ARF8, and miR5658-APL might be key factors to biological phenomenon of taproot type variation, and a putative regulatory model of taproot thickening and development was proposed. Furthermore, several critical genes including SUS1, EXPB3, and CDC5 were characterized and profiled by RT-qPCR analysis.

CONCLUSION

This integrated study on the transcriptional and post-transcriptional profiles could provide new insights into comprehensive understanding of the molecular regulatory mechanism underlying taproot thickening in root vegetable crops.

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.

Controlling Apomixis: Shared Features and Distinct Characteristics of Gene Regulation

In higher plants, sexual and asexual reproduction through seeds (apomixis) have evolved as alternative strategies. As apomixis leads to the formation of clonal offspring, its great potential for agricultural applications has long been recognized. However, the genetic basis and the molecular control underlying apomixis and its evolutionary origin are to date not fully understood. Both in sexual and apomictic plants, reproduction is tightly controlled by versatile mechanisms regulating gene expression, translation, and protein abundance and activity. Increasing evidence suggests that interrelated pathways including epigenetic regulation, cell-cycle control, hormonal pathways, and signal transduction processes are relevant for apomixis. Additional molecular mechanisms are being identified that involve the activity of DNA- and RNA-binding proteins, such as RNA helicases which are increasingly recognized as important regulators of reproduction. Together with other factors including non-coding RNAs, their association with ribosomes is likely to be relevant for the formation and specification of the apomictic reproductive lineage. Subsequent seed formation appears to involve an interplay of transcriptional activation and repression of developmental programs by epigenetic regulatory mechanisms. In this review, insights into the genetic basis and molecular control of apomixis are presented, also taking into account potential relations to environmental stress, and considering aspects of evolution.

Small RNA and mRNA Sequencing Reveal the Roles of microRNAs Involved in Pomegranate Female Sterility

Female sterility is a key factor restricting plant reproduction. Our previous studies have revealed that pomegranate female sterility mainly arose from the abnormality of ovule development. MicroRNAs (miRNAs) play important roles in ovule development. However, little is known about the roles of miRNAs in female sterility. In this study, a combined high-throughput sequencing approach was used to investigate the miRNAs and their targeted transcripts involved in female development. A total of 103 conserved and 58 novel miRNAs were identified. Comparative profiling indicated that the expression of 43 known miRNAs and 14 novel miRNAs were differentially expressed between functional male flowers (FMFs) and bisexual flowers (BFs), 30 known miRNAs and nine novel miRNAs showed significant differences among different stages of BFs, and 20 known miRNAs and 18 novel miRNAs exhibited remarkable expression differences among different stages of FMFs. Gene ontology (GO) analyses of 144 predicted targets of differentially expressed miRNAs indicated that the “reproduction process” and “floral whorl development” processes were significantly enriched. The miRNA–mRNA interaction analyses revealed six pairs of candidate miRNAs and their targets associated with female sterility. Interestingly, pg-miR166a-3p was accumulated, whereas its predicted targets (Gglean012177.1 and Gglean013966.1) were repressed in functional male flowers (FMFs), and the interaction between pg-miR166a-3p and its targets (Gglean012177.1 and Gglean013966.1) were confirmed by transient assay. A. thaliana transformed with 35S-pre-pg-miR166a-3p verified the role of pg-miR166a-3p in ovule development, which indicated pg-miR166a-3p’s potential role in pomegranate female sterility. The results provide new insights into molecular mechanisms underlying the female sterility at the miRNA level.

LINC00459 sponging miR-218 to elevate DKK3 inhibits proliferation and invasion in melanoma

The lncRNA biomarkers in melanoma remain to be further explored. The lncRNAs with different expression levels in melanoma tissue were identified by microarray analysis. To investigate the biological functions of target lncRNA, several in-vivo and in-vitro studies were performed. Potential mechanisms of competitive endogenous RNAs (ceRNAs) were predicted by using bioinformatics analysis and explored by western blot assay, fluorescence in situ hybridization assay, real-time quantitative PCR (RT-qPCR) array, RNA pull-down analysis, AGO2-RIP assay, and dual-luciferase reporter assay. The results demonstrated decreased LINC00459 in melanoma cell lines and tissues. According to the in-vitro and in-vivo experiments, up-regulated LINC00459 had inhibitory effect on cell proliferation and invasion. Bioinformatics analyses suggested that miR-218 could be a direct target of LINC00459. In addition, miR-218 was proved to be able to directly target the dickkopf-related protein 3 (DKK3) gene. In conclusion, our analysis suggested that the LINC00459 could sponge miR-218 and increase the expression of DKK3 gene, thus inhibiting the invasion and proliferation of melanoma cells, which indicated that the LINC00459 could be an effective biomarker for melanoma and its potential as the therapeutic target.

Integration of Transcriptomes, Small RNAs, and Degradome Sequencing to Identify Putative miRNAs and Their Targets Related to Eu-Rubber Biosynthesis in Eucommia ulmoides

Eucommia ulmoides has attracted much attention as a valuable natural rubber (Eu-rubber) production tree. As a strategic material, Eu-rubber plays a vital role in general and defence industries. However, the study of Eu-rubber biosynthesis at a molecular level is scarce, and the regulatory network between microRNAs (miRNAs) and messenger RNAs (mRNAs) in Eu-rubber biosynthesis has not been assessed. In this study, we comprehensively analyzed the transcriptomes, small RNAs (sRNAs) and degradome to reveal the regulatory network of Eu-rubber biosynthesis in E. ulmoides. A total of 82,065 unigenes and 221 miRNAs were identified using high-throughput sequencing; 20,815 targets were predicted using psRNATarget software. Of these targets, 779 miRNA-target pairs were identified via degradome sequencing. Thirty-one miRNAs were differentially expressed; 22 targets of 34 miRNAs were annotated in the terpenoid backbone biosynthesis pathway (ko00900) based on the Kyoto Encyclopedia of Genes and Genomes (KEGG). These miRNAs were putatively related to Eu-rubber biosynthesis. A regulatory network was constructed according to the expression profiles of miRNAs and their targets. These results provide a comprehensive analysis of transcriptomics, sRNAs and degradome to reveal the Eu-rubber accumulation, and provide new insights into genetic engineering techniques which may improve the content of Eu-rubber in E. ulmoides.

Small RNA-seq reveals novel regulatory components for apomixis in Paspalum notatum

Background

Apomixis is considered an evolutionary deviation of the sexual reproductive pathway leading to the generation of clonal maternal progenies by seeds. Recent evidence from model and non-model species suggested that this trait could be modulated by epigenetic mechanisms involving small RNAs (sRNAs). Here we profiled floral sRNAs originated from apomictic and sexual Paspalum notatum genotypes in order to identify molecular pathways under epigenetic control that might be involved in the transition from sexuality to agamospermy.

Results

The mining of genes participating in sRNA-directed pathways from floral Paspalum transcriptomic resources showed these routes are functional during reproductive development, with several members differentially expressed in apomictic and sexual plants. Triplicate floral sRNA libraries derived from apomictic and a sexual genotypes were characterized by using high-throughput sequencing technology. EdgeR was apply to compare the number of sRNA reads between sexual and apomictic libraries that map over all Paspalum floral transcripts. A total of 1525 transcripts showed differential sRNA representation, including genes related to meiosis, plant hormone signaling, biomolecules transport, transcription control and cell cycle. Survey for miRNA precursors on transcriptome and genome references allowed the discovery of 124 entities, including 40 conserved and 8 novel ones. Fifty-six clusters were differentially represented in apomictic and sexual plants. All differentially expressed miRNAs were up-regulated in apomictic libraries but miR2275, which showed different family members with opposed representation. Examination of predicted miRNAs targets detected 374 potential candidates. Considering sRNA, miRNAs and target surveys together, 14 genes previously described as related with auxin metabolism, transport and signaling were detected, including AMINO ACID/AUXIN PERMEASE 15, IAA-AMIDO SYNTHETASE GH3–8, IAA30, miR160, miR167, miR164, miR319, ARF2, ARF8, ARF10, ARF12, AFB2, PROLIFERATING CELL FACTOR 6 and NITRATE TRANSPORTER 1.1.

Conclusions

This work provides a comprehensive survey of the sRNA differential representation in flowers of sexual and apomictic Paspalum notatum plants. An integration of the small RNA profiling data presented here and previous transcriptomic information suggests that sRNA-mediated regulation of auxin pathways is pivotal in promoting apomixis. These results will underlie future functional characterization of the molecular components mediating the switch from sexuality to apomixis.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5881-0) contains supplementary material, which is available to authorized users.

LINC00612 enhances the proliferation and invasion ability of bladder cancer cells as ceRNA by sponging miR-590 to elevate expression of PHF14

Background

Bladder cancer (BC) is a common type of cancer that involves tumors of the urinary system and poses a serious threat to human health. Long noncoding RNAs (lncRNAs) have emerged as crucial biomarkers and regulators in many cancers. Novel lncRNA biomarkers in BC urgently need to be investigated in regard to its function and regulatory mechanisms.

Methods

Identification of differentially expressed lncRNAs in BC tissue was performed via microarray analysis. To investigate the biological functions of LINC00612, loss-of-function and gain-of-function experiments were performed in vitro and in vivo. Bioinformatics analysis, dual-luciferase reporter assays, AGO2-RIP assays, RNA pull-down assays, real-time quantitative PCR (RT-qPCR) arrays, fluorescence in situ hybridization assays, and western blot assays were conducted to explore the underlying mechanisms of competitive endogenous RNAs (ceRNAs).

Results

LINC00612 was upregulated in BC tissues and cell lines. Functionally, downregulation of LINC00612 inhibited cell proliferation and invasion in vitro and in vivo, whereas overexpression of LINC00612 resulted in the opposite effects. Bioinformatics analysis and luciferase assays revealed that miR-590 was a direct target of LINC0061, which was validated by dual-luciferase reporter assays, AGO2-RIP assays, RNA pull-down assays, RT-qPCR arrays, and rescue experiments. Additionally, miR-590 was shown to directly target the PHD finger protein 14 (PHF14) gene. LNIC00612 modulated the expression of E-cadherin and vimentin by competitively sponging miR-590 to elevate the expression of PHF14, thus affecting BC cellular epithelial-mesenchymal transition (EMT).

Conclusions

Our results indicate that LINC00612 enhances the proliferation and invasion ability of BC cells by sponging miR-590 to upregulate PHF14 expression and promote BC cellular EMT, suggesting that LINC00612 may act as a potential biomarker and therapeutic target for BC.

Lnc SMAD5-AS1 as ceRNA inhibit proliferation of diffuse large B cell lymphoma via Wnt/β-catenin pathway by sponging miR-135b-5p to elevate expression of APC

Diffuse large B cell lymphoma (DLBCL) is a common and fatal hematological malignancy. Long noncoding RNAs (lncRNAs) have emerged as crucial biomarkers and regulators in many cancers. Novel lncRNA biomarker in DLBCL needs to be investigated badly, as well as its function and molecular mechanism. To further explore, microarray analysis was performed to identify the differentially expressed lncRNAs in DLBCL tissues. To investigate the biological functions of SMAD5-AS1, we performed gain- and loss-of-function experiments in vitro and in vivo. Furthermore, bioinformatics analysis, dual-luciferase reporter assays, Argonaute 2-RNA immunoprecipitation (AGO2-RIP), RNA pull-down assay, quantitative PCR arrays, western blot assay, TOPFlash/FOPFlash reporter assay, and rescue experiments were conducted to explore the underlying mechanisms of competitive endogenous RNAs (ceRNAs). We found that SMAD5-AS1 was down-regulated in DLBCL tissues and cell lines. Functionally, SMAD5-AS1 downregulation promoted cell proliferation in vitro and in vivo, whereas SMAD5-AS1 overexpression could lead to the opposite effects in vitro and in vivo. Bioinformatics analysis and luciferase assays revealed that miR-135b-5p was a direct target of SMAD5-AS1, which was validated by dual-luciferase reporter assays, AGO2-RIP, RNA pull-down assay, and rescue experiments. Also, dual-luciferase reporter assays and rescue experiments demonstrated that miR-135b-5p targeted the adenomatous polyposis coli (APC) gene directly. SMAD5-AS1/miR-135b-5p inhibits the cell proliferation via inactivating the classic Wnt/β-catenin pathway in the form of APC dependency. Our results indicated that SMAD5-AS1 inhibits DLBCL proliferation by sponging miR-135b-5p to up-regulate APC expression and inactivate classic Wnt/β-catenin pathway, suggesting that SMAD5-AS1 may act as a potential biomarker and therapeutic target for DLBCL.

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.

MITE insertion-dependent expression of CitRKD1 with a RWP-RK domain regulates somatic embryogenesis in citrus nucellar tissues

BACKGROUND

Somatic embryogenesis in nucellar tissues is widely recognized to induce polyembryony in major citrus varieties such as sweet oranges, satsuma mandarins and lemons. This capability for apomixis is attractive in agricultural production systems using hybrid seeds, and many studies have been performed to elucidate the molecular mechanisms of various types of apomixis. To identify the gene responsible for somatic embryogenesis in citrus, a custom oligo-DNA microarray including predicted genes in the citrus polyembryonic locus was used to compare the expression profiles in reproductive tissues between monoembryonic and polyembryonic varieties. The full length of CitRKD1, which was identified as a candidate gene responsible for citrus somatic embryogenesis, was isolated from satsuma mandarin and its molecular function was investigated using transgenic 'Hamlin' sweet orange by antisense-overexpression.

RESULTS

The candidate gene CitRKD1, predominantly transcribed in reproductive tissues of polyembryonic varieties, is a member of the plant RWP-RK domain-containing protein. CitRKD1 of satsuma mandarin comprised two alleles (CitRKD1-mg1 and CitRKD1-mg2) at the polyembryonic locus controlling embryonic type (mono/polyembryony) that were structurally divided into two types with or without a miniature inverted-repeat transposable element (MITE)-like insertion in the upstream region. CitRKD1-mg2 with the MITE insertion was the predominant transcript in flowers and young fruits where somatic embryogenesis of nucellar cells occurred. Loss of CitRKD1 function by antisense-overexpression abolished somatic embryogenesis in transgenic sweet orange and the transgenic T₁ plants were confirmed to derive from zygotic embryos produced by self-pollination by DNA diagnosis. Genotyping PCR analysis of 95 citrus traditional and breeding varieties revealed that the CitRKD1 allele with the MITE insertion (polyembryonic allele) was dominant and major citrus varieties with the polyembryonic allele produced polyembryonic seeds.

CONCLUSION

CitRKD1 at the polyembryonic locus plays a principal role in regulating citrus somatic embryogenesis. CitRKD1 comprised multiple alleles that were divided into two types, polyembryonic alleles with a MITE insertion in the upstream region and monoembryonic alleles without it. CitRKD1 was transcribed in reproductive tissues of polyembryonic varieties with the polyembryonic allele. The MITE insertion in the upstream region of CitRKD1 might be involved in regulating the transcription of CitRKD1.

Reproduction in woody perennial Citrus: an update on nucellar embryony and self-incompatibility

Review on citrus reproduction. Citrus is one of the most important and widely grown fruit crops. It possesses several special reproductive characteristics, such as nucellar embryony and self-incompatibility. The special phenomenon of nucellar embryony in citrus, also known as the polyembryony, is a kind of sporophytic apomixis. During the past decade, the emergence of novel technologies and the construction of multiple citrus reference genomes have facilitated rapid advances to our understanding of nucellar embryony. Indeed, several research teams have preliminarily determined the genetic basis of citrus apomixis. On the other hand, the phenomenon of self-incompatibility that promotes genetic diversity by rejecting self-pollen and accepting non-self-pollen is difficult to study in citrus because the long juvenile period of citrus presents challenges to identifying candidate genes that control this phenomenon. In this review, we focus on advances to our understanding of reproduction in citrus from the last decade and discuss priorities for the coming decade.

Comparative transcriptome analysis during early fruit development between three seedy citrus genotypes and their seedless mutants

Identification of genes with differential transcript abundance (GDTA) in seedless mutants may enhance understanding of seedless citrus development. Transcriptome analysis was conducted at three time points during early fruit development (Phase 1) of three seedy citrus genotypes: Fallglo (Bower citrus hybrid (Citrus reticulata×C. reticulata×C. paradisi)×Temple (C. reticulata×C. sinensis)), grapefruit (C. paradisi), Pineapple sweet orange (C. sinensis), and their seedless mutants. Seed abortion in seedless mutants was observed at 26 days post anthesis (Time point 2). Affymetrix transcriptomic analysis revealed 359 to 1077 probe sets with differential transcript abundance in the comparison of seedless versus seedy fruits for each citrus genotypes and time points. The GDTA identified by 18 microarray probe sets were validated by qPCR. Hierarchical clustering analysis revealed a range of GDTA associated with development, hormone and protein metabolism, all of which may reflect genes associated with seedless fruit development. There were 14, 9 and 12 genes found exhibiting similar abundance ratios in all three seedless versus seedy genotype comparisons at time point 1, 2 and 3, respectively. Among those genes were genes coding for an aspartic protease and a cysteine protease, which may play important roles in seedless fruit development. New insights into seedless citrus fruit development may contribute to biotech approaches to create seedless cultivars.