Investigating circulating miRNA in transition dairy cows: What miRNAomics tells about metabolic adaptation.

In the current study, we investigated dairy cows' circulating microRNA (miRNA) expression signature during several key time points around calving, to get insights into different aspects of metabolic adaptation. In a trial with 32 dairy cows, plasma samples were collected on days -21, 1, 28, and 63 relative to calving. Individually extracted total RNA was subjected to RNA sequencing using NovaSeq 6,000 (Illumina, CA) on the respective platform of IGA Technology Services, Udine, Italy. MiRDeep2 was used to identify known and novel miRNA according to the miRbase collection. Differentially expressed miRNA (DEM) were assessed at a threshold of fold-change > 1.5 and false discovery rate < 0.05 using the edgeR package. The MiRWalk database was used to predict DEM targets and their associated KEGG pathways. Among a total of 1,692 identified miRNA, 445 known miRNA were included for statistical analysis, of which 84, 59, and 61 DEM were found between days -21 to 1, 1 to 28, and 28 to 63, respectively. These miRNA were annotated to KEGG pathways targeting the insulin, MAPK, Ras, Wnt, Hippo, sphingolipid, T cell receptor, and mTOR signaling pathways. MiRNA-mRNA network analysis identified miRNA as master regulators of the biological process including miR-138, miR-149-5p, miR-2466-3p, miR-214, miR-504, and miR-6523a. This study provided new insights into the miRNA signatures of transition to the lactation period. Calving emerged as a critical time point when miRNA were most affected, while the following period appeared to be recovering from massive parturition changes. The primarily affected pathways were key signaling pathways related to establishing metabolic and immune adaptations.

Identification of a Prognostic Microenvironment-Related Gene Signature in Glioblastoma Patients Treated with Carmustine Wafers.

Despite the state-of-the-art treatment, patients diagnosed with glioblastoma (GBM) have a median overall survival (OS) of 14 months. The insertion of carmustine wafers (CWs) into the resection cavity as adjuvant treatment represents a promising option, although its use has been limited due to contrasting clinical results. Our retrospective evaluation of CW efficacy showed a significant improvement in terms of OS in a subgroup of patients. Given the crucial role of the tumor microenvironment (TME) in GBM progression and response to therapy, we hypothesized that the TME of patients who benefited from CW could have different properties compared to that of patients who did not show any advantage. Using an in vitro model of the glioma microenvironment, represented by glioma-associated-stem cells (GASC), we performed a transcriptomic analysis of GASC isolated from tumors of patients responsive and not responsive to CW to identify differentially expressed genes. We found different transcriptomic profiles, and we identified four genes, specifically down-regulated in GASC isolated from long-term survivors, correlated with clinical data deposited in the TCGA-GBM dataset. Our results highlight that studying the in vitro properties of patient-specific glioma microenvironments can help to identify molecular determinants potentially prognostic for patients treated with CW.

Genotypic and tissue-specific variation of Populus nigra transcriptome profiles in response to drought.

Climate change is one of the most important challenges for mankind in the far and near future. In this regard, sustainable production of woody crops on marginal land with low water availability is a major challenge to tackle. This dataset is part of an experiment, in which we exposed three genetically differentiated genotypes of Populus nigra originating from contrasting natural habitats to gradually increasing moderate drought. RNA sequencing was performed on fine roots, developing xylem and leaves of those three genotypes under control and moderate drought conditions in order to get a comprehensive dataset on the transcriptional changes at the whole plant level under water limiting conditions. This dataset has already provided insight in the transcriptional control of saccharification potential of the three Populus genotypes under drought conditions and we suggest that our data will be valuable for further in-depth analysis regarding candidate gene identification or, on a bigger scale, for meta-transcriptome analysis.

First whole genome characterization of porcine astrovirus detected in swine faeces in Italy.

Porcine astroviruses (PoAstV) are found in the gastrointestinal tract of healthy and diseased pigs worldwide. However, their role in causing enteric disease in pigs and other animals has not been elucidated. In the present report, we describe for the first time in Italy, the identification and genetic characterization, through whole genome sequencing, of a PoAstV2 in pigs in Northeast Italy in 2015. This instance is the first detection of PoAstV2 in pigs in Italy. The phylogenetic analysis of the complete ORF2 segment highlights the high similarity of this virus to those circulating that same year in Japan. There are very few full astrovirus genomes available, and the present data represent an important contribution towards a better understanding of the characteristics and evolution of these viruses.

Alkaloid chemodiversity in Mandragora spp. is associated with loss-of-functionality of MoH6H, a hyoscyamine 6ß-hydroxylase gene.

Mandrakes (Mandragora spp., Solanaceae) are known to contain tropane alkaloids and have been used since antiquity in traditional medicine. Tropane alkaloids such as scopolamine and hyoscyamine are used in modern medicine to treat pain, motion sickness, as eye pupil dilators and antidotes against organo-phosphate poisoning. Hyoscyamine is converted to 6ß-hydroxyhyoscyamine (anisodamine) and scopolamine by hyoscyamine 6ß-hydroxylase (H6H), a 2-oxoglutarate dependent dioxygenase. We describe here a marked chemo-diversity in the tropane alkaloid content in Mandragora spp. M. officinarum and M. turcomanica lack anisodamine and scopolamine but display up to 10 fold higher hyoscyamine levels as compared with M. autumnalis. Transcriptomic analyses revealed that H6H is highly conserved among scopolamine-producing Solanaceae. MoH6H present in M. officinarum differs in several amino acid residues including a homozygotic mutation in the substrate binding region of the protein and its prevalence among accessions was confirmed by Cleaved-Amplified-Polymorphic-Sequence analyses. Functional expression revealed that MaH6H, a gene isolated from M. autumnalis encodes an active H6H enzyme while the MoH6H sequence isolated from M. officinarum was functionally inactive. A single G to T mutation in nucleotide 663 of MoH6H is associated with the lack of anisodamine and scopolamine in M. officinalis.

Genomic tools for durum wheat breeding: de novo assembly of Svevo transcriptome and SNP discovery in elite germplasm.

BACKGROUND: The tetraploid durum wheat (Triticum turgidum L. ssp. durum Desf. Husnot) is an important crop which provides the raw material for pasta production and a valuable source of genetic diversity for breeding hexaploid wheat (Triticum aestivum L.). Future breeding efforts to enhance yield potential and climate resilience will increasingly rely on genomics-based approaches to identify and select beneficial alleles. A deeper characterisation of the molecular and functional diversity of the durum wheat transcriptome will be instrumental to more effectively harness its genetic diversity. RESULTS: We report on the de novo transcriptome assembly of durum wheat cultivar 'Svevo'. The transcriptome of four tissues/organs (shoots and roots at the seedling stage, reproductive organs and developing grains) was assembled de novo, yielding 180,108 contigs, with a N50 length of 1121 bp and mean contig length of 883 bp. Alignment against the transcriptome of nine plant species identified 43% of transcripts with homology to at least one reference transcriptome. The functional annotation was completed by means of a combination of complementary software. The presence of differential expression between the A- and B-homoeolog copies of the durum wheat tetraploid genome was ascertained by phase reconstruction of polymorphic sites based on the T. urartu transcripts and inferring homoeolog-specific sequences. We observed greater expression divergence between A and B homoeologs in grains rather than in leaves and roots. The transcriptomes of 13 durum wheat cultivars spanning the breeding period from 1969 to 2005 were analysed for SNP diversity, leading to 95,358 non-rare, hemi-SNPs shared among two or more cultivars and 33,747 locus-specific (diploid inheritance) SNPs. CONCLUSIONS: Our study updates and expands the de novo transcriptome reference assembly available for durum wheat. Out of 180,108 assembled transcripts, 13,636 were specific to the Svevo cultivar as compared to the only other reference transcriptome available for durum, thus contributing to the identification of the tetraploid wheat pan-transcriptome. Additionally, the analysis of 13 historically relevant hallmark varieties produced a SNP dataset that could successfully validate the genotyping in tetraploid wheat and provide a valuable resource for genomics-assisted breeding of both tetraploid and hexaploid wheats.

Genes and gene clusters related to genotype and drought-induced variation in saccharification potential, lignin content and wood anatomical traits in Populus nigra.

Wood is a renewable resource that can be employed for the production of second generation biofuels by enzymatic saccharification and subsequent fermentation. Knowledge on how the saccharification potential is affected by genotype-related variation of wood traits and drought is scarce. Here, we used three Populus nigra L. genotypes from habitats differing in water availability to (i) investigate the relationships between wood anatomy, lignin content and saccharification and (ii) identify genes and co-expressed gene clusters related to genotype and drought-induced variation in wood traits and saccharification potential. The three poplar genotypes differed in wood anatomy, lignin content and saccharification potential. Drought resulted in reduced cambial activity, decreased vessel and fiber lumina, and increased the saccharification potential. The saccharification potential was unrelated to lignin content as well as to most wood anatomical traits. RNA sequencing of the developing xylem revealed that 1.5% of the analyzed genes were differentially expressed in response to drought, while 67% differed among the genotypes. Weighted gene correlation network analysis identified modules of co-expressed genes correlated with saccharification potential. These modules were enriched in gene ontology terms related to cell wall polysaccharide biosynthesis and modification and vesicle transport, but not to lignin biosynthesis. Among the most strongly saccharification-correlated genes, those with regulatory functions, especially kinases, were prominent. We further identified transcription factors whose transcript abundances differed among genotypes, and which were co-regulated with genes for biosynthesis and modifications of hemicelluloses and pectin. Overall, our study suggests that the regulation of pectin and hemicellulose metabolism is a promising target for improving wood quality of second generation bioenergy crops. The causal relationship of the identified genes and pathways with saccharification potential needs to be validated in further experiments.

The build-up of osmotic stress responses within the growing root apex using kinematics and RNA-sequencing.

Molecular regulation of growth must include spatial and temporal coupling of cell production and cell expansion. The underlying mechanisms, especially under environmental challenge, remain obscure. Spatial patterns of cell processes make the root apex well suited to deciphering stress signaling pathways, and to investigating both processes. Kinematics and RNA-sequencing were used to analyze the immediate growth response of hydroponically grown Populus nigra cuttings submitted to osmotic stress. About 7400 genes and unannotated transcriptionally active regions were differentially expressed between the division and elongation zones. Following the onset of stress, growth decreased sharply, probably due to mechanical effects, before recovering partially. Stress impaired cell expansion over the apex, progressively shortened the elongation zone, and reduced the cell production rate. Changes in gene expression revealed that growth reduction was mediated by a shift in hormone homeostasis. Osmotic stress rapidly elicited auxin, ethylene, and abscisic acid. When growth restabilized, transcriptome remodeling became complex and zone specific, with the deployment of hormone signaling cascades, transcriptional regulators, and stress-responsive genes. Most transcriptional regulations fit growth reduction, but stress also promoted expression of some growth effectors, including aquaporins and expansins Together, osmotic stress interfered with growth by activating regulatory proteins rather than by repressing the machinery of expansive growth.

A knowledge base for Vitis vinifera functional analysis.

BACKGROUND: Vitis vinifera (Grapevine) is the most important fruit species in the modern world. Wine and table grapes sales contribute significantly to the economy of major wine producing countries. The most relevant goals in wine production concern quality and safety. In order to significantly improve the achievement of these objectives and to gain biological knowledge about cultivars, a genomic approach is the most reliable strategy. The recent grapevine genome sequencing offers the opportunity to study the potential roles of genes and microRNAs in fruit maturation and other physiological and pathological processes. Although several systems allowing the analysis of plant genomes have been reported, none of them has been designed specifically for the functional analysis of grapevine genomes of cultivars under environmental stress in connection with microRNA data. DESCRIPTION: Here we introduce a novel knowledge base, called BIOWINE, designed for the functional analysis of Vitis vinifera genomes of cultivars present in Sicily. The system allows the analysis of RNA-seq experiments of two different cultivars, namely Nero d'Avola and Nerello Mascalese. Samples were taken under different climatic conditions of phenological phases, diseases, and geographic locations. The BIOWINE web interface is equipped with data analysis modules for grapevine genomes. In particular users may analyze the current genome assembly together with the RNA-seq data through a customized version of GBrowse. The web interface allows users to perform gene set enrichment by exploiting third-party databases. CONCLUSIONS: BIOWINE is a knowledge base implementing a set of bioinformatics tools for the analysis of grapevine genomes. The system aims to increase our understanding of the grapevine varieties and species of Sicilian products focusing on adaptability to different climatic conditions, phenological phases, diseases, and geographic locations.

A high-density, SNP-based consensus map of tetraploid wheat as a bridge to integrate durum and bread wheat genomics and breeding.

Consensus linkage maps are important tools in crop genomics. We have assembled a high-density tetraploid wheat consensus map by integrating 13 data sets from independent biparental populations involving durum wheat cultivars (Triticum turgidum ssp. durum), cultivated emmer (T. turgidum ssp. dicoccum) and their ancestor (wild emmer, T. turgidum ssp. dicoccoides). The consensus map harboured 30 144 markers (including 26 626 SNPs and 791 SSRs) half of which were present in at least two component maps. The final map spanned 2631 cM of all 14 durum wheat chromosomes and, differently from the individual component maps, all markers fell within the 14 linkage groups. Marker density per genetic distance unit peaked at centromeric regions, likely due to a combination of low recombination rate in the centromeric regions and even gene distribution along the chromosomes. Comparisons with bread wheat indicated fewer regions with recombination suppression, making this consensus map valuable for mapping in the A and B genomes of both durum and bread wheat. Sequence similarity analysis allowed us to relate mapped gene-derived SNPs to chromosome-specific transcripts. Dense patterns of homeologous relationships have been established between the A- and B-genome maps and between nonsyntenic homeologous chromosome regions as well, the latter tracing to ancient translocation events. The gene-based homeologous relationships are valuable to infer the map location of homeologs of target loci/QTLs. Because most SNP and SSR markers were previously mapped in bread wheat, this consensus map will facilitate a more effective integration and exploitation of genes and QTL for wheat breeding purposes.

Genetic dissection of barley morphology and development.

Since the early 20th century, barley (Hordeum vulgare) has been a model for investigating the effects of physical and chemical mutagens and for exploring the potential of mutation breeding in crop improvement. As a consequence, extensive and well-characterized collections of morphological and developmental mutants have been assembled that represent a valuable resource for exploring a wide range of complex and fundamental biological processes. We constructed a collection of 881 backcrossed lines containing mutant alleles that induce a majority of the morphological and developmental variation described in this species. After genotyping these lines with up to 3,072 single nucleotide polymorphisms, comparison to their recurrent parent defined the genetic location of 426 mutant alleles to chromosomal segments, each representing on average <3% of the barley genetic map. We show how the gene content in these segments can be predicted through conservation of synteny with model cereal genomes, providing a route to rapid gene identification.

Exploiting induced variation to dissect quantitative traits in barley.

The identification of genes underlying complex quantitative traits such as grain yield by means of conventional genetic analysis (positional cloning) requires the development of several large mapping populations. However, it is possible that phenotypically related, but more extreme, allelic variants generated by mutational studies could provide a means for more efficient cloning of QTLs (quantitative trait loci). In barley (Hordeum vulgare), with the development of high-throughput genome analysis tools, efficient genome-wide identification of genetic loci harbouring mutant alleles has recently become possible. Genotypic data from NILs (near-isogenic lines) that carry induced or natural variants of genes that control aspects of plant development can be compared with the location of QTLs to potentially identify candidate genes for development--related traits such as grain yield. As yield itself can be divided into a number of allometric component traits such as tillers per plant, kernels per spike and kernel size, mutant alleles that both affect these traits and are located within the confidence intervals for major yield QTLs may represent extreme variants of the underlying genes. In addition, the development of detailed comparative genomic models based on the alignment of a high-density barley gene map with the rice and sorghum physical maps, has enabled an informed prioritization of 'known function' genes as candidates for both QTLs and induced mutant genes.