Nasopharyngeal airway long noncoding RNAs of infants with bronchiolitis and subsequent risk of developing childhood asthma.

BACKGROUND: Severe bronchiolitis (ie, bronchiolitis requiring hospitalization) during infancy is a major risk factor for developing childhood asthma. However, the biological mechanisms linking these 2 conditions remain unclear. OBJECTIVE: We sought to investigate the longitudinal relationship between nasopharyngeal airway long noncoding RNA (lncRNA) in infants with severe bronchiolitis and subsequent asthma development. METHODS: In this multicenter prospective cohort study of infants with severe bronchiolitis, we performed RNA sequencing of nasopharyngeal airway lncRNAs at index hospitalization. First, we identified differentially expressed lncRNAs (DE-lncRNAs) associated with asthma development by age 6 years. Second, we investigated the associations of DE-lncRNAs with asthma-related clinical characteristics. Third, to characterize the function of DE-lncRNAs, we performed pathway analysis for mRNA targeted by DE-lncRNAs. Finally, we examined the associations of DE-lncRNAs with nasal cytokines at index hospitalization. RESULTS: Among 343 infants with severe bronchiolitis (median age, 3 months), we identified 190 DE-lncRNAs (false-discovery rate [FDR] < 0.05) associated with asthma development (eg, LINC02145, RAMP2-AS1, and PVT1). These DE-lncRNAs were associated with asthma-related clinical characteristics (FDR < 0.05), for example, respiratory syncytial virus or rhinovirus infection, infant eczema, and IgE sensitization. Furthermore, DE-lncRNAs were characterized by asthma-related pathways, including mitogen-activated protein kinase, FcɛR, and phosphatidylinositol 3-kinase (PI3K)-protein kinase B signaling pathways (FDR < 0.05). These DE-lncRNAs were also associated with nasal cytokines (eg, IL-1ß, IL-4, and IL-13; FDR < 0.05). CONCLUSIONS: In a multicenter cohort study of infants with severe bronchiolitis, we identified nasopharyngeal airway lncRNAs associated with childhood asthma development, characterized by asthma-related clinical characteristics, asthma-related pathways, and nasal cytokines. Our approach identifies lncRNAs underlying the bronchiolitis-asthma link and facilitates the early identification of infants at high risk of subsequent asthma development.

Integrated nasopharyngeal airway metagenome and asthma genetic risk endotyping of severe bronchiolitis in infancy and risk of childhood asthma.

The severe bronchiolitis endotype characterized by a high abundance of H. influenzae, high proportion of RV-A and RV-C infections, and high asthma genetic risk had a significantly higher risk for developing asthma. BACKGROUND: Infants with bronchiolitis are at increased risk for developing asthma. Growing evidence suggests bronchiolitis is a heterogeneous condition. However, little is known about its biologically distinct subgroups based on the integrated metagenome and asthma genetic risk signature and their longitudinal relationships with asthma development. METHODS: In a multi-center prospective cohort study of infants with severe bronchiolitis (i.e., bronchiolitis requiring hospitalization), we profiled nasopharyngeal airway metagenome and virus at hospitalization, and calculated the polygenic risk score of asthma. Using similarity network fusion clustering approach, we identified integrated metagenome-asthma genetic risk endotypes. We also examined their longitudinal association with the risk of developing asthma by age six years. RESULTS: Of 450 infants with bronchiolitis (median age, 3 months), we identified five distinct endotypes-characterized by their nasopharyngeal metagenome, virus, and asthma genetic risk profiles. Compared with endotype A infants (who clinically resembled "classic" bronchiolitis), endotype E infants (characterized by a high abundance of H. influenzae, high proportion of RV-A and RV-C infections, and high asthma genetic risk) had a significantly higher risk for developing asthma (35.9% versus 16.7%; ORadj, 2.24; 95%CI, 1.02-4.97; p=0.046). The pathway analysis showed that endotype E had enriched microbial pathways (e.g., glycolysis, L-lysine, arginine metabolism) and host pathways (e.g., IFNs, IL-6/JAK/STAT3, fatty acids, MHC, and immunoglobin-related) (FDR<0.05). Additionally, endotype E had a significantly higher proportion of neutrophils (FDR<0.05). CONCLUSION: In this multi-center prospective cohort study of infant bronchiolitis, the clustering analysis of integrated-omics data identified biologically distinct endotypes with differential risks for developing asthma.

Race-Neutral Equations and Pulmonary Function Test Interpretation in Two Pediatric Cohorts.

OBJECTIVE: To investigate the changes in predicted lung function measurements when using race-neutral equations in children, based upon the new Global Lung Initiative (GLI) reference equations, utilizing a race-neutral approach in interpreting spirometry results compared with the 2012 race-specific GLI equations. STUDY DESIGN: We analyzed data from 2 multicenter prospective cohorts comprised of healthy children and children with history of severe (requiring hospitalization) bronchiolitis. Spirometry testing was done at the 6-year physical exam, and 677 tests were analyzed using new GLI Global and 2012 GLI equations. We used multivariable logistic regression, adjusted for age, height, and sex, to examine the association of race with the development of new impairment or increased severity (forced expiratory volume in the first second (FEV1) z-score ≤ -1.645) as per 2022 American Thoracic Society (ATS) guidelines. RESULTS: Compared with the race-specific GLI, the race-neutral equation yielded increases in the median forced expiratory volume in the first second and forced vital capacity (FVC) percent predicted in White children but decreases in these two measures in Black children. The prevalence of obstruction increased in White children by 21%, and the prevalence of possible restriction increased in Black children by 222%. Compared with White race, Black race was associated with increased prevalence of new impairments (aOR 7.59; 95%CI, 3.00-19.67; P < .001) and increased severity (aOR 35.40; 95%CI, 4.70-266.40; P = .001). Results were similar across both cohorts. CONCLUSIONS: As there are no biological justifications for the inclusion of race in spirometry interpretation, use of race-neutral spirometry reference equations led to an increase in both the prevalence and severity of respiratory impairments among Black children.

Early-life exposure to residential greenness and risk of asthma in a U.S. bronchiolitis cohort.

INTRODUCTION: Severe bronchiolitis (i.e., bronchiolitis requiring hospitalization) is linked to childhood asthma development. Despite a growing understanding of risk factors for developing post-bronchiolitis asthma, protective factors remain unclear. In this study, we aimed to investigate whether exposure to residential greenness between birth and bronchiolitis hospitalization is associated with asthma and atopic asthma development by age 6 years. METHODS: We analyzed a US severe bronchiolitis cohort from hospitalization to age 6 years, investigating how the normalized difference vegetation index (NDVI) and chlorophyll index green (CI green), measured in small (100 m) and large (500 m) radiuses around homes, relate to asthma and atopic asthma by age 6 years. We also explored whether maternal antibiotic use, daycare attendance, and respiratory virus type during hospitalization act as effect modifiers. RESULTS: The study cohort included 861 infants, with 239 (28%) developing asthma by age 6 years-152 atopic, 17 nonatopic, and 70 unclassified. Early life residential exposure to high NDVI and CI green levels was associated with lower odds of asthma (ORAdj for NDVI within a 100 m radius, 0.18; 95% CI, 0.05-0.78; and ORAdj for CI green levels within a 100 m radius, 0.53; 95% CI, 0.31-0.90). Associations also were significant for the development of atopic asthma (ORAdj 0.16; 95% CI, 0.03-0.96; and ORAdj 0.46; 95% CI, 0.25-0.92; respectively). Results were similar for the 500 m radius exposures. No effect modification was noted. CONCLUSION: In a U.S. bronchiolitis cohort, exposure to residential greenness between birth and bronchiolitis hospitalization is linked to lower asthma and atopic asthma risk by age 6 years.

The bZIP transcription factor BATF3/ZIP-10 suppresses innate immunity by attenuating PMK-1/p38 signaling.

Innate immunity is the first line of host defense against pathogenic invasion in metazoans. The transcription factor basic leucine zipper transcriptional factor ATF-like 3 (BATF3) plays a crucial role in the development of conventional dendritic cells and the program of CD8 + T cell survival and memory, but the role of BATF3 in innate immune responses remains unclear. Here, we show an evolutionarily conserved basic-region leucine zipper (bZIP) transcription factor BATF3/ZIP-10 suppresses innate immune response through repressing the p38/PMK-1 mitogen-activated protein kinase (MAPK) pathway in vitro and in vivo. The worm mutant lacking the Caenorhabditis elegans homolog BATF3, ZIP-10, exhibited enhanced resistance to PA14 infection, which was completely rescued by transgenic expression of either endogenous zip-10 or mouse or human Batf3 cDNA driven by the worm zip-10 promoter. ZIP-10 expression was inhibited by a microRNA miR-60 that was downregulated upon PA14 infection. Moreover, the level of phosphorylated but not total PMK-1/p38 was attenuated by ZIP-10 and stimulated by miR-60. The human HEK293 cells with Batf3 overexpression or RNA-interference knockdown exhibited a reduction or increase of the cell viability upon Pseudomonas aeruginosa PA14 infection, respectively. The overexpression of either worm ZIP-10 or human BATF3 abolished the activation of p38 and inhibited the expression of antimicrobial peptides and cytokine genes in HEK293 cells. Our findings indicate that the genetic transcriptional program of the evolutionally conserved bZIP transcription factor BATF3/ZIP-10 suppresses innate immunity by attenuating the p38 MAPK signaling activity, which expands our understanding of the pathological mechanisms underlying relevant infectious diseases.

Respiratory Virus-Specific Nasopharyngeal Lipidome Signatures and Severity in Infants With Bronchiolitis: A Prospective Multicenter Study.

BACKGROUND: In infant bronchiolitis, recent evidence indicates that respiratory viruses (eg, respiratory syncytial virus [RSV], rhinovirus [RV]) contribute to the heterogeneity of disease severity. Of the potential pathobiological molecules, lipids serve as signaling molecules in airway inflammation. However, little is known about the role of the airway lipidome in between-virus heterogeneity and disease severity. METHODS: In this multicenter prospective study of 800 infants hospitalized for RSV or RV bronchiolitis, we analyzed nasopharyngeal lipidome data. We examined discriminatory lipids between RSV and RV infection and the association of the discriminatory lipids with bronchiolitis severity, defined by positive pressure ventilation (PPV) use. RESULTS: We identified 30 discriminatory nasopharyngeal lipid species and 8 fatty acids between RSV and RV infection. In the multivariable models adjusting for patient-level confounders, 8 lipid species-for example, phosphatidylcholine (18:2/18:2) (adjusted odds ratio [aOR], 0.23 [95% confidence interval {CI}, .11-.44]; false discovery rate [FDR] = 0.0004) and dihydroceramide (16:0) (aOR, 2.17 [95% CI, 1.12-3.96]; FDR = 0.04)-were significantly associated with the risk of PPV use. Additionally, 6 fatty acids-for example, eicosapentaenoic acid (aOR, 0.27 [95% CI, .11-.57]; FDR = 0.01)-were also significantly associated with the risk of PPV use. CONCLUSIONS: In infants hospitalized for bronchiolitis, the nasopharyngeal lipidome plays an important role in the pathophysiology of between-virus heterogeneity and disease severity.

Nasal airway microRNA profiling of infants with severe bronchiolitis and risk of childhood asthma: a multicentre prospective study.

BACKGROUND: Severe bronchiolitis (i.e. bronchiolitis requiring hospitalisation) during infancy is a major risk factor for childhood asthma. However, the exact mechanism linking these common conditions remains unclear. We examined the longitudinal relationship between nasal airway miRNAs during severe bronchiolitis and the risk of developing asthma. METHODS: In a 17-centre prospective cohort study of infants with severe bronchiolitis, we sequenced their nasal microRNA at hospitalisation. First, we identified differentially expressed microRNAs (DEmiRNAs) associated with the risk of developing asthma by age 6 years. Second, we characterised the DEmiRNAs based on their association with asthma-related clinical features, and expression level by tissue and cell types. Third, we conducted pathway and network analyses by integrating DEmiRNAs and their mRNA targets. Finally, we investigated the association of DEmiRNAs and nasal cytokines. RESULTS: In 575 infants (median age 3 months), we identified 23 DEmiRNAs associated with asthma development (e.g. hsa-miR-29a-3p; false discovery rate (FDR) <0.10), particularly in infants with respiratory syncytial virus infection (FDR for the interaction <0.05). These DEmiRNAs were associated with 16 asthma-related clinical features (FDR <0.05), e.g. infant eczema and corticosteroid use during hospitalisation. In addition, these DEmiRNAs were highly expressed in lung tissue and immune cells (e.g. T-helper cells, neutrophils). Third, DEmiRNAs were negatively correlated with their mRNA targets (e.g. hsa-miR-324-3p/IL13), which were enriched in asthma-related pathways (FDR <0.05), e.g. toll-like receptor, PI3K-Akt and FcɛR signalling pathways, and validated by cytokine data. CONCLUSION: In a multicentre cohort of infants with severe bronchiolitis, we identified nasal miRNAs during illness that were associated with major asthma-related clinical features, immune response, and risk of asthma development.

VirusCircBase: a database of virus circular RNAs.

Circular RNAs (circRNAs) are covalently closed long noncoding RNAs critical in diverse cellular activities and multiple human diseases. Several cancer-related viral circRNAs have been identified in double-stranded DNA viruses (dsDNA), yet no systematic study about the viral circRNAs has been reported. Herein, we have performed a systematic survey of 11 924 circRNAs from 23 viral species by computational prediction of viral circRNAs from viral-infection-related RNA sequencing data. Besides the dsDNA viruses, our study has also revealed lots of circRNAs in single-stranded RNA viruses and retro-transcribing viruses, such as the Zika virus, the Influenza A virus, the Zaire ebolavirus, and the Human immunodeficiency virus 1. Most viral circRNAs had reverse complementary sequences or repeated sequences at the flanking sequences of the back-splice sites. Most viral circRNAs only expressed in a specific cell line or tissue in a specific species. Functional enrichment analysis indicated that the viral circRNAs from dsDNA viruses were involved in KEGG pathways associated with cancer. All viral circRNAs presented in the current study were stored and organized in VirusCircBase, which is freely available at http://www.computationalbiology.cn/ViruscircBase/home.html and is the first virus circRNA database. VirusCircBase forms the fundamental atlas for the further exploration and investigation of viral circRNAs in the context of public health.

Classification and characterization of multigene family proteins of African swine fever viruses.

African swine fever virus (ASFV) poses serious threats to the pig industry. The multigene family (MGF) proteins are extensively distributed in ASFVs and are generally classified into five families, including MGF-100, MGF-110, MGF-300, MGF-360 and MGF-505. Most MGF proteins, however, have not been well characterized and classified within each family. To bridge this gap, this study first classified MGF proteins into 31 groups based on protein sequence homology and network clustering. A web server for classifying MGF proteins was established and kept available for free at http://www.computationalbiology.cn/MGF/home.html. Results showed that MGF groups of the same family were most similar to each other and had conserved sequence motifs; the genetic diversity of MGF groups varied widely, mainly due to the occurrence of indels. In addition, the MGF proteins were predicted to have large structural and functional diversity, and MGF proteins of the same MGF family tended to have similar structure, location and function. Reconstruction of the ancestral states of MGF groups along the ASFV phylogeny showed that most MGF groups experienced either the copy number variations or the gain-or-loss changes, and most of these changes happened within strains of the same genotype. It is found that the copy number decrease and the loss of MGF groups were much larger than the copy number increase and the gain of MGF groups, respectively, suggesting the ASFV tended to lose MGF proteins in the evolution. Overall, the work provides a detailed classification for MGF proteins and would facilitate further research on MGF proteins.

Identification and characterization of circRNAs encoded by MERS-CoV, SARS-CoV-1 and SARS-CoV-2.

The life-threatening coronaviruses MERS-CoV, SARS-CoV-1 and SARS-CoV-2 (SARS-CoV-1/2) have caused and will continue to cause enormous morbidity and mortality to humans. Virus-encoded noncoding RNAs are poorly understood in coronaviruses. Data mining of viral-infection-related RNA-sequencing data has resulted in the identification of 28 754, 720 and 3437 circRNAs encoded by MERS-CoV, SARS-CoV-1 and SARS-CoV-2, respectively. MERS-CoV exhibits much more prominent ability to encode circRNAs in all genomic regions than those of SARS-CoV-1/2. Viral circRNAs typically exhibit low expression levels. Moreover, majority of the viral circRNAs exhibit expressions only in the late stage of viral infection. Analysis of the competitive interactions of viral circRNAs, human miRNAs and mRNAs in MERS-CoV infections reveals that viral circRNAs up-regulated genes related to mRNA splicing and processing in the early stage of viral infection, and regulated genes involved in diverse functions including cancer, metabolism, autophagy, viral infection in the late stage of viral infection. Similar analysis in SARS-CoV-2 infections reveals that its viral circRNAs down-regulated genes associated with metabolic processes of cholesterol, alcohol, fatty acid and up-regulated genes associated with cellular responses to oxidative stress in the late stage of viral infection. A few genes regulated by viral circRNAs from both MERS-CoV and SARS-CoV-2 were enriched in several biological processes such as response to reactive oxygen and centrosome localization. This study provides the first glimpse into viral circRNAs in three deadly coronaviruses and would serve as a valuable resource for further studies of circRNAs in coronaviruses.

Severe bronchiolitis profiling as the first step towards prevention of asthma.

Bronchiolitis is the most common respiratory infection leading to hospitalization and constitutes a significant healthcare burden. The two main viral agents causing bronchiolitis, respiratory syncytial virus (RSV) and rhinovirus (RV), have distinct cytopathic, immune response, and clinical characteristics. Different approaches have been suggested for subtyping bronchiolitis based on viral etiology, atopic status, transcriptome profiles in blood, airway metabolome, lipidomic data, and airway microbiota. The highest risk of asthma at school age has been in a subgroup of bronchiolitis characterized by older age, high prevalence of RV infection, previous breathing problems, and/or eczema. Regarding solely viral etiology, RV-bronchiolitis in infancy has been linked to a nearly three times higher risk of developing asthma than RSV-bronchiolitis. Although treatment with betamimetics and systemic corticosteroids has been found ineffective in bronchiolitis overall, it can be beneficial for infants with severe RV bronchiolitis. Thus, there is a need to develop a more individualized therapeutic approach for bronchiolitis and follow-up strategies for infants at higher risk of asthma in the future perspective.

Metabolome subtyping of severe bronchiolitis in infancy and risk of childhood asthma.

BACKGROUND: Infants with bronchiolitis are at increased risk for developing asthma. Growing evidence suggests bronchiolitis is a heterogeneous condition. OBJECTIVES: We sought to identify biologically distinct subgroups based on the metabolome signatures (metabotypes) in infants with severe bronchiolitis and to examine the longitudinal relationships of metabotypes with asthma development. METHODS: In a multicenter prospective cohort study of infants (age, <12 months) hospitalized for bronchiolitis, the nasopharyngeal airway metabolome was profiled at hospitalization. Using a clustering approach, this study identified mutually exclusive metabotypes. This study also examined their longitudinal association with the risk of developing asthma by 5 years of age. RESULTS: Of 918 infants hospitalized for bronchiolitis (median age, 3 months), this study identified 5 distinct metabotypes-characterized by their nasopharyngeal metabolome profile: A, glycerophosphocholine-high; B, amino acid-high, polyunsaturated fatty acid-low; C, amino acid-high, glycerophospholipid-low; D, glycerophospholipid-high; and E, mixed. Compared with infants with metabotype A (who clinically resembled "classic" bronchiolitis), infants with metabotype B had a significantly higher risk for developing asthma (23% vs 41%; adjusted odds ratio, 2.22; 95% CI, 1.07-4.69). The pathway analysis showed that metabotype B had enriched amino acid (eg, methionine, histidine, glutathione) and α-linolenic/linoleic acid metabolism pathways (false discovery rate, <5 × 10-14 for all). Finally, the transcriptome analysis revealed that infants with metabotype B had upregulated IFN-α and IL-6/JAK/STAT3 pathways and downregulated fatty acid metabolism pathways (false discovery rate, <0.05 for both). CONCLUSIONS: In this multicenter prospective cohort study of infants with severe bronchiolitis, the clustering analysis of metabolome data identified biologically distinct metabotypes, including a metabotype characterized by high inflammatory amino acids and low polyunsaturated fatty acids that is at significantly increased risk for developing asthma.

Nasopharyngeal airway dual-transcriptome of infants with severe bronchiolitis and risk of childhood asthma: A multicenter prospective study.

BACKGROUND: Severe bronchiolitis (ie, bronchiolitis requiring hospitalization) during infancy is a major risk factor for childhood asthma. However, the exact mechanism linking these common conditions remains unclear. OBJECTIVES: This study sought to examine the integrated role of airway microbiome (both taxonomy and function) and host response in asthma development in this high-risk population. METHODS: This multicenter prospective cohort study of 244 infants with severe bronchiolitis (median age, 3 months) examined the infants' nasopharyngeal metatranscriptomes (microbiomes) and transcriptomes (hosts), as well as metabolomes at hospitalization. The longitudinal relationships investigated include (1) major bacterial species (Streptococcus pneumoniae, Haemophilus influenzae, and Moraxella catarrhalis), (2) microbial function, and (3) host response with risks of developing asthma by age 6 years. RESULTS: First, the abundance of S pneumoniae was associated with greater risks of asthma (P = .01), particularly in infants with nonrhinovirus infection (Pinteraction = .04). Second, of 328 microbial functional pathways that are differentially enriched by asthma development, the top pathways (eg, fatty acid and glycolysis pathways; false discovery rate [FDR] < 1 × 10-12) were driven by these 3 major species (eg, positive association of S pneumoniae with glycolysis; FDR < 0.001). These microbial functional pathways were validated with the parallel metabolome data. Third, 104 transcriptome pathways were differentially enriched (FDR < .05)-for example, downregulated interferon-α and -γ and upregulated T-cell activation pathways. S pneumoniae was associated with most differentially expressed transcripts (eg, DAGLB; FDR < 0.05). CONCLUSIONS: By applying metatranscriptomic, transcriptomic, and metabolomic approaches to a multicenter cohort of infants with bronchiolitis, this study found an interplay between major bacterial species, their function, and host response in the airway, and their longitudinal relationship with asthma development.

A genome-wide cross-trait analysis identifies shared loci and causal relationships of type 2 diabetes and glycaemic traits with polycystic ovary syndrome.

AIMS/HYPOTHESIS: The link underlying abnormal glucose metabolism, type 2 diabetes and polycystic ovary syndrome (PCOS) that is independent of BMI remains unclear in observational studies. We aimed to clarify this association using a genome-wide cross-trait approach. METHODS: Summary statistics from the hitherto largest genome-wide association studies conducted for type 2 diabetes, type 2 diabetes mellitus adjusted for BMI (T2DMadjBMI), fasting glucose, fasting insulin, 2h glucose after an oral glucose challenge (all adjusted for BMI), HbA1c and PCOS, all in populations of European ancestry, were used. We quantified overall and local genetic correlations, identified pleiotropic loci and expression-trait associations, and made causal inferences across traits. RESULTS: A positive overall genetic correlation between type 2 diabetes and PCOS was observed, largely influenced by BMI (rg=0.31, p=1.63×10-8) but also independent of BMI (T2DMadjBMI-PCOS: rg=0.12, p=0.03). Sixteen pleiotropic loci affecting type 2 diabetes, glycaemic traits and PCOS were identified, suggesting mechanisms of association that are independent of BMI. Two shared expression-trait associations were found for type 2 diabetes/T2DMadjBMI and PCOS targeting tissues of the cardiovascular, exocrine/endocrine and digestive systems. A putative causal effect of fasting insulin adjusted for BMI and type 2 diabetes on PCOS was demonstrated. CONCLUSIONS/INTERPRETATION: We found a genetic link underlying type 2 diabetes, glycaemic traits and PCOS, driven by both biological pleiotropy and causal mediation, some of which is independent of BMI. Our findings highlight the importance of controlling fasting insulin levels to mitigate the risk of PCOS, as well as screening for and long-term monitoring of type 2 diabetes in all women with PCOS, irrespective of BMI.

Nasopharyngeal lipidomic endotypes of infants with bronchiolitis and risk of childhood asthma: a multicentre prospective study.

BACKGROUND: Bronchiolitis is the leading cause of hospitalisation of US infants and an important risk factor for childhood asthma. Recent evidence suggests that bronchiolitis is clinically heterogeneous. We sought to derive bronchiolitis endotypes by integrating clinical, virus and lipidomics data and to examine their relationship with subsequent asthma risk. METHODS: This is a multicentre prospective cohort study of infants (age <12 months) hospitalised for bronchiolitis. We identified endotypes by applying clustering approaches to clinical, virus and nasopharyngeal airway lipidomic data measured at hospitalisation. We then determined their longitudinal association with the risk for developing asthma by age 6 years by fitting a mixed-effects logistic regression model. To account for multiple comparisons of the lipidomics data, we computed the false discovery rate (FDR). To understand the underlying biological mechanism of the endotypes, we also applied pathway analyses to the lipidomics data. RESULTS: Of 917 infants with bronchiolitis (median age, 3 months), we identified clinically and biologically meaningful lipidomic endotypes: (A) cinicalclassiclipidmixed (n=263), (B) clinicalseverelipidsphingolipids-high (n=281), (C) clinicalmoderatelipidphospholipids-high (n=212) and (D) clinicalatopiclipidsphingolipids-low (n=161). Endotype A infants were characterised by 'classic' clinical presentation of bronchiolitis. Profile D infants were characterised by a higher proportion of parental asthma, IgE sensitisation and rhinovirus infection and low sphingolipids (eg, sphingomyelins, ceramides). Compared with endotype A, profile D infants had a significantly higher risk of asthma (22% vs 50%; unadjusted OR, 3.60; 95% CI 2.31 to 5.62; p<0.001). Additionally, endotype D had a significantly lower abundance of polyunsaturated fatty acids (eg, docosahexaenoic acid; FDR=0.01). The pathway analysis revealed that sphingolipid metabolism pathway was differentially expressed in endotype D (FDR=0.048). CONCLUSIONS: In this multicentre prospective cohort study of infants with bronchiolitis, integrated clustering of clinical, virus and lipidomic data identified clinically and biologically distinct endotypes that have a significantly differential risk for developing asthma.Delete.

Genomic correlation, shared loci, and causal relationship between obesity and polycystic ovary syndrome: a large-scale genome-wide cross-trait analysis.

BACKGROUND: The comorbidity between polycystic ovary syndrome (PCOS) and obesity has long been observed in clinical settings, but their shared genetic basis remains unclear. METHODS: Leveraging summary statistics of large-scale GWAS(s) conducted in European-ancestry populations on body mass index (adult BMI, Nfemale=434,794; childhood BMI, N=39,620), waist-to-hip ratio (WHR, Nfemale=381,152), WHR adjusted for BMI (WHRadjBMI, Nfemale=379,501), and PCOS (Ncase=10,074, Ncontrol=103,164), we performed a large-scale genome-wide cross-trait analysis to quantify overall and local genetic correlation, to identify shared loci, and to infer causal relationship. RESULTS: We found positive genetic correlations between PCOS and adult BMI (rg=0.47, P=2.19×10-16), childhood BMI (rg=0.31, P=6.72×10-5), and WHR (rg=0.32, P=1.34×10-10), all withstanding Bonferroni correction. A suggestive significant genetic correlation was found between PCOS and WHRadjBMI (rg=0.09, P=0.04). Partitioning the whole genome into 1703 nearly independent regions, we observed a significant local genetic correlation for adult BMI and PCOS at chromosome 18: 57630483-59020751. We identified 16 shared loci underlying PCOS and obesity-related traits via cross-trait meta-analysis including 9 loci shared between BMI and PCOS (adult BMI and PCOS: 5 loci; childhood BMI and PCOS: 4 loci), 6 loci shared between WHR and PCOS, and 5 loci shared between WHRadjBMI and PCOS. Mendelian randomization (MR) supported the causal roles of both adult BMI (OR=2.92, 95% CI=2.33-3.67) and childhood BMI (OR=2.76, 95% CI=2.09-3.66) in PCOS, but not WHR (OR=1.19, 95% CI=0.93-1.52) or WHRadjBMI (OR=1.03, 95% CI=0.87-1.22). Genetic predisposition to PCOS did not seem to influence the risk of obesity-related traits. CONCLUSIONS: Our cross-trait analysis suggests a shared genetic basis underlying obesity and PCOS and provides novel insights into the biological mechanisms underlying these complex traits. Our work informs public health intervention by confirming the important role of weight management in PCOS prevention.

Nasopharyngeal metatranscriptome profiles of infants with bronchiolitis and risk of childhood asthma: a multicentre prospective study.

BACKGROUND: Bronchiolitis is not only the leading cause of hospitalisation in US infants but also a major risk factor for asthma development. Growing evidence supports clinical heterogeneity within bronchiolitis. Our objectives were to identify metatranscriptome profiles of infant bronchiolitis, and to examine their relationship with the host transcriptome and subsequent asthma development. METHODS: As part of a multicentre prospective cohort study of infants (age <1 year) hospitalised for bronchiolitis, we integrated virus and nasopharyngeal metatranscriptome (species-level taxonomy and function) data measured at hospitalisation. We applied network-based clustering approaches to identify metatranscriptome profiles. We then examined their association with the host transcriptome at hospitalisation and risk for developing asthma. RESULTS: We identified five metatranscriptome profiles of bronchiolitis (n=244): profile A: virusRSVmicrobiomecommensals; profile B: virusRSV/RV-Amicrobiome H.influenzae ; profile C: virusRSVmicrobiome S.pneumoniae ; profile D: virusRSVmicrobiome M.nonliquefaciens ; and profile E: virusRSV/RV-Cmicrobiome M.catarrhalis . Compared with profile A, profile B infants were characterised by a high proportion of eczema, Haemophilus influenzae abundance and enriched virulence related to antibiotic resistance. These profile B infants also had upregulated T-helper 17 and downregulated type I interferon pathways (false discovery rate (FDR) <0.005), and significantly higher risk for developing asthma (17.9% versus 38.9%; adjusted OR 2.81, 95% CI 1.11-7.26). Likewise, profile C infants were characterised by a high proportion of parental asthma, Streptococcus pneumoniae dominance, and enriched glycerolipid and glycerophospholipid metabolism of the microbiome. These profile C infants had an upregulated RAGE signalling pathway (FDR <0.005) and higher risk of asthma (17.9% versus 35.6%; adjusted OR 2.49, 95% CI 1.10-5.87). CONCLUSIONS: Metatranscriptome and clustering analysis identified biologically distinct metatranscriptome profiles that have differential risks of asthma.

Proteomics endotyping of infants with severe bronchiolitis and risk of childhood asthma.

BACKGROUND: Bronchiolitis is the leading cause of hospitalization in U.S. infants and a major risk factor for childhood asthma. Growing evidence supports clinical heterogeneity within bronchiolitis. We aimed to identify endotypes of infant bronchiolitis by integrating clinical, virus, and serum proteome data, and examine their relationships with asthma development. METHODS: This is a multicenter prospective cohort study of infants hospitalized for physician-diagnosis of bronchiolitis. We identified bronchiolitis endotypes by applying unsupervised machine learning (clustering) approaches to integrated clinical, virus (respiratory syncytial virus [RSV], rhinovirus [RV]), and serum proteome data measured at hospitalization. We then examined their longitudinal association with the risk for developing asthma by age 6 years. RESULTS: In 140 infants hospitalized with bronchiolitis, we identified three endotypes: (1) clinicalatopic virusRV proteomeNFκB-dysregulated , (2) clinicalnon-atopic virusRSV/RV proteomeTNF-dysregulated , and (3) clinicalclassic virusRSV proteomeNFκB/TNF-regulated endotypes. Endotype 1 infants were characterized by high proportion of IgE sensitization and RV infection. These endotype 1 infants also had dysregulated NFκB pathways (FDR < 0.001) and significantly higher risks for developing asthma (53% vs. 22%; adjOR 4.04; 95% CI, 1.49-11.0; p = 0.006), compared with endotype 3 (clinically resembling "classic" bronchiolitis). Likewise, endotype 2 infants were characterized by low proportion of IgE sensitization and high proportion of RSV or RV infection. These endotype 2 infants had dysregulated tumor necrosis factor (TNF)-mediated signaling pathway (FDR <0.001) and significantly higher risks for developing asthma (44% vs. 22%; adjOR 2.71; 95% CI, 1.03-7.11, p = 0.04). CONCLUSION: In this multicenter cohort, integrated clustering of clinical, virus, and proteome data identified biologically distinct endotypes of bronchiolitis that have differential risks of asthma development.

Systematic identification and characterization of repeat sequences in African swine fever virus genomes.

African swine fever virus (ASFV) is a large DNA virus that infects domestic pigs with high morbidity and mortality rates. Repeat sequences, which are DNA sequence elements that are repeated more than twice in the genome, play an important role in the ASFV genome. The majority of repeat sequences, however, have not been identified and characterized in a systematic manner. In this study, three types of repeat sequences, including microsatellites, minisatellites and short interspersed nuclear elements (SINEs), were identified in the ASFV genome, and their distribution, structure, function, and evolutionary history were investigated. Most repeat sequences were observed in noncoding regions and at the 5' end of the genome. Noncoding repeat sequences tended to form enhancers, whereas coding repeat sequences had a lower ratio of alpha-helix and beta-sheet and a higher ratio of loop structure and surface amino acids than nonrepeat sequences. In addition, the repeat sequences tended to encode penetrating and antimicrobial peptides. Further analysis of the evolution of repeat sequences revealed that the pan-repeat sequences presented an open state, showing the diversity of repeat sequences. Finally, CpG islands were observed to be negatively correlated with repeat sequence occurrences, suggesting that they may affect the generation of repeat sequences. Overall, this study emphasizes the importance of repeat sequences in ASFVs, and these results can aid in understanding the virus's function and evolution.

Prokaryotic virus host predictor: a Gaussian model for host prediction of prokaryotic viruses in metagenomics.

BACKGROUND: Viruses are ubiquitous biological entities, estimated to be the largest reservoirs of unexplored genetic diversity on Earth. Full functional characterization and annotation of newly discovered viruses requires tools to enable taxonomic assignment, the range of hosts, and biological properties of the virus. Here we focus on prokaryotic viruses, which include phages and archaeal viruses, and for which identifying the viral host is an essential step in characterizing the virus, as the virus relies on the host for survival. Currently, the method for determining the viral host is either to culture the virus, which is low-throughput, time-consuming, and expensive, or to computationally predict the viral hosts, which needs improvements at both accuracy and usability. Here we develop a Gaussian model to predict hosts for prokaryotic viruses with better performances than previous computational methods. RESULTS: We present here Prokaryotic virus Host Predictor (PHP), a software tool using a Gaussian model, to predict hosts for prokaryotic viruses using the differences of k-mer frequencies between viral and host genomic sequences as features. PHP gave a host prediction accuracy of 34% (genus level) on the VirHostMatcher benchmark dataset and a host prediction accuracy of 35% (genus level) on a new dataset containing 671 viruses and 60,105 prokaryotic genomes. The prediction accuracy exceeded that of two alignment-free methods (VirHostMatcher and WIsH, 28-34%, genus level). PHP also outperformed these two alignment-free methods much (24-38% vs 18-20%, genus level) when predicting hosts for prokaryotic viruses which cannot be predicted by the BLAST-based or the CRISPR-spacer-based methods alone. Requiring a minimal score for making predictions (thresholding) and taking the consensus of the top 30 predictions further improved the host prediction accuracy of PHP. CONCLUSIONS: The Prokaryotic virus Host Predictor software tool provides an intuitive and user-friendly API for the Gaussian model described herein. This work will facilitate the rapid identification of hosts for newly identified prokaryotic viruses in metagenomic studies.