Search

Although γδ T cells are known to participate in immune dysregulation in solid tumors, their relevance to human microsatellite-stable (MSS) colorectal cancer (CRC) is still undefined. Here, using integrated gene expression analysis and T cell receptor sequencing, we characterized γδ T cells in MSS CRC, with a focus on Vδ1 + T cells. We identified Vδ1 + T cells with shared motifs in the third complementarity-determining region of the δ-chain, reflective of antigen recognition. Changes in gene and protein expression levels suggested a dysfunctional effector state of Vδ1 + T cells in MSS CRC, distinct from Vδ1 + T cells in microsatellite-instable (MSI). Interaction analysis highlighted an immunosuppressive role of fibroblasts in the dysregulation of Vδ1 + T cells in MSS CRC via the TIGIT-NECTIN2 axis. Blocking this pathway with a TIGIT antibody partially restored cytotoxicity of the dysfunctional Vδ1 phenotype. These results define an operative pathway in γδ T cells in MSS CRC., (© 2024. The Author(s).)

Background: Aging in postmitotic tissues is associated with clonal expansion of somatic mitochondrial deletions, the origin of which is not well understood. Such deletions are often flanked by direct nucleotide repeats, but this alone does not fully explain their distribution. Here, we hypothesized that the close proximity of direct repeats on single-stranded mitochondrial DNA (mtDNA) might play a role in the formation of deletions., Results: By analyzing human mtDNA deletions in the major arc of mtDNA, which is single-stranded during replication and is characterized by a high number of deletions, we found a non-uniform distribution with a "hot spot" where one deletion breakpoint occurred within the region of 6-9 kb and another within 13-16 kb of the mtDNA. This distribution was not explained by the presence of direct repeats, suggesting that other factors, such as the spatial proximity of these two regions, can be the cause. In silico analyses revealed that the single-stranded major arc may be organized as a large-scale hairpin-like loop with a center close to 11 kb and contacting regions between 6-9 kb and 13-16 kb, which would explain the high deletion activity in this contact zone. The direct repeats located within the contact zone, such as the well-known common repeat with a first arm at 8470-8482 bp (base pair) and a second arm at 13,447-13,459 bp, are three times more likely to cause deletions compared to direct repeats located outside of the contact zone. A comparison of age- and disease-associated deletions demonstrated that the contact zone plays a crucial role in explaining the age-associated deletions, emphasizing its importance in the rate of healthy aging., Conclusions: Overall, we provide topological insights into the mechanism of age-associated deletion formation in human mtDNA, which could be used to predict somatic deletion burden and maximum lifespan in different human haplogroups and mammalian species., (© 2023. The Author(s).)

Cellular development is tightly regulated as mature cells with aberrant functions may initiate pathogenic processes. The endometrium is a highly regenerative tissue, shedding and regenerating each month. Endometrial stromal fibroblasts are regenerated each cycle from mesenchymal stem cells and play a pivotal role in endometriosis, a disease characterised by endometrial cells that grow outside the uterus. Why the cells of some women are more capable of developing into endometriosis lesions is not clear. Using isolated, purified and cultured endometrial cells of mesenchymal origin from 19 women with (n = 10) and without (n = 9) endometriosis we analysed the transcriptome of 33,758 individual cells and compared these to clinical characteristics and in vitro growth profiles. We show purified mesenchymal cell cultures include a mix of mesenchymal stem cells and two endometrial stromal fibroblast subtypes with distinct transcriptomic signatures indicative of varied progression through the differentiation processes. The fibroblast subgroup characterised by incomplete differentiation was predominantly (81%) derived from women with endometriosis and exhibited an altered in vitro growth profile. These results uncover an inherent difference in endometrial cells of women with endometriosis and highlight the relevance of cellular differentiation and its potential to contribute to disease susceptibility., (© 2022. The Author(s).)

Squamous cell carcinoma (SCC) is a common type of skin cancer that typically arises from premalignant precursor lesions named actinic keratoses (AK). Chronic inflammation is a well-known promoter of skin cancer progression. AK and SCC have been associated with an overabundance of the bacterium Staphylococcus aureus ( S. aureus ). Certain secreted products from S. aureus are known to promote cutaneous pro-inflammatory responses; however, not all S. aureus strains produce these. As inflammation plays a key role in SCC development, we investigated the pro-inflammatory potential and toxin secretion profiles of skin-cancer associated S. aureus . Sterile culture supernatants ("secretomes") of S. aureus clinical strains isolated from AK and SCC were applied to human keratinocytes in vitro. Some S. aureus secretomes induced keratinocytes to overexpress inflammatory mediators that have been linked to skin carcinogenesis, including IL-6, IL-8, and TNFα. A large phenotypic variation between the tested clinical strains was observed. Strains that are highly pro-inflammatory in vitro also caused more pronounced skin inflammation in mice. Proteomic characterization of S. aureus secretomes using mass spectrometry established that specific S. aureus enzymes and cytolytic toxins, including hemolysins, phenol-soluble modulins, and serine proteases, as well as currently uncharacterized proteins, correlate with the pro-inflammatory S. aureus phenotype. This study is the first to describe the toxin secretion profiles of AK and SCC-associated S. aureus , and their potential to induce a pro-inflammatory environment in the skin. Further studies are needed to establish whether these S. aureus products promote SCC development by mediating chronic inflammation., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Krueger, Zaugg, Chisholm, Linedale, Lachner, Teoh, Tuong, Lukowski, Morrison, Soyer, Hugenholtz, Hill and Frazer.)

Langerhans cells (LC) are skin-resident antigen-presenting cells that regulate immune responses to epithelial microorganisms. Human papillomavirus (HPV) infection can promote malignant epithelial transformation. As LCs are considered important for controlling HPV infection, we compared the transcriptome of murine LCs from skin transformed by K14E7 oncoprotein and from healthy skin. We identified transcriptome heterogeneity at the single cell level amongst LCs in normal skin, associated with ontogeny, cell cycle, and maturation. We identified a balanced co-existence of immune-stimulatory and immune-inhibitory LC cell states in normal skin that was significantly disturbed in HPV16 E7-transformed skin. Hyperplastic skin was depleted of immune-stimulatory LCs and enriched for LCs with an immune-inhibitory gene signature, and LC-keratinocyte crosstalk was dysregulated. We identified reduced expression of interleukin (IL)-34, a critical molecule for LC homeostasis. Enrichment of an immune-inhibitory LC gene signature and reduced levels of epithelial IL-34 were also found in human HPV-associated cervical epithelial cancers., Competing Interests: The authors declare no competing interests., (© 2021 The Authors.)

Here, we outline detailed protocols to isolate and profile murine splenic dendritic cells (DCs) through advanced flow cytometry of the myeloid compartment and single-cell transcriptomic profiling with integrated cell surface protein expression through CITE-seq. This protocol provides a general transferrable road map for different tissues and species. For complete details on the use and execution of this protocol, please refer to Lukowski et al. (2021)., Competing Interests: The authors declare no competing interests. S.L. is currently an employee of Boehringer-Ingelheim., (© 2021 The Author(s).)

Conventional dendritic cells (cDCs) are traditionally subdivided into cDC1 and cDC2 lineages. Batf3 is a cDC1-required transcription factor, and we observed that Batf3 -/- mice harbor a population of cDC1-like cells co-expressing cDC2-associated surface molecules. Using single-cell RNA sequencing with integrated cell surface protein expression (CITE-seq), we found that Batf3 -/- mitotic immature cDC1-like cells showed reduced expression of cDC1 features and increased levels of cDC2 features. In wild type, we also observed a proportion of mature cDC1 cells expressing surface features characteristic to cDC2 and found that overall cDC cell state heterogeneity was mainly driven by developmental stage, proliferation, and maturity. We detected population diversity within Sirpa + cDC2 cells, including a Cd33 + cell state expressing high levels of Sox4 and lineage-mixed features characteristic to cDC1, cDC2, pDCs, and monocytes. In conclusion, these data suggest that multiple cDC cell states can co-express lineage-overlapping features, revealing a level of previously unappreciated cDC plasticity., Competing Interests: The authors declare no competing interests. S.L. is currently an employee of Boehringer-Ingelheim., (© 2021 The Author(s).)

Thrombosis is a leading cause of morbidity and mortality. Fibrinogen, the soluble substrate for fibrin-based clotting, has a central role in haemostasis and thrombosis and its plasma concentration correlates with cardiovascular disease event risk and a prothrombotic state in experimental models. We aimed to identify chemical entities capable of changing fibrinogen production and test their impact on experimental thrombosis. A total of 1,280 bioactive compounds were screened for their ability to alter fibrinogen production by hepatocyte-derived cancer cells and a selected panel was tested in zebrafish larvae. Anthralin and all- trans retinoic acid (RA) were identified as fibrinogen-lowering and fibrinogen-increasing moieties, respectively. In zebrafish larvae, anthralin prolonged laser-induced venous- occlusion times and reduced thrombocyte accumulation at injury sites. RA had opposite effects. Treatment with RA, a nuclear receptor ligand, increased fibrinogen mRNA levels. Using an antisense morpholino oligonucleotide to deplete zebrafish fibrinogen, we correlated a shortening of laser-induced venous thrombosis times with RA treatment and fibrinogen protein levels. Anthralin had little effect on fibrinogen mRNA in zebrafish larvae, despite leading to lower detectable fibrinogen. Therefore, we made a proteomic scan of anthralin-treated cells and larvae. A reduced representation of proteins linked to the canonical secretory pathway was detected, suggesting that anthralin affects protein secretion. In summary, we found that chemical modulation of fibrinogen levels correlates with measured effects on experimental venous thrombosis and could be investigated as a therapeutic avenue for thrombosis prevention., Competing Interests: R. J. Fish has received travel support from Novo Nordisk. All other authors have no conflict of interest to declare., (Thieme. All rights reserved.)

Background: The discovery that somatic cells can be reprogrammed to induced pluripotent stem cells (iPSCs) has provided a foundation for in vitro human disease modelling, drug development and population genetics studies. Gene expression plays a critical role in complex disease risk and therapeutic response. However, while the genetic background of reprogrammed cell lines has been shown to strongly influence gene expression, the effect has not been evaluated at the level of individual cells which would provide significant resolution. By integrating single cell RNA-sequencing (scRNA-seq) and population genetics, we apply a framework in which to evaluate cell type-specific effects of genetic variation on gene expression., Results: Here, we perform scRNA-seq on 64,018 fibroblasts from 79 donors and map expression quantitative trait loci (eQTLs) at the level of individual cell types. We demonstrate that the majority of eQTLs detected in fibroblasts are specific to an individual cell subtype. To address if the allelic effects on gene expression are maintained following cell reprogramming, we generate scRNA-seq data in 19,967 iPSCs from 31 reprogramed donor lines. We again identify highly cell type-specific eQTLs in iPSCs and show that the eQTLs in fibroblasts almost entirely disappear during reprogramming., Conclusions: This work provides an atlas of how genetic variation influences gene expression across cell subtypes and provides evidence for patterns of genetic architecture that lead to cell type-specific eQTL effects.

Regulatory T cells (Tregs) are recruited to nonlymphoid tissues in chronic disease, including cancer, and the tissue environment is held to shape the Treg phenotype diversity. Using single-cell RNA sequencing, we examined the transcriptomic and TCR profile of Tregs recruited to hyperproliferative HPV16 E7-expressing transgenic and control nontransgenic murine skin grafts. Tregs were more abundant in E7 transgenic skin grafts than control grafts, without evidence of E7 specificity. E7 transgenic grafts attracted both Klrg1 + Tregs and Il1r2 + Tregs, which were phenotypically distinct but shared a core gene signature with previously described tumor-infiltrating Tregs. Pseudotime trajectory analysis of Tregs of defined TCR clonotypes predicted phenotypic plasticity within the skin and between the skin and draining lymph nodes. Thus, oncogene-induced hyperproliferative skin expressing a single defined non-self-antigen can attract and induce non-Ag-specific Tregs that acquire distinct regulatory phenotypes characterized by specific effector gene signatures., (Copyright © 2021 The Authors.)

Human papillomavirus (HPV) infection is associated with a range of malignancies that affect anogenital and oropharyngeal sites. α-HPVs dominantly infect basal epithelial cells of mucosal tissues, where they dysregulate cell division and local immunity. The cervix is one of the mucosal sites most susceptible to HPV infections. It consists of anatomically diverse regions, and the majority of cervical intraepithelial neoplasia and cancers arise within the cervical squamo-columnar junction where undifferentiated basal progenitor cells with stem cell properties are found. The cancer stem cell theory particularly associates tumorigenesis, invasion, dissemination, and metastasis with cancer cells exhibiting stem cell properties. In this perspective, we discuss evidence of a cervical cancer stem cell niche and explore the association of stemness related genes with 5-year survival using a publicly available transcriptomic dataset of a cervical cancer cohort. We report that poor prognosis in this cohort correlates with overexpression of a subset of stemness pathway genes, a majority of which regulate the central Focal Adhesion pathway, and are also found to be enriched in the HPV infection pathway. These observations support therapeutic targeting of stemness genes overexpressed by mucosal cells infected with high-risk HPVs., (Copyright © 2020 Budhwani, Lukowski, Porceddu, Frazer and Chandra.)

The libraries generated by high-throughput single cell RNA-sequencing (scRNA-seq) platforms such as the Chromium from 10× Genomics require considerable amounts of sequencing, typically due to the large number of cells. The ability to use these data to address biological questions is directly impacted by the quality of the sequence data. Here we have compared the performance of the Illumina NextSeq 500 and NovaSeq 6000 against the BGI MGISEQ-2000 platform using identical Single Cell 3' libraries consisting of over 70 000 cells generated on the 10× Genomics Chromium platform. Our results demonstrate a highly comparable performance between the NovaSeq 6000 and MGISEQ-2000 in sequencing quality, and the detection of genes, cell barcodes, Unique Molecular Identifiers. The performance of the NextSeq 500 was also similarly comparable to the MGISEQ-2000 based on the same metrics. Data generated by both sequencing platforms yielded similar analytical outcomes for general single-cell analysis. The performance of the NextSeq 500 and MGISEQ-2000 were also comparable for the deconvolution of multiplexed cell pools via variant calling, and detection of guide RNA (gRNA) from a pooled CRISPR single-cell screen. Our study provides a benchmark for high-capacity sequencing platforms applied to high-throughput scRNA-seq libraries., (© The Author(s) 2019. Published by Oxford University Press on behalf of NAR Genomics and Bioinformatics.)

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 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., Competing Interests: Disclosures: Dr. Molyneaux reported grants from Boehringer Ingelheim, personal fees from Hoffman-La Roche, and grants from AstraZeneca outside the submitted work. Dr. Chambers reported "other" from Seqbio Pty Ltd outside the submitted work. No other disclosures were reported., (© 2020 Byrne et al.)

Study Question: Are genetic effects on endometrial gene expression tissue specific and/or associated with reproductive traits and diseases?, Summary Answer: Analyses of RNA-sequence data and individual genotype data from the endometrium identified novel and disease associated, genetic mechanisms regulating gene expression in the endometrium and showed evidence that these mechanisms are shared across biologically similar tissues., What Is Known Already: The endometrium is a complex tissue vital for female reproduction and is a hypothesized source of cells initiating endometriosis. Understanding genetic regulation specific to, and shared between, tissue types can aid the identification of genes involved in complex genetic diseases., Study Design, Size, Duration: RNA-sequence and genotype data from 206 individuals was analysed and results were compared with large publicly available datasets., Participants/materials, Setting, Methods: RNA-sequencing and genotype data from 206 endometrial samples was used to identify the influence of genetic variants on gene expression, via expression quantitative trait loci (eQTL) analysis and to compare these endometrial eQTLs with those in other tissues. To investigate the association between endometrial gene expression regulation and reproductive traits and diseases, we conducted a tissue enrichment analysis, transcriptome-wide association study (TWAS) and summary data-based Mendelian randomisation (SMR) analyses. Transcriptomic data was used to test differential gene expression between women with and without endometriosis., Main Results and the Role of Chance: A tissue enrichment analysis with endometriosis genome-wide association study summary statistics showed that genes surrounding endometriosis risk loci were significantly enriched in reproductive tissues. A total of 444 sentinel cis-eQTLs (P < 2.57 × 10-9) and 30 trans-eQTLs (P < 4.65 × 10-13) were detected, including 327 novel cis-eQTLs in endometrium. A large proportion (85%) of endometrial eQTLs are present in other tissues. Genetic effects on endometrial gene expression were highly correlated with the genetic effects on reproductive (e.g. uterus, ovary) and digestive tissues (e.g. salivary gland, stomach), supporting a shared genetic regulation of gene expression in biologically similar tissues. The TWAS analysis indicated that gene expression at 39 loci is associated with endometriosis, including five known endometriosis risk loci. SMR analyses identified potential target genes pleiotropically or causally associated with reproductive traits and diseases including endometriosis. However, without taking account of genetic variants, a direct comparison between women with and without endometriosis showed no significant difference in endometrial gene expression., Large Scale Data: The eQTL dataset generated in this study is available at http://reproductivegenomics.com.au/shiny/endo&#95;eqtl&#95;rna/. Additional datasets supporting the conclusions of this article are included within the article and the supplementary information files, or are available on reasonable request., Limitations, Reasons for Caution: Data are derived from fresh tissue samples and expression levels are an average of expression from different cell types within the endometrium. Subtle cell-specifc expression changes may not be detected and differences in cell composition between samples and across the menstrual cycle will contribute to sample variability. Power to detect tissue specific eQTLs and differences between women with and without endometriosis was limited by the sample size in this study. The statistical approaches used in this study identify the likely gene targets for specific genetic risk factors, but not the functional mechanism by which changes in gene expression may influence disease risk., Wider Implications of the Findings: Our results identify novel genetic variants that regulate gene expression in endometrium and the majority of these are shared across tissues. This allows analysis with large publicly available datasets to identify targets for female reproductive traits and diseases. Much larger studies will be required to identify genetic regulation of gene expression that will be specific to endometrium., Study Funding/competing Interest(s): This work was supported by the National Health and Medical Research Council (NHMRC) under project grants GNT1026033, GNT1049472, GNT1046880, GNT1050208, GNT1105321, GNT1083405 and GNT1107258. G.W.M is supported by a NHMRC Fellowship (GNT1078399). J.Y is supported by an ARC Fellowship (FT180100186). There are no competing interests., (© The Author(s) 2020. Published by Oxford University Press on behalf of the European Society of Human Reproduction and Embryology.)

Human Papillomavirus infection is highly prevalent worldwide. While most types of HPV cause benign warts, some high-risk types are known to cause cervical cancer, as well as cancer of the oral cavity and head and neck. Persistent cutaneous HPV infection can be particularly problematic in patients with chronic immunosuppression, for example following organ transplantation. Due to unknown mechanisms, these patients may develop numerous warts, as well as present with a dramatically increased skin cancer prevalence. Despite an association between HPV persistence in the epidermis and excessive wart or squamous cancer development, the molecular mechanisms linking immunosuppression, HPV expression and excessive epidermal proliferation have not been determined, largely due to low-sensitivity methodology to capture rare viral transcription events. Here, we use single-cell RNA sequencing to profile HPV-positive skin lesions from an immunosuppressed patient that were found to express the alphapapillomavirus HPV78 in basal keratinocytes, suprabasal keratinocytes and hair follicle stem cells. This method can be applied to detect and investigate HPV transcripts in cutaneous lesions, allowing mechanistic links between immunosuppression-induced HPV life cycle and epidermal hyperproliferation to be uncovered., (Copyright © 2019 Elsevier Inc. All rights reserved.)

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, Müller 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 Müller 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., (© 2019 The Authors.)

Background: Recent developments in single-cell RNA sequencing (scRNA-seq) platforms have vastly increased the number of cells typically assayed in an experiment. Analysis of scRNA-seq data is multidisciplinary in nature, requiring careful consideration of the application of statistical methods with respect to the underlying biology. Few analysis packages exist that are at once robust, are computationally fast, and allow flexible integration with other bioinformatics tools and methods., Findings: ascend is an R package comprising tools designed to simplify and streamline the preliminary analysis of scRNA-seq data, while addressing the statistical challenges of scRNA-seq analysis and enabling flexible integration with genomics packages and native R functions, including fast parallel computation and efficient memory management. The package incorporates both novel and established methods to provide a framework to perform cell and gene filtering, quality control, normalization, dimension reduction, clustering, differential expression, and a wide range of visualization functions., Conclusions: ascend is designed to work with scRNA-seq data generated by any high-throughput platform and includes functions to convert data objects between software packages. The ascend workflow is simple and interactive, as well as suitable for implementation by a broad range of users, including those with little programming experience., (© The Author(s) 2019. Published by Oxford University Press.)

The cellular and molecular profiles that govern the endothelial heterogeneity of the circulatory system have yet to be elucidated. Using a data-driven approach to study the endothelial compartment via single-cell RNA sequencing, we characterized cell subpopulations within and assigned them to a defined endothelial hierarchy. We show that two transcriptionally distinct endothelial populations exist within the aorta and, using two independent trajectory analysis methods, confirm that they represent transitioning cells rather than discrete cell types. Gene co-expression analysis revealed crucial regulatory networks underlying each population, including significant metabolic gene networks in progenitor cells. Using mitochondrial activity assays and phenotyping, we confirm that endovascular progenitors display higher mitochondrial content compared to differentiated endothelial cells. The identities of these populations were further validated against bulk RNA sequencing (RNA-seq) data obtained from normal and tumor-derived vasculature. Our findings validate the heterogeneity of the aortic endothelium and previously suggested hierarchy between progenitor and differentiated cells., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

Background: Major challenges in understanding the functional consequences of genetic risk factors for human disease are which tissues and cell types are affected and the limited availability of suitable tissue. The aim of this study was to evaluate tissue-specific genotype-epigenetic characteristics in DNA samples from both endometrium and blood collected from women at different stages of the menstrual cycle and relate results to genetic risk factors for reproductive traits and diseases., Results: We analysed DNA methylation (DNAm) data from endometrium and blood samples from 66 European women. Methylation profiles were compared between stages of the menstrual cycle, and changes in methylation overlaid with changes in transcription and genotypes. We observed large changes in methylation (27,262 DNAm probes) across the menstrual cycle in endometrium that were not observed in blood. Individual genotype data was tested for association with methylation at 443,016 and 443,101 DNAm probes in endometrium and blood respectively to identify methylation quantitative trait loci (mQTLs). A total of 4546 sentinel cis-mQTLs (P < 1.13 × 10 -10 ) and 434 sentinel trans-mQTLs (P < 2.29 × 10 -12 ) were detected in endometrium and 6615 sentinel cis-mQTLs (P < 1.13 × 10 -10 ) and 590 sentinel trans-mQTLs (P < 2.29 × 10 -12 ) were detected in blood. Following secondary analyses, conducted to test for overlap between mQTLs in the two tissues, we found that 62% of endometrial cis-mQTLs were also observed in blood and the genetic effects between tissues were highly correlated. A number of mQTL SNPs were associated with reproductive traits and diseases, including one mQTL located in a known risk region for endometriosis (near GREB1)., Conclusions: We report novel findings characterising genetic regulation of methylation in endometrium and the association of endometrial mQTLs with endometriosis risk and other reproductive traits and diseases. The high correlation of genetic effects between tissues highlights the potential to exploit the power of large mQTL datasets in endometrial research and identify target genes for functional studies. However, tissue-specific methylation profiles and genetic effects also highlight the importance of also using disease-relevant tissues when investigating molecular mechanisms of disease risk.

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., (© 2018 Asian Pacific Society of Respirology.)

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., (Copyright © 2018 The Authors. Published by Elsevier Inc. All rights reserved.)

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., (Copyright © 2018 Elsevier Inc. All rights reserved.)

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.

Background: Giant cell arteritis (GCA) is the most common form of vasculitis affecting elderly people. It is one of the few true ophthalmic emergencies but symptoms and signs are variable thereby making it a challenging disease to diagnose. A temporal artery biopsy is the gold standard to confirm GCA, but there are currently no specific biochemical markers to aid diagnosis. 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: We recruited 16 patients with histological evidence of GCA at the Royal Victorian Eye and Ear Hospital, 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. 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., Results: 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, two transcripts remained differentially expressed after 12 months; SGTB, associated with neuronal apoptosis, and FCGR3A, associatied with Takayasu arteritis. We detected genes that correlate with both symptoms and biochemical markers used 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 including Polymyalgia Rheumatica, bilateral blindness and death within 12 months., Conclusions: This 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. In the acute phase, the gene-phenotype relationships we have identified could provide insight to potential disease severity and as such guide in initiating appropriate patient management.

Heterogeneity of cell states represented in pluripotent cultures has 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 composed of 165 unique genes that denote the specific pluripotency states; 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 threefold) 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., (© 2018 Nguyen et al.; Published by Cold Spring Harbor Laboratory Press.)

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.

The majority of aneuploid fetuses are spontaneously miscarried. Nevertheless, some aneuploid individuals survive despite the strong genetic insult. Here, we investigate if the survival probability of aneuploid fetuses is affected by the genome-wide burden of slightly deleterious variants. We analyzed two cohorts of live-born Down syndrome individuals (388 genotyped samples and 16 fibroblast transcriptomes) and observed a deficit of slightly deleterious variants on Chromosome 21 and decreased transcriptome-wide variation in the expression level of highly constrained genes. We interpret these results as signatures of embryonic selection, and propose a genetic handicap model whereby an individual bearing an extremely severe deleterious variant (such as aneuploidy) could escape embryonic lethality if the genome-wide burden of slightly deleterious variants is sufficiently low. This approach can be used to study the composition and effect of the numerous slightly deleterious variants in humans and model organisms., (© 2018 Popadin et al.; Published by Cold Spring Harbor Laboratory Press.)

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.

Study Question: Do genetic effects regulate gene expression in human endometrium?, Summary Answer: 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., What Is Known Already: 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., Study Design, Size, Duration: We analysed genome-wide genotype and gene expression data to map cis expression quantitative trait loci (eQTL) in endometrium., Participants/materials, Setting, Methods: 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., Main Results and the Role of Chance: We identified a total of 18 595 cis SNP-probe associations at a study-wide level of significance (P < 1 × 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., Large Scale Data: The summary eQTL results are publicly available to browse or download., Limitations, Reasons for Caution: 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., Wider Implications of the Findings: Identification of endometrial eQTLs provides a platform for better understanding genetic effects on endometriosis risk and other endometrial-related pathologies., Study Funding/competing Interest(s): Funding for this work was provided by NHMRC Project Grants GNT1026033, GNT1049472, GNT1046880, GNT1050208, GNT1105321 and APP1083405. There are no competing interests., (© The Author 2017. Published by Oxford University Press on behalf of the European Society of Human Reproduction and Embryology. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com)

Introduction: Fibrinogen storage disease (FSD) is characterized by hypofibrinogenemia and hepatic inclusions due to impaired release of mutant fibrinogen which accumulates and aggregates in the hepatocellular endoplasmic reticulum. Liver disease is variable., Aim: We studied a new Swiss family with fibrinogen Aguadilla. In order to understand the molecular peculiarity of FSD mutations, fibrinogen Aguadilla and the three other causative mutations, all located in the γD domain, were modelled., Method: The proband is a Swiss girl aged 4 investigated because of fatigue and elevated liver enzymes. Protein structure models were prepared using the Swiss-PdbViewer and POV-Ray software., Results: The proband was found to be heterozygous for fibrinogen Aguadilla: FGG Arg375Trp. Familial screening revealed that her mother and maternal grandmother were also affected and, in addition, respectively heterozygous and homozygous for the hereditary haemochromatosis mutation HFE C282Y. Models of backbone and side-chain interactions for fibrinogen Aguadilla in a 10-angstrom region revealed the loss of five H-bonds and the gain of one H-bond between structurally important amino acids. The structure predicted for fibrinogen Angers showed a novel helical structure in place of hole 'a' on the outer edge of γD likely to have a negative impact on fibrinogen assembly and secretion., Conclusion: The mechanism by which FSD mutations generate hepatic intracellular inclusions is still not clearly established although the promotion of aberrant intermolecular strand insertions is emerging as a likely cause. Reporting new cases is essential in the light of novel opportunities of treatment offered by increasing knowledge of the degradation pathway and autophagy., (© 2015 John Wiley & Sons Ltd.)

The expression of plasma proteins changes dramatically as a result of cytokine induction, particularly interleukin-6, and their levels are used as clinical markers of inflammation. miRNAs are important regulators of gene expression and play significant roles in many inflammatory diseases and processes. The interactions between miRNAs and the genes that they regulate during the acute phase response have not been investigated. We examined the effects of IL-6 stimulation on the transcriptome and miRNome of human and mouse primary hepatocytes and the HepG2 cell line. Using an integrated analysis, we identified differentially expressed miRNAs whose seed sequences are significantly enriched in the 3' untranslated regions of differentially expressed genes, many of which are involved in inflammation-related pathways. Our finding that certain miRNAs may de-repress critical acute phase proteins within acute timeframes has important biological and clinical implications., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Post-transcriptional regulation of gene expression through 5' untranslated region (5'UTR)-encoded cis-acting elements is an important mechanism for the control of protein expression levels. Through controlling specific aspects of translation initiation, expression can be tightly regulated while remaining responsive to cellular requirements. With respect to cystic fibrosis (CF), the overexpression of cystic fibrosis transmembrane conductance regulator (CFTR) protein trafficking mutants, such as delta-F508, is of great biological and clinical interest. By understanding the post-transcriptional mechanisms that regulate CFTR expression, new procedures can be developed to enhance CFTR expression in homozygous delta-F508 CF patients. We have identified the key elements of a complex negative regulatory mechanism that is encoded within the human CFTR 5'UTR and show how these elements act in combination to restrict CFTR gene expression to a consistently low level in a transcript-specific manner. This study shows, for the first time, that endogenous human CFTR expression is post-transcriptionally regulated through a 5'UTR-mediated mechanism. We show that the very low levels of endogenous CFTR expression, compared with other low expression genes, are maintained through the co-operative inhibitory effects of an upstream open reading frame and a thermodynamically stable RNA secondary structure., (© The Author 2014. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Mammals express the sialic acids N-acetylneuraminic acid (Neu5Ac) and N-glycolylneuraminic acid (Neu5Gc) on cell surfaces, where they act as receptors for pathogens, including influenza A virus (IAV). Neu5Gc is synthesized from Neu5Ac by the enzyme cytidine monophosphate-N-acetylneuraminic acid hydroxylase (CMAH). In humans, this enzyme is inactive and only Neu5Ac is produced. Ferrets are susceptible to human-adapted IAV strains and have been the dominant animal model for IAV studies. Here we show that ferrets, like humans, do not synthesize Neu5Gc. Genomic analysis reveals an ancient, nine-exon deletion in the ferret CMAH gene that is shared by the Pinnipedia and Musteloidia members of the Carnivora. Interactions between two human strains of IAV with the sialyllactose receptor (sialic acid--α2,6Gal) confirm that the type of terminal sialic acid contributes significantly to IAV receptor specificity. Our results indicate that exclusive expression of Neu5Ac contributes to the susceptibility of ferrets to human-adapted IAV strains.

Gene expression levels can be subject to selection. We hypothesized that the age of gene origin is associated with expression constraints, given that it affects the level of gene integration into the functional cellular environment. By studying the genetic variation affecting gene expression levels (cis expression quantitative trait loci [cis-eQTLs]) and protein levels (cis protein QTLs [cis-pQTLs]), we determined that young, primate-specific genes are enriched in cis-eQTLs and cis-pQTLs. Compared to cis-eQTLs of old genes originating before the zebrafish divergence, cis-eQTLs of young genes have a higher effect size, are located closer to the transcription start site, are more significant, and tend to influence genes in multiple tissues and populations. These results suggest that the expression constraint of each gene increases throughout its lifespan. We also detected a positive correlation between expression constraints (approximated by cis-eQTL properties) and coding constraints (approximated by Ka/Ks) and observed that this correlation might be driven by gene age. To uncover factors associated with the increase in gene-age-related expression constraints, we demonstrated that gene connectivity, gene involvement in complex regulatory networks, gene haploinsufficiency, and the strength of posttranscriptional regulation increase with gene age. We also observed an increase in heritability of gene expression levels with age, implying a reduction of the environmental component. In summary, we show that gene age shapes key gene properties during evolution and is therefore an important component of genome function., (Copyright © 2014 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Escherichia coli ST131 is a globally disseminated, multidrug resistant clone responsible for a high proportion of urinary tract and bloodstream infections. The rapid emergence and successful spread of E. coli ST131 is strongly associated with antibiotic resistance; however, this phenotype alone is unlikely to explain its dominance amongst multidrug resistant uropathogens circulating worldwide in hospitals and the community. Thus, a greater understanding of the molecular mechanisms that underpin the fitness of E. coli ST131 is required. In this study, we employed hyper-saturated transposon mutagenesis in combination with multiplexed transposon directed insertion-site sequencing to define the essential genes required for in vitro growth and the serum resistome (i.e. genes required for resistance to human serum) of E. coli EC958, a representative of the predominant E. coli ST131 clonal lineage. We identified 315 essential genes in E. coli EC958, 231 (73%) of which were also essential in E. coli K-12. The serum resistome comprised 56 genes, the majority of which encode membrane proteins or factors involved in lipopolysaccharide (LPS) biosynthesis. Targeted mutagenesis confirmed a role in serum resistance for 46 (82%) of these genes. The murein lipoprotein Lpp, along with two lipid A-core biosynthesis enzymes WaaP and WaaG, were most strongly associated with serum resistance. While LPS was the main resistance mechanism defined for E. coli EC958 in serum, the enterobacterial common antigen and colanic acid also impacted on this phenotype. Our analysis also identified a novel function for two genes, hyxA and hyxR, as minor regulators of O-antigen chain length. This study offers novel insight into the genetic make-up of E. coli ST131, and provides a framework for future research on E. coli and other Gram-negative pathogens to define their essential gene repertoire and to dissect the molecular mechanisms that enable them to survive in the bloodstream and cause disease., Competing Interests: The authors have declared that no competing interests exist.

Cystic fibrosis (CF) is characterized as a single-gene disorder with a simple, autosomal recessive mode of inheritance. However, translation of cystic fibrosis transmembrane conductance regulator (CFTR) genotype into CF phenotype is influenced by nucleotide sequence variations at multiple genetic loci, and individuals heterozygous for CFTR mutations are predisposed to a range of CFTR-related conditions, such as disseminated bronchiectasis. CF disease severity and CFTR-related conditions are more akin to complex, multifactorial traits, which are increasingly being associated with mutations that perturb gene expression. We have identified a patient with disseminated bronchiectasis, who is heterozygous for a single-nucleotide substitution in the CFTR 5' untranslated region (UTR) (c.-34C>T). The c.-34C>T mutation creates an upstream AUG codon and upstream open reading frame that overlaps, and is out of frame with, the CFTR protein coding sequence. Using luciferase reporter constructs, we have shown that the c.-34C>T mutation decreases gene expression by 85-99%, by reducing translation efficiency and mRNA stability. This is the first CFTR regulatory mutation shown to act at a posttranscriptional level that reduces the synthesis of normal CFTR (Class V), and reaffirms the importance of regulatory mutations as a genetic basis of multifactorial phenotypes., (© 2011 Wiley-Liss, Inc.)