Genomic Footprinting Analyses from DNase-seq Data to Construct Gene Regulatory Networks.

Chromatin accessibility is directly linked with transcription in eukaryotes. Accessible regions associated with regulatory proteins are highly sensitive to DNase I digestion and are termed DNase I hypersensitive sites (DHSs). DHSs can be identified by DNase I digestion, followed by high-throughput DNA sequencing (DNase-seq). The single-base-pair resolution digestion patterns from DNase-seq allows identifying transcription factor (TF) footprints of local DNA protection that predict TF-DNA binding. The identification of differential footprinting between two conditions allows mapping relevant TF regulatory interactions. Here, we provide step-by-step instructions to build gene regulatory networks from DNase-seq data. Our pipeline includes steps for DHSs calling, identification of differential TF footprints between treatment and control conditions, and construction of gene regulatory networks. Even though the data we used in this example was obtained from Arabidopsis thaliana, the workflow developed in this guide can be adapted to work with DNase-seq data from any organism with a sequenced genome.

RhizoBindingSites, a Database of DNA-Binding Motifs in Nitrogen-Fixing Bacteria Inferred Using a Footprint Discovery Approach.

Basic knowledge of transcriptional regulation is needed to understand the mechanisms governing biological processes, i.e., nitrogen fixation by Rhizobiales bacteria in symbiosis with leguminous plants. The RhizoBindingSites database is a computer-assisted framework providing motif-gene-associated conserved sequences potentially implicated in transcriptional regulation in nine symbiotic species. A dyad analysis algorithm was used to deduce motifs in the upstream regulatory region of orthologous genes, and only motifs also located in the gene seed promoter with a p-value of 1e-4 were accepted. A genomic scan analysis of the upstoream sequences with these motifs was performed. These predicted binding sites were categorized according to low, medium and high homology between the matrix and the upstream regulatory sequence. On average, 62.7% of the genes had a motif, accounting for 80.44% of the genes per genome, with 19613 matrices (a matrix is a representation of a motif). The RhizoBindingSites database provides motif and gene information, motif conservation in the order Rhizobiales, matrices, motif logos, regulatory networks constructed from theoretical or experimental data, a criterion for selecting motifs and a guide for users. The RhizoBindingSites database is freely available online at rhizobindingsites.ccg.unam.mx.

Local Changes in Chromatin Accessibility and Transcriptional Networks Underlying the Nitrate Response in Arabidopsis Roots.

Transcriptional regulation, determined by the chromatin structure and regulatory elements interacting at promoter regions, is a key step in plant responses to environmental cues. Nitrate (NO3-) is a nutrient signal that regulates the expression of hundreds of genes in Arabidopsis thaliana. Here, we integrate mRNA sequencing, genome-wide RNA polymerase II (RNPII), chromatin immunoprecipitation sequencing, and DNase sequencing datasets to establish the relationship between RNPII occupancy and chromatin accessibility in response to NO3- treatments in Arabidopsis roots. Genomic footprinting allowed us to identify in vivo regulatory elements controlling gene expression in response to NO3- treatments. NO3--modulated transcription factor (TF) footprints are important for a rapid increase in RNPII occupancy and transcript accumulation over time. We mapped key TF regulatory interactions and functionally validated the role of NAP, an NAC-domain containing TF, as a new regulatory factor in NO3- transport. Taken together, our study provides a comprehensive view of transcriptional networks in response to a nutrient signal in Arabidopsis roots.

Metabolic Footprinting of a Clear Cell Renal Cell Carcinoma in Vitro Model for Human Kidney Cancer Detection.

A protocol for harvesting and extracting extracellular metabolites from an in vitro model of human renal cell lines was developed to profile the exometabolome by means of a discovery-based metabolomics approach using ultraperformance liquid chromatography coupled to quadrupole-time-of-flight mass spectrometry. Metabolic footprints provided by conditioned media (CM) samples ( n = 66) of two clear cell Renal Cell Carcinoma (ccRCC) cell lines with different genetic backgrounds and a nontumor renal cell line, were compared with the human serum metabolic profile of a pilot cohort ( n = 10) comprised of stage IV ccRCC patients and healthy individuals. Using a cross-validated orthogonal projection to latent structures-discriminant analysis model, a panel of 21 discriminant features selected by iterative multivariate classification, allowed differentiating control from tumor cell lines with 100% specificity, sensitivity, and accuracy. Isoleucine/leucine, phenylalanine, N-lactoyl-leucine, and N-acetyl-phenylalanine, and cysteinegluthatione disulfide (CYSSG) were identified by chemical standards, and hydroxyprolyl-valine was identified with MS and MS/MS experiments. A subset of 9 discriminant features, including the identified metabolites except for CYSSG, produced a fingerprint of classification value that enabled discerning ccRCC patients from healthy individuals. To our knowledge, this is the first time that N-lactoyl-leucine is associated with ccRCC. Results from this study provide a proof of concept that CM can be used as a serum proxy to obtain disease-related metabolic signatures.

Estimating 20-year land-use change and derived CO2 emissions associated with crops, pasture and forestry in Brazil and each of its 27 states.

Land-use change (LUC) in Brazil has important implications on global climate change, ecosystem services and biodiversity, and agricultural expansion plays a critical role in this process. Concerns over these issues have led to the need for estimating the magnitude and impacts associated with that, which are increasingly reported in the environmental assessment of products. Currently, there is an extensive debate on which methods are more appropriate for estimating LUC and related emissions and regionalized estimates are lacking for Brazil, which is a world leader in agricultural production (e.g. food, fibres and bioenergy). We developed a method for estimating scenarios of past 20-year LUC and derived CO2 emission rates associated with 64 crops, pasture and forestry in Brazil as whole and in each of its 27 states, based on time-series statistics and in accordance with most used carbon-footprinting standards. The scenarios adopted provide a range between minimum and maximum rates of CO2 emissions from LUC according to different possibilities of land-use transitions, which can have large impacts in the results. Specificities of Brazil, like multiple cropping and highly heterogeneous carbon stocks, are also addressed. The highest CO2 emission rates are observed in the Amazon biome states and crops with the highest rates are those that have undergone expansion in this region. Some states and crops showing large agricultural areas have low emissions associated, especially in southern and eastern Brazil. Native carbon stocks and time of agricultural expansion are the most decisive factors to the patterns of emissions. Some implications on LUC estimation methods and standards and on agri-environmental policies are discussed.

Screening of Transcription Factors Involved in Fetal Hemoglobin Regulation Using Phylogenetic Footprinting.

Fetal hemoglobin (Hb F) is an important genetic modulator of the beta-hemoglobinopathies. The regulation of Hb F levels is influenced by transcription factors. We used phylogenetic footprinting to screen transcription factors that have binding sites in HBG1 and HBG2 genes' noncoding regions in order to know the genetic determinants of the Hb F expression. Our analysis showed 354 conserved motifs in the noncoding regions of HBG1 gene and 231 motifs in the HBG2 gene between the analyzed species. Of these motifs, 13 showed relation to Hb F regulation: cell division cycle-5 (CDC5), myelo-blastosis viral oncogene homolog (c-MYB), transcription factor CP2 (TFCP2), GATA binding protein 1 (GATA-1), GATA binding protein 2 (GATA-2), nuclear factor erythroid 2 (NF-E2), nuclear transcription factor Y (NF-Y), runt-related transcription factor 1 (RUNX-1), T-cell acute lymphocytic leukemia 1 (TAL-1), YY1 transcription factor (YY1), beta protein 1 (BP1), chicken ovalbumin upstream promoter-transcription factor II (COUP-TFII), and paired box 1 (PAX-1). The last three motifs were conserved only in the noncoding regions of the HBG1 gene. The understanding of genetic elements involved in the maintenance of high Hb F levels may provide new efficient therapeutic strategies in the beta-hemoglobinopathies treatment, promoting reduction in clinical complications of these genetic disorders.