Summary 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 mRNA and miRNA from lymphoblastoid cell lines of 462 individuals from the 1000 Genomes Project – the first uniformly processed RNA-seq data from multiple human populations with high-quality genome sequences. We discovered extremely widespread genetic variation affecting regulation of the majority of genes, with transcript structure and expression level variation being equally common but genetically largely independent. Our characterization of causal regulatory variation sheds light on cellular mechanisms of regulatory and loss-of-function variation, and allowed 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.

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.

The availability of complete genome sequence data from both bacteria and eukaryotes provides information about the contribution of bacterial genes to the origin and evolution of mitochondria. Phylogenetic analyses based on genes located in the mitochondrial genome indicate that these genes originated from within the a-proteobacteria. A number of ancestral bacterial genes have also been transferred from the mitochondrial to the nuclear genome, as evidenced by the presence of orthologous genes in the mitochondrial genome in some species and in the nuclear genome of other species. However, a multitude of mitochondrial proteins encoded in the nucleus display no homology to bacterial proteins, indicating that these originated within the eukaryotic cell subsequent to the acquisition of the endosymbiont. An analysis of the expression patterns of yeast nuclear genes coding for mitochondrial proteins has shown that genes predicted to be of eukaryotic origin are mainly translated on polysomes that are free in the cytosol whereas those of putative bacterial origin are translated on polysomes attached to the mitochondrion. The strong relationship with a-proteobacterial genes observed for some mitochondrial genes, combined with the lack of such a relationship for others, indicates that the modern mitochondrial proteome is the product of both reductive and expansive processes.

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.

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 variablenumber 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.

Cell-free DNA analysis is becoming adopted for first line aneuploidy screening, however for most healthcare programs, cost and workflow complexity is limiting adoption of the test. We report a novel cost effective method, the Vanadis NIPT assay, designed for high precision digitally-enabled measurement of chromosomal aneuploidies in maternal plasma. Reducing NIPT assay complexity is achieved by using novel molecular probe technology that specifically label target chromosomes combined with a new readout format using a nanofilter to enrich single molecules for imaging and counting without DNA amplification, microarrays or sequencing. The primary objective of this study was to assess the Vanadis NIPT assay with respect to analytical precision and clinical feasibility. Analysis of reference DNA samples indicate that samples which are challenging to analyze with low fetal-fraction can be readily detected with a limit of detection determined at <2% fetal-fraction. In total of 286 clinical samples were analysed and 30 out of 30 pregnancies affected by trisomy 21 were classified correctly. This method has the potential to make cost effective NIPT more widely available with more women benefiting from superior detection and false positive rates.

Bartonella henselae is an emerging zoonotic pathogen causing a wide range of disease manifestations in humans. In this study, we report on the analysis of the sarcosine-insoluble outer membrane fraction of B. henselae ATCC 49882 Houston-1 by one-dimensional sodium dodecyl sulfate-polyacrylamide gel electrophoresis (1-D SDS-PAGE) and two-dimensional nonequilibrium pH gradient polyacrylamide gel electrophoresis (2-D NEPHGE). Protein species were identified by matrix-assisted laser desorption/ionization-time of flight-mass spectrometry (MALDI-TOF-MS) and subsequent database query against the B. henselae genome sequence. Subcellular fractionation, application of the ionic detergent lauryl sarcosine, assessment of trypsin sensitivity, and heat modifiability of surface-exposed proteins represented valuable tools for the analysis of the outer membrane subproteome of B. henselae. 2-D NEPHGE was applied to display and catalogue a substantial number of proteins associated with the B. henselae sarcosine-insoluble outer membrane fraction, resulting in the establishment of a first 2-D reference map of this compartment. Thus, 53 distinct protein species associated with the outer membrane subproteome fraction were identified. This study provides novel insights into the membrane biology and the associated putative virulence factors of this pathogen of increasing medical importance.

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.