A score system for quality evaluation of RNA sequence tags: an improvement for gene expression profiling

Background: High-throughput molecular approaches for gene expression profiling, such as Serial Analysis of Gene Expression (SAGE), Massively Parallel Signature Sequencing (MPSS) or Sequencing-by-Synthesis (SBS) represent powerful techniques that provide global transcription profiles of different cell types through sequencing of short fragments of transcripts, denominated sequence tags. These techniques have improved our understanding about the relationships between these expression profiles and cellular phenotypes. Despite this, more reliable datasets are still necessary. In this work, we present a web-based tool named S3T: Score System for Sequence Tags, to index sequenced tags in accordance with their reliability. This is made through a series of evaluations based on a defined rule set. S3T allows the identification/selection of tags, considered more reliable for further gene expression analysis. Results: This methodology was applied to a public SAGE dataset. In order to compare data before and after filtering, a hierarchical clustering analysis was performed in samples from the same type of tissue, in distinct biological conditions, using these two datasets. Our results provide evidences suggesting that it is possible to find more congruous clusters after using S3T scoring system. Conclusion: These results substantiate the proposed application to generate more reliable data. This is a significant contribution for determination of global gene expression profiles. The library analysis with S3T is freely available at http://gdm.fmrp.usp.br/s3t/.S3T source code and datasets can also be downloaded from the aforementioned website.

Chagas Cardiomyopathy Manifestations and Trypanosoma cruzi Genotypes Circulating in Chronic Chagasic Patients

Chagas disease caused by Trypanosoma cruzi is a complex disease that is endemic and an important problem in public health in Latin America. The T. cruzi parasite is classified into six discrete taxonomic units (DTUs) based on the recently proposed nomenclature (TcI, TcII, TcIII, TcIV, TcV and TcVI). The discovery of genetic variability within TcI showed the presence of five genotypes (Ia, Ib, Ic, Id and Ie) related to the transmission cycle of Chagas disease. In Colombia, TcI is more prevalent but TcII has also been reported, as has mixed infection by both TcI and TcII in the same Chagasic patient. The objectives of this study were to determine the T. cruzi DTUs that are circulating in Colombian chronic Chagasic patients and to obtain more information about the molecular epidemiology of Chagas disease in Colombia. We also assessed the presence of electrocardiographic, radiologic and echocardiographic abnormalities with the purpose of correlating T. cruzi genetic variability and cardiac disease. Molecular characterization was performed in Colombian adult chronic Chagasic patients based on the intergenic region of the mini-exon gene, the 24S alpha and 18S regions of rDNA and the variable region of satellite DNA, whereby the presence of T. cruzi I, II, III and IV was detected. In our population, mixed infections also occurred, with TcI-TcII, TcI-TcIII and TcI-TcIV, as well as the existence of the TcI genotypes showing the presence of genotypes Ia and Id. Patients infected with TcI demonstrated a higher prevalence of cardiac alterations than those infected with TcII. These results corroborate the predominance of TcI in Colombia and show the first report of TcIII and TcIV in Colombian Chagasic patients. Findings also indicate that Chagas cardiomyopathy manifestations are more correlated with TcI than with TcII in Colombia.

Selection of suitable housekeeping genes for expression analysis in glioblastoma using quantitative RT-PCR

Background: Considering the broad variation in the expression of housekeeping genes among tissues and experimental situations, studies using quantitative RT-PCR require strict definition of adequate endogenous controls. For glioblastoma, the most common type of tumor in the central nervous system, there was no previous report regarding this issue. Results: Here we show that amongst seven frequently used housekeeping genes TBP and HPRT1 are adequate references for glioblastoma gene expression analysis. Evaluation of the expression levels of 12 target genes utilizing different endogenous controls revealed that the normalization method applied might introduce errors in the estimation of relative quantities. Genes presenting expression levels which do not significantly differ between tumor and normal tissues can be considered either increased or decreased if unsuitable reference genes are applied. Most importantly, genes showing significant differences in expression levels between tumor and normal tissues can be missed. We also demonstrated that the Holliday Junction Recognizing Protein, a novel DNA repair protein over expressed in lung cancer, is extremely over-expressed in glioblastoma, with a median change of about 134 fold. Conclusion: Altogether, our data show the relevance of previous validation of candidate control genes for each experimental model and indicate TBP plus HPRT1 as suitable references for studies on glioblastoma gene expression.

RNA interference-mediated knockdown of CD49e (alpha 5 integrin chain) in human thymic epithelial cells modulates the expression of multiple genes and decreases thymocyte adhesion

Background: The thymus is a central lymphoid organ, in which bone marrow-derived T cell precursors undergo a complex process of maturation. Developing thymocytes interact with thymic microenvironment in a defined spatial order. A component of thymic microenvironment, the thymic epithelial cells, is crucial for the maturation of T-lymphocytes through cell-cell contact, cell matrix interactions and secretory of cytokines/chemokines. There is evidence that extracellular matrix molecules play a fundamental role in guiding differentiating thymocytes in both cortical and medullary regions of the thymic lobules. The interaction between the integrin alpha 5 beta 1 (CD49e/CD29; VLA-5) and fibronectin is relevant for thymocyte adhesion and migration within the thymic tissue. Our previous results have shown that adhesion of thymocytes to cultured TEC line is enhanced in the presence of fibronectin, and can be blocked with anti-VLA-5 antibody. Results: Herein, we studied the role of CD49e expressed by the human thymic epithelium. For this purpose we knocked down the CD49e by means of RNA interference. This procedure resulted in the modulation of more than 100 genes, some of them coding for other proteins also involved in adhesion of thymocytes; others related to signaling pathways triggered after integrin activation, or even involved in the control of F-actin stress fiber formation. Functionally, we demonstrated that disruption of VLA-5 in human TEC by CD49e-siRNA-induced gene knockdown decreased the ability of TEC to promote thymocyte adhesion. Such a decrease comprised all CD4/CD8-defined thymocyte subsets. Conclusion: Conceptually, our findings unravel the complexity of gene regulation, as regards key genes involved in the heterocellular cell adhesion between developing thymocytes and the major component of the thymic microenvironment, an interaction that is a mandatory event for proper intrathymic T cell differentiation.

Group I introns and associated homing endonuclease genes reveals a clinal structure for Porphyra spiralis var. amplifolia (Bangiales, Rhodophyta) along the Eastern coast of South America

Background: Group I introns are found in the nuclear small subunit ribosomal RNA gene (SSU rDNA) of some species of the genus Porphyra (Bangiales, Rhodophyta). Size polymorphisms in group I introns has been interpreted as the result of the degeneration of homing endonuclease genes (HEG) inserted in peripheral loops of intron paired elements. In this study, intron size polymorphisms were characterized for different Porphyra spiralis var. amplifolia (PSA) populations on the Southern Brazilian coast, and were used to infer genetic relationships and genetic structure of these PSA populations, in addition to cox2-3 and rbcL-S regions. Introns of different sizes were tested qualitatively for in vitro self-splicing. Results: Five intron size polymorphisms within 17 haplotypes were obtained from 80 individuals representing eight localities along the distribution of PSA in the Eastern coast of South America. In order to infer genetic structure and genetic relationships of PSA, these polymorphisms and haplotypes were used as markers for pairwise Fst analyses, Mantel's test and median joining network. The five cox2-3 haplotypes and the unique rbcL-S haplotype were used as markers for summary statistics, neutrality tests Tajima's D and Fu's Fs and for median joining network analyses. An event of demographic expansion from a population with low effective number, followed by a pattern of isolation by distance was obtained for PSA populations with the three analyses. In vitro experiments have shown that introns of different lengths were able to self-splice from pre-RNA transcripts. Conclusion: The findings indicated that degenerated HEGs are reminiscent of the presence of a full-length and functional HEG, once fixed for PSA populations. The cline of HEG degeneration determined the pattern of isolation by distance. Analyses with the other markers indicated an event of demographic expansion from a population with low effective number. The different degrees of degeneration of the HEG do not refrain intron self-splicing. To our knowledge, this was the first study to address intraspecific evolutionary history of a nuclear group I intron; to use nuclear, mitochondrial and chloroplast DNA for population level analyses of Porphyra; and intron size polymorphism as a marker for population genetics.

The taxonomic status of the endangered thin-spined porcupine, Chaetomys subspinosus (Olfers, 1818), based on molecular and karyologic data

Background: The thin-spined porcupine, also known as the bristle-spined rat, Chaetomys subspinosus (Olfers, 1818), the only member of its genus, figures among Brazilian endangered species. In addition to being threatened, it is poorly known, and even its taxonomic status at the family level has long been controversial. The genus Chaetomys was originally regarded as a porcupine in the family Erethizontidae, but some authors classified it as a spiny-rat in the family Echimyidae. Although the dispute seems to be settled in favor of the erethizontid advocates, further discussion of its affinities should be based on a phylogenetic framework. In the present study, we used nucleotide-sequence data from the complete mitochondrial cytochrome b gene and karyotypic information to address this issue. Our molecular analyses included one individual of Chaetomys subspinosus from the state of Bahia in northeastern Brazil, and other hystricognaths. Results: All topologies recovered in our molecular phylogenetic analyses strongly supported Chaetomys subspinosus as a sister clade of the erethizontids. Cytogenetically, Chaetomys subspinosus showed 2n = 52 and FN = 76. Although the sexual pair could not be identified, we assumed that the X chromosome is biarmed. The karyotype included 13 large to medium metacentric and submetacentric chromosome pairs, one small subtelocentric pair, and 12 small acrocentric pairs. The subtelocentric pair 14 had a terminal secondary constriction in the short arm, corresponding to the nucleolar organizer region (Ag-NOR), similar to the erethizontid Sphiggurus villosus, 2n = 42 and FN = 76, and different from the echimyids, in which the secondary constriction is interstitial. Conclusion: Both molecular phylogenies and karyotypical evidence indicated that Chaetomys is closely related to the Erethizontidae rather than to the Echimyidae, although in a basal position relative to the rest of the Erethizontidae. The high levels of molecular and morphological divergence suggest that Chaetomys belongs to an early radiation of the Erethizontidae that may have occurred in the Early Miocene, and should be assigned to its own subfamily, the Chaetomyinae.

Phylogenetic position of the acariform mites: sensitivity to homology assessment under total evidence

Background: Mites (Acari) have traditionally been treated as monophyletic, albeit composed of two major lineages: Acariformes and Parasitiformes. Yet recent studies based on morphology, molecular data, or combinations thereof, have increasingly drawn their monophyly into question. Furthermore, the usually basal (molecular) position of one or both mite lineages among the chelicerates is in conflict to their morphology, and to the widely accepted view that mites are close relatives of Ricinulei. Results: The phylogenetic position of the acariform mites is examined through employing SSU, partial LSU sequences, and morphology from 91 chelicerate extant terminals (forty Acariformes). In a static homology framework, molecular sequences were aligned using their secondary structure as guide, whereby regions of ambiguous alignment were discarded, and pre-aligned sequences analyzed under parsimony and different mixed models in a Bayesian inference. Parsimony and Bayesian analyses led to trees largely congruent concerning infraordinal, well-supported branches, but with low support for inter-ordinal relationships. An exception is Solifugae + Acariformes (P. P = 100%, J. = 0.91). In a dynamic homology framework, two analyses were run: a standard POY analysis and an analysis constrained by secondary structure. Both analyses led to largely congruent trees; supporting a (Palpigradi (Solifugae Acariformes)) clade and Ricinulei as sister group of Tetrapulmonata with the topology (Ricinulei (Amblypygi (Uropygi Araneae))). Combined analysis with two different morphological data matrices were run in order to evaluate the impact of constraining the analysis on the recovered topology when employing secondary structure as a guide for homology establishment. The constrained combined analysis yielded two topologies similar to the exclusively molecular analysis for both morphological matrices, except for the recovery of Pedipalpi instead of the (Uropygi Araneae) clade. The standard (direct optimization) POY analysis, however, led to the recovery of trees differing in the absence of the otherwise well-supported group Solifugae + Acariformes. Conclusions: Previous studies combining ribosomal sequences and morphology often recovered topologies similar to purely morphological analyses of Chelicerata. The apparent stability of certain clades not recovered here, like Haplocnemata and Acari, is regarded as a byproduct of the way the molecular homology was previously established using the instrumentalist approach implemented in POY. Constraining the analysis by a priori homology assessment is defended here as a way of maintaining the severity of the test when adding new data to the analysis. Although the strength of the method advocated here is keeping phylogenetic information from regions usually discarded in an exclusively static homology framework; it still has the inconvenience of being uninformative on the effect of alignment ambiguity on resampling methods of clade support estimation. Finally, putative morphological apomorphies of Solifugae + Acariformes are the reduction of the proximal cheliceral podomere, medial abutting of the leg coxae, loss of sperm nuclear membrane, and presence of differentiated germinative and secretory regions in the testis delivering their products into a common lumen.

Complete Genome Sequence of Mycoplasma suis and Insights into Its Biology and Adaption to an Erythrocyte Niche

Mycoplasma suis, the causative agent of porcine infectious anemia, has never been cultured in vitro and mechanisms by which it causes disease are poorly understood. Thus, the objective herein was to use whole genome sequencing and analysis of M. suis to define pathogenicity mechanisms and biochemical pathways. M. suis was harvested from the blood of an experimentally infected pig. Following DNA extraction and construction of a paired end library, whole-genome sequencing was performed using GS-FLX (454) and Titanium chemistry. Reads on paired-end constructs were assembled using GS De Novo Assembler and gaps closed by primer walking; assembly was validated by PFGE. Glimmer and Manatee Annotation Engine were used to predict and annotate protein-coding sequences (CDS). The M. suis genome consists of a single, 742,431 bp chromosome with low G+C content of 31.1%. A total of 844 CDS, 3 single copies, unlinked rRNA genes and 32 tRNAs were identified. Gene homologies and GC skew graph show that M. suis has a typical Mollicutes oriC. The predicted metabolic pathway is concise, showing evidence of adaptation to blood environment. M. suis is a glycolytic species, obtaining energy through sugars fermentation and ATP-synthase. The pentose-phosphate pathway, metabolism of cofactors and vitamins, pyruvate dehydrogenase and NAD(+) kinase are missing. Thus, ribose, NADH, NADPH and coenzyme A are possibly essential for its growth. M. suis can generate purines from hypoxanthine, which is secreted by RBCs, and cytidine nucleotides from uracil. Toxins orthologs were not identified. We suggest that M. suis may cause disease by scavenging and competing for host nutrients, leading to decreased life-span of RBCs. In summary, genome analysis shows that M. suis is dependent on host cell metabolism and this characteristic is likely to be linked to its pathogenicity. The prediction of essential nutrients will aid the development of in vitro cultivation systems.

Identification of protein-coding and non-coding RNA expression profiles in CD34(+) and in stromal cells in refractory anemia with ringed sideroblasts

Background: Myelodysplastic syndromes (MDS) are a group of clonal hematological disorders characterized by ineffective hematopoiesis with morphological evidence of marrow cell dysplasia resulting in peripheral blood cytopenia. Microarray technology has permitted a refined high-throughput mapping of the transcriptional activity in the human genome. Non-coding RNAs (ncRNAs) transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression, and in the regulation of exon-skipping and intron retention. Characterization of ncRNAs in progenitor cells and stromal cells of MDS patients could be strategic for understanding gene expression regulation in this disease. Methods: In this study, gene expression profiles of CD34(+) cells of 4 patients with MDS of refractory anemia with ringed sideroblasts (RARS) subgroup and stromal cells of 3 patients with MDS-RARS were compared with healthy individuals using 44 k combined intron-exon oligoarrays, which included probes for exons of protein-coding genes, and for non-coding RNAs transcribed from intronic regions in either the sense or antisense strands. Real-time RT-PCR was performed to confirm the expression levels of selected transcripts. Results: In CD34(+) cells of MDS-RARS patients, 216 genes were significantly differentially expressed (q-value <= 0.01) in comparison to healthy individuals, of which 65 (30%) were non-coding transcripts. In stromal cells of MDS-RARS, 12 genes were significantly differentially expressed (q-value <= 0.05) in comparison to healthy individuals, of which 3 (25%) were non-coding transcripts. Conclusions: These results demonstrated, for the first time, the differential ncRNA expression profile between MDS-RARS and healthy individuals, in CD34(+) cells and stromal cells, suggesting that ncRNAs may play an important role during the development of myelodysplastic syndromes.

Identification of archaeal proteins that affect the exosome function in vitro

Background: The archaeal exosome is formed by a hexameric RNase PH ring and three RNA binding subunits and has been shown to bind and degrade RNA in vitro. Despite extensive studies on the eukaryotic exosome and on the proteins interacting with this complex, little information is yet available on the identification and function of archaeal exosome regulatory factors. Results: Here, we show that the proteins PaSBDS and PaNip7, which bind preferentially to poly-A and AU-rich RNAs, respectively, affect the Pyrococcus abyssi exosome activity in vitro. PaSBDS inhibits slightly degradation of a poly-rA substrate, while PaNip7 strongly inhibits the degradation of poly-A and poly-AU by the exosome. The exosome inhibition by PaNip7 appears to depend at least partially on its interaction with RNA, since mutants of PaNip7 that no longer bind RNA, inhibit the exosome less strongly. We also show that FITC-labeled PaNip7 associates with the exosome in the absence of substrate RNA. Conclusions: Given the high structural homology between the archaeal and eukaryotic proteins, the effect of archaeal Nip7 and SBDS on the exosome provides a model for an evolutionarily conserved exosome control mechanism.

Haplotype distribution of five nuclear genes based on network genealogies and Bayesian inference indicates that Trypanosoma cruzi hybrid strains are polyphyletic

Chagas disease is still a major public health problem in Latin America. Its causative agent, Trypanosoma cruzi, can be typed into three major groups, T. cruzi I, T. cruzi II and hybrids. These groups each have specific genetic characteristics and epidemiological distributions. Several highly virulent strains are found in the hybrid group; their origin is still a matter of debate. The null hypothesis is that the hybrids are of polyphyletic origin, evolving independently from various hybridization events. The alternative hypothesis is that all extant hybrid strains originated from a single hybridization event. We sequenced both alleles of genes encoding EF-1 alpha, actin and SSU rDNA of 26 T. cruzi strains and DHFR-TS and TR of 12 strains. This information was used for network genealogy analysis and Bayesian phylogenies. We found T. cruzi I and T. cruzi II to be monophyletic and that all hybrids had different combinations of T. cruzi I and T. cruzi II haplotypes plus hybrid-specific haplotypes. Bootstrap values (networks) and posterior probabilities (Bayesian phylogenies) of clades supporting the monophyly of hybrids were far below the 95% confidence interval, indicating that the hybrid group is polyphyletic. We hypothesize that T. cruzi I and T. cruzi II are two different species and that the hybrids are extant representatives of independent events of genome hybridization, which sporadically have sufficient fitness to impact on the epidemiology of Chagas disease.

Molecular and morphological characterization of Amathia distans Busk and Amathia brasiliensis Busk (Bryozoa: Ctenostomata) from the tropical and subtropical Western Atlantic

Morphological and molecular analyses have proven to be complementary tools of taxonomic information for the redescription of the ctenostome bryozoans Amathia brasiliensis Busk, 1886 and Amathia distans Busk, 1886. The two species, originally described from material collected by the `Challenger` expedition but synonymized by later authors, now have their status fixed by means of the selection of lectotypes, morphological observations and analyses of DNA sequences described here. The morphological characters allowing the identification of living and/or preserved specimens are (1) A. brasiliensis: whitish-pale pigment spots in the frontal surface of stolons and zooids, and a wide stolon with biserial zooid clusters growing in clockwise and anti-clockwise spirals along it, the spirality direction being maintained from maternal to daughter stolons; and (2) A. distans: bright yellow pigment spots in stolonal and zooidal surfaces including lophophores, and a slender stolon, thickly cuticularized, with biserial zooid clusters growing in clockwise and anti-clockwise spirals along it and the spirality direction not maintained from maternal to daughter stolons. Pairwise comparisons of DNA sequences of the mitochondrial genes cytochrome c oxidase subunit I and large ribosomal RNA subunit revealed deep genetic divergence between A. brasiliensis and A. distans. Finally, analyses of those sequences within a Bayesian phylogenetic context recovered their genealogical species status.

Diversity of free-living marine nematodes (Enoplida) from Baja California assessed by integrative taxonomy

We used morphological and molecular approaches to evaluate the diversity of free-living marine nematodes (order Enoplida) at four coastal sites in the Gulf of California and three on the Pacific coast of Baja California, Mexico. We identified 22 morphological species belonging to six families, of which Thoracostomopsidae and Oncholaimidae were the most diverse. The genus Mesacanthion (Thoracostomopsidae) was the most widespread and diverse. Five allopatric species, genetically and morphologically differentiated, were found in two localities in the Gulf of California (M. sp1 and M. sp2) and three in the Pacific coast (M. sp3, M. sp4 and M. sp5). Overall, we produced 19 and 20 sequences for the 18S and 28S genes, respectively. Neither gene displayed intraspecific polymorphisms, which allowed us to establish that some morphological variation was likely either ontogenetic or due to phenotypic plasticity. Although 18S and 28S phylogenies were topologically congruent (incongruence length difference test, P > 0.05), divergences between species were much higher in the 28S gene. Moreover, this gene possessed a stronger phylogenetic signal to resolve relationships involving Rhabdodemania and Bathylaimus. On the other hand, the close relationship of Pareurystomina (Enchilidiidae) with oncholaimids warrants further study. The 28S sequences (D2D3 domain) may be better suited for DNA barcoding of marine nematodes than those from the 18S rDNA, particularly for differentiating closely related or cryptic species. Finally, our results underline the relevance of adopting an integrative approach encompassing morphological and molecular analyses to improve the assessment of marine nematode diversity and advance their taxonomy.

Highly Toxic Microcystis aeruginosa Strain, Isolated from So Paulo-Brazil, Produce Hepatotoxins and Paralytic Shellfish Poison Neurotoxins

While evaluating several laboratory-cultured cyanobacteria strains for the presence of paralytic shellfish poison neurotoxins, the hydrophilic extract of Microcystis aeruginosa strain SPC777-isolated from Billings`s reservoir, So Paulo, Brazil-was found to exhibit lethal neurotoxic effect in mouse bioassay. The in vivo test showed symptoms that unambiguously were those produced by PSP. In order to identify the presence of neurotoxins, cells were lyophilized, and the extracts were analyzed by HPLC-FLD and HPLC-MS. HPLC-FLD analysis revealed four main Gonyautoxins: GTX4(47.6%), GTX2(29.5%), GTX1(21.9%), and GTX3(1.0%). HPLC-MS analysis, on other hand, confirmed both epimers, with positive Zwitterions M(+) 395.9 m/z for GTX3/GTX2 and M(+) 411 m/z for GTX4/GTX1 epimers. The hepatotoxins (Microcystins) were also evaluated by ELISA and HPLC-MS analyses. Positive immunoreaction was observed by ELISA assay. Alongside, the HPLC-MS analyses revealed the presence of [l-ser(7)] MCYST-RR. The N-methyltransferase (NMT) domain of the microcystin synthetase gene mcyA was chosen as the target sequence to detect the presence of the mcy gene cluster. PCR amplification of the NMT domain, using the genomic DNA of the SPC777 strain and the MSF/MSR primer set, resulted in the expected 1,369 bp product. The phylogenetic analyses grouped the NMT sequence with the NMT sequences of other known Microcystis with high bootstrap support. The taxonomical position of M. aeruginosa SPC777 was confirmed by a detailed morphological description and a phylogenetic analysis of 16S rRNA gene sequence. Therefore, co-production of PSP neurotoxins and microcystins by an isolated M. aeruginosa strain is hereby reported for the first time.

Diversity and identification of methanogenic archaea and sulphate-reducing bacteria in sediments from a pristine tropical mangrove

Mangrove sediments are anaerobic ecosystems rich in organic matter. This environment is optimal for anaerobic microorganisms, such as sulphate-reducing bacteria and methanogenic archaea, which are responsible for nutrient cycling. In this study, the diversity of these two functional guilds was evaluated in a pristine mangrove forest using denaturing gradient gel electrophoresis (DGGE) and clone library sequencing in a 50 cm vertical profile sampled every 5.0 cm. DGGE profiles indicated that both groups presented higher richness in shallow samples (0-30 cm) with a steep decrease in richness beyond that depth. According to redundancy analysis, this alteration significantly correlated with a decrease in the amount of organic matter. Clone library sequencing indicated that depth had a strong effect on the selection of dissimilatory sulphate reductase (dsrB) operational taxonomic units (OTUs), as indicated by the small number of shared OTUs found in shallow (0.0 cm) and deep (40.0 cm) libraries. On the other hand, methyl coenzyme-M reductase (mcrA) libraries indicated that most of the OTUs found in the shallow library were present in the deep library. These results show that these two guilds co-exist in these mangrove sediments and indicate important roles for these organisms in nutrient cycling within this ecosystem.