SAMHD1-deficient fibroblasts from Aicardi-Goutières Syndrome patients can escape senescence and accumulate mutations.

In mammalian cells, the catabolic activity of the dNTP triphosphohydrolase SAMHD1 sets the balance and concentration of the four dNTPs. Deficiency of SAMHD1 leads to unequally increased pools and marked dNTP imbalance. Imbalanced dNTP pools increase mutation frequency in cancer cells, but it is not known if the SAMHD1-induced dNTP imbalance favors accumulation of somatic mutations in non-transformed cells. Here, we have investigated how fibroblasts from Aicardi-Goutières Syndrome (AGS) patients with mutated SAMHD1 react to the constitutive pool imbalance characterized by a huge dGTP pool. We focused on the effects on dNTP pools, cell cycle progression, dynamics and fidelity of DNA replication, and efficiency of UV-induced DNA repair. AGS fibroblasts entered senescence prematurely or upregulated genes involved in G1/S transition and DNA replication. The normally growing AGS cells exhibited unchanged DNA replication dynamics and, when quiescent, faster rate of excision repair of UV-induced DNA damages. To investigate whether the lack of SAMHD1 affects DNA replication fidelity, we compared de novo mutations in AGS and WT cells by exome next-generation sequencing. Somatic variant analysis indicated a mutator phenotype suggesting that SAMHD1 is a caretaker gene whose deficiency is per se mutagenic, promoting genome instability in non-transformed cells.

Optimizing Microbiota Profiles for Athletes.

Gut microbiome influences athletes' physiology, but because of the complexity of sport performance and the great intervariability of microbiome features, it is not reasonable to define a single healthy microbiota profile for athletes. We suggest the use of specific meta-omics analysis coupled with innovative computational systems to uncover the hidden association between microbes and athlete's physiology and predict personalized recommendation.

Revealing metabolic mechanisms of interaction in the anaerobic digestion microbiome by flux balance analysis.

Anaerobic digestion is a key biological process for renewable energy, yet the mechanistic knowledge on its hidden microbial dynamics is still limited. The present work charted the interaction network in the anaerobic digestion microbiome via the full characterization of pairwise interactions and the associated metabolite exchanges. To this goal, a novel collection of 836 genome-scale metabolic models was built to represent the functional capabilities of bacteria and archaea species derived from genome-centric metagenomics. Dominant microbes were shown to prefer mutualistic, parasitic and commensalistic interactions over neutralism, amensalism and competition, and are more likely to behave as metabolite importers and profiteers of the coexistence. Additionally, external hydrogen injection positively influences microbiome dynamics by promoting commensalism over amensalism. Finally, exchanges of glucogenic amino acids were shown to overcome auxotrophies caused by an incomplete tricarboxylic acid cycle. Our novel strategy predicted the most favourable growth conditions for the microbes, overall suggesting strategies to increasing the biogas production efficiency. In principle, this approach could also be applied to microbial populations of biomedical importance, such as the gut microbiome, to allow a broad inspection of the microbial interplays.

Genetic Analyses in Dent Disease and Characterization of CLCN5 Mutations in Kidney Biopsies.

Dent disease (DD), an X-linked renal tubulopathy, is mainly caused by loss-of-function mutations in CLCN5 (DD1) and OCRL genes. CLCN5 encodes the ClC-5 antiporter that in proximal tubules (PT) participates in the receptor-mediated endocytosis of low molecular weight proteins. Few studies have analyzed the PT expression of ClC-5 and of megalin and cubilin receptors in DD1 kidney biopsies. About 25% of DD cases lack mutations in either CLCN5 or OCRL genes (DD3), and no other disease genes have been discovered so far. Sanger sequencing was used for CLCN5 gene analysis in 158 unrelated males clinically suspected of having DD. The tubular expression of ClC-5, megalin, and cubilin was assessed by immunolabeling in 10 DD1 kidney biopsies. Whole exome sequencing (WES) was performed in eight DD3 patients. Twenty-three novel CLCN5 mutations were identified. ClC-5, megalin, and cubilin were significantly lower in DD1 than in control biopsies. The tubular expression of ClC-5 when detected was irrespective of the type of mutation. In four DD3 patients, WES revealed 12 potentially pathogenic variants in three novel genes (SLC17A1, SLC9A3, and PDZK1), and in three genes known to be associated with monogenic forms of renal proximal tubulopathies (SLC3A, LRP2, and CUBN). The supposed third Dent disease-causing gene was not discovered.

Scuba: scalable kernel-based gene prioritization.

BACKGROUND: The uncovering of genes linked to human diseases is a pressing challenge in molecular biology and precision medicine. This task is often hindered by the large number of candidate genes and by the heterogeneity of the available information. Computational methods for the prioritization of candidate genes can help to cope with these problems. In particular, kernel-based methods are a powerful resource for the integration of heterogeneous biological knowledge, however, their practical implementation is often precluded by their limited scalability. RESULTS: We propose Scuba, a scalable kernel-based method for gene prioritization. It implements a novel multiple kernel learning approach, based on a semi-supervised perspective and on the optimization of the margin distribution. Scuba is optimized to cope with strongly unbalanced settings where known disease genes are few and large scale predictions are required. Importantly, it is able to efficiently deal both with a large amount of candidate genes and with an arbitrary number of data sources. As a direct consequence of scalability, Scuba integrates also a new efficient strategy to select optimal kernel parameters for each data source. We performed cross-validation experiments and simulated a realistic usage setting, showing that Scuba outperforms a wide range of state-of-the-art methods. CONCLUSIONS: Scuba achieves state-of-the-art performance and has enhanced scalability compared to existing kernel-based approaches for genomic data. This method can be useful to prioritize candidate genes, particularly when their number is large or when input data is highly heterogeneous. The code is freely available at https://github.com/gzampieri/Scuba .

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.

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.

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.

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.

Transcriptome analysis of the regenerating tail vs. the scarring limb in lizard reveals pathways leading to successful vs. unsuccessful organ regeneration in amniotes.

BACKGROUND: Lizards are amniotes regenerating the tail but not the limb, and no information on their different gene expression is available. RESULTS: Transcriptomes of regenerating tail and limb blastemas show differences in gene expression between the two organs. In tail blastemal, snoRNAs and Wnt signals appear up-regulated probably in association with the apical epidermal peg (AEP), an epithelial region that sustains tail regeneration but is absent in the limb. A balance between pro-oncogenes and tumor suppressors is likely present in tail blastema allowing a regulated proliferation. Small collagens, protease inhibitors, embryonic keratins are up-regulated in the regenerating tail blastema but not in the limb where Wnt inhibitors, inflammation-immune and extracellular matrix proteins depress cell growth. CONCLUSIONS: The AEP and the spinal cord in the tail maintains Wnt and fibroblast growth signaling that stimulate blastema cell proliferation and growth while these signals are absent in the limb as a consequence of the intense inflammation. Regeneration of amniote appendages requires a control of cell proliferation and inflammatory-immune reactions to form an apical epidermal cap. Genes that control cell proliferation and inflammation, addressing regeneration and not tumor formation in the tail and scarring in the limb are discussed for future studies. Developmental Dynamics 246:116-134, 2017. © 2016 Wiley Periodicals, Inc.

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/ .

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.

Brain RNA-Seq Profiling of the Mucopolysaccharidosis Type II Mouse Model.

Lysosomal storage disorders (LSDs) are a group of about 50 genetic metabolic disorders, mainly affecting children, sharing the inability to degrade specific endolysosomal substrates. This results in failure of cellular functions in many organs, including brain that in most patients may go through progressive neurodegeneration. In this study, we analyzed the brain of the mouse model for Hunter syndrome, a LSD mostly presenting with neurological involvement. Whole transcriptome analysis of the cerebral cortex and midbrain/diencephalon/hippocampus areas was performed through RNA-seq. Genes known to be involved in several neurological functions showed a significant differential expression in the animal model for the disease compared to wild type. Among the pathways altered in both areas, axon guidance, calcium homeostasis, synapse and neuroactive ligand-receptor interaction, circadian rhythm, neuroinflammation and Wnt signaling were the most significant. Application of RNA sequencing to dissect pathogenic alterations of complex syndromes allows to photograph perturbations, both determining and determined by these disorders, which could simultaneously occur in several metabolic and biochemical pathways. Results also emphasize the common, altered pathways between neurodegenerative disorders affecting elderly and those associated with pediatric diseases of genetic origin, perhaps pointing out a general common course for neurodegeneration, independent from the primary triggering cause.

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.

Differential expression and localization of Ankrd2 isoforms in human skeletal and cardiac muscles.

Four human Ankrd2 transcripts, reported in the Ensembl database, code for distinct protein isoforms (360, 333, 327 and 300 aa), and so far, their existence, specific expression and localization patterns have not been studied in detail. Ankrd2 is preferentially expressed in the slow fibers of skeletal muscle. It is found in both the nuclei and the cytoplasm of skeletal muscle cells, and its localization is prone to change during differentiation and upon stress. Ankrd2 has also been detected in the heart, in ventricular cardiomyocytes and in the intercalated disks (ICDs). The main objective of this study was to distinguish between the Ankrd2 isoforms and to determine the contribution of each one to the general profile of Ankrd2 expression in striated muscles. We demonstrated that the known expression and localization pattern of Ankrd2 in striated muscle can be attributed to the isoform of 333 aa which is dominant in both tissues, while the designated cardiac and canonical isoform of 360 aa was less expressed in both tissues. The 360 aa isoform has a distinct nuclear localization in human skeletal muscle, as well as in primary myoblasts and myotubes. In contrast to the isoform of 333 aa, it was not preferentially expressed in slow fibers and not localized to the ICDs of human cardiomyocytes. Regulation of the expression of both isoforms is achieved at the transcriptional level. Our results set the stage for investigation of the specific functions and interactions of the Ankrd2 isoforms in healthy and diseased human striated muscles.

Profiling of skeletal muscle Ankrd2 protein in human cardiac tissue and neonatal rat cardiomyocytes.

Muscle-specific mechanosensors Ankrd2/Arpp (ankyrin repeat protein 2) and Ankrd1/CARP (cardiac ankyrin repeat protein) have an important role in transcriptional regulation, myofibrillar assembly, cardiogenesis and myogenesis. In skeletal muscle myofibrils, Ankrd2 has a structural role as a component of a titin associated stretch-sensing complex, while in the nucleus it exerts regulatory function as transcriptional co-factor. It is also involved in myogenic differentiation and coordination of myoblast proliferation. Although expressed in the heart, the role of Ankrd2 in the cardiac muscle is completely unknown. Recently, we have shown that hypertrophic and dilated cardiomyopathy pathways are altered upon Ankrd2 silencing suggesting the importance of this protein in cardiac tissue. Here we provide the underlying basis for the functional investigation of Ankrd2 in the heart. We confirmed reduced Ankrd2 expression levels in human heart in comparison with Ankrd1 using RNAseq and Western blot. For the first time we demonstrated that, apart from the sarcomere and nucleus, both proteins are localized to the intercalated disks of human cardiomyocytes. We further tested the expression and localization of endogenous Ankrd2 in rat neonatal cardiomyocytes, a well-established model for studying cardiac-specific proteins. Ankrd2 was found to be expressed in both the cytoplasm and nucleus, independently from maturation status of cardiomyocytes. In contrast to Ankrd1, it is not responsive to the cardiotoxic drug Doxorubicin, suggesting that different mechanisms govern their expression in cardiac cells.

Comprehensive transcript profiling of two grapevine rootstock genotypes contrasting in drought susceptibility links the phenylpropanoid pathway to enhanced tolerance.

In light of ongoing climate changes in wine-growing regions, the selection of drought-tolerant rootstocks is becoming a crucial factor for developing a sustainable viticulture. In this study, M4, a new rootstock genotype that shows tolerance to drought, was compared from a genomic and transcriptomic point of view with the less drought-tolerant genotype 101.14. The root and leaf transcriptome of both 101.14 and the M4 rootstock genotype was analysed, following exposure to progressive drought conditions. Multifactorial analyses indicated that stress treatment represents the main factor driving differential gene expression in roots, whereas in leaves the genotype is the prominent factor. Upon stress, M4 roots and leaves showed a higher induction of resveratrol and flavonoid biosynthetic genes, respectively. The higher expression of VvSTS genes in M4, confirmed by the accumulation of higher levels of resveratrol in M4 roots compared with 101.14, was coupled to an up-regulation of several VvWRKY transcription factors. Interestingly, VvSTS promoter analyses performed on both the resequenced genomes highlighted a significantly higher number of W-BOX elements in the tolerant genotype. It is proposed that the elevated synthesis of resveratrol in M4 roots upon water stress could enhance the plant's ability to cope with the oxidative stress usually associated with water deficit.

A deep survey of alternative splicing in grape reveals changes in the splicing machinery related to tissue, stress condition and genotype.

BACKGROUND: Alternative splicing (AS) significantly enhances transcriptome complexity. It is differentially regulated in a wide variety of cell types and plays a role in several cellular processes. Here we describe a detailed survey of alternative splicing in grape based on 124 SOLiD RNAseq analyses from different tissues, stress conditions and genotypes. RESULTS: We used the RNAseq data to update the existing grape gene prediction with 2,258 new coding genes and 3,336 putative long non-coding RNAs. Several gene structures have been improved and alternative splicing was described for about 30% of the genes. A link between AS and miRNAs was shown in 139 genes where we found that AS affects the miRNA target site. A quantitative analysis of the isoforms indicated that most of the spliced genes have one major isoform and tend to simultaneously co-express a low number of isoforms, typically two, with intron retention being the most frequent alternative splicing event. CONCLUSIONS: As described in Arabidopsis, also grape displays a marked AS tissue-specificity, while stress conditions produce splicing changes to a minor extent. Surprisingly, some distinctive splicing features were also observed between genotypes. This was further supported by the observation that the panel of Serine/Arginine-rich splicing factors show a few, but very marked differences between genotypes. The finding that a part the splicing machinery can change in closely related organisms can lead to some interesting hypotheses for evolutionary adaptation, that could be particularly relevant in the response to sudden and strong selective pressures.

PASS-bis: a bisulfite aligner suitable for whole methylome analysis of Illumina and SOLiD reads.

SUMMARY: The sequencing of bisulfite-treated DNA (Bi-Seq) is becoming a gold standard for methylation studies. The mapping of Bi-Seq reads is complex and requires special alignment algorithms. This problem is particularly relevant for SOLiD color space, where the bisulfite conversion C/T changes two adjacent colors into 16 possible combinations. Here, we present an algorithm that efficiently aligns Bi-Seq reads obtained either from SOLiD or Illumina. An accompanying methylation-caller program creates a genomic view of methylated and unmethylated Cs on both DNA strands. AVAILABILITY AND IMPLEMENTATION: The algorithm has been implemented as an option of the program PASS, freely available at http://pass.cribi.unipd.it. CONTACT: pass@cribi.unipd.it SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Oxidative stress response and nitrogen utilization are strongly variable in Saccharomyces cerevisiae wine strains with different fermentation performances.

We used RNA-sequencing (RNA-seq) to analyze the expression profile of four vineyard strains of Saccharomyces cerevisiae having different fermentation performances. The expression profiles obtained in two steps of the fermentation process were compared with those obtained for the industrial wine strain EC1118 and for the laboratory strain S288c. The two strains with low fermentation efficiency, namely, S288c and the vineyard strain R103, exhibited markedly different expression profiles when compared to the other four strains. We also found that the vineyard strains P283 and P301 are characterized by a high expression of the transcription factor Met32p in the first step of the fermentation. Met32p, in coordination with the Hap4p transcription factor, determined the over-expression of the genes involved in the respiration processes, in the response to oxidative stress and in the sulfur amino acids biosynthesis. These combined actions are likely to increase the level of antioxidants whose protective effect could contribute to improve the fermentation process. Gene expression and phenotypic data revealed that the vineyard strain P301 has low nitrogen utilization in comparison to the other wine strains, combined with high fermentation efficiency. Analysis of the genes involved in fermentation stress response revealed a lower expression in strains characterized by low fermentation efficiency, particularly in the first fermentation phase. These findings evidenced the high variability of transcriptional profiles among different wine yeast strains and clarify their connection with complex phenotypic traits, such as the fermentation efficiency and the nitrogen sources utilization.