Phylogenetic tree-informed microRNAome analysis uncovers conserved and lineage-specific miRNAs in Camellia during floral organ development

Small RNA profiling reveals that an ovule-specific microRNA, cja-miR5179, targets a B-class MADS-box gene in Camellia japonica

BACKGROUND AND AIMS

The functional specialization of microRNA and its target genes is often an important factor in the establishment of spatiotemporal patterns of gene expression that are essential to plant development and growth. In different plant lineages, understanding the functional conservation and divergence of microRNAs remains to be explored.

METHODS

To identify small regulatory RNAs underlying floral patterning, we performed a tissue-specific profiling of small RNAs in various floral organs from single and double flower varieties (flowers characterized by multiple layers of petals) in Camellia japonica. We identified cja-miR5179, which belongs to a deeply conserved microRNA family that is conserved between angiosperms and basal plants but frequently lost in eudicots. We characterized the molecular function of cja-miR5179 and its target - a B-function MADS-box gene - through gene expression analysis and transient expression assays.

KEY RESULTS

We showed that cja-miR5179 is exclusively expressed in ovule tissues at the early stage of floral development. We found that cja-miR5179 targets the coding sequences of a DEFICIENS-like B-class gene (CjDEF) mRNA, which is located in the K motif of the MADS-box domain; and the target sites of miR5179/MADS-box were consistent in Camellia and orchids. Furthermore, through a petal transient-expression assay, we showed that the BASIC PENTACYSTEINE proteins bind to the GA-rich motifs in the cja-miR5179 promoter region and suppresses its expression.

CONCLUSIONS

We propose that the regulation between miR5179 and a B-class MADS-box gene in C. japonica has a deep evolutionary origin before the separation of monocots and dicots. During floral development of C. japonica, cja-miR5179 is specifically expressed in the ovule, which may be required for the inhibition of CjDEF function. This work highlights the evolutionary conservation as well as functional divergence of small RNAs in floral development.

Functional role of microRNA in the regulation of biotic and abiotic stress in agronomic plants

The increasing demand for food is the result of an increasing population. It is crucial to enhance crop yield for sustainable production. Recently, microRNAs (miRNAs) have gained importance because of their involvement in crop productivity by regulating gene transcription in numerous biological processes, such as growth, development and abiotic and biotic stresses. miRNAs are small, non-coding RNA involved in numerous other biological functions in a plant that range from genomic integrity, metabolism, growth, and development to environmental stress response, which collectively influence the agronomic traits of the crop species. Additionally, miRNA families associated with various agronomic properties are conserved across diverse plant species. The miRNA adaptive responses enhance the plants to survive environmental stresses, such as drought, salinity, cold, and heat conditions, as well as biotic stresses, such as pathogens and insect pests. Thus, understanding the detailed mechanism of the potential response of miRNAs during stress response is necessary to promote the agronomic traits of crops. In this review, we updated the details of the functional aspects of miRNAs as potential regulators of various stress-related responses in agronomic plants.

Transgressive and parental dominant gene expression and cytosine methylation during seed development in Brassica napus hybrids

KEY MESSAGE

Transcriptomic and epigenomic profiling of gene expression and small RNAs during seed and seedling development reveals expression and methylation dominance levels with implications on early stage heterosis in oilseed rape. The enhanced performance of hybrids through heterosis remains a key aspect in plant breeding; however, the underlying mechanisms are still not fully elucidated. To investigate the potential role of transcriptomic and epigenomic patterns in early expression of hybrid vigor, we investigated gene expression, small RNA abundance and genome-wide methylation in hybrids from two distant Brassica napus ecotypes during seed and seedling developmental stages using next-generation sequencing. A total of 31117, 344, 36229 and 7399 differentially expressed genes, microRNAs, small interfering RNAs and differentially methylated regions were identified, respectively. Approximately 70% of the differentially expressed or methylated features displayed parental dominance levels where the hybrid followed the same patterns as the parents. Via gene ontology enrichment and microRNA-target association analyses during seed development, we found copies of reproductive, developmental and meiotic genes with transgressive and paternal dominance patterns. Interestingly, maternal dominance was more prominent in hypermethylated and downregulated features during seed formation, contrasting to the general maternal gamete demethylation reported during gametogenesis in angiosperms. Associations between methylation and gene expression allowed identification of putative epialleles with diverse pivotal biological functions during seed formation. Furthermore, most differentially methylated regions, differentially expressed siRNAs and transposable elements were in regions that flanked genes without differential expression. This suggests that differential expression and methylation of epigenomic features may help maintain expression of pivotal genes in a hybrid context. Differential expression and methylation patterns during seed formation in an F₁ hybrid provide novel insights into genes and mechanisms with potential roles in early heterosis.

Genome-wide identification and development of miniature inverted-repeat transposable elements and intron length polymorphic markers in tea plant (Camellia sinensis)

Marker-assisted breeding and tagging of important quantitative trait loci for beneficial traits are two important strategies for the genetic improvement of plants. However, the scarcity of diverse and informative genetic markers covering the entire tea genome limits our ability to achieve such goals. In the present study, we used a comparative genomic approach to mine the tea genomes of Camellia sinensis var. assamica (CSA) and C. sinensis var. sinensis (CSS) to identify the markers to differentiate tea genotypes. In our study, 43 and 60 Camellia sinensis miniature inverted-repeat transposable element (CsMITE) families were identified in these two sequenced tea genomes, with 23,170 and 37,958 putative CsMITE sequences, respectively. In addition, we identified 4912 non-redundant, Camellia sinensis intron length polymorphic (CsILP) markers, 85.8% of which were shared by both the CSS and CSA genomes. To validate, a subset of randomly chosen 10 CsMITE markers and 15 CsILP markers were tested and found to be polymorphic among the 36 highly diverse tea genotypes. These genome-wide markers, which were identified for the first time in tea plants, will be a valuable resource for genetic diversity analysis as well as marker-assisted breeding of tea genotypes for quality improvement.

Transcriptome and Degradome Profiling Reveals a Role of miR530 in the Circadian Regulation of Gene Expression in Kalanchoë marnieriana

Crassulacean acid metabolism (CAM) is an important photosynthetic pathway for plant adaptation to dry environments. CAM plants feature a coordinated interaction between mesophyll and epidermis functions that involves refined regulations of gene expression. Plant microRNAs (miRNAs) are crucial post-transcription regulators of gene expression, however, their roles underlying the CAM pathway remain poorly investigated. Here, we present a study characterizing the expression of miRNAs in an obligate CAM species Kalanchoë marnieriana. Through sequencing of transcriptome and degradome in mesophyll and epidermal tissues under the drought treatments, we identified differentially expressed miRNAs that were potentially involved in the regulation of CAM. In total, we obtained 84 miRNA genes, and eight of them were determined to be Kalanchoë-specific miRNAs. It is widely accepted that CAM pathway is regulated by circadian clock. We showed that miR530 was substantially downregulated in epidermal peels under drought conditions; miR530 targeted two tandem zinc knuckle/PLU3 domain encoding genes (TZPs) that were potentially involved in light signaling and circadian clock pathways. Our work suggests that the miR530-TZPs module might play a role of regulating CAM-related gene expression in Kalanchoë.

Relationship between Secondary Metabolism and miRNA for Important Flavor Compounds in Different Tissues of Tea Plant (Camellia sinensis) As Revealed by Genome-Wide miRNA Analysis

This study investigated the regulatory relationship between important flavor compounds and microRNA (miRNA) in nine different tissues of tea plant by analyzing the related metabolites, small RNAs (sRNAs), degradome, and coexpression network. A total of 272 differential expressed miRNAs (DEmiRNAs) were obtained, including 198 conserved miRNAs and 74 novel miRNAs. Meanwhile, the expression patterns of miR159-GAMYB, miR167-ARF, and miR396-GRF pairs were investigated by quantitative real-time polymerase chain reaction (qRT-PCR) and the target sites were verified by 5'RNA ligase-mediated RACE (5' RLM-RACE). Further coexpression analysis showed that the content of gallated catechins was significantly and negatively correlated with the expression of miR156, but positively correlated with the expression of miR166 and miR172. Additionally, the expression of miR169a, miR169l, and miR319h was shown to be positively correlated with the content of nongallated catechins and the experssion levels of ANRa, ANRb, and LARb. Moreover, important volatile compounds, such as linalool, geraniol, and 2-phenylethanol, were found to be highly positively correlated with the expression of miR171o, miRN71a, miRN71b, miRN71c, and miRN71d. Our data indicate that these miRNAs may play important roles in regulating the biosynthesis of flavor compounds in different tissues of tea plant.

Independent recurrent evolution of MICRORNA genes converging onto similar non-canonical organisation across green plant lineages is driven by local and segmental duplication events in species, family and lineages

The relationship between evolutionary history, organisation and transcriptional regulation of genes are intrinsically linked. These have been well studied in canonically organised protein-coding genes but not of MIRNAs. In the present study, we investigated the non-canonical arrangement of MIRNAs across taxonomic boundaries from algae to angiosperms employing a combination of genome organization, phylogeny and synteny. We retrieved the complete dataset of MIRNA from twenty-five species to identify and classify based on organisational patterns. The median size of cluster was between 2-5 kb and between 1-20 % of all MIRNAs are organized in head-to-head (with bidirectional promoter), head-to-tail (tandem), and overlapping manner. Although majority of the clusters are composed of MIRNA homologs, 25% of all clusters comprises of non-homologous genes with a potential of generating functional and regulatory complexity. A comparison of phylogeny and organizational patterns revealed that multiple independent events, some of which are species-specific, and some ancient, in different lineages, are responsible for non-canonical organization. Detailed investigation of MIR395 family across the plants revealed a complex origin of non-canonical arrangement through ancient and recent, segmental and local duplications; analysis of MIR399 family revealed major expansion occurred prior to monocot-dicot split, with few lineage-specific events. Evolution of "convergent" organization pattern of non-canonical arrangement originating from independent loci through recurrent event highlights our poor understanding of evolutionary process of MIRNA genes. The present investigation thus paves way for comparative functional genomics to understand the role of non-canonical organization on transcriptional regulation and regulatory diversity in MIRNA gene families.

Anti-NMDA Receptor Encephalitis, Vaccination and Virus

Anti-N-methyl-d-aspartate (Anti-NMDA) receptor encephalitis is an acute autoimmune disorder. The symptoms range from psychiatric symptoms, movement disorders, cognitive impairment, and autonomic dysfunction. Previous studies revealed that vaccination might induce this disease. A few cases were reported to be related to H1N1 vaccine, tetanus/diphtheria/pertussis and polio vaccine, and Japanese encephalitis vaccine. Although vaccination is a useful strategy to prevent infectious diseases, in a low risk, it may trigger serious neurological symptoms. In addition to anti-NMDA receptor encephalitis, other neurological diseases were reported to be associated with a number of vaccines. In this paper, the anti-NMDA receptor encephalitis cases related to a number of vaccines and other neurological symptoms that might be induced by these vaccines were reviewed. In addition, anti-NMDA receptor encephalitis cases that were induced by virus infection were also reviewed.

Dual-mode ultrasound radiomics and intrinsic imaging phenotypes for diagnosis of lymph node lesions

Background

The ultrasonic diagnosis of lymph node lesions is usually based on a small number of subjective visual features from a single ultrasonic modality, which limits diagnostic accuracy. Therefore, our study aimed to propose a computerized method for using dual-mode ultrasound radiomics and the intrinsic imaging phenotypes for accurately differentiating benign, lymphomatous, and metastatic lymph nodes.

Methods

A total of 543 lymph nodes from 538 patients were examined with both B-mode ultrasonography and elastography. The data set was randomly divided into a training set of 407 nodes and a validation set of 136 nodes. First, we extracted 430 radiomic features from dual-mode images. Then, we combined the least absolute shrinkage and selection operator with the analysis of variance to select several typical features. We retrieved the intrinsic imaging phenotypes by using a hierarchical clustering of all radiomics features, and we integrated the phenotypes with the selected features for the classification of benign, lymphomatous, and metastatic nodes.

Results

The areas under the receiver operating characteristic curves (AUCs) on the validation set were 0.960 for benign vs. lymphomatous, 0.716 for benign vs. metastatic, 0.933 for lymphomatous vs. metastatic, and 0.856 for benign vs. malignant.

Conclusions

The radiomics features and intrinsic imaging phenotypes derived from the dual-mode ultrasound can capture the distinctions between benign, lymphomatous, and metastatic nodes and are valuable in node differentiation.

Identification of vital candidate microRNA/mRNA pairs regulating ovule development using high-throughput sequencing in hazel

BACKGROUND

Hazels (Corylus spp.) are economically important nut-producing species in which ovule development determines seed plumpness, one of the key parameters reflecting nut quality. microRNAs (miRNAs) play important roles in RNA silencing and the post-transcriptional regulation of gene expression. However, very little is currently known regarding the miRNAs involved in regulating ovule growth and development.

RESULTS

In this study, we accordingly sought to determine the important miRNAs involved in ovule development and growth in hazel. We examined ovules at four developmental stages, namely, ovule formation (Ov1), early ovule growth (Ov2), rapid ovule growth (Ov3), and ovule maturity (Ov4). On the basis of small RNA and mRNA sequencing using the Illumina sequencing platform, we identified 970 miRNAs in hazel, of which 766 and 204 were known and novel miRNAs, respectively. In Ov1-vs-Ov2, Ov1-vs-Ov3, Ov1-vs-Ov4, Ov2-vs-Ov3, Ov2-vs-Ov4, and Ov3-vs-Ov4 paired comparisons, 471 differentially expressed microRNAs (DEmiRNAs) and their 3117 target differentially expressed messenger RNAs (DEmRNAs) formed 11,199 DEmiRNA/DEmRNA pairs, with each DEmiRNA changing the expression of an average of 6.62 target mRNAs. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of all DEmRNAs revealed 29 significantly enriched KEGG pathways in the six paired comparisons, including protein export (ko03060), fatty acid elongation (ko00062), starch and sucrose metabolism (ko00500), fatty acid biosynthesis (ko00061), and amino sugar and nucleotide sugar metabolism (ko00520). Our results indicate that DEmiRNA/DEmRNA pairs showing opposite change trends were related to stress tolerance, embryo and seed development, cell proliferation, auxin transduction, and the biosynthesis of proteins, starch, and fats may participate in ovule growth and development.

CONCLUSIONS

These findings contribute to a better understanding of ovule development at the level of post-transcriptional regulation, and lay the foundation for further functional analyses of hazelnut ovule growth and development.

Identification of alternatively spliced gene isoforms and novel noncoding RNAs by single-molecule long-read sequencing in Camellia

Direct single-molecule sequencing of full-length transcripts allows efficient identification of gene isoforms, which is apt to alternative splicing (AS), polyadenylation, and long non-coding RNA analyses. However, the identification of gene isoforms and long non-coding RNAs with novel regulatory functions remains challenging, especially for species without a reference genome. Here, we present a comprehensive analysis of a combined long-read and short-read transcriptome sequencing in Camellia japonica. Through a novel bioinformatic pipeline of reverse-tracing the split-sites, we have uncovered 257,692 AS sites from 61,838 transcripts; and 13,068 AS isoforms have been validated by aligning the short reads. We have identified the tissue-specific AS isoforms along with 6,373 AS events that were found in all tissues. Furthermore, we have analysed the polyadenylation (polyA) patterns of transcripts, and found that the preference for polyA signals was different between the AS and non-AS transcripts. Moreover, we have predicted the phased small interfering RNA (phasiRNA) loci through integrative analyses of transcriptome and small RNA sequencing. We have shown that a newly evolved phasiRNA locus from lipoxygenases generated 12 consecutive 21 bp secondary RNAs, which were responsive to cold and heat stress in Camellia. Our studies of the isoform transcriptome provide insights into gene splicing and functions that may facilitate the mechanistic understanding of plants.

Integrative expression network analysis of microRNA and gene isoforms in sacred lotus

Background

Gene expression is complex and regulated by multiple molecular mechanisms, such as miRNA-mediated gene inhibition and alternative-splicing of pre-mRNAs. However, the coordination of interaction between miRNAs with different splicing isoforms, and the change of splicing isoform in response to different cellular environments are largely unexplored in plants. In this study, we analyzed the miRNA and mRNA transcriptome from lotus (Nelumbo nucifera), an economically important flowering plant.

Results

Through RNA-seq analyses on miRNAs and their target genes (isoforms) among six lotus tissues, expression of most miRNAs seem to be negatively correlated with their targets and tend to be tissue-specific. Further, our results showed that preferential interactions between miRNAs and hub gene isoforms in one coexpression module which is highly correlated with leaf. Intriguingly, for many genes, their corresponding isoforms were assigned to different co-expressed modules, and they exhibited more divergent mRNA structures including presence and absence of miRNA binding sites, suggesting functional divergence for many isoforms is escalated by both structural and expression divergence. Further detailed functional enrichment analysis of miRNA targets revealed that miRNAs are involved in the regulation of lotus growth and development by regulating plant hormone-related pathway genes.

Conclusions

Taken together, our comprehensive analyses of miRNA and mRNA transcriptome elucidate the coordination of interaction between miRNAs and different splicing isoforms, and highlight the functional divergence of many transcript isoforms from the same locus in lotus.

Genome-Wide microRNA Profiling Using Oligonucleotide Microarray Reveals Regulatory Networks of microRNAs in Nicotiana benthamiana During Beet Necrotic Yellow Vein Virus Infection

Beet necrotic yellow vein virus (BNYVV) infections induce stunting and leaf curling, as well as root and floral developmental defects and leaf senescence in Nicotiana benthamiana. A microarray analysis with probes capable of detecting 1596 candidate microRNAs (miRNAs) was conducted to investigate differentially expressed miRNAs and their targets upon BNYVV infection of N. benthamiana plants. Eight species-specific miRNAs of N. benthamiana were identified. Comprehensive characterization of the N. benthamiana microRNA profile in response to the BNYVV infection revealed that 129 miRNAs were altered, including four species-specific miRNAs. The targets of the differentially expressed miRNAs were predicted accordingly. The expressions of miR164, 160, and 393 were up-regulated by BNYVV infection, and those of their target genes, NAC21/22, ARF17/18, and TIR, were down-regulated. GRF1, which is a target of miR396, was also down-regulated. Further genetic analysis of GRF1, by Tobacco rattle virus-induced gene silencing, assay confirmed the involvement of GRF1 in the symptom development during BNYVV infection. BNYVV infection also induced the up-regulation of miR168 and miR398. The miR398 was predicted to target umecyanin, and silencing of umecyanin could enhance plant resistance against viruses, suggesting the activation of primary defense response to BNYVV infection in N. benthamiana. These results provide a global profile of miRNA changes induced by BNYVV infection and enhance our understanding of the mechanisms underlying BNYVV pathogenesis.

Identification of Populus Small RNAs Responsive to Mutualistic Interactions With Mycorrhizal Fungi, Laccaria bicolor and Rhizophagus irregularis

Ecto- and endo-mycorrhizal colonization of Populus roots have a positive impact on the overall tree health and growth. A complete molecular understanding of these interactions will have important implications for increasing agricultural or forestry sustainability using plant:microbe-based strategies. These beneficial associations entail extensive morphological changes orchestrated by the genetic reprogramming in both organisms. In this study, we performed a comparative analysis of two Populus species (Populus deltoides and P. trichocarpa) that were colonized by either an arbuscular mycorrhizal fungus (AmF), Rhizophagus irregularis or an ectomycorrhizal fungus (EmF), Laccaria bicolor, to describe the small RNA (sRNA) landscape including small open reading frames (sORFs) and micro RNAs (miRNAs) involved in these mutualistic interactions. We identified differential expression of sRNAs that were, to a large extent, (1) within the genomic regions lacking annotated genes in the Populus genome and (2) distinct for each fungal interaction. These sRNAs may be a source of novel sORFs within a genome, and in this regard, we identified potential sORFs encoded by the sRNAs. We predicted a higher number of differentially-expressed miRNAs in P. trichocarpa (4 times more) than in P. deltoides (conserved and novel). In addition, 44 miRNAs were common in P. trichocarpa between the EmF and AmF treatments, and only 4 miRNAs were common in P. deltoides between the treatments. Root colonization by either fungus was more effective in P. trichocarpa than in P. deltoides, thus the relatively few differentially-expressed miRNAs predicted in P. deltoides might reflect the extent of the symbiosis. Finally, we predicted several genes targets for the plant miRNAs identified here, including potential fungal gene targets. Our findings shed light on additional molecular tiers with a role in Populus-fungal mutualistic associations and provides a set of potential molecular targets for future enhancement.

Identifying High Confidence microRNAs in the Developing Seeds of Jatropha curcas

MicroRNAs (miRNAs) are endogenously short noncoding regulatory RNAs implicated in plant development and physiology. Nine small RNA (sRNA) libraries from three typical seed developmental stages (young, intermediate, and mature) were generated by deep sequencing to identify the miRNAs of J. curcas, a potential oilseed crop for the production of renewable oil. Strict criteria were adopted to identify 93 high confidence miRNAs including 48 conserved miRNAs and 45 novel miRNAs. Target genes of these miRNAs were involved in a broad range of physiological functions, including gene expression regulation, primary & secondary metabolism, growth & development, signal transduction, and stress response. About one third (29 out of 93) miRNAs showed significant changes in expression levels during the seed developmental process, indicating that the miRNAs might regulate its targets by their changes of transcription levels in seed development. However, most miRNAs were found differentially expressed in the late stage of seed development, suggesting that miRNAs play more important roles in the stage when seed accumulating organic matters and suffering dehydration stress. This study presents the first large scale identification of high confidence miRNAs in the developing seeds of J. curcas.

Unraveling the Roles of Regulatory Genes during Domestication of Cultivated Camellia: Evidence and Insights from Comparative and Evolutionary Genomics

With the increasing power of DNA sequencing, the genomics-based approach is becoming a promising resolution to dissect the molecular mechanism of domestication of complex traits in trees. Genus Camellia possesses rich resources with a substantial value for producing beverage, ornaments, edible oil and more. Currently, a vast number of genetic and genomic research studies in Camellia plants have emerged and provided an unprecedented opportunity to expedite the molecular breeding program. In this paper, we summarize the recent advances of gene expression and genomic resources in Camellia species and focus on identifying genes related to key economic traits such as flower and fruit development and stress tolerances. We investigate the genetic alterations and genomic impacts under different selection programs in closely related species. We discuss future directions of integrating large-scale population and quantitative genetics and multiple omics to identify key candidates to accelerate the breeding process. We propose that future work of exploiting the genomic data can provide insights related to the targets of domestication during breeding and the evolution of natural trait adaptations in genus Camellia.

Comparison of multi-lineage differentiation of hiPSCs reveals novel miRNAs that regulate lineage specification

MicroRNAs (miRNAs) are known to be crucial players in governing the differentiation of human induced pluripotent stem cells (hiPSCs). Despite their utter importance, identifying key lineage specifiers among the myriads of expressed miRNAs remains challenging. We believe that the current practice in mining miRNA specifiers via delineating dynamic fold-changes only is inadequate. Our study, therefore, provides evidence to pronounce "lineage specificity" as another important attribute to qualify for these lineage specifiers. Adopted hiPSCs were differentiated into representative lineages (hepatic, nephric and neuronal) over all three germ layers whilst the depicted miRNA expression changes compiled into an integrated atlas. We demonstrated inter-lineage analysis shall aid in the identification of key miRNAs with lineage-specificity, while these shortlisted candidates were collectively known as "lineage-specific miRNAs". Subsequently, we followed through the fold-changes along differentiation via computational analysis to identify miR-192 and miR-372-3p, respectively, as representative candidate key miRNAs for the hepatic and nephric lineages. Indeed, functional characterization validated that miR-192 and miR-372-3p regulate lineage differentiation via modulation of the expressions of lineage-specific genes. In summary, our presented miRNA atlas is a resourceful ore for the mining of key miRNAs responsible for lineage specification.

Comparative genomics analysis reveals gene family expansion and changes of expression patterns associated with natural adaptations of flowering time and secondary metabolism in yellow Camellia

Yellow-flowering species are unique in the genus Camellia not only for their bright yellow pigments but also the health-improving substances in petals. However, little is known regarding the biosynthesis pathways of pigments and secondary metabolites. Here, we performed comparative genomics studies in two yellow-flowered species of the genus Camellia with distinctive flowering periods. We obtained 112,190 and 89,609 unigenes from Camellia nitidissima and Camellia chuongtsoensis, respectively, and identified 9547 gene family clusters shared with various plant species and 3414 single-copy gene families. Global gene expression analysis revealed six comparisons of differentially expressed gene sets in different developmental stages of floral bud. Through the identification of orthologous pairs, conserved and specific differentially expressed genes (DEGs) between species were compared. Functional enrichment analysis suggested that the gibberellin (GA) biosynthesis pathway might be related to the alteration of flowering responses. Furthermore, the expression patterns of secondary metabolism pathway genes were analyzed between yellow- and red-flowered Camellias. We showed that the key enzymes involved in glycosylation of flavonoids displayed differential expression patterns, indicating that the direct glycosylation of flavonols rather than anthocyanins was pivotal to coloration and health-improving metabolites in the yellow Camellia petals. Finally, the gene family analysis of UDP-glycosyltransferases revealed an expansion of group C members in C. nitidissima. Through comparative genomics analysis, we demonstrate that changes of gene expression and gene family members are critical to the variation of natural traits. This work provides valuable insights into the molecular regulation of trait adaptations of floral pigmentation and flowering timing.

Characterization of microRNAs of Beta macrocarpa and their responses to Beet necrotic yellow vein virus infection

Plant microRNAs (miRNAs) are a class of non-coding RNAs that play important roles in plant development, defense, and symptom development. Here, 547 known miRNAs representing 129 miRNA families, and 282 potential novel miRNAs were identified in Beta macrocarpa using small RNA deep sequencing. A phylogenetic analysis was performed, and 8 Beta lineage-specific miRNAs were identified. Through a differential expression analysis, miRNAs associated with Beet necrotic yellow vein virus (BNYVV) infection were identified and confirmed using a microarray analysis and stem-loop RT-qPCR. In total, 103 known miRNAs representing 38 miRNA families, and 45 potential novel miRNAs were differentially regulated, with at least a two-fold change, in BNYVV-infected plants compared with that of the mock-inoculated control. Targets of these differentially expressed miRNAs were also predicted by degradome sequencing. These differentially expressed miRNAs were involved in hormone biosynthesis and signal transduction pathways, and enhanced axillary bud development and plant defenses. This work is the first to describe miRNAs of the plant genus Beta and may offer a reference for miRNA research in other species in the genus. It provides valuable information on the pathogenicity mechanisms of BNYVV.

Global gene expression defines faded whorl specification of double flower domestication in Camellia

Double flowers in cultivated camellias are divergent in floral patterns which present a rich resource for demonstrating molecular modifications influenced by the human demands. Despite the key principle of ABCE model in whorl specification, the underlying mechanism of fine-tuning double flower formation remains largely unclear. Here a comprehensive comparative transcriptomics interrogation of gene expression among floral organs of wild type and “formal double” and “anemone double” is presented. Through a combination of transcriptome, small RNA and “degradome” sequencing, we studied the regulatory gene expression network underlying the double flower formation. We obtained the differentially expressed genes between whorls in wild and cultivated Camellia. We showed that the formation of double flowers tends to demolish gene expression canalization of key functions; the faded whorl specification mechanism was fundamental under the diverse patterns of double flowers. Furthermore, we identified conserved miRNA-targets regulations in the control of double flowers, and we found that miR172-AP2, miR156-SPLs were critical regulatory nodes contributing to the diversity of double flower forms. This work highlights the hierarchical patterning of global gene expression in floral development, and supports the roles of “faded ABC model” mechanism and miRNA-targets regulations underlying the double flower domestication.