Functional characterization of the Hyles euphorbiae hawkmoth transcriptome reveals strong expression of phorbol ester detoxification and seasonal cold hardiness genes.

BACKGROUND: The European spurge hawkmoth, Hyles euphorbiae (Lepidoptera, Sphingidae), has been intensively studied as a model organism for insect chemical ecology, cold hardiness and evolution of species delineation. To understand species isolation mechanisms at a molecular level, this study aims at determining genetic factors underlying two adaptive ecological trait candidates, phorbol ester (TPA) detoxification and seasonal cold acclimation. METHOD: A draft transcriptome of H. euphorbiae was generated using Illumina sequencing, providing the first genomic resource for the hawkmoth subfamily Macroglossinae. RNA expression levels in tissues of experimental TPA feeding larvae and cooled pupae was compared to levels in control larvae and pupae using 26 bp RNA sequence tag libraries (DeepSuperSAGE). Differential gene expression was assessed by homology searches of the tags in the transcriptome. RESULTS: In total, 389 and 605 differentially expressed transcripts for detoxification and cold hardiness, respectively, could be identified and annotated with proteins. The majority (22 of 28) of differentially expressed detox transcripts of the four 'drug metabolism' enzyme groups (cytochrome P450 (CYP), carboxylesterases (CES), glutathione S-transferases (GST) and lipases) are up-regulated. Triacylglycerol lipase was significantly over proportionally annotated among up-regulated detox transcripts. We record several up-regulated lipases, GSTe2, two CESs, CYP9A21, CYP6BD6 and CYP9A17 as candidate genes for further H. euphorbiae TPA detoxification analyses. Differential gene expression of the cold acclimation treatment is marked by metabolic depression with enriched Gene Ontology terms among down-regulated transcripts almost exclusively comprising metabolism, aerobic respiration and dissimilative functions. Down-regulated transcripts include energy expensive respiratory proteins like NADH dehydrogenase, cytochrome oxidase and ATP synthase. Gene expression patterns show shifts in carbohydrate metabolism towards cryoprotectant production. The Glycolysis enzymes, G1Pase, A1e, Gpi and an Akr isoform are up-regulated. Glycerol, an osmolyte which lowers the body liquid supercooling point, appears to be the predominant polyol cryoprotectant in H. euphorbiae diapause pupae. Several protein candidates involved in glucose, glycerol, myo-inositol and potentially sorbitol and trehalose synthesis were identified. CONCLUSIONS: A majority of differently expressed transcripts unique for either detoxification or cold hardiness indicates highly specialized functional adaptation which may have evolved from general cell metabolism and stress response.The transcriptome and extracted candidate biomarkers provide a basis for further gene expression studies of physiological processes and adaptive traits in H. euphorbiae.

Deep Super-SAGE transcriptomic analysis of cold acclimation in lentil (Lens culinaris Medik.).

BACKGROUND: Frost is one of the main abiotic stresses limiting plant distribution and crop production. To cope with the stress, plants evolved adaptations known as cold acclimation or chilling tolerance to maximize frost tolerance. Cold acclimation is a progressive acquisition of freezing tolerance by plants subjected to low non-freezing temperatures which subsequently allows them to survive exposure to frost. Lentil is a cool season grain legume that is challenged by winter frost in some areas of its cultivation. RESULTS: To better understand the genetic base of frost tolerance differential gene expression in response to cold acclimation was investigated. Recombinant inbred lines (RILs) from the cross Precoz x WA8649041 were first classified as cold tolerant or cold susceptible according to their response to temperatures between -3 to -15 °C. Then, RILs from both extremes of the response curve were cold acclimated and the leaf transcriptomes of two bulks each of eight frost tolerant and seven cold susceptible RILs were investigated by Deep Super-SAGE transcriptome profiling. Thus, four RNA bulks were analysed: the acclimated susceptible, the acclimated tolerant and the respective controls (non-acclimated susceptible and non-acclimated tolerant). Approximately 16.5 million 26 nucleotide long Super-SAGE tags were sequenced in the four sets (between ~3 and 5.4 millions). In total, 133,077 different unitags, each representing a particular transcript isoform, were identified in these four sets. Tags which showed a significantly different abundance in any of the bulks (fold change ≥4.0 and a significant p-value <0.001) were selected and used to identify the corresponding lentil gene sequence. Three hundred of such lentil sequences were identified. Most of their known homologs coded for glycine-rich, cold and drought-regulated proteins, dormancy-associated proteins, proline-rich proteins (PRPs) and other membrane proteins. These were generally but not exclusively over-expressed in the acclimated tolerant lines. CONCLUSIONS: This set of candidate genes implicated in the response to frost in lentil represents an useful base for deeper and more detailed investigations into this important agronomic trait in future.

Massive analysis of cDNA ends (MACE) reveals a co-segregating candidate gene for LpPg1 stem rust resistance in perennial ryegrass (Lolium perenne).

KEY MESSAGE: Molecular markers including a potential resistance gene co-segregating with the LpPg1 stem rust resistance locus in perennial ryegrass were identified by massive analysis of cDNA ends (MACE) transcriptome profiling. Stem rust caused by Puccinia graminis subsp. graminicola is a severe fungal disease in the forage crop perennial ryegrass and other grasses. The previously identified LpPg1 locus confers efficient resistance against the pathogen. The aim of this study was to identify candidate genes involved in rust resistance and to use them as a resource for the development of molecular markers for LpPg1. To identify such candidates, bulked segregant analysis was combined with NGS-based massive analysis of cDNA ends (MACE) transcriptome profiling. Total RNA was isolated from bulks of infected and non-infected leaf segments from susceptible and resistant genotypes of a full-sibling mapping population and their respective parental lines and MACE was performed. Bioinformatic analysis detected 330 resistance-specific SNPs in 178 transcripts and 341 transcripts that were exclusively expressed in the resistant bulk. The sequences of many of these transcripts were homologous to genes in distinct regions of chromosomes one and four of the model grass Brachypodium distachyon. Of these, 30 were genetically mapped to a 50.8 cM spanning region surrounding the LpPg1 locus. One candidate NBS-LRR gene co-segregated with the resistance locus. Quantitative analysis of gene expression suggests that LpPg1 mediates an efficient resistance mechanism characterized by early recognition of the pathogen, fast defense signaling and rapid induction of antifungal proteins. We demonstrate here that MACE is a cost-efficient, fast and reliable tool that detects polymorphisms for genetic mapping of candidate resistance genes and simultaneously reveals deep insight into the molecular and genetic base of resistance.

Identification of novel small ncRNAs in pollen of tomato.

BACKGROUND: The unprecedented role of sncRNAs in the regulation of pollen biogenesis on both transcriptional and epigenetic levels has been experimentally proven. However, little is known about their global regulation, especially under stress conditions. We used tomato pollen in order to identify pollen stage-specific sncRNAs and their target mRNAs. We further deployed elevated temperatures to discern stress responsive sncRNAs. For this purpose high throughput sncRNA-sequencing as well as Massive Analysis of cDNA Ends (MACE) were performed for three-replicated sncRNAs libraries derived from tomato tetrad, post-meiotic, and mature pollen under control and heat stress conditions. RESULTS: Using the omiRas analysis pipeline we identified known and predicted novel miRNAs as well as sncRNAs from other classes, responsive or not to heat. Differential expression analysis revealed that post-meiotic and mature pollen react most strongly by regulation of the expression of coding and non-coding genomic regions in response to heat. To gain insight to the function of these miRNAs, we predicted targets and annotated them to Gene Ontology terms. This approach revealed that most of them belong to protein binding, transcription, and Serine/Threonine kinase activity GO categories. Beside miRNAs, we observed differential expression of both tRNAs and snoRNAs in tetrad, post-meiotic, and mature pollen when comparing normal and heat stress conditions. CONCLUSIONS: Thus, we describe a global spectrum of sncRNAs expressed in pollen as well as unveiled those which are regulated at specific time-points during pollen biogenesis. We integrated the small RNAs into the regulatory network of tomato heat stress response in pollen.

Characterization and mapping of LanrBo: a locus conferring anthracnose resistance in narrow-leafed lupin (Lupinus angustifolius L.).

KEY MESSAGE: A novel and highly effective source of anthracnose resistance in narrow-leafed lupin was identified. Resistance was shown to be governed by a single dominant locus. Molecular markers have been developed, which can be used for selecting resistant genotypes in lupin breeding. A screening for anthracnose resistance of a set of plant genetic resources of narrow-leafed lupin (Lupinus angustifolius L.) identified the breeding line Bo7212 as being highly resistant to anthracnose (Colletotrichum lupini). Segregation analysis indicated that the resistance of Bo7212 is inherited by a single dominant locus. The corresponding resistance gene was given the designation LanrBo. Previously published molecular anchor markers allowed us to locate LanrBo on linkage group NLL-11 of narrow-leafed lupin. Using information from RNAseq data obtained with inoculated resistant vs. susceptible lupin entries as well as EST-sequence information from the model genome Lotus japonicus, additional SNP and EST markers linked to LanrBo were derived. A bracket of two LanrBo-flanking markers allows for precise marker-assisted selection of the novel resistance gene in narrow-leafed lupin breeding programs.

Allelic diversity in the transcriptomes of contrasting rust-infected genotypes of Lathyrus sativus, a lasting resource for smart breeding.

BACKGROUND: Grass pea (Lathyrus sativus L.) is a valuable resource for potentially durable partial resistance to rust. To gain insight into the resistance mechanism and identify potential resistance genes, we generated the first comprehensive transcriptome assemblies from control and Uromyces pisi inoculated leafs of a susceptible and a partially rust-resistant grass pea genotype by RNA-seq. RESULTS: 134,914 contigs, shared by both libraries, were used to analyse their differential expression in response to rust infection. Functional annotation grouped 60.4% of the contigs present in plant databases (37.8% of total) to 33 main functional categories, being "protein", "RNA", "signalling", "transport" and "stress" the most represented. Transcription profiles revealed considerable differences in regulation of major phytohormone signalling pathways: whereas Salicylic and Abscisic Acid pathways were up-regulated in the resistant genotype, Jasmonate and Ethylene pathways were down-regulated in the susceptible one. As potential Resistance-genes we identified a mildew resistance locus O (MLO)-like gene, and MLO-related transcripts. Also, several pathogenesis-related genes were up-regulated in the resistant and exclusively down regulated in the susceptible genotype. Pathogen effectors identified in both inoculated libraries, as e.g. the rust Rtp1 transcript, may be responsible for the down-regulation of defence-related transcripts. The two genotypes contained 4,892 polymorphic contigs with SNPs unevenly distributed between different functional categories. Protein degradation (29.7%) and signalling receptor kinases (8.2%) were the most diverged, illustrating evolutionary adaptation of grass pea to the host/pathogens arms race. CONCLUSIONS: The vast array of novel, resistance-related genomic information we present here provides a highly valuable resource for future smart breeding approaches in this hitherto under-researched, valuable legume crop.

SymGRASS: a database of sugarcane orthologous genes involved in arbuscular mycorrhiza and root nodule symbiosis.

BACKGROUND: The rationale for gathering information from plants procuring nitrogen through symbiotic interactions controlled by a common genetic program for a sustainable biofuel production is the high energy demanding application of synthetic nitrogen fertilizers. We curated sequence information publicly available for the biofuel plant sugarcane, performed an analysis of the common SYM pathway known to control symbiosis in other plants, and provide results, sequences and literature links as an online database. METHODS: Sugarcane sequences and informations were downloaded from the nucEST database, cleaned and trimmed with seqclean, assembled with TGICL plus translating mapping method, and annotated. The annotation is based on BLAST searches against a local formatted plant Uniprot90 generated with CD-HIT for functional assignment, rpsBLAST to CDD database for conserved domain analysis, and BLAST search to sorghum's for Gene Ontology (GO) assignment. Gene expression was normalized according the Unigene standard, presented as ESTs/100 kb. Protein sequences known in the SYM pathway were used as queries to search the SymGRASS sequence database. Additionally, antimicrobial peptides described in the PhytAMP database served as queries to retrieve and generate expression profiles of these defense genes in the libraries compared to the libraries obtained under symbiotic interactions. RESULTS: We describe the SymGRASS, a database of sugarcane orthologous genes involved in arbuscular mycorrhiza (AM) and root nodule (RN) symbiosis. The database aggregates knowledge about sequences, tissues, organ, developmental stages and experimental conditions, and provides annotation and level of gene expression for sugarcane transcripts and SYM orthologous genes in sugarcane through a web interface. Several candidate genes were found for all nodes in the pathway, and interestingly a set of symbiosis specific genes was found. CONCLUSIONS: The knowledge integrated in SymGRASS may guide studies on molecular, cellular and physiological mechanisms by which sugarcane controls the establishment and efficiency of endophytic associations. We believe that the candidate sequences for the SYM pathway together with the pool of exclusively expressed tentative consensus (TC) sequences are crucial for the design of molecular studies to unravel the mechanisms controlling the establishment of symbioses in sugarcane, ultimately serving as a basis for the improvement of grass crops.

Bulked segregant transcriptome analysis in pea identifies key expression markers for resistance to Peyronellaea pinodes.

Peyronellaea pinodes is a devastating pathogen of pea crop. Quantitative trait loci (QTL) associated with resistance have been identified, as well as genes differentially expressed between resistant and susceptible pea lines. The key question is which of these many genes located into these QTLs, or differentially expressed, are the key genes that distinguish resistant from susceptible plants and could be used as markers. To identify these key genes, in the present study we applied MACE (Massive Analysis of cDNA Ends) -Seq to a whole Recombinant Inbred Line population segregating for resistance to this disease and their parental lines and identified those genes which expression was more correlated with the level of resistance. We also compared gene expression profiles between the most resistant and the most susceptible families of the RIL population. A total of 6780 transcripts were differentially expressed between the parental lines after inoculation. Of them, 803 showed the same expression pattern in the bulks formed by the most resistant and most susceptible RIL families. These genes, showing a consistent expression pattern, could be used as expression markers to distinguish resistant from susceptible plants. The analysis of these genes also discovered the crucial mechanisms acting against P. pinodes.

Short read Illumina data for the de novo assembly of a non-model snail species transcriptome (Radix balthica, Basommatophora, Pulmonata), and a comparison of assembler performance.

BACKGROUND: Until recently, read lengths on the Solexa/Illumina system were too short to reliably assemble transcriptomes without a reference sequence, especially for non-model organisms. However, with read lengths up to 100 nucleotides available in the current version, an assembly without reference genome should be possible. For this study we created an EST data set for the common pond snail Radix balthica by Illumina sequencing of a normalized transcriptome. Performance of three different short read assemblers was compared with respect to: the number of contigs, their length, depth of coverage, their quality in various BLAST searches and the alignment to mitochondrial genes. RESULTS: A single sequencing run of a normalized RNA pool resulted in 16,923,850 paired end reads with median read length of 61 bases. The assemblies generated by VELVET, OASES, and SeqMan NGEN differed in the total number of contigs, contig length, the number and quality of gene hits obtained by BLAST searches against various databases, and contig performance in the mt genome comparison. While VELVET produced the highest overall number of contigs, a large fraction of these were of small size (< 200bp), and gave redundant hits in BLAST searches and the mt genome alignment. The best overall contig performance resulted from the NGEN assembly. It produced the second largest number of contigs, which on average were comparable to the OASES contigs but gave the highest number of gene hits in two out of four BLAST searches against different reference databases. A subsequent meta-assembly of the four contig sets resulted in larger contigs, less redundancy and a higher number of BLAST hits. CONCLUSION: Our results document the first de novo transcriptome assembly of a non-model species using Illumina sequencing data. We show that de novo transcriptome assembly using this approach yields results useful for downstream applications, in particular if a meta-assembly of contig sets is used to increase contig quality. These results highlight the ongoing need for improvements in assembly methodology.

Identification of Tick Ixodes ricinus Midgut Genes Differentially Expressed During the Transmission of Borrelia afzelii Spirochetes Using a Transcriptomic Approach.

Lyme borreliosis is an emerging tick-borne disease caused by spirochetes Borrelia burgdorferi sensu lato. In Europe, Lyme borreliosis is predominantly caused by Borrelia afzelii and transmitted by Ixodes ricinus. Although Borrelia behavior throughout tick development is quite well documented, specific molecular interactions between Borrelia and the tick have not been satisfactorily examined. Here, we present the first transcriptomic study focused on the expression of tick midgut genes regulated by Borrelia. By using massive analysis of cDNA ends (MACE), we searched for tick transcripts expressed differentially in the midgut of unfed, 24h-fed, and fully fed I. ricinus nymphs infected with B. afzelii. In total, we identified 553 upregulated and 530 downregulated tick genes and demonstrated that B. afzelii interacts intensively with the tick. Technical and biological validations confirmed the accuracy of the transcriptome. The expression of five validated tick genes was silenced by RNA interference. Silencing of the uncharacterized protein (GXP_Contig_30818) delayed the infection progress and decreased infection prevalence in the target mice tissues. Silencing of other genes did not significantly affect tick feeding nor the transmission of B. afzelii, suggesting a possible role of these genes rather in Borrelia acquisition or persistence in ticks. Identification of genes and proteins exploited by Borrelia during transmission and establishment in a tick could help the development of novel preventive strategies for Lyme borreliosis.

Epigenetics in non-classical monocytes support their pro-inflammatory gene expression.

Non-classical human monocytes are characterized by high-level expression of cytokines like TNF, but the mechanisms involved are elusive. We have identified miRNAs and CpG-methylation sites that are unique to non-classical monocytes, defined via CD14 and CD16 expression levels. For down-regulated miRNAs that are linked to up-regulated mRNAs the dominant gene ontology term was intracellular signal transduction. This included down-regulated miRNA-20a-5p and miRNA-106b-5p, which both are linked to increased mRNA for the TRIM8 signaling molecule. Methylation analysis revealed 16 hypo-methylated CpG sites upstream of 14 differentially increased mRNAs including 2 sites upstream of TRIM8. Consistent with a positive role in signal transduction, high TRIM8 levels went along with high basal TNF mRNA levels in non-classical monocytes. Since cytokine expression levels in monocytes strongly increase after stimulation with toll-like-receptor ligands, we have analyzed non-classical monocytes (defined via slan expression) after stimulation with lipopolysaccharide (LPS). LPS-stimulated cells continued to have low miRNA-20a and miRNA-106b and high TRIM8 mRNA levels and they showed a 10-fold increase in TNF mRNA. These data suggest that decreased miRNAs and CpG hypo-methylation is linked to enhanced expression of TRIM8 and that this can contribute to the increased TNF levels in non-classical human monocytes.

A combined transcriptomic approach to identify candidates for an anti-tick vaccine blocking B. afzelii transmission.

Ixodes ricinus is the vector for Borrelia afzelii, the predominant cause of Lyme borreliosis in Europe, whereas Ixodes scapularis is the vector for Borrelia burgdorferi in the USA. Transcription of several I. scapularis genes changes in the presence of B. burgdorferi and contributes to successful infection. To what extend B. afzelii influences gene expression in I. ricinus salivary glands is largely unknown. Therefore, we measured expression of uninfected vs. infected tick salivary gland genes during tick feeding using Massive Analysis of cDNA Ends (MACE) and RNAseq, quantifying 26.179 unique transcripts. While tick feeding was the main differentiator, B. afzelii infection significantly affected expression of hundreds of transcripts, including 465 transcripts after 24 h of tick feeding. Validation of the top-20 B. afzelii-upregulated transcripts at 24 h of tick feeding in ten biological genetic distinct replicates showed that expression varied extensively. Three transcripts could be validated, a basic tail protein, a lipocalin and an ixodegrin, and might be involved in B. afzelii transmission. However, vaccination with recombinant forms of these proteins only marginally altered B. afzelii infection in I. ricinus-challenged mice for one of the proteins. Collectively, our data show that identification of tick salivary genes upregulated in the presence of pathogens could serve to identify potential pathogen-blocking vaccine candidates.

Biofilm Formation of Listeria monocytogenes Strains Under Food Processing Environments and Pan-Genome-Wide Association Study.

Concerns about food contamination by Listeria monocytogenes are on the rise with increasing consumption of ready-to-eat foods. Biofilm production of L. monocytogenes is presumed to be one of the ways that confer its increased resistance and persistence in the food chain. In this study, a collection of isolates from foods and food processing environments (FPEs) representing persistent, prevalent, and rarely detected genotypes was evaluated for biofilm forming capacities including adhesion and sessile biomass production under diverse environmental conditions. The quantity of sessile biomass varied according to growth conditions, lineage, serotype as well as genotype but association of clonal complex (CC) 26 genotype with biofilm production was evidenced under cold temperature. In general, relative biofilm productivity of each strain varied inconsistently across growth conditions. Under our experimental conditions, there were no clear associations between biofilm formation efficiency and persistent or prevalent genotypes. Distinct extrinsic factors affected specific steps of biofilm formation. Sudden nutrient deprivation enhanced cellular adhesion while a prolonged nutrient deficiency impeded biofilm maturation. Salt addition increased biofilm production, moreover, nutrient limitation supplemented by salt significantly stimulated biofilm formation. Pan-genome-wide association study (Pan-GWAS) assessed genetic composition with regard to biofilm phenotypes for the first time. The number of reported genes differed depending on the growth conditions and the number of common genes was low. However, a broad overview of the ontology contents revealed similar patterns regardless of the conditions. Functional analysis showed that functions related to transformation/competence and surface proteins including Internalins were highly enriched.

Gene Expression Profiling and Fine Mapping Identifies a Gibberellin 2-Oxidase Gene Co-segregating With the Dominant Dwarfing Gene Ddw1 in Rye (Secale cereale L.).

The gibberellin (GA)-sensitive dwarfing gene Ddw1 provides an opportunity to genetically reduce plant height in rye. Genetic analysis in a population of recombinant inbred lines confirmed a monogenetic dominant inheritance of Ddw1. Significant phenotypic differences in PH between homo- and heterozygotic genotypes indicate an incomplete dominance of Ddw1. De novo transcriptome sequencing of Ddw1 mutant as well as tall genotypes resulted in 113,547 contigs with an average length of 318 bp covering 36.18 Mbp rye DNA. A hierarchical cluster analysis based on individual groups of rye homologs of functionally characterized rice genes controlling morphological or physiological traits including plant height, flowering time, and source activity identified the gene expression profile of stems at the begin of heading to most comprehensively mirror effects of Ddw1. Genome-wide expression profiling identified 186 transcripts differentially expressed between semi-dwarf and tall genotypes in stems. In total, 29 novel markers have been established and mapped to a 27.2 cM segment in the distal part of the long arm of chromosome 5R. Ddw1 could be mapped within a 0.4 cM interval co-segregating with a marker representing the C20-GA2-oxidase gene ScGA2ox12, that is up-regulated in stems of Ddw1 genotypes. The increased expression of ScGA2ox12 observed in semi-dwarf rye as well as structural alterations in transcript sequences associated with the ScGA2ox12 gene implicate, that Ddw1 is a dominant gain-of-function mutant. Integration of the target interval in the wheat reference genome sequence indicated perfect micro-colinearity between the Ddw1 locus and a 831 kb segment on chromosome 5A, which resides inside of a 11.21 Mb interval carrying the GA-sensitive dwarfing gene Rht12 in wheat. The potential of Ddw1 as a breeder's option to improve lodging tolerance in rye is discussed.

First genetic linkage map of Lathyrus cicera based on RNA sequencing-derived markers: Key tool for genetic mapping of disease resistance.

The Lathyrus cicera transcriptome was analysed in response to rust (Uromyces pisi) infection to develop novel molecular breeding tools with potential for genetic mapping of resistance in this robust orphan legume species. One RNA-seq library each was generated from control and rust-inoculated leaves from two L. cicera genotypes with contrasting quantitative resistance, de novo assembled into contigs and sequence polymorphisms were identified. In toto, 19,224 SNPs differentiate the susceptible from the partially resistant genotype's transcriptome. In addition, we developed and tested 341 expressed E-SSR markers from the contigs, of which 60.7% varied between the two L. cicera genotypes. A first L. cicera linkage map was created using part of the developed markers in a RIL population from the cross of the two genotypes. This map contains 307 markers, covered 724.2 cM and is organised in 7 major and 2 minor linkage groups, with an average mapping interval of 2.4 cM. The genic markers also enabled us to compare their position in L. cicera map with the physical position of the same markers mapped on Medicago truncatula genome, highlighting a high macrosyntenic conservation between both species. This study provides a large new set of genic polymorphic molecular markers with potential for mapping rust resistances. It represents the first step towards genomics-assisted precision breeding in L. cicera.

Loss of the Chr16p11.2 ASD candidate gene QPRT leads to aberrant neuronal differentiation in the SH-SY5Y neuronal cell model.

Background: Altered neuronal development is discussed as the underlying pathogenic mechanism of autism spectrum disorders (ASD). Copy number variations of 16p11.2 have recurrently been identified in individuals with ASD. Of the 29 genes within this region, quinolinate phosphoribosyltransferase (QPRT) showed the strongest regulation during neuronal differentiation of SH-SY5Y neuroblastoma cells. We hypothesized a causal relation between this tryptophan metabolism-related enzyme and neuronal differentiation. We thus analyzed the effect of QPRT on the differentiation of SH-SY5Y and specifically focused on neuronal morphology, metabolites of the tryptophan pathway, and the neurodevelopmental transcriptome. Methods: The gene dosage-dependent change of QPRT expression following Chr16p11.2 deletion was investigated in a lymphoblastoid cell line (LCL) of a deletion carrier and compared to his non-carrier parents. Expression of QPRT was tested for correlation with neuromorphology in SH-SY5Y cells. QPRT function was inhibited in SH-SY5Y neuroblastoma cells using (i) siRNA knockdown (KD), (ii) chemical mimicking of loss of QPRT, and (iii) complete CRISPR/Cas9-mediated knock out (KO). QPRT-KD cells underwent morphological analysis. Chemically inhibited and QPRT-KO cells were characterized using viability assays. Additionally, QPRT-KO cells underwent metabolite and whole transcriptome analyses. Genes differentially expressed upon KO of QPRT were tested for enrichment in biological processes and co-regulated gene-networks of the human brain. Results: QPRT expression was reduced in the LCL of the deletion carrier and significantly correlated with the neuritic complexity of SH-SY5Y. The reduction of QPRT altered neuronal morphology of differentiated SH-SY5Y cells. Chemical inhibition as well as complete KO of the gene were lethal upon induction of neuronal differentiation, but not proliferation. The QPRT-associated tryptophan pathway was not affected by KO. At the transcriptome level, genes linked to neurodevelopmental processes and synaptic structures were affected. Differentially regulated genes were enriched for ASD candidates, and co-regulated gene networks were implicated in the development of the dorsolateral prefrontal cortex, the hippocampus, and the amygdala. Conclusions: In this study, QPRT was causally related to in vitro neuronal differentiation of SH-SY5Y cells and affected the regulation of genes and gene networks previously implicated in ASD. Thus, our data suggest that QPRT may play an important role in the pathogenesis of ASD in Chr16p11.2 deletion carriers.

Reprint of: MicroRNA profiling of human intermediate monocytes.

Among the three human monocyte subsets, intermediate CD14++CD16+ monocytes have been characterized as particularly proinflammatory cells in experimental studies and as potential biomarkers of cardiovascular risk in clinical cohorts. To further substantiate the distinct role of intermediate monocytes within human monocyte heterogeneity, we assessed subset-specific expression of miRNAs as central epigenetic regulators of gene expression. We hypothesized that intermediate monocytes have a distinct miRNA profile compared to classical and non-classical monocytes. By using small RNA-seq we analyzed 662 miRNAs in the three monocyte subsets. We identified 38 miRNAs that are differentially expressed in intermediate monocytes compared to both classical and non-classical monocytes with a p value of <10-10, of which two miRNAs - miR-6087 (upregulated) and miR-150-5p (downregulated) - differed in their expression more than ten-fold. Pathway analysis of the 38 differentially expressed miRNAs linked intermediate monocytes to distinct biological processes such as gene regulation, cell differentiation, toll-like receptor signaling as well as antigen processing and presentation. Moreover, differentially expressed miRNAs were connected to those genes that we previously identified as markers of intermediate monocytes. In aggregation, we provide first genome-wide miRNA data in the context of monocyte heterogeneity, which substantiate the concept of monocyte trichotomy in human immunity. The identification of miRNAs that are specific for intermediate monocytes may allow to develop strategies, which particularly target this cell population while sparing the other two subsets.

MicroRNA profiling of human intermediate monocytes.

Among the three human monocyte subsets, intermediate CD14++CD16+ monocytes have been characterized as particularly proinflammatory cells in experimental studies and as potential biomarkers of cardiovascular risk in clinical cohorts. To further substantiate the distinct role of intermediate monocytes within human monocyte heterogeneity, we assessed subset-specific expression of miRNAs as central epigenetic regulators of gene expression. We hypothesized that intermediate monocytes have a distinct miRNA profile compared to classical and non-classical monocytes. By using small RNA-seq we analyzed 662 miRNAs in the three monocyte subsets. We identified 38 miRNAs that are differentially expressed in intermediate monocytes compared to both classical and non-classical monocytes with a p value of <10-10, of which two miRNAs - miR-6087 (upregulated) and miR-150-5p (downregulated) - differed in their expression more than ten-fold. Pathway analysis of the 38 differentially expressed miRNAs linked intermediate monocytes to distinct biological processes such as gene regulation, cell differentiation, toll-like receptor signaling as well as antigen processing and presentation. Moreover, differentially expressed miRNAs were connected to those genes that we previously identified as markers of intermediate monocytes. In aggregation, we provide first genome-wide miRNA data in the context of monocyte heterogeneity, which substantiate the concept of monocyte trichotomy in human immunity. The identification of miRNAs that are specific for intermediate monocytes may allow to develop strategies, which particularly target this cell population while sparing the other two subsets.

Precision medicine for hepatocelluar carcinoma using molecular pattern diagnostics: results from a preclinical pilot study.

The aim of this study was to design a road map for personalizing cancer therapy in hepatocellular carcinoma (HCC) by using molecular pattern diagnostics. As an exploratory study, we investigated molecular patterns of tissues of two tumors from individual HCC patients, which in previous experiments had shown contrasting reactions to the phase 2 transforming growth factor beta receptor 1 inhibitor galunisertib. Cancer-driving molecular patterns encompass - inter alias - altered transcription profiles and somatic mutations in coding regions differentiating tumors from their respective peritumoral tissues and from each other. Massive analysis of cDNA ends and all-exome sequencing demonstrate a highly divergent transcriptional and mutational landscape, respectively, for the two tumors, that offers potential explanations for the tumors contrasting responses to galunisertib. Molecular pattern diagnostics (MPDs) suggest alternative, individual-tumor-specific therapies, which in both cases deviate from the standard sorafenib treatment and from each other. Suggested personalized therapies use kinase inhibitors and immune-focused drugs as well as low-toxicity natural compounds identified using an advanced bioinformatics routine included in the MPD protocol. The MPD pipeline we describe here for the prediction of suitable drugs for treatment of two contrasting HCCs may serve as a blueprint for the design of therapies for various types of cancer.

A Combined Comparative Transcriptomic, Metabolomic, and Anatomical Analyses of Two Key Domestication Traits: Pod Dehiscence and Seed Dormancy in Pea (Pisum sp.).

The origin of the agriculture was one of the turning points in human history, and a central part of this was the evolution of new plant forms, domesticated crops. Seed dispersal and germination are two key traits which have been selected to facilitate cultivation and harvesting of crops. The objective of this study was to analyze anatomical structure of seed coat and pod, identify metabolic compounds associated with water-impermeable seed coat and differentially expressed genes involved in pea seed dormancy and pod dehiscence. Comparative anatomical, metabolomics, and transcriptomic analyses were carried out on wild dormant, dehiscent Pisum elatius (JI64, VIR320) and cultivated, indehiscent Pisum sativum non-dormant (JI92, Cameor) and recombinant inbred lines (RILs). Considerable differences were found in texture of testa surface, length of macrosclereids, and seed coat thickness. Histochemical and biochemical analyses indicated genotype related variation in composition and heterogeneity of seed coat cell walls within macrosclereids. Liquid chromatography-electrospray ionization/mass spectrometry and Laser desorption/ionization-mass spectrometry of separated seed coats revealed significantly higher contents of proanthocyanidins (dimer and trimer of gallocatechin), quercetin, and myricetin rhamnosides and hydroxylated fatty acids in dormant compared to non-dormant genotypes. Bulk Segregant Analysis coupled to high throughput RNA sequencing resulted in identification of 770 and 148 differentially expressed genes between dormant and non-dormant seeds or dehiscent and indehiscent pods, respectively. The expression of 14 selected dormancy-related genes was studied by qRT-PCR. Of these, expression pattern of four genes: porin (MACE-S082), peroxisomal membrane PEX14-like protein (MACE-S108), 4-coumarate CoA ligase (MACE-S131), and UDP-glucosyl transferase (MACE-S139) was in agreement in all four genotypes with Massive analysis of cDNA Ends (MACE) data. In case of pod dehiscence, the analysis of two candidate genes (SHATTERING and SHATTERPROOF) and three out of 20 MACE identified genes (MACE-P004, MACE-P013, MACE-P015) showed down-expression in dorsal and ventral pod suture of indehiscent genotypes. Moreover, MACE-P015, the homolog of peptidoglycan-binding domain or proline-rich extensin-like protein mapped correctly to predicted Dpo1 locus on PsLGIII. This integrated analysis of the seed coat in wild and cultivated pea provides new insight as well as raises new questions associated with domestication and seed dormancy and pod dehiscence.