Complete Sructure of Nuclear rDNA of the Obligate Plant Parasite Plasmodiophora brassicae: Intraspecific Polymorphisms in the Exon and Group I Intron of the Large Subunit rDNA

Genome-Wide Mining of CULLIN E3 Ubiquitin Ligase Genes from Uncaria rhynchophylla

CULLIN (CUL) protein is a subtype of E3 ubiquitin ligase that is involved in a variety of biological processes and responses to stress in plants. In Uncaria rhynchophylla, the CUL gene family has not been identified and its role in plant development, stress response and secondary metabolite synthesis has not been studied. In this study, 12 UrCUL gene members all contained the typical N-terminal domain and C-terminal domain identified from the U. rhynchophylla genome and were classified into four subfamilies based on the phylogenetic relationship with CULs in Arabidopsis thaliana. They were unevenly distributed on eight chromosomes but had a similar structural composition in the same subfamily, indicating that they were relatively conserved and potentially had similar gene functions. An interspecific and intraspecific collinearity analysis showed that fragment duplication played an important role in the evolution of the CUL gene family. The analysis of the cis-acting elements suggests that the UrCULs may play an important role in various biological processes, including the abscisic acid (ABA) response. To investigate this hypothesis, we treated the roots of U. rhynchophylla tissue-cultured seedlings with ABA. The expression pattern analysis showed that all the UrCUL genes were widely expressed in roots with various expression patterns. The co-expression association analysis of the UrCULs and key enzyme genes in the terpenoid indole alkaloid (TIA) synthesis pathway revealed the complex expression patterns of 12 UrCUL genes and some key TIA enzyme genes, especially UrCUL1, UrCUL1-likeA, UrCUL2-likeA and UrCUL2-likeB, which might be involved in the biosynthesis of TIAs. The results showed that the UrCULs were involved in the response to ABA hormones, providing important information for elucidating the function of UrCULs in U. rhynchophylla. The mining of UrCULs in the whole genome of U. rhynchophylla provided new information for understanding the CUL gene and its function in plant secondary metabolites, growth and development.

Hidden Fungi: Combining Culture-Dependent and -Independent DNA Barcoding Reveals Inter-Plant Variation in Species Richness of Endophytic Root Fungi in Elymus repens

The root endophyte community of the grass species Elymus repens was investigated using both a culture-dependent approach and a direct amplicon sequencing method across five sites and from individual plants. There was much heterogeneity across the five sites and among individual plants. Focusing on one site, 349 OTUs were identified by direct amplicon sequencing but only 66 OTUs were cultured. The two approaches shared ten OTUs and the majority of cultured endophytes do not overlap with the amplicon dataset. Media influenced the cultured species richness and without the inclusion of 2% MEA and full-strength MEA, approximately half of the unique OTUs would not have been isolated using only PDA. Combining both culture-dependent and -independent methods for the most accurate determination of root fungal species richness is therefore recommended. High inter-plant variation in fungal species richness was demonstrated, which highlights the need to rethink the scale at which we describe endophyte communities.

Comparative analysis of the ribosomal DNA repeat unit (rDNA) of Perna viridis (Linnaeus, 1758) and Perna canaliculus (Gmelin, 1791)

Perna viridis and P. canaliculus are economically and ecologically important species of shellfish. In this study, the complete ribosomal DNA (rDNA) unit sequences of these species were determined for the first time. The gene order, 18S rRNA-internal transcribed spacer (ITS) 1-5.8S rRNA-ITS2-28S rRNA-intergenic spacer (IGS), was similar to that observed in other eukaryotes. The lengths of the P. viridis and P. canaliculus rDNA sequences ranged from 8,432 to 8,616 bp and from 7,597 to 7,610 bp, respectively, this variability was mainly attributable to the IGS region. The putative transcription termination site and initiation site were confirmed. Perna viridis and P. canaliculus rDNA contained two (length: 93 and 40 bp) and one (length: 131 bp) repeat motifs, respectively. Individual intra-species differences mainly involved the copy number of repeat units. In P. viridis, three cytosine-guanine (CpG) sites with sizes of 440, 1,075 and 537 bp were found to cover nearly the entire IGS sequence, whereas in P. canaliculus, two CpG islands with sizes of 361 and 484 bp were identified. The phylogenetic trees constructed with maximum likelihood and neighbour-joining methods and based on ITS sequences were identical and included three major clusters. Species of the same genus were easily clustered together.

Infection of Plasmodiophora brassicae changes the fungal endophyte community of tumourous stem mustard roots as revealed by high-throughput sequencing and culture-dependent methods

Diverse fungal endophytes live in plants and are shaped by some abiotic and biotic stresses. Plant disease as particular biotic stress possibly gives an impact on the communities of fungal endophytes. In this study, clubroot disease caused by an obligate biotroph protist, Plasmodiophora brassicae, was considered to analyze its influence on the fungal endophyte community using an internal transcribed spacer (ITS) through high-throughput sequencing and culture-dependent methods. The results showed that the diversity of the endophyte community in the healthy roots was much higher than the clubroots. Ascomycota was the dominant group of endophytes (Phoma, Mortierella, Penicillium, etc.) in the healthy roots while P. brassicae was the dominant taxon in the clubroots. Hierarchical clustering, principal component analysis (PCA), principal coordinates analysis (PCoA) and analysis of similarities (ANOSIM) indicated significant differences between the endophyte communities in the healthy roots and clubroots. Linear discriminant analysis effect size (LefSe) analysis showed that the dominant genera could be regarded as potential biomarkers. The endophyte community in the healthy roots had a more complex network compared with the clubroots. Also, many plant pathogenic Fusarium were isolated from the clubroots by the culture-dependent method. The outcome of this study illustrates that P. brassicae infection may change the fungal endophyte community associated with the roots of tumourous stem mustard and facilitates the entry of soil pathogen into the roots.

Combinations of Independent Dominant Loci Conferring Clubroot Resistance in All Four Turnip Accessions (Brassica rapa) From the European Clubroot Differential Set

Clubroot disease is devastating to Brassica crop production when susceptible cultivars are planted in infected fields. European turnips are the most resistant sources and their resistance genes have been introduced into other crops such oilseed rape (Brassica napus L.), Chinese cabbage and other Brassica vegetables. The European clubroot differential (ECD) set contains four turnip accessions (ECD1-4). These ECD turnips exhibited high levels of resistance to clubroot when they were tested under controlled environmental conditions with Canadian field isolates. Gene mapping of the clubroot resistance genes in ECD1-4 were performed and three independent dominant resistance loci were identified. Two resistance loci were mapped on chromosome A03 and the third on chromosome A08. Each ECD turnip accession contained two of these three resistance loci. Some resistance loci were homozygous in ECD accessions while others showed heterozygosity based on the segregation of clubroot resistance in 20 BC₁ families derived from ECD1 to 4. Molecular markers were developed linked to each clubroot resistance loci for the resistance gene introgression in different germplasm.

Combinations of Independent Dominant Loci Conferring Clubroot Resistance in All Four Turnip Accessions (Brassica rapa) From the European Clubroot Differential Set

Clubroot disease is devastating to Brassica crop production when susceptible cultivars are planted in infected fields. European turnips are the most resistant sources and their resistance genes have been introduced into other crops such oilseed rape (Brassica napus L.), Chinese cabbage and other Brassica vegetables. The European clubroot differential (ECD) set contains four turnip accessions (ECD1-4). These ECD turnips exhibited high levels of resistance to clubroot when they were tested under controlled environmental conditions with Canadian field isolates. Gene mapping of the clubroot resistance genes in ECD1-4 were performed and three independent dominant resistance loci were identified. Two resistance loci were mapped on chromosome A03 and the third on chromosome A08. Each ECD turnip accession contained two of these three resistance loci. Some resistance loci were homozygous in ECD accessions while others showed heterozygosity based on the segregation of clubroot resistance in 20 BC₁ families derived from ECD1 to 4. Molecular markers were developed linked to each clubroot resistance loci for the resistance gene introgression in different germplasm.

Detection of Ribosomal DNA Sequence Polymorphisms in the Protist Plasmodiophora brassicae for the Identification of Geographical Isolates

Clubroot is a soil-borne disease caused by the protist Plasmodiophora brassicae (P. brassicae). It is one of the most economically important diseases of Brassica rapa and other cruciferous crops as it can cause remarkable yield reductions. Understanding P. brassicae genetics, and developing efficient molecular markers, is essential for effective detection of harmful races of this pathogen. Samples from 11 Korean field populations of P. brassicae (geographic isolates), collected from nine different locations in South Korea, were used in this study. Genomic DNA was extracted from the clubroot-infected samples to sequence the ribosomal DNA. Primers and probes for P. brassicae were designed using a ribosomal DNA gene sequence from a Japanese strain available in GenBank (accession number AB526843; isolate NGY). The nuclear ribosomal DNA (rDNA) sequence of P. brassicae, comprising 6932 base pairs (bp), was cloned and sequenced and found to include the small subunits (SSUs) and a large subunit (LSU), internal transcribed spacers (ITS1 and ITS2), and a 5.8s. Sequence variation was observed in both the SSU and LSU. Four markers showed useful differences in high-resolution melting analysis to identify nucleotide polymorphisms including single- nucleotide polymorphisms (SNPs), oligonucleotide polymorphisms, and insertions/deletions (InDels). A combination of three markers was able to distinguish the geographical isolates into two groups.

The Large Subunit rDNA Sequence of Plasmodiophora brassicae Does not Contain Intra-species Polymorphism

Clubroot disease caused by Plasmodiophora brassicae is one of the most important diseases of cultivated brassicas. P. brassicae occurs in pathotypes which differ in the aggressiveness towards their Brassica host plants. To date no DNA based method to distinguish these pathotypes has been described. In 2011 polymorphism within the 28S rDNA of P. brassicae was reported which potentially could allow to distinguish pathotypes without the need of time-consuming bioassays. However, isolates of P. brassicae from around the world analysed in this study do not show polymorphism in their LSU rDNA sequences. The previously described polymorphism most likely derived from soil inhabiting Cercozoa more specifically Neoheteromita-like glissomonads. Here we correct the LSU rDNA sequence of P. brassicae. By using FISH we demonstrate that our newly generated sequence belongs to the causal agent of clubroot disease.

History of oilseed rape cropping and geographic origin affect the genetic structure of Plasmodiophora brassicae populations

The soilborne pathogen Plasmodiophora brassicae causes clubroot on Brassica crops, a common disease in many oilseed rape growing regions. Here, we investigate genetic diversity and geographic differentiation of P. brassicae populations from different regions in Germany. We compared three regions that differ in oilseed rape cropping history, oilseed rape acreage, and incidence of clubroot. These regions were either spatially separated or separated by the former inner German border. Plasmodiophora isolates were collected from 59 fields (29, 17, and 13 fields per region, respectively) and 174 amplified fragment length polymorphism (AFLP) markers were analyzed. Every field isolate showed a unique genotype pattern; that is, no genotype was shared among the regions and different fields. The mean gene diversity was 0.27, suggesting that P. brassicae is a genetically diverse species. The comparison of indexes (gene diversity, genotypic diversity, and linkage disequilibrium) between the regions does not support our hypotheses that cropping history, oilseed rape acreage, and incidence of clubroot affect these estimates. Principal component analysis (PCA), fixation index (FST), and generalized linear model (GLM) were suitable to specify regional differences. PCA revealed two clusters of isolates based on the geographic origin of the isolates and FST showed that these clusters were highly differentiated. Hypotheses about association of genotypes with different spatial scales were tested with GLM: the region, reflecting the cropping history, and the individual field had a significant effect on the AFLP pattern. We propose that individual field isolates represent a discrete population and that geographic differentiation results from low levels of gene flow due to the limited dispersal of this soilborne pathogen and from localized selection pressure as unifying force on the genotypes.

Identification, molecular characterization, and evolution of group I introns at the expansion segment D11 of 28S rDNA in Rhizoctonia species

The nuclear ribosomal DNA of Rhizoctonia species is polymorphic in terms of the nucleotide composition and length. Insertions of 349-410 nucleotides in length with characteristics of group I introns were detected at a single insertion point at the expansion segment D11 of 28S rDNA in 12 out of 64 isolates. Eleven corresponded to Rhizoctonia solani (teleomorph: Thanatephorous) and one (AG-Q) to Rhizoctonia spp. (teleomorph: Ceratobasidium). Sequence data showed that all but AG-Q contained conserved DNA catalytic core regions (P, Q, R, and S) essential for selfsplicing. The predicted secondary structure revealed that base-paired helices corresponded to subgroup IC1. Isolates from same anastomosis group and even subgroups within R. solani were variable with regard to possession of introns. Phylogenetic analyses indicated that introns were vertically transmitted. Unfortunately, sequence data from the conserved region from all 64 isolates were not useful for delimiting species. Analyses with IC1 introns at same insertion point, of both Ascomycota and Basidiomycota indicated the possibility of horizontal transfer at this site. The present study uncovered new questions on evolutionary pattern of change of these introns within Rhizoctonia species.

Identification and characterization of Crr1a, a gene for resistance to clubroot disease (Plasmodiophora brassicae Woronin) in Brassica rapa L

Clubroot disease, caused by the obligate biotrophic protist Plasmodiophora brassicae Woronin, is one of the most economically important diseases of Brassica crops in the world. Although many clubroot resistance (CR) loci have been identified through genetic analysis and QTL mapping, the molecular mechanisms of defense responses against P. brassicae remain unknown. Fine mapping of the Crr1 locus, which was originally identified as a single locus, revealed that it comprises two gene loci, Crr1a and Crr1b. Here we report the map-based cloning and characterization of Crr1a, which confers resistance to clubroot in Brassica rapa. Crr1a(G004), cloned from the resistant line G004, encodes a Toll-Interleukin-1 receptor/nucleotide-binding site/leucine-rich repeat (TIR-NB-LRR) protein expressed in the stele and cortex of hypocotyl and roots, where secondary infection of the pathogen occurs, but not in root hairs, where primary infection occurs. Gain-of-function analysis proved that Crr1a(G004) alone conferred resistance to isolate Ano-01 in susceptible Arabidopsis and B. rapa. In comparison, the susceptible allele Crr1a(A9709) encodes a truncated NB-LRR protein, which lacked more than half of the TIR domain on account of the insertion of a solo-long terminal repeat (LTR) in exon 1 and included several substitutions and insertion-deletions in the LRR domain. This study provides a basis for further molecular analysis of defense mechanisms against P. brassicae and will contribute to the breeding of resistant cultivars of Brassica vegetables by marker-assisted selection.Data deposition The sequence reported in this paper has been deposited in the GenBank database (accession no. AB605024).