Widespread amplification of amplified fragment length polymorphisms (AFLPs) in marine Antarctic animals

Although recent years have witnessed a rapid growth in the number of genetic studies of Antarctic organisms, relatively few studies have so far used nuclear markers, possibly due to the perceived cost and difficulty of isolating markers such as microsatellites. However, an often overlooked alternative is to use amplified fragment length polymorphisms (AFLPs), a versatile and low-cost method capable of generating large numbers of predominantly nuclear loci in virtually any organism. We conducted a literature review of population genetic studies of Antarctic organisms, finding that fewer than 10% used AFLPs. Moreover, a strong taxonomic bias was found, with studies employing mitochondrial DNA or microsatellites focussing predominantly on animals, while those using AFLPs were mostly of plants or lower organisms. Consequently, we explored the extent to which AFLPs amplify across a range of Antarctic marine animal taxa by genotyping eight individuals each of twelve different species, ranging from echinoderms through soft corals to pelagic fish, at four selective primer combinations. AFLPs readily amplified across all of the taxa, generating between 32 and 84 loci per species, with on average 56.5% of these being polymorphic. In general, levels of polymorphism bore little relationship with expectations based on larger populations of broadcast-spawning species being more variable, though we did find a tentative positive correlation between the number of AFLP bands amplified and a measure of effective population size. Our study lends further support for the utility and ease of use of AFLPs and their suitability for studies of Antarctic species across a wide range of taxa.     

Gene genealogies, cryptic species, and molecular evolution in the human pathogen Coccidioides immitis and relatives (Ascomycota, Onygenales).

Previous genealogical analyses of population structure in Coccidioides immitis revealed the presence of two cryptic and sexual species in this pathogenic fungus but did not clarify their origin and relationships with respect to other taxa. By combining the C. immitis data with those of two of its closest relatives, the free-living saprophytes Auxarthron zuffianum and Uncinocarpus reesii, we show that the C. immitis species complex is monophyletic, indicating a single origin of pathogenicity. Cryptic species also were found in both A. zuffianum and U. reesii, indicating that they can be found in both pathogenic and free-living fungi. Our study, together with a few others, indicates that the current list of known fungal species might be augmented by a factor of at least two. However, at least in the C. immitis, A. zuffianum, and U. reesii complexes, cryptic species represent subdivisions at the tips of deep monophyletic clades and thus well within the existing framework of generic classification. An analysis of silent and expressed divergence and polymorphism values between and within the taxa identified by genealogical concordance did not reveal faster evolution in C. immitis as a consequence of adaptation to the pathogenic habit, nor did it show positive Darwinian evolution in a region of a dioxygenase gene (tcrP gene coding for 4-HPPD) known to cause antigenic responses in humans. Instead, the data suggested relative stasis, indicative of purifying selection against mostly deleterious mutations. Two introns in the same gene fragment were considerably more divergent than exons and were unalignable between species complexes but had very low polymorphism within taxa.     

New genes for phylogenetic studies of lichenized fungi: glyceraldehyde-3-phosphate dehydrogenase and beta-tubulin genes

Abstract:Primers for amplification and sequencing of partial glyceraldehyde-3-phosphate dehydrogenase (gpd) gene were designed for lichenized fungi. The 5′ gpd primer is most probably fungal specific, since a BLAST search in GenBank found identical sequences only from ascomycetous taxa, whereas the 3′ gpd primer was more universal. Utility of the gpd primers and previously designed beta-tubulin primers was tested in nine lichen taxa. Both the gpd and beta-tubulin primer pairs amplified in most of the taxa examined: the gpd primers generated a c. 1100 nucleotide fragment, whereas the PCR product obtained from the beta-tubulin primers was c. 900 nucleotides long. The gpd amplification products of Cladonia arbuscula and C. rangiferina were sequenced and both were found to contain three introns, the length of which varied between 49 to 83 nucleotides. To examine the applicability of gpd sequences in resolving relationships within Ascomycota, trees were calculated from 22 fungal gpd sequences obtained from GenBank together with the twoCladonia sequences using parsimony jackknifing. The gpd tree was compared with the SSU rDNA tree of the respective species (or genera). A similar analysis of the beta-tubulin gene was not performed, because only a few beta-tublin sequences from the same taxa were available in GenBank. The gpd tree was well resolved but in conflict with the SSU rDNA tree. In contrast to the SSU rDNA tree, the gpd tree did not support the monophyly of the Ascomycota. Analysis of the combined data set produced a tree very similar to that of the SSU rDNA data. However, the relationship of Lecanorales to the other orders remained unresolved. Even though gpd and beta-tubulin are highly conserved proteins, the third codon positions and introns are variable and both genes have the potential for inferring phylogenetic relationships at the lower taxonomic levels in the lichenized fungi. The two genes may be useful even below species level, depending on the species investigated.     

True and false gharials: a nuclear gene phylogeny of crocodylia

The phylogeny of Crocodylia offers an unusual twist on the usual molecules versus morphology story. The true gharial (Gavialis gangeticus) and the false gharial (Tomistoma schlegelii), as their common names imply, have appeared in all cladistic morphological analyses as distantly related species, convergent upon a similar morphology. In contrast, all previous molecular studies have shown them to be sister taxa. We present the first phylogenetic study of Crocodylia using a nuclear gene. We cloned and sequenced the c-myc proto-oncogene from Alligator mississippiensis to facilitate primer design and then sequenced an 1,100-base pair fragment that includes both coding and noncoding regions and informative indels for one species in each extant crocodylian genus and six avian outgroups. Phylogenetic analyses using parsimony, maximum likelihood, and Bayesian inference all strongly agreed on the same tree, which is identical to the tree found in previous molecular analyses: Gavialis and Tomistoma are sister taxa and together are the sister group of Crocodylidae. Kishino-Hasegawa tests rejected the morphological tree in favor of the molecular tree. We excluded long-branch attraction and variation in base composition among taxa as explanations for this topology. To explore the causes of discrepancy between molecular and morphological estimates of crocodylian phylogeny, we examined puzzling features of the morphological data using a priori partitions of the data based on anatomical regions and investigated the effects of different coding schemes for two obvious morphological similarities of the two gharials.     

Amplified fragment length polymorphism (AFLP) markers reveal extra-pair parentage in a bird species: the bluethroat (Luscinia svecica)

We tested the use of amplified fragment length polymorphism (AFLP) to assess the frequency of extra-pair parentage in a bluethroat (Luscinia svecica namnetum) population. Thirty-six families totalling 162 nestlings were analysed. Using a combination of three primer pairs, we reached an exclusion probability of 93% for the population. This probability can reach 99% considering families independently. We revealed that extra-pair fertilizations are very common: 63.8% of all broods contain at least one extra-pair young, totalling 41.9% of all young analysed. However, with the technique and the three primer pairs used it was not possible to attribute the parentage exclusions to extra-pair paternity, maternity or both. As brood parasitism has never been reported in this species, it seems likely that the exclusions are due to extra-pair males. This study shows that dominant AFLP markers can be useful for studying the mating system of taxa for which no microsatellite primers are available. This technique allows the approximate estimation of parentage exclusions despite the fact that it is not possible to know which parent has to be excluded.     

Assessment of genetic diversity in 31 species of mangroves and their associates through RAPD and AFLP markers

Random amplified polymorphic DNA (RAPD) and amplified fragment length polymorphism (AFLP) markers were used to assess the genetic diversity in 31 species of mangroves and mangrove associates. Four AFLP primer combinations resulted in the amplification of 840 bands with an average of 210 bands per primer combination and 11 RAPD primers yielded 319 bands with an average of 29 bands per primer. The percentage of polymorphism detected was too high indicating the high degree of genetic variability in mangrove taxa both at inter- and intra-generic levels. In the dendrogram, species belonging to a particular family/ genus, taxa inhabiting similar habitats or having similar adaptations tended to be together. There were exceptions too; as many unrelated species of mangroves form clusters. The intrafamilial classification and inter-relationships of genera in the family Rhizophoraceae could be confirmed by molecular analysis. Both the markers RAPD and AFLP were found equally informative and useful for a better understanding of the genetic variability and genome relationships among mangroves and their associated species.     

The usefulness of amplified fragment length polymorphism markers for taxon discrimination across graduated fine evolutionary levels in Caribbean Anolis lizards

Fine-level taxon discrimination is important in biodiversity assessment and ecogeographical research. Genomic markers are often required for studies on closely related taxa, however, most existing mitochondrial and nuclear markers require prior knowledge of the genome and are impractical for use in small conservation projects. This study describes the application of amplified fragment length polymorphism (AFLP) to discriminate at four progressively finer evolutionary levels of Caribbean Anolis lizards from the central Lesser Antilles. AFLP is shown to be a rapid and effective method for discriminating between species. Separation increases with primer pair number and choice of primer combination appears to be noncritical. Initial population-level results show markedly less discriminatory power. A screening technique for the identification of population informative markers combining principal component and principal coordinate analyses is presented and assessed. Subsequent results show selected conspecific AFLP data to be remarkably congruent with those of mitochondrial DNA, microsatellite and morphological markers. The use of AFLP as a low-cost nuclear marker in species-level taxon discrimination is supported, whereas population level application demands further consideration.     

The pollination ecology of Grevillea beadleana McGillivray,an endangered shrub from northern New South Wales,Australia

Grevillea beadleana (Proteaceae) is an endangered species known from five populations in northern New South Wales, Australia. The reproductive ecology of G. beadleana was compared in two populations with a ten-fold difference in the number of plants. Grevillea beadleana was found to be self-compatible in both populations and an examination of pollen viability and stigma maturation revealed that the species is protandrous. Flowering within inflorescences is acropetallous. In the first season plants in the largest population produced approx. ten-fold more inflorescences than those in the smaller population and, although the number of flowers per inflorescence did not vary significantly between populations the first season, the larger population produced more fruit per inflorescence than the smaller population. However, fruit to flower ratios were less than 0.2 in both seasons and populations. In both populations the number of fruit was significantly greater at the proximal end of the inflorescence, where flowers open first, compared with medial and distal positions. Several bird species were observed visiting flowers, although few birds were recorded foraging at plants in the smaller population. Within both populations, birds tended to make more within- than between-plant visits. Self-compatibility, acropetally and proximal fruit-set, combined with the predominantly within-plant movement of honeyeaters, suggests inbreeding may be common within both populations of G. beadleana. Pollination and fruiting success are discussed for G. beadleana and breeding systems among rare and common taxa in Grevillea are reviewed.     

Substitution bias, rapid saturation, and the use of mtDNA for nematode systematics

Only relatively recently have researchers turned to molecular methods for nematode phylogeny reconstruction. Thus, we lack the extensive literature on evolutionary patterns and phylogenetic usefulness of different DNA regions for nematodes that exists for other taxa. Here, we examine the usefulness of mtDNA for nematode phylogeny reconstruction and provide data that can be used for a priori character weighting or for parameter specification in models of sequence evolution. We estimated the substitution pattern for the mitochondrial ND4 gene from intraspecific comparisons in four species of parasitic nematodes from the family Trichostrongylidae (38-50 sequences per species). The resulting pattern suggests a strong mutational bias toward A and T, and a lower transition/transversion ratio than is typically observed in other taxa. We also present information on the relative rates of substitution at first, second, and third codon positions and on relative rates of saturation of different types of substitutions in comparisons ranging from intraspecific to interordinal. Silent sites saturate extremely quickly, presumably owing to the substitution bias and, perhaps, to an accelerated mutation rate. Results emphasize the importance of using only the most closely related sequences in order to infer patterns of substitution accurately for nematodes or for other taxa having strongly composition-biased DNA. ND4 also shows high amino acid polymorphism at both the intra- and interspecific levels, and in higher level comparisons, there is evidence of saturation at variable amino acid sites. In general, we recommend using mtDNA coding genes only for phylogenetics of relatively closely related nematode species and, even then, using only nonsynonymous substitutions and the more conserved mitochondrial genes (e.g., cytochrome oxidases). On the other hand, the high substitution rate in genes such as ND4 should make them excellent for population genetics studies, identifying cryptic species, and resolving relationships among closely related congeners when other markers show insufficient variation.      

Molecular analyses of Old World Leishmania RAPD markers and development of a PCR assay selective for parasites of the L. donovani species Complex

Three amplicons, appearing in a species-specific manner on the electrophoregrams of RAPD reactions that were obtained with primer OPA1, OPA1-800, OPA1- 900, and OPA1-1200, are analyzed in this study. The study revealed that each of these products is composed of one Leishmania DNA band, taxonomically conserved among the different Old World species studied. Subsequently, only the electrophoretic position of the RAPD products can be considered species-specific. In addition, sequence data, genomic organization, and chromosomal location have proved that these fragments are different and physically independent. However, they possess common features related to the presence of different kinds of short DNA repeats, more particularly microsatellites and a CCCTTC motive, corresponding to the 3' half of the OPA1 primer. These results suggest that the OPA1 primer has initiated amplification from different priming sites, having a species-specific location. This corresponds to sequence micro-heterogeneity of DNA fragments present within the different species and leading eventually to a selective amplification of different RAPD products. This characteristic has been used to develop an original selective PCR test based on the sequence of the OPA1-800 product, in which only DNAs from the L. donovani species complex are amplified. Restriction site polymorphisms and sequence variations are identified within the PCR fragment amplified from these parasite DNAs. In fact, the OPA1-800 fragment proved to be a useful DNA marker either as a DNA probe or as a target for PCR-based assays. This tool can therefore be recommended for the control of Old World Leishmania parasites, such as species discrimination, molecular tracking of isolates, or study of polymorphisms within the L. donovani species complex. Moreover, the molecular bases underlying the amplification of the RAPD fragments studied correspond to mechanisms already described. Although they do not account for the amplification of all Leishmania RAPD products, such mechanisms stress some of the pitfalls of the technique, which need to be taken into consideration. We have identified at least misleading observations of DNA bands amplified in a species-specific manner, in spite of their presence in the genome of the other taxa, and relatedness between bands within the amplification profiles. Therefore, recommendations for careful interpretation of RAPD data in population genetics or phylogenetic analyses are reiterated. Molecular analyses are essential to validate conclusions.     

The phylogenetic relationships of flies in the family Drosophilidae deduced from mtDNA sequences.

The phylogenetic relationships of several taxa from representative genera, subgenera, groups, and subgroups in the Drosophilidae were examined using sequences from a 905-bp mtDNA fragment. Conventional cloning and sequencing techniques were used to obtain nucleotide sequences. In addition, polymerase chain reaction primers were designed for the rapid amplification and sequencing of this region for the species examined in the Drosophilidae. Phylogenetic analysis was done by cladistic techniques. Because of the coding nature of the 905-bp mtDNA fragment, several separate analyses of these sequences were performed. The genera Scaptomyza and Hirtodrosophila occupy ancestral branching positions in the molecular phylogeny. The genera Chymomyza and Zaprionus have intermediate branching positions, while the subgenera Drosophila and Sophophora are in the most derived position in the molecular phylogeny. Within the subgenus Sophophora, there is little resolution using these sequences, while within the subgenus Drosophila, D. melanica, D. funebris, and D. pinicola form a clade in a derived part of the phylogeny, with D. robusta and D. immigrans branching in an intermediate position in the phylogeny. D. mercatorum, a member of the repleta species group, occupies an ancestral position in the molecular phylogeny.     

Phylogenetic relationships among anchovies, sardines, herrings and their relatives (Clupeiformes), inferred from whole mitogenome sequences

The relationships among and within the main lineages of the order Clupeiformes have been explored in few morphological studies and still remain poorly understood. Using whole mitogenome sequences, we inferred the relationships among 25 clupeiform species, sampled from each clupeiform family and subfamily, and a large selection of non-clupeiform teleosts. Our character sets, including unambiguously aligned, concatenated mitogenome sequences that we have divided into four (1st and 2nd codon positions, tRNA genes, and rRNA genes) or five partitions (same as before plus the transversions at 3rd codon positions, using 'RY' coding), were analyzed by the partitioned Bayesian method. The result strongly supported the monophyly of the Clupeiformes within the Otocephala, with Denticeps clupeoides as the sister group of a clade comprising all the remaining clupeiforms species (= suborder Clupeoidei). Within the Clupeoidei, the family Engraulidae was the sister group of the remaining taxa, comprising members of Sundasalangidae, Pristigasteridae, Clupeidae and Chirocentridae. Relationships among the latter four families remained ambiguous. In particular, the position of the Chirocentridae was difficult to estimate possibly owing to its higher molecular evolutionary rate. Of the five subfamilies in the family Clupeidae, monophylies of three (Alosinae, Clupeinae and Dorosomatinae) were statistically rejected. Instead, our mitogenomic data provide strong support for new clades within the Clupeidae, some of which are composed of members of more than one of the previously accepted subfamilies.     

Lower level relationships in the mushroom genus Cortinarius (Basidiomycota, Agaricales): a comparison of RPB1, RPB2, and ITS phylogenies

We sampled and analyzed approximately 2900bp across the three loci from 54 taxa belonging to a taxonomically difficult group of Cortinarius subgenus Phlegmacium. The combined analyses of ITS and variable regions of RPB1 and RPB2 greatly increase the resolution and nodal support for phylogenies of these closely related species belonging to clades that until now have proven very difficult to resolve with the ribosomal markers, nLSU and ITS. We present the first study of the utility of variable regions of the genes encoding the two largest subunits of RNA polymerase II (RPB1 and RPB2) for inferring the phylogeny of mushroom-forming fungi in combination with and compared to the widely used ribosomal marker ITS. The studied region of RPB1 contains an intron of the size and variability of ITS along with many variable positions in coding regions. Though almost entirely coding, the studied region of RPB2 is more variable than ITS. Both RNA polymerase II genes were alignable across all taxa. Our results indicate that several sections of Cortinarius need redefinition, and that several taxa treated at subspecific and varietal level should be treated at specific level. We suggest a new section for the two species, C. caesiocortinatus and C. prasinocyaneus, which constitute a well-supported separate lineage. We speculate that sequence information from RNA polymerase II genes have the potential for resolving phylogenetic problems at several levels of the diverse and taxonomically very challenging genus Cortinarius.     

Multiplex-terminal restriction fragment length polymorphism

A novel method called "multiplex-terminal restriction fragment length polymorphism (M-TRFLP)" has been recently developed which can be used for simultaneous analysis of the community composition of two or more microbial taxa (up to four). This method can also be used for microbial diagnostic purposes. For M-TRFLP analysis, primers specific to different target genes are used for multiplex-PCR, with one primer for each target being labeled with a unique fluorescent dye at its 5' end. Restriction digestion of the amplified products followed by fragment size analysis on a DNA sequencer produces profiles for targeted genes, which can be distinguished from each other by the color of the terminal fragments imparted by the unique fluorescent dye used for primer labeling. In contrast to current protocols, M-TRFLP allows multiple communities or multiple targets (genes) data to be obtained in just one reaction and therefore saves time, cost and labor. This protocol can be completed in 5-8 h.     

Evolution of functionally conserved enhancers can be accelerated in large populations: a population-genetic model

The evolution of cis-regulatory elements (or enhancers) appears to proceed at dramatically different rates in different taxa. Vertebrate enhancers are often very highly conserved in their sequences, and relative positions, across distantly related taxa. In contrast, functionally equivalent enhancers in closely related Drosophila species can differ greatly in their sequences and spatial organization. We present a population-genetic model to explain this difference. The model examines the dynamics of fixation of pairs of individually deleterious, but compensating, mutations. As expected, small populations are predicted to have a high rate of evolution, and the rate decreases with increasing population size. In contrast to previous models, however, this model predicts that the rate of evolution by pairs of compensatory mutations increases dramatically for population sizes above several thousand individuals, to the point of greatly exceeding the neutral rate. Application of this model predicts that species with moderate population sizes will have relatively conserved enhancers, whereas species with larger populations will be expected to evolve their enhancers at much higher rates. We propose that the different degree of conservation seen in vertebrate and Drosophila enhancers may be explained solely by differences in their population sizes and generation times.     

Characterization of a genetic resource collection for Miscanthus (Saccharinae,Andropogoneae, Poaceae) using AFLP and ISSR PCR

Amplified fragment length polymorphism (AFLP) and inter-simple sequence repeat markers were employed to characterize a genetic resource collection of Miscanthus, a grass under trial in Europe as a biomass crop. The 26 polymorphic markers produced by two ISSR fingerprinting primers were able to discriminate taxa and identify putative clones. AFLP fingerprints were fully reproducible and produced a larger number of markers for the three primer pairs tested, of which 998 were polymorphic (representing 79.3% of all bands). AFLP markers distinguished species, infra-specific taxa (varieties and cultivars) and putatively clonal material. They were also used to assess the inter-relationships of the taxa, to investigate the origin of important hybrid plants and to estimate the overall level of genetic variation in the collection. They were useful for assessing the species status of certain taxa such as M. transmorrisonensis, an endemic from Taiwan that was clearly distinct from M. sinensis; whereas other taxa of disputed species status, such as M. condensatus and M. yakushimanum were not genetically distinct from M. sinensis. The AFLP markers detected a high degree of infra-specific variation and allowed subdivisions of the genetic resource collection to be made, particularly within M. sinensis.     

Parallel declines in species and genetic diversity in tropical forest fragments

The potential for parallel impacts of habitat change on multiple biodiversity levels has important conservation implications. We report on the first empirical test of the 'species-genetic diversity correlation' across co-distributed taxa with contrasting ecological traits in the context of habitat fragmentation. In a rainforest landscape undergoing conversion to oil palm, we show that depauperate species richness in fragments is mirrored by concomitant declines in population genetic diversity in the taxon predicted to be most susceptible to fragmentation. This association, not seen in the other species, relates to fragment area rather than isolation. While highlighting the over-simplification of extrapolating across taxa, we show that fragmentation presents a double jeopardy for some species. For these, conserving genetic diversity at levels of pristine forest could require sites 15-fold larger than those needed to safeguard species numbers. Importantly, however, each fragment contributes to regional species richness, with larger ones tending to contain more species.     

The accuracy of methods for coding and sampling higher-level taxa for phylogenetic analysis: a simulation study

Many phylogenetic analyses, particularly morphological studies, use higher taxa (e.g., genera, families) rather than species as terminal taxa. This general approach requires dealing with interspecific variation among the species that make up the higher taxon. In this paper, I review different parsimony methods for coding and sampling higher taxa and compare their relative accuracies using computer simulations. Despite their widespread use, methods that involve coding higher taxa as terminals perform poorly in simulations, relative to splitting up the higher taxa and using species as terminals. Among the methods that use higher taxa as terminals, coding a taxon based on the most common condition among the included species (majority or modal coding) is generally more accurate than other coding methods, such as coding taxa as missing or polymorphic. The success of the majority method, and results of further simulations, suggest that in many cases "common equals primitive" within variable taxa, at least for low and intermediate rates of character change. The fixed-only method (excluding variable characters) performs very poorly, a result that is indirectly supported by analyses of published data for squamate reptiles. Sampling only a single species per higher taxon also yields low accuracy under many conditions. Along with recent studies of intraspecific polymorphism, the results of this study show the general importance of (1) including characters despite variation within taxa and (2) using methods that incorporate detailed information on the distribution of states within variable taxa.     

Exploring species limits in two closely related Chinese oaks

BACKGROUND: The species status of two closely related Chinese oaks, Quercus liaotungensis and Q. mongolica, has been called into question. The objective of this study was to investigate the species status and to estimate the degree of introgression between the two taxa using different approaches. METHODOLOGY/PRINCIPAL FINDINGS: Using SSR (simple sequence repeat) and AFLP (amplified fragment length polymorphism) markers, we found that interspecific genetic differentiation is significant and higher than the differentiation among populations within taxa. Bayesian clusters, principal coordinate analysis and population genetic distance trees all classified the oaks into two main groups consistent with the morphological differentiation of the two taxa rather than with geographic locations using both types of markers. Nevertheless, a few individuals in Northeast China and many individuals in North China have hybrid ancestry according to Bayesian assignment. One SSR locus and five AFLPs are significant outliers against neutral expectations in the interspecific FST simulation analysis, suggesting a role for divergent selection in differentiating species. MAIN CONCLUSIONS/SIGNIFICANCE: All results based on SSRS and AFLPS reached the same conclusion: Q. liaotungensis and Q. mongolica maintain distinct gene pools in most areas of sympatry. They should therefore be considered as discrete taxonomic units. Yet, the degree of introgression varies between the two species in different contact zones, which might be caused by different population history or by local environmental factors.