Allelic mapping bias in RNA-sequencing is not a major confounder in eQTL studies

Background

RNA sequencing (RNA-seq) is the current gold-standard method to quantify gene expression for expression quantitative trait locus (eQTL) studies. However, a potential caveat in these studies is that RNA-seq reads carrying the non-reference allele of variant loci can have lower probability to map correctly to the reference genome, which could bias gene quantifications and cause false positive eQTL associations. In this study, we analyze the effect of this allelic mapping bias in eQTL discovery.

Results

We simulate RNA-seq read mapping over 9.5 M common SNPs and indels, with 15.6% of variants showing biased mapping rate for reference versus non-reference reads. However, removing potentially biased RNA-seq reads from an eQTL dataset of 185 individuals has a very small effect on gene and exon quantifications and eQTL discovery. We detect only a handful of likely false positive eQTLs, and overall eQTL SNPs show no significant enrichment for high mapping bias.

Conclusion

Our results suggest that RNA-seq quantifications are generally robust against allelic mapping bias, and that this does not have a severe effect on eQTL discovery. Nevertheless, we provide our catalog of putatively biased loci to allow better controlling for mapping bias to obtain more accurate results in future RNA-seq studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0467-2) contains supplementary material, which is available to authorized users.     

Characterization of low-density granulocytes in COVID-19

Severe COVID-19 is characterized by extensive pulmonary complications, to which host immune responses are believed to play a role. As the major arm of innate immunity, neutrophils are one of the first cells recruited to the site of infection where their excessive activation can contribute to lung pathology. Low-density granulocytes (LDGs) are circulating neutrophils, whose numbers increase in some autoimmune diseases and cancer, but are poorly characterized in acute viral infections. Using flow cytometry, we detected a significant increase of LDGs in the blood of acute COVID-19 patients, compared to healthy controls. Based on their surface marker expression, COVID-19-related LDGs exhibit four different populations, which display distinctive stages of granulocytic development and most likely reflect emergency myelopoiesis. Moreover, COVID-19 LDGs show a link with an elevated recruitment and activation of neutrophils. Functional assays demonstrated the immunosuppressive capacities of these cells, which might contribute to impaired lymphocyte responses during acute disease. Taken together, our data confirms a significant granulocyte activation during COVID-19 and suggests that granulocytes of lower density play a role in disease progression.     

Ectomycorrhizal community structure varies among Norway spruce (Picea abies) clones

In northern boreal forests, the diversity of ectomycorrhizal (ECM) species is much greater than that of their host trees. This field study investigated the role of individual trees in shaping the ECM community. We compared ECM communities of eight Norway spruce (Picea abies) clones planted in a clear-cut area in 1994 with a randomized block design. In 2003, the ECM fungi were identified from randomly sampled root tips using denaturing gradient gel electrophoresis (DGGE) and rDNA internal transcribed spacer (ITS) sequence similarity. ECM diversity varied among clone groups, showing twofold growth differences. Moreover, according to detrended correspondence analysis (DCA), ECM community structure varied not only among but also within slow-growing or fast-growing clones. Results suggest that ECM diversity and community structure are related to the growth rate or size of the host. A direct or indirect influence of host genotype was also observed, and we therefore suggest that individual trees are partly responsible for the high diversity and patchy distribution of ECM communities in boreal forests.     

Cytovillin, a microvillar Mr 75,000 protein. cDNA sequence, prokaryotic expression, and chromosomal localization

Cytovillin is a microvillar cytoplasmic peripheral membrane protein, with prominent expression in vivo in placental syncytiotrophoblasts and certain human tumors. Cytovillin cDNA was cloned from a human placental lambda gt11 library using affinity purified antibodies. The identity of cytovillin cDNA clones was confirmed by expression of cytovillin in Escherichia coli and using antibodies raised against the expressed fusion protein in comparison with antibodies against cytovillin purified from cultured human choriocarcinoma cells. In these cells Northern blotting analysis identified a major 3.5-kilobase cytovillin mRNA. The cDNA encodes a protein of 575 amino acids corresponding to a molecular weight of 68,084. According to secondary structure prediction, cytovillin is a hydrophilic protein with an extensive internal alpha-helical region ending in a sequence of 7 consecutive prolines. The predicted alpha-helical region showed limited homology to alpha-helical regions of cytoskeletal proteins and certain other proteins, but no extensive homologies were found in the cytovillin cDNA or the deduced amino acid sequence to other registered DNA or protein sequences. Southern blot analysis of a DNA panel of human mouse somatic cell hybrids localized the cytovillin gene to the end of the long arm of chromosome 6 (6q22-q27). Our results show that cytovillin is representative of a novel class of microvillar proteins.     

The Arrival of Siberian Ancestry Connecting the Eastern Baltic to Uralic Speakers further East

Download : Download video (141MB)     

Production of pectin methylesterase from Erwinia chrysanthemi B374 in Bacillus subtilis

Summary The gene coding for pectin methylesterase (PME) of Erwinia chrysanthemi B374 (pme) was cloned by a polymerase chain reaction. The pme gene was expressed in Bacillus subtilis using a secretion vector based on the promoter and signal sequence of the -amylase gene from B. amyloliquefaciens. The cultivation of B. subtilis cells carrying the cloned pme resulted in efficient secretion of PME into the culture medium based on enzymatic and sodium dodecyl sulphate-polyacrylamide gel electrophoresis characterizations. The NH₂-terminal sequence analysis of the secreted PME revealed two different NH₂-termini. Heterologous processing was probably due to a second putative signal peptidase cleavage site at the joint region between the PME and -amylase signal peptide. Offprint requests to: R. Heikinheimo     

Collagen prolyl 4-hydroxylase tetramers and dimers show identical decreases in Km values for peptide substrates with increasing chain length: mutation of one of the two catalytic sites in the tetramer inactivates the enzyme by more than half

The collagen prolyl 4-hydroxylases (collagen P4Hs, EC 1.14.11.2) play a key role in the synthesis of the extracellular matrix. The vertebrate enzymes are alpha(2)beta(2) tetramers, the beta subunit being identical to protein disulfide isomerase (PDI). The main Caenorhabditis elegans collagen P4H form is an unusual PHY-1/PHY-2/(PDI)(2) mixed tetramer consisting of two types of catalytic alpha subunit, but the PHY-1 and PHY-2 polypeptides also form active PHY/PDI dimers. The lengths of peptide substrates have a major effect on their interaction with the P4H tetramers, the K(m) values decreasing markedly with increasing chain length. This phenomenon has been explained in terms of processive binding of the two catalytic subunits to long peptides. We determined here the K(m) values of a collagen P4H having two catalytic sites, the C. elegans mixed tetramer, and a form having only one such site, the PHY-1/PDI dimer, for peptides of varying lengths. All the K(m) values of the PHY-1/PDI dimer were found to be about 1.5-2.5 times those of the tetramer, but increasing peptide length led to identical decreases in the values of both enzyme forms. The K(m) for a nonhydroxylated collagen fragment with 33 -X-Y-Gly-triplets but only 11 -X-Pro-Gly-triplets was found to correspond to the number of the former rather than the latter. To study the individual roles of the two catalytic sites in a tetramer, we produced mutant PHY-1/PHY-2/(PDI)(2) tetramers in which binding of the Fe(2+) ion or 2-oxoglutarate to one of the two catalytic sites was prevented. The activities of the mutant tetramers decreased to markedly less than 50% of that of the wild type, being about 5-10% and 20-30% with the enzymes having one of the two Fe(2+)-binding sites or 2-oxoglutarate-binding sites inactivated, respectively, while the K(m) values for these cosubstrates or peptide substrates were not affected. Our data thus indicate that although collagen P4Hs do not act on peptide substrates by a processive mechanism, prevention of hydroxylation at one of the two catalytic sites in the tetramer impairs the function of the other catalytic site.     

Cytovillin and other microvillar proteins of human choriocarcinoma cells

Microvilli were isolated from cultured human JEG-3 choriocarcinoma cells using a gentle shearing method. The protein components of the isolated microvilli were examined by sodium dodecylsulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and immunoblotting. The major Mr 42,000 and Mr 100,000 polypeptide bands reacted with anti-actin and anti-alpha-actinin antisera, respectively. Extraction of the isolated JEG-3 microvilli with Triton X-100 left an insoluble cytoskeletal residue containing mainly actin, alpha-actin, and polypeptides of Mr 200,000, 55,000 and 35,000. The Mr 35,000 polypeptide remained insoluble only at high concentrations of free Ca2+. Immunoblotting analysis of the JEG-3 microvilli indicated that they were devoid of tropomyosin, although the total JEG-3 protein lysates gave a strong positive reaction with anti-tropomyosin antiserum. The different subcellular localization of cytovillin and tropomyosin was also shown by indirect immunofluorescence microscopy. Cytovillin, an Mr 75,000 microvillus-specific membrane protein of JEG-3 cells, existed in an oligomeric form (dimer or trimer) as shown by gel filtration of Triton X-100 solubilized microvillar proteins and by native polyacrylamide gel electrophoresis of purified cytovillin. Disulfide bridges were not involved in the aggregation, because the mobility of cytovillin was similar under reducing and nonreducing conditions in SDS-PAGE. Cytovillin was shown to be closely related to ezrin, a minor component of chicken intestinal brush border microvilli.     

Comparative analysis of some aspects of mitochondrial metabolism in differentiated and undifferentiated neuroblastoma cells

The aim of the present study is to clarify some aspects of the mechanisms of regulation of mitochondrial metabolism in neuroblastoma (NB) cells. Experiments were performed on murine Neuro-2a (N2a) cell line, and the same cells differentiated by all-trans-retinoic acid (dN2a) served as in vitro model of normal neurons. Oxygraphy and Metabolic Control Analysis (MCA) were applied to characterize the function of mitochondrial oxidative phosphorylation (OXPHOS) in NB cells. Flux control coefficients (FCCs) for components of the OXPHOS system were determined using titration studies with specific non-competitive inhibitors in the presence of exogenously added ADP. Respiration rates of undifferentiated Neuro-2a cells (uN2a) and the FCC of Complex-II in these cells were found to be considerably lower than those in dN2a cells. Our results show that NB is not an exclusively glycolytic tumor and could produce a considerable part of ATP via OXPHOS. Two important enzymes - hexokinase-2 and adenylate kinase-2 can play a role in the generation of ATP in NB cells. MCA has shown that in uN2a cells the key sites in the regulation of OXPHOS are complexes I, II and IV, whereas in dN2a cells complexes II and IV. Results obtained for the phosphate and adenine nucleotide carriers showed that in dN2a cells these carriers exerted lower control over the OXPHOS than in undifferentiated cells. The sum of FCCs for both types of NB cells was found to exceed significantly that for normal cells suggesting that in these cells the respiratory chain was somehow reorganized or assembled into large supercomplexes.     

Direct measurement of energy fluxes from mitochondria into cytoplasm in permeabilized cardiac cells <Emphasis Type="Italic">in situ:</Emphasis> some evidence for mitochondrial interactosome

The aim of this study was to measure energy fluxes from mitochondria in isolated permeabilized cardiomyocytes. Respiration of permeabilized cardiomyocytes and mitochondrial membrane potential were measured in presence of MgATP, pyruvate kinase – phosphoenolpyruvate and creatine. ATP and phosphocreatine concentrations in medium surrounding cardiomyocytes were determined. While ATP concentration did not change in time, mitochondria effectively produced phosphocreatine (PCr) with PCr/O₂ ratio equal to 5.68 ± 0.14. Addition of heterodimeric tubulin to isolated mitochondria was found to increase apparent Km for exogenous ADP from 11 ± 2 μM to 330 ± 47 μM, but creatine again decreased it to 23 ± 6 μM. These results show directly that under physiological conditions the major energy carrier from mitochondria into cytoplasm is PCr, produced by mitochondrial creatine kinase (MtCK), which functional coupling to adenine nucleotide translocase is enhanced by selective limitation of permeability of mitochondrial outer membrane within supercomplex ATP Synthasome-MtCK-VDAC-tubulin, Mitochondrial Interactosome.