ALLOZYME DIVERGENCE WITHIN THE BOVIDAE

We describe a phylogeny of the Bovidae based on 40 allozyme loci in 27 species, representing 10 of the 14 bovid tribes described by Vrba (1985). Giraffe represented a related family (Giraffidae). A phenogram was derived using the unweighted pair-group method with arithmetic means (UPGMA), based on Nei's genetic distances (ND) between species. A tree was also derived using the neighbor-joining technique, also based on ND. To provide a cladistic interpretation, the data were analyzed by a maximum parsimony method (phylogenetic analysis using parsimony, PAUP). We found marked divergence within the Bovidae, consistent with the appearance of the family in the early Miocene. Unexpectedly, the most divergent species was the impala, which occupied a basal position in all trees. Species in the tribe Alcelaphini were the most derived taxa in all trees. These patterns conflict strongly with the previous taxonomic alliance, based on immuno-distance and anatomical evidence, of the impala as a sister group of the Alcelaphini. All trees agreed that tribes described by Vrba (1985) are monophyletic, except the Neotragini, which was polyphyletic, with suni occupying a long branch by itself. The dikdik and klipspringer were consistently placed as sister taxa to species in the Antilopini. Three tribes (Aepycerotini, Tragelaphini and Cephalophini), whose fossils have not been found outside Africa, were basal in all trees, suggesting that bovids originated in Africa. Nodes connecting the remaining tribes were closely clustered, a pattern that agrees with fossil evidence of rapid divergence within the Bovidae in the mid-Miocene (about 15 mybp). The allozyme data suggested a second phase of rapid divergence within tribes during the Plio-Pleistocene, a pattern that also agrees with fossil evidence. Rates of bovid divergence have therefore been far from constant. However, the clustering of nodes imparts considerable uncertainty to the branching order leading to the derived tribes, and to a lesser extent, species within tribes. The classical division of the Bovidae into the Boodontia and Aegeodontia does not agree with the phylogenetic grouping of tribes presented in this analysis. However, the maximum parsimony tree derived using ‘local’ branch swapping clustered all grazing species into a derived, monophyletic group, suggesting that grazing may have evolved only once in bovid evolution.     

Multiple data sets, high homoplasy, and the phylogeny of softshell turtles (Testudines: Trionychidae)

We present a phylogenetic hypothesis and novel, rank-free classification for all extant species of softshell turtles (Testudines:Trionychidae). Our data set included DNA sequence data from two mitochondrial protein-coding genes and a approximately 1-kb nuclear intron for 23 of 26 recognized species, and 59 previously published morphological characters for a complimentary set of 24 species. The combined data set provided complete taxonomic coverage for this globally distributed clade of turtles, with incomplete data for a few taxa. Although our taxonomic sampling is complete, most of the modern taxa are representatives of old and very divergent lineages. Thus, due to biological realities, our sampling consists of one or a few representatives of several ancient lineages across a relatively deep phylogenetic tree. Our analyses of the combined data set converge on a set of well-supported relationships, which is in accord with many aspects of traditional softshell systematics including the monophyly of the Cyclanorbinae and Trionychinae. However, our results conflict with other aspects of current taxonomy and indicate that most of the currently recognized tribes are not monophyletic. We use this strong estimate of the phylogeny of softshell turtles for two purposes: (1) as the basis for a novel rank-free classification, and (2) to retrospectively examine strategies for analyzing highly homoplasious mtDNA data in deep phylogenetic problems where increased taxon sampling is not an option. Weeded and weighted parsimony, and model-based techniques, generally improved the phylogenetic performance of highly homoplasious mtDNA sequences, but no single strategy completely mitigated the problems of associated with these highly homoplasious data. Many deep nodes in the softshell turtle phylogeny were confidently recovered only after the addition of largely nonhomoplasious data from the nuclear intron.     

A Cladistic Analysis of Mitochondrial Ribosomal DNA from the Bovidae

There is a huge data base of genetic information for the domestic artiodactyl speciesBos taurus(cow),Ovis aries(sheep), andCapra hircus(goat). However, the phylogenetic relationships of these economically critical taxa and their close relatives, family Bovidae, remain for the most part unresolved. In this report, we aligned new mitochondrial (mt) 12S and 16S ribosomal (r) DNA sequences from 26 bovid taxa with published sequences. Phylogenetic analyses of the more than 64 kilobases of mt rDNA from 57 taxa support a basal division in the Bovidae that separatesBosand its close relatives fromCapra, Ovis,and their kin. As suggested by previous molecular and morphological studies, “antelopes” are a paraphyletic assemblage. Caprinae (sheep, goats, goat antelopes, and musk oxen) groups consistently with hippotragine and alcelaphine antelopes, while Bovini (cattle and buffaloes) clusters with tragelaphine and boselaphine antelopes. The traditional tribal subdivisions of Bovidae are supported in most cases, but there are exceptions within Caprinae and Antilopinae (gazelles and close relatives). The rDNA data consistently place the enigmatic generaPelea, Pantholops,andSaiga,but the origin ofAepyceros,the impala, remains obscure. Combined phylogenetic analyses of the rDNA data with the skeletal characters of Gentry (1992) were used to assess the stability of the molecular results.     

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.     

Diversification of the North American shrub genus Ceanothus (Rhamnaceae): conflicting phylogenies from nuclear ribosomal DNA and chloroplast DNA

Ceanothus comprises ～55 morphologically and ecologically diverse species of woody perennials endemic to North America. Interpretations of the natural history of Ceanothus have served as a general model of evolution for woody perennials with simple entomophilous pollination systems, but these interpretations lacked explicit phylogenetic context. We used cladistic analysis of sequences of the chloroplast-encoded matK and the internal transcribed spacers (ITS) and 5.8S coding region of nuclear ribosomal DNA (nrDNA) to reconstruct the phylogeny of Ceanothus. The nuclear and organellar phylogenies exhibited very low levels of both topological and character congruence. Subgenera Ceanothus and Cerastes are monophyletic sister taxa in both phylogenies, but both data sets suffer from a lack of resolution below the level of subgenus. Lack of taxonomic congruence between the two data sets may be a result of introgression and/or lineage sorting. The ITS tree was accepted as the better estimate of a species phylogeny for Ceanothus, on the assumption that nuclear markers are less prone to introgression. Three of five polytypic species in the ITS data set were paraphyletic, and four of six polytypic species in the matK data set were paraphyletic. This study demonstrates the degree to which matched independent data sets can produce conflicting summaries of evolutionary history.     

SuperTRI: A new approach based on branch support analyses of multiple independent data sets for assessing reliability of phylogenetic inferences

Supermatrix and supertree are two methods for constructing a phylogenetic tree by using multiple data sets. However, these methods are not a panacea, as conflicting signals between data sets can lead to misinterpret the evolutionary history of taxa. In particular, the supermatrix approach is expected to be misleading if the species-tree signal is not dominant after the combination of the data sets. Moreover, most current supertree methods suffer from two limitations: (i) they ignore or misinterpret secondary (non-dominant) phylogenetic signals of the different data sets; and (ii) the logical basis of node robustness measures is unclear.To overcome these limitations, we propose a new approach, called SuperTRI, which is based on the branch support analyses of the independent data sets, and where the reliability of the nodes is assessed using three measures: the supertree Bootstrap percentage and two other values calculated from the separate analyses: the mean branch support (mean Bootstrap percentage or mean posterior probability) and the reproducibility index.The SuperTRI approach is tested on a data matrix including seven genes for 82 taxa of the family Bovidae (Mammalia, Ruminantia), and the results are compared to those found with the supermatrix approach. The phylogenetic analyses of the supermatrix and independent data sets were done using four methods of tree reconstruction: Bayesian inference, maximum likelihood, and unweighted and weighted maximum parsimony. The results indicate, firstly, that the SuperTRI approach shows less sensitivity to the four phylogenetic methods, secondly, that it is more accurate to interpret the relationships among taxa, and thirdly, that interesting conclusions on introgression and radiation can be drawn from the comparisons between SuperTRI and supermatrix analyses. To cite this article: A. Ropiquet et al., C. R. Biologies 332 (2009).     

Effects of Asymmetric Nuclear Introgression,Introgressive Mitochondrial Sweep,and Purifying Selection on Phylogenetic Reconstruction and Divergence Estimates in the Pacific Clade of Locustella Warblers

When isolated but reproductively compatible populations expand geographically and meet, simulations predict asymmetric introgression of neutral loci from a local to invading taxon. Genetic introgression may affect phylogenetic reconstruction by obscuring topology and divergence estimates. We combined phylogenetic analysis of sequences from one mtDNA and 12 nuDNA loci with analysis of gene flow among 5 species of Pacific Locustella warblers to test for presence of genetic introgression and its effects on tree topology and divergence estimates. Our data showed that nuDNA introgression was substantial and asymmetrical among all members of superspecies groups whereas mtDNA showed no introgression except a single species pair where the invader''s mtDNA was swept by mtDNA of the local species. This introgressive sweep of mtDNA had the opposite direction of the nuDNA introgression and resulted in the paraphyly of the local species'' mtDNA haplotypes with respect to those of the invader. Тhe multilocus nuDNA species tree resolved all inter- and intraspecific relationships despite substantial introgression. However, the node ages on the species tree may be underestimated as suggested by the differences in node age estimates based on non-introgressing mtDNA and introgressing nuDNA. In turn, the introgressive sweep and strong purifying selection appear to elongate internal branches in the mtDNA gene tree.     

Physical organization of the 1.709 satellite IV DNA family in Bovini and Tragelaphini tribes of the Bovidae: sequence and chromosomal evolution

Repetitive DNA in the mammalian genome is a valuable record and marker for evolution, providing information about the order and driving forces related to evolutionary events. The evolutionarily young 1.709 satellite IV DNA family is present near the centromeres of many chromosomes in the Bovidae. Here, we isolated 1.709 satellite DNA sequences from five Bovidae species belonging to Bovini: Bos taurus (BTA, cattle), Bos indicus (BIN, zebu), Bubalus bubalis (BBU, water buffalo) and Tragelaphini tribes: Taurotragus oryx (TOR, eland) and Tragelaphus euryceros (TEU, bongo). Its presence in both tribes shows the sequence predates the evolutionary separation of the two tribes (more than 10 million years ago), and primary sequence shows increasing divergence with evolutionary distance. Genome organization (Southern hybridization) and physical distribution (in situ hybridization) revealed differences in the molecular organization of these satellite DNA sequences. The data suggest that the sequences on the sex chromosomes and the autosomes evolve as relatively independent groups, with the repetitive sequences suggesting that Bovini autosomes and the Tragelaphini sex chromosomes represent the more primitive chromosome forms.     

Combining immunological and molecular data to assess phylogenetic relations of some Greek Podarcis species

Most recent molecular studies revealed the phylogeny of Greek Podarcis species, which for years remained elusive, due to discordant data produced from various chromosomal, complement fixation and protein studies. In this report, we analyzed cellular immune responses of spleen-derived lymphocytes from six allopatric Podarcis species encountered in Greece, by assessing two-way mixed lymphocyte reaction (MLR)-induced proliferation. On the basis of stimulation indices (S.I.) as determined from cultures set up from xenogeneic splenocytes coincubated in pairs, we generated a phylogenetic tree, fully consistent with the phylogenetic relationships of Podarcis as determined by parallel analyses based on partial mitochondrial (mt) DNA sequences. Although the exact mechanisms triggering lymphocyte responses in lizard two-way xenogeneic MLR are not fully understood, our results show the potential use of cell-mediated immune responses as an additional approach to mtDNA analysis, for species delimitation within specific lizard taxa.     

Molecular phylogeny of the subfamily Tephritinae (Diptera: Tephritidae) based on mitochondrial 16S rDNA sequences

The phylogeny of the subfamily Tephritinae (Diptera: Tephritidae) was reconstructed from mitochondrial 16S ribosomal RNA gene sequences using 53 species representing 11 currently recognized tribes of the Tephritinae and 10 outgroup species. The minimum evolution and Bayesian trees suggested the following phylogenetic relationships: (1) monophyly of the Tephritinae was strongly supported; (2) a sister group relationship between the Tephritinae and Plioreocepta was supported by the Bayesian tree; (3) the tribes Tephrellini, Myopitini, and Terelliini (excluding Neaspilota) were supported as monophyletic groups; (4) the non-monophyletic nature of the tribes Dithrycini, Eutretini, Noeetini, Tephritini, Cecidocharini, and Xyphosiini; and (5) recognition of 10 putative tribal groups, most of which were supported strongly by the statistical tests of the interior branches. Our results, therefore, convincingly suggest that an extensive rearrangement of the tribal classification of the Tephritinae is necessary. Since our sampling of taxa heavily relied on the current accepted classification, some lineages identified by the present study were severely under-sampled and other possible major lineages of the Tephritinae were probably not even represented in our dataset. We believe that our results provide baseline information for a more rigorous sampling of additional taxa representing all possible major lineages of the subfamily, which is essential for a comprehensive revision of the tephritine tribal classification.     

Phylogenetic position of the saola (Pseudoryx nghetinhensis) inferred from cytogenetic analysis of eleven species of Bovidae

Previous morphological and molecular analyses failed to resolve the phylogenetic position of the critically endangered saola (Pseudoryx nghetinhensis) with respect to its placement in Bovina (cattle, bison, and yak) or Bubalina (Asian and African buffaloes). In the present study, G- and C-banding, Ag-staining and FISH with 28S and telomeric probes was undertaken for 17 bovid species. An analysis of these data allowed us to identify 49 structural rearrangements that included autosomes, gonosomes and 17 different NOR sites. The combined data set was subjected to a cladistic analysis aimed at: (i) providing new insights on phylogenetic relationships of the saola and other species within the subfamily Bovinae, and (ii) testing the suitability of different classes of chromosomal characters for phylogenetic reconstruction of the family Bovidae. The study revealed that nucleolar organizing regions (NORs) are phylogenetically informative. It was shown that at least one, or sometimes two of these characters punctuate divergences that include nodes that are the most basal in the tree, to those that are the most recent. In this context, the shared presence of three NORs in saola and species of Syncerus and Bubalus strongly suggests the saola's placement within the subtribe Bubalina. This contrasts with Robertsonian rearrangements which are informative only at the generic level. These findings suggest that NORs are an important and frequently overlooked source of additional phylogenetic information within the Bovidae that may also have applicability at higher taxonomic levels, possibly even for Pecora.     

The major opsin in bees (Insecta: Hymenoptera): A promising nuclear gene for higher level phylogenetics.

We report the phylogenetic utility of the nuclear gene encoding the long-wavelength opsin (LW Rh) for tribes of bees. Aligned nucleotide sequences were examined in multiple taxa from the four tribes comprising the corbiculate bees within the subfamily Apinae. Phylogenetic analyses of sequence variation in a 502-bp fragment (approx 40% of the coding region) strongly supported the monophyly of each of the four tribes, which are well established from previous studies of morphology and DNA. Trees estimated from parsimony and maximum likelihood analyses of LW Rh sequences show a strongly supported relationship between the tribes Meliponini and Bombini, a relationship that has been found uniformly in studies of other genes (28S, 16S, and cytochrome b). All of the tribal clades as well as relationships among the tribes are supported by high bootstrap values, suggesting the utility of LW Rh in estimating tribal and subfamily rank for these bees. The sequences exhibit minimal base composition bias. Both 1st + 2nd and 3rd position sites provide information for estimating a reliable tree topology. These results suggest that LW Rh, which has not been reported previously in studies of organismal phylogenetics, could provide important new data from the nuclear genome for phylogeny reconstruction.     

Phylogeny of the Bovidae (Artiodactyla, Mammalia), based on mitochondrial ribosomal DNA sequences.

Portions of the 12S and 16S mitochondrial ribosomal genes for 16 species representing nine tribes in the mammal family Bovidae were compared with six previously published orthologous sequences. Phylogenetic analysis of variable nucleotide positions under different constraints and weighting schemes revealed no robust groupings among tribes. Consensus trees support previous hypotheses of monophyly for four clades, including the traditional subfamily Bovinae. However, the basal diversification of bovid tribes, which was largely unresolved by morphological, immunodiffusion, allozyme, and protein sequence data, remains unresolved with the addition of DNA sequence data. The intractability of this systematic problem is consistent with a rapid radiation of the major bovid groups. Several analyses of our data show that monophyly of the Bovidae, which was weakly supported by previous morphological and molecular work, is questionable.     

Phylogenetic evidence of mitochondrial DNA introgression among pocket gophers in New Mexico (family Geomyidae)

Mitochondrial DNA (mtDNA) variation in the cytochrome b gene was determined for two divergent taxa of pocket gophers, Thomomys bottae actuosus and T. b. ruidosae . These two taxa hybridize in a narrow contact zone, but introgression of nuclear markers such as allozymes or chromosomes does not extend much beyond the hybrid zone (Patton et al. 1979). We found that despite their distinctness, the two subspecies shared very similar mtDNA haplotypes. By a comparison of phylogenetic histories derived from nuclear markers (allozymes) and from mtDNA haplotypes sampled in different populations of T. bottae from New Mexico, we show that apparent similarity is due to an introgression of T. b. ruidosae mtDNA into T. b. actuosus nuclear background. Evidence of introgression is not limited to the present-day contact zone between these two taxa, but extends at least 75 km away from it. The actuosus haplotype coexists along with the ruidosae mtDNA in the Gallinas Mts., which are inhabited by otherwise pure T. b. actuosus , while further north only typical actuosus haplotypes were detected. Of several potential mechanisms which could lead to such a geographical pattern of variation, we argue that a combination of range shifts due to climatic fluctuations, and genetic drift are most likely. Horizontal gene transfers due to hybridization are historical events which seem rather common among pocket gophers. Although they can be identified with careful phylogenetic study using independent data sets, the potential for misinterpreting a gene tree as an organismal tree is great in this and other groups of animals.     

The tribal radiation of the family Bovidae (Artiodactyla) and the evolution of the mitochondrial cytochrome b gene

The nucleotide sequence of the complete mitochondrial cytochrome b gene has been determined and compared for 51 species of the family Bovidae and 10 potential pecoran and tragulid outgroups. A detailed saturation analysis at each codon position relative to the maximum parsimony procedure indicates that all transitions on third codon positions do not accumulate in a similar fashion: C-T are more saturated than A-G substitutions. The same trend is observed for second positions but not for first positions where A-G and C-T transitions exhibit roughly the same levels of saturation. Maximum parsimony reconstructions were weighted according to these observations. Maximum parsimony, maximum likelihood, and distance phylogenetic reconstructions all depict a major split within Bovidae. The subfamily Bovinae includes four multifurcating tribes and subtribes: Boselaphini, Tragelaphini, cattle-Bovini (Bos and Bison), and buffalo-Bovini (Bubalus and Syncerus). Its sister group is the subfamily Antilopinae, i.e., all non-Bovinae taxa, represented by seven lineages: Antilopini (including Saiga), Caprini sensu lato (i. e., Caprinae including Pantholops), Hippotragini, Alcelaphini, Reduncini (including Pelea), Aepyceros possibly linked to Neotragus, and Cephalophini possibly linked to Oreotragus (the suni and the klipspringer being members of a polyphyletic Neotragini). These various tribes and major lineages were produced by two noteworthy explosive radiations, which occurred simultaneously between 12.0 and 15.3 MY (Middle Miocene) in the subfamilies Bovinae and Antilopinae.     

A biogeographic–ecological approach to disentangle reticulate evolution in the Triatoma phyllosoma species group (Heteroptera: Triatominae), vectors of Chagas disease

Introgression and incomplete lineage sorting (ILS) are two of the main sources of gene‐tree incongruence; both can confound the assessment of phylogenetic relationships among closely related species. The Triatoma phyllosoma species group is a clade of partially co‐distributed and cross‐fertile Chagas disease vectors. Despite previous efforts, the phylogeny of this group remains unresolved, largely because of substantial gene‐tree incongruence. Here, we sequentially address introgression and ILS to provide a robust phylogenetic hypothesis for the T. phyllosoma species group. To identify likely instances of introgression prior to molecular scrutiny, we assessed biogeographic data and information on fertility of inter‐specific crosses. We first derived a few explicit hybridization hypotheses by considering the degree of spatial overlap within each species pair. Then, we assessed the plausibility of these hypotheses in the light of each species pair's cross‐fertility. Using this contextual information, we evaluated mito‐nuclear (cyt b, ITS‐2) gene‐tree incongruence and found evidence suggesting introgression within two species pairs. Finally, we modeled ILS using a Bayesian multispecies coalescent approach and either (a) a “complete” dataset with all the specimens in our sample, or (b) a “filtered” dataset without putatively introgressed specimens. The “filtered tree” had higher posterior‐probability support, as well as more plausible topology and divergence times, than the “complete tree.” Detecting and filtering out introgression and modeling ILS allowed us to derive an improved phylogenetic hypothesis for the T. phyllosoma species group. Our results illustrate how biogeographic and ecological‐reproductive contextual information can help clarify the systematics and evolution of recently diverged taxa prone to introgression and ILS.     

Phylogenetic analysis of RAD‐seq data: examining the influence of gene genealogy conflict on analysis of concatenated data

One of the major issues in phylogenetic analysis is that gene genealogies from different gene regions may not reflect the true species tree or history of speciation. This has led to considerable debate about whether concatenation of loci is the best approach for phylogenetic analysis. The application of Next‐generation sequencing techniques such as RAD‐seq generates thousands of relatively short sequence reads from across the genomes of the sampled taxa. These data sets are typically concatenated for phylogenetic analysis leading to data sets that contain millions of base pairs per taxon. The influence of gene region conflict among so many loci in determining the phylogenetic relationships among taxa is unclear. We simulated RAD‐seq data by sampling 100 and 500 base pairs from alignments of over 6000 coding regions that each produce one of three highly supported alternative phylogenies of seven species of Drosophila. We conducted phylogenetic analyses on different sets of these regions to vary the sampling of loci with alternative gene trees to examine the effect on detecting the species tree. Irrespective of sequence length sampled per region and which subset of regions was used, phylogenetic analyses of the concatenated data always recovered the species tree. The results suggest that concatenated alignments of Next‐generation data that consist of many short sequences are robust to gene tree/species tree conflict when the goal is to determine the phylogenetic relationships among taxa.     

A molecular phylogeny of the Cephinae (Hymenoptera, Cephidae) based on mtDNA COI gene: a test of traditional classification

Cephinae is traditionally divided into three tribes and about 24 genera based on morphology and host utilization. There has been no study testing the monophyly of taxa under a strict phylogenetic criterion. A molecular phylogeny of Cephinae based on a total of 68 sequences of mtDNA COI gene, representing seven genera of Cephinae, is reconstructed to test the traditional limits and relationships of taxa. Monophyly of the traditional tribes is not supported. Monophyly of the genera are largely supported except for Pachycephus. A few host shift events are suggested based on phylogenetic relationships among taxa. These results indicate that a more robust phylogeny is required for a more plausible conclusion. We also report two species of Cephus for the first time from Turkey.     

Brassicaceae phylogeny inferred from phytochrome A and ndhF sequence data: tribes and trichomes revisited

The family Brassicaceae comprises 3710 species in 338 genera, 25 recently delimited tribes, and three major lineages based on phylogenetic results from the chloroplast gene ndhF. To assess the credibility of the lineages and newly delimited tribes, we sequenced an approximately 1.8-kb region of the nuclear phytochrome A (PHYA) gene for taxa previously sampled for the chloroplast gene ndhF. Using parsimony, likelihood, and Bayesian methods, we reconstructed the phylogeny of the gene and used the approximately unbiased (AU) test to compare phylogenetic results from PHYA with findings from ndhF. We also combined ndhF and PHYA data and used a Bayesian mixed model approach to infer phylogeny. PHYA and combined analyses recovered the same three large lineages as those recovered in ndhF trees, increasing confidence in these lineages. The combined tree confirms the monophyly of most of the recently delimited tribes (only Alysseae, Anchonieae, and Descurainieae are not monophyletic), while 13 of the 23 sampled tribes are monophyletic in PHYA trees. In addition to phylogenetic results, we documented the trichome branching morphology of species across the phylogeny and explored the evolution of different trichome morphologies using the AU test. Our results indicate that dendritic, medifixed, and stellate trichomes likely evolved independently several times in the Brassicaceae.