Whole genome sequence comparisons and "full-length" cDNA sequences: a combined approach to evaluate and improve Arabidopsis genome annotation

Linking discoveries, mechanisms, and technologies to develop a clearer perspective on plant long noncoding RNAs

Long noncoding RNAs (lncRNAs) are a large and diverse class of genes in eukaryotic genomes that contribute to a variety of regulatory processes. Functionally characterized lncRNAs play critical roles in plants, ranging from regulating flowering to controlling lateral root formation. However, findings from the past decade have revealed that thousands of lncRNAs are present in plant transcriptomes, and characterization has lagged far behind identification. In this setting, distinguishing function from noise is challenging. However, the plant community has been at the forefront of discovery in lncRNA biology, providing many functional and mechanistic insights that have increased our understanding of this gene class. In this review, we examine the key discoveries and insights made in plant lncRNA biology over the past two and a half decades. We describe how discoveries made in the pregenomics era have informed efforts to identify and functionally characterize lncRNAs in the subsequent decades. We provide an overview of the functional archetypes into which characterized plant lncRNAs fit and speculate on new avenues of research that may uncover yet more archetypes. Finally, this review discusses the challenges facing the field and some exciting new molecular and computational approaches that may help inform lncRNA comparative and functional analyses.

Transcriptome profiling based on Illumina- and SMRT-based RNA-seq reveals circadian regulation of key pathways in flower bud development in walnut

Flower bud development is a defining feature of walnut, which contributes to the kernel yield, yield stability, fruit quality and commodity value. However, little is known about the mechanism of the flower bud development in walnut. Here, the stages of walnut female flower bud development were divided into five period (P01-05) by using histological observation. They were further studied through PacBio Iso-Seq and RNA-seq analysis. Accordingly, we obtained 52,875 full-length transcripts, where 4,579 were new transcripts, 3,065 were novel genes, 1,437 were consensus lncRNAs and 20,813 were alternatively spliced isoforms. These transcripts greatly improved the current genome annotation and enhanced our understanding of the walnut transcriptome. Next, RNA sequencing of female flower buds at five periods revealed that circadian rhythm-plant was commonly enriched along with the flower bud developmental gradient. A total of 14 differentially expressed genes (DEGs) were identified, and six of them were confirmed by real-time quantitative analysis. Additionally, six and two differentially expressed clock genes were detected to be regulated by AS events and lncRNAs, respectively. All these detected plant circadian genes form a complex interconnected network to regulate the flower bud development. Thus, investigation of key genes associated with the circadian clock could clarify the process of flower bud development in walnut.

The mutation nrpb1-A325V in the largest subunit of RNA polymerase II suppresses compromised growth of Arabidopsis plants deficient in a function of the general transcription factor IIF

The evolutionarily conserved 12-subunit RNA polymerase II (Pol II) is a central catalytic component that drives RNA synthesis during the transcription cycle that consists of transcription initiation, elongation, and termination. A diverse set of general transcription factors, including a multifunctional TFIIF, govern Pol II selectivity, kinetic properties, and transcription coupling with posttranscriptional processes. Here, we show that TFIIF of Arabidopsis (Arabidopsis thaliana) resembles the metazoan complex that is composed of the TFIIFα and TFIIFβ polypeptides. Arabidopsis has two TFIIFβ subunits, of which TFIIFβ1/MAN1 is essential and TFIIFβ2/MAN2 is not. In the partial loss-of-function mutant allele man1-1, the winged helix domain of Arabidopsis TFIIFβ1/MAN1 was dispensable for plant viability, whereas the cellular organization of the shoot and root apical meristems were abnormal. Forward genetic screening identified an epistatic interaction between the largest Pol II subunit nrpb1-A325V variant and the man1-1 mutation. The suppression of the man1-1 mutant developmental defects by a mutation in Pol II suggests a link between TFIIF functions in Arabidopsis transcription cycle and the maintenance of cellular organization in the shoot and root apical meristems.

A galactosyltransferase acting on arabinogalactan protein glycans is essential for embryo development in Arabidopsis

Arabinogalactan proteins (AGPs) are a complex family of cell-wall proteoglycans that are thought to play major roles in plant growth and development. Genetic approaches to studying AGP function have met limited success so far, presumably due to redundancy within the large gene families encoding AGP backbones. Here we used an alternative approach for genetic dissection of the role of AGPs in development by modifying their glycan side chains. We have identified an Arabidopsis glycosyltransferase of CAZY family GT31 (AtGALT31A) that galactosylates AGP side chains. A mutation in the AtGALT31A gene caused the arrest of embryo development at the globular stage. The presence of the transcript in the suspensor of globular-stage embryos is consistent with a role for AtGALT31A in progression of embryo development beyond the globular stage. The first observable defect in the mutant is perturbation of the formative asymmetric division of the hypophysis, indicating an essential role for AGP proteoglycans in either specification of the hypophysis or orientation of the asymmetric division plane.

Extraplastidial cytidinediphosphate diacylglycerol synthase activity is required for vegetative development in Arabidopsis thaliana

Cytidinediphosphate diacylglycerol synthase (CDS) catalyzes the activation of phosphatidic acid to cytidinediphosphate (CDP)-diacylglycerol, a central intermediate in glycerolipid biosynthesis in prokaryotic and eukaryotic organisms. Cytidinediphosphate-diacylglycerol is the precursor to phosphatidylinositol, phosphatidylglycerol (PG) and cardiolipin of eukaryotic phospholipids that are essential for various cellular functions. Isoforms of CDS are located in plastids, mitochondria and the endomembrane system of plants and are encoded by five genes in Arabidopsis. Two genes have previously been shown to code for the plastidial isoforms which are indispensable for the biosynthesis of plastidial PG, and thus biogenesis and function of thylakoid membranes. Here we have focused on the extraplastidial CDS isoforms, encoded by CDS1 and CDS2 which are constitutively expressed contrary to CDS3. We provide evidence that these closely related CDS genes code for membrane proteins located in the endoplasmic reticulum and possess very similar enzymatic properties. Development and analysis of Arabidopsis mutants lacking either one or both CDS1 and CDS2 genes clearly shows that these two genes have redundant functions. As reflected in the seedling lethal phenotype of the cds1cds2 double mutant, plant cells require at least one catalytically active microsomal CDS isoform for cell division and expansion. According to the altered glycerolipid composition of the double mutant in comparison with wild-type seedlings, it is likely that the drastic decrease in the level of phosphatidylinositol and the increase in phosphatidic acid cause defects in cell division and expansion.

Two interacting proteins are necessary for the editing of the NdhD-1 site in Arabidopsis plastids

After transcription, mRNA editing in angiosperm chloroplasts and mitochondria results in the conversion of cytidine to uridine by deamination. Analysis of Arabidopsis thaliana mutants affected in RNA editing have shown that many pentatricopeptide repeat proteins (PPRs) are required for specific cytidine deamination events. PPR proteins have been shown to be sequence-specific RNA binding proteins allowing the recognition of the C to be edited. The C-terminal DYW domain present in many editing factors has been proposed to catalyze C deamination, as it shows sequence similarities with cytidine deaminases in other organisms. However, many editing factors, such as the first to be discovered, CHLORORESPIRATORY REDUCTION4 (CRR4), lack this domain, so its importance has been unclear. Using a reverse genetic approach, we identified DYW1, an RNA editing factor acting specifically on the plastid ndhD-1 editing site recognized by CRR4. Unlike other known editing factors, DYW1 contains no identifiable PPR motifs but does contain a clear DYW domain. We were able to show interaction between CRR4 and DYW1 by bimolecular fluorescence complementation and to reconstitute a functional chimeric CRR4-DYW1 protein complementing the crr4 dyw1double mutant. We propose that CRR4 and DYW1 act together to edit the ndhD-1 site.

Identification of a DNA methylation-independent imprinting control region at the Arabidopsis MEDEA locus

Genomic imprinting is exclusive to mammals and seed plants and refers to parent-of-origin-dependent, differential transcription. As previously shown in mammals, studies in Arabidopsis have implicated DNA methylation as an important hallmark of imprinting. The current model suggests that maternally expressed imprinted genes, such as MEDEA (MEA), are activated by the DNA glycosylase DEMETER (DME), which removes DNA methylation established by the DNA methyltransferase MET1. We report the systematic functional dissection of the MEA cis-regulatory region, resulting in the identification of a 200-bp fragment that is necessary and sufficient to mediate MEA activation and imprinted expression, thus containing the imprinting control region (ICR). Notably, imprinted MEA expression mediated by this ICR is independent of DME and MET1, consistent with the lack of any significant DNA methylation in this region. This is the first example of an ICR without differential DNA methylation, suggesting that factors other than DME and MET1 are required for imprinting at the MEA locus.

Characterization of the Arabidopsis augmin complex uncovers its critical function in the assembly of the acentrosomal spindle and phragmoplast microtubule arrays

Plant cells assemble the bipolar spindle and phragmoplast microtubule (MT) arrays in the absence of the centrosome structure. Our recent findings in Arabidopsis thaliana indicated that AUGMIN subunit3 (AUG3), a homolog of animal dim γ-tubulin 3, plays a critical role in γ-tubulin-dependent MT nucleation and amplification during mitosis. Here, we report the isolation of the entire plant augmin complex that contains eight subunits. Among them, AUG1 to AUG6 share low sequence similarity with their animal counterparts, but AUG7 and AUG8 share homology only with proteins of plant origin. Genetic analyses indicate that the AUG1, AUG2, AUG4, and AUG5 genes are essential, as stable mutations in these genes could only be transmitted to heterozygous plants. The sterile aug7-1 homozygous mutant in which AUG7 expression is significantly reduced exhibited pleiotropic phenotypes of seriously retarded vegetative and reproductive growth. The aug7-1 mutation caused delocalization of γ-tubulin in the mitotic spindle and phragmoplast. Consequently, spindles were abnormally elongated, and their poles failed to converge, as MTs were splayed to discrete positions rendering deformed arrays. In addition, the mutant phragmoplasts often had disorganized MT bundles with uneven edges. We conclude that assembly of MT arrays during plant mitosis depends on the augmin complex, which includes two plant-specific subunits.

Construction and analysis of full-length and normalized cDNA libraries from citrus

We have developed an integrated method to generate a normalized cDNA collection enriched in full-length and rare transcripts from citrus, using different species and multiple tissues and developmental stages. Interpretation of ever-increasing raw sequence information generated by modern genome sequencing technologies faces multiple challenges, such as gene function analysis and genome annotation. In this regard, the availability of full-length cDNA clones facilitates functional analysis of the corresponding genes enabling manipulation of their expression and the generation of a variety of tagged versions of the native protein. The development of full-length cDNA sequences has the power to improve the quality of genome annotation, as well as provide tools for functional characterization of genes.

Augmin plays a critical role in organizing the spindle and phragmoplast microtubule arrays in Arabidopsis

In higher plant cells, microtubules (MTs) are nucleated and organized in a centrosome-independent manner. It is unclear whether augmin-dependent mechanisms underlie spindle MT organization in plant cells as they do in animal cells. When AUGMIN subunit3 (AUG3), which encodes a homolog of animal dim γ-tubulin 3/human augmin-like complex, subunit 3, was disrupted in Arabidopsis thaliana, gametogenesis frequently failed due to defects in cell division. Compared with the control microspores, which formed bipolar spindles at the cell periphery, the mutant cells often formed peripheral half spindles that only attached to condensed chromosomes or formed elongated spindles with unfocused interior poles. In addition, defective cells exhibited disorganized phragmoplast MT arrays, which caused aborted cytokinesis. The resulting pollen grains were either shrunken or contained two nuclei in an undivided cytoplasm. AUG3 was localized along MTs in the spindle and phragmoplast, and its signal was pronounced in anaphase spindle poles. An AUG3-green fluorescent protein fusion exhibited a dynamic distribution pattern, similar to that of the γ-tubulin complex protein2. When AUG3 was enriched from seedlings by affinity chromatography, AUG1 was detected by immunoblotting, suggesting an augmin-like complex was present in vivo. We conclude that augmin plays a critical role in MT organization during plant cell division.

Infection of Arabidopsis thaliana by Phytophthora parasitica and identification of variation in host specificity

Oomycete pathogens cause severe damage to a wide range of agriculturally important crops and natural ecosystems. They represent a unique group of plant pathogens that are evolutionarily distant from true fungi. In this study, we established a new plant-oomycete pathosystem in which the broad host range pathogen Phytophthora parasitica was demonstrated to be capable of interacting compatibly with the model plant Arabidopsis thaliana. Water-soaked lesions developed on leaves within 3 days and numerous sporangia formed within 5 days post-inoculation of P. parasitica zoospores. Cytological characterization showed that P. parasitica developed appressoria-like swellings and penetrated epidermal cells directly and preferably at the junction between anticlinal host cell walls. Multiple haustoria-like structures formed in both epidermal cells and mesophyll cells 1 day post-inoculation of zoospores. Pathogenicity assays of 25 A. thaliana ecotypes with six P. parasitica strains indicated the presence of a natural variation in host specificity between A. thaliana and P. parasitica. Most ecotypes were highly susceptible to P. parasitica strains Pp014, Pp016 and Pp025, but resistant to strains Pp008 and Pp009, with the frequent appearance of cell wall deposition and active defence response-based cell necrosis. Gene expression and comparative transcriptomic analysis further confirmed the compatible interaction by the identification of up-regulated genes in A. thaliana which were characteristic of biotic stress. The established A. thaliana-P. parasitica pathosystem expands the model systems investigating oomycete-plant interactions, and will facilitate a full understanding of Phytophthora biology and pathology.

Construction and characterization of a cDNA library from wheat infected with Fusarium graminearum Fg 2

Total RNA from wheat spikes infected with F. graminearum Fg2 was extracted and the mRNA was purified. Switching Mechanism at 5′ end of the RNA Transcript (SMART) technique and CDS Ill/3′ primer were used for first-strand cDNA synthesis using reverse transcriptase by RT-PCR. Primer extension polymerase chain reaction was used to construct the double-strand cDNA that was digested by proteinase K, then by Sfi I and fractionated. cDNAs longer than 0.5 kb were collected and ligated to λTriplEx2 vector followed λ phage packaging reaction and library amplification. The qualities of both unamplified and amplified cDNA libraries were strictly checked by conventional titer determination. One hundred and sixty five plaques were randomly picked and tested using PCR with universal primers derived from the sequence flanking the vector. A high quality cDNA library from wheat spikes that have been infected by F. graminearum was successfully constructed.

Targeted SNP discovery in Atlantic salmon (Salmo salar) genes using a 3'UTR-primed SNP detection approach

Background

Single nucleotide polymorphisms (SNPs) represent the most widespread type of DNA variation in vertebrates and may be used as genetic markers for a range of applications. This has led to an increased interest in identification of SNP markers in non-model species and farmed animals. The in silico SNP mining method used for discovery of most known SNPs in Atlantic salmon (Salmo salar) has applied a global (genome-wide) approach. In this study we present a targeted 3'UTR-primed SNP discovery strategy that utilizes sequence data from Salmo salar full length sequenced cDNAs (FLIcs). We compare the efficiency of this new strategy to the in silico SNP mining method when using both methods for targeted SNP discovery.

Results

The SNP discovery efficiency of the two methods was tested in a set of FLIc target genes. The 3'UTR-primed SNP discovery method detected novel SNPs in 35% of the target genes while the in silico SNP mining method detected novel SNPs in 15% of the target genes. Furthermore, the 3'UTR-primed SNP discovery strategy was the less labor intensive one and revealed a higher success rate than the in silico SNP mining method in the initial amplification step. When testing the methods we discovered 112 novel bi-allelic polymorphisms (type I markers) in 88 salmon genes [dbSNP: ss179319972-179320081, ss250608647-250608648], and three of the SNPs discovered were missense substitutions.

Conclusions

Full length insert cDNAs (FLIcs) are important genomic resources that have been developed in many farmed animals. The 3'UTR-primed SNP discovery strategy successfully utilized FLIc data to detect novel SNPs in the partially tetraploid Atlantic salmon. This strategy may therefore be useful for targeted SNP discovery in several species, and particularly useful in species that, like salmonids, have duplicated genomes.

Comparative genomic analysis of 1047 completely sequenced cDNAs from an Arabidopsis-related model halophyte, Thellungiella halophila

BACKGROUND

Thellungiella halophila (also known as T. salsuginea) is a model halophyte with a small size, short life cycle, and small genome. Thellungiella genes exhibit a high degree of sequence identity with Arabidopsis genes (90% at the cDNA level). We previously generated a full-length enriched cDNA library of T. halophila from various tissues and from whole plants treated with salinity, chilling, freezing stress, or ABA. We determined the DNA sequences of 20 000 cDNAs at both the 5'- and 3' ends, and identified 9569 distinct genes.

RESULTS

Here, we completely sequenced 1047 Thellungiella full-length cDNAs representing abiotic-stress-related genes, transcription factor genes, and protein phosphatase 2C genes. The predicted coding sequences, 5'-UTRs, and 3'-UTRs were compared with those of orthologous genes from Arabidopsis for length, sequence similarity, and structure. The 5'-UTR sequences of Thellungiella and Arabidopsis orthologs shared a significant level of similarity, although the motifs were rearranged. While examining the stress-related Thellungiella coding sequences, we found a short splicing variant of T. halophila salt overly sensitive 1 (ThSOS1), designated ThSOS1S. ThSOS1S contains the transmembrane domain of ThSOS1 but lacks the C-terminal hydrophilic region. The expression level of ThSOS1S under normal growth conditions was higher than that of ThSOS1. We also compared the expression levels of Na+-transport-system genes between Thellungiella and Arabidopsis by using full-length cDNAs from each species as probes. Several genes that play essential roles in Na+ excretion, compartmentation, and diffusion (SOS1, SOS2, NHX1, and HKT1) were expressed at higher levels in Thellungiella than in Arabidopsis.

CONCLUSIONS

The full-length cDNA sequences obtained in this study will be essential for the ongoing annotation of the Thellungiella genome, especially for further improvement of gene prediction. Moreover, they will enable us to find splicing variants such as ThSOS1S (AB562331).

Identification and characterization of full-length cDNAs in channel catfish (Ictalurus punctatus) and blue catfish (Ictalurus furcatus)

BACKGROUND

Genome annotation projects, gene functional studies, and phylogenetic analyses for a given organism all greatly benefit from access to a validated full-length cDNA resource. While increasingly common in model species, full-length cDNA resources in aquaculture species are scarce.

METHODOLOGY AND PRINCIPAL FINDINGS

Through in silico analysis of catfish (Ictalurus spp.) ESTs, a total of 10,037 channel catfish and 7,382 blue catfish cDNA clones were identified as potentially encoding full-length cDNAs. Of this set, a total of 1,169 channel catfish and 933 blue catfish full-length cDNA clones were selected for re-sequencing to provide additional coverage and ensure sequence accuracy. A total of 1,745 unique gene transcripts were identified from the full-length cDNA set, including 1,064 gene transcripts from channel catfish and 681 gene transcripts from blue catfish, with 416 transcripts shared between the two closely related species. Full-length sequence characteristics (ortholog conservation, UTR length, Kozak sequence, and conserved motifs) of the channel and blue catfish were examined in detail. Comparison of gene ontology composition between full-length cDNAs and all catfish ESTs revealed that the full-length cDNA set is representative of the gene diversity encoded in the catfish transcriptome.

CONCLUSIONS

This study describes the first catfish full-length cDNA set constructed from several cDNA libraries. The catfish full-length cDNA sequences, and data gleaned from sequence characteristics analysis, will be a valuable resource for ongoing catfish whole-genome sequencing and future gene-based studies of function and evolution in teleost fishes.

Two closely related genes of Arabidopsis encode plastidial cytidinediphosphate diacylglycerol synthases essential for photoautotrophic growth

Cytidinediphosphate diacylglycerol synthase (CDS) catalyzes the formation of cytidinediphosphate diacylglycerol, an essential precursor of anionic phosphoglycerolipids like phosphatidylglycerol or -inositol. In plant cells, CDS isozymes are located in plastids, mitochondria, and microsomes. Here, we show that these isozymes are encoded by five genes in Arabidopsis (Arabidopsis thaliana). Alternative translation initiation or alternative splicing of CDS2 and CDS4 transcripts can result in up to 10 isoforms. Most of the cDNAs encoding the various plant isoforms were functionally expressed in yeast and rescued the nonviable phenotype of the mutant strain lacking CDS activity. The closely related genes CDS4 and CDS5 were found to encode plastidial isozymes with similar catalytic properties. Inactivation of both genes was required to obtain Arabidopsis mutant lines with a visible phenotype, suggesting that the genes have redundant functions. Analysis of these Arabidopsis mutants provided further independent evidence for the importance of plastidial phosphatidylglycerol for structure and function of thylakoid membranes and, hence, for photoautotrophic growth.

Ins and outs of systems biology vis-à-vis molecular biology: continuation or clear cut?

The comprehension of living organisms in all their complexity poses a major challenge to the biological sciences. Recently, systems biology has been proposed as a new candidate in the development of such a comprehension. The main objective of this paper is to address what systems biology is and how it is practised. To this end, the basic tools of a systems biological approach are explored and illustrated. In addition, it is questioned whether systems biology 'revolutionizes' molecular biology and 'transcends' its assumed reductionism. The strength of this claim appears to depend on how molecular and systems biology are characterised and on how reductionism is interpreted. Doing credit to molecular biology and to methodological reductionism, it is argued that the distinction between molecular and systems biology is gradual rather than sharp. As such, the classical challenge in biology to manage, interpret and integrate biological data into functional wholes is further intensified by systems biology's use of modelling and bioinformatics, and by its scale enlargement.

Identification and analysis of in planta expressed genes of Magnaporthe oryzae

Background

Infection of plants by pathogens and the subsequent disease development involves substantial changes in the biochemistry and physiology of both partners. Analysis of genes that are expressed during these interactions represents a powerful strategy to obtain insights into the molecular events underlying these changes. We have employed expressed sequence tag (EST) analysis to identify rice genes involved in defense responses against infection by the blast fungus Magnaporthe oryzae and fungal genes involved in infectious growth within the host during a compatible interaction.

Results

A cDNA library was constructed with RNA from rice leaves (Oryza sativa cv. Hwacheong) infected with M. oryzae strain KJ201. To enrich for fungal genes, subtraction library using PCR-based suppression subtractive hybridization was constructed with RNA from infected rice leaves as a tester and that from uninfected rice leaves as the driver. A total of 4,148 clones from two libraries were sequenced to generate 2,302 non-redundant ESTs. Of these, 712 and 1,562 ESTs could be identified to encode fungal and rice genes, respectively. To predict gene function, Gene Ontology (GO) analysis was applied, with 31% and 32% of rice and fungal ESTs being assigned to GO terms, respectively. One hundred uniESTs were found to be specific to fungal infection EST. More than 80 full-length fungal cDNA sequences were used to validate ab initio annotated gene model of M. oryzae genome sequence.

Conclusion

This study shows the power of ESTs to refine genome annotation and functional characterization. Results of this work have advanced our understanding of the molecular mechanisms underpinning fungal-plant interactions and formed the basis for new hypothesis.

Characterization of full-length sequenced cDNA inserts (FLIcs) from Atlantic salmon (Salmo salar)

Background

Sequencing of the Atlantic salmon genome is now being planned by an international research consortium. Full-length sequenced inserts from cDNAs (FLIcs) are an important tool for correct annotation and clustering of the genomic sequence in any species. The large amount of highly similar duplicate sequences caused by the relatively recent genome duplication in the salmonid ancestor represents a particular challenge for the genome project. FLIcs will therefore be an extremely useful resource for the Atlantic salmon sequencing project. In addition to be helpful in order to distinguish between duplicate genome regions and in determining correct gene structures, FLIcs are an important resource for functional genomic studies and for investigation of regulatory elements controlling gene expression. In contrast to the large number of ESTs available, including the ESTs from 23 developmental and tissue specific cDNA libraries contributed by the Salmon Genome Project (SGP), the number of sequences where the full-length of the cDNA insert has been determined has been small.

Results

High quality full-length insert sequences from 560 pre-smolt white muscle tissue specific cDNAs were generated, accession numbers [GenBank: BT043497 - BT044056]. Five hundred and ten (91%) of the transcripts were annotated using Gene Ontology (GO) terms and 440 of the FLIcs are likely to contain a complete coding sequence (cCDS). The sequence information was used to identify putative paralogs, characterize salmon Kozak motifs, polyadenylation signal variation and to identify motifs likely to be involved in the regulation of particular genes. Finally, conserved 7-mers in the 3'UTRs were identified, of which some were identical to miRNA target sequences.

Conclusion

This paper describes the first Atlantic salmon FLIcs from a tissue and developmental stage specific cDNA library. We have demonstrated that many FLIcs contained a complete coding sequence (cCDS). This suggests that the remaining cDNA libraries generated by SGP represent a valuable cCDS FLIc source. The conservation of 7-mers in 3'UTRs indicates that these motifs are functionally important. Identity between some of these 7-mers and miRNA target sequences suggests that they are miRNA targets in Salmo salar transcripts as well.

A new set of ESTs and cDNA clones from full-length and normalized libraries for gene discovery and functional characterization in citrus

Background

Interpretation of ever-increasing raw sequence information generated by modern genome sequencing technologies faces multiple challenges, such as gene function analysis and genome annotation. Indeed, nearly 40% of genes in plants encode proteins of unknown function. Functional characterization of these genes is one of the main challenges in modern biology. In this regard, the availability of full-length cDNA clones may fill in the gap created between sequence information and biological knowledge. Full-length cDNA clones facilitate functional analysis of the corresponding genes enabling manipulation of their expression in heterologous systems and the generation of a variety of tagged versions of the native protein. In addition, the development of full-length cDNA sequences has the power to improve the quality of genome annotation.

Results

We developed an integrated method to generate a new normalized EST collection enriched in full-length and rare transcripts of different citrus species from multiple tissues and developmental stages. We constructed a total of 15 cDNA libraries, from which we isolated 10,898 high-quality ESTs representing 6142 different genes. Percentages of redundancy and proportion of full-length clones range from 8 to 33, and 67 to 85, respectively, indicating good efficiency of the approach employed. The new EST collection adds 2113 new citrus ESTs, representing 1831 unigenes, to the collection of citrus genes available in the public databases. To facilitate functional analysis, cDNAs were introduced in a Gateway-based cloning vector for high-throughput functional analysis of genes in planta. Herein, we describe the technical methods used in the library construction, sequence analysis of clones and the overexpression of CitrSEP, a citrus homolog to the Arabidopsis SEP3 gene, in Arabidopsis as an example of a practical application of the engineered Gateway vector for functional analysis.

Conclusion

The new EST collection denotes an important step towards the identification of all genes in the citrus genome. Furthermore, public availability of the cDNA clones generated in this study, and not only their sequence, enables testing of the biological function of the genes represented in the collection. Expression of the citrus SEP3 homologue, CitrSEP, in Arabidopsis results in early flowering, along with other phenotypes resembling the over-expression of the Arabidopsis SEPALLATA genes. Our findings suggest that the members of the SEP gene family play similar roles in these quite distant plant species.

Clone-based functional genomics

Annotated genomes have provided a wealth of information about gene structure and gene catalogs in a wide range of species. Taking advantage of these developments, novel techniques have been implemented to investigate systematically diverse aspects of gene and protein functions underpinning biology processes. Here, we review functional genomics applications that require the mass production of cloned sequence repertoires, including ORFeomes and silencing tag collections. We discuss the techniques employed in large-scale cloning projects and we provide an up-to-date overview of the clone resources available for model plant species and of the current applications that may be scaled up for systematic plant gene studies.

Functional genomics using RIKEN Arabidopsis thaliana full-length cDNAs

Full-length cDNAs are essential for the correct annotation of genomic sequences as well as for the functional analysis of genes and their products. We have isolated about 240,000 RIKEN Arabidopsis full-length (RAFL) cDNA clones. These clones were clustered into about 17,000 non-redundant cDNA groups, i.e., about 60% of all Arabidopsis predicted genes. The sequence information of the RAFL cDNAs is useful for promoter analysis, and for the correct annotation of predicted transcriptional units and gene products. We prepared cDNA microarrays containing independent full-length cDNA groups and studied the expression profiles of genes under various stress- and hormone-treatment conditions, and in various mutants and transgenic plants. These expression profiling studies have shown the expression levels of many genes as a detailed snapshot describing the state of a biological system in planta under various conditions. We have applied RAFL cDNAs to the functional analysis of proteins using the full-length cDNA over-expressing (FOX) gene hunting system and the wheat germ cell-free protein synthesis system. The RAFL cDNA collection was also used for determination of the domain structure of proteins by NMR. In this review, we summarize the present state and perspectives of functional genomics using RAFL cDNAs.

A Plasmodium falciparum FcB1-schizont-EST collection providing clues to schizont specific gene structure and polymorphism

BACKGROUND

The Plasmodium falciparum genome (3D7 strain) published in 2002, revealed ~5,400 genes, mostly based on in silico predictions. Experimental data is therefore required for structural and functional assessments of P. falciparum genes and expression, and polymorphic data are further necessary to exploit genomic information to further qualify therapeutic target candidates. Here, we undertook a large scale analysis of a P. falciparum FcB1-schizont-EST library previously constructed by suppression subtractive hybridization (SSH) to study genes expressed during merozoite morphogenesis, with the aim of: 1) obtaining an exhaustive collection of schizont specific ESTs, 2) experimentally validating or correcting P. falciparum gene models and 3) pinpointing genes displaying protein polymorphism between the FcB1 and 3D7 strains.

RESULTS

A total of 22,125 clones randomly picked from the SSH library were sequenced, yielding 21,805 usable ESTs that were then clustered on the P. falciparum genome. This allowed identification of 243 protein coding genes, including 121 previously annotated as hypothetical. Statistical analysis of GO terms, when available, indicated significant enrichment in genes involved in "entry into host-cells" and "actin cytoskeleton". Although most ESTs do not span full-length gene reading frames, detailed sequence comparison of FcB1-ESTs versus 3D7 genomic sequences allowed the confirmation of exon/intron boundaries in 29 genes, the detection of new boundaries in 14 genes and identification of protein polymorphism for 21 genes. In addition, a large number of non-protein coding ESTs were identified, mainly matching with the two A-type rRNA units (on chromosomes 5 and 7) and to a lower extent, two atypical rRNA loci (on chromosomes 1 and 8), TARE subtelomeric regions (several chromosomes) and the recently described telomerase RNA gene (chromosome 9).

CONCLUSION

This FcB1-schizont-EST analysis confirmed the actual expression of 243 protein coding genes, allowing the correction of structural annotations for a quarter of these sequences. In addition, this analysis demonstrated the actual transcription of several remarkable non-protein coding loci: 2 atypical rRNA, TARE region and telomerase RNA gene. Together with other collections of P. falciparum ESTs, usually generated from mixed parasite stages, this collection of FcB1-schizont-ESTs provides valuable data to gain further insight into the P. falciparum gene structure, polymorphism and expression.

Possible diversifying selection in the imprinted gene, MEDEA, in Arabidopsis

Coevolutionary conflict among imprinted genes that influence traits such as offspring growth may arise when maternal and paternal genomes have different evolutionary optima. This conflict is expected in outcrossing taxa with multiple paternity, but not self-fertilizing taxa. MEDEA (MEA) is an imprinted plant gene that influences seed growth. Disagreement exists regarding the type of selection acting on this gene. We present new data and analyses of sequence diversity of MEA in self-fertilizing and outcrossing Arabidopsis and its relatives, to help clarify the form of selection acting on this gene. Codon-based branch analysis among taxa (PAML) suggests that selection on the coding region is changing over time, and nonsynonymous substitution is elevated in at least one outcrossing branch. Codon-based analysis of diversity within outcrossing Arabidopsis lyrata ssp. petraea (OmegaMap) suggests that diversifying selection is acting on a portion of the gene, to cause elevated nonsynonymous polymorphism. Providing further support for balancing selection in A. lyrata, Hudson, Kreitman and Aguadé analysis indicates that diversity/divergence at silent sites in the MEA promoter and genic region is elevated relative to reference genes, and there are deviations from the neutral frequency spectrum. This combination of positive selection as well as balancing and diversifying selection in outcrossing lineages is consistent with other genes influence by evolutionary conflict, such as disease resistance genes. Consistent with predictions that conflict would be eliminated in self-fertilizing taxa, we found no evidence of positive, balancing, or diversifying selection in A. thaliana promoter or genic region.

BTB and TAZ domain scaffold proteins perform a crucial function in Arabidopsis development

In Arabidopsis, bric-a-brac, tramtrack and broad (BTB) domain scaffold proteins form a family of 80 proteins that have involvement in various signaling pathways. The five members of the subfamily of BTB AND TAZ DOMAIN proteins (BT1-BT5) have a typical domain structure that is only observed in land plants. Here, we present a functional analysis of the BT family, of which at least four members are encoded by auxin-responsive genes. BT1 is a short-lived protein that is characteristically targeted for degradation by the 26S proteasome. Expression pattern, gene structure and sequence analyses indicate that BT1 and BT2 are closely related. They both localize to the nucleus and the cytosol, whereas the remaining BT proteins were determined as cytosolic proteins. Detailed molecular and phenotypic analysis of plants segregating for null mutations in the BT family revealed substantial redundancy among the BT members, and highlighted that BT proteins perform crucial roles in both male and female gametophyte development. BT2 seems to be the predominant gene in this process, in which it is functionally replaced by BT3 and BT1 through reciprocal transcription regulation. Compensational expression alters the steady-state mRNA levels among the remaining BT family members when other BT members are lost, and this contributes towards functional redundancy. Our data provide a surprising example of functional redundancy among genes required during gametophyte development, something that could not be detected in the current screens for gametophyte mutants.

Annotating genomes with massive-scale RNA sequencing

A method for de novo genome annotation using high-throughput cDNA sequencing data.

Next generation technologies enable massive-scale cDNA sequencing (so-called RNA-Seq). Mainly because of the difficulty of aligning short reads on exon-exon junctions, no attempts have been made so far to use RNA-Seq for building gene models de novo, that is, in the absence of a set of known genes and/or splicing events. We present G-Mo.R-Se (Gene Modelling using RNA-Seq), an approach aimed at building gene models directly from RNA-Seq and demonstrate its utility on the grapevine genome.

Novel long non-protein coding RNAs involved in Arabidopsis differentiation and stress responses

Long non-protein coding RNAs (npcRNA) represent an emerging class of riboregulators, which either act directly in this long form or are processed to shorter miRNA and siRNA. Genome-wide bioinformatic analysis of full-length cDNA databases identified 76 Arabidopsis npcRNAs. Fourteen npcRNAs were antisense to protein-coding mRNAs, suggesting cis-regulatory roles. Numerous 24-nt siRNA matched to five different npcRNAs, suggesting that these npcRNAs are precursors of this type of siRNA. Expression analyses of the 76 npcRNAs identified a novel npcRNA that accumulates in a dcl1 mutant but does not appear to produce trans-acting siRNA or miRNA. Additionally, another npcRNA was the precursor of miR869 and shown to be up-regulated in dcl4 but not in dcl1 mutants, indicative of a young miRNA gene. Abiotic stress altered the accumulation of 22 npcRNAs among the 76, a fraction significantly higher than that observed for the RNA binding protein-coding fraction of the transcriptome. Overexpression analyses in Arabidopsis identified two npcRNAs as regulators of root growth during salt stress and leaf morphology, respectively. Hence, together with small RNAs, long npcRNAs encompass a sensitive component of the transcriptome that have diverse roles during growth and differentiation.

The Brazilian contribution to the study of the Schistosoma mansoni transcriptome

Brazilian scientists have made definitive contributions towards the understanding of the transcriptome of Schistosoma mansoni. The transcriptome, as a collection of expressed genes, was studied with two basic approaches, the sequencing of cloned cDNA ends and the production and sequencing of random fragments of cDNA. The data that was generated covers nearly the entire set of transcribed genes. These data have been important for gene discovery and annotation of gene function, annotation of the genome, construction of microarrays, identification of peptides in the study of the proteome and can be mined for splice variants, polymorphisms such as single nucleotide polymorphisms and microsatellites, for example. This review will describe the advances achieved with the contribution of Brazilian researchers to the understanding of the transcribed regions of the S. mansoni genome.

Web-based Arabidopsis functional and structural genomics resources

As plant research moves to a "post-genomic" era, many diverse internet resources become available to the international research community. Arabidopsis thaliana, because of its small size, rapid life cycle and simple genome, has been a model system for decades, with much research funding and many projects devoted to creation of functional and structural genomics resources. Different types of data, including genome, transcriptome, proteome, phenome, metabolome and ionome are stored in these resources. In this chapter, a variety of genomics resources are introduced, with simple descriptions of how some can be accessed by laboratory researchers via the internet.

A conifer genomics resource of 200,000 spruce (Picea spp.) ESTs and 6,464 high-quality, sequence-finished full-length cDNAs for Sitka spruce (Picea sitchensis)

Background

Members of the pine family (Pinaceae), especially species of spruce (Picea spp.) and pine (Pinus spp.), dominate many of the world's temperate and boreal forests. These conifer forests are of critical importance for global ecosystem stability and biodiversity. They also provide the majority of the world's wood and fiber supply and serve as a renewable resource for other industrial biomaterials. In contrast to angiosperms, functional and comparative genomics research on conifers, or other gymnosperms, is limited by the lack of a relevant reference genome sequence. Sequence-finished full-length (FL)cDNAs and large collections of expressed sequence tags (ESTs) are essential for gene discovery, functional genomics, and for future efforts of conifer genome annotation.

Results

As part of a conifer genomics program to characterize defense against insects and adaptation to local environments, and to discover genes for the production of biomaterials, we developed 20 standard, normalized or full-length enriched cDNA libraries from Sitka spruce (P. sitchensis), white spruce (P. glauca), and interior spruce (P. glauca-engelmannii complex). We sequenced and analyzed 206,875 3'- or 5'-end ESTs from these libraries, and developed a resource of 6,464 high-quality sequence-finished FLcDNAs from Sitka spruce. Clustering and assembly of 147,146 3'-end ESTs resulted in 19,941 contigs and 26,804 singletons, representing 46,745 putative unique transcripts (PUTs). The 6,464 FLcDNAs were all obtained from a single Sitka spruce genotype and represent 5,718 PUTs.

Conclusion

This paper provides detailed annotation and quality assessment of a large EST and FLcDNA resource for spruce. The 6,464 Sitka spruce FLcDNAs represent the third largest sequence-verified FLcDNA resource for any plant species, behind only rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana), and the only substantial FLcDNA resource for a gymnosperm. Our emphasis on capturing FLcDNAs and ESTs from cDNA libraries representing herbivore-, wound- or elicitor-treated induced spruce tissues, along with incorporating normalization to capture rare transcripts, resulted in a rich resource for functional genomics and proteomics studies. Sequence comparisons against five plant genomes and the non-redundant GenBank protein database revealed that a substantial number of spruce transcripts have no obvious similarity to known angiosperm gene sequences. Opportunities for future applications of the sequence and clone resources for comparative and functional genomics are discussed.

The genome sequence of the model ascomycete fungus Podospora anserina

A 10X draft sequence of Podospora anserina genome shows highly dynamic evolution since its divergence from Neurospora crassa.

Background

The dung-inhabiting ascomycete fungus Podospora anserina is a model used to study various aspects of eukaryotic and fungal biology, such as ageing, prions and sexual development.

Results

We present a 10X draft sequence of P. anserina genome, linked to the sequences of a large expressed sequence tag collection. Similar to higher eukaryotes, the P. anserina transcription/splicing machinery generates numerous non-conventional transcripts. Comparison of the P. anserina genome and orthologous gene set with the one of its close relatives, Neurospora crassa, shows that synteny is poorly conserved, the main result of evolution being gene shuffling in the same chromosome. The P. anserina genome contains fewer repeated sequences and has evolved new genes by duplication since its separation from N. crassa, despite the presence of the repeat induced point mutation mechanism that mutates duplicated sequences. We also provide evidence that frequent gene loss took place in the lineages leading to P. anserina and N. crassa. P. anserina contains a large and highly specialized set of genes involved in utilization of natural carbon sources commonly found in its natural biotope. It includes genes potentially involved in lignin degradation and efficient cellulose breakdown.

Conclusion

The features of the P. anserina genome indicate a highly dynamic evolution since the divergence of P. anserina and N. crassa, leading to the ability of the former to use specific complex carbon sources that match its needs in its natural biotope.

Riboregulators in plant development

npcRNA (non-protein-coding RNAs) are an emerging class of regulators, so-called riboregulators, and include a large diversity of small RNAs [miRNAs (microRNAs)/siRNAs (small interfering RNAs)] that are involved in various developmental processes in plants and animals. In addition, several other npcRNAs encompassing various transcript sizes (up to several kilobases) have been identified using different genomic approaches. Much less is known about the mechanism of action of these other classes of riboregulators also present in the cell. The organogenesis of nitrogen-fixing nodules in legume plants is initiated in specific root cortical cells that express the npcRNA MtENOD40 (Medicago truncatula early nodulin 40). We have identified a novel RBP (RNA-binding protein), MtRBP1 (M. truncatula RBP 1), which interacts with the MtENOD40 RNA, and is exported into the cytoplasm during legume nodule development in the region expressing MtENOD40. A direct involvement of the MtENOD40 RNA in the relocalization of this RBP into cytoplasmic granules could be demonstrated, revealing a new RNA function in the cell. To extend these results, we searched for npcRNAs in the model plant Arabidopsis thaliana whose genome is completely known. We have identified 86 novel npcRNAs from which 27 corresponded to antisense RNAs of known coding regions. Using a dedicated 'macroarray' containing these npcRNAs and a collection of RBPs, we characterized their regulation in different tissues and plants subjected to environmental stresses. Most of the npcRNAs showed high variations in gene expression in contrast with the RBP genes. Recent large-scale analysis of the sRNA component of the transcriptome revealed an enormous diversity of siRNAs/miRNAs in the Arabidopsis genome. Bioinformatic analysis revealed that 34 large npcRNAs are precursors of siRNAs/miRNAs. npcRNAs, which are a sensitive component of the transcriptome, may reveal novel riboregulatory mechanisms involved in post-transcriptional control of differentiation or environmental responses.

Analysis of 4,664 high-quality sequence-finished poplar full-length cDNA clones and their utility for the discovery of genes responding to insect feeding

Background

The genus Populus includes poplars, aspens and cottonwoods, which will be collectively referred to as poplars hereafter unless otherwise specified. Poplars are the dominant tree species in many forest ecosystems in the Northern Hemisphere and are of substantial economic value in plantation forestry. Poplar has been established as a model system for genomics studies of growth, development, and adaptation of woody perennial plants including secondary xylem formation, dormancy, adaptation to local environments, and biotic interactions.

Results

As part of the poplar genome sequencing project and the development of genomic resources for poplar, we have generated a full-length (FL)-cDNA collection using the biotinylated CAP trapper method. We constructed four FLcDNA libraries using RNA from xylem, phloem and cambium, and green shoot tips and leaves from the P. trichocarpa Nisqually-1 genotype, as well as insect-attacked leaves of the P. trichocarpa × P. deltoides hybrid. Following careful selection of candidate cDNA clones, we used a combined strategy of paired end reads and primer walking to generate a set of 4,664 high-accuracy, sequence-verified FLcDNAs, which clustered into 3,990 putative unique genes. Mapping FLcDNAs to the poplar genome sequence combined with BLAST comparisons to previously predicted protein coding sequences in the poplar genome identified 39 FLcDNAs that likely localize to gaps in the current genome sequence assembly. Another 173 FLcDNAs mapped to the genome sequence but were not included among the previously predicted genes in the poplar genome. Comparative sequence analysis against Arabidopsis thaliana and other species in the non-redundant database of GenBank revealed that 11.5% of the poplar FLcDNAs display no significant sequence similarity to other plant proteins. By mapping the poplar FLcDNAs against transcriptome data previously obtained with a 15.5 K cDNA microarray, we identified 153 FLcDNA clones for genes that were differentially expressed in poplar leaves attacked by forest tent caterpillars.

Conclusion

This study has generated a high-quality FLcDNA resource for poplar and the third largest FLcDNA collection published to date for any plant species. We successfully used the FLcDNA sequences to reassess gene prediction in the poplar genome sequence, perform comparative sequence annotation, and identify differentially expressed transcripts associated with defense against insects. The FLcDNA sequences will be essential to the ongoing curation and annotation of the poplar genome, in particular for targeting gaps in the current genome assembly and further improvement of gene predictions. The physical FLcDNA clones will serve as useful reagents for functional genomics research in areas such as analysis of gene functions in defense against insects and perennial growth. Sequences from this study have been deposited in NCBI GenBank under the accession numbers EF144175 to EF148838.

A bifunctional locus (BIO3-BIO1) required for biotin biosynthesis in Arabidopsis

We identify here the Arabidopsis (Arabidopsis thaliana) gene encoding the third enzyme in the biotin biosynthetic pathway, dethiobiotin synthetase (BIO3; At5g57600). This gene is positioned immediately upstream of BIO1, which is known to be associated with the second reaction in the pathway. Reverse genetic analysis demonstrates that bio3 insertion mutants have a similar phenotype to the bio1 and bio2 auxotrophs identified using forward genetic screens for arrested embryos rescued on enriched nutrient medium. Unexpectedly, bio3 and bio1 mutants define a single genetic complementation group. Reverse transcription-polymerase chain reaction analysis demonstrates that separate BIO3 and BIO1 transcripts and two different types of chimeric BIO3-BIO1 transcripts are produced. Consistent with genetic data, one of the fused transcripts is monocistronic and encodes a bifunctional fusion protein. A splice variant is bicistronic, with distinct but overlapping reading frames. The dual functionality of the monocistronic transcript was confirmed by complementing the orthologous auxotrophs of Escherichia coli (bioD and bioA). BIO3-BIO1 transcripts from other plants provide further evidence for differential splicing, existence of a fusion protein, and localization of both enzymatic reactions to mitochondria. In contrast to most biosynthetic enzymes in eukaryotes, which are encoded by genes dispersed throughout the genome, biotin biosynthesis in Arabidopsis provides an intriguing example of a bifunctional locus that catalyzes two sequential reactions in the same metabolic pathway. This complex locus exhibits several unusual features that distinguish it from biotin operons in bacteria and from other genes known to encode bifunctional enzymes in plants.

A collection of 10,096 indica rice full-length cDNAs reveals highly expressed sequence divergence between Oryza sativa indica and japonica subspecies

Relatively few indica rice full-length cDNAs were available to aid in the annotation of rice genes. The data presented here described the sequencing and analysis of 10,096 full-length cDNAs from Oryza sativa subspecies indica Guangluai 4. Of them, 9,029 matched rice genomic sequences in publicly-available databases, and 1,200 were identified as new rice genes. Comparison with the knowledge-based Oryza Molecular Biological Encyclopedia japonica cDNA collection indicated that 3,316 (41.6%) of the 7,965 indica-japonica cDNA pairs showed no distinct variations at protein level (2,117 indica-japonica cDNA pairs showed fully identical and 1,199 indica-japonica cDNA pairs showed no frame shift). Moreover, 3,645 (45.8%) of the indica-japonica pairs showed substantial differences at the protein level due to single nucleotide polymorphisms (SNPs), insertions or deletions, and sequence-segment variations between indica and japonica subspecies. Further experimental verifications using PCR screening and quantitative reverse transcriptional PCR revealed unique transcripts for indica subspecies. Comparative analysis also showed that most of rice genes were evolved under purifying selection. These variations might distinguish the phenotypic changes of the two cultivated rice subspecies indica and japonica. Analysis of these cDNAs extends known rice genes and identifies new ones in rice.

A cluster of disease resistance genes in Arabidopsis is coordinately regulated by transcriptional activation and RNA silencing

The RPP5 (for recognition of Peronospora parasitica 5) locus in the Arabidopsis thaliana Columbia strain contains a cluster of paralogous disease Resistance (R) genes that play important roles in innate immunity. Among the R genes in this locus, RPP4 confers resistance to two races of the fungal pathogen Hyaloperonospora parasitica, while activation of SNC1 (for suppressor of npr1-1, constitutive 1) results in the resistance to another race of H. parasitica and to pathovars of the bacterial pathogen Pseudomonas syringae through the accumulation of salicylic acid (SA). Here, we demonstrate that other Columbia RPP5 locus R genes can be induced by transgenic overexpression of SNC1, which itself is regulated by a positive amplification loop involving SA accumulation. We also show that small RNA species that can target RPP5 locus R genes are produced in wild-type plants and that these R genes can be cosuppressed in transgenic plants overexpressing SNC1. Steady state expression levels of SNC1 increase in some mutants (dcl4-4, ago1-36, and upf1-5) defective in RNA silencing as well as in transgenic plants expressing the P1/Helper Component-Protease viral suppressor of RNA silencing. However, steady state levels of small RNA species do not change in mutants that upregulate SNC1. These data indicate many Columbia RPP5 locus R genes can be coordinately regulated both positively and negatively and suggest that the RPP5 locus is poised to respond to pathogens that disturb RNA silencing.

Identification of novel conserved peptide uORF homology groups in Arabidopsis and rice reveals ancient eukaryotic origin of select groups and preferential association with transcription factor-encoding genes

BACKGROUND

Upstream open reading frames (uORFs) can mediate translational control over the largest, or major ORF (mORF) in response to starvation, polyamine concentrations, and sucrose concentrations. One plant uORF with conserved peptide sequences has been shown to exert this control in an amino acid sequence-dependent manner but generally it is not clear what kinds of genes are regulated, or how extensively this mechanism is invoked in a given genome.

RESULTS

By comparing full-length cDNA sequences from Arabidopsis and rice we identified 26 distinct homology groups of conserved peptide uORFs, only three of which have been reported previously. Pairwise Ka/Ks analysis showed that purifying selection had acted on nearly all conserved peptide uORFs and their associated mORFs. Functions of predicted mORF proteins could be inferred for 16 homology groups and many of these proteins appear to have a regulatory function, including 6 transcription factors, 5 signal transduction factors, 3 developmental signal molecules, a homolog of translation initiation factor eIF5, and a RING finger protein. Transcription factors are clearly overrepresented in this data set when compared to the frequency calculated for the entire genome (p = 1.2 x 10(-7)). Duplicate gene pairs arising from a whole genome duplication (ohnologs) with a conserved uORF are much more likely to have been retained in Arabidopsis (Arabidopsis thaliana) than are ohnologs of other genes (39% vs 14% of ancestral genes, p = 5 x 10(-3)). Two uORF groups were found in animals, indicating an ancient origin of these putative regulatory elements.

CONCLUSION

Conservation of uORF amino acid sequence, association with homologous mORFs over long evolutionary time periods, preferential retention after whole genome duplications, and preferential association with mORFs coding for transcription factors suggest that the conserved peptide uORFs identified in this study are strong candidates for translational controllers of regulatory genes.

The Schistosoma mansoni transcriptome: an update

Large scale EST sequencing projects have been carried out for Schistosoma mansoni and Schistosoma japonicum. This update will briefly review the most recent accomplishments in the area and discuss the use of EST data for the purposes of gene discovery, gene model development, genome annotation and SNP analysis. In addition, the use of ESTs for studying other features of the transcriptome such as splice site and transcription initiation variants will be discussed as well as approaches to assigning function to unknown transcripts. Although EST sequencing has contributed much for schistosome research, other data mining possibilities exist, including the identification of putative drug and vaccine targets.

Gene models from ESTs (GeneModelEST): an application on the Solanum lycopersicum genome

BACKGROUND

The structure annotation of a genome is based either on ab initio methodologies or on similaritiy searches versus molecules that have been already annotated. Ab initio gene predictions in a genome are based on a priori knowledge of species-specific features of genes. The training of ab initio gene finders is based on the definition of a data-set of gene models. To accomplish this task the common approach is to align species-specific full length cDNA and EST sequences along the genomic sequences in order to define exon/intron structure of mRNA coding genes.

RESULTS

GeneModelEST is the software here proposed for defining a data-set of candidate gene models using exclusively evidence derived from cDNA/EST sequences.GeneModelEST requires the genome coordinates of the spliced-alignments of ESTs and of contigs (tentative consensus sequences) generated by an EST clustering/assembling procedure to be formatted in a General Feature Format (GFF) standard file. Moreover, the alignments of the contigs versus a protein database are required as an NCBI BLAST formatted report file. The GeneModelEST analysis aims to i) evaluate each exon as defined from contig spliced alignments onto the genome sequence; ii) classify the contigs according to quality levels in order to select candidate gene models; iii) assign to the candidate gene models preliminary functional annotations. We discuss the application of the proposed methodology to build a data-set of gene models of Solanum lycopersicum, whose genome sequencing is an ongoing effort by the International Tomato Genome Sequencing Consortium.

CONCLUSION

The contig classification procedure used by GeneModelEST supports the detection of candidate gene models, the identification of potential alternative transcripts and it is useful to filter out ambiguous information. An automated procedure, such as the one proposed here, is fundamental to support large scale analysis in order to provide species-specific gene models, that could be useful as a training data-set for ab initio gene finders and/or as a reference gene list for a human curated annotation.

PIC1, an ancient permease in Arabidopsis chloroplasts, mediates iron transport

In chloroplasts, the transition metals iron and copper play an essential role in photosynthetic electron transport and act as cofactors for superoxide dismutases. Iron is essential for chlorophyll biosynthesis, and ferritin clusters in plastids store iron during germination, development, and iron stress. Thus, plastidic homeostasis of transition metals, in particular of iron, is crucial for chloroplast as well as plant development. However, very little is known about iron uptake by chloroplasts. Arabidopsis thaliana PERMEASE IN CHLOROPLASTS1 (PIC1), identified in a screen for metal transporters in plastids, contains four predicted alpha-helices, is targeted to the inner envelope, and displays homology with cyanobacterial permease-like proteins. Knockout mutants of PIC1 grew only heterotrophically and were characterized by a chlorotic and dwarfish phenotype reminiscent of iron-deficient plants. Ultrastructural analysis of plastids revealed severely impaired chloroplast development and a striking increase in ferritin clusters. Besides upregulation of ferritin, pic1 mutants showed differential regulation of genes and proteins related to iron stress or transport, photosynthesis, and Fe-S cluster biogenesis. Furthermore, PIC1 and its cyanobacterial homolog mediated iron accumulation in an iron uptake-defective yeast mutant. These observations suggest that PIC1 functions in iron transport across the inner envelope of chloroplasts and hence in cellular metal homeostasis.

PIC1, an ancient permease in Arabidopsis chloroplasts, mediates iron transport

In chloroplasts, the transition metals iron and copper play an essential role in photosynthetic electron transport and act as cofactors for superoxide dismutases. Iron is essential for chlorophyll biosynthesis, and ferritin clusters in plastids store iron during germination, development, and iron stress. Thus, plastidic homeostasis of transition metals, in particular of iron, is crucial for chloroplast as well as plant development. However, very little is known about iron uptake by chloroplasts. Arabidopsis thaliana PERMEASE IN CHLOROPLASTS1 (PIC1), identified in a screen for metal transporters in plastids, contains four predicted alpha-helices, is targeted to the inner envelope, and displays homology with cyanobacterial permease-like proteins. Knockout mutants of PIC1 grew only heterotrophically and were characterized by a chlorotic and dwarfish phenotype reminiscent of iron-deficient plants. Ultrastructural analysis of plastids revealed severely impaired chloroplast development and a striking increase in ferritin clusters. Besides upregulation of ferritin, pic1 mutants showed differential regulation of genes and proteins related to iron stress or transport, photosynthesis, and Fe-S cluster biogenesis. Furthermore, PIC1 and its cyanobacterial homolog mediated iron accumulation in an iron uptake-defective yeast mutant. These observations suggest that PIC1 functions in iron transport across the inner envelope of chloroplasts and hence in cellular metal homeostasis.

Reaping the Benefits of SAGE

Serial analysis of gene expression (SAGE) is a powerful technique which yields a digital measure of gene expression through the sequencing of libraries of specific mRNA-derived fragments, namely SAGE tags. This chapter introduces the methods and software tools that are available for researchers to analyze gene expression through SAGE analysis. A detailed examination of SAGE analysis in Arabidopsis thaliana using the publicly available analysis tool, SaskSAGE, is provided. The use of this software allows the user to maximize the information gained from SAGE experiments in a model system with a fully sequenced genome.

Annotation and expression profile analysis of 2073 full-length cDNAs from stress-induced maize (Zea mays L.) seedlings

Full-length cDNAs are very important for genome annotation and functional analysis of genes. The number of full-length cDNAs from maize (Zea mays L.) remains limited. Here we report the construction of a full-length enriched cDNA library from osmotically stressed maize seedlings by using the modified CAP trapper method. From this library, 2073 full-length cDNAs were collected and further analyzed by sequencing from both the 5'- and 3'-ends. A total of 1728 (83.4%) sequences did not match known maize mRNA and full-length cDNA sequences in the GenBank database and represent new full-length genes. After alignment of the 2073 full-length cDNAs with 448 maize BAC sequences, it was found that 84 full-length cDNAs could be mapped to the BACs. Of these, 43 genes (51.2%) have been correctly annotated from the BAC clones, 37 genes (44.0%) have been annotated with a different exon-intron structure from our cDNA, and four genes (4.76%) had no annotations in the TIGR database. Expression analysis of 2073 full-length maize cDNAs using a cDNA macroarray led to the identification of 79 genes upregulated by stress treatments and 329 downregulated genes. Of the 79 stress-inducible genes, 30 genes contain ABRE, DRE, MYB, MYC core sequences or other abiotic-responsive cis-acting elements in their promoters. These results suggest that these cis-acting elements and the corresponding transcription factors take part in plant responses to osmotic stress either cooperatively or independently. Additionally, the data suggest that an ethylene signaling pathway may be involved in the maize response to drought stress.

ASEtrap: a biological method for speeding up the exploration of spliceomes

Alternative splicing (AS) of pre-messenger RNA is a major mechanism for generating protein diversity from a limited number of genes in higher eukaryotes, and it constitutes a central mode of genetic regulation. Thus, efficient methods are needed to systematically identify new AS events at a genomic scale across different tissues, stages of development, and physiological or pathological conditions in order to better understand gene expression. To fulfill this goal, we have designed the ASEtrap, which is a cloning procedure for producing AS libraries that is based on a single-stranded trap consisting of an ssDNA-binding protein. In this paper, we have applied our approach to the construction of an AS library and a Control library from human placenta. By analyzing 9226 and 9999 sequences of the AS and Control libraries, respectively, we show that internal AS events (events that can be identified by the sole resources provided by either the AS or the Control library) and the discovery rate of new AS events measured at early stages of sequencing were nine to 10 times higher in the former than in the latter. Moreover, by performing a search for new AS events within a group of 162 known drug target genes, we identified six new events in six genes, and we observed that they all were discovered exclusively through the AS library. Thus, it appears that ASEtrap has the potential to greatly facilitate the determination of the total complement of splice variants expressed in human, as well as other organisms.

Simultaneous high-throughput recombinational cloning of open reading frames in closed and open configurations

Comprehensive open reading frame (ORF) clone collections, ORFeomes, are key components of functional genomics projects. When recombinational cloning systems are used to capture ORFs in master clones, these DNA sequences can be easily transferred into a variety of expression plasmids, each designed for a specific assay. Depending on downstream applications, an ORF is cloned either with or without a stop codon at its original position, referred to as closed or open configuration, respectively. The former is preferred when the encoded protein is produced in its native form or with an amino-terminal tag; the latter is obligatory when the protein is produced as a fusion with a carboxyl-terminal tag. We developed a streamlined protocol for high-throughput, simultaneous cloning of both open and closed ORF entry clones with the Gateway recombinational cloning system. The protocol is straightforward to set up in large-scale ORF cloning projects, and is cost-effective, because the initial ORF amplification and the cloning in a pDONR vector are performed only once to obtain the two ORF configurations. We illustrated its implementation for the isolation and validation of 346 Arabidopsis ORF entry clones.

Comparative genomics as a tool for gene discovery

With the increasing availability of data from multiple eukaryotic genome sequencing projects, attention has focused on interspecific comparisons to discover novel genes and transcribed genomic sequences. Generally, these extrinsic strategies combine ab initio gene prediction with expression and/or homology data to identify conserved gene candidates between two or more genomes. Interspecific sequence analyses have proven invaluable for the improvement of existing annotations, automation of annotation, and identification of novel coding regions and splice variants. Further, comparative genomic approaches hold the promise of improved prediction of terminal or small exons, microRNA precursors, and small peptide-encoding open reading frames--sequence elements that are difficult to identify through purely intrinsic methodologies in the absence of experimental data.

Cloned sequence repertoires for small- and large-scale biology

How to assign function to the tens of thousands of genes discovered in the chromosomes of a few model species? How to complement the classical genetic approaches that are not always ideally suited to decode complex mechanisms? The solutions to these pressing questions are not simple and rely on the development of novel resources and technologies. Here I critically review what clone collections are available and how they can be exploited for the systematic analysis of gene functions in plants.

Fuzzy logic in medicine and bioinformatics

The purpose of this paper is to present a general view of the current applications of fuzzy logic in medicine and bioinformatics. We particularly review the medical literature using fuzzy logic. We then recall the geometrical interpretation of fuzzy sets as points in a fuzzy hypercube and present two concrete illustrations in medicine (drug addictions) and in bioinformatics (comparison of genomes).

Computational and comparative analyses of 150 full-length cDNA sequences from the oomycete plant pathogen Phytophthora infestans

Phytophthora infestans is a devastating phytopathogenic oomycete that causes late blight on tomato and potato. Recent genome sequencing efforts of P. infestans and other Phytophthora species are generating vast amounts of sequence data providing opportunities to unlock the complex nature of pathogenesis. However, accurate annotation of Phytophthora genomes will be a significant challenge. Most of the information about gene structure in these species was gathered from a handful of genes resulting in significant limitations for development of ab initio gene-calling programs. In this study, we collected a total of 150 bioinformatically determined near full-length cDNA (FLcDNA) sequences of P. infestans that were predicted to contain full open reading frame sequences. We performed detailed computational analyses of these FLcDNA sequences to obtain a snapshot of P. infestans gene structure, gauge the degree of sequence conservation between P. infestans genes and those of Phytophthora sojae and Phytophthora ramorum, and identify patterns of gene conservation between P. infestans and various eukaryotes, particularly fungi, for which genome-wide translated protein sequences are available. These analyses helped us to define the structural characteristics of P. infestans genes using a validated data set. We also determined the degree of sequence conservation within the genus Phytophthora and identified a set of fast evolving genes. Finally, we identified a set of genes that are shared between Phytophthora and fungal phytopathogens but absent in animal fungal pathogens. These results confirm that plant pathogenic oomycetes and fungi share virulence components, and suggest that eukaryotic microbial pathogens that share similar lifestyles also share a similar set of genes independently of their phylogenetic relatedness.