A comprehensive evolutionary analysis of aquaporins, a family of intrinsic membrane proteins that function as water channels, was conducted to establish groups of homology (i.e., to identify orthologues and paralogues) within the family and to gain insights into the functional constraints acting on the structure of the aquaporin molecule structure. Aquaporins are present in all living organisms, and therefore, they provide an excellent opportunity to further our understanding of the broader biological significance of molecular evolution by gene duplication followed by functional and structural specialization. Based on the resulting phylogeny, the 153 channel proteins analyzed were classified into six major paralogous groups: (1) GLPs, or glycerol-transporting channel proteins, which include mammalian AQP3, AQP7, and AQP9, several nematode paralogues, a yeast paralogue, and Escherichia coli GLP; (2) AQPs, or aquaporins, which include metazoan AQP0, AQP1, AQP2, AQP4, AQP5, and AQP6; (3) PIPs, or plasma membrane intrinsic proteins of plants, which include PIP1 and PIP2; (4) TIPs, or tonoplast intrinsic proteins of plants, which include alphaTIP, gammaTIP, and deltaTIP; (5) NODs, or nodulins of plants; and (6) AQP8s, or metazoan aquaporin 8 proteins. Of these groups, AQPs, PIPs, and TIPs cluster together. According to the results, the capacity to transport glycerol shown by several members of the family was acquired only early in the history of the family. The new phylogeny reveals that several water channel proteins are misclassified and require reassignment, whereas several previously undetermined ones can now be classified with confidence. The deduced phylogenetic framework was used to characterize the molecular features of water channel proteins. Three motifs are common to all family members: AEF (Ala-Glu-Phe), which is located in the N-terminal domain; and two NPA (Asp-Pro-Ala) boxes, which are located in the center and C-terminal domains, respectively. Other residues are found to be conserved within the major groups but not among them. Overall, the PIP subfamily showed the least variation. In general, no radical amino acid replacements affecting tertiary structure were identified, with the exception of Ala-->Ser in the TIP subfamily. Constancy of rates of evolution was demonstrated within the different paralogues but rejected among several of them (GLP and NOD).