The first draft genome of Picrorhiza kurrooa, an endangered medicinal herb from Himalayas.

Picrorhiza kurrooa is an endangered medicinal herb which is distributed across the Himalayan region at an altitude between 3000-5000 m above mean sea level. The medicinal properties of P. kurrooa are attributed to monoterpenoid picrosides present in leaf, rhizome and root of the plant. However, no genomic information is currently available for P. kurrooa, which limits our understanding about its molecular systems and associated responses. The present study brings the first assembled draft genome of P. kurrooa by using 227 Gb of raw data generated by Illumina and PacBio RS II sequencing platforms. The assembled genome has a size of n = ~ 1.7 Gb with 12,924 scaffolds. Four pronged assembly quality validations studies, including experimentally reported ESTs mapping and directed sequencing of the assembled contigs, confirmed high reliability of the assembly. About 76% of the genome is covered by complex repeats alone. Annotation revealed 24,798 protein coding and 9789 non-coding genes. Using the assembled genome, a total of 710 miRNAs were discovered, many of which were found responsible for molecular response against temperature changes. The miRNAs and targets were validated experimentally. The availability of draft genome sequence will aid in genetic improvement and conservation of P. kurrooa. Also, this study provided an efficient approach for assembling complex genomes while dealing with repeats when regular assemblers failed to progress due to repeats.

Chromosome-level genome assembly of Ophiorrhiza pumila reveals the evolution of camptothecin biosynthesis.

Plant genomes remain highly fragmented and are often characterized by hundreds to thousands of assembly gaps. Here, we report chromosome-level reference and phased genome assembly of Ophiorrhiza pumila, a camptothecin-producing medicinal plant, through an ordered multi-scaffolding and experimental validation approach. With 21 assembly gaps and a contig N50 of 18.49 Mb, Ophiorrhiza genome is one of the most complete plant genomes assembled to date. We also report 273 nitrogen-containing metabolites, including diverse monoterpene indole alkaloids (MIAs). A comparative genomics approach identifies strictosidine biogenesis as the origin of MIA evolution. The emergence of strictosidine biosynthesis-catalyzing enzymes precede downstream enzymes' evolution post γ whole-genome triplication, which occurred approximately 110 Mya in O. pumila, and before the whole-genome duplication in Camptotheca acuminata identified here. Combining comparative genome analysis, multi-omics analysis, and metabolic gene-cluster analysis, we propose a working model for MIA evolution, and a pangenome for MIA biosynthesis, which will help in establishing a sustainable supply of camptothecin.

Identification of ancient viruses from metagenomic data of the Jomon people.

Ancient DNA studies provide genomic information about the origins, population structures, and physical characteristics of ancient humans that cannot be solely examined by archeological studies. The DNAs extracted from ancient human bones, teeth, or tissues are often contaminated with coexisting bacterial and viral genomes that contain DNA from ancient microbes infecting those of ancient humans. Information on ancient viral genomes is useful in making inferences about the viral evolution. Here, we have utilized metagenomic sequencing data from the dental pulp of five Jomon individuals, who lived on the Japanese archipelago more than 3000 years ago; this is to detect ancient viral genomes. We conducted de novo assembly of the non-human reads where we have obtained 277,387 contigs that were longer than 1000 bp. These contigs were subjected to homology searches against a collection of modern viral genome sequences. We were able to detect eleven putative ancient viral genomes. Among them, we reconstructed the complete sequence of the Siphovirus contig89 (CT89) viral genome. The Jomon CT89-like sequence was determined to contain 59 open reading frames, among which five genes known to encode phage proteins were under strong purifying selection. The host of CT89 was predicted to be Schaalia meyeri, a bacterium residing in the human oral cavity. Finally, the CT89 phylogenetic tree showed two clusters, from both of which the Jomon sequence was separated. Our results suggest that metagenomic information from the dental pulp of the Jomon people is essential in retrieving ancient viral genomes used to examine their evolution.

A draft genome for Spatholobus suberectus.

Spatholobus suberectus Dunn (S. suberectus), which belongs to the Leguminosae, is an important medicinal plant in China. Owing to its long growth cycle and increased use in human medicine, wild resources of S. suberectus have decreased rapidly and may be on the verge of extinction. De novo assembly of the whole S. suberectus genome provides us a critical potential resource towards biosynthesis of the main bioactive components and seed development regulation mechanism of this plant. Utilizing several sequencing technologies such as Illumina HiSeq X Ten, single-molecule real-time sequencing, 10x Genomics, as well as new assembly techniques such as FALCON and chromatin interaction mapping (Hi-C), we assembled a chromosome-scale genome about 798 Mb in size. In total, 748 Mb (93.73%) of the contig sequences were anchored onto nine chromosomes with the longest scaffold being 103.57 Mb. Further annotation analyses predicted 31,634 protein-coding genes, of which 93.9% have been functionally annotated. All data generated in this study is available in public databases.

Molecular characterization and gene expression patterns of retinoid receptors, in normal and regenerating tissues of the sea cucumber, Holothuria glaberrima.

Retinoic acid receptors (RAR) and retinoid X receptors (RXR) are ligand-mediated transcription factors that synchronize intricate signaling networks in metazoans. Dimer formation between these two nuclear receptors mediates the recruitment of co-regulatory complexes coordinating the progression of signaling cascades during developmental and regenerative events. In the present study we identified and characterized the receptors for retinoic acid in the sea cucumber Holothuria glaberrima; a model system capable of regenerative organogenesis during adulthood. Molecular characterizations revealed the presence of three isoforms of RAR and two of RXR as a consequence of alternative splicing events. Various analyses including: primary structure sequencing, phylogenetic analysis, protein domain prediction, and multiple sequence alignment further confirmed their identity. Semiquantitative reverse transcription PCR analysis of each receptor isoform herein identified showed that the retinoid receptors are expressed in all tissues sampled: the mesenteries, respiratory trees, muscles, gonads, and the digestive tract. During regenerative organogenesis two of the receptors (RAR-L and RXR-T) showed differential expression in the posterior segment while RAR-S is differentially expressed in the anterior segment of the intestine. This work presents the first description of the components relaying the signaling for retinoic acid within this model system.

Crop wild relative populations of Beta vulgaris allow direct mapping of agronomically important genes.

Rapid identification of agronomically important genes is of pivotal interest for crop breeding. One source of such genes are crop wild relative (CWR) populations. Here we used a CWR population of <200 wild beets (B. vulgaris ssp. maritima), sampled in their natural habitat, to identify the sugar beet (Beta vulgaris ssp. vulgaris) resistance gene Rz2 with a modified version of mapping-by-sequencing (MBS). For that, we generated a draft genome sequence of the wild beet. Our results show the importance of preserving CWR in situ and demonstrate the great potential of CWR for rapid discovery of causal genes relevant for crop improvement. The candidate gene for Rz2 was identified by MBS and subsequently corroborated via RNA interference (RNAi). Rz2 encodes a CC-NB-LRR protein. Access to the DNA sequence of Rz2 opens the path to improvement of resistance towards rhizomania not only by marker-assisted breeding but also by genome editing.

Hybrid de novo genome assembly of the Chinese herbal fleabane Erigeron breviscapus.

Background: The plants in the Erigeron genus of the Compositae (Asteraceae) family are commonly called fleabanes, possibly due to the belief that certain chemicals in these plants repel fleas. In the traditional Chinese medicine, Erigeron breviscapus , which is native to China, was widely used in the treatment of cerebrovascular disease. A handful of bioactive compounds, including scutellarin, 3,5-dicaffeoylquinic acid, and 3,4-dicaffeoylquinic acid, have been isolated from the plant. With the purpose of finding novel medicinal compounds and understanding their biosynthetic pathways, we propose to sequence the genome of E. breviscapus . We assembled the highly heterozygous E. breviscapus genome using a combination of PacBio single-molecular real-time sequencing and next-generation sequencing methods on the Illumina HiSeq platform. The final draft genome is approximately 1.2 Gb, with contig and scaffold N50 sizes of 18.8 kb and 31.5 kb, respectively. Further analyses predicted 37 504 protein-coding genes in the E. breviscapus genome and 8172 shared gene families among Compositae species. The E. breviscapus genome provides a valuable resource for the investigation of novel bioactive compounds in this Chinese herb.

Draft genome of the lined seahorse, Hippocampus erectus.

Background: The lined seahorse, Hippocampus erectus , is an Atlantic species and mainly inhabits shallow sea beds or coral reefs. It has become very popular in China for its wide use in traditional Chinese medicine. In order to improve the aquaculture yield of this valuable fish species, we are trying to develop genomic resources for assistant selection in genetic breeding. Here, we provide whole genome sequencing, assembly, and gene annotation of the lined seahorse, which can enrich genome resource and further application for its molecular breeding. A total of 174.6 Gb (Gigabase) raw DNA sequences were generated by the Illumina Hiseq2500 platform. The final assembly of the lined seahorse genome is around 458 Mb, representing 94% of the estimated genome size (489 Mb by k-mer analysis). The contig N50 and scaffold N50 reached 14.57 kb and 1.97 Mb, respectively. Quality of the assembled genome was assessed by BUSCO with prediction of 85% of the known vertebrate genes and evaluated using the de novo assembled RNA-seq transcripts to prove a high mapping ratio (more than 99% transcripts could be mapped to the assembly). Using homology-based, de novo and transcriptome-based prediction methods, we predicted 20 788 protein-coding genes in the generated assembly, which is less than our previously reported gene number (23 458) of the tiger tail seahorse ( H. comes ). We report a draft genome of the lined seahorse. These generated genomic data are going to enrich genome resource of this economically important fish, and also provide insights into the genetic mechanisms of its iconic morphology and male pregnancy behavior.

From Pine Cones to Read Clouds: Rescaffolding the Megagenome of Sugar Pine (Pinus lambertiana).

We investigate the utility and scalability of new read cloud technologies to improve the draft genome assemblies of the colossal, and largely repetitive, genomes of conifers. Synthetic long read technologies have existed in various forms as a means of reducing complexity and resolving repeats since the outset of genome assembly. Recently, technologies that combine subhaploid pools of high molecular weight DNA with barcoding on a massive scale have brought new efficiencies to sample preparation and data generation. When combined with inexpensive light shotgun sequencing, the resulting data can be used to scaffold large genomes. The protocol is efficient enough to consider routinely for even the largest genomes. Conifers represent the largest reference genome projects executed to date. The largest of these is that of the conifer Pinus lambertiana (sugar pine), with a genome size of 31 billion bp. In this paper, we report on the molecular and computational protocols for scaffolding the P. lambertiana genome using the library technology from 10× Genomics. At 247,000 bp, the NG50 of the existing reference sequence is the highest scaffold contiguity among the currently published conifer assemblies; this new assembly's NG50 is 1.94 million bp, an eightfold increase.

Transcriptome characterization and gene expression analysis related to sexual dimorphism in the ghost moth, Thitarodes pui, a host of Ophiocordyceps sinensis.

Thitarodes pui is one of the host species of the Chinese caterpillar fungus Ophiocordyceps sinensis as a traditional Chinese medicine with economic and medical importance. The pupal and adult stages of T. pui are sexually dimorphic. In order to elucidate the molecular mechanisms involved in the sexually dimorphic development of T. pui, we compared the transcriptomes of female and male pupae and adults. We obtained 15,881,734, 16,962,086, 17,514,743, and 17,770,904 clean reads from female pupae, male pupae, female adults, and male adults, respectively. The reads obtained from the four samples were pooled and assembled into 65,165 unigenes, 23,597 of which were annotated. Candidate genes involved in sexual development were identified and analysed. Gene expression analysis revealed that 1406 genes were differentially expressed in male and female pupae, 448 of which were up-regulated in males and 958 were up-regulated in females. A total of 2025 genes were differentially expressed in male and females adults, 1304 of which were up-regulated in males and 721 were up-regulated in females. The functional enrichment of the differentially expressed genes indicated that reproduction and cuticle synthesis were regulated differently between the sexes. The transcriptome data obtained provide significant information regarding the genes involved in sexually dimorphic development, which will improve our understanding of the molecular mechanisms related to sexual dimorphism and helpful for the moth mass rearing which would provide enough host insects for the sustainable utilization of O. sinensis.

Assembly of the draft genome of buckwheat and its applications in identifying agronomically useful genes.

Buckwheat (Fagopyrum esculentum Moench; 2n = 2x = 16) is a nutritionally dense annual crop widely grown in temperate zones. To accelerate molecular breeding programmes of this important crop, we generated a draft assembly of the buckwheat genome using short reads obtained by next-generation sequencing (NGS), and constructed the Buckwheat Genome DataBase. After assembling short reads, we determined 387,594 scaffolds as the draft genome sequence (FES_r1.0). The total length of FES_r1.0 was 1,177,687,305 bp, and the N50 of the scaffolds was 25,109 bp. Gene prediction analysis revealed 286,768 coding sequences (CDSs; FES_r1.0_cds) including those related to transposable elements. The total length of FES_r1.0_cds was 212,917,911 bp, and the N50 was 1,101 bp. Of these, the functions of 35,816 CDSs excluding those for transposable elements were annotated by BLAST analysis. To demonstrate the utility of the database, we conducted several test analyses using BLAST and keyword searches. Furthermore, we used the draft genome as a reference sequence for NGS-based markers, and successfully identified novel candidate genes controlling heteromorphic self-incompatibility of buckwheat. The database and draft genome sequence provide a valuable resource that can be used in efforts to develop buckwheat cultivars with superior agronomic traits.

Transfer of energy pathway genes in microbial enhanced biological phosphorus removal communities.

BACKGROUND: Lateral gene transfer (LGT) is an important evolutionary process in microbial evolution. In sewage treatment plants, LGT of antibiotic resistance and xenobiotic degradation-related proteins has been suggested, but the role of LGT outside these processes is unknown. Microbial communities involved in Enhanced Biological Phosphorus Removal (EBPR) have been used to treat wastewater in the last 50 years and may provide insights into adaptation to an engineered environment. We introduce two different types of analysis to identify LGT in EBPR sewage communities, based on identifying assembled sequences with more than one strong taxonomic match, and on unusual phylogenetic patterns. We applied these methods to investigate the role of LGT in six energy-related metabolic pathways. RESULTS: The analyses identified overlapping but non-identical sets of transferred enzymes. All of these were homologous with sequences from known mobile genetic elements, and many were also in close proximity to transposases and integrases in the EBPR data set. The taxonomic method had higher sensitivity than the phylogenetic method, identifying more potential LGTs. Both analyses identified the putative transfer of five enzymes within an Australian community, two in a Danish community, and none in a US-derived culture. CONCLUSIONS: Our methods were able to identify sequences with unusual phylogenetic or compositional properties as candidate LGT events. The association of these candidates with known mobile elements supports the hypothesis of transfer. The results of our analysis strongly suggest that LGT has influenced the development of functionally important energy-related pathways in EBPR systems, but transfers may be unique to each community due to different operating conditions or taxonomic composition.

Identification of candidate genes associated with fertility restoration of cytoplasmic male-sterility in onion (Allium cepa L.) using a combination of bulked segregant analysis and RNA-seq.

KEY MESSAGE: A combination of BSA and RNA-seq was performed to identify candidates for the restorer-of-fertility gene in onion. The AcPMS1 involved in DNA mismatch repair was identified as the best candidate. To identify candidate genes of the restorer-of-fertility gene (Ms) responsible for fertility restoration of onion cytoplasmic male-sterility, a combined approach of bulked segregant analysis and RNA-seq was employed. From 32,674 de novo assembled contigs, 430 perfectly homozygous SNPs between male-fertile (MF) and male-sterile (MS) bulks were identified in 141 contigs. After verifying the homozygosity of the SNPs by PCR amplification and sequencing, the SNPs on 139 of the contigs were genotypes for the two recombinants which contained crossover events between the Ms locus and two tightly linked molecular markers. As a result, 30 contigs showing perfect linkage with the Ms locus in the large-sized segregating population were identified. Among them, 14 showed perfect linkage disequilibrium (LD) with the Ms locus, as determined by genotyping 251 domestic breeding lines. Furthermore, molecular markers tagging the 14 contigs also showed almost perfect LD with each other in 124 exotic accessions introduced from 21 countries, except for one accession which contained a crossover event by which the 14 markers were divided into two groups. After sequencing of the full-length cDNA of the 14 contigs showing perfect LD, the deduced amino acids sequences of the MF and MS alleles were compared. Four genes were shown to harbor putative critical amino acid changes in the known domains. Among them, the gene encoding PMS1, involved in the DNA mismatch repair pathway, was assumed to be the best candidate gene responsible for fertility restoration of male-sterility in onion.

Genome reconstruction and gene expression of "Candidatus Accumulibacter phosphatis" Clade IB performing biological phosphorus removal.

We report the first integrated metatranscriptomic and metagenomic analysis of enhanced biological phosphorus removal (EBPR) sludge. A draft genome of Candidatus Accumulibacter spp. strain HKU-1, a member of Clade IB, was retrieved. It was estimated to be ∼90% complete and shared average nucleotide identities of 83% and 88% with the finished genome CAP IIA UW-1 and the draft genome CAP IA UW-2, respectively. Different from CAP IIA UW-1, the phosphotransferase (pap) in polyphosphate metabolism and V-ATPase in orthophosphate transport were absent from CAP IB HKU-1. Additionally, unlike CAP IA UW-2, CAP IB HKU-1 carried the genes for carbon fixation and nitrogen fixation. Despite these differences, the key genes required for acetate uptake, glycolysis and polyhydroxyalkanoate (PHA) synthesis were conserved in all these Accumulibacter genomes. The preliminary metatranscriptomic results revealed that the most significantly up-regulated genes of CAP IB HKU-1 from the anaerobic to the aerobic phase were responsible for assimilatory sulfate reduction, genetic information processing and phosphorus absorption, while the down-regulated genes were related to N2O reduction, PHA synthesis and acetyl-CoA formation. This study yielded another important Accumulibacter genome, revealed the functional difference within the Accumulibacter Type I, and uncovered the genetic responses to EBPR stimuli at a higher resolution.

Transcriptome analysis of Catharanthus roseus for gene discovery and expression profiling.

The medicinal plant, Catharanthus roseus, accumulates wide range of terpenoid indole alkaloids, which are well documented therapeutic agents. In this study, deep transcriptome sequencing of C. roseus was carried out to identify the pathways and enzymes (genes) involved in biosynthesis of these compounds. About 343 million reads were generated from different tissues (leaf, flower and root) of C. roseus using Illumina platform. Optimization of de novo assembly involving a two-step process resulted in a total of 59,220 unique transcripts with an average length of 1284 bp. Comprehensive functional annotation and gene ontology (GO) analysis revealed the representation of many genes involved in different biological processes and molecular functions. In total, 65% of C. roseus transcripts showed homology with sequences available in various public repositories, while remaining 35% unigenes may be considered as C. roseus specific. In silico analysis revealed presence of 11,620 genic simple sequence repeats (excluding mono-nucleotide repeats) and 1820 transcription factor encoding genes in C. roseus transcriptome. Expression analysis showed roots and leaves to be actively participating in bisindole alkaloid production with clear indication that enzymes involved in pathway of vindoline and vinblastine biosynthesis are restricted to aerial tissues. Such large-scale transcriptome study provides a rich source for understanding plant-specialized metabolism, and is expected to promote research towards production of plant-derived pharmaceuticals.

The large-scale investigation of gene expression in Leymus chinensis stigmas provides a valuable resource for understanding the mechanisms of poaceae self-incompatibility.

BACKGROUND: Many Poaceae species show a gametophytic self-incompatibility (GSI) system, which is controlled by at least two independent and multiallelic loci, S and Z. Until currently, the gene products for S and Z were unknown. Grass SI plant stigmas discriminate between pollen grains that land on its surface and support compatible pollen tube growth and penetration into the stigma, whereas recognizing incompatible pollen and thus inhibiting pollination behaviors. Leymus chinensis (Trin.) Tzvel. (sheepgrass) is a Poaceae SI species. A comprehensive analysis of sheepgrass stigma transcriptome may provide valuable information for understanding the mechanism of pollen-stigma interactions and grass SI. RESULTS: The transcript abundance profiles of mature stigmas, mature ovaries and leaves were examined using high-throughput next generation sequencing technology. A comparative transcriptomic analysis of these tissues identified 1,025 specifically or preferentially expressed genes in sheepgrass stigmas. These genes contained a significant proportion of genes predicted to function in cell-cell communication and signal transduction. We identified 111 putative transcription factors (TFs) genes and the most abundant groups were MYB, C2H2, C3H, FAR1, MADS. Comparative analysis of the sheepgrass, rice and Arabidopsis stigma-specific or preferential datasets showed broad similarities and some differences in the proportion of genes in the Gene Ontology (GO) functional categories. Potential SI candidate genes identified in other grasses were also detected in the sheepgrass stigma-specific or preferential dataset. Quantitative real-time PCR experiments validated the expression pattern of stigma preferential genes including homologous grass SI candidate genes. CONCLUSIONS: This study represents the first large-scale investigation of gene expression in the stigmas of an SI grass species. We uncovered many notable genes that are potentially involved in pollen-stigma interactions and SI mechanisms, including genes encoding receptor-like protein kinases (RLK), CBL (calcineurin B-like proteins) interacting protein kinases, calcium-dependent protein kinase, expansins, pectinesterase, peroxidases and various transcription factors. The availability of a pool of stigma-specific or preferential genes for L. chinensis offers an opportunity to elucidate the mechanisms of SI in Poaceae.

Antimicrobial peptides in the centipede Scolopendra subspinipes mutilans.

The centipede Scolopendra subspinipes mutilans is an environmentally beneficial and medically important arthropod species. Although this species is increasingly applied as a reliable source of new antimicrobial peptides, the transcriptome of this species is a prerequisite for more rational selection of antimicrobial peptides. In this report, we isolated total RNA from the whole body of adult centipedes, S. subspinipes mutilans, that were nonimmunized and immunized against Escherichia coli, and we generated a total of 77,063 pooled contigs and singletons using high-throughput sequencing. To screen putative antimicrobial peptides, in silico analyses of the S. subspinipes mutilans transcriptome were performed based on the physicochemical evidence of length, charge, isoelectric point, and in vitro and in vivo aggregation scores together with the existence of continuous antimicrobial peptide stretches. Moreover, we excluded some transcripts that showed similarity with both previously known antimicrobial peptides and the human proteome, had a proteolytic cleavage site, and had downregulated expression compared with the nonimmunized sample. As a result, we selected 17 transcripts and tested their antimicrobial activity with a radial diffusion assay. Among them, ten synthetic peptides experimentally showed antimicrobial activity against microbes and no toxicity to mouse erythrocytes. Our results provide not only a useful set of antimicrobial peptide candidates and an efficient strategy for novel antimicrobial peptide development but also the transcriptome data of a big centipede as a valuable resource.

Identification of phenylpropanoid biosynthetic genes and phenylpropanoid accumulation by transcriptome analysis of Lycium chinense.

BACKGROUND: Lycium chinense is well known in traditional Chinese herbal medicine for its medicinal value and composition, which have been widely studied for decades. However, further research on Lycium chinense is limited due to the lack of transcriptome and genomic information. RESULTS: The transcriptome of L. chinense was constructed by using an Illumina HiSeq 2000 sequencing platform. All 56,526 unigenes with an average length of 611 nt and an N50 equaling 848 nt were generated from 58,192,350 total raw reads after filtering and assembly. Unigenes were assembled by BLAST similarity searches and annotated with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) orthology identifiers. Using these transcriptome data, the majority of genes that are associated with phenylpropanoid biosynthesis in L. chinense were identified. In addition, phenylpropanoid biosynthesis-related gene expression and compound content in different organs were analyzed. We found that most phenylpropanoid genes were highly expressed in the red fruits, leaves, and flowers. An important phenylpropanoid, chlorogenic acid, was also found to be extremely abundant in leaves. CONCLUSIONS: Using Illumina sequencing technology, we have identified the function of novel homologous genes that regulate metabolic pathways in Lycium chinense.

Identification of expressed sequences in the coffee genome potentially associated with somatic embryogenesis.

Brazil possesses the most modern and productive coffee growing farms in the world, but technological development is desired to cope with the increasing world demand. One way to increase Brazilian coffee growing productivity is wide scale production of clones with superior genotypes, which can be obtained with in vitro propagation technique, or from tissue culture. These procedures can generate thousands of clones. However, the methodologies for in vitro cultivation are genotype-dependent, which leads to an almost empirical development of specific protocols for each species. Therefore, molecular markers linked to the biochemical events of somatic embryogenesis would greatly facilitate the development of such protocols. In this context, sequences potentially involved in embryogenesis processes in the coffee plant were identified in silico from libraries generated by the Brazilian Coffee Genome Project. Through these in silico analyses, we identified 15 EST-contigs related to the embryogenesis process. Among these, 5 EST-contigs (3605, 9850, 13686, 17240, and 17265) could readily be associated with plant embryogenesis. Sequence analysis of EST-contig 3605, 9850, and 17265 revealed similarity to a polygalacturonase, to a cysteine-proteinase, and to an allergenine, respectively. Results also show that EST-contig 17265 sequences presented similarity to an expansin. Finally, analysis of EST-contig 17240 revealed similarity to a protein of unknown function, but it grouped in the similarity dendrogram with the WUSCHEL transcription factor. The data suggest that these EST-contigs are related to the embryogenic process and have potential as molecular markers to increase methodological efficiency in obtaining coffee plant embryogenic materials.

Identification of genes and pathways related to phenol degradation in metagenomic libraries from petroleum refinery wastewater.

Two fosmid libraries, totaling 13,200 clones, were obtained from bioreactor sludge of petroleum refinery wastewater treatment system. The library screening based on PCR and biological activity assays revealed more than 400 positive clones for phenol degradation. From these, 100 clones were randomly selected for pyrosequencing in order to evaluate the genetic potential of the microorganisms present in wastewater treatment plant for biodegradation, focusing mainly on novel genes and pathways of phenol and aromatic compound degradation. The sequence analysis of selected clones yielded 129,635 reads at an estimated 17-fold coverage. The phylogenetic analysis showed Burkholderiales and Rhodocyclales as the most abundant orders among the selected fosmid clones. The MG-RAST analysis revealed a broad metabolic profile with important functions for wastewater treatment, including metabolism of aromatic compounds, nitrogen, sulphur and phosphorus. The predicted 2,276 proteins included phenol hydroxylases and cathecol 2,3- dioxygenases, involved in the catabolism of aromatic compounds, such as phenol, byphenol, benzoate and phenylpropanoid. The sequencing of one fosmid insert of 33 kb unraveled the gene that permitted the host, Escherichia coli EPI300, to grow in the presence of aromatic compounds. Additionally, the comparison of the whole fosmid sequence against bacterial genomes deposited in GenBank showed that about 90% of sequence showed no identity to known sequences of Proteobacteria deposited in the NCBI database. This study surveyed the functional potential of fosmid clones for aromatic compound degradation and contributed to our knowledge of the biodegradative capacity and pathways of microbial assemblages present in refinery wastewater treatment system.