Microbiota changes induced by microencapsulated sodium butyrate in patients with inflammatory bowel disease

BACKGROUND

Butyrate has shown anti-inflammatory and regenerative properties, providing symptomatic relief when orally supplemented in patients suffering from various colonic diseases. We investigated the effect of a colonic-delivery formulation of butyrate on the fecal microbiota of patients with inflammatory bowel diseases (IBDs).

METHODS

In this double-blind, placebo-controlled, pilot study, 49 IBD patients (n = 19 Crohn's disease, CD and n = 30 ulcerative colitis, UC) were randomized to oral administration of microencapsulated-sodium-butyrate (BLM) or placebo for 2 months, in addition to conventional therapy. Eighteen healthy volunteers (HVs) were recruited to provide a healthy microbiota model of the local people. Fecal microbiota from stool samples was assessed by 16S sequencing. Clinical disease activity and quality of life (QoL) were evaluated before and after treatment.

KEY RESULTS

At baseline, HVs showed a different microbiota composition compared with IBD patients. Sodium-butyrate altered the gut microbiota of IBD patients by increasing bacteria able to produce SCFA in UC patients (Lachnospiraceae spp.) and the butyrogenic colonic bacteria in CD patients (Butyricicoccus). In UC patients, QoL was positively affected by treatment.

CONCLUSIONS AND INFERENCES

Sodium-butyrate supplementation increases the growth of bacteria able to produce SCFA with potentially anti-inflammatory action. The clinical impact of this finding requires further investigation.

Assessment of statistical methods from single cell, bulk RNA-seq, and metagenomics applied to microbiome data

BACKGROUND

The correct identification of differentially abundant microbial taxa between experimental conditions is a methodological and computational challenge. Recent work has produced methods to deal with the high sparsity and compositionality characteristic of microbiome data, but independent benchmarks comparing these to alternatives developed for RNA-seq data analysis are lacking.

RESULTS

We compare methods developed for single-cell and bulk RNA-seq, and specifically for microbiome data, in terms of suitability of distributional assumptions, ability to control false discoveries, concordance, power, and correct identification of differentially abundant genera. We benchmark these methods using 100 manually curated datasets from 16S and whole metagenome shotgun sequencing.

CONCLUSIONS

The multivariate and compositional methods developed specifically for microbiome analysis did not outperform univariate methods developed for differential expression analysis of RNA-seq data. We recommend a careful exploratory data analysis prior to application of any inferential model and we present a framework to help scientists make an informed choice of analysis methods in a dataset-specific manner.

A draft genome of sweet cherry (Prunus avium L.) reveals genome-wide and local effects of domestication

Sweet cherry (Prunus avium L.) trees are both economically important fruit crops but also important components of natural forest ecosystems in Europe, Asia and Africa. Wild and domesticated trees currently coexist in the same geographic areas with important questions arising on their historical relationships. Little is known about the effects of the domestication process on the evolution of the sweet cherry genome. We assembled and annotated the genome of the cultivated variety "Big Star*" and assessed the genetic diversity among 97 sweet cherry accessions representing three different stages in the domestication and breeding process (wild trees, landraces and modern varieties). The genetic diversity analysis revealed significant genome-wide losses of variation among the three stages and supports a clear distinction between wild and domesticated trees, with only limited gene flow being detected between wild trees and domesticated landraces. We identified 11 domestication sweeps and five breeding sweeps covering, respectively, 11.0 and 2.4 Mb of the P. avium genome. A considerable fraction of the domestication sweeps overlaps with those detected in the related species, Prunus persica (peach), indicating that artificial selection during domestication may have acted independently on the same regions and genes in the two species. We detected 104 candidate genes in sweep regions involved in different processes, such as the determination of fruit texture, the regulation of flowering and fruit ripening and the resistance to pathogens. The signatures of selection identified will enable future evolutionary studies and provide a valuable resource for genetic improvement and conservation programs in sweet cherry.

Microbial metabarcoding highlights different bacterial and fungal populations in honey samples from local beekeepers and market in north-eastern Italy

Due to its chemical properties, honey does not foster the growth of microorganisms, however it may contain a rich microbial community, including viable, stressed, and not viable microbes. In order to characterize honey microbiota focusing on the difference between products from beekeepers and large retail in the present study a culture-independent approach based on DNA metabarcoding was applied. Honey samples were collected from Local Beekeepers (LB) and Market sales (M) during four years with the aim to investigate the microbiological quality in the honey market. Extraction and amplification of DNA from honey samples showed reduced efficiency with increasing age of honey, with the loss of 50-80% of samples four years old (2014). For this reason, only samples of similar age were compared and the analysis of microbial communities focused on year 2017, for a total of 75 samples. Differences in alpha and beta-diversity were evidenced comparing microbial communities between LB and M samples. In particular, contaminant bacteria dominated the microbiota in M samples while LB samples were enriched in Lactic Acid Bacteria (LAB) that cannot be isolated with culture-dependent approaches.

The Mitochondrial Genome Assembly of Fennel (Foeniculum vulgare) Reveals Two Different atp6 Gene Sequences in Cytoplasmic Male Sterile Accessions

Cytoplasmic male sterility (CMS) has always aroused interest among researchers and breeders, being a valuable resource widely exploited not only to breed F1 hybrid varieties but also to investigate genes that control stamen and pollen development. With the aim of identifying candidate genes for CMS in fennel, we adopted an effective strategy relying on the comparison between mitochondrial genomes (mtDNA) of both fertile and sterile genotypes. mtDNA raw reads derived from a CMS genotype were assembled in a single molecule (296,483 bp), while a draft mtDNA assembly (166,124 nucleotides, 94 contigs) was performed using male fertile sample (MF) sequences. From their annotation and alignment, two atp6-like sequences were identified. atp6-, the putative mutant copy with a 300 bp truncation at the 5'-end, was found only in the mtDNA of CMS samples, while the wild type copy (atp6+) was detected only in the MF mtDNA. Further analyses (i.e., reads mapping and Sanger sequencing), revealed an atp6+ copy also in CMS samples, probably in the nuclear DNA. However, qPCRs performed on different tissues proved that, despite its availability, atp6+ is expressed only in MF samples, while apt6- mRNA was always detected in CMS individuals. In the light of these findings, the energy deficiency model could explain the pollen deficiency observed in male sterile flower. atp6- could represent a gene whose mRNA is translated into a not-fully functional protein leading to suboptimal ATP production that guarantees essential cellular processes but not a high energy demand process such as pollen development. Our study provides novel insights into the fennel mtDNA genome and its atp6 genes, and paves the way for further studies aimed at understanding their functional roles in the determination of male sterility.

Engineering a 3D in vitro model of human skeletal muscle at the single fiber scale

The reproduction of reliable in vitro models of human skeletal muscle is made harder by the intrinsic 3D structural complexity of this tissue. Here we coupled engineered hydrogel with 3D structural cues and specific mechanical properties to derive human 3D muscle constructs (“myobundles”) at the scale of single fibers, by using primary myoblasts or myoblasts derived from embryonic stem cells. To this aim, cell culture was performed in confined, laminin-coated micrometric channels obtained inside a 3D hydrogel characterized by the optimal stiffness for skeletal muscle myogenesis. Primary myoblasts cultured in our 3D culture system were able to undergo myotube differentiation and maturation, as demonstrated by the proper expression and localization of key components of the sarcomere and sarcolemma. Such approach allowed the generation of human myobundles of ~10 mm in length and ~120 μm in diameter, showing spontaneous contraction 7 days after cell seeding. Transcriptome analyses showed higher similarity between 3D myobundles and skeletal signature, compared to that found between 2D myotubes and skeletal muscle, mainly resulting from expression in 3D myobundles of categories of genes involved in skeletal muscle maturation, including extracellular matrix organization. Moreover, imaging analyses confirmed that structured 3D culture system was conducive to differentiation/maturation also when using myoblasts derived from embryonic stem cells. In conclusion, our structured 3D model is a promising tool for modelling human skeletal muscle in healthy and diseases conditions.

A single polyploidization event at the origin of the tetraploid genome of Coffea arabica is responsible for the extremely low genetic variation in wild and cultivated germplasm

The genome of the allotetraploid species Coffea arabica L. was sequenced to assemble independently the two component subgenomes (putatively deriving from C. canephora and C. eugenioides) and to perform a genome-wide analysis of the genetic diversity in cultivated coffee germplasm and in wild populations growing in the center of origin of the species. We assembled a total length of 1.536 Gbp, 444 Mb and 527 Mb of which were assigned to the canephora and eugenioides subgenomes, respectively, and predicted 46,562 gene models, 21,254 and 22,888 of which were assigned to the canephora and to the eugeniodes subgenome, respectively. Through a genome-wide SNP genotyping of 736 C. arabica accessions, we analyzed the genetic diversity in the species and its relationship with geographic distribution and historical records. We observed a weak population structure due to low-frequency derived alleles and highly negative values of Taijma’s D, suggesting a recent and severe bottleneck, most likely resulting from a single event of polyploidization, not only for the cultivated germplasm but also for the entire species. This conclusion is strongly supported by forward simulations of mutation accumulation. However, PCA revealed a cline of genetic diversity reflecting a west-to-east geographical distribution from the center of origin in East Africa to the Arabian Peninsula. The extremely low levels of variation observed in the species, as a consequence of the polyploidization event, make the exploitation of diversity within the species for breeding purposes less interesting than in most crop species and stress the need for introgression of new variability from the diploid progenitors.

Draft genome assembly and transcriptome data of the icefish Chionodraco myersi reveal the key role of mitochondria for a life without hemoglobin at subzero temperatures

Antarctic fish belonging to Notothenioidei represent an extraordinary example of radiation in the cold. In addition to the absence of hemoglobin, icefish show a number of other striking peculiarities including large-diameter blood vessels, high vascular densities, mitochondria-rich muscle cells, and unusual mitochondrial architecture. In order to investigate the bases of icefish adaptation to the extreme Southern Ocean conditions we sequenced the complete genome of the icefish Chionodraco myersi. Comparative analyses of the icefish genome with those of other teleost species, including two additional white-blooded and five red-blooded notothenioids, provided a new perspective on the evolutionary loss of globin genes. Muscle transcriptome comparative analyses against red-blooded notothenioids as well as temperate fish revealed the peculiar regulation of genes involved in mitochondrial function in icefish. Gene duplication and promoter sequence divergence were identified as genome-wide patterns that likely contributed to the broad transcriptional program underlying the unique features of icefish mitochondria.

Inhibition of Virulence-Related Traits in Pseudomonas syringae pv. actinidiae by Gunpowder Green Tea Extracts

Green tea is a widely-consumed healthy drink produced from the leaves of Camellia sinensis. It is renowned for its antioxidant and anticarcinogenic properties, but also displays significant antimicrobial activity against numerous human pathogens. Here we analyzed the antimicrobial activity of Gunpowder green tea against Pseudomonas syringae pv. actinidiae (Psa), the agent that causes kiwifruit bacterial canker. At the phenotypic level, tea extracts strongly inhibited Psa growth and swimming motility, suggesting it could reduce Psa epiphytic survival during plant colonization. The loss of bacterial virulence-related traits following treatment with tea extracts was also investigated by large-scale transcriptome analysis, which confirmed the in vitro phenotypes and revealed the induction of adaptive responses in the treated bacteria allowing them to cope with iron deficiency and oxidative stress. Such molecular changes may account for the ability of Gunpowder green tea to protect kiwifruit against Psa infection.

Bark and Grape Microbiome of Vitis vinifera: Influence of Geographic Patterns and Agronomic Management on Bacterial Diversity

In recent years, the concept of “microbial terroir” has been introduced in the frame of the more renowned notion of “vitivinicultural terroir,’ since several studies demonstrated that wine characteristics are related to regional microbial community compositions. Most of the existing research focused on grape berries microbiota, since it can directly impact wine quality. In this work we studied, for the first time through next-generation sequencing, the epiphytic bacterial community of vine bark and its relationships with grape microbiota. The study was carried out in two Italian wine appellations (situated in different regions) to explore the impact of biogeography, and the influence of two agronomical practices (biodynamic and conventional) was evaluated as well. Overall, our results show that grapevine bark harbors a rich epiphytic microbiota and displays a higher microbial biodiversity than grape berry. Moreover, this study suggests that geographic and anthropogenic factors impact both bark and grape bacteriomes, but to a different extent. The evidence of a “microbial terroir” seems to be even more marked in bark than in berries, possibly due to its permanence over time and to its physical proximity with soil. The importance of vine trunk bark, as potential source of inoculum for grapes and as interesting bacterial diversity habitat, is evidenced. This opens new fields of investigation, not only for researchers that aim at describing this little-known habitat within the vineyard, but also for stakeholders from the wine industry that want to understand the roles of microorganisms on the entire winemaking process, from vineyard to cellar.

Multilevel comparative bioinformatics to investigate evolutionary relationships and specificities in gene annotations: an example for tomato and grapevine

Background

“Omics” approaches may provide useful information for a deeper understanding of speciation events, diversification and function innovation. This can be achieved by investigating the molecular similarities at sequence level between species, allowing the definition of ortholog and paralog genes. However, the spreading of sequenced genome, often endowed with still preliminary annotations, requires suitable bioinformatics to be appropriately exploited in this framework.

Results

We presented here a multilevel comparative approach to investigate on genome evolutionary relationships and peculiarities of two fleshy fruit species of relevant agronomic interest, Solanum lycopersicum (tomato) and Vitis vinifera (grapevine). We defined 17,823 orthology relationships between tomato and grapevine reference gene annotations. The resulting orthologs are associated with the detected paralogs in each species, permitting the definition of gene networks, useful to investigate the different relationships. The reconciliation of the compared collections in terms of an updating of the functional descriptions was also exploited. All the results were made accessible in ComParaLogs, a dedicated bioinformatics platform available at http://biosrv.cab.unina.it/comparalogs/gene/search.

Conclusions

The aim of the work was to suggest a reliable approach to detect all similarities of gene loci between two species based on the integration of results from different levels of information, such as the gene, the transcript and the protein sequences, overcoming possible limits due to exclusive protein versus protein comparisons. This to define reliable ortholog and paralog genes, as well as species specific gene loci in the two species, overcoming limits due to the possible draft nature of preliminary gene annotations. Moreover, reconciled functional descriptions, as well as common or peculiar enzymatic classes and protein domains from tomato and grapevine, together with the definition of species-specific gene sets after the pairwise comparisons, contributed a comprehensive set of information useful to comparatively exploit the two species gene annotations and investigate on differences between species with climacteric and non-climacteric fruits. In addition, the definition of networks of ortholog genes and of associated paralogs, and the organization of web-based interfaces for the exploration of the results, defined a friendly computational bench-work in support of comparative analyses between two species.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2420-y) contains supplementary material, which is available to authorized users.

Timing and Order of the Molecular Events Marking the Onset of Berry Ripening in Grapevine

Grapevine (Vitis vinifera) is a model for the investigation of physiological and biochemical changes during the formation and ripening of nonclimacteric fleshy fruits. However, the order and complexity of the molecular events during fruit development remain poorly understood. To identify the key molecular events controlling berry formation and ripening, we created a highly detailed transcriptomic and metabolomic map of berry development, based on samples collected every week from fruit set to maturity in two grapevine genotypes for three consecutive years, resulting in 219 samples. Major transcriptomic changes were represented by coordinated waves of gene expression associated with early development, veraison (onset of ripening)/midripening, and late-ripening and were consistent across vintages. The two genotypes were clearly distinguished by metabolite profiles and transcriptional changes occurring primarily at the veraison/midripening phase. Coexpression analysis identified a core network of transcripts as well as variations in the within-module connections representing varietal differences. By focusing on transcriptome rearrangements close to veraison, we identified two rapid and successive shared transitions involving genes whose expression profiles precisely locate the timing of the molecular reprogramming of berry development. Functional analyses of two transcription factors, markers of the first transition, suggested that they participate in a hierarchical cascade of gene activation at the onset of ripening. This study defined the initial transcriptional events that mark and trigger the onset of ripening and the molecular network that characterizes the whole process of berry development, providing a framework to model fruit development and maturation in grapevine.

Biocontrol traits of Bacillus licheniformis GL174, a culturable endophyte of Vitis vinifera cv. Glera

Background

Bacillus licheniformis GL174 is a culturable endophytic strain isolated from Vitis vinifera cultivar Glera, the grapevine mainly cultivated for the Prosecco wine production. This strain was previously demonstrated to possess some specific plant growth promoting traits but its endophytic attitude and its role in biocontrol was only partially explored. In this study, the potential biocontrol action of the strain was investigated in vitro and in vivo and, by genome sequence analyses, putative functions involved in biocontrol and plant-bacteria interaction were assessed.

Results

Firstly, to confirm the endophytic behavior of the strain, its ability to colonize grapevine tissues was demonstrated and its biocontrol properties were analyzed. Antagonism test results showed that the strain could reduce and inhibit the mycelium growth of diverse plant pathogens in vitro and in vivo. The strain was demonstrated to produce different molecules of the lipopeptide class; moreover, its genome was sequenced, and analysis of the sequences revealed the presence of many protein-coding genes involved in the biocontrol process, such as transporters, plant-cell lytic enzymes, siderophores and other secondary metabolites.

Conclusions

This step-by-step analysis shows that Bacillus licheniformis GL174 may be a good biocontrol agent candidate, and describes some distinguished traits and possible key elements involved in this process. The use of this strain could potentially help grapevine plants to cope with pathogen attacks and reduce the amount of chemicals used in the vineyard.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1306-5) contains supplementary material, which is available to authorized users.

Genomic analysis of Sparus aurata reveals the evolutionary dynamics of sex-biased genes in a sequential hermaphrodite fish

Sexual dimorphism is a fascinating subject in evolutionary biology and mostly results from sex-biased expression of genes, which have been shown to evolve faster in gonochoristic species. We report here genome and sex-specific transcriptome sequencing of Sparus aurata, a sequential hermaphrodite fish. Evolutionary comparative analysis reveals that sex-biased genes in S. aurata are similar in number and function, but evolved following strikingly divergent patterns compared with gonochoristic species, showing overall slower rates because of stronger functional constraints. Fast evolution is observed only for highly ovary-biased genes due to female-specific patterns of selection that are related to the peculiar reproduction mode of S. aurata, first maturing as male, then as female. To our knowledge, these findings represent the first genome-wide analysis on sex-biased loci in a hermaphrodite vertebrate species, demonstrating how having two sexes in the same individual profoundly affects the fate of a large set of evolutionarily relevant genes.

The leaf transcriptome of fennel (Foeniculum vulgare Mill.) enables characterization of the t-anethole pathway and the discovery of microsatellites and single-nucleotide variants

Fennel is a plant species of both agronomic and pharmaceutical interest that is characterized by a shortage of genetic and molecular data. Taking advantage of NGS technology, we sequenced and annotated the first fennel leaf transcriptome using material from four different lines and two different bioinformatic approaches: de novo and genome-guided transcriptome assembly. A reference transcriptome for assembly was produced by combining these two approaches. Among the 79,263 transcripts obtained, 47,775 were annotated using BLASTX analysis performed against the NR protein database subset with 11,853 transcripts representing putative full-length CDS. Bioinformatic analyses revealed 1,011 transcripts encoding transcription factors, mainly from the BHLH, MYB-related, C2H2, MYB, and ERF families, and 6,411 EST-SSR regions. Single-nucleotide variants of SNPs and indels were identified among the 8 samples at a frequency of 0.5 and 0.04 variants per Kb, respectively. Finally, the assembled transcripts were screened to identify genes related to the biosynthesis of t-anethole, a compound well-known for its nutraceutical and medical properties. For each of the 11 genes encoding structural enzymes in the t-anethole biosynthetic pathway, we identified at least one transcript showing a significant match. Overall, our work represents a treasure trove of information exploitable both for marker-assisted breeding and for in-depth studies on thousands of genes, including those involved in t-anethole biosynthesis.

High-Throughput Sequencing of microRNAs in Glucocorticoid Sensitive Paediatric Inflammatory Bowel Disease Patients

The aim of this research was the identification of novel pharmacogenomic biomarkers for better understanding the complex gene regulation mechanisms underpinning glucocorticoid (GC) action in paediatric inflammatory bowel disease (IBD). This goal was achieved by evaluating high-throughput microRNA (miRNA) profiles during GC treatment, integrated with the assessment of expression changes in GC receptor (GR) heterocomplex genes. Furthermore, we tested the hypothesis that differentially expressed miRNAs could be directly regulated by GCs through investigating the presence of GC responsive elements (GREs) in their gene promoters. Ten IBD paediatric patients responding to GCs were enrolled. Peripheral blood was obtained at diagnosis (T0) and after four weeks of steroid treatment (T4). MicroRNA profiles were analyzed using next generation sequencing, and selected significantly differentially expressed miRNAs were validated by quantitative reverse transcription-polymerase chain reaction. In detail, 18 miRNAs were differentially expressed from T0 to T4, 16 of which were upregulated and 2 of which were downregulated. Out of these, three miRNAs (miR-144, miR-142, and miR-96) could putatively recognize the 3’UTR of the GR gene and three miRNAs (miR-363, miR-96, miR-142) contained GREs sequences, thereby potentially enabling direct regulation by the GR. In conclusion, we identified miRNAs differently expressed during GC treatment and miRNAs which could be directly regulated by GCs in blood cells of young IBD patients. These results could represent a first step towards their translation as pharmacogenomic biomarkers.

Genome-Wide Transcriptional Changes and Lipid Profile Modifications Induced by Medicago truncatula N5 Overexpression at an Early Stage of the Symbiotic Interaction with Sinorhizobium meliloti

Plant lipid-transfer proteins (LTPs) are small basic secreted proteins, which are characterized by lipid-binding capacity and are putatively involved in lipid trafficking. LTPs play a role in several biological processes, including the root nodule symbiosis. In this regard, the Medicago truncatula nodulin 5 (MtN5) LTP has been proved to positively regulate the nodulation capacity, controlling rhizobial infection and nodule primordia invasion. To better define the lipid transfer protein MtN5 function during the symbiosis, we produced MtN5-downregulated and -overexpressing plants, and we analysed the transcriptomic changes occurring in the roots at an early stage of Sinorhizobium meliloti infection. We also carried out the lipid profile analysis of wild type (WT) and MtN5-overexpressing roots after rhizobia infection. The downregulation of MtN5 increased the root hair curling, an early event of rhizobia infection, and concomitantly induced changes in the expression of defence-related genes. On the other hand, MtN5 overexpression favoured the invasion of the nodules by rhizobia and determined in the roots the modulation of genes that are involved in lipid transport and metabolism as well as an increased content of lipids, especially galactolipids that characterize the symbiosome membranes. Our findings suggest the potential participation of LTPs in the synthesis and rearrangement of membranes occurring during the formation of the infection threads and the symbiosome membrane.

Characterization and Evolution of the Cell Cycle-Associated Mob Domain-Containing Proteins in Eukaryotes

The MOB family includes a group of cell cycle-associated proteins highly conserved throughout eukaryotes, whose founding members are implicated in mitotic exit and co-ordination of cell cycle progression with cell polarity and morphogenesis. Here we report the characterization and evolution of the MOB domain-containing proteins as inferred from the 43 eukaryotic genomes so far sequenced. We show that genes for Mob-like proteins are present in at least 41 of these genomes, confirming the universal distribution of this protein family and suggesting its prominent biological function. The phylogenetic analysis reveals five distinct MOB domain classes, showing a progressive expansion of this family from unicellular to multicellular organisms, reaching the highest number in mammals. Plant Mob genes appear to have evolved from a single ancestor, most likely after the loss of one or more genes during the early stage of Viridiplantae evolutionary history. Three of the Mob classes are widespread among most of the analyzed organisms. The possible biological and molecular function of Mob proteins and their role in conserved signaling pathways related to cell proliferation, cell death and cell polarity are also presented and critically discussed.

Downregulation of lizard immuno-genes in the regenerating tail and myogenes in the scarring limb suggests that tail regeneration occurs in an immuno-privileged organ

Amputated tails of lizards regenerate while limbs form scars which histological structure is very different from the original organs. Lizards provide useful information for regenerative medicine and some hypotheses on the loss of regeneration in terrestrial vertebrates. Analysis of tail and limb transcriptomes shows strong downregulation in the tail blastema for immunoglobulins and surface B and T receptors, cell function, and metabolism. In contrast, in the limb blastema genes for myogenesis, muscle and cell function, and extracellular matrix deposition but not immunity are variably downregulated. The upregulated genes show that the regenerating tail is an embryonic organ driven by the Wnt pathway and non-coding RNAs. The strong inflammation following amputation, the non-activation of the Wnt pathway, and the upregulation of inflammatory genes with no downregulation of immune genes indicate that the amputated limb does not activate an embryonic program. Intense inflammation in limbs influences in particular the activity of genes coding for muscle proteins, cell functions, and stimulates the deposition of dense extracellular matrix proteins resulting in scarring limb outgrowths devoid of muscles. The present study complements that on upregulated genes, and indicates that the regenerating tail requires immune suppression to maintain this embryonic organ connected to the rest of the tail without be rejected or turned into a scar. It is hypothesized that the evolution of the adaptive immune system determined scarring instead of organ regeneration in terrestrial vertebrates and that lizards evolved the process of tail regeneration through a mechanism of immuno-evasion.

QueryOR: a comprehensive web platform for genetic variant analysis and prioritization

Background

Whole genome and exome sequencing are contributing to the extraordinary progress in the study of human genetic variants. In this fast developing field, appropriate and easily accessible tools are required to facilitate data analysis.

Results

Here we describe QueryOR, a web platform suitable for searching among known candidate genes as well as for finding novel gene-disease associations. QueryOR combines several innovative features that make it comprehensive, flexible and easy to use. Instead of being designed on specific datasets, it works on a general XML schema specifying formats and criteria of each data source. Thanks to this flexibility, new criteria can be easily added for future expansion. Currently, up to 70 user-selectable criteria are available, including a wide range of gene and variant features. Moreover, rather than progressively discarding variants taking one criterion at a time, the prioritization is achieved by a global positive selection process that considers all transcript isoforms, thus producing reliable results. QueryOR is easy to use and its intuitive interface allows to handle different kinds of inheritance as well as features related to sharing variants in different patients. QueryOR is suitable for investigating single patients, families or cohorts.

Conclusions

QueryOR is a comprehensive and flexible web platform eligible for an easy user-driven variant prioritization. It is freely available for academic institutions at http://queryor.cribi.unipd.it/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1654-4) contains supplementary material, which is available to authorized users.

The Accumulation of miRNAs Differentially Modulated by Drought Stress Is Affected by Grafting in Grapevine

Grapevine (Vitis vinifera) is routinely grafted, and rootstocks inducing drought tolerance represent a source for adapting vineyards to climate change in temperate areas. Our goal was to investigate drought stress effects on microRNA (miRNA) abundance in a drought-resistant grapevine rootstock, M4 (Vitis vinifera × Vitis berlandieri), compared with a commercial cultivar, Cabernet Sauvignon, using their autografts and reciprocal grafts. RNA extracted from roots and leaves of droughted and irrigated plants of different graft combinations was used to prepare cDNA libraries for small RNA sequencing and to analyze miRNAs by quantitative real-time polymerase chain reaction (RT-qPCR). Measurements of leaf water potential, leaf gas exchange, and root hydraulic conductance attested that, under irrigation, M4 reduced water loss in comparison with cultivar Cabernet Sauvignon mostly through nonhydraulic, root-specific mechanisms. Under drought, stomatal conductance decreased at similar levels in the two genotypes. Small RNA sequencing allowed the identification of 70 conserved miRNAs and the prediction of 28 novel miRNAs. Different accumulation trends of miRNAs, observed upon drought and in different genotypes and organs, were confirmed by RT-qPCR Corresponding target transcripts, predicted in silico and validated by RT-qPCR, often showed opposite expression profiles than the related miRNAs. Drought effects on miRNA abundance differed between the two genotypes. Furthermore, the concentration of drought-responsive miRNAs in each genotype was affected by reciprocal grafting, suggesting either the movement of signals inducing miRNA expression in the graft partner or, possibly, miRNA transport between scion and rootstock. These results open new perspectives in the selection of rootstocks for improving grapevine adaptation to drought.

Transcriptional Characterization of a Widely-Used Grapevine Rootstock Genotype under Different Iron-Limited Conditions

Iron chlorosis is a serious deficiency that affects orchards and vineyards reducing quality and yield production. Chlorotic plants show abnormal photosynthesis and yellowing shoots. In grapevine iron uptake and homeostasis are most likely controlled by a mechanism known as "Strategy I," characteristic of non-graminaceous plants and based on a system of soil acidification, iron reduction and transporter-mediated uptake. Nowadays, grafting of varieties of economic interest on tolerant rootstocks is widely used practice against many biotic and abiotic stresses. Nevertheless, many interspecific rootstocks, and in particular those obtained by crossing exclusively non-vinifera genotypes, can show limited nutrient uptake and transport, in particular for what concerns iron. In the present study, 101.14, a commonly used rootstock characterized by susceptibility to iron chlorosis was subjected to both Fe-absence and Fe-limiting conditions. Grapevine plantlets were grown in control, Fe-deprived, and bicarbonate-supplemented hydroponic solutions. Whole transcriptome analyses, via mRNA-Seq, were performed on root apices of stressed and unstressed plants. Analysis of differentially expressed genes (DEGs) confirmed that Strategy I is the mechanism responsible for iron uptake in grapevine, since many orthologs genes to the Arabidopsis "ferrome" were differentially regulated in stressed plant. Molecular differences in the plant responses to Fe absence and presence of bicarbonate were also identified indicating the two treatments are able to induce response-mechanisms only partially overlapping. Finally, we measured the expression of a subset of genes differentially expressed in 101.14 (such as IRT1, FERRITIN1, bHLH38/39) or known to be fundamental in the "strategy I" mechanism (AHA2 and FRO2) also in a tolerant rootstock (M1) finding important differences which could be responsible for the different degrees of tolerance observed.

Light Remodels Lipid Biosynthesis in Nannochloropsis gaditana by Modulating Carbon Partitioning between Organelles

The seawater microalga Nannochloropsis gaditana is capable of accumulating a large fraction of reduced carbon as lipids. To clarify the molecular bases of this metabolic feature, we investigated light-driven lipid biosynthesis in Nannochloropsis gaditana cultures combining the analysis of photosynthetic functionality with transcriptomic, lipidomic and metabolomic approaches. Light-dependent alterations are observed in amino acid, isoprenoid, nucleic acid, and vitamin biosynthesis, suggesting a deep remodeling in the microalgal metabolism triggered by photoadaptation. In particular, high light intensity is shown to affect lipid biosynthesis, inducing the accumulation of diacylglyceryl-N,N,N-trimethylhomo-Ser and triacylglycerols, together with the up-regulation of genes involved in their biosynthesis. Chloroplast polar lipids are instead decreased. This situation correlates with the induction of genes coding for a putative cytosolic fatty acid synthase of type 1 (FAS1) and polyketide synthase (PKS) and the down-regulation of the chloroplast fatty acid synthase of type 2 (FAS2). Lipid accumulation is accompanied by the regulation of triose phosphate/inorganic phosphate transport across the chloroplast membranes, tuning the carbon metabolic allocation between cell compartments, favoring the cytoplasm, mitochondrion, and endoplasmic reticulum at the expense of the chloroplast. These results highlight the high flexibility of lipid biosynthesis in N. gaditana and lay the foundations for a hypothetical mechanism of regulation of primary carbon partitioning by controlling metabolite allocation at the subcellular level.

Transcriptome dynamics in the asexual cycle of the chordate Botryllus schlosseri

Background

We performed an analysis of the transcriptome during the blastogenesis of the chordate Botryllus schlosseri, focusing in particular on genes involved in cell death by apoptosis. The tunicate B. schlosseri is an ascidian forming colonies characterized by the coexistence of three blastogenetic generations: filter-feeding adults, buds on adults, and budlets on buds. Cyclically, adult tissues undergo apoptosis and are progressively resorbed and replaced by their buds originated by asexual reproduction. This is a feature of colonial tunicates, the only known chordates that can reproduce asexually.

Results

Thanks to a newly developed web-based platform (http://botryllus.cribi.unipd.it), we compared the transcriptomes of the mid-cycle, the pre-take-over, and the take-over phases of the colonial blastogenetic cycle. The platform is equipped with programs for comparative analysis and allows to select the statistical stringency. We enriched the genome annotation with 11,337 new genes; 581 transcripts were resolved as complete open reading frames, translated in silico into amino acid sequences and then aligned onto the non-redundant sequence database. Significant differentially expressed genes were classified within the gene ontology categories. Among them, we recognized genes involved in apoptosis activation, de-activation, and regulation.

Conclusions

With the current work, we contributed to the improvement of the first released B. schlosseri genome assembly and offer an overview of the transcriptome changes during the blastogenetic cycle, showing up- and down-regulated genes. These results are important for the comprehension of the events underlying colony growth and regression, cell proliferation, colony homeostasis, and competition among different generations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2598-1) contains supplementary material, which is available to authorized users.