EVPsort: An Atlas of Small ncRNA Profiling and Sorting in Extracellular Vesicles and Particles.

Extracellular vesicles and particles (EVPs) play a crucial role in mediating cell-to-cell communication by transporting various molecular cargos, with small non-coding RNAs (ncRNAs) holding particular significance. A thorough investigation into the abundance and sorting mechanisms of ncRNA within EVPs is imperative for advancing their clinical applications. We have developed EVPsort, which not only provides an extensive overview of ncRNA profiling in 3,162 samples across various biofluids, cell lines, and disease contexts but also seamlessly integrates 19 external databases and tools. This integration encompasses information on associations between ncRNAs and RNA-binding proteins (RBPs), motifs, targets, pathways, diseases, and drugs. With its rich resources and powerful analysis tools, EVPsort extends its profiling capabilities to investigate ncRNA sorting, identify relevant RBPs and motifs, and assess functional implications. EVPsort stands as a pioneering database dedicated to comprehensively addressing both the abundance and sorting of ncRNA within EVPs. It is freely accessible at https://bioinfo.vanderbilt.edu/evpsort/.

FindAdapt: A python package for fast and accurate adapter detection in small RNA sequencing.

Adapter trimming is an essential step for analyzing small RNA sequencing data, where reads are generally longer than target RNAs ranging from 18 to 30 bp. Most adapter trimming tools require adapter information as input. However, adapter information is hard to access, specified incorrectly, or not provided with publicly available datasets, hampering their reproducibility and reusability. Manual identification of adapter patterns from raw reads is labor-intensive and error-prone. Moreover, the use of randomized adapters to reduce ligation biases during library preparation makes adapter detection even more challenging. Here, we present FindAdapt, a Python package for fast and accurate detection of adapter patterns without relying on prior information. We demonstrated that FindAdapt was far superior to existing approaches. It identified adapters successfully in 180 simulation datasets with diverse read structures and 3,184 real datasets covering a variety of commercial and customized small RNA library preparation kits. FindAdapt is stand-alone software that can be easily integrated into small RNA sequencing analysis pipelines.

Probable digenic inheritance of Diamond-Blackfan anemia.

A 26-year-old female proband with a clinical diagnosis and consistent phenotype of Diamond-Blackfan anemia (DBA, OMIM 105650) without an identified genotype was referred to the Undiagnosed Diseases Network. DBA is classically associated with monoallelic variants that have an autosomal-dominant or -recessive mode of inheritance. Intriguingly, her case was solved by a detection of a digenic interaction between non-allelic RPS19 and RPL27 variants. This was confirmed with a machine learning structural model, co-segregation analysis, and RNA sequencing. This is the first report of DBA caused by a digenic effect of two non-allelic variants demonstrated by machine learning structural model. This case suggests that atypical phenotypic presentations of DBA may be caused by digenic inheritance in some individuals. We also conclude that a machine learning structural model can be useful in detecting digenic models of possible interactions between products encoded by alleles of different genes inherited from non-affected carrier parents that can result in DBA with an unrealized 25% recurrence risk.

Mendelian Randomization Analysis of Genetic Proxies of Thiazide Diuretics and the Reduction of Kidney Stone Risk.

Importance: Clinical trial data have called into question the efficacy of thiazide diuretics for the prevention of kidney stones. Objective: To identify whether there is an association between genetic proxies of thiazide diuretics and the risk of kidney stones. Design, Setting, and Participants: This genetic association study undertook a mendelian randomization analysis of derived exposures and outcomes from genome-wide association study summary statistics. Genetic proxies of thiazide diuretics were derived from the International Consortium for Blood Pressure. Kidney stone cases and controls were derived from the Million Veteran Program, UK Biobank, and the FinnGen study. These cross-sectional designs do not report a duration of follow-up. Data analysis was performed in May 2023. Exposure: Genetic proxies of thiazide diuretics were genetic variants in the thiazide-sensitive sodium chloride cotransporter gene associated with systolic blood pressure. Genetic proxies of ß-blockers and systolic blood pressure served as negative controls. Main Outcomes and Measures: The main outcome was the odds of kidney stones. The secondary outcomes were serum laboratory values relevant to the treatment of kidney stones. Results: The main analysis included up to 1 079 657 individuals, including 50 832 kidney stone cases and 1 028 825 controls. In a meta-analysis of all cohorts, genetic proxies of thiazide diuretics were associated with a lower odds of kidney stones (OR, 0.85; 95% CI, 0.81-0.89; P < .001). Genetic proxies of ß-blockers (OR, 1.02; 95% CI, 0.96-1.07; P = .52) and systolic blood pressure (OR, 1.00; 95% CI, 1.00-1.01; P = .49) were not associated with kidney stones. Genetic proxies of thiazide diuretics were associated with higher serum calcium (ß [SE], 0.051 [0.0092]; P < .001) and total cholesterol (ß [SE], 0.065 [0.015]; P < .001), but lower serum potassium (ß [SE], -0.073 [0.022]; P < .001). Conclusions and Relevance: In this genetic association study, genetic proxies of thiazide diuretics were associated with reduced kidney stone risk. This finding reflects a drug effect over the course of a lifetime, unconstrained by the limited follow-up period of clinical trials.

Genetic Inhibition of APOL1 Pore-Forming Function Prevents APOL1-Mediated Kidney Disease.

SIGNIFICANCE STATEMENT: African Americans are at increased risk of CKD in part due to high-risk (HR) variants in the apolipoprotein L1 ( APOL1 ) gene, termed G1/G2. A different APOL1 variant, p.N264K , reduced the risk of CKD and ESKD among carriers of APOL1 HR variants to levels comparable with individuals with APOL1 low-risk variants in an analysis of 121,492 participants of African ancestry from the Million Veteran Program (MVP). Functional genetic studies in cell models showed that APOL1 p.N264K blocked APOL1 pore-forming function and ion channel conduction and reduced toxicity of APOL1 HR mutations. Pharmacologic inhibitors that mimic this mutation blocking APOL1 -mediated pore formation may be able to prevent and/or treat APOL1 -associated kidney disease. BACKGROUND: African Americans are at increased risk for nondiabetic CKD in part due to HR variants in the APOL1 gene. METHODS: We tested whether a different APOL1 variant, p.N264K , modified the association between APOL1 HR genotypes (two copies of G1/G2) and CKD in a cross-sectional analysis of 121,492 participants of African ancestry from the MVP. We replicated our findings in the Vanderbilt University Biobank ( n =14,386) and National Institutes of Health All of Us ( n =14,704). Primary outcome was CKD and secondary outcome was ESKD among nondiabetic patients. Primary analysis compared APOL1 HR genotypes with and without p.N264K . Secondary analyses included APOL1 low-risk genotypes and tested for interaction. In MVP, we performed sequential logistic regression models adjusting for demographics, comorbidities, medications, and ten principal components of ancestry. Functional genomic studies expressed APOL1 HR variants with and without APOL1 p.N264K in cell models. RESULTS: In the MVP cohort, 15,604 (12.8%) had two APOL1 HR variants, of which 582 (0.5%) also had APOL1 p.N264K . In MVP, 18,831 (15%) had CKD, 4177 (3%) had ESKD, and 34% had diabetes. MVP APOL1 HR, without p.N264K , was associated with increased odds of CKD (odds ratio [OR], 1.72; 95% confidence interval [CI], 1.60 to 1.85) and ESKD (OR, 3.94; 95% CI, 3.52 to 4.41). In MVP, APOL1 p.N264K mitigated the renal risk of APOL1 HR, in CKD (OR, 0.43; 95% CI, 0.28 to 0.65) and ESKD (OR, 0.19; CI 0.07 to 0.51). In the replication cohorts meta-analysis, APOL1 p.N264K mitigated the renal risk of APOL1 HR in CKD (OR, 0.40; 95% CI, 0.18 to 0.92) and ESKD (OR, 0.19; 95% CI, 0.05 to 0.79). In the mechanistic studies, APOL1 p.N264K blocked APOL1 pore-forming function and ion channel conduction and reduced toxicity of APOL1 HR variants. CONCLUSIONS: APOL1 p.N264K is associated with reduced risk of CKD and ESKD among carriers of APOL1 HR to levels comparable with individuals with APOL1 low-risk genotypes.

scDemultiplex: An iterative beta-binomial model-based method for accurate demultiplexing with hashtag oligos.

Single-cell sequencing have been widely used to characterize cellular heterogeneity. Sample multiplexing where multiple samples are pooled together for single-cell experiments, attracts wide attention due to its benefits of increasing capacity, reducing costs, and minimizing batch effects. To analyze multiplexed data, the first crucial step is to demultiplex, the process of assigning cells to individual samples. Inaccurate demultiplexing will create false cell types and result in misleading characterization. We propose scDemultiplex, which models hashtag oligo (HTO) counts with beta-binomial distribution and uses an iterative strategy for further refinement. Compared with seven existing demultiplexing approaches, scDemultiplex achieved great performance in both high-quality and low-quality data. Additionally, scDemultiplex can be combined with other approaches to improve their performance.

Systematic Assessment of Small RNA Profiling in Human Extracellular Vesicles.

MOTIVATION: Extracellular vesicles (EVs) are produced and released by most cells and are now recognized to play a role in intercellular communication through the delivery of molecular cargo, including proteins, lipids, and RNA. Small RNA sequencing (small RNA-seq) has been widely used to characterize the small RNA content in EVs. However, there is a lack of a systematic assessment of the quality, technical biases, RNA composition, and RNA biotypes enrichment for small RNA profiling of EVs across cell types, biofluids, and conditions. METHODS: We collected and reanalyzed small RNA-seq datasets for 2756 samples from 83 studies involving 55 with EVs only and 28 with both EVs and matched donor cells. We assessed their quality by the total number of reads after adapter trimming, the overall alignment rate to the host and non-host genomes, and the proportional abundance of total small RNA and specific biotypes, such as miRNA, tRNA, rRNA, and Y RNA. RESULTS: We found that EV extraction methods varied in their reproducibility in isolating small RNAs, with effects on small RNA composition. Comparing proportional abundances of RNA biotypes between EVs and matched donor cells, we discovered that rRNA and tRNA fragments were relatively enriched, but miRNAs and snoRNA were depleted in EVs. Except for the export of eight miRNAs being context-independent, the selective release of most miRNAs into EVs was study-specific. CONCLUSION: This work guides quality control and the selection of EV isolation methods and enhances the interpretation of small RNA contents and preferential loading in EVs.

Phenome-Wide Association Study of UMOD Gene Variants and Differential Associations With Clinical Outcomes Across Populations in the Million Veteran Program a Multiethnic Biobank.

Introduction: Common variants in the UMOD gene are considered an evolutionary adaptation against urinary tract infections (UTIs) and have been implicated in kidney stone formation, chronic kidney disease (CKD), and hypertension. However, differences in UMOD variant-phenotype associations across population groups are unclear. Methods: We tested associations between UMOD/PDILT variants and up to 1528 clinical diagnosis codes mapped to phenotype groups in the Million Veteran Program (MVP), using published phenome-wide association study (PheWAS) methodology. Associations were tested using logistic regression adjusted for age, sex, and 10 principal components of ancestry. Bonferroni correction for multiple comparisons was applied. Results: Among 648,593 veterans, mean (SD) age was 62 (14) years; 9% were female, 19% Black, and 8% Hispanic. In White patients, the rs4293393 UMOD risk variant associated with increased uromodulin was associated with increased odds of CKD (odds ratio [OR]: 1.22, 95% CI: 1.20-1.24, P = 5.90 × 10-111), end-stage kidney disease (OR: 1.17, 95% CI: 1.11-1.24, P = 2.40 × 10-09), and hypertension (OR: 1.03, 95% CI: 1.05-1.05, P = 2.11 × 10-06) and significantly lower odds of UTIs (OR: 0.94, 95% CI: 0.92-0.96, P = 1.21 × 10-10) and kidney calculus (OR: 0.85, 95% CI: 0.83-0.86, P = 4.27 × 10-69). Similar findings were observed across UMOD/PDILT variants. The rs77924615 PDILT variant had stronger associations with acute cystitis in White female (OR: 0.73, 95% CI: 0.59-0.91, P = 4.98 × 10-03) versus male (OR: 0.99, 95% CI: 0.89-1.11, P = 8.80 × 10-01) (P interaction = 0.01) patients. In Black patients, the rs77924615 PDILT variant was significantly associated with pyelonephritis (OR: 0.65, 95% CI: 0.54-0.79, P = 1.05 × 10-05), whereas associations with UMOD promoter variants were attenuated. Conclusion: Robust associations were observed between UMOD/PDILT variants linked with increased uromodulin expression and lower odds of UTIs and calculus and increased odds of CKD and hypertension. However, these associations varied significantly across ancestry groups and sex.

Pharmacogenetic variants and risk of remdesivir-associated liver enzyme elevations in Million Veteran Program participants hospitalized with COVID-19.

Remdesivir is the first US Food and Drug Administration (FDA)-approved drug for the treatment of coronavirus disease 2019 (COVID-19). We conducted a retrospective pharmacogenetic study to examine remdesivir-associated liver enzyme elevation among Million Veteran Program participants hospitalized with COVID-19 between March 15, 2020, and June 30, 2021. Pharmacogene phenotypes were assigned using Stargazer. Linear regression was performed on peak log-transformed enzyme values, stratified by population, adjusted for age, sex, baseline liver enzymes, comorbidities, and 10 population-specific principal components. Patients on remdesivir had higher peak alanine aminotransferase (ALT) values following treatment initiation compared with patients not receiving remdesivir. Remdesivir administration was associated with a 33% and 24% higher peak ALT in non-Hispanic White (NHW) and non-Hispanic Black (NHB) participants (p < 0.001), respectively. In a multivariable model, NHW CYP2C19 intermediate/poor metabolizers had a 9% increased peak ALT compared with NHW normal/rapid/ultrarapid metabolizers (p = 0.015); this association was not observed in NHB participants. In summary, remdesivir-associated ALT elevations appear to be multifactorial, and further studies are needed.

Genetic Variants Associated With Mineral Metabolism Traits in Chronic Kidney Disease.

CONTEXT: Chronic kidney disease (CKD) causes multiple interrelated disturbances in mineral metabolism. Genetic studies in the general population have identified common genetic variants associated with circulating phosphate, calcium, parathyroid hormone (PTH), and fibroblast growth factor 23 (FGF23). OBJECTIVE: In this study we aimed to discover genetic variants associated with circulating mineral markers in CKD. METHODS: We conducted candidate single-nucleotide variation (SNV) analysis in 3027 participants in the multiethnic Chronic Renal Insufficiency Cohort (CRIC) to determine the associations between SNVs and circulating levels of mineral markers. RESULTS: SNVs adjacent to or within genes encoding the regulator of G protein-coupled signaling 14 (RGS14) and the calcium-sensing receptor (CASR) were associated with levels of mineral metabolites. The strongest associations (P < .001) were at rs4074995 (RGS14) for phosphate (0.09 mg/dL lower per minor allele) and FGF23 (8.6% lower), and at rs1801725 (CASR) for calcium (0.12 mg/dL higher). In addition, the prevalence of hyperparathyroidism differed by rs4074995 (RGS14) genotype (chi-square P < .0001). Differential inheritance by race was noted for the minor allele of RGS14. Expression quantitative loci (eQTL) analysis showed that rs4074995 was associated with lower RGS14 gene expression in glomeruli (P = 1.03 × 10-11) and tubules (P = 4.0 × 10-4). CONCLUSION: We evaluated genetic variants associated with mineral metabolism markers in a CKD population. Participants with CKD and the minor allele of rs4074995 (RGS14) had lower phosphorus, lower plasma FGF23, and lower prevalence of hyperparathyroidism. The minor allele of RGS14 was also associated with lower gene expression in the kidney. Further studies are needed to elucidate the effect of rs4074995 on the pathogenesis of disordered mineral metabolism in CKD.

APOL1 Risk Variants, Acute Kidney Injury, and Death in Participants With African Ancestry Hospitalized With COVID-19 From the Million Veteran Program.

IMPORTANCE: Coronavirus disease 2019 (COVID-19) confers significant risk of acute kidney injury (AKI). Patients with COVID-19 with AKI have high mortality rates. OBJECTIVE: Individuals with African ancestry with 2 copies of apolipoprotein L1 (APOL1) variants G1 or G2 (high-risk group) have significantly increased rates of kidney disease. We tested the hypothesis that the APOL1 high-risk group is associated with a higher-risk of COVID-19-associated AKI and death. DESIGN, SETTING, AND PARTICIPANTS: This retrospective cohort study included 990 participants with African ancestry enrolled in the Million Veteran Program who were hospitalized with COVID-19 between March 2020 and January 2021 with available genetic information. EXPOSURES: The primary exposure was having 2 APOL1 risk variants (RV) (APOL1 high-risk group), compared with having 1 or 0 risk variants (APOL1 low-risk group). MAIN OUTCOMES AND MEASURES: The primary outcome was AKI. The secondary outcomes were stages of AKI severity and death. Multivariable logistic regression analyses adjusted for preexisting comorbidities, medications, and inpatient AKI risk factors; 10 principal components of ancestry were performed to study these associations. We performed a subgroup analysis in individuals with normal kidney function prior to hospitalization (estimated glomerular filtration rate ≥60 mL/min/1.73 m2). RESULTS: Of the 990 participants with African ancestry, 905 (91.4%) were male with a median (IQR) age of 68 (60-73) years. Overall, 392 (39.6%) patients developed AKI, 141 (14%) developed stages 2 or 3 AKI, 28 (3%) required dialysis, and 122 (12.3%) died. One hundred twenty-five (12.6%) of the participants were in the APOL1 high-risk group. Patients categorized as APOL1 high-risk group had significantly higher odds of AKI (adjusted odds ratio [OR], 1.95; 95% CI, 1.27-3.02; P = .002), higher AKI severity stages (OR, 2.03; 95% CI, 1.37-2.99; P < .001), and death (OR, 2.15; 95% CI, 1.22-3.72; P = .007). The association with AKI persisted in the subgroup with normal kidney function (OR, 1.93; 95% CI, 1.15-3.26; P = .01). Data analysis was conducted between February 2021 and April 2021. CONCLUSIONS AND RELEVANCE: In this cohort study of veterans with African ancestry hospitalized with COVID-19 infection, APOL1 kidney risk variants were associated with higher odds of AKI, AKI severity, and death, even among individuals with prior normal kidney function.

APOL1 risk variants in individuals of African genetic ancestry drive endothelial cell defects that exacerbate sepsis.

The incidence and severity of sepsis is higher among individuals of African versus European ancestry. We found that genetic risk variants (RVs) in the trypanolytic factor apolipoprotein L1 (APOL1), present only in individuals of African ancestry, were associated with increased sepsis incidence and severity. Serum APOL1 levels correlated with sepsis and COVID-19 severity, and single-cell sequencing in human kidneys revealed high expression of APOL1 in endothelial cells. Analysis of mice with endothelial-specific expression of RV APOL1 and in vitro studies demonstrated that RV APOL1 interfered with mitophagy, leading to cytosolic release of mitochondrial DNA and activation of the inflammasome (NLRP3) and the cytosolic nucleotide sensing pathways (STING). Genetic deletion or pharmacological inhibition of NLRP3 and STING protected mice from RV APOL1-induced permeability defects and proinflammatory endothelial changes in sepsis. Our studies identify the inflammasome and STING pathways as potential targets to reduce APOL1-associated health disparities in sepsis and COVID-19.

A domain damage index to prioritizing the pathogenicity of missense variants.

Prioritizing causal variants is one major challenge for the clinical application of sequencing data. Prompted by the observation that 74.3% of missense pathogenic variants locate in protein domains, we developed an approach named domain damage index (DDI). DDI identifies protein domains depleted of rare missense variations in the general population, which can be further used as a metric to prioritize variants. DDI is significantly correlated with phylogenetic conservation, variant-level metrics, and reported pathogenicity. DDI achieved great performance for distinguishing pathogenic variants from benign ones in three benchmark datasets. The combination of DDI with the other two best approaches improved the performance of each individual method considerably, suggesting DDI provides a powerful and complementary way of variant prioritization.

Beyond Programmed Death-Ligand 1: B7-H6 Emerges as a Potential Immunotherapy Target in SCLC.

INTRODUCTION: The programmed death-ligand 1 (PD-L1) immune checkpoint inhibitors, atezolizumab and durvalumab, have received regulatory approval for the first-line treatment of patients with extensive-stage SCLC. Nevertheless, when used in combination with platinum-based chemotherapy, these PD-L1 inhibitors only improve overall survival by 2 to 3 months. This may be due to the observation that less than 20% of SCLC tumors express PD-L1 at greater than 1%. Evaluating the composition and abundance of checkpoint molecules in SCLC may identify molecules beyond PD-L1 that are amenable to therapeutic targeting. METHODS: We analyzed RNA-sequencing data from SCLC cell lines (n = 108) and primary tumor specimens (n = 81) for expression of 39 functionally validated inhibitory checkpoint ligands. Furthermore, we generated tissue microarrays containing SCLC cell lines and patient with SCLC specimens to confirm expression of these molecules by immunohistochemistry. We annotated patient outcomes data, including treatment response and overall survival. RESULTS: The checkpoint protein B7-H6 (NCR3LG1) exhibited increased protein expression relative to PD-L1 in cell lines and tumors (p < 0.05). Higher B7-H6 protein expression correlated with longer progression-free survival (p = 0.0368) and increased total immune infiltrates (CD45+) in patients. Furthermore, increased B7-H6 gene expression in SCLC tumors correlated with a decreased activated natural killer cell gene signature, suggesting a complex interplay between B7-H6 expression and immune signature in SCLC. CONCLUSIONS: We investigated 39 inhibitory checkpoint molecules in SCLC and found that B7-H6 is highly expressed and associated with progression-free survival. In addition, 26 of 39 immune checkpoint proteins in SCLC tumors were more abundantly expressed than PD-L1, indicating an urgent need to investigate additional checkpoint targets for therapy in addition to PD-L1.

A SAC Phosphoinositide Phosphatase Controls Rice Development via Hydrolyzing PI4P and PI(4,5)P2.

Phosphoinositides (PIs) as regulatory membrane lipids play essential roles in multiple cellular processes. Although the exact molecular targets of PI-dependent modulation remain largely elusive, the effects of disturbed PI metabolism could be employed to identify regulatory modules associated with particular downstream targets of PIs. Here, we identified the role of GRAIN NUMBER AND PLANT HEIGHT1 (GH1), which encodes a suppressor of actin (SAC) domain-containing phosphatase with unknown function in rice (Oryza sativa). Endoplasmic reticulum-localized GH1 specifically dephosphorylated and hydrolyzed phosphatidylinositol 4-phosphate (PI4P) and phosphatidylinositol 4,5-bisphosphate [PI(4,5)P2]. Inactivation of GH1 resulted in massive accumulation of both PI4P and PI(4,5)P2, while excessive GH1 caused their depletion. Notably, superabundant PI4P and PI(4,5)P2 could both disrupt actin cytoskeleton organization and suppress cell elongation. Interestingly, both PI4P and PI(4,5)P2 inhibited actin-related protein2 and -3 (Arp2/3) complex-nucleated actin-branching networks in vitro, whereas PI(4,5)P2 showed more dramatic effects in a dose-dependent manner. Overall, the overaccumulation of PI(4,5)P2 resulting from dysfunction of SAC phosphatase possibly perturbs Arp2/3 complex-mediated actin polymerization, thereby disordering cell development. These findings imply that the Arp2/3 complex might be the potential molecular target of PI(4,5)P2-dependent modulation in eukaryotes, thereby providing insights into the relationship between PI homeostasis and plant growth and development.

Natural alleles of a proteasome α2 subunit gene contribute to thermotolerance and adaptation of African rice.

Global warming threatens many aspects of human life, for example, by reducing crop yields. Breeding heat-tolerant crops using genes conferring thermotolerance is a fundamental way to help deal with this challenge. Here we identify a major quantitative trait locus (QTL) for thermotolerance in African rice (Oryza glaberrima), Thermo-tolerance 1 (TT1), which encodes an α2 subunit of the 26S proteasome involved in the degradation of ubiquitinated proteins. Ubiquitylome analysis indicated that OgTT1 protects cells from heat stress through more efficient elimination of cytotoxic denatured proteins and more effective maintenance of heat-response processes than achieved with OsTT1. Variation in TT1 has been selected for on the basis of climatic temperature and has had an important role in local adaptation during rice evolution. In addition, we found that overexpression of OgTT1 was associated with markedly enhanced thermotolerance in rice, Arabidopsis and Festuca elata. This discovery may lead to an increase in crop security in the face of the ongoing threat of global warming.