Limited Phylogeographic Signal in Sex-Linked and Autosomal Loci Despite Geographically,Ecologically, and Phenotypically Concordant Structure of mtDNA Variation in the Holarctic Avian Genus Eremophila

Phylogeographic studies of Holarctic birds are challenging because they involve vast geographic scale, complex glacial history, extensive phenotypic variation, and heterogeneous taxonomic treatment across countries, all of which require large sample sizes. Knowledge about the quality of phylogeographic information provided by different loci is crucial for study design. We use sequences of one mtDNA gene, one sex-linked intron, and one autosomal intron to elucidate large scale phylogeographic patterns in the Holarctic lark genus Eremophila. The mtDNA ND2 gene identified six geographically, ecologically, and phenotypically concordant clades in the Palearctic that diverged in the Early - Middle Pleistocene and suggested paraphyly of the horned lark (E. alpestris) with respect to the Temminck''s lark (E. bilopha). In the Nearctic, ND2 identified five subclades which diverged in the Late Pleistocene. They overlapped geographically and were not concordant phenotypically or ecologically. Nuclear alleles provided little information on geographic structuring of genetic variation in horned larks beyond supporting the monophyly of Eremophila and paraphyly of the horned lark. Multilocus species trees based on two nuclear or all three loci provided poor support for haplogroups identified by mtDNA. The node ages calculated using mtDNA were consistent with the available paleontological data, whereas individual nuclear loci and multilocus species trees appeared to underestimate node ages. We argue that mtDNA is capable of discovering independent evolutionary units within avian taxa and can provide a reasonable phylogeographic hypothesis when geographic scale, geologic history, and phenotypic variation in the study system are too complex for proposing reasonable a priori hypotheses required for multilocus methods. Finally, we suggest splitting the currently recognized horned lark into five Palearctic and one Nearctic species.     

Is homoplasy or lineage sorting the source of incongruent mtdna and nuclear gene trees in the stiff-tailed ducks (Nomonyx-Oxyura)?

We evaluated the potential effects of homoplasy, ancestral polymorphism, and hybridization as obstacles to resolving phylogenetic relationships within Nomonyx-Oxyura stiff-tailed ducks (Oxyurinae; subtribe Oxyurina). Mitochondrial DNA (mtDNA) control region sequences from 94 individuals supported monophyly of mtDNA haplotypes for each of the six species and provided no evidence of extant incomplete lineage sorting or inter-specific hybridization. The ruddy ducks (O. j. jamaicensis,O. j. andina, O. j. ferruginea) are each others' closest relatives, but the lack of shared haplotypes between O. j. jamaicensis and O. j. ferruginea suggests long-standing historical isolation. In contrast, O. j. andina shares haplotypes with O. j. jamaicensis and O. j. ferruginea, which supports Todd's (1979) and Fjelds?'s (1986) hypothesis that O. j. andina is an intergrade or hybrid subspecies of O. j. jamaicensis and O. j. ferruginea. Control region data and a much larger data set composed of approximately 8800 base pairs of mitochondrial and nuclear sequence for each species indicate that the two New World species, O. vittata and O. jamaicensis, branch basally within Oxyura. A clade of three Old World species (O. australis, O. maccoa, O. leucocephala) is well supported, but different loci and also different characters within the mtDNA data support three different resolutions of the Old World clade, yielding an essentially unresolved trichotomy. Fundamentally different factors limited the resolution of the mtDNA and nuclear gene trees. Gene trees for most nuclear loci were unresolved due to slow rates of mutation and a lack of informative variation, whereas uncertain resolution of the mtDNA gene tree was due to homoplasy. Within the mtDNA, approximately equal numbers of characters supported each of three possible resolutions. Parametric and nonparametric bootstrap analyses suggest that resolution of the mtDNA tree based on ~4300 bp per taxon is uncertain but that complete mtDNA sequences would yield a fully resolved gene tree. A short internode separating O. leucocephala from (O. australis, O. maccoa) in the best mtDNA tree combined with long terminal branches and substantial rate variation among nucleotide sites allowed the small number of changes occurring on the internode to be obscured by homoplasy in a significant portion of simulated data sets. Although most nuclear loci were uninformative, two loci supported a resolution of the Old World clade (O. maccoa, O. leucocephala) that is incongruent with the best mtDNA tree. Thus, incongruence between nuclear and mtDNA trees may be due to random sorting of ancestral lineages during the short internode, homoplasy in the mtDNA data, or both. The Oxyura trichotomy represents a difficult though likely common problem in molecular systematics. Given a short internode, the mtDNA tree has a greater chance of being congruent with the history of speciation because its effective population size (N(e)) is one-quarter that of any nuclear locus, but its resolution is more likely to be obscured by homoplasy. In contrast, gene trees for more slowly evolving nuclear loci will be difficult to resolve due to a lack of substitutions during the internode, and when resolved are more likely to be incongruent with the species history due to the stochastic effects of lineage sorting. We suggest that researchers consider first whether independent gene trees are adequately resolved and then whether those trees are congruent with the species history. In the case of Oxyura, the answer to both questions may be no. Complete mtDNA sequences combined with data from a very large number of nuclear loci may be the only way to resolve such trichotomies.     

The impact of microsatellite electromorph size homoplasy on multilocus population structure estimates in a tropical tree (Corythophora alta) and an anadromous fish (Morone saxatilis)

Microsatellite allelic states are determined by electrophoretic sizing of polymerase chain reaction fragments to define electromorphs. Numerous studies have documented that identical microsatellite electromorphs are potentially heterogeneous at the DNA sequence level, a phenomenon called electromorph size homoplasy. Few studies have examined the impact of electromorph size homoplasy on estimates of population genetic parameters. We investigated the frequency of microsatellite electromorph size homoplasy for 12 loci in the tropical tree Corythophora alta and 11 loci in the anadromous fish Morone saxatilis by sequencing 14-23 homozygotes per locus sampled from multiple populations for a total of 453 sequences. Sequencing revealed no homoplasy for M. saxatilis loci. Seven C. alta loci exhibited homoplasy, including single and compound repeat motifs both with and without interruptions. Between 12.5 and 42.9% of electromorphs sampled per locus showed size homoplasy. Two methods of correction for homoplasy in C. alta generally produced little or no change in single-locus estimates of RST, except for two loci in which some additional differentiation among populations was revealed. Twelve-locus estimates of RST (including the seven loci corrected for homoplasy) were slightly greater than estimates from uncorrected data, although the 95% confidence intervals overlapped. The frequency of methodological errors such as clerical mistakes or sample mislabelling per genotype scored was estimated at 5.4 and 7.3% for C. alta and M. saxatilis, respectively. Simulations showed that the increase in RST produced by homoplasy correction was only slightly larger than variation in RST estimates expected to be caused by methodological errors.     

Impact of amplified fragment length polymorphism size homoplasy on the estimation of population genetic diversity and the detection of selective loci

AFLP markers are becoming one of the most popular tools for genetic analysis in the fields of evolutionary genetics and ecology and conservation of genetic resources. The technique combines a high-information content and fidelity with the possibility of carrying out genomewide scans. However, a potential problem with this technique is the lack of homology of bands with the same electrophoretic mobility, what is known as fragment-size homoplasy. We carried out a theoretical analysis aimed at quantifying the impact of AFLP homoplasy on the estimation of within- and between-neutral population genetic diversity in a model of a structured finite population with migration among subpopulations. We also investigated the performance of a currently used method (DFDIST software) to detect selective loci from the comparison between genetic differentiation and heterozygosis of dominant molecular markers, as well as the impact of AFLP homoplasy on its effectiveness. The results indicate that the biases produced by homoplasy are: (1) an overestimation of the frequency of the allele determining the presence of the band, (2) an underestimation of the degree of differentiation between subpopulations, and (3) an overestimation or underestimation of the heterozygosis, depending on the allele frequency of the markers. The impact of homoplasy is quickly diminished by reducing the number of fragments analyzed per primer combination. However, substantial biases on the expected heterozygosity (up to 15-25%) may occur with approximately 50-100 fragments per primer combination. The performance of the DFDIST software to detect selective loci from dominant markers is highly dependent on the number of selective loci in the genome and their average effects, the estimate of genetic differentiation chosen to be used in the analysis, and the critical bound probability used to detect outliers. Overall, the results indicate that the software should be used with caution. AFLP homoplasy can produce a reduction of up to 15% in the power to detect selective loci.     

Inferring phylogeny despite incomplete lineage sorting

It is now well known that incomplete lineage sorting can cause serious difficulties for phylogenetic inference, but little attention has been paid to methods that attempt to overcome these difficulties by explicitly considering the processes that produce them. Here we explore approaches to phylogenetic inference designed to consider retention and sorting of ancestral polymorphism. We examine how the reconstructability of a species (or population) phylogeny is affected by (a) the number of loci used to estimate the phylogeny and (b) the number of individuals sampled per species. Even in difficult cases with considerable incomplete lineage sorting (times between divergences less than 1 N(e) generations), we found the reconstructed species trees matched the "true" species trees in at least three out of five partitions, as long as a reasonable number of individuals per species were sampled. We also studied the tradeoff between sampling more loci versus more individuals. Although increasing the number of loci gives more accurate trees for a given sampling effort with deeper species trees (e.g., total depth of 10 N(e) generations), sampling more individuals often gives better results than sampling more loci with shallower species trees (e.g., depth = 1 N(e)). Taken together, these results demonstrate that gene sequences retain enough signal to achieve an accurate estimate of phylogeny despite widespread incomplete lineage sorting. Continued improvement in our methods to reconstruct phylogeny near the species level will require a shift to a compound model that considers not only nucleotide or character state substitutions, but also the population genetics processes of lineage sorting. [Coalescence; divergence; population; speciation.].     

Molecular resolution of the systematics of a problematic group of fishes (Teleostei: Osmeridae) and evidence for morphological homoplasy

Relationships among the species of Northern Hemisphere smelts (family Osmeridae) have long been debated in the fish systematics literature. Eight independent studies based on morphological characters failed to reach any consensus on osmerid interrelationships. We reconstruct the osmerid phylogeny based on DNA sequence data from three mitochondrial (cytb, 16S, 12S) and three nuclear (ITS2, S71, RAG1) gene regions from multiple individuals of the 14 species in 6 genera, using the Japanese ayu (Plecoglossus altivelis) as the outgroup. Analyses with different combinations of nuclear and mitochondrial datasets yielded a generally well-resolved phylogeny of the genera that conflicts with previous hypotheses of osmerid interrelationships, and Shimodaira–Hasegawa tests suggest our topology with the current molecular dataset is significantly better than earlier reconstructions. In addition, mapping 114 morphological characters used in previous studies onto our phylogeny shows widespread homoplasy, which is likely the source of the systematic disagreement produced in earlier works.     

Distinguishing the effects of selection from demographic history in the genetic variation of two sister passerines based on mitochondrial–nuclear comparison

Determining the mechanisms responsible for the distribution of genetic diversity in natural populations has occupied a central role in molecular evolution. Our study was motivated by the unprecedented observation that a widespread Eurasian flycatcher, Ficedula albicilla, exhibited no variation at the mitochondrial DNA (mtDNA) ND2 gene in 75 individuals sampled over a 5000-km distance. In contrast, its sister species, F. parva, had low but considerably higher levels of mtDNA variation. We assessed whether natural selection or demographic factors could explain the absence of mtDNA variation in F. albicilla. Eighteen nuclear genes were sequenced to estimate the two species'' phylogeographic histories, and for comparison to the mtDNA data. Multilocus coalescence analyses suggested that F. albicilla experienced a population expansion perhaps following a population bottleneck. Simulations based on this demographic history, however, did not replicate the extremely low level of mtDNA variation. Historical range changes based on ecological niche models also failed to explain the observed mtDNA patterns. Neutrality tests (DHEW and ML-HKA) suggested a non-neutral pattern in the mtDNA of F. albicilla. We found a transmembrane-skewed distribution of nonsynonymous substitutions between the two species, three of which caused functional change; the results implied that positive selection could have targeted mtDNA. Several lines of evidence support selection rather than demographic history as the main force influencing the patterns of mtDNA variation. Despite the influence of natural selection, many of the phylogeographic inferences derived from mtDNA were robust, including species limits and a high level of gene flow among populations within species.     

PARSIMONY,MOLECULAR EVOLUTION,AND BIOGEOGRAPHY: THE CASE OF THE NORTH AMERICAN GIANT SALAMANDER

To draw biogeographic conclusions about the Central Highlands region of the United States, we reconstructed the phylogeny of hellbender (Cryptobranchus alleganiensis) populations from restriction-site variation in mtDNA. We were unable to root the phylogeny using an outgroup and therefore could not weight restriction-site gains more heavily than site losses. As a result, maximum parsimony results in low phylogenetic resolution because of high levels of homoplasy in the data set. Use of a recently published algorithm based on an explicit model of molecular evolution yielded much greater resolution of the mtDNA relationships. This phylogeny indicates the two subspecies of hellbenders are paraphyletic with respect to one another. Hellbenders found in the southern Ozarks (C. a. bishopi) are either most closely related to populations of C. a. alleganiensis inhabiting the Tennessee River drainage or are so divergent that phylogenetic affinities are undetectable. Extremely low levels of divergence among mtDNA haplotypes found in populations from Pennsylvania, Indiana, Illinois, and the northern Missouri Ozarks suggest a recent, probably post-Pleistocene, invasion of this region from a refugium in one of these areas. Biogeographic hypotheses of the causes and timing of hellbender distributions differ significantly from those postulated from analyses of fish species relationships. Possible reasons for the discrepancy are discussed.     

Amplified fragment length polymorphism analysis identifies hybrids between two subspecies of warblers

In the western Pyrenees (Southwest France and Northwest Spain), a narrow hybrid zone exists between the common chiffchaff Phylloscopus (collybita) collybita and the Iberian chiffchaff Phylloscopus (c.) brehmii. In this zone, which is approximately 20 km wide, mixed matings and individuals singing the songs of both taxa occur at substantial frequencies (24 and 8.6%, respectively), suggesting frequent hybridization. Previous studies have shown very weak mitochondrial gene flow (Nm = 0.065), whereas four microsatellites suggested much higher nuclear gene flow (Nm = 4.9). In this study we used the amplified fragment length polymorphism (AFLP) method in order to identify hybrids and early backcrosses. We typed 91 birds from both allopatric and sympatric areas for 12 informative AFLP markers (of > 141 polymorphic fragments), obtained by screening 13 AFLP primer combinations. These individuals were previously typed for song (brehmii, collybita or mixed singers), mitochondrial DNA (mtDNA) haplotype and allelic genotypes at four microsatellite loci. Assignment tests demonstrated that in the zone of sympatry, a substantial number of intermediate genotypes existed among the birds previously believed to be pure collybita and brehmii, based on song and mtDNA haplotype. The majority of the mixed singers had intermediate genotypes. Our data suggest that the fraction of the adult population having a hybrid origin (hybrids or backcrosses) is in the order of 10%. With such a frequency of genetic hybrids, there would have been much more mtDNA introgression than observed, had female hybrids been perfectly fertile/viable. This result is consistent with male-biased gene flow and Haldane's rule.     

An efficient method for screening faecal DNA genotypes and detecting new individuals and hybrids in the red wolf (<Emphasis Type="Italic">Canis rufus</Emphasis>) experimental population area

Previously, sequencing of mitochondrial DNA (mtDNA) from non-invasively collected faecal material (scat) has been used to help manage hybridization in the wild red wolf (Canis rufus) population. This method is limited by the maternal inheritance of mtDNA and the inability to obtain individual identification. Here, we optimize the use of nuclear DNA microsatellite markers on red wolf scat DNA to distinguish between individuals and detect hybrids. We develop a data filtering method in which scat genotypes are compared to known blood genotypes to reduce the number of PCR amplifications needed. We apply our data filtering method and the more conservative maximum likelihood ratio method (MLR) of Miller et al. (2002 Genetics 160:357–366) to a scat dataset previously screened for hybrids by sequencing of mtDNA. Using seven microsatellite loci, we obtained genotypes for 105 scats, which were matched to 17 individuals. The PCR amplification success rate was 50% and genotyping error rates ranged from 6.6% to 52.1% per locus. Our data filtering method produced comparable results to the MLR method, and decreased the time and cost of analysis by 25%. Analysis of this dataset using our data filtering method verified that no hybrid individuals were present in the Alligator River National Wildlife Refuge, North Carolina in 2000. Our results demonstrate that nuclear DNA microsatellite analysis of red wolf scats provides an efficient and accurate approach to screen for new individuals and hybrids.     

Combining multiple autosomal introns for studying shallow phylogeny and taxonomy of Laurasiatherian mammals: Application to the tribe Bovini (Cetartiodactyla,Bovidae)

Mitochondrial sequences are widely used for species identification and for studying phylogenetic relationships among closely related species or populations of the same species. However, many studies of mammals have shown that the maternal history of the mitochondrial genome can be discordant with the true evolutionary history of the taxa. In such cases, the analyses of multiple nuclear genes can be more powerful for deciphering interspecific relationships.Here, we designed primers for amplifying 13 new exon-primed intron-crossing (EPIC) autosomal loci for studying shallow phylogeny and taxonomy of Laurasiatherian mammals. Three criteria were used for the selection of the markers: gene orthology, a PCR product length between 600 and 1200 nucleotides, and different chromosomal locations in the bovine genome. Positive PCRs were obtained from different species representing the orders Carnivora, Cetartiodactyla, Chiroptera, Perissodactyla and Pholidota.The newly developed markers were analyzed in a phylogenetic study of the tribe Bovini (the group containing domestic and wild cattle, bison, yak, African buffalo, Asian buffalo, and saola) based on 17 taxa and 18 nuclear genes, representing a total alignment of 13,095 nucleotides. The phylogenetic results were compared to those obtained from analyses of the complete mitochondrial genome and Y chromosomal genes. Our analyses support a basal divergence of the saola (Pseudoryx) and a sister-group relationship between yak and bison. These results contrast with recent molecular studies but are in better agreement with morphology. The comparison of pairwise nucleotide distances shows that our nuDNA dataset provides a good signal for identifying taxonomic levels, such as species, genera, subtribes, tribes and subfamilies, whereas the mtDNA genome fails because of mtDNA introgression and higher levels of homoplasy. Accordingly, we conclude that the genus Bison should be regarded as a synonym of Bos, with the European bison relegated to a subspecies rank within Bos bison. We compared our molecular dating estimates to the fossil record in order to propose a biogeographic scenario for the evolution of Bovini during the Neogene.     

Genetic introgression between an introduced babbler, the Chinese hwamei Leucodioptron c. canorum, and the endemic Taiwan hwamei L. taewanus: a multiple marker systems analysis

Identification of interspecific hybrids is often a subject of primary concern in the development of conservation strategies. Here we performed a genetic analysis combining mitochondrial DNA (mtDNA), microsatellites and single nucleotide polymorphic sites (SNPs) to assay the level of hybridization and introgression between an introduced babbler, Chinese hwamei Leucodioptron canorus, and its close relative, the endemic Taiwan hwamei L. taewanus in Taiwan. Fifty‐five Chinese hwameis from the Asian mainland and 69 Taiwan hwameis, including nine morphological hybrids, were sampled and analyzed. Evidences of mitochondrial introgression were found in three hybrids and one Taiwan hwamei. Five unlinked interspecific SNPs were identified at nine anonymous nuclear loci with interspecific differentiation (total Fst=0.77) that was much higher than that at seven highly polymorphic microsatellite loci combined (total Fst=0.1). Bayesian cluster analysis based on five interspecific SNP loci and two highly differentiated microsatellite loci (Fst>0.08) suggested that twelve individuals sampled in Taiwan were likely F2 or backcross hybrids, among which eight were morphological intermediates. A total of 20.3% (14/69) individuals sampled in Taiwan were suggested to be hybrids, suggesting that fitness reduction in hybrids might be negligible. These results imply that without an effective management strategy, the entire Taiwan hwamei population could easily become an admixed with Chinese hwamei and loose its evolutionary integrity. To reduce introgressive hybridization, illegal trade of Chinese hwamei should be strictly regulated and only the expensive male Chinese hwameis should be legally imported to minimize the chance for Chinese hwameis being released into the field. In our study we also found interspecific SNP markers to outperform microsatellite loci in detecting hybridization and introgression between two closely related species, which may be ascribed to the lower level of homoplasy of SNP loci.     

Regional structure despite limited mtDNA sequence diversity found in the endangered Huchen, <Emphasis Type="Italic">Hucho hucho</Emphasis> (Linnaeus, 1758)

We evaluate the hypothesis of no geographic structure in the Huchen (Hucho hucho), a large, predatory salmonid endemic to the Danube basin. Forty-seven individuals sampled from throughout the Huchen’s native range were genetically characterized. Extremely limited sequence diversity across 1,800 bases of mtDNA (the complete control region and partial NADH-1 subunit) evidenced by four closely related mtDNA haplotypes was found. Nonetheless, the geographic distribution of mtDNA repeats (5–10, 82-bp long copies per individual) as well as allelic diversity across two microsatellite loci indicated large-scale geographic structure between the north-western (Austria and Slovenia) distribution area and eastern (Slovakia and Ukraine) or southern (Bosnia-Herzegovina and Montenegro) sample sites. An extremely slow rate of substitution for the H. hucho mtDNA is considered along with the alternative hypotheses to explain the limited mtDNA diversity. Considering the regional genetic structure implied by our data, we advocate restrictions on the transport of brood fish or yearlings across the range of the species distribution and sale of Huchen across international boundaries. Future genetic analysis to support local conservation and monitoring efforts must focus on developing a high-resolution screen that may be applied to identify hatchery versus naturally reproduced individuals in the wild.     

Molecular systematics of Xenocyprinae (teleostei: cyprinidae): taxonomy, biogeography, and coevolution of a special group restricted in East Asia

We surveyed mitochondrial DNA (mtDNA) sequence variation in the subfamily Xenocyprinae from China and used these data to estimate intraspecific, interspecific, and intergeneric phylogeny and assess biogeographic scenarios underlying the geographic structure of lineages. We sequenced 1140 bp of cytochrome b from 30 individuals of Xenocyprinae and one putative outgroup (Myxocypris asiaticus) and also sequenced 297 bp of ND4L, 1380 bp of ND4, 68 bp of tRNA(His), and 69 bp of tRNA(Ser) from 17 individuals of Xenocyprinae and the outgroup (M. asiaticus). We detected high levels of nucleotide variation among populations, species, and genera. The phylogenetic analysis suggested that Distoechodon hupeinensis might be transferred to the genus Xenocypris, the taxonomic status of the genus Plagiognathops might be preserved, and species of Xenocypris and Plagiognathops form a monophyletic group that is sister to the genus Distoechodon and Pseudobrama. The introgressive hybridization might occur among the populations of X. argentea and X. davidi, causing the two species to not be separated by mtDNA patterns according to their species identification, and the process and direction of hybridization are discussed. The spatial distributions of mtDNA lineages among populations of Xenocypris were compatible with the major geographic region, which indicated that the relationship between Hubei + Hunan and Fujian is closer than that between Hubei + Hunan and Sichuan. From a perspective of parasite investigation, our data suggested that the fauna of Hexamita in Xenocyprinae could be used to infer the phylogeny of their hosts.     

Local scale population genetic structure of Entomophthora muscae epidemics

Isolates of the entomopathogenic fungus Entomophthora muscae were obtained from houseflies sampled at five stables on Zealand during an epidemic in fall 2005. DNA fingerprints were generated from single conidia from 40 E. muscae isolates using the PCR-based method of Inter-Simple Sequence Repeats (ISSR). This resulted in fingerprint patterns consisting of about 50 fragments, of which, 14 were polymorphic. From the polymorphic loci we estimated the reproductive mode, genetic differentiation and gene diversity of E. muscae populations using the statistics of index of association, Weir & Cockerhams theta and Nei's analysis of gene diversity in subdivided populations. Our results revealed no significant differences in allele frequencies among the five populations. The index of association test rejected the null hypothesis of random mating, but the test for paired locus compatibility showed weak or no linkage between loci indicating recombination or homoplasy in the dataset. From this study we cannot exclude the possibility that the population genetic structure underlying E. muscae epidemics could be panmictic consisting of several lineages with a high level of reciprocal migration.     

Population stock structure of leatherback turtles (Dermochelys coriacea) in the Atlantic revealed using mtDNA and microsatellite markers

This study presents a comprehensive genetic analysis of stock structure for leatherback turtles (Dermochelys coriacea), combining 17 microsatellite loci and 763 bp of the mtDNA control region. Recently discovered eastern Atlantic nesting populations of this critically endangered species were absent in a previous survey that found little ocean-wide mtDNA variation. We added rookeries in West Africa and Brazil and generated longer sequences for previously analyzed samples. A total of 1,417 individuals were sampled from nine nesting sites in the Atlantic and SW Indian Ocean. We detected additional mtDNA variation with the longer sequences, identifying ten polymorphic sites that resolved a total of ten haplotypes, including three new variants of haplotypes previously described by shorter sequences. Population differentiation was substantial between all but two adjacent rookery pairs, and F _ST values ranged from 0.034 to 0.676 and 0.004 to 0.205 for mtDNA and microsatellite data respectively, suggesting that male-mediated gene flow is not as widespread as previously assumed. We detected weak (F _ST = 0.008 and 0.006) but significant differentiation with microsatellites between the two population pairs that were indistinguishable with mtDNA data. POWSIM analysis showed that our mtDNA marker had very low statistical power to detect weak structure (F _ST < 0.005), while our microsatellite marker array had high power. We conclude that the weak differentiation detected with microsatellites reflects a fine scale level of demographic independence that warrants recognition, and that all nine of the nesting colonies should be considered as demographically independent populations for conservation. Our findings illustrate the importance of evaluating the power of specific genetic markers to detect structure in order to correctly identify the appropriate population units to conserve.     

Molecular phylogeny of Salmo of the western Balkans,based upon multiple nuclear loci

Background

Classification of species within the genus Salmo is still a matter of discussion due to their high level of diversity and to the low power of resolution of mitochondrial (mt)DNA-based phylogeny analyses that have been traditionally used in evolutionary studies of the genus. We apply a new marker system based on nuclear (n)DNA loci to present a novel view of the phylogeny of Salmo representatives and we compare it with the mtDNA-based phylogeny.

Methods

Twenty-two nDNA loci were sequenced for 76 individuals of the brown trout complex: Salmo trutta (Danubian, Atlantic, Adriatic, Mediterranean and Duero mtDNA lineages), Salmo marmoratus (marble trout), Salmo obtusirostris (softmouth trout), and Salmo ohridanus (Ohrid belvica or belushka). Sequences were phylogenetically analyzed using maximum-likelihood and Bayesian Inference methods. The divergence time of the major clades was estimated using the program BEAST.

Results

The existence of five genetic units i.e. S. salar, S. ohridanus, S. obtusirostris, S. marmoratus and the S. trutta complex, including its major phylogenetic lineages was confirmed. Contrary to previous observations, S. obtusirostris was found to be sister to the S. trutta complex and the S. marmoratus clade rather than to the S. ohridanus clade. Reticulate evolution of S. obtusirostris was confirmed and a time for its pre-glacial origin suggested. S. marmoratus was found to be a separate species as S. trutta and S. obtusirostris. Relationships among lineages within the S. trutta complex were weakly supported and remain largely unresolved.

Conclusions

Nuclear DNA-based results showed a fairly good match with the phylogeny of Salmo inferred from mtDNA analyses. The comparison of nDNA and mtDNA data revealed at least four cases of mitochondrial–nuclear DNA discordance observed that were all confined to the Adriatic basin of the Western Balkans. Together with the well-known extensive morphological and genetic variability of Balkan trouts, this observation highlights an interesting and variegated evolutionary history of Salmo in this area.     

Lack of phylogeographic structure in three widespread Australian birds reinforces emerging challenges in Australian historical biogeography

Aims Building on molecular studies of widespread Australian vertebrates, we tested whether each of three widespread Australian bird species, namely the singing honeyeater, Lichenostomus virescens, spiny‐cheeked honeyeater, Acanthagenys rufogularis (Passeriformes: Meliphagidae), and black‐faced woodswallow, Artamus cinereus (Passeriformes: Artamidae), has undergone a recent (Pleistocene) range expansion across the Australian continent. We related the findings to the presence or absence of geographic variation in each species’ external phenotype and whether historical or non‐historical factors have been involved in generating variation. Methods A total of 92 specimens of the three species were collected from, as far as possible, the same localities across Australia. They were sampled for mitochondrial DNA (mtDNA) diversity in the 1041 base pairs of the ND2 gene, and these data were analysed with nucleotide diversity statistics, unrooted networks, nested clade analysis, and tests of range expansion or stability. Results Range expansions could not be rejected in any of the species in our study. Each had low, geographically unstructured nucleotide diversity. Patterns of geographic variation in the singing honeyeater's and, to a lesser extent, the black‐faced woodswallow's external phenotypes are not correlated with mtDNA diversity in ND2. Main conclusions Our study adds to the increasing number of data sets suggesting the apparent prevalence of Pleistocene population expansions in widespread Australian birds. Furthermore, it shows that observable geographic structure may evolve very quickly, in response either to environmental gradients or to historical factors that operated too recently to be detected by ND2 sequences (e.g. in the singing honeyeater). Conversely, we have shown that a species that has had a recent population expansion need not necessarily be geographically invariant. To understand fully the interplay between vicariance and dispersal in the history of widespread Australian arid‐zone birds, or between the historical and non‐historical origins of their differentiation, carefully conducted case‐by‐case molecular studies will be necessary. Only then will biogeographical patterns and the processes that led to them emerge. Study of the historical biogeography and the more recent population history of Australian arid‐zone birds has reached a point where mtDNA‐based studies, while still informative and contributing to a growing data base of such work, should be complemented with data from multiple, rapidly evolving nuclear loci.     

Comprehensive phylogenomic time tree of bryophytes reveals deep relationships and uncovers gene incongruences in the last 500 million years of diversification

Premise

Bryophytes form a major component of terrestrial plant biomass, structuring ecological communities in all biomes. Our understanding of the evolutionary history of hornworts, liverworts, and mosses has been significantly reshaped by inferences from molecular data, which have highlighted extensive homoplasy in various traits and repeated bursts of diversification. However, the timing of key events in the phylogeny, patterns, and processes of diversification across bryophytes remain unclear.

Methods

Using the GoFlag probe set, we sequenced 405 exons representing 228 nuclear genes for 531 species from 52 of the 54 orders of bryophytes. We inferred the species phylogeny from gene tree analyses using concatenated and coalescence approaches, assessed gene conflict, and estimated the timing of divergences based on 29 fossil calibrations.

Results

The phylogeny resolves many relationships across the bryophytes, enabling us to resurrect five liverwort orders and recognize three more and propose 10 new orders of mosses. Most orders originated in the Jurassic and diversified in the Cretaceous or later. The phylogenomic data also highlight topological conflict in parts of the tree, suggesting complex processes of diversification that cannot be adequately captured in a single gene-tree topology.

Conclusions

We sampled hundreds of loci across a broad phylogenetic spectrum spanning at least 450 Ma of evolution; these data resolved many of the critical nodes of the diversification of bryophytes. The data also highlight the need to explore the mechanisms underlying the phylogenetic ambiguity at specific nodes. The phylogenomic data provide an expandable framework toward reconstructing a comprehensive phylogeny of this important group of plants.     

A review of the genetic relationships of Atlantic walrus (Odobenus rosmarus rosmarus) east and west of Greenland

Studies of the genetic variation involving allozymes, mitochondrial and nuclear DNA (microsatellites) in walruses (Odobenus rosmarus) were reviewed. In addition, the genetic relationships of a total of 211 Atlantic walruses, O. r. rosmarus, from 5 sampling areas west and east of Greenland were studied using 12 nuclear DNA-microsatellite loci and restriction fragment length polymorphism obtained from the ND1, ND2 and ND3/4 segments of the mitochondrial DNA (mtDNA). At the mtDNA level, no divergence was observed among the three sampling areas east of Greenland (i.e. East Greenland, Svalbard and Franz Josef Land), whereas areas west of Greenland (i.e. Northwest and West Greenland) showed some differentiation. The genetic variation at the microsatellite loci grouped the individuals into four sub-populations: Northwest Greenland, West Greenland, East Greenland and a common Svalbard-Franz Josef Land sub-population. A significant correlation between genetic distance and geographic distance between the sampling areas (isolation-by-distance effect) was detected, especially at the mtDNA level. At a small-scale phylo-geographical level, the mtDNA data indicated that Atlantic walruses have diverged into two major groups: one northwest (i.e. in the North Water) and one east of Greenland (i.e. an East Greenland-Svalbard-Franz Josef Land group), whereas the haplotype distribution in the West Greenland sample reflected a mixture of both these groups. The microsatellite data supported a general grouping of walruses to the west and east of Greenland.