Genetic control of the dynamic transcriptional response to immune stimuli and glucocorticoids at single-cell resolution.

Synthetic glucocorticoids, such as dexamethasone, have been used as a treatment for many immune conditions, such as asthma and, more recently, severe COVID-19. Single-cell data can capture more fine-grained details on transcriptional variability and dynamics to gain a better understanding of the molecular underpinnings of inter-individual variation in drug response. Here, we used single-cell RNA-seq to study the dynamics of the transcriptional response to glucocorticoids in activated peripheral blood mononuclear cells from 96 African American children. We used novel statistical approaches to calculate a mean-independent measure of gene expression variability and a measure of transcriptional response pseudotime. Using these approaches, we showed that glucocorticoids reverse the effects of immune stimulation on both gene expression mean and variability. Our novel measure of gene expression response dynamics, based on the diagonal linear discriminant analysis, separated individual cells by response status on the basis of their transcriptional profiles and allowed us to identify different dynamic patterns of gene expression along the response pseudotime. We identified genetic variants regulating gene expression mean and variability, including treatment-specific effects, and showed widespread genetic regulation of the transcriptional dynamics of the gene expression response.

Extended application of genomic selection to screen multiomics data for prognostic signatures of prostate cancer.

Prognostic tests using expression profiles of several dozen genes help provide treatment choices for prostate cancer (PCa). However, these tests require improvement to meet the clinical need for resolving overtreatment, which continues to be a pervasive problem in PCa management. Genomic selection (GS) methodology, which utilizes whole-genome markers to predict agronomic traits, was adopted in this study for PCa prognosis. We leveraged The Cancer Genome Atlas (TCGA) database to evaluate the prediction performance of six GS methods and seven omics data combinations, which showed that the Best Linear Unbiased Prediction (BLUP) model outperformed the other methods regarding predictability and computational efficiency. Leveraging the BLUP-HAT method, an accelerated version of BLUP, we demonstrated that using expression data of a large number of disease-relevant genes and with an integration of other omics data (i.e. miRNAs) significantly increased outcome predictability when compared with panels consisting of a small number of genes. Finally, we developed a novel stepwise forward selection BLUP-HAT method to facilitate searching multiomics data for predictor variables with prognostic potential. The new method was applied to the TCGA data to derive mRNA and miRNA expression signatures for predicting relapse-free survival of PCa, which were validated in six independent cohorts. This is a transdisciplinary adoption of the highly efficient BLUP-HAT method and its derived algorithms to analyze multiomics data for PCa prognosis. The results demonstrated the efficacy and robustness of the new methodology in developing prognostic models in PCa, suggesting a potential utility in managing other types of cancer.

Inference of Chromosome-Length Haplotypes Using Genomic Data of Three or a Few More Single Gametes.

Compared with genomic data of individual markers, haplotype data provide higher resolution for DNA variants, advancing our knowledge in genetics and evolution. Although many computational and experimental phasing methods have been developed for analyzing diploid genomes, it remains challenging to reconstruct chromosome-scale haplotypes at low cost, which constrains the utility of this valuable genetic resource. Gamete cells, the natural packaging of haploid complements, are ideal materials for phasing entire chromosomes because the majority of the haplotypic allele combinations has been preserved. Therefore, compared with the current diploid-based phasing methods, using haploid genomic data of single gametes may substantially reduce the complexity in inferring the donor's chromosomal haplotypes. In this study, we developed the first easy-to-use R package, Hapi, for inferring chromosome-length haplotypes of individual diploid genomes with only a few gametes. Hapi outperformed other phasing methods when analyzing both simulated and real single gamete cell sequencing data sets. The results also suggested that chromosome-scale haplotypes may be inferred by using as few as three gametes, which has pushed the boundary to its possible limit. The single gamete cell sequencing technology allied with the cost-effective Hapi method will make large-scale haplotype-based genetic studies feasible and affordable, promoting the use of haplotype data in a wide range of research.

Efficient multivariate analysis algorithms for longitudinal genome-wide association studies.

MOTIVATION: Current dynamic phenotyping system introduces time as an extra dimension to genome-wide association studies (GWAS), which helps to explore the mechanism of dynamical genetic control for complex longitudinal traits. However, existing methods for longitudinal GWAS either ignore the covariance among observations of different time points or encounter computational efficiency issues. RESULTS: We herein developed efficient genome-wide multivariate association algorithms for longitudinal data. In contrast to existing univariate linear mixed model analyses, the proposed method has improved statistic power for association detection and computational speed. In addition, the new method can analyze unbalanced longitudinal data with thousands of individuals and more than ten thousand records within a few hours. The corresponding time for balanced longitudinal data is just a few minutes. AVAILABILITY AND IMPLEMENTATION: A software package to implement the efficient algorithm named GMA (https://github.com/chaoning/GMA) is available freely for interested users in relevant fields. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Analysis of trait heritability in functionally partitioned rice genomes.

Knowledge of the genetic architecture of importantly agronomical traits can speed up genetic improvement in cultivated rice (Oryza sativa L.). Many recent investigations have leveraged genome-wide association studies (GWAS) to identify single nucleotide polymorphisms (SNPs), associated with agronomic traits in various rice populations. The reported trait-relevant SNPs appear to be arbitrarily distributed along the genome, including genic and nongenic regions. Whether the SNPs in different genomic regions play different roles in trait heritability and which region is more responsible for phenotypic variation remains opaque. We analyzed a natural rice population of 524 accessions with 3,616,597 SNPs to compare the genetic contributions of functionally distinct genomic regions for five agronomic traits, i.e., yield, heading date, plant height, grain length, and grain width. An analysis of heritability in the functionally partitioned rice genome showed that regulatory or intergenic regions account for the most trait heritability. A close look at the trait-associated SNPs (TASs) indicated that the majority of the TASs are located in nongenic regions, and the genetic effects of the TASs in nongenic regions are generally greater than those in genic regions. We further compared the predictabilities using the genetic variants from genic regions with those using nongenic regions. The results revealed that nongenic regions play a more important role than genic regions in trait heritability in rice, which is consistent with findings in humans and maize. This conclusion not only offers clues for basic research to disclose genetics behind these agronomic traits, but also provides a new perspective to facilitate genomic selection in rice.

GDCRNATools: an R/Bioconductor package for integrative analysis of lncRNA, miRNA and mRNA data in GDC.

Motivation: The large-scale multidimensional omics data in the Genomic Data Commons (GDC) provides opportunities to investigate the crosstalk among different RNA species and their regulatory mechanisms in cancers. Easy-to-use bioinformatics pipelines are needed to facilitate such studies. Results: We have developed a user-friendly R/Bioconductor package, named GDCRNATools, for downloading, organizing and analyzing RNA data in GDC with an emphasis on deciphering the lncRNA-mRNA related competing endogenous RNAs regulatory network in cancers. Many widely used bioinformatics tools and databases are utilized in our package. Users can easily pack preferred downstream analysis pipelines or integrate their own pipelines into the workflow. Interactive shiny web apps built in GDCRNATools greatly improve visualization of results from the analysis. Availability and implementation: GDCRNATools is an R/Bioconductor package that is freely available at Bioconductor (http://bioconductor.org/packages/devel/bioc/html/GDCRNATools.html). Detailed instructions, manual and example code are also available in Github (https://github.com/Jialab-UCR/GDCRNATools).

Identification of optimal prediction models using multi-omic data for selecting hybrid rice.

Genomic prediction benefits hybrid rice breeding by increasing selection intensity and accelerating breeding cycles. With the rapid advancement of technology, other omic data, such as metabolomic data and transcriptomic data, are readily available for predicting breeding values for agronomically important traits. In this study, the best prediction strategies were determined for yield, 1000 grain weight, number of grains per panicle, and number of tillers per plant of hybrid rice (derived from recombinant inbred lines) by comprehensively evaluating all possible combinations of omic datasets with different prediction methods. It was demonstrated that, in rice, the predictions using a combination of genomic and metabolomic data generally produce better results than single-omics predictions or predictions based on other combined omic data. Best linear unbiased prediction (BLUP) appears to be the most efficient prediction method compared to the other commonly used approaches, including least absolute shrinkage and selection operator (LASSO), stochastic search variable selection (SSVS), support vector machines with radial basis function and epsilon regression (SVM-R(EPS)), support vector machines with radial basis function and nu regression (SVM-R(NU)), support vector machines with polynomial kernel and epsilon regression (SVM-P(EPS)), support vector machines with polynomial kernel and nu regression (SVM-P(NU)) and partial least squares regression (PLS). This study has provided guidelines for selection of hybrid rice in terms of which types of omic datasets and which method should be used to achieve higher trait predictability. The answer to these questions will benefit academic research and will also greatly reduce the operative cost for the industry which specializes in breeding and selection.

Genomic selection using principal component regression.

Many statistical methods are available for genomic selection (GS) through which genetic values of quantitative traits are predicted for plants and animals using whole-genome SNP data. A large number of predictors with much fewer subjects become a major computational challenge in GS. Principal components regression (PCR) and its derivative, i.e., partial least squares regression (PLSR), provide a solution through dimensionality reduction. In this study, we show that PCR can perform better than PLSR in cross validation. PCR often requires extracting more components to achieve the maximum predictive ability than PLSR and thus may be associated with a higher computational cost. However, application of the HAT method (a strategy of describing the relationship between the fitted and observed response variables with a hat matrix) to PCR circumvents conventional cross validation in testing predictive ability, resulting in substantially improved computational efficiency over PLSR where cross validation is mandatory. Advantages of PCR over PLSR are illustrated with a simulated trait of a hypothetical population and four agronomical traits of a rice population. The benefit of using PCR in genomic selection is further demonstrated in an effort to predict 1000 metabolomic traits and 24,973 transcriptomic traits in the same rice population.

Metabolome-wide association studies for agronomic traits of rice.

Identification of trait-associated metabolites will advance the knowledge and understanding of the biosynthetic and catabolic pathways that are relevant to the complex traits of interest. In the past, the association between metabolites (treated as quantitative traits) and genetic variants (e.g., SNPs) has been extensively studied using metabolomic quantitative trait locus (mQTL) mapping. Nevertheless, the research on the association between metabolites with agronomic traits has been inadequate. In practice, the regular approaches for QTL mapping analysis may be adopted for metabolites-phenotypes association analysis due to the similarity in data structure of these two types of researches. In the study, we compared four regular QTL mapping approaches, i.e., simple linear regression (LR), linear mixed model (LMM), Bayesian analysis with spike-slab priors (Bayes B) and least absolute shrinkage and selection operator (LASSO), by testing their performances on the analysis of metabolome-phenotype associations. Simulation studies showed that LASSO had the higher power and lower false positive rate than the other three methods. We investigated the associations of 839 metobolites with five agronomic traits in a collection of 533 rice varieties. The results implied that a total of 25 metabolites were significantly associated with five agronomic traits. Literature search and bioinformatics analysis indicated that the identified 25 metabolites are significantly involved in some growth and development processes potentially related to agronomic traits. We also explored the predictability of agronomic traits based on the 839 metabolites through cross-validation, which showed that metabolomic prediction was efficient and its application in plant breeding has been justified.

Differential gene expression in uterine endometrium during implantation in pigs.

Embryonic mortality during the implantation period strongly affects litter size in pigs. To analyze the differentially expressed genes (DEGs) in the endometrium during implantation and further to identify candidate genes for litter size, tissues of endometrial attachment sites and intersites were collected from nine pregnant sows on Days 13, 18, and 24 of pregnancy. Endometrium tissue was also collected from another three nonpregnant sows. Samples were hybridized to the porcine Agilent GeneChip microarray. The analysis of gene expression patterns over the implantation period revealed 858 DEGs at endometrial attachment sites. Comparisons of the gene files of attachment sites and intersites revealed 12, 51, and 89 DEGs on Days 13, 18, and 24 of pregnancy, respectively. Annotated function was used to identify overrepresented genetic processes, and several biological processes were considered as the most enriched. Genes related to vascular development, proteolysis, RNA metabolism and translation, protein modification, immune response, and hormone-related are discussed in detail. Then we combined microarray technology and linkage analysis to select powerful candidate genes for quantitative trait loci affecting pig litter size. Eighty-seven DEGs were located in quantitative trait loci related to litter size, that is, total number born and number born alive. Those candidate genes were thought to affect litter size by influencing embryonic implantation. Furthermore, single nucleotide polymorphism of VEGFA was shown to be associated with litter size in pigs. This study identified candidate genes for litter size that were related to embryonic implantation and could be a resource for target studies of genetic markers for litter size in pigs.

Characterization of caffeine response regulatory variants in vascular endothelial cells.

Genetic variants in gene regulatory sequences can modify gene expression and mediate the molecular response to environmental stimuli. In addition, genotype-environment interactions (GxE) contribute to complex traits such as cardiovascular disease. Caffeine is the most widely consumed stimulant and is known to produce a vascular response. To investigate GxE for caffeine, we treated vascular endothelial cells with caffeine and used a massively parallel reporter assay to measure allelic effects on gene regulation for over 43,000 genetic variants. We identified 665 variants with allelic effects on gene regulation and 6 variants that regulate the gene expression response to caffeine (GxE, false discovery rate [FDR] < 5%). When overlapping our GxE results with expression quantitative trait loci colocalized with coronary artery disease and hypertension, we dissected their regulatory mechanisms and showed a modulatory role for caffeine. Our results demonstrate that massively parallel reporter assay is a powerful approach to identify and molecularly characterize GxE in the specific context of caffeine consumption.

Integrated of multi-omics and molecular docking reveal PHGDH, PSAT1 and PSPH in the serine synthetic pathway as potential targets of T-2 toxin exposure in pig intestinal tract.

T-2 toxin (T-2) with a molecular weight of 466.52 g/mol is an inevitable mycotoxin in food products and feeds, posing a significant threat to human and animal health. However, the underlying molecular mechanisms of the cytotoxic effects of T-2 exposure on porcine intestinal epithelial cells (IPEC-J2) remain unclear. Here, we investigated the cytotoxic effects of T-2 exposure on IPEC-J2 through the detection of cell viability, cell morphology, mitochondrial membrane potential, ROS, apoptosis and autophagy. Further transcriptomic and proteomic analyses of IPEC-J2 upon T-2 exposure were performed by using RNA-seq and TMT techniques. A total of 546 differential expressed genes (DEGs) and 269 differentially expressed proteins (DEPs) were detected. Among these, 24 common DEGs/DEPs were involved in IPEC-J2 upon T-2 exposure. Interestingly, molecular docking analysis revealed potential interactions between T-2 and three key enzymes (PHGDP, PSAT1, and PSPH) in the serine biosynthesis pathway. Besides, further experimental showed that PSAT1 knockdown exacerbated T-2-induced oxidative damage. Together, our findings indicated that the serine biosynthesis pathway including PHGDP, PSAT1, PSPH genes probably acts critical roles in the regulation of T-2-induced cell damage. This study provided new insights into the global molecular effects of T-2 exposure and identified the serine biosynthesis pathway as molecular targets and potential treatment strategies against T-2.

Analysis of transcriptional changes in the immune system associated with pubertal development in a longitudinal cohort of children with asthma.

Puberty is an important developmental period marked by hormonal, metabolic and immune changes. Puberty also marks a shift in sex differences in susceptibility to asthma. Yet, little is known about the gene expression changes in immune cells that occur during pubertal development. Here we assess pubertal development and leukocyte gene expression in a longitudinal cohort of 251 children with asthma. We identify substantial gene expression changes associated with age and pubertal development. Gene expression changes between pre- and post-menarcheal females suggest a shift from predominantly innate to adaptive immunity. We show that genetic effects on gene expression change dynamically during pubertal development. Gene expression changes during puberty are correlated with gene expression changes associated with asthma and may explain sex differences in prevalence. Our results show that molecular data used to study the genetics of early onset diseases should consider pubertal development as an important factor that modifies the transcriptome.

Single Nucleus Transcriptomics Reveals Pervasive Glial Activation in Opioid Overdose Cases.

Dynamic interactions of neurons and glia in the ventral midbrain (VM) mediate reward and addiction behavior. We studied gene expression in 212,713 VM single nuclei from 95 human opioid overdose cases and drug-free controls. Chronic exposure to opioids left numerical proportions of VM glial and neuronal subtypes unaltered, while broadly affecting glial transcriptomes, involving 9.5 - 6.2% of expressed genes within microglia, oligodendrocytes, and astrocytes, with prominent activation of the immune response including interferon, NFkB signaling, and cell motility pathways, sharply contrasting with down-regulated expression of synaptic signaling and plasticity genes in VM non-dopaminergic neurons. VM transcriptomic reprogramming in the context of opioid exposure and overdose included 325 genes with genetic variation linked to substance use traits in the broader population, thereby pointing to heritable risk architectures in the genomic organization of the brain's reward circuitry.

Single nucleus transcriptomics of ventral midbrain identifies glial activation associated with chronic opioid use disorder.

Dynamic interactions of neurons and glia in the ventral midbrain mediate reward and addiction behavior. We studied gene expression in 212,713 ventral midbrain single nuclei from 95 individuals with history of opioid misuse, and individuals without drug exposure. Chronic exposure to opioids was not associated with change in proportions of glial and neuronal subtypes, however glial transcriptomes were broadly altered, involving 9.5 - 6.2% of expressed genes within microglia, oligodendrocytes, and astrocytes. Genes associated with activation of the immune response including interferon, NFkB signaling, and cell motility pathways were upregulated, contrasting with down-regulated expression of synaptic signaling and plasticity genes in ventral midbrain non-dopaminergic neurons. Ventral midbrain transcriptomic reprogramming in the context of chronic opioid exposure included 325 genes that previous genome-wide studies had linked to risk of substance use traits in the broader population, thereby pointing to heritable risk architectures in the genomic organization of the brain's reward circuitry.

Transcriptome-wide association analyses reveal the impact of regulatory variants on rice panicle architecture and causal gene regulatory networks.

Panicle architecture is a key determinant of rice grain yield and is mainly determined at the 1-2 mm young panicle stage. Here, we investigated the transcriptome of the 1-2 mm young panicles from 275 rice varieties and identified thousands of genes whose expression levels were associated with panicle traits. Multimodel association studies suggested that many small-effect genetic loci determine spikelet per panicle (SPP) by regulating the expression of genes associated with panicle traits. We found that alleles at cis-expression quantitative trait loci of SPP-associated genes underwent positive selection, with a strong preference for alleles increasing SPP. We further developed a method that integrates the associations of cis- and trans-expression components of genes with traits to identify causal genes at even small-effect loci and construct regulatory networks. We identified 36 putative causal genes of SPP, including SDT (MIR156j) and OsMADS17, and inferred that OsMADS17 regulates SDT expression, which was experimentally validated. Our study reveals the impact of regulatory variants on rice panicle architecture and provides new insights into the gene regulatory networks of panicle traits.

Photocatalytic Microbial Fuel Cells and Performance Applications: A Review.

In recent years, photocatalytic microbial fuel cells have gradually become a hot research topic in pollutant treatment, using either in situ or indirectly the oxidation of organic pollutants by catalytic materials under light and the biodegradation and mineralization of various components in wastewater by microorganisms, or through the generation of electricity by the microbial fuel cell (MFC) system to promote the photogeneration and separation of electrons and holes by the catalytic materials of the photocatalytic cell (PC) system. This study aims to provide new ideas for the development of environmentally friendly wastewater treatment technologies by investigating the use of photocatalytic cells for the efficient degradation and resource utilization of target pollutants. This study aims to raise awareness of the use of photocatalytic microbial fuel cells for pollutant degradation by providing an overview of the practical status of photocatalytic microbial fuel cells. This is achieved by reviewing the key cathode development, production capacity, and progress in the degradation of pollutants in photocatalytic microbial fuel cells. The issues facing future developments are also discussed in terms of how photocatalytic microbial fuel cells work and how they degrade pollutants. This study shows that photocatalytic microbial fuel cells are beneficial for achieving renewable energy (bioenergy, photovoltaic, etc.) capacity and dealing with environmental pollution and that this is a novel technology that deserves to be promoted to achieve the current dual carbon targets.

Genetic control of the dynamic transcriptional response to immune stimuli and glucocorticoids at single cell resolution.

Synthetic glucocorticoids, such as dexamethasone, have been used as treatment for many immune conditions, such as asthma and more recently severe COVID-19. Single cell data can capture more fine-grained details on transcriptional variability and dynamics to gain a better understanding of the molecular underpinnings of inter-individual variation in drug response. Here, we used single cell RNA-seq to study the dynamics of the transcriptional response to glucocorticoids in activated Peripheral Blood Mononuclear Cells from 96 African American children. We employed novel statistical approaches to calculate a mean-independent measure of gene expression variability and a measure of transcriptional response pseudotime. Using these approaches, we demonstrated that glucocorticoids reverse the effects of immune stimulation on both gene expression mean and variability. Our novel measure of gene expression response dynamics, based on the diagonal linear discriminant analysis, separated individual cells by response status on the basis of their transcriptional profiles and allowed us to identify different dynamic patterns of gene expression along the response pseudotime. We identified genetic variants regulating gene expression mean and variability, including treatment-specific effects, and demonstrated widespread genetic regulation of the transcriptional dynamics of the gene expression response.

Leveraging Single-Cell RNA-seq Data to Uncover the Association Between Cell Type and Chronic Liver Diseases.

BACKGROUND: Components of liver microenvironment is complex, which makes it difficult to clarify pathogenesis of chronic liver diseases (CLD). Genome-wide association studies (GWASs) have greatly revealed the role of host genetic background in CLD pathogenesis and prognosis, while single-cell RNA sequencing (scRNA-seq) enables interrogation of the cellular diversity and function of liver tissue at unprecedented resolution. Here, we made integrative analysis on the GWAS and scRNA-seq data of CLD to uncover CLD-related cell types and provide clues for understanding on the pathogenesis. METHODS: We downloaded three GWAS summary data and three scRNA-seq data on CLD. After defining the cell types for each scRNA-seq data, we used RolyPoly and LDSC-cts to integrate the GWAS and scRNA-seq. In addition, we analyzed one scRNA-seq data without association to CLD to validate the specificity of our findings. RESULTS: After processing the scRNA-seq data, we obtain about 19,002-32,200 cells and identified 10-17 cell types. For the HCC analysis, we identified the association between B cell and HCC in two datasets. RolyPoly also identified the association, when we integrated the two scRNA-seq datasets. In addition, we also identified natural killer (NK) cell as HCC-associated cell type in one dataset. In specificity analysis, we identified no significant cell type associated with HCC. As for the cirrhosis analysis, we obtained no significant related cell type. CONCLUSION: In this integrative analysis, we identified B cell and NK cell as HCC-related cell type. More attention and verification should be paid to them in future research.

A Fast and Powerful Empirical Bayes Method for Genome-Wide Association Studies.

Linear mixed model (LMM) is an efficient method for GWAS. There are numerous forms of LMM-based GWAS methods. However, improving statistical power and computing efficiency have always been the research hotspots of the LMM-based GWAS methods. Here, we proposed a fast empirical Bayes method, which is based on linear mixed models. We call it Fast-EB-LMM in short. The novelty of this method is that it uses a modified kinship matrix accounting for individual relatedness to avoid competition between the locus of interest and its counterpart in the polygene. This property has increased statistical power. We adopted two special algorithms to ease the computational burden: Eigenvalue decomposition and Woodbury matrix identity. Simulation studies showed that Fast-EB-LMM has significantly increased statistical power of marker detection and improved computational efficiency compared with two widely used GWAS methods, EMMA and EB. Real data analyses for two carcass traits in a Chinese Simmental beef cattle population showed that the significant single-nucleotide polymorphisms (SNPs) and candidate genes identified by Fast-EB-LMM are highly consistent with results of previous studies. We therefore believe that the Fast-EB-LMM method is a reliable and efficient method for GWAS.