Olof Karlberg

Transcriptome and genome sequencing uncovers functional variation in humans.

Genome sequencing projects are discovering millions of genetic variants in humans, and interpretation of their functional effects is essential for understanding the genetic basis of variation in human traits. Here we report sequencing and deep analysis of messenger RNA and microRNA from lymphoblastoid cell lines of 462 individuals from the 1000 Genomes Project—the first uniformly processed high-throughput RNA-sequencing data from multiple human populations with high-quality genome sequences. We discover extremely widespread genetic variation affecting the regulation of most genes, with transcript structure and expression level variation being equally common but genetically largely independent. Our characterization of causal regulatory variation sheds light on the cellular mechanisms of regulatory and loss-of-function variation, and allows us to infer putative causal variants for dozens of disease-associated loci. Altogether, this study provides a deep understanding of the cellular mechanisms of transcriptome variation and of the landscape of functional variants in the human genome.

The Dual Origin of the Yeast Mitochondrial Proteome

We propose a scheme for the origin of mitochondria based on phylogenetic reconstructions with more than 400 yeast nuclear genes that encode mitochondrial proteins. Half of the yeast mitochondrial proteins have no discernable bacterial homologues, while one‐tenth are unequivocally of α‐proteobacterial origin. These data suggest that the majority of genes encoding yeast mitochondrial proteins are descendants of two different genomic lineages that have evolved in different modes. First, the ancestral free‐living α‐proteobacterium evolved into an endosymbiont of an anaerobic host. Most of the ancestral bacterial genes were lost, but a small fraction of genes supporting bioenergetic and translational processes were retained and eventually transferred to what became the host nuclear genome. In a second, parallel mode, a larger number of novel mitochondrial genes were recruited from the nuclear genome to complement the remaining genes from the bacterial ancestor. These eukaryotic genes, which are primarily involved in transport and regulatory functions, transformed the endosymbiont into an ATP‐exporting organelle. Copyright

Reproducibility of high-throughput mRNA and small RNA sequencing across laboratories

RNA sequencing is an increasingly popular technology for genome-wide analysis of transcript sequence and abundance. However, understanding of the sources of technical and interlaboratory variation is still limited. To address this, the GEUVADIS consortium sequenced mRNAs and small RNAs of lymphoblastoid cell lines of 465 individuals in seven sequencing centers, with a large number of replicates. The variation between laboratories appeared to be considerably smaller than the already limited biological variation. Laboratory effects were mainly seen in differences in insert size and GC content and could be adequately corrected for. In small-RNA sequencing, the microRNA (miRNA) content differed widely between samples owing to competitive sequencing of rRNA fragments. This did not affect relative quantification of miRNAs. We conclude that distributing RNA sequencing among different laboratories is feasible, given proper standardization and randomization procedures. We provide a set of quality measures and guidelines for assessing technical biases in RNA-seq data.

Deep origin of plastid/parasite ATP/ADP translocases.

Membrane proteins that transport ATP and ADP have been identified in mitochondria, plastids, and obligate intracellular parasites. The mitochondrial ATP/ADP transporters are derived from a broad-specificity transport family of eukaryotic origin, whereas the origin of the plastid/parasite ATP/ADP translocase is more elusive. Here we present the sequences of five genes coding for ATP/ADP translocases from four species of Rickettsia. The results are consistent with an early duplication and divergence of the five ATP/ADP translocases within the rickettsial lineage. A comparison of the phylogenetic depths of the mitochondrial and the plastid/parasite ATP/ADP translocases indicates a deep origin for both transporters. The results provide no evidence for a recent acquisition of the ATP/ADP transporters in Rickettsia via horizontal gene transfer, as previously suggested. A possible function of the two types of ATP/ADP translocases was to allow switches between glycolysis and aerobic respiration in the early eukaryotic cell and its endosymbiont.

Five-Member Gene Family of Bartonella quintana

ABSTRACT Bartonella quintana, the agent of trench fever and an etiologic agent of bacillary angiomatosis, has an extraordinarily high hemin requirement for growth compared to other bacterial pathogens. We previously identified the major hemin receptor of the pathogen as a 30-kDa surface protein, termed HbpA. This report describes four additional homologues that share approximately 48% amino acid sequence identity with hbpA. Three of the genes form a paralagous cluster, termed hbpCAB, whereas the other members, hbpD and hbpE, are unlinked. Secondary structure predictions and other evidence suggest that Hbp family members are β-barrels located in the outer membrane and contain eight transmembrane domains plus four extracellular loops. Homologs from a variety of gram-negative pathogens were identified, including Bartonella henselae Pap31, Brucella Omp31, Agrobacterium tumefaciens Omp25, and neisserial opacity proteins (Opa). Family members expressed in vitro-synthesized proteins ranging from ca. 26.5 to 35.1 kDa, with the exception of HbpB, an ∼55.9-kDa protein whose respective gene has been disrupted by a ∼510 GC-rich element containing variable-number tandem repeats. Transcription analysis by quantitative reverse transcriptase-PCR (RT-PCR) indicates that all family members are expressed under normal culture conditions, with hbpD and hbpB transcripts being the most abundant and the rarest, respectively. Mutagenesis of hbpA by allelic exchange produced a strain that exhibited an enhanced hemin-binding phenotype relative to the parental strain, and analysis by quantitative RT-PCR showed elevated transcript levels for the other hbp family members, suggesting that compensatory expression occurs.

Performance of Microarray and Liquid Based Capture Methods for Target Enrichment for Massively Parallel Sequencing and SNP Discovery

Targeted sequencing is a cost-efficient way to obtain answers to biological questions in many projects, but the choice of the enrichment method to use can be difficult. In this study we compared two hybridization methods for target enrichment for massively parallel sequencing and single nucleotide polymorphism (SNP) discovery, namely Nimblegen sequence capture arrays and the SureSelect liquid-based hybrid capture system. We prepared sequencing libraries from three HapMap samples using both methods, sequenced the libraries on the Illumina Genome Analyzer, mapped the sequencing reads back to the genome, and called variants in the sequences. 74–75% of the sequence reads originated from the targeted region in the SureSelect libraries and 41–67% in the Nimblegen libraries. We could sequence up to 99.9% and 99.5% of the regions targeted by capture probes from the SureSelect libraries and from the Nimblegen libraries, respectively. The Nimblegen probes covered 0.6 Mb more of the original 3.1 Mb target region than the SureSelect probes. In each sample, we called more SNPs and detected more novel SNPs from the libraries that were prepared using the Nimblegen method. Thus the Nimblegen method gave better results when judged by the number of SNPs called, but this came at the cost of more over-sampling.

Genome-wide profiling of target genes for the systemic lupus erythematosus-associated transcription factors IRF5 and STAT4

Background The transcription factors interferon regulatory factor 5 (IRF5) and signal transducer and activator of transcription 4 (STAT4) are encoded by two of the strongest susceptibility genes for systemic lupus erythematosus (SLE). Objective To investigate the target genes and functional roles of IRF5 and STAT4 in human peripheral blood mononuclear cells (PBMCs). Methods Chromatin immunoprecipitation-sequencing (ChIP-seq) was performed in PBMCs stimulated to activate IRF5 and STAT4. The expression of the target genes of IRF5 and STAT4 was investigated in a publicly available dataset generated from PBMCs from patients with SLE and healthy controls. The genomic regions bound by the transcription complexes mediated by IRF5 and STAT4 were examined for transcription factor binding motifs and SLE-associated sequence variants. Results More than 7000 target genes for IRF5 and STAT4 were identified in stimulated PBMCs. These genes were enriched to functional pathways in the type I interferon system, and have key roles in the inflammatory response. The expression patterns of the target genes were characteristic for patients with SLE. The transcription factors high mobility group-I/Y, specificity protein 1, and paired box 4 may function cooperatively with IRF5 and STAT4 in transcriptional regulation. Eight of the target regions for IRF5 and STAT4 contain SLE-associated sequence variants. Conclusions By participating in transcription complex with other co-factors, IRF5 and STAT4 harbour the potential of regulating a large number of target genes, which may contribute to their strong association with SLE.

Efficient application of next-generation sequencing for the diagnosis of rare genetic syndromes

Aims The causes of intellectual disability, which affects 1%–3% of the general population, are highly heterogeneous and the genetic defect remains unknown in around 40% of patients. The application of next-generation sequencing is changing the nature of biomedical diagnosis. This technology has quickly become the method of choice for searching for pathogenic mutations in rare uncharacterised genetic diseases. Methods Whole-exome sequencing was applied to a series of families affected with intellectual disability in order to identify variants underlying disease phenotypes. Results We present data of three families in which we identified the disease-causing mutations and which benefited from receiving a clinical diagnosis: Cornelia de Lange, Cohen syndrome and Dent-2 disease. The genetic heterogeneity and the variability in clinical presentation of these disorders could explain why these patients are difficult to diagnose. Conclusions The accessibility to next-generation sequencing allows clinicians to save much time and cost in identifying the aetiology of rare diseases. The presented cases are excellent examples that demonstrate the efficacy of next-generation sequencing in rare disease diagnosis.

A novel splicing mutation in the IQSEC2 gene that modulates the phenotype severity in a family with intellectual disability.

The IQSEC2 gene is located on chromosome Xp11.22 and encodes a guanine nucleotide exchange factor for the ADP-ribosylation factor family of small GTPases. This gene is known to have a significant role in cytoskeletal organization, dendritic spine morphology and synaptic organization. Variants in IQSEC2 cause moderate to severe intellectual disability in males and a variable phenotype in females because this gene escapes from X-chromosome inactivation. Here we report on the first splicing variant in IQSEC2 (g.88032_88033del; NG_021296.1) that co-segregates in a family diagnosed with an X-linked form of ID. In a percentage of the cells, the variant activates an intraexonic splice acceptor site that abolishes 26 amino acids from the highly conserved PH domain of IQSEC2 and creates a premature stop codon 36 amino acids later in exon 13. Interestingly, the percentage of aberrant splicing seems to correlate with the severity of the disease in each patient. The impact of this variant in the target tissue is unknown, but we can hypothesize that these differences may be related to the amount of abnormal IQSEC2 transcript. To our knowledge, we are reporting a novel mechanism of IQSEC2 involvement in ID. Variants that affect splicing are related to many genetic diseases and the understanding of their role in disease expands potential opportunities for gene therapy. Modulation of aberrant splicing transcripts can become a potent therapeutic approach for many of these diseases.