Chromatin accessibility and microRNA expression in nephron progenitor cells during kidney development.

Mammalian nephrons originate from a population of nephron progenitor cells, and changes in these cells' transcriptomes contribute to the cessation of nephrogenesis, an important determinant of nephron number. To characterize microRNA (miRNA) expression and identify putative cis-regulatory regions, we collected nephron progenitor cells from mouse kidneys at embryonic day 14.5 and postnatal day zero and assayed small RNA expression and transposase-accessible chromatin. We detect expression of 1104 miRNA (114 with expression changes), and 46,374 chromatin accessible regions (2103 with changes in accessibility). Genome-wide, our data highlight processes like cellular differentiation, cell migration, extracellular matrix interactions, and developmental signaling pathways. Furthermore, they identify new candidate cis-regulatory elements for Eya1 and Pax8, both genes with a role in nephron progenitor cell differentiation. Finally, we associate expression-changing miRNAs, including let-7-5p, miR-125b-5p, miR-181a-2-3p, and miR-9-3p, with candidate cis-regulatory elements and target genes. These analyses highlight new putative cis-regulatory loci for miRNA in nephron progenitors.

Unusual sequence characteristics of human chromosome 19 are conserved across 11 nonhuman primates.

BACKGROUND: Human chromosome 19 has many unique characteristics including gene density more than double the genome-wide average and 20 large tandemly clustered gene families. It also has the highest GC content of any chromosome, especially outside gene clusters. The high GC content and concomitant high content of hypermutable CpG sites raises the possibility chromosome 19 exhibits higher levels of nucleotide diversity both within and between species, and may possess greater variation in DNA methylation that regulates gene expression. RESULTS: We examined GC and CpG content of chromosome 19 orthologs across representatives of the primate order. In all 12 primate species with suitable genome assemblies, chromosome 19 orthologs have the highest GC content of any chromosome. CpG dinucleotides and CpG islands are also more prevalent in chromosome 19 orthologs than other chromosomes. GC and CpG content are generally higher outside the gene clusters. Intra-species variation based on SNPs in human common dbSNP, rhesus, crab eating macaque, baboon and marmoset datasets is most prevalent on chromosome 19 and its orthologs. Inter-species comparisons based on phyloP conservation show accelerated nucleotide evolution for chromosome 19 promoter flanking and enhancer regions. These same regulatory regions show the highest CpG density of any chromosome suggesting they possess considerable methylome regulatory potential. CONCLUSIONS: The pattern of high GC and CpG content in chromosome 19 orthologs, particularly outside gene clusters, is present from human to mouse lemur representing 74 million years of primate evolution. Much CpG variation exists both within and between primate species with a portion of this variation occurring in regulatory regions.

Regulatory architecture of cell identity genes and housekeeping genes.

Gene regulation and chromosome architecture are intimately linked. Genes with prominent roles in cell identity are often regulated by clusters of enhancer elements. By contrast, a recent study shows housekeeping genes are often regulated through clustering of promoters. We discuss here new regulatory insights for these two types of genes.

Riboswitches, from cognition to transformation.

Riboswitches are functional RNA elements that regulate gene expression by directly detecting metabolites. Twenty years have passed since it was first discovered, researches on riboswitches are becoming increasingly standardized and refined, which could significantly promote people's cognition of RNA function as well. Here, we focus on some representative orphan riboswitches, enumerate the structural and functional transformation and artificial design of riboswitches including the coupling with ribozymes, hoping to attain a comprehensive understanding of riboswitch research.

Convergent and lineage-specific genomic differences in limb regulatory elements in limbless reptile lineages.

Loss of limbs evolved many times in squamate reptiles. Here we investigated the genomic basis of convergent limb loss in reptiles. We sequenced the genomes of a closely related pair of limbless-limbed gymnophthalmid lizards and performed a comparative genomic analysis including five snakes and the limbless glass lizard. Our analysis of these three independent limbless lineages revealed that signatures of shared sequence or transcription factor binding site divergence in individual limb regulatory elements are generally rare. Instead, shared divergence occurs more often at the level of signaling pathways, involving different regulatory elements associated with the same limb genes (such as Hand2 or Hox) and/or patterning mechanisms (such as Shh signaling). Interestingly, although snakes are known to have mutations in the Shh ZRS limb enhancer, this enhancer lacks relevant mutations in limbless lizards. Thus, different mechanisms could contribute to limb loss, and there are likely multiple evolutionary paths to limblessness in reptiles.

Methods of massive parallel reporter assays for investigation of enhancers.

The correct deployment of genetic programs for development and differentiation relies on finely coordinated regulation of specific gene sets. Genomic regulatory elements play an exceptional role in this process. There are few types of gene regulatory elements, including promoters, enhancers, insulators and silencers. Alterations of gene regulatory elements may cause various pathologies, including cancer, congenital disorders and autoimmune diseases. The development of high-throughput genomic assays has made it possible to significantly accelerate the accumulation of information about the characteristic epigenetic properties of regulatory elements. In combination with high-throughput studies focused on the genome-wide distribution of epigenetic marks, regulatory proteins and the spatial structure of chromatin, this significantly expands the understanding of the principles of epigenetic regulation of genes and allows potential regulatory elements to be searched for in silico. However, common experimental approaches used to study the local characteristics of chromatin have a number of technical limitations that may reduce the reliability of computational identification of genomic regulatory sequences. Taking into account the variability of the functions of epigenetic determinants and complex multicomponent regulation of genomic elements activity, their functional verification is often required. A plethora of methods have been developed to study the functional role of regulatory elements on the genome scale. Common experimental approaches for in silico identification of regulatory elements and their inherent technical limitations will be described. The present review is focused on original high-throughput methods of enhancer activity reporter analysis that are currently used to validate predicted regulatory elements and to perform de novo searches. The methods described allow assessing the functional role of the nucleotide sequence of a regulatory element, to determine its exact boundaries and to assess the influence of the local state of chromatin on the activity of enhancers and gene expression. These approaches have contributed substantially to the understanding of the fundamental principles of gene regulation.

Comprehensive multi-omics integration identifies differentially active enhancers during human brain development with clinical relevance.

BACKGROUND: Non-coding regulatory elements (NCREs), such as enhancers, play a crucial role in gene regulation, and genetic aberrations in NCREs can lead to human disease, including brain disorders. The human brain is a complex organ that is susceptible to numerous disorders; many of these are caused by genetic changes, but a multitude remain currently unexplained. Understanding NCREs acting during brain development has the potential to shed light on previously unrecognized genetic causes of human brain disease. Despite immense community-wide efforts to understand the role of the non-coding genome and NCREs, annotating functional NCREs remains challenging. METHODS: Here we performed an integrative computational analysis of virtually all currently available epigenome data sets related to human fetal brain. RESULTS: Our in-depth analysis unravels 39,709 differentially active enhancers (DAEs) that show dynamic epigenomic rearrangement during early stages of human brain development, indicating likely biological function. Many of these DAEs are linked to clinically relevant genes, and functional validation of selected DAEs in cell models and zebrafish confirms their role in gene regulation. Compared to enhancers without dynamic epigenomic rearrangement, DAEs are subjected to higher sequence constraints in humans, have distinct sequence characteristics and are bound by a distinct transcription factor landscape. DAEs are enriched for GWAS loci for brain-related traits and for genetic variation found in individuals with neurodevelopmental disorders, including autism. CONCLUSION: This compendium of high-confidence enhancers will assist in deciphering the mechanism behind developmental genetics of human brain and will be relevant to uncover missing heritability in human genetic brain disorders.

Parallel genetics of regulatory sequences using scalable genome editing in vivo.

How regulatory sequences control gene expression is fundamental for explaining phenotypes in health and disease. Regulatory elements must ultimately be understood within their genomic environment and development- or tissue-specific contexts. Because this is technically challenging, few regulatory elements have been characterized in vivo. Here, we use inducible Cas9 and multiplexed guide RNAs to create hundreds of mutations in enhancers/promoters and 3' UTRs of 16 genes in C. elegans. Our software crispr-DART analyzes indel mutations in targeted DNA sequencing. We quantify the impact of mutations on expression and fitness by targeted RNA sequencing and DNA sampling. When applying our approach to the lin-41 3' UTR, generating hundreds of mutants, we find that the two adjacent binding sites for the miRNA let-7 can regulate lin-41 expression independently of each other. Finally, we map regulatory genotypes to phenotypic traits for several genes. Our approach enables parallel analysis of regulatory sequences directly in animals.

AML Subtype Is a Major Determinant of the Association between Prognostic Gene Expression Signatures and Their Clinical Significance.

Relapse in acute myeloid leukemia (AML) may result from variable genetic origins or convergence on common biological processes. Exploiting the specificity and sensitivity of regulatory DNA, we analyze patient samples of multiple clinical outcomes covering various AML molecular subtypes. We uncover regulatory variation among patients translating into a transcriptional signature that predicts relapse risk. In addition, we find clusters of coexpressed genes within this signature selectively link to relapse risk in distinct patient subgroups defined by molecular subtype or AML maturation. Analyzing these gene clusters and the AML subtypes separately enhances their prognostic value substantially and provides insight in the mechanisms underlying relapse risk across the distinct patient subgroups. We propose that prognostic gene expression signatures in AML are valid only within patient subgroups and do not transcend these subgroups.

Global donor and acceptor splicing site kinetics in human cells.

RNA splicing is an essential part of eukaryotic gene expression. Although the mechanism of splicing has been extensively studied in vitro, in vivo kinetics for the two-step splicing reaction remain poorly understood. Here, we combine transient transcriptome sequencing (TT-seq) and mathematical modeling to quantify RNA metabolic rates at donor and acceptor splice sites across the human genome. Splicing occurs in the range of minutes and is limited by the speed of RNA polymerase elongation. Splicing kinetics strongly depends on the position and nature of nucleotides flanking splice sites, and on structural interactions between unspliced RNA and small nuclear RNAs in spliceosomal intermediates. Finally, we introduce the 'yield' of splicing as the efficiency of converting unspliced to spliced RNA and show that it is highest for mRNAs and independent of splicing kinetics. These results lead to quantitative models describing how splicing rates and yield are encoded in the human genome.

Proviruses with Long-Term Stable Expression Accumulate in Transcriptionally Active Chromatin Close to the Gene Regulatory Elements: Comparison of ASLV-, HIV- and MLV-Derived Vectors.

Individual groups of retroviruses and retroviral vectors differ in their integration site preference and interaction with the host genome. Hence, immediately after infection genome-wide distribution of integrated proviruses is non-random. During long-term in vitro or persistent in vivo infection, the genomic position and chromatin environment of the provirus affects its transcriptional activity. Thus, a selection of long-term stably expressed proviruses and elimination of proviruses, which have been gradually silenced by epigenetic mechanisms, helps in the identification of genomic compartments permissive for proviral transcription. We compare here the extent and time course of provirus silencing in single cell clones of the K562 human myeloid lymphoblastoma cell line that have been infected with retroviral reporter vectors derived from avian sarcoma/leukosis virus (ASLV), human immunodeficiency virus type 1 (HIV) and murine leukaemia virus (MLV). While MLV proviruses remain transcriptionally active, ASLV proviruses are prone to rapid silencing. The HIV provirus displays gradual silencing only after an extended time period in culture. The analysis of integration sites of long-term stably expressed proviruses shows a strong bias for some genomic features-especially integration close to the transcription start sites of active transcription units. Furthermore, complex analysis of histone modifications enriched at the site of integration points to the accumulation of proviruses of all three groups in gene regulatory segments, particularly close to the enhancer loci. We conclude that the proximity to active regulatory chromatin segments correlates with stable provirus expression in various retroviral species.

Development of a Pseudomonas putida cell-free protein synthesis platform for rapid screening of gene regulatory elements.

Cell-free protein synthesis (CFPS) systems enable the production of protein without the use of living, intact cells. An emerging area of interest is to use CFPS systems to characterize individual elements for genetic programs [e.g. promoters, ribosome binding sites (RBS)]. To enable this research area, robust CFPS systems must be developed from new chassis organisms. One such chassis is the Gram-negative Pseudomonas bacteria, which have been studied extensively for their diverse metabolism with promises in the field of bioremediation and biosynthesis. Here, we report the development and optimization of a high-yielding (198 ± 5.9 µg/ml) batch CFPS system from Pseudomonas putida ATCC 12633. Importantly, both circular and linear DNA templates can be applied directly to the CFPS reaction to program protein synthesis. Therefore, it is possible to prepare hundreds or even thousands of DNA templates without time-consuming cloning work. This opens the possibility to rapidly assess and validate genetic part performance in vitro before performing experiments in cells. To validate the P. putida CFPS system as a platform for prototyping genetic parts, we designed and constructed a library consisting of 15 different RBSs upstream of the reporter protein sfGFP, which covered an order of magnitude range in expression. Looking forward, our P. putida CFPS platform will not only expand the protein synthesis toolkit for synthetic biology but also serve as a platform in expediting the screening and prototyping of gene regulatory elements.

Uncovering drug-responsive regulatory elements.

Nucleotide changes in gene regulatory elements can have a major effect on interindividual differences in drug response. For example, by reviewing all published pharmacogenomic genome-wide association studies, we show here that 96.4% of the associated single nucleotide polymorphisms reside in noncoding regions. We discuss how sequencing technologies are improving our ability to identify drug response-associated regulatory elements genome-wide and to annotate nucleotide variants within them. We highlight specific examples of how nucleotide changes in these elements can affect drug response and illustrate the techniques used to find them and functionally characterize them. Finally, we also discuss challenges in the field of drug-responsive regulatory elements that need to be considered in order to translate these findings into the clinic.

Integrating multi-omic features exploiting Chromosome Conformation Capture data.

The representation, integration, and interpretation of omic data is a complex task, in particular considering the huge amount of information that is daily produced in molecular biology laboratories all around the world. The reason is that sequencing data regarding expression profiles, methylation patterns, and chromatin domains is difficult to harmonize in a systems biology view, since genome browsers only allow coordinate-based representations, discarding functional clusters created by the spatial conformation of the DNA in the nucleus. In this context, recent progresses in high throughput molecular biology techniques and bioinformatics have provided insights into chromatin interactions on a larger scale and offer a formidable support for the interpretation of multi-omic data. In particular, a novel sequencing technique called Chromosome Conformation Capture allows the analysis of the chromosome organization in the cell's natural state. While performed genome wide, this technique is usually called Hi-C. Inspired by service applications such as Google Maps, we developed NuChart, an R package that integrates Hi-C data to describe the chromosomal neighborhood starting from the information about gene positions, with the possibility of mapping on the achieved graphs genomic features such as methylation patterns and histone modifications, along with expression profiles. In this paper we show the importance of the NuChart application for the integration of multi-omic data in a systems biology fashion, with particular interest in cytogenetic applications of these techniques. Moreover, we demonstrate how the integration of multi-omic data can provide useful information in understanding why genes are in certain specific positions inside the nucleus and how epigenetic patterns correlate with their expression.

IntergenicDB: a database for intergenic sequences.

UNLABELLED: A whole genome contains not only coding regions, but also non-coding regions. These are located between the end of a given coding region and the beginning of the following coding region. For this reason, the information about gene regulation process underlies in intergenic regions. There is no easy way to obtain intergenic regions from current available databases. IntergenicDB was developed to integrate data of intergenic regions and their gene related information from NCBI databases. The main goal of INTERGENICDB is to offer friendly database for intergenic sequences of bacterial genomes. AVAILABILITY: http://intergenicdb.bioinfoucs.com/