Colin A. Semple

Genome-wide analysis of mammalian promoter architecture and evolution

Mammalian promoters can be separated into two classes, conserved TATA box–enriched promoters, which initiate at a well-defined site, and more plastic, broad and evolvable CpG-rich promoters. We have sequenced tags corresponding to several hundred thousand transcription start sites (TSSs) in the mouse and human genomes, allowing precise analysis of the sequence architecture and evolution of distinct promoter classes. Different tissues and families of genes differentially use distinct types of promoters. Our tagging methods allow quantitative analysis of promoter usage in different tissues and show that differentially regulated alternative TSSs are a common feature in protein-coding genes and commonly generate alternative N termini. Among the TSSs, we identified new start sites associated with the majority of exons and with 3′ UTRs. These data permit genome-scale identification of tissue-specific promoters and analysis of the cis-acting elements associated with them.

Genome-wide association scan identifies a colorectal cancer susceptibility locus on 11q23 and replicates risk loci at 8q24 and 18q21.

In a genome-wide association study to identify loci associated with colorectal cancer (CRC) risk, we genotyped 555,510 SNPs in 1,012 early-onset Scottish CRC cases and 1,012 controls (phase 1). In phase 2, we genotyped the 15,008 highest-ranked SNPs in 2,057 Scottish cases and 2,111 controls. We then genotyped the five highest-ranked SNPs from the joint phase 1 and 2 analysis in 14,500 cases and 13,294 controls from seven populations, and identified a previously unreported association, rs3802842 on 11q23 (OR = 1.1; P = 5.8 × 10−10), showing population differences in risk. We also replicated and fine-mapped associations at 8q24 (rs7014346; OR = 1.19; P = 8.6 × 10−26) and 18q21 (rs4939827; OR = 1.2; P = 7.8 × 10−28). Risk was greater for rectal than for colon cancer for rs3802842 (P < 0.008) and rs4939827 (P < 0.009). Carrying all six possible risk alleles yielded OR = 2.6 (95% CI = 1.75–3.89) for CRC. These findings extend our understanding of the role of common genetic variation in CRC etiology.

Mutations in genes encoding ribonuclease H2 subunits cause Aicardi-Goutières syndrome and mimic congenital viral brain infection

Aicardi-Goutières syndrome (AGS) is an autosomal recessive neurological disorder, the clinical and immunological features of which parallel those of congenital viral infection. Here we define the composition of the human ribonuclease H2 enzyme complex and show that AGS can result from mutations in the genes encoding any one of its three subunits. Our findings demonstrate a role for ribonuclease H in human neurological disease and suggest an unanticipated relationship between ribonuclease H2 and the antiviral immune response that warrants further investigation.

Sushi gets serious: the draft genome sequence of the pufferfish Fugu rubripes

The publication of the Fugu rubripes draft genome sequence will take this fish from culinary delicacy to potent tool in deciphering the mysteries of human genome function.

Proteins containing the UBA domain are able to bind to multi-ubiquitin chains

The UBA domain is a motif found in a variety of proteins, some of which are associated with the ubiquitin–proteasome system. We describe the isolation of a fission-yeast gene, mud1+, which encodes a UBA domain containing protein that is able to bind multi-ubiquitin chains. We show that the UBA domain is responsible for this activity. Two other proteins containing this motif, the fission-yeast homologues of Rad23 and Dsk2, are also shown to bind multi-ubiquitin chains via their UBA domains. These two proteins are implicated, along with the fission-yeast Pus1(S5a/Rpn10) subunit of the 26 S proteasome, in the recognition and turnover of substrates by this proteolytic complex.

POCUS: mining genomic sequence annotation to predict disease genes

Here we present POCUS (prioritization of candidate genes using statistics), a novel computational approach to prioritize candidate disease genes that is based on over-representation of functional annotation between loci for the same disease. We show that POCUS can provide high (up to 81-fold) enrichment of real disease genes in the candidate-gene shortlists it produces compared with the original large sets of positional candidates. In contrast to existing methods, POCUS can also suggest counterintuitive candidates.

Common variation near CDKN1A , POLD3 and SHROOM2 influences colorectal cancer risk

We performed a meta-analysis of five genome-wide association studies to identify common variants influencing colorectal cancer (CRC) risk comprising 8,682 cases and 9,649 controls. Replication analysis was performed in case-control sets totaling 21,096 cases and 19,555 controls. We identified three new CRC risk loci at 6p21 (rs1321311, near CDKN1A; P = 1.14 × 10−10), 11q13.4 (rs3824999, intronic to POLD3; P = 3.65 × 10−10) and Xp22.2 (rs5934683, near SHROOM2; P = 7.30 × 10−10) This brings the number of independent loci associated with CRC risk to 20 and provides further insight into the genetic architecture of inherited susceptibility to CRC.

Conservation and divergence in Toll-like receptor 4-regulated gene expression in primary human versus mouse macrophages

Evolutionary change in gene expression is generally considered to be a major driver of phenotypic differences between species. We investigated innate immune diversification by analyzing interspecies differences in the transcriptional responses of primary human and mouse macrophages to the Toll-like receptor (TLR)–4 agonist lipopolysaccharide (LPS). By using a custom platform permitting cross-species interrogation coupled with deep sequencing of mRNA 5′ ends, we identified extensive divergence in LPS-regulated orthologous gene expression between humans and mice (24% of orthologues were identified as “divergently regulated”). We further demonstrate concordant regulation of human-specific LPS target genes in primary pig macrophages. Divergently regulated orthologues were enriched for genes encoding cellular “inputs” such as cell surface receptors (e.g., TLR6, IL-7Rα) and functional “outputs” such as inflammatory cytokines/chemokines (e.g., CCL20, CXCL13). Conversely, intracellular signaling components linking inputs to outputs were typically concordantly regulated. Functional consequences of divergent gene regulation were confirmed by showing LPS pretreatment boosts subsequent TLR6 responses in mouse but not human macrophages, in keeping with mouse-specific TLR6 induction. Divergently regulated genes were associated with a large dynamic range of gene expression, and specific promoter architectural features (TATA box enrichment, CpG island depletion). Surprisingly, regulatory divergence was also associated with enhanced interspecies promoter conservation. Thus, the genes controlled by complex, highly conserved promoters that facilitate dynamic regulation are also the most susceptible to evolutionary change.

Duplication and selection in the evolution of primate β-defensin genes

BackgroundInnate immunity is the first line of defense against microorganisms in vertebrates and acts by providing an initial barrier to microorganisms and triggering adaptive immune responses. Peptides such as β-defensins are an important component of this defense, providing a broad spectrum of antimicrobial activity against bacteria, fungi, mycobacteria and several enveloped viruses. β-defensins are small cationic peptides that vary in their expression patterns and spectrum of pathogen specificity. Disruptions in β-defensin function have been implicated in human diseases, including cystic fibrosis, and a fuller understanding of the variety, function and evolution of human β-defensins might form the basis for novel therapies. Here we use a combination of laboratory and computational techniques to characterize the main human β-defensin locus on chromosome 8p22-p23.ResultsIn addition to known genes in the region we report the genomic structures and expression patterns of four novel human β-defensin genes and a related pseudogene. These genes show an unusual pattern of evolution, with rapid divergence between second exon sequences that encode the mature β-defensin peptides matched by relative stasis in first exons that encode signal peptides.ConclusionsWe conclude that the 8p22-p23 locus has evolved by successive rounds of duplication followed by substantial divergence involving positive selection, to produce a diverse cluster of paralogous genes established before the human-baboon divergence more than 23 million years ago. Positive selection, disproportionately favoring alterations in the charge of amino-acid residues, is implicated as driving second exon divergence in these genes.

Two novel proteins recruited by synaptonemal complex protein 1 (SYCP1) are at the centre of meiosis

Completion of meiosis in mammals depends on the formation of the synaptonemal complex, a tripartite structure that physically links homologous chromosomes during prophase I. Several components of the synaptonemal complex are known, including constituents of the cohesin core, the axial/lateral element and the transverse filaments. No protein has previously been identified as an exclusive component of the central element. Mutations in some synaptonemal-complex proteins results in impaired meiosis. In humans, cases of male infertility have been associated with failure to build the synaptonemal complex. To search for new components of the meiotic machinery, we have used data from microarray expression profiling and found two proteins localising solely to the central element of the mammalian synaptonemal complex. These new proteins, SYCE1 and CESC1, interact with the transverse filament protein SYCP1, and their localisation to the central element appears to depend on recruitment by SYCP1. This suggests a role for SYCE1 and CESC1 in synaptonemal-complex assembly, and perhaps also stability and recombination.