豌豆蚜基因组P450基因家族的分析

细胞色素P450在各种内源和外源物质代谢中起着非常重要的作用。利用豌豆蚜Acyrthosiphon pisum基因组mRNA和氨基酸数据库研究P450基因功能的进化规律。通过生物信息学方法对豌豆蚜全基因组P450进行分析, 结果显示: 在豌豆蚜基因组中发现69个P450基因, 它们分别属于13个P450家族和18个亚家族, 是一个典型的多基因家族。进一步将这些基因与豌豆蚜ESTs数据库进行了比对分析, 其中39个候选基因有EST证据, 证明了这些P450基因在转录水平的真实性。以氨基酸相似度大于60%为标准对豌豆蚜基因组中P450基因进行分组, 69个P450基因中, 除18个基因序列因差异太大, 不能被归入任何一组, 其余51个可归入10个组, 其中8个组(包含47条序列)适合于正选择和基因转换分析。正选择和基因转换分析结果表明: 仅有1个组(含9个基因)显著受到正选择压力作用, 正选择概率大于95%的氨基酸位点分别是20T和27N, 20T位于底物识别位点SRS1, 27N位于D. helix; 有3个组(包含8个基因）显示显著的基因转换事件。 参与基因转换的基因均为CYP4家族成员, 分别是CYP4C, CYP4G和CYP4V亚家族。 参与基因转换的成员之间的蛋白相似度较高（70％~95％）, 且XM_001944991与XM_001951794同存在于SCAFFOLD12542上, XM_001945510与XM_001944057同存在于SCAFFOLD7010上。这可能暗示豌豆蚜P450基因通过基因复制, 然后通过基因转换使P450获得新的功能, 以适应多变的生存环境。此外, 鉴定出20个不同的基序, 其中有5条基序在90％以上的基因中出现。      

Comparison of Cytochrome P450 Genes from Six Plant Genomes

Plants depend on cytochrome P450 (CYP) enzymes for nearly every aspect of their biology. In several sequenced angiosperms, CYP genes constitute up to 1% of the protein coding genes. The angiosperm sequence diversity is encapsulated by 59 CYP families, of which 52 families form a widely distributed core set. In the 20 years since the first plant P450 was sequenced, 3,387 P450 sequences have been identified and annotated in plant databases. As no new angiosperm CYP families have been discovered since 2004, it is now apparent that the sampling of CYP diversity is beginning to plateau. This review presents a comparison of 1,415 cytochrome P450 sequences from the six sequenced genomes of Vitis vinifera (grape), Carica papaya (papaya), Populus trichocarpa (poplar), Oryza sativa (rice), Arabidopsis thaliana (Arabidopsis or mouse ear’s cress) and Physcomitrella patens (moss). An evolutionary analysis is presented that tracks land plant P450 innovation over time from the most ancient and conserved sequences to the newest dicot-specific families. The earliest or oldest P450 families are devoted to the essential biochemistries of sterol and carotenoid synthesis. The next evolutionary radiation of P450 families appears to mediate crucial adaptations to a land environment. And, the newest CYP families appear to have driven the diversity of angiosperms in mediating the synthesis of pigments, odorants, flavors and order-/genus-specific secondary metabolites. Family-by-family comparisons allow the visualization of plant genome plasticity by whole genome duplications and massive gene family expansions via tandem duplications. Molecular evidence of human domestication is quite apparent in the repeated P450 gene duplications occurring in the grape genome.     

Genome Sequence of the Pea Aphid Acyrthosiphon pisum

Aphids are important agricultural pests and also biological models for studies of insect-plant interactions, symbiosis, virus vectoring, and the developmental causes of extreme phenotypic plasticity. Here we present the 464 Mb draft genome assembly of the pea aphid Acyrthosiphon pisum. This first published whole genome sequence of a basal hemimetabolous insect provides an outgroup to the multiple published genomes of holometabolous insects. Pea aphids are host-plant specialists, they can reproduce both sexually and asexually, and they have coevolved with an obligate bacterial symbiont. Here we highlight findings from whole genome analysis that may be related to these unusual biological features. These findings include discovery of extensive gene duplication in more than 2000 gene families as well as loss of evolutionarily conserved genes. Gene family expansions relative to other published genomes include genes involved in chromatin modification, miRNA synthesis, and sugar transport. Gene losses include genes central to the IMD immune pathway, selenoprotein utilization, purine salvage, and the entire urea cycle. The pea aphid genome reveals that only a limited number of genes have been acquired from bacteria; thus the reduced gene count of Buchnera does not reflect gene transfer to the host genome. The inventory of metabolic genes in the pea aphid genome suggests that there is extensive metabolite exchange between the aphid and Buchnera, including sharing of amino acid biosynthesis between the aphid and Buchnera. The pea aphid genome provides a foundation for post-genomic studies of fundamental biological questions and applied agricultural problems.     

Intron-exon organization and phylogeny in a large superfamily, the paralogous cytochrome P450 genes of Arabidopsis thaliana

The cytochrome P450 gene superfamily is represented by 80 genes in animal genomes and perhaps more than 300 genes in plant genomes. We analyzed about half of all Arabidopsis P450 genes, a very large dataset of truly paralogous genes. Sequence alignments were used to draw phylogenetic trees, and this information was compared with the intron-exon organization of each P450 gene. We found 60 unique intron positions, of which 37 were phase 0 introns. Our results confirm the polyphyletic origin of plant P450 genes. One group of these genes, the A-type P450s, are plant specific and characterized by a simple organization, with one highly conserved intron. Closely related A-type P450 genes are often clustered in the genome with as many as a dozen genes (e.g., of the CYP71 subfamily) on a short stretch of chromosome. The other P450 genes (non-A-type) form several distinct clades and are characterized by numerous introns. One such clade contains the two CYP51 genes, which are thought to encode obtusifoliol 14a demethylase. The two CYP51 genes have a single intron that is not shared with CYP51 genes from vertebrates or fungi, or with any other Arabidopsis P450 gene. Only a few of the Arabidopsis P450 genes are intronless (e.g., the CYP710A and CYP96A subfamilies). There was a relatively good correlation between intron conservation and phylogenetic relationships between members of the P450 subfamilies. Gene organization appears to be a useful tool in establishing the evolutionary relatedness of P450 genes, which may help in predictions of P450 function.     

Prediction of signal recognition particle RNA genes

We describe a method for prediction of genes that encode the RNA component of the signal recognition particle (SRP). A heuristic search for the strongly conserved helix 8 motif of SRP RNA is combined with covariance models that are based on previously known SRP RNA sequences. By screening available genomic sequences we have identified a large number of novel SRP RNA genes and we can account for at least one gene in every genome that has been completely sequenced. Novel bacterial RNAs include that of Thermotoga maritima, which, unlike all other non-gram-positive eubacteria, is predicted to have an Alu domain. We have also found the RNAs of Lactococcus lactis and Staphylococcus to have an unusual UGAC tetraloop in helix 8 instead of the normal GNRA sequence. An investigation of yeast RNAs reveals conserved sequence elements of the Alu domain that aid in the analysis of these RNAs. Analysis of the human genome reveals only two likely genes, both on chromosome 14. Our method for SRP RNA gene prediction is the first convenient tool for this task and should be useful in genome annotation.     

The cytochrome P450 genes of channel catfish: Their involvement in disease defense responses as revealed by meta-analysis of RNA-Seq data sets

Background

Cytochrome P450s (CYPs) encode one of the most diverse enzyme superfamily in nature. They catalyze oxidative reactions of endogenous molecules and exogenous chemicals.

Methods

We identified CYPs genes through in silico analysis using EST, RNA-Seq and genome databases of channel catfish. Phylogenetic analyses and conserved syntenic analyses were conducted to determine their identities and orthologies. Meta-analysis of RNA-Seq databases was conducted to analyze expression profile of CYP genes following bacterial infection.

Results

A full set of 61 CYP genes was identified and characterized in channel catfish. Phylogenetic tree and conserved synteny provided strong evidence of their identities and orthorlogy. Lineage-specific gene duplication was evident in a number of clans in channel catfish. CYP46A1 is missing in the catfish genome as observed with syntenic analysis and RT-PCR analysis. Thirty CYPs were found up- or down-regulated in liver, while seven and eight CYPs were observed regulated in intestine and gill following bacterial infection.

Conclusion

We systematically identified and characterized a full set of 61 CYP genes in channel catfish and studied their expression profiles after bacterial infection. While bacterial challenge altered the expression of large numbers of CYP genes, the mechanisms and significance of these changes are not known.

General significance

This work provides an example to systematically study CYP genes in non-model species. Moreover, it provides a basis for further toxicological and physiological studies in channel catfish.     

Structural organization and classification of cytochrome P450 genes in flax (Linum usitatissimum L.)

Flax CYPome analysis resulted in the identification of 334 putative cytochrome P450 (CYP450) genes in the cultivated flax genome. Classification of flax CYP450 genes based on the sequence similarity with Arabidopsis orthologs and CYP450 nomenclature, revealed 10 clans representing 44 families and 98 subfamilies. CYP80, CYP83, CYP92, CYP702, CYP705, CYP708, CYP728, CYP729, CYP733 and CYP736 families are absent in the flax genome. The subfamily members exhibited conserved sequences, length of exons and phasing of introns. Similarity search of the genomic resources of wild flax species Linum bienne with CYP450 coding sequences of the cultivated flax, revealed the presence of 127 CYP450 gene orthologs, indicating amplification of novel CYP450 genes in the cultivated flax. Seven families CYP73, 74, 75, 76, 77, 84 and 709, coding for enzymes associated with phenylpropanoid/fatty acid metabolism, showed extensive gene amplification in the flax. About 59% of the flax CYP450 genes were present in the EST libraries.     

Genome-wide analysis of cytochrome P450 monooxygenase genes in the silkworm, Bombyx mori

Based on the advances in the silkworm genome project, a new genome-wide analysis of cytochrome P450 genes was performed. A total of 84 CYP-related sequences were identified and could be classified into 26 families and 47 subfamilies according to standard nomenclature. Seventy eight of the eighty four genes appear to be functional and six are probable pseudogenes. The distribution of Bombyx mori P450s in the genome shows that most of them are tandem arranged on chromosomes, only 34 genes are present as singletons, with 8 clusters including 3 or more than 3 genes. Sequence alignments were used to reconstruct phylogenetic trees and to analyze the intron-exon organizations of the functional genes. The conserved intron positioning agrees perfectly with their common grouping on the tree. The presence of three extremely ancient introns which are conserved across different clans indicates that a few introns are still highly conserved after they have undergone extensive evolutionary changes of B. mori P450 duplication and divergence. Comparison of the P450s from B. mori to the P450s from Drosophila melanogaster shows that the expansion is not uniform across the gene families. Remarkably, two mitochondrial families, the B. mori CYP333 and D. melanogaster Cyp12, formed two orthologous groups in the phylogenetic tree. All CYP333s can be proposed to be related to xenobiotic metabolism in accordance with the D. melanogaster Cyp12s. The characterization and evolutionary analysis of P450s from B. mori in the current study provide useful information for understanding the characteristics and diversity of P450s from B. mori and the baseline for functional analyses of individual P450s in this model Lepidopteran insect.     

The Complete Mitochondrial Genome of Gossypium hirsutum and Evolutionary Analysis of Higher Plant Mitochondrial Genomes

Background

Mitochondria are the main manufacturers of cellular ATP in eukaryotes. The plant mitochondrial genome contains large number of foreign DNA and repeated sequences undergone frequently intramolecular recombination. Upland Cotton (Gossypium hirsutum L.) is one of the main natural fiber crops and also an important oil-producing plant in the world. Sequencing of the cotton mitochondrial (mt) genome could be helpful for the evolution research of plant mt genomes.

Methodology/Principal Findings

We utilized 454 technology for sequencing and combined with Fosmid library of the Gossypium hirsutum mt genome screening and positive clones sequencing and conducted a series of evolutionary analysis on Cycas taitungensis and 24 angiosperms mt genomes. After data assembling and contigs joining, the complete mitochondrial genome sequence of G. hirsutum was obtained. The completed G.hirsutum mt genome is 621,884 bp in length, and contained 68 genes, including 35 protein genes, four rRNA genes and 29 tRNA genes. Five gene clusters are found conserved in all plant mt genomes; one and four clusters are specifically conserved in monocots and dicots, respectively. Homologous sequences are distributed along the plant mt genomes and species closely related share the most homologous sequences. For species that have both mt and chloroplast genome sequences available, we checked the location of cp-like migration and found several fragments closely linked with mitochondrial genes.

Conclusion

The G. hirsutum mt genome possesses most of the common characters of higher plant mt genomes. The existence of syntenic gene clusters, as well as the conservation of some intergenic sequences and genic content among the plant mt genomes suggest that evolution of mt genomes is consistent with plant taxonomy but independent among different species.     

Genome-wide identification and characterization of putative cytochrome P450 genes in the model legume <Emphasis Type="Italic">Medicago truncatula</Emphasis>

In plants, cytochrome P450 is a group of monooxygenases existing as a gene superfamily and plays important roles in metabolizing physiologically important compounds. However, to date only a limited number of P450s have been identified and characterized in legumes. In this study, data mining methods were used, and 151 putative P450 genes in the model legume Medicago truncatula were identified, including 135 novel sequences. These genes were classified into 9 clans and 44 families by sequence similarity, and among those 4 new clans and 21 new families not reported previously in legumes. By comparison of these genes with P450 genes in Arabidopsis and rice, it was found that most of the known P450 families in dicot species exist in M. truncatula. The representative protein sequences of putative P450s were aligned, and the secondary elements were assigned based on the known structure P450BM3. Putative substrate recognition sites (SRSs) and substrate binding sites were also identified in these sequences. In addition, the ESTs-derived expression profiles (digital Northern) of the putative P450 genes were analyzed, which was confirmed by semi-quantitative RT-PCR analyses of several selected P450 genes. These results will provide a base for catalogue information on P450 genes in M. truncatula and for further functional analysis of P450 superfamily genes in legumes.     

Cytochrome P450 Genes from the Sacred Lotus Genome

Cytochrome P450 monooxygenases (P450s) in the sacred lotus (Nelumbo nucifera) genome have been identified and named according to systematic P450 nomenclatures. Comparisons of these sequences with those in the papaya and grape CYPomes have indicated that gene blooms exist in the CYP89, CYP94, CYP96 and CYP714 families and that less dramatic expansions exist in the CYP71 and CYP72 families. Expansions in the CYP94 and CYP96 families may be associated with generation of the extremely hydrophobic leaf surfaces associated with the “lotus effect” in this water-adapted species, since these families are known to hydroxylate fatty acids and alkanes in the wax biosynthetic pathways of other plant species. Evolution of the CYP719 and CYP80 families may be associated with production of a number of benzylisoquinoline and aporphine alkaloids. Structures for anonaine and roemerine, two of the most abundant aporphine alkaloids in lotus leaves and seeds, contain methylenedioxy bridges that are known to be generated by members of the CYP719 family. With only one CYP719A22 gene existing in the lotus genome, it is likely that it is involved in making aporphine alkaloids. The fact that CYP719 has not previously been seen in angiosperm phylogeny below the order of Ranunculales suggests that its presence in lotus (in the Proteales) presents an evolutionary terminus prior to its loss in more recent eudicot species. With several CYP80 family genes existing in the lotus genome, there are multiple candidates for those involved in conducting benzylisoquinoline alkaloid synthesis.     

Origins of P450 diversity

The P450 enzymes maintain a conserved P450 fold despite a considerable variation in sequence. The P450 family even includes proteins that lack the single conserved cysteine and are therefore no longer haem-thiolate proteins. The mechanisms of successive gene duplications leading to large families in plants and animals are well established. Comparisons of P450 CYP gene clusters in related species illustrate the rapid changes in CYPome sizes. Examples of CYP copy number variation with effects on fitness are emerging, and these provide an opportunity to study the proximal causes of duplication or pseudogenization. Birth and death models can explain the proliferation of CYP genes that is amply illustrated by the sequence of every new genome. Thus, the distribution of P450 diversity within the CYPome of plants and animals, a few families with many genes (P450 blooms) and many families with few genes, follows similar power laws in both groups. A closer look at some families with few genes shows that these, often single member families, are not stable during evolution. The enzymatic prowess of P450 may predispose them to switch back and forth between metabolism of critical structural or signal molecules and metabolism dedicated to environmental response.     

棉花DUR3基因的鉴定及进化分析

植物DUR3同源蛋白属于钠离子/溶质共运蛋白家族的尿素高亲和力运输蛋白，在植物体对外源尿素的主动吸收及内源尿素的再分配过程中具有重要作用。为明确棉花DUR3基因的结构和进化情况，基于生物信息学的方法，从全基因组水平鉴定陆地棉和雷蒙德氏棉的DUR3基因，并对基因结构、跨膜结构域、基序分布、进化关系等进行分析。结果表明：(1)从陆地棉A亚组和D亚组染色体各鉴定出1个DUR3基因，从雷蒙德氏棉基因组鉴定出1个DUR3基因。这3个棉花DUR3同源蛋白同其他植物DUR3同源蛋白一样，具有15个跨膜结构域，具有3个位置一致、高度保守的基序。(2)基因结构分析表明，双子叶植物DUR3基因的外显子个数明显多于单子叶植物，这3个棉花DUR3基因的外显子个数亦是如此。(3)根据物种间种属亲缘关系，对不同物种DUR3氨基酸序列构建的进化树显示，棉花的同双子叶植物的聚在一起。(4)DUR3直系同源基因和旁系同源基因的K_a/K_s比值普遍均大于1,说明这些基因在进化过程中主要受到正向选择的作用。该研究结果为深入研究棉花DUR3同源蛋白提供了理论基础。     

Genome-wide analysis of simple sequence repeats in the model medicinal mushroom Ganoderma lucidum

Simple sequence repeats (SSRs) or microsatellites are one of the most popular sources of genetic markers and play a significant role in gene function and genome organization. We identified SSRs in the genome of Ganoderma lucidum and analyzed their frequency and distribution in different genomic regions. We also compared the SSRs in G. lucidum with six other Agaricomycetes genomes: Coprinopsis cinerea, Laccaria bicolor, Phanerochaete chrysosporium, Postia placenta, Schizophyllum commune and Serpula lacrymans. Based on our search criteria, the total number of SSRs found ranged from 1206 to 6104 and covered from 0.04% to 0.15% of the fungal genomes. The SSR abundance was not correlated with the genome size, and mono- to tri-nucleotide repeats outnumbered other SSR categories in all of the species examined. In G. lucidum, a repertoire of 2674 SSRs was detected, with mono-nucleotides being the most abundant. SSRs were found in all genomic regions and were more abundant in non-coding regions than coding regions. The highest SSR relative abundance was found in introns (108 SSRs/Mb), followed by intergenic regions (84 SSRs/Mb). A total of 684 SSRs were found in the protein-coding sequences (CDSs) of 588 gene models, with 81.4% of them being tri- or hexa-nucleotides. After scanning for InterPro domains, 280 of these genes were successfully annotated, and 215 of them could be assigned to Gene Ontology (GO) terms. SSRs were also identified in 28 bioactive compound synthesis-related gene models, including one 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR), three polysaccharide biosynthesis genes and 24 cytochrome P450 monooxygenases (CYPs). Primers were designed for the identified SSR loci, providing the basis for the future development of SSR markers of this medicinal fungus.     

CYP52 (cytochrome P450alk) multigene family in Candida maltosa: molecular cloning and nucleotide sequence of the two tandemly arranged genes

Southern blot analysis under low-stringency conditions using a previously isolated n-alkane-inducible cytochrome P450 (P450alk) gene as a probe revealed the presence of multiple P450alk-related genes in the genome of Candida maltosa. Nine P450alk-related genes (one reported previously and eight in the present report) were isolated from a genomic library constructed from this strain, and these were classified on the basis of sequence similarities into three pairs of putative allelic genes and three nonallelic genes. Two pairs of these alleles were tandemly arranged in the genome. The complete nucleotide sequences of one of these pairs were determined and compared to other members of this P450 family (CYP52) in C. maltosa and C. tropicalis. Northern blot analysis further showed that these genes were regulated by carbon sources. These results provide evidence for a P450alk (CYP52) multigene family in C. maltosa.