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.

Genome Size, Karyotype Polymorphism and Chromosomal Evolution in Trypanosoma cruzi

Background: The Trypanosoma cruzi genome was sequenced from a hybrid strain (CL Brener). However, high allelic variation and the repetitive nature of the genome have prevented the complete linear sequence of chromosomes being determined. Determining the full complement of chromosomes and establishing syntenic groups will be important in defining the structure of T. cruzi chromosomes. A large amount of information is now available for T. cruzi and Trypanosoma brucei, providing the opportunity to compare and describe the overall patterns of chromosomal evolution in these parasites. Methodology/Principal Findings: The genome sizes, repetitive DNA contents, and the numbers and sizes of chromosomes of nine strains of T. cruzi from four lineages (TcI, TcII, TcV and TcVI) were determined. The genome of the TcI group was statistically smaller than other lineages, with the exception of the TcI isolate Tc1161 (Jose-IMT). Satellite DNA content was correlated with genome size for all isolates, but this was not accompanied by simultaneous amplification of retrotransposons. Regardless of chromosomal polymorphism, large syntenic groups are conserved among T. cruzi lineages. Duplicated chromosome-sized regions were identified and could be retained as paralogous loci, increasing the dosage of several genes. By comparing T. cruzi and T. brucei chromosomes, homologous chromosomal regions in T. brucei were identified. Chromosomes Tb9 and Tb11 of T. brucei share regions of syntenic homology with three and six T. cruzi chromosomal bands, respectively. Conclusions: Despite genome size variation and karyotype polymorphism, T. cruzi lineages exhibit conservation of chromosome structure. Several syntenic groups are conserved among all isolates analyzed in this study. The syntenic regions are larger than expected if rearrangements occur randomly, suggesting that they are conserved owing to positive selection. Mapping of the syntenic regions on T. cruzi chromosomal bands provides evidence for the occurrence of fusion and split events involving T. brucei and T. cruzi chromosomes.

Genome-wide SNP genotyping highlights the role of natural selection in Plasmodium falciparum population divergence

Background: The malaria parasite Plasmodium falciparum exhibits abundant genetic diversity, and this diversity is key to its success as a pathogen. Previous efforts to study genetic diversity in P. falciparum have begun to elucidate the demographic history of the species, as well as patterns of population structure and patterns of linkage disequilibrium within its genome. Such studies will be greatly enhanced by new genomic tools and recent large-scale efforts to map genomic variation. To that end, we have developed a high throughput single nucleotide polymorphism (SNP) genotyping platform for P. falciparum. Results: Using an Affymetrix 3,000 SNP assay array, we found roughly half the assays (1,638) yielded high quality, 100% accurate genotyping calls for both major and minor SNP alleles. Genotype data from 76 global isolates confirm significant genetic differentiation among continental populations and varying levels of SNP diversity and linkage disequilibrium according to geographic location and local epidemiological factors. We further discovered that nonsynonymous and silent (synonymous or noncoding) SNPs differ with respect to within-population diversity, interpopulation differentiation, and the degree to which allele frequencies are correlated between populations. Conclusions: The distinct population profile of nonsynonymous variants indicates that natural selection has a significant influence on genomic diversity in P. falciparum, and that many of these changes may reflect functional variants deserving of follow-up study. Our analysis demonstrates the potential for new high-throughput genotyping technologies to enhance studies of population structure, natural selection, and ultimately enable genome-wide association studies in P. falciparum to find genes underlying key phenotypic traits.

Leishmania-Specific Surface Antigens Show Sub-Genus Sequence Variation and Immune Recognition

Background: A family of hydrophilic acylated surface (HASP) proteins, containing extensive and variant amino acid repeats, is expressed at the plasma membrane in infective extracellular (metacyclic) and intracellular (amastigote) stages of Old World Leishmania species. While HASPs are antigenic in the host and can induce protective immune responses, the biological functions of these Leishmania-specific proteins remain unresolved. Previous genome analysis has suggested that parasites of the sub-genus Leishmania (Viannia) have lost HASP genes from their genomes. Methods/Principal Findings: We have used molecular and cellular methods to analyse HASP expression in New World Leishmania mexicana complex species and show that, unlike in L. major, these proteins are expressed predominantly following differentiation into amastigotes within macrophages. Further genome analysis has revealed that the L. (Viannia) species, L. (V.) braziliensis, does express HASP-like proteins of low amino acid similarity but with similar biochemical characteristics, from genes present on a region of chromosome 23 that is syntenic with the HASP/SHERP locus in Old World Leishmania species and the L. (L.) mexicana complex. A related gene is also present in Leptomonas seymouri and this may represent the ancestral copy of these Leishmania-genus specific sequences. The L. braziliensis HASP-like proteins (named the orthologous (o) HASPs) are predominantly expressed on the plasma membrane in amastigotes and are recognised by immune sera taken from 4 out of 6 leishmaniasis patients tested in an endemic region of Brazil. Analysis of the repetitive domains of the oHASPs has shown considerable genetic variation in parasite isolates taken from the same patients, suggesting that antigenic change may play a role in immune recognition of this protein family. Conclusions/Significance: These findings confirm that antigenic hydrophilic acylated proteins are expressed from genes in the same chromosomal region in species across the genus Leishmania. These proteins are surface-exposed on amastigotes (although L. (L.) major parasites also express HASPB on the metacyclic plasma membrane). The central repetitive domains of the HASPs are highly variant in their amino acid sequences, both within and between species, consistent with a role in immune recognition in the host.

Three-dimensional visualization of human hemoglobin phenotypes with HPLC

Hemoglobinopathies were included in the Brazilian Neonatal Screening Program on June 6, 2001. Automated high-performance liquid chromatography (HPLC) was indicated as one of the diagnostic methods. The amount of information generated by these systems is immense, and the behavior of groups cannot always be observed in individual analyses. Three-dimensional (3-D) visualization techniques can be applied to extract this information, for extracting patterns, trends or relations from the results stored in databases. We applied the 3-D visualization tool to analyze patterns in the results of hemoglobinopathy based on neonatal diagnosis by HPLC. The laboratory results of 2520 newborn analyses carried out in 2001 and 2002 were used. The ""Fast"", ""F1"", ""F"" and ""A"" peaks, which were detected by the analytical system, were chosen as attributes for mapping. To establish a behavior pattern, the results were classified into groups according to hemoglobin phenotype: normal (N = 2169), variant (N = 73) and thalassemia (N = 279). 3-D visualization was made with the FastMap DB tool; there were two distribution patterns in the normal group, due to variation in the amplitude of the values obtained by HPLC for the F1 window. It allowed separation of the samples with normal Hb from those with alpha thalassemia, based on a significant difference (P < 0.05) between the mean values of the ""Fast"" and ""A"" peaks, demonstrating the need for better evaluation of chromatograms; this method could be used to help diagnose alpha thalassemia in newborns.

Gene Expression Noise in Spatial Patterning: hunchback Promoter Structure Affects Noise Amplitude and Distribution in Drosophila Segmentation

Positional information in developing embryos is specified by spatial gradients of transcriptional regulators. One of the classic systems for studying this is the activation of the hunchback (hb) gene in early fruit fly (Drosophila) segmentation by the maternally-derived gradient of the Bicoid (Bcd) protein. Gene regulation is subject to intrinsic noise which can produce variable expression. This variability must be constrained in the highly reproducible and coordinated events of development. We identify means by which noise is controlled during gene expression by characterizing the dependence of hb mRNA and protein output noise on hb promoter structure and transcriptional dynamics. We use a stochastic model of the hb promoter in which the number and strength of Bcd and Hb (self-regulatory) binding sites can be varied. Model parameters are fit to data from WT embryos, the self-regulation mutant hb(14F), and lacZ reporter constructs using different portions of the hb promoter. We have corroborated model noise predictions experimentally. The results indicate that WT (self-regulatory) Hb output noise is predominantly dependent on the transcription and translation dynamics of its own expression, rather than on Bcd fluctuations. The constructs and mutant, which lack self-regulation, indicate that the multiple Bcd binding sites in the hb promoter (and their strengths) also play a role in buffering noise. The model is robust to the variation in Bcd binding site number across a number of fly species. This study identifies particular ways in which promoter structure and regulatory dynamics reduce hb output noise. Insofar as many of these are common features of genes (e. g. multiple regulatory sites, cooperativity, self-feedback), the current results contribute to the general understanding of the reproducibility and determinacy of spatial patterning in early development.

VARIATION IN STOMATAL NUMBERS OF GLOSSOPTERIS LEAVES FROM THE LOWER PERMIAN OF PARANA BASIN, BRAZIL

The stomatal density and index in compressed leaves of Glossopteris communis from two different roof shales from the Lower Permian in Parana Basin, Brazil (Western Gondwana) have been investigated to test the possible relationship with modeled global changes in atmospheric CO(2) during the Phanerozoic. The obtained parameters show that the genus Glossopteris from the Cool Temperate biome can be used as CO(2) -proxy, despite the impossibility of being compared with living relatives or equivalents. When confronted with already published data for the Tropical Summer Wet biome, the present results confirm the detection of low levels of atmospheric CO(2) during the Early Permian, as predicted by the modeled curve. Nevertheless, the lower stomatal numbers detected at the climax of the coal interval (Faxinal Coalfield, Sakmarian) when compared to the higher ones obtained in leaves from a younger interval (Figueira Coalfield, Artinskian) could be attributed to temporarily high levels of atmospheric CO(2). Therefore, the occurrence of an extensive peat generating event at the southern part of the basin and subsequent greenhouse gases emissions from this environment may have been enough to reverse regionally and temporarily the reduction trend in atmospheric CO(2). Additionally, the Faxinal flora is preserved in a tonstein layer, which is a record of volcanic activity that could also cause a rise in atmospheric CO(2). During the Artinskian, the scarce generation of peat mires, as revealed by the occurrence of thin and discontinuous coal layers, and the lack of volcanism evidence would be insufficient to affect the general low CO(2) trend.

Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

Inference of gene regulatory networks from time series by Tsallis entropy

Background: The inference of gene regulatory networks (GRNs) from large-scale expression profiles is one of the most challenging problems of Systems Biology nowadays. Many techniques and models have been proposed for this task. However, it is not generally possible to recover the original topology with great accuracy, mainly due to the short time series data in face of the high complexity of the networks and the intrinsic noise of the expression measurements. In order to improve the accuracy of GRNs inference methods based on entropy (mutual information), a new criterion function is here proposed. Results: In this paper we introduce the use of generalized entropy proposed by Tsallis, for the inference of GRNs from time series expression profiles. The inference process is based on a feature selection approach and the conditional entropy is applied as criterion function. In order to assess the proposed methodology, the algorithm is applied to recover the network topology from temporal expressions generated by an artificial gene network (AGN) model as well as from the DREAM challenge. The adopted AGN is based on theoretical models of complex networks and its gene transference function is obtained from random drawing on the set of possible Boolean functions, thus creating its dynamics. On the other hand, DREAM time series data presents variation of network size and its topologies are based on real networks. The dynamics are generated by continuous differential equations with noise and perturbation. By adopting both data sources, it is possible to estimate the average quality of the inference with respect to different network topologies, transfer functions and network sizes. Conclusions: A remarkable improvement of accuracy was observed in the experimental results by reducing the number of false connections in the inferred topology by the non-Shannon entropy. The obtained best free parameter of the Tsallis entropy was on average in the range 2.5 <= q <= 3.5 (hence, subextensive entropy), which opens new perspectives for GRNs inference methods based on information theory and for investigation of the nonextensivity of such networks. The inference algorithm and criterion function proposed here were implemented and included in the DimReduction software, which is freely available at http://sourceforge.net/projects/dimreduction and http://code.google.com/p/dimreduction/.

Descriptions and phylogenetic significance of the fronto-lateral gland pores and dorsal lattice organs of cyprid larvae of seven species of barnacles (Cirripedia : Thoracica : Pedunculata)

The paired fronto-lateral gland pores and lattice organs (LO1, 2, 3, 4, and 5) of seven species of pedunculate barnacles belonging to two thoracican suborders, Heteralepadomorpha (family Heteralepadidae: Heteralepas sp. 1 and 2) and Lepadomorpha (families Poecilasmatidae: Poecilasma inaequilaterale and Octolasmis aymonini geryonophila and Lepadidae: Lepas pacifica, Dosima fascicularis, and Conchoderma virgatum), were investigated by scanning electron microscopy (SEM). While the fronto-lateral gland pores exhibit slight variation among species, with only L. pacifica showing a different morphology, the variations in the arrangement of LOs are phylogenetically instructive. The lattice organs in the foregoing species correspond in general to the inferred advanced type (Type C), but the distinct keel in the pore field in P. inaequilaterale and L. pacifica is reminiscent of, but not necessarily identical with the less advanced Type B. The arrangement of the anterior LOs (1-2) is rhomboidal in the two heteralepadomorph species, the two poecilasmatid species, and two of the three lepadid species, as it is in all previously and presently known lepadomorph cyprids except D. fascicularis. In this last species, they are deployed linearly along the hinge line. A linear arrangement of all the lattice organs is presumably the plesiomorphic condition for the Thoracica; an obvious exception being the pattern seen in Ibla cumingi. The arrangement of the first two pairs of posterior LOs (3-4) in O. a. geryonophila and C. virgatum differs from that of all previously described Lepadomorpha in being rhomboidal rather than aligned linearly along the hinge line. This same arrangement of LOs 3 and 4 in the two heteralepadomorph species is notable since it is not known in other thoracicans. Our results concerning variation in lattice organs of the lower Pedunculata are more or less consistent with current phylogenetic speculations and genetic information that ally Heteralepadomorpha with Lepadomorpha. Significance of this variation at lower taxonomic levels is also evident in the two similar forms of Heteralepas.

INTERANNUAL VARIABILITY IN TROPICAL ESTUARINE COPEPOD ASSEMBLAGES OFF NORTHEASTERN BRAZIL

We studied copepod assemblage variability among years, seasons, and tidal states in the Mucuri River estuary (Bahia State, Brazil). Zooplankton samples were collected seasonally through five years (2002-2006) at three sampling stations, one of which was sampled over a complete tidal cycle (two ebb and two flood tides). Temperature, salinity, river flux, and rainfall data were collected. Winter and summer represented dry and wet seasons, respectively. Copepod abundances ranged from 40 to 63% of the total zooplankton assemblage and comprised 46 taxa, among which, common estuarine species such as Temora turbinata (first record for the studied area), Parvocalanus crassirostris, Acartia lilljeborgi, Oithona hebes were the most abundant (euryhaline species). Interannual and seasonal variations were most marked in stenohaline species, e.g.. Notodiaptomus sp. and Thermocyclops minutus; density variations of euryhaline species, which made up the majority of the abundant taxa, were most closely related to tides. Diversity and richness also followed an intertidal pattern of variation.

Retinoic Acid-Treated Pluripotent Stem Cells Undergoing Neurogenesis Present Increased Aneuploidy and Micronuclei Formation

The existence of loss and gain of chromosomes, known as aneuploidy, has been previously described within the central nervous system. During development, at least one-third of neural progenitor cells (NPCs) are aneuploid. Notably, aneuploid NPCs may survive and functionally integrate into the mature neural circuitry. Given the unanswered significance of this phenomenon, we tested the hypothesis that neural differentiation induced by all-trans retinoic acid (RA) in pluripotent stem cells is accompanied by increased levels of aneuploidy, as previously described for cortical NPCs in vivo. In this work we used embryonal carcinoma (EC) cells, embryonic stem (ES) cells and induced pluripotent stem (iPS) cells undergoing differentiation into NPCs. Ploidy analysis revealed a 2-fold increase in the rate of aneuploidy, with the prevalence of chromosome loss in RA primed stem cells when compared to naive cells. In an attempt to understand the basis of neurogenic aneuploidy, micronuclei formation and survivin expression was assessed in pluripotent stem cells exposed to RA. RA increased micronuclei occurrence by almost 2-fold while decreased survivin expression by 50%, indicating possible mechanisms by which stem cells lose their chromosomes during neural differentiation. DNA fragmentation analysis demonstrated no increase in apoptosis on embryoid bodies treated with RA, indicating that cell death is not the mandatory fate of aneuploid NPCs derived from pluripotent cells. In order to exclude that the increase in aneuploidy was a spurious consequence of RA treatment, not related to neurogenesis, mouse embryonic fibroblasts were treated with RA under the same conditions and no alterations in chromosome gain or loss were observed. These findings indicate a correlation amongst neural differentiation, aneuploidy, micronuclei formation and survivin downregulation in pluripotent stem cells exposed to RA, providing evidence that somatically generated chromosomal variation accompanies neurogenesis in vitro.

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.

Retention-favorable timing of propagule release in barnacles and periwinkles

Recent studies on the ecology of marine larvae suggest that retention near coastal areas and self-recruitment are probably much more common than previously thought. In light of this, the advective potential of pelagic stages can be partially determined by the timing of propagule release. We sampled the upper-shore levels of a subtropical coastline in southeastern Brazil to examine the temporal patterns of propagule release for the common barnacle Chthamalus bisinuatus and the periwinkle Nodilittorina lineolata. The release timing in both species was very consistent between sites a few kilometers apart. Barnacles released nauplii in a rhythmic pattern, mostly coinciding with neap tides, when the speed of tidal currents was lowest. There was no variation in propagule release in relation to diel or flood-ebb tidal cycles. Periwinkles released propagules in a very irregular pattern, which remarkably matched a time series of wave heights. Egg capsules were released during periods of rough seas, when onshore surface currents were expected to prevail. Eggs were released mostly during ebb tides, and there was no diel variation. Propagule release rhythms are usually viewed as a means to ensure fast offshore advection, thus avoiding presumed hazardous conditions for larvae. However, propagule release strategies in these upper-shore species would more likely contribute to the retention of offspring close to parental stocks.