Simple sequence repeats reveal uneven distribution of genetic diversity in chloroplast genomes of Brassica oleracea L. and (n=9) wild relatives

Diversity in the chloroplast genome of 171 accessions representing the Brassica 'C' (n = 9) genome, including domesticated and wild B. oleracea and nine inter-fertile related wild species, was investigated using six chloroplast SSR (microsatellite) markers. The lack of diversity detected among 105 cultivated and wild accessions of B. oleracea contrasted starkly with that found within its wild relatives. The vast majority of B. oleracea accessions shared a single haplotype, whereas as many as six haplotypes were detected in two wild species, B. villosa Biv. and B. cretica Lam.. The SSRs proved to be highly polymorphic across haplotypes, with calculated genetic diversity values (H) of 0.23-0.87. In total, 23 different haplotypes were detected in C genome species, with an additional five haplotypes detected in B. rapa L. (A genome n = 10) and another in B. nigra L. (B genome, n = 8). The low chloroplast diversity of B. oleracea is not suggestive of multiple domestication events. The predominant B. oleracea haplotype was also common in B. incana Ten. and present in low frequencies in B. villosa, B. macrocarpa Guss, B. rupestris Raf. and B. cretica. The chloroplast SSRs reveal a wealth of diversity within wild Brassica species that will facilitate further evolutionary and phylogeographic studies of this important crop genus.

Inter-simple sequence repeat and aggressiveness analyses revealed high genetic diversity, recombination and long-range dispersal in Fusarium culmorum

Inter-simple sequence repeat (ISSR) analysis and aggressiveness assays were used to investigate genetic variability within a global collection of Fusarium culmorum isolates. A set of four ISSR primers were tested, of which three primers amplified a total of 37 bands out of which 30 (81%) were polymorphic. The intraspecific diversity was high, ranging from four to 28 different ISSR genotypes for F. culmorum depending on the primer. The combined analysis of ISSR data revealed 59 different genotypes clustered into seven distinct clades amongst 75 isolates of F. culmorum examined. All the isolates were assayed to test their aggressiveness on a winter wheat cv. 'Armada'. A significant quantitative variation for aggressiveness was found among the isolates. The ISSR and aggressiveness variation existed on a macro- as well as micro-geographical scale. The data suggested a long-range dispersal of F. culmorum and indicated that this fungus may have been introduced into Canada from Europe. In addition to the high level of intraspecific diversity observed in F. culmorum, the index of multilocus association calculated using ISSR data indicated that reproduction in F. culmorum cannot be exclusively clonal and recombination is likely to occur.

A genetic linkage map of an apple rootstock progeny anchored to the Malus genome sequence

An apple rootstock progeny raised from the cross between the very dwarfing ‘M.27’ and the more vigorous ‘M.116’ (‘M.M.106’ × ‘M.27’) was used for the construction of a linkage map comprising a total of 324 loci: 252 previously mapped SSRs, 71 newly characterised or previously unmapped SSR loci (including 36 amplified by 33 out of the 35 novel markers reported here), and the self-incompatibility locus. The map spanned the 17 linkage groups (LG) expected for apple covering a genetic distance of 1,229.5 cM, an estimated 91% of the Malus genome. Linkage groups were well populated and, although marker density ranged from 2.3 to 6.2 cM/SSR, just 15 gaps of more than 15 cM were observed. Moreover, only 17.5% of markers displayed segregation distortion and, unsurprisingly in a semi-compatible backcross, distortion was particularly pronounced surrounding the self-incompatibility locus (S) at the bottom of LG17. DNA sequences of 273 SSR markers and the S locus, representing a total of 314 loci in this investigation, were used to anchor to the ‘Golden Delicious’ genome sequence. More than 260 of these loci were located on the expected pseudo-chromosome on the ‘Golden Delicious’ genome or on its homeologous pseudo-chromosome. In total, 282.4 Mbp of sequence from 142 genome sequence scaffolds of the Malus genome were anchored to the ‘M.27’ × ‘M.116’ map, providing an interface between the marker data and the underlying genome sequence. This will be exploited for the identification of genes responsible for traits of agronomic importance such as dwarfing and water use efficiency.

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.

Molecular genetic differentiation in earthworms inhabiting a heterogeneous Pb-polluted landscape

A Pb-mine site situated on acidic soil, but comprising of Ca-enriched islands around derelict buildings was used to study the spatial pattern of genetic diversity in Lumbricus rubellus. Two distinct genetic lineages ('A' and 'B'), differentiated at both the mitochondrial (mtDNA COII) and nuclear level (AFLPs) were revealed with a mean inter-lineage mtDNA sequence divergence of approximately 13%, indicative of a cryptic species complex. AFLP analysis indicates that lineage A individuals within one central 'ecological island' site are uniquely clustered, with little genetic overlap with lineage A individuals at the two peripheral sites. FTIR microspectroscopy of Pb-sequestering chloragocytes revealed different phosphate profiles in residents of adjacent acidic and calcareous islands. Bioinformatics found over-representation of Ca pathway genes in ESTPb libraries. Subsequent sequencing of a Ca-transport gene, SERCA, revealed mutations in the protein's cytosolic domain. We recommend the mandatory genotyping of all individuals prior to field-based ecotoxicological assays, particularly those using discriminating genomic technologies.

Genetic structure and systematic relationships within the Ophrys fuciflora aggregate (Orchidaceae: Orchidinae): high diversity in Kent and a wind-induced discontinuity bisecting the Adriatic

A recent phylogenetic study based on multiple datasets is used as the framework for a more detailed examination of one of the ten molecularly circumscribed groups identified, the Ophrys fuciflora aggregate. The group is highly morphologically variable, prone to phenotypic convergence, shows low levels of sequence divergence and contains an unusually large proportion of threatened taxa, including the rarest Ophrys species in the UK. The aims of this study were to (a) circumscribe minimum resolvable genetically distinct entities within the O. fuciflora aggregate, and (b) assess the likelihood of gene flow between genetically and geographically distinct entities at the species and population levels. Fifty-five accessions sampled in Europe and Asia Minor from the O. fuciflora aggregate were studied using the AFLP genetic fingerprinting technique to evaluate levels of infraspecific and interspecific genetic variation and to assess genetic relationships between UK populations of O. fuciflora s.s. in Kent and in their continental European and Mediterranean counterparts. The two genetically and geographically distinct groups recovered, one located in England and central Europe and one in south-eastern Europe, are incongruent with current species delimitation within the aggregate as a whole and also within O. fuciflora s.s. Genetic diversity is higher in Kent than in the rest of western and central Europe. Gene flow is more likely to occur between populations in closer geographical proximity than those that are morphologically more similar. Little if any gene flow occurs between populations located in the south-eastern Mediterranean and those dispersed throughout the remainder of the distribution, revealing a genetic discontinuity that runs north-south through the Adriatic. This discontinuity is also evident in other clades of Ophrys and is tentatively attributed to the long-term influence of prevailing winds on the long-distance distribution of pollinia and especially seeds. A cline of gene flow connects populations from Kent and central and southern Europe; these individuals should therefore be considered part of an extensive meta-population. Gene flow is also evident among populations from Kent, which appear to constitute a single metapopulation. They show some evidence of hybridization, and possibly also introgression, with O. apifera.

Genetic variation in a chilean endangered endemic: Gomortega keule (Molina) Baillon

Gomortega keule (Molina) Baillon is an endangered, rare species, the only representative of its genus, and endemic to Central Chile. Populations of this tree are now fragmented and few individuals can be found in any of them. Genetic diversity was studied in 33 individuals from three populations in Cauquenes, a coastal mountain area (35°58'S-72°41'W). Fifteen InterSimple Sequence Repeat primers were used to determine the degree of similarity between and within populations. This revealed that 30% of the variation exhibited was between populations while 70% was within; nevertheless individuals were clearly clustered in a pattern which reflected a narrow base of diversity. Three other species from the Laurales order were used in order to provide an external reference as to the degree of diversity. In addition, an external wild population from the native species, Peumus boldus, was used to verify the utility of the markers. We show that the primers are effective in quickly giving an estimate of the degree of diversity of a population, thus giving important topical information relevant to preserving endangered species. Aspects of the conservation and management policy for the species in order to maintain the remaining populations and to preserve the genetic resources are discussed.

A genetic linkage map of microsatellite, gene-specific and morphological markers in diploid Fragaria

Diploid Fragaria provide a potential model for genomic studies in the Rosaceae. To develop a genetic linkage map of diploid Fragaria, we scored 78 markers (68 microsatellites, one sequence-characterised amplified region, six gene-specific markers and three morphological traits) in an interspecific F2 population of 94 plants generated from a cross of F.vesca f. semperflorens × F. nubicola. Co-segregation analysis arranged 76 markers into seven discrete linkage groups covering 448 cM, with linkage group sizes ranging from 100.3 cM to 22.9 cM. Marker coverage was generally good; however some clustering of markers was observed on six of the seven linkage groups. Segregation distortion was observed at a high proportion of loci (54%), which could reflect the interspecific nature of the progeny and, in some cases, the self-incompatibility of F. nubicola. Such distortion may also account for some of the marker clustering observed in the map. One of the morphological markers, pale-green leaf (pg) has not previously been mapped in Fragaria and was located to the mid-point of linkage group VI. The transferable nature of the markers used in this study means that the map will be ideal for use as a framework for additional marker incorporation aimed at enhancing and resolving map coverage of the diploid Fragaria genome. The map also provides a sound basis for linkage map transfer to the cultivated octoploid strawberry.

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.

Progressive erosion of genetic and epigenetic variation in callus-derived cocoa (Theobroma cacao) plants

Relatively little is known about the timing of genetic and epigenetic forms of somaclonal variation arising from callus growth. We surveyed for both types of change in cocoa (Theobroma cacao) plants regenerated from calli of various ages, and also between tissues from the source trees. For genetic change, we used 15 single sequence repeat (SSR) markers from four source trees and from 233 regenerated plants. For epigenetic change, we used 386 methylation-sensitive amplified polymorphism (MSAP) markers on leaf and explant (staminode) DNA from two source trees and on leaf DNA from 114 regenerants. Genetic variation within source trees was limited to one slippage mutation in one leaf. Regenerants were far more variable, with 35% exhibiting at least one mutation. Genetic variation initially accumulated with culture age but subsequently declined. MSAP (epigenetic) profiles diverged between leaf and staminode samples from source trees. Multivariate analysis revealed that leaves from regenerants occupied intermediate eigenspace between leaves and staminodes of source plants but became progressively more similar to source tree leaves with culture age. Statistical analysis confirmed this rather counterintuitive finding that leaves of ‘late regenerants’ exhibited significantly less genetic and epigenetic divergence from source leaves than those exposed to short periods of callus growth.

Escherichia coli isolates from extraintestinal organs of livestock animals harbour diverse virulence genes and belong to multiple genetic lineages

Escherichia coli, the most common cause of bacteraemia in humans in the UK, can also cause serious diseases in animals. However the population structure, virulence and antimicrobial resistance genes of those from extraintestinal organs of livestock animals are poorly characterised. The aims of this study were to investigate the diversity of these isolates from livestock animals and to understand if there was any correlation between the virulence and antimicrobial resistance genes and the genetic backbone of the bacteria and if these isolates were similar to those isolated from humans. Here 39 E. coli isolates from liver (n=31), spleen (n=5) and blood (n=3) of cattle (n=34), sheep (n=3), chicken (n=1) and pig (n=1) were assigned to 19 serogroups with O8 being the most common (n=7), followed by O101, O20 (both n=3) and O153 (n=2). They belong to 29 multi-locus sequence types, 20 clonal complexes with ST23 (n=7), ST10 (n=6), ST117 and ST155 (both n=3) being most common and were distributed among phylogenetic group A (n=16), B1 (n=12), B2 (n=2) and D (n=9). The pattern of a subset of putative virulence genes was different in almost all isolates. No correlation between serogroups, animal hosts, MLST types, virulence and antimicrobial resistance genes was identified. The distributions of clonal complexes and virulence genes were similar to other extraintestinal or commensal E. coli from humans and other animals, suggesting a zoonotic potential. The diverse and various combinations of virulence genes implied that the infections were caused by different mechanisms and infection control will be challenging.

Comparative genomics of Brachyspira pilosicoli strains: genome rearrangements, reductions and correlation of genetic compliment with phenotypic diversity.

Background. The anaerobic spirochaete Brachyspira pilosicoli causes enteric disease in avian, porcine and human hosts, amongst others. To date, the only available genome sequence of B. pilosicoli is that of strain 95/1000, a porcine isolate. In the first intra-species genome comparison within the Brachyspira genus, we report the whole genome sequence of B. pilosicoli B2904, an avian isolate, the incomplete genome sequence of B. pilosicoli WesB, a human isolate, and the comparisons with B. pilosicoli 95/1000. We also draw on incomplete genome sequences from three other Brachyspira species. Finally we report the first application of the high-throughput Biolog phenotype screening tool on the B. pilosicoli strains for detailed comparisons between genotype and phenotype. Results. Feature and sequence genome comparisons revealed a high degree of similarity between the three B. pilosicoli strains, although the genomes of B2904 and WesB were larger than that of 95/1000 (~2,765, 2.890 and 2.596 Mb, respectively). Genome rearrangements were observed which correlated largely with the positions of mobile genetic elements. Through comparison of the B2904 and WesB genomes with the 95/1000 genome, features that we propose are non-essential due to their absence from 95/1000 include a peptidase, glycine reductase complex components and transposases. Novel bacteriophages were detected in the newly-sequenced genomes, which appeared to have involvement in intra- and inter-species horizontal gene transfer. Phenotypic differences predicted from genome analysis, such as the lack of genes for glucuronate catabolism in 95/1000, were confirmed by phenotyping. Conclusions. The availability of multiple B. pilosicoli genome sequences has allowed us to demonstrate the substantial genomic variation that exists between these strains, and provides an insight into genetic events that are shaping the species. In addition, phenotype screening allowed determination of how genotypic differences translated to phenotype. Further application of such comparisons will improve understanding of the metabolic capabilities of Brachyspira species.

Distribution, gene sequence and expression in vivo of the plasmid encoded fimbrial antigen of Salmonella serotype Enteritidis

The pefA gene which encoded the serotype associated plasmid (SAP) mediated fimbrial major subunit antigen of Salmonella enterica serotype Typhimurium shared genetic identity with 128 of 706 salmonella isolates as demonstrated by dot (colony) hybridization. Seventy-seven of 113 isolates of Typhimurium and individual isolates of serotypes Bovis-morbificans, Cholerae-suis and Enteritidis phage type 9b hybridized pefA strongly, whereas 48 isolates of Enteritidis hybridized pefA weakly and one Enteritidis isolate of phage type 14b failed to hybridize. Individual isolates of 294 serotypes and 247 individual isolates of serotype Dublin did not hybridize pefA. Southern hybridization of plasmids extracted from Enteritidis demonstrated that the pefA gene probe hybridized strongly an atypical SAP of 80 kb in size harboured by one Enteritidis isolate of phage-type 9b, whereas the typical SAP of 58 kb in size harboured by 48 Enteritidis isolates hybridized weakly. One Enteritidis isolate of phage type 14b which failed to hybridize pefA in dot (colony) hybridization experiments was demonstrated to be plasmid free. A cosmid library of Enteritidis phage type 4 expressed in Escherichia coli K12 was screened by hybridization for the presence of pef sequences. Recombinant clones which were deduced to harbour the entire pef operon elaborated a PEF-like fimbrial structure at the cell surface. The PEF-like fimbrial antigen was purified from one cosmid clone and used in western blot experiments with sera from chickens infected with Enteritidis phage-type 4. Seroconversion to the fimbrial antigen was observed which indicated that the Enteritidis PEF-like fimbrial structure was expressed at some stage during infection. Nucleotide sequence analysis demonstrated that the pefA alleles of Typhimurium and Enteritidis phage-type 4 shared 76% DNA nucleotide and 82% deduced amino acid sequence identity.

Genetic relatedness and phenotypic characteristics of Treponema associated with human periodontal tissues and ruminant foot disease

Treponema have been implicated recently in the pathogenesis of digital dermatitis (DID) and contagious ovine digital dermatitis (CODD) that are infectious diseases of bovine and ovine foot tissues, respectively. Previous analyses of treponemal 16S rDNA sequences, PCR-amplified directly from DID or CODD lesions, have suggested relatedness of animal Treponema to some human oral Treponema species isolated from periodontal tissues. In this study a range of adhesion and virulence-related properties of three animal Treponema isolates have been compared with representative human oral strains of Treponema denticola and Treponema vincentii. In adhesion assays using biotinylated treponemal cells, T denticola cells bound in consistently higher numbers to fibronectin, laminin, collagen type 1, gelatin, keratin and lactoferrin than did T. vincentii or animal Treponema isolates. However, animal DID strains adhered to fibrinogen at equivalent or greater levels than T denticola. All Treponema strains bound to the amino-terminal heparin l/fibrin I domain of fibronectin. 16S rDNA sequence analyses placed ovine strain UB1090 and bovine strain UB1467 within a cluster that was phylogenetically related to T vincentii, while ovine strain UB1466 appeared more closely related to T denticola. These observations correlated with phenotypic properties. Thus, T denticola ATCC 35405, GM-1, and Treponema UB1466 had similar outer-membrane protein profiles, produced chymotrypsin-like protease (CTLP), trypsin-like protease and high levels of proline iminopeptidase, and co-aggregated with human oral bacteria Porphyromonas gingivalis and Streptococcus crista. Conversely, T vincentii ATCC 35580, D2A-2, and animal strains UB1090 and UB1467 did not express CTLP or trypsin-like protease and did not co-aggregate with P. gingivalis or S. crista. Taken collectively, these results suggest that human oral-related Treponema have broad host specificity and that similar control or preventive strategies might be developed for human and animal Treponema-associated infections.

Genetic structure and linkage disequilibrium in landrace populations of barley in Sardinia

Multilocus digenic linkage disequilibria (LD) and their population structure were investigated in eleven landrace populations of barley (Hordeum vulgare ssp. vulgare L.) in Sardinia, using 134 dominant simple-sequence amplified polymorphism markers. The analysis of molecular variance for these markers indicated that the populations were partially differentiated (F ST = 0.18), and clustered into three geographic areas. Consistent with this population pattern, STRUCTURE analysis allocated individuals from a bulk of all populations into four genetic groups, and these groups also showed geographic patterns. In agreement with other molecular studies in barley, the general level of LD was low (13 % of locus pairs, with P < 0.01) in the bulk of 337 lines, and decayed steeply with map distance between markers. The partitioning of multilocus associations into various components indicated that genetic drift and founder effects played a major role in determining the overall genetic makeup of the diversity in these landrace populations, but that epistatic homogenising or diversifying selection was also present. Notably, the variance of the disequilibrium component was relatively high, which implies caution in the pooling of barley lines for association studies. Finally, we compared the analyses of multilocus structure in barley landrace populations with parallel analyses in both composite crosses of barley on the one hand and in natural populations of wild barley on the other. Neither of these serves as suitable mimics of landraces in barley, which require their own study. Overall, the results suggest that these populations can be exploited for LD mapping if population structure is controlled.