LOXL2-mediated chromatin compaction is required to maintain the oncogenic properties of triple-negative breast cancer cells

Oxidation of histone H3 at lysine 4 (H3K4ox) is catalyzed by lysyl oxidase homolog 2 (LOXL2). This histone modification is enriched in heterochromatin in triple-negative breast cancer (TNBC) cells and has been linked to the maintenance of compacted chromatin. However, the molecular mechanism underlying this maintenance is still unknown. Here, we show that LOXL2 interacts with RuvB-Like 1 (RUVBL1), RuvB-Like 2 (RUVBL2), Actin-like protein 6A (ACTL6A), and DNA methyltransferase 1associated protein 1 (DMAP1), a complex involved in the incorporation of the histone variant H2A.Z. Our experiments indicate that this interaction and the active form of RUVBL2 are required to maintain LOXL2-dependent chromatin compaction. Genome-wide experiments showed that H2A.Z, RUVBL2, and H3K4ox colocalize in heterochromatin regions. In the absence of LOXL2 or RUVBL2, global levels of the heterochromatin histone mark H3K9me3 were strongly reduced, and the ATAC-seq signal in the H3K9me3 regions was increased. Finally, we observed that the interplay between these series of events is required to maintain H3K4ox-enriched heterochromatin regions, which in turn is key for maintaining the oncogenic properties of the TNBC cell line tested (MDA-MB-231).

Distinct structural variants and repeat landscape shape the genomes of the ancient grapes Aglianico and Falanghina

Mounting evidence recognizes structural variations (SVs) and repetitive DNA sequences as crucial players in shaping the existing grape phenotypic diversity at intra- and inter-species levels. To deepen our understanding on the abundance, diversity, and distribution of SVs and repetitive DNAs, including transposable elements (TEs) and tandemly repeated satellite DNA (satDNAs), we re-sequenced the genomes of the ancient grapes Aglianico and Falanghina. The analysis of large copy number variants (CNVs) detected candidate polymorphic genes that are involved in the enological features of these varieties. In a comparative analysis of Aglianico and Falanghina sequences with 21 publicly available genomes of cultivated grapes, we provided a genome-wide annotation of grape TEs at the lineage level. We disclosed that at least two main clusters of grape cultivars could be identified based on the TEs content. Multiple TEs families appeared either significantly enriched or depleted. In addition, in silico and cytological analyses provided evidence for a diverse chromosomal distribution of several satellite repeats between Aglianico, Falanghina, and other grapes. Overall, our data further improved our understanding of the intricate grape diversity held by two Italian traditional varieties, unveiling a pool of unique candidate genes never so far exploited in breeding for improved fruit quality.

The Chloroplast Genome of Endive (Cichorium endivia L.): Cultivar Structural Variants and Transcriptome Responses to Stress Due to Rain Extreme Events

The chloroplast (cp) genome diversity has been used in phylogeny studies, breeding, and variety protection, and its expression has been shown to play a role in stress response. Smooth- and curly-leafed endives (Cichorium endivia var. latifolium and var. crispum) are of nutritional and economic importance and are the target of ever-changing breeding programmes. A reference cp genome sequence was assembled and annotated (cultivar 'Confiance'), which was 152,809 base pairs long, organized into the angiosperm-typical quadripartite structure, harboring two inverted repeats separated by the large- and short- single copy regions. The annotation included 136 genes, 90 protein-coding genes, 38 transfer, and 8 ribosomal RNAs and the sequence generated a distinct phyletic group within Asteraceae with the well-separated C. endivia and intybus species. SSR variants within the reference genome were mostly of tri-nucleotide type, and the cytosine to uracil (C/U) RNA editing recurred. The cp genome was nearly fully transcribed, hence sequence polymorphism was investigated by RNA-Seq of seven cultivars, and the SNP number was higher in smooth- than curly-leafed ones. All cultivars maintained C/U changes in identical positions, suggesting that RNA editing patterns were conserved; most cultivars shared SNPs of moderate impact on protein changes in the ndhD, ndhA, and psbF genes, suggesting that their variability may have a potential role in adaptive response. The cp transcriptome expression was investigated in leaves of plants affected by pre-harvest rainfall and rainfall excess plus waterlogging events characterized by production loss, compared to those of a cycle not affected by extreme rainfall. Overall, the analyses evidenced stress- and cultivar-specific responses, and further revealed that genes of the Cytochrome b6/f, and PSI-PSII systems were commonly affected and likely to be among major targets of extreme rain-related stress.

Transcriptomic landscape of tomato traditional long shelf-life landraces under low water regimes

'Corbarino' (COR) and 'Lucariello' (LUC) belong to the family of Mediterranean long shelf-life tomato landraces, producing high quality fruits under low water input cultivation regime in their traditional cultivation area. Understanding the morpho-physiological and molecular details of the peculiar drought stress tolerance of these two genotypes may be key to their valorization as breeding material. RNA sequencing of leaf samples of COR and LUC subjected to drought stress by water withholding in a semi-controlled greenhouse identified 3089 and 2135 differentially expressed genes respectively. These included COR- and LUC-specific annotated genes, as well as genes containing single nucleotide polymorphisms as compared to reference genome. Enriched Gene Ontology categories showed that categories such as response to water, oxidoreductase activity, nucleotide salvation and lipid biosynthesis-related processes were enriched among up-regulated DEGs. By contrast, growth and photosynthesis related genes were down-regulated after drought stress, consistent with leaf gas exchange and biomass accumulation measurements. Genes encoding cell wall degrading enzymes of the pectinase family were also down-regulated in drought stress conditions and upregulated in rewatering, indicating that cell wall composition/hardness is important for drought stress responses. Globally our results contribute to understanding the transcriptomic and physiological responses of representative tomato genotypes from Southern Italy, highlighting a promising set of genes to be investigated to improve tomato tolerance to drought.

Mitochondrial DNA is a target of HBV integration

Hepatitis B virus (HBV) may integrate into the genome of infected cells and contribute to hepatocarcinogenesis. However, the role of HBV integration in hepatocellular carcinoma (HCC) development remains unclear. In this study, we apply a high-throughput HBV integration sequencing approach that allows sensitive identification of HBV integration sites and enumeration of integration clones. We identify 3339 HBV integration sites in paired tumour and non-tumour tissue samples from 7 patients with HCC. We detect 2107 clonally expanded integrations (1817 in tumour and 290 in non-tumour tissues), and a significant enrichment of clonal HBV integrations in mitochondrial DNA (mtDNA) preferentially occurring in the oxidative phosphorylation genes (OXPHOS) and D-loop region. We also find that HBV RNA sequences are imported into the mitochondria of hepatoma cells with the involvement of polynucleotide phosphorylase (PNPASE), and that HBV RNA might have a role in the process of HBV integration into mtDNA. Our results suggest a potential mechanism by which HBV integration may contribute to HCC development.

Effect of a Co-Feed Liquid Whey-Integrated Diet on Crossbred Pigs' Fecal Microbiota

This study assessed the potential effect of a co-feed liquid whey-integrated diet on the fecal microbiota of 14 crossbred pigs. The experimental design was as follows: seven pigs were in the control group, fed with a control feed, and seven were in the experimental group, fed with the same control feed supplemented daily with liquid whey. The collection of fecal samples was conducted on each animal before the dietary treatment (T0) and one (T1), and two (T2) months after the beginning of the co-feed integration. In addition, blood samples were collected from each pig at the same time points in order to evaluate the physiological parameters. Taxonomic analysis showed a bacterial community dominated by Firmicutes, Bacteroidetes, Spirochaetes, and Proteobacteria phyla that populated the crossbred pig feces. The diversity metrics suggested that the co-feed supplementation affected some alpha diversity indexes of the fecal microbiota. In addition, the differential abundance analysis at the genus level revealed significant differences for various genera, suggesting that the liquid whey supplementation potentially influenced a part of the bacterial community over time. Spearman's correlations revealed that the differential abundant genera identified are positively or negatively correlated with the physiological parameters.

Transcriptomic profiles of the ruminal wall in Italian Mediterranean dairy buffaloes fed green forage

BACKGROUND

Green feed diet in ruminants exerts a beneficial effect on rumen metabolism and enhances the content of milk nutraceutical quality. At present, a comprehensive analysis focused on the identification of genes, and therefore, biological processes modulated by the green feed in buffalo rumen has never been reported. We performed RNA-sequencing in the rumen of buffaloes fed a total mixed ration (TMR) + the inclusion of 30% of ryegrass green feed (treated) or TMR (control), and identified differentially expressed genes (DEGs) using EdgeR and NOISeq tools.

RESULTS

We found 155 DEGs using EdgeR (p-values < 0.05) and 61 DEGs using NOISeq (prob ≥0.8), 30 of which are shared. The rt-qPCR validation suggested a higher reliability of EdgeR results as compared with NOISeq data, in our biological context. Gene Ontology analysis of DEGs identified using EdgeR revealed that green feed modulates biological processes relevant for the rumen physiology and, then, health and well-being of buffaloes, such as lipid metabolism, response to the oxidative stress, immune response, and muscle structure and function. Accordingly, we found: (i) up-regulation of HSD17B13, LOC102410803 (or PSAT1) and HYKK, and down-regulation of CDO1, SELENBP1 and PEMT, encoding factors involved in energy, lipid and amino acid metabolism; (ii) enhanced expression of SIM2 and TRIM14, whose products are implicated in the immune response and defense against infections, and reduced expression of LOC112585166 (or SAAL1), ROR2, SMOC2, and S100A11, encoding pro-inflammatory factors; (iii) up-regulation of NUDT18, DNAJA4 and HSF4, whose products counteract stressful conditions, and down-regulation of LOC102396388 (or UGT1A9) and LOC102413340 (or MRP4/ABCC4), encoding detoxifying factors; (iv) increased expression of KCNK10, CACNG4, and ATP2B4, encoding proteins modulating Ca²⁺ homeostasis, and reduced expression of the cytoskeleton-related MYH11 and DES.

CONCLUSION

Although statistically unpowered, this study suggests that green feed modulates the expression of genes involved in biological processes relevant for rumen functionality and physiology, and thus, for welfare and quality production in Italian Mediterranean dairy buffaloes. These findings, that need to be further confirmed through the validation of additional DEGs, allow to speculate a role of green feed in the production of nutraceutical molecules, whose levels might be enhanced also in milk.

A transcriptomic insight into the human sperm microbiome through next-generation sequencing

The purpose of this study is to provide novel information through Next Generation Sequencing (NGS) for the characterization of viral and bacterial RNA cargo of human sperm cells from healthy fertile donors. For this, RNA-seq raw data of poly(A) RNA from 12 sperm samples from fertile donors were aligned to microbiome databases using the GAIA software. Species of viruses and bacteria were quantified in Operational Taxonomic Units (OTU) and filtered by minimal expression level (>1% OTU in at least one sample). Mean expression values (and their standard deviation) of each species were estimated. A Hierarchical Cluster Analysis (HCA) and a Principal Component Analysis (PCA) were performed to detect common microbiome patterns among samples. Sixteen microbiome species, families, domains, and orders surpassed the established expression threshold. Of the 16 categories, nine corresponded to viruses (23.07% OTU) and seven to bacteria (2.77% OTU), among which the Herperviriales order and Escherichia coli were the most abundant, respectively. HCA and PCA displayed four clusters of samples with a differentiated microbiome fingerprint. This work represents a pilot study into the viruses and bacteria that make up the human sperm microbiome. Despite the high variability observed, some patterns of similarity among individuals were identified. Further NGS studies under standardized methodological procedures are necessary to achieve a deep knowledge of the semen microbiome and its implications in male fertility.

Phylogenetic and Structural Analyses of VPS13 Proteins in Archaeplastida Reveal Their Complex Evolutionary History in Viridiplantae

VPS13 is a lipid transfer protein family conserved among Eukaryotes and playing roles in fundamental processes involving vesicular transport and membrane expansion including autophagy and organelle biogenesis. VPS13 folds into a long hydrophobic tunnel, allowing lipid transport, decorated by distinct domains involved in protein localization and regulation. Whereas VPS13 organization and function have been extensively studied in yeast and mammals, information in organisms originating from primary endosymbiosis is scarce. In the higher plant Arabidopsis thaliana, four paralogs, AtVPS13S, X, M1, and M2, were identified, AtVPS13S playing a role in the regulation of root growth, cell patterning, and reproduction. In this work, we performed phylogenetic, as well as domain and structural modeling of VPS13 proteins in Archaeplastida in order to understand their general organization and evolutionary history. We confirmed the presence of human VPS13B orthologues in some phyla and described two new VPS13 families presenting a particular domain arrangement: VPS13R in Rhodophytes and VPS13Y in Chlorophytes and Streptophytes. By focusing on Viridiplantae, we were able to draw the evolutionary history of these proteins made by multiple gene gains and duplications as well as domain rearrangements. We showed that some Chlorophytes have only three (AtVPS13M, S, Y) whereas some Charophytes have up to six VPS13 paralogs (AtVPS13M1, M2, S, Y, X, B). We also highlighted specific structural features of VPS13M and X paralogs. This study reveals the complex evolution of VPS13 family and opens important perspectives for their functional characterization in photosynthetic organisms.

Exploitation of epigenetic variation of crop wild relatives for crop improvement and agrobiodiversity preservation

Crop wild relatives (CWRs) are recognized as the best potential source of traits for crop improvement. However, successful crop improvement using CWR relies on identifying variation in genes controlling desired traits in plant germplasms and subsequently incorporating them into cultivars. Epigenetic diversity may provide an additional layer of variation within CWR and can contribute novel epialleles for key traits for crop improvement. There is emerging evidence that epigenetic variants of functional and/or agronomic importance exist in CWR gene pools. This provides a rationale for the conservation of epigenotypes of interest, thus contributing to agrobiodiversity preservation through conservation and (epi)genetic monitoring. Concepts and techniques of classical and modern breeding should consider integrating recent progress in epigenetics, initially by identifying their association with phenotypic variations and then by assessing their heritability and stability in subsequent generations. New tools available for epigenomic analysis offer the opportunity to capture epigenetic variation and integrate it into advanced (epi)breeding programmes. Advances in -omics have provided new insights into the sources and inheritance of epigenetic variation and enabled the efficient introduction of epi-traits from CWR into crops using epigenetic molecular markers, such as epiQTLs.

P421 Whole-transcriptome analysis of Sporothrix brasiliensis grown in mold- and yeast-inducing conditions

Poster session 3, September 23, 2022, 12:30 PM - 1:30 PM

 

Objectives

Sporothrix brasiliensis is an emerging Sporothrix species limited to Brazil capable of causing sporotrichosis in humans and animals, especially in cats. Like other pathogenic Sporothrix species, S. brasiliensis exhibits a temperature-dependent dimorphic switch and is therefore, able to undergo a reversible morphological transition (mold and yeast), in response to environmental thermal stimuli.

While dimorphism appears to be essential for virulence in Sporothrix spp, the molecular mechanisms involved in this phenomenon have not yet been fully elucidated.

In this study, we used the strand-specific RNA-Seq technique and bioinformatics analysis to investigate the transcriptomic signatures associated with mold and yeast phases of S. brasiliensis. Furthermore, we generated an accurate version of the S. brasiliensis genome annotation in order to perform high-quality gene expression analysis and other functional or structural genomic studies.

Methods

The whole transcriptome of S. brasiliensis ATCC-MYA-4823, grown in both yeast-inducing (YPD medium at 37°C) and mold-inducing (YPD medium at 25°C) conditions, was sequenced in this study. High-quality RNA was used to prepare Illumina TruSeq Stranded mRNA-paired-end sequencing libraries (2 × 150 bp) that were sequenced using the HiSeq-2500 platform. A total of three biological replicates were sequenced for each condition.

Before transcriptome assembly, adapters and low-quality reads (Phred-score <25) were removed. The StringTie software was used to assemble the transcriptomes imported into the Apollo webtool to manually curate the genome annotation. Transcripts were investigated using TransDecoder and CPC2 programs to determine whether a gene was potentially protein-coding or non-coding. Finally, differential gene expression analysis between yeast and mold forms of S. brasiliensis was conducted using the edgeR package.

Results

Illumina sequencing resulted in a total of ∼217 million raw reads. After quality filtering and trimming, ∼99.5% of reads were used for downstream bioinformatics analysis. The updated S. brasiliensis genome annotation consisted of a total of 14 664 genes of which 10 243 protein-coding genes, 4259 lncRNAs, 140 tRNAs, and 22 rRNAs.

Gene expression analysis revealed a total of 13 838 and 13 938 transcripts expressed in mold- and yeast-form, respectively. Of these, 192 and 292 were expressed exclusively in the mold and yeast-phase, respectively. Moreover, a total of 6802 genes (FDR <0.05) were differentially expressed between the two examined conditions. In particular, 3420 of these genes were up-regulated in the yeast-form (2450 coding, 970 non-coding), and 3382 genes in the filamentous form (2507 coding, 875 non-coding). The raw reads have been deposited into the SRA database and are available under BioProjectID: PRJNA646214.

Conclusions

The characterization of the whole-transcriptome of S. brasiliensis mycelial and yeast-like forms represents an essential starting point for investigating the molecular pathways and regulatory frameworks associated with these two morphological stages. Our results provide new insight into global gene expression profiles of S. brasiliensis, emphasizing the role of non-coding RNAs in its complex transcriptional network.

All transcriptomic data have also been integrated into the ‘Sporothrix Genome DataBase’ (www. sporothrixgenomedatabase.unime.it) in order to expand the current knowledge of Sporothrix genomics and to allow a more in-depth structural exploration of S. brasiliensis gene models, including gene expression patterns related to its saprophytic and pathogenic lifestyle.

LINE-1 RNA causes heterochromatin erosion and is a target for amelioration of senescent phenotypes in progeroid syndromes

Constitutive heterochromatin is responsible for genome repression of DNA enriched in repetitive sequences, telomeres, and centromeres. During physiological and pathological premature aging, heterochromatin homeostasis is profoundly compromised. Here, we showed that LINE-1 (Long Interspersed Nuclear Element-1; L1) RNA accumulation was an early event in both typical and atypical human progeroid syndromes. L1 RNA negatively regulated the enzymatic activity of the histone-lysine N-methyltransferase SUV39H1 (suppression of variegation 3-9 homolog 1), resulting in heterochromatin loss and onset of senescent phenotypes in vitro. Depletion of L1 RNA in dermal fibroblast cells from patients with different progeroid syndromes using specific antisense oligonucleotides (ASOs) restored heterochromatin histone 3 lysine 9 and histone 3 lysine 27 trimethylation marks, reversed DNA methylation age, and counteracted the expression of senescence-associated secretory phenotype genes such as p16, p21, activating transcription factor 3 (ATF3), matrix metallopeptidase 13 (MMP13), interleukin 1a (IL1a), BTG anti-proliferation factor 2 (BTG2), and growth arrest and DNA damage inducible beta (GADD45b). Moreover, systemic delivery of ASOs rescued the histophysiology of tissues and increased the life span of a Hutchinson-Gilford progeria syndrome mouse model. Transcriptional profiling of human and mouse samples after L1 RNA depletion demonstrated that pathways associated with nuclear chromatin organization, cell proliferation, and transcription regulation were enriched. Similarly, pathways associated with aging, inflammatory response, innate immune response, and DNA damage were down-regulated. Our results highlight the role of L1 RNA in heterochromatin homeostasis in progeroid syndromes and identify a possible therapeutic approach to treat premature aging and related syndromes.

Pseudoautosomal Region 1 Overdosage Affects the Global Transcriptome in iPSCs From Patients With Klinefelter Syndrome and High-Grade X Chromosome Aneuploidies

Klinefelter syndrome (KS) is the most prevalent aneuploidy in males and is characterized by a 47,XXY karyotype. Less frequently, higher grade sex chromosome aneuploidies (HGAs) can also occur. Here, using a paradigmatic cohort of KS and HGA induced pluripotent stem cells (iPSCs) carrying 49,XXXXY, 48,XXXY, and 47,XXY karyotypes, we identified the genes within the pseudoautosomal region 1 (PAR1) as the most susceptible to dosage-dependent transcriptional dysregulation and therefore potentially responsible for the progressively worsening phenotype in higher grade X aneuploidies. By contrast, the biallelically expressed non-PAR escape genes displayed high interclonal and interpatient variability in iPSCs and differentiated derivatives, suggesting that these genes could be associated with variable KS traits. By interrogating KS and HGA iPSCs at the single-cell resolution we showed that PAR1 and non-PAR escape genes are not only resilient to the X-inactive specific transcript (XIST)-mediated inactivation but also that their transcriptional regulation is disjointed from the absolute XIST expression level. Finally, we explored the transcriptional effects of X chromosome overdosage on autosomes and identified the nuclear respiratory factor 1 (NRF1) as a key regulator of the zinc finger protein X-linked (ZFX). Our study provides the first evidence of an X-dosage-sensitive autosomal transcription factor regulating an X-linked gene in low- and high-grade X aneuploidies.

Multi-omics data integration provides insights into the post-harvest biology of a long shelf-life tomato landrace

In this study we investigated the transcriptome and epigenome dynamics of the tomato fruit during post-harvest in a landrace belonging to a group of tomatoes (Solanum lycopersicum L.) collectively known as "Piennolo del Vesuvio", all characterized by a long shelf-life. Expression of protein-coding genes and microRNAs as well as DNA methylation patterns and histone modifications were analysed in distinct post-harvest phases. Multi-omics data integration contributed to the elucidation of the molecular mechanisms underlying processes leading to long shelf-life. We unveiled global changes in transcriptome and epigenome. DNA methylation increased and the repressive histone mark H3K27me3 was lost as the fruit progressed from red ripe to 150 days post-harvest. Thousands of genes were differentially expressed, about half of which were potentially epi-regulated as they were engaged in at least one epi-mark change in addition to being microRNA targets in ~5% of cases. Down-regulation of the ripening regulator MADS-RIN and of genes involved in ethylene response and cell wall degradation was consistent with the delayed fruit softening. Large-scale epigenome reprogramming that occurred in the fruit during post-harvest likely contributed to delayed fruit senescence.

Complete genome sequence of the biocontrol yeast Papiliotrema terrestris strain LS28

Papiliotrema terrestris strain LS28 is a biocontrol agent selected for its antagonistic activity against several plant pathogens both in the field and postharvest. The availability of a genome sequencing sets the foundation for the identification of the genetic mechanisms of its antagonistic activity. The genome size is 21.29 Mbp with a G+C content of 58.65%, and genome annotation predicts 8,626 protein-encoding genes. Phylogenetic analysis based on whole-genome data confirms that P. terrestris is a Tremellomycetes more closely related to Papiliotrema flavescens than Papiliotrema laurentii.

PRGdb 4.0: an updated database dedicated to genes involved in plant disease resistance process

The Plant Resistance Genes database (PRGdb; http://prgdb.org/prgdb4/) has been greatly expanded, keeping pace with the increasing amount of available knowledge and data (sequenced proteomes, cloned genes, public analysis data, etc.). The easy-to-use style of the database website has been maintained, while an updated prediction tool, more data and a new section have been added. This new section will contain plant resistance transcriptomic experiments, providing additional easy-to-access experimental information. DRAGO3, the tool for automatic annotation and prediction of plant resistance genes behind PRGdb, has been improved in both accuracy and sensitivity, leading to more reliable predictions. PRGdb offers 199 reference resistance genes and 586.652 putative resistance genes from 182 sequenced proteomes. Compared to the previous release, PRGdb 4.0 has increased the number of reference resistance genes from 153 to 199, the number of putative resistance genes from 177K from 76 proteomes to 586K from 182 sequenced proteomes. A new section has been created that collects plant-pathogen transcriptomic data for five species of agricultural interest. Thereby, with these improvements and data expansions, PRGdb 4.0 aims to serve as a reference to the plant scientific community and breeders worldwide, helping to further study plant resistance mechanisms that contribute to fighting pathogens.

isoCNV: in silico optimization of copy number variant detection from targeted or exome sequencing data

Background

Accurate copy number variant (CNV) detection is especially challenging for both targeted sequencing (TS) and whole‐exome sequencing (WES) data. To maximize the performance, the parameters of the CNV calling algorithms should be optimized for each specific dataset. This requires obtaining validated CNV information using either multiplex ligation-dependent probe amplification (MLPA) or array comparative genomic hybridization (aCGH). They are gold standard but time-consuming and costly approaches.

Results

We present isoCNV which optimizes the parameters of DECoN algorithm using only NGS data. The parameter optimization process is performed using an in silico CNV validated dataset obtained from the overlapping calls of three algorithms: CNVkit, panelcn.MOPS and DECoN. We evaluated the performance of our tool and showed that increases the sensitivity in both TS and WES real datasets.

Conclusions

isoCNV provides an easy-to-use pipeline to optimize DECoN that allows the detection of analysis-ready CNV from a set of DNA alignments obtained under the same conditions. It increases the sensitivity of DECoN without the need for orthogonal methods. isoCNV is available at https://gitlab.com/sequentiateampublic/isocnv.

An Oil Hyper-Accumulator Mutant Highlights Peroxisomal ATP Import as a Regulatory Step for Fatty Acid Metabolism in Aurantiochytrium limacinum

Thraustochytrids are marine protists that naturally accumulate triacylglycerol with long chains of polyunsaturated fatty acids, such as ω3-docosahexaenoic acid (DHA). They represent a sustainable response to the increasing demand for these "essential" fatty acids (FAs). Following an attempt to transform a strain of Aurantiochytrium limacinum, we serendipitously isolated a clone that did not incorporate any recombinant DNA but contained two to three times more DHA than the original strain. Metabolic analyses indicated a deficit in FA catabolism. However, whole transcriptome analysis did not show down-regulation of genes involved in FA catabolism. Genome sequencing revealed extensive DNA deletion in one allele encoding a putative peroxisomal adenylate transporter. Phylogenetic analyses and yeast complementation experiments confirmed the gene as a peroxisomal adenylate nucleotide transporter (AlANT1), homologous to yeast ScANT1 and plant peroxisomal adenylate nucleotide carrier AtPNC genes. In yeast and plants, a deletion of the peroxisomal adenylate transporter inhibits FA breakdown and induces FA accumulation, a phenotype similar to that described here. In response to this metabolic event, several compensatory mechanisms were observed. In particular, genes involved in FA biosynthesis were upregulated, also contributing to the high FA accumulation. These results support AlANT1 as a promising target for enhancing DHA production in Thraustochytrids.

Early Alpine occupation backdates westward human migration in Late Glacial Europe

Before the end of the Last Glacial Maximum (LGM, ∼16.5 ka ago)¹ set in motion major shifts in human culture and population structure,² a consistent change in lithic technology, material culture, settlement pattern, and adaptive strategies is recorded in Southern Europe at ∼18-17 ka ago. In this time frame, the landscape of Northeastern Italy changed considerably, and the retreat of glaciers allowed hunter-gatherers to gradually recolonize the Alps.^3-6 Change within this renewed cultural frame (i.e., during the Late Epigravettian phase) is currently associated with migrations favored by warmer climate linked to the Bølling-Allerød onset (14.7 ka ago),^7-11 which replaced earlier genetic lineages with ancestry found in an individual who lived ∼14 ka ago at Riparo Villabruna, Italy, and shared among different contexts (Villabruna Cluster).⁹ Nevertheless, these dynamics and their chronology are still far from being disentangled due to fragmentary evidence for long-distance interactions across Europe.¹² Here, we generate new genomic data from a human mandible uncovered at Riparo Tagliente (Veneto, Italy), which we directly dated to 16,980-16,510 cal BP (2σ). This individual, affected by focal osseous dysplasia, is genetically affine to the Villabruna Cluster. Our results therefore backdate by at least 3 ka the diffusion in Southern Europe of a genetic component linked to Balkan/Anatolian refugia, previously believed to have spread during the later Bølling/Allerød event. In light of the new genetic evidence, this population replacement chronologically coincides with the very emergence of major cultural transitions in Southern and Western Europe.

Meiocyte Isolation by INTACT and Meiotic Transcriptome Analysis in Arabidopsis

Isolation of nuclei tagged in specific cell types (INTACT) is a method developed to isolate cell-type-specific nuclei that are tagged through in vivo biotin labeling of a nuclear targeting fusion (NTF) protein. In our work, INTACT was used to capture nuclei of meiocytes and to generate a meiotic transcriptome in Arabidopsis. Using the promoter of AtDMC1 recombinase to label meiotic nuclei, we generated transgenic plants carrying AtDMC1:NTF along with biotin ligase enzyme (BirA) under the constitutive ACTIN2 (ACT2) promoter. AtDMC1-driven expression of biotin-labeled NTF allowed us to collect nuclei of meiocytes by streptavidin-coated magnetic beads. The nuclear meiotic transcriptome was obtained by RNA-seq using low-quantity input RNA. Transcripts grouped into different categories according to their expression levels were investigated by gene ontology enrichment analysis (GOEA). The most enriched GO term "DNA demethylation" in mid/high-expression classes suggests that this biological process is particularly relevant to meiosis onset. The majority of genes with established roles in meiosis were distributed in the classes of mid/high and high expression. Meiotic transcriptome was compared with public available transcriptomes from other tissues in Arabidopsis. Bioinformatics analysis by expression network identified a core of more than 1,500 genes related to meiosis landmarks.

Transcriptome-wide expression profiling of Sporothrix schenckii yeast and mycelial forms and the establishment of the Sporothrix Genome DataBase

Sporothrix schenckii is a dimorphic fungus existing as mould in the environment and as yeast in the host. The morphological shift between mycelial/yeast phases is crucial for its virulence, but the transcriptional networks implicated in dimorphic transition are still not fully understood. Here, we report the global transcriptomic differences occurring between mould and yeast phases of S. schenckii, including changes in gene expression profiles associated with these distinct cellular phenotypes. Moreover, we also propose a new genome annotation, which reveals a more complex transcriptional architecture than previously assumed. Using RNA-seq, we identified a total of 17 307 genes, of which 11 217 were classified as protein-encoding genes, whereas 6090 were designated as non-coding RNAs (ncRNAs). Approximately ~71 % of all annotated genes were found to overlap and the different-strand overlapping type was the most common. Gene expression analysis revealed that 8795 genes were differentially regulated among yeast and mould forms. Differential gene expression was also observed for antisense ncRNAs overlapping neighbouring protein-encoding genes. The release of transcriptome-wide data and the establishment of the Sporothrix Genome DataBase (http://sporothrixgenomedatabase.unime.it) represent an important milestone for Sporothrix research, because they provide a strong basis for future studies on the molecular pathways involved in numerous biological processes.

Fruit Development in Ficus carica L.: Morphological and Genetic Approaches to Fig Buds for an Evolution From Monoecy Toward Dioecy

The mechanism behind the bud evolution towards breba or main crop in Ficus carica L. is uncertain. Anatomical and genetic studies may put a light on the possible similarities/differences between the two types of fruits. For this reason, we collected complimentary data from anatomical, X-ray imaging, and genetic techniques. The RNA seq together with structural genome annotation allowed the prediction of 34,629 known genes and 938 novel protein-coding genes. Transcriptome analysis of genes during bud differentiation revealed differentially expressed genes in two fig varieties (Dottato and Petrelli) and in breba and main crop. We chose Dottato and Petrelli because the first variety does not require pollination to set main crop and the latter does; moreover, Petrelli yields many brebas whereas Dottato few. Of the 1,615 and 1,904 loci expressed in Dottato and Petrelli, specifically in breba or main crop, respectively, only 256 genes appeared to be transcripts in both varieties. The buds of the two fig varieties were observed under optical microscope and using 3D X-ray tomography, highlighting differences mainly related to the stage of development. The X-ray images of buds showed a great structural similarity between breba and main crop during the initial stages of development. Analysis at the microscope indicated that inflorescence differentiation of breba was split in two seasons whereas that of main crop started at the end of winter of season 2 and was completed within 2 to 3 months. The higher expression of floral homeotic protein AGAMOUS in breba with respect to main crop, since this protein is required for normal development of stamens and carpels in the flower, may indicate an original role of these fruits for staminate flowers production for pollination of the main crop, as profichi in the caprifig. Several genes related to auxin (auxin efflux carrier, auxin response factor, auxin binding protein, auxin responsive protein) and to GA synthesis (GA20ox) were highly expressed in brebas with respect to main crop for the development of this parthenocarpic fruit.

Illumina and PacBio DNA sequencing data, de novo assembly and annotation of the genome of Aurantiochytrium limacinum strain CCAP_4062/1

The complete genome of the thraustochytrid Aurantiochytrium limacinum strain CCAP_4062/1 was sequenced using both Illumina Novaseq 6000 and third generation sequencing technology PacBio RSII in order to obtain trustworthy assembly and annotation. The reads from both platforms were combined at multiple levels in order to obtain a reliable assembly, then compared to the A. limacinum ATCCⓇ MYA1381™ reference genome. The final assembly was annotated with the help of strain CCAP_4062/1 RNAseq data. A. limacinum strain CCAP_4062/1 is an industrial strain used for the production of very long chain polyunsaturated fatty acids, like the docosahexaenoic acid that is an essential fatty acid synthesised only at very low pace in humans and vertebrates . Thraustochytrids in general and Aurantiochytrium more specifically, are used for carotenoid and squalene production as well. Beside their biotechnological interest, thraustochytrids play a crucial role in both inshore and oceanic basins ecosystems. Genome sequences will foster biotechnological as well as ecological studies.

Transcriptome analysis of non-ochratoxigenic Aspergillus carbonarius strains and interactions between some black aspergilli species

Aspergillus carbonarius consistently produces large amounts of ochratoxin A (OTA), a mycotoxin with nephrotoxic effects on animals and humans. In the present study, we analyzed the transcriptional changes associated to OTA production in three atypical non-ochratoxigenic strains of A. carbonarius. In addition, in vitro interactions between ochratoxigenic strains of A. carbonarius and A. niger and non-ochratoxigenic strains of A. carbonarius and A. tubingensis were studied in order to evaluate their potential for controlling OTA production. RNA-seq analysis revealed that there are 696 differentially expressed genes identified in the three non-OTA-producing strains, including 280 up-regulated and 333 down-regulated genes. A functional and gene ontology enrichment analysis revealed that the processes related to metabolic and oxidation processes, associated with functions such as oxidoreductase and hydrolase activity were down regulated. All the genes related with OTA biosynthesis in A. carbonarius were the most down-regulated genes in non-ochratoxigenic strains. We also showed that these strains possess a deleterious mutation in the AcOTApks gene required for OTA biosynthesis. Moreover, one of these strains gave the best control of OTA production resulting in an OTA reduction of 98-100% in co-inoculation with an ochratoxigenic strain of A. niger and an OTA reduction of 79-89% with an ochratoxigenic strain of A. carbonarius. Results of this study provided novel insights into the knowledge of the OTA biosynthetic pathway in these non-ochratoxigenic wild strains, and showed the biocontrol potential of these strains.

Recombination suppression in heterozygotes for a pericentric inversion induces the interchromosomal effect on crossovers in Arabidopsis

During meiosis, recombination ensures allelic exchanges through crossovers (COs) between the homologous chromosomes. Advances in our understanding of the rules of COs have come from studies of mutations including structural chromosomal rearrangements that, when heterozygous, are known to impair COs in various organisms. In this work, we investigate the effect of a large heterozygous pericentric inversion on male and female recombination in Arabidopsis. The inversion was discovered in the Atmcc1 mutant background and was characterized through genetic and next-generation sequencing analysis. Reciprocal backcross populations, each consisting of over 400 individuals, obtained from the mutant and the wild type, both crossed with Landsberg erecta, were analyzed genome-wide by 143 single-nucleotide polymorphisms. The negative impact of inversion became evident in terms of CO loss in the rearranged chromosome in both male and female meiosis. No single-CO event was detected within the inversion, consistent with a post-meiotic selection operating against unbalanced gametes. Cytological analysis of chiasmata in F₁ plants confirmed that COs were reduced in male meiosis in the chromosome with inversion. Crossover suppression on the rearranged chromosome is associated with a significant increase of COs in the other chromosomes, thereby maintaining unchanged the number of COs per cell. The CO pattern observed in our study is consistent with the interchromosomal (IC) effect as first described in Drosophila. In contrast to male meiosis, in female meiosis no IC effect is visible. This may be related to the greater strength of interference that constrains the CO number in excess of the minimum value imposed by CO assurance in Arabidopsis female meiosis.

Whole genome SNPs discovery in Nero Siciliano pig

Autochthonous pig breeds represent an important genetic reserve to be utilized mainly for the production of typical products. To explore its genetic variability, here we present for the first time whole genome sequencing data and SNPs discovered in a male domestic Nero Siciliano pig compared to the last pig reference genome Sus scrofa11.1.A total of 346.8 million paired reads were generated by sequencing. After quality control, 99.03% of the reads were mapped to the reference genome, and over 11 million variants were detected.Additionally, we evaluated sequence diversity in 21 fitness-related loci selected based on their biological function and/or their proximity to relevant QTLs. We focused on genes that have been related to environmental adaptation and reproductive traits in previous studies regarding local breeds. A total of 6,747 variants were identified resulting in a rate of 1 variant every ~276 bases. Among these variants 1,132 were novel to the dbSNP151 database. This study represents a first step in the genetic characterization of Nero Siciliano pig and also provides a platform for future comparative studies between this and other swine breeds.

Epigenetic remodelling licences adult cholangiocytes for organoid formation and liver regeneration

Following severe or chronic liver injury, adult ductal cells (cholangiocytes) contribute to regeneration by restoring both hepatocytes and cholangiocytes. We recently showed that ductal cells clonally expand as self-renewing liver organoids that retain their differentiation capacity into both hepatocytes and ductal cells. However, the molecular mechanisms by which adult ductal-committed cells acquire cellular plasticity, initiate organoids and regenerate the damaged tissue remain largely unknown. Here, we describe that ductal cells undergo a transient, genome-wide, remodelling of their transcriptome and epigenome during organoid initiation and in vivo following tissue damage. TET1-mediated hydroxymethylation licences differentiated ductal cells to initiate organoids and activate the regenerative programme through the transcriptional regulation of stem-cell genes and regenerative pathways including the YAP-Hippo signalling. Our results argue in favour of the remodelling of genomic methylome/hydroxymethylome landscapes as a general mechanism by which differentiated cells exit a committed state in response to tissue damage.

The lipid metabolism in thraustochytrids

Thraustochytrids are unicellular heterotrophic marine protists of the Stramenopile group, often considered as non-photosynthetic microalgae. They have been isolated from a wide range of habitats including deep sea, but are mostly present in waters rich in sediments and organic materials. They are abundant in mangrove forests where they are major colonizers, feeding on decaying leaves and initiating the mangrove food web. Discovered 80 years ago, they have recently attracted considerable attention due to their biotechnological potential. This interest arises from their fast growth, their specific lipid metabolism and the improvement of the genetic tools and transformation techniques. These organisms are particularly rich in ω3-docosahexaenoic acid (DHA), an 'essential' fatty acid poorly encountered in land plants and animals but required for human health. To produce their DHA, thraustochytrids use a sophisticated system different from the classical fatty acid synthase system. They are also a potential source of squalene and carotenoids. Here we review our current knowledge about the life cycle, ecophysiology, and metabolism of these organisms, with a particular focus on lipid dynamics. We describe the different pathways involved in lipid and fatty acid syntheses, emphasizing their specificity, and we report on the recent efforts aimed to engineer their lipid metabolism.

Whole-Genome Sequencing of an Uncommon Cryptococcus neoformans MLST43 Genotype Isolated in Nigeria

Cryptococcosis is a human infection caused mainly by two species of the Cryptococcus genus, Cryptococcus neoformans and Cryptococcus gattii, whose populations contain several phylogenetically related haploid (VN/VG-types) and hybrid genotypes. Here, we report the whole-genome sequencing of a Nigerian C. neoformans VNII, Mat-α, strain with a rare multilocus-sequence-type (MLST) genotype (ST43).

LOXL2-mediated H3K4 oxidation reduces chromatin accessibility in triple-negative breast cancer cells

Oxidation of H3 at lysine 4 (H3K4ox) by lysyl oxidase-like 2 (LOXL2) generates an H3 modification with an unknown physiological function. We find that LOXL2 and H3K4ox are higher in triple-negative breast cancer (TNBC) cell lines and patient-derived xenografts (PDXs) than those from other breast cancer subtypes. ChIP-seq revealed that H3K4ox is located primarily in heterochromatin, where it is involved in chromatin compaction. Knocking down LOXL2 reduces H3K4ox levels and causes chromatin decompaction, resulting in a sustained activation of the DNA damage response (DDR) and increased susceptibility to anticancer agents. This critical role that LOXL2 and oxidized H3 play in chromatin compaction and DDR suggests that functionally targeting LOXL2 could be a way to sensitize TNBC cells to conventional therapy.

PRGdb 3.0: a comprehensive platform for prediction and analysis of plant disease resistance genes

The Plant Resistance Genes database (PRGdb; http://prgdb.org) has been redesigned with a new user interface, new sections, new tools and new data for genetic improvement, allowing easy access not only to the plant science research community but also to breeders who want to improve plant disease resistance. The home page offers an overview of easy-to-read search boxes that streamline data queries and directly show plant species for which data from candidate or cloned genes have been collected. Bulk data files and curated resistance gene annotations are made available for each plant species hosted. The new Gene Model view offers detailed information on each cloned resistance gene structure to highlight shared attributes with other genes. PRGdb 3.0 offers 153 reference resistance genes and 177 072 annotated candidate Pathogen Receptor Genes (PRGs). Compared to the previous release, the number of putative genes has been increased from 106 to 177 K from 76 sequenced Viridiplantae and algae genomes. The DRAGO 2 tool, which automatically annotates and predicts (PRGs) from DNA and amino acid with high accuracy and sensitivity, has been added. BLAST search has been implemented to offer users the opportunity to annotate and compare their own sequences. The improved section on plant diseases displays useful information linked to genes and genomes to connect complementary data and better address specific needs. Through, a revised and enlarged collection of data, the development of new tools and a renewed portal, PRGdb 3.0 engages the plant science community in developing a consensus plan to improve knowledge and strategies to fight diseases that afflict main crops and other plants.

Three-Dimensional Genomic Structure and Cohesin Occupancy Correlate with Transcriptional Activity during Spermatogenesis

Mammalian gametogenesis involves dramatic and tightly regulated chromatin remodeling, whose regulatory pathways remain largely unexplored. Here, we generate a comprehensive high-resolution structural and functional atlas of mouse spermatogenesis by combining in situ chromosome conformation capture sequencing (Hi-C), RNA sequencing (RNA-seq), and chromatin immunoprecipitation sequencing (ChIP-seq) of CCCTC-binding factor (CTCF) and meiotic cohesins, coupled with confocal and super-resolution microscopy. Spermatogonia presents well-defined compartment patterns and topological domains. However, chromosome occupancy and compartmentalization are highly re-arranged during prophase I, with cohesins bound to active promoters in DNA loops out of the chromosomal axes. Compartment patterns re-emerge in round spermatids, where cohesin occupancy correlates with transcriptional activity of key developmental genes. The compact sperm genome contains compartments with actively transcribed genes but no fine-scale topological domains, concomitant with the presence of protamines. Overall, we demonstrate how genome-wide cohesin occupancy and transcriptional activity is associated with three-dimensional (3D) remodeling during spermatogenesis, ultimately reprogramming the genome for the next generation.

A de novo transcriptome assembly approach elucidates the dynamics of ovarian maturation in the swordfish (Xiphias gladius)

The Mediterranean swordfish (Xiphias gladius) has been recently classified as overfished and in 2016, the International Commission for the Conservation of the Atlantic Tunas (ICCAT) established a multi-annual management plan to recover this stock. To successfully achieve this goal, knowledge about swordfish biology is needed. To date, few studies on swordfish have been performed and none of them has provided useful insights into the reproductive biology at molecular level. Here we set to characterise the molecular dynamics underlying ovarian maturation by employing a de novo transcriptome assembly approach. Differential gene expression analysis in mature and immature ovaries identified a number of differentially expressed genes associated with biological processes driving ovarian maturation. Focusing on ovarian steroidogenesis and vitellogenin uptake, we depict the molecular dynamics characterizing these processes while a phylogenetic analysis let us identify a candidate vitellogenin receptor. This is the first swordfish transcriptome assembly and these findings provide in-depth understanding of molecular processes describing ovarian maturation. Moreover, the establishment of a publicly available database containing information on the swordfish transcriptome aims to boost research on this species with the long-term of developing more comprehensive and successful stock management plans.

Assessment of genomic changes in a CRISPR/Cas9 Phaeodactylum tricornutum mutant through whole genome resequencing

The clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9 system, co-opted from a bacterial defense natural mechanism, is the cutting edge technology to carry out genome editing in a revolutionary fashion. It has been shown to work in many different model organisms, from human to microbes, including two diatom species, Phaeodactylum tricornutum and Thalassiosira pseudonana. Transforming P. tricornutum by bacterial conjugation, we have performed CRISPR/Cas9-based mutagenesis delivering the nuclease as an episome; this allowed for avoiding unwanted perturbations due to random integration in the genome and for excluding the Cas9 activity when it was no longer required, reducing the probability of obtaining off-target mutations, a major drawback of the technology. Since there are no reports on off-target occurrence at the genome level in microalgae, we performed whole-genome Illumina sequencing and found a number of different unspecific changes in both the wild type and mutant strains, while we did not observe any preferential mutation in the genomic regions in which off-targets were predicted. Our results confirm that the CRISPR/Cas9 technology can be efficiently applied to diatoms, showing that the choice of the conjugation method is advantageous for minimizing unwanted changes in the genome of P. tricornutum.

Whole-genome re-sequencing of two Italian tomato landraces reveals sequence variations in genes associated with stress tolerance, fruit quality and long shelf-life traits

Tomato is a high value crop and the primary model for fleshy fruit development and ripening. Breeding priorities include increased fruit quality, shelf life and tolerance to stresses. To contribute towards this goal, we re-sequenced the genomes of Corbarino (COR) and Lucariello (LUC) landraces, which both possess the traits of plant adaptation to water deficit, prolonged fruit shelf-life and good fruit quality. Through the newly developed pipeline Reconstructor, we generated the genome sequences of COR and LUC using datasets of 65.8 M and 56.4 M of 30-150 bp paired-end reads, respectively. New contigs including reads that could not be mapped to the tomato reference genome were assembled, and a total of 43, 054 and 44, 579 gene loci were annotated in COR and LUC. Both genomes showed novel regions with similarity to Solanum pimpinellifolium and Solanum pennellii. In addition to small deletions and insertions, 2, 000 and 1, 700 single nucleotide polymorphisms (SNPs) could exert potentially disruptive effects on 1, 371 and 1, 201 genes in COR and LUC, respectively. A detailed survey of the SNPs occurring in fruit quality, shelf life and stress tolerance related-genes identified several candidates of potential relevance. Variations in ethylene response components may concur in determining peculiar phenotypes of COR and LUC.

Maize RNA PolIV affects the expression of genes with nearby TE insertions and has a genome-wide repressive impact on transcription

BACKGROUND

RNA-directed DNA methylation (RdDM) is a plant-specific epigenetic process that relies on the RNA polymerase IV (Pol IV) for the production of 24 nucleotide small interfering RNAs (siRNA) that guide the cytosine methylation and silencing of genes and transposons. Zea mays RPD1/RMR6 gene encodes the largest subunit of Pol IV and is required for normal plant development, paramutation, transcriptional repression of certain transposable elements (TEs) and transcriptional regulation of specific alleles.

RESULTS

In this study we applied a total RNA-Seq approach to compare the B73 and rpd1/rmr6 leaf transcriptomes. Although previous studies indicated that loss of siRNAs production in RdDM mutants provokes a strong loss of CHH DNA methylation but not massive gene or TEs transcriptional activation in both Arabidopsis and maize, our total RNA-Seq analysis of rpd1/rmr6 transcriptome reveals that loss of Pol IV activity causes a global increase in the transcribed fraction of the maize genome. Our results point to the genes with nearby TE insertions as being the most strongly affected by Pol IV-mediated gene silencing. TEs modulation of nearby gene expression is linked to alternative methylation profiles on gene flanking regions, and these profiles are strictly dependent on specific characteristics of the TE member inserted. Although Pol IV is essential for the biogenesis of siRNAs, the genes with associated siRNA loci are less affected by the pol IV mutation.

CONCLUSIONS

This deep and integrated analysis of gene expression, TEs distribution, smallRNA targeting and DNA methylation levels, reveals that loss of Pol IV activity globally affects genome regulation, pointing at TEs as modulator of nearby gene expression and indicating the existence of multiple level epigenetic silencing mechanisms. Our results also suggest a predominant role of the Pol IV-mediated RdDM pathway in genome dominance regulation, and subgenome stability and evolution in maize.

What the human sperm methylome tells us

Aim: To characterize the sperm methylome in semen samples from 19 donors with proven fertility. Materials & methods: Bisulfite-converted sperm DNA was hybridized on the HumanMethylation450 Infinium BeadChip platform. CpG fluorescence intensities were extracted and converted to β-values. Results: The sperm methylome is highly homogeneous and hypomethylated. Genes with hypomethylated promoters are ontologically associated to biological functions related to spermatogenesis and embryogenesis. Sex chromosomes are the most hypomethylated chromosomes, supporting data that indicated their essential role in spermatogenesis. A total of 94 genes are resistant to demethylation, being strong candidates for transgenerational inheritance. Conclusion: Spermatozoa carry a homogeneous methylation profile that is a footprint of past events (spermatogenesis), is designed to facilitate future events (embryogenesis) and has a possible influence in the adult life (transgenerational effects).

The hypothesis that Helicobacter pylori predisposes to Alzheimer's disease is biologically plausible

There is epidemiological evidence that H. pylori might predispose to Alzheimer's disease. To understand the cellular processes potentially linking such unrelated events, we incubated the human gastric cells MNK-28 with the H. pylori peptide Hp(2-20). We then monitored the activated genes by global gene expression. The peptide modulated 77 genes, of which 65 are listed in the AlzBase database and include the hallmarks of Alzheimer's disease: APP, APOE, PSEN1, and PSEN2. A large fraction of modulated genes (30 out of 77) belong to the inflammation pathway. Remarkably, the pathways dis-regulated in Alzheimer's and Leasch-Nyhan diseases result dis-regulated also in this study. The unsuspected links between such different diseases - though still awaiting formal validation - suggest new directions for the study of neurological diseases.

The hypothesis that Helicobacter pylori predisposes to Alzheimer’s disease is biologically plausible

There is epidemiological evidence that H. pylori might predispose to Alzheimer’s disease. To understand the cellular processes potentially linking such unrelated events, we incubated the human gastric cells MNK-28 with the H. pylori peptide Hp(2-20). We then monitored the activated genes by global gene expression. The peptide modulated 77 genes, of which 65 are listed in the AlzBase database and include the hallmarks of Alzheimer’s disease: APP, APOE, PSEN1, and PSEN2. A large fraction of modulated genes (30 out of 77) belong to the inflammation pathway. Remarkably, the pathways dis-regulated in Alzheimer’s and Leasch-Nyhan diseases result dis-regulated also in this study. The unsuspected links between such different diseases – though still awaiting formal validation – suggest new directions for the study of neurological diseases.

Exploiting the great potential of Sequence Capture data by a new tool, SUPER-CAP

The recent development of Sequence Capture methodology represents a powerful strategy for enhancing data generation to assess genetic variation of targeted genomic regions. Here, we present SUPER-CAP, a bioinformatics web tool aimed at handling Sequence Capture data, fine calculating the allele frequency of variations and building genotype-specific sequence of captured genes. The dataset used to develop this in silico strategy consists of 378 loci and related regulative regions in a collection of 44 tomato landraces. About 14,000 high-quality variants were identified. The high depth (>40×) of coverage and adopting the correct filtering criteria allowed identification of about 4,000 rare variants and 10 genes with a different copy number variation. We also show that the tool is capable to reconstruct genotype-specific sequences for each genotype by using the detected variants. This allows evaluating the combined effect of multiple variants in the same protein. The architecture and functionality of SUPER-CAP makes the software appropriate for a broad set of analyses including SNP discovery and mining. Its functionality, together with the capability to process large data sets and efficient detection of sequence variation, makes SUPER-CAP a valuable bioinformatics tool for genomics and breeding purposes.

Stress-induced and epigenetic-mediated maize transcriptome regulation study by means of transcriptome reannotation and differential expression analysis

Plant's response and adaptation to abiotic stresses involve sophisticated genetic and epigenetic regulatory systems. To obtain a global view of molecular response to osmotic stresses, including the non-coding portion of genome, we conducted a total leaf transcriptome analysis on maize plants subjected to prolonged drought and salt stresses. Stress application to both B73 wild type and the epiregulator mutant rpd1-1/rmr6 allowed dissection of the epigenetic component of stress response. Coupling total RNA-Seq and transcriptome re-assembly we annotated thousands of new maize transcripts, together with 13,387 lncRNAs that may play critical roles in regulating gene expression. Differential expression analysis revealed hundreds of genes modulated by long-term stress application, including also many lncRNAs and transposons specifically induced by stresses. The amplitude and dynamic of the stress-modulated gene sets are very different between B73 and rpd1-1/rmr6 mutant plants, as result of stress-like effect on genome regulation caused by the mutation itself, which activates many stress-related genes even in control condition. The analyzed extensive set of total RNA-Seq data, together with the improvement of the transcriptome and the identification of the non-coding portion of the transcriptome give a revealing insight into the genetic and epigenetic mechanism responsible for maize molecular response to abiotic stresses.

Flowering and trichome development share hormonal and transcription factor regulation

Gibberellins (GAs) and cytokinins (CKs) are plant hormones that act either synergistically or antagonistically during the regulation of different developmental processes. In Arabidopsis thaliana, GAs and CKs overlap in the positive regulation of processes such as the transition from the vegetative to the reproductive phase and the development of epidermal adaxial trichomes. Despite the fact that both developmental processes originate in the rosette leaves, they occur separately in time and space. Here we review how, as genetic and molecular mechanisms are being unraveled, both processes might be closely related. Additionally, this shared genetic network is not only dependent on GA and CK hormone signaling but is also strictly controlled by specific clades of transcription factor families. Some key flowering genes also control other rosette leaf developmental processes such as adaxial trichome formation. Conversely, most of the trichome activator genes, which belong to the MYB, bHLH and C2H2 families, were found to positively control the floral transition. Furthermore, three MADS floral organ identity genes, which are able to convert leaves into floral structures, are also able to induce trichome proliferation in the flower. These data lead us to propose that the spatio-temporal regulation and integration of diverse signals control different developmental processes, such as floral induction and trichome formation, which are intimately connected through similar genetic pathways.

Effects of Mecp2 loss of function in embryonic cortical neurons: a bioinformatics strategy to sort out non-neuronal cells variability from transcriptome profiling

Background

Mecp2 null mice model Rett syndrome (RTT) a human neurological disorder affecting females after apparent normal pre- and peri-natal developmental periods. Neuroanatomical studies in cerebral cortex of RTT mouse models revealed delayed maturation of neuronal morphology and autonomous as well as non-cell autonomous reduction in dendritic complexity of postnatal cortical neurons. However, both morphometric parameters and high-resolution expression profile of cortical neurons at embryonic developmental stage have not yet been studied. Here we address these topics by using embryonic neuronal primary cultures from Mecp2 loss of function mouse model.

Results

We show that embryonic primary cortical neurons of Mecp2 null mice display reduced neurite complexity possibly reflecting transcriptional changes. We used RNA-sequencing coupled with a bioinformatics comparative approach to identify and remove the contribution of variable and hard to quantify non-neuronal brain cells present in our in vitro cell cultures.

Conclusions

Our results support the need to investigate both Mecp2 morphological as well as molecular effect in neurons since prenatal developmental stage, long time before onset of Rett symptoms.

Electronic supplementary material

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

GREENC: a Wiki-based database of plant lncRNAs

Long non-coding RNAs (lncRNAs) are functional non-translated molecules greater than 200 nt. Their roles are diverse and they are usually involved in transcriptional regulation. LncRNAs still remain largely uninvestigated in plants with few exceptions. Experimentally validated plant lncRNAs have been shown to regulate important agronomic traits such as phosphate starvation response, flowering time and interaction with symbiotic organisms, making them of great interest in plant biology and in breeding. There is still a lack of lncRNAs in most sequenced plant species, and in those where they have been annotated, different methods have been used, so making the lncRNAs less useful in comparisons within and between species. We developed a pipeline to annotate lncRNAs and applied it to 37 plant species and six algae, resulting in the annotation of more than 120 000 lncRNAs. To facilitate the study of lncRNAs for the plant research community, the information gathered is organised in the Green Non-Coding Database (GreeNC, http://greenc.sciencedesigners.com/).

The Solanum commersonii Genome Sequence Provides Insights into Adaptation to Stress Conditions and Genome Evolution of Wild Potato Relatives

Here, we report the draft genome sequence of Solanum commersonii, which consists of ∼830 megabases with an N50 of 44,303 bp anchored to 12 chromosomes, using the potato (Solanum tuberosum) genome sequence as a reference. Compared with potato, S. commersonii shows a striking reduction in heterozygosity (1.5% versus 53 to 59%), and differences in genome sizes were mainly due to variations in intergenic sequence length. Gene annotation by ab initio prediction supported by RNA-seq data produced a catalog of 1703 predicted microRNAs, 18,882 long noncoding RNAs of which 20% are shown to target cold-responsive genes, and 39,290 protein-coding genes with a significant repertoire of nonredundant nucleotide binding site-encoding genes and 126 cold-related genes that are lacking in S. tuberosum. Phylogenetic analyses indicate that domesticated potato and S. commersonii lineages diverged ∼2.3 million years ago. Three duplication periods corresponding to genome enrichment for particular gene families related to response to salt stress, water transport, growth, and defense response were discovered. The draft genome sequence of S. commersonii substantially increases our understanding of the domesticated germplasm, facilitating translation of acquired knowledge into advances in crop stability in light of global climate and environmental changes.