Construction of genetic linkage map and identification of a novel major locus for resistance to pine wood nematode in Japanese black pine (Pinus thunbergii)

Haplotype-resolved genome assembly and resequencing analysis provide insights into genome evolution and allelic imbalance in Pinus densiflora

Haplotype-level allelic characterization facilitates research on the functional, evolutionary and breeding-related features of extremely large and complex plant genomes. We report a 21.7-Gb chromosome-level haplotype-resolved assembly in Pinus densiflora. We found genome rearrangements involving translocations and inversions between chromosomes 1 and 3 of Pinus species and a proliferation of specific long terminal repeat (LTR) retrotransposons (LTR-RTs) in P. densiflora. Evolutionary analyses illustrated that tandem and LTR-RT-mediated duplications led to an increment of transcription factor (TF) genes in P. densiflora. The haplotype sequence comparison showed allelic imbalances, including presence-absence variations of genes (PAV genes) and their functional contributions to flowering and abiotic stress-related traits in P. densiflora. Allele-aware resequencing analysis revealed PAV gene diversity across P. densiflora accessions. Our study provides insights into key mechanisms underlying the evolution of genome structure, LTR-RTs and TFs within the Pinus lineage as well as allelic imbalances and diversity across P. densiflora.

Transcription Factor and Protein Regulatory Network of PmACRE1 in Pinus massoniana Response to Pine Wilt Nematode Infection

Pine wilt disease, caused by Bursaphelenchus xylophilus, is a highly destructive and contagious forest affliction. Often termed the "cancer" of pine trees, it severely impacts the growth of Masson pine (Pinus massoniana). Previous studies have demonstrated that ectopic expression of the PmACRE1 gene from P. massoniana in Arabidopsis thaliana notably enhances resistance to pine wilt nematode infection. To further elucidate the transcriptional regulation and protein interactions of the PmACRE1 in P. massoniana in response to pine wilt nematode infection, we cloned a 1984 bp promoter fragment of the PmACRE1 gene, a transient expression vector was constructed by fusing this promoter with the reporter GFP gene, which successfully activated the GFP expression. DNA pull-down assays identified PmMYB8 as a trans-acting factor regulating PmACRE1 gene expression. Subsequently, we found that the PmACRE1 protein interacts with several proteins, including the ATP synthase CF1 α subunit, ATP synthase CF1 β subunit, extracellular calcium-sensing receptor (PmCAS), caffeoyl-CoA 3-O-methyltransferase (PmCCoAOMT), glutathione peroxidase, NAD+-dependent glyceraldehyde-3-phosphate dehydrogenase, phosphoglycerate kinase 1, cinnamyl alcohol dehydrogenase, auxin response factor 16, and dehydrin 1 protein. Bimolecular fluorescence complementation (BiFC) assays confirmed the interactions between PmACRE1 and PmCCoAOMT, as well as PmCAS proteins in vitro. These findings provide preliminary insights into the regulatory role of PmACRE1 in P. massoniana's defense against pine wilt nematode infection.

Biomass-Based, Dual Enzyme-Responsive Nanopesticides: Eco-friendly and Efficient Control of Pine Wood Nematode Disease

Pine wood nematode (PWN) disease is a globally devastating forest disease caused by infestation with PWN, Bursaphelenchus xylophilus, which mainly occurs through the vector insect Japanese pine sawyer (JPS), Monochamus alternatus. PWN disease is notoriously difficult to manage effectively and is known as the "cancer of pine trees." In this study, dual enzyme-responsive nanopesticides (AVM@EC@Pectin) were prepared using nanocoating avermectin (AVM) after modification with natural polymers. The proposed treatment can respond to the cell wall-degrading enzymes secreted by PWNs and vector insects during pine tree infestation to intelligently release pesticides to cut off the transmission and infestation pathways and realize the integrated control of PWN disease. The LC50 value of AVM@EC@Pectin was 11.19 mg/L for PWN and 26.31 mg/L for JPS. The insecticidal activity of AVM@EC@Pectin was higher than that of the commercial emulsifiable concentrate (AVM-EC), and the photostability, adhesion, and target penetration were improved. The half-life (t1/2) of AVM@EC@Pectin was 133.7 min, which is approximately twice that of AVM-EC (68.2 min). Sprayed and injected applications showed that nanopesticides had superior bidirectional transportation, with five-times higher AVM contents detected in the roots relative to those of AVM-EC when sprayed at the top. The safety experiment showed that the proposed treatment had lower toxicity and higher safety for nontarget organisms in the application environment and human cells. This study presents a green, safe, and effective strategy for the integrated management of PWN disease.

Pinus thunbergii Parl. Somatic Plants' Resistance to Bursaphelenchus xylophilus Depends on Pathogen-Induced Differential Transcriptomic Responses

Pinus thunbergii Parl. is an economically and medicinally important plant, as well as a world-renowned horticultural species of the Pinus genus. Pine wilt disease is a dangerous condition that affects P. thunbergii. However, understanding of the genetics underlying resistance to this disease is poor. Our findings reveal that P. thunbergii's resistance mechanism is based on differential transcriptome responses generated by the early presence of the pathogen Bursaphelenchus xylophilus, also known as the pine wood nematode. A transcriptome analysis (RNA-seq) was performed to examine gene expression in shoot tissues from resistant and susceptible P. thunbergii trees. RNA samples were collected from the shoots of inoculated pines throughout the infection phases by the virulent Bursaphelenchus xylophilus AMA3 strain. The photosynthesis and plant-pathogen interaction pathways were significantly enriched in the first and third days after infection. Flavonoid biosynthesis was induced in response to late infestation (7 and 14 days post-infestation). Calmodulin, RBOH, HLC protein, RPS, PR1, and genes implicated in phytohormone crosstalk (e.g., SGT1, MYC2, PP2C, and ERF1) showed significant alterations between resistant and susceptible trees. Furthermore, salicylic acid was found to aid pine wood nematodes tolerate adverse conditions and boost reproduction, which may be significant for pine wood nematode colonization within pines. These findings provide new insights into how host defenses overcame pine wood nematode infection in the early stage, which could potentially contribute to the development of novel strategies for the control of pine wilt disease.

Resistant and Susceptible Pinus thunbergii ParL. Show Highly Divergent Patterns of Differentially Expressed Genes during the Process of Infection by Bursaphelenchus xylophilus

Pine wilt disease (PWD) is a devastating disease that threatens pine forests worldwide, and breeding resistant pines is an important management strategy used to reduce its impact. A batch of resistant seeds of P. thunbergii was introduced from Japan. Based on the resistant materials, we obtained somatic plants through somatic embryogenesis. In this study, we performed transcriptome analysis to further understand the defense response of resistant somatic plants of P. thunbergii to PWD. The results showed that, after pine wood nematode (PWN) infection, resistant P. thunbergii stimulated more differential expression genes (DEGs) and involved more regulatory pathways than did susceptible P. thunbergii. For the first time, the alpha-linolenic acid metabolism and linoleic acid metabolism were intensively observed in pines resisting PWN infection. The related genes disease resistance protein RPS2 (SUMM2) and pathogenesis-related genes (PR1), as well as reactive oxygen species (ROS)-related genes were significantly up-expressed in order to contribute to protection against PWN inoculation in P. thunbergii. In addition, the diterpenoid biosynthesis pathway was significantly enriched only in resistant P. thunbergii. These findings provided valuable genetic information for future breeding of resistant conifers, and could contribute to the development of new diagnostic tools for early screening of resistant pine seedlings based on specific PWN-tolerance-related markers.

Transcriptomic analysis reveals differentially expressed genes associated with pine wood nematode resistance in resistant Pinus thunbergii

Pine wilt disease (caused by the nematode Bursaphelenchus xylophilus) is extremely harmful to pine forests in East Asia. As a low-resistance pine species, Pinus thunbergii is more vulnerable to pine wood nematode (PWN) than Pinus densiflora and Pinus massoniana. Field inoculation experiments were conducted on PWN-resistant and -susceptible P. thunbergii, and the difference in transcription profiles 24 h after inoculation was analyzed. We identified 2603 differentially expressed genes (DEGs) in PWN-susceptible P. thunbergii, while 2559 DEGs were identified in PWN-resistant P. thunbergii. Before inoculation, DEGs between PWN-resistant and PWN-susceptible P. thunbergii were enriched in the REDOX (Oxidation-Reduction) activity pathway (152 DEGs), followed by the oxidoreductase activity pathway (106 DEGs). After inoculation with PWN, however, the opposite was observed; DEGs were enriched in the oxidoreductase activity pathway (119 DEGs), followed by the REDOX activity pathway (84 DEGs). Before inoculation, according to the metabolic pathway analysis results, we found more genes upregulated in phenylpropanoid metabolic pathways and enriched in lignin synthesis pathways; cinnamoyl-CoA reductase-coding genes related to lignin synthesis were upregulated in PWN-resistant P. thunbergii and downregulated in PWN-susceptible P. thunbergii, and the lignin content was always higher in resistant than in susceptible P. thunbergii. These results reveal distinctive strategies of resistant and susceptible P. thunbergii in dealing with PWN infections.

Promoting the application of Pinus thunbergii Parl. to enhance the growth and survival rates of post-germination somatic plantlets

OBJECTIVE

There is a growing need for nematode resistant Pinaceae species plantlets to cope with the global scale degradation of coniferous forests, due to the prevalence of pine wilt disease. One of the bottlenecks that limits the commercialization of Pinaceae species plantlets is regeneration following their transfer from controlled sterile environments to the field while maintaining high survival rates.

METHODS

The growth factors of somatic plantlets (SPs), such as sucrose, media, culture substrate, brassinolide and spectrum were investigated to promote the application of somatic nematode-resistant P. thunbergii plants in afforestation.

RESULTS

The 1/2 WPM liquid medium, culture substrate (perlite and vermiculite =1:1), and carbohydrate (20 g/L sucrose) were effective in stimulating the growth of rooted SPs. While for unrooted SPs, 1 ug/L of brassinolide enhanced plantlet growth and rooting. And blue light (B) significantly promoted the longitudinal growth of shoots, while red light (R) was beneficial for root growth during the laboratory domestication stage. High quality SPs were obtained at a R/B ratio of 8:2. Following this acclimatization protocol, the P. thunbergii SPs could be directly transplanted to the field with a higher survival rate (85.20 %) in a forcing house.

CONCLUSION

this acclimatization protocol extremely improved the survival rate of P. thunbergii SPs. Moreover, this work will contribute to enhancing the possibilities for somatic plant afforestation with Pinus species.

Soil Fungal Community Structure and Function Shift during a Disease-Driven Forest Succession

Forest succession is important for sustainable forest management in terrestrial ecosystems. However, knowledge about the response of soil microbes to forest disease-driven succession is limited. In this study, we investigated the soil fungal biomass, soil enzyme activity, and fungal community structure and function in forests suffering succession processes produced by pine wilt disease from conifer to broadleaved forests using Illumina Miseq sequencing coupled with FUNGuild analysis. The results showed that the broadleaved forest had the highest fungal biomass and soil enzyme activities in C, N, and S cycles, whereas the conifer forest had the highest enzyme activity in the P cycle. Along the succession, the fungal diversity and richness significantly increased (P < 0.05). The fungal communities were dominated by Ascomycota (42.0%), Basidiomycota (38.0%), and Mortierellomycota (9.5%), among which the abundance of Ascomycota significantly increased (P < 0.05), whereas that of Basidiomycota and Mortierellomycota decreased (P < 0.05). The abundance of species Mortierella humilis, Lactarius salmonicolor, and Russula sanguinea decreased, whereas that of Mortierella minutissima increased (P < 0.05). The forests in different succession stages formed distinct fungal communities and functional structures (P < 0.05). Functionally, the saprotrophs, symbiotrophs, and pathotrophs were the dominant groups in the conifer, mixed, and broadleaved forests, respectively. Soil pH and soil organic carbon were the key factors influencing the fungal community and functional structures during the succession. These findings provide useful information for better understanding the plant-microbe interaction during forest succession caused by forest disease. IMPORTANCE The studies on soil fungal communities in disease-driven forest succession are rare. This study showed that during the disease-driven forest succession, the soil enzyme activity, soil fungal diversity, and biomass increased along succession. The disease-driven forest succession changed the soil fungal community structure and function, in which the symbiotrophs were the most dominant group along the succession. These findings provide useful information for better understanding the plant-microbe interaction during forest succession caused by forest disease.

Comprehensive collection of genes and comparative analysis of full-length transcriptome sequences from Japanese larch (Larix kaempferi) and Kuril larch (Larix gmelinii var. japonica)

BACKGROUND

Japanese larch (Larix kaempferi) is an economically important deciduous conifer species that grows in cool-temperate forests and is endemic to Japan. Kuril larch (L. gmelinii var. japonica) is a variety of Dahurian larch that is naturally distributed in the Kuril Islands and Sakhalin. The hybrid larch (L. gmelinii var. japonica × L. kaempferi) exhibits heterosis, which manifests as rapid juvenile growth and high resistance to vole grazing. Since these superior characteristics have been valued by forestry managers, the hybrid larch is one of the most important plantation species in Hokkaido. To accelerate molecular breeding in these species, we collected and compared full-length cDNA isoforms (Iso-Seq) and RNA-Seq short-read, and merged them to construct candidate gene as reference for both Larix species. To validate the results, candidate protein-coding genes (ORFs) related to some flowering signal-related genes ​were screened from the reference sequences, and the phylogenetic relationship with closely related species was elucidated.

RESULTS

Using the isoform sequencing of PacBio RS ll and the de novo assembly of RNA-Seq short-read sequences, we identified 50,690 and 38,684 ORFs in Japanese larch and Kuril larch, respectively. BUSCO completeness values were 90.5% and 92.1% in the Japanese and Kuril larches, respectively. After comparing the collected ORFs from the two larch species, a total of 19,813 clusters, comprising 22,571 Japanese larch ORFs and 22,667 Kuril larch ORFs, were contained in the intersection of the Venn diagram. In addition, we screened several ORFs related to flowering signals (SUPPRESSER OF OVEREXPRESSION OF CO1: SOC1, LEAFY: LFY, FLOWERING Locus T: FT, CONSTANCE: CO) from both reference sequences, and very similar found in other species.

CONCLUSIONS

The collected ORFs will be useful as reference sequences for molecular breeding of Japanese and Kuril larches, and also for clarifying the evolution of the conifer genome and investigating functional genomics.

SNP Detection in Pinus pinaster Transcriptome and Association with Resistance to Pinewood Nematode

Pinewood nematode (PWN, Bursaphelenchus xylophilus) is the causal agent of pine wilt disease (PWD), which severely affects Pinus pinaster stands in southwestern Europe. Despite the high susceptibility of P. pinaster, individuals of selected half-sib families have shown genetic variability in survival after PWN inoculation, indicating that breeding for resistance can be a valuable strategy to control PWD. In this work, RNA-seq data from susceptible and resistant plants inoculated with PWN were used for SNP discovery and analysis. A total of 186,506 SNPs were identified, of which 31 were highly differentiated between resistant and susceptible plants, including SNPs in genes involved in cell wall lignification, a process previously linked to PWN resistance. Fifteen of these SNPs were selected for validation through Sanger sequencing and 14 were validated. To evaluate SNP-phenotype associations, 40 half-sib plants were genotyped for six validated SNPs. Associations with phenotype after PWN inoculation were found for two SNPs in two different genes (MEE12 and PCMP-E91), as well as two haplotypes of HIPP41, although significance was not maintained following Bonferroni correction. SNPs here detected may be useful for the development of molecular markers for PWD resistance and should be further investigated in future association studies.

Quantitative Trait Loci Analysis Based on High-Density Mapping of Single-Nucleotide Polymorphisms by Genotyping-by-Sequencing Against Pine Wilt Disease in Japanese Black Pine (Pinus thunbergii)

Identifying genes/loci for resistance to pine wilt disease (PWD) caused by the pine wood nematode (PWN) is beneficial for improving resistance breeding in Pinus thunbergii, but to date, genetic information using molecular markers has been limited. Here, we constructed a high-density linkage map using genotyping-by-sequencing (GBS) and conducted quantitative trait loci (QTL) analysis for PWD resistance for the self-pollinated progeny of "Namikata 73," which is the most resistant variety among resistant varieties of P. thunbergii, following inoculation tests with PWN. An S1 mapping population consisting of the 116 progenies derived from self-pollination of the resistant variety, "Namikata 73" (resistance rank 5 to PWN), was inoculated with PWN isolate Ka-4 and evaluated for disease symptoms. To construct a high-density linkage map, we used single-nucleotide polymorphisms (SNPs) identified by GBS based on next-generation sequencing technology and some anchor DNA markers, expressed sequence tag (EST)-derived SNP markers and EST-derived simple sequence repeat (SSR) markers, and genomic SSR markers. The linkage map had 13 linkage groups (LGs) consisting of 2,365 markers including 2,243 GBS-SNP markers over a total map distance of 1968.4 centimorgans (cM). Results from QTL analysis using phenotype data and the linkage map indicated that PWD resistance is controlled by a single locus located on LG-3, as identified in a previous study. This locus showed overdominant genetic action in the present study. With the confirmation of PWD1 in two different mapping populations (present study and a previous study), the locus associated with this region is thought to be a good target for marker-assisted selection in P. thunbergii breeding programs in order to obtain high levels of resistance to PWD caused by PWN.

QTL mapping for growth-related traits by constructing the first genetic linkage map in Simao pine

BACKGROUND

Simao pine is one of the primary economic tree species for resin and timber production in southwest China. The exploitation and utilization of Simao pine are constrained by the relatively lacking of genetic information. Construction a fine genetic linkage map and detecting quantitative trait locis (QTLs) for growth-related traits is a prerequisite section of Simao Pine's molecular breeding program.

RESULTS

In our study, a high-resolution Simao pine genetic map employed specific locus amplified fragment sequencing (SLAF-seq) technology and based on an F₁ pseudo-testcross population has been constructed. There were 11,544 SNPs assigned to 12 linkage groups (LGs), and the total length of the map was 2,062.85 cM with a mean distance of 0.37 cM between markers. According to the phenotypic variation analysis for three consecutive years, a total of seventeen QTLs for four traits were detected. Among 17 QTLs, there were six for plant height (Dh.16.1, Dh16.2, Dh17.1, Dh18.1-3), five for basal diameter (Dbd.17.1-5), four for needle length (Dnl17.1-3, Dnl18.1) and two for needle diameter (Dnd17.1 and Dnd18.1) respectively. These QTLs individually explained phenotypic variance from 11.0-16.3%, and the logarithm of odds (LOD) value ranged from 2.52 to 3.87.

CONCLUSIONS

In our study, a fine genetic map of Simao pine applied the technology of SLAF-seq has been constructed for the first time. Based on the map, a total of 17 QTLs for four growth-related traits were identified. It provides helpful information for genomic studies and marker-assisted selection (MAS) in Simao pine.

Optimization of factors affecting the rooting of pine wilt disease resistant Masson pine (Pinus massoniana) stem cuttings

Pine wilt disease (PWD) is a devastating disease affecting trees belonging to the genus Pinus. To control the spread of PWD in the Masson pine forest in China, PWD resistant Masson pine clones have been selected by the Anhui Academy of Forestry. However, because Masson pine is a difficult-to-root species, producing seedlings is challenging, especially from trees older than 5 years of age, which impedes the application of PWD resistant clones. In this study, we investigated the factors affecting rooting of PWD resistant clones and established a cheap, reliable, and simple method that promotes rooting. We tested the effects of three management methods, four substrates, two cutting materials, two cutting treatments, and three collection times on the rooting of cuttings obtained from 9-year-old PWD resistant clones. Rooting was observed only in stem cuttings treated with the full-light automatic spray management method. Additionally, stem cuttings showed a significantly higher rooting rate and root quality than needles cuttings. Compared with other substrates, stem cuttings planted in perlite produced the longest adventitious root and the highest total root length and lateral root number. Moreover, stem cuttings of PWD resistant clones collected in May showed a significantly higher rooting rate and root quality than those collected in June and July. Moreover, stem cuttings prepared with a horizontal cut while retaining the needles showed significantly higher rooting rate and root quality than those prepared with a diagonal cut while partly removing the needles. This study promotes the reproduction of seedlings of PWD-resistant Masson pine clones which helps control the spread of PWD, meanwhile, provides a technical reference for the propagation of mature pine trees via cuttings.

SNP Genotyping with Target Amplicon Sequencing Using a Multiplexed Primer Panel and Its Application to Genomic Prediction in Japanese Cedar, Cryptomeria japonica (L.f.) D.Don

Along with progress in sequencing technology and accumulating knowledge of genome and gene sequences, molecular breeding techniques have been developed for predicting the genetic potential of individual genotypes and for selecting superior individuals. For Japanese cedar (Cryptomeria japonica (L.f.) D.Don), which is the most common coniferous species in Japanese forestry, we constructed a custom primer panel for target amplicon sequencing in order to simultaneously determine 3034 informative single nucleotide polymorphisms (SNPs). We performed primary evaluation of the custom primer panel with actual sequencing and in silico PCR. Genotyped SNPs had a distribution over almost the entire region of the C. japonica linkage map and verified the high reproducibility of genotype calls compared to SNPs obtained by genotyping arrays. Genotyping was performed for 576 individuals of the F1 population, and genomic prediction models were constructed for growth and wood property-related traits using the genotypes. Amplicon sequencing with the custom primer panel enables efficient obtaining genotype data in order to perform genomic prediction, manage clones, and advance forest tree breeding.