Protein Binding Microarray Research Articles

Rice protein-binding microarrays: a tool to detect cis-acting elements near promoter regions in rice

Main conclusionThe present study showed that a rice (Oryza sativa)-specific protein-binding microarray (RPBM) can be applied to analyze DNA-binding motifs with a TF where binding is evaluated in extended natural promoter regions. The analysis may facilitate identifying TFs and their downstream genes and constructing gene networks through cis-elements.Transcription factors (TFs) regulate gene expression at the transcriptional level by binding a specific DNA sequence. Thus, predicting the DNA-binding motifs of TFs is one of the most important areas in the functional analysis of TFs in the postgenomic era. Although many methods have been developed to address this challenge, many TFs still have unknown DNA-binding motifs. In this study, we designed RPBM with 40-bp probes and 20-bp of overlap, yielding 49 probes spanning the 1-kb upstream region before the translation start site of each gene in the entire genome. To confirm the efficiency of RPBM technology, we selected two previously studied TFs, OsWOX13 and OsSMF1, and an uncharacterized TF, OsWRKY34. We identified the ATTGATTG and CCACGTCA DNA-binding sequences of OsWOX13 and OsSMF1, respectively. In total, 635 and 932 putative feature genes were identified for OsWOX13 and OsSMF1, respectively. We discovered the CGTTGACTTT DNA-binding sequence and 195 putative feature genes of OsWRKY34. RPBM could be applicable in the analysis of DNA-binding motifs for TFs where binding is evaluated in the promoter and 5′ upstream CDS regions. The analysis may facilitate identifying TFs and their downstream genes and constructing gene networks through cis-elements.

Adaptive gene misregulation.

Because gene expression is important for evolutionary adaptation, its misregulation is an important cause of maladaptation. A misregulated gene can be incorrectly silent ("off") when a transcription factor (TF) that is required for its activation does not binds its regulatory region. Conversely, a misregulated gene can be incorrectly active ("on") when a TF not normally involved in its activation binds its regulatory region, a phenomenon also known as regulatory crosstalk. DNA mutations that destroy or create TF binding sites on DNA are an important source of misregulation and crosstalk. Although misregulation reduces fitness in an environment to which an organism is well-adapted, it may become adaptive in a new environment. Here, I derive simple yet general mathematical expressions that delimit the conditions under which misregulation can be adaptive. These expressions depend on the strength of selection against misregulation, on the fraction of DNA sequence space filled with TF binding sites, and on the fraction of genes that must be expressed for optimal adaptation. I then use empirical data from RNA sequencing, protein-binding microarrays, and genome evolution, together with population genetic simulations to ask when these conditions are likely to be met. I show that they can be met under realistic circumstances, but these circumstances may vary among organisms and environments. My analysis provides a framework in which improved theory and data collection can help us demonstrate the role of misregulation in adaptation. It also shows that misregulation, like DNA mutation, is one of life's many imperfections that can help propel Darwinian evolution.

Zinc Finger Protein SALL4 Functions through an AT-Rich Motif to Regulate Gene Expression

SUMMARYThe zinc finger transcription factor SALL4 is highly expressed in embryonic stem cells, downregulated in most adult tissues, but reactivated in many aggressive cancers. This unique expression pattern makes SALL4 an attractive therapeutic target. However, whether SALL4 binds DNA directly to regulate gene expression is unclear, and many of its targets in cancer cells remain elusive. Here, through an unbiased screen of protein binding microarray (PBM) and cleavage under targets and release using nuclease (CUT&RUN) experiments, we identify and validate the DNA binding domain of SALL4 and its consensus binding sequence. Combined with RNA sequencing (RNA-seq) analyses after SALL4 knockdown, we discover hundreds of new SALL4 target genes that it directly regulates in aggressive liver cancer cells, including genes encoding a family of histone 3 lysine 9-specific demethylases (KDMs). Taken together, these results elucidate the mechanism of SALL4 DNA binding and reveal pathways and molecules to target in SALL4-dependent tumors.

DAP-Seq Identification of Transcription Factor-Binding Sites in Potato.

Plant growth and adaptation to environmental fluctuations involve a tight control of cellular processes which, to a great extent, are mediated by changes at the transcriptional level. This regulation is exerted by transcription factors (TFs), a group of regulatory proteins that control gene expression by directly binding to the gene promoter regions via their cognate TF-binding sites (TFBS). The nature of TFBS defines the pattern of expression of the various plant loci, the precise combinatorial assembly of these elements being key in conferring plant's adaptation ability and in domestication. As such, TFs are main potential targets for biotechnological interventions, prompting in the last decade notable protein-DNA interaction efforts toward definition of their TFBS. Distinct methods based on in vivo or in vitro approaches defined the TFBS for many TFs, mainly in Arabidopsis, but comprehensive information on the transcriptional networks for many regulators is still lacking, especially in crops. In this chapter, detailed protocols for DAP-seq studies to unbiased identification of TFBS in potato are provided. This methodology relies on the affinity purification of genomic DNA-protein complexes in vitro, and high-throughput sequencing of the eluted DNA fragments. DAP-seq outperforms other in vitro DNA-motif definition strategies, such as protein-binding microarrays and SELEX-seq, since the protein of interest is directly bound to the genomic DNA extracted from plants yielding all the potential sites bound by the TF in the genome. Actually, data generated from DAP-seq experiments are highly similar to those out of ChIP-seq methods, but are generated much faster. We also provide a standard procedure to the analysis of the DAP-seq data, addressed to non-experienced users, that involves two consecutive steps: (1) processing of raw data (trimming, filtering, and read alignment) and (2) peak calling and identification of enriched motifs. This method allows identification of the binding profiles of dozens of TFs in crops, in a timely manner.

High-Throughput Analysis of the Cell and DNA Site-Specific Binding of Native NF-κB Dimers Using Nuclear Extract Protein-Binding Microarrays (NextPBMs).

Nuclear factor-kappa B (NF-κB) transcription factors coordinate gene expression in response to a broad array of cellular signals. In vertebrates, there are five NF-κB proteins (c-Rel, RelA/p65, RelB, p50, and p52) that can form various dimeric combinations exhibiting both common and dimer-specific DNA-binding specificity. In this chapter, we describe the use of the nuclear extract protein-binding microarray (nextPBM), a high-throughput method to characterize the DNA binding of transcription factors present in cell nuclear extracts. NextPBMs allow for sensitive analysis of the DNA binding of NF-κB dimers and their interactions with cell-specific cofactors.

Transcription Factor Binding Affinities and DNA Shape Readout.

SummaryAn essential event in gene regulation is the binding of a transcription factor (TF) to its target DNA. Models considering the interactions between the TF and the DNA geometry proved to be successful approaches to describe this binding event, while conserving data interpretability. However, a direct characterization of the DNA shape contribution to binding is still missing due to the lack of accurate and large-scale binding affinity data. Here, we use a binding assay we recently established to measure with high sensitivity the binding specificities of 13 Drosophila TFs, including dinucleotide dependencies to capture non-independent amino acid-base interactions. Correlating the binding affinities with all DNA shape features, we find that shape readout is widely used by these factors. A shape readout/TF-DNA complex structure analysis validates our approach while providing biological insights such as positively charged or highly polar amino acids often contact nucleotides that exhibit strong shape readout.

BZIP Dimers CREB1, ATF2, Zta, ATF3|cJun, and cFos|cJun Prefer to Bind to Some Double-Stranded DNA Sequences Containing 5-Formylcytosine and 5-Carboxylcytosine.

In mammalian cells, 5-methylcytosine (5mC) occurs in genomic double-stranded DNA (dsDNA) and is enzymatically oxidized to 5-hydroxymethylcytosine (5hmC), then to 5-formylcytosine (5fC), and finally to 5-carboxylcytosine (5caC). These cytosine modifications are enriched in regulatory regions of the genome. The effect of these oxidative products on five bZIP dimers (CREB1, ATF2, Zta, ATF3|cJun, and cFos|cJun) binding to five types of dsDNA was measured using protein binding microarrays. The five dsDNAs contain either cytosine in both DNA strands or cytosine in one strand and either 5mC, 5hmC, 5fC, or 5caC in the second strand. Some sequences containing the CEBP half-site GCAA are bound more strongly by all five bZIP domains when dsDNA contains 5mC, 5hmC, or 5fC. dsDNA containing 5caC in some TRE (AP-1)-like sequences, e.g., TGACTAA, is better bound by Zta, ATF3|cJun, and cFos|cJun.

Zinc Finger Protein SALL4 Functions Through an AT-Rich Motif to Regulate Gene Expression

The zinc finger transcription factor SALL4 is highly expressed in embryonic stem cells, down-regulated in most adult tissues, but reactivated in many aggressive cancers. This unique expression pattern makes SALL4 an attractive target for designing therapeutic strategies. However, whether SALL4 binds DNA directly to regulate gene expression is unclear and many of its targets in cancer cells remain elusive. Here, through an unbiased screen of protein binding microarray (PBM) and Cleavage Under Targets and Release Using Nuclease (CUT&RUN) experiments, we identified and validated the DNA binding domain of SALL4 and its consensus binding sequence. Combined with RNA-seq analyses after SALL4 knockdown, we discovered hundreds of new SALL4 target genes that it directly regulates in aggressive liver cancer cells, including genes encoding a family of Histone 3 Lysine 9-specific Demethylases (KDMs). Taken together, these results elucidated the mechanism of SALL4 DNA binding and revealed novel pathways and molecules to target in SALL4-dependent tumors.

79: Genetic variants at AGTR2 locus increase risk of preterm birth

Preterm birth (PTB) complicates approximately 10% of pregnancies in the U.S. While genetic factors contribute to PTB, there has been limited replication of specific genetic risk loci associated with PTB. Recently, we conducted the first large-scale, maternal genome-wide association study (GWAS), utilizing 43,568 women from 23andMe, that identified a locus of variants near the angiotensin II receptor type 2 receptor (AGTR2) gene associated with both PTB and gestational length. We hypothesized that a causative variant in this locus influences expression of AGTR2. We assessed genetic variants previously identified in our GWAS located near AGTR2 (combined analysis p-value of 7.2x10-16 in association with gestational length) to identify genotype-dependent regulation of AGTR2 expression. We computationally prioritized disease variants using available functional genomics data. Implicated variants were then evaluated for possible effects on transcription factor (TF) binding using protein-binding microarray data. Selected variants predicted to alter TF binding were analyzed using electrophoretic mobility shift assays (EMSAs). Luciferase reporter assays were performed to identify rs7889204 genotype-dependent regulation of the AGTR2 promoter. Our computational analyses identified rs7889204 as the top functional candidate variant. EMSA analysis confirmed predicted allele-dependent binding for the transcription factors CEBPB and HOXA10. Luciferase reporter assays revealed allele-dependent reporter activity for rs7889204, with decreased expression for the non-reference “C” allele. We identified candidate mechanisms in which rs7889204 at the AGTR2 locus differentially binds CEBPB and HOXA10. Reporter assays demonstrate decreased enhancer activity for the non-reference allele at rs7889204. These results are consistent with a previously established eQTL for AGTR2 in the uterus, suggesting decreased binding of CEBPB and HOXA10 to the preterm birth rs7889204 risk allele results in decreased AGTR2 expression and possible downstream decidual implantation abnormalities leading to PTB.View Large Image Figure ViewerDownload Hi-res image Download (PPT)

Context and number of noncanonical repeat variable diresidues impede the design of TALE proteins with improved DNA targeting.

Transcription activator-like effector (TALE) proteins have been used extensively for targeted binding of fusion proteins to loci of interest in (epi)genome engineering. Such approaches typically utilize four canonical TALE repeat variable diresidue (RVD) types, corresponding to the identities of two key amino acids, to target each nucleotide. Alternate RVDs with improved specificity are desired. Here, we focused on seven noncanonical RVDs that have been suggested to have improved specificity for their target nucleotides. We used custom protein binding microarrays to characterize the DNA-binding activity of 65 TALEs containing these alternate or corresponding canonical RVDs at multiple positions to ~5,000 unique DNA sequences per protein. We found that none of the noncanonical thymine-targeting RVDs displayed stronger preference for thymine than did the canonical RVD. Of the noncanonical RVDs with putatively improved specificity for guanine, only EN and NH showed greater discrimination of guanine over adenine. This improved specificity, however, comes at a cost: more substitutions of a noncanonical RVD for a canonical RVD generally decreased the protein's DNA-binding activity. Our results highlight the need to investigate RVD-nucleotide specificities in multiple protein contexts and suggest that a balance between canonical and noncanonical RVDs is needed to build TALEs with improved specificity.

Predicting in-vitro Transcription Factor Binding Sites Using DNA Sequence + Shape.

Discovery of transcription factor binding sites (TFBSs) is essential for understanding the underlying binding mechanisms and cellular functions. Recently, Convolutional neural network (CNN) has succeeded in predicting TFBSs from the primary DNA sequences. In addition to DNA sequences, several evidences suggest that protein-DNA binding is partly mediated by properties of DNA shape. Although many methods have been proposed to jointly account for DNA sequences and shape properties in predicting TFBSs, they ignore the power of the combination of deep learning and DNA sequence + shape. Therefore we develop a deep-learning-based sequence + shape framework (DLBSS) in this paper, which appropriately integrates DNA sequences and shape properties, to better understand protein-DNA binding preference. This method uses a shared CNN to find their common patterns from DNA sequences and their corresponding shape features, which are then concatenated to compute a predicted value. Using 66 in-vitro datasets derived from universal protein binding microarrays (uPBMs), we show that our proposed method DLBSS significantly improves the performance of predicting TFBSs. In addition, we explain the reason why we should use the shared CNN, and explore the performance of DLBSS when using a deeper CNN, through a series of experiments.

A map of direct TF-DNA interactions in the human genome.

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is the most popular assay to identify genomic regions, called ChIP-seq peaks, that are bound in vivo by transcription factors (TFs). These regions are derived from direct TF–DNA interactions, indirect binding of the TF to the DNA (through a co-binding partner), nonspecific binding to the DNA, and noise/bias/artifacts. Delineating the bona fide direct TF–DNA interactions within the ChIP-seq peaks remains challenging. We developed a dedicated software, ChIP-eat, that combines computational TF binding models and ChIP-seq peaks to automatically predict direct TF–DNA interactions. Our work culminated with predicted interactions covering >2% of the human genome, obtained by uniformly processing 1983 ChIP-seq peak data sets from the ReMap database for 232 unique TFs. The predictions were a posteriori assessed using protein binding microarray and ChIP-exo data, and were predominantly found in high quality ChIP-seq peaks. The set of predicted direct TF–DNA interactions suggested that high-occupancy target regions are likely not derived from direct binding of the TFs to the DNA. Our predictions derived co-binding TFs supported by protein-protein interaction data and defined cis-regulatory modules enriched for disease- and trait-associated SNPs. We provide this collection of direct TF–DNA interactions and cis-regulatory modules through the UniBind web-interface (http://unibind.uio.no).

Comprehensive study of nuclear receptor DNA binding provides a revised framework for understanding receptor specificity

The type II nuclear receptors (NRs) function as heterodimeric transcription factors with the retinoid X receptor (RXR) to regulate diverse biological processes in response to endogenous ligands and therapeutic drugs. DNA-binding specificity has been proposed as a primary mechanism for NR gene regulatory specificity. Here we use protein-binding microarrays (PBMs) to comprehensively analyze the DNA binding of 12 NR:RXRα dimers. We find more promiscuous NR-DNA binding than has been reported, challenging the view that NR binding specificity is defined by half-site spacing. We show that NRs bind DNA using two distinct modes, explaining widespread NR binding to half-sites in vivo. Finally, we show that the current models of NR specificity better reflect binding-site activity rather than binding-site affinity. Our rich dataset and revised NR binding models provide a framework for understanding NR regulatory specificity and will facilitate more accurate analyses of genomic datasets.

The bZIP mutant CEBPB (V285A) has sequence specific DNA binding propensities similar to CREB1

The bZIP homodimers CEBPB and CREB1 bind DNA containing methylated cytosines differently. CREB1 binds stronger to the C/EBP half-site GCAA when the cytosine is methylated. For CEBPB, methylation of the same cytosine does not affect DNA binding. The X-ray structure of CREB1 binding the half site GTCA identifies an alanine in the DNA binding region interacting with the methyl group of T, structurally analogous to the methyl group of methylated C. This alanine is replaced with a valine in CEBPB. To explore the contribution of this amino acid to binding with methylated cytosine of the GCAA half-site, we made the reciprocal mutants CEBPB(V285A) and CREB1(A297V) and used protein binding microarrays (PBM) to examine binding to four types of double-stranded DNA (dsDNA): 1) DNA with cytosine in both strands (DNA(C|C)), 2) DNA with 5-methylcytosine (M) in one strand and cytosine in the second strand (DNA(M|C)), 3) DNA with 5-hydroxymethylcytosine (H) in one strand and cytosine in the second strand (DNA(H|C)), and 4) DNA with both cytosines in all CG dinucleotides containing 5-methylcytosine (DNA(5mCG)). When binding to DNA(C|C), CEBPB (V285A) preferentially binds the CRE consensus motif (TGACGTCA), similar to CREB1. The reciprocal mutant, CREB1(A297V) binds DNA with some similarity to CEBPB, with strongest binding to the methylated PAR site 8-mer TTACGTAA. These data demonstrate that V285 residue inhibits CEBPB binding to methylated cytosine of the GCAA half-site.

NextPBM: a platform to study cell-specific transcription factor binding and cooperativity.

High-throughput (HT) in vitro methods for measuring protein-DNA binding have become invaluable for characterizing transcription factor (TF) complexes and modeling gene regulation. However, current methods do not utilize endogenous proteins and, therefore, do not quantify the impact of cell-specific post-translational modifications (PTMs) and cooperative cofactors. We introduce the HT nextPBM (nuclear extract protein-binding microarray) approach to study DNA binding of native cellular TFs that accounts for PTMs and cell-specific cofactors. We integrate immune-depletion and phosphatase treatment steps into our nextPBM pipeline to characterize the impact of cofactors and phosphorylation on TF binding. We analyze binding of PU.1/SPI1 and IRF8 from human monocytes, delineate DNA-sequence determinants for their cooperativity, and show how PU.1 affinity correlates with enhancer status and the presence of cooperative and collaborative cofactors. We describe how nextPBMs, and our accompanying computational framework, can be used to discover cell-specific cofactors, screen for synthetic cooperative DNA elements, and characterize TF cooperativity.

A map of direct TF-DNA interactions in the human genome.

Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is the most popular assay to identify genomic regions, called ChIP-seq peaks, that are bound in vivo by transcription factors (TFs). These regions are derived from direct TF–DNA interactions, indirect binding of the TF to the DNA (through a co-binding partner), nonspecific binding to the DNA, and noise/bias/artifacts. Delineating the bona fide direct TF–DNA interactions within the ChIP-seq peaks remains challenging. We developed a dedicated software, ChIP-eat, that combines computational TF binding models and ChIP-seq peaks to automatically predict direct TF–DNA interactions. Our work culminated with predicted interactions covering >2% of the human genome, obtained by uniformly processing 1983 ChIP-seq peak data sets from the ReMap database for 232 unique TFs. The predictions were a posteriori assessed using protein binding microarray and ChIP-exo data, and were predominantly found in high quality ChIP-seq peaks. The set of predicted direct TF–DNA interactions suggested that high-occupancy target regions are likely not derived from direct binding of the TFs to the DNA. Our predictions derived co-binding TFs supported by protein-protein interaction data and defined cis-regulatory modules enriched for disease- and trait-associated SNPs. We provide this collection of direct TF–DNA interactions and cis-regulatory modules through the UniBind web-interface (http://unibind.uio.no).

OsTGA2 confers disease resistance to rice against leaf blight by regulating expression levels of disease related genes via interaction with NH1.

How plants defend themselves from microbial infection is one of the most critical issues for sustainable crop production. Some TGA transcription factors belonging to bZIP superfamily can regulate disease resistance through NPR1-mediated immunity mechanisms in Arabidopsis. Here, we examined biological roles of OsTGA2 (grouped into the same subclade as Arabidopsis TGAs) in bacterial leaf blight resistance. Transcriptional level of OsTGA2 was accumulated after treatment with salicylic acid, methyl jasmonate, and Xathomonas oryzae pv. Oryzae (Xoo), a bacterium causing serious blight of rice. OsTGA2 formed homo- and hetero-dimer with OsTGA3 and OsTGA5 and interacted with rice NPR1 homologs 1 (NH1) in rice. Results of quadruple 9-mer protein-binding microarray analysis indicated that OsTGA2 could bind to TGACGT DNA sequence. Overexpression of OsTGA2 increased resistance of rice to bacterial leaf blight, although overexpression of OsTGA3 resulted in disease symptoms similar to wild type plant upon Xoo infection. Overexpression of OsTGA2 enhanced the expression of defense related genes containing TGA binding cis-element in the promoter such as AP2/EREBP 129, ERD1, and HOP1. These results suggest that OsTGA2 can directly regulate the expression of defense related genes and increase the resistance of rice against bacterial leaf blight disease.

How bacterial xenogeneic silencer rok distinguishes foreign from self DNA in its resident genome.

Bacterial xenogeneic silencers play important roles in bacterial evolution by recognizing and inhibiting expression from foreign genes acquired through horizontal gene transfer, thereby buffering against potential fitness consequences of their misregulated expression. Here, the detailed DNA binding properties of Rok, a xenogeneic silencer in Bacillus subtilis, was studied using protein binding microarray, and the solution structure of its C-terminal DNA binding domain was determined in complex with DNA. The C-terminal domain of Rok adopts a typical winged helix fold, with a novel DNA recognition mechanism different from other winged helix proteins or xenogeneic silencers. Rok binds the DNA minor groove by forming hydrogen bonds to bases through N154, T156 at the N-terminal of α3 helix and R174 of wing W1, assisted by four lysine residues interacting electrostatically with DNA backbone phosphate groups. These structural features endow Rok with preference towards DNA sequences harboring AACTA, TACTA, and flexible multiple TpA steps, while rigid A-tracts are disfavored. Correspondingly, the Bacillus genomes containing Rok are rich in A-tracts and show a dramatic underrepresentation of AACTA and TACTA, which are significantly enriched in Rok binding regions. These observations suggest that the xenogeneic silencing protein and its resident genome may have evolved cooperatively.

A WUSCHEL Homeobox Transcription Factor, OsWOX13, Enhances Drought Tolerance and Triggers Early Flowering in Rice.

Plants have evolved strategies to cope with drought stress by maximizing physiological capacity and adjusting developmental processes such as flowering time. The WOX13 orthologous group is the most conserved among the clade of WOX homeodomain-containing proteins and is found to function in both drought stress and flower development. In this study, we isolated and characterized OsWOX13 from rice. OsWOX13 was regulated spatially in vegetative organs but temporally in flowers and seeds. Overexpression of OsWOX13 (OsWOX13-ov) in rice under the rab21 promoter resulted in drought resistance and early flowering by 7–10 days. Screening of gene expression profiles in mature leaf and panicles of OsWOX13-ov showed a broad spectrum of effects on biological processes, such as abiotic and biotic stresses, exerting a cross-talk between responses. Protein binding microarray and electrophoretic mobility shift assay analyses supported ATTGATTG as the putative cis-element binding of OsWOX13. OsDREB1A and OsDREB1F, drought stress response transcription factors, contain ATTGATTG motif(s) in their promoters and are preferentially expressed in OsWOX13-ov. In addition, Heading date 3a and OsMADS14, regulators in the flowering pathway and development, were enhanced in OsWOX13-ov. These results suggest that OsWOX13 mediates the stress response and early flowering and, thus, may be a regulator of genes involved in drought escape.

Coordination of Chloroplast Development through the Action of the GNC and GLK Transcription Factor Families.

Fundamental questions regarding how chloroplasts develop from proplastids remain poorly understood despite their central importance to plant life. Two families of nuclear transcription factors, the GATA NITRATE-INDUCIBLE CARBON-METABOLISM-INVOLVED (GNC) and GOLDEN TWO-LIKE (GLK) families, have been implicated in directly and positively regulating chloroplast development. Here, we determined the degree of functional overlap between the two transcription factor families in Arabidopsis (Arabidopsis thaliana), characterizing their ability to regulate chloroplast biogenesis both alone and in concert. We determined the DNA-binding motifs for GNC and GLK2 using protein-binding microarrays; the enrichment of these motifs in transcriptome datasets indicates that GNC and GLK2 are repressors and activators of gene expression, respectively. ChIP-seq analysis of GNC identified PHYTOCHROME INTERACTING FACTOR and brassinosteroid activity genes as targets whose repression by GNC facilitates chloroplast biogenesis. In addition, GNC targets and represses genes involved in ERECTA signaling and thereby facilitates stomatal development. Our results define key regulatory features of the GNC and GLK transcription factor families that contribute to the control of chloroplast biogenesis and photosynthetic activity, including areas of independence and cross talk.