IgG is an aging factor that drives adipose tissue fibrosis and metabolic decline.

Aging is underpinned by pronounced metabolic decline; however, the drivers remain obscure. Here, we report that IgG accumulates during aging, particularly in white adipose tissue (WAT), to impair adipose tissue function and metabolic health. Caloric restriction (CR) decreases IgG accumulation in WAT, whereas replenishing IgG counteracts CR's metabolic benefits. IgG activates macrophages via Ras signaling and consequently induces fibrosis in WAT through the TGF-ß/SMAD pathway. Consistently, B cell null mice are protected from aging-associated WAT fibrosis, inflammation, and insulin resistance, unless exposed to IgG. Conditional ablation of the IgG recycling receptor, neonatal Fc receptor (FcRn), in macrophages prevents IgG accumulation in aging, resulting in prolonged healthspan and lifespan. Further, targeting FcRn by antisense oligonucleotide restores WAT integrity and metabolic health in aged mice. These findings pinpoint IgG as a hidden culprit in aging and enlighten a novel strategy to rejuvenate metabolic health.

Deciphering the decline of metabolic elasticity in aging and obesity.

Organisms must adapt to fluctuating nutrient availability to maintain energy homeostasis. Here, we term the capacity for such adaptation and restoration "metabolic elasticity" and model it through ad libitum-fasting-refeeding cycles. Metabolic elasticity is achieved by coordinate versatility in gene expression, which we call "gene elasticity." We have developed the gene elasticity score as a systematic method to quantify the elasticity of the transcriptome across metabolically active tissues in mice and non-human primates. Genes involved in lipid and carbohydrate metabolism show high gene elasticity, and their elasticity declines with age, particularly with PPARγ dysregulation in adipose tissue. Synchronizing PPARγ activity with nutrient conditions through feeding-timed agonism optimizes their metabolic benefits and safety. We further broaden the conceptual scope of metabolic and gene elasticity to dietary challenges, revealing declines in diet-induced obesity similar to those in aging. Altogether, our findings provide a dynamic perspective on the dysmetabolic consequences of aging and obesity.

Multidimensional conservation analysis decodes the expression of conserved long noncoding RNAs.

Although long noncoding RNAs (lncRNAs) experience weaker evolutionary constraints and exhibit lower sequence conservation than coding genes, they can still conserve their features in various aspects. Here, we used multiple approaches to systemically evaluate the conservation between human and mouse lncRNAs from various dimensions including sequences, promoter, global synteny, and local synteny, which led to the identification of 1,731 conserved lncRNAs with 427 high-confidence ones meeting multiple criteria. Conserved lncRNAs, compared with non-conserved ones, generally have longer gene bodies, more exons and transcripts, stronger connections with human diseases, and are more abundant and widespread across different tissues. Transcription factor (TF) profile analysis revealed a significant enrichment of TF types and numbers in the promoters of conserved lncRNAs. We further identified a set of TFs that preferentially bind to conserved lncRNAs and exert stronger regulation on conserved than non-conserved lncRNAs. Our study has reconciled some discrepant interpretations of lncRNA conservation and revealed a new set of transcriptional factors ruling the expression of conserved lncRNAs.

PPARγ (Peroxisome Proliferator-Activated Receptor γ) Deacetylation Suppresses Aging-Associated Atherosclerosis and Hypercholesterolemia.

BACKGROUND: Atherosclerosis is a medical urgency manifesting at the onset of hypercholesterolemia and is associated with aging. Activation of PPARγ (peroxisome proliferator-activated receptor γ) counteracts metabolic dysfunction influenced by aging, and its deacetylation displays an atheroprotective property. Despite the marked increase of PPARγ acetylation during aging, it is unknown whether PPARγ acetylation is a pathogenic contributor to aging-associated atherosclerosis. METHODS: Mice with constitutive deacetylation-mimetic PPARγ mutations on lysine residues K268 and K293 (2KR) in an LDL (low-density lipoprotein)-receptor knockout (Ldlr-/-) background (2KR:Ldlr-/-) were aged for 18 months on a standard laboratory diet to examine the cardiometabolic phenotype, which was confirmed in Western-type diet-fed 2KR:Ldlr+/- mice. Whole-liver RNA-sequencing and in vitro studies in bone marrow-derived macrophages were conducted to decipher the mechanism. RESULTS: In contrast to severe atherosclerosis in WT:Ldlr-/- mice, aged 2KR:Ldlr-/- mice developed little to no plaque, which was underlain by a significantly improved plasma lipid profile, with particular reductions in circulating LDL. The protection from hypercholesterolemia was recapitulated in Western-type diet-fed 2KR:Ldlr+/- mice. Liver RNA-sequencing analysis revealed suppression of liver inflammation rather than changes in cholesterol metabolism. This anti-inflammatory effect of 2KR was attributed to polarized M2 activation of macrophages. Additionally, the upregulation of core circadian component Bmal1 (brain and muscle ARNT-like 1), perceived to be involved in anti-inflammatory immunity, was observed in the liver and bone marrow-derived macrophages. CONCLUSIONS: PPARγ deacetylation in mice prevents the development of aging-associated atherosclerosis and hypercholesterolemia, in association with the anti-inflammatory phenotype of 2KR macrophages.

TOX4, an insulin receptor-independent regulator of hepatic glucose production, is activated in diabetic liver.

Increased hepatic glucose production (HGP) contributes to hyperglycemia in type 2 diabetes. Hormonal regulation of this process is primarily, but not exclusively, mediated by the AKT-FoxO1 pathway. Here, we show that cAMP and dexamethasone regulate the high-mobility group superfamily member TOX4 to mediate HGP, independent of the insulin receptor/FoxO1 pathway. TOX4 inhibition decreases glucose production in primary hepatocytes and liver and increases glucose tolerance. Combined genetic ablation of TOX4 and FoxO1 in liver has additive effects on glucose tolerance and gluconeogenesis. Moreover, TOX4 ablation fails to reverse the metabolic derangement brought by insulin receptor knockout. TOX4 expression is increased in livers of patients with steatosis and diabetes and in diet-induced obese and db/db mice. In the latter two murine models, knockdown Tox4 decreases glycemia and improves glucose tolerance. We conclude that TOX4 is an insulin receptor-independent regulator of HGP and a candidate contributor to the pathophysiology of diabetes.

SC-MEB: spatial clustering with hidden Markov random field using empirical Bayes.

Spatial transcriptomics has been emerging as a powerful technique for resolving gene expression profiles while retaining tissue spatial information. These spatially resolved transcriptomics make it feasible to examine the complex multicellular systems of different microenvironments. To answer scientific questions with spatial transcriptomics and expand our understanding of how cell types and states are regulated by microenvironment, the first step is to identify cell clusters by integrating the available spatial information. Here, we introduce SC-MEB, an empirical Bayes approach for spatial clustering analysis using a hidden Markov random field. We have also derived an efficient expectation-maximization algorithm based on an iterative conditional mode for SC-MEB. In contrast to BayesSpace, a recently developed method, SC-MEB is not only computationally efficient and scalable to large sample sizes but is also capable of choosing the smoothness parameter and the number of clusters. We performed comprehensive simulation studies to demonstrate the superiority of SC-MEB over some existing methods. We applied SC-MEB to analyze the spatial transcriptome of human dorsolateral prefrontal cortex tissues and mouse hypothalamic preoptic region. Our analysis results showed that SC-MEB can achieve a similar or better clustering performance to BayesSpace, which uses the true number of clusters and a fixed smoothness parameter. Moreover, SC-MEB is scalable to large 'sample sizes'. We then employed SC-MEB to analyze a colon dataset from a patient with colorectal cancer (CRC) and COVID-19, and further performed differential expression analysis to identify signature genes related to the clustering results. The heatmap of identified signature genes showed that the clusters identified using SC-MEB were more separable than those obtained with BayesSpace. Using pathway analysis, we identified three immune-related clusters, and in a further comparison, found the mean expression of COVID-19 signature genes was greater in immune than non-immune regions of colon tissue. SC-MEB provides a valuable computational tool for investigating the structural organizations of tissues from spatial transcriptomic data.

Adipsin promotes bone marrow adiposity by priming mesenchymal stem cells.

Background: Marrow adipose tissue (MAT) has been shown to be vital for regulating metabolism and maintaining skeletal homeostasis in the bone marrow (BM) niche. As a reflection of BM remodeling, MAT is highly responsive to nutrient fluctuations, hormonal changes, and metabolic disturbances such as obesity and diabetes mellitus. Expansion of MAT has also been strongly associated with bone loss in mice and humans. However, the regulation of BM plasticity remains poorly understood, as does the mechanism that links changes in marrow adiposity with bone remodeling. Methods: We studied deletion of Adipsin, and its downstream effector, C3, in C57BL/6 mice as well as the bone-protected PPARγ constitutive deacetylation 2KR mice to assess BM plasticity. The mice were challenged with thiazolidinedione treatment, calorie restriction, or aging to induce bone loss and MAT expansion. Analysis of bone mineral density and marrow adiposity was performed using a µCT scanner and by RNA analysis to assess adipocyte and osteoblast markers. For in vitro studies, primary bone marrow stromal cells were isolated and subjected to osteoblastogenic or adipogenic differentiation or chemical treatment followed by morphological and molecular analyses. Clinical data was obtained from samples of a previous clinical trial of fasting and high-calorie diet in healthy human volunteers. Results: We show that Adipsin is the most upregulated adipokine during MAT expansion in mice and humans in a PPARγ acetylation-dependent manner. Genetic ablation of Adipsin in mice specifically inhibited MAT expansion but not peripheral adipose depots, and improved bone mass during calorie restriction, thiazolidinedione treatment, and aging. These effects were mediated through its downstream effector, complement component C3, to prime common progenitor cells toward adipogenesis rather than osteoblastogenesis through inhibiting Wnt/ß-catenin signaling. Conclusions: Adipsin promotes new adipocyte formation and affects skeletal remodeling in the BM niche. Our study reveals a novel mechanism whereby the BM sustains its own plasticity through paracrine and endocrine actions of a unique adipokine. Funding: This work was supported by the National Institutes of Health T32DK007328 (NA), F31DK124926 (NA), R01DK121140 (JCL), R01AR068970 (BZ), R01AR071463 (BZ), R01DK112943 (LQ), R24DK092759 (CJR), and P01HL087123 (LQ).

Protocol for comprehensive RNA sequencing analysis of murine long non-coding RNAs during aging.

Comprehensive analyses of lncRNAs in aging have been lacking because previous studies have mainly focused on the protein-coding genes during aging. Here, we describe a protocol for the organism-wide analysis of murine lncRNAs during aging. We provide step-by-step instructions to identify lncRNAs that contribute to aging and to determine their underlying functions in each tissue. We further describe methods to compare the lncRNA expression patterns and dynamic changes among multiple tissues. For complete details on the use and execution of this protocol, please refer to Zhou et al. (2020).

Metabolic Health Status Contributes to Transcriptome Alternation in Human Visceral Adipose Tissue During Obesity.

OBJECTIVE: BMI is a well-established factor affecting the transcriptome profile of adipose tissue, but there are few reports on the relationship between the metabolic health status of people with obesity and the transcriptional changes, particularly in visceral adipose tissue. METHODS: Visceral adipose tissue was collected from three subgroups of patients, lean (n = 11), metabolically healthy obesity (MHO; n = 22), and metabolically unhealthy obesity (MUO; n = 26), and RNA sequencing was conducted to profile the transcriptome changes between these groups in a pairwise manner. RESULTS: Comparing MUO with lean and comparing MHO with lean revealed similar patterns in gene expression and pathway changes: obesity, regardless of metabolic health, was associated with upregulated inflammatory pathways. However, the inflammatory signature in MUO was stronger than in MHO. Pairwise comparisons among MUO, MHO, and lean samples identified 34 common differentially expressed genes; 12 out of 34 genes were associated with inflammatory pathways and exhibited a gradually increased expression pattern in the order of lean, MHO, and MUO. CONCLUSIONS: This study reveals not only that BMI plays an important role in determining the gene expression profile in visceral adipose tissue but also that a metabolically healthy condition is associated with a less inflammatory transcriptional change during obesity.

A Landscape of Murine Long Non-Coding RNAs Reveals the Leading Transcriptome Alterations in Adipose Tissue during Aging.

Aging is an inevitable process that involves profound physiological changes. Long non-coding RNAs (lncRNAs) are emerging as important regulators in various biological processes but are not systemically studied in aging. To provide an organism-wide lncRNA landscape during aging, we conduct comprehensive RNA sequencing (RNA-seq) analyses across the mouse lifespan. Of the 1,675 aging-regulated lncRNAs (AR-lncRNAs) identified, the majority are connected to inflammation-related biological pathways. AR-lncRNAs exhibit high tissue specificity; conversely, those with higher tissue specificity are preferentially regulated during aging. White adipose tissue (WAT) displays the highest number of AR-lncRNAs and develops the most dynamic crosstalk between AR-lncRNA and AR-mRNA during aging. An adipose-enriched AR-lncRNA, lnc-adipoAR1, is negatively correlated with aging, and knocking it down inhibits adipogenesis, phenocopying the compromised adipogenic capacity of aged fat. Our works together reveal AR-lncRNAs as essential components in aging and suggest that although each tissue ages in a distinct manner, WAT is a leading contributor to aging-related health decline.

A Comparison of Co-expression Networks in Silk Gland Reveals the Causes of Silk Yield Increase During Silkworm Domestication.

Long-term domestication and selective breeding have increased the silk yield of the domestic silkworm (Bombyx mori) by several times the amount of the silk yield of its wild ancestor (Bombyx mandarina). However, little is known about the molecular mechanisms behind the increase in silk yield during domestication. Based on dynamic patterns of functional divergence in the silk gland between domestic and wild silkworms, we found that at early and intermediate stages of silk gland development, the up-regulated genes of the domestic silkworm were mainly involved in DNA integration, nucleic acid binding, and transporter activity, which are related to the division and growth of cells. This has led to the posterior silk gland (PSG) of the domestic silkworm having significantly more cells ("factories" of fibroin protein synthesis) than that of the wild silkworm. At the late stage of silk gland development, the up-regulated genes in the domestic silkworm was enriched in protein processing and ribosome pathways, suggesting protein synthesis efficiency is greatly improved during silkworm domestication. While there was an increase in fibroin protein synthesis, the production of sericin protein was simultaneously reduced in the silk gland of the domestic silkworm. This reflects that domestic and wild silkworms have been under different selection pressures. Importantly, we found that the network co-expressed with the silk-coding genes of the domestic silkworm was larger than that of the wild silkworm. Furthermore, many more genes co-expressed with silk-coding genes in the domestic silkworm were subjected to artificial selection than those in the wild silkworm. Our results revealed that the increase of silk yield during silkworm domestication is involved in improvement of a biological system which includes not only expansion of "factories" (cells of PSG) of protein synthesis, but also a high expression of silk-coding genes and silk production-related genes such as biological energy, transport, and ribosome pathway genes.

Genetic and genomic analysis for cocoon yield traits in silkworm.

Domestic species provides a powerful model for examining genetic mechanisms in the evolution of yield traits. The domestic silkworm (Bombyx mori) is an important livestock species in sericulture. While the mechanisms controlling cocoon yield are largely unknown. Here, using B. mori and its wild relative B. mandarina as intercross parents, 100 BC1 individuals were sequenced by restriction site-associated DNA sequencing (RAD-Seq). The linkage map contained 9,632 markers was constructed. We performed high-resolution quantitative trait locus (QTL) mapping for four cocoon yield traits. A total of 11 QTLs were identified, including one yield-enhancing QTL from wild silkworm. By integrating population genomics and transcriptomic analysis with QTLs, some favourable genes were revealed, including 14 domestication-related genes and 71 differentially expressed genes (DEGs) in the fifth-instar larval silk gland transcriptome between B. mori and B. mandarina. The relationships between the expression of two important candidate genes (KWMTBOMO04917 and KWMTBOMO12906) and cocoon yield were supported by quantitative real-time PCR (qPCR). Our results provide some new insights into the molecular mechanisms of complex yield traits in silkworm. The combined method might be an efficient approach for identifying putative causal genes in domestic livestock and wild relatives.

Evidence of peripheral olfactory impairment in the domestic silkworms: insight from the comparative transcriptome and population genetics.

BACKGROUND: The insect olfactory system is a highly specific and sensitive chemical detector, which plays important roles in feeding, mating and finding an appropriate oviposition site. The ecological niche of Bombyx mori has changed greatly since domestication from B. mandarina, and its olfactory response to environmental odorants clearly decreased. However, the mechanisms that result in the olfactory impairment are largely unknown. RESULTS: The antennal transcriptomes were compared between the domestic and wild silkworms. Comparison of the same sex between the domestic and wild silkworms revealed 1410 and 1173 differentially expressed genes (DEGs) in males and females, respectively. To understand the olfactory impairment, we mainly focused on the olfactory-related genes. In total, 30 olfactory genes and 19 odorant-degrading enzymes (ODEs) showed differential expression in the two comparisons, in which 19 and 14 were down-regulated in the domestic silkworm, respectively. Based on population genomic data, the down-regulated odorant receptors (ORs) showed a higher ratio of unique non-synonymous polymorphisms to synonymous polymorphisms (N/S ratio) in the domestic populations than that in the wild silkworms. Furthermore, one deleterious mutation was found in OR30 of the domestic population, which was located in transmembrane helix 6 (TM6). CONCLUSIONS: Our results suggested that down-regulation of the olfactory-related genes and relaxed selection might be the major reasons for olfactory impairment of the domestic silkworm reared completely indoor environment. Reversely, wild silkworm may increase expression and remove deleterious polymorphisms of olfactory-related genes to retain sensitive olfaction.

Unexpected invasion of miniature inverted-repeat transposable elements in viral genomes.

BACKGROUND: Transposable elements (TEs) are common and often present with high copy numbers in cellular genomes. Unlike in cellular organisms, TEs were previously thought to be either rare or absent in viruses. Almost all reported TEs display only one or two copies per viral genome. In addition, the discovery of pandoraviruses with genomes up to 2.5-Mb emphasizes the need for biologists to rethink the fundamental nature of the relationship between viruses and cellular life. RESULTS: Herein, we performed the first comprehensive analysis of miniature inverted-repeat transposable elements (MITEs) in the 5170 viral genomes for which sequences are currently available. Four hundred and fifty one copies of ten miniature inverted-repeat transposable elements (MITEs) were found and each MITE had reached relatively large copy numbers (some up to 90) in viruses. Eight MITEs belonging to two DNA superfamilies (hobo/Activator/Tam3 and Chapaev-Mirage-CACTA) were for the first time identified in viruses, further expanding the organismal range of these two superfamilies. TEs may play important roles in shaping the evolution of pandoravirus genomes, which were here found to be very rich in MITEs. We also show that putative autonomous partners of seven MITEs are present in the genomes of viral hosts, suggesting that viruses may borrow the transpositional machinery of their cellular hosts' autonomous elements to spread MITEs and colonize their own genomes. The presence of seven similar MITEs in viral hosts, suggesting horizontal transfers (HTs) as the major mechanism for MITEs propagation. CONCLUSIONS: Our discovery highlights that TEs contribute to shape genome evolution of pandoraviruses. We concluded that as for cellular organisms, TEs are part of the pandoraviruses' diverse mobilome.

Identification and comparison of long non-coding RNAs in the silk gland between domestic and wild silkworms.

Under long-term artificial selection, the domestic silkworm (Bombyx mori) has increased its silk yield tremendously in comparison with its wild progenitor, Bombyx mandarina. However, the molecular mechanism of silk yield increase is still unknown. Comparative analysis of long non-coding RNAs (lncRNAs) may provide some insights into understanding this phenotypic variation. In this study, using RNA sequencing technology data of silk gland in domestic and wild silkworms, we identified 599 lncRNAs in the silk gland of the silkworm. Compared with protein-coding genes, the silk gland lncRNA genes tend to have fewer exon numbers, shorter transcript length and lower GC-content. Moreover, we found that three lncRNA genes are significantly and differentially expressed between domestic and wild silkworms. The potential targets of two differentially expressed lncRNAs (DELs) (dw4sg_0040 and dw4sg_0483) and the expression-correlated genes with the two DELs are mainly enriched in the related processes of silk protein translation. This implies that these DELs may affect the phenotypic variation in silk yield between the domestic and wild silkworms through the post-transcriptional regulation of silk protein.

iMITEdb: the genome-wide landscape of miniature inverted-repeat transposable elements in insects.

Miniature inverted-repeat transposable elements (MITEs) have attracted much attention due to their widespread occurrence and high copy numbers in eukaryotic genomes. However, the systematic knowledge about MITEs in insects and other animals is still lacking. In this study, we identified 6012 MITE families from 98 insect species genomes. Comparison of these MITEs with known MITEs in the NCBI non-redundant database and Repbase showed that 5701(∼95%) of 6012 MITE families are novel. The abundance of MITEs varies drastically among different insect species, and significantly correlates with genome size. In general, larger genomes contain more MITEs than small genomes. Furthermore, all identified MITEs were included in a newly constructed database (iMITEdb) (http://gene.cqu.edu.cn/iMITEdb/), which has functions such as browse, search, BLAST and download. Overall, our results not only provide insight on insect MITEs but will also improve assembly and annotation of insect genomes. More importantly, the results presented in this study will promote studies of MITEs function, evolution and application in insects. DATABASE URL: http://gene.cqu.edu.cn/iMITEdb/.

BmncRNAdb: a comprehensive database of non-coding RNAs in the silkworm, Bombyx mori.

BACKGROUND: Long non-coding RNAs (lncRNAs) may play critical roles in a wide range of developmental processes of higher organisms. Recently, lncRNAs have been widely identified across eukaryotes and many databases of lncRNAs have been developed for human, mouse, fruit fly, etc. However, there is rare information about them in the only completely domesticated insect, silkworm (Bombyx mori). DESCRIPTION: In this study, we systematically scanned lncRNAs using the available silkworm RNA-seq data and public unigenes. Finally, we identified and collected 6281 lncRNAs in the silkworm. Besides, we also collected 1986 microRNAs (miRNAs) from previous studies. Then, we organized them into a comprehensive and web-based database, BmncRNAdb. This database offers a user-friendly interface for data browse and online analysis as well as the three online tools for users to predict the target genes of lncRNA or miRNA. CONCLUSIONS: We have systematically identified and collected the silkworm lncRNAs and constructed a comprehensive database of the silkworm lncRNAs and miRNAs. This work gives a glimpse into lncRNAs of the silkworm and lays foundations for the ncRNAs study of the silkworm and other insects in the future. The BmncRNAdb is freely available at http://gene.cqu.edu.cn/BmncRNAdb/index.php .

Comparative analysis of the silk gland transcriptomes between the domestic and wild silkworms.

BACKGROUND: Bombyx mori was domesticated from the Chinese wild silkworm, Bombyx mandarina. Wild and domestic silkworms are good models in which to investigate genes related to silk protein synthesis that may be differentially expressed in silk glands, because their silk productions are very different. Here we used the mRNA deep sequencing (RNA-seq) approach to identify the differentially expressed genes (DEGs) in the transcriptomes of the median/posterior silk glands of two domestic and two wild silkworms. RESULTS: The results indicated that about 58% of the total genes were expressed (reads per kilo bases per million reads (RPKM) ≥ 1) in each silkworm. Comparisons of the domestic and wild silkworm transcriptomes revealed 32 DEGs, of which 16 were up-regulated in the domestic silkworms compared with in the wild silkworms, and the other 16 were up-regulated in the wild silkworms compared with in the domestic silkworms. Quantitative real-time polymerase chain reaction (qPCR) was performed for 15 randomly selected DEGs in domestic versus wild silkworms. The qPCR results were mostly consistent with the expression levels determined from the RNA-seq data. Based on a Gene Ontology (GO) enrichment analysis and manual annotation, five of the up-regulated DEGs in the wild silkworms were predicted to be involved in immune response, and seven of the up-regulated DEGs were related to the GO term "oxidoreductase activity", which is associated with antioxidant systems. In the domestic silkworms, the up-regulated DEGs were related mainly to tissue development, secretion of proteins and metabolism. CONCLUSIONS: The up-regulated DEGs in the two domestic silkworms may be involved mainly in the highly efficient biosynthesis and secretion of silk proteins, while the up-regulated DEGs in the two wild silkworms may play more important roles in tolerance to pathogens and environment adaptation. Our results provide a foundation for understanding the molecular mechanisms of the silk production difference between domestic and wild silkworms.