Phylogenetically diverse wild plant species use common biochemical strategies to thrive in the Atacama Desert.

The best ideotypes are under mounting pressure due to increased aridity. Understanding the conserved molecular mechanisms that evolve in wild plants adapted to harsh environments is crucial in developing new strategies for agriculture. Yet our knowledge of such mechanisms in wild species is scant. We performed metabolic pathway reconstruction using transcriptome information from 32 Atacama and phylogenetically related species that do not live in Atacama (Sisters species). We analyzed reaction enrichment to understand the commonalities and differences of Atacama plants. To gain insights into the mechanisms that ensure survival, we compared expressed gene isoform numbers and gene expression patterns between the annotated biochemical reactions from 32 Atacama and Sister species. We found biochemical convergences characterized by reactions enriched in at least 50% of the Atacama species, pointing to potential advantages against drought and nitrogen starvation, for instance. These findings suggest that the adaptation in the Atacama Desert may result in part from shared genetic legacies governing the expression of key metabolic pathways to face harsh conditions. Enriched reactions corresponded to ubiquitous compounds common to extreme and agronomic species and were congruent with our previous metabolomic analyses. Convergent adaptive traits offer promising candidates for improving abiotic stress resilience in crop species.

The Botrytis cinerea Gene Expression Browser.

For comprehensive gene expression analyses of the phytopathogenic fungus Botrytis cinerea, which infects a number of plant taxa and is a cause of substantial agricultural losses worldwide, we developed BEB, a web-based B. cinerea gene Expression Browser. This computationally inexpensive web-based application and its associated database contain manually curated RNA-Seq data for B. cinerea. BEB enables expression analyses of genes of interest under different culture conditions by providing publication-ready heatmaps depicting transcript levels, without requiring advanced computational skills. BEB also provides details of each experiment and user-defined gene expression clustering and visualization options. If needed, tables of gene expression values can be downloaded for further exploration, including, for instance, the determination of differentially expressed genes. The BEB implementation is based on open-source computational technologies that can be deployed for other organisms. In this case, the new implementation will be limited only by the number of transcriptomic experiments that are incorporated into the platform. To demonstrate the usability and value of BEB, we analyzed gene expression patterns across different conditions, with a focus on secondary metabolite gene clusters, chromosome-wide gene expression, previously described virulence factors, and reference genes, providing the first comprehensive expression overview of these groups of genes in this relevant fungal phytopathogen. We expect this tool to be broadly useful in B. cinerea research, providing a basis for comparative transcriptomics and candidate gene identification for functional assays.

Spatiotemporal analysis identifies ABF2 and ABF3 as key hubs of endodermal response to nitrate.

Nitrate is a nutrient and a potent signal that impacts global gene expression in plants. However, the regulatory factors controlling temporal and cell type-specific nitrate responses remain largely unknown. We assayed nitrate-responsive transcriptome changes in five major root cell types of the Arabidopsis thaliana root as a function of time. We found that gene-expression response to nitrate is dynamic and highly localized and predicted cell type-specific transcription factor (TF)-target interactions. Among cell types, the endodermis stands out as having the largest and most connected nitrate-regulatory gene network. ABF2 and ABF3 are major hubs for transcriptional responses in the endodermis cell layer. We experimentally validated TF-target interactions for ABF2 and ABF3 by chromatin immunoprecipitation followed by sequencing and a cell-based system to detect TF regulation genome-wide. Validated targets of ABF2 and ABF3 account for more than 50% of the nitrate-responsive transcriptome in the endodermis. Moreover, ABF2 and ABF3 are involved in nitrate-induced lateral root growth. Our approach offers an unprecedented spatiotemporal resolution of the root response to nitrate and identifies important components of cell-specific gene regulatory networks.

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.

TGA class II transcription factors are essential to restrict oxidative stress in response to UV-B stress in Arabidopsis.

Plants possess a robust metabolic network for sensing and controlling reactive oxygen species (ROS) levels upon stress conditions. Evidence shown here supports a role for TGA class II transcription factors as critical regulators of genes controlling ROS levels in the tolerance response to UV-B stress in Arabidopsis. First, tga256 mutant plants showed reduced capacity to scavenge H2O2 and restrict oxidative damage in response to UV-B, and also to methylviologen-induced photooxidative stress. The TGA2 transgene (tga256/TGA2 plants) complemented these phenotypes. Second, RNAseq followed by clustering and Gene Ontology term analyses indicate that TGA2/5/6 positively control the UV-B-induced expression of a group of genes with oxidoreductase, glutathione transferase, and glucosyltransferase activities, such as members of the glutathione S-transferase Tau subfamily (GSTU), which encodes peroxide-scavenging enzymes. Accordingly, increased glutathione peroxidase activity triggered by UV-B was impaired in tga256 mutants. Third, the function of TGA2/5/6 as transcriptional activators of GSTU genes in the UV-B response was confirmed for GSTU7, GSTU8, and GSTU25, using quantitative reverse transcription-PCR and ChIP analyses. Fourth, expression of the GSTU7 transgene complemented the UV-B-susceptible phenotype of tga256 mutant plants. Together, this evidence indicates that TGA2/5/6 factors are key regulators of the antioxidant/detoxifying response to an abiotic stress such as UV-B light overexposure.

Transcriptomic profiles reveal differences in zinc metabolism, inflammation, and tight junction proteins in duodenum from cholesterol gallstone subjects.

Cholesterol Gallstone Disease (GSD) is a common multifactorial disorder characterized by crystallization and aggregation of biliary cholesterol in the gallbladder. The global prevalence of GSD is ~10-20% in the adult population but rises to 28% in Chile (17% among men and 30% among women). The small intestine may play a role in GSD pathogenesis, but the molecular mechanisms have not been clarified. Our aim was to identify the role of the small intestine in GSD pathogenesis. Duodenal biopsy samples were obtained from patients with GSD and healthy volunteers. GSD status was defined by abdominal ultrasonography. We performed a transcriptome study in a discovery cohort using Illumina HiSeq. 2500, and qPCR, immunohistochemistry and immunofluorescence were used to validate differentially expressed genes among additional case-control cohorts. 548 differentially expressed genes between GSD and control subjects were identified. Enriched biological processes related to cellular response to zinc, and immune and antimicrobial responses were observed in GSD patients. We validated lower transcript levels of metallothionein, NPC1L1 and tight junction genes and higher transcript levels of genes involved in immune and antimicrobial pathways in GSD patients. Interestingly, serum zinc and phytosterol to cholesterol precursor ratios were lower in GSD patients. A significant association was observed between serum zinc and phytosterol levels. Our results support a model where proximal small intestine plays a key role in GSD pathogenesis. Zinc supplementation, modulation of proximal microbiota and/or intestinal barrier may be novel targets for strategies to prevent GSD.

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.

Whole Genome Sequence, Variant Discovery and Annotation in Mapuche-Huilliche Native South Americans.

Whole human genome sequencing initiatives help us understand population history and the basis of genetic diseases. Current data mostly focuses on Old World populations, and the information of the genomic structure of Native Americans, especially those from the Southern Cone is scant. Here we present annotation and variant discovery from high-quality complete genome sequences of a cohort of 11 Mapuche-Huilliche individuals (HUI) from Southern Chile. We found approximately 3.1 × 106 single nucleotide variants (SNVs) per individual and identified 403,383 (6.9%) of novel SNVs events. Analyses of large-scale genomic events detected 680 copy number variants (CNVs) and 4,514 structural variants (SVs), including 398 and 1,910 novel events, respectively. Global ancestry composition of HUI genomes revealed that the cohort represents a sample from a marginally admixed population from the Southern Cone, whose main genetic component derives from Native American ancestors. Additionally, we found that HUI genomes contain variants in genes associated with 5 of the 6 leading causes of noncommunicable diseases in Chile, which may have an impact on the risk of prevalent diseases in Chilean and Amerindian populations. Our data represents a useful resource that can contribute to population-based studies and for the design of early diagnostics or prevention tools for Native and admixed Latin American populations.

MicroRNA Profiling and Bioinformatics Target Analysis in Dorsal Hippocampus of Chronically Stressed Rats: Relevance to Depression Pathophysiology.

Studies conducted in rodents subjected to chronic stress and some observations in humans after psychosocial stress, have allowed to establish a link between stress and the susceptibility to many complex diseases, including mood disorders. The studies in rodents have revealed that chronic exposure to stress negatively affects synaptic plasticity by triggering changes in the production of trophic factors, subunit levels of glutamate ionotropic receptors, neuron morphology, and neurogenesis in the adult hippocampus. These modifications may account for the impairment in learning and memory processes observed in chronically stressed animals. It is plausible then, that stress modifies the interplay between signal transduction cascades and gene expression regulation in the hippocampus, therefore leading to altered neuroplasticity and functioning of neural circuits. Considering that miRNAs play an important role in post-transcriptional-regulation of gene expression and participate in several hippocampus-dependent functions; we evaluated the consequences of chronic stress on the expression of miRNAs in dorsal (anterior) portion of the hippocampus, which participates in memory formation in rodents. Here, we show that male rats exposed to daily restraint stress (2.5 h/day) during 7 and 14 days display a differential profile of miRNA levels in dorsal hippocampus and remarkably, we found that some of these miRNAs belong to the miR-379-410 cluster. We confirmed a rise in miR-92a and miR-485 levels after 14 days of stress by qPCR, an effect that was not mimicked by chronic administration of corticosterone (14 days). Our in silico study identified the top-10 biological functions influenced by miR-92a, nine of which were shared with miR-485: Nervous system development and function, Tissue development, Behavior, Embryonic development, Organ development, Organismal development, Organismal survival, Tissue morphology, and Organ morphology. Furthermore, our in silico study provided a landscape of potential miRNA-92a and miR-485 targets, along with relevant canonical pathways related to axonal guidance signaling and cAMP signaling, which may influence the functioning of several neuroplastic substrates in dorsal hippocampus. Additionally, the combined effect of miR-92a and miR-485 on transcription factors, along with histone-modifying enzymes, may have a functional relevance by producing changes in gene regulatory networks that modify the neuroplastic capacity of the adult dorsal hippocampus under stress.

Step-by-Step Construction of Gene Co-expression Networks from High-Throughput Arabidopsis RNA Sequencing Data.

The rapid increase in the availability of transcriptomics data generated by RNA sequencing represents both a challenge and an opportunity for biologists without bioinformatics training. The challenge is handling, integrating, and interpreting these data sets. The opportunity is to use this information to generate testable hypothesis to understand molecular mechanisms controlling gene expression and biological processes (Fig. 1). A successful strategy to generate tractable hypotheses from transcriptomics data has been to build undirected network graphs based on patterns of gene co-expression. Many examples of new hypothesis derived from network analyses can be found in the literature, spanning different organisms including plants and specific fields such as root developmental biology.In order to make the process of constructing a gene co-expression network more accessible to biologists, here we provide step-by-step instructions using published RNA-seq experimental data obtained from a public database. Similar strategies have been used in previous studies to advance root developmental biology. This guide includes basic instructions for the operation of widely used open source platforms such as Bio-Linux, R, and Cytoscape. Even though the data we used in this example was obtained from Arabidopsis thaliana, the workflow developed in this guide can be easily adapted to work with RNA-seq data from any organism.

Use of transcriptomics and co-expression networks to analyze the interconnections between nitrogen assimilation and photorespiratory metabolism.

Nitrogen is one of the most important nutrients for plants and, in natural soils, its availability is often a major limiting factor for plant growth. Here we examine the effect of different forms of nitrogen nutrition and of photorespiration on gene expression in the model legume Lotus japonicus with the aim of identifying regulatory candidate genes co-ordinating primary nitrogen assimilation and photorespiration. The transcriptomic changes produced by the use of different nitrogen sources in leaves of L. japonicus plants combined with the transcriptomic changes produced in the same tissue by different photorespiratory conditions were examined. The results obtained provide novel information on the possible role of plastidic glutamine synthetase in the response to different nitrogen sources and in the C/N balance of L. japonicus plants. The use of gene co-expression networks establishes a clear relationship between photorespiration and primary nitrogen assimilation and identifies possible transcription factors connected to the genes of both routes.

Transcriptome analysis reveals regulatory networks underlying differential susceptibility to Botrytis cinerea in response to nitrogen availability in Solanum lycopersicum.

Nitrogen (N) is one of the main limiting nutrients for plant growth and crop yield. It is well documented that changes in nitrate availability, the main N source found in agricultural soils, influences a myriad of developmental programs and processes including the plant defense response. Indeed, many agronomical reports indicate that the plant N nutritional status influences their ability to respond effectively when challenged by different pathogens. However, the molecular mechanisms involved in N-modulation of plant susceptibility to pathogens are poorly characterized. In this work, we show that Solanum lycopersicum defense response to the necrotrophic fungus Botrytis cinerea is affected by plant N availability, with higher susceptibility in nitrate-limiting conditions. Global gene expression responses of tomato against B. cinerea under contrasting nitrate conditions reveals that plant primary metabolism is affected by the fungal infection regardless of N regimes. This result suggests that differential susceptibility to pathogen attack under contrasting N conditions is not only explained by a metabolic alteration. We used a systems biology approach to identify the transcriptional regulatory network implicated in plant response to the fungus infection under contrasting nitrate conditions. Interestingly, hub genes in this network are known key transcription factors involved in ethylene and jasmonic acid signaling. This result positions these hormones as key integrators of nitrate and defense against B. cinerea in tomato plants. Our results provide insights into potential crosstalk mechanisms between necrotrophic defense response and N status in plants.

Transcriptional networks in the nitrate response of Arabidopsis thaliana.

Nitrogen is an essential macronutrient for plants and its availability is a key determinant of plant growth and development and crop yield. Besides their nutritional role, N nutrients and metabolites are signals that activate signaling pathways that modulate many plant processes. Because the most abundant inorganic N source for plants in agronomic soils is nitrate, much of the work to understand plant N-signaling has focused on this nutrient. Over the last years, several studies defined a comprehensive catalog of nitrate-responsive genes, involved in nitrate transport, metabolism and a variety of other processes. Despite significant progress in recent years, primarily using Arabidopsis thaliana as a model system, the molecular mechanisms by which nitrate elicits changes in transcript abundance are still not fully understood. Here we highlight recent advancements in identifying key transcription factors and transcriptional mechanisms that orchestrate the gene expression response to changes in nitrate availability in A. thaliana.

Constructing simple biological networks for understanding complex high-throughput data in plants.

Technological advances in the last decade have enabled biologists to produce increasing amounts of information for the transcriptome, proteome, interactome, and other -omics data sets in many model organisms. A major challenge is integration and biological interpretation of these massive data sets in order to generate testable hypotheses about gene regulatory networks or molecular mechanisms that govern system behaviors. Constructing gene networks requires bioinformatics skills to adequately manage, integrate, analyze and productively use the data to generate biological insights. In this chapter, we provide detailed methods for users without prior knowledge of bioinformatics to construct gene networks and derive hypotheses that can be experimentally verified. Step-by-step instructions for acquiring, integrating, analyzing, and visualizing genome-wide data are provided for two widely used open source platforms, R and Cytoscape platforms. The examples provided are based on Arabidopsis data, but the protocols presented should be readily applicable to any organism for which similar data can be obtained.

Nitrogen control of developmental phase transitions in Arabidopsis thaliana.

Nitrogen (N) is an essential macronutrient and a key structural component of macromolecules in plants. N nutrients and metabolites can act as signals that impact on many aspects of plant biology. The plant life cycle involves a series of developmental phase transitions that must be tightly coordinated to external and internal cues in order to ensure plant survival and reproduction. N availability is one of the factors controlling phase changes. In this review, we integrate and summarize the known effects of N over different developmental stages in plants. Substantial advances have been made in our understanding of signalling and N-responsive gene regulatory networks. We focus on the molecular mechanisms underlying N regulation of developmental transitions and the role of putative new regulators that might link N availability to pathways controlling Arabidopsis growth and development from seed germination through the plant reproductive transition.

Systems analysis of transcriptome data provides new hypotheses about Arabidopsis root response to nitrate treatments.

Nitrogen (N) is an essential macronutrient for plant growth and development. Plants adapt to changes in N availability partly by changes in global gene expression. We integrated publicly available root microarray data under contrasting nitrate conditions to identify new genes and functions important for adaptive nitrate responses in Arabidopsis thaliana roots. Overall, more than 2000 genes exhibited changes in expression in response to nitrate treatments in Arabidopsis thaliana root organs. Global regulation of gene expression by nitrate depends largely on the experimental context. However, despite significant differences from experiment to experiment in the identity of regulated genes, there is a robust nitrate response of specific biological functions. Integrative gene network analysis uncovered relationships between nitrate-responsive genes and 11 highly co-expressed gene clusters (modules). Four of these gene network modules have robust nitrate responsive functions such as transport, signaling, and metabolism. Network analysis hypothesized G2-like transcription factors are key regulatory factors controlling transport and signaling functions. Our meta-analysis highlights the role of biological processes not studied before in the context of the nitrate response such as root hair development and provides testable hypothesis to advance our understanding of nitrate responses in plants.

Integrated RNA-seq and sRNA-seq analysis identifies novel nitrate-responsive genes in Arabidopsis thaliana roots.

BACKGROUND: Nitrate and other nitrogen metabolites can act as signals that regulate global gene expression in plants. Adaptive changes in plant morphology and physiology triggered by changes in nitrate availability are partly explained by these changes in gene expression. Despite several genome-wide efforts to identify nitrate-regulated genes, no comprehensive study of the Arabidopsis root transcriptome under contrasting nitrate conditions has been carried out. RESULTS: In this work, we employed the Illumina high throughput sequencing technology to perform an integrated analysis of the poly-A + enriched and the small RNA fractions of the Arabidopsis thaliana root transcriptome in response to nitrate treatments. Our sequencing strategy identified new nitrate-regulated genes including 40 genes not represented in the ATH1 Affymetrix GeneChip, a novel nitrate-responsive antisense transcript and a new nitrate responsive miRNA/TARGET module consisting of a novel microRNA, miR5640 and its target, AtPPC3. CONCLUSIONS: Sequencing of small RNAs and mRNAs uncovered new genes, and enabled us to develop new hypotheses for nitrate regulation and coordination of carbon and nitrogen metabolism.

Systems approaches map regulatory networks downstream of the auxin receptor AFB3 in the nitrate response of Arabidopsis thaliana roots.

Auxin is a key phytohormone regulating central processes in plants. Although the mechanism by which auxin triggers changes in gene expression is well understood, little is known about the specific role of the individual members of the TIR1/AFB auxin receptors, Aux/IAA repressors, and ARF transcription factors and/or molecular pathways acting downstream leading to plant responses to the environment. We previously reported a role for AFB3 in coordinating primary and lateral root growth to nitrate availability. In this work, we used an integrated genomics, bioinformatics, and molecular genetics approach to dissect regulatory networks acting downstream of AFB3 that are activated by nitrate in roots. We found that the NAC4 transcription factor is a key regulatory element controlling a nitrate-responsive network, and that nac4 mutants have altered lateral root growth but normal primary root growth in response to nitrate. This finding suggests that AFB3 is able to activate two independent pathways to control root system architecture. Our systems approach has unraveled key components of the AFB3 regulatory network leading to changes in lateral root growth in response to nitrate.