Exon structure conservation despite low sequence similarity: a relic of dramatic events in evolution?

In Silico identification and characterization of SOS gene family in soybean: Potential of calcium in salinity stress mitigation

Leguminous crops are usually sensitive to saline stress during germination and plant growth stages. The Salt Overly Sensitive (SOS) pathway is one of the key signaling pathways involved in salt translocation and tolerance in plants however, it is obscure in soybean. The current study describes the potential of calcium application on the mitigation of salinity stress and its impact on seed germination, morphological, physiological and biochemical attributes of soybean. The seeds from previously reported salt-tolerant and salt-susceptible soybean varieties were primed with water, calcium (10 and 20 mM), and stressed under 60, 80 and 100 mM NaCl and evaluated in various combinations. Results show that germination increased by 7% in calcium primed non-stressed seeds under non-stressing, whereas an improvement of 15%-25% was observed in germination under NaCl stress. Likewise, improvement in seedling length (3%-8%), plant height (9%-18%), number of nodes (3%-14%), SOD activity (20%) and Na+/K+ concentration (3%-5% reduction) in calcium primed plants, indicates alleviation of salinity-induced negative effects. In addition, this study also included in silico identification and confirmation of presence of Arabidopsis thaliana SOS genes orthologs in soybean. The research of amino acid sequences of SOS proteins from Arabidopsis thaliana (AtSOSs) within Glycine max genome displayed protein identity (60-80%) thus these identified homologs were called as GmSOS. Further phylogeny and in silico analyses showed that GmSOS orthologs contain similar gene structures, close evolutionary relationship, and same conserved motifs, reinforcing that GmSOSs belong to SOS family and they share many common features with orthologs from other species thus may perform similar functions. This is the first study that reports role of SOSs in salt-stress mitigation in soybean.

Unlocking the role of novel primary/di-amine oxidases in crop improvement: Tissue specificity leads to specific roles connected to abiotic stress, hormone responses and sensing nitrogen

Genetic improvements of solanaceous crops for quality and stress responsive traits are needed because of the central role vegetables and fruits have in providing nutrients to human diets. Copper amine oxidase (CuAO) encoding genes involved in metabolism of primary/di-amine nitrogenous compounds, play a role in balancing internal nitrogen (N) pools especially when external N supply fluctuates during growth, development and environmental stresses. In the present study, we investigated the occurrence, molecular evolution and possible role(s) of these unknown genes in tomato crops. Multiple genome-wide bioinformatics approaches led to the identification of eight bona fide CuAO genes (SlCuAO1-SlCuAO8) in the tomato genome with gene numbers like those in Arabidopsis and rice indicating their conserved functional relevance with a tandemly duplicated SlCuAO6-SlCuAO7 pair at chr.9. A conserved intron-exon size and phase distribution for SlCuAO2, 3, 4 pairs are similar to a recently identified single duckweed SpCuAO1 orthologue gene indicating its evolutionary conservation. Synteny analysis showed their closest association to Arabidopsis and but not with rice. Transcriptome data indicated that gene expression for about six genes (SlCuAO1, 2, 3, 4, 6, 7) is root specific, fruit specific for SlCuAO5 and flower specific for SlCuAO8 thus indicating amine oxidation is variable across tissues with a prominance in the root tissue. The majority of CuAO genes are negatively regulated by methyl jasmonate. Positive regulation, however, involves CuAO3/8. Transcript analysis of the ethylene-deficient transgenic lines indicated that ethylene is required for activation of SlCuAO4. CuAO4 and CuAO5 exhibited most significant tissues-independent gene expression responses across various nitrogen regimes. Drought, heat and N stress identified CuAO5 as an overlapping highly expressed gene that corroborates with putrescine accumulation for free and conjugated forms with an opposite abundance of bound forms. Taken together our study highlights new insights into the roles of copper amine oxidation genes and identifies CuAO5 as a multiple stress induced gene that can be used in genetic improvement programs for combining heat, drought and nitrogen use efficiency related traits.

Genome-Wide Analysis of C/S1-bZIP Subfamilies in Populus tomentosa and Unraveling the Role of PtobZIP55/21 in Response to Low Energy

C/S1 basic leucine zipper (bZIP) transcription factors are essential for plant survival under energy deficiency. However, studies on the responses of C/S1-bZIPs to low energy in woody plants have not yet been reported. In this study, members of C/S1-bZIP subfamilies in Populus tomentosa were systematically analyzed using bioinformatic approaches. Four C-bZIPs and 10 S1-bZIPs were identified, and their protein properties, phylogenetic relationships, gene structures, conserved motifs, and uORFs were systematically investigated. In yeast two-hybrid assays, direct physical interactions between C-bZIP and S1-bZIP members were observed, highlighting their potential functional synergy. Moreover, expression profile analyses revealed that low energy induced transcription levels of most C/S1-bZIP members, with bZIP55 and bZIP21 (a homolog of bZIP55) exhibiting particularly significant upregulation. When the expression of bZIP55 and bZIP21 was co-suppressed using artificial microRNA mediated gene silencing in transgenic poplars, root growth was promoted. Further analyses revealed that bZIP55/21 negatively regulated the root development of P. tomentosa in response to low energy. These findings provide insights into the molecular mechanisms by which C/S1-bZIPs regulate poplar growth and development in response to energy deprivation.

Genomic survey and expression analysis of cellulose synthase superfamily and COBRA-like gene family in Zanthoxylum bungeanum stipule thorns

The Cellulose Synthase gene (CS) superfamily and COBRA-like (COBL) gene family are essential for synthesizing cellulose and hemicellulose, which play a crucial role in cell wall biosynthesis and the hardening of plant tissues. Our study identified 126 ZbCS and 31 ZbCOBL genes from the Zanthoxylum bungeanum (Zb) genome. Phylogenetic analysis and conservative domain analysis unfolded that ZbCS and ZbCOBL genes were divided into seven and two subfamilies, respectively. Gene duplication data suggested that more than 75% of these genes had tandem and fragment duplications. Codon usage patterns analysis indicated that the ZbCS and ZbCOBL genes prefer ending with A/T base, with weak codon preference. Furthermore, seven key ZbCS and five key ZbCOBL genes were identified based on the content of cellulose and hemicellulose and the expression characteristics of ZbCS and ZbCOBL genes in various stages of stipule thorns. Altogether, these results improve the understanding of CS and COBL genes and provide valuable reference data for cultivating Zb with soft thorns.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01432-x.

Genome-wide characterization of LBD transcription factors in switchgrass (Panicum virgatum L.) and the involvement of PvLBD12 in salt tolerance

PvLBD12 enhanced the salt tolerance by increasing proline accumulation, improving K⁺ accumulation, and decreasing reactive oxygen species level in switchgrass. Abiotic stresses are the serious factors which limit plant development and productivity and restrict the agricultural economy. It is important, therefore, to understand the mechanism of abiotic tolerance in plants. Lateral organ boundaries domain (LBD) proteins as plant-specific transcription factors play important function in plant lateral organ development, plant regeneration, and abiotic stress. In our study, we identify 69 LBD members from switchgrass genome-wide sequences and classify them based on their homology with LBD proteins in Arabidopsis. RT-qPCR showed that PvLBD genes had different expression patterns under abiotic stress conditions, indicating that they play important roles in various stress. PvLBD12 was selected as a candidate gene for further functional analysis because it had the highest expression level under salt stress. Overexpression of PvLBD12 enhanced salt tolerance by altering a wide range of physiological responses (like increased proline accumulation, reduced malondialdehyde production, improved K⁺ accumulation, and reduced Na⁺ absorption) in switchgrass. Some stress response genes such as proline biosynthesis gene PvP5CS1, vacuolar Na⁺(K⁺)/H⁺ antiporter gene PvNHX1, two key ROS-scavenging enzyme genes PvCAT and PvSOD were all upregulated in PvLBD12 overexpression lines. Taken together, PvLBD12 plays a pivotal role in response to salt stress by increasing proline accumulation, improving K⁺ accumulation, reducing Na⁺ absorption, and decreasing reactive oxygen species level. It will be better to understand the potential biological functions of LBD genes in other plants.

Genome-wide identification glutathione-S-transferase gene superfamily in Daphnia pulex and its transcriptional response to nanoplastics

Glutathione S-transferases (GSTs) are key multifunctional phase II detoxification enzymes involved in the regulation of growth, development, and stress responses. However, the knowledge of GSTs in the model invertebrate organism Daphnia pulex at the genomic level remains limited. In the present study, 35 GST genes were identified in D. pulex (Dp-GST), belonging to eight subfamilies, with the sigma, mu, and delta/epsilon subfamilies constituting approximately 29 %, 20 %, and 20 % of the GST superfamily, respectively. Chromosome tandem duplication of genes within the same subfamily was observed, which may be the main force driving GST expansion in D. pulex. DpGST genes showed different expression patterns in response to nanoplastic exposure for 96 h and 21 days. Some homologous GST genes in D. pulex showed similar expression patterns in response to nanoplastic exposure, likely owing to their unique motifs. For example, motif 9 is found in all delta/epsilon GST genes, whereas motifs 1, 2, 3, 5, and 7 are highly conserved in sigma GST genes. The characterization of D. pulex GSTs extending the knowledge of GST-mediated environmental contaminants, especially nanoplastics.

Gene Structure Evolution of the Short-Chain Dehydrogenase/Reductase (SDR) Family

SDR (Short-chain Dehydrogenases/Reductases) are one of the oldest and heterogeneous superfamily of proteins, whose classification is problematic because of the low percent identity, even within families. To get clearer insights into SDR molecular evolution, we explored the splicing site organization of the 75 human SDR genes across their vertebrate and invertebrate orthologs. We found anomalous gene structures in members of the human SDR7C and SDR42E families that provide clues of retrogene properties and independent evolutionary trajectories from a common invertebrate ancestor. The same analyses revealed that the identity value between human and invertebrate non-allelic variants is not necessarily associated with the homologous gene structure. Accordingly, a revision of the SDR nomenclature is proposed by including the human SDR40C1 and SDR7C gene in the same family.

Evolution of γ chain cytokines: Mechanisms, methods and applications

The common γ chain family of cytokines and their receptors play fundamental roles in the immune system. Evolutionary studies of γ chain cytokines have elegantly illustrated how the immune system adapts to ever-changing environmental conditions. Indeed, these studies have revealed the uniqueness of cytokine evolution, which exhibits strong positive selection pressure needed to adapt to rapidly evolving threats whilst still conserving their receptor binding capabilities. In this review, we summarise the evolutionary mechanisms that gave rise to the characteristically diverse family of γ chain cytokines. We also speculate on the benefits of studying cytokine evolution, which may provide alternative ways to design novel cytokine therapeutic strategies. Additionally, we discuss current evolutionary models that elucidate the emergence of distinct cytokines (IL-4 and IL-13) and cytokine receptors (IL-2Rα and IL-15Rα). Finally, we address and reflect on the difficulties associated with evolutionary studies of rapidly evolving genes and describe a variety of computational methods that have revealed numerous aspects of cytokine evolution.

CCL25 Inhibition Alleviates Sepsis-Induced Acute Lung Injury and Inflammation

Purpose

Acute lung injury (ALI) is a common clinical syndrome with high mortality. The chemokine ligand 25 (CCL25) is involved in inflammation, leukocyte trafficking and immunoregulation. However, the role and mechanism of CCL25 in ALI are not fully understood yet. The aim of this study was to explore the relationship between acute lung injury and CCL25.

Patients and Methods

In this study, we first examined chemokine expression in sepsis patients and found that serum CCL25 expression levels were relatively high in sepsis patients compared to healthy individuals. Based on this, we designed in vitro and in vivo experiments to verify the validity of the theory. In vitro, we used lipopolysaccharide-stimulated human pulmonary microvascular endothelial cells (HPMECs). In vivo, we established male C57BL/6 mice cecal ligation puncture (CLP) model of sepsis.

Results

In vitro, we used lipopolysaccharide-stimulated human pulmonary microvascular endothelial cells (HPMECs) and found significantly higher expression of CCL25 by enzyme-linked immunosorbent assay. Inhibition of CCL25 resulted in a significant decrease in the expression of inflammatory cytokines in HPMECs. In addition, we found that CCL25 promoted increased endothelial permeability by reducing the expression of tight junction proteins and was associated with activation of the P38 MAPK pathway by measuring the transepithelial electrical resistance and fluorescence intensity of fluorescein isothiocyanate. Results from luciferase assays and chromatin immunoprecipitation assays showed that inhibition of NF-κB activity in HPMECs decreased CCL25 expression, but addition of recombinant CCL25 increased cell permeability and inflammatory cytokine expression. In vivo, we established male C57BL/6 mice cecal ligation puncture (CLP) model of sepsis. We found that inhibition of CCL25 significantly reduced inflammatory cytokine expression in a CLP-induced sepsis model, thereby alleviating lung tissue damage in mice.

Conclusion

Our study suggests that CCL25 contributed to the development of ALI by modulating the functions of microvascular endothelial cells.

Genome-wide identification and characterization of the Prunus persica ferredoxin gene family and its role in improving heat tolerance

Ferredoxin is involved in many biological processes, such as carbon fixation, nitrogen assimilation, chlorophyll metabolism, and fatty acid synthesis, and it plays a role in plant resistance to stress. However, the functions of Fds in peach during stress are unclear. In this study, 11 members of the peach Fd gene family were identified and divided into six groups (I- VI). We carried out bioinformatics analysis on these sequences, analyzed the physical and chemical properties of PpFd protein and the cis-elements in its promoter region, and predicted and compared the differences in gene structure and conserved protein motifs among groups. The results showed that the PpFd protein was highly conserved in plant species. In addition, overexpression of PpFd08 significantly increased the tolerance of transgenic tomato to high-temperature stress. The transcriptome analysis and qRT-PCR results of PpFd08 transgenic apple calli showed that PpFd08 might enhance heat resistance by modulating the expression of heat tolerance related genes. The results of this study provide a new understanding for the further study of the function of PpFd protein in peach and a candidate gene for improving the heat resistance of peach.

Genome-wide bioinformatics analysis of Cellulose Synthase gene family in common bean (Phaseolus vulgaris L.) and the expression in the pod development

Background

CesA and Csl gene families, which belong to the cellulose synthase gene superfamily, plays an important role in the biosynthesis of the plant cell wall. Although researchers have investigated this gene superfamily in several model plants, to date, no comprehensive analysis has been conducted in the common bean.

Results

In this study, we identified 39 putative cellulose synthase genes from the common bean genome sequence. Then, we performed a bioinformatics analysis of this gene family involving sequence alignment, phylogenetic analysis, gene structure, collinearity analysis and chromosome location. We found all members possess a cellulose_synt domain. Phylogenetic analysis revealed that these cellulose synthase genes may be classified into five subfamilies, and that members in the same subfamily share conserved exon-intron distribution and motif compositions. Abundant and distinct cis-acting elements in the 2 k basepairs upstream regulatory regions indicate that the cellulose synthase gene family may plays a vital role in the growth and development of common bean. Moreover, the 39 cellulose synthase genes are distributed on 10 of the 11 chromosomes. Additionally expression analysis shows that all CesA/Csl genes selected are constitutively expressed in the pod development.

Conclusions

This research reveals both the putative biochemical and physiological functions of cellulose synthase genes in common bean and implies the importance of studying non-model plants to understand the breadth and diversity of cellulose synthase genes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01026-0.

OUP accepted manuscript

Homology of highly divergent genes often cannot be determined from sequence similarity alone. For example, we recently identified in the aphid Hormaphis cornu a family of rapidly evolving bicycle genes, which encode novel proteins implicated as plant gall effectors, and sequence similarity search methods yielded few putative bicycle homologs in other species. Coding sequence-independent features of genes, such as intron-exon boundaries, often evolve more slowly than coding sequences, however, and can provide complementary evidence for homology. We found that a linear logistic regression classifier using only structural features of bicycle genes identified many putative bicycle homologs in other species. Independent evidence from sequence features and intron locations supported homology assignments. To test the potential roles of bicycle genes in other aphids, we sequenced the genome of a second gall-forming aphid, Tetraneura nigriabdominalis and found that many bicycle genes are strongly expressed in the salivary glands of the gall forming foundress. In addition, bicycle genes are strongly overexpressed in the salivary glands of a non-gall forming aphid, Acyrthosiphon pisum, and in the non-gall forming generations of H. cornu. These observations suggest that Bicycle proteins may be used by multiple aphid species to manipulate plants in diverse ways. Incorporation of gene structural features into sequence search algorithms may aid identification of deeply divergent homologs, especially of rapidly evolving genes involved in host-parasite interactions.

Evolutionary history of histamine receptors: Early vertebrate origin and expansion of the H3-H4 subtypes

Histamine receptors belonging to the superfamily of G protein-coupled receptors (GPCRs) mediate the diverse biological effects of biogenic histamine. They are classified into four phylogenetically distinct subtypes H₁-H₄, each with a different binding affinity for histamine and divergent downstream signaling pathways. Here we present the evolutionary history of the histamine receptors using a phylogenetic approach complemented with comparative genomics analyses of the sequences, gene structures, and synteny of gene neighborhoods. The data indicate the earliest emergence of histamine-mediated GPCR signaling by a H₂ in a prebilaterian ancestor. The analyses support a revised classification of the vertebrate H₃-H₄ receptor subtypes. We demonstrate the presence of the H₄ across vertebrates, contradicting the currently held notion that H₄ is restricted to mammals. These non-mammalian vertebrate H₄ orthologs have been mistaken for H₃. We also identify the presence of a new H₃ subtype (H₃B), distinct from the canonical H₃ (H₃A), and propose that the H₃A, H₃B, and H₄ likely emerged from a H₃ progenitor through the 1R/2R whole genome duplications in an ancestor of the vertebrates. It is apparent that the ability of the H₁, H₂, and H_3-4 to bind histamine was acquired convergently. We identified genomic signatures suggesting that the H₁ and H₃-H₄ shared a last common ancestor with the muscarinic receptor in a bilaterian predecessor whereas, the H₂ and the α-adrenoreceptor shared a progenitor in a prebilaterian ancestor. Furthermore, site-specific analysis of the vertebrate subtypes revealed potential residues that may account for the functional divergence between them.

Genome-wide analysis of the abiotic stress-related bZIP family in switchgrass

The large basic leucine zipper (bZIP) transcription factor family is conserved in plants. These proteins regulate growth, development, and stress response. Here, we conducted a genome-wide analysis to identify the bZIP genes associated with stress resistance in switchgrass (Panicum virgatum L.). We identified 178 PvbZIPs unevenly distributed on 18 switchgrass chromosomes. An evolutionary analysis segregated them into 10 subfamilies. Gene structure and conserved motif analyses indicated that the same subfamily members shared similar intron-exon modes and motif compositions. This finding corroborated the proposed PvbZIP family grouping. A promoter analysis showed that PvbZIP genes participate in various stress responses. Phylogenetic and synteny analyses characterized 111 switchgrass bZIPs as orthologs of 70 rice bZIPs. A protein interaction network analysis revealed that 22 proteins are involved in salt and drought tolerance. An expression atlas disclosed that the expression patterns of several PvbZIPs differ among various tissues and developmental stages. Online data demonstrated that 16 PvbZIPs were significantly downregulated and five were significantly upregulated in response to heat stress. Other PvbZIPs participated in responses to abiotic stress such as salt, drought, cold, and heat. Our genome-wide analysis and identification of the switchgrass bZIP family characterized multiple candidate PvbZIPs that regulate growth and stress response. This study lays theoretical and empirical foundations for future functional investigations into other transcription factors.

Genome-wide analysis of glycerol-3-phosphate O-acyltransferase gene family and functional characterization of two cutin group GPATs in Brassica napus

Genome-wide identification, spatio-temporal expression analysis and functional characterization of selected Brassica napus GPATs highlight their roles in cuticular wax biosynthesis and defense against fungal pathogens. Glycerol-3-phosphate 1-O-acyltransferase (GPAT) is a key enzyme in the biosynthesis of glycerolipids, a major component of cellular membranes and extracellular protective layers, such as cuticles in plants. Brassica napus is an economically important crop and cultivated worldwide mostly for its edible oil. The B. napus GPATs (BnGPATs) are insufficiently characterized. Here, we performed genome-wide analysis to identify putative GPATs in B. napus and its diploid progenitors B. rapa and B oleracea. The 32 B. napus BnGPATs are phylogenetically divided into three major groups, cutin, suberin, and diverse ancient groups. Analysis of transcriptomes of different tissues and seeds at different developmental stages revealed the spatial and temporal expression profiles of BnGPATs. The yield and oil quality of B. napus are adversely affected by the necrotrophic fungus, Sclerotinia sclerotiorum. We showed that several BnGPATs, including cutin-related BnGPAT19 and 21, were upregulated in the S. sclerotiorum resistant line. RNAi-mediated suppression of BnGPAT19 and 21 in B. napus resulted in thinner cuticle, leading to rapid water and chlorophyll loss in toluidine blue staining and leaf bleaching assays. In addition, the RNAi plants also developed severe necrotic lesions following fungal inoculation compared to the wild-type plants, indicating that BnGPAT19 and 21 are likely involved in cuticular wax biosynthesis that is critical for initial pathogen defense. Taken together, we provided a comprehensive account of GPATs B. napus and characterized BnGPAT19 and 21 for their potential roles in cuticular wax biosynthesis and defense against fungal pathogens.

Genome-wide characterization of the cellulose synthase gene superfamily in Solanum lycopersicum

The cellulose synthase gene superfamily, which includes the cellulose synthase (CesA) and cellulose synthase-like (Csl) gene families, plays a vital role in the biosynthesis of cellulose and hemicellulose in plants. However, these genes have not been extensively studied in tomato (Solanum lycopersicum), a model for Solanaceae plants and for fleshy fruit development. Here, we identified and systematically analyzed 38 CesA/Csl family members that contained cellulose synthase domain regions, and categorized their encoded proteins into 6 subfamilies (CesA, CslA, CslB, CslD, CslE, and CslG) based on phylogenetic analysis. Most CesA/Csl genes from tomato are closely related to those from Arabidopsis, but the families have distinct features regarding gene structure, chromosome distribution and localization, phylogeny, and deduced protein sequence, indicating that they arose via different evolutionary process. Furthermore, expression analysis of CesA/Csl genes in different tissues at various developmental stages showed that most CesAs were constitutively expressed with differential expression levels in various organs; three CslD genes were expressed specifically in flowers, and four CesA and five Csl putative genes were preferentially expressed in fruits. Our results provide insight into the general characteristics of the CesA/Csl genes in tomato, and lay the foundation for further functional studies of CesA/Csl genes in tomato and other Solanaceae species.

Characterization of Cellulose Synthase A (CESA) Gene Family in Eudicots

Cellulose synthase A (CESA) is a key enzyme involved in the complex process of plant cell wall biosynthesis, and it remains a productive subject for research. We employed systems biology approaches to explore structural diversity of eudicot CESAs by exon-intron organization, mode of duplication, synteny, and splice site analyses. Using a combined phylogenetics and comparative genomics approach coupled with co-expression networks we reconciled the evolution of cellulose synthase gene family in eudicots and found that the basic forms of CESA proteins are retained in angiosperms. Duplications have played an important role in expansion of CESA gene family members in eudicots. Co-expression networks showed that primary and secondary cell wall modules are duplicated in eudicots. We also identified 230 simple sequence repeat markers in 103 eudicot CESAs. The 13 identified conserved motifs in eudicots will provide a basis for gene identification and functional characterization in other plants. Furthermore, we characterized (in silico) eudicot CESAs against senescence and found that expression levels of CESAs decreased during leaf senescence.

OMGene: mutual improvement of gene models through optimisation of evolutionary conservation

Background

The accurate determination of the genomic coordinates for a given gene – its gene model – is of vital importance to the utility of its annotation, and the accuracy of bioinformatic analyses derived from it. Currently-available methods of computational gene prediction, while on the whole successful, frequently disagree on the model for a given predicted gene, with some or all of the variant gene models often failing to match the biologically observed structure. Many prediction methods can be bolstered by using experimental data such as RNA-seq. However, these resources are not always available, and rarely give a comprehensive portrait of an organism’s transcriptome due to temporal and tissue-specific expression profiles.

Results

Orthology between genes provides evolutionary evidence to guide the construction of gene models. OMGene (Optimise My Gene) aims to improve gene model accuracy in the absence of experimental data by optimising the consistency of multiple sequence alignments of orthologous genes from multiple species. Using RNA-seq data sets from plants, mammals, and fungi, considering intron/exon junction representation and exon coverage, and assessing the intra-orthogroup consistency of subcellular localisation predictions, we demonstrate the utility of OMGene for improving gene models in annotated genomes.

Conclusions

We show that significant improvements in the accuracy of gene model annotations can be made, both in established and in de novo annotated genomes, by leveraging information from multiple species.

Systems Identification and Characterization of Cell Wall Reassembly and Degradation Related Genes in Glycine max (L.) Merill, a Bioenergy Legume

Soybean is a promising biomass resource for generation of second-generation biofuels. Despite the utility of soybean cellulosic biomass and post-processing residues in biofuel generation, there is no comprehensive information available on cell wall loosening and degradation related gene families. In order to achieve enhanced lignocellulosic biomass with softened cell walls and reduced recalcitrance, it is important to identify genes involved in cell wall polymer loosening and degrading. Comprehensive genome-wide analysis of gene families involved in cell wall modifications is an efficient stratagem to find new candidate genes for soybean breeding for expanding biofuel industry. We report the identification of 505 genes distributed among 12 gene families related to cell wall loosening and degradation. 1262 tandem duplication events contributed towards expansion and diversification of studied gene families. We identified 687 Simple Sequence Repeat markers and 5 miRNA families distributed on 316 and 10 genes, respectively. Publically available microarray datasets were used to explore expression potential of identified genes in soybean plant developmental stages, 68 anatomical parts, abiotic and biotic stresses. Co-expression networks revealed transcriptional coordination of different gene families involved in cell wall loosening and degradation process.

Genome and transcriptome-wide analyses of cellulose synthase gene superfamily in soybean

The plant cellulose synthase gene superfamily belongs to the category of type-2 glycosyltransferases, and is involved in cellulose and hemicellulose biosynthesis. These enzymes are vital for maintaining cell-wall structural integrity throughout plant life. Here, we identified 78 putative cellulose synthases (CS) in the soybean genome. Phylogenetic analysis against 40 reference Arabidopsis CS genes clustered soybean CSs into seven major groups (CESA, CSL A, B, C, D, E and G), located on 19 chromosomes (except chromosome 18). Soybean CS expansion occurred in 66 duplication events. Additionally, we identified 95 simple sequence repeat makers related to 44 CSs. We next performed digital expression analysis using publically available datasets to understand potential CS functions in soybean. We found that CSs were highly expressed during soybean seed development, a pattern confirmed with an Affymatrix soybean IVT array and validated with RNA-seq profiles. Within CS groups, CESAs had higher relative expression than CSLs. Soybean CS models were designed based on maximum average RPKM values. Gene co-expression networks were developed to explore which CSs could work together in soybean. Finally, RT-PCR analysis confirmed the expression of 15 selected CSs during all four seed developmental stages.

Genome-Wide Identification of bZIP Family Genes Involved in Drought and Heat Stresses in Strawberry (Fragaria vesca)

Basic leucine zipper (bZIP) genes are known to play a crucial role in response to various processes in plant as well as abiotic or biotic stress challenges. We have performed an identification and characterization of 50 bZIP genes across the woodland strawberry (Fragaria vesca) genome, which were divided into 10 clades according to the phylogenetic relationship of the strawberry bZIP proteins with those in Arabidopsis and rice. Five categories of intron patterns were observed within basic and hinge regions of the bZIP domains. Some additional conserved motifs have been found with the group specificity. Further, we predicted DNA-binding specificity of the basic and hinge regions as well as dimerization properties of leucine zipper regions, which was consistent with our phylogenetic clade and classified into 20 subfamilies. Across the different developmental stages of 15 organs and two types of fruits, the clade A bZIP members showed different tissue-specific expression patterns and the duplicated genes were differentially regulated, indicating a functional diversification coupled with the expansion of this gene family in strawberry. Under normal growth conditions, mrna11837 and mrna30280 of clade A showed very weak expression levels in organs and fruits, respectively; but higher expression was observed with different set of genes following drought and heat treatment, which may be caused by the separate response pathway between drought and heat treatments.

New fungal defensin-like peptides provide evidence for fold change of proteins in evolution

Defensins containing a consensus cystine framework, Cys_[1]…Cys_[2]X₃Cys_[3]…Cys_[4]… Cys_[5]X₁Cys_[6] (X, any amino acid except Cys; …, variable residue numbers), are extensively distributed in a variety of multicellular organisms (plants, fungi and invertebrates) and essentially involved in immunity as microbicidal agents. This framework is a prerequisite for forming the cysteine-stabilized α-helix and β-sheet (CSαβ) fold, in which the two invariant motifs, Cys_[2]X₃Cys_[3]/Cys_[5]X₁Cys_[6], are key determinants of fold formation. By using a computational genomics approach, we identified a large superfamily of fungal defensin-like peptides (fDLPs) in the phytopathogenic fungal genus – Zymoseptoria, which includes 132 structurally typical and 63 atypical members. These atypical fDLPs exhibit an altered cystine framework and accompanying fold change associated with their secondary structure elements and disulfide bridge patterns, as identified by protein structure modelling. Despite this, they definitely are homologous with the typical fDLPs in view of their precise gene structure conservation and identical precursor organization. Sequence and structural analyses combined with functional data suggest that most of Zymoseptoria fDLPs might have lost their antimicrobial activity. The present study provides a clear example of fold change in the evolution of proteins and is valuable in establishing remote homology among peptide superfamily members with different folds.

Annotating Protein Functional Residues by Coupling High-Throughput Fitness Profile and Homologous-Structure Analysis

Identification and annotation of functional residues are fundamental questions in protein sequence analysis. Sequence and structure conservation provides valuable information to tackle these questions. It is, however, limited by the incomplete sampling of sequence space in natural evolution. Moreover, proteins often have multiple functions, with overlapping sequences that present challenges to accurate annotation of the exact functions of individual residues by conservation-based methods. Using the influenza A virus PB1 protein as an example, we developed a method to systematically identify and annotate functional residues. We used saturation mutagenesis and high-throughput sequencing to measure the replication capacity of single nucleotide mutations across the entire PB1 protein. After predicting protein stability upon mutations, we identified functional PB1 residues that are essential for viral replication. To further annotate the functional residues important to the canonical or noncanonical functions of viral RNA-dependent RNA polymerase (vRdRp), we performed a homologous-structure analysis with 16 different vRdRp structures. We achieved high sensitivity in annotating the known canonical polymerase functional residues. Moreover, we identified a cluster of noncanonical functional residues located in the loop region of the PB1 β-ribbon. We further demonstrated that these residues were important for PB1 protein nuclear import through the interaction with Ran-binding protein 5. In summary, we developed a systematic and sensitive method to identify and annotate functional residues that are not restrained by sequence conservation. Importantly, this method is generally applicable to other proteins about which homologous-structure information is available.

IMPORTANCE

To fully comprehend the diverse functions of a protein, it is essential to understand the functionality of individual residues. Current methods are highly dependent on evolutionary sequence conservation, which is usually limited by sampling size. Sequence conservation-based methods are further confounded by structural constraints and multifunctionality of proteins. Here we present a method that can systematically identify and annotate functional residues of a given protein. We used a high-throughput functional profiling platform to identify essential residues. Coupling it with homologous-structure comparison, we were able to annotate multiple functions of proteins. We demonstrated the method with the PB1 protein of influenza A virus and identified novel functional residues in addition to its canonical function as an RNA-dependent RNA polymerase. Not limited to virology, this method is generally applicable to other proteins that can be functionally selected and about which homologous-structure information is available.

Analysis of basic leucine zipper genes and their expression during bud dormancy in peach (Prunus persica)

Dormancy is a biological characteristic developed to resist the cold conditions in winter. The bZIP transcription factors are present exclusively in eukaryotes and have been identified and classified in many species. bZIP proteins are known to regulate numerous biological processes, however, the role of bZIP in bud dodormancy has not been studied extensively. In total, 50 PpbZIP transcription factor-encoding genes were identified and categorized them into 10 groups (A-I and S). Similar intron/exon structures, additional conserved motifs, and DNA-binding site specificity supported our classification scheme. Additionally, chromosomal distribution and collinearity analyses suggested that expansion of the PpbZIP transcription factor family was due to segment/chromosomal duplications. We also predicted the dimerization properties based on characteristic features of the leucine zipper and classified PpbZIP proteins into 23 subfamilies. Furthermore, qRT-PCR results indicated that PpbZIPs genes may be involved in regulating dormancy. The same gene of different species might participate in different regulating networks through interactions with specific partners. Our expression profiling results complemented the microarray data, suggesting that co-expression patterns of bZIP transcription factors during dormancy differed among deciduous fruit trees. Our findings further clarify the molecular characteristics of the PpbZIP transcription factor family, including potential gene functions during dormancy. This information may facilitate further research on the evolutionary history and biological functions of bZIP proteins in peach and other rosaceae plants.

CAB-Align: A Flexible Protein Structure Alignment Method Based on the Residue-Residue Contact Area

Proteins are flexible, and this flexibility has an essential functional role. Flexibility can be observed in loop regions, rearrangements between secondary structure elements, and conformational changes between entire domains. However, most protein structure alignment methods treat protein structures as rigid bodies. Thus, these methods fail to identify the equivalences of residue pairs in regions with flexibility. In this study, we considered that the evolutionary relationship between proteins corresponds directly to the residue–residue physical contacts rather than the three-dimensional (3D) coordinates of proteins. Thus, we developed a new protein structure alignment method, contact area-based alignment (CAB-align), which uses the residue–residue contact area to identify regions of similarity. The main purpose of CAB-align is to identify homologous relationships at the residue level between related protein structures. The CAB-align procedure comprises two main steps: First, a rigid-body alignment method based on local and global 3D structure superposition is employed to generate a sufficient number of initial alignments. Then, iterative dynamic programming is executed to find the optimal alignment. We evaluated the performance and advantages of CAB-align based on four main points: (1) agreement with the gold standard alignment, (2) alignment quality based on an evolutionary relationship without 3D coordinate superposition, (3) consistency of the multiple alignments, and (4) classification agreement with the gold standard classification. Comparisons of CAB-align with other state-of-the-art protein structure alignment methods (TM-align, FATCAT, and DaliLite) using our benchmark dataset showed that CAB-align performed robustly in obtaining high-quality alignments and generating consistent multiple alignments with high coverage and accuracy rates, and it performed extremely well when discriminating between homologous and nonhomologous pairs of proteins in both single and multi-domain comparisons. The CAB-align software is freely available to academic users as stand-alone software at http://www.pharm.kitasato-u.ac.jp/bmd/bmd/Publications.html.

Genome-wide systematic characterization of the bZIP transcriptional factor family in tomato (Solanum lycopersicum L.)

BACKGROUND

Transcription factors of the basic leucine zipper (bZIP) family represent exclusively in eukaryotes and have been shown to regulate diverse biological processes in plant growth and development as well as in abiotic and biotic stress responses. However, little is known about the bZIP family in tomato (Solanum lycopersicum L.).

METHODS

The SlbZIP genes were identified using local BLAST and hidden Markov model profile searches. The phylogenetic trees, conserved motifs and gene structures were generated by MEGA6.06, MEME tool and gene Structure Display Server, respectively. The syntenic block diagrams were generated by the Circos software. The transcriptional gene expression profiles were obtained using Genevestigator tool and quantitative RT-PCR.

RESULTS

In the present study, we carried out a genome-wide identification and systematic analyses of 69 SlbZIP genes that distributes unevenly on the tomato chromosomes. This family can be divided into 9 groups according to the phylogenetic relationship among the SlbZIP proteins. Six kinds of intron patterns (a-f) within the basic and hinge regions are defined. The additional conserved motifs and their presence of the group specificity were also identified. Further, we predicted the DNA-binding patterns and the dimerization property on the basis of the characteristic features in the basic and hinge regions and the leucine zipper, respectively, which supports our classification greatly and helps to classify 24 distinct subfamilies. Within the SlbZIP family, a total of 40 SlbZIP genes are located in the segmental duplicate regions in the tomato genome, suggesting that the segment chromosomal duplications contribute greatly to the expansion of the tomato SlbZIP family. Expression profiling analyses of 59 SlbZIP genes using quantitative RT-PCR and publicly available microarray data indicate that the tomato SlbZIP genes have distinct and diverse expression patterns in different tissues and developmental stages and many of the tomato bZIP genes might be involved in responses to various abiotic and biotic stresses as well as in response to light.

CONCLUSIONS

This genome-wide systematic characterization identified a total of 69 members in the SlbZIP family and the analyses of the protein features and gene expression patterns provide useful clues for further functional characterization of the bZIP transcription factors in tomato.

Genome-wide systematic characterization of the bZIP transcriptional factor family in tomato (Solanum lycopersicum L.)

BACKGROUND

Transcription factors of the basic leucine zipper (bZIP) family represent exclusively in eukaryotes and have been shown to regulate diverse biological processes in plant growth and development as well as in abiotic and biotic stress responses. However, little is known about the bZIP family in tomato (Solanum lycopersicum L.).

METHODS

The SlbZIP genes were identified using local BLAST and hidden Markov model profile searches. The phylogenetic trees, conserved motifs and gene structures were generated by MEGA6.06, MEME tool and gene Structure Display Server, respectively. The syntenic block diagrams were generated by the Circos software. The transcriptional gene expression profiles were obtained using Genevestigator tool and quantitative RT-PCR.

RESULTS

In the present study, we carried out a genome-wide identification and systematic analyses of 69 SlbZIP genes that distributes unevenly on the tomato chromosomes. This family can be divided into 9 groups according to the phylogenetic relationship among the SlbZIP proteins. Six kinds of intron patterns (a-f) within the basic and hinge regions are defined. The additional conserved motifs and their presence of the group specificity were also identified. Further, we predicted the DNA-binding patterns and the dimerization property on the basis of the characteristic features in the basic and hinge regions and the leucine zipper, respectively, which supports our classification greatly and helps to classify 24 distinct subfamilies. Within the SlbZIP family, a total of 40 SlbZIP genes are located in the segmental duplicate regions in the tomato genome, suggesting that the segment chromosomal duplications contribute greatly to the expansion of the tomato SlbZIP family. Expression profiling analyses of 59 SlbZIP genes using quantitative RT-PCR and publicly available microarray data indicate that the tomato SlbZIP genes have distinct and diverse expression patterns in different tissues and developmental stages and many of the tomato bZIP genes might be involved in responses to various abiotic and biotic stresses as well as in response to light.

CONCLUSIONS

This genome-wide systematic characterization identified a total of 69 members in the SlbZIP family and the analyses of the protein features and gene expression patterns provide useful clues for further functional characterization of the bZIP transcription factors in tomato.

Evolutionary history of the neuropeptide S receptor/neuropeptide S system

The neuropeptide S receptor (NPSR) belongs to the G protein-coupled receptor (GPCR) superfamily and is activated by the neuropeptide S (NPS). Although recently discovered, the vertebrate NPSR-NPS system has been established as an important signaling system in the central nervous system and is involved in physiological processes such as locomotor activity, wakefulness, asthma pathogenesis, anxiety and food intake. The availability of a large number of genome sequences from multiple bilaterian lineages has provided an opportunity to establish the evolutionary history of the system. This review describes the origin and the molecular evolution of the NPSR-NPS system using data derived primarily from comparative genomic analyses. These analyses indicate that the NPSR-NPS system and the vasopressin-like receptor-vasopressin/oxytocin peptide (VPR-VP/OT) system originated from a single system in an ancestral bilaterian. Multiple duplications of this ancestral system gave rise to the bilaterian VPR-VP/OT system and to the protostomian cardioacceleratory peptide receptor-cardioacceleratory peptide (CCAPR-CCAP) system and to the NPSR-NPS system in the deuterostomes. Gene structure features of the receptors were consistent with the orthology annotations derived from phylogenetic analyses. The orthology of the peptide precursors closely paralleled that of the receptors suggesting an ancient coevolution of the receptor-peptide pair. An important challenge for the coevolution hypothesis will be to establish the molecular and structural basis of the divergence between orthologous receptor-ligand pairs in this system.

Pierced Lasso Bundles are a new class of knot-like motifs

A four-helix bundle is a well-characterized motif often used as a target for designed pharmaceutical therapeutics and nutritional supplements. Recently, we discovered a new structural complexity within this motif created by a disulphide bridge in the long-chain helical bundle cytokine leptin. When oxidized, leptin contains a disulphide bridge creating a covalent-loop through which part of the polypeptide chain is threaded (as seen in knotted proteins). We explored whether other proteins contain a similar intriguing knot-like structure as in leptin and discovered 11 structurally homologous proteins in the PDB. We call this new helical family class the Pierced Lasso Bundle (PLB) and the knot-like threaded structural motif a Pierced Lasso (PL). In the current study, we use structure-based simulation to investigate the threading/folding mechanisms for all the PLBs along with three unthreaded homologs as the covalent loop (or lasso) in leptin is important in folding dynamics and activity. We find that the presence of a small covalent loop leads to a mechanism where structural elements slipknot to thread through the covalent loop. Larger loops use a piercing mechanism where the free terminal plugs through the covalent loop. Remarkably, the position of the loop as well as its size influences the native state dynamics, which can impact receptor binding and biological activity. This previously unrecognized complexity of knot-like proteins within the helical bundle family comprises a completely new class within the knot family, and the hidden complexity we unraveled in the PLBs is expected to be found in other protein structures outside the four-helix bundles. The insights gained here provide critical new elements for future investigation of this emerging class of proteins, where function and the energetic landscape can be controlled by hidden topology, and should be take into account in ab initio predictions of newly identified protein targets.

Genomic surveys and expression analysis of bZIP gene family in castor bean (Ricinus communis L.)

The basic leucine zipper (bZIP) transcription factors comprise a family of transcriptional regulators present extensively in plants, involved in regulating diverse biological processes such as flower and vascular development, seed maturation, stress signaling and pathogen defense. Castor bean (Ricinus communis L. Euphorbiaceae) is one of the most important non-edible oilseed crops and its seed oil is broadly used for industrial applications. We performed a comprehensive genome-wide identification and analysis of the bZIP transcription factors that exist in the castor bean genome in this study. In total, 49 RcbZIP transcription factors were identified, characterized and categorized into 11 groups (I-XI) based on their gene structure, DNA-binding sites, conserved motifs, and phylogenetic relationships. The dimerization properties of 49 RcbZIP proteins were predicted on the basis of the characteristic features in the leucine zipper. Global expression profiles of 49 RcbZIP genes among different tissues were examined using high-throughput sequencing of digital gene expression profiles, and resulted in diverse expression patterns that may provide basic information to further reveal the function of the 49 RcbZIP genes in castor bean. The results obtained from this study would provide valuable information in understanding the molecular basis of the RcbZIP transcription factor family and their potential function in regulating the growth and development, particularly in seed filling of castor bean.

CEPiNS: Conserved Exon Prediction in Novel Species

Exon structure is relatively well conserved among orthologs in several large clades of species (e.g. Mammalia, Diptera, Lepidoptera) across evolutionary distances of up to 80 million years. Thus, it should be straightforward to predict the exon structures in novel species based upon the known exon structures of species that have had their genomes sequenced and well assembled. Being able to predict the exon boundaries in the genes of novel species is important given the quickly growing numbers of transcriptome sequencing projects. CEPiNS is a new pipeline for mining exon boundaries of predicted gene sets from model species and then using this information to identify the exon boundaries in a novel species through codon based alignment. The pipeline uses the freeware SPIDEY, an exon boundary prediction tool, and BLAST (BLASTN, BLASTP, TBLASTX), both of which are part of NCBI's toolkit. CEPiNS provides an important tool to analyze the transcriptome of novel species.

Evolutionary origin of β-defensins

β-Defensins are a group of vertebrate-specific antimicrobial peptides (AMPs) with microbicidal and immune regulatory functions. In spite of their conservation across the vertebrate lineage ranging from bony fish to human, the evolutionary origin of these molecules remains unsolved. We addressed this issue by comparing three-dimensional (3D) structure and genomic organization of β-defensins with those of big defensins, a family of invertebrate-derived β-defensin-related peptides with two distinct structural and functional domains. β-Defensins and the carboxyl-terminal domain of big defensins adopt a conserved β-sheet topology stabilized by three identical disulfide bridges. Genomic organization analysis revealed that the defensin domain of these two classes of molecules is encoded by a single exon with a positionally conserved phase-1 intron in its upstream. The genomic and 3D structural conservation provides convincing evidence for their evolutionary relationship, in which β-defensins emerged from an ancestral big defensin through exon shuffling or intronization of exonic sequences. The phylogenetic distribution of big defensins in Arthropoda, Mollusca and Cephalochordata suggests an early origin of the β-defensin domain, which can be traced to the common ancestor of bilateral metazoans.

The C. elegans rab family: identification, classification and toolkit construction

Rab monomeric GTPases regulate specific aspects of vesicle transport in eukaryotes including coat recruitment, uncoating, fission, motility, target selection and fusion. Moreover, individual Rab proteins function at specific sites within the cell, for example the ER, golgi and early endosome. Importantly, the localization and function of individual Rab subfamily members are often conserved underscoring the significant contributions that model organisms such as Caenorhabditis elegans can make towards a better understanding of human disease caused by Rab and vesicle trafficking malfunction. With this in mind, a bioinformatics approach was first taken to identify and classify the complete C. elegans Rab family placing individual Rabs into specific subfamilies based on molecular phylogenetics. For genes that were difficult to classify by sequence similarity alone, we did a comparative analysis of intron position among specific subfamilies from yeast to humans. This two-pronged approach allowed the classification of 30 out of 31 C. elegans Rab proteins identified here including Rab31/Rab50, a likely member of the last eukaryotic common ancestor (LECA). Second, a molecular toolset was created to facilitate research on biological processes that involve Rab proteins. Specifically, we used Gateway-compatible C. elegans ORFeome clones as starting material to create 44 full-length, sequence-verified, dominant-negative (DN) and constitutive active (CA) rab open reading frames (ORFs). Development of this toolset provided independent research projects for students enrolled in a research-based molecular techniques course at California State University, East Bay (CSUEB).

3D profile-based approach to proteome-wide discovery of novel human chemokines

Chemokines are small secreted proteins with important roles in immune responses. They consist of a conserved three-dimensional (3D) structure, so-called IL8-like chemokine fold, which is supported by disulfide bridges characteristic of this protein family. Sequence- and profile-based computational methods have been proficient in discovering novel chemokines by making use of their sequence-conserved cysteine patterns. However, it has been recently shown that some chemokines escaped annotation by these methods due to low sequence similarity to known chemokines and to different arrangement of cysteines in sequence and in 3D. Innovative methods overcoming the limitations of current techniques may allow the discovery of new remote homologs in the still functionally uncharacterized fraction of the human genome. We report a novel computational approach for proteome-wide identification of remote homologs of the chemokine family that uses fold recognition techniques in combination with a scaffold-based automatic mapping of disulfide bonds to define a 3D profile of the chemokine protein family. By applying our methodology to all currently uncharacterized human protein sequences, we have discovered two novel proteins that, without having significant sequence similarity to known chemokines or characteristic cysteine patterns, show strong structural resemblance to known anti-HIV chemokines. Detailed computational analysis and experimental structural investigations based on mass spectrometry and circular dichroism support our structural predictions and highlight several other chemokine-like features. The results obtained support their functional annotation as putative novel chemokines and encourage further experimental characterization. The identification of remote homologs of human chemokines may provide new insights into the molecular mechanisms causing pathologies such as cancer or AIDS, and may contribute to the development of novel treatments. Besides, the genome-wide applicability of our methodology based on 3D protein family profiles may open up new possibilities for improving and accelerating protein function annotation processes.

LB3D: a protein three-dimensional substructure search program based on the lower bound of a root mean square deviation value

Searching for protein structure-function relationships using three-dimensional (3D) structural coordinates represents a fundamental approach for determining the function of proteins with unknown functions. Since protein structure databases are rapidly growing in size, the development of a fast search method to find similar protein substructures by comparison of protein 3D structures is essential. In this article, we present a novel protein 3D structure search method to find all substructures with root mean square deviations (RMSDs) to the query structure that are lower than a given threshold value. Our new algorithm runs in O(m + N/m(0.5)) time, after O(N log N) preprocessing, where N is the database size and m is the query length. The new method is 1.8-41.6 times faster than the practically best known O(N) algorithm, according to computational experiments using a huge database (i.e., >20,000,000 C-alpha coordinates).

Structural basis for the dual recognition of helical cytokines IL-34 and CSF-1 by CSF-1R

Lacking any discernible sequence similarity, interleukin-34 (IL-34) and colony stimulating factor 1 (CSF-1) signal through a common receptor CSF-1R on cells of mononuclear phagocyte lineage. Here, the crystal structure of dimeric IL-34 reveals a helical cytokine fold homologous to CSF-1, and we further show that the complex architecture of IL-34 bound to the N-terminal immunoglobulin domains of CSF-1R is similar to the CSF-1/CSF-1R assembly. However, unique conformational adaptations in the receptor domain geometry and intermolecular interface explain the cross-reactivity of CSF-1R for two such distantly related ligands. The docking adaptations of the IL-34 and CSF-1 quaternary complexes, when compared to the stem cell factor assembly, draw a common evolutionary theme for transmembrane signaling. In addition, the structure of IL-34 engaged by a Fab fragment reveals the mechanism of a neutralizing antibody that can help deconvolute IL-34 from CSF-1 biology, with implications for therapeutic intervention in diseases with myeloid pathogenic mechanisms.

Molecular evolution of the neuropeptide S receptor

The neuropeptide S receptor (NPSR) is a recently deorphanized member of the G protein-coupled receptor (GPCR) superfamily and is activated by the neuropeptide S (NPS). NPSR and NPS are widely expressed in central nervous system and are known to have crucial roles in asthma pathogenesis, locomotor activity, wakefulness, anxiety and food intake. The NPS-NPSR system was previously thought to have first evolved in the tetrapods. Here we examine the origin and the molecular evolution of the NPSR using in-silico comparative analyses and document the molecular basis of divergence of the NPSR from its closest vertebrate paralogs. In this study, NPSR-like sequences have been identified in a hemichordate and a cephalochordate, suggesting an earlier emergence of a NPSR-like sequence in the metazoan lineage. Phylogenetic analyses revealed that the NPSR is most closely related to the invertebrate cardioacceleratory peptide receptor (CCAPR) and the group of vasopressin-like receptors. Gene structure features were congruent with the phylogenetic clustering and supported the orthology of NPSR to the invertebrate NPSR-like and CCAPR. A site-specific analysis between the vertebrate NPSR and the well studied paralogous vasopressin-like receptor subtypes revealed several putative amino acid sites that may account for the observed functional divergence between them. The data can facilitate experimental studies aiming at deciphering the common features as well as those related to ligand binding and signal transduction processes specific to the NPSR.

Identification of gene-oriented exon orthology between human and mouse

BACKGROUND

Gene orthology has been well studied in the evolutionary area and is thought to be an important implication to functional genome annotations. As the accumulation of transcriptomic data, alternative splicing is taken into account in the assignments of gene orthologs and the orthology is suggested to be further considered at transcript level. Whether gene or transcript orthology, exons are the basic units that represent the whole gene structure; however, there is no any reported study on how to build exon level orthology in a whole genome scale. Therefore, it is essential to establish a gene-oriented exon orthology dataset.

RESULTS

Using a customized pipeline, we first build exon orthologous relationships from assigned gene orthologs pairs in two well-annotated genomes: human and mouse. More than 92% of non-overlapping exons have at least one ortholog between human and mouse and only a small portion of them own more than one ortholog. The exons located in the coding region are more conserved in terms of finding their ortholog counterparts. Within the untranslated region, the 5' UTR seems to have more diversity than the 3' UTR according to exon orthology designations. Interestingly, most exons located in the coding region are also conserved in length but this conservation phenomenon dramatically drops down in untranslated regions. In addition, we allowed multiple assignments in exon orthologs and a subset of exons with possible fusion/split events were defined here after a thorough analysis procedure.

CONCLUSIONS

Identification of orthologs at the exon level is essential to provide a detailed way to interrogate gene orthology and splicing analysis. It could be used to extend the genome annotation as well. Besides examining the one-to-one orthologous relationship, we manage the one-to-multi exon pairs to represent complicated exon generation behavior. Our results can be further applied in many research fields studying intron-exon structure and alternative/constitutive exons in functional genomic areas.

In silico characterization and comparative genomic analysis of the Culex quinquefasciatus glutathione S-transferase (GST) supergene family

The glutathione S-transferases (GSTs) are phase II class of detoxification enzymes that are involved both directly and indirectly in insecticide resistance mechanisms. The Culex quinquefasciatus GST superfamily was analyzed by utilizing the public domain Culex genome sequence. In total, 35 cytosolic (seven classes) and 5 microsomal putatively active GSTs were retrieved, classified, and annotated. The study revealed the presence of three unclassified GSTs. Of 35 cytosolic GSTs, 65% contributed by insect specific Delta-Epsilon classes. Gene cluster analysis revealed that most of the genes of Delta, Epsilon, and Theta classes were organized into gene clusters. The gene organization analysis revealed the dominance of phase "0" introns in the Culex GST family. The studies on intron loss and gain events revealed that the Delta GSTs have experienced a higher number of loss and gains during their evolution. A positive correlation was observed between the phylogenetic relationship of members of the GST superfamily and their corresponding exon-intron organization. In addition, the genes within the gene clusters revealed the monophyletic phylogenetic relationship implying the importance of gene duplication events in the gene families' evolution. Finally, the comparative genomic analysis has shown a complex evolutionary scenario associated with the GST supergene family evolution in insects.

Genome-wide identification and evolutionary analysis of Arabidopsis sm genes family

Sm proteins are members of a family of small proteins that are widespread in biosphere and found associated with RNA metabolism. To date, to our knowledge, only Arabidopsis SAD1 gene has been studied functionally in plant. In this study, 42 Sm genes are identified through comprehensive analysis in Arabidopsis. And a complete overview of this gene family is presented, including the gene structures, phylogeny, chromosome locations, selection pressure and expression. The results reveal that gene duplication contributes to the expansion of the Sm gene family in Arabidopsis genome, diverse expression patterns suggest their functional differentiation and divergence analysis indicates purifying selection as a key role in evolution. Our comparative genomics analysis of Sm genes will provide the first step towards the future experimental research on determining the functions of these genes.

Powdery mildew fungal effector candidates share N-terminal Y/F/WxC-motif

Background

Powdery mildew and rust fungi are widespread, serious pathogens that depend on developing haustoria in the living plant cells. Haustoria are separated from the host cytoplasm by a plant cell-derived extrahaustorial membrane. They secrete effector proteins, some of which are subsequently transferred across this membrane to the plant cell to suppress defense.

Results

In a cDNA library from barley epidermis containing powdery mildew haustoria, two-thirds of the sequenced ESTs were fungal and represented ~3,000 genes. Many of the most highly expressed genes encoded small proteins with N-terminal signal peptides. While these proteins are novel and poorly related, they do share a three-amino acid motif, which we named "Y/F/WxC", in the N-terminal of the mature proteins. The first amino acid of this motif is aromatic: tyrosine, phenylalanine or tryptophan, and the last is always cysteine. In total, we identified 107 such proteins, for which the ESTs represent 19% of the fungal clones in our library, suggesting fundamental roles in haustoria function. While overall sequence similarity between the powdery mildew Y/F/WxC-proteins is low, they do have a highly similar exon-intron structure, suggesting they have a common origin. Interestingly, searches of public fungal genome and EST databases revealed that haustoria-producing rust fungi also encode large numbers of novel, short proteins with signal peptides and the Y/F/WxC-motif. No significant numbers of such proteins were identified from genome and EST sequences from either fungi which do not produce haustoria or from haustoria-producing Oomycetes.

Conclusion

In total, we identified 107, 178 and 57 such Y/F/WxC-proteins from the barley powdery mildew, the wheat stem rust and the wheat leaf rust fungi, respectively. All together, our findings suggest the Y/F/WxC-proteins to be a new class of effectors from haustoria-producing pathogenic fungi.

Antimicrobial peptide-like genes in Nasonia vitripennis: a genomic perspective

BACKGROUND

Antimicrobial peptides (AMPs) are an essential component of innate immunity which can rapidly respond to diverse microbial pathogens. Insects, as a rich source of AMPs, attract great attention of scientists in both understanding of the basic biology of the immune system and searching molecular templates for anti-infective drug design. Despite a large number of AMPs have been identified from different insect species, little information in terms of these peptides is available from parasitic insects.

RESULTS

By using integrated computational approaches to systemically mining the Hymenopteran parasitic wasp Nasonia vitripennis genome, we establish the first AMP repertoire whose members exhibit extensive sequence and structural diversity and can be distinguished into multiple molecular types, including insect and fungal defensin-like peptides (DLPs) with the cysteine-stabilized alpha-helical and beta-sheet (CSalphabeta) fold; Pro- or Gly-rich abaecins and hymenoptaecins; horseshoe crab tachystatin-type AMPs with the inhibitor cystine knot (ICK) fold; and a linear alpha-helical peptide. Inducible expression pattern of seven N. vitripennis AMP genes were verified, and two representative peptides were synthesized and functionally identified to be antibacterial. In comparison with Apis mellifera (Hymenoptera) and several non-Hymenopteran model insects, N. vitripennis has evolved a complex antimicrobial immune system with more genes and larger protein precursors. Three classical strategies that are likely responsible for the complexity increase have been recognized: 1) Gene duplication; 2) Exon duplication; and 3) Exon-shuffling.

CONCLUSION

The present study established the N. vitripennis peptidome associated with antimicrobial immunity by using a combined computational and experimental strategy. As the first AMP repertoire of a parasitic wasp, our results offer a basic platform for further studying the immunological and evolutionary significances of these newly discovered AMP-like genes in this class of insects.

Did cathelicidins, a family of multifunctional host-defense peptides, arise from a cysteine protease inhibitor?

Cystatins, the cysteine protease inhibitors, and the cathelin-like domain (CLD) of the antimicrobial cathelicidins are classified into the same superfamily because of their overall structural similarity. However, their evolutionary relationship has remained obscure owing to low sequence similarity. Structural similarity of two proteins often provides evidence for divergent evolution; however, structural convergence can not be completely ruled out in this case. Conserved gene structure and related function provide new evidence in favor of a common ancestral origin for cystatins and CLDs. On the basis of two observations, the C-terminal location of the cathelicidin antimicrobial domain and evolutionary gain of one 3' intron, I propose a gradual evolution model to explain how the AMD evolved from the ancestral cystatin scaffold.

Exalign: a new method for comparative analysis of exon-intron gene structures

The evolution of genes is usually studied and reconstructed at the sequence level, that is, by comparing and aligning their genomic, transcript or protein sequences. However, including the exon–intron structure of genes in the analysis can provide further and useful information, for example to draw reliable phylogenetic relationships left unsolved by traditional sequence-based evolutionary studies, or to shed further light on patterns of intron gain and loss. In spite of this, no tool especially devised for this task is currently available. In this work we present Exalign, an algorithm designed to retrieve, compare and search for the exon-intron structure of existing gene annotations, that has been implemented in a software tool freely accessible through a web interface as well as available for download. We present different applications of our method, from the reconstruction of the evolutionary history of homologous gene families to the detection of as of today unknown cases of intron loss in human and rodents, and, remarkably, two never reported intron gain events in human and mouse. The web interface for accessing Exalign is available at http://www.pesolelab.it/exalign/ or http://www.beacon.unimi.it/exalign/

The primate-specific protein TBC1D3 is required for optimal macropinocytosis in a novel ARF6-dependent pathway

The generation of novel genes and proteins throughout evolution has been proposed to occur as a result of whole genome and gene duplications, exon shuffling, and retrotransposition events. The analysis of such genes might thus shed light into the functional complexity associated with highly evolved species. One such case is represented by TBC1D3, a primate-specific gene, harboring a TBC domain. Because TBC domains encode Rab-specific GAP activities, TBC-containing proteins are predicted to play a major role in endocytosis and intracellular traffic. Here, we show that the TBC1D3 gene originated late in evolution, likely through a duplication of the RNTRE locus, and underwent gene amplification during primate speciation. Despite possessing a TBC domain, TBC1D3 is apparently devoid of Rab-GAP activity. However, TBC1D3 regulates the optimal rate of epidermal growth factor-mediated macropinocytosis by participating in a novel pathway involving ARF6 and RAB5. In addition, TBC1D3 binds and colocalize to GGA3, an ARF6-effector, in an ARF6-dependent manner, and synergize with it in promoting macropinocytosis, suggesting that the two proteins act together in this process. Accordingly, GGA3 siRNA-mediated ablation impaired TBC1D3-induced macropinocytosis. We thus uncover a novel signaling pathway that appeared after primate speciation. Within this pathway, a TBC1D3:GGA3 complex contributes to optimal propagation of signals, ultimately facilitating the macropinocytic process.

Phylogenetic analysis of vertebrate kininogen genes

Kininogens, the precursors of bradykinins, vary extremely in both structure and function among different taxa of animals, in particular between mammals and amphibians. This includes even the most conserved bradykinin domain in terms of biosynthesis mode and structure. To elucidate the evolutionary dynamics of kininogen genes, we have identified 19 novel amino acid sequences from EST and genomic databases (for mammals, birds, and fishes) and explored their phylogenetic relationships using combined amino acid sequence and gene structure as markers. Our results show that there were initially two paralogous kininogen genes in vertebrates. During their evolution, the original gene was saved with frequent multiplication in amphibians, but lost in fishes, birds, and mammals, while the novel gene was saved with multiple functions in fishes, birds, and mammals, but became a pseudogene in amphibians. We also propose that the defense mechanism against specific predators in amphibian skin secretions has been bradykinin receptor dependent. Our findings may provide a foundation for identification and structural, functional, and evolutionary analyses of more kininogen genes and other gene families.

Protein arginine methyltransferases: Evolution and assessment of their pharmacological and therapeutic potential

Protein arginine N-methylation is a post-translational modification whose influence on cell function is becoming widely appreciated. Protein arginine methyltransferases (PRMT) catalyze the methylation of terminal nitrogen atoms of guanidinium side chains within arginine residues of proteins. Recently, several new members of the PRMT family have been cloned and their catalytic function determined. In this report, we present a review and phylogenetic analysis of the PRMT found so far in genomes. PRMT are found in nearly all groups of eukaryotes. Many human PRMT originated early in eukaryote evolution. Homologs of PRMT1 and PRMT5 are found in nearly every eukaryote studied. The gene structure of PRMT vary: most introns appear to be inserted randomly into the open reading frame. The change in catalytic specificity of some PRMT occurred with changes in the arginine binding pocket within the active site. Because of the high degree of conservation of sequence among the family throughout evolution, creation of specific PRMT inhibitors in pathogenic organisms may be difficult, but could be very effective if developed. Furthermore, because of the intricate involvement of several PRMT in cellular physiology, their inhibition may be fraught with unwanted side effects. Nevertheless, development of pharmaceutical agents to control PRMT functions could lead to significant new targets.