A Structural In Silico Analysis of the Immunogenicity of L-Asparaginase from Penicillium cerradense.

L-asparaginase is an essential drug used to treat acute lymphoid leukemia (ALL), a cancer of high prevalence in children. Several adverse reactions associated with L-asparaginase have been observed, mainly caused by immunogenicity and allergenicity. Some strategies have been adopted, such as searching for new microorganisms that produce the enzyme and applying protein engineering. Therefore, this work aimed to elucidate the molecular structure and predict the immunogenic profile of L-asparaginase from Penicillium cerradense, recently revealed as a new fungus of the genus Penicillium and producer of the enzyme, as a motivation to search for alternatives to bacterial L-asparaginase. In the evolutionary relationship, L-asparaginase from P. cerradense closely matches Aspergillus species. Using in silico tools, we characterized the enzyme as a protein fragment of 378 amino acids (39 kDa), including a signal peptide containing 17 amino acids, and the isoelectric point at 5.13. The oligomeric state was predicted to be a homotetramer. Also, this L-asparaginase presented a similar immunogenicity response (T- and B-cell epitopes) compared to Escherichia coli and Dickeya chrysanthemi enzymes. These results suggest a potentially useful L-asparaginase, with insights that can drive strategies to improve enzyme production.

Scalable and versatile container-based pipelines for de novo genome assembly and bacterial annotation.

Background: Advancements in DNA sequencing technology have transformed the field of bacterial genomics, allowing for faster and more cost effective chromosome level assemblies compared to a decade ago. However, transforming raw reads into a complete genome model is a significant computational challenge due to the varying quality and quantity of data obtained from different sequencing instruments, as well as intrinsic characteristics of the genome and desired analyses. To address this issue, we have developed a set of container-based pipelines using Nextflow, offering both common workflows for inexperienced users and high levels of customization for experienced ones. Their processing strategies are adaptable based on the sequencing data type, and their modularity enables the incorporation of new components to address the community's evolving needs. Methods: These pipelines consist of three parts: quality control, de novo genome assembly, and bacterial genome annotation. In particular, the genome annotation pipeline provides a comprehensive overview of the genome, including standard gene prediction and functional inference, as well as predictions relevant to clinical applications such as virulence and resistance gene annotation, secondary metabolite detection, prophage and plasmid prediction, and more. Results: The annotation results are presented in reports, genome browsers, and a web-based application that enables users to explore and interact with the genome annotation results. Conclusions: Overall, our user-friendly pipelines offer a seamless integration of computational tools to facilitate routine bacterial genomics research. The effectiveness of these is illustrated by examining the sequencing data of a clinical sample of Klebsiella pneumoniae.

Complete Genome Sequence of Pantoea stewartii RON18713 from Brazil Nut Tree Phyllosphere Reveals Genes Involved in Plant Growth Promotion.

The Amazonian rainforest is a hyper-diverse ecosystem in the number of species and the myriad of intertaxon relationships that are mostly understudied. In order to characterize a dominant and economically important Amazonian species, the Brazil nut tree (Bertholletia excelsa Bonpl.), at the genome level, wegenerated high-coverage long-read sequencing data from the leaves of a single individual. The genome assembly revealed an unexpected discovery: two circular contigs that could be assigned to the chromosome and a plasmid of a Pantoea stewartii strain. Comparative genomics revealed that this strain belongs to the indologenes subspecies and displays high synteny with other strains isolated from diseased leaves of the neotropical palm Bactris gasipaes Kunth. Investigation of pathogenicity-related genes revealed the absence of the entire type III secretion system gene cluster in the plasmid, which was otherwise highly similar to a plasmid from an isolate known to cause disease in Dracaena sanderiana Mast. In contrast, several genes associated with plant-growth promoting traits were detected, including genes involved in indole-3-acetic acid (IAA) production, phosphate solubilization, and biosynthesis of siderophores. In summary, we report the genome of an uncultivated P. stewartii subsp. indologenes strain associated with the Brazil nut tree and potentially a plant growth-promoting bacteria.

Genomic and physiological characterization of Novosphingobium terrae sp. nov., an alphaproteobacterium isolated from Cerrado soil containing a mega-sized chromid.

A novel bacterial strain, designated GeG2T, was isolated from soils of the native Cerrado, a highly biodiverse savanna-like Brazilian biome. 16S rRNA gene analysis of GeG2T revealed high sequence identity (100%) to the alphaproteobacterium Novosphingobium rosa; however, comparisons with N. rosa DSM 7285T showed several distinctive features, prompting a full characterization of the new strain in terms of physiology, morphology, and, ultimately, its genome. GeG2T cells were Gram-stain-negative bacilli, facultatively anaerobic, motile, positive for catalase and oxidase activities, and starch hydrolysis. Strain GeG2T presented planktonic-sessile dimorphism and cell aggregates surrounded by extracellular matrix and nanometric spherical structures were observed, suggesting the production of exopolysaccharides (EPS) and outer membrane vesicles (OMVs). Despite high 16S rDNA identity, strain GeG2T showed 90.38% average nucleotide identity and 42.60% digital DNA-DNA hybridization identity with N. rosa, below species threshold. Whole-genome assembly revealed four circular replicons: a 4.1 Mb chromosome, a 2.7 Mb extrachromosomal megareplicon, and two plasmids (212.7 and 68.6 kb). The megareplicon contains a few core genes and plasmid-type replication/maintenance systems, consistent with its classification as a chromid. Genome annotation shows a vast repertoire of carbohydrate-active enzymes and genes involved in the degradation of aromatic compounds, highlighting the biotechnological potential of the new isolate. Chemotaxonomic features, including polar lipid and fatty acid profiles, as well as physiological, molecular, and whole-genome comparisons showed significant differences between strain GeG2T and N. rosa, indicating that it represents a novel species, for which the name Novosphingobium terrae is proposed. The type strain is GeG2T (= CBMAI 2313T = CBAS 753 T).

A Novel Multidrug Resistant, Non-Tn4401 Genetic Element-Bearing, Strain of Klebsiella pneumoniae Isolated From an Urban Lake With Drinking and Recreational Water Reuse.

Antimicrobial resistance (AMR) is an increasing and urgent issue for human health worldwide, as it leads to the reduction of available antibiotics to treat bacterial infections, in turn increasing hospital stays and lethality. Therefore, the study and genomic surveillance of bacterial carriers of resistance in and outside of clinical settings is of utter importance. A colony of multidrug resistant (MDR) bacteria identified as Klebsiella spp., by 16S rDNA amplicon sequencing, has been isolated from an urban lake in Brazil, during a drug-degrading bacterial prospection. Genomic analyses revealed the bacteria as Klebsiella pneumoniae species. Furthermore, the in silico Multilocus Sequence Typing (MLST) identified the genome as a new sequence type, ST5236. The search for antimicrobial resistance genes (ARGs) detected the presence of genes against beta-lactams, fosfomycin, acriflavine and efflux pumps, as well as genes for heavy metal resistance. Of particular note, an extended-spectrum beta-lactamase gene (blaCTX-M-15) has been detected in close proximity to siphoviridae genes, while a carbapenemase gene (KPC-2) has been found in an extrachromosomal contig, within a novel non-Tn4401 genetic element (NTEKPC). An extrachromosomal contig found in the V3 isolate is identical to a contig of a K. pneumoniae isolate from a nearby hospital, which indicates a putative gene flow from the hospital network into Paranoá lake. The discovery of a MDR isolate in this lake is worrisome, as the region has recently undergone periods of water scarcity causing the lake, which receives treated wastewater effluent, and is already used for recreational purposes, to be used as an environmental buffer for drinking water reuse. Altogether, our results indicate an underrepresentation of environmental K. pneumoniae among available genomes, which may hamper the understanding of the population dynamics of the species in the environment and its consequences in the spread of ARGs and virulence genes.

Multidrug-Resistant (MDR) Klebsiella variicola Strains Isolated in a Brazilian Hospital Belong to New Clones.

Klebsiella variicola is mainly associated with opportunistic infections and frequently identified as Klebsiella pneumoniae. This misidentification implies a wrong epidemiology result as well as incorrect attribution to K. pneumoniae as the etiology of some severe infections. Recently, huge efforts have been made to study K. variicola, however, the biological aspects of this species are still unclear. Here we characterized five K. variicola strains initially identified as K. pneumoniae, with a Vitek-2 System and 16S rRNA sequencing. One-step multiplex polymerase chain reaction and Whole Genome Sequencing (WGS) identified them as K. variicola. Additionally, WGS analysis showed that all the strains are closely related with K. variicola genomes, forming a clustered group, apart from K. pneumoniae and K. quasipneumoniae. Multilocus sequence typing analysis showed four different sequence types (STs) among the strains and for two of them (Kv97 and Kv104) the same ST was assigned. All strains were multidrug-resistant (MDR) and three showed virulence phenotypes including invasion capacity to epithelial cells, and survival in human blood and serum. These results showed the emergence of new K. variicola clones with pathogenic potential to colonize and cause infection in different tissues. These characteristics associated with MDR strains raise great concern for human health.

A cost-effective approach to DNA methylation detection by Methyl Sensitive DArT sequencing.

Several studies suggest the relation of DNA methylation to diseases in humans and important phenotypes in plants drawing attention to this epigenetic mark as an important source of variability. In the last decades, several methodologies were developed to assess the methylation state of a genome. However, there is still a lack of affordable and precise methods for genome wide analysis in large sample size studies. Methyl sensitive double digestion MS-DArT sequencing method emerges as a promising alternative for methylation profiling. We developed a computational pipeline for the identification of DNA methylation using MS-DArT-seq data and carried out a pilot study using the Eucalyptus grandis tree sequenced for the species reference genome. Using a statistic framework as in differential expression analysis, 72,515 genomic sites were investigated and 5,846 methylated sites identified, several tissue specific, distributed along the species 11 chromosomes. We highlight a bias towards identification of DNA methylation in genic regions and the identification of 2,783 genes and 842 transposons containing methylated sites. Comparison with WGBS, DNA sequencing after treatment with bisulfite, data demonstrated a precision rate higher than 95% for our approach. The availability of a reference genome is useful for determining the genomic context of methylated sites but not imperative, making this approach suitable for any species. Our approach provides a cost effective, broad and reliable examination of DNA methylation profile on MspI/HpaII restriction sites, is fully reproducible and the source code is available on GitHub (https://github.com/wendelljpereira/ms-dart-seq).

Exploring the Brazilian diversity of Aspergillus sp. strains for lovastatin and itaconic acid production.

Filamentous fungi are well known for producing secondary metabolites applied in various industrial segments. Among these, lovastatin and itaconic acid, produced by Aspergillus terreus, have applications in the pharmaceutical and chemical industries. Lovastatin is primarily used for the control of hypercholesterolemia, while itaconic acid is a building block for the production of synthetic fibers, coating adhesives, among others. In this study, for the first time, 35 strains of Aspergillus sp. from four Brazilian culture collections were evaluated for lovastatin and itaconic acid production and compared to a reference strain, ATCC 20542. From an initial screening, the strains ATCC 20542, URM 224, URM1876, URM 5061, URM 5254, URM 5256, URM 5650, and URM 5961 were selected for genomic comparison. Among tested strains, the locus corresponding to the lovastatin genomic cluster was assembled, showing that all genes essential for lovastatin biosynthesis were present in producing URM 5961 and URM 5650 strains, with 100% and 98.5% similarity to ATCC 20542, respectively. However, in the no producing URM 1876, URM 224, URM 5254, URM 5061, and URM 5256 strains, this cluster was either fragmented or missing. Among the 35 strains evaluated for itaconic acid production in this study, only three strains had titers above 0.5 g/L, 16 strains had production below 0.5 g/L, and the remaining 18 strains had no production, with the highest production of itaconic acid observed in the URM 5254 strain with 2.2 g/L. The essential genes for itaconic acid production, mttA, cadA msfA were also mapped, where all three genes linked to itaconic acid production were found in a single contig in the assembly of each strain. In contrast to lovastatin loci, there is no correlation between the level of itaconic acid production and genetic polymorphisms in the genes associated with its biosynthesis.

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.

A Wor1-Like Transcription Factor Is Essential for Virulence of Cryptococcus neoformans.

Gti1/Pac2 transcription factors occur exclusively in fungi and their roles vary according to species, including regulating morphological transition and virulence, mating and secondary metabolism. Many of these functions are important for fungal pathogenesis. We therefore hypothesized that one of the two proteins of this family in Cryptococcus neoformans, a major pathogen of humans, would also control virulence-associated cellular processes. Elimination of this protein in C. neoformans results in reduced polysaccharide capsule expression and defective cytokinesis and growth at 37°C. The mutant loses virulence in a mouse model of cryptococcal infection and retains only partial virulence in the Galleria mellonella alternative model at 30°C. We performed RNA-Seq experiments on the mutant and found abolished transcription of genes that, in combination, are known to account for all the observed phenotypes. The protein has been named Required for cytokinesis and virulence 1 (Rcv1).

Draft Genome Sequence of the Antimicrobial-Producing Strain Paenibacillus elgii AC13.

A Paenibacillus elgii strain isolated from soil samples from Cerrado, Brazil, showed antimicrobial activity. Its genome sequence was acquired (GS20 FLX Titanium 454 platform) and comprises 108 contigs (N50, 198,427 bp) and 6,810 predicted sequences. Here, we shed some light on the antimicrobial genes of the strain, including a nonribosomal peptide synthetase (NRPS) module identified as part of a pelgipeptin gene cluster.

Transcriptome analyses of the Dof-like gene family in grapevine reveal its involvement in berry, flower and seed development.

The Dof (DNA-binding with one finger) protein family spans a group of plant transcription factors involved in the regulation of several functions, such as plant responses to stress, hormones and light, phytochrome signaling and seed germination. Here we describe the Dof-like gene family in grapevine (Vitis vinifera L.), which consists of 25 genes coding for Dof. An extensive in silico characterization of the VviDofL gene family was performed. Additionally, the expression of the entire gene family was assessed in 54 grapevine tissues and organs using an integrated approach with microarray (cv Corvina) and real-time PCR (cv Pinot Noir) analyses. The phylogenetic analysis comparing grapevine sequences with those of Arabidopsis, tomato, poplar and already described Dof genes in other species allowed us to identify several duplicated genes. The diversification of grapevine DofL genes during evolution likely resulted in a broader range of biological roles. Furthermore, distinct expression patterns were identified between samples analyzed, corroborating such hypothesis. Our expression results indicate that several VviDofL genes perform their functional roles mainly during flower, berry and seed development, highlighting their importance for grapevine growth and production. The identification of similar expression profiles between both approaches strongly suggests that these genes have important regulatory roles that are evolutionally conserved between grapevine cvs Corvina and Pinot Noir.

SNP discovery in common bean by restriction-associated DNA (RAD) sequencing for genetic diversity and population structure analysis.

Researchers have made great advances into the development and application of genomic approaches for common beans, creating opportunities to driving more real and applicable strategies for sustainable management of the genetic resource towards plant breeding. This work provides useful polymorphic single-nucleotide polymorphisms (SNPs) for high-throughput common bean genotyping developed by RAD (restriction site-associated DNA) sequencing. The RAD tags were generated from DNA pooled from 12 common bean genotypes, including breeding lines of different gene pools and market classes. The aligned sequences identified 23,748 putative RAD-SNPs, of which 3357 were adequate for genotyping; 1032 RAD-SNPs with the highest ADT (assay design tool) score are presented in this article. The RAD-SNPs were structurally annotated in different coding (47.00 %) and non-coding (53.00 %) sequence components of genes. A subset of 384 RAD-SNPs with broad genome distribution was used to genotype a diverse panel of 95 common bean germplasms and revealed a successful amplification rate of 96.6 %, showing 73 % of polymorphic SNPs within the Andean group and 83 % in the Mesoamerican group. A slightly increased He (0.161, n = 21) value was estimated for the Andean gene pool, compared to the Mesoamerican group (0.156, n = 74). For the linkage disequilibrium (LD) analysis, from a group of 580 SNPs (289 RAD-SNPs and 291 BARC-SNPs) genotyped for the same set of genotypes, 70.2 % were in LD, decreasing to 0.10 %in the Andean group and 0.77 % in the Mesoamerican group. Haplotype patterns spanning 310 Mb of the genome (60 %) were characterized in samples from different origins. However, the haplotype frameworks were under-represented for the Andean (7.85 %) and Mesoamerican (5.55 %) gene pools separately. In conclusion, RAD sequencing allowed the discovery of hundreds of useful SNPs for broad genetic analysis of common bean germplasm. From now, this approach provides an excellent panel of molecular tools for whole genome analysis, allowing integrating and better exploring the common bean breeding practices.

Genome-wide discovery and validation of Eucalyptus small RNAs reveals variable patterns of conservation and diversity across species of Myrtaceae.

BACKGROUND: Micro RNAs are a class of small non coding RNAs of 20-24 nucleotides transcribed as single stranded precursors from MIR gene loci. Initially described as post-transcriptional regulators involved in development, two decades ago, miRNAs have been proven to regulate a wide range of processes in plants such as germination, morphology and responses to biotic and abiotic stress. Despite wide conservation in plants, a number of miRNAs are lineage specific. We describe the first genome wide survey of Eucalyptus miRNAs based on high throughput sequencing. RESULTS: In addition to discovering small RNA sequences, MIR loci were mapped onto the reference genome and interspecific variability investigated. Sequencing was carried out for the two most world widely planted species, E. grandis and E. globulus. To maximize discovery, E. grandis samples were from BRASUZ1, the same tree whose genome provided the reference sequence. Interspecific analysis reinforces the variability in small RNA repertoire even between closely related species. Characterization of Eucalyptus small RNA sequences showed 95 orthologous to conserved miRNAs and 193 novel miRNAs. In silico target prediction confirmed 163 novel miRNAs and degradome sequencing experimentally confirmed several hundred targets. Experimental evidence based on the exclusive expression of a set of small RNAs across 16 species within Myrtaceae further highlighted variable patterns of conservation and diversity of these regulatory elements. CONCLUSIONS: The description of miRNAs in Eucalyptus contributes to scientific knowledge of this vast genre, which is the most widely planted hardwood crop in the tropical and subtropical world, adding another important element to the annotation of Eucalyptus grandis reference genome.

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.

Transcription profile of soybean-root-knot nematode interaction reveals a key role of phythormones in the resistance reaction.

BACKGROUND: Root-knot nematodes (RKN- Meloidogyne genus) present extensive challenges to soybean crop. The soybean line (PI 595099) is known to be resistant against specific strains and races of nematode species, thus its differential gene expression analysis can lead to a comprehensive gene expression profiling in the incompatible soybean-RKN interaction. Even though many disease resistance genes have been studied, little has been reported about phytohormone crosstalk on modulation of ROS signaling during soybean-RKN interaction. RESULTS: Using 454 technology to explore the common aspects of resistance reaction during both parasitism and resistance phases it was verified that hormone, carbohydrate metabolism and stress related genes were consistently expressed at high levels in infected roots as compared to mock control. Most noteworthy genes include those encoding glycosyltransferases, peroxidases, auxin-responsive proteins and gibberellin-regulated genes. Our data analysis suggests the key role of glycosyltransferases, auxins and components of gibberellin signal transduction, biosynthesis and deactivation pathways in the resistance reaction and their participation in jasmonate signaling and redox homeostasis in mediating aspects of plant growth and responses to biotic stress. CONCLUSIONS: Based on this study we suggest a reasonable model regarding to the complex mechanisms of crosstalk between plant hormones, mainly gibberellins and auxins, which can be crucial to modulate the levels of ROS in the resistance reaction to nematode invasion. The model also includes recent findings concerning to the participation of DELLA-like proteins and ROS signaling controlling plant immune or stress responses. Furthermore, this study provides a dataset of potential candidate genes involved in both nematode parasitism and resistance, which can be tested further for their role in this biological process using functional genomics approaches.

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.

Development of expressed sequence tag and expressed sequence tag-simple sequence repeat marker resources for Musa acuminata.

BACKGROUND AND AIMS: Banana (Musa acuminata) is a crop contributing to global food security. Many varieties lack resistance to biotic stresses, due to sterility and narrow genetic background. The objective of this study was to develop an expressed sequence tag (EST) database of transcripts expressed during compatible and incompatible banana-Mycosphaerella fijiensis (Mf) interactions. Black leaf streak disease (BLSD), caused by Mf, is a destructive disease of banana. Microsatellite markers were developed as a resource for crop improvement. METHODOLOGY: cDNA libraries were constructed from in vitro-infected leaves from BLSD-resistant M. acuminata ssp. burmaniccoides Calcutta 4 (MAC4) and susceptible M. acuminata cv. Cavendish Grande Naine (MACV). Clones were 5'-end Sanger sequenced, ESTs assembled with TGICL and unigenes annotated using BLAST, Blast2GO and InterProScan. Mreps was used to screen for simple sequence repeats (SSRs), with markers evaluated for polymorphism using 20 diploid (AA) M. acuminata accessions contrasting in resistance to Mycosphaerella leaf spot diseases. PRINCIPAL RESULTS: A total of 9333 high-quality ESTs were obtained for MAC4 and 3964 for MACV, which assembled into 3995 unigenes. Of these, 2592 displayed homology to genes encoding proteins with known or putative function, and 266 to genes encoding proteins with unknown function. Gene ontology (GO) classification identified 543 GO terms, 2300 unigenes were assigned to EuKaryotic orthologous group categories and 312 mapped to Kyoto Encyclopedia of Genes and Genomes pathways. A total of 624 SSR loci were identified, with trinucleotide repeat motifs the most abundant in MAC4 (54.1 %) and MACV (57.6 %). Polymorphism across M. acuminata accessions was observed with 75 markers. Alleles per polymorphic locus ranged from 2 to 8, totalling 289. The polymorphism information content ranged from 0.08 to 0.81. CONCLUSIONS: This EST collection offers a resource for studying functional genes, including transcripts expressed in banana-Mf interactions. Markers are applicable for genetic mapping, diversity characterization and marker-assisted breeding.

Genomic characterization of DArT markers based on high-density linkage analysis and physical mapping to the Eucalyptus genome.

Diversity Arrays Technology (DArT) provides a robust, high throughput, cost-effective method to query thousands of sequence polymorphisms in a single assay. Despite the extensive use of this genotyping platform for numerous plant species, little is known regarding the sequence attributes and genome-wide distribution of DArT markers. We investigated the genomic properties of the 7,680 DArT marker probes of a Eucalyptus array, by sequencing them, constructing a high density linkage map and carrying out detailed physical mapping analyses to the Eucalyptus grandis reference genome. A consensus linkage map with 2,274 DArT markers anchored to 210 microsatellites and a framework map, with improved support for ordering, displayed extensive collinearity with the genome sequence. Only 1.4 Mbp of the 75 Mbp of still unplaced scaffold sequence was captured by 45 linkage mapped but physically unaligned markers to the 11 main Eucalyptus pseudochromosomes, providing compelling evidence for the quality and completeness of the current Eucalyptus genome assembly. A highly significant correspondence was found between the locations of DArT markers and predicted gene models, while most of the 89 DArT probes unaligned to the genome correspond to sequences likely absent in E. grandis, consistent with the pan-genomic feature of this multi-Eucalyptus species DArT array. These comprehensive linkage-to-physical mapping analyses provide novel data regarding the genomic attributes of DArT markers in plant genomes in general and for Eucalyptus in particular. DArT markers preferentially target the gene space and display a largely homogeneous distribution across the genome, thereby providing superb coverage for mapping and genome-wide applications in breeding and diversity studies. Data reported on these ubiquitous properties of DArT markers will be particularly valuable to researchers working on less-studied crop species who already count on DArT genotyping arrays but for which no reference genome is yet available to allow such detailed characterization.

Global transcriptome analysis of two wild relatives of peanut under drought and fungi infection.

BACKGROUND: Cultivated peanut (Arachis hypogaea) is one of the most widely grown grain legumes in the world, being valued for its high protein and unsaturated oil contents. Worldwide, the major constraints to peanut production are drought and fungal diseases. Wild Arachis species, which are exclusively South American in origin, have high genetic diversity and have been selected during evolution in a range of environments and biotic stresses, constituting a rich source of allele diversity. Arachis stenosperma harbors resistances to a number of pests, including fungal diseases, whilst A. duranensis has shown improved tolerance to water limited stress. In this study, these species were used for the creation of an extensive databank of wild Arachis transcripts under stress which will constitute a rich source for gene discovery and molecular markers development. RESULTS: Transcriptome analysis of cDNA collections from A. stenosperma challenged with Cercosporidium personatum (Berk. and M.A. Curtis) Deighton, and A. duranensis submitted to gradual water limited stress was conducted using 454 GS FLX Titanium generating a total of 7.4 x 10(5) raw sequence reads covering 211 Mbp of both genomes. High quality reads were assembled to 7,723 contigs for A. stenosperma and 12,792 for A. duranensis and functional annotation indicated that 95% of the contigs in both species could be appointed to GO annotation categories. A number of transcription factors families and defense related genes were identified in both species. Additionally, the expression of five A. stenosperma Resistance Gene Analogs (RGAs) and four retrotransposon (FIDEL-related) sequences were analyzed by qRT-PCR. This data set was used to design a total of 2,325 EST-SSRs, of which a subset of 584 amplified in both species and 214 were shown to be polymorphic using ePCR. CONCLUSIONS: This study comprises one of the largest unigene dataset for wild Arachis species and will help to elucidate genes involved in responses to biological processes such as fungal diseases and water limited stress. Moreover, it will also facilitate basic and applied research on the genetics of peanut through the development of new molecular markers and the study of adaptive variation across the genus.