Data size sufficiency analyses of haplotype inference algortihms

We present experimental and theoretical analyses of data requirements for haplotype inference algorithms. Our experiments include a broad range of problem sizes under two standard models of tree distribution and were designed to yield statistically robust results despite the size of the sample space. Our results validate Gusfield's conjecture that a population size of n log n is required to give (with high probability) sufficient information to deduce the n haplotypes and their complete evolutionary history. The experimental results inspired our experimental finding with theoretical bounds on the population size. We also analyze the population size required to deduce some fixed fraction of the evolutionary history of a set of n haplotypes and establish linear bounds on the required sample size. These linear bounds are also shown theoretically.

Multilocus sequence analysis reveals the genetic diversity of European fruit tree phytoplasmas and supports the existence of inter-species recombination

The genetic diversity of three temperate fruit tree phytoplasmas ‘Candidatus Phytoplasma prunorum’, ‘Ca. P. mali’ and ‘Ca. P. pyri’ has been established by multilocus sequence analysis. Among the four genetic loci used, the genes imp and aceF distinguished 30 and 24 genotypes, respectively, and showed the highest variability. Percentage of substitution for imp ranged from 50 to 68% according to species. Percentage of substitution varied between 9 and 12% for aceF, whereas it was between 5 and 6% for pnp and secY. In the case of ‘Ca P. prunorum’ the three most prevalent aceF genotypes were detected in both plants and insect vectors, confirming that the prevalent isolates are propagated by insects. The four isolates known to be hypo-virulent had the same aceF sequence, indicating a possible monophyletic origin. Haplotype network reconstructed by eBURST revealed that among the 34 haplotypes of ‘Ca. P. prunorum’, the four hypo-virulent isolates also grouped together in the same clade. Genotyping of some Spanish and Azerbaijanese ‘Ca. P. pyri’ isolates showed that they shared some alleles with ‘Ca. P. prunorum’, supporting for the first time to our knowledge, the existence of inter-species recombination between these two species.

On the use of principal components in contemporany population genetics: a case study

In human Population Genetics, routine applications of principal component techniques are oftenrequired. Population biologists make widespread use of certain discrete classifications of humansamples into haplotypes, the monophyletic units of phylogenetic trees constructed from severalsingle nucleotide bimorphisms hierarchically ordered. Compositional frequencies of the haplotypesare recorded within the different samples. Principal component techniques are then required as adimension-reducing strategy to bring the dimension of the problem to a manageable level, say two,to allow for graphical analysis.Population biologists at large are not aware of the special features of compositional data and normally make use of the crude covariance of compositional relative frequencies to construct principalcomponents. In this short note we present our experience with using traditional linear principalcomponents or compositional principal components based on logratios, with reference to a specificdataset

Minimizing recombinations in consensus networks for phylogeographic studies

Background: We address the problem of studying recombinational variations in (human) populations. In this paper, our focus is on one computational aspect of the general task: Given two networks G1 and G2, with both mutation and recombination events, defined on overlapping sets of extant units the objective is to compute a consensus network G3 with minimum number of additional recombinations. We describe a polynomial time algorithm with a guarantee that the number of computed new recombination events is within ϵ = sz(G1, G2) (function sz is a well-behaved function of the sizes and topologies of G1 and G2) of the optimal number of recombinations. To date, this is the best known result for a network consensus problem.Results: Although the network consensus problem can be applied to a variety of domains, here we focus on structure of human populations. With our preliminary analysis on a segment of the human Chromosome X data we are able to infer ancient recombinations, population-specific recombinations and more, which also support the widely accepted 'Out of Africa' model. These results have been verified independently using traditional manual procedures. To the best of our knowledge, this is the first recombinations-based characterization of human populations. Conclusion: We show that our mathematical model identifies recombination spots in the individual haplotypes; the aggregate of these spots over a set of haplotypes defines a recombinational landscape that has enough signal to detect continental as well as population divide based on a short segment of Chromosome X. In particular, we are able to infer ancient recombinations, population-specific recombinations and more, which also support the widely accepted 'Out of Africa' model. The agreement with mutation-based analysis can be viewed as an indirect validation of our results and the model. Since the model in principle gives us more information embedded in the networks, in our future work, we plan to investigate more non-traditional questions via these structures computed by our methodology.

African signatures of recent positive selection in human FOXI1

Background: The human FOXI1 gene codes for a transcription factor involved in the physiology of the inner ear, testis, and kidney. Using three interspecies comparisons, it has been suggested that this may be a gene underhuman-specific selection. We sought to confirm this finding by using an extended set of orthologous sequences.Additionally, we explored for signals of natural selection within humans by sequencing the gene in 20 Europeans,20 East Asians and 20 Yorubas and by analysing SNP variation in a 2 Mb region centered on FOXI1 in 39worldwide human populations from the HGDP-CEPH diversity panel.Results: The genome sequences recently available from other primate and non-primate species showed that FOXI1divergence patterns are compatible with neutral evolution. Sequence-based neutrality tests were not significant inEuropeans, East Asians or Yorubas. However, the Long Range Haplotype (LRH) test, as well as the iHS and XP-Rsbstatistics revealed significantly extended tracks of homozygosity around FOXI1 in Africa, suggesting a recentepisode of positive selection acting on this gene. A functionally relevant SNP, as well as several SNPs either on theputatively selected core haplotypes or with significant iHS or XP-Rsb values, displayed allele frequencies stronglycorrelated with the absolute geographical latitude of the populations sampled.Conclusions: We present evidence for recent positive selection in the FOXI1 gene region in Africa. Climate mightbe related to this recent adaptive event in humans. Of the multiple functions of FOXI1, its role in kidney-mediatedwater-electrolyte homeostasis is the most obvious candidate for explaining a climate-related adaptation.

Genetic evidence of heterogeneity in intrahepatic cholestasis of pregnancy

Background and aims: The aim of this study was to investigate the genetic aetiology of intrahepatic cholestasis of pregnancy (ICP) and the impact of known cholestasis genes (BSEP, FIC1, and MDR3) on the development of this disease. Patients and methods: Sixty nine Finnish ICP patients were prospectively interviewed for a family history of ICP, and clinical features were compared in patients with familial ICP (patients with a positive family history, n=11) and sporadic patients (patients with no known family history of ICP, n=58). For molecular genetic analysis, 16 individuals from two independently ascertained Finnish ICP families were genotyped for the flanking markers for BSEP, FIC1, and MDR3. Results: The pedigree structures in 16% (11/69) of patients suggested dominant inheritance. Patients with familial ICP had higher serum aminotransferase levels and a higher recurrence risk (92% v 40%). Both segregation of haplotypes and multipoint linkage analysis excluded BSEP, FIC1, and MDR3 genes in the studied pedigrees. Additionally, the MDR3 gene, previously shown to harbour mutations in ICP patients, was negative for mutations when sequenced in four affected individuals from the two families. Conclusions: These results support the hypothesis that the aetiology of ICP is heterogeneous and that ICP is due to a genetic predisposition in a proportion of patients. The results of molecular genetic analysis further suggest that the previously identified three cholestasis genes are not likely to be implicated in these Finnish ICP families with dominant inheritance.

Nitric oxide synthase (NOS) inhibition for one week improves renal sodium and water excretion in cirrhotic rats with ascites

Normalization of the increased vascular nitric oxide (NO) generation with low doses of NG-nitro-L-arginine methyl ester (L-NAME) corrects the hemodynamic abnormalities of cirrhotic rats with ascites. We have undertaken this study to investigate the effect of the normalization of vascular NO production, as estimated by aortic cyclic guanosine monophosphate (cGMP) concentration and endothelial nitric oxide synthase (eNOS) protein expression in the aorta and mesenteric artery, on sodium and water excretion. Rats with carbon tetrachloride-induced cirrhosis and ascites were investigated using balance studies. The cirrhotic rats were separated into two groups, one receiving 0.5 mg/kg per day of L-NAME (CIR-NAME) during 7 d, whereas the other group (CIR) was administrated the same volume of vehicle. Two other groups of rats were used as controls, one group treated with L-NAME and another group receiving the same volume of vehicle. Sodium and water excretion was measured on days 0 and 7. On day 8, blood samples were collected for electrolyte and hormone measurements, and aorta and mesenteric arteries were harvested for cGMP determination and nitric oxide synthase (NOS) immunoblotting. Aortic cGMP and eNOS protein expression in the aorta and mesenteric artery were increased in CIR as compared with CIR-NAME. Both cirrhotic groups had a similar decrease in sodium excretion on day 0 (0.7 versus 0.6 mmol per day, NS) and a positive sodium balance (+0.9 versus +1.2 mmol per day, NS). On day 7, CIR-NAME rats had an increase in sodium excretion as compared with the CIR rats (sodium excretion: 2.4 versus 0.7 mmol per day, P < 0.001) and a negative sodium balance (-0.5 versus +0.8 mmol per day, P < 0.001). The excretion of a water load was also increased after L-NAME administration (from 28+/-5% to 65+/-7, P < 0.05). Plasma renin activity, aldosterone and arginine vasopressin were also significantly decreased in the CIR-NAME, as compared with the CIR rats. The results thus indicate that normalization of aortic cGMP and eNOS protein expression in vascular tissue is associated with increased sodium and water excretion in cirrhotic rats with ascites.

Nitric oxide and the release of lipoprotein lipase from white adipose tissue

Background/Aim: Lipoprotein lipase (LPL) is the main enzyme responsible for the distribution of circulating triacylglycerides in tissues. Its regulation via release from active sites in the vascular endothelium is poorly understood. In a previous study we reported that in response to acute immobilization (IMMO), LPL activity rapidly increases in plasma and decreases in white adipose tissue (WAT) in rats. In other stress situations IMMO triggers a generalized increase in nitric oxide (NO) production. Methods/Results: Here we demonstrate that in rats: 1) in vivo acute IMMO rapidly increases NO concentrations in plasma 2) during acute IMMO the WAT probably produces NO via the endothelial isoform of nitric oxide synthase (eNOS) from vessels, and 3) epididymal WAT perfused in situ with an NO donor rapidly releases LPL from the endothelium. Conclusion: We propose the following chain of events: stress stimulus / rapid increase of NO production in WAT (by eNOS) / release of LPL from the endothelium in WAT vessels. This chain of events could be a new mechanism that promotes the rapid decrease of WAT LPL activity in response to a physiological stimulus.

Dugesia sicula (Platyhelminthes, Tricladida): the colonizing success of an asexual Planarian

Background Dugesia sicula is the only species of its genus not presenting an endemic or restricted distribution within the Mediterranean area. It mostly comprises fissiparous populations (asexual reproduction by body division and regeneration), most likely sexually sterile, and characterized by an extremely low genetic diversity interpreted as the consequence of a recent anthropic expansion. However, its fissiparous reproduction can result in an apparent lack of diversity within the species, since genetic variation within individuals can be as large as between them because most individuals within a population are clones. We have estimated haplotype and nucleotide diversity of cytochrome oxidase I within and among individuals along the species distribution of a broad sample of D. sicula, including asexual and the two only sexual populations known today; and predicted its potential distribution based on climatic variables. Our aim was to determine the centre of colonisation origin, whether the populations are recent, and whether the species is expanding. Results The species presents 3 most frequent haplotypes, differing in a maximum of 11 base pairs. As expected from their fissiparous mode of reproduction, in half of all the analysed localities many individuals have multiple heteroplasmic haplotypes. The distribution of haplotypes is not geographically structured; however, the distribution of haplotypes and heteroplasmic populations shows higher diversity in the central Mediterranean region. The potential distribution predicted by climatic variables based modelling shows a preference for coastal areas and fits well with the observed data. Conclusions The distribution and frequency of the most frequent haplotypes and the presence of heteroplasmic individuals allow us to gain an understanding of the recent history of the species, together with previous knowledge on its phylogenetic relationships and age: The species most probably originated in Africa and dispersed through the central Mediterranean. After one or multiple populations became triploid and fissiparous, the species colonized the Mediterranean basin, likely both by its own means and helped by human activities. Its present distribution practically fulfils its potential distribution as modelled with climatic variables. Its prevalence in coastal regions with higher water temperatures predicts a likely future expansion to northern and more interior areas following the increase in temperatures due to climate change.

Polymorphisms in DNA-repair genes in a cohort of prostate cancer patients from different areas in Spain: heterogeneity between populations as a confounding factor in association studies

Background: Differences in the distribution of genotypes between individuals of the same ethnicity are an important confounder factor commonly undervalued in typical association studies conducted in radiogenomics. Objective: To evaluate the genotypic distribution of SNPs in a wide set of Spanish prostate cancer patients for determine the homogeneity of the population and to disclose potential bias. Design, Setting, and Participants: A total of 601 prostate cancer patients from Andalusia, Basque Country, Canary and Catalonia were genotyped for 10 SNPs located in 6 different genes associated to DNA repair: XRCC1 (rs25487, rs25489, rs1799782), ERCC2 (rs13181), ERCC1 (rs11615), LIG4 (rs1805388, rs1805386), ATM (rs17503908, rs1800057) and P53 (rs1042522). The SNP genotyping was made in a Biotrove OpenArrayH NT Cycler. Outcome Measurements and Statistical Analysis: Comparisons of genotypic and allelic frequencies among populations, as well as haplotype analyses were determined using the web-based environment SNPator. Principal component analysis was made using the SnpMatrix and XSnpMatrix classes and methods implemented as an R package. Non-supervised hierarchical cluster of SNP was made using MultiExperiment Viewer. Results and Limitations: We observed that genotype distribution of 4 out 10 SNPs was statistically different among the studied populations, showing the greatest differences between Andalusia and Catalonia. These observations were confirmed in cluster analysis, principal component analysis and in the differential distribution of haplotypes among the populations. Because tumor characteristics have not been taken into account, it is possible that some polymorphisms may influence tumor characteristics in the same way that it may pose a risk factor for other disease characteristics. Conclusion: Differences in distribution of genotypes within different populations of the same ethnicity could be an important confounding factor responsible for the lack of validation of SNPs associated with radiation-induced toxicity, especially when extensive meta-analysis with subjects from different countries are carried out.

Interrogating eleven fast-evolving genes for signatures of recent positive selection in worldwide human populations

Different signatures of natural selection persist over varying time scales in our genome, revealing possible episodes of adaptative evolution during human history. Here, we identify genes showing signatures of ancestral positive selection in the human lineage and investigate whether some of those genes have been evolving adaptatively in extant human populations. Specifically, we compared more than 11,000 human genes with their orthologs inchimpanzee, mouse, rat and dog and applied a branch-site likelihood method to test for positive selection on the human lineage. Among the significant cases, a robust set of 11 genes were then further explored for signatures of recent positive selection using SNP data. We genotyped 223 SNPs in 39 worldwide populations from the HGDP Diversity panel and supplemented this information with available genotypes for up to 4,814 SNPs distributed along 2 Mb centered on each gene. After exploring the allele frequency spectrum, population differentiation and the maintainance of long unbroken haplotypes, we found signals of recent adaptative phenomena in only one of the 11 candidate gene regions. However, the signal ofrecent selection in this region may come from a different, neighbouring gene (CD5) ratherthan from the candidate gene itself (VPS37C). For this set of positively-selected genes in thehuman lineage, we find no indication that these genes maintained their rapid evolutionarypace among human populations. Based on these data, it therefore appears that adaptation forhuman-specific and for population-specific traits may have involved different genes.

CX3CR1 is a modifying gene of survival and progression in Amyotrophic Lateral Sclerosis

The objective of this study was to investigate the association of functional variants of the human CX3CR1 gene (Fractalkine receptor) with the risk of Amyotrophic Lateral Sclerosis disease (ALS), the survival and the progression rate of the disease symptoms in a Spanish ALS cohort. 187 ALS patients (142sporadic [sALS] and 45 familial) and 378 controls were recruited. We investigated CX3CR1 V249I (rs3732379) and T280M (rs3732378) genotypes and their haplotypes as predictors of survival, the progression rate of the symptoms (as measured by ALSFRS-R and FVC decline) and the risk of suffering ALS disease. The results indicated that sALS patients with CX3CR1 249I/I or 249V/I genotypes presented a shorter survival time (42.21±4.82) than patients with 249V/V genotype (67.48±7.28; diff=-25.27 months 95%CI [-42.09,-8.45]; p=0.003). The survival time was shorter in sALS patients with spinal topography and CX3CR1 249I alleles (diff -29.78 months; 95%CI [-49.42,-10.14]; p=0.003). The same effects were also observed in the spinal sALS patients with 249I-280M haplotype (diff=-27.23 months; 95%CI [-49.83,-4.64]; p=0.018). In the sALS group, the CX3CR1 249I variant was associated with a faster progression of the disease symptoms (OR= 2.32; 95IC%[1.24, 4.33]; p=0.007). There was no evidence for association of these two CX3CR1 variants with ALS disease risk. The association evidenced herein is clinically relevant and indicates that CX3CR1 could be a disease-modifying gene in sALS. The progression rate of the disease"s symptoms and the survival time is affected in patients with one or two copies of the CX3CR1 249I allele. The CX3CR1 is the most potent ALS survival genetic factor reported to date. These results reinforce the role of the immune system in ALS pathogenesis.

Association of MMP-2 polymorphisms with severe and very severe COPD: a case control study of MMPs-1, 9 and 12 in a European population.

BACKGROUND: Genetic factors play a role in chronic obstructive pulmonary disease (COPD) but are poorly understood. A number of candidate genes have been proposed on the basis of the pathogenesis of COPD. These include the matrix metalloproteinase (MMP) genes which play a role in tissue remodelling and fit in with the protease--antiprotease imbalance theory for the cause of COPD. Previous genetic studies of MMPs in COPD have had inadequate coverage of the genes, and have reported conflicting associations of both single nucleotide polymorphisms (SNPs) and SNP haplotypes, plausibly due to under-powered studies. METHODS: To address these issues we genotyped 26 SNPs, providing comprehensive coverage of reported SNP variation, in MMPs- 1, 9 and 12 from 977 COPD patients and 876 non-diseased smokers of European descent and evaluated their association with disease singly and in haplotype combinations. We used logistic regression to adjust for age, gender, centre and smoking history. RESULTS: Haplotypes of two SNPs in MMP-12 (rs652438 and rs2276109), showed an association with severe/very severe disease, corresponding to GOLD Stages III and IV. CONCLUSIONS: Those with the common A-A haplotype for these two SNPs were at greater risk of developing severe/very severe disease (p = 0.0039) while possession of the minor G variants at either SNP locus had a protective effect (adjusted odds ratio of 0.76; 95% CI 0.61 - 0.94). The A-A haplotype was also associated with significantly lower predicted FEV1 (42.62% versus 44.79%; p = 0.0129). This implicates haplotypes of MMP-12 as modifiers of disease severity.

Natural or naturalized? phylogeography suggests that the abundant sea urchin arbacia lixula Is a recent Colonizer of the Mediterranean

We present the global phylogeography of the black sea urchin Arbacia lixula, an amphi-Atlantic echinoid with potential to strongly impact shallow rocky ecosystems. Sequences of the mitochondrial cytochrome c oxidase gene of 604 specimens from 24 localities were obtained, covering most of the distribution area of the species, including the Mediterranean and both shores of the Atlantic. Genetic diversity measures, phylogeographic patterns, demographic parameters and population differentiation were analysed. We found high haplotype diversity but relatively low nucleotide diversity, with 176 haplotypes grouped within three haplogroups: one is shared between Eastern Atlantic (including Mediterranean) and Brazilian populations, the second is found in Eastern Atlantic and the Mediterranean and the third is exclusively from Brazil. Significant genetic differentiation was found between Brazilian, Eastern Atlantic and Mediterranean regions, but no differentiation was found among Mediterranean sub-basins or among Eastern Atlantic sub-regions. The star-shaped topology of the haplotype network and the unimodal mismatch distributions of Mediterranean and Eastern Atlantic samples suggest that these populations have suffered very recent demographic expansions. These expansions could be dated 94-205 kya in the Mediterranean, and 31-67 kya in the Eastern Atlantic. In contrast, Brazilian populations did not show any signature of population expansion. Our results indicate that all populations of A. lixula constitute a single species. The Brazilian populations probably diverged from an Eastern Atlantic stock. The present-day genetic structure of the species in Eastern Atlantic and the Mediterranean is shaped by very recent demographic processes. Our results support the view (backed by the lack of fossil record) that A. lixula is a recent thermophilous colonizer which spread throughout the Mediterranean during a warm period of the Pleistocene, probably during the last interglacial. Implications for the possible future impact of A. lixula on shallow Mediterranean ecosystems in the context of global warming trends must be considered.

Time and space: genetic structure of the common sea urchin paracentrotus lividus in the south-eastern iberian coast

Spatio-temporal variability in settlement and recruitment, high mortality during the first life-history stages, and selection may determine the genetic structure of cohorts of long-lived marine invertebrates at small scales. We conducted a spatial and temporal analysis of the common Mediterranean Sea urchin Paracentrotus lividus to determine the genetic structure of cohorts at different scales. In Tossa de Mar (NW Mediterranean), recruitment was followed over 5 consecutive springs (2006-2010). In spring 2008, recruits and two-year-old individuals were collected at 6 locations along East and South Iberian coasts separated from 200 to over 1,100 km. All cohorts presented a high genetic diversity based on a fragment of mtCOI. Our results showed a marked genetic homogeneity in the temporal monitoring and a low degree of spatial structure in 2006. In 2008, coupled with an abnormality in the usual circulation patterns in the area, the genetic structure of the southern populations studied changed markedly, with arrival of many private haplotypes. This fact highlights the importance of point events in renewing the genetic makeup of populations, which can only be detected through analysis of the cohort structure coupling temporal and spatial perspectives.