Search

Transposable elements (TEs) are a major force in the production of new alleles during domestication; nevertheless, their use in association studies has been limited because of their complexity. We have developed a TE genotyping pipeline (TEmarker) and applied it to whole-genome genome-wide association study (GWAS) data from 176 Oryza sativa subsp. japonica accessions to identify genetic elements associated with specific agronomic traits. TE markers recovered a large proportion (69%) of single-nucleotide polymorphism (SNP)-based GWAS peaks, and these TE peaks retained ca. 25% of the SNPs. The use of TEs in GWASs may reduce false positives associated with linkage disequilibrium (LD) among SNP markers. A genome scan revealed positive selection on TEs associated with agronomic traits. We found several cases of insertion and deletion variants that potentially resulted from the direct action of TEs, including an allele of LOC_Os11g08410 associated with plant height and panicle length traits. Together, these findings reveal the utility of TE markers for connecting genotype to phenotype and suggest a potential role for TEs in influencing phenotypic variations in rice that impact agronomic traits. School of Plant and Environmental Sciences at Virginia Tech Published version The authors want to acknowledge the School of Plant and Environmental Sciences at Virginia Tech for the support of H.Y. They also want to acknowledge the Advanced Research Computing unit (www.arc.vt.edu) for the use of computational resources. Finally, the authors would like to acknowledge the reviewers for their useful suggestions that helped to improve the article.

Background K-seq, a new genotyping methodology based on the amplification of genomic regions using two steps of Klenow amplification with short oligonucleotides, followed by standard PCR and Illumina sequencing, is presented. The protocol was accompanied by software developed to aid with primer set design. Results As the first examples, K-seq in species as diverse as tomato, dog and wheat was developed. K-seq provided genetic distances similar to those based on WGS in dogs. Experiments comparing K-seq and GBS in tomato showed similar genetic results, although K-seq had the advantage of finding more SNPs for the same number of Illumina reads. The technology reproducibility was tested with two independent runs of the tomato samples, and the correlation coefficient of the SNP coverages between samples was 0.8 and the genotype match was above 94%. K-seq also proved to be useful in polyploid species. The wheat samples generated specific markers for all subgenomes, and the SNPs generated from the diploid ancestors were located in the expected subgenome with accuracies greater than 80%. Conclusion K-seq is an open, patent-unencumbered, easy-to-set-up, cost-effective and reliable technology ready to be used by any molecular biology laboratory without special equipment in many genetic studies. Published version

Background Mango, Mangifera indica L., an important tropical fruit crop, is grown for its sweet and aromatic fruits. Past improvement of this species has predominantly relied on chance seedlings derived from over 1000 cultivars in the Indian sub-continent with a large variation for fruit size, yield, biotic and abiotic stress resistance, and fruit quality among other traits. Historically, mango has been an orphan crop with very limited molecular information. Only recently have molecular and genomics-based analyses enabled the creation of linkage maps, transcriptomes, and diversity analysis of large collections. Additionally, the combined analysis of genomic and phenotypic information is poised to improve mango breeding efficiency. Results This study sequenced, de novo assembled, analyzed, and annotated the genome of the monoembryonic mango cultivar ‘Tommy Atkins’. The draft genome sequence was generated using NRGene de-novo Magic on high molecular weight DNA of ‘Tommy Atkins’, supplemented by 10X Genomics long read sequencing to improve the initial assembly. A hybrid population between ‘Tommy Atkins’ x ‘Kensington Pride’ was used to generate phased haplotype chromosomes and a highly resolved phased SNP map. The final ‘Tommy Atkins’ genome assembly was a consensus sequence that included 20 pseudomolecules representing the 20 chromosomes of mango and included ~ 86% of the ~ 439 Mb haploid mango genome. Skim sequencing identified ~ 3.3 M SNPs using the ‘Tommy Atkins’ x ‘Kensington Pride’ mapping population. Repeat masking identified 26,616 genes with a median length of 3348 bp. A whole genome duplication analysis revealed an ancestral 65 MYA polyploidization event shared with Anacardium occidentale. Two regions, one on LG4 and one on LG7 containing 28 candidate genes, were associated with the commercially important fruit size characteristic in the mapping population. Conclusions The availability of the complete ‘Tommy Atkins’ mango genome will aid global initiatives to study mango genetics. Published version

Additional file 3: Supplemental Data S3. Fruit weight LG4 QTL report. Supplemental Data S4. Fruit weight LG7 QTL report. Supplemental Table 1. Summary of the M. indica genome assembly by chromosome. Supplemental Table 2. BUSCO analysis for the consensus diploid genome assembly TA4. Supplemental Table 3. Occurrence and distribution of repetitive DNA sequences in the mango genome. Supplemental Table 4. Protein domain content. Supplemental Table 5. Number of homologs pairs with Ks between 0.1 and 0.6 per chromosome. Supplemental Table 6. Repetitive elements comparative analysis. Supplemental Table 7. Comparison of the SNP variants and variants rate for ‘Tommy Atkins’ and ‘Kensington Pride’ by pseudomolecule. Supplemental Table S8. Concurrence of linkage groups between ‘Tommy Atkins’ and ‘Alphonso’. Supplementary Figure 1. Repeat composition for M. indica, P. vera and C. sinensis

Background Understanding the DNA methylome and its relationship with non-coding RNAs, including microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), is essential for elucidating the molecular mechanisms underlying key biological processes in plants. Few studies have examined the functional roles of the DNA methylome in grass species with highly heterozygous polyploid genomes. Results We performed genome-wide DNA methylation profiling in the tetraploid switchgrass (Panicum virgatum L.) cultivar ‘Alamo’ using bisulfite sequencing. Single-base-resolution methylation patterns were observed in switchgrass leaf and root tissues, which allowed for characterization of the relationship between DNA methylation and mRNA, miRNA, and lncRNA populations. The results of this study revealed that siRNAs positively regulate DNA methylation of the mCHH sites surrounding genes, and that DNA methylation interferes with gene and lncRNA expression in switchgrass. Ninety-six genes covered by differentially methylated regions (DMRs) were annotated by GO analysis as being involved in stimulus-related processes. Functionally, 82% (79/96) of these genes were found to be hypomethylated in switchgrass root tissue. Sequencing analysis of lncRNAs identified two lncRNAs that are potential precursors of miRNAs, which are predicted to target genes that function in cellulose biosynthesis, stress regulation, and stem and root development. Conclusions This study characterized the DNA methylome in switchgrass and elucidated its relevance to gene and non-coding RNAs. These results provide valuable genomic resources and references that will aid further epigenetic research in this important biofuel crop. Published version

Additional file 8: Table S3. Comparison of methylation levels between class I and class II transposons in leaf and root tissues.

Additional file 16: Table S7. Correlation analysis between methylation levels and gene expression in different genic regions.

Background Tobacco (Nicotiana tabacum) is an important plant model system that has played a key role in the early development of molecular plant biology. The tobacco genome is large and its characterisation challenging because it is an allotetraploid, likely arising from hybridisation between diploid N. sylvestris and N. tomentosiformis ancestors. A draft assembly was recently published for N. tabacum, but because of the aforementioned genome complexities it was of limited utility due to a high level of fragmentation. Results Here we report an improved tobacco genome assembly, which, aided by the application of optical mapping, achieves an N50 size of 2.17 Mb and enables anchoring of 64% of the genome to pseudomolecules; a significant increase from the previous value of 19%. We use this assembly to identify two homeologous genes that explain the differentiation of the burley tobacco market class, with potential for greater understanding of Nitrogen Utilization Efficiency and Nitrogen Use Efficiency in plants; an important trait for future sustainability of agricultural production. Conclusions Development of an improved genome assembly for N. tabacum enables what we believe to be the first successful map-based gene discovery for the species, and demonstrates the value of an improved assembly for future research in this model and commercially-important species. Published version

Plant forms display a wide variety of architectures, depending on the number of lateral branches, internode elongation and phyllotaxy. These are in turn determined by the number, the position and the fate of the Axillary Meristems (AMs). Mutants that affect AM determination during the vegetative phase have been isolated in several model plants. Among these genes, the GRAS transcription factor LATERAL SUPPRESSOR (Ls) plays a pivotal role in AM determination during the vegetative phase. Hereby we characterize the phylogenetic orthologue of Ls in Antirrhinum, ERAMOSA (ERA). Our data supported ERA control of AM formation during both the vegetative and the reproductive phase in snapdragon. A phylogenetic analysis combined with an analysis of the synteny of Ls in several species strongly supported the hypothesis that ERA is a phylogenetic orthologue of Ls, although it plays a broader role. During the reproductive phase ERA promotes the establishment of the stem niche at the bract axis but, after the reproductive transition, it is antagonized by the MADS box transcription factor SQUAMOSA (SQUA). Surprisingly double mutant era squa plants display a squa phenotype developing axillary meristems, which can eventually turn into inflorescences or flowers. Cariplo Foundation [2011-2257]; FP7-PEOPLE-2013-IRSES [FRUIT LOOK] We thank Xue Yongbiao (Beijing Institute of Genomics, The Chinese Academy of Sciences) for sharing data in advance of publication. Edward Kiegle (University of Milan) is acknowledged for editing the paper. This work was supported by Cariplo Foundation [grant number 2011-2257 to C.M. and S.M.] and by FP7-PEOPLE-2013-IRSES [FRUIT LOOK to S.M.].

BACKGROUND: Switchgrass (Panicum virgatum L.) is a warm-season perennial grass that can be used as a second generation bioenergy crop. However, foliar fungal pathogens, like switchgrass rust, have the potential to significantly reduce switchgrass biomass yield. Despite its importance as a prominent bioenergy crop, a genome-wide comprehensive analysis of NB-LRR disease resistance genes has yet to be performed in switchgrass. RESULTS: In this study, we used a homology-based computational approach to identify 1011 potential NB-LRR resistance gene homologs (RGHs) in the switchgrass genome (v 1.1). In addition, we identified 40 RGHs that potentially contain unique domains including major sperm protein domain, jacalin-like binding domain, calmodulin-like binding, and thioredoxin. RNA-sequencing analysis of leaf tissue from 'Alamo', a rust-resistant switchgrass cultivar, and 'Dacotah', a rust-susceptible switchgrass cultivar, identified 2634 high quality variants in the RGHs between the two cultivars. RNA-sequencing data from field-grown cultivar 'Summer' plants indicated that the expression of some of these RGHs was developmentally regulated. CONCLUSIONS: Our results provide useful insight into the molecular structure, distribution, and expression patterns of members of the NB-LRR gene family in switchgrass. These results also provide a foundation for future work aimed at elucidating the molecular mechanisms underlying disease resistance in this important bioenergy crop. Published version

RNA-seq has been used to perform global expression analysis of the achene and the receptacle at four stages of fruit ripening, and of the roots and leaves of strawberry (Fragaria × ananassa). About 967 million reads and 191 Gb of sequence were produced, using Illumina sequencing. Mapping the reads in the related genome of the wild diploid Fragaria vesca revealed differences between the achene and receptacle development program, and reinforced the role played by ethylene in the ripening receptacle. For the strawberry transcriptome assembly, a de novo strategy was followed, generating separate assemblies for each of the ten tissues and stages sampled. The Trinity program was used for these assemblies, resulting in over 1.4 M isoforms. Filtering by a threshold of 0.3 FPKM, and doing Blastx (E-value < 1 e-30) against the UniProt database of plants reduced the number to 472,476 isoforms. Their assembly with the MIRA program (90% homology) resulted in 26,087 contigs. From these, 91.34 percent showed high homology to Fragaria vesca genes and 87.30 percent Fragaria iinumae (BlastN E-value < 1 e-100). Mapping back the reads on the MIRA contigs identified polymorphisms at nucleotide level, using FREEBAYES, as well as estimate their relative abundance in each sample.

Background Switchgrass (Panicum virgatum L.) is a warm-season perennial grass that can be used as a second generation bioenergy crop. However, foliar fungal pathogens, like switchgrass rust, have the potential to significantly reduce switchgrass biomass yield. Despite its importance as a prominent bioenergy crop, a genome-wide comprehensive analysis of NB-LRR disease resistance genes has yet to be performed in switchgrass. Results In this study, we used a homology-based computational approach to identify 1011 potential NB-LRR resistance gene homologs (RGHs) in the switchgrass genome (v 1.1). In addition, we identified 40 RGHs that potentially contain unique domains including major sperm protein domain, jacalin-like binding domain, calmodulin-like binding, and thioredoxin. RNA-sequencing analysis of leaf tissue from ‘Alamo’, a rust-resistant switchgrass cultivar, and ‘Dacotah’, a rust-susceptible switchgrass cultivar, identified 2634 high quality variants in the RGHs between the two cultivars. RNA-sequencing data from field-grown cultivar ‘Summer’ plants indicated that the expression of some of these RGHs was developmentally regulated. Conclusions Our results provide useful insight into the molecular structure, distribution, and expression patterns of members of the NB-LRR gene family in switchgrass. These results also provide a foundation for future work aimed at elucidating the molecular mechanisms underlying disease resistance in this important bioenergy crop. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-3201-5) contains supplementary material, which is available to authorized users.

Background Solanum dulcamara (bittersweet, climbing nightshade) is one of the few species of the Solanaceae family native to Europe. As a common weed it is adapted to a wide range of ecological niches and it has long been recognized as one of the alternative hosts for pathogens and pests responsible for many important diseases in potato, such as Phytophthora. At the same time, it may represent an alternative source of resistance genes against these diseases. Despite its unique ecology and potential as a genetic resource, genomic research tools are lacking for S. dulcamara. We have taken advantage of next-generation sequencing to speed up research on and use of this non-model species. Results In this work, we present the first large-scale characterization of the S. dulcamara transcriptome. Through comparison of RNAseq reads from two different accessions, we were able to predict transcript-based SNP and SSR markers. Using the SNP markers in combination with genomic AFLP and CAPS markers, the first genome-wide genetic linkage map of bittersweet was generated. Based on gene orthology, the markers were anchored to the genome of related Solanum species (tomato, potato and eggplant), revealing both conserved and novel chromosomal rearrangements. This allowed a better estimation of the evolutionary moment of rearrangements in a number of cases and showed that chromosomal breakpoints are regularly re-used. Conclusion Knowledge and tools developed as part of this study pave the way for future genomic research and exploitation of this wild Solanum species. The transcriptome assembly represents a resource for functional analysis of genes underlying interesting biological and agronomical traits and, in the absence of the full genome, provides a reference for RNAseq gene expression profiling aimed at understanding the unique biology of S. dulcamara. Cross-species orthology-based marker selection is shown to be a powerful tool to quickly generate a comparative genetic map, which may speed up gene mapping and contribute to the understanding of genome evolution within the Solanaceae family.