Functionally specified protein signatures distinctive for each of the different blue copper proteins

Molecular evolution of Phytocyanin gene and analysis of expression at different coloring periods in apple (Malus domestica)

BACKGROUND

PC (phytocyanin) is a class of copper-containing electron transfer proteins closely related to plant photosynthesis, abiotic stress responses growth and development in plants, and regulation of the expression of some flavonoids and phenylpropanoids, etc., however, compared with other plants, the PC gene family has not been systematically characterized in apple.

RESULTS

A total of 59 MdPC gene members unevenly distributed across 12 chromosomes were identified at the genome-wide level. The proteins of the MdPC family were classified into four subfamilies based on differences in copper binding sites and glycosylation sites: Apple Early nodulin-like proteins (MdENODLs), Apple Uclacyanin-like proteins (MdUCLs), Apple Stellacyanin-like proteins (MdSCLs), and Apple Plantacyanin-like proteins (MdPLCLs). Some MdPC members with similar gene structures and conserved motifs belong to the same group or subfamily. The internal collinearity analysis revealed 14 collinearity gene pairs among members of the apple MdPC gene. Interspecific collinearity analysis showed that apple had 31 and 35 homologous gene pairs with strawberry and grape, respectively. Selection pressure analysis indicated that the MdPC gene was under purifying selection. Prediction of protein interactions showed that MdPC family members interacted strongly with the Nad3 protein. GO annotation results indicated that the MdPC gene also regulated the biosynthesis of phenylpropanoids. Chip data analysis showed that (MdSCL3, MdSCL7 and MdENODL27) were highly expressed in mature fruits and peels. Many cis-regulatory elements related to light response, phytohormones, abiotic stresses and flavonoid biosynthetic genes regulation were identified 2000 bp upstream of the promoter of the MdPC gene, and qRT-PCR results showed that gene members in Group IV (MdSCL1/3, MdENODL27) were up-regulated at all five stages of apple coloring, but the highest expression was observed at the DAF13 (day after fruit bag removal) stage. The gene members in Group II (MdUCL9, MdPLCL3) showed down-regulated or lower expression in the first four stages of apple coloring but up-regulated and highest expression in the DAF 21 stage.

CONCLUSION

Herein, one objective of these findings is to provide valuable information for understanding the structure, molecular evolution, and expression pattern of the MdPC gene, another major objective in this study was designed to lay the groundwork for further research on the molecular mechanism of PC gene regulation of apple fruit coloration.

The PIF1-miR408-PLANTACYANIN repression cascade regulates light-dependent seed germination

Light-dependent seed germination is a vital process for many seed plants. A decisive event in light-induced germination is degradation of the central repressor PHYTOCHROME INTERACTING FACTOR 1 (PIF1). The balance between gibberellic acid (GA) and abscisic acid (ABA) helps to control germination. However, the cellular mechanisms linking PIF1 turnover to hormonal balancing remain elusive. Here, employing far-red light-induced Arabidopsis thaliana seed germination as the experimental system, we identified PLANTACYANIN (PCY) as an inhibitor of germination. It is a blue copper protein associated with the vacuole that is both highly expressed in mature seeds and rapidly silenced during germination. Molecular analyses showed that PIF1 binds to the miR408 promoter and represses miR408 accumulation. This in turn posttranscriptionally modulates PCY abundance, forming the PIF1-miR408-PCY repression cascade for translating PIF1 turnover to PCY turnover during early germination. Genetic analysis, RNA-sequencing, and hormone quantification revealed that PCY is necessary and sufficient to maintain the PIF1-mediated seed transcriptome and the low-GA-high-ABA state. Furthermore, we found that PCY domain organization and regulation by miR408 are conserved features in seed plants. These results revealed a cellular mechanism whereby PIF1-relayed external light signals are converted through PCY turnover to internal hormonal profiles for controlling seed germination.

Genome-wide identification, classification, and expression of phytocyanins in Populus trichocarpa

74 phytocyanin genes were identified in the Populus trichocarpa genome. Phylogenetic analysis grouped the PC proteins into four subfamilies (UCs, PLCs, SCs, and ENODLs). Closely related PC proteins share similar motifs, implying similar functions. Expression profiles of PtPC genes were analyzed in response to drought and salt-stress. Phytocyanins (PCs) are blue copper proteins associated with electron carrier activity that have a large influence on plant growth and resistance. The majority of PCs are chimeric arabinogalactan proteins (AGPs). In this work, we identified 74 PC genes in Populus trichocarpa and analyzed them comprehensively. Based on the ligands composition of copper-binding sites, glycosylation state, the domain structure and spectral characteristics of PC genes, PCs were divided into four subfamilies [uclacyanins (UCs), plantacyanins (PLCs), stellacyanins (SCs) and early nodulin-like proteins (ENODLs)], and phylogenetic relationship analysis classified them into seven groups. All PtPCs are randomly distributed on 17 of the 19 poplar chromosomes, and they appear to have undergone expansion via segmental duplication. Eight PtPCs do not contain introns, and each group has a similar conserved motif structure. Promoter analysis revealed cis-elements related to growth, development and stress responses, and established orthology relationships of PCs between Arabidopsis and poplar by synteny analysis. Expression profile analysis and qRT-PCR analysis showed that PtPCs were expressed widely in various tissues. Quantitative real-time RT-PCR analysis of PC genes expression in response to salt and drought stress revealed their stress-responses profiles. This work provides a theoretical basis for a further study of stress resistance mechanisms and the function of PC genes in poplar growth and development.

Comparative analysis of the phytocyanin gene family in 10 plant species: a focus on Zea mays

Phytocyanins (PCs) are plant-specific blue copper proteins, which play essential roles in electron transport. While the origin and expansion of this gene family is not well-investigated in plants. Here, we investigated their evolution by undertaking a genome-wide identification and comparison in 10 plants: Arabidopsis, rice, poplar, tomato, soybean, grape, maize, Selaginella moellendorffii, Physcomitrella patens, and Chlamydomonas reinhardtii. We found an expansion process of this gene family in evolution. Except PCs in Arabidopsis and rice, which have described in previous researches, a structural analysis of PCs in other eight plants indicated that 292 PCs contained N-terminal secretion signals and 217 PCs were expected to have glycosylphosphatidylinositol-anchor signals. Moreover, 281 PCs had putative arabinogalactan glycomodules and might be AGPs. Chromosomal distribution and duplication patterns indicated that tandem and segmental duplication played dominant roles for the expansion of PC genes. In addition, gene organization and motif compositions are highly conserved in each clade. Furthermore, expression profiles of maize PC genes revealed diversity in various stages of development. Moreover, all nine detected maize PC genes (ZmUC10, ZmUC16, ZmUC19, ZmSC2, ZmUC21, ZmENODL10, ZmUC22, ZmENODL13, and ZmENODL15) were down-regulated under salt treatment, and five PCs (ZmUC19, ZmSC2, ZmENODL10, ZmUC22, and ZmENODL13) were down-regulated under drought treatment. ZmUC16 was strongly expressed after drought treatment. This study will provide a basis for future understanding the characterization of this family.

Comparative genomics in acid mine drainage biofilm communities reveals metabolic and structural differentiation of co-occurring archaea

Background

Metal sulfide mineral dissolution during bioleaching and acid mine drainage (AMD) formation creates an environment that is inhospitable to most life. Despite dominance by a small number of bacteria, AMD microbial biofilm communities contain a notable variety of coexisting and closely related Euryarchaea, most of which have defied cultivation efforts. For this reason, we used metagenomics to analyze variation in gene content that may contribute to niche differentiation among co-occurring AMD archaea. Our analyses targeted members of the Thermoplasmatales and related archaea. These results greatly expand genomic information available for this archaeal order.

Results

We reconstructed near-complete genomes for uncultivated, relatively low abundance organisms A-, E-, and Gplasma, members of Thermoplasmatales order, and for a novel organism, Iplasma. Genomic analyses of these organisms, as well as Ferroplasma type I and II, reveal that all are facultative aerobic heterotrophs with the ability to use many of the same carbon substrates, including methanol. Most of the genomes share genes for toxic metal resistance and surface-layer production. Only Aplasma and Eplasma have a full suite of flagellar genes whereas all but the Ferroplasma spp. have genes for pili production. Cryogenic-electron microscopy (cryo-EM) and tomography (cryo-ET) strengthen these metagenomics-based ultrastructural predictions. Notably, only Aplasma, Gplasma and the Ferroplasma spp. have predicted iron oxidation genes and Eplasma and Iplasma lack most genes for cobalamin, valine, (iso)leucine and histidine synthesis.

Conclusion

The Thermoplasmatales AMD archaea share a large number of metabolic capabilities. All of the uncultivated organisms studied here (A-, E-, G-, and Iplasma) are metabolically very similar to characterized Ferroplasma spp., differentiating themselves mainly in their genetic capabilities for biosynthesis, motility, and possibly iron oxidation. These results indicate that subtle, but important genomic differences, coupled with unknown differences in gene expression, distinguish these organisms enough to allow for co-existence. Overall this study reveals shared features of organisms from the Thermoplasmatales lineage and provides new insights into the functioning of AMD communities.

DLocalMotif: a discriminative approach for discovering local motifs in protein sequences

MOTIVATION

Local motifs are patterns of DNA or protein sequences that occur within a sequence interval relative to a biologically defined anchor or landmark. Current protein motif discovery methods do not adequately consider such constraints to identify biologically significant motifs that are only weakly over-represented but spatially confined. Using negatives, i.e. sequences known to not contain a local motif, can further increase the specificity of their discovery.

RESULTS

This article introduces the method DLocalMotif that makes use of positional information and negative data for local motif discovery in protein sequences. DLocalMotif combines three scoring functions, measuring degrees of motif over-representation, entropy and spatial confinement, specifically designed to discriminatively exploit the availability of negative data. The method is shown to outperform current methods that use only a subset of these motif characteristics. We apply the method to several biological datasets. The analysis of peroxisomal targeting signals uncovers several novel motifs that occur immediately upstream of the dominant peroxisomal targeting signal-1 signal. The analysis of proline-tyrosine nuclear localization signals uncovers multiple novel motifs that overlap with C2H2 zinc finger domains. We also evaluate the method on classical nuclear localization signals and endoplasmic reticulum retention signals and find that DLocalMotif successfully recovers biologically relevant sequence properties.

AVAILABILITY

http://bioinf.scmb.uq.edu.au/dlocalmotif/

Terminal oxidase diversity and function in "Metallosphaera yellowstonensis": gene expression and protein modeling suggest mechanisms of Fe(II) oxidation in the sulfolobales

"Metallosphaera yellowstonensis" is a thermoacidophilic archaeon isolated from Yellowstone National Park that is capable of autotrophic growth using Fe(II), elemental S, or pyrite as electron donors. Analysis of the draft genome sequence from M. yellowstonensis strain MK1 revealed seven different copies of heme copper oxidases (subunit I) in a total of five different terminal oxidase complexes, including doxBCEF, foxABCDEFGHIJ, soxABC, and the soxM supercomplex, as well as a novel hypothetical two-protein doxB-like polyferredoxin complex. Other genes found in M. yellowstonensis with possible roles in S and or Fe cycling include a thiosulfate oxidase (tqoAB), a sulfite oxidase (som), a cbsA cytochrome b(558/566), several small blue copper proteins, and a novel gene sequence coding for a putative multicopper oxidase (Mco). Results from gene expression studies, including reverse transcriptase (RT) quantitative PCR (qPCR) of cultures grown autotrophically on either Fe(II), pyrite, or elemental S showed that the fox gene cluster and mco are highly expressed under conditions where Fe(II) is an electron donor. Metagenome sequence and gene expression studies of Fe-oxide mats confirmed the importance of fox genes (e.g., foxA and foxC) and mco under Fe(II)-oxidizing conditions. Protein modeling of FoxC suggests a novel lysine-lysine or lysine-arginine heme B binding domain, indicating that it is likely the cytochrome component of a heterodimer complex with foxG as a ferredoxin subunit. Analysis of mco shows that it encodes a novel multicopper blue protein with two plastocyanin type I copper domains that may play a role in the transfer of electrons within the Fox protein complex. An understanding of metabolic pathways involved in aerobic iron and sulfur oxidation in Sulfolobales has broad implications for understanding the evolution and niche diversification of these thermophiles as well as practical applications in fields such as bioleaching of trace metals from pyritic ores.

Cotton BCP genes encoding putative blue copper-binding proteins are functionally expressed in fiber development and involved in response to high-salinity and heavy metal stresses

Copper is vitally required for plants at low concentrations but extremely toxic for plants at elevated concentrations. Plants have evolved a series of mechanisms to prevent the consequences of the excess or deficit of copper. These mechanisms require copper-interacting proteins involved in copper trafficking. Blue copper-binding proteins (BCPs) are a class of copper proteins containing one blue copper-binding domain binding a single type I copper. To investigate the role of BCPs in plant development and in response to stresses, we isolated nine cDNAs encoding the putative blue copper-binding proteins (GhBCPs) from cotton (Gossypium hirsutum). Meanwhile, four corresponding genes (including GhBCP1-GhBCP4), which contain a single intron inserted in their conserved position, were isolated from cotton genome. Quantitative reverse transcriptase-polymerase chain reaction (RT-PCR) analysis indicated that the nine GhBCP genes are differentially expressed in cotton tissues. Among them, GhBCP1 and GhBCP4 were predominantly expressed in fibers, while the transcripts of GhBCP2 and GhBCP3 were accumulated at relatively high levels in fibers. These four genes were strongly expressed in early fiber elongation, but dramatically declined with further fiber development. In addition, these GhBCP genes were upregulated in fibers by Cu(2+) , Zn(2+) , high-salinity and drought stresses, but downregulated in fibers by Al(3+) treatment. Overexpression of GhBCP1 and GhBCP4 in yeast (Schizosaccharomyces pombe) significantly increased the cell growth rate under Cu(2+) , Zn(2+) and high-salinity stresses. These results suggested that these GhBCPs may participate in the regulation of fiber development and in response to high-salinity and heavy metal stresses in cotton.

The 1.4 A resolution structure of Paracoccus pantotrophus pseudoazurin

Pseudoazurins are small type 1 copper proteins that are involved in the flow of electrons between various electron donors and acceptors in the bacterial periplasm, mostly under denitrifying conditions. The previously determined structure of Paracoccus pantotrophus pseudoazurin in the oxidized form was improved to a nominal resolution of 1.4 A, with R and R(free) values of 0.188 and 0.206, respectively. This high-resolution structure makes it possible to analyze the interactions between the monomers and the solvent structure in detail. Analysis of the high-resolution structure revealed the structural regions that are responsible for monomer-monomer recognition during dimer formation and for protein-protein interaction and that are important for partner recognition. The pseudoazurin structure was compared with other structures of various type 1 copper proteins and these were grouped into families according to similarities in their secondary structure; this may be useful in the annotation of copper proteins in newly sequenced genomes and in the identification of novel copper proteins.

Homologues of nitrite reductases in ammonia-oxidizing archaea: diversity and genomic context

Ammonia-oxidizing archaea are frequent and ubiquitous inhabitants of terrestrial and marine environments. As they have only recently been detected, most aspects of their metabolism are yet unknown. Here we report on the occurrence of genes encoding potential homologues of copper-dependent nitrite reductases (NirK) in ammonia-oxidizing archaea of soils and other environments using metagenomic approaches and PCR amplification. Two pairs of highly overlapping 40 kb genome fragments, each containing nirK genes of archaea, were isolated from a metagenomic soil library. Between 68% and 85% of the open reading frames on these genome fragments had homologues in the genomes of the marine archaeal ammonia oxidizers Nitrosopumilus maritimus and Cenarchaeum symbiosum. Extensions of NirK homologues with C-terminal fused amicyanin domains were deduced from two of the four fosmids indicating structural variation of these multicopper proteins in archaea. Phylogenetic analyses including all major groups of currently known NirK homologues revealed that the deduced protein sequences of marine and soil archaea were separated into two highly divergent lineages that did not contain bacterial homologues. In contrast, another separated lineage contained potential multicopper oxidases of both domains, archaea and bacteria. More nirK gene variants directly amplified by PCR from several environments indicated further diversity of the gene and a widespread occurrence in archaea. Transcription of the potential archaeal nirK in soil was demonstrated at different water contents, but no significant increase in transcript copy number was observed with increased denitrifying activity.

The genome sequence of the metal-mobilizing, extremely thermoacidophilic archaeon Metallosphaera sedula provides insights into bioleaching-associated metabolism

Despite their taxonomic description, not all members of the order Sulfolobales are capable of oxidizing reduced sulfur species, which, in addition to iron oxidation, is a desirable trait of biomining microorganisms. However, the complete genome sequence of the extremely thermoacidophilic archaeon Metallosphaera sedula DSM 5348 (2.2 Mb, approximately 2,300 open reading frames [ORFs]) provides insights into biologically catalyzed metal sulfide oxidation. Comparative genomics was used to identify pathways and proteins involved (directly or indirectly) with bioleaching. As expected, the M. sedula genome contains genes related to autotrophic carbon fixation, metal tolerance, and adhesion. Also, terminal oxidase cluster organization indicates the presence of hybrid quinol-cytochrome oxidase complexes. Comparisons with the mesophilic biomining bacterium Acidithiobacillus ferrooxidans ATCC 23270 indicate that the M. sedula genome encodes at least one putative rusticyanin, involved in iron oxidation, and a putative tetrathionate hydrolase, implicated in sulfur oxidation. The fox gene cluster, involved in iron oxidation in the thermoacidophilic archaeon Sulfolobus metallicus, was also identified. These iron- and sulfur-oxidizing components are missing from genomes of nonleaching members of the Sulfolobales, such as Sulfolobus solfataricus P2 and Sulfolobus acidocaldarius DSM 639. Whole-genome transcriptional response analysis showed that 88 ORFs were up-regulated twofold or more in M. sedula upon addition of ferrous sulfate to yeast extract-based medium; these included genes for components of terminal oxidase clusters predicted to be involved with iron oxidation, as well as genes predicted to be involved with sulfur metabolism. Many hypothetical proteins were also differentially transcribed, indicating that aspects of the iron and sulfur metabolism of M. sedula remain to be identified and characterized.

A Solid State 13C NMR, Crystallographic, and Quantum Chemical Investigation of Chemical Shifts and Hydrogen Bonding in Histidine Dipeptides

We report the first solid-state NMR, crystallographic, and quantum chemical investigation of the origins of the 13C NMR chemical shifts of the imidazole group in histidine-containing dipeptides. The chemical shift ranges for Cgamma and Cdelta2 seen in eight crystalline dipeptides were very large (12.7-13.8 ppm); the shifts were highly correlated (R2= 0.90) and were dominated by ring tautomer effects and intermolecular interactions. A similar correlation was found in proteins, but only for buried residues. The imidazole 13C NMR chemical shifts were predicted with an overall rms error of 1.6-1.9 ppm over a 26 ppm range, by using quantum chemical methods. Incorporation of hydrogen bond partner molecules was found to be essential in order to reproduce the chemical shifts seen experimentally. Using AIM (atoms in molecules) theory we found that essentially all interactions were of a closed shell nature and the hydrogen bond critical point properties were highly correlated with the N...H...O (average R2= 0.93) and Nepsilon2...H...N (average R2= 0.98) hydrogen bond lengths. For Cepsilon1, the 13C chemical shifts were also highly correlated with each of these properties (at the Nepsilon2 site), indicating the dominance of intermolecular interactions for Cepsilon1. These results open up the way to analyzing 13C NMR chemical shifts, tautomer states (from Cdelta2, Cepsilon1 shifts), and hydrogen bond properties (from Cepsilon1 shifts) of histidine residue in proteins and should be applicable to imidazole-containing drug molecules bound to proteins, as well.