Single-Cell RNA-Seq Reveals Lineage and X Chromosome Dynamics in Human Preimplantation Embryos.

Mouse studies have been instrumental in forming our current understanding of early cell-lineage decisions; however, similar insights into the early human development are severely limited. Here, we present a comprehensive transcriptional map of human embryo development, including the sequenced transcriptomes of 1,529 individual cells from 88 human preimplantation embryos. These data show that cells undergo an intermediate state of co-expression of lineage-specific genes, followed by a concurrent establishment of the trophectoderm, epiblast, and primitive endoderm lineages, which coincide with blastocyst formation. Female cells of all three lineages achieve dosage compensation of X chromosome RNA levels prior to implantation. However, in contrast to the mouse, XIST is transcribed from both alleles throughout the progression of this expression dampening, and X chromosome genes maintain biallelic expression while dosage compensation proceeds. We envision broad utility of this transcriptional atlas in future studies on human development as well as in stem cell research.

Identification of spatial expression trends in single-cell gene expression data.

As methods for measuring spatial gene expression at single-cell resolution become available, there is a need for computational analysis strategies. We present trendsceek, a method based on marked point processes that identifies genes with statistically significant spatial expression trends. trendsceek finds these genes in spatial transcriptomic and sequential fluorescence in situ hybridization data, and also reveals significant gene expression gradients and hot spots in low-dimensional projections of dissociated single-cell RNA-seq data.

scphaser: haplotype inference using single-cell RNA-seq data.

UNLABELLED: Determination of haplotypes is important for modelling the phenotypic consequences of genetic variation in diploid organisms, including cis-regulatory control and compound heterozygosity. We realized that single-cell RNA-seq (scRNA-seq) data are well suited for phasing genetic variants, since both transcriptional bursts and technical bottlenecks cause pronounced allelic fluctuations in individual single cells. Here we present scphaser, an R package that phases alleles at heterozygous variants to reconstruct haplotypes within transcribed regions of the genome using scRNA-seq data. The devised method efficiently and accurately reconstructed the known haplotype for ≥93% of phasable genes in both human and mouse. It also enables phasing of rare and de novo variants and variants far apart within genes, which is hard to attain with population-based computational inference. AVAILABILITY AND IMPLEMENTATION: scphaser is implemented as an R package. Tutorial and code are available at https://github.com/edsgard/scphaser CONTACT: rickard.sandberg@ki.se SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The role of germline alterations in the DNA damage response genes BRIP1 and BRCA2 in melanoma susceptibility.

We applied a targeted sequencing approach to identify germline mutations conferring a moderately to highly increased risk of cutaneous and uveal melanoma. Ninety-two high-risk melanoma patients were screened for inherited variation in 120 melanoma candidate genes. Observed gene variants were filtered based on frequency in reference populations, cosegregation with melanoma in families and predicted functional effect. Several novel or rare genetic variants in genes involved in DNA damage response, cell-cycle regulation and transcriptional control were identified in melanoma patients. Among identified genetic alterations was an extremely rare variant (minor allele frequency of 0.00008) in the BRIP1 gene that was found to cosegregate with the melanoma phenotype. We also found a rare nonsense variant in the BRCA2 gene (rs11571833), previously associated with cancer susceptibility but not with melanoma, which showed weak association with melanoma susceptibility in the Swedish population. Our results add to the growing knowledge about genetic factors associated with melanoma susceptibility and also emphasize the role of DNA damage response as an important factor in melanoma etiology. © 2016 Wiley Periodicals, Inc.

HExpoChem: a systems biology resource to explore human exposure to chemicals.

SUMMARY: Humans are exposed to diverse hazardous chemicals daily. Although an exposure to these chemicals is suspected to have adverse effects on human health, mechanistic insights into how they interact with the human body are still limited. Therefore, acquisition of curated data and development of computational biology approaches are needed to assess the health risks of chemical exposure. Here we present HExpoChem, a tool based on environmental chemicals and their bioactivities on human proteins with the objective of aiding the qualitative exploration of human exposure to chemicals. The chemical-protein interactions have been enriched with a quality-scored human protein-protein interaction network, a protein-protein association network and a chemical-chemical interaction network, thus allowing the study of environmental chemicals through formation of protein complexes and phenotypic outcomes enrichment. AVAILABILITY: HExpoChem is available at http://www.cbs.dtu.dk/services/HExpoChem-1.0/.

Hereditary uveal melanoma: a report of a germline mutation in BAP1.

Melanoma of the eye is a rare and distinct subtype of melanoma, which only rarely are familial. However, cases of uveal melanoma (UM) have been found in families with mixed cancer syndromes. Here, we describe a comprehensive search for inherited genetic variation in a family with multiple cases of UM but no aggregation of other cancer diagnoses. The proband is a woman diagnosed with UM at 16 years who within 6 months developed liver metastases. We also identified two older paternal relatives of the proband who had died from UM. We performed exome sequencing of germline DNA from members of the affected family. Exome-wide analysis identified a novel loss-of-function mutation in the BAP1 gene, previously suggested as a tumor suppressor. The mutation segregated with the UM phenotype in this family, and we detected a loss of the wild-type allele in the UM tumor of the proband, strongly supporting a causative association with UM. Screening of BAP1 germline mutations in families predisposed for UM may be used to identify individuals at increased risk of disease. Such individuals may then be enrolled in preventive programs and regular screenings to facilitate early detection and thereby improve prognosis.

A genome-wide association study of men with symptoms of testicular dysgenesis syndrome and its network biology interpretation.

BACKGROUND: Testicular dysgenesis syndrome (TDS) is a common disease that links testicular germ cell cancer, cryptorchidism and some cases of hypospadias and male infertility with impaired development of the testis. The incidence of these disorders has increased over the last few decades, and testicular cancer now affects 1% of the Danish and Norwegian male population. METHODS: To identify genetic variants that span the four TDS phenotypes, the authors performed a genome-wide association study (GWAS) using Affymetrix Human SNP Array 6.0 to screen 488 patients with symptoms of TDS and 439 selected controls with excellent reproductive health. Furthermore, they developed a novel integrative method that combines GWAS data with other TDS-relevant data types and identified additional TDS markers. The most significant findings were replicated in an independent cohort of 671 Nordic men. RESULTS: Markers located in the region of TGFBR3 and BMP7 showed association with all TDS phenotypes in both the discovery and replication cohorts. An immunohistochemistry investigation confirmed the presence of transforming growth factor ß receptor type III (TGFBR3) in peritubular and Leydig cells, in both fetal and adult testis. Single-nucleotide polymorphisms in the KITLG gene showed significant associations, but only with testicular cancer. CONCLUSIONS: The association of single-nucleotide polymorphisms in the TGFBR3 and BMP7 genes, which belong to the transforming growth factor ß signalling pathway, suggests a role for this pathway in the pathogenesis of TDS. Integrating data from multiple layers can highlight findings in GWAS that are biologically relevant despite having border significance at currently accepted statistical levels.

ChemProt: a disease chemical biology database.

Systems pharmacology is an emergent area that studies drug action across multiple scales of complexity, from molecular and cellular to tissue and organism levels. There is a critical need to develop network-based approaches to integrate the growing body of chemical biology knowledge with network biology. Here, we report ChemProt, a disease chemical biology database, which is based on a compilation of multiple chemical-protein annotation resources, as well as disease-associated protein-protein interactions (PPIs). We assembled more than 700,000 unique chemicals with biological annotation for 30,578 proteins. We gathered over 2-million chemical-protein interactions, which were integrated in a quality scored human PPI network of 428,429 interactions. The PPI network layer allows for studying disease and tissue specificity through each protein complex. ChemProt can assist in the in silico evaluation of environmental chemicals, natural products and approved drugs, as well as the selection of new compounds based on their activity profile against most known biological targets, including those related to adverse drug events. Results from the disease chemical biology database associate citalopram, an antidepressant, with osteogenesis imperfect and leukemia and bisphenol A, an endocrine disruptor, with certain types of cancer, respectively. The server can be accessed at http://www.cbs.dtu.dk/services/ChemProt/.

Zfp697 is an RNA-binding protein that regulates skeletal muscle inflammation and regeneration.

Muscular atrophy is a mortality risk factor that happens with disuse, chronic disease, and aging. Recovery from atrophy requires changes in several cell types including muscle fibers, and satellite and immune cells. Here we show that Zfp697/ZNF697 is a damage-induced regulator of muscle regeneration, during which its expression is transiently elevated. Conversely, sustained Zfp697 expression in mouse muscle leads to a gene expression signature of chemokine secretion, immune cell recruitment, and extracellular matrix remodeling. Myofiber-specific Zfp697 ablation hinders the inflammatory and regenerative response to muscle injury, compromising functional recovery. We uncover Zfp697 as an essential interferon gamma mediator in muscle cells, interacting primarily with ncRNAs such as the pro-regenerative miR-206. In sum, we identify Zfp697 as an integrator of cell-cell communication necessary for tissue regeneration.

Transcriptomes of the desiccation-tolerant resurrection plant Craterostigma plantagineum.

Studies of the resurrection plant Craterostigma plantagineum have revealed some of the mechanisms which these desiccation-tolerant plants use to survive environments with extreme dehydration and restricted seasonal water. Most resurrection plants are polyploid with large genomes, which has hindered efforts to obtain whole genome sequences and perform mutational analysis. However, the application of deep sequencing technologies to transcriptomics now permits large-scale analyses of gene expression patterns despite the lack of a reference genome. Here we use pyro-sequencing to characterize the transcriptomes of C. plantagineum leaves at four stages of dehydration and rehydration. This reveals that genes involved in several pathways, such as those required for vitamin K and thiamin biosynthesis, are tightly regulated at the level of gene expression. Our analysis also provides a comprehensive picture of the array of cellular responses controlled by gene expression that allow resurrection plants to survive desiccation.

GeneiASE: Detection of condition-dependent and static allele-specific expression from RNA-seq data without haplotype information.

Allele-specific expression (ASE) is the imbalance in transcription between maternal and paternal alleles at a locus and can be probed in single individuals using massively parallel DNA sequencing technology. Assessing ASE within a single sample provides a static picture of the ASE, but the magnitude of ASE for a given transcript may vary between different biological conditions in an individual. Such condition-dependent ASE could indicate a genetic variation with a functional role in the phenotypic difference. We investigated ASE through RNA-sequencing of primary white blood cells from eight human individuals before and after the controlled induction of an inflammatory response, and detected condition-dependent and static ASE at 211 and 13021 variants, respectively. We developed a method, GeneiASE, to detect genes exhibiting static or condition-dependent ASE in single individuals. GeneiASE performed consistently over a range of read depths and ASE effect sizes, and did not require phasing of variants to estimate haplotypes. We observed condition-dependent ASE related to the inflammatory response in 19 genes, and static ASE in 1389 genes. Allele-specific expression was confirmed by validation of variants through real-time quantitative RT-PCR, with RNA-seq and RT-PCR ASE effect-size correlations r = 0.67 and r = 0.94 for static and condition-dependent ASE, respectively.

Using Whole-Exome Sequencing to Identify Genetic Markers for Carboplatin and Gemcitabine-Induced Toxicities.

PURPOSE: Chemotherapies are associated with significant interindividual variability in therapeutic effect and adverse drug reactions. In lung cancer, the use of gemcitabine and carboplatin induces grade 3 or 4 myelosuppression in about a quarter of the patients, while an equal fraction of patients is basically unaffected in terms of myelosuppressive side effects. We therefore set out to identify genetic markers for gemcitabine/carboplatin-induced myelosuppression. EXPERIMENTAL DESIGN: We exome sequenced 32 patients that suffered extremely high neutropenia and thrombocytopenia (grade 3 or 4 after first chemotherapy cycle) or were virtually unaffected (grade 0 or 1). The genetic differences/polymorphism between the groups were compared using six different bioinformatics strategies: (i) whole-exome nonsynonymous single-nucleotide variants association analysis, (ii) deviation from Hardy-Weinberg equilibrium, (iii) analysis of genes selected by a priori biologic knowledge, (iv) analysis of genes selected from gene expression meta-analysis of toxicity datasets, (v) Ingenuity Pathway Analysis, and (vi) FunCoup network enrichment analysis. RESULTS: A total of 53 genetic variants that differed among these groups were validated in an additional 291 patients and were correlated to the patients' myelosuppression. In the validation, we identified rs1453542 in OR4D6 (P = 0.0008; OR, 5.2; 95% CI, 1.8-18) as a marker for gemcitabine/carboplatin-induced neutropenia and rs5925720 in DDX53 (P = 0.0015; OR, 0.36; 95% CI, 0.17-0.71) as a marker for thrombocytopenia. Patients homozygous for the minor allele of rs1453542 had a higher risk of neutropenia, and for rs5925720 the minor allele was associated with a lower risk for thrombocytopenia. CONCLUSIONS: We have identified two new genetic markers with the potential to predict myelosuppression induced by gemcitabine/carboplatin chemotherapy.

Over Expression of NANOS3 and DAZL in Human Embryonic Stem Cells.

The mechanisms underlying human germ cell development are largely unknown, partly due to the scarcity of primordial germ cells and the inaccessibility of the human germline to genetic analysis. Human embryonic stem cells can differentiate to germ cells in vitro and can be genetically modified to study the genetic requirements for germ cell development. Here, we studied NANOS3 and DAZL, which have critical roles in germ cell development in several species, via their over expression in human embryonic stem cells using global transcriptional analysis, in vitro germ cell differentiation, and in vivo germ cell formation assay by xenotransplantation. We found that NANOS3 over expression prolonged pluripotency and delayed differentiation. In addition, we observed a possible connection of NANOS3 with inhibition of apoptosis. For DAZL, our results suggest a post-transcriptional regulation mechanism in hES cells. In addition, we found that DAZL suppressed the translation of OCT4, and affected the transcription of several genes associated with germ cells, cell cycle arrest, and cell migration. Furthermore, DAZL over expressed cells formed spermatogonia-like colonies in a rare instance upon xenotransplantation. These data can be used to further elucidate the role of NANOS3 and DAZL in germ cell development both in vitro and in vivo.

Genome-wide assessment of the association of rare and common copy number variations to testicular germ cell cancer.

Testicular germ cell cancer (TGCC) is one of the most heritable forms of cancer. Previous genome-wide association studies have focused on single nucleotide polymorphisms, largely ignoring the influence of copy number variants (CNVs). Here we present a genome-wide study of CNV on a cohort of 212 cases and 437 controls from Denmark, which was genotyped at ∼1.8 million markers, half of which were non-polymorphic copy number markers. No association of common variants were found, whereas analysis of rare variants (present in less than 1% of the samples) initially indicated a single gene with significantly higher accumulation of rare CNVs in cases as compared to controls, at the gene PTPN1 (P = 3.8 × 10(-2), 0.9% of cases and 0% of controls). However, the CNV could not be verified by qPCR in the affected samples. Further, the CNV calling of the array-data was validated by sequencing of the GSTM1 gene, which showed that the CNV frequency was in complete agreement between the two platforms. This study therefore disconfirms the hypothesis that there exists a single CNV locus with a major effect size that predisposes to TGCC. Genome-wide pathway association analysis indicated a weak association of rare CNVs related to cell migration (false-discovery rate = 0.021, 1.8% of cases and 1.1% of controls). Dysregulation during migration of primordial germ cells has previously been suspected to be a part of TGCC development and this set of multiple rare variants may thereby have a minor contribution to an increased susceptibility of TGCCs.

Analysis of gene expression profiles of microdissected cell populations indicates that testicular carcinoma in situ is an arrested gonocyte.

Testicular germ cell cancers in young adult men derive from a precursor lesion called carcinoma in situ (CIS) of the testis. CIS cells were suggested to arise from primordial germ cells or gonocytes. However, direct studies on purified samples of CIS cells are lacking. To overcome this problem, we performed laser microdissection of CIS cells. Highly enriched cell populations were obtained and subjected to gene expression analysis. The expression profile of CIS cells was compared with microdissected gonocytes, oogonia, and cultured embryonic stem cells with and without genomic aberrations. Three samples of each tissue type were used for the analyses. Unique expression patterns for these developmentally very related cell types revealed that CIS cells were very similar to gonocytes because only five genes distinguished these two cell types. We did not find indications that CIS was derived from a meiotic cell, and the similarity to embryonic stem cells was modest compared with gonocytes. Thus, we provide new evidence that the molecular phenotype of CIS cells is similar to that of gonocytes. Our data are in line with the idea that CIS cells may be gonocytes that survived in the postnatal testis. We speculate that disturbed development of somatic cells in the fetal testis may play a role in allowing undifferentiated cells to survive in the postnatal testes. The further development of CIS into invasive germ cell tumors may depend on signals from their postpubertal niche of somatic cells, including hormones and growth factors from Leydig and Sertoli cells.