Genetic architecture of skewed X inactivation in the laboratory mouse

Distributed X chromosome inactivation in brain circuitry is associated with X-linked disease penetrance of behavior

The precise anatomical degree of brain X chromosome inactivation (XCI) that is sufficient to alter X-linked disorders in females is unclear. Here, we quantify whole-brain XCI at single-cell resolution to discover a prevalent activation ratio of maternal to paternal X at 60:40 across all divisions of the adult brain. This modest, non-random XCI influences X-linked disease penetrance: maternal transmission of the fragile X mental retardation 1 (Fmr1)-knockout (KO) allele confers 55% of total brain cells with mutant X-active, which is sufficient for behavioral penetrance, while 40% produced from paternal transmission is tolerated. Local XCI mosaicism within affected maternal Fmr1-KO mice further specifies sensorimotor versus social anxiety phenotypes depending on which distinct brain circuitry is most affected, with only a 50%-55% mutant X-active threshold determining penetrance. Thus, our results define a model of X-linked disease penetrance in females whereby distributed XCI among single cells populating brain circuitries can regulate the behavioral penetrance of an X-linked mutation.

Population variability in X-chromosome inactivation across 9 mammalian species

One of the two X chromosomes in female mammals is epigenetically silenced in embryonic stem cells by X chromosome inactivation (XCI). This creates a mosaic of cells expressing either the maternal or the paternal X allele. The XCI ratio, the proportion of inactivated parental alleles, varies widely among individuals, representing the largest instance of epigenetic variability within mammalian populations. While various contributing factors to XCI variability are recognized, namely stochastic and/or genetic effects, their relative contributions are poorly understood. This is due in part to limited cross-species analysis, making it difficult to distinguish between generalizable or species-specific mechanisms for XCI ratio variability. To address this gap, we measured XCI ratios in nine mammalian species (9,143 individual samples), ranging from rodents to primates, and compared the strength of stochastic models or genetic factors for explaining XCI variability. Our results demonstrate the embryonic stochasticity of XCI is a general explanatory model for population XCI variability in mammals, while genetic factors play a minor role.

Conservation and divergence of canonical and non-canonical imprinting in murids

BACKGROUND

Genomic imprinting affects gene expression in a parent-of-origin manner and has a profound impact on complex traits including growth and behavior. While the rat is widely used to model human pathophysiology, few imprinted genes have been identified in this murid. To systematically identify imprinted genes and genomic imprints in the rat, we use low input methods for genome-wide analyses of gene expression and DNA methylation to profile embryonic and extraembryonic tissues at allele-specific resolution.

RESULTS

We identify 14 and 26 imprinted genes in these tissues, respectively, with 10 of these genes imprinted in both tissues. Comparative analyses with mouse reveal that orthologous imprinted gene expression and associated canonical DNA methylation imprints are conserved in the embryo proper of the Muridae family. However, only 3 paternally expressed imprinted genes are conserved in the extraembryonic tissue of murids, all of which are associated with non-canonical H3K27me3 imprints. The discovery of 8 novel non-canonical imprinted genes unique to the rat is consistent with more rapid evolution of extraembryonic imprinting. Meta-analysis of novel imprinted genes reveals multiple mechanisms by which species-specific imprinted expression may be established, including H3K27me3 deposition in the oocyte, the appearance of ZFP57 binding motifs, and the insertion of endogenous retroviral promoters.

CONCLUSIONS

In summary, we provide an expanded list of imprinted loci in the rat, reveal the extent of conservation of imprinted gene expression, and identify potential mechanisms responsible for the evolution of species-specific imprinting.

A long noncoding RNA influences the choice of the X chromosome to be inactivated

X chromosome inactivation (XCI) is the process of silencing one of the X chromosomes in cells of the female mammal which ensures dosage compensation between the sexes. Although theoretically random in somatic tissues, the choice of which X chromosome is chosen to be inactivated can be biased in mice by genetic element(s) associated with the so-called X-controlling element (Xce). Although the Xce was first described and genetically localized nearly 40 y ago, its mode of action remains elusive. In the approach presented here, we identify a single long noncoding RNA (lncRNA) within the Xce locus, Lppnx, which may be the driving factor in the choice of which X chromosome will be inactivated in the developing female mouse embryo. Comparing weak and strong Xce alleles we show that Lppnx modulates the expression of Xist lncRNA, one of the key factors in XCI, by controlling the occupancy of pluripotency factors at Intron1 of Xist. This effect is counteracted by enhanced binding of Rex1 in DxPas34, another key element in XCI regulating the activity of Tsix lncRNA, the main antagonist of Xist, in the strong but not in the weak Xce allele. These results suggest that the different susceptibility for XCI observed in weak and strong Xce alleles results from differential transcription factor binding of Xist Intron 1 and DxPas34, and that Lppnx represents a decisive factor in explaining the action of the Xce.

Simultaneous brain cell type and lineage determined by scRNA-seq reveals stereotyped cortical development

Mutations are acquired frequently, such that each cell's genome inscribes its history of cell divisions. Common genomic alterations involve loss of heterozygosity (LOH). LOH accumulates throughout the genome, offering large encoding capacity for inferring cell lineage. Using only single-cell RNA sequencing (scRNA-seq) of mouse brain cells, we found that LOH events spanning multiple genes are revealed as tracts of monoallelically expressed, constitutionally heterozygous single-nucleotide variants (SNVs). We simultaneously inferred cell lineage and marked developmental time points based on X chromosome inactivation and the total number of LOH events while identifying cell types from gene expression patterns. Our results are consistent with progenitor cells giving rise to multiple cortical cell types through stereotyped expansion and distinct waves of neurogenesis. This type of retrospective analysis could be incorporated into scRNA-seq pipelines and, compared with experimental approaches for determining lineage in model organisms, is applicable where genetic engineering is prohibited, such as humans.

Elastic dosage compensation by X-chromosome upregulation

X-chromosome inactivation and X-upregulation are the fundamental modes of chromosome-wide gene regulation that collectively achieve dosage compensation in mammals, but the regulatory link between the two remains elusive and the X-upregulation dynamics are unknown. Here, we use allele-resolved single-cell RNA-seq combined with chromatin accessibility profiling and finely dissect their separate effects on RNA levels during mouse development. Surprisingly, we uncover that X-upregulation elastically tunes expression dosage in a sex- and lineage-specific manner, and moreover along varying degrees of X-inactivation progression. Male blastomeres achieve X-upregulation upon zygotic genome activation while females experience two distinct waves of upregulation, upon imprinted and random X-inactivation; and ablation of Xist impedes female X-upregulation. Female cells carrying two active X chromosomes lack upregulation, yet their collective RNA output exceeds that of a single hyperactive allele. Importantly, this conflicts the conventional dosage compensation model in which naïve female cells are initially subject to biallelic X-upregulation followed by X-inactivation of one allele to correct the X dosage. Together, our study provides key insights to the chain of events of dosage compensation, explaining how transcript copy numbers can remain remarkably stable across developmental windows wherein severe dose imbalance would otherwise be experienced by the cell.

Mechanisms of Choice in X-Chromosome Inactivation

Early in development, placental and marsupial mammals harbouring at least two X chromosomes per nucleus are faced with a choice that affects the rest of their lives: which of those X chromosomes to transcriptionally inactivate. This choice underlies phenotypical diversity in the composition of tissues and organs and in their response to the environment, and can determine whether an individual will be healthy or affected by an X-linked disease. Here, we review our current understanding of the process of choice during X-chromosome inactivation and its implications, focusing on the strategies evolved by different mammalian lineages and on the known and unknown molecular mechanisms and players involved.

Diverse epigenetic mechanisms maintain parental imprints within the embryonic and extraembryonic lineages

Genomic imprinting and X chromosome inactivation (XCI) require epigenetic mechanisms to encode allele-specific expression, but how these specific tasks are accomplished at single loci or across chromosomal scales remains incompletely understood. Here, we systematically disrupt essential epigenetic pathways within polymorphic embryos in order to examine canonical and non-canonical genomic imprinting as well as XCI. We find that DNA methylation and Polycomb group repressors are indispensable for autosomal imprinting, albeit at distinct gene sets. Moreover, the extraembryonic ectoderm relies on a broader spectrum of imprinting mechanisms, including non-canonical targeting of maternal endogenous retrovirus (ERV)-driven promoters by the H3K9 methyltransferase G9a. We further identify Polycomb-dependent and -independent gene clusters on the imprinted X chromosome, which appear to reflect distinct domains of Xist-mediated suppression. From our data, we assemble a comprehensive inventory of the epigenetic pathways that maintain parent-specific imprinting in eutherian mammals, including an expanded view of the placental lineage.

Systematic Analysis of Monoallelic Gene Expression and Chromatin Accessibility Across Multiple Tissues in Hybrid Mice

In diploid eukaryotic organisms, both alleles of each autosomal gene are usually assumed to be simultaneously expressed at similar levels. However, some genes can be expressed preferentially or strictly from a single allele, a process known as monoallelic expression. Classic monoallelic expression of X-chromosome-linked genes, olfactory receptor genes and developmentally imprinted genes is the result of epigenetic modifications. Genetic-origin-dependent monoallelic expression, however, is caused by cis-regulatory differences between the alleles. There is a paucity of systematic study to investigate these phenomena across multiple tissues, and the mechanisms underlying such monoallelic expression are not yet fully understood. Here we provide a detailed portrait of monoallelic gene expression across multiple tissues/cell lines in a hybrid mouse cross between the Mus musculus strain C57BL/6J and the Mus spretus strain SPRET/EiJ. We observed pervasive tissue-dependent allele-specific gene expression: in total, 1,839 genes exhibited monoallelic expression in at least one tissue, and 410 genes in at least two tissues. Among these 88 are monoallelic genes with different active alleles between tissues, probably representing genetic-origin-dependent monoallelic expression. We also identified six autosomal monoallelic genes with the active allele being identical in all eight tissues, which are likely novel candidates of imprinted genes. To depict the underlying regulatory mechanisms at the chromatin layer, we performed ATAC-seq in two different cell lines derived from the F1 mouse. Consistent with the global expression pattern, cell-type dependent monoallelic peaks were found, and a higher proportion of C57BL/6J-active peaks were observed in both cell types, implying possible species-specific regulation. Finally, only a small part of monoallelic gene expression could be explained by allelic differences in chromatin organization in promoter regions, suggesting that other distal elements may play important roles in shaping the patterns of allelic gene expression across tissues.

Uncovering Modifier Genes of X-Linked Alport Syndrome Using a Novel Multiparent Mouse Model

BACKGROUND

Mutations in COL4A5 are responsible for 80% of cases of X-linked Alport Syndrome (XLAS). Although genes that cause AS are well characterized, people with AS who have similar genetic mutations present with a wide variation in the extent of kidney impairment and age of onset, suggesting the activities of modifier genes.

METHODS

We created a cohort of genetically diverse XLAS male and female mice using the Diversity Outbred mouse resource and measured albuminuria, GFR, and gene expression. Using a quantitative trait locus approach, we mapped modifier genes that can best explain the underlying phenotypic variation measured in our diverse population.

RESULTS

Genetic analysis identified several loci associated with the variation in albuminuria and GFR, including a locus on the X chromosome associated with X inactivation and a locus on chromosome 2 containing Fmn1. Subsequent analysis of genetically reduced Fmn1 expression in Col4a5 knockout mice showed a decrease in albuminuria, podocyte effacement, and podocyte protrusions in the glomerular basement membrane, which support the candidacy of Fmn1 as a modifier gene for AS.

CONCLUSION

With this novel approach, we emulated the variability in the severity of kidney phenotypes found in human patients with Alport Syndrome through albuminuria and GFR measurements. This approach can identify modifier genes in kidney disease that can be used as novel therapeutic targets.

Simpson-Golabi-Behmel syndrome: One family, same mutation, different outcome

Simpson-Golabi-Behmel syndrome (SGBS) is a rare X-linked condition characterized by pre and postnatal overgrowth with visceral and skeletal abnormalities. The syndrome is caused mainly by mutations in the X-linked gene GPC3. Clinical presentation of SGBS in affected males is well defined, but there is a lack of knowledge about affected females, with very few reported cases. In total, eight female carriers with clinical expression of SGBS have been reported to date. In the present report, we describe the ninth patient and her family history. The interesting features of our female patient are the Wilms' tumor and the transfontanelar ultrasound findings. The patient's older sister, carrier of the same mutation, has minor facial dysmorphisms but no congenital anomalies and so far, no further clinical findings, as well as her mother and grandmother. There is a lesson to be learned from these rare cases, namely that SGBS may have a significant clinical expression in females, and therefore, screening should be considered in all patients with SGBS regardless of the sex or phenotypic severity.

Bayesian modeling of skewed X inactivation in genetically diverse mice identifies a novel Xce allele associated with copy number changes

Female mammals are functional mosaics of their parental X-linked gene expression due to X chromosome inactivation (XCI). This process inactivates one copy of the X chromosome in each cell during embryogenesis and that state is maintained clonally through mitosis. In mice, the choice of which parental X chromosome remains active is determined by the X chromosome controlling element (Xce), which has been mapped to a 176-kb candidate interval. A series of functional Xce alleles has been characterized or inferred for classical inbred strains based on biased, or skewed, inactivation of the parental X chromosomes in crosses between strains. To further explore the function structure basis and location of the Xce, we measured allele-specific expression of X-linked genes in a large population of F1 females generated from Collaborative Cross (CC) strains. Using published sequence data and applying a Bayesian "Pólya urn" model of XCI skew, we report two major findings. First, inter-individual variability in XCI suggests mouse epiblasts contain on average 20-30 cells contributing to brain. Second, CC founder strain NOD/ShiLtJ has a novel and unique functional allele, Xceg, that is the weakest in the Xce allelic series. Despite phylogenetic analysis confirming that NOD/ShiLtJ carries a haplotype almost identical to the well-characterized C57BL/6J (Xceb), we observed unexpected patterns of XCI skewing in females carrying the NOD/ShiLtJ haplotype within the Xce. Copy number variation is common at the Xce locus and we conclude that the observed allelic series is a product of independent and recurring duplications shared between weak Xce alleles.

In vivo Firre and Dxz4 deletion elucidates roles for autosomal gene regulation

Recent evidence has determined that the conserved X chromosome mega-structures controlled by the Firre and Dxz4 loci are not required for X chromosome inactivation (XCI) in cell lines. Here, we examined the in vivo contribution of these loci by generating mice carrying a single or double deletion of Firre and Dxz4. We found that these mutants are viable, fertile and show no defect in random or imprinted XCI. However, the lack of these elements results in many dysregulated genes on autosomes in an organ-specific manner. By comparing the dysregulated genes between the single and double deletion, we identified superloop, megadomain, and Firre locus-dependent gene sets. The largest transcriptional effect was observed in all strains lacking the Firre locus, indicating that this locus is the main driver for these autosomal expression signatures. Collectively, these findings suggest that these X-linked loci are involved in autosomal gene regulation rather than XCI biology.

A Diallel of the Mouse Collaborative Cross Founders Reveals Strong Strain-Specific Maternal Effects on Litter Size

Reproductive success in the eight founder strains of the Collaborative Cross (CC) was measured using a diallel-mating scheme. Over a 48-month period we generated 4,448 litters, and provided 24,782 weaned pups for use in 16 different published experiments. We identified factors that affect the average litter size in a cross by estimating the overall contribution of parent-of-origin, heterosis, inbred, and epistatic effects using a Bayesian zero-truncated overdispersed Poisson mixed model. The phenotypic variance of litter size has a substantial contribution (82%) from unexplained and environmental sources, but no detectable effect of seasonality. Most of the explained variance was due to additive effects (9.2%) and parental sex (maternal vs. paternal strain; 5.8%), with epistasis accounting for 3.4%. Within the parental effects, the effect of the dam's strain explained more than the sire's strain (13.2% vs. 1.8%), and the dam's strain effects account for 74.2% of total variation explained. Dams from strains C57BL/6J and NOD/ShiLtJ increased the expected litter size by a mean of 1.66 and 1.79 pups, whereas dams from strains WSB/EiJ, PWK/PhJ, and CAST/EiJ reduced expected litter size by a mean of 1.51, 0.81, and 0.90 pups. Finally, there was no strong evidence for strain-specific effects on sex ratio distortion. Overall, these results demonstrate that strains vary substantially in their reproductive ability depending on their genetic background, and that litter size is largely determined by dam's strain rather than sire's strain effects, as expected. This analysis adds to our understanding of factors that influence litter size in mammals, and also helps to explain breeding successes and failures in the extinct lines and surviving CC strains.

Genome-wide stability of the DNA replication program in single mammalian cells

Here, we report a single-cell DNA replication sequencing method, scRepli-seq, a genome-wide methodology that measures copy number differences between replicated and unreplicated DNA. Using scRepli-seq, we demonstrate that replication-domain organization is conserved among individual mouse embryonic stem cells (mESCs). Differentiated mESCs exhibited distinct profiles, which were also conserved among cells. Haplotype-resolved scRepli-seq revealed similar replication profiles of homologous autosomes, while the inactive X chromosome was clearly replicated later than its active counterpart. However, a small degree of cell-to-cell replication-timing heterogeneity was present, which was smallest at the beginning and the end of S phase. In addition, developmentally regulated domains were found to deviate from others and showed a higher degree of heterogeneity, thus suggesting a link to developmental plasticity. Moreover, allelic expression imbalance was found to strongly associate with replication-timing asynchrony. Our results form a foundation for single-cell-level understanding of DNA replication regulation and provide insights into three-dimensional genome organization.

Short communication: Identification of the pseudoautosomal region in the Hereford bovine reference genome assembly ARS-UCD1.2

In cattle, the X chromosome accounts for approximately 3 and 6% of the genome in bulls and cows, respectively. In spite of the large size of this chromosome, very few studies report analysis of the X chromosome in genome-wide association studies and genomic selection. This lack of genetic interrogation is likely due to the complexities of undertaking these studies given the hemizygous state of some, but not all, of the X chromosome in males. The first step in facilitating analysis of this gene-rich chromosome is to accurately identify coordinates for the pseudoautosomal boundary (PAB) to split the chromosome into a region that may be treated as autosomal sequence (pseudoautosomal region) and a region that requires more complex statistical models. With the recent release of ARS-UCD1.2, a more complete and accurate assembly of the cattle genome than was previously available, it is timely to fine map the PAB for the first time. Here we report the use of SNP chip genotypes, short-read sequences, and long-read sequences to fine map the PAB (X chromosome:133,300,518) and simultaneously determine the neighboring regions of reduced homology and true pseudoautosomal region. These results greatly facilitate the inclusion of the X chromosome in genome-wide association studies, genomic selection, and other genetic analysis undertaken on this reference genome.

Epigenetic and Cellular Diversity in the Brain through Allele-Specific Effects

The benefits of diploidy are considered to involve masking partially recessive mutations and increasing genetic diversity. Here, we review new studies showing evidence for diverse allele-specific expression and epigenetic states in mammalian brain cells, which suggest that diploidy expands the landscape of gene regulatory and expression programs in cells. Allele-specific expression has been thought to be restricted to a few specific classes of genes. However, new studies show novel genomic imprinting effects that are brain-region-, cell-type- and age-dependent. In addition, novel forms of random monoallelic expression that impact many autosomal genes have been described in vitro and in vivo. We discuss the implications for understanding the benefits of diploidy, and the mechanisms shaping brain development, function, and disease.

Intragenomic Conflict and Immune Tolerance: Do Selfish X-Linked Alleles Drive Skewed X Chromosome Inactivation?

In mammalian females, diploid somatic cells contain two X chromosomes, one of which is transcriptionally silenced, in a process termed X chromosome inactivation (XCI). Whereas XCI is largely random in placental females, many women exhibit skewed XCI (SXCI), in which the vast majority cells have the same X chromosome inactivated. SXCI has serious health consequences, associated with conditions ranging from Alzheimer’s to various autoimmune disorders. SXCI is also associated with outcomes of pregnancies, with higher rates of recurrent spontaneous abortion in women with SXCI. Here, I suggest that SXCI could be driven by selfish X-linked alleles. Consistent with the association of SXCI with autoimmunity, I first note the possibility that recurrent spontaneous abortion could reflect immune rejection of fetuses inheriting alleles from the largely silenced maternal X chromosome. Preferential abortion of fetuses carrying silenced X-linked alleles implies a transmission advantage for X-linked alleles on the largely expressed chromosome, which could drive the emergence of X-linked alleles that make the chromosome resistant to XCI. I discuss the evolutionary dynamics, fitness tradeoffs and implications of this hypothesis, and suggest future directions.

The genomic ancestry, landscape genetics and invasion history of introduced mice in New Zealand

The house mouse (Mus musculus) provides a fascinating system for studying both the genomic basis of reproductive isolation, and the patterns of human-mediated dispersal. New Zealand has a complex history of mouse invasions, and the living descendants of these invaders have genetic ancestry from all three subspecies, although most are primarily descended from M. m. domesticus. We used the GigaMUGA genotyping array (approximately 135 000 loci) to describe the genomic ancestry of 161 mice, sampled from 34 locations from across New Zealand (and one Australian city-Sydney). Of these, two populations, one in the south of the South Island, and one on Chatham Island, showed complete mitochondrial lineage capture, featuring two different lineages of M. m. castaneus mitochondrial DNA but with only M. m. domesticus nuclear ancestry detectable. Mice in the northern and southern parts of the North Island had small traces (approx. 2-3%) of M. m. castaneus nuclear ancestry, and mice in the upper South Island had approximately 7-8% M. m. musculus nuclear ancestry including some Y-chromosomal ancestry-though no detectable M. m. musculus mitochondrial ancestry. This is the most thorough genomic study of introduced populations of house mice yet conducted, and will have relevance to studies of the isolation mechanisms separating subspecies of mice.

Mapping the mouse Allelome reveals tissue-specific regulation of allelic expression

To determine the dynamics of allelic-specific expression during mouse development, we analyzed RNA-seq data from 23 F1 tissues from different developmental stages, including 19 female tissues allowing X chromosome inactivation (XCI) escapers to also be detected. We demonstrate that allelic expression arising from genetic or epigenetic differences is highly tissue-specific. We find that tissue-specific strain-biased gene expression may be regulated by tissue-specific enhancers or by post-transcriptional differences in stability between the alleles. We also find that escape from X-inactivation is tissue-specific, with leg muscle showing an unexpectedly high rate of XCI escapers. By surveying a range of tissues during development, and performing extensive validation, we are able to provide a high confidence list of mouse imprinted genes including 18 novel genes. This shows that cluster size varies dynamically during development and can be substantially larger than previously thought, with the Igf2r cluster extending over 10 Mb in placenta.

DOI:http://dx.doi.org/10.7554/eLife.25125.001

Bovine mammary gland X chromosome inactivation

X chromosome inactivation (XCI) is a process by which 1 of the 2 copies of the X chromosomes present in female mammals is inactivated. The transcriptional silencing of one X chromosome achieves dosage compensation between XX females and XY males and ensures equal expression of X-linked genes in both sexes. Although all mammals use this form of dosage compensation, the complex mechanisms that regulate XCI vary between species, tissues, and development. These mechanisms include not only varying levels of inactivation, but also the nature of inactivation, which can range from being random in nature to driven by parent of origin. To date, no data describing XCI in calves or adult cattle have been reported and we are reliant on data from mice to infer potential mechanisms and timings for this process. In the context of dairy cattle breeding and genomic prediction, the implications of X chromosome inheritance and XCI in the mammary gland are particularly important where a relatively small number of bulls pass their single X chromosome on to all of their daughters. We describe here the use of RNA-seq, whole genome sequencing and Illumina BovineHD BeadChip (Illumina, San Diego, CA) genotypes to assess XCI in lactating mammary glands of dairy cattle. At a population level, maternally and paternally inherited copies of the X chromosome are expressed equally in the lactating mammary gland consistent with random inactivation of the X chromosome. However, average expression of the paternal chromosome ranged from 10 to 90% depending on the individual animal. These results suggest that either the mammary gland arises from 1 or 2 stem cells, or a nongenetic mechanism that skews XCI exists. Although a considerable amount of future work is required to fully understand XCI in cattle, the data reported here represent an initial step in ensuring that X chromosome variation is captured and used in an appropriate manner for future genomic selection.

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.

Developmental regulation of X-chromosome inactivation

With the emergence of sex-determination by sex chromosomes, which differ in composition and number between males and females, appeared the need to equalize X-chromosomal gene dosage between the sexes. Mammals have devised the strategy of X-chromosome inactivation (XCI), in which one of the two X-chromosomes is rendered transcriptionally silent in females. In the mouse, the best-studied model organism with respect to XCI, this inactivation process occurs in different forms, imprinted and random, interspersed by periods of X-chromosome reactivation (XCR), which is needed to switch between the different modes of XCI. In this review, I describe the recent advances with respect to the developmental control of XCI and XCR and in particular their link to differentiation and pluripotency. Furthermore, I review the mechanisms, which influence the timing and choice, with which one of the two X-chromosomes is chosen for inactivation during random XCI. This has an impact on how females are mosaics with regard to which X-chromosome is active in different cells, which has implications on the severity of diseases caused by X-linked mutations.

The Mouse Universal Genotyping Array: From Substrains to Subspecies

Genotyping microarrays are an important resource for genetic mapping, population genetics, and monitoring of the genetic integrity of laboratory stocks. We have developed the third generation of the Mouse Universal Genotyping Array (MUGA) series, GigaMUGA, a 143,259-probe Illumina Infinium II array for the house mouse (Mus musculus). The bulk of the content of GigaMUGA is optimized for genetic mapping in the Collaborative Cross and Diversity Outbred populations, and for substrain-level identification of laboratory mice. In addition to 141,090 single nucleotide polymorphism probes, GigaMUGA contains 2006 probes for copy number concentrated in structurally polymorphic regions of the mouse genome. The performance of the array is characterized in a set of 500 high-quality reference samples spanning laboratory inbred strains, recombinant inbred lines, outbred stocks, and wild-caught mice. GigaMUGA is highly informative across a wide range of genetically diverse samples, from laboratory substrains to other Mus species. In addition to describing the content and performance of the array, we provide detailed probe-level annotation and recommendations for quality control.

Noncanonical Genomic Imprinting Effects in Offspring

Here, we describe an RNA-sequencing (RNA-seq)-based approach that accurately detects even modest maternal or paternal allele expression biases at the tissue level, which we call noncanonical genomic imprinting effects. We profile imprinting in the arcuate nucleus (ARN) and dorsal raphe nucleus of the female mouse brain as well as skeletal muscle (mesodermal) and liver (endodermal). Our study uncovers hundreds of noncanonical autosomal and X-linked imprinting effects. Noncanonical imprinting is highly tissue-specific and enriched in the ARN, but rare in the liver. These effects are reproducible across different genetic backgrounds and associated with allele-specific chromatin. Using in situ hybridization for nascent RNAs, we discover that autosomal noncanonical imprinted genes with a tissue-level allele bias exhibit allele-specific expression effects in subpopulations of neurons in the brain in vivo. We define noncanonical imprinted genes that regulate monoamine signaling and determine that these effects influence the impact of inherited mutations on offspring behavior.

Allelome.PRO, a pipeline to define allele-specific genomic features from high-throughput sequencing data

Detecting allelic biases from high-throughput sequencing data requires an approach that maximises sensitivity while minimizing false positives. Here, we present Allelome.PRO, an automated user-friendly bioinformatics pipeline, which uses high-throughput sequencing data from reciprocal crosses of two genetically distinct mouse strains to detect allele-specific expression and chromatin modifications. Allelome.PRO extends approaches used in previous studies that exclusively analyzed imprinted expression to give a complete picture of the ‘allelome’ by automatically categorising the allelic expression of all genes in a given cell type into imprinted, strain-biased, biallelic or non-informative. Allelome.PRO offers increased sensitivity to analyze lowly expressed transcripts, together with a robust false discovery rate empirically calculated from variation in the sequencing data. We used RNA-seq data from mouse embryonic fibroblasts from F1 reciprocal crosses to determine a biologically relevant allelic ratio cutoff, and define for the first time an entire allelome. Furthermore, we show that Allelome.PRO detects differential enrichment of H3K4me3 over promoters from ChIP-seq data validating the RNA-seq results. This approach can be easily extended to analyze histone marks of active enhancers, or transcription factor binding sites and therefore provides a powerful tool to identify candidate cis regulatory elements genome wide.

Epigenetic events regulating monoallelic gene expression

In mammals, generally it is assumed that the genes inherited from each parent are expressed to similar levels. However, it is now apparent that in non-sex chromosomes, 6-10% of genes are selected for monoallelic expression. Monoallelic expression or allelic exclusion is established either in an imprinted (parent-of-origin) or a stochastic manner. The stochastic model explains random selection while the imprinted model describes parent-of-origin specific selection of alleles for expression. Allelic exclusion occurs during X chromosome inactivation, parent-of-origin expression of imprinted genes and stochastic monoallelic expression of cell surface molecules, clustered protocadherin (PCDH) genes. Mis-regulation or loss of allelic exclusion contributes to developmental diseases. Epigenetic mechanisms are fundamental players that determine this type of expression despite a homogenous genetic background. DNA methylation and histone modifications are two mediators of the epigenetic phenomena. The majority of DNA methylation is found on cytosines of the CpG dinucleotide in mammals. Several covalent modifications of histones change the electrostatic forces between DNA and histones modifying gene expression. Long-range chromatin interactions organize chromatin into transcriptionally permissive and prohibitive regions leading to simultaneous regulation of gene expression and repression. Non-coding RNAs (ncRNAs) are also players in regulating gene expression. Together, these epigenetic mechanisms fine-tune gene expression levels essential for normal development and survival. In this review, first we discuss what is known about monoallelic gene expression. Then, we focus on the molecular mechanisms that regulate expression of three monoallelically expressed gene classes: the X-linked genes, selected imprinted genes and PCDH genes.

X-chromosome inactivation: new insights into cis and trans regulation

X-chromosome inactivation (XCI) is a developmentally associated process that evolved in mammals to enable gene dosage compensation between XX and XY individuals. In placental mammals, it is triggered by the long noncoding RNA Xist, which is produced from a complex regulatory locus, the X-inactivation centre (Xic). Recent insights into the regulatory landscape of the Xic, including its partitioning into topological associating domains (TADs) and its genetic dissection, have important implications for the monoallelic regulation of Xist. Here, we present some of the latest studies on X inactivation with a special focus on the regulation of Xist, its various functions and the putative role of chromosome conformation in regulating the dynamics of this locus during development and differentiation.

Analyses of allele-specific gene expression in highly divergent mouse crosses identifies pervasive allelic imbalance

Complex human traits are influenced by variation in regulatory DNA through mechanisms that are not fully understood. Since regulatory elements are conserved between humans and mice, a thorough annotation of cis regulatory variants in mice could aid in this process. Here we provide a detailed portrait of mouse gene expression across multiple tissues in a three-way diallel. Greater than 80% of mouse genes have cis regulatory variation. These effects influence complex traits and usually extend to the human ortholog. Further, we estimate that at least one in every thousand SNPs creates a cis regulatory effect. We also observe two types of parent-of-origin effects, including classical imprinting and a novel, global allelic imbalance in favor of the paternal allele. We conclude that, as with humans, pervasive regulatory variation influences complex genetic traits in mice and provide a new resource toward understanding the genetic control of transcription in mammals.

A multi-megabase copy number gain causes maternal transmission ratio distortion on mouse chromosome 2

Significant departures from expected Mendelian inheritance ratios (transmission ratio distortion, TRD) are frequently observed in both experimental crosses and natural populations. TRD on mouse Chromosome (Chr) 2 has been reported in multiple experimental crosses, including the Collaborative Cross (CC). Among the eight CC founder inbred strains, we found that Chr 2 TRD was exclusive to females that were heterozygous for the WSB/EiJ allele within a 9.3 Mb region (Chr 2 76.9 - 86.2 Mb). A copy number gain of a 127 kb-long DNA segment (designated as responder to drive, R2d) emerged as the strongest candidate for the causative allele. We mapped R2d sequences to two loci within the candidate interval. R2d1 is located near the proximal boundary, and contains a single copy of R2d in all strains tested. R2d2 maps to a 900 kb interval, and the number of R2d copies varies from zero in classical strains (including the mouse reference genome) to more than 30 in wild-derived strains. Using real-time PCR assays for the copy number, we identified a mutation (R2d2WSBdel1) that eliminates the majority of the R2d2WSB copies without apparent alterations of the surrounding WSB/EiJ haplotype. In a three-generation pedigree segregating for R2d2WSBdel1, the mutation is transmitted to the progeny and Mendelian segregation is restored in females heterozygous for R2d2WSBdel1, thus providing direct evidence that the copy number gain is causal for maternal TRD. We found that transmission ratios in R2d2WSB heterozygous females vary between Mendelian segregation and complete distortion depending on the genetic background, and that TRD is under genetic control of unlinked distorter loci. Although the R2d2WSB transmission ratio was inversely correlated with average litter size, several independent lines of evidence support the contention that female meiotic drive is the cause of the distortion. We discuss the implications and potential applications of this novel meiotic drive system.

Known unknowns for allele-specific expression and genomic imprinting effects

Recent studies have provided evidence for non-canonical imprinting effects that are associated with allele-specific expression biases at the tissue level in mice. These imprinting effects have features that are distinct from canonical imprinting effects that involve allele silencing. Here, I discuss some of the evidence for non-canonical imprinting effects in the context of random X-inactivation and epigenetic allele-specific expression effects on the autosomes. I propose several mechanisms that may underlie non-canonical imprinting effects and outline future directions and approaches to study these effects at the cellular level in vivo. The growing evidence for complex allele-specific expression effects that are cell- and developmental stage-specific has opened a new frontier for study. Currently, the function of these effects and the underlying regulatory mechanisms are largely unknown.