Standardization and Interpretation of RNA-sequencing for Transplantation.

RNA-sequencing (RNA-seq) is a technique to determine the order of nucleotides in an RNA segment. Modern sequencing platforms simultaneously sequence millions of RNA molecules. Advances in bioinformatics have allowed us to collect, store, analyze, and disseminate data from RNA-seq experiments and decipher biological insights from large sequencing datasets. Although bulk RNA-seq has significantly advanced our understanding of tissue-specific gene expression and regulation, recent advances in single-cell RNA-seq have allowed such information to be mapped to individual cells, thus remarkably enhancing our insight into discrete cellular functions within a biospecimen. These different RNA-seq experimental approaches require specialized computational tools. Herein, we will first review the RNA-seq experimental workflow, discuss the common terminologies used in RNA-seq, and suggest approaches for standardization across multiple studies. Next, we will provide an up-to-date appraisal of the applications of bulk RNA-seq and single-cell/nucleus RNA-seq in preclinical and clinical research on kidney transplantation, as well as typical bioinformatic workflows utilized in such analysis. Lastly, we will deliberate on the limitations of this technology in transplantation research and briefly summarize newer technologies that could be combined with RNA-seq to permit more powerful dissections of biological functions. Because each step in RNA-seq workflow has numerous variations and could potentially impact the results, as conscientious citizens of the research community, we must strive to continuously modernize our analytical pipelines and exhaustively report their technical details.

Detection of infiltrating fibroblasts by single-cell transcriptomics in human kidney allografts.

We tested the hypothesis that single-cell RNA-sequencing (scRNA-seq) analysis of human kidney allograft biopsies will reveal distinct cell types and states and yield insights to decipher the complex heterogeneity of alloimmune injury. We selected 3 biopsies of kidney cortex from 3 individuals for scRNA-seq and processed them fresh using an identical protocol on the 10x Chromium platform; (i) HK: native kidney biopsy from a living donor, (ii) AK1: allograft kidney with transplant glomerulopathy, tubulointerstitial fibrosis, and worsening graft function, and (iii) AK2: allograft kidney after successful treatment of active antibody-mediated rejection. We did not study T-cell-mediated rejections. We generated 7217 high-quality single cell transcriptomes. Taking advantage of the recipient-donor sex mismatches revealed by X and Y chromosome autosomal gene expression, we determined that in AK1 with fibrosis, 42 months after transplantation, more than half of the kidney allograft fibroblasts were recipient-derived and therefore likely migratory and graft infiltrative, whereas in AK2 without fibrosis, 84 months after transplantation, most fibroblasts were donor-organ-derived. Furthermore, AK1 was enriched for tubular progenitor cells overexpressing profibrotic extracellular matrix genes. AK2, eight months after successful treatment of rejection, contained plasmablast cells with high expression of immunoglobulins, endothelial cell elaboration of T cell chemoattractant cytokines, and persistent presence of cytotoxic T cells. In addition to these key findings, our analysis revealed unique cell types and states in the kidney. Altogether, single-cell transcriptomics yielded novel mechanistic insights, which could pave the way for individualizing the care of transplant recipients.

Dynamic genome-wide gene expression and immune cell composition in the developing human placenta.

Despite the central role of the placenta in supporting a pregnancy, relatively little is known about transcriptomic and immune-cell changes that occur across gestation. To generate a reference gene expression map of first (T1), second (T2) and third (T3) trimester human placenta, and assess differences in transcriptome in maternal versus fetal side tissues sections of full-term placenta, we performed RNA-Seq analysis on 64 biopsy samples from 18 placentas across all three gestations. We identified 1120 differentially expressed genes in placenta tissues from T1 and T3 samples using a generalized linear model within DESeq2. In total, 411 and 709 genes were positively associated with T1 and T3 placenta, respectively. Comparison of the top 200 differentially expressed genes in the T1 placenta with T3 showed that most of the top enriched biological processes were related to cell division and proliferation. T1 and T2 tissues shared expression of fibroblast-specific COL6A2, HGF, and SPP1 genes. In T3 samples, the expression of genes relating to vasculature development and regulation were highly enriched. Monocytes and NK cell population increased in T3 compared to T1 and T2, whereas Hofbauer cell proportion expanded significantly in T2 and then decreased in T3 samples. There were no significant gene expression differences in the maternal vs. fetal side in T3 placentas. Gene expression patterns shift temporally across trimesters but not spatially across the placenta, at least at the resolution of the biopsy samples. The genes and gene set we identified here represent a valuable resource for studying pathology in pregnancy-related disorders.

Invasive phenotype in triple negative breast cancer is inhibited by blocking SIN3A-PF1 interaction through KLF9 mediated repression of ITGA6 and ITGB1.

SIN3A, a scaffold protein has regulatory functions in tumor biology. Through its Paired amphipathic helix (PAH2) domain, SIN3A interacts with PHF12 (PF1), a protein with SIN3 interaction domain (SID) that forms a complex with MRG15 and KDM5A/B. These components are often overexpressed in cancer. In the present study, we evaluated the role of SIN3A and its interacting partner PF1 in mediating inhibition of tumor growth and invasion in triple negative breast cancer (TNBC). We found profound inhibition of invasion, migration, and induction of cellular senescence by specific disruption of the PF1/SIN3A PAH2 domain interaction in TNBC cells expressing PF1-SID transcript or peptide treatment. Genome-wide transcriptomic analysis by RNA-seq revealed that PF1-SID downregulates several gene sets and pathways linked to invasion and migration. Integrin α6 (ITGA6) and integrin ß1 (ITGB1) and their downstream target proteins were downregulated in PF1-SID cells. We further determined increased presence of SIN3A and transcriptional repressor, KLF9, on promoters of ITGA6 and ITGB1 in PF1-SID cells. Knockdown of KLF9 leads to re-expression of ITGA6 and ITGB1 and restoration of the invasive phenotype, functionally linking KLF9 to this process. Overall, these data demonstrate that specific disruption of PF1/SIN3A, inhibits tumor growth, migration, and invasion. Also, PF1-SID not only inhibits tumor growth by senescence induction and reduced proliferation, but it also targets cancer stem cell gene expression and blocks mammosphere formation. Overall, these data demonstrate a mechanism whereby invasion and metastasis of TNBC can be suppressed by inhibiting SIN3A-PF1 interaction and enhancing KLF9 mediated suppression of ITGA6 and ITGB1.

Plasma microRNA Interindividual Variability in Healthy Individuals, Pregnant Women, and an Individual with a Stably Altered Neuroendocrine Phenotype.

BACKGROUND: Extracellular RNAs (exRNAs) in biofluids are amenable to quantitative analysis and proposed as noninvasive biomarkers for monitoring organ function. Cell-lineage-specific microRNAs (miRNAs) are present in plasma as soluble ribonucleoproteins or enclosed in exRNA carriers and transported through the vasculature. However, more extensive studies of healthy individuals are needed to gain insights into the variability of plasma miRNA abundance and composition. METHODS: The exRNA composition of platelet-depleted plasma collected twice from 236 healthy individuals was characterized by small RNA sequencing. Plasma of pregnant women featuring dramatically increased placental miRNAs and samples from subject P12 with noticeably increased epithelial- and neuroendocrine-origin miRNAs were included for comparison. The miRNA content of 10 000g and 100 000g pellet fractions of plasma generated by ultracentrifugation was also determined. Data analysis methods included Pearson correlation, differential gene expression, and unsupervised clustering. RESULTS: The abundance changes for more variable miRNAs in plasma of normal individuals correlated between coexpressed cell-lineage-specific miRNAs of the liver, neuroendocrine organs, epithelial cells, and muscle. ExRNA of pellet fractions contained <2% of total plasma miRNA with modest enrichment of lineage-specific and variable miRNAs compared to supernatant. The abundance fold changes of miRNAs observed in pregnancy and P12 compared to normal exceeded interquartile variability of healthy individuals. The neuroendocrine miRNA signature of P12 persisted for more than 4 years and was absent in other individuals. CONCLUSIONS: This study defines the framework and effect size for screening of extensive plasma collections for miRNA phenotypes and biomarker discovery.

Field-deployable, rapid diagnostic testing of saliva for SARS-CoV-2.

To safely re-open economies and prevent future outbreaks, rapid, frequent, point-of-need, SARS-CoV-2 diagnostic testing is necessary. However, existing field-deployable COVID-19 testing methods require the use of uncomfortable swabs and trained providers in PPE, while saliva-based methods must be transported to high complexity laboratories for testing. Here, we report the development and clinical validation of High-Performance Loop-mediated isothermal Amplification (HP-LAMP), a rapid, saliva-based, SARS-CoV-2 test with a limit of detection of 1.4 copies of virus per µl of saliva and a sensitivity and specificity with clinical samples of > 96%, on par with traditional RT-PCR based methods using swabs, but can deliver results using only a single fluid transfer step and simple heat block. Testing of 120 patient samples in 40 pools comprised of 5 patient samples each with either all negative or a single positive patient sample was 100% accurate. Thus, HP-LAMP may enable rapid and accurate results in the field using saliva, without need of a high-complexity laboratory.

Author Correction: Comprehensive aptamer-based screening identifies a spectrum of urinary biomarkers of lupus nephritis across ethnicities.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Comprehensive aptamer-based screening identifies a spectrum of urinary biomarkers of lupus nephritis across ethnicities.

Emerging urinary biomarkers continue to show promise in evaluating lupus nephritis (LN). Here, we screen urine from active LN patients for 1129 proteins using an aptamer-based platform, followed by ELISA validation in two independent cohorts comprised of 127 inactive lupus, 107 active LN, 67 active non-renal lupus patients and 74 healthy controls, of three different ethnicities. Urine proteins that best distinguish active LN from inactive disease are ALCAM, PF-4, properdin, and VCAM-1 among African-Americans, sE-selectin, VCAM-1, BFL-1 and Hemopexin among Caucasians, and ALCAM, VCAM-1, TFPI and PF-4 among Asians. Most of these correlate significantly with disease activity indices in the respective ethnic groups, and surpass conventional metrics in identifying active LN, with better sensitivity, and negative/positive predictive values. Several elevated urinary molecules are also expressed within the kidneys in LN, based on single-cell RNAseq analysis. Longitudinal studies are warranted to assess the utility of these biomarkers in tracking lupus nephritis.

Single-cell transcriptome analysis of human skin identifies novel fibroblast subpopulation and enrichment of immune subsets in atopic dermatitis.

BACKGROUND: Atopic dermatitis (AD) is a prevalent inflammatory skin disease with a complex pathogenesis involving immune cell and epidermal abnormalities. Despite whole tissue biopsy studies that have advanced the mechanistic understanding of AD, single cell-based molecular alterations are largely unknown. OBJECTIVE: Our aims were to construct a detailed, high-resolution atlas of cell populations and assess variability in cell composition and cell-specific gene expression in the skin of patients with AD versus in controls. METHODS: We performed single-cell RNA sequencing on skin biopsy specimens from 5 patients with AD (4 lesional samples and 5 nonlesional samples) and 7 healthy control subjects, using 10× Genomics. RESULTS: We created transcriptomic profiles for 39,042 AD (lesional and nonlesional) and healthy skin cells. Fibroblasts demonstrated a novel COL6A5+COL18A1+ subpopulation that was unique to lesional AD and expressed CCL2 and CCL19 cytokines. A corresponding LAMP3+ dendritic cell (DC) population that expressed the CCL19 receptor CCR7 was also unique to AD lesions, illustrating a potential role for fibroblast signaling to immune cells. The lesional AD samples were characterized by expansion of inflammatory DCs (CD1A+FCER1A+) and tissue-resident memory T cells (CD69+CD103+). The frequencies of type 2 (IL13+)/type 22 (IL22+) T cells were higher than those of type 1 (IFNG+) in lesional AD, whereas this ratio was slightly diminished in nonlesional AD and further diminished in controls. CONCLUSION: AD lesions were characterized by expanded type 2/type 22 T cells and inflammatory DCs, and by a unique inflammatory fibroblast that may interact with immune cells to regulate lymphoid cell organization and type 2 inflammation.

Cell atlas of the foetal human heart and implications for autoimmune-mediated congenital heart block.

AIMS: Investigating human heart development and applying this to deviations resulting in disease is incomplete without molecular characterization of the cell types required for normal functioning. We investigated foetal human heart single-cell transcriptomes from mid-gestational healthy and anti-SSA/Ro associated congenital heart block (CHB) samples. METHODS AND RESULTS: Three healthy foetal human hearts (19th to 22nd week of gestation) and one foetal heart affected by autoimmune-associated CHB (21st week of gestation) were subjected to enzymatic dissociation using the Langendorff preparation to obtain single-cell suspensions followed by 10× Genomics- and Illumina-based single-cell RNA-sequencing (scRNA-seq). In addition to the myocytes, fibroblasts, immune cells, and other minor cell types, previously uncharacterized diverse sub-populations of endothelial cells were identified in the human heart. Differential gene expression analysis revealed increased and heterogeneous interferon responses in varied cell types of the CHB heart compared with the healthy controls. In addition, we also identified matrisome transcripts enriched in CHB stromal cells that potentially contribute to extracellular matrix deposition and subsequent fibrosis. CONCLUSION: These data provide an information-rich resource to further our understanding of human heart development, which, as illustrated by comparison to a heart exposed to a maternal autoimmune environment, can be leveraged to provide insight into the pathogenesis of disease.

Author Correction: Tubular cell and keratinocyte single-cell transcriptomics applied to lupus nephritis reveal type I IFN and fibrosis relevant pathways.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Single-cell RNA sequencing for the study of lupus nephritis.

Single-cell RNA sequencing (scRNA-seq) has recently undergone rapid advances in the development of this technology, leading to high throughput and accelerating discovery in many biological systems and diseases. The single-cell resolution of the technique allows for the investigation of heterogeneity in cell populations, and the pinpointing of pathological populations contributing to disease. Here we review the development of scRNA-seq technology and the analysis that has evolved with the ever-increasing throughput. Finally, we highlight recent applications of scRNA-seq to understand the molecular pathogenesis of lupus and lupus nephritis.

Tubular cell and keratinocyte single-cell transcriptomics applied to lupus nephritis reveal type I IFN and fibrosis relevant pathways.

The molecular and cellular processes that lead to renal damage and to the heterogeneity of lupus nephritis (LN) are not well understood. We applied single-cell RNA sequencing (scRNA-seq) to renal biopsies from patients with LN and evaluated skin biopsies as a potential source of diagnostic and prognostic markers of renal disease. Type I interferon (IFN)-response signatures in tubular cells and keratinocytes distinguished patients with LN from healthy control subjects. Moreover, a high IFN-response signature and fibrotic signature in tubular cells were each associated with failure to respond to treatment. Analysis of tubular cells from patients with proliferative, membranous and mixed LN indicated pathways relevant to inflammation and fibrosis, which offer insight into their histologic differences. In summary, we applied scRNA-seq to LN to deconstruct its heterogeneity and identify novel targets for personalized approaches to therapy.

A single-cell survey of the human first-trimester placenta and decidua.

The placenta and decidua interact dynamically to enable embryonic and fetal development. Here, we report single-cell RNA sequencing of 14,341 and 6754 cells from first-trimester human placental villous and decidual tissues, respectively. Bioinformatic analysis identified major cell types, many known and some subtypes previously unknown in placental villi and decidual context. Further detailed analysis revealed proliferating subpopulations, enrichment of cell type-specific transcription factors, and putative intercellular communication in the fetomaternal microenvironment. This study provides a blueprint to further the understanding of the roles of these cells in the placenta and decidua for maintenance of early gestation as well as pathogenesis in pregnancy-related disorders.

Single cell RNA sequencing to dissect the molecular heterogeneity in lupus nephritis.

Lupus nephritis is a leading cause of mortality among systemic lupus erythematosus (SLE) patients, and its heterogeneous nature poses a significant challenge to the development of effective diagnostics and treatments. Single cell RNA sequencing (scRNA-seq) offers a potential solution to dissect the heterogeneity of the disease and enables the study of similar cell types distant from the site of renal injury to identify novel biomarkers. We applied scRNA-seq to human renal and skin biopsy tissues and demonstrated that scRNA-seq can be performed on samples obtained during routine care. Chronicity index, IgG deposition, and quantity of proteinuria correlated with a transcriptomic-based score composed of IFN-inducible genes in renal tubular cells. Furthermore, analysis of cumulative expression profiles of single cell keratinocytes dissociated from nonlesional, non-sun-exposed skin of patients with lupus nephritis also revealed upregulation of IFN-inducible genes compared with keratinocytes isolated from healthy controls. This indicates the possible use of scRNA-seq analysis of skin biopsies as a biomarker of renal disease. These data support the potential utility of scRNA-seq to provide new insights into the pathogenesis of lupus nephritis and pave the way for exploiting a readily accessible tissue to reflect injury in the kidney.

DND1 maintains germline stem cells via recruitment of the CCR4-NOT complex to target mRNAs.

The vertebrate-conserved RNA-binding protein DND1 is required for the survival of primordial germ cells (PGCs), as well as the suppression of germ cell tumours in mice. Here we show that in mice DND1 binds a UU(A/U) trinucleotide motif predominantly in the 3' untranslated regions of mRNA, and destabilizes target mRNAs through direct recruitment of the CCR4-NOT deadenylase complex. Transcriptomic analysis reveals that the extent of suppression is dependent on the number of DND1-binding sites. This DND1-dependent mRNA destabilization is required for the survival of mouse PGCs and spermatogonial stem cells by suppressing apoptosis. The spectrum of target RNAs includes positive regulators of apoptosis and inflammation, and modulators of signalling pathways that regulate stem-cell pluripotency, including the TGFß superfamily, all of which are aberrantly elevated in DND1-deficient PGCs. We propose that the induction of the post-transcriptional suppressor DND1 synergizes with concurrent transcriptional changes to ensure precise developmental transitions during cellular differentiation and maintenance of the germ line.

A Simple Alternative to Stereotactic Injection for Brain Specific Knockdown of miRNA.

MicroRNAs (miRNAs) are key regulators of gene expression. In the brain, vital processes like neurodevelopment and neuronal functions depend on the correct expression of microRNAs. Perturbation of microRNAs in the brain can be used to model neurodegenerative diseases by modulating neuronal cell death. Currently, stereotactic injection is used to deliver miRNA knockdown agents to specific location in the brain. Here, we discuss strategies to design antagomirs against miRNA with locked nucleotide modifications (LNA). Subsequently describe a method for brain specific delivery of antagomirs, uniformly across different regions of the brain. This method is simple and widely applicable since it overcomes the surgery, associated injury and limitation of local delivery in stereotactic injections. We prepared a complex of neurotropic, cell-penetrating peptide Rabies Virus Glycoprotein (RVG) with antagomir against miRNA-29 and injected through tail vein, to specifically deliver in the brain. The antagomir design incorporated features that allow specific targeting of the miRNA and formation of non-covalent complexes with the peptide. The knock-down of the miRNA in neuronal cells, resulted in apoptotic cell death and associated behavioural defects. Thus, the method can be used for acute models of neuro-degeneration through the perturbation of miRNAs.

Mammalian miRNA curation through next-generation sequencing.

Characteristic small RNA biogenesis processing patterns are used for the discovery of novel microRNAs (miRNAs) from next-generation sequencing data. Here, we highlight and discuss key criteria for mammalian - specifically human - miRNA database curation based on small RNA sequencing data. Sequence reads obtained from small RNA cDNA libraries are aligned to reference genomic regions, and miRNA genes are revealed by their distinct read length and bimodal read frequency distribution, the predicted secondary structure of the deduced miRNA stem-loop precursor molecule, and, to a lesser degree, based on evolutionary conservation of small RNAs from other vertebrates. Properly curated miRNA databases are an important resource for investigators interested in miRNA biology, diagnostics, and therapeutics.

Antagonism of microRNA function in zebrafish embryos by using locked nucleic acid enzymes (LNAzymes).

MicroRNAs (miRNAs) have crucial functions in many cellular processes, such as differentiation, proliferation and apoptosis; aberrant expression of miRNAs has been linked to human diseases, including cancer. Tools that allow specific and efficient knockdown of miRNAs would be of immense importance for exploring miRNA function. Zebrafish serves as an excellent vertebrate model system to understand the functions of miRNAs involved in a variety of biological processes. We designed and employed a strategy based on locked nucleic acid enzymes (LNAzymes) for in vivo knockdown of miRNA in zebrafish embryos. We demonstrate that LNAzyme can efficiently knockdown miRNAs with minimal toxicity to the zebrafish embryos.