From hormones to secondary metabolism: the emergence of metabolic gene clusters in plants

Gene clusters for the synthesis of secondary metabolites are a common feature of microbial genomes. Well-known examples include clusters for the synthesis of antibiotics in actinomycetes, and also for the synthesis of antibiotics and toxins in filamentous fungi. Until recently it was thought that genes for plant metabolic pathways were not clustered, and this is certainly true in many cases; however, five plant secondary metabolic gene clusters have now been discovered, all of them implicated in synthesis of defence compounds. An obvious assumption might be that these eukaryotic gene clusters have arisen by horizontal gene transfer from microbes, but there is compelling evidence to indicate that this is not the case. This raises intriguing questions about how widespread such clusters are, what the significance of clustering is, why genes for some metabolic pathways are clustered and those for others are not, and how these clusters form. In answering these questions we may hope to learn more about mechanisms of genome plasticity and adaptive evolution in plants. It is noteworthy that for the five plant secondary metabolic gene clusters reported so far, the enzymes for the first committed steps all appear to have been recruited directly or indirectly from primary metabolic pathways involved in hormone synthesis. This may or may not turn out to be a common feature of plant secondary metabolic gene clusters as new clusters emerge.     

Origin and distribution of epipolythiodioxopiperazine (ETP) gene clusters in filamentous ascomycetes

Background  

Genes responsible for biosynthesis of fungal secondary metabolites are usually tightly clustered in the genome and co-regulated with metabolite production. Epipolythiodioxopiperazines (ETPs) are a class of secondary metabolite toxins produced by disparate ascomycete fungi and implicated in several animal and plant diseases. Gene clusters responsible for their production have previously been defined in only two fungi. Fungal genome sequence data have been surveyed for the presence of putative ETP clusters and cluster data have been generated from several fungal taxa where genome sequences are not available. Phylogenetic analysis of cluster genes has been used to investigate the assembly and heredity of these gene clusters.     

Gene duplications and the time thereafter – examples from plant secondary metabolism

Gene duplications are regarded as one of the central mechanisms for the origin of new genes. Recent studies in plant secondary metabolism have provided several examples of genes that originated by duplication with successive diversification. In this review, the mechanisms of gene duplication are explained and several models discussed that suggest the way that gene duplicates develop into genes with new functions. Signatures of gene duplication and diversification processes are discussed using the biosynthesis of benzoxazinones and of pyrrolizidine alkaloids as examples.     

Evolution of Chemical Diversity in Echinocandin Lipopeptide Antifungal Metabolites

The echinocandins are a class of antifungal drugs that includes caspofungin, micafungin, and anidulafungin. Gene clusters encoding most of the structural complexity of the echinocandins provided a framework for hypotheses about the evolutionary history and chemical logic of echinocandin biosynthesis. Gene orthologs among echinocandin-producing fungi were identified. Pathway genes, including the nonribosomal peptide synthetases (NRPSs), were analyzed phylogenetically to address the hypothesis that these pathways represent descent from a common ancestor. The clusters share cooperative gene contents and linkages among the different strains. Individual pathway genes analyzed in the context of similar genes formed unique echinocandin-exclusive phylogenetic lineages. The echinocandin NRPSs, along with the NRPS from the inp gene cluster in Aspergillus nidulans and its orthologs, comprise a novel lineage among fungal NRPSs. NRPS adenylation domains from different species exhibited a one-to-one correspondence between modules and amino acid specificity that is consistent with models of tandem duplication and subfunctionalization. Pathway gene trees and Ascomycota phylogenies are congruent and consistent with the hypothesis that the echinocandin gene clusters have a common origin. The disjunct Eurotiomycete-Leotiomycete distribution appears to be consistent with a scenario of vertical descent accompanied by incomplete lineage sorting and loss of the clusters from most lineages of the Ascomycota. We present evidence for a single evolutionary origin of the echinocandin family of gene clusters and a progression of structural diversification in two fungal classes that diverged approximately 290 to 390 million years ago. Lineage-specific gene cluster evolution driven by selection of new chemotypes contributed to diversification of the molecular functionalities.     

Localization of retina/pineal-expressed sequences: identification of novel candidate genes for inherited retinal disorders.

More than 100 genes causing inherited retinal diseases have been mapped to chromosomal locations, but less than half of these genes have been cloned. Mutations in many retina/pineal-specific genes are known to cause inherited retinal diseases. Examples include mutations in arrestin, rhodopsin kinase, and the cone-rod homeobox gene, CRX. To identify additional candidate genes for inherited retinal disorders, novel retina/pineal-expressed EST clusters were identified from the TIGR Human Gene Index database and mapped to specific chromosomal sites. After known human gene sequences were excluded, and repeat sequences were masked, 26 novel retina and pineal gland cDNA clusters were identified. The retinal expression of each novel EST cluster was confirmed by PCR assay of a retinal cDNA library, and each cluster was localized in the genome using the GeneBridge 4.0 radiation hybrid panel. In silico expression data from the TIGR database suggest that these EST clusters are retina/pineal-specific or predominantly expressed in these tissues. This combination of database analysis and laboratory investigation has localized several EST clusters that are potential candidates for genes causing inherited retinopathy.     

Chloroplast genomes of the diatoms <Emphasis Type="Italic">Phaeodactylum tricornutum</Emphasis> and <Emphasis Type="Italic">Thalassiosira pseudonana</Emphasis>: comparison with other plastid genomes of the red lineage

The chloroplast genomes of the pennate diatom Phaeodactylum tricornutum and the centric diatom Thalassiosira pseudonana have been completely sequenced and are compared with those of other secondary plastids of the red lineage: the centric diatom Odontella sinensis, the haptophyte Emiliania huxleyi, and the cryptophyte Guillardia theta. All five chromist genomes are compact, with small intergenic regions and no introns. The three diatom genomes are similar in gene content with 127-130 protein-coding genes, and genes for 27 tRNAs, three ribosomal RNAs and two small RNAs (tmRNA and signal recognition particle RNA). All three genomes have open-reading frames corresponding to ORFs148, 355 and 380 of O. sinensis, which have been assigned the names ycf88, ycf89 and ycf90. Gene order is not strictly conserved, but there are a number of conserved gene clusters showing remnants of red algal origin. The acpP, tsf and psb28 genes appear to be on the way from the plastid to the host nucleus, indicating that endosymbiotic gene transfer is a continuing process.     

Current strategies to induce secondary metabolites from microbial biosynthetic cryptic gene clusters

The genome of actinomycetes and several other microorganisms are endowed with many cryptic gene clusters that can code for previously undetected, a plethora of complex secondary metabolites. Under standard laboratory controlled conditions, the genes regulating these biosynthetic clusters are expressed at very low levels or remain phenotypically cryptic (silent). Over the past several decades, multi-drug-resistant bacteria have been observed with increased frequency, posing a significant threat to human health worldwide. The present alarming situation urgently calls for concerted global efforts for the discovery of new antimicrobials. The present situation, if not controlled, will take us again to the pre-antibiotic era. Today, in the post-genomic era, various new strategies such as the activation of cryptic gene clusters in microorganisms rejuvenate a new conviction in the field of natural product research that may lead to the identification of yet-unidentified novel secondary metabolites of therapeutic and other use. Decryptification of this versatile endogenous genetic reservoir may provide in the near future the more concrete rationale for antibiotic discovery. The present review is an attempt to provide a comprehensive detail, outlining current strategies that have been shown successful to activate cryptic biosynthetic gene clusters in microorganisms.     

Connected gene neighborhoods in prokaryotic genomes

A computational method was developed for delineating connected gene neighborhoods in bacterial and archaeal genomes. These gene neighborhoods are not typically present, in their entirety, in any single genome, but are held together by overlapping, partially conserved gene arrays. The procedure was applied to comparing the orders of orthologous genes, which were extracted from the database of Clusters of Orthologous Groups of proteins (COGs), in 31 prokaryotic genomes and resulted in the identification of 188 clusters of gene arrays, which included 1001 of 2890 COGs. These clusters were projected onto actual genomes to produce extended neighborhoods including additional genes, which are adjacent to the genes from the clusters and are transcribed in the same direction, which resulted in a total of 2387 COGs being included in the neighborhoods. Most of the neighborhoods consist predominantly of genes united by a coherent functional theme, but also include a minority of genes without an obvious functional connection to the main theme. We hypothesize that although some of the latter genes might have unsuspected roles, others are maintained within gene arrays because of the advantage of expression at a level that is typical of the given neighborhood. We designate this phenomenon ‘genomic hitchhiking’. The largest neighborhood includes 79 genes (COGs) and consists of overlapping, rearranged ribosomal protein superoperons; apparent genome hitchhiking is particularly typical of this neighborhood and other neighborhoods that consist of genes coding for translation machinery components. Several neighborhoods involve previously undetected connections between genes, allowing new functional predictions. Gene neighborhoods appear to evolve via complex rearrangement, with different combinations of genes from a neighborhood fixed in different lineages.     

Beyond aflatoxin: four distinct expression patterns and functional roles associated with Aspergillus flavus secondary metabolism gene clusters

Species of Aspergillus produce a diverse array of secondary metabolites, and recent genomic analysis has predicted that these species have the capacity to synthesize many more compounds. It has been possible to infer the presence of 55 gene clusters associated with secondary metabolism in Aspergillus flavus ; however, only three metabolic pathways—aflatoxin, cyclopiazonic acid (CPA) and aflatrem—have been assigned to these clusters. To gain an insight into the regulation of and to infer the ecological significance of the 55 secondary metabolite gene clusters predicted in A. flavus, we examined their expression over 28 diverse conditions. Variables included culture medium and temperature, fungal development, colonization of developing maize seeds and misexpression of laeA , a global regulator of secondary metabolism. Hierarchical clustering analysis of expression profiles allowed us to categorize the gene clusters into four distinct clades. Gene clusters for the production of aflatoxins, CPA and seven other unknown compound(s) were identified as belonging to one clade. To further explore the relationships found by gene expression analysis, aflatoxin and CPA production were quantified under five different cell culture environments known to be conducive or nonconducive for aflatoxin biosynthesis and during the colonization of developing maize seeds. Results from these studies showed that secondary metabolism gene clusters have distinctive gene expression profiles. Aflatoxin and CPA were found to have unique regulation, but are sufficiently similar that they would be expected to co-occur in substrates colonized with A. flavus .     

Analysis of homologous gene clusters in Caenorhabditis elegans reveals striking regional cluster domains

An algorithm for detecting local clusters of homologous genes was applied to the genome of Caenorhabditis elegans. Clusters of two or more homologous genes are abundant, totaling 1391 clusters containing 4607 genes, over one-fifth of all genes in C. elegans. Cluster genes are distributed unevenly in the genome, with the large majority located on autosomal chromosome arms, regions characterized by higher genetic recombination and more repeat sequences than autosomal centers and the X chromosome. Cluster genes are transcribed at much lower levels than average and very few have gross phenotypes as assayed by RNAi-mediated reduction of function. The molecular identity of cluster genes is unusual, with a preponderance of nematode-specific gene families that encode putative secreted and transmembrane proteins, and enrichment for genes implicated in xenobiotic detoxification and innate immunity. Gene clustering in Drosophila melanogaster is also substantial and the molecular identity of clustered genes follows a similar pattern. I hypothesize that autosomal chromosome arms in C. elegans undergo frequent local gene duplication and that these duplications support gene diversification and rapid evolution in response to environmental challenges. Although specific gene clusters have been documented in C. elegans, their abundance, genomic distribution, and unusual molecular identities were previously unrecognized.     

Complete genome sequence of Nitrosospira multiformis, an ammonia-oxidizing bacterium from the soil environment

The complete genome of the ammonia-oxidizing bacterium Nitrosospira multiformis (ATCC 25196(T)) consists of a circular chromosome and three small plasmids totaling 3,234,309 bp and encoding 2,827 putative proteins. Of the 2,827 putative proteins, 2,026 proteins have predicted functions and 801 are without conserved functional domains, yet 747 of these have similarity to other predicted proteins in databases. Gene homologs from Nitrosomonas europaea and Nitrosomonas eutropha were the best match for 42% of the predicted genes in N. multiformis. The N. multiformis genome contains three nearly identical copies of amo and hao gene clusters as large repeats. The features of N. multiformis that distinguish it from N. europaea include the presence of gene clusters encoding urease and hydrogenase, a ribulose-bisphosphate carboxylase/oxygenase-encoding operon of distinctive structure and phylogeny, and a relatively small complement of genes related to Fe acquisition. Systems for synthesis of a pyoverdine-like siderophore and for acyl-homoserine lactone were unique to N. multiformis among the sequenced genomes of ammonia-oxidizing bacteria. Gene clusters encoding proteins associated with outer membrane and cell envelope functions, including transporters, porins, exopolysaccharide synthesis, capsule formation, and protein sorting/export, were abundant. Numerous sensory transduction and response regulator gene systems directed toward sensing of the extracellular environment are described. Gene clusters for glycogen, polyphosphate, and cyanophycin storage and utilization were identified, providing mechanisms for meeting energy requirements under substrate-limited conditions. The genome of N. multiformis encodes the core pathways for chemolithoautotrophy along with adaptations for surface growth and survival in soil environments.     

Phylogenomic analyses of KCNA gene clusters in vertebrates: why do gene clusters stay intact?

Background  

Gene clusters are of interest for the understanding of genome evolution since they provide insight in large-scale duplications events as well as patterns of individual gene losses. Vertebrates tend to have multiple copies of gene clusters that typically are only single clusters or are not present at all in genomes of invertebrates. We investigated the genomic architecture and conserved non-coding sequences of vertebrate KCNA gene clusters. KCNA genes encode shaker-related voltage-gated potassium channels and are arranged in two three-gene clusters in tetrapods. Teleost fish are found to possess four clusters. The two tetrapod KNCA clusters are of approximately the same age as the Hox gene clusters that arose through duplications early in vertebrate evolution. For some genes, their conserved retention and arrangement in clusters are thought to be related to regulatory elements in the intergenic regions, which might prevent rearrangements and gene loss. Interestingly, this hypothesis does not appear to apply to the KCNA clusters, as too few conserved putative regulatory elements are retained.     

Two plus two does not equal three: statistical tests for multiple genome comparison

Gene clusters that span three or more chromosomal regions are of increasing importance, yet statistical tests to validate such clusters are in their infancy. Current approaches either conduct several pairwise comparisons or consider only the number of genes that occur in all of the regions. In this paper, we provide statistical tests for clusters spanning exactly three regions based on genome models of typical comparative genomics problems, including analysis of conserved linkage within multiple species and identification of large-scale duplications. Our tests are the first to combine evidence from genes shared among all three regions and genes shared between pairs of regions. We show that our tests of clusters spanning three regions are more sensitive than existing approaches, and can thus be used to identify more diverged homologous regions.     

Identifying clusters of functionally related genes in genomes

MOTIVATION: An increasing body of literature shows that genomes of eukaryotes can contain clusters of functionally related genes. Most approaches to identify gene clusters utilize microarray data or metabolic pathway databases to find groups of genes on chromosomes that are linked by common attributes. A generalized method that can find gene clusters regardless of the mechanism of origin would provide researchers with an unbiased method for finding clusters and studying the evolutionary forces that give rise to them. RESULTS: We present an algorithm to identify gene clusters in eukaryotic genomes that utilizes functional categories defined in graph-based vocabularies such as the Gene Ontology (GO). Clusters identified in this manner need only have a common function and are not constrained by gene expression or other properties. We tested the algorithm by analyzing genomes of a representative set of species. We identified species-specific variation in percentage of clustered genes as well as in properties of gene clusters including size distribution and functional annotation. These properties may be diagnostic of the evolutionary forces that lead to the formation of gene clusters. AVAILABILITY: A software implementation of the algorithm and example output files are available at http://fcg.tamu.edu/C_Hunter/.