Pine Genome Research Articles

Transcription and splicing variations of SR genes accompany with genome-wide accumulation of long-introns in pine

Most of mRNAs in Eukaryote were matured after the removal of introns in their pre-mRNA transcripts. Serine/arginine-rich (SR) proteins are a group of splicing regulators regulating the splicing processes globally. Expressions of SR proteins themselves were extensively regulated, at both transcription and splicing levels, under different environmental conditions, specially heat stress conditions. The pine genome is characterized by super-long and easily methylated introns in a large number of genes that derived from the extensive accumulation of transposons (TEs). Here, we identified and analyzed the phylogenetic characteristics of 24 SR proteins and their encoding genes from the pine genome. Then we explored transcription and pre-mRNA splicing expression patterns of SR genes in P. massoniana seedlings under normal and heat stress temperature conditions. Our results showed that the transcription patterns of SR genes in pine exhibited significant changes compared to other plant species, and these changes were not strictly correlated with the intron length and DNA methylation intensity of the SR genes. Interestingly, none of the long introns of SR genes underwent alternative splicing (AS) in our experiment. Furthermore, the intensity of AS regulation may be related to the potential DNA methylation intensity of SR genes. Taken together, this study explores for the first time the characteristics of significant variations in the transcription and splicing patterns of SR proteins in a plant species with an over-accumulation of super-long introns.

Stone Pine (Pinus pinea L.) High-Added-Value Genetics: An Overview.

Stone pine (Pinus pinea L.) has received limited attention in terms of genetic research. However, genomic techniques hold promise for decoding the stone pine genome and contributing to developing a more resilient bioeconomy. Retrotransposon and specific genetic markers are effective tools for determining population-specific genomic diversity. Studies on the transcriptome and proteome have identified differentially expressed genes PAS1, CLV1, ATAF1, and ACBF involved in shoot bud formation. The stone pine proteome shows variation among populations and shows the industrial potential of the enzyme pinosylvin. Microsatellite studies have revealed low levels of polymorphism and a unique genetic diversity in stone pine, which may contribute to its environmental adaptation. Transcriptomic and proteomic analyses uncover the genetic and molecular responses of stone pine to fungal infections and nematode infestations, elucidating the defense activation, gene regulation, and the potential role of terpenes in pathogen resistance. Transcriptomics associated with carbohydrate metabolism, dehydrins, and transcription factors show promise as targets for improving stone pine's drought stress response and water retention capabilities. Stone pine presents itself as an important model tree for studying climate change adaptation due to its characteristics. While knowledge gaps exist, stone pine's genetic resources hold significant potential, and ongoing advancements in techniques offer prospects for future exploration.

The methylation landscape of giga-genome and the epigenetic timer of age in Chinese pine

Epigenetics has been revealed to play a crucial role in the long-term memory in plants. However, little is known about whether the epigenetic modifications occur with age progressively in conifers. Here, we present the single-base resolution DNA methylation landscapes of the 25-gigabase Chinese pine (Pinus tabuliformis) genome at different ages. The result shows that DNA methylation is closely coupled with the regulation of gene transcription. The age-dependent methylation profile with a linearly increasing trend is the most significant pattern of DMRs between ages. Two segments at the five-prime end of the first ultra-long intron in DAL1, a conservative age biomarker in conifers, shows a gradual decline of CHG methylation as the age increased, which is highly correlated with its expression profile. Similar high correlation is also observed in nine other age marker genes. Our results suggest that DNA methylation serves as an important epigenetic signature of developmental age in conifers.

The Chinese pine genome and methylome unveil key features of conifer evolution

Conifers dominate the world's forest ecosystems and are the most widely planted tree species. Their giant and complex genomes present great challenges for assembling a complete reference genome for evolutionary and genomic studies. We present a 25.4-Gb chromosome-level assembly of Chinese pine (Pinus tabuliformis) and revealed that its genome size is mostly attributable to huge intergenic regions and long introns with high transposable element (TE) content. Large genes with long introns exhibited higher expressions levels. Despite a lack of recent whole-genome duplication, 91.2% of genes were duplicated through dispersed duplication, and expanded gene families are mainly related to stress responses, which may underpin conifers' adaptation, particularly in cold and/or arid conditions. The reproductive regulation network is distinct compared with angiosperms. Slow removal of TEs with high-level methylation may have contributed to genomic expansion. This study provides insights into conifer evolution and resources for advancing research on conifer adaptation and development.

NLR diversity and candidate fusiform rust resistance genes in loblolly pine

Resistance to fusiform rust disease in loblolly pine (Pinus taeda) is a classic gene-for-gene system. Early resistance gene mapping in the P. taeda family 10-5 identified RAPD markers for a major fusiform rust resistance gene, Fr1. More recently, single nucleotide polymorphism (SNP) markers associated with resistance were mapped to a full-length gene model in the loblolly pine genome encoding for a nucleotide-binding site leucine-rich repeat (NLR) protein. NLR genes are one of the most abundant gene families in plant genomes and are involved in effector-triggered immunity. Inter- and intraspecies studies of NLR gene diversity and expression have resulted in improved disease resistance. To characterize NLR gene diversity and discover potential resistance genes, we assembled de novo transcriptomes from 92 loblolly genotypes from across the natural range of the species. In these transcriptomes, we identified novel NLR transcripts that are not present in the loblolly pine reference genome and found significant geographic diversity of NLR genes providing evidence of gene family evolution. We designed capture probes for these NLRs to identify and map SNPs that stably cosegregate with resistance to the SC20-21 isolate of Cronartium quercuum f.sp. fusiforme (Cqf) in half-sib progeny of the 10-5 family. We identified 10 SNPs and 2 quantitative trait loci associated with resistance to SC20-21 Cqf. The geographic diversity of NLR genes provides evidence of NLR gene family evolution in loblolly pine. The SNPs associated with rust resistance provide a resource to enhance breeding and deployment of resistant pine seedlings.

Fine dissection of limber pine resistance to Cronartium ribicola using targeted sequencing of the NLR family

BackgroundProteins with nucleotide binding site (NBS) and leucine-rich repeat (LRR) domains (NLR) make up one of most important resistance (R) families for plants to resist attacks from various pathogens and pests. The available transcriptomes of limber pine (Pinus flexilis) allow us to characterize NLR genes and related resistance gene analogs (RGAs) in host resistance against Cronartium ribicola, the causal fungal pathogen of white pine blister rust (WPBR) on five-needle pines throughout the world. We previously mapped a limber pine major gene locus (Cr4) that confers complete resistance to C. ribicola on the Pinus consensus linkage group 8 (LG-8). However, genetic distribution of NLR genes as well as their divergence between resistant and susceptible alleles are still unknown.ResultsTo identify NLR genes at the Cr4 locus, the present study re-sequenced a total of 480 RGAs using targeted sequencing in a Cr4-segregated seed family. Following a call of single nucleotide polymorphisms (SNPs) and genetic mapping, a total of 541 SNPs from 155 genes were mapped across 12 LGs. Three putative NLR genes were newly mapped in the Cr4 region, including one that co-segregated with Cr4. The tight linkage of NLRs with Cr4-controlled phenotypes was further confirmed by bulked segregation analysis (BSA) using extreme-phenotype genome-wide association study (XP-GWAS) for significance test. Local tandem duplication in the Cr4 region was further supported by syntenic analysis using the sugar pine genome sequence. Significant gene divergences have been observed in the NLR family, revealing that diversifying selection pressures are relatively higher in local duplicated genes. Most genes showed similar expression patterns at low levels, but some were affected by genetic background related to disease resistance. Evidence from fine genetic dissection, evolutionary analysis, and expression profiling suggests that two NLR genes are the most promising candidates for Cr4 against WPBR.ConclusionThis study provides fundamental insights into genetic architecture of the Cr4 locus as well as a set of NLR variants for marker-assisted selection in limber pine breeding. Novel NLR genes were identified at the Cr4 locus and the Cr4 candidates will aid deployment of this R gene in combination with other major/minor genes in the limber pine breeding program.

Structural and Functional Characteristics of Two Molecular Variants of the Nitrogen Sensor PII in Maritime Pine.

High levels of nitrogen are stored as arginine during the last stages of seed formation in maritime pine (Pinus pinaster Aiton). The protein sensor PII regulates the feedback inhibition of arginine biosynthesis through interaction with the key enzyme N-acetylglutamate kinase (NAGK). In this study, the structural and functional characteristics of PII have been investigated in maritime pine to get insights into the regulation of arginine metabolism. Two different forms of PII have been identified, PpPIIa and PpPIIb, which differ in their amino acid sequence and most likely correspond to splicing variants of a single gene in the pine genome. Two PII variants are also present in other pine species but not in other conifers such as spruces. PpPIIa and PpPIIb are trimeric proteins for which structural modeling predicts similar tridimensional protein core structures. Both are located in the chloroplast, where the PII-target enzyme PpNAGK is also found. PpPIIa, PpPIIb, and PpNAGK have been recombinantly produced to investigate the formation of NAGK-PII complexes. The interaction of PpPIIa/PpPIIb and PpNAGK may be enhanced by glutamine and contribute to relieve the feedback inhibition of PpNAGK by arginine. Expression analysis of PpPII genes revealed that PpIIa transcripts were predominant during embryogenesis and germination. The potential roles of PpPIIa and PpPIIb in the regulation of arginine metabolism of maritime pine are discussed.

Involvement of VNS NAC-domain transcription factors in tracheid formation in Pinus taeda.

Vascular plants have two types of water-conducting cells, xylem vessel cells (in angiosperms) and tracheid cells (in ferns and gymnosperms). These cells are commonly characterized by secondary cell wall (SCW) formation and programmed cell death (PCD), which increase the efficiency of water conduction. The differentiation of xylem vessel cells is regulated by a set of NAC (NAM, ATAF1/2 and CUC2) transcription factors, called the VASCULAR-RELATED NAC-DOMAIN (VND) family, in Arabidopsis thaliana Linne. The VNDs regulate the transcriptional induction of genes required for SCW formation and PCD. However, information on the transcriptional regulation of tracheid cell differentiation is still limited. Here, we performed functional analysis of loblolly pine (Pinus taeda Linne) VND homologs (PtaVNS, for VND, NST/SND, SMB-related protein). We identified five PtaVNS genes in the loblolly pine genome, and four of these PtaVNS genes were highly expressed in tissues with tracheid cells, such as shoot apices and developing xylem. Transient overexpression of PtaVNS genes induced xylem vessel cell-like patterning of SCW deposition in tobacco (Nicotiana benthamiana Domin) leaves, and up-regulated the promoter activities of loblolly pine genes homologous to SCW-related MYB transcription factor genes and cellulose synthase genes, as well as to cysteine protease genes for PCD. Collectively, our data indicated that PtaVNS proteins possess transcriptional activity to induce the molecular programs required for tracheid formation, i.e., SCW formation and PCD. Moreover, these findings suggest that the VNS-MYB-based transcriptional network regulating water-conducting cell differentiation in angiosperm and moss plants is conserved in gymnosperms.

Association mapping of ectomycorrhizal traits in loblolly pine (Pinus taeda L.).

To understand how diverse mutualisms coevolve and how species adapt to complex environments, a description of the underlying genetic basis of the traits involved must be provided. For example, in diverse coevolving mutualisms, such as the interaction of host plants with a suite of symbiotic mycorrhizal fungi, a key question is whether host plants can coevolve independently with multiple species of symbionts, which depends on whether those interactions are governed independently by separate genes or pleiotropically by shared genes. To provide insight into this question, we employed an association mapping approach in a clonally replicated field experiment of loblolly pine (Pinus taeda L.) to identify genetic components of host traits governing ectomycorrhizal (EM) symbioses (mycorrhizal traits). The relative abundances of different EM fungi as well as the total number of root tips per cm root colonized by EM fungi were analyzed as separate mycorrhizal traits of loblolly pine. Single-nucleotide polymorphisms (SNPs) within candidate genes of loblolly pine were associated with loblolly pine mycorrhizal traits, mapped to the loblolly pine genome, and their putative protein function obtained when available. The results support the hypothesis that ectomycorrhiza formation is governed by host genes of large effect that apparently have independent influences on host interactions with different symbiont species.

Resequencing of massive pine genomes helps to unlock the genetic underpinning of quantitative traits in conifer trees.

This article is a Commentary on De La Torre et al., 221: 1789–1801.

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.

Glutamate synthases from conifers: gene structure and phylogenetic studies

BackgroundPlants synthesize glutamate from ammonium by the combined activity of the enzymes glutamine synthetase (GS) and glutamate synthase (GOGAT) through the glutamate synthase cycle. In plants, there are two forms of glutamate synthases that differ in their electron donors, NADH-GOGAT (EC 1.4.1.14) and Fd-GOGAT (EC 1.4.7.1), which have differential roles either in primary ammonia assimilation or in the reassimilation of ammonium from different catabolic processes. Glutamate synthases are complex iron-sulfur flavoproteins containing functional domains involved in the control and coordination of their catalytic activities in annual plants. In conifers, partial cDNA sequences for GOGATs have been isolated and used for gene expression studies. However, knowledge of the gene structure and of phylogenetic relationships with other plant enzymes is quite scant.ResultsTechnological advances in conifer megagenomes sequencing have made it possible to obtain full-length cDNA sequences encoding Fd- and NADH-GOGAT from maritime pine, as well as BAC clones containing sequences for NADH-GOGAT and Fd-GOGAT genes. In the current study, we studied the genomic organization of pine GOGAT genes, the size of their exons/introns, copy numbers in the pine genome and relationships with other plant genes. Phylogenetic analysis was performed, and the degree of preservation and dissimilarity of key domains for the catalytic activities of these enzymes in different taxa were determined.ConclusionsFd- and NADH-GOGAT are encoded by single-copy genes in the maritime pine genome. The Fd-GOGAT gene is extremely large spanning more than 330 kb and the presence of very long introns highlights the important contribution of LTR retrotransposons to the gene size in conifers. In contrast, the structure of the NADH-GOGAT gene is similar to the orthologous genes in angiosperms. Our phylogenetic analysis indicates that these two genes had different origins during plant evolution. The results provide new insights into the structure and molecular evolution of these essential genes.

Assessment of gene copy number variation of Scots pine thaumatin-like protein gene using real-time PCR based methods

The importance and impact of gene copy number variations (CNVs) as a source of polymorphism in the human and other genomes is being increasingly recognized. Less information is available about CNVs in forest tree species, mainly due to the relative lack of genomic resources. In this study, several methods—quantitative polymerase chain reaction, comparative high-resolution melting curve analysis (C-HRM), and digital polymerase chain reaction (dPCR)—were used to investigate CNV of the Scots pine thaumatin-like protein gene (PsTLP). The obtained results were supported by transcriptome analysis of a single Pinus sylvestris individual and publically available pine genome sequences. Although estimations of gene copy number (CN) varied, depending on the region of the PsTLP gene investigated and the endogenous control utilized, our results revealed the existence of copy number variations of the PsTLP gene between Scots pine individuals. Of 23 individuals analyzed, two had an increased calculated relative CN regardless of the analyzed gene region and endogenous control used, while several samples had increased copy numbers of regions of the PsTLP gene. C-HRM results were highly correlated with qPCR data (R 2 TLP3′ = 0.88; R 2 TLPc = 0.92), but interpretation of gene CN from C-HRM results proved to be difficult. The results from selected samples analyzed by digital PCR also were highly correlated with qPCR results (R 2 = 0.90).

Development of microsatellite genetic markers in Siberian stone pine (Pinus sibirica Du Tour) based on the de novo whole genome sequencing

Siberian stone pine, Pinus sibirica Du Tour is one of the most economically and environmentally important forest-forming species of conifers in Russia. To study these forests a large number of highly polymorphic molecular genetic markers, such as microsatellite loci, are required. Prior to the new high-throughput next generation sequencing (NGS) methods, discovery of microsatellite loci and development of micro-satellite markers were very time consuming and laborious. The recently developed draft assembly of the Siberian stone pine genome, sequenced using the NGS methods, allowed us to identify a large number of microsatellite loci in the Siberian stone pine genome and to develop species-specific PCR primers for amplification and genotyping of 70 microsatellite loci. The primers were designed using contigs containing short simple sequence tandem repeats from the Siberian stone pine whole genome draft assembly. Based on the testing of primers for 70 microsatellite loci with tri-, tetra- or pentanucleotide repeats, 18 most promising, reliable and polymorphic loci were selected that can be used further as molecular genetic markers in population genetic studies of Siberian stone pine.

Complete Chloroplast Genome of the Wollemi Pine (Wollemia nobilis): Structure and Evolution.

The Wollemi pine (Wollemia nobilis) is a rare Southern conifer with striking morphological similarity to fossil pines. A small population of W. nobilis was discovered in 1994 in a remote canyon system in the Wollemi National Park (near Sydney, Australia). This population contains fewer than 100 individuals and is critically endangered. Previous genetic studies of the Wollemi pine have investigated its evolutionary relationship with other pines in the family Araucariaceae, and have suggested that the Wollemi pine genome contains little or no variation. However, these studies were performed prior to the widespread use of genome sequencing, and their conclusions were based on a limited fraction of the Wollemi pine genome. In this study, we address this problem by determining the entire sequence of the W. nobilis chloroplast genome. A detailed analysis of the structure of the genome is presented, and the evolution of the genome is inferred by comparison with the chloroplast sequences of other members of the Araucariaceae and the related family Podocarpaceae. Pairwise alignments of whole genome sequences, and the presence of unique pseudogenes, gene duplications and insertions in W. nobilis and Araucariaceae, indicate that the W. nobilis chloroplast genome is most similar to that of its sister taxon Agathis. However, the W. nobilis genome contains an unusually high number of repetitive sequences, and these could be used in future studies to investigate and conserve any remnant genetic diversity in the Wollemi pine.

The genetic architecture of local adaptation I: the genomic landscape of foxtail pine (Pinus balfouriana Grev. & Balf.) as revealed from a high-density linkage map

Explaining the origin and evolutionary dynamics of the genetic architecture of adaptation is a major research goal of evolutionary genetics. Despite controversy surrounding success of the attempts to accomplish this goal, a full understanding of adaptive genetic variation necessitates knowledge about the genomic location and patterns of dispersion for the genetic components affecting fitness-related phenotypic traits. Even with advances in next-generation sequencing technologies, the production of full genome sequences for non-model species is often cost-prohibitive, especially for tree species such as pines where genome size often exceeds 20 to 30 Gbp. We address this need by constructing a dense linkage map for foxtail pine (Pinus balfouriana Grev. & Balf.), with the ultimate goal of uncovering and explaining the origin and evolutionary dynamics of adaptive genetic variation in natural populations of this forest tree species. We utilized megagametophyte arrays (n = 76–95 megagametophytes/tree) from four maternal trees in combination with double digest restriction-site associated DNA sequencing (ddRADseq) to produce a consensus linkage map covering 98.58 % of the foxtail pine genome, which was estimated to be 1276 cM in length (95 % CI, 1174 to 1378 cM). A novel bioinformatic approach using iterative rounds of marker ordering and imputation was employed to produce single-tree linkage maps (507–17,066 contigs/map; lengths, 1037.40–1572.80 cM). These linkage maps were collinear across maternal trees, with highly correlated marker orderings (Spearman’s ρ>0.95). A consensus linkage map derived from these single-tree linkage maps contained 12 linkage groups along which 20,655 contigs were non-randomly distributed across 901 unique positions (n=23 contigs/position), with an average spacing of 1.34 cM between adjacent positions. Of the 20,655 contigs positioned on the consensus linkage map, 5627 had enough sequence similarity to contigs contained within the most recent build of the loblolly pine (Pinus taeda L.) genome to identify them as putative homologues containing both genic and non-genic loci. Importantly, all 901 unique positions on the consensus linkage map had at least one contig with putative homology to loblolly pine. When combined with the other biological signals that predominate in our data (e.g., correlations of recombination fractions across single trees), we show that dense linkage maps for non-model forest tree species can be efficiently constructed using next-generation sequencing technologies. We subsequently discuss the usefulness of these maps as community-wide resources and as tools with which to test hypotheses about the genetic architecture of local adaptation.

Optimization of the genotyping-by-sequencing strategy for population genomic analysis in conifers.

Flexibility and low cost make genotyping-by-sequencing (GBS) an ideal tool for population genomic studies of nonmodel species. However, to utilize the potential of the method fully, many parameters affecting library quality and single nucleotide polymorphism (SNP) discovery require optimization, especially for conifer genomes with a high repetitive DNA content. In this study, we explored strategies for effective GBS analysis in pine species. We constructed GBS libraries using HpaII, PstI and EcoRI-MseI digestions with different multiplexing levels and examined the effect of restriction enzymes on library complexity and the impact of sequencing depth and size selection of restriction fragments on sequence coverage bias. We tested and compared UNEAK, Stacks and GATK pipelines for the GBS data, and then developed a reference-free SNP calling strategy for haploid pine genomes. Our GBS procedure proved to be effective in SNP discovery, producing 7000-11 000 and 14 751 SNPs within and among three pine species, respectively, from a PstI library. This investigation provides guidance for the design and analysis of GBS experiments, particularly for organisms for which genomic information is lacking.

Stress-induced transcriptional activation of retrotransposon-like sequences in the Scots pine (Pinus sylvestris L.) genome

Transposition of mobile elements has been implicated in genome instability, rearrangements and therefore also adaptation to changing environmental conditions. Transposons could influence gene activity directly by transposition inside or close to coding regions by their disruption or by addition of regulative sequences. Further, class I transposable elements, which are the most abundant in plant genomes, utilize a RNA intermediate in their life cycle, therefore retrotransposons could act by producing non-coding RNAs that could affect other transcripts by RNA interference. Transposition activity is suppressed by chromatin modifications, and both classes of transposons have been shown to be activated in plants under various stress conditions and developmental stages. Using a nonspecific amplification approach, we demonstrate the differential transcriptional activation of sequences with homology to transposable elements and other associated sequences in the complex genome of Scots pine (Pinus sylvestris L.) after exposure to heat stress, infestation with pine woolly aphids, and salicylic acid and abscisic acid treatment. Sequences with homology to several retrotransposon classes and families were identified, as well as several chimeric transcript types. Some of them represent chloroplast sequence insertions into the pine nuclear genome and these sequences are highly represented in EST databases of a wide range of species. In this study, we identified several retrotransposon classes and families with differing levels of similarity with known transposable elements from other plant species, and which are differentially expressed under various stress conditions in Scots pine.

Selection of haploid cell lines from megagametophyte cultures of maritime pine as a DNA source for massive sequencing of the species

The potential of haploid tissues for genetic studies in conifers is hampered by the lack of abundant and homogeneous plant material suitable for DNA isolation. In this work we have determined factors promoting haploid callus induction and proliferation from megagametophytes of Oria 6, a genotype of Pinus pinaster Aiton (maritime pine) from the natural population Sierra de Oria (Almeria, Spain), selected based on its response to extreme drought conditions. The generation of haploid callus was restricted to megagametophytes isolated from light-brown cones with no dehydrated seeds collected in September. Culture medium composition did not significantly affect callus induction, but a modified Murashige and Skoog medium with 2,4-dichlorophenoxyacetic acid and 6-benzyladenine favored further multiplication. The ploidy status of the callus lines was determined by flow cytometry and seven polymorphic microsatellites. A total of sixteen haploid callus lines were established and one of these is being used as a source of DNA for massive sequencing of maritime pine genome.

Sequencing and Assembly of the 22-Gb Loblolly Pine Genome

Conifers are the predominant gymnosperm. The size and complexity of their genomes has presented formidable technical challenges for whole-genome shotgun sequencing and assembly. We employed novel strategies that allowed us to determine the loblolly pine (Pinus taeda) reference genome sequence, the largest genome assembled to date. Most of the sequence data were derived from whole-genome shotgun sequencing of a single megagametophyte, the haploid tissue of a single pine seed. Although that constrained the quantity of available DNA, the resulting haploid sequence data were well-suited for assembly. The haploid sequence was augmented with multiple linking long-fragment mate pair libraries from the parental diploid DNA. For the longest fragments, we used novel fosmid DiTag libraries. Sequences from the linking libraries that did not match the megagametophyte were identified and removed. Assembly of the sequence data were aided by condensing the enormous number of paired-end reads into a much smaller set of longer “super-reads,” rendering subsequent assembly with an overlap-based assembly algorithm computationally feasible. To further improve the contiguity and biological utility of the genome sequence, additional scaffolding methods utilizing independent genome and transcriptome assemblies were implemented. The combination of these strategies resulted in a draft genome sequence of 20.15 billion bases, with an N50 scaffold size of 66.9 kbp.