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.

Inherited human c-Rel deficiency disrupts myeloid and lymphoid immunity to multiple infectious agents.

We studied a child with severe viral, bacterial, fungal, and parasitic diseases, who was homozygous for a loss-of-function mutation of REL, encoding c-Rel, which is selectively expressed in lymphoid and myeloid cells. The patient had low frequencies of NK, effector memory cells reexpressing CD45RA (Temra) CD8+ T cells, memory CD4+ T cells, including Th1 and Th1*, Tregs, and memory B cells, whereas the counts and proportions of other leukocyte subsets were normal. Functional deficits of myeloid cells included the abolition of IL-12 and IL-23 production by conventional DC1s (cDC1s) and monocytes, but not cDC2s. c-Rel was also required for induction of CD86 expression on, and thus antigen-presenting cell function of, cDCs. Functional deficits of lymphoid cells included reduced IL-2 production by naive T cells, correlating with low proliferation and survival rates and poor production of Th1, Th2, and Th17 cytokines by memory CD4+ T cells. In naive CD4+ T cells, c-Rel is dispensable for early IL2 induction but contributes to later phases of IL2 expression. The patient's naive B cells displayed impaired MYC and BCL2L1 induction, compromising B cell survival and proliferation and preventing their differentiation into Ig-secreting plasmablasts. Inherited c-Rel deficiency disrupts the development and function of multiple myeloid and lymphoid cells, compromising innate and adaptive immunity to multiple infectious agents.

Deep sequencing of DNA from urine of kidney allograft recipients to estimate donor/recipient-specific DNA fractions.

Kidney transplantation is the treatment of choice for patients with end-stage kidney failure, but transplanted allograft could be affected by viral and bacterial infections and by immune rejection. The standard test for the diagnosis of acute pathologies in kidney transplants is kidney biopsy. However, noninvasive tests would be desirable. Various methods using different techniques have been developed by the transplantation community. But these methods require improvements. We present here a cost-effective method for kidney rejection diagnosis that estimates donor/recipient-specific DNA fraction in recipient urine by sequencing urinary cell DNA. We hypothesized that in the no-pathology stage, the largest tissue types present in recipient urine are donor kidney cells, and in case of rejection, a larger number of recipient immune cells would be observed. Extensive in-silico simulation was used to tune the sequencing parameters: number of variants and depth of coverage. Sequencing of DNA mixture from 2 healthy individuals showed the method is highly predictive (maximum error < 0.04). We then demonstrated the insignificant impact of familial relationship and ethnicity using an in-house and public database. Lastly, we performed deep DNA sequencing of urinary cell pellets from 32 biopsy-matched samples representing two pathology groups: acute rejection (AR, 11 samples) and acute tubular injury (ATI, 12 samples) and 9 samples with no pathology. We found a significant association between the donor/recipient-specific DNA fraction in the two pathology groups compared to no pathology (P = 0.0064 for AR and P = 0.026 for ATI). We conclude that deep DNA sequencing of urinary cells from kidney allograft recipients offers a noninvasive means of diagnosing acute pathologies in the human kidney allograft.

Identification of genetic variants controlling RNA editing and their effect on RNA structure stabilization.

Post-transcriptional modification of RNA (RNA editing, RNAe) results in differences between the RNA transcript and the genomic DNA sequence (RDD). Enzymatic modification of adenosine to inosine (A2I) by ADAR is the most studied type of RNAe. However, few genetic association studies with A2I RNAe events have been conducted. Some studies have analyzed the inter-population RNAe-QTL diversity in humans, but the sample size of these studies was limited. Other types of RNA and DNA differences have been reported but are largely understudied. Here, we report a comprehensive analysis of all types of RDD, based on two independent datasets. We found that A2I was by far the most observed type of RDD. Moreover, manual curation suggests that A2I is likely the only enzymatically driven RNAe type observed in blood derived DNA, all other non-A2I RDD could either be attributed to sequencing and processing artifacts, or are a result of somatic DNA rearrangements. We then conducted an in-cis genetic association study and identified 472 genetic associations (RNAe-QTL), that were replicated in both datasets. We confirm the potential effect of the RNAe-QTL on RNA structure by showing that allele specific RNAe occurs in heterozygotes. Although the generally assumed function of RNAe is to destabilize double stranded RNA structure, we found clear evidence for the potential additional involvement of RNAe in maintaining RNA hairpin that has been altered by the RNAe-QTL. Our study confirms, in two independent datasets, the potential role of RNAe in maintaining RNA structure in the presence of genetic variation.

A recessive form of hyper-IgE syndrome by disruption of ZNF341-dependent STAT3 transcription and activity.

Heterozygosity for human signal transducer and activator of transcription 3 (STAT3) dominant-negative (DN) mutations underlies an autosomal dominant form of hyper-immunoglobulin E syndrome (HIES). We describe patients with an autosomal recessive form of HIES due to loss-of-function mutations of a previously uncharacterized gene, ZNF341 ZNF341 is a transcription factor that resides in the nucleus, where it binds a specific DNA motif present in various genes, including the STAT3 promoter. The patients' cells have low basal levels of STAT3 mRNA and protein. The autoinduction of STAT3 production, activation, and function by STAT3-activating cytokines is strongly impaired. Like patients with STAT3 DN mutations, ZNF341-deficient patients lack T helper 17 (TH17) cells, have an excess of TH2 cells, and have low memory B cells due to the tight dependence of STAT3 activity on ZNF341 in lymphocytes. Their milder extra-hematopoietic manifestations and stronger inflammatory responses reflect the lower ZNF341 dependence of STAT3 activity in other cell types. Human ZNF341 is essential for the STAT3 transcription-dependent autoinduction and sustained activity of STAT3.

Single nucleotide variant counts computed from RNA sequencing and cellular traffic into human kidney allografts.

Advances in bioinformatics allow identification of single nucleotide polymorphisms (variants) from RNA sequence data. In an allograft biopsy, 2 genomes contribute to the RNA pool, 1 from the donor organ and the other from the infiltrating recipient's cells. We hypothesize that imbalances in genetic variants of RNA sequence data of kidney allograft biopsies provide an objective measure of cellular infiltration of the allograft. We performed mRNA sequencing of 40 kidney allograft biopsies, selected to represent a comprehensive range of diagnostic categories. We analyzed the sequencing reads of these biopsies and of 462 lymphoblastoid cell lines from the 1000 Genomes Project, for RNA variants. The ratio of heterozygous to nonreference genome homozygous variants (Het/Hom ratio) on all autosomes was determined for each sample, and the estimation of stromal and immune cells in malignant tumors using expression data (ESTIMATE) score was computed as a complementary estimate of the degree of cellular infiltration into biopsies. The Het/Hom ratios (P = .02) and the ESTIMATE scores (P < .001) were associated with the biopsy diagnosis. Both measures correlated significantly (r = .67, P < .0001), even though the Het/Hom ratio is based on mRNA sequence variation, while the ESTIMATE score uses mRNA expression. Het/Hom ratio and the ESTIMATE score may offer unbiased and quantitative parameters for characterizing cellular traffic into human kidney allografts.

Recurrent rhinovirus infections in a child with inherited MDA5 deficiency.

MDA5 is a cytosolic sensor of double-stranded RNA (ds)RNA including viral byproducts and intermediates. We studied a child with life-threatening, recurrent respiratory tract infections, caused by viruses including human rhinovirus (HRV), influenza virus, and respiratory syncytial virus (RSV). We identified in her a homozygous missense mutation in IFIH1 that encodes MDA5. Mutant MDA5 was expressed but did not recognize the synthetic MDA5 agonist/(ds)RNA mimic polyinosinic-polycytidylic acid. When overexpressed, mutant MDA5 failed to drive luciferase activity from the IFNB1 promoter or promoters containing ISRE or NF-κB sequence motifs. In respiratory epithelial cells or fibroblasts, wild-type but not knockdown of MDA5 restricted HRV infection while increasing IFN-stimulated gene expression and IFN-ß/λ. However, wild-type MDA5 did not restrict influenza virus or RSV replication. Moreover, nasal epithelial cells from the patient, or fibroblasts gene-edited to express mutant MDA5, showed increased replication of HRV but not influenza or RSV. Thus, human MDA5 deficiency is a novel inborn error of innate and/or intrinsic immunity that causes impaired (ds)RNA sensing, reduced IFN induction, and susceptibility to the common cold virus.

Human Adaptive Immunity Rescues an Inborn Error of Innate Immunity.

The molecular basis of the incomplete penetrance of monogenic disorders is unclear. We describe here eight related individuals with autosomal recessive TIRAP deficiency. Life-threatening staphylococcal disease occurred during childhood in the proband, but not in the other seven homozygotes. Responses to all Toll-like receptor 1/2 (TLR1/2), TLR2/6, and TLR4 agonists were impaired in the fibroblasts and leukocytes of all TIRAP-deficient individuals. However, the whole-blood response to the TLR2/6 agonist staphylococcal lipoteichoic acid (LTA) was abolished only in the index case individual, the only family member lacking LTA-specific antibodies (Abs). This defective response was reversed in the patient, but not in interleukin-1 receptor-associated kinase 4 (IRAK-4)-deficient individuals, by anti-LTA monoclonal antibody (mAb). Anti-LTA mAb also rescued the macrophage response in mice lacking TIRAP, but not TLR2 or MyD88. Thus, acquired anti-LTA Abs rescue TLR2-dependent immunity to staphylococcal LTA in individuals with inherited TIRAP deficiency, accounting for incomplete penetrance. Combined TIRAP and anti-LTA Ab deficiencies underlie staphylococcal disease in this patient.

Human HOIP and LUBAC deficiency underlies autoinflammation, immunodeficiency, amylopectinosis, and lymphangiectasia.

Inherited, complete deficiency of human HOIL-1, a component of the linear ubiquitination chain assembly complex (LUBAC), underlies autoinflammation, infections, and amylopectinosis. We report the clinical description and molecular analysis of a novel inherited disorder of the human LUBAC complex. A patient with multiorgan autoinflammation, combined immunodeficiency, subclinical amylopectinosis, and systemic lymphangiectasia, is homozygous for a mutation in HOIP, the gene encoding the catalytic component of LUBAC. The missense allele (L72P, in the PUB domain) is at least severely hypomorphic, as it impairs HOIP expression and destabilizes the whole LUBAC complex. Linear ubiquitination and NF-κB activation are impaired in the patient's fibroblasts stimulated by IL-1ß or TNF. In contrast, the patient's monocytes respond to IL-1ß more vigorously than control monocytes. However, the activation and differentiation of the patient's B cells are impaired in response to CD40 engagement. These cellular and clinical phenotypes largely overlap those of HOIL-1-deficient patients. Clinical differences between HOIL-1- and HOIP-mutated patients may result from differences between the mutations, the loci, or other factors. Our findings show that human HOIP is essential for the assembly and function of LUBAC and for various processes governing inflammation and immunity in both hematopoietic and nonhematopoietic cells.

Inherited IL-17RC deficiency in patients with chronic mucocutaneous candidiasis.

Chronic mucocutaneous candidiasis (CMC) is characterized by recurrent or persistent infections of the skin, nail, oral, and genital mucosae with Candida species, mainly C. albicans. Autosomal-recessive (AR) IL-17RA and ACT1 deficiencies and autosomal-dominant IL-17F deficiency, each reported in a single kindred, underlie CMC in otherwise healthy patients. We report three patients from unrelated kindreds, aged 8, 12, and 37 yr with isolated CMC, who display AR IL-17RC deficiency. The patients are homozygous for different nonsense alleles that prevent the expression of IL-17RC on the cell surface. The defect is complete, abolishing cellular responses to IL-17A and IL-17F homo- and heterodimers. However, in contrast to what is observed for the IL-17RA- and ACT1-deficient patients tested, the response to IL-17E (IL-25) is maintained in these IL-17RC-deficient patients. These experiments of nature indicate that human IL-17RC is essential for mucocutaneous immunity to C. albicans but is otherwise largely redundant.

An ACT1 mutation selectively abolishes interleukin-17 responses in humans with chronic mucocutaneous candidiasis.

Patients with inborn errors of interleukin-17F (IL-17F) or IL-17RA display chronic mucocutaneous candidiasis (CMC). We report a biallelic missense mutation (T536I) in the adaptor molecule ACT1 in two siblings with CMC. The mutation, located in the SEFIR domain, abolished the homotypic interaction of ACT1 with IL-17 receptors, with no effect on homodimerization. The patients' fibroblasts failed to respond to IL-17A and IL-17F, and their T cells to IL-17E. By contrast, healthy individuals homozygous for the common variant D10N, located in the ACT1 tumor necrosis factor receptor-associated factor-interacting domain and previously associated with psoriasis, had impaired, but not abolished, responses to IL-17 cytokines. SEFIR-independent interactions of ACT1 with other proteins, such as CD40, heat shock protein 70 (HSP70) and HSP90, were not affected by the T536I mutation. Overall, human IL-17A and IL-17F depend on ACT1 to mediate protective mucocutaneous immunity. Moreover, other ACT1-dependent IL-17 cytokines seem to be largely redundant in host defense.