Analysis of Y-chromosome STRs in Chile confirms an extensive introgression of European male lineages in urban populations

We analyzed the Y chromosome haplotypes (Yfiler) of 978 non-related Chilean males grouped in five sampling regions (Iquique, Santiago de Chile, Concepción, Temuco and Punta Arenas) covering main geo-political regions. Overall, 803 different haplotypes and 688 singletons were observed. Molecular diversity was moderately lower than in other neighboring countries (e.g. Argentina); and AMOVA analysis on Y-STR haplotypes showed that among variation within Chile accounted for only 0.25% of the total variation. Punta Arenas, in the southern cone, showed the lowest haplotype diversity, and discrimination capacity, and also the highest matching probability of the five Chilean samples, probably reflecting its more marked geographic isolation compared to the other regions. Multidimensional scaling (MDS) analysis based on RST genetic distances suggested a close proximity of Chilean Y-chromosome profiles to European ones. Consistently, haplogroups inferred from Y-STR profiles revealed that the Native American component constituted only 8% of all the haplotypes, and this component ranged from 5% in the Centre of the country to 9-10% in the South and 13% in the North, which is in good agreement with the distribution of Native American communities in these regions. AMOVA computed on inferred haplogroups confirmed the very low among variation observed in Chilean populations. The present project provides the first Chilean dataset to the international Y-chromosome STR Haplotype Reference Database (YHRD) and it is also the first reference database for Y-chromosome forensic casework of the country.

The complete mitogenome of a 500-year-old Inca child mummy

In 1985, a frozen mummy was found in Cerro Aconcagua (Argentina). Archaeological studies identified the mummy as a seven-year-old Inca sacrifice victim who lived >500 years ago, at the time of the expansion of the Inca Empire towards the southern cone. The sequence of its entire mitogenome was obtained. After querying a large worldwide database of mitogenomes (>28,000) we found that the Inca haplotype belonged to a branch of haplogroup C1b (C1bi) that has not yet been identified in modern Native Americans. The expansion of C1b into the Americas, as estimated using 203 C1b mitogenomes, dates to the initial Paleoindian settlements (~18.3 thousand years ago [kya]); however, its internal variation differs between Mesoamerica and South America. By querying large databases of control region haplotypes (>150,000), we found only a few C1bi members in Peru and Bolivia (e.g. Aymaras), including one haplotype retrieved from ancient DNA of an individual belonging to the Wari Empire (Peruvian Andes). Overall, the results suggest that the profile of the mummy represents a very rare sub-clade that arose 14.3 (5–23.6) kya and could have been more frequent in the past. A Peruvian Inca origin for present-day C1bi haplotypes would satisfy both the genetic and paleo-anthropological findings.

Viral Co-Infections in Pediatric Patients Hospitalized with Lower Tract Acute Respiratory Infections

BACKGROUND

Molecular techniques can often reveal a broader range of pathogens in respiratory infections. We aim to investigate the prevalence and age pattern of viral co-infection in children hospitalized with lower tract acute respiratory infection (LT-ARI), using molecular techniques.

METHODS

A nested polymerase chain reaction approach was used to detect Influenza (A, B), metapneumovirus, respiratory syncytial virus (RSV), parainfluenza (1-4), rhinovirus, adenovirus (A-F), bocavirus and coronaviruses (NL63, 229E, OC43) in respiratory samples of children with acute respiratory infection prospectively admitted to any of the GENDRES network hospitals between 2011-2013. The results were corroborated in an independent cohort collected in the UK.

RESULTS

A total of 204 and 97 nasopharyngeal samples were collected in the GENDRES and UK cohorts, respectively. In both cohorts, RSV was the most frequent pathogen (52.9% and 36.1% of the cohorts, respectively). Co-infection with multiple viruses was found in 92 samples (45.1%) and 29 samples (29.9%), respectively; this was most frequent in the 12-24 months age group. The most frequently observed co-infection patterns were RSV-Rhinovirus (23 patients, 11.3%, GENDRES cohort) and RSV-bocavirus / bocavirus-influenza (5 patients, 5.2%, UK cohort).

CONCLUSION

The presence of more than one virus in pediatric patients admitted to hospital with LT-ARI is very frequent and seems to peak at 12-24 months of age. The clinical significance of these findings is unclear but should warrant further analysis.

Nonspecific (Heterologous) Protection of Neonatal BCG Vaccination Against Hospitalization Due to Respiratory Infection and Sepsis

BACKGROUND

Bacille Calmette-Guerin (BCG) vaccination has been suggested to have nonspecific beneficial effects in children from developing countries, reducing morbidity and mortality caused by unrelated pathogens.

OBJECTIVE

We aimed to assess the heterologous protective effects of BCG vaccination against respiratory infection (RI) and sepsis not attributable to tuberculosis in children born in Spain.

METHODS

We conducted a retrospective epidemiological study using data from the Official Spanish Registry of Hospitalizations (CMBD-HA) to identify differences in hospitalization rates (HR) in BCG-vaccinated children (Basque Country, where neonatal BCG is part of the immunization schedule and has a 100% coverage) as compared to non-BCG-vaccinated children (from the rest of Spain, where BCG is not used).

RESULTS

A total of 464 611 hospitalization episodes from 1992 to 2011 were analyzed. The HR due to RI not attributable to tuberculosis in BCG-vaccinated children was significant lower compared to non-BCG-vaccinated children for all age groups, with a total preventive fraction (PF) of 41.4% (95% confidence interval: 40.3-42.5; P-value <.001). According to age group, PF was 32.4% (30.9-33.9; P-value <.001) for children under 1 year old, 60.1% (58.5-61.7; P-value <.001) for children between 1 and 4 years old, 66.6% (62.8-70.2; P-value <.001) for children between 5 and 9 years old, and 69.6% (63.3-75.0; P-value <.001) for children between 10 and 14 years old. The HR due to sepsis not attributable to tuberculosis in BCG-vaccinated children under 1 year of age was also significantly lower, with a PF of 52.8% (43.8-60.7; P-value <.001).

CONCLUSIONS

BCG vaccination at birth may decrease hospitalization due to RI and sepsis not related to tuberculosis through heterologous protection.

Evaluating the accuracy of AIM panels at quantifying genome ancestry

Background

There is a growing interest among geneticists in developing panels of Ancestry Informative Markers (AIMs) aimed at measuring the biogeographical ancestry of individual genomes. The efficiency of these panels is commonly tested empirically by contrasting self-reported ancestry with the ancestry estimated from these panels.

Results

Using SNP data from HapMap we carried out a simulation-based study aimed at measuring the effect of SNP coverage on the estimation of genome ancestry. For three of the main continental groups (Africans, East Asians, Europeans) ancestry was first estimated using the whole HapMap SNP database as a proxy for global genome ancestry; these estimates were subsequently compared to those obtained from pre-designed AIM panels. Panels that consider >400 AIMs capture genome ancestry reasonably well, while those containing a few dozen AIMs show a large variability in ancestry estimates. Curiously, 500-1,000 SNPs selected at random from the genome provide an unbiased estimate of genome ancestry and perform as well as any AIM panel of similar size. In simulated scenarios of population admixture, panels containing few AIMs also show important deficiencies to measure genome ancestry.

Conclusions

The results indicate that the ability to estimate genome ancestry is strongly dependent on the number of AIMs used, and not primarily on their individual informativeness. Caution should be taken when making individual (medical, forensic, or anthropological) inferences based on AIMs.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-543) contains supplementary material, which is available to authorized users.

A generalized model to estimate the statistical power in mitochondrial disease studies involving 2×k tables

Background

Mitochondrial DNA (mtDNA) variation (i.e. haplogroups) has been analyzed in regards to a number of multifactorial diseases. The statistical power of a case-control study determines the a priori probability to reject the null hypothesis of homogeneity between cases and controls.

Methods/Principal Findings

We critically review previous approaches to the estimation of the statistical power based on the restricted scenario where the number of cases equals the number of controls, and propose a methodology that broadens procedures to more general situations. We developed statistical procedures that consider different disease scenarios, variable sample sizes in cases and controls, and variable number of haplogroups and effect sizes. The results indicate that the statistical power of a particular study can improve substantially by increasing the number of controls with respect to cases. In the opposite direction, the power decreases substantially when testing a growing number of haplogroups. We developed mitPower (http://bioinformatics.cesga.es/mitpower/), a web-based interface that implements the new statistical procedures and allows for the computation of the a priori statistical power in variable scenarios of case-control study designs, or e.g. the number of controls needed to reach fixed effect sizes.

Conclusions/Significance

The present study provides with statistical procedures for the computation of statistical power in common as well as complex case-control study designs involving 2×k tables, with special application (but not exclusive) to mtDNA studies. In order to reach a wide range of researchers, we also provide a friendly web-based tool – mitPower – that can be used in both retrospective and prospective case-control disease studies.

Ancestry analysis reveals a predominant Native American component with moderate European admixture in Bolivians

We have genotyped 46 Ancestry Informative Markers (AIMs) in two of the most populated areas in Bolivia, namely, La Paz (Andean region; n=105), and Chuquisaca (Sub-Andean region; n=73). Using different analytical tools, we inferred admixture proportions of these two American communities by comparing the genetic profiles with those publicly available from the CEPH (Centre d'Etude du Polymorphisme Humain) panel representing three main continental groups (Africa, Europe, and America). By way of simulations, we first evaluated the minimum sample size needed in order to obtain accurate estimates of ancestry proportions. The results indicated that sample sizes above 30 individuals could be large enough to estimate main continental ancestry proportions using the 46 AIMs panel. With the exception of a few individuals, the results also indicated that Bolivians showed a predominantly Native American ancestry with variable levels of European admixture. The proportions of ancestry were statistically different in La Paz and Chuquisaca: the Native American component was 86% and 77% (Mann-Whitney U-test: un-adjusted P-value=2.1×10(-5)), while the European ancestry was 13% and 21% (Mann-Whitney U-test: un-adjusted P-value=3.6×10(-5)), respectively. The African ancestry in Bolivians captured by the AIMs analyzed in the present study was below 2%. The inferred ancestry of Bolivians fits well with previous studies undertaken on haplotype data, indicating a major proportion of Native American lineages. The genetic differences observed in these two groups suggest that forensic genetic analysis should be better performed based on local databases built in the main Bolivian areas.

The genetic legacy of the pre-colonial period in contemporary Bolivians

Only a few genetic studies have been carried out to date in Bolivia. However, some of the most important (pre)historical enclaves of South America were located in these territories. Thus, the (sub)-Andean region of Bolivia was part of the Inca Empire, the largest state in Pre-Columbian America. We have genotyped the first hypervariable region (HVS-I) of 720 samples representing the main regions in Bolivia, and these data have been analyzed in the context of other pan-American samples (>19,000 HVS-I mtDNAs). Entire mtDNA genome sequencing was also undertaken on selected Native American lineages. Additionally, a panel of 46 Ancestry Informative Markers (AIMs) was genotyped in a sub-set of samples. The vast majority of the Bolivian mtDNAs (98.4%) were found to belong to the main Native American haplogroups (A: 14.3%, B: 52.6%, C: 21.9%, D: 9.6%), with little indication of sub-Saharan and/or European lineages; however, marked patterns of haplogroup frequencies between main regions exist (e.g. haplogroup B: Andean [71%], Sub-Andean [61%], Llanos [32%]). Analysis of entire genomes unraveled the phylogenetic characteristics of three Native haplogroups: the pan-American haplogroup B2b (originated ∼21.4 thousand years ago [kya]), A2ah (∼5.2 kya), and B2o (∼2.6 kya). The data suggest that B2b could have arisen in North California (an origin even in the north most region of the American continent cannot be disregarded), moved southward following the Pacific coastline and crossed Meso-America. Then, it most likely spread into South America following two routes: the Pacific path towards Peru and Bolivia (arriving here at about ∼15.2 kya), and the Amazonian route of Venezuela and Brazil southwards. In contrast to the mtDNA, Ancestry Informative Markers (AIMs) reveal a higher (although geographically variable) European introgression in Bolivians (25%). Bolivia shows a decreasing autosomal molecular diversity pattern along the longitudinal axis, from the Altiplano to the lowlands. Both autosomes and mtDNA revealed a low impact (1-2%) of a sub-Saharan component in Bolivians.