Project description:For decades, policies regarding generic medicines have sought to provide patients with economical access to safe and effective drugs, while encouraging the development of new therapies. This balance is becoming more challenging for physicians and regulators as biologics and non-biological complex drugs (NBCDs) such as glatiramer acetate demonstrate remarkable efficacy, because generics for these medicines are more difficult to assess. We sought to develop computational methods that use transcriptional profiles to compare branded medicines to generics, robustly characterizing differences in biological impact. We combined multiple computational methods to determine whether differentially expressed genes result from random variation, or point to consistent differences in biological impact of the generic compared to the branded medicine. We applied these methods to analyze gene expression data from mouse splenocytes exposed to either branded glatiramer acetate or a generic. The computational methods identified extensive evidence that branded glatiramer acetate has a more consistent biological impact across batches than the generic, and has a distinct impact on regulatory T cells and myeloid lineage cells. In summary, we developed a computational pipeline that integrates multiple methods to compare two medicines in an innovative way. This pipeline, and the specific findings distinguishing branded glatiramer acetate from a generic, can help physicians and regulators take appropriate steps to ensure safety and efficacy.
Project description:For decades, policies regarding generic medicines have sought to provide patients with economical access to safe and effective drugs, while encouraging the development of new therapies. This balance is becoming more challenging for physicians and regulators as biologics and non-biological complex drugs (NBCDs) such as glatiramer acetate demonstrate remarkable efficacy, because generics for these medicines are more difficult to assess. We sought to develop computational methods that use transcriptional profiles to compare branded medicines to generics, robustly characterizing differences in biological impact. We combined multiple computational methods to determine whether differentially expressed genes result from random variation, or point to consistent differences in biological impact of the generic compared to the branded medicine. We applied these methods to analyze gene expression data from mouse splenocytes exposed to either branded glatiramer acetate or a generic. The computational methods identified extensive evidence that branded glatiramer acetate has a more consistent biological impact across batches than the generic, and has a distinct impact on regulatory T cells and myeloid lineage cells. In summary, we developed a computational pipeline that integrates multiple methods to compare two medicines in an innovative way. This pipeline, and the specific findings distinguishing branded glatiramer acetate from a generic, can help physicians and regulators take appropriate steps to ensure safety and efficacy. Glatiramoid samples for SPL cell activation were grouped into four categories: 1) Verified GA, which included GA-RS (22 biological samples) and GA drug product (GA-DP, 34 biological samples from 30 batches) manufactured by Teva; 2) Deliberately Modified GA (DM-GA; 9 biological samples), which included glatiramoids made by Teva that were similar to GA but modified in a variety of ways: prepared with different ingredients (e.g., missing a constituent amino acid); prepared with the same amino acids in the same molar ratio as GA, but with defined amino acid sequences and different molecular weights (referred to as peptide markers TV-35 and TV-66); synthesized by a different process (e.g., changing acetolytic cleavage conditions, alternating polymerization initiator); exposed to destabilizing conditions (e.g., degradation by acid, base, and heat); 3) Unverified Glatiramoids, which included 4 biological samples, TV-5010 and 3 glatiramoids synthesized to be similar to GA but not manufactured using the Teva-patented manufacturing process; and 4) Unverified Generic GA, which included samples from 2 glatiramoids (M-bM-^@M-^\GA-QM-bM-^@M-^] 11 biological samples from 5 different batches, and 2 M-bM-^@M-^\GA-CM-bM-^@M-^] samples from a single batch) marketed as generic GA manufactured by companies other than Teva.
Project description:Introgressed variants from other species can be an important source of genetic variation because they may arise rapidly, can include multiple mutations on a single haplotype, and have often been pretested by selection in the species of origin. Although introgressed alleles are generally deleterious, several studies have reported introgression as the source of adaptive alleles-including the rodenticide-resistant variant of Vkorc1 that introgressed from Mus spretus into European populations of Mus musculus domesticus. Here, we conducted bidirectional genome scans to characterize introgressed regions into one wild population of M. spretus from Spain and three wild populations of M. m. domesticus from France, Germany, and Iran. Despite the fact that these species show considerable intrinsic postzygotic reproductive isolation, introgression was observed in all individuals, including in the M. musculus reference genome (GRCm38). Mus spretus individuals had a greater proportion of introgression compared with M. m. domesticus, and within M. m. domesticus, the proportion of introgression decreased with geographic distance from the area of sympatry. Introgression was observed on all autosomes for both species, but not on the X-chromosome in M. m. domesticus, consistent with known X-linked hybrid sterility and inviability genes that have been mapped to the M. spretus X-chromosome. Tract lengths were generally short with a few outliers of up to 2.7 Mb. Interestingly, the longest introgressed tracts were in olfactory receptor regions, and introgressed tracts were significantly enriched for olfactory receptor genes in both species, suggesting that introgression may be a source of functional novelty even between species with high barriers to gene flow.
Project description:Translational research is commonly performed in the C57B6/J mouse strain, chosen for its genetic homogeneity and phenotypic uniformity. Here, we evaluate the suitability of the white-footed deer mouse (Peromyscus leucopus) as a model organism for aging research, offering a comparative analysis against C57B6/J and diversity outbred (DO) Mus musculus strains. Our study includes comparisons of body composition, skeletal muscle function, and cardiovascular parameters, shedding light on potential applications and limitations of P. leucopus in aging studies. Notably, P. leucopus exhibits distinct body composition characteristics, emphasizing reduced muscle force exertion and a unique metabolism, particularly in fat mass. Cardiovascular assessments showed changes in arterial stiffness, challenging conventional assumptions and highlighting the need for a nuanced interpretation of aging-related phenotypes. Our study also highlights inherent challenges associated with maintaining and phenotyping P. leucopus cohorts. Behavioral considerations, including anxiety-induced responses during handling and phenotyping assessment, pose obstacles in acquiring meaningful data. Moreover, the unique anatomy of P. leucopus necessitates careful adaptation of protocols designed for Mus musculus. While showcasing potential benefits, further extensive analyses across broader age ranges and larger cohorts are necessary to establish the reliability of P. leucopus as a robust and translatable model for aging studies.
Project description:BackgroundCopy number variation is an important dimension of genetic diversity and has implications in development and disease. As an important model organism, the mouse is a prime candidate for copy number variant (CNV) characterization, but this has yet to be completed for a large sample size. Here we report CNV analysis of publicly available, high-density microarray data files for 351 mouse tail samples, including 290 mice that had not been characterized for CNVs previously.ResultsWe found 9634 putative autosomal CNVs across the samples affecting 6.87% of the mouse reference genome. We find significant differences in the degree of CNV uniqueness (single sample occurrence) and the nature of CNV-gene overlap between wild-caught mice and classical laboratory strains. CNV-gene overlap was associated with lipid metabolism, pheromone response and olfaction compared to immunity, carbohydrate metabolism and amino-acid metabolism for wild-caught mice and classical laboratory strains, respectively. Using two subspecies of wild-caught Mus musculus, we identified putative CNVs unique to those subspecies and show this diversity is better captured by wild-derived laboratory strains than by the classical laboratory strains. A total of 9 genic copy number variable regions (CNVRs) were selected for experimental confirmation by droplet digital PCR (ddPCR).ConclusionThe analysis we present is a comprehensive, genome-wide analysis of CNVs in Mus musculus, which increases the number of known variants in the species and will accelerate the identification of novel variants in future studies.
Project description:BackgroundLong terminal repeat (LTR) retrotransposons make up a large fraction of the typical mammalian genome. They comprise about 8% of the human genome and approximately 10% of the mouse genome. On account of their abundance, LTR retrotransposons are believed to hold major significance for genome structure and function. Recent advances in genome sequencing of a variety of model organisms has provided an unprecedented opportunity to evaluate better the diversity of LTR retrotransposons resident in eukaryotic genomes.ResultsUsing a new data-mining program, LTR_STRUC, in conjunction with conventional techniques, we have mined the GenBank mouse (Mus musculus) database and the more complete Ensembl mouse dataset for LTR retrotransposons. We report here that the M. musculus genome contains at least 21 separate families of LTR retrotransposons; 13 of these families are described here for the first time.ConclusionsAll families of mouse LTR retrotransposons are members of the gypsy-like superfamily of retroviral-like elements. Several different families of unrelated non-autonomous elements were identified, suggesting that the evolution of non-autonomy may be a common event. High sequence similarity between several LTR retrotransposons identified in this study and those found in distantly-related species suggests that horizontal transfer has been a significant factor in the evolution of mouse LTR retrotransposons.
Project description:House mice (Mus musculus) emit ultrasonic vocalizations (USVs), which are surprisingly complex and have features of bird song, but their functions are not well understood. Previous studies have reported mixed evidence on whether there are sex differences in USV emission, though vocalization rate or other features may depend upon whether potential receivers are of the same or opposite sex. We recorded the USVs of wild-derived adult house mice (F1 of wild-caught Mus musculus musculus), and we compared the vocalizations of males and females in response to a stimulus mouse of the same- or opposite-sex. To detect and quantify vocalizations, we used an algorithm that automatically detects USVs (Automatic Mouse Ultrasound Detector or A-MUD). We found high individual variation in USV emission rates (4 to 2083 elements/10 min trial) and a skewed distribution, with most mice (60%) emitting few (≤50) elements. We found no differences in the rates of calling between the sexes overall, but mice of both sexes emitted vocalizations at a higher rate and higher frequencies during opposite- compared to same-sex interactions. We also observed a trend toward higher amplitudes by males when presented with a male compared to a female stimulus. Our results suggest that mice modulate the rate and frequency of vocalizations depending upon the sex of potential receivers.
Project description:To characterize the genetic basis of hybrid male sterility in detail, we used a systems genetics approach, integrating mapping of gene expression traits with sterility phenotypes and QTL. We measured genome-wide testis expression in 305 male F2s from a cross between wild-derived inbred strains of M. musculus musculus and M. m. domesticus. We identified several thousand cis- and trans-acting QTL contributing to expression variation (eQTL). Many trans eQTL cluster into eleven ‘hotspots,’ seven of which co-localize with QTL for sterility phenotypes identified in the cross. The number and clustering of trans eQTL - but not cis eQTL - were substantially lower when mapping was restricted to a ‘fertile’ subset of mice, providing evidence that trans eQTL hotspots are related to sterility. Functional annotation of transcripts with eQTL provides insights into the biological processes disrupted by sterility loci and guides prioritization of candidate genes. Using a conditional mapping approach, we identified eQTL dependent on interactions between loci, revealing a complex system of epistasis. Our results illuminate established patterns, including the role of the X chromosome in hybrid sterility.
Project description:The mammalian vomeronasal organ (VNO) expresses two G-protein coupled receptor gene families that mediate pheromone responses, the V1R and V2R receptor genes. In rodents, there are ~150 V1R genes comprising 12 subfamilies organized in gene clusters at multiple chromosomal locations. Previously, we showed that several of these subfamilies had been extensively modulated by gene duplications, deletions, and gene conversions around the time of the evolutionary split of the mouse and rat lineages, consistent with the hypothesis that V1R repertoires might be involved in reinforcing speciation events. Here, we generated genome sequence for one large cluster containing two V1R subfamilies in Mus spretus, a closely related and sympatric species to Mus musculus, and investigated evolutionary change in these repertoires along the two mouse lineages.We describe a comparison of spretus and musculus with respect to genome organization and synteny, as well as V1R gene content and phylogeny, with reference to previous observations made between mouse and rat. Unlike the mouse-rat comparisons, synteny seems to be largely conserved between the two mouse species. Disruption of local synteny is generally associated with differences in repeat content, although these differences appear to arise more from deletion than new integrations. Even though unambiguous V1R orthology is evident, we observe dynamic modulation of the functional repertoires, with two of seven V1Rb and one of eleven V1Ra genes lost in spretus, two V1Ra genes becoming pseudogenes in musculus, two additional orthologous pairs apparently subject to strong adaptive selection, and another divergent orthologous pair that apparently was subjected to gene conversion.Therefore, eight of the 18 (~44%) presumptive V1Ra/V1Rb genes in the musculus-spretus ancestor appear to have undergone functional modulation since these two species diverged. As compared to the rat-mouse split, where modulation is evident by independent expansions of these two V1R subfamilies, divergence between musculus and spretus has arisen more by mutations within coding sequences. These results support the hypothesis that adaptive changes in functional V1R repertoires contribute to the delineation of very closely related species.
Project description:BACKGROUND:Identification of protein-protein interactions is an important first step to understand living systems. High-throughput experimental approaches have accumulated large amount of information on protein-protein interactions in human and other model organisms. Such interaction information has been successfully transferred to other species, in which the experimental data are limited. However, the annotation transfer method could yield false positive interologs due to the lack of conservation of interactions when applied to phylogenetically distant organisms. RESULTS:To address this issue, we used phylogenetic profile method to filter false positives in interologs based on the notion that evolutionary conserved interactions show similar patterns of occurrence along the genomes. The approach was applied to Mus musculus, in which the experimentally identified interactions are limited. We first inferred the protein-protein interactions in Mus musculus by using two approaches: i) identifying mouse orthologs of interacting proteins (interologs) based on the experimental protein-protein interaction data from other organisms; and ii) analyzing frequency of mouse ortholog co-occurrence in predicted operons of bacteria. We then filtered possible false-positives in the predicted interactions using the phylogenetic profiles. We found that this filtering method significantly increased the frequency of interacting protein-pairs coexpressed in the same cells/tissues in gene expression omnibus (GEO) database as well as the frequency of interacting protein-pairs shared the similar Gene Ontology (GO) terms for biological processes and cellular localizations. The data supports the notion that phylogenetic profile helps to reduce the number of false positives in interologs. CONCLUSION:We have developed protein-protein interaction database in mouse, which contains 41109 interologs. We have also developed a web interface to facilitate the use of database http://lgsun.grc.nia.nih.gov/mppi/.