Long-read proteogenomics to connect disease-associated sQTLs to the protein isoform effectors of disease.

A major fraction of loci identified by genome-wide association studies (GWASs) mediate alternative splicing, but mechanistic interpretation is hindered by the technical limitations of short-read RNA sequencing (RNA-seq), which cannot directly link splicing events to full-length protein isoforms. Long-read RNA-seq represents a powerful tool to characterize transcript isoforms, and recently, infer protein isoform existence. Here, we present an approach that integrates information from GWASs, splicing quantitative trait loci (sQTLs), and PacBio long-read RNA-seq in a disease-relevant model to infer the effects of sQTLs on the ultimate protein isoform products they encode. We demonstrate the utility of our approach using bone mineral density (BMD) GWAS data. We identified 1,863 sQTLs from the Genotype-Tissue Expression (GTEx) project in 732 protein-coding genes that colocalized with BMD associations (H4PP ≥ 0.75). We generated PacBio Iso-Seq data (N = ∼22 million full-length reads) on human osteoblasts, identifying 68,326 protein-coding isoforms, of which 17,375 (25%) were unannotated. By casting the sQTLs onto protein isoforms, we connected 809 sQTLs to 2,029 protein isoforms from 441 genes expressed in osteoblasts. Overall, we found that 74 sQTLs influenced isoforms likely impacted by nonsense-mediated decay and 190 that potentially resulted in the expression of unannotated protein isoforms. Finally, we functionally validated colocalizing sQTLs in TPM2, in which siRNA-mediated knockdown in osteoblasts showed two TPM2 isoforms with opposing effects on mineralization but exhibited no effect upon knockdown of the entire gene. Our approach should be to generalize across diverse clinical traits and to provide insights into protein isoform activities modulated by GWAS loci.

Bridging the splicing gap in human genetics with long-read RNA sequencing: finding the protein isoform drivers of disease.

Aberrant splicing underlies many human diseases, including cancer, cardiovascular diseases and neurological disorders. Genome-wide mapping of splicing quantitative trait loci (sQTLs) has shown that genetic regulation of alternative splicing is widespread. However, identification of the corresponding isoform or protein products associated with disease-associated sQTLs is challenging with short-read RNA-seq, which cannot precisely characterize full-length transcript isoforms. Furthermore, contemporary sQTL interpretation often relies on reference transcript annotations, which are incomplete. Solutions to these issues may be found through integration of newly emerging long-read sequencing technologies. Long-read sequencing offers the capability to sequence full-length mRNA transcripts and, in some cases, to link sQTLs to transcript isoforms containing disease-relevant protein alterations. Here, we provide an overview of sQTL mapping approaches, the use of long-read sequencing to characterize sQTL effects on isoforms, the linkage of RNA isoforms to protein-level functions and comment on future directions in the field. Based on recent progress, long-read RNA sequencing promises to be part of the human disease genetics toolkit to discover and treat protein isoforms causing rare and complex diseases.

Differential Transcriptome Responses in Human THP-1 Macrophages Following Exposure to T98G and LN-18 Human Glioblastoma Secretions: A Simplified Bioinformatics Approach to Understanding Patient-Glioma-Specific Effects on Tumor-Associated Macrophages.

A common theme in glioma disease progression is robust infiltration of immune cells within the tumor microenvironment, resulting in a state of chronic inflammation. This disease state is characterized by an abundance of CD68+ microglia and CD163+ bone marrow-derived macrophages with the greater the percentage of CD163+ cells, the poorer the prognosis. These macrophages are "cold," in that their phenotype is of an alternatively activated state (M0-M2-like) supporting tumor growth rather than being engaged with classically activated, pro-inflammatory, and anti-tumor activities, referred to as "hot", or M1-like. Herein, we have developed an in vitro approach that uses two human glioma cell lines, T98G and LN-18, which exhibit a variety of differing mutations and characteristics, to demonstrate their disparate effects on differentiated THP-1 macrophages. We first developed an approach to differentiating THP-1 monocytes to macrophages with mixed transcriptomic phenotypes we regard as M0-like macrophages. We then found that supernatants from the two different glioma cell lines induced different gene expression profiles in THP-1 macrophages, suggesting that from patient to patient, gliomas may be considered as different diseases. This study suggests that in addition to standard glioma treatment modalities, transcriptome profiling of the effects of cultured glioma cells on a standard THP-1 macrophage in vitro model may lead to future druggable targets that aim to reprogram tumor-associated macrophages towards an anti-tumor phenotype.

Using "-omics" Data to Inform Genome-wide Association Studies (GWASs) in the Osteoporosis Field.

PURPOSE OF REVIEW: Osteoporosis constitutes a major societal health problem. Genome-wide association studies (GWASs) have identified over 1100 loci influencing bone mineral density (BMD); however, few of the causal genes have been identified. Here, we review approaches that use "-omics" data and genetic- and systems genetics-based analytical strategies to facilitate causal gene discovery. RECENT FINDINGS: The bone field is beginning to adopt approaches that are commonplace in other disease disciplines. The slower progress has been due in part to the lack of large-scale "omics" data on bone and bone cells. This is however changing, and approaches such as eQTL colocalization, transcriptome-wide association studies (TWASs), network, and integrative approaches are beginning to provide significant insight into the genes responsible for BMD GWAS associations. The use of "-omics" data to inform BMD GWASs has increased in recent years, leading to the identification of novel regulators of BMD in humans. The ultimate goal will be to use this information to develop more effective therapies to treat and ultimately prevent osteoporosis.

Long-read proteogenomics to connect disease-associated sQTLs to the protein isoform effectors of disease.

A major fraction of loci identified by genome-wide association studies (GWASs) lead to alterations in alternative splicing, but interpretation of how such alterations impact proteins is hindered by the technical limitations of short-read RNA-seq, which cannot directly link splicing events to full-length transcript or protein isoforms. Long-read RNA-seq represents a powerful tool to define and quantify transcript isoforms, and recently, infer protein isoform existence. Here we present a novel approach that integrates information from GWAS, splicing QTL (sQTL), and PacBio long-read RNA-seq in a disease-relevant model to infer the effects of sQTLs on the ultimate protein isoform products they encode. We demonstrate the utility of our approach using bone mineral density (BMD) GWAS data. We identified 1,863 sQTLs from the Genotype-Tissue Expression (GTEx) project in 732 protein-coding genes which colocalized with BMD associations (H 4 PP ≥ 0.75). We generated deep coverage PacBio long-read RNA-seq data (N=∼22 million full-length reads) on human osteoblasts, identifying 68,326 protein-coding isoforms, of which 17,375 (25%) were novel. By casting the colocalized sQTLs directly onto protein isoforms, we connected 809 sQTLs to 2,029 protein isoforms from 441 genes expressed in osteoblasts. Using these data, we created one of the first proteome-scale resources defining full-length isoforms impacted by colocalized sQTLs. Overall, we found that 74 sQTLs influenced isoforms likely impacted by nonsense mediated decay (NMD) and 190 that potentially resulted in the expression of new protein isoforms. Finally, we identified colocalizing sQTLs in TPM2 for splice junctions between two mutually exclusive exons, and two different transcript termination sites, making it impossible to interpret without long-read RNA-seq data. siRNA mediated knockdown in osteoblasts showed two TPM2 isoforms with opposing effects on mineralization. We expect our approach to be widely generalizable across diverse clinical traits and accelerate system-scale analyses of protein isoform activities modulated by GWAS loci.

Mitochondrial ß-oxidation of adipose-derived fatty acids by osteoblasts fuels parathyroid hormone-induced bone formation.

The energetic costs of bone formation require osteoblasts to coordinate their activities with tissues, like adipose, that can supply energy-dense macronutrients. In the case of intermittent parathyroid hormone (PTH) treatment, a strategy used to reduce fracture risk, bone formation is preceded by a change in systemic lipid homeostasis. To investigate the requirement for fatty acid oxidation by osteoblasts during PTH-induced bone formation, we subjected mice with osteoblast-specific deficiency of mitochondrial long-chain ß-oxidation as well as mice with adipocyte-specific deficiency for the PTH receptor or adipose triglyceride lipase to an anabolic treatment regimen. PTH increased the release of fatty acids from adipocytes and ß-oxidation by osteoblasts, while the genetic mouse models were resistant to the hormone's anabolic effect. Collectively, these data suggest that PTH's anabolic actions require coordinated signaling between bone and adipose, wherein a lipolytic response liberates fatty acids that are oxidized by osteoblasts to fuel bone formation.

Single-Cell Transcriptomics of Bone Marrow Stromal Cells in Diversity Outbred Mice: A Model for Population-Level scRNA-Seq Studies.

Genome-wide association studies (GWASs) have advanced our understanding of the genetics of osteoporosis; however, the challenge has been converting associations to causal genes. Studies have utilized transcriptomics data to link disease-associated variants to genes, but few population transcriptomics data sets have been generated on bone at the single-cell level. To address this challenge, we profiled the transcriptomes of bone marrow-derived stromal cells (BMSCs) cultured under osteogenic conditions from five diversity outbred (DO) mice using single-cell RNA-seq (scRNA-seq). The goal of the study was to determine if BMSCs could serve as a model to generate cell type-specific transcriptomic profiles of mesenchymal lineage cells from large populations of mice to inform genetic studies. By enriching for mesenchymal lineage cells in vitro, coupled with pooling of multiple samples and downstream genotype deconvolution, we demonstrate the scalability of this model for population-level studies. We demonstrate that dissociation of BMSCs from a heavily mineralized matrix had little effect on viability or their transcriptomic signatures. Furthermore, we show that BMSCs cultured under osteogenic conditions are diverse and consist of cells with characteristics of mesenchymal progenitors, marrow adipogenic lineage precursors (MALPs), osteoblasts, osteocyte-like cells, and immune cells. Importantly, all cells were similar from a transcriptomic perspective to cells isolated in vivo. We employed scRNA-seq analytical tools to confirm the biological identity of profiled cell types. SCENIC was used to reconstruct gene regulatory networks (GRNs), and we observed that cell types show GRNs expected of osteogenic and pre-adipogenic lineage cells. Further, CELLECT analysis showed that osteoblasts, osteocyte-like cells, and MALPs captured a significant component of bone mineral density (BMD) heritability. Together, these data suggest that BMSCs cultured under osteogenic conditions coupled with scRNA-seq can be used as a scalable and biologically informative model to generate cell type-specific transcriptomic profiles of mesenchymal lineage cells in large populations. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Identification of Known and Novel Long Noncoding RNAs Potentially Responsible for the Effects of Bone Mineral Density (BMD) Genomewide Association Study (GWAS) Loci.

Osteoporosis, characterized by low bone mineral density (BMD), is the most common complex disease affecting bone and constitutes a major societal health problem. Genome-wide association studies (GWASs) have identified over 1100 associations influencing BMD. It has been shown that perturbations to long noncoding RNAs (lncRNAs) influence BMD and the activities of bone cells; however, the extent to which lncRNAs are involved in the genetic regulation of BMD is unknown. Here, we combined the analysis of allelic imbalance (AI) in human acetabular bone fragments with a transcriptome-wide association study (TWAS) and expression quantitative trait loci (eQTL) colocalization analysis using data from the Genotype-Tissue Expression (GTEx) project to identify lncRNAs potentially responsible for GWAS associations. We identified 27 lncRNAs in bone that are located in proximity to a BMD GWAS association and harbor single-nucleotide polymorphisms (SNPs) demonstrating AI. Using GTEx data we identified an additional 31 lncRNAs whose expression was associated (false discovery rate [FDR] correction < 0.05) with BMD through TWAS and had a colocalizing eQTL (regional colocalization probability [RCP] > 0.1). The 58 lncRNAs are located in 43 BMD associations. To further support a causal role for the identified lncRNAs, we show that 23 of the 58 lncRNAs are differentially expressed as a function of osteoblast differentiation. Our approach identifies lncRNAs that are potentially responsible for BMD GWAS associations and suggest that lncRNAs play a role in the genetics of osteoporosis. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Genetic variability affects the skeletal response to immobilization in founder strains of the diversity outbred mouse population.

Mechanical unloading decreases bone volume and strength. In humans and mice, bone mineral density is highly heritable, and in mice the response to changes in loading varies with genetic background. Thus, genetic variability may affect the response of bone to unloading. As a first step to identify genes involved in bones' response to unloading, we evaluated the effects of unloading in eight inbred mouse strains: C57BL/6J, PWK/PhJ, WSB/EiJ, A/J, 129S1/SvImJ, NOD/ShiLtJ, NZO/HlLtJ, and CAST/EiJ. C57BL/6J and NOD/ShiLtJ mice had the greatest unloading-induced loss of diaphyseal cortical bone volume and strength. NZO/HlLtJ mice had the greatest metaphyseal trabecular bone loss, and C57BL/6J, WSB/EiJ, NOD/ShiLtJ, and CAST/EiJ mice had the greatest epiphyseal trabecular bone loss. Bone loss in the epiphyses displayed the highest heritability. With immobilization, mineral:matrix was reduced, and carbonate:phosphate and crystallinity were increased. A/J mice displayed the greatest unloading-induced loss of mineral:matrix. Changes in gene expression in response to unloading were greatest in NOD/ShiLtJ and CAST/EiJ mice. The most upregulated genes in response to unloading were associated with increased collagen synthesis and extracellular matrix formation. Our results demonstrate a strong differential response to unloading as a function of strain. Diversity outbred (DO) mice are a high-resolution mapping population derived from these eight inbred founder strains. These results suggest DO mice will be highly suited for examining the genetic basis of the skeletal response to unloading.

Impact of off-road vehicles on soil and vegetation in a desert rangeland in Saudi Arabia.

Off-road vehicle driving is considered as main contributor to land degradation in arid regions. This study examined the impact of off-road vehicles (ORV) on soil and vegetation in a natural recreational desert meadow of Raudhat Khuraim, Saudi Arabia. Vegetation canopy cover and plant height away from road tracks were assessed. Also, species density and canopy cover, bare ground cover and soil attributes were assessed in four microhabitats; tracks, inter-tracks, verges, and away from vehicle tracks (undisturbed natural areas). Results show that the cover of forbs and grasses was negatively associated with distance from road verges. It was observed that the height of woody species responded negatively to distance away from tracks. Cover of native species decreased under verge, inter-track and track microhabitats giving more opportunity for weeds to flourish. Bare ground was highest (60.7%) in tracks. ORV impact on soil bulk density was clear with an increase of 38% under tracks compared to soils of undisturbed natural vegetation and a similar decrease in porosity was observed. On the other hand, soil electrical conductivity was significantly higher (5.45 mS cm-1) under disturbance compared to 1.32 mS cm-1 in undisturbed natural vegetation. Organic matter and nitrogen were not affected significantly by ORV disturbance. The results emphasize that managing off-road vehicle driving is essential for conserving native vegetation.

Lentil ( Lens culinaris Medikus) Diet Affects the Gut Microbiome and Obesity Markers in Rat.

Lentil, a moderate-energy high-protein pulse crop, provides significant amounts of essential nutrients for healthy living. The objective of this study was to determine if a lentil-based diet affects food and energy intake, body weight, percent body fat, liver weight, and body plasma triacylglycerols (TGs) as well as the composition of fecal microbiota in rats. A total of 36 Sprague-Dawley rats were treated with either a standard diet, a 3.5% high amylose corn starch diet, or a 70.8% red lentil diet for 6 weeks. By week 6, rats fed the lentil diet had significantly lower mean body weight (443 ± 47 g/rat) than those fed the control (511 ± 51 g/rat) or corn (502 ± 38 g/rat) diets. Further, mean percent body fat and TG concentration were lower, and lean body mass was higher in rats fed the lentil diet than those fed the corn diet. Fecal abundance of Actinobacteria and Bacteriodetes were greater in rats fed the lentil or corn starch diets than those fed the control diet. Fecal abundance of Firmicutes, a bacterial phylum comprising multiple pathogenic species, decreased in rats fed the lentil and high-amylose corn starch diets vs the control diet. The lentil-based diet decreased body weight, percent body fat, and plasma triacylglycerols in rats and suppressed intestinal colonization by pathogens.

Molecular and phytochemical analysis of wild type and olive cultivars grown under Saudi Arabian environment.

This study aimed to assess genetic variability at molecular and phytochemical levels among the four most commonly grown olive cultivars and the wild-type olive of Saudi Arabia. Sixty-six and 80 amplicons were generated from 9 random amplified polymorphic DNA (RAPD) and inter simple sequence repeats (ISSR) primers, each, producing an average of 95.9 and 86.44% polymorphism for the two markers, respectively. The PIC values were 82.2% for the RAPD and 85.4% for the ISSR markers and the discrimination power for both the markers was 11.1%. The UPGMA cluster analysis based on the RAPD and ISSR data resulted in the aggregation of cultivars and wild accession with a good bootstrapping value according to their origin. Furthermore, a total of 199 compounds were identified in the cultivars based on peak area, retention time, and molecular formula using GC-MS analyses of methanolic and ethanolic extracts. These compounds were classified according to their chemical class; most of them were fatty acids, alcoholic compounds, carboxylic acids, aldehydes, heterocyclic compounds, ketones, alkanes, and phenols. Genetic and phytochemical distances were significantly correlated, based on the Mantel test. The Saudi wild accession also had high numbers of fatty acids and their esters, and can be used in breeding programs for generating new genotypes with interesting characters.