Search

Pluripotent stem cells underpin a growing sector that leverages their differentiation potential for research, industry, and clinical applications. This review evaluates the landscape of methods in single-cell transcriptomics that are enabling accelerated discovery in stem cell science. We focus on strategies for scaling stem cell differentiation through multiplexed single-cell analyses, for evaluating molecular regulation of cell differentiation using new analysis algorithms, and methods for integration and projection analysis to classify and benchmark stem cell derivatives against in vivo cell types. By discussing the available methods, comparing their strengths, and illustrating strategies for developing integrated analysis pipelines, we provide user considerations to inform their implementation and interpretation.

Human pluripotent stem cell (hPSC)-derived progenies are immature versions of cells, presenting a potential limitation to the accurate modelling of disease associated with maturity or age. Hence, it is important to characterise how closely cells used in culture resemble their native counterparts. In order to select appropriate points in time for RPE cultures to reflect native counterparts, we characterised the transcriptomic profiles of hPSC-derived retinal pigment epithelium (RPE) cells from 1- and 12-month cultures. We differentiated the human embryonic stem cell line H9 into RPE cells, performed single cell RNA-sequencing of a total of 16,576 cells, and analysed the resulting data to assess the molecular changes of RPE cells across these two culture time points. Our results indicate the stability of the RPE transcriptomic signature, with no evidence of an epithelial – mesenchymal transition, and with maturing populations of RPE observed with time in culture. Assessment of gene ontology pathways revealed that as cultures age, RPE cells upregulate expression of genes involved in metal binding and antioxidant functions. This might reflect an increased ability to handle oxidative stress as cells mature. Comparison with native human RPE data confirmed a maturing transcriptional profile of RPE cells in culture. These results suggest that in vitro long-term culture of RPE cells allow the modelling of specific phenotypes observed in native mature tissue. Our work highlights the transcriptional landscape of hPSC-derived RPE as they age in culture, which provides a reference for native and patient-samples to be benchmarked against.

The human immune system displays substantial variation between individuals, leading to differences in susceptibility to autoimmune disease. We present single-cell RNA sequencing (scRNA-seq) data from 1,267,758 peripheral blood mononuclear cells from 982 healthy human subjects. For 14 cell types, we identified 26,597 independent cis-expression quantitative trait loci (eQTLs) and 990 trans-eQTLs, with most showing cell type–specific effects on gene expression. We subsequently show how eQTLs have dynamic allelic effects in B cells that are transitioning from naïve to memory states and demonstrate how commonly segregating alleles lead to interindividual variation in immune function. Finally, using a Mendelian randomization approach, we identify the causal route by which 305 risk loci contribute to autoimmune disease at the cellular level. This work brings together genetic epidemiology with scRNA-seq to uncover drivers of interindividual variation in the immune system.

Population health research is increasingly focused on the genetic determinants of healthy ageing, but there is no public resource of whole genome sequences and phenotype data from healthy elderly individuals. Here we describe the first release of the Medical Genome Reference Bank (MGRB), comprising whole genome sequence and phenotype of 2570 elderly Australians depleted for cancer, cardiovascular disease, and dementia. We analyse the MGRB for single-nucleotide, indel and structural variation in the nuclear and mitochondrial genomes. MGRB individuals have fewer disease-associated common and rare germline variants, relative to both cancer cases and the gnomAD and UK Biobank cohorts, consistent with risk depletion. Age-related somatic changes are correlated with grip strength in men, suggesting blood-derived whole genomes may also provide a biologic measure of age-related functional deterioration. The MGRB provides a broadly applicable reference cohort for clinical genetics and genomic association studies, and for understanding the genetics of healthy ageing., Healthspan and healthy aging are areas of research with potential socioeconomic impact. Here, the authors present the Medical Genome Reference Bank (MGRB) which consist of over 4,000 individuals aged 70 years and older without a history of the major age-related diseases and report on results from whole-genome sequencing and association analyses.

Previous genetic association studies have failed to identify loci robustly associated with sepsis, and there have been no published genetic association studies or polygenic risk score analyses of patients with septic shock, despite evidence suggesting genetic factors may be involved. We systematically collected genotype and clinical outcome data in the context of a randomized controlled trial from patients with septic shock to enrich the presence of disease-associated genetic variants. We performed genomewide association studies of susceptibility and mortality in septic shock using 493 patients with septic shock and 2442 population controls, and polygenic risk score analysis to assess genetic overlap between septic shock risk/mortality with clinically relevant traits. One variant, rs9489328, located in AL589740.1 noncoding RNA, was significantly associated with septic shock (p = 1.05 × 10–10); however, it is likely a false-positive. We were unable to replicate variants previously reported to be associated (p < 1.00 × 10–6 in previous scans) with susceptibility to and mortality from sepsis. Polygenic risk scores for hematocrit and granulocyte count were negatively associated with 28-day mortality (p = 3.04 × 10–3; p = 2.29 × 10–3), and scores for C-reactive protein levels were positively associated with susceptibility to septic shock (p = 1.44 × 10–3). Results suggest that common variants of large effect do not influence septic shock susceptibility, mortality and resolution; however, genetic predispositions to clinically relevant traits are significantly associated with increased susceptibility and mortality in septic individuals.

Background Advances in droplet-based single-cell RNA-sequencing (scRNA-seq) have dramatically increased throughput, allowing tens of thousands of cells to be routinely sequenced in a single experiment. In addition to cells, droplets capture cell-free “ambient” RNA predominantly caused by lysis of cells during sample preparation. Samples with high ambient RNA concentration can create challenges in accurately distinguishing cell-containing droplets and droplets containing ambient RNA. Current methods to separate these groups often retain a significant number of droplets that do not contain cells or empty droplets. Additionally, there are currently no methods available to detect droplets containing damaged cells, which comprise partially lysed cells, the original source of the ambient RNA. Results Here, we describe DropletQC, a new method that is able to detect empty droplets, damaged, and intact cells, and accurately distinguish them from one another. This approach is based on a novel quality control metric, the nuclear fraction, which quantifies for each droplet the fraction of RNA originating from unspliced, nuclear pre-mRNA. We demonstrate how DropletQC provides a powerful extension to existing computational methods for identifying empty droplets such as EmptyDrops. Conclusions We implement DropletQC as an R package, which can be easily integrated into existing single-cell analysis workflows.

The tumor immune microenvironment is a main contributor to cancer progression and a promising therapeutic target for oncology. However, immune microenvironments vary profoundly between patients, and biomarkers for prognosis and treatment response lack precision. A comprehensive compendium of tumor immune cells is required to pinpoint predictive cellular states and their spatial localization. We generated a single-cell tumor immune atlas, jointly analyzing published data sets of >500,000 cells from 217 patients and 13 cancer types, providing the basis for a patient stratification based on immune cell compositions. Projecting immune cells from external tumors onto the atlas facilitated an automated cell annotation system. To enable in situ mapping of immune populations for digital pathology, we applied SPOTlight, combining single-cell and spatial transcriptomics data and identifying colocalization patterns of immune, stromal, and cancer cells in tumor sections. We expect the tumor immune cell atlas, together with our versatile toolbox for precision oncology, to advance currently applied stratification approaches for prognosis and immunotherapy. This publication is part of a project (BCLLATLAS) that has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 810287). This work has received funding from the Ministerio de Ciencia, Innovación y Universidades (SAF2017-89109-P; AEI/FEDER, UE) and the Fundació La Marató de TV3 (201903-30-31-32).We further acknowledge funding from the St. Vincent’s Clinic Foundation (V.T.C.) and the National Health and Medical Research Council Investigator Grant (APP1175781, J.E.P.), the Fundación Asociación Española contra el Cáncer (AECC), FERO (EDM), Ramón Areces Foundation, Cellex Foundation, BBVA (CAIMI), the ISCIII, FIS (PI16/01278), Juan de la Cierva formación fellowship (C.R.-P.) and Sara Borrell fellowship (E.P.-R.). Core funding is from the ISCIII and the Generalitat de Catalunya. We acknowledge support of the Spanish Ministry of Economy, Industry and Competitiveness (MEIC) to the EMBL partnership, the Centro de Excelencia Severo Ochoa, the CERCA Programme / Generalitat de Catalunya, the Spanish Ministry of Economy, Industry and Competitiveness (MEIC) through the nstituto de Salud Carlos III and the Generalitat de Catalunya through Departament de Salut and Departament d’Empresa i Coneixement. We also acknowledge the cofinancing by the Spanish Ministry of Economy, Industry and Competitiveness (MEIC) with funds from the European Regional Development Fund (ERDF) corresponding to the 2014–2020 Smart Growth Operating Program

Trait-associated genetic variants affect complex phenotypes primarily via regulatory mechanisms on the transcriptome. To investigate the genetics of gene expression, we performed cis- and trans-expression quantitative trait locus (eQTL) analyses using blood-derived expression from 31,684 individuals through the eQTLGen Consortium. We detected cis-eQTL for 88% of genes, and these were replicable in numerous tissues. Distal trans-eQTL (detected for 37% of 10,317 trait-associated variants tested) showed lower replication rates, partially due to low replication power and confounding by cell type composition. However, replication analyses in single-cell RNA-seq data prioritized intracellular trans-eQTL. Trans-eQTL exerted their effects via several mechanisms, primarily through regulation by transcription factors. Expression of 13% of the genes correlated with polygenic scores for 1,263 phenotypes, pinpointing potential drivers for those traits. In summary, this work represents a large eQTL resource, and its results serve as a starting point for in-depth interpretation of complex phenotypes.Analyses of expression profiles from whole blood of 31,684 individuals identify cis-expression quantitative trait loci (eQTL) effects for 88% of genes and trans-eQTL effects for 37% of trait-associated variants.

Induced pluripotent stem cells generated from patients with geographic atrophy as well as healthy individuals were differentiated to retinal pigment epithelium (RPE) cells. By integrating transcriptional profiles of 127,659 RPE cells generated from 43 individuals with geographic atrophy and 36 controls with genotype data, we identified 439 expression Quantitative Trait (eQTL) loci in cis that were associated with disease status and specific to subpopulations of RPE cells. We identified loci linked to two genes with known associations with geographic atrophy - PILRB and PRPH2, in addition to 43 genes with significant genotype x disease interactions that are candidates for novel genetic associations for geographic atrophy. On a transcriptome-only level, we identified molecular pathways significantly upregulated in geographic atrophy-RPE including in extracellular cellular matrix reorganisation, neurodegeneration, and mitochondrial functions. We subsequently implemented a large-scale proteomics analysis, confirming modification in proteins associated with these pathways. We also identified six significant protein (p) QTL that regulate protein expression in the RPE cells and in geographic atrophy - two of which share variants with cis-eQTL. Transcriptome-wide association analysis identified genes at loci previously associated with age-related macular degeneration. Further analysis conditional on disease status, implicated statistically significant RPE-specific eQTL. This study uncovers important differences in RPE homeostasis associated with geographic atrophy.

Epistasis has rarely been shown among natural polymorphisms in human traits; this research using advanced computation and gene expression data reveals many instances of epistasis between common single nucleotide polymorphisms in humans, with epistasis and the direction of its effect replicating in independent cohorts. Although frequently demonstrated in model organisms and domesticated species, very few examples of epistasis — where the effect of one polymorphism on a trait depends on other polymorphisms present in the genome — have been demonstrated in humans. Using advanced computation and a gene expression study design, these authors reveal hundreds of instances of epistasis between common single nucleotide polymorphisms (SNPs) in humans. They show that epistasis and direction of its effect replicate in independent cohorts. Epistatic networks of three or more SNPs are shown to influence the expression levels of many genes, whereby one cis-acting SNP is modulated by several trans-acting SNPs. Epistasis is the phenomenon whereby one polymorphism’s effect on a trait depends on other polymorphisms present in the genome. The extent to which epistasis influences complex traits1 and contributes to their variation2,3 is a fundamental question in evolution and human genetics. Although often demonstrated in artificial gene manipulation studies in model organisms4,5, and some examples have been reported in other species6, few examples exist for epistasis among natural polymorphisms in human traits7,8. Its absence from empirical findings may simply be due to low incidence in the genetic control of complex traits2,3, but an alternative view is that it has previously been too technically challenging to detect owing to statistical and computational issues9. Here we show, using advanced computation10 and a gene expression study design, that many instances of epistasis are found between common single nucleotide polymorphisms (SNPs). In a cohort of 846 individuals with 7,339 gene expression levels measured in peripheral blood, we found 501 significant pairwise interactions between common SNPs influencing the expression of 238 genes (P

Background High throughput single-cell RNA sequencing (scRNA-Seq) has emerged as a powerful tool for exploring cellular heterogeneity among complex human cancers. scRNA-Seq studies using fresh human surgical tissue are logistically difficult, preclude histopathological triage of samples, and limit the ability to perform batch processing. This hindrance can often introduce technical biases when integrating patient datasets and increase experimental costs. Although tissue preservation methods have been previously explored to address such issues, it is yet to be examined on complex human tissues, such as solid cancers and on high throughput scRNA-Seq platforms. Methods Using the Chromium 10X platform, we sequenced a total of ~ 120,000 cells from fresh and cryopreserved replicates across three primary breast cancers, two primary prostate cancers and a cutaneous melanoma. We performed detailed analyses between cells from each condition to assess the effects of cryopreservation on cellular heterogeneity, cell quality, clustering and the identification of gene ontologies. In addition, we performed single-cell immunophenotyping using CITE-Seq on a single breast cancer sample cryopreserved as solid tissue fragments. Results Tumour heterogeneity identified from fresh tissues was largely conserved in cryopreserved replicates. We show that sequencing of single cells prepared from cryopreserved tissue fragments or from cryopreserved cell suspensions is comparable to sequenced cells prepared from fresh tissue, with cryopreserved cell suspensions displaying higher correlations with fresh tissue in gene expression. We showed that cryopreservation had minimal impacts on the results of downstream analyses such as biological pathway enrichment. For some tumours, cryopreservation modestly increased cell stress signatures compared to freshly analysed tissue. Further, we demonstrate the advantage of cryopreserving whole-cells for detecting cell-surface proteins using CITE-Seq, which is impossible using other preservation methods such as single nuclei-sequencing. Conclusions We show that the viable cryopreservation of human cancers provides high-quality single-cells for multi-omics analysis. Our study guides new experimental designs for tissue biobanking for future clinical single-cell RNA sequencing studies.

East Coast fever, a tick-borne cattle disease caused by theTheileria parvaparasite, is among the biggest natural killers of cattle in East Africa, leading to over 1 million deaths annually. Here we report on the genetic analysis of a cohort of Boran cattle demonstrating heritable tolerance to infection byT. parva(h2= 0.65, s.e. 0.57). Through a linkage analysis we identify a 6 Mb genomic region onBos tauruschromosome 15 that is significantly associated with survival outcome followingT. parvaexposure. Testing this locus in an independent cohort of animals replicates this association with survival followingT. parvainfection. A stop gained polymorphism in this region was found to be highly associated with survival across both related and unrelated animals, with only one of the 20 homozygote carriers (T/T) of this change succumbing to the disease in contrast to 44 out of 97 animals homozygote for the reference allele (C/C). Consequently, we present a genetic locus linked to tolerance of one of Africa’s most important cattle diseases, raising the promise of marker-assisted selection for cattle that are less susceptible to infection byT. parva.Author SummaryMore than a million cattle die of East Coast fever in Africa each year, the impact of which disproportionately falls onto low-income, smallholder farmers. The lack of a widely accessible vaccine, heavy reliance on chemicals to control the tick vector and inadequate drug treatments means that new approaches for controlling the disease are urgently required. Through a genetic study of an extended pedigree of Boran cattle that are more than three times less likely to succumb to the disease than matched controls, we identify a region on chromosome 15 of the cattle genome associated with a high level of tolerance to the disease. We show that a variant in this region is also associated with survival in an independent cohort, and is linked to rates of cell expansion during infection. This genetic variant can therefore support marker-assisted selection, allowing farmers to breed tolerant cattle and offers a route to introduce this beneficial DNA to non-native breeds, enabling reduced disease incidence and increased productivity, which would be of benefit to millions of rural smallholder farmers across Africa.

Optical Coherence Tomography (OCT) enables non-invasive imaging of the retina and is used to diagnose and manage ophthalmic diseases including glaucoma. We present the first large-scale genome-wide association study of inner retinal morphology using phenotypes derived from OCT images of 31,434 UK Biobank participants. We identify 46 loci associated with thickness of the retinal nerve fibre layer or ganglion cell inner plexiform layer. Only one of these loci has been associated with glaucoma, and despite its clear role as a biomarker for the disease, Mendelian randomisation does not support inner retinal thickness being on the same genetic causal pathway as glaucoma. We extracted overall retinal thickness at the fovea, representative of foveal hypoplasia, with which three of the 46 SNPs were associated. We additionally associate these three loci with visual acuity. In contrast to the Mendelian causes of severe foveal hypoplasia, our results suggest a spectrum of foveal hypoplasia, in part genetically determined, with consequences on visual function., Author summary The thickness of the inner retinal layers is one of the biomarkers for glaucoma, the leading cause of irreversible blindness globally. Here we utilised the large-scale of the UK Biobank and the images of the retina it contains to look for genetic variants that effect the thickness of the inner retina. We find many variants associated with this variable, but surprisingly only one that also affects glaucoma. Further analysis shows that glaucoma and genetically determined inner retinal thickness are not on the same genetic pathway and it is rather the change of thickness over time that is indicative of the disease. This is important as it invites the potential for the discovered variants to be used as a representation of baseline thickness in the clinic in the future. We also show that foveal hypoplasia, the lack of the normal valley-like shape of the central retina, is present at a population level in a mild form and is affected by three variants that also affect visual acuity. This is an interesting discovery as foveal hypoplasia was previously thought of as an outcome of rare Mendelian disease.

Coronary artery disease remains the leading cause of death globally and is a major burden to every health system in the world. There have been significant improvements in risk modification, treatments, and mortality; however, our ability to detect asymptomatic disease for early intervention remains limited. Recent discoveries regarding the inflammatory nature of atherosclerosis have prompted investigation into new methods of diagnosis and treatment of coronary artery disease. This article reviews some of the highlights of the important developments in cardioimmunology and summarizes the clinical evidence linking the immune system and atherosclerosis. It provides an overview of the major serological biomarkers that have been associated with atherosclerosis, noting the limitations of these markers attributable to low specificity, and then contrasts these serological markers with the circulating immune cell subtypes that have been found to be altered in coronary artery disease. This review then outlines the technique of mass cytometry and its ability to provide high‐dimensional single‐cell data and explores how this high‐resolution quantification of specific immune cell subpopulations may assist in the diagnosis of early atherosclerosis in combination with other complimentary techniques such as single‐cell RNA sequencing. We propose that this improved specificity has the potential to transform the detection of coronary artery disease in its early phases, facilitating targeted preventative approaches in the precision medicine era.

Breast cancers are complex cellular ecosystems where heterotypic interactions play central roles in disease progression and response to therapy. However, our knowledge of their cellular composition and organization is limited. Here we present a single-cell and spatially resolved transcriptomics analysis of human breast cancers. We developed a single-cell method of intrinsic subtype classification (SCSubtype) to reveal recurrent neoplastic cell heterogeneity. Immunophenotyping using cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) provides high-resolution immune profiles, including new PD-L1/PD-L2+ macrophage populations associated with clinical outcome. Mesenchymal cells displayed diverse functions and cell-surface protein expression through differentiation within three major lineages. Stromal-immune niches were spatially organized in tumors, offering insights into antitumor immune regulation. Using single-cell signatures, we deconvoluted large breast cancer cohorts to stratify them into nine clusters, termed 'ecotypes', with unique cellular compositions and clinical outcomes. This study provides a comprehensive transcriptional atlas of the cellular architecture of breast cancer.

The tumour stroma regulates nearly all stages of carcinogenesis. Stromal heterogeneity in human triple-negative breast cancers (TNBCs) remains poorly understood, limiting the development of stromal-targeted therapies. Single-cell RNA sequencing of five TNBCs revealed two cancer-associated fibroblast (CAF) and two perivascular-like (PVL) subpopulations. CAFs clustered into two states: the first with features of myofibroblasts and the second characterised by high expression of growth factors and immunomodulatory molecules. PVL cells clustered into two states consistent with a differentiated and immature phenotype. We showed that these stromal states have distinct morphologies, spatial relationships and functional properties in regulating the extracellular matrix. Using cell signalling predictions, we provide evidence that stromal-immune crosstalk acts via a diverse array of immunoregulatory molecules. Importantly, the investigation of gene signatures from inflammatory-CAFs and differentiated-PVL cells in independent TNBC patient cohorts revealed strong associations with cytotoxic T-cell dysfunction and exclusion, respectively. Such insights present promising candidates to further investigate for new therapeutic strategies in the treatment of TNBCs.

BACKGROUND Chronic lung allograft dysfunction (CLAD) is the leading cause of mortality in lung transplant recipients. CLAD is characterized by respiratory failure owing to the accumulation of fibrotic cells in small airways and alveoli, inducing tissue contraction and architectural destruction. However, the source of the fibroblastic cells and the mechanism(s) underlying the accumulation and activation remain unexplained. Mesenchymal stromal cells (MSCs) are multipotent progenitors that are normally located in the lung tissue but can be isolated from the alveolar space in lung transplant recipients, where they have a profibrotic phenotype. Our objective was to identify the mediator(s) inducing migration and contractile differentiation of lung tissue MSCs. METHODS Bronchoalveolar lavage (BAL) (7 healthy controls and 21 lung transplant recipients), CCL2, HGF, TGFB, EGF, and PDGF-BB and autotaxin were measured by enzyme-linked immunosorbent assay. BAL (7 healthy controls and 31 lung transplant recipients) lysophosphatidic acid (LPA) (16:0, 18:0, 18:1, 22:4) was measured by liquid chromatography with tandem mass spectrometry. The effect of inhibition of candidate mediators on BAL-mediated chemoattraction of MSCs and contraction of MSC-spiked collagen gel assays was assessed. BAL cells from a lung transplant recipient with CLAD were analyzed by single-cell RNA sequencing. RESULTS We first demonstrate that BAL fluid from lung transplant recipients and particularly those with CLAD is potently chemoattractive to human lung tissue‒derived MSCs and induces a contractile phenotype. After excluding several candidate mediators, we show that LPA blockade completely abrogated transplant recipient BAL‒mediated chemoattraction of MSCs and contraction of MSC-spiked collagen gels. Furthermore, LPA levels were enriched in transplant recipient BAL, and LPA replicated the observed in vitro profibrotic effects of transplant recipient BAL. Finally, we identify BAL monocyte‒derived macrophages with autotaxin (ENPP2) and fibrotic transcriptional signature. CONCLUSIONS Autotaxin-expressing alveolar macrophages are present in CLAD BAL. These cells potentially provide a local source of autotaxin/LPA that drives MSC recruitment and tissue contraction in CLAD. These cells are analogous to an aberrant macrophage population recently identified in idiopathic pulmonary fibrosis, suggesting an overlap in pathogenesis between CLAD and other forms of lung fibrosis.

The discovery that somatic cells can be reprogrammed to induced pluripotent stem cells (iPSCs) - cells that can be differentiated into any cell type of the three germ layers - has provided a foundation forin vitrohuman disease modelling1,2, drug development1,2, and population genetics studies3,4. In the majority of instances, the expression levels of genes, plays a critical role in contributing to disease risk, or the ability to identify therapeutic targets. However, while the effect of the genetic background of cell lines has been shown to strongly influence gene expression, the effect has not been evaluated at the level of individual cells. Differences in the effect of genetic variation on the gene expression of different cell-types, would provide significant resolution for in vitro research using preprogramed cells. By bringing together single cell RNA sequencing15–21and population genetics, we now have a framework in which to evaluate the cell-types specific effects of genetic variation on gene expression. Here, we performed single cell RNA-sequencing on 64,018 fibroblasts from 79 donors and we mapped expression quantitative trait loci (eQTL) at the level of individual cell types. We demonstrate that the large majority of eQTL detected in fibroblasts are specific to an individual sub-type of cells. To address if the allelic effects on gene expression are dynamic across cell reprogramming, we generated scRNA-seq data in 19,967 iPSCs from 31 reprogramed donor lines. We again identify highly cell type specific eQTL in iPSCs, and show that that the eQTL in fibroblasts are almost entirely disappear during reprogramming. This work provides an atlas of how genetic variation influences gene expression across cell subtypes, and provided evidence for patterns of genetic architecture that lead to cell-types specific eQTL effects.

Idiopathic pulmonary fibrosis (IPF) is a fatal lung disease in which airway macrophages (AMs) play a key role. Itaconate has emerged as a mediator of macrophage function, but its role during fibrosis is unknown. Here, we reveal that itaconate is an endogenous antifibrotic factor in the lung. Itaconate levels are reduced in bronchoalveolar lavage, and itaconate-synthesizing cis-aconitate decarboxylase expression (ACOD1) is reduced in AMs from patients with IPF compared with controls. In the murine bleomycin model of pulmonary fibrosis, Acod1-/- mice develop persistent fibrosis, unlike wild-type (WT) littermates. Profibrotic gene expression is increased in Acod1-/- tissue-resident AMs compared with WT, and adoptive transfer of WT monocyte-recruited AMs rescued mice from disease phenotype. Culture of lung fibroblasts with itaconate decreased proliferation and wound healing capacity, and inhaled itaconate was protective in mice in vivo. Collectively, these data identify itaconate as critical for controlling the severity of lung fibrosis, and targeting this pathway may be a viable therapeutic strategy.

In recent years, functional genomics approaches combining genetic information with bulk RNA-sequencing data have identified the downstream expression effects of disease-associated genetic risk factors through so-called expression quantitative trait locus (eQTL) analysis. Single-cell RNA-sequencing creates enormous opportunities for mapping eQTLs across different cell types and in dynamic processes, many of which are obscured when using bulk methods. Rapid increase in throughput and reduction in cost per cell now allow this technology to be applied to large-scale population genetics studies. To fully leverage these emerging data resources, we have founded the single-cell eQTLGen consortium (sc-eQTLGen), aimed at pinpointing the cellular contexts in which disease-causing genetic variants affect gene expression. Here, we outline the goals, approach and potential utility of the sc-eQTLGen consortium. We also provide a set of study design considerations for future single-cell eQTL studies.

The past four decades have yielded advances in molecular biology allowing detailed characterization of the cellular genome and the transcriptome: the complete set of RNA species transcribed by a cell or tissue. Through transcriptomics and next-generation sequencing, we can now attain an unprecedented level of detail in understanding cellular phenotypes through examining the genes expressed in specific physiological and pathological states. In this review, we provide an overview of transcriptomics and RNA-sequencing in the analysis of whole tissue and single cells. We describe the techniques and pitfalls involved in the isolation and sequencing of single cells, and what additional benefits this application can provide. Finally, we look to how these technologies are being applied in pulmonary research, and how they may translate in the near future into clinical practice.

A variety of methods have been developed to demultiplex pooled samples in a single cell RNA sequencing (scRNA-seq) experiment which either require hashtag barcodes or sample genotypes prior to pooling. We introduce scSplit which utilizes genetic differences inferred from scRNA-seq data alone to demultiplex pooled samples. scSplit also enables mapping clusters to original samples. Using simulated, merged, and pooled multi-individual datasets, we show that scSplit prediction is highly concordant with demuxlet predictions and is highly consistent with the known truth in cell-hashing dataset. scSplit is ideally suited to samples without external genotype information and is available at: https://github.com/jon-xu/scSplit

We characterized monocyte/macrophage populations from infants to older adults and lineage traced these cells using sex-mismatched lung transplant recipients. Classical monocytes peak in adulthood and decline in older adults. Importantly, our data show that after transplant, donor human airway macrophages (AMs) are readily replaced by recipient-derived monocytes., The ontogeny of airway macrophages (AMs) in human lung and their contribution to disease are poorly mapped out. In mice, aging is associated with an increasing proportion of peripherally, as opposed to perinatally derived AMs. We sought to understand AM ontogeny in human lung during healthy aging and after transplant. We characterized monocyte/macrophage populations from the peripheral blood and airways of healthy volunteers across infancy/childhood (2–12 yr), maturity (20–50 yr), and older adulthood (>50 yr). Single-cell RNA sequencing (scRNA-seq) was performed on airway inflammatory cells isolated from sex-mismatched lung transplant recipients. During healthy aging, the proportions of blood and bronchoalveolar lavage (BAL) classical monocytes peak in adulthood and decline in older adults. scRNA-seq of BAL cells from lung transplant recipients indicates that after transplant, the majority of AMs are recipient derived. These data show that during aging, the peripheral monocyte phenotype is consistent with that found in the airways and, furthermore, that the majority of human AMs after transplant are derived from circulating monocytes.

Trans-eQTLs have been implicated in complex traits and common diseases, but many were initially identified on the basis of having an effect in cis, and there has been no assessment of the significance of the overlap in relation to chance expectations. Here, we investigated whether trans-expression quantitative trait loci (eQTL) associations identified in whole blood contribute to variance in complex traits by determining (1) whether genome-wide significant (GWS) single-nucleotide polymorphisms (SNPs) were enriched for trans-eQTL (including trans-only eQTL), and (2) whether the genomic regions surrounding associated trans-genes were enriched for statistical associations in the relevant GWAS. On average for a given phenotype, we identify 4.8% of GWS SNPs overlapping with trans-eQTL present in blood, and show that for the majority of these phenotypes, this observation does not exceed that expected by chance. Likewise, we observe no enrichment for genetic associations with the GWAS phenotype in the regions surrounding the linked trans-genes, with the exception of rheumatoid arthritis. Interestingly, the GWS SNPs for each phenotype were consistently more enriched for unique trans-eQTL SNPs than trans-eQTL SNP-probe pairs (p = 4 × 10−7), with schizophrenia the only exception. This relative enrichment for trans-eQTL SNPs over trans-eQTL SNP-probe pairs implies that trait-associated trans-eQTL SNPs in whole blood are less likely to be 'master regulators' than random trans-eQTL SNPs. Taken together, these results suggest little evidence for the role of blood-based trans-eQTL in complex traits and disease, although this may reflect the finite size of currently available data sets and our findings may not hold for trans-eQTLs in more trait-relevant tissues. All software is publically available at https://github.com/IMB-Computational-Genomics-Lab/eqtlOverlapper .

Heterogeneity of cell states represented in pluripotent cultures have not been described at the transcriptional level. Since gene expression is highly heterogeneous between cells, single-cell RNA sequencing can be used to identify how individual pluripotent cells function. Here, we present results from the analysis of single-cell RNA sequencing data from 18,787 individual WTC CRISPRi human induced pluripotent stem cells. We developed an unsupervised clustering method, and through this identified four subpopulations distinguishable on the basis of their pluripotent state including: a core pluripotent population (48.3%), proliferative (47.8%), early-primed for differentiation (2.8%) and late-primed for differentiation (1.1%). For each subpopulation we were able to identify the genes and pathways that define differences in pluripotent cell states. Our method identified four transcriptionally distinct predictor gene sets comprised of 165 unique genes that denote the specific pluripotency states; and using these sets, we developed a multigenic machine learning prediction method to accurately classify single cells into each of the subpopulations. Compared against a set of established pluripotency markers, our method increases prediction accuracy by 10%, specificity by 20%, and explains a substantially larger proportion of deviance (up to 3-fold) from the prediction model. Finally, we developed an innovative method to predict cells transitioning between subpopulations, and support our conclusions with results from two orthogonal pseudotime trajectory methods.

Summary: DNA replication timing and three-dimensional (3D) genome organization are associated with distinct epigenome patterns across large domains. However, whether alterations in the epigenome, in particular cancer-related DNA hypomethylation, affects higher-order levels of genome architecture is still unclear. Here, using Repli-Seq, single-cell Repli-Seq, and Hi-C, we show that genome-wide methylation loss is associated with both concordant loss of replication timing precision and deregulation of 3D genome organization. Notably, we find distinct disruption in 3D genome compartmentalization, striking gains in cell-to-cell replication timing heterogeneity and loss of allelic replication timing in cancer hypomethylation models, potentially through the gene deregulation of DNA replication and genome organization pathways. Finally, we identify ectopic H3K4me3-H3K9me3 domains from across large hypomethylated domains, where late replication is maintained, which we purport serves to protect against catastrophic genome reorganization and aberrant gene transcription. Our results highlight a potential role for the methylome in the maintenance of 3D genome regulation.

The evolutionary theories of mutation accumulation (MA) and disposable soma (DS) provide possible explanations for the existence of human aging. To better understand the relative importance of these theories, we devised a test to identify MA- and DS-consistent sites across the genome using familial DNA methylation data. Two key characteristics of DNA methylation allowed us to do so. First, DNA methylation exhibits distinct and widespread changes with age, with numerous age-differentially-methylated sites observed across the genome. Second, many sites show heritable DNA methylation patterns within families. We extended heritability predictions of MA and DS to DNA methylation, predicting that MA-consistent age-differentially-methylated sites will show increasing heritability with age, while DS-consistent sites will show the opposite. Variance components models were used to test for changing heritability of methylation with age at 48,601 age-differentially-methylated sites across the genome in 610 individuals from 176 families. Of these, 102 sites showed significant MA-consistent increases in heritability with age, while 2266 showed significant DS-consistent decreases in heritability. These results suggest that both MA and DS play a role in explaining aging and aging-related changes, and that while the majority of DNA methylation changes observed in aging are consistent with epigenetic drift, targeted changes exist and may mediate effects of aging-related genes.

Expression quantitative trait locus (eQTL) detection has emerged as an important tool for unraveling of the relationship between genetic risk factors and disease or clinical phenotypes. Most studies use single marker linear regression to discover primary signals, followed by sequential conditional modeling to detect secondary genetic variants affecting gene expression. However, this approach assumes that functional variants are sparsely distributed and that close linkage between them has little impact on estimation of their precise location and the magnitude of effects. We describe a series of simulation studies designed to evaluate the impact of linkage disequilibrium (LD) on the fine mapping of causal variants with typical eQTL effect sizes. In the presence of multisite regulation, even though between 80 and 90% of modeled eSNPs associate with normally distributed traits, up to 10% of all secondary signals could be statistical artifacts, and at least 5% but up to one-quarter of credible intervals of SNPs within r2 > 0.8 of the peak may not even include a causal site. The Bayesian methods eCAVIAR and DAP (Deterministic Approximation of Posteriors) provide only modest improvement in resolution. Given the strong empirical evidence that gene expression is commonly regulated by more than one variant, we conclude that the fine mapping of causal variants needs to be adjusted for multisite influences, as conditional estimates can be highly biased by interference among linked sites, but ultimately experimental verification of individual effects is needed. Presumably similar conclusions apply not just to eQTL mapping, but to multisite influences on fine mapping of most types of quantitative trait.

Background Typical migraine is a frequent, debilitating and painful headache disorder with an estimated heritability of about 50%. Although genome-wide association (GWA) studies have identified over 40 single nucleotide polymorphisms associated with migraine, further research is required to determine their biological role in migraine susceptibility. Therefore, we performed a study of genome-wide gene expression in a large sample of 83 migraine cases and 83 non-migraine controls to determine whether altered expression levels of genes and pathways could provide insights into the biological mechanisms underlying migraine. Methods We assessed whole blood gene expression data for 17994 expression probes measured using IlluminaHT-12 v4.0 BeadChips. Differential expression was assessed using multivariable logistic regression. Gene expression probes with a nominal p value Results Association analyses between migraine and probe expression levels, adjusted for age and gender, revealed an excess of small p values, but there was no significant single-probe association after correction for multiple testing. Network analysis of pooled expression data identified 10 modules of co-expressed genes. One module harboured a significant number of differentially expressed genes and was strongly enriched with immune-inflammatory pathways, including multiple pathways expressed in microglial cells. Conclusions These data suggest immune-inflammatory pathways play an important role in the pathogenesis, manifestation, and/or progression of migraine in some patients. Furthermore, gene-expression associations are measurable in whole blood, suggesting the analysis of blood gene expression can inform our understanding of the biological mechanisms underlying migraine, identify biomarkers, and facilitate the discovery of novel pathways and thus determine new targets for drug therapy.

We show that genotype-by-environment interaction can be inferred from an analysis without environmental data in a large sample., Genotype-by-environment interaction (GEI) is a fundamental component in understanding complex trait variation. However, it remains challenging to identify genetic variants with GEI effects in humans largely because of the small effect sizes and the difficulty of monitoring environmental fluctuations. Here, we demonstrate that GEI can be inferred from genetic variants associated with phenotypic variability in a large sample without the need of measuring environmental factors. We performed a genome-wide variance quantitative trait locus (vQTL) analysis of ~5.6 million variants on 348,501 unrelated individuals of European ancestry for 13 quantitative traits in the UK Biobank and identified 75 significant vQTLs with P < 2.0 × 10−9 for 9 traits, especially for those related to obesity. Direct GEI analysis with five environmental factors showed that the vQTLs were strongly enriched with GEI effects. Our results indicate pervasive GEI effects for obesity-related traits and demonstrate the detection of GEI without environmental data.

Expression QTL (eQTL) detection has emerged as an important tool for unraveling the relationship between genetic risk factors and disease or clinical phenotypes. Most studies are predicated on the assumption that only a single causal variant explains the association signal in each interval. This greatly simplifies the statistical modeling, but is liable to biases in scenarios where multiple local causal-variants are responsible. Here, our primary goal was to address the prevalence of secondary cis-eQTL signals regulating peripheral blood gene expression locally, utilizing two large human cohort studies, each >2500 samples with accompanying whole genome genotypes. The CAGE (Consortium for the Architecture of Gene Expression) dataset is a compendium of Illumina microarray studies, and the Framingham Heart Study is a two-generation Affymetrix dataset. We also describe Bayesian colocalization analysis of the extent of sharing of cis-eQTL detected in both studies as well as with the BIOS RNAseq dataset. Stepwise conditional modeling demonstrates that multiple eQTL signals are present for ∼40% of over 3500 eGenes in both microarray datasets, and that the number of loci with additional signals reduces by approximately two-thirds with each conditioning step. Although

Many computational methods have been developed to infer cell type proportions from bulk transcriptomics data. However, an evaluation of the impact of data transformation, pre-processing, marker selection, cell type composition and choice of methodology on the deconvolution results is still lacking. Using five single-cell RNA-sequencing (scRNA-seq) datasets, we generate pseudo-bulk mixtures to evaluate the combined impact of these factors. Both bulk deconvolution methodologies and those that use scRNA-seq data as reference perform best when applied to data in linear scale and the choice of normalization has a dramatic impact on some, but not all methods. Overall, methods that use scRNA-seq data have comparable performance to the best performing bulk methods whereas semi-supervised approaches show higher error values. Moreover, failure to include cell types in the reference that are present in a mixture leads to substantially worse results, regardless of the previous choices. Altogether, we evaluate the combined impact of factors affecting the deconvolution task across different datasets and propose general guidelines to maximize its performance.

Genome-wide association studies (GWAS) have identified thousands of genetic variants associated with human complex traits. However, the genes or functional DNA elements through which these variants exert their effects on the traits are often unknown. We propose a method (called SMR) that integrates summary-level data from GWAS with data from expression quantitative trait locus (eQTL) studies to identify genes whose expression levels are associated with a complex trait because of pleiotropy. We apply the method to five human complex traits using GWAS data on up to 339,224 individuals and eQTL data on 5,311 individuals, and we prioritize 126 genes (for example, TRAF1 and ANKRD55 for rheumatoid arthritis and SNX19 and NMRAL1 for schizophrenia), of which 25 genes are new candidates; 77 genes are not the nearest annotated gene to the top associated GWAS SNP. These genes provide important leads to design future functional studies to understand the mechanism whereby DNA variation leads to complex trait variation.

Summary There is pressing urgency to understand the pathogenesis of the severe acute respiratory syndrome coronavirus clade 2 (SARS-CoV-2), which causes the disease COVID-19. SARS-CoV-2 spike (S) protein binds angiotensin-converting enzyme 2 (ACE2), and in concert with host proteases, principally transmembrane serine protease 2 (TMPRSS2), promotes cellular entry. The cell subsets targeted by SARS-CoV-2 in host tissues and the factors that regulate ACE2 expression remain unknown. Here, we leverage human, non-human primate, and mouse single-cell RNA-sequencing (scRNA-seq) datasets across health and disease to uncover putative targets of SARS-CoV-2 among tissue-resident cell subsets. We identify ACE2 and TMPRSS2 co-expressing cells within lung type II pneumocytes, ileal absorptive enterocytes, and nasal goblet secretory cells. Strikingly, we discovered that ACE2 is a human interferon-stimulated gene (ISG) in vitro using airway epithelial cells and extend our findings to in vivo viral infections. Our data suggest that SARS-CoV-2 could exploit species-specific interferon-driven upregulation of ACE2, a tissue-protective mediator during lung injury, to enhance infection., Graphical Abstract, Highlights • Meta-analysis of human, non-human primate, and mouse single-cell RNA-seq datasets for putative SARS-CoV-2 targets • Type II pneumocytes, nasal secretory cells, and absorptive enterocytes are ACE2 + TMPRSS2 + • Interferon and influenza increase ACE2 in human nasal epithelia and lung tissue • Mouse Ace2 is not upregulated by interferon, raising implications for disease modeling, Analysis of single-cell RNA-seq datasets from human, non-human primate, and mouse barrier tissues identifies putative cellular targets of SARS-CoV-2 on the basis of ACE2 and TMPRSS2 expression. ACE2 represents a previously unappreciated interferon-stimulated gene in human, but not mouse, epithelial tissues, identifying anti-viral induction of a host tissue-protective mechanism, but also a potential means for viral exploitation of the host response.

The retina is a specialized neural tissue that senses light and initiates image processing. Although the functional organization of specific retina cells has been well studied, the molecular profile of many cell types remains unclear in humans. To comprehensively profile the human retina, we performed single-cell RNA sequencing on 20,009 cells from three donors and compiled a reference transcriptome atlas. Using unsupervised clustering analysis, we identified 18 transcriptionally distinct cell populations representing all known neural retinal cells: rod photoreceptors, cone photoreceptors, Muller glia, bipolar cells, amacrine cells, retinal ganglion cells, horizontal cells, astrocytes, and microglia. Our data captured molecular profiles for healthy and putative early degenerating rod photoreceptors, and revealed the loss of MALAT1 expression with longer post-mortem time, which potentially suggested a novel role of MALAT1 in rod photoreceptor degeneration. We have demonstrated the use of this retina transcriptome atlas to benchmark pluripotent stem cell-derived cone photoreceptors and an adult Muller glia cell line. This work provides an important reference with unprecedented insights into the transcriptional landscape of human retinal cells, which is fundamental to understanding retinal biology and disease.

Gene expression varies markedly across the menstrual cycle and expression levels for many genes are under genetic control. We analyzed gene expression and mapped expression quantitative trait loci (eQTLs) in endometrial tissue samples from 229 women and then analyzed the overlap of endometrial eQTL signals with genomic regions associated with endometriosis and other reproductive traits. We observed a total of 45,923 cis-eQTLs for 417 unique genes and 2,968 trans-eQTLs affecting 82 unique genes. Two eQTLs were located in known risk regions for endometriosis including LINC00339 on chromosome 1 and VEZT on chromosome 12 and there was evidence for eQTLs that may be target genes in genomic regions associated with other reproductive diseases. Dynamic changes in expression of individual genes across cycle include alterations in both mean expression and transcriptional silencing. Significant effects of cycle stage on mean expression levels were observed for (2,427/15,262) probes with detectable expression in at least 90% of samples and for (2,877/9,626) probes expressed in some, but not all samples. Pathway analysis supports similar biological control of both altered expression levels and transcriptional silencing. Taken together, these data identify strong genetic effects on genes with diverse functions in human endometrium and provide a platform for better understanding genetic effects on endometrial-related pathologies.

Chronic liver disease (CLD) is associated with tissue-destructive fibrosis. Considering that common mechanisms drive fibrosis across etiologies, and that steatosis is an important cofactor for pathology, we performed RNA sequencing on liver biopsies of patients with different fibrosis stages, resulting from infection with hepatitis C virus (HCV) (with or without steatosis) or fatty liver disease. In combination with enhanced liver fibrosis score correlation analysis, we reveal a common set of genes associated with advanced fibrosis, as exemplified by those encoding the transcription factor ETS-homologous factor (EHF) and the extracellular matrix protein versican (VCAN). We identified 17 fibrosis-associated genes as candidate EHF targets and demonstrated that EHF regulates multiple fibrosis-associated genes, including VCAN, in hepatic stellate cells. Serum VCAN levels were also elevated in advanced fibrosis patients. Comparing biopsies from patients with HCV with or without steatosis, we identified a steatosis-enriched gene set associated with advanced fibrosis, validating follistatin-like protein 1 (FSTL1) as an exemplar of this profile. In patients with advanced fibrosis, serum FSTL1 levels were elevated in those with steatosis (versus those without). Liver Fstl1 mRNA levels were also elevated in murine CLD models. We thus reveal a common gene signature for CLD-associated liver fibrosis and potential biomarkers and/or targets for steatosis-associated liver fibrosis.

Persistent human papillomavirus (HPV) infection is responsible for at least 5% of human malignancies. Most HPV-associated cancers are initiated by the HPV16 genotype, as confirmed by detection of integrated HPV DNA in cells of oral and anogenital epithelial cancers. However, single-cell RNA sequencing may enable prediction of HPV involvement in carcinogenesis at other sites. We conducted single-cell RNA sequencing on keratinocytes from a mouse transgenic for the E7 gene of HPV16 and showed sensitive and specific detection of HPV16-E7 mRNA, predominantly in basal keratinocytes. We showed that increased E7 mRNA copy number per cell was associated with increased expression of E7 induced genes. This technique enhances detection of active viral transcription in solid tissue and may clarify possible linkage of HPV infection to development of squamous cell carcinoma.

BackgroundGiant cell arteritis (GCA) is the most common form of vasculitis affecting elderly people. It is one of the few true ophthalmic emergencies. GCA is a heterogenous disease, symptoms and signs are variable thereby making it challenging to diagnose and often delaying diagnosis. A temporal artery biopsy is the gold standard to test for GCA, and there are currently no specific biochemical markers to categorize or aid diagnosis of the disease. We aimed to identify a less invasive method to confirm the diagnosis of GCA, as well as to ascertain clinically relevant predictive biomarkers by studying the transcriptome of purified peripheral CD4+ and CD8+ T lymphocytes in patients with GCA.Methods and FindingsWe recruited 16 patients with histological evidence of GCA at the Royal Victorian Eye and Ear Hospital (RVEEH), Melbourne, Australia, and aimed to collect blood samples at six time points: acute phase, 2–3 weeks, 6–8 weeks, 3 months, 6 months and 12 months after clinical diagnosis. CD4+ and CD8+ T-cells were positively selected at each time point through magnetic-assisted cell sorting (MACS). RNA was extracted from all 195 collected samples for subsequent RNA sequencing. The expression profiles of patients were compared to those of 16 age-matched controls. Over the 12-month study period, polynomial modelling analyses identified 179 and 4 statistically significant transcripts with altered expression profiles (FDR < 0.05) between cases and controls in CD4+ and CD8+ populations, respectively. In CD8+ cells, we identified two transcripts that remained differentially expressed after 12 months, namely SGTB, associated with neuronal apoptosis, and FCGR3A, which has been found in association with Takayasu arteritis (TA), another large vessel vasculitis. We detected genes that correlate with both symptoms and biochemical markers used in the acute setting for predicting long-term prognosis. 15 genes were shared across 3 phenotypes in CD4 and 16 across CD8 cells. In CD8, IL32 was common to 5 phenotypes: a history of Polymyalgia Rheumatica, both visual disturbance and raised neutrophils at the time of presentation, bilateral blindness and death within 12 months. Altered IL32 gene expression could provide risk evaluation of GCA diagnosis at the time of presentation and give an indication of prognosis, which may influence management.ConclusionsThis is the first longitudinal gene expression study undertaken to identify robust transcriptomic biomarkers of GCA. Our results show cell type-specific transcript expression profiles, novel gene-phenotype associations, and uncover important biological pathways for this disease. These data significantly enhance the current knowledge of relevant biomarkers, their association with clinical prognostic markers, as well as potential candidates for detecting disease activity in whole blood samples. In the acute phase, the gene-phenotype relationships we have identified could provide insight to potential disease severity and as such guide us in initiating appropriate patient management.

Transcript co-expression is regulated by a combination of shared genetic and environmental factors. Here, we estimate the proportion of co-expression that is due to shared genetic variance. To do so, we estimated the genetic correlations between each pairwise combination of 2469 transcripts that are highly heritable and expressed in whole blood in 1748 unrelated individuals of European ancestry. We identify 556 pairs with a significant genetic correlation of which 77% are located on different chromosomes, and report 934 expression quantitative trait loci, identified in an independent cohort, with significant effects on both transcripts in a genetically correlated pair. We show significant enrichment for transcription factor control and physical proximity through chromatin interactions as possible mechanisms of shared genetic control. Finally, we construct networks of interconnected transcripts and identify their underlying biological functions. Using genetic correlations to investigate transcriptional co-regulation provides valuable insight into the nature of the underlying genetic architecture of gene regulation., Covariance of gene expression pairs is due to a combination of shared genetic and environmental factors. Here the authors estimate the genetic correlation between highly heritable pairs and identify transcription factor control and chromatin interactions as possible mechanisms of correlation.

Cardiac differentiation of human pluripotent stem cells (hPSCs) requires orchestration of dynamic gene regulatory networks during stepwise fate transitions, but often generates immature cell types that do not fully recapitulate properties of their adult counterparts, suggesting incomplete activation of key transcriptional networks. We performed extensive single-cell transcriptomic analyses to map fate choices and gene expression programs during cardiac differentiation of hPSCs, and identified strategies to improve in vitro cardiomyocyte differentiation. Utilizing genetic gain- and loss-of-function approaches, we found that hypertrophic signaling is not effectively activated during monolayer-based cardiac differentiation, thereby preventing expression of HOPX and its activation of downstream genes that govern late stages of cardiomyocyte maturation. This study therefore provides a key transcriptional roadmap of in vitro cardiac differentiation at single-cell resolution, revealing fundamental mechanisms underlying heart development and differentiation of hPSC-derived cardiomyocytes.

We used single cell sequencing technology to characterize the transcriptomes of 1,174 human embryonic stem cell-derived retinal ganglion cells (RGCs) at the single cell level. The human embryonic stem cell line BRN3B-mCherry (A81-H7), was differentiated to RGCs using a guided differentiation approach. Cells were harvested at day 36 and prepared for single cell RNA sequencing. Our data indicates the presence of three distinct subpopulations of cells, with various degrees of maturity. One cluster of 288 cells showed increased expression of genes involved in axon guidance together with semaphorin interactions, cell-extracellular matrix interactions and ECM proteoglycans, suggestive of a more mature RGC phenotype.

We develop a Bayesian mixed linear model that simultaneously estimates single-nucleotide polymorphism (SNP)-based heritability, polygenicity (proportion of SNPs with nonzero effects), and the relationship between SNP effect size and minor allele frequency for complex traits in conventionally unrelated individuals using genome-wide SNP data. We apply the method to 28 complex traits in the UK Biobank data (N = 126,752) and show that on average, 6% of SNPs have nonzero effects, which in total explain 22% of phenotypic variance. We detect significant (P 0.05/28) signatures of natural selection in the genetic architecture of 23 traits, including reproductive, cardiovascular, and anthropometric traits, as well as educational attainment. The significant estimates of the relationship between effect size and minor allele frequency in complex traits are consistent with a model of negative (or purifying) selection, as confirmed by forward simulation. We conclude that negative selection acts pervasively on the genetic variants associated with human complex traits.

Do genetic effects regulate gene expression in human endometrium?This study demonstrated strong genetic effects on endometrial gene expression and some evidence for genetic regulation of gene expression in a menstrual cycle stage-specific manner.Genetic effects on expression levels for many genes are tissue specific. Endometrial gene expression varies across menstrual cycle stages and between individuals, but there are limited data on genetic control of expression in endometrium.We analysed genome-wide genotype and gene expression data to map cis expression quantitative trait loci (eQTL) in endometrium.We recruited 123 women of European ancestry. DNA samples from blood were genotyped on Illumina HumanCoreExome chips. Total RNA was extracted from endometrial tissues. Whole-transcriptome profiles were characterized using Illumina Human HT-12 v4.0 Expression Beadchips. We performed eQTL mapping with ~8 000 000 genotyped and imputed single nucleotide polymorphisms (SNPs) and 12 329 genes.We identified a total of 18 595 cis SNP-probe associations at a study-wide level of significance (P1 × 10-7), which correspond to independent eQTLs for 198 unique genes. The eQTLs with the largest effect in endometrial tissue were rs4902335 for CHURC1 (P = 1.05 × 10-32) and rs147253019 for ZP3 (P = 8.22 × 10-30). We further performed a context-specific eQTL analysis to investigate if genetic effects on gene expression regulation act in a menstrual cycle-specific manner. Interestingly, five cis-eQTLs were identified with a significant stage-by-genotype interaction. The strongest stage interaction was the eQTL for C10ORF33 (PYROXD2) with SNP rs2296438 (P = 2.0 × 10-4), where we observe a 2-fold difference in the average expression levels of heterozygous samples depending on the stage of the menstrual cycle.The summary eQTL results are publicly available to browse or download.A limitation of the present study was the relatively modest sample size. It was not powered to identify trans-eQTLs and larger sample sizes will also be needed to provide better power to detect cis-eQTLs and cycle stage-specific effects, given the substantial changes in expression across the menstrual cycle for many genes.Identification of endometrial eQTLs provides a platform for better understanding genetic effects on endometriosis risk and other endometrial-related pathologies.Funding for this work was provided by NHMRC Project Grants GNT1026033, GNT1049472, GNT1046880, GNT1050208, GNT1105321 and APP1083405. There are no competing interests.

The mitochondrial and nuclear genomes coordinate and co-evolve in eukaryotes in order to adapt to environmental changes. Variation in the mitochondrial genome is capable of affecting expression of genes on the nuclear genome. Sex-specific mitochondrial genetic control of gene expression has been demonstrated in Drosophila melanogaster, where males were found to drive most of the total variation in gene expression. This has potential implications for male-related health and disease resulting from variation in mtDNA solely inherited from the mother. We used a family-based study comprised of 47,323 gene expression probes and 78 mitochondrial SNPs (mtSNPs) from n = 846 individuals to examine the extent of mitochondrial genetic control of gene expression in humans. This identified 15 significant probe-mtSNP associations (P

Background: Multiple Sclerosis (MS) is an immune-mediated disease of the central nervous system which responds to therapies targeting circulating immune cells. Objective: Our aim was to test if the T-cell activation gene expression pattern (TCAGE) we had previously described from whole blood was replicated in an independent cohort. Methods: We used RNA-seq to interrogate the whole blood transcriptomes of 72 individuals (40 healthy controls, 32 untreated MS). A cohort of 862 control individuals from the Brisbane Systems Genetics Study (BSGS) was used to assess heritability and seasonal expression. The effect of interferon beta (IFNB) therapy on expression was evaluated. Results: The MS/TCAGE association was replicated and rationalized to a single marker, ribosomal protein S6 (RPS6). Expression of RPS6 was higher in MS than controls ( pConclusions: These data support investigation of RPS6 as a potential therapeutic target and candidate biomarker for measuring clinical response to IFNB and other MS therapies, and of MS disease heterogeneity.

There is evidence across several species for genetic control of phenotypic variation of complex traits, such that the variance among phenotypes is genotype dependent. Understanding genetic control of variability is important in evolutionary biology, agricultural selection programmes and human medicine, yet for complex traits, no individual genetic variants associated with variance, as opposed to the mean, have been identified. Here we perform a meta-analysis of genome-wide association studies of phenotypic variation using ∼170,000 samples on height and body mass index (BMI) in human populations. We report evidence that the single nucleotide polymorphism (SNP) rs7202116 at the FTO gene locus, which is known to be associated with obesity (as measured by mean BMI for each rs7202116 genotype), is also associated with phenotypic variability. We show that the results are not due to scale effects or other artefacts, and find no other experiment-wise significant evidence for effects on variability, either at loci other than FTO for BMI or at any locus for height. The difference in variance for BMI among individuals with opposite homozygous genotypes at the FTO locus is approximately 7%, corresponding to a difference of ∼0.5 kilograms in the standard deviation of weight. Our results indicate that genetic variants can be discovered that are associated with variability, and that between-person variability in obesity can partly be explained by the genotype at the FTO locus. The results are consistent with reported FTO by environment interactions for BMI, possibly mediated by DNA methylation. Our BMI results for other SNPs and our height results for all SNPs suggest that most genetic variants, including those that influence mean height or mean BMI, are not associated with phenotypic variance, or that their effects on variability are too small to detect even with samples sizes greater than 100,000. ispartof: Nature vol:490 issue:7419 pages:267-272 ispartof: location:England status: published

Genome-wide association studies (GWAS) have identified markers within the WNT4 region on chromosome 1p36.12 showing consistent and strong association with increasing endometriosis risk. Fine mapping using sequence and imputed genotype data has revealed strong candidates for the causal SNPs within these critical regions; however, the molecular pathogenesis of these SNPs is currently unknown. We used gene expression data collected from whole blood from 862 individuals and endometrial tissue from 136 individuals from independent populations of European descent to examine the mechanism underlying endometriosis susceptibility. Association mapping results from 7,090 individuals (2,594 cases and 4,496 controls) supported rs3820282 as the SNP with the strongest association for endometriosis risk (P = 1.84 × 10−5, OR = 1.244 (1.126-1.375)). SNP rs3820282 is a significant eQTL in whole blood decreasing expression of LINC00339 (also known as HSPC157) and increasing expression of CDC42 (P = 2.0 ×10−54 and 4.5x10−4 respectively). The largest effects were for two LINC00339 probes (P = 2.0 ×10−54; 1.0 × 10−34). The eQTL for LINC00339 was also observed in endometrial tissue (P = 2.4 ×10−8) with the same direction of effect for both whole blood and endometrial tissue. There was no evidence for eQTL effects for WNT4. Chromatin conformation capture provides evidence for risk SNPs interacting with the promoters of both LINC00339 and CDC4 and luciferase reporter assays suggest the risk SNP rs12038474 is located in a transcriptional silencer for CDC42 and the risk allele increases expression of CDC42. However, no effect of rs3820282 was observed in the LINC00339 expression in Ishikawa cells. Taken together, our results suggest that SNPs increasing endometriosis risk in this region act through CDC42, but further functional studies are required to rule out inverse regulation of both LINC00339 and CDC42.

Do endometriosis risk-associated single nucleotide polymorphisms (SNPs) found at the 12q22 locus have effects on vezatin ( ITALIC! VEZT) expression?The original genome-wide association study (GWAS) SNP (rs10859871), and other newly identified association signals, demonstrate strong evidence for ITALIC! cis-expression quantitative trait loci (eQTL) effects on ITALIC! VEZT expression.GWAS have identified several disease-risk loci (SNPs) associated with endometriosis. The SNP rs10859871 is located within the ITALIC! VEZT gene. ITALIC! VEZT expression is altered in the endometrium of endometriosis patients and is an excellent candidate for having a causal role in endometriosis. Most of the SNPs identified from GWAS are not located within the coding region of the genome. However, they are likely to have an effect on the regulation of gene expression. Genetic variants that affect levels of gene expression are called expression quantitative trait loci (eQTL).Samples for genotyping and ITALIC! VEZT variant screening were drawn from women recruited for genetic studies in Australia/New Zealand and women undergoing surgery in a tertiary care centre. Coding variants for ITALIC! VEZT were screened in blood from 100 unrelated individuals (endometriosis-dense families) from the QIMR Berghofer Medical Research Institute dataset. SNPs at the 12q22 locus were imputed and reanalysed for their association with endometriosis. Reanalysis of endometriosis risk-association was performed on a final combined Australian dataset of 2594 cases and 4496 controls. Gene expression was performed on 136 endometrial samples. eQTL analysis in whole blood was performed on 862 individuals from the Brisbane Systems Genetics Study. Endometrial tissue-specific eQTL analysis was performed on 122 samples (eutopic endometrium) collected following laparoscopic surgery. VEZT protein expression studies employed ITALIC! n = 56 (western blotting) and ITALIC! n = 42 (immunohistochemistry) endometrial samples.The women recruited for this study provided blood and/or endometrial tissue samples in a hospital setting. Genomic DNA was screened for common and coding variants. SNPs of interest in the 12q22 region were genotyped using Agena MassARRAY technology or Taqman SNP genotyping assay. Gene expression profiles from RNA extracted from blood and endometrial tissue samples were generated using Illumina whole-genome expression chips (Human HT-12 v4.0). Whole protein extracted from endometrium was used for VEZT western blots, and paraffin sections of endometrium were employed for VEZT immunohistochemistry semi-quantitative analysis.A total of 11 coding variants of ITALIC! VEZT (including one novel variant) were identified from an endometriosis-dense cohort. Polymorphic coding and imputed SNPs were combined with previous GWAS data to reanalyse the endometriosis risk association of the 12q22 region. The disease association signal at 12q22 was due to coding variants in ITALIC! VEZT or ITALIC! FGD6 (FYVE, RhoGEF and PH domain-containing 6) and SNPs with the strongest signals were either intronic or intergenic. We found strong evidence for ITALIC! VEZT cis-eQTLs with the sentinel SNP (rs10859871) in blood and endometrium, where the endometriosis risk allele (C) was associated with an increase in ITALIC! VEZT expression. We could not demonstrate this genotype-specific effect on VEZT protein expression in endometrium. However, we did observe a menstrual cycle stage specific increase in VEZT protein expression in endometrial glands, specific to the secretory phase ( ITALIC! P = 2.0 × 10(-4)).In comparison to the blood sample datasets, the study numbers of endometrial tissues were substantially reduced. Protein studies failed to complement RNA results, also likely a reflection of the low study numbers in these experiments. ITALIC! In silico prediction tools used in this investigation are typically based on cell lines different to our tissues of interest, thus any functional annotations drawn from these approaches should be considered carefully. Therefore, functional studies on VEZT and related pathway components are still warranted to unequivocally implicate a causal role for VEZT in endometriosis pathophysiology.GWAS have proven to be very valuable tools for deciphering complex diseases. Endometriosis is a text-book example of a complex disease, involving genetic, lifestyle and environmental influences. Our focused investigation of the 12q22 region validates an association with increased endometriosis risk. Endometriosis risk SNPs (including rs10859871) located within this locus demonstrated evidence for ITALIC! cis-eQTLs on ITALIC! VEZT expression. By examining women who possess an enhanced genetic risk of developing endometriosis, we have identified an effect on ITALIC! VEZT expression and therefore a potential gene/gene pathway in endometriosis disease establishment and development.Funding for this work was provided by NHMRC Project Grants GNT1012245, GNT1026033, GNT1049472 and GNT1046880. G.W.M. is supported by the NHMRC Fellowship scheme (GNT1078399). S.J.H.-C. is supported by the J.N. Peters Bequest Fellowship. The authors declare no competing interests.N/A.

Comparison between groups of monozygotic (MZ) and dizygotic (DZ) twins enables an estimation of the relative contribution of genetic and shared and nonshared environmental factors to phenotypic variability. Using DNA methylation profiling of ∼20,000 CpG sites as a phenotype, we have examined discordance levels in three neonatal tissues from 22 MZ and 12 DZ twin pairs. MZ twins exhibit a wide range of within-pair differences at birth, but show discordance levels generally lower than DZ pairs. Within-pair methylation discordance was lowest in CpG islands in all twins and increased as a function of distance from islands. Variance component decomposition analysis of DNA methylation in MZ and DZ pairs revealed a low mean heritability across all tissues, although a wide range of heritabilities was detected for specific genomic CpG sites. The largest component of variation was attributed to the combined effects of nonshared intrauterine environment and stochastic factors. Regression analysis of methylation on birth weight revealed a general association between methylation of genes involved in metabolism and biosynthesis, providing further support for epigenetic change in the previously described link between low birth weight and increasing risk for cardiovascular, metabolic, and other complex diseases. Finally, comparison of our data with that of several older twins revealed little evidence for genome-wide epigenetic drift with increasing age. This is the first study to analyze DNA methylation on a genome scale in twins at birth, further highlighting the importance of the intrauterine environment on shaping the neonatal epigenome.