Multisystem inflammatory syndrome in children (MIS-C) is a severe complication of SARS-CoV-2 infections that occurs in the pediatric population. We sought to characterize T cell responses in MIS-C ...compared to COVID-19 and pediatric hyperinflammatory syndromes.
MIS-C was distinct from COVID-19 and hyperinflammatory syndromes due to an expansion of T cells expressing TRBV11–2 that was not associated with HLA genotype. Children diagnosed with MIS-C, but who were negative for SARS-CoV-2 by PCR and serology, did not display Vβ skewing. There was no difference in the proportion of T cells that became activated after stimulation with SARS-CoV-2 peptides in children with MIS-C compared to convalescent COVID-19. The frequency of SARS-CoV-2-specific TCRs and the antigens recognized by these TCRs were comparable in MIS-C and COVID-19.
Expansion of Vβ11–2+ T cells was a specific biomarker of MIS-C patients with laboratory confirmed SARS-CoV-2 infections. Children with MIS-C had robust antigen-specific T cell responses to SARS-CoV-2.
Repetitive exposure to antigen in chronic infection and cancer drives T cell exhaustion, limiting adaptive immunity. In contrast, aberrant, sustained T cell responses can persist over decades in ...human allergic disease. To understand these divergent outcomes, we employed bioinformatic, immunophenotyping and functional approaches with human diseased tissues, identifying an abundant population of type 2 helper T (T
2) cells with co-expression of TCF7 and LEF1, and features of chronic activation. These cells, which we termed T
2-multipotent progenitors (T
2-MPP) could self-renew and differentiate into cytokine-producing effector cells, regulatory T (T
) cells and follicular helper T (T
) cells. Single-cell T-cell-receptor lineage tracing confirmed lineage relationships between T
2-MPP, T
2 effectors, T
cells and T
cells. T
2-MPP persisted despite in vivo IL-4 receptor blockade, while thymic stromal lymphopoietin (TSLP) drove selective expansion of progenitor cells and rendered them insensitive to glucocorticoid-induced apoptosis in vitro. Together, our data identify T
2-MPP as an aberrant T cell population with the potential to sustain type 2 inflammation and support the paradigm that chronic T cell responses can be coordinated over time by progenitor cells.
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.
Secondary lymphedema occurs in up to 20% of patients after lymphadenectomy performed for the surgical management of tumors involving the breast, prostate, uterus, and skin. Patients develop ...progressive edema of the affected extremity due to retention of protein-rich lymphatic fluid. Despite compression therapy, patients progress to chronic lymphedema in which noncompressible fibrosis and adipose tissue are deposited within the extremity. The presence of fibrosis led to our hypothesis that rosiglitazone, a PPARγ agonist that inhibits fibrosis, would reduce fibrosis in a mouse model of secondary lymphedema after hind limb lymphadenectomy. In vivo, rosiglitazone reduced fibrosis in the hind limb after lymphadenectomy. Our findings verified that rosiglitazone reestablished the adipogenic features of TGF-β1-treated mesenchymal cells in vitro. Despite this, rosiglitazone led to a reduction in adipose tissue deposition. Single-cell RNA-Seq data obtained from human tissues and flow cytometric and histological evaluation of mouse tissues demonstrated increased presence of PDGFRα+ cells in lymphedema; human tissue analysis verified these cells have the capacity for adipogenic and fibrogenic differentiation. Upon treatment with rosiglitazone, we noted a reduction in the overall quantity of PDGFRα+ cells and LipidTOX+ cells. Our findings provide a framework for treating secondary lymphedema as a condition of fibrosis and adipose tissue deposition, both of which, paradoxically, can be prevented with a pro-adipogenic agent.
DNA sequence variation has been associated with quantitative changes in molecular phenotypes such as gene expression, but its impact on chromatin states is poorly characterized. To understand the ...interplay between chromatin and genetic control of gene regulation, we quantified allelic variability in transcription factor binding, histone modifications, and gene expression within humans. We found abundant allelic specificity in chromatin and extensive local, short-range, and long-range allelic coordination among the studied molecular phenotypes. We observed genetic influence on most of these phenotypes, with histone modifications exhibiting strong context-dependent behavior. Our results implicate transcription factors as primary mediators of sequence-specific regulation of gene expression programs, with histone modifications frequently reflecting the primary regulatory event.
RNA sequencing (RNA-seq) is the current gold-standard method to quantify gene expression for expression quantitative trait locus (eQTL) studies. However, a potential caveat in these studies is that ...RNA-seq reads carrying the non-reference allele of variant loci can have lower probability to map correctly to the reference genome, which could bias gene quantifications and cause false positive eQTL associations. In this study, we analyze the effect of this allelic mapping bias in eQTL discovery.
We simulate RNA-seq read mapping over 9.5 M common SNPs and indels, with 15.6% of variants showing biased mapping rate for reference versus non-reference reads. However, removing potentially biased RNA-seq reads from an eQTL dataset of 185 individuals has a very small effect on gene and exon quantifications and eQTL discovery. We detect only a handful of likely false positive eQTLs, and overall eQTL SNPs show no significant enrichment for high mapping bias.
Our results suggest that RNA-seq quantifications are generally robust against allelic mapping bias, and that this does not have a severe effect on eQTL discovery. Nevertheless, we provide our catalog of putatively biased loci to allow better controlling for mapping bias to obtain more accurate results in future RNA-seq studies.
was among the first loci shown to confer risk for inflammatory arthritis in the absence of an associated coding variant, but its genetic mechanism remains undefined. Using Immunochip data from 3,939 ...patients with juvenile idiopathic arthritis (JIA) and 14,412 control individuals, we identified 132 plausible common non-coding variants, reduced serially by single-nucleotide polymorphism sequencing (SNP-seq), electrophoretic mobility shift, and luciferase studies to the single variant rs7034653 in the third intron of
. Genetically manipulated experimental cells and primary monocytes from genotyped donors establish that the risk G allele reduces binding of Fos-related antigen 2 (FRA2), encoded by
, resulting in reduced TRAF1 expression and enhanced tumor necrosis factor (TNF) production. Conditioning on this JIA variant eliminated attributable risk for rheumatoid arthritis, implicating a mechanism shared across the arthritis spectrum. These findings reveal that rs7034653, FRA2, and TRAF1 mediate a pathway through which a non-coding functional variant drives risk of inflammatory arthritis in children and adults.
Identical sequences with a minimal length of about 300 base pairs (bp) have been involved in the generation of various meiotic/mitotic genomic rearrangements through non-allelic homologous ...recombination (NAHR) events. Genomic disorders and structural variation, together with gene remodelling processes have been associated with many of these rearrangements. Based on these observations, we identified and integrated all the 100% identical repeats of at least 300 bp in the NCBI version 36.2 human genome reference assembly into non-overlapping regions, thus defining the Identical Repeated Backbone (IRB) of the reference human genome.
The IRB sequences are distributed all over the genome in 66,600 regions, which correspond to approximately 2% of the total NCBI human genome reference assembly. Important structural and functional elements such as common repeats, segmental duplications, and genes are contained in the IRB. About 80% of the IRB bp overlap with known copy-number variants (CNVs). By analyzing the genes embedded in the IRB, we were able to detect some identical genes not previously included in the Ensembl release 50 annotation of human genes. In addition, we found evidence of IRB gene copy-number polymorphisms in raw sequence reads of two diploid sequenced genomes.
In general, the IRB offers new insight into the complex organization of the identical repeated sequences of the human genome. It provides an accurate map of potential NAHR sites which could be used in targeting the study of novel CNVs, predicting DNA copy-number variation in newly sequenced genomes, and improve genome annotation.
Population-scale genome sequencing allows the characterization of functional effects of a broad spectrum of genetic variants underlying human phenotypic variation. Here, we investigate the influence ...of rare and common genetic variants on gene expression patterns, using variants identified from sequencing data from the 1000 genomes project in an African and European population sample and gene expression data from lymphoblastoid cell lines. We detect comparable numbers of expression quantitative trait loci (eQTLs) when compared to genotypes obtained from HapMap 3, but as many as 80% of the top expression quantitative trait variants (eQTVs) discovered from 1000 genomes data are novel. The properties of the newly discovered variants suggest that mapping common causal regulatory variants is challenging even with full resequencing data; however, we observe significant enrichment of regulatory effects in splice-site and nonsense variants. Using RNA sequencing data, we show that 46.2% of nonsynonymous variants are differentially expressed in at least one individual in our sample, creating widespread potential for interactions between functional protein-coding and regulatory variants. We also use allele-specific expression to identify putative rare causal regulatory variants. Furthermore, we demonstrate that outlier expression values can be due to rare variant effects, and we approximate the number of such effects harboured in an individual by effect size. Our results demonstrate that integration of genomic and RNA sequencing analyses allows for the joint assessment of genome sequence and genome function.