Application of deep learning in radiation therapy for cancer

In recent years, with the development of artificial intelligence, deep learning has been gradually applied to clinical treatment and research. It has also found its way into the applications in radiotherapy, a crucial method for cancer treatment. This study summarizes the commonly used and latest deep learning algorithms (including transformer, and diffusion models), introduces the workflow of different radiotherapy, and illustrates the application of different algorithms in different radiotherapy modules, as well as the defects and challenges of deep learning in the field of radiotherapy, so as to provide some help for the development of automatic radiotherapy for cancer.

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.

Myokines and Biomarkers of Frailty in Older Inpatients with Undernutrition: A Prospective Study

BACKGROUND

Population aging might increase the prevalence of undernutrition in older people, which increases the risk of frailty. Numerous studies have indicated that myokines are released by skeletal myocytes in response to muscular contractions and might be associated with frailty. This study aimed to evaluate whether myokines are biomarkers of frailty in older inpatients with undernutrition.

METHODS

The frailty biomarkers were extracted from the Gene Expression Omnibus and Genecards datasets. Relevant myokines and health-related variables were assessed in 55 inpatients aged ≥ 65 years from the Peking Union Medical College Hospital prospective longitudinal frailty study. Serum was prepared for enzyme-linked immunosorbent assay using the appropriate kits. Correlations between biomarkers and frailty status were calculated by Spearman's correlation analysis. Multiple linear regression was performed to investigate the association between factors and frailty scores.

RESULTS

The prevalence of frailty was 13.21%. The bioinformatics analysis indicated that leptin, adenosine 5'-monophosphate-activated protein kinase (AMPK), irisin, decorin, and myostatin were potential biomarkers of frailty. The frailty group had significantly higher concentrations of leptin, AMPK, and MSTN than the robust group (p < 0.05). AMPK was significantly positively correlated with frailty (p < 0.05). The pre-frailty and frailty groups had significantly lower concentrations of irisin than the robust group (p < 0.05), whereas the DCN concentration did not differ among the groups. Multiple linear regression suggested that the 15 factors influencing the coefficients of association, the top 50% were the ADL score, MNA-SF score, serum albumin concentration, urination function, hearing function, leptin concentration, GDS-15 score, and MSTN concentration.

CONCLUSIONS

Proinflammatory myokines, particularly leptin, myostatin, and AMPK, negatively affect muscle mass and strength in older adults. ADL and nutritional status play major roles in the development of frailty. Our results confirm that identification of frailty relies upon clinical variables, myokine concentrations, and functional parameters, which might enable the identification and monitoring of frailty.

Use of TSAR, Thermal Shift Analysis in R, to identify Folic Acid as a Molecule that Interacts with HIV-1 Capsid

Thermal shift assay (TSA) is a versatile biophysical technique for studying protein interactions. Here, we report a free, open-source software tool TSAR (Thermal Shift Analysis in R) to expedite and automate the analysis of thermal shift data derived either from individual experiments or large screens of chemical libraries. The TSAR package incorporates multiple, dynamic workflows to facilitate the analysis of TSA data and returns publication-ready graphics or processed results. Further, the package includes a graphic user interface (GUI) that enables easy use by non-programmers, aiming to simplify TSA analysis while diversifying visualization. To exemplify the utility of TSAR we screened a chemical library of vitamins to identify molecules that interact with the capsid protein (CA) of human immunodeficiency virus type 1 (HIV-1). Our data show that hexameric CA interacts with folic acid in vitro.

GWAS meta-analysis of psoriasis identifies new susceptibility alleles impacting disease mechanisms and therapeutic targets

Psoriasis is a common, debilitating immune-mediated skin disease. Genetic studies have identified biological mechanisms of psoriasis risk, including those targeted by effective therapies. However, the genetic liability to psoriasis is not fully explained by variation at robustly identified risk loci. To move towards a saturation map of psoriasis susceptibility we meta-analysed 18 GWAS comprising 36,466 cases and 458,078 controls and identified 109 distinct psoriasis susceptibility loci, including 45 that have not been previously reported. These include susceptibility variants at loci in which the therapeutic targets IL17RA and AHR are encoded, and deleterious coding variants supporting potential new drug targets (including in STAP2, CPVL and POU2F3). We conducted a transcriptome-wide association study to identify regulatory effects of psoriasis susceptibility variants and cross-referenced these against single cell expression profiles in psoriasis-affected skin, highlighting roles for the transcriptional regulation of haematopoietic cell development and epigenetic modulation of interferon signalling in psoriasis pathobiology.

Metabolome-wide Mendelian randomization characterizes heterogeneous and shared causal effects of metabolites on human health

Metabolites are small molecules that are useful for estimating disease risk and elucidating disease biology. Nevertheless, their causal effects on human diseases have not been evaluated comprehensively. We performed two-sample Mendelian randomization to systematically infer the causal effects of 1,099 plasma metabolites measured in 6,136 Finnish men from the METSIM study on risk of 2,099 binary disease endpoints measured in 309,154 Finnish individuals from FinnGen. We identified evidence for 282 causal effects of 70 metabolites on 183 disease endpoints (FDR<1%). We found 25 metabolites with potential causal effects across multiple disease domains, including ascorbic acid 2-sulfate affecting 26 disease endpoints in 12 disease domains. Our study suggests that N-acetyl-2-aminooctanoate and glycocholenate sulfate affect risk of atrial fibrillation through two distinct metabolic pathways and that N-methylpipecolate may mediate the causal effect of N6, N6-dimethyllysine on anxious personality disorder. This study highlights the broad causal impact of plasma metabolites and widespread metabolic connections across diseases.

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.

JUMP: replicability analysis of high-throughput experiments with applications to spatial transcriptomic studies

MOTIVATION

Replicability is the cornerstone of scientific research. The current statistical method for high-dimensional replicability analysis either cannot control the false discovery rate (FDR) or is too conservative.

RESULTS

We propose a statistical method, JUMP, for the high-dimensional replicability analysis of two studies. The input is a high-dimensional paired sequence of p-values from two studies and the test statistic is the maximum of p-values of the pair. JUMP uses four states of the p-value pairs to indicate whether they are null or non-null. Conditional on the hidden states, JUMP computes the cumulative distribution function of the maximum of p-values for each state to conservatively approximate the probability of rejection under the composite null of replicability. JUMP estimates unknown parameters and uses a step-up procedure to control FDR. By incorporating different states of composite null, JUMP achieves a substantial power gain over existing methods while controlling the FDR. Analyzing two pairs of spatially resolved transcriptomic datasets, JUMP makes biological discoveries that otherwise cannot be obtained by using existing methods.

AVAILABILITY AND IMPLEMENTATION

An R package JUMP implementing the JUMP method is available on CRAN (https://CRAN.R-project.org/package=JUMP).

[Oncologic outcomes of surgical treatments of advanced sinonasal malignancies]

Objective: To investigate the prognoses of advanced (T3-T4) sinonasal malignancies (SNM). Methods: The clinical data of 229 patients with advanced (T3-4) SNM who underwent surgical treatments in the First Affiliated Hospital of Sun Yat-sen University from 2000 to 2018 were retrospectively analyzed, including 162 males and 67 females, aged (46.8±18.5) years old. Among them, 167 cases received endoscopic surgery alone, 30 cases received assisted incision endoscopic surgery, and 32 cases received open surgery. The Kaplan-Meier method was used to estimate the 3-year and 5-year overall survival (OS) and event-free survival (EFS). Univariate and multivariate Cox regression analyses were performed to explore significant prognostic factors. Results: The 3-year and 5-year OS were respectively 69.7% and 64.0%. The median OS time was 43 months. The 3-year and 5-year EFS were respectively 57.8% and 47.4%. The median EFS time was 34 months. The 5-year OS of the patients with epithelial-derived tumors was better than that of the patients with mesenchymal-derived tumors and malignant melanoma (5-year OS was respectively 72.3%, 47.8% and 30.0%, χ²=36.01, P<0.001). Patients with microscopically margin-negative resection (R0 resection) had the best prognosis, followed by macroscopically margin-negative resection (R1 resection), and debulking surgery was the worst (5-year OS was respectively 78.4%, 55.1% and 37.4%, χ²=24.63, P<0.001). There was no significant difference in 5-year OS between the endoscopic surgery group and the open surgery group (65.8% vs. 53.4%, χ²=2.66, P=0.102). Older patients had worse OS (HR=1.02, P=0.011) and EFS (HR=1.01, P=0.027). Patients receiving adjuvant therapy had a lower risk of death (HR=0.62, P=0.038). Patients with a history of nasal radiotherapy had a higher risk of recurrence (HR=2.48, P=0.002) and a higher risk of death (HR=2.03, P=0.020). Conclusion: For patients with advanced SNM, the efficacy of endoscopic surgery can be comparable to that of open surgery when presence of safe surgical margins, and a treatment plan based on transnasal endoscopic surgery as the main comprehensive treatment is recommended.

Mapping genomic regulation of kidney disease and traits through high-resolution and interpretable eQTLs

Expression quantitative trait locus (eQTL) studies illuminate genomic variants that regulate specific genes and contribute to fine-mapped loci discovered via genome-wide association studies (GWAS). Efforts to maximize their accuracy are ongoing. Using 240 glomerular (GLOM) and 311 tubulointerstitial (TUBE) micro-dissected samples from human kidney biopsies, we discovered 5371 GLOM and 9787 TUBE genes with at least one variant significantly associated with expression (eGene) by incorporating kidney single-nucleus open chromatin data and transcription start site distance as an "integrative prior" for Bayesian statistical fine-mapping. The use of an integrative prior resulted in higher resolution eQTLs illustrated by (1) smaller numbers of variants in credible sets with greater confidence, (2) increased enrichment of partitioned heritability for GWAS of two kidney traits, (3) an increased number of variants colocalized with the GWAS loci, and (4) enrichment of computationally predicted functional regulatory variants. A subset of variants and genes were validated experimentally in vitro and using a Drosophila nephrocyte model. More broadly, this study demonstrates that tissue-specific eQTL maps informed by single-nucleus open chromatin data have enhanced utility for diverse downstream analyses.

Association of pectoralis muscle area on computed tomography with airflow limitation severity and respiratory outcomes in COPD: A population-based prospective cohort study

BACKGROUND

Previous studies have shown that patients with chronic obstructive pulmonary disease (COPD) of severe or very severe airflow limitation have a reduced pectoralis muscle area (PMA), which is associated with mortality. However, whether patients with COPD of mild or moderate airflow limitation also have a reduced PMA remains unclear. Additionally, limited evidence is available regarding the associations between PMA and respiratory symptoms, lung function, computed tomography (CT) imaging, lung function decline, and exacerbations. Therefore, we conducted this study to evaluate the presence of PMA reduction in COPD and to clarify its associations with the referred variables.

METHODS

This study was based on the subjects enrolled from July 2019 to December 2020 in the Early Chronic Obstructive Pulmonary Disease (ECOPD) study. Data including questionnaire, lung function, and CT imaging were collected. The PMA was quantified on full-inspiratory CT at the aortic arch level using predefined -50 and 90 Hounsfield unit attenuation ranges. Multivariate linear regression analyses were performed to assess the association between the PMA and airflow limitation severity, respiratory symptoms, lung function, emphysema, air trapping, and the annual decline in lung function. Cox proportional hazards analysis and Poisson regression analysis were used to evaluate the PMA and exacerbations after adjustment.

RESULTS

We included 1352 subjects at baseline (667 with normal spirometry, 685 with spirometry-defined COPD). The PMA was monotonically lower with progressive airflow limitation severity of COPD after adjusting for confounders (vs. normal spirometry; Global Initiative for Chronic Obstructive Lung Disease [GOLD] 1: β=-1.27, P=0.028; GOLD 2: β=-2.29, P<0.001; GOLD 3: β=-4.88, P<0.001; GOLD 4: β=-6.47, P=0.014). The PMA was negatively associated with the modified British Medical Research Council dyspnea scale (β=-0.005, P=0.026), COPD Assessment Test score (β=-0.06, P=0.001), emphysema (β=-0.07, P<0.001), and air trapping (β=-0.24, P<0.001) after adjustment. The PMA was positively associated with lung function (all P<0.05). Similar associations were discovered for the pectoralis major muscle area and pectoralis minor muscle area. After the 1-year follow-up, the PMA was associated with the annual decline in the post-bronchodilator forced expiratory volume in 1 s percent of predicted value (β=0.022, P=0.002) but not with the annual rate of exacerbations or the time to first exacerbation.

CONCLUSION

Patients with mild or moderate airflow limitation exhibit a reduced PMA. The PMA is associated with airflow limitation severity, respiratory symptoms, lung function, emphysema, and air trapping, suggesting that PMA measurement can assist with COPD assessment.

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.

Probabilistic integration of transcriptome-wide association studies and colocalization analysis identifies key molecular pathways of complex traits

Integrative genetic association methods have shown great promise in post-GWAS (genome-wide association study) analyses, in which one of the most challenging tasks is identifying putative causal genes and uncovering molecular mechanisms of complex traits. Recent studies suggest that prevailing computational approaches, including transcriptome-wide association studies (TWASs) and colocalization analysis, are individually imperfect, but their joint usage can yield robust and powerful inference results. This paper presents INTACT, a computational framework to integrate probabilistic evidence from these distinct types of analyses and implicate putative causal genes. This procedure is flexible and can work with a wide range of existing integrative analysis approaches. It has the unique ability to quantify the uncertainty of implicated genes, enabling rigorous control of false-positive discoveries. Taking advantage of this highly desirable feature, we further propose an efficient algorithm, INTACT-GSE, for gene set enrichment analysis based on the integrated probabilistic evidence. We examine the proposed computational methods and illustrate their improved performance over the existing approaches through simulation studies. We apply the proposed methods to analyze the multi-tissue eQTL data from the GTEx project and eight large-scale complex- and molecular-trait GWAS datasets from multiple consortia and the UK Biobank. Overall, we find that the proposed methods markedly improve the existing putative gene implication methods and are particularly advantageous in evaluating and identifying key gene sets and biological pathways underlying complex traits.

[Effects of insulin glargine at different times on organ oxidative stress in burned rats with delayed resuscitation]

Objective: To examine the antioxidant effect of low dosage insulin glargine intervention at different time in rats with delayed resuscitation after burn, in order to acquire a better time of antioxidant intervention during delayed resuscitation following burn injury. Methods: With 10 rats in each group, 50 male SD rats were assigned to sham injury group, delayed resuscitation group, immediate post-burn insulin glargine treatment group (immediate treatment group), 2 hours post-burn insulin glargine treatment group(2 h treatment group), and 6 hours post-burn insulin treatment group(6 h treatment group) with random number table. Each treatment group received subcutaneous injections of insulin glargine (1.0 U·kg^-1·d^-1) immediately, two hours and six hours after the burn, while the delayed resuscitation group received the same amount of normal saline six hours after the burn. To imitate delayed fluid resuscitation, the delayed resuscitation group and each therapy group were intraperitoneally injected with normal saline (40 ml/kg) 6 hours after injury. No medicine and fluid resuscitation was administered to the sham injury group. Rats in the sham injury group had their abdominal aortic blood, hearts, and kidney tissues collected immediately after injury, while rats in the other groups had their blood and tissues collected 24 hours later. To analyze the timing of antioxidant intervention, the activities of CuZn-superoxide dismutase (CuZn-SOD), catalase (CAT), glutathione peroxidase (GSH-Px), total antioxidant capacity (T-AOC), xanthine oxidase (XOD) and myeloperoxidase (MPO) in blood glucose and myocardial and renal tissues were measured by spectrophotometry. Results: Compared with the sham group, blood glucose levels in the delayed resuscitation group increased [(10.72±0.80) vs (6.57±0.82)mmol/L,P<0.001], while in the myocardium and kidney, the activities of CuZn-SOD, CAT, GSH-Px and T-AOC decreased (all P<0.05) and the activities of XOD and MPO increased (all P<0.05). Compared with the delayed resuscitation group, blood glucose decreased in the immediate, 2 h, and 6 h treatment groups (all P<0.05). In the immediate and 2 h treatment group, the activities of CuZn-SOD, CAT, GSH-Px and T-AOC in the myocardium and kidney increased(all P<0.05). In the 6 h treatment group, only the activities of GSH-Px in myocardium, CAT and GSH-Px in kidney increased (all P<0.05). Compared with the delayed resuscitation group, in the immediate treatment group, the activities of MPO and XOD in myocardial tissue and XOD in renal tissue decreased (all P<0.05). The activities of MPO and XOD in myocardial and renal tissues of the 2 h treatment group both decreased (all P<0.05). In the 6 h treatment group, the activities of MPO in myocardial tissue and XOD in renal tissue both decreased (all P<0.05). Compared with the immediate treatment group, the activity of GSH-Px in myocardial tissue increased (P<0.05), and the activities of CuZn-SOD, CAT, GSH-Px and T-AOC in renal tissue increased in the 2 h treatment group (all P<0.05). The activities of CuZn-SOD, CAT, GSH-Px and T-AOC in myocardium of 6 h treatment group decreased (all P<0.05). Compared with the immediate treatment group, the activities of XOD and MPO in myocardial tissue and XOD in renal tissue of the 2 h treatment group had no significant difference (all P>0.05), but the activity of MPO in renal tissue decreased (P<0.05). The activities of XOD and MPO in myocardial tissue of the 6 h treatment group increased (all P<0.05). Compared with the 2 h treatment group, the activities of CuZn-SOD, CAT and GSH-Px and T-AOC in myocardium and kidney tissues in the 6 h treatment group decreased (all P<0.05), while the activities of XOD and MPO in myocardium and kidney tissues increased [myocardium: (374±8) vs (290±19) U/g, (0.021 8±0.003 9) vs (0.010 7±0.002 4) U/g, kidney: (157±6) vs (128±9) U/g, (0.026 8±0.004 3) vs (0.013 4±0.003 1) U/g, all P<0.05]. Conclusions: The timing of the intervention is related to the antioxidant effect of insulin glargine during delayed burn resuscitation. The intervention immediately and 2 hours after burn could have a better antioxidant effect compared to the intervention at 6 hours after burn.

Integrating transcriptomics, metabolomics, and GWAS helps reveal molecular mechanisms for metabolite levels and disease risk

Transcriptomics data have been integrated with genome-wide association studies (GWASs) to help understand disease/trait molecular mechanisms. The utility of metabolomics, integrated with transcriptomics and disease GWASs, to understand molecular mechanisms for metabolite levels or diseases has not been thoroughly evaluated. We performed probabilistic transcriptome-wide association and locus-level colocalization analyses to integrate transcriptomics results for 49 tissues in 706 individuals from the GTEx project, metabolomics results for 1,391 plasma metabolites in 6,136 Finnish men from the METSIM study, and GWAS results for 2,861 disease traits in 260,405 Finnish individuals from the FinnGen study. We found that genetic variants that regulate metabolite levels were more likely to influence gene expression and disease risk compared to the ones that do not. Integrating transcriptomics with metabolomics results prioritized 397 genes for 521 metabolites, including 496 previously identified gene-metabolite pairs with strong functional connections and suggested 33.3% of such gene-metabolite pairs shared the same causal variants with genetic associations of gene expression. Integrating transcriptomics and metabolomics individually with FinnGen GWAS results identified 1,597 genes for 790 disease traits. Integrating transcriptomics and metabolomics jointly with FinnGen GWAS results helped pinpoint metabolic pathways from genes to diseases. We identified putative causal effects of UGT1A1/UGT1A4 expression on gallbladder disorders through regulating plasma (E,E)-bilirubin levels, of SLC22A5 expression on nasal polyps and plasma carnitine levels through distinct pathways, and of LIPC expression on age-related macular degeneration through glycerophospholipid metabolic pathways. Our study highlights the power of integrating multiple sets of molecular traits and GWAS results to deepen understanding of disease pathophysiology.

Toxicity Profiles and Protective Effects of Antifreeze Proteins From Insect in Mammalian Models

Antifreeze proteins (AFPs), found in many cold-adapted organisms, can protect them from cold and freezing damages and have thus been considered as additional protectants in current cold tissue preservation solutions that generally include electrolytes, osmotic agents, colloids and antioxidants, to reduce the loss of tissue viability associated with cold-preservation. Due to the lack of toxicity profile studies on AFPs, their inclusion in cold preservation solutions has been a trial-and-error process limiting the development of AFPs' application in cold preservation. To assess the feasibility of translating the technology of AFPs for mammalian cell cold or cryopreservation, we determined the toxicity profile of two highly active beetle AFPs, DAFP1 and TmAFP, from Dendroides canadensis and Tenebrio molitor in this study. Toxicity was examined on a panel of representative mammalian cell lines including testicular spermatogonial stem cells and Leydig cells, macrophages, and hepatocytes. Treatments with DAFP1 and TmAFP at up to 500μg/mL for 48 and 72hours were safe in three of the cell lines, except for a 20% decrease in spermatogonia treated with TmAFP. However, both AFPs at 500μg/mL or below reduced hepatocyte viability by 20 to 40% at 48 and 72h. At 1000μg/mL, DAFP1 and TmAFP reduced viability in most cell lines. While spermatogonia and Leydig cell functions were not affected by 1000μg/mL DAFP1, this treatment induced inflammatory responses in macrophages. Adding 1000μg/ml DAFP1 to rat kidneys stored at 4°C for 48hours protected the tissues from cold-related damage, based on tissue morphology and gene and protein expression of two markers of kidney function. However, DAFP1 and TmAFP did not prevent the adverse effects of cold on kidneys over 72hours. Overall, DAFP1 is less toxic at high dose than TmAFP, and has potential for use in tissue preservation at doses up to 500μg/mL. However, careful consideration must be taken due to the proinflammatory potential of DAFP1 on macrophages at higher doses and the heighten susceptibility of hepatocytes to both AFPs.

Analyzing and reconciling colocalization and transcriptome-wide association studies from the perspective of inferential reproducibility

Transcriptome-wide association studies and colocalization analysis are popular computational approaches for integrating genetic-association data from molecular and complex traits. They show the unique ability to go beyond variant-level genetic-association evidence and implicate critical functional units, e.g., genes, in disease etiology. However, in practice, when the two approaches are applied to the same molecular and complex-trait data, the inference results can be markedly different. This paper systematically investigates the inferential reproducibility between the two approaches through theoretical derivation, numerical experiments, and analyses of four complex trait GWAS and GTEx eQTL data. We identify two classes of inconsistent inference results. We find that the first class of inconsistent results (i.e., genes with strong colocalization but weak transcriptome-wide association study [TWAS] signals) might suggest an interesting biological phenomenon, i.e., horizontal pleiotropy; thus, the two approaches are truly complementary. The inconsistency in the second class (i.e., genes with weak colocalization but strong TWAS signals) can be understood and effectively reconciled. To this end, we propose a computational approach for locus-level colocalization analysis. We demonstrate that the joint TWAS and locus-level colocalization analysis improves specificity and sensitivity for implicating biologically relevant genes.

Genome-wide association studies of metabolites in Finnish men identify disease-relevant loci

Few studies have explored the impact of rare variants (minor allele frequency < 1%) on highly heritable plasma metabolites identified in metabolomic screens. The Finnish population provides an ideal opportunity for such explorations, given the multiple bottlenecks and expansions that have shaped its history, and the enrichment for many otherwise rare alleles that has resulted. Here, we report genetic associations for 1391 plasma metabolites in 6136 men from the late-settlement region of Finland. We identify 303 novel association signals, more than one third at variants rare or enriched in Finns. Many of these signals identify genes not previously implicated in metabolite genome-wide association studies and suggest mechanisms for diseases and disease-related traits.

Interspecies variation in hominid gut microbiota controls host gene regulation

The gut microbiome exhibits extreme compositional variation between hominid hosts. However, it is unclear how this variation impacts host physiology across species and whether this effect can be mediated through microbial regulation of host gene expression in interacting epithelial cells. Here, we characterize the transcriptional response of human colonic epithelial cells in vitro to live microbial communities extracted from humans, chimpanzees, gorillas, and orangutans. We find that most host genes exhibit a conserved response, whereby they respond similarly to the four hominid microbiomes. However, hundreds of host genes exhibit a divergent response, whereby they respond only to microbiomes from specific host species. Such genes are associated with intestinal diseases in humans, including inflammatory bowel disease and Crohn’s disease. Last, we find that inflammation-associated microbial species regulate the expression of host genes previously associated with inflammatory bowel disease, suggesting health-related consequences for species-specific host-microbiome interactions across hominids.

Muehlbauer et al. investigate how variation between different hominid microbiomes drives host gene expression in colonic epithelial cell cultures. They find that host genes that respond only to microbiomes from a specific hominid species are linked to gastrointestinal diseases, suggesting implications for understanding how the microbiome can impact human health.

[Mage-D1 binding to activated p75NTR positively regulates mineralization of rat ectomesenchymal stem cells in vitro]

OBJECTIVE

To detect the binding of Mage-D1 with activated p75NTR and explore their role in regulating mineralization of ectomesenchymal stem cells (EMSCs).

METHODS

EMSCs were isolated from the tooth germs of embryonic SD rats (19.5 days of gestation) by tissue explant culture and were identified for surface markers using flow cytometry. The cultured cells were divided into blank control group, 100 ng/mL nerve growth factor (NGF) stimulation group, and lentivirus-mediated Mage-D1 interference (SH-Mage-D1) group. Proximity ligation assay was used to detect the binding of Mage-D1 with activated p75NTR in the EMSCs, and the binding strength was compared among the 3 groups. Alizarin red staining and ALP staining were used to observe mineralization of the induced cells. The expressions of ALP, Runx2, OCN, BSP, OPN, Msx1 and Dlx1 at both the mRNA and protein levels were detected using RT-PCR and Western blotting.

RESULTS

The isolated EMSCs expressed high levels of cell surface markers CD44, CD90, CD29, CD146, and CD105 with a low expression of CD45. The results of proximity ligation assay showed that the binding of Mage-D1 with activated p75NTR in the cells increased over time, and the binding strength was significantly greater in NFG-treated cells than in the cells in the other two groups (P < 0.05). Alizarin red staining and ALP staining of the induced cells showed that the changes in the mineralization nodules were consistent with those of ALP activity. The cells treated with 100 ng/mL NGF exhibited significantly increased expressions of ALP, Runx2, OCN, BSP, OPN, Col1, Msx1 and Dlx1 as compared with the cells in the other two groups (P < 0.05).

CONCLUSION

Mage-D1 directly binds to activated p75NTR in embryonic rat EMSCs to positively regulate the mineralization of the EMSCs.

Psychosocial experiences modulate asthma-associated genes through gene-environment interactions

Social interactions and the overall psychosocial environment have a demonstrated impact on health, particularly for people living in disadvantaged urban areas. Here, we investigated the effect of psychosocial experiences on gene expression in peripheral blood immune cells of children with asthma in Metro Detroit. Using RNA-sequencing and a new machine learning approach, we identified transcriptional signatures of 19 variables including psychosocial factors, blood cell composition, and asthma symptoms. Importantly, we found 169 genes associated with asthma or allergic disease that are regulated by psychosocial factors and 344 significant gene-environment interactions for gene expression levels. These results demonstrate that immune gene expression mediates the link between negative psychosocial experiences and asthma risk.

Functional dynamic genetic effects on gene regulation are specific to particular cell types and environmental conditions

Genetic effects on gene expression and splicing can be modulated by cellular and environmental factors; yet interactions between genotypes, cell type, and treatment have not been comprehensively studied together. We used an induced pluripotent stem cell system to study multiple cell types derived from the same individuals and exposed them to a large panel of treatments. Cellular responses involved different genes and pathways for gene expression and splicing and were highly variable across contexts. For thousands of genes, we identified variable allelic expression across contexts and characterized different types of gene-environment interactions, many of which are associated with complex traits. Promoter functional and evolutionary features distinguished genes with elevated allelic imbalance mean and variance. On average, half of the genes with dynamic regulatory interactions were missed by large eQTL mapping studies, indicating the importance of exploring multiple treatments to reveal previously unrecognized regulatory loci that may be important for disease.

Prognostic analysis of Behçet's disease with aortic regurgitation or involvement

Background

Aortic regurgitation is the most common cardiovascular damage in Chinese patients with Behçet’s disease (BD) and is usually associated with aortic disease. These patients are easily misdiagnosed, and their prognosis is poor, even after surgical treatment. This study aimed to analyse potential factors that can improve the prognosis of BD patients with aortic regurgitation and/or aortic involvement.

Methods

Twenty-two patients with diagnosed or suspected BD as well as aortic regurgitation and/or aortic involvement in our hospital from 2012 through 2017 were collected in this study. Their clinical characteristics were listed, and the diagnosis of BD was evaluated by two different criteria sets. The influences of surgical treatment and immunosuppressive therapy (IST) on their prognosis were also explored.

Results

The diagnostic positive rate of the International Criteria for Behçet’s Disease was higher than that of the International Study Group criteria (kappa value 0.31, p < 0.05), indicating that the diagnostic consistency between the criteria sets was poor. There was no significant difference in survival between patients who had undergone ≤ 1 operation and those with ≥ 2 operations. Aortic valve replacement alone or in combination with aortic root replacement had no significant effect on the incidence of reoperation or death, but IST did significantly reduce this incidence (p < 0.05). However, there was no significant difference in the occurrence of reoperation or death between preoperative and postoperative IST versus postoperative IST only.

Conclusion

IST significantly improved the prognosis of BD patients with aortic regurgitation and/or aortic involvement.

Supplementary Information

The online version of this article (10.1007/s12471-021-01567-6) contains supplementary material, which is available to authorized users.

Population-scale tissue transcriptomics maps long non-coding RNAs to complex disease

Long non-coding RNA (lncRNA) genes have well-established and important impacts on molecular and cellular functions. However, among the thousands of lncRNA genes, it is still a major challenge to identify the subset with disease or trait relevance. To systematically characterize these lncRNA genes, we used Genotype Tissue Expression (GTEx) project v8 genetic and multi-tissue transcriptomic data to profile the expression, genetic regulation, cellular contexts, and trait associations of 14,100 lncRNA genes across 49 tissues for 101 distinct complex genetic traits. Using these approaches, we identified 1,432 lncRNA gene-trait associations, 800 of which were not explained by stronger effects of neighboring protein-coding genes. This included associations between lncRNA quantitative trait loci and inflammatory bowel disease, type 1 and type 2 diabetes, and coronary artery disease, as well as rare variant associations to body mass index.

The long non-coding RNA DDX11-AS1 facilitates cell progression and oxaliplatin resistance via regulating miR-326/IRS1 axis in gastric cancer

OBJECTIVE

The long non-coding RNA DDX11 antisense RNA 1 (DDX11-AS1) was found to be highly expressed in gastric cancer (GC). This study was to explore the role and molecular mechanism in oxaliplatin (OXA) resistance.

PATIENTS AND METHODS

The levels of DDX11-AS1, microRNA-326 (miR-326) and insulin receptor substrate 1 (IRS1) were measured by quantitative Real-time polymerase chain reaction (qRT-PCR). Cell proliferation, migration, invasion and apoptosis were examined by methylthiazolyldiphenyl-tetrazolium bromide (MTT), transwell and flow cytometry assays, respectively. Levels of all protein were detected using Western blot. The correlation between miR-326 and DDX11-AS1/IRS1 was confirmed by Dual-Luciferase reporter and RNA immunoprecipitation (RIP) assays. The xenograft model was constructed to explore the effect of DDX11-AS1 in vivo.

RESULTS

DDX11-AS1 was overexpressed in OXA-resistant GC tissues and cells, and DDX11-AS1 knockdown inhibited cell proliferation, migration, invasion and OXA resistance, and promoted apoptosis in OXA-resistant GC cells. Mechanically, DDX11-AS1 directly targeted miR-326 and miR-326 could bind to IRS1 in OXA-resistant GC cells. Functionally, silencing DDX11-AS1 repressed the progression and OXA resistance in OXA-resistant GC cells by down-modulating IRS1 expression via sponging miR-326 in vitro and in vivo.

CONCLUSIONS

DDX11-AS1 accelerated the progression and OXA chemoresistance of GC cells in vitro and in vivo by sponging miR-326, thus increasing the expression of IRS1, suggesting DDX11-AS1 might be a promising prognostic biomarker and therapeutic target in GC.

Revisiting the genome-wide significance threshold for common variant GWAS

Over the last decade, GWAS meta-analyses have used a strict P-value threshold of 5 × 10-8 to classify associations as significant. Here, we use our current understanding of frequently studied traits including lipid levels, height, and BMI to revisit this genome-wide significance threshold. We compare the performance of studies using the P = 5 × 10-8 threshold in terms of true and false positive rate to other multiple testing strategies: (1) less stringent P-value thresholds, (2) controlling the FDR with the Benjamini-Hochberg and Benjamini-Yekutieli procedure, and (3) controlling the Bayesian FDR with posterior probabilities. We applied these procedures to re-analyze results from the Global Lipids and GIANT GWAS meta-analysis consortia and supported them with extensive simulation that mimics the empirical data. We observe in simulated studies with sample sizes ∼20,000 and >120,000 that relaxing the P-value threshold to 5 × 10-7 increased discovery at the cost of 18% and 8% of additional loci being false positive results, respectively. FDR and Bayesian FDR are well controlled for both sample sizes with a few exceptions that disappear under a less stringent definition of true positives and the two approaches yield similar results. Our work quantifies the value of using a relaxed P-value threshold in large studies to increase their true positive discovery but also show the excess false positive rates due to such actions in modest-sized studies. These results may guide investigators considering different thresholds in replication studies and downstream work such as gene-set enrichment or pathway analysis. Finally, we demonstrate the viability of FDR-controlling procedures in GWAS.

Exploiting the GTEx resources to decipher the mechanisms at GWAS loci

The resources generated by the GTEx consortium offer unprecedented opportunities to advance our understanding of the biology of human diseases. Here, we present an in-depth examination of the phenotypic consequences of transcriptome regulation and a blueprint for the functional interpretation of genome-wide association study-discovered loci. Across a broad set of complex traits and diseases, we demonstrate widespread dose-dependent effects of RNA expression and splicing. We develop a data-driven framework to benchmark methods that prioritize causal genes and find no single approach outperforms the combination of multiple approaches. Using colocalization and association approaches that take into account the observed allelic heterogeneity of gene expression, we propose potential target genes for 47% (2519 out of 5385) of the GWAS loci examined.

Probabilistic colocalization of genetic variants from complex and molecular traits: promise and limitations

Colocalization analysis has emerged as a powerful tool to uncover the overlapping of causal variants responsible for both molecular and complex disease phenotypes. The findings from colocalization analysis yield insights into the molecular pathways of complex diseases. In this paper, we conduct an in-depth investigation of the promise and limitations of the available colocalization analysis approaches. Focusing on variant-level colocalization approaches, we first establish the connections between various existing methods. We proceed to discuss the impacts of various controllable analytical factors and uncontrollable practical factors on outcomes of colocalization analysis through realistic simulations and real data examples. We identify a single analytical factor, the specification of prior enrichment levels, which can lead to severe inflation of false-positive colocalization findings. Meanwhile, the combination of many other analytical and practical factors all lead to diminished power. Consequently, we recommend the following strategies for the best practice of colocalization analysis: (1) estimating prior enrichment level from the observed data and (2) separating fine-mapping and colocalization analysis. Our analysis of 4,091 complex traits and the multi-tissue expression quantitative trait loci (eQTL) data from the GTEx (v.8) suggests that colocalizations of molecular QTLs and causal complex trait associations are widespread. However, only a small proportion can be confidently identified from currently available data due to a lack of power. Our findings set a benchmark for current and future integrative genetic association analysis applications.

Effect of culture conditions on conidia production and enhancement of environmental stress resistance of Esteya vermicola in solid-state fermentation

AIMS

Esteya vermicola is an endoparasitic fungus producing lunate conidia, which kill pine wood nematode (PWN), and PWN could cause pine wilt disease (PWD). The aims of this study were to increase production and confirm the resistance (temperature and UV irradiation) of lunate conidia, and further determine the effective concentrations of conidia infecting PWN.

METHODS AND RESULTS

In this study, rice was used as a carrier to absorb conidial suspension to propagate conidia. The optimal conditions for lunate conidia production were 25°C temperature, 9 days of culture time, 2 : 1 rice/distilled water ratio and 10% inoculum size. The germination rate of E. vermicola cultured on potato dextrose agar was influenced by UV irradiation, similar to growth on rice. Esteya vermicola cultured on rice under heat stress might be more suitable for application in the field. The concentration (1 × 10⁸ conidia per ml) to kill PWN had the highest infectivity among the four conidia concentrations tested after 3 days of inoculation.

CONCLUSIONS

This study showed a rice substrate-supported high-quality conidia production and the optimal infectivity concentration of E. vermicola.

SIGNIFICANCE AND IMPACT OF THE STUDY

These results provide the necessary process of an economical and efficient biological control strategy against PWD.