Genetic associations with resistance to Meloidogyne enterolobii in guava (Psidium sp.) using cross-genera SNPs and comparative genomics to Eucalyptus highlight evolutionary conservation across the Myrtaceae.

Tropical fruit tree species constitute a yet untapped supply of outstanding diversity of taste and nutritional value, barely developed from the genetics standpoint, with scarce or no genomic resources to tackle the challenges arising in modern breeding practice. We generated a de novo genome assembly of the Psidium guajava, the super fruit "apple of the tropics", and successfully transferred 14,268 SNP probesets from Eucalyptus to Psidium at the nucleotide level, to detect genomic loci linked to resistance to the root knot nematode (RKN) Meloidogyne enterolobii derived from the wild relative P. guineense. Significantly associated loci with resistance across alternative analytical frameworks, were detected at two SNPs on chromosome 3 in a pseudo-assembly of Psidium guajava genome built using a syntenic path approach with the Eucalyptus grandis genome to determine the order and orientation of the contigs. The P. guineense-derived resistance response to RKN and disease onset is conceivably triggered by mineral nutrients and phytohormone homeostasis or signaling with the involvement of the miRNA pathway. Hotspots of mapped resistance quantitative trait loci and functional annotation in the same genomic region of Eucalyptus provide further indirect support to our results, highlighting the evolutionary conservation of genomes across genera of Myrtaceae in the adaptation to pathogens. Marker assisted introgression of the resistance loci mapped should accelerate the development of improved guava cultivars and hybrid rootstocks.

Transcriptome Profiling-Based Analysis of Carbohydrate-Active Enzymes in Aspergillus terreus Involved in Plant Biomass Degradation.

Given the global abundance of plant biomass residues, potential exists in biorefinery-based applications with lignocellulolytic fungi. Frequently isolated from agricultural cellulosic materials, Aspergillus terreus is a fungus efficient in secretion of commercial enzymes such as cellulases, xylanases and phytases. In the context of biomass saccharification, lignocellulolytic enzyme secretion was analyzed in a strain of A. terreus following liquid culture with sugarcane bagasse (SB) (1% w/v) and soybean hulls (SH) (1% w/v) as sole carbon source, in comparison to glucose (G) (1% w/v). Analysis of the fungal secretome revealed a maximum of 1.017 UI.mL-1 xylanases after growth in minimal medium with SB, and 1.019 UI.mL-1 after incubation with SH as carbon source. The fungal transcriptome was characterized on SB and SH, with gene expression examined in comparison to equivalent growth on G as carbon source. Over 8000 genes were identified, including numerous encoding enzymes and transcription factors involved in the degradation of the plant cell wall, with significant expression modulation according to carbon source. Eighty-nine carbohydrate-active enzyme (CAZyme)-encoding genes were identified following growth on SB, of which 77 were differentially expressed. These comprised 78% glycoside hydrolases, 8% carbohydrate esterases, 2.5% polysaccharide lyases, and 11.5% auxiliary activities. Analysis of the glycoside hydrolase family revealed significant up-regulation for genes encoding 25 different GH family proteins, with predominance for families GH3, 5, 7, 10, and 43. For SH, from a total of 91 CAZyme-encoding genes, 83 were also significantly up-regulated in comparison to G. These comprised 80% glycoside hydrolases, 7% carbohydrate esterases, 5% polysaccharide lyases, 7% auxiliary activities (AA), and 1% glycosyltransferases. Similarly, within the glycoside hydrolases, significant up-regulation was observed for genes encoding 26 different GH family proteins, with predominance again for families GH3, 5, 10, 31, and 43. A. terreus is a promising species for production of enzymes involved in the degradation of plant biomass. Given that this fungus is also able to produce thermophilic enzymes, this first global analysis of the transcriptome following cultivation on lignocellulosic carbon sources offers considerable potential for the application of candidate genes in biorefinery applications.

Spring Is Coming: Genetic Analyses of the Bud Break Date Locus Reveal Candidate Genes From the Cold Perception Pathway to Dormancy Release in Apple (Malus × domestica Borkh.).

Chilling requirement (CR) for bud dormancy completion determines the time of bud break in apple (Malus × domestica Borkh.). The molecular control of bud dormancy is highly heritable, suggesting a strong genetic control of the trait. An available Infinium II SNP platform for genotyping containing 8,788 single nucleotide polymorphic markers was employed, and linkage maps were constructed in a F1 cross from the low CR M13/91 and the moderate CR cv. Fred Hough. These maps were used to identify quantitative trait loci (QTL) for bud break date as a trait related to dormancy release. A major QTL for bud break was detected at the beginning of linkage group 9 (LG9). This QTL remained stable during seven seasons in two different growing sites. To increase mapping efficiency in detecting contributing genes underlying this QTL, 182 additional SNP markers located at the locus for bud break were used. Combining linkage mapping and structural characterization of the region, the high proportion of the phenotypic variance in the trait explained by the QTL is related to the coincident positioning of Arabidopsis orthologs for ICE1, FLC, and PRE1 protein-coding genes. The proximity of these genes from the most explanatory markers of this QTL for bud break suggests potential genetic additive effects, reinforcing the hypothesis of inter-dependent mechanisms controlling dormancy induction and release in apple trees.

Analysis of the Transcriptome in Aspergillus tamarii During Enzymatic Degradation of Sugarcane Bagasse.

The production of bioethanol from non-food agricultural residues represents an alternative energy source to fossil fuels for incorporation into the world's economy. Within the context of bioconversion of plant biomass into renewable energy using improved enzymatic cocktails, Illumina RNA-seq transcriptome profiling was conducted on a strain of Aspergillus tamarii, efficient in biomass polysaccharide degradation, in order to identify genes encoding proteins involved in plant biomass saccharification. Enzyme production and gene expression was compared following growth in liquid and semi-solid culture with steam-exploded sugarcane bagasse (SB) (1% w/v) and glucose (1% w/v) employed as contrasting sole carbon sources. Enzyme production following growth in liquid minimum medium supplemented with SB resulted in 0.626 and 0.711 UI.mL-1 xylanases after 24 and 48 h incubation, respectively. Transcriptome profiling revealed expression of over 7120 genes, with groups of genes modulated according to solid or semi-solid culture, as well as according to carbon source. Gene ontology analysis of genes expressed following SB hydrolysis revealed enrichment in xyloglucan metabolic process and xylan, pectin and glucan catabolic process, indicating up-regulation of genes involved in xylanase secretion. According to carbohydrate-active enzyme (CAZy) classification, 209 CAZyme-encoding genes were identified with significant differential expression on liquid or semi-solid SB, in comparison to equivalent growth on glucose as carbon source. Up-regulated CAZyme-encoding genes related to cellulases (CelA, CelB, CelC, CelD) and hemicellulases (XynG1, XynG2, XynF1, XylA, AxeA, arabinofuranosidase) showed up to a 10-fold log2FoldChange in expression levels. Five genes from the AA9 (GH61) family, related to lytic polysaccharide monooxygenase (LPMO), were also identified with significant expression up-regulation. The transcription factor gene XlnR, involved in induction of hemicellulases, showed up-regulation on liquid and semi-solid SB culture. Similarly, the gene ClrA, responsible for regulation of cellulases, showed increased expression on liquid SB culture. Over 150 potential transporter genes were also identified with increased expression on liquid and semi-solid SB culture. This first comprehensive analysis of the transcriptome of A. tamarii contributes to our understanding of genes and regulatory systems involved in cellulose and hemicellulose degradation in this fungus, offering potential for application in improved enzymatic cocktail development for plant biomass degradation in biorefinery applications.

Design and evaluation of a sequence capture system for genome-wide SNP genotyping in highly heterozygous plant genomes: a case study with a keystone Neotropical hardwood tree genome.

Targeted sequence capture coupled to high-throughput sequencing has become a powerful method for the study of genome-wide sequence variation. Following our recent development of a genome assembly for the Pink Ipê tree (Handroanthus impetiginosus), a widely distributed Neotropical timber species, we now report the development of a set of 24,751 capture probes for single-nucleotide polymorphisms (SNPs) characterization and genotyping across 18,216 distinct loci, sampling more than 10 Mbp of the species genome. This system identifies nearly 200,000 SNPs located inside or in close proximity to almost 14,000 annotated protein-coding genes, generating quality genotypic data in populations spanning wide geographic distances across the species native range. To provide recommendations for future developments of similar systems for highly heterozygous plant genomes we investigated issues such as probe design, sequencing coverage and bioinformatics, including the evaluation of the capture efficiency and a reassessment of the technical reproducibility of the assay for SNPs recall and genotyping precision. Our results highlight the value of a detailed probe screening on a preliminary genome assembly to produce reliable data for downstream genetic studies. This work should inspire and assist the development of similar genomic resources for other orphan crops and forest trees with highly heterozygous genomes.

Genome assembly of the Pink Ipê (Handroanthus impetiginosus, Bignoniaceae), a highly valued, ecologically keystone Neotropical timber forest tree.

Background: Handroanthus impetiginosus (Mart. ex DC.) Mattos is a keystone Neotropical hardwood tree widely distributed in seasonally dry tropical forests of South and Mesoamerica. Regarded as the "new mahogany," it is the second most expensive timber, the most logged species in Brazil, and currently under significant illegal trading pressure. The plant produces large amounts of quinoids, specialized metabolites with documented antitumorous and antibiotic effects. The development of genomic resources is needed to better understand and conserve the diversity of the species, to empower forensic identification of the origin of timber, and to identify genes for important metabolic compounds. Findings: The genome assembly covers 503.7 Mb (N50 = 81 316 bp), 90.4% of the 557-Mbp genome, with 13 206 scaffolds. A repeat database with 1508 sequences was developed, allowing masking of ∼31% of the assembly. Depth of coverage indicated that consensus determination adequately removed haplotypes assembled separately due to the extensive heterozygosity of the species. Automatic gene prediction provided 31 688 structures and 35 479 messenger RNA transcripts, while external evidence supported a well-curated set of 28 603 high-confidence models (90% of total). Finally, we used the genomic sequence and the comprehensive gene content annotation to identify genes related to the production of specialized metabolites. Conclusions: This genome assembly is the first well-curated resource for a Neotropical forest tree and the first one for a member of the Bignoniaceae family, opening exceptional opportunities to empower molecular, phytochemical, and breeding studies. This work should inspire the development of similar genomic resources for the largely neglected forest trees of the mega-diverse tropical biomes.

Regional heritability mapping and genome-wide association identify loci for complex growth, wood and disease resistance traits in Eucalyptus.

Although genome-wide association studies (GWAS) have provided valuable insights into the decoding of the relationships between sequence variation and complex phenotypes, they have explained little heritability. Regional heritability mapping (RHM) provides heritability estimates for genomic segments containing both common and rare allelic effects that individually contribute too little variance to be detected by GWAS. We carried out GWAS and RHM for seven growth, wood and disease resistance traits in a breeding population of 768 Eucalyptus hybrid trees using EuCHIP60K. Total genomic heritabilities accounted for large proportions (64-89%) of pedigree-based trait heritabilities, providing additional evidence that complex traits in eucalypts are controlled by many sequence variants across the frequency spectrum, each with small contributions to the phenotypic variance. RHM detected 26 quantitative trait loci (QTLs) encompassing 2191 single nucleotide polymorphisms (SNPs), whereas GWAS detected 13 single SNP-trait associations. RHM and GWAS QTLs individually explained 5-15% and 4-6% of the genomic heritability, respectively. RHM was superior to GWAS in capturing larger proportions of genomic heritability. Equated to previously mapped QTLs, our results highlighted genomic regions for further examination towards gene discovery. RHM-QTLs bearing a combination of common and rare variants could be useful enhancements to incorporate prior knowledge of the underlying genetic architecture in genomic prediction models.

Sugarcane giant borer transcriptome analysis and identification of genes related to digestion.

Sugarcane is a widely cultivated plant that serves primarily as a source of sugar and ethanol. Its annual yield can be significantly reduced by the action of several insect pests including the sugarcane giant borer (Telchin licus licus), a lepidopteran that presents a long life cycle and which efforts to control it using pesticides have been inefficient. Although its economical relevance, only a few DNA sequences are available for this species in the GenBank. Pyrosequencing technology was used to investigate the transcriptome of several developmental stages of the insect. To maximize transcript diversity, a pool of total RNA was extracted from whole body insects and used to construct a normalized cDNA database. Sequencing produced over 650,000 reads, which were de novo assembled to generate a reference library of 23,824 contigs. After quality score and annotation, 43% of the contigs had at least one BLAST hit against the NCBI non-redundant database, and 40% showed similarities with the lepidopteran Bombyx mori. In a further analysis, we conducted a comparison with Manduca sexta midgut sequences to identify transcripts of genes involved in digestion. Of these transcripts, many presented an expansion or depletion in gene number, compared to B. mori genome. From the sugarcane giant borer (SGB) transcriptome, a number of aminopeptidase N (APN) cDNAs were characterized based on homology to those reported as Cry toxin receptors. This is the first report that provides a large-scale EST database for the species. Transcriptome analysis will certainly be useful to identify novel developmental genes, to better understand the insect's biology and to guide the development of new strategies for insect-pest control.

Transcriptome analysis in cotton boll weevil (Anthonomus grandis) and RNA interference in insect pests.

Cotton plants are subjected to the attack of several insect pests. In Brazil, the cotton boll weevil, Anthonomus grandis, is the most important cotton pest. The use of insecticidal proteins and gene silencing by interference RNA (RNAi) as techniques for insect control are promising strategies, which has been applied in the last few years. For this insect, there are not much available molecular information on databases. Using 454-pyrosequencing methodology, the transcriptome of all developmental stages of the insect pest, A. grandis, was analyzed. The A. grandis transcriptome analysis resulted in more than 500.000 reads and a data set of high quality 20,841 contigs. After sequence assembly and annotation, around 10,600 contigs had at least one BLAST hit against NCBI non-redundant protein database and 65.7% was similar to Tribolium castaneum sequences. A comparison of A. grandis, Drosophila melanogaster and Bombyx mori protein families' data showed higher similarity to dipteran than to lepidopteran sequences. Several contigs of genes encoding proteins involved in RNAi mechanism were found. PAZ Domains sequences extracted from the transcriptome showed high similarity and conservation for the most important functional and structural motifs when compared to PAZ Domains from 5 species. Two SID-like contigs were phylogenetically analyzed and grouped with T. castaneum SID-like proteins. No RdRP gene was found. A contig matching chitin synthase 1 was mined from the transcriptome. dsRNA microinjection of a chitin synthase gene to A. grandis female adults resulted in normal oviposition of unviable eggs and malformed alive larvae that were unable to develop in artificial diet. This is the first study that characterizes the transcriptome of the coleopteran, A. grandis. A new and representative transcriptome database for this insect pest is now available. All data support the state of the art of RNAi mechanism in insects.