Differential allelic representation (DAR) identifies candidate eQTLs and improves transcriptome analysis.

In comparisons between mutant and wild-type genotypes, transcriptome analysis can reveal the direct impacts of a mutation, together with the homeostatic responses of the biological system. Recent studies have highlighted that, when the effects of homozygosity for recessive mutations are studied in non-isogenic backgrounds, genes located proximal to the mutation on the same chromosome often appear over-represented among those genes identified as differentially expressed (DE). One hypothesis suggests that DE genes chromosomally linked to a mutation may not reflect functional responses to the mutation but, instead, result from an unequal distribution of expression quantitative trait loci (eQTLs) between sample groups of mutant or wild-type genotypes. This is problematic because eQTL expression differences are difficult to distinguish from genes that are DE due to functional responses to a mutation. Here we show that chromosomally co-located differentially expressed genes (CC-DEGs) are also observed in analyses of dominant mutations in heterozygotes. We define a method and a metric to quantify, in RNA-sequencing data, localised differential allelic representation (DAR) between those sample groups subjected to differential expression analysis. We show how the DAR metric can predict regions prone to eQTL-driven differential expression, and how it can improve functional enrichment analyses through gene exclusion or weighting-based approaches. Advantageously, this improved ability to identify probable eQTLs also reveals examples of CC-DEGs that are likely to be functionally related to a mutant phenotype. This supports a long-standing prediction that selection for advantageous linkage disequilibrium influences chromosome evolution. By comparing the genomes of zebrafish (Danio rerio) and medaka (Oryzias latipes), a teleost with a conserved ancestral karyotype, we find possible examples of chromosomal aggregation of CC-DEGs during evolution of the zebrafish lineage. Our method for DAR analysis requires only RNA-sequencing data, facilitating its application across new and existing datasets.

Genome-Wide CRISPR Screen Identifies RACK1 as a Critical Host Factor for Flavivirus Replication.

Cellular factors have important roles in all facets of the flavivirus replication cycle. Deciphering viral-host protein interactions is essential for understanding the flavivirus life cycle as well as development of effective antiviral strategies. To uncover novel host factors that are co-opted by multiple flaviviruses, a CRISPR/Cas9 genome wide knockout (KO) screen was employed to identify genes required for replication of Zika virus (ZIKV). Receptor for Activated Protein C Kinase 1 (RACK1) was identified as a novel host factor required for ZIKV replication, which was confirmed via complementary experiments. Depletion of RACK1 via siRNA demonstrated that RACK1 is important for replication of a wide range of mosquito- and tick-borne flaviviruses, including West Nile Virus (WNV), Dengue Virus (DENV), Powassan Virus (POWV) and Langat Virus (LGTV) as well as the coronavirus SARS-CoV-2, but not for YFV, EBOV, VSV or HSV. Notably, flavivirus replication was only abrogated when RACK1 expression was dampened prior to infection. Utilising a non-replicative flavivirus model, we show altered morphology of viral replication factories and reduced formation of vesicle packets (VPs) in cells lacking RACK1 expression. In addition, RACK1 interacted with NS1 protein from multiple flaviviruses; a key protein for replication complex formation. Overall, these findings reveal RACK1's crucial role to the biogenesis of pan-flavivirus replication organelles. IMPORTANCE Cellular factors are critical in all facets of viral lifecycles, where overlapping interactions between the virus and host can be exploited as possible avenues for the development of antiviral therapeutics. Using a genome-wide CRISPR knockout screening approach to identify novel cellular factors important for flavivirus replication we identified RACK1 as a pro-viral host factor for both mosquito- and tick-borne flaviviruses in addition to SARS-CoV-2. Using an innovative flavivirus protein expression system, we demonstrate for the first time the impact of the loss of RACK1 on the formation of viral replication factories known as 'vesicle packets' (VPs). In addition, we show that RACK1 can interact with numerous flavivirus NS1 proteins as a potential mechanism by which VP formation can be induced by the former.

Placental Transcription Profiling in 6-23 Weeks' Gestation Reveals Differential Transcript Usage in Early Development.

The human placenta is a rapidly developing transient organ that is key to pregnancy success. Early development of the conceptus occurs in a low oxygen environment before oxygenated maternal blood begins to flow into the placenta at ~10-12 weeks' gestation. This process is likely to substantially affect overall placental gene expression. Transcript variability underlying gene expression has yet to be profiled. In this study, accurate transcript expression profiles were identified for 84 human placental chorionic villus tissue samples collected across 6-23 weeks' gestation. Differential gene expression (DGE), differential transcript expression (DTE) and differential transcript usage (DTU) between 6-10 weeks' and 11-23 weeks' gestation groups were assessed. In total, 229 genes had significant DTE yet no significant DGE. Integration of DGE and DTE analyses found that differential expression patterns of individual transcripts were commonly masked upon aggregation to the gene-level. Of the 611 genes that exhibited DTU, 534 had no significant DGE or DTE. The four most significant DTU genes ADAM10, VMP1, GPR126, and ASAH1, were associated with hypoxia-responsive pathways. Transcript usage is a likely regulatory mechanism in early placentation. Identification of functional roles will facilitate new insight in understanding the origins of pregnancy complications.

Transcriptome analyses of 7-day-old zebrafish larvae possessing a familial Alzheimer's disease-like mutation in psen1 indicate effects on oxidative phosphorylation, ECM and MCM functions, and iron homeostasis.

BACKGROUND: Early-onset familial Alzheimer's disease (EOfAD) is promoted by dominant mutations, enabling the study of Alzheimer's disease (AD) pathogenic mechanisms through generation of EOfAD-like mutations in animal models. In a previous study, we generated an EOfAD-like mutation, psen1Q96_K97del, in zebrafish and performed transcriptome analysis comparing entire brains from 6-month-old wild type and heterozygous mutant fish. We identified predicted effects on mitochondrial function and endolysosomal acidification. Here we aimed to determine whether similar effects occur in 7 day post fertilization (dpf) zebrafish larvae that might be exploited in screening of chemical libraries to find ameliorative drugs. RESULTS: We generated clutches of wild type and heterozygous psen1Q96_K97del 7 dpf larvae using a paired-mating strategy to reduce extraneous genetic variation before performing a comparative transcriptome analysis. We identified 228 differentially expressed genes and performed various bioinformatics analyses to predict cellular functions. CONCLUSIONS: Our analyses predicted a significant effect on oxidative phosphorylation, consistent with our earlier observations of predicted effects on ATP synthesis in adult heterozygous psen1Q96_K97del brains. The dysregulation of minichromosome maintenance protein complex (MCM) genes strongly contributed to predicted effects on DNA replication and the cell cycle and may explain earlier observations of genome instability due to PSEN1 mutation. The upregulation of crystallin gene expression may be a response to defective activity of mutant Psen1 protein in endolysosomal acidification. Genes related to extracellular matrix (ECM) were downregulated, consistent with previous studies of EOfAD mutant iPSC neurons and postmortem late onset AD brains. Also, changes in expression of genes controlling iron ion transport were observed without identifiable changes in the prevalence of transcripts containing iron responsive elements (IREs) in their 3' untranslated regions (UTRs). These changes may, therefore, predispose to the apparent iron dyshomeostasis previously observed in 6-month-old heterozygous psen1Q96_K97del EOfAD-like mutant brains.

ngsReports: a Bioconductor package for managing FastQC reports and other NGS related log files.

MOTIVATION: High throughput next generation sequencing (NGS) has become exceedingly cheap, facilitating studies to be undertaken containing large sample numbers. Quality control (QC) is an essential stage during analytic pipelines and the outputs of popular bioinformatics tools such as FastQC and Picard can provide information on individual samples. Although these tools provide considerable power when carrying out QC, large sample numbers can make inspection of all samples and identification of systemic bias a challenge. RESULTS: We present ngsReports, an R package designed for the management and visualization of NGS reports from within an R environment. The available methods allow direct import into R of FastQC reports along with outputs from other tools. Visualization can be carried out across many samples using default, highly customizable plots with options to perform hierarchical clustering to quickly identify outlier libraries. Moreover, these can be displayed in an interactive shiny app or HTML report for ease of analysis. AVAILABILITY AND IMPLEMENTATION: The ngsReports package is available on Bioconductor and the GUI shiny app is available at https://github.com/UofABioinformaticsHub/shinyNgsreports. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Peptidase inhibitor 16 identifies a human regulatory T-cell subset with reduced FOXP3 expression over the first year of recent onset type 1 diabetes.

CD4+ T-cell subsets play a major role in the host response to infection, and a healthy immune system requires a fine balance between reactivity and tolerance. This balance is in part maintained by regulatory T cells (Treg), which promote tolerance, and loss of immune tolerance contributes to autoimmunity. As the T cells which drive immunity are diverse, identifying and understanding how these subsets function requires specific biomarkers. From a human CD4 Tconv/Treg cell genome wide analysis we identified peptidase inhibitor 16 (PI16) as a CD4 subset biomarker and we now show detailed analysis of its distribution, phenotype and links to Treg function in type 1 diabetes. To determine the clinical relevance of Pi16 Treg, we analysed PI16+ Treg cells from type 1 diabetes patient samples. We observed that FOXP3 expression levels declined with disease progression, suggesting loss of functional fitness in these Treg cells in Type 1 diabetes, and in particular the rate of loss of FOXP3 expression was greatest in the PI16+ve Treg. We propose that PI16 has utility as a biomarker of functional human Treg subsets and may be useful for tracking loss of immune function in vivo. The ability to stratify at risk patients so that tailored interventions can be applied would open the door to personalised medicine for Type 1 diabetes.

Adaptive changes in the genomes of wild rabbits after 16 years of viral epidemics.

Since its introduction to control overabundant invasive European rabbits (Oryctolagus cuniculus), the highly virulent rabbit haemorrhagic disease virus (RHDV) has caused regular annual disease outbreaks in Australian rabbit populations. Although initially reducing rabbit abundance by 60%, continent-wide, experimental evidence has since indicated increased genetic resistance in wild rabbits that have experienced RHDV-driven selection. To identify genetic adaptations, which explain the increased resistance to this biocontrol virus, we investigated genome-wide SNP (single nucleotide polymorphism) allele frequency changes in a South Australian rabbit population that was sampled in 1996 (pre-RHD genomes) and after 16 years of RHDV outbreaks. We identified several SNPs with changed allele frequencies within or close to genes potentially important for increased RHD resistance. The identified genes are known to be involved in virus infections and immune reactions or had previously been identified as being differentially expressed in healthy versus acutely RHDV-infected rabbits. Furthermore, we show in a simulation study that the allele/genotype frequency changes cannot be explained by drift alone and that several candidate genes had also been identified as being associated with surviving RHD in a different Australian rabbit population. Our unique data set allowed us to identify candidate genes for RHDV resistance that have evolved under natural conditions, and over a time span that would not have been feasible in an experimental setting. Moreover, it provides a rare example of host genetic adaptations to virus-driven selection in response to a suddenly emerging infectious disease.

Short-term environmental enrichment, and not physical exercise, alleviate cognitive decline and anxiety from middle age onwards without affecting hippocampal gene expression.

Physical exercise (PE) and environmental enrichment (EE) have consistently been shown to modulate behavior and neurobiological mechanisms. The current literature lacks evidence to confirm the relationship between PE and EE, if any, and whether short-term treatment with PE, EE, or PE+EE could be considered to correct age-related behavioral deficits. Three-, 8-, and 13-month-old C57BL/6 mice were assigned to either PE, EE, or PE+EE treatment groups (n = 12-16/group) for 4 weeks before behavioral testing and were compared to controls. Differential effects of the treatments on various behaviors and hippocampal gene expression were measured using an established behavioral battery and high-throughput qPCR respectively. Short-term EE enhanced locomotor activity at 9 and 14 months of age, whereas the combination of PE and EE reduced locomotor activity in the home cage at 14 months. Short-term EE also was found to reverse the age-related increase in anxiety at 9 months and spatial memory deficits at 14 months of age. Conversely, short-term PE induced spatial learning impairment and depressive-like behavior at four months but showed no effects in 9- and 14-month-old mice. PE and PE+EE, but not EE, modified the expression of several hippocampal genes at 9 months of age compared with control mice. In conclusion, short-term EE may help to alleviate age-related cognitive decline and increase in anxiety, without altering hippocampal gene expression. On the contrary, PE is detrimental at a young age for both affective-like behaviors and spatial learning and memory but showed no effects at middle and late middle age despite hippocampal gene expression alterations.

PCDH19 regulation of neural progenitor cell differentiation suggests asynchrony of neurogenesis as a mechanism contributing to PCDH19 Girls Clustering Epilepsy.

PCDH19-Girls Clustering Epilepsy (PCDH19-GCE) is a childhood epileptic encephalopathy characterised by a spectrum of neurodevelopmental problems. PCDH19-GCE is caused by heterozygous loss-of-function mutations in the X-chromosome gene, Protocadherin 19 (PCDH19) encoding a cell-cell adhesion molecule. Intriguingly, hemizygous males are generally unaffected. As PCDH19 is subjected to random X-inactivation, heterozygous females are comprised of a mosaic of cells expressing either the normal or mutant allele, which is thought to drive pathology. Despite being the second most prevalent monogeneic cause of epilepsy, little is known about the role of PCDH19 in brain development. In this study we show that PCDH19 is highly expressed in human neural stem and progenitor cells (NSPCs) and investigate its function in vitro in these cells of both mouse and human origin. Transcriptomic analysis of mouse NSPCs lacking Pcdh19 revealed changes to genes involved in regulation of neuronal differentiation, and we subsequently show that loss of Pcdh19 causes increased NSPC neurogenesis. We reprogramed human fibroblast cells harbouring a pathogenic PCDH19 mutation into human induced pluripotent stem cells (hiPSC) and employed neural differentiation of these to extend our studies into human NSPCs. As in mouse, loss of PCDH19 function caused increased neurogenesis, and furthermore, we show this is associated with a loss of human NSPC polarity. Overall our data suggests a conserved role for PCDH19 in regulating mammalian cortical neurogenesis and has implications for the pathogenesis of PCDH19-GCE. We propose that the difference in timing or "heterochrony" of neuronal cell production originating from PCDH19 wildtype and mutant NSPCs within the same individual may lead to downstream asynchronies and abnormalities in neuronal network formation, which in-part predispose the individual to network dysfunction and epileptic activity.

Resistance to RHD virus in wild Australian rabbits: Comparison of susceptible and resistant individuals using a genomewide approach.

Deciphering the genes involved in disease resistance is essential if we are to understand host-pathogen coevolutionary processes. The rabbit haemorrhagic disease virus (RHDV) was imported into Australia in 1995 as a biocontrol agent to manage one of the most successful and devastating invasive species, the European rabbit (Oryctolagus cuniculus). During the first outbreaks of the disease, RHDV caused mortality rates of up to 97%. Recently, however, increased genetic resistance to RHDV has been reported. Here, we have aimed to identify genomic differences between rabbits that survived a natural infection with RHDV and those that died in the field using a genomewide next-generation sequencing (NGS) approach. We detected 72 SNPs corresponding to 133 genes associated with survival of a RHD infection. Most of the identified genes have known functions in virus infections and replication, immune responses or apoptosis, or have previously been found to be regulated during RHD. Some of the genes identified in experimental studies, however, did not seem to play a role under natural selection regimes, highlighting the importance of field studies to complement the genomic background of wildlife diseases. Our study provides a set of candidate markers as a tool for the future scanning of wild rabbits for their resistance to RHDV. This is important both for wild rabbit populations in southern Europe where RHD is regarded as a serious problem decimating the prey of endangered predator species and for assessing the success of currently planned RHDV variant biocontrol releases in Australia.

The androgen receptor interacts with GATA3 to transcriptionally regulate a luminal epithelial cell phenotype in breast cancer.

BACKGROUND: The androgen receptor (AR) is a tumor suppressor in estrogen receptor (ER) positive breast cancer, a role sustained in some ER negative breast cancers. Key factors dictating AR genomic activity in a breast context are largely unknown. Herein, we employ an unbiased chromatin immunoprecipitation-based proteomic technique to identify endogenous AR interacting co-regulatory proteins in ER positive and negative models of breast cancer to gain new insight into mechanisms of AR signaling in this disease. RESULTS: The DNA-binding factor GATA3 is identified and validated as a novel AR interacting protein in breast cancer cells irrespective of ER status. AR activation by the natural ligand 5α-dihydrotestosterone (DHT) increases nuclear AR-GATA3 interactions, resulting in AR-dependent enrichment of GATA3 chromatin binding at a sub-set of genomic loci. Silencing GATA3 reduces but does not prevent AR DNA binding and transactivation of genes associated with AR/GATA3 co-occupied loci, indicating a co-regulatory role for GATA3 in AR signaling. DHT-induced AR/GATA3 binding coincides with upregulation of luminal differentiation genes, including EHF and KDM4B, established master regulators of a breast epithelial cell lineage. These findings are validated in a patient-derived xenograft model of breast cancer. Interaction between AR and GATA3 is also associated with AR-mediated growth inhibition in ER positive and ER negative breast cancer. CONCLUSIONS: AR and GATA3 interact to transcriptionally regulate luminal epithelial cell differentiation in breast cancer regardless of ER status. This interaction facilitates the tumor suppressor function of AR and mechanistically explains why AR expression is associated with less proliferative, more differentiated breast tumors and better overall survival in breast cancer.

Differential allelic representation (DAR) identifies candidate eQTLs and improves transcriptome analysis.

In comparisons between mutant and wild-type genotypes, transcriptome analysis can reveal the direct impacts of a mutation, together with the homeostatic responses of the biological system. Recent studies have highlighted that, when homozygous mutations are studied in non-isogenic backgrounds, genes from the same chromosome as a mutation often appear over-represented among differentially expressed (DE) genes. One hypothesis suggests that DE genes chromosomally linked to a mutation may not reflect true biological responses to the mutation but, instead, result from differences in representation of expression quantitative trait loci (eQTLs) between sample groups selected on the basis of mutant or wild-type genotype. This is problematic when inclusion of spurious DE genes in a functional enrichment study results in incorrect inferences of mutation effect. Here we show that chromosomally co-located differentially expressed genes (CC-DEGs) can also be observed in analyses of dominant mutations in heterozygotes. We define a method and a metric to quantify, in RNA-sequencing data, localised differential allelic representation (DAR) between groups of samples subject to differential expression analysis. We show how the DAR metric can predict regions prone to eQTL-driven differential expression, and how it can improve functional enrichment analyses through gene exclusion or weighting of gene-level rankings. Advantageously, this improved ability to identify probable eQTLs also reveals examples of CC-DEGs that are likely to be functionally related to a mutant phenotype. This supports a long-standing prediction that selection for advantageous linkage disequilibrium influences chromosome evolution. By comparing the genomes of zebrafish (Danio rerio) and medaka (Oryzias latipes), a teleost with a conserved ancestral karyotype, we find possible examples of chromosomal aggregation of CC-DEGs during evolution of the zebrafish lineage. The DAR metric provides a solid foundation for addressing the eQTL issue in new and existing datasets because it relies solely on RNA-sequencing data.

The Use of Zebrafish in Transcriptome Analysis of the Early Effects of Mutations Causing Early Onset Familial Alzheimer's Disease and Other Inherited Neurodegenerative Conditions.

The degree to which non-human animals can be used to model Alzheimer's disease is a contentious issue, particularly as there is still widespread disagreement regarding the pathogenesis of this neurodegenerative dementia. The currently popular transgenic models are based on artificial expression of genes mutated in early onset forms of familial Alzheimer's disease (EOfAD). Uncertainty regarding the veracity of these models led us to focus on heterozygous, single mutations of endogenous genes (knock-in models) as these most closely resemble the genetic state of humans with EOfAD, and so incorporate the fewest assumptions regarding pathological mechanism. We have generated a number of lines of zebrafish bearing EOfAD-like and non-EOfAD-like mutations in genes equivalent to human PSEN1, PSEN2, and SORL1. To analyze the young adult brain transcriptomes of these mutants, we exploited the ability of zebrafish to produce very large families of simultaneous siblings composed of a variety of genotypes and raised in a uniform environment. This "intra-family" analysis strategy greatly reduced genetic and environmental "noise" thereby allowing detection of subtle changes in gene sets after bulk RNA sequencing of entire brains. Changes to oxidative phosphorylation were predicted for all EOfAD-like mutations in the three genes studied. Here we describe some of the analytical lessons learned in our program combining zebrafish genome editing with transcriptomics to understand the molecular pathologies of neurodegenerative disease.

Unbiased gene expression analysis of the delayed fracture healing observed in Zucker diabetic fatty rats.

Aims: Impaired fracture repair in patients with type 2 diabetes mellitus (T2DM) is not fully understood. In this study, we aimed to characterize the local changes in gene expression (GE) associated with diabetic fracture. We used an unbiased approach to compare GE in the fracture callus of Zucker diabetic fatty (ZDF) rats relative to wild-type (WT) littermates at three weeks following femoral osteotomy. Methods: Zucker rats, WT and homozygous for leptin receptor mutation (ZDF), were fed a moderately high-fat diet to induce T2DM only in the ZDF animals. At ten weeks of age, open femoral fractures were simulated using a unilateral osteotomy stabilized with an external fixator. At three weeks post-surgery, the fractured femur from each animal was retrieved for analysis. Callus formation and the extent of healing were assessed by radiograph and histology. Bone tissue was processed for total RNA extraction and messenger RNA (mRNA) sequencing (mRNA-Seq). Results: Radiographs and histology demonstrated impaired fracture healing in ZDF rats with incomplete bony bridge formation and an influx of intramedullary inflammatory tissue. In comparison, near-complete bridging between cortices was observed in Sham WT animals. Of 13,160 genes, mRNA-Seq analysis identified 13 that were differentially expressed in ZDF rat callus, using a false discovery rate (FDR) threshold of 10%. Seven genes were upregulated with high confidence (FDR = 0.05) in ZDF fracture callus, most with known roles in inflammation. Conclusion: These findings suggest that elevated or prolonged inflammation contributes to delayed fracture healing in T2DM. The identified genes may be used as biomarkers to monitor and treat delayed fracture healing in diabetic patients.

Determination of Tr1 cell populations correlating with distinct activation states in acute IAV infection.

Type I regulatory (Tr1) cells are defined as FOXP3-IL-10-secreting clusters of differentiation (CD4+) T cells that contribute to immune suppression and typically express the markers LAG-3 and CD49b and other co-inhibitory receptors. These cells have not been studied in detail in the context of the resolution of acute infection in the lung. Here, we identify FOXP3- interleukin (IL)-10+ CD4+ T cells transiently accumulating in the lung parenchyma during resolution of the response to sublethal influenza A virus (IAV) infection in mice. These cells were dependent on IL-27Rα, which was required for timely recovery from IAV-induced weight loss. LAG-3 and CD49b were not generally co-expressed by FOXP3- IL-10+ CD4+ T cells in this model and four populations of these cells based on LAG-3 and CD49b co-expression were apparent [LAG-3-CD49b- (double negative), LAG-3+CD49b+ (double positive), LAG-3+CD49b- (LAG-3+), LAG-3-CD49b+ (CD49b+)]. However, each population exhibited suppressive potential consistent with the definition of Tr1 cells. Notably, differences between these populations of Tr1 cells were apparent including differential dependence on IL-10 to mediate suppression and expression of markers indicative of different activation states and terminal differentiation. Sort-transfer experiments indicated that LAG-3+ Tr1 cells exhibited the capacity to convert to double negative and double positive Tr1 cells, indicative of plasticity between these populations. Together, these data determine the features and suppressive potential of Tr1 cells in the resolution of IAV infection and identify four populations delineated by LAG-3 and CD49b, which likely correspond to different Tr1 cell activation states.

Parallel recovery of chromatin accessibility and gene expression dynamics from frozen human regulatory T cells.

Epigenetic features such as DNA accessibility dictate transcriptional regulation in a cell type- and cell state- specific manner, and mapping this in health vs. disease in clinically relevant material is opening the door to new mechanistic insights and new targets for therapy. Assay for Transposase Accessible Chromatin Sequencing (ATAC-seq) allows chromatin accessibility profiling from low cell input, making it tractable on rare cell populations, such as regulatory T (Treg) cells. However, little is known about the compatibility of the assay with cryopreserved rare cell populations. Here we demonstrate the robustness of an ATAC-seq protocol comparing primary Treg cells recovered from fresh or cryopreserved PBMC samples, in the steady state and in response to stimulation. We extend this method to explore the feasibility of conducting simultaneous quantitation of chromatin accessibility and transcriptome from a single aliquot of 50,000 cryopreserved Treg cells. Profiling of chromatin accessibility and gene expression in parallel within the same pool of cells controls for cellular heterogeneity and is particularly beneficial when constrained by limited input material. Overall, we observed a high correlation of accessibility patterns and transcription factor dynamics between fresh and cryopreserved samples. Furthermore, highly similar transcriptomic profiles were obtained from whole cells and from the supernatants recovered from ATAC-seq reactions. We highlight the feasibility of applying these techniques to profile the epigenomic landscape of cells recovered from cryopreservation biorepositories.

Genome-wide identification of human FOXP3 target genes in natural regulatory T cells.

The transcription factor FOXP3 is essential for the formation and function of regulatory T cells (Tregs), and Tregs are essential for maintaining immune homeostasis and tolerance. This is demonstrated by a lethal autoimmune defect in mice lacking Foxp3 and in immunodysregulation polyendocrinopathy enteropathy X-linked syndrome patients. However, little is known about the molecular basis of human FOXP3 function or the relationship between direct and indirect targets of FOXP3 in human Tregs. To investigate this, we have performed a comprehensive genome-wide analysis for human FOXP3 target genes from cord blood Tregs using chromatin immunoprecipitation array profiling and expression profiling. We have identified 5579 human FOXP3 target genes and derived a core Treg gene signature conserved across species using mouse chromatin immunoprecipitation data sets. A total of 739 of the 5579 FOXP3 target genes were differentially regulated in Tregs compared with Th cells, thus allowing the identification of a number of pathways and biological functions overrepresented in Tregs. We have identified gene families including cell surface molecules and microRNAs that are differentially expressed in FOXP3(+) Tregs. In particular, we have identified a novel role for peptidase inhibitor 16, which is expressed on the cell surface of >80% of resting human CD25(+)FOXP3(+) Tregs, suggesting that in conjunction with CD25 peptidase inhibitor 16 may be a surrogate surface marker for Tregs with potential clinical application.

3DFAACTS-SNP: using regulatory T cell-specific epigenomics data to uncover candidate mechanisms of type 1 diabetes (T1D) risk.

BACKGROUND: Genome-wide association studies (GWAS) have enabled the discovery of single nucleotide polymorphisms (SNPs) that are significantly associated with many autoimmune diseases including type 1 diabetes (T1D). However, many of the identified variants lie in non-coding regions, limiting the identification of mechanisms that contribute to autoimmune disease progression. To address this problem, we developed a variant filtering workflow called 3DFAACTS-SNP to link genetic variants to target genes in a cell-specific manner. Here, we use 3DFAACTS-SNP to identify candidate SNPs and target genes associated with the loss of immune tolerance in regulatory T cells (Treg) in T1D. RESULTS: Using 3DFAACTS-SNP, we identified from a list of 1228 previously fine-mapped variants, 36 SNPs with plausible Treg-specific mechanisms of action. The integration of cell type-specific chromosome conformation capture data in 3DFAACTS-SNP identified 266 regulatory regions and 47 candidate target genes that interact with these variant-containing regions in Treg cells. We further demonstrated the utility of the workflow by applying it to three other SNP autoimmune datasets, identifying 16 Treg-centric candidate variants and 60 interacting genes. Finally, we demonstrate the broad utility of 3DFAACTS-SNP for functional annotation of all known common (> 10% allele frequency) variants from the Genome Aggregation Database (gnomAD). We identified 9376 candidate variants and 4968 candidate target genes, generating a list of potential sites for future T1D or other autoimmune disease research. CONCLUSIONS: We demonstrate that it is possible to further prioritise variants that contribute to T1D based on regulatory function, and illustrate the power of using cell type-specific multi-omics datasets to determine disease mechanisms. Our workflow can be customised to any cell type for which the individual datasets for functional annotation have been generated, giving broad applicability and utility.

Transcriptional targets of senataxin and E2 promoter binding factors are associated with neuro-degenerative pathways during increased autophagic flux.

Autophagy is an intracellular recycling process that degrades harmful molecules and enables survival during starvation, with implications for diseases including dementia, cancer and atherosclerosis. Previous studies demonstrate how a limited number of transcription factors (TFs) can increase autophagy. However, this knowledge has not resulted in translation into therapy, thus, to gain understanding of more suitable targets, we utilized a systems biology approach. We induced autophagy by amino acid starvation and mTOR inhibition in HeLa, HEK 293 and SH-SY5Y cells and measured temporal gene expression using RNA-seq. We observed 456 differentially expressed genes due to starvation and 285 genes due to mTOR inhibition (PFDR < 0.05 in every cell line). Pathway analyses implicated Alzheimer's and Parkinson's diseases (PFDR ≤ 0.024 in SH-SY5Y and HeLa) and amyotrophic lateral sclerosis (ALS, PFDR < 0.05 in mTOR inhibition experiments). Differential expression of the Senataxin (SETX) target gene set was predicted to activate multiple neurodegenerative pathways (PFDR ≤ 0.04). In the SH-SY5Y cells of neuronal origin, the E2F transcription family was predicted to activate Alzheimer's disease pathway (PFDR ≤ 0.0065). These exploratory analyses suggest that SETX and E2F may mediate transcriptional regulation of autophagy and further investigations into their possible role in neuro-degeneration are warranted.

Transient Polycomb activity represses developmental genes in growing oocytes.

BACKGROUND: Non-genetic disease inheritance and offspring phenotype are substantially influenced by germline epigenetic programming, including genomic imprinting. Loss of Polycomb Repressive Complex 2 (PRC2) function in oocytes causes non-genetically inherited effects on offspring, including embryonic growth restriction followed by post-natal offspring overgrowth. While PRC2-dependent non-canonical imprinting is likely to contribute, less is known about germline epigenetic programming of non-imprinted genes during oocyte growth. In addition, de novo germline mutations in genes encoding PRC2 lead to overgrowth syndromes in human patients, but the extent to which PRC2 activity is conserved in human oocytes is poorly understood. RESULTS: In this study, we identify a discrete period of early oocyte growth during which PRC2 is expressed in mouse growing oocytes. Deletion of Eed during this window led to the de-repression of 343 genes. A high proportion of these were developmental regulators, and the vast majority were not imprinted genes. Many of the de-repressed genes were also marked by the PRC2-dependent epigenetic modification histone 3 lysine 27 trimethylation (H3K27me3) in primary-secondary mouse oocytes, at a time concurrent with PRC2 expression. In addition, we found H3K27me3 was also enriched on many of these genes by the germinal vesicle (GV) stage in human oocytes, strongly indicating that this PRC2 function is conserved in the human germline. However, while the 343 genes were de-repressed in mouse oocytes lacking EED, they were not de-repressed in pre-implantation embryos and lost H3K27me3 during pre-implantation development. This implies that H3K27me3 is a transient feature that represses a wide range of genes in oocytes. CONCLUSIONS: Together, these data indicate that EED has spatially and temporally distinct functions in the female germline to repress a wide range of developmentally important genes and that this activity is conserved in the mouse and human germlines.