Molecular evolution of the cadherin superfamily

This review deals with the large and pleiotropic superfamily of cadherins and its molecular evolution. We compiled literature data and an in-depth phylogenetic analysis of more than 350 members of this superfamily from about 30 species, covering several but not all representative branches within metazoan evolution. We analyzed the sequence homology between either ectodomains or cytoplasmic domains, and we reviewed protein structural data and genomic architecture. Cadherins and cadherin-related molecules are defined by having an ectodomain in which at least two consecutive calcium-binding cadherin repeats are present. There are usually 5 or 6 domains, but in some cases as many as 34. Additional protein modules in the ectodomains point at adaptive evolution. Despite the occurrence of several conserved motifs in subsets of cytoplasmic domains, these domains are even more diverse than ectodomains and most likely have evolved separately from the ectodomains. By fine tuning molecular classifications, we reduced the number of solitary superfamily members. We propose a cadherin major branch, subdivided in two families and 8 subfamilies, and a cadherin-related major branch, subdivided in four families and 11 subfamilies. Accordingly, we propose a more appropriate nomenclature. Although still fragmentary, our insight into the molecular evolution of these remarkable proteins is steadily growing. Consequently, we can start to propose testable hypotheses for structure-function relationships with impact on our models of molecular evolution. An emerging concept is that the ever evolving diversity of cadherin structures is serving dual and important functions: specific cell adhesion and intricate cell signaling.     

Searching for natural lineages within the Cerylonid Series (Coleoptera: Cucujoidea)

Phylogenetic relationships within the diverse beetle superfamily Cucujoidea are poorly known. The Cerylonid Series (C.S.) is the largest of all proposed superfamilial cucujoid groups, comprising eight families and representing most of the known cucujoid species diversity. The monophyly of the C.S., however, has never been formally tested and the higher-level relationships among and within the constituent families remain equivocal. Here we present a phylogenetic study based on 18S and 28S rDNA for 16 outgroup taxa and 61 C.S. ingroup taxa, representing seven of the eight C.S. families and 20 of 39 subfamilies. We test the monophyly of the C.S., investigate the relationships among the C.S. families, and test the monophyly of the constituent families and subfamilies. Phylogenetic reconstruction of the combined data was achieved via standard static alignment parsimony analyses, Direct Optimization using parsimony, and partitioned Bayesian analysis. All three analyses support the paraphyly of Cucujoidea with respect to Tenebrionoidea and confirm the monophyly of the C.S. The C.S. families Bothrideridae, Cerylonidae, Discolomatidae, Coccinellidae and Corylophidae are supported as monophyletic in all analyses. Only the Bayesian analysis recovers a monophyletic Latridiidae. Endomychidae is recovered as polyphyletic in all analyses. Of the 14 subfamilies with multiple terminals in this study, 11 were supported as monophyletic. The corylophid subfamily Corylophinae and the coccinellid subfamilies Chilocorinae and Scymninae are recovered as paraphyletic. A sister grouping of Anamorphinae+Corylophidae is supported in all analyses. Other taxonomic implications are discussed in light of our results.     

Probing the S100 protein family through genomic and functional analysis

The EF-hand superfamily of calcium binding proteins includes the S100, calcium binding protein, and troponin subfamilies. This study represents a genome, structure, and expression analysis of the S100 protein family, in mouse, human, and rat. We confirm the high level of conservation between mammalian sequences but show that four members, including S100A12, are present only in the human genome. We describe three new members of the S100 family in the three species and their locations within the S100 genomic clusters and propose a revised nomenclature and phylogenetic relationship between members of the EF-hand superfamily. Two of the three new genes were induced in bone-marrow-derived macrophages activated with bacterial lipopolysaccharide, suggesting a role in inflammation. Normal human and murine tissue distribution profiles indicate that some members of the family are expressed in a specific manner, whereas others are more ubiquitous. Structure-function analysis of the chemotactic properties of murine S100A8 and human S100A12, particularly within the active hinge domain, suggests that the human protein is the functional homolog of the murine protein. Strong similarities between the promoter regions of human S100A12 and murine S100A8 support this possibility. This study provides insights into the possible processes of evolution of the EF-hand protein superfamily. Evolution of the S100 proteins appears to have occurred in a modular fashion, also seen in other protein families such as the C2H2-type zinc-finger family.     

Structural Variations in Protein Superfamilies: Actin and Tubulin

Structures of homologous proteins are usually conserved during evolution, as are critical active site residues. This is the case for actin and tubulin, the two most important cytoskeleton proteins in eukaryotes. Actins and their related proteins (Arps) constitute a large superfamily whereas the tubulin family has fewer members. Unaligned sequences of these two protein families were analysed by searching for short groups of family-specific amino acid residues, that we call motifs, and by counting the number of residues from one motif to the next. For each sequence, the set of motif-to-motif residue counts forms a subfamily-specific pattern (landmark pattern) allowing actin and tubulin superfamily members to be identified and sorted into subfamilies. The differences between patterns of individual subfamilies are due to inserts and deletions (indels). Inserts appear to have arisen at an early stage in eukaryote evolution as suggested by the small but consistent kingdom-dependent differences found within many Arp subfamilies and in γ-tubulins. Inserts tend to be in surface loops where they can influence subfamily-specific function without disturbing the core structure of the protein. The relatively few indels found for tubulins have similar positions to established results, whereas we find many previously unreported indel positions and lengths for the metazoan Arps.     

Expanding the nitrogen regulatory protein superfamily: Homology detection at below random sequence identity

Nitrogen regulatory (PII) proteins are signal transduction molecules involved in controlling nitrogen metabolism in prokaryots. PII proteins integrate the signals of intracellular nitrogen and carbon status into the control of enzymes involved in nitrogen assimilation. Using elaborate sequence similarity detection schemes, we show that five clusters of orthologs (COGs) and several small divergent protein groups belong to the PII superfamily and predict their structure to be a (betaalphabeta)(2) ferredoxin-like fold. Proteins from the newly emerged PII superfamily are present in all major phylogenetic lineages. The PII homologs are quite diverse, with below random (as low as 1%) pairwise sequence identities between some members of distant groups. Despite this sequence diversity, evidence suggests that the different subfamilies retain the PII trimeric structure important for ligand-binding site formation and maintain a conservation of conservations at residue positions important for PII function. Because most of the orthologous groups within the PII superfamily are composed entirely of hypothetical proteins, our remote homology-based structure prediction provides the only information about them. Analogous to structural genomics efforts, such prediction gives clues to the biological roles of these proteins and allows us to hypothesize about locations of functional sites on model structures or rationalize about available experimental information. For instance, conserved residues in one of the families map in close proximity to each other on PII structure, allowing for a possible metal-binding site in the proteins coded by the locus known to affect sensitivity to divalent metal ions. Presented analysis pushes the limits of sequence similarity searches and exemplifies one of the extreme cases of reliable sequence-based structure prediction. In conjunction with structural genomics efforts to shed light on protein function, our strategies make it possible to detect homology between highly diverse sequences and are aimed at understanding the most remote evolutionary connections in the protein world.     

Phylogenetic analysis of the chromate ion transporter (CHR) superfamily

ChrA is a membrane protein that confers resistance to the toxic ion chromate through the energy-dependent chromate efflux from the cytoplasm. In the protein databases, ChrA is a member of the chromate ion transporter (CHR) superfamily, composed of at least several dozens of members, distributed in the three domains of life. The aim of this work was to perform a phylogenetic analysis of the CHR superfamily. An exhaustive search for ChrA homologous proteins was carried out at the National Center for Biotechnology Information database. One hundred and thirty-five sequences were identified as members of the CHR superfamily [77 long-chain sequences, or bidomains (LCHR), and 58 short-chain sequences, or monodomains (SCHR)], organized mainly as tandem pairs of genes whose resultant proteins probably possess oppositely oriented membrane topology. LCHR sequences were split into amino and carboxyl domains, and the resultant domains were aligned with the SCHR proteins. A phylogenetic tree was reconstructed using four different methods, obtaining similar results. The domains were grouped into three clusters: the SCHR proteins cluster, the amino domain cluster of LCHR proteins and the carboxyl domain cluster of LCHR proteins. These results, as well as differences in the genomic context of CHR proteins, enabled the proteins to be sorted into two families (SCHR and LCHR), and 10 subfamilies. Evidence was found suggesting an ancient origin of LCHR proteins from the fusion of two SCHR protein-encoding genes; however, some secondary events of fusion and fission may have occurred later. The separate distribution of the LCHR and SCHR proteins, differences in the genomic context in both groups and the fact that chromate transport has been demonstrated only in LCHR proteins suggest that the CHR proteins comprise two or more paralogous groups in the CHR superfamily.     

Phylogeny of Myrmeleontiformia based on larval morphology (Neuropterida: Neuroptera)

The suborder Myrmeleontiformia is a derived lineage of lacewings (Insecta: Neuroptera) including the families Psychopsidae, Nemopteridae, Nymphidae, Ascalaphidae and Myrmeleontidae. In particular, Myrmeleontidae (antlions) are the most diverse neuropteran family, representing a conspicuous component of the insect fauna of xeric environments. We present the first detailed quantitative phylogenetic analysis of Myrmeleontiformia, based on 107 larval morphological and behavioural characters for 36 genera whose larvae are known (including at least one representative of all the subfamilies of the suborder). Four related families were used as outgroups to polarize character states. Phylogenetic analyses were conducted using both parsimony and Bayesian methods. The reconstructions resulting from our analyses corroborate the monophyly of Myrmeleontiformia. Within this clade, Psychopsidae are recovered as the sister family to all the remaining taxa. Nemopteridae (including both subfamilies Nemopterinae and Crocinae) are recovered as monophyletic and sister to the clade comprising Nymphidae + (Myrmeleontidae + Ascalaphidae). Nymphidae consist of two well‐supported clades corresponding to the subfamilies Nymphinae and Myiodactylinae. Our results suggest that Ascalaphidae may not be monophyletic, as they collapse into an unresolved polytomy under the Bayesian analysis. In addition, the recovered phylogenetic relationships diverge from the traditional classification scheme for ascalaphids. Myrmeleontidae are reconstructed as monophyletic, with the subfamilies Stilbopteryginae, Palparinae and Myrmeleontinae. We retrieved a strongly supported clade comprising taxa with a fossorial habit of the preimaginal instars, which represents a major antlion radiation, also including the monophyletic pit‐trap building species.     

Diversity and dispersal of a ubiquitous protein family: acyl-CoA dehydrogenases

Acyl-CoA dehydrogenases (ACADs), which are key enzymes in fatty acid and amino acid catabolism, form a large, pan-taxonomic protein family with at least 13 distinct subfamilies. Yet most reported ACAD members have no subfamily assigned, and little is known about the taxonomic distribution and evolution of the subfamilies. In completely sequenced genomes from approximately 210 species (eukaryotes, bacteria and archaea), we detect ACAD subfamilies by rigorous ortholog identification combining sequence similarity search with phylogeny. We then construct taxonomic subfamily-distribution profiles and build phylogenetic trees with orthologous proteins. Subfamily profiles provide unparalleled insight into the organisms’ energy sources based on genome sequence alone and further predict enzyme substrate specificity, thus generating explicit working hypotheses for targeted biochemical experimentation. Eukaryotic ACAD subfamilies are traditionally considered as mitochondrial proteins, but we found evidence that in fungi one subfamily is located in peroxisomes and participates in a distinct β-oxidation pathway. Finally, we discern horizontal transfer, duplication, loss and secondary acquisition of ACAD genes during evolution of this family. Through these unorthodox expansion strategies, the ACAD family is proficient in utilizing a large range of fatty acids and amino acids—strategies that could have shaped the evolutionary history of many other ancient protein families.     

Post-Messinian evolutionary relationships across the Sicilian channel: Mitochondrial and nuclear markers link a new green toad from Sicily to African relatives

Background

The major birch pollen allergen, Bet v 1, is a member of the ubiquitous PR-10 family of plant pathogenesis-related proteins. In recent years, a number of diverse plant proteins with low sequence similarity to Bet v 1 was identified. In addition, determination of the Bet v 1 structure revealed the existence of a large superfamily of structurally related proteins. In this study, we aimed to identify and classify all Bet v 1-related structures from the Protein Data Bank and all Bet v 1-related sequences from the Uniprot database.

Results

Structural comparisons of representative members of already known protein families structurally related to Bet v 1 with all entries of the Protein Data Bank yielded 47 structures with non-identical sequences. They were classified into eleven families, five of which were newly identified and not included in the Structural Classification of Proteins database release 1.71. The taxonomic distribution of these families extracted from the Pfam protein family database showed that members of the polyketide cyclase family and the activator of Hsp90 ATPase homologue 1 family were distributed among all three superkingdoms, while members of some bacterial families were confined to a small number of species. Comparison of ligand binding activities of Bet v 1-like superfamily members revealed that their functions were related to binding and metabolism of large, hydrophobic compounds such as lipids, hormones, and antibiotics. Phylogenetic relationships within the Bet v 1 family, defined as the group of proteins with significant sequence similarity to Bet v 1, were determined by aligning 264 Bet v 1-related sequences. A distance-based phylogenetic tree yielded a classification into 11 subfamilies, nine exclusively containing plant sequences and two subfamilies of bacterial proteins. Plant sequences included the pathogenesis-related proteins 10, the major latex proteins/ripening-related proteins subfamily, and polyketide cyclase-like sequences.

Conclusion

The ubiquitous distribution of Bet v 1-related proteins among all superkingdoms suggests that a Bet v 1-like protein was already present in the last universal common ancestor. During evolution, this protein diversified into numerous families with low sequence similarity but with a common fold that succeeded as a versatile scaffold for binding of bulky ligands.     

Molecular phylogenetics of the family Cyprinidae (Actinopterygii: Cypriniformes) as evidenced by sequence variation in the first intron of S7 ribosomal protein-coding gene: further evidence from a nuclear gene of the systematic chaos in the family

The family Cyprinidae is the largest freshwater fish group in the world, including over 200 genera and 2100 species. The phylogenetic relationships of major clades within this family are simply poorly understood, largely because of the overwhelming diversity of the group; however, several investigators have advanced different hypotheses of relationships that pre- and post-date the use of shared-derived characters as advocated through phylogenetic systematics. As expected, most previous investigations used morphological characters. Recently, mitochondrial DNA (mtDNA) sequences and combined morphological and mtDNA investigations have been used to explore and advance our understanding of species relationships and test monophyletic groupings. Limitations of these studies include limited taxon sampling and a strict reliance upon maternally inherited mtDNA variation. The present study is the first endeavor to recover the phylogenetic relationships of the 12 previously recognized monophyletic subfamilies within the Cyprinidae using newly sequenced nuclear DNA (nDNA) for over 50 species representing members of the different previously hypothesized subfamily and family groupings within the Cyprinidae and from other cypriniform families as outgroup taxa. Hypothesized phylogenetic relationships are constructed using maximum parsimony and Basyesian analyses of 1042 sites, of which 971 sites were variable and 790 were phylogenetically informative. Using other appropriate cypriniform taxa of the families Catostomidae (Myxocyprinus asiaticus), Gyrinocheilidae (Gyrinocheilus aymonieri), and Balitoridae (Nemacheilus sp. and Beaufortia kweichowensis) as outgroups, the Cyprinidae is resolved as a monophyletic group. Within the family the genera Raiamas, Barilius, Danio, and Rasbora, representing many of the tropical cyprinids, represent basal members of the family. All other species can be classified into variably supported and resolved monophyletic lineages, depending upon analysis, that are consistent with or correspond to Barbini and Leuciscini. The Barbini includes taxa traditionally aligned with the subfamily Cyprininae sensu previous morphological revisionary studies by Howes (Barbinae, Labeoninae, Cyprininae and Schizothoracinae). The Leuciscini includes six other subfamilies that are mainly divided into three separate lineages. The relationships among genera and subfamilies are discussed as well as the possible origins of major lineages.     

Bioinformatic and enzymatic characterization of the MAPEG superfamily

The membrane associated proteins in eicosanoid and glutathione metabolism (MAPEG) superfamily includes structurally related membrane proteins with diverse functions of widespread origin. A total of 136 proteins belonging to the MAPEG superfamily were found in database and genome screenings. The members were found in prokaryotes and eukaryotes, but not in any archaeal organism. Multiple sequence alignments and calculations of evolutionary trees revealed a clear subdivision of the eukaryotic MAPEG members, corresponding to the six families of microsomal glutathione transferases (MGST) 1, 2 and 3, leukotriene C4 synthase (LTC4), 5-lipoxygenase activating protein (FLAP), and prostaglandin E synthase. Prokaryotes contain at least two distinct potential ancestral subfamilies, of which one is unique, whereas the other most closely resembles enzymes that belong to the MGST2/FLAP/LTC4 synthase families. The insect members are most similar to MGST1/prostaglandin E synthase. With the new data available, we observe that fish enzymes are present in all six families, showing an early origin for MAPEG family differentiation. Thus, the evolutionary origins and relationships of the MAPEG superfamily can be defined, including distinct sequence patterns characteristic for each of the subfamilies. We have further investigated and functionally characterized representative gene products from Escherichia coli, Synechocystis sp., Arabidopsis thaliana and Drosophila melanogaster, and the fish liver enzyme, purified from pike (Esox lucius). Protein overexpression and enzyme activity analysis demonstrated that all proteins catalyzed the conjugation of 1-chloro-2,4-dinitrobenzene with reduced glutathione. The E. coli protein displayed glutathione transferase activity of 0.11 micromol.min(-1).mg(-1) in the membrane fraction from bacteria overexpressing the protein. Partial purification of the Synechocystis sp. protein yielded an enzyme of the expected molecular mass and an N-terminal amino acid sequence that was at least 50% pure, with a specific activity towards 1-chloro-2,4-dinitrobenzene of 11 micromol.min(-1).mg(-1). Yeast microsomes expressing the Arabidopsis enzyme showed an activity of 0.02 micromol.min(-1).mg(-1), whereas the Drosophila enzyme expressed in E. coli was highly active at 3.6 micromol.min(-1).mg(-1). The purified pike enzyme is the most active MGST described so far with a specific activity of 285 micromol.min(-1).mg(-1). Drosophila and pike enzymes also displayed glutathione peroxidase activity towards cumene hydroperoxide (0.4 and 2.2 micromol.min(-1).mg(-1), respectively). Glutathione transferase activity can thus be regarded as a common denominator for a majority of MAPEG members throughout the kingdoms of life whereas glutathione peroxidase activity occurs in representatives from the MGST1, 2 and 3 and PGES subfamilies.     

Phylogenetic analysis of alpha-galactosidases of the GH27 family

Amino acid sequence analysis of alpha-galactosidases and other proteins of glycoside hydrolase family 27 (GH27) allowed isolation of three major subfamilies, 27a-27c. Unique isomalto-dextranase of Arthrobacter globiformis clustered separately. Eukaryotic proteins formed five clusters on a phylogenetic tree of the family. Bacterial GH27 proteins, which are relatively few, did not form stable clusters. A monophyletic origin of the GH27 family was demonstrated with the use of related proteins of the GH36 family. The structure of the active center and evolution of alpha-galactosidases are discussed.     

Medium-chain dehydrogenases/reductases (MDR). Family characterizations including genome comparisons and active site modeling.

Completed eukaryotic genomes were screened for medium-chain dehydrogenases/reductases (MDR). In the human genome, 23 MDR forms were found, a number that probably will increase, because the genome is not yet fully interpreted. Partial sequences already indicate that at least three further members exist. Within the MDR superfamily, at least eight families were distinguished. Three families are formed by dimeric alcohol dehydrogenases (ADH; originally detected in animals/plants), cinnamyl alcohol dehydrogenases (originally detected in plants) and tetrameric alcohol dehydrogenases (originally detected in yeast). Three further families are centred around forms initially detected as mitochondrial respiratory function proteins, acetyl-CoA reductases of fatty acid synthases, and leukotriene B4 dehydrogenases. The two remaining families with polyol dehydrogenases (originally detected as sorbitol dehydrogenase) and quinone reductases (originally detected as zeta-crystallin) are also distinct but with variable sequences. The most abundant families in the human genome are the dimeric ADH forms and the quinone oxidoreductases. The eukaryotic patterns are different from those of Escherichia coli. The different families were further evaluated by molecular modelling of their active sites as to geometry, hydrophobicity and volume of substrate-binding pockets. Finally, sequence patterns were derived that are diagnostic for the different families and can be used in genome annotations.     

Structure-Based Phylogenetic Analysis of the Lipocalin Superfamily

Lipocalins constitute a superfamily of extracellular proteins that are found in all three kingdoms of life. Although very divergent in their sequences and functions, they show remarkable similarity in 3-D structures. Lipocalins bind and transport small hydrophobic molecules. Earlier sequence-based phylogenetic studies of lipocalins highlighted that they have a long evolutionary history. However the molecular and structural basis of their functional diversity is not completely understood. The main objective of the present study is to understand functional diversity of the lipocalins using a structure-based phylogenetic approach. The present study with 39 protein domains from the lipocalin superfamily suggests that the clusters of lipocalins obtained by structure-based phylogeny correspond well with the functional diversity. The detailed analysis on each of the clusters and sub-clusters reveals that the 39 lipocalin domains cluster based on their mode of ligand binding though the clustering was performed on the basis of gross domain structure. The outliers in the phylogenetic tree are often from single member families. Also structure-based phylogenetic approach has provided pointers to assign putative function for the domains of unknown function in lipocalin family. The approach employed in the present study can be used in the future for the functional identification of new lipocalin proteins and may be extended to other protein families where members show poor sequence similarity but high structural similarity.     

Evolution of function in protein superfamilies, from a structural perspective

The recent growth in protein databases has revealed the functional diversity of many protein superfamilies. We have assessed the functional variation of homologous enzyme superfamilies containing two or more enzymes, as defined by the CATH protein structure classification, by way of the Enzyme Commission (EC) scheme. Combining sequence and structure information to identify relatives, the majority of superfamilies display variation in enzyme function, with 25 % of superfamilies in the PDB having members of different enzyme types. We determined the extent of functional similarity at different levels of sequence identity for 486,000 homologous pairs (enzyme/enzyme and enzyme/non-enzyme), with structural and sequence relatives included. For single and multi-domain proteins, variation in EC number is rare above 40 % sequence identity, and above 30 %, the first three digits may be predicted with an accuracy of at least 90 %. For more distantly related proteins sharing less than 30 % sequence identity, functional variation is significant, and below this threshold, structural data are essential for understanding the molecular basis of observed functional differences. To explore the mechanisms for generating functional diversity during evolution, we have studied in detail 31 diverse structural enzyme superfamilies for which structural data are available. A large number of variations and peculiarities are observed, at the atomic level through to gross structural rearrangements. Almost all superfamilies exhibit functional diversity generated by local sequence variation and domain shuffling. Commonly, substrate specificity is diverse across a superfamily, whilst the reaction chemistry is maintained. In many superfamilies, the position of catalytic residues may vary despite playing equivalent functional roles in related proteins. The implications of functional diversity within supefamilies for the structural genomics projects are discussed. More detailed information on these superfamilies is available at http://www.biochem.ucl.ac.uk/bsm/FAM-EC/.     

Phylogenetic diversification of the globin gene superfamily in chordates

Phylogenetic reconstructions provide a means of inferring the branching relationships among members of multigene families that have diversified via successive rounds of gene duplication and divergence. Such reconstructions can illuminate the pathways by which particular expression patterns and protein functions evolved. For example, phylogenetic analyses can reveal cases in which similar expression patterns or functional properties evolved independently in different lineages, either through convergence, parallelism, or evolutionary reversals. The purpose of this article is to provide a robust phylogenetic framework for interpreting experimental data and for generating hypotheses about the functional evolution of globin proteins in chordate animals. To do this, we present a consensus phylogeny of the chordate globin gene superfamily. We document the relative roles of gene duplication and whole-genome duplication in fueling the functional diversification of vertebrate globins, and we unravel patterns of shared ancestry among globin genes from representatives of the three chordate subphyla (Craniata, Urochordata, and Cephalochordata). Our results demonstrate the value of integrating phylogenetic analyses with genomic analyses of conserved synteny to infer the duplicative origins and evolutionary histories of globin genes. We also discuss a number of case studies that illustrate the importance of phylogenetic information when making inferences about the evolution of globin gene expression and protein function. Finally, we discuss why the globin gene superfamily presents special challenges for phylogenetic analysis, and we describe methodological approaches that can be used to meet those challenges.     

A new and unified nomenclature for male fertility restorer (RF) proteins in higher plants

The male fertility restorer (RF) proteins belong to extended protein families associated with the cytoplasmic male sterility in higher plants. Up till now, there is no devised nomenclature for naming the RF proteins. The systematic sequencing of new plant species in recent years has uncovered the existence of several novel RF genes and their encoded proteins. Their naming has been simply arbitrary and could not be adequately handled in the context of comparative functional genomics. We propose in this study a unified nomenclature for the RF extended protein families across all plant species. This new and unified nomenclature relies upon previously developed nomenclature for the first ever characterized RF gene, RF2A/ALDH2B2, a member of ALDH gene superfamily, and adheres to the guidelines issued by the ALDH Genome Nomenclature Committees. The proposed nomenclature reveals that RF gene superfamily encodes currently members of 51 families. This unified nomenclature accommodates functional RF genes and pseudogenes, and offers the flexibility needed to incorporate additional RFs as they become available in future. In addition, we provide a phylogenetic relationship between the RF extended families and use computational protein modeling to demonstrate the high divergence of RF functional specializations through specific structural features of selected members of RF superfamily.