Virulence related sequences; insights provided by comparative genomics of Streptococcus uberis of differing virulence

Background

Streptococcus uberis, a Gram-positive, catalase-negative member of the family Streptococcaceae is an important environmental pathogen responsible for a significant proportion of subclinical and clinical bovine intramammary infections. Currently, the genome of only a single reference strain (0140J) has been described. Here we present a comparative analysis of complete draft genome sequences of an additional twelve S. uberis strains.

Results

Pan and core genome analysis revealed the core genome common to all strains to be 1,550 genes in 1,509 orthologous clusters, complemented by 115-246 accessory genes present in one or more S. uberis strains but absent in the reference strain 0140J. Most of the previously predicted virulent genes were present in the core genome of all 13 strains but gene gain/loss was observed between the isolates in CDS associated with clustered regularly interspaced short palindromic repeats (CRISPRs), prophage and bacteriocin production. Experimental challenge experiments confirmed strain EF20 as non-virulent; only able to infect in a transient manner that did not result in clinical mastitis. Comparison of the genome sequence of EF20 with the validated virulent strain 0140J identified genes associated with virulence, however these did not relate clearly with clinical/non-clinical status of infection.

Conclusion

The gain/loss of mobile genetic elements such as CRISPRs and prophage are a potential driving force for evolutionary change. This first “whole-genome” comparison of strains isolated from clinical vs non-clinical intramammary infections including the type virulent vs non-virulent strains did not identify simple gene gain/loss rules that readily explain, or be confidently associated with, differences in virulence. This suggests that a more complex dynamic determines infection potential and clinical outcome not simply gene content.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1512-6) contains supplementary material, which is available to authorized users.     

Identification of three extra-chromosomal replicons in Leptospira pathogenic strain and development of new shuttle vectors

Background

The genome of pathogenic Leptospira interrogans contains two chromosomes. Plasmids and prophages are known to play specific roles in gene transfer in bacteria and can potentially serve as efficient genetic tools in these organisms. Although plasmids and prophage remnants have recently been reported in Leptospira species, their characteristics and potential applications in leptospiral genetic transformation systems have not been fully evaluated.

Results

Three extrachromosomal replicons designated lcp1 (65,732 bp), lcp2 (56,757 bp), and lcp3 (54,986 bp) in the L. interrogans serovar Linhai strain 56609 were identified through whole genome sequencing. All three replicons were stable outside of the bacterial chromosomes. Phage particles were observed in the culture supernatant of 56609 after mitomycin C induction, and lcp3, which contained phage-related genes, was considered to be an inducible prophage. L. interrogans–Escherichia coli shuttle vectors, constructed with the predicted replication elements of single rep or rep combined with parAB loci from the three plasmids were shown to successfully transform into both saprophytic and pathogenic Leptospira species, suggesting an essential function for rep genes in supporting auto-replication of the plasmids. Additionally, a wide distribution of homologs of the three rep genes was identified in L. interrogans isolates, and correlation tests showed that the transformability of the shuttle vectors in L. interrogans isolates depended, to certain extent, on genetic compatibility between the rep sequences of both plasmid and host.

Conclusions

Three extrachromosomal replicons co-exist in L. interrogans, one of which we consider to be an inducible prophage. The vectors constructed with the rep genes of the three replicons successfully transformed into saprophytic and pathogenic Leptospira species alike, but this was partly dependent on genetic compatibility between the rep sequences of both plasmid and host.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1321-y) contains supplementary material, which is available to authorized users.     

High resolution assembly and characterization of genomes of Canadian isolates of Salmonella Enteritidis

Background

There is a need to characterize genomes of the foodborne pathogen, Salmonella enterica serovar Enteritidis (SE) and identify genetic information that could be ultimately deployed for differentiating strains of the organism, a need that is yet to be addressed mainly because of the high degree of clonality of the organism. In an effort to achieve the first characterization of the genomes of SE of Canadian origin, we carried out massively parallel sequencing of the nucleotide sequence of 11 SE isolates obtained from poultry production environments (n = 9), a clam and a chicken, assembled finished genomes and investigated diversity of the SE genome.

Results

The median genome size was 4,678,683 bp. A total of 4,833 chromosomal genes defined the pan genome of our field SE isolates consisting of 4,600 genes present in all the genomes, i.e., core genome, and 233 genes absent in at least one genome (accessory genome). Genome diversity was demonstrable by the presence of 1,360 loci showing single nucleotide polymorphism (SNP) in the core genome which was used to portray the genetic distances by means of a phylogenetic tree for the SE isolates. The accessory genome consisted mostly of previously identified SE prophage sequences as well as two, apparently full- sized, novel prophages namely a 28 kb sequence provisionally designated as SE-OLF-10058 (3) prophage and a 43 kb sequence provisionally designated as SE-OLF-10012 prophage.

Conclusions

The number of SNPs identified in the relatively large core genome of SE is a reflection of substantial diversity that could be exploited for strain differentiation as shown by the development of an informative phylogenetic tree. Prophage sequences can also be exploited for SE strain differentiation and lineage tracking. This work has laid the ground work for further studies to develop a readily adoptable laboratory test for the subtyping of SE.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-713) contains supplementary material, which is available to authorized users.     

Prophage-like elements present in Mycobacterium genomes

Background

Prophages, integral components of many bacterial genomes, play significant roles in cognate host bacteria, such as virulence, toxin biosynthesis and secretion, fitness cost, genomic variations, and evolution. Many prophages and prophage-like elements present in sequenced bacterial genomes, such as Bifidobacteria, Lactococcus and Streptococcus, have been described. However, information for the prophage of Mycobacterium remains poorly defined.

Results

In this study, based on the search of the complete genome database from GenBank, the Whole Genome Shotgun (WGS) databases, and some published literatures, thirty-three prophages were described in detail. Eleven of them were full-length prophages, and others were prophage-like elements. Eleven prophages were firstly revealed. They were phiMAV_1, phiMAV_2, phiMmcs_1, phiMmcs_2, phiMkms_1, phiMkms_2, phiBN42_1, phiBN44_1, phiMCAN_1, phiMycsm_1, and phiW7S_1. Their genomes and gene contents were firstly analyzed. Furthermore, comparative genomics analyses among mycobacterioprophages showed that full-length prophage phi172_2 belonged to mycobacteriophage Cluster A and the phiMmcs_1, phiMkms_1, phiBN44_1, and phiMCAN_1 shared high homology and could be classified into one group.

Conclusions

To our knowledge, this is the first systematic characterization of mycobacterioprophages, their genomic organization and phylogeny. This information will afford more understanding of the biology of Mycobacterium.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-243) contains supplementary material, which is available to authorized users.     

CRISPR-Cas systems target a diverse collection of invasive mobile genetic elements in human microbiomes

Background

Bacteria and archaea develop immunity against invading genomes by incorporating pieces of the invaders'' sequences, called spacers, into a clustered regularly interspaced short palindromic repeats (CRISPR) locus between repeats, forming arrays of repeat-spacer units. When spacers are expressed, they direct CRISPR-associated (Cas) proteins to silence complementary invading DNA. In order to characterize the invaders of human microbiomes, we use spacers from CRISPR arrays that we had previously assembled from shotgun metagenomic datasets, and identify contigs that contain these spacers'' targets.

Results

We discover 95,000 contigs that are putative invasive mobile genetic elements, some targeted by hundreds of CRISPR spacers. We find that oral sites in healthy human populations have a much greater variety of mobile genetic elements than stool samples. Mobile genetic elements carry genes encoding diverse functions: only 7% of the mobile genetic elements are similar to known phages or plasmids, although a much greater proportion contain phage- or plasmid-related genes. A small number of contigs share similarity with known integrative and conjugative elements, providing the first examples of CRISPR defenses against this class of element. We provide detailed analyses of a few large mobile genetic elements of various types, and a relative abundance analysis of mobile genetic elements and putative hosts, exploring the dynamic activities of mobile genetic elements in human microbiomes. A joint analysis of mobile genetic elements and CRISPRs shows that protospacer-adjacent motifs drive their interaction network; however, some CRISPR-Cas systems target mobile genetic elements lacking motifs.

Conclusions

We identify a large collection of invasive mobile genetic elements in human microbiomes, an important resource for further study of the interaction between the CRISPR-Cas immune system and invaders.     

Reticulate evolution in North American black-fruited hawthorns (Crataegus section Douglasia; Rosaceae): evidence from nuclear ITS2 and plastid sequences

Background and Aims

The taxonomic complexity of Crataegus (hawthorn; Rosaceae, Maleae), especially in North America, has been attributed by some to hybridization in combination with gametophytic apomixis and polyploidization, whereas others have considered the roles of hybridization and apomixis to be minimal. Study of the chemical composition and therapeutic value of hawthorn extracts requires reproducible differentiation of entities that may be difficult to distinguish by morphology alone. This study sought to address this by using the nuclear ribosomal spacer region ITS2 as a supplementary DNA barcode; however, a lack of success prompted an investigation to discover why this locus gave unsatisfactory results.

Methods

ITS2 was extensively cloned so as to document inter- and intraindividual variation in this locus, using hawthorns of western North America where the genus Crataegus is represented by only two widely divergent groups, the red-fruited section Coccineae and the black-fruited section Douglasia. Additional sequence data from selected loci on the plastid genome were obtained to enhance further the interpretation of the ITS2 results.

Key Results

In the ITS2 gene tree, ribotypes from western North American hawthorns are found in two clades. Ribotypes from diploid members of section Douglasia occur in one clade (with representatives of the east-Asian section Sanguineae). The other clade comprises those from diploid and polyploid members of section Coccineae. Both clades contribute ribotypes to polyploid Douglasia. Data from four plastid-derived intergenic spacers demonstrate the maternal parentage of these allopolyploids.

Conclusions

Repeated hybridization between species of section Douglasia and western North American members of section Coccineae involving the fertilization of unreduced female gametes explains the observed distribution of ribotypes and accounts for the phenetic intermediacy of many members of section Douglasia.     

Identification of prophages in bacterial genomes by dinucleotide relative abundance difference

Background

Prophages are integrated viral forms in bacterial genomes that have been found to contribute to interstrain genetic variability. Many virulence-associated genes are reported to be prophage encoded. Present computational methods to detect prophages are either by identifying possible essential proteins such as integrases or by an extension of this technique, which involves identifying a region containing proteins similar to those occurring in prophages. These methods suffer due to the problem of low sequence similarity at the protein level, which suggests that a nucleotide based approach could be useful.

Methodology

Earlier dinucleotide relative abundance (DRA) have been used to identify regions, which deviate from the neighborhood areas, in genomes. We have used the difference in the dinucleotide relative abundance (DRAD) between the bacterial and prophage DNA to aid location of DNA stretches that could be of prophage origin in bacterial genomes. Prophage sequences which deviate from bacterial regions in their dinucleotide frequencies are detected by scanning bacterial genome sequences. The method was validated using a subset of genomes with prophage data from literature reports. A web interface for prophage scan based on this method is available at http://bicmku.in:8082/prophagedb/dra.html. Two hundred bacterial genomes which do not have annotated prophages have been scanned for prophage regions using this method.

Conclusions

The relative dinucleotide distribution difference helps detect prophage regions in genome sequences. The usefulness of this method is seen in the identification of 461 highly probable loci pertaining to prophages which have not been annotated so earlier. This work emphasizes the need to extend the efforts to detect and annotate prophage elements in genome sequences.     

Single-molecule sequencing reveals the molecular basis of multidrug-resistance in ST772 methicillin-resistant Staphylococcus aureus

Background

Methicillin-resistant Staphylococcus aureus (MRSA) is a major cause of hospital-associated infection, but there is growing awareness of the emergence of multidrug-resistant lineages in community settings around the world. One such lineage is ST772-MRSA-V, which has disseminated globally and is increasingly prevalent in India. Here, we present the complete genome sequence of DAR4145, a strain of the ST772-MRSA-V lineage from India, and investigate its genomic characteristics in regards to antibiotic resistance and virulence factors.

Results

Sequencing using single-molecule real-time technology resulted in the assembly of a single continuous chromosomal sequence, which was error-corrected, annotated and compared to nine draft genome assemblies of ST772-MRSA-V from Australia, Malaysia and India. We discovered numerous and redundant resistance genes associated with mobile genetic elements (MGEs) and known core genome mutations that explain the highly antibiotic resistant phenotype of DAR4145. Staphylococcal toxins and superantigens, including the leukotoxin Panton-Valentinin Leukocidin, were predominantly associated with genomic islands and the phage φ-IND772PVL. Some of these mobile resistance and virulence factors were variably present in other strains of the ST772-MRSA-V lineage.

Conclusions

The genomic characteristics presented here emphasize the contribution of MGEs to the emergence of multidrug-resistant and highly virulent strains of community-associated MRSA. Antibiotic resistance was further augmented by chromosomal mutations and redundancy of resistance genes. The complete genome of DAR4145 provides a valuable resource for future investigations into the global dissemination and phylogeography of ST772-MRSA-V.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1599-9) contains supplementary material, which is available to authorized users.     

Five vicariant genera from Gondwana: the Velloziaceae as shown by molecules and morphology [corrected

Background and Aims

The amount of data collected previously for Velloziaceae neither clarified relationships within the family nor helped determine an appropriate classification, which has led to huge discordance among treatment by different authors. To achieve an acceptable phylogenetic result and understand the evolution and roles of characters in supporting groups, a total evidence analysis was developed which included approx. 20 % of the species and all recognized genera and sections of Velloziaceae, plus outgroups representatives of related families within Pandanales.

Methods

Analyses were undertaken with 48 species of Velloziaceae, representing all ten genera, with DNA sequences from the atpB-rbcL spacer, trnL-trnF spacer, trnL intron, trnH-psbA spacer, ITS ribosomal DNA spacers and morphology.

Key Results

Four groups consistently emerge from the analyses. Persistent leaves, two phloem strands, stem cortex divided in three regions and violet tepals support Acanthochlamys as sister to Velloziaceae s.s., which are supported mainly by leaves with marginal bundles, transfusion tracheids and inflorescence without axis. Within Velloziaceae s.s., an African Xerophyta + Talbotia clade is uniquely supported by basal loculicidal capsules; an American clade, Barbacenia s.l. + Barbaceniopsis + Nanuza + Vellozia, is supported by only homoplastic characters. Barbacenia s.l. (= Aylthonia + Barbacenia + Burlemarxia + Pleurostima) is supported by a double sheath in leaf vascular bundles and a corona; Barbaceniopsis + Nanuza + Vellozia is not supported by an unambiguous character, but Barbaceniopsis is supported by five characters, including diclinous flowers, Nanuza + Vellozia is supported mainly by horizontal stigma lobes and stem inner cortex cells with secondary walls, and Vellozia alone is supported mainly by pollen in tetrads.

Conclusions

The results imply recognition of five genera (Acanthochlamys (Xerophyta (Barbacenia (Barbaceniopsis, Vellozia)))), solving the long-standing controversies among recent classifications of the family. They also suggest a Gondwanan origin for Velloziaceae, with a vicariant pattern of distribution.     

Why Africa matters: evolution of Old World Salvia (Lamiaceae) in Africa

Background and Aims

Salvia is the largest genus in Lamiaceae and it has recently been found to be non-monophyletic. Molecular data on Old World Salvia are largely lacking. In this study, we present data concerning Salvia in Africa. The focus is on the colonization of the continent, character evolution and the switch of pollination systems in the genus.

Methods

Maximum likelihood and Bayesian inference were used for phylogenetic reconstruction. Analyses were based on two nuclear markers [internal transcribed spacer (ITS) and external transcribed spacer (ETS)] and one plastid marker (rpl32-trnL). Sequence data were generated for 41 of the 62 African taxa (66 %). Mesquite was used to reconstruct ancestral character states for distribution, life form, calyx shape, stamen type and pollination syndrome.

Key Results

Salvia in Africa is non-monophyletic. Each of the five major regions in Africa, except Madagascar, was colonized at least twice, and floristic links between North African, south-west Asian and European species are strongly supported. The large radiation in Sub-Saharan Africa (23 species) can be traced back to dispersal from North Africa via East Africa to the Cape Region. Adaptation to bird pollination in southern Africa and Madagascar reflects parallel evolution.

Conclusions

The phenotypic diversity in African Salvia is associated with repeated introductions to the continent. Many important evolutionary processes, such as colonization, adaptation, parallelism and character transformation, are reflected in this comparatively small group. The data presented in this study can help to understand the evolution of Salvia sensu lato and other large genera.     

Evolutionary history of the Afro-Madagascan Ixora species (Rubiaceae): species diversification and distribution of key morphological traits inferred from dated molecular phylogenetic trees

Background and Aims

Previous work on the pantropical genus Ixora has revealed an Afro-Madagascan clade, but as yet no study has focused in detail on the evolutionary history and morphological trends in this group. Here the evolutionary history of Afro-Madagascan Ixora spp. (a clade of approx. 80 taxa) is investigated and the phylogenetic trees compared with several key morphological traits in taxa occurring in Madagascar.

Methods

Phylogenetic relationships of Afro-Madagascan Ixora are assessed using sequence data from four plastid regions (petD, rps16, rpoB-trnC and trnL-trnF) and nuclear ribosomal external transcribed spacer (ETS) and internal transcribed spacer (ITS) regions. The phylogenetic distribution of key morphological characters is assessed. Bayesian inference (implemented in BEAST) is used to estimate the temporal origin of Ixora based on fossil evidence.

Key Results

Two separate lineages of Madagascan taxa are recovered, one of which is nested in a group of East African taxa. Divergence in Ixora is estimated to have commenced during the mid Miocene, with extensive cladogenesis occurring in the Afro-Madagascan clade during the Pliocene onwards.

Conclusions

Both lineages of Madagascan Ixora exhibit morphological innovations that are rare throughout the rest of the genus, including a trend towards pauciflorous inflorescences and a trend towards extreme corolla tube length, suggesting that the same ecological and selective pressures are acting upon taxa from both Madagascan lineages. Novel ecological opportunities resulting from climate-induced habitat fragmentation and corolla tube length diversification are likely to have facilitated species radiation on Madagascar.     

Multispacer sequence typing relapsing fever Borreliae in Africa

Background

In Africa, relapsing fevers are neglected arthropod-borne infections caused by closely related Borrelia species. They cause mild to deadly undifferentiated fever particularly severe in pregnant women. Lack of a tool to genotype these Borrelia organisms limits knowledge regarding their reservoirs and their epidemiology.

Methodology/Principal Findings

Genome sequence analysis of Borrelia crocidurae, Borrelia duttonii and Borrelia recurrentis yielded 5 intergenic spacers scattered between 10 chromosomal genes that were incorporated into a multispacer sequence typing (MST) approach. Sequencing these spacers directly from human blood specimens previously found to be infected by B. recurrentis (30 specimens), B. duttonii (17 specimens) and B. crocidurae (13 specimens) resolved these 60 strains and the 3 type strains into 13 species-specific spacer types in the presence of negative controls. B. crocidurae comprised of 8 spacer types, B. duttonii of 3 spacer types and B. recurrentis of 2 spacer types.

Conclusions/Significance

Phylogenetic analyses of MST data suggested that B. duttonii, B. crocidurae and B. recurrentis are variants of a unique ancestral Borrelia species. MST proved to be a suitable approach for identifying and genotyping relapsing fever borreliae in Africa. It could be applied to both vectors and clinical specimens.     

Comparative genomic analysis of clinical and environmental strains provides insight into the pathogenicity and evolution of Vibrio parahaemolyticus

Background

Vibrio parahaemolyticus is a Gram-negative halophilic bacterium. Infections with the bacterium could become systemic and can be life-threatening to immunocompromised individuals. Genome sequences of a few clinical isolates of V. parahaemolyticus are currently available, but the genome dynamics across the species and virulence potential of environmental strains on a genome-scale have not been described before.

Results

Here we present genome sequences of four V. parahaemolyticus clinical strains from stool samples of patients and five environmental strains in Hong Kong. Phylogenomics analysis based on single nucleotide polymorphisms revealed a clear distinction between the clinical and environmental isolates. A new gene cluster belonging to the biofilm associated proteins of V. parahaemolyticus was found in clincial strains. In addition, a novel small genomic island frequently found among clinical isolates was reported. A few environmental strains were found harboring virulence genes and prophage elements, indicating their virulence potential. A unique biphenyl degradation pathway was also reported. A database for V. parahaemolyticus (http://kwanlab.bio.cuhk.edu.hk/vp) was constructed here as a platform to access and analyze genome sequences and annotations of the bacterium.

Conclusions

We have performed a comparative genomics analysis of clinical and environmental strains of V. parahaemolyticus. Our analyses could facilitate understanding of the phylogenetic diversity and niche adaptation of this bacterium.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1135) contains supplementary material, which is available to authorized users.     

Extreme phenotypic variation in Cetraria aculeata (lichenized Ascomycota): adaptation or incidental modification?

Background and Aims

Phenotypic variability is a successful strategy in lichens for colonizing different habitats. Vagrancy has been reported as a specific adaptation for lichens living in steppe habitats around the world. Among the facultatively vagrant species, the cosmopolitan Cetraria aculeata apparently forms extremely modified vagrant thalli in steppe habitats of Central Spain. The aim of this study was to investigate whether these changes are phenotypic plasticity (a single genotype producing different phenotypes), by characterizing the anatomical and ultrastructural changes observed in vagrant morphs, and measuring differences in ecophysiological performance.

Methods

Specimens of vagrant and attached populations of C. aculeata were collected on the steppes of Central Spain. The fungal internal transcribed spacer (ITS), glyceraldehyde-3-phosphate dehydrogenase (GPD) and the large sub-unit of the mitochondrial ribosomal DNA (mtLSUm), and the algal ITS and actin were studied within a population genetics framework. Semi-thin and ultrathin sections were analysed by means of optical, scanning electron and transmission electron microscopy. Gas exchange and chlorophyll fluorescence were used to compare the physiological performance of both morphs.

Key Results and Conclusions

Vagrant and attached morphs share multilocus haplotypes which may indicate that they belong to the same species in spite of their completely different anatomy. However, differentiation tests suggested that vagrant specimens do not represent a random sub-set of the surrounding population. The morphological differences were related to anatomical and ultrastructural differences. Large intercalary growth rates of thalli after the loss of the basal–apical thallus polarity may be the cause of the increased growth shown by vagrant specimens. The anatomical and morphological changes lead to greater duration of ecophysiological activity in vagrant specimens. Although the anatomical and physiological changes could be chance effects, the genetic differentiation between vagrant and attached sub-populations and the higher biomass of the former show fitness effects and adaptation to dry environmental conditions in steppe habitats.     

Plastid DNA sequencing and nuclear SNP genotyping help resolve the puzzle of central American Platanus

Background and Aims

Recent research on the history of Platanus reveals that hybridization phenomena occurred in the central American species. This study has two goals: to help resolve the evolutive puzzle of central American Platanus, and to test the potential of real-time polymerase chain reaction (PCR) for detecting ancient hybridization.

Methods

Sequencing of a uniparental plastid DNA marker [psbA-trnH^(GUG) intergenic spacer] and qualitative and quantitative single nucleotide polymorphism (SNP) genotyping of biparental nuclear ribosomal DNA (nrDNA) markers [LEAFY intron 2 (LFY-i2) and internal transcribed spacer 2 (ITS2)] were used.

Key Results

Based on the SNP genotyping results, several Platanus accessions show the presence of hybridization/introgression, including some accessions of P. rzedowskii and of P. mexicana var. interior and one of P. mexicana var. mexicana from Oaxaca (= P. oaxacana). Based on haplotype analyses of the psbA-trnH spacer, five haplotypes were detected. The most common of these is present in taxa belonging to P. orientalis, P. racemosa sensu lato, some accessions of P. occidentalis sensu stricto (s.s.) from Texas, P. occidentalis var. palmeri, P. mexicana s.s. and P. rzedowskii. This is highly relevant to genetic relationships with the haplotypes present in P. occidentalis s.s. and P. mexicana var. interior.

Conclusions

Hybridization and introgression events between lineages ancestral to modern central and eastern North American Platanus species occurred. Plastid haplotypes and qualitative and quantitative SNP genotyping provide information critical for understanding the complex history of Mexican Platanus. Compared with the usual molecular techniques of sub-cloning, sequencing and genotyping, real-time PCR assay is a quick and sensitive technique for analysing complex evolutionary patterns.     

Lessons from Plectocephalus (Compositae, Cardueae-Centaureinae): ITS disorientation in annuals and Beringian dispersal as revealed by molecular analyses

Background and Aims

The geographic distribution of the genus Plectocephalus comprises a single species in Ethiopia, two in North America and possibly four more in South America, in a striking disjunction that is exceptional for genera of the tribe Cardueae. The enormity of this disjunction cast doubts on the precise taxonomic delineation of the genus, which is not unanimously recognized as a natural entity. The aims of this study were to define the generic boundaries of Plectocephalus and to formulate a hypothesis that would explain its natural range.

Methods

A combined molecular approach, using nuclear internal transcribed spacers (ITS) and external transcribed spacers (ETS), and plastid trnL-trnL-F, rpl32-trnL^UAG and ndhF markers, was chosen for phylogenetic reconstruction by maximum parsimony and Bayesian inference.

Key Results

Phylogenetic analysis shows that Plectocephalus is a natural genus that includes the African species P. varians, together with all the native South American species, currently classified as Centaurea, C. cachinalensis, C. floccosa and C. tweediei. The recognition of Centaurodendron as an independent genus, which we consider appropriate, would make Plectocephalus paraphyletic. Affinities of Plectocephalus should lie with eastern representatives of Centaureinae. Geographic disjunction is explained as a consequence of dispersal via the Bering Land Bridge during the Miocene–Pliocene. The phylogeny of the basal grade of Centaureinae differs from previous phylogenies, and artefacts resulting from differences in mutation rates of annual and perennial taxa are confirmed. Sensitivity of ITS to these differences was the highest observed for all DNA regions used in this study.

Conclusions

The natural status of the genus Plectocephalus is confirmed and several nomenclatural combinations are proposed. New evidence contributes to the debate concerning problems posed by the use of ITS in the phylogenetic reconstruction of groups that differ in terms of their life cycles. Dispersal from Caucasus and Anatolia along the Siberian route and then across the Bering Land Bridge follows a route previously proposed for other taxonomic groups.