Ambiguous genes due to aligners and their impact on RNA-seq data analysis.

The main scope of the study is ambiguous genes, i.e. genes whose expression is difficult to estimate from the data produced by next-generation sequencing technologies. We focused on the RNA sequencing (RNA-Seq) type of experiment performed on the Illumina platform. It is crucial to identify such genes and understand the cause of their difficulty, as these genes may be involved in some diseases. By giving misleading results, they could contribute to a misunderstanding of the cause of certain diseases, which could lead to inappropriate treatment. We thought that the ambiguous genes would be difficult to map because of their complex structure. So we looked at RNA-seq analysis using different mappers to find genes that would have different measurements from the aligners. We were able to identify such genes using a generalized linear model with two factors: mappers and groups introduced by the experiment. A large proportion of ambiguous genes are pseudogenes. High sequence similarity of pseudogenes to functional genes may indicate problems in alignment procedures. In addition, predictive analysis verified the performance of difficult genes in classification. The effectiveness of classifying samples into specific groups was compared, including the expression of difficult and not difficult genes as covariates. In almost all cases considered, ambiguous genes have less predictive power.

Temporal and sequential transcriptional dynamics define lineage shifts in corticogenesis.

The cerebral cortex contains billions of neurons, and their disorganization or misspecification leads to neurodevelopmental disorders. Understanding how the plethora of projection neuron subtypes are generated by cortical neural stem cells (NSCs) is a major challenge. Here, we focused on elucidating the transcriptional landscape of murine embryonic NSCs, basal progenitors (BPs), and newborn neurons (NBNs) throughout cortical development. We uncover dynamic shifts in transcriptional space over time and heterogeneity within each progenitor population. We identified signature hallmarks of NSC, BP, and NBN clusters and predict active transcriptional nodes and networks that contribute to neural fate specification. We find that the expression of receptors, ligands, and downstream pathway components is highly dynamic over time and throughout the lineage implying differential responsiveness to signals. Thus, we provide an expansive compendium of gene expression during cortical development that will be an invaluable resource for studying neural developmental processes and neurodevelopmental disorders.

Chronic Leptin Treatment Induces Epithelial-Mesenchymal Transition in MCF10A Mammary Epithelial Cells.

Leptin is a cytokine-like hormone that functions as a link between obesity and breast cancer (BC). Leptin treatment induces Epithelial to Mesenchymal Transition (EMT) in BC cell lines. In non-tumoral breast epithelial MCF10A cells, acute leptin treatment induces partial EMT. However, the effect of chronic leptin treatment on EMT in non-tumorigenic breast cells has not been fully explored. This study aimed to evaluate the effect of chronic leptin treatment on the induction of EMT in MCF10A cells. We found that chronic leptin treatment induces a switch from an epithelial to a mesenchymal morphology, partial loss of E-cadherin and gain of vimentin expression. Immunolocalization experiments showed a partial loss of E-cadherin at cell junctions and increased cytoplasmic localization of vimentin in leptin-treated cells. Moreover, chronic leptin treatment increased collective cell migration and invasion. Furthermore, when cultured in non-adherent conditions leptin treated cells exhibited reduced cell aggregation, increased survival, and decreased apoptosis, which correlates with increased FAK and AKT phosphorylation. Finally, bioinformatic analysis in two publicly available RNAseq datasets from normal breast tissue shows that high levels of leptin mRNA correlate positively with the expression of mesenchymal markers, and negatively with epithelial markers. Thus, our results demonstrate that chronic leptin treatment induces EMT in non-tumorigenic MCF10A cells and suggest that high leptin expression in normal breast tissue may induce EMT and contribute to increased risk of breast cancer.

Heterogeneous expression of ACE2 and TMPRRS2 in mesenchymal stromal cells.

The outbreak of COVID-19 has become a serious public health emergency. The virus targets cells by binding the ACE2 receptor. After infection, the virus triggers in some humans an immune storm containing the release of proinflammatory cytokines and chemokines followed by multiple organ failure. Several vaccines are enrolled, but an effective treatment is still missing. Mesenchymal stem cells (MSCs) have shown to secrete immunomodulatory factors that suppress this cytokine storm. Therefore, MSCs have been suggested as a potential treatment option for COVID-19. We report here that the ACE2 expression is minimal or nonexistent in MSC derived from three different human tissue sources (adipose tissue, umbilical cord Wharton`s jelly and bone marrow). In contrast, TMPRSS2 that is implicated in SARS-CoV-2 entry has been detected in all MSC samples. These results are of particular importance for future MSC-based cell therapies to treat severe cases after COVID-19 infection.

SwissPKcdw - A clinical data warehouse for the optimization of pediatric dosing regimens.

Clinical trials have been performed mainly in adults and accordingly the necessary information is lacking for pediatric patients, especially regarding dosage recommendation for approved drugs. This gap in information could be filled with results from pharmacokinetic (PK) modeling, based on data collected in daily clinical routine. In order to make this data accessible and usable for research, the Swiss Pharmacokinetics Clinical Data Warehouse (SwissPKcdw ) project has been set up, including a clinical data warehouse (CDW) and the regulatory framework for data transfer and use within. Embedded into the secure BioMedIT network, the CDW can connect to various data providers and researchers in order to collaborate on the data securely. Due to its modularity, partially containerized deployment and open-source software, each of the components can be extended, modified, and re-used for similar projects that require integrated data management, data analysis, and web tools in a secure scientific data and information technology (IT) environment. Here, we describe a collaborative and interprofessional effort to implement the aforementioned infrastructure between several partners from medical health care and academia. Furthermore, we describe a real-world use case where blood samples from pediatric patients were analyzed for the presence of genetic polymorphisms and the results were aggregated and further analyzed together with the health-related patient data in the SwissPKcdw .

Estradiol Inhibits Human Brain Vascular Pericyte Migration Activity: A Functional and Transcriptomic Analysis.

Stroke is the third leading cause of mortality in women and it kills twice as many women as breast cancer. A key role in the pathophysiology of stroke plays the disruption of the blood-brain barrier (BBB) within the neurovascular unit. While estrogen induces vascular protective actions, its influence on stroke remains unclear. Moreover, experiments assessing its impact on endothelial cells to induce barrier integrity are non-conclusive. Since pericytes play an active role in regulating BBB integrity and function, we hypothesize that estradiol may influence BBB by regulating their activity. In this study using human brain vascular pericytes (HBVPs) we investigated the impact of estradiol on key pericyte functions known to influence BBB integrity. HBVPs expressed estrogen receptors (ER-α, ER-ß and GPER) and treatment with estradiol (10 nM) inhibited basal cell migration but not proliferation. Since pericyte migration is a hallmark for BBB disruption following injury, infection and inflammation, we investigated the effects of estradiol on TNFα-induced PC migration. Importantly, estradiol prevented TNFα-induced pericyte migration and this effect was mimicked by PPT (ER-α agonist) and DPN (ER-ß agonist), but not by G1 (GPR30 agonist). The modulatory effects of estradiol were abrogated by MPP and PHTPP, selective ER-α and ER-ß antagonists, respectively, confirming the role of ER-α and ER-ß in mediating the anti-migratory actions of estrogen. To delineate the intracellular mechanisms mediating the inhibitory actions of estradiol on PC migration, we investigated the role of AKT and MAPK activation. While estradiol consistently reduced the TNFα-induced MAPK and Akt phosphorylation, only the inhibition of MAPK, but not Akt, significantly abrogated the migratory actions of TNFα. In transendothelial electrical resistance measurements, estradiol induced barrier function (TEER) in human brain microvascular endothelial cells co-cultured with pericytes, but not in HBMECs cultured alone. Importantly, transcriptomics analysis of genes modulated by estradiol in pericytes showed downregulation of genes known to increase cell migration and upregulation of genes known to inhibit cell migration. Taken together, our findings provide the first evidence that estradiol modulates pericyte activity and thereby improves endothelial integrity.

Transcryptomic Analysis of Human Brain -Microvascular Endothelial Cell Driven Changes in -Vascular Pericytes.

Many pathological conditions of the brain are associated with structural abnormalities within the neurovascular system and linked to pericyte (PC) loss and/or dysfunction. Since crosstalk between endothelial cells (ECs) and PCs greatly impacts the function of the blood-brain barrier (BBB), effects of PCs on endothelial integrity and function have been investigated extensively. However, the impact of ECs on the function and activity of PCs remains largely unknown. Hence, using co-cultures of human brain vascular PCs with human cerebral microvascular ECs on opposite sides of porous Transwell inserts which facilitates direct EC-PC contact and improves EC barrier function, we analyzed EC-driven transcriptomic changes in PCs using microarrays and changes in cytokines/chemokines using proteome arrays. Gene expression analysis (GEA) in PCs co-cultured with ECs versus PCs cultured alone showed significant upregulation of 1'334 genes and downregulation of 964 genes. GEA in co-cultured PCs revealed increased expression of five prominent PC markers as well as soluble factors, such as transforming growth factor beta, fibroblast growth factor, angiopoietin 1, brain-derived neurotrophic factor, all of which are involved in EC-PC crosstalk and BBB induction. Pathway enrichment analysis of modulated genes showed a strong impact on many inflammatory and extracellular matrix (ECM) pathways including interferon and interleukin signaling, TGF-ß and interleukin-1 regulation of ECM, as well as on the mRNA processing pathway. Interestingly, while co-culture induced the mRNA expression of many chemokines and cytokines, including several CCL- and CXC-motif ligands and interleukins, we observed a decreased expression of the same inflammatory mediators on the protein level. Importantly, in PCs, ECs significantly induced interferon associated proteins (IFIT1, IFI44L, IF127, IFIT3, IFI6, IFI44) with anti-viral actions; downregulated prostaglandin E receptor 2 (prevent COX-2 mediated BBB damage); upregulated fibulin-3 and connective tissue growth factor essential for BBB integrity; and multiple ECMs (collagens and integrins) that inhibit cell migration. Our findings suggest that via direct contact, ECs prime PCs to induce molecules to promote BBB integrity and cell survival during infection and inflammatory insult. Taken together, we provide first evidence that interaction with ECs though porous membranes induces major changes in the transcriptomic and proteomic profile of PCs. ECs influence genes involved in diverse aspects of PC function including PC maturation, cell survival, anti-viral defense, blood flow regulation, immuno-modulation and ECM deposition.

Transcryptomic Analysis of Human Brain-Microvascular Endothelial Response to -Pericytes: Cell Orientation Defines Barrier Function.

Pericytes facilitate blood-brain barrier (BBB) integrity; however, the mechanisms involved remain unclear. Hence, using co-cultures of human cerebral microvascular endothelial cells (ECs) and vascular pericytes (PCs) in different spatial arrangements, as well as PC conditioned media, we investigated the impact of PC-EC orientation and PC-derived soluble factors on EC barrier function. We provide the first evidence that barrier-inducing properties of PCs require basolateral contact with ECs. Gene expression analysis (GEA) in ECs co-cultured with PCs versus ECs alone showed significant upregulation of 38 genes and downregulation of 122 genes. Pathway enrichment analysis of modulated genes showed significant regulation of several pathways, including transforming growth factor-ß and interleukin-1 regulated extracellular matrix, interferon and interleukin signaling, immune system signaling, receptor of advanced glycation end products (RAGE), and cytokine-cytokine receptor interaction. Transcriptomic analysis showed a reduction in molecules such as pro-inflammatory cytokines and chemokines, which are known to be induced during BBB disruption. Moreover, cytokine proteome array confirmed the downregulation of key pro-inflammatory cytokines and chemokines on the protein level. Other molecules which influence BBB and were favorably modulated upon EC-PC co-culture include IL-18 binding protein, kallikrein-3, CSF2 CSF3, CXCL10, CXCL11 (downregulated) and IL-1-R4; HGF, PDGF-AB/BB, PECAM, SERPIN E1 (upregulated). In conclusion, we provide the first evidence that (1) basolateral contact between ECs and PCs is essential for EC barrier function and integrity; (2) in ECs co-cultured with PCs, the profile of BBB disrupting pro-inflammatory molecules and cytokines/chemokines is downregulated; (3) PCs significantly modulate EC mechanisms known to improve barrier function, including TGF-ß regulated ECM pathway, anti-inflammatory cytokines, growth factors and matrix proteins. This human PC-EC co-culture may serve as a viable in vitro model for investigating BBB function and drug transport.

Structure of SRSF1 RRM1 bound to RNA reveals an unexpected bimodal mode of interaction and explains its involvement in SMN1 exon7 splicing.

The human prototypical SR protein SRSF1 is an oncoprotein that contains two RRMs and plays a pivotal role in RNA metabolism. We determined the structure of the RRM1 bound to RNA and found that the domain binds preferentially to a CN motif (N is for any nucleotide). Based on this solution structure, we engineered a protein containing a single glutamate to asparagine mutation (E87N), which gains the ability to bind to uridines and thereby activates SMN exon7 inclusion, a strategy that is used to cure spinal muscular atrophy. Finally, we revealed that the flexible inter-RRM linker of SRSF1 allows RRM1 to bind RNA on both sides of RRM2 binding site. Besides revealing an unexpected bimodal mode of interaction of SRSF1 with RNA, which will be of interest to design new therapeutic strategies, this study brings a new perspective on the mode of action of SRSF1 in cells.

Temporal activation of WNT/ß-catenin signaling is sufficient to inhibit SOX10 expression and block melanoma growth.

Despite advances in the systemic treatment of patients with metastatic melanoma using immune checkpoint and tyrosine kinase inhibitors (TKI), the majority of stage IV melanoma patients eventually succumb to the disease. We have previously identified the transcription factor Sox10 as a crucial player in melanoma, yet the underlying molecular mechanisms mediating Sox10-dependent tumorigenesis remain largely uncharacterized. Here, we show that MEK and RAF inhibitors do not suppress levels of SOX10 protein in patient-derived cells in vitro, as well as in melanoma patients in vivo. In a search for pharmacological inhibitors of SOX10, we performed a mass spectrometry-based screen in human melanoma cells. Subsequent analysis revealed that SOX10 directly interacts with ß-catenin, which is a key mediator of canonical Wnt/ß-catenin signaling. We demonstrate that inhibitors of glycogen synthase kinase 3 alpha/beta (GSK3α/ß) efficiently abrogate SOX10 protein in human melanoma cells in vitro and in melanoma mouse models in vivo. The mechanism of action of GSK3-mediated SOX10 suppression is transcription-independent and relies on the presence of a proteasome degradable form of ß-catenin. Taken together, we provide evidence that activation of canonical Wnt signaling has a profound effect on melanoma growth and is able to counteract Sox10-dependent melanoma maintenance both in vitro and in vivo.

Recording transcriptional histories using Record-seq.

It is difficult to elucidate the transcriptional history of a cell using current experimental approaches, as they are destructive in nature and therefore describe only a moment in time. To overcome these limitations, we recently established Record-seq, a technology that enables transcriptional recording by CRISPR spacer acquisition from RNA. The recorded transcriptomes are recovered by SENECA, a method that selectively amplifies expanded CRISPR arrays, followed by deep sequencing. The resulting CRISPR spacers are aligned to the host genome, thereby enabling transcript quantification and associated analyses. Here, we describe the experimental procedures of the Record-seq workflow as well as subsequent data analysis. Beginning with the experimental design, Record-seq data can be obtained and analyzed within 1-2 weeks.

Medullary Respiratory Circuit Is Reorganized by a Seasonally-Induced Program in Preparation for Hibernation.

Deep hibernators go through several cycles of profound drops in body temperature during the winter season, with core temperatures sometimes reaching near freezing. Yet unlike non-hibernating mammals, they can sustain breathing rhythms. The physiological processes that make this possible are still not understood. In this study, we focused on the medullary Ventral Respiratory Column of a facultative hibernator, the Syrian hamster. Using shortened day-lengths, we induced a "winter-adapted" physiological state, which is a prerequisite for hibernation. When recording electrophysiological signals from acute slices in the winter-adapted pre-Bötzinger complex (preBötC), spike trains showed higher spike rates, amplitudes, complexity, as well as higher temperature sensitivity, suggesting an increase in connectivity and/or synaptic strength during the winter season. We further examined action potential waveforms and found that the depolarization integral, as measured by the area under the curve, is selectively enhanced in winter-adapted animals. This suggests that a shift in the ion handling kinetics is also being induced by the winter-adaptation program. RNA sequencing of respiratory pre-motor neurons, followed by gene set enrichment analysis, revealed differential regulation and splicing in structural, synaptic, and ion handling genes. Splice junction analysis suggested that differential exon usage is occurring in a select subset of ion handling subunits (ATP1A3, KCNC3, SCN1B), and synaptic structure genes (SNCB, SNCG, RAB3A). Our findings show that the hamster respiratory center undergoes a seasonally-cued alteration in electrophysiological properties, likely protecting against respiratory failure at low temperatures.

Gene expression in chronic granulomatous disease and interferon-γ receptor-deficient cells treated in vitro with interferon-γ.

Interferon-γ (IFN-γ) plays an important role in innate and adaptive immunity against intracellular infections and is used clinically for the prevention and control of infections in chronic granulomatous disease (CGD) and inborn defects in the IFN-γ/interleukin (IL)-12 axis. Using transcriptome profiling (RNA-seq), we sought to identify differentially expressed genes, transcripts and exons in Epstein-Barr virus-transformed B lymphocytes (B-EBV) cells from CGD patients, IFN-γ receptor deficiency patients, and normal controls, treated in vitro with IFN-γ for 48 hours. Our results show that IFN-γ increased the expression of a diverse array of genes related to different cellular programs. In cells from normal controls and CGD patients, IFN-γ-induced expression of genes relevant to oxidative killing, nitric oxide synthase pathway, proteasome-mediated degradation, antigen presentation, chemoattraction, and cell adhesion. IFN-γ also upregulated genes involved in diverse stages of messenger RNA (mRNA) processing including pre-mRNA splicing, as well as others implicated in the folding, transport, and assembly of proteins. In particular, differential exon expression of WARS (encoding tryptophanyl-transfer RNA synthetase, which has an essential function in protein synthesis) induced by IFN-γ in normal and CGD cells suggests that this gene may have an important contribution to the benefits of IFN-γ treatment for CGD. Upregulation of mRNA and protein processing related genes in CGD and IFNRD cells could mediate some of the effects of IFN-γ treatment. These data support the concept that IFN-γ treatment may contribute to increased immune responses against pathogens through regulation of genes important for mRNA and protein processing.

SeQuiLa: an elastic, fast and scalable SQL-oriented solution for processing and querying genomic intervals.

SUMMARY: Efficient processing of large-scale genomic datasets has recently become possible due to the application of 'big data' technologies in bioinformatics pipelines. We present SeQuiLa-a distributed, ANSI SQL-compliant solution for speedy querying and processing of genomic intervals that is available as an Apache Spark package. Proposed range join strategy is significantly (∼22×) faster than the default Apache Spark implementation and outperforms other state-of-the-art tools for genomic intervals processing. AVAILABILITY AND IMPLEMENTATION: The project is available at http://biodatageeks.org/sequila/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Metrics for Assessing the Similarity of Microwave Breast Imaging Scans of Healthy Volunteers.

Microwave radar imaging is promising as a complementary medical imaging modality. However, the unique nature of the images means interpretation can be difficult. As a result, it is important to understand the sources of image differences, and how much variability is inherent in the imaging system itself. To address this issue, we compare the effectiveness of six different measures of image similarity for quantifying the similarity (or difference) between two microwave radar images. The structural similarity index has become the de facto standard for image comparison, but we propose that useful information can be acquired from a measure known as the Modified Hausdorff Distance. We apply each measure to image pairs from sequential scans of both phantoms and volunteers. We find that rather than using a single value to quantify the image similarity, by computing a number of values that are designed to capture different image aspects, we can better assess the ways in which the images differ.

Genes sharing the protein family domain decrease the performance of classification with RNA-seq genomic signatures.

BACKGROUND: The experience with running various types of classification on the CAMDA neuroblastoma dataset have led us to the conclusion that the results are not always obvious and may differ depending on type of analysis and selection of genes used for classification. This paper aims in pointing out several factors that may influence the downstream machine learning analysis. In particular those factors are: type of the primary analysis, type of the classifier and increased correlation between the genes sharing a protein domain. They influence the analysis directly, but also interplay between them may be important. We have compiled the gene-domain database and used it for analysis to see the differences between the genes that share a domain versus the rest of the genes in the datasets. RESULTS: The major findings are: pairs of genes that share a domain have an increased Spearman's correlation coefficients of counts; genes sharing a domain are expected to have a lower predictive power due to increased correlation. For most of the cases it can be seen with the higher number of misclassified samples; classifiers performance may vary depending on a method, still in most cases using genes sharing a domain in the training set results in a higher misclassification rate; increased correlation in genes sharing a domain results most often in worse performance of the classifiers regardless of the primary analysis tools used, even if the primary analysis alignment yield varies. CONCLUSIONS: The effect of sharing a domain is likely more a results of real biological co-expression than just sequence similarity and artifacts of mapping and counting. Still, this is more difficult to conclude and needs further research. The effect is interesting itself, but we also point out some practical aspects in which it may influence the RNA sequencing analysis and RNA biomarker use. In particular it means that a gene signature biomarker set build out of RNA-sequencing results should be depleted for genes sharing common domains. It may cause to perform better when applying classification. REVIEWERS: This article was reviewed by Dimitar Vassiliev and Susmita Datta.

RNA-Seq analysis during the life cycle of Cryptosporidium parvum reveals significant differential gene expression between proliferating stages in the intestine and infectious sporozoites.

Cryptosporidium parvum is a major cause of diarrhoea in humans and animals. There are no vaccines and few drugs available to control C. parvum. In this study, we used RNA-Seq to compare gene expression in sporozoites and intracellular stages of C. parvum to identify genes likely to be important for successful completion of the parasite's life cycle and, thereby, possible targets for drugs or vaccines. We identified 3774 protein-encoding transcripts in C. parvum. Applying a stringent cut-off of eight fold for determination of differential expression, we identified 173 genes (26 coding for predicted secreted proteins) upregulated in sporozoites. On the other hand, expression of 1259 genes was upregulated in intestinal stages (merozoites/gamonts) with a gene ontology enrichment for 63 biological processes and upregulation of 117 genes in 23 metabolic pathways. There was no clear stage specificity of expression of AP2-domain containing transcription factors, although sporozoites had a relatively small repertoire of these important regulators. Our RNA-Seq analysis revealed a new calcium-dependent protein kinase, bringing the total number of known calcium-dependent protein kinases (CDPKs) in C. parvum to 11. One of these, CDPK1, was expressed in all stages, strengthening the notion that it is a valid drug target. By comparing parasites grown in vivo (which produce bona fide thick-walled oocysts) and in vitro (which are arrested in sexual development prior to oocyst generation) we were able to confirm that genes encoding oocyst wall proteins are expressed in gametocytes and that the proteins are stockpiled rather than generated de novo in zygotes. RNA-Seq analysis of C. parvum revealed genes expressed in a stage-specific manner and others whose expression is required at all stages of development. The functional significance of these can now be addressed through recent advances in transgenics for C. parvum, and may lead to the identification of viable drug and vaccine targets.

Molecular insights into Cassava brown streak virus susceptibility and resistance by profiling of the early host response.

Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV) are responsible for significant cassava yield losses in eastern sub-Saharan Africa. To study the possible mechanisms of plant resistance to CBSVs, we inoculated CBSV-susceptible and CBSV-resistant cassava varieties with a mixed infection of CBSVs using top-cleft grafting. Transcriptome profiling of the two cassava varieties was performed at the earliest time point of full infection (28 days after grafting) in the susceptible scions. The expression of genes encoding proteins in RNA silencing, salicylic acid pathways and callose deposition was altered in the susceptible cassava variety, but transcriptional changes were limited in the resistant variety. In total, the expression of 585 genes was altered in the resistant variety and 1292 in the susceptible variety. Transcriptional changes led to the activation of ß-1,3-glucanase enzymatic activity and a reduction in callose deposition in the susceptible cassava variety. Time course analysis also showed that CBSV replication in susceptible cassava induced a strong up-regulation of RDR1, a gene previously reported to be a susceptibility factor in other potyvirus-host pathosystems. The differences in the transcriptional responses to CBSV infection indicated that susceptibility involves the restriction of callose deposition at plasmodesmata. Aniline blue staining of callose deposits also indicated that the resistant variety displays a moderate, but significant, increase in callose deposition at the plasmodesmata. Transcriptome data suggested that resistance does not involve typical antiviral defence responses (i.e. RNA silencing and salicylic acid). A meta-analysis of the current RNA-sequencing (RNA-seq) dataset and selected potyvirus-host and virus-cassava RNA-seq datasets revealed that the conservation of the host response across pathosystems is restricted to genes involved in developmental processes.

Benchmarking distributed data warehouse solutions for storing genomic variant information.

Database URL: https://github.com/ZSI-Bio/variantsdwh.

Corrigendum: Pramel7 mediates ground-state pluripotency through proteasomal-epigenetic combined pathways.

This corrects the article DOI: 10.1038/ncb3554.