Taming the massive genome of Scots pine with PiSy50k, a new genotyping array for conifer research.

Pinus sylvestris (Scots pine) is the most widespread coniferous tree in the boreal forests of Eurasia, with major economic and ecological importance. However, its large and repetitive genome presents a challenge for conducting genome-wide analyses such as association studies, genetic mapping and genomic selection. We present a new 50K single-nucleotide polymorphism (SNP) genotyping array for Scots pine research, breeding and other applications. To select the SNP set, we first genotyped 480 Scots pine samples on a 407 540 SNP screening array and identified 47 712 high-quality SNPs for the final array (called 'PiSy50k'). Here, we provide details of the design and testing, as well as allele frequency estimates from the discovery panel, functional annotation, tissue-specific expression patterns and expression level information for the SNPs or corresponding genes, when available. We validated the performance of the PiSy50k array using samples from Finland and Scotland. Overall, 39 678 (83.2%) SNPs showed low error rates (mean = 0.9%). Relatedness estimates based on array genotypes were consistent with the expected pedigrees, and the level of Mendelian error was negligible. In addition, array genotypes successfully discriminate between Scots pine populations of Finnish and Scottish origins. The PiSy50k SNP array will be a valuable tool for a wide variety of future genetic studies and forestry applications.

A key role for sex chromosomes in the regulation of parthenogenesis in the brown alga Ectocarpus.

Although evolutionary transitions from sexual to asexual reproduction are frequent in eukaryotes, the genetic bases of these shifts remain largely elusive. Here, we used classic quantitative trait analysis, combined with genomic and transcriptomic information to dissect the genetic basis of asexual, parthenogenetic reproduction in the brown alga Ectocarpus. We found that parthenogenesis is controlled by the sex locus, together with two additional autosomal loci, highlighting the key role of the sex chromosome as a major regulator of asexual reproduction. We identify several negative effects of parthenogenesis on male fitness, and different fitness effects of parthenogenetic capacity depending on the life cycle generation. Although allele frequencies in natural populations are currently unknown, we discuss the possibility that parthenogenesis may be under both sex-specific selection and generation/ploidally-antagonistic selection, and/or that the action of fluctuating selection on this trait may contribute to the maintenance of polymorphisms in populations. Importantly, our data provide the first empirical illustration, to our knowledge, of a trade-off between the haploid and diploid stages of the life cycle, where distinct parthenogenesis alleles have opposing effects on sexual and asexual reproduction and may help maintain genetic variation. These types of fitness trade-offs have profound evolutionary implications in natural populations and may structure life history evolution in organisms with haploid-diploid life cycles.

DISTAG/TBCCd1 Is Required for Basal Cell Fate Determination in Ectocarpus.

Brown algae are one of the most developmentally complex groups within the eukaryotes. As in many land plants and animals, their main body axis is established early in development, when the initial cell gives rise to two daughter cells that have apical and basal identities, equivalent to shoot and root identities in land plants, respectively. We show here that mutations in the Ectocarpus DISTAG (DIS) gene lead to loss of basal structures during both the gametophyte and the sporophyte generations. Several abnormalities were observed in the germinating initial cell in dis mutants, including increased cell size, disorganization of the Golgi apparatus, disruption of the microtubule network, and aberrant positioning of the nucleus. DIS encodes a TBCCd1 protein, which has a role in internal cell organization in animals, Chlamydomonas reinhardtii, and trypanosomes. Our study highlights the key role of subcellular events within the germinating initial cell in the determination of apical/basal cell identities in a brown alga and emphasizes the remarkable functional conservation of TBCCd1 in regulating internal cell organization across extremely distant eukaryotic groups.

Scots pine aminopropyltransferases shed new light on evolution of the polyamine biosynthesis pathway in seed plants.

Background and Aims: Polyamines are small metabolites present in all living cells and play fundamental roles in numerous physiological events in plants. The aminopropyltransferases (APTs), spermidine synthase (SPDS), spermine synthase (SPMS) and thermospermine synthase (ACL5), are essential enzymes in the polyamine biosynthesis pathway. In angiosperms, SPMS has evolved from SPDS via gene duplication, whereas in gymnosperms APTs are mostly unexplored and no SPMS gene has been reported. The present study aimed to investigate the functional properties of the SPDS and ACL5 proteins of Scots pine (Pinus sylvestris L.) in order to elucidate the role and evolution of APTs in higher plants. Methods: Germinating Scots pine seeds and seedlings were analysed for polyamines by high-performance liquid chromatography (HPLC) and the expression of PsSPDS and PsACL5 genes by in situ hybridization. Recombinant proteins of PsSPDS and PsACL5 were produced and investigated for functional properties. Also gene structures, promoter regions and phylogenetic relationships of PsSPDS and PsACL5 genes were analysed. Key Results: Scots pine tissues were found to contain spermidine, spermine and thermospermine. PsSPDS enzyme catalysed synthesis of both spermidine and spermine. PsACL5 was found to produce thermospermine, and PsACL5 gene expression was localized in the developing procambium in embryos and tracheary elements in seedlings. Conclusions: Contrary to previous views, our results demonstrate that SPMS activity is not a novel feature developed solely in the angiosperm lineage of seed plants but also exists as a secondary property in the Scots pine SPDS enzyme. The discovery of bifunctional SPDS from an evolutionarily old conifer reveals the missing link in the evolution of the polyamine biosynthesis pathway. The finding emphasizes the importance of pre-existing secondary functions in the evolution of new enzyme activities via gene duplication. Our results also associate PsACL5 with the development of vascular structures in Scots pine.

Re-annotation, improved large-scale assembly and establishment of a catalogue of noncoding loci for the genome of the model brown alga Ectocarpus.

The genome of the filamentous brown alga Ectocarpus was the first to be completely sequenced from within the brown algal group and has served as a key reference genome both for this lineage and for the stramenopiles. We present a complete structural and functional reannotation of the Ectocarpus genome. The large-scale assembly of the Ectocarpus genome was significantly improved and genome-wide gene re-annotation using extensive RNA-seq data improved the structure of 11 108 existing protein-coding genes and added 2030 new loci. A genome-wide analysis of splicing isoforms identified an average of 1.6 transcripts per locus. A large number of previously undescribed noncoding genes were identified and annotated, including 717 loci that produce long noncoding RNAs. Conservation of lncRNAs between Ectocarpus and another brown alga, the kelp Saccharina japonica, suggests that at least a proportion of these loci serve a function. Finally, a large collection of single nucleotide polymorphism-based markers was developed for genetic analyses. These resources are available through an updated and improved genome database. This study significantly improves the utility of the Ectocarpus genome as a high-quality reference for the study of many important aspects of brown algal biology and as a reference for genomic analyses across the stramenopiles.

High rate of adaptive evolution in two widespread European pines.

Comparing related organisms with differing ecological requirements and evolutionary histories can shed light on the mechanisms and drivers underlying genetic adaptation. Here, by examining a common set of hundreds of loci, we compare patterns of nucleotide diversity and molecular adaptation of two European conifers (Scots pine and maritime pine) living in contrasted environments and characterized by distinct population genetic structure (low and clinal in Scots pine, high and ecotypic in maritime pine) and demographic histories. We found higher nucleotide diversity in Scots pine than in maritime pine, whereas rates of new adaptive substitutions (ωa ), as estimated from the distribution of fitness effects, were similar across species and among the highest found in plants. Sample size and population genetic structure did not appear to have resulted in significant bias in estimates of ωa . Moreover, population contraction-expansion dynamics for each species did not affect differentially the rate of adaptive substitution in these two pines. Several methodological and biological factors may underlie the unusually high rate of adaptive evolution of Scots pine and maritime pine. By providing two new case studies with contrasting evolutionary histories, we contribute to disentangling the multiple factors potentially affecting adaptive evolution in natural plant populations.

Association of FLOWERING LOCUS T/TERMINAL FLOWER 1-like gene FTL2 expression with growth rhythm in Scots pine (Pinus sylvestris).

Understanding the genetic basis of the timing of bud set, an important trait in conifers, is relevant for adaptation and forestry practice. In common garden experiments, both Scots pine (Pinus sylvestris) and Norway spruce (Picea abies) show a latitudinal cline in the trait. We compared the regulation of their bud set biology by examining the expression of PsFTL2, a Pinus sylvestris homolog to PaFTL2, a FLOWERING LOCUS T/TERMINAL FLOWER 1 (FT/TFL1)-like gene, the expression levels of which have been found previously to be associated with the timing of bud set in Norway spruce. In a common garden study, we analyzed the relationship of bud phenology under natural and artificial photoperiods and the expression of PsFTL2 in a set of Scots pine populations from different latitudes. The expression of PsFTL2 increased in the needles preceding bud set and decreased during bud burst. In the northernmost population, even short night periods were efficient to trigger this expression, which also increased earlier under all photoperiodic regimes compared with the southern populations. Despite the different biology, with few limitations, the two conifers that diverged 140 million yr ago probably share an association of FTL2 with bud set, pointing to a common mechanism for the timing of growth cessation in conifers.

Within- and between-population comparisons suggest independently acting selection maintaining parallel clines in Scots pine (Pinus sylvestris).

Parallel clines in traits related to adaptation in a species can be due to independent selection on a pair of traits, or due to selection in one trait resulting in a parallel cline in a correlated trait. To distinguish between the mechanisms giving rise to parallel adaptive population divergence of multiple traits along an environmental gradient we need to study variation, correlations, and selective forces within individual populations along the gradient. In many tree species, budset timing (BST) forms a latitudinal cline, and parallel clinal variation is also found in other seedling traits, such as first-year height (FYH) and fall frost injury (FFI). In this study, we set up a common garden experiment with open pollinated progeny from natural populations of Scots pine (Pinus sylvestris), with one large sample from single population (500 families) and smaller samples from across a latitudinal gradient. BST, FYH and induced FFI were first measured in a greenhouse. The seedlings were then planted in the field, where survival and height were measured at the age of 9 years as fitness proxies. We compared between- and within-population variation and genetic correlations of these three seedling traits, and estimated selection gradients at the family level in our main population, taking into account the potential effects of seed weight. Between-population genetic correlations between seedling traits were high (0.76-0.95). Within-population genetic correlations in the main population were lower (0.14-0.35), as in other populations (0.10-0.39). Within population, extensive adaptive variation persists in the seedling traits, in line with rather weak selection gradients, yet maintaining the clines. Although our sampling does not cover the whole cline equally, the results suggest that the individual clines in these traits are maintained by largely independently acting selection, which results in fewer constraints in adaptation under changing climate.

An initial assessment of linkage disequilibrium (LD) in coffee trees: LD patterns in groups of Coffea canephora Pierre using microsatellite analysis.

BACKGROUND: A reciprocal recurrent selection program has been under way for the Coffea canephora coffee tree for approximately thirty years in the Ivory Coast. Association genetics would help to speed up this program by more rapidly selecting zones of interest in the genome. However, prior to any such studies, the linkage disequilibrium (LD) needs to be assessed between the markers on the genome. These data are essential for guiding association studies. RESULTS: This article describes the first results of an LD assessment in a coffee tree species. Guinean and Congolese breeding populations of C. canephora have been used for this work, with the goal of identifying ways of using these populations in association genetics. We identified changes in the LD along the genome within the different C. canephora diversity groups. In the different diversity groups studied, the LD was variable. Some diversity groups displayed disequilibria over long distances (up to 25 cM), whereas others had disequilibria not exceeding 1 cM. We also discovered a fine structure within the Guinean group. CONCLUSIONS: Given these results, association studies can be used within the species C. canephora. The coffee recurrent selection scheme being implemented in the Ivory Coast can thus be optimized. Lastly, our results could be used to improve C. arabica because one of its parents is closely related to C. canephora.

Genetic variability and QTL mapping of freezing tolerance and related traits in Medicago truncatula.

Freezing is a major environmental limitation to crop productivity for a number of species including legumes. We investigated the genetic determinism of freezing tolerance in the model legume Medicago truncatula Gaertn (M. truncatula). After having observed a large variation for freezing tolerance among 15 M. truncatula accessions, the progeny of a F6 recombinant inbred line population, derived from a cross between two accessions, was acclimated to low above-freezing temperatures and assessed for: (a) number of leaves (NOL), leaf area (LA), chlorophyll content index (CCI), shoot and root dry weights (SDW and RDW) at the end of the acclimation period and (b) visual freezing damage (FD) during the freezing treatment and 2 weeks after regrowth and foliar electrolyte leakage (EL) 2 weeks after regrowth. Consistent QTL positions with additive effects for FD were found on LG1, LG4 and LG6, the latter being the most explanatory (R (2) ≈ 40 %). QTL for NOL, QTL for EL, NOL and RDW, and QTL for EL and CCI colocalized with FD QTL on LG1, LG4 and LG6, respectively. Favorable alleles for these additive effects were brought by the same parent suggesting that this accession contributes to superior freezing tolerance by affecting plants' capacity to maintain growth at low above-freezing temperatures. No epistatic effects were found between FD QTL, but for each of the studied traits, 3-6 epistatic effects were detected between loci not detected directly as QTL. These results open the way to the assessment of syntenic relationships between QTL for frost tolerance in M. truncatula and cultivated legume species.

LiDAR Is Effective in Characterizing Vine Growth and Detecting Associated Genetic Loci.

The strong societal demand to reduce pesticide use and adaptation to climate change challenges the capacities of phenotyping new varieties in the vineyard. High-throughput phenotyping is a way to obtain meaningful and reliable information on hundreds of genotypes in a limited period. We evaluated traits related to growth in 209 genotypes from an interspecific grapevine biparental cross, between IJ119, a local genitor, and Divona, both in summer and in winter, using several methods: fresh pruning wood weight, exposed leaf area calculated from digital images, leaf chlorophyll concentration, and LiDAR-derived apparent volumes. Using high-density genetic information obtained by the genotyping by sequencing technology (GBS), we detected 6 regions of the grapevine genome [quantitative trait loci (QTL)] associated with the variations of the traits in the progeny. The detection of statistically significant QTLs, as well as correlations (R2) with traditional methods above 0.46, shows that LiDAR technology is effective in characterizing the growth features of the grapevine. Heritabilities calculated with LiDAR-derived total canopy and pruning wood volumes were high, above 0.66, and stable between growing seasons. These variables provided genetic models explaining up to 47% of the phenotypic variance, which were better than models obtained with the exposed leaf area estimated from images and the destructive pruning weight measurements. Our results highlight the relevance of LiDAR-derived traits for characterizing genetically induced differences in grapevine growth and open new perspectives for high-throughput phenotyping of grapevines in the vineyard.

The complete reference genome for grapevine (Vitis vinifera L.) genetics and breeding.

Grapevine is one of the most economically important crops worldwide. However, the previous versions of the grapevine reference genome tipically consist of thousands of fragments with missing centromeres and telomeres, limiting the accessibility of the repetitive sequences, the centromeric and telomeric regions, and the study of inheritance of important agronomic traits in these regions. Here, we assembled a telomere-to-telomere (T2T) gap-free reference genome for the cultivar PN40024 using PacBio HiFi long reads. The T2T reference genome (PN_T2T) is 69 Mb longer with 9018 more genes identified than the 12X.v0 version. We annotated 67% repetitive sequences, 19 centromeres and 36 telomeres, and incorporated gene annotations of previous versions into the PN_T2T assembly. We detected a total of 377 gene clusters, which showed associations with complex traits, such as aroma and disease resistance. Even though PN40024 derives from nine generations of selfing, we still found nine genomic hotspots of heterozygous sites associated with biological processes, such as the oxidation-reduction process and protein phosphorylation. The fully annotated complete reference genome therefore constitutes an important resource for grapevine genetic studies and breeding programs.

Adaptive potential of northernmost tree populations to climate change, with emphasis on Scots pine (Pinus sylvestris L.).

The adaptive potential of the northernmost Pinus sylvestris L. (and other northern tree) populations is considered by examining first the current patterns of quantitative genetic adaptive traits, which show high population differentiation and clines. We then consider the postglacial history of the populations using both paleobiological and genetic data. The current patterns of diversity at nuclear genes suggest that the traces of admixture are mostly visible in mitochondrial DNA variation patterns. There is little evidence of increased diversity due to admixture between an eastern and western colonization lineage, but no signal of reduced diversity (due to sequential bottlenecks) either. Quantitative trait variation in the north is not associated with the colonizing lineages. The current clines arose rapidly and may be based on standing genetic variation. The initial phenotypic response of Scots pine in the north is predicted to be increased survival and growth. The genetic responses are examined based on quantitative genetic predictions of sustained selection response and compared with earlier simulation results that have aimed at more ecological realism. The phenotypic responses of increased growth and survival reduce the opportunity for selection and delay the evolutionary responses. The lengthening of the thermal growing period also causes selection on the critical photoperiod in the different populations. Future studies should aim at including multiple ecological and genetic factors in evaluating potential responses.

instaGRAAL: chromosome-level quality scaffolding of genomes using a proximity ligation-based scaffolder.

Hi-C exploits contact frequencies between pairs of loci to bridge and order contigs during genome assembly, resulting in chromosome-level assemblies. Because few robust programs are available for this type of data, we developed instaGRAAL, a complete overhaul of the GRAAL program, which has adapted the latter to allow efficient assembly of large genomes. instaGRAAL features a number of improvements over GRAAL, including a modular correction approach that optionally integrates independent data. We validate the program using data for two brown algae, and human, to generate near-complete assemblies with minimal human intervention.

The genome of Ectocarpus subulatus - A highly stress-tolerant brown alga.

Brown algae are multicellular photosynthetic stramenopiles that colonize marine rocky shores worldwide. Ectocarpus sp. Ec32 has been established as a genomic model for brown algae. Here we present the genome and metabolic network of the closely related species, Ectocarpus subulatus Kützing, which is characterized by high abiotic stress tolerance. Since their separation, both strains show new traces of viral sequences and the activity of large retrotransposons, which may also be related to the expansion of a family of chlorophyll-binding proteins. Further features suspected to contribute to stress tolerance include an expanded family of heat shock proteins, the reduction of genes involved in the production of halogenated defence compounds, and the presence of fewer cell wall polysaccharide-modifying enzymes. Overall, E. subulatus has mainly lost members of gene families down-regulated in low salinities, and conserved those that were up-regulated in the same condition. However, 96% of genes that differed between the two examined Ectocarpus species, as well as all genes under positive selection, were found to encode proteins of unknown function. This underlines the uniqueness of brown algal stress tolerance mechanisms as well as the significance of establishing E. subulatus as a comparative model for future functional studies.

Association of sugar content QTL and PQL with physiological traits relevant to frost damage resistance in pea under field and controlled conditions.

To increase yield in pea (Pisum sativum L.), autumn sowing would be preferable. Hence, frost tolerance of pea became a major trait of interest for breeders. In order to better understand the cold acclimation in pea, Champagne a frost tolerant line and Terese, a frost sensitive line, and their recombinant inbred lines (RIL) were studied. RIL frost tolerance was evaluated by a frost damage scale under field as well as controlled conditions. A quantitative trait loci (QTL) approach was used to identify chromosomal regions linked to frost tolerance. The detected QTL explained from 6.5 to 46.5% of the phenotypic variance. Amongst them, those located on linkage groups 5 and 6 were consistent with over all experiments, in field as well as in controlled environments. In order to improve the understanding of the frost tolerance mechanisms, several cold acclimation key characters such as concentration of sugars, electrolyte leakage, osmotic pressure, and activity of RuBisCO were assessed. Some of these physiological QTL colocalised with QTL for frost damage, in particular two raffinose QTL on LG5 and LG6 and one RuBisCO activity QTL on LG6, explaining 8.8 to 27.0% of the phenotypic variance. In addition, protein quantitative loci were mapped; some of them colocalised with frost damage and physiological QTL on LG5 and LG6, explaining 16.0-43.6% of the phenotypic variance. Raffinose metabolism and RuBisCO activity and its effect on photosynthesis might play a major role in cold acclimation of pea.

Thermospermine Synthase (ACL5) and Diamine Oxidase (DAO) Expression Is Needed for Zygotic Embryogenesis and Vascular Development in Scots Pine.

Unlike in flowering plants, the detailed roles of the enzymes in the polyamine (PA) pathway in conifers are poorly known. We explored the sequence conservation of the PA biosynthetic genes and diamine oxidase (DAO) in conifers and flowering plants to reveal the potential functional diversification of the enzymes between the plant lineages. The expression of the genes showing different selective constraints was studied in Scots pine zygotic embryogenesis and early seedling development. We found that the arginine decarboxylase pathway is strongly preferred in putrescine production in the Scots pine as well as generally in conifers and that the reduced use of ornithine decarboxylase (ODC) has led to relaxed purifying selection in ODC genes. Thermospermine synthase (ACL5) genes evolve under strong purifying selection in conifers and the DAO gene is also highly conserved in pines. In developing Scots pine seeds, the expression of both ACL5 and DAO increased as embryogenesis proceeded. Strong ACL5 expression was present in the procambial cells of the embryo and in the megagametophyte cells destined to die via morphologically necrotic cell death. Thus, the high sequence conservation of ACL5 genes in conifers may indicate the necessity of ACL5 for both embryogenesis and vascular development. Moreover, the result suggests the involvement of ACL5 in morphologically necrotic cell death and supports the view of the genetic regulation of necrosis in Scots pine embryogenesis and in plant development. DAO transcripts were located close to the cell walls and between the walls of adjacent cells in Scots pine zygotic embryos and in the roots of young seedlings. We propose that DAO, in addition to the role in Put oxidation for providing H2O2 during the cell-wall structural processes, may also participate in cell-to-cell communication at the mRNA level. To conclude, our findings indicate that the PA pathway of Scots pines possesses several special functional characteristics which differ from those of flowering plants.

Convergent recruitment of TALE homeodomain life cycle regulators to direct sporophyte development in land plants and brown algae.

Three amino acid loop extension homeodomain transcription factors (TALE HD TFs) act as life cycle regulators in green algae and land plants. In mosses these regulators are required for the deployment of the sporophyte developmental program. We demonstrate that mutations in either of two TALE HD TF genes, OUROBOROS or SAMSARA, in the brown alga Ectocarpus result in conversion of the sporophyte generation into a gametophyte. The OUROBOROS and SAMSARA proteins heterodimerise in a similar manner to TALE HD TF life cycle regulators in the green lineage. These observations demonstrate that TALE-HD-TF-based life cycle regulation systems have an extremely ancient origin, and that these systems have been independently recruited to regulate sporophyte developmental programs in at least two different complex multicellular eukaryotic supergroups, Archaeplastida and Chromalveolata.

Genetic Diversity in the UV Sex Chromosomes of the Brown Alga Ectocarpus.

Three types of sex chromosome system exist in nature: diploid XY and ZW systems and haploid UV systems. For many years, research has focused exclusively on XY and ZW systems, leaving UV chromosomes and haploid sex determination largely neglected. Here, we perform a detailed analysis of DNA sequence neutral diversity levels across the U and V sex chromosomes of the model brown alga Ectocarpus using a large population dataset. We show that the U and V non-recombining regions of the sex chromosomes (SDR) exhibit about half as much neutral diversity as the autosomes. This difference is consistent with the reduced effective population size of these regions compared with the rest of the genome, suggesting that the influence of additional factors such as background selection or selective sweeps is minimal. The pseudoautosomal region (PAR) of this UV system, in contrast, exhibited surprisingly high neutral diversity and there were several indications that genes in this region may be under balancing selection. The PAR of Ectocarpus is known to exhibit unusual genomic features and our results lay the foundation for further work aimed at understanding whether, and to what extent, these structural features underlie the high level of genetic diversity. Overall, this study fills a gap between available information on genetic diversity in XY/ZW systems and UV systems and significantly contributes to advancing our knowledge of the evolution of UV sex chromosomes.

High-density genetic map and identification of QTLs for responses to temperature and salinity stresses in the model brown alga Ectocarpus.

Deciphering the genetic architecture of adaptation of brown algae to environmental stresses such as temperature and salinity is of evolutionary as well as of practical interest. The filamentous brown alga Ectocarpus sp. is a model for the brown algae and its genome has been sequenced. As sessile organisms, brown algae need to be capable of resisting the various abiotic stressors that act in the intertidal zone (e.g. osmotic pressure, temperature, salinity, UV radiation) and previous studies have shown that an important proportion of the expressed genes is regulated in response to hyposaline, hypersaline or oxidative stress conditions. Using the double digest RAD sequencing method, we constructed a dense genetic map with 3,588 SNP markers and identified 39 QTLs for growth-related traits and their plasticity under different temperature and salinity conditions (tolerance to high temperature and low salinity). GO enrichment tests within QTL intervals highlighted membrane transport processes such as ion transporters. Our study represents a significant step towards deciphering the genetic basis of adaptation of Ectocarpus sp. to stress conditions and provides a substantial resource to the increasing list of tools generated for the species.