Search

The CRISPR/Cas genome editing approach in non-model organisms poses challenges that remain to be resolved. Here, we demonstrated a generalized roadmap for a de novo genome annotation approach applied to the non-model organism Macrobrachium rosenbergii . We also addressed the typical genome editing challenges arising from genetic variations, such as a high frequency of single nucleotide polymorphisms, differences in sex chromosomes, and repetitive sequences that can lead to off-target events. For the genome editing of M. rosenbergii , our laboratory recently adapted the CRISPR/Cas genome editing approach to embryos and the embryonic primary cell culture. In this continuation study, an annotation pipeline was trained to predict the gene models by leveraging the available genomic, transcriptomic, and proteomic data, and enabling accurate gene prediction and guide design for knock-outs. A next-generation sequencing analysis demonstrated a high frequency of genetic variations in genes on both autosomal and sex chromosomes, which have been shown to affect the accuracy of editing analyses. To enable future applications based on the CRISPR/Cas tool in non-model organisms, we also verified the reliability of editing efficiency and tracked off-target frequencies. Despite the lack of comprehensive information on non-model organisms, this study provides an example of the feasibility of selecting and editing specific genes with a high degree of certainty.

Pathogenic variants underlying Mendelian diseases often disrupt the normal physiology of a few tissues and organs. However, variant effect prediction tools that aim to identify pathogenic variants are typically oblivious to tissue contexts. Here we report a machine-learning framework, denoted "Tissue Risk Assessment of Causality by Expression for variants" (TRACEvar, https://netbio.bgu.ac.il/TRACEvar/ ), that offers two advancements. First, TRACEvar predicts pathogenic variants that disrupt the normal physiology of specific tissues. This was achieved by creating 14 tissue-specific models that were trained on over 14,000 variants and combined 84 attributes of genetic variants with 495 attributes derived from tissue omics. TRACEvar outperformed 10 well-established and tissue-oblivious variant effect prediction tools. Second, the resulting models are interpretable, thereby illuminating variants' mode of action. Application of TRACEvar to variants of 52 rare-disease patients highlighted pathogenicity mechanisms and relevant disease processes. Lastly, the interpretation of all tissue models revealed that top-ranking determinants of pathogenicity included attributes of disease-affected tissues, particularly cellular process activities. Collectively, these results show that tissue contexts and interpretable machine-learning models can greatly enhance the etiology of rare diseases., (© 2024. The Author(s).)

Deconvolution algorithms mostly rely on single-cell RNA-sequencing (scRNA-seq) data applied onto bulk RNA-sequencing (bulk RNA-seq) to estimate tissues' cell-type composition, with performance accuracy validated on deposited databases. Adipose tissues' cellular composition is highly variable, and adipocytes can only be captured by single-nucleus RNA-sequencing (snRNA-seq). Here we report the development of sNucConv, a Scaden deep-learning-based deconvolution tool, trained using 5 hSAT and 7 hVAT snRNA-seq-based data corrected by (i) snRNA-seq/bulk RNA-seq highly correlated genes and (ii) individual cell-type regression models. Applying sNucConv on our bulk RNA-seq data resulted in cell-type proportion estimation of 15 and 13 cell types, with accuracy of R = 0.93 (range: 0.76-0.97) and R = 0.95 (range: 0.92-0.98) for hVAT and hSAT, respectively. This performance level was further validated on an independent set of samples (5 hSAT; 5 hVAT). The resulting model was depot specific, reflecting depot differences in gene expression patterns. Jointly, sNucConv provides proof-of-concept for producing validated deconvolution models for tissues un-amenable to scRNA-seq., Competing Interests: The authors declare no competing interests., (© 2024 The Authors.)

Background: Organisms from many distinct evolutionary lineages acquired the capacity to enter a dormant state in response to environmental conditions incompatible with maintaining normal life activities. Most studied organisms exhibit seasonal or annual episodes of dormancy, but numerous less studied organisms enter long-term dormancy, lasting decades or even centuries. Intriguingly, many planktonic animals produce encased embryos known as resting eggs or cysts that, like plant seeds, may remain dormant for decades. Herein, we studied a rotifer Brachionus plicatilis as a model planktonic species that forms encased dormant embryos via sexual reproduction and non-dormant embryos via asexual reproduction and raised the following questions: Which genes are expressed at which time points during embryogenesis? How do temporal transcript abundance profiles differ between the two types of embryos? When does the cell cycle arrest? How do dormant embryos manage energy?, Results: As the molecular developmental kinetics of encased embryos remain unknown, we employed single embryo RNA sequencing (CEL-seq) of samples collected during dormant and non-dormant embryogenesis. We identified comprehensive and temporal transcript abundance patterns of genes and their associated enriched functional pathways. Striking differences were uncovered between dormant and non-dormant embryos. In early development, the cell cycle-associated pathways were enriched in both embryo types but terminated with fewer nuclei in dormant embryos. As development progressed, the gene transcript abundance profiles became increasingly divergent between dormant and non-dormant embryos. Organogenesis was suspended in dormant embryos, concomitant with low transcript abundance of homeobox genes, and was replaced with an ATP-poor preparatory phase characterized by very high transcript abundance of genes encoding for hallmark dormancy proteins (e.g., LEA proteins, sHSP, and anti-ROS proteins, also found in plant seeds) and proteins involved in dormancy exit. Surprisingly, this period appeared analogous to the late maturation phase of plant seeds., Conclusions: The study highlights novel divergent temporal transcript abundance patterns between dormant and non-dormant embryos. Remarkably, several convergent functional solutions appear during the development of resting eggs and plant seeds, suggesting a similar preparatory phase for long-term dormancy. This study accentuated the broad novel molecular features of long-term dormancy in encased animal embryos that behave like "animal seeds"., (© 2024. The Author(s).)

Cell death is physiologically induced by specific mediators. However, our power to trigger the process in selected cells is quite limited. The protandric shrimp Hippolyte inermis offers a possible answer. Here, we analyse a de novo transcriptome of shrimp post-larvae fed on diatoms. The sex ratio of diatom-fed shrimps versus shrimps fed on control diets was dramatically altered, demonstrating the disruption of the androgenic gland, and their transcriptome revealed key modifications in gene expression. A wide transcriptomic analysis, validated by real-time qPCR, revealed that ferroptosis represents the primary factor to re-shape the body of this invertebrate, followed by further apoptotic events, and our findings open biotechnological perspectives for controlling the destiny of selected tissues. Ferroptosis was detected here for the first time in a crustacean. In addition, this is the first demonstration of a noticeable effect prompted by an ingested food, deeply impacting the gene networks of a young metazoan, definitely determining its future physiology and sexual differentiation.

Background: Alzheimer's disease (AD) exhibits mitochondrial dysfunctions associated with dysregulated metabolism, brain inflammation, synaptic loss, and neuronal cell death. As a key protein serving as the mitochondrial gatekeeper, the voltage-dependent anion channel-1 (VDAC1) that controls metabolism and Ca 2+ homeostasis is positioned at a convergence point for various cell survival and death signals. Here, we targeted VDAC1 with VBIT-4, a newly developed inhibitor of VDAC1 that prevents its pro-apoptotic activity, and mitochondria dysfunction., Methods: To address the multiple pathways involved in AD, neuronal cultures and a 5 × FAD mouse model of AD were treated with VBIT-4. We addressed multiple topics related to the disease and its molecular mechanisms using immunoblotting, immunofluorescence, q-RT-PCR, 3-D structural analysis and several behavioral tests., Results: In neuronal cultures, amyloid-beta (Aβ)-induced VDAC1 and p53 overexpression and apoptotic cell death were prevented by VBIT-4. Using an AD-like 5 × FAD mouse model, we showed that VDAC1 was overexpressed in neurons surrounding Aβ plaques, but not in astrocytes and microglia, and this was associated with neuronal cell death. VBIT-4 prevented the associated pathophysiological changes including neuronal cell death, neuroinflammation, and neuro-metabolic dysfunctions. VBIT-4 also switched astrocytes and microglia from being pro-inflammatory/neurotoxic to neuroprotective phenotype. Moreover, VBIT-4 prevented cognitive decline in the 5 × FAD mice as evaluated using several behavioral assessments of cognitive function. Interestingly, VBIT-4 protected against AD pathology, with no significant change in phosphorylated Tau and only a slight decrease in Aβ-plaque load., Conclusions: The study suggests that mitochondrial dysfunction with its gatekeeper VDAC1 is a promising target for AD therapeutic intervention, and VBIT-4 is a promising drug candidate for AD treatment., (© 2022. The Author(s).)

Plant adaptation to a desert environment and its endemic heat stress is poorly understood at the molecular level. The naturally heat-tolerant Brassicaceae species Anastatica hierochuntica is an ideal extremophyte model to identify genetic adaptations that have evolved to allow plants to tolerate heat stress and thrive in deserts. We generated an A. hierochuntica reference transcriptome and identified extremophyte adaptations by comparing Arabidopsis thaliana and A. hierochuntica transcriptome responses to heat, and detecting positively selected genes in A. hierochuntica. The two species exhibit similar transcriptome adjustment in response to heat and the A. hierochuntica transcriptome does not exist in a constitutive heat 'stress-ready' state. Furthermore, the A. hierochuntica global transcriptome as well as heat-responsive orthologs, display a lower basal and higher heat-induced expression than in A. thaliana. Genes positively selected in multiple extremophytes are associated with stomatal opening, nutrient acquisition, and UV-B induced DNA repair while those unique to A. hierochuntica are consistent with its photoperiod-insensitive, early-flowering phenotype. We suggest that evolution of a flexible transcriptome confers the ability to quickly react to extreme diurnal temperature fluctuations characteristic of a desert environment while positive selection of genes involved in stress tolerance and early flowering could facilitate an opportunistic desert lifestyle., (© 2022 The Authors. New Phytologist © 2022 New Phytologist Foundation.)

Objective: Up-regulated expression of transcription-factor E2F1 in human visceral adipose tissue (VAT) characterizes a dysmetabolic obesity sub-phenotype. An E2F1-miRNA network has been described in multiple cancers. Here we investigated whether elevated VAT-E2F1 in obesity is associated with VAT-miRNA alterations similar to, or distinct from, those described in cancer. Furthermore, we assessed if E2F1-associated miRNA changes may contribute to the link between high- VAT-E2F1 and a dysmetabolic obesity phenotype. Methods: We assembled a cohort of patients with obesity and high-VAT-E2F1, matched by age, sex, ±BMI to patients with low-VAT-E2F1, with and without obesity (8 patients/groupX3 groups). We performed Nanostring©-based miRNA profiling of VAT samples from all 24 patients. Candidate E2F1-related miRNAs were validated by qPCR in an independent cohort of patients with extreme obesity, with or without type-2-diabetes (T2DM) (n = 20). Bioinformatic tools and manipulation of E2F1 expression in cells were used to establish the plausibility of the functional VAT-E2F1-miRNA network in obesity. Results: Among n = 798 identified miRNAs, 17 were differentially expressed in relation to E2F1 and not to obesity itself. No evidence for the cancer-related E2F1-miRNA network was identified in human VAT in obesity. In HEK293-cells, overexpression/downregulation of E2F1 correspondingly altered the expression of miRNA-206 and miRNA-210-5p, two miRNAs with reported metabolic functions consistent with those of E2F1. In VAT from both cohorts, the expression of both miRNA-206 and 210-5p intercorrelated, and correlated with the expression of E2F1. In cohort 1 we did not detect significant associations with biochemical parameters. In cohort 2 of patients with extreme obesity, all those with high VAT-E2F1 showed a diabetes-complicated obesity phenotype and higher expression of miRNA-206 and miRNA-210-5p, which also correlated with fasting glucose levels (both miRNAs) and fasting insulin (miRNA-210-5p). Conclusions: Whilst the previously described cancer-related E2F1-miRNA network does not appear to operate in VAT in obesity, miRNAs-206 and 210-5p may link high-E2F1 expression in VAT with diabetes-complicated extreme obesity phenotype.

Mitochondrial SMAC/Diablo induces apoptosis by binding the inhibitor of apoptosis proteins (IAPs), thereby activating caspases and, subsequently, apoptosis. Previously, we found that despite its pro-apoptotic activity, SMAC/Diablo is overexpressed in cancer, and demonstrated that in cancer it possesses new essential and non-apoptotic functions that are associated with regulating phospholipid synthesis including modulating mitochondrial phosphatidylserine decarboxylase activity. Here, we demonstrate additional functions for SMAC/Diablo associated with inflammation and immunity. CRISPR/Cas9 SMAC/Diablo-depleted A549 lung cancer cells displayed inhibited cell proliferation and migration. Proteomics analysis of these cells revealed altered expression of proteins associated with lipids synthesis and signaling, vesicular transport and trafficking, metabolism, epigenetics, the extracellular matrix, cell signaling, and neutrophil-mediated immunity. SMAC-KO A549 cell-showed inhibited tumor growth and proliferation and activated apoptosis. The small SMAC-depleted "tumor" showed a morphology of alveoli-like structures, reversed epithelial-mesenchymal transition, and altered tumor microenvironment. The SMAC-lacking tumor showed reduced expression of inflammation-related proteins such as NF-kB and TNF-α, and of the PD-L1, associated with immune system suppression. These results suggest that SMAC is involved in multiple processes that are essential for tumor growth and progression. Thus, targeting SMAC's non-canonical function is a potential strategy to treat cancer., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Pandey, Shteinfer-Kuzmine, Chalifa-Caspi and Shoshan-Barmatz.)

Studies in rodent models suggest that calls emitted by isolated pups serve as an early behavioral manifestation of communication deficits and autistic like behavior. Previous studies in our labs showed that gestational exposure to the pesticide chlorpyrifos (CPF) and the Mthfr-knock-out mice are associated with impaired social preference and restricted or repetitive behavior. To extend these studies, we examine how pup communication via ultrasonic vocalizations is altered in these ASD models. We implemented an unsupervised hierarchical clustering method based on the spectral properties of the syllables in order to exploit syllable classification to homogeneous categories while avoiding over-categorization. Comparative exploration of the spectral and temporal aspects of syllables emitted by pups in two ASD models point to the following: (1) Most clusters showed a significant effect of the ASD factor on the start and end frequencies and bandwidth and (2) The highest percent change due to the ASD factor was on the bandwidth and duration. In addition, we found sex differences in the spectral and temporal properties of the calls in both control groups as well as an interaction between sex and the gene/environment factor. Considering the basal differences in the characteristics of syllables emitted by pups of the C57Bl/6 and Balb/c strains used as a background in the two models, we suggest that the above spectral-temporal parameters start frequency, bandwidth, and duration are the most sensitive USV features that may represent developmental changes in ASD models., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Shekel, Giladi, Raykin, Weiner, Chalifa-Caspi, Lederman, Kofman and Golan.)

In chronic obesity, activated adipose tissue proinflammatory cascades are tightly linked to metabolic dysfunction. Yet, close temporal analyses of the responses to obesogenic environment such as high-fat feeding (HFF) in susceptible mouse strains question the causal relationship between inflammation and metabolic dysfunction, and/or raises the possibility that certain inflammatory cascades play adaptive/homeostatic, rather than pathogenic roles. Here, we hypothesized that CTRP6, a C1QTNF family member, may constitute an early responder to acute nutritional changes in adipose tissue, with potential physiological roles. Both 3-days high-fat feeding (3dHFF) and acute obesity reversal [2-wk switch to low-fat diet after 8-wk HFF (8wHFF)] already induced marked changes in whole body fuel utilization. Although adipose tissue expression of classical proinflammatory cytokines ( Tnf-α , Ccl2 , and Il1b ) exhibited no, or only minor, change, C1qtnf6 uniquely increased, and decreased, in response to 3dHFF and acute obesity reversal, respectively. CTRP6 knockout (KO) mouse embryonic fibroblasts (MEFs) exhibited increased adipogenic gene expression ( Pparg , Fabp4 , and Adipoq ) and markedly reduced inflammatory genes ( Tnf-α , Ccl2 , and Il6 ) compared with wild-type MEFs, and recombinant CTRP6 induced the opposite gene expression signature, as assessed by RNA sequencing. Consistently, 3dHFF of CTRP6-KO mice induced a greater whole body and adipose tissue weight gain compared with wild-type littermates. Collectively, we propose CTRP6 as a gene that rapidly responds to acute changes in caloric intake, acting in acute overnutrition to induce a "physiological inflammatory response" that limits adipose tissue expansion. NEW & NOTEWORTHY CTRP6 (C1qTNF6), a member of adiponectin gene family, regulates inflammation and metabolism in established obesity. Here, short-term high-fat feeding in mice is shown to increase adipose tissue expression of CTRP6 before changes in the expression of classical inflammatory genes occur. Conversely, CTRP6 expression in adipose tissue decreases early in the course of obesity reversal. Gain- and loss-of-function models suggest CTRP6 as a positive regulator of inflammatory cascades, and a negative regulator of adipogenesis and adipose tissue expansion.

The tumor microenvironment (TME) plays an important role in cell growth, proliferation, migration, immunity, malignant transformation, and apoptosis. Thus, better insight into tumor-host interactions is required. Most of these processes involve the metabolic reprogramming of cells. Here, we focused on this reprogramming in cancerous cells and its effect on the TME. A major limitation in the study of tumor-host interactions is the difficulty in separating cancerous from non-cancerous signaling pathways within a tumor. Our strategy involved specifically silencing the expression of VDAC1 in the mitochondria of human-derived A549 lung cancer xenografts in mice, but not in the mouse-derived cells of the TME. Next-generation sequencing (NGS) analysis allows distinguishing the human or mouse origin of genes, thus enabling the separation of the bidirectional cross-talk between the TME and malignant cells. We demonstrate that depleting VDAC1 in cancer cells led to metabolic reprogramming, tumor regression, and the disruption of tumor-host interactions. This was reflected in the altered expression of a battery of genes associated with TME, including those involved in extracellular matrix organization and structure, matrix-related peptidases, angiogenesis, intercellular interacting proteins, integrins, and growth factors associated with stromal activities. We show that metabolic rewiring upon mitochondrial VDAC1 silencing in cancer cells affected several components of the TME, such as structural protein matrix metalloproteinases and Lox, and elicited a stromal response resembling the reaction to a foreign body in wound healing. As tumor progression requires a cooperative interplay between the host and cancer cells, and the ECM is intensively remodeled during cancer progression, VDAC1 depletion induced metabolic reprogramming that targeted both tumor cells and resulted in the alteration of the whole spectrum of TME-related genes, affecting the reciprocal feedback between ECM molecules, host cells, and cancer cells. Thus, VDAC1 depletion using si-VDAC1 represents therapeutic potential, inhibiting cancer cell proliferation and also inducing the modulation of TME components, which influences cancer progression, migration, and invasion.

The transcriptional coactivator BRD4 has a fundamental role in transcription regulation and thus became a promising epigenetic therapeutic candidate to target diverse pathologies. However, the regulation of BRD4 by posttranslational modifications has been largely unexplored. Here, we show that BRD4 is methylated on chromatin at lysine-99 by the protein lysine methyltransferase SETD6. BRD4 methylation negatively regulates the expression of genes that are involved in translation and inhibits total mRNA translation in cells. Mechanistically, we provide evidence that supports a model where BRD4 methylation by SETD6 does not have a direct role in the association with acetylated histone H4 at chromatin. However, this methylation specifically determines the recruitment of the transcription factor E2F1 to selected target genes that are involved in mRNA translation. Together, our findings reveal a previously unknown molecular mechanism for BRD4 methylation-dependent gene-specific targeting, which may serve as a new direction for the development of therapeutic applications., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).)

RNA-Seq is nowadays an indispensable approach for comparative transcriptome profiling in model and nonmodel organisms. Analyzing RNA-Seq data from nonmodel organisms poses unique challenges, due to unavailability of a high-quality genome reference and to relative sparsity of tools for downstream functional analyses. In this chapter, we provide an overview of the analysis steps in RNA-Seq projects of nonmodel organisms, while elaborating on aspects that are unique to this analysis. These will include (1) strategic decisions that have to be made in advance, regarding sequencing technology and reference to use; (2) how to search for available draft genomes, and, if necessary, how to improve their gene prediction and annotation; (3) how to clean raw reads before de novo assembly; (4) how to separate the reads in RNA-Seq projects of symbiont organisms; (5) how to design and carry out a de novo transcriptome assembly that will be comprehensive and reliable; (6) how to assess transcriptome quality; (7) when and how to reduce redundancy in the transcriptome; (8) techniques and considerations in transcriptome functional annotation; (9) quantitating transcript abundance in the face of high transcriptome redundancy; and, most importantly, (10) how to achieve functional enrichment testing using available tools which either support a large range of species or enable a universal, non-species-specific analysis.Throughout the chapter, we will refer to a variety of useful software tools. For the initial analysis steps involving high-volume data, these will include Linux-based programs. For the later steps, we will describe both Linux and R packages for advanced users, as well as many user-friendly tools for nonprogrammers. Finally, we will present a full workflow for RNA-Seq analysis of nonmodel organisms using the NeatSeq-Flow platform, which can be used locally through a user-friendly interface.

One fundamental character common to pancrustaceans (Crustacea and Hexapoda) is a mineralized rigid exoskeleton whose principal organic components are chitin and proteins. In contrast to traditional research in the field that has been devoted to the structural and physicochemical aspects of biomineralization, the present study explores transcriptomic aspects of biomineralization as a first step towards adding a complementary molecular layer to this field. The rigidity of the exoskeleton in pancrustaceans dictates essential molt cycles enabling morphological changes and growth. Thus, formation and mineralization of the exoskeleton are concomitant to the timeline of the molt cycle. Skeletal proteinaceous toolkit elements have been discovered in previous studies using innovative molt-related binary gene expression patterns derived from transcriptomic libraries representing the major stages comprising the molt cycle of the decapod crustacean Cherax quadricarinatus. Here, we revisited some prominent exoskeleton-related structural proteins encoding and, using the above molt-related binary pattern methodology, enlarged the transcriptomic database of C. quadricarinatus. The latter was done by establishing a new transcriptomic library of the cuticle forming epithelium and molar tooth at four different molt stages (i.e., inter-molt, early pre-molt, late pre-molt and post-molt) and incorporating it to a previous transcriptome derived from the gastroliths and mandible. The wider multigenic approach facilitated by the newly expanded transcriptomic database not only revisited single genes of the molecular toolkit, but also provided both scattered and specific information that broaden the overview of proteins and gene clusters which are involved in the construction and biomineralization of the exoskeleton in decapod crustaceans., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

Elevated expression of E2F1 in adipocyte fraction of human visceral adipose tissue (hVAT) associates with a poor cardiometabolic profile. We hypothesized that beyond directly activating autophagy and MAP3K5 (ASK)-MAP kinase signaling, E2F1 governs a distinct transcriptome that contributes to adipose tissue and metabolic dysfunction in obesity. We performed RNA sequencing of hVAT samples from age-, sex-, and BMI-matched patients, all obese, whose visceral E2F1 protein expression was either high (E2F1 high ) or low (E2F1 low ). Tumor necrosis factor superfamily (TNFSF) members, including TRAIL ( TNFSF10 ), TL1A ( TNFSF15 ), and their receptors, were enriched in E2F1 high While TRAIL was equally expressed in adipocytes and stromal vascular fraction (SVF), TL1A was mainly expressed in SVF, and TRAIL-induced TL1A was attributed to CD4 + and CD8 + subclasses of hVAT T cells. In human adipocytes, TL1A enhanced basal and impaired insulin-inhibitable lipolysis and altered adipokine secretion, and in human macrophages it induced foam cell biogenesis and M1 polarization. Two independent human cohorts confirmed associations between TL1A and TRAIL expression in hVAT and higher leptin and IL6 serum concentrations, diabetes status, and hVAT-macrophage lipid content. Jointly, we propose an intra-adipose tissue E2F1-associated TNFSF paracrine loop engaging lymphocytes, macrophages, and adipocytes, ultimately contributing to adipose tissue dysfunction in obesity., (© 2020 by the American Diabetes Association.)

In the Australian redclaw crayfish, Cherax quadricarinatus (WZ/ZZ system), intersexuals, although exhibiting both male and female gonopores, are functional males bearing a female genotype (WZ males). Therefore, the occurrence of the unusual homogametic WW females in nature is plausible. We developed W/Z genomic sex markers and used them to investigate the genotypic structure of experimental and native C. quadricarinatus populations in Australia. We discovered, for the first time, the natural occurrence of WW females in crustacean populations. By modeling population dynamics, we found that intersexuals contribute to the growth rate of crayfish populations in the short term. Given the vastly fragmented C. quadricarinatus habitat, which is characterized by drought-flood cycles, we speculate that intersexuals contribute to the fitness of this species since they lead to occasional increment in the population growth rate which potentially supports crayfish population restoration and establishment under extinction threats or colonization events., Competing Interests: A patent regarding sex-specific genomic markers in the Australian redclaw crayfish is pending (International application number: PCT/IL2018/051046, International publication number: WO/2019/058371)., (© 2020 The Authors.)

Motivation: Differential network analysis, designed to highlight network changes between conditions, is an important paradigm in network biology. However, differential network analysis methods have been typically designed to compare between two conditions and were rarely applied to multiple protein interaction networks (interactomes). Importantly, large-scale benchmarks for their evaluation have been lacking., Results: Here, we present a framework for assessing the ability of differential network analysis of multiple human tissue interactomes to highlight tissue-selective processes and disorders. For this, we created a benchmark of 6499 curated tissue-specific Gene Ontology biological processes. We applied five methods, including four differential network analysis methods, to construct weighted interactomes for 34 tissues. Rigorous assessment of this benchmark revealed that differential analysis methods perform well in revealing tissue-selective processes (AUCs of 0.82-0.9). Next, we applied differential network analysis to illuminate the genes underlying tissue-selective hereditary disorders. For this, we curated a dataset of 1305 tissue-specific hereditary disorders and their manifesting tissues. Focusing on subnetworks containing the top 1% differential interactions in disease-relevant tissue interactomes revealed significant enrichment for disorder-causing genes in 18.6% of the cases, with a significantly high success rate for blood, nerve, muscle and heart diseases., Summary: Altogether, we offer a framework that includes expansive manually curated datasets of tissue-selective processes and disorders to be used as benchmarks or to illuminate tissue-selective processes and genes. Our results demonstrate that differential analysis of multiple human tissue interactomes is a powerful tool for highlighting processes and genes with tissue-selective functionality and clinical impact., Availability and Implementation: Datasets are available as part of the Supplementary data., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2020. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Carcinogenesis is a complicated process that involves the deregulation of epigenetics, resulting in cellular transformational events, such as proliferation, differentiation, and metastasis. Most chromatin-modifying enzymes utilize metabolites as co-factors or substrates and thus are directly dependent on such metabolites as acetyl-coenzyme A, S-adenosylmethionine, and NAD+. Here, we show that using specific siRNA to deplete a tumor of VDAC1 not only led to reprograming of the cancer cell metabolism but also altered several epigenetic-related enzymes and factors. VDAC1, in the outer mitochondrial membrane, controls metabolic cross-talk between the mitochondria and the rest of the cell, thus regulating the metabolic and energetic functions of mitochondria, and has been implicated in apoptotic-relevant events. We previously demonstrated that silencing VDAC1 expression in glioblastoma (GBM) U-87MG cell-derived tumors, resulted in reprogramed metabolism leading to inhibited tumor growth, angiogenesis, epithelial-mesenchymal transition and invasiveness, and elimination of cancer stem cells, while promoting the differentiation of residual tumor cells into neuronal-like cells. These VDAC1 depletion-mediated effects involved alterations in transcription factors regulating signaling pathways associated with cancer hallmarks. As the epigenome is sensitive to cellular metabolism, this study was designed to assess whether depleting VDAC1 affects the metabolism-epigenetics axis. Using DNA microarrays, q-PCR, and specific antibodies, we analyzed the effects of si-VDAC1 treatment of U-87MG-derived tumors on histone modifications and epigenetic-related enzyme expression levels, as well as the methylation and acetylation state, to uncover any alterations in epigenetic properties. Our results demonstrate that metabolic rewiring of GBM via VDAC1 depletion affects epigenetic modifications, and strongly support the presence of an interplay between metabolism and epigenetics.

Reprograming of the metabolism of cancer cells is an event recognized as a hallmark of the disease. The mitochondrial gatekeeper, voltage-dependent anion channel 1 (VDAC1), mediates transport of metabolites and ions in and out of mitochondria, and is involved in mitochondria-mediated apoptosis. Here, we compared the effects of reducing hVDAC1 expression in a glioblastoma xenograft using human-specific si-RNA (si-hVDAC1) for a short (19 days) and a long term (40 days). Tumors underwent reprograming, reflected in rewired metabolism, eradication of cancer stem cells (CSCs) and differentiation. Short- and long-term treatments of the tumors with si-hVDAC1 similarly reduced the expression of metabolism-related enzymes, and translocator protein (TSPO) and CSCs markers. In contrast, differentiation into cells expressing astrocyte or neuronal markers was noted only after a long period during which the tumor cells were hVDAC1-depleted. This suggests that tumor cell differentiation is a prolonged process that precedes metabolic reprograming and the "disappearance" of CSCs. Tumor proteomics analysis revealing global changes in the expression levels of proteins associated with signaling, synthesis and degradation of proteins, DNA structure and replication and epigenetic changes, all of which were highly altered after a long period of si-hVDAC1 tumor treatment. The depletion of hVDAC1 greatly reduced the levels of the multifunctional translocator protein TSPO, which is overexpressed in both the mitochondria and the nucleus of the tumor. The results thus show that VDAC1 depletion-mediated cancer cell metabolic reprograming involves a chain of events occurring in a sequential manner leading to a reversal of the unique properties of the tumor, indicative of the interplay between metabolism and oncogenic signaling networks., Competing Interests: The authors declare no conflict of interest.

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., (Copyright © 2019 Mulla, Hajaj, Elyagon, Mor, Gillis, Murninkas, Klapper-Goldstein, Plaschkes, Chalifa-Caspi, Etzion and Etzion.)

The cultivation of monosex populations is common in animal husbandry. However, preselecting the desired gender remains a major biotechnological and ethical challenge. To achieve an efficient biotechnology for all-female aquaculture in the economically important prawn (Macrobrachium rosenbergii), we achieved - for the first time - WW males using androgenic gland cells transplantation which caused full sex-reversal of WW females to functional males. Crossing the WW males with WW females yielded all-female progeny lacking the Z chromosome. We now have the ability to manipulate - by non-genomic means - all possible genotype combinations (ZZ, WZ and WW) to retain either male or female phenotypes and hence to produce monosex populations of either gender. This calls for a study of the genomic basis underlying this striking sexual plasticity, questioning the content of the W and Z chromosomes. Here, we report on the sequencing of a high-quality genome exhibiting distinguishable paternal and maternal sequences. This assembly covers ~ 87.5% of the genome and yielded a remarkable N50 value of ~ 20 × 10 6  bp. Genomic sex markers were used to initiate the identification and validation of parts of the W and Z chromosomes for the first time in arthropods.

Age-associated changes in CD4 T-cell functionality have been linked to chronic inflammation and decreased immunity. However, a detailed characterization of CD4 T cell phenotypes that could explain these dysregulated functional properties is lacking. We used single-cell RNA sequencing and multidimensional protein analyses to profile thousands of CD4 T cells obtained from young and old mice. We found that the landscape of CD4 T cell subsets differs markedly between young and old mice, such that three cell subsets-exhausted, cytotoxic, and activated regulatory T cells (aT regs )-appear rarely in young mice but gradually accumulate with age. Most unexpected were the extreme pro- and anti-inflammatory phenotypes of cytotoxic CD4 T cells and aT regs , respectively. These findings provide a comprehensive view of the dynamic reorganization of the CD4 T cell milieu with age and illuminate dominant subsets associated with chronic inflammation and immunity decline, suggesting new therapeutic avenues for age-related diseases.

Numerous aquatic invertebrates survive harsh environments by displaying dormancy as encysted embryos. This study aimed at determining whether metabolomics could provide molecular insight to explain the "dormancy syndrome" by highlighting functional pathways and metabolites, hence offering a novel comprehensive molecular view of dormancy. We compared the metabolome of morphologically distinct dormant encysted embryos (resting eggs) and non-dormant embryos (amictic eggs) of a rotifer (Brachionus plicatilis). Metabolome profiling revealed ~5,000 features, 1,079 of which were annotated. Most of the features were represented at significantly higher levels in non-dormant than dormant embryos. A large number of features was assigned to putative functional pathways indicating novel differences between dormant and non-dormant states. These include features associated with glycolysis, the TCA and urea cycles, amino acid, purine and pyrimidine metabolism. Interestingly, ATP, nucleobases, cyclic nucleotides, thymidine and uracil, were not detected in dormant resting eggs, suggesting an impairment of response to environmental and internal cues, cessation of DNA synthesis, transcription and plausibly translation in the dormant embryos. The levels of trehalose or its analogues, with a role in survival under desiccation conditions, were higher in resting eggs. In conclusion, the current study highlights metabolomics as a major analytical tool to functionally compare dormancy across species.

Motivation: Identification of conserved syntenic blocks across microbial genomes is important for several problems in comparative genomics such as gene annotation, study of genome organization and evolution and prediction of gene interactions. Current tools for syntenic block discovery do not scale up to the large quantity of prokaryotic genomes available today., Results: We present a novel methodology for the discovery, ranking and taxonomic distribution analysis of colinear syntenic blocks (CSBs)-groups of genes that are consistently located close to each other, in the same order, across a wide range of taxa. We present an efficient algorithm that identifies CSBs in large genomic datasets. The algorithm is implemented and incorporated in a novel tool with a graphical user interface, denoted CSBFinder, that ranks the discovered CSBs according to a probabilistic score and clusters them to families according to their gene content similarity. We apply CSBFinder to data mine 1487 prokaryotic genomes including chromosomes and plasmids. For post-processing analysis, we generate heatmaps for visualizing the distribution of CSB family members across various taxa. We exemplify the utility of CSBFinder in operon prediction, in deciphering unknown gene function and in taxonomic analysis of colinear syntenic blocks., Availability and Implementation: CSBFinder software and code are publicly available at https://github.com/dinasv/CSBFinder., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2018. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Under natural conditions, plants are regularly exposed to combinations of stress factors. A common example is the conjunction between nitrogen (N) deficiency and excess light. The combined effect of stress factors is often ignored in studies using controlled conditions, possibly resulting in misleading conclusions. To address this issue, the present study examined the physiological behavior of Arabidopsis thaliana under the effect of varying nitrogen levels and light intensities. The joint influence of low N and excess light had an adverse effect on plant growth, chlorophyll and anthocyanin concentrations, photochemical capacity and the abundance of proteins involved in carbon assimilation and antioxidative metabolism. In contrast, no adverse physiological responses were observed for plants under either nitrogen limitation or high light (HL) intensity conditions (i.e. single stress). The underlying mechanisms for the increased growth in conditions of HL and sufficient nitrogen were a combination of chlorophyll accumulation and an increased number of proteins involved in C3 carbon assimilation, amino acids biosynthesis and chloroplast development. In contrast, combined stress conditions shifts plants from growth to survival by displaying anthocyanin accumulation and an increased number of proteins involved in catabolism of lipids and amino acids as energy substrates. Ultimately switching plants development from growth to survival. Our results suggest that an assessment of the physiological response to the combined effect of multiple stresses cannot be directly extrapolated from the physiological response to a single stress. Specifically, the synergistic interaction between N deficiency and saturating light in Arabidopsis plants could not have been modeled via only one of the stress factors., (© 2018 Scandinavian Plant Physiology Society.)

Trigger factor (TF) has a known cytoplasmic function as a chaperone. In a previous study we showed that pneumococcal TF is also cell-wall localized and this finding combined with the immunogenic characteristic of TF, has led us to determine the vaccine potential of TF and decipher its involvement in pneumococcal pathogenesis. Bioinformatic analysis revealed that TF is conserved among pneumococci and has no human homologue. Immunization of mice with recombinant (r)TF elicited a protective immune response against a pneumococcal challenge, suggesting that TF contributes to pneumococcal pathogenesis. Indeed, rTF and an anti-rTF antiserum inhibited bacterial adhesion to human lung derived epithelial cells, indicating that TF contributes to the bacterial adhesion to the host. Moreover, bacteria lacking TF demonstrated reduced adhesion, in vitro, to lung-derived epithelial cells, neural cells and glial cells. The reduced adhesion could be restored by chromosomal complementation. Furthermore, bacteria lacking TF demonstrated significantly reduced virulence in a mouse model. Taken together, the ability of rTF to elicit a protective immune response, involvement of TF in bacterial adhesion, conservation of the protein among pneumococcal strains and the lack of human homologue, all suggest that rTF can be considered as a future candidate vaccine with a much broader coverage as compared to the currently available pneumococcal vaccines.

Asteriscus graveolens ( A. graveolens ) plants contain among other metabolites, sesquiterpene lactone asteriscunolide isomers (AS). The crude extract and its fractions affected the viability of mouse BS-24-1 lymphoma cells (BS-24-1 cells) with an IC50 of 3 μg/mL. The fraction was cytotoxic to cancer cells but not to non-cancerous cells (human induced pluripotent stem cells); its activity was accompanied by a concentration- and time-dependent appearance of apoptosis as determined by DNA fragmentation and caspase-3 activity. High levels of Reactive Oxygen Species (ROS) were rapidly observed (less than 1 min) after addition of the fraction followed by an increase in caspase-3 activity three hours later. Comparison of RNA-seq transcriptome profiles from pre-and post-treatment of BS-24-1 cells with crude extract of A. graveolens yielded a list of 2293 genes whose expression was significantly affected. This gene set included genes encoding proteins involved in cell cycle arrest, protection against ROS, and activation of the tumor suppressor P53 pathway, supporting the biochemical findings on ROS species-dependent apoptosis induced by A. graveolens fraction. Interestingly, several of the pathways and genes affected by A. graveolens extract are expressed following treatment of human cancer cells with chemotherapy drugs. We suggest, that A. graveolens extracts maybe further developed into selective chemotherapy.

We formalize a new problem variant in gene-block discovery, denoted Reference-Anchored Gene Blocks (RAGB), given a query sequence Q of length n, representing the gene array of a DNA element, a window size bound d on the length of a substring of interest in Q, and a set of target gene sequences [Formula: see text]. Our objective is to identify gene blocks in [Formula: see text] that are centered in a subset q of co-localized genes from Q, and contain genomes from [Formula: see text] in which the corresponding orthologs of the genes from q are also co-localized. We cast RAGB as a variant of a (colored) biclique problem in bipartite graphs, and analyze its parameterized complexity, as well as the parameterized complexity of other related problems. We give an [Formula: see text] time algorithm for the uncolored variant of our biclique problem, where m is the number of areas of interest that are parsed from the target sequences, and n and d are as defined earlier. Our algorithm can be adapted to compute all maximal bicliques in the graph within the same time complexity, and to handle edge weights with a slight [Formula: see text] increase to its time complexity. For the colored version of the problem, our algorithm has a time complexity of [Formula: see text]. We implement the algorithm and exemplify its application to the data mining of proteobacterial gene blocks that are centered in predicted proteobacterial genomic islands, leading to the identification of putatively mobilized clusters of virulence, pathogenicity, and resistance genes.

Pneumococcal flavin reductase (FlaR) is known to be cell-wall associated and possess age dependent antigenicity in children. This study aimed at characterizing FlaR and elucidating its involvement in pneumococcal physiology and virulence. Bioinformatic analysis of FlaR sequence identified three-conserved cysteine residues, suggesting a transition metal-binding capacity. Recombinant FlaR (rFlaR) bound Fe 2+ and exhibited FAD-dependent NADP-reductase activity, which increased in the presence of cysteine or excess Fe 2+ and inhibited by divalent-chelating agents. flaR mutant was highly susceptible to H 2 O 2 compared to its wild type (WT) and complemented strains, suggesting a role for FlaR in pneumococcal oxidative stress resistance. Additionally, flaR mutant demonstrated significantly decreased mice mortality following intraperitoneal infection. Interestingly, lack of FlaR did not affect the extent of phagocytosis by primary mouse peritoneal macrophages but reduced adhesion to A549 cells compared to the WT and complemented strains. Noteworthy are the findings that immunization with rFlaR elicited protection in mice against intraperitoneal lethal challenge and anti-FlaR antisera neutralized bacterial virulence. Taken together, FlaR's roles in pneumococcal physiology and virulence, combined with its lack of significant homology to human proteins, point towards rFlaR as a vaccine candidate.

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 plicatilis A 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., (© 2017 by The American Society for Biochemistry and Molecular Biology, Inc.)

It is commonly assumed that the dead, hardened floral bracts of the dispersal unit of grasses have been evolved to protect seeds from predation and / or assist in fruit/caryopsis dispersal. While these structures have important agronomical and economical implications, their adaptive value has not been fully explored. We investigated the hypothesis that the maternally derived hardened floral bracts have been evolved not just as a means for caryopsis protection and dispersal, but also as storage for substances that might affect seed germination and seedling vigor. Dead glumes as well as lemmas and paleas of wild emmer wheat (Triticum turgidum var dicoccoides) were found to store and release upon hydration active hydrolases including nucleases and chitinases. High nuclease activity was released upon hydration from glumes derived from wild strains of wheat including Triticum urartu and wild emmer wheat, while very low nuclease activity was detected in glumes derived from domesticated, free-threshing wheat cultivars (e.g., durum wheat). Germination from the intact dispersal unit of wild emmer wheat was delayed, but post germination growth was better than those of separated caryopses. Most notable was a significant increase in lateral root production on seedlings germinated from the intact dispersal unit. Proteome analysis of wild emmer wheat glumes revealed many proteins stored and released upon hydration including S1-type nucleases, peptidases, antifungal hydrolases such as chitinases and β-1,3-glucanase as well as pectin acetylesterase, a protein involved in cell wall degradation and remodeling. Also, reactive oxygen species (ROS)-detoxifying enzymes such as superoxide dismutase and ascorbate peroxidase were overrepresented in dead glumes of wild emmer wheat. Thus our study highlighted previously unknown features of the dispersal unit in wild wheat in which the dead, hardened floral bracts enclosing the caryopsis store active hydrolases and nutritional elements and probably growth promoting substances that facilitate seed longevity and germination and increase seedling vigor.

Aim: The transmission dynamics of Acinetobacter baumannii in endemic settings, and the relation between microbial properties and patients' clinical outcomes, are yet obscure and hampered by insufficient metadata., Methods & Results: Of 20 consecutive patients with A. baumannii bloodstream infection that were thoroughly analyzed at a single center, at least one transmission opportunity was evident for 85% of patients. This implies that patient-to-patient transmission is the major mode of A. baumannii acquisitions in health facilities. Moreover, all patients who died immediately (<24 h of admission) were infected with a single clone (ST457; relative risk = 1.6; p = 0.05)., Conclusion: This preliminary analysis should prompt further investigation by mapping genomic virulence determinants among A. baumannii ST457 lineage compared with other strains.

Background: The ability of regeneration is essential for the homeostasis of all animals as it allows the repair and renewal of tissues and body parts upon normal turnover or injury. The extent of this ability varies greatly in different animals with the sea anemone Nematostella vectensis, a basal cnidarian model animal, displaying remarkable whole-body regeneration competence., Results: In order to study this process in Nematostella we performed an RNA-Seq screen wherein we analyzed and compared the transcriptional response to bisection in the wound-proximal body parts undergoing oral (head) or aboral (tail) regeneration at several time points up to the initial restoration of the basic body shape. The transcriptional profiles of regeneration responsive genes were analyzed so as to define the temporal pattern of differential gene expression associated with the tissue-specific oral and aboral regeneration. The identified genes were characterized according to their GO (gene ontology) assignations revealing groups that were enriched in the regeneration process with particular attention to their affiliation to the major developmental signaling pathways. While some of the genes and gene groups thus analyzed were previously known to be active in regeneration, we have also revealed novel and surprising candidate genes such as cilia-associated genes that likely participate in this important developmental program., Conclusions: This work highlighted the main groups of genes which showed polarization upon regeneration, notably the proteinases, multiple transcription factors and the Wnt pathway genes that were highly represented, all displaying an intricate temporal balance between the two sides. In addition, the evolutionary comparison performed between regeneration in different animal model systems may reveal the basic mechanisms playing a role in this fascinating process.

In Streptococcus pneumonia, phosphoenolpyruvate protein phosphotransferase (PtsA) is an intracellular protein of the monosaccharide phosphotransferase systems. Biochemical and immunostaining methods were applied to show that PtsA also localizes to the bacterial cell-wall. Thus, it was suspected that PtsA has functions other than its main cytoplasmic enzymatic role. Indeed, recombinant PtsA and anti-rPtsA antiserum were shown to inhibit adhesion of S. pneumoniae to cultured human lung adenocarcinoma A549 cells. Screening of a combinatorial peptide library expressed in a filamentous phage with rPtsA identified epitopes that were capable of inhibiting S. pneumoniae adhesion to A549 cells. The insert peptides in the phages were sequenced, and homologous sequences were found in human BMPER, multimerin1, protocadherin19, integrinβ4, epsin1 and collagen type VIIα1 proteins, all of which can be found in A549 cells except the latter. Six peptides, synthesized according to the homologous sequences in the human proteins, specifically bound rPtsA in the micromolar range and significantly inhibited pneumococcal adhesion in vitro to lung- and tracheal-derived cell lines. In addition, the tested peptides inhibited lung colonization after intranasal inoculation of mice with S. pneumoniae. Immunization with rPtsA protected the mice against a sublethal intranasal and a lethal intravenous pneumococcal challenge. In addition, mouse anti rPtsA antiserum reduced bacterial virulence in the intravenous inoculation mouse model. These findings showed that the surface-localized PtsA functions as an adhesin, PtsA binding peptides derived from its putative target molecules can be considered for future development of therapeutics, and rPtsA should be regarded as a candidate for vaccine development.

The emerging microbial source tracking (MST) methodologies aim to identify fecal contamination originating from domestic and wild animals, and from humans. Avian MST is especially challenging, primarily because the Aves class includes both domesticated and wild species with highly diverse habitats and dietary characteristics. The quest for specific fecal bacterial MST markers can be difficult with respect to attaining sufficient assay sensitivity and specificity. The present study utilizes high throughput sequencing (HTS) to screen bacterial 16S rRNA genes from fecal samples collected from both domestic and wild avian species. Operational taxonomic unit (OTU) analysis was then performed, from which sequences were retained for downstream quantitative polymerase chain reaction (qPCR) marker development. Identification of unique avian host DNA sequences, absent in non-avian hosts, was then carried out using a dedicated database of bacterial 16S rRNA gene taken from the Ribosomal Database Project. Six qPCR assays were developed targeting the 16S rRNA gene of Lactobacillus, Gallibacterium, Firmicutes, Fusobacteriaceae, and other bacteria. Two assays (Av4143 and Av163) identified most of the avian fecal samples and demonstrated sensitivity values of 91 and 70%, respectively. The Av43 assay only identified droppings from battery hens and poultry, whereas each of the other three assays (Av24, Av13, and Av216) identified waterfowl species with lower sensitivities values. The development of an MST assay-panel, which includes both domestic and wild avian species, expands the currently known MST analysis capabilities for decoding fecal contamination.

As challenges to food security increase, the demand for lead genes for improving crop production is growing. However, genetic screens of plant mutants typically yield very low frequencies of desired phenotypes. Here, we present a powerful computational approach for selecting candidate genes for screening insertion mutants. We combined ranking of Arabidopsis thaliana regulatory genes according to their expression in response to multiple abiotic stresses (Multiple Stress [MST] score), with stress-responsive RNA co-expression network analysis to select candidate multiple stress regulatory (MSTR) genes. Screening of 62 T-DNA insertion mutants defective in candidate MSTR genes, for abiotic stress germination phenotypes yielded a remarkable hit rate of up to 62%; this gene discovery rate is 48-fold greater than that of other large-scale insertional mutant screens. Moreover, the MST score of these genes could be used to prioritize them for screening. To evaluate the contribution of the co-expression analysis, we screened 64 additional mutant lines of MST-scored genes that did not appear in the RNA co-expression network. The screening of these MST-scored genes yielded a gene discovery rate of 36%, which is much higher than that of classic mutant screens but not as high as when picking candidate genes from the co-expression network. The MSTR co-expression network that we created, AraSTressRegNet is publicly available at http://netbio.bgu.ac.il/arnet. This systems biology-based screening approach combining gene ranking and network analysis could be generally applicable to enhancing identification of genes regulating additional processes in plants and other organisms provided that suitable transcriptome data are available., (© 2014 Society for Experimental Biology, Association of Applied Biologists and John Wiley & Sons Ltd.)

In crustaceans, like all arthropods, growth is accompanied by a molting cycle. This cycle comprises major physiological events in which mineralized chitinous structures are built and degraded. These events are in turn governed by genes whose patterns of expression are presumably linked to the molting cycle. To study these genes we performed next generation sequencing and constructed a molt-related transcriptomic library from two exoskeletal-forming tissues of the crayfish Cherax quadricarinatus, namely the gastrolith and the mandible cuticle-forming epithelium. To simplify the study of such a complex process as molting, a novel approach, binary patterning of gene expression, was employed. This approach revealed that key genes involved in the synthesis and breakdown of chitin exhibit a molt-related pattern in the gastrolith-forming epithelium. On the other hand, the same genes in the mandible cuticle-forming epithelium showed a molt-independent pattern of expression. Genes related to the metabolism of glucosamine-6-phosphate, a chitin precursor synthesized from simple sugars, showed a molt-related pattern of expression in both tissues. The binary patterning approach unfolds typical patterns of gene expression during the molt cycle of a crustacean. The use of such a simplifying integrative tool for assessing gene patterning seems appropriate for the study of complex biological processes.