Loss of EIF4G2 mediates aggressiveness in distinct human endometrial cancer subpopulations with poor survival outcome in patients

The non-canonical translation initiation factor EIF4G2 plays essential roles in cellular stress responses via translation of selective mRNA cohorts. Currently there is limited and conflicting information regarding its involvement in cancer development and progression. Here we assessed its role in endometrial cancer (EC), in a cohort of 280 EC patients across different types, grades, and stages, and found that low EIF4G2 expression highly correlated with poor overall- and recurrence-free survival in Grade 2 EC patients, monitored over a period of up to 12 years. To establish a causative connection between low EIF4G2 expression and cancer progression, we stably knocked-down EIF4G2 in two human EC cell lines in parallel. EIF4G2 depletion resulted in increased resistance to conventional therapies and increased the prevalence of molecular markers for aggressive cell subsets, altering their transcriptional and proteomic landscapes. Prominent among the proteins with decreased abundance were Kinesin-1 motor proteins, KIF5B and KLC1, 2, 3. Multiplexed imaging of the EC patient tumor cohort showed a correlation between decreased expression of the kinesin proteins, and poor survival in patients with tumors of certain grades and stages. These findings reveal potential novel biomarkers for Grade 2 EC with ramifications for patient stratification and therapeutic interventions.

PRDM16 co-operates with LHX2 to shape the human brain

PRDM16 is a dynamic transcriptional regulator of various stem cell niches, including adipocytic, hematopoietic, cardiac progenitors, and neural stem cells. PRDM16 has been suggested to contribute to 1p36 deletion syndrome, one of the most prevalent subtelomeric microdeletion syndromes. We report a patient with a de novo nonsense mutation in the PRDM16 coding sequence, accompanied by lissencephaly and microcephaly features. Human stem cells were genetically modified to mimic this mutation, generating cortical organoids that exhibited altered cell cycle dynamics. RNA sequencing of cortical organoids at day 32 unveiled changes in cell adhesion and WNT-signaling pathways. ChIP-seq of PRDM16 identified binding sites in postmortem human fetal cortex, indicating the conservation of PRDM16 binding to developmental genes in mice and humans, potentially at enhancer sites. A shared motif between PRDM16 and LHX2 was identified and further examined through comparison with LHX2 ChIP-seq data from mice. These results suggested a collaborative partnership between PRDM16 and LHX2 in regulating a common set of genes and pathways in cortical radial glia cells, possibly via their synergistic involvement in cortical development.

Metabolic inputs in the probiotic bacterium Lacticaseibacillus rhamnosus contribute to cell-wall remodeling and increased fitness

Lacticaseibacillus rhamnosus GG (LGG) is a Gram-positive beneficial bacterium that resides in the human intestinal tract and belongs to the family of lactic acid bacteria (LAB). This bacterium is a widely used probiotic and was suggested to provide numerous benefits for human health. However, as in most LAB strains, the molecular mechanisms that mediate the competitiveness of probiotics under different diets remain unknown. Fermentation is a fundamental process in LAB, allowing the oxidation of simple carbohydrates (e.g., glucose, mannose) for energy production under oxygen limitation, as in the human gut. Our results indicate that fermentation reshapes the metabolome, volatilome, and proteome architecture of LGG. Furthermore, fermentation alters cell envelope remodeling and peptidoglycan biosynthesis, which leads to altered cell wall thickness, aggregation properties, and cell wall composition. In addition, fermentable sugars induced the secretion of known and novel metabolites and proteins targeting the enteric pathogens Enterococcus faecalis and Salmonella enterica Serovar Typhimurium. Overall, our results link simple carbohydrates with cell wall remodeling, aggregation to host tissues, and biofilm formation in probiotic strains and connect them with the production of broad-spectrum antimicrobial effectors.

Sedentary Behavior Impacts on the Epigenome and Transcriptome: Lessons from Muscle Inactivation in Drosophila Larvae

The biological mechanisms linking sedentary lifestyles and metabolic derangements are incompletely understood. In this study, temporal muscle inactivation in Drosophila larvae carrying a temperature-sensitive mutation in the shibire (shi1) gene was induced to mimic sedentary behavior during early life and study its transcriptional outcome. Our findings indicated a significant change in the epigenetic profile, as well as the genomic profile, of RNA Pol II binding in the inactive muscles relative to control, within a relatively short time period. Whole-genome analysis of RNA-Pol II binding to DNA by muscle-specific targeted DamID (TaDa) protocol revealed that muscle inactivity altered Pol II binding in 121 out of 2010 genes (6%), with a three-fold enrichment of genes coding for lncRNAs. The suppressed protein-coding genes included genes associated with longevity, DNA repair, muscle function, and ubiquitin-dependent proteostasis. Moreover, inducing muscle inactivation exerted a multi-level impact upon chromatin modifications, triggering an altered epigenetic balance of active versus inactive marks. The downregulated genes in the inactive muscles included genes essential for muscle structure and function, carbohydrate metabolism, longevity, and others. Given the multiple analogous genes in Drosophila for many human genes, extrapolating our findings to humans may hold promise for establishing a molecular link between sedentary behavior and metabolic diseases.

PRDM16 co-operates with LHX2 to shape the human brain

PRDM16 is a dynamic transcriptional regulator of various stem cell niches, including adipocytic, hematopoietic, cardiac progenitors, and neural stem cells. PRDM16 has been suggested to contribute to 1p36 deletion syndrome, one of the most prevalent subtelomeric microdeletion syndromes. We report a patient with a de novo nonsense mutation in the PRDM16 coding sequence, accompanied by lissencephaly and microcephaly features. Human stem cells were genetically modified to mimic this mutation, generating cortical organoids that exhibited altered cell cycle dynamics. RNA sequencing of cortical organoids at day 32 unveiled changes in cell adhesion and WNT-signaling pathways. ChIP-seq of PRDM16 identified binding sites in postmortem human fetal cortex, indicating the conservation of PRDM16 binding to developmental genes in mice and humans, potentially at enhancer sites. A shared motif between PRDM16 and LHX2 was identified and further examined through comparison with LHX2 ChIP-seq data from mice. These results suggested a collaborative partnership between PRDM16 and LHX2 in regulating a common set of genes and pathways in cortical radial glia cells, possibly via their synergistic involvement in cortical development.

LIS1 RNA-binding orchestrates the mechanosensitive properties of embryonic stem cells in AGO2-dependent and independent ways

Lissencephaly-1 (LIS1) is associated with neurodevelopmental diseases and is known to regulate the molecular motor cytoplasmic dynein activity. Here we show that LIS1 is essential for the viability of mouse embryonic stem cells (mESCs), and it governs the physical properties of these cells. LIS1 dosage substantially affects gene expression, and we uncovered an unexpected interaction of LIS1 with RNA and RNA-binding proteins, most prominently the Argonaute complex. We demonstrate that LIS1 overexpression partially rescued the extracellular matrix (ECM) expression and mechanosensitive genes conferring stiffness to Argonaute null mESCs. Collectively, our data transforms the current perspective on the roles of LIS1 in post-transcriptional regulation underlying development and mechanosensitive processes.

Global analysis of contact-dependent human-to-mouse intercellular mRNA and lncRNA transfer in cell culture

Full-length mRNAs transfer between adjacent mammalian cells via direct cell-to-cell connections called tunneling nanotubes (TNTs). However, the extent of mRNA transfer at the transcriptome-wide level (the 'transferome') is unknown. Here, we analyzed the transferome in an in vitro human-mouse cell co-culture model using RNA-sequencing. We found that mRNA transfer is non-selective, prevalent across the human transcriptome, and that the amount of transfer to mouse embryonic fibroblasts (MEFs) strongly correlates with the endogenous level of gene expression in donor human breast cancer cells. Typically, <1% of endogenous mRNAs undergo transfer. Non-selective, expression-dependent RNA transfer was further validated using synthetic reporters. RNA transfer appears contact-dependent via TNTs, as exemplified for several mRNAs. Notably, significant differential changes in the native MEF transcriptome were observed in response to co-culture, including the upregulation of multiple cancer and cancer-associated fibroblast-related genes and pathways. Together, these results lead us to suggest that TNT-mediated RNA transfer could be a phenomenon of physiological importance under both normal and pathogenic conditions.

The LINC Complex Inhibits Excessive Chromatin Repression

The Linker of Nucleoskeleton and Cytoskeleton (LINC) complex transduces nuclear mechanical inputs suggested to control chromatin organization and gene expression; however, the underlying mechanism is currently unclear. We show here that the LINC complex is needed to minimize chromatin repression in muscle tissue, where the nuclei are exposed to significant mechanical inputs during muscle contraction. To this end, the genomic binding profiles of Polycomb, Heterochromatin Protein1 (HP1a) repressors, and of RNA-Pol II were studied in Drosophila larval muscles lacking functional LINC complex. A significant increase in the binding of Polycomb and parallel reduction of RNA-Pol-II binding to a set of muscle genes was observed. Consistently, enhanced tri-methylated H3K9 and H3K27 repressive modifications and reduced chromatin activation by H3K9 acetylation were found. Furthermore, larger tri-methylated H3K27me3 repressive clusters, and chromatin redistribution from the nuclear periphery towards nuclear center, were detected in live LINC mutant larval muscles. Computer simulation indicated that the observed dissociation of the chromatin from the nuclear envelope promotes growth of tri-methylated H3K27 repressive clusters. Thus, we suggest that by promoting chromatin-nuclear envelope binding, the LINC complex restricts the size of repressive H3K27 tri-methylated clusters, thereby limiting the binding of Polycomb transcription repressor, directing robust transcription in muscle fibers.

rG4detector, a novel RNA G-quadruplex predictor, uncovers their impact on stress granule formation

RNA G-quadruplexes (rG4s) are RNA secondary structures, which are formed by guanine-rich sequences and have important cellular functions. Existing computational tools for rG4 prediction rely on specific sequence features and/or were trained on small datasets, without considering rG4 stability information, and are therefore sub-optimal. Here, we developed rG4detector, a convolutional neural network to identify potential rG4s in transcriptomics data. rG4detector outperforms existing methods in both predicting rG4 stability and in detecting rG4-forming sequences. To demonstrate the biological-relevance of rG4detector, we employed it to study RNAs that are bound by the RNA-binding protein G3BP1. G3BP1 is central to the induction of stress granules (SGs), which are cytoplasmic biomolecular condensates that form in response to a variety of cellular stresses. Unexpectedly, rG4detector revealed a dynamic enrichment of rG4s bound by G3BP1 in response to cellular stress. In addition, we experimentally characterized G3BP1 cross-talk with rG4s, demonstrating that G3BP1 is a bona fide rG4-binding protein and that endogenous rG4s are enriched within SGs. Furthermore, we found that reduced rG4 availability impairs SG formation. Hence, we conclude that rG4s play a direct role in SG biology via their interactions with RNA-binding proteins and that rG4detector is a novel useful tool for rG4 transcriptomics data analyses.

DAP5 drives translation of specific mRNA targets with upstream ORFs in human embryonic stem cells

Death associated protein 5 (DAP5/eIF4G2/NAT1) is a member of the eIF4G translation initiation factors that has been shown to mediate noncanonical and/or cap-independent translation. It is essential for embryonic development and for differentiation of embryonic stem cells (ESCs), specifically its ability to drive translation of specific target mRNAs. In order to expand the repertoire of DAP5 target mRNAs, we compared ribosome profiles in control and DAP5 knockdown (KD) human ESCs (hESCs) to identify mRNAs with decreased ribosomal occupancy upon DAP5 silencing. A cohort of 68 genes showed decreased translation efficiency in DAP5 KD cells. Mass spectrometry confirmed decreased protein abundance of a significant portion of these targets. Among these was KMT2D, a histone methylase previously shown to be essential for ESC differentiation and embryonic development. We found that nearly half of the cohort of DAP5 target mRNAs displaying reduced translation efficiency of their main coding sequences upon DAP5 KD contained upstream open reading frames (uORFs) that are actively translated independently of DAP5. This is consistent with previously suggested mechanisms by which DAP5 mediates leaky scanning through uORFs and/or reinitiation at the main coding sequence. Crosslinking protein-RNA immunoprecipitation experiments indicated that a significant subset of DAP5 mRNA targets bound DAP5, indicating that direct binding between DAP5 protein and its target mRNAs is a frequent but not absolute requirement for DAP5-dependent translation of the main coding sequence. Thus, we have extended DAP5's function in translation of specific mRNAs in hESCs by a mechanism allowing translation of the main coding sequence following upstream translation of short ORFs.

Heterogeneous nuclear ribonucleoprotein U (HNRNPU) safeguards the developing mouse cortex

HNRNPU encodes the heterogeneous nuclear ribonucleoprotein U, which participates in RNA splicing and chromatin organization. Microdeletions in the 1q44 locus encompassing HNRNPU and other genes and point mutations in HNRNPU cause brain disorders, including early-onset seizures and severe intellectual disability. We aimed to understand HNRNPU’s roles in the developing brain. Our work revealed that HNRNPU loss of function leads to rapid cell death of both postmitotic neurons and neural progenitors, with an apparent higher sensitivity of the latter. Further, expression and alternative splicing of multiple genes involved in cell survival, cell motility, and synapse formation are affected following Hnrnpu’s conditional truncation. Finally, we identified pharmaceutical and genetic agents that can partially reverse the loss of cortical structures in Hnrnpu mutated embryonic brains, ameliorate radial neuronal migration defects and rescue cultured neural progenitors’ cell death.

HNRNPU is an RNA splicing protein associated with brain disorders such as early onset seizures. Here they show that HNRNPU functions to maintain neural progenitors and their progeny by regulating splicing of key neuronal genes.

The roles of intracellular and extracellular calcium in Bacillus subtilis biofilms

In nature, bacteria reside in biofilms– multicellular differentiated communities held together by an extracellular matrix. This work identified a novel subpopulation—mineral-forming cells—that is essential for biofilm formation in Bacillus subtilis biofilms. This subpopulation contains an intracellular calcium-accumulating niche, in which the formation of a calcium carbonate mineral is initiated. As the biofilm colony develops, this mineral grows in a controlled manner, forming a functional macrostructure that serves the entire community. Consistently, biofilm development is prevented by the inhibition of calcium uptake. Our results provide a clear demonstration of the orchestrated production of calcite exoskeleton, critical to morphogenesis in simple prokaryotes.

•

The orchestrated formation of calcite scaffolds supports the morphogenesis of microbial biofilms

•

A novel subpopulation—mineral-forming cells—is essential for biofilm formation

•

This subpopulation contains an intracellular calcium-accumulating niche, supporting the formation of calcium carbonate

•

Intracellular calcium homeostasis and calcium export are associated with a functional biofilm macrostructure

Whole-genome sequencing reveals that variants in the Interleukin 18 Receptor Accessory Protein 3'UTR protect against ALS

The noncoding genome is substantially larger than the protein-coding genome but has been largely unexplored by genetic association studies. Here, we performed region-based rare variant association analysis of >25,000 variants in untranslated regions of 6,139 amyotrophic lateral sclerosis (ALS) whole genomes and the whole genomes of 70,403 non-ALS controls. We identified interleukin-18 receptor accessory protein (IL18RAP) 3' untranslated region (3'UTR) variants as significantly enriched in non-ALS genomes and associated with a fivefold reduced risk of developing ALS, and this was replicated in an independent cohort. These variants in the IL18RAP 3'UTR reduce mRNA stability and the binding of double-stranded RNA (dsRNA)-binding proteins. Finally, the variants of the IL18RAP 3'UTR confer a survival advantage for motor neurons because they dampen neurotoxicity of human induced pluripotent stem cell (iPSC)-derived microglia bearing an ALS-associated expansion in C9orf72, and this depends on NF-κB signaling. This study reveals genetic variants that protect against ALS by reducing neuroinflammation and emphasizes the importance of noncoding genetic association studies.

Bi-fated tendon-to-bone attachment cells are regulated by shared enhancers and KLF transcription factors

The mechanical challenge of attaching elastic tendons to stiff bones is solved by the formation of a unique transitional tissue. Here, we show that murine tendon-to-bone attachment cells are bi-fated, activating a mixture of chondrocyte and tenocyte transcriptomes, under regulation of shared regulatory elements and Krüppel-like factors (KLFs) transcription factors. High-throughput bulk and single-cell RNA sequencing of humeral attachment cells revealed expression of hundreds of chondrogenic and tenogenic genes, which was validated by in situ hybridization and single-molecule ISH. ATAC sequencing showed that attachment cells share accessible intergenic chromatin areas with either tenocytes or chondrocytes. Epigenomic analysis revealed enhancer signatures for most of these regions. Transgenic mouse enhancer reporter assays verified the shared activity of some of these enhancers. Finally, integrative chromatin and motif analyses and transcriptomic data implicated KLFs as regulators of attachment cells. Indeed, blocking expression of both Klf2 and Klf4 in developing limb mesenchyme impaired their differentiation.

Spatiotemporal Proteomic Analysis of Stress Granule Disassembly Using APEX Reveals Regulation by SUMOylation and Links to ALS Pathogenesis

Stress granules (SGs) are cytoplasmic assemblies of proteins and non-translating mRNAs. Whereas much has been learned about SG formation, a major gap remains in understanding the compositional changes SGs undergo during normal disassembly and under disease conditions. Here, we address this gap by proteomic dissection of the SG temporal disassembly sequence using multi-bait APEX proximity proteomics. We discover 109 novel SG proteins and characterize distinct SG substructures. We reveal dozens of disassembly-engaged proteins (DEPs), some of which play functional roles in SG disassembly, including small ubiquitin-like modifier (SUMO) conjugating enzymes. We further demonstrate that SUMOylation regulates SG disassembly and SG formation. Parallel proteomics with amyotrophic lateral sclerosis (ALS)-associated C9ORF72 dipeptides uncovered attenuated DEP recruitment during SG disassembly and impaired SUMOylation. Accordingly, SUMO activity ameliorated C9ORF72-ALS-related neurodegeneration in Drosophila. By dissecting the SG spatiotemporal proteomic landscape, we provide an in-depth resource for future work on SG function and reveal basic and disease-relevant mechanisms of SG disassembly.

Human genetics and neuropathology suggest a link between miR-218 and amyotrophic lateral sclerosis pathophysiology

Motor neuron-specific microRNA-218 (miR-218) has recently received attention because of its roles in mouse development. However, miR-218 relevance to human motor neuron disease was not yet explored. Here, we demonstrate by neuropathology that miR-218 is abundant in healthy human motor neurons. However, in amyotrophic lateral sclerosis (ALS) motor neurons, miR-218 is down-regulated and its mRNA targets are reciprocally up-regulated (derepressed). We further identify the potassium channel Kv10.1 as a new miR-218 direct target that controls neuronal activity. In addition, we screened thousands of ALS genomes and identified six rare variants in the human miR-218-2 sequence. miR-218 gene variants fail to regulate neuron activity, suggesting the importance of this small endogenous RNA for neuronal robustness. The underlying mechanisms involve inhibition of miR-218 biogenesis and reduced processing by DICER. Therefore, miR-218 activity in motor neurons may be susceptible to failure in human ALS, suggesting that miR-218 may be a potential therapeutic target in motor neuron disease.

The extracellular matrix protein TasA is a developmental cue that maintains a motile subpopulation within Bacillus subtilis biofilms

In nature, bacteria form biofilms-differentiated multicellular communities attached to surfaces. Within these generally sessile biofilms, a subset of cells continues to express motility genes. We found that this subpopulation enabled Bacillus subtilis biofilms to expand on high-friction surfaces. The extracellular matrix (ECM) protein TasA was required for the expression of flagellar genes. In addition to its structural role as an adhesive fiber for cell attachment, TasA acted as a developmental signal stimulating a subset of biofilm cells to revert to a motile phenotype. Transcriptomic analysis revealed that TasA stimulated the expression of a specific subset of genes whose products promote motility and repress ECM production. Spontaneous suppressor mutations that restored motility in the absence of TasA revealed that activation of the biofilm-motility switch by the two-component system CssR/CssS antagonized the TasA-mediated reversion to motility in biofilm cells. Our results suggest that although mostly sessile, biofilms retain a degree of motility by actively maintaining a motile subpopulation.

A unified nomenclature for vertebrate olfactory receptors

BACKGROUND

Olfactory receptors (ORs) are G protein-coupled receptors with a crucial role in odor detection. A typical mammalian genome harbors ~ 1000 OR genes and pseudogenes; however, different gene duplication/deletion events have occurred in each species, resulting in complex orthology relationships. While the human OR nomenclature is widely accepted and based on phylogenetic classification into 18 families and further into subfamilies, for other mammals different and multiple nomenclature systems are currently in use, thus concealing important evolutionary and functional insights.

RESULTS

Here, we describe the Mutual Maximum Similarity (MMS) algorithm, a systematic classifier for assigning a human-centric nomenclature to any OR gene based on inter-species hierarchical pairwise similarities. MMS was applied to the OR repertoires of seven mammals and zebrafish. Altogether, we assigned symbols to 10,249 ORs. This nomenclature is supported by both phylogenetic and synteny analyses. The availability of a unified nomenclature provides a framework for diverse studies, where textual symbol comparison allows immediate identification of potential ortholog groups as well as species-specific expansions/deletions; for example, Or52e5 and Or52e5b represent a rat-specific duplication of OR52E5. Another example is the complete absence of OR subfamily OR6Z among primate OR symbols. In other mammals, OR6Z members are located in one genomic cluster, suggesting a large deletion in the great ape lineage. An additional 14 mammalian OR subfamilies are missing from the primate genomes. While in chimpanzee 87% of the symbols were identical to human symbols, this number decreased to ~ 50% in dog and cow and to ~ 30% in rodents, reflecting the adaptive changes of the OR gene superfamily across diverse ecological niches. Application of the proposed nomenclature to zebrafish revealed similarity to mammalian ORs that could not be detected from the current zebrafish olfactory receptor gene nomenclature.

CONCLUSIONS

We have consolidated a unified standard nomenclature system for the vertebrate OR superfamily. The new nomenclature system will be applied to cow, horse, dog and chimpanzee by the Vertebrate Gene Nomenclature Committee and its implementation is currently under consideration by other relevant species-specific nomenclature committees.

Interplay of LIS1 and MeCP2: Interactions and Implications With the Neurodevelopmental Disorders Lissencephaly and Rett Syndrome

LIS1 is the main causative gene for lissencephaly, while MeCP2 is the main causative gene for Rett syndrome, both of which are neurodevelopmental diseases. Here we report nuclear functions for LIS1 and identify previously unrecognized physical and genetic interactions between the products of these two genes in the cell nucleus, that has implications on MeCP2 organization, neuronal gene expression and mouse behavior. Reduced LIS1 levels affect the association of MeCP2 with chromatin. Transcriptome analysis of primary cortical neurons derived from wild type, Lis1±, MeCP2−/y, or double mutants mice revealed a large overlap in the differentially expressed (DE) genes between the various mutants. Overall, our findings provide insights on molecular mechanisms involved in the neurodevelopmental disorders lissencephaly and Rett syndrome caused by dysfunction of LIS1 and MeCP2, respectively.

A secretion-enhancing cis regulatory targeting element (SECReTE) involved in mRNA localization and protein synthesis

The localization of mRNAs encoding secreted/membrane proteins (mSMPs) to the endoplasmic reticulum (ER) likely facilitates the co-translational translocation of secreted proteins. However, studies have shown that mSMP recruitment to the ER in eukaryotes can occur in a manner that is independent of the ribosome, translational control, and the signal recognition particle, although the mechanism remains largely unknown. Here, we identify a cis-acting RNA sequence motif that enhances mSMP localization to the ER and appears to increase mRNA stability, and both the synthesis and secretion of secretome proteins. Termed SECReTE, for secretion-enhancing cisregulatory targeting element, this motif is enriched in mRNAs encoding secretome proteins translated on the ER in eukaryotes and on the inner membrane of prokaryotes. SECReTE consists of ≥10 nucleotide triplet repeats enriched with pyrimidine (C/U) every third base (i.e. NNY, where N = any nucleotide, Y = pyrimidine) and can be present in the untranslated as well as the coding regions of the mRNA. Synonymous mutations that elevate the SECReTE count in a given mRNA (e.g. SUC2, HSP150, and CCW12) lead to an increase in protein secretion in yeast, while a reduction in count led to less secretion and physiological defects. Moreover, the addition of SECReTE to the 3’UTR of an mRNA for an exogenously expressed protein (e.g. GFP) led to its increased secretion from yeast cells. Thus, SECReTE constitutes a novel RNA motif that facilitates ER-localized mRNA translation and protein secretion.

Proteins destined for secretion from the cell, including soluble secreted and membrane proteins (SMPs), are translocated into the endoplasmic reticulum (ER) either directly upon translation on the ER surface or post-translationally, as in the case of type II membrane proteins. Interestingly, several studies have demonstrated that mRNAs encoding SMPs (mSMPs) are also enriched on the ER, yet how they target this organelle is less clear. The signal recognition particle (SRP), which recognizes N-terminal hydrophobic signals of nascent polypeptides and targets them to an ER-localized receptor, was proposed to mediate the co-translational ER targeting of mRNA. However, more recent studies show that SRP inactivation, as well as the inhibition of translation, do not prevent targeting. Thus, how mSMPs reach the ER and whether the process is translation-independent remain open. Here we identify a cis-acting sequence element in mSMPs that appears to facilitate mRNA stability and localization to the ER and, more importantly, enhances protein secretion. This motif, entitled “SECReTE” (secretion-enhancing cisregulatory targeting element) is enriched in nearly all mSMPs in eukaryotes and its addition or removal from mRNAs results in either enhanced or reduced protein secretion, respectively. Thus, SECReTE is a RNA sequence motif that regulates protein translation and secretion.

An active β-lactamase is a part of an orchestrated cell wall stress resistance network of Bacillus subtilis and related rhizosphere species

A hallmark of the Gram-positive bacteria, such as the soil-dwelling bacterium Bacillus subtilis, is their cell wall. Here, we report that d-leucine and flavomycin, biofilm inhibitors targeting the cell wall, activate the β-lactamase PenP. This β-lactamase contributes to ampicillin resistance in B. subtilis under all conditions tested. In contrast, both Spo0A, a master regulator of nutritional stress, and the general cell wall stress response, differentially contribute to β-lactam resistance under different conditions. To test whether β-lactam resistance and β-lactamase genes are widespread in other Bacilli, we isolated Bacillus species from undisturbed soils, and found that their genomes can encode up to five β-lactamases with differentiated activity spectra. Surprisingly, the activity of environmental β-lactamases and PenP, as well as the general stress response, resulted in a similarly reduced lag phase of the culture in the presence of β-lactam antibiotics, with little or no impact on the logarithmic growth rate. The length of the lag phase may determine the outcome of the competition between β-lactams and β-lactamases producers. Overall, our work suggests that antibiotic resistance genes in B. subtilis and related species are ancient and widespread, and could be selected by interspecies competition in undisturbed soils.

The Disordered Landscape of the 20S Proteasome Substrates Reveals Tight Association with Phase Separated Granules

Proteasomal degradation is the main route of regulated proteostasis. The 20S proteasome is the core particle (CP) responsible for the catalytic activity of all proteasome complexes. Structural constraints mean that only unfolded, extended polypeptide chains may enter the catalytic core of the 20S proteasome. It has been previously shown that the 20S CP is active in degradation of certain intrinsically disordered proteins (IDP) lacking structural constrains. Here, a comprehensive analysis of the 20S CP substrates in vitro is conducted. It is revealed that the 20S CP substrates are highly disordered. However, not all the IDPs are 20S CP substrates. The group of the IDPs that are 20S CP substrates, termed 20S-IDPome are characterized by having significantly more protein binding partners, more posttranslational modification sites, and are highly enriched for RNA binding proteins. The vast majority of them are involved in splicing, mRNA processing, and translation. Remarkably, it is found that low complexity proteins with prion-like domain (PrLD), which interact with GR or PR di-peptide repeats, are the most preferential 20S CP substrates. The finding suggests roles of the 20S CP in gene transcription and formation of phase-separated granules.

Developmental activities of the complement pathway in migrating neurons

In recent years the notion that malfunctioning of the immune system may result in developmental brain diseases has emerged. However, the role of immune molecules in the developing brain has not been well explored. The complement pathway converges to cleave C3. Here we show that key proteins in the lectin arm of this pathway, MASP1, MASP2 and C3, are expressed in the developing cortex and that neuronal migration is impaired in knockout and knockdown mice. Molecular mimics of C3 cleavage products rescue the migration defects that have been seen following knockdown of C3 or Masp2. Pharmacological activation of the downstream receptors rescue Masp2 and C3 knockdown as well as C3 knockout. Therefore, we propose that the complement pathway is functionally important in migrating neurons of the developing cortex.

Characterization of proteome dynamics during growth in oleate reveals a new peroxisome-targeting receptor

To optimally perform the diversity of metabolic functions that occur within peroxisomes, cells must dynamically regulate peroxisome size, number and content in response to the cell state and the environment. Except for transcriptional regulation little is known about the mechanisms used to perform this complicated feat. Focusing on the yeast Saccharomyces cerevisiae, we used complementary high-content screens to follow changes in localization of most proteins during growth in oleate. We found extensive changes in cellular architecture and identified several proteins that colocalized with peroxisomes that had not previously been considered peroxisomal proteins. One of the newly identified peroxisomal proteins, Ymr018w, is a protein with an unknown function that is similar to the yeast and human peroxisomal targeting receptor Pex5. We demonstrate that Ymr018w is a new peroxisomal-targeting receptor that targets a subset of matrix proteins to peroxisomes. We, therefore, renamed Ymr018w, Pex9, and suggest that Pex9 is a condition-specific targeting receptor that enables the dynamic rewiring of peroxisomes in response to metabolic needs. Moreover, we suggest that Pex5-like receptors might also exist in vertebrates.

ORDB, HORDE, ODORactor and other on-line knowledge resources of olfactory receptor-odorant interactions

We present here an exploration of the evolution of three well-established, web-based resources dedicated to the dissemination of information related to olfactory receptors (ORs) and their functional ligands, odorants. These resources are: the Olfactory Receptor Database (ORDB), the Human Olfactory Data Explorer (HORDE) and ODORactor. ORDB is a repository of genomic and proteomic information related to ORs and other chemosensory receptors, such as taste and pheromone receptors. Three companion databases closely integrated with ORDB are OdorDB, ORModelDB and OdorMapDB; these resources are part of the SenseLab suite of databases (http://senselab.med.yale.edu). HORDE (http://genome.weizmann.ac.il/horde/) is a semi-automatically populated database of the OR repertoires of human and several mammals. ODORactor (http://mdl.shsmu.edu.cn/ODORactor/) provides information related to OR-odorant interactions from the perspective of the odorant. All three resources are connected to each other via web-links.

Database URL: http://senselab.med.yale.edu; http://genome.weizmann.ac.il/horde/; http://mdl.shsmu.edu.cn/ODORactor/

The human olfactory transcriptome

BACKGROUND

Olfaction is a versatile sensory mechanism for detecting thousands of volatile odorants. Although molecular basis of odorant signaling is relatively well understood considerable gaps remain in the complete charting of all relevant gene products. To address this challenge, we applied RNAseq to four well-characterized human olfactory epithelial samples and compared the results to novel and published mouse olfactory epithelium as well as 16 human control tissues.

RESULTS

We identified 194 non-olfactory receptor (OR) genes that are overexpressed in human olfactory tissues vs.

CONTROLS

The highest overexpression is seen for lipocalins and bactericidal/permeability-increasing (BPI)-fold proteins, which in other species include secreted odorant carriers. Mouse-human discordance in orthologous lipocalin expression suggests different mammalian evolutionary paths in this family. Of the overexpressed genes 36 have documented olfactory function while for 158 there is little or no previous such functional evidence. The latter group includes GPCRs, neuropeptides, solute carriers, transcription factors and biotransformation enzymes. Many of them may be indirectly implicated in sensory function, and ~70 % are over expressed also in mouse olfactory epithelium, corroborating their olfactory role. Nearly 90 % of the intact OR repertoire, and ~60 % of the OR pseudogenes are expressed in the olfactory epithelium, with the latter showing a 3-fold lower expression. ORs transcription levels show a 1000-fold inter-paralog variation, as well as significant inter-individual differences. We assembled 160 transcripts representing 100 intact OR genes. These include 1-4 short 5' non-coding exons with considerable alternative splicing and long last exons that contain the coding region and 3' untranslated region of highly variable length. Notably, we identified 10 ORs with an intact open reading frame but with seemingly non-functional transcripts, suggesting a yet unreported OR pseudogenization mechanism. Analysis of the OR upstream regions indicated an enrichment of the homeobox family transcription factor binding sites and a consensus localization of a specific transcription factor binding site subfamily (Olf/EBF).

CONCLUSIONS

We provide an overview of expression levels of ORs and auxiliary genes in human olfactory epithelium. This forms a transcriptomic view of the entire OR repertoire, and reveals a large number of over-expressed uncharacterized human non-receptor genes, providing a platform for future discovery.

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.

A role for TENM1 mutations in congenital general anosmia

Congenital general anosmia (CGA) is a neurological disorder entailing a complete innate inability to sense odors. While the mechanisms underlying vertebrate olfaction have been studied in detail, there are still gaps in our understanding of the molecular genetic basis of innate olfactory disorders. Applying whole-exome sequencing to a family multiply affected with CGA, we identified three members with a rare X-linked missense mutation in the TENM1 (teneurin 1) gene (ENST00000422452:c.C4829T). In Drosophila melanogaster, TENM1 functions in synaptic-partner-matching between axons of olfactory sensory neurons and target projection neurons and is involved in synapse organization in the olfactory system. We used CRISPR-Cas9 system to generate a Tenm1 disrupted mouse model. Tenm1(-/-) and point-mutated Tenm1(A) (/A) adult mice were shown to have an altered ability to locate a buried food pellet. Tenm1(A) (/A) mice also displayed an altered ability to sense aversive odors. Results of our study, that describes a new Tenm1 mouse, agree with the hypothesis that TENM1 has a role in olfaction. However, additional studies should be done in larger CGA cohorts, to provide statistical evidence that loss-of-function mutations in TENM1 can solely cause the disease in our and other CGA cases.

Identification of a Functional Risk Variant for Pemphigus Vulgaris in the ST18 Gene

Pemphigus vulgaris (PV) is a life-threatening autoimmune mucocutaneous blistering disease caused by disruption of intercellular adhesion due to auto-antibodies directed against epithelial components. Treatment is limited to immunosuppressive agents, which are associated with serious adverse effects. The propensity to develop the disease is in part genetically determined. We therefore reasoned that the delineation of PV genetic basis may point to novel therapeutic strategies. Using a genome-wide association approach, we recently found that genetic variants in the vicinity of the ST18 gene confer a significant risk for the disease. Here, using targeted deep sequencing, we identified a PV-associated variant residing within the ST18 promoter region (p<0.0002; odds ratio = 2.03). This variant was found to drive increased gene transcription in a p53/p63-dependent manner, which may explain the fact that ST18 is up-regulated in the skin of PV patients. We then discovered that when overexpressed, ST18 stimulates PV serum-induced secretion of key inflammatory molecules and contributes to PV serum-induced disruption of keratinocyte cell-cell adhesion, two processes previously implicated in the pathogenesis of PV. Thus, the present findings indicate that ST18 may play a direct role in PV and consequently represents a potential target for the treatment of this disease.

Genic insights from integrated human proteomics in GeneCards

GeneCards is a one-stop shop for searchable human gene annotations (http://www.genecards.org/). Data are automatically mined from ∼120 sources and presented in an integrated web card for every human gene. We report the application of recent advances in proteomics to enhance gene annotation and classification in GeneCards. First, we constructed the Human Integrated Protein Expression Database (HIPED), a unified database of protein abundance in human tissues, based on the publically available mass spectrometry (MS)-based proteomics sources ProteomicsDB, Multi-Omics Profiling Expression Database, Protein Abundance Across Organisms and The MaxQuant DataBase. The integrated database, residing within GeneCards, compares favourably with its individual sources, covering nearly 90% of human protein-coding genes. For gene annotation and comparisons, we first defined a protein expression vector for each gene, based on normalized abundances in 69 normal human tissues. This vector is portrayed in the GeneCards expression section as a bar graph, allowing visual inspection and comparison. These data are juxtaposed with transcriptome bar graphs. Using the protein expression vectors, we further defined a pairwise metric that helps assess expression-based pairwise proximity. This new metric for finding functional partners complements eight others, including sharing of pathways, gene ontology (GO) terms and domains, implemented in the GeneCards Suite. In parallel, we calculated proteome-based differential expression, highlighting a subset of tissues that overexpress a gene and subserving gene classification. This textual annotation allows users of VarElect, the suite’s next-generation phenotyper, to more effectively discover causative disease variants. Finally, we define the protein–RNA expression ratio and correlation as yet another attribute of every gene in each tissue, adding further annotative information. The results constitute a significant enhancement of several GeneCards sections and help promote and organize the genome-wide structural and functional knowledge of the human proteome.

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

Deficiency of asparagine synthetase causes congenital microcephaly and a progressive form of encephalopathy

We analyzed four families that presented with a similar condition characterized by congenital microcephaly, intellectual disability, progressive cerebral atrophy, and intractable seizures. We show that recessive mutations in the ASNS gene are responsible for this syndrome. Two of the identified missense mutations dramatically reduce ASNS protein abundance, suggesting that the mutations cause loss of function. Hypomorphic Asns mutant mice have structural brain abnormalities, including enlarged ventricles and reduced cortical thickness, and show deficits in learning and memory mimicking aspects of the patient phenotype. ASNS encodes asparagine synthetase, which catalyzes the synthesis of asparagine from glutamine and aspartate. The neurological impairment resulting from ASNS deficiency may be explained by asparagine depletion in the brain or by accumulation of aspartate/glutamate leading to enhanced excitability and neuronal damage. Our study thus indicates that asparagine synthesis is essential for the development and function of the brain but not for that of other organs.

HORDE: comprehensive resource for olfactory receptor genomics

Olfactory receptors (ORs) constitute the largest gene family in the mammalian genome. The existence of these proteins underlies the nature of, and variability in, odorant perception. The Human Olfactory Receptor Data Explorer (HORDE, http://genome.weizmann.ac.il/horde/ ) is a free online resource, which presents a complete compendium of all OR genes and pseudogenes in the genome of human and four other vertebrates. HORDE includes three parts: (1) an automated pipeline, which mines OR gene and pseudogene sequences out of complete genomes, and generates gene symbols based on sequence similarity; (2) a card generator that obtains and displays annotative information on individual ORs retrieved from external databases and relevant studies; and (3) a search engine that allows user retrieval of OR information. For human ORs, HORDE specifically addresses the universe of interindividual variation, as obtained from several sources, including whole genome sequences made possible by next-generation sequencing. This encompasses single nucleotide polymorphisms (SNP) and copy number variation (CNV), including deleterious mutational events. HORDE also hosts a number of tools designed specifically to assist in the study of OR evolution and function. In this chapter, we describe the status of HORDE (build #43). We also discuss plans for future enhancements and a road map for HORDE to become a better community-based bioinformatics tool. We highlight HORDE's role as a major research tool in the study of an expanding cohort of OR repertoires.

Mutation in TECPR2 reveals a role for autophagy in hereditary spastic paraparesis

We studied five individuals from three Jewish Bukharian families affected by an apparently autosomal-recessive form of hereditary spastic paraparesis accompanied by severe intellectual disability, fluctuating central hypoventilation, gastresophageal reflux disease, wake apnea, areflexia, and unique dysmorphic features. Exome sequencing identified one homozygous variant shared among all affected individuals and absent in controls: a 1 bp frameshift TECPR2 deletion leading to a premature stop codon and predicting significant degradation of the protein. TECPR2 has been reported as a positive regulator of autophagy. We thus examined the autophagy-related fate of two key autophagic proteins, SQSTM1 (p62) and MAP1LC3B (LC3), in skin fibroblasts of an affected individual, as compared to a healthy control, and found that both protein levels were decreased and that there was a more pronounced decrease in the lipidated form of LC3 (LC3II). siRNA knockdown of TECPR2 showed similar changes, consistent with aberrant autophagy. Our results are strengthened by the fact that autophagy dysfunction has been implicated in a number of other neurodegenerative diseases. The discovered TECPR2 mutation implicates autophagy, a central intracellular mechanism, in spastic paraparesis.

General olfactory sensitivity database (GOSdb): candidate genes and their genomic variations

Genetic variations in olfactory receptors likely contribute to the diversity of odorant-specific sensitivity phenotypes. Our working hypothesis is that genetic variations in auxiliary olfactory genes, including those mediating transduction and sensory neuronal development, may constitute the genetic basis for general olfactory sensitivity (GOS) and congenital general anosmia (CGA). We thus performed a systematic exploration for auxiliary olfactory genes and their documented variation. This included a literature survey, seeking relevant functional in vitro studies, mouse gene knockouts and human disorders with olfactory phenotypes, as well as data mining in published transcriptome and proteome data for genes expressed in olfactory tissues. In addition, we performed next-generation transcriptome sequencing (RNA-seq) of human olfactory epithelium and mouse olfactory epithelium and bulb, so as to identify sensory-enriched transcripts. Employing a global score system based on attributes of the 11 data sources utilized, we identified a list of 1,680 candidate auxiliary olfactory genes, of which 450 are shortlisted as having higher probability of a functional role. For the top-scoring 136 genes, we identified genomic variants (probably damaging single nucleotide polymorphisms, indels, and copy number deletions) gleaned from public variation repositories. This database of genes and their variants should assist in rationalizing the great interindividual variation in human overall olfactory sensitivity (http://genome.weizmann.ac.il/GOSdb).

Personal receptor repertoires: olfaction as a model

Background

Information on nucleotide diversity along completely sequenced human genomes has increased tremendously over the last few years. This makes it possible to reassess the diversity status of distinct receptor proteins in different human individuals. To this end, we focused on the complete inventory of human olfactory receptor coding regions as a model for personal receptor repertoires.

Results

By performing data-mining from public and private sources we scored genetic variations in 413 intact OR loci, for which one or more individuals had an intact open reading frame. Using 1000 Genomes Project haplotypes, we identified a total of 4069 full-length polypeptide variants encoded by these OR loci, average of ~10 per locus, constituting a lower limit for the effective human OR repertoire. Each individual is found to harbor as many as 600 OR allelic variants, ~50% higher than the locus count. Because OR neuronal expression is allelically excluded, this has direct effect on smell perception diversity of the species. We further identified 244 OR segregating pseudogenes (SPGs), loci showing both intact and pseudogene forms in the population, twenty-six of which are annotatively “resurrected” from a pseudogene status in the reference genome. Using a custom SNP microarray we validated 150 SPGs in a cohort of 468 individuals, with every individual genome averaging 36 disrupted sequence variations, 15 in homozygote form. Finally, we generated a multi-source compendium of 63 OR loci harboring deletion Copy Number Variations (CNVs). Our combined data suggest that 271 of the 413 intact OR loci (66%) are affected by nonfunctional SNPs/indels and/or CNVs.

Conclusions

These results portray a case of unusually high genetic diversity, and suggest that individual humans have a highly personalized inventory of functional olfactory receptors, a conclusion that might apply to other receptor multigene families.

Evolutionary grass roots for odor recognition

Considerable evidence supports the idea that odorant recognition depends on specific sequence variations in olfactory receptor (OR) proteins. Much of this emerges from in vitro screens in heterogenous expression systems. However, the ultimate proof should arise from measurements of odorant thresholds in human individuals harboring different OR genetic variants, a research vein that has so far been only scantly explored. The study of McRae et al., published in this issue of Chemical Senses, shows how the recognition of a grassy odorant depends on specific OR interindividual sequence changes. It provides a clear relevant example for the impact of genetics on olfaction and is an excellent portrayal of the power of human genomics to decipher olfactory perception.

In-silico human genomics with GeneCards

Since 1998, the bioinformatics, systems biology, genomics and medical communities have enjoyed a synergistic relationship with the GeneCards database of human genes (http://www.genecards.org). This human gene compendium was created to help to introduce order into the increasing chaos of information flow. As a consequence of viewing details and deep links related to specific genes, users have often requested enhanced capabilities, such that, over time, GeneCards has blossomed into a suite of tools (including GeneDecks, GeneALaCart, GeneLoc, GeneNote and GeneAnnot) for a variety of analyses of both single human genes and sets thereof. In this paper, we focus on inhouse and external research activities which have been enabled, enhanced, complemented and, in some cases, motivated by GeneCards. In turn, such interactions have often inspired and propelled improvements in GeneCards. We describe here the evolution and architecture of this project, including examples of synergistic applications in diverse areas such as synthetic lethality in cancer, the annotation of genetic variations in disease, omics integration in a systems biology approach to kidney disease, and bioinformatics tools.

Genome sequence of the pattern forming Paenibacillus vortex bacterium reveals potential for thriving in complex environments

Background

The pattern-forming bacterium Paenibacillus vortex is notable for its advanced social behavior, which is reflected in development of colonies with highly intricate architectures. Prior to this study, only two other Paenibacillus species (Paenibacillus sp. JDR-2 and Paenibacillus larvae) have been sequenced. However, no genomic data is available on the Paenibacillus species with pattern-forming and complex social motility. Here we report the de novo genome sequence of this Gram-positive, soil-dwelling, sporulating bacterium.

Results

The complete P. vortex genome was sequenced by a hybrid approach using 454 Life Sciences and Illumina, achieving a total of 289× coverage, with 99.8% sequence identity between the two methods. The sequencing results were validated using a custom designed Agilent microarray expression chip which represented the coding and the non-coding regions. Analysis of the P. vortex genome revealed 6,437 open reading frames (ORFs) and 73 non-coding RNA genes. Comparative genomic analysis with 500 complete bacterial genomes revealed exceptionally high number of two-component system (TCS) genes, transcription factors (TFs), transport and defense related genes. Additionally, we have identified genes involved in the production of antimicrobial compounds and extracellular degrading enzymes.

Conclusions

These findings suggest that P. vortex has advanced faculties to perceive and react to a wide range of signaling molecules and environmental conditions, which could be associated with its ability to reconfigure and replicate complex colony architectures. Additionally, P. vortex is likely to serve as a rich source of genes important for agricultural, medical and industrial applications and it has the potential to advance the study of social microbiology within Gram-positive bacteria.

Systematic inference of copy-number genotypes from personal genome sequencing data reveals extensive olfactory receptor gene content diversity

Copy-number variations (CNVs) are widespread in the human genome, but comprehensive assignments of integer locus copy-numbers (i.e., copy-number genotypes) that, for example, enable discrimination of homozygous from heterozygous CNVs, have remained challenging. Here we present CopySeq, a novel computational approach with an underlying statistical framework that analyzes the depth-of-coverage of high-throughput DNA sequencing reads, and can incorporate paired-end and breakpoint junction analysis based CNV-analysis approaches, to infer locus copy-number genotypes. We benchmarked CopySeq by genotyping 500 chromosome 1 CNV regions in 150 personal genomes sequenced at low-coverage. The assessed copy-number genotypes were highly concordant with our performed qPCR experiments (Pearson correlation coefficient 0.94), and with the published results of two microarray platforms (95–99% concordance). We further demonstrated the utility of CopySeq for analyzing gene regions enriched for segmental duplications by comprehensively inferring copy-number genotypes in the CNV-enriched >800 olfactory receptor (OR) human gene and pseudogene loci. CopySeq revealed that OR loci display an extensive range of locus copy-numbers across individuals, with zero to two copies in some OR loci, and two to nine copies in others. Among genetic variants affecting OR loci we identified deleterious variants including CNVs and SNPs affecting ∼15% and ∼20% of the human OR gene repertoire, respectively, implying that genetic variants with a possible impact on smell perception are widespread. Finally, we found that for several OR loci the reference genome appears to represent a minor-frequency variant, implying a necessary revision of the OR repertoire for future functional studies. CopySeq can ascertain genomic structural variation in specific gene families as well as at a genome-wide scale, where it may enable the quantitative evaluation of CNVs in genome-wide association studies involving high-throughput sequencing.

Human individual genome sequencing has recently become affordable, enabling highly detailed genetic sequence comparisons. While the identification and genotyping of single-nucleotide polymorphisms has already been successfully established for different sequencing platforms, the detection, quantification and genotyping of large-scale copy-number variants (CNVs), i.e., losses or gains of long genomic segments, has remained challenging. We present a computational approach that enables detecting CNVs in sequencing data and accurately identifies the actual copy-number at which DNA segments of interest occur in an individual genome. This approach enabled us to obtain novel insights into the largest human gene family – the olfactory receptors (ORs) – involved in smell perception. While previous studies reported an abundance of CNVs in ORs, our approach enabled us to globally identify absolute differences in OR gene counts that exist between humans. While several OR genes have very high gene counts, other ORs are found only once or are missing entirely in some individuals. The latter have a particularly high probability of influencing individual differences in the perception of smell, a question that future experimental efforts can now address. Furthermore, we observed differences in OR gene counts between populations, pointing at ORs that might contribute to population-specific differences in smell.

Synthetic lethal hubs associated with vincristine resistant neuroblastoma

Chemotherapy of cancer experiences a number of shortcomings including development of drug resistance. This fact also holds true for neuroblastoma utilizing chemotherapeutics as vincristine. We performed a comparative analysis of molecular and cellular mechanisms associated with vincristine resistance utilizing cell line as well as human tissue data. Differential gene expression analysis revealed molecular features, processes and pathways afflicted with drug resistance mechanisms in general, and specifically with vincristine significantly involving actin associated features. However, specific mode of resistance as well as underlying genotype of parental, vincristine sensitive cells apparently exhibited significant heterogeneity. No consensus profile for vincristine resistance could be derived, but resistance-associated changes on the level of individual neuroblastoma cell lines as well as individual patient profiles became clearly evident. Based on these prerequisites we utilized the concept of synthetic lethality aimed at identifying hub proteins which when inhibited promise to induce cell death due to a synthetic lethal interaction with down-regulated, chemoresistance associated features. Our screening procedure identified synthetic lethal hub proteins afflicted with actin associated processes holding synthetic lethal interactions to down-regulated features individually found in all chemoresistant cell lines tested, therefore promising an improved therapeutic window. Verification of such synthetic lethal hub candidates in human neuroblastoma tissue expression profiles indicated the feasibility of this screening approach for addressing vincristine resistance in neuroblastoma.

GeneCards Version 3: the human gene integrator

GeneCards (www.genecards.org) is a comprehensive, authoritative compendium of annotative information about human genes, widely used for nearly 15 years. Its gene-centric content is automatically mined and integrated from over 80 digital sources, resulting in a web-based deep-linked card for each of >73 000 human gene entries, encompassing the following categories: protein coding, pseudogene, RNA gene, genetic locus, cluster and uncategorized. We now introduce GeneCards Version 3, featuring a speedy and sophisticated search engine and a revamped, technologically enabling infrastructure, catering to the expanding needs of biomedical researchers. A key focus is on gene-set analyses, which leverage GeneCards’ unique wealth of combinatorial annotations. These include the GeneALaCart batch query facility, which tabulates user-selected annotations for multiple genes and GeneDecks, which identifies similar genes with shared annotations, and finds set-shared annotations by descriptor enrichment analysis. Such set-centric features address a host of applications, including microarray data analysis, cross-database annotation mapping and gene-disorder associations for drug targeting. We highlight the new Version 3 database architecture, its multi-faceted search engine, and its semi-automated quality assurance system. Data enhancements include an expanded visualization of gene expression patterns in normal and cancer tissues, an integrated alternative splicing pattern display, and augmented multi-source SNPs and pathways sections. GeneCards now provides direct links to gene-related research reagents such as antibodies, recombinant proteins, DNA clones and inhibitory RNAs and features gene-related drugs and compounds lists. We also portray the GeneCards Inferred Functionality Score annotation landscape tool for scoring a gene’s functional information status. Finally, we delineate examples of applications and collaborations that have benefited from the GeneCards suite.

Database URL:www.genecards.org

The genome of a songbird

The zebra finch is an important model organism in several fields with unique relevance to human neuroscience. Like other songbirds, the zebra finch communicates through learned vocalizations, an ability otherwise documented only in humans and a few other animals and lacking in the chicken-the only bird with a sequenced genome until now. Here we present a structural, functional and comparative analysis of the genome sequence of the zebra finch (Taeniopygia guttata), which is a songbird belonging to the large avian order Passeriformes. We find that the overall structures of the genomes are similar in zebra finch and chicken, but they differ in many intrachromosomal rearrangements, lineage-specific gene family expansions, the number of long-terminal-repeat-based retrotransposons, and mechanisms of sex chromosome dosage compensation. We show that song behaviour engages gene regulatory networks in the zebra finch brain, altering the expression of long non-coding RNAs, microRNAs, transcription factors and their targets. We also show evidence for rapid molecular evolution in the songbird lineage of genes that are regulated during song experience. These results indicate an active involvement of the genome in neural processes underlying vocal communication and identify potential genetic substrates for the evolution and regulation of this behaviour.

GeneDecks: paralog hunting and gene-set distillation with GeneCards annotation

Sophisticated genomic navigation strongly benefits from a capacity to establish a similarity metric among genes. GeneDecks is a novel analysis tool that provides such a metric by highlighting shared descriptors between pairs of genes, based on the rich annotation within the GeneCards compendium of human genes. The current implementation addresses information about pathways, protein domains, Gene Ontology (GO) terms, mouse phenotypes, mRNA expression patterns, disorders, drug relationships, and sequence-based paralogy. GeneDecks has two modes: (1) Paralog Hunter, which seeks functional paralogs based on combinatorial similarity of attributes; and (2) Set Distiller, which ranks descriptors by their degree of sharing within a given gene set. GeneDecks enables the elucidation of unsuspected putative functional paralogs, and a refined scrutiny of various gene-sets (e.g., from high-throughput experiments) for discovering relevant biological patterns.

Evidence for an interaction of schizophrenia susceptibility loci on chromosome 6q23.3 and 10q24.33-q26.13 in Arab Israeli families

A genome scan for schizophrenia related loci in Arab Israeli families by Lerer et al. [Lerer et al. (2003); Mol Psychiatry 8:488-498] detected significant evidence for linkage at chromosome 6q23. Subsequent fine mapping [Levi et al. (2005); Eur J Hum Genet 13:763-771], association [Amann-Zalcenstein et al. (2006); Eur J Hum Genet 14:1111-1119] and replication studies [Ingason et al. (2007); Eur J Hum Genet 15:988-991] identified AHI1 as a putative susceptibility gene. The same genome scan revealed suggestive evidence for a schizophrenia susceptibility locus in the 10q23-26 region. Genes at these two loci may act independently in the pathogenesis of the disease in our homogeneous sample of Arab Israeli families or may interact with each other and with other factors in a common biological pathway. The purpose of our current study was to test the hypothesis of genetic interaction between these two loci and to identify the type of interaction between them. The initial stage of our study focused on the 10q23-q26 region which has not been explored further in our sample. The second stage of the study included a test for possible genetic interaction between the 6q23.3 locus and the refined 10q24.33-q26.13 locus. A final candidate region of 19.9 Mb between markers D10S222 (105.3 Mb) and D10S587 (125.2 Mb) was found on chromosome 10 by non-parametric and parametric linkage analyses. These linkage findings are consistent with previous reports in the same chromosomal region. Two-locus multipoint linkage analysis under three complex disease inheritance models (heterogeneity, multiplicative, and additive models) yielded a best maximum LOD score of 7.45 under the multiplicative model suggesting overlapping function of the 6q23.3 and 10q24.33-q26.13 loci.

Mutations and lethality in simulated prebiotic networks

The Graded Autocatalysis Replication Domain (GARD) model describes an origin of life scenario which involves non-covalent compositional assemblies, made of monomeric mutually catalytic molecules. GARD constitutes an alternative to informational biopolymers as a mechanism of primordial inheritance. In the present work, we examined the effect of mutations, one of the most fundamental mechanisms for evolution, in the context of the networks of mutual interaction within GARD prebiotic assemblies. We performed a systematic analysis analogous to single and double gene deletions within GARD. While most deletions have only a small effect on both growth rate and molecular composition of the assemblies, ~10% of the deletions caused lethality, or sometimes showed enhanced fitness. Analysis of 14 different network properties on 2,000 different GARD networks indicated that lethality usually takes place when the deleted node has a high molecular count, or when it is a catalyst for such node. A correlation was also found between lethality and node degree centrality, similar to what is seen in real biological networks. Addressing double knockout mutations, our results demonstrate the occurrence of both synthetic lethality and extragenic suppression within GARD networks, and convey an attempt to correlate synthetic lethality to network node-pair properties. The analyses presented help establish GARD as a workable alternative prebiotic scenario, suggesting that life may have begun with large molecular networks of low fidelity, that later underwent evolutionary compaction and fidelity augmentation.

Common peptides shed light on evolution of Olfactory Receptors

BACKGROUND

Olfactory Receptors (ORs) form the largest multigene family in vertebrates. Their evolution and their expansion in the vertebrate genomes was the subject of many studies. In this paper we apply a motif-based approach to this problem in order to uncover evolutionary characteristics.

RESULTS

We extract deterministic motifs from ORs belonging to ten species using the MEX (Motif Extraction) algorithm, thus defining Common Peptides (CPs) characteristic to ORs. We identify species-specific CPs and show that their relative abundance is high only in fish and frog, suggesting relevance to water-soluble odorants. We estimate the origins of CPs according to the tree of life and track the gains and losses of CPs through evolution. We identify major CP gain in tetrapods and major losses in reptiles. Although the number of human ORs is less than half of the number of ORs in other mammals, the fraction of lost CPs is only 11%.By examining the positions of CPs along the OR sequence, we find two regions that expanded only in tetrapods. Using CPs we are able to establish remote homology relations between ORs and non-OR GPCRs.Selecting CPs according to their evolutionary age, we bicluster ORs and CPs for each species. Clean biclustering emerges when using relatively novel CPs. Evolutionary age is used to track the history of CP acquisition in the collection of mammalian OR families within HORDE (Human Olfactory Receptor Data Explorer).

CONCLUSION

The CP method provides a novel perspective that reveals interesting traits in the evolution of olfactory receptors. It is consistent with previous knowledge, and provides finer details. Using available phylogenetic trees, evolution can be rephrased in terms of CP origins.