Stephen Pederson

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.

FOXP3 and FOXP3-regulated microRNAs suppress SATB1 in breast cancer cells

The transcription factor FOXP3 has been identified as a tumour suppressor in the breast and prostate epithelia, but little is known about its specific mechanism of action. We have identified a feed-forward regulatory loop in which FOXP3 suppresses the expression of the oncogene SATB1. In particular, we demonstrate that SATB1 is not only a direct target of FOXP3 repression, but that FOXP3 also induces two miRs, miR-7 and miR-155, which specifically target the 3′-UTR of SATB1 to further regulate its expression. We conclude that FOXP3-regulated miRs form part of the mechanism by which FOXP3 prevents the transformation of the healthy breast epithelium to a cancerous phenotype. Approaches aimed at restoring FOXP3 function and the miRs it regulates could help provide new approaches to target breast cancer.

Unravelling the molecular basis for regulatory T-cell plasticity and loss of function in disease

Regulatory T cells (Treg) are critical for preventing autoimmunity and curtailing responses of conventional effector T cells (Tconv). The reprogramming of T‐cell fate and function to generate Treg requires switching on and off of key gene regulatory networks, which may be initiated by a subtle shift in expression levels of specific genes. This can be achieved by intermediary regulatory processes that include microRNA and long noncoding RNA‐based regulation of gene expression. There are well‐documented microRNA profiles in Treg and Tconv, and these can operate to either reinforce or reduce expression of a specific set of target genes, including FOXP3 itself. This type of feedforward/feedback regulatory loop is normally stable in the steady state, but can alter in response to local cues or genetic risk. This may go some way to explaining T‐cell plasticity. In addition, in chronic inflammation or autoimmunity, altered Treg/Tconv function may be influenced by changes in enhancer–promoter interactions, which are highly cell type‐specific. These interactions are impacted by genetic risk based on genome‐wide association studies and may cause subtle alterations to the gene regulatory networks controlled by or controlling FOXP3 and its target genes. Recent insights into the 3D organisation of chromatin and the mapping of noncoding regulatory regions to the genes they control are shedding new light on the direct impact of genetic risk on T‐cell function and susceptibility to inflammatory and autoimmune conditions.

Resistance to RHD virus in wild Australian rabbits: comparison of susceptible and resistant individuals using a genome-wide 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.

Accelerated brain aging towards transcriptional inversion in a zebrafish model of familial Alzheimer's disease

Alzheimer’s disease (AD) develops silently over decades. We cannot easily access and analyse pre-symptomatic brains, so the earliest molecular changes that initiate AD remain unclear. Previously, we demonstrated that the genes mutated in early-onset, dominantly-inherited familial forms of AD (fAD) are evolving particularly rapidly in mice and rats. Fortunately, some non-mammalian vertebrates such as the zebrafish preserve fAD-relevant transcript isoforms of the PRESENILIN (PSEN1 and PSEN2) genes that these rodents have lost. Zebrafish are powerful vertebrate genetic models for many human diseases, but no genetic model of fAD in zebrafish currently exists. We edited the zebrafish genome to model the unique, protein-truncating fAD mutation of human PSEN2, K115fs. Analysing the brain transcriptome and proteome of young (6-month-old) and aged, infertile (24-month-old) wild type and heterozygous fAD-like mutant female sibling zebrafish supports accelerated brain aging and increased glucocorticoid signalling in young fAD-like fish, leading to a transcriptional ‘inversion’ into glucocorticoid resistance and vast changes in biological pathways in aged, infertile fAD-like fish. Notably, one of these changes involving microglia-associated immune responses regulated by the ETS transcription factor family is preserved between our zebrafish fAD model and human early-onset AD. Importantly, these changes occur before obvious histopathology and likely in the absence of Aβ. Our results support the contributions of early metabolic and oxidative stresses to immune and stress responses favouring AD pathogenesis and highlight the value of our fAD-like zebrafish genetic model for elucidating early changes in the brain that promote AD pathogenesis. The success of our approach has important implications for future modelling of AD.

Genome-Wide SNP Discovery In Field And Laboratory Colonies Of Australian Plutella Species

Understanding dispersal and gene flow is an important focus of evolutionary biology, conservation biology and pest management. The diamondback moth, Plutella xylostella, is a worldwide pest of Brassica vegetable and oilseed cropping systems. This insect has high dispersal ability, which has important consequences for population dynamics and the potential spread of insecticide resistance genes. Population genetic studies of the diamondback moth have found little evidence of population structure, suggesting that frequent intermixing occurs within regions, however the patterns of local and regional dispersal remain to be identified. For this and many other pest species, understanding dispersal is crucial for developing integrated management tactics such as forecasting systems and insecticide resistance management plans. In recent years, next generation sequencing (NGS) methods have provided previously unparalleled resolution for population genetic studies in a wide range of species. Here, we assessed the potential of NGS-derived molecular markers to provide new insights about population structure in the diamondback moth. We use restriction-site-associated DNA sequencing (RAD-Seq) to discover hundreds to thousands of single nucleotide polymorphism (SNP) markers in nine field and laboratory-reared populations collected from Australia. Genotypic data from RAD-Seq markers identified a cryptic species, P. australiana, among individuals collected from a wild host, Diplotaxis sp., indicating strong divergence in the nuclear genomes of two Australian Plutella lineages. Significant genetic differentiation was detected among populations of P. xylostella used in our study, however this could be explained by reduced heterozogosity and genetic drift in laboratory-reared populations founded by relatively few individuals. This study demonstrates that RAD-Seq is a powerful method for generating SNP markers for population genetic studies in this species.

ngsReports: An R 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 can be found in the outputs of popular bioinformatics tools such as FastQC and Picard. Although these tools provide considerable power when carrying out QC, large sample numbers can make 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 output as well as that from aligners such as HISAT2, STAR and Bowtie2. Visualization can be carried out across many samples using heatmaps rendered using ggplot2 and plotly. Moreover, these can be displayed in an interactive shiny app or a HTML report. We also provide methods to assess observed GC content in an organism dependent manner for both transcriptomic and genomic datasets. Importantly, hierarchical clustering can be carried out on heatmaps with large sample sizes to quickly identify outliers and batch effects. Availability and Implementation ngsReports is available at https://github.com/UofABioinformaticsHub/ngsReports.

FOXP3 and miR-155 cooperate to control the invasive potential of human breast cancer cells by down regulating ZEB2 independently of ZEB1

Control of oncogenes, including ZEB1 and ZEB2, is a major checkpoint for preventing cancer, and loss of this control contributes to many cancers, including breast cancer. Thus tumour suppressors, such as FOXP3, which is mutated or lost in many cancer tissues, play an important role in maintaining normal tissue homeostasis. Here we show for the first time that ZEB2 is selectively down regulated by FOXP3 and also by the FOXP3 induced microRNA, miR-155. Interestingly, neither FOXP3 nor miR-155 directly altered the expression of ZEB1. In breast cancer cells repression of ZEB2, independently of ZEB1, resulted in reduced expression of a mesenchymal marker, Vimentin and reduced invasion. However, there was no de-repression of E-cadherin and migration was enhanced. Small interfering RNAs targeting ZEB2 suggest that this was a direct effect of ZEB2 and not FOXP3/miR-155. In normal human mammary epithelial cells, depletion of endogenous FOXP3 resulted in de-repression of ZEB2, accompanied by upregulated expression of vimentin, increased E-cadherin expression and cell morphological changes. We suggest that FOXP3 may help maintain normal breast epithelial characteristics through regulation of ZEB2, and loss of FOXP3 in breast cancer cells results in deregulation of ZEB2.

Transposable elements and gene expression during the evolution of amniotes

BackgroundTransposable elements (TEs) are primarily responsible for the DNA losses and gains in genome sequences that occur over time within and between species. TEs themselves evolve, with clade specific LTR/ERV, LINEs and SINEs responsible for the bulk of species-specific genomic features. Because TEs can contain regulatory motifs, they can be exapted as regulators of gene expression. While TE insertions can provide evolutionary novelty for the regulation of gene expression, their overall impact on the evolution of gene expression is unclear. Previous investigators have shown that tissue specific gene expression in amniotes is more similar across species than within species, supporting the existence of conserved developmental gene regulation. In order to understand how species-specific TE insertions might affect the evolution/conservation of gene expression, we have looked at the association of gene expression in six tissues with TE insertions in six representative amniote genomes.ResultsA novel bootstrapping approach has been used to minimise the conflation of effects of repeat types on gene expression. We compared the expression of orthologs containing recent TE insertions to orthologs that contained older TE insertions, and the expression of non-orthologs containing recent TE insertions to non-orthologs with older TE insertions. Both orthologs and non-orthologs showed significant differences in gene expression associated with TE insertions. TEs were found associated with species-specific changes in gene expression, and the magnitude and direction of expression changes were noteworthy. Overall, orthologs containing species-specific TEs were associated with lower gene expression, while in non-orthologs, non-species specific TEs were associated with higher gene expression. Exceptions were SINE elements in human and chicken, which had an opposite association with gene expression compared to other species.ConclusionsOur observed species-specific associations of TEs with gene expression support a role for TEs in speciation/response to selection by species. TEs do not exhibit consistent associations with gene expression and observed associations can vary depending on the age of TE insertions. Based on these observations, it would be prudent to refrain from extrapolating these and previously reported associations to distantly related species.

Specific growth conditions induce a Streptococcus pneumoniae non-mucoidal, small colony variant and determine the outcome of its co-culture with Haemophilus influenzae

Haemophilus influenzae and Streptococcus pneumoniae are known aetiologic agents of chronic otitis media, frequently as a multispecies infection. In this study, we show that the outcome of H. influenzae/S. pneumoniae interactions is dependent on the nutrient source. In continuous culture containing chemically defined media with lactose, S. pneumoniae was non-viable in mono-culture, and in co-culture remained non-viable until 288 h. With glucose, S. pneumoniae became non-viable in mono-culture, but uniquely existed in 3 distinct states in co-culture: parental cells (until 24 h), a dormant state until 336 h and its re-emergence as a non-mucoidal, small colony variant (SCV). The S. pneumoniae SCV was stable and whole genome sequencing showed three major single nucleotide polymorphisms in the SCV cells-cap3A (capsule biosynthesis pathway), fpg (DNA glycosylase of the DNA repair mechanism) and glutamate-5-kinase. Previously, fpg mutants have shown increased mutator rates, permitting bacterial survival against host-generated stresses. Transcriptomics showed these SCV cells up-regulated sugar transporters and toxin/antitoxin systems. An animal model revealed a reduced survival in the lungs and ear by SCV cells. This is the first study documenting the effect of carbon source and the development of a distinct S. pneumoniae cell type during H. influenzae/S. pneumoniae interactions.