Eukaryotic transcription factors

Evolutionary analysis of the OSCA gene family in sunflower (Helianthus annuus L) and expression analysis under NaCl stress

Hyperosmolality-gated calcium-permeable channels (OSCA) are Ca^2 + nonselective cation channels that contain the calcium-dependent DUF221 domain, which plays an important role in plant response to stress and growth. However, the OSCA gene has not been fully identified and analyzed in sunflowers. In this study, we comprehensively analyzed the number, structure, collinearity, and phylogeny of the OSCA gene family in the sunflower, six Compositae species (Arctium lappa, Chrysanthemum morifolium, Cichorium endivia, Cichorium intybus, Lactuca sativa var. Angustata, and Carthamus tinctorius), and six other plants (soybean, Arabidopsis thaliana, rice, grape, and maize). The expression of the sunflower OSCA gene in nine different tissues, six different hormones, and NaCl stress conditions were analyzed based on transcriptome data and qRT-PCR. A total of 15 OSCA proteins, distributed on 10 chromosomes, were identified in the sunflower, and all of them were located in the endoplasmic reticulum. Using the phylogenetic tree, collinearity, gene structure, and motif analysis of the six Compositae species and six other plants, we found that the sunflower OSCA protein had only three subfamilies and lacked the Group 4 subfamily, which is conserved in the evolution of Compositae and subject to purification selection. The OSCA gene structure and motif analysis of the sunflower and six Compositae showed that there was a positive correlation between the number of motifs of most genes and the length of the gene, different subfamilies had different motifs, and the Group 4 subfamily had the smallest number of genes and the simplest gene structure. RNA-seq and qRT-PCR analysis showed that the expression levels of most OSCA genes in the sunflower changed to varying degrees under salt stress, and HaOSCA2.6 and HaOSCA3.1 were the most important in the sunflower's response to salt stress. The coexpression network of the sunflower genes under salt stress was constructed based on weighted gene co-expression network analysis (WGCNA). In conclusion, our findings suggest that the OSCA gene family is conserved during the sunflower's evolution and plays an important role in salt tolerance. These results will deepen our understanding of the evolutionary relationship of the sunflower OSCA gene family and provide a basis for their functional studies under salt stress.

Transcriptome Analysis and Functional Validation Identify a Putative bZIP Transcription Factor, Fpkapc, that Regulates Development, Stress Responses, and Virulence in Fusarium pseudograminearum

Fusarium pseudograminearum is a soilborne, hemibiotrophic phytopathogenic fungus that causes Fusarium crown rot and Fusarium head blight in wheat. The basic leucine zipper proteins (bZIPs) are evolutionarily conserved transcription factors that play crucial roles in a range of growth and developmental processes and the responses to biotic and abiotic stresses. However, the roles of bZIP transcription factors remains unknown in F. pseudograminearum. In this study, a bZIP transcription factor Fpkapc was identified to localize to the nucleus in F. pseudograminearum. A mutant strain (Δfpkapc) was constructed to determine the role of Fpkapc in growth and pathogenicity of F. pseudograminearum. Transcriptomic analyses revealed that many genes involved in basic metabolism and oxidation-reduction processes were downregulated, whereas many genes involved in metal iron binding were upregulated in the Δfpkapc strain, compared with the wild type (WT). Correspondingly, the mutant had severe growth defects and displayed abnormal hyphal tips. Conidiation in the Fpkapc mutant was reduced, with more conidia in smaller size and fewer septa than in the WT. Also, relative to WT, the Δfpkapc strain showed greater tolerance to ion stress, but decreased tolerance to H₂O₂. The mutant caused smaller disease lesions on wheat and barley plants, but significantly increased TRI gene expression, compared with the WT. In summary, Fpkapc plays multiple roles in governing growth, development, stress responses, and virulence in F. pseudograminearum.

Genome-wide identification and expression analysis of the WRKY transcription factor family in flax (Linum usitatissimum L.)

BACKGROUND

Members of the WRKY protein family, one of the largest transcription factor families in plants, are involved in plant growth and development, signal transduction, senescence, and stress resistance. However, little information is available about WRKY transcription factors in flax (Linum usitatissimum L.).

RESULTS

In this study, comprehensive genome-wide characterization of the flax WRKY gene family was conducted that led to prediction of 102 LuWRKY genes. Based on bioinformatics-based predictions of structural and phylogenetic features of encoded LuWRKY proteins, 95 LuWRKYs were classified into three main groups (Group I, II, and III); Group II LuWRKYs were further assigned to five subgroups (IIa-e), while seven unique LuWRKYs (LuWRKYs 96-102) could not be assigned to any group. Most LuWRKY proteins within a given subgroup shared similar motif compositions, while a high degree of motif composition variability was apparent between subgroups. Using RNA-seq data, expression patterns of the 102 predicted LuWRKY genes were also investigated. Expression profiling data demonstrated that most genes associated with cellulose, hemicellulose, or lignin content were predominantly expressed in stems, roots, and less in leaves. However, most genes associated with stress responses were predominantly expressed in leaves and exhibited distinctly higher expression levels in developmental stages 1 and 8 than during other stages.

CONCLUSIONS

Ultimately, the present study provides a comprehensive analysis of predicted flax WRKY family genes to guide future investigations to reveal functions of LuWRKY proteins during plant growth, development, and stress responses.

Time-resolved mRNA and miRNA expression profiling reveals crucial coregulation of molecular pathways involved in epithelial-pneumococcal interactions

Streptococcus pneumoniae is a major causative agent of pneumonia worldwide and its complex interaction with the lung epithelium has not been thoroughly characterized. In this study, we exploited both RNA‐sequencing and microRNA (miRNA)‐sequencing approaches to monitor the transcriptional changes in human lung alveolar epithelial cells infected by S. pneumoniae in a time‐resolved manner. A total of 1330 differentially expressed (DE) genes and 45 DE miRNAs were identified in all comparisons during the infection process. Clustering analysis showed that all DE genes were grouped into six clusters, several of which were primarily involved in inflammatory or immune responses. In addition, target gene enrichment analyses identified 11 transcription factors that were predicted to link at least one of four clusters, revealing transcriptional coregulation of multiple processes or pathways by common transcription factors. Notably, pharmacological treatment suggested that phosphorylation of p65 is important for optimal transcriptional regulation of target genes in epithelial cells exposed to pathogens. Furthermore, network‐based clustering analysis separated the DE genes negatively regulated by DE miRNAs into two functional modules (M1 and M2), with an enrichment in immune responses and apoptotic signaling pathways for M1. Integrated network analyses of potential regulatory interactions in M1 revealed that multiple DE genes related to immunity and apoptosis were regulated by multiple miRNAs, indicating the coordinated regulation of multiple genes by multiple miRNAs. In conclusion, time‐series expression profiling of messenger RNA and miRNA provides a wealth of information for global transcriptional changes, and offers comprehensive insight into the molecular mechanisms underlying host–pathogen interactions.

The bZIP transcription factor FpAda1 is essential for fungal growth and conidiation in Fusarium pseudograminearum

Fusarium pseudograminearum is an important pathogen of Fusarium crown rot and Fusarium head blight, which is able to infect wheat and barley worldwide, causing great economic losses. Transcription factors (TFs) of the basic leucine zipper (bZIP) protein family control important processes in all eukaryotes. In this study, we identified a gene, designated FpAda1, encoding a bZIP TF in F. pseudograminearum. The homolog of FpAda1 is also known to affect hyphal growth in Neurospora crassa. Deletion of FpAda1 in F. pseudograminearum resulted in defects in hyphal growth, mycelial branching and conidia formation. Pathogenicity assays showed that virulence of the Δfpada1 mutant was dramatically decreased on wheat coleoptiles and barley leaves. However, wheat coleoptile inoculation assay showed that Δfpada1 could penetrate and proliferate in wheat cells. Moreover, the FpAda1 was required for abnormal nuclear morphology in conidia and transcription of FpCdc2 and FpCdc42. Taken together, these results indicate that FpAda1 is an important transcription factor involved in growth and development in F. pseudograminearum.

Genome-wide identification, phylogeny, evolutionary expansion and expression analyses of bZIP transcription factor family in tartaty buckwheat

Background

In reported plants, the bZIP family is one of the largest transcription factor families. bZIP genes play roles in the light signal, seed maturation, flower development, cell elongation, seed accumulation protein, abiotic and biological stress and other biological processes. While, no detailed identification and genome-wide analysis of bZIP family genes in Fagopyum talaricum (tartary buckwheat) has previously been published. The recently reported genome sequence of tartary buckwheat provides theoretical basis for us to study and discuss the characteristics and expression of bZIP genes in tartary buckwheat based on the whole genome.

Results

In this study, 96 FtbZIP genes named from FtbZIP1 to FtbZIP96 were identified and divided into 11 subfamilies according to their genetic relationship with 70 bZIPs of A. thaliana. FtbZIP genes are not evenly distributed on the chromosomes, and we found tandem and segmental duplication events of FtbZIP genes on 8 tartary buckwheat chromosomes. According to the results of gene and motif composition, FtbZIP located in the same group contained analogous intron/exon organizations and motif composition. By qRT-PCR, we quantified the expression of FtbZIP members in stem, root, leaf, fruit, and flower and during fruit development. Exogenous ABA treatment increased the weight of tartary buckwheat fruit and changed the expressions of FtbZIP genes in group A.

Conclusions

Through our study, we identified 96 FtbZIP genes in tartary buckwheat and synthetically further analyzed the structure composition, evolution analysis and expression pattern of FtbZIP proteins. The expression pattern indicates that FtbZIP is important in the course of plant growth and development of tartary buckwheat. Through comprehensively analyzing fruit weight and FtbZIP genes expression after ABA treatment and endogenous ABA content of tartary buckwheat fruit, ABA may regulate downstream gene expression by regulating the expression of FtPinG0003523300.01 and FtPinG0003196200.01, thus indirectly affecting the fruit development of tartary buckwheat. This will help us to further study the function of FtbZIP genes in the tartary buckwheat growth and improve the fruit of tartary buckwheat.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5882-z) contains supplementary material, which is available to authorized users.

Stochasticity in transcriptional expression of a negative regulator of Arabidopsis ABA network

Stably heritable spatiotemporal co/over-expression of distinct transcriptional regulators across generations is a desired target as they signal traffic in the cell. Here, the stability and expression pattern of AtHB7 (Arabidopsis homeodomain-leucine zipper class I) cDNA was characterized in 220 random population of transformed tomato clones where no AtHB7 orthologous has been identified in to date. Integration of p35S::AtHB7 casette was tested by the amplification of the stretches (700/425 bp) in the target by NPT II/AtHB7 oligos. Transcriptional expression pattern for the amplicons of the specific transcripts in the leaf tissues of transformants were determined by qRT-PCR. Transgene copy number was negatively correlated with transgene expression level, yet a majority of transformants (78%) carried single-copy of transgene. About 1:3 of the lines containing two-copy inserts showed less transcript expression. Heterologous CaMV 35S promoter drove AtHB7, illuminated no penalty on transgene expression levels, stability or plant phenotype under drought stress. Integration and expression analysis of transcription factors is of great significance for reliable prediction of gene dosing/functions in plant genomes so as to sustain breeding under abiotic stress to guarantee food security.

Phylogenetic analysis of F-bZIP transcription factors indicates conservation of the zinc deficiency response across land plants

Basic leucine zipper (bZIP) transcription factors control important developmental and physiological processes in plants. In Arabidopsis thaliana, the three gene F-bZIP subfamily has been associated with zinc deficiency and salt stress response. Benefiting from the present abundance of plant genomic data, we performed an evolutionary and structural characterization of plant F-bZIPs. We observed divergence during seed plant evolution, into two groups and inferred different selective pressures for each. Group 1 contains AtbZIP19 and AtbZIP23 and appears more conserved, whereas Group 2, containing AtbZIP24, is more prone to gene loss and expansion events. Transcriptomic and experimental data reinforced AtbZIP19/23 as pivotal regulators of the zinc deficiency response, mostly via the activation of genes from the ZIP metal transporter family, and revealed that they are the main regulatory switch of AtZIP4. A survey of AtZIP4 orthologs promoters across different plant taxa revealed an enrichment of the Zinc Deficiency Response Element (ZDRE) to which both AtbZIP19/23 bind. Overall, our results indicate that while the AtbZIP24 function in the regulation of the salt stress response may be the result of neo-functionalization, the AtbZIP19/23 function in the regulation of the zinc deficiency response may be conserved in land plants (Embryophytes).

Promoter Polymorphism of Toll-Like Receptor 4 is Associated with a Decreased Risk of Coronary Artery Disease: A Case-Control Study in the Chinese Han Population

Background

Coronary artery disease (CAD) is considered a chronic inflammatory disease of the blood vessels. Toll-like receptor 4 (TLR4) is a transmembrane receptor involved in inflammatory reactions. The aim of this study was to determine the association between polymorphisms in the promoter region and 3′-untranslated region (3′-UTR) of TLR4, and the associated CAD risk.

Material/Methods

This study enrolled 424 participants with CAD and 424 controls without CAD. The polymorphisms in the promoter region and 3′-UTR of TLR4 were identified from the HapMap database, including rs10116253, rs10983755, and rs11536889. Genomic DNA was extracted from peripheral blood. Polymerase chain reaction-restriction fragment length polymorphism was performed to identify genotype polymorphisms. Relative luciferase activity was measured using the dual-luciferase reporter assay system.

Results

TLR4 rs10116253 in the promoter region was associated with CAD risk. The variant (CC+TC) genotypes of rs10116253 were associated with a decreased CAD risk (OR 95% CI 0.73 (0.54–0.98), p=0.034). In the stratification analyses, the variant (CC+TC) genotypes of rs10116253 were observed to have a relationship with decreased CAD risk in the male subgroup (OR: 95% CI 0.68 (0.48–0.98), p=0.041). Moreover, the variant CC and (CC+TC) genotypes of rs10116253 were correlated with a decreased CAD risk in participants younger than 60-year-old (TC: OR (95% CI 0.62 (0.39–0.98), p=0.042; TC+CC: OR 95% CI 0.63 (0.41–0.98), p=0.039). Regarding rs10116253, the luciferase activity of the mutant C allele construct was lower than that of the wild T allele construct (5.215±0.009 vs. 5.304±0.041; p=0.087).

Conclusions

The results provided evidence of an association between the TLR4 rs10116253 in the promoter region and a reduced risk of CAD.

Genome-wide analyses of the bZIP family reveal their involvement in the development, ripening and abiotic stress response in banana

The leucine zipper (bZIP) transcription factors play important roles in multiple biological processes. However, less information is available regarding the bZIP family in the important fruit crop banana. In this study, 121 bZIP transcription factor genes were identified in the banana genome. Phylogenetic analysis showed that MabZIPs were classified into 11 subfamilies. The majority of MabZIP genes in the same subfamily shared similar gene structures and conserved motifs. The comprehensive transcriptome analysis of two banana genotypes revealed the differential expression patterns of MabZIP genes in different organs, in various stages of fruit development and ripening, and in responses to abiotic stresses, including drought, cold, and salt. Interaction networks and co-expression assays showed that group A MabZIP-mediated networks participated in various stress signaling, which was strongly activated in Musa ABB Pisang Awak. This study provided new insights into the complicated transcriptional control of MabZIP genes and provided robust tissue-specific, development-dependent, and abiotic stress-responsive candidate MabZIP genes for potential applications in the genetic improvement of banana cultivars.

A novel wheat bZIP transcription factor, TabZIP60, confers multiple abiotic stress tolerances in transgenic Arabidopsis

The basic region/leucine zipper (bZIP) transcription factors (TFs) play vital roles in the response to abiotic stress. However, little is known about the function of bZIP genes in wheat abiotic stress. In this study, we report the isolation and functional characterization of the TabZIP60 gene. Three homologous genome sequences of TabZIP60 were isolated from hexaploid wheat and mapped to the wheat homoeologous group 6. A subcellular localization analysis indicated that TabZIP60 is a nuclear-localized protein that activates transcription. Furthermore, TabZIP60 gene transcripts were strongly induced by polyethylene glycol, salt, cold and exogenous abscisic acid (ABA) treatments. Further analysis showed that the overexpression of TabZIP60 in Arabidopsis resulted in significantly improved tolerances to drought, salt, freezing stresses and increased plant sensitivity to ABA in seedling growth. Meanwhile, the TabZIP60 was capable of binding ABA-responsive cis-elements that are present in promoters of many known ABA-responsive genes. A subsequent analysis showed that the overexpression of TabZIP60 led to enhanced expression levels of some stress-responsive genes and changes in several physiological parameters. Taken together, these results suggest that TabZIP60 enhances multiple abiotic stresses through the ABA signaling pathway and that modifications of its expression may improve multiple stress tolerances in crop plants.

Genetic variants in ABCA1 promoter affect transcription activity and plasma HDL level in pigs

Excess accumulation of cholesterol in plasma may result in coronary artery disease. Numerous studies have demonstrated that ATP-binding cassette protein A1 (ABCA1) mediates the efflux of cholesterol and phospholipids to apolipoproteins, a process necessary for plasma high density lipoprotein (HDL) formation. Higher plasma levels of HDL are associated with lower risk for cardiovascular disease. Studies of human disease and animal models had shown that an increased hepatic ABCA1 activity relates to an enhanced plasma HDL level. In this study, we hypothesized that functional mutations in the ABCA1 promoter in pigs may affect gene transcription activity, and consequently the HDL level in plasma. The promoter region of ABCA1 was comparatively scanned by direct sequencing with pool DNA of high- and low-HDL groups (n=30 for each group). Two polymorphisms, c. - 608A>G and c. - 418T>A, were revealed with reverse allele distribution in the two groups. The two polymorphisms were completely linked and formed only G-A or A-T haplotypes when genotyped in a larger population (n=526). Furthermore, we found that the G-A/G-A genotype was associated with higher HDL and ABCA1 mRNA level than A-T/A-T genotype. Luciferase assay also revealed that G-A haplotype promoter had higher activity than A-T haplotype. Single-nucleotide mutant assay showed that c.-418T>A was the causal mutation for ABCA1 transcription activity alteration. Conclusively, we identified two completely linked SNPs in porcine ABCA1 promoter region which have influence on the plasma HDL level by altering ABCA1 gene transcriptional activity.

Genomic organization and promoter analysis of a transcriptional repressor gene from Fenneropenaeus chinensis

In this study, we cloned and sequenced genomic sequences from a Fenneropenaeus chinensis transcriptional repressor gene, FcTR. The FcTR gene is 2,671 bp in length and has four exons and three introns. The 873 bp promoter contains several transcription factor binding sites, including a TATA box and a cyclic AMP-responsive element. Promoter deletion analysis using a luciferase reporter gene identified regulatory elements. Challenge with white spot syndrome virus increased expression from the promoter-deletion constructs. These results suggest that FcTR might play an important role in the shrimp immune response.

Genome-wide analysis of the bZIP transcription factors in cucumber

bZIP proteins are one of the largest transcriptional regulators playing crucial roles in plant development, physiological processes, and biotic/abiotic stress responses. Despite the availability of recently published draft genome sequence of Cucumis sativus, no comprehensive investigation of these family members has been presented for cucumber. We have identified 64 bZIP transcription factor-encoding genes in the cucumber genome. Based on structural features of their encoded proteins, CsbZIP genes could be classified into 6 groups. Cucumber bZIP genes were expanded mainly by segmental duplication rather than tandem duplication. Although segmental duplication rate of the CsbZIP genes was lower than that of Arabidopsis, rice and sorghum, it was observed as a common expansion mechanism. Some orthologous relationships and chromosomal rearrangements were observed according to comparative mapping analysis with other species. Genome-wide expression analysis of bZIP genes indicated that 64 CsbZIP genes were differentially expressed in at least one of the ten sampled tissues. A total of 4 CsbZIP genes displayed higher expression values in leaf, flowers and root tissues. The in silico micro-RNA (miRNA) and target transcript analyses identified that a total of 21 CsbZIP genes were targeted by 38 plant miRNAs. CsbZIP20 and CsbZIP22 are the most targeted by miR165 and miR166 family members, respectively. We also analyzed the expression of ten CsbZIP genes in the root and leaf tissues of drought-stressed cucumber using quantitative RT-PCR. All of the selected CsbZIP genes were measured as increased in root tissue at 24th h upon PEG treatment. Contrarily, the down-regulation was observed in leaf tissues of all analyzed CsbZIP genes. CsbZIP12 and CsbZIP44 genes showed gradual induction of expression in root tissues during time points. This genome-wide identification and expression profiling provides new opportunities for cloning and functional analyses, which may be used in further studies for improving stress tolerance in plants.

TdIF1 recognizes a specific DNA sequence through its Helix-Turn-Helix and AT-hook motifs to regulate gene transcription

TdIF1 was originally identified as a protein that directly binds to DNA polymerase TdT. TdIF1 is also thought to function in transcription regulation, because it binds directly to the transcriptional factor TReP-132, and to histone deacetylases HDAC1 and HDAC2. Here we show that TdIF1 recognizes a specific DNA sequence and regulates gene transcription. By constructing TdIF1 mutants, we identify amino acid residues essential for its interaction with DNA. An in vitro DNA selection assay, SELEX, reveals that TdIF1 preferentially binds to the sequence 5′-GNTGCATG-3′ following an AT-tract, through its Helix-Turn-Helix and AT-hook motifs. We show that four repeats of this recognition sequence allow TdIF1 to regulate gene transcription in a plasmid-based luciferase reporter assay. We demonstrate that TdIF1 associates with the RAB20 promoter, and RAB20 gene transcription is reduced in TdIF1-knocked-down cells, suggesting that TdIF1 stimulates RAB20 gene transcription.

Genome-wide expansion and expression divergence of the basic leucine zipper transcription factors in higher plants with an emphasis on sorghum

Plant bZIP transcription factors play crucial roles in multiple biological processes. However, little is known about the sorghum bZIP gene family although the sorghum genome has been completely sequenced. In this study, we have carried out a genome-wide identification and characterization of this gene family in sorghum. Our data show that the genome encodes at least 92 bZIP transcription factors. These bZIP genes have been expanded mainly by segmental duplication. Such an expansion mechanism has also been observed in rice, arabidopsis and many other plant organisms, suggesting a common expansion mode of this gene family in plants. Further investigation shows that most of the bZIP members have been present in the most recent common ancestor of sorghum and rice and the major expansion would occur before the sorghum-rice split era. Although these bZIP genes have been duplicated with a long history, they exhibited limited functional divergence as shown by nonsynonymous substitutions (Ka)/synonymous substitutions (Ks) analyses. Their retention was mainly due to the high percentages of expression divergence. Our data also showed that this gene family might play a role in multiple developmental stages and tissues and might be regarded as important regulators of various abiotic stresses and sugar signaling.

Assessing energetic contributions to binding from a disordered region in a protein-protein interaction

Many functional proteins are at least partially disordered prior to binding. Although the structural transitions upon binding of disordered protein regions can influence the affinity and specificity of protein complexes, their precise energetic contributions to binding are unknown. Here, we use a model protein-protein interaction system in which a locally disordered region has been modified by directed evolution to quantitatively assess the thermodynamic and structural contributions to binding of disorder-to-order transitions. Through X-ray structure determination of the protein binding partners before and after complex formation and isothermal titration calorimetry of the interactions, we observe a correlation between protein ordering and binding affinity for complexes along this affinity maturation pathway. Additionally, we show that discrepancies between observed and calculated heat capacities based on buried surface area changes in the protein complexes can be explained largely by heat capacity changes that would result solely from folding the locally disordered region. Previously developed algorithms for predicting binding energies of protein-protein interactions, however, are unable to correctly model the energetic contributions of the structural transitions in our model system. While this highlights the shortcomings of current computational methods in modeling conformational flexibility, it suggests that the experimental methods used here could provide training sets of molecular interactions for improving these algorithms and further rationalizing molecular recognition in protein-protein interactions.

A strained DNA binding helix is conserved for site recognition, folding nucleation, and conformational modulation

Nucleic acid recognition is often mediated by alpha-helices or disordered regions that fold into alpha-helix on binding. A peptide bearing the DNA recognition helix of HPV16 E2 displays type II polyproline (PII) structure as judged by pH, temperature, and solvent effects on the CD spectra. NMR experiments indicate that the canonical alpha-helix is stabilized at the N-terminus, while the PII forms at the C-terminus half of the peptide. Re-examination of the dihedral angles of the DNA binding helix in the crystal structure and analysis of the NMR chemical shift indexes confirm that the N-terminus half is a canonical alpha-helix, while the C-terminal half adopts a 3(10) helix structure. These regions precisely match two locally driven folding nucleii, which partake in the native hydrophobic core and modulate a conformational switch in the DNA binding helix. The peptide shows only weak and unspecific residual DNA binding, 10(4)-fold lower affinity, and 500-fold lower discrimination capacity compared with the domain. Thus, the precise side chain conformation required for modulated and tight physiological binding by HPV E2 is largely determined by the noncanonical strained alpha-helix conformation, "presented" by this unique architecture. (c) 2009 Wiley Periodicals, Inc. Biopolymers 91: 432-443, 2009.

The role of bZIP transcription factors in green plant evolution: adaptive features emerging from four founder genes

Background

Transcription factors of the basic leucine zipper (bZIP) family control important processes in all eukaryotes. In plants, bZIPs are regulators of many central developmental and physiological processes including photomorphogenesis, leaf and seed formation, energy homeostasis, and abiotic and biotic stress responses. Here we performed a comprehensive phylogenetic analysis of bZIP genes from algae, mosses, ferns, gymnosperms and angiosperms.

Methodology/Principal Findings

We identified 13 groups of bZIP homologues in angiosperms, three more than known before, that represent 34 Possible Groups of Orthologues (PoGOs). The 34 PoGOs may correspond to the complete set of ancestral angiosperm bZIP genes that participated in the diversification of flowering plants. Homologous genes dedicated to seed-related processes and ABA-mediated stress responses originated in the common ancestor of seed plants, and three groups of homologues emerged in the angiosperm lineage, of which one group plays a role in optimizing the use of energy.

Conclusions/Significance

Our data suggest that the ancestor of green plants possessed four bZIP genes functionally involved in oxidative stress and unfolded protein responses that are bZIP-mediated processes in all eukaryotes, but also in light-dependent regulations. The four founder genes amplified and diverged significantly, generating traits that benefited the colonization of new environments.

Safety assessment considerations for food and feed derived from plants with genetic modifications that modulate endogenous gene expression and pathways

The current globally recognized comparative food and feed safety assessment paradigm for biotechnology-derived crops is a robust and comprehensive approach for evaluating the safety of both the inserted gene product and the resulting crop. Incorporating many basic concepts from food safety, toxicology, nutrition, molecular biology, and plant breeding, this approach has been used effectively by scientists and regulatory agencies for 10-15 years. Current and future challenges in agriculture include the need for improved yields, tolerance to biotic and abiotic stresses, and improved nutrition. The next generation of biotechnology-derived crops may utilize regulatory proteins, such as transcription factors that modulate gene expression and/or endogenous plant pathways. In this review, we discuss the applicability of the current safety assessment paradigm to biotechnology-derived crops developed using modifications involving regulatory proteins. The growing literature describing the molecular biology underlying plant domestication and conventional breeding demonstrates the naturally occurring genetic variation found in plants, including significant variation in the classes, expression, and activity of regulatory proteins. Specific examples of plant modifications involving insertion or altered expression of regulatory proteins are discussed as illustrative case studies supporting the conclusion that the current comparative safety assessment process is appropriate for these types of biotechnology-developed crops.

Molecular mechanism of interaction between oligopeptides and double-stranded DNA

Formation of Ala-Glu-Asp-Gly peptide complex with duplex DNA [poly(dA-dT):poly(dA-dT)] at neutral pH and ambient temperature was studied experimentally. Peptide binding to duplex DNA was associated with hyperchromic effect indicating local separation of the double helix. Energies of paired interactions between amino acid residues and nucleotide bases were compared with the energy of bond between two bases in the nucleotide pair (dA-dT). A new concept of interaction between two types of information-carrying molecules is put forward. This interaction underlies triggering of protein synthesis mechanism and can explain the emergence of life on the Earth.

Exploring and designing protein function with restricted diversity

Combinatorial libraries with restricted diversity can be used to rapidly map binding energetics across protein interfaces. Shotgun scanning strategies have been used for alanine scanning and for alternative mutagenesis schemes that provide high-resolution functional views of binding interfaces. In addition, synthetic antibodies have been derived from naïve libraries restricted to a binary code to explore the minimal requirements for molecular recognition. These studies shed light on the underlying principles governing molecular recognition, and provide rapid yet quantitative alternatives to conventional biophysical methods for exploring protein structure and function.

Long-range cooperative binding effects in a T cell receptor variable domain

Although cellular processes depend on protein-protein interactions, our understanding of molecular recognition between proteins remains far from comprehensive. Protein-protein interfaces are structural and energetic mosaics in which a subset of interfacial residues, called hot spots, contributes disproportionately to the affinity of the complex. These hot-spot residues can be further clustered into hot regions. It has been proposed that binding energetics between residues within a hot region are cooperative, whereas those between hot regions are strictly additive. If this idea held true for all protein-protein interactions, then energetically significant long-range conformational effects would be unlikely to occur. In the present study, we show cooperative binding energetics between distinct hot regions that are separated by >20 A. Using combinatorial mutagenesis and surface plasmon resonance binding analysis to dissect additivity and cooperativity in a complex formed between a variable domain of a T cell receptor and a bacterial superantigen, we find that combinations of mutations from each of two hot regions exhibited significant cooperative energetics. Their connecting sequence is composed primarily of a single beta-strand of the T cell receptor variable Ig domain, which has been observed to undergo a strand-switching event and does not form an integral part of the stabilizing core of this Ig domain. We propose that these cooperative effects are propagated through a dynamic structural network. Cooperativity between hot regions has significant implications for the prediction and inhibition of protein-protein interactions.

Human transcription factors contain a high fraction of intrinsically disordered regions essential for transcriptional regulation

Human transcriptional regulation factors, such as activators, repressors, and enhancer-binding factors are quite different from their prokaryotic counterparts in two respects: the average sequence in human is more than twice as long as that in prokaryotes, while the fraction of sequence aligned to domains of known structure is 31% in human transcription factors (TFs), less than half of that in bacterial TFs (72%). Intrinsically disordered (ID) regions were identified by a disorder-prediction program, and were found to be in good agreement with available experimental data. Analysis of 401 human TFs with experimental evidence from the Swiss-Prot database showed that as high as 49% of the entire sequence of human TFs is occupied by ID regions. More than half of the human TFs consist of a small DNA binding domain (DBD) and long ID regions frequently sandwiching unassigned regions. The remaining TFs have structural domains in addition to DBDs and ID regions. Experimental studies, particularly those with NMR, revealed that the transactivation domains in unbound TFs are usually unstructured, but become structured upon binding to their partners. The sequences of human and mouse TF orthologues are 90.5% identical despite a high incidence of ID regions, probably reflecting important functional roles played by ID regions. In general ID regions occupy a high fraction in TFs of eukaryotes, but not in prokaryotes. Implications of this dichotomy are discussed in connection with their functional roles in transcriptional regulation and evolution.

Structural basis of affinity maturation and intramolecular cooperativity in a protein-protein interaction

Although protein-protein interactions are involved in nearly all cellular processes, general rules for describing affinity and selectivity in protein-protein complexes are lacking, primarily because correlations between changes in protein structure and binding energetics have not been well determined. Here, we establish the structural basis of affinity maturation for a protein-protein interaction system that we had previously characterized energetically. This model system exhibits a 1500-fold affinity increase. Also, its affinity maturation is restricted by negative intramolecular cooperativity. With three complex and six unliganded variant X-ray crystal structures, we provide molecular snapshots of protein interface remodeling events that span the breadth of the affinity maturation process and present a comprehensive structural view of affinity maturation. Correlating crystallographically observed structural changes with measured energetic changes reveals molecular bases for affinity maturation, intramolecular cooperativity, and context-dependent binding.

Alternative views of functional protein binding epitopes obtained by combinatorial shotgun scanning mutagenesis

Combinatorial shotgun scanning mutagenesis was used to analyze two large, related protein binding sites to assess the specificity and importance of individual side chain contributions to binding affinity. The strategy allowed for cost-effective generation of a plethora of functional data. The ease of the technology promoted comprehensive investigations, in which the classic alanine-scanning approach was expanded with two additional strategies, serine- and homolog-scanning. Binding of human growth hormone (hGH) to the hGH receptor served as the model system. The entire high affinity receptor-binding sites (site 1) of wild-type hGH (hGHwt) and of an affinity-improved variant (hGHv) were investigated and the results were compared. The contributions that 35 residue positions make to binding were assessed on each hormone molecule by both serine- and homolog-scanning. The hormone molecules were displayed on the surfaces of bacteriophage, and the 35 positions were randomized simultaneously to allow equal starting frequencies of the wild-type residue and either serine or a homologous mutation in separate libraries. Functional selections for binding to the hGH receptor shifted the relative wild-type/mutant frequencies at each position to an extent characteristic of the functional importance of the side chain. Functional epitope maps were created and compared to previous maps obtained by alanine-scanning. Comparisons between the different scans provide insights into the affinity maturation process that produced hGHv. The serine and homolog-scanning results expand upon and complement the alanine-scanning results and provide additional data on the robustness of the high affinity receptor-binding site of hGH.

Intramolecular cooperativity in a protein binding site assessed by combinatorial shotgun scanning mutagenesis

Combinatorial shotgun alanine-scanning was used to assess intramolecular cooperativity in the high affinity site (site 1) of human growth hormone (hGH) for binding to its receptor. A total of 19 side-chains were analyzed and statistically significant data were obtained for 145 of the 171 side-chain pairs. The analysis revealed that 90% of the side-chain pairs exhibited no statistically significant pair interactions, and the remaining 10% of side-chain pairs exhibited only small interactions corresponding to cooperative interaction energies with magnitudes less than 0.4 kcal/mol. The statistical predictions were tested by measuring affinities for purified mutant proteins and were found to be accurate for five of six side-chain pairs tested. The results reveal that hGH site 1 behaves in a highly additive manner and suggest that shotgun scanning should be useful for assessing cooperative effects in other protein-protein interactions.

Identification of the 'NORE' (N-Oct-3 responsive element), a novel structural motif and composite element

N-Oct-3 is a neuronal transcription factor widely expressed in the developing mammalian central nervous system, and necessary to maintain neural cell differentiation. The key role of N-Oct-3 in the transcriptional regulation of a multiplicity of genes is primarily due to the structural plasticity of its so-called ‘POU’ (acronym of Pit, Oct, Unc) DNA-binding domain. We have recently reported about the unusual dual neuro-specific transcriptional regulation displayed by N-Oct-3 [Blaud,M., Vossen,C., Joseph,G., Alazard,R., Erard,M. and Nieto,L. (2004) J. Mol. Biol., 339, 1049–1058]. To elucidate the underlying molecular mechanisms, we have now made use of molecular modeling, DNA footprinting and electrophoretic mobility shift assay techniques. This combined approach has allowed us to uncover a novel mode of homodimerization adopted by the N-Oct-3 POU domain bound to the neuronal aromatic amino acids de-carboxylase and corticotropin-releasing hormone gene promoters and to demonstrate that this pattern is induced by a structural motif that we have termed ‘NORE’ (N-Oct-3 responsive element), comprising the 14 bp sequence element TNNRTAAATAATRN. In addition, we have been able to explain how the same structural motif can also induce the formation of a heterodimer in association with hepatocyte nuclear factor 3β(/Forkhead box a2). Finally, we discuss the possible role of the NORE motif in relation to neuroendocrine lung tumor formation, and in particular the development of small cell lung cancer.

A genome-wide and nonredundant mouse transcription factor database

Here we describe the development of a genome-wide and nonredundant mouse transcription factor database and its viewer (http://genome.gsc.riken.gp/TFdb/). We systematically selected transcription factors with DNA-binding properties and their regulators on the basis of their LocusLink and Gene Ontology annotations. We also incorporated into our database information regarding the corresponding available cDNA clones and their structural properties. Because of these features, our database is unique and may provide useful information for systematic genome-wide studies of transcriptional regulation.

Unique gene expression program of human germinal center T helper cells

Gene expression profiling was used to compare the gene expression patterns of human germinal center (GC) T helper (Th) cells with other CD4+ T-cell subsets (naive, central, and effector memory T cells). GC-Th cells, specifically localized in germinal centers to help B cells, are distantly related to central and effector memory T cells in global gene expression profiles. GC-Th cells displayed substantial differences in mRNA for adhesion molecules, chemoattractant receptors, and cytokines compared with other populations. Distinct expression of transcriptional factors by GC-Th cells is consistent with the hypothesis that they may be different from other T cells in cell lineage. Interestingly, CXCL13, a critical chemokine for B-cell entry to lymphoid follicles, is one of the most highly up-regulated genes in GC-Th cells. GC-Th cells (but not other T cells) produce and secrete large amounts of functional CXCL13 upon T-cell receptor activation, a process that is dependent on costimulation, requires translation and transcription, and is dramatically enhanced by activation in the presence of GC-B cells. This study revealed for the first time the unique gene expression program of GC-Th cells.

Structure and DNA-binding Sites of the SWI1 AT-rich Interaction Domain (ARID) Suggest Determinants for Sequence-specific DNA Recognition

ARID (AT-rich interaction domain) is a homologous family of DNA-binding domains that occur in DNA-binding proteins from a wide variety of species, ranging from yeast to nematodes, insects, mammals, and plants. SWI1, a member of the SWI/SNF protein complex that is involved in chromatin remodeling during transcription, contains the ARID motif. The ARID domain of human SWI1 (also known as p270) does not select for a specific DNA sequence from a random sequence pool. The lack of sequence specificity shown by the SWI1 ARID domain stands in contrast to the other characterized ARID domains, which recognize specific AT-rich sequences. We have solved the three-dimensional structure of human SWI1 ARID using solution NMR methods. In addition, we have characterized nonspecific DNA binding by the SWI1 ARID domain. Results from this study indicate that a flexible, long, internal loop in the ARID motif is likely to be important for sequence-specific DNA recognition. The structure of the human SWI1 ARID domain also represents a distinct structural subfamily. Studies of ARID indicate that the boundary of DNA binding structural and functional domains can extend beyond the sequence homologous region in a homologous family of proteins. Structural studies of homologous domains such as the ARID family of DNA-binding domains should provide information to better predict the boundary of structural and functional domains in structural genomic studies.

Magnetic resonance image-guided proteomics of human glioblastoma multiforme

PURPOSE

To investigate the correlation between gadolinium contrast-enhancement patterns on T1-weighted magnetic resonance (MR) images and spatial changes in protein expression profiles in human glioblastoma multiforme (GBM) and the use of imaging as a noninvasive technique to evaluate the heterogeneity of solid tumors prior to microarray analysis.

MATERIALS AND METHODS

Four patients with MR images and confirmed diagnosis of GBM were enrolled in the study. Intraoperative stereotaxy was used in conjunction with MR images to identify contrast-enhanced (CE) and nonenhanced (NE) regions of the tumor during surgical resection. Total protein was extracted from resected tumor samples using standard techniques and subjected to proteomic analysis using surface enhanced laser desorption/ionization time of flight mass spectrometry (SELDI-TOF-MS).

RESULTS

We found that protein profiles from CE and NE regions within a given tumor have qualitative and semiquantitative proteomic pattern differences, suggesting an altered gene expression profile that correlates with detectable tissue imaging parameters. We also found that CE regions within the same tumor exhibited distinct differences in protein expression profiles, despite similar histological features. In addition, there were marked similarities in the proteomic patterns among the NE regions across all patients, while the CE regions were distinct, suggesting that the CE regions have complex protein profiles unique to individuals.

CONCLUSION

The results demonstrate that major differences in protein expression patterns within a tumor can be correlated to radiographic findings. Image-guided proteomics holds promise for characterizing tissue prior to microarray analysis designed to identify specific diagnostic markers and therapeutic targets within solid tumors.

Small Molecules that Mimic Components of Bioactive Protein Surfaces

Small molecules designed to mimic specific structural components of a protein (peptide strands, sheets, turns, helices, or amino acids) can be expected to display agonist or antagonist biological responses by virtue of interacting with the same receptors that recognize the protein. Here we describe some minimalist approaches to structural mimetics of amino acids and of strand, turn, or helix segments of proteins. The designed molecules show potent and selective inhibition of protease, transferase, and phospholipase enzymes, or antagonism of G-protein coupled or transcriptional receptors, and have potent anti-tumour, anti-inflammatory, or antiviral activity.

Common domains in the initiators of DNA replication in Bacteria, Archaea and Eukarya: combined structural, functional and phylogenetic perspectives

Although DNA replication is the universal process for the transmission of genetic information in all living organisms, until very recently evidence was lacking for a related structure and function in the proteins (initiators) that trigger replication in the three 'Life Domains' (Bacteria, Archaea and Eukarya). In this article new data concerning the presence of common features in the initiators of chromosomal replication in bacteria, archaea and eukaryotes are reviewed. Initiators are discussed in the light of: (i) The structure and function of their conserved ATPases Associated with various cellular Activities (AAA+) and winged-helix domains. (ii) The nature of the macromolecular assemblies that they constitute at the replication origins. (iii) Their possible phylogenetic relationship, attempting to sketch the essentials of a hypothetical DNA replication initiator in the micro-organism proposed to be the ancestor of all living cells.