The Amplification and Application of Ribosomal RNA (rDNA) Gene Sequences of Blidingia minima (Chlorophyta, Blidingia)

The sequence of the ribosomal RNA gene (rDNA) plays 
an important role in species identification and phylogenetic analysis. However, 
the only published full-length sequence of a 
ribosomal gene of green algae is that of Ulva mutabilis. In this 
study, we amplified the full-length sequence of each 
ribosomal gene unit of the ribosomal gene of Blidingia minima. The 
full-length sequence of the ribosomal gene in Blidingia minima is 8676 
bp, including the 1759 bp 18S rDNA, 576 bp internal transcribed spacer (ITS) + 5.8S 
rDNA, 3282 bp 28S rDNA, and 3059 bp 
intergenic spacer (IGS) region. We then carried out a series of genetic 
analyses based on the ITS and IGS sequences, to verify whether IGS sequences are useful for studying the genetic 
diversity of green algae from different locations. We amplified the IGS 
sequences of Blidingia minima from 10 different locations in the Yellow 
Sea. Multiple alignments of the IGS sequences of samples from these 10 different sites revealed varying degrees of base differences, and comparative analysis 
of the ITS sequences revealed that our amplified species was classified as Blidingia 
minima and distinct from other green algae. In conclusion, our full-length 
amplified ribosomal gene provides useful information to enrich the data on 
green algae ribosomal genes and provides an effective molecular marker 
for the analysis of the interspecies and intraspecies relationships of Blidingia minima.


Introduction
combined with increased temperatures [1] [2], and such events have a negative impact on the local economy and the environment. In addition, green tides can be dispersed by wind and currents to other coastal areas [3]. The world's largest green tide events were mainly caused by Ulva prolifera O. F. Müller (Chlorophyta, Ulvophyceae), which erupted in 2008 in Qingdao on the Yellow Sea coast in north-eastern China. Based on satellite remote sensing images, it has been postulated that the origin of the green tide was closely related to the expansion of Pyropia aquaculture regions along the coastline of Jiangsu Province [4] [5].
There was molecular evidence that the species attached on the raft in Jiangsu Province and the species that caused the green tide were the same [3]. In addition to U. prolifera, a variety of other green algae were also identified on the Pyropia farming raft. Another of the fixed green algae species growing on the Pyropia farming rafts throughout the year is a Blidingia species with extensive biomass ( Figure 1). Blidingia has similar habits to U. prolifera on the Pyropia farming rafts. The growth of large quantities of green algae on the Pyropia farming raft, especially Ulva prolifera and Blidingia minima, causes large economic losses to the aquaculture industries in China. During the outbreak of green tide, the biomass of Ulva prolifera is predominant. The biomass of other green algae, such as Ulva flexuosa and Ulva compressa, gradually decreased, but the amount of Blidingia minima remained almost unchanged. Few studies have shown why Blidingia minima can grow on Pyropia farming rafts throughout the year. Therefore, our amplified ribosomal gene cluster can provide some molecular information for subsequent research.
In eukaryotes, nuclear ribosomal DNA (nrDNA) is a cluster structure composed of multiple tandem transcription units, and each rDNA transcriptional unit consists of three rRNA coding sequences (the 18S rDNA gene, 5.8S rDNA gene and 28S rDNA gene), two internal transcribed spacers (ITSs) and one intergenic spacer (IGS) [6] [7] [8]. Due to different evolutionary rates, different regions of ribosomal transcription units can be regarded as molecular markers for species identification and phylogenetic studies at different levels in marine algal research. Reference [9] combined 18S rDNA and rbcL to analyse the phylogeny of Chloromonas and Chlamydomonas (Chlorophyceae, Volvocales), with an emphasis on snow and other cold-temperature habitats. Reference [10] demonstrated that Phytophthora de Bary and Peronosporales were a common natural group in phylogenetic studies based on nuclear large subunit ribosomal DNA sequences. Reference [11] combined ITS data, including the 5.8S rDNA sequence, with morphological data and found that the bloom-forming algae in the Yellow Sea in 2009 and 2010 were the same species. Reference [12] used the intergenic spacer (IGS) region to analyse the genetic variability of Pyropia haitanensis in different areas.
Similarity analysis and multisequencing alignment of the sequences indicated that the IGS sequences of different regions of P. haitanensis had notable variability.
Thus, the sequences of ribosomal transcription units are of great significance in phylogenetic research on marine algae.
This article provides additional information on the classification and identification of green algae. In addition, the methods outlined in this article can be used to obtain the entire ribosomal RNA gene (rDNA) sequence in other green algae taxa.

Sample Collection
Samples of Blidingia minima were collected from ten different sites in the Yellow Sea of China (Table 1). Once the samples were obtained at each site, they were immediately washed with seawater and dried in the shade to a moisture content of 30% -40%. Finally, they were transported to the laboratory in an insulated specimen box at 4˚C. All samples used in this experiment were stored in the Algae Laboratory of the School of Basic Medical and Biological Sciences, Soochow University, and when fresh algae were obtained, the primary axis, diverged branches, colours and basal part (holdfasts) were recorded and photographed.
The samples were sectioned into thin slices with a knife for observation under a microscope. Then, the cell shapes, sizes and arrangements were observed horizontally and vertically; the cell interior structures, such as chromatophores and pyrenoids, were also measured. A microscope (Nikon N90i) and a stereomicroscope (NiKon SMZ 1500) with an image collection system were applied in this research to obtain the main characteristics of the samples ( Figure 2).

Total RNA Extraction
Fresh algae samples were ground to a powder using an RNase-free mortar and pestle with liquid nitrogen, followed by the addition of 1 ml of Total RNA Extractor (Sangon Biotech, Shanghai) to completely cover the powder, according to the manufacturer's instructions. The concentration and quality of the extracted RNA were preliminarily tested by electrophoresis on 1% agarose gels. RNA samples were stored at −80˚C until use.

DNA Extraction
The algae were removed from the low-temperature refrigerator and placed in sterile double-distilled water for several hours. The algal surface was constantly brushed to remove the debris on the surface of the leaf thallus. The algae were then soaked in 0.7% KI for 10 min and rinsed several times with sterile seawater. Total genomic DNA was extracted by using a Plant Genomic DNA Kit (TIANGEN, BIOTECH (BEIJING) CO., Ltd) according to the instruction manual.

3' RACE cDNA Preparation
First, a poly(A) tail was added to the 3' end of the RNA using poly(A) polymerase (Takara Biotechnology (Dalian) Co., Ltd). Then, 11 μL of RNA was combined with 1 μL of 3'-CDS primer A, heated to 72˚C for 3 min and cooled to 42˚C for 2 min. Next SMARTScribe Reverse Transcriptase were added to each reaction tube according to the instructions of the SMART RACE Kit (Clontech). The 20 μL reaction was incubated at 42˚C for 90 min and then cooled to 70˚C for 10 min. The cDNA products were diluted with 90 μL of EDTA buffer.

PCR Amplification and Sequencing
An overview of the sequence cloning method of the ribosomal RNA gene in Blidingia minima is presented in Figure 3. The primers ( Table 2)  The PCR amplification cycle is presented in Table 3. The colony PCR conditions were as follows: 3 min initial denaturation at 95˚C; 30 cycles of denaturation at

Sequence Alignment and Analysis
Nucleotide BLAST was used for sequence alignment with the database to confirm the amplification results. Tandem repeats within the IGS sequence were identified using the Tandem Repeat Finder [14]. Multiple sequence alignment analysis was carried out using the DNAMAN software. Phylogenetic analysis was conducted using MEGA version 5 software [15]. Sequences were aligned using the CLUSTALW program. Bootstrap values obtained after 1000 replications are shown on the branches. The numbers around the branches indicate bootstrap supports.

Amplification of the 18S rDNA, ITS + 5.8S rDNA, 28S rDNA and IGS Sequences
Total genomic DNA or cDNA was used to amplify the fragments of ribosomal RNA in Blidingia minima with the corresponding primers. The 18S rDNA was amplified with two pairs of primers (Figure 4

Base Composition of the 18S rDNA, ITS + 5.8S rDNA, 28S rDNA and IGS Sequences
The base composition and GC content of each rDNA fragment were analysed, and the results showed that the C-base contents of the ITS and IGS sequences were the highest, and the G-base contents of the 18S rDNA and 28S rDNA sequences were the highest. The GC contents of the 18S rDNA, ITS + 5.8S rDNA, 28S rDNA and IGS sequences were 49.57%, 62.84%, 51.55% and 53.08%, respectively (Table 4).

Phylogenetic Tree Analysis
Maximum likelihood phylogenetic trees were established based on the ITS sequences ( Figure 5). Comparative analysis of the ITS sequences revealed that our amplified species were classified as Blidingia minima and distinguished from other green algae.
Maximum likelihood phylogenetic trees were established based on the IGS sequences ( Figure 6) showed that the strains were divided into two big clades. The ZQ1 and ZQ2 strains formed one clade, other eight strains formed another big clade. This result shows that IGS is suitable for intraspecies relationship analysis.   Table 1. Open Journal of Marine Science

Analysis of Special Structures in the IGS
Two tandem repeats of different lengths at the 5' end of the IGS sequence were found by the Tandem Repeat Finder. In addition, short dyad and palindromic sequences were found in the 3' end and in the middle of the IGS sequence. The locations of specific structures in the IGS sequence are shown in Figure 7, and the detailed characteristics of these special structures are presented in Table 5.

IGS Sequence Similarity Analysis
The lengths of the IGS sequences of 10 strains were all about the same ( Figure   8), and the sequences were relatively conserved. Multiple alignment of the IGS sequences of Blidingia minima from 10 different strains (Figure 9) showed multiple base differences, which indicated that the IGS sequences of the strains in different provinces were different, while the strains of the same province had only a few base differences.     Table 1.

Discussion
In eukaryotic organisms, the 18S rDNA, 5.8S rDNA, and 28S rDNA cluster as one transcript and are separated from the ITS and IGS. Comparing their full-length ribosomal RNA genes revealed that the ribosomal RNA gene of Blidingia minima was much shorter than those of Pyropia yezoensis and Bangia, and the sequence length of each region was also very different. In the process of amplifying the ribosomal RNA gene, we could search for only a partial sequence of the 28S 5' end of the nucleic acid sequence in the NCBI library, and it was very difficult to amplify the complete 28S sequence and IGS sequence with only one pair of primers. Therefore, we amplified the known partial sequence by PCR and then amplified the 28S 3' end sequence by rapid amplifica-Open Journal of Marine Science tion using cDNA ends technology. Finally, through sequencing and splicing, a 3282 bp 28S sequence was obtained. In the process of IGS sequence amplification, due to the presence of some specific structures in the IGS sequence, more specific primers and stringent reaction conditions were required in the PCR. Fortunately, we finally obtained a pair of efficient primers and identified the best reaction conditions. The GC content of each unit of the ribosomal RNA sequence unit was analysed and found to be highest in the ITS regions and lowest in the 18S rDNA genes. The total GC content of the ribosomal RNA gene was 54.26%, which was in the range of 35% -56% [19]. Given their high evolutionary rate, ITS sequences are expected to be suitable for detecting interspecific variation [20] [21]. We have made a phylogenetic tree of the amplified ITS sequences, and the results show that our amplified species are clustered together with Blidingia minima, whereas other algae are clustered into another branch, indicating that our expanded ITS can be used as an effective molecular marker to distinguish interspecies relations.
These regions can also contain functional sequences, such as spacers, promoters and enhancers, and transcription initiation and termination sites [28] [29]. Indeed, the IGS structure is attractive to study because of its importance in the transcription regulation of rRNA genes through the regulation of stop and start signals for transcription of the rDNA units that are located within the region [30] [31]. We compared the full sequence of IGS from 10 strains from different provinces. Strains in the same province were significantly different from IGS sequences in different provinces, which is consistent with the conclusion that the IGS is a highly variable region with a rapid evolutionary rate in Blidingia minima. Therefore, it can provide powerful classification and identification support at the species level and below.

Conclusion
In this paper, we first amplified the complete nuclear ribosomal DNA sequence of Blidingia minima and carried out a series of genetic analyses based on the ITS and IGS sequences. ITS can be effectively used as a molecular marker for interspecies identification. The phylogenetic tree shows that Blidingia minima can be distinguished from other green algae. We also used IGS to compare the homology of 10 different strains in the Yellow Sea. The results also show that IGS can provide important information for studying the intraspecific relationships and genetic diversity of Blidingia minima.