Population subdivision and molecular sequence variation: Theory and analysis of Drosophila ananassae data

Population subdivision complicates analysis of molecular variation. Even if neutrality is assumed, three evolutionary forces need to be considered: migration, mutation, and drift. Simplification can be achieved by assuming that the process of migration among and drift within subpopulations is occurring fast compared to Mutation and drift in the entire population. This allows a two-step approach in the analysis: (i) analysis of population subdivision and (ii) analysis of molecular variation in the migrant pool. We model population subdivision using an infinite island model, where we allow the migration/drift parameter Theta to vary among populations. Thus, central and peripheral populations can be differentiated. For inference of Theta, we use a coalescence approach, implemented via a Markov chain Monte Carlo (MCMC) integration method that allows estimation of allele frequencies in the migrant pool. The second step of this approach (analysis of molecular variation in the migrant pool) uses the estimated allele frequencies in the migrant pool for the study of molecular variation. We apply this method to a Drosophila ananassae sequence data set. We find little indication of isolation by distance, but large differences in the migration parameter among populations. The population as a whole seems to be expanding. A population from Bogor (Java, Indonesia) shows the highest variation and seems closest to the species center.

PPARGC1A sequence variation and cardiovascular risk-factor levels: a study of the main genetic effects and gene x environment interactions in children from the European Youth Heart Study.

AIMS/HYPOTHESIS: The PPARGC1A gene coactivates multiple nuclear transcription factors involved in cellular energy metabolism and vascular stasis. In the present study, we genotyped 35 tagging polymorphisms to capture all common PPARGC1A nucleotide sequence variations and tested for association with metabolic and cardiovascular traits in 2,101 Danish and Estonian boys and girls from the European Youth Heart Study, a multicentre school-based cross-sectional cohort study. METHODS: Fasting plasma glucose concentrations, anthropometric variables and blood pressure were measured. Habitual physical activity and aerobic fitness were objectively assessed using uniaxial accelerometry and a maximal aerobic exercise stress test on a bicycle ergometer, respectively. RESULTS: In adjusted models, nominally significant associations were observed for BMI (rs10018239, p = 0.039), waist circumference (rs7656250, p = 0.012; rs8192678 [Gly482Ser], p = 0.015; rs3755863, p = 0.02; rs10018239, beta = -0.01 cm per minor allele copy, p = 0.043), systolic blood pressure (rs2970869, p = 0.018) and fasting glucose concentrations (rs11724368, p = 0.045). Stronger associations were observed for aerobic fitness (rs7656250, p = 0.005; rs13117172, p = 0.008) and fasting glucose concentrations (rs7657071, p = 0.002). None remained significant after correcting for the number of statistical comparisons. We proceeded by testing for gene x physical activity interactions for the polymorphisms that showed nominal evidence of association in the main effect models. None of these tests was statistically significant. CONCLUSIONS/INTERPRETATION: Variants at PPARGC1A may influence several metabolic traits in this European paediatric cohort. However, variation at PPARGC1A is unlikely to have a major impact on cardiovascular or metabolic health in these children.

Effects of a sex-ratio distorting endosymbiont on mtDNA variation in a global insect pest

Background: Patterns of mtDNA variation within a species reflect long-term population structure, but may also be influenced by maternally inherited endosymbionts, such as Wolbachia. These bacteria often alter host reproductive biology and can drive particular mtDNA haplotypes through populations. We investigated the impacts of Wolbachia infection and geography on mtDNA variation in the diamondback moth, a major global pest whose geographic distribution reflects both natural processes and transport via human agricultural activities. Results: The mtDNA phylogeny of 95 individuals sampled from 10 countries on four continents revealed two major clades. One contained only Wolbachia-infected individuals from Malaysia and Kenya, while the other contained only uninfected individuals, from all countries including Malaysia and Kenya. Within the uninfected group was a further clade containing all individuals from Australasia and displaying very limited sequence variation. In contrast, a biparental nuclear gene phylogeny did not have infected and uninfected clades, supporting the notion that maternally-inherited Wolbachia are responsible for the mtDNA pattern. Only about 5% (15/306) of our global sample of individuals was infected with the plutWBI isolate and even within infected local populations, many insects were uninfected. Comparisons of infected and uninfected isofemale lines revealed that plutWBI is associated with sex ratio distortion. Uninfected lines have a 1:1 sex ratio, while infected ones show a 2:1 female bias. Conclusion: The main correlate of mtDNA variation in P. xylostella is presence or absence of the plutWBI infection. This is associated with substantial sex ratio distortion and the underlying mechanisms deserve further study. In contrast, geographic origin is a poor predictor of moth mtDNA sequences, reflecting human activity in moving the insects around the globe. The exception is a clade of Australasian individuals, which may reflect a bottleneck during their recent introduction to this region.

DNA variation in a 13-Mb region including the F9 gene: inferring the genealogical history and causal role of a hemophilia B mutation (IVS 5+13 A -> G)

About 5.5% of all UK hemophilia B patients have the base substitution IVS 5+13 A-->G as the only change in their factor (F)IX gene (F9). This generates a novel donor splice site which fits the consensus better than the normal intron 5 donor splice. Use of the novel splice site should result in a missense mutation followed by the abnormal addition of four amino acids to the patients' FIX. In order to explain the prevalence of this mutation, its genealogical history is examined. Analysis of restriction fragment length polymorphism in the 21 reference UK individuals (from different families) with the above mutation showed identical haplotypes in 19 while two differed from the rest and from each other. In order to investigate the history of the mutation and to verify that it had occurred independently more than once, the sequence variation in 1.5-kb segments scattered over a 13-Mb region including F9 was examined in 18 patients and 15 controls. This variation was then analyzed with a recently developed Bayesian approach that reconstructs the genealogy of the gene investigated while providing evidence of independent mutations that contribute disconnected branches to the genealogical tree. The method also provides minimum estimates of the age of the mutation inherited by the members of coherent trees. This revealed that 17 or 18 mutant genes descend from a founder who probably lived 450 years ago, while one patient carries an independent mutation. The independent recurrence of the IVS5+13 A-->G mutation strongly supports the conclusion that it is the cause of these patients' mild hemophilia.

Selection of candidate coding DNA barcoding regions for use on land plants

An in silico screen of 41 of the 81 coding regions of the Nicotiana plastid genome generated a shortlist of 12 candidates as DNA barcoding loci for land plants. These loci were evaluated for amplification and sequence variation against a reference set of 98 land plant taxa. The deployment of multiple primers and a modified multiplexed tandem polymerase chain reaction yielded 85–94% amplification across taxa, and mean sequence differences between sister taxa of 6.1 from 156 bases of accD to 22 from 493 bases of matK. We conclude that loci should be combined for effective diagnosis, and recommend further investigation of the following six loci: matK, rpoB, rpoC1, ndhJ, ycf5 and accD.

Inbreeding depression and founder diversity among captive and free-living populations of the endangered pink pigeon Columba mayeri

The endemic pink pigeon has recovered from less than 20 birds in the mid-1970s to 355 free-living individuals in 2003. A major concern for the species' recovery has been the potential genetic problem of inbreeding. Captive pink pigeons bred for reintroduction were managed to maximise founder representation and minimise inbreeding. In this paper, we quantify the effect of inbreeding on survival and reproductive parameters in captive and wild populations and quantify DNA sequence variation in the mitochondrial d-loop region for pink pigeon founders. Inbreeding affected egg fertility, squab, juvenile and adult survival, but effects were strongest in highly inbred birds (F≥0.25). Inbreeding depression was more apparent in free-living birds where even moderate levels of inbreeding affected survival, although highly inbred birds were equally compromised in both captive and wild populations. Mitochondrial DNA haplotypic diversity in pink pigeon founders is low, suggesting that background inbreeding is contributing to low fertility and depressed productivity in this species, as well as comparable survival of some groups of non-inbred and nominally inbred birds. Management of wild populations has boosted population growth and may be required long-term to offset the negative effects of inbreeding depression and enhance the species' survival.

A functional genomics approach reveals novel quantitative trait loci associated with platelet signalling pathways

Platelet response to activation varies widely between individuals but shows interindividual consistency and strong heritability. The genetic basis of this variation has not been properly explored. We therefore systematically measured the effect on function of sequence variation in 97 candidate genes in the collagen and adenosine-diphosphate (ADP) signaling pathways. Resequencing of the genes in 48 European DNA samples nearly doubled the number of known single nucleotide polymorphisms (SNPs) and informed the selection of 1327 SNPs for genotyping in 500 healthy Northern European subjects with known platelet responses to collagen-related peptide (CRP-XL) and ADP. This identified 17 novel associations with platelet function (P < .005) accounting for approximately 46% of the variation in response. Further investigations with platelets of known genotype explored the mechanisms behind some of the associations. SNPs in PEAR1 associated with increased platelet response to CRP-XL and increased PEAR1 protein expression after platelet degranulation. The minor allele of a 3' untranslated region (UTR) SNP (rs2769668) in VAV3 was associated with higher protein expression (P = .03) and increased P-selectin exposure after ADP activation (P = .004). Furthermore the minor allele of the intronic SNP rs17786144 in ITPR1 modified Ca2+ levels after activation with ADP (P < .004). These data provide novel insights into key hubs within platelet signaling networks.

Friends or relatives? Phylogenetics and species delimitation in the controversial European orchid Genus ophrys

Background and Aims Highly variable, yet possibly convergent, morphology and lack of sequence variation have severely hindered production of a robust phylogenetic framework for the genus Ophrys. The aim of this study is to produce this framework as a basis for more rigorous species delimitation and conservation recommendations. Methods Nuclear and plastid DNA sequencing and amplified fragment length polymorphism (AFLP) were performed on 85 accessions of Ophrys, spanning the full range of species aggregates currently recognized. Data were analysed using a combination of parsimony and Bayesian tree-building techniques and by principal coordinates analysis. Key Results Complementary phylogenetic analyses and ordinations using nuclear, plastid and AFLP datasets identify ten genetically distinct groups (six robust) within the genus that may in turn be grouped into three sections (treated as subgenera by some authors). Additionally, genetic evidence is provided for a close relationship between the O. tenthredinifera, O. bombyliflora and O. speculum groups. The combination of these analytical techniques provides new insights into Ophrys systematics, notably recognition of the novel O. umbilicata group. Conclusions Heterogeneous copies of the nuclear ITS region show that some putative Ophrys species arose through hybridization rather than divergent speciation. The supposedly highly specific pseudocopulatory pollination syndrome of Ophrys is demonstrably 'leaky', suggesting that the genus has been substantially over-divided at the species level.

The evolution of mobile DNAs: when will transposons create phylogenies that look as if there is a master gene?

Some families of mammalian interspersed repetitive DNA, such as the Alu SINE sequence, appear to have evolved by the serial replacement of one active sequence with another, consistent with there being a single source of transposition: the "master gene." Alternative models, in which multiple source sequences are simultaneously active, have been called "transposon models." Transposon models differ in the proportion of elements that are active and in whether inactivation occurs at the moment of transposition or later. Here we examine the predictions of various types of transposon model regarding the patterns of sequence variation expected at an equilibrium between transposition, inactivation, and deletion. Under the master gene model, all bifurcations in the true tree of elements occur in a single lineage. We show that this property will also hold approximately for transposon models in which most elements are inactive and where at least some of the inactivation events occur after transposition. Such tree shapes are therefore not conclusive evidence for a single source of transposition.

The impact of genome defense on mobile elements in Microbotryum

Repeat induced point mutation (RIP), a mechanism causing hypermutation of repetitive DNA sequences in fungi, has been described as a ‘genome defense’ which functions to inactivate mobile elements and inhibit their deleterious effects on genome stability. Here we address the interactions between RIP and transposable elements in the Microbotryum violaceum species complex. Ten strains of M. violaceum, most of which belong to different species of the fungus, were all found to contain intragenomic populations of copia-like retrotransposons. Intragenomic DNA sequence variation among the copia-like elements was analyzed for evidence of RIP. Among species with RIP, there was no significant correlation between the frequency of RIP-induced mutations and inferred transposition rate based on diversity. Two strains of M. violaceum, from two different plant species but belonging to the same fungal lineage, contained copia-like elements with very low diversity, as would result from a high transposition rate, and these were also unique in showing no evidence of the hypermutation patterns indicative of the RIP genome defense. In this species, evidence of RIP was also absent from a Class II helitron-like transposable element. However, unexpectedly the absolute repetitive element load was lower than in other strains.

Detection of gyrA mutations in quinolone-resistant Salmonella enterica by denaturing high-performance liquid chromatography

Denaturing high-performance liquid chromatography (DHPLC) was evaluated as a rapid screening and identification method for DNA sequence variation detection in the quinolone resistance-determining region of gyrA from Salmonella serovars. A total of 203 isolates of Salmonella were screened using this method. DHPLC analysis of 14 isolates representing each type of novel or multiple mutations and the wild type were compared with LightCycler-based PCR-gyrA hybridization mutation assay (GAMA) and single-strand conformational polymorphism (SSCP) analyses. The 14 isolates gave seven different SSCP patterns, and LightCycler detected four different mutations. DHPLC detected 11 DNA sequence variants at eight different codons, including those detected by LightCycler or SSCP. One of these mutations was silent. Five isolates contained multiple mutations, and four of these could be distinguished from the composite sequence variants by their DHPLC profile. Seven novel mutations were identified at five different loci not previously described in quinolone-resistant salmonella. DHPLC analysis proved advantageous for the detection of novel and multiple mutations. DHPLC also provides a rapid, high-throughput alternative to LightCycler and SSCP for screening frequently occurring mutations.

Mapping the platelet profile for functional genomic studies and demonstration of the effect size of the GP6 locus

BACKGROUND: Evidence suggests the wide variation in platelet response within the population is genetically controlled. Unraveling the complex relationship between sequence variation and platelet phenotype requires accurate and reproducible measurement of platelet response. OBJECTIVE: To develop a methodology suitable for measuring signaling pathway-specific platelet phenotype, to use this to measure platelet response in a large cohort, and to demonstrate the effect size of sequence variation in a relevant model gene. METHODS: Three established platelet assays were evaluated: mobilization of [Ca(2+)](i), aggregometry and flow cytometry, each in response to adenosine 5'-diphosphate (ADP) or the glycoprotein (GP) VI-specific crosslinked collagen-related peptide (CRP). Flow cytometric measurement of fibrinogen binding and P-selectin expression in response to a single, intermediate dose of each agonist gave the best combination of reproducibility and inter-individual variability and was used to measure the platelet response in 506 healthy volunteers. Pathway specificity was ensured by blocking the main subsidiary signaling pathways. RESULTS: Individuals were identified who were hypo- or hyper-responders for both pathways, or who had differential responses to the two agonists, or between outcomes. 89 individuals, retested three months later using the same methodology, showed high concordance between the two visits in all four assays (r(2) = 0.872, 0.868, 0.766 and 0.549); all subjects retaining their phenotype at recall. The effect of sequence variation at the GP6 locus accounted for approximately 35% of the variation in the CRP-XL response. CONCLUSION: Genotyping-phenotype association studies in a well-characterized, large cohort provides a powerful strategy to measure the effect of sequence variation in genes regulating the platelet response.

Detection of somaclonal variation during cocoa somatic embryogenesis characterised using cleaved amplified polymorphic sequence and the new freeware Artbio

The scarcity and stochastic nature of genetic mutations presents a significant challenge for scientists seeking to characterise de novo mutation frequency at specific loci. Such mutations can be particularly numerous during regeneration of plants from in vitro culture and can undermine the value of germplasm conservation efforts. We used cleaved amplified polymorphic sequence (CAPS) analysis to characterise new mutations amongst a clonal population of cocoa plants regenerated via a somatic embryogenesis protocol used previously for cocoa cryopreservation. Efficacy of the CAPS system for mutation detection was greatly improved after an ‘a priori’ in silico screen of reference target sequences for actual and potential restriction enzyme recognition sites using a new freely available software called Artbio. Artbio surveys known sequences for existing restriction enzyme recognition sites but also identifies all single nucleotide polymorphism (SNP) deviations from such motifs. Using this software, we performed an in silico screen of seven loci for restriction sites and their potential mutant SNP variants that were possible from 21 restriction enzymes. The four most informative locus-enzyme combinations were then used to survey the regenerant populations for de novo mutants. We characterised the pattern of point mutations and, using the outputs of Artbio, calculated the ratio of base substitution in 114 somatic embryo-derived cocoa regenerants originating from two explant genotypes. We found 49 polymorphisms, comprising 26.3% of the samples screened, with an inferred rate of 2.8 × 10−3 substitutions/screened base. This elevated rate is of a similar order of magnitude to previous reports of de novo microsatellite length mutations arising in the crop and suggests caution should be exercised when applying somatic embryogenesis for the conservation of plant germplasm.

C. colchicum - variation measured by DNA sequence

Cyclamen colchicum has a mixed history in the hands of botanists. This paper examines the genetic identity of a group of wild Cyclamen populations from the Caucasus to discover whther they are Cyclamen colchicum, C. purpurascens or a mixture of the two. The collections supplemented by material collected at the type locality for C. colchicum, proved to be a single but variable genetic group of C. colchicum that was distinct from C. purpurascens.

Controlled variation of monomer sequence-distribution in the synthesis of aromatic poly(ether ketone)s

The effects of varying the alkali metal cation in the high-temperature nucleophilic synthesis of a semi-crystalline, aromatic poly(ether ketone) have been systematically investigated, and striking variations in the sequence-distributions and thermal characteristics of the resulting polymers were found. Polycondensation of 4,4'-dihydroxybenzophenone with 1,3-bis(4-fluorobenzoyl)benzene in diphenylsulfone as solvent, in the presence of an alkali metal carbonate M2CO3 (M= Li, Na, K, or Rb) as base, affords a range of different polymers that vary in the distribution pattern of 2-ring and 3-ring monomer units along the chain. Lithium carbonate gives an essentially alternating and highly crystalline polymer, but the degree of sequence-randomisation increases progressively as the alkali metal series is descended, with rubidium carbonate giving a fully random and non-thermally-crystallisable polymer. Randomisation during polycondensation is shown to result from reversible cleavage of the ether linkages in the polymer by fluoride ions, and an isolated sample of alternating-sequence polymer is thus converted to a fully randomised material on heating with rubidium fluoride.