Search

Comparison of the androgen-binding protein (Abp) gene regions of six Mus genomes provides insights into the evolutionary history of this large murid rodent gene family. We identified 206 unique Abp sequences and mapped their physical relationships. At least 48 are duplicated and thus present in more than two identical copies. All six taxa have substantially elevated LINE1 densities in Abp regions compared with flanking regions, similar to levels in mouse and rat genomes, although nonallelic homologous recombination seems to have only occurred in Mus musculus domesticus. Phylogenetic and structural relationships support the hypothesis that the extensive Abp expansion began in an ancestor of the genus Mus. We also found duplicated Abpa27's in two taxa, suggesting that previously reported selection on a27 alleles may have actually detected selection on haplotypes wherein different paralogs were lost in each. Other studies reported that a27 gene and species trees were incongruent, likely because of homoplasy. However, L1MC3 phylogenies, supposed to be homoplasy-free compared with coding regions, support our paralog hypothesis because the L1MC3 phylogeny was congruent with the a27 topology. This paralog hypothesis provides an alternative explanation for the origin of the a27 gene that is suggested to be fixed in the three different subspecies of Mus musculus and to mediate sexual selection and incipient reinforcement between at least two of them. Finally, we ask why there are so many Abp genes, especially given the high frequency of pseudogenes and suggest that relaxed selection operates over a large part of the gene clusters., (© The Author(s) 2021. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Understanding the mechanisms driving lineage-specific evolution in both primates and rodents has been hindered by the lack of sister clades with a similar phylogenetic structure having high-quality genome assemblies. Here, we have created chromosome-level assemblies of the Mus caroli and Mus pahari genomes. Together with the Mus musculus and Rattus norvegicus genomes, this set of rodent genomes is similar in divergence times to the Hominidae (human-chimpanzee-gorilla-orangutan). By comparing the evolutionary dynamics between the Muridae and Hominidae, we identified punctate events of chromosome reshuffling that shaped the ancestral karyotype of Mus musculus and Mus caroli between 3 and 6 million yr ago, but that are absent in the Hominidae. Hominidae show between four- and sevenfold lower rates of nucleotide change and feature turnover in both neutral and functional sequences, suggesting an underlying coherence to the Muridae acceleration. Our system of matched, high-quality genome assemblies revealed how specific classes of repeats can play lineage-specific roles in related species. Recent LINE activity has remodeled protein-coding loci to a greater extent across the Muridae than the Hominidae, with functional consequences at the species level such as reproductive isolation. Furthermore, we charted a Muridae-specific retrotransposon expansion at unprecedented resolution, revealing how a single nucleotide mutation transformed a specific SINE element into an active CTCF binding site carrier specifically in Mus caroli , which resulted in thousands of novel, species-specific CTCF binding sites. Our results show that the comparison of matched phylogenetic sets of genomes will be an increasingly powerful strategy for understanding mammalian biology., (© 2018 Thybert et al.; Published by Cold Spring Harbor Laboratory Press.)

The Androgen-binding protein ( Abp ) gene region of the mouse genome contains 64 genes, some encoding pheromones that influence assortative mating between mice from different subspecies. Using CNVnator and quantitative PCR, we explored copy number variation in this gene family in natural populations of Mus musculus domesticus ( Mmd ) and Mus musculus musculus ( Mmm ), two subspecies of house mice that form a narrow hybrid zone in Central Europe. We found that copy number variation in the center of the Abp gene region is very common in wild Mmd , primarily representing the presence/absence of the final duplications described for the mouse genome. Clustering of Mmd individuals based on this variation did not reflect their geographical origin, suggesting no population divergence in the Abp gene cluster. However, copy number variation patterns differ substantially between Mmd and other mouse taxa. Large blocks of Abp genes are absent in Mmm , Mus musculus castaneus and an outgroup, Mus spretus , although with differences in variation and breakpoint locations. Our analysis calls into question the reliance on a reference genome for interpreting the detailed organization of genes in taxa more distant from the Mmd reference genome. The polymorphic nature of the gene family expansion in all four taxa suggests that the number of Abp genes, especially in the central gene region, is not critical to the survival and reproduction of the mouse. However, Abp haplotypes of variable length may serve as a source of raw genetic material for new signals influencing reproductive communication and thus speciation of mice., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

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., (Copyright © 2017 by the Genetics Society of America.)

Retrotransposons comprise a large portion of mammalian genomes. They contribute to structural changes and more importantly to gene regulation. The expansion and diversification of gene families have been implicated as sources of evolutionary novelties. Given the roles retrotransposons play in genomes, their contribution to the evolution of gene families warrants further exploration. In this study, we found a significant association between two major retrotransposon classes, LINEs and LTRs, and lineage-specific gene family expansions in both the human and mouse genomes. The distribution and diversity differ between LINEs and LTRs, suggesting that each has a distinct involvement in gene family expansion. LTRs are associated with open chromatin sites surrounding the gene families, supporting their involvement in gene regulation, whereas LINEs may play a structural role promoting gene duplication. Our findings also suggest that gene family expansions, especially in the mouse genome, undergo two phases. The first phase is characterized by elevated deposition of LTRs and their utilization in reshaping gene regulatory networks. The second phase is characterized by rapid gene family expansion due to continuous accumulation of LINEs and it appears that, in some instances at least, this could become a runaway process. We provide an example in which this has happened and we present a simulation supporting the possibility of the runaway process. Altogether we provide evidence of the contribution of retrotransposons to the expansion and evolution of gene families. Our findings emphasize the putative importance of these elements in diversification and adaptation in the human and mouse lineages., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

We produced a tear proteome of the genome mouse, C57BL/6, that contained 139 different protein identifications: 110 from a two-dimensional (2D) gel with subsequent trypsin digestion, 19 from a one-dimensional (1D) gel with subsequent trypsin digestion and ten from a 1D gel with subsequent Asp-N digestion. We compared this tear proteome with a C57BL/6 mouse saliva proteome produced previously. Sixteen of the 139 tear proteins are shared between the two proteomes, including six proteins that combat microbial growth. Among the 123 other tear proteins, were members of four large protein families that have no counterparts in humans: Androgen-binding proteins (ABPs) with different members expressed in the two proteomes, Exocrine secreted peptides (ESPs) expressed exclusively in the tear proteome, major urinary proteins (MUPs) expressed in one or both proteomes and the mouse-specific Kallikreins (subfamily b KLKs) expressed exclusively in the saliva proteome. All four families have members with suggested roles in mouse communication, which may influence some aspect of reproductive behavior. We discuss this in the context of functional adaptation involving tear and saliva proteins in the secretions of mouse lacrimal and salivary glands, respectively.

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.

The overall goal of our study was to compare the proteins found in the saliva proteomes of three mammals: human, mouse and rat. Our first objective was to compare two human proteomes with very different analysis depths. The 89 shared proteins in this comparison apparently represent a core of highly-expressed human salivary proteins. Of the proteins unique to each proteome, one-half to 2/3 lack signal peptides and probably are contaminants instead of less highly-represented salivary proteins. We recently published the first rodent saliva proteomes with salivas collected from the genome mouse (C57BL/6) and the genome rat (BN/SsNHsd/Mcwi). Our second objective was to compare the proteins in the human proteome with those we identified in the genome mouse and rat to determine those common to all three mammals as well as the specialized rodent subset. We also identified proteins unique to each of the three mammals because differences in the secreted protein constitutions can provide clues to differences in the evolutionary adaptation of the secretions in the three different mammals., Competing Interests: The authors declare no conflict of interest.

Background: Retrotransposons have been suggested to provide a substrate for non-allelic homologous recombination (NAHR) and thereby promote gene family expansion. Their precise role, however, is controversial. Here we ask whether retrotransposons contributed to the recent expansions of the Androgen-binding protein (Abp) gene families that occurred independently in the mouse and rat genomes., Results: Using dot plot analysis, we found that the most recent duplication in the Abp region of the mouse genome is flanked by L1Md_T elements. Analysis of the sequence of these elements revealed breakpoints that are the relicts of the recombination that caused the duplication, confirming that the duplication arose as a result of NAHR using L1 elements as substrates. L1 and ERVII retrotransposons are considerably denser in the Abp regions than in one Mb flanking regions, while other repeat types are depleted in the Abp regions compared to flanking regions. L1 retrotransposons preferentially accumulated in the Abp gene regions after lineage separation and roughly followed the pattern of Abp gene expansion. By contrast, the proportion of shared vs. lineage-specific ERVII repeats in the Abp region resembles the rest of the genome., Conclusions: We confirmed the role of L1 repeats in Abp gene duplication with the identification of recombinant L1Md_T elements at the edges of the most recent mouse Abp gene duplication. High densities of L1 and ERVII repeats were found in the Abp gene region with abrupt transitions at the region boundaries, suggesting that their higher densities are tightly associated with Abp gene duplication. We observed that the major accumulation of L1 elements occurred after the split of the mouse and rat lineages and that there is a striking overlap between the timing of L1 accumulation and expansion of the Abp gene family in the mouse genome. Establishing a link between the accumulation of L1 elements and the expansion of the Abp gene family and identification of an NAHR-related breakpoint in the most recent duplication are the main contributions of our study.

Two decades ago, we developed a congenic strain of Mus musculus, called b-congenic, by replacing the androgen-binding protein Abpa27(a) allele in the C3H/HeJ genome with the Abpa27(b) allele from DBA/2J. We and other researchers used this b-congenic strain and its C3H counterpart, the a-congenic strain, to test the hypothesis that, given the choice between signals from two strains with different a27 alleles on the same genetic background, test subjects would prefer the homosubspecific one. It was our purpose in undertaking this study to characterize the segment transferred from DBA to the C3H background in producing the b-congenic strain on which a role for ABPA27 in behavior has been predicated. We determined the size of the chromosome 7 segment transferred from DBA and the genes it contains that might influence preference. We found that the "functional" DBA segment is about 1% the size of the mouse haploid genome and contains at least 29 genes expressed in salivary glands, however, only three of these encode proteins identified in the mouse salivary proteome. At least two of the three genes Abpa27, Abpbg26 and Abpbg27 encoding the subunits of androgen-binding protein ABP dimers evolved under positive selection and the third one may have also. In the sense that they are subunits of the same two functional entities, the ABP dimers, we propose that their evolutionary histories might not be independent of each other.

Three proteinaceous pheromone families, the androgen-binding proteins (ABPs), the exocrine-gland secreting peptides (ESPs) and the major urinary proteins (MUPs) are encoded by large gene families in the genomes of Mus musculus and Rattus norvegicus. We studied the evolutionary histories of the Mup and Esp genes and compared them with what is known about the Abp genes. Apparently gene conversion has played little if any role in the expansion of the mouse Class A and Class B Mup genes and pseudogenes, and the rat Mups. By contrast, we found evidence of extensive gene conversion in many Esp genes although not in all of them. Our studies of selection identified at least two amino acid sites in β-sheets as having evolved under positive selection in the mouse Class A and Class B MUPs and in rat MUPs. We show that selection may have acted on the ESPs by determining K(a)/K(s) for Exon 3 sequences with and without the converted sequence segment. While it appears that purifying selection acted on the ESP signal peptides, the secreted portions of the ESPs probably have undergone much more rapid evolution. When the inner gene converted fragment sequences were removed, eleven Esp paralogs were present in two or more pairs with K(a)/K(s) >1.0 and thus we propose that positive selection is detectable by this means in at least some mouse Esp paralogs. We compare and contrast the evolutionary histories of all three mouse pheromone gene families in light of their proposed functions in mouse communication.

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.

Assortative mating, a potentially efficient prezygotic reproductive barrier, may prevent loss of genetic potential by avoiding the production of unfit hybrids (i.e., because of hybrid infertility or hybrid breakdown) that occur at regions of secondary contact between incipient species. In the case of the mouse hybrid zone, where two subspecies of Mus musculus (M. m. domesticus and M. m. musculus) meet and exchange genes to a limited extent, assortative mating requires a means of subspecies recognition. We based the work reported here on the hypothesis that, if there is a pheromone sufficiently diverged between M. m. domesticus and M. m. musculus to mediate subspecies recognition, then that process must also require a specific receptor(s), also sufficiently diverged between the subspecies, to receive the signal and elicit an assortative mating response. We studied the mouse V1R genes, which encode a large family of receptors in the vomeronasal organ (VNO), by screening Perlegen SNP data and identified one, Vmn1r67, with 24 fixed SNP differences most of which (15/24) are nonsynonymous nucleotide substitutions between M. m. domesticus and M. m. musculus. We observed substantial linkage disequilibrium (LD) between Vmn1r67 and Abpa27, a mouse salivary androgen-binding protein gene that encodes a proteinaceous pheromone (ABP) capable of mediating assortative mating, perhaps in conjunction with its bound small lipophilic ligand. The LD we observed is likely a case of association rather than residual physical linkage from a very recent selective sweep, because an intervening gene, Vmn1r71, shows significant intra(sub)specific polymorphism but no inter(sub)specific divergence in its nucleotide sequence. We discuss alternative explanations of these observations, for example that Abpa27 and Vmn1r67 are coevolving as signal and receptor to reinforce subspecies hybridization barriers or that the unusually divergent Vmn1r67 allele was not a product of fast positive selection, but was derived from an introgressed allele, possibly from Mus spretus.

Three families of proteinaceous pheromones have been described in the house mouse: androgen-binding proteins (ABPs), exocrine gland-secreting peptides (ESPs), and major urinary proteins (MUPs), each of which is thought to communicate different information. All three are encoded by large gene clusters in different regions of the mouse genome, clusters that have expanded dramatically during mouse evolutionary history. We report copy number variation among the most recently duplicated Abp genes, which suggests substantial volatility in this gene region. It appears that groups of these genes behave as low copy repeats (LCRs), duplicating as relatively large blocks of genes by nonallelic homologous recombination. An analysis of gene conversion suggested that it did not contribute to the very low or absent divergence among the paralogs duplicated in this way. We evaluated the ESP and MUP gene regions for signs of the LCR pattern but could find no compelling evidence for duplication of gene blocks of any significant size. Assessment of the entire Abp gene region with the Mouse Paralogy Browser supported the conclusion that substantial volatility has occurred there. This was especially evident when comparing strains with all or part of the Mus musculus musculus or Mus musculus castaneus Abp region. No particularly remarkable volatility was observed in the other two gene families, and we discuss the significance of this in light of the various roles proposed for the three families of mouse proteinaceous pheromones.

Male reproductive fitness is strongly affected by seminal fluid. In addition to interacting with the female environment, seminal fluid mediates important physiological characteristics of sperm, including capacitation and motility. In mammals, the male reproductive tract shows a striking degree of compartmentalization, with at least six distinct tissue types contributing material that is combined with sperm in an ejaculate. Although studies of whole ejaculates have been undertaken in some species, we lack a comprehensive picture of the specific proteins produced by different accessory tissues. Here, we perform proteomic investigations of six regions of the male reproductive tract in mice -- seminal vesicles, anterior prostate, dorsolateral prostate, ventral prostate, bulbourethral gland, and bulbourethral diverticulum. We identify 766 proteins that could be mapped to 506 unique genes and compare them with a high-quality human seminal fluid data set. We find that Gene Ontology functions of seminal proteins are largely conserved between mice and humans. By placing these data in an evolutionary framework, we show that seminal vesicle proteins have experienced a significantly higher rate of nonsynonymous substitution compared with the genome, which could be the result of adaptive evolution. In contrast, proteins from the other five tissues showed significantly lower nonsynonymous substitution, revealing a previously unappreciated level of evolutionary constraint acting on the majority of male reproductive proteins.

Proteins involved in reproductive fitness have evolved unusually rapidly across diverse groups of organisms. These reproductive proteins show unusually high rates of amino acid substitutions, suggesting that the proteins have been subject to positive selection. We sought to identify seminal fluid proteins experiencing adaptive evolution because such proteins are often involved in sperm competition, host immunity to pathogens, and manipulation of female reproductive physiology and behavior. We performed an evolutionary screen of the mouse prostate transcriptome for genes with elevated evolutionary rates between mouse and rat. We observed that secreted rodent prostate proteins evolve approximately twice as fast as nonsecreted proteins, remarkably similar to findings in the primate prostate and in the Drosophila male accessory gland. Our screen led us to identify and characterize a group of seminal vesicle secretion (Svs) proteins and to show that the gene Svs7 is evolving very rapidly, with many amino acid sites under positive selection. Another gene in this group, Svs5, showed evidence of branch-specific selection in the rat. We also found that Svs7 is under selection in primates and, by using three-dimensional models, demonstrated that the same regions have been under selection in both groups. Svs7 has been identified as mouse caltrin, a protein involved in sperm capacitation, the process responsible for the timing of changes in sperm activity and behavior, following ejaculation. We propose that the most likely explanation of the adaptive evolution of Svs7 that we have observed in rodents and primates stems from an important function in sperm competition.

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.

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.

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., (Copyright 2004 Cold Spring Harbor Laboratory Press ISSN)

Previous behavioral experiments showed that mouse salivary androgen-binding protein (ABP) was involved in interindividual recognition and might play a role in sexual isolation between house mouse (Mus musculus) subspecies. The pattern of evolution of Abpa, the gene for the alpha subunit of ABP, was found to be consistent with this hypothesis. Abpa apparently diverged rapidly between species and subspecies with a large excess of nonsynonymous substitutions, a lack of exon polymorphism within each of the three subspecies, and a lack of intron polymorphism in the one subspecies studied (M. musculus domesticus). Here we characterized the intron and exon sequence variations of this gene in house mouse populations from central Eurasia, a region yet unsampled and thought to be close to the cradle of the radiation of the subspecies. We also determined the intron and exon sequences in seven other species of the genus Mus. We confirmed the general pattern of rapid evolution by essentially nonsynonymous substitutions, both inter- and intraspecifically, supporting the idea that Darwinian selection has driven the evolution of this gene. We also observed a uniform intron sequence in five samples of M. musculus musculus, suggesting that a selective sweep might have occurred for that allele. In contrast to previous results, however, we found extensive intron and exon polymorphism in some house mouse populations from central Eurasia. We also found evidence for secondary admixture of the subspecies-specific alleles in regions of transition between the subspecies in central Eurasia. Furthermore, an abnormal intron phylogeny suggested that interspecific exchanges had occurred between the house mouse subspecies and three other Palearctic species. These observations appear to be at variance with the simple hypothesis that Abpa is involved in reproductive isolation. Although we do not rule out a role in recognition, the situation appears to be more complex than previously thought. Thus the selective mechanism behind the evolution of Abpa remains to be resolved, and we suggest that it may have changed during the recent colonization history of the house mouse.

We studied the effects of a single genetic change on a complex mammalian behavior using animals congenic for two variants of Abpa, the gene for the alpha subunit of mouse salivary androgen-binding protein (ABP), in two-way preference tests. Females exhibited a preference for investigating salivas of males of their own genetic type of ABP but not for urines of either type of male. This preference behavior is consistent for samples of mice from geographically diverse populations of Mus musculus domesticus and M. m. musculus. These findings provide an explanation for the observation that this gene is evolving under strong selection.

Previous work has shown that the gene for the alpha subunit of androgen-binding protein, Abpa, may be involved in premating isolation between different subspecies of the house mouse, Mus musculus. We investigated patterns of DNA sequence variation at Abpa within and between species of mice to test several predictions of a model of neutral molecular evolution. Intraspecific variation among 10 Mus musculus domesticus alleles was compared with divergence between M. m. domesticus and M. caroli for Abpa and two X-linked genes, Glra2 and Amg. No variation was observed at Abpa within M. m. domesticus. The ratio of polymorphism to divergence was significantly lower at Abpa than at Glra2 and Amg, despite the fact that all three genes experience similar rates of recombination. Interspecific comparisons among M. m. domesticus, Mus musculus musculus, Mus musculus castaneus, Mus spretus, Mus spicilegus, and Mus caroli revealed that the ratio of nonsynonymous substitutions to synonymous substitutions on a per-site basis (Ka/Ks) was generally greater than one. The combined observations of no variation at Abpa within M. m: domesticus and uniformly high Ka/Ks values between species suggest that positive directional selection has acted recently at this locus.

Mouse salivary androgen-binding protein (ABP) is a family of dimeric proteins that may play a pheromonal role in Mus musculus. The protein dimer consists of a common alpha subunit disulfide-bonded to a variable (beta or gamma) subunit. Here we report N-terminal sequences of the beta and gamma subunits, showing that they are very similar to each other while being quite different from the alpha subunit. We demonstrate differential androgen binding by the two dimers. Both bind dihydrotestosterone to about the same extent but the alpha:beta dimer binds significantly more testosterone than the alpha:gamma dimer. We discuss the possible significance of this diversity of androgen binding with respect to the possibility that androgen binding is related to a putative pheromonal role for the protein.

The goals of this study were to determine the steroid-binding specificity of the mouse salivary androgen-binding protein (ABP) family and to ascertain whether there might be other proteins in mouse saliva capable of binding steroids. The optimal conditions for testosterone binding by mouse salivary proteins were determined using a small-scale chromatography system to separate bound from unbound steroid. Testosterone binding appeared to be biphasic but was directly proportional to saliva concentration, with an optimum temperature of 37 degrees C in the second phase. These results were used to develop a steroid-binding protocol to study the steroid specificity of salivary proteins separated by electrophoresis. The ABP family bound testosterone and progesterone well and HO-progesterone and DHT to a lesser extent but did not bind either cholesterol or estradiol. Steroid structural comparisons suggest that binding by ABP is governed by the A ring of the steroid. Another protein that is not a member of the ABP family bound cholesterol specifically but no protein that specifically bound estradiol was observed.

We wanted to determine whether the microevolution of the mouse salivary androgen-binding protein (ABP) Alpha subunit gene (Abpa) could mediate sexual selection and thereby have a potential role in maintaining gene pool integrity where radiating mouse subspecies make secondary contact. This hypothesis is based upon previous work in this laboratory, which has shown that each subspecies apparently has its own allele and that these alleles have a 25-fold excess of nonsynonymous/synonymous base substitutions compared to an average protein under purifying selection. We provide direct evidence for ABP-assortative mate selection in a laboratory setting: Mus musculus domesticus and M. m. musculus female mice recognize and discriminate between the territories of male mice that essentially differ solely in their Abpa genotype and, when the males are present, the female prefers to mate with the one of her own ABP type. The observation that females could differentiate between the territories of the two males when those mice were absent suggests that the males marked their territories with ABP. In this study, we also detected ABP on the pelts of male mice and in their environment. It is likely that the animals apply the protein to their pelts by licking and that it is then deposited in their surroundings. We suggest that females of the two subspecies are able to discriminate between males of those subspecies on the basis of this protein molecule. Mouse salivary ABP might present a worthwhile system with which to study a prezygotic isolation mechanism in a mammal., (© 1997 The Society for the Study of Evolution.)

Mouse salivary androgen-binding protein (ABP) is a major secretory product of the submaxillary glands. Although it is a common salivary protein among rodents generally, the function of ABP has yet to be determined. Here we report a comparison of the DNA coding sequences and putative amino acid sequences they determine for the three common alleles of the Alpha subunit gene (Abpa), alleles that appear to be diagnostic for the three subspecies of Mus musculus. Three other unique sequences were found in the species M. caroli, M. spretus, and M. spicilegus. Comparison of the six sequences shows that 8 of the 20 base substitution sites produce a high degree of variability in amino acids 32, 33, 36, and 39, a variability that creates unique sequence combinations in each species and subspecies. We compare the possibilities that selection or genetic drift caused this unusual microevolution and argue that selection is the more likely explanation. We speculate on the potential significance of this with respect to the proposal that ABP is involved in assortive mate kin selection.

Androgen-binding protein (ABP) is found in the salivas of a wide variety of rodents and it has been proposed that ABP functions in sex and/or subspecies recognition (Karn and Dlouhy, J. Hered. 82, 453, 1991). This is a report of significant identity between the alpha subunit of mouse salivary ABP and Chain 1 of cat allergen Fel dI (50% identity), as well as with two other proteins that share identity with Chain 1 of Fel dI, rabbit uteroglobin (27% identity with ABP alpha) and human lung Clara 10 (27% identity with ABP alpha). The secondary structure predicted for the mouse ABP alpha subunit is a very good fit with the secondary structure determined by X-ray crystallography for rabbit uteroglobin, a protein that shares with mouse ABP the capability of binding steroid. Fel dI is found in cat saliva, sebaceous glands, and pelt. Its function is not known but it has been proposed to be involved in protecting dry epithelia, a parallel to uteroglobin protecting wet epithelia. Since mice, like cats, lick themselves and each other extensively, coating their pelts with ABP may be part of this or another biological function.

It has been suggested that three distinct genes, Abpa, Abpb, and Abpg, determine the three subunits of mouse salivary androgen-binding protein (ABP) (Dlouhy, S. R., et al., Genetics 115, 535, 1987). We report the putative amino acid sequence of the subunit common to all forms of ABP, the Alpha subunit, and the partial amino acid sequence of the Beta subunit. These sequences have little in common, supporting the notion of at least two distinct genes coding for the subunits of the most common form of salivary ABP, the A:B dimer. A search of GenBank showed that these sequences have not been reported previously. The Beta subunit shows significant homology with helospectin, a member of the glucagon superfamily, but not enough homology to assign it to the family. No homology exists between ABP subunits and members of the ligand-binding carrier family of proteins nor does ABP show homology with other androgen-binding proteins. Particularly interesting is the observation that there is no relationship to rat prostatic steroid binding protein (PBP), given the similarities in protein tertiary structure, the numbers of subunits and their genes, and the earlier observation of ABP cross-reactive material in mouse prostate.

We have shown that genes for at least six human parotid proteins, parotid acidic protein (Pa), proline-rich protein (Pr), double-banded protein (Db), glycoprotein (Gl), parotid middle-band protein (Pm), and parotid-size variant (Ps) are linked. We have designated this complex of genes as the salivary protein complex (SPC). Several of the genes in this complex show marked associations that are most likely the result of linkage disequilibrium. It seems likely that the SPC arose through the process of gene duplication. This hypothesis is supported by the results of our present study that demonstrate the biochemical similarity of the protein products of several SPC genes. The amino acid compositions of the major SPC proteins are compared, including several (Ps 1 and 2, and Db) that have not been published. All of these proteins are quite similar and consist to a large extent of the amino acids, proline, glycine, and gix (glutamine and/or glutamic acid).

We demonstrate that the previously described gene Androgen binding protein (Abp; Dlouhy and Karn, 1984) codes for the Alpha subunit of ABP and rename the locus Androgen binding protein alpha (Abpa). A study of recombinant inbred strains demonstrates that Abpa is located on chromosome 7 near Glucose phosphate isomerase-1 (Gpi-1). Biochemical and genetic evidence indicates the existence of another ABP subunit, Gamma, and its locus, Androgen binding protein gamma (Abpg), that is closely linked to Abpa. Although no polymorphism has yet been found for the previously described Beta subunit of ABP (Dlouhy and Karn, 1983; 1984), we suggest that it represents a third locus. Androgen binding protein beta (Abpb). ABP subunits appear to dimerize randomly and a model is presented demonstrating the origin of six ABP dimers in the salivas of AbpaaAbpga/AbpabAbpgb heterozygous mice. The results of cell-free translation studies in which the pre-ABP subunits are identified specifically by immunoprecipitation with anti-ABP antibody supports the idea that independent mRNAs code for the Alpha, Beta and Gamma subunits.

This report describes a gene which influences the electrophoretic mobility of a protein in the salivas of adult mice. Three categories of phenotype have been observed: the two single-banded types, F (Fast) and S (Slow), and the two-banded type, SF (Slow-Fast), with the two bands represented in varying proportions. All females, regardless of age or strain, and all males before puberty show only the F phenotype. Males of the BALB/c and C57BL/6J strains show the F phenotype throughout puberty and adult life, whereas males of the C3H/St and C57BL/KsJ strains show the SF phenotype in puberty and the S phenotype in adult life. We have designated this variation the sex-limited saliva pattern (Ssp). The results from genetic crosses indicate that the variation among the strains is determined by an autosomal locus, Ssp, with two alleles, SspS and SspF, where SspS is dominant to SspF. Testosterone treatment can accelerate the acquisition of the S type in males of the strains C3H/St and C57BL/KsJ and also induces that phenotype in C3H/St females and C57BL/6J males. Thus it appears that the observed strain-specific differences reflect a genetic variation in androgen levels and/or androgen sensitivity rather than variation in a structural gene.

The allelically determined human salivary proteins, Ps 1 and 2, were purified on sodium dodecyl sulfate gels, eluted, and compared by limited proteolysis with Streptomyces griseus protease VI, Bacillus subtilis protease VII, and Staphylococcus aureus protease V8. Prior dansylation of the Ps isoproteins facilitated visualization of the peptides. Digestion patterns indicate considerable homology between the Ps isoproteins and support the conclusion [Azen, A. E., and Denniston, C. (1980). Biochem. Genet. 18:483] that there is an actual molecular weight difference between them. Further, the results suggest that this difference owes to an extension of the Ps 2 chain at one of its ends.

Human proline-rich proteins (PRPs) are polymorphic, homologous in sequence, and linked in a cluster called the human salivary protein complex (SPC). Recently this complex was localized to human chromosome band 12p13.2 (Mamula et al., Cytogenet. Cell Genet. 39:279, 1985). We have isolated a PRP cDNA, EO27, from a human parotid gland library, identified it by DNA sequencing, and used it to study the molecular and cellular biology of PRP production. Cell-free translation and mRNA characterization with EO27 indicate that the numerous PRPs seen in saliva are produced from relatively few, large precursors, probably by posttranslational cleavage. This supports an hypothesis originally proposed by Friedman and Karn in 1977 (Am. J. Hum. Genet. 29:44 A; Biochem. Genet. 15:549) and later supported by biochemical studies (Karn et al., Biochem Genet. 17:1061, 1979) and molecular studies (Mamula et al., Fed. Proc. 43:1522, 1984; Maeda et al., J. Biol. Chem. 260:1123, 1985). EO27 was also used in this study to localize PRP mRNA production to the acinar cells of the parotid gland by in situ hybridization.

The autosomal recessive gene diabetes (db) produces a condition similar to human insulin-dependent diabetes mellitus in certain strains of inbred mice. In this investigation, the effects of expression of the db gene on the development of the submandibular glands, electrophoretic protein patterns in salivas, fasting blood glucose levels, and glycosylated hemoglobin levels were evaluated in mice undergoing puberty. Three sex-limited effects of the db gene were observed in diabetic male mice: (1) a compromise of the development of the specialized submandibular glands with the extensive tubular portion normally found in males. (2) failure to develop a salivary protein pattern unique to male mice, and (3) attainment of higher levels of fasting blood glucose than found in female diabetic mice. Since it has been documented that homozygous mice fail to develop functional gonads, apparently due to insufficient production of gonadotropin, it is likely that the compromised development of the specialized submandibular glands, and, consequently, the male salivary protein pattern, is a result of decreased testosterone production. Experiments in which diabetic mice were treated with testosterone support that conclusion, since testosterone caused transformation of the salivary protein pattern to one identical with that of normal male littermate controls and increased the tubular portion of the submandibular glands.

A novel variation in electrophoretic phenotype is described for a mouse salivary androgen binding protein (Abp). Crosses show that the variation is inherited in an autosomal codominant manner and protein characterization studies show that the variant Abp differs in isoelectric point from the common form of the protein. Those observations suggest that the variation involves the structural gene for the mouse salivary Abp. The genetic studies also show that the electrophoretic mobility of the variant Abp can be influenced by the sex-limited saliva pattern (Ssp) gene. The SspS allele alters the electrophoretic mobility of Abp in males at puberty or in females which have received exogenous testosterone [Karn, R.C., Dlouhy, S.R., Springer, K.R., Hjorth, J.P., and Nielsen, J.T. (1982). Biochem Genet. 20:493]. This study shows that Abp and Ssp are distinct genes which are not closely linked and that SspS is trans active in F1 (Abpa/Abpb, SspS/SspF) males.

Based on data from 76 informative families, linkages between Pa and Pr and between Pr and Db have been established by two-point linkage analysis. In both pairs of loci, there were no significant sex heterogeneity in recombination fractions. Linkage between Pa and Db cannot be established based on two-point analyses, but a significant sex difference in the recombination fraction between Pa and Db was observed. Strong confirming evidence was obtained from three-point analysis to place Pa, Pr, and Db in one linkage group. The most likely order is Pa-Pr-Db, but the relative odds over second order Pr-Pa-Db are small. Haplotype frequencies of Pr, Pa, and Db were obtained based on the phenotypes of the 685 random Caucasians, providing evidence for marked linkage disequilibrium among the three loci.

A new polymorphism, Pc, has been identified in human saliva. Two proteins, Pc 1 and Pc 2, are determined by alleles Pc1 and Pc2, respectively, which show autosomal codominant inheritance. No null phenotype has been encountered in 225 randomly collected salivas. The frequencies of the two alleles differ in the Black and White American populations, with Pc1 and Pc2 being 0.670 and 0.330 in the Black (N = 47) and 0.461 and 0.539 in the White (N = 178) populations, respectively. The alleles are in equilibrium in the two populations and segregation analyses (30 families) do not suggest the existence of a null allele in either population. Of seven polymorphic human salivary proteins determined by genes in the salivary protein complex (SPC), Pc phenotypes show association only with Ps phenotypes. Based on that association, our linkage studies, and the biochemical similarities with other SPC proteins, we tentatively conclude that Pc is a member of the SPC, bringing the total number of genes in that group to 13.

Expression of the dominant allele, Ssps, at the Sex-limited saliva pattern locus in mice results in an alteration of the electrophoretic mobility of a mouse salivary protein from the F (fast) to the S (slow) type [Karn, R.C., et al. (1980). Genetics 94:s52; Karn, R.C., et al. (1982). Biochem. Genet. 20:493]. We now demonstrate that the protein affected binds androgen and has a basic heterodimeric structure with the subunits connected by disulfide bridging. It is produced in the submaxillary gland where the alteration takes place. The change in electrophoretic mobility from F to S appears to be primarily the result of an increase in the molecular weight of the larger subunit since the isoelectric point changes very little. We discuss the possible causes of the increase in molecular weight.

The acidic proline-rich proteins (Pr) showing genetic polymorphism were purified from human parotid salivas by gel filtration and ion exchange chromatography. Molecular weight determinations, amino acid composition analyses, and polypeptide mapping experiments indicate that the Pr 3 protein is a fragment of the Pr 1 protein. Studies of a parotid saliva factor capable of converting Pr 1 to Pr 3 and Pr 2 to Pr 4 indicate that Pr 3 and Pr 4 are generated from Pr 1 and Pr 2, respectively. Evidence suggests that the converting factor is a protease capable of posttranslationally cleaving Pr 1 and Pr 2, the primary or derived products of alleles Pr1 and Pr2.

A genetic analysis of human urinary pepsinogen isozymes is presented. Nine discrete phenotypes were identified in a population survey of 215 unrelated Caucasian individuals. The phenotypes were characterized by differences among the staining intensities of the activated group I pepsinogens, Pg 5, Pg 4, Pg 3, and Pg 2. The genetic studies demonstrated that the codominant expression of four alleles, PgA, PgB, PgC and PgD, at a single genetic locus determined the nine phenotypes identified. Linkage analysis excluded close linkage of the Pg locus with the chromosome 6 markers HLA, GLO1, and Bf.

Isoelectric focusing was performed on parotid salivas selected for their electrophoretic phenotypes of proline-rich acidic salivary proteins. Fractions encompassing narrow pH regions were pooled and examined by polyacrylamide gel electrophoresis. Isoelectric focusing yielded partial purification of major and minor acidic proline-rich proteins which were subsequently compared by immunoelectrophoresis and double immunodiffusion against goat anti-human parotid saliva. Cross-reactivity without spurring between all fractions containing major Pr proteins in both immunoelectrophoresis and double immunodiffusion suggests that these proteins are immunologically very similar or identical.

We previously demonstrated polymorphism of a mouse salivary protein which, because of its ability to bind androgen, we designated androgen binding protein (ABP) and its structural gene, Androgen binding protein (Abp). This report describes the purification of salivary ABP and presents the amino acid composition of its subunits. Using an antibody raised against the purified protein, we demonstrated the presence of cross-reactive material (CRM) in mouse parotid, submaxillary, sublingual, and prostate glands by double immunodiffusion. Immunohistochemical detection of proteins on electroblots of polyacrylamide electrophoresis and isoelectric focusing gels, however shows that the prostate CRM is a protein with electrophoretic and isoelectric focusing behavior distinct from that of salivary ABP isoproteins. Because the DBA/2J strain has a variant of salivary ABP, that strain was analyzed to determine if a prostate ABP-CRM variant was present. The approach was hampered by an inability to detect CRM in the prostates of DBA/2J mice. Prostate CRM was detected, however, in some progeny from repeated backcrosses of DBA/2J X C3H/St hybrids to the C3H/St and DBA/2J parental strains. The prostate CRM detected in samples from animals heterozygous for salivary Abp appears to be identical to C3H/St prostate CRM, suggesting that the gene controlling prostate ABP-CRM is related to, but distinct from, Abp. The reason for reduced or absent CRM in the prostates of DBA/2J males is unknown but this finding suggests that there are strain-related differences in the expression of this prostate protein.

The age-dependent development of alpha-amylase expression in utero and during the first two years of life is reported. Separation of salivary and pancreatic amylase isozymes in a discontinuous buffered sheet polyacrylamide electrophoretic system, with subsequent densitometry, provides a reliable semiquantitative method of estimating the proportions of salivary and pancreatic amylases in urine and amniotic fluid samples. In the newborn the predominant amylase isozymes seen in the urine are of salivary origin. As the child ages the level of amylase in the urine rises and an increase in the proportion of pancreatic amylase isozymes occurs. Amniotic fluids of late first and early second trimester pregnancies contain salivary isozymes. None of the amniotic fluid samples examined has pancreatic amylase isozymes. These data reflect a differential development of the expression of the two amylase approaches adult levels by 16 months of age. Conversely, the salivary (Amy1) locus is expressed as early as 18 weeks of gestation and remains relatively constant with but a small increase in salivary amylase (units/ml) activity during early development, as the total amylase activity approaches adult values.