Ontogeny tends to recapitulate phylogeny in digital organisms

Abstract Biologists have long debated whether ontogeny recapitulates phylogeny and, if so, why. Two plausible explanations are that (i) changes to early developmental stages are selected against because they tend to disrupt later development and (ii) simpler structures often precede more complex ones in both ontogeny and phylogeny if the former serve as building blocks for the latter. It is difficult to test these hypotheses experimentally in natural systems, so we used a computational system that exhibits evolutionary dynamics. We observed that ontogeny does indeed recapitulate phylogeny; traits that arose earlier in a lineage's history also tended to be expressed earlier in the development of individuals. The relative complexity of traits contributed substantially to this correlation, but a significant tendency toward recapitulation remained even after accounting for trait complexity. This additional effect provides evidence that selection against developmental disruption also contributed to the conservation of early stages in development.     

Pollen flow and post-pollination barriers in two varieties of Dactylorhiza incarnata s.l. (Orchidaceae)

Dactylorhiza orchids are known for their high variation in morphology and distinct varieties have been named in D. incarnata s.l. However, it is not known how these varieties interact in mixed populations and why they remain stable. We conducted three field experiments in West-Estonian populations of D. incarnata to examine if the two most common varieties co-occurring these are separated from each other either by pre-pollination or post-pollination reproductive barriers. We found that pollinia were far more frequently transferred between the purple-flowered var. incarnata and the pale-flowered var. ochroleuca than between plants of the same variety. Furthermore, in hand-pollination and germination experiments concerning the same two varieties, we found that pollen source (self-pollination, within- or between-varieties pollination) did not affect seed production or probability of fungus infection of the germinating seeds. These two varieties of D. incarnata thus had no pre-pollination or early functioning post-pollination reproductive barriers. Post-pollination barriers may, however, act later in seedling or adult stage.     

Armed and accurate: engineering cytotoxic T cells for eradication of leukemia

After publication we discovered an error in the identification of the origin of the cell line reported in our article in BMC Biotechnology (2010, 10:50), entitled "Ao38, a new cell line from eggs of the black witch moth, Ascalapha odorata (Lepidoptera: Noctuidae), is permissive for AcMNPV infection and produces high levels of recombinant proteins". Upon analysis of primary A. odorata cultures, we found that they were contaminated with cells of Trichoplusia ni origin. The origin of the Ao38 cell line was determined as T. ni using three marker genes and the Ao38 cell line was renamed BTI-Tnao38. References to the origin of the cell line as Ascalapha odorata should be replaced with "a cell line of Trichoplusia ni origin". The absence of TNCL virus detection in the BTI-Tnao38 (Ao38) cell line was confirmed using a highly sensitive RT-PCR protocol capable of detecting TNCL virus RNA at approximately 0.018 copies/cell. Because of these observations, we have revised the title of the original article to "Correction: BTI-Tnao38, a new cell line derived from Trichoplusia ni, is permissive for AcMNPV infection and produces high levels of recombinant proteins" and two additional authors were added to reflect their contributions to the analysis of this cell line.     

The Quantitative Methods Boot Camp: Teaching Quantitative Thinking and Computing Skills to Graduate Students in the Life Sciences

The past decade has seen a rapid increase in the ability of biologists to collect large amounts of data. It is therefore vital that research biologists acquire the necessary skills during their training to visualize, analyze, and interpret such data. To begin to meet this need, we have developed a “boot camp” in quantitative methods for biology graduate students at Harvard Medical School. The goal of this short, intensive course is to enable students to use computational tools to visualize and analyze data, to strengthen their computational thinking skills, and to simulate and thus extend their intuition about the behavior of complex biological systems. The boot camp teaches basic programming using biological examples from statistics, image processing, and data analysis. This integrative approach to teaching programming and quantitative reasoning motivates students’ engagement by demonstrating the relevance of these skills to their work in life science laboratories. Students also have the opportunity to analyze their own data or explore a topic of interest in more detail. The class is taught with a mixture of short lectures, Socratic discussion, and in-class exercises. Students spend approximately 40% of their class time working through both short and long problems. A high instructor-to-student ratio allows students to get assistance or additional challenges when needed, thus enhancing the experience for students at all levels of mastery. Data collected from end-of-course surveys from the last five offerings of the course (between 2012 and 2014) show that students report high learning gains and feel that the course prepares them for solving quantitative and computational problems they will encounter in their research. We outline our course here which, together with the course materials freely available online under a Creative Commons License, should help to facilitate similar efforts by others.This is part of the PLOS Computational Biology Education collection.     

Field Metabolic Rate and PCB Adipose Tissue Deposition Efficiency in East Greenland Polar Bears Derived from Contaminant Monitoring Data

Climate change will increasingly affect the natural habitat and diet of polar bears (Ursus maritimus). Understanding the energetic needs of polar bears is therefore important. We developed a theoretical method for estimating polar bear food consumption based on using the highly recalcitrant polychlorinated biphenyl (PCB) congener, 2,2′,4,4′,55-hexaCB (CB153) in bear adipose tissue as an indicator of food intake. By comparing the CB153 tissue concentrations in wild polar bears with estimates from a purposely designed individual-based model, we identified the possible combinations of field metabolic rates (FMR) and CB153 deposition efficiencies in East Greenland polar bears. Our simulations indicate that if 30% of the CB153 consumed by polar bear individuals were deposited into their adipose tissue, the corresponding FMR would be only two times the basal metabolic rate. In contrast, if the modelled CB153 deposition efficiency were 10%, adult polar bears would require six times more energy than that needed to cover basal metabolism. This is considerably higher than what has been assumed for polar bears in previous studies though it is similar to FMRs found in other marine mammals. An implication of this result is that even relatively small reductions in future feeding opportunities could impact the survival of East Greenland polar bears.     

Mitochondrial DNA Variant in COX1 Subunit Significantly Alters Energy Metabolism of Geographically Divergent Wild Isolates in Caenorhabditis elegans

Mitochondrial DNA (mtDNA) sequence variation can influence the penetrance of complex diseases and climatic adaptation. While studies in geographically defined human populations suggest that mtDNA mutations become fixed when they have conferred metabolic capabilities optimally suited for a specific environment, it has been challenging to definitively assign adaptive functions to specific mtDNA sequence variants in mammals. We investigated whether mtDNA genome variation functionally influences Caenorhabditis elegans wild isolates of distinct mtDNA lineages and geographic origins. We found that, relative to N2 (England) wild-type nematodes, CB4856 wild isolates from a warmer native climate (Hawaii) had a unique p.A12S amino acid substitution in the mtDNA-encoded COX1 core catalytic subunit of mitochondrial complex IV (CIV). Relative to N2, CB4856 worms grown at 20 °C had significantly increased CIV enzyme activity, mitochondrial matrix oxidant burden, and sensitivity to oxidative stress but had significantly reduced lifespan and mitochondrial membrane potential. Interestingly, mitochondrial membrane potential was significantly increased in CB4856 grown at its native temperature of 25 °C. A transmitochondrial cybrid worm strain, chpIR (M, CB4856 > N2), was bred as homoplasmic for the CB4856 mtDNA genome in the N2 nuclear background. The cybrid strain also displayed significantly increased CIV activity, demonstrating that this difference results from the mtDNA-encoded p.A12S variant. However, chpIR (M, CB4856 > N2) worms had significantly reduced median and maximal lifespan relative to CB4856, which may relate to their nuclear–mtDNA genome mismatch. Overall, these data suggest that C. elegans wild isolates of varying geographic origins may adapt to environmental challenges through mtDNA variation to modulate critical aspects of mitochondrial energy metabolism.     

Individual small in‐stream barriers contribute little to strong local population genetic structure five strictly aquatic macroinvertebrate taxa

Water flow in river networks is frequently regulated by man‐made in‐stream barriers. These obstacles can hinder dispersal of aquatic organisms and isolate populations leading to the loss of genetic diversity. Although millions of small in‐stream barriers exist worldwide, their impact on dispersal of macroinvertebrates remains unclear. Therefore, we, therefore, assessed the effects of such barriers on the population structure and effective dispersal of five macroinvertebrate species with strictly aquatic life cycles: the amphipod crustacean Gammarus fossarum (clade 11), three snail species of the Ancylus fluviatilis species complex and the flatworm Dugesia gonocephala. We studied populations at nine weirs and eight culverts (3 pipes, 5 tunnels), built 33–109 years ago, mainly in the heavily fragmented catchment of the river Ruhr (Sauerland, Germany). To assess fragmentation and barrier effects, we generated genome‐wide SNP data using ddRAD sequencing and evaluated clustering, differentiation between populations up‐ and downstream of each barrier and effective migration rates among sites and across barriers. Additionally, we applied population genomic simulations to assess expected differentiation patterns under different gene flow scenarios. Our data show that populations of all species are highly isolated at regional and local scales within few kilometers. While the regional population structure likely results from historical processes, the strong local differentiation suggests that contemporary dispersal barriers exist. However, we identified significant barrier effects only for pipes (for A. fluviatilis II and III) and few larger weirs (>1.3 m; for D. gonocephala). Therefore, our data suggest that most small in‐stream barriers can probably be overcome by all studied taxa frequently enough to prevent fragmentation. However, it remains to be tested if the strong local differentiation is a result of a cumulative effect of small barriers, or if larger in‐stream barriers, land use, chemical pollution, urbanization, or a combination of these factors impede gene flow.     

Compound design guidelines for evading the efflux and permeation barriers of Escherichia coli with the oxazolidinone class of antibacterials: Test case for a general approach to improving whole cell Gram-negative activity

Previously we reported the results from an effort to improve Gram-negative antibacterial activity in the oxazolidinone class of antibiotics via a systematic medicinal chemistry campaign focused entirely on C-ring modifications. In that series we set about testing if the efflux and permeation barriers intrinsic to the outer membrane of Escherichia coli could be rationally overcome by designing analogs to reside in specific property limits associated with Gram-negative activity: i) low MW (<400), ii) high polarity (clogD_7.4 <1), and iii) zwitterionic character at pH 7.4. Indeed, we observed that only analogs residing within these limits were able to overcome these barriers. Herein we report the results from a parallel effort where we explored structural changes throughout all three rings in the scaffold for the same purpose. Compounds were tested against a diagnostic MIC panel of Escherichia coli and Staphylococcus aureus strains to determine the impact of combining structural modifications in overcoming the OM barriers and in bridging the potency gap between the species. The results demonstrated that distributing the charge-carrying moieties across two rings was also beneficial for avoidance of the outer membrane barriers. Importantly, analysis of the structure-permeation relationship (SPR) obtained from this and the prior study indicated that in addition to MW, polarity, and zwitterionic character, having ≤4 rotatable bonds is also associated with evasion of the OM barriers. These combined results provide the medicinal chemist with a framework and strategy for overcoming the OM barriers in GNB in antibacterial drug discovery efforts.     

Genetic Barriers to Historical Gene Flow between Cryptic Species of Alpine Bumblebees Revealed by Comparative Population Genomics

Evidence is accumulating that gene flow commonly occurs between recently diverged species, despite the existence of barriers to gene flow in their genomes. However, we still know little about what regions of the genome become barriers to gene flow and how such barriers form. Here, we compare genetic differentiation across the genomes of bumblebee species living in sympatry and allopatry to reveal the potential impact of gene flow during species divergence and uncover genetic barrier loci. We first compared the genomes of the alpine bumblebee Bombus sylvicola and a previously unidentified sister species living in sympatry in the Rocky Mountains, revealing prominent islands of elevated genetic divergence in the genome that colocalize with centromeres and regions of low recombination. This same pattern is observed between the genomes of another pair of closely related species living in allopatry (B. bifarius and B. vancouverensis). Strikingly however, the genomic islands exhibit significantly elevated absolute divergence (d_XY) in the sympatric, but not the allopatric, comparison indicating that they contain loci that have acted as barriers to historical gene flow in sympatry. Our results suggest that intrinsic barriers to gene flow between species may often accumulate in regions of low recombination and near centromeres through processes such as genetic hitchhiking, and that divergence in these regions is accentuated in the presence of gene flow.     

An Affinity–Effect Relationship for Microbial Communities in Plant–Soil Feedback Loops

Feedback loops involving soil microorganisms can regulate plant populations. Here, we hypothesize that microorganisms are most likely to play a role in plant–soil feedback loops when they possess an affinity for a particular plant and the capacity to consistently affect the growth of that plant for good or ill. We characterized microbial communities using whole-community DNA fingerprinting from multiple "home-and-away" experiments involving giant ragweed (Ambrosia trifida L.) and common sunflower (Helianthus annuus L.), and we looked for affinity–effect relationships in these microbial communities. Using canonical ordination and partial least squares regression, we developed indices expressing each microorganism's affinity for ragweed or sunflower and its putative effect on plant biomass, and we used linear regression to analyze the relationship between microbial affinity and effect. Significant linear affinity–effect relationships were found in 75 % of cases. Affinity–effect relationships were stronger for ragweed than for sunflower, and ragweed affinity–effect relationships showed consistent potential for negative feedback loops. The ragweed feedback relationships indicated the potential involvement of multiple microbial taxa, resulting in strong, consistent affinity–effect relationships in spite of large-scale microbial variability between trials. In contrast, sunflower plant–soil feedback may involve just a few key players, making it more sensitive to underlying microbial variation. We propose that affinity–effect relationship can be used to determine key microbial players in plant–soil feedback against a low "signal-to-noise" background of complex microbial datasets.     

Interlinked nonlinear subnetworks underlie the formation of robust cellular patterns in Arabidopsis epidermis: a dynamic spatial model

Background

Analysis of microarray data has been used for the inference of gene-gene interactions. If, however, the aim is the discovery of disease-related biological mechanisms, then the criterion for defining such interactions must be specifically linked to disease.

Results

Here we present a computational methodology that jointly analyzes two sets of microarray data, one in the presence and one in the absence of a disease, identifying gene pairs whose correlation with disease is due to cooperative, rather than independent, contributions of genes, using the recently developed information theoretic measure of synergy. High levels of synergy in gene pairs indicates possible membership of the two genes in a shared pathway and leads to a graphical representation of inferred gene-gene interactions associated with disease, in the form of a "synergy network." We apply this technique on a set of publicly available prostate cancer expression data and successfully validate our results, confirming that they cannot be due to pure chance and providing a biological explanation for gene pairs with exceptionally high synergy.

Conclusion

Thus, synergy networks provide a computational methodology helpful for deriving "disease interactomes" from biological data. When coupled with additional biological knowledge, they can also be helpful for deciphering biological mechanisms responsible for disease.     

Teaching Bioinformatics in Concert

Can biology students without programming skills solve problems that require computational solutions? They can if they learn to cooperate effectively with computer science students. The goal of the in-concert teaching approach is to introduce biology students to computational thinking by engaging them in collaborative projects structured around the software development process. Our approach emphasizes development of interdisciplinary communication and collaboration skills for both life science and computer science students.     

Adaptation in Toxic Environments: Arsenic Genomic Islands in the Bacterial Genus Thiomonas

Acid mine drainage (AMD) is a highly toxic environment for most living organisms due to the presence of many lethal elements including arsenic (As). Thiomonas (Tm.) bacteria are found ubiquitously in AMD and can withstand these extreme conditions, in part because they are able to oxidize arsenite. In order to further improve our knowledge concerning the adaptive capacities of these bacteria, we sequenced and assembled the genome of six isolates derived from the Carnoulès AMD, and compared them to the genomes of Tm. arsenitoxydans 3As (isolated from the same site) and Tm. intermedia K12 (isolated from a sewage pipe). A detailed analysis of the Tm. sp. CB2 genome revealed various rearrangements had occurred in comparison to what was observed in 3As and K12 and over 20 genomic islands (GEIs) were found in each of these three genomes. We performed a detailed comparison of the two arsenic-related islands found in CB2, carrying the genes required for arsenite oxidation and As resistance, with those found in K12, 3As, and five other Thiomonas strains also isolated from Carnoulès (CB1, CB3, CB6, ACO3 and ACO7). Our results suggest that these arsenic-related islands have evolved differentially in these closely related Thiomonas strains, leading to divergent capacities to survive in As rich environments.     

The dynamics of Brazilian protozoology over the past century

Brazilian scientists have been contributing to the protozoology field for more than 100 years with important discoveries of new species such asTrypanosoma cruzi and Leishmania spp. In this work, we used a Brazilian thesis database (Coordination for the Improvement of Higher Education Personnel) covering the period from 1987-2011 to identify researchers who contributed substantially to protozoology. We selected 248 advisors by filtering to obtain researchers who supervised at least 10 theses. Based on a computational analysis of the thesis databases, we found students who were supervised by these scientists. A computational procedure was developed to determine the advisors’ scientific ancestors using the Lattes Platform. These analyses provided a list of 1,997 researchers who were inspected through Lattes CV examination and allowed the identification of the pioneers of Brazilian protozoology. Moreover, we investigated the areas in which researchers who earned PhDs in protozoology are now working. We found that 68.4% of them are still in protozoology, while 16.7% have migrated to other fields. We observed that support for protozoology by national or international agencies is clearly correlated with the increase of scientists in the field. Finally, we described the academic genealogy of Brazilian protozoology by formalising the “forest” of Brazilian scientists involved in the study of protozoa and their vectors over the past century.