The gut microbiota plays a crucial role in the bi-directional gut–brain axis, a communication that integrates the gut and central nervous system (CNS) activities. Animal studies reveal that gut bacteria influence behaviour, Brain-Derived Neurotrophic Factor (BDNF) levels and serotonin metabolism. In the present study, we report for the first time an analysis of the microbiota–gut–brain axis in zebrafish (Danio rerio). After 28 days of dietary administration with the probiotic Lactobacillus rhamnosus IMC 501, we found differences in shoaling behaviour, brain expression levels of bdnf and of genes involved in serotonin signalling/metabolism between control and treated zebrafish group. In addition, in microbiota we found a significant increase of Firmicutes and a trending reduction of Proteobacteria. This study demonstrates that selected microbes can be used to modulate endogenous neuroactive molecules in zebrafish.

Bilateral symmetry of flowers is a relevant novelty that has occurred many times throughout the evolution of flowering plants. In Antirrhinum majus, establishment of flower dorso-ventral asymmetry is mainly due to interaction of TCP (CYC and DICH) and MYB (DIV, RAD, and DRIF) transcription factors. In the present study, we characterized 8 DIV-, 4 RAD-, and 2 DRIF-like genes from the transcriptome of Orchis italica, an orchid species with bilaterally symmetric and resupinate flowers. We found a similar number of DIV- and RAD-like genes within the genomes of Phalaenopsis equestris and Dendrobium catenatum orchids. Orchid DIV- and RAD-like proteins share conserved motifs whose distribution reflects their phylogeny and analysis of the genomic organization revealed a single intron containing many traces of transposable elements. Evolutionary analysis has shown that purifying selection acts on the DIV- and RAD-like coding regions in orchids, with relaxation of selective constraints in a branch of the DIV-like genes. Analysis of the expression patterns of DIV- and RAD-like genes in O. italica revealed possible redundant functions for some of them. In the perianth of O. italica, the ortholog of DIV and DRIF of A. majus are expressed in all tissues, whereas RAD is mainly expressed in the outer tepals and lip. These data allow for proposal of an evolutionary conserved model in which the expression of the orthologs of the DIV, RAD, and DRIF genes might be related to establishment of flower bilateral symmetry in the nonmodel orchid species O. italica.

Plant microRNAs (miRNAs) are small, regulatory non-coding RNAs involved in a wide range of biological processes, from organ development to response to stimuli. In recent years, an increasing number of studies on model plant species have highlighted the evolutionary conservation of a high number of miRNA families and the existence of taxon-specific ones. However, few studies have examined miRNAs in non-model species such as orchids, which are characterized by highly diversified floral structures and pollination strategies. Therefore, we analysed a small RNA library of inflorescence tissue of the Mediterranean orchid Orchis italica to increase the knowledge on miRNAs in a non-model plant species. The high-throughput sequencing and analysis of a small RNA library of inflorescence of O. italica revealed 23 conserved and 161 putative novel miRNA families. Among the putative miRNA targets, experimental validation demonstrated that a DEF-like MADS-box transcript is cleaved by the homolog of miR5179 of O. italica. The presence of conserved miRNA families in the inflorescence of O. italica indicates that the basic developmental flower regulatory mechanisms mediated by miRNAs are maintained through evolution. Because, according to the “orchid code” theory, DEF-like genes exert a key function in the diversification of tepals and lip, the cleavage-mediated inhibitory activity of miR5179 on a OitaDEF-like transcript suggests that, in orchids, miRNAs play an important role in the diversification of the perianth organs.

TCP proteins are plant-specific transcription factors involved in many different processes. Because of their involvement in a large number of developmental pathways, their roles have been investigated in various plant species. However, there are almost no studies of this transcription factor family in orchids. Based on the available transcriptome of the inflorescence of the orchid Orchis italica, in the present study we identified 12 transcripts encoding TCP proteins. The phylogenetic analysis showed that they belong to different TCP classes (I and II) and groups (PCF, CIN and CYC/TB1), and that they display a number of conserved motifs when compared with the TCPs of Arabidopsis and Oryza. The presence of a specific cleavage site for the microRNA miRNA319, an important post-transcriptional regulator of several TCP genes in other species, was demonstrated for one transcript of O. italica, and the analysis of the expression pattern of the TCP transcripts in different inflorescence organs and in leaf tissue suggests that some TCP transcripts of O. italica exert their role only in specific tissues, while others may play multiple roles in different tissues. In addition, the evolutionary analysis showed a general purifying selection acting on the coding region of these transcripts.

The floral transcriptome of Orchis italica, a wild orchid species, was obtained using Illumina RNA-seq technology and specific de novo assembly and analysis tools. More than 100 million raw reads were processed resulting in 132,565 assembled transcripts and 86,079 unigenes with an average length of 606 bp and N50 of 956 bp. Functional annotation assigned 38,984 of the unigenes to records present in the NCBI non-redundant protein database, 32,161 of them to Gene Ontology terms, 15,775 of them to Eukaryotic Orthologous Groups (KOG) and 7,143 of them to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The in silico expression analysis based on the Fragments Per Kilobase of transcript per Million mapped reads (FPKM) was confirmed by real-time RT-PCR experiments on 10 selected unigenes, which showed high and statistically significant positive correlation with the RNA-seq based expression data. The prediction of putative long non-coding RNAs was assessed using two different software packages, CPC and Portrait, resulting in 7,779 unannotated unigenes that matched the threshold values for both of the analyses. Among the predicted long non-coding RNAs, one is the homologue of TAS3, a long non-coding RNA precursor of trans-acting small interfering RNAs (ta-siRNAs). The differential expression pattern observed for the selected putative long non-coding RNAs suggests their possible functional role in different floral tissues.