Jennifer M. Taylor

QuantiSNP: an Objective Bayes Hidden-Markov Model to detect and accurately map copy number variation using SNP genotyping data

Array-based technologies have been used to detect chromosomal copy number changes (aneuploidies) in the human genome. Recent studies identified numerous copy number variants (CNV) and some are common polymorphisms that may contribute to disease susceptibility. We developed, and experimentally validated, a novel computational framework (QuantiSNP) for detecting regions of copy number variation from BeadArray™ SNP genotyping data using an Objective Bayes Hidden-Markov Model (OB-HMM). Objective Bayes measures are used to set certain hyperparameters in the priors using a novel re-sampling framework to calibrate the model to a fixed Type I (false positive) error rate. Other parameters are set via maximum marginal likelihood to prior training data of known structure. QuantiSNP provides probabilistic quantification of state classifications and significantly improves the accuracy of segmental aneuploidy identification and mapping, relative to existing analytical tools (Beadstudio, Illumina), as demonstrated by validation of breakpoint boundaries. QuantiSNP identified both novel and validated CNVs. QuantiSNP was developed using BeadArray™ SNP data but it can be adapted to other platforms and we believe that the OB-HMM framework has widespread applicability in genomic research. In conclusion, QuantiSNP is a novel algorithm for high-resolution CNV/aneuploidy detection with application to clinical genetics, cancer and disease association studies.

Global microRNA expression profiles in insulin target tissues in a spontaneous rat model of type 2 diabetes

Aims/hypothesisMicroRNAs regulate a broad range of biological mechanisms. To investigate the relationship between microRNA expression and type 2 diabetes, we compared global microRNA expression in insulin target tissues from three inbred rat strains that differ in diabetes susceptibility.MethodsUsing microarrays, we measured the expression of 283 microRNAs in adipose, liver and muscle tissue from hyperglycaemic (Goto–Kakizaki), intermediate glycaemic (Wistar Kyoto) and normoglycaemic (Brown Norway) rats (n = 5 for each strain). Expression was compared across strains and validated using quantitative RT-PCR. Furthermore, microRNA expression variation in adipose tissue was investigated in 3T3-L1 adipocytes exposed to hyperglycaemic conditions.ResultsWe found 29 significantly differentiated microRNAs (padjusted < 0.05): nine in adipose tissue, 18 in liver and two in muscle. Of these, five microRNAs had expression patterns that correlated with the strain-specific glycaemic phenotype. MiR-222 (padjusted = 0.0005) and miR-27a (padjusted = 0.006) were upregulated in adipose tissue; miR-195 (padjusted = 0.006) and miR-103 (padjusted = 0.04) were upregulated in liver; and miR-10b (padjusted = 0.004) was downregulated in muscle. Exposure of 3T3-L1 adipocytes to increased glucose concentration upregulated the expression of miR-222 (p = 0.008), miR-27a (p = 0.02) and the previously reported miR-29a (p = 0.02). Predicted target genes of these differentially expressed microRNAs are involved in pathways relevant to type 2 diabetes.ConclusionThe expression patterns of miR-222, miR-27a, miR-195, miR-103 and miR-10b varied with hyperglycaemia, suggesting a role for these microRNAs in the pathophysiology of type 2 diabetes, as modelled by the Gyoto–Kakizaki rat. We observed similar patterns of expression of miR-222, miR-27a and miR-29a in adipocytes as a response to increased glucose levels, which supports our hypothesis that altered expression of microRNAs accompanies primary events related to the pathogenesis of type 2 diabetes.

High-Throughput Analysis of Promoter Occupancy Reveals Direct Neural Targets of FOXP2, a Gene Mutated in Speech and Language Disorders

We previously discovered that mutations of the human FOXP2 gene cause a monogenic communication disorder, primarily characterized by difficulties in learning to make coordinated sequences of articulatory gestures that underlie speech. Affected people have deficits in expressive and receptive linguistic processing and display structural and/or functional abnormalities in cortical and subcortical brain regions. FOXP2 provides a unique window into neural processes involved in speech and language. In particular, its role as a transcription factor gene offers powerful functional genomic routes for dissecting critical neurogenetic mechanisms. Here, we employ chromatin immunoprecipitation coupled with promoter microarrays (ChIP-chip) to successfully identify genomic sites that are directly bound by FOXP2 protein in native chromatin of human neuron-like cells. We focus on a subset of downstream targets identified by this approach, showing that altered FOXP2 levels yield significant changes in expression in our cell-based models and that FOXP2 binds in a specific manner to consensus sites within the relevant promoters. Moreover, we demonstrate significant quantitative differences in target expression in embryonic brains of mutant mice, mediated by specific in vivo Foxp2-chromatin interactions. This work represents the first identification and in vivo verification of neural targets regulated by FOXP2. Our data indicate that FOXP2 has dual functionality, acting to either repress or activate gene expression at occupied promoters. The identified targets suggest roles in modulating synaptic plasticity, neurodevelopment, neurotransmission, and axon guidance and represent novel entry points into in vivo pathways that may be disturbed in speech and language disorders.

Foxp2 Regulates Gene Networks Implicated in Neurite Outgrowth in the Developing Brain

Forkhead-box protein P2 is a transcription factor that has been associated with intriguing aspects of cognitive function in humans, non-human mammals, and song-learning birds. Heterozygous mutations of the human FOXP2 gene cause a monogenic speech and language disorder. Reduced functional dosage of the mouse version (Foxp2) causes deficient cortico-striatal synaptic plasticity and impairs motor-skill learning. Moreover, the songbird orthologue appears critically important for vocal learning. Across diverse vertebrate species, this well-conserved transcription factor is highly expressed in the developing and adult central nervous system. Very little is known about the mechanisms regulated by Foxp2 during brain development. We used an integrated functional genomics strategy to robustly define Foxp2-dependent pathways, both direct and indirect targets, in the embryonic brain. Specifically, we performed genome-wide in vivo ChIP–chip screens for Foxp2-binding and thereby identified a set of 264 high-confidence neural targets under strict, empirically derived significance thresholds. The findings, coupled to expression profiling and in situ hybridization of brain tissue from wild-type and mutant mouse embryos, strongly highlighted gene networks linked to neurite development. We followed up our genomics data with functional experiments, showing that Foxp2 impacts on neurite outgrowth in primary neurons and in neuronal cell models. Our data indicate that Foxp2 modulates neuronal network formation, by directly and indirectly regulating mRNAs involved in the development and plasticity of neuronal connections.

FLOWERING LOCUS C (FLC) regulates development pathways throughout the life cycle of arabidopsis.

FLOWERING LOCUS C (FLC) has a key role in the timing of the initiation of flowering in Arabidopsis. FLC binds and represses two genes that promote flowering, FT and SOC1. We show that FLC binds to many other genes, indicating that it has regulatory roles other than the repression of flowering. We identified 505 FLC binding sites, mostly located in the promoter regions of genes and containing at least one CArG box, the motif known to be associated with MADS-box proteins such as FLC. We examined 40 of the target genes, and 20 showed increased transcript levels in an flc mutant compared with the wild type. Five genes showed decreased expression in the mutant, indicating that FLC binding can result in either transcriptional repression or activation. The genes we identified as FLC targets are involved in developmental pathways throughout the life history of the plant, many of which are associated with reproductive development. FLC is also involved in vegetative development, as evidenced by its binding to SPL15, delaying the progression from juvenile to adult phase. Some of the FLC target genes are also bound by two other MADS-box proteins, AP1 and SEP3, suggesting that MADS-box genes may operate in a network of control at different stages of the life cycle, many ultimately contributing to the development of the reproductive phase of the plant.

A genome-wide survey of imprinted genes in rice seeds reveals imprinting primarily occurs in the endosperm.

Genomic imprinting causes the expression of an allele depending on its parental origin. In plants, most imprinted genes have been identified in Arabidopsis endosperm, a transient structure consumed by the embryo during seed formation. We identified imprinted genes in rice seed where both the endosperm and embryo are present at seed maturity. RNA was extracted from embryos and endosperm of seeds obtained from reciprocal crosses between two subspecies Nipponbare (Japonica rice) and 93-11 (Indica rice). Sequenced reads from cDNA libraries were aligned to their respective parental genomes using single-nucleotide polymorphisms (SNPs). Reads across SNPs enabled derivation of parental expression bias ratios. A continuum of parental expression bias states was observed. Statistical analyses indicated 262 candidate imprinted loci in the endosperm and three in the embryo (168 genic and 97 non-genic). Fifty-six of the 67 loci investigated were confirmed to be imprinted in the seed. Imprinted loci are not clustered in the rice genome as found in mammals. All of these imprinted loci were expressed in the endosperm, and one of these was also imprinted in the embryo, confirming that in both rice and Arabidopsis imprinted expression is primarily confined to the endosperm. Some rice imprinted genes were also expressed in vegetative tissues, indicating that they have additional roles in plant growth. Comparison of candidate imprinted genes found in rice with imprinted candidate loci obtained from genome-wide surveys of imprinted genes in Arabidopsis to date shows a low degree of conservation, suggesting that imprinting has evolved independently in eudicots and monocots.

Efficient experimental design and analysis strategies for the detection of differential expression using RNA-Sequencing

BackgroundRNA sequencing (RNA-Seq) has emerged as a powerful approach for the detection of differential gene expression with both high-throughput and high resolution capabilities possible depending upon the experimental design chosen. Multiplex experimental designs are now readily available, these can be utilised to increase the numbers of samples or replicates profiled at the cost of decreased sequencing depth generated per sample. These strategies impact on the power of the approach to accurately identify differential expression. This study presents a detailed analysis of the power to detect differential expression in a range of scenarios including simulated null and differential expression distributions with varying numbers of biological or technical replicates, sequencing depths and analysis methods.ResultsDifferential and non-differential expression datasets were simulated using a combination of negative binomial and exponential distributions derived from real RNA-Seq data. These datasets were used to evaluate the performance of three commonly used differential expression analysis algorithms and to quantify the changes in power with respect to true and false positive rates when simulating variations in sequencing depth, biological replication and multiplex experimental design choices.ConclusionsThis work quantitatively explores comparisons between contemporary analysis tools and experimental design choices for the detection of differential expression using RNA-Seq. We found that the DESeq algorithm performs more conservatively than edgeR and NBPSeq. With regard to testing of various experimental designs, this work strongly suggests that greater power is gained through the use of biological replicates relative to library (technical) replicates and sequencing depth. Strikingly, sequencing depth could be reduced as low as 15% without substantial impacts on false positive or true positive rates.

Aberrant succination of proteins in fumarate hydratase-deficient mice and HLRCC patients is a robust biomarker of mutation status

Germline mutations in the FH gene encoding the Krebs cycle enzyme fumarate hydratase predispose to hereditary leiomyomatosis and renal cell cancer (HLRCC) syndrome. FH‐deficient cells and tissues accumulate high levels of fumarate, which may act as an oncometabolite and contribute to tumourigenesis. A recently proposed role for fumarate in the covalent modification of cysteine residues to S‐(2‐succinyl) cysteine (2SC) (termed protein succination) prompted us to assess 2SC levels in our existing models of HLRCC. Herein, using a previously characterized antibody against 2SC, we show that genetic ablation of FH causes high levels of protein succination. We next hypothesized that immunohistochemistry for 2SC would serve as a metabolic biomarker for the in situ detection of FH‐deficient tissues. Robust detection of 2SC was observed in Fh1 (murine FH)‐deficient renal cysts and in a retrospective series of HLRCC tumours (n = 16) with established FH mutations. Importantly, 2SC was undetectable in normal tissues (n = 200) and tumour types not associated with HLRCC (n = 1342). In a prospective evaluation of cases referred for genetic testing for HLRCC, the presence of 2SC‐modified proteins (2SCP) correctly predicted genetic alterations in FH in every case. In two series of unselected type II papillary renal cancer (PRCC), prospectively analysed by 2SCP staining followed by genetic analysis, the biomarker accurately identified previously unsuspected FH mutations (2/33 and 1/36). The investigation of whether metabolites in other tumour types produce protein modification signature(s) that can be assayed using similar strategies will be of interest in future studies of cancer. Copyright

Trans Chromosomal Methylation in Arabidopsis hybrids

The heterotic hybrid offspring of Arabidopsis accessions C24 and Landsberg erecta have altered methylomes. Changes occur most frequently at loci where parental methylation levels are different. There are context-specific biases in the nonadditive methylation patterns with mCG generally increased and mCHH decreased relative to the parents. These changes are a result of two main mechanisms, Trans Chromosomal Methylation and Trans Chromosomal deMethylation, where the methylation level of one parental allele alters to resemble that of the other parent. Regions of altered methylation are enriched around genic regions and are often correlated with changes in siRNA levels. We identified examples of genes with altered expression likely to be due to methylation changes and suggest that in crosses between the C24 and Ler accessions, epigenetic controls can be important in the generation of altered transcription levels that may contribute to the increased biomass of the hybrids.

MicroRNA-125a is over-expressed in insulin target tissues in a spontaneous rat model of Type 2 Diabetes

BackgroundMicroRNAs (miRNAs) are non-coding RNA molecules involved in post-transcriptional control of gene expression of a wide number of genes, including those involved in glucose homeostasis. Type 2 diabetes (T2D) is characterized by hyperglycaemia and defects in insulin secretion and action at target tissues. We sought to establish differences in global miRNA expression in two insulin-target tissues from inbred rats of spontaneously diabetic and normoglycaemic strains.MethodsWe used a miRNA microarray platform to measure global miRNA expression in two insulin-target tissues: liver and adipose tissue from inbred rats of spontaneously diabetic (Goto-Kakizaki [GK]) and normoglycaemic (Brown-Norway [BN]) strains which are extensively used in genetic studies of T2D. MiRNA data were integrated with gene expression data from the same rats to investigate how differentially expressed miRNAs affect the expression of predicted target gene transcripts.ResultsThe expression of 170 miRNAs was measured in liver and adipose tissue of GK and BN rats. Based on a p-value for differential expression between GK and BN, the most significant change in expression was observed for miR-125a in liver (FC = 5.61, P = 0.001, Padjusted= 0.10); this overexpression was validated using quantitative RT-PCR (FC = 13.15, P = 0.0005). MiR-125a also showed over-expression in the GK vs. BN analysis within adipose tissue (FC = 1.97, P = 0.078, Padjusted= 0.99), as did the previously reported miR-29a (FC = 1.51, P = 0.05, Padjusted= 0.99). In-silico tools assessing the biological role of predicted miR-125a target genes suggest an over-representation of genes involved in the MAPK signaling pathway. Gene expression analysis identified 1308 genes with significantly different expression between GK and BN rats (Padjusted < 0.05): 233 in liver and 1075 in adipose tissue. Pathways related to glucose and lipid metabolism were significantly over-represented among these genes. Enrichment analysis suggested that differentially expressed genes in GK compared to BN included more predicted miR-125a target genes than would be expected by chance in adipose tissue (FDR = 0.006 for up-regulated genes; FDR = 0.036 for down-regulated genes) but not in liver (FDR = 0.074 for up-regulated genes; FDR = 0.248 for down-regulated genes).ConclusionMiR-125a is over-expressed in liver in hyperglycaemic GK rats relative to normoglycaemic BN rats, and our array data also suggest miR-125a is over-expressed in adipose tissue. We demonstrate the use of in-silico tools to provide the basis for further investigation of the potential role of miR-125a in T2D. In particular, the enrichment of predicted miR-125a target genes among differentially expressed genes has identified likely target genes and indicates that integrating global miRNA and mRNA expression data may give further insights into miRNA-mediated regulation of gene expression.