A new method to reconstruct recombination events at a genomic scale

Recombination is one of the main forces shaping genome diversity, but the information it generates is often overlooked. A recombination event creates a junction between two parental sequences that may be transmitted to the subsequent generations. Just like mutations, these junctions carry evidence of the shared past of the sequences. We present the IRiS algorithm, which detects past recombination events from extant sequences and specifies the place of each recombination and which are the recombinants sequences. We have validated and calibrated IRiS for the human genome using coalescent simulations replicating standard human demographic history and a variable recombination rate model, and we have fine-tuned IRiS parameters to simultaneously optimize for false discovery rate, sensitivity, and accuracy in placing the recombination events in the sequence. Newer recombinations overwrite traces of past ones and our results indicate more recent recombinations are detected by IRiS with greater sensitivity. IRiS analysis of the MS32 region, previously studied using sperm typing, showed good concordance with estimated recombination rates. We also applied IRiS to haplotypes for 18 X-chromosome regions in HapMap Phase 3 populations. Recombination events detected for each individual were recoded as binary allelic states and combined into recotypes. Principal component analysis and multidimensional scaling based on recotypes reproduced the relationships between the eleven HapMap Phase III populations that can be expected from known human population history, thus further validating IRiS. We believe that our new method will contribute to the study of the distribution of recombination events across the genomes and, for the first time, it will allow the use of recombination as genetic marker to study human genetic variation.     

Human F7 sequence is split into three deep clades that are related to FVII plasma levels

It is widely accepted that FVII levels are strongly, consistently, and independently related to cardiovascular risk. These levels are influenced by genetic and environmental factors. Among the genetic factors, only a limited number of polymorphisms in the F7 gene have been reported, and they explain only a small proportion of the genetic variability. Recently, we have accomplished the complete dissection of the F7 quantitative trait locus responsible for all of the genetic variability observed in FVII levels. Now, we present the thorough study of the haplotype organization of F7 DNA sequence variation among individuals and the evolutionary processes that produced this variation, by sequencing 15 kb of genomic DNA sequence from the F7 locus in 40 unrelated individual (80 chromosomes) from the genetic analysis of idiopathic thrombophilia (GAIT) project as well as four non-human primate species. Our study revealed 49 polymorphisms, of which 39 SNPs were further considered. Genotyping of these DNA variations in the whole family-based GAIT sample helped resolve linkage phases, and a total of 37 distinct haplotypes were identified. Tajima’s D was significantly positive in this sample, suggesting balancing selection. This parameter was a reflection of the phylogenetic structure of F7 haplotype, which was deeply split into three well-supported clades or haplogroups, suggesting that functional differences among F7 variants do not depend on a few single-site variations. Moreover, haplogroup 2 was associated with high FVII levels and haplogroup 3 with low levels. In this study, we have for the first time established a clear relation between genotypic variability structure and phenotypic variability of a particular quantitative trait involved in a complex disease. Electronic Supplementary Material Supplementary material is available for this article at     

Evolutionary dynamics of the human <Emphasis Type="Italic">ABO</Emphasis> gene

The ABO polymorphism has long been suspected to be under balancing selection. To explore this possibility, we analyzed two datasets: (1) a set of 94 23-Kb sequences in European- and African-Americans produced by the Seattle SNPs project, and (2) a set of 814 2-Kb sequences in O alleles from seven worldwide populations. A phylogenetic analysis of the Seattle sequences showed a complex pattern in which the action of recombination and gene conversion are evident, and in which four main lineages could be individuated. The sequence patterns could be linked to the expected blood group phenotype; in particular, the main mutation giving rise to the null O allele is likely to have appeared at least three times in human evolution, giving rise to allele lineages O02, O01, and O09. However, the genealogy changes along the gene and variations of both numbers of branches and of their time depth were observed, which could result from a combined action of recombination and selection. Several neutrality tests clearly demonstrated deviations compatible with balancing selection, peaking at several locations along the gene. The time depth of the genealogy was also incompatible with neutral evolution, particularly in the region from exons 6 to 7, which codes for most of the catalytic domain. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

FABSIM: a software for generating FST distributions with various ascertainment biases

We have developed a software that applies ascertainment bias on simulated DNA sequences and calculates F(ST) on them, so they can be used to generate neutral distributions that are appropriate to test whether the genetic differentiation of a particular gene between populations is compatible with neutral evolution, or, on the contrary, suggests local adaptation by natural selection. AVAILABILITY: FABSIM is available from http://www.snpator.com/public/downloads/aRamirez/FABSIM/.     

Gene flow in environmental Legionella pneumophila leads to genetic and pathogenic heterogeneity within a Legionnaires’ disease outbreak

Background

Legionnaires’ disease is a severe form of pneumonia caused by the environmental bacterium Legionella pneumophila. Outbreaks commonly affect people with known risk factors, but the genetic and pathogenic complexity of L. pneumophila within an outbreak is not well understood. Here, we investigate the etiology of the major Legionnaires’ disease outbreak that occurred in Edinburgh, UK, in 2012, by examining the evolutionary history, genome content, and virulence of L. pneumophila clinical isolates.

Results

Our high resolution genomic approach reveals that the outbreak was caused by multiple genetic subtypes of L. pneumophila, the majority of which had diversified from a single progenitor through mutation, recombination, and horizontal gene transfer within an environmental reservoir prior to release. In addition, we discover that some patients were infected with multiple L. pneumophila subtypes, a finding which can affect the certainty of source attribution. Importantly, variation in the complement of type IV secretion systems encoded by different genetic subtypes correlates with virulence in a Galleria mellonella model of infection, revealing variation in pathogenic potential among the outbreak source population of L. pneumophila.

Conclusions

Taken together, our study indicates previously cryptic levels of pathogen heterogeneity within a Legionnaires’ disease outbreak, a discovery that impacts on source attribution for future outbreak investigations. Furthermore, our data suggest that in addition to host immune status, pathogen diversity may be an important influence on the clinical outcome of individual outbreak infections.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0504-1) contains supplementary material, which is available to authorized users.     

Repeated bouts of eccentrically biased endurance exercise stimulate salivary IgA secretion rate

To determine the salivary secretory immunoglobulin A (sIgA) response to repeated bouts of unaccustomed, downhill running (eccentrically biased) and examine potential protective immunological adaption from a repeated bout effect. Eleven active but untrained males (age: 19.7±0.4 years; VO_2peak: 47.8± 3.6 ml · kg⁻¹ · min ⁻¹) performed two 60 min bouts (Run 1 and Run 2) of downhill running (−13.5% gradient), separated by 14 days, at a speed eliciting 75% of their VO_2peak on a level grade. Saliva samples were collected before (baseline), immediately post exercise (IPE), and every hour for 12 h and every 24 h for 6 days after each run. Salivary sIgA concentration was measured and sIgA secretion rate was calculated. Results were analysed using repeated measures ANOVA (12 h period: 2x14; 24 h intervals: 2x7; p ≤ 0.05) with Tukey post-hoc tests where appropriate. Results are reported as means ± SE. There was a significant (p < 0.0001) interaction effect for sIgA secretion rate, IPE, with higher values after Run 2, as well as a significant (p < 0.01) time effect with elevated levels IPE and between 24 h and 144 h. There was a run effect (p < 0.0001), with the sIgA secretion rate significantly higher after Run 2. Repeated bouts of unaccustomed, eccentrically biased exercise induced alterations in the salivary sIgA secretion rate. This may serve as a protective mucosal adaptation to exercise-induced tissue damage.     

Nicaraguan population data on LDLR,GYPA, D7S8, HBGG,GC and HLA-DQA1 loci

Nicaraguans have become the most numerous and fastest increasing minority in Costa Rica: at present they represent around 6% of the total population of the country. We have analyzed the allele and genotype frequencies of six PCR-based genetic markers (LDLR, GYPA, HBGG, D7S8, GC, and HLA-DQA1) in 100 unrelated Nicaraguans living in Costa Rica. All loci studied were in Hardy-Weinberg equilibrium. Some statistical parameters of forensic interest were also calculated (h, PD and CE). Allele frequencies of the markers HLA-DQA1 and GYPA were found to be significantly different between the populations of Nicaragua and Costa Rica. Nevertheless, genetic distances showed that Nicaragua is close to other Hispanic-admixed populations like those from Argentina, Chile, Colombia, Costa Rica, and USA Hispanics. The loci set was assessed to be useful for paternity testing and individual identification in the Nicaraguan population residing in Costa Rica.     

Georgian and kurd mtDNA sequence analysis shows a lack of correlation between languages and female genetic lineages

Mitochondrial DNA sequences from Georgians and Kurds were analyzed in order to test the possible correlation between female lineages and languages in these two neighboring West Eurasian groups. Mitochondrial sequence pools in both populations are very similar despite their different linguistic and prehistoric backgrounds. Both populations present mtDNA lineages that clearly belong to the European gene pool, as shown by 1) similar nucleotide and sequence diversities; 2) a large number of sequences shared with the rest of European samples; 3) nonsignificant genetic distances; and 4) classification of the present lineages into the major European mtDNA haplogroups already described. The outlier position of the populations from the Caucasus according to classical genetic markers is not recognized in the present Georgian mtDNA sequence pool. This result suggests that the differentiation of mtDNA sequences in West Eurasia and the outlier features of Caucasian populations should be attributed to different processes. Moreover, the putative linguistic relationship between Caucasian groups and the Basques, another outlier population within Europe for classical genetic markers, is not detected by the analysis of mtDNA sequences.