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

Bibliometric analysis on cannibalism/infanticide and maternal aggression towards pups in laboratory rodents

Animal welfare has evolved during the past decades to improve not only the quality of life of laboratory rodents but also the quality and reproducibility of scientific investigations. Bibliometric analysis has become an important tool to complete the current knowledge with academic databases. Our objective was to investigate whether scientific research on cannibalism/infanticide is connected with maternal aggression towards the offspring in laboratory rodents. To carry out our research, we performed a specific search for published articles on each concept. Results were analyzed in the open-source environment RStudio with the package Bibliometrix. We obtained 253 and 134 articles for the first search (cannibalism/infanticide) and the second search (maternal aggression towards the pups) respectively. We observed that the interest in infanticide/cannibalism started in the 1950s, while researchers started showing interest in maternal aggression towards the pups 30 years later. Our analyses indicated that maternal aggression had better citations in scientific literature. In addition, although our results showed some common features (e.g. oxytocin or medial preoptic area in the brain), we observed a gap between cannibalism/infanticide and maternal aggression towards the pups with only 14 published articles in common for both the searches. Therefore, we recommend researchers to combine both concepts in further investigations in the context of cannibalism for better dissemination and higher impact in laboratory rodents' welfare research.

Systematic evaluation of single-cell RNA-seq analyses performance based on long-read sequencing platforms

INTRODUCTION

The rapid development of next-generation sequencing (NGS)-based single-cell RNA sequencing (scRNA-seq) allows for detecting and quantifying gene expression in a high-throughput manner, providing a powerful tool for comprehensively understanding cellular function in various biological processes. However, the NGS-based scRNA-seq only quantifies gene expression and cannot reveal the exact transcript structures (isoforms) of each gene due to the limited read length. On the other hand, the long read length of third-generation sequencing (TGS) technologies, including Oxford Nanopore Technologies (ONT) and Pacific Biosciences (PacBio), enable direct reading of intact cDNA molecules.

OBJECTIVES

Both ONT and PacBio have been used in conjunction with scRNA-seq, but their performance in single-cell analyses has not been systematically evaluated.

METHODS

To address this, we generated ONT and PacBio data from the same single-cell cDNA libraries containing different amount of cells.

RESULTS

Using NGS as a control, we assessed the performance of each platform in cell type identification. Additionally, the reliability in identifying novel isoforms and allele-specific gene/isoform expression by both platforms was verified, providing a systematic evaluation to design the sequencing strategies in single-cell transcriptome studies.

CONCLUSION

Beyond gene expression analysis, which the NGS-based scRNA-seq only affords, TGS-based scRNA-seq achieved gene splicing analyses, identifying novel isoforms. Attribute to higher sequencing quality of PacBio, it outperforms ONT in accuracy of novel transcripts identification and allele-specific gene/isoform expression.

The pivotal role of the X-chromosome in the genetic architecture of the human brain

Genes on the X-chromosome are extensively expressed in the human brain. However, little is known for the X-chromosome's impact on the brain anatomy, microstructure, and functional network. We examined 1,045 complex brain imaging traits from 38,529 participants in the UK Biobank. We unveiled potential autosome-X-chromosome interactions, while proposing an atlas outlining dosage compensation (DC) for brain imaging traits. Through extensive association studies, we identified 72 genome-wide significant trait-locus pairs (including 29 new associations) that share genetic architectures with brain-related disorders, notably schizophrenia. Furthermore, we discovered unique sex-specific associations and assessed variations in genetic effects between sexes. Our research offers critical insights into the X-chromosome's role in the human brain, underscoring its contribution to the differences observed in brain structure and functionality between sexes.

Extensive folding variability between homologous chromosomes in mammalian cells

Genetic variation and 3D chromatin structure have major roles in gene regulation. Due to challenges in mapping chromatin conformation with haplotype-specific resolution, the effects of genetic sequence variation on 3D genome structure and gene expression imbalance remain understudied. Here, we applied Genome Architecture Mapping (GAM) to a hybrid mouse embryonic stem cell (mESC) line with high density of single nucleotide polymorphisms (SNPs). GAM resolved haplotype-specific 3D genome structures with high sensitivity, revealing extensive allelic differences in chromatin compartments, topologically associating domains (TADs), long-range enhancer-promoter contacts, and CTCF loops. Architectural differences often coincide with allele-specific differences in gene expression, mediated by Polycomb repression. We show that histone genes are expressed with allelic imbalance in mESCs, are involved in haplotype-specific chromatin contact marked by H3K27me3, and are targets of Polycomb repression through conditional knockouts of Ezh2 or Ring1b. Our work reveals highly distinct 3D folding structures between homologous chromosomes, and highlights their intricate connections with allelic gene expression.

The genetic basis of novel trait gain in walking fish

A major goal in biology is to understand how organisms evolve novel traits. Multiple studies have identified genes contributing to regressive evolution, the loss of structures that existed in a recent ancestor. However, fewer examples exist for genes underlying constructive evolution, the gain of novel structures and capabilities in lineages that previously lacked them. Sea robins are fish that have evolved enlarged pectoral fins, six mobile locomotory fin rays (legs) and six novel macroscopic lobes in the central nervous system (CNS) that innervate the corresponding legs. Here, we establish successful husbandry and use a combination of transcriptomics, CRISPR-Cas9 editing, and behavioral assays to identify key transcription factors that are required for leg formation and function in sea robins. We also generate hybrids between two sea robin species with distinct leg morphologies and use allele-specific expression analysis and gene editing to explore the genetic basis of species-specific trait diversity, including a novel sensory gain of function. Collectively, our study establishes sea robins as a new model for studying the genetic basis of novel organ formation, and demonstrates a crucial role for the conserved limb gene tbx3a in the evolution of chemosensory legs in walking fish.

Environmentally robust cis-regulatory changes underlie rapid climatic adaptation

Changes in gene expression are thought to play a major role in adaptive evolution. While it is known that gene expression is highly sensitive to the environment, very few studies have determined the influence of genetic and environmental effects on adaptive gene expression differences in natural populations. Here, we utilize allele-specific expression to characterize cis and trans gene regulatory divergence in temperate and tropical house mice in two metabolic tissues under two thermal conditions. First, we show that gene expression divergence is pervasive between populations and across thermal conditions, with roughly 5 to 10% of genes exhibiting genotype-by-environment interactions. Second, we found that most expression divergence was due to cis-regulatory changes that were stable across temperatures. In contrast, patterns of expression plasticity were largely attributable to trans-effects, which showed greater sensitivity to temperature. Nonetheless, we found a small subset of temperature-dependent cis-regulatory changes, thereby identifying loci underlying expression plasticity. Finally, we performed scans for selection in wild house mice to identify genomic signatures of rapid adaptation. Genomic outliers were enriched in genes with evidence for cis-regulatory divergence. Notably, these genes were associated with phenotypes that affected body weight and metabolism, suggesting that cis-regulatory changes are a possible mechanism for adaptive body size evolution between populations. Our results show that gene expression plasticity, largely controlled in trans, may facilitate the colonization of new environments, but that evolved changes in gene expression are largely controlled in cis, illustrating the genetic and nongenetic mechanisms underlying the establishment of populations in new environments.

SEESAW: detecting isoform-level allelic imbalance accounting for inferential uncertainty

Detecting allelic imbalance at the isoform level requires accounting for inferential uncertainty, caused by multi-mapping of RNA-seq reads. Our proposed method, SEESAW, uses Salmon and Swish to offer analysis at various levels of resolution, including gene, isoform, and aggregating isoforms to groups by transcription start site. The aggregation strategies strengthen the signal for transcripts with high uncertainty. The SEESAW suite of methods is shown to have higher power than other allelic imbalance methods when there is isoform-level allelic imbalance. We also introduce a new test for detecting imbalance that varies across a covariate, such as time.

Inference of putative cell-type-specific imprinted regulatory elements and genes during human neuronal differentiation

Genomic imprinting results in gene expression bias caused by parental chromosome of origin and occurs in genes with important roles during human brain development. However, the cell-type and temporal specificity of imprinting during human neurogenesis is generally unknown. By detecting within-donor allelic biases in chromatin accessibility and gene expression that are unrelated to cross-donor genotype, we inferred imprinting in both primary human neural progenitor cells and their differentiated neuronal progeny from up to 85 donors. We identified 43/20 putatively imprinted regulatory elements (IREs) in neurons/progenitors, and 133/79 putatively imprinted genes in neurons/progenitors. Although 10 IREs and 42 genes were shared between neurons and progenitors, most putative imprinting was only detected within specific cell types. In addition to well-known imprinted genes and their promoters, we inferred novel putative IREs and imprinted genes. Consistent with both DNA methylation-based and H3K27me3-based regulation of imprinted expression, some putative IREs also overlapped with differentially methylated or histone-marked regions. Finally, we identified a progenitor-specific putatively imprinted gene overlapping with copy number variation that is associated with uniparental disomy-like phenotypes. Our results can therefore be useful in interpreting the function of variants identified in future parent-of-origin association studies.

Allele-specific regulatory effects on the pig transcriptome

BACKGROUND

Allele-specific expression (ASE) refers to the preferential expression of one allele over the other and contributes to adaptive phenotypic plasticity. Here, we used a reciprocal cross-model between phenotypically divergent European Berkshire and Asian Tibetan pigs to characterize 2 ASE classes: imprinting (i.e., the unequal expression between parental alleles) and sequence dependent (i.e., unequal expression between breed-specific alleles). We examined 3 transcript types, including protein-coding genes (PCGs), long noncoding RNAs, and transcripts of unknown coding potential, across 7 representative somatic tissues from hybrid pigs generated by reciprocal crosses.

RESULTS

We identified a total of 92 putative imprinted transcripts, 69 (75.00%) of which are described here for the first time. By combining the transcriptome from purebred Berkshire and Tibetan pigs, we found ∼6.59% of PCGs are differentially expressed between breeds that are regulated by trans-elements (e.g., transcriptional factors), while only ∼1.35% are attributable to cis (e.g., promoters). The higher prevalence of trans-PCGs indicates the dominated effects of trans-regulation in driving expression differences and shaping adaptive phenotypic plasticity between breeds, which were supported by functional enrichment analysis. We also found strong evidence that expression changes mediated by cis-effects were associated with accumulated variants in promoters.

CONCLUSIONS

Our study provides a comprehensive map of expression regulation that constitutes a valuable resource for the agricultural improvement of pig breeds.

Variable allelic expression of imprinted genes at the Peg13, Trappc9, Ago2 cluster in single neural cells

Genomic imprinting is an epigenetic process through which genes are expressed in a parent-of-origin specific manner resulting in mono-allelic or strongly biased expression of one allele. For some genes, imprinted expression may be tissue-specific and reliant on CTCF-influenced enhancer-promoter interactions. The Peg13 imprinting cluster is associated with neurodevelopmental disorders and comprises canonical imprinted genes, which are conserved between mouse and human, as well as brain-specific imprinted genes in mouse. The latter consist of Trappc9, Chrac1 and Ago2, which have a maternal allelic expression bias of ∼75% in brain. Findings of such allelic expression biases on the tissue level raise the question of how they are reflected in individual cells and whether there is variability and mosaicism in allelic expression between individual cells of the tissue. Here we show that Trappc9 and Ago2 are not imprinted in hippocampus-derived neural stem cells (neurospheres), while Peg13 retains its strong bias of paternal allele expression. Upon analysis of single neural stem cells and in vitro differentiated neurons, we find not uniform, but variable states of allelic expression, especially for Trappc9 and Ago2. These ranged from mono-allelic paternal to equal bi-allelic to mono-allelic maternal, including biased bi-allelic transcriptional states. Even Peg13 expression deviated from its expected paternal allele bias in a small number of cells. Although the cell populations consisted of a mosaic of cells with different allelic expression states, as a whole they reflected bulk tissue data. Furthermore, in an attempt to identify potential brain-specific regulatory elements across the Trappc9 locus, we demonstrate tissue-specific and general silencer activities, which might contribute to the regulation of its imprinted expression bias.

Sex- and age-dependent genetics of longevity in a heterogeneous mouse population

DNA variants that modulate life span provide insight into determinants of health, disease, and aging. Through analyses in the UM-HET3 mice of the Interventions Testing Program (ITP), we detected a sex-independent quantitative trait locus (QTL) on chromosome 12 and identified sex-specific QTLs, some of which we detected only in older mice. Similar relations between life history and longevity were uncovered in mice and humans, underscoring the importance of early access to nutrients and early growth. We identified common age- and sex-specific genetic effects on gene expression that we integrated with model organism and human data to create a hypothesis-building interactive resource of prioritized longevity and body weight genes. Finally, we validated Hipk1, Ddost, Hspg2, Fgd6, and Pdk1 as conserved longevity genes using Caenorhabditis elegans life-span experiments.

Evolutionary changes in the chromatin landscape contribute to reorganization of a developmental gene network during rapid life history evolution in sea urchins

Chromatin configuration is highly dynamic during embryonic development in animals, exerting an important point of control in transcriptional regulation. Yet there exists remarkably little information about the role of evolutionary changes in chromatin configuration to the evolution of gene expression and organismal traits. Genome-wide assays of chromatin configuration, coupled with whole-genome alignments, can help address this gap in knowledge in several ways. In this study we present a comparative analysis of regulatory element sequences and accessibility throughout embryogenesis in three sea urchin species with divergent life histories: a lecithotroph Heliocidaris erythrogramma, a closely related planktotroph H. tuberculata, and a distantly related planktotroph Lytechinus variegatus. We identified distinct epigenetic and mutational signatures of evolutionary modifications to the function of putative cis-regulatory elements in H. erythrogramma that have accumulated nonuniformly throughout the genome, suggesting selection, rather than drift, underlies many modifications associated with the derived life history. Specifically, regulatory elements composing the sea urchin developmental gene regulatory network are enriched for signatures of positive selection and accessibility changes which may function to alter binding affinity and access of developmental transcription factors to these sites. Furthermore, regulatory element changes often correlate with divergent expression patterns of genes involved in cell type specification, morphogenesis, and development of other derived traits, suggesting these evolutionary modifications have been consequential for phenotypic evolution in H. erythrogramma. Collectively, our results demonstrate that selective pressures imposed by changes in developmental life history rapidly reshape the cis-regulatory landscape of core developmental genes to generate novel traits and embryonic programs.

Favoring Expression of Yak Alleles in Interspecies F1 Hybrids of Cattle and Yak Under High-Altitude Environments

Both cis- and trans-regulation could cause differential expression between the parental alleles in diploid species that might have broad biological implications. Due to the relatively distant genetic divergence between cattle and yak, as well as their differential adaptation to high-altitude environments, we investigated genome-wide allelic differential expression (ADE) in their F1 hybrids using Nanopore long-read RNA-seq technology. From adult F1 hybrids raised in high-altitude, ten lung and liver tissues were individually sequenced for producing 31.6 M full-length transcript sequences. Mapping against autosomal homologous regions between cattle and yak, we detected 17,744 and 14,542 protein-encoding genes expressed in lung and liver tissues, respectively. According to the parental assignments of transcript sequences, a total of 3,381 genes were detected to show ADE in at least one sample. There were 186 genes showing ubiquitous ADE in all the studied animals, and among them 135 and 37 genes had consistent higher expression of yak and cattle alleles, respectively. Functional analyses revealed that the genes with favoring expression of yak alleles have been involved in the biological progresses related with hypoxia adaptation and immune response. In contrast, the genes with favoring expression of cattle alleles have been enriched into different biological progresses, such as secretion of endocrine hormones and lipid metabolism. Our results would support unequal contribution of parental genes to environmental adaptation in the F1 hybrids of cattle and yak.

Gene expression changes following chronic antipsychotic exposure in single cells from mouse striatum

Schizophrenia is an idiopathic psychiatric disorder with a high degree of polygenicity. Evidence from genetics, single-cell transcriptomics, and pharmacological studies suggest an important, but untested, overlap between genes involved in the etiology of schizophrenia and the cellular mechanisms of action of antipsychotics. To directly compare genes with antipsychotic-induced differential expression to genes involved in schizophrenia, we applied single-cell RNA-sequencing to striatal samples from male C57BL/6 J mice chronically exposed to a typical antipsychotic (haloperidol), an atypical antipsychotic (olanzapine), or placebo. We identified differentially expressed genes in three cell populations identified from the single-cell RNA-sequencing (medium spiny neurons [MSNs], microglia, and astrocytes) and applied multiple analysis pipelines to contextualize these findings, including comparison to GWAS results for schizophrenia. In MSNs in particular, differential expression analysis showed that there was a larger share of differentially expressed genes (DEGs) from mice treated with olanzapine compared with haloperidol. DEGs were enriched in loci implicated by genetic studies of schizophrenia, and we highlighted nine genes with convergent evidence. Pathway analyses of gene expression in MSNs highlighted neuron/synapse development, alternative splicing, and mitochondrial function as particularly engaged by antipsychotics. In microglia, we identified pathways involved in microglial activation and inflammation as part of the antipsychotic response. In conclusion, single-cell RNA sequencing may provide important insights into antipsychotic mechanisms of action and links to findings from psychiatric genomic studies.

Comparative 3D genome architecture in vertebrates

Background

The three-dimensional (3D) architecture of the genome has a highly ordered and hierarchical nature, which influences the regulation of essential nuclear processes at the basis of gene expression, such as gene transcription. While the hierarchical organization of heterochromatin and euchromatin can underlie differences in gene expression that determine evolutionary differences among species, the way 3D genome architecture is affected by evolutionary forces within major lineages remains unclear. Here, we report a comprehensive comparison of 3D genomes, using high resolution Hi-C data in fibroblast cells of fish, chickens, and 10 mammalian species.

Results

This analysis shows a correlation between genome size and chromosome length that affects chromosome territory (CT) organization in the upper hierarchy of genome architecture, whereas lower hierarchical features, including local transcriptional availability of DNA, are selected through the evolution of vertebrates. Furthermore, conservation of topologically associating domains (TADs) appears strongly associated with the modularity of expression profiles across species. Additionally, LINE and SINE transposable elements likely contribute to heterochromatin and euchromatin organization, respectively, during the evolution of genome architecture.

Conclusions

Our analysis uncovers organizational features that appear to determine the conservation and transcriptional regulation of functional genes across species. These findings can guide ongoing investigations of genome evolution by extending our understanding of the mechanisms shaping genome architecture.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01301-7.

Genetic variants of the hypoxia-inducible factor 3 alpha subunit (Hif3a) gene in the Fat and Lean mouse selection lines

BACKGROUND

Adipose tissue hypoxia and members of the hypoxia-inducible factor alpha (HIFA) are involved in development of obesity. However, the mechanism and functions of HIF3A, one of three HIFA paralogs, in fat deposition have not been sufficiently studied.

METHODS AND RESULTS

In the present study, we investigated whether Hif3a sequence variants are associated with divergent fat deposition in mouse selection lines for fatness and leanness. Sequencing and RFLP were used to analyse sequence variants within Hif3a. To identify candidate regulatory variants, we performed literature screening and used databases and bioinformatics tools like Ensembl, MethPrimer, TargetScanMouse, miRDB, PolyAsite, RISE, LncRRIsearch, RNAfold, PredictProtein, CAIcal, and switches.ELM Resource. There are 90 sequence variants in Hif3a between the two mouse lines. While most Fat line variants locate within intronic regions, Lean line variants are mainly in 3' UTR. We constructed a map of Hif3a potential regulatory regions and identified 39 regulatory variants by integrating data on constrained and regulatory elements, CpGs, and miRNAs and lncRNAs binding sites. Moreover, 3' UTR and two exonic variants may influence mRNA stability, translation rate and protein functionality. We propose as priority candidates for further functional studies a missense (rs37398126) and synonymous (rs37739792) variants, and intronic (rs47471302) variant that overlap conserved element in promoter region and predicted lncRNAs binding site.

CONCLUSION

The results indicate a potential involvement of Hif3a in fat deposition. Additionally, approach used in the present study may serve as a general guideline for constructing an integrative gene map for prioritizing candidate gene variants with phenotypic effects.

Noncanonical genomic imprinting in the monoamine system determines naturalistic foraging and brain-adrenal axis functions

Noncanonical genomic imprinting can cause biased expression of one parental allele in a tissue; however, the functional relevance of such biases is unclear. To investigate ethological roles for noncanonical imprinting in dopa decarboxylase (Ddc) and tyrosine hydroxylase (Th), we use machine learning to decompose naturalistic foraging in maternal and paternal allele mutant heterozygous mice. We uncover distinct roles for the maternal versus paternal alleles on foraging, where maternal alleles affect sons while daughters are under paternal allelic control. Each parental allele controls specific action sequences reflecting decisions in naive or familiar contexts. The maternal Ddc allele is preferentially expressed in subsets of hypothalamic GABAergic neurons, while the paternal allele predominates in subsets of adrenal cells. Each Ddc allele affects distinct molecular and endocrine components of the brain-adrenal axis. Thus, monoaminergic noncanonical imprinting has ethological roles in foraging and endocrine functions and operates by affecting discrete subsets of cells.

Transcriptional kinetics and molecular functions of long noncoding RNAs

An increasing number of long noncoding RNAs (lncRNAs) have experimentally confirmed functions, yet little is known about their transcriptional dynamics and it is challenging to determine their regulatory effects. Here, we used allele-sensitive single-cell RNA sequencing to demonstrate that, compared to messenger RNAs, lncRNAs have twice as long duration between two transcriptional bursts. Additionally, we observed increased cell-to-cell variability in lncRNA expression due to lower frequency bursting producing larger numbers of RNA molecules. Exploiting heterogeneity in asynchronously growing cells, we identified and experimentally validated lncRNAs with cell state-specific functions involved in cell cycle progression and apoptosis. Finally, we identified cis-functioning lncRNAs and showed that knockdown of these lncRNAs modulated the nearby protein-coding gene’s transcriptional burst frequency or size. In summary, we identified distinct transcriptional regulation of lncRNAs and demonstrated a role for lncRNAs in the regulation of mRNA transcriptional bursting.

Allele-sensitive single-cell RNA sequencing analysis of long noncoding RNA (lncRNA) transcriptional kinetics shows that their lower expression compared to mRNA is due to lower burst frequencies and highlights cell-state-specific functions for several lncRNAs.

Extreme phenotypic divergence and the evolution of development

As analyses of developmental mechanisms extend to ever more species, it becomes important to understand not just what is conserved or altered during evolution, but why. Closely related species that exhibit extreme phenotypic divergence can be uniquely informative in this regard. A case in point is the sea urchin genus Heliocidaris, which contains species that recently evolved a life history involving nonfeeding larvae following nearly half a billion years of prior evolution with feeding larvae. The resulting shift in selective regimes produced rapid and surprisingly extensive changes in developmental mechanisms that are otherwise highly conserved among echinoderm species. The magnitude and extent of these changes challenges the notion that conservation of early development in echinoderms is largely due to internal constraints that prohibit modification and instead suggests that natural selection actively maintains stability of inherently malleable trait developmental mechanisms over immense time periods. Knowing how and why natural selection changed during the evolution of nonfeeding larvae can also reveal why developmental mechanisms do and do not change in particular ways.

cis-Regulatory changes in locomotor genes are associated with the evolution of burrowing behavior

How evolution modifies complex, innate behaviors is largely unknown. Divergence in many morphological traits, and some behaviors, is linked to cis-regulatory changes in gene expression. Given this, we compare brain gene expression of two interfertile sister species of Peromyscus mice that show large and heritable differences in burrowing behavior. Species-level differential expression and allele-specific expression in F₁ hybrids indicate a preponderance of cis-regulatory divergence, including many genes whose cis-regulation is affected by burrowing behavior. Genes related to locomotor coordination show the strongest signals of lineage-specific selection on burrowing-induced cis-regulatory changes. Furthermore, genetic markers closest to these candidate genes associate with variation in burrow shape in a genetic cross, suggesting an enrichment for loci affecting burrowing behavior near these candidate locomotor genes. Our results provide insight into how cis-regulated gene expression can depend on behavioral context and how this dynamic regulatory divergence between species may contribute to behavioral evolution.

Power calculator for detecting allelic imbalance using hierarchical Bayesian model

OBJECTIVE

Allelic imbalance (AI) is the differential expression of the two alleles in a diploid. AI can vary between tissues, treatments, and environments. Methods for testing AI exist, but methods are needed to estimate type I error and power for detecting AI and difference of AI between conditions. As the costs of the technology plummet, what is more important: reads or replicates?

RESULTS

We find that a minimum of 2400, 480, and 240 allele specific reads divided equally among 12, 5, and 3 replicates is needed to detect a 10, 20, and 30%, respectively, deviation from allelic balance in a condition with power > 80%. A minimum of 960 and 240 allele specific reads divided equally among 8 replicates is needed to detect a 20 or 30% difference in AI between conditions with comparable power. Higher numbers of replicates increase power more than adding coverage without affecting type I error. We provide a Python package that enables simulation of AI scenarios and enables individuals to estimate type I error and power in detecting AI and differences in AI between conditions.

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.

Single-nuclei chromatin profiling of ventral midbrain reveals cell identity transcription factors and cell-type-specific gene regulatory variation

BACKGROUND

Cell types in ventral midbrain are involved in diseases with variable genetic susceptibility, such as Parkinson's disease and schizophrenia. Many genetic variants affect regulatory regions and alter gene expression in a cell-type-specific manner depending on the chromatin structure and accessibility.

RESULTS

We report 20,658 single-nuclei chromatin accessibility profiles of ventral midbrain from two genetically and phenotypically distinct mouse strains. We distinguish ten cell types based on chromatin profiles and analysis of accessible regions controlling cell identity genes highlights cell-type-specific key transcription factors. Regulatory variation segregating the mouse strains manifests more on transcriptome than chromatin level. However, cell-type-level data reveals changes not captured at tissue level. To discover the scope and cell-type specificity of cis-acting variation in midbrain gene expression, we identify putative regulatory variants and show them to be enriched at differentially expressed loci. Finally, we find TCF7L2 to mediate trans-acting variation selectively in midbrain neurons.

CONCLUSIONS

Our data set provides an extensive resource to study gene regulation in mesencephalon and provides insights into control of cell identity in the midbrain and identifies cell-type-specific regulatory variation possibly underlying phenotypic and behavioural differences between mouse strains.

The integrative biology of genetic dominance

Dominance is a basic property of inheritance systems describing the link between a diploid genotype at a single locus and the resulting phenotype. Models for the evolution of dominance have long been framed as an opposition between the irreconcilable views of Fisher in 1928 supporting the role of largely elusive dominance modifiers and Wright in 1929, who viewed dominance as an emerging property of the structure of enzymatic pathways. Recent theoretical and empirical advances however suggest that these opposing views can be reconciled, notably using models investigating the regulation of gene expression and developmental processes. In this more comprehensive framework, phenotypic dominance emerges from departures from linearity between any levels of integration in the genotype‐to‐phenotype map. Here, we review how these different models illuminate the emergence and evolution of dominance. We then detail recent empirical studies shedding new light on the diversity of molecular and physiological mechanisms underlying dominance and its evolution. By reconciling population genetics and functional biology, we hope our review will facilitate cross‐talk among research fields in the integrative study of dominance evolution.

Perspectives on Allele-Specific Expression

Diploidy has profound implications for population genetics and susceptibility to genetic diseases. Although two copies are present for most genes in the human genome, they are not necessarily both active or active at the same level in a given individual. Genomic imprinting, resulting in exclusive or biased expression in favor of the allele of paternal or maternal origin, is now believed to affect hundreds of human genes. A far greater number of genes display unequal expression of gene copies due to cis-acting genetic variants that perturb gene expression. The availability of data generated by RNA sequencing applied to large numbers of individuals and tissue types has generated unprecedented opportunities to assess the contribution of genetic variation to allelic imbalance in gene expression. Here we review the insights gained through the analysis of these data about the extent of the genetic contribution to allelic expression imbalance, the tools and statistical models for gene expression imbalance, and what the results obtained reveal about the contribution of genetic variants that alter gene expression to complex human diseases and phenotypes.

Replicate sequencing libraries are important for quantification of allelic imbalance

A sensitive approach to quantitative analysis of transcriptional regulation in diploid organisms is analysis of allelic imbalance (AI) in RNA sequencing (RNA-seq) data. A near-universal practice in such studies is to prepare and sequence only one library per RNA sample. We present theoretical and experimental evidence that data from a single RNA-seq library is insufficient for reliable quantification of the contribution of technical noise to the observed AI signal; consequently, reliance on one-replicate experimental design can lead to unaccounted-for variation in error rates in allele-specific analysis. We develop a computational approach, Qllelic, that accurately accounts for technical noise by making use of replicate RNA-seq libraries. Testing on new and existing datasets shows that application of Qllelic greatly decreases false positive rate in allele-specific analysis while conserving appropriate signal, and thus greatly improves reproducibility of AI estimates. We explore sources of technical overdispersion in observed AI signal and conclude by discussing design of RNA-seq studies addressing two biologically important questions: quantification of transcriptome-wide AI in one sample, and differential analysis of allele-specific expression between samples.

A unified atlas of CD8 T cell dysfunctional states in cancer and infection

CD8 T cells play an essential role in defense against viral and bacterial infections and in tumor immunity. Deciphering T cell loss of functionality is complicated by the conspicuous heterogeneity of CD8 T cell states described across experimental and clinical settings. By carrying out a unified analysis of over 300 assay for transposase-accessible chromatin sequencing (ATAC-seq) and RNA sequencing (RNA-seq) experiments from 12 studies of CD8 T cells in cancer and infection, we defined a shared differentiation trajectory toward dysfunction and its underlying transcriptional drivers and revealed a universal early bifurcation of functional and dysfunctional T cell states across models. Experimental dissection of acute and chronic viral infection using single-cell ATAC (scATAC)-seq and allele-specific single-cell RNA (scRNA)-seq identified state-specific drivers and captured the emergence of similar TCF1+ progenitor-like populations at an early branch point, at which functional and dysfunctional T cells diverge. Our atlas of CD8 T cell states will facilitate mechanistic studies of T cell immunity and translational efforts.

Males That Silence Their Father's Genes: Genomic Imprinting of a Complete Haploid Genome

Genetic conflict is considered a key driver in the evolution of reproductive systems with non-Mendelian inheritance, where parents do not contribute equally to the genetic makeup of their offspring. One of the most extraordinary examples of non-Mendelian inheritance is paternal genome elimination (PGE), a form of haplodiploidy which has evolved repeatedly across arthropods. Under PGE, males are diploid but only transmit maternally inherited chromosomes, while the paternally inherited homologues are excluded from sperm. This asymmetric inheritance is thought to have evolved through an evolutionary arms race between the paternal and maternal genomes over transmission to future generations. In several PGE clades, such as the mealybugs (Hemiptera: Pseudococcidae), paternal chromosomes are not only eliminated from sperm, but also heterochromatinized early in development and thought to remain inactive, which could result from genetic conflict between parental genomes. Here, we present a parent-of-origin allele-specific transcriptome analysis in male mealybugs showing that expression is globally biased toward the maternal genome. However, up to 70% of somatically expressed genes are to some degree paternally expressed, while paternal genome expression is much more restricted in the male reproductive tract, with only 20% of genes showing paternal contribution. We also show that parent-of-origin-specific gene expression patterns are remarkably similar across genotypes, and that genes with completely biparental expression show elevated rates of molecular evolution. Our results provide the clearest example yet of genome-wide genomic imprinting in insects and enhance our understanding of PGE, which will aid future empirical tests of evolutionary theory regarding the origin of this unusual reproductive strategy.

Tissue context determines the penetrance of regulatory DNA variation

Functional assessment of disease-associated sequence variation at non-coding regulatory elements is complicated by their high degree of context sensitivity to both the local chromatin and nuclear environments. Allelic profiling of DNA accessibility across individuals has shown that only a select minority of sequence variation affects transcription factor (TF) occupancy, yet low sequence diversity in human populations means that no experimental assessment is available for the majority of disease-associated variants. Here we describe high-resolution in vivo maps of allelic DNA accessibility in liver, kidney, lung and B cells from 5 increasingly diverged strains of F1 hybrid mice. The high density of heterozygous sites in these hybrids enables precise quantification of effect size and cell-type specificity for hundreds of thousands of variants throughout the mouse genome. We show that chromatin-altering variants delineate characteristic sensitivity profiles for hundreds of TF motifs. We develop a compendium of TF-specific sensitivity profiles accounting for genomic context effects. Finally, we link maps of allelic accessibility to allelic transcript levels in the same samples. This work provides a foundation for quantitative prediction of cell-type specific effects of non-coding variation on TF activity, which will facilitate both fine-mapping and systems-level analyses of common disease-associated variation in human genomes.

Testcrosses are an efficient strategy for identifying cis-regulatory variation: Bayesian analysis of allele-specific expression (BayesASE)

Allelic imbalance (AI) occurs when alleles in a diploid individual are differentially expressed and indicates cis acting regulatory variation. What is the distribution of allelic effects in a natural population? Are all alleles the same? Are all alleles distinct? The approach described applies to any technology generating allele-specific sequence counts, for example for chromatin accessibility and can be applied generally including to comparisons between tissues or environments for the same genotype. Tests of allelic effect are generally performed by crossing individuals and comparing expression between alleles directly in the F1. However, a crossing scheme that compares alleles pairwise is a prohibitive cost for more than a handful of alleles as the number of crosses is at least (n2-n)/2 where n is the number of alleles. We show here that a testcross design followed by a hypothesis test of AI between testcrosses can be used to infer differences between nontester alleles, allowing n alleles to be compared with n crosses. Using a mouse data set where both testcrosses and direct comparisons have been performed, we show that the predicted differences between nontester alleles are validated at levels of over 90% when a parent-of-origin effect is present and of 60%-80% overall. Power considerations for a testcross, are similar to those in a reciprocal cross. In all applications, the testing for AI involves several complex bioinformatics steps. BayesASE is a complete bioinformatics pipeline that incorporates state-of-the-art error reduction techniques and a flexible Bayesian approach to estimating AI and formally comparing levels of AI between conditions. The modular structure of BayesASE has been packaged in Galaxy, made available in Nextflow and as a collection of scripts for the SLURM workload manager on github (https://github.com/McIntyre-Lab/BayesASE).

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.

aScan: A Novel Method for the Study of Allele Specific Expression in Single Individuals

In diploid organisms, two copies of each allele are normally inherited from parents. Paternal and maternal alleles can be regulated and expressed unequally, which is referred to as allele-specific expression (ASE). In this work, we present aScan, a novel method for the identification of ASE from the analysis of matched individual genomic and RNA sequencing data. By performing extensive analyses of both real and simulated data, we demonstrate that aScan can correctly identify ASE with high accuracy and sensitivity in different experimental settings. Additionally, by applying our method to a small cohort of individuals that are not included in publicly available databases of human genetic variation, we outline the value of possible applications of ASE analysis in single individuals for deriving a more accurate annotation of "private" low-frequency genetic variants associated with regulatory effects on transcription. All in all, we believe that aScan will represent a beneficial addition to the set of bioinformatics tools for the analysis of ASE. Finally, while our method was initially conceived for the analysis of RNA-seq data, it can in principle be applied to any quantitative NGS assay for which matched genotypic and expression data are available. AVAILABILITY: aScan is currently available in the form of an open source standalone software package at: https://github.com/Federico77z/aScan/. aScan version 1.0.3, available at https://github.com/Federico77z/aScan/releases/tag/1.0.3, has been used for all the analyses included in this manuscript. A Docker image of the tool has also been made available at https://github.com/pmandreoli/aScanDocker.

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.

Antipsychotic Behavioral Phenotypes in the Mouse Collaborative Cross Recombinant Inbred Inter-Crosses (RIX)

Schizophrenia is an idiopathic disorder that affects approximately 1% of the human population, and presents with persistent delusions, hallucinations, and disorganized behaviors. Antipsychotics are the standard pharmacological treatment for schizophrenia, but are frequently discontinued by patients due to inefficacy and/or side effects. Chronic treatment with the typical antipsychotic haloperidol causes tardive dyskinesia (TD), which manifests as involuntary and often irreversible orofacial movements in around 30% of patients. Mice treated with haloperidol develop many of the features of TD, including jaw tremors, tongue protrusions, and vacuous chewing movements (VCMs). In this study, we used genetically diverse Collaborative Cross (CC) recombinant inbred inter-cross (RIX) mice to elucidate the genetic basis of antipsychotic-induced adverse drug reactions (ADRs). We performed a battery of behavioral tests in 840 mice from 73 RIX lines (derived from 62 CC strains) treated with haloperidol or placebo in order to monitor the development of ADRs. We used linear mixed models to test for strain and treatment effects. We observed highly significant strain effects for almost all behavioral measurements investigated (P < 0.001). Further, we observed strong strain-by-treatment interactions for most phenotypes, particularly for changes in distance traveled, vertical activity, and extrapyramidal symptoms (EPS). Estimates of overall heritability ranged from 0.21 (change in body weight) to 0.4 (VCMs and change in distance traveled) while the portion attributable to the interactions of treatment and strain ranged from 0.01 (for change in body weight) to 0.15 (for change in EPS). Interestingly, close to 30% of RIX mice exhibited VCMs, a sensitivity to haloperidol exposure, approximately similar to the rate of TD in humans chronically exposed to haloperidol. Understanding the genetic basis for the susceptibility to antipsychotic ADRs may be possible in mouse, and extrapolation to humans could lead to safer therapeutic approaches for schizophrenia.

Trappc9 deficiency causes parent-of-origin dependent microcephaly and obesity

Some imprinted genes exhibit parental origin specific expression bias rather than being transcribed exclusively from one copy. The physiological relevance of this remains poorly understood. In an analysis of brain-specific allele-biased expression, we identified that Trappc9, a cellular trafficking factor, was expressed predominantly (~70%) from the maternally inherited allele. Loss-of-function mutations in human TRAPPC9 cause a rare neurodevelopmental syndrome characterized by microcephaly and obesity. By studying Trappc9 null mice we discovered that homozygous mutant mice showed a reduction in brain size, exploratory activity and social memory, as well as a marked increase in body weight. A role for Trappc9 in energy balance was further supported by increased ad libitum food intake in a child with TRAPPC9 deficiency. Strikingly, heterozygous mice lacking the maternal allele (70% reduced expression) had pathology similar to homozygous mutants, whereas mice lacking the paternal allele (30% reduction) were phenotypically normal. Taken together, we conclude that Trappc9 deficient mice recapitulate key pathological features of TRAPPC9 mutations in humans and identify a role for Trappc9 and its imprinting in controlling brain development and metabolism.

Omics Application in Animal Science-A Special Emphasis on Stress Response and Damaging Behaviour in Pigs

Increasing stress resilience of livestock is important for ethical and profitable meat and dairy production. Susceptibility to stress can entail damaging behaviours, a common problem in pig production. Breeding animals with increased stress resilience is difficult for various reasons. First, studies on neuroendocrine and behavioural stress responses in farm animals are scarce, as it is difficult to record adequate phenotypes under field conditions. Second, damaging behaviours and stress susceptibility are complex traits, and their biology is not yet well understood. Dissecting complex traits into biologically better defined, heritable and easily measurable proxy traits and developing biomarkers will facilitate recording these traits in large numbers. High-throughput molecular technologies (“omics”) study the entirety of molecules and their interactions in a single analysis step. They can help to decipher the contributions of different physiological systems and identify candidate molecules that are representative of different physiological pathways. Here, we provide a general overview of different omics approaches and we give examples of how these techniques could be applied to discover biomarkers. We discuss the genetic dissection of the stress response by different omics techniques and we provide examples and outline potential applications of omics tools to understand and prevent outbreaks of damaging behaviours.

The role of polygenic susceptibility to obesity among carriers of pathogenic mutations in MC4R in the UK Biobank population

BACKGROUND

Melanocortin 4 receptor (MC4R) deficiency, caused by mutations in MC4R, is the most common cause of monogenic forms of obesity. However, these mutations have often been identified in small-scale, case-focused studies. Here, we assess the penetrance of previously reported MC4R mutations at a population level. Furthermore, we examine why some carriers of pathogenic mutations remain of normal weight, to gain insight into the mechanisms that control body weight.

METHODS AND FINDINGS

We identified 59 known obesity-increasing mutations in MC4R from the Human Gene Mutation Database (HGMD) and Clinvar. We assessed their penetrance and effect on obesity (body mass index [BMI] ≥ 30 kg/m2) in >450,000 individuals (age 40-69 years) of the UK Biobank, a population-based cohort study. Of these 59 mutations, only 11 had moderate-to-high penetrance and increased the odds of obesity by more than 2-fold. We subsequently focused on these 11 mutations and examined differences between carriers of normal weight and carriers with obesity. Twenty-eight of the 182 carriers of these 11 mutations were of normal weight. Body composition of carriers of normal weight was similar to noncarriers of normal weight, whereas among individuals with obesity, carriers had a somewhat higher BMI than noncarriers (1.44 ± 0.07 standard deviation scores [SDSs] ± standard error [SE] versus 1.29 ± 0.001, P = 0.03), because of greater lean mass (1.44 ± 0.09 versus 1.15 ± 0.002, P = 0.002). Carriers of normal weight more often reported that, already at age 10 years, their body size was below average or average (72%) compared with carriers with obesity (48%) (P = 0.01). To assess the polygenic contribution to body weight in carriers of normal weight and carriers with obesity, we calculated a genome-wide polygenic risk score for BMI (PRSBMI). The PRSBMI of carriers of normal weight (PRSBMI = -0.64 ± 0.18) was significantly lower than of carriers with obesity (0.40 ± 0.11; P = 1.7 × 10-6), and tended to be lower than that of noncarriers of normal weight (-0.29 ± 0.003; P = 0.05). Among carriers, those with a low PRSBMI (bottom quartile) have an approximately 5-kg/m2 lower BMI (approximately 14 kg of body weight for a 1.7-m-tall person) than those with a high PRS (top quartile). Because the UK Biobank population is healthier than the general population in the United Kingdom, penetrance may have been somewhat underestimated.

CONCLUSIONS

We showed that large-scale data are needed to validate the impact of mutations observed in small-scale and case-focused studies. Furthermore, we observed that despite the key role of MC4R in obesity, the effects of pathogenic MC4R mutations may be countered, at least in part, by a low polygenic risk potentially representing other innate mechanisms implicated in body weight regulation.

Allele-specific expression is widespread in Bos indicus muscle and affects meat quality candidate genes

Differences between the expression of the two alleles of a gene are known as allele-specific expression (ASE), a common event in the transcriptome of mammals. Despite ASE being a source of phenotypic variation, its occurrence and effects on genetic prediction of economically relevant traits are still unexplored in bovines. Furthermore, as ASE events are likely driven by cis-regulatory mutations, scanning them throughout the bovine genome represents a significant step to elucidate the mechanisms underlying gene expression regulation. To address this question in a Bos indicus population, we built the ASE profile of the skeletal muscle tissue of 190 Nelore steers, using RNA sequencing data and SNPs genotypes from the Illumina BovineHD BeadChip (770 K bp). After quality control, 820 SNPs showed at least one sample with ASE. These SNPs were widespread among all autosomal chromosomes, being 32.01% found in 3'UTR and 31.41% in coding regions. We observed a considerable variation of ASE profile among individuals, which highlighted the need for biological replicates in ASE studies. Functional analysis revealed that ASE genes play critical biological functions in the development and maintenance of muscle tissue. Additionally, some of these genes were previously reported as associated with beef production and quality traits in livestock, thus indicating a possible source of bias on genomic predictions for these traits.

Cigarette Smoke Exposure and Inflammatory Signaling Increase the Expression of the SARS-CoV-2 Receptor ACE2 in the Respiratory Tract

The factors mediating fatal SARS-CoV-2 infections are poorly understood. Here, we show that cigarette smoke causes a dose-dependent upregulation of angiotensin converting enzyme 2 (ACE2), the SARS-CoV-2 receptor, in rodent and human lungs. Using single-cell sequencing data, we demonstrate that ACE2 is expressed in a subset of secretory cells in the respiratory tract. Chronic smoke exposure triggers the expansion of this cell population and a concomitant increase in ACE2 expression. In contrast, quitting smoking decreases the abundance of these secretory cells and reduces ACE2 levels. Finally, we demonstrate that ACE2 expression is responsive to inflammatory signaling and can be upregulated by viral infections or interferon treatment. Taken together, these results may partially explain why smokers are particularly susceptible to severe SARS-CoV-2 infections. Furthermore, our work identifies ACE2 as an interferon-stimulated gene in lung cells, suggesting that SARS-CoV-2 infections could create positive feedback loops that increase ACE2 levels and facilitate viral dissemination.

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.

Short paired-end reads trump long single-end reads for expression analysis

Background

Typical experimental design advice for expression analyses using RNA-seq generally assumes that single-end reads provide robust gene-level expression estimates in a cost-effective manner, and that the additional benefits obtained from paired-end sequencing are not worth the additional cost. However, in many cases (e.g., with Illumina NextSeq and NovaSeq instruments), shorter paired-end reads and longer single-end reads can be generated for the same cost, and it is not obvious which strategy should be preferred. Using publicly available data, we test whether short-paired end reads can achieve more robust expression estimates and differential expression results than single-end reads of approximately the same total number of sequenced bases.

Results

At both the transcript and gene levels, 2 × 40 paired-end reads unequivocally provide expression estimates that are more highly correlated with 2 × 125 than 1 × 75 reads; in nearly all cases, those correlations are also greater than for 1 × 125, despite the greater total number of sequenced bases for the latter. Across an array of metrics, differential expression tests based upon 2 × 40 consistently outperform those using 1 × 75.

Conclusion

Researchers seeking a cost-effective approach for gene-level expression analysis should prefer short paired-end reads over a longer single-end strategy. Short paired-end reads will also give reasonably robust expression estimates and differential expression results at the isoform level.

Genome-wide analysis of spatiotemporal allele-specific expression in F1 hybrids of meat- and egg-type chickens

In diploid organisms, each gene locus is composed of two parental alleles, which would interact with each other for determining the phenotypic variation. Better understanding of the allele-specific expression (ASE) in farm animals is much important to explore the genetic basis underlying economically important traits, which have been poorly understood yet. In this study, genome-wide analysis was applied to explore the spatiotemporal pattern of ASE in the F1 hybrids of chicken. First, meat- and egg-type chickens were selected for producing a full-sib F1 hybrid population (n = 57). Then, genome resequencing of two parents and 38 offspring were performed and liver and breast muscle samples (n = 38) were subjected to strand-specific RNA sequencing (ssRNA-seq) for ASE detection at 1, 28, and 56 days of age, respectively. The results accurately identified a total of 465 informative genes that could be distinguished with respect to their parental origins. There were 0.4% - 4.1% of informative genes showing ASE, and 57 of them were found across different tissues and time points. Besides, most ASE genes in chickens were tissue-specific, and no matter what the time-point pattern of one ASE gene, the same parental allele of this gene almost showed consistently higher or lower expression across all time points in the same type tissue. In conclusion, this study indicated that most of ASE genes were tissue-specific and time-dependent.

Genetic basis for divergence in developmental gene expression in two closely related sea urchins

The genetic basis for divergence in developmental gene expression among species is poorly understood, despite growing evidence that such changes underlie many interesting traits. Here we quantify transcription in hybrids of Heliocidaris tuberculata and Heliocidaris erythrogramma, two closely related sea urchins with highly divergent developmental gene expression and life histories. We find that most expression differences between species result from genetic influences that affect one stage of development, indicating limited pleiotropic consequences for most mutations that contribute to divergence in gene expression. Activation of zygotic transcription is broadly delayed in H. erythrogramma, the species with the derived life history, despite its overall faster premetamorphic development. Altered expression of several terminal differentiation genes associated with the derived larval morphology of H. erythrogramma is based largely on differences in the expression or function of their upstream regulators, providing insights into the genetic basis for the evolution of key life history traits.

Role of the imprinted allele of the Cdkn1c gene in mouse neocortical development

Imprinted genes are expressed from only one allele in a parent of origin-specific manner. The cyclin-dependent kinase inhibitor p57kip2 is encoded by an imprinted gene Cdkn1c, with the paternal allele being silenced. The possible expression and function of the paternal allele of Cdkn1c have remained little studied, however. We now show that the paternal allele of the Cdkn1c gene is expressed at a low level in the developing mouse neocortex. Surprisingly, the central nervous system-specific conditional deletion of the paternal allele (pat cKO) at the Cdkn1c locus resulted in a marked reduction in brain size. Furthermore, pat cKO gradually reduced the number of neural stem-progenitor cells (NPCs) during neocortical development, and thus reduced the number of upper-layer neurons, which were derived from late-stage NPCs. Our results thus show that the paternal allele of the Cdkn1c locus plays a key role in maintenance of NPCs during neocortical development.

Integrative QTL analysis of gene expression and chromatin accessibility identifies multi-tissue patterns of genetic regulation

Gene transcription profiles across tissues are largely defined by the activity of regulatory elements, most of which correspond to regions of accessible chromatin. Regulatory element activity is in turn modulated by genetic variation, resulting in variable transcription rates across individuals. The interplay of these factors, however, is poorly understood. Here we characterize expression and chromatin state dynamics across three tissues-liver, lung, and kidney-in 47 strains of the Collaborative Cross (CC) mouse population, examining the regulation of these dynamics by expression quantitative trait loci (eQTL) and chromatin QTL (cQTL). QTL whose allelic effects were consistent across tissues were detected for 1,101 genes and 133 chromatin regions. Also detected were eQTL and cQTL whose allelic effects differed across tissues, including local-eQTL for Pik3c2g detected in all three tissues but with distinct allelic effects. Leveraging overlapping measurements of gene expression and chromatin accessibility on the same mice from multiple tissues, we used mediation analysis to identify chromatin and gene expression intermediates of eQTL effects. Based on QTL and mediation analyses over multiple tissues, we propose a causal model for the distal genetic regulation of Akr1e1, a gene involved in glycogen metabolism, through the zinc finger transcription factor Zfp985 and chromatin intermediates. This analysis demonstrates the complexity of transcriptional and chromatin dynamics and their regulation over multiple tissues, as well as the value of the CC and related genetic resource populations for identifying specific regulatory mechanisms within cells and tissues.

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.

Evolution of cis- and trans-regulatory divergence in the chicken genome between two contrasting breeds analyzed using three tissue types at one-day-old

Background

Gene expression variation is a key underlying factor influencing phenotypic variation, and can occur via cis- or trans-regulation. To understand the role of cis- and trans-regulatory variation on population divergence in chicken, we developed reciprocal crosses of two chicken breeds, White Leghorn and Cornish Game, which exhibit major differences in body size and reproductive traits, and used them to determine the degree of cis versus trans variation in the brain, liver, and muscle tissue of male and female 1-day-old specimens.

Results

We provided an overview of how transcriptomes are regulated in hybrid progenies of two contrasting breeds based on allele specific expression analysis. Compared with cis-regulatory divergence, trans-acting genes were more extensive in the chicken genome. In addition, considerable compensatory cis- and trans-regulatory changes exist in the chicken genome. Most importantly, stronger purifying selection was observed on genes regulated by trans-variations than in genes regulated by the cis elements.

Conclusions

We present a pipeline to explore allele-specific expression in hybrid progenies of inbred lines without a specific reference genome. Our research is the first study to describe the regulatory divergence between two contrasting breeds. The results suggest that artificial selection associated with domestication in chicken could have acted more on trans-regulatory divergence than on cis-regulatory divergence.

Fast effect size shrinkage software for beta-binomial models of allelic imbalance

Allelic imbalance occurs when the two alleles of a gene are differentially expressed within a diploid organism, and can indicate important differences in cis-regulation and epigenetic state across the two chromosomes. Because of this, the ability to accurately quantify the proportion at which each allele of a gene is expressed is of great interest to researchers. This becomes challenging in the presence of small read counts and/or sample sizes, which can cause estimates for allelic expression proportions to have high variance. Investigators have traditionally dealt with this problem by filtering out genes with small counts and samples. However, this may inadvertently remove important genes that have truly large allelic imbalances. Another option is to use Bayesian estimators to reduce the variance. To this end, we evaluated the accuracy of three different estimators, the latter two of which are Bayesian shrinkage estimators: maximum likelihood, approximate posterior estimation of GLM coefficients (apeglm) and adaptive shrinkage (ash). We also wrote C++ code to quickly calculate ML and apeglm estimates, and integrated it into the apeglm package. The three methods were evaluated on both simulated and real data. Apeglm consistently performed better than ML according to a variety of criteria, including mean absolute error and concordance at the top. While ash had lower error and greater concordance than ML on the simulations, it also had a tendency to over-shrink large effects, and performed worse on the real data according to error and concordance. Furthermore, when compared to five other packages that also fit beta-binomial models, the apeglm package was substantially faster, making our package useful for quick and reliable analyses of allelic imbalance. Apeglm is available as an R/Bioconductor package at http://bioconductor.org/packages/apeglm.