PIRLs: a novel class of plant intracellular leucine-rich repeat proteins

Leucine-rich repeat (LRR) proteins feature tandem leucine-rich motifs that form a protein-protein interaction domain. Plants contain diverse classes of LRR proteins, many of which take part in signal transduction. We have identified a novel family of nine Arabidopsis LRR proteins that, based on predicted intracellular location and LRR motif consensus sequence, are related to Ras-binding LRR proteins found in signaling complexes in animals and yeast. This new class has been named plant intracellular Ras group-related LRR proteins (PIRLs). We have characterized PIRL cDNAs, rigorously defined gene and protein annotations, investigated gene family evolution and surveyed mRNA expression. While LRR regions suggested a relationship to Ras group LRR proteins, outside of their LRR domains PIRLs differed from Ras group proteins, exhibiting N- and C-terminal regions containing low complexity stretches and clusters of charged amino acids. PIRL genes grouped into three subfamilies based on sequence relationships and gene structures. Related gene pairs and dispersed chromosomal locations suggested family expansion by ancestral genomic or segmental duplications. Expression surveys revealed that all PIRL mRNAs are actively transcribed, with three expressed differentially in leaves, roots or flowers. These results define PIRLs as a distinct, plant-specific class of intracellular LRR proteins that probably mediate protein interactions, possibly in the context of signal transduction. T-DNA knock-out mutants have been isolated as a starting point for systematic functional analysis of this intriguing family.     

Highly dynamic exon shuffling in candidate pathogen receptors ... what if brown algae were capable of adaptive immunity?

Pathogen recognition is the first step of immune reactions. In animals and plants, direct or indirect pathogen recognition is often mediated by a wealth of fast-evolving receptors, many of which contain ligand-binding and signal transduction domains, such as leucine-rich or tetratricopeptide repeat (LRR/TPR) and NB-ARC domains, respectively. In order to identify candidates potentially involved in algal defense, we mined the genome of the brown alga Ectocarpus siliculosus for homologues of these genes and assessed the evolutionary pressures acting upon them. We thus annotated all Ectocarpus LRR-containing genes, in particular an original group of LRR-containing GTPases of the ROCO family, and 24 NB-ARC-TPR proteins. They exhibit high birth and death rates, while a diversifying selection is acting on their LRR (respectively TPR) domain, probably affecting the ligand-binding specificities. Remarkably, each repeat is encoded by an exon, and the intense exon shuffling underpins the variability of LRR and TPR domains. We conclude that the Ectocarpus ROCO and NB-ARC-TPR families are excellent candidates for being involved in recognition/transduction events linked to immunity. We further hypothesize that brown algae may generate their immune repertoire via controlled somatic recombination, so far only known from the vertebrate adaptive immune systems.     

Leucine‐rich repeat receptor‐like kinase II phylogenetics reveals five main clades throughout the plant kingdom

Receptor‐like kinases (RLKs) represent the largest group of cell surface receptors in plants. The monophyletic leucine‐rich repeat (LRR)‐RLK subfamily II is considered to contain the somatic embryogenesis receptor kinases (SERKs) and NSP‐interacting kinases known to be involved in developmental processes and cellular immunity in plants. There are only a few published studies on the phylogenetics of LRR‐RLKII; unfortunately these suffer from poor taxon/gene sampling. Hence, it is not clear how many and what main clades this family contains, let alone what structure–function relationships exist. We used 1342 protein sequences annotated as ‘SERK’ and ‘SERK‐like’ plus related sequences in order to estimate phylogeny within the LRR‐RLKII clade, using the nematode protein kinase Pelle as an outgroup. We reconstruct five main clades (LRR‐RLKII 1–5), in each of which the main pattern of land plant relationships re‐occurs, confirming previous hypotheses that duplication events happened in this gene subfamily prior to divergence among land plant lineages. We show that domain structures and intron–exon boundaries within the five clades are well conserved in evolution. Furthermore, phylogenetic patterns based on the separate LRR and kinase parts of LRR‐RLKs are incongruent: whereas the LRR part supports a LRR‐RLKII 2/3 sister group relationship, the kinase part supports clades 1/2. We infer that the kinase part includes few ‘radical’ amino acid changes compared with the LRR part. Finally, our results confirm that amino acids involved in each LRR‐RLKII–receptor complex interaction are located at N‐capping residues, and that the short amino acid motifs of this interaction domain are highly conserved throughout evolution within the five LRR‐RLKII clades.     

Correlation of gene and protein structures in the FXYD family proteins

The FXYD family proteins are auxiliary subunits of the Na,K-ATPase, expressed primarily in tissues that specialize in fluid or solute transport, or that are electrically excitable. These proteins range in size from about 60 to 160 amino acid residues, and share a core homology of 35 amino acid residues in and around a single transmembrane segment. Despite their relatively small sizes, they are all encoded by genes with six to nine small exons. We show that the helical secondary structures of three FXYD family members, FXYD1, FXYD3, and FXYD4, determined in micelles by NMR spectroscopy, reflect the structures of their corresponding genes. The coincidence of helical regions, and connecting segments, with the positions of intron-exon junctions in the genes, support the hypothesis that the FXYD proteins may have been assembled from discrete structural modules through exon shuffling.     

Gene expansion and retention leads to a diverse tyrosine kinase superfamily in amphioxus

Tyrosine kinase (TK) proteins play a central role in cellular behavior and development of animals. The expansion of this superfamily is regarded as a key event in the evolution of the complex signaling pathways and gene networks of metazoans and is a prominent example of how shuffling of protein modules may generate molecular novelties. Using the intron/exon structure within the TK domain (TK intron code) as a complementary tool for the assignment of orthology and paralogy, we identified and studied the 118 TK proteins of the amphioxus Branchiostoma floridae genome to elucidate TK gene family evolution in metazoans and chordates in particular. Unlike all characterized metazoans to date, amphioxus has members of all known widespread TK families, with not a single loss. Putting amphioxus TKs in an evolutionary context, including new data from the cnidarian Nematostella vectensis, the echinoderm Strongylocentrotus purpuratus, and the ascidian Ciona intestinalis, we suggest new evolutionary histories for different TK families and draw a new global picture of gene loss/gain in the different phyla. Surprisingly, our survey also detected an unprecedented expansion of a group of closely related TK families, including TIE, FGFR, PDGFR, and RET, due most probably to massive gene duplication and exon shuffling. Based on their highly similar intron/exon structure at the TK domain, we suggest that this group of TK families constitute a superfamily of TK proteins, which we termed EXpanding TK, after their seemingly unique propensity to gene duplication and exon shuffling, not only in amphioxus but also across all metazoan groups. Due to this extreme tendency to both retention and expansion of TK genes, amphioxus harbors the richest and most diverse TK repertoire among all metazoans studied so far, retaining most of the gene complement of its ancestors, but having evolved its own repertoire of genetic novelties.     

Isolation of a novel ras gene from Trichomonas vaginalis: a possible evolutionary ancestor of the Ras and Rap genes of higher eukaryotes.

The Ras subfamily proteins are small, monomeric GTP-binding proteins with vital roles in regulating eukaryotic signal transduction pathways. Gene duplication and divergence have been postulated as the mechanism by which such family members have evolved their specific functions. A cDNA clone of TvRsp was isolated and sequenced from a cDNA expression library of the primitive eukaryote Trichomonas vaginalis. The genomic DNA corresponding to the cDNA sequence was amplified by PCR and sequenced. Sequence analysis suggested that TvRsp was an intronless gene. This gene encoded a protein of 181 amino acids and contained the 5 conserved G domains that designated it as a Ras or Rap subfamily member. However, the deduced amino acid sequence shared only 34%-37% overall identity with other Ras subfamily members of different species, and the presence of motifs characteristic of both the Ras and Rap families of GTPase confused the familial classification of this gene. Phylogenetic analysis showed its origins at the divergence point of the Ras/Rap families and suggested that TvRsp was a possible evolutionary ancestral gene of the ras/rap genes of higher eukaryotes. This information was of importance not only from the perspective of understanding the evolution and diversity of eukaryotic signal transduction pathways but also in providing a framework by which to understand protein processing in the growth and differentiation of single-celled microorganisms.     

Computational identification of novel chitinase-like proteins in the Drosophila melanogaster genome

MOTIVATION: Multiple chitinases as well as lectins closely related to them have been characterized previously from many insect species and the corresponding genes/cDNAs have been cloned. However, the identification of the entire assortment of genes for chitinase family proteins and their differences in biochemical properties have not been carried out in any individual insect species. The completion of the entire DNA sequence of Drosophila melanogaster (fruit fly) genome and identification of open reading frames presents an opportunity to study the structures and functions of chitinase-like proteins, and also to identify new members of this family in DROSOPHILA: We are, therefore, interested in studying the functional genomics of chitinase-like gene families in insects. METHODS: We searched the Drosophila protein sequences database using fully characterized insect chitinase sequences and BLASTP software, identified all the putative chitinase-like proteins encoded in Drosophila genome, and predicted their structures using domain analysis tools. A phylogenetic analysis of the chitinase-like proteins from Drosophila and several other insect species was carried out. The structures of these chitinases were modeled using homology modeling software. RESULTS: Our analysis revealed the presence of 18 chitinase-like proteins in the Drosophila protein database. Among these are seven novel chitinase-like proteins that contain four signature amino acid sequences of chitinases belonging to family 18 glycosylhydrolases, including both acidic and hydrophobic amino acid residues critical for enzyme activity. All the proteins contain at least one catalytic domain with one having four catalytic domains. Phylogenetic analysis of chitinase-like proteins from Drosophila and other insects revealed an evolutionary relationship among all these proteins, which indicated gene duplication and domain shuffling to generate the observed diversity in the encoded proteins. Homology modeling showed that all the Drosophila chitinase-like proteins contain one or more catalytic domains with a (alpha/beta)8 barrel-like structure. Our results suggest that insects utilize multiple family 18 chitinolytic enzymes and also non-enzymatic chitinase-like proteins for degrading/remodeling/binding to chitin in different insect anatomical extracellular structures, such as the cuticle, peritrophic membrane, trachea and mouth parts during insect development, and possibly for other roles including chitin synthesis. AVAILABILITY: Perl program and supplementary material are available at http://www.ksu.edu/bioinformatics/supplementary.htm     

The diversification of plant cytosolic small heat shock proteins preceded the divergence of mosses.

A cDNA library was constructed with mRNA isolated from heat-stressed cell cultures of Funaria hygrometrica (Bryophyta, Musci, Funariaceae). cDNA clones encoding six cytosolic small heat shock proteins (sHSPs) were identified using differential screening. Phylogenetic analysis of these sHSP sequences with other known sHSPs identified them as members of the previously described higher plant cytosolic class I and II families. Four of the F. hygrometrica sHSPs are members of the cytosolic class I family, and the other two are members of the cytosolic class II family. The presence of members of the cytosolic I and II sHSP families in a bryophyte indicates that these gene families are ancient, and evolved at least 450 MYA. This result also indicates that the plant sHSP gene families duplicated much earlier than did the well-studied phytochrome gene family. Members of the cytosolic I and II sHSP families are developmentally regulated in seeds and flowers in higher plants. Our findings show that the two cytosolic sHSP families evolved before the appearance of these specialized structures. Previous analysis of angiosperm sHSPs had identified class- or family-specific amino acid consensus regions and determined that rate heterogeneity exists among the different sHSP families. The analysis of the F. hygrometrica sHSP sequences reveals patterns and rates of evolution distinct from those seen among angiosperm sHSPs. Some, but not all, of the amino acid consensus regions identified in seed plants are conserved in the F. hygrometrica sHSPs. Taken together, the results of this study illuminate the evolution of the sHSP gene families and illustrate the importance of including representatives of basal land plant lineages in plant molecular evolutionary studies.     

Fliih, a gelsolin-related cytoskeletal regulator essential for early mammalian embryonic development

The Drosophila melanogaster flightless I gene is required for normal cellularization of the syncytial blastoderm. Highly conserved homologues of flightless I are present in Caenorhabditis elegans, mouse, and human. We have disrupted the mouse homologue Fliih by homologous recombination in embryonic stem cells. Heterozygous Fliih mutant mice develop normally, although the level of Fliih protein is reduced. Cultured homozygous Fliih mutant blastocysts hatch, attach, and form an outgrowing trophoblast cell layer, but egg cylinder formation fails and the embryos degenerate. Similarly, Fliih mutant embryos initiate implantation in vivo but then rapidly degenerate. We have constructed a transgenic mouse carrying the complete human FLII gene and shown that the FLII transgene is capable of rescuing the embryonic lethality of the homozygous targeted Fliih mutation. These results confirm the specific inactivation of the Fliih gene and establish that the human FLII gene and its gene product are functional in the mouse. The Fliih mouse mutant phenotype is much more severe than in the case of the related gelsolin family members gelsolin, villin, and CapG, where the homozygous mutant mice are viable and fertile but display alterations in cytoskeletal actin regulation.     

Polygalacturonase-inhibiting proteins--leucine-rich repeat proteins in plant defence

Plant polygalacturonase-inhibiting proteins (PGIPs) belong to the leucine-rich repeat (LRR) family and are known to prevent pathogen invasion by inhibiting the plant cell wall degrading enzyme, polygalacturonase. Our study reveals that these multigene-encoded defence proteins found in flowering plants only exhibit identical domain architecture with 10 tandemly-arranged LRRs. This implies that variations of PGIP inhibitory properties are not associated with the number of the repeats but with subtle changes in the sequence content of the repeats. The first and eighth repeat contain more mutations compared to the strict conservation of the plant-specific LRRs or any repeat at other positions. Each of these repeats forms a separate cluster in the phylogenetic tree, both within and across plant families, thus suggesting uniqueness with respect to their position. A study of the genes encoding PGIPs, shows the existence of two categories (i) single exon and hence no intron; and (ii) two exons with an intron in between. Analyses of the intron phase and correlation of the exon-intron structure with the compact structural modules in PGIPs support insertion of introns in the pre-existing single exon genes and thus the intron late model. Lack of conservation of phase across families and formation of individual clusters for each family in the phylogenetic tree drawn with the intron sequences illustrate the event of insertion that took place separately in each of these families.     

A screen of cell-surface molecules identifies leucine-rich repeat proteins as key mediators of synaptic target selection

In Drosophila embryos and larvae, a small number of identified motor neurons innervate body wall muscles in a highly stereotyped pattern. Although genetic screens have identified many proteins that are required for axon guidance and synaptogenesis in this system, little is known about the mechanisms by which muscle fibers are defined as targets for specific motor axons. To identify potential target labels, we screened 410 genes encoding cell-surface and secreted proteins, searching for those whose overexpression on all muscle fibers causes motor axons to make targeting errors. Thirty such genes were identified, and a number of these were members of a large gene family encoding proteins whose extracellular domains contain leucine-rich repeat (LRR) sequences, which are protein interaction modules. By manipulating gene expression in muscle 12, we showed that four LRR proteins participate in the selection of this muscle as the appropriate synaptic target for the RP5 motor neuron.     

SCAN domain-containing 2 gene (SCAND2) is a novel nuclear protein derived from the zinc finger family by exon shuffling

The SCAN domain is a recently recognized protein domain that characterizes a subfamily of the Krüppel-like zinc finger proteins. We have previously described a novel SCAN domain-containing 2 gene (SCAND2) that does not belong to the zinc finger family. We report structural and sequence analyzes of all known members of the SCAN family and use these data to illustrate a model of gene family evolution. Most of the SCAN containing genes share common gene organization features that support the proposed origin for SCAND2 by disruption of an ancestral SCAN-zinc finger gene by a retroposition event and subsequent exon shuffling.     

Designing repeat proteins: modular leucine-rich repeat protein libraries based on the mammalian ribonuclease inhibitor family

We present a novel approach to design repeat proteins of the leucine-rich repeat (LRR) family for the generation of libraries of intracellular binding molecules. From an analysis of naturally occurring LRR proteins, we derived the concept to assemble repeat proteins with randomized surface positions from libraries of consensus repeat modules. As a guiding principle, we used the mammalian ribonuclease inhibitor (RI) family, which comprises cytosolic LRR proteins known for their extraordinary affinities to many RNases. By aligning the amino acid sequences of the internal repeats of human, pig, rat, and mouse RI, we derived a first consensus sequence for the characteristic alternating 28 and 29 amino acid residue A-type and B-type repeats. Structural considerations were used to replace all conserved cysteine residues, to define less conserved positions, and to decide where to introduce randomized amino acid residues. The so devised consensus RI repeat library was generated at the DNA level and assembled by stepwise ligation to give libraries of 2-12 repeats. Terminal capping repeats, known to shield the continuous hydrophobic core of the LRR domain from the surrounding solvent, were adapted from human RI. In this way, designed LRR protein libraries of 4-14 LRRs (equivalent to 130-415 amino acid residues) were obtained. The biophysical analysis of randomly chosen library members showed high levels of soluble expression in the Escherichia coli cytosol, monomeric behavior as characterized by gel-filtration, and alpha-helical CD spectra, confirming the success of our design approach.     

Phylogeny and domain evolution in the APETALA2-like gene family

The combined processes of gene duplication, nucleotide substitution, domain duplication, and intron/exon shuffling can generate a complex set of related genes that may differ substantially in their expression patterns and functions. The APETALA2-like (AP2-like) gene family exhibits patterns of both gene and domain duplication, coupled with changes in sequence, exon arrangement, and expression. In angiosperms, these genes perform an array of functions including the establishment of the floral meristem, the specification of floral organ identity, the regulation of floral homeotic gene expression, the regulation of ovule development, and the growth of floral organs. To determine patterns of gene diversification, we conducted a series of broad phylogenetic analyses of AP2-like sequences from green plants. These studies indicate that the AP2 domain was duplicated prior to the divergence of the two major lineages of AP2-like genes, euAP2 and AINTEGUMENTA (ANT). Structural features of the AP2-like genes as well as phylogenetic analyses of nucleotide and amino acid (aa) sequences of the AP2-like gene family support the presence of the two major lineages. The ANT lineage is supported by a 10-aa insertion in the AP2-R1 domain and a 1-aa insertion in the AP2-R2 domain, relative to all other members of the AP2-like family. MicroRNA172-binding sequences, the function of which has been studied in some of the AP2-like genes in Arabidopsis, are restricted to the euAP2 lineage. Within the ANT lineage, the euANT lineage is characterized by four conserved motifs: one in the 10-aa insertion in the AP2-R1 domain (euANT1) and three in the predomain region (euANT2, euANT3, and euANT4). Our expression studies show that the euAP2 homologue from Amborella trichopoda, the putative sister to all other angiosperms, is expressed in all floral organs as well as leaves.     

Receptor-like genes in the major resistance locus of lettuce are subject to divergent selection.

Disease resistance genes in plants are often found in complex multigene families. The largest known cluster of disease resistance specificities in lettuce contains the RGC2 family of genes. We compared the sequences of nine full-length genomic copies of RGC2 representing the diversity in the cluster to determine the structure of genes within this family and to examine the evolution of its members. The transcribed regions range from at least 7.0 to 13.1 kb, and the cDNAs contain deduced open reading frames of approximately 5. 5 kb. The predicted RGC2 proteins contain a nucleotide binding site and irregular leucine-rich repeats (LRRs) that are characteristic of resistance genes cloned from other species. Unique features of the RGC2 gene products include a bipartite LRR region with >40 repeats. At least eight members of this family are transcribed. The level of sequence diversity between family members varied in different regions of the gene. The ratio of nonsynonymous (Ka) to synonymous (Ks) nucleotide substitutions was lowest in the region encoding the nucleotide binding site, which is the presumed effector domain of the protein. The LRR-encoding region showed an alternating pattern of conservation and hypervariability. This alternating pattern of variation was also found in all comparisons within families of resistance genes cloned from other species. The Ka /Ks ratios indicate that diversifying selection has resulted in increased variation at these codons. The patterns of variation support the predicted structure of LRR regions with solvent-exposed hypervariable residues that are potentially involved in binding pathogen-derived ligands.