Timber! Felling the loblolly pine genome

Conventional short read sequences derived from haploid DNA were extended into long super-reads enabling assembly of the massive 22 Gbp loblolly pine, Pinus taeda, genome.See related research http://genomebiology.com/2014/15/3/R59     

Decoding the oak genome: public release of sequence data,assembly, annotation and publication strategies

The 1.5 Gbp/2C genome of pedunculate oak (Quercus robur) has been sequenced. A strategy was established for dealing with the challenges imposed by the sequencing of such a large, complex and highly heterozygous genome by a whole‐genome shotgun (WGS) approach, without the use of costly and time‐consuming methods, such as fosmid or BAC clone‐based hierarchical sequencing methods. The sequencing strategy combined short and long reads. Over 49 million reads provided by Roche 454 GS‐FLX technology were assembled into contigs and combined with shorter Illumina sequence reads from paired‐end and mate‐pair libraries of different insert sizes, to build scaffolds. Errors were corrected and gaps filled with Illumina paired‐end reads and contaminants detected, resulting in a total of 17 910 scaffolds (>2 kb) corresponding to 1.34 Gb. Fifty per cent of the assembly was accounted for by 1468 scaffolds (N50 of 260 kb). Initial comparison with the phylogenetically related Prunus persica gene model indicated that genes for 84.6% of the proteins present in peach (mean protein coverage of 90.5%) were present in our assembly. The second and third steps in this project are genome annotation and the assignment of scaffolds to the oak genetic linkage map. In accordance with the Bermuda and Fort Lauderdale agreements and the more recent Toronto Statement, the oak genome data have been released into public sequence repositories in advance of publication. In this presubmission paper, the oak genome consortium describes its principal lines of work and future directions for analyses of the nature, function and evolution of the oak genome.     

Finished sequence and assembly of the DUF1220-rich 1q21 region using a haploid human genome

Background

Although the reference human genome sequence was declared finished in 2003, some regions of the genome remain incomplete due to their complex architecture. One such region, 1q21.1-q21.2, is of increasing interest due to its relevance to human disease and evolution. Elucidation of the exact variants behind these associations has been hampered by the repetitive nature of the region and its incomplete assembly. This region also contains 238 of the 270 human DUF1220 protein domains, which are implicated in human brain evolution and neurodevelopment. Additionally, examinations of this protein domain have been challenging due to the incomplete 1q21 build. To address these problems, a single-haplotype hydatidiform mole BAC library (CHORI-17) was used to produce the first complete sequence of the 1q21.1-q21.2 region.

Results

We found and addressed several inaccuracies in the GRCh37sequence of the 1q21 region on large and small scales, including genomic rearrangements and inversions, and incorrect gene copy number estimates and assemblies. The DUF1220-encoding NBPF genes required the most corrections, with 3 genes removed, 2 genes reassigned to the 1p11.2 region, 8 genes requiring assembly corrections for DUF1220 domains (~91 DUF1220 domains were misassigned), and multiple instances of nucleotide changes that reassigned the domain to a different DUF1220 subtype. These corrections resulted in an overall increase in DUF1220 copy number, yielding a haploid total of 289 copies. Approximately 20 of these new DUF1220 copies were the result of a segmental duplication from 1q21.2 to 1p11.2 that included two NBPF genes. Interestingly, this duplication may have been the catalyst for the evolutionarily important human lineage-specific chromosome 1 pericentric inversion.

Conclusions

Through the hydatidiform mole genome sequencing effort, the 1q21.1-q21.2 region is complete and misassemblies involving inter- and intra-region duplications have been resolved. The availability of this single haploid sequence path will aid in the investigation of many genetic diseases linked to 1q21, including several associated with DUF1220 copy number variations. Finally, the corrected sequence identified a recent segmental duplication that added 20 additional DUF1220 copies to the human genome, and may have facilitated the chromosome 1 pericentric inversion that is among the most notable human-specific genomic landmarks.     

Paternal inheritance of chloroplast DNA and maternal inheritance of mitochondrial DNA in loblolly pine

Summary The inheritance of organelle DNAs in loblolly pine was studied by using restriction fragment length polymorphisms. Chloroplast DNA from loblolly pine is paternally inherited in pitch pine x loblolly pine hybrids. Mitochondrial DNA is maternally inherited in loblolly pine crosses. The uniparental inheritance of organelle genomes from opposite sexes within the same plant appears to be unique among those higher plants that have been tested and indicates that loblolly pine, and possibly other conifers, must have special mechanisms for organelle exclusion or degradation or both. This genetic system creates an exceptional opportunity for the study of maternal and paternal genetic lineages within a single species.     

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 (Almería, 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.     

Bottom-up genome assembly using the Bacillus subtilis genome vector

We established a protocol to construct complete recombinant genomes from their small contiguous DNA pieces and obtained the genomes of mouse mitochondrion and rice chloroplast using a B. subtilis genome (BGM) vector. This method allows the design of any recombinant genomes, valuable not only for fundamental research in systems biology and synthetic biology but also for various applications in the life sciences.     

Gas exchange and canopy structure of 9-year-old loblolly pine,pitch pine and pitch x loblolly hybrids

Summary Seasonal gas exchange and canopy structure were compared among 9-year-old loblolly pine (Pinus taeda L.), pitch pine (Pinus rigida Mill.), and pitch x loblolly hybrids (Pinus rigida x taeda) growing in an F2 plantation located in Critz, Va., USA. Leaf net photosynthesis, conductance, internal CO₂ concentration (c_i), water use efficiency (WUE; photosynthesis/conductance), dark respiration and the ratio of net photosynthesis/respiration did not vary among or within the three taxa. Significant differences in volume production, crown length, total crown leaf surface area and the silhouette area of shade shoots among the taxa were observed. The loblolly-South Carolina source had greater volume and crown surface area than the pitch pine, and the hybrid taxa were intermediate between the two. Although the silhouette area ratio of shade foliage varied among taxa, it was not related to volume. A strong relationship between total leaf surface area and volume was observed. Leaf conductance, c_i, WUE and leaf water potential were the physiological parameters significantly and positively correlated with volume. This study suggests that the amount of needle surface in the canopy is more important in early stand volume growth than the leaf carbon exchange rate and the degree of needle self-shading in the lower canopy.     

Comparative mapping in loblolly and radiata pine using RFLP and microsatellite markers

Genetic linkage maps were constructed for loblolly pine (Pinus taeda L.) and radiata pine (P. radiata D. Don) using a common set of RFLP and microsatellite markers. The map for loblolly pine combined data from two full-sib families and consisted of 20 linkage groups covering 1281 cM. The map for radiata pine had 14 linkage groups and covered 1223 cM. All of the RFLP probes readily hybridise between loblolly and radiata pine often producing similar hybridisation patterns. There were in total 60 homologous RFLP loci mapped in both species which could be used for comparative purposes. A set of 20 microsatellite markers derived from radiata pine were also assayed; however, only 9 amplified and revealed polymorphic loci in both species. Single-locus RFLP and microsatellite markers were used to match up linkage groups and compare order between species. Twelve syntenic groups were obtained each consisting of from 3 to 9 homologous loci. The order of homologous loci was colinear in most cases, suggesting no major chromosomal rearrangements in the evolution of these species. Comparative mapping between loblolly and radiata pine should facilitate genetic research in both species and provide a framework for mapping in other pine species. Received: 25 November 1998 / /Accepted: 19 December 1998     

Phenylalanine ammonia-lyase from loblolly pine : purification of the enzyme and isolation of complementary DNA clones

Phenylalanine ammonia-lyase (EC 4.3.1.5) has been purified from differentiating secondary xylem of loblolly pine (Pinus taeda L.). Native molecular weight of the enzyme was estimated to be 280,000, with a subunit molecular weight of 74,000; isoelectric point, 5.8; and Michaelis constant for i-phenylalanine, 27 micromolar. No evidence was obtained for the existence of isoforms of the enzyme, nor for negative cooperativity of substrate binding. Polyclonal antibodies were raised against the phenylalanine ammonia-lyase subunit and used to identify a pal clone in an expression library of xylem complementary DNA (cDNA). Polymerase chain reaction, using oligonucleotide primers made from N-terminal amino acid sequence and from the 5′ end of the clone isolated from the expression library, was also used to isolate cDNA clones. These methods yielded cDNA clones covering the protein coding region of the pal messenger RNA. Comparisons of nucleotide sequence of pal cDNAs from pine, bean, sweet potato, and rice showed 60 to 62% identity between the pine clone and the angiosperm clones.     

Characterization of EST-SSRs in loblolly pine and spruce

In the first large study of conifer expressed sequence tag-simple sequence repeats (EST-SSRs), two large conifer EST databases were characterized for EST-SSRs. One database was from “interior spruce” (white and Engelmann spruce in Southern British Columbia) and Sitka spruce, while the other was from loblolly pine. We found 475 and 629 unique EST-SSRs in loblolly pine and spruce, respectively. 3′ ESTs contained 14% more SSRs than 5′ EST reads in loblolly pine and 41% more in spruce. Conifer EST-SSRs differed conspicuously from angiosperm EST-SSRs in several aspects. EST-SSRs were considerably less frequent in conifers (one EST-SSR every ∼50 kb) than in angiosperms (one EST-SSR every ∼20 kb). Dinucleotide repeats were the most abundant repeat class in conifers, while in angiosperms, trinucleotides were most common. Finally, the AT motif was the dominant motif recovered in both conifer species, whereas AG was the most common dinucleotide repeat in angiosperms. Also, as these EST-SSRs in conifers could be developed into useful genetic markers, our work demonstrates the value of large-scale EST sequencing projects for in-silico approaches for marker development.     

In vitro regeneration of loblolly pine and random amplified polymorphic DNA analyses of regenerated plantlets

 Adventitious buds were induced from organogenic callus derived from mature zygotic embryos of three lines (E-311, E-440, and E-822) of loblolly pine (Pinus taeda L.) within 27 weeks of culture. The influence of cytokinins, silver nitrate, and low-temperature treatment on the differentiation of adventitious buds was analyzed. Elongation of adventitious buds was achieved on TE medium supplemented with 0.5 mg/l indole-3-butyric acid (IBA) and 1 mg/l 6-benzyladenine (BA). After adventitious shoots had rooted on TE medium supplemented with 0.5 mg/l IBA, 2 mg/l BA, and 0.5 mg/l gibberellic acid, 498 regenerated plantlets were transferred to a perlite:peatmoss:vermiculite (1 :>: 1) soil mixture; 351 of these survived in the field. Total DNA was extracted from 21 regenerated plantlets randomly chosen from the 151 regenerated plantlets of line E-822. Random amplified polymorphic DNA (RAPD) analysis using 80 arbitrary oligonucleotide 10-mers showed that 21 primers gave 107 clear reproducible bands, with the amplification products being monomorphic for all of the plantlets of line E-822 tested. A total of 2,247 bands obtained from these studies exhibited no aberration in RAPD banding patterns among the tested plantlets. These results suggest that somatic organogenesis can be used for clonal micropropagation of some lines of loblolly pine without the fear of the appearance of unwanted somaclonal variants. Received: 5 August 2000 / Revision received: 5 September 2000 / Accepted: 10 October 2000     

Rare instances of haploid inducer DNA in potato dihaploids and ploidy-dependent genome instability

In cultivated tetraploid potato (Solanum tuberosum), reduction to diploidy (dihaploidy) allows for hybridization to diploids and introgression breeding and may facilitate the production of inbreds. Pollination with haploid inducers (HIs) yields maternal dihaploids, as well as triploid and tetraploid hybrids. Dihaploids may result from parthenogenesis, entailing the development of embryos from unfertilized eggs, or genome elimination, entailing missegregation and the loss of paternal chromosomes. A sign of genome elimination is the occasional persistence of HI DNA in some dihaploids. We characterized the genomes of 919 putative dihaploids and 134 hybrids produced by pollinating tetraploid clones with three HIs: IVP35, IVP101, and PL-4. Whole-chromosome or segmental aneuploidy was observed in 76 dihaploids, with karyotypes ranging from 2n = 2x − 1 = 23 to 2n = 2x + 3 = 27. Of the additional chromosomes in 74 aneuploids, 66 were from the non-inducer parent and 8 from the inducer parent. Overall, we detected full or partial chromosomes from the HI parent in 0.87% of the dihaploids, irrespective of parental genotypes. Chromosomal breaks commonly affected the paternal genome in the dihaploid and tetraploid progeny, but not in the triploid progeny, correlating instability to sperm ploidy and to haploid induction. The residual HI DNA discovered in the progeny is consistent with genome elimination as the mechanism of haploid induction.

A large potato progeny population produced by crossing tetraploid cultivated clones to diploid Phureja lines displays rare instances of haploid inducer chromosomes, which are frequently damaged.     

Adventures in the enormous: a 1.8 million clone BAC library for the 21.7 Gb genome of loblolly pine

Loblolly pine (LP; Pinus taeda L.) is the most economically important tree in the U.S. and a cornerstone species in southeastern forests. However, genomics research on LP and other conifers has lagged behind studies on flowering plants due, in part, to the large size of conifer genomes. As a means to accelerate conifer genome research, we constructed a BAC library for the LP genotype 7-56. The LP BAC library consists of 1,824,768 individually-archived clones making it the largest single BAC library constructed to date, has a mean insert size of 96 kb, and affords 7.6X coverage of the 21.7 Gb LP genome. To demonstrate the efficacy of the library in gene isolation, we screened macroarrays with overgos designed from a pine EST anchored on LP chromosome 10. A positive BAC was sequenced and found to contain the expected full-length target gene, several gene-like regions, and both known and novel repeats. Macroarray analysis using the retrotransposon IFG-7 (the most abundant repeat in the sequenced BAC) as a probe indicates that IFG-7 is found in roughly 210,557 copies and constitutes about 5.8% or 1.26 Gb of LP nuclear DNA; this DNA quantity is eight times the Arabidopsis genome. In addition to its use in genome characterization and gene isolation as demonstrated herein, the BAC library should hasten whole genome sequencing of LP via next-generation sequencing strategies/technologies and facilitate improvement of trees through molecular breeding and genetic engineering. The library and associated products are distributed by the Clemson University Genomics Institute (www.genome.clemson.edu).     

Development of DNA methylation-based epigenetic age predictors in loblolly pine (Pinus taeda)

Biological ageing is connected to life history variation across ecological scales and informs a basic understanding of age-related declines in organismal function. Altered DNA methylation dynamics are a conserved aspect of biological ageing and have recently been modelled to predict chronological age among vertebrate species. In addition to their utility in estimating individual age, differences between chronological and predicted ages arise due to acceleration or deceleration of epigenetic ageing, and these discrepancies are linked to disease risk and multiple life history traits. Although evidence suggests that patterns of DNA methylation can describe ageing in plants, predictions with epigenetic clocks have yet to be performed. Here, we resolve the DNA methylome across CpG, CHG, and CHH-methylation contexts in the loblolly pine tree (Pinus taeda) and construct epigenetic clocks capable of predicting ages in this species within 6% of its maximum lifespan. Although patterns of CHH-methylation showed little association with age, both CpG and CHG-methylation contexts were strongly associated with ageing, largely becoming hypomethylated with age. Among age-associated loci were those in close proximity to malate dehydrogenase, NADH dehydrogenase, and 18S and 26S ribosomal RNA genes. This study reports one of the first epigenetic clocks in plants and demonstrates the universality of age-associated DNA methylation dynamics which can inform conservation and management practices, as well as our ecological and evolutionary understanding of biological ageing in plants.     

Morphological and anatomical relationships of loblolly pine fine roots

 Suberized or brown roots have been traditionally considered secondary or woody tissues. The validity of using morphological features such as color to infer root anatomy for southern pines is questionable and unproven. The objectives of this study were (i) to establish relationships between root color, diameter, and developmental stage (i.e., primary or secondary tissues) for loblolly pine, (ii) to determine the percentages of primary and secondary brown roots by diameter class, and (iii) to use these percentages to make first order estimates of the amount of brown root length and surface area that is in the primary and secondary developmental stages for sampled roots of a semi-mature loblolly pine stand. ”Unsectioned” roots were collected by coring to a 25-cm depth 3 times a year and measuring roots for length and surface area by diameter class. ”Sectioned” roots were sampled from a one-time core and from periodic grab samples. These roots were sectioned and characterized by their color, diameter and developmental stage. Diameters of sectioned roots (n=353) ranged from 0.21 to 8.24 mm. White and orange roots ranged from 0.23 to 2.50 mm, while brown roots spanned the range. White roots were developmentally primary, whereas orange/brown roots were either primary (from 0.21 to 2.50 mm), secondary (from 0.33 to 8.24 mm), or in transition (from 0.27 to 0.76). Total live root length of the sampled stands was estimated to be composed of 38% primary tissue, 58% secondary tissue, and 4% transition tissue. Lastly, neither root color nor diameter was a reliable predictor of developmental stage unless roots were white (primary), or orange/brown and >2.5 mm in diameter (secondary). Received: 30 June 1997 / Accepted: 28 January 1998     

The mapping of the barley genome by RAPD analysis using double haploid strains

A map of the barley genome consisting of 107 loci was constructed using double haploid lines. 33 RAPD loci marks all barley chromosomes. The total length of the obtained map is 1047 centiMorgans (cM) with a 9.8 cM average distance between markers. 7 linkage groups was identified as certain chromosomes. The distribution of the markers was analyzed, comparison with other published maps was made. Some statistical considerations are presented. Alternative strategies of gene mapping are considered.