Propionigenium and Vibrio species identified as possible component causes of shrimp white feces syndrome (WFS) associated with the microsporidian Enterocytozoon hepatopenaei.

White feces syndrome (WFS) in cultivated shrimp is characterized by white shrimp midguts (intestines) and white fecal strings that float as mats on pond surfaces. The etiology of WFS is complex, but one type called EHP-WFS is associated with the microsporidian Enterocytozoon hepatopenaei (EHP). The hepatopancreas (HP), midgut and fecal strings of EHP-WFS shrimp exhibit massive quantities of EHP spores together with mixed, unidentified bacteria. In EHP-WFS ponds, some EHP-infected shrimp show white midguts (WG) and produce white feces while other EHP-infected shrimp in the same pond show grossly normal midguts (NG) and produce no white feces. We hypothesized that comparison of the microbial flora between WG and NG shrimp would reveal probable combinations of microbes significantly associated with EHP-WFS. To test this, we selected a Penaeus vannamei cultivation pond exhibiting severe WFS and used microscopic and microbial profiling analyses to compare WG and NG samples. Histologically, EHP was confirmed in the HP and midgut of both WG and NG shrimp, but EHP burdens were higher and EHP tissue damage was more severe in WG shrimp. Further, intestinal microbiomes in WG shrimp were less diverse and had higher abundance of bacteria from the genera Vibrio and Propionigenium. Propionigenium burden in the HP of WG shrimp (9364 copies/100 ng DNA) was significantly higher (P = 1.1 × 10-5) than in NG shrimp (12 copies/100 ng DNA). These findings supported our hypothesis by revealing two candidate bacterial genera that should be tested in combination with EHP as potential component causes of EHP-WFS in P. vannamei.

Decay of the glycolytic pathway and adaptation to intranuclear parasitism within Enterocytozoonidae microsporidia.

Glycolysis and oxidative phosphorylation are the fundamental pathways of ATP generation in eukaryotes. Yet in microsporidia, endoparasitic fungi living at the limits of cellular streamlining, oxidative phosphorylation has been lost: energy is obtained directly from the host or, during the dispersive spore stage, via glycolysis. It was therefore surprising when the first sequenced genome from the Enterocytozoonidae - a major family of human and animal-infecting microsporidians - appeared to have lost genes for glycolysis. Here, we sequence and analyse genomes from additional members of this family, shedding new light on their unusual biology. Our survey includes the genome of Enterocytozoon hepatopenaei, a major aquacultural parasite currently causing substantial economic losses in shrimp farming, and Enterospora canceri, a pathogen that lives exclusively inside epithelial cell nuclei of its crab host. Our analysis of gene content across the clade suggests that Ent. canceri's adaptation to intranuclear life is underpinned by the expansion of transporter families. We demonstrate that this entire lineage of pathogens has lost glycolysis and, uniquely amongst eukaryotes, lacks any obvious intrinsic means of generating energy. Our study provides an important resource for the investigation of host-pathogen interactions and reductive evolution in one of the most medically and economically important microsporidian lineages.

ShrimpGPAT: a gene and protein annotation tool for knowledge sharing and gene discovery in shrimp.

BACKGROUND: Although captured and cultivated marine shrimp constitute highly important seafood in terms of both economic value and production quantity, biologists have little knowledge of the shrimp genome and this partly hinders their ability to improve shrimp aquaculture. To help improve this situation, the Shrimp Gene and Protein Annotation Tool (ShrimpGPAT) was conceived as a community-based annotation platform for the acquisition and updating of full-length complementary DNAs (cDNAs), Expressed Sequence Tags (ESTs), transcript contigs and protein sequences of penaeid shrimp and their decapod relatives and for in-silico functional annotation and sequence analysis. DESCRIPTION: ShrimpGPAT currently holds quality-filtered, molecular sequences of 14 decapod species (~500,000 records for six penaeid shrimp and eight other decapods). The database predominantly comprises transcript sequences derived by both traditional EST Sanger sequencing and more recently by massive-parallel sequencing technologies. The analysis pipeline provides putative functions in terms of sequence homologs, gene ontologies and protein-protein interactions. Data retrieval can be conducted easily either by a keyword text search or by a sequence query via BLAST, and users can save records of interest for later investigation using tools such as multiple sequence alignment and BLAST searches against pre-defined databases. In addition, ShrimpGPAT provides space for community insights by allowing functional annotation with tags and comments on sequences. Community-contributed information will allow for continuous database enrichment, for improvement of functions and for other aspects of sequence analysis. CONCLUSIONS: ShrimpGPAT is a new, free and easily accessed service for the shrimp research community that provides a comprehensive and up-to-date database of quality-filtered decapod gene and protein sequences together with putative functional prediction and sequence analysis tools. An important feature is its community-based functional annotation capability that allows the research community to contribute knowledge and insights about the properties of molecular sequences for better, shared, functional characterization of shrimp genes. Regularly updated and expanded with data on more decapods, ShrimpGPAT is publicly available at http://shrimpgpat.sc.mahidol.ac.th/.

Shrimp injection with dsRNA targeting the microsporidian EHP polar tube protein reduces internal and external parasite amplification.

The microsporidian Enterocytozoon hepatopenaei (EHP) is a major threat to shrimp health worldwide. Severe EHP infections in shrimp cause growth retardation and increase susceptibility to opportunistic infections. EHP produces spores with a chitin wall that enables them to survive prolonged environmental exposure. Previous studies showed that polar tube extrusion is a prerequisite for EHP infection, such that inhibiting extrusion should prevent infection. Using a proteomic approach, polar tube protein 2 of EHP (EhPTP2) was found abundantly in protein extracts obtained from extruded spores. Using an immunofluorescent antibody against EhPTP2 for immunohistochemistry, extruded spores were found in the shrimp hepatopancreas (HP) and intestine, but not in the stomach. We hypothesized that presence of EhPTP2 might be required for successful EHP spore extrusion. To test this hypothesis, we injected EhPTP2-specific double-stranded RNA (dsRNA) and found that it significantly diminished EHP copy numbers in infected shrimp. This indicated reduced amplification of EHP-infected cells in the HP by spores released from previously infected cells. In addition, injection of the dsRNA into EHP-infected shrimp prior to their use in cohabitation with naïve shrimp significantly (p < 0.05) reduced the rate of EHP transmission to naïve shrimp. The results revealed that EhPTP2 plays a crucial role in the life cycle of EHP and that dsRNA targeting EHP mRNA can effectively reach the parasite developing in host cells. This approach is a model for future investigations to identify critical genes for EHP survival and spread as potential targets for preventative and therapeutic measures in shrimp.

Molecular and cellular characterization of four putative nucleotide transporters from the shrimp microsporidian Enterocytozoon hepatopenaei (EHP).

Microsporidia are obligate intracellular parasites that lost several enzymes required in energy production. The expansion of transporter families in these organisms enables them to hijack ATP from hosts. In this study, nucleotide transporters of the microsporidian Enterocytozoon hepatopenaei (EHP), which causes slow growth in economically valuable Penaeus shrimp, were characterized. Analysis of the EHP genome suggested the presence of four putative nucleotide transporter genes, namely EhNTT1, EhNTT2, EhNTT3, and EhNTT4. Sequence alignment revealed four charged amino acids that are conserved in previously characterized nucleotide transporters. Phylogenetic analysis suggested that EhNTT1, 3, and 4 were derived from one horizontal gene transfer event, which was independent from that of EhNTT2. Localization of EhNTT1 and EhNTT2 using immunofluorescence analysis revealed positive signals within the envelope of developing plasmodia and on mature spores. Knockdown of EhNTT2 by double administration of sequence specific double-stranded RNA resulted in a significant reduction in EHP copy numbers, suggesting that EhNTT2 is crucial for EHP replication in shrimp. Taken together, the insight into the roles of NTTs in microsporidian proliferation can provide the biological basis for the development of alternative control strategies for microsporidian infection in shrimp.

Incidence, genetic diversity, and antimicrobial resistance profiles of Vibrio parahaemolyticus in seafood in Bangkok and eastern Thailand.

Background: Emergence of Vibrio parahaemolyticus pandemic strain O3:K6 was first documented in 1996. Since then it has been accounted for large outbreaks of diarrhea globally. In Thailand, prior studies on pandemic and non-pandemic V. parahaemolyticus had mostly been done in the south. The incidence and molecular characterization of pandemic and non-pandemic strains in other parts of Thailand have not been fully characterized. This study examined the incidence of V. parahaemolyticus in seafood samples purchased in Bangkok and collected in eastern Thailand and characterized V. parahaemolyticus isolates. Potential virulence genes, VPaI-7, T3SS2, and biofilm were examined. Antimicrobial resistance (AMR) profiles and AMR genes (ARGs) were determined. Methods: V. parahaemolyticus was isolated from 190 marketed and farmed seafood samples by a culture method and confirmed by polymerase chain reaction (PCR). The incidence of pandemic and non-pandemic V. parahaemolyticus and VPaI-7, T3SS2, and biofilm genes was examined by PCR. AMR profiles were verified by a broth microdilution technique. The presence of ARGs was verified by genome analysis. V. parahaemolyticus characterization was done by multilocus sequence typing (MLST). A phylogenomic tree was built from nucleotide sequences by UBCG2.0 and RAxML softwares. Results: All 50 V. parahaemolyticus isolates including 21 pathogenic and 29 non-pathogenic strains from 190 samples had the toxRS/old sequence, indicating non-pandemic strains. All isolates had biofilm genes (VP0950, VP0952, and VP0962). None carried T3SS2 genes (VP1346 and VP1367), while VPaI-7 gene (VP1321) was seen in two isolates. Antimicrobial susceptibility profiles obtained from 36 V. parahaemolyticus isolates revealed high frequency of resistance to colistin (100%, 36/36) and ampicillin (83%, 30/36), but susceptibility to amoxicillin/clavulanic acid and piperacillin/tazobactam (100%, 36/36). Multidrug resistance (MDR) was seen in 11 isolates (31%, 11/36). Genome analysis revealed ARGs including blaCARB (100%, 36/36), tet(34) (83%, 30/36), tet(35) (42%, 15/36), qnrC (6%, 2/36), dfrA6 (3%, 1/36), and blaCTX-M-55 (3%, 1/36). Phylogenomic and MLST analyses classified 36 V. parahaemolyticus isolates into 5 clades, with 12 known and 13 novel sequence types (STs), suggesting high genetic variation among the isolates. Conclusions: Although none V. parahaemolyticus strains isolated from seafood samples purchased in Bangkok and collected in eastern Thailand were pandemic strains, around one third of isolates were MDR V. parahaemolyticus strains. The presence of resistance genes of the first-line antibiotics for V. parahaemolyticus infection raises a major concern for clinical treatment outcome since these resistance genes could be highly expressed under suitable circumstances.

Complete genome sequence of virulence-enhancing Siphophage VHS1 from Vibrio harveyi.

Vibrio harveyi siphophage 1 (VHS1) is a tailed phage with an icosahedral head of approximately 66 nm in diameter and an unornamented, flexible tail of approximately 153 nm in length. When Vibrio harveyi 1114GL is lysogenized with VHS1, its virulence for the black tiger shrimp (Penaeus monodon) increases by more than 100 times, and this coincides with production of a toxin(s) associated with shrimp hemocyte agglutination. Curiously, the lysogen does not show increased virulence for the whiteleg shrimp (Penaeus [Litopenaeus] vannamei). Here we present and annotate the complete, circular genome of VHS1 (81,509 kbp; GenBank accession number JF713456). By software analysis, the genome contains 125 putative open reading frames (ORFs), all of which appear to be located on the same DNA strand, similar to the case for many other bacteriophages. Most of the putative ORFs show no significant homology to known sequences in GenBank. Notable exceptions are ORFs for a putative DNA polymerase and putative phage structural proteins, including a portal protein, a phage tail tape measure protein, and a phage head protein. The last protein was identified as a component of the species-specific toxin mixture described above as being associated with agglutination of hemocytes from P. monodon.

Molecular characterization and expression profiling of transformer 2 and fruitless-like homologs in the black tiger shrimp, Penaeus monodon.

Transformer 2 (tra 2) and fruitless (fru) genes have been proven to play a key role in sex determination pathways in many Arthropods, including insects and crustaceans. In this study, a paralog of P. monodon tra 2 (Pmtra 2), P. monodon ovarian associated transformer 2 (PmOvtra 2) and 2 isoforms of P. monodon fruitless-like gene (Pmfru-1 and Pmfru-2) were identified and characterized. The full cDNA sequence of PmOvtra 2 consisted of 1,774 bp with the longest open reading frame (ORF) of 744 bp encoding for 247 amino acids. The PmOvtra 2 exhibited a predicted RNA-recognition motif (RRM) domain and two arginine-serine (RS) regions, suggesting its function in RNA splicing. The full cDNA sequence of Pmfru-1 consisted of 1,306 bp with 1,182 bp ORF encoding for 393 amino acids, whereas the full cDNA sequence of Pmfru-2 consisted of 1,858 bp with 1,437 bp ORF encoding 478 amino acids. The deduced amino acid sequences of Pmfru-1 and Pmfru-2 exhibited highly conserved domains of Fru proteins, including Broad-complex, Tramtrack and Bric-a-brac (BTB), and zinc finger (ZF) domains. In addition, Pmfru-1 and Pmfru-2 were suggestively originated from the same single genomic locus by genomic sequence analysis. Specifically, Pmfru pre-mRNA was alternatively spliced for Pmfru-1 and Pmfru-2 to include mutually exclusive exon 7 and exon 6, respectively. Temporal and spatial expression of PmOvtra 2, Pmfru-1, and Pmfru-2 were also investigated by qPCR. The results showed that all were expressed in early developmental stages with undifferentiated gonads starting from nauplius until postlarvae. The expression of PmOvtra 2 started at nauplius stage and gradually increased from mysis to postlarvae (PL) 1. However, the expression of Pmfru-1 was low at the nauplii stage and slightly increased from protozoea to PL5, whereas the expression of Pmfru-2 maintained a low level from nauplius to mysis and then gradually increased at the PL stages. Expressions of PmOvtra 2, Pmfru-1, and Pmfru-2 were detected in various tissues including nervous tissue, gill, heart, hepatopancreas, gut, and gonads. Interestingly, the sexually dimorphic expression of PmOvtra 2, Pmfru-1, and Pmfru-2 was demonstrated in fully developed gonads in which the ovary showed significantly higher expressions than the testis. The great difference in the expression pattern of PmOvtra 2, Pmfru-1, and Pmfru-2 in the ovary and testis suggested their roles in the female sex determination in P. monodon.

Applying Modified VP53A Recombinant Protein as an Anti-White Spot Syndrome Virus Biological Agent in Litopenaeus vannamei Farming.

Shrimp farming is an important economic activity. However, due to the spread of pathogens, shrimp aquaculture is becoming increasingly difficult. Many studies have confirmed that white spot syndrome virus (WSSV) recombinant proteins can inhibit viral infection. Among them, VP53 recombinant protein has been found to reduce mortality upon WSSV challenge. This study was conducted in Kaohsiung, Taiwan and reports the first field feeding trial to demonstrate that WSSV recombinant proteins can improve shrimp survival rates at a farming scale. Prior to the feeding trial, the shrimp were confirmed to be slightly infected with WSSV, Vibrio parahaemolyticus strains causing acute hepatopancreatic necrosis disease (AHPND), non-AHPND V. parahaemolyticus strains, and Enterocytozoon hepatopenaei (EHP), which are common pathogens that shrimp farmers often face. The shrimp were then divided into two groups: a control group (C group) fed with a commercial diet and a protein group (P group) fed with the same commercial feed with VP53 recombinant protein. Our findings indicated that the survival rate and expression of immune genes of the P group were higher than those of the C group. The intestinal microbiota of the two groups were also analysed. Collectively, our results confirmed that the recombinant WSSV envelope protein derivative can be used as an effective anti-virus biological agent in shrimp farms.

CRISPR-Cas fluorescent cleavage assay coupled with recombinase polymerase amplification for sensitive and specific detection of Enterocytozoon hepatopenaei.

Enterocytozoon hepatopenaei (EHP) is a parasite that infects pacific whiteleg shrimp, Penaeus vannamei, causing growth retardation and uneven size distributions that lead to severe losses in shrimp productivity. Routine monitoring is crucial to timely prevention and management of EHP, but field-deployable diagnostic kits for EHP are still scarce. Here, we proposed the use of recombinase polymerase amplification (RPA) and CRISPR-Cas12a fluorescence assay, henceforth RPA-Cas12a, for detection of EHP. Targeting ptp2 gene, RPA-Cas12a could detect as few as 50 copies of DNA and showed no reactivity with closely related microsporidia. The entire procedure could be performed at a temperature close to 37 °C within 1 h. Naked eye visualization was possible with UV/blue-light excitation or lateral flow detection. Thus, RPA-Cas12a is a rapid, sensitive and specific detection platform that requires no sophisticated equipment and shows promise for on-site surveillance of EHP.

Shewanella khirikhana sp. nov. - a shrimp pathogen isolated from a cultivation pond exhibiting early mortality syndrome.

Early mortality syndrome (EMS) in cultivated shrimp is of complex aetiology. One of the causes is acute hepatopancreatic necrosis disease (AHPND) caused by unique Vibrio isolates that carry two Pirvp toxin genes, but other causes of EMS remain mostly unexplained. Here, we describe the discovery of a Shewanella isolate TH2012T from an EMS/AHPND outbreak pond and demonstrate its virulence for shrimp (the mean lethal concentration of 105 colony-forming units per millilitre by immersion challenge) accompanied by distinctive histopathology, particularly of the ventral nerve cord and lymphoid organ but also including the digestive tract. On the basis of its complete genome sequence, multilocus phylogenetic trees, digital DNA-DNA hybridization analysis and differential phenotypic characteristics, we propose that Shewanella isolate TH2012T represents a novel species, separated sufficiently from the type strains S. litorisediminis and S. amazonensis to justify naming it Shewanella khirikhana sp. nov. Analysis of the TH2012T genome revealed no homologues of the Pirvp toxin genes but revealed a number of other potential virulence factors. It constitutes the first Shewanella isolate reported to be pathogenic to shrimp.

Complete Genome Sequence of Shewanella sp. Strain TH2012, Isolated from Shrimp in a Cultivation Pond Exhibiting Early Mortality Syndrome.

Here, we present the complete genome sequence of a Shewanella isolate, TH2012, from a shrimp pond in which shrimp exhibited early mortality syndrome (EMS)/acute hepatopancreatic necrosis disease (AHPND). The complete genome of TH2012 has a prophage-like element and a number of potential virulence factors, making TH2012 a possible contributing factor to EMS outbreaks.

Comparative genomics of Vibrio campbellii strains and core species of the Vibrio Harveyi clade.

The core of the Vibrio Harveyi clade contains V. harveyi, V. campbellii, V. owensii, V. jasicida, and V. rotiferianus. They are well recognized aquatic animal pathogens, but misclassification has been common due to similarities in their rDNA sequences and phenotypes. To better understand their evolutionary relationships and functional features, we sequenced a shrimp pathogen strain V. harveyi 1114GL, reclassified it as V. campbellii and compared this and 47 other sequenced Vibrio genomes in the Harveryi clade. A phylogeny based on 1,775 genes revealed that both V. owensii and V. jasicida were closer to V. campbellii than to V. harveyi and that V. campbellii strains can be divided into two distinct groups. Species-specific genes such as intimin and iron acquisition genes were identified in V. campbellii. In particular, the 1114GL strain contains two bacterial immunoglobulin-like genes for cell adhesion with 22 Big_2 domains that have been extensively reshuffled and are by far the most expanded among all species surveyed in this study. The 1114GL strain differed from ATCC BAA-1116 by ~9% at the synonymous sites, indicating high diversity within V. campbellii. Our study revealed the characteristics of V. campbellii in the Harveyi clade and the genetic basis for their wide-spread pathogenicity.

Intron size correlates positively with recombination rate in Caenorhabditis elegans.

A negative correlation between intron size and recombination rate has been reported for the Drosophila melanogaster and human genomes. Population-genetic models suggest that this pattern could be caused by an interaction between recombination rate and the efficacy of natural selection. To test this idea, we examined variation in intron size and recombination rate across the genome of the nematode Caenorhabditis elegans. Interestingly, we found that intron size correlated positively with recombination rate in this species.

A novel lectin domain-containing protein (LvCTLD) associated with response of the whiteleg shrimp Penaeus (Litopenaeus) vannamei to yellow head virus (YHV).

When using mRNA from gills of normal whiteleg shrimp Penaeus (Litopenaeus) vannamei as the tester and mRNA from yellow head virus (YHV)-infected shrimp as the driver, subtractive suppression hybridization (SSH) revealed that a novel EST clone of 198 bp with a putative C-type lectin-like domain (CTLD) was downregulated in YHV-infected shrimp. The clone nucleotide sequence had 99% identity with one contig MGID1052359 (1,380 bp) reported in an EST database of P. vannamei, and the presence of this target in normal shrimp was confirmed by RT-PCR using primers designed from the MGID1052359 sequence. Analysis of the primary structure of the deduced amino acid (a.a.) sequence of the contig revealed a short portion (40 a.a. residues) at its N-terminus with high similarity to a low density lipoprotein receptor (LDLR) class A domain and another 152 a.a. residues at its C-terminus with high similarity to a C-type lectin domain. Thus, the clone was named LvCTLD and three recombinant proteins (LvCTLD, the LDLR domain and the CTLD domain) were synthesized in a bacterial system based on its sequence. An in vitro encapsulation assay revealed that Sepharose 4B beads coated with rLvCTLD were encapsulated by shrimp hemocytes and that melanization followed by 24 h post-encapsulation. The encapsulation activity of rLvCTLD was inhibited by 100 mM galactose, but not mannose or EDTA. In vivo injection of rLvCTLD or rLvCTLD plus YHV resulted in a significant elevation of PO activity in the hemolymph of the challenged shrimp when compared to shrimp injected with buffer, suggesting that rLvCTLD could activate the proPO system. An ELISA test revealed that rLvCTLD could bind to YHV particles in the presence of shrimp hemolymph. Phylogenetic analysis suggested that the LvCTLD sequence was more closely related to an antiviral gene found in Penaeus monodon (PmAV) than to other reported shrimp lectins. Taken together, we conclude that a novel shrimp LvCTLD is a host recognition molecule involved in the shrimp defense mechanism against YHV via recruitment of hemocytes, probably at the site of viral infection, and via activation of the proPO system.

Gene number expansion and contraction in vertebrate genomes with respect to invertebrate genomes.

Where did vertebrate genes come from? Here we address this question by analyzing eight completely sequenced land vertebrate genomes and six completely sequenced invertebrate genomes. Approximately 70% of the vertebrate genes can be found in the six invertebrate genomes with the standard homology search criteria (denoted as V.MCL), another approximately 6% can be found with relaxed search criteria, and an additional approximately 2% can be found in sequenced fungal and bacterial genomes. Thus, a substantial proportion of vertebrate genes (approximately 22%) cannot be found in the nonvertebrate genomes studied (denoted as Vonly). Interestingly, genes in Vonly are predominantly singletons, while the majority of genes in the other three groups belong to gene families. The proteins of Vonly tend to evolve faster than those of V.MCL. Surprisingly, in many cases the family sizes in V.MCL are only as large as or even smaller than their counterparts in the invertebrates, contrary to the general perception of a larger family size in vertebrates. Interestingly, in comparison with the family size in invertebrates, vertebrate gene families involved in regulation, signal transduction, transcription, protein transport, and protein modification tend to be expanded, whereas those involved in metabolic processes tend to be contracted. Furthermore, for almost all of the functional categories with family size expansion in vertebrates, the number of gene types (i.e., the number of singletons plus the number of gene families) tends to be over-represented in Vonly, but under-represented in V.MCL. Our study suggests that gene function is a major determinant of gene family size.

Protein function, connectivity, and duplicability in yeast.

Protein-protein interaction networks have evolved mainly through connectivity rewiring and gene duplication. However, how protein function influences these processes and how a network grows in time have not been well studied. Using protein-protein interaction data and genomic data from the budding yeast, we first examined whether there is a correlation between the age and connectivity of yeast proteins. A steady increase in connectivity with protein age is observed for yeast proteins except for those that can be traced back to Eubacteria. Second, we investigated whether protein connectivity and duplicability vary with gene function. We found a higher average duplicability for proteins interacting with external environments than for proteins localized within intracellular compartments. For example, proteins that function in the cell periphery (mainly transporters) show a high duplicability but are lowly connected. Conversely, proteins that function within the nucleus (e.g., transcription, RNA and DNA metabolisms, and ribosome biogenesis and assembly) are highly connected but have a low duplicability. Finally, we found a negative correlation between protein connectivity and duplicability.

Higher gene duplicabilities for metabolic proteins than for nonmetabolic proteins in yeast and E. coli.

Although the evolutionary significance of gene duplication has long been appreciated, it remains unclear what factors determine gene duplicability. In this study we investigated whether metabolism is an important determinant of gene duplicability because cellular metabolism is crucial for the survival and reproduction of an organism. Using genomic data and metabolic pathway data from the yeast (Saccharomyces cerevisiae) and Escherichia coli, we found that metabolic proteins indeed tend to have higher gene duplicability than nonmetabolic proteins. Moreover, a detailed analysis of metabolic pathways in these two organisms revealed that genes in the central metabolic pathways and the catabolic pathways have, on average, higher gene duplicability than do other genes and that most genes in anabolic pathways are single-copy genes.