CBRG, Oxford University

The Computational Biology Research Group (CBRG) provides computing support for bioinformatics analysis at the University of Oxford. We have expertise in many aspects of bioinformatics (sequence analysis, microarrays, proteomics and integration). We especially encourage collaborations that require writing custom software, bioinformatics tools and databases.

An account with the CBRG has many benefits and gives automatic access to a large number of molecular biology computing packages and to numerous biological databases.

Bioinformatics analysis tools online:

A wide range of bioinformatics programs are available online via EMBOSS Explorer.
You will need a molbiol username and password to use these tools.

Other web tools available include: BLAST, MASCOT and BASE (see full list).

We have combined the circular chromosome conformation capture protocol with high-throughput, genome-wide sequence analysis to characterize the cis-acting regulatory network at a single locus. In contrast to methods which identify large interacting regions (10-1000 kb), the 4C approach provides a comprehensive, high-resolution analysis of a specific locus with the aim of defining, in detail, the cis-regulatory elements controlling a single gene or gene cluster. Using the human Î±-globin locus as a model, we detected all known local and long-range interactions with this gene cluster. In addition, we identified two interactions with genes located 300 kb (NME4) and 625 kb (FAM173a) from the Î±-globin cluster.

Congenital myasthenic syndromes are a heterogeneous group of inherited disorders that arise from impaired signal transmission at the neuromuscular synapse. They are characterized by fatigable muscle weakness. We performed linkage analysis, whole-exome and whole-genome sequencing to determine the underlying defect in patients with an inherited limb-girdle pattern of myasthenic weakness. We identify ALG14 and ALG2 as novel genes in which mutations cause a congenital myasthenic syndrome. Through analogy with yeast, ALG14 is thought to form a multiglycosyltransferase complex with ALG13 and DPAGT1 that catalyses the first two committed steps of asparagine-linked protein glycosylation. We show that ALG14 is concentrated at the muscle motor endplates and small interfering RNA silencing of ALG14 results in reduced cell-surface expression of muscle acetylcholine receptor expressed in human embryonic kidney 293 cells. ALG2 is an alpha-1,3-mannosyltransferase that also catalyses early steps in the asparagine-linked glycosylation pathway. Mutations were identified in two kinships, with mutation ALG2p.Val68Gly found to severely reduce ALG2 expression both in patient muscle, and in cell cultures. Identification of DPAGT1, ALG14 and ALG2 mutations as a cause of congenital myasthenic syndrome underscores the importance of asparagine-linked protein glycosylation for proper functioning of the neuromuscular junction. These syndromes form part of the wider spectrum of congenital disorders of glycosylation caused by impaired asparagine-linked glycosylation. It is likely that further genes encoding components of this pathway will be associated with congenital myasthenic syndromes or impaired neuromuscular transmission as part of a more severe multisystem disorder. Our findings suggest that treatment with cholinesterase inhibitors may improve muscle function in many of the congenital disorders of glycosylation.

The SNP rs12252-C allele alters the function of interferon-induced transmembrane protein-3 increasing the disease severity of influenza virus infection in Caucasians, but the allele is rare. However, rs12252-C is much more common in Han Chinese. Here we report that the CC genotype is found in 69% of Chinese patients with severe pandemic influenza A H1N1/09 virus infection compared with 25% in those with mild infection. Specifically, the CC genotype was estimated to confer a sixfold greater risk for severe infection than the CT and TT genotypes. More importantly, because the risk genotype occurs with such a high frequency, its effect translates to a large population-attributable risk of 54.3% for severe infection in the Chinese population studied compared with 5.4% in Northern Europeans. Interferon-induced transmembrane protein-3 genetic variants could, therefore, have a strong effect of the epidemiology of influenza in China and in people of Chinese descent.

The extracellular signal-related kinases 1 and 2 (ERK1/2) are key proteins mediating mitogen-activated protein kinase signaling downstream of RAS: phosphorylation of ERK1/2 leads to nuclear uptake and modulation of multiple targets. Here, we show that reduced dosage of ERF, which encodes an inhibitory ETS transcription factor directly bound by ERK1/2 (refs. 2,3,4,5,6,7), causes complex craniosynostosis (premature fusion of the cranial sutures) in humans and mice. Features of this newly recognized clinical disorder include multiple-suture synostosis, craniofacial dysmorphism, Chiari malformation and language delay. Mice with functional Erf levels reduced to âˆ¼30% of normal exhibit postnatal multiple-suture synostosis; by contrast, embryonic calvarial development appears mildly delayed. Using chromatin immunoprecipitation in mouse embryonic fibroblasts and high-throughput sequencing, we find that ERF binds preferentially to elements away from promoters that contain RUNX or AP-1 motifs. This work identifies ERF as a novel regulator of osteogenic stimulation by RAS-ERK signaling, potentially by competing with activating ETS factors in multifactor transcriptional complexes.

Craniosynostosis, the premature fusion of the cranial sutures, is a heterogeneous disorder with a prevalence of âˆ¼1 in 2,200 (refs. 1,2). A specific genetic etiology can be identified in âˆ¼21% of cases, including mutations of TWIST1, which encodes a class II basic helix-loop-helix (bHLH) transcription factor, and causes Saethre-Chotzen syndrome, typically associated with coronal synostosis. Using exome sequencing, we identified 38 heterozygous TCF12 mutations in 347 samples from unrelated individuals with craniosynostosis. The mutations predominantly occurred in individuals with coronal synostosis and accounted for 32% and 10% of subjects with bilateral and unilateral pathology, respectively. TCF12 encodes one of three class I E proteins that heterodimerize with class II bHLH proteins such as TWIST1. We show that TCF12 and TWIST1 act synergistically in a transactivation assay and that mice doubly heterozygous for loss-of-function mutations in Tcf12 and Twist1 have severe coronal synostosis. Hence, the dosage of TCF12-TWIST1 heterodimers is critical for normal coronal suture development.

Craniofrontonasal syndrome (CFNS), an X-linked disorder caused by loss-of-function mutations of EFNB1, exhibits a paradoxical sex reversal in phenotypic severity: females characteristically have frontonasal dysplasia, craniosynostosis and additional minor malformations, but males are usually more mildly affected with hypertelorism as the only feature. X-inactivation is proposed to explain the more severe outcome in heterozygous females, as this leads to functional mosaicism for cells with differing expression of EPHRIN-B1, generating abnormal tissue boundaries-a process that cannot occur in hemizygous males. Apparently challenging this model, males occasionally present with a more severe female-like CFNS phenotype. We hypothesized that such individuals might be mosaic for EFNB1 mutations and investigated this possibility in multiple tissue samples from six sporadically presenting males. Using denaturing high performance liquid chromatography, massively parallel sequencing and multiplex-ligation-dependent probe amplification (MLPA) to increase sensitivity above standard dideoxy sequencing, we identified mosaic mutations of EFNB1 in all cases, comprising three missense changes, two gene deletions and a novel point mutation within the 5' untranslated region (UTR). Quantification by Pyrosequencing and MLPA demonstrated levels of mutant cells between 15 and 69%. The 5' UTR variant mutates the stop codon of a small upstream open reading frame that, using a dual-luciferase reporter construct, was demonstrated to exacerbate interference with translation of the wild-type protein. These results demonstrate a more severe outcome in mosaic than in constitutionally deficient males in an X-linked dominant disorder and provide further support for the cellular interference mechanism, normally related to X-inactivation in females.

Human U1 small nuclear (sn)RNA, required for splicing of pre-mRNA, is encoded by genes on chromosome 1 (1p36). Imperfect copies of these U1 snRNA genes, also located on chromosome 1 (1q12-21), were thought to be pseudogenes. However, many of these "variant" (v)U1 snRNA genes produce fully processed transcripts. Using antisense oligonucleotides to block the activity of a specific vU1 snRNA in HeLa cells, we have identified global transcriptome changes following interrogation of the Affymetrix Human Exon ST 1.0 array. Our results indicate that this vU1 snRNA regulates expression of a subset of target genes at the level of pre-mRNA processing. This is the first indication that variant U1 snRNAs have a biological function in vivo. Furthermore, some vU1 snRNAs are packaged into unique ribonucleoproteins (RNPs), and many vU1 snRNA genes are differentially expressed in human embryonic stem cells (hESCs) and HeLa cells, suggesting developmental control of RNA processing through expression of different sets of vU1 snRNPs.

Natural killer (NK) cells play an important role in the containment of HIV replication during primary infection, though their functions are impaired during chronic HIV infection. Infants experience more rapid HIV disease progression than adults, but contributions of infant NK cells to containing HIV infection are unknown. The aim of this study was to determine the impact of HIV infection on infant NK cell phenotype by evaluating samples and data from a cohort study of women and their infants, conducted in Nairobi, Kenya between 1999 and 2003. The percentage and phenotype of NK cells was evaluated longitudinally by multi-parameter flow cytometry over the first year of life in HIV-infected (HIV+, = 16), HIV-exposed uninfected (HIV-EU, n = 6), and healthy unexposed controls (HIV-, n = 4). At birth, NK subset distributions based on expression of CD56 and CD16 did not differ between HIV+, HIV-EU, or HIV- infants. However, HIV infection was associated with a subsequent decline in NK cells as a percentage of total lymphocytes (p < 0.001), and an expanding proportion of CD56-CD16+ NK cells (p < 0.001). Activated CD38(bright)CD69+ NK cells were more frequent in the HIV+ infants, followed by HIV-EU and HIV- infants, in both CD56(dim) (p = 0.005) and CD56(bright) compartments (p = 0.03). HIV infection and exposure was also associated with a significant decline in the percentage of perforin-expressing NK cells in the CD56(dim) compartment over the first year of life, with HIV+ infants losing approximately 2.5% (p < 0.001) and HIV-EU infants losing 3.0% (p = 0.01) of perforin+ cells per month. Thus, infant HIV infection is associated with alterations in NK cell subsets, activation, and cytolytic potential that could contribute to their poor control over HIV infection. Furthermore, exposure to HIV infection in infants who escaped infection is also associated with alterations in NK cells that may contribute to the reduced ability to fight infections that is observed in HIV-EU infants.

Tumour necrosis factor alpha (TNFÎ±) is a potent cytokine that signals through nuclear factor kappa B (NFÎºB) to activate a subset of human genes. It is usually assumed that this involves RNA polymerases transcribing responsive genes wherever they might be in the nucleus. Using primary human endothelial cells, variants of chromosome conformation capture (including 4C and chromatin interaction analysis with paired-end tag sequencing), and fluorescence in situ hybridization to detect single nascent transcripts, we show that TNFÎ± induces responsive genes to congregate in discrete 'NFÎºB factories'. Some factories further specialize in transcribing responsive genes encoding micro-RNAs that target downregulated mRNAs. We expect all signalling pathways to contain this extra leg, where responding genes are transcribed in analogous specialized factories.

Carpenter syndrome is an autosomal-recessive multiple-congenital-malformation disorder characterized by multisuture craniosynostosis and polysyndactyly of the hands and feet; many other clinical features occur, and the most frequent include obesity, umbilical hernia, cryptorchidism, and congenital heart disease. Mutations of RAB23, encoding a small GTPase that regulates vesicular transport, are present in the majority of cases. Here, we describe a disorder caused by mutations in multiple epidermal-growth-factor-like-domains 8 (MEGF8), which exhibits substantial clinical overlap with Carpenter syndrome but is frequently associated with abnormal left-right patterning. We describe five affected individuals with similar dysmorphic facies, and three of them had either complete situs inversus, dextrocardia, or transposition of the great arteries; similar cardiac abnormalities were previously identified in a mouse mutant for the orthologous Megf8. The mutant alleles comprise one nonsense, three missense, and two splice-site mutations; we demonstrate in zebrafish that, in contrast to the wild-type protein, the proteins containing all three missense alterations provide only weak rescue of an early gastrulation phenotype induced by Megf8 knockdown. We conclude that mutations in MEGF8 cause a Carpenter syndrome subtype frequently associated with defective left-right patterning, probably through perturbation of signaling by hedgehog and nodal family members. We did not observe any subject with biallelic loss-of function mutations, suggesting that some residual MEGF8 function might be necessary for survival and might influence the phenotypes observed.

The dominant congenital disorders Apert syndrome, achondroplasia and multiple endocrine neoplasia-caused by specific missense mutations in the FGFR2, FGFR3 and RET proteins respectively-represent classical examples of paternal age-effect mutation, a class that arises at particularly high frequencies in the sperm of older men. Previous analyses of DNA from randomly selected cadaveric testes showed that the levels of the corresponding FGFR2, FGFR3 and RET mutations exhibit very uneven spatial distributions, with localised hotspots surrounded by large mutation-negative areas. These studies imply that normal testes are mosaic for clusters of mutant cells: these clusters are predicted to have altered growth and signalling properties leading to their clonal expansion (selfish spermatogonial selection), but DNA extraction eliminates the possibility to study such processes at a tissue level. Using a panel of antibodies optimised for the detection of spermatocytic seminoma, a rare tumour of spermatogonial origin, we demonstrate that putative clonal events are frequent within normal testes of elderly men (mean age: 73.3 yrs) and can be classed into two broad categories. We found numerous small (less than 200 cells) cellular aggregations with distinct immunohistochemical characteristics, localised to a portion of the seminiferous tubule, which are of uncertain significance. However more infrequently we identified additional regions where entire seminiferous tubules had a circumferentially altered immunohistochemical appearance that extended through multiple serial sections that were physically contiguous (up to 1 mm in length), and exhibited enhanced staining for antibodies both to FGFR3 and a marker of downstream signal activation, pAKT. These findings support the concept that populations of spermatogonia in individual seminiferous tubules in the testes of older men are clonal mosaics with regard to their signalling properties and activation, thus fulfilling one of the specific predictions of selfish spermatogonial selection.

X chromosome inactivation involves multiple levels of chromatin modification, established progressively and in a stepwise manner during early development. The chromosomal protein Smchd1 was recently shown to play an important role in DNA methylation of CpG islands (CGIs), a late step in the X inactivation pathway that is required for long-term maintenance of gene silencing. Here we show that inactive X chromosome (Xi) CGI methylation can occur via either Smchd1-dependent or -independent pathways. Smchd1-dependent CGI methylation, the primary pathway, is acquired gradually over an extended period, whereas Smchd1-independent CGI methylation occurs rapidly after the onset of X inactivation. The de novo methyltransferase Dnmt3b is required for methylation of both classes of CGI, whereas Dnmt3a and Dnmt3L are dispensable. Xi CGIs methylated by these distinct pathways differ with respect to their sequence characteristics and immediate chromosomal environment. We discuss the implications of these results for understanding CGI methylation during development.

Congenital myasthenic syndromes are a heterogeneous group of inherited disorders that arise from impaired signal transmission at the neuromuscular synapse. They are characterized by fatigable muscle weakness. We performed whole-exome sequencing to determine the underlying defect in a group of individuals with an inherited limb-girdle pattern of myasthenic weakness. We identify DPAGT1 as a gene in which mutations cause a congenital myasthenic syndrome. We describe seven different mutations found in five individuals with DPAGT1 mutations. The affected individuals share a number of common clinical features, including involvement of proximal limb muscles, response to treatment with cholinesterase inhibitors and 3,4-diaminopyridine, and the presence of tubular aggregates in muscle biopsies. Analyses of motor endplates from two of the individuals demonstrate a severe reduction of endplate acetylcholine receptors. DPAGT1 is an essential enzyme catalyzing the first committed step of N-linked protein glycosylation. Our findings underscore the importance of N-linked protein glycosylation for proper functioning of the neuromuscular junction. Using the DPAGT1-specific inhibitor tunicamycin, we show that DPAGT1 is required for efficient glycosylation of acetylcholine-receptor subunits and for efficient export of acetylcholine receptors to the cell surface. We suggest that the primary pathogenic mechanism of DPAGT1 mutations is reduced levels of acetylcholine receptors at the endplate region. These individuals share clinical features similar to those of congenital myasthenic syndrome due to GFPT1 mutations, and their disorder might be part of a larger subgroup comprising the congenital myasthenic syndromes that result from defects in the N-linked glycosylation pathway and that manifest through impaired neuromuscular transmission.

Delayed HIV-1 disease progression is associated with a single nucleotide polymorphism upstream of the HLA-C gene that correlates with differential expression of the HLA-C Ag. This polymorphism was recently shown to be a marker for a protective variant in the 3'UTR of HLA-C that disrupts a microRNA binding site, resulting in enhanced HLA-C expression at the cell surface. Whether individuals with "high" HLA-C expression show a stronger HLA-C-restricted immune response exerting better viral control than that of their counterparts has not been established. We hypothesized that the magnitude of the HLA-C-restricted immune pressure on HIV would be greater in subjects with highly expressed HLA-C alleles. Using a cohort derived from a unique narrow source epidemic in China, we identified mutations in HIV proviral DNA exclusively associated with HLA-C, which were used as markers for the intensity of the immune pressure exerted on the virus. We found an increased frequency of mutations in individuals with highly expressed HLA-C alleles, which also correlated with IFN-Î³ production by HLA-C-restricted CD8(+) T cells. These findings show that immune pressure on HIV is stronger in subjects with the protective genotype and highlight the potential role of HLA-C-restricted responses in HIV control. This is, to our knowledge, the first in vivo evidence supporting the protective role of HLA-C-restricted responses in nonwhites during HIV infection.

Determining previous infecting dengue virus (DENV) serotypes has been difficult due to highly cross-reactive immune responses from previous DENV infections. Determining the correlates of serotype-specific immune responses would be crucial in understanding dengue transmission in the community and would also help to determine the correlates of protective immune responses. Therefore, we set out to define highly conserved, serotype-specific regions of the DENVs. Serotype-specific and highly conserved regions of the four DENV serotypes were identified using Basic Local Alignment Search Tool (BLAST) searches and custom perl scripts. Using ex-vivo and cultured enzyme-linked immunospot (ELISPOT) assays, we identified serotype-specific T cell epitopes within the four DENV serotypes in healthy adult donors from Sri Lanka. We identified T cell responses to 19 regions of the four DENV serotypes. Six peptides were from the NS2A region and four peptides were from the NS4A region. All immune donors responded to peptides of at least two DENV serotypes, suggesting that heterologous infection is common in Sri Lanka. Eight of 20 individuals responded to at least two peptides of DENV-4, despite this serotype not being implicated previously in any of the epidemics in Sri Lanka. The use of these regions to determine past and current infecting DENV serotypes will be of value to characterize further the dynamics of silent dengue transmission in the community. In addition, these T cell responses to these regions could be used to characterize DENV serotype-specific immune responses and thus possibly help us to understand the immune correlates of a protective immune response.

CBRG support
This file last modified Tuesday May 22, 2012
© Computational Biology Research Group, 2004-10