15 years of GDR: New data and functionality in the Genome Database for Rosaceae.

The Genome Database for Rosaceae (GDR, https://www.rosaceae.org) is an integrated web-based community database resource providing access to publicly available genomics, genetics and breeding data and data-mining tools to facilitate basic, translational and applied research in Rosaceae. The volume of data in GDR has increased greatly over the last 5 years. The GDR now houses multiple versions of whole genome assembly and annotation data from 14 species, made available by recent advances in sequencing technology. Annotated and searchable reference transcriptomes, RefTrans, combining peer-reviewed published RNA-Seq as well as EST datasets, are newly available for major crop species. Significantly more quantitative trait loci, genetic maps and markers are available in MapViewer, a new visualization tool that better integrates with other pages in GDR. Pathways can be accessed through the new GDR Cyc Pathways databases, and synteny among the newest genome assemblies from eight species can be viewed through the new synteny browser, SynView. Collated single-nucleotide polymorphism diversity data and phenotypic data from publicly available breeding datasets are integrated with other relevant data. Also, the new Breeding Information Management System allows breeders to upload, manage and analyze their private breeding data within the secure GDR server with an option to release data publicly.

Soil indigenous microbiome and plant genotypes cooperatively modify soybean rhizosphere microbiome assembly.

BACKGROUND: Plants have evolved intimate interactions with soil microbes for a range of beneficial functions including nutrient acquisition, pathogen resistance and stress tolerance. Further understanding of this system is a promising way to advance sustainable agriculture by exploiting the versatile benefits offered by the plant microbiome. The rhizosphere is the interface between plant and soil, and functions as the first step of plant defense and root microbiome recruitment. It features a specialized microbial community, intensive microbe-plant and microbe-microbe interactions, and complex signal communication. To decipher the rhizosphere microbiome assembly of soybean (Glycine max), we comprehensively characterized the soybean rhizosphere microbial community using 16S rRNA gene sequencing and evaluated the structuring influence from both host genotype and soil source. RESULTS: Comparison of the soybean rhizosphere to bulk soil revealed significantly different microbiome composition, microbe-microbe interactions and metabolic capacity. Soil type and soybean genotype cooperatively modulated microbiome assembly with soil type predominantly shaping rhizosphere microbiome assembly while host genotype slightly tuned this recruitment process. The undomesticated progenitor species, Glycine soja, had higher rhizosphere diversity in both soil types tested in comparison to the domesticated soybean genotypes. Rhizobium, Novosphingobium, Phenylobacterium, Streptomyces, Nocardioides, etc. were robustly enriched in soybean rhizosphere irrespective of the soil tested. Co-occurrence network analysis revealed dominant soil type effects and genotype specific preferences for key microbe-microbe interactions. Functional prediction results demonstrated converged metabolic capacity in the soybean rhizosphere between soil types and among genotypes, with pathways related to xenobiotic degradation, plant-microbe interactions and nutrient transport being greatly enriched in the rhizosphere. CONCLUSION: This comprehensive comparison of the soybean microbiome between soil types and genotypes expands our understanding of rhizosphere microbe assembly in general and provides foundational information for soybean as a legume crop for this assembly process. The cooperative modulating role of the soil type and host genotype emphasizes the importance of integrated consideration of soil condition and plant genetic variability for future development and application of synthetic microbiomes. Additionally, the detection of the tuning role by soybean genotype in rhizosphere microbiome assembly provides a promising way for future breeding programs to integrate host traits participating in beneficial microbiota assembly.

The essential chloroplast ribosomal protein uL15c interacts with the chloroplast RNA helicase ISE2 and affects intercellular trafficking through plasmodesmata.

Chloroplasts retain part of their ancestral genomes and the machinery for expression of those genomes. The nucleus-encoded chloroplast RNA helicase INCREASED SIZE EXCLUSION LIMIT2 (ISE2) is required for chloroplast ribosomal RNA processing and chloro-ribosome assembly. To further elucidate ISE2's role in chloroplast translation, two independent approaches were used to identify its potential protein partners. Both a yeast two-hybrid screen and a pull-down assay identified plastid ribosomal protein L15, uL15c (formerly RPL15), as interacting with ISE2. The interaction was confirmed in vivo by co-immunoprecipitation. Interestingly, we found that rpl15 null mutants do not complete embryogenesis, indicating that RPL15 is an essential gene for autotrophic growth of Arabidopsis thaliana. Arabidopsis and Nicotiana benthamiana plants with reduced expression of RPL15 developed chlorotic leaves, had reduced photosynthetic capacity and exhibited defective chloroplast development. Processing of chloroplast ribosomal RNAs and assembly of ribosomal subunits were disrupted by reduced expression of RPL15. Chloroplast translation was also decreased, reducing accumulation of chloroplast-encoded proteins, in such plants compared to wild-type plants. Notably, knockdown of RPL15 expression increased intercellular trafficking, a phenotype also observed in plants with reduced ISE2 expression. This finding provides further evidence for chloroplast function in modulating intercellular trafficking via plasmodesmata.

Development of Genomic Resources for the Powdery Mildew, Erysiphe pulchra.

Powdery mildews (PMs) are important plant pathogens causing widespread damage. Here, we report the first draft genome of Erysiphe pulchra, the causative agent of PM of flowering dogwood, Cornus florida. The assembled genome was 63.5 Mbp and resulted in formation of 19,442 contigs (N50 = 11,686 bp) that contained an estimated 6,860 genes with a genome coverage of 62×. We found 102 candidate secreted effector proteins (CSEPs) in E. pulchra similar to E. necator genes that are potentially involved in disease development. This draft genome is an initial step for understanding the evolutionary history of the PMs and will also provide insight into evolutionary strategies that led to the wide host expansion and environmental adaptations so effectively employed by the PM lineages.

A novel picornavirus-like genome from transcriptome sequencing of sugar beet cyst nematode represents a new putative genus.

Analysis of transcriptome sequence data from eggs and second-stage juveniles (J2s) of sugar beet cyst nematode (SBCN, Heterodera schachtii) identified the full-length genome of a positive-sense single-stranded RNA virus, provisionally named sugar beet cyst nematode virus 1 (SBCNV1). The SBCNV1 sequence was detected in both eggs and J2s, indicating its possible vertical transmission. The 9503-nucleotide genome sequence contains a single long open reading frame, which was predicted to encode a polyprotein with conserved domains for picornaviral structural proteins proximal to its amino terminus and RNA helicase, cysteine proteinase and RNA-dependent RNA polymerase (RdRp) conserved domains proximal to its carboxyl terminus, hallmarks of viruses belonging to the order Picornavirales. Phylogenetic analysis of the predicted SBCNV1 RdRp amino acid sequence indicated that the SBCNV1 sequence is most closely related to members of the family Secoviridae, which includes genera of nematode-transmitted plant-infecting viruses. SBCNV1 represents the first fully sequenced viral genome from SBCN.

Genome-wide association study identifies loci associated with milk leukocyte phenotypes following experimental challenge with Streptococcus uberis.

Mastitis is a detrimental disease in the dairy industry that decreases milk quality and costs upwards of $2 billion annually. Often, mastitis results from bacteria entering the gland through the teat opening. Streptococcus uberis is responsible for a high percentage of subclinical and clinical mastitis. Following an intramammary experimental challenge with S. uberis on Holstein cows (n = 40), milk samples were collected and somatic cell counts (SCC) were determined by the Dairy Herd Improvement Association Laboratory. Traditional genome-wide association studies (GWAS) have utilized test day SCC or SCC lactation averages to identify loci of interest. Our approach utilizes SCC collected following a S. uberis experimental challenge to generate three novel phenotypes: (1) area under the curve (AUC) of SCC for 0-7 days and (2) 0-28 days post-challenge; and (3) when SCC returned to below 200,000 cells/mL post-challenge (< 21 days, 21-28 days, or > 28 days). Polymorphisms were identified using Illumina's BovineSNP50 v2 DNA BeadChip. Associations were tested using Plink software and identified 16 significant (p < 1.0 × 10-4) single-nucleotide polymorphisms (SNPs) across the phenotypes. Most significant SNPs were in genes linked to cell signaling, migration, and apoptosis. Several have been recognized in relation to infectious processes (ATF7, SGK1, and PACRG), but others less so (TRIO, GLRA1, CELSR2, TIAM2, CPE). Further investigation of these genes and their roles in inflammation (e.g., SCC) can provide potential targets that influence resolution of mammary gland infection. Likewise, further investigation of the identified SNP with mastitis and other disease phenotypes can provide greater insight to the potential of these SNP as genetic markers.

Genotypic and phenotypic evaluation of off-type grasses in hybrid Bermudagrass [Cynodon dactylon (L.) Pers. x C. transvaalensis Burtt-Davy] putting greens using genotyping-by-sequencing and morphological characterization.

BACKGROUND: Interspecific hybrid bermudagrass [Cynodon dactylon (L.) Pers. x C. transvaalensis Burtt-Davy] is one of the most widely used grasses on golf courses, with cultivars derived from 'Tifgreen' or 'Tifdwarf' particularly used for putting greens. Many bermudagrass cultivars established for putting greens can be genetically unstable and lead to the occurrence of undesirable off-type grasses that vary in phenotype. The objective of this research was to genetically and phenotypically differentiate off-type grasses and hybrid cultivars. Beginning in 2013, off-type and desirable hybrid bermudagrass samples were collected from golf course putting greens in the southeastern United States and genetically and phenotypically characterized using genotyping-by-sequencing and morphology. RESULTS: Genotyping-by-sequencing determined that 11% (5) of off-type and desirable samples from putting greens were genetically divergent from standard cultivars such as Champion, MiniVerde, Tifdwarf, TifEagle, and Tifgreen. In addition, genotyping-by-sequencing was unable to genetically distinguish all standard cultivars from one another due to their similar origin and clonal propagation; however, over 90,000 potentially informative nucleotide variants were identified among the triploid hybrid cultivars. CONCLUSIONS: Although few genetic differences were found in this research, samples harvested from golf course putting greens had variable morphology and were clustered into three distinct phenotypic groups. The majority of off-type grasses in hybrid bermudagrass putting greens were genetically similar with variable morphological traits. Off-type grasses within golf course putting greens have the potential to compromise putting surface functionality and aesthetics.

High-quality genetic mapping with ddRADseq in the non-model tree Quercus rubra.

BACKGROUND: Restriction site associated DNA sequencing (RADseq) has the potential to be a broadly applicable, low-cost approach for high-quality genetic linkage mapping in forest trees lacking a reference genome. The statistical inference of linear order must be as accurate as possible for the correct ordering of sequence scaffolds and contigs to chromosomal locations. Accurate maps also facilitate the discovery of chromosome segments containing allelic variants conferring resistance to the biotic and abiotic stresses that threaten forest trees worldwide. We used ddRADseq for genetic mapping in the tree Quercus rubra, with an approach optimized to produce a high-quality map. Our study design also enabled us to model the results we would have obtained with less depth of coverage. RESULTS: Our sequencing design produced a high sequencing depth in the parents (248×) and a moderate sequencing depth (15×) in the progeny. The digital normalization method of generating a de novo reference and the SAMtools SNP variant caller yielded the most SNP calls (78,725). The major drivers of map inflation were multiple SNPs located within the same sequence (77% of SNPs called). The highest quality map was generated with a low level of missing data (5%) and a genome-wide threshold of 0.025 for deviation from Mendelian expectation. The final map included 849 SNP markers (1.8% of the 78,725 SNPs called). Downsampling the individual FASTQ files to model lower depth of coverage revealed that sequencing the progeny using 96 samples per lane would have yielded too few SNP markers to generate a map, even if we had sequenced the parents at depth 248×. CONCLUSIONS: The ddRADseq technology produced enough high-quality SNP markers to make a moderately dense, high-quality map. The success of this project was due to high depth of coverage of the parents, moderate depth of coverage of the progeny, a good framework map, an optimized bioinformatics pipeline, and rigorous premapping filters. The ddRADseq approach is useful for the construction of high-quality genetic maps in organisms lacking a reference genome if the parents and progeny are sequenced at sufficient depth. Technical improvements in reduced representation sequencing (RRS) approaches are needed to reduce the amount of missing data.

The green ash transcriptome and identification of genes responding to abiotic and biotic stresses.

BACKGROUND: To develop a set of transcriptome sequences to support research on environmental stress responses in green ash (Fraxinus pennsylvanica), we undertook deep RNA sequencing of green ash tissues under various stress treatments. The treatments, including emerald ash borer (EAB) feeding, heat, drought, cold and ozone, were selected to mimic the increasing threats of climate change and invasive pests faced by green ash across its native habitat. RESULTS: We report the generation and assembly of RNA sequences from 55 green ash samples into 107,611 putative unique transcripts (PUTs). 52,899 open reading frames were identified. Functional annotation of the PUTs by comparison to the Uniprot protein database identified matches for 63 % of transcripts and for 98 % of transcripts with ORFs. Further functional annotation identified conserved protein domains and assigned gene ontology terms to the PUTs. Examination of transcript expression across different RNA libraries revealed that expression patterns clustered based on tissues regardless of stress treatment. The transcripts from stress treatments were further examined to identify differential expression. Tens to hundreds of differentially expressed PUTs were identified for each stress treatment. A set of 109 PUTs were found to be consistently up or down regulated across three or more different stress treatments, representing basal stress response candidate genes in green ash. In addition, 1956 simple sequence repeats were identified in the PUTs, of which we identified 465 high quality DNA markers and designed flanking PCR primers. CONCLUSIONS: North American native ash trees have suffered extensive mortality due to EAB infestation, creating a need to breed or select for resistant green ash genotypes. Stress from climate change is an additional concern for longevity of native ash populations. The use of genomics could accelerate management efforts. The green ash transcriptome we have developed provides important sequence information, genetic markers and stress-response candidate genes.

Diversity of ABC transporter genes across the plant kingdom and their potential utility in biotechnology.

BACKGROUND: The ATP-binding cassette (ABC) transporter gene superfamily is ubiquitous among extant organisms and prominently represented in plants. ABC transporters act to transport compounds across cellular membranes and are involved in a diverse range of biological processes. Thus, the applicability to biotechnology is vast, including cancer resistance in humans, drug resistance among vertebrates, and herbicide and other xenobiotic resistance in plants. In addition, plants appear to harbor the highest diversity of ABC transporter genes compared with any other group of organisms. This study applied transcriptome analysis to survey the kingdom-wide ABC transporter diversity in plants and suggest biotechnology applications of this diversity. RESULTS: We utilized sequence similarity-based informatics techniques to infer the identity of ABC transporter gene candidates from 1295 phylogenetically-diverse plant transcriptomes. A total of 97,149 putative (approximately 25 % were full-length) ABC transporter gene members were identified; each RNA-Seq library (plant sample) had 88 ± 30 gene members. As expected, simpler organisms, such as algae, had fewer unique members than vascular land plants. Differences were also noted in the richness of certain ABC transporter subfamilies. Land plants had more unique ABCB, ABCC, and ABCG transporter gene members on average (p < 0.005), and green algae, red algae, and bryophytes had significantly more ABCF transporter gene members (p < 0.005). Ferns had significantly fewer ABCA transporter gene members than all other plant groups (p < 0.005). CONCLUSIONS: We present a transcriptomic overview of ABC transporter gene members across all major plant groups. An increase in the number of gene family members present in the ABCB, ABCC, and ABCD transporter subfamilies may indicate an expansion of the ABC transporter superfamily among green land plants, which include all crop species. The striking difference between the number of ABCA subfamily transporter gene members between ferns and other plant taxa is surprising and merits further investigation. Discussed is the potential exploitation of ABC transporters in plant biotechnology, with an emphasis on crops.

Developing expressed sequence tag libraries and the discovery of simple sequence repeat markers for two species of raspberry (Rubus L.).

BACKGROUND: Due to a relatively high level of codominant inheritance and transferability within and among taxonomic groups, simple sequence repeat (SSR) markers are important elements in comparative mapping and delineation of genomic regions associated with traits of economic importance. Expressed sequence tags (ESTs) are a source of SSRs that can be used to develop markers to facilitate plant breeding and for more basic research across genera and higher plant orders. METHODS: Leaf and meristem tissue from 'Heritage' red raspberry (Rubus idaeus) and 'Bristol' black raspberry (R. occidentalis) were utilized for RNA extraction. After conversion to cDNA and library construction, ESTs were sequenced, quality verified, assembled and scanned for SSRs.  Primers flanking the SSRs were designed and a subset tested for amplification, polymorphism and transferability across species. ESTs containing SSRs were functionally annotated using the GenBank non-redundant (nr) database and further classified using the gene ontology database. RESULTS: To accelerate development of EST-SSRs in the genus Rubus (Rosaceae), 1149 and 2358 cDNA sequences were generated from red raspberry and black raspberry, respectively. The cDNA sequences were screened using rigorous filtering criteria which resulted in the identification of 121 and 257 SSR loci for red and black raspberry, respectively. Primers were designed from the surrounding sequences resulting in 131 and 288 primer pairs, respectively, as some sequences contained more than one SSR locus. Sequence analysis revealed that the SSR-containing genes span a diversity of functions and share more sequence identity with strawberry genes than with other Rosaceous species. CONCLUSION: This resource of Rubus-specific, gene-derived markers will facilitate the construction of linkage maps composed of transferable markers for studying and manipulating important traits in this economically important genus.

Mapping QTLs for blight resistance and morpho-phenological traits in inter-species hybrid families of chestnut (Castanea spp.).

Chestnut blight (caused by Cryphonectria parasitica), together with Phytophthora root rot (caused by Phytophthora cinnamomi), has nearly extirpated American chestnut (Castanea dentata) from its native range. In contrast to the susceptibility of American chestnut, many Chinese chestnut (C. mollissima) genotypes are resistant to blight. In this research, we performed a series of genome-wide association studies for blight resistance originating from three unrelated Chinese chestnut trees (Mahogany, Nanking and M16) and a Quantitative Trait Locus (QTL) study on a Mahogany-derived inter-species F2 family. We evaluated trees for resistance to blight after artificial inoculation with two fungal strains and scored nine morpho-phenological traits that are the hallmarks of species differentiation between American and Chinese chestnuts. Results support a moderately complex genetic architecture for blight resistance, as 31 QTLs were found on 12 chromosomes across all studies. Additionally, although most morpho-phenological trait QTLs overlap or are adjacent to blight resistance QTLs, they tend to aggregate in a few genomic regions. Finally, comparison between QTL intervals for blight resistance and those previously published for Phytophthora root rot resistance, revealed five common disease resistance regions on chromosomes 1, 5, and 11. Our results suggest that it will be difficult, but still possible to eliminate Chinese chestnut alleles for the morpho-phenological traits while achieving relatively high blight resistance in a backcross hybrid tree. We see potential for a breeding scheme that utilizes marker-assisted selection early for relatively large effect QTLs followed by genome selection in later generations for smaller effect genomic regions.

Chemical profile and analysis of biosynthetic pathways and genes of volatile terpenes in Pityopsis ruthii, a rare and endangered flowering plant.

It is critical to gather biological information about rare and endangered plants to incorporate into conservation efforts. The secondary metabolism of Pityopsis ruthii, an endangered flowering plant that only occurs along limited sections of two rivers (Ocoee and Hiwassee) in Tennessee, USA was studied. Our long-term goal is to understand the mechanisms behind P. ruthii's adaptation to restricted areas in Tennessee. Here, we profiled the secondary metabolites, specifically in flowers, with a focus on terpenes, aiming to uncover the genomic and molecular basis of terpene biosynthesis in P. ruthii flowers using transcriptomic and biochemical approaches. By comparative profiling of the nonpolar portion of metabolites from various tissues, P. ruthii flowers were rich in terpenes, which included 4 monoterpenes and 10 sesquiterpenes. These terpenes were emitted from flowers as volatiles with monoterpenes and sesquiterpenes accounting for almost 68% and 32% of total emission of terpenes, respectively. These findings suggested that floral terpenes play important roles for the biology and adaptation of P. ruthii to its limited range. To investigate the biosynthesis of floral terpenes, transcriptome data for flowers were produced and analyzed. Genes involved in the terpene biosynthetic pathway were identified and their relative expressions determined. Using this approach, 67 putative terpene synthase (TPS) contigs were detected. TPSs in general are critical for terpene biosynthesis. Seven full-length TPS genes encoding putative monoterpene and sesquiterpene synthases were cloned and functionally characterized. Three catalyzed the biosynthesis of sesquiterpenes and four catalyzed the biosynthesis of monoterpenes. In conclusion, P. ruthii plants employ multiple TPS genes for the biosynthesis of a mixture of floral monoterpenes and sesquiterpenes, which probably play roles in chemical defense and attracting insect pollinators alike.

Chloroplast genome of the invasive Pyrus calleryana underscores the high molecular diversity of the species.

The complete chloroplast genome of Pyrus calleryana (GenBank OM541581.1) was developed by de novo assembly from whole-genome sequencing data. Reference-guided (P. phaeocarpa) read mapping and assembly were followed by annotation and phylogenetic comparisons. The 159,965 bp P. calleryana chloroplast genome represented 36.56% GC content with a classical quadripartite architecture and two inverted repeats regions (IRs; each 26,392 bp) separating the large single-copy region (LSC; 87,942 bp) and the small single-copy region (SSC; 19.239 bp). In total, 125 unique features were annotated in that genome, including 83 protein coding genes, 38 tRNA coding genes, and 4 rRNA coding genes. Phylogenetic analyses based on the whole chloroplast genome sequences placed the P. calleryana among other Rosaceae plants, specifically among the Asian species of Pyrus.

Characterization and microsatellite marker development for a common bark and ambrosia beetle associate, Geosmithia obscura.

Symbioses between Geosmithia fungi and wood-boring and bark beetles seldom result in disease induction within the plant host. Yet, exceptions exist such as Geosmithia morbida, the causal agent of Thousand Cankers Disease (TCD) of walnuts and wingnuts, and Geosmithia sp. 41, the causal agent of Foamy Bark Canker disease of oaks. Isolates of G. obscura were recovered from black walnut trees in eastern Tennessee and at least one isolate induced cankers following artificial inoculation. Due to the putative pathogenicity and lack of recovery of G. obscura from natural lesions, a molecular diagnostic screening tool was developed using microsatellite markers mined from the G. obscura genome. A total of 3256 candidate microsatellite markers were identified (2236, 789, 137 di-, tri-, and tetranucleotide motifs, respectively), with 2011, 703, 101 di-, tri-, and tetranucleotide motifs, respectively, containing markers with primers. From these, 75 microsatellite markers were randomly selected, screened, and optimized, resulting in 28 polymorphic markers that yielded single, consistently recovered bands, which were used in downstream analyses. Five of these microsatellite markers were found to be specific to G. obscura and did not cross-amplify into other, closely related species. Although the remaining tested markers could be useful, they cross-amplified within different Geosmithia species, making them not reliable for G. obscura detection. Five novel microsatellite markers (GOBS9, GOBS10, GOBS41, GOBS43, and GOBS50) were developed based on the G. obscura genome. These species-specific microsatellite markers are available as a tool for use in molecular diagnostics and can assist future surveillance studies.

Development and Characterization of 20 Genomic SSR Markers for Ornamental Cultivars of Weigela.

Weigela (Caprifoliaceae) is a genus of ornamental plants popular for its phenotypic variation and hardiness, that includes species hybridized to produce the commercially available cultivars. Despite its popularity, limited genetic resources exist for the genus. Twenty genomic simple sequence repeat (gSSR) markers distributed across the genome were developed using low coverage whole-genome sequencing data of Weigela Spilled Wine®. A cross-amplification evaluation with these 20 gSSR markers on a collection of 18 Weigela cultivars revealed a total of 111 unique alleles, including 36 private alleles. A diagrammatic key was constructed to identify cultivars using only six of the gSSR markers, demonstrating the newly developed gSSR markers are immediately useful for cultivar identification. Future uses could include breeding with marker-assisted selection, determining the history of hybridization of the current cultivated lines, aiding in the construction of genetic maps, and assessing the patterns of population genetic structure of Weigela spp.

A genetically anchored physical framework for Theobroma cacao cv. Matina 1-6.

BACKGROUND: The fermented dried seeds of Theobroma cacao (cacao tree) are the main ingredient in chocolate. World cocoa production was estimated to be 3 million tons in 2010 with an annual estimated average growth rate of 2.2%. The cacao bean production industry is currently under threat from a rise in fungal diseases including black pod, frosty pod, and witches' broom. In order to address these issues, genome-sequencing efforts have been initiated recently to facilitate identification of genetic markers and genes that could be utilized to accelerate the release of robust T. cacao cultivars. However, problems inherent with assembly and resolution of distal regions of complex eukaryotic genomes, such as gaps, chimeric joins, and unresolvable repeat-induced compressions, have been unavoidably encountered with the sequencing strategies selected. RESULTS: Here, we describe the construction of a BAC-based integrated genetic-physical map of the T. cacao cultivar Matina 1-6 which is designed to augment and enhance these sequencing efforts. Three BAC libraries, each comprised of 10× coverage, were constructed and fingerprinted. 230 genetic markers from a high-resolution genetic recombination map and 96 Arabidopsis-derived conserved ortholog set (COS) II markers were anchored using pooled overgo hybridization. A dense tile path consisting of 29,383 BACs was selected and end-sequenced. The physical map consists of 154 contigs and 4,268 singletons. Forty-nine contigs are genetically anchored and ordered to chromosomes for a total span of 307.2 Mbp. The unanchored contigs (105) span 67.4 Mbp and therefore the estimated genome size of T. cacao is 374.6 Mbp. A comparative analysis with A. thaliana, V. vinifera, and P. trichocarpa suggests that comparisons of the genome assemblies of these distantly related species could provide insights into genome structure, evolutionary history, conservation of functional sites, and improvements in physical map assembly. A comparison between the two T. cacao cultivars Matina 1-6 and Criollo indicates a high degree of collinearity in their genomes, yet rearrangements were also observed. CONCLUSIONS: The results presented in this study are a stand-alone resource for functional exploitation and enhancement of Theobroma cacao but are also expected to complement and augment ongoing genome-sequencing efforts. This resource will serve as a template for refinement of the T. cacao genome through gap-filling, targeted re-sequencing, and resolution of repetitive DNA arrays.

Sequencing of a QTL-rich region of the Theobroma cacao genome using pooled BACs and the identification of trait specific candidate genes.

BACKGROUND: BAC-based physical maps provide for sequencing across an entire genome or a selected sub-genomic region of biological interest. Such a region can be approached with next-generation whole-genome sequencing and assembly as if it were an independent small genome. Using the minimum tiling path as a guide, specific BAC clones representing the prioritized genomic interval are selected, pooled, and used to prepare a sequencing library. RESULTS: This pooled BAC approach was taken to sequence and assemble a QTL-rich region, of ~3 Mbp and represented by twenty-seven BACs, on linkage group 5 of the Theobroma cacao cv. Matina 1-6 genome. Using various mixtures of read coverages from paired-end and linear 454 libraries, multiple assemblies of varied quality were generated. Quality was assessed by comparing the assembly of 454 reads with a subset of ten BACs individually sequenced and assembled using Sanger reads. A mixture of reads optimal for assembly was identified. We found, furthermore, that a quality assembly suitable for serving as a reference genome template could be obtained even with a reduced depth of sequencing coverage. Annotation of the resulting assembly revealed several genes potentially responsible for three T. cacao traits: black pod disease resistance, bean shape index, and pod weight. CONCLUSIONS: Our results, as with other pooled BAC sequencing reports, suggest that pooling portions of a minimum tiling path derived from a BAC-based physical map is an effective method to target sub-genomic regions for sequencing. While we focused on a single QTL region, other QTL regions of importance could be similarly sequenced allowing for biological discovery to take place before a high quality whole-genome assembly is completed.

Development and Characterization of 15 Novel Genomic SSRs for Viburnum farreri.

The Viburnum genus is of particular interest to horticulturalists, phylogeneticists, and biogeographers. Despite its popularity, there are few existing molecular markers to investigate genetic diversity in this large genus, which includes over 160 species. There are also few polymorphic molecular tools that can delineate closely related species within the genus. Viburnum farreri, a member of the Solenotinus subclade and one of the centers of diversity for Viburnum, was selected for DNA sequencing and development of genomic simple sequence repeats (gSSRs). In this study, 15 polymorphic gSSRs were developed and characterized for a collection of 19 V. farreri samples. Number of alleles per locus ranged from two- to- eight and nine loci had four or more alleles. Observed heterozygosity ranged from 0 to 0.84 and expected heterozygosity ranged from 0.10 to 0.80 for the 15 loci. Shannon diversity index values across these loci ranged from 0.21 to 1.62. The markers developed in this study add to the existing molecular toolkit for the genus and will be used in future studies investigating cross-transferability, genetic variation, and species and cultivar delimitation in the Viburnum genus and closely allied genera in the Adoxaceae and Caprifoliaceae.

Genome-Wide Changes of Regulatory Non-Coding RNAs Reveal Pollen Development Initiated at Ecodormancy in Peach.

Bud dormancy is under the regulation of complex mechanisms including genetic and epigenetic factors. To study the function of regulatory non-coding RNAs in winter dormancy release, we analyzed the small RNA and long non-coding RNA (lncRNA) expression from peach (Prunus persica) floral buds in endodormancy, ecodormancy and bud break stages. Small RNAs underwent a major shift in expression primarily between dormancy and flowering with specific pairs of microRNAs and their mRNA target genes undergoing coordinated differential expression. From endodormancy to ecodormancy, ppe-miR6285 was significantly upregulated while its target gene, an ASPARAGINE-RICH PROTEIN involved in the regulation of abscisic acid signaling, was downregulated. At ecodormancy, ppe-miR2275, a homolog of meiosis-specific miR2275 across angiosperms, was significantly upregulated, supporting microsporogenesis in anthers at a late stage of dormancy. The expression of 785 lncRNAs, unlike the overall expression pattern in the small RNAs, demonstrated distinctive expression signatures across all dormancy and flowering stages. We predicted that a subset of lncRNAs were targets of microRNAs and found 18 lncRNA/microRNA target pairs with both differentially expressed across time points. The genome-wide differential expression and network analysis of non-coding RNAs and mRNAs from the same tissues provide new candidate loci for dormancy regulation and suggest complex noncoding RNA interactions control transcriptional regulation across these key developmental time points.