An improved assembly of the loblolly pine mega-genome using long-read single-molecule sequencing

Experimental Induction of Extreme Indented Growth Rings (Hazel Wood) in Pinus halepensis Miller by Wide and Long Parallel Bark and Vascular Cambium Woundings

Indented growth rings were found long ago to be experimentally induced in Pinus halepensis Miller by thin parallel axial scratching of the bark up to the vascular cambium with a sharp blade. Here, we show that when the bark and vascular cambium of P. halepensis are wounded by wide and long parallel axial wounds ("windows") rather than by thin scratches, the induced indented growth rings become dramatically more indented. All ten trees that were wounded by long parallel "windows" responded with very strong growth (especially in the first two years) that resulted in the formation of very conspicuous, extremely indented growth rings in the wood formed in between the long and wide woundings. This is true for both the trunks that were wounded all around their circumference and those that were wounded only in part of their circumference. We also suggest further lines of research.

Genome-wide identification and characterization of the AP2/ERF gene family in loblolly pine (Pinus taeda L.)

The loblolly pine (Pinus taeda L.) is one of the most profitable forest species worldwide owing to its quick growth, high wood yields, and strong adaptability. The AP2/ERF gene family plays a widespread role in the physiological processes of plant defense responses and the biosynthesis of metabolites. Nevertheless, there are no reports on this gene family in loblolly pine (P. taeda). In this study, a total of 303 members of the AP2/ERF gene family were identified. Through multiple sequence alignment and phylogenetic analysis, they were classified into four subfamilies, including AP2 (34), RAV (17), ERF (251), and Soloist (1). An analysis of the conservation domains, conserved motifs, and gene structure revealed that every PtAP2/ERF transcription factor (TF) had at least one AP2 domain. While evolutionary conservation was displayed within the same subfamilies, the distribution of conserved domains, conserved motifs, and gene architectures varied between subfamilies. Cis-element analysis revealed abundant light-responsive elements, phytohormone-responsive elements, and stress-responsive elements in the promoter of the PtAP2/ERF genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of potential target genes showed that the AP2/ERF gene family might play a critical role in plant growth and development, the response to environmental stresses, and metabolite biosynthesis. Utilizing quantitative real-time PCR (qRT-PCR), we examined the expression patterns of 10 randomly selected genes from Group IX after 6 h of treatments with mechanical injury, ethephon (Eth), and methyl jasmonate (MeJA). The AP2/ERF gene family in the loblolly pine was systematically analyzed for the first time in this study, offering a theoretical basis for exploring the functions and applications of AP2/ERF genes.

A genome sequence for the threatened whitebark pine

Whitebark pine (WBP, Pinus albicaulis) is a white pine of subalpine regions in the Western contiguous United States and Canada. WBP has become critically threatened throughout a significant part of its natural range due to mortality from the introduced fungal pathogen white pine blister rust (WPBR, Cronartium ribicola) and additional threats from mountain pine beetle (Dendroctonus ponderosae), wildfire, and maladaptation due to changing climate. Vast acreages of WBP have suffered nearly complete mortality. Genomic technologies can contribute to a faster, more cost-effective approach to the traditional practices of identifying disease-resistant, climate-adapted seed sources for restoration. With deep-coverage Illumina short reads of haploid megagametophyte tissue and Oxford Nanopore long reads of diploid needle tissue, followed by a hybrid, multistep assembly approach, we produced a final assembly containing 27.6 Gb of sequence in 92,740 contigs (N50 537,007 bp) and 34,716 scaffolds (N50 2.0 Gb). Approximately 87.2% (24.0 Gb) of total sequence was placed on the 12 WBP chromosomes. Annotation yielded 25,362 protein-coding genes, and over 77% of the genome was characterized as repeats. WBP has demonstrated the greatest variation in resistance to WPBR among the North American white pines. Candidate genes for quantitative resistance include disease resistance genes known as nucleotide-binding leucine-rich repeat receptors (NLRs). A combination of protein domain alignments and direct genome scanning was employed to fully describe the 3 subclasses of NLRs. Our high-quality reference sequence and annotation provide a marked improvement in NLR identification compared to previous assessments that leveraged de novo-assembled transcriptomes.

Assembly and annotation of the black spruce genome provide insights on spruce phylogeny and evolution of stress response

Black spruce (Picea mariana [Mill.] B.S.P.) is a dominant conifer species in the North American boreal forest that plays important ecological and economic roles. Here, we present the first genome assembly of P. mariana with a reconstructed genome size of 18.3 Gbp and NG50 scaffold length of 36.0 kbp. A total of 66,332 protein-coding sequences were predicted in silico and annotated based on sequence homology. We analyzed the evolutionary relationships between P. mariana and 5 other spruces for which complete nuclear and organelle genome sequences were available. The phylogenetic tree estimated from mitochondrial genome sequences agrees with biogeography; specifically, P. mariana was strongly supported as a sister lineage to P. glauca and 3 other taxa found in western North America, followed by the European Picea abies. We obtained mixed topologies with weaker statistical support in phylogenetic trees estimated from nuclear and chloroplast genome sequences, indicative of ancient reticulate evolution affecting these 2 genomes. Clustering of protein-coding sequences from the 6 Picea taxa and 2 Pinus species resulted in 34,776 orthogroups, 560 of which appeared to be specific to P. mariana. Analysis of these specific orthogroups and dN/dS analysis of positive selection signatures for 497 single-copy orthogroups identified gene functions mostly related to plant development and stress response. The P. mariana genome assembly and annotation provides a valuable resource for forest genetics research and applications in this broadly distributed species, especially in relation to climate adaptation.

Identification and high-throughput genotyping of single nucleotide polymorphism markers in a non-model conifer (Abies nordmanniana (Steven) Spach)

Single nucleotide polymorphism (SNP) markers are powerful tools for investigating population structures, linkage analysis, and genome-wide association studies, as well as for breeding and population management. The availability of SNP markers has been limited to the most commercially important timber species, primarily due to the cost of genome sequencing required for SNP discovery. In this study, a combination of reference-based and reference-free approaches were used to identify SNPs in Nordmann fir (Abies nordmanniana), a species previously lacking genomic sequence information. Using a combination of a genome assembly of the closely related Silver fir (Abies alba) species and a de novo assembly of low-copy regions of the Nordmann fir genome, we identified a high density of reliable SNPs. Reference-based approaches identified two million SNPs in common between the Silver fir genome and low-copy regions of Nordmann fir. A combination of one reference-free and two reference-based approaches identified 250 shared SNPs. A subset of 200 SNPs were used to genotype 342 individuals and thereby tested and validated in the context of identity analysis and/or clone identification. The tested SNPs successfully identified all ramets per clone and five mislabeled individuals via identity and genomic relatedness analysis. The identified SNPs will be used in ad hoc breeding of Nordmann fir in Denmark.

Identifying genetic variation associated with environmental gradients and drought-tolerance phenotypes in ponderosa pine

As climate changes, understanding the genetic basis of local adaptation in plants becomes an ever more pressing issue. Combining genotype-environment association (GEA) with genotype-phenotype association (GPA) analysis has an exciting potential to uncover the genetic basis of environmental responses. We use these approaches to identify genetic variants linked to local adaptation to drought in Pinus ponderosa. Over 4 million Single Nucleotide Polymorphisms (SNPs) were identified using 223 individuals from across the Sierra Nevada of California. 927,740 (22.3%) SNPs were retained after filtering for proximity to genes and used in our association analyses. We found 1374 associated with five major climate variables, with the largest number (1151) associated with April 1st snowpack. We also conducted a greenhouse study with various drought-tolerance traits measured in first-year seedlings of a subset of the genotyped trees grown in the greenhouse. 796 SNPs were associated with control-condition trait values, while 1149 were associated with responsiveness of these traits to drought. While no individual SNPs were associated with both the environmental variables and the measured traits, several annotated genes were associated with both, particularly those involved in cell wall formation, biotic and abiotic stress responses, and ubiquitination. However, the functions of many of the associated genes have not yet been determined due to the lack of gene annotation information for conifers. Future studies are needed to assess the developmental roles and ecological significance of these unknown genes.

The genome of the Wollemi pine, a critically endangered "living fossil" unchanged since the Cretaceous, reveals extensive ancient transposon activity

We present the genome of the living fossil, Wollemia nobilis, a southern hemisphere conifer morphologically unchanged since the Cretaceous. Presumed extinct until rediscovery in 1994, the Wollemi pine is critically endangered with less than 60 wild adults threatened by intensifying bushfires in the Blue Mountains of Australia. The 12 Gb genome is among the most contiguous large plant genomes assembled, with extremely low heterozygosity and unusual abundance of DNA transposons. Reduced representation and genome re-sequencing of individuals confirms a relictual population since the last major glacial/drying period in Australia, 120 ky BP. Small RNA and methylome sequencing reveal conservation of ancient silencing mechanisms despite the presence of thousands of active and abundant transposons, including some transferred horizontally to conifers from arthropods in the Jurassic. A retrotransposon burst 8-6 my BP coincided with population decline, possibly as an adaptation enhancing epigenetic diversity. Wollemia, like other conifers, is susceptible to Phytophthora, and a suite of defense genes, similar to those in loblolly pine, are targeted for silencing by sRNAs in leaves. The genome provides insight into the earliest seed plants, while enabling conservation efforts.

Strong Purifying Selection in Haploid Tissue-Specific Genes of Scots Pine Supports the Masking Theory

The masking theory states that genes expressed in a haploid stage will be under more efficient selection. In contrast, selection will be less efficient in genes expressed in a diploid stage, where the fitness effects of recessive deleterious or beneficial mutations can be hidden from selection in heterozygous form. This difference can influence several evolutionary processes such as the maintenance of genetic variation, adaptation rate, and genetic load. Masking theory expectations have been confirmed in single-cell haploid and diploid organisms. However, in multicellular organisms, such as plants, the effects of haploid selection are not clear-cut. In plants, the great majority of studies indicating haploid selection have been carried out using male haploid tissues in angiosperms. Hence, evidence in these systems is confounded with the effects of sexual selection and intraspecific competition. Evidence from other plant groups is scarce, and results show no support for the masking theory. Here, we have used a gymnosperm Scots pine megagametophyte, a maternally derived seed haploid tissue, and four diploid tissues to test the strength of purifying selection on a set of genes with tissue-specific expression. By using targeted resequencing data of those genes, we obtained estimates of genetic diversity, the site frequency spectrum of 0-fold and 4-fold sites, and inferred the distribution of fitness effects of new mutations in haploid and diploid tissue-specific genes. Our results show that purifying selection is stronger for tissue-specific genes expressed in the haploid megagametophyte tissue and that this signal of strong selection is not an artifact driven by high expression levels.

A long-read and short-read transcriptomics approach provides the first high-quality reference transcriptome and genome annotation for Pseudotsuga menziesii (Douglas-fir)

Douglas-fir (Pseudotsuga menziesii) is native to western North America. It grows in a wide range of environmental conditions and is an important timber tree. Although there are several studies on the gene expression responses of Douglas-fir to abiotic cues, the absence of high-quality transcriptome and genome data is a barrier to further investigation. Like for most conifers, the available transcriptome and genome reference dataset for Douglas-fir remains fragmented and requires refinement. We aimed to generate a highly accurate, and complete reference transcriptome and genome annotation. We deep-sequenced the transcriptome of Douglas-fir needles from seedlings that were grown under nonstress control conditions or a combination of heat and drought stress conditions using long-read (LR) and short-read (SR) sequencing platforms. We used 2 computational approaches, namely de novo and genome-guided LR transcriptome assembly. Using the LR de novo assembly, we identified 1.3X more high-quality transcripts, 1.85X more "complete" genes, and 2.7X more functionally annotated genes compared to the genome-guided assembly approach. We predicted 666 long noncoding RNAs and 12,778 unique protein-coding transcripts including 2,016 putative transcription factors. We leveraged the LR de novo assembled transcriptome with paired-end SR and a published single-end SR transcriptome to generate an improved genome annotation. This was conducted with BRAKER2 and refined based on functional annotation, repetitive content, and transcriptome alignment. This high-quality genome annotation has 51,419 unique gene models derived from 322,631 initial predictions. Overall, our informatics approach provides a new reference Douglas-fir transcriptome assembly and genome annotation with considerably improved completeness and functional annotation.

Targeted and untargeted metabolomics reveals deep analysis of drought stress responses in needles and roots of Pinus taeda seedlings

Drought stress is one of major environmental stresses affecting plant growth and yield. Although Pinus taeda trees are planted in rainy southern China, local drought sometime occurs and can last several months, further affecting their growth and resin production. In this study, P. taeda seedlings were treated with long-term drought (42 d), and then targeted and untargeted metabolomics analysis were carried out to evaluate drought tolerance of P. taeda. Targeted metabolomics analysis showed that levels of some sugars, phytohormones, and amino acids significantly increased in the roots and needles of water-stressed (WS) P. taeda seedlings, compared with well-watered (WW) pine seedlings. These metabolites included sucrose in pine roots, the phytohormones abscisic acid and sacylic acid in pine needles, the phytohormone gibberellin (GA4) and the two amino acids, glycine and asparagine, in WS pine roots. Compared with WW pine seedlings, the neurotransmitter acetylcholine significantly increased in needles of WS pine seedlings, but significantly reduced in their roots. The neurotransmitters L-glutamine and hydroxytyramine significantly increased in roots and needles of WS pine seedlings, respectively, compared with WW pine seedlings, but the neurotransmitter noradrenaline significantly reduced in needles of WS pine seedlings. Levels of some unsaturated fatty acids significantly reduced in roots or needles of WS pine seedlings, compared with WW pine seedlings, such as linoleic acid, oleic acid, myristelaidic acid, myristoleic acid in WS pine roots, and palmitelaidic acid, erucic acid, and alpha-linolenic acid in WS pine needles. However, three saturated fatty acids significantly increased in WS pine seedlings, i.e., dodecanoic acid in WS pine needles, tricosanoic acid and heptadecanoic acid in WS pine roots. Untargeted metabolomics analysis showed that levels of some metabolites increased in WS pine seedlings, especially sugars, long-chain lipids, flavonoids, and terpenoids. A few of specific metabolites increased greatly, such as androsin, piceatanol, and panaxatriol in roots and needles of WS pine seedlings. Comparing with WW pine seedlings, it was found that the most enriched pathways in WS pine needles included flavone and flavonol biosynthesis, ABC transporters, diterpenoid biosynthesis, plant hormone signal transduction, and flavonoid biosynthesis; in WS pine roots, the most enriched pathways included tryptophan metabolism, caffeine metabolism, sesquiterpenoid and triterpenoid biosynthesis, plant hormone signal transduction, biosynthesis of phenylalanine, tyrosine, and tryptophan. Under long-term drought stress, P. taeda seedlings showed their own metabolomics characteristics, and some new metabolites and biosynthesis pathways were found, providing a guideline for breeding drought-tolerant cultivars of P. taeda.

Genetic diversity and structure of the 4th cycle breeding population of Chinese fir (Cunninghamia lanceolata (lamb.) hook)

Studying population genetic structure and diversity is crucial for the marker-assisted selection and breeding of coniferous tree species. In this study, using RAD-seq technology, we developed 343,644 high-quality single nucleotide polymorphism (SNP) markers to resolve the genetic diversity and population genetic structure of 233 Chinese fir selected individuals from the 4^th cycle breeding program, representing different breeding generations and provenances. The genetic diversity of the 4^th cycle breeding population was high with nucleotide diversity (P_i ) of 0.003, and H_o and H_e of 0.215 and 0.233, respectively, indicating that the breeding population has a broad genetic base. The genetic differentiation level between the different breeding generations and different provenances was low (F_st < 0.05), with population structure analysis results dividing the 233 individuals into four subgroups. Each subgroup has a mixed branch with interpenetration and weak population structure, which might be related to breeding rather than provenance, with aggregation from the same source only being in the local branches. Our results provide a reference for further research on the marker-assisted selective breeding of Chinese fir and other coniferous trees.

A draft Diabrotica virgifera virgifera genome: insights into control and host plant adaption by a major maize pest insect

BACKGROUND

Adaptations by arthropod pests to host plant defenses of crops determine their impacts on agricultural production. The larval host range of western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is restricted to maize and a few grasses. Resistance of D. v. virgifera to crop rotation practices and multiple insecticides contributes to its status as the most damaging pest of cultivated maize in North America and Europe. The extent to which adaptations by this pest contributes to host plant specialization remains unknown.

RESULTS

A 2.42 Gb draft D. v. virgifera genome, Dvir_v2.0, was assembled from short shotgun reads and scaffolded using long-insert mate-pair, transcriptome and linked read data. K-mer analysis predicted a repeat content of ≥ 61.5%. Ortholog assignments for Dvir_2.0 RefSeq models predict a greater number of species-specific gene duplications, including expansions in ATP binding cassette transporter and chemosensory gene families, than in other Coleoptera. A majority of annotated D. v. virgifera cytochrome P450s belong to CYP4, 6, and 9 clades. A total of 5,404 transcripts were differentially-expressed between D. v. virgifera larvae fed maize roots compared to alternative host (Miscanthus), a marginal host (Panicum virgatum), a poor host (Sorghum bicolor) and starvation treatments; Among differentially-expressed transcripts, 1,908 were shared across treatments and the least number were between Miscanthus compared to maize. Differentially-expressed transcripts were enriched for putative spliceosome, proteosome, and intracellular transport functions. General stress pathway functions were unique and enriched among up-regulated transcripts in marginal host, poor host, and starvation responses compared to responses on primary (maize) and alternate hosts.

CONCLUSIONS

Manual annotation of D. v. virgifera Dvir_2.0 RefSeq models predicted expansion of paralogs with gene families putatively involved in insecticide resistance and chemosensory perception. Our study also suggests that adaptations of D. v. virgifera larvae to feeding on an alternate host plant invoke fewer transcriptional changes compared to marginal or poor hosts. The shared up-regulation of stress response pathways between marginal host and poor host, and starvation treatments may reflect nutrient deprivation. This study provides insight into transcriptomic responses of larval feeding on different host plants and resources for genomic research on this economically significant pest of maize.

Proteomics research in forest trees: A 2012-2022 update

This review is a compilation of proteomic studies on forest tree species published in the last decade (2012-2022), mostly focused on the most investigated species, including Eucalyptus, Pinus, and Quercus. Improvements in equipment, platforms, and methods in addition to the increasing availability of genomic data have favored the biological knowledge of these species at the molecular, organismal, and community levels. Integration of proteomics with physiological, biochemical and other large-scale omics in the direction of the Systems Biology, will provide a comprehensive understanding of different biological processes, from growth and development to responses to biotic and abiotic stresses. As main issue we envisage that proteomics in long-living plants will thrive light on the plant responses and resilience to global climate change, contributing to climate mitigation strategies and molecular breeding programs. Proteomics not only will provide a molecular knowledge of the mechanisms of resilience to either biotic or abiotic stresses, but also will allow the identification on key gene products and its interaction. Proteomics research has also a translational character being applied to the characterization of the variability and biodiversity, as well as to wood and non-wood derived products, traceability, allergen and bioactive peptides identification, among others. Even thought, the full potential of proteomics is far from being fully exploited in forest tree research, with PTMs and interactomics being reserved to plant model systems. The most outstanding achievements in forest tree proteomics in the last decade as well as prospects are discussed.

Genome assembly and isoform analysis of a highly heterozygous New Zealand fisheries species, the tarakihi (Nemadactylus macropterus)

Although being some of the most valuable and heavily exploited wild organisms, few fisheries species have been studied at the whole-genome level. This is especially the case in New Zealand, where genomics resources are urgently needed to assist fisheries management. Here, we generated 55 Gb of short Illumina reads (92× coverage) and 73 Gb of long Nanopore reads (122×) to produce the first genome assembly of the marine teleost tarakihi [Nemadactylus macropterus (Forster, 1801)], a highly valuable fisheries species in New Zealand. An additional 300 Mb of Iso-Seq reads were obtained to assist in gene annotation. The final genome assembly was 568 Mb long with an N50 of 3.37 Mb. The genome completeness was high, with 97.8% of complete Actinopterygii Benchmarking Universal Single-Copy Orthologs. Heterozygosity values estimated through k-mer counting (1.00%) and bi-allelic SNPs (0.64%) were high compared with the same values reported for other fishes. Iso-Seq analysis recovered 91,313 unique transcripts from 15,515 genes (mean ratio of 5.89 transcripts per gene), and the most common alternative splicing event was intron retention. This highly contiguous genome assembly and the isoform-resolved transcriptome will provide a useful resource to assist the study of population genomics and comparative eco-evolutionary studies in teleosts and related organisms.

Distinct composition and amplification dynamics of transposable elements in sacred lotus (Nelumbo nucifera Gaertn.)

Sacred lotus (Nelumbo nucifera Gaertn.) is a basal eudicot plant with a unique lifestyle, physiological features, and evolutionary characteristics. Here we report the unique profile of transposable elements (TEs) in the genome, using a manually curated repeat library. TEs account for 59% of the genome, and hAT (Ac/Ds) elements alone represent 8%, more than in any other known plant genome. About 18% of the lotus genome is comprised of Copia LTR retrotransposons, and over 25% of them are associated with non-canonical termini (non-TGCA). Such high abundance of non-canonical LTR retrotransposons has not been reported for any other organism. TEs are very abundant in genic regions, with retrotransposons enriched in introns and DNA transposons primarily in flanking regions of genes. The recent insertion of TEs in introns has led to significant intron size expansion, with a total of 200 Mb in the 28 455 genes. This is accompanied by declining TE activity in intergenic regions, suggesting distinct control efficacy of TE amplification in different genomic compartments. Despite the prevalence of TEs in genic regions, some genes are associated with fewer TEs, such as those involved in fruit ripening and stress responses. Other genes are enriched with TEs, and genes in epigenetic pathways are the most associated with TEs in introns, indicating a dynamic interaction between TEs and the host surveillance machinery. The dramatic differential abundance of TEs with genes involved in different biological processes as well as the variation of target preference of different TEs suggests the composition and activity of TEs influence the path of evolution.

The genome of homosporous maidenhair fern sheds light on the euphyllophyte evolution and defences

Euphyllophytes encompass almost all extant plants, including two sister clades, ferns and seed plants. Decoding genomes of ferns is the key to deep insight into the origin of euphyllophytes and the evolution of seed plants. Here we report a chromosome-level genome assembly of Adiantum capillus-veneris L., a model homosporous fern. This fern genome comprises 30 pseudochromosomes with a size of 4.8-gigabase and a contig N50 length of 16.22 Mb. Gene co-expression network analysis uncovered that homospore development in ferns has relatively high genetic similarities with that of the pollen in seed plants. Analysing fern defence response expands understanding of evolution and diversity in endogenous bioactive jasmonates in plants. Moreover, comparing fern genomes with those of other land plants reveals changes in gene families important for the evolutionary novelties within the euphyllophyte clade. These results lay a foundation for studies on fern genome evolution and function, as well as the origin and evolution of euphyllophytes.

Conifer Biotechnology: An Overview

The peculiar characteristics of conifers determine the difficulty of their study and their great importance from various points of view. However, their study faces numerous important scientific, methodological, cultural, economic, social, and legal challenges. This paper presents an approach to several of those challenges and proposes a multidisciplinary scientific perspective that leads to a holistic understanding of conifers from the perspective of the latest technical, computer, and scientific advances. This review highlights the deep connection that all scientific contributions to conifers can have in each other as fully interrelated communicating vessels.

The Larix kaempferi genome reveals new insights into wood properties

Here, through single-molecule real-time sequencing, we present a high-quality genome sequence of the Japanese larch (Larix kaempferi), a conifer species with great value for wood production and ecological afforestation. The assembled genome is 10.97 Gb in size, harboring 45,828 protein-coding genes. Of the genome, 66.8% consists of repeat sequences, of which long terminal repeat retrotransposons are dominant and make up 69.86%. We find that tandem duplications have been responsible for the expansion of genes involved in transcriptional regulation and stress responses, unveiling their crucial roles in adaptive evolution. Population transcriptome analysis reveals that lignin content in L. kaempferi is mainly determined by the process of monolignol polymerization. The expression values of six genes (LkCOMT7, LkCOMT8, LkLAC23, LkLAC102, LkPRX148, and LkPRX166) have significantly positive correlations with lignin content. These results indicated that the increased expression of these six genes might be responsible for the high lignin content of the larches' wood. Overall, this study provides new genome resources for investigating the evolution and biological function of conifer trees, and also offers new insights into wood properties of larches.

Evaluating the accuracy of variant calling methods using the frequency of parent-offspring genotype mismatch

The use of next‐generation sequencing (NGS) data sets has increased dramatically over the last decade, but there have been few systematic analyses quantifying the accuracy of the commonly used variant caller programs. Here we used a familial design consisting of diploid tissue from a single lodgepole pine (Pinus contorta) parent and the maternally derived haploid tissue from 106 full‐sibling offspring, where mismatches could only arise due to mutation or bioinformatic error. Given the rarity of mutation, we used the rate of mismatches between parent and offspring genotype calls to infer the single nucleotide polymorphism (SNP) genotyping error rates of FreeBayes, HaplotypeCaller, SAMtools, UnifiedGenotyper, and VarScan. With baseline filtering HaplotypeCaller and UnifiedGenotyper yielded more SNPs and higher error rates by one to two orders of magnitude, whereas FreeBayes, SAMtools and VarScan yielded lower numbers of SNPs and more modest error rates. To facilitate comparison between variant callers we standardized each SNP set to the same number of SNPs using additional filtering, where UnifiedGenotyper consistently produced the smallest proportion of genotype errors, followed by HaplotypeCaller, VarScan, SAMtools, and FreeBayes. Additionally, we found that error rates were minimized for SNPs called by more than one variant caller. Finally, we evaluated the performance of various commonly used filtering metrics on SNP calling. Our analysis provides a quantitative assessment of the accuracy of five widely used variant calling programs and offers valuable insights into both the choice of variant caller program and the choice of filtering metrics, especially for researchers using non‐model study systems.

Genomic Prediction of Complex Traits in Perennial Plants: A Case for Forest Trees

This chapter provides an overview of the genomic selection progress in long-lived forest tree species. Factors affecting the prediction accuracy in genomic prediction are assessed with examples from empirical studies. Infrastructure and resources required for the implementation of genomic selection are evaluated. Some general guidelines are provided for the successful application of genomic selection in forest tree breeding programs.

MicroRNA-mediated post-transcriptional regulation of Pinus pinaster response and resistance to pinewood nematode

Pine wilt disease (PWD), caused by the parasitic nematode Bursaphelenchus xylophilus, or pinewood nematode (PWN), is a serious threat to pine forests in Europe. Pinus pinaster is highly susceptible to the disease and it is currently the most affected European pine species. In this work, we investigated the role of small RNAs (sRNAs) in regulating P. pinaster–PWN interaction in an early stage of infection. After performing an artificial PWN inoculation assay, we have identified 105 plant microRNAs (miRNAs) responsive to PWN. Based on their predicted targets, part of these miRNAs was associated with roles in jasmonate-response pathway, ROS detoxification, and terpenoid biosynthesis. Furthermore, by comparing resistant and susceptible plants, eight miRNAs with putative functions in plant defence and resistance to PWN have been identified. Finally, we explored the possibility of bidirectional trans-kingdom RNA silencing, identifying several P. pinaster genes putatively targeted by PWN miRNAs, which was supported by degradome analysis. Targets for P. pinaster miRNAs were also predicted in PWN, suggesting a role for trans-kingdom miRNA transfer and gene silencing both in PWN parasitism as in P. pinaster resistance to PWD. Our results provide new insights into previously unexplored roles of sRNA post-transcriptional regulation in P. pinaster response and resistance to PWN.

Identification and Tissue-Specific Expression Analysis of CYP720B Subfamily Genes in Slash Pine and Loblolly Pine

Diterpene resin acids (DRAs) are major components of pine oleoresin that can effectively resist the invasion of insects and pathogenic microorganisms. The subfamily of cytochrome P450s, CYP720B, catalyzes diterpene products into DRAs. Identifying CYP720B subfamily members and revealing the characteristics of tissue-specific expression would help understand diterpene-rich structures and diverse types. Slash pine and loblolly pine are important pines that provide oleoresin products. In this study, we identified CYP720B candidate genes based on the Pinus taeda V2.0 genome and full-length transcriptome of slash pine by PacBio. A total of 17 genes in slash pine and 19 in loblolly pine were identified and classified into four main clades by phylogenetic analysis. An analysis of cis-acting elements showed that CYP720B genes were closely related to adversity resistance. The gene expression of these candidates in different tissues was quantified by real-time quantitative PCR (RT–qPCR) analysis. Most of the genes showed relatively higher expression levels in roots and stems than in the other tissues, corresponding with the results of DRA component detection by gas chromatography–mass spectrometry (GC–MS), which indicated that stems and roots might be important tissues in oleoresin biosynthesis. These results provide a valuable resource for a better understanding of the biological role of individual CYP720Bs in slash pine and loblolly pine.

Genome-Wide Prediction of Transcription Start Sites in Conifers

The identification of promoters is an essential step in the genome annotation process, providing a framework for gene regulatory networks and their role in transcription regulation. Despite considerable advances in the high-throughput determination of transcription start sites (TSSs) and transcription factor binding sites (TFBSs), experimental methods are still time-consuming and expensive. Instead, several computational approaches have been developed to provide fast and reliable means for predicting the location of TSSs and regulatory motifs on a genome-wide scale. Numerous studies have been carried out on the regulatory elements of mammalian genomes, but plant promoters, especially in gymnosperms, have been left out of the limelight and, therefore, have been poorly investigated. The aim of this study was to enhance and expand the existing genome annotations using computational approaches for genome-wide prediction of TSSs in the four conifer species: loblolly pine, white spruce, Norway spruce, and Siberian larch. Our pipeline will be useful for TSS predictions in other genomes, especially for draft assemblies, where reliable TSS predictions are not usually available. We also explored some of the features of the nucleotide composition of the predicted promoters and compared the GC properties of conifer genes with model monocot and dicot plants. Here, we demonstrate that even incomplete genome assemblies and partial annotations can be a reliable starting point for TSS annotation. The results of the TSS prediction in four conifer species have been deposited in the Persephone genome browser, which allows smooth visualization and is optimized for large data sets. This work provides the initial basis for future experimental validation and the study of the regulatory regions to understand gene regulation in gymnosperms.

Identifying and testing marker-trait associations for growth and phenology in three pine species: Implications for genomic prediction

In tree species, genomic prediction offers the potential to forecast mature trait values in early growth stages, if robust marker-trait associations can be identified. Here we apply a novel multispecies approach using genotypes from a new genotyping array, based on 20,795 single nucleotide polymorphisms (SNPs) from three closely related pine species (Pinus sylvestris, Pinus uncinata and Pinus mugo), to test for associations with growth and phenology data from a common garden study. Predictive models constructed using significantly associated SNPs were then tested and applied to an independent multisite field trial of P. sylvestris and the capability to predict trait values was evaluated. One hundred and eighteen SNPs showed significant associations with the traits in the pine species. Common SNPs (MAF > 0.05) associated with bud set were only found in genes putatively involved in growth and development, whereas those associated with growth and budburst were also located in genes putatively involved in response to environment and, to a lesser extent, reproduction. At one of the two independent sites, the model we developed produced highly significant correlations between predicted values and observed height data (YA, height 2020: r = 0.376, p < 0.001). Predicted values estimated with our budburst model were weakly but positively correlated with duration of budburst at one of the sites (GS, 2015: r = 0.204, p = 0.034; 2018: r = 0.205, p = 0.034-0.037) and negatively associated with budburst timing at the other (YA: r = -0.202, p = 0.046). Genomic prediction resulted in the selection of sets of trees whose mean height was taller than the average for each site. Our results provide tentative support for the capability of prediction models to forecast trait values in trees, while highlighting the need for caution in applying them to trees grown in different environments.

Achievements and Challenges of Genomics-Assisted Breeding in Forest Trees: From Marker-Assisted Selection to Genome Editing

Forest tree breeding efforts have focused mainly on improving traits of economic importance, selecting trees suited to new environments or generating trees that are more resilient to biotic and abiotic stressors. This review describes various methods of forest tree selection assisted by genomics and the main technological challenges and achievements in research at the genomic level. Due to the long rotation time of a forest plantation and the resulting long generation times necessary to complete a breeding cycle, the use of advanced techniques with traditional breeding have been necessary, allowing the use of more precise methods for determining the genetic architecture of traits of interest, such as genome-wide association studies (GWASs) and genomic selection (GS). In this sense, main factors that determine the accuracy of genomic prediction models are also addressed. In turn, the introduction of genome editing opens the door to new possibilities in forest trees and especially clustered regularly interspaced short palindromic repeats and CRISPR-associated protein 9 (CRISPR/Cas9). It is a highly efficient and effective genome editing technique that has been used to effectively implement targetable changes at specific places in the genome of a forest tree. In this sense, forest trees still lack a transformation method and an inefficient number of genotypes for CRISPR/Cas9. This challenge could be addressed with the use of the newly developing technique GRF-GIF with speed breeding.

A genome-wide SNP genotyping resource for tropical pine tree species

We performed gene and genome targeted SNP discovery towards the development of a genome-wide, multi-species genotyping array for tropical pines. Pooled RNA-seq data from shoots of seedlings from five tropical pine species was used to identify transcript-based SNPs resulting in 1.3 million candidate Affymetrix SNP probe sets. In addition, we used a custom 40K probe set to perform capture-seq in pooled DNA from 81 provenances representing the natural ranges of six tropical pine species in Mexico and Central America resulting in 563K candidate SNP probe sets. Altogether, 300K RNA-seq (72%) and 120K capture-seq (28%) derived SNP probe sets were tiled on a 420K screening array that was used to genotype 576 trees representing the 81 provenances and commercial breeding material. Based on the screening array results, 50K SNPs were selected for commercial SNP array production (Axiom 384 format, Thermo Fisher Scientific) including 20K polymorphic SNPs for P. patula, P. tecunumanii, P. oocarpa and P. caribaea, 15K for P. greggii and P. maximinoi, 13K for P. elliottii and 8K for P. pseudostrobus. We included 9.7K ancestry informative SNPs that will be valuable for species and hybrid discrimination. Of the 50K SNP markers, 25% are polymorphic in only one species, while 75% are shared by two or more species. The Pitro50K SNP chip will be useful for population genomics and molecular breeding in this group of pine species that, together with their hybrids, represent the majority of fast-growing tropical and subtropical pine plantations globally.

Sequencing and Chromosome-Scale Assembly of Plant Genomes, Brassica rapa as a Use Case

With the rise of long-read sequencers and long-range technologies, delivering high-quality plant genome assemblies is no longer reserved to large consortia. Not only sequencing techniques, but also computer algorithms have reached a point where the reconstruction of assemblies at the chromosome scale is now feasible at the laboratory scale. Current technologies, in particular long-range technologies, are numerous, and selecting the most promising one for the genome of interest is crucial to obtain optimal results. In this study, we resequenced the genome of the yellow sarson, Brassica rapa cv. Z1, using the Oxford Nanopore PromethION sequencer and assembled the sequenced data using current assemblers. To reconstruct complete chromosomes, we used and compared three long-range scaffolding techniques, optical mapping, Omni-C, and Pore-C sequencing libraries, commercialized by Bionano Genomics, Dovetail Genomics, and Oxford Nanopore Technologies, respectively, or a combination of the three, in order to evaluate the capability of each technology.

Evolutionary history of the mediterranean Pinus halepensis-brutia species complex using gene-resequencing and transcriptomic approaches

Complementary gene-resequencing and transcriptomic approaches reveal contrasted evolutionary histories in a species complex. Pinus halepensis and Pinus brutia are closely related species that can intercross, but occupy different geographical ranges and bioclimates. To study the evolution of this species complex and to provide genomic resources for further research, we produce and analyze two new complementary sets of genetic resources: (i) a set of 172 re-sequenced genomic target loci analyzed in 45 individuals, and (ii) a set of 11 transcriptome assemblies. These two datasets provide insights congruent with previous studies: P. brutia displays high level of genetic diversity and no genetic sub-structure, while P. halepensis shows three main genetic clusters, the western Mediterranean and North African clusters displaying much lower genetic diversity than the eastern Mediterranean cluster, the latter cluster having similar genetic diversity to P. brutia. In addition, these datasets provide new insights on the timing of the species-complex history: the two species would have split at the end of the tertiary, and the changing climatic conditions of the Mediterranean region at the end of the Tertiary-beginning of the Quaternary, together with the distinct species tolerance to harsh climatic conditions would have resulted in different geographic distributions, demographic histories and genetic patterns of the two pines. The multiple glacial-interglacial cycles during the Quaternary would have led to the expansion of P. brutia in the Middle East, while P. halepensis would have been through bottlenecks. The last glaciations, from 0.6 Mya on, would have affected further the Western genetic pool of P. halepensis.

Toward genomic selection in Pinus taeda: Integrating resources to support array design in a complex conifer genome

PREMISE

An informatics approach was used for the construction of an Axiom genotyping array from heterogeneous, high-throughput sequence data to assess the complex genome of loblolly pine (Pinus taeda).

METHODS

High-throughput sequence data, sourced from exome capture and whole genome reduced-representation approaches from 2698 trees across five sequence populations, were analyzed with the improved genome assembly and annotation for the loblolly pine. A variant detection, filtering, and probe design pipeline was developed to detect true variants across and within populations. From 8.27 million variants, a total of 642,275 were evaluated and 423,695 of those were screened across a range-wide population.

RESULTS

The final informatics and screening approach delivered an Axiom array representing 46,439 high-confidence variants to the forest tree breeding and genetics community. Based on the annotated reference genome, 34% were located in or directly upstream or downstream of genic regions.

DISCUSSION

The Pita50K array represents a genome-wide resource developed from sequence data for an economically important conifer, loblolly pine. It uniquely integrates independent projects that assessed trees sampled across the native range. The challenges associated with the large and repetitive genome are addressed in the development of this resource.

Breeding for Climate Change Resilience: A Case Study of Loblolly Pine (Pinus taeda L.) in North America

Earth's atmosphere is warming and the effects of climate change are becoming evident. A key observation is that both the average levels and the variability of temperature and precipitation are changing. Information and data from new technologies are developing in parallel to provide multidisciplinary opportunities to address and overcome the consequences of these changes in forest ecosystems. Changes in temperature and water availability impose multidimensional environmental constraints that trigger changes from the molecular to the forest stand level. These can represent a threat for the normal development of the tree from early seedling recruitment to adulthood both through direct mortality, and by increasing susceptibility to pathogens, insect attack, and fire damage. This review summarizes the strengths and shortcomings of previous work in the areas of genetic variation related to cold and drought stress in forest species with particular emphasis on loblolly pine (Pinus taeda L.), the most-planted tree species in North America. We describe and discuss the implementation of management and breeding strategies to increase resilience and adaptation, and discuss how new technologies in the areas of engineering and genomics are shaping the future of phenotype-genotype studies. Lessons learned from the study of species important in intensively-managed forest ecosystems may also prove to be of value in helping less-intensively managed forest ecosystems adapt to climate change, thereby increasing the sustainability and resilience of forestlands for the future.

Candidate Genes for the High-Altitude Adaptations of Two Mountain Pine Taxa

Mountain plants, challenged by vegetation time contractions and dynamic changes in environmental conditions, developed adaptations that help them to balance their growth, reproduction, survival, and regeneration. However, knowledge regarding the genetic basis of species adaptation to higher altitudes remain scarce for most plant species. Here, we attempted to identify such corresponding genomic regions of high evolutionary importance in two closely related European pines, Pinus mugo and P. uncinata, contrasting them with a reference lowland relative—P. sylvestris. We genotyped 438 samples at thousands of single nucleotide polymorphism (SNP) markers, tested their genetic differentiation and population structure followed by outlier detection and gene ontology annotations. Markers clearly differentiated the species and uncovered patterns of population structure in two of them. In P. uncinata three Pyrenean sites were grouped together, while two outlying populations constituted a separate cluster. In P. sylvestris, Spanish population appeared distinct from the remaining four European sites. Between mountain pines and the reference species, 35 candidate genes for altitude-dependent selection were identified, including such encoding proteins responsible for photosynthesis, photorespiration and cell redox homeostasis, regulation of transcription, and mRNA processing. In comparison between two mountain pines, 75 outlier SNPs were found in proteins involved mainly in the gene expression and metabolism.

Mass spectrometry-based forest tree metabolomics

Research in forest tree species has advanced slowly when compared with other agricultural crops and model organisms, mainly due to the long-life cycles, large genome sizes, and lack of genomic tools. Additionally, trees are complex matrices, and the presence of interferents (e.g., oleoresins and cellulose) challenges the analysis of tree tissues with mass spectrometry (MS)-based analytical platforms. In this review, advances in MS-based forest tree metabolomics are discussed. Given their economic and ecological significance, particular focus is given to Pinus, Quercus, and Eucalyptus forest tree species to better understand their metabolite responses to abiotic and biotic stresses in the current climate change scenario. Furthermore, MS-based metabolomics technologies produce large and complex datasets that require expertize to adequately manage, process, analyze, and store the data in dedicated repositories. To ensure that the full potential of forest tree metabolomics data are translated into new knowledge, these data should comply with the FAIR principles (i.e., Findable, Accessible, Interoperable, and Re-usable). It is essential that adequate standards are implemented to annotate metadata from forest tree metabolomics studies as is already required by many science and governmental agencies and some major scientific publishers. © 2019 John Wiley & Sons Ltd. Mass Spec Rev 40:126-157, 2021.

Transcript Dynamics in Wounded and Inoculated Scots Pine

Comparative transcriptome analysis provides a useful tool for the exploration of plant-pathogen interaction by allowing in-depth comparison of gene expression between unaffected, inoculated and wounded organisms. Here we present the results of comparative transcriptome analysis in genetically identical one-year-old Scots pine ramets after wounding and inoculation with Heterobasidion annosum. We identified 230 genes that were more than 2-fold upregulated in inoculated samples (compared to controls) and 116 downregulated genes. Comparison of inoculated samp les with wounded samples identified 32 differentially expressed genes (30 were upregulated after inoculation). Several of the genes upregulated after inoculation are involved in protection from oxidative stress, while genes involved in photosynthesis, water transport and drought stress tolerance were downregulated. An NRT3 family protein was the most upregulated transcript in response to both inoculation and wounding, while a U-box domain-containing protein gene was the most upregulated gene comparing inoculation to wounding. The observed transcriptome dynamics suggest involvement of auxin, ethylene, jasmonate, gibberellin and reactive oxygen species pathways and cell wall modification regulation in response to H. annosum infection. The results are compared to methyl jasmonate induced transcriptome dynamics.

Comparative Study of Pine Reference Genomes Reveals Transposable Element Interconnected Gene Networks

Sequencing the giga-genomes of several pine species has enabled comparative genomic analyses of these outcrossing tree species. Previous studies have revealed the wide distribution and extraordinary diversity of transposable elements (TEs) that occupy the large intergenic spaces in conifer genomes. In this study, we analyzed the distribution of TEs in gene regions of the assembled genomes of Pinus taeda and Pinus lambertiana using high-performance computing resources. The quality of draft genomes and the genome annotation have significant consequences for the investigation of TEs and these aspects are discussed. Several TE families frequently inserted into genes or their flanks were identified in both species’ genomes. Potentially important sequence motifs were identified in TEs that could bind additional regulatory factors, promoting gene network formation with faster or enhanced transcription initiation. Node genes that contain many TEs were observed in multiple potential transposable element-associated networks. This study demonstrated the increased accumulation of TEs in the introns of stress-responsive genes of pines and suggests the possibility of rewiring them into responsive networks and sub-networks interconnected with node genes containing multiple TEs. Many such regulatory influences could lead to the adaptive environmental response clines that are characteristic of naturally spread pine populations.

Molecular study of drought response in the Mediterranean conifer Pinus pinaster Ait.: Differential transcriptomic profiling reveals constitutive water deficit-independent drought tolerance mechanisms

Adaptation of long-living forest trees to respond to environmental changes is essential to secure their performance under adverse conditions. Water deficit is one of the most significant stress factors determining tree growth and survival. Maritime pine (Pinus pinaster Ait.), the main source of softwood in southwestern Europe, is subjected to recurrent drought periods which, according to climate change predictions for the years to come, will progressively increase in the Mediterranean region. The mechanisms regulating pine adaptive responses to environment are still largely unknown. The aim of this work was to go a step further in understanding the molecular mechanisms underlying maritime pine response to water stress and drought tolerance at the whole plant level. A global transcriptomic profiling of roots, stems, and needles was conducted to analyze the performance of siblings showing contrasted responses to water deficit from an ad hoc designed full-sib family. Although P. pinaster is considered a recalcitrant species for vegetative propagation in adult phase, the analysis was conducted using vegetatively propagated trees exposed to two treatments: well-watered and moderate water stress. The comparative analyses led us to identify organ-specific genes, constitutively expressed as well as differentially expressed when comparing control versus water stress conditions, in drought-sensitive and drought-tolerant genotypes. Different response strategies can point out, with tolerant individuals being pre-adapted for coping with drought by constitutively expressing stress-related genes that are detected only in latter stages on sensitive individuals subjected to drought.

Comparative Transcriptome Analysis of Pine Trees Treated with Resistance-Inducing Substances against the Nematode Bursaphelenchus xylophilus

The pinewood nematode (PWN) Bursaphelenchus xylophilus causes pine wilt disease, which results in substantial economic and environmental losses across pine forests worldwide. Although systemic acquired resistance (SAR) is effective in controlling PWN, the detailed mechanisms underlying the resistance to PWN are unclear. Here, we treated pine samples with two SAR elicitors, acibenzolar-S-methyl (ASM) and methyl salicylic acid (MeSA) and constructed an in vivo transcriptome of PWN-infected pines under SAR conditions. A total of 252 million clean reads were obtained and mapped onto the reference genome. Compared with untreated pines, 1091 and 1139 genes were differentially upregulated following the ASM and MeSA treatments, respectively. Among these, 650 genes showed co-expression patterns in response to both SAR elicitors. Analysis of these patterns indicated a functional linkage among photorespiration, peroxisome, and glycine metabolism, which may play a protective role against PWN infection-induced oxidative stress. Further, the biosynthesis of flavonoids, known to directly control parasitic nematodes, was commonly upregulated under SAR conditions. The ASM- and MeSA-specific expression patterns revealed functional branches for myricetin and quercetin production in flavonol biosynthesis. This study will enhance the understanding of the dynamic interactions between pine hosts and PWN under SAR conditions.

Development of a single nucleotide polymorphism array for population genomic studies in four European pine species

Pines are some of the most ecologically and economically important tree species in the world, and many have enormous natural distributions or have been extensively planted. However, a lack of rapid genotyping capability is hampering progress in understanding the molecular basis of genetic variation in these species. Here, we deliver an efficient tool for genotyping thousands of single nucleotide polymorphism (SNP) markers across the genome that can be applied to genetic studies in pines. Polymorphisms from resequenced candidate genes and transcriptome sequences of P. sylvestris, P. mugo, P. uncinata, P. uliginosa and P. radiata were used to design a 49,829 SNP array (Axiom_PineGAP, Thermo Fisher). Over a third (34.68%) of the unigenes identified from the P. sylvestris transcriptome were represented on the array, which was used to screen samples of four pine species. The conversion rate for the array on all samples was 42% (N = 20,795 SNPs) and was similar for SNPs sourced from resequenced candidate gene and transcriptome sequences. The broad representation of gene ontology terms by unigenes containing converted SNPs reflected their coverage across the full transcriptome. Over a quarter of successfully converted SNPs were polymorphic among all species, and the data were successful in discriminating among the species and some individual populations. The SNP array provides a valuable new tool to advance genetic studies in these species and demonstrates the effectiveness of the technology for rapid genotyping in species with large and complex genomes.

Two Independent Plastid accD Transfers to the Nuclear Genome of Gnetum and Other Insights on Acetyl-CoA Carboxylase Evolution in Gymnosperms

Acetyl-CoA carboxylase (ACCase) is the key regulator of fatty acid biosynthesis. In most plants, ACCase exists in two locations (cytosol and plastids) and in two forms (homomeric and heteromeric). Heteromeric ACCase comprises four subunits, three of them (ACCA-C) are nuclear encoded (nr) and the fourth (ACCD) is usually plastid encoded. Homomeric ACCase is encoded by a single nr-gene (ACC). We investigated the ACCase gene evolution in gymnosperms by examining the transcriptomes of newly sequenced Gnetum ula, combined with 75 transcriptomes and 110 plastomes of other gymnosperms. AccD-coding sequences are elongated through the insertion of repetitive DNA in four out of five cupressophyte families (except Sciadopityaceae) and were functionally transferred to the nucleus of gnetophytes and Sciadopitys. We discovered that, among the three genera of gnetophytes, only Gnetum has two copies of nr-accD. Furthermore, using protoplast transient expression assays, we experimentally verified that the nr-accD precursor proteins in Gnetum and Sciadopitys can be delivered to the plastids. Of the two nr-accD copies of Gnetum, one dually targets plastids and mitochondria, whereas the other potentially targets plastoglobuli. The distinct transit peptides, gene architectures, and flanking sequences between the two Gnetum accDs suggest that they have independent origins. Our findings are the first account of two distinctly targeted nr-accDs of any green plants and the most comprehensive analyses of ACCase evolution in gymnosperms to date.

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

BACKGROUND

Pine wilt disease (PWD), which is caused by the pine wood nematode (PWN) Bursaphelenchus xylophilus, is currently the greatest threat to pine forests in Europe and East Asian countries including Japan. Constructing a detailed linkage map of DNA markers and identifying PWD resistance genes/loci lead to improved resistance in Pinus thunbergii, as well as other Pinus species that are also susceptible to PWD.

RESULTS

A total F₁ mapping population of 188 individuals derived from a cross between the PWD-resistant P. thunbergii varieties 'Tanabe 54' (resistant rank 2 to PWD) and 'Tosashimizu 63' (resistant rank 4 to PWD) was inoculated with PWN, and was evaluated for disease symptoms. To perform linkage analysis for PWN resistance, a set of three maps was constructed; two parental maps generated using the integrated two-way pseudo-testcross method, and a consensus map with population-type cross-pollination. The linkage map of 'Tanabe 54' consisted of 167 loci, and covered 14 linkage groups (LGs), with a total genetic distance of 1214.6 cM. The linkage map of 'Tosashimizu 63' consisted of 252 loci, and covered 14 LGs, with a total genetic distance of 1422.1 cM. The integrated consensus map comprised 12 LGs with the basic chromosome number of P. thunbergii, and a total genetic distance of 1403.6 cM. Results from quantitative trait loci (QTL) analysis using phenotype data and linkage maps indicated that PWN resistance is controlled by a single dominant allele, which was derived from the 'Tanabe 54' female parent. This major QTL was located on linkage group 3 and was designated PWD1 for PINE WILT DISEASE 1.

CONCLUSIONS

The PWD1 locus is a major resistance QTL located on the Pinus consensus LG03 that acts in a dominant manner to confer pine wood nematode resistance. Information from the present study will be useful for P. thunbergii breeding programs to improve resistance to PWD, and also to help identify susceptibility genes in Pinus species.

Leger J, Le J, Natarajan A, Rivera M, Guillergan M, Jaeger E, Steffy B, Zimin A. New de novo assembly of the Atlantic bottlenose dolphin (Tursiops truncatus) improves genome completeness and provides haplotype phasing

High-quality genomes are essential to resolve challenges in breeding, comparative biology, medicine, and conservation planning. New library preparation techniques along with better assembly algorithms result in continued improvements in assemblies for non-model organisms, moving them toward reference-quality genomes. We report on the latest genome assembly of the Atlantic bottlenose dolphin, leveraging Illumina sequencing data coupled with a combination of several library preparation techniques. These include Linked-Reads (Chromium, 10x Genomics), mate pairs (MP), long insert paired ends, and standard paired end. Data were assembled with the commercial DeNovoMAGIC assembly software, resulting in two assemblies, a traditional "haploid" assembly (Tur_tru_Illumina_hap_v1) that is a mosaic of the two parental haplotypes and a phased assembly (Tur_tru_Illumina_phased_v1) where each scaffold has sequence from a single homologous chromosome. We show that Tur_tru_Illumina_hap_v1 is more complete and more accurate compared to the current best reference based on the amount and composition of sequence, the consistency of the MP alignments to the assembled scaffolds, and on the analysis of conserved single-copy mammalian orthologs. The phased de novo assembly Tur_tru_Illumina_phased_v1 is the first publicly available for this species and provides the community with novel and accurate ways to explore the heterozygous nature of the dolphin genome.

Stepwise large genome assembly approach: a case of Siberian larch (Larix sibirica Ledeb)

Background

De novo assembling of large genomes, such as in conifers (~ 12–30 Gbp), which also consist of ~ 80% of repetitive DNA, is a very complex and computationally intense endeavor. One of the main problems in assembling such genomes lays in computing limitations of nucleotide sequence assembly programs (DNA assemblers). As a rule, modern assemblers are usually designed to assemble genomes with a length not exceeding the length of the human genome (3.24 Gbp). Most assemblers cannot handle the amount of input sequence data required to provide sufficient coverage needed for a high-quality assembly.

Results

An original stepwise method of de novo assembly by parts (sets), which allows to bypass the limitations of modern assemblers associated with a huge amount of data being processed, is presented in this paper. The results of numerical assembling experiments conducted using the model plant Arabidopsis thaliana, Prunus persica (peach) and four most popular assemblers, ABySS, SOAPdenovo, SPAdes, and CLC Assembly Cell, showed the validity and effectiveness of the proposed stepwise assembling method.

Conclusion

Using the new stepwise de novo assembling method presented in the paper, the genome of Siberian larch, Larix sibirica Ledeb. (12.34 Gbp) was completely assembled de novo by the CLC Assembly Cell assembler. It is the first genome assembly for larch species in addition to only five other conifer genomes sequenced and assembled for Picea abies, Picea glauca, Pinus taeda, Pinus lambertiana, and Pseudotsuga menziesii var. menziesii.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2570-y) contains supplementary material, which is available to authorized users.

Finding Nemo: hybrid assembly with Oxford Nanopore and Illumina reads greatly improves the clownfish (Amphiprion ocellaris) genome assembly

Background

Some of the most widely recognized coral reef fishes are clownfish or anemonefish, members of the family Pomacentridae (subfamily: Amphiprioninae). They are popular aquarium species due to their bright colours, adaptability to captivity, and fascinating behavior. Their breeding biology (sequential hermaphrodites) and symbiotic mutualism with sea anemones have attracted much scientific interest. Moreover, there are some curious geographic-based phenotypes that warrant investigation. Leveraging on the advancement in Nanopore long read technology, we report the first hybrid assembly of the clown anemonefish (Amphiprion ocellaris) genome utilizing Illumina and Nanopore reads, further demonstrating the substantial impact of modest long read sequencing data sets on improving genome assembly statistics.

Results

We generated 43 Gb of short Illumina reads and 9 Gb of long Nanopore reads, representing approximate genome coverage of 54× and 11×, respectively, based on the range of estimated k-mer-predicted genome sizes of between 791 and 967 Mbp. The final assembled genome is contained in 6404 scaffolds with an accumulated length of 880 Mb (96.3% BUSCO-calculated genome completeness). Compared with the Illumina-only assembly, the hybrid approach generated 94% fewer scaffolds with an 18-fold increase in N₅₀ length (401 kb) and increased the genome completeness by an additional 16%. A total of 27 240 high-quality protein-coding genes were predicted from the clown anemonefish, 26 211 (96%) of which were annotated functionally with information from either sequence homology or protein signature searches.

Conclusions

We present the first genome of any anemonefish and demonstrate the value of low coverage (∼11×) long Nanopore read sequencing in improving both genome assembly contiguity and completeness. The near-complete assembly of the A. ocellaris genome will be an invaluable molecular resource for supporting a range of genetic, genomic, and phylogenetic studies specifically for clownfish and more generally for other related fish species of the family Pomacentridae.

Independent assessment and improvement of wheat genome sequence assemblies using Fosill jumping libraries

Background

The accurate sequencing and assembly of very large, often polyploid, genomes remains a challenging task, limiting long-range sequence information and phased sequence variation for applications such as plant breeding. The 15-Gb hexaploid bread wheat (Triticum aestivum) genome has been particularly challenging to sequence, and several different approaches have recently generated long-range assemblies. Mapping and understanding the types of assembly errors are important for optimising future sequencing and assembly approaches and for comparative genomics.

Results

Here we use a Fosill 38-kb jumping library to assess medium and longer–range order of different publicly available wheat genome assemblies. Modifications to the Fosill protocol generated longer Illumina sequences and enabled comprehensive genome coverage. Analyses of two independent Bacterial Artificial Chromosome (BAC)-based chromosome-scale assemblies, two independent Illumina whole genome shotgun assemblies, and a hybrid Single Molecule Real Time (SMRT-PacBio) and short read (Illumina) assembly were carried out. We revealed a surprising scale and variety of discrepancies using Fosill mate-pair mapping and validated several of each class. In addition, Fosill mate-pairs were used to scaffold a whole genome Illumina assembly, leading to a 3-fold increase in N50 values.

Conclusions

Our analyses, using an independent means to validate different wheat genome assemblies, show that whole genome shotgun assemblies based solely on Illumina sequences are significantly more accurate by all measures compared to BAC-based chromosome-scale assemblies and hybrid SMRT-Illumina approaches. Although current whole genome assemblies are reasonably accurate and useful, additional improvements will be needed to generate complete assemblies of wheat genomes using open-source, computationally efficient, and cost-effective methods.

Improved de novo genome assembly and analysis of the Chinese cucurbit Siraitia grosvenorii, also known as monk fruit or luo-han-guo

Background

Luo-han-guo (Siraitia grosvenorii), also called monk fruit, is a member of the Cucurbitaceae family. Monk fruit has become an important area for research because of the pharmacological and economic potential of its noncaloric, extremely sweet components (mogrosides). It is also commonly used in traditional Chinese medicine for the treatment of lung congestion, sore throat, and constipation. Recently, a single reference genome became available for monk fruit, assembled from 36.9x genome coverage reads via Illumina sequencing platforms. This genome assembly has a relatively short (34.2 kb) contig N50 length and lacks integrated annotations. These drawbacks make it difficult to use as a reference in assembling transcriptomes and discovering novel functional genes.

Findings

Here, we offer a new high-quality draft of the S. grosvenorii genome assembled using 31 Gb (∼73.8x) long single molecule real time sequencing reads and polished with ∼50 Gb Illumina paired-end reads. The final genome assembly is approximately 469.5 Mb, with a contig N50 length of 432,384 bp, representing a 12.6-fold improvement. We further annotated 237.3 Mb of repetitive sequence and 30,565 consensus protein coding genes with combined evidence. Phylogenetic analysis showed that S. grosvenorii diverged from members of the Cucurbitaceae family approximately 40.9 million years ago. With comprehensive transcriptomic analysis and differential expression testing, we identified 4,606 up-regulated genes in the early fruit compared to the leaf, a number of which were linked to metabolic pathways regulating fruit development and ripening.

Conclusions

The availability of this new monk fruit genome assembly, as well as the annotations, will facilitate the discovery of new functional genes and the genetic improvement of monk fruit.

Exploring the genome and transcriptome of the cave nectar bat Eonycteris spelaea with PacBio long-read sequencing

Background

In the past two decades, bats have emerged as an important model system to study host-pathogen interactions. More recently, it has been shown that bats may also serve as a new and excellent model to study aging, inflammation, and cancer, among other important biological processes. The cave nectar bat or lesser dawn bat (Eonycteris spelaea) is known to be a reservoir for several viruses and intracellular bacteria. It is widely distributed throughout the tropics and subtropics from India to Southeast Asia and pollinates several plant species, including the culturally and economically important durian in the region. Here, we report the whole-genome and transcriptome sequencing, followed by subsequent de novo assembly, of the E. spelaea genome solely using the Pacific Biosciences (PacBio) long-read sequencing platform.

Findings

The newly assembled E. spelaea genome is 1.97 Gb in length and consists of 4,470 sequences with a contig N50 of 8.0 Mb. Identified repeat elements covered 34.65% of the genome, and 20,640 unique protein-coding genes with 39,526 transcripts were annotated.

Conclusions

We demonstrated that the PacBio long-read sequencing platform alone is sufficient to generate a comprehensive de novo assembled genome and transcriptome of an important bat species. These results will provide useful insights and act as a resource to expand our understanding of bat evolution, ecology, physiology, immunology, viral infection, and transmission dynamics.

Genomic Variation Among and Within Six Juglans Species

Genomic analysis in Juglans (walnuts) is expected to transform the breeding and agricultural production of both nuts and lumber. To that end, we report here the determination of reference sequences for six additional relatives of Juglans regia: Juglans sigillata (also from section Dioscaryon), Juglans nigra, Juglans microcarpa, Juglans hindsii (from section Rhysocaryon), Juglans cathayensis (from section Cardiocaryon), and the closely related Pterocarya stenoptera. While these are ‘draft’ genomes, ranging in size between 640Mbp and 990Mbp, their contiguities and accuracies can support powerful annotations of genomic variation that are often the foundation of new avenues of research and breeding. We annotated nucleotide divergence and synteny by creating complete pairwise alignments of each reference genome to the remaining six. In addition, we have re-sequenced a sample of accessions from four Juglans species (including regia). The variation discovered in these surveys comprises a critical resource for experimentation and breeding, as well as a solid complementary annotation. To demonstrate the potential of these resources the structural and sequence variation in and around the polyphenol oxidase loci, PPO1 and PPO2 were investigated. As reported for other seed crops variation in this gene is implicated in the domestication of walnuts. The apparently Juglandaceae specific PPO1 duplicate shows accelerated divergence and an excess of amino acid replacement on the lineage leading to accessions of the domesticated nut crop species, Juglans regia and sigillata.

Exploring the loblolly pine (Pinus taeda L.) genome by BAC sequencing and Cot analysis

Loblolly pine (LP; Pinus taeda L.) is an economically and ecologically important tree in the southeastern U.S. To advance understanding of the loblolly pine (LP; Pinus taeda L.) genome, we sequenced and analyzed 100 BAC clones and performed a Cot analysis. The Cot analysis indicates that the genome is composed of 57, 24, and 10% highly-repetitive, moderately-repetitive, and single/low-copy sequences, respectively (the remaining 9% of the genome is a combination of fold back and damaged DNA). Although single/low-copy DNA only accounts for 10% of the LP genome, the amount of single/low-copy DNA in LP is still 14 times the size of the Arabidopsis genome. Since gene numbers in LP are similar to those in Arabidopsis, much of the single/low-copy DNA of LP would appear to be composed of DNA that is both gene- and repeat-poor. Macroarrays prepared from a LP bacterial artificial chromosome (BAC) library were hybridized with probes designed from cell wall synthesis/wood development cDNAs, and 50 of the "targeted" clones were selected for further analysis. An additional 25 clones were selected because they contained few repeats, while 25 more clones were selected at random. The 100 BAC clones were Sanger sequenced and assembled. Of the targeted BACs, 80% contained all or part of the cDNA used to target them. One targeted BAC was found to contain fungal DNA and was eliminated from further analysis. Combinations of similarity-based and ab initio gene prediction approaches were utilized to identify and characterize potential coding regions in the 99 BACs containing LP DNA. From this analysis, we identified 154 gene models (GMs) representing both putative protein-coding genes and likely pseudogenes. Ten of the GMs (all of which were specifically targeted) had enough support to be classified as intact genes. Interestingly, the 154 GMs had statistically indistinguishable (α = 0.05) distributions in the targeted and random BAC clones (15.18 and 12.61 GM/Mb, respectively), whereas the low-repeat BACs contained significantly fewer GMs (7.08 GM/Mb). However, when GM length was considered, the targeted BACs had a significantly greater percentage of their length in GMs (3.26%) when compared to random (1.63%) and low-repeat (0.62%) BACs. The results of our study provide insight into LP evolution and inform ongoing efforts to produce a reference genome sequence for LP, while characterization of genes involved in cell wall production highlights carbon metabolism pathways that can be leveraged for increasing wood production.

A guide to sequence your favorite plant genomes

With the rapid development of sequencing technology and the plummeting cost, assembling whole genomes from non-model plants will soon become routine for plant systematists and evolutionary biologists. Here we summarize and compare several of the latest genome sequencing and assembly approaches, offering a practical guide on how to approach a genome project. We also highlight certain precautions that need to be taken before investing time and money into a genome project.

ConTEdb: a comprehensive database of transposable elements in conifers

Conifers are the largest and most ubiquitous group of gymnosperms and have significant ecological significance and economic importance. However, the huge and complex genomes have hindered the sequencing and mining of conifer genomes. In this study, we identified 413 423 transposable elements (TEs) from Picea abies, Picea glauca and Pinus taeda using a combination of multiple approaches and classified them into 11 133 families. A comprehensive web-based database, ConTEdb, was constructed and served for researchers. ConTEdb enables users to browse, retrieve and download the TE sequences from the database. Several analysis tools are integrated into ConTEdb to help users mine the TE data easily and effectively. In summary, ConTEdb provides a platform to study TE biology and functional genomics in conifers.