Global, regional, and national incidence and mortality burden of non-COVID-19 lower respiratory infections and aetiologies, 1990-2021: a systematic analysis from the Global Burden of Disease Study 2021

BACKGROUND

Lower respiratory infections (LRIs) are a major global contributor to morbidity and mortality. In 2020-21, non-pharmaceutical interventions associated with the COVID-19 pandemic reduced not only the transmission of SARS-CoV-2, but also the transmission of other LRI pathogens. Tracking LRI incidence and mortality, as well as the pathogens responsible, can guide health-system responses and funding priorities to reduce future burden. We present estimates from the Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2021 of the burden of non-COVID-19 LRIs and corresponding aetiologies from 1990 to 2021, inclusive of pandemic effects on the incidence and mortality of select respiratory viruses, globally, regionally, and for 204 countries and territories.

METHODS

We estimated mortality, incidence, and aetiology attribution for LRI, defined by the GBD as pneumonia or bronchiolitis, not inclusive of COVID-19. We analysed 26 259 site-years of mortality data using the Cause of Death Ensemble model to estimate LRI mortality rates. We analysed all available age-specific and sex-specific data sources, including published literature identified by a systematic review, as well as household surveys, hospital admissions, health insurance claims, and LRI mortality estimates, to generate internally consistent estimates of incidence and prevalence using DisMod-MR 2.1. For aetiology estimation, we analysed multiple causes of death, vital registration, hospital discharge, microbial laboratory, and literature data using a network analysis model to produce the proportion of LRI deaths and episodes attributable to the following pathogens: Acinetobacter baumannii, Chlamydia spp, Enterobacter spp, Escherichia coli, fungi, group B streptococcus, Haemophilus influenzae, influenza viruses, Klebsiella pneumoniae, Legionella spp, Mycoplasma spp, polymicrobial infections, Pseudomonas aeruginosa, respiratory syncytial virus (RSV), Staphylococcus aureus, Streptococcus pneumoniae, and other viruses (ie, the aggregate of all viruses studied except influenza and RSV), as well as a residual category of other bacterial pathogens.

FINDINGS

Globally, in 2021, we estimated 344 million (95% uncertainty interval [UI] 325-364) incident episodes of LRI, or 4350 episodes (4120-4610) per 100 000 population, and 2·18 million deaths (1·98-2·36), or 27·7 deaths (25·1-29·9) per 100 000. 502 000 deaths (406 000-611 000) were in children younger than 5 years, among which 254 000 deaths (197 000-320 000) occurred in countries with a low Socio-demographic Index. Of the 18 modelled pathogen categories in 2021, S pneumoniae was responsible for the highest proportions of LRI episodes and deaths, with an estimated 97·9 million (92·1-104·0) episodes and 505 000 deaths (454 000-555 000) globally. The pathogens responsible for the second and third highest episode counts globally were other viral aetiologies (46·4 million [43·6-49·3] episodes) and Mycoplasma spp (25·3 million [23·5-27·2]), while those responsible for the second and third highest death counts were S aureus (424 000 [380 000-459 000]) and K pneumoniae (176 000 [158 000-194 000]). From 1990 to 2019, the global all-age non-COVID-19 LRI mortality rate declined by 41·7% (35·9-46·9), from 56·5 deaths (51·3-61·9) to 32·9 deaths (29·9-35·4) per 100 000. From 2019 to 2021, during the COVID-19 pandemic and implementation of associated non-pharmaceutical interventions, we estimated a 16·0% (13·1-18·6) decline in the global all-age non-COVID-19 LRI mortality rate, largely accounted for by a 71·8% (63·8-78·9) decline in the number of influenza deaths and a 66·7% (56·6-75·3) decline in the number of RSV deaths.

INTERPRETATION

Substantial progress has been made in reducing LRI mortality, but the burden remains high, especially in low-income and middle-income countries. During the COVID-19 pandemic, with its associated non-pharmaceutical interventions, global incident LRI cases and mortality attributable to influenza and RSV declined substantially. Expanding access to health-care services and vaccines, including S pneumoniae, H influenzae type B, and novel RSV vaccines, along with new low-cost interventions against S aureus, could mitigate the LRI burden and prevent transmission of LRI-causing pathogens.

FUNDING

Bill & Melinda Gates Foundation, Wellcome Trust, and Department of Health and Social Care (UK).

Global, regional, and national burden of disorders affecting the nervous system, 1990-2021: a systematic analysis for the Global Burden of Disease Study 2021

BACKGROUND

Disorders affecting the nervous system are diverse and include neurodevelopmental disorders, late-life neurodegeneration, and newly emergent conditions, such as cognitive impairment following COVID-19. Previous publications from the Global Burden of Disease, Injuries, and Risk Factor Study estimated the burden of 15 neurological conditions in 2015 and 2016, but these analyses did not include neurodevelopmental disorders, as defined by the International Classification of Diseases (ICD)-11, or a subset of cases of congenital, neonatal, and infectious conditions that cause neurological damage. Here, we estimate nervous system health loss caused by 37 unique conditions and their associated risk factors globally, regionally, and nationally from 1990 to 2021.

METHODS

We estimated mortality, prevalence, years lived with disability (YLDs), years of life lost (YLLs), and disability-adjusted life-years (DALYs), with corresponding 95% uncertainty intervals (UIs), by age and sex in 204 countries and territories, from 1990 to 2021. We included morbidity and deaths due to neurological conditions, for which health loss is directly due to damage to the CNS or peripheral nervous system. We also isolated neurological health loss from conditions for which nervous system morbidity is a consequence, but not the primary feature, including a subset of congenital conditions (ie, chromosomal anomalies and congenital birth defects), neonatal conditions (ie, jaundice, preterm birth, and sepsis), infectious diseases (ie, COVID-19, cystic echinococcosis, malaria, syphilis, and Zika virus disease), and diabetic neuropathy. By conducting a sequela-level analysis of the health outcomes for these conditions, only cases where nervous system damage occurred were included, and YLDs were recalculated to isolate the non-fatal burden directly attributable to nervous system health loss. A comorbidity correction was used to calculate total prevalence of all conditions that affect the nervous system combined.

FINDINGS

Globally, the 37 conditions affecting the nervous system were collectively ranked as the leading group cause of DALYs in 2021 (443 million, 95% UI 378-521), affecting 3·40 billion (3·20-3·62) individuals (43·1%, 40·5-45·9 of the global population); global DALY counts attributed to these conditions increased by 18·2% (8·7-26·7) between 1990 and 2021. Age-standardised rates of deaths per 100 000 people attributed to these conditions decreased from 1990 to 2021 by 33·6% (27·6-38·8), and age-standardised rates of DALYs attributed to these conditions decreased by 27·0% (21·5-32·4). Age-standardised prevalence was almost stable, with a change of 1·5% (0·7-2·4). The ten conditions with the highest age-standardised DALYs in 2021 were stroke, neonatal encephalopathy, migraine, Alzheimer's disease and other dementias, diabetic neuropathy, meningitis, epilepsy, neurological complications due to preterm birth, autism spectrum disorder, and nervous system cancer.

INTERPRETATION

As the leading cause of overall disease burden in the world, with increasing global DALY counts, effective prevention, treatment, and rehabilitation strategies for disorders affecting the nervous system are needed.

FUNDING

Bill & Melinda Gates Foundation.

Global age-sex-specific mortality, life expectancy, and population estimates in 204 countries and territories and 811 subnational locations, 1950-2021, and the impact of the COVID-19 pandemic: a comprehensive demographic analysis for the Global Burden of Disease Study 2021

BACKGROUND

Estimates of demographic metrics are crucial to assess levels and trends of population health outcomes. The profound impact of the COVID-19 pandemic on populations worldwide has underscored the need for timely estimates to understand this unprecedented event within the context of long-term population health trends. The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) 2021 provides new demographic estimates for 204 countries and territories and 811 additional subnational locations from 1950 to 2021, with a particular emphasis on changes in mortality and life expectancy that occurred during the 2020-21 COVID-19 pandemic period.

METHODS

22 223 data sources from vital registration, sample registration, surveys, censuses, and other sources were used to estimate mortality, with a subset of these sources used exclusively to estimate excess mortality due to the COVID-19 pandemic. 2026 data sources were used for population estimation. Additional sources were used to estimate migration; the effects of the HIV epidemic; and demographic discontinuities due to conflicts, famines, natural disasters, and pandemics, which are used as inputs for estimating mortality and population. Spatiotemporal Gaussian process regression (ST-GPR) was used to generate under-5 mortality rates, which synthesised 30 763 location-years of vital registration and sample registration data, 1365 surveys and censuses, and 80 other sources. ST-GPR was also used to estimate adult mortality (between ages 15 and 59 years) based on information from 31 642 location-years of vital registration and sample registration data, 355 surveys and censuses, and 24 other sources. Estimates of child and adult mortality rates were then used to generate life tables with a relational model life table system. For countries with large HIV epidemics, life tables were adjusted using independent estimates of HIV-specific mortality generated via an epidemiological analysis of HIV prevalence surveys, antenatal clinic serosurveillance, and other data sources. Excess mortality due to the COVID-19 pandemic in 2020 and 2021 was determined by subtracting observed all-cause mortality (adjusted for late registration and mortality anomalies) from the mortality expected in the absence of the pandemic. Expected mortality was calculated based on historical trends using an ensemble of models. In location-years where all-cause mortality data were unavailable, we estimated excess mortality rates using a regression model with covariates pertaining to the pandemic. Population size was computed using a Bayesian hierarchical cohort component model. Life expectancy was calculated using age-specific mortality rates and standard demographic methods. Uncertainty intervals (UIs) were calculated for every metric using the 25th and 975th ordered values from a 1000-draw posterior distribution.

FINDINGS

Global all-cause mortality followed two distinct patterns over the study period: age-standardised mortality rates declined between 1950 and 2019 (a 62·8% [95% UI 60·5-65·1] decline), and increased during the COVID-19 pandemic period (2020-21; 5·1% [0·9-9·6] increase). In contrast with the overall reverse in mortality trends during the pandemic period, child mortality continued to decline, with 4·66 million (3·98-5·50) global deaths in children younger than 5 years in 2021 compared with 5·21 million (4·50-6·01) in 2019. An estimated 131 million (126-137) people died globally from all causes in 2020 and 2021 combined, of which 15·9 million (14·7-17·2) were due to the COVID-19 pandemic (measured by excess mortality, which includes deaths directly due to SARS-CoV-2 infection and those indirectly due to other social, economic, or behavioural changes associated with the pandemic). Excess mortality rates exceeded 150 deaths per 100 000 population during at least one year of the pandemic in 80 countries and territories, whereas 20 nations had a negative excess mortality rate in 2020 or 2021, indicating that all-cause mortality in these countries was lower during the pandemic than expected based on historical trends. Between 1950 and 2021, global life expectancy at birth increased by 22·7 years (20·8-24·8), from 49·0 years (46·7-51·3) to 71·7 years (70·9-72·5). Global life expectancy at birth declined by 1·6 years (1·0-2·2) between 2019 and 2021, reversing historical trends. An increase in life expectancy was only observed in 32 (15·7%) of 204 countries and territories between 2019 and 2021. The global population reached 7·89 billion (7·67-8·13) people in 2021, by which time 56 of 204 countries and territories had peaked and subsequently populations have declined. The largest proportion of population growth between 2020 and 2021 was in sub-Saharan Africa (39·5% [28·4-52·7]) and south Asia (26·3% [9·0-44·7]). From 2000 to 2021, the ratio of the population aged 65 years and older to the population aged younger than 15 years increased in 188 (92·2%) of 204 nations.

INTERPRETATION

Global adult mortality rates markedly increased during the COVID-19 pandemic in 2020 and 2021, reversing past decreasing trends, while child mortality rates continued to decline, albeit more slowly than in earlier years. Although COVID-19 had a substantial impact on many demographic indicators during the first 2 years of the pandemic, overall global health progress over the 72 years evaluated has been profound, with considerable improvements in mortality and life expectancy. Additionally, we observed a deceleration of global population growth since 2017, despite steady or increasing growth in lower-income countries, combined with a continued global shift of population age structures towards older ages. These demographic changes will likely present future challenges to health systems, economies, and societies. The comprehensive demographic estimates reported here will enable researchers, policy makers, health practitioners, and other key stakeholders to better understand and address the profound changes that have occurred in the global health landscape following the first 2 years of the COVID-19 pandemic, and longer-term trends beyond the pandemic.

FUNDING

Bill & Melinda Gates Foundation.

Integrated genomic network analysis revealed potential of a druggable target for hemorrhoid treatment

Hemorrhoids are a prevalent medical condition that necessitates effective treatment options. The current options for treatment consist of oral medications, topical applications, or surgery, yet a scarcity of highly effective drugs still exists. Genetic markers provide promising avenues for investigating the treatment of hemorrhoids, as they may reveal intricate biological mechanisms and targeted drug therapies, ultimately enhancing more precise treatment tailored to the patient. This study aims to identify new drug candidates for treating hemorrhoids through a meticulous bioinformatics approach and integrated with genomic network analysis. After extracting 21 druggable target genes using DrugBank from 293 genes connected to hemorrhoids, 87 possible drugs were selected. Three of these drugs (ketamine, methylene blue, and fulvestrant) hold potential in addressing issues associated with hemorrhoids and have been supported by clinical or preclinical studies. Eighty-four compounds present new therapeutic possibilities for managing hemorrhoids. We highlight that our findings indicate that NOX1 and NOS3 genes are promising biomarkers, with NOS3 gaining significance owing to its robust systemic functional annotations. Sapropterin, an existing drug, is closely associated with NOS3, providing a clear target for biomarker-driven interventions. This study illustrates the potential of combining genomic network analysis with bioinformatics to repurpose drugs for treating hemorrhoids. Subsequent research will explore the mechanisms for utilizing NOS3 targeting in the treatment of hemorrhoids.

A genomic and bioinformatic-based approach to identify genetic variants for liver cancer across multiple continents

Liver cancer is the fourth leading cause of death worldwide. Well-known risk factors include hepatitis B virus and hepatitis C virus, along with exposure to aflatoxins, excessive alcohol consumption, obesity, and type 2 diabetes. Genomic variants play a crucial role in mediating the associations between these risk factors and liver cancer. However, the specific variants involved in this process remain under-explored. This study utilized a bioinformatics approach to identify genetic variants associated with liver cancer from various continents. Single-nucleotide polymorphisms associated with liver cancer were retrieved from the genome-wide association studies catalog. Prioritization was then performed using functional annotation with HaploReg v4.1 and the Ensembl database. The prevalence and allele frequencies of each variant were evaluated using Pearson correlation coefficients. Two variants, rs2294915 and rs2896019, encoded by the PNPLA3 gene, were found to be highly expressed in the liver tissue, as well as in the skin, cell-cultured fibroblasts, and adipose-subcutaneous tissue, all of which contribute to the risk of liver cancer. We further found that these two SNPs (rs2294915 and rs2896019) were positively correlated with the prevalence rate. Positive associations with the prevalence rate were more frequent in East Asian and African populations. We highlight the utility of this population-specific PNPLA3 genetic variant for genetic association studies and for the early prognosis and treatment of liver cancer. This study highlights the potential of integrating genomic databases with bioinformatic analysis to identify genetic variations involved in the pathogenesis of liver cancer. The genetic variants investigated in this study are likely to predispose to liver cancer and could affect its progression and aggressiveness. We recommend future research prioritizing the validation of these variations in clinical settings.

Identification of druggable genes for multiple myeloma based on genomic information

Multiple myeloma (MM) is a hematological malignancy. It is widely believed that genetic factors play a significant role in the development of MM, as investigated in numerous studies. However, the application of genomic information for clinical purposes, including diagnostic and prognostic biomarkers, remains largely confined to research. In this study, we utilized genetic information from the Genomic-Driven Clinical Implementation for Multiple Myeloma database, which is dedicated to clinical trial studies on MM. This genetic information was sourced from the genome-wide association studies catalog database. We prioritized genes with the potential to cause MM based on established annotations, as well as biological risk genes for MM, as potential drug target candidates. The DrugBank database was employed to identify drug candidates targeting these genes. Our research led to the discovery of 14 MM biological risk genes and the identification of 10 drugs that target three of these genes. Notably, only one of these 10 drugs, panobinostat, has been approved for use in MM. The two most promising genes, calcium signal-modulating cyclophilin ligand (CAMLG) and histone deacetylase 2 (HDAC2), were targeted by four drugs (cyclosporine, belinostat, vorinostat, and romidepsin), all of which have clinical evidence supporting their use in the treatment of MM. Interestingly, five of the 10 drugs have been approved for other indications than MM, but they may also be effective in treating MM. Therefore, this study aimed to clarify the genomic variants involved in the pathogenesis of MM and highlight the potential benefits of these genomic variants in drug discovery.

Leveraging Genomic and Bioinformatic Analysis to Enhance Drug Repositioning for Dermatomyositis

Dermatomyositis (DM) is an autoimmune disease that is classified as a type of idiopathic inflammatory myopathy, which affects human skin and muscles. The most common clinical symptoms of DM are muscle weakness, rash, and scaly skin. There is currently no cure for DM. Genetic factors are known to play a pivotal role in DM progression, but few have utilized this information geared toward drug discovery for the disease. Here, we exploited genomic variation associated with DM and integrated this with genomic and bioinformatic analyses to discover new drug candidates. We first integrated genome-wide association study (GWAS) and phenome-wide association study (PheWAS) catalogs to identify disease-associated genomic variants. Biological risk genes for DM were prioritized using strict functional annotations, further identifying candidate drug targets based on druggable genes from databases. Overall, we analyzed 1239 variants associated with DM and obtained 43 drugs that overlapped with 13 target genes (JAK2, FCGR3B, CD4, CD3D, LCK, CD2, CD3E, FCGR3A, CD3G, IFNAR1, CD247, JAK1, IFNAR2). Six drugs clinically investigated for DM, as well as eight drugs under pre-clinical investigation, are candidate drugs that could be repositioned for DM. Further studies are necessary to validate potential biomarkers for novel DM therapeutics from our findings.

Transcriptomics-driven drug repositioning for the treatment of diabetic foot ulcer

Diabetic foot ulcers (DFUs) are a common complication of diabetes and can lead to severe disability and even amputation. Despite advances in treatment, there is currently no cure for DFUs and available drugs for treatment are limited. This study aimed to identify new candidate drugs and repurpose existing drugs to treat DFUs based on transcriptomics analysis. A total of 31 differentially expressed genes (DEGs) were identified and used to prioritize the biological risk genes for DFUs. Further investigation using the database DGIdb revealed 12 druggable target genes among 50 biological DFU risk genes, corresponding to 31 drugs. Interestingly, we highlighted that two drugs (urokinase and lidocaine) are under clinical investigation for DFU and 29 drugs are potential candidates to be repurposed for DFU therapy. The top 5 potential biomarkers for DFU from our findings are IL6ST, CXCL9, IL1R1, CXCR2, and IL10. This study highlights IL1R1 as a highly promising biomarker for DFU due to its high systemic score in functional annotations, that can be targeted with an existing drug, Anakinra. Our study proposed that the integration of transcriptomic and bioinformatic-based approaches has the potential to drive drug repurposing for DFUs. Further research will further examine the mechanisms by which targeting IL1R1 can be used to treat DFU.

A bioinformatic approach to identify pathogenic variants for Stevens-Johnson syndrome

Stevens-Johnson syndrome (SJS) produces a severe hypersensitivity reaction caused by Herpes simplex virus or mycoplasma infection, vaccination, systemic disease, or other agents. Several studies have investigated the genetic susceptibility involved in SJS. To provide further genetic insights into the pathogenesis of SJS, this study prioritized high-impact, SJS-associated pathogenic variants through integrating bioinformatic and population genetic data. First, we identified SJS-associated single nucleotide polymorphisms from the genome-wide association studies catalog, followed by genome annotation with HaploReg and variant validation with Ensembl. Subsequently, expression quantitative trait locus (eQTL) from GTEx identified human genetic variants with differential gene expression across human tissues. Our results indicate that two variants, namely rs2074494 and rs5010528, which are encoded by the HLA-C (human leukocyte antigen C) gene, were found to be differentially expressed in skin. The allele frequencies for rs2074494 and rs5010528 also appear to significantly differ across continents. We highlight the utility of these population-specific HLA-C genetic variants for genetic association studies, and aid in early prognosis and disease treatment of SJS.

Bioinformatics Analysis to Uncover the Potential Drug Targets Responsible for Mycobacterium tuberculosis Peptidoglycan and Lysine Biosynthesis

Drug-resistant tuberculosis (TB), which results mainly from the selection of naturally resistant strains of Mycobacterium tuberculosis (MTB) due to mismanaged treatment, poses a severe challenge to the global control of TB. Therefore, screening novel and unique drug targets against this pathogen is urgently needed. The metabolic pathways of Homo sapiens and MTB were compared using the Kyoto Encyclopedia of Genes and Genomes tool, and further, the proteins that are involved in the metabolic pathways of MTB were subtracted and proceeded to protein-protein interaction network analysis, subcellular localization, drug ability testing, and gene ontology. The study aims to identify enzymes for the unique pathways for further screening to determine the feasibility of the therapeutic targets. The qualitative characteristics of 28 proteins identified as drug target candidates were studied. The results showed that 12 were cytoplasmic, 2 were extracellular, 12 were transmembrane, and 3 were unknown. Furthermore, druggability analysis revealed 14 druggable proteins, of which 12 were novel and responsible for MTB peptidoglycan and lysine biosynthesis. The novel targets obtained in this study are used to develop antimicrobial treatments against pathogenic bacteria. Future studies should further shed light on the clinical implementation to identify antimicrobial therapies against MTB.

Investigation of susceptibility genes for chickenpox disease across multiple continents

Chickenpox (varicella) is caused by infection with the varicella-zoster virus (VZV), a neurotropic alpha herpes virus with a double-stranded DNA genome. Chickenpox can cause life-threatening complications, including subsequent bacterial infections, central nervous system symptoms, and even death without any risk factors. Few studies have been reported to investigate genetic susceptibility implicated in chickenpox. Herein, our study identified global genetic variants that potentially contributed to chickenpox susceptibility by utilizing the established bioinformatic-based approach. We integrated several databases, such as genome-wide association studies (GWAS) catalog, GTEx portal, HaploReg version 4.1, and Ensembl databases analyses to investigate susceptibility genes associated with chickenpox. Notably, increased expression of HLA-S, HCG4P5, and ABHD16A genes underlie enhanced chickenpox susceptibility in the European, American, and African populations. As compared to the Asian population, Europeans, Americans, and Africans have higher allele frequencies of the extant variants rs9266089, rs10947050, and rs79501286 from the susceptibility genes. Our study suggested that these susceptibility genes and associated genetic variants might play a critical role in chickenpox progression based on host genetics with clinical implications.

Integrated genomic analysis to identify druggable targets for pancreatic cancer

According to the National Comprehensive Cancer Network and the American Society of Clinical Oncology, the standard treatment for pancreatic cancer (PC) is gemcitabine and fluorouracil. Other chemotherapeutic agents have been widely combined. However, drug resistance remains a huge challenge, leading to the ineffectiveness of cancer therapy. Therefore, we are trying to discover new treatments for PC by utilizing genomic information to identify PC-associated genes as well as drug target genes for drug repurposing. Genomic information from a public database, the cBio Cancer Genomics Portal, was employed to retrieve the somatic mutation genes of PC. Five functional annotations were applied to prioritize the PC risk genes: Kyoto Encyclopedia of Genes and Genomes; biological process; knockout mouse; Gene List Automatically Derived For You; and Gene Expression Omnibus Dataset. DrugBank database was utilized to extract PC drug targets. To narrow down the most promising drugs for PC, CMap Touchstone analysis was applied. Finally, ClinicalTrials.gov and a literature review were used to screen the potential drugs under clinical and preclinical investigation. Here, we extracted 895 PC-associated genes according to the cBioPortal database and prioritized them by using five functional annotations; 318 genes were assigned as biological PC risk genes. Further, 216 genes were druggable according to the DrugBank database. CMap Touchstone analysis indicated 13 candidate drugs for PC. Among those 13 drugs, 8 drugs are in the clinical trials, 2 drugs were supported by the preclinical studies, and 3 drugs are with no evidence status for PC. Importantly, we found that midostaurin (targeted PRKA) and fulvestrant (targeted ESR1) are promising candidate drugs for PC treatment based on the genomic-driven drug repurposing pipelines. In short, integrated analysis using a genomic information database demonstrated the viability for drug repurposing. We proposed two drugs (midostaurin and fulvestrant) as promising drugs for PC.

Integration of genomic variants and bioinformatic-based approach to drive drug repurposing for multiple sclerosis

Multiple sclerosis (MS) is a chronic autoimmune disease in the central nervous system (CNS) marked by inflammation, demyelination, and axonal loss. Currently available MS medication is limited, thereby calling for a strategy to accelerate new drug discovery. One of the strategies to discover new drugs is to utilize old drugs for new indications, an approach known as drug repurposing. Herein, we first identified 421 MS-associated SNPs from the Genome-Wide Association Study (GWAS) catalog (p-value < 5 × 10-8), and a total of 427 risk genes associated with MS using HaploReg version 4.1 under the criterion r 2  > 0.8. MS risk genes were then prioritized using bioinformatics analysis to identify biological MS risk genes. The prioritization was performed based on six defined categories of functional annotations, namely missense mutation, cis-expression quantitative trait locus (cis-eQTL), molecular pathway analysis, protein-protein interaction (PPI), genes overlap with knockout mouse phenotype, and primary immunodeficiency (PID). A total of 144 biological MS risk genes were found and mapped into 194 genes within an expanded PPI network. According to the DrugBank and the Therapeutic Target Database, 27 genes within the list targeted by 68 new candidate drugs were identified. Importantly, the power of our approach is confirmed with the identification of a known approved drug (dimethyl fumarate) for MS. Based on additional data from ClinicalTrials.gov, eight drugs targeting eight distinct genes are prioritized with clinical evidence for MS disease treatment. Notably, CD80 and CD86 pathways are promising targets for MS drug repurposing. Using in silico drug repurposing, we identified belatacept as a promising MS drug candidate. Overall, this study emphasized the integration of functional genomic variants and bioinformatic-based approach that reveal important biological insights for MS and drive drug repurposing efforts for the treatment of this devastating disease.

Identification of Hub Genes and Potential Biomarkers for Childhood Asthma by Utilizing an Established Bioinformatic Analysis Approach

Childhood asthma represents a heterogeneous disease resulting from the interaction between genetic factors and environmental exposures. Currently, finding reliable biomarkers is necessary for the clinical management of childhood asthma. However, only a few biomarkers are being used in clinical practice in the pediatric population. In the long run, new biomarkers for asthma in children are required and would help direct therapy approaches. This study aims to identify potential childhood asthma biomarkers using a genetic-driven biomarkers approach. Herein, childhood asthma-associated Single Nucleotide Polymorphisms (SNPs) were utilized from the GWAS database to drive and facilitate the biomarker of childhood asthma. We uncovered 466 childhood asthma-associated loci by extending to proximal SNPs based on r2 > 0.8 in Asian populations and utilizing HaploReg version 4.1 to determine 393 childhood asthma risk genes. Next, the functional roles of these genes were subsequently investigated using Gene Ontology (GO) term enrichment analysis, a Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, and a protein−protein interaction (PPI) network. MCODE and CytoHubba are two Cytoscape plugins utilized to find biomarker genes from functional networks created using childhood asthma risk genes. Intriguingly, 10 hub genes (IL6, IL4, IL2, IL13, PTPRC, IL5, IL33, TBX21, IL2RA, and STAT6) were successfully identified and may have been identified to play a potential role in the pathogenesis of childhood asthma. Among 10 hub genes, we strongly suggest IL6 and IL4 as prospective childhood asthma biomarkers since both of these biomarkers achieved a high systemic score in Cytohubba’s MCC algorithm. In summary, this study offers a valuable genetic-driven biomarker approach to facilitate the potential biomarkers for asthma in children.

The use of genomic variants to drive drug repurposing for chronic hepatitis B

Background

One of the main challenges in personalized medicine is to establish and apply a large number of variants from genomic databases into clinical diagnostics and further facilitate genome-driven drug repurposing. By utilizing biological chronic hepatitis B infection (CHB) risk genes, our study proposed a systematic approach to use genomic variants to drive drug repurposing for CHB.

Method

The genomic variants were retrieved from the Genome-Wide Association Study (GWAS) and Phenome-Wide Association Study (PheWAS) databases. Then, the biological CHB risk genes crucial for CHB progression were prioritized based on the scoring system devised with five strict functional annotation criteria. A score of ≥ 2 were categorized as the biological CHB risk genes and further shed light on drug target genes for CHB treatments. Overlapping druggable targets were identified using two drug databases (DrugBank and Drug-Gene Interaction Database (DGIdb)).

Results

A total of 44 biological CHB risk genes were screened based on the scoring system from five functional annotation criteria. Interestingly, we found 6 druggable targets that overlapped with 18 drugs with status of undergoing clinical trials for CHB, and 9 druggable targets that overlapped with 20 drugs undergoing preclinical investigations for CHB. Eight druggable targets were identified, overlapping with 25 drugs that can potentially be repurposed for CHB. Notably, CD40 and HLA-DPB1 were identified as promising targets for CHB drug repurposing based on the target scores.

Conclusion

Through the integration of genomic variants and a bioinformatic approach, our findings suggested the plausibility of CHB genomic variant-driven drug repurposing for CHB.

Identification of Druggable Genes for Asthma by Integrated Genomic Network Analysis

Asthma is a common and heterogeneous disease characterized by chronic airway inflammation. Currently, the two main types of asthma medicines are inhaled corticosteroids and long-acting β2-adrenoceptor agonists (LABAs). In addition, biological drugs provide another therapeutic option, especially for patients with severe asthma. However, these drugs were less effective in preventing severe asthma exacerbation, and other drug options are still limited. Herein, we extracted asthma-associated single nucleotide polymorphisms (SNPs) from the genome-wide association studies (GWAS) and phenome-wide association studies (PheWAS) catalog and prioritized candidate genes through five functional annotations. Genes enriched in more than two categories were defined as “biological asthma risk genes.” Then, DrugBank was used to match target genes with FDA-approved medications and identify candidate drugs for asthma. We discovered 139 biological asthma risk genes and identified 64 drugs targeting 22 of these genes. Seven of them were approved for asthma, including reslizumab, mepolizumab, theophylline, dyphylline, aminophylline, oxtriphylline, and enprofylline. We also found 17 drugs with clinical or preclinical evidence in treating asthma. In addition, eleven of the 40 candidate drugs were further identified as promising asthma therapy. Noteworthy, IL6R is considered a target for asthma drug repurposing based on its high target scores. Through in silico drug repurposing approach, we identified sarilumab and satralizumab as the most promising drug for asthma treatment.

Drug Repurposing for Atopic Dermatitis by Integration of Gene Networking and Genomic Information

Atopic Dermatitis (AD) is a chronic and relapsing skin disease. The medications for treating AD are still limited, most of them are topical corticosteroid creams or antibiotics. The current study attempted to discover potential AD treatments by integrating a gene network and genomic analytic approaches. Herein, the Single Nucleotide Polymorphism (SNPs) associated with AD were extracted from the GWAS catalog. We identified 70 AD-associated loci, and then 94 AD risk genes were found by extending to proximal SNPs based on r2 > 0.8 in Asian populations using HaploReg v4.1. Next, we prioritized the AD risk genes using in silico pipelines of bioinformatic analysis based on six functional annotations to identify biological AD risk genes. Finally, we expanded them according to the molecular interactions using the STRING database to find the drug target genes. Our analysis showed 27 biological AD risk genes, and they were mapped to 76 drug target genes. According to DrugBank and Therapeutic Target Database, 25 drug target genes overlapping with 53 drugs were identified. Importantly, dupilumab, which is approved for AD, was successfully identified in this bioinformatic analysis. Furthermore, ten drugs were found to be potentially useful for AD with clinical or preclinical evidence. In particular, we identified filgotinub and fedratinib, targeting gene JAK1, as potential drugs for AD. Furthermore, four monoclonal antibody drugs (lebrikizumab, tralokinumab, tocilizumab, and canakinumab) were successfully identified as promising for AD repurposing. In sum, the results showed the feasibility of gene networking and genomic information as a potential drug discovery resource.

Genetic variants that influence SARS-CoV-2 receptor TMPRSS2 expression among population cohorts from multiple continents

The World Health Organization recently announced that pandemic status has been achieved for coronavirus disease 2019 (COVID-19), which is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Exponential increases in patient numbers have been reported around the world, along with proportional increases in the number of COVID-19-related deaths. The SARS-CoV-2 infection rate in a population is expected to be influenced by social practices, availability of vaccines or prophylactics, and the prevalence of susceptibility genes in the population. Previous work revealed that cellular uptake of SARS-CoV-2 requires Angiotensin Converting Enzyme 2 (ACE-2) and a cellular protease. The spike (S) protein on SARS-CoV-2 binds ACE-2, which functions as an entry receptor. Following receptor binding, transmembrane protease serine 2 (encoded by TMPRSS2) primes the S protein to allow cellular uptake. Therefore, individual expression of TMPRSS2 may be a crucial determinant of SARS-CoV-2 infection susceptibility. Here, we utilized multiple large genome databases, including the GTEx portal, SNP nexus, and Ensembl genome project, to identify gene expression profiles for TMPRSS2 and its important expression quantitative trait loci. Our results show that four variants (rs464397, rs469390, rs2070788 and rs383510) affect expression of TMPRSS2 in lung tissue. The allele frequency of each variant was then assessed in regional populations, including African, American, European, and three Asian cohorts (China, Japan and Taiwan). Interestingly, our data shows that TMPRSS2-upregulating variants are at higher frequencies in European and American populations than in the Asian populations, which implies that these populations might be relatively susceptible to SARS-CoV-2 infection.

The Effect of Education through Short Message Service (SMS) Messages on Diabetic Patients Adherence

Poor adherence and a lack of understanding of medication instructions for oral antidiabetic use are key factors that inhibit the control of glycemic levels. The aforementioned situation needs intervention to improve medication adherence and the therapy. This study was conducted with a quasi-experimental design with prospective data collection. The subjects of this study were 50 outpatients with type 2 diabetes melitus (T2DM) who had received oral antidiabetic medicine therapy at least six months prior to adherence measurement. The patients were classified into two groups-the control group and the intervention group. The intervention group received Short Message Service (SMS) messages of diabetes education, while the control group did not. Data collection was conducted by doing interviews and administering the Morisky Medication Adherence Scale (MMAS) questionnaire. The results showed the increase in adherence in the intervention group as 1.15 ± 1.04 and that in the control group as 0.72 ± 0.90. These results indicated that there were significant differences in MMAS score between the control and intervention groups (p < 0.05). The decrease in fasting blood glucose and glucose measured 2 h postprandially was greater in the intervention group than that in the control group. It was concluded that the provision of education through SMS had a positive effect on medication adherence and glycemic levels.