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

Evidence of Nonrandom Patterns of Functional Chromosome Organization in Danaus plexippus

Our understanding on the interplay between gene functionality and gene arrangement at different chromosome scales relies on a few Diptera and the honeybee, species with quality reference genome assemblies, accurate gene annotations, and abundant transcriptome data. Using recently generated 'omic resources in the monarch butterfly Danaus plexippus, a species with many more and smaller chromosomes relative to Drosophila species and the honeybee, we examined the organization of genes preferentially expressed at broadly defined developmental stages (larva, pupa, adult males, and adult females) at both fine and whole-chromosome scales. We found that developmental stage-regulated genes do not form more clusters, but do form larger clusters, than expected by chance, a pattern consistent across the gene categories examined. Notably, out of the 30 chromosomes in the monarch genome, 12 of them, plus the fraction of the chromosome Z that corresponds to the ancestral Z in other Lepidoptera, were found enriched for developmental stage-regulated genes. These two levels of nonrandom gene organization are not independent as enriched chromosomes for developmental stage-regulated genes tend to harbor disproportionately large clusters of these genes. Further, although paralogous genes were overrepresented in gene clusters, their presence is not enough to explain two-thirds of the documented cases of whole-chromosome enrichment. The composition of the largest clusters often included paralogs from more than one multigene family as well as unrelated single-copy genes. Our results reveal intriguing patterns at the whole-chromosome scale in D. plexippus while shedding light on the interplay between gene expression and chromosome organization beyond Diptera and Hymenoptera.

Androgen-Binding Protein (Abp) Evolutionary History: Has Positive Selection Caused Fixation of Different Paralogs in Different Taxa of the Genus Mus?

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.

Lacrimal androgen-binding proteins protect against Aspergillus fumigatus keratitis in mice

AIM

To clarify the regulatory mechanisms of lacrimal androgen-binding proteins (ABPs) in mice with keratitis caused by Aspergillus fumigatus (A. fumigatus).

METHODS

Mouse models of A. fumigatus keratitis were established. Lacrimal glands were removed after 24 h for general and histological comparison. Lacrimal ABPs were detected by qRT-PCR and quantitative proteomic analysis, or were detected by qRT-PCR after subconjunctival or lacrimal gland injection with dexamethasone. Unique inflammatory factors were detected by qRT-PCR, Western blot and/or immunofluorescence. Interleukin-1β (IL-1β) was injected into the lacrimal gland to explore the relationship between IL-1β and lacrimal ABPs.

RESULTS

The lacrimal glands of mice with fungal keratitis were larger than normal mice and these structures became disorganized. Moreover, the expression of ABP ε and ABP δ were increased. Subconjunctival injection with dexamethasone could reduce the size of the lacrimal gland and increase the expression of ABP ε and ABP δ, while lacrimal gland injection with dexamethasone had no obvious effects. The expression of IL-1β in the lacrimal gland of mice with A. fumigatus keratitis was increased. When IL-1β was injected into the lacrimal gland, the lacrimal gland enlarged and the expression of ABP ε and ABP δ decreased.

CONCLUSION

Lacrimal glands contributed to protection in fungal keratitis, which was not due to the involvement of inflammatory cells in mice. ABP δ and ABP ε of mice were involved in reducing the severity of corneal damage in mice with A. fumigatus keratitis. Moreover, the expression of IL-1β and ABP δ and ABP ε were intrinsically linked.

Functional evolution of vertebrate sensory receptors

Sensory receptors enable animals to perceive their external world, and functional properties of receptors evolve to detect the specific cues relevant for an organism's survival. Changes in sensory receptor function or tuning can directly impact an organism's behavior. Functional tests of receptors from multiple species and the generation of chimeric receptors between orthologs with different properties allow for the dissection of the molecular basis of receptor function and identification of the key residues that impart functional changes in different species. Knowledge of these functionally important sites facilitates investigation into questions regarding the role of epistasis and the extent of convergence, as well as the timing of sensory shifts relative to other phenotypic changes. However, as receptors can also play roles in non-sensory tissues, and receptor responses can be modulated by numerous other factors including varying expression levels, alternative splicing, and morphological features of the sensory cell, behavioral validation can be instrumental in confirming that responses observed in heterologous systems play a sensory role. Expression profiling of sensory cells and comparative genomics approaches can shed light on cell-type specific modifications and identify other proteins that may affect receptor function and can provide insight into the correlated evolution of complex suites of traits. Here we review the evolutionary history and diversity of functional responses of the major classes of sensory receptors in vertebrates, including opsins, chemosensory receptors, and ion channels involved in temperature-sensing, mechanosensation and electroreception.

The expression of lacrimal androgen-binding proteins in mice Pseudomonas aeruginosa keratitis

AIM

To investigate the expression of lacrimal androgen-binding proteins (ABPs) in mice Pseudomonas aeruginosa (P. aeruginosa) keratitis.

METHODS

P. aeruginosa mice model from different gender was developed by intra-stromal injection. The expression of lacrimal ABPs in lacrimal gland specimens from P. aeruginosa keratitis mice was detected by the quantitative polymerase chain reaction (qRT-PCR). Corneal virulence was evaluated based on clinical scores. To study the mechanism of lacrimal ABPs' expression, experimental subjects were pre-treated with 4E-BP1 inhibitor, and were used to evaluate the expression levels by qRT-PCR.

RESULTS

Compared with control groups, the expression of ABPα, ABPη and ABPζ in lacrimal gland from P. aeruginosa keratitis mice had no meaningful changes, while ABPε and ABPδ were significantly higher at 1d after infection. The expression of ABPδ in lacrimal gland of male mice was higher than female mice, regardless of whether or not P. aeruginosa keratitis occurred. After 4E-BP1 inhibitor subconjunctival injection or lacrimal injection, the expression of ABPδ and ABPε has no significant change compared with the control group.

CONCLUSION

ABPδ and ABPε secreted by mice lacrimal gland may involve in the progress of alleviating the severity of corneal damage in P. aeruginosa keratitis. The expression of ABPδ and ABPε upon P. aeruginosa infection is independent of cap-dependent mRNA translation activated by 4E-BP1.

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.

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.

Repeat associated mechanisms of genome evolution and function revealed by the Mus caroli and Mus pahari genomes

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.

Analysis of Copy Number Variation in the Abp Gene Regions of Two House Mouse Subspecies Suggests Divergence during the Gene Family Expansions

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.

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.

Immune and biochemical responses in skin differ between bovine hosts genetically susceptible and resistant to the cattle tick Rhipicephalus microplus

Background

Ticks attach to and penetrate their hosts’ skin and inactivate multiple components of host responses in order to acquire a blood meal. Infestation loads with the cattle tick, Rhipicephalus microplus, are heritable: some breeds carry high loads of reproductively successful ticks, whereas in others, few ticks feed and reproduce efficiently.

Methods

In order to elucidate the mechanisms that result in the different outcomes of infestations with cattle ticks, we examined global gene expression and inflammation induced by tick bites in skins from one resistant and one susceptible breed of cattle that underwent primary infestations with larvae and nymphs of R. microplus. We also examined the expression profiles of genes encoding secreted tick proteins that mediate parasitism in larvae and nymphs feeding on these breeds.

Results

Functional analyses of differentially expressed genes in the skin suggest that allergic contact-like dermatitis develops with ensuing production of IL-6, CXCL-8 and CCL-2 and is sustained by HMGB1, ISG15 and PKR, leading to expression of pro-inflammatory chemokines and cytokines that recruit granulocytes and T lymphocytes. Importantly, this response is delayed in susceptible hosts. Histopathological analyses of infested skins showed inflammatory reactions surrounding tick cement cones that enable attachment in both breeds, but in genetically tick-resistant bovines they destabilized the cone. The transcription data provided insights into tick-mediated activation of basophils, which have previously been shown to be a key to host resistance in model systems. Skin from tick-susceptible bovines expressed more transcripts encoding enzymes that detoxify tissues. Interestingly, these enzymes also produce volatile odoriferous compounds and, accordingly, skin rubbings from tick-susceptible bovines attracted significantly more tick larvae than rubbings from resistant hosts. Moreover, transcripts encoding secreted modulatory molecules by the tick were significantly more abundant in larval and in nymphal salivary glands from ticks feeding on susceptible bovines.

Conclusions

Compared with tick-susceptible hosts, genes encoding enzymes producing volatile compounds exhibit significantly lower expression in resistant hosts, which may render them less attractive to larvae; resistant hosts expose ticks to an earlier inflammatory response, which in ticks is associated with significantly lower expression of genes encoding salivary proteins that suppress host immunity, inflammation and coagulation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1945-z) contains supplementary material, which is available to authorized users.

Whole Genome Sequence of Two Wild-Derived Mus musculus domesticus Inbred Strains, LEWES/EiJ and ZALENDE/EiJ, with Different Diploid Numbers

Wild-derived mouse inbred strains are becoming increasingly popular for complex traits analysis, evolutionary studies, and systems genetics. Here, we report the whole-genome sequencing of two wild-derived mouse inbred strains, LEWES/EiJ and ZALENDE/EiJ, of Mus musculus domesticus origin. These two inbred strains were selected based on their geographic origin, karyotype, and use in ongoing research. We generated 14× and 18× coverage sequence, respectively, and discovered over 1.1 million novel variants, most of which are private to one of these strains. This report expands the number of wild-derived inbred genomes in the Mus genus from six to eight. The sequence variation can be accessed via an online query tool; variant calls (VCF format) and alignments (BAM format) are available for download from a dedicated ftp site. Finally, the sequencing data have also been stored in a lossless, compressed, and indexed format using the multi-string Burrows-Wheeler transform. All data can be used without restriction.

Comparative Proteomics of Mouse Tears and Saliva: Evidence from Large Protein Families for Functional Adaptation

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.

Molecular and phenotypic distinction of the very recently evolved insular subspecies Mus musculus helgolandicus ZIMMERMANN, 1953

Background

Populations and subspecies of the house mouse Mus musculus were able to invade new regions worldwide in the wake of human expansion. Here we investigate the origin and colonization history of the house mouse inhabiting the small island of Heligoland on the German Bight - Mus musculus helgolandicus. It was first described by Zimmermann in 1953, based on morphological descriptions which were considered to be a mosaic between the subspecies M. m. domesticus and M. m. musculus. Since mice on islands are excellent evolutionary model systems, we have focused here on a molecular characterization and an extended phenotype analysis.

Results

The molecular data show that the mice from Heligoland are derived from M. m. domesticus based on mitochondrial D-loop sequences as well as on four nuclear diagnostic markers, including one each from the sex-chromosomes. STRUCTURE analysis based on 21 microsatellite markers assigns Heligoland mice to a distinct population and D-loop network analysis suggests that they are derived from a single colonization event. In spite of mice from the mainland arriving by ships, they are apparently genetically refractory against further immigration. Mutation frequencies in complete mitochondrial genome sequences date the colonization age to approximately 400 years ago. Complete genome sequences from three animals revealed a genomic admixture with M. m. musculus genomic regions with at least 6.5 % of the genome affected. Geometric morphometric analysis of mandible shapes including skull samples from two time points during the last century suggest specific adaptations to a more carnivorous diet.

Conclusions

The molecular and morphological analyses confirm that M. m. helgolandicus consists of a distinct evolutionary lineage with specific adaptations. It shows a remarkable resilience against genetic mixture with mainland populations of M. m. domesticus despite major disturbances in the past century and a high ship traffic. The genomic admixture with M. m. musculus genetic material may have contributed to the genomic distinction of the Heligoland mice. In spite of its young age, M. m. helgolandicus may thus be considered as a true subspecies of Mus, whose evolution was triggered through fast divergence on a small island.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0439-5) contains supplementary material, which is available to authorized users.

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.

The role of retrotransposons in gene family expansions: insights from the mouse Abp gene family

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.

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.

The roles of gene duplication, gene conversion and positive selection in rodent Esp and Mup pheromone gene families with comparison to the Abp family

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.

Update of the human secretoglobin (SCGB) gene superfamily and an example of 'evolutionary bloom' of androgen-binding protein genes within the mouse Scgb gene superfamily

The secretoglobins (SCGBs) comprise a family of small, secreted proteins found in animals exclusively of mammalian lineage. There are 11 human SCGB genes and five pseudogenes. Interestingly, mice have 68 Scgb genes, four of which are highly orthologous to human SCGB genes; the remainder represent an 'evolutionary bloom' and make up a large gene family represented by only six counterparts in humans. SCGBs are found in high concentrations in many mammalian secretions, including fluids of the lung, lacrimal gland, salivary gland, prostate and uterus. Whereas the biological activities of most individual SCGBs have not been fully characterised, what already has been discovered suggests that this family has an important role in the modulation of inflammation, tissue repair and tumorigenesis. In mice, the large Scgb1b and Scgb2b gene families encode the androgen-binding proteins, which have been shown to play a role in mate selection. Although much has been learned about SCGBs in recent years, clearly more research remains to be done to allow a better understanding of the roles of these proteins in human health and disease. Such information is predicted to reveal valuable novel drug targets for the treatment of inflammation, as well as designing biomarkers that might identify tissue damage or cancer.

Congenic strain analysis reveals genes that are rapidly evolving components of a prezygotic isolation mechanism mediating incipient reinforcement

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.

Isolation and gene flow: inferring the speciation history of European house mice

Inferring the history of isolation and gene flow during species differentiation can inform us on the processes underlying their formation. Following their recent expansion in Europe, two subspecies of the house mouse (Mus musculus domesticus and Mus musculus musculus) have formed a hybrid zone maintained by hybrid incompatibilities and possibly behavioural reinforcement, offering a good model of incipient speciation. We reconstruct the history of their divergence using an approximate Bayesian computation framework and sequence variation at 57 autosomal loci. We find support for a long isolation period preceding the advent of gene flow around 200,000 generations ago, much before the formation of the European hybrid zone a few thousand years ago. The duration of the allopatric episode appears long enough (74% of divergence time) to explain the accumulation of many post-zygotic incompatibilities expressed in the present hybrid zone. The ancient contact inferred could have played a role in mating behaviour divergence and laid the ground for further reinforcement. We suggest that both subspecies originally colonized the Middle East from the northern Indian subcontinent, domesticus settling on the shores of the Persian Gulf and musculus on those of the Caspian Sea. Range expansions during interglacials would have induced secondary contacts, presumably in Iran, where they must have also interacted with Mus musculus castaneus. Future studies should incorporate this possibility, and we point to Iran and its surroundings as a hot spot for house mouse diversity and speciation studies.

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.

Reinforcement selection acting on the European house mouse hybrid zone

Behavioural isolation may lead to complete speciation when partial postzygotic isolation acts in the presence of divergent-specific mate-recognition systems. These conditions exist where Mus musculus musculus and M. m. domesticus come into contact and hybridize. We studied two mate-recognition signal systems, based on urinary and salivary proteins, across a Central European portion of the mouse hybrid zone. Introgression of the genomic regions responsible for these signals: the major urinary proteins (MUPs) and androgen binding proteins (ABPs), respectively, was compared to introgression at loci assumed to be nearly neutral and those under selection against hybridization. The preference of individuals taken from across the zone regarding these signals was measured in Y mazes, and we develop a model for the analysis of the transition of such traits under reinforcement selection. The strongest assortative preferences were found in males for urine and females for ABP. Clinal analyses confirm nearly neutral introgression of an Abp locus and two loci closely linked to the Abp gene cluster, whereas two markers flanking the Mup gene region reveal unexpected introgression. Geographic change in the preference traits matches our reinforcement selection model significantly better than standard cline models. Our study confirms that behavioural barriers are important components of reproductive isolation between the house mouse subspecies.

A novel defensive mechanism against acetaminophen toxicity in the mouse lateral nasal gland: role of CYP2A5-mediated regulation of testosterone homeostasis and salivary androgen-binding protein expression

To identify novel factors or mechanisms that are important for the resistance of tissues to chemical toxicity, we have determined the mechanisms underlying the previously observed increases in resistance to acetaminophen (APAP) toxicity in the lateral nasal gland (LNG) of the male Cyp2g1-null/Cyp2a5-low mouse. Initial studies established that Cyp2a5-null mice, but not a newly generated strain of Cyp2g1-null mice, were resistant to APAP toxicity in the LNG; therefore, subsequent studies were focused on the Cyp2a5-null mice. Compared with the wild-type (WT) male mouse, the Cyp2a5-null male mouse had intact capability to metabolize APAP to reactive intermediates in the LNG, as well as unaltered circulating levels of APAP, APAP-GSH, APAP-glucuronide, and APAP-sulfate. However, it displayed reduced tissue levels of APAP and APAP-GSH and increased tissue levels of testosterone and salivary androgen-binding protein (ABP) in the LNG. Furthermore, we found that ABP was able to compete with GSH and cellular proteins for adduction with reactive metabolites of APAP in vitro. The amounts of APAP-ABP adducts formed in vivo were greater, whereas the amounts of APAP adducts formed with other cellular proteins were substantially lower, in the LNG of APAP-treated male Cyp2a5-null mice compared with the LNG of APAP-treated male WT mice. We propose that through its critical role in testosterone metabolism, CYP2A5 regulates 1) the bioavailability of APAP and APAP-GSH (presumably through modulation of the rates of xenobiotic excretion from the LNG) and 2) the expression of ABP, which can quench reactive APAP metabolites and thereby spare critical cellular proteins from inactivation.

A candidate subspecies discrimination system involving a vomeronasal receptor gene with different alleles fixed in M. m. domesticus and M. m. musculus

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.

The mechanism of expansion and the volatility it created in three pheromone gene clusters in the mouse (Mus musculus) genome

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.

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.

Separating derived from ancestral features of mouse and human genomes

To take full advantage of the mouse as a model organism, it is essential to distinguish lineage-specific biology from what is shared between human and mouse. Investigations into shared genetic elements common to both have been well served by the draft human and mouse genome sequences. More recently, the virtually complete euchromatic sequences of the two reference genomes have been finished. These reveal a high ( approximately 5%) level of sequence duplications that had previously been recalcitrant to sequencing and assembly. Within these duplications lie large numbers of rodent- or primate-specific genes. In the present paper, we review the sequence properties of the two genomes, dwelling most on the duplications, deletions and insertions that separate each of them from their most recent common ancestor, approx. 90 million years ago. We consider the differences in gene numbers and repertoires between the two species, and speculate on their contributions to lineage-specific biology. Loss of ancient single-copy genes are rare, as are gains of new functional genes through retrotransposition. Instead, most changes to the gene repertoire have occurred in large multicopy families. It has been proposed that numbers of such 'environmental genes' rise and fall, and their sequences change, as adaptive responses to infection and other environmental pressures, including conspecific competition. Nevertheless, many such genes may be under little or no selection.

Lineage-specific biology revealed by a finished genome assembly of the mouse

A finished clone-based assembly of the mouse genome reveals extensive recent sequence duplication during recent evolution and rodent-specific expansion of certain gene families. Newly assembled duplications contain protein-coding genes that are mostly involved in reproductive function.

The mouse (Mus musculus) is the premier animal model for understanding human disease and development. Here we show that a comprehensive understanding of mouse biology is only possible with the availability of a finished, high-quality genome assembly. The finished clone-based assembly of the mouse strain C57BL/6J reported here has over 175,000 fewer gaps and over 139 Mb more of novel sequence, compared with the earlier MGSCv3 draft genome assembly. In a comprehensive analysis of this revised genome sequence, we are now able to define 20,210 protein-coding genes, over a thousand more than predicted in the human genome (19,042 genes). In addition, we identified 439 long, non–protein-coding RNAs with evidence for transcribed orthologs in human. We analyzed the complex and repetitive landscape of 267 Mb of sequence that was missing or misassembled in the previously published assembly, and we provide insights into the reasons for its resistance to sequencing and assembly by whole-genome shotgun approaches. Duplicated regions within newly assembled sequence tend to be of more recent ancestry than duplicates in the published draft, correcting our initial understanding of recent evolution on the mouse lineage. These duplicates appear to be largely composed of sequence regions containing transposable elements and duplicated protein-coding genes; of these, some may be fixed in the mouse population, but at least 40% of segmentally duplicated sequences are copy number variable even among laboratory mouse strains. Mouse lineage-specific regions contain 3,767 genes drawn mainly from rapidly-changing gene families associated with reproductive functions. The finished mouse genome assembly, therefore, greatly improves our understanding of rodent-specific biology and allows the delineation of ancestral biological functions that are shared with human from derived functions that are not.

The availability of an accurate genome sequence provides the bedrock upon which modern biomedical research is based. Here we describe a high-quality assembly, Build 36, of the mouse genome. This assembly was put together by aligning overlapping individual clones representing parts of the genome, and it provides a more complete picture than previous assemblies, because it adds much rodent-specific sequence that was previously unavailable. The addition of these sequences provides insight into both the genomic architecture and the gene complement of the mouse. In particular, it highlights recent gene duplications and the expansion of certain gene families during rodent evolution. An improved understanding of the mouse genome and thus mouse biology will enhance the utility of the mouse as a model for human disease.

Species specificity in major urinary proteins by parallel evolution

Species-specific chemosignals, pheromones, regulate social behaviors such as aggression, mating, pup-suckling, territory establishment, and dominance. The identity of these cues remains mostly undetermined and few mammalian pheromones have been identified. Genetically-encoded pheromones are expected to exhibit several different mechanisms for coding 1) diversity, to enable the signaling of multiple behaviors, 2) dynamic regulation, to indicate age and dominance, and 3) species-specificity. Recently, the major urinary proteins (Mups) have been shown to function themselves as genetically-encoded pheromones to regulate species-specific behavior. Mups are multiple highly related proteins expressed in combinatorial patterns that differ between individuals, gender, and age; which are sufficient to fulfill the first two criteria. We have now characterized and fully annotated the mouse Mup gene content in detail. This has enabled us to further analyze the extent of Mup coding diversity and determine their potential to encode species-specific cues.

Our results show that the mouse Mup gene cluster is composed of two subgroups: an older, more divergent class of genes and pseudogenes, and a second class with high sequence identity formed by recent sequential duplications of a single gene/pseudogene pair. Previous work suggests that truncated Mup pseudogenes may encode a family of functional hexapeptides with the potential for pheromone activity. Sequence comparison, however, reveals that they have limited coding potential. Similar analyses of nine other completed genomes find Mup gene expansions in divergent lineages, including those of rat, horse and grey mouse lemur, occurring independently from a single ancestral Mup present in other placental mammals. Our findings illustrate that increasing genomic complexity of the Mup gene family is not evolutionarily isolated, but is instead a recurring mechanism of generating coding diversity consistent with a species-specific function in mammals.

Adaptive evolution in rodent seminal vesicle secretion 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.