Concerted evolution in the ribosomal RNA cistron

High-resolution melting curve analysis: A detection assay for Ceratocystis eucalypticola and C. manginecans in infected Eucalyptus

Eucalyptus spp. in plantations are negatively affected by canker and wilt diseases caused by several species of Ceratocystis, particularly those in the Latin American Clade (LAC). Ceratocystis eucalypticola and Ceratocystis manginecans are of particular concern where disease epidemics are reported globally, with recent outbreaks emerging in South African and Indonesian Eucalyptus plantations. Consequently, a rapid screening protocol is required for these pathogens. In this study, a high-resolution melting curve analysis (HRMA) was developed to detect C. eucalypticola and C. manginecans that bypasses time-consuming isolation and post-PCR procedures. Primers targeting a 172 bp region of the cerato-platanin (CP) gene were designed. Using these primers, the accuracy of HRMA to detect and distinguish between these two LAC species was assessed using pure fungal DNA, and DNA extracted directly from Eucalyptus samples naturally infected with C. eucalypticola. The assay accurately detected the presence of C. eucalypticola and C. manginecans and quantifies their DNA, both from cultures, and directly from wood samples. HRMA further differentiated these two species from all other tested LAC individuals. This assay was also able to detect the presence of all the tested LAC species and distinguish seven of these, including C. fimbriata, to species level. Ceratocystis polyconidia was the only non-LAC off-target species detected. Based on these results, the developed assay can be used to rapidly identify C. eucalypticola and C. manginecans directly from infected plant material or fungal cultures, with the potential to also screen for several other LAC species.

Anatomy of a pest control failure: introgression of cytochrome P450 337B3 alleles from invasive old-world bollworm into native corn earworm (Lepidoptera: Noctuidae)

The establishment of invasive species populations can threaten the ecological balance in naïve habitats and impact agricultural production practices. Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) (old-world bollworm, OWBW) and Helicoverpa zea (corn earworm, CEW) were geographically separated prior to the 2013 report of OWBW invasion into South America. Introgression of OWBW-specific cytochrome P450 337B3 (CYP337B3) gene into CEW was repeatedly detected across South America and the Caribbean. Two hybrids were documented from Texas in 2019. In this study, screening insects collected in Olathe, CO, USA, where a failure of pyrethroids to control CEW damage to conventional sweetcorn in 2023 detected 28.6% of insects with the OWBW-specific CYP337B3 marker. Nucleotide sequencing of the CYP337B3 gene identified 73.1% and 26.9% of insects carried CYP337B3v2 and CYP337B3v6 alleles, respectively, and 0.15 overall frequency of CYP337B3 alleles. Based on prior data for distinct phylogeographic origins of CYP337B3v2 and v6 alleles, our results indicate Olathe samples were derived from 2 different introductions: An uncertain source of the v6 allele that was initially reported in West Africa and possibly South American or Caribbean origin of the globally distributed v2 allele. One of the 1618 individuals screened also carried a ribosomal RNA internal transcribed spacer 1 derived from OWBW. Local selection pressures at the Olathe location imposed by repeated pyrethroid exposures are likely attributed to the prevalence of CYP337B3, where control practices hasten the accumulation of phenotypic resistance by adaptive introgression. Pyrethroid and other resistance factors carried by invasive OWBW may continue to impact CEW management tactics across the Americas.

A taxonomic revision of the genus Angelica (Apiaceae) in Taiwan with a new species A. aliensis

BACKGROUND

Angelica L. sensu lato is a taxonomically complex genus, and many studies have utilized morphological and molecular features to resolve its classification issues. In Taiwan, there are six taxa within Angelica, and their taxonomic treatments have been a subject of controversy. In this study, we conducted a comprehensive analysis incorporating morphological and molecular (cpDNA and nrDNA) characteristics to revise the taxonomic treatments of Angelica in Taiwan.

RESULTS

As a result of our research, we have revised the classification between A. dahurica var. formosana and A. pubescens and merged two varieties of A. morrisonicola into a single taxon. A new taxon, A. aliensis, has been identified and found to share a close relationship with A. tarokoensis. Based on the morphological and molecular characteristics data, it has been determined that the former three taxa should be grouped into the Eurasian Angelica clade, while the remaining four taxa should belong to the littoral Angelica clade. Furthermore, Angelica species in Taiwan distributed at higher altitudes displayed higher genetic diversity, implying that the central mountain range of Taiwan serves as a significant reservoir of plant biodiversity. Genetic drift, such as bottlenecks, has been identified as a potential factor leading to the fixation or reduction of genetic diversity of populations in most Angelica species. We provide key to taxa, synopsis, phenology, and distribution for each taxon of Taiwan.

CONCLUSIONS

Our comprehensive analysis of morphological and molecular features has shed light on the taxonomic complexities within Angelica in Taiwan, resolving taxonomic issues and providing valuable insights into the phylogenetic relationships of Angelica in Taiwan.

A Custom Regional DNA Barcode Reference Library for Lichen-Forming Fungi of the Intermountain West, USA, Increases Successful Specimen Identification

DNA barcoding approaches provide powerful tools for characterizing fungal diversity. However, DNA barcoding is limited by poor representation of species-level diversity in fungal sequence databases. Can the development of custom, regionally focused DNA reference libraries improve species-level identification rates for lichen-forming fungi? To explore this question, we created a regional ITS database for lichen-forming fungi (LFF) in the Intermountain West of the United States. The custom database comprised over 4800 sequences and represented over 600 formally described and provisional species. Lichen communities were sampled at 11 sites throughout the Intermountain West, and LFF diversity was characterized using high-throughput ITS2 amplicon sequencing. We compared the species-level identification success rates from our bulk community samples using our regional ITS database and the widely used UNITE database. The custom regional database resulted in significantly higher species-level assignments (72.3%) of candidate species than the UNITE database (28.3-34.2%). Within each site, identification of candidate species ranged from 72.3-82.1% using the custom database; and 31.5-55.4% using the UNITE database. These results highlight that developing regional databases may accelerate a wide range of LFF research by improving our ability to characterize species-level diversity using DNA barcoding.

Genomic Resources for the First Federally Endangered Lichen: The Florida Perforate Cladonia (Cladonia perforata)

Thirty years after its designation as a federally endangered species, the Florida Perforate Cladonia (FPC) remains imperiled in isolated populations in the Florida scrub in the southeastern USA. For threatened and endangered species, such as FPC, reference genomes provide critical insight into genomic diversity, local adaptations, landscape-level genetics, and phylogenomics. Using high-throughput sequencing, we assemble the first draft nuclear and mitochondrial genomes for the FPC mycobiont-Cladonia perforata. We also assess genetic diversity within and among populations in southeastern Florida using genome-scale data and investigate diversity across the entire nuclear ribosomal cistron, including the standard DNA barcoding marker for fungi. The draft nuclear genome spanned 33.6 Mb, and the complete, circular mitochondrial genome was 59 Kb. We also generated the first chloroplast genome, to our knowledge, for the photobiont genus associated with FPC, an undescribed Asterochloris species. We inferred the presence of multiple, distinct mycobiont parental genotypes (genets) occurring at local scales in southeastern Florida, and strikingly, no genets were shared among even the closest sample sites. All sampled thalli shared identical mitochondrial genomes, while the nuclear ribosomal cistron showed limited variability-highlighting the genetic resolution provided by nuclear genome-scale datasets. The genomic resources generated here provide critical resources for informed conservation efforts for the FPC.

Evidence of Biparental Mitochondrial Inheritance from Self-Fertile Crosses between Closely Related Species of Ceratocystis

Hybridization is recognized as a notable driver of evolution and adaptation, which closely related species may exploit in the form of incomplete reproductive barriers. Three closely related species of Ceratocystis (i.e., C. fimbriata, C. manginecans and C. eucalypticola) have previously been shown to hybridize. In such studies, naturally occurring self-sterile strains were mated with an unusual laboratory-generated sterile isolate type, which could have impacted conclusions regarding the prevalence of hybridization and inheritance of mitochondria. In the current study, we investigated whether interspecific crosses between fertile isolates of these three species are possible and, if so, how mitochondria are inherited by the progeny. For this purpose, a PCR-RFLP method and a mitochondrial DNA-specific PCR technique were custom-made. These were applied in a novel approach of typing complete ascospore drops collected from the fruiting bodies in each cross to distinguish between self-fertilizations and potential hybridization. These markers showed hybridization between C. fimbriata and C. eucalypticola and between C. fimbriata and C. manginecans, while no hybridization was detected in the crosses involving C. manginecans and C. eucalypticola. In both sets of hybrid progeny, we detected biparental inheritance of mitochondria. This study was the first to successfully produce hybrids from a cross involving self-fertile isolates of Ceratocystis and also provided the first direct evidence of biparental mitochondrial inheritance in the Ceratocystidaceae. This work lays the foundation for further research focused on investigating the role of hybridization in the speciation of Ceratocystis species and if mitochondrial conflict could have influenced the process.

Structure-Based GC Investigation Sheds New Light on ITS2 Evolution in Corydalis Species

Guanine and cytosine (GC) content is a fundamental component of genetic diversity and essential for phylogenetic analyses. However, the GC content of the ribosomal internal transcribed spacer 2 (ITS2) remains unknown, despite the fact that ITS2 is a widely used phylogenetic marker. Here, the ITS2 was high-throughput sequenced from 29 Corydalis species, and their GC contents were comparatively investigated in the context of ITS2's characteristic secondary structure and concerted evolution. Our results showed that the GC contents of ITS2 were 131% higher than those of their adjacent 5.8S regions, suggesting that ITS2 underwent GC-biased evolution. These GCs were distributed in a heterogeneous manner in the ITS2 secondary structure, with the paired regions being 130% larger than the unpaired regions, indicating that GC is chosen for thermodynamic stability. In addition, species with homogeneous ITS2 sequences were always GC-rich, supporting GC-biased gene conversion (gBGC), which occurred with ITS2's concerted evolution. The RNA substitution model inferred also showed a GC preference among base pair transformations, which again supports gBGC. Overall, structurally based GC investigation reveals that ITS2 evolves under structural stability and gBGC selection, significantly increasing its GC content.

Galactomyces candidus diversity in the complex mycobiota of cow-milk bryndza cheese comprising antagonistic and sensitive strains

Galactomyces candidus (orthographic variant: Galactomyces candidum) is a heterogeneous species of Saccharomycetales that comprises dimorphic yeasts described previously under various names (e.g. Geotrichum, Dipodascus). Its strains are common components of the cheese surface mycobiota. This study identified genetically and physiologically heterogeneous G. candidus strains in the complex mycobiota of artisanal cow-milk bryndza samples. The traditional Slovak bryndza is a cheese produced from ewe's milk in cooler mountainous regions and from cow's milk in warmer low-land regions. The taxonomic analysis of the culturable yeasts of the latter version carried out in this study revealed considerable differences from the yeast biota previously described for ovine bryndza. However, the conventional D1/D2- and ITS-based barcode analyses could not assign unanimously all isolates to species because of the intragenomic barcode diversity in certain groups and the discordance between the D1/D2 and ITS results in other groups. The identified species and groups of isolates had different abilities to utilise the carbon and energy sources (lactose, lactate, lipids and proteins) available in milk and ripening cheese. The G. candidus strains did not metabolise lactose and lactate, hydrolysed milk proteins with diverse, usually moderate efficiency and only could grow on certain amino acids as only energy sources. Their preferred substrate was lipid. Under aerobic conditions, its hyphae penetrated the lipid droplets and degraded their content from inside by developing a dense internal mycelium. Sporulation and different MLST (multilocus sequence typing) patterns indicated that the Galactomyces strains could sexually interact and their genomes could recombine. The Galactomyces and Kluyveromyces isolates had antagonistic effects against other members of the mycobiota.

Beyond Nuclear Ribosomal DNA Sequences: Evolution, Taxonomy, and Closest Known Saprobic Relatives of Powdery Mildew Fungi (Erysiphaceae) Inferred From Their First Comprehensive Genome-Scale Phylogenetic Analyses

Powdery mildew fungi (Erysiphaceae), common obligate biotrophic pathogens of many plants, including important agricultural and horticultural crops, represent a monophyletic lineage within the Ascomycota. Within the Erysiphaceae, molecular phylogenetic relationships and DNA-based species and genera delimitations were up to now mostly based on nuclear ribosomal DNA (nrDNA) phylogenies. This is the first comprehensive genome-scale phylogenetic analysis of this group using 751 single-copy orthologous sequences extracted from 24 selected powdery mildew genomes and 14 additional genomes from Helotiales, the fungal order that includes the Erysiphaceae. Representative genomes of all powdery mildew species with publicly available whole-genome sequencing (WGS) data that were of sufficient quality were included in the analyses. The 24 powdery mildew genomes included in the analysis represented 17 species belonging to eight out of 19 genera recognized within the Erysiphaceae. The epiphytic genera, all but one represented by multiple genomes, belonged each to distinct, well-supported lineages. Three hemiendophytic genera, each represented by a single genome, together formed the hemiendophytic lineage. Out of the 14 other taxa from the Helotiales, Arachnopeziza araneosa, a saprobic species, was the only taxon that grouped together with the 24 genome-sequenced powdery mildew fungi in a monophyletic clade. The close phylogenetic relationship between the Erysiphaceae and Arachnopeziza was revealed earlier by a phylogenomic study of the Leotiomycetes. Further analyses of powdery mildew and Arachnopeziza genomes may discover signatures of the evolutionary processes that have led to obligate biotrophy from a saprobic way of life. A separate phylogeny was produced using the 18S, 5.8S, and 28S nrDNA sequences of the same set of powdery mildew specimens and compared to the genome-scale phylogeny. The nrDNA phylogeny was largely congruent to the phylogeny produced using 751 orthologs. This part of the study has revealed multiple contamination and other quality issues in some powdery mildew genomes. We recommend that the presence of 28S, internal transcribed spacer (ITS), and 18S nrDNA sequences in powdery mildew WGS datasets that are identical to those determined by Sanger sequencing should be used to assess the quality of assemblies, in addition to the commonly used Benchmarking Universal Single-Copy Orthologs (BUSCO) values.

Population Diversity and Genetic Structure Reveal Patterns of Host Association and Anthropogenic Impact for the Globally Important Fungal Tree Pathogen Ceratocystis manginecans

Species in the Ceratocystis manginecans complex are important fungal pathogens of plantation trees globally. The most important hosts include species of Eucalyptus, Acacia, Mangifera, and Punica. Despite their relevance and widespread occurrence, little is known regarding their population genetics and how this might relate to their host associations or geographic regions in which they occur. A global collection of 491 isolates representing the C. manginecans complex, from four different plant hosts and nine countries, were genotyped using microsatellite markers. Population genetic analyses using numerous tools were conducted to interrogate how their genetic diversity and structure might be affected by host or areas of occurrence. Results of genetic diversity studies showed that when grouping isolates into populations based on their host associations, the population on Eucalyptus was most diverse, and it also has a broad global distribution. When considering countries of origin as a basis for defining populations, the gene and genotypic diversity were highest in populations from China, Indonesia, and Brazil. In contrast, populations from Oman and Pakistan collected from Mangifera had the lowest genetic diversity and were clonal. Molecular variance, population differentiation, and network and structure analyses showed that the genetic structure of isolates in the C. manginecans complex is influenced by both host association as well as geographical isolation. Furthermore, the results reflected the movement of genotypes between plant hosts and geographic regions that have implications regarding the broad global distribution of this pathogen.

The role of internal transcribed spacer 2 secondary structures in classifying mycoparasitic Ampelomyces

Many fungi require specific growth conditions before they can be identified. Direct environmental DNA sequencing is advantageous, although for some taxa, specific primers need to be used for successful amplification of molecular markers. The internal transcribed spacer region is the preferred DNA barcode for fungi. However, inter- and intra-specific distances in ITS sequences highly vary among some fungal groups; consequently, it is not a solely reliable tool for species delineation. Ampelomyces, mycoparasites of the fungal phytopathogen order Erysiphales, can have ITS genetic differences up to 15%; this may lead to misidentification with other closely related unknown fungi. Indeed, Ampelomyces were initially misidentified as other pycnidial mycoparasites, but subsequent research showed that they differ in pycnidia morphology and culture characteristics. We investigated whether the ITS2 nucleotide content and secondary structure was different between Ampelomyces ITS2 sequences and those unrelated to this genus. To this end, we retrieved all ITS sequences referred to as Ampelomyces from the GenBank database. This analysis revealed that fungal ITS environmental DNA sequences are still being deposited in the database under the name Ampelomyces, but they do not belong to this genus. We also detected variations in the conserved hybridization model of the ITS2 proximal 5.8S and 28S stem from two Ampelomyces strains. Moreover, we suggested for the first time that pseudogenes form in the ITS region of this mycoparasite. A phylogenetic analysis based on ITS2 sequences-structures grouped the environmental sequences of putative Ampelomyces into a different clade from the Ampelomyces-containing clades. Indeed, when conducting ITS2 analysis, resolution of genetic distances between Ampelomyces and those putative Ampelomyces improved. Each clade represented a distinct consensus ITS2 S2, which suggested that different pre-ribosomal RNA (pre-rRNA) processes occur across different lineages. This study recommends the use of ITS2 S2s as an important tool to analyse environmental sequencing and unveiling the underlying evolutionary processes.

Identification and Characterization of Ceratocystis fimbriata Causing Lethal Wilt on the Lansium Tree in Indonesia

Bark canker, wood discoloration, and wilting of the duku tree (Lansium domesticum) along the watershed of Komering River, South Sumatra Province, Indonesia first appeared in 2013. The incidence of tree mortality was 100% within 3 years in badly infected orchards. A Ceratocystis species was consistently isolated from the diseased tissue and identified by morphological and sequence analyses of the internal transcribed spacer (ITS) and β-tubulin regions. Pathogenicity tests were conducted and Koch's postulates were confirmed. The fungus was also pathogenic on Acacia mangium, but was less pathogenic on mango. Partial flooding was unfavourable for disease development. Two described isolates (WRC and WBC) had minor variation in morphology and DNA sequences, but the former exhibited a more pathogenic on both duku and acacia. The ITS phylogenies grouped the most pathogenic isolate (WRC) causing wilting of the duku tree within the aggressive and widely distributed ITS5 haplotype of C. fimbriata.

What Is the Best Lens? Comparing the Resolution Power of Genome-Derived Markers and Standard Barcodes

Fungal species delimitation was traditionally carried out with multicopy ribosomal RNA (rRNA) genes, principally for their ease of amplification. Since the efficacy of these markers has been questioned, single-copy protein-encoding genes have been proposed alone or in combination for Multi-Locus Sequence Typing (MLST). In this context, the role of the many sequences obtained with Next-Generation Sequencing (NGS) techniques, in both genomics and metagenomics, further pushes toward an analysis of the efficacy of NGS-derived markers and of the metrics to evaluate the marker efficacy in discriminating fungal species. This paper aims at proposing MeTRe (Mean Taxonomic Resolution), a novel index that could be used both for measuring marker efficacy and for assessing the actual resolution (i.e., the level of separation) between species obtained with different markers or their combinations. In this paper, we described and then employed this index to compare the efficacy of two rRNAs and four single-copy markers obtained from public databases as both an amplicon-based approach and genome-derived sequences. Two different groups of species were used, one with a pathogenic species of Candida that was characterized by relatively well-separated taxa, whereas the other, comprising some relevant species of the sensu stricto group of the genus Saccharomyces, included close species and interspecific hybrids. The results showed the ability of MeTRe to evaluate marker efficacy in general and genome-derived markers specifically.

Molecular evidence of the hybrid origin of Cryptocoryne ×purpurea Ridl. nothovar. purpurea (Araceae)

Natural hybridization has been considered a source of taxonomic complexity in Cryptocoryne. A combined study of DNA sequencing data from the internal transcribed spacer (ITS) of nuclear ribosomal DNA and the trnK-matK region of chloroplast DNA was used to identify the parents of Cryptocoryne putative hybrids from Peninsular Malaysia. Based on the intermediate morphology and sympatric distribution area, the plants were tentatively identified as the hybrid Cryptocoryne ×purpurea nothovar. purpurea. The plants were pollen sterile and had long been considered as hybrids, supposedly between two related and co-existing species, C. cordata var. cordata and C. griffithii. The status of C. ×purpurea nothovar. purpurea was independently confirmed by the presence of an additive ITS sequence pattern from these two parental species in hybrid individuals. An analysis of the chloroplast trnK-matK sequences showed that the hybridization is bidirectional with the putative hybrids sharing identical sequences from C. cordata var. cordata and C. griffithii, indicating that both putative parental species had been the maternal parent in different accessions.

Empirical evidence for concerted evolution in the 18S rDNA region of the planktonic diatom genus Chaetoceros

Concerted evolution is a process of homogenisation of repetitive sequences within a genome through unequal crossing over and gene conversion. This homogenisation is never fully achieved because mutations always create new variants. Classically, concerted evolution has been detected as "noise" in electropherograms and these variants have been characterised through cloning and sequencing of subsamples of amplified products. However, this approach limits the number of detectable variants and provides no information about the abundance of each variant. In this study, we investigated concerted evolution by using environmental time-series metabarcoding data, single strain high-throughput sequencing (HTS) and a collection of Sanger reference barcode sequences. We used six species of the marine planktonic diatom genus Chaetoceros as study system. Abundance plots obtained from environmental metabarcoding and single strain HTS showed the presence of a haplotype far more abundant than all the others (the "dominant" haplotype) and identical to the reference sequences of that species obtained with Sanger sequencing. This distribution fitted best with Zipf's law among the rank abundance/ dominance models tested. Furthermore, in each strain 99% of reads showed a similarity of 99% with the dominant haplotype, confirming the efficiency of the homogenisation mechanism of concerted evolution. We also demonstrated that minor haplotypes found in the environmental samples are not only technical artefacts, but mostly intragenomic variation generated by incomplete homogenisation. Finally, we showed that concerted evolution can be visualised inferring phylogenetic networks from environmental data. In conclusion, our study provides an important contribution to the understanding of concerted evolution and to the interpretation of DNA barcoding and metabarcoding data based on multigene family markers.

Pseudodiaptomus marinus Sato, 1913 in the Black Sea: morphology, genetic analysis, and variability in seasonal and interannual abundance

Calanoid copepod Peudodiaptomus marinus Sato, 1913 was first recorded in Sevastopol Bay in the northern Black Sea in September 2016. We performed regular observations of this new invasive species between October 2016 and December 2018. We conducted bi-weekly plankton sampling at three stations located within or adjacent to Sevastopol Bay. This is the first paper to combine (i) a detailed morphological study, (ii) molecular genetic analysis, and (iii) an investigation of P. marinus seasonal dynamics and interannual abundance variability in the coastal Black Sea. Our morphological research showed similarities between Pseudodiaptomus specimens and existing P. marinus illustrations and descriptions. Our morphological analysis results were confirmed using molecular genetic studies. Based on the genetic variability of ITS2 and cytb, we found that all investigated specimens from Sevastopol Bay belonged to P. marinus. Investigations of P. marinus seasonal and interannual abundance variability showed the same seasonal patterns throughout the studied period, with a higher seasonal abundance from October to November and one pronounced density peak in autumn. The highest abundances (2,000 ind m^–3 at the mouth of the bay and more than 5,000 ind m^–3at its centre) were recorded in November 2018. In the samples, we found adults of both sexes, including ovigerous females, copepodites I–V, and nauplii, suggesting that the species reproduce in Sevastopol Bay. Our research indicated that P. marinus is a new non-indigenous species (NIS) in the Black Sea, and we will discuss a possible vector of its introduction into this basin.

The fate of 35S rRNA genes in the allotetraploid grass Brachypodium hybridum

Nucleolar dominance (ND) consists of the reversible silencing of 35S/45S rDNA loci inherited from one of the ancestors of an allopolyploid. The molecular mechanisms by which one ancestral rDNA set is selected for silencing remain unclear. We applied a combination of molecular (Southern blot hybridization and reverse-transcription cleaved amplified polymorphic sequence analysis), genomic (analysis of variants) and cytogenetic (fluorescence in situ hybridization) approaches to study the structure, expression and epigenetic landscape of 35S rDNA in an allotetraploid grass that exhibits ND, Brachypodium hybridum (genome composition DDSS), and its putative progenitors, Brachypodium distachyon (DD) and Brachypodium stacei (SS). In progenitor genomes, B. stacei showed a higher intragenomic heterogeneity of rDNA compared with B. distachyon. In all studied accessions of B. hybridum, there was a reduction in the copy number of S homoeologues, which was accompanied by their inactive transcriptional status. The involvement of DNA methylation in CG and CHG contexts in the silencing of the S-genome rDNA loci was revealed. In the B. hybridum allotetraploid, ND is stabilized towards the D-genome units, irrespective of the polyphyletic origin of the species, and does not seem to be influenced by homoeologous 35S rDNA ratios and developmental stage.

Metschnikowia pulcherrima and Related Pulcherrimin-Producing Yeasts: Fuzzy Species Boundaries and Complex Antimicrobial Antagonism

Yeasts affiliated with the Metschnikowia pulcherrima clade (subclade) of the large ascomycetous genus Metschnikowia frequently turn out to produce the characteristic maroon-red pulcherrimin when tested for pigment production and prove to exert antagonistic effects on many types of microorganisms. The determination of the exact taxonomic position of the strains is hampered by the shortage of distinctive morphological and physiological properties of the species of the clade and the lack of rDNA barcode gaps. The rDNA repeats of the type strains of the species are not homogenized and are assumed to evolve by a birth-and-death mechanism combined with reticulation. The taxonomic division is further hampered by the incomplete biological (reproductive) isolation of the species: certain type strains can be hybridized and genome sequencing revealed chimeric genome structures in certain strains that might have evolved from interspecies hybrids (alloploid genome duplication). Various mechanisms have been proposed for the antimicrobial antagonism. One is related to pulcherrimin production. The diffusible precursor of pulcherrimin, the pulcherriminic acid is secreted by the cells into the environment where it forms the insoluble pulcherrimin with the ferric ions. The lack of free iron caused by the immobilization of ferric ions inhibits the growth of many microorganisms. Recent results of research into the complexity of the taxonomic division of the pulcherrimin-producing Metschnikowia yeasts and the mechanism(s) underlying their antimicrobial antagonism are discussed in this review.

Species, Sequence Types and Alleles: Dissecting Genetic Variation in Acanthamoeba

Species designations within Acanthamoeba are problematic because of pleomorphic morphology. Molecular approaches, including DNA sequencing, hinted at a resolution that has yet to be fully achieved. Alternative approaches were required. In 1996, the Byers/Fuerst lab introduced the concept of sequence types. Differences between isolates of Acanthamoeba could be quantitatively assessed by comparing sequences of the nuclear 18S rRNA gene, ultimately producing 22 sequence types, designated T1 through T22. The concept of sequence types helps our understanding of Acanthamoeba evolution. Nevertheless, substantial variation in the 18S rRNA gene differentiates many isolates within each sequence type. Because the majority of isolates with sequences in the international DNA databases have been studied for only a small segment of the gene, designated ASA.S1, genetic variation within this hypervariable region of the 18S rRNA gene has been scrutinized. In 2002, we first categorized variation in this region in a sample of T3 and T4 isolates from Hong Kong, observing ten “alleles” within type T4 and five “alleles” within T3. Subsequently, confusion occurred when different labs applied redundant numerical labels to identify different alleles. A more unified approach was required. We have tabulated alleles occurring in the sequences submitted to the international DNA databases, and determined their frequencies. Over 150 alleles have occurred more than once within 3500+ isolates of sequence type T4. Results from smaller samples of other sequence types (T3, T5, T11 and T15, and supergroup T2/6) have also been obtained. Our results provide new insights into the evolutionary history of Acanthamoeba, further illuminating the degree of genetic separation between significant taxonomic units within the genus, perhaps eventually elucidating what constitutes a species of Acanthamoeba.

Phylogenomic incongruence in Ceratocystis: a clue to speciation?

Background

The taxonomic history of Ceratocystis, a genus in the Ceratocystidaceae, has been beset with questions and debate. This is due to many of the commonly used species recognition concepts (e.g., morphological and biological species concepts) providing different bases for interpretation of taxonomic boundaries. Species delineation in Ceratocystis primarily relied on genealogical concordance phylogenetic species recognition (GCPSR) using multiple standard molecular markers.

Results

Questions have arisen regarding the utility of these markers e.g., ITS, BT and TEF1-α due to evidence of intragenomic variation in the ITS, as well as genealogical incongruence, especially for isolates residing in a group referred to as the Latin-American clade (LAC) of the species. This study applied a phylogenomics approach to investigate the extent of phylogenetic incongruence in Ceratocystis. Phylogenomic analyses of a total of 1121 shared BUSCO genes revealed widespread incongruence within Ceratocystis, particularly within the LAC, which was typified by three equally represented topologies. Comparative analyses of the individual gene trees revealed evolutionary patterns indicative of hybridization. The maximum likelihood phylogenetic tree generated from the concatenated dataset comprised of 1069 shared BUSCO genes provided improved phylogenetic resolution suggesting the need for multiple gene markers in the phylogeny of Ceratocystis.

Conclusion

The incongruence observed among single gene phylogenies in this study call into question the utility of single or a few molecular markers for species delineation. Although this study provides evidence of interspecific hybridization, the role of hybridization as the source of discordance will require further research because the results could also be explained by high levels of shared ancestral polymorphism in this recently diverged lineage. This study also highlights the utility of BUSCO genes as a set of multiple orthologous genes for phylogenomic studies.

Lost in translation: Population genomics and long-read sequencing reveals relaxation of concerted evolution of the ribosomal DNA cistron

Concerted evolution of the ribosomal DNA array has been studied in numerous eukaryotic taxa, yet is still poorly understood. rDNA genes are repeated dozens to hundreds of times in the eukaryotic genome (Eickbush and Eickbush, 2007) and it is believed that these arrays are homogenized through concerted evolution (Zimmer et al., 1980; Dover, 1993) preventing the accumulation of intragenomic, and intraspecific, variation. However, numerous studies have reported rampant intragenomic and intraspecific variation in the rDNA array (Ganley and Kobayashi, 2011; Naidoo et al., 2013; Hughes and Petersen, 2001; Lindner and Banik, 2011; Li et al., 2013; Lindner et al., 2013; Hughes et al., 2018), contradicting our current understanding of concerted evolution. The internal transcribed spacers (ITS) of the rDNA cistron are the most commonly used DNA barcoding region in Fungi (Schoch et al., 2012), and rely on concerted evolution to homogenize the rDNA array leading to a "barcode gap" (Puillandre et al., 2012). Here we show that in Boletus edulis Bull., ITS intragenomic variation persists at low allele frequencies throughout the rDNA array, this variation does not correlate with genomic relatedness between populations, and rDNA genes may not evolve in a strictly concerted fashion despite the presence of unequal recombination and gene conversion. Under normal assumptions, heterozygous positions found in ITS sequences represent hybridization between populations, yet through allelic mapping of the rDNA array we found numerous heterozygous alleles to be stochastically introgressed throughout, presenting a dishonest signal of gene flow. Moreover, despite the signal of gene flow in ITS, our organisms were highly inbred, indicating a disconnect between true gene flow and barcoding signals. In addition, we show that while the mechanisms of concerted evolution are ongoing in pseudo-heterozygous individuals, they are not fully homogenizing the ITS array. Concerted evolution of the rDNA array may insufficiently homogenize the ITS gene, allowing for misleading signals of gene flow to persist, vastly complicating the use of the ITS locus for DNA barcoding in Fungi.

Characterizing the ribosomal tandem repeat and its utility as a DNA barcode in lichen-forming fungi

Background

Regions within the nuclear ribosomal operon are a major tool for inferring evolutionary relationships and investigating diversity in fungi. In spite of the prevalent use of ribosomal markers in fungal research, central features of nuclear ribosomal DNA (nrDNA) evolution are poorly characterized for fungi in general, including lichenized fungi. The internal transcribed spacer (ITS) region of the nrDNA has been adopted as the primary DNA barcode identification marker for fungi. However, little is known about intragenomic variation in the nrDNA in symbiotic fungi. In order to better understand evolution of nrDNA and the utility of the ITS region for barcode identification of lichen-forming fungal species, we generated nearly complete nuclear ribosomal operon sequences from nine species in the Rhizoplaca melanophthalma species complex using short reads from high-throughput sequencing.

Results

We estimated copy numbers for the nrDNA operon, ranging from nine to 48 copies for members of this complex, and found low levels of intragenomic variation in the standard barcode region (ITS). Monophyly of currently described species in this complex was supported in phylogenetic inferences based on the ITS, 28S, intergenic spacer region, and some intronic regions, independently; however, a phylogenetic inference based on the 18S provided much lower resolution. Phylogenetic analysis of concatenated ITS and intergenic spacer sequence data generated from 496 specimens collected worldwide revealed previously unrecognized lineages in the nrDNA phylogeny.

Conclusions

The results from our study support the general assumption that the ITS region of the nrDNA is an effective barcoding marker for fungi. For the R. melanophthalma group, the limited amount of potential intragenomic variability in the ITS region did not correspond to fixed diagnostic nucleotide position characters separating taxa within this species complex. Previously unrecognized lineages inferred from ITS sequence data may represent undescribed species-level lineages or reflect uncharacterized aspects of nrDNA evolution in the R. melanophthalma species complex.

IGS sequences in Cestrum present AT- and GC-rich conserved domains, with strong regulatory potential for 5S rDNA

The 35S and 5S ribosomal DNA (rDNA) organized in thousands of copies in genomes, have been widely used in numerous comparative cytogenetic studies. Nevertheless, several questions related to the diversity and organization of regulatory motifs in 5S rDNA remain to be addressed. The 5S rDNA unit is composed of a conserved 120 bp length coding region and an intergenic spacer (IGS) containing potential regulatory motifs (Poly-T, AT-rich and GC-rich) differing in number, redundancy and position along the IGS. The Cestrum species (Solanaceae) have large genomes (about 10 pg/1C) and conserved 2n = 16 karyotypes. Strikingly, these genomes show high diversity of heterochromatin distribution, variability in 35S rDNA loci and the occurrence of B chromosomes. However, the 5S rDNA loci are highly conserved in the proximal region of chromosome 8. Comparison of seventy-one IGS sequences in plants revealed several conserved motifs with potential regulatory function. The AT- and GC-rich domains appeared highly conserved in Cestrum chromosomes. The 5S genic and the GC-rich IGS probe produced FISH signals in both A (pair 8) and B chromosomes. The GC-rich domain presented a strong potential for regulation because it may be associated with CpG islands organization, as well as to hairpin and loop organization. Another interesting aspect was the ability of AT- and GC-rich motifs to produce non-heterochromatic CMA/DAPI signals. While the length of the 5S rDNA IGS region varied in size between the Cestrum species, the individual sequence motifs seem to be conserved suggesting their regulatory function. The most striking feature was the conserved GC-rich domain in Cestrum, which is recognized as a signature trait of the proximal region of chromosome pair 8.

The extremely halotolerant black yeast Hortaea werneckii - a model for intraspecific hybridization in clonal fungi

The polymorphic black yeast Hortaea werneckii (Capnodiales, Ascomycota) is extremely halotolerant (growth from 0 to 30% [w/v] NaCl) and has been extensively studied as a model for halotolerance in Eukaryotes for over two decades. Its most frequent sources are hypersaline environments and adjacent sea-water habitats in temperate, subtropical and tropical climates. Although typically saprobic, H. werneckii can also act as a commensal coloniser on human skin, causing tinea nigra on hands and soles. Here, we report that addition of NaCl to culture media expands the growth range of H. werneckii to 37 °C, which explains its colonisation of human skin, with its increased salinity. The morphological and physiological plasticity/ versatility of H. werneckii indicate that a species complex might be involved. This was investigated in this polyphasic taxonomic analysis based on the global diversity of H. werneckii strains collected from hypersaline environments, and from humans and animals. Analysis of D1/D2domains of 28S and internal transcribed spacer rDNA revealed 10 and 17 genotypes, respectively, that were not always compliant. The genotypes have global distributions. Human and environmental strains with the same genotypes are intermingled. Due to the limited number of phylogenetically informative characters in the ribosomal DNA dataset, the partial genes encoding for β-tubulin (BTB) and mini-chromosome maintenance protein (MCM7) were also sequenced. The use of these genes was hampered by ambiguous sequences obtained by Sanger sequencing, as a consequence of the diploid and highly heterozygous genome of many H. werneckii strains. Analysis of the BTB and MCM7 genes showed that in some cases two copies of the gene from the same genome are positioned in distant phylogenetic clusters of the intraspecific gene tree. Analysis of whole-genome sequences of selected H. werneckii strains generally confirmed the phylogenetic distances estimated on the basis of ribosomal genes, but also showed substantial reticulation within the phylogenetic history of the strains. This is in line with the hypothesis that the diploid genomes of H. werneckii were formed by hybridizations, which have sometimes occurred between relatively divergent strains.

Early Ongoing Speciation of Ogataea uvarum Sp. Nov. Within the Grape Ecosystem Revealed by the Internal Variability Among the rDNA Operon Repeats

A yeast strain was isolated during a study on vineyard-associated yeast strains from Apulia in Southern Italy. ITS and LSU D1/D2 rDNA sequences showed this strain not to belong to any known species and was described as the type strain of Ogataea uvarum sp. nov., a close relative of O. philodendri. Several secondary peaks appeared in the sequences, suggesting internal heterogeneity among the copies of the rDNA. This hypothesis was tested by sequencing single clones of the marker region. The analyses showed different levels of variability throughout the operon with differences between the rRNA encoding genes and the internally transcribed regions. O. uvarum and O. philodendri share high frequency variants, i.e., variants frequently found in many clones, whereas there is a large variability of the low frequency polymorphisms, suggesting that the mechanism of homogenization is more active with the former than with the latter type of variation. These findings indicate that low frequency variants are detected in Sanger sequencing as secondary peaks whereas in Next Generation Sequencing (NGS) of metagenomics DNA would lead to an overestimate of the alpha diversity. For the first time in our knowledge, this investigation shed light on the variation of the copy number of the rDNA cistron during the yeast speciation process. These polymorphisms can be used to investigate on the processes occurring in these taxonomic markers during the separation of fungal species, it being a genetic process highly frequent in the complex microbial ecosystem existing in grape, must and wine.

NGS barcode sequencing in taxonomy and diagnostics, an application in "Candida" pathogenic yeasts with a metagenomic perspective

Species identification of yeasts and other Fungi is currently carried out with Sanger sequences of selected molecular markers, mainly from the ribosomal DNA operon, characterized by hundreds of tandem repeats of the 18S, ITS1, 5.8S, ITS2 and LSU loci. The ITS region has been recently proposed as a primary barcode marker making this region the most used one in taxonomy, phylogeny and diagnostics. The introduction of NGS is providing tools of high efficacy and relatively low cost to amplify two or more markers simultaneously with great sequencing depth. However, the presence of intra-genomic variability between the repeats requires specific analytical procedures and pipelines. In this study, 286 strains belonging to 11 pathogenic yeasts species were analysed with NGS of the region spanning from ITS1 to the D1/D2 domain of the LSU encoding ribosomal DNA. Results showed that relatively high heterogeneity can hamper the use of these sequences for the identification of single strains and even more of complex microbial mixtures. These observations point out that the metagenomics studies could be affected by species inflection at levels higher than currently expected.

New Ceratocystis species associated with rapid death of Metrosideros polymorpha in Hawai'i

The native ‘ōhi’a lehua (Metrosideros polymorpha) has cultural, biological and ecological significance to Hawai’i, but it is seriously threatened by a disease commonly referred to as rapid ‘ōhi’a death (ROD). Preliminary investigations showed that a Ceratocystis species similar to C. fimbriata s.lat. was the cause of the disease. In this study, we used a combination of the phylogenetic, morphological and biological species concepts, as well as pathogenicity tests and microsatellite analyses, to characterise isolates collected from diseased ‘ōhi’a trees across Hawai’i Island. Two distinct lineages, representing new species of Ceratocystis, were evident based on multigene phylogenetic analyses. These are described here as C. lukuohia and C. huliohia. Ceratocystis lukuohia forms part of the Latin American clade (LAC) and was most closely associated with isolates from Syngonium and Xanthosoma from the Caribbean and elsewhere, including Hawai’i, and C. platani, which is native to eastern USA. Ceratocystis huliohia resides in the Asian-Australian clade (AAC) and is most closely related to C. uchidae, C. changhui and C. cercfabiensis, which are thought to be native to Asia. Morphology and interfertility tests support the delineation of these two new species and pathogenicity tests show that both species are aggressive pathogens on seedlings of M. polymorpha. Characterisation of isolates using microsatellite markers suggest that both species are clonal and likely represent recently-introduced strains. Intensive research is underway to develop rapid screening protocols for early detection of the pathogens and management strategies in an attempt to prevent the spread of the pathogens to the other islands of Hawai’i, which are currently disease free.

Marked intra-genomic variation and pseudogenes in the ITS1-5.8S-ITS2 rDNA of Symphurus plagiusa (Pleuronectiformes: Cynoglossidae)

The eukaryotic ribosomal DNA (rDNA) cluster consists of multiple copies of three genes (18S, 5.8S, and 28S rDNA) and two internal transcribed spacers (ITS1 and ITS2). In recent years, an increasing number of rDNA sequence polymorphisms have been identified in numerous species. In the present study, we provide 33 complete ITS (ITS1-5.8S-ITS2) sequences from two Symphurus plagiusa individuals. To the best of our knowledge, these sequences are the first detailed information on ITS sequences in Pleuronectiformes. Here, two divergent types (Type A and B) of the ITS1-5.8S-ITS2 rDNA sequence were found, which mainly differ in sequence length, GC content, nucleotide diversity (π), secondary structure and minimum free energy. The ITS1-5.8S-ITS2 rDNA sequence of Type B was speculated to be a putative pseudogene according to pseudogene identification criteria. Cluster analysis showed that sequences from the same type clustered into one group and two major groups were formed. The high degree of ITS1-5.8S-ITS2 sequence polymorphism at the intra-specific level indicated that the S. plagiusa genome has evolved in a non-concerted evolutionary manner. These results not only provide useful data for ribosomal pseudogene identification, but also further contribute to the study of rDNA evolution in teleostean genomes.

SMORE: Synteny Modulator of Repetitive Elements

Several families of multicopy genes, such as transfer ribonucleic acids (tRNAs) and ribosomal RNAs (rRNAs), are subject to concerted evolution, an effect that keeps sequences of paralogous genes effectively identical. Under these circumstances, it is impossible to distinguish orthologs from paralogs on the basis of sequence similarity alone. Synteny, the preservation of relative genomic locations, however, also remains informative for the disambiguation of evolutionary relationships in this situation. In this contribution, we describe an automatic pipeline for the evolutionary analysis of such cases that use genome-wide alignments as a starting point to assign orthology relationships determined by synteny. The evolution of tRNAs in primates as well as the history of the Y RNA family in vertebrates and nematodes are used to showcase the method. The pipeline is freely available.

Evolution of ribosomal internal transcribed spacers in Deuterostomia

Sequences of ribosomal internal transcribed spacers (ITSs) are of great importance to molecular phylogenetics and DNA barcoding, but remain unstudied in some large taxa of Deuterostomia. We have analyzed complete ITS1 and ITS2 sequences in 62 species from 16 Deuterostomia classes, with ITS sequences in 24 species from 11 classes initially obtained using unannotated contigs and raw read sequences. A general tendency for both ITS length and GC-content increase from interior to superior Deuterostomia taxa, a uniform GC-content in both ITSs within the same species, thymine content decrease in sense DNA sequences of both ITSs are shown. A possible role of GC-based gene conversion in Deuterostomia ITS evolutionary changes is hypothesized. The first example of non-LTR retrotransposon insertion into ITS sequence in Deuterostomia is described in turtle Geochelone nigra. The roles of mobile genetic element insertions in the evolution of ITS sequences in some Sauropsida taxa are discussed as well.

Evolutionary and polymorphism analyses reveal the central role of BTN3A2 in the concerted evolution of the BTN3 gene family

The butyrophilin 3 (BTN3) receptors are implicated in the T lymphocytes regulation and present a wide plasticity in mammals. In order to understand how these genes have been diversified, we studied their evolution and show that the three human BTN3 are the result of two successive duplications in Primates and that the three genes are present in Hominoids and the Old World Monkey groups. A thorough phylogenetic analysis reveals a concerted evolution of BTN3 characterized by a strong and recurrent homogenization of the region encoding the signal peptide and the immunoglobulin variable (IgV) domain in Hominoids, where the sequences of BTN3A1 or BTN3A3 are replaced by BTN3A2 sequence. In human, the analysis of the diversity of these genes in 1683 individuals representing 26 worldwide populations shows that the three genes are polymorphic, with more than 46 alleles for each gene, and marked by extreme homogenization of the IgV sequences. The same analysis performed for the BTN2 genes shows also a concerted evolution; however, it is not as strong and recurrent as for BTN3. This study shows that BTN3 receptors are marked by extreme concerted evolution at the IgV domain and that BTN3A2 plays a central role in this evolution.

Genera of phytopathogenic fungi: GOPHY 1

Genera of Phytopathogenic Fungi (GOPHY) is introduced as a new series of publications in order to provide a stable platform for the taxonomy of phytopathogenic fungi. This first paper focuses on 21 genera of phytopathogenic fungi: Bipolaris, Boeremia, Calonectria, Ceratocystis, Cladosporium, Colletotrichum, Coniella, Curvularia, Monilinia, Neofabraea, Neofusicoccum, Pilidium, Pleiochaeta, Plenodomus, Protostegia, Pseudopyricularia, Puccinia, Saccharata, Thyrostroma, Venturia and Wilsonomyces. For each genus, a morphological description and information about its pathology, distribution, hosts and disease symptoms are provided. In addition, this information is linked to primary and secondary DNA barcodes of the presently accepted species, and relevant literature. Moreover, several novelties are introduced, i.e. new genera, species and combinations, and neo-, lecto- and epitypes designated to provide a stable taxonomy. This first paper includes one new genus, 26 new species, ten new combinations, and four typifications of older names.

Orthologs, turn-over, and remolding of tRNAs in primates and fruit flies

Background

Transfer RNAs (tRNAs) are ubiquitous in all living organism. They implement the genetic code so that most genomes contain distinct tRNAs for almost all 61 codons. They behave similar to mobile elements and proliferate in genomes spawning both local and non-local copies. Most tRNA families are therefore typically present as multicopy genes. The members of the individual tRNA families evolve under concerted or rapid birth-death evolution, so that paralogous copies maintain almost identical sequences over long evolutionary time-scales. To a good approximation these are functionally equivalent. Individual tRNA copies thus are evolutionary unstable and easily turn into pseudogenes and disappear. This leads to a rapid turnover of tRNAs and often large differences in the tRNA complements of closely related species. Since tRNA paralogs are not distinguished by sequence, common methods cannot not be used to establish orthology between tRNA genes.

Results

In this contribution we introduce a general framework to distinguish orthologs and paralogs in gene families that are subject to concerted evolution. It is based on the use of uniquely aligned adjacent sequence elements as anchors to establish syntenic conservation of sequence intervals. In practice, anchors and intervals can be extracted from genome-wide multiple sequence alignments. Syntenic clusters of concertedly evolving genes of different families can then be subdivided by list alignments, leading to usually small clusters of candidate co-orthologs. On the basis of recent advances in phylogenetic combinatorics, these candidate clusters can be further processed by cograph editing to recover their duplication histories. We developed a workflow that can be conceptualized as stepwise refinement of a graph of homologous genes. We apply this analysis strategy with different types of synteny anchors to investigate the evolution of tRNAs in primates and fruit flies. We identified a large number of tRNA remolding events concentrated at the tips of the phylogeny. With one notable exception all phylogenetically old tRNA remoldings do not change the isoacceptor class.

Conclusions

Gene families evolving under concerted evolution are not amenable to classical phylogenetic analyses since paralogs maintain identical, species-specific sequences, precluding the estimation of correct gene trees from sequence differences. This leaves conservation of syntenic arrangements with respect to “anchor elements” that are not subject to concerted evolution as the only viable source of phylogenetic information. We have demonstrated here that a purely synteny-based analysis of tRNA gene histories is indeed feasible. Although the choice of synteny anchors influences the resolution in particular when tight gene clusters are present, and the quality of sequence alignments, genome assemblies, and genome rearrangements limits the scope of the analysis, largely coherent results can be obtained for tRNAs. In particular, we conclude that a large fraction of the tRNAs are recent copies. This proliferation is compensated by rapid pseudogenization as exemplified by many very recent alloacceptor remoldings.

Electronic supplementary material

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

Contrasting Patterns of rDNA Homogenization within the Zygosaccharomyces rouxii Species Complex

Arrays of repetitive ribosomal DNA (rDNA) sequences are generally expected to evolve as a coherent family, where repeats within such a family are more similar to each other than to orthologs in related species. The continuous homogenization of repeats within individual genomes is a recombination process termed concerted evolution. Here, we investigated the extent and the direction of concerted evolution in 43 yeast strains of the Zygosaccharomyces rouxii species complex (Z. rouxii, Z. sapae, Z. mellis), by analyzing two portions of the 35S rDNA cistron, namely the D1/D2 domains at the 5’ end of the 26S rRNA gene and the segment including the internal transcribed spacers (ITS) 1 and 2 (ITS regions). We demonstrate that intra-genomic rDNA sequence variation is unusually frequent in this clade and that rDNA arrays in single genomes consist of an intermixing of Z. rouxii, Z. sapae and Z. mellis-like sequences, putatively evolved by reticulate evolutionary events that involved repeated hybridization between lineages. The levels and distribution of sequence polymorphisms vary across rDNA repeats in different individuals, reflecting four patterns of rDNA evolution: I) rDNA repeats that are homogeneous within a genome but are chimeras derived from two parental lineages via recombination: Z. rouxii in the ITS region and Z. sapae in the D1/D2 region; II) intra-genomic rDNA repeats that retain polymorphisms only in ITS regions; III) rDNA repeats that vary only in their D1/D2 domains; IV) heterogeneous rDNA arrays that have both polymorphic ITS and D1/D2 regions. We argue that an ongoing process of homogenization following allodiplodization or incomplete lineage sorting gave rise to divergent evolutionary trajectories in different strains, depending upon temporal, structural and functional constraints. We discuss the consequences of these findings for Zygosaccharomyces species delineation and, more in general, for yeast barcoding.

A possible centre of diversity in South East Asia for the tree pathogen, Ceratocystis manginecans

The fungal pathogen, Ceratocystis manginecans, has caused serious canker and wilt disease on Mangifera indica (mango), on legume tree species in Oman and Pakistan and on Acacia spp. in Indonesia. A Ceratocystis species, with similar morphology to C. manginecans, has recently been reported in Vietnam, causing severe disease of Acacia trees. Previous population genetic studies on isolates from M. indica in Oman and Pakistan have shown that the pathogen represents a single clonal haplotype, indicative of an introduced pathogen. The aim of this study was to investigate the genetic diversity and population structure of 160 C. manginecans isolates, from four host-associated populations from Oman, Pakistan, Indonesia and Vietnam. This was done by applying a combination of 14 previously developed microsatellite markers and a new set, designed in this study from two different C. manginecans genomes. Sequence data confirmed that the isolates in Vietnam are the same species as those in Indonesia and were thus identified as C. manginecans. Unlike the populations in Oman and Pakistan, relatively high levels of genetic variation were found for the isolates from Indonesia and Vietnam. The Vietnam population was significantly differentiated from the other populations and isolates from this area had the highest level of genetic diversity thus far encountered for the pathogen.

Intra-Genomic Internal Transcribed Spacer Region Sequence Heterogeneity and Molecular Diagnosis in Clinical Microbiology

Internal transcribed spacer region (ITS) sequencing is the most extensively used technology for accurate molecular identification of fungal pathogens in clinical microbiology laboratories. Intra-genomic ITS sequence heterogeneity, which makes fungal identification based on direct sequencing of PCR products difficult, has rarely been reported in pathogenic fungi. During the process of performing ITS sequencing on 71 yeast strains isolated from various clinical specimens, direct sequencing of the PCR products showed ambiguous sequences in six of them. After cloning the PCR products into plasmids for sequencing, interpretable sequencing electropherograms could be obtained. For each of the six isolates, 10-49 clones were selected for sequencing and two to seven intra-genomic ITS copies were detected. The identities of these six isolates were confirmed to be Candida glabrata (n=2), Pichia (Candida) norvegensis (n=2), Candida tropicalis (n=1) and Saccharomyces cerevisiae (n=1). Multiple sequence alignment revealed that one to four intra-genomic ITS polymorphic sites were present in the six isolates, and all these polymorphic sites were located in the ITS1 and/or ITS2 regions. We report and describe the first evidence of intra-genomic ITS sequence heterogeneity in four different pathogenic yeasts, which occurred exclusively in the ITS1 and ITS2 spacer regions for the six isolates in this study.

Species or Genotypes? Reassessment of Four Recently Described Species of the Ceratocystis Wilt Pathogen, Ceratocystis fimbriata, on Mangifera indica

Ceratocystis wilt is among the most important diseases on mango (Mangifera indica) in Brazil, Oman, and Pakistan. The causal agent was originally identified in Brazil as Ceratocystis fimbriata, which is considered by some as a complex of many cryptic species, and four new species on mango trees were distinguished from C. fimbriata based on variation in internal transcribed spacer sequences. In the present study, phylogenetic analyses using DNA sequences of mating type genes, TEF-1α, and β-tubulin failed to identify lineages corresponding to the four new species names. Further, mating experiments found that the mango isolates representing the new species were interfertile with each other and a tester strain from sweet potato (Ipomoea batatas), on which the name C. fimbriata is based, and there was little morphological variation among the mango isolates. Microsatellite markers found substantial differentiation among mango isolates at the regional and population levels, but certain microsatellite genotypes were commonly found in multiple populations, suggesting that these genotypes had been disseminated in infected nursery stock. The most common microsatellite genotypes corresponded to the four recently named species (C. manginecans, C. acaciivora, C. mangicola, and C. mangivora), which are considered synonyms of C. fimbriata. This study points to the potential problems of naming new species based on introduced genotypes of a pathogen, the value of an understanding of natural variation within and among populations, and the importance of phenotype in delimiting species.

New Ceratocystis species from Eucalyptus and Cunninghamia in South China

During routine surveys for possible fungal pathogens in the rapidly expanding plantations of Eucalyptus and Cunninghamia lanceolata in China, numerous isolates of unknown species in the genus Ceratocystis (Microascales) were obtained from tree wounds. In this study we identified the Ceratocystis isolates from Eucalyptus and Cunninghamia in the GuangDong, GuangXi, FuJian and HaiNan Provinces of South China based on morphology and through comparisons of DNA sequence data for the ITS, partial β-tubulin and TEF-1α gene regions. Morphological and DNA sequence comparisons revealed two previously unknown species residing in the Indo-Pacific Clade. These are described here as Ceratocystis cercfabiensis sp. nov. and Ceratocystis collisensis sp. nov. Isolates of Ceratocystis cercfabiensis showed intragenomic variation in their ITS sequences and four strains were selected for cloning of the ITS gene region. Twelve ITS haplotypes were obtained from 17 clones selected for sequencing, differing in up to seven base positions and representing two separate phylogenetic groups. This is the first evidence of multiple ITS types in isolates of Ceratocystis residing in the Indo-Pacific Clade. Caution should thus be exercised when using the ITS gene region as a barcoding marker for Ceratocystis species in this clade. This study also represents the first record of a species of Ceratocystis from Cunninghamia.

Redefining Ceratocystis and allied genera

The genus Ceratocystis was established in 1890 and accommodates many important fungi. These include serious plant pathogens, significant insect symbionts and agents of timber degradation that result in substantial economic losses. Virtually since its type was described from sweet potatoes, the taxonomy of Ceratocystis has been confused and vigorously debated. In recent years, particulary during the last two decades, it has become very obvious that this genus includes a wide diversity of very different fungi. These have been roughly lumped together due to their similar morphological structures that have clearly evolved through convergent evolution linked to an insect-associated ecology. As has been true for many other groups of fungi, the emergence of DNA-based sequence data and associated phylogenetic inferences, have made it possible to robustly support very distinct boundaries defined by morphological characters and ecological differences. In this study, DNA-sequence data for three carefully selected gene regions (60S, LSU, MCM7) were generated for 79 species residing in the aggregate genus Ceratocystis sensu lato and these data were subjected to rigorous phylogenetic analyses. The results made it possible to distinguish seven major groups for which generic names have been chosen and descriptions either provided or emended. The emended genera included Ceratocystis sensu stricto, Chalaropsis, Endoconidiophora, Thielaviopsis, and Ambrosiella, while two new genera, Davidsoniella and Huntiella, were described. In total, 30 new combinations have been made. This major revision of the generic boundaries in the Ceratocystidaceae will simplify future treatments and work with an important group of fungi including distantly related species illogically aggregated under a single name.

ITS polymorphisms shed light on hybrid evolution in apomictic plants: a case study on the Ranunculus auricomus complex

The reconstruction of reticulate evolutionary histories in plants is still a major methodological challenge. Sequences of the ITS nrDNA are a popular marker to analyze hybrid relationships, but variation of this multicopy spacer region is affected by concerted evolution, high intraindividual polymorphism, and shifts in mode of reproduction. The relevance of changes in secondary structure is still under dispute. We aim to shed light on the extent of polymorphism within and between sexual species and their putative natural as well as synthetic hybrid derivatives in the Ranunculus auricomus complex to test morphology-based hypotheses of hybrid origin and parentage of taxa. We employed direct sequencing of ITS nrDNA from 68 individuals representing three sexuals, their synthetic hybrids and one sympatric natural apomict, as well as cloning of ITS copies in four representative individuals, RNA secondary structure analysis, and landmark geometric morphometric analysis on leaves. Phylogenetic network analyses indicate additivity of parental ITS variants in both synthetic and natural hybrids. The triploid synthetic hybrids are genetically much closer to their maternal progenitors, probably due to ploidy dosage effects, although exhibiting a paternal-like leaf morphology. The natural hybrids are genetically and morphologically closer to the putative paternal progenitor species. Secondary structures of ITS1-5.8S-ITS2 were rather conserved in all taxa. The observed similarities in ITS polymorphisms suggest that the natural apomict R. variabilis is an ancient hybrid of the diploid sexual species R. notabilis and the sexual species R. cassubicifolius. The additivity pattern shared by R. variabilis and the synthetic hybrids supports an evolutionary and biogeographical scenario that R. variabilis originated from ancient hybridization. Concerted evolution of ITS copies in R. variabilis is incomplete, probably due to a shift to asexual reproduction. Under the condition of comprehensive inter- and intraspecific sampling, ITS polymorphisms are powerful for elucidating reticulate evolutionary histories.