The dynamic epigenetic regulation of the inactive X chromosome in healthy human B cells is dysregulated in lupus patients.

Systemic lupus erythematous (SLE) is a female-predominant disease characterized by autoimmune B cells and pathogenic autoantibody production. Individuals with two or more X chromosomes are at increased risk for SLE, suggesting that X-linked genes contribute to the observed sex bias of this disease. To normalize X-linked gene expression between sexes, one X in female cells is randomly selected for transcriptional silencing through X-chromosome inactivation (XCI), resulting in allele-specific enrichment of epigenetic modifications, including histone methylation and the long noncoding RNA XIST/Xist on the inactive X (Xi). As we have previously shown that epigenetic regulation of the Xi in female lymphocytes from mice is unexpectedly dynamic, we used RNA fluorescence in situ hybridization and immunofluorescence to profile epigenetic features of the Xi at the single-cell level in human B cell subsets from pediatric and adult SLE patients and healthy controls. Our data reveal that abnormal XCI maintenance in B cells is a feature of SLE. Using single-cell and bulk-cell RNA sequencing datasets, we found that X-linked immunity genes escape XCI in specific healthy human B cell subsets and that human SLE B cells exhibit aberrant expression of X-linked genes and XIST RNA interactome genes. Our data reveal that mislocalized XIST RNA, coupled with a dramatic reduction in heterochromatic modifications at the Xi in SLE, predispose for aberrant X-linked gene expression from the Xi, thus defining a genetic and epigenetic pathway that affects X-linked gene expression in human SLE B cells and likely contributes to the female bias in SLE.

Myeloid-derived miR-6236 potentiates adipocyte insulin signaling and prevents hyperglycemia during obesity.

Adipose tissue macrophages (ATMs) influence obesity-associated metabolic dysfunction, but the mechanisms by which they do so are not well understood. We show that miR-6236 is a bona fide miRNA that is secreted by ATMs during obesity. Global or myeloid cell-specific deletion of miR-6236 aggravates obesity-associated adipose tissue insulin resistance, hyperglycemia, hyperinsulinemia, and hyperlipidemia. miR-6236 augments adipocyte insulin sensitivity by inhibiting translation of negative regulators of insulin signaling, including PTEN. The human genome harbors a miR-6236 homolog that is highly expressed in the serum and adipose tissue of obese people. hsa-MIR-6236 expression negatively correlates with hyperglycemia and glucose intolerance, and positively correlates with insulin sensitivity. Together, our findings establish miR-6236 as an ATM-secreted miRNA that potentiates adipocyte insulin signaling and protects against metabolic dysfunction during obesity.

The role of extracellular vesicle-derived miRNAs in adipose tissue function and metabolic health.

Extracellular vesicles (EVs) are nanometer size lipid particles that are released from virtually every cell type. Recent studies have shown that miRNAs carried by EVs play important roles in intercellular and interorgan communication. In the context of obesity and insulin resistance, EV-derived miRNAs functionally bridge major metabolic organs, including the adipose tissue, skeletal muscle, liver, and pancreas, to regulate insulin secretion and signaling. As a result, many of these EV-derived miRNAs have been proposed as potential disease biomarkers and/or therapeutic agents. However, the field's knowledge of EV miRNA-mediated regulation of mammalian metabolism is still in its infancy. Here, we review the evidence indicating that EV-derived miRNAs provide cell-to-cell and organ-to-organ communication to support metabolic health, highlight the potential medical relevance of these discoveries, and discuss the most important knowledge gaps and future directions for this field.

Altered adipose tissue macrophage populations in people with HIV on integrase inhibitor-containing antiretroviral therapy.

OBJECTIVE: Antiretroviral therapy (ART) extends the life of people with HIV (PWH), but these individuals are at increased risk for obesity, dyslipidemia, diabetes, and cardiovascular disease. These comorbidities may be a consequence of HIV-related chronic inflammation and/or adverse effects of ART on tissue regulatory adipose tissue macrophages (ATMs). We sought to determine the effects of HIV/ART on metabolically beneficial ATM populations and functions. DESIGN: We examined subcutaneous ATMs from PWH on integrase inhibitor-containing ART ( n  = 5) and uninfected persons ( n  = 9). We complemented these studies with ex vivo and in vitro analyses of peripheral blood mononuclear cell (PBMC) and murine macrophage lipid metabolism and fatty acid oxidation gene expression. METHODS: ATM populations were examined by flow cytometry. Macrophage lipid metabolism and fatty acid oxidation gene expression were examined by Seahorse assay and quantitative PCR. RESULTS: Adipose tissue from PWH had reduced populations of metabolically activated CD9 + ATMs compared to that of uninfected controls ( P  < 0.001). PBMCs of PWH had lower fatty acid metabolism compared to those of uninfected controls ( P  < 0.01). Analysis of murine macrophages revealed that dolutegravir reduced lipid metabolism ( P  < 0.001) and increased expression of the fatty acid beta-oxidation enzyme enoyl-CoA hydratase, short chain 1 ( P  < 0.05). CONCLUSIONS: We report the loss of metabolically beneficial ATM populations in PWH on ART, altered fatty acid metabolism of blood immune cells, and evidence that dolutegravir alters macrophage fatty acid metabolism. Future studies should examine direct or indirect effects and mechanisms of dolutegravir, and other integrase inhibitors and ART classes, on fatty acid beta-oxidation.

Altered X-chromosome inactivation in T cells may promote sex-biased autoimmune diseases.

Systemic lupus erythematosus (SLE) is an autoimmune disorder that predominantly affects women and is driven by autoreactive T cell-mediated inflammation. It is known that individuals with multiple X-chromosomes are at increased risk for developing SLE; however, the mechanisms underlying this genetic basis are unclear. Here, we use single cell imaging to determine the epigenetic features of the inactive X (Xi) in developing thymocytes, mature T cell subsets, and T cells from SLE patients and mice. We show that Xist RNA and heterochromatin modifications transiently reappear at the Xi and are missing in mature single positive T cells. Activation of mature T cells restores Xist RNA and heterochromatin marks simultaneously back to the Xi. Notably, X-chromosome inactivation (XCI) maintenance is altered in T cells of SLE patients and late-stage-disease NZB/W F1 female mice, and we show that X-linked genes are abnormally upregulated in SLE patient T cells. SLE T cells also have altered expression of XIST RNA interactome genes, accounting for perturbations of Xi epigenetic features. Thus, abnormal XCI maintenance is a feature of SLE disease, and we propose that Xist RNA localization at the Xi could be an important factor for maintaining dosage compensation of X-linked genes in T cells.

Crosstalk among lncRNAs, microRNAs and mRNAs in the muscle 'degradome' of rainbow trout.

In fish, protein-coding and noncoding genes involved in muscle atrophy are not fully characterized. In this study, we characterized coding and noncoding genes involved in gonadogenesis-associated muscle atrophy, and investigated the potential functional interplay between these genes. Using RNA-Seq, we compared expression pattern of mRNAs, long noncoding RNAs (lncRNAs) and microRNAs of atrophying skeletal muscle from gravid females and control skeletal muscle from age-matched sterile individuals. A total of 852 mRNAs, 1,160 lncRNAs and 28 microRNAs were differentially expressed (DE) between the two groups. Muscle atrophy appears to be mediated by many genes encoding ubiquitin-proteasome system, autophagy related proteases, lysosomal proteases and transcription factors. Transcripts encoding atrogin-1 and mir-29 showed exceptional high expression in atrophying muscle, suggesting an important role in bulk muscle proteolysis. DE genes were co-localized in the genome with strong expression correlation, and they exhibited extensive 'lncRNA-mRNA', 'lncRNA-microRNA', 'mRNA-microRNA' and 'lncRNA-protein' physical interactions. DE genes exhibiting potential functional interactions comprised the highly correlated 'lncRNA-mRNA-microRNA' gene network described as 'degradome'. This study pinpoints extensive coding and noncoding RNA interactions during muscle atrophy in fish, and provides valuable resources for future mechanistic studies.

RNA-Seq reveals MicroRNA expression signature and genetic polymorphism associated with growth and muscle quality traits in rainbow trout.

The role of microRNA expression and genetic variation in microRNA-binding sites of target genes on growth and muscle quality traits is poorly characterized. We used RNA-Seq approach to investigate their importance on 5 growth and muscle quality traits: whole body weight (WBW), muscle yield, muscle crude-fat content, muscle shear force and whiteness. Phenotypic data were collected from 471 fish, representing 98 families (~5 fish/family) from a growth-selected line. Muscle microRNAs and mRNAs were sequenced from 22 families showing divergent phenotypes. Ninety microRNAs showed differential expression between families with divergent phenotypes, and their expression was strongly associated with variation in phenotypes. A total of 204 single nucleotide polymorphisms (SNPs) present in 3' UTR of target genes either destroyed or created novel illegitimate microRNA target sites; of them, 78 SNPs explained significant variation in the aforementioned 5 muscle traits. Majority of the phenotype-associated SNPs were present in microRNA-binding sites of genes involved in energy metabolism and muscle structure. These findings suggest that variation in microRNA expression and/or sequence variation in microRNA binding sites in target genes play an important role in mediating differences in fish growth and muscle quality phenotypes.

Genome-Wide Discovery of Long Non-Coding RNAs in Rainbow Trout.

The ENCODE project revealed that ~70% of the human genome is transcribed. While only 1-2% of the RNAs encode for proteins, the rest are non-coding RNAs. Long non-coding RNAs (lncRNAs) form a diverse class of non-coding RNAs that are longer than 200 nt. Emerging evidence indicates that lncRNAs play critical roles in various cellular processes including regulation of gene expression. LncRNAs show low levels of gene expression and sequence conservation, which make their computational identification in genomes difficult. In this study, more than two billion Illumina sequence reads were mapped to the genome reference using the TopHat and Cufflinks software. Transcripts shorter than 200 nt, with more than 83-100 amino acids ORF, or with significant homologies to the NCBI nr-protein database were removed. In addition, a computational pipeline was used to filter the remaining transcripts based on a protein-coding-score test. Depending on the filtering stringency conditions, between 31,195 and 54,503 lncRNAs were identified, with only 421 matching known lncRNAs in other species. A digital gene expression atlas revealed 2,935 tissue-specific and 3,269 ubiquitously-expressed lncRNAs. This study annotates the lncRNA rainbow trout genome and provides a valuable resource for functional genomics research in salmonids.

Differential expression of long non-coding RNAs in three genetic lines of rainbow trout in response to infection with Flavobacterium psychrophilum.

Bacterial cold-water disease caused by Flavobacterium psychrophilum is one of the major causes of mortality of salmonids. Three genetic lines of rainbow trout designated as ARS-Fp-R (resistant), ARS-Fp-C (control) and ARS-Fp-S (susceptible) have significant differences in survival rate following F. psychrophilum infection. Previous study identified transcriptome differences of immune-relevant protein-coding genes at basal and post infection levels among these genetic lines. Using RNA-Seq approach, we quantified differentially expressed (DE) long non-coding RNAs (lncRNAs) in response to F. psychrophilum challenge in these genetic lines. Pairwise comparison between genetic lines and different infection statuses identified 556 DE lncRNAs. A positive correlation existed between the number of the differentially regulated lncRNAs and that of the protein-coding genes. Several lncRNAs showed strong positive and negative expression correlation with their overlapped, neighboring and distant immune related protein-coding genes including complement components, cytokines, chemokines and several signaling molecules involved in immunity. The correlated expressions and genome-wide co-localization suggested that some lncRNAs may be involved in regulating immune-relevant protein-coding genes. This study provides the first evidence of lncRNA-mediated regulation of the anti-bacterial immune response in a commercially important aquaculture species and will likely help developing new genetic markers for rainbow trout disease resistance.

Transcriptome assembly, gene annotation and tissue gene expression atlas of the rainbow trout.

Efforts to obtain a comprehensive genome sequence for rainbow trout are ongoing and will be complemented by transcriptome information that will enhance genome assembly and annotation. Previously, transcriptome reference sequences were reported using data from different sources. Although the previous work added a great wealth of sequences, a complete and well-annotated transcriptome is still needed. In addition, gene expression in different tissues was not completely addressed in the previous studies. In this study, non-normalized cDNA libraries were sequenced from 13 different tissues of a single doubled haploid rainbow trout from the same source used for the rainbow trout genome sequence. A total of ~1.167 billion paired-end reads were de novo assembled using the Trinity RNA-Seq assembler yielding 474,524 contigs > 500 base-pairs. Of them, 287,593 had homologies to the NCBI non-redundant protein database. The longest contig of each cluster was selected as a reference, yielding 44,990 representative contigs. A total of 4,146 contigs (9.2%), including 710 full-length sequences, did not match any mRNA sequences in the current rainbow trout genome reference. Mapping reads to the reference genome identified an additional 11,843 transcripts not annotated in the genome. A digital gene expression atlas revealed 7,678 housekeeping and 4,021 tissue-specific genes. Expression of about 16,000-32,000 genes (35-71% of the identified genes) accounted for basic and specialized functions of each tissue. White muscle and stomach had the least complex transcriptomes, with high percentages of their total mRNA contributed by a small number of genes. Brain, testis and intestine, in contrast, had complex transcriptomes, with a large numbers of genes involved in their expression patterns. This study provides comprehensive de novo transcriptome information that is suitable for functional and comparative genomics studies in rainbow trout, including annotation of the genome.

Whole-body transcriptome of selectively bred, resistant-, control-, and susceptible-line rainbow trout following experimental challenge with Flavobacterium psychrophilum.

Genetic improvement for enhanced disease resistance in fish is an increasingly utilized approach to mitigate endemic infectious disease in aquaculture. In domesticated salmonid populations, large phenotypic variation in disease resistance has been identified but the genetic basis for altered responsiveness remains unclear. We previously reported three generations of selection and phenotypic validation of a bacterial cold water disease (BCWD) resistant line of rainbow trout, designated ARS-Fp-R. This line has higher survival after infection by either standardized laboratory challenge or natural challenge as compared to two reference lines, designated ARS-Fp-C (control) and ARS-Fp-S (susceptible). In this study, we utilized 1.1 g fry from the three genetic lines and performed RNA-seq to measure transcript abundance from the whole body of naive and Flavobacterium psychrophilum infected fish at day 1 (early time-point) and at day 5 post-challenge (onset of mortality). Sequences from 24 libraries were mapped onto the rainbow trout genome reference transcriptome of 46,585 predicted protein coding mRNAs that included 2633 putative immune-relevant gene transcripts. A total of 1884 genes (4.0% genome) exhibited differential transcript abundance between infected and mock-challenged fish (FDR < 0.05) that included chemokines, complement components, tnf receptor superfamily members, interleukins, nod-like receptor family members, and genes involved in metabolism and wound healing. The largest number of differentially expressed genes occurred on day 5 post-infection between naive and challenged ARS-Fp-S line fish correlating with high bacterial load. After excluding the effect of infection, we identified 21 differentially expressed genes between the three genetic lines. In summary, these data indicate global transcriptome differences between genetic lines of naive animals as well as differentially regulated transcriptional responses to infection.