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.