Assessment of effectiveness and impact of universal prophylaxis with nirsevimab for prevention of hospitalizations due to respiratory syncytial virus in infants. The NIRSE-GAL study protocol

Nirsevimab has been recently licensed for universal RSV prophylaxis in infants. NIRSE-GAL is a three-year population-based study initiated in Galicia in September 2023. It aims to evaluate nirsevimab effectiveness against RSV-related hospitalizations lower respiratory tract infections (LRTI), severe RSV, all-cause LRTI, and all-cause hospitalization. NIRSE-GAL also aims to estimate nirsevimab impact on primary healthcare use in the short and mid-term, children's wheezing and asthma, and medical prescriptions for RSV. The immunization campaigns will be scheduled based on the expected start week for the RSV season and will last the whole season. Immunization will be offered to: i) infants born during the campaign (seasonal), ii) infants < 6 months at the start of the campaign (catch-up), and iii) infants with high-risk factors, aged 6-24 months at the start of the campaign (high-risk). The follow-up period will start: i) the immunization date for all immunized infants, ii) the start of the campaign, for the non-immunized catch-up or high-risk groups, or iii) the birthdate for the non-immunized seasonal group. Infants will be followed up until outcome occurrence, death, or end of study. Nirsevimab effectiveness will be estimated using Poisson and Cox regression models. Sensitivity and stratified analyses will be undertaken. The number of averted cases and the number needed to immunize will be estimated. Immunization failure and nirsevimab safety will be monitored. NIRSE-GAL was approved by the ethics committee of Galicia (CEIC 2023-377) and registered in ClinicalTrials.gov (ID: NCT06180993). Findings will be mainly shared via peer-reviewed publications and scientific conferences.

COVID-19 InfoVaccines: A WHO-supported educational project to promote COVID-19 vaccination information among professionals and the general population

COVID-19 vaccine uptake varied across countries, in part due to vaccine hesitancy fueled by a lack of trustworthy information. To help health workers provide evidence-based answers to common questions about COVID-19 vaccines and vaccination, and thereby, assist individuals´ decisions on vaccine acceptance, COVID-19 InfoVaccines, a joint WHO-EU project, was launched in February 2021 to support COVID-19 vaccine rollout in 6 Eastern European countries. COVID-19 InfoVaccines was made available in seven languages and shared on social media networks. A total of 262,592 users accessed COVID-19 InfoVaccines.com between February 11, 2021, and January 31st, 2023. The users were most interested in: general questions; vaccine efficacy and duration of protection; vaccine safety; vaccine co-administration, and dose-interval and interchangeability; though the interest in a specific theme varied in function of the epidemiological situation. A total of 118,510 (45.1%) and 46,644 (17.7%) users scrolled up to 35% and 75% of the COVID-19 InfoVaccines webpage, respectively. The average engagement rate was 71.61%. The users accessed COVID-19 InfoVaccines from 231 countries and territories, but the majority were in Ukraine (N = 38,404; 14.6%), Spain (N = 23,327; 8.9%), and Argentina (N = 21,167; 8.1%). Older Facebook users were more interested in COVID-19 information than younger individuals (X2 p-value < .0001). Two hundred twenty-eight videos were shared on YouTube. The average Click-Through-Rate on Facebook was 7.82%, and that on YouTube was 4.4%, with 60 videos having a Click-Through-Rate >5%, falling in the range of average YouTube video Click-Through-Rate (2% - 10%). As misinformation about vaccines and vaccination spreads easily and can negatively impact health-related decisions, initiatives like COVID-19 InfoVaccines are crucial to facilitate access to reliable information.

Short- and mid-term morbidity and primary-care burden due to infant respiratory syncytial virus infection: A Spanish 6-year population-based longitudinal study

BACKGROUND

The morbidity burden of respiratory syncytial virus (RSV) in infants extends beyond hospitalization. Defining the RSV burden before implementing prophylaxis programs is essential for evaluating any potential impact on short- to mid-term morbidity and the utilization of primary healthcare (PHC) and emergency services (ES). We established this reference data using a population-based cohort approach.

METHODS

Infants hospitalized for RSV from January 2016 to March 2023 were matched with non-hospitalized ones based on birthdate and sex. We defined the exposure as severe RSV hospitalization. The main study outcomes were as follows: (1) PHC and ES visits for RSV, categorized using the International Classification of Primary Care codes, (2) prescriptions for respiratory airway obstructive disease, and (3) antibacterial prescriptions. Participants were followed up from 30 days before hospitalization for severe RSV until the outcome occurrence or end of the study. Adjusted incidence rate ratios (IRRs) of the outcomes along with their 95% confidence intervals (CI) were estimated using Poisson regression models. Stratified analyses by type of PHC visit (nurse, pediatrician, or pharmacy) and follow-up period were undertaken. We defined mid-term outcomes as those taking place up to 24 months of follow-up period.

RESULTS

The study included 6626 children (3313 RSV-hospitalized; 3313 non-hospitalized) with a median follow-up of 53.7 months (IQR = 27.9, 69.4). After a 3-month follow-up, severe RSV was associated with a considerable increase in PHC visits for wheezing/asthma (IRR = 4.31, 95% CI: 3.84-4.84), lower respiratory infections (IRR = 4.91, 95% CI: 4.34-5.58), and bronchiolitis (IRR = 4.68, 95% CI: 2.93-7.65). Severe RSV was also associated with more PHC visits for the pediatrician (IRR = 2.00, 95% CI: 1.96-2.05), nurse (IRR = 1.89, 95% CI: 1.75-1.92), hospital emergency (IRR = 2.39, 95% CI: 2.17-2.63), primary healthcare emergency (IRR: 1.54, 95% CI: 1.31-1.82), as well as with important increase in prescriptions for obstructive airway diseases (IRR = 5.98, 95% CI: 5.43-6.60) and antibacterials (IRR = 4.02, 95% CI: 3.38-4.81). All findings remained substantial until 2 years of post-infection.

CONCLUSIONS

Severe RSV infection in infants significantly increases short- to mid-term respiratory morbidity leading to an escalation in healthcare utilization (PHC/ES attendance) and medication prescriptions for up to 2 years afterward. Our approach could be useful in assessing the impact and cost-effectiveness of RSV prevention programs.

Effectiveness and impact of universal prophylaxis with nirsevimab in infants against hospitalisation for respiratory syncytial virus in Galicia, Spain: initial results of a population-based longitudinal study

BACKGROUND

Galicia (Spain) was one of the first regions worldwide to incorporate nirsevimab for universal respiratory syncytial virus (RSV) prophylaxis in infants into its immunisation programme. The NIRSE-GAL longitudinal population-based study aimed to assess nirsevimab effectiveness in preventing hospitalisations (ie, admittance to hospital).

METHODS

The 2023-24 immunisation campaign with nirsevimab in Galicia began on Sept 25, 2023, and concluded on March 31, 2024. The campaign targeted three groups: infants born during the campaign (seasonal group), infants younger than 6 months at the start of the campaign (catch-up group), and infants aged 6-24 months with high-risk factors at the start of the campaign (high-risk group). Infants in the seasonal group were offered immunisation on the first day of life before discharge from hospital. Infants in the catch-up and high-risk groups received electronic appointments to attend a public hospital or health-care centre for nirsevimab administration. For this interim analysis, we used data collected from Sept 25 to Dec 31, 2023, from children born up to Dec 15, 2023. Data were retrieved from public health registries. Nirsevimab effectiveness in preventing RSV-associated lower respiratory tract infection (LRTI) hospitalisations; severe RSV-related LRTI requiring intensive care unit admission, mechanical ventilation, or oxygen support; all-cause LRTI hospitalisations; and all-cause hospitalisations was estimated using adjusted Poisson regression models. Data from five past RSV seasons (2016-17, 2017-18, 2018-19, 2019-20, and 2022-23), excluding the COVID-19 pandemic period, were used to estimate the number of RSV-related LRTI hospitalisations averted along with its IQR. The number needed to immunise to avoid one case in the 2023-24 season was then estimated from the averted cases. Nirsevimab safety was routinely monitored. The NIRSE-GAL study protocol was registered on ClinicalTrials.gov (NCT06180993), and follow-up of participants is ongoing.

FINDINGS

9408 (91·7%) of 10 259 eligible infants in the seasonal and catch-up groups received nirsevimab, including 6220 (89·9%) of 6919 in the seasonal group and 3188 (95·4%) of 3340 in the catch-up group. 360 in the high-risk group were offered nirsevimab, 348 (97%) of whom received it. Only infants in the seasonal and catch-up groups were included in analyses to estimate nirsevimab effectiveness and impact because there were too few events in the high-risk group. In the catch-up and seasonal groups combined, 30 (0·3%) of 9408 infants who received nirsevimab and 16 (1·9%) of 851 who did not receive nirsevimab were hospitalised for RSV-related LRTI, corresponding to an effectiveness of 82·0% (95% CI 65·6-90·2). Effectiveness was 86·9% (69·1-94·2) against severe RSV-related LRTI requiring oxygen support, 69·2% (55·9-78·0) against all-cause LRTI hospitalisations, and 66·2% (56·0-73·7) against all-cause hospitalisations. Nirsevimab effectiveness against other endpoints of severe RSV-related LRTI could not be estimated because of too few events. RSV-related LRTI hospitalisations were reduced by 89·8% (IQR 87·5-90·3), and the number needed to immunise to avoid one RSV-related LRTI hospitalisation was 25 (IQR 24-32). No severe adverse events related to nirsevimab were registered.

INTERPRETATION

Nirsevimab substantially reduced infant hospitalisations for RSV-associated LRTI, severe RSV-associated LRTI requiring oxygen, and all-cause LRTI when given in real-world conditions. These findings offer policy makers and health authorities robust, real-world, population-based evidence to guide the development of strategies for RSV prevention.

FUNDING

Sanofi and AstraZeneca.

TRANSLATION

For the Spanish translation of the abstract see Supplementary Materials section.

Music compensates for altered gene expression in age-related cognitive disorders

Extensive literature has explored the beneficial effects of music in age-related cognitive disorders (ACD), but limited knowledge exists regarding its impact on gene expression. We analyzed transcriptomes of ACD patients and healthy controls, pre-post a music session (n = 60), and main genes/pathways were compared to those dysregulated in mild cognitive impairment (MCI) and Alzheimer's disease (AD) as revealed by a multi-cohort study (n = 1269 MCI/AD and controls). Music was associated with 2.3 times more whole-genome gene expression, particularly on neurodegeneration-related genes, in ACD than in controls. Co-expressed gene-modules and pathways analysis demonstrated that music impacted autophagy, vesicle and endosome organization, biological processes commonly dysregulated in MCI/AD. Notably, the data indicated a strong negative correlation between musically-modified genes/pathways in ACD and those dysregulated in MCI/AD. These findings highlight the compensatory effect of music on genes/biological processes affected in MCI/AD, providing insights into the molecular mechanisms underlying the benefits of music on these disorders.

COVID-19 vaccine effectiveness in children by age groups. A population-based study in Galicia, Spain

BACKGROUND

Studies on vaccine effectiveness (VE) against COVID-19 in the pediatric population are outgoing. We aimed to quantify VE against SARS-CoV-2 in two pediatric age groups, 5-11 and 12-17-year-old, while considering vaccine type, SARS-CoV-2 variant, and duration of protection.

METHODS

A population-based test-negative control study was undertaken in Galicia, Spain. Children 5-11-year-old received the Comirnaty® (Pfizer, US) vaccine, while those aged 12-17-year-old received the Comirnaty® (Pfizer, US) or SpikeVax® (ModernaTX, Inc) vaccine. Participants were categorized into unvaccinated (0 doses or one dose with <14 days since vaccination), partially vaccinated (only one dose with ≥14 days, or two doses with <14 days after the second dose administration), and fully vaccinated (two doses with ≥14 days after the second injection). Adjusted odds ratios (OR) and their 95% confidence intervals (CI) were estimated using multiple logistic regression models. VE was calculated as (1-OR) * 100. Stratified and sensitivity analyses were performed.

RESULTS

In the fully vaccinated 5-11-year-old children, VE against the Omicron variant was 44.1% (95% CI: 38.2%-49.4%). In the fully vaccinated 12-17-year-old individuals, VE was 83.4% (95% CI: 81.2%-85.3%) against Delta and 74.8% (95% CI: 58.5%-84.9%) against Omicron. Comirnaty® and SpikeVax® vaccines showed a similar magnitude of VE against Delta [Comirnaty® VE: 81.9% (95% CI: 79.3%-84.1%) and SpikeVax® VE: 85.3% (95% CI: 81.9%-88.1%)]. Comirnaty® (Pfizer, US; VE: 79.7%; 95% CI: 50.7%-92.4%) showed a slightly higher magnitude of protection against Omicron than SpikeVax® (ModernaTX, Inc), yet with an overlapping CI (VE: 74.3%; 95% CI: 56.6%-84.9%). VE was maintained in all age subgroups in both pediatric populations, but it declined over time.

CONCLUSIONS

In Galicia, mRNA VE was moderate against SARS-CoV-2 infections in the 5-11-year-old populations, but high in older children. VE declined over time, suggesting a potential need for booster dose schedules.

Sex-biased expression of the TLR7 gene in severe COVID-19 patients: Insights from transcriptomics and epigenomics

There is abundant epidemiological data indicating that the incidence of severe cases of coronavirus disease (COVID-19) is significantly higher in males than females worldwide. Moreover, genetic variation at the X-chromosome linked TLR7 gene has been associated with COVID-19 severity. It has been suggested that the sex-biased incidence of COVID-19 might be related to the fact that TLR7 escapes X-chromosome inactivation during early embryogenesis in females, thus encoding a doble dose of its gene product compared to males. We analyzed TLR7 expression in two acute phase cohorts of COVID-19 patients that used two different technological platforms, one of them in a multi-tissue context including saliva, nasal, and blood samples, and a third cohort that included different post-infection timepoints of long-COVID-19 patients. We additionally explored methylation patterns of TLR7 using epigenomic data from an independent cohort of COVID-19 patients stratified by severity and sex. In line with genome-wide association studies, we provide supportive evidence indicating that TLR7 has altered CpG methylation patterns and it is consistently downregulated in males compared to females in the most severe cases of COVID-19.

Effectiveness of COVID-19 vaccine booster in the general population and in subjects with comorbidities. A population-based study in Spain

BACKGROUND

Research on the effectiveness of COVID-19 booster-based vaccine schedule is ongoing and real-world data on vaccine effectiveness (VE) in comorbid patients are limited. We aimed to estimate booster dose VE against SARS-CoV-2 infection and COVID-19 severity in the general population and in comorbid patients.

METHOD

A retrospective test-negative control study was undertaken in Galicia-Spain (December 2020-November 2021). VE and 95% confidence interval (CI) were estimated using multivariate logistic regression models.

RESULTS

1,512,415 (94.13%) negative and 94,334 (5.87%) positive SARS-CoV-2 test results were included. A booster dose of COVID-19 vaccine is associated with substantially higher protection against SARS-CoV-2 infection than vaccination without a booster [VE_boosted = 87% (95%CI: 83%; 89%); VE_non-boosted = 66% (95%CI: 65%; 67%)]. The high VE was observed in all ages, but was more pronounced in subjects older than 65 years. VE against COVID-19 severity was analyzed in a mixed population of boosted and non-boosted individuals and considerable protection was obtained [VE: hospitalization = 72% (95%CI: 68%; 75%); intensive care unit administration = 83% (95%CI: 78%; 88%), in-hospital mortality = 66% (95%CI: 53%; 75%)]. Boosted comorbid patients are more protected against SARS-CoV-2 infection than those who were non-boosted. This was observed in a wide range of major diseases including cancer (81% versus 54%), chronic obstructive pulmonary disease (84% versus 61%), diabetes (84% versus 65%), hypertension (82% versus 65%) and obesity (91% versus 67%), among others.

CONCLUSIONS

A booster dose of COVID-19 vaccine increases the protection against SARS-CoV-2 infection and COVID-19 severity in the general population and in comorbid patients.

Exome-wide association study to identify rare variants influencing COVID-19 outcomes: Results from the Host Genetics Initiative

Host genetics is a key determinant of COVID-19 outcomes. Previously, the COVID-19 Host Genetics Initiative genome-wide association study used common variants to identify multiple loci associated with COVID-19 outcomes. However, variants with the largest impact on COVID-19 outcomes are expected to be rare in the population. Hence, studying rare variants may provide additional insights into disease susceptibility and pathogenesis, thereby informing therapeutics development. Here, we combined whole-exome and whole-genome sequencing from 21 cohorts across 12 countries and performed rare variant exome-wide burden analyses for COVID-19 outcomes. In an analysis of 5,085 severe disease cases and 571,737 controls, we observed that carrying a rare deleterious variant in the SARS-CoV-2 sensor toll-like receptor TLR7 (on chromosome X) was associated with a 5.3-fold increase in severe disease (95% CI: 2.75-10.05, p = 5.41x10-7). This association was consistent across sexes. These results further support TLR7 as a genetic determinant of severe disease and suggest that larger studies on rare variants influencing COVID-19 outcomes could provide additional insights.

A Timeframe for SARS-CoV-2 Genomes: A Proof of Concept for Postmortem Interval Estimations

Establishing the timeframe when a particular virus was circulating in a population could be useful in several areas of biomedical research, including microbiology and legal medicine. Using simulations, we demonstrate that the circulation timeframe of an unknown SARS-CoV-2 genome in a population (hereafter, estimated time of a queried genome [QG]; tE-QG) can be easily predicted using a phylogenetic model based on a robust reference genome database of the virus, and information on their sampling dates. We evaluate several phylogeny-based approaches, including modeling evolutionary (substitution) rates of the SARS-CoV-2 genome (~10-3 substitutions/nucleotide/year) and the mutational (substitutions) differences separating the QGs from the reference genomes (RGs) in the database. Owing to the mutational characteristics of the virus, the present Viral Molecular Clock Dating (VMCD) method covers timeframes going backwards from about a month in the past. The method has very low errors associated to the tE-QG estimates and narrow intervals of tE-QG, both ranging from a few days to a few weeks regardless of the mathematical model used. The SARS-CoV-2 model represents a proof of concept that can be extrapolated to any other microorganism, provided that a robust genome sequence database is available. Besides obvious applications in epidemiology and microbiology investigations, there are several contexts in forensic casework where estimating tE-QG could be useful, including estimation of the postmortem intervals (PMI) and the dating of samples stored in hospital settings.

A multi-tissue study of immune gene expression profiling highlights the key role of the nasal epithelium in COVID-19 severity

Coronavirus Disease-19 (COVID-19) symptoms range from mild to severe illness; the cause for this differential response to infection remains unknown. Unravelling the immune mechanisms acting at different levels of the colonization process might be key to understand these differences. We carried out a multi-tissue (nasal, buccal and blood; n = 156) gene expression analysis of immune-related genes from patients affected by different COVID-19 severities, and healthy controls through the nCounter technology. Mild and asymptomatic cases showed a powerful innate antiviral response in nasal epithelium, characterized by activation of interferon (IFN) pathway and downstream cascades, successfully controlling the infection at local level. In contrast, weak macrophage/monocyte driven innate antiviral response and lack of IFN signalling activity were present in severe cases. Consequently, oral mucosa from severe patients showed signals of viral activity, cell arresting and viral dissemination to the lower respiratory tract, which ultimately could explain the exacerbated innate immune response and impaired adaptative immune responses observed at systemic level. Results from saliva transcriptome suggest that the buccal cavity might play a key role in SARS-CoV-2 infection and dissemination in patients with worse prognosis. Co-expression network analysis adds further support to these findings, by detecting modules specifically correlated with severity involved in the abovementioned biological routes; this analysis also provides new candidate genes that might be tested as biomarkers in future studies. We also found tissue specific severity-related signatures mainly represented by genes involved in the innate immune system and cytokine/chemokine signalling. Local immune response could be key to determine the course of the systemic response and thus COVID-19 severity. Our findings provide a framework to investigate severity host gene biomarkers and pathways that might be relevant to diagnosis, prognosis, and therapy.

Role and Diagnostic Performance of Host Epigenome in Respiratory Morbidity after RSV Infection: The EPIRESVi Study

Background

Respiratory syncytial virus (RSV) infection has been associated with the subsequent development of recurrent wheezing and asthma, although the mechanisms involved are still unknown. We investigate the role of epigenetics in the respiratory morbidity after infection by comparing methylation patterns from children who develop recurrent wheezing (RW-RSV), subsequent asthma (AS-RVS), and those experiencing complete recovery (CR-RSV).

Methods

Prospective, observational study of infants aged < 2 years with RSV respiratory infection admitted to hospital and followed-up after discharge for at least three years. According to their clinical course, patients were categorized into subgroups: RW-RSV (n = 36), AS-RSV (n = 9), and CR-RSV (n = 32). The DNA genome-wide methylation pattern was analyzed in whole blood samples, collected during the acute phase of the infection, using the Illumina Infinium Methylation EPIC BeadChip (850K CpG sites). Differences in methylation were determined through a linear regression model adjusted for age, gender and cell composition.

Results

Patients who developed respiratory sequelae showed a statistically significant higher proportion of NK and CD8T cells (inferred through a deconvolution approach) than those with complete recovery. We identified 5,097 significant differentially methylated positions (DMPs) when comparing RW-RSV and AS-RVS together against CR-RSV. Methylation profiles affect several genes involved in airway inflammation processes. The most significant DMPs were found to be hypomethylated in cases and therefore generally leading to overexpression of affected genes. The lead CpG position (cg24509398) falls at the gene body of EYA3 (P-value = 2.77×10^-10), a tyrosine phosphatase connected with pulmonary vascular remodeling, a key process in the asthma pathology. Logistic regression analysis resulted in a diagnostic epigenetic signature of 3-DMPs (involving genes ZNF2698, LOC102723354 and RPL15/NKIRAS1) that allows to efficiently differentiate sequelae cases from CR-RSV patients (AUC = 1.00). Enrichment pathway analysis reveals the role of the cell cycle checkpoint (FDR P-value = 4.71×10^-2), DNA damage (FDP-value = 2.53×10^-2), and DNA integrity checkpoint (FDR P-value = 2.56×10^-2) in differentiating sequelae from CR-RSV patients.

Conclusions

Epigenetic mechanisms might play a fundamental role in the long-term sequelae after RSV infection, contributing to explain the different phenotypes observed.

CovidPhy: A tool for phylogeographic analysis of SARS-CoV-2 variation

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the pathogen responsible for the coronavirus disease 2019 (COVID-19) pandemic. SARS-CoV-2 genomes have been sequenced massively and worldwide and are now available in different public genome repositories. There is much interest in generating bioinformatic tools capable to analyze and interpret SARS-CoV-2 variation. We have designed CovidPhy (http://covidphy.eu), a web interface that can process SARS-CoV-2 genome sequences in plain fasta text format or provided through identity codes from the Global Initiative on Sharing Avian Influenza Data (GISAID) or GenBank. CovidPhy aggregates information available on the large GISAID database (>1.49 M genomes). Sequences are first aligned against the reference sequence and the interface provides different sources of information, including automatic classification of genomes into a pre-computed phylogeny and phylogeographic information, haplogroup/lineage frequencies, and sequencing variation, indicating also if the genome contains known variants of concern (VOC). Additionally, CovidPhy allows searching for variants and haplotypes introduced by the user and includes a list of genomes that are good candidates for being responsible for large outbreaks worldwide, most likely mediated by important superspreading events, indicating their possible geographic epicenters and their relative impact as recorded in the GISAID database.

Case Report: Everolimus reduced bone turnover markers but showed no clinical benefit in a patient with severe progressive osseous heteroplasia

BACKGROUND

Progressive osseous heteroplasia (POH) is an ultrarare genetic disorder characterized by an inactivating mutation in the GNAS gene that causes heterotopic ossification. Inhibition of the mammalian target of the rapamycin (mTOR) signalling pathway has been proposed as a therapy for progressive bone fibrodysplasia and non-genetic forms of bone heteroplasia. Herein, we describe the impact of using Everolimus as a rescue therapy for an identical twin girl exhibiting an aggressive clinical phenotype of POH.

METHODS

Clinical evaluation of the progression of the disease during Everolimus treatment was performed periodically. Cytokine markers involved in bone metabolism and protein markers related to bone activity were analyzed to explore bone turnover activity.

RESULTS

The patient received Everolimus therapy for 36 weeks. During treatment, no clinical improvement of the disease was perceived. Analysis of biochemical parameters, namely, β-CTX (r ² = -0.576, P-value = 0.016) and PNIP (r ² = -0.598, P-value = 0.011), indicated that bone turnover activity was significantly reduced. Additionally, bone metabolism-related biomarkers showed only a significant positive correlation with PTH levels.

CONCLUSIONS

Everolimus treatment did not modify the clinical progression of the disease in an aggressive form of POH, although an impact on the protein markers studied was observed.

Superspreading in the emergence of COVID-19 variants

Superspreading and variants of concern (VOC) of the human pathogen severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are the main catalyzers of the coronavirus disease 2019 (COVID-19) pandemic. However, measuring their individual impact is challenging. By examining the largest database of SARS-CoV-2 genomes The Global Initiative on Sharing Avian Influenza Data [GISAID; n >1.2 million high-quality (HQ) sequences], we present evidence suggesting that superspreading has had a key role in the epidemiological predominance of VOC. There are clear signatures in the database compatible with large superspreading events (SSEs) coinciding chronologically with the worst epidemiological scenarios triggered by VOC. The data suggest that, without the randomness effect of the genetic drift facilitated by superspreading, new VOC of SARS-CoV-2 would have had more limited chance of success.

Development and Evaluation of the Ancestry Informative Marker Panel of the VISAGE Basic Tool

We detail the development of the ancestry informative single nucleotide polymorphisms (SNPs) panel forming part of the VISAGE Basic Tool (BT), which combines 41 appearance predictive SNPs and 112 ancestry predictive SNPs (three SNPs shared between sets) in one massively parallel sequencing (MPS) multiplex, whereas blood-based age analysis using methylation markers is run in a parallel MPS analysis pipeline. The selection of SNPs for the BT ancestry panel focused on established forensic markers that already have a proven track record of good sequencing performance in MPS, and the overall SNP multiplex scale closely matched that of existing forensic MPS assays. SNPs were chosen to differentiate individuals from the five main continental population groups of Africa, Europe, East Asia, America, and Oceania, extended to include differentiation of individuals from South Asia. From analysis of 1000 Genomes and HGDP-CEPH samples from these six population groups, the BT ancestry panel was shown to have no classification error using the Bayes likelihood calculators of the Snipper online analysis portal. The differentiation power of the component ancestry SNPs of BT was balanced as far as possible to avoid bias in the estimation of co-ancestry proportions in individuals with admixed backgrounds. The balancing process led to very similar cumulative population-specific divergence values for Africa, Europe, America, and Oceania, with East Asia being slightly below average, and South Asia an outlier from the other groups. Comparisons were made of the African, European, and Native American estimated co-ancestry proportions in the six admixed 1000 Genomes populations, using the BT ancestry panel SNPs and 572,000 Affymetrix Human Origins array SNPs. Very similar co-ancestry proportions were observed down to a minimum value of 10%, below which, low-level co-ancestry was not always reliably detected by BT SNPs. The Snipper analysis portal provides a comprehensive population dataset for the BT ancestry panel SNPs, comprising a 520-sample standardised reference dataset; 3445 additional samples from 1000 Genomes, HGDP-CEPH, Simons Foundation and Estonian Biocentre genome diversity projects; and 167 samples of six populations from in-house genotyping of individuals from Middle East, North and East African regions complementing those of the sampling regimes of the other diversity projects.

Identification of a Minimal 3-Transcript Signature to Differentiate Viral from Bacterial Infection from Best Genome-Wide Host RNA Biomarkers: A Multi-Cohort Analysis

The fight against the spread of antibiotic resistance is one of the most important challenges facing health systems worldwide. Given the limitations of current diagnostic methods, the development of fast and accurate tests for the diagnosis of viral and bacterial infections would improve patient management and treatment, as well as contribute to reducing antibiotic misuse in clinical settings. In this scenario, analysis of host transcriptomics constitutes a promising target to develop new diagnostic tests based on the host-specific response to infections. We carried out a multi-cohort meta-analysis of blood transcriptomic data available in public databases, including 11 different studies and 1209 samples from virus- (n = 695) and bacteria- (n = 514) infected patients. We applied a Parallel Regularized Regression Model Search (PReMS) on a set of previously reported genes that distinguished viral from bacterial infection to find a minimum gene expression bio-signature. This strategy allowed us to detect three genes, namely BAFT, ISG15 and DNMT1, that clearly differentiate groups of infection with high accuracy (training set: area under the curve (AUC) 0.86 (sensitivity: 0.81; specificity: 0.87); testing set: AUC 0.87 (sensitivity: 0.82; specificity: 0.86)). BAFT and ISG15 are involved in processes related to immune response, while DNMT1 is related to the preservation of methylation patterns, and its expression is modulated by pathogen infections. We successfully tested this three-transcript signature in the 11 independent studies, demonstrating its high performance under different scenarios. The main advantage of this three-gene signature is the low number of genes needed to differentiate both groups of patient categories.

BCG vaccination improves DTaP immune responses in mice and is associated with lower pertussis incidence in ecological epidemiological studies

BACKGROUND

The Bacillus Calmette-Guérin (BCG), the only vaccine against tuberculosis (TB) currently in use, has shown beneficial effects against unrelated infections and to enhance immune responses to vaccines. However, there is little evidence regarding the influence of BCG vaccination on pertussis.

METHODS

Here, we studied the ability of BCG to improve the immune responses to diphtheria, tetanus, and acellular (DTaP) or whole-cell pertussis (DTwP) vaccination in a mouse model. We included MTBVAC, an experimental live-attenuated vaccine derived from Mycobacterium tuberculosis, in our studies to explore if it presents similar heterologous immunity as BCG. Furthermore, we explored the potential effect of routine BCG vaccination on pertussis incidence worldwide.

FINDINGS

We found that both BCG and MTBVAC when administered before DTaP, triggered Th1 immune responses against diphtheria, tetanus, and pertussis in mice. Immunization with DTaP alone failed to trigger a Th1 response, as measured by the production of IFN-γ. Humoral responses against DTaP antigens were also enhanced by previous immunization with BCG or MTBVAC. Furthermore, exploration of human epidemiological data showed that pertussis incidence was 10-fold lower in countries that use DTaP and BCG compared to countries that use only DTaP.

INTERPRETATION

BCG vaccination may have a beneficial impact on the protection against pertussis conferred by DTaP. Further randomized controlled trials are needed to properly define the impact of BCG on pertussis incidence in a controlled setting. This could be a major finding that would support changes in immunization policies.

FUNDING

This work was supported by the Ministry of "Economía y Competitividad"; European Commission H2020 program, "Gobierno de Aragón"; CIBERES; "Fundação Butantan"; Instituto de Salud Carlos III and "Fondo FEDER".

Pitfalls of barcodes in the study of worldwide SARS-CoV-2 variation and phylodynamics

Analysis of SARS-CoV-2 genome variation using a minimal number of selected informative sites conforming a genetic barcode presents several drawbacks. We show that purely mathematical procedures for site selection should be supervised by known phylogeny (i) to ensure that solid tree branches are represented instead of mutational hotspots with poor phylogeographic proprieties, and (ii) to avoid phylogenetic redundancy. We propose a procedure that prevents information redundancy in site selection by considering the cumulative informativeness of previously selected sites (as a proxy for phylogenetic-based criteria). This procedure demonstrates that, for short barcodes (e.g., 11 sites), there are thousands of informative site combinations that improve previous proposals. We also show that barcodes based on worldwide databases inevitably prioritize variants located at the basal nodes of the phylogeny, such that most representative genomes in these ancestral nodes are no longer in circulation. Consequently, coronavirus phylodynamics cannot be properly captured by universal genomic barcodes because most SARS-CoV-2 variation is generated in geographically restricted areas by the continuous introduction of domestic variants.

Pitfalls of barcodes in the study of worldwide SARS-CoV-2 variation and phylodynamics

Analysis of SARS-CoV-2 genome variation using a minimal number of selected informative sites conforming a genetic barcode presents several drawbacks. We show that purely mathematical procedures for site selection should be supervised by known phylogeny (i) to ensure that solid tree branches are represented instead of mutational hotspots with poor phylogeographic proprieties, and (ii) to avoid phylogenetic redundancy. We propose a procedure that prevents information redundancy in site selection by considering the cumulative informativeness of previously selected sites (as a proxy for phylogenetic-based criteria). This procedure demonstrates that, for short barcodes (e.g., 11 sites), there are thousands of informative site combinations that improve previous proposals. We also show that barcodes based on worldwide databases inevitably prioritize variants located at the basal nodes of the phylogeny, such that most representative genomes in these ancestral nodes are no longer in circulation. Consequently, coronavirus phylodynamics cannot be properly captured by universal genomic barcodes because most SARS-CoV-2 variation is generated in geographically restricted areas by the continuous introduction of domestic variants.

Corrigendum: Case Report: Two Monochorionic Twins With a Critically Different Course of Progressive Osseous Heteroplasia

[This corrects the article DOI: 10.3389/fped.2021.662669.].

Case Report: Two Monochorionic Twins With a Critically Different Course of Progressive Osseus Heteroplasia

Progressive osseous heteroplasia (POH; OMIM 166350) is a rare autosomal-dominant genetic disorder in which extra-skeletal bone forms within skin and muscle tissue. POH is one of the clinical manifestations of an inactivating mutation in the GNAS gene. GNAS gene alterations are difficult matter to address, as GNAS alleles show genetic imprinting and produce several transcript products, and the same mutation may lead to strikingly different phenotypes. Also, most of the publications concerning POH patients are either clinical depictions of a case (or a case series), descriptions of their genetic background, or a tentative correlation of both clinical and molecular findings. Treatment for POH is rarely addressed, and POH still lacks therapeutic options. We describe a unique case of POH in two monochorionic twins, who presented an almost asymptomatic vs. the severe clinical course, despite sharing the same mutation and genetic background. We also report the results of the therapeutic interventions currently available for heterotopic ossification in the patient with the severe course. This article not only critically supports the assumption that the POH course is strongly influenced by factors beyond genetic background but also remarks the lack of options for patients suffering an orphan disease, even after testing drugs with promising in vitro results.

Phylogeography of SARS-CoV-2 pandemic in Spain: a story of multiple introductions, micro-geographic stratification, founder effects, and super-spreaders

Spain has been one of the main global pandemic epicenters for coronavirus disease 2019 (COVID-19). Here, we analyzed >41 000 genomes (including >26 000 high-quality (HQ) genomes) downloaded from the GISAID repository, including 1 245 (922 HQ) sampled in Spain. The aim of this study was to investigate genome variation of novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and reconstruct phylogeographic and transmission patterns in Spain. Phylogeographic analysis suggested at least 34 independent introductions of SARS-CoV-2 to Spain at the beginning of the outbreak. Six lineages spread very successfully in the country, probably favored by super-spreaders, namely, A2a4 (7.8%), A2a5 (38.4%), A2a10 (2.8%), B3a (30.1%), and B9 (8.7%), which accounted for 87.9% of all genomes in the Spanish database. One distinct feature of the Spanish SARS-CoV-2 genomes was the higher frequency of B lineages (39.3%, mainly B3a+B9) than found in any other European country. While B3a, B9, (and an important sub-lineage of A2a5, namely, A2a5c) most likely originated in Spain, the other three haplogroups were imported from other European locations. The B3a strain may have originated in the Basque Country from a B3 ancestor of uncertain geographic origin, whereas B9 likely emerged in Madrid. The time of the most recent common ancestor (TMRCA) of SARS-CoV-2 suggested that the first coronavirus entered the country around 11 February 2020, as estimated from the TMRCA of B3a, the first lineage detected in the country. Moreover, earlier claims that the D614G mutation is associated to higher transmissibility is not consistent with the very high prevalence of COVID-19 in Spain when compared to other countries with lower disease incidence but much higher frequency of this mutation (56.4% in Spain vs. 82.4% in rest of Europe). Instead, the data support a major role of genetic drift in modeling the micro-geographic stratification of virus strains across the country as well as the role of SARS-CoV-2 super-spreaders.

Phylogeography of SARS-CoV-2 pandemic in Spain: a story of multiple introductions, micro-geographic stratification, founder effects, and super-spreaders

Spain has been one of the main global pandemic epicenters for coronavirus disease 2019 (COVID-19). Here, we analyzed >41 000 genomes (including >26 000 high-quality (HQ) genomes) downloaded from the GISAID repository, including 1 245 (922 HQ) sampled in Spain. The aim of this study was to investigate genome variation of novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and reconstruct phylogeographic and transmission patterns in Spain. Phylogeographic analysis suggested at least 34 independent introductions of SARS-CoV-2 to Spain at the beginning of the outbreak. Six lineages spread very successfully in the country, probably favored by super-spreaders, namely, A2a4 (7.8%), A2a5 (38.4%), A2a10 (2.8%), B3a (30.1%), and B9 (8.7%), which accounted for 87.9% of all genomes in the Spanish database. One distinct feature of the Spanish SARS-CoV-2 genomes was the higher frequency of B lineages (39.3%, mainly B3a+B9) than found in any other European country. While B3a, B9, (and an important sub-lineage of A2a5, namely, A2a5c) most likely originated in Spain, the other three haplogroups were imported from other European locations. The B3a strain may have originated in the Basque Country from a B3 ancestor of uncertain geographic origin, whereas B9 likely emerged in Madrid. The time of the most recent common ancestor (TMRCA) of SARS-CoV-2 suggested that the first coronavirus entered the country around 11 February 2020, as estimated from the TMRCA of B3a, the first lineage detected in the country. Moreover, earlier claims that the D614G mutation is associated to higher transmissibility is not consistent with the very high prevalence of COVID-19 in Spain when compared to other countries with lower disease incidence but much higher frequency of this mutation (56.4% in Spain vs. 82.4% in rest of Europe). Instead, the data support a major role of genetic drift in modeling the micro-geographic stratification of virus strains across the country as well as the role of SARS-CoV-2 super-spreaders.

Mapping genome variation of SARS-CoV-2 worldwide highlights the impact of COVID-19 super-spreaders

The human pathogen severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the major pandemic of the twenty-first century. We analyzed more than 4700 SARS-CoV-2 genomes and associated metadata retrieved from public repositories. SARS-CoV-2 sequences have a high sequence identity (>99.9%), which drops to >96% when compared to bat coronavirus genome. We built a mutation-annotated reference SARS-CoV-2 phylogeny with two main macro-haplogroups, A and B, both of Asian origin, and more than 160 sub-branches representing virus strains of variable geographical origins worldwide, revealing a rather uniform mutation occurrence along branches that could have implications for diagnostics and the design of future vaccines. Identification of the root of SARS-CoV-2 genomes is not without problems, owing to conflicting interpretations derived from either using the bat coronavirus genomes as an outgroup or relying on the sampling chronology of the SARS-CoV-2 genomes and TMRCA estimates; however, the overall scenario favors haplogroup A as the ancestral node. Phylogenetic analysis indicates a TMRCA for SARS-CoV-2 genomes dating to November 12, 2019, thus matching epidemiological records. Sub-haplogroup A2 most likely originated in Europe from an Asian ancestor and gave rise to subclade A2a, which represents the major non-Asian outbreak, especially in Africa and Europe. Multiple founder effect episodes, most likely associated with super-spreader hosts, might explain COVID-19 pandemic to a large extent.

Seroprevalence of SARS-CoV-2 Among Pediatric Healthcare Workers in Spain

Spain is one of the countries most severely affected by the SARS-CoV-2 pandemic, with almost 190,000 cases as of April 18, 2020. As healthcare workers (HCW) are one of the groups hardest hit by the infection, it is important to know the seroprevalence of antibodies against SARS-CoV-2 in pediatric departments. We performed 175 immunoglobulin (Ig)M and IgG immunochromatographic rapid tests in the personnel working at the Pediatric Department of the Hospital Clínico Universitario of Santiago de Compostela (Spain), including pediatricians, residents, nurses, and other staff, on days 31-33 since the lockdown started. Seven out of the 175 tests were positive, including four for IgM and three for IgG, leading to a seroprevalence of 4.0% (95% CI: 1.1-6.9%). Only one of them had symptoms at the time of testing (sore throat). All seropositive cases yielded negative RT-PCR of the upper and lower respiratory tract. This is the first SARS-CoV-2 serological survey among HCWs reported in Spain. Notwithstanding the test limitations, our results reveal that personal protection policy and lockdown measures have been effective to limit population exposure. The low seroprevalence rate poses a significant challenge for the next strategic steps of pandemic control.

A qPCR expression assay of IFI44L gene differentiates viral from bacterial infections in febrile children

The diagnosis of bacterial infections in hospital settings is currently performed using bacterial culture from sterile site, but they are lengthy and limited. Transcriptomic biomarkers are becoming promising tools for diagnosis with potential applicability in clinical settings. We evaluated a RT-qPCR assay for a 2-transcript host expression signature (FAM89A and IFI44L genes) inferred from microarray data that allow to differentiate between viral and bacterial infection in febrile children. This assay was able to discriminate viral from bacterial infections (P-value = 1.04 × 10-4; AUC = 92.2%; sensitivity = 90.9%; specificity = 85.7%) and showed very high reproducibility regardless of the reference gene(s) used to normalize the data. Unexpectedly, the monogenic IFI44L expression signature yielded better results than those obtained from the 2-transcript test (P-value = 3.59 × 10-5; AUC = 94.1%; sensitivity = 90.9%; specificity = 92.8%). We validated this IFI44L signature in previously published microarray and whole-transcriptome data from patients affected by different types of viral and bacterial infections, confirming that this gene alone differentiates between both groups, thus saving time, effort, and costs. Herein, we demonstrate that host expression microarray data can be successfully translated into a fast, highly accurate and relatively inexpensive in vitro assay that could be implemented in the clinical routine.

Ancestry patterns inferred from massive RNA-seq data

There is a growing body of evidence suggesting that patterns of gene expression vary within and between human populations. However, the impact of this variation in human diseases has been poorly explored, in part owing to the lack of a standardized protocol to estimate biogeographical ancestry from gene expression studies. Here we examine several studies that provide new solid evidence indicating that the ancestral background of individuals impacts gene expression patterns. Next, we test a procedure to infer genetic ancestry from RNA-seq data in 25 data sets where information on ethnicity was reported. Genome data of reference continental populations retrieved from The 1000 Genomes Project were used for comparisons. Remarkably, only eight out of 25 data sets passed FastQC default filters. We demonstrate that, for these eight population sets, the ancestral background of donors could be inferred very efficiently, even in data sets including samples with complex patterns of admixture (e.g., American-admixed populations). For most of the gene expression data sets of suboptimal quality, ancestral inference yielded odd patterns. The present study thus brings a cautionary note for gene expression studies highlighting the importance to control for the potential confounding effect of ancestral genetic background.

Lifestyle and comorbid conditions as risk factors for community-acquired pneumonia in outpatient adults (NEUMO-ES-RISK project)

Introduction

Information about community-acquired pneumonia (CAP) risk in primary care is limited. We assess different lifestyle and comorbid conditions as risk factors (RF) for CAP in adults in primary care.

Methods

A retrospective-observational-controlled study was designed. Adult CAP cases diagnosed at primary care in Spain between 2009 and 2013 were retrieved using the National Surveillance System of Primary Care Data (BiFAP). Age-matched and sex-matched controls were selected by incidence density sampling (ratio 2:1). Associations are presented as percentages and OR. Binomial regression models were constructed to avoid bias effects.

Results

51 139 patients and 102 372 controls were compared. Mean age (SD) was 61.4 (19.9) years. RF more significantly linked to CAP were: HIV (OR [95% CI]: 5.21 [4.35 to 6.27]), chronic obstructive pulmonary disease (COPD) (2.97 [2.84 to 3.12]), asthma (2.16 [2.07,2.26]), smoking (1.96 [1.91 to 2.02]) and poor dental hygiene (1.45 [1.41 to 1.49]). Average prevalence of any RF was 82.2% in cases and 69.2% in controls (2.05 [2.00 to 2.10]). CAP rate increased with the accumulation of RF and age: risk associated with 1RF was 1.42 (1.37 to 1.47) in 18–60-year-old individuals vs 1.57 (1.49 to 1.66) in >60 years of age, with 2RF 1.88 (1.80 to 1.97) vs 2.35 (2.23, 2.48) and with ≥ 3 RF 3.11 (2.95, 3.30) vs 4.34 (4.13 to 4.57).

Discussion

Prevalence of RF in adult CAP in primary care is high. Main RFs associated are HIV, COPD, asthma, smoking and poor dental hygiene. Our risk stacking results could help clinicians identify patients at higher risk of pneumonia.

Rotavirus infection beyond the gut

The landscape of rotavirus (RV) infection has changed substantially in recent years. Autoimmune triggering has been added to clinical spectrum of this pathology, which is now known to be much broader than diarrhea. The impact of RV vaccines in these other conditions is becoming a growing field of research. The importance of host genetic background in RV susceptibility has been revealed, therefore increasing our understanding of vaccine effectiveness and giving some clues about the limited efficacy of RV vaccines in low-income settings. Also, interaction of RV with intestinal microbiota seems to play a key role in the process of infection vaccine effect. This article reviews current findings on the extraintestinal impact of RV infection and their widening clinical picture, and the recently described mechanisms of host susceptibility to infection and vaccine effectiveness. RV infection is a systemic disease with clinical and pathophysiological implications beyond the gut. We propose an “iceberg” model for this pathology with almost hidden clinical implications away from the gastrointestinal tract and eventually triggering the development of autoimmune diseases. Impact of current vaccines is being influenced by host genetics and gut microbiota interactions and these factors must be taken into account in the development of public health programs.

The peopling of South America and the trans-Andean gene flow of the first settlers

Genetic and archaeological data indicate that the initial Paleoindian settlers of South America followed two entry routes separated by the Andes and the Amazon rainforest. The interactions between these paths and their impact on the peopling of South America remain unclear. Analysis of genetic variation in the Peruvian Andes and regions located south of the Amazon River might provide clues on this issue. We analyzed mitochondrial DNA variation at different Andean locations and >360,000 autosomal SNPs from 28 Native American ethnic groups to evaluate different trans-Andean demographic scenarios. Our data reveal that the Peruvian Altiplano was an important enclave for early Paleoindian expansions and point to a genetic continuity in the Andes until recent times, which was only marginally affected by gene flow from the Amazonian lowlands. Genomic variation shows a good fit with the archaeological evidence, indicating that the genetic interactions between the descendants of the settlers that followed the Pacific and Atlantic routes were extremely limited.

Life-threatening infections in children in Europe (the EUCLIDS Project): a prospective cohort study

BACKGROUND

Sepsis and severe focal infections represent a substantial disease burden in children admitted to hospital. We aimed to understand the burden of disease and outcomes in children with life-threatening bacterial infections in Europe.

METHODS

The European Union Childhood Life-threatening Infectious Disease Study (EUCLIDS) was a prospective, multicentre, cohort study done in six countries in Europe. Patients aged 1 month to 18 years with sepsis (or suspected sepsis) or severe focal infections, admitted to 98 participating hospitals in the UK, Austria, Germany, Lithuania, Spain, and the Netherlands were prospectively recruited between July 1, 2012, and Dec 31, 2015. To assess disease burden and outcomes, we collected demographic and clinical data using a secured web-based platform and obtained microbiological data using locally available clinical diagnostic procedures.

FINDINGS

2844 patients were recruited and included in the analysis. 1512 (53·2%) of 2841 patients were male and median age was 39·1 months (IQR 12·4-93·9). 1229 (43·2%) patients had sepsis and 1615 (56·8%) had severe focal infections. Patients diagnosed with sepsis had a median age of 27·6 months (IQR 9·0-80·2), whereas those diagnosed with severe focal infections had a median age of 46·5 months (15·8-100·4; p<0·0001). Of 2844 patients in the entire cohort, the main clinical syndromes were pneumonia (511 [18·0%] patients), CNS infection (469 [16·5%]), and skin and soft tissue infection (247 [8·7%]). The causal microorganism was identified in 1359 (47·8%) children, with the most prevalent ones being Neisseria meningitidis (in 259 [9·1%] patients), followed by Staphylococcus aureus (in 222 [7·8%]), Streptococcus pneumoniae (in 219 [7·7%]), and group A streptococcus (in 162 [5·7%]). 1070 (37·6%) patients required admission to a paediatric intensive care unit. Of 2469 patients with outcome data, 57 (2·2%) deaths occurred: seven were in patients with severe focal infections and 50 in those with sepsis.

INTERPRETATION

Mortality in children admitted to hospital for sepsis or severe focal infections is low in Europe. The disease burden is mainly in children younger than 5 years and is largely due to vaccine-preventable meningococcal and pneumococcal infections. Despite the availability and application of clinical procedures for microbiological diagnosis, the causative organism remained unidentified in approximately 50% of patients.

FUNDING

European Union's Seventh Framework programme.

Rotavirus intestinal infection induces an oral mucosa cytokine response

INTRODUCTION

Salivary glands are known immune effector sites and considered to be part of the whole mucosal immune system. The aim of the present study was to assess the salivary immune response to rotavirus (RV) infection through the analysis of the cytokine immune profile in saliva.

MATERIAL AND METHODS

A prospective comparative study of serial saliva samples from 27 RV-infected patients (sampled upon admission to the hospital during acute phase and at convalescence-i.e. at least three months after recovery) and 36 healthy controls was performed. Concentrations of 11 salivary cytokines (IFN-γ, IFN-α2, IL-1β, IL-6, IL-8, IL-10, IL-15, IL12p70, TNF-α, IFN-λ1, IL-22) were determined. Cytokine levels were compared between healthy controls acute infection and convalescence. The correlation between clinical data and salivary cytokine profile in infected children was assessed.

RESULTS

The salivary cytokine profile changes significantly in response to acute RV infection. In RV-infected patients, IL-22 levels were increased in the acute phase with respect to convalescence (P-value < 0.001). Comparisons between infected and control group showed significant differences in salivary IFN-α2, IL-1β, IL-6, IL-8, IL-10 and IL-22. Although acute-phase levels of IL-12, IL-10, IL-6 and IFN-γ showed nominal association with Vesikari's severity, this trend did not reach statistical significance after multiple test adjustment.

CONCLUSIONS

RV infection induces a host salivary immune response, indicating that immune mucosal response to RV infection is not confined to the intestinal mucosa. Our data point to a whole mucosal implication in the RV infection as a result of the integrative mucosal immune response, and suggest the salivary gland as effector site for RV infection.

The geographic mosaic of Ecuadorian Y-chromosome ancestry

Ecuadorians originated from a complex mixture of Native American indigenous people with Europeans and Africans. We analyzed Y-chromosome STRs (Y-STRs) in a sample of 415 Ecuadorians (145 using the AmpFlSTR^® Yfiler™ system [Life Technologies, USA] and 270 using the PowerPlex^®Y23 system [Promega Corp., USA]; hereafter Yfiler and PPY23, respectively) representing three main ecological continental regions of the country, namely Amazon rainforest, Andes, and Pacific coast. Diversity values are high in the three regions, and the PPY23 exhibits higher discrimination power than the Yfiler set. While summary statistics, AMOVA, and R_ST distances show low to moderate levels of population stratification, inferred ancestry derived from Y-STRs reveal clear patterns of geographic variation. The major ancestry in Ecuadorian males is European (61%), followed by an important Native American component (34%); whereas the African ancestry (5%) is mainly concentrated in the Northwest corner of the country. We conclude that classical procedures for measuring population stratification do not have the desirable sensitivity. Statistical inference of ancestry from Y-STRS is a satisfactory alternative for revealing patterns of spatial variation that would pass unnoticed when using popular statistical summary indices.

Whole Exome Sequencing reveals new candidate genes in host genomic susceptibility to Respiratory Syncytial Virus Disease

Respiratory syncytial virus (RSV) is an important cause of serious lower respiratory tract disease in infants. Several studies have shown evidence pointing to the genome of the host as an important factor determining susceptibility to respiratory disease caused by RSV. We sequenced the complete exomes of 54 patients infected by RSV that needed hospitalization due to development of severe bronchiolitis. The Iberian sample (IBS) from The 1000 Genomes Project (1000G) was used as control group; all the association results were pseudo-replicated using other 1000G-European controls and Spanish controls. The study points to SNP rs199665292 in the olfactory receptor (OR) gene OR13C5 as the best candidate variant (P-value = 1.16 × 10-12; OR = 5.56). Genetic variants at HLA genes (HLA-DQA1, HLA-DPB1), and in the mucin 4 gene (MUC4) also emerge as susceptibility candidates. By collapsing rare variants in genes and weighing by pathogenicity, we obtained confirmatory signals of association in the OR gene OR8U1/OR8U8, the taste receptor TAS2R19, and another mucin gene (MUC6). Overall, we identified new predisposition variants and genes related to RSV infection. Of special interest is the association of RSV to olfactory and taste receptors; this finding is in line with recent evidence pointing to their role in viral infectious diseases.

Phylogeographic and genome-wide investigations of Vietnam ethnic groups reveal signatures of complex historical demographic movements

The territory of present-day Vietnam was the cradle of one of the world’s earliest civilizations, and one of the first world regions to develop agriculture. We analyzed the mitochondrial DNA (mtDNA) complete control region of six ethnic groups and the mitogenomes from Vietnamese in The 1000 Genomes Project (1000G). Genome-wide data from 1000G (~55k SNPs) were also investigated to explore different demographic scenarios. All Vietnamese carry South East Asian (SEA) haplotypes, which show a moderate geographic and ethnic stratification, with the Mong constituting the most distinctive group. Two new mtDNA clades (M7b1a1f1 and F1f1) point to historical gene flow between the Vietnamese and other neighboring countries. Bayesian-based inferences indicate a time-deep and continuous population growth of Vietnamese, although with some exceptions. The dramatic population decrease experienced by the Cham 700 years ago (ya) fits well with the Nam tiến (“southern expansion”) southwards from their original heartland in the Red River Delta. Autosomal SNPs consistently point to important historical gene flow within mainland SEA, and add support to a main admixture event occurring between Chinese and a southern Asian ancestral composite (mainly represented by the Malay). This admixture event occurred ~800 ya, again coinciding with the Nam tiến.

Salivary epidermal growth factor correlates with hospitalization length in rotavirus infection

BACKGROUND

The IFI27 interferon gene expression has been found to be largely increased in rotavirus (RV)-infected patients. IFI27 gene encodes for a protein of unknown function, very recently linked to epidermal proliferation and related to the epidermal growth factor (EGF) protein. The EGF is a low-molecular-weight polypeptide that is mainly produced by submandibular and parotid glands, and it plays an important physiological role in the maintenance of oro-esophageal and gastric tissue integrity. Our aim was to determine salivary EGF levels in RV-infected patients in order to establish its potential relationship with IFI27 increased expression and EGF-mediated mucosal protection in RV infection.

METHODS

We conducted a prospective comparative study using saliva samples from 27 infants infected with RV (sampled at recruitment during hospital admission and at convalescence, i.e. at least 3 months after recovery) and from 36 healthy control children.

RESULTS

Median (SD) EGF salivary concentration was 777 (529) pg/ml in RV-infected group at acute phase and 356 (242) pg/m at convalescence, while it was 337 (119) pg/ml in the healthy control group. A significant association was found between EGF levels and hospitalization length of stay (P-value = 0.022; r2 = -0.63).

CONCLUSIONS

The salivary levels of EGF are significantly increased during the acute phase of natural RV infection, and relate to length of hospitalization. Further assessment of this non-invasive biomarker in RV disease is warranted.

Genomic continuity of Argentinean Mennonites

Mennonites are Anabaptist communities that originated in Central Europe about 500 years ago. They initially migrated to different European countries, and in the early 18^th century they established their first communities in North America, from where they moved to other American regions. We aimed to analyze an Argentinean Mennonite congregation from a genome-wide perspective by way of investigating >580.000 autosomal SNPs. Several analyses show that Argentinean Mennonites have European ancestry without signatures of admixture with other non-European American populations. Among the worldwide datasets used for population comparison, the CEU, which is the best-subrogated Central European population existing in The 1000 Genome Project, is the dataset showing the closest genome affinity to the Mennonites. When compared to other European population samples, the Mennonites show higher inbreeding coefficient values. Argentinean Mennonites show signatures of genetic continuity with no evidence of admixture with Americans of Native American or sub-Saharan African ancestry. Their genome indicates the existence of an increased endogamy compared to other Europeans most likely mirroring their lifestyle that involve small communities and historical consanguineous marriages.

Incidence and risk factor prevalence of community-acquired pneumonia in adults in primary care in Spain (NEUMO-ES-RISK project)

Background

Community-acquired pneumonia (CAP) is a major cause of morbidity and mortality in adults even in developed countries. Several lifestyle factors and comorbidities have been linked to an increased risk, although their prevalence has not been well documented in the primary care setting. The aim of this study is to assess the incidence, risk factor and comorbid conditions distribution of CAP in adults in primary care in Spain.

Methods

Retrospective observational study in adults (>18 years-old) with CAP diagnosed and attended at primary care in Spain between 2009 and 2013, using the Computerized Database for Pharmacoepidemiological Studies in Primary Care (BIFAP).

Results

Twenty-eight thousand four hundred thirteen patient records were retrieved and analyzed. Mean age (standard deviation): 60.5 (20.3) years, 51.7 % males. Global incidence of CAP in adults was estimated at 4.63 per 1000 persons/year. CAP incidence increased progressively with age, ranging from a 1.98 at 18–20 years of age to 23.74 in patients over 90 years of age. According to sex, global CAP incidence was slightly higher in males (5.04) than females (4.26); CAP incidence from 18 to 65 year-olds up was comparable between males (range: 2.18–5.75) and females (range: 1.47–5.21), whereas from 65 years of age, CAP incidence was noticeable higher in males (range: 7.06–36.93) than in females (range: 5.43–19.62). Average prevalence of risk factors was 71.3 %, which increased with age, doubling the risk in males by the age of 75 (females 20 % vs males 40 %). From 55 years of age, at least one risk factor was identified in 85.7 % of cases: one risk factor (23.8 %), two risk factors (23.4 %), three or more risk factors (38.5 %). Major risk factors were: metabolic disease (27.4 %), cardiovascular disease (17.8 %) and diabetes (15.5 %).

Conclusions

The annual incidence of CAP in primary care adults in Spain is high, comparable between males and females up to 65 years of age, but clearly increasing in males from that age. CAP risk increases with age and doubles in males older than 75 years. The majority of CAP cases in patients over 55 years of age is associated to at least one risk factor. The main risk factors associated were metabolic disease, cardiovascular disease, and diabetes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-016-1974-4) contains supplementary material, which is available to authorized users.

Development and Validation of a New Clinical Scale for Infants with Acute Respiratory Infection: The ReSVinet Scale

Background and Aims

A properly validated scoring system allowing objective categorization of infants with acute respiratory infections (ARIs), avoiding the need for in-person assessment and that could also be used by non-health professionals is currently not available. We aimed to develop a new clinical assessment scale meeting these specifications.

Methods

We designed a clinical scale (ReSVinet scale) based on seven parameters (feeding intolerance, medical intervention, respiratory difficulty, respiratory frequency, apnoea, general condition, fever) that were assigned different values (from 0 to 3) for a total of 20 points.170 children under two years of age with ARI were assessed independently by three pediatricians using this scale. Parents also evaluated their offspring with an adapted version of the scale in a subset of 61 cases. The scale was tested for internal consistency (Cronbach’s alpha), Pearson correlation coefficient for the items in the scale, inter-observer reliability (kappa index) and floor-ceiling effect.

Results

Internal consistency was good for all the observers, with the lowest Cronbach’s alpha being 0.72. There was a strong correlation between the investigators (r-value ranged 0.76–0.83) and also between the results obtained by the parents and the investigators(r = 0.73). Light’s kappa for the observations of the three investigators was 0.74. Weighted kappa in the group evaluated by the parents was 0.73. The final score was correlated with length of hospital stay, PICU admission and Wood-Downes Score.

Conclusions

The ReSVinet scale may be useful and reliable in the evaluation of infants with ARI, particularly acute bronchiolitis, even with data obtained from medical records and when employed by parents. Although further studies are necessary, ReSVinet scale already complies with more score validation criteria than the vast majority of the alternatives currently available and used in the clinical practice.

Does Viral Co-Infection Influence the Severity of Acute Respiratory Infection in Children?

Background

Multiple viruses are often detected in children with respiratory infection but the significance of co-infection in pathogenesis, severity and outcome is unclear.

Objectives

To correlate the presence of viral co-infection with clinical phenotype in children admitted with acute respiratory infections (ARI).

Methods

We collected detailed clinical information on severity for children admitted with ARI as part of a Spanish prospective multicenter study (GENDRES network) between 2011–2013. A nested polymerase chain reaction (PCR) approach was used to detect respiratory viruses in respiratory secretions. Findings were compared to an independent cohort collected in the UK.

Results

204 children were recruited in the main cohort and 97 in the replication cohort. The number of detected viruses did not correlate with any markers of severity. However, bacterial superinfection was associated with increased severity (OR: 4.356; P-value = 0.005), PICU admission (OR: 3.342; P-value = 0.006), higher clinical score (1.988; P-value = 0.002) respiratory support requirement (OR: 7.484; P-value < 0.001) and longer hospital length of stay (OR: 1.468; P-value < 0.001). In addition, pneumococcal vaccination was found to be a protective factor in terms of degree of respiratory distress (OR: 2.917; P-value = 0.035), PICU admission (OR: 0.301; P-value = 0.011), lower clinical score (-1.499; P-value = 0.021) respiratory support requirement (OR: 0.324; P-value = 0.016) and oxygen necessity (OR: 0.328; P-value = 0.001). All these findings were replicated in the UK cohort.

Conclusion

The presence of more than one virus in hospitalized children with ARI is very frequent but it does not seem to have a major clinical impact in terms of severity. However bacterial superinfection increases the severity of the disease course. On the contrary, pneumococcal vaccination plays a protective role.

Mapping the genomic mosaic of two 'Afro-Bolivians' from the isolated Yungas valleys

Background

Unraveling the ancestry of ‘Afro-American’ communities is hampered by the complex demographic processes that took place during the Transatlantic Slave Trade (TAST) and the (post-)colonization periods. ‘Afro-Bolivians’ from the subtropical Yungas valleys constitute small and isolated communities that live surrounded by the predominant Native American community of Bolivia. By genotyping >580,000 SNPs in two ‘Afro-Bolivians’, and comparing these genomic profiles with data compiled from more than 57 African groups and other reference ancestral populations (n = 1,161 in total), we aimed to disentangle the complex admixture processes undergone by ‘Afro-Bolivians’.

Results

The data indicate that these two genomes constitute a complex mosaic of ancestries that is approximately 80 % of recent African origin; the remaining ~20 % being European and Native American. West-Central Africa contributed most of the African ancestry to ‘Afro-Bolivians’, and this component is related to populations living along the Atlantic coast (i.e. Senegal, Ghana, Nigeria). Using tract length distribution of genomic segments attributable to distinct ancestries, we could date the time of admixture in about 400 years ago. This time coincides with the maximum importation of slaves to Bolivia to compensate the diminishing indigenous labor force needed for the development of the National Mint of Potosí.

Conclusions

Overall, the data indicate that the genome of ‘Afro-Bolivians’ was shaped by a complex process of admixture occurring in America among individuals originating in different West-Central African populations; their genomic mosaics received additional contributions of Europeans and local Native Americans (e.g. Aymaras).

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2520-x) contains supplementary material, which is available to authorized users.

Bacteremia in Children Hospitalized with Respiratory Syncytial Virus Infection

BACKGROUND

The risk of bacteremia is considered low in children with acute bronchiolitis. However the rate of occult bacteremia in infants with RSV infection is not well established. The aim was to determine the actual rate and predictive factors of bacteremia in children admitted to hospital due to confirmed RSV acute respiratory illness (ARI), using both conventional culture and molecular techniques.

METHODS

A prospective multicenter study (GENDRES-network) was conducted between 2011-2013 in children under the age of two admitted to hospital because of an ARI. Among those RSV-positive, bacterial presence in blood was assessed using PCR for Meningococcus, Streptococcus pneumoniae, Haemophilus influenzae, Streptococcus pyogenes, Klebsiella pneumoniae, Pseudomonas aeruginosa, Escherichia coli, and Staphylococcus aureus, in addition to conventional cultures.

RESULTS

66 children with positive RSV respiratory illness were included. In 10.6% patients, bacterial presence was detected: H. influenzae (n = 4) and S. pneumoniae (n = 2). In those patients with bacteremia, there was a previous suspicion of bacterial superinfection and had received empirical antibiotic treatment 6 out of 7 (85.7%) patients. There were significant differences in terms of severity between children with positive bacterial PCR and those with negative results: PICU admission (100% vs. 50%, P-value = 0.015); respiratory support necessity (100% vs. 18.6%, P-value < 0.001); Wood-Downes score (mean = 8.7 vs. 4.8 points, P-value < 0.001); GENVIP scale (mean = 17 vs. 10.1, P-value < 0.001); and length of hospitalization (mean = 12.1 vs. 7.5 days, P-value = 0.007).

CONCLUSION

Bacteremia is not frequent in infants hospitalized with RSV respiratory infection, however, it should be considered in the most severe cases.

Charting the Y-chromosome ancestry of present-day Argentinean Mennonites

Old Order Mennonite communities initially arose in Northern Europe (centered in the Netherlands) and derived from the Anabaptist movement of the 16th century. Mennonites migrated to the New World in the early 18th century, first to North America, and more recently to Mesoamerica and South America. We analyzed Y-chromosome short tandem repeats (STRs) and single nucleotide polymorphisms in males from a community of Mennonites, 'La Nueva Esperanza', which arrived to Argentina in 1985 from colonies in Bolivia and Mexico. Molecular diversity indices coupled with demographic simulations show that Mennonites have a reduced variability when compared with local Argentinean populations and 69 European population samples. Mennonite Y-STR haplotypes were mainly observed in Central Europe. In agreement, multidimensional scaling analyses based on RST genetic distances indicate that Mennonite Y-chromosomes are closely related to Central/Northern Europeans (the Netherlands, Switzerland and Denmark). In addition, statistical inferences made on the most likely geographic origin of Y-chromosome haplotypes point more specifically to the Netherlands as the populations that best represent the majority of the Mennonite Y-chromosomes. Overall, Y-chromosome variation of Mennonites shows the signatures of moderate reduction of variability when compared with source populations, which is in good agreement with their lifestyle in small endogamous demes. These genetic singularities could also help to understand disease conditions that are more prevalent among Mennonites.