Species trees for the tree swallows (Genus Tachycineta): An alternative phylogenetic hypothesis to the mitochondrial gene tree

The New World swallow genus Tachycineta comprises nine species that collectively have a wide geographic distribution and remarkable variation both within- and among-species in ecologically important traits. Existing phylogenetic hypotheses for Tachycineta are based on mitochondrial DNA sequences, thus they provide estimates of a single gene tree. In this study we sequenced multiple individuals from each species at 16 nuclear intron loci. We used gene concatenated approaches (Bayesian and maximum likelihood) as well as coalescent-based species tree inference to reconstruct phylogenetic relationships of the genus. We examined the concordance and conflict between the nuclear and mitochondrial trees and between concatenated and coalescent-based inferences. Our results provide an alternative phylogenetic hypothesis to the existing mitochondrial DNA estimate of phylogeny. This new hypothesis provides a more accurate framework in which to explore trait evolution and examine the evolution of the mitochondrial genome in this group.     

Mitochondrial genome primers for Lake Malawi cichlids

Resolving the evolutionary history of rapidly diversifying lineages like the Lake Malawi Cichlid Flock demands powerful phylogenetic tools. Although this clade of over 500 species of fish likely diversified in less than two million years, the availability of extensive sequence data sets, such as complete mitochondrial genomes, could help resolve evolutionary patterns in this group. Using a large number of newly developed primers, we generated whole mitochondrial genome sequences for 14 Lake Malawi cichlids. We compared sequence divergence across protein‐coding regions of the mitochondrial genome and also compared divergence in the mitochondrial loci to divergence at two nuclear protein‐coding loci, Mitfb and Dlx2. Despite the widespread sharing of haplotypes of identical sequences at individual loci, the combined use of all protein‐coding mitochondrial loci provided a bifurcating phylogenetic hypothesis for the exemplars of major lineages within the Lake Malawi cichlid radiation. The primers presented here could have substantial utility for evolutionary analyses of mitochondrial evolution and hybridization within this diverse clade.     

LOW METABOLIC RATES IN SALAMANDERS ARE CORRELATED WITH WEAK SELECTIVE CONSTRAINTS ON MITOCHONDRIAL GENES

Mitochondria are the site for the citric acid cycle and oxidative phosphorylation (OXPHOS), the final steps of ATP synthesis via cellular respiration. Each mitochondrion contains its own genome; in vertebrates, this is a small, circular DNA molecule that encodes 13 subunits of the multiprotein OXPHOS electron transport complexes. Vertebrate lineages vary dramatically in metabolic rates; thus, functional constraints on mitochondrial‐encoded proteins likely differ, potentially impacting mitochondrial genome evolution. Here, we examine mitochondrial genome evolution in salamanders, which have the lowest metabolic requirements among tetrapods. We show that salamanders experience weaker purifying selection on protein‐coding sequences than do frogs, a comparable amphibian clade with higher metabolic rates. In contrast, we find no evidence for weaker selection against mitochondrial genome expansion in salamanders. Together, these results suggest that different aspects of mitochondrial genome evolution (i.e., nucleotide substitution, accumulation of noncoding sequences) are differently affected by metabolic variation across tetrapod lineages.     

The evolution of the adenylate-forming protein family in beetles: multiple luciferase gene paralogues in fireflies and glow-worms

Bioluminescence in beetles is dependent upon the enzyme luciferase. It has been hypothesised luciferase evolved from a fatty acyl-CoA synthetase gene deriving a novel bioluminescent function (neofunctionalization) after a gene duplication event. We evaluated this hypothesis within a phylogenetic framework using independent evidence obtained from the genome of Tribolium castaneum, published luciferase genes and novel luciferase and luciferase-like sequences. This phylogenetic study provides evidence for a large gene family of luciferase and luciferase-like paralogues in bioluminescent and non-bioluminescent beetles. All luciferase sequences formed a clade supporting a protoluciferase existing prior to the divergence of the Lampyridae, Elateridae and Phengodidae (Elateroidea). Multiple luciferase genes were identified from members of the Photurinae and the Luciolinae indicating complex gene duplication events within lampyrid genomes. The majority of luciferase residues were identified to be under purifying selection as opposed to positive selection. We conclude that beetle luciferase may have arisen from a process of subfunctionalization as opposed to neofunctionalization early on in the evolution of the Elateroidea.     

Genealogical concordance between the mating type locus and seven other nuclear genes supports formal recognition of nine phylogenetically distinct species within the Fusarium graminearum clade

Species limits were investigated within the Fusarium graminearum clade (Fg clade) through phylogenetic analyses of DNA sequences from portions of 11 nuclear genes including the mating-type (MAT) locus. Nine phylogenetically distinct species were resolved within the Fg clade, and they all possess contiguous MAT1-1 and MAT1-2 idiomorphs consistent with a homothallic reproductive mode. In contrast, only one of the two MAT idiomorphs was found in five other species, four of which were putatively asexual, and the other was heterothallic. Molecular evolutionary analyses indicate the MAT genes are under strong purifying selection and that they are functionally constrained, even in species for which a sexual state is unknown. The phylogeny supports a monophyletic and apomorphic origin of homothallism within this clade. Morphological analyses demonstrate that a combination of conidial characters could be used to differentiate three species and three species pairs. Species rank is formally proposed for the eight unnamed species within the Fg clade using fixed nucleotide characters.     

Molecular evolution of a family of resistance gene analogs of nucleotide-binding site sequences in <Emphasis Type="Italic">Solanum lycopersicum</Emphasis>

Nucleotide-binding site-leucine-rich repeats (NBS-LRR) gene families are one of the major plant resistance genes. Genomic NBS evolution was studied in many plant species for diverse arrays of NBS gene families. In this study, we focused on one family of NBS sequences in an attempt to understand how closely related NBS sequences evolved in the light of selection in domesticated plant species. A phylogenetic analysis revealed five major clades (A–E) and five subclades (A1–A5) within clade A of cloned NBS sequences. Positive selection was only detected in newly evolved NBS lineages in subclades of clade A. Positively selected codon sites were found among NBS sequences of clade A. A sliding-window analysis revealed that regions with Ka/Ks ratios of >1 were in the inter-motifs when paired clades were compared, but regions with Ka/Ks ratios of >1 were found across NBS sequences when subclades of clade A were compared. Our results based on a family of closely related NBS sequences showed that positive selection was first exerted on specific lineages across all NBS sequences after selective constraints. Subsequently, sequences with mutations in commonly conserved motifs were scrutinized by purifying selection. In the long term, conserved high frequency alleles in commonly conserved motifs and changes in inter-motifs were maintained in the investigated family of NBS sequences. Moreover, codons identified to be under positive selection in the inter-motifs were mainly located in regions involved in functions of ATP binding or hydrolysis.     

Two mitochondrial genes under episodic positive selection in subterranean octodontoid rodents

Tuco-tucos (Ctenomys) and related coruros (Spalacopus) are South American subterranean rodents. An energetically demanding lifestyle within the hypoxic, underground atmosphere may change the selective regime on oxidative phosphorylation. We examined whether weak and/or episodic positive directional selection affected the evolution of two mitochondrial genes (COX2, CytB), in a background of purifying selection in these lineages. We estimated rates of synonymous (dS) and non-synonymous (dN) substitutions and found: 1) significantly higher dN/dS ratio in subterranean groups relative to non-subterranean related species, and 2) two codons in each gene under episodic selection: 94 and 277 of COX2 and 269 and 307 of CytB.     

High altitude adaptation of the schizothoracine fishes (Cyprinidae) revealed by the mitochondrial genome analyses

The schizothoracine fishes, also known as “mountain carps” are widely distributed in the Qinghai-Tibetan Plateau and its peripheral regions. Although they provide a prime example of high altitude adaptation, the phylogenetic relationships and the divergence times among these carp lineages are still controversial. Moreover, the genetic basis for high altitude adaptation is also poorly understood. In this study, we determined the mitochondrial genomes from two species of the schizothoracine fishes, representing a “morphologically primitive” clade and “morphologically specialized” clade, respectively. The phylogenetic tree and the divergence times were estimated within the evolutionary framework of the entire order Cypriniformes. Our results indicate a polyphylyetic relationship of the schizothoracine fishes and suggest two independent migration events into the Qinghai-Tibetan Plateau: one by the “morphologically primitive” clade in the Late Miocene and another by the “morphologically specialized” clade in the Eocene. Rapid speciation events of each clade from the Late Miocene to the Pliocene correspond to the timing of the geologic acceleration of the Qinghai-Tibetan Plateau. Interestingly, we found evidence for positive selection acting on the protein coding genes in the mitochondrial genomes of the “morphologically specialized” clade, implying a possible genetic basis for high altitude adaptation in this derived lineage of cypriniform fishes.     

Ancient origin of elicitin gene clusters in Phytophthora genomes

The genus Phytophthora belongs to the oomycetes in the eukaryotic stramenopile lineage and is comprised of over 65 species that are all destructive plant pathogens on a wide range of dicotyledons. Phytophthora produces elicitins (ELIs), a group of extracellular elicitor proteins that cause a hypersensitive response in tobacco. Database mining revealed several new classes of elicitin-like (ELL) sequences with diverse elicitin domains in Phytophthora infestans, Phytophthora sojae, Phytophthora brassicae, and Phytophthora ramorum. ELIs and ELLs were shown to be unique to Phytophthora and Pythium species. They are ubiquitous among Phytophthora species and belong to one of the most highly conserved and complex protein families in the Phytophthora genus. Phylogeny construction with elicitin domains derived from 156 ELIs and ELLs showed that most of the diversified family members existed prior to divergence of Phytophthora species from a common ancestor. Analysis to discriminate diversifying and purifying selection showed that all 17 ELI and ELL clades are under purifying selection. Within highly similar ELI groups there was no evidence for positively selected amino acids suggesting that purifying selection contributes to the continued existence of this diverse protein family. Characteristic cysteine spacing patterns were found for each phylogenetic clade. Except for the canonical clade ELI-1, ELIs and ELLs possess C-terminal domains of variable length, many of which have a high threonine, serine, or proline content suggesting an association with the cell wall. In addition, some ELIs and ELLs have a predicted glycosylphosphatidylinositol site suggesting anchoring of the C-terminal domain to the cell membrane. The eli and ell genes belonging to different clades are clustered in the genomes. Overall, eli and ell genes are expressed at different levels and in different life cycle stages but those sharing the same phylogenetic clade appear to have similar expression patterns.     

Complete sequence of the mitochondrial genome of the tapeworm Hymenolepis diminuta: gene arrangements indicate that Platyhelminths are Eutrochozoans

Using "long-PCR," we amplified in overlapping fragments the complete mitochondrial genome of the tapeworm Hymenolepis diminuta (Platyhelminthes: Cestoda) and determined its 13,900-nt sequence. The gene content is the same as that typically found for animal mitochondrial DNA (mtDNA) except that atp8 appears to be lacking, a condition found previously for several other animals. Despite the small size of this mtDNA, there are two large noncoding regions, one of which contains 13 repeats of a 31-nt sequence and a potential stem-loop structure of 25 bp with an 11-member loop. Large potential secondary structures were identified also for the noncoding regions of two other cestode mtDNAS: Comparison of the mitochondrial gene arrangement of H. diminuta with those previously published supports a phylogenetic position of flatworms as members of the Eutrochozoa, rather than placing them basal to either a clade of protostomes or a clade of coelomates.     

Systematics of Mukdenia and Oresitrophe (Saxifragaceae): Insights from genome skimming data

Oresitrophe and Mukdenia (Saxifragaceae) are epilithic sister genera used in traditional Chinese medicine. The taxonomy of Mukdenia, especially of M. acanthifolia, has been controversial. To address this, we produced plastid and mitochondrial data using genome skimming for Mukdenia acanthifolia and Mukdenia rossii, including three individuals of each species. We assembled complete plastomes, mitochondrial CDS and nuclear ribosomal ETS/ITS sequences using these data. Comparative analysis shows that the plastomes of Mukdenia and Oresitrophe are relatively conservative in terms of genome size, structure, gene content, RNA editing sites and codon usage. Five plastid regions that represent hotspots of change (trnH-psbA, psbC-trnS, trnM-atpE, petA-psbJ and ccsA-ndhD) are identified within Mukdenia, and six regions (trnH-psbA, petN-psbM, trnM-atpE, rps16-trnQ, ycf1 and ndhF) contain a higher number of species-specific parsimony-informative sites that may serve as potential DNA barcodes for species identification. To infer phylogenetic relationships between Mukdenia and Oresitrophe, we combined our data with published data based on three different datasets. The monophyly of each species (Oresitrophe rupifraga, M. acanthifolia and M. rossii) and the inferred topology ((M. rossii, M. acanthifolia), O. rupifraga) are well supported in trees reconstructed using the complete plastome sequences, but M. acanthifolia and M. rossii did not form a separate clade in the trees based on ETS + ITS data, while the mitochondrial CDS trees are not well-resolved. We found low recovery of genes in the Angiosperms353 target enrichment panel from our unenriched genome skimming data. Hybridization or incomplete lineage sorting may be the cause of discordance between trees reconstructed from organellar and nuclear data. Considering its morphological distinctiveness and our molecular phylogenetic results, we strongly recommend that M. acanthifolia be treated as a distinct species.     

Lineage-specific selection in human mtDNA: lack of polymorphisms in a segment of MTND5 gene in haplogroup J

Human mitochondrial DNA (mtDNA) is a nonrecombining genome that codes for 13 subunits of the mitochondrial oxidative phosphorylation system, 2 rRNAs, and 22 tRNAs. Mutations have accumulated sequentially in mtDNA lineages that diverged tens of thousands of years ago. The genes in mtDNA are subject to different functional constraints and are therefore expected to evolve at different rates, but the rank order of these rates should be the same in all lineages of a phylogeny. Previous studies have indicated, however, that specific regions of mtDNA may have experienced different histories of selection in different lineages, possibly because of lineage-specific interactions or environmental factors such as climate. We report here on a survey for lineage-specific patterns of nucleotide polymorphism in human mtDNA. We calculated molecular polymorphism indices and neutrality tests for classes of functional sites and genes in 837 human mtDNA sequences, compared the results between continent-specific mtDNA lineages, and used two sliding window methods to identify differences in the patterns of polymorphism between haplogroups. A general correlation between nucleotide position and the level of nucleotide polymorphism was identified in the coding region of the mitochondrial genome. Nucleotide diversity in the protein-coding sequence of mtDNA was generally not much higher than that found for many genes in nuclear DNA. A comparison of nonsynonymous/synonymous rate ratios in the 13 protein-coding genes suggested differences in the relative levels of selection between haplogroups, including the European haplogroup clusters. Interestingly, a segment of the MTND5 gene was found to be almost void of segregating sites and nonsynonymous mutations in haplogroup J, which has been associated with susceptibility to certain complex diseases. Our results suggest that there are haplogroup-specific differences in the intensity of selection against particular regions of the mitochondrial genome, indicating that some mutations may be non-neutral within specific phylogenetic lineages but neutral within others.     

Complete mitochondrial genome sequence from an endangered Indian snake, Python molurus molurus (Serpentes, Pythonidae)

This paper reports the complete mitochondrial genome sequence of an endangered Indian snake, Python molurus molurus (Indian Rock Python). A typical snake mitochondrial (mt) genome of 17258 bp length comprising of 37 genes including the 13 protein coding genes, 22 tRNA genes, and 2 ribosomal RNA genes along with duplicate control regions is described herein. The P. molurus molurus mt. genome is relatively similar to other snake mt. genomes with respect to gene arrangement, composition, tRNA structures and skews of AT/GC bases. The nucleotide composition of the genome shows that there are more A-C % than T-G% on the positive strand as revealed by positive AT and CG skews. Comparison of individual protein coding genes, with other snake genomes suggests that ATP8 and NADH3 genes have high divergence rates. Codon usage analysis reveals a preference of NNC codons over NNG codons in the mt. genome of P. molurus. Also, the synonymous and non-synonymous substitution rates (ka/ks) suggest that most of the protein coding genes are under purifying selection pressure. The phylogenetic analyses involving the concatenated 13 protein coding genes of P. molurus molurus conformed to the previously established snake phylogeny.     

The mastodon mitochondrial genome: a mammoth accomplishment

The mitochondrial genome of an American mastodon was recently sequenced and used to root a phylogenetic analysis that included full mitochondrial genome sequences from woolly mammoths and the two living elephant genera. The study definitively established that mammoth and Asian elephant mitochondrial DNA lineages are more closely related than either is to African elephants. However, it also suggests that a complex evolutionary picture could ultimately emerge and points to similarities between the early evolution of the Elephantidae and that of the gorilla-human-chimpanzee clade.     

Analysis of nuclear copies of mitochondrial sequences in honeybee (Apis mellifera) genome

At least 0.08% of the Apis mellifera nuclear genome contains sequences that originated from mitochondria. These nuclear copies of mitochondrial sequences (numts) are scattered all over the honeybee chromosomes and have originated by multiple independent insertions of mitochondrial DNA (mtDNA) as evident by phylogenetic analysis. Apart from original insertions, moderate duplications of numts also contributed to the present pattern and distribution of mitochondrial sequences in honeybee chromosomes. Assimilation of mitochondrial genes in the nuclear genome is mediated by extensive fragmentations of the original inserts. Replication slippage seems to be a major mechanism by which small sequences are inserted or deleted from mtDNA destined to nucleus. Most of the honeybee numts (84%) are located in the nongenic regions. The majority (94%) of the numts that are located in predicted nuclear genes have originated from mitochondrial genes coding for cytochrome oxidase and NADH dehydrogenase subunits. On the other hand, the mitochondrial rRNA or tRNA gene sequences are predominantly (88%) located in nongenic regions of the genome. Evidences also support for exertion of purifying selection on numts located in specific genes. Comparative analysis of numts of European, African, and Africanized honeybees suggests that numt evolution in A. mellifera is probably not demarked by speciation time frame but may be a continuous and dynamic process.     

High-throughput sequencing of six bamboo chloroplast genomes: phylogenetic implications for temperate woody bamboos (Poaceae: Bambusoideae)

Background

Bambusoideae is the only subfamily that contains woody members in the grass family, Poaceae. In phylogenetic analyses, Bambusoideae, Pooideae and Ehrhartoideae formed the BEP clade, yet the internal relationships of this clade are controversial. The distinctive life history (infrequent flowering and predominance of asexual reproduction) of woody bamboos makes them an interesting but taxonomically difficult group. Phylogenetic analyses based on large DNA fragments could only provide a moderate resolution of woody bamboo relationships, although a robust phylogenetic tree is needed to elucidate their evolutionary history. Phylogenomics is an alternative choice for resolving difficult phylogenies.

Methodology/Principal Findings

Here we present the complete nucleotide sequences of six woody bamboo chloroplast (cp) genomes using Illumina sequencing. These genomes are similar to those of other grasses and rather conservative in evolution. We constructed a phylogeny of Poaceae from 24 complete cp genomes including 21 grass species. Within the BEP clade, we found strong support for a sister relationship between Bambusoideae and Pooideae. In a substantial improvement over prior studies, all six nodes within Bambusoideae were supported with ≥0.95 posterior probability from Bayesian inference and 5/6 nodes resolved with 100% bootstrap support in maximum parsimony and maximum likelihood analyses. We found that repeats in the cp genome could provide phylogenetic information, while caution is needed when using indels in phylogenetic analyses based on few selected genes. We also identified relatively rapidly evolving cp genome regions that have the potential to be used for further phylogenetic study in Bambusoideae.

Conclusions/Significance

The cp genome of Bambusoideae evolved slowly, and phylogenomics based on whole cp genome could be used to resolve major relationships within the subfamily. The difficulty in resolving the diversification among three clades of temperate woody bamboos, even with complete cp genome sequences, suggests that these lineages may have diverged very rapidly.     

Signatures of natural selection in the mitochondrial genomes of Tachycineta swallows and their implications for latitudinal patterns of the ‘pace of life’

Latitudinal variation in avian life histories can be summarized as a slow–fast continuum, termed the ‘pace of life’, that encompasses patterns in life span, reproduction, and rates of development among tropical and temperate species. Much of the variation in avian pace of life is tied to differences in rates of long-term metabolic energy expenditure. Given the vital role of the mitochondrion in metabolic processes, studies of variation in the mitochondrial genome may offer opportunities to establish mechanistic links between genetic variation and latitudinal ‘pace of life’ patterns. Using comparative genomic analyses, we examined complete mitochondrial genome sequences obtained from nine, broadly distributed Tachycineta swallow species to test for signatures of natural selection across the mitogenome within a phylogenetic framework. Our results show that although purifying selection is the dominant selective force acting on the mitochondrial genome in Tachycineta, three mitochondrial genes (ND2, ND5, and CYTB) contain regions that exhibit signatures of diversifying selection. Two of these genes (ND2 and ND5) encode interacting subunits of NADH dehydrogenase, and amino residues that were inferred to be targets of positive selection were disproportionately concentrated in these genes. Moreover, the positively selected sites exhibited a phylogenetic pattern that could be indicative of adaptive divergence between “fast” and “slow” lineages. These results suggest that functional variation in cytochrome b and NADH dehydrogenase could mechanistically contribute to latitudinal ‘pace of life’ patterns in Tachycineta.     

Multilocus phylogeny of the widely distributed South American lizard clade Eulaemus (Liolaemini,Liolaemus)

The lizard genus Liolaemus and different clades within it have been the focus of several recent phylogenetic studies mainly based on morphology and mtDNA. Although there is general consensus for recognizing two clades (subgenera) within the genus, [Liolaemus (sensu stricto) and Eulaemus], phylogenetic relationships within each subgenus remain difficult to elucidate, given incomplete taxonomic sampling and large discordance between published studies. Here, new phylogenetic relationships for the Eulaemus subgenus are proposed based on the largest molecular data set ever used for this clade, which includes 188 individuals and 14 loci representing different parts of the genome (mtDNA, anonymous nuclear loci and nuclear protein‐coding loci). This data set was analysed using two species tree approaches (*beast and MDC). Levels of discordance among methods were found, and with previously published studies, but results are robust enough to propose new phylogenetic hypotheses for the Eulaemus clade. Specifically well‐resolved and well‐supported novel hypotheses are provided within the lineomaculatus section, and we formally recognize the zullyae clade, the sarmientoi clade and the hatcheri group. We also resolve species relationships within the montanus section, and particularly within the melanops series. We found discordance between mitochondrial and nuclear trees and discussed alternative hypotheses for the lineomaculatus and montanus sections, as well as the challenge in resolving phylogenetic relationships for large clades in general.     

Phylogenetic relationships among middle American pocket gophers (genus Orthogeomys) based on mitochondrial DNA sequences.

Relationships among members representing each of the three subgenera of the Middle American rodent genus Orthogeomys (Rodentia: Geomyidae) were studied by comparing DNA sequence data from two regions of the mitochondrial genome. Results from 527 bp from the 16 S rDNA region and a 402-bp fragment of the cytochrome b gene indicate that the three subgenera are well differentiated genetically, with the subgenus Orthogeomys being distantly related to Macrogeomys and Heterogeomys, and Macrogeomys appearing as the most derived. Within the subgenus Macrogeomys, O. heterodus and O. cherriei form a distinct clade, as do O. dariensis and O. cavator. As with previous protein-electrophoretic studies, the placement of O. underwoodi could not be determined definitively within the subgenus Macrogeomys. We interpret our inability to determine phylogenetic relationships among these three clades as evidence for a rapid phyletic radiation within this subgenus. Sequence divergence estimates indicate that the Macrogeomys radiation took place following the time of completion of the Panamanian land bridge (1.9-2.9 mya). Additionally, the near identity of sequences of a newly described species, O. thaeleri, with those of O. dariensis (percentage sequence divergence = 0.3%) suggests that the two may be conspecific.     

Molecular systematics of North American phoxinin genera (Actinopterygii: Cyprinidae) inferred from mitochondrial 12S and 16S ribosomal RNA sequences

The cyprinid fish fauna of North America is relatively large, with approximately 300 species, and all but one of these are considered phoxinins. The phylogenetic relationships of the North American phoxinins continue to pose difficulties for systematists. Results of morphological analyses are not consistent owing to differences interpreting and coding characters. Herein, we present phylogenetic analyses of mitochondrial 12S and 16S ribosomal RNA sequence data for representatives of nearly all genera of North American phoxinins. The data were analysed using parsimony, weighted parsimony, maximum likelihood and bayesian analyses. Results from weighted parsimony, likelihood and the bayesian analysis are largely consistent as they all account for differing substitution rates between transitions and transversions. Several major clades within the fauna can be recognized and are strongly supported by all analyses. These include the western clade, creek chub–plagopterin clade and the open posterior myodome clade. The shiner clade is nested in the open posterior myodome clade and is the most species-rich clade of North American phoxinins. Relationships within this clade were not well resolved by our analyses. This may reflect the inability of the mitochondrial RNA genes to resolve recent speciation events or taxon sampling within the shiner clade.  © 2003 The Linnean Society of London, Zoological Journal of the Linnean Society , 2003, 139 , 63–80.