RESUMEN

BACKGROUND: Rhizobium giardinii has been demonstrated to colonize the roots of a variety of legume species, including common beans, and to increase nitrogen fixation. This suggests that Rhizobium giardinii might be a beneficial tool for sustainable agriculture by lowering dependency on synthetic nitrogen fertilizers and enhancing soil fertility. Understanding the regulatory components in the R. giardinii A3AY_RS01 genes might also lead to the creation of innovative ways for increasing the effectiveness of nitrogen fixation in other agriculturally important bacteria. Therefore, this study was aimed to predict regulatory element of R. giardinii DNA-binding response regulator A3AY_RS01 genes. RESULTS: The locations for 19 % of the Transcriptional start site (TSSs) were within -300 bp relative to the start codon and ten candidate motifs were identified that are shared by at least 50 % of the R. giardinii A3AY_RS01 promoter input sequences from both strands. Motif 1 was revealed as the common promoter motif for all of R. giardinii A3AY_RS01 genes that serves as binding sites for TFs involved in the expression regulation of these genes. Hence, it was revealed that Motif 1 may serve as the binding site chiefly for Ferric uptake regulator (Fur) transcription factor family to regulate expression of A3AY_RS01 genes. High CpG density in the promoter than body regions were observed for most of the genes except for A3AY_RS0102950, A3AY_RS0120195 and A3AY_RS0131150 genes. Nonetheless, promoter areas were richer than body regions in both techniques. CONCLUSIONS: MV1 motif can serve as a binding site for the Fur transcription factor gene family in R. giardinii to regulate the expression of R. giardinii A3AY_RS01 genes. R. giardinii A3AY_RS01 genes are rich in CpG Islands, and play an important role in the regulation of the gene expression of nitrogen fixing in this bacterium.

RESUMEN

Long-read RNA sequencing has shed light on transcriptomic complexity, but questions remain about the functionality of downstream protein products. We introduce Biosurfer, a computational approach for comparing protein isoforms, while systematically tracking the transcriptional, splicing, and translational variations that underlie differences in the sequences of the protein products. Using Biosurfer, we analyzed the differences in 32,799 pairs of GENCODE annotated protein isoforms, finding a majority (70%) of variable N-termini are due to the alternative transcription start sites, while only 9% arise from 5' UTR alternative splicing. Biosurfer's detailed tracking of nucleotide-to-residue relationships helped reveal an uncommonly tracked source of single amino acid residue changes arising from the codon splits at junctions. For 17% of internal sequence changes, such split codon patterns lead to single residue differences, termed "ragged codons". Of variable C-termini, 72% involve splice- or intron retention-induced reading frameshifts. We found an unusual pattern of reading frame changes, in which the first frameshift is closely followed by a distinct second frameshift that restores the original frame, which we term a "snapback" frameshift. We analyzed long read RNA-seq-predicted proteome of a human cell line and found similar trends as compared to our GENCODE analysis, with the exception of a higher proportion of isoforms predicted to undergo nonsense-mediated decay. Biosurfer's comprehensive characterization of long-read RNA-seq datasets should accelerate insights of the functional role of protein isoforms, providing mechanistic explanation of the origins of the proteomic diversity driven by the alternative splicing. Biosurfer is available as a Python package at https://github.com/sheynkman-lab/biosurfer.

RESUMEN

BACKGROUND: Genome-wide landscape of alternative promoter use remains unknown. We determined expression profiles of promoters in 26 lung adenocarcinoma cell lines using the transcriptional start site-sequencing data and proposed an index 'canonical promoter usage' to quantify the diversity of alternative promoter usage. METHODS: Transcriptional start site-sequencing and other datasets were obtained from the DataBase of Transcriptional Start Sites. Transcriptional start site-sequencing read clusters were mapped onto RefGene to determine the promoters. Commonly used promoters were designated as canonical promoters. The sequence logos, CpG islands, DNA methylation and histone modifications of canonical and non-canonical promoters were examined. Canonical promoter usage was calculated by dividing 'read counts of a canonical promoter' by 'read counts of all the units of promoters' on each gene. The expressed genes were subjected to hierarchical clustering according to their canonical promoter usage. RESULTS: Among 104 455 promoters for 14 297 genes, 8659 canonical and 68 197 non-canonical promoters were identified. Corresponding to higher expression, canonical promoters showed core promoter sequences, higher CpG island positivity, less DNA methylation and higher transcription-promoting histone modifications. Gene ontology enrichment analysis revealed that the clusters with lower canonical promoter usage were related to signalling pathways, whereas clusters of tightly regulated genes with higher canonical promoter usage were related to housekeeping genes. CONCLUSION: Canonical promoters were regulated by conventional transcriptional machinery, while non-canonical promoters would be targets of 'leaky' expression. Further investigation is warranted to analyse the correlation between alternative promoter usage and biological characteristics contributing to carcinogenesis.

Asunto(s)

RESUMEN

OBJECTIVE: To establish a quantitative method to evaluate the DNA methylation level of an immediate upstream region of major BRCA1 transcriptional start sites (TSSs), and to investigate whether methylation of the region is a prognostic factor in high-grade serous ovarian cancer patients after neoadjuvant chemotherapy. METHODS: Ninety-two FFPE samples of advanced high-grade serous ovarian cancers after neoadjuvant chemotherapy between 2011 and 2018 were used for mutation and methylation analysis. DNA methylation levels were assessed by pyrosequencing and DNA methylation microarray. An association between methylation level (or a mutation) and progression-free survival was assessed by Kaplan-Meier analysis. RESULT: Major BRCA1 transcripts and CpG sites immediately upstream of their TSSs were identified, and a pyrosequencing method was developed. BRCA1 methylation, BRCA1/2 mutations, and a RAD51C mutation were detected in 17/79 (21.5%), 17/92 (18.5%), and 1/92 (1.1%) high-grade serious ovarian cancer samples. In univariate analysis, BRCA1 methylation and no residual tumor were associated with progression-free survival (BRCA1 methylation: P = 0.025, no residual tumor: P = 0.0026). Multivariate analysis showed that both BRCA1 methylation (P = 0.038, HR = 0.47, 95% CI: 0.21-0.96) and no residual tumor (P = 0.012, HR = 0.49, 95% CI: 0.28-0.85) were significant favorable prognostic factors. CONCLUSION: A quantitative method to estimate the methylation level of the immediate upstream region of major BRCA1 TSSs was established. Methylation of the region of was an independent favorable prognostic factor in high-grade serous ovarian cancer patients.

Asunto(s)

RESUMEN

Dichloromethane (DCM, methylene chloride) is a toxic halogenated volatile organic compound massively used for industrial applications, and consequently often detected in the environment as a major pollutant. DCM biotransformation suggests a sustainable decontamination strategy of polluted sites. Among methylotrophic bacteria able to use DCM as a sole source of carbon and energy for growth, Methylorubrum extorquens DM4 is a longstanding reference strain. Here, the primary 5'-ends of transcripts were obtained using a differential RNA-seq (dRNA-seq) approach to provide the first transcription start site (TSS) genome-wide landscape of a methylotroph using DCM or methanol. In total, 7231 putative TSSs were annotated and classified with respect to their localization to coding sequences (CDSs). TSSs on the opposite strand of CDS (antisense TSS) account for 31% of all identified TSSs. One-third of the detected TSSs were located at a distance to the start codon inferior to 250 nt (average of 84 nt) with 7% of leaderless mRNA. Taken together, the global TSS map for bacterial growth using DCM or methanol will facilitate future studies in which transcriptional regulation is crucial, and efficient DCM removal at polluted sites is limited by regulatory processes.

RESUMEN

Dysferlin (DYSF) has drawn much attention due to its involvement in dysferlinopathy and was reported to affect monocyte functions in recent studies. However, the role of DYSF in the pathogenesis of atherosclerotic cardiovascular diseases (ASCVD) and the regulation mechanism of DYSF expression have not been fully studied. In this study, Gene Expression Omnibus (GEO) database and epigenome-wide association study (EWAS) literatures were searched to find the DNA methylation-driven genes (including DYSF) of ASCVD. The hub genes related to DYSF were also identified through weighted correlation network analysis (WGCNA). Regulation of DYSF expression through its promoter methylation status was verified using peripheral blood leucocytes (PBLs) from ASCVD patients and normal controls, and experiments on THP1 cells and Apoe-/- mice. Similarly, the expressions of DYSF related hub genes, mainly contained SELL, STAT3 and TMX1, were also validated. DYSF functions were then evaluated by phagocytosis, transwell and adhesion assays in DYSF knock-down and overexpressed THP1 cells. The results showed that DYSF promoter hypermethylation up-regulated its expression in clinical samples, THP1 cells and Apoe-/- mice, confirming DYSF as a DNA methylation-driven gene. The combination of DYSF expression and methylation status in PBLs had a considerable prediction value for ASCVD. Besides, DYSF could enhance the phagocytosis, migration and adhesion ability of THP1 cells. Among DYSF related hub genes, SELL was proven to be the downstream target of DYSF by wet experiments. In conclusion, DYSF promoter hypermethylation upregulated its expression and promoted monocytes activation, which further participated in the pathogenesis of ASCVD.

Asunto(s)

RESUMEN

The sigma H (σΗ) and sigma E (σE) subunits of Corynebacterium glutamicum RNA polymerase belong to Group 4 of sigma factors, also called extracytoplasmic function (ECF) sigma factors. Genes of the C. glutamicum σΗ regulon that are involved in heat and oxidative stress response have already been defined, whereas the genes of the σE regulon, which is involved in cell surface stress response, have not been explored until now. Using the C. glutamicum RES167 strain and its derivative C. glutamicum ΔcseE with a deletion in the anti-σΕ gene, differential gene expression was analyzed by RNA sequencing. We found 296 upregulated and 398 downregulated genes in C. glutamicum ΔcseE compared to C. glutamicum RES167. To confirm the functional link between σΕ and the corresponding promoters, we tested selected promoters using the in vivo two-plasmid system with gfpuv as a reporter gene and by in vitro transcription. Analyses with RNAP+σΗ and RNAP+σΕ, which were previously shown to recognize similar promoters, proved that the σΗ and σE regulons significantly overlap. The σE-controlled genes were found to be involved for example in protein quality control (dnaK, dnaJ2, clpB, and clpC), the regulation of Clp proteases (clgR), and membrane integrity maintenance. The single-promoter analyses with σΗ and σΕ revealed that there are two groups of promoters: those which are exclusively σΗ-specific, and the other group of promoters, which are σΗ/σE-dependent. No exclusively σE-dependent promoter was detected. We defined the consensus sequences of exclusively σΗ-regulated promotors to be -35 GGAAt and - 10 GTT and σΗ/σE-regulated promoters to be -35 GGAAC and - 10 cGTT. Fifteen genes were found to belong to the σΗ/σΕ regulon. Homology modeling showed that there is a specific interaction between Met170 in σΗ and the nucleotides -31 and - 30 within the non-coding strand (AT or CT) of the σΗ-dependent promoters. In σE, Arg185 was found to interact with the nucleotides GA at the same positions in the σE-dependent promoters.

RESUMEN

Chromodomain-Helicase-DNA-binding protein 8 (CHD8) is a chromatin remodeler that is central to regulation of gene expression pathways during brain development. Many loss-of-function mutations in transcribed regions of the chd8 gene have been identified in autism spectrum disorder patients. Nothing is known about transcription of the human chd8 gene. Defects in expression of this gene could represent another mechanism leading to reduced amount of CHD8. We identify two major promoters for the human chd8 gene, both of which are located many thousand base pairs upstream of the coding region. Each proximal promoter directs a similar transcriptional efficiency in transfected cells. At least two elements within 200bp of the 5'flanking regions of these promoters are important to drive highest transcriptional levels in transient transfection experiments. RNA polymerase II occupancy levels at each promoter are roughly equivalent. Lastly, each promoter directs a dispersed set of start sites in a cultured cell line. This work could provide the framework for future studies to investigate the importance of chd8 gene expression for diseases associated with brain and neuronal development.

Asunto(s)

RESUMEN

Gene transcription in bacteria often starts some nucleotides upstream of the start codon. Identifying the specific Transcriptional Start Site (TSS) is essential for genetic manipulation, as in many cases upstream of the start codon there are sequence elements that are involved in gene expression regulation. Taken into account the classical gene structure, we are able to identify two kinds of transcriptional start site: primary and secondary. A primary transcriptional start site is located some nucleotides upstream of the translational start site, while a secondary transcriptional start site is located within the gene encoding sequence. Here, we present a step by step protocol for genome-wide transcriptional start sites determination by differential RNA-sequencing (dRNA-seq) using the enteric pathogen Shigella flexneri serotype 5a strain M90T as model. However, this method can be employed in any other bacterial species of choice. In the first steps, total RNA is purified from bacterial cultures using the hot phenol method. Ribosomal RNA (rRNA) is specifically depleted via hybridization probes using a commercial kit. A 5'-monophosphate-dependent exonuclease (TEX)-treated RNA library enriched in primary transcripts is then prepared for comparison with a library that has not undergone TEX-treatment, followed by ligation of an RNA linker adaptor of known sequence allowing the determination of TSS with single nucleotide precision. Finally, the RNA is processed for Illumina sequencing library preparation and sequenced as purchased service. TSS are identified by in-house bioinformatic analysis. Our protocol is cost-effective as it minimizes the use of commercial kits and employs freely available software.

RESUMEN

Eukaryotes typically utilize two distinct aminoacyl-tRNA synthetase isoforms, one for cytosolic and one for mitochondrial protein synthesis. However, the genome of budding yeast (Saccharomyces cerevisiae) contains only one cysteinyl-tRNA synthetase gene (YNL247W, also known as CRS1). In this study, we report that CRS1 encodes both cytosolic and mitochondrial isoforms. The 5' complementary DNA end method and GFP reporter gene analyses indicated that yeast CRS1 expression yields two classes of mRNAs through alternative transcription starts: a long mRNA containing a mitochondrial targeting sequence and a short mRNA lacking this targeting sequence. We found that the mitochondrial Crs1 is the product of translation from the first initiation AUG codon on the long mRNA, whereas the cytosolic Crs1 is produced from the second in-frame AUG codon on the short mRNA. Genetic analysis and a ChIP assay revealed that the transcription factor heme activator protein (Hap) complex, which is involved in mitochondrial biogenesis, determines the transcription start sites of the CRS1 gene. We also noted that Hap complex-dependent initiation is regulated according to the needs of mitochondrial energy production. The results of our study indicate energy-dependent initiation of alternative transcription of CRS1 that results in production of two Crs1 isoforms, a finding that suggests Crs1's potential involvement in mitochondrial energy metabolism in yeast.

Asunto(s)

RESUMEN

In recent years, the knowledge generated by decoding the human genome has allowed groundbreaking genetic research to better understand genomic architecture and heritability in healthy and disease states. The vast amount of data generated over time and yet to be generated provides the basis for translational research towards the development of preventive and therapeutic strategies for many conditions. In this special issue, we highlight the discoveries of disease-associated and protective DNA variations in common human diseases and developmental disorders.

Asunto(s)

RESUMEN

Muscle-enriched lamin-interacting protein (Mlip) is an alternatively spliced gene whose splicing specificity is dictated by tissue type. MLIP is most abundantly expressed in brain, cardiac, and skeletal muscle. In the present study, we systematically mapped the transcriptional start and stop sites of murine Mlip Rapid amplification of cDNA ends (RACE) of Mlip transcripts from the brain, heart, and skeletal muscle revealed two transcriptional start sites (TSSs), exon 1a and exon 1b, and only one transcriptional termination site. RT-PCR analysis of the usage of the two identified TSSs revealed that the heart utilizes only exon 1a for MLIP expression, whereas the brain exclusively uses exon 1b TSS. Loss of Mlip exon 1a in mice resulted in a 7-fold increase in the prevalence of centralized nuclei in muscle fibers with the Mlip exon1a-deficient satellite cells on single fibers exhibiting a significant delay in commitment to a MYOD-positive phenotype. Furthermore, we demonstrate that the A-type lamin-binding domain in MLIP is encoded in exon 1a, indicating that MLIP isoforms generated with exon 1b TSS lack the A-type lamin-binding domain. Collectively these findings suggest that Mlip tissue-specific expression and alternative splicing play a critical role in determining MLIP's functions in mice.

Asunto(s)

RESUMEN

Micro (mi)RNAs regulate gene expression in many eukaryotic organisms where they control diverse biological processes. Their biogenesis, from primary transcripts to mature miRNAs, have been extensively characterized in animals and plants, showing distinct differences between these phylogenetically distant groups of organisms. However, comparably little is known about miRNA biogenesis in organisms whose evolutionary position is placed in between plants and animals and/or in unicellular organisms. Here, we investigate miRNA maturation in the unicellular amoeba Dictyostelium discoideum, belonging to Amoebozoa, which branched out after plants but before animals. High-throughput sequencing of small RNAs and poly(A)-selected RNAs demonstrated that the Dicer-like protein DrnB is required, and essentially specific, for global miRNA maturation in D. discoideum. Our RNA-seq data also showed that longer miRNA transcripts, generally preceded by a T-rich putative promoter motif, accumulate in a drnB knock-out strain. For two model miRNAs we defined the transcriptional start sites (TSSs) of primary (pri)-miRNAs and showed that they carry the RNA polymerase II specific m7G-cap. The generation of the 3'-ends of these pri-miRNAs differs, with pri-mir-1177 reading into the downstream gene, and pri-mir-1176 displaying a distinct end. This 3´-end is processed to shorter intermediates, stabilized in DrnB-depleted cells, of which some carry a short oligo(A)-tail. Furthermore, we identified 10 new miRNAs, all DrnB dependent and developmentally regulated. Thus, the miRNA machinery in D. discoideum shares features with both plants and animals, which is in agreement with its evolutionary position and perhaps also an adaptation to its complex lifestyle: unicellular growth and multicellular development.

Asunto(s)

RESUMEN

The non-pathogenic bacterium Mycobacterium smegmatis mc2155 has been widely used as a model organism in mycobacterial research, yet a detailed study about its transcription landscape remains to be established. Here we report the transcriptome, expression profiles and transcriptional structures through growth-phase-dependent RNA sequencing (RNA-seq) as well as other related experiments. We found: (1) 2,139 transcriptional start sites (TSSs) in the genome-wide scale, of which eight samples were randomly selected and further verified by 5'-RACE; (2) 2,233 independent monocistronic or polycistronic mRNAs in the transcriptome within the operon/sub-operon structures which are classified into five groups; (3) 47.50% (1016/2139) genes were transcribed into leaderless mRNAs, with the TSSs of 41.3% (883/2139) mRNAs overlapping with the first base of the annotated start codon. Initial amino acids of MSMEG_4921 and MSMEG_6422 proteins were identified by Edman degradation, indicating the presence of distinctive widespread leaderless features in M. smegmatis mc2155. (4) 150 genes with potentially wrong structural annotation, of which 124 proposed genes have been corrected; (5) eight highly active promoters, with their activities further determined by ß-galactosidase assays. These data integrated the transcriptional landscape to genome information of model organism mc2155 and lay a solid foundation for further works in Mycobacterium.

RESUMEN

A major determinant in the efficiency of ribosome loading onto mRNAs is the 5' TL (transcript leader or 5' UTR). In addition, elements within this region also impact on start site selection demonstrating that it can modulate the protein readout at both quantitative and qualitative levels. With the increasing wealth of data generated by the mining of the mammalian transcriptome, it has become evident that a genes 5' TL is not homogeneous but actually exhibits significant heterogeneity. This arises due to the utilization of alternative promoters, and is further compounded by significant variability with regards to the precise transcriptional start sites of each (not to mention alternative splicing). Consequently, the transcript for a protein coding gene is not a unique mRNA, but in-fact a complexed quasi-species of variants whose composition may respond to the changing physiological environment of the cell. Here we examine the potential impact of these events with regards to the protein readout.

RESUMEN

The SRY-related high-mobility box 9 (SOX9) gene is expressed in many different tissues. To better understand the DNA elements that control tissue-specific expression, we cloned and sequenced a 2.5 kb fragment lying 5' to the zebrafish sox9b gene transcriptional start site. Three regions of this clone contained stable secondary structures that hindered cloning, sequencing, and amplification. This segment and smaller fragmentswere inserted 5' of an EGFP reporter and transgenic fish were raised with the different reporters. Reporter expression was also observed in embryos directly injected with the constructs to transiently express the reporter. Heart expression required only a very short 5' sequence, as a 0.6 kb sox9b fragment produced reporter expression in heart in transgenic zebrafish, and transient experiments showed heart expression from a minimal sox9b promoter region containing a conserved TATA box and an EGR2 element (-74/+29 bp). Reporter expression in transgenic skeletal muscle was consistently lower than in other tissues. Jaw, brain, and notochord expression was strong with the full-length clone, but was dramatically reduced as the size of the fragment driving the reporter decreased from approximately 1.8 to 0.9 kb. The 2.5 kb region 5' of the sox9b contained 7 conserved non-coding elements (CNEs) that included putative hypoxia inducible factor 1α (HIF1α), CAAT box (CCAAT), early growth response protein 2 (EGR2), and core promoter elements. While a synthetic fragment containing all 7 CNEs produced some degree of reporter expression in muscle, jaw, heart and brain, the degree of reporter expression was considerably lower than that produced by the full length clone. These results can account for the tissue-specific expression of sox9b in the developing zebrafish.

Asunto(s)

RESUMEN

Rapid amplification of cDNA ends (RACE) is a technique that was developed to swiftly and efficiently amplify full-length RNA molecules in which the terminal ends have not been characterized. Current usage of this procedure has been more focused on sequencing and characterizing RNA 5' and 3' untranslated regions. Herein is described an adapted RACE protocol to amplify bacterial RNA transcripts.

Asunto(s)

RESUMEN

Heterochromatin protein 1α (HP1α) encoded from the CBX5-gene is an evolutionary conserved protein that binds histone H3 di- or tri-methylated at position lysine 9 (H3K9me2/3), a hallmark for heterochromatin, and has an essential role in forming higher order chromatin structures. HP1α has diverse functions in heterochromatin formation, gene regulation, and mitotic progression, and forms complex networks of gene, RNA, and protein interactions. Emerging evidence has shown that HP1α serves a unique biological role in breast cancer related processes and in particular for epigenetic control mechanisms involved in aberrant cell proliferation and metastasis. However, how HP1α deregulation plays dual mechanistic functions for cancer cell proliferation and metastasis suppression and the underlying cellular mechanisms are not yet comprehensively described. In this paper we provide an overview of the role of HP1α as a new sight of epigenetics in proliferation and metastasis of human breast cancer. This highlights the importance of addressing HP1α in breast cancer diagnostics and therapeutics.

Asunto(s)

RESUMEN

Deep-sequencing technologies applied to RNA have tremendous potential to identify novel transcripts with single-nucleotide resolution. By combining whole-transcript cDNA sequencing (RNA-seq) and genome-wide identification of transcription start sites (dRNA-seq), it is possible to characterize long 5'-untranslated regions potentially endowed with regulatory capacities and to detect premature termination of transcription. This can be used to identify new potential riboswitches. In this chapter, we provide a detailed protocol of the dRNA-seq method based on differential pretreatment of RNAs with tobacco acid pyrophosphatase to differentiate between 5'-ends of primary and processed RNAs. We also give a briefer protocol of the preparation of RNA-seq libraries and of how to go through data bioinformatics analysis and data visualization using genome browsers. This approach is powerful to identify novel riboswitches and to demonstrate the functionality of riboswitches predicted in silico.

Asunto(s)