A chromosome-scale assembly reveals chromosomal aberrations and exchanges generating genetic diversity in Coffea arabica germplasm

In order to better understand the mechanisms generating genetic diversity in the recent allotetraploid species Coffea arabica, here we present a chromosome-level assembly obtained with long read technology. Two genomic compartments with different structural and functional properties are identified in the two homoeologous genomes. The resequencing data from a large set of accessions reveals low intraspecific diversity in the center of origin of the species. Across a limited number of genomic regions, diversity increases in some cultivated genotypes to levels similar to those observed within one of the progenitor species, Coffea canephora, presumably as a consequence of introgressions deriving from the so-called Timor hybrid. It also reveals that, in addition to few, early-occurring exchanges between homoeologous chromosomes, there are numerous recent chromosomal aberrations including aneuploidies, deletions, duplications and exchanges. These events are still polymorphic in the germplasm and could represent a fundamental source of genetic variation in such a lowly variable species.

Lignocellulolytic Potential of Microbial Consortia Isolated from a Local Biogas Plant: The Case of Thermostable Xylanases Secreted by Mesophilic Bacteria

Lignocellulose biomasses (LCB), including spent mushroom substrate (SMS), pose environmental challenges if not properly managed. At the same time, these renewable resources hold immense potential for biofuel and chemicals production. With the mushroom market growth expected to amplify SMS quantities, repurposing or disposal strategies are critical. This study explores the use of SMS for cultivating microbial communities to produce carbohydrate-active enzymes (CAZymes). Addressing a research gap in using anaerobic digesters for enriching microbiomes feeding on SMS, this study investigates microbial diversity and secreted CAZymes under varied temperatures (37 °C, 50 °C, and 70 °C) and substrates (SMS as well as pure carboxymethylcellulose, and xylan). Enriched microbiomes demonstrated temperature-dependent preferences for cellulose, hemicellulose, and lignin degradation, supported by thermal and elemental analyses. Enzyme assays confirmed lignocellulolytic enzyme secretion correlating with substrate degradation trends. Notably, thermogravimetric analysis (TGA), coupled with differential scanning calorimetry (TGA-DSC), emerged as a rapid approach for saccharification potential determination of LCB. Microbiomes isolated at mesophilic temperature secreted thermophilic hemicellulases exhibiting robust stability and superior enzymatic activity compared to commercial enzymes, aligning with biorefinery conditions. PCR-DGGE and metagenomic analyses showcased dynamic shifts in microbiome composition and functional potential based on environmental conditions, impacting CAZyme abundance and diversity. The meta-functional analysis emphasised the role of CAZymes in biomass transformation, indicating microbial strategies for lignocellulose degradation. Temperature and substrate specificity influenced the degradative potential, highlighting the complexity of environmental-microbial interactions. This study demonstrates a temperature-driven microbial selection for lignocellulose degradation, unveiling thermophilic xylanases with industrial promise. Insights gained contribute to optimizing enzyme production and formulating efficient biomass conversion strategies. Understanding microbial consortia responses to temperature and substrate variations elucidates bioconversion dynamics, emphasizing tailored strategies for harnessing their biotechnological potential.

CK2β Regulates Hematopoietic Stem Cell Biology and Erythropoiesis

The Ser-Thr kinase CK2 plays important roles in sustaining cell survival and resistance to stress and these functions are exploited by different types of blood tumors. Yet, the physiological involvement of CK2 in normal blood cell development is poorly known. Here, we discovered that the β regulatory subunit of CK2 is critical for normal hematopoiesis in the mouse. Fetal livers of conditional CK2β knockout embryos showed increased numbers of hematopoietic stem cells associated to a higher proliferation rate compared to control animals. Both hematopoietic stem and progenitor cells (HSPCs) displayed alterations in the expression of transcription factors involved in cell quiescence, self-renewal, and lineage commitment. HSPCs lacking CK2β were functionally impaired in supporting both in vitro and in vivo hematopoiesis as demonstrated by transplantation assays. Furthermore, KO mice developed anemia due to a reduced number of mature erythroid cells. This compartment was characterized by dysplasia, proliferative defects at early precursor stage, and apoptosis at late-stage erythroblasts. Erythroid cells exhibited a marked compromise of signaling cascades downstream of the cKit and erythropoietin receptor, with a defective activation of ERK/JNK, JAK/STAT5, and PI3K/AKT pathways and perturbations of several transcriptional programs as demonstrated by RNA-Seq analysis. Moreover, we unraveled an unforeseen molecular mechanism whereby CK2 sustains GATA1 stability and transcriptional proficiency. Thus, our work demonstrates new and crucial functions of CK2 in HSPC biology and in erythropoiesis.

Exploring the Effects of Probiotic Treatment on Urinary and Serum Metabolic Profiles in Healthy Individuals

Probiotics are live microorganisms that confer health benefits when administered in adequate amounts. They are used to promote gut health and alleviate various disorders. Recently, there has been an increasing interest in the potential effects of probiotics on human physiology. In the presented study, the effects of probiotic treatment on the metabolic profiles of human urine and serum using a nuclear magnetic resonance (NMR)-based metabonomic approach were investigated. Twenty-one healthy volunteers were enrolled in the study, and they received two different dosages of probiotics for 8 weeks. During the study, urine and serum samples were collected from volunteers before and during probiotic supplementation. The results showed that probiotics had a significant impact on the urinary and serum metabolic profiles without altering their phenotypes. This study demonstrated the effects of probiotics in terms of variations of metabolite levels resulting also from the different probiotic posology. Overall, the results suggest that probiotic administration may affect both urine and serum metabolomes, although more research is needed to understand the mechanisms and clinical implications of these effects. NMR-based metabonomic analysis of biofluids is a powerful tool for monitoring host-gut microflora dynamic interaction as well as for assessing the individual response to probiotic treatment.

MicroRNAs associated with AGL6 and IAA9 function in tomato fruit set

OBJECTIVE

Fruit set is triggered after ovule fertilization, as a consequence of the downregulation of ovary growth repressors, such as the tomato transcription factors Auxin/indole-3-acetic acid 9 (IAA9) and Agamous-like 6 (AGL6). In a recent work, we developed a method to silence IAA9 and AGL6 in tomato ovaries using exogenous dsRNAs. We also produced small RNA libraries from IAA9- and AGL6-silenced ovaries to confirm the presence of siRNAs, derived from exogenous dsRNA, targeting IAA9 and AGL6. The objective of this work is to exploit these sRNA libraries to identify miRNAs differentially expressed in IAA9- and AGL6-silenced ovaries as compared with unpollinated control ovaries.

RESULTS

We identified by RNA sequencing 125 and 104 known and 509 and 516 novel miRNAs from reads mapped to mature or hairpin sequences, respectively. Of the known miRNAs, 7 and 45 were differentially expressed in IAA9- and AGL6-silenced ovaries compared to control ones, respectively. Six miRNAs were common to both datasets, suggesting their importance in the fruit set process. The expression pattern of two of these (miR393 and miR482e-5p) was verified by stem-loop qRT-PCR. The identified miRNAs represent a pool of regulatory sRNAs potentially involved in tomato fruit initiation.

An improved reference of the grapevine genome reasserts the origin of the PN40024 highly-homozygous genotype

The genome sequence of the diploid and highly homozygous V. vinifera genotype PN40024 serves as the reference for many grapevine studies. Despite several improvements to the PN40024 genome assembly, its current version PN12X.v2 is quite fragmented and only represents the haploid state of the genome with mixed haplotypes. In fact, being nearly homozygous, this genome contains several heterozygous regions that are yet to be resolved. Taking the opportunity of improvements that long-read sequencing technologies offer to fully discriminate haplotype sequences, an improved version of the reference, called PN40024.v4, was generated. Through incorporating long genomic sequencing reads to the assembly, the continuity of the 12X.v2 scaffolds was highly increased with a total number decreasing from 2,059 to 640 and a reduction in N bases of 88%. Additionally, the full alternative haplotype sequence was built for the first time, the chromosome anchoring was improved and the number of unplaced scaffolds was reduced by half. To obtain a high-quality gene annotation that outperforms previous versions, a liftover approach was complemented with an optimized annotation workflow for Vitis. Integration of the gene reference catalogue and its manual curation have also assisted in improving the annotation, while defining the most reliable estimation of 35,230 genes to date. Finally, we demonstrated that PN40024 resulted from nine selfings of cv. 'Helfensteiner' (cross of cv. 'Pinot noir' and 'Schiava grossa') instead of a single 'Pinot noir'. These advances will help maintain the PN40024 genome as a gold-standard reference, also contributing towards the eventual elaboration of the grapevine pangenome.

CK2β-regulated signaling controls B cell differentiation and function

Serine-Threonine kinase CK2 supports malignant B-lymphocyte growth but its role in B-cell development and activation is largely unknown. Here, we describe the first B-cell specific knockout (KO) mouse model of the β regulatory subunit of CK2. CK2β^KO mice present an increase in marginal zone (MZ) and a reduction in follicular B cells, suggesting a role for CK2 in the regulation of the B cell receptor (BCR) and NOTCH2 signaling pathways. Biochemical analyses demonstrate an increased activation of the NOTCH2 pathway in CK2β^KO animals, which sustains MZ B-cell development. Transcriptomic analyses indicate alterations in biological processes involved in immune response and B-cell activation. Upon sheep red blood cells (SRBC) immunization CK2β^KO mice exhibit enlarged germinal centers (GCs) but display a limited capacity to generate class-switched GC B cells and immunoglobulins. In vitro assays highlight that B cells lacking CK2β have an impaired signaling downstream of BCR, Toll-like receptor, CD40, and IL-4R all crucial for B-cell activation and antigen presenting efficiency. Somatic hypermutations analysis upon 4-Hydroxy-3-nitrophenylacetyl hapten conjugated to Chicken Gamma Globulin (NP-CGG) evidences a reduced NP-specific W33L mutation frequency in CK2β^KO mice suggesting the importance of the β subunit in sustaining antibody affinity maturation. Lastly, since diffuse large B cell lymphoma (DLBCL) cells derive from GC or post-GC B cells and rely on CK2 for their survival, we sought to investigate the consequences of CK2 inhibition on B cell signaling in DLBCL cells. In line with the observations in our murine model, CK2 inactivation leads to signaling defects in pathways that are essential for malignant B-lymphocyte activation.

benchdamic: benchmarking of differential abundance methods for microbiome data

SUMMARY

Recently, an increasing number of methodological approaches have been proposed to tackle the complexity of metagenomics and microbiome data. In this scenario, reproducibility and replicability have become two critical issues, and the development of computational frameworks for the comparative evaluations of such methods is of utmost importance. Here, we present benchdamic, a Bioconductor package to benchmark methods for the identification of differentially abundant taxa.

AVAILABILITY AND IMPLEMENTATION

benchdamic is available as an open-source R package through the Bioconductor project at https://bioconductor.org/packages/benchdamic/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A red wine intervention does not modify plasma trimethylamine N-oxide but is associated with broad shifts in the plasma metabolome and gut microbiota composition

BACKGROUND

Gut microbiota profiles are closely related to cardiovascular diseases through mechanisms that include the reported deleterious effects of metabolites, such as trimethylamine N-oxide (TMAO), which have been studied as diagnostic and therapeutic targets. Moderate red wine (RW) consumption is reportedly cardioprotective, possibly by affecting the gut microbiota.

OBJECTIVES

To investigate the effects of RW consumption on the gut microbiota, plasma TMAO, and the plasma metabolome in men with documented coronary artery disease (CAD) using a multiomics assessment in a crossover trial.

METHODS

We conducted a randomized, crossover, controlled trial involving 42 men (average age, 60 y) with documented CAD comparing 3-wk RW consumption (250 mL/d, 5 d/wk) with an equal period of alcohol abstention, both preceded by a 2-wk washout period. The gut microbiota was analyzed via 16S rRNA high-throughput sequencing. Plasma TMAO was evaluated by LC-MS/MS. The plasma metabolome of 20 randomly selected participants was evaluated by ultra-high-performance LC-MS/MS. The effect of RW consumption was assessed by individual comparisons using paired tests during the abstention and RW periods.

RESULTS

Plasma TMAO did not differ between RW intervention and alcohol abstention, and TMAO concentrations showed low intraindividual concordance over time, with an intraclass correlation coefficient of 0.049 during the control period. After RW consumption, there was significant remodeling of the gut microbiota, with a difference in β diversity and predominance of Parasutterella, Ruminococcaceae, several Bacteroides species, and Prevotella. Plasma metabolomic analysis revealed significant changes in metabolites after RW consumption, consistent with improved redox homeostasis.

CONCLUSIONS

Modulation of the gut microbiota may contribute to the putative cardiovascular benefits of moderate RW consumption. The low intraindividual concordance of TMAO presents challenges regarding its role as a cardiovascular risk biomarker at the individual level. This study was registered at clinical trials.gov as NCT03232099.

MetaPhage: an Automated Pipeline for Analyzing, Annotating, and Classifying Bacteriophages in Metagenomics Sequencing Data

Phages are the most abundant biological entities on the planet, and they play an important role in controlling density, diversity, and network interactions among bacterial communities through predation and gene transfer. To date, a variety of bacteriophage identification tools have been developed that differ in the phage mining strategies used, input files requested, and results produced. However, new users attempting bacteriophage analysis can struggle to select the best methods and interpret the variety of results produced. Here, we present MetaPhage, a comprehensive reads-to-report pipeline that streamlines the use of multiple phage miners and generates an exhaustive report. The report both summarizes and visualizes the key findings and enables further exploration of key results via interactive filterable tables. The pipeline is implemented in Nextflow, a widely adopted workflow manager that enables an optimized parallelization of tasks in different locations, from local server to the cloud; this ensures reproducible results from containerized packages. MetaPhage is designed to enable scalability and reproducibility; also, it can be easily expanded to include new miners and methods as they are developed in this continuously growing field. MetaPhage is freely available under a GPL-3.0 license at https://github.com/MattiaPandolfoVR/MetaPhage. IMPORTANCE Bacteriophages (viruses that infect bacteria) are the most abundant biological entities on earth and are increasingly studied as members of the resident microbiota community in many environments, from oceans to soils and the human gut. Their identification is of great importance to better understand complex bacterial dynamics and microbial ecosystem function. A variety of metagenome bacteriophage identification tools have been developed that differ in the phage mining strategies used, input files requested, and results produced. To facilitate the management and the execution of such a complex workflow, we developed MetaPhage (MP), a comprehensive reads-to-report pipeline that streamlines the use of multiple phage miners and generates an exhaustive report. The pipeline is implemented in Nextflow, a widely adopted workflow manager that enables an optimized parallelization of tasks. MetaPhage is designed to enable scalability and reproducibility and offers an installation-free, dependency-free, and conflict-free workflow execution.

Nanovector-mediated exogenous delivery of dsRNA induces silencing of target genes in very young tomato flower buds

RNA interference (RNAi) is a post-translational regulatory mechanism that controls gene expression in plants. This process can be artificially induced by double-stranded RNA (dsRNA) molecules with sequence homology to target mRNAs. Exogenously applied dsRNA on leaves has been shown to silence virulence genes of fungi and viruses, conferring protection to plants. Coupling dsRNA to nanoparticles has been demonstrated to prolong the silencing effect. The ability of exogenous dsRNA to silence endogenous genes in plants is currently under debate, mainly due to the difficulty in delivering dsRNA into plant tissues and organs. Our study aims to develop a method based on the exogenous application of dsRNA on tomato flowers for silencing endogenous genes controlling ovary growth. Two methods of dsRNA delivery into tomato flower buds (i.e., pedicel soaking and injection) were compared to test their efficacy in silencing the tomato Aux/IAA9 (SlIAA9) gene, which encodes for a known repressor of ovary growth. We examined the silencing effect of dsRNA alone and coupled to layered double hydroxide (LDHs) nanoparticles. We found that injection into the pedicel led to the silencing of SlIAA9 and the efficacy of the method was confirmed by choosing a different ovary growth repressor gene (SlAGAMOUS-like 6; SlAGL6). The coupling of dsRNA to LDHs increased the silencing effect in the case of SlIAA9. Silencing of the two repressors caused an increase in ovary size only when flower buds were treated with dsRNA coupled to LDHs. RNA-Seq of small RNAs showed that induction of RNAi was caused by the processing of injected dsRNA. In this work, we demonstrate for the first time that exogenous dsRNA coupled to LDHs can induce post-transcriptional gene silencing in the young tomato ovary by injection into the flower pedicel. This method represents a silencing tool for the study of the molecular changes occurring during the early stages of ovary/fruit growth as a consequence of downregulation of target genes, without the need to produce transgenic plants stably expressing RNAi constructs.

Randomized Double-Blind Crossover Study for Evaluating a Probiotic Mixture on Gastrointestinal and Behavioral Symptoms of Autistic Children

Autism spectrum disorders (ASDs) represent a diagnostic challenge with a still partially uncertain etiology, in which genetic and environmental factors have now been assessed. Among the hypotheses underlying the involvement of biological and environmental factors, the gut–brain axis is of particular interest in autism spectrum disorders. Several studies have highlighted the related incidence of particular gastrointestinal symptoms (GISs) in children suffering from ASDs. Probiotics have shown success in treating several gastrointestinal dysbiotic disorders; therefore, it is plausible to investigate whether they can alleviate behavioral symptoms as well. On these bases, a randomized double-blind crossover study with a placebo was conducted, evaluating the effects of a mixture of probiotics in a group of 61 subjects aged between 24 months and 16 years old with a diagnosis of ASD. Behavioral evaluation was performed through the administration of a questionnaire including a Parenting Stress Index (PSI) test and the Vineland Adaptive Behavior Scale (VABS). The Psycho-Educational Profile and the Autism Spectrum Rating Scale (ASRS) were also evaluated. Microbial composition analyses of fecal samples of the two groups was also performed. The study showed significant improvements in GISs, communication skills, maladaptive behaviors, and perceived parental stress level after the administration of probiotics. Microbiome alpha diversity was comparable between treatment arms and no significant differences were found, although beta diversity results were significantly different in the treatment group between T0 and T1 time points. Streptococcus thermophilus, Bifidobacterium longum, Limosilactobacillus fermentum, and Ligilactobacillus salivarius species were identified as some of the most discriminant taxa positively associated with T1 samples. This preliminary study corroborates the relationship between intestinal microbiota and ASD recently described in the literature.

Salivary microbiota composition may discriminate between patients with eosinophilic oesophagitis (EoE) and non-EoE subjects

BACKGROUND

Data on the role of the microbiome in adult patients with eosinophilic oesophagitis (EoE) are limited.

AIMS

To prospectively collect and characterise the salivary, oesophageal and gastric microbiome in patients with EoE, further correlating the findings with disease activity.

METHODS

Adult patients with symptoms of oesophageal dysfunction undergoing upper endoscopy were consecutively enrolled. Patients were classified as EoE patients, in case of more than 15 eosinophils per high-power field, or non-EoE controls, in case of lack of eosinophilic infiltration. Before and during endoscopy, saliva, oesophageal and gastric fundus biopsies were collected. Microbiota assessment was performed by 16 s rRNA analysis. A Sparse Partial Least Squares Discriminant Analysis (sPLS-DA) was implemented to identify biomarkers.

RESULTS

Saliva samples were collected from 29 EoE patients and 20 non-EoE controls;, biopsies from 25 EoE and 5 non-EoE controls. In saliva samples, 23 Amplicon Sequence Variants (ASVs) were positively associated with EoE and 27 ASVs with controls, making it possible to discriminate between EoE and non-EoE patients with a classification error (CE) of 24%. In a validation cohort, the accuracy, sensitivity, specificity, positive predictive value and negative predictive value of this model were 78.6%, 80%, 75%, 80% and 60%, respectively. Moreover, the analysis of oesophageal microbiota samples observed a clear microbial pattern able to discriminate between active and inactive EoE (CE = 8%).

CONCLUSION

Our preliminary data suggest that salivary metabarcoding analysis in combination with machine learning approaches could become a valid, cheap, non-invasive test to segregate between EoE and non-EoE patients.

Comparative analysis identifies micro-RNA associated with nutrient homeostasis, development and stress response in Arabidopsis thaliana upon high Zn and metal hyperaccumulator Arabidopsis halleri

miRNAs have been found to be key players in mineral homeostasis, both in the control of nutrient balance and in the response to toxic trace elements. However, the effect of Zn excess on miRNAs has not been elucidated; moreover, no data are present regarding miRNAs in hyperaccumulator species, where metal homeostasis is tightly regulated. Therefore, expression levels of mature miRNAs were measured by RNA-Seq in Zn-sensitive Arabidopsis thaliana grown in control conditions and upon high Zn, in soil and in Zn-hyperaccumulator Arabidopsis halleri grown in control conditions. Differential expression of notable miRNAs and their targets was confirmed by real-time RT-PCR. The comparison in A. thaliana revealed a small subset modulated upon Zn treatment that is associated with stress response and nutrient homeostasis. On the other hand, a more consistent group of miRNAs was differentially expressed in A. halleri compared with A. thaliana, reflecting inherent differences in nutritional requirements and response to stresses and plant growth and development. Overall, these results confirm the involvement of miRNAs in Zn homeostasis and support the hypothesis of distinct regulatory pathways in hyperaccumulator species.

A novel Penicillium sumatraense isolate reveals an arsenal of degrading enzymes exploitable in algal bio-refinery processes

BACKGROUND

Microalgae are coming to the spotlight due to their potential applications in a wide number of fields ranging from the biofuel to the pharmaceutical sector. However, several factors such as low productivity, expensive harvesting procedures and difficult metabolite extractability limit their full utilization at industrial scale. Similarly to the successful employment of enzymatic arsenals from lignocellulolytic fungi to convert lignocellulose into fermentable sugars for bioethanol production, specific algalytic formulations could be used to improve the extractability of lipids from microalgae to produce biodiesel. Currently, the research areas related to algivorous organisms, algal saprophytes and the enzymes responsible for the hydrolysis of algal cell wall are still little explored.

RESULTS

Here, an algal trap method for capturing actively growing microorganisms was successfully used to isolate a filamentous fungus, that was identified by whole-genome sequencing, assembly and annotation as a novel Penicillium sumatraense isolate. The fungus, classified as P. sumatraense AQ67100, was able to assimilate heat-killed Chlorella vulgaris cells by an enzymatic arsenal composed of proteases such as dipeptidyl- and amino-peptidases, β-1,3-glucanases and glycosidases including α- and β-glucosidases, β-glucuronidase, α-mannosidases and β-galactosidases. The treatment of C. vulgaris with the filtrate from P. sumatraense AQ67100 increased the release of chlorophylls and lipids from the algal cells by 42.6 and 48.9%, respectively.

CONCLUSIONS

The improved lipid extractability from C. vulgaris biomass treated with the fungal filtrate highlighted the potential of algal saprophytes in the bioprocessing of microalgae, posing the basis for the sustainable transformation of algal metabolites into biofuel-related compounds.

VviNAC33 promotes organ de-greening and represses vegetative growth during the vegetative-to-mature phase transition in grapevine

Plants undergo several developmental transitions during their life cycle. In grapevine, a perennial woody fruit crop, the transition from vegetative/green-to-mature/woody growth involves transcriptomic reprogramming orchestrated by a small group of genes encoding regulators, but the underlying molecular mechanisms are not fully understood. We investigated the function of the transcriptional regulator VviNAC33 by generating and characterizing transgenic overexpressing grapevine lines and a chimeric repressor, and by exploring its putative targets through a DNA affinity purification sequencing (DAP-seq) approach combined with transcriptomic data. We demonstrated that VviNAC33 induces leaf de-greening, inhibits organ growth and directly activates the expression of STAY-GREEN PROTEIN 1 (SGR1), which is involved in Chl and photosystem degradation, and AUTOPHAGY 8f (ATG8f), which is involved in the maturation of autophagosomes. Furthermore, we show that VviNAC33 directly inhibits AUXIN EFFLUX FACILITATOR PIN1, RopGEF1 and ATP SYNTHASE GAMMA CHAIN 1T (ATPC1), which are involved in photosystem II integrity and activity. Our results show that VviNAC33 plays a major role in terminating photosynthetic activity and organ growth as part of a regulatory network governing the vegetative-to-mature phase transition.

Supplementation with Bifidobacterium breve BR03 and B632 strains improved insulin sensitivity in children and adolescents with obesity in a cross-over, randomized double-blind placebo-controlled trial

BACKGROUND & AIMS

Variations in gut microbiota might impact metabolism leading to body weight excess. We assessed the impact of a probiotic supplementation in pediatric obesity on weight, metabolic alterations, selected gut microbial groups, and functionality.

METHODS

Cross-over, double-blind, randomized control trial (BIFI-OBESE trial; NCT03261466). 101 youths (6-18 years, Tanner stage ≥2) with obesity and insulin-resistance on diet were randomized to 2 × 10⁹ CFU/AFU/day of Bifidobacterium breve BR03 (DSM 16604) and B. breve B632 (DSM 24706) (51) or placebo (50) for 8 weeks with a 4-weeks wash-out period.

RESULTS

All subjects (M/F 54/47) completed the first 8 weeks, and 82 (M/F 43/39) the last part without adverse events. Mixed-effects models revealed a carry-over effect on many variables in the entire study, narrowing the analysis to the first 8 weeks before the wash-out periods. All subjects improved metabolic parameters, and decreased weight and Escherichia coli counts. Probiotics improved insulin sensitivity at fasting (QUICKI, 0.013 CI95%0.0-0.03) and during OGTT (ISI, 0.654 CI95%-0.11-1.41). Cytokines, GLP1, and target microbial counts did not vary. Of 25 SCFAs, acetic acid and acetic acid pentyl-ester relative abundance remained stable in the probiotics, while increased in the placebo (p < 0.02). A signature of five butanoic esters identified three clusters, one of them had better glucose responses during probiotics.

CONCLUSION

An 8 weeks treatment with B. breve BR03 and B632 had beneficial effects on insulin sensitivity in youths with obesity. Microbiota functionality could influence metabolic answers to probiotics. Long-term studies to confirm and enrich our findings are justified. Tailored probiotic treatments could be an additional strategy for obesity.

TRIAL REGISTRATION

NCT03261466.

Natural contaminants in bee pollen: DNA metabarcoding as a tool to identify floral sources of pyrrolizidine alkaloids and fungal diversity

The use of bee pollen as a food supplement has increased in recent years as it contains several nutrients and phytochemicals. However, depending on floral composition, bee pollen can be contaminated by pyrrolizidine alkaloids (PAs), PA N-oxides (PANOs) and toxigenic fungi found in plants, which may pose a potential health risk for consumers. Thus, a DNA metabarcoding approach based on internal transcribed spacer 2 (ITS2) region was used to identify the plant sources of 17 PAs/PANOs detected by a validated method in liquid chromatography coupled to mass spectrometry (LC-MS/MS), as well as floral and fungal diversity in 61 bee pollen samples. According to LC-MS/MS analysis, 67% of the samples contained PAs/PANOs with mean concentration of 339 µg/kg. The contamination pattern was characterised by lycopsamine- and senecionine-type PAs/PANOs. PA/PANO-producing plants were identified in 54% of the PA/PANO-contaminated samples analysed by DNA metabarcoding, which also allowed identifying the overall floral and fungal composition of 56 samples. To evaluate the performance of the molecular approach, a subset of 25 samples was analysed by classical palynology. Palynological analysis partially confirmed the results of DNA metabarcoding, which had a better performance in distinguishing pollens of different genera from Asteraceae (76%) and Brassicaceae (88%). However, the molecular analysis did not identify pollens from Castanea, Eucalyptus, Hedera and Salix, which were abundant in 11 samples according to palynology. On the other hand, the molecular analysis allowed identifying several fungal genera in 33 samples, including the toxigenic fungi Alternaria and Aspergillus, which were positively correlated to the plant genus Hypericum. Despite limitations in identifying some pollen types, these preliminary results suggest that the DNA metabarcoding could be applied in a multidisciplinary approach to give a picture of floral and fungal diversity, which can be sources of natural contaminants in bee pollen and would help to control its safety.

High Carotenoid Mutants of Chlorella vulgaris Show Enhanced Biomass Yield under High Irradiance

Microalgae represent a carbon-neutral source of bulk biomass, for extraction of high-value compounds and production of renewable fuels. Due to their high metabolic activity and reproduction rates, species of the genus Chlorella are highly productive when cultivated in photobioreactors. However, wild-type strains show biological limitations making algal bioproducts expensive compared to those extracted from other feedstocks. Such constraints include inhomogeneous light distribution due to high optical density of the culture, and photoinhibition of the surface-exposed cells. Thus, the domestication of algal strains for industry makes it increasingly important to select traits aimed at enhancing light-use efficiency while withstanding excess light stress. Carotenoids have a crucial role in protecting against photooxidative damage and, thus, represent a promising target for algal domestication. We applied chemical mutagenesis to Chlorella vulgaris and selected for enhanced tolerance to the carotenoid biosynthesis inhibitor norflurazon. The NFR (norflurazon-resistant) strains showed an increased carotenoid pool size and enhanced tolerance towards photooxidative stress. Growth under excess light revealed an improved carbon assimilation rate of NFR strains with respect to WT. We conclude that domestication of Chlorella vulgaris, by optimizing both carotenoid/chlorophyll ratio and resistance to photooxidative stress, boosted light-to-biomass conversion efficiency under high light conditions typical of photobioreactors. Comparison with strains previously reported for enhanced tolerance to singlet oxygen, reveals that ROS resistance in Chlorella is promoted by at least two independent mechanisms, only one of which is carotenoid-dependent.

Pangenome analyses of LuxS-coding genes and enzymatic repertoires in cocoa-related lactic acid bacteria

Lactobacillaceae presents potential for interspecific Quorum Sensing (QS) in spontaneous cocoa fermentation, correlated with high abundance of luxS. Three Brazilian isolates from cocoa fermentation were characterized by Whole Genome Sequencing and luxS gene was surveyed in their genomes, in comparison with public databases. They were classified as Lactiplantibacillus plantarum, Limosilactobacillus fermentum and Pediococcus acidilactici. LuxS genes were conserved in core genomes of the novel isolates, but in some non-cocoa related Lactic Acid Bacteria (LAB) it was accessory and plasmid-borne. The conservation and horizontal acquisition of luxS reinforces that QS is determinant for bacterial adaptation in several environments, especially taking into account the luxS has been correlated with modulation of bacteriocin production, stress tolerance and biofilm formation. Therefore, in this paper, new clade and species-specific primers were designed for future application for screening of luxS gene in LAB to evaluate the adaptive potential to diverse food fermentations.

Metabarcoding analysis of gut microbiota of healthy individuals reveals impact of probiotic and maltodextrin consumption

In a previously published double-blind, placebo-controlled study, we showed that probiotics intake exerted a positive effect on sleep quality and a general improvement across time in different aspects of the profile of mood state, like sadness, anger, and fatigue in 33 healthy individuals. The present work investigates the impact of the probiotic product, constituted of Limosilactobacillus fermentum LF16, Lacticaseibacillus rhamnosus LR06, Lactiplantibacillus plantarum LP01 (all former members of Lactobacillus genus), and Bifidobacterium longum 04, on the gut microbiota composition of the same cohort through a metabarcoding analysis. Both the placebo and probiotic treatments had a significant impact on the microbiota composition. Statistical analysis showed that the microbiota of the individuals could be clustered into three groups, or bacteriotypes, at the baseline, and, inherently, bacterial compositions were linked to different responses to probiotic and placebo intakes. Interestingly, L. rhamnosus and L. fermentum were retrieved in the probiotic-treated cohort, while a bifidogenic effect of maltodextrin, used as placebo, was observed. The present study shed light on the importance of defining bacteriotypes to assess the impact of interventions on the gut microbiota and allowed to reveal microbial components which could be related to positive effects (i.e. sleep quality improvement) to be verified in further studies.

The Grapevine E3 Ubiquitin Ligase VriATL156 Confers Resistance against the Downy Mildew Pathogen Plasmopara viticola

Downy mildew, caused by Plasmopara viticola, is one of the most severe diseases of grapevine (Vitis vinifera L.). Genetic resistance is an effective and sustainable control strategy, but major resistance genes (encoding receptors for specific pathogen effectors) introgressed from wild Vitis species, although effective, may be non-durable because the pathogen can evolve to avoid specific recognition. Previous transcriptomic studies in the resistant species Vitis riparia highlighted the activation of signal transduction components during infection. The transfer of such components to V. vinifera might confer less specific and therefore more durable resistance. Here, we describe the generation of transgenic V. vinifera lines constitutively expressing the V. riparia E3 ubiquitin ligase gene VriATL156. Phenotypic and molecular analysis revealed that the transgenic plants were less susceptible to P. viticola than vector-only controls, confirming the role of this E3 ubiquitin ligase in the innate immune response. Two independent transgenic lines were selected for detailed analysis of the resistance phenotype by RNA-Seq and microscopy, revealing the profound reprogramming of transcription to achieve resistance that operates from the earliest stages of pathogen infection. The introduction of VriATL156 into elite grapevine cultivars could therefore provide an effective and sustainable control measure against downy mildew.

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.