Ensembl 2007

The mechanism of enterogenous toxin methylmalonic acid aggravating calcium-phosphorus metabolic disorder in uremic rats by regulating the Wnt/β-catenin pathway

BACKGROUND

Uremia (UR) is caused by increased UR-related toxins in the bloodstream. We explored the mechanism of enterogenous toxin methylmalonic acid (MMA) in calcium-phosphorus metabolic disorder in UR rats via the Wnt/β-catenin pathway.

METHODS

The UR rat model was established by 5/6 nephrectomy. The fecal bacteria of UR rats were transplanted into Sham rats. Sham rats were injected with exogenous MMA or Salinomycin (SAL). Pathological changes in renal/colon tissues were analyzed. MMA concentration, levels of renal function indicators, serum inflammatory factors, Ca2+/P3+, and parathyroid hormone, intestinal flora structure, fecal metabolic profile, intestinal permeability, and glomerular filtration rate (GFR) were assessed. Additionally, rat glomerular podocytes were cultured, with cell viability and apoptosis measured.

RESULTS

Intestinal flora richness and diversity in UR rats were decreased, along with unbalanced flora structure. Among the screened 133 secondary differential metabolites, the MMA concentration rose, showing the most significant difference. UR rat fecal transplantation caused elevated MMA concentration in the serum and renal tissues of Sham rats. The intestinal flora metabolite MMA or exogenous MMA promoted intestinal barrier impairment, increased intestinal permeability, induced glomerular podocyte loss, and reduced GFR, causing calcium-phosphorus metabolic disorder. The intestinal flora metabolite MMA or exogenous MMA induced inflammatory responses and facilitated glomerular podocyte apoptosis by activating the Wnt/β-catenin pathway, which could be counteracted by repressing the Wnt/β-catenin pathway.

CONCLUSIONS

Enterogenous toxin MMA impelled intestinal barrier impairment in UR rats, enhanced intestinal permeability, and activated the Wnt/β-catenin pathway to induce glomerular podocyte loss and reduce GFR, thus aggravating calcium-phosphorus metabolic disorder.

Regulation of eukaryotic transcription initiation in response to cellular stress

Transcription initiation is a vital step in the regulation of eukaryotic gene expression. It can be dysregulated in response to various cellular stressors which is associated with numerous human diseases including cancer. Transcription initiation is facilitated via many gene-specific trans-regulatory elements such as transcription factors, activators, and coactivators through their interactions with transcription pre-initiation complex (PIC). These trans-regulatory elements can uniquely facilitate PIC formation (hence, transcription initiation) in response to cellular nutrient stress. Cellular nutrient stress also regulates the activity of other pathways such as target of rapamycin (TOR) pathway. TOR pathway exhibits distinct regulatory mechanisms of transcriptional activation in response to stress. Like TOR pathway, the cell cycle regulatory pathway is also found to be linked to transcriptional regulation in response to cellular stress. Several transcription factors such as p53, C/EBP Homologous Protein (CHOP), activating transcription factor 6 (ATF6α), E2F, transforming growth factor (TGF)-β, Adenomatous polyposis coli (APC), SMAD, and MYC have been implicated in regulation of transcription of target genes involved in cell cycle progression, apoptosis, and DNA damage repair pathways. Additionally, cellular metabolic and oxidative stressors have been found to regulate the activity of long non-coding RNAs (lncRNA). LncRNA regulates transcription by upregulating or downregulating the transcription regulatory proteins involved in metabolic and cell signaling pathways. Numerous human diseases, triggered by chronic cellular stressors, are associated with abnormal regulation of transcription. Hence, understanding these mechanisms would help unravel the molecular regulatory insights with potential therapeutic interventions. Therefore, here we emphasize the recent advances of regulation of eukaryotic transcription initiation in response to cellular stress.

Highly variable biological effects of statins on cancer, non-cancer, and stem cells in vitro

Statins, the drugs used for the treatment of hypercholesterolemia, have come into the spotlight not only as chemoadjuvants, but also as potential stem cell modulators in the context of regenerative therapy. In our study, we compared the in vitro effects of all clinically used statins on the viability of human pancreatic cancer (MiaPaCa-2) cells, non-cancerous human embryonic kidney (HEK 293) cells and adipose-derived mesenchymal stem cells (ADMSC). Additionally, the effect of statins on viability of MiaPaCa-2 and ADMSC cells spheroids was tested. Furthermore, we performed a microarray analysis on ADMSCs treated with individual statins (12 μM) and compared the importance of the effects of statins on gene expression between stem cells and pancreatic cancer cells. Concentrations of statins that significantly affected cancer cells viability (< 40 μM) did not affect stem cells viability after 24 h. Moreover, statins that didn´t affect viability of cancer cells grown in a monolayer, induce the disintegration of cancer cell spheroids. The effect of statins on gene expression was significantly less pronounced in stem cells compared to pancreatic cancer cells. In conclusion, the low efficacy of statins on non-tumor and stem cells at concentrations sufficient for cancer cells growth inhibition, support their applicability in chemoadjuvant tumor therapy.

Nitric oxide inhibits FTO demethylase activity to regulate N6-methyladenosine mRNA methylation

N6-methyladenosine (m6A) is the most abundant internal modification on eukaryotic mRNAs. Demethylation of m6A on mRNA is catalyzed by the enzyme fat mass and obesity-associated protein (FTO), a member of the nonheme Fe(II) and 2-oxoglutarate (2-OG)-dependent family of dioxygenases. FTO activity and m6A-mRNA are dysregulated in multiple diseases including cancers, yet endogenous signaling molecules that modulate FTO activity have not been identified. Here we show that nitric oxide (NO) is a potent inhibitor of FTO demethylase activity by directly binding to the catalytic iron center, which causes global m6A hypermethylation of mRNA in cells and results in gene-specific enrichment of m6A on mRNA of NO-regulated transcripts. Both cell culture and tumor xenograft models demonstrated that endogenous NO synthesis can regulate m6A-mRNA levels and transcriptional changes of m6A-associated genes. These results build a direct link between NO and m6A-mRNA regulation and reveal a novel signaling mechanism of NO as an endogenous regulator of the epitranscriptome.

A review from biological mapping to computation-based subcellular localization

Subcellular localization is crucial to the study of virus and diseases. Specifically, research on protein subcellular localization can help identify clues between virus and host cells that can aid in the design of targeted drugs. Research on RNA subcellular localization is significant for human diseases (such as Alzheimer's disease, colon cancer, etc.). To date, only reviews addressing subcellular localization of proteins have been published, which are outdated for reference, and reviews of RNA subcellular localization are not comprehensive. Therefore, we collated (the most up-to-date) literature on protein and RNA subcellular localization to help researchers understand changes in the field of protein and RNA subcellular localization. Extensive and complete methods for constructing subcellular localization models have also been summarized, which can help readers understand the changes in application of biotechnology and computer science in subcellular localization research and explore how to use biological data to construct improved subcellular localization models. This paper is the first review to cover both protein subcellular localization and RNA subcellular localization. We urge researchers from biology and computational biology to jointly pay attention to transformation patterns, interrelationships, differences, and causality of protein subcellular localization and RNA subcellular localization.

Pharmacogenetics of Praziquantel Metabolism: Evaluating the Cytochrome P450 Genes of Zimbabwean Patients During a Schistosomiasis Treatment

Schistosomiasis is a parasitic disease infecting over 236 million people annually, with the majority affected residing on the African continent. Control of this disease is reliant on the drug praziquantel (PZQ), with treatment success dependent on an individual reaching PZQ concentrations lethal to schistosomes. Despite the complete reliance on PZQ to treat schistosomiasis in Africa, the characterization of the pharmacogenetics associated with PZQ metabolism in African populations has been sparse. We aimed to characterize genetic variation in the drug-metabolising cytochrome P450 enzymes (CYPs) and determine the association between each variant and the efficacy of PZQ treatment in Zimbabwean patients exposed to Schistosoma haematobium infection. Genomic DNA from blood samples of 114 case-control Zimbabweans infected with schistosomes were sequenced using the CYP1A2, CYP2C9, CYP2C19, CYP2D6, CYP3A4, and CYP3A5 genes as targets. Bioinformatic tools were used to identify and predict functional effects of detected single nucleotide polymorphisms (SNPs). A random forest (RF) model was then used to assess SNPs most predictive of PZQ efficacy, with a misclassification rate of 29%. SNPs were detected across all six genes, with 70 SNPs identified and multiple functional changes to the CYP enzymes predicted. Only four SNPs were significantly associated with PZQ efficacy using χ² tests, with rs951840747 (OR: 3.61, p = 0.01) in the CYP1A2 gene having the highest odds of an individual possessing this SNP clearing infection, and rs6976017 (OR: 2.19, p = 0.045) of CYP3A5 determined to be the most predictive of PZQ efficacy via the RF. Only the rs28371702 (CC) genotype (OR: 2.36, p = 0.024) of CYP2D6 was significantly associated with an unsuccessful PZQ treatment. This study adds to the genomic characterization of the diverse populations in Africa and identifies variants relevant to other pharmacogenetic studies crucial for the development and usage of drugs in these populations.

PANTHER: Making genome-scale phylogenetics accessible to all

Phylogenetics is a powerful tool for analyzing protein sequences, by inferring their evolutionary relationships to other proteins. However, phylogenetics analyses can be challenging: they are computationally expensive and must be performed carefully in order to avoid systematic errors and artifacts. Protein Analysis THrough Evolutionary Relationships (PANTHER; http://pantherdb.org) is a publicly available, user-focused knowledgebase that stores the results of an extensive phylogenetic reconstruction pipeline that includes computational and manual processes and quality control steps. First, fully reconciled phylogenetic trees (including ancestral protein sequences) are reconstructed for a set of "reference" protein sequences obtained from fully sequenced genomes of organisms across the tree of life. Second, the resulting phylogenetic trees are manually reviewed and annotated with function evolution events: inferred gains and losses of protein function along branches of the phylogenetic tree. Here, we describe in detail the current contents of PANTHER, how those contents are generated, and how they can be used in a variety of applications. The PANTHER knowledgebase can be downloaded or accessed via an extensive API. In addition, PANTHER provides software tools to facilitate the application of the knowledgebase to common protein sequence analysis tasks: exploring an annotated genome by gene function; performing "enrichment analysis" of lists of genes; annotating a single sequence or large batch of sequences by homology; and assessing the likelihood that a genetic variant at a particular site in a protein will have deleterious effects.

Co-transcriptional RNA cleavage by Drosha homolog Pac1 triggers transcription termination in fission yeast

Transcription termination of protein-coding genes in eukaryotic cells usually relies on a tight coordination between the cleavage and polyadenylation of the pre-mRNA, and 5′-3′ degradation of the downstream nascent transcript. Here we investigated the contribution of the essential fission yeast endonuclease Pac1, a homolog of human Drosha that cleaves hairpin RNA structures, in triggering polyadenylation-independent transcription termination. Using ChIP-sequencing in Pac1-deficient cells, we found that Pac1 triggers transcription termination at snRNA and snoRNA genes as well as at specific protein-coding genes. Notably, we found that Pac1-dependent premature termination occurred at two genes encoding conserved transmembrane transporters whose expression were strongly repressed by Pac1. Analysis by genome editing indicated that a stem-loop structure in the nascent transcript directs Pac1-mediated cleavage and that the regions upstream and downstream of the Pac1 cleavage site in the targeted mRNAs were stabilized by mutation of nuclear 3′-5′ and 5′-3′ exonucleases, respectively. Our findings unveil a premature transcription termination pathway that uncouples co-transcriptional RNA cleavage from polyadenylation, triggering rapid nuclear RNA degradation.

Comparison of Transcriptomic Profiles of MiaPaCa-2 Pancreatic Cancer Cells Treated with Different Statins

Statins have been widely used for the treatment of hypercholesterolemia due to their ability to inhibit HMG-CoA reductase, the rate-limiting enzyme of de novo cholesterol synthesis, via the so-called mevalonate pathway. However, their inhibitory action also causes depletion of downstream intermediates of the pathway, resulting in the pleiotropic effects of statins, including the beneficial impact in the treatment of cancer. In our study, we compared the effect of all eight existing statins on the expression of genes, the products of which are implicated in cancer inhibition and suggested the molecular mechanisms of their action in epigenetic and posttranslational regulation, and in cell-cycle arrest, death, migration, or invasion of the cancer cells.

ANO1 Expression Orchestrates p27Kip1/MCL1-Mediated Signaling in Head and Neck Squamous Cell Carcinoma

Simple Summary

Our aim was to elucidate the molecular mechanisms of how ANO1 contributes to oncogenic processes in squamous cell carcinoma of the head and neck (HNSCC). We explored transcriptional programs influenced by ANO1 knockdown in patient-derived UT-SCC cell lines with 11q13 amplification and ANO1 overexpression. ANO1 depletion led to downregulation of broad pro-survival BCL2 family protein members, including MCL1, and simultaneously induced upregulation of the cell cycle inhibitor p27^Kip1 and its redistribution from the cytoplasm into the nucleus in the studied HNSCC cells. Gene set enrichment analysis highlighted pathways associated with perturbed cell cycle and apoptosis in the ANO1-depleted samples. Silencing of ANO1 and application of an ANO1-targeting small-molecule inhibitor led to ANO1 degradation and reduction of cell viability. These findings suggest that ANO1 has drug target potential that deserves further evaluation in preclinical in vivo models.

Abstract

Head and neck squamous cell carcinoma (HNSCC) is a heterogeneous group of tumors that derive from the mucosal epithelium of the upper aerodigestive tract and present high mortality rate. Lack of efficient targeted-therapies and biomarkers towards patients’ stratification are caveats in the disease treatment. Anoctamin 1 (ANO1) gene is amplified in 30% of HNSCC cases. Evidence suggests involvement of ANO1 in proliferation, migration, and evasion of apoptosis; however, the exact mechanisms remain elusive. Aim of this study was to unravel the ANO1-dependent transcriptional programs and expand the existing knowledge of ANO1 contribution to oncogenesis and drug response in HNSCC. We cultured two HNSCC cell lines established from primary tumors harboring amplification and high expression of ANO1 in three-dimensional collagen. Differential expression analysis of ANO1-depleted HNSCC cells demonstrated downregulation of MCL1 and simultaneous upregulation of p27^Kip1 expression. Suppressing ANO1 expression led to redistribution of p27^Kip1 from the cytoplasm to the nucleus and associated with a cell cycle arrested phenotype. ANO1 silencing or pharmacological inhibition resulted in reduction of cell viability and ANO1 protein levels, as well as suppression of pro-survival BCL2 family proteins. Collectively, these data provide insights of ANO1 involvement in HNSCC carcinogenesis and support the rationale that ANO1 is an actionable drug target.

SARS-CoV-2 Proteins Exploit Host's Genetic and Epigenetic Mediators for the Annexation of Key Host Signaling Pathways

The constant rise of the death toll and cases of COVID-19 has made this pandemic a serious threat to human civilization. Understanding of host-SARS-CoV-2 interaction in viral pathogenesis is still in its infancy. In this study, we utilized a blend of computational and knowledgebase approaches to model the putative virus-host interplay in host signaling pathways by integrating the experimentally validated host interactome proteins and differentially expressed host genes in SARS-CoV-2 infection. While searching for the pathways in which viral proteins interact with host proteins, we discovered various antiviral immune response pathways such as hypoxia-inducible factor 1 (HIF-1) signaling, autophagy, retinoic acid-inducible gene I (RIG-I) signaling, Toll-like receptor signaling, fatty acid oxidation/degradation, and IL-17 signaling. All these pathways can be either hijacked or suppressed by the viral proteins, leading to improved viral survival and life cycle. Aberration in pathways such as HIF-1 signaling and relaxin signaling in the lungs suggests the pathogenic lung pathophysiology in COVID-19. From enrichment analysis, it was evident that the deregulated genes in SARS-CoV-2 infection might also be involved in heart development, kidney development, and AGE-RAGE signaling pathway in diabetic complications. Anomalies in these pathways might suggest the increased vulnerability of COVID-19 patients with comorbidities. Moreover, we noticed several presumed infection-induced differentially expressed transcription factors and epigenetic factors, such as miRNAs and several histone modifiers, which can modulate different immune signaling pathways, helping both host and virus. Our modeling suggests that SARS-CoV-2 integrates its proteins in different immune signaling pathways and other cellular signaling pathways for developing efficient immune evasion mechanisms while leading the host to a more complicated disease condition. Our findings would help in designing more targeted therapeutic interventions against SARS-CoV-2.

Transcriptome of nasopharyngeal samples from COVID-19 patients and a comparative analysis with other SARS-CoV-2 infection models reveal disparate host responses against SARS-CoV-2

BACKGROUND

Although it is becoming evident that individual's immune system has a decisive influence on SARS-CoV-2 disease progression, pathogenesis is largely unknown. In this study, we aimed to profile the host transcriptome of COVID-19 patients from nasopharyngeal samples along with virus genomic features isolated from respective host, and a comparative analyses of differential host responses in various SARS-CoV-2 infection systems.

RESULTS

Unique and rare missense mutations in 3C-like protease observed in all of our reported isolates. Functional enrichment analyses exhibited that the host induced responses are mediated by innate immunity, interferon, and cytokine stimulation. Surprisingly, induction of apoptosis, phagosome, antigen presentation, hypoxia response was lacking within these patients. Upregulation of immune and cytokine signaling genes such as CCL4, TNFA, IL6, IL1A, CCL2, CXCL2, IFN, and CCR1 were observed in lungs. Lungs lacked the overexpression of ACE2 as suspected, however, high ACE2 but low DPP4 expression was observed in nasopharyngeal cells. Interestingly, directly or indirectly, viral proteins specially non-structural protein mediated overexpression of integrins such as ITGAV, ITGA6, ITGB7, ITGB3, ITGA2B, ITGA5, ITGA6, ITGA9, ITGA4, ITGAE, and ITGA8 in lungs compared to nasopharyngeal samples suggesting the possible way of enhanced invasion. Furthermore, we found comparatively highly expressed transcription factors such as CBP, CEBP, NFAT, ATF3, GATA6, HDAC2, TCF12 which have pivotal roles in lung injury.

CONCLUSIONS

Even though this study incorporates a limited number of cases, our data will provide valuable insights in developing potential studies to elucidate the differential host responses on the viral pathogenesis in COVID-19, and incorporation of further data will enrich the search of an effective therapeutics.

Vertebrate Chromosome Evolution

The study of chromosome evolution is undergoing a resurgence of interest owing to advances in DNA sequencing technology that facilitate the production of chromosome-scale whole-genome assemblies de novo. This review focuses on the history, methods, discoveries, and current challenges facing the field, with an emphasis on vertebrate genomes. A detailed examination of the literature on the biology of chromosome rearrangements is presented, specifically the relationship between chromosome rearrangements and phenotypic evolution, adaptation, and speciation. A critical review of the methods for identifying, characterizing, and visualizing chromosome rearrangements and computationally reconstructing ancestral karyotypes is presented. We conclude by looking to the future, identifying the enormous technical and scientific challenges presented by the accumulation of hundreds and eventually thousands of chromosome-scale assemblies.

Lung transcriptome of a COVID-19 patient and systems biology predictions suggest impaired surfactant production which may be druggable by surfactant therapy

An incomplete understanding of the molecular mechanisms behind impairment of lung pathobiology by COVID-19 complicates its clinical management. In this study, we analyzed the gene expression pattern of cells obtained from biopsies of COVID-19-affected patient and compared to the effects observed in typical SARS-CoV-2 and SARS-CoV-infected cell-lines. We then compared gene expression patterns of COVID-19-affected lung tissues and SARS-CoV-2-infected cell-lines and mapped those to known lung-related molecular networks, including hypoxia induced responses, lung development, respiratory processes, cholesterol biosynthesis and surfactant metabolism; all of which are suspected to be downregulated following SARS-CoV-2 infection based on the observed symptomatic impairments. Network analyses suggest that SARS-CoV-2 infection might lead to acute lung injury in COVID-19 by affecting surfactant proteins and their regulators SPD, SPC, and TTF1 through NSP5 and NSP12; thrombosis regulators PLAT, and EGR1 by ORF8 and NSP12; and mitochondrial NDUFA10, NDUFAF5, and SAMM50 through NSP12. Furthermore, hypoxia response through HIF-1 signaling might also be targeted by SARS-CoV-2 proteins. Drug enrichment analysis of dysregulated genes has allowed us to propose novel therapies, including lung surfactants, respiratory stimulants, sargramostim, and oseltamivir. Our study presents a distinct mechanism of probable virus induced lung damage apart from cytokine storm.

Khademul Islam AB. Perversely expressed long noncoding RNAs can alter host response and viral proliferation in SARS-CoV-2 infection

BACKGROUND

Regulatory roles of long noncoding RNAs (lncRNAs) during viral infection has become more evident in last decade, but are yet to be explored for SARS-CoV-2.

MATERIALS & METHODS

We analyzed RNA-seq dataset of SARS-CoV-2 infected lung epithelial cells to identify differentially expressed genes.

RESULTS

Our analyses uncover 21 differentially expressed lncRNAs broadly involved in cell survival and regulation of gene expression. These lncRNAs can directly interact with six differentially expressed protein-coding genes, and ten host genes that interact with SARS-CoV-2 proteins. Also, they can block the suppressive effect of nine microRNAs induced in viral infections.

CONCLUSION

Our investigation determines that deregulated lncRNAs in SARS-CoV-2 infection are involved in viral proliferation, cellular survival, and immune response, ultimately determining disease outcome.

Epigenetic Regulator miRNA Pattern Differences Among SARS-CoV, SARS-CoV-2, and SARS-CoV-2 World-Wide Isolates Delineated the Mystery Behind the Epic Pathogenicity and Distinct Clinical Characteristics of Pandemic COVID-19

A detailed understanding of the molecular mechanism of SARS-CoV-2 pathogenesis is still elusive, and there is a need to address its deadly nature and to design effective therapeutics. Here, we present a study that elucidates the interplay between the SARS-CoV and SARS-CoV-2 viruses' and host's miRNAs, an epigenetic regulator, as a mode of pathogenesis; and we explored how the SARS-CoV and SARS-CoV-2 infections differ in terms of their miRNA-mediated interactions with the host and the implications this has in terms of disease complexity. We have utilized computational approaches to predict potential host and viral miRNAs and their possible roles in different important functional pathways. We have identified several putative host antiviral miRNAs that can target the SARS viruses and also predicted SARS viruses-encoded miRNAs targeting host genes. In silico predicted targets were also integrated with SARS-infected human cell microarray and RNA-seq gene expression data. A comparison between the host miRNA binding profiles on 67 different SARS-CoV-2 genomes from 24 different countries with respective country's normalized death count surprisingly uncovered some miRNA clusters, which are associated with increased death rates. We have found that induced cellular miRNAs can be both a boon and a bane to the host immunity, as they have possible roles in neutralizing the viral threat; conversely, they can also function as proviral factors. On the other hand, from over representation analysis, our study revealed that although both SARS-CoV and SARS-CoV-2 viral miRNAs could target broad immune-signaling pathways; only some of the SARS-CoV-2 miRNAs are found to uniquely target some immune-signaling pathways, such as autophagy, IFN-I signaling, etc., which might suggest their immune-escape mechanisms for prolonged latency inside some hosts without any symptoms of COVID-19. Furthermore, SARS-CoV-2 can modulate several important cellular pathways that might lead to the increased anomalies in patients with comorbidities like cardiovascular diseases, diabetes, breathing complications, etc. This might suggest that miRNAs can be a key epigenetic modulator behind the overcomplications amongst the COVID-19 patients. Our results support that miRNAs of host and SARS-CoV-2 can indeed play a role in the pathogenesis which can be further concluded with more experiments. These results will also be useful in designing RNA therapeutics to alleviate the complications from COVID-19.

Rbf Activates the Myogenic Transcriptional Program to Promote Skeletal Muscle Differentiation

The importance of the retinoblastoma tumor suppressor protein pRB in cell cycle control is well established. However, less is known about its role in differentiation during animal development. Here, we investigated the role of Rbf, the Drosophila pRB homolog, in adult skeletal muscles. We found that the depletion of Rbf severely reduced muscle growth and altered myofibrillogenesis but only minimally affected myoblast proliferation. We identified an Rbf-dependent transcriptional program in late muscle development that is distinct from the canonical role of Rbf in cell cycle control. Unexpectedly, Rbf acts as a transcriptional activator of the myogenic and metabolic genes in the growing muscles. The genomic regions bound by Rbf contained the binding sites of several factors that genetically interacted with Rbf by modulating Rbf-dependent phenotype. Thus, our results reveal a distinctive role for Rbf as a direct activator of the myogenic transcriptional program that drives late muscle differentiation.

Inactivation of the tumor suppressor RB, an obligatory step in most cancers, results in unrestrained cell cycle progression. Zappia et al. show that Rbf, the RB Drosophila ortholog, directly activates the metabolic program that accompanies muscle development. This work expands the understanding of the plethora of Rbf functions.

Exploratory analysis of genetic variants influencing molecular traits in cerebral cortex of suicide completers

Genetic factors have been implicated in suicidal behavior. It has been suggested that one of the roles of genetic factors in suicide could be represented by the effect of genetic variants on gene expression regulation. Alteration in the expression of genes participating in multiple biological systems in the suicidal brain has been demonstrated, so it is imperative to identify genetic variants that could influence gene expression or its regulatory mechanisms. In this study, we integrated DNA methylation, gene expression, and genotype data from the prefrontal cortex of suicides to identify genetic variants that could be factors in the regulation of gene expression, generally called quantitative trait locus (xQTLs). We identify 6,224 methylation quantitative trait loci and 2,239 expression quantitative trait loci (eQTLs) in the prefrontal cortex of suicide completers. The xQTLs identified influence the expression of genes involved in neurodevelopment and cell organization. Two of the eQTLs identified (rs8065311 and rs1019238) were previously associated with cannabis dependence, highlighting a candidate genetic variant for the increased suicide risk in subjects with substance use disorders. Our findings suggest that genetic variants may regulate gene expression in the prefrontal cortex of suicides through the modulation of promoter and enhancer activity, and to a lesser extent, binding transcription factors.

Deep Intact Proteoform Characterization in Human Cell Lysate Using High-pH and Low-pH Reversed-Phase Liquid Chromatography

Post-translational modifications (PTMs) play critical roles in biological processes and have significant effects on the structures and dynamics of proteins. Top-down proteomics methods were developed for and applied to the study of intact proteins and their PTMs in human samples. However, the large dynamic range and complexity of human samples makes the study of human proteins challenging. To address these challenges, we developed a 2D pH RP/RPLC-MS/MS technique that fuses high-resolution separation and intact protein characterization to study the human proteins in HeLa cell lysate. Our results provide a deep coverage of soluble proteins in human cancer cells. Compared to 225 proteoforms from 124 proteins identified when 1D separation was used, 2778 proteoforms from 628 proteins were detected and characterized using our 2D separation method. Many proteoforms with critically functional PTMs including phosphorylation were characterized. Additionally, we present the first detection of intact human GcvH proteoforms with rare modifications such as octanoylation and lipoylation. Overall, the increase in the number of proteoforms identified using 2DLC separation is largely due to the reduction in sample complexity through improved separation resolution, which enables the detection of low-abundance PTM-modified proteoforms. We demonstrate here that 2D pH RP/RPLC is an effective technique to analyze complex protein samples using top-down proteomics.

Immunogenetics of marsupial B-cells

Marsupials and eutherians are mammals that differ in their physiological traits, predominately their reproductive and developmental strategies; eutherians give birth to well-developed young, while marsupials are born highly altricial after a much shorter gestation. These developmental traits also result in differences in the development of the immune system of eutherian and marsupial species. In eutherians, B-cells are the key to humoral immunity as they are found in multiple lymphoid organs and have the unique ability to mediate the production of antigen-specific antibodies in the presence of extracellular pathogens. The development of B-cells in marsupials has been reported and hypothesised to be similar to that of eutherians, except that haematopoiesis occurs in the liver, postpartum, until the bone marrow fully matures. In eutherians, specific genes are linked to specific stages in B-cell development, maturation, and differentiation processes, and have been identified including immunoglobulins (heavy and light chains), cluster of differentiation markers (CD10, 19, 34 and CD79α/β), signal transduction molecules (BTK, Lyn and Syk) and transcriptional regulators (EBF1, E2A, and Pax5). This review aims to discuss the known similarities and differences between marsupial and eutherian B-cells, in regards to their genetic presence, homology, and developmental stages, as well as to highlight the areas requiring further investigation. By enhancing our understanding of the genes that are involved with B-cells in the marsupial lineage, it will, in turn, aid our understanding of the marsupial immune system and support the development of specific immunological reagents for research and wildlife conservation purposes.

SwiftOrtho: A fast, memory-efficient, multiple genome orthology classifier

BACKGROUND

Gene homology type classification is required for many types of genome analyses, including comparative genomics, phylogenetics, and protein function annotation. Consequently, a large variety of tools have been developed to perform homology classification across genomes of different species. However, when applied to large genomic data sets, these tools require high memory and CPU usage, typically available only in computational clusters.

FINDINGS

Here we present a new graph-based orthology analysis tool, SwiftOrtho, which is optimized for speed and memory usage when applied to large-scale data. SwiftOrtho uses long k-mers to speed up homology search, while using a reduced amino acid alphabet and spaced seeds to compensate for the loss of sensitivity due to long k-mers. In addition, it uses an affinity propagation algorithm to reduce the memory usage when clustering large-scale orthology relationships into orthologous groups. In our tests, SwiftOrtho was the only tool that completed orthology analysis of proteins from 1,760 bacterial genomes on a computer with only 4 GB RAM. Using various standard orthology data sets, we also show that SwiftOrtho has a high accuracy.

CONCLUSIONS

SwiftOrtho enables the accurate comparative genomic analyses of thousands of genomes using low-memory computers. SwiftOrtho is available at https://github.com/Rinoahu/SwiftOrtho.

UniProt genomic mapping for deciphering functional effects of missense variants

Understanding the association of genetic variation with its functional consequences in proteins is essential for the interpretation of genomic data and identifying causal variants in diseases. Integration of protein function knowledge with genome annotation can assist in rapidly comprehending genetic variation within complex biological processes. Here, we describe mapping UniProtKB human sequences and positional annotations, such as active sites, binding sites, and variants to the human genome (GRCh38) and the release of a public genome track hub for genome browsers. To demonstrate the power of combining protein annotations with genome annotations for functional interpretation of variants, we present specific biological examples in disease-related genes and proteins. Computational comparisons of UniProtKB annotations and protein variants with ClinVar clinically annotated single nucleotide polymorphism (SNP) data show that 32% of UniProtKB variants colocate with 8% of ClinVar SNPs. The majority of colocated UniProtKB disease-associated variants (86%) map to 'pathogenic' ClinVar SNPs. UniProt and ClinVar are collaborating to provide a unified clinical variant annotation for genomic, protein, and clinical researchers. The genome track hubs, and related UniProtKB files, are downloadable from the UniProt FTP site and discoverable as public track hubs at the UCSC and Ensembl genome browsers.

GenTree, an integrated resource for analyzing the evolution and function of primate-specific coding genes

The origination of new genes contributes to phenotypic evolution in humans. Two major challenges in the study of new genes are the inference of gene ages and annotation of their protein-coding potential. To tackle these challenges, we created GenTree, an integrated online database that compiles age inferences from three major methods together with functional genomic data for new genes. Genome-wide comparison of the age inference methods revealed that the synteny-based pipeline (SBP) is most suited for recently duplicated genes, whereas the protein-family–based methods are useful for ancient genes. For SBP-dated primate-specific protein-coding genes (PSGs), we performed manual evaluation based on published PSG lists and showed that SBP generated a conservative data set of PSGs by masking less reliable syntenic regions. After assessing the coding potential based on evolutionary constraint and peptide evidence from proteomic data, we curated a list of 254 PSGs with different levels of protein evidence. This list also includes 41 candidate misannotated pseudogenes that encode primate-specific short proteins. Coexpression analysis showed that PSGs are preferentially recruited into organs with rapidly evolving pathways such as spermatogenesis, immune response, mother–fetus interaction, and brain development. For brain development, primate-specific KRAB zinc-finger proteins (KZNFs) are specifically up-regulated in the mid-fetal stage, which may have contributed to the evolution of this critical stage. Altogether, hundreds of PSGs are either recruited to processes under strong selection pressure or to processes supporting an evolving novel organ.

Relevance of Titin Missense and Non-Frameshifting Insertions/Deletions Variants in Dilated Cardiomyopathy

Recent advancements in next generation sequencing (NGS) technology have led to the identification of the giant sarcomere gene, titin (TTN), as a major human disease gene. Truncating variants of TTN (TTNtv) especially in the A-band region account for 20% of dilated cardiomyopathy (DCM) cases. Much attention has been focused on assessment and interpretation of TTNtv in human disease; however, missense and non-frameshifting insertions/deletions (NFS-INDELs) are difficult to assess and interpret in clinical diagnostic workflow. Targeted sequencing covering all exons of TTN was performed on a cohort of 530 primary DCM patients from three cardiogenetic centres across Europe. Using stringent bioinformatic filtering, twenty-nine and two rare TTN missense and NFS-INDELs variants predicted deleterious were identified in 6.98% and 0.38% of DCM patients, respectively. However, when compared with those identified in the largest available reference population database, no significant enrichment of such variants was identified in DCM patients. Moreover, DCM patients and reference individuals had comparable frequencies of splice-region missense variants with predicted splicing alteration. DCM patients and reference populations had comparable frequencies of rare predicted deleterious TTN missense variants including splice-region missense variants suggesting that these variants are not independently causative for DCM. Hence, these variants should be classified as likely benign in the clinical diagnostic workflow, although a modifier effect cannot be excluded at this stage.

Predictions of novel Schistosoma mansoni - human protein interactions consistent with experimental data

Infection by the human blood fluke, Schistosoma mansoni involves a variety of cross-species protein- protein interactions. The pathogen expresses a diverse arsenal of proteins that facilitate the breach of physical and biochemical barriers present in skin evasion of the immune system, and digestion of human plasma proteins including albumin and hemoglobin, allowing schistosomes to reside in the host for years. However, only a small number of specific interactions between S. mansoni and human proteins have been identified. We present and apply a protocol that generates testable predictions of S. mansoni-human protein interactions. In this study, we have preliminary predictions of novel interactions between schistosome and human proteins relevant to infection and the ability of the parasite to evade the immune system. We applied a computational whole-genome comparative approach to predict potential S. mansoni-human protein interactions based on similarity to known protein complexes. We first predict S. mansoni -human protein interactions based on similarity to known protein complexes. Putative interactions were then scored and assessed using several contextual filters, including the use of annotation automatically derived from literature using a simple natural language processing methodology. Next, in vitro experiments were carried out between schistosome and host proteins to validate several prospective predictions. Our method predicted 7 out of the 10 previously known cross-species interactions involved in pathogenesis between S. mansoni and its human host. Interestingly, two novel putative interactions involving Schistosoma proteins, the cercarial elastase SmCE, and the adult tegument surface protein Sm29, were also predicted and experimentally characterized. Preliminary data suggest that elafin, a host endogenous serine protease inhibitor, may be a novel substrate for SmCE. Additionally, CD59, an inhibitor of the membrane attack complex, could interact with Sm29. Furthermore, the application framework provides an integrated methodology for investigation of host-pathogen interactions and an extensive source of orthogonal data for experimental analysis. We have made the predictions available for community perusal.

Simultaneous sequencing of coding and noncoding RNA reveals a human transcriptome dominated by a small number of highly expressed noncoding genes

Comparing the abundance of one RNA molecule to another is crucial for understanding cellular functions but most sequencing techniques can target only specific subsets of RNA. In this study, we used a new fragmented ribodepleted TGIRT sequencing method that uses a thermostable group II intron reverse transcriptase (TGIRT) to generate a portrait of the human transcriptome depicting the quantitative relationship of all classes of nonribosomal RNA longer than 60 nt. Comparison between different sequencing methods indicated that FRT is more accurate in ranking both mRNA and noncoding RNA than viral reverse transcriptase-based sequencing methods, even those that specifically target these species. Measurements of RNA abundance in different cell lines using this method correlate with biochemical estimates, confirming tRNA as the most abundant nonribosomal RNA biotype. However, the single most abundant transcript is 7SL RNA, a component of the signal recognition particle. Structured noncoding RNAs (sncRNAs) associated with the same biological process are expressed at similar levels, with the exception of RNAs with multiple functions like U1 snRNA. In general, sncRNAs forming RNPs are hundreds to thousands of times more abundant than their mRNA counterparts. Surprisingly, only 50 sncRNA genes produce half of the non-rRNA transcripts detected in two different cell lines. Together the results indicate that the human transcriptome is dominated by a small number of highly expressed sncRNAs specializing in functions related to translation and splicing.

Proteomic Analysis Shows Constitutive Secretion of MIF and p53-associated Activity of COX-2-/- Lung Fibroblasts

The differential expression of two closelyassociated cyclooxygenase isozymes, COX-1 and COX-2, exhibited functions beyond eicosanoid metabolism. We hypothesized that COX-1 or COX-2 knockout lung fibroblasts may display altered protein profiles which may allow us to further differentiate the functional roles of these isozymes at the molecular level. Proteomic analysis shows constitutive production of macrophage migration inhibitory factor (MIF) in lung fibroblasts derived from COX-2^−/− but not wild-type (WT) or COX-1^−/− mice. MIF was spontaneously released in high levels into the extracellular milieu of COX2^−/− fibroblasts seemingly from the preformed intracellular stores, with no change in the basal gene expression of MIF. The secretion and regulation of MIF in COX-2^−/− was “prostaglandin-independent.” GO analysis showed that concurrent with upregulation of MIF, there is a significant surge in expression of genes related to fibroblast growth, FK506 binding proteins, and isomerase activity in COX-2^−/− cells. Furthermore, COX-2^−/− fibroblasts also exhibit a significant increase in transcriptional activity of various regulators, antagonists, and co-modulators of p53, as well as in the expression of oncogenes and related transcripts. Integrative Oncogenomics Cancer Browser (IntroGen) analysis shows downregulation of COX-2 and amplification of MIF and/or p53 activity during development of glioblastomas, ependymoma, and colon adenomas. These data indicate the functional role of the MIF-COX-p53 axis in inflammation and cancer at the genomic and proteomic levels in COX-2-ablated cells. This systematic analysis not only shows the proinflammatory state but also unveils a molecular signature of a pro-oncogenic state of COX-1 in COX-2 ablated cells.

ChIP-Seq analysis identifies p27(Kip1)-target genes involved in cell adhesion and cell signalling in mouse embryonic fibroblasts

The protein p27^Kip1 (p27), a member of the Cip-Kip family of cyclin-dependent kinase inhibitors, is involved in tumorigenesis and a correlation between reduced levels of this protein in human tumours and a worse prognosis has been established. Recent reports revealed that p27 also behaves as a transcriptional regulator. Thus, it has been postulated that the development of tumours with low amounts of p27 could be propitiated by deregulation of transcriptional programs under the control of p27. However, these programs still remain mostly unknown. The aim of this study has been to define the transcriptional programs regulated by p27 by first identifying the p27-binding sites (p27-BSs) on the whole chromatin of quiescent mouse embryonic fibroblasts. The chromatin regions associated to p27 have been annotated to the most proximal genes and it has been considered that the expression of these genes could by regulated by p27. The identification of the chromatin p27-BSs has been performed by Chromatin Immunoprecipitation Sequencing (ChIP-seq). Results revealed that p27 associated with 1839 sites that were annotated to 1417 different genes being 852 of them protein coding genes. Interestingly, most of the p27-BSs were in distal intergenic regions and introns whereas, in contrast, its association with promoter regions was very low. Gene ontology analysis of the protein coding genes revealed a number of relevant transcriptional programs regulated by p27 as cell adhesion, intracellular signalling and neuron differentiation among others. We validated the interaction of p27 with different chromatin regions by ChIP followed by qPCR and demonstrated that the expressions of several genes belonging to these programs are actually regulated by p27. Finally, cell adhesion assays revealed that the adhesion of p27^-/- cells to the plates was much higher that controls, revealing a role of p27 in the regulation of a transcriptional program involved in cell adhesion.

Two Distinct Regulatory Mechanisms of Transcriptional Initiation in Response to Nutrient Signaling

SAGA (Spt-Ada-Gcn5-Acetyltransferase) and TFIID (transcription factor IID) have been previously shown to facilitate the formation of the PIC (pre-initiation complex) at the promoters of two distinct sets of genes. Here, we demonstrate that TFIID and SAGA differentially participate in the stimulation of PIC formation (and hence transcriptional initiation) at the promoter of PHO84, a gene for the high-affinity inorganic phosphate (P_i) transporter for crucial cellular functions, in response to nutrient signaling. We show that transcriptional initiation of PHO84 occurs predominantly in a TFIID-dependent manner in the absence of P_i in the growth medium. Such TFIID dependency is mediated via the NuA4 (nucleosome acetyltransferase of H4) histone acetyltransferase (HAT). Intriguingly, transcriptional initiation of PHO84 also occurs in the presence of P_i in the growth medium, predominantly via the SAGA complex, but independently of NuA4 HAT. Thus, P_i in the growth medium switches transcriptional initiation of PHO84 from NuA4-TFIID to SAGA dependency. Further, we find that both NuA4-TFIID- and SAGA-dependent transcriptional initiations of PHO84 are facilitated by the 19S proteasome subcomplex or regulatory particle (RP) via enhanced recruitment of the coactivators SAGA and NuA4 HAT, which promote TFIID-independent and -dependent PIC formation for transcriptional initiation, respectively. NuA4 HAT does not regulate activator binding to PHO84, but rather facilitates PIC formation for transcriptional initiation in the absence of Pi in the growth medium. On the other hand, SAGA promotes activator recruitment to PHO84 for transcriptional initiation in the growth medium containing Pi. Collectively, our results demonstrate two distinct stimulatory pathways for PIC formation (and hence transcriptional initiation) at PHO84 by TFIID, SAGA, NuA4, and 19S RP in the presence and absence of an essential nutrient, P_i, in the growth media, thus providing new regulatory mechanisms of transcriptional initiation in response to nutrient signaling.

Identification of Novel Retroid Agents in Danio rerio, Oryzias latipes, Gasterosteus aculeatus and Tetraodon nigroviridis

Retroid agents are genomes that encode a reverse transcriptase (RT) and replicate or transpose by way of an RNA intermediate. The Genome Parsing Suite (GPS) is software created to identify and characterize Retroid agents in any genome database (McClure et al. 2005). The detailed analysis of all Retroid agents found by the GPS in Danio rerio (zebrafish), Oryzias latipes (medaka), Gasterosteus aculeatus (stickleback) and Tetraodon nigroviridis (spotted green pufferfish) reveals extensive Retroid agent diversity in the compact genomes of all four fish. Novel Retroid agents were identified by the GPS software: the telomerase reverse transcriptase (TERT) in O. latipes, G aculeatus and T. nigroviridis and a potential TERT in D. rerio, a retrotransposon in D. rerio, and multiple lineages of endogenous retroviruses (ERVs) in D. rerio, O. latipes and G aculeatus.

Public Databases and Software for the Pathway Analysis of Cancer Genomes

The study of pathway disruption is key to understanding cancer biology. Advances in high throughput technologies have led to the rapid accumulation of genomic data. The explosion in available data has generated opportunities for investigation of concerted changes that disrupt biological functions, this in turns created a need for computational tools for pathway analysis. In this review, we discuss approaches to the analysis of genomic data and describe the publicly available resources for studying biological pathways.

p27Kip1, PCAF and PAX5 cooperate in the transcriptional regulation of specific target genes

The cyclin-dependent kinase inhibitor p27^Kip1 (p27) also behaves as a transcriptional repressor. Data showing that the p300/CBP-associated factor (PCAF) acetylates p27 inducing its degradation suggested that PCAF and p27 could collaborate in the regulation of transcription. However, this possibility remained to be explored. We analyzed here the transcriptional programs regulated by PCAF and p27 in the colon cancer cell line HCT116 by chromatin immunoprecipitation sequencing (ChIP-seq). We identified 269 protein-encoding genes that contain both p27 and PCAF binding sites being the majority of these sites different for PCAF and p27. PCAF or p27 knock down revealed that both regulate the expression of these genes, PCAF as an activator and p27 as a repressor. The double knock down of PCAF and p27 strongly reduced their expression indicating that the activating role of PCAF overrides the repressive effect of p27. We also observed that the transcription factor Pax5 interacts with both p27 and PCAF and that the knock down of Pax5 induces the expression of p27/PCAF target genes indicating that it also participates in the transcriptional regulation mediated by p27/PCAF. In summary, we report here a previously unknown mechanism of transcriptional regulation mediated by p27, Pax5 and PCAF.

Hook, Line and Infection: A Guide to Culturing Parasites, Establishing Infections and Assessing Immune Responses in the Three-Spined Stickleback

The three-spined stickleback (Gasterosteus aculeatus) is a model organism with an extremely well-characterized ecology, evolutionary history, behavioural repertoire and parasitology that is coupled with published genomic data. These small temperate zone fish therefore provide an ideal experimental system to study common diseases of coldwater fish, including those of aquacultural importance. However, detailed information on the culture of stickleback parasites, the establishment and maintenance of infections and the quantification of host responses is scattered between primary and grey literature resources, some of which is not readily accessible. Our aim is to lay out a framework of techniques based on our experience to inform new and established laboratories about culture techniques and recent advances in the field. Here, essential knowledge on the biology, capture and laboratory maintenance of sticklebacks, and their commonly studied parasites is drawn together, highlighting recent advances in our understanding of the associated immune responses. In compiling this guide on the maintenance of sticklebacks and a range of common, taxonomically diverse parasites in the laboratory, we aim to engage a broader interdisciplinary community to consider this highly tractable model when addressing pressing questions in evolution, infection and aquaculture.

RNA splicing in human disease and in the clinic

Defects at the level of the pre-mRNA splicing process represent a major cause of human disease. Approximately 15-50% of all human disease mutations have been shown to alter functioning of basic and auxiliary splicing elements. These elements are required to ensure proper processing of pre-mRNA splicing molecules, with their disruption leading to misprocessing of the pre-mRNA molecule and disease. The splicing process is a complex process, with much still to be uncovered before we are able to accurately predict whether a reported genomic sequence variant (GV) represents a splicing-associated disease mutation or a harmless polymorphism. Furthermore, even when a mutation is correctly identified as affecting the splicing process, there still remains the difficulty of providing an exact evaluation of the potential impact on disease onset, severity and duration. In this review, we provide a brief overview of splicing diagnostic methodologies, from in silico bioinformatics approaches to wet lab in vitro and in vivo systems to evaluate splicing efficiencies. In particular, we provide an overview of how the latest developments in high-throughput sequencing can be applied to the clinic, and are already changing clinical approaches.

Mutation in cpsf6/CFIm68 (Cleavage and Polyadenylation Specificity Factor Subunit 6) causes short 3'UTRs and disturbs gene expression in developing embryos, as revealed by an analysis of primordial germ cell migration using the medaka mutant naruto

Our previous studies analyzing medaka mutants defective in primordial germ cell (PGC) migration identified cxcr4b and cxcr7, which are both receptors of the chemokine sdf1/cxcl12, as key regulators of PGC migration. Among PGC migration mutants, naruto (nar) is unique in that the mutant phenotype includes gross morphological abnormalities of embryos, suggesting that the mutation affects a broader range of processes. A fine genetic linkage mapping and genome sequencing showed the nar gene encodes Cleavage and Polyadenylation Specificity Factor subunit 6 (CPSF6/CFIm68). CPSF6 is a component of the Cleavage Factor Im complex (CFIm) which plays a key role in pre-mRNA 3'-cleavage and polyadenylation. 3'RACE of sdf1a/b and cxcr7 transcripts in the mutant embryos indicated shorter 3'UTRs with poly A additions occurring at more upstream positions than wild-type embryos, suggesting CPSF6 functions to prevent premature 3'UTR cleavage. In addition, expression of the coding region sequences of sdf1a/b in nar mutants was more anteriorly extended in somites than wild-type embryos, accounting for the abnormally extended distribution of PGCs in nar mutants. An expected consequence of shortening 3'UTR is the escape from the degradation mechanism mediated by microRNAs interacting with distal 3'UTR sequence. The abnormal expression pattern of sdf1a coding sequence may be at least partially accounted for by this mechanism. Given the pleiotropic effects of nar mutation, further analysis using the nar mutant will reveal processes in which CPSF6 plays essential regulatory roles in poly A site selection and involvement of 3'UTRs in posttranscriptional gene regulation in various genes in vivo.

Prediction and Validation of Transcription Factors Modulating the Expression of Sestrin3 Gene Using an Integrated Computational and Experimental Approach

SESN3 has been implicated in multiple biological processes including protection against oxidative stress, regulation of glucose and lipid metabolism. However, little is known about the factors and mechanisms controlling its gene expression at the transcriptional level. We performed in silico phylogenetic footprinting analysis of 5 kb upstream regions of a diverse set of human SESN3 orthologs for the identification of high confidence conserved binding motifs (BMo). We further analyzed the predicted BMo by a motif comparison tool to identify the TFs likely to bind these discovered motifs. Predicted TFs were then integrated with experimentally known protein-protein interactions and experimentally validated to delineate the important transcriptional regulators of SESN3. Our study revealed high confidence set of BMos (integrated with DNase I hypersensitivity sites) in the upstream regulatory regions of SESN3 that could be bound by transcription factors from multiple families including FOXOs, SMADs, SOXs, TCFs and HNF4A. TF-TF network analysis established hubs of interaction that include SMAD3, TCF3, SMAD2, HDAC2, SOX2, TAL1 and TCF12 as well as the likely protein complexes formed between them. We show using ChIP-PCR as well as over-expression and knock out studies that FOXO3 and SOX2 transcriptionally regulate the expression of SESN3 gene. Our findings provide an important roadmap to further our understanding on the regulation of SESN3.

Investigating polymorphisms by bioinformatics is a potential cost-effective method to screen for germline mutations in Chinese familial adenomatous polyposis patients

The aim of this study was to investigate germline mutations of the APC, MUTYH and AXIN2 genes in Chinese patients with familial adenomatous polyposis (FAP), and further assess the value of bioinformatics in screening the pathogenic changes predisposing to FAP. APC genes from 11 unrelated FAP patients in Yunnan province in China were firstly examined by exon-specific DNA sequencing. For samples without already known pathogenic changes predisposing to FAP in the APC gene, whole-gene sequencing of MUTYH and AXIN2 was performed. Mutational analysis of each gene was performed by bioinformatics. Eleven different types of APC polymorphisms were observed in the cohort of families analyzed. Of these polymorphisms, four were missense substitutions (V1822D, V1173G, P1760H and K2057), one was a nonsense substitution (S1196X), and six were silent substitutions (Y486Y, T449T, T1493T, G1678G, S1756S and P1960P). One missense mutation (Q335H) and two intronic substitutions (c.264+11G>A and c.420+35A>G) were detected in the MUTYH gene, and four synonymous mutations (I144I, P455P, P462P and L688L) and three intonic mutations (c.1060–77G>T, c.1060–287A>G and c.1060–282 A>G) of the AXIN2 gene were observed. In addition to the already reported pathogenic mutations, by using function assessment tools and databases, the synonymous substitutions observed in the APC gene of our samples were predicted to affect splicing regulation in the translation of mRNA, while the missense mutations observed in the APC gene and MUTYH gene were predicted to be disease-related polymorphisms; however, no functional effect of the mutations was observed in the AXIN2 gene. Comprehensive screening for germline mutations in APC, MUTYH and AXIN2 genes followed by prediction of pathogenicity using bioinformatic tools contributes to a cost-effective way of screening germline mutations in Chinese familial adenomatous polyposis patients.

Inflammatory stimuli acutely modulate peripheral taste function

Inflammation-mediated changes in taste perception can affect health outcomes in patients, but little is known about the underlying mechanisms. In the present work, we hypothesized that proinflammatory cytokines directly modulate Na(+) transport in taste buds. To test this, we measured acute changes in Na(+) flux in polarized fungiform taste buds loaded with a Na(+) indicator dye. IL-1β elicited an amiloride-sensitive increase in Na(+) transport in taste buds. In contrast, TNF-α dramatically and reversibly decreased Na(+) flux in polarized taste buds via amiloride-sensitive and amiloride-insensitive Na(+) transport systems. The speed and partial amiloride sensitivity of these changes in Na(+) flux indicate that IL-1β and TNF-α modulate epithelial Na(+) channel (ENaC) function. A portion of the TNF-mediated decrease in Na(+) flux is also blocked by the TRPV1 antagonist capsazepine, although TNF-α further reduced Na(+) transport independently of both amiloride and capsazepine. We also assessed taste function in vivo in a model of infection and inflammation that elevates these and additional cytokines. In rats administered systemic lipopolysaccharide (LPS), CT responses to Na(+) were significantly elevated between 1 and 2 h after LPS treatment. Low, normally preferred concentrations of NaCl and sodium acetate elicited high response magnitudes. Consistent with this outcome, codelivery of IL-1β and TNF-α enhanced Na(+) flux in polarized taste buds. These results demonstrate that inflammation elicits swift changes in Na(+) taste function, which may limit salt consumption during illness.

Bulk Segregant Analysis Reveals the Genetic Basis of a Natural Trait Variation in Fission Yeast

Although the fission yeast Schizosaccharomyces pombe is a well-established model organism, studies of natural trait variations in this species remain limited. To assess the feasibility of segregant-pool-based mapping of phenotype-causing genes in natural strains of fission yeast, we investigated the cause of a maltose utilization defect (Mal(-)) of the S. pombe strain CBS5557 (originally known as Schizosaccharomyces malidevorans). Analyzing the genome sequence of CBS5557 revealed 955 nonconservative missense substitutions, and 61 potential loss-of-function variants including 47 frameshift indels, 13 early stop codons, and 1 splice site mutation. As a side benefit, our analysis confirmed 146 sequence errors in the reference genome and improved annotations of 27 genes. We applied bulk segregant analysis to map the causal locus of the Mal(-) phenotype. Through sequencing the segregant pools derived from a cross between CBS5557 and the laboratory strain, we located the locus to within a 2.23-Mb chromosome I inversion found in most S. pombe isolates including CBS5557. To map genes within the inversion region that occupies 18% of the genome, we created a laboratory strain containing the same inversion. Analyzing segregants from a cross between CBS5557 and the inversion-containing laboratory strain narrowed down the locus to a 200-kb interval and led us to identify agl1, which suffers a 5-bp deletion in CBS5557, as the causal gene. Interestingly, loss of agl1 through a 34-kb deletion underlies the Mal(-) phenotype of another S. pombe strain CGMCC2.1628. This work adapts and validates the bulk segregant analysis method for uncovering trait-gene relationship in natural fission yeast strains.

Investigation of glycan evolution based on a comprehensive analysis of glycosyltransferases using phylogenetic profiling

Glycans play important roles in such cell-cell interactions as signaling and adhesion, including processes involved in pathogenic infections, cancers, and neurological diseases. Glycans are biosynthesized by multiple glycosyltransferases (GTs), which function sequentially. Excluding mucin-type O-glycosylation, the non-reducing terminus of glycans is biosynthesized in the Golgi apparatus after the reducing terminus is biosynthesized in the ER. In the present study, we performed genome-wide analyses of human GTs by investigating the degree of conservation of homologues in other organisms, as well as by elucidating the phylogenetic relationship between cephalochordates and urochordates, which has long been controversial in deuterostome phylogeny. We analyzed 173 human GTs and functionally linked glycan synthesis enzymes by phylogenetic profiling and clustering, compiled orthologous genes from the genomes of other organisms, and converted them into a binary sequence based on the presence (1) or absence (0) of orthologous genes in the genomes. Our results suggest that the non-reducing terminus of glycans is biosynthesized by newly evolved GTs. According to our analysis, the phylogenetic profiles of GTs resemble the phylogenetic tree of life, where deuterostomes, metazoans, and eukaryotes are resolved into separate branches. Lineage-specific GTs appear to play essential roles in the divergence of these particular lineages. We suggest that urochordates lose several genes that are conserved among metazoans, such as those expressing sialyltransferases, and that the Golgi apparatus acquires the ability to synthesize glycans after the ER acquires this function.

Whole exome sequencing in extended families with autism spectrum disorder implicates four candidate genes

Autism spectrum disorders (ASDs) are a group of neurodevelopmental disorders, characterized by impairment in communication and social interactions, and by repetitive behaviors. ASDs are highly heritable, and estimates of the number of risk loci range from hundreds to >1000. We considered 7 extended families (size 12-47 individuals), each with ≥3 individuals affected by ASD. All individuals were genotyped with dense SNP panels. A small subset of each family was typed with whole exome sequence (WES). We used a 3-step approach for variant identification. First, we used family-specific parametric linkage analysis of the SNP data to identify regions of interest. Second, we filtered variants in these regions based on frequency and function, obtaining exactly 200 candidates. Third, we compared two approaches to narrowing this list further. We used information from the SNP data to impute exome variant dosages into those without WES. We regressed affected status on variant allele dosage, using pedigree-based kinship matrices to account for relationships. The p value for the test of the null hypothesis that variant allele dosage is unrelated to phenotype was used to indicate strength of evidence supporting the variant. A cutoff of p = 0.05 gave 28 variants. As an alternative third filter, we required Mendelian inheritance in those with WES, resulting in 70 variants. The imputation- and association-based approach was effective. We identified four strong candidate genes for ASD (SEZ6L, HISPPD1, FEZF1, SAMD11), all of which have been previously implicated in other studies, or have a strong biological argument for their relevance.

Genome-wide ultraconserved elements exhibit higher phylogenetic informativeness than traditional gene markers in percomorph fishes

Ultraconserved elements (UCEs) have become popular markers in phylogenomic studies because of their cost effectiveness and their potential to resolve problematic phylogenetic relationships. Although UCE datasets typically contain a much larger number of loci and sites than more traditional datasets of PCR-amplified, single-copy, protein coding genes, a fraction of UCE sites are expected to be part of a nearly invariant core, and the relative performance of UCE datasets versus protein coding gene datasets is poorly understood. Here we use phylogenetic informativeness (PI) to compare the resolving power of multi-locus and UCE datasets in a sample of percomorph fishes with sequenced genomes (genome-enabled). We compare three data sets: UCE core regions, flanking sequence adjacent to the UCE core and a set of ten protein coding genes commonly used in fish systematics. We found the net informativeness of UCE core and flank regions to be roughly ten-fold and 100-fold more informative than that of the protein coding genes. On a per locus basis UCEs and protein coding genes exhibited similar levels of phylogenetic informativeness. Our results suggest that UCEs offer enormous potential for resolving relationships across the percomorph tree of life.

VEGF mRNA and protein concentrations in the developing human eye

BACKGROUND

Vascular endothelial growth factor (VEGF), a well-characterized regulator of angiogenesis, has been mechanistically implicated in retinal neovascularization and in the pathogenesis of retinopathy of prematurity. However, the ontogeny of VEGF expression in the human fetal retina is not well known. Because retinal vasculature grows with gestational maturation, we hypothesized that VEGF expression also increases in the midgestation human fetal eye as a function of gestational age.

METHODS

To identify changes in VEGF gene expression during normal human development, we measured VEGF mRNA by quantitative PCR and measured VEGF protein by enzyme-linked immunosorbent assay and western blots in 10-24 wk gestation fetal vitreous, retina, and serum.

RESULTS

VEGF mRNA expression in the retina increased with gestational age. VEGF isoform A, particularly its VEGF121 splice variant, contributed to this positive correlation. Consistent with these findings, we detected increasing VEGF121 protein concentrations in vitreous humor from fetuses of 10-24 wk gestation, while VEGF concentrations decreased in fetal serum.

CONCLUSION

VEGF121 mRNA and protein concentrations increase with increasing gestational age in the developing human retina. We speculate that VEGF plays an important role in normal retinal vascular development, and that preterm delivery affects production of this vascular growth factor.

Unraveling the association between mRNA expressions and mutant phenotypes in a genome-wide assessment of mice

High-throughput gene expression profiling has revealed substantial leaky and extraneous transcription of eukaryotic genes, challenging the perceptions that transcription is strictly regulated and that changes in transcription have phenotypic consequences. To assess the functional implications of mRNA transcription directly, we analyzed mRNA expression data derived from microarrays, RNA-sequencing, and in situ hybridization, together with phenotype data of mouse mutants as a proxy of gene function at the tissue level. The results indicated that despite the presence of widespread ectopic transcription, mRNA expression and mutant phenotypes of mammalian genes or tissues remain associated. The expression-phenotype association at the gene level was particularly strong for tissue-specific genes, and the association could be underestimated due to data insufficiency and incomprehensive phenotyping of mouse mutants; the strength of expression-phenotype association at the tissue level depended on tissue functions. Mutations on genes expressed at higher levels or expressed at earlier embryonic stages more often result in abnormal phenotypes in the tissues where they are expressed. The mRNA expression profiles that have stronger associations with their phenotype profiles tend to be more evolutionarily conserved, indicating that the evolution of transcriptome and the evolution of phenome are coupled. Therefore, mutations resulting in phenotypic aberrations in expressed tissues are more likely to occur in highly transcribed genes, tissue-specific genes, genes expressed during early embryonic stages, or genes with evolutionarily conserved mRNA expression profiles.

Identification of Cdca7 as a novel Notch transcriptional target involved in hematopoietic stem cell emergence

Guiu et al. use ChIP-on-chip analysis for the Notch partner RBPj, using embryonic tissue from the aorta-gonad-mesonephros region to identify potential novel Notch target genes involved in HSC emergence. They show that c-MYC–responsive gene Cdca7 is expressed in different HSC and progenitor subpopulations and that CDCA7 is important for maintaining the undifferentiated phenotype. Cdca7 acts downstream of Notch in HSCs in zebrafish, mouse, and human, indicating a highly conserved Notch/RBPj/Cdca7 axis in hematopoietic development.

Hematopoietic stem cell (HSC) specification occurs in the embryonic aorta and requires Notch activation; however, most of the Notch-regulated elements controlling de novo HSC generation are still unknown. Here, we identify putative direct Notch targets in the aorta-gonad-mesonephros (AGM) embryonic tissue by chromatin precipitation using antibodies against the Notch partner RBPj. By ChIP-on-chip analysis of the precipitated DNA, we identified 701 promoter regions that were candidates to be regulated by Notch in the AGM. One of the most enriched regions corresponded to the Cdca7 gene, which was subsequently confirmed to recruit the RBPj factor but also Notch1 in AGM cells. We found that during embryonic hematopoietic development, expression of Cdca7 is restricted to the hematopoietic clusters of the aorta, and it is strongly up-regulated in the hemogenic population during human embryonic stem cell hematopoietic differentiation in a Notch-dependent manner. Down-regulation of Cdca7 mRNA in cultured AGM cells significantly induces hematopoietic differentiation and loss of the progenitor population. Finally, using loss-of-function experiments in zebrafish, we demonstrate that CDCA7 contributes to HSC emergence in vivo during embryonic development. Thus, our study identifies Cdca7 as an evolutionary conserved Notch target involved in HSC emergence.