Search

Background Colorectal cancer (CRC) is a multifactorial disease resulting from both genetic predisposition and environmental factors including the gut microbiota (GM), but deciphering the influence of genetic variants, environmental variables, and interactions with the GM is exceedingly difficult. We previously observed significant differences in intestinal adenoma multiplicity between C57BL/6 J-ApcMin (B6-Min/J) from The Jackson Laboratory (JAX), and original founder strain C57BL/6JD-ApcMin (B6-Min/D) from the University of Wisconsin. Methods To resolve genetic and environmental interactions and determine their contributions we utilized two genetically inbred, independently isolated ApcMin mouse colonies that have been separated for over 20 generations. Whole genome sequencing was used to identify genetic variants unique to the two substrains. To determine the influence of genetic variants and the impact of differences in the GM on phenotypic variability, we used complex microbiota targeted rederivation to generate two Apc mutant mouse colonies harboring complex GMs from two different sources (GMJAX originally from JAX or GMHSD originally from Envigo), creating four ApcMin groups. Untargeted metabolomics were used to characterize shifts in the fecal metabolite profile based on genetic variation and differences in the GM. Results WGS revealed several thousand high quality variants unique to the two substrains. No homozygous variants were present in coding regions, with the vast majority of variants residing in noncoding regions. Host genetic divergence between Min/J and Min/D and the complex GM additively determined differential adenoma susceptibility. Untargeted metabolomics revealed that both genetic lineage and the GM collectively determined the fecal metabolite profile, and that each differentially regulates bile acid (BA) metabolism. Metabolomics pathway analysis facilitated identification of a functionally relevant private noncoding variant associated with the bile acid transporter Fatty acid binding protein 6 (Fabp6). Expression studies demonstrated differential expression of Fabp6 between Min/J and Min/D, and the variant correlates with adenoma multiplicity in backcrossed mice. Conclusions We found that both genetic variation and differences in microbiota influences the quantitiative adenoma phenotype in ApcMin mice. These findings demonstrate how the use of metabolomics datasets can aid as a functional genomic tool, and furthermore illustrate the power of a multi-omics approach to dissect complex disease susceptibility of noncoding variants.

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.

Recent analyses of the genetics of peripheral T-cell lymphoma (PTCL) have shown that a large proportion of cases are derived from normal follicular helper (Tfh) T-cells. The sanroque mouse strain bears a mutation that increases Tfh cell number and heterozygous animals (Roquinsan/+) develop lymphomas similar to human Tfh lymphoma. Here we demonstrate the usefulness of Roquinsan/+ animals as a pre-clinical model of Tfh lymphoma. Long latency of development and incomplete penetrance in this strain suggests the lymphomas are genetically diverse. We carried out preliminary genetic characterisation by whole exome sequencing and detected tumor specific mutations in Hsp90ab1, Ccnb3 and RhoA. Interleukin-2-inducible kinase (ITK) is expressed in Tfh lymphoma and is a potential therapeutic agent. A preclinical study of ibrutinib, a small molecule inhibitor of mouse and human ITK, in established lymphoma was carried out and showed lymphoma regression in 8/12 (67%) mice. Using T2-weighted MRI to assess lymph node volume and diffusion weighted MRI scanning as a measure of function, we showed that treatment increased mean apparent diffusion coefficient (ADC) suggesting cell death, and that change in ADC following treatment correlated with change in lymphoma volume. We suggest that heterozygous sanroque mice are a useful model of Tfh cell derived lymphomas in an immunocompetent animal.

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.

The most commonly employed mammalian model organism is the laboratory mouse. A wide variety of genetically diverse inbred mouse strains, representing distinct physiological states, disease susceptibilities, and biological mechanisms have been developed over the last century. We report full length draft de novo genome assemblies for 16 of the most widely used inbred strains and reveal for the first time extensive strain-specific haplotype variation. We identify and characterise 2,567 regions on the current Genome Reference Consortium mouse reference genome exhibiting the greatest sequence diversity between strains. These regions are enriched for genes involved in defence and immunity, and exhibit enrichment of transposable elements and signatures of recent retrotransposition events. Combinations of alleles and genes unique to an individual strain are commonly observed at these loci, reflecting distinct strain phenotypes. Several immune related loci, some in previously identified QTLs for disease response have novel haplotypes not present in the reference that may explain the phenotype. We used these genomes to improve the mouse reference genome resulting in the completion of 10 new gene structures, and 62 new coding loci were added to the reference genome annotation. Notably this high quality collection of genomes revealed a previously unannotated gene (Efcab3-like) encoding 5,874 amino acids, one of the largest known in the rodent lineage. Interestingly, Efcab3-like−/− mice exhibit severe size anomalies in four regions of the brain suggesting a mechanism of Efcab3-like regulating brain development.

We report full-length draft de novo genome assemblies for 16 widely used inbred mouse strains and find extensive strain-specific haplotype variation. We identify and characterize 2,567 regions on the current mouse reference genome exhibiting the greatest sequence diversity. These regions are enriched for genes involved in pathogen defence and immunity and exhibit enrichment of transposable elements and signatures of recent retrotransposition events. Combinations of alleles and genes unique to an individual strain are commonly observed at these loci, reflecting distinct strain phenotypes. We used these genomes to improve the mouse reference genome, resulting in the completion of 10 new gene structures. Also, 62 new coding loci were added to the reference genome annotation. These genomes identified a large, previously unannotated, gene (Efcab3-like) encoding 5,874 amino acids. Mutant Efcab3-like mice display anomalies in multiple brain regions, suggesting a possible role for this gene in the regulation of brain development. Medical Research Council and the Wellcome Trust

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 to 6 MYA, but that are absent in the Hominidae. In fact, Hominidae show between four-and seven-fold 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. For example, 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. This process resulted in thousands of novel, species-specific CTCF binding sites. Our results demonstrate that the comparison of matched phylogenetic sets of genomes will be an increasingly powerful strategy for understanding mammalian biology.

The Mouse Genomes Project was initiated in 2009 with the goal of using next-generation sequencing technologies to catalogue molecular variation in the common laboratory mouse strains, and a selected set of wild-derived inbred strains. The initial sequencing and survey of sequence variation in 17 inbred strains was completed in 2011 and included comprehensive catalogue of single nucleotide polymorphisms, short insertion/deletions, larger structural variants including their fine scale architecture and landscape of transposable element variation, and genomic sites subject to post-transcriptional alteration of RNA. From this beginning, the resource has expanded significantly to include 36 fully sequenced inbred laboratory mouse strains, a refined and updated data processing pipeline, and new variation querying and data visualisation tools which are available on the project's website ( http://www.sanger.ac.uk/resources/mouse/genomes/ ). The focus of the project is now the completion of de novo assembled chromosome sequences and strain-specific gene structures for the core strains. We discuss how the assembled chromosomes will power comparative analysis, data access tools and future directions of mouse genetics.

peer-reviewed peer-reviewed Background: Single nucleotide polymorphisms (SNPs) are the most abundant genetic variant found in vertebrates and invertebrates. SNP discovery has become a highly automated, robust and relatively inexpensive process allowing the identification of many thousands of mutations for model and non-model organisms. Annotating large numbers of SNPs can be a difficult and complex process. Many tools available are optimised for use with organisms densely sampled for SNPs, such as humans. There are currently few tools available that are species non-specific or support non-model organism data. Results: Here we present SNPdat, a high throughput analysis tool that can provide a comprehensive annotation of both novel and known SNPs for any organism with a draft sequence and annotation. Using a dataset of 4,566 SNPs identified in cattle using high-throughput DNA sequencing we demonstrate the annotations performed and the statistics that can be generated by SNPdat. Conclusions: SNPdat provides users with a simple tool for annotation of genomes that are either not supported by other tools or have a small number of annotated SNPs available. SNPdat can also be used to analyse datasets from organisms which are densely sampled for SNPs. As a command line tool it can easily be incorporated into existing SNP discovery pipelines and fills a niche for analyses involving non-model organisms that are not supported by many available SNP annotation tools. SNPdat will be of great interest to scientists involved in SNP discovery and analysis projects, particularly those with limited bioinformatics experience. Teagasc Walsh Fellowship Scheme; Science Foundation Ireland (SFI) Stokes lecturer scheme (number: 07/SK/B1236A)

Metastasis remains the primary cause of patient morbidity and mortality in solid tumors and is due to the action of a large number of tumor-autonomous and non-autonomous factors. Here we report the results of a genome-wide integrated strategy to identify novel metastasis susceptibility candidate genes and molecular pathways in breast cancer metastasis. This analysis implicates a number of transcriptional regulators and suggests cell-mediated immunity is an important determinant. Moreover, the analysis identified novel or FDA-approved drugs as potentially useful for anti-metastatic therapy. Further explorations implementing this strategy may therefore provide a variety of information for clinical applications in the control and treatment of advanced neoplastic disease., Author Summary Metastasis, the spread and growth of tumor cells from the original tumor to secondary sites throughout the body, is the primary cause of cancer-related death for most solid tumor types. The process of metastasis is very complex, requiring multiple individual steps and the cooperation of different cell types during the dissemination and proliferation steps. Many genes are involved in this process, but at present few have been identified and characterized. In this study, we have integrated multiple genome-wide analysis methods to try to identify large numbers of candidate metastasis-associated genes and pathways based on a highly metastatic mouse model. Using this strategy, we have identified a number of genes that predict outcome of human breast cancer. These genes implicate specific molecular and cellular pathways in the metastatic process that might be used to intervene in the process. Furthermore, this integrated analysis implicates pre-existing drugs that might be re-purposed to help prevent or reduce metastatic burden in patients. The combined results obtained from this analytical strategy therefore provide an important platform for further genome-wide analysis into the etiology of metastatic disease.

s: AACR Special Conference on Tumor Metastasis; November 30-December 3, 2015; Austin, TX Metastasis remains the primary cause of patient morbidity and mortality in solid tumors and is due to the action of a large number of tumor-autonomous and non-autonomous factors. Hundreds or thousands of genes are thought to be associated with metastasis however how many of these genes contribute etiologically to tumor progression is not currently known. Identification of the genes contributing mechanistically to the metastatic processes will not only deepen our understanding of how metastases occur, but also provide novel targets for either preventing their formation or combatting their life-threatening effects. Previously our laboratory demonstrated that inbred mice of distinct phylogenetic lineages possess distinct propensity for metastatic disease, indicating that inherited susceptibility for metastasis exists and that polymorphisms can both mark and functionally effect genes within the metastatic cascade. Here we report the results of a genome-wide integrated strategy to identify novel metastasis susceptibility candidate genes and molecular pathways which implicates a number of transcriptional regulators and suggests cell-mediated immunity is an important determinant. The strategy described integrates meiotic genetic screens with epigenetic control of gene expression and three dimensional chromatin structure analyses to identify genes, molecular and cellular pathways likely to be important in metastatic disease. Moreover, the analysis identified novel or FDA-approved drugs as potentially useful for anti-metastatic therapy. Further explorations implementing this strategy may therefore provide a variety of information for clinical applications in the control and treatment of advanced neoplastic disease. Citation Format: Kent W. Hunter, Ling Bai, Howard H. Yang, Ying Hu, Anjali Shukla, Ngoc-Han Ha, Anthony Doran, Farhoud Faraji, Natalie Goldberger, Maxwell P. Lee, Thomas Keane. An integrated genome-wide systems genetics screen for breast cancer metastasis susceptibility genes. [abstract]. In: Proceedings of the AACR Special Conference on Tumor Metastasis; 2015 Nov 30-Dec 3; Austin, TX. Philadelphia (PA): AACR; Cancer Res 2016;76(7 Suppl):Abstract nr IA21.

Breast cancer is the most frequent cancer and the second leading cause of cancer mortality in women. The vast majority of breast cancer mortality is due to metastatic disease since the primary tumor can be relatively easily surgically resected. More comprehensive understanding of biology of metastasis is therefore clearly warranted to unveil novel metastasis-associated molecules and cellular processes that might be targeted for clinical intervention. Previously we demonstrated that like cancer incidence, metastatic progression has a significant inherited component. Using mouse complex trait mapping strategies we have been isolating metastasis susceptibility genes from backcross analysis. However, studies suggest that many of these QTL peaks are the result of the contribution of multiple genes, suggesting that causal genes are not being identified. To increase our ability to detect causal genes we have implemented a new integrated strategy. Whole genome sequencing of appropriate high or low metastatic mouse strains is performed to identify variants, which are then filtered for those within DNAse hypersensitive sites (DHS), based on the hypothesis that most inherited phenotypes are due to expression differences rather than missense variants. The genes associated with polymorphic DHS are then screened through mouse and human tumor expression databases to identify those genes associated with development of metastatic disease. This strategy was able to identify genes previously associated with metastatic breast cancer in both mouse and humans and new genes have also been validated. The results further suggest that inherited variation in cellular mediated immune response may be an important contributor to metastatic disease. Citation Format: Anjali Shukla, Ling Bai, Howard Yang, Anthony Doran, Ying Hu, Thomas Geiger, Maxwell Lee, Thomas Keane, Kent W. Hunter. Integrating SNPs, epigenetics and transcriptomics to better understand the inherited predisposition to breast cancer metastasis. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4138. doi:10.1158/1538-7445.AM2015-4138

Background Four traits related to carcass performance have been identified as economically important in beef production: carcass weight, carcass fat, carcass conformation of progeny and cull cow carcass weight. Although Holstein-Friesian cattle are primarily utilized for milk production, they are also an important source of meat for beef production and export. Because of this, there is great interest in understanding the underlying genomic structure influencing these traits. Several genome-wide association studies have identified regions of the bovine genome associated with growth or carcass traits, however, little is known about the mechanisms or underlying biological pathways involved. This study aims to detect regions of the bovine genome associated with carcass performance traits (employing a panel of 54,001 SNPs) using measures of genetic merit (as predicted transmitting abilities) for 5,705 Irish Holstein-Friesian animals. Candidate genes and biological pathways were then identified for each trait under investigation. Results Following adjustment for false discovery (q-value 0.5) with at least one of the four traits. In total, 557 unique bovine genes, which mapped to 426 human orthologs, were within 500kbs of QTL found associated with a trait using the Bayesian approach. Using this information, 24 significantly over-represented pathways were identified across all traits. The most significantly over-represented biological pathway was the peroxisome proliferator-activated receptor (PPAR) signaling pathway. Conclusions A large number of genomic regions putatively associated with bovine carcass traits were detected using two different statistical approaches. Notably, several significant associations were detected in close proximity to genes with a known role in animal growth such as glucagon and leptin. Several biological pathways, including PPAR signaling, were shown to be involved in various aspects of bovine carcass performance. These core genes and biological processes may form the foundation for further investigation to identify causative mutations involved in each trait. Results reported here support previous findings suggesting conservation of key biological processes involved in growth and metabolism.