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.

Nakazawaea atacamensis f.a., sp. nov. a novel nonconventional fermentative ascomycetous yeast species from the Atacama Desert.

In this study, we describe Nakazawaea atacamensis f. a., sp. nov., a novel species obtained from Neltuma chilensis plant samples in Chile's hyperarid Atacama Desert. In total, three strains of N. atacamensis were obtained from independent N. chilensis samples (synonym Prosopis chilensis, Algarrobo). Two strains were obtained from bark samples, while the third strain was obtained from bark-exuded gum from another tree. The novel species was defined using molecular characteristics and subsequently characterized with respect to morphological, physiological, and biochemical properties. A neighbor-joining analysis using the sequences of the D1/D2 domains of the large subunit ribosomal RNA gene revealed that N. atacamensis clustered with Nakazawaea pomicola. The sequence of N. atacamensis differed from closely related species by 1.3%-5.2% in the D1/D2 domains. A phylogenomic analysis based on single-nucleotide polymorphism's data confirms that the novel species belongs to the genus Nakazawaea, where N. atacamensis clustered with N. peltata. Phenotypic comparisons demonstrated that N. atacamensis exhibited distinct carbon assimilation patterns compared to its related species. Genome sequencing of the strain ATA-11A-BT revealed a genome size of approximately 12.4 Mbp, similar to other Nakazawaea species, with 5116 protein-coding genes annotated using InterProScan. In addition, N. atacamensis exhibited the capacity to ferment synthetic wine must, representing a potential new yeast for mono or co-culture wine fermentations. This comprehensive study expands our understanding of the genus Nakazawaea and highlights the ecological and industrial potential of N. atacamensis in fermentation processes. The holotype of N. atacamensis sp. nov. is CBS 18375T . The Mycobank number is MB 849680.

Characterization of the signalling pathways involved in the repression of root nitrate uptake by nitrate in Arabidopsis thaliana.

In Arabidopsis thaliana, root high-affinity nitrate (NO3-) uptake depends mainly on NRT2.1, 2.4, and 2.5, which are repressed by high NO3- supply at the transcript level. For NRT2.1, this regulation is due to the action of (i) feedback down-regulation by N metabolites and (ii) repression by NO3- itself mediated by the transceptor NRT1.1(NPF6.3). However, for NRT2.4 and NRT2.5, the signalling pathway(s) remain unknown as do the molecular elements involved. Here we show that unlike NRT2.1, NRT2.4 and NRT2.5 are not induced in an NO3- reductase mutant but are up-regulated following replacement of NO3- by ammonium (NH4+) as the N source. Moreover, increasing the NO3- concentration in a mixed nutrient solution with constant NH4+ concentration results in a gradual repression of NRT2.4 and NRT2.5, which is suppressed in an nrt1.1 mutant. This indicates that NRT2.4 and NRT2.5 are subjected to repression by NRT1.1-mediated NO3- sensing, and not to feedback repression by reduced N metabolites. We further show that key regulators of NRT2 transporters, such as HHO1, HRS1, PP2C, LBD39, BT1, and BT2, are also regulated by NRT1.1-mediated NO3- sensing, and that several of them are involved in NO3- repression of NRT2.1, NRT2.4, and NRT2.5. Finally, we provide evidence that it is the phosphorylated form of NRT1.1 at the T101 residue, which is most active in triggering the NRT1.1-mediated NO3- regulation of all these genes. Altogether, these data led us to propose a regulatory model for high-affinity NO3- uptake in Arabidopsis, highlighting several NO3- transduction cascades downstream of the phosphorylated form of the NRT1.1 transceptor.

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.

Plant ecological genomics at the limits of life in the Atacama Desert.

The Atacama Desert in Chile-hyperarid and with high-ultraviolet irradiance levels-is one of the harshest environments on Earth. Yet, dozens of species grow there, including Atacama-endemic plants. Herein, we establish the Talabre-Lejía transect (TLT) in the Atacama as an unparalleled natural laboratory to study plant adaptation to extreme environmental conditions. We characterized climate, soil, plant, and soil-microbe diversity at 22 sites (every 100 m of altitude) along the TLT over a 10-y period. We quantified drought, nutrient deficiencies, large diurnal temperature oscillations, and pH gradients that define three distinct vegetational belts along the altitudinal cline. We deep-sequenced transcriptomes of 32 dominant plant species spanning the major plant clades, and assessed soil microbes by metabarcoding sequencing. The top-expressed genes in the 32 Atacama species are enriched in stress responses, metabolism, and energy production. Moreover, their root-associated soils are enriched in growth-promoting bacteria, including nitrogen fixers. To identify genes associated with plant adaptation to harsh environments, we compared 32 Atacama species with the 32 closest sequenced species, comprising 70 taxa and 1,686,950 proteins. To perform phylogenomic reconstruction, we concatenated 15,972 ortholog groups into a supermatrix of 8,599,764 amino acids. Using two codon-based methods, we identified 265 candidate positively selected genes (PSGs) in the Atacama plants, 64% of which are located in Pfam domains, supporting their functional relevance. For 59/184 PSGs with an Arabidopsis ortholog, we uncovered functional evidence linking them to plant resilience. As some Atacama plants are closely related to staple crops, these candidate PSGs are a "genetic goldmine" to engineer crop resilience to face climate change.

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.

c-Abl regulates a synaptic plasticity-related transcriptional program involved in memory and learning.

Memory consolidation requires activation of a gene expression program that allows de novo protein synthesis. But the molecular mechanisms that favour or restrict that program are poorly understood. The kinase c-Abl can modulate gene expression through transcription factors and chromatin modifiers. Here, we show that c-Abl ablation in the brain improves learning acquisition and memory consolidation in mice. Its absence also affects gene expression profiles in the mouse hippocampus. We found that genes involved in synaptic plasticity and actin cytoskeleton dynamics, such as Arp2 and Thorase, are up-regulated at the mRNA and protein levels in trained c-Abl KO mice and by a chemical-LTP stimulus. Trained c-Abl KO mice also show that dendritic spines are larger than in wild-type mice and present at a higher density. These results indicate that c-Abl kinase is an important part of the mechanism that limits or restricts signalling of relevant gene programs involved in morphological and functional spine changes upon neuronal stimulation.

Genome-wide analysis in response to nitrogen and carbon identifies regulators for root AtNRT2 transporters.

In Arabidopsis (Arabidopsis thaliana), the High-Affinity Transport System (HATS) for root nitrate (NO3-) uptake depends mainly on four NRT2 NO3- transporters, namely NRT2.1, NRT2.2, NRT2.4, and NRT2.5. The HATS is the target of many regulations to coordinate nitrogen (N) acquisition with the N status of the plant and with carbon (C) assimilation through photosynthesis. At the molecular level, C and N signaling pathways control gene expression of the NRT2 transporters. Although several regulators of these transporters have been identified in response to either N or C signals, the response of NRT2 gene expression to the interaction of these signals has never been specifically investigated, and the underlying molecular mechanisms remain largely unknown. To address this question we used an original systems biology approach to model a regulatory gene network targeting NRT2.1, NRT2.2, NRT2.4, and NRT2.5 in response to N/C signals. Our systems analysis of the data identified three transcription factors, TGA3, MYC1, and bHLH093. Functional analysis of mutants combined with yeast one-hybrid experiments confirmed that all three transcription factors are regulators of NRT2.4 or NRT2.5 in response to N or C signals. These results reveal a role for TGA3, MYC1, and bHLH093 in controlling the expression of root NRT2 transporter genes.

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.

Loop assembly: a simple and open system for recursive fabrication of DNA circuits.

High-efficiency methods for DNA assembly have enabled the routine assembly of synthetic DNAs of increased size and complexity. However, these techniques require customization, elaborate vector sets or serial manipulations for the different stages of assembly. We have developed Loop assembly based on a recursive approach to DNA fabrication. The system makes use of two Type IIS restriction endonucleases and corresponding vector sets for efficient and parallel assembly of large DNA circuits. Standardized level 0 parts can be assembled into circuits containing 1, 4, 16 or more genes by looping between the two vector sets. The vectors also contain modular sites for hybrid assembly using sequence overlap methods. Loop assembly enables efficient and versatile DNA fabrication for plant transformation. We show the construction of plasmids up to 16 genes and 38 kb with high efficiency (> 80%). We have characterized Loop assembly on over 200 different DNA constructs and validated the fidelity of the method by high-throughput Illumina plasmid sequencing. Our method provides a simple generalized solution for DNA construction with standardized parts. The cloning system is provided under an OpenMTA license for unrestricted sharing and open access.

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.