Phylogenetic analysis, structural evolution and functional divergence of the 12-oxo-phytodienoate acid reductase gene family in plants

Background  

The 12-oxo-phytodienoic acid reductases (OPRs) are enzymes that catalyze the reduction of double-bonds in α, β -unsaturated aldehydes or ketones and are part of the octadecanoid pathway that converts linolenic acid to jasmonic acid. In plants, OPRs belong to the old yellow enzyme family and form multigene families. Although discoveries about this family in Arabidopsis and other species have been reported in some studies, the evolution and function of multiple OPRs in plants are not clearly understood.     

Phylogenetic analysis, structural evolution and functional divergence of the 12-oxo-phytodienoate acid reductase gene family in plants

Background

Whenever different data sets arrive at conflicting phylogenetic hypotheses, only testable causal explanations of sources of errors in at least one of the data sets allow us to critically choose among the conflicting hypotheses of relationships. The large (28S) and small (18S) subunit rRNAs are among the most popular markers for studies of deep phylogenies. However, some nodes supported by this data are suspected of being artifacts caused by peculiarities of the evolution of these molecules. Arthropod phylogeny is an especially controversial subject dotted with conflicting hypotheses which are dependent on data set and method of reconstruction. We assume that phylogenetic analyses based on these genes can be improved further i) by enlarging the taxon sample and ii) employing more realistic models of sequence evolution incorporating non-stationary substitution processes and iii) considering covariation and pairing of sites in rRNA-genes.

Results

We analyzed a large set of arthropod sequences, applied new tools for quality control of data prior to tree reconstruction, and increased the biological realism of substitution models. Although the split-decomposition network indicated a high noise content in the data set, our measures were able to both improve the analyses and give causal explanations for some incongruities mentioned from analyses of rRNA sequences. However, misleading effects did not completely disappear.

Conclusion

Analyses of data sets that result in ambiguous phylogenetic hypotheses demand for methods, which do not only filter stochastic noise, but likewise allow to differentiate phylogenetic signal from systematic biases. Such methods can only rely on our findings regarding the evolution of the analyzed data. Analyses on independent data sets then are crucial to test the plausibility of the results. Our approach can easily be extended to genomic data, as well, whereby layers of quality assessment are set up applicable to phylogenetic reconstructions in general.     

Structural and functional analysis reveals that human OASL binds dsRNA to enhance RIG-I signaling

The oligoadenylate synthetase (OAS) enzymes are cytoplasmic dsRNA sensors belonging to the antiviral innate immune system. Upon binding to viral dsRNA, the OAS enzymes synthesize 2′-5′ linked oligoadenylates (2-5As) that initiate an RNA decay pathway to impair viral replication. The human OAS-like (OASL) protein, however, does not harbor the catalytic activity required for synthesizing 2-5As and differs from the other human OAS family members by having two C-terminal ubiquitin-like domains. In spite of its lack of enzymatic activity, human OASL possesses antiviral activity. It was recently demonstrated that the ubiquitin-like domains of OASL could substitute for K63-linked poly-ubiquitin and interact with the CARDs of RIG-I and thereby enhance RIG-I signaling. However, the role of the OAS-like domain of OASL remains unclear. Here we present the crystal structure of the OAS-like domain, which shows a striking similarity with activated OAS1. Furthermore, the structure of the OAS-like domain shows that OASL has a dsRNA binding groove. We demonstrate that the OAS-like domain can bind dsRNA and that mutating key residues in the dsRNA binding site is detrimental to the RIG-I signaling enhancement. Hence, binding to dsRNA is an important feature of OASL that is required for enhancing RIG-I signaling.     

Phylogenetic and transcriptional analysis of a strictosidine synthase-like gene family in Arabidopsis thaliana reveals involvement in plant defence responses

Protein domains with similarity to plant strictosidine synthase-like (SSL) sequences have been uncovered in the genomes of all multicellular organisms sequenced so far and are known to play a role in animal immune responses. Among several distinct groups of Arabidopsis thaliana SSL sequences, four genes ( AtSSL4–AtSSL7 ) arranged in tandem on chromosome 3 show more similarity to SSL genes from Drosophila melanogaster and Caenorhabditis elegans than to other Arabidopsis SSL genes. To examine whether any of the four AtSSL genes are immune-inducible, we analysed the expression of each of the four AtSSL genes after exposure to microbial pathogens, wounding and plant defence elicitors using real-time quantitative RT-PCR, Northern blot hybridisation and Western blot analysis with antibodies raised against recombinant At SSL proteins. While the AtSSL4 gene was constitutively expressed and not significantly induced by any treatment, the other three AtSSL genes were induced to various degrees by plant defence signalling compounds, such as salicylic acid, methyl jasmonate and ethylene, as well as by wounding and exposure to the plant pathogens Alternaria brassicicola and cucumber mosaic virus . Our data demonstrate that the four SSL-coding genes are regulated individually, suggesting specific roles in basal ( SSL4 ) and inducible ( SSL5-7 ) plant defence mechanisms.     

Identification of novel members reveals the structural and functional divergence of lepidopteran-specific Lipoprotein_11 family

30K proteins (30KPs) are classified into the lepidopteran-specific Lipoprotein_11 family. They are involved in various physiological processes such as energy storage, embryonic development, and immune response in the silkworm. To date, 30KPs were only found in Bombyx mori and Manduca sexta. Moreover, the C-termini of ENF peptide binding proteins (ENF-BPs) show similarity to 30KPs. ENF peptides are multifunctional insect cytokines and involved in growth regulation and defense reaction, whereas ENF-BPs act as active regulators of ENF peptides. In order to get insights into this gene family in Lepidoptera, we performed an extensive survey of lepidopteran-derived genome and EST datasets. We identified 73 30KP homologous genes in 12 lepidopteran species, of which 56 are novel members. The structural and phylogenetic analyses revealed that these genes could be classified into three groups: ENF-BP genes, typical 30KP genes, and serine/threonine-rich 30KP (S/T-rich 30KP) genes. The C-terminal regions are common to all the three subfamilies, but the N-termini are highly variable. We found a novel subfamily of Lipoprotein_11 and named it S/T-rich 30KP according to its exclusive S/T-rich domain in the N terminus. ENF-BP was also found to contain a special domain in the N terminus, which is homologous to Pp-0912 of Pseudomonas putida. Microarray data and semi-quantitative RT-PCR showed that the three groups have their respective temporal–spatial expression patterns. S/T-rich 30KP genes have enriched expression in the mature testis and might be involved in spermiogenesis or fertilization. Typical 30KP genes are expressed mainly in the fat body and integument at the larvae and pupae stages. ENF-BP genes are expressed predominantly in the hemocyte. The differential spatial–temporal expression profiles revealed the functional divergence of three Lipoprotein_11 subfamilies.     

Comparative transcriptional profiling and evolutionary analysis of the GRAM domain family in eukaryotes

The GRAM domain was found in glucosyltransferases, myotubularins and other membrane-associated proteins. So far, functions for majority of these proteins are yet to be uncovered. In order to address the evolutionary and functional significance of this family members, we have performed a comprehensive investigation on their genome-wide identification, phylogenetic relationship and expression divergence in five different organisms representing monocot/dicot plants, vertebrate/invertebrate animals and yeast, namely, Oryza sativa, Arabidopsis thaliana, Mus musculus, Drosophila melanogaster and Saccharomyces cerevisiae, respectively. We have identified 65 members of GRAM domain family from these organisms. Our data revealed that this family was an ancient group and various organisms had evolved into different family sizes. Large-scale genome duplication and divergence in both expression patterns and functions were significantly contributed to the expansion and retention of this family. Mouse and Drosophila members showed higher divergences in their proteins as indicated by higher Ka/Ks ratios and possessed multiple domains in various combinations. However, in plants, their protein functions were possibly retained with a relatively low divergence as signified by lower Ka/Ks ratios and only one additional domain was combined during evolution. On the other hand, this family in all five organisms exhibited high divergence in their expression patterns both at tissue level and under various biotic and abiotic stresses. These highly divergent expression patterns unraveled the complexity of functions of GRAM domain family. Each member may play specialized roles in a specific tissue or stress condition and may function as regulators of environmental and hormonal signaling.     

Phylogenetic footprinting reveals a nuclear protein which binds to silencer sequences in the human gamma and epsilon globin genes.

Tissue- and developmental stage-specific expression of the human beta-like globin genes is regulated by a combination of ubiquitous and erythroid-restricted trans factors that bind to cis elements near each of the five active genes. Additional interactions of these cis and trans factors with sequences located in the far 5' end of the cluster occur by as yet obscure mechanisms. Because of the complexity of this regulatory puzzle, precise identification of the determinants that control hemoglobin switching has proven difficult. Phylogenetic footprinting is an evolutionary approach to this problem which is based on the supposition that the basic mechanisms of switching are conserved throughout mammalian phylogeny. Alignment of the 5' flanking regions of the gamma genes of several species allows the identification of footprints of 100% conserved sequence. We have now tested oligomers spanning 13 such phylogenetic footprints and find that 12 are bound by nuclear proteins. One conserved element located at -1086 from the gamma genes exhibits repressor activity in transient transfection studies. The protein that binds this element, CSBP-1 (conserved sequence-binding protein 1), also binds at three sites within a silencer element upstream from the epsilon globin gene. Further analysis reveals that the CSBP-1 binding activity is identical to that of a recently cloned zinc finger protein that has been shown to act as a repressor in other systems. The binding of CSPB-1 to silencer sequences in the epsilon and gamma globin genes may be important in the stage-specific silencing of these genes.     

Molecular evolution and functional divergence of the bestrophin protein family

Background  

Mutations in human bestrophin 1 are associated with at least three autosomal-dominant macular dystrophies including Best disease, adult onset vitelliform macular dystrophy and autosomal dominant vitreo-retinochoroidopathy. The protein is integral to the membrane and is likely involved in Ca²⁺-dependent transport of chloride ions across cellular membranes. Bestrophin 1 together with its three homologues forms a phylogenetically highly conserved family of proteins.     

Structural and functional analysis of the human HDAC4 catalytic domain reveals a regulatory structural zinc-binding domain

Histone deacetylases (HDACs) regulate chromatin status and gene expression, and their inhibition is of significant therapeutic interest. To date, no biological substrate for class IIa HDACs has been identified, and only low activity on acetylated lysines has been demonstrated. Here, we describe inhibitor-bound and inhibitor-free structures of the histone deacetylase-4 catalytic domain (HDAC4cd) and of an HDAC4cd active site mutant with enhanced enzymatic activity toward acetylated lysines. The structures presented, coupled with activity data, provide the molecular basis for the intrinsically low enzymatic activity of class IIa HDACs toward acetylated lysines and reveal active site features that may guide the design of class-specific inhibitors. In addition, these structures reveal a conformationally flexible structural zinc-binding domain conserved in all class IIa enzymes. Importantly, either the mutation of residues coordinating the structural zinc ion or the binding of a class IIa selective inhibitor prevented the association of HDAC4 with the N-CoR.HDAC3 repressor complex. Together, these data suggest a key role of the structural zinc-binding domain in the regulation of class IIa HDAC functions.     

Phylogenetic analysis of the vertebrate galectin family

Galectins form a family of structurally related carbohydrate binding proteins (lectins) that have been identified in a large variety of metazoan phyla. They are involved in many biological processes such as morphogenesis, control of cell death, immunological response, and cancer. To elucidate the evolutionary history of galectins and galectin-like proteins in chordates, we have exploited three independent lines of evidence: (i) location of galectin encoding genes (LGALS) in the human genome; (ii) exon-intron organization of galectin encoding genes; and (iii) sequence comparison of carbohydrate recognition domains (CRDs) of chordate galectins. Our results suggest that a duplication of a mono-CRD galectin gene gave rise to an original bi-CRD galectin gene, before or early in chordate evolution. The N-terminal and C-terminal CRDs of this original galectin subsequently diverged into two different subtypes, defined by exon-intron structure (F4-CRD and F3-CRD). We show that all vertebrate mono-CRD galectins known to date belong to either the F3- or F4- subtype. A sequence of duplication and divergence events of the different galectins in chordates is proposed.     

Identification of a domain that mediates vesicle aggregation reveals functional diversity of annexin repeats.

Annexins are structurally related proteins that bind phospholipids in a Ca2(+)-dependent manner and possess at least four conserved 70-amino acid repeat domains. The ability of certain annexins to promote contact between vesicle membranes in vitro has prompted the suggestion that these proteins regulate membrane traffic in exocytosis. We have previously found that annexins I and II promote contact between vesicles whereas annexin V does not. In order to understand the mechanism of annexin I-mediated vesicle-vesicle contact, we prepared a monoclonal antibody that specifically inhibits annexin I-mediated vesicle aggregation. We identified the domain of annexin I recognized by this monoclonal antibody by using it to screen an expression library containing random fragments of annexin I cDNA. The antibody identified a fragment encoding amino acids 41-118 (the first repeat plus 8 residues of the amino-terminal tail). We constructed a chimeric protein containing these amino acids of annexin I fused to the second, third, and fourth repeats of annexin V. Transfer of this domain conferred the ability to promote vesicle aggregation, confirming that this domain participates directly in mediating contact between vesicle membranes.