Showing papers in "BMC Biology in 2007"

Ab initio modeling of small proteins by iterative TASSER simulations

Abstract: Predicting 3-dimensional protein structures from amino-acid sequences is an important unsolved problem in computational structural biology. The problem becomes relatively easier if close homologous proteins have been solved, as high-resolution models can be built by aligning target sequences to the solved homologous structures. However, for sequences without similar folds in the Protein Data Bank (PDB) library, the models have to be predicted from scratch. Progress in the ab initio structure modeling is slow. The aim of this study was to extend the TASSER (threading/assembly/refinement) method for the ab initio modeling and examine systemically its ability to fold small single-domain proteins. We developed I-TASSER by iteratively implementing the TASSER method, which is used in the folding test of three benchmarks of small proteins. First, data on 16 small proteins (< 90 residues) were used to generate I-TASSER models, which had an average Cα-root mean square deviation (RMSD) of 3.8A, with 6 of them having a Cα-RMSD < 2.5A. The overall result was comparable with the all-atomic ROSETTA simulation, but the central processing unit (CPU) time by I-TASSER was much shorter (150 CPU days vs. 5 CPU hours). Second, data on 20 small proteins (< 120 residues) were used. I-TASSER folded four of them with a Cα-RMSD < 2.5A. The average Cα-RMSD of the I-TASSER models was 3.9A, whereas it was 5.9A using TOUCHSTONE-II software. Finally, 20 non-homologous small proteins (< 120 residues) were taken from the PDB library. An average Cα-RMSD of 3.9A was obtained for the third benchmark, with seven cases having a Cα-RMSD < 2.5A. Our simulation results show that I-TASSER can consistently predict the correct folds and sometimes high-resolution models for small single-domain proteins. Compared with other ab initio modeling methods such as ROSETTA and TOUCHSTONE II, the average performance of I-TASSER is either much better or is similar within a lower computational time. These data, together with the significant performance of automated I-TASSER server (the Zhang-Server) in the 'free modeling' section of the recent Critical Assessment of Structure Prediction (CASP)7 experiment, demonstrate new progresses in automated ab initio model generation. The I-TASSER server is freely available for academic users http://zhang.bioinformatics.ku.edu/I-TASSER .

Late-acting dominant lethal genetic systems and mosquito control

Abstract: Reduction or elimination of vector populations will tend to reduce or eliminate transmission of vector-borne diseases. One potential method for environmentally-friendly, species-specific population control is the Sterile Insect Technique (SIT). SIT has not been widely used against insect disease vectors such as mosquitoes, in part because of various practical difficulties in rearing, sterilization and distribution. Additionally, vector populations with strong density-dependent effects will tend to be resistant to SIT-based control as the population-reducing effect of induced sterility will tend to be offset by reduced density-dependent mortality. We investigated by mathematical modeling the effect of manipulating the stage of development at which death occurs (lethal phase) in an SIT program against a density-dependence-limited insect population. We found late-acting lethality to be considerably more effective than early-acting lethality. No such strains of a vector insect have been described, so as a proof-of-principle we constructed a strain of the principal vector of the dengue and yellow fever viruses, Aedes (Stegomyia) aegypti, with the necessary properties of dominant, repressible, highly penetrant, late-acting lethality. Conventional SIT induces early-acting (embryonic) lethality, but genetic methods potentially allow the lethal phase to be tailored to the program. For insects with strong density-dependence, we show that lethality after the density-dependent phase would be a considerable improvement over conventional methods. For density-dependent parameters estimated from field data for Aedes aegypti, the critical release ratio for population elimination is modeled to be 27% to 540% greater for early-acting rather than late-acting lethality. Our success in developing a mosquito strain with the key features that the modeling indicated were desirable demonstrates the feasibility of this approach for improved SIT for disease control.

Classification and nomenclature of all human homeobox genes

Abstract: The homeobox genes are a large and diverse group of genes, many of which play important roles in the embryonic development of animals. Increasingly, homeobox genes are being compared between genomes in an attempt to understand the evolution of animal development. Despite their importance, the full diversity of human homeobox genes has not previously been described. We have identified all homeobox genes and pseudogenes in the euchromatic regions of the human genome, finding many unannotated, incorrectly annotated, unnamed, misnamed or misclassified genes and pseudogenes. We describe 300 human homeobox loci, which we divide into 235 probable functional genes and 65 probable pseudogenes. These totals include 3 genes with partial homeoboxes and 13 pseudogenes that lack homeoboxes but are clearly derived from homeobox genes. These figures exclude the repetitive DUX1 to DUX5 homeobox sequences of which we identified 35 probable pseudogenes, with many more expected in heterochromatic regions. Nomenclature is established for approximately 40 formerly unnamed loci, reflecting their evolutionary relationships to other loci in human and other species, and nomenclature revisions are proposed for around 30 other loci. We use a classification that recognizes 11 homeobox gene 'classes' subdivided into 102 homeobox gene 'families'. We have conducted a comprehensive survey of homeobox genes and pseudogenes in the human genome, described many new loci, and revised the classification and nomenclature of homeobox genes. The classification scheme may be widely applicable to homeobox genes in other animal genomes and will facilitate comparative genomics of this important gene superclass.

Broadening the horizon – level 2.5 of the HUPO-PSI format for molecular interactions

Abstract: Molecular interaction Information is a key resource in modern biomedical research. Publicly available data have previously been provided in a broad array of diverse formats, making access to this very difficult. The publication and wide implementation of the Human Proteome Organisation Proteomics Standards Initiative Molecular Interactions (HUPO PSI-MI) format in 2004 was a major step towards the establishment of a single, unified format by which molecular interactions should be presented, but focused purely on protein-protein interactions. The HUPO-PSI has further developed the PSI-MI XML schema to enable the description of interactions between a wider range of molecular types, for example nucleic acids, chemical entities, and molecular complexes. Extensive details about each supported molecular interaction can now be captured, including the biological role of each molecule within that interaction, detailed description of interacting domains, and the kinetic parameters of the interaction. The format is supported by data management and analysis tools and has been adopted by major interaction data providers. Additionally, a simpler, tab-delimited format MITAB2.5 has been developed for the benefit of users who require only minimal information in an easy to access configuration. The PSI-MI XML2.5 and MITAB2.5 formats have been jointly developed by interaction data producers and providers from both the academic and commercial sector, and are already widely implemented and well supported by an active development community. PSI-MI XML2.5 enables the description of highly detailed molecular interaction data and facilitates data exchange between databases and users without loss of information. MITAB2.5 is a simpler format appropriate for fast Perl parsing or loading into Microsoft Excel.

A 100%-complete sequence reveals unusually simple genomic features in the hot-spring red alga Cyanidioschyzon merolae

Abstract: All previously reported eukaryotic nuclear genome sequences have been incomplete, especially in highly repeated units and chromosomal ends. Because repetitive DNA is important for many aspects of biology, complete chromosomal structures are fundamental for understanding eukaryotic cells. Our earlier, nearly complete genome sequence of the hot-spring red alga Cyanidioschyzon merolae revealed several unique features, including just three ribosomal DNA copies, very few introns, and a small total number of genes. However, because the exact structures of certain functionally important repeated elements remained ambiguous, that sequence was not complete. Obviously, those ambiguities needed to be resolved before the unique features of the C. merolae genome could be summarized, and the ambiguities could only be resolved by completing the sequence. Therefore, we aimed to complete all previous gaps and sequence all remaining chromosomal ends, and now report the first nuclear-genome sequence for any eukaryote that is 100% complete. Our present complete sequence consists of 16546747 nucleotides covering 100% of the 20 linear chromosomes from telomere to telomere, representing the simple and unique chromosomal structures of the eukaryotic cell. We have unambiguously established that the C. merolae genome contains the smallest known histone-gene cluster, a unique telomeric repeat for all chromosomal ends, and an extremely low number of transposons. By virtue of these attributes and others that we had discovered previously, C. merolae appears to have the simplest nuclear genome of the non-symbiotic eukaryotes. These unusually simple genomic features in the 100% complete genome sequence of C. merolae are extremely useful for further studies of eukaryotic cells.

Nutritional upgrading for omnivorous carpenter ants by the endosymbiont Blochmannia

Abstract: Carpenter ants (genus Camponotus) are considered to be omnivores. Nonetheless, the genome sequence of Blochmannia floridanus, the obligate intracellular endosymbiont of Camponotus floridanus, suggests a function in nutritional upgrading of host resources by the bacterium. Thus, the strongly reduced genome of the endosymbiont retains genes for all subunits of a functional urease, as well as those for biosynthetic pathways for all but one (arginine) of the amino acids essential to the host. Nutritional upgrading by Blochmannia was tested in 90-day feeding experiments with brood-raising in worker-groups on chemically defined diets with and without essential amino acids and treated or not with antibiotics. Control groups were fed with cockroaches, honey water and Bhatkar agar. Worker-groups were provided with brood collected from the queenright mother-colonies (45 eggs and 45 first instar larvae each). Brood production did not differ significantly between groups of symbiotic workers on diets with and without essential amino acids. However, aposymbiotic worker groups raised significantly less brood on a diet lacking essential amino acids. Reduced brood production by aposymbiotic workers was compensated when those groups were provided with essential amino acids in their diet. Decrease of endosymbionts due to treatment with antibiotic was monitored by qRT-PCR and FISH after the 90-day experimental period. Urease function was confirmed by feeding experiments using 15N-labelled urea. GC-MS analysis of 15N-enrichment of free amino acids in workers revealed significant labelling of the non-essential amino acids alanine, glycine, aspartic acid, and glutamic acid, as well as of the essential amino acids methionine and phenylalanine. Our results show that endosymbiotic Blochmannia nutritionally upgrade the diet of C. floridanus hosts to provide essential amino acids, and that it may also play a role in nitrogen recycling via its functional urease. Blochmannia may confer a significant fitness advantage via nutritional upgrading by enhancing competitive ability of Camponotus with other ant species lacking such an endosymbiont. Domestication of the endosymbiont may have facilitated the evolutionary success of the genus Camponotus.

Faced with inequality: chicken do not have a general dosage compensation of sex-linked genes

Abstract: The contrasting dose of sex chromosomes in males and females potentially introduces a large-scale imbalance in levels of gene expression between sexes, and between sex chromosomes and autosomes. In many organisms, dosage compensation has thus evolved to equalize sex-linked gene expression in males and females. In mammals this is achieved by X chromosome inactivation and in flies and worms by up- or down-regulation of X-linked expression, respectively. While otherwise widespread in systems with heteromorphic sex chromosomes, the case of dosage compensation in birds (males ZZ, females ZW) remains an unsolved enigma.

Extensive population genetic structure in the giraffe

Abstract: A central question in the evolutionary diversification of large, widespread, mobile mammals is how substantial differentiation can arise, particularly in the absence of topographic or habitat barriers to dispersal. All extant giraffes (Giraffa camelopardalis) are currently considered to represent a single species classified into multiple subspecies. However, geographic variation in traits such as pelage pattern is clearly evident across the range in sub-Saharan Africa and abrupt transition zones between different pelage types are typically not associated with extrinsic barriers to gene flow, suggesting reproductive isolation. By analyzing mitochondrial DNA sequences and nuclear microsatellite loci, we show that there are at least six genealogically distinct lineages of giraffe in Africa, with little evidence of interbreeding between them. Some of these lineages appear to be maintained in the absence of contemporary barriers to gene flow, possibly by differences in reproductive timing or pelage-based assortative mating, suggesting that populations usually recognized as subspecies have a long history of reproductive isolation. Further, five of the six putative lineages also contain genetically discrete populations, yielding at least 11 genetically distinct populations. Such extreme genetic subdivision within a large vertebrate with high dispersal capabilities is unprecedented and exceeds that of any other large African mammal. Our results have significant implications for giraffe conservation, and imply separate in situ and ex situ management, not only of pelage morphs, but also of local populations.

The vertebrate phylotypic stage and an early bilaterian-related stage in mouse embryogenesis defined by genomic information

Abstract: Embryos of taxonomically different vertebrates are thought to pass through a stage in which they resemble one another morphologically. This "vertebrate phylotypic stage" may represent the basic vertebrate body plan that was established in the common ancestor of vertebrates. However, much controversy remains about when the phylotypic stage appears, and whether it even exists. To overcome the limitations of studies based on morphological comparison, we explored a comprehensive quantitative method for defining the constrained stage using expressed sequence tag (EST) data, gene ontologies (GO), and available genomes of various animals. If strong developmental constraints occur during the phylotypic stage of vertebrate embryos, then genes conserved among vertebrates would be highly expressed at this stage. We established a novel method for evaluating the ancestral nature of mouse embryonic stages that does not depend on comparative morphology. The numerical "ancestor index" revealed that the mouse indeed has a highly conserved embryonic period at embryonic day 8.0–8.5, the time of appearance of the pharyngeal arch and somites. During this period, the mouse prominently expresses GO-determined developmental genes shared among vertebrates. Similar analyses revealed the existence of a bilaterian-related period, during which GO-determined developmental genes shared among bilaterians are markedly expressed at the cleavage-to-gastrulation period. The genes associated with the phylotypic stage identified by our method are essential in embryogenesis. Our results demonstrate that the mid-embryonic stage of the mouse is indeed highly constrained, supporting the existence of the phylotypic stage. Furthermore, this candidate stage is preceded by a putative bilaterian ancestor-related period. These results not only support the developmental hourglass model, but also highlight the hierarchical aspect of embryogenesis proposed by von Baer. Identification of conserved stages and tissues by this method in various animals would be a powerful tool to examine the phylotypic stage hypothesis, and to understand which kinds of developmental events and gene sets are evolutionarily constrained and how they limit the possible variations of animal basic body plans.

Immune system responses and fitness costs associated with consumption of bacteria in larvae of Trichoplusia ni.

Abstract: Insects helped pioneer, and persist as model organisms for, the study of specific aspects of immunity Although they lack an adaptive immune system, insects possess an innate immune system that recognizes and destroys intruding microorganisms Its operation under natural conditions has not been well studied, as most studies have introduced microbes to laboratory-reared insects via artificial mechanical wounding One of the most common routes of natural exposure and infection, however, is via food; thus, the role of dietary microbial communities in herbivorous insect immune system evolution invites study Here, we examine the immune system response and consequences of exposing a lepidopteran agricultural pest to non-infectious microorganisms via simple oral consumption Immune system response was compared between Trichoplusia ni larvae reared on diets with or without non-pathogenic bacteria (Escherichia coli and Micrococcus luteus) Two major immune response-related enzymatic activities responded to diets differently – phenoloxidase activity was inhibited in the bacteria-fed larvae, whereas general antibacterial activity was enhanced Eight proteins were highly expressed in the hemolymph of the bacteria fed larvae, among them immune response related proteins arylphorin, apolipophorin III and gloverin Expression response among 25 putative immune response-related genes were assayed via RT-qPCR Seven showed more than fivefold up regulation in the presence of bacterial diet, with 22 in total being differentially expressed, among them apolipophorin III, cecropin, gallerimycin, gloverin, lysozyme, and phenoloxidase inhibiting enzyme Finally, potential life-history trade-offs were studied, with pupation time and pupal mass being negatively affected in bacteria fed larvae The presence of bacteria in food, even if non-pathogenic, can trigger an immune response cascade with life history tradeoffs Trichoplusia ni larvae are able to detect and respond to environmental microbes encountered in the diet, possibly even using midgut epithelial tissue as a sensing organ Potential benefits of this immune system priming may outweigh the observed tradeoffs, as priming based on environmentally sensed bacterial may decrease risk of serious infection These results show that food plant microbial communities represent a dynamic and unstudied part of the coevolutionary interactions between plants and their insect herbivores

Design and applicability of DNA arrays and DNA barcodes in biodiversity monitoring

Abstract: The rapid and accurate identification of species is a critical component of large-scale biodiversity monitoring programs. DNA arrays (micro and macro) and DNA barcodes are two molecular approaches that have recently garnered much attention. Here, we compare these two platforms for identification of an important group, the mammals. Our analyses, based on the two commonly used mitochondrial genes cytochrome c oxidase I (the standard DNA barcode for animal species) and cytochrome b (a common species-level marker), suggest that both arrays and barcodes are capable of discriminating mammalian species with high accuracy. We used three different datasets of mammalian species, comprising different sampling strategies. For DNA arrays we designed three probes for each species to address intraspecific variation. As for DNA barcoding, our analyses show that both cytochrome c oxidase I and cytochrome b genes, and even smaller fragments of them (mini-barcodes) can successfully discriminate species in a wide variety of specimens. This study showed that DNA arrays and DNA barcodes are valuable molecular methods for biodiversity monitoring programs. Both approaches were capable of discriminating among mammalian species in our test assemblages. However, because designing DNA arrays require advance knowledge of target sequences, the use of this approach could be limited in large scale monitoring programs where unknown haplotypes might be encountered. DNA barcodes, by contrast, are sequencing-based and therefore could provide more flexibility in large-scale studies.

Rise of oceanographic barriers in continuous populations of a cetacean: the genetic structure of harbour porpoises in Old World waters.

Abstract: Understanding the role of seascape in shaping genetic and demographic population structure is highly challenging for marine pelagic species such as cetaceans for which there is generally little evidence of what could effectively restrict their dispersal. In the present work, we applied a combination of recent individual-based landscape genetic approaches to investigate the population genetic structure of a highly mobile extensive range cetacean, the harbour porpoise in the eastern North Atlantic, with regards to oceanographic characteristics that could constrain its dispersal. Analyses of 10 microsatellite loci for 752 individuals revealed that most of the sampled range in the eastern North Atlantic behaves as a 'continuous' population that widely extends over thousands of kilometres with significant isolation by distance (IBD). However, strong barriers to gene flow were detected in the south-eastern part of the range. These barriers coincided with profound changes in environmental characteristics and isolated, on a relatively small scale, porpoises from Iberian waters and on a larger scale porpoises from the Black Sea. The presence of these barriers to gene flow that coincide with profound changes in oceanographic features, together with the spatial variation in IBD strength, provide for the first time strong evidence that physical processes have a major impact on the demographic and genetic structure of a cetacean. This genetic pattern further suggests habitat-related fragmentation of the porpoise range that is likely to intensify with predicted surface ocean warming.

Pair of lice lost or parasites regained: the evolutionary history of anthropoid primate lice

Abstract: The parasitic sucking lice of primates are known to have undergone at least 25 million years of coevolution with their hosts. For example, chimpanzee lice and human head/body lice last shared a common ancestor roughly six million years ago, a divergence that is contemporaneous with their hosts. In an assemblage where lice are often highly host specific, humans host two different genera of lice, one that is shared with chimpanzees and another that is shared with gorillas. In this study, we reconstruct the evolutionary history of primate lice and infer the historical events that explain the current distribution of these lice on their primate hosts. Phylogenetic and cophylogenetic analyses suggest that the louse genera Pediculus and Pthirus are each monophyletic, and are sister taxa to one another. The age of the most recent common ancestor of the two Pediculus species studied matches the age predicted by host divergence (ca. 6 million years), whereas the age of the ancestor of Pthirus does not. The two species of Pthirus (Pthirus gorillae and Pthirus pubis) last shared a common ancestor ca. 3–4 million years ago, which is considerably younger than the divergence between their hosts (gorillas and humans, respectively), of approximately 7 million years ago. Reconciliation analysis determines that there are two alternative explanations that account for the current distribution of anthropoid primate lice. The more parsimonious of the two solutions suggests that a Pthirus species switched from gorillas to humans. This analysis assumes that the divergence between Pediculus and Pthirus was contemporaneous with the split (i.e., a node of cospeciation) between gorillas and the lineage leading to chimpanzees and humans. Divergence date estimates, however, show that the nodes in the host and parasite trees are not contemporaneous. Rather, the shared coevolutionary history of the anthropoid primates and their lice contains a mixture of evolutionary events including cospeciation, parasite duplication, parasite extinction, and host switching. Based on these data, the coevolutionary history of primates and their lice has been anything but parsimonious.

Identification of novel conserved peptide uORF homology groups in Arabidopsis and rice reveals ancient eukaryotic origin of select groups and preferential association with transcription factor-encoding genes

Abstract: Upstream open reading frames (uORFs) can mediate translational control over the largest, or major ORF (mORF) in response to starvation, polyamine concentrations, and sucrose concentrations. One plant uORF with conserved peptide sequences has been shown to exert this control in an amino acid sequence-dependent manner but generally it is not clear what kinds of genes are regulated, or how extensively this mechanism is invoked in a given genome. By comparing full-length cDNA sequences from Arabidopsis and rice we identified 26 distinct homology groups of conserved peptide uORFs, only three of which have been reported previously. Pairwise K a /K s analysis showed that purifying selection had acted on nearly all conserved peptide uORFs and their associated mORFs. Functions of predicted mORF proteins could be inferred for 16 homology groups and many of these proteins appear to have a regulatory function, including 6 transcription factors, 5 signal transduction factors, 3 developmental signal molecules, a homolog of translation initiation factor eIF5, and a RING finger protein. Transcription factors are clearly overrepresented in this data set when compared to the frequency calculated for the entire genome (p = 1.2 × 10-7). Duplicate gene pairs arising from a whole genome duplication (ohnologs) with a conserved uORF are much more likely to have been retained in Arabidopsis (Arabidopsis thaliana) than are ohnologs of other genes (39% vs 14% of ancestral genes, p = 5 × 10-3). Two uORF groups were found in animals, indicating an ancient origin of these putative regulatory elements. Conservation of uORF amino acid sequence, association with homologous mORFs over long evolutionary time periods, preferential retention after whole genome duplications, and preferential association with mORFs coding for transcription factors suggest that the conserved peptide uORFs identified in this study are strong candidates for translational controllers of regulatory genes.

Global and regional brain metabolic scaling and its functional consequences

Abstract: Background Information processing in the brain requires large amounts of metabolic energy, the spatial distribution of which is highly heterogeneous, reflecting the complex activity patterns in the mammalian brain.

Real-time phase-contrast x-ray imaging: a new technique for the study of animal form and function

Abstract: Despite advances in imaging techniques, real-time visualization of the structure and dynamics of tissues and organs inside small living animals has remained elusive. Recently, we have been using synchrotron x-rays to visualize the internal anatomy of millimeter-sized opaque, living animals. This technique takes advantage of partially-coherent x-rays and diffraction to enable clear visualization of internal soft tissue not viewable via conventional absorption radiography. However, because higher quality images require greater x-ray fluxes, there exists an inherent tradeoff between image quality and tissue damage. We evaluated the tradeoff between image quality and harm to the animal by determining the impact of targeted synchrotron x-rays on insect physiology, behavior and survival. Using 25 keV x-rays at a flux density of 80 μW/mm-2, high quality video-rate images can be obtained without major detrimental effects on the insects for multiple minutes, a duration sufficient for many physiological studies. At this setting, insects do not heat up. Additionally, we demonstrate the range of uses of synchrotron phase-contrast imaging by showing high-resolution images of internal anatomy and observations of labeled food movement during ingestion and digestion. Synchrotron x-ray phase contrast imaging has the potential to revolutionize the study of physiology and internal biomechanics in small animals. This is the only generally applicable technique that has the necessary spatial and temporal resolutions, penetrating power, and sensitivity to soft tissue that is required to visualize the internal physiology of living animals on the scale from millimeters to microns.

Accelerated gene evolution and subfunctionalization in the pseudotetraploid frog Xenopus laevis

Abstract: Ancient whole genome duplications have been implicated in the vertebrate and teleost radiations, and in the emergence of diverse angiosperm lineages, but the evolutionary response to such a perturbation is still poorly understood. The African clawed frog Xenopus laevis experienced a relatively recent tetraploidization ~40 million years ago. Analysis of the considerable amount of EST sequence available for this species together with the genome sequence of the related diploid Xenopus tropicalis provides a unique opportunity to study the genomic response to whole genome duplication. We identified 2218 gene triplets in which a single gene in X. tropicalis corresponds to precisely two co-orthologous genes in X. laevis – the largest such collection published from any duplication event in animals. Analysis of these triplets reveals accelerated evolution or relaxation of constraint in the peptides of the X. laevis pairs compared with the orthologous sequences in X. tropicalis and other vertebrates. In contrast, single-copy X. laevis genes do not show this acceleration. Duplicated genes can differ substantially in expression levels and patterns. We find no significant difference in gene content in the duplicated set, versus the single-copy set based on molecular and biological function ontologies. These results support a scenario in which duplicate genes are retained through a process of subfunctionalization and/or relaxation of constraint on both copies of an ancestral gene.

Male urine signals social rank in the Mozambique tilapia (Oreochromis mossambicus).

Abstract: Background: The urine of freshwater fish species investigated so far acts as a vehicle for reproductive pheromones affecting the behaviour and physiology of the opposite sex. However, the role of urinary pheromones in intra-sexual competition has received less attention. This is particularly relevant in lek-breeding species, such as the Mozambique tilapia (Oreochromis mossambicus), where males establish dominance hierarchies and there is the possibility for chemical communication in the modulation of aggression among males. To investigate whether males use urine during aggressive interactions, we measured urination frequency of dye-injected males during paired interactions between size-matched males. Furthermore, we assessed urinary volume stored in the bladder of males in a stable social hierarchy and the olfactory potency of their urine by recording of the electro-olfactogram. Results: Males released urine in pulses of short duration (about one second) and markedly increased urination frequency during aggressive behaviour, but did not release urine whilst submissive. In the stable hierarchy, subordinate males stored less urine than males of higher social rank; the olfactory potency of the urine was positively correlated with the rank of the male donor. Conclusion: Dominant males store urine and use it as a vehicle for odorants actively released during aggressive disputes. The olfactory potency of the urine is positively correlated with the social status of the male. We suggest that males actively advertise their dominant status through urinary odorants which may act as a 'dominance' pheromone to modulate aggression in rivals, thereby contributing to social stability within the lek.

Evolution favors protein mutational robustness in sufficiently large populations

Abstract: An important question is whether evolution favors properties such as mutational robustness or evolvability that do not directly benefit any individual, but can influence the course of future evolution. Functionally similar proteins can differ substantially in their robustness to mutations and capacity to evolve new functions, but it has remained unclear whether any of these differences might be due to evolutionary selection for these properties. Here we use laboratory experiments to demonstrate that evolution favors protein mutational robustness if the evolving population is sufficiently large. We neutrally evolve cytochrome P450 proteins under identical selection pressures and mutation rates in populations of different sizes, and show that proteins from the larger and thus more polymorphic population tend towards higher mutational robustness. Proteins from the larger population also evolve greater stability, a biophysical property that is known to enhance both mutational robustness and evolvability. The excess mutational robustness and stability is well described by mathematical theory, and can be quantitatively related to the way that the proteins occupy their neutral network. Our work is the first experimental demonstration of the general tendency of evolution to favor mutational robustness and protein stability in highly polymorphic populations. We suggest that this phenomenon could contribute to the mutational robustness and evolvability of viruses and bacteria that exist in large populations.

Adaptive sequence evolution in a color gene involved in the formation of the characteristic egg-dummies of male haplochromine cichlid fishes

Abstract: The exceptionally diverse species flocks of cichlid fishes in East Africa are prime examples of parallel adaptive radiations. About 80% of East Africa's more than 1 800 endemic cichlid species, and all species of the flocks of Lakes Victoria and Malawi, belong to a particularly rapidly evolving lineage, the haplochromines. One characteristic feature of the haplochromines is their possession of egg-dummies on the males' anal fins. These egg-spots mimic real eggs and play an important role in the mating system of these maternal mouthbrooding fish. Here, we show that the egg-spots of haplochromines are made up of yellow pigment cells, xanthophores, and that a gene coding for a type III receptor tyrosine kinase, colony-stimulating factor 1 receptor a (csf1ra), is expressed in egg-spot tissue. Molecular evolutionary analyses reveal that the extracellular ligand-binding and receptor-interacting domain of csf1ra underwent adaptive sequence evolution in the ancestral lineage of the haplochromines, coinciding with the emergence of egg-dummies. We also find that csf1ra is expressed in the egg-dummies of a distantly related cichlid species, the ectodine cichlid Ophthalmotilapia ventralis, in which markings with similar functions evolved on the pelvic fin in convergence to those of the haplochromines. We conclude that modifications of existing signal transduction mechanisms might have evolved in the haplochromine lineage in association with the origination of anal fin egg-dummies. That positive selection has acted during the evolution of a color gene that seems to be involved in the morphogenesis of a sexually selected trait, the egg-dummies, highlights the importance of further investigations of the comparative genomic basis of the phenotypic diversification of cichlid fishes.

A clade uniting the green algae Mesostigma viride and Chlorokybus atmophyticus represents the deepest branch of the Streptophyta in chloroplast genome-based phylogenies

Abstract: The Viridiplantae comprise two major phyla: the Streptophyta, containing the charophycean green algae and all land plants, and the Chlorophyta, containing the remaining green algae. Despite recent progress in unravelling phylogenetic relationships among major green plant lineages, problematic nodes still remain in the green tree of life. One of the major issues concerns the scaly biflagellate Mesostigma viride, which is either regarded as representing the earliest divergence of the Streptophyta or a separate lineage that diverged before the Chlorophyta and Streptophyta. Phylogenies based on chloroplast and mitochondrial genomes support the latter view. Because some green plant lineages are not represented in these phylogenies, sparse taxon sampling has been suspected to yield misleading topologies. Here, we describe the complete chloroplast DNA (cpDNA) sequence of the early-diverging charophycean alga Chlorokybus atmophyticus and present chloroplast genome-based phylogenies with an expanded taxon sampling. The 152,254 bp Chlorokybus cpDNA closely resembles its Mesostigma homologue at the gene content and gene order levels. Using various methods of phylogenetic inference, we analyzed amino acid and nucleotide data sets that were derived from 45 protein-coding genes common to the cpDNAs of 37 green algal/land plant taxa and eight non-green algae. Unexpectedly, all best trees recovered a robust clade uniting Chlorokybus and Mesostigma. In protein trees, this clade was sister to all streptophytes and chlorophytes and this placement received moderate support. In contrast, gene trees provided unequivocal support to the notion that the Mesostigma + Chlorokybus clade represents the earliest-diverging branch of the Streptophyta. Independent analyses of structural data (gene content and/or gene order) and of subsets of amino acid data progressively enriched in slow-evolving sites led us to conclude that the latter topology reflects the true organismal relationships. In disclosing a sister relationship between the Mesostigmatales and Chlorokybales, our study resolves the long-standing debate about the nature of the unicellular flagellated ancestors of land plants and alters significantly our concepts regarding the evolution of streptophyte algae. Moreover, in predicting a richer chloroplast gene repertoire than previously inferred for the common ancestor of all streptophytes, our study has contributed to a better understanding of chloroplast genome evolution in the Viridiplantae.

Sleep in Kcna2 knockout mice

Abstract: Shaker codes for a Drosophila voltage-dependent potassium channel. Flies carrying Shaker null or hypomorphic mutations sleep 3–4 h/day instead of 8–14 h/day as their wild-type siblings do. Shaker-like channels are conserved across species but it is unknown whether they affect sleep in mammals. To address this issue, we studied sleep in Kcna2 knockout (KO) mice. Kcna2 codes for Kv1.2, the alpha subunit of a Shaker-like voltage-dependent potassium channel with high expression in the mammalian thalamocortical system. Continuous (24 h) electroencephalograph (EEG), electromyogram (EMG), and video recordings were used to measure sleep and waking in Kcna2 KO, heterozygous (HZ) and wild-type (WT) pups (P17) and HZ and WT adult mice (P67). Sleep stages were scored visually based on 4-s epochs. EEG power spectra (0–20 Hz) were calculated on consecutive 4-s epochs. KO pups die by P28 due to generalized seizures. At P17 seizures are either absent or very rare in KO pups (< 1% of the 24-h recording time), and abnormal EEG activity is only present during the seizure. KO pups have significantly less non-rapid eye movement (NREM) sleep (-23%) and significantly more waking (+21%) than HZ and WT siblings with no change in rapid eye movement (REM) sleep time. The decrease in NREM sleep is due to an increase in the number of waking episodes, with no change in number or duration of sleep episodes. Sleep patterns, daily amounts of sleep and waking, and the response to 6 h sleep deprivation are similar in HZ and WT adult mice. Kv1.2, a mammalian homologue of Shaker, regulates neuronal excitability and affects NREM sleep.

Species status of Neisseria gonorrhoeae: evolutionary and epidemiological inferences from multilocus sequence typing

Abstract: Various typing methods have been developed for Neisseria gonorrhoeae, but none provide the combination of discrimination, reproducibility, portability, and genetic inference that allows the analysis of all aspects of the epidemiology of this pathogen from a single data set. Multilocus sequence typing (MLST) has been used successfully to characterize the related organisms Neisseria meningitidis and Neisseria lactamica. Here, the same seven locus Neisseria scheme was used to characterize a diverse collection of N. gonorrhoeae isolates to investigate whether this method would allow differentiation among isolates, and to distinguish these three species. A total of 149 gonococcal isolates were typed and submitted to the Neisseria MLST database. Although relatively few (27) polymorphisms were detected among the seven MLST loci, a total of 66 unique allele combinations (sequence types, STs), were observed, a number comparable to that seen among isolate collections of the more diverse meningococcus. Patterns of genetic variation were consistent with high levels of recombination generating this diversity. There was no evidence for geographical structuring among the isolates examined, with isolates collected in Liverpool, UK, showing levels of diversity similar to a global collection of isolates. There was, however, evidence that populations of N. meningitidis, N. gonorrhoeae and N. lactamica were distinct, with little support for frequent genetic recombination among these species, with the sequences from the gdh locus alone grouping the species into distinct clusters. The seven loci Neisseria MLST scheme was readily adapted to N. gonorrhoeae isolates, providing a highly discriminatory typing method. In addition, these data permitted phylogenetic and population genetic inferences to be made, including direct comparisons with N. meningitidis and N. lactamica. Examination of these data demonstrated that alleles were rarely shared among the three species. Analysis of variation at a single locus, gdh, provided a rapid means of identifying misclassified isolates and determining whether mixed cultures were present.

Systematics and plastid genome evolution of the cryptically photosynthetic parasitic plant genus Cuscuta (Convolvulaceae)

Abstract: The genus Cuscuta L. (Convolvulaceae), commonly known as dodders, are epiphytic vines that invade the stems of their host with haustorial feeding structures at the points of contact. Although they lack expanded leaves, some species are noticeably chlorophyllous, especially as seedlings and in maturing fruits. Some species are reported as crop pests of worldwide distribution, whereas others are extremely rare and have local distributions and apparent niche specificity. A strong phylogenetic framework for this large genus is essential to understand the interesting ecological, morphological and molecular phenomena that occur within these parasites in an evolutionary context. Here we present a well-supported phylogeny of Cuscuta using sequences of the nuclear ribosomal internal transcribed spacer and plastid rps2, rbcL and matK from representatives across most of the taxonomic diversity of the genus. We use the phylogeny to interpret morphological and plastid genome evolution within the genus. At least three currently recognized taxonomic sections are not monophyletic and subgenus Cuscuta is unequivocally paraphyletic. Plastid genes are extremely variable with regards to evolutionary constraint, with rbcL exhibiting even higher levels of purifying selection in Cuscuta than photosynthetic relatives. Nuclear genome size is highly variable within Cuscuta, particularly within subgenus Grammica, and in some cases may indicate the existence of cryptic species in this large clade of morphologically similar species. Some morphological characters traditionally used to define major taxonomic splits within Cuscuta are homoplastic and are of limited use in defining true evolutionary groups. Chloroplast genome evolution seems to have evolved in a punctuated fashion, with episodes of loss involving suites of genes or tRNAs followed by stabilization of gene content in major clades. Nearly all species of Cuscuta retain some photosynthetic ability, most likely for nutrient apportionment to their seeds, while complete loss of photosynthesis and possible loss of the entire chloroplast genome is limited to a single small clade of outcrossing species found primarily in western South America.

Pygo1 and Pygo2 roles in Wnt signaling in mammalian kidney development

Abstract: The pygopus gene of Drosophila encodes an essential component of the Armadillo (β-catenin) transcription factor complex of canonical Wnt signaling. To better understand the functions of Pygopus-mediated canonical Wnt signaling in kidney development, targeted mutations were made in the two mammalian orthologs, Pygo1 and Pygo2. Each mutation deleted >80% of the coding sequence, including the critical PHD domain, and almost certainly resulted in null function. Pygo2 homozygous mutants, with rare exception, died shortly after birth, with a phenotype including lens agenesis, growth retardation, altered kidney development, and in some cases exencephaly and cleft palate. Pygo1 homozygous mutants, however, were viable and fertile, with no detectable developmental defects. Double Pygo1/Pygo2 homozygous mutants showed no apparent synergy in phenotype severity. The BAT-gal transgene reporter of canonical Wnt signaling showed reduced levels of expression in Pygo1 -/-/Pygo2 -/- mutants, with tissue-specific variation in degree of diminution. The Pygo1 and Pygo2 genes both showed widespread expression in the developing kidney, with raised levels in the stromal cell compartment. Confocal analysis of the double mutant kidneys showed disturbance of both the ureteric bud and metanephric mesenchyme-derived compartments. Branching morphogenesis of the ureteric bud was altered, with expanded tips and reduced tip density, probably contributing to the smaller size of the mutant kidney. In addition, there was an expansion of the zone of condensed mesenchyme capping the ureteric bud. Nephron formation, however, proceeded normally. Microarray analysis showed changed expression of several genes, including Cxcl13, Slc5a2, Klk5, Ren2 and Timeless, which represent candidate Wnt targets in kidney development. The mammalian Pygopus genes are required for normal branching morphogenesis of the ureteric bud during kidney development. Nevertheless, the relatively mild phenotype observed in the kidney, as well as other organ systems, indicates a striking evolutionary divergence of Pygopus function between mammals and Drosophila. In mammals, the Pygo1/Pygo2 genes are not absolutely required for canonical Wnt signaling in most developing systems, but rather function as quantitative transducers, or modulators, of Wnt signal intensity.

Broad genomic and transcriptional analysis reveals a highly derived genome in dinoflagellate mitochondria

Abstract: Dinoflagellates comprise an ecologically significant and diverse eukaryotic phylum that is sister to the phylum containing apicomplexan endoparasites. The mitochondrial genome of apicomplexans is uniquely reduced in gene content and size, encoding only three proteins and two ribosomal RNAs (rRNAs) within a highly compacted 6 kb DNA. Dinoflagellate mitochondrial genomes have been comparatively poorly studied: limited available data suggest some similarities with apicomplexan mitochondrial genomes but an even more radical type of genomic organization. Here, we investigate structure, content and expression of dinoflagellate mitochondrial genomes. From two dinoflagellates, Crypthecodinium cohnii and Karlodinium micrum, we generated over 42 kb of mitochondrial genomic data that indicate a reduced gene content paralleling that of mitochondrial genomes in apicomplexans, i.e., only three protein-encoding genes and at least eight conserved components of the highly fragmented large and small subunit rRNAs. Unlike in apicomplexans, dinoflagellate mitochondrial genes occur in multiple copies, often as gene fragments, and in numerous genomic contexts. Analysis of cDNAs suggests several novel aspects of dinoflagellate mitochondrial gene expression. Polycistronic transcripts were found, standard start codons are absent, and oligoadenylation occurs upstream of stop codons, resulting in the absence of termination codons. Transcripts of at least one gene, cox3, are apparently trans-spliced to generate full-length mRNAs. RNA substitutional editing, a process previously identified for mRNAs in dinoflagellate mitochondria, is also implicated in rRNA expression. The dinoflagellate mitochondrial genome shares the same gene complement and fragmentation of rRNA genes with its apicomplexan counterpart. However, it also exhibits several unique characteristics. Most notable are the expansion of gene copy numbers and their arrangements within the genome, RNA editing, loss of stop codons, and use of trans-splicing.

Ankyrin repeat domain-encoding genes in the wPip strain of Wolbachia from the Culex pipiens group

Abstract: Background: Wolbachia are obligate endosymbiotic bacteria maternally transmitted through the egg cytoplasm that are responsible for several reproductive disorders in their insect hosts, such as cytoplasmic incompatibility (CI) in infected mosquitoes. Species in the Culex pipiens complex display an unusually high number of Wolbachia-induced crossing types, and based on present data, only the wPip strain is present. Results: The sequencing of the wPip strain of Wolbachia revealed the presence of 60 ankyrin repeat domain (ANK) encoding genes and expression studies of these genes were carried out in adult mosquitoes. One of these ANK genes, pk2, is shown to be part of an operon of three prophageassociated genes with sex-specific expression, and is present in two identical copies in the genome. Another homolog of pk2 is also present that is differentially expressed in different Cx. pipiens group strains. A further two ANK genes showed sex-specific regulation in wPip-infected Cx. pipiens group adults. Conclusion: The high number, variability and differential expression of ANK genes in wPip suggest an important role in Wolbachia biology, and the gene family provides both markers and promising candidates for the study of reproductive manipulation.

Food-associated cues alter forebrain functional connectivity as assessed with immediate early gene and proenkephalin expression

Abstract: Cues predictive of food availability are powerful modulators of appetite as well as food-seeking and ingestive behaviors. The neurobiological underpinnings of these conditioned responses are not well understood. Monitoring regional immediate early gene expression is a method used to assess alterations in neuronal metabolism resulting from upstream intracellular and extracellular signaling. Furthermore, assessing the expression of multiple immediate early genes offers a window onto the possible sequelae of exposure to food cues, since the function of each gene differs. We used immediate early gene and proenkephalin expression as a means of assessing food cue-elicited regional activation and alterations in functional connectivity within the forebrain. Contextual cues associated with palatable food elicited conditioned motor activation and corticosterone release in rats. This motivational state was associated with increased transcription of the activity-regulated genes homer1a, arc, zif268, ngfi-b and c-fos in corticolimbic, thalamic and hypothalamic areas and of proenkephalin within striatal regions. Furthermore, the functional connectivity elicited by food cues, as assessed by an inter-regional multigene-expression correlation method, differed substantially from that elicited by neutral cues. Specifically, food cues increased cortical engagement of the striatum, and within the nucleus accumbens, shifted correlations away from the shell towards the core. Exposure to the food-associated context also induced correlated gene expression between corticostriatal networks and the basolateral amygdala, an area critical for learning and responding to the incentive value of sensory stimuli. This increased corticostriatal-amygdalar functional connectivity was absent in the control group exposed to innocuous cues. The results implicate correlated activity between the cortex and the striatum, especially the nucleus accumbens core and the basolateral amygdala, in the generation of a conditioned motivated state that may promote excessive food intake. The upregulation of a number of genes in unique patterns within corticostriatal, thalamic, and hypothalamic networks suggests that food cues are capable of powerfully altering neuronal processing in areas mediating the integration of emotion, cognition, arousal, and the regulation of energy balance. As many of these genes play a role in plasticity, their upregulation within these circuits may also indicate the neuroanatomic and transcriptional correlates of extinction learning.

Primate brain architecture and selection in relation to sex

Abstract: Social and competitive demands often differ between the sexes in mammals. These differing demands should be expected to produce variation in the relative sizes of various brain structures. Sexual selection on males can be predicted to influence brain components handling sensory-motor skills that are important for physical competition or neural pathways involving aggression. Conversely, because female fitness is more closely linked to ecological factors and social interactions that enable better acquisition of resources, social selection on females should select for brain components important for navigating social networks. Sexual and social selection acting on one sex could produce sexual dimorphism in brain structures, which would result in larger species averages for those same brain structures. Alternatively, sex-specific selection pressures could produce correlated effects in the other sex, resulting in larger brain structures for both males and females of a species. Data are presently unavailable for the sex-specific sizes of brain structures for anthropoid primates, but under either scenario, the effects of sexual and social selection should leave a detectable signal in average sizes of brain structures for different species. The degree of male intra-sexual selection was positively correlated with several structures involved in autonomic functions and sensory-motor skills, and in pathways relating to aggression and aggression control. The degree of male intra-sexual selection was not correlated with relative neocortex size, which instead was significantly positively correlated with female social group size, but negatively correlated with male group size. Sexual selection on males and social selection on females have exerted different effects on primate brain architecture. Species with a higher degree of male intra-sexual selection carry a neural signature of an evolutionary history centered on physical conflicts, but no traces of increased demands on sociocognitive tasks. Conversely, female sociality is indicated to have driven the evolution of socio-cognitive skills. Primate brain architecture is therefore likely to be a product of ecological and species-specific social factors as well as different sex-specific selection pressures. Our results also highlight the need for acquisition and analysis of sex-specific brain components in mammals.

Spatial chemical conservation of hot spot interactions in protein-protein complexes

Abstract: Conservation of the spatial binding organizations at the level of physico-chemical interactions is important for the formation and stability of protein-protein complexes as well as protein and drug design. Due to the lack of computational tools for recognition of spatial patterns of interactions shared by a set of protein-protein complexes, the conservation of such interactions has not been addressed previously. We performed extensive spatial comparisons of physico-chemical interactions common to different types of protein-protein complexes. We observed that 80% of these interactions correspond to known hot spots. Moreover, we show that spatially conserved interactions allow prediction of hot spots with a success rate higher than obtained by methods based on sequence or backbone similarity. Detection of spatially conserved interaction patterns was performed by our novel MAPPIS algorithm. MAPPIS performs multiple alignments of the physico-chemical interactions and the binding properties in three dimensional space. It is independent of the overall similarity in the protein sequences, folds or amino acid identities. We present examples of interactions shared between complexes of colicins with immunity proteins, serine proteases with inhibitors and T-cell receptors with superantigens. We unravel previously overlooked similarities, such as the interactions shared by the structurally different RNase-inhibitor families. The key contribution of MAPPIS is in discovering the 3D patterns of physico-chemical interactions. The detected patterns describe the conserved binding organizations that involve energetically important hot spot residues and are crucial for the protein-protein associations.