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.     

Effect of polyamines on in vitro anther culture of Citrus clementina Hort. ex Tan

The improvement of the induction rate in Citrus anther culture is important for taking practical advantage of the haploid potential in breeding. The influence of polyamines on anther culture of Citrus clementina, cv Nules, with particular attention to the free, soluble and insoluble-conjugated polyamine levels, has been investigated. Putrescine, spermidine and putrescine plus spermidine, were added to the standard induction medium. Before culture, spermidine was the most abundant among the free polyamines detected in anthers. The exogenous supply of either putrescine or spermidine, either independently or combined, effected greater uptake and accumulation of polyamines. The addition of 2 mM spermidine to the medium stimulated gametic embryogenesis in clementine Nules, whereas putrescine did not influence embryo production. Regenerants were mostly tri-haploids; a few doubled-haploids and no haploid plants were obtained.     

Gynogenetic haploids of Citrus after in vitro pollination with triploid pollen grains

This study reports haploid plantlet regeneration through gynogenesis in Citrus clementina Hort. ex Tan., cv. Nules, induced by in vitro pollination with pollen grains of Oroblanco, a triploid cultivar of grapefruit. It indicates that parthenogenesis induced in vitro by triploid pollen can be an alternative method to obtain haploids in monoembryonic cultivars of Citrus. Actually, despite considerable efforts, androgenesis has not been yet successful in many genotypes of Citrus. Pollination and mature stage of pistils was necessary for gynogenic embryo regeneration. Fourteen haploid gynogenic embryos of Nules clementine were obtained. Embryo conversion was high (85.7%) and embryos vigorously germinated producing twelve plantlets. Chromosome counting, performed in the roots of obtained embryos, showed the haploid level (n=x=9). Isozyme analyses confirmed the expected homozygous state of embryos and plantlets. This revised version was published online in June 2006 with corrections to the Cover Date.     

Evidence and evolutionary analysis of ancient whole-genome duplication in barley predating the divergence from rice

Background  

Well preserved genomic colinearity among agronomically important grass species such as rice, maize, Sorghum, wheat and barley provides access to whole-genome structure information even in species lacking a reference genome sequence. We investigated footprints of whole-genome duplication (WGD) in barley that shaped the cereal ancestor genome by analyzing shared synteny with rice using a ~2000 gene-based barley genetic map and the rice genome reference sequence.     

Accelerated gene evolution and subfunctionalization in the pseudotetraploid frog Xenopus laevis

Background  

Ancient whole genome duplications have been implicated in the vertebrate and teleost radiations, and in the emergence of diverse angiosperm lineages, but the evolutionary response to such a perturbation is still poorly understood. The African clawed frog Xenopus laevis experienced a relatively recent tetraploidization ~40 million years ago. Analysis of the considerable amount of EST sequence available for this species together with the genome sequence of the related diploid Xenopus tropicalis provides a unique opportunity to study the genomic response to whole genome duplication.     

<Emphasis Type="Italic">Medicago truncatula</Emphasis> contains a second gene encoding a plastid located glutamine synthetase exclusively expressed in developing seeds

Background  

Nitrogen is a crucial nutrient that is both essential and rate limiting for plant growth and seed production. Glutamine synthetase (GS), occupies a central position in nitrogen assimilation and recycling, justifying the extensive number of studies that have been dedicated to this enzyme from several plant sources. All plants species studied to date have been reported as containing a single, nuclear gene encoding a plastid located GS isoenzyme per haploid genome. This study reports the existence of a second nuclear gene encoding a plastid located GS in Medicago truncatula.     

Chaos game representation for comparison of whole genomes

Background  

Chaos game representation of genome sequences has been used for visual representation of genome sequence patterns as well as alignment-free comparisons of sequences based on oligonucleotide frequencies. However the potential of this representation for making alignment-based comparisons of whole genome sequences has not been exploited.     

Genome-wide SNP identification by high-throughput sequencing and selective mapping allows sequence assembly positioning using a framework genetic linkage map

Background  

Determining the position and order of contigs and scaffolds from a genome assembly within an organism's genome remains a technical challenge in a majority of sequencing projects. In order to exploit contemporary technologies for DNA sequencing, we developed a strategy for whole genome single nucleotide polymorphism sequencing allowing the positioning of sequence contigs onto a linkage map using the bin mapping method.     

Machine learning methods for metabolic pathway prediction

Background  

A key challenge in systems biology is the reconstruction of an organism's metabolic network from its genome sequence. One strategy for addressing this problem is to predict which metabolic pathways, from a reference database of known pathways, are present in the organism, based on the annotated genome of the organism.     

Rapid attainment of a doubled haploid line from transgenic maize (<Emphasis Type="Italic">Zea mays</Emphasis> L.) plants by means of anther culture

Summary We present a strategy for establishing a transgenic doubled haploid maize line from heterozygous transgenic material by means of anther culture. Compared to conventional inbreeding, the in vitro androgenesis technique enables a faster generation of virtually fully homozygous lines. Since the androgenic response is highly genotype-dependent, we crossed transgenic, non-androgenic plants carrying a herbicide resistance marker gene (pat, encoding for phosphinothricin acetyl transferase) with a highly androgenic genotype. The transgenic progenies were used as donor plants for anther culture. One transgenic and three non-transgenic doubled haploid lines have been established within approximately 1 yr. The homozygosity of all four doubled haploid lines was tested by analysis of simple sequence repeat (SSR) markers at 19 different loci. Polymorphisms were found between the lines but not within the lines indicating the homozygous nature of the entire plant genome gained by anther culture. Southern blot analysis revealed that the transgenic donor plants and their doubled haploid progeny exhibited the same integration pattern of the pat gene. No segregation of the herbicide resistance trait has been observed among the progeny of the transgenic doubled haploid line.     

Rapid and sensitive detection of Citrus Bacterial Canker by loop-mediated isothermal amplification combined with simple visual evaluation methods

Background  

Citrus Bacterial Canker (CBC) is a major, highly contagious disease of citrus plants present in many countries in Asia, Africa and America, but not in the Mediterranean area. There are three types of Citrus Bacterial Canker, named A, B, and C that have different genotypes and posses variation in host range within citrus species. The causative agent for type A CBC is Xanthomonas citri subsp. citri, while Xanthomonas fuscans subsp. aurantifolii, strain B causes type B CBC and Xanthomonas fuscans subsp. aurantifolii strain C causes CBC type C. The early and accurate identification of those bacteria is essential for the protection of the citrus industry. Detection methods based on bacterial isolation, antibodies or polymerase chain reaction (PCR) have been developed previously; however, these approaches may be time consuming, laborious and, in the case of PCR, it requires expensive laboratory equipment. Loop-mediated isothermal amplification (LAMP), which is a novel isothermal DNA amplification technique, is sensitive, specific, fast and requires no specialized laboratory equipment.     

Inheritance in doubled-diploid clementine and comparative study with SDR unreduced gametes of diploid clementine

Key message

Tetraploid clementine displays mainly tetrasomic inheritance. Genetic structures of 2n SDR and 2 × gametes from DD clementine are complementary and will guides triploids citrus breeding strategies.

Abstract

Triploid breeding is developed worldwide to create new seedless cultivars. Citrus triploid hybrids can be recovered from 2x × 2x sexual hybridizations as a consequence of the formation of unreduced gametes (2n), or from 4x × 2x interploid hybridizations in which tetraploid parents used are most often doubled-diploid (DD). Here we have analyzed the inheritance in doubled-diploid clementine and compared the genetic structures of gametes of DD clementine with SDR unreduced gametes of diploid clementine. Parental heterozygosity restitution (PHR) with DD parents depends on the rate of preferential chromosome pairing and thus the proportion of disomic versus tetrasomic segregations. Doubled-diploid clementine largely exhibited tetrasomic segregation. However, three linkage groups had intermediate segregation and one had a tendency for disomy. Significant doubled reduction rates (DR) rates were observed in six of the nine LGs. Differences of PHR between 2n SDR and 2x DD gametes were highest in the centromeric region and progressively decreased toward the distal regions where they were not significant. Over all markers, PHR was lower (two-thirds) in SDR 2n gametes than in DD-derived diploid gametes. The two strategies appear complementary in terms of genotypic variability. Interploid 4x × 2x hybridization is potentially more efficient for developing new cultivars that are phenotypically closer to the diploid parent of the DD than sexual hybridization through SDR 2n gametes. Conversely, 2x × 2x triploidisation has the potential to produce novel products with characteristics for market segmentation strategies.

     

Construction of a virtual Mycobacterium tuberculosis consensus genome and its application to data from a next generation sequencer

Background

Although Mycobacterium tuberculosis isolates are consisted of several different lineages and the epidemiology analyses are usually assessed relative to a particular reference genome, M. tuberculosis H37Rv, which might introduce some biased results. Those analyses are essentially based genome sequence information of M. tuberculosis and could be performed in sillico in theory, with whole genome sequence (WGS) data available in the databases and obtained by next generation sequencers (NGSs). As an approach to establish higher resolution methods for such analyses, whole genome sequences of the M. tuberculosis complexes (MTBCs) strains available on databases were aligned to construct virtual reference genome sequences called the consensus sequence (CS), and evaluated its feasibility in in sillico epidemiological analyses.

Results

The consensus sequence (CS) was successfully constructed and utilized to perform phylogenetic analysis, evaluation of read mapping efficacy, which is crucial for detecting single nucleotide polymorphisms (SNPs), and various MTBC typing methods virtually including spoligotyping, VNTR, Long sequence polymorphism and Beijing typing. SNPs detected based on CS, in comparison with H37Rv, were utilized in concatemer-based phylogenetic analysis to determine their reliability relative to a phylogenetic tree based on whole genome alignment as the gold standard. Statistical comparison of phylogenic trees based on CS with that of H37Rv indicated the former showed always better results that that of later. SNP detection and concatenation with CS was advantageous because the frequency of crucial SNPs distinguishing among strain lineages was higher than those of H37Rv. The number of SNPs detected was lower with the consensus than with the H37Rv sequence, resulting in a significant reduction in computational time. Performance of each virtual typing was satisfactory and accorded with those published when those are available.

Conclusions

These results indicated that virtual CS constructed from genome sequence data is an ideal approach as a reference for MTBC studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1368-9) contains supplementary material, which is available to authorized users.     

Barcodes for genomes and applications

Background  

Each genome has a stable distribution of the combined frequency for each k-mer and its reverse complement measured in sequence fragments as short as 1000 bps across the whole genome, for 1<k<6. The collection of these k-mer frequency distributions is unique to each genome and termed the genome's barcode.