Genome-wide association analysis of Monilinia fructicola lesion in a collection of Spanish peach landraces.

Brown rot, caused by the Monilinia spp., is the disease that causes the greatest losses in stone fruit worldwide. Currently, M. fructicola has become the dominant species in the main peach production area in Spain. The fruit cuticle is the first barrier of protection against external aggressions and may have a key role in the susceptibility to brown rot. However, information on the role of skin fruit on the resistance to brown rot in peach is scarce. Previous genetic analyses in peach have demonstrated that brown rot resistance is a complex and quantitative trait in which different fruit parts and resistance mechanisms are involved. To search for genomic areas involved in the control of the cultivar susceptibility to brown rot and to elucidate the role of fruit skin against this infection, we have studied, for two consecutive seasons (2019 and 2020), the fruit susceptibility to M. fructicola, together with fruit cuticle thickness (CT) and density (CD), in a collection of 80 Spanish and 5 foreign peach cultivars from the National Peach Collection at CITA (Zaragoza, Spain). Brown rot incidence, lesion diameter, and severity index were calculated after 5 days of inoculation on non-wounded fruit. The peach collection has also been genotyped using the new peach SNP chip (9 + 9K). Genotypic and phenotypic data have been used to perform a genome-wide association analysis (GWAS). Phenotyping has shown a wide variability on the brown rot susceptibility within the Spanish germplasm as well as on CD and CT. The GWAS results have identified several significant SNPs associated with disease severity index (DSI), CD, and CT, five of which were considered as reliable SNP-trait associations. A wide protein network analysis, using 127 genes within the regions of the reliable SNPs and previously identified candidate genes (169) associated with Monilinia spp. resistance, highlighted several genes involved in classical hypersensitive response (HR), genes related to wax layers as ceramidases and lignin precursors catalyzers, and a possible role of autophagy during brown rot infection. This work adds relevant information on the complexity resistance mechanisms to brown rot infection in peach fruits and the genetics behind them.

Orchard Management and Incorporation of Biochemical and Molecular Strategies for Improving Drought Tolerance in Fruit Tree Crops.

Water scarcity is one of the greatest concerns for agronomy worldwide. In recent years, many water resources have been depleted due to multiple factors, especially mismanagement. Water resource shortages lead to cropland expansion, which likely influences climate change and affects global agriculture, especially horticultural crops. Fruit yield is the final aim in commercial orchards; however, drought can slow tree growth and/or decrease fruit yield and quality. It is therefore necessary to find approaches to solve this problem. The main objective of this review is to discuss the most recent horticultural, biochemical, and molecular strategies adopted to improve the response of temperate fruit crops to water stress. We also address the viability of cultivating fruit trees in dry areas and provide precise protection methods for planting fruit trees in arid lands. We review the main factors involved in planting fruit trees in dry areas, including plant material selection, regulated deficit irrigation (DI) strategies, rainwater harvesting (RWH), and anti-water stress materials. We also provide a detailed analysis of the molecular strategies developed to combat drought, such as Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR) through gene overexpression or gene silencing. Finally, we look at the molecular mechanisms associated with the contribution of the microbiome to improving plant responses to drought.

Development and Evaluation of an AxiomTM 60K SNP Array for Almond (Prunus dulcis).

A high-density single nucleotide polymorphism (SNP) array is essential to enable faster progress in plant breeding for new cultivar development. In this regard, we have developed an Axiom 60K almond SNP array by resequencing 81 almond accessions. For the validation of the array, a set of 210 accessions were genotyped and 82.8% of the SNPs were classified in the best recommended SNPs. The rate of missing data was between 0.4% and 2.7% for the almond accessions and less than 15.5% for the few peach and wild accessions, suggesting that this array can be used for peach and interspecific peach × almond genetic studies. The values of the two SNPs linked to the RMja (nematode resistance) and SK (bitterness) genes were consistent. We also genotyped 49 hybrids from an almond F2 progeny and could build a genetic map with a set of 1159 SNPs. Error rates, less than 1%, were evaluated by comparing replicates and by detection of departures from Mendelian inheritance in the F2 progeny. This almond array is commercially available and should be a cost-effective genotyping tool useful in the search for new genes and quantitative traits loci (QTL) involved in the control of agronomic traits.

Spontaneous, Artificial, and Genome Editing-Mediated Mutations in Prunus.

Mutation is a source of genetic diversity widely used in breeding programs for the acquisition of agronomically interesting characters in commercial varieties of the Prunus species, as well as in the rest of crop species. Mutation can occur in nature at a very low frequency or can be induced artificially. Spontaneous or bud sport mutations in somatic cells can be vegetatively propagated to get an individual with the mutant phenotype. Unlike animals, plants have unlimited growth and totipotent cells that let somatic mutations to be transmitted to the progeny. On the other hand, in vitro tissue culture makes it possible to induce mutation in plant material and perform large screenings for mutant's selection and cleaning of chimeras. Finally, targeted mutagenesis has been boosted by the application of CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)/Cas9 and Transcription activator-like effector nuclease (TALEN) editing technologies. Over the last few decades, environmental stressors such as global warming have been threatening the supply of global demand for food based on population growth in the near future. For this purpose, the release of new varieties adapted to such changes is a requisite, and selected or generated Prunus mutants by properly regulated mechanisms could be helpful to this task. In this work, we reviewed the most relevant mutations for breeding traits in Prunus species such as flowering time, self-compatibility, fruit quality, and disease tolerance, including new molecular perspectives in the present postgenomic era including CRISPR/Cas9 and TALEN editing technologies.

Whole Transcriptome Analyses of Apricots and Japanese Plum Fruits after 1-MCP (Ethylene-Inhibitor) and Ethrel (Ethylene-Precursor) Treatments Reveal New Insights into the Physiology of the Ripening Process.

The physiology of Prunus fruit ripening is a complex and not completely understood process. To improve this knowledge, postharvest behavior during the shelf-life period at the transcriptomic level has been studied using high-throughput sequencing analysis (RNA-Seq). Monitoring of fruits has been analyzed after different ethylene regulator treatments, including 1-MCP (ethylene-inhibitor) and Ethrel (ethylene-precursor) in two contrasting selected apricot (Prunus armeniaca L.) and Japanese plum (P. salicina L.) cultivars, 'Goldrich' and 'Santa Rosa'. KEEG and protein-protein interaction network analysis unveiled that the most significant metabolic pathways involved in the ripening process were photosynthesis and plant hormone signal transduction. In addition, previously discovered genes linked to fruit ripening, such as pectinesterase or auxin-responsive protein, have been confirmed as the main genes involved in this process. Genes encoding pectinesterase in the pentose and glucuronate interconversions pathway were the most overexpressed in both species, being upregulated by Ethrel. On the other hand, auxin-responsive protein IAA and aquaporin PIP were both upregulated by 1-MCP in 'Goldrich' and 'Santa Rosa', respectively. Results also showed the upregulation of chitinase and glutaredoxin 3 after Ethrel treatment in 'Goldrich' and 'Santa Rosa', respectively, while photosystem I subunit V psaG (photosynthesis) was upregulated after 1-MCP in both species. Furthermore, the overexpression of genes encoding GDP-L-galactose and ferredoxin in the ascorbate and aldarate metabolism and photosynthesis pathways caused by 1-MCP favored antioxidant activity and therefore slowed down the fruit senescence process.

Genetic Diversity and Genome-Wide Association Study of Morphological and Quality Traits in Peach Using Two Spanish Peach Germplasm Collections.

Peach [Prunus persica (L.) Batsch] is one of the most important stone fruits species in world production. Spanish peach production is currently the second largest in the world and the available cultivars in Spain includes a great source of genetic diversity with variability in fruit quality traits and postharvest disorders tolerance. In order to explore the genetic diversity and single nucleotide polymorphism (SNP)-trait associations in the Spanish germplasm, the new peach 18K SNP v2 array was used to genotype 287 accessions belonging to the two National Peach Germplasm Collections placed at the Agrifood Research and Technology Centre of Aragon (CITA) and at the Experimental Station of Aula Dei (EEAD)-CSIC. The high density of the new SNP array allowed the identification of 30 groups of synonymies, which had not been identified before using low-density markers. In addition, a possible large-scale molecular event in 'Starcrest', a sport of 'Springcrest', was detected showing a possible chromosome replacement of a 13.5 Mb region. Previous suggestions about Spanish diversification regions agreed with our genetic diversity and linkage disequilibrium (LD) decay results using high-density markers. A genome-wide association study (GWAS) detected 34 significant SNP-trait association with the type of leaf glands (TLG), fruit hairiness (FH), and flesh texture (FT). The impact of the significant SNPs was studied with SnpEff. Candidate genes encode several important family proteins involved in trichome formation and powdery mildew resistance (linked to TLG in peach). The genetic distance among cultivars obtained, together with SNP-trait associations found, provide new knowledge for marker-assisted selection and crossing approaches in peach breeding programmes.

sRNA Analysis Evidenced the Involvement of Different Plant Viruses in the Activation of RNA Silencing-Related Genes and the Defensive Response Against Plum pox virus of 'GF305' Peach Grafted with 'Garrigues' Almond.

Plum pox virus (PPV) causes sharka disease in Prunus trees. Peach (P. persica) trees are severely affected by PPV, and no definitive source of genetic resistance has been identified. However, previous results showed that PPV-resistant 'Garrigues' almond (P. dulcis) was able to transfer its resistance to 'GF305' peach through grafting, reducing symptoms and viral load in PPV-infected plants. A recent study tried to identify genes responsible for this effect by studying messenger RNA expression through RNA sequencing in peach and almond plants, before and after grafting and before and after PPV infection. In this work, we used the same peach and almond samples but focused the high-throughput analyses on small RNA (sRNA) expression. We studied massive sequencing data and found an interesting pattern of sRNA overexpression linked to antiviral defense genes that suggested activation of these genes followed by downregulation to basal levels. We also discovered that 'Garrigues' almond plants were infected by different plant viruses that were transferred to peach plants. The large amounts of viral sRNA found in grafted peaches indicated a strong RNA silencing antiviral response and led us to postulate that these plant viruses could be collaborating in the observed "Garrigues effect."

Gene Expression Analysis of Induced Plum pox virus (Sharka) Resistance in Peach (Prunus persica) by Almond (P. dulcis) Grafting.

No natural sources of resistance to Plum pox virus (PPV, sharka disease) have been identified in peach. However, previous studies have demonstrated that grafting a "Garrigues" almond scion onto "GF305" peach rootstock seedlings heavily infected with PPV can progressively reduce disease symptoms and virus accumulation. Furthermore, grafting a "Garrigues" scion onto the "GF305" rootstock has been shown to completely prevent virus infection. This study aims to analyse the rewiring of gene expression associated with this resistance to PPV transmitted by grafting through the phloem using RNA-Seq and RT-qPCR analysis. A total of 18 candidate genes were differentially expressed after grafting "Garrigues" almond onto healthy "GF305" peach. Among the up-regulated genes, a HEN1 homolog stands out, which, together with the differential expression of RDR- and DCL2-homologs, suggests that the RNA silencing machinery is activated by PPV infection and can contribute to the resistance induced by "Garrigues" almond. Glucan endo-1,3-beta D-glucosidase could be also relevant for the "Garrigues"-induced response, since its expression is much higher in "Garrigues" than in "GF305". We also discuss the potential relevance of the following in PPV infection and "Garrigues"-induced resistance: several pathogenesis-related proteins; no apical meristem proteins; the transcription initiation factor, TFIIB; the speckle-type POZ protein; in addition to a number of proteins involved in phytohormone signalling.

Pedigree analysis of 220 almond genotypes reveals two world mainstream breeding lines based on only three different cultivars.

Loss of genetic variability is an increasing challenge in tree breeding programs due to the repeated use of a reduced number of founder genotypes. However, in almond, little is known about the genetic variability in current breeding stocks, although several cases of inbreeding depression have been reported. To gain insights into the genetic structure in modern breeding programs worldwide, marker-verified pedigree data of 220 almond cultivars and breeding selections were analyzed. Inbreeding coefficients, pairwise relatedness, and genetic contribution were calculated for these genotypes. The results reveal two mainstream breeding lines based on three cultivars: "Tuono", "Cristomorto", and "Nonpareil". Descendants from "Tuono" or "Cristomorto" number 76 (sharing 34 descendants), while "Nonpareil" has 71 descendants. The mean inbreeding coefficient of the analyzed genotypes was 0.041, with 14 genotypes presenting a high inbreeding coefficient, over 0.250. Breeding programs from France, the USA, and Spain showed inbreeding coefficients of 0.075, 0.070, and 0.037, respectively. According to their genetic contribution, modern cultivars from Israel, France, the USA, Spain, and Australia trace back to a maximum of six main founding genotypes. Among the group of 65 genotypes carrying the Sf allele for self-compatibility, the mean relatedness coefficient was 0.125, with "Tuono" as the main founding genotype (24.7% of total genetic contribution). The results broaden our understanding about the tendencies followed in almond breeding over the last 50 years and will have a large impact into breeding decision-making process worldwide. Increasing current genetic variability is required in almond breeding programs to assure genetic gain and continuing breeding progress.

Development of a highly efficient Axiom™ 70 K SNP array for Pyrus and evaluation for high-density mapping and germplasm characterization.

BACKGROUND: Both a source of diversity and the development of genomic tools, such as reference genomes and molecular markers, are equally important to enable faster progress in plant breeding. Pear (Pyrus spp.) lags far behind other fruit and nut crops in terms of employment of available genetic resources for new cultivar development. To address this gap, we designed a high-density, high-efficiency and robust single nucleotide polymorphism (SNP) array for pear, with the main objectives of conducting genetic diversity and genome-wide association studies. RESULTS: By applying a two-step design process, which consisted of the construction of a first 'draft' array for the screening of a small subset of samples, we were able to identify the most robust and informative SNPs to include in the Applied Biosystems™ Axiom™ Pear 70 K Genotyping Array, currently the densest SNP array for pear. Preliminary evaluation of this 70 K array in 1416 diverse pear accessions from the USDA National Clonal Germplasm Repository (NCGR) in Corvallis, OR identified 66,616 SNPs (93% of all the tiled SNPs) as high quality and polymorphic (PolyHighResolution). We further used the Axiom Pear 70 K Genotyping Array to construct high-density linkage maps in a bi-parental population, and to make a direct comparison with available genotyping-by-sequencing (GBS) data, which suggested that the SNP array is a more robust method of screening for SNPs than restriction enzyme reduced representation sequence-based genotyping. CONCLUSIONS: The Axiom Pear 70 K Genotyping Array, with its high efficiency in a widely diverse panel of Pyrus species and cultivars, represents a valuable resource for a multitude of molecular studies in pear. The characterization of the USDA-NCGR collection with this array will provide important information for pear geneticists and breeders, as well as for the optimization of conservation strategies for Pyrus.

A new genomic tool for walnut (Juglans regia L.): development and validation of the high-density Axiom™ J. regia 700K SNP genotyping array.

Over the last 20 years, global production of Persian walnut (Juglans regia L.) has grown enormously, likely reflecting increased consumption due to its numerous benefits to human health. However, advances in genome-wide association (GWA) studies and genomic selection (GS) for agronomically important traits in walnut remain limited due to the lack of powerful genomic tools. Here, we present the development and validation of a high-density 700K single nucleotide polymorphism (SNP) array in Persian walnut. Over 609K high-quality SNPs have been thoroughly selected from a set of 9.6 m genome-wide variants, previously identified from the high-depth re-sequencing of 27 founders of the Walnut Improvement Program (WIP) of University of California, Davis. To validate the effectiveness of the array, we genotyped a collection of 1284 walnut trees, including 1167 progeny of 48 WIP families and 26 walnut cultivars. More than half of the SNPs (55.7%) fell in the highest quality class of 'Poly High Resolution' (PHR) polymorphisms, which were used to assess the WIP pedigree integrity. We identified 151 new parent-offspring relationships, all confirmed with the Mendelian inheritance test. In addition, we explored the genetic variability among cultivars of different origin, revealing how the varieties from Europe and California were differentiated from Asian accessions. Both the reconstruction of the WIP pedigree and population structure analysis confirmed the effectiveness of the Applied Biosystems™ Axiom™ J. regia 700K SNP array, which initiates a novel genomic and advanced phase in walnut genetics and breeding.

DNA Methylation Analysis of Dormancy Release in Almond (Prunus dulcis) Flower Buds Using Epi-Genotyping by Sequencing.

DNA methylation and histone post-translational modifications have been described as epigenetic regulation mechanisms involved in developmental transitions in plants, including seasonal changes in fruit trees. In species like almond (Prunus dulcis (Mill.) D.A: Webb), prolonged exposure to cold temperatures is required for dormancy release and flowering. Aiming to identify genomic regions with differential methylation states in response to chill accumulation, we carried out Illumina reduced-representation genome sequencing on bisulfite-treated DNA from floral buds. To do this, we analyzed almond genotypes with different chilling requirements and flowering times both before and after dormancy release for two consecutive years. The study was performed using epi-Genotyping by Sequencing (epi-GBS). A total of 7317 fragments were sequenced and the samples compared. Out of these fragments, 677 were identified as differentially methylated between the almond genotypes. Mapping these fragments using the Prunus persica (L.) Batsch v.2 genome as reference provided information about coding regions linked to early and late flowering methylation markers. Additionally, the methylation state of ten gene-coding sequences was found to be linked to the dormancy release process.

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.

Novel Insights into Tree Biology and Genome Evolution as Revealed Through Genomics.

Reference genome sequences are the key to the discovery of genes and gene families that determine traits of interest. Recent progress in sequencing technologies has enabled a rapid increase in genome sequencing of tree species, allowing the dissection of complex characters of economic importance, such as fruit and wood quality and resistance to biotic and abiotic stresses. Although the number of reference genome sequences for trees lags behind those for other plant species, it is not too early to gain insight into the unique features that distinguish trees from nontree plants. Our review of the published data suggests that, although many gene families are conserved among herbaceous and tree species, some gene families, such as those involved in resistance to biotic and abiotic stresses and in the synthesis and transport of sugars, are often expanded in tree genomes. As the genomes of more tree species are sequenced, comparative genomics will further elucidate the complexity of tree genomes and how this relates to traits unique to trees.

Assessing the Gene Content of the Megagenome: Sugar Pine (Pinus lambertiana).

Sugar pine (Pinus lambertiana Douglas) is within the subgenus Strobus with an estimated genome size of 31 Gbp. Transcriptomic resources are of particular interest in conifers due to the challenges presented in their megagenomes for gene identification. In this study, we present the first comprehensive survey of the P. lambertiana transcriptome through deep sequencing of a variety of tissue types to generate more than 2.5 billion short reads. Third generation, long reads generated through PacBio Iso-Seq have been included for the first time in conifers to combat the challenges associated with de novo transcriptome assembly. A technology comparison is provided here to contribute to the otherwise scarce comparisons of second and third generation transcriptome sequencing approaches in plant species. In addition, the transcriptome reference was essential for gene model identification and quality assessment in the parallel project responsible for sequencing and assembly of the entire genome. In this study, the transcriptomic data were also used to address questions surrounding lineage-specific Dicer-like proteins in conifers. These proteins play a role in the control of transposable element proliferation and the related genome expansion in conifers.

Sequence of the Sugar Pine Megagenome.

Until very recently, complete characterization of the megagenomes of conifers has remained elusive. The diploid genome of sugar pine (Pinus lambertiana Dougl.) has a highly repetitive, 31 billion bp genome. It is the largest genome sequenced and assembled to date, and the first from the subgenus Strobus, or white pines, a group that is notable for having the largest genomes among the pines. The genome represents a unique opportunity to investigate genome "obesity" in conifers and white pines. Comparative analysis of P. lambertiana and P. taeda L. reveals new insights on the conservation, age, and diversity of the highly abundant transposable elements, the primary factor determining genome size. Like most North American white pines, the principal pathogen of P. lambertiana is white pine blister rust (Cronartium ribicola J.C. Fischer ex Raben.). Identification of candidate genes for resistance to this pathogen is of great ecological importance. The genome sequence afforded us the opportunity to make substantial progress on locating the major dominant gene for simple resistance hypersensitive response, Cr1 We describe new markers and gene annotation that are both tightly linked to Cr1 in a mapping population, and associated with Cr1 in unrelated sugar pine individuals sampled throughout the species' range, creating a solid foundation for future mapping. This genomic variation and annotated candidate genes characterized in our study of the Cr1 region are resources for future marker-assisted breeding efforts as well as for investigations of fundamental mechanisms of invasive disease and evolutionary response.

The walnut (Juglans regia) genome sequence reveals diversity in genes coding for the biosynthesis of non-structural polyphenols.

The Persian walnut (Juglans regia L.), a diploid species native to the mountainous regions of Central Asia, is the major walnut species cultivated for nut production and is one of the most widespread tree nut species in the world. The high nutritional value of J. regia nuts is associated with a rich array of polyphenolic compounds, whose complete biosynthetic pathways are still unknown. A J. regia genome sequence was obtained from the cultivar 'Chandler' to discover target genes and additional unknown genes. The 667-Mbp genome was assembled using two different methods (SOAPdenovo2 and MaSuRCA), with an N50 scaffold size of 464 955 bp (based on a genome size of 606 Mbp), 221 640 contigs and a GC content of 37%. Annotation with MAKER-P and other genomic resources yielded 32 498 gene models. Previous studies in walnut relying on tissue-specific methods have only identified a single polyphenol oxidase (PPO) gene (JrPPO1). Enabled by the J. regia genome sequence, a second homolog of PPO (JrPPO2) was discovered. In addition, about 130 genes in the large gallate 1-ß-glucosyltransferase (GGT) superfamily were detected. Specifically, two genes, JrGGT1 and JrGGT2, were significantly homologous to the GGT from Quercus robur (QrGGT), which is involved in the synthesis of 1-O-galloyl-ß-d-glucose, a precursor for the synthesis of hydrolysable tannins. The reference genome for J. regia provides meaningful insight into the complex pathways required for the synthesis of polyphenols. The walnut genome sequence provides important tools and methods to accelerate breeding and to facilitate the genetic dissection of complex traits.

YeATSAM analysis of the walnut and chickpea transcriptome reveals key genes undetected by current annotation tools.

Background: The transcriptome, a treasure trove of gene space information, remains severely under-used by current genome annotation methods.  Methods: Here, we present an annotation method in the YeATS suite (YeATSAM), based on information encoded by the transcriptome, that demonstrates artifacts of the assembler, which must be addressed to achieve proper annotation.  Results and Discussion: YeATSAM was applied to the transcriptome obtained from twenty walnut tissues and compared to MAKER-P annotation of the recently published walnut genome sequence (WGS). MAKER-P and YeATSAM both failed to annotate several hundred proteins found by the other. Although many of these unannotated proteins have repetitive sequences (possibly transposable elements), other crucial proteins were excluded by each method. An egg cell-secreted protein and a homer protein were undetected by YeATSAM, although these did not produce any transcripts. Importantly, MAKER-P failed to classify key photosynthesis-related proteins, which we show emanated from Trinity assembly artifacts potentially not handled by MAKER-P. Also, no proteins from the large berberine bridge enzyme (BBE) family were annotated by MAKER-P. BBE is implicated in biosynthesis of several alkaloids metabolites, like anti-microbial berberine. As further validation, YeATSAM identified ~1000 genes that are not annotated in the NCBI database by Gnomon. YeATSAM used a RNA-seq derived chickpea ( Cicer arietinum L.) transcriptome assembled using Newbler v2.3.  Conclusions: Since the current version of YeATSAM does not have an ab initio module, we suggest a combined annotation scheme using both MAKER-P and YeATSAM to comprehensively and accurately annotate the WGS.

A Consensus Genetic Map for Pinus taeda and Pinus elliottii and Extent of Linkage Disequilibrium in Two Genotype-Phenotype Discovery Populations of Pinus taeda.

A consensus genetic map for Pinus taeda (loblolly pine) and Pinus elliottii (slash pine) was constructed by merging three previously published P. taeda maps with a map from a pseudo-backcross between P. elliottii and P. taeda. The consensus map positioned 3856 markers via genotyping of 1251 individuals from four pedigrees. It is the densest linkage map for a conifer to date. Average marker spacing was 0.6 cM and total map length was 2305 cM. Functional predictions of mapped genes were improved by aligning expressed sequence tags used for marker discovery to full-length P. taeda transcripts. Alignments to the P. taeda genome mapped 3305 scaffold sequences onto 12 linkage groups. The consensus genetic map was used to compare the genome-wide linkage disequilibrium in a population of distantly related P. taeda individuals (ADEPT2) used for association genetic studies and a multiple-family pedigree used for genomic selection (CCLONES). The prevalence and extent of LD was greater in CCLONES as compared to ADEPT2; however, extended LD with LGs or between LGs was rare in both populations. The average squared correlations, r(2), between SNP alleles less than 1 cM apart were less than 0.05 in both populations and r(2) did not decay substantially with genetic distance. The consensus map and analysis of linkage disequilibrium establish a foundation for comparative association mapping and genomic selection in P. taeda and P. elliottii.

Decoding the massive genome of loblolly pine using haploid DNA and novel assembly strategies.

BACKGROUND: The size and complexity of conifer genomes has, until now, prevented full genome sequencing and assembly. The large research community and economic importance of loblolly pine, Pinus taeda L., made it an early candidate for reference sequence determination. RESULTS: We develop a novel strategy to sequence the genome of loblolly pine that combines unique aspects of pine reproductive biology and genome assembly methodology. We use a whole genome shotgun approach relying primarily on next generation sequence generated from a single haploid seed megagametophyte from a loblolly pine tree, 20-1010, that has been used in industrial forest tree breeding. The resulting sequence and assembly was used to generate a draft genome spanning 23.2 Gbp and containing 20.1 Gbp with an N50 scaffold size of 66.9 kbp, making it a significant improvement over available conifer genomes. The long scaffold lengths allow the annotation of 50,172 gene models with intron lengths averaging over 2.7 kbp and sometimes exceeding 100 kbp in length. Analysis of orthologous gene sets identifies gene families that may be unique to conifers. We further characterize and expand the existing repeat library based on the de novo analysis of the repetitive content, estimated to encompass 82% of the genome. CONCLUSIONS: In addition to its value as a resource for researchers and breeders, the loblolly pine genome sequence and assembly reported here demonstrates a novel approach to sequencing the large and complex genomes of this important group of plants that can now be widely applied.