SET-PP2A complex as a new therapeutic target in KMT2A (MLL) rearranged AML

KMT2A-rearranged (KMT2A-R) is an aggressive and chemo-refractory acute leukemia which mostly affects children. Transcriptomics-based characterization and chemical interrogation identified kinases as key drivers of survival and drug resistance in KMT2A-R leukemia. In contrast, the contribution and regulation of phosphatases is unknown. In this study we uncover the essential role and underlying mechanisms of SET, the endogenous inhibitor of Ser/Thr phosphatase PP2A, in KMT2A-R-leukemia. Investigation of SET expression in acute myeloid leukemia (AML) samples demonstrated that SET is overexpressed, and elevated expression of SET is correlated with poor prognosis and with the expression of MEIS and HOXA genes in AML patients. Silencing SET specifically abolished the clonogenic ability of KMT2A-R leukemic cells and the transcription of KMT2A targets genes HOXA9 and HOXA10. Subsequent mechanistic investigations showed that SET interacts with both KMT2A wild type and fusion proteins, and it is recruited to the HOXA10 promoter. Pharmacological inhibition of SET by FTY720 disrupted SET-PP2A interaction leading to cell cycle arrest and increased sensitivity to chemotherapy in KMT2A-R-leukemic models. Phospho-proteomic analyses revealed that FTY720 reduced the activity of kinases regulated by PP2A, including ERK1, GSK3β, AURB and PLK1 and led to suppression of MYC, supporting the hypothesis of a feedback loop among PP2A, AURB, PLK1, MYC, and SET. Our findings illustrate that SET is a novel player in KMT2A-R leukemia and they provide evidence that SET antagonism could serve as a novel strategy to treat this aggressive leukemia.

microRNA sequencing for biomarker detection in the diagnosis, classification and prognosis of Diffuse Large B Cell Lymphoma

Despite being considered a single disease, Diffuse Large B Cell Lymphoma (DLBCL) presents with variable backgrounds, which results in heterogeneous outcomes among patients, with 40% of them still having primary refractory disease or relapse. Thus, novel biomarkers are needed. In addition, multiple factors regarding its pathogenesis remain unclear. In this context, recent investigations point to the relevance of microRNAs (miRNAs) in cancer. However, regarding DLBCL, there is inconsistency in the data reported. Therefore, in this work, the main goals were to determine a miRNA set with utility as biomarkers for DLBCL diagnosis, classification, prognosis and treatment response, as well as to decipher the mechanism of action of deregulated miRNAs in the origin of the disease. We analyzed miRNA expression in a cohort of 78 DLBCL patients and 17 controls using small RNA sequencing and performed a miRNA-mRNA interaction network analysis. This way, we were able to define new miRNA expression signatures for diagnosis, classification, treatment response and prognosis, and we identified plausible mechanisms of action by which deregulated miRNAs could be involved in DLBCL pathogenesis. In summary, our study remarks that miRNAs could play an important role in DLBCL.

Nanopore Direct RNA Sequencing Data Processing and Analysis Using MasterOfPores

This chapter describes MasterOfPores v.2 (MoP2), an open-source suite of pipelines for processing and analyzing direct RNA Oxford Nanopore sequencing data. The MoP2 relies on the Nextflow DSL2 framework and Linux containers, thus enabling reproducible data analysis in transcriptomic and epitranscriptomic studies. We introduce the key concepts of MoP2 and provide a step-by-step fully reproducible and complete example of how to use the workflow for the analysis of S. cerevisiae total RNA samples sequenced using MinION flowcells. The workflow starts with the pre-processing of raw FAST5 files, which includes basecalling, read quality control, demultiplexing, filtering, mapping, estimation of per-gene/transcript abundances, and transcriptome assembly, with support of the GPU computing for the basecalling and read demultiplexing steps. The secondary analyses of the workflow focus on the estimation of RNA poly(A) tail lengths and the identification of RNA modifications. The MoP2 code is available at https://github.com/biocorecrg/MOP2 and is distributed under the MIT license.

Genome-wide transcriptomic changes reveal the genetic pathways involved in insect migration

Insects are capable of extraordinary feats of long‐distance movement that have profound impacts on the function of terrestrial ecosystems. The ability to undertake these movements arose multiple times through the evolution of a suite of traits that make up the migratory syndrome, however the underlying genetic pathways involved remain poorly understood. Migratory hoverflies (Diptera: Syrphidae) are an emerging model group for studies of migration. They undertake seasonal movements in huge numbers across large parts of the globe and are important pollinators, biological control agents and decomposers. Here, we assembled a high‐quality draft genome of the marmalade hoverfly (Episyrphus balteatus). We leveraged this genomic resource to undertake a genome‐wide transcriptomic comparison of actively migrating Episyrphus, captured from a high mountain pass as they flew south to overwinter, with the transcriptomes of summer forms which were non‐migratory. We identified 1543 genes with very strong evidence for differential expression. Interrogation of this gene set reveals a remarkable range of roles in metabolism, muscle structure and function, hormonal regulation, immunity, stress resistance, flight and feeding behaviour, longevity, reproductive diapause and sensory perception. These features of the migrant phenotype have arisen by the integration and modification of pathways such as insulin signalling for diapause and longevity, JAK/SAT for immunity, and those leading to octopamine production and fuelling to boost flight capabilities. Our results provide a powerful genomic resource for future research, and paint a comprehensive picture of global expression changes in an actively migrating insect, identifying key genomic components involved in this important life‐history strategy.

Controlled X-chromosome dynamics defines meiotic potential of female mouse in vitro germ cells

The mammalian germline is characterized by extensive epigenetic reprogramming during its development into functional eggs and sperm. Specifically, the epigenome requires resetting before parental marks can be established and transmitted to the next generation. In the female germline, X‐chromosome inactivation and reactivation are among the most prominent epigenetic reprogramming events, yet very little is known about their kinetics and biological function. Here, we investigate X‐inactivation and reactivation dynamics using a tailor‐made in vitro system of primordial germ cell‐like cell (PGCLC) differentiation from mouse embryonic stem cells. We find that X‐inactivation in PGCLCs in vitro and in germ cell‐competent epiblast cells in vivo is moderate compared to somatic cells, and frequently characterized by escaping genes. X‐inactivation is followed by step‐wise X‐reactivation, which is mostly completed during meiotic prophase I. Furthermore, we find that PGCLCs which fail to undergo X‐inactivation or reactivate too rapidly display impaired meiotic potential. Thus, our data reveal fine‐tuned X‐chromosome remodelling as a critical feature of female germ cell development towards meiosis and oogenesis.

Selective cleavage of ncRNA and antiviral activity by RNase2/EDN in THP1-induced macrophages

RNase2 is the member of the RNaseA family most abundant in macrophages. Here, we knocked out RNase2 in THP-1 cells and analysed the response to Respiratory Syncytial Virus (RSV). RSV induced RNase2 expression, which significantly enhanced cell survival. Next, by cP-RNAseq sequencing, which amplifies the cyclic-phosphate endonuclease products, we analysed the ncRNA population. Among the ncRNAs accumulated in WT vs KO cells, we found mostly tRNA-derived fragments (tRFs) and second miRNAs. Differential sequence coverage identified tRFs from only few parental tRNAs, revealing a predominant cleavage at anticodon and d-loops at U/C (B1) and A (B2) sites. Selective tRNA cleavage was confirmed in vitro using the recombinant protein. Likewise, only few miRNAs were significantly more abundant in WT vs RNase2-KO cells. Complementarily, by screening of a tRF & tiRNA array, we identified an enriched population associated to RNase2 expression and RSV exposure. The results confirm the protein antiviral action and provide the first evidence of its cleavage selectivity on ncRNAs.

Graphical abstract

Cross talk between the upstream exon-intron junction and Prp2 facilitates splicing of non-consensus introns

Splicing of mRNA precursors is essential in the regulation of gene expression. U2AF65 recognizes the poly-pyrimidine tract and helps in the recognition of the branch point. Inactivation of fission yeast U2AF65 (Prp2) blocks splicing of most, but not all, pre-mRNAs, for reasons that are not understood. Here, we have determined genome-wide the splicing efficiency of fission yeast cells as they progress into synchronous meiosis in the presence or absence of functional Prp2. Our data indicate that in addition to the splicing elements at the 3' end of any intron, the nucleotides immediately upstream the intron will determine whether Prp2 is required or dispensable for splicing. By changing those nucleotides in any given intron, we regulate its Prp2 dependency. Our results suggest a model in which Prp2 is required for the coordinated recognition of both intronic ends, placing Prp2 as a key regulatory element in the determination of the exon-intron boundaries.

ExOrthist: a tool to infer exon orthologies at any evolutionary distance

Several bioinformatic tools have been developed for genome-wide identification of orthologous and paralogous genes. However, no corresponding tool allows the detection of exon homology relationships. Here, we present ExOrthist, a fully reproducible Nextflow-based software enabling inference of exon homologs and orthogroups, visualization of evolution of exon-intron structures, and assessment of conservation of alternative splicing patterns. ExOrthist evaluates exon sequence conservation and considers the surrounding exon-intron context to derive genome-wide multi-species exon homologies at any evolutionary distance. We demonstrate its use in different evolutionary scenarios: whole genome duplication in frogs and convergence of Nova-regulated splicing networks (https://github.com/biocorecrg/ExOrthist).

NFAT5 Amplifies Antipathogen Responses by Enhancing Chromatin Accessibility, H3K27 Demethylation, and Transcription Factor Recruitment

The ability of innate immune cells to respond to pathogen-associated molecular patterns across a wide range of intensities is fundamental to limit the spreading of infections. Studies on transcription responses to pathogen-activated TLRs have often used relatively high TLR ligand concentrations, and less is known about their regulation under mild stimulatory conditions. We had shown that the transcription factor NFAT5 facilitates expression of antipathogen genes under TLR stimulation conditions corresponding to low pathogen loads. In this study, we analyze how NFAT5 optimizes TLR-activated responses in mouse macrophages. We show that NFAT5 was required for effective recruitment of central effectors p65/NF-κB and c-Fos to specific proinflammatory target genes, such as Nos2, Il6, and Tnf in primary macrophages responding to low doses of the TLR4 ligand LPS. By contrast, NFAT5 was not required for p65/NF-κB recruitment in response to high LPS doses. Using the transposase-accessible chromatin with high-throughput sequencing assay, we show that NFAT5 facilitated chromatin accessibility mainly at promoter regions of multiple TLR4-responsive genes. Analysis of various histone marks that regulate gene expression in response to pathogens identified H3K27me3 demethylation as an early NFAT5-dependent mechanism that facilitates p65 recruitment to promoters of various TLR4-induced genes. Altogether, these results advance our understanding about specific mechanisms that optimize antipathogen responses to limit infections.

Citizen-science based study of the oral microbiome in Cystic fibrosis and matched controls reveals major differences in diversity and abundance of bacterial and fungal species

Introduction: Cystic fibrosis (CF) is an autosomal genetic disease, associated with the production of excessively thick mucosa and with life-threatening chronic lung infections. The microbiota of the oral cavity can act as a reservoir or as a barrier for infectious microorganisms that can colonize the lungs. However, the specific composition of the oral microbiome in CF is poorly understood.Methods: In collaboration with CF associations in Spain, we collected oral rinse samples from 31 CF persons (age range 7-47) and matched controls, and then performed 16S rRNA metabarcoding and high-throughput sequencing, combined with culture and proteomics-based identification of fungi to survey the bacterial and fungal oral microbiome.Results: We found that CF is associated with less diverse oral microbiomes, which were characterized by higher prevalence of Candida albicans and differential abundances of a number of bacterial taxa that have implications in both the connection to lung infections in CF, as well as potential oral health concerns, particularly periodontitis and dental caries.Conclusion: Overall, our study provides a first global snapshot of the oral microbiome in CF. Future studies are required to establish the relationships between the composition of the oral and lung microbiomes in CF.

QCloud2: An Improved Cloud-based Quality-Control System for Mass-Spectrometry-based Proteomics Laboratories

QCloud is a cloud-based system to support proteomics laboratories in daily quality assessment using a user-friendly interface, easy setup, and automated data processing. Since its release, QCloud has facilitated automated quality control for proteomics experiments in many laboratories. QCloud provides a quick and effortless evaluation of instrument performance that helps to overcome many analytical challenges derived from clinical and translational research. Here we present an improved version of the system, QCloud2. This new version includes enhancements in the scalability and reproducibility of the quality-control pipelines, and it features an improved front end for data visualization, user management, and chart annotation. The QCloud2 system also includes programmatic access and a standalone local version.

Oral microbiome in down syndrome and its implications on oral health

Introduction: The oral cavity harbors an abundant and diverse microbial community (i.e. the microbiome), whose composition and roles in health and disease have been the focus of intense research. Down syndrome (DS) is associated with particular characteristics in the oral cavity, and with a lower incidence of caries and higher incidence of periodontitis and gingivitis compared to control populations. However, the overall composition of the oral microbiome in DS and how it varies with diverse factors like host age or the pH within the mouth are still poorly understood. Methods: Using a Citizen-Science approach in collaboration with DS associations in Spain, we performed 16S rRNA metabarcoding and high-throughput sequencing, combined with culture and proteomics-based identification of fungi to survey the bacterial and fungal oral microbiome in 27 DS persons (age range 7–55) and control samples matched by geographical distribution, age range, and gender.

Results: We found that DS is associated with low salivary pH and less diverse oral microbiomes, which were characterized by lower levels of Alloprevotella, Atopobium, Candidatus Saccharimonas, and higher amounts of Kingella, Staphylococcus, Gemella, Cardiobacterium, Rothia, Actinobacillus, and greater prevalence of Candida.

Conclusion: Altogether, our study provides a first global snapshot of the oral microbiome in DS. Future studies are required to establish whether the observed differences are related to differential pathology in the oral cavity in DS.

Determination of primary microRNA processing in clinical samples by targeted pri-miR-sequencing

MicroRNA expression is important for gene regulation and deregulated microRNA expression is often observed in diseases such as cancer. The processing of primary microRNA transcripts is an important regulatory step in microRNA biogenesis. Due to low expression level and association with chromatin, primary microRNAs are challenging to study in clinical samples where input material is limited. Here, we present a high-sensitivity targeted method to determine processing efficiency of several hundred primary microRNAs from total RNA that requires relatively few RNA sequencing reads. We validate the method using RNA from HeLa cells and show the applicability to clinical samples by analyzing RNA from normal liver and hepatocellular carcinoma. We identify 24 primary microRNAs with significant changes in processing efficiency from normal liver to hepatocellular carcinoma, among those the highly expressed miRNA-122 and miRNA-21, demonstrating that differential processing of primary microRNAs is occurring and could be involved in disease. With our method presented here we provide means to study pri-miRNA processing in disease from clinical samples.

The altered transcriptome of pediatric myelodysplastic syndrome revealed by RNA sequencing

Pediatric myelodysplastic syndrome (PMDS) is a very rare and still poorly characterized disorder. In this work, we identified novel potential targets of PMDS by determining genes with aberrant expression, which can be correlated with PMDS pathogenesis. We identified 291 differentially expressed genes (DEGs) in PMDS patients, comprising genes involved in the regulation of apoptosis and the cell cycle, ribosome biogenesis, inflammation and adaptive immunity. Ten selected DEGs were then validated, confirming the sequencing data. These DEGs will potentially represent new molecular biomarkers and therapeutic targets for PMDS.

The penile microbiota of Black South African men: relationship with human papillomavirus and HIV infection

BACKGROUND

To date, the microbiota of the human penis has been studied mostly in connection with circumcision, HIV risk and female partner bacterial vaginosis (BV). These studies have shown that male circumcision reduces penile anaerobic bacteria, that greater abundance of penile anaerobic bacteria is correlated with increased cytokine levels and greater risk of HIV infection, and that the penile microbiota is an important harbour for BV-associated bacteria. While circumcision has been shown to significantly reduce the risk of acquiring human papillomavirus (HPV) infection, the relationship of the penile microbiota with HPV is still unknown. In this study, we examined the penile microbiota of HPV-infected men as well as the impact of HIV status.

RESULTS

The penile skin microbiota of 238 men from Cape Town (South Africa) were profiled using Illumina sequencing of the V3-V4 hypervariable regions of the 16S rRNA gene. Corynebacterium and Prevotella were found to be the most abundant genera. Six distinct community state types (CSTs) were identified. CST-1, dominated by Corynebacterium, corresponded to less infections with high-risk HPV (HR-HPV) relative to CSTs 2-6. Men in CST-5 had greater relative abundances of Prevotella, Clostridiales, and Porphyromonas and a lower relative abundance of Corynebacterium. Moreover, they were significantly more likely to have HPV or HR-HPV infections than men in CST-1. Using a machine learning approach, we identified greater relative abundances of the anaerobic BV-associated bacteria (Prevotella, Peptinophilus, and Dialister) and lower relative abundance of Corynebacterium in HR-HPV-infected men compared to HR-HPV-uninfected men. No association was observed between HIV and CST, although the penile microbiota of HIV-infected men had greater relative abundances of Staphylococcus compared to HIV-uninfected men.

CONCLUSIONS

We found significant differences in the penile microbiota composition of men with and without HPV and HIV infections. HIV and HR-HPV infections were strongly associated with greater relative abundances of Staphylococcus and BV-associated bacterial taxa (notably Prevotella, Peptinophilus and Dialister), respectively. It is possible that these taxa could increase susceptibility to HIV and HR-HPV acquisition, in addition to creating conditions in which infections persist. Further longitudinal studies are required to establish causal relationships and to determine the extent of the effect.

MasterOfPores: A Workflow for the Analysis of Oxford Nanopore Direct RNA Sequencing Datasets

The direct RNA sequencing platform offered by Oxford Nanopore Technologies allows for direct measurement of RNA molecules without the need of conversion to complementary DNA, fragmentation or amplification. As such, it is virtually capable of detecting any given RNA modification present in the molecule that is being sequenced, as well as provide polyA tail length estimations at the level of individual RNA molecules. Although this technology has been publicly available since 2017, the complexity of the raw Nanopore data, together with the lack of systematic and reproducible pipelines, have greatly hindered the access of this technology to the general user. Here we address this problem by providing a fully benchmarked workflow for the analysis of direct RNA sequencing reads, termed MasterOfPores. The pipeline starts with a pre-processing module, which converts raw current intensities into multiple types of processed data including FASTQ and BAM, providing metrics of the quality of the run, quality-filtering, demultiplexing, base-calling and mapping. In a second step, the pipeline performs downstream analyses of the mapped reads, including prediction of RNA modifications and estimation of polyA tail lengths. Four direct RNA MinION sequencing runs can be fully processed and analyzed in 10 h on 100 CPUs. The pipeline can also be executed in GPU locally or in the cloud, decreasing the run time fourfold. The software is written using the NextFlow framework for parallelization and portability, and relies on Linux containers such as Docker and Singularity for achieving better reproducibility. The MasterOfPores workflow can be executed on any Unix-compatible OS on a computer, cluster or cloud without the need of installing any additional software or dependencies, and is freely available in Github (https://github.com/biocorecrg/master_of_pores). This workflow simplifies direct RNA sequencing data analyses, facilitating the study of the (epi)transcriptome at single molecule resolution.

Fine tuning the extracellular environment accelerates the derivation of kidney organoids from human pluripotent stem cells

The generation of organoids is one of the biggest scientific advances in regenerative medicine. Here, by lengthening the time that human pluripotent stem cells (hPSCs) were exposed to a three-dimensional microenvironment, and by applying defined renal inductive signals, we generated kidney organoids that transcriptomically matched second-trimester human fetal kidneys. We validated these results using ex vivo and in vitro assays that model renal development. Furthermore, we developed a transplantation method that utilizes the chick chorioallantoic membrane. This approach created a soft in vivo microenvironment that promoted the growth and differentiation of implanted kidney organoids, as well as providing a vascular component. The stiffness of the in ovo chorioallantoic membrane microenvironment was recapitulated in vitro by fabricating compliant hydrogels. These biomaterials promoted the efficient generation of renal vesicles and nephron structures, demonstrating that a soft environment accelerates the differentiation of hPSC-derived kidney organoids.

Identification of a Novel Transcription Factor Required for Osteogenic Differentiation of Mesenchymal Stem Cells

Osteogenic differentiation is a complex and still poorly understood biological process regulated by intrinsic cellular signals and extrinsic microenvironmental cues. Following appropriate stimuli, mesenchymal stem cells (MSCs) differentiate into osteoblasts through a tightly regulated multistep process driven by several transcription factors and characterized by the expression of a number of bone-specific proteins. In this study, we describe a novel transcription factor that we named osteoblast inducer (ObI)-1, involved in MSC differentiation toward the osteogenic lineage. ObI-1 encodes for a nuclear protein subjected to proteasomal degradation and expressed during osteoblast differentiation both in a murine multipotent mesenchymal cell line (W20-17) and in primary murine MSCs. RNA interference-mediated knockdown of ObI-1 expression significantly impairs osteoblast differentiation and matrix mineralization with reduced expression of the osteogenic markers, Runt-related transcription factor 2 (Runx2) and osteopontin. Conversely, ObI-1 overexpression enhances osteogenic differentiation and bone-specific markers expression. ObI-1 stimulates bone morphogenetic protein (BMP)-4 expression and the consequent activation of the Smad pathway; treatment with a BMP receptor type I antagonist completely abolishes ObI-1-mediated stimulation of osteogenic differentiation. Collectively, our findings suggest that ObI-1 modulates osteogenic differentiation, at least in part, through the BMP signaling pathway, increasing Runx2 activation and leading to osteoblast commitment and maturation.

Citizen science charts two major "stomatotypes" in the oral microbiome of adolescents and reveals links with habits and drinking water composition

BACKGROUND

The oral cavity comprises a rich and diverse microbiome, which plays important roles in health and disease. Previous studies have mostly focused on adult populations or in very young children, whereas the adolescent oral microbiome remains poorly studied. Here, we used a citizen science approach and 16S profiling to assess the oral microbiome of 1500 adolescents around Spain and its relationships with lifestyle, diet, hygiene, and socioeconomic and environmental parameters.

RESULTS

Our results provide a detailed snapshot of the adolescent oral microbiome and how it varies with lifestyle and other factors. In addition to hygiene and dietary habits, we found that the composition of tap water was related to important changes in the abundance of several bacterial genera. This points to an important role of drinking water in shaping the oral microbiota, which has been so far poorly explored. Overall, the microbiome samples of our study can be clustered into two broad compositional patterns (stomatotypes), driven mostly by Neisseria and Prevotella, respectively. These patterns show striking similarities with those found in unrelated populations.

CONCLUSIONS

We hypothesize that these stomatotypes represent two possible global optimal equilibria in the oral microbiome that reflect underlying constraints of the human oral niche. As such, they should be found across a variety of geographical regions, lifestyles, and ages.

An atlas of alternative splicing profiles and functional associations reveals new regulatory programs and genes that simultaneously express multiple major isoforms

Alternative splicing (AS) generates remarkable regulatory and proteomic complexity in metazoans. However, the functions of most AS events are not known, and programs of regulated splicing remain to be identified. To address these challenges, we describe the Vertebrate Alternative Splicing and Transcription Database (VastDB), the largest resource of genome-wide, quantitative profiles of AS events assembled to date. VastDB provides readily accessible quantitative information on the inclusion levels and functional associations of AS events detected in RNA-seq data from diverse vertebrate cell and tissue types, as well as developmental stages. The VastDB profiles reveal extensive new intergenic and intragenic regulatory relationships among different classes of AS and previously unknown and conserved landscapes of tissue-regulated exons. Contrary to recent reports concluding that nearly all human genes express a single major isoform, VastDB provides evidence that at least 48% of multiexonic protein-coding genes express multiple splice variants that are highly regulated in a cell/tissue-specific manner, and that >18% of genes simultaneously express multiple major isoforms across diverse cell and tissue types. Isoforms encoded by the latter set of genes are generally coexpressed in the same cells and are often engaged by translating ribosomes. Moreover, they are encoded by genes that are significantly enriched in functions associated with transcriptional control, implying they may have an important and wide-ranging role in controlling cellular activities. VastDB thus provides an unprecedented resource for investigations of AS function and regulation.

A De Novo Nonsense Mutation in MAGEL2 in a Patient Initially Diagnosed as Opitz-C: Similarities Between Schaaf-Yang and Opitz-C Syndromes

Opitz trigonocephaly C syndrome (OTCS) is a rare genetic disorder characterized by craniofacial anomalies, variable intellectual and psychomotor disability, and variable cardiac defects with a high mortality rate. Different patterns of inheritance and genetic heterogeneity are known in this syndrome. Whole exome and genome sequencing of a 19-year-old girl (P7), initially diagnosed with OTCS, revealed a de novo nonsense mutation, p.Q638*, in the MAGEL2 gene. MAGEL2 is an imprinted, maternally silenced, gene located at 15q11-13, within the Prader-Willi region. Patient P7 carried the mutation in the paternal chromosome. Recently, mutations in MAGEL2 have been described in Schaaf-Yang syndrome (SHFYNG) and in severe arthrogryposis. Patient P7 bears resemblances with SHFYNG cases but has other findings not described in this syndrome and common in OTCS. We sequenced MAGEL2 in nine additional OTCS patients and no mutations were found. This study provides the first clear molecular genetic basis for an OTCS case, indicates that there is overlap between OTCS and SHFYNG syndromes, and confirms that OTCS is genetically heterogeneous. Genes encoding MAGEL2 partners, either in the retrograde transport or in the ubiquitination-deubiquitination complexes, are promising candidates as OTCS disease-causing genes.

The complete 12 Mb genome and transcriptome of Nonomuraea gerenzanensis with new insights into its duplicated "magic" RNA polymerase

In contrast to the widely accepted consensus of the existence of a single RNA polymerase in bacteria, several actinomycetes have been recently shown to possess two forms of RNA polymerases due the to co-existence of two rpoB paralogs in their genome. However, the biological significance of the rpoB duplication is obscure. In this study we have determined the genome sequence of the lipoglycopeptide antibiotic A40926 producer Nonomuraea gerenzanensis ATCC 39727, an actinomycete with a large genome and two rpoB genes, i.e. rpoB(S) (the wild-type gene) and rpoB(R) (the mutant-type gene). We next analyzed the transcriptional and metabolite profiles in the wild-type gene and in two derivative strains over-expressing either rpoB(R) or a mutated form of this gene to explore the physiological role and biotechnological potential of the "mutant-type" RNA polymerase. We show that rpoB(R) controls antibiotic production and a wide range of metabolic adaptive behaviors in response to environmental pH. This may give interesting perspectives also with regard to biotechnological applications.

Extreme genomic erosion after recurrent demographic bottlenecks in the highly endangered Iberian lynx

BACKGROUND

Genomic studies of endangered species provide insights into their evolution and demographic history, reveal patterns of genomic erosion that might limit their viability, and offer tools for their effective conservation. The Iberian lynx (Lynx pardinus) is the most endangered felid and a unique example of a species on the brink of extinction.

RESULTS

We generate the first annotated draft of the Iberian lynx genome and carry out genome-based analyses of lynx demography, evolution, and population genetics. We identify a series of severe population bottlenecks in the history of the Iberian lynx that predate its known demographic decline during the 20th century and have greatly impacted its genome evolution. We observe drastically reduced rates of weak-to-strong substitutions associated with GC-biased gene conversion and increased rates of fixation of transposable elements. We also find multiple signatures of genetic erosion in the two remnant Iberian lynx populations, including a high frequency of potentially deleterious variants and substitutions, as well as the lowest genome-wide genetic diversity reported so far in any species.

CONCLUSIONS

The genomic features observed in the Iberian lynx genome may hamper short- and long-term viability through reduced fitness and adaptive potential. The knowledge and resources developed in this study will boost the research on felid evolution and conservation genomics and will benefit the ongoing conservation and management of this emblematic species.

PHF13 is a molecular reader and transcriptional co-regulator of H3K4me2/3

PHF13 is a chromatin affiliated protein with a functional role in differentiation, cell division, DNA damage response and higher chromatin order. To gain insight into PHF13's ability to modulate these processes, we elucidate the mechanisms targeting PHF13 to chromatin, its genome wide localization and its molecular chromatin context. Size exclusion chromatography, mass spectrometry, X-ray crystallography and ChIP sequencing demonstrate that PHF13 binds chromatin in a multivalent fashion via direct interactions with H3K4me2/3 and DNA, and indirectly via interactions with PRC2 and RNA PolII. Furthermore, PHF13 depletion disrupted the interactions between PRC2, RNA PolII S5P, H3K4me3 and H3K27me3 and resulted in the up and down regulation of genes functionally enriched in transcriptional regulation, DNA binding, cell cycle, differentiation and chromatin organization. Together our findings argue that PHF13 is an H3K4me2/3 molecular reader and transcriptional co-regulator, affording it the ability to impact different chromatin processes.

DOI:http://dx.doi.org/10.7554/eLife.10607.001

In human and other eukaryotic cells, DNA is packaged around proteins called histones to form a structure known as chromatin. Chemical tags added to the histones alter how the DNA is packaged and the activity of the genes encoded by that DNA. For example, many active genes are packaged around histone H3 proteins that have “Lysine 4 tri-methyl” tags attached to them. Another protein that is associated with chromatin is called PHF13 and it has several roles, including repairing damaged DNA. However, it was not known whether PHF13 binds to chromatin via the chemical tags, or in another way.

Ho-Ryun, Xu, Fuchs et al. used several biochemical techniques in mouse and human cells to explore how PHF13 specifically interacts with chromatin. These experiments showed that PHF13 binds specifically to DNA and to two types of methyl tags (lysine 4-tri-methyl or lysine 4-di-methyl). These chemical tags are predominantly found at active promoters as well as at a small subset of less active promoters known as bivalent promoters. PHF13 interacted with other proteins on the chromatin that are known to either drive or repress gene activity and it’s depletion affected the activity of many genes. Whether PHF13 increased or decreased gene activity depended on whether it was bound to active or bivalent promoters. The active promoters targeted by PHF13 had higher numbers of the tri-methyl tags whereas the di-methyl tags were more common on the bivalent promoters.

These findings provide preliminary evidence that a protein binding to different methyl tags in the same place on histone H3 can have opposite effects on gene activity. Ho-Ryun, Xu, Fuchs et al. now intend to find out more about the other proteins that interact with PHF13 on chromatin.

DOI:http://dx.doi.org/10.7554/eLife.10607.002

Genome and transcriptome analysis of the Mesoamerican common bean and the role of gene duplications in establishing tissue and temporal specialization of genes

BACKGROUND

Legumes are the third largest family of angiosperms and the second most important crop class. Legume genomes have been shaped by extensive large-scale gene duplications, including an approximately 58 million year old whole genome duplication shared by most crop legumes.

RESULTS

We report the genome and the transcription atlas of coding and non-coding genes of a Mesoamerican genotype of common bean (Phaseolus vulgaris L., BAT93). Using a comprehensive phylogenomics analysis, we assessed the past and recent evolution of common bean, and traced the diversification of patterns of gene expression following duplication. We find that successive rounds of gene duplications in legumes have shaped tissue and developmental expression, leading to increased levels of specialization in larger gene families. We also find that many long non-coding RNAs are preferentially expressed in germ-line-related tissues (pods and seeds), suggesting that they play a significant role in fruit development. Our results also suggest that most bean-specific gene family expansions, including resistance gene clusters, predate the split of the Mesoamerican and Andean gene pools.

CONCLUSIONS

The genome and transcriptome data herein generated for a Mesoamerican genotype represent a counterpart to the genomic resources already available for the Andean gene pool. Altogether, this information will allow the genetic dissection of the characters involved in the domestication and adaptation of the crop, and their further implementation in breeding strategies for this important crop.

Prp4 Kinase Grants the License to Splice: Control of Weak Splice Sites during Spliceosome Activation

The genome of the fission yeast Schizosaccharomyces pombe encodes 17 kinases that are essential for cell growth. These include the cell-cycle regulator Cdc2, as well as several kinases that coordinate cell growth, polarity, and morphogenesis during the cell cycle. In this study, we further characterized another of these essential kinases, Prp4, and showed that the splicing of many introns is dependent on Prp4 kinase activity. For detailed characterization, we chose the genes res1 and ppk8, each of which contains one intron of typical size and position. Splicing of the res1 intron was dependent on Prp4 kinase activity, whereas splicing of the ppk8 intron was not. Extensive mutational analyses of the 5’ splice site of both genes revealed that proper transient interaction with the 5’ end of snRNA U1 governs the dependence of splicing on Prp4 kinase activity. Proper transient interaction between the branch sequence and snRNA U2 was also important. Therefore, the Prp4 kinase is required for recognition and efficient splicing of introns displaying weak exon1/5’ splice sites and weak branch sequences.

Prp4 is an essential protein kinase that is involved in the splicing of some introns. Using a conditional mutant of Prp4, we showed that a subset of genes, including several cell cycle–regulatory genes, are dependent on Prp4 for splicing. Furthermore, we could convert genes between Prp4-dependent and -independent states by introducing single-nucleotide mutations in the exon1/5’ splice sites and branch sequence of introns. This work shows that Prp4 activity is required for splicing surveillance in a subset of mRNAs.

Comparative "-omics" in Mycoplasma pneumoniae Clinical Isolates Reveals Key Virulence Factors

The human respiratory tract pathogen M. pneumoniae is one of the best characterized minimal bacterium. Until now, two main groups of clinical isolates of this bacterium have been described (types 1 and 2), differing in the sequence of the P1 adhesin gene. Here, we have sequenced the genomes of 23 clinical isolates of M. pneumoniae. Studying SNPs, non-synonymous mutations, indels and genome rearrangements of these 23 strains and 4 previously sequenced ones, has revealed new subclasses in the two main groups, some of them being associated with the country of isolation. Integrative analysis of in vitro gene essentiality and mutation rates enabled the identification of several putative virulence factors and antigenic proteins; revealing recombination machinery, glycerol metabolism and peroxide production as possible factors in the genetics and physiology of these pathogenic strains. Additionally, the transcriptomes and proteomes of two representative strains, one from each of the two main groups, have been characterized to evaluate the impact of mutations on RNA and proteins levels. This study has revealed that type 2 strains show higher expression levels of CARDS toxin, a protein recently shown to be one of the major factors of inflammation. Thus, we propose that type 2 strains could be more toxigenic than type 1 strains of M. pneumoniae.

MyMpn: a database for the systems biology model organism Mycoplasma pneumoniae

MyMpn (http://mympn.crg.eu) is an online resource devoted to studying the human pathogen Mycoplasma pneumoniae, a minimal bacterium causing lower respiratory tract infections. Due to its small size, its ability to grow in vitro, and the amount of data produced over the past decades, M. pneumoniae is an interesting model organisms for the development of systems biology approaches for unicellular organisms. Our database hosts a wealth of omics-scale datasets generated by hundreds of experimental and computational analyses. These include data obtained from gene expression profiling experiments, gene essentiality studies, protein abundance profiling, protein complex analysis, metabolic reactions and network modeling, cell growth experiments, comparative genomics and 3D tomography. In addition, the intuitive web interface provides access to several visualization and analysis tools as well as to different data search options. The availability and—even more relevant—the accessibility of properly structured and organized data are of up-most importance when aiming to understand the biology of an organism on a global scale. Therefore, MyMpn constitutes a unique and valuable new resource for the large systems biology and microbiology community.

Complete sequencing of Novosphingobium sp. PP1Y reveals a biotechnologically meaningful metabolic pattern

Background

Novosphingobium sp. strain PP1Y is a marine α-proteobacterium adapted to grow at the water/fuel oil interface. It exploits the aromatic fraction of fuel oils as a carbon and energy source. PP1Y is able to grow on a wide range of mono-, poly- and heterocyclic aromatic hydrocarbons. Here, we report the complete functional annotation of the whole Novosphingobium genome.

Results

PP1Y genome analysis and its comparison with other Sphingomonadal genomes has yielded novel insights into the molecular basis of PP1Y’s phenotypic traits, such as its peculiar ability to encapsulate and degrade the aromatic fraction of fuel oils. In particular, we have identified and dissected several highly specialized metabolic pathways involved in: (i) aromatic hydrocarbon degradation; (ii) resistance to toxic compounds; and (iii) the quorum sensing mechanism.

Conclusions

In summary, the unraveling of the entire PP1Y genome sequence has provided important insight into PP1Y metabolism and, most importantly, has opened new perspectives about the possibility of its manipulation for bioremediation purposes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-384) contains supplementary material, which is available to authorized users.

Zrf1 is required to establish and maintain neural progenitor identity

The molecular mechanisms underlying specification from embryonic stem cells (ESCs) and maintenance of neural progenitor cells (NPCs) are largely unknown. Recently, we reported that the Zuotin-related factor 1 (Zrf1) is necessary for chromatin displacement of the Polycomb-repressive complex 1 (PRC1). We found that Zrf1 is required for NPC specification from ESCs and that it promotes the expression of NPC markers, including the key regulator Pax6. Moreover, Zrf1 is essential to establish and maintain Wnt ligand expression levels, which are necessary for NPC self-renewal. Reactivation of proper Wnt signaling in Zrf1-depleted NPCs restores Pax6 expression and the self-renewal capacity. ESC-derived NPCs in vitro resemble most of the characteristics of the self-renewing NPCs located in the developing embryonic cortex, which are termed radial glial cells (RGCs). Depletion of Zrf1 in vivo impairs the expression of key self-renewal regulators and Wnt ligand genes in RGCs. Thus, we demonstrate that Zrf1 plays an essential role in NPC generation and maintenance.

Comparative transcriptomics of early dipteran development

Background

Modern sequencing technologies have massively increased the amount of data available for comparative genomics. Whole-transcriptome shotgun sequencing (RNA-seq) provides a powerful basis for comparative studies. In particular, this approach holds great promise for emerging model species in fields such as evolutionary developmental biology (evo-devo).

Results

We have sequenced early embryonic transcriptomes of two non-drosophilid dipteran species: the moth midge Clogmia albipunctata, and the scuttle fly Megaselia abdita. Our analysis includes a third, published, transcriptome for the hoverfly Episyrphus balteatus. These emerging models for comparative developmental studies close an important phylogenetic gap between Drosophila melanogaster and other insect model systems. In this paper, we provide a comparative analysis of early embryonic transcriptomes across species, and use our data for a phylogenomic re-evaluation of dipteran phylogenetic relationships.

Conclusions

We show how comparative transcriptomics can be used to create useful resources for evo-devo, and to investigate phylogenetic relationships. Our results demonstrate that de novo assembly of short (Illumina) reads yields high-quality, high-coverage transcriptomic data sets. We use these data to investigate deep dipteran phylogenetic relationships. Our results, based on a concatenation of 160 orthologous genes, provide support for the traditional view of Clogmia being the sister group of Brachycera (Megaselia, Episyrphus, Drosophila), rather than that of Culicomorpha (which includes mosquitoes and blackflies).

RYBP and Cbx7 define specific biological functions of polycomb complexes in mouse embryonic stem cells

The Polycomb repressive complex 1 (PRC1) is required for decisions of stem cell fate. In mouse embryonic stem cells (ESCs), two major variations of PRC1 complex, defined by the mutually exclusive presence of Cbx7 or RYBP, have been identified. Here, we show that although the genomic localization of the Cbx7- and RYBP-containing PRC1 complexes overlaps in certain genes, it can also be mutually exclusive. At the molecular level, Cbx7 is necessary for recruitment of Ring1B to chromatin, whereas RYBP enhances the PRC1 enzymatic activity. Genes occupied by RYBP show lower levels of Ring1B and H2AK119ub and are consequently more highly transcribed than those bound by Cbx7. At the functional level, we show that genes occupied by RYBP are primarily involved in the regulation of metabolism and cell-cycle progression, whereas those bound by Cbx7 predominantly control early-lineage commitment of ESCs. Altogether, our results indicate that different PRC1 subtypes establish a complex pattern of gene regulation that regulates common and nonoverlapping aspects of ESC pluripotency and differentiation.

Age-associated inflammation inhibits epidermal stem cell function

Altered stem cell homeostasis is linked to organismal aging. However, the mechanisms involved remain poorly understood. Here we report novel alterations in hair follicle stem cells during skin aging, including increased numbers, decreased function, and an inability to tolerate stress. Performing high-throughput RNA sequencing on aging stem cells, cytokine arrays, and functional assays, we identify an age-associated imbalance in epidermal Jak-Stat signaling that inhibits stem cell function. Collectively, this study reveals a role for the aging epidermis in the disruption of cytokine and stem cell homeostasis, suggesting that stem cell decline during aging may be part of broader tumor-suppressive mechanisms.

Direct targets of Klf5 transcription factor contribute to the maintenance of mouse embryonic stem cell undifferentiated state

Background

A growing body of evidence has shown that Krüppel-like transcription factors play a crucial role in maintaining embryonic stem cell (ESC) pluripotency and in governing ESC fate decisions. Krüppel-like factor 5 (Klf5) appears to play a critical role in these processes, but detailed knowledge of the molecular mechanisms of this function is still not completely addressed.

Results

By combining genome-wide chromatin immunoprecipitation and microarray analysis, we have identified 161 putative primary targets of Klf5 in ESCs. We address three main points: (1) the relevance of the pathways governed by Klf5, demonstrating that suppression or constitutive expression of single Klf5 targets robustly affect the ESC undifferentiated phenotype; (2) the specificity of Klf5 compared to factors belonging to the same family, demonstrating that many Klf5 targets are not regulated by Klf2 and Klf4; and (3) the specificity of Klf5 function in ESCs, demonstrated by the significant differences between Klf5 targets in ESCs compared to adult cells, such as keratinocytes.

Conclusions

Taken together, these results, through the definition of a detailed list of Klf5 transcriptional targets in mouse ESCs, support the important and specific functional role of Klf5 in the maintenance of the undifferentiated ESC phenotype.

See: http://www.biomedcental.com/1741-7007/8/125

miRNA 34a, 100, and 137 modulate differentiation of mouse embryonic stem cells

MicroRNAs (miRNAs) play an important role in proper function and differentiation of mouse embryonic stem cells (ESCs). We performed a systematic comparison of miRNA expression in undifferentiated vs. differentiating ESCs. We report that 138 miRNAs are increased on the induction of differentiation. We compared the entire list of candidate mRNA targets of up-regulated miRNAs with that of mRNA down-regulated in ESCs on induction of differentiation. Among the candidate targets emerging from this analysis, we found three genes, Smarca5, Jarid1b, and Sirt1, previously demonstrated to be involved in sustaining the undifferentiated phenotype in ESCs. On this basis, we first demonstrated that Smarca5 is a direct target of miR-100, Jarid1b of miR-137, and we also confirmed previously published data demonstrating that Sirt1 is a direct target of miR-34a in a different context. The suppression of these three miRNAs by anti-miRs caused the block of ESC differentiation induced by LIF withdrawal. On the other hand, the overexpression of the three miRNAs resulted in an altered expression of differentiation markers. These results demonstrate that miR-34a, miR-100, and miR-137 are required for proper differentiation of mouse ESCs, and that they function in part by targeting Sirt1, Smarca5, and Jarid1b mRNAs.

Systematic identification of stem-loop containing sequence families in bacterial genomes

BACKGROUND

Analysis of non-coding sequences in several bacterial genomes brought to the identification of families of repeated sequences, able to fold as secondary structures. These sequences have often been claimed to be transcribed and fulfill a functional role. A previous systematic analysis of a representative set of 40 bacterial genomes produced a large collection of sequences, potentially able to fold as stem-loop structures (SLS). Computational analysis of these sequences was carried out by searching for families of repetitive nucleic acid elements sharing a common secondary structure.

RESULTS

The initial clustering procedure identified clusters of similar sequences in 29 genomes, corresponding to about 1% of the whole population. Sequences selected in this way have a substantially higher aptitude to fold into a stable secondary structure than the initial set. Removal of redundancies and regrouping of the selected sequences resulted in a final set of 92 families, defined by HMM analysis. 25 of them include all well-known SLS containing repeats and others reported in literature, but not analyzed in detail. The remaining 67 families have not been previously described. Two thirds of the families share a common predicted secondary structure and are located within intergenic regions.

CONCLUSION

Systematic analysis of 40 bacterial genomes revealed a large number of repeated sequence families, including known and novel ones. Their predicted structure and genomic location suggest that, even in compact bacterial genomes, a relatively large fraction of the genome consists of non-protein-coding sequences, possibly functioning at the RNA level.

Effects of a 8-oxoadenosine incorporation on quadruplex structures: thermal stabilities and structural studies

The effects of incorporation of 8-oxoadenosine in two different truncations of human telomeric sequence forming quadruplex structures are reported. In order to characterise their structures, a combination of NMR and UV spectroscopy and computational techniques were used. Both oligonucleotides have been found to form fourfold symmetric quadruplex structures. As a tautomeric equilibrium between keto and enol forms of 8-oxoadenosine may establish in solution and intrinsic stabilities effects, such as internal H-bonds, for example, may determine the predominance of some particular tautomer, molecular modelling studies were performed on quadruplex structures containing both the tautomeric forms. Both molecules resulted to be thermally less stable than the natural.