Evolution of the secretoglobins: a genomic and proteomic view

Hello, kitty: could cat allergy be a form of intoxication?

Background

The relationship between slow loris (Nycticebus spp.) venom (BGE protein) and the major cat allergen (Fel d 1) from domestic cat (Felis catus) is known for about two decades. Along this time, evidence was accumulated regarding convergences between them, including their almost identical mode of action.

Methods

Large-scale database mining for Fel d 1 and BGE proteins in Felidae and Nycticebus spp., alignment, phylogeny proposition and molecular modelling, associated with directed literature review were assessed.

Results

Fel d 1 sequences for 28 non-domestic felids were identified, along with two additional loris BGE protein sequences. Dimer interfaces are less conserved among sequences, and the chain 1 shows more sequence similarity than chain 2. Post-translational modification similarities are highly probable.

Conclusions

Fel d 1 functions beyond allergy are discussed, considering the great conservation of felid orthologs of this protein. Reasons for toxicity being found only in domestic cats are proposed in the context of domestication. The combination of the literature review, genome-derived sequence data, and comparisons with the venomous primate slow loris may point to domestic cats as potentially poisonous mammals.

Another cat and mouse game: Deciphering the evolution of the SCGB superfamily and exploring the molecular similarity of major cat allergen Fel d 1 and mouse ABP using computational approaches

The mammalian secretoglobin (SCGB) superfamily contains functionally diverse members, among which the major cat allergen Fel d 1 and mouse salivary androgen-binding protein (ABP) display similar subunits. We searched for molecular similarities between Fel d 1 and ABP to examine the possibility that they play similar roles. We aimed to i) cluster the evolutionary relationships of the SCGB superfamily; ii) identify divergence patterns, structural overlap, and protein-protein docking between Fel d 1 and ABP dimers; and iii) explore the residual interaction between ABP dimers and steroid binding in chemical communication using computational approaches. We also report that the evolutionary tree of the SCGB superfamily comprises seven unique palm-like clusters, showing the evolutionary pattern and divergence time tree of Fel d 1 with 28 ABP paralogs. Three ABP subunits (A27, BG27, and BG26) share phylogenetic relationships with Fel d 1 chains. The Fel d 1 and ABP subunits show similarities in terms of sequence conservation, identical motifs and binding site clefts. Topologically equivalent positions were visualized through superimposition of ABP A27:BG27 (AB) and ABP A27:BG26 (AG) dimers on a heterodimeric Fel d 1 model. In docking, Fel d 1-ABP dimers exhibit the maximum surface binding ability of AG compared with that of AB dimers and the several polar interactions between ABP dimers with steroids. Hence, cat Fel d 1 is an ABP-like molecule in which monomeric chains 1 and 2 are the equivalent of the ABPA and ABPBG monomers, respectively. These findings suggest that the biological and molecular function of Fel d 1 is similar to that of ABP in chemical communication, possibly via pheromone and/or steroid binding.

Studies of an Androgen-Binding Protein Knockout Corroborate a Role for Salivary ABP in Mouse Communication

The house mouse Androgen-binding protein (Abp) gene family is comprised of 64 paralogs, 30 Abpa and 34 Abpbg, encoding the alpha (ABPA) and beta-gamma (ABPBG) protein subunits that are disulfide-bridged to form dimers in secretions. Only 14 Abp genes are expressed in distinct patterns in the lacrimal (11) and submandibular glands (3). We created a knockout mouse line lacking two of the three genes expressed in submandibular glands, Abpa27 and Abpbg27, by replacing them with the neomycin resistance gene. The knockout genotype (-/-) showed no Abpa27 or Abpbg27 transcripts in submandibular gland complementary DNA (cDNA) libraries and there was a concomitant lack of protein expression of ABPA27 and ABPBG27 in the -/- genotype saliva, shown by elimination of these two proteins from the saliva proteome and the loss of cross-reactive material in the acinar cells of the submandibular glands. We also observed a decrease in BG26 protein in the -/- animals, suggesting monomer instability. Overall, we observed no major phenotypic changes in the -/- genotype, compared with their +/+ and +/- siblings raised in a laboratory setting, including normal growth curves, tissue histology, fecundity, and longevity. The only difference is that male and female C57BL/6 mice preferred saliva of the opposite sex containing ABP statistically significantly more than saliva of the opposite sex without ABP in a Y-maze test. These results show for the first time that mice can sense the presence of ABP between saliva targets with and without ABPs, and that they spend more time investigating the target containing ABP.

Selection shaped the evolution of mouse androgen-binding protein (ABP) function and promoted the duplication of Abp genes

In the present article, we summarize two aspects of our work on mouse ABP (androgen-binding protein): (i) the sexual selection function producing incipient reinforcement on the European house mouse hybrid zone, and (ii) the mechanism behind the dramatic expansion of the Abp gene region in the mouse genome. Selection unifies these two components, although the ways in which selection has acted differ. At the functional level, strong positive selection has acted on key sites on the surface of one face of the ABP dimer, possibly to influence binding to a receptor. A different kind of selection has apparently driven the recent and rapid expansion of the gene region, probably by increasing the amount of Abp transcript, in one or both of two ways. We have shown previously that groups of Abp genes behave as LCRs (low-copy repeats), duplicating as relatively large blocks of genes by NAHR (non-allelic homologous recombination). The second type of selection involves the close link between the accumulation of L1 elements and the expansion of the Abp gene family by NAHR. It is probably predicated on an initial selection for increased transcription of existing Abp genes and/or an increase in Abp gene number providing more transcriptional sites. Either or both could increase initial transcript production, a quantitative change similar to increasing the volume of a radio transmission. In closing, we also provide a note on Abp gene nomenclature.

Did androgen-binding protein paralogs undergo neo- and/or Subfunctionalization as the Abp gene region expanded in the mouse genome?

The Androgen-binding protein (Abp) region of the mouse genome contains 30 Abpa genes encoding alpha subunits and 34 Abpbg genes encoding betagamma subunits, their products forming dimers composed of an alpha and a betagamma subunit. We endeavored to determine how many Abp genes are expressed as proteins in tears and saliva, and as transcripts in the exocrine glands producing them. Using standard PCR, we amplified Abp transcripts from cDNA libraries of C57BL/6 mice and found fifteen Abp gene transcripts in the lacrimal gland and five in the submandibular gland. Proteomic analyses identified proteins corresponding to eleven of the lacrimal gland transcripts, all of them different from the three salivary ABPs reported previously. Our qPCR results showed that five of the six transcripts that lacked corresponding proteins are expressed at very low levels compared to those transcripts with proteins. We found 1) no overlap in the repertoires of expressed Abp paralogs in lacrimal gland/tears and salivary glands/saliva; 2) substantial sex-limited expression of lacrimal gland/tear expressed-paralogs in males but no sex-limited expression in females; and 3) that the lacrimal gland/tear expressed-paralogs are found exclusively in ancestral clades 1, 2 and 3 of the five clades described previously while the salivary glands/saliva expressed-paralogs are found only in clade 5. The number of instances of extremely low levels of transcription without corresponding protein production in paralogs specific to tears and saliva suggested the role of subfunctionalization, a derived condition wherein genes that may have been expressed highly in both glands ancestrally were down-regulated subsequent to duplication. Thus, evidence for subfunctionalization can be seen in our data and we argue that the partitioning of paralog expression between lacrimal and salivary glands that we report here occurred as the result of adaptive evolution.

Evolution of the ABPA subunit of androgen-binding protein expressed in the submaxillary glands in New and Old World rodent taxa

The salivary androgen-binding proteins (ABPs) are members of the secretoglobin gene family present in mammals. Each ABP is a heterodimer assembled as an ABPA subunit encoded by an Abpa gene and linked by disulfide bridges to an ABPBG subunit encoded by an Abpbg gene. The ABP dimers are secreted into the saliva of mice and then transferred to the pelage after grooming and subsequently to the environment allowing an animal to mark territory with a biochemical signal. The putative role of the mouse salivary ABPs is that of pheromones mediating mate selection resulting in assortative mating in the Mus musculus species complex. We focused on comparing patterns of molecular evolution between the Abpa genes expressed in the submaxillary glands of species of New World and Old World muroids. We found that in both sets of rodents the Abpa genes expressed in the submaxillary glands appear to be evolving under a similar evolutionary regime, with relatively high nonsynonymous substitution rates, suggesting that ABP might play a similar biological role in both systems. Thus, ABP could be involved with mate recognition and species isolation in New World as well as Old World muroids.

Positive selection shaped the convergent evolution of independently expanded kallikrein subfamilies expressed in mouse and rat saliva proteomes

We performed proteomics studies of salivas from the genome mouse (C57BL/6 strain) and the genome rat (BN/SsNHsd/Mcwi strain). Our goal was to identify salivary proteins with one or more of three characteristics that may indicate that they have been involved in adaptation: 1) rapid expansion of their gene families; 2) footprints of positive selection; and/or 3) sex-limited expression. The results of our proteomics studies allow direct comparison of the proteins expressed and their levels between the sexes of the two rodent species. Twelve members of the Mus musculus species-specific kallikrein subfamily Klk1b showed sex-limited expression in the mouse saliva proteomes. By contrast, we did not find any of the Rattus norvegicus species-specific kallikrein subfamily Klk1c proteins in male or female genome rat, nor transcripts in their submandibular glands. On the other hand, we detected expression of this family as transcripts in the submandibular glands of both sexes of Sprague-Dawley rats. Using the CODEML program in the PAML package, we demonstrate that the two rodent kallikrein subfamilies have apparently evolved rapidly under the influence of positive selection that continually remodeled the amino acid sites on the same face in the members of the subfamilies. Thus, although their kallikrein subfamily expansions were independent, this evolutionary pattern has occurred in parallel in the two rodent species, suggesting a form of convergent evolution at the molecular level. On the basis of this new data, we suggest that the previous speculative function of the species-specific rodent kallikreins as important solely in wound healing in males be investigated further. In addition to or instead of that function, we propose that their sex-limited expression, coupled with their rapid evolution may be clues to an as-yet-undetermined interaction between the sexes.

High abundance of testosterone and salivary androgen-binding protein in the lateral nasal gland of male mice

To better understand androgen function in the mammalian nose, we have determined the levels of testosterone (T) in the olfactory mucosa (OM, which harbors the olfactory receptor neurons) and the lateral nasal gland (LNG, which is the largest anterior nasal gland) of C57BL/6 mice. The results indicated that, in adult male mice, T levels in the LNG were substantially higher than those in the OM and other non-reproductive or non-endocrine tissues examined, including liver, kidney, and brain. Furthermore, in the LNG, the high T levels were accompanied by high levels of salivary androgen-binding protein (sABP) and low microsomal T-hydroxylase activities. The high abundance of T and sABP in the LNG suggests not only that the LNG is a storage site for androgen, but also the possibility that unusually high T levels may occur in other organs that have abundant expression of sABP but low expression of steroid-metabolizing enzymes. Our findings suggest a critical need to determine androgen levels in various organs, as well as to establish the functional significance of an unusually high T level in the LNG, a gland known for its secretion of biologically active molecules, such as odorant binding proteins and immunoglobulin A, to the nasal cavity.

Signalling components of the house mouse mate recognition system

Subspecies-specific mate recognition may represent significant barrier to gene flow between diverged genomes potentially leading to speciation. In the house mouse, assortative mating involves the coevolution of several signals and receptors. We compared signalling ability of bedding material, faeces, urine, saliva, salivary androgen binding proteins (ABP) and combinations of urine with saliva and urine with ABP in mate choice in two wild-derived inbred strains (one of Mus musculus musculus and one of Mus musculus domesticus origin). We observed high levels of variation in assortative preferences between the two strains and sexes. The strongest preferences were observed in M. m. musculus-derived individuals in tests where urine was present either alone or as part of a composite signal target. M. m. domesticus-derived mice displayed strain-specific preferences for faeces. Saliva was the least preferred stimulus in both strains and sexes. No effect of two-compound cues was detected. We conclude that there is divergence across both the stimulus and preference parts of the recognition system for both house mouse strains. Of the tested stimuli, those that have the capacity to carry a signal for extended periods under natural conditions (such as urine and faeces) seem to be the most important substances in strain-specific recognition.

The chicken egg white proteome

Using 1-D PAGE and LC-MS/MS and MS(3) we identified 78 chicken egg white proteins, 54 of which were identified in egg white for the first time. All proteins were quantitated by calculating their exponentially modified protein abundance index (emPAI). Some previously known egg white components not characterized by amino acid sequences before, such as alpha-2-macroglobulin, were associated to a sequence for the first time. The predicted sequence was confirmed by MS-sequenced peptides covering 42% of the entire sequence. alpha-2-Macroglobulin occurred in egg white at the same concentration as ovostatin with which it showed 35% identity. For other proteins, which were previously only characterized by partial sequences, such as beta-ovomucin or ovalbumin X, we identified and confirmed predicted complete sequences with a high coverage by MS-sequenced peptides. New proteins included a 7 kDa protein consisting of a single secretoglobin sequence (ovosecretoglobin), a 7 kDa protein with similarity to black swan cygnin and turkey meleagrin (gallin) and proteins involved in binding, modification, and possibly detoxification, of bacterial lipopolysaccaride. The list of egg white proteins provided is by far the most comprehensive at present and is intended to serve as a starting point for the isolation and functional characterization of interesting new proteins.

Rapid bursts of androgen-binding protein (Abp) gene duplication occurred independently in diverse mammals

BACKGROUND

The draft mouse (Mus musculus) genome sequence revealed an unexpected proliferation of gene duplicates encoding a family of secretoglobin proteins including the androgen-binding protein (ABP) alpha, beta and gamma subunits. Further investigation of 14 alpha-like (Abpa) and 13 beta- or gamma-like (Abpbg) undisrupted gene sequences revealed a rich diversity of developmental stage-, sex- and tissue-specific expression. Despite these studies, our understanding of the evolution of this gene family remains incomplete. Questions arise from imperfections in the initial mouse genome assembly and a dearth of information about the gene family structure in other rodents and mammals.

RESULTS

Here, we interrogate the latest 'finished' mouse (Mus musculus) genome sequence assembly to show that the Abp gene repertoire is, in fact, twice as large as reported previously, with 30 Abpa and 34 Abpbg genes and pseudogenes. All of these have arisen since the last common ancestor with rat (Rattus norvegicus). We then demonstrate, by sequencing homologs from species within the Mus genus, that this burst of gene duplication occurred very recently, within the past seven million years. Finally, we survey Abp orthologs in genomes from across the mammalian clade and show that bursts of Abp gene duplications are not specific to the murid rodents; they also occurred recently in the lagomorph (rabbit, Oryctolagus cuniculus) and ruminant (cattle, Bos taurus) lineages, although not in other mammalian taxa.

CONCLUSION

We conclude that Abp genes have undergone repeated bursts of gene duplication and adaptive sequence diversification driven by these genes' participation in chemosensation and/or sexual identification.

Diverse spatial, temporal, and sexual expression of recently duplicated androgen-binding protein genes in Mus musculus

BACKGROUND

The genes for salivary androgen-binding protein (ABP) subunits have been evolving rapidly in ancestors of the house mouse Mus musculus, as evidenced both by recent and extensive gene duplication and by high ratios of nonsynonymous to synonymous nucleotide substitution rates. This makes ABP an appropriate model system with which to investigate how recent adaptive evolution of paralogous genes results in functional innovation (neofunctionalization).

RESULTS

It was our goal to find evidence for the expression of as many of the Abp paralogues in the mouse genome as possible. We observed expression of six Abpa paralogues and five Abpbg paralogues in ten glands and other organs located predominantly in the head and neck (olfactory lobe of the brain, three salivary glands, lacrimal gland, Harderian gland, vomeronasal organ, and major olfactory epithelium). These Abp paralogues differed dramatically in their specific expression in these different glands and in their sexual dimorphism of expression. We also studied the appearance of expression in both late-stage embryos and postnatal animals prior to puberty and found significantly different timing of the onset of expression among the various paralogues.

CONCLUSION

The multiple changes in the spatial expression profile of these genes resulting in various combinations of expression in glands and other organs in the head and face of the mouse strongly suggest that neofunctionalization of these genes, driven by adaptive evolution, has occurred following duplication. The extensive diversification in expression of this family of proteins provides two lines of evidence for a pheromonal role for ABP: 1) different patterns of Abpa/Abpbg expression in different glands; and 2) sexual dimorphism in the expression of the paralogues in a subset of those glands. These expression patterns differ dramatically among various glands that are located almost exclusively in the head and neck, where the sensory organs are located. Since mice are nocturnal, it is expected that they will make extensive use of olfactory as opposed to visual cues. The glands expressing Abp paralogues produce secretions (lacrimal and salivary) or detect odors (MOE and VNO) and thus it appears highly likely that ABP proteins play a role in olfactory communication.

Comparative evolutionary genomics of androgen-binding protein genes

Allelic variation within the mouse androgen-binding protein (ABP) alpha subunit gene (Abpa) has been suggested to promote assortative mating and thus prezygotic isolation. This is consistent with the elevated evolutionary rates observed for the Abpa gene, and the Abpb and Abpg genes whose products (ABPbeta and ABPgamma) form heterodimers with ABPalpha. We have investigated the mouse sequence that contains the three Abpa/b/g genes, and orthologous regions in rat, human, and chimpanzee genomes. Our studies reveal extensive "remodeling" of this region: Duplication rates of Abpa-like and Abpbg-like genes in mouse are >2 orders of magnitude higher than the average rate for all mouse genes; synonymous nucleotide substitution rates are twofold higher; and the Abpabg genomic region has expanded nearly threefold since divergence of the rodents. During this time, one in six amino acid sites in ABPbetagamma-like proteins appear to have been subject to positive selection; these may constitute a site of interaction with receptors or ligands. Greater adaptive variation among Abpbg-like sequences than among Abpa-like sequences suggests that assortative mating preferences are more influenced by variation in Abpbg-like genes. We propose a role for ABPalpha/beta/gamma proteins as pheromones, or in modulating odorant detection. This would account for the extraordinary adaptive evolution of these genes, and surrounding genomic regions, in murid rodents.