Molecular evolution of nitrate reductase genes

To understand the evolutionary mechanisms and relationships of nitrate reductases (NRs), the nucleotide sequences encoding 19 nitrate reductase (NR) genes from 16 species of fungi, algae, and higher plants were analyzed. The NR genes examined show substantial sequence similarity, particularly within functional domains, and large variations in GC content at the third codon position and intron number. The intron positions were different between the fungi and plants, but conserved within these groups. The overall and nonsynonymous substitution rates among fungi, algae, and higher plants were estimated to be 4.33 × 10⁻¹⁰ and 3.29 × 10⁻¹⁰ substitutions per site per year. The three functional domains of NR genes evolved at about one-third of the rate of the N-terminal and the two hinge regions connecting the functional domains. Relative rate tests suggested that the nonsynonymous substitution rates were constant among different lineages, while the overall nucleotide substitution rates varied between some lineages. The phylogenetic trees based on NR genes correspond well with the phylogeny of the organisms determined from systematics and other molecular studies. Based on the nonsynonymous substitution rate, the divergence time of monocots and dicots was estimated to be about 340 Myr when the fungi–plant or algae–higher plant divergence times were used as reference points and 191 Myr when the rice–barley divergence time was used as a reference point. These two estimates are consistent with other estimates of divergence times based on these reference points. The lack of consistency between these two values appears to be due to the uncertainty of the reference times. Received: 10 April 1995 / Accepted: 10 September 1995     

Molecular evolution of calmodulin-like domain protein kinases (CDPKs) in plants and protists

Many genes for calmodulin-like domain protein kinases (CDPKs) have been identified in plants and Alveolate protists. To study the molecular evolution of the CDPK gene family, we performed a phylogenetic analysis of CDPK genomic sequences. Analysis of introns supports the phylogenetic analysis; CDPK genes with similar intron/exon structure are grouped together on the phylogenetic tree. Conserved introns support a monophyletic origin for plant CDPKs, CDPK-related kinases, and phosphoenolpyruvate carboxylase kinases. Plant CDPKs divide into two major branches. Plant CDPK genes on one branch share common intron positions with protist CDPK genes. The introns shared between protist and plant CDPKs presumably originated before the divergence of plants from Alveolates. Additionally, the calmodulin-like domains of protist CDPKs have intron positions in common with animal and fungal calmodulin genes. These results, together with the presence of a highly conserved phase zero intron located precisely at the beginning of the calmodulin-like domain, suggest that the ancestral CDPK gene could have originated from the fusion of protein kinase and calmodulin genes facilitated by recombination of ancient introns. Received: 11 July 2000 / Accepted: 18 April 2001     

Gene Conversion and Evolution of Daphniid Hemoglobins (Crustacea, Cladocera)

The extracellular hemoglobins of cladocerans derive from the aggregation of 12 two-domain globin subunits that are apparently encoded by four genes. This study establishes that at least some of these genes occur as a tandem array in both Daphnia magna and Daphnia exilis. The genes share a uniform structure; a bridge intron separates two globin domains which each include three exons and two introns. Introns are small, averaging just 77 bp, but a longer sequence (2.2–3.2 kb) separates adjacent globin genes. A survey of structural diversity in globin genes from other daphniids revealed three independent cases of intron loss, but exon lengths were identical, excepting a 3-bp insertion in exon 5 of Simocephalus. Heterogeneity in the extent of nucleotide divergence was marked among exons, largely as a result of the pronounced diversification of the terminal exon. This variation reflected, in part, varying exposure to concerted evolution. Conversion events were frequent in exons 1–4 but were absent from exons 5 and 6. Because of this difference, the results of phylogenetic analyses were strongly affected by the sequences employed in this construction. Phylogenies based on total nucleotide divergence in exons 1–4 revealed affinities among all genes isolated from a single species, reflecting the impact of gene conversion events. In contrast, phylogenies based on total nucleotide divergence in exons 5 and 6 revealed affinities among orthologous genes from different taxa. Received: 8 March 1999 / Accepted: 14 July 1999     

Molecular evolution of the Metazoan protein kinase C multigene family

Protein kinases C (PKCs) comprise closely related Ser/Thr kinases, ubiquitously present in animal tissues; they respond to second messengers, e.g., Ca²⁺ and/or diacylglycerol, to express their activities. Two PKCs have been sequenced from Geodia cydonium, a member of the lowest multicellular animals, the sponges (Porifera). One sponge G. cydonium PKC, GCPKC1, belongs to the ``novel' (Ca<sup>2+</sup>-independent) PKC (nPKC) subfamily while the second one, GCPKC2, has the hallmarks of the ``conventional' (Ca²⁺-dependent) PKC (cPKC) subfamily. The alignment of the Ser/Thr catalytic kinase domains, of the predicted aa sequences for these cDNAs with respective segments from previously reported sequences, revealed highest homology to PKCs from animals but also distant relationships to Ser/Thr kinases from protozoa, plants, and bacteria. However, a comparison of the complete structures of the sponge PKCs, which are—already—identical to those of nPKCs and cPKCs from higher metazoa, with the structures of protozoan, plant, and bacterial Ser/Thr kinases indicates that the metazoan PKCs have to be distinguished from the nonmetazoan enzymes. These data indicate that metazoan PKCs have a universal common ancestor which they share with the nonmetazoan Ser/Thr kinases with respect to the kinase domain, but they differ from them in overall structural composition. Received: 10 January 1996 / Accepted: 12 March 1996     

The Molecular Evolution of Visual Pigments of Freshwater Crayfishes (Decapoda: Cambaridae)

This study examines the diverse maximum wavelength absorption (λ_max) found in crayfishes (Decapoda: Cambaridae and Parastacidae) and the associated genetic variation in their opsin locus. We measured the wavelength absorption in the photoreceptors of six species that inhabit environments of different light intensities (i.e., burrows, streams, standing waters, and subterranean waters). Our results indicate that there is relatively little variation in λ_max (522–530 nm) among species from different genera and families. The existing variation did not correlate with the habitat differences of the crayfishes studied. We simultaneously sequenced the rhodopsin gene to identify the amino acid replacements that affect shifts in maximum wavelength absorption. We then related these to changes that correlated with shifts in λ_max by reconstructing ancestral character states using a maximum-likelihood approach. Using amino acid sequences obtained from five species (all were 301 amino acids in length), we identified a number of candidates for producing shifts of 4 to 8 nm in λ_max. These amino acid replacements occurred in similar regions to those involved in spectral shifts in vertebrates. Received: 12 March 1997 / Accepted: 3 June 1997     

Evolution of trophic types in emperor fishes ( Lethrinus,Lethrinidae, Percoidei) based on cytochrome B gene sequence variation

Three trophic categories exist within emperor fishes, genus Lethrinus, relating to body form and dentition type. One group contains low-bodied, high speed, stalking predators with conical teeth. Another group comprises high-bodied, slow speed carnivores with molariform teeth capable of crushing hard-shelled benthic prey. A third group is also high-bodied but with conical teeth feeding mostly on small or soft-shelled benthic prey. Inferring the evolution of these trophic types within Lethrinus using morphology is problematic since these characters are typically correlated with feeding mode and are potentially homoplasious. We use mitochondrial DNA sequences, to independently determine a phylogenetic hypothesis for Lethrinus, which are not dependent on morphological characters relating to trophic categories. We analyzed complete cytochrome b gene sequences (1140 bp) for 20 species of Lethrinus, representing the three trophic types, and for 13 outgroup species, including four other representatives of the Lethrinidae. A monophyletic Lethrinidae did not resolve, but the monophyly of Lethrinus is well supported. In addition, two major clades within Lethrinus are well supported. One of these clades exclusively contains low-bodied species with conical teeth while the other clade only comprises the high-bodied species with molariform teeth. A high-bodied species with conical teeth, Lethrinus miniatus, appears most ancestral and sister to all other Lethrinus species. We hypothesize that this generalist trophic type was the evolutionary precursor to both of the other primary trophic types.     

The tryptophan biosynthetic pathway of aphid endosymbionts (Buchnera): Genetics and evolution of plasmid-associated anthranilate synthase (trpEG) within the aphididae

The bacterial endosymbionts (Buchnera) from the aphids Rhopalosiphum padi, R. maidis, Schizaphis graminum, and Acyrthosiphon pisum contain the genes for anthranilate synthase (trpEG) on plasmids made up of one or more 3.6-kb units. Anthranilate synthase is the first as well as the rate-limiting enzyme in the tryptophan biosynthetic pathway. The amplification of trpEG on plasmids may result in an increase of enzyme protein and overproduction of this essential amino acid, which is required by the aphid host. The nucleotide sequence of trpEG from endosymbionts of different species of aphids is highly conserved, as is an approximately 500-bp upstream DNA segment which has the characteristics of an origin of replication. Phylogenetic analyses were performed using trpE and trpG from the endosymbionts of these four aphids as well as from the endosymbiont of Schlechtendalia chinensis, in which trpEG occurs on the chromosome. The resulting phylogeny was congruent with trees derived from sequences of two chromosome-located bacterial genes (part of trpB and 16S ribosomal DNA). In turn, trees obtained from plasmid-borne and bacterial chromosome-borne sequences were congruent with the tree resulting from phylogenetic analysis of three aphid mitochondrial regions (portions of the small and large ribosomal DNA subunits, as well as cytochrome oxidase II). Congruence of trees based on genes from host mitochondria and from bacteria adds to previous support for exclusively vertical transmission of the endosymbionts within aphid lineages. Congruence with trees based on plasmid-borne genes supports the origin of the plasmid-borne trpEG from the chromosomal genes of the same lineage and the absence of subsequent plasmid exchange among endosymbionts of different species of aphids. Received: 22 August 1995 / Accepted: 6 September 1995     

Rapid evolution of a young L1 (LINE-1) clade in recently speciated rattus taxa

L1 elements are retrotransposons that have been replicating and evolving in mammalian genomes since before the mammalian radiation. Rattus norvegicus shares the young L1_mlvi2 clade only with its sister taxon, Rattus cf moluccarius. Here we compared the L1_mlvi2 clade in these recently diverged species and found that it evolved rapidly into closely related but distinct clades: the L1_mlvi2-rm clade (or subfamily), characterized here from R. cf moluccarius, and the L1_mlvi2-rn clade, originally described in R. norvegicus. In addition to other differences, these clades are distinguished by a cluster of amino acid replacement substitutions in ORF I. Both rat species contain the L1_mlvi2-rm clade, but the L1_mlvi2-rn clade is restricted to R. norvegicus. Therefore, the L1_mlvi2-rm clade arose prior to the divergence of R. norvegicus and R. cf moluccarius, and the L1_mlvi2-rn clade amplified after their divergence. The total number of L1_mlvi2-rm elements in R. cf moluccarius is about the same as the sum of the L1_mlvi2-rm and L1_mlvi2-rn elements in R. norvegicus. The possibility that L1 amplification is in some way limited so that the two clades compete for replicative supremacy as well as the implications of the other distinguishing characteristic of the L1_mlvi2-rn and L1_mlvi2-rm clades are discussed. Received: 6 November 1996 / Accepted: 12 April 1997     

Genetic Algorithm-Based Maximum-Likelihood Analysis for Molecular Phylogeny

A heuristic approach to search for the maximum-likelihood (ML) phylogenetic tree based on a genetic algorithm (GA) has been developed. It outputs the best tree as well as multiple alternative trees that are not significantly worse than the best one on the basis of the likelihood criterion. These near-optimum trees are subjected to further statistical tests. This approach enables ones to infer phylogenetic trees of over 20 taxa taking account of the rate heterogeneity among sites on practical time scales on a PC cluster. Computer simulations were conducted to compare the efficiency of the present approach with that of several likelihood-based methods and distance-based methods, using amino acid sequence data of relatively large (5–24) taxa. The superiority of the ML method over distance-based methods increases as the condition of simulations becomes more realistic (an incorrect model is assumed or many taxa are involved). This approach was applied to the inference of the universal tree based on the concatenated amino acid sequences of vertically descendent genes that are shared among all genomes whose complete sequences have been reported. The inferred tree strongly supports that Archaea is paraphyletic and Eukarya is specifically related to Crenarchaeota. Apart from the paraphyly of Archaea and some minor disagreements, the universal tree based on these genes is largely consistent with the universal tree based on SSU rRNA. Received: 4 January 2001 / Accepted: 16 May 2001     

Muscle metabolism and growth in Antarctic fishes (suborder Notothenioidei): evolution in a cold environment

The radiation of notothenioid fishes (order Perciformes) in the Southern Ocean provides a model system for investigating evolution and adaptation to a low temperature environment. The Notothenioid fishes comprising eight families, 43 genera and 122 species dominate the fish fauna in Antarctica. The diversification of the clade probably began 15–20 million years ago after the formation of the Antarctic Polar Front. The radiation was, therefore, associated with climatic cooling down to the present day temperature of −1.86 °C. Origins and Evolution of the Antarctic Biota Geological Society Special Publication No. 47, Geological Society of London. pp. 253–268). The success of the group has been closely linked with the evolution of glycopeptide and peptide antifreezes, which are amongst the most abundant proteins in blood and interstitial fluid. The radiation of the clade has been associated with disaptation (evolutionary loss of function) and recovery. For example, it is thought that the icefishes (Channichyidae) lost haemoglobin through a single mutational event leading to the deletion of the entire β-globin gene and the 5′ end of the linked α-globin gene, resulting in compensatory adaptations of the cardiovascular system. Phylogenetically based statistical methods also indicate a progressive and dramatic reduction in the number of skeletal muscle fibres (FN_max) at the end of the recruitment phase of growth in basal compared to derived families. The reduction in FN_max is associated with a compensatory increase in the maximum fibre diameter, which can reach 100 μm in slow and 600 μm in fast muscle fibres. At −1 to 0 °C, the oxygen consumption of isolated mitochondria per mg mitochondrial protein shows no evidence of up-regulation relative to mitochondria from temperate and tropical Perciform fishes. The mitochondria content of slow muscle fibres in Antarctic notothenioids is towards the upper end of the range reported for teleosts with similar lifestyles, reaching 50% in Channichthyids. High mitochondrial densities facilitate ATP production and oxygen diffusion through the membrane lipid compartment of the fibre. Modelling studies suggest that adequate oxygen flux in the large diameter muscle fibres of notothenioids is possible because of the reduced metabolic demand and enhanced solubility of oxygen associated with low temperature. At the whole animal level size-corrected resting metabolic rate fits on the same temperature relationship as for Perciformes from warmer climates. It seems likely that the additional energetic costs associated with antifreeze synthesis and high mitochondrial densities are compensated for by reductions in other energy requiring processes: a hypothesis that could be tested with detailed energy budget studies. One plausible candidate is a reduction in membrane leak pathways linked to the loss of muscle fibres, which would serve to minimise the cost of maintaining ionic gradients.     

Molecular Data from the 16S rRNA Gene for the Phylogeny of Pectinidae (Mollusca: Bivalvia)

The phylogenetic relationships among the species belonging to the family Pectinidae are still an issue of debate. The mitochondrial DNA sequences from the large ribosomal RNA gene may be of great value for systematic and phylogenetic studies within families. Partial sequences of the 16S rRNA gene were obtained for the scallop species Adamussium colbecki, Aequipecten opercularis, Chlamys glabra, C. islandica, C. varia, and Pecten jacobeus and compared with the published sequence of Pecten maximus. The present molecular data show that Chlamys are polyphyletic and do not support the assignment of these species to the two subfamilies Chlamydinae and Pectininae. Moreover, the minimal genetic distance between P. maximus and P. jacobeus suggests that they could belong to the same species. Received: 24 May 1999 / Accepted: 1 September 1999     

Molecular evolution of a portion of the mitochondrial 16S ribosomal gene region in scleractinian corals

Relationships among families and suborders of scleractinian corals are poorly understood because of difficulties 1) in making inferences about the evolution of the morphological characters used in coral taxonomy and 2) in interpreting their 240-million-year fossil record. Here we describe patterns of molecular evolution in a segment of the mitochondrial (mt) 16S ribosomal gene from taxa of 14 families of corals and the use of this gene segment in a phylogenetic analysis of relationships within the order. We show that sequences obtained from scleractinians are homologous to other metazoan 16S ribosomal sequences and fall into two distinct clades defined by size of the amplified gene product. Comparisons of sequences from the two clades demonstrate that both sets of sequences are evolving under similar evolutionary constraints: they do not differ in nucleotide composition, numbers of transition and transversion substitutions, spatial patterns of substitutions, or in rates of divergence. The characteristics and patterns observed in these sequences as well as the secondary structures, are similar to those observed in mt 16S ribosomal DNA sequences from other taxa. Phylogenetic analysis of these sequences shows that they are useful for evaluating relationships within the order. The hypothesis generated from this analysis differs from traditional hypotheses for evolutionary relationships among the Scleractinia and suggests that a reevaluation of evolutionary affinities in the order is needed. Received: 4 September 1996 / Accepted: 7 April 1997     

Molecular dissection of the large mechanosensitive ion channel (MscL) of E. coli: Mutants with altered channel gating and pressure sensitivity

In the search for the essential functional domains of the large mechanosensitive ion channel (MscL) of E. coli, we have cloned several mutants of the mscL gene into a glutathione S-transferase fusion protein expression system. The resulting mutated MscL proteins had either amino acid additions, substitutions or deletions in the amphipathic N-terminal region, and/or deletions in the amphipathic central or hydrophilic C-terminal regions. Proteolytic digestion of the isolated fusion proteins by thrombin yielded virtually pure recombinant MscL proteins that were reconstituted into artificial liposomes and examined for function by the patch-clamp technique. The addition of amino acid residues to the N-terminus of the MscL did not affect channel activity, whereas N-terminal deletions or changes to the N-terminal amino acid sequence were poorly tolerated and resulted in channels exhibiting altered pressure sensitivity and gating. Deletion of 27 amino acids from the C-terminus resulted in MscL protein that formed channels similar to the wild-type, while deletion of 33 C-terminal amino acids extinguished channel activity. Similarly, deletion of the internal amphipathic region of the MscL abolished activity. In accordance with a recently proposed spatial model of the MscL, our results suggest that (i) the N-terminal portion participates in the channel activation by pressure, and (ii) the essential channel functions are associated with both, the putative central amphipathic α-helical portion of the protein and the six C-terminal residues RKKEEP forming a charge cluster following the putative M2 membrane spanning α-helix. Received: 25 September 1996/Revised: 21 November 1996     

Molecular Phylogeny of Metazoan Intermediate Filament Proteins

We have cloned cytoplasmic intermediate filament (IF) proteins from a large number of invertebrate phyla using cDNA probes, the monoclonal antibody IFA, peptide sequence information, and various RT-PCR procedures. Novel IF protein sequences reported here include the urochordata and nine protostomic phyla, i.e., Annelida, Brachiopoda, Chaetognatha, Echiura, Nematomorpha, Nemertea, Platyhelminthes, Phoronida, and Sipuncula. Taken together with the wealth of data on IF proteins of vertebrates and the results on IF proteins of Cephalochordata, Mollusca, Annelida, and Nematoda, two IF prototypes emerge. The L-type, which includes 35 sequences from 11 protostomic phyla, shares with the nuclear lamins the long version of the coil 1b subdomain and, in most cases, a homology segment of some 120 residues in the carboxyterminal tail domain. The S-type, which includes all four subfamilies (types I to IV) of vertebrate IF proteins, lacks 42 residues in the coil 1b subdomain and the carboxyterminal lamin homology segment. Since IF proteins from all three phyla of the chordates have the 42-residue deletion, this deletion arose in a progenitor prior to the divergence of the chordates into the urochordate, cephalochordate, and vertebrate lineages, possibly already at the origin of the deuterostomic branch. Four phyla recently placed into the protostomia on grounds of their 18S rDNA sequences (Brachiopoda, Nemertea, Phoronida, and Platyhelminthes) show IF proteins of the L-type and fit by sequence identity criteria into the lophotrochozoic branch of the protostomia. Received: 2 April 1998 / Accepted: 19 June 1998     

Origin and evolution of GBV-C/hepatitis G virus and relationships with ancient human migrations

The GB virus C/hepatitis G virus (GBV-C/HGV) is a newly identified human RNA virus, belonging to the Flaviviridae family. Persistent infection by GBV-C/HGV is common in humans, and genetically divergent isolates have been identified in different parts of the world. Due to the absence of a real pathogenic role of GBV-C/HGV in liver disease and its extremely low mutation rate, this virus is a potential marker to trace prehistoric links between human populations. In this study, origin and evolution of GBV-C/HGV were examined using a set of fully sequenced strains of worldwide origin. A first phylogenetic analysis, addressed to the short (255 nucleotides) NS5A overlapping coding region by the neighbor-joining method, suggested an ancient African origin of GBV-C/HGV. This notion was confirmed when the same analysis was applied to the genomic regions showing the lowest rate of synonymous substitutions, covering one-fourth (2184 nucleotides) of the total coding potential of the virus genome. By using a multivariate statistical method and extending the analysis to the complete coding region, fine details of the evolutionary history of GBV-C/HGV were further elucidated. By this approach, isolates from Southeast Asia appeared to be the most closely related to those of African origin, consistent with a major route of ancient human migrations from Africa to southeastern parts of the Asian continent. Received: 26 October 2000 / Accepted: 28 February 2001     

Molecular Regulation of Renal Phosphate Transport

Received: 18 April 1996/Revised: 26 June 1996     

Simultaneous Molecular and Morphological Analysis of Braconid Relationships (Insecta: Hymenoptera: Braconidae) Indicates Independent mt-tRNA Gene Inversions Within a Single Wasp Family

We investigated the phylogeny of the Braconidae (Insecta: Hymenoptera) with a much expanded data set compared with that of previous attempts, employing 16S and 28S rDNA gene fragments, together with a suite of morphological characters, from 74 ingroup taxa. Most notably, parsimony analyses under a range of models recovered the Aphidiinae as sister group to the cyclostomes and the Ichneutinae as sister group to the microgastroids. The cyclostomes were recovered as a natural group only if certain, putatively misplaced genera (Mesostoa, Aspilodemon) were excluded from them. Further, mapping of rearrangement characters onto this phylogeny of the Braconidae indicated parallel inversions of the mt-tRNA^D gene, with the two instances of inversion distinguishable by the presence or absence of an additional tRNA gene (tRNA^H). This is the first report of a parallel inversion of a mt-tRNA gene and makes the Braconidae the first metazoan family to display both parallel inversions and translocations. Received: 6 April 2001 / Accepted: 9 July 2001     

Multiple Lineages of the Mitochondrial Gene NADH Dehydrogenase Subunit 1 (ND1) in Parasitic Helminths: Implications for Molecular Evolutionary Studies of Facultatively Anaerobic Eukaryotes

The task of using partial ND1 sequences to infer a phylogeny for species of the genus Paragonimus (Trematoda: Digenea) was complicated by the discovery of at least two ND1 lineages within individual worms. The divergence of the ND1 lineages is shown by phylogenetic analysis not only to predate the divergence of the three Paragonimus species or species groups investigated but also the divergence of some trematode families. Some sequences are clearly pseudogenes as they contain single base deletions and/or premature termination codons. The presence of both pseudogenes and/or mitochondrial heteroplasmy are invoked to explain the presence of multiple and divergent ND1 lineages in these trematodes, which have two distinct cytochemical types of mitochondria. The implications for phylogenetic studies generally and of parasitic helminths specifically, using ND1 sequence data, are discussed. The ability of these organisms to adapt their metabolic processes to the variable availability of oxygen as an electron acceptor are proposed to explain some of the molecular diversity observed in parasitic helminths and possibly also in other anaerobically adapted eukaryotes. Received: 18 October 1999 / Accepted: 23 June 2000     

Na+:HCO(3-) cotransporters (NBC): cloning and characterization

The sodium bicarbonate cotransporter (NBC1) is essential for bicarbonate transport across plasma membranes in epithelial and nonepithelial cells. The direction of the NaHCO₃ movement in secretory epithelia is opposite to that in reabsorptive epithelia. In secretory epithelia (such as pancreatic duct cells) NBC is responsible for the transport of bicarbonate from blood to the cell for eventual secretion at the apical membrane. In reabsorptive epithelia (such as kidney proximal tubule cells) NBC is responsible for the reabsorption of bicarbonate from cell to the blood. In nonepithelial cells this transporter is mainly involved with cell pH regulation. Recent molecular cloning experiments have identified the existence of four NBC isoforms (NBC1, 2, 3 and 4) and two NBC-related proteins AE4 and NCBE (Anion Exchanger 4 and Na-dependent Chloride-Bicarbonate Exchanger). All but AE4 are presumed to mediate the cotransport of Na⁺ and HCO₃ ⁻ under normal conditions and may be functionally altered in certain pathologic states. NBC1 shows a limited tissue expression pattern, is electrogenic and plays an important role in bicarbonate reabsorption in kidney proximal tubule. In addition to the kidney, NBC1 is expressed in pancreatic duct cells, is activated by cystic fibrosis transmembrane conductance regulator (CFTR) and plays an important role in HCO₃ ⁻ secretion. NBC2 and NBC3 have a wider tissue distribution than NBC1, are electroneutral, and are involved with cell pH regulation. The characterization of NBC4 is incomplete. The NBC-related protein called NCBE mediates Na-dependent, Cl⁻/Bicarbonate Exchange. The purpose of this review is to summarize recent advances on the cloning of NBC isoforms and related proteins and their role and regulation in physiologic and pathologic states. Received: 26 February 2001/Revised: 14 May 2001