An approximate Bayesian computation approach to overcome biases that arise when using amplified fragment length polymorphism markers to study population structure

There is great interest in using amplified fragment length polymorphism (AFLP) markers because they are inexpensive and easy to produce. It is, therefore, possible to generate a large number of markers that have a wide coverage of species genotnes. Several statistical methods have been proposed to study the genetic structure using AFLP's but they assume Hardy-Weinberg equilibrium and do not estimate the inbreeding coefficient, F-IS. A Bayesian method has been proposed by Holsinger and colleagues that relaxes these simplifying assumptions but we have identified two sources of bias that can influence estimates based on these markers: (i) the use of a uniform prior on ancestral allele frequencies and (ii) the ascertainment bias of AFLP markers. We present a new Bayesian method that avoids these biases by using an implementation based on the approximate Bayesian computation (ABC) algorithm. This new method estimates population-specific F-IS and F-ST values and offers users the possibility of taking into account the criteria for selecting the markers that are used in the analyses. The software is available at our web site (http://www-leca.uif-grenoble.fi-/logiciels.htm). Finally, we provide advice on how to avoid the effects of ascertainment bias.

Half tetrad analysis in alfalfa using multiple restriction fragment length polymorphism markers.

A maximum likelihood approach of half tetrad analysis (HTA) based on multiple restriction fragment length polymorphism (RFLP) markers was developed. This procedure estimates the relative frequencies of 2n gametes produced by mechanisms genetically equivalent to first division restitution (FDR) or second division restitution and simultaneously locates the centromere within a linkage group of RFLP marker loci. The method was applied to the diploid alfalfa clone PG-F9 (2n = 2x = 16) previously selected because of its high frequency of 2n egg production. HTA was based on four RFLP loci for which PG-F9 was heterozygous with codominant alleles that were absent in the tetraploid tester. Models including three linked and one unlinked RFLP loci were developed and tested. Results of the HTA showed that PG-F9 produced 6% FDR and 94% second division restitution 2n eggs. Information from a marker locus belonging to one linkage group was used to more precisely locate the centromere on a different linkage group. HTA, together with previous cytological analysis, indicated that in PG-F9, FDR 2n eggs are likely produced by diplospory, a mechanism common among apomictic species. The occurrence of FDR 2n eggs in plant species and their importance for crop evolution and breeding is discussed together with the potential applicability of multilocus HTA in the study of reproductive mutants.

A population genetics study of Anopheles darlingi (Diptera: Culicidae) from Colombia based on random amplified polymorphic DNA-polymerase chain reaction and amplified fragment lenght polymorphism markers

The genetic variation and population structure of three populations of Anopheles darlingi from Colombia were studied using random amplified polymorphic markers (RAPDs) and amplified fragment length polymorphism markers (AFLPs). Six RAPD primers produced 46 polymorphic fragments, while two AFLP primer combinations produced 197 polymorphic fragments from 71 DNA samples. Both of the evaluated genetic markers showed the presence of gene flow, suggesting that Colombian An. darlingi populations are in panmixia. Average genetic diversity, estimated from observed heterozygosity, was 0.374 (RAPD) and 0.309 (AFLP). RAPD and AFLP markers showed little evidence of geographic separation between eastern and western populations; however, the F ST values showed high gene flow between the two western populations (RAPD: F ST = 0.029; Nm: 8.5; AFLP: F ST = 0.051; Nm: 4.7). According to molecular variance analysis (AMOVA), the genetic distance between populations was significant (RAPD:phiST = 0.084; AFLP:phiST = 0.229, P < 0.001). The F ST distances and AMOVAs using AFLP loci support the differentiation of the Guyana biogeographic province population from those of the Chocó-Magdalena. In this last region, Chocó and Córdoba populations showed the highest genetic flow.

Restriction fragment length polymorphism of 195 bp repeated satellite dna of Trypanosoma cruzi supports the existence of two phylogenetic groups

The restriction fragment length polymorphism of the 195 bp repeated DNA sequence of Trypanosoma cruzi was analyzed among 23 T. cruzi stocks giving a reliable picture of the whole phylogenetic variability of the species. The profiles observed with the enzymes Hinf I and Hae III were linked together and supported the existence of two groups. Group 1 shows a 195 bp repeated unit (Hinf I) and high molecular weight DNA (Hae III), while group 2 presents a ladder profile for each enzyme, which is a characteristic of tandemly repeated DNA. The two groups, respectively, clustered stocks pertaining to the two principal lineages evidenced by isoenzyme and RAPD markers. The congruence among these three independent genomic markers corroborates the existence of two real phylogenetic lineages in T. cruzi. The specific monomorphic profiles for each major phylogenetic lineage suggest the existence of ancient sexuality and cryptic biological speciation.

Characterization of six rat strains (Rattus norvegicus) by mitochondrial DNA restriction fragment length polymorphism

Restriction fragment length polymorphism (RFLP) was used to examine the extent of mtDNA polymorphism among six strains of rats (Rattus norvegicus) - Wistar, Wistar Munich, Brown Norway, Wistar Kyoto, SHR and SHR-SP. A survey of 26 restriction enzymes has revealed a low level of genetic divergence among strains. The sites of cleavage by EcoRI, NcoI and XmnI were shown to be polymorphic. The use of these three enzymes allows the 6 strains to be classified into 4 haplotypes and identifies specific markers for each one. The percentage of sequence divergence among all pairs of haplotypes ranged from 0.035 to 0.33%, which is the result of a severe population constriction undergone by the strains. These haplotypes are easily demonstrable and therefore RFLP analysis can be employed for genetic monitoring of rats within animal facilities or among different laboratories.

Detection of the TLR4 1196C > T Polymorphism by Mismatched-Polymerase Chain Reaction Using Plasmid DNA as Internal Control in Restriction Fragment Length Polymorphism Assays

Background: Restriction fragment length polymorphism (RFLP) is a common molecular assay used for genotyping, and it requires validated quality control procedures to prevent mistyping caused by impaired endonuclease activity. We have evaluated the usefulness of a plasmid-based internal control in RFLP assays. Results: Blood samples were collected from 102 individuals with acute myocardial infarction (AMI) and 108 non-AMI individuals (controls) for DNA extraction and laboratory analyses. The 1196C> T polymorphism in the toll-like receptor 4 (TLR4) gene was amplified by mismatched-polymerase chain reaction (PCR). Amplicons and pBluescript II SK-plasmid were simultaneously digested with endonuclease HincII. Fragments were separated on 2% agarose gels. Plasmid was completely digested using up to 55.2 nmL/L DNA solutions and 1 mu L PCR product. Nevertheless, plasmid DNA with 41.4 nM or higher concentrations was incompletely digested in the presence of 7 mL PCR product. In standardized conditions, TLR4 1196C> T variant was accurately genotyped. TLR4 1196T allele frequency was similar between AMI (3.1%) and controls (2.0%, p = 0.948). TLR4 SNP was not associated with AMI in this sample population. In conclusion, the plasmid-based control is a useful approach to prevent mistyping in RFLP assays, and it is validate for genetic association studies such as TLR4 1196C> T.

Genotyping of Brazilian Erysipelothrix spp. strains by amplified fragment length polymorphism

One hundred and fifty-one Erysipelothrix spp. isolates from diseased and carrier swine from Brazil were identified by PCR, submitted to serotyping and analyzed by amplified fragment length polymorphism with a single enzyme (AFLP). Reference strains from Australia and the United Kingdom were also examined. The 151 strains were classified into 18 different serotypes (1a, 1b, 2a, 2b, 4, 5, 6, 7, 8, 10, 11, 12, 15, 17, 19, 21, 24 and 25), being serotype 2b the most frequent (39.7%). By associating serotyping and PCR results, it was possible to identify 146 strains as E. rhusiopathiae and five strains as E. tonsillarum. Despite the fact that for this genus AFLP did not cluster all isolates according to serotype, origin, disease or isolation data, the execution of the technique was easy and fast, demonstrating high discriminatory power. The results produced by the AFLP analysis of Erysipelothrix spp. could also support its use as a discriminatory tool for E. rhusiopathiae and E. tonsillarum species. (C) 2010 Elsevier B.V. All rights reserved.

Characterization by restriction fragment length polymorphism and sequence analysis of field and vaccine strains of infectious laryngotracheitis virus involved in severe outbreaks

At the end of 2002 and throughout 2003, there was a severe outbreak of infectious laryngotracheitis (ILT) in an intensive production area of commercial hens in the Sao Paulo State of Brazil. ILT virus was isolated from 28 flocks, and 21 isolates were genotyped by polymerase chain reaction and restriction fragment length polymorphism (PCR-RFLP) using four genes and eight restriction enzymes, and by partial sequencing of the infected cell protein 4 (ICP4) and thymidine kinase (TK) genes. Three groups resulted from the combinations of PCR-RFLP patterns: 19 field isolates formed Group I, and the remaining two isolates together with the chicken embryo origin (CEO) vaccine strains formed Group II. Group III comprised the tissue-culture origin (TCO) vaccine strain by itself. The PCR-RFLP results agreed with the sequencing results of two ICP4 gene fragments. The ICP4 gene sequence analysis showed that the 19 field isolates classified into Group I by RFLP-PCR were identical among themselves, but were different to the TCO and CEO vaccines. The two Group II isolates could not be distinguished from one of the CEO vaccines. The nucleotide and amino acid sequence analyses discriminated between the Brazilian and non-Brazilian isolates, as well as between the TCO and CEO vaccines. Sequence analysis of the TK gene enabled classification of the field isolates (Group I) as virulent and non-vaccine. This work shows that the severe ILT outbreak was caused by a highly virulent, non-vaccine strain.

Development and evaluation of a species diagnostic polymerase chain reaction-restriction fragment-length polymorphism procedure for cryptic members of the Culex sitiens (Diptera : Culicidae) subgroup in Australia and the southwest Pacific

Members of the Culex sitiens subgroup are important vectors of arboviruses, including Japanese encephalitis virus, Murray Valley encephalitis virus and Ross River virus. Of the eight described species, Cx. annulirostris Skuse, Cx. sitiens Wiedemann, and Cx. palpalis Taylor appear to be the most abundant and widespread throughout northern Australia and Papua New Guinea (PNG). Recent investigations using allozymes have shown this subgroup to contain cryptic species that possess overlapping adult morphology. We report the development of a polymerase chain reaction-restriction fragment-length polymorphism (PCR-RFLP) procedure that reliably separates these three species. This procedure utilizes the sequence variation in the ribosomal DNA ITS1 and demonstrates species-specific PCR-RFLP profiles from both colony and field collected material. Assessment of the consistency of this procedure was undertaken on mosquitoes sampled from a wide geographic area including Australia, PNG, and the Solomon Islands. Overlapping adult morphology was observed for Cx. annulirostris and Cx. palpalis in both northern Queensland and PNG and for all three species at one site in northwest Queensland.

The use of 16s rDNA restriction fragment length polymorphism analysis for the identification of non-fermenting gram negative bacilli in a cystic fibrosis microbiology referral service

The utility of 16s rDNA restriction fragment length polymorphism (RFLP) analysis for the partial genomovar differentiation of Burkholderia cepacia complex bacterium is well documented. We compared the 16s rDNA RFLP signatures for a number of non-fermenting gram negative bacilli (NF GNB) LMG control strains and clinical isolates pertaining to the genera Burkholderia, Pseudomonas, Achromobacter (Alcaligenes), Ralstonia, Stenotrophomonas and Pandoraea. A collection of 24 control strain (LMG) and 25 clinical isolates were included in the study. Using conventional PCR, a 1.2 kbp 16s rDNA fragment was generated for each organism. Following restriction digestion and electrophoresis, each clinical isolate RFLP signature was compared to those of the control strain panel. Nineteen different RFLP signatures were detected from the 28 control strains included in the study. TwentyoneyTwenty- five of the clinical isolates could be classified by RFLP analysis into a single genus and species when compared to the patterns produced by the control strain panel. Four clinical B. pseudomallei isolates produced RFLP signatures which were indistinguishable from B. cepacia genomovars I, III and VIII. The identity of these four isolates were confirmed using B. pseudomallei specific PCR. 16s rDNA RFLP analysis can be a useful identification strategy when applied to NF GNB, particularly for those which exhibit colistin sulfate resistance. The use of this molecular based methodology has proved very useful in the setting of a CF referral laboratory particularly when utilised in conjunction with B. cepacia complex and genomovar specific PCR techniques. Species specific PCR or sequence analysis should be considered for selected isolates; especially where discrepancies between epidemiology, phenotypic and genotypic characteristics occur.

Genetic relationships of Corynebacterium diphtheriae strains isolated from a diphtheria case and carriers by restriction fragment length polymorphism of rRNA genes

In the present study we report the results of an analysis, based on ribotyping of Corynebacterium diphtheriae intermedius strains isolated from a 9 years old child with clinical diphtheria and his 5 contacts. Quantitative analysis of RFLPs of rRNA was used to determine relatedness of these 7 C.diphtheriae strains providing support data in the diphtheria epidemiology. We have also tested those strains for toxigenicity in vitro by using the Elek's gel diffusion method and in vivo by using cell culture method on cultured monkey kidney cell (VERO cells). The hybridization results revealed that the 5 C.diphtheriae strains isolated from contacts and one isolated from the clinical case (nose case strain) had identical RFLP patterns with all 4 restriction endonucleases used, ribotype B. The genetic distance from this ribotype and ribotype A (throat case strain), that we initially assumed to be responsible for the illness of the patient, was of 0.450 showing poor genetic correlation among these two ribotypes. We found no significant differences concerned to the toxin production by using the cell culture method. In conclusion, the use of RFLPs of rRNA gene was successful in detecting minor differences in closely related toxigenic C.diphtheriae intermedius strains and providing information about genetic relationships among them.

Polymerase chain reaction and restriction fragment length polymorphism analysis of the ITS2 region for differatiation of Brazilian Biomphalaria intermediate hosts of the Schistosoma mansoni

We sequenced the internal transcribed spacer 2 of the ribosomal DNA (ITS2-DNAr) from the three Schistosoma mansoni intermediate hosts in Brazil: Biomphalaria glabrata, Biomphalaria tenagophila and Biomphalaria straminea. Analysis of a restriction map from those sequences allowed us to select putative restriction enzymes able to identify the snail species under study. Four restriction enzymes were used and HpaII provided simple species-specific profiles easily visualized in polyacrylamide gels. The use of ITS2 is advantageous as it provides a small fragment of 460 bp which may be easily amplified by PCR. In the current work, we showed that the amplification of ITS2-DNAr together with HpaII enzyme restriction is an auxiliary molecular tool for the morphological identification of such snails as well as for taxonomic and phylogenetic studies of neotropical planorbids.

Molecular characterization of Mycobacterium tuberculosis isolates from Tehran, Iran by restriction fragment length polymorphism analysis and spoligotyping

Abstract: INTRODUCTION Characterization of Mycobacterium tuberculosis (MTB) isolates by DNA fingerprinting has contributed to tuberculosis (TB) control. The aim of this study was to determine the genetic diversity of MTB isolates from Tehran province in Iran. METHODS MTB isolates from 60 Iranian and 10 Afghan TB patients were fingerprinted by standard IS6110-restriction fragment length polymorphism (RFLP) analysis and spoligotyping. RESULTS The copy number of IS6110 ranged from 10-24 per isolate. The isolates were classified into 22 clusters showing ≥ 80% similarity by RFLP analysis. Fourteen multidrug-resistant (MDR) isolates were grouped into 4 IS6110-RFLP clusters, with 10 isolates [71% (95% CI: 45-89%)] in 1 cluster, suggesting a possible epidemiological linkage. Eighteen Iranian isolates showed ≥ 80% similarity with Afghan isolates. There were no strains with identical fingerprints. Spoligotyping of 70 isolates produced 23 distinct patterns. Sixty (85.7%) isolates were grouped into 13 clusters, while the remaining 10 isolates (14.2%) were not clustered. Ural (formerly Haarlem4) (n = 22, 31.4%) was the most common family followed by Central Asian strain (CAS) (n = 18, 25.7%) and T (n = 9, 12.8%) families. Only 1strain was characterized as having the Beijing genotype. Among 60 Iranian and 10 Afghan MTB isolates, 25% (95% CI: 16-37) and 70% (95% CI: 39-89) were categorized as Ural lineage, respectively. CONCLUSIONS A higher prevalence of Ural family MTB isolates among Afghan patients than among Iranian patients suggests the possible transmission of this lineage following the immigration of Afghans to Iran.

Species identification in mammals from mixed biological samples based on mitochondrial DNA control region length polymorphism.

The ability to identify the species origin of an unknown biological sample is relevant in the fields of human and wildlife forensics. However, the detection of several species mixed in the same sample still remains a challenge. We developed and tested a new approach for mammal DNA identification in mixtures of two or three species, based on the analysis of mitochondrial DNA control region interspecific length polymorphism followed by direct sequencing. Contrary to other published methods dealing with species mixtures, our protocol requires a single universal primer pair and is not based on a pre-defined panel of species. Amplicons can be separated either on agarose gels or using CE. The advantages and limitations of the assay are discussed under different conditions, such as variable template concentration, amplicon sizes and size difference among the amplicons present in the mixture. For the first time, this protocol provides a simple, reliable and flexible method for simultaneous identification of multiple mammalian species from mixtures, without any prior knowledge of the species involved.

Molecular identification of similar species of the genus Biomphalaria (Mollusca: Planorbidae) determined by a polymerase chain reaction-restriction fragment length polymorphism

The freshwater snails Biomphalaria straminea, B. intermedia, B. kuhniana and B. peregrina, are morphologically similar; based on this similarity the first three species were therefore grouped in the complex B. straminea. The morphological identification of these species is based on characters such as vaginal wrinkling, relation between prepuce: penial sheath:deferens vas and number of muscle layers in the penis wall. In this study the polymerase chain reaction restriction fragment length polymorphism technique was used for molecular identification of these molluscs. This technique is based on the amplification of the internal transcribed spacer regions ITS1 e ITS2 of the ribosomal RNA gene and subsequent digestion of these fragments by restriction enzymes. Six enzymes were tested: Dde I, Mnl I, Hae III, Rsa I, Hpa II e Alu I. The restriction patterns obtained with DdeI presented the best profile for separation of the four species of Biomphalaria. The profiles obtained with all the enzymes were used to estimate the genetic distances among the species through analysis of common banding patterns.