Journal of Oceanology and Limnology   2021, Vol. 39 issue(3): 1103-1117     PDF       
http://dx.doi.org/10.1007/s00343-020-0184-5
Institute of Oceanology, Chinese Academy of Sciences
0

Article Information

WANG Lu, XU Kuipeng, TANG Xianghai, WANG Junhao, KONG Fanna, MAO Yunxiang
Construction of high-density genetic linkage map of Pyropia yezoensis (Bangiales, Rhodophyta) and identification of red color trait QTLs in the thalli
Journal of Oceanology and Limnology, 39(3): 1103-1117
http://dx.doi.org/10.1007/s00343-020-0184-5

Article History

Received May. 4, 2020
accepted in principle Jun. 29, 2020
accepted for publication Jul. 6, 2020
Construction of high-density genetic linkage map of Pyropia yezoensis (Bangiales, Rhodophyta) and identification of red color trait QTLs in the thalli
Lu WANG1,3#, Kuipeng XU1#, Xianghai TANG1, Junhao WANG1, Fanna KONG1, Yunxiang MAO2,3     
1 Key Laboratory of Marine Genetics and Breeding(Ministry of Education), College of Marine Life Sciences, Ocean University of China, Qingdao 266003, China;
2 Key Laboratory of Utilization and Conservation of Tropical Marine Bioresource(Ministry of Education), College of Fisheries and Life Science, Hainan Tropical Ocean University, Sanya 572022, China;
3 Marine Biology and Biotechnology Laboratory, Pilot National Laboratory for Marine Science and Technology(Qingdao), Qingdao 266071, China
Abstract: Pyropia yezoensis is an important macroalga because of its extensive global distribution and economic importance. Color is an important trait in the thalli of P. yezoensis, it is also an effective marker to identify the hybridization in genetic breeding. In this study, a high-density genetic linkage map was constructed based on high-throughput single nucleotide polymorphism (SNP) markers, and used for analyzing the quantitative trait loci (QTLs) of red color trait in the thalli of P. yezoensis. The conchospore undergoes meiosis to develop into an ordered tetrad, and each cell has a haploid phenotype and can grow into a single individual. Based on this theory, F1 haploid population was used as the mapping population. The map included 531 SNP markers, 394.57 cM long on average distance of 0.74 cM. Collinear analysis of the genetic linkage map and the physical map indicated that the coverage between the two maps was 79.42%. Furthermore, QTL mapping identified six QTLs for the chromosomal regions associated with the red color trait of the thalli. The value of phenotypic variance explained (PVE) by an individual QTL ranged from 4.71%-63.11%. And QTL qRed-1-1, with a PVE of 63.11%, was considered the major QTL. Thus, these data may provide a platform for gene and QTL fine mapping, and marker-assisted breeding in P. yezoensis in the future.
Keywords: Pyropia yezoensis    high-density genetic linkage map    quantitative trait loci (QTL) mapping    F1 haploid population    red pigment variant    
1 INTRODUCTION

The genus Pyropia / Porphyra (Bangiales, Rhodophyta), popularly known as nori and laver, is one of the most economically important seaweeds in the world (Sutherland et al., 2011). According to FAO's statistics, in 2017, the global nori production was about 2.56 million tons in fresh weight, with a commercial value of more than $2.3 billion (http://www.fao.org/fishery/factsheets/en). In China, Pyropia yezoensis and Pyropia haitanensis are the two most commercially valuable Pyropia species. P. haitanensis is mainly cultivated in the Southern coastline of China. P. yezoensis is mainly cultivated in the Northern coastline of China and also cultivated in Korea and Japan. The life cycle of P. yezoensis consists of two alternative phases, macroscopic leafy thallus and microscopic filamentous conchocelis. The conchocelis, also called the sporophyte, is diploid and can release conchospores, which undergo meiosis and germinate into haploid thalli. It has been reported that P. yezoensis has the chromosome numbers of 2n=6 (Yabu and Tokida, 1963; Migita, 1967; Yabu, 1969) and the nuclear genome size is fewer than 100 Mb (Nakamura et al., 2013).

The selection and breeding of new varieties is important for developing economical seaweeds for the seaweed industry. The cultivation and breeding of laver originated in the 1950s (Fei et al., 1998). Traditional breeding of laver in China was based on selection, mutation, and crossing (Hu2006, 2006). Many chemical and physical techniques including application of chemical mutagens (Yan et al., 2000, 2004), γ rays (Wang et al., 2000; Yan et al., 2005), and heavy-ion beam (Niwa et al., 2011) was used to produce mutants of Pyropia, most of them were pigment mutations. In recent years, crossbreeding has been widely performed to produce new varieties with new traits. For example, Pyropia sp. from India was crossed with P. haitanensis to produce an improved strain (HR-5), which was characterized by fast growth, high quality, and production of large amount of conchospores (Ding et al., 2018). The intraspecific hybridization of the greentype pigmentation mutant and red-type pigmentation mutant of P. yezoensis was carried out, and then, an improved strain (A-18) with fast growth rate and wildtype color was isolated from the recombinants (Jiang et al., 2018). However, these traditional methods were limiting because of the time-consuming, inefficient, and difficult to distinguish phenotype for the blade plasticity for environmental factors. They were the bottlenecks in the development of laver industry. With the development of molecular breeding, molecular marker-assisted selection (MAS), and whole genome selection breeding, which have become the main tools providing direction for plant breeding enabling plant breeders to directly select the genotypes, thereby accelerating crop improvement (Johnson, 2003; Moose and Mumm, 2008; Ribaut et al., 2010; Yu et al., 2011; Qing et al., 2019; Randhawa et al., 2019; Wang et al., 2019). High resolution genetic linkage map has been the basis for these breeding strategies.

Genetic linkage maps have been essential for the genetics studies, genomic evolution, and mapping quantitative trait locus (QTL) (Cheng et al., 2019). To date, genetic linkage maps have been successively constructed for dozens of different species of plants and animals. In many crops, genetic linkage maps were also used to identify the QTLs associated with economic traits, for example plant height (PH), leaf length (LL), maize chlorotic mottle virus (MCMV) resistance, flag leaf-related traits, grain shape and salinity tolerance (De Leon et al., 2016; Liu et al., 2018; Chen et al., 2019; Choi et al., 2019; Sitonik et al., 2019). As compared to the land crops, the development of genetic linkage maps in seaweeds is relatively lacking. The genetic linkage maps have been reported for only few species of seaweeds. A tentative amplified fragment length polymorphismsimple sequence repeat (AFLP-SSR) linkage map of Laminaria was constructed with an average distance of 7.91 cM (Yang et al., 2009). In Saccharina japonica, based on the specific length amplified fragment sequencing (SLAF-seq), 7 627 single nucleotide polymorphism (SNP) loci were mapped to 31 linkage groups with an average spacing of 0.69 cM, and the QTLs for the blade length and blade width phenotypes were detected (Wang et al., 2018b). A genetic map of Ectocarpus siliculosus was constructed containing 406 AFLP markers resulting in 34 linkage groups, which spanned a total distance of 2 459.3 cM and an average distance of 6.6 cM (Heesch et al., 2010). In Undaria pinnatifida, a highdensity genetic linkage map was constructed using SLAF-seq technique, which spanned a distance of 1 816.28 cM and consisted of 30 linkage groups, with an average distance of 0.39 cM (Shan et al., 2015). In Pyropia, the first high-density genetic linkage map was constructed in P. haitanensis using SLAF-seq, which included 4 550 length polymorphic markers. These markers were combined into 740 bins on five linkage groups, with a length of 874.33 cM and an average distance of 1.18 cM between the adjacent bins (Xu et al., 2015). Recently, a genetic map of P. yezoensis was constructed using only 92 sequence-related amplified polymorphism (SRAP) loci to produce three linkage groups (LGs), with a mean interlocus space of 6.23 cM (Huang and Yan, 2019).

Although two genetic linkage maps of Pyropia have been constructed, some problems still exist. Thus, it is still necessary to construct a genetic linkage map of Pyropia with high quality and density, especially for P. yezoensis. It was difficult to apply the molecular markers technique in laver breeding because the thallus of F1 gametophyte was chimera. The color mutants of the thalli were widespread used in the genetic studies of Pyropia/ Porphyra (Yan et al., 2000). In nature or mutant breeding of P. yezoensis, it is common to identify red pigment variants or mutants. A red pigment variant was found previously in the laboratory, and the contents of phycoerythrin and phycocyanin were significantly higher than that in the wild type. And we found that the red color trait probably was a qualitative trait, which could be used to distinguish genotypes of the tetrad. So far, there is no research on the red color trait of thalli yet. Therefore, in order to apply this red color trait as a marker in molecular pyramiding breeding, the QTL mapping for the red color trait was imperative. In the present study, we constructed a high-density genetic linkage map of P. yezoensis using SNP markers identified by genome resequencing and F1 haploid population, and further detected the red color trait of thalli using the QTL method. Our study has provided a means to correlate the traits and genotypes in P. yezoensis and may contribute towards the ultimate goal of molecular marker-assisted breeding.

2 MATERIAL AND METHOD 2.1 Plant material

Two pure lines, RZ and HT, with different blade colors were used as the parent strains to construct the mapping populations in this study. Pure RZ line with a wild-type, established by our laboratory was the maternal parent (Fig. 1). Single somatic cells were enzymatically isolated from the thallus of RZ under the microscope with glass capillaries. Pure HT line with a red-type was the paternal parent and a spontaneous mutant which was isolated from a pure line PYL-349 (Fig. 1). Healthy and immature thalli of parents were separately cultured in 500-mL flasks containing PES medium (UTEX Culture Collection of Algae) and agitated by aeration at 15 ℃ under 80 μmol photons/(m2·s) (10-h L꞉14-h D) (Niwa, 2010).

Fig.1 The thallus of the male parental strain (RZ), the female parental strain (HT), and the F1 thallus selected from the cross between the RZ and HT strains
2.2 Construction of map population

The crossing was performed according to the method of Niwa (2010). The diagrammatic sketch is shown in Fig. 2. Two pure lines, HT line and RZ line, were used for the cross to construct the mapping population. The selected male and female thalli were co-cultured at a ratio of 3꞉1 for about 10 days until carposporangia appeared in the maternal piece. Then the maternal blade pieces were washed and gently rubbed in a petri dish with sterilized seawater which were both filtered and boiled. The cleaned pieces were chopped with a blade, then, the carpospores were sucked out with glass micropipettes (inside diameter, 0.5 mm) and cultured in 96-well plates to grow into F1 conchocelis colonies at 20 ℃ under 30 μmol photons/(m2·s) (12-h L꞉12-h D). When the conchocelis colonies grew to a certain size, they were fragmented to inoculate the shells, and continually cultured at 24 ℃ under 30 μmol photons/(m2·s) (10-h L꞉14-h D) for about a month and a half until the conchosporangia formed. Once the conchospores were released and attached to the cotton strings at 20±1 ℃ under 135 μmol photons/(m2·s) (12-h L: 12-h D), they were transferred to a flask containing 800 mL of PES culture medium and cultured with aeration at 20±1 ℃ under 135 μmol photons/(m2·s) (12-h L꞉12-h D) to develop into thalli. The culture medium was renewed every 3 days. After 2 weeks, the temperature was reduced to 10±1 ℃ and the culture was continued for one and a half months under the conditions of 160 μmol photons/(m2·s) (12-h L꞉12-h D) until the thalli were obtained. The F1 thalli with only four color-sectors were selected and each color-sector blade was singly developed into an individual. These thalli were obtained to form the F1 haploid population.

Fig.2 Schematic view of the cross between RZ and HT, and construction of F1 haploid population of P. yezoensis
2.3 DNA extraction and Illumina library construction

DNA was isolated from the thallus of each parental line and F1 haploid lines. The collected blades were ground into powder in liquid nitrogen, and the DNA was extracted and purified using the Plant Genomic DNA Kit (DP305-02, TIANGEN) following the manufacturer's method. The DNA quality was determined using a Nanophotometer N60 (Implen, Germany) and the concentration was determined based on the Qubit2.0. The size and integrity of DNA were assessed by 1.0% agarose gel electrophoresis using lambda DNA as the standard. After qualified detection, DNA samples were used for Illumina library construction using the VAHTS Universal DNA Library Prep Kit for Illumina® V3 ND607 (Vazyme Biotech Co., Ltd; Nanjing, China) according to the manufacturer's protocol.

2.4 High-throughput sequencing and SNP genotyping

The eligible libraries were diluted for pair-end sequencing on Illumina HiseqTM Xten sequencing platform (Illumina, Inc; SanDiego, CA, USA) at Novogene Bioinformatics Technology Co. Ltd. in Tianjin, China. Data quality control was performed based on the obtained resequencing data considering the following four criteria: i) removal of adapter contaminated reads; ii) the low quality (bases with a mass value less than 5 and accounting for more than 20% of the total reads) reads were dropped; iii) removal of reads with a high ratio of N (the base measured as N accounts for more than 10% of the total reads); and iv) avoiding duplication.

The clean reads of the parents and F1 individuals were first mapped to the reference genome (WMLA00000000) with BWA backtrack (Li and Durbin, 2009). The reference genome was indexed and the command 'aln' was used with the parameters '- o 1 - e 10 - t 4 - l 32 - i 15 - q 10' to find the suffix array coordinates of good matches for each read (Xu et al., 2019). The best alignments were then converted into the Sequence Alignment/Map (SAM) format using the 'sample' command (Li and Durbin, 2009). We improved the accuracy of sequence alignment using the following filtering criteria: i) removal of the number of mismatched bases less than 5 and the alignment result with a quality of 0; ii) use of Picard software (McKenna et al., 2010) with the command 'AddOrReplaceReadGroups' to get a BAM file, and used command 'FixMateInformation' to ensure each read and its mate pair were synchronized; iii) Removal of potential PCR duplicates.

After alignment, used the SAMtools program (Li et al., 2009) with Bayesian algorithm, SNP calling was performed on a population scale. The command 'mpileup' was used with the parameters '- q 1 - C 50 - S - D - m 2 - F 0.002 - u' (Xu et al., 2019). The subsequent analysis filtered the raw SNPs using Variant Call Format (VCF) tools (Danecek et al., 2011) and (Genome Analysis Toolkit) GATK (McKenna et al., 2010) using the following criteria: i) coverage depth ≥4, ≤1 000; ii) mapping quality ≥20; iii) distance of adjacent SNPs ≥5 bp; iv) retention of only those SNPs that occurred in more than 50% of the individuals; v) minor allele frequency (MAF) ≥0.05.

During data processing for genotyping, since the genotype of F1 individual was haploid, they were artificially doubled into diploid. The SNPs between the parents were considered as the original markers for genotyping of 131 individuals, and only the individual with the segregation patterns "aabb" were retained for genetic linkage map construction. Genotyping was defined by covering at least 80% of all individuals. It means that at least 105 individuals out of 131 individuals can be detected genotypes (a/b) at a site.

2.5 Genetic linkage map construction

The genetic linkage map was constructed by Joinmap 4.1 (Cheema and Dicks, 2009). Marker segregation ratios were calculated using the chisquare test, and the markers showing significant segregation distortion (P < 0.001) were excluded. Logarithm of the odds (LOD) value ≥5.0 was used as thresholds for groups. The map distances were calculated by Kosambi mapping function. Genetic maps were drawn using lepmap3.0 software.

2.6 Evaluation of the genetic map

In our group, high-throughput chromosome conformation capture (Hi-C) technique was used to assemble the genome of P. yezoensis at the chromosome level. Three scaffolds were obtained that were equivalent to the three chromosomes observed in haploid P. yezoensis (Ohme and Miura, 1988). The quality of the genetic map was evaluated with heat maps generated for each LG, and the collinear analysis between the genetic linkage map and the physical map was performed. The pair-wise recombination values of all the mapped SNP markers were used to assess the quality of the genetic linkage map in heat maps. Furthermore, heat maps were also to reflect the relationship of recombination between the markers from one single LG (Shan et al., 2015). Collinear analysis reflected the positional relationship of markers between the genetic linkage map and the physical map. Both drawings were implemented in R language.

2.7 Quantitative trait locus (QTL) analyses of red color trait

The phenotypes of F1 individuals for color as the criterion were statistically analyzed; the red color was denoted as R, and the normal color was denoted as W. MapQTL6.0 program was used to determine the phenotypic thresholds, and the significant LOD thresholds for the red color trait were calculated by the permutation test of α < 0.05 and n=1 000 for significant linkages (Van Ooijen and Kyazma, 2009). Based on these permutations, QTL analysis was carried out using WinQTLcart2.5 (Wang et al., 2007), and the composite interval mapping (CIM) method (Zeng, 1994) was performed to detect any significant associations between the red color trait and marker loci. A LOD score of 3.0 at genome level was used as a minimum cutoff to declare the presence of a QTL. The value of phenotypic variance explained (PVE) of a major QTL is defined as more than 15% (Kumawat et al., 2012; Liu et al., 2014; Wang et al., 2018a).

3 RESULT 3.1 The construction of F1 haploid population

Based on the scheme in Fig. 2, pure line RZ (♀) was crossed with HT (♂) and carpospores were released from fertilized carposporangia on the maternal blade piece, which developed into conchocelis. Then the conchocelis colonies were further released conchospores, which grew to F1 thalli. Among the F1 thalli that were obtained, only the chimeras of color-sectors with two or three or four color-sectors (Fig. 2) showed that the crossing process was successfully. In this research, only the thalli with tetrad chimera were selected for constructing population. Finally, 152 color-sectors from 38 F1 thalli were obtained as mapping population.

3.2 SNP calling and genotyping

A total of 152 color-sectors were used to build the DNA libraries, and only 131 were successful. Two parents and 131 individuals were used to generate SNP data for subsequent analysis, and 96% of the sequences were correctly aligned to the reference genome of P. yezoensis. In total, 24.24 Gb clean data were generated for the parents, 8.81 Gb for the female parent (RZ), and 15.43 Gb for the male parent (HT). The average sequencing depth was 75.89× in RZ and 144.21× in HT. For 131 individuals, a total of 676.59 Gb clean data were obtained containing 43.57 hundred million reads, on average of 150 bp in length for each read. The number of reads per single individual ranged from 0.99–75.35 million, with an average of 33 million reads and the sequencing depth of 46.27×.

A total of 131 F1 haploid individuals from colorsectors were used to construct the mapping population. For each blade color-sector, its genotype was haploid. It could develop into double haploid conchocelis by self-cross, with a homozygous genotype. Thus, in this study, to reduce the time and to improve the efficiency, we did not complete the self-crossing of the F1 blade color-sectors. Before genotyping, the genotypes of F1 color-sector individuals were artificially directed as doubled into diploid. Then a F1 haploid population of 131 individuals was constructed and used for genotyping and constructing the genetic linkage map. The F1 haploid population was genotyped by 24 438 SNP polymorphic markers between the two parents, and finally 15 822 polymorphic marker loci were successfully genotyped in 131 individuals (Supplementary Table S1). These markers showed the aa×bb segregation pattern. No abnormal bases were found, indicating that the accuracy of genotyping was high. After filtering out the markers with less than 80% genotypic integrity coverage, 668 markers were finally obtained. The segregation distortion threshold p was set to be 0.001, and only three distorted segregation markers were filtered out. Ultimately, a total of 665 valid markers were entered into the linkage analysis.

3.3 High-resolution genetic linkage map construction

After linkage analysis, 531 of the 665 markers were successfully mapped onto the genetic linkage map, while other 134 markers failed to link to any group. The 531 markers were distributed on the three linkage groups (Supplementary Table S2). The final genetic linkage map was 394.57 cM in length with an average distance of 0.74 cM between the adjacent markers (Fig. 3, Table 1). The number of LGs in the genetic linkage map was same as the chromosome number of the haploid genome (3 chromosomes) of P. yezoensis (Yabu and Tokida, 1963; Migita, 1967; Yabu, 1969). As shown in Table 1, the length of the three LGs was 176.22, 141.35, and 77.00 cM, respectively (mean, 131.52 cM). The degree of linkage between the markers was reflected by "Gaps < 5 cM", which ranged between 97.436% and 99.425%, with an average value of 98.542%. Among the three LGs, LG1 was the largest LG, including 162 SNP markers and spanning 176.22 cM with an average marker interval of 1.09 cM. The second largest group was LG2 with 195 SNP markers and spanning 141.35 cM. The average inter-marker distance for LG2 was 0.72 cM. LG3 was the most saturated and the shortest, containing 174 SNP markers spanning 77.00 cM with an average marker interval of 0.44 cM.

Fig.3 High-density linkage map for P. yezoensis Marker loci are shown on the right and locations are shown on the left. The markers labelled with red color are the anchor markers. The specific details of the linkage map are provided in Supplementary Table S2.
Table 1 Summary statistics of the genetic linkage map of P. yezoensis
3.4 Evaluation of the genetic map

Collinearity analysis was performed by mapping the markers on the genetic linkage map against the scaffolds, as shown in Fig. 4. In general, three linkage groups showed high collinearity corresponding to the scaffolds. For LG1 and scaffold_1, the region before marker mk37 at the 120 cM could correspond to the physical region from 1 to 21 617 320 bp on scaffold_1 in turn; while an inversion occurred behind the region of 120 cM on LG1 and the region from 25 184 629 bp to the end of scaffold_1. And the inverted region in LG1 contained a total of 35 markers. For LG2 and scaffold_2, most of the marker positions except the region from 66 cM to 91 cM which contained 76 markers could correspond to scaffold_2 as well. Similarly, the marker positions of LG3 could correspond to scaffold_3 consistently. Overall, the coverage between the physical and genetic linkage maps was 79.42%.

Fig.4 Collinearity analysis of genetic map and physical map in P. yezoensis Genetic maps are shown on the left (red color) and physical maps are shown on the right (dark blue color). The light blue lines in the middle reflect the corresponding position relationship between the two maps.

The heat map reflects the recombination relationship between the markers from one single linkage group, which can be used to identify errors in the order. As shown in Supplementary Figs.S1–S3, the markers on each linkage group were used to calculate pair-wise recombination values using three heat maps. If the order between the markers was correct, the color on the diagonal was the darkest and reddest, indicating that the linkage relationship was the strongest. The linkage relationship between the markers became weaker as the distance increased. In general, the visualization of the heat map indicated that the LGs performed well.

3.5 Identification of red color related QTL

The gametophytic phenotype of parent RZ is naturally brown. The paternal parent HT is the natural variant from RZ, with a red color gametophytic phenotype. The phenotypes of blade color trait, of the parents and the 131 F1 singles were recorded as the phenotype data of the QTLs (Supplementary Table S3). Among the 131 F1 singles, there were 70 wildtype (R) and 61 red-type (H) (Fig. 5). The statistical significance of the QTL effects was evaluated by 1 000 permutations, and the threshold at genomewide significance level 0.05 of the LOD scores was determined as 3.0. According to the threshold at genome level, 6 QTLs associated with the red color trait were identified on LG1 and LG2; no QTLs were found on LG3 (Table 2). Among the 6 QTLs, the minimum and maximum LOD scores were 3.194 and 27.157 respectively, and the value of PVE ranged from 4.71%–63.11% (Table 2).

Fig.5 The distribution of blade color trait in the 131 F1 singles R: represents the blade color trait is wild-type; H: represents the blade color trait is red-type.
Table 2 QTLs detected for red color trait based on CIM analyses in the F1 haploid population

The QTL qRed-1-1 was located at 49.9–50.3 cM between mk159 and mk162 on LG1, with the maximum PVE value of 63.11%. This locus explained 63.11% of the phenotypic variance, indicating that it may be a major locus for the red color trait. Based on the genome data of this region, from 37 302 208 bp to 41 082 413 bp on scaffold_1, a total of 501 genes were identified, 282 were annotated in the nonredundant (NR) database (Supplementary Table S4) and 219 genes with unknown functions. The Gene Ontology (GO) enrichment analysis of these 282 genes showed that their functions including enzymatic activity, basal metabolic process, photomorphogenesis, and so on (Supplementary Table S5). Subsequent fine mapping and gene editing techniques are needed for further verification. The QTL qRed-1-2 was located at 59.1–60.4 cM, with the physical position between 21 327 428 bp to 21 420 62 bp of scaffold_1, explaining 14.83% of the phenotypic variance. The QTL qRed-2-1 was located on LG2 at position 22.6– 23.6 cM, with the physical position between 32 143 876 bp to 32 430 593 bp of scaffold_2, explaining 10.22% of the phenotypic variance. The QTL qRed-2-2 was located on LG2 at position 34.4– 36.3 cM, with the physical position between 30 634 713 bp to 32 430 699 bp of scaffold_2, explaining 14.09% of the phenotypic variance. In addition, the other two QTLs (qRed-2-3 and qRed-1-3) were located at 40.0–45.2 cM and 74.0–79.1 cM along LG2 and LG1, with PVE of 8.12% and 4.71%, respectively. The CI (95% confidence interval for QTL length) ranged from 0.4–5.1 cM among the 3 QTLs on LG1, and the CI ranged from 1.0–5.2 cM among the 3 QTLs on LG2, and half of them were less than 1.0 cM. Therefore, 3 QTLs with a distance of less than 1.0 cM from the nearest markers could be used for further studies.

4 DISCUSSION 4.1 Mapping population of P. yezoensis

The construction of population map is the first step in constructing a genetic linkage map. Currently, among the commonly used mapping populations, such as backcross 1 (BC1), F2 generation population, double haploid (DH) population, and recombinant inbred lines (RIL), the DH population is considered as the ideal mapping population because of its homozygous genotype. In plants, there are various methods to construct DH population, such as anther culture, haploid induction, and natural doubling (Zhu et al., 2004; Chai et al., 2007; Ma et al., 2018; Shen et al., 2018). For P. haitanensis, its gametophytic single cell in vitro can directly develop into conchocelis and the genotype can be doubled into diploid. Thus, it is easy to construct a DH population for P. haitanensis after obtaining F1 thallus. A DH population of 100 lines was built to construct the genetic linkage map of P. haitanensis (Xu et al., 2015). However, P. yezoensis does not own this character. Its thalli are monoecious and could be self-fertilized. Thus, it can complete double haploid by self-fertilization.

In addition, F1 conchospores of Pyropia divided into ordered tetrads via meiosis, which were chimeras. If there were no obvious markers, it was difficult to select the heterozygous F1 thallus to ensure if the crossing occurred. The color sector embedded in the F1 heterozygous thallus is an ideal selection marker to resolve this problem (Yan and Aruga, 2000). The ordered tetrad continued to divide and develop into 2–4 pigment containing chimeras. In this study, we used red color pigment variant to make a cross with the wild type, and then F1 thalli with color-sectors were used to construct the mapping population (Fig. 2).

Since each color-sector of the thalli is haploid, after the thalli self-fertilized to form carpogonium and developed into conchocelis, their genotypes were doubled. Huang and Yan (2019) used this strategy to construct a genetic linkage map of P. yezoensis, but the process was time-consuming. Conchospore germination of P. yezoensis developed in a linear pattern, and cell divisions were not synchronous (He et al., 2019). In this study, although we selected F1 haploid population as the mapping population to construct the genetic map, yet during genotyping, the genotypes of F1(a/b) individuals were automatically doubled and changed into aa/bb. We confirmed that this strategy is very scientifically feasible.

4.2 Features of the high-density map of P. yezoensis

As important macroalgae in aquaculture, a highdensity genetic linkage map is essential for molecular marker-assisted selection and QTL mapping of P. yezoensis. The acquisition of high-density markers is critical to increase the density of genetic linkage maps. High-throughput sequencing provides support for the development of high-throughput SNP molecular markers at the whole genome level, including SLAF, Restriction-site associated DNA sequencing (RAD-seq), IIb restriction site-associated DNA (2b-RAD) and re-sequencing (Davey et al., 2011; Sun et al., 2013). As a result, reducedrepresentation genome sequencing (RRGS) technology has also emerged. Although SLAF-seq was used to develop long polymorphic markers in P. haitanensis, but most of them were found to be redundant (Xu et al., 2015). This may be due to the peculiarity of Pyropia genome as lots of microorganisms and Pyropia thalli showed serious symbiosis relationship. When isolating the genome of Pyropia, most of the nucleic acids are derived from the microorganisms. In addition, the genome of plastids and mitochondria also account for a very high ratio. Among the 1748 polymorphic markers from SLAF-seq, only two markers could be used for constructing the genetic map (Xu et al., 2015). Then, the bin-marker approach was further utilized to construct the final genetic map, but it still lacked genomic information for reference. Huang and Yan (2019) used fluorescent SRAP markers to construct the genetic linkage map of P. yezoensis; this technique is based on PCR technology, which inevitably limits the density of the map and the throughput of the molecular markers. The map included 92 SRAP markers (Huang and Yan, 2019); however, the SRAP primers were designed from higher plants (Dendrobium loddigesii Rolfe, Eucommia ulmoides Oliver and Malus Mill), and most of the loci were located in the open reading frame (ORF) regions of the gene. In addition, each marker listed on their map was one of the bands of SRAP amplification products and its sequence was not provided in this paper. Therefore, it was difficult to determine the locations of SRAP markers on the genome and integrate our map with their map.

In this study, we applied re-sequencing based on the reference genome (WMLA00000000) with superior quality to develop SNP markers at the genomic level. All the contaminating sequences from the bacteria, plastids, and mitochondrial genome were removed from the genome sequences. Using this high-quality reference genome, more effective SNP markers were mapped to the genome to increase the density of the markers. As compared to the published genetic map of Pyropia, such as P. haitanensis, with an average density of 1.18 cM and P. yezoensis with an average density of 6.23 cM (Xu et al., 2015; Huang and Yan, 2019), the mapping in our study with an average density of 0.74 cM is the highest (Fig. 3, Table 1). For mapping, 668 markers out of 15 822 polymorphic markers were finally obtained based on 80% genotypic integrity coverage. After filtering out three distorted segregation markers, a total of 665 valid markers were entered into the linkage analysis. The rest of the 15 154 markers showed no recombination and were found to be redundant. Although our map has the highest density as compared with other Pyropia, the density of our map can be further improved by increasing the parental polymorphism and the number of population individuals.

High-resolution genetic mapping with many markers is helpful in improving the genome sequence assembly in plants (Ganal et al., 2011). The number of LG in our study was 3 (Table 1), which is equivalent to the haploid chromosome number of P. yezoensis (Ohme and Miura, 1988). This indicated that the mapping strategy in this research was scientific. Most regions of the genetic linkage map showed high linearity with the scaffolds of the reference sequence of P. yezoensis (Fig. 4). Some short inversion and translocation phenomena existed in the partial region (Fig. 4). Given that the reference genome comes from RZ pure line, the occurrence of the inversion and translocation may be because the genome of pure line HT was used in crossing. Thus, this phenomenon has also provided a guide to confirm and improve the assembly of genome sequence data.

4.3 QTLs of red color traits

In nature or mutant breeding, it is easy to identify color variants or mutants. Color variants or mutants play an important role in the genetic studies of Pyropia / Porphyra (Miura and Shin, 1989). Ohme et al. (1986) used the mutant pigment color as the marker to clarify meiosis as the dividing process for the conchospore. Ohme and Miura (1988) also used tetrad analysis to determine the distance between the centromere and the loci of four color mutants that were assigned to three different linkage groups. Color variants or mutants have also been used as genetic markers for intraspecific or interspecific crossbreeding of Pyropia to produce improved lines (Yan and Aruga, 2000; Liu et al., 2015). Although color mutants have also been used to construct genetic linkage maps (Huang and Yan, 2019), there are no reports regarding their genetic location on the chromosomes.

We first identified 6 QTLs distributed on 3 LGs that associated with red color trait, including one major QTL and five minor QTLs (Table 2). The PVE of a major QTL is normally defined as more than 15% (Kumawat et al., 2012; Liu et al., 2014; Wang et al., 2018a). Therefore, in this study, the QTLs with PVE larger than 15% were defined as the major QTLs. It was noted that major QTL qRed-1-1 showed a PVE of up to 63.11%, with the highest LOD value of 27.157 and the 95% CI of 49.9–50.3 cM (Table 2). The PVE of the remaining five QTLs were all less than 15%. Therefore, we believe that the contribution of these five minor QTLs to red color trait could be negligible, and QTL qRed-1-1 made an absolute contribution to the red color trait. Higher LOD values could better control the occurrence of false QTLs and are suitable for fine mapping of the major QTLs (Li et al., 2010; Tao et al., 2013). Thus, QTL qRed-1-1 should be further analyzed to identify the genes controlling the red color trait.

The color-related traits in higher plants showed that they were controlled by the single-gene or QTL regions through QTL analysis. The monogenic mutations lead to changes in the fruit color of pepper (Brand et al., 2012). The red color leaf of peach and the yellow color of Lophopyrum ponticum are also controlled by single gene, the Gr gene, and the yellow pigment gene Y, respectively (Dirlewanger et al., 2004; Zhang et al., 2005). Yet, the leaf color in Chinese cabbage (Brassica rapa L. subsp. pekinensis) and switchgrass (Panicum virgatum) have been detected to be controlled by the QTLs (Xu et al., 2007; Lowry et al., 2015). In pot azalea (Rhododendron simsii hybrids), pink petal color has been confirmed to result from a gene-dosage effect, and two major QTLs were found to be associated with the petal color (De Keyser et al., 2013). Candidate genes involved in skin and flesh coloration in peach (PprMYB10), cherry (PavMYB10), and apple (MdMYB1 / MdMYBA / MdMYB10) have also been reported (Kui et al., 2010; Sooriyapathirana et al., 2010; Frett et al., 2014). In the present study, the major QTL for the red color trait of blade, qRed-1-1 on LG1, existed on scaffold_1 from 37 302 208 bp to 41 082 413 bp (Supplementary Table S2), and a total of 501 genes were identified. According to the NR functional annotation, 282 genes were annotated with known functions (Supplementary Table S4). GO enrich analysis showed that the genes related to organelle organization accounted for the largest rate, and the genes related to catalytic activity and cellular components ranked second. To prove the authenticity of a QTL, it is necessary to convert a QTL from a statistical to a genetic trait (Mo, 2003). Further reducing the distances between the QTLs with red color trait by fine-mapping is necessary to identify the genes controlling this trait. The gene knock-out or genome editing technologies can also be used to validate the genes controlling the red color trait of Pyropia.

5 CONCLUSION

In this study, a F1 haploid mapping population consisting of 131 individuals was constructed based on the color-sectors F1 thallus from the cross between the pure line RZ and the red variant HT. By genome resquencing, a high-density genetic linkage map was constructed, containing three linkage groups with a total length of 394.57 cM and an average genetic distance of 0.74 cM. This map showed better collearity with the scaffold of P. yezoensis generally. Six QTLs related to blade red color trait were mapped to linkage groups 1 and 2. QTL qRed-1-1 with PVE of 63.11% was considered as the major QTL. These results will not only facilitate marker-assisted selection breeding of P. yezoensis, but also lay a foundation for further isolating of the genes control red color trait.

6 DATA AVAILABILITY STATEMENT

The data generated and/or analyzed during the current study are available from the corresponding author on reasonable request.

Electronic supplementary material

Supplementary material (Supplementary Tables S1–S5 and Figs.S1–S3) is available in the online version of this article at https://doi.org/10.1007/s00343-020-0184-5.

References
Brand A, Borovsky Y, Meir S, Rogachev I, Aharoni A, Paran I. 2012. pc8.1, a major QTL for pigment content in pepper fruit, is associated with variation in plastid compartment size. Planta, 235(3): 579-588. DOI:10.1007/s00425-011-1530-9
Chai L G, Zhang J J, Lin G P, Wang Y, Xu F S. 2007. Construction of two DH populations and identification of chromosome ploidy in burley tobacco. Acta Tabacaria Sinica, 13(2): 33-37. (in Chinese with English abstract) DOI:10.3321/j.issn:1004-5708.2007.02.007
Cheema J, Dicks J. 2009. Computational approaches and software tools for genetic linkage map estimation in plants. Briefings in Bioinformatics, 10(6): 595-608. DOI:10.1093/bib/bbp045
Chen J X, Zhou H, Gu Y, Xia D, Wu B, Gao G J, Zhang Q L, He Y Q. 2019. Mapping and verification of grain shape QTLs based on high-throughput SNP markers in rice. Molecular Breeding, 39(3): 42. DOI:10.1007/s11032-019-0955-x
Cheng L R, Chen X C, Jiang C H, Ma B, Ren M, Cheng Y Z, Liu D, Geng R M, Yang A G. 2019. High-density SNP genetic linkage map construction and quantitative trait locus mapping for resistance to cucumber mosaic virus in tobacco (Nicotiana tabacum L. ). The Crop Journal, 7(4): 539-547. DOI:10.1016/j.cj.2018.11.010
Choi J K, Sa K J, Park D H, Lim S E, Ryu S H, Park J Y, Park K J, Rhee H I, Lee M, Lee J K. 2019. Construction of genetic linkage map and identification of QTLs related to agronomic traits in DH population of maize (Zea mays L. ) using SSR markers. Genes & Genomics, 41(6): 667-678. DOI:10.1007/s13258-019-00813-x
Danecek P, Auton A, Abecasis G, Albers C A, Banks E, DePristo M A, Handsaker R E, Lunter G, Marth G T, Sherry S T, McVean G, Durbin R, 1000 Genomes Project Analysis Group. 2011. The variant call format and VCFtools. Bioinformatics, 27(15): 2 156-2 158. DOI:10.1093/bioinformatics/btr330
Davey J W, Hohenlohe P A, Etter P D, Boone J Q, Catchen J M, Blaxter M L. 2011. Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nature Reviews Genetics, 12(7): 499-510. DOI:10.1038/nrg3012
De Keyser E, Lootens P, Van Bockstaele E, De Riek J. 2013. Image analysis for QTL mapping of flower colour and leaf characteristics in pot azalea (Rhododendron simsii hybrids). Euphytica, 189(3): 445-460. DOI:10.1007/s10681-012-0809-7
De Leon T B, Linscombe S, Subudhi P K. 2016. Molecular dissection of seedling salinity tolerance in rice (Oryza sativa L. ) using a high-density GBS-based SNP linkage map. Rice, 9(1): 52. DOI:10.1186/s12284-016-0125-2
Ding H C, Lv F, Wu H X, Yan X H. 2018. Selection and characterization of an improved strain produced by inter-species hybridization between Pyropia sp. from India and Pyropia haitanensis from China. Aquaculture, 491: 177-183. DOI:10.1016/j.aquaculture.2018.03.015
Dirlewanger E, Cosson P, Howad W, Capdeville G, Bosselut N, Claverie M, Voisin R, Poizat C, Lafargue B, Baron O, Laigret F, Kleinhentz M, Arús P, Esmenjaud D. 2004. Microsatellite genetic linkage maps of myrobalan plum and an almond-peach hybrid-location of root-knot nematode resistance genes. Theoretical and Applied Genetics, 109(4): 827-838. DOI:10.1007/s00122-004-1694-9
Fei X G, Wang S J, Lu S. 1998. Red Algae Cell Engineering Breeding and Seedling. Marine Biotechnology. Shandong Science and Technology Press, Jinan, China. p. 7-26. (in Chinese)
Frett T J, Reighard G L, Okie W R, Gasic K. 2014. Mapping quantitative trait loci associated with blush in peach[Prunuspersica (L.) Batsch]. Tree Genetics & Genomes, 10(2): 367-381. DOI:10.1007/s11295-013-0692-y
Ganal M W, Durstewitz G, Polley A, Bérard A, Buckler E S, Charcosset A, Clarke J D, Graner E M, Hansen M, Joets J, Le Paslier M C, McMullen M D, Montalent P, Rose M, Schön C C, Sun Q, Walter H, Martin O C, Falque M. 2011. A large maize (Zea mays L. ) SNP genotyping array: development and germplasm genotyping, and genetic mapping to compare with the B73 reference genome. PLoS One, 6(12): e28334. DOI:10.1371/journal.pone.0028334
He Y, Zhou W, Li J M, Shen S D. 2019. Study of mosaicism of Pyropia yezoensis. Marine Biology Research, 15(2): 150-158. DOI:10.1080/17451000.2019.1605183
Heesch S, Cho G Y, Peters A F, Le Corguillé G, Falentin C, Boutet G, Coëdel S, Jubin C, Samson G, Corre E, Coelho S M, Cock J M. 2010. A sequence-tagged genetic map for the brown alga Ectocarpus siliculosus provides large-scale assembly of the genome sequence. New Phytologist, 188(1): 42-51. DOI:10.1111/j.1469-8137.2010.03273.x
Hu Y M. 2006. The history and actuality of Porphyra's breeding. Journal of Ningbo Institute of Education, 8(4): 22-26. (in Chinese with English abstract) DOI:10.3969/j.issn.1009-2560.2006.04.007
Huang L B, Yan X H. 2019. Construction of a genetic linkage map in Pyropia yezoensis (Bangiales, Rhodophyta) and QTL analysis of several economic traits of blades. PLoS One, 14(3): e0209128. DOI:10.1371/journal.pone.0209128
Jiang H, Ding H C, Yan X H. 2018. Selection and characterization of an improved strain (A-18) by hybridization recombinant in Pyropia yezoensis (Bangiales, Rhodophyta). Acta Oceanologica Sinica, 40(2): 95-103. (in Chinese with English abstract) DOI:10.3969/j.issn.0253-4193.2018.02.010
Johnson R. 2003. Marker-assisted selection. In: Janick J ed. Plant Breeding Reviews. John Wiley & Sons, New York, p. 293-309, https://doi.org/10.1002/9780470650240.ch13.
Kui L W, Bolitho K, Grafton K, Kortstee A, Karunairetnam S, McGhie T K, Espley R V, Hellens R P, Allan A C. 2010. An R2R3 MYB transcription factor associated with regulation of the anthocyanin biosynthetic pathway in Rosaceae. BMC Plant Biology, 10(1): 50. DOI:10.1186/1471-2229-10-50
Kumawat G, Raje R S, Bhutani S, Pal J K, Mithra A S V C R, Gaikwad K, Sharma T R, Singh N K. 2012. Molecular mapping of QTLs for plant type and earliness traits in pigeonpea (Cajanus cajan L. Millsp. ). BMC Genetics, 13(1): 84. DOI:10.1186/1471-2156-13-84
Li H, Durbin R. 2009. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics, 25(14): 1 754-1 760. DOI:10.1093/bioinformatics/btp324
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, 1000 Genome Project Data Processing Subgroup. 2009. The sequence alignment/map format and SAMtools. Bioinformatics, 25(16): 2 078-2 079. DOI:10.1093/bioinformatics/btp352
Li H, Hearne S, Bänziger M, Li Z, Wang J. 2010. Statistical properties of QTL linkage mapping in biparental genetic populations. Heredity, 105(3): 257-267. DOI:10.1038/hdy.2010.56
Liu J, Shikano T, Leinonen T, Cano J M, Li M H, Merilä J. 2014. Identification of major and minor QTL for ecologically important morphological traits in three-spined sticklebacks (Gasterosteus aculeatus). Genes, Genomes, Genetics, 4(4): 595-604. DOI:10.1534/g3.114.010389
Liu K Y, Xu H, Liu G, Guan P F, Zhou X Y, Peng H R, Yao Y Y, Ni Z F, Sun Q X, Du J K. 2018. QTL mapping of flag leaf-related traits in wheat (Triticum aestivum L. ). Theoretical and Applied Genetics, 131(4): 839-849. DOI:10.1007/s00122-017-3040-z
Liu M J, Huang L B, Yan X H. 2015. Isolation and characterization of the improved strain HW-4 by intraspecific hybridization in Pyropia yezoensis. Journal of Fishery Sciences of China, 22(1): 33-43. (in Chinese with English abstract) DOI:10.3724/SP.J.1118.2015.00213
Lowry D B, Taylor S H, Bonnette J, Aspinwall M J, Asmus A L, Keitt T H, Tobias C M, Juenger T E. 2015. QTLs for biomass and developmental traits in Switchgrass (Panicum virgatum). Bioenergy Research, 8(4): 1 856-1 867. DOI:10.1007/s12155-015-9629-7
Ma L L, Guan Z R, Zhang Z T, Zhang X X, Zhang Y L, Zou C Y, Peng H W, Pan G T, Lee M, Shen Y O, Lübberstedt T. 2018. Identification of quantitative trait loci for leaf-related traits in an IBM Syn10 DH maize population across three environments. Plant Breeding, 137(2): 127-138. DOI:10.1111/pbr.12566
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo M A. 2010. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Research, 20(9): 1 297-1 303. DOI:10.1101/gr.107524.110
Migita S. 1967. Cytological studies on Porphyra yezoensis Ueda. Nagasaki University Bulletin of the Faculty of Fisheries, 24: 55-64.
Miura A, Shin J A. 1989. Crossbreeding in cultivars of Porphyra yezoensis (Bangiales, Rhodophyta)-Preliminary report. The Korean Journal of Phycology, 4(2): 207-211.
Mo H D. 2003. Look back and reflect on genetic researches of variation for quantitative traits: a challenge for quantitative genetics in post-genome ear. Journal of Yangzhou University (Agricultural and Life Sciences Edition), 24(2): 24-31. (in Chinese with English abstract)
Moose S P, Mumm R H. 2008. Molecular plant breeding as the foundation for 21st century crop improvement. Plant Physiology, 147(3): 969-977. DOI:10.1104/pp.108.118232
Nakamura Y, Sasaki N, Kobayashi M, Ojima N, Yasuike M, Shigenobu Y, Satomi M, Fukuma Y, Shiwaku K, Tsujimoto A, Kobayashi T, Nakayama I, Ito F, Nakajima K, Sano M, Wada T, Kuhara S, Inouye K, Gojobori T, Ikeo K. 2013. The first symbiont-free genome sequence of marine red alga, Susabi-nori (Pyropia yezoensis). PLoS One, 8(3): e57122. DOI:10.1371/journal.pone.0057122
Niwa K, Yamamoto T, Furuita H, Abe T. 2011. Mutation breeding in the marine crop Porphyra yezoensis (Bangiales, Rhodophyta): cultivation experiment of the artificial red mutant isolated by heavy-ion beam mutagenesis. Aquaculture, 314(1-4): 182-187. DOI:10.1016/j.aquaculture.2011.02.007
Niwa K. 2010. Genetic analysis of artificial green and red mutants of Porphyra yezoensis Ueda (Bangiales, Rhodophyta). Aquaculture, 308(1-2): 6-12. DOI:10.1016/j.aquaculture.2010.08.007
Ohme M, Kunifuji Y, Miura A. 1986. Cross experiments of the color mutants in Porphyra yezoensis Ueda. Journal of Applied Phycology, 34: 101-106.
Ohme M, Miura A. 1988. Tetrad analysis in conchospore germlings of Porphyra yezoensis (Rhodophyta, Bangiales). Plant Science, 57(2): 135-140. DOI:10.1016/0168-9452(88)90079-9
Qing D J, Dai G X, Zhou W Y, Huang S S, Liang H F, Gao L J, Gao J, Huang J, Zhou M, Chen R T, Chen W W, Huang F K, Deng G F. 2019. Development of molecular marker and introgression of Bph3 into elite rice cultivars by marker-assisted selection. Breeding Science, 69(1): 40-46. DOI:10.1270/jsbbs.18080
Randhawa M S, Bains N S, Sohu V S, Chhuneja P, Trethowan R M, Bariana H S, Bansal U. 2019. Marker assisted transfer of stripe rust and stem rust resistance genes into four wheat cultivars. Agronomy, 9(9): 497. DOI:10.3390/agronomy9090497
Ribaut J M, De Vicente M C, Delannay X. 2010. Molecular breeding in developing countries: challenges and perspectives. Current Opinion in Plant Biology, 13(2): 213-218. DOI:10.1016/j.pbi.2009.12.011
Shan T F, Pang S J, Li J, Li X, Su L. 2015. Construction of a high-density genetic map and mapping of a sex-linked locus for the brown alga Undaria pinnatifida (Phaeophyceae) based on large scale marker development by specific length amplified fragment (SLAF) sequencing. BMC Genomics, 16(1): 902. DOI:10.1186/s12864-015-2184-y
Shen Y S, Yang Y, Xu E S, Ge X H, Xiang Y, Li Z Y. 2018. Novel and major QTL for branch angle detected by using DH population from an exotic introgression in rapeseed (Brassica napus L. ). . Theoretical and Applied Genetics, 131(1): 67-78. DOI:10.1007/s00122-017-2986-1
Sitonik C A, Suresh L M, Beyene Y, Olsen M S, Makumbi D, Oliver K, Das B, Bright J M, Mugo S, Crossa J, Tarekegne A, Prasanna B M, Gowda M. 2019. Genetic architecture of maize chlorotic mottle virus and maize lethal necrosis through GWAS, linkage analysis and genomic prediction in tropical maize germplasm. Theoretical and Applied Genetics, 132(8): 2 381-2 399. DOI:10.1007/s00122-019-03360-x
Sooriyapathirana S S, Khan A, Sebolt A M, Wang D C, Bushakra J M, Kui L W, Allan A C, Gardiner S E, Chagné D, Iezzoni A F. 2010. QTL analysis and candidate gene mapping for skin and flesh color in sweet cherry fruit (Prunus avium L. ). . Tree Genetics & Genomes, 6(6): 821-832. DOI:10.1007/s11295-010-0294-x
Sun X W, Liu D Y, Zhang X F, Li W B, Liu H, Hong W G, Jiang C B, Guan N, Ma C X, Zeng H P, Xu C H, Song J, Huang L, Wang C M, Shi J J, Wang R, Zheng X H, Lu C Y, Wang X W, Zheng H K. 2013. SLAF-seq: an efficient method of large-scale de novo SNP discovery and genotyping using high-throughput sequencing. PLoS One, 8(3): e58700. DOI:10.1371/journal.pone.0058700
Sutherland J E, Lindstrom S C, Nelson W A, Brodie J, Lynch M D J, Hwang M S, Choi H G, Miyata M, Kikuchi N, Oliveira M C, Farr T, Neefus C, Mols-Mortensen A, Milstein D, Müller K M. 2011. A new look at an ancient order: generic revision of the Bangiales (Rhodophyta). Journal of Phycology, 47(5): 1 131-1 151. DOI:10.1111/j.1529-8817.2011.01052.x
Tao Y F, Liu Q C, Wang H H, Zhang Y J, Huang X Y, Wang B B, Lai J S, Ye J R, Liu B S, Xu M L. 2013. Identification and fine-mapping of a QTL, qMrdd1, that confers recessive resistance to maize rough dwarf disease. BMC Plant Biology, 13(1): 145. DOI:10.1186/1471-2229-13-145
Van Ooijen J W, Kyazma B V. 2009. MapQTL 6: Software for the Mapping of Quantitative Trait Loci in Experimental Populations of Diploid Species. Kyazma BV, Wageningen, Netherlands.
Wang B B, Liu H, Liu Z P, Dong X M, Guo J J, Li W, Chen J, Gao C, Zhu Y B, Zheng X M, Chen Z L, Chen J, Song W B, Hauck A, Lai J S. 2018a. Identification of minor effect QTLs for plant architecture related traits using super high density genotyping and large recombinant inbred population in maize (Zea mays). BMC Plant Biology, 18(1): 17. DOI:10.1186/s12870-018-1233-5
Wang C N, Qiao A H, Fang X F, Sun L, Gao P, Davis A R, Liu S, Luan F S. 2019. Fine mapping of lycopene content and flesh color related gene and development of molecular marker-assisted selection for flesh color in watermelon (Citrullus lanatus). Frontiers in Plant Science, 10: 1 240. DOI:10.3389/fpls.2019.01240
Wang S J, Zheng Y Z, Ma L B, Xu P, Zhu J Y. 2000. Gamma-rays induction of mutation in conchocelis of Porphyra yezoensis. Chinese Journal of Oceanology and Limnology, 18(1): 47. DOI:10.1007/BF02842541
Wang S, Basten C, Zeng Z. 2007. Windows QTL Cartographer 2.5. Department of Statistics, North Carolina State University, Raleigh, NC.
Wang X L, Chen Z H, Li Q Y, Zhang J, Liu S, Duan D L. 2018b. High-density SNP-based QTL mapping and candidate gene screening for yield-related blade length and width in Saccharina japonica (Laminariales, Phaeophyta). Scientific Reports, 8(1): 13 591. DOI:10.1038/s41598-018-32015-y
Xu D H, Sun R F, Zhang Y G, Yuan Y X, Kang J G, Wu J, Zhang H, Song X F, Li X N, Song Y M, Wang X W. 2007. Mapping and analysis of QTL related to leaf color in Chinese cabbage (Brassica rapa L. ssp. pekinensis). Acta Horticulturae Sinica, 34(1): 99-104. (in Chinese with English abstract) DOI:10.3321/j.issn:0513-353X.2007.01.020
Xu K P, Yu X Z, Tang X H, Kong F N, Mao Y X. 2019. Organellar genome variation and genetic diversity of Chinese Pyropia yezoensis. Frontiers in Marine Science, 6: 756. DOI:10.3389/fmars.2019.00756
Xu Y, Huang L, Ji D H, Chen C S, Zheng H K, Xie C T. 2015. Construction of a dense genetic linkage map and mapping quantitative trait loci for economic traits of a doubled haploid population of Pyropia haitanensis (Bangiales, Rhodophyta). BMC Plant Biology, 15(1): 228. DOI:10.1186/s12870-015-0604-4
Yabu H. 1969. Observation on chromosomes in some species of Porphyra. Hokkaido University Bulletin of the Faculty of Fisheries, 19(4): 239-243.
Yan X H, Aruga Y. 2000. Genetic analysis of artificial pigmentation mutants in Porphyra yezoensis Ueda (Bangiales, Rhodophyta). Phycological Research, 48(3): 177-187. DOI:10.1046/j.1440-1835.2000.00203.x
Yan X H, Fujita Y, Aruga Y. 2000. Induction and characterization of pigmentation mutants in Porphyra yezoensis (Bangiales, Rhodophyta). Journal of Applied Phycology, 12(1): 69-81. DOI:10.1023/A:1008129119065
Yan X H, Fujita Y, Aruga Y. 2004. High monospore-producing mutants obtained by treatment with MNNG in Porphyra yezoensis Ueda (Bangiales, Rhodophyta). Hydrobiologia, 512(1-3): 133-140. DOI:10.1023/B:Hydr.0000020318.00749.2d
Yan X H, Liang Z Q, Song W L, Huang J, Ma P, Aruga Y. 2005. Induction and isolation of artificial pigmentation mutants in Porphyra haitanensis Chang et Zheng (Bangiales, Rhodophyta).. Journal of Fisheries of China, 29(2): 166-172. (in Chinese with English abstract) DOI:10.3321/j.issn:1000-0615.2005.02.005
Yang G P, Sun Y, Shi Y Y, Zhang L N, Guo S S, Li B J, Li X J, Li Z L, Cong Y Z, Zhao Y S, Wang W Q. 2009. Construction and characterization of a tentative amplified fragment length polymorphism-simple sequence repeat linkage map of Laminaria (Laminariales, Phaeophyta). Journal of Phycology, 45(4): 873-878. DOI:10.1111/j.1529-8817.2009.00720.x
Yu Y, Zhang X J, Li F H, Xiang J H. 2011. Strategy of whole genomic selection breeding and its application prospect in aquaculture. Journal of Fishery Sciences of China, 18(4): 936-943. (in Chinese with English abstract) DOI:10.3724/SP.J.1118.2011.00935
Zeng Z B. 1994. Precision mapping of quantitative trait loci. Genetics, 136(4): 1 457-1 468. DOI:10.1093/genetics/136.4.1457
Zhang W J, Lukaszewski A J, Kolmer J, Soria M A, Goyal S, Dubcovsky J. 2005. Molecular characterization of durum and common wheat recombinant lines carrying leaf rust resistance (Lr19) and yellow pigment (Y) genes from Lophopyrum ponticum. Theoretical and Applied Genetics, 111(3): 573-582. DOI:10.1007/s00122-005-2048-y
Zhu H Q, Zhang X G, Xue Q Z. 2004. Developing practicable system of chromosomal doubling for the construction of doubled haploid population in tobacco. Molecular Plant Breeding, 2(5): 643-648. (in Chinese with English abstract) DOI:10.3969/j.issn.1672-416X.2004.05.007