A novel multi-alignment pipeline for high-throughput sequencing data

Phylogeography of introgression: Spatial and temporal analyses identify two introgression events between brown and American black bears

Brown bears (Ursus arctos) colonized North America from Eurasia in two distinct and temporally separated waves. Once in North America they encountered endemic American black bears (U. americanus) during range expansions from eastern Beringia southwards into the interior of the continent. The establishment of sympatry between these species provided the opportunity for hybridization and introgression, which was previously identified at the species level using D-statistics. Both species have broad spatial ranges that should limit the extent of introgression, such that it is found primarily between sympatric populations. Here, we used range-wide sampling and whole genome sequencing of both bear species to test for spatial variability in introgression. We identified two pulses of introgression between brown and American black bears, and demonstrate the introgressed segments occur across spatially structured lineages in both species. The first pulse occurred 270-120 kya, near the initiation of intraspecific divergence, approximately 99-93 kya, within each species. This pulse occurred as sympatry was established in western North America. The second pulse occurred between western American black bears and North American brown bears and lasted to 9 kya. Introgression was bidirectional and sympatric lineages had more introgressed tracts and a larger proportion of the genome introgressed from the other species. This study advances our phylogeographic understanding of both iconic bear species through investigating the timing of divergence and gene flow as bears expanded and contracted their ranges across North America.

Different complex regulatory phenotypes underlie hybrid male sterility in divergent rodent crosses

Hybrid incompatibilities are a critical component of species barriers and may arise due to negative interactions between divergent regulatory elements in parental species. We used a comparative approach to identify common themes in the regulatory phenotypes associated with hybrid male sterility in two divergent rodent crosses, dwarf hamsters and house mice. We investigated three potential characteristic gene expression phenotypes in hybrids including the propensity of transgressive differentially expressed genes toward over or underexpression, the influence of developmental stage on patterns of misexpression, and the role of the sex chromosomes on misexpression phenotypes. In contrast to near pervasive overexpression in hybrid house mice, we found that misexpression in hybrid dwarf hamsters was dependent on developmental stage. In both house mouse and dwarf hamster hybrids, however, misexpression increased with the progression of spermatogenesis, although to varying extents and with potentially different consequences. In both systems, we detected sex chromosome-specific overexpression in stages of spermatogenesis where inactivated X chromosome expression was expected, but the hybrid overexpression phenotypes were fundamentally different. Importantly, misexpression phenotypes support the presence of multiple developmental blocks to spermatogenesis in dwarf hamster hybrids, including a potential role of meiotic stalling or breakdown early in spermatogenesis. Collectively, we demonstrate that while there are some similarities in hybrid expression phenotypes of house mice and dwarf hamsters, there are also clear differences that point toward unique mechanisms underlying hybrid male sterility. Our results highlight the potential of comparative approaches in helping to understand the causes and consequences of disrupted gene expression in speciation.

Different complex regulatory phenotypes underlie hybrid male sterility in divergent rodent crosses

Hybrid incompatibilities are a critical component of species barriers and may arise due to negative interactions between divergent regulatory elements in parental species. We used a comparative approach to identify common themes in the regulatory phenotypes associated with hybrid male sterility in two divergent rodent crosses, dwarf hamsters and house mice. We investigated three potential characteristic gene expression phenotypes in hybrids including the propensity of transgressive differentially expressed genes towards over or underexpression, the influence of developmental stage on patterns of misexpression, and the role of the sex chromosomes on misexpression phenotypes. In contrast to near pervasive overexpression in hybrid house mice, we found that misexpression in hybrid dwarf hamsters was dependent on developmental stage. In both house mouse and dwarf hamster hybrids, however, misexpression increased with the progression of spermatogenesis, although to varying extents and with potentially different consequences. In both systems, we detected sex-chromosome specific overexpression in stages of spermatogenesis where inactivated X chromosome expression was expected, but the hybrid overexpression phenotypes were fundamentally different. Importantly, misexpression phenotypes support the presence of multiple developmental blocks to spermatogenesis in dwarf hamster hybrids, including a potential role of meiotic stalling or breakdown early in spermatogenesis. Collectively, we demonstrate that while there are some similarities in hybrid expression phenotypes of house mice and dwarf hamsters, there are also clear differences that point towards unique mechanisms underlying hybrid male sterility. Our results highlight the potential of comparative approaches in helping to understand the causes and consequences of disrupted gene expression in speciation.

Biparental graph strategy to represent and analyze hybrid plant genomes

Hybrid plants are found extensively in the wild, and they often demonstrate superior performance of complex traits over their parents and other selfing plants. This phenomenon, known as heterosis, has been extensively applied in plant breeding for decades. However, the process of decoding hybrid plant genomes has seriously lagged due to the challenges associated with genome assembly and the lack of appropriate methodologies for their subsequent representation and analysis. Here, we present the assembly and analysis of 2 hybrids, an intraspecific hybrid between 2 maize (Zea mays ssp. mays) inbred lines and an interspecific hybrid between maize and its wild relative teosinte (Z. mays ssp. parviglumis), utilizing a combination of PacBio High Fidelity sequencing and chromatin conformation capture sequencing data. The haplotypic assemblies are well phased at chromosomal scale, successfully resolving the complex loci with extensive parental structural variations (SVs). By integrating into a biparental genome graph, the haplotypic assemblies can facilitate downstream short-read-based SV calling and allele-specific gene expression analysis, demonstrating outstanding advantages over a single linear genome. Our work offers a comprehensive workflow that aims to facilitate the decoding of numerous hybrid plant genomes, particularly those with unknown or inaccessible parentage, thereby enhancing our understanding of genome evolution and heterosis.

An integrated transcription factor framework for Treg identity and diversity

Vertebrate cell identity depends on the combined activity of scores of transcription factors (TF). While TFs have often been studied in isolation, a systematic perspective on their integration has been missing. Focusing on FoxP3+ regulatory T cells (Tregs), key guardians of immune tolerance, we combined single-cell chromatin accessibility, machine learning, and high-density genetic variation, to resolve a validated framework of diverse Treg chromatin programs, each shaped by multi-TF inputs. This framework identified previously unrecognized Treg controllers (Smarcc1) and illuminated the mechanism of action of FoxP3, which amplified a pre-existing Treg identity, diversely activating or repressing distinct programs, dependent on different regulatory partners. Treg subpopulations in the colon relied variably on FoxP3, Helios+ Tregs being completely dependent, but RORγ+ Tregs largely independent. These differences were rooted in intrinsic biases decoded by the integrated framework. Moving beyond master regulators, this work unravels how overlapping TF activities coalesce into Treg identity and diversity.

Reference Genome Choice and Filtering Thresholds Jointly Influence Phylogenomic Analyses

Molecular phylogenies are a cornerstone of modern comparative biology and are commonly employed to investigate a range of biological phenomena, such as diversification rates, patterns in trait evolution, biogeography, and community assembly. Recent work has demonstrated that significant biases may be introduced into downstream phylogenetic analyses from processing genomic data; however, it remains unclear whether there are interactions among bioinformatic parameters or biases introduced through the choice of reference genome for sequence alignment and variant calling. We address these knowledge gaps by employing a combination of simulated and empirical data sets to investigate the extent to which the choice of reference genome in upstream bioinformatic processing of genomic data influences phylogenetic inference, as well as the way that reference genome choice interacts with bioinformatic filtering choices and phylogenetic inference method. We demonstrate that more stringent minor allele filters bias inferred trees away from the true species tree topology, and that these biased trees tend to be more imbalanced and have a higher center of gravity than the true trees. We find the greatest topological accuracy when filtering sites for minor allele count (MAC) >3-4 in our 51-taxa data sets, while tree center of gravity was closest to the true value when filtering for sites with MAC >1-2. In contrast, filtering for missing data increased accuracy in the inferred topologies; however, this effect was small in comparison to the effect of minor allele filters and may be undesirable due to a subsequent mutation spectrum distortion. The bias introduced by these filters differs based on the reference genome used in short read alignment, providing further support that choosing a reference genome for alignment is an important bioinformatic decision with implications for downstream analyses. These results demonstrate that attributes of the study system and dataset (and their interaction) add important nuance for how best to assemble and filter short-read genomic data for phylogenetic inference.

AIRE relies on Z-DNA to flag gene targets for thymic T cell tolerization

AIRE is an unconventional transcription factor that enhances the expression of thousands of genes in medullary thymic epithelial cells and promotes clonal deletion or phenotypic diversion of self-reactive T cells1-4. The biological logic of AIRE's target specificity remains largely unclear as, in contrast to many transcription factors, it does not bind to a particular DNA sequence motif. Here we implemented two orthogonal approaches to investigate AIRE's cis-regulatory mechanisms: construction of a convolutional neural network and leveraging natural genetic variation through analysis of F1 hybrid mice5. Both approaches nominated Z-DNA and NFE2-MAF as putative positive influences on AIRE's target choices. Genome-wide mapping studies revealed that Z-DNA-forming and NFE2L2-binding motifs were positively associated with the inherent ability of a gene's promoter to generate DNA double-stranded breaks, and promoters showing strong double-stranded break generation were more likely to enter a poised state with accessible chromatin and already-assembled transcriptional machinery. Consequently, AIRE preferentially targets genes with poised promoters. We propose a model in which Z-DNA anchors the AIRE-mediated transcriptional program by enhancing double-stranded break generation and promoter poising. Beyond resolving a long-standing mechanistic conundrum, these findings suggest routes for manipulating T cell tolerance.

Evolution of parent-of-origin effects on placental gene expression in house mice

The mammalian placenta is a hotspot for the evolution of genomic imprinting, a form of gene regulation that involves the parent-specific epigenetic silencing of one allele. Imprinted genes are central to placental development and are thought to contribute to the evolution of reproductive barriers between species. However, it is unclear how rapidly imprinting evolves or how functional specialization among placental tissues influences the evolution of imprinted expression. We compared parent-of-origin expression bias across functionally distinct placental layers sampled from reciprocal crosses within three closely related lineages of mice ( Mus ). Using genome-wide gene expression and DNA methylation data from fetal and maternal tissues, we developed an analytical strategy to minimize pervasive bias introduced by maternal contamination of placenta samples. We corroborated imprinted expression at 42 known imprinted genes and identified five candidate imprinted genes showing parent-of-origin specific expression and DNA methylation. Paternally-biased expression was enriched in the labyrinth zone, a layer specialized in nutrient transfer, and maternally-biased genes were enriched in the junctional zone, which specializes in modulation of maternal physiology. Differentially methylated regions were predominantly determined through epigenetic modification of the maternal genome and were associated with both maternally- and paternally-biased gene expression. Lastly, comparisons between lineages revealed a small set of co-regulated genes showing rapid divergence in expression levels and imprinted status in the M. m. domesticus lineage. Together, our results reveal important links between core functional elements of placental biology and the evolution of imprinted gene expression among closely related rodent species.

Cystathionine β-synthase as novel endogenous regulator of lymphangiogenesis via modulating VEGF receptor 2 and 3

Lymphangiogenesis is a key player in several diseases such as tumor metastasis, obesity, and graft rejection. Endogenous regulation of lymphangiogenesis is only partly understood. Here we use the normally avascular cornea as a model to identify endogenous regulators of lymphangiogenesis. Quantitative trait locus analysis of a large low-lymphangiogenic BALB/cN x high-lymphangiogenic C57BL/6 N intercross and prioritization by whole-transcriptome sequencing identify a novel gene responsible for differences in lymphatic vessel architecture on chromosome 17, the cystathionine β-synthase (Cbs). Inhibition of CBS in lymphatic endothelial cells results in reduce proliferation, migration, altered tube-formation, and decrease expression of vascular endothelial growth factor (VEGF) receptor 2 (VEGF-R2) and VEGF-R3, but not their ligands VEGF-C and VEGF-D. Also in vivo inflammation-induced lymphangiogenesis is significantly reduce in C57BL/6 N mice after pharmacological inhibition of CBS. The results confirm CBS as a novel endogenous regulator of lymphangiogenesis acting via VEGF receptor 2 and 3-regulation and open new treatment avenues in diseases associated with pathologic lymphangiogenesis.

Characterization of sequence determinants of enhancer function using natural genetic variation

Sequence variation in enhancers that control cell-type-specific gene transcription contributes significantly to phenotypic variation within human populations. However, it remains difficult to predict precisely the effect of any given sequence variant on enhancer function due to the complexity of DNA sequence motifs that determine transcription factor (TF) binding to enhancers in their native genomic context. Using F₁-hybrid cells derived from crosses between distantly related inbred strains of mice, we identified thousands of enhancers with allele-specific TF binding and/or activity. We find that genetic variants located within the central region of enhancers are most likely to alter TF binding and enhancer activity. We observe that the AP-1 family of TFs (Fos/Jun) are frequently required for binding of TEAD TFs and for enhancer function. However, many sequence variants outside of core motifs for AP-1 and TEAD also impact enhancer function, including sequences flanking core TF motifs and AP-1 half sites. Taken together, these data represent one of the most comprehensive assessments of allele-specific TF binding and enhancer function to date and reveal how sequence changes at enhancers alter their function across evolutionary timescales.

Stage-specific disruption of X chromosome expression during spermatogenesis in sterile house mouse hybrids

Hybrid sterility is a complex phenotype that can result from the breakdown of spermatogenesis at multiple developmental stages. Here, we disentangle two proposed hybrid male sterility mechanisms in the house mice, Mus musculus domesticus and M. m. musculus, by comparing patterns of gene expression in sterile F1 hybrids from a reciprocal cross. We found that hybrid males from both cross directions showed disrupted X chromosome expression during prophase of meiosis I consistent with a loss of meiotic sex chromosome inactivation (MSCI) and Prdm9-associated sterility, but that the degree of disruption was greater in mice with an M. m. musculus X chromosome consistent with previous studies. During postmeiotic development, gene expression on the X chromosome was only disrupted in one cross direction, suggesting that misexpression at this later stage was genotype-specific and not a simple downstream consequence of MSCI disruption which was observed in both reciprocal crosses. Instead, disrupted postmeiotic expression may depend on the magnitude of earlier disrupted MSCI, or the disruption of particular X-linked genes or gene networks. Alternatively, only hybrids with a potential deficit of Sly copies, a Y-linked ampliconic gene family, showed overexpression in postmeiotic cells, consistent with a previously proposed model of antagonistic coevolution between the X- and Y-linked ampliconic genes contributing to disrupted expression late in spermatogenesis. The relative contributions of these two regulatory mechanisms and their impact on sterility phenotypes await further study. Our results further support the hypothesis that X-linked hybrid sterility in house mice has a variable genetic basis, and that genotype-specific disruption of gene regulation contributes to overexpression of the X chromosome at different stages of development. Overall, these findings underscore the critical role of epigenetic regulation of the X chromosome during spermatogenesis and suggest that these processes are prone to disruption in hybrids.

Unraveling patterns of disrupted gene expression across a complex tissue

Whole tissue RNASeq is the standard approach for studying gene expression divergence in evolutionary biology and provides a snapshot of the comprehensive transcriptome for a given tissue. However, whole tissues consist of diverse cell types differing in expression profiles, and the cellular composition of these tissues can evolve across species. Here, we investigate the effects of different cellular composition on whole tissue expression profiles. We compared gene expression from whole testes and enriched spermatogenesis populations in two species of house mice, Mus musculus musculus and M. m. domesticus, and their sterile and fertile F1 hybrids, which differ in both cellular composition and regulatory dynamics. We found that cellular composition differences skewed expression profiles and differential gene expression in whole testes samples. Importantly, both approaches were able to detect large-scale patterns such as disrupted X chromosome expression, although whole testes sampling resulted in decreased power to detect differentially expressed genes. We encourage researchers to account for histology in RNASeq and consider methods that reduce sample complexity whenever feasible. Ultimately, we show that differences in cellular composition between tissues can modify expression profiles, potentially altering inferred gene ontological processes, insights into gene network evolution, and processes governing gene expression evolution.

Molecular Evolution across Mouse Spermatogenesis

Genes involved in spermatogenesis tend to evolve rapidly, but we lack a clear understanding of how protein sequences and patterns of gene expression evolve across this complex developmental process. We used fluorescence-activated cell sorting (FACS) to generate expression data for early (meiotic) and late (postmeiotic) cell types across 13 inbred strains of mice (Mus) spanning ∼7 My of evolution. We used these comparative developmental data to investigate the evolution of lineage-specific expression, protein-coding sequences, and expression levels. We found increased lineage specificity and more rapid protein-coding and expression divergence during late spermatogenesis, suggesting that signatures of rapid testis molecular evolution are punctuated across sperm development. Despite strong overall developmental parallels in these components of molecular evolution, protein and expression divergences were only weakly correlated across genes. We detected more rapid protein evolution on the X chromosome relative to the autosomes, whereas X-linked gene expression tended to be relatively more conserved likely reflecting chromosome-specific regulatory constraints. Using allele-specific FACS expression data from crosses between four strains, we found that the relative contributions of different regulatory mechanisms also differed between cell types. Genes showing cis-regulatory changes were more common late in spermatogenesis, and tended to be associated with larger differences in expression levels and greater expression divergence between species. In contrast, genes with trans-acting changes were more common early and tended to be more conserved across species. Our findings advance understanding of gene evolution across spermatogenesis and underscore the fundamental importance of developmental context in molecular evolutionary studies.

The Transcription Factor Foxp3 Shapes Regulatory T Cell Identity by Tuning the Activity of trans-Acting Intermediaries

Regulatory T (Treg) cell identity is defined by the lineage-specifying transcription factor (TF) Foxp3. Here we examined mechanisms of Foxp3 function by leveraging naturally occurring genetic variation in wild-derived inbred mice, which enables the identification of DNA sequence motifs driving epigenetic features. Chromatin accessibility, TF binding, and gene expression patterns in resting and activated subsets of Treg cells, conventional CD4 T cells, and cells expressing a Foxp3 reporter null allele revealed that the majority of Foxp3-dependent changes occurred at sites not bound by Foxp3. Chromatin accessibility of these indirect Foxp3 targets depended on the presence of DNA binding motifs for other TFs, including TCF1. Foxp3 expression correlated with decreased TCF1 and reduced accessibility of TCF1-bound chromatin regions. Deleting one copy of the Tcf7 gene recapitulated Foxp3-dependent negative regulation of chromatin accessibility. Thus, Foxp3 defines Treg cell identity in a largely indirect manner by fine-tuning the activity of other major chromatin remodeling TFs such as TCF1.

Mutation of NEKL-4/NEK10 and TTLL genes suppress neuronal ciliary degeneration caused by loss of CCPP-1 deglutamylase function

Ciliary microtubules are subject to post-translational modifications that act as a "Tubulin Code" to regulate motor traffic, binding proteins and stability. In humans, loss of CCP1, a cytosolic carboxypeptidase and tubulin deglutamylating enzyme, causes infantile-onset neurodegeneration. In C. elegans, mutations in ccpp-1, the homolog of CCP1, result in progressive degeneration of neuronal cilia and loss of neuronal function. To identify genes that regulate microtubule glutamylation and ciliary integrity, we performed a forward genetic screen for suppressors of ciliary degeneration in ccpp-1 mutants. We isolated the ttll-5(my38) suppressor, a mutation in a tubulin tyrosine ligase-like glutamylase gene. We show that mutation in the ttll-4, ttll-5, or ttll-11 gene suppressed the hyperglutamylation-induced loss of ciliary dye filling and kinesin-2 mislocalization in ccpp-1 cilia. We also identified the nekl-4(my31) suppressor, an allele affecting the NIMA (Never in Mitosis A)-related kinase NEKL-4/NEK10. In humans, NEK10 mutation causes bronchiectasis, an airway and mucociliary transport disorder caused by defective motile cilia. C. elegans NEKL-4 localizes to the ciliary base but does not localize to cilia, suggesting an indirect role in ciliary processes. This work defines a pathway in which glutamylation, a component of the Tubulin Code, is written by TTLL-4, TTLL-5, and TTLL-11; is erased by CCPP-1; is read by ciliary kinesins; and its downstream effects are modulated by NEKL-4 activity. Identification of regulators of microtubule glutamylation in diverse cellular contexts is important to the development of effective therapies for disorders characterized by changes in microtubule glutamylation. By identifying C. elegans genes important for neuronal and ciliary stability, our work may inform research into the roles of the tubulin code in human ciliopathies and neurodegenerative diseases.

KSP: an integrated method for predicting catalyzing kinases of phosphorylation sites in proteins

BACKGROUND

Protein phosphorylation by kinases plays crucial roles in various biological processes including signal transduction and tumorigenesis, thus a better understanding of protein phosphorylation events in cells is fundamental for studying protein functions and designing drugs to treat diseases caused by the malfunction of phosphorylation. Although a large number of phosphorylation sites in proteins have been identified using high-throughput phosphoproteomic technologies, their specific catalyzing kinases remain largely unknown. Therefore, computational methods are urgently needed to predict the kinases that catalyze the phosphorylation of these sites.

RESULTS

We developed KSP, a new algorithm for predicting catalyzing kinases for experimentally identified phosphorylation sites in human proteins. KSP constructs a network based on known protein-protein interactions and kinase-substrate relationships. Based on the network, it computes an affinity score between a phosphorylation site and kinases, and returns the top-ranked kinases of the score as candidate catalyzing kinases. When tested on known kinase-substrate pairs, KSP outperforms existing methods including NetworKIN, iGPS, and PKIS.

CONCLUSIONS

We developed a novel accurate tool for predicting catalyzing kinases of known phosphorylation sites. It can work as a complementary network approach for sequence-based phosphorylation site predictors.

An Activity-Mediated Transition in Transcription in Early Postnatal Neurons

The maturation of the mammalian brain occurs after birth, and this stage of neuronal development is frequently impaired in neurological disorders, such as autism and schizophrenia. However, the mechanisms that regulate postnatal brain maturation are poorly defined. By purifying neuronal subpopulations across brain development in mice, we identify a postnatal switch in the transcriptional regulatory circuits that operates in the maturing mammalian brain. We show that this developmental transition includes the formation of hundreds of cell-type-specific neuronal enhancers that appear to be modulated by neuronal activity. Once selected, these enhancers are active throughout adulthood, suggesting that their formation in early life shapes neuronal identity and regulates mature brain function.

Landscape of genomic imprinting and its functions in the mouse mammary gland

Genomic imprinting is an epigenetic modification of DNA, whereby gene expression is restricted to either maternally or paternally inherited alleles. Imprinted genes (IGs) in the placenta and embryo are essential for growth regulation and nutrient supply. However, despite being an important nutrition delivery organ, studies on mammary gland genomic imprinting remain limited. In this study, we found that both the number of IGs and their expression levels decreased during development of the mouse mammary gland. IG expression was lineage-specific and related to mammary gland development and lactation. Meta-analysis of single-cell RNA sequencing data revealed that mammary gland IGs were co-expressed in a network that regulated cell stemness and differentiation, which was confirmed by our functional studies. Accordingly, our data indicated that IGs were essential for the self-renewal of mammary gland stem cells and IG decline was correlated with mammary gland maturity. Taken together, our findings revealed the importance of IGs in a poorly studied nutrition-related organ, i.e. the mammary gland, thus providing a reference for further studies on genomic imprinting.

Developmental Analysis of Mimulus Seed Transcriptomes Reveals Functional Gene Expression Clusters and Four Imprinted, Endosperm-Expressed Genes

The double fertilization of the female gametophyte initiates embryogenesis and endosperm development in seeds via the activation of genes involved in cell differentiation, organ patterning, and growth. A subset of genes expressed in endosperm exhibit imprinted expression, and the correct balance of gene expression between parental alleles is critical for proper endosperm and seed development. We use a transcriptional time series analysis to identify genes that are associated with key shifts in seed development, including genes associated with secondary cell wall synthesis, mitotic cell cycle, chromatin organization, auxin synthesis, fatty acid metabolism, and seed maturation. We relate these genes to morphological changes in Mimulus seeds. We also identify four endosperm-expressed transcripts that display imprinted (paternal) expression bias. The imprinted status of these four genes is conserved in other flowering plants, suggesting that they are functionally important in endosperm development. Our study explores gene regulatory dynamics in a species with ab initio cellular endosperm development, broadening the taxonomic focus of the literature on gene expression in seeds. Moreover, it is the first to validate genes with imprinted endosperm expression in Mimulus guttatus, and will inform future studies on the genetic causes of seed failure in this model system.

Novel Neuroprotective Loci Modulating Ischemic Stroke Volume in Wild-Derived Inbred Mouse Strains

To identify genes involved in cerebral infarction, we have employed a forward genetic approach in inbred mouse strains, using quantitative trait loci (QTL) mapping for cerebral infarct volume after middle cerebral artery occlusion. We had previously observed that infarct volume is inversely correlated with cerebral collateral vessel density in most strains. In this study, we expanded the pool of allelic variation among classical inbred mouse strains by utilizing the eight founder strains of the Collaborative Cross and found a wild-derived strain, WSB/EiJ, that breaks this general rule that collateral vessel density inversely correlates with infarct volume. WSB/EiJ and another wild-derived strain, CAST/EiJ, show the highest collateral vessel densities of any inbred strain, but infarct volume of WSB/EiJ mice is 8.7-fold larger than that of CAST/EiJ mice. QTL mapping between these strains identified four new neuroprotective loci modulating cerebral infarct volume while not affecting collateral vessel phenotypes. To identify causative variants in genes, we surveyed nonsynonymous coding SNPs between CAST/EiJ and WSB/EiJ and found 96 genes harboring coding SNPs predicted to be damaging and mapping within one of the four intervals. In addition, we performed RNA-sequencing for brain tissue of CAST/EiJ and WSB/EiJ mice and identified 79 candidate genes mapping in one of the four intervals showing strain-specific differences in expression. The identification of the genes underlying these neuroprotective loci will provide new understanding of genetic risk factors of ischemic stroke, which may provide novel targets for future therapeutic intervention of human ischemic stroke.

Natural Genetic Variation Reveals Key Features of Epigenetic and Transcriptional Memory in Virus-Specific CD8 T Cells

Stable changes in chromatin states and gene expression in cells of the immune system form the basis for memory of infections and other challenges. Here, we used naturally occurring cis-regulatory variation in wild-derived inbred mouse strains to explore the mechanisms underlying long-lasting versus transient gene regulation in CD8 T cells responding to acute viral infection. Stably responsive DNA elements were characterized by dramatic and congruent chromatin remodeling events affecting multiple neighboring sites and required distinct transcription factor (TF) binding motifs for their accessibility. Specifically, we found that cooperative recruitment of T-box and Runx family transcription factors to shared targets mediated stable chromatin remodeling upon T cell activation. Our observations provide insights into the molecular mechanisms driving virus-specific CD8 T cell responses and suggest a general mechanism for the formation of transcriptional and epigenetic memory applicable to other immune and non-immune cells.

Investigation and identification of functional post-translational modification sites associated with drug binding and protein-protein interactions

Background

Protein post-translational modification (PTM) plays an essential role in various cellular processes that modulates the physical and chemical properties, folding, conformation, stability and activity of proteins, thereby modifying the functions of proteins. The improved throughput of mass spectrometry (MS) or MS/MS technology has not only brought about a surge in proteome-scale studies, but also contributed to a fruitful list of identified PTMs. However, with the increase in the number of identified PTMs, perhaps the more crucial question is what kind of biological mechanisms these PTMs are involved in. This is particularly important in light of the fact that most protein-based pharmaceuticals deliver their therapeutic effects through some form of PTM. Yet, our understanding is still limited with respect to the local effects and frequency of PTM sites near pharmaceutical binding sites and the interfaces of protein-protein interaction (PPI). Understanding PTM’s function is critical to our ability to manipulate the biological mechanisms of protein.

Results

In this study, to understand the regulation of protein functions by PTMs, we mapped 25,835 PTM sites to proteins with available three-dimensional (3D) structural information in the Protein Data Bank (PDB), including 1785 modified PTM sites on the 3D structure. Based on the acquired structural PTM sites, we proposed to use five properties for the structural characterization of PTM substrate sites: the spatial composition of amino acids, residues and side-chain orientations surrounding the PTM substrate sites, as well as the secondary structure, division of acidity and alkaline residues, and solvent-accessible surface area. We further mapped the structural PTM sites to the structures of drug binding and PPI sites, identifying a total of 1917 PTM sites that may affect PPI and 3951 PTM sites associated with drug-target binding. An integrated analytical platform (CruxPTM), with a variety of methods and online molecular docking tools for exploring the structural characteristics of PTMs, is presented. In addition, all tertiary structures of PTM sites on proteins can be visualized using the JSmol program.

Conclusion

Resolving the function of PTM sites is important for understanding the role that proteins play in biological mechanisms. Our work attempted to delineate the structural correlation between PTM sites and PPI or drug-target binding. CurxPTM could help scientists narrow the scope of their PTM research and enhance the efficiency of PTM identification in the face of big proteome data. CruxPTM is now available at http://csb.cse.yzu.edu.tw/CruxPTM/.

Electronic supplementary material

The online version of this article (10.1186/s12918-017-0506-1) contains supplementary material, which is available to authorized users.

Genomes of the Mouse Collaborative Cross

The Collaborative Cross (CC) is a multiparent panel of recombinant inbred (RI) mouse strains derived from eight founder laboratory strains. RI panels are popular because of their long-term genetic stability, which enhances reproducibility and integration of data collected across time and conditions. Characterization of their genomes can be a community effort, reducing the burden on individual users. Here we present the genomes of the CC strains using two complementary approaches as a resource to improve power and interpretation of genetic experiments. Our study also provides a cautionary tale regarding the limitations imposed by such basic biological processes as mutation and selection. A distinct advantage of inbred panels is that genotyping only needs to be performed on the panel, not on each individual mouse. The initial CC genome data were haplotype reconstructions based on dense genotyping of the most recent common ancestors (MRCAs) of each strain followed by imputation from the genome sequence of the corresponding founder inbred strain. The MRCA resource captured segregating regions in strains that were not fully inbred, but it had limited resolution in the transition regions between founder haplotypes, and there was uncertainty about founder assignment in regions of limited diversity. Here we report the whole genome sequence of 69 CC strains generated by paired-end short reads at 30× coverage of a single male per strain. Sequencing leads to a substantial improvement in the fine structure and completeness of the genomes of the CC. Both MRCAs and sequenced samples show a significant reduction in the genome-wide haplotype frequencies from two wild-derived strains, CAST/EiJ and PWK/PhJ. In addition, analysis of the evolution of the patterns of heterozygosity indicates that selection against three wild-derived founder strains played a significant role in shaping the genomes of the CC. The sequencing resource provides the first description of tens of thousands of new genetic variants introduced by mutation and drift in the CC genomes. We estimate that new SNP mutations are accumulating in each CC strain at a rate of 2.4 ± 0.4 per gigabase per generation. The fixation of new mutations by genetic drift has introduced thousands of new variants into the CC strains. The majority of these mutations are novel compared to currently sequenced laboratory stocks and wild mice, and some are predicted to alter gene function. Approximately one-third of the CC inbred strains have acquired large deletions (>10 kb) many of which overlap known coding genes and functional elements. The sequence of these mice is a critical resource to CC users, increases threefold the number of mouse inbred strain genomes available publicly, and provides insight into the effect of mutation and drift on common resources.

Phylogenomic Insights into Mouse Evolution Using a Pseudoreference Approach

Comparative genomic studies are now possible across a broad range of evolutionary timescales, but the generation and analysis of genomic data across many different species still present a number of challenges. The most sophisticated genotyping and down-stream analytical frameworks are still predominantly based on comparisons to high-quality reference genomes. However, established genomic resources are often limited within a given group of species, necessitating comparisons to divergent reference genomes that could restrict or bias comparisons across a phylogenetic sample. Here, we develop a scalable pseudoreference approach to iteratively incorporate sample-specific variation into a genome reference and reduce the effects of systematic mapping bias in downstream analyses. To characterize this framework, we used targeted capture to sequence whole exomes (∼54 Mbp) in 12 lineages (ten species) of mice spanning the Mus radiation. We generated whole exome pseudoreferences for all species and show that this iterative reference-based approach improved basic genomic analyses that depend on mapping accuracy while preserving the associated annotations of the mouse reference genome. We then use these pseudoreferences to resolve evolutionary relationships among these lineages while accounting for phylogenetic discordance across the genome, contributing an important resource for comparative studies in the mouse system. We also describe patterns of genomic introgression among lineages and compare our results to previous studies. Our general approach can be applied to whole or partitioned genomic data and is easily portable to any system with sufficient genomic resources, providing a useful framework for phylogenomic studies in mice and other taxa.

The evolutionary capacitor HSP90 buffers the regulatory effects of mammalian endogenous retroviruses

Understanding how genotypes are linked to phenotypes is important in biomedical and evolutionary studies. The chaperone heat-shock protein 90 (HSP90) buffers genetic variation by stabilizing proteins with variant sequences, thereby uncoupling phenotypes from genotypes. Here we report an unexpected role of HSP90 in buffering cis-regulatory variation affecting gene expression. By using the tripartite-motif-containing 28 (TRIM28; also known as KAP1)-mediated epigenetic pathway, HSP90 represses the regulatory influence of endogenous retroviruses (ERVs) on neighboring genes that are critical for mouse development. Our data based on natural variations in the mouse genome show that genes respond to HSP90 inhibition in a manner dependent on their genomic location with regard to strain-specific ERV-insertion sites. The evolutionary-capacitor function of HSP90 may thus have facilitated the exaptation of ERVs as key modifiers of gene expression and morphological diversification. Our findings add a new regulatory layer through which HSP90 uncouples phenotypic outcomes from individual genotypes.

The Composite Regulatory Basis of the Large X-Effect in Mouse Speciation

The disruption of meiotic sex chromosome inactivation (MSCI) has been proposed to be a major developmental mechanism underlying the rapid evolution of hybrid male sterility. We tested this idea by analyzing cell-specific gene expression across spermatogenesis in two lineages of house mice and their sterile and fertile reciprocal hybrids. We found pervasive disruption of sex chromosome gene expression in sterile hybrids at every stage of spermatogenesis. Failure of MSCI was developmentally preceded by increased silencing of autosomal genes, supporting the hypothesis that divergence at the hybrid incompatibility gene, Prdm9, results in increased rates of autosomal asynapsis which in turn triggers widespread silencing of unsynapsed chromatin. We also detected opposite patterns of postmeiotic overexpression or hyper-repression of the sex chromosomes in reciprocal hybrids, supporting the hypothesis that genomic conflict has driven functional divergence that leads to deleterious X-Y dosage imbalances in hybrids. Our developmental timeline also exposed more subtle patterns of mitotic misregulation on the X chromosome, a previously undocumented stage of spermatogenic disruption in this cross. These results indicate that multiple hybrid incompatibilities have converged on a common regulatory phenotype, the disrupted expression of the sex chromosomes during spermatogenesis. Collectively, these data reveal a composite regulatory basis to hybrid male sterility in mice that helps resolve the mechanistic underpinnings of the well-documented large X-effect in mice speciation. We propose that the inherent sensitivity of spermatogenesis to X-linked regulatory disruption has the potential to be a major driver of reproductive isolation in species with chromosomal sex determination.

Genomic imprinting, disrupted placental expression, and speciation

The importance of regulatory incompatibilities to the early stages of speciation remains unclear. Hybrid mammals often show extreme parent-of-origin growth effects that are thought to be a consequence of disrupted genetic imprinting (parent-specific epigenetic gene silencing) during early development. Here, we test the long-standing hypothesis that abnormal hybrid growth reflects disrupted gene expression due to loss of imprinting (LOI) in hybrid placentas, resulting in dosage imbalances between paternal growth factors and maternal growth repressors. We analyzed placental gene expression in reciprocal dwarf hamster hybrids that show extreme parent-of-origin growth effects relative to their parental species. In massively enlarged hybrid placentas, we observed both extensive transgressive expression of growth-related genes and biallelic expression of many genes that were paternally silenced in normal sized hybrids. However, the apparent widespread disruption of paternal silencing was coupled with reduced gene expression levels overall. These patterns are contrary to the predictions of the LOI model and indicate that hybrid misexpression of dosage-sensitive genes is caused by other regulatory mechanisms in this system. Collectively, our results support a central role for disrupted gene expression and imprinting in the evolution of mammalian hybrid inviability, but call into question the generality of the widely invoked LOI model.

Local regulation of gene expression by lncRNA promoters, transcription and splicing

Mammalian genomes are pervasively transcribed to produce thousands of long non-coding RNAs (lncRNAs). A few of these lncRNAs have been shown to recruit regulatory complexes through RNA-protein interactions to influence the expression of nearby genes, and it has been suggested that many other lncRNAs can also act as local regulators. Such local functions could explain the observation that lncRNA expression is often correlated with the expression of nearby genes. However, these correlations have been challenging to dissect and could alternatively result from processes that are not mediated by the lncRNA transcripts themselves. For example, some gene promoters have been proposed to have dual functions as enhancers, and the process of transcription itself may contribute to gene regulation by recruiting activating factors or remodelling nucleosomes. Here we use genetic manipulation in mouse cell lines to dissect 12 genomic loci that produce lncRNAs and find that 5 of these loci influence the expression of a neighbouring gene in cis. Notably, none of these effects requires the specific lncRNA transcripts themselves and instead involves general processes associated with their production, including enhancer-like activity of gene promoters, the process of transcription, and the splicing of the transcript. Furthermore, such effects are not limited to lncRNA loci: we find that four out of six protein-coding loci also influence the expression of a neighbour. These results demonstrate that cross-talk among neighbouring genes is a prevalent phenomenon that can involve multiple mechanisms and cis-regulatory signals, including a role for RNA splice sites. These mechanisms may explain the function and evolution of some genomic loci that produce lncRNAs and broadly contribute to the regulation of both coding and non-coding genes.

Contrasting Levels of Molecular Evolution on the Mouse X Chromosome

The mammalian X chromosome has unusual evolutionary dynamics compared to autosomes. Faster-X evolution of spermatogenic protein-coding genes is known to be most pronounced for genes expressed late in spermatogenesis, but it is unclear if these patterns extend to other forms of molecular divergence. We tested for faster-X evolution in mice spanning three different forms of molecular evolution-divergence in protein sequence, gene expression, and DNA methylation-across different developmental stages of spermatogenesis. We used FACS to isolate individual cell populations and then generated cell-specific transcriptome profiles across different stages of spermatogenesis in two subspecies of house mice (Mus musculus), thereby overcoming a fundamental limitation of previous studies on whole tissues. We found faster-X protein evolution at all stages of spermatogenesis and faster-late protein evolution for both X-linked and autosomal genes. In contrast, there was less expression divergence late in spermatogenesis (slower late) on the X chromosome and for autosomal genes expressed primarily in testis (testis-biased). We argue that slower-late expression divergence reflects strong regulatory constraints imposed during this critical stage of sperm development and that these constraints are particularly acute on the tightly regulated sex chromosomes. We also found slower-X DNA methylation divergence based on genome-wide bisulfite sequencing of sperm from two species of mice (M. musculus and M. spretus), although it is unclear whether slower-X DNA methylation reflects development constraints in sperm or other X-linked phenomena. Our study clarifies key differences in patterns of regulatory and protein evolution across spermatogenesis that are likely to have important consequences for mammalian sex chromosome evolution, male fertility, and speciation.

Genetic Architectures of Quantitative Variation in RNA Editing Pathways

RNA editing refers to post-transcriptional processes that alter the base sequence of RNA. Recently, hundreds of new RNA editing targets have been reported. However, the mechanisms that determine the specificity and degree of editing are not well understood. We examined quantitative variation of site-specific editing in a genetically diverse multiparent population, Diversity Outbred mice, and mapped polymorphic loci that alter editing ratios globally for C-to-U editing and at specific sites for A-to-I editing. An allelic series in the C-to-U editing enzyme Apobec1 influences the editing efficiency of Apob and 58 additional C-to-U editing targets. We identified 49 A-to-I editing sites with polymorphisms in the edited transcript that alter editing efficiency. In contrast to the shared genetic control of C-to-U editing, most of the variable A-to-I editing sites were determined by local nucleotide polymorphisms in proximity to the editing site in the RNA secondary structure. Our results indicate that RNA editing is a quantitative trait subject to genetic variation and that evolutionary constraints have given rise to distinct genetic architectures in the two canonical types of RNA editing.

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.

PRDM9 drives evolutionary erosion of hotspots in Mus musculus through haplotype-specific initiation of meiotic recombination

Meiotic recombination generates new genetic variation and assures the proper segregation of chromosomes in gametes. PRDM9, a zinc finger protein with histone methyltransferase activity, initiates meiotic recombination by binding DNA at recombination hotspots and directing the position of DNA double-strand breaks (DSB). The DSB repair mechanism suggests that hotspots should eventually self-destruct, yet genome-wide recombination levels remain constant, a conundrum known as the hotspot paradox. To test if PRDM9 drives this evolutionary erosion, we measured activity of the Prdm9^Cst allele in two Mus musculus subspecies, M.m. castaneus, in which Prdm9^Cst arose, and M.m. domesticus, into which Prdm9^Cst was introduced experimentally. Comparing these two strains, we find that haplotype differences at hotspots lead to qualitative and quantitative changes in PRDM9 binding and activity. Using Mus spretus as an outlier, we found most variants affecting PRDM9^Cst binding arose and were fixed in M.m. castaneus, suppressing hotspot activity. Furthermore, M.m. castaneus×M.m. domesticus F1 hybrids exhibit novel hotspots, with large haplotype biases in both PRDM9 binding and chromatin modification. These novel hotspots represent sites of historic evolutionary erosion that become activated in hybrids due to crosstalk between one parent's Prdm9 allele and the opposite parent's chromosome. Together these data support a model where haplotype-specific PRDM9 binding directs biased gene conversion at hotspots, ultimately leading to hotspot erosion.

Sexually reproducing creatures need to produce germ cells, notably sperm and egg, and do so using a specialized cell division, termed meiosis. A hallmark of meiosis is the process of recombination, in which pieces of maternal and paternal genetic material are exchanged, creating new combinations that are inherited by their progeny. Recombination increases diversity in subsequent generations, facilitating evolution. However, if recombination goes awry, it can lead to genetic disorders and spontaneous abortions in humans. In most mammals the sites of recombination are directed by the enzyme PRDM9 to specific regions on DNA, termed hotspots. For several decades it has been speculated that the process of recombination should lead to the eventual inactivation of hotspots, resulting in the loss of ability to reproduce. The discovery of PRDM9 provided a potential solution to this dilemma when the appearance of new PRDM9 alleles with altered DNA binding specificity would immediately create a new set of hotspots. We have now used the power of mouse genetics and large scale measurements of PRDM9 location and activity to show that this cycle of hotspot loss and recovery does indeed occur over the course of hundreds of thousands of years, and is directed by PRDM9.

The role of cis regulatory evolution in maize domestication

Gene expression differences between divergent lineages caused by modification of cis regulatory elements are thought to be important in evolution. We assayed genome-wide cis and trans regulatory differences between maize and its wild progenitor, teosinte, using deep RNA sequencing in F1 hybrid and parent inbred lines for three tissue types (ear, leaf and stem). Pervasive regulatory variation was observed with approximately 70% of ∼17,000 genes showing evidence of regulatory divergence between maize and teosinte. However, many fewer genes (1,079 genes) show consistent cis differences with all sampled maize and teosinte lines. For ∼70% of these 1,079 genes, the cis differences are specific to a single tissue. The number of genes with cis regulatory differences is greatest for ear tissue, which underwent a drastic transformation in form during domestication. As expected from the domestication bottleneck, maize possesses less cis regulatory variation than teosinte with this deficit greatest for genes showing maize-teosinte cis regulatory divergence, suggesting selection on cis regulatory differences during domestication. Consistent with selection on cis regulatory elements, genes with cis effects correlated strongly with genes under positive selection during maize domestication and improvement, while genes with trans regulatory effects did not. We observed a directional bias such that genes with cis differences showed higher expression of the maize allele more often than the teosinte allele, suggesting domestication favored up-regulation of gene expression. Finally, this work documents the cis and trans regulatory changes between maize and teosinte in over 17,000 genes for three tissues.