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.

Nucleotide, cytogenetic and expression impact of the human chromosome 8p23.1 inversion polymorphism

Background: The human chromosome 8p23.1 region contains a 3.8–4.5 Mb segment which can be found in different orientations (defined as genomic inversion) among individuals. The identification of single nucleotide polymorphisms (SNPs) tightly linked to the genomic orientation of a given region should be useful to indirectly evaluate the genotypes of large genomic orientations in the individuals. Results: We have identified 16 SNPs, which are in linkage disequilibrium (LD) with the 8p23.1 inversion as detected by fluorescent in situ hybridization (FISH). The variability of the 8p23.1 orientation in 150 HapMap samples was predicted using this set of SNPs and was verified by FISH in a subset of samples. Four genes (NEIL2, MSRA, CTSB and BLK) were found differentially expressed (p<0.0005) according to the orientation of the 8p23.1 region. Finally, we have found variable levels of mosaicism for the orientation of the 8p23.1 as determined by FISH. Conclusion: By means of dense SNP genotyping of the region, haplotype-based computational analyses and FISH experiments we could infer and verify the orientation status of alleles in the 8p23.1 region by detecting two short haplotype stretches at both ends of the inverted region, which are likely the relic of the chromosome in which the original inversion occurred. Moreover, an impact of 8p23.1 inversion on gene expression levels cannot be ruled out, since four genes from this region have statistically significant different expression levels depending on the inversion status. FISH results in lymphoblastoid cell lines suggest the presence of mosaicism regarding the 8p23.1 inversion.

Altered brain-derived neurotrophic factor blood levels and gene variability are associated with anorexia and bulimia

Murine models and association studies in eating disorder (ED) patients have shown a role for the brain-derived neurotrophic factor (BDNF) in eating behavior. Some studies have shown association of BDNF -270C/T single-nucleotide polymorphism (SNP) with bulimia nervosa (BN), while BDNF Val66Met variant has been shown to be associated with both BN and anorexia nervosa (AN). To further test the role of this neurotrophin in humans, we screened 36 SNPs in the BDNF gene and tested for their association with ED and plasma BDNF levels as a quantitative trait. We performed a family-based association study in 106 ED nuclear families and analyzed BDNF blood levels in 110 ED patients and in 50 sib pairs discordant for ED. The rs7124442T/rs11030102C/rs11030119G haplotype was found associated with high BDNF levels (mean BDNF TCG haplotype carriers = 43.6 ng/ml vs. mean others 23.0 ng/ml, P = 0.016) and BN (Z = 2.64; P recessive = 0.008), and the rs7934165A/270T haplotype was associated with AN (Z =-2.64; P additive = 0.008). The comparison of BDNF levels in 50 ED discordant sib pairs showed elevated plasma BDNF levels for the ED group (mean controls = 41.0 vs. mean ED = 52.7; P = 0.004). Our data strongly suggest that altered BDNF levels modulated by BDNF gene variability are associated with the susceptibility to ED, providing physiological evidence that BDNF plays a role in the development of AN and BN, and strongly arguing for its involvement in eating behavior and body weight regulation.

A variant in the gene FUT9 is associated with susceptibility to placental malaria infection

Malaria in pregnancy forms a substantial part of the worldwide burden of malaria, with an estimated annual death toll of up to 200,000 infants, as well as increased maternal morbidity and mortality. Studies of genetic susceptibility to malaria have so far focused on infant malaria, with only a few studies investigating the genetic basis of placental malaria, focusing only on a limited number of candidate genes. The aim of this study therefore was to identify novel host genetic factors involved in placental malaria infection. To this end we carried out a nested case-control study on 180 Mozambican pregnant women with placental malaria infection, and 180 controls within an intervention trial of malaria prevention. We genotyped 880 SNPs in a set of 64 functionally related genes involved in glycosylation and innate immunity. A SNP located in the gene FUT9, rs3811070, was significantly associated with placental malaria infection (OR = 2.31, permutation p-value = 0.028). Haplotypic analysis revealed a similarly strong association of a common haplotype of four SNPs including rs3811070. FUT9 codes for a fucosyl-transferase that is catalyzing the last step in the biosynthesis of the Lewis-x antigen, which forms part of the Lewis blood group-related antigens. These results therefore suggest an involvement of this antigen in the pathogenesis of placental malaria infection.

Evolutionary analysis of genes of two pathways involved in placental malaria infection

Placental malaria is a special form of malaria that causes up to 200,000 maternal and infant deaths every year. Previous studies show that two receptor molecules, hyaluronic acid and chondroitin sulphate A, are mediating the adhesion of parasite-infected erythrocytes in the placenta of patients, which is believed to be a key step in the pathogenesis of the disease. In this study, we aimed at identifying sites of malaria-induced adaptation by scanning for signatures of natural selection in 24 genes in the complete biosynthesis pathway of these two receptor molecules. We analyzed a total of 24 Mb of publicly available polymorphism data from the International HapMap project for three human populations with European, Asian and African ancestry, with the African population from a region of presently and historically high malaria prevalence. Using the methods based on allele frequency distributions, genetic differentiation between populations, and on long-range haplotype structure, we found only limited evidence for malaria-induced genetic adaptation in this set of genes in the African population; however, we identified one candidate gene with clear evidence of selection in the Asian population. Although historical exposure to malaria in this population cannot be ruled out, we speculate that it might be caused by other pathogens, as there is growing evidence that these molecules are important receptors in a variety of host-pathogen interactions. We propose to use the present methods in a systematic way to help identify candidate regions under positive selection as a consequence of malaria.

Genetic characterization of the ABO blood group in Neandertals

Background: The high polymorphism rate in the human ABO blood group gene seems to be related to susceptibility to different pathogens. It has been estimated that all genetic variation underlying the human ABO alleles appeared along the human lineage, after the divergence from the chimpanzee lineage. A paleogenetic analysis of the ABO blood group gene in Neandertals allows us to directly test for the presence of the ABO alleles in these extinct humans. Results: We have analysed two male Neandertals that were retrieved under controlled conditions at the El Sidron site in Asturias (Spain) and that appeared to be almost free of modern human DNA contamination. We find a human specific diagnostic deletion for blood group O (O01 haplotype) in both Neandertal individuals. Conclusion: These results suggest that the genetic change responsible for the O blood group in humans predates the human and Neandertal divergence. A potential selective event associated with the emergence of the O allele may have therefore occurred after humans separated from their common ancestor with chimpanzees and before the human-Neandertal population divergence.

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.

Schmidtea mediterranea phylogeography: an old species surviving on a few mediterranean islands?

Schmidtea mediterranea (Platyhelminthes, Tricladida, Continenticola) is found in scattered localities on a few islands and in coastal areas of the western Mediterranean. Although S. mediterranea is the object of many regeneration studies, little is known about its evolutionary history. Its present distribution has been proposed to stem from the fragmentation and migration of the Corsica-Sardinia microplate during the formation of the western Mediterranean basin, which implies an ancient origin for the species. To test this hypothesis, we obtained a large number of samples from across its distribution area. Using known and new molecular markers and, for the first time in planarians, a molecular clock, we analysed the genetic variability and demographic parameters within the species and between its sexual and asexual populations to estimate when they diverged. Results: A total of 2 kb from three markers (COI, CYB and a nuclear intron N13) was amplified from ~200 specimens. Molecular data clustered the studied populations into three groups that correspond to the west, central and southeastern geographical locations of the current distribution of S. mediterranea. Mitochondrial genes show low haplotype and nucleotide diversity within populations but demonstrate higher values when all individuals are considered. The nuclear marker shows higher values of genetic diversity than the mitochondrial genes at the population level, but asexual populations present lower variability than the sexual ones. Neutrality tests are significant for some populations. Phylogenetic and dating analyses show the three groups to be monophyletic, with the west group being the basal group. The time when the diversification of the species occurred is between ~20 and ~4 mya, although the asexual nature of the western populations could have affected the dating analyses. Conclusions: S. mediterranea is an old species that is sparsely distributed in a harsh habitat, which is probably the consequence of the migration of the Corsica-Sardinia block. This species probably adapted to temperate climates in the middle of a changing Mediterranean climate that eventually became dry and hot. These data also suggest that in the mainland localities of Europe and Africa, sexual individuals of S. mediterranea are being replaced by asexual individuals that are either conspecific or are from other species that are better adapted to the Mediterranean climate.

Statistical properties of new neutrality tests against population growth.

A number of statistical tests for detecting population growth are described. We compared the statistical power of these tests with that of others available in the literature. The tests evaluated fall into three categories: those tests based on the distribution of the mutation frequencies, on the haplotype distribution, and on the mismatch distribution. We found that, for an extensive variety of cases, the most powerful tests for detecting population growth are Fu"s FS test and the newly developed R2 test. The behavior of the R2 test is superior for small sample sizes, whereas FS is better for large sample sizes. We also show that some popular statistics based on the mismatch distribution are very conservative. Key words: population growth, population expansion, coalescent simulations, neutrality tests

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.

DnaSP v5: A software for comprehensive analysis of DNA polymorphism data

DnaSP is a software package for a comprehensive analysis of DNA polymorphism data. Version 5 implements a number of new features and analytical methods allowing extensive DNA polymorphism analyses on large datasets. Among other features, the newly implemented methods allow for: (i) analyses on multiple data files; (ii) haplotype phasing; (iii) analyses on insertion/deletion polymorphism data; (iv) visualizing sliding window results integrated with available genome annotations in the UCSC browser.