Oral bacteria accelerate pancreatic cancer development in mice

OBJECTIVE

Epidemiological studies highlight an association between pancreatic ductal adenocarcinoma (PDAC) and oral carriage of the anaerobic bacterium Porphyromonas gingivalis, a species highly linked to periodontal disease. We analysed the potential for P. gingivalis to promote pancreatic cancer development in an animal model and probed underlying mechanisms.

DESIGN

We tracked P. gingivalis bacterial translocation from the oral cavity to the pancreas following administration to mice. To dissect the role of P. gingivalis in PDAC development, we administered bacteria to a genetically engineered mouse PDAC model consisting of inducible acinar cell expression of mutant Kras (Kras +/LSL-G12D; Ptf1a-CreER, iKC mice). These mice were used to study the cooperative effects of Kras mutation and P. gingivalis on the progression of pancreatic intraepithelial neoplasia (PanIN) to PDAC. The direct effects of P. gingivalis on acinar cells and PDAC cell lines were studied in vitro.

RESULTS

P. gingivalis migrated from the oral cavity to the pancreas in mice and can be detected in human PanIN lesions. Repetitive P. gingivalis administration to wild-type mice induced pancreatic acinar-to-ductal metaplasia (ADM), and altered the composition of the intrapancreatic microbiome. In iKC mice, P. gingivalis accelerated PanIN to PDAC progression. In vitro, P. gingivalis infection induced acinar cell ADM markers SOX9 and CK19, and intracellular bacteria protected PDAC cells from reactive oxygen species-mediated cell death resulting from nutrient stress.

CONCLUSION

Taken together, our findings demonstrate a causal role for P. gingivalis in pancreatic cancer development in mice.

Asynchronous Pattern of MAPKs' Activity during Aging of Different Tissues and of Distinct Types of Skeletal Muscle

The MAPK p38α was proposed to be a prominent promoter of skeletal muscle aging. The skeletal muscle tissue is composed of various muscle types, and it is not known if p38α is associated with aging in all of them. It is also not known if p38α is associated with aging of other tissues. JNK and ERK were also proposed to be associated with aging of several tissues. Nevertheless, the pattern of p38α, JNK, and ERK activity during aging was not documented. Here, we documented the levels of phosphorylated/active p38α, Erk1/2, and JNKs in several organs as well as the soleus, tibialis anterior, quadriceps, gastrocnemius, and EDL muscles of 1-, 3-, 6-, 13-, 18-, and 24-month-old mice. We report that in most tissues and skeletal muscles, the MAPKs' activity does not change in the course of aging. In most tissues and muscles, p38α is in fact active at younger ages. The quadriceps and the lungs are exceptions, where p38α is significantly active only in mice 13 months old or older. Curiously, levels of active JNK and ERKs are also elevated in aged lungs and quadriceps. RNA-seq analysis of the quadriceps during aging revealed downregulation of proteins related to the extra-cellular matrix (ECM) and ERK signaling. A panel of mRNAs encoding cell cycle inhibitors and senescence-associated proteins, considered to be aging markers, was not found to be elevated. It seems that the pattern of MAPKs' activation in aging, as well as expression of known 'aging' components, are tissue- and muscle type-specific, supporting a notion that the process of aging is tissue- and even cell-specific.

A short peptide protects from age-onset proteotoxicity

Aberrant protein aggregation jeopardizes cellular functionality and underlies the development of a myriad of late-onset maladies including Alzheimer's disease (AD) and Huntington's disease (HD). Accordingly, molecules that mitigate the toxicity of hazardous protein aggregates are of great interest as potential future therapeutics. Here we asked whether a small peptide, composed of five amino acids (5MER peptide) that was derived from the human pro-inflammatory CD44 protein, could protect model nematodes from the toxicity of aggregative proteins that underlie the development of neurodegenerative disorders in humans. We found that the 5MER peptide mitigates the toxicity that stems from both; the AD-causing Aβ peptide and a stretch of poly-glutamine that is accountable for the development of several disorders including HD, while minimally affecting lifespan. This protection was dependent on the activity of aging-regulating transcription factors and associated with enhanced Aβ and polyQ35-YFP aggregation. A transcriptomic analysis unveiled that the peptide modifies signaling pathways, thereby modulating the expression of various genes, including these, which are known as protein homeostasis (proteostasis) regulators such as txt-13 and modifiers of proteasome activity. The knockdown of txt-13 protects worms from proteotoxicity to the same extent as the 5MER peptide, suggesting that the peptide activates the transcellular chaperone signaling to promote proteostasis. Together, our results propose that the 5MER peptide should be considered as a component of future therapeutic cocktails for the treatment of neurodegenerative maladies.

Systematic identification of gene combinations to target in innate immune cells to enhance T cell activation

Genetic engineering of immune cells has opened new avenues for improving their functionality but it remains a challenge to pinpoint which genes or combination of genes are the most beneficial to target. Here, we conduct High Multiplicity of Perturbations and Cellular Indexing of Transcriptomes and Epitopes (HMPCITE-seq) to find combinations of genes whose joint targeting improves antigen-presenting cell activity and enhances their ability to activate T cells. Specifically, we perform two genome-wide CRISPR screens in bone marrow dendritic cells and identify negative regulators of CD86, that participate in the co-stimulation programs, including Chd4, Stat5b, Egr2, Med12, and positive regulators of PD-L1, that participate in the co-inhibitory programs, including Sptlc2, Nckap1l, and Pi4kb. To identify the genetic interactions between top-ranked genes and find superior combinations to target, we perform high-order Perturb-Seq experiments and we show that targeting both Cebpb and Med12 results in a better phenotype compared to the single perturbations or other combinations of perturbations.

Isoforms of the TAL1 transcription factor have different roles in hematopoiesis and cell growth

T-cell acute lymphoblastic leukemia (T-ALL) protein 1 (TAL1) is a central transcription factor in hematopoiesis. The timing and level of TAL1 expression orchestrate the differentiation to specialized blood cells and its overexpression is a common cause of T-ALL. Here, we studied the 2 protein isoforms of TAL1, short and long, which are generated by the use of alternative promoters as well as by alternative splicing. We analyzed the expression of each isoform by deleting an enhancer or insulator, or by opening chromatin at the enhancer location. Our results show that each enhancer promotes expression from a specific TAL1 promoter. Expression from a specific promoter gives rise to a unique 5' UTR with differential regulation of translation. Moreover, our study suggests that the enhancers regulate TAL1 exon 3 alternative splicing by inducing changes in the chromatin at the splice site, which we demonstrate is mediated by KMT2B. Furthermore, our results indicate that TAL1-short binds more strongly to TAL1 E-protein partners and functions as a stronger transcription factor than TAL1-long. Specifically TAL1-short has a unique transcription signature promoting apoptosis. Finally, when we expressed both isoforms in mice bone marrow, we found that while overexpression of both isoforms prevents lymphoid differentiation, expression of TAL-short alone leads to hematopoietic stem cell exhaustion. Furthermore, we found that TAL1-short promoted erythropoiesis and reduced cell survival in the CML cell line K562. While TAL1 and its partners are considered promising therapeutic targets in the treatment of T-ALL, our results show that TAL1-short could act as a tumor suppressor and suggest that altering TAL1 isoform's ratio could be a preferred therapeutic approach.

UMSBP2 is chromatin remodeler that functions in regulation of gene expression and suppression of antigenic variation in trypanosomes

Universal Minicircle Sequence binding proteins (UMSBPs) are CCHC-type zinc-finger proteins that bind the single-stranded G-rich UMS sequence, conserved at the replication origins of minicircles in the kinetoplast DNA, the mitochondrial genome of kinetoplastids. Trypanosoma brucei UMSBP2 has been recently shown to colocalize with telomeres and to play an essential role in chromosome end protection. Here we report that TbUMSBP2 decondenses in vitro DNA molecules, which were condensed by core histones H2B, H4 or linker histone H1. DNA decondensation is mediated via protein-protein interactions between TbUMSBP2 and these histones, independently of its previously described DNA binding activity. Silencing of the TbUMSBP2 gene resulted in a significant decrease in the disassembly of nucleosomes in T. brucei chromatin, a phenotype that could be reverted, by supplementing the knockdown cells with TbUMSBP2. Transcriptome analysis revealed that silencing of TbUMSBP2 affects the expression of multiple genes in T. brucei, with a most significant effect on the upregulation of the subtelomeric variant surface glycoproteins (VSG) genes, which mediate the antigenic variation in African trypanosomes. These observations suggest that UMSBP2 is a chromatin remodeling protein that functions in the regulation of gene expression and plays a role in the control of antigenic variation in T. brucei.

A reversible mutation in a genomic hotspot saves bacterial swarms from extinction

Microbial adaptation to changing environmental conditions is frequently mediated by hypermutable sequences. Here we demonstrate that such a hypermutable hotspot within a gene encoding a flagellar unit of Paenibacillus glucanolyticus generated spontaneous non-swarming mutants with increased stress resistance. These mutants, which survived conditions that eliminated wild-type cultures, could be carried by their swarming siblings when the colony spread, consequently increasing their numbers at the spreading edge. Of interest, the hypermutable nature of the aforementioned sequence enabled the non-swarming mutants to serve as "seeds" for a new generation of wild-type cells through reversion of the mutation. Using a mathematical model, we examined the survival dynamics of P. glucanolyticus colonies under fluctuating environments. Our experimental and theoretical results suggest that the non-swarming, stress-resistant mutants can save the colony from extinction. Notably, we identified this hypermutable sequence in flagellar genes of additional Paenibacillus species, suggesting that this phenomenon could be wide-spread and ecologically important.

Chronic Intermittent Hypoxia during Sleep Causes Browning of Interscapular Adipose Tissue Accompanied by Local Insulin Resistance in Mice

Obstructive sleep apnea (OSA) is a highly prevalent condition, characterized by intermittent hypoxia (IH), sleep disruption, and altered autonomic nervous system function. OSA has been independently associated with dyslipidemia, insulin resistance, and metabolic syndrome. Brown adipose tissue (BAT) has been suggested as a modulator of systemic glucose tolerance through adaptive thermogenesis. Reductions in BAT mass have been associated with obesity and metabolic syndrome. No studies have systematically characterized the effects of chronic IH on BAT. Thus, we aimed to delineate IH effects on BAT and concomitant metabolic changes. C57BL/6J 8-week-old male mice were randomly assigned to IH during sleep (alternating 90 s cycles of 6.5% F_IO₂ followed by 21% F_IO₂) or normoxia (room air, RA) for 10 weeks. Mice were subjected to glucose tolerance testing and ¹⁸F-FDG PET-MRI towards the end of the exposures followed by BAT tissues analyses for morphological and global transcriptomic changes. Animals exposed to IH were glucose intolerant despite lower total body weight and adiposity. BAT tissues in IH-exposed mice demonstrated characteristic changes associated with "browning"-smaller lipids, increased vascularity, and a trend towards higher protein levels of UCP1. Conversely, mitochondrial DNA content and protein levels of respiratory chain complex III were reduced. Pro-inflammatory macrophages were more abundant in IH-exposed BAT. Transcriptomic analysis revealed increases in fatty acid oxidation and oxidative stress pathways in IH-exposed BAT, along with a reduction in pathways related to myogenesis, hypoxia, and IL-4 anti-inflammatory response. Functionally, IH-exposed BAT demonstrated reduced absorption of glucose on PET scans and reduced phosphorylation of AKT in response to insulin. Current studies provide initial evidence for the presence of a maladaptive response of interscapular BAT in response to chronic IH mimicking OSA, resulting in a paradoxical divergence, namely, BAT browning but tissue-specific and systemic insulin resistance. We postulate that oxidative stress, mitochondrial dysfunction, and inflammation may underlie these dichotomous outcomes in BAT.

P450 oxidoreductase regulates barrier maturation by mediating retinoic acid metabolism in a model of the human BBB

The blood-brain barrier (BBB) selectively regulates the entry of molecules into the central nervous system (CNS). A crosstalk between brain microvascular endothelial cells (BMECs) and resident CNS cells promotes the acquisition of functional tight junctions (TJs). Retinoic acid (RA), a key signaling molecule during embryonic development, is used to enhance in vitro BBB models’ functional barrier properties. However, its physiological relevance and affected pathways are not fully understood. P450 oxidoreductase (POR) regulates the enzymatic activity of microsomal cytochromes. POR-deficient (PORD) patients display impaired steroid homeostasis and cognitive disabilities. Here, we used both patient-specific POR-deficient and CRISPR-Cas9-mediated POR-depleted induced pluripotent stem cell (iPSC)-derived BMECs (iBMECs) to study the role of POR in the acquisition of functional barrier properties. We demonstrate that POR regulates cellular RA homeostasis and that POR deficiency leads to the accumulation of RA within iBMECs, resulting in the impaired acquisition of TJs and, consequently, to dysfunctional development of barrier properties.

•

Retinoic acid (RA) promotes functional barrier properties

•

POR-deficient iPS-brain endothelial-like cells display impaired barrier development

•

POR mediates CYP26-dependent cellular RA catabolism

•

RA accumulation induces a pro-inflammatory response

Single-cell transcriptomics reveals a senescence-associated IL-6/CCR6 axis driving radiodermatitis

Irradiation‐induced alopecia and dermatitis (IRIAD) are two of the most visually recognized complications of radiotherapy, of which the molecular and cellular basis remains largely unclear. By combining scRNA‐seq analysis of whole skin‐derived irradiated cells with genetic ablation and molecular inhibition studies, we show that senescence‐associated IL‐6 and IL‐1 signaling, together with IL‐17 upregulation and CCR6⁺‐mediated immune cell migration, are crucial drivers of IRIAD. Bioinformatics analysis colocalized irradiation‐induced IL‐6 signaling with senescence pathway upregulation largely within epidermal hair follicles, basal keratinocytes, and dermal fibroblasts. Loss of cytokine signaling by genetic ablation in IL‐6^−/− or IL‐1R^−/− mice, or by molecular blockade, strongly ameliorated IRIAD, as did deficiency of CCL20/CCR6‐mediated immune cell migration in CCR6^−/− mice. Moreover, IL‐6 deficiency strongly reduced IL‐17, IL‐22, CCL20, and CCR6 upregulation, whereas CCR6 deficiency reciprocally diminished IL‐6, IL‐17, CCL3, and MHC upregulation, suggesting that proximity‐dependent cellular cross talk promotes IRIAD. Therapeutically, topical application of Janus kinase blockers or inhibition of T‐cell activation by cyclosporine effectively reduced IRIAD, suggesting the potential of targeted approaches for the treatment of dermal side effects in radiotherapy patients.

Rapid Atrial Pacing Promotes Atrial Fibrillation Substrate in Unanesthetized Instrumented Rats

Aim

The self-perpetuating nature of atrial fibrillation (AF) has been a subject of intense research in large mammalian models exposed to rapid atrial pacing (RAP). Recently, rodents are increasingly used to gain insight into the pathophysiology of AF. However, little is known regarding the effects of RAP on the atria of rats and mice. Using an implantable device for electrophysiological studies in rodents, we examined on a daily basis, the effects of continuous RAP on the developed AF substrate of unanesthetized rats and mice.

Methods and Results

Aggressive burst pacing did not induce AF at baseline in the large majority of rodents, but repeatedly induced AF episodes in rats exposed to RAP for more than 2 days. A microarray study of left atrial tissue from rats exposed to RAP for 2 days vs. control pacing identified 304 differentially expressed genes. Enrichment analysis and comparison with a dataset of atrial tissue from AF patients revealed indications of increased carbohydrate metabolism and changes in pathways that are thought to play critical roles in human AF, including TGF-beta and IL-6 signaling. Among 19 commonly affected genes in comparison with human AF, downregulation of FOXP1 and upregulation of the KCNK2 gene encoding the Kir2.1 potassium channel were conspicuous findings, suggesting NFAT activation. Further results included reduced expression of MIR-26 and MIR-101, which is in line with NFAT activation.

Conclusion

Our results demonstrate electrophysiological evidence for AF promoting effects of RAP in rats and several molecular similarities between the effects of RAP in large and small mammalian models.

Complete Assembly of the Genome of an Acidovorax citrulli Strain Reveals a Naturally Occurring Plasmid in This Species

Acidovorax citrulli is the causal agent of bacterial fruit blotch (BFB), a serious threat to cucurbit crop production worldwide. Based on genetic and phenotypic properties, A. citrulli strains are divided into two major groups: group I strains have been generally isolated from melon and other non-watermelon cucurbits, while group II strains are closely associated with watermelon. In a previous study, we reported the genome of the group I model strain, M6. At that time, the M6 genome was sequenced by MiSeq Illumina technology, with reads assembled into 139 contigs. Here, we report the assembly of the M6 genome following sequencing with PacBio technology. This approach not only allowed full assembly of the M6 genome, but it also revealed the occurrence of a ∼53 kb plasmid. The M6 plasmid, named pACM6, was further confirmed by plasmid extraction, Southern-blot analysis of restricted fragments and obtention of M6-derivative cured strains. pACM6 occurs at low copy numbers (average of ∼4.1 ± 1.3 chromosome equivalents) in A. citrulli M6 and contains 63 open reading frames (ORFs), most of which (55.6%) encoding hypothetical proteins. The plasmid contains several genes encoding type IV secretion components, and typical plasmid-borne genes involved in plasmid maintenance, replication and transfer. The plasmid also carries an operon encoding homologs of a Fic-VbhA toxin-antitoxin (TA) module. Transcriptome data from A. citrulli M6 revealed that, under the tested conditions, the genes encoding the components of this TA system are among the highest expressed genes in pACM6. Whether this TA module plays a role in pACM6 maintenance is still to be determined. Leaf infiltration and seed transmission assays revealed that, under tested conditions, the loss of pACM6 did not affect the virulence of A. citrulli M6. We also show that pACM6 or similar plasmids are present in several group I strains, but absent in all tested group II strains of A. citrulli.

GeneHancer: genome-wide integration of enhancers and target genes in GeneCards

A major challenge in understanding gene regulation is the unequivocal identification of enhancer elements and uncovering their connections to genes. We present GeneHancer, a novel database of human enhancers and their inferred target genes, in the framework of GeneCards. First, we integrated a total of 434 000 reported enhancers from four different genome-wide databases: the Encyclopedia of DNA Elements (ENCODE), the Ensembl regulatory build, the functional annotation of the mammalian genome (FANTOM) project and the VISTA Enhancer Browser. Employing an integration algorithm that aims to remove redundancy, GeneHancer portrays 285 000 integrated candidate enhancers (covering 12.4% of the genome), 94 000 of which are derived from more than one source, and each assigned an annotation-derived confidence score. GeneHancer subsequently links enhancers to genes, using: tissue co-expression correlation between genes and enhancer RNAs, as well as enhancer-targeted transcription factor genes; expression quantitative trait loci for variants within enhancers; and capture Hi-C, a promoter-specific genome conformation assay. The individual scores based on each of these four methods, along with gene–enhancer genomic distances, form the basis for GeneHancer’s combinatorial likelihood-based scores for enhancer–gene pairing. Finally, we define ‘elite’ enhancer–gene relations reflecting both a high-likelihood enhancer definition and a strong enhancer–gene association.

GeneHancer predictions are fully integrated in the widely used GeneCards Suite, whereby candidate enhancers and their annotations are displayed on every relevant GeneCard. This assists in the mapping of non-coding variants to enhancers, and via the linked genes, forms a basis for variant–phenotype interpretation of whole-genome sequences in health and disease.

Database URL:http://www.genecards.org/

Rational confederation of genes and diseases: NGS interpretation via GeneCards, MalaCards and VarElect

Background

A key challenge in the realm of human disease research is next generation sequencing (NGS) interpretation, whereby identified filtered variant-harboring genes are associated with a patient’s disease phenotypes. This necessitates bioinformatics tools linked to comprehensive knowledgebases. The GeneCards suite databases, which include GeneCards (human genes), MalaCards (human diseases) and PathCards (human pathways) together with additional tools, are presented with the focus on MalaCards utility for NGS interpretation as well as for large scale bioinformatic analyses.

Results

VarElect, our NGS interpretation tool, leverages the broad information in the GeneCards suite databases. MalaCards algorithms unify disease-related terms and annotations from 69 sources. Further, MalaCards defines hierarchical relatedness—aliases, disease families, a related diseases network, categories and ontological classifications. GeneCards and MalaCards delineate and share a multi-tiered, scored gene-disease network, with stringency levels, including the definition of elite status—high quality gene-disease pairs, coming from manually curated trustworthy sources, that includes 4500 genes for 8000 diseases. This unique resource is key to NGS interpretation by VarElect. VarElect, a comprehensive search tool that helps infer both direct and indirect links between genes and user-supplied disease/phenotype terms, is robustly strengthened by the information found in MalaCards. The indirect mode benefits from GeneCards’ diverse gene-to-gene relationships, including SuperPaths—integrated biological pathways from 12 information sources. We are currently adding an important information layer in the form of “disease SuperPaths”, generated from the gene-disease matrix by an algorithm similar to that previously employed for biological pathway unification. This allows the discovery of novel gene-disease and disease–disease relationships. The advent of whole genome sequencing necessitates capacities to go beyond protein coding genes. GeneCards is highly useful in this respect, as it also addresses 101,976 non-protein-coding RNA genes. In a more recent development, we are currently adding an inclusive map of regulatory elements and their inferred target genes, generated by integration from 4 resources.

Conclusions

MalaCards provides a rich big-data scaffold for in silico biomedical discovery within the gene-disease universe. VarElect, which depends significantly on both GeneCards and MalaCards power, is a potent tool for supporting the interpretation of wet-lab experiments, notably NGS analyses of disease. The GeneCards suite has thus transcended its 2-decade role in biomedical research, maturing into a key player in clinical investigation.

Dormancy in Embryos: Insight from Hydrated Encysted Embryos of an Aquatic Invertebrate

Numerous aquatic invertebrates remain dormant for decades in a hydrated state as encysted embryos. In search for functional pathways associated with this form of dormancy, we used label-free quantitative proteomics to compare the proteomes of hydrated encysted dormant embryos (resting eggs; RE) with nondormant embryos (amictic eggs; AM) of the rotifer Brachionus plicatilisA total of 2631 proteins were identified in rotifer eggs. About 62% proteins showed higher abundance in AM relative to RE (Fold Change>3; p = 0.05). Proteins belonging to numerous putative functional pathways showed dramatic changes during dormancy. Most striking were changes in the mitochondria indicating an impeded metabolism. A comparison between the abundance of proteins and their corresponding transcript levels, revealed higher concordance for RE than for AM. Surprisingly, numerous highly abundant dormancy related proteins show corresponding high mRNA levels in metabolically inactive RE. As these mRNAs and proteins degrade at the time of exit from dormancy they may serve as a source of nucleotides and amino acids during the exit from dormancy. Because proteome analyses point to a similarity in functional pathways of hydrated RE and desiccated life forms, REs were dried. Similar hatching and reproductive rates were found for wet and dried REs, suggesting analogous pathways for long-term survival in wet or dry forms. Analysis by KEGG pathways revealed a few general strategies for dormancy, proposing an explanation for the low transcriptional similarity among dormancies across species, despite the resemblance in physiological phenotypes.

GeneHancer: genome-wide integration of enhancers and target genes in GeneCards

A major challenge in understanding gene regulation is the unequivocal identification of enhancer elements and uncovering their connections to genes. We present GeneHancer, a novel database of human enhancers and their inferred target genes, in the framework of GeneCards. First, we integrated a total of 434 000 reported enhancers from four different genome-wide databases: the Encyclopedia of DNA Elements (ENCODE), the Ensembl regulatory build, the functional annotation of the mammalian genome (FANTOM) project and the VISTA Enhancer Browser. Employing an integration algorithm that aims to remove redundancy, GeneHancer portrays 285 000 integrated candidate enhancers (covering 12.4% of the genome), 94 000 of which are derived from more than one source, and each assigned an annotation-derived confidence score. GeneHancer subsequently links enhancers to genes, using: tissue co-expression correlation between genes and enhancer RNAs, as well as enhancer-targeted transcription factor genes; expression quantitative trait loci for variants within enhancers; and capture Hi-C, a promoter-specific genome conformation assay. The individual scores based on each of these four methods, along with gene–enhancer genomic distances, form the basis for GeneHancer’s combinatorial likelihood-based scores for enhancer–gene pairing. Finally, we define ‘elite’ enhancer–gene relations reflecting both a high-likelihood enhancer definition and a strong enhancer–gene association. GeneHancer predictions are fully integrated in the widely used GeneCards Suite, whereby candidate enhancers and their annotations are displayed on every relevant GeneCard. This assists in the mapping of non-coding variants to enhancers, and via the linked genes, forms a basis for variant–phenotype interpretation of whole-genome sequences in health and disease.

DATABASE URL

http://www.genecards.org/.

GeneAnalytics: An Integrative Gene Set Analysis Tool for Next Generation Sequencing, RNAseq and Microarray Data

Postgenomics data are produced in large volumes by life sciences and clinical applications of novel omics diagnostics and therapeutics for precision medicine. To move from “data-to-knowledge-to-innovation,” a crucial missing step in the current era is, however, our limited understanding of biological and clinical contexts associated with data. Prominent among the emerging remedies to this challenge are the gene set enrichment tools. This study reports on GeneAnalytics™ (geneanalytics.genecards.org), a comprehensive and easy-to-apply gene set analysis tool for rapid contextualization of expression patterns and functional signatures embedded in the postgenomics Big Data domains, such as Next Generation Sequencing (NGS), RNAseq, and microarray experiments. GeneAnalytics' differentiating features include in-depth evidence-based scoring algorithms, an intuitive user interface and proprietary unified data. GeneAnalytics employs the LifeMap Science's GeneCards suite, including the GeneCards®—the human gene database; the MalaCards—the human diseases database; and the PathCards—the biological pathways database. Expression-based analysis in GeneAnalytics relies on the LifeMap Discovery®—the embryonic development and stem cells database, which includes manually curated expression data for normal and diseased tissues, enabling advanced matching algorithm for gene–tissue association. This assists in evaluating differentiation protocols and discovering biomarkers for tissues and cells. Results are directly linked to gene, disease, or cell “cards” in the GeneCards suite. Future developments aim to enhance the GeneAnalytics algorithm as well as visualizations, employing varied graphical display items. Such attributes make GeneAnalytics a broadly applicable postgenomics data analyses and interpretation tool for translation of data to knowledge-based innovation in various Big Data fields such as precision medicine, ecogenomics, nutrigenomics, pharmacogenomics, vaccinomics, and others yet to emerge on the postgenomics horizon.

MalaCards: an amalgamated human disease compendium with diverse clinical and genetic annotation and structured search

The MalaCards human disease database (http://www.malacards.org/) is an integrated compendium of annotated diseases mined from 68 data sources. MalaCards has a web card for each of ∼20 000 disease entries, in six global categories. It portrays a broad array of annotation topics in 15 sections, including Summaries, Symptoms, Anatomical Context, Drugs, Genetic Tests, Variations and Publications. The Aliases and Classifications section reflects an algorithm for disease name integration across often-conflicting sources, providing effective annotation consolidation. A central feature is a balanced Genes section, with scores reflecting the strength of disease-gene associations. This is accompanied by other gene-related disease information such as pathways, mouse phenotypes and GO-terms, stemming from MalaCards’ affiliation with the GeneCards Suite of databases. MalaCards’ capacity to inter-link information from complementary sources, along with its elaborate search function, relational database infrastructure and convenient data dumps, allows it to tackle its rich disease annotation landscape, and facilitates systems analyses and genome sequence interpretation. MalaCards adopts a ‘flat’ disease-card approach, but each card is mapped to popular hierarchical ontologies (e.g. International Classification of Diseases, Human Phenotype Ontology and Unified Medical Language System) and also contains information about multi-level relations among diseases, thereby providing an optimal tool for disease representation and scrutiny.

VarElect: the phenotype-based variation prioritizer of the GeneCards Suite

Background

Next generation sequencing (NGS) provides a key technology for deciphering the genetic underpinnings of human diseases. Typical NGS analyses of a patient depict tens of thousands non-reference coding variants, but only one or very few are expected to be significant for the relevant disorder. In a filtering stage, one employs family segregation, rarity in the population, predicted protein impact and evolutionary conservation as a means for shortening the variation list. However, narrowing down further towards culprit disease genes usually entails laborious seeking of gene-phenotype relationships, consulting numerous separate databases. Thus, a major challenge is to transition from the few hundred shortlisted genes to the most viable disease-causing candidates.

Results

We describe a novel tool, VarElect (http://ve.genecards.org), a comprehensive phenotype-dependent variant/gene prioritizer, based on the widely-used GeneCards, which helps rapidly identify causal mutations with extensive evidence. The GeneCards suite offers an effective and speedy alternative, whereby >120 gene-centric automatically-mined data sources are jointly available for the task. VarElect cashes on this wealth of information, as well as on GeneCards’ powerful free-text Boolean search and scoring capabilities, proficiently matching variant-containing genes to submitted disease/symptom keywords. The tool also leverages the rich disease and pathway information of MalaCards, the human disease database, and PathCards, the unified pathway (SuperPaths) database, both within the GeneCards Suite. The VarElect algorithm infers direct as well as indirect links between genes and phenotypes, the latter benefitting from GeneCards’ diverse gene-to-gene data links in GenesLikeMe. Finally, our tool offers an extensive gene-phenotype evidence portrayal (“MiniCards”) and hyperlinks to the parent databases.

Conclusions

We demonstrate that VarElect compares favorably with several often-used NGS phenotyping tools, thus providing a robust facility for ranking genes, pointing out their likelihood to be related to a patient’s disease. VarElect’s capacity to automatically process numerous NGS cases, either in stand-alone format or in VCF-analyzer mode (TGex and VarAnnot), is indispensable for emerging clinical projects that involve thousands of whole exome/genome NGS analyses.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2722-2) contains supplementary material, which is available to authorized users.

The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses

GeneCards, the human gene compendium, enables researchers to effectively navigate and inter-relate the wide universe of human genes, diseases, variants, proteins, cells, and biological pathways. Our recently launched Version 4 has a revamped infrastructure facilitating faster data updates, better-targeted data queries, and friendlier user experience. It also provides a stronger foundation for the GeneCards suite of companion databases and analysis tools. Improved data unification includes gene-disease links via MalaCards and merged biological pathways via PathCards, as well as drug information and proteome expression. VarElect, another suite member, is a phenotype prioritizer for next-generation sequencing, leveraging the GeneCards and MalaCards knowledgebase. It automatically infers direct and indirect scored associations between hundreds or even thousands of variant-containing genes and disease phenotype terms. VarElect's capabilities, either independently or within TGex, our comprehensive variant analysis pipeline, help prepare for the challenge of clinical projects that involve thousands of exome/genome NGS analyses. © 2016 by John Wiley & Sons, Inc.

A crayfish molar tooth protein with putative mineralized exoskeletal chitinous matrix properties

Some crustaceans possess exoskeletons that are reinforced with calcium carbonate. In the crayfish Cherax quadricarinatus, the molar tooth, which is part of the mandibular exoskeleton, contains an unusual crystalline enamel-like apatite layer. As this layer resembles vertebrate enamel in composition and function, it offers an interesting example of convergent evolution. Unlike other parts of the crayfish exoskeleton, which is periodically shed and regenerated during the molt cycle, molar mineral deposition takes place during the pre-molt stage. The molar mineral composition transforms continuously from fluorapatite through amorphous calcium phosphate to amorphous calcium carbonate and is mounted on chitin. The process of crayfish molar formation is entirely extracellular and presumably controlled by proteins, lipids, polysaccharides, low-molecular weight molecules and calcium salts. We have identified a novel molar protein termed Cq-M15 from C. quadricarinatus and cloned its transcript from the molar-forming epithelium. Its transcript and differential expression were confirmed by a next-generation sequencing library. The predicted acidic pI of Cq-M15 suggests its possible involvement in mineral arrangement. Cq-M15 is expressed in several exoskeletal tissues at pre-molt and its silencing is lethal. Like other arthropod cuticular proteins, Cq-M15 possesses a chitin-binding Rebers-Riddiford domain, with a recombinant version of the protein found to bind chitin. Cq-M15 was also found to interact with calcium ions in a concentration-dependent manner. This latter property might make Cq-M15 useful for bone and dental regenerative efforts. We suggest that, in the molar tooth, this protein might be involved in calcium phosphate and/or carbonate precipitation.

Early and late response of Nematostella vectensis transcriptome to heavy metals

Environmental contamination from heavy metals poses a global concern for the marine environment, as heavy metals are passed up the food chain and persist in the environment long after the pollution source is contained. Cnidarians play an important role in shaping marine ecosystems, but environmental pollution profoundly affects their vitality. Among the cnidarians, the sea anemone Nematostella vectensis is an advantageous model for addressing questions in molecular ecology and toxicology as it tolerates extreme environments and its genome has been published. Here, we employed a transcriptome-wide RNA-Seq approach to analyse N. vectensis molecular defence mechanisms against four heavy metals: Hg, Cu, Cd and Zn. Altogether, more than 4800 transcripts showed significant changes in gene expression. Hg had the greatest impact on up-regulating transcripts, followed by Cu, Zn and Cd. We identified, for the first time in Cnidaria, co-up-regulation of immediate-early transcription factors such as Egr1, AP1 and NF-κB. Time-course analysis of these genes revealed their early expression as rapidly as one hour after exposure to heavy metals, suggesting that they may complement or substitute for the roles of the metal-mediating Mtf1 transcription factor. We further characterized the regulation of a large array of stress-response gene families, including Hsp, ABC, CYP members and phytochelatin synthase, that may regulate synthesis of the metal-binding phytochelatins instead of the metallothioneins that are absent from Cnidaria genome. This study provides mechanistic insight into heavy metal toxicity in N. vectensis and sheds light on ancestral stress adaptations.

Metabolite and transcript profiling of berry skin during fruit development elucidates differential regulation between Cabernet Sauvignon and Shiraz cultivars at branching points in the polyphenol pathway

BACKGROUND

Grapevine berries undergo complex biochemical changes during fruit maturation, many of which are dependent upon the variety and its environment. In order to elucidate the varietal dependent developmental regulation of primary and specialized metabolism, berry skins of Cabernet Sauvignon and Shiraz were subjected to gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-mass spectrometry (LC-MS) based metabolite profiling from pre-veraison to harvest. The generated dataset was augmented with transcript profiling using RNAseq.

RESULTS

The analysis of the metabolite data revealed similar developmental patterns of change in primary metabolites between the two cultivars. Nevertheless, towards maturity the extent of change in the major organic acid and sugars (i.e. sucrose, trehalose, malate) and precursors of aromatic and phenolic compounds such as quinate and shikimate was greater in Shiraz compared to Cabernet Sauvignon. In contrast, distinct directional projections on the PCA plot of the two cultivars samples towards maturation when using the specialized metabolite profiles were apparent, suggesting a cultivar-dependent regulation of the specialized metabolism. Generally, Shiraz displayed greater upregulation of the entire polyphenol pathway and specifically higher accumulation of piceid and coumaroyl anthocyanin forms than Cabernet Sauvignon from veraison onwards. Transcript profiling revealed coordinated increased transcript abundance for genes encoding enzymes of committing steps in the phenylpropanoid pathway. The anthocyanin metabolite profile showed F3'5'H-mediated delphinidin-type anthocyanin enrichment in both varieties towards maturation, consistent with the transcript data, indicating that the F3'5'H-governed branching step dominates the anthocyanin profile at late berry development. Correlation analysis confirmed the tightly coordinated metabolic changes during development, and suggested a source-sink relation between the central and specialized metabolism, stronger in Shiraz than Cabernet Sauvignon. RNAseq analysis also revealed that the two cultivars exhibited distinct pattern of changes in genes related to abscisic acid (ABA) biosynthesis enzymes.

CONCLUSIONS

Compared with CS, Shiraz showed higher number of significant correlations between metabolites, which together with the relatively higher expression of flavonoid genes supports the evidence of increased accumulation of coumaroyl anthocyanins in that cultivar. Enhanced stress related metabolism, e.g. trehalose, stilbene and ABA in Shiraz berry-skin are consistent with its relatively higher susceptibility to environmental cues.

Inositol-related gene knockouts mimic lithium's effect on mitochondrial function

The inositol-depletion hypothesis proposes that lithium attenuates phosphatidylinositol signaling. Knockout (KO) mice of two genes (IMPA1 or Slc5a3), each encoding for a protein related to inositol metabolism, were studied in comparison with lithium-treated mice. Since we previously demonstrated that these KO mice exhibit a lithium-like neurochemical and behavioral phenotype, here we searched for pathways that may mediate lithium's/the KO effects. We performed a DNA-microarray study searching for pathways affected both by chronic lithium treatment and by the KO of each of the genes. The data were analyzed using three different bioinformatics approaches. We found upregulation of mitochondria-related genes in frontal cortex of lithium-treated, IMPA1 and Slc5a3 KO mice. Three out of seven genes differentially expressed in all three models, Cox5a, Ndufs7, and Ndufab, all members of the mitochondrial electron transfer chain, have previously been associated with bipolar disorder and/or lithium treatment. Upregulation of the expression of these genes was verified by real-time PCR. To further support the link between mitochondrial function and lithium's effect on behavior, we determined the capacity of chronic low-dose rotenone, a mitochondrial respiratory chain complex I inhibitor, to alter lithium-induced behavior as measured by the forced-swim and the amphetamine-induced hyperlocomotion paradigms. Rontenone treatment counteracted lithium's effect on behavior, supporting the proposition suggested by the bioinformatics analysis for a mitochondrial function involvement in behavioral effects of lithium mediated by inositol metabolism alterations.The results provide support for the notion that mitochondrial dysfunction is linked to bipolar disorder and can be ameliorated by lithium. The phenotypic similarities between lithium-treated wild-type mice and the two KO models suggest that lithium may affect behavior by altering inositol metabolism.

Gene expression analysis reveals functional pathways of glatiramer acetate activation

BACKGROUND

Glatiramer acetate (GA, Copaxone®), a mixture of polymers comprising four amino acids, is approved for treatment of relapsing-remitting multiple sclerosis and clinically isolated syndrome. GA mediates its activity by induction of GA-specific T cells that shift the T cell balance from a dominant proinflammatory phenotype (Th1/Th17) to an anti-inflammatory phenotype (Th2/Treg).

OBJECTIVE

To characterize the functional pathways by which GA acts on immune cells, the authors conducted gene expression profiling using glatiramoid-stimulated splenocytes.

METHODS

Mice were immunized with GA and harvested splenocytes were reactivated ex vivo with GA or a purported generic GA. Gene expression profiles and functional pathways were evaluated in reactivated splenocytes.

RESULTS

Overall, 1,474 genes were significantly upregulated or downregulated by GA. The main functional pathways induced by GA were: increased proliferation and activation of immune cells including T and B lymphocytes, stimulation of antigen presenting cells and differentiation of effector T lymphocytes. T-helper cell differentiation was the most significant canonical pathway associated with gene transcripts altered by GA. These expression patterns were not observed when splenocytes were activated with generic GA.

CONCLUSION

GA-induced functional pathways coincide with known mechanisms of GA activity in MS patients and further support the unique therapeutic effect of this drug.

Plant response to stress meets dedifferentiation

Plant response to various stress conditions often results in expression of common genes, known as stress-responsive/inducible genes. Accumulating data point to a common, yet elusive process underlying the response of plant cells to stress. Evidence derived from transcriptome profiling of shoot apical meristem stem cells, dedifferentiating protoplast cells as well as from senescing cells lends support to a model in which a common response of cells to certain biotic and abiotic stresses converges on cellular dedifferentiation whereby cells first acquire a stem cell-like state before assuming a new fate.

A lectin array-based methodology for the analysis of protein glycosylation

Glycosylation is the most versatile and one of the most abundant protein modifications. It has a structural role as well as diverse functional roles in many specific biological functions, including cancer development, viral and bacterial infections, and autoimmunity. The diverse roles of glycosylation in biological processes are rapidly growing areas of research, however, Glycobiology research is limited by the lack of a technology for rapid analysis of glycan composition of glycoproteins. Currently used methods for glycoanalysis are complex, typically requiring high levels of expertise and days to provide answers, and are not readily available to all researcher. We have developed a lectin array-based method, Qproteome GlycoArray kits, for rapid analysis of glycosylation profiles of glycoproteins. Glycoanalysis is performed on intact glycoproteins, requiring only 4-6 h for most analysis types. The method, demonstrated in this manuscript by several examples, is based on binding of an intact glycoprotein to the arrayed lectins, resulting in a characteristic fingerprint that is highly sensitive to changes in the protein's glycan composition. The large number of lectins, each with its specific recognition pattern, ensures high sensitivity to changes in the glycosylation pattern. A set of proprietary algorithms automatically interpret the fingerprint signals to provide a comprehensive glycan profile output.

DNA topoisomerase I in the mouse central nervous system: Age and sex dependence

Topoisomerase I (topo I) is a nuclear enzyme responsible for the topological state of DNA and therefore participates in most DNA transactions, particularly in transcription. Topo I, a ubiquitous enzyme, was identified and characterized in various cell types and tissues; however, the characterization of topo I in the intact central nervous system was not performed. Here we investigated, for the first time, the activity, level, and distribution pattern of topo I in the various selected brain regions in the mouse. In the visual cortex, cerebellum, and striatum the activity of topo I was 3-4-fold higher compared to that found in the hippocampus and hypothalamus. Immunohistochemical and immunofluorescence analyses revealed specific distribution patterns of topo I protein in neurons of each of the areas examined. The highest topo I levels were observed in inhibitory neurons. In addition to the expected nuclear localization of this protein, some neurons exhibited significant cytoplasmic content as well. The activity and level of topo I is age- and gender-dependent. It increases from birth to maturity and decreases, more significantly in males, with senescence. These results point to a possible importance and involvement of topo I activity and regulation in various brain functions.