A field guide to cultivating computational biology

Evolving in sync with the computation revolution over the past 30 years, computational biology has emerged as a mature scientific field. While the field has made major contributions toward improving scientific knowledge and human health, individual computational biology practitioners at various institutions often languish in career development. As optimistic biologists passionate about the future of our field, we propose solutions for both eager and reluctant individual scientists, institutions, publishers, funding agencies, and educators to fully embrace computational biology. We believe that in order to pave the way for the next generation of discoveries, we need to improve recognition for computational biologists and better align pathways of career success with pathways of scientific progress. With 10 outlined steps, we call on all adjacent fields to move away from the traditional individual, single-discipline investigator research model and embrace multidisciplinary, data-driven, team science.

Do you want to attract computational biologists to your project or to your department? Despite the major contributions of computational biology, those attempting to bridge the interdisciplinary gap often languish in career advancement, publication, and grant review. Here, sixteen computational biologists around the globe present "A field guide to cultivating computational biology," focusing on solutions.

Biology in the digital era requires computation and collaboration. A modern research project may include multiple model systems, use multiple assay technologies, collect varying data types, and require complex computational strategies, which together make effective design and execution difficult or impossible for any individual scientist. While some labs, institutions, funding bodies, publishers, and other educators have already embraced a team science model in computational biology and thrived [1–7], others who have not yet fully adopted it risk severely lagging behind the cutting edge. We propose a general solution: “deep integration” between biology and the computational sciences. Many different collaborative models can yield deep integration, and different problems require different approaches (Fig 1).Open in a separate windowFig 1Supporting interdisciplinary team science will accelerate biological discoveries.Scientists who have little exposure to different fields build silos, in which they perform science without external input. To solve hard problems and to extend your impact, collaborate with diverse scientists, communicate effectively, recognize the importance of core facilities, and embrace research parasitism. In biologically focused parasitism, wet lab biologists use existing computational tools to solve problems; in computationally focused parasitism, primarily dry lab biologists analyze publicly available data. Both strategies maximize the use and societal benefit of scientific data.In this article, we define computational science extremely broadly to include all quantitative approaches such as computer science, statistics, machine learning, and mathematics. We also define biology broadly, including any scientific inquiry pertaining to life and its many complications. A harmonious deep integration between biology and computer science requires action—we outline 10 immediate calls to action in this article and aim our speech directly at individual scientists, institutions, funding agencies, and publishers in an attempt to shift perspectives and enable action toward accepting and embracing computational biology as a mature, necessary, and inevitable discipline (Box 1).Box 1. Ten calls to action for individual scientists, funding bodies, publishers, and institutions to cultivate computational biology. Many actions require increased funding support, while others require a perspective shift. For those actions that require funding, we believe convincing the community of need is the first step toward agencies and systems allocating sufficient support

1. Respect collaborators’ specific research interests and motivationsProblem: Researchers face conflicts when their goals do not align with collaborators. For example, projects with routine analyses provide little benefit for computational biologists.Solution: Explicit discussion about interests/expertise/goals at project onset.Opportunity: Clearly defined expectations identify gaps, provide commitment to mutual benefit.
2. Seek necessary input during project design and throughout the project life cycleProblem: Modern research projects require multiple experts spanning the project’s complexity.Solution: Engage complementary scientists with necessary expertise throughout the entire project life cycle.Opportunity: Better designed and controlled studies with higher likelihood for success.
3. Provide and preserve budgets for computational biologists’ workProblem: The perception that analysis is “free” leads to collaborator budget cuts.Solution: When budget cuts are necessary, ensure that they are spread evenly.Opportunity: More accurate, reproducible, and trustworthy computational analyses.
4. Downplay publication author order as an evaluation metric for computational biologistsProblem: Computational biologist roles on publications are poorly understood and undervalued.Solution: Journals provide more equitable opportunities, funding bodies and institutions improve understanding of the importance of team science, scientists educate each other.Opportunity: Engage more computational biologist collaborators, provide opportunities for more high-impact work.
5. Value software as an academic productProblem: Software is relatively undervalued and can end up poorly maintained and supported, wasting the time put into its creation.Solution: Scientists cite software, and funding bodies provide more software funding opportunities.Opportunity: More high-quality maintainable biology software will save time, reduce reimplementation, and increase analysis reproducibility.
6. Establish academic structures and review panels that specifically reward team scienceProblem: Current mechanisms do not consistently reward multidisciplinary work.Solution: Separate evaluation structures to better align peer review to reward indicators of team science.Opportunity: More collaboration to attack complex multidisciplinary problems.
7. Develop and reward cross-disciplinary training and mentoringProblem: Academic labs and institutions are often insufficiently equipped to provide training to tackle the next generation of biological problems, which require computational skills.Solution: Create better training programs aligned to necessary on-the-job skills with an emphasis on communication, encourage wet/dry co-mentorship, and engage younger students to pursue computational biology.Opportunity: Interdisciplinary students uncover important insights in their own data.
8. Support computing and experimental infrastructure to empower computational biologistsProblem: Individual computational labs often fund suboptimal cluster computing systems and lack access to data generation facilities.Solution: Institutions can support centralized compute and engage core facilities to provide data services.Opportunity: Time and cost savings for often overlooked administrative tasks.
9. Provide incentives and mechanisms to share open data to empower discovery through reanalysisProblem: Data are often siloed and have untapped potential.Solution: Provide institutional data storage with standardized identifiers and provide separate funding mechanisms and publishing venues for data reuse.Opportunity: Foster new breed of researchers, “research parasites,” who will integrate multimodal data and enhance mechanistic insights.
10. Consider infrastructural, ethical, and cultural barriers to clinical data accessProblem: Identifiable health data, which include sensitive information that must be kept hidden, are distributed and disorganized, and thus underutilized.Solution: Leadership must enforce policies to share deidentifiable data with interoperable metadata identifiers.Opportunity: Derive new insights from multimodal data integration and build datasets with increased power to make biological discoveries.

     

Evaluating the Ribosomal Internal Transcribed Spacer (ITS) as a Candidate Dinoflagellate Barcode Marker

Background

DNA barcoding offers an efficient way to determine species identification and to measure biodiversity. For dinoflagellates, an ancient alveolate group of about 2000 described extant species, DNA barcoding studies have revealed large amounts of unrecognized species diversity, most of which is not represented in culture collections. To date, two mitochondrial gene markers, Cytochrome Oxidase I (COI) and Cytochrome b oxidase (COB), have been used to assess DNA barcoding in dinoflagellates, and both failed to amplify all taxa and suffered from low resolution. Nevertheless, both genes yielded many examples of morphospecies showing cryptic speciation and morphologically distinct named species being genetically similar, highlighting the need for a common marker. For example, a large number of cultured Symbiodinium strains have neither taxonomic identification, nor a common measure of diversity that can be used to compare this genus to other dinoflagellates.

Methodology/Principal Findings

The purpose of this study was to evaluate the Internal Transcribed Spacer units 1 and 2 (ITS) of the rDNA operon, as a high resolution marker for distinguishing species dinoflagellates in culture. In our study, from 78 different species, the ITS barcode clearly differentiated species from genera and could identify 96% of strains to a known species or sub-genus grouping. 8.3% showed evidence of being cryptic species. A quarter of strains identified had no previous species identification. The greatest levels of hidden biodiversity came from Scrippsiella and the Pfiesteriaceae family, whilst Heterocapsa strains showed a high level of mismatch to their given species name.

Conclusions/Significance

The ITS marker was successful in confirming species, revealing hidden diversity in culture collections. This marker, however, may have limited use for environmental barcoding due to paralogues, the potential for unidentifiable chimaeras and priming across taxa. In these cases ITS would serve well in combination with other markers or for specific taxon studies.     

Immunochemical studies on blood groups: The internal structure and immunological properties of water-soluble human blood group A substance studied by Smith degradation,liberation, and fractionation of oligosaccharides and reaction with lectins

Carbohydrate structures in the interior of a blood group A active substance (MSS) were exposed by one and by two Smith degradations. Reactivities of the original glycoprotein and its Smith degraded products with 13 different lectins and with anti-I Ma were studied by quantitative precipitin assay. MSS and its first Smith degraded product completely precipitated Ricinus communis hemagglutinin with five times less of the first Smith degraded glycoprotein being required for 50% precipitation. The second Smith degraded material precipitated only 90% of the lectin. MSS did not precipitate peanut lectin, whereas its first and second Smith degraded products completely precipitated the lectin. The first Smith degraded glycoprotein also reacted well with Wistaria floribunda, Maclura pomifera, Bauhinia purpurea alba, and Geodia lectins indicating that its carbohydrate moiety could contain dGalNAc, dGalβ1 → 3dGalNAc, dGalβ1 → 4dGlcNAc, dGalβ1 → 3dGlcNAcβ1 → 3dGal and/or dGalβ1 → 4dGlcNAcβ1 → 6dGal and/or dGalβ1 → 4dGlcNAcβ1 → 6dGalNAc determinants at nonreducing ends. The second Smith degraded material precipitated well with Ricinus communis hemagglutinin, Arachis hypogaea, Geodia cydonium, Maclura pomifera, and Helix pomatia lectins showing that dGalNAc, dGalβ1 → 3dGalNAc, dGalβ1 → 4dGlcNAc residues at terminal nonreducing ends could be involved. Monoclonal anti-I Ma (group 1) serum reacted strongly with the first Smith degraded product indicating large numbers of anti-I Ma determinants, dGalβ1 → 4dGlcNAcβ1 → d 6dGal and/or dGalβ1 → 4dGlcNAcβ1 → 6dGalNAc at nonreducing ends. The comparable activities of the native and Smith degraded products with wheat germ lectin indicate capacity to react with DGlcNAc residues at nonreducing ends and/or at positions in the interior of the chain. The totality of lectin reactivities indicates heterogeneity of the carbohydrate side chains. Oligosaccharides with ³H at their reducing ends released from the protein core of the first and second Smith degraded products were obtained by treatment with 0.05 m NaOH and 1 M NaB³H₄ at 50 °C for 16 h (Carlson degradation). The liberated reduced oligosaccharides were fractionated by dialysis, followed by retardion, Bio-Gel P-2, P-4, and P-6 columns. They were further purified on charcoal-celite columns, and by preparative paper chromatography and high-pressure liquid chromatography. Their distribution by size was estimated by the yields on dialysis, Bio-Gel P-2, and Bio-Gel P-6 chromatography, and from the radioactivity of the reduced sugars. Of the oligosaccharide fractions from the first Smith degraded product, about 77% of the carbohydrate side chain residues contained from 1 to 6 sugars, 13% from 7 to perhaps 12 sugars, and 10% was nondialyzable (polysaccharides and glycopeptide fragments). Of the second Smith degraded product, approximately 82% of carbohydrate residues had from 1 to 6 sugars, 14% from 7 to perhaps 20 sugars and 4% was nondialyzable. The biological activity profile of the two Smith degraded products together with the size distributions of the oligosaccharides indicated that their carbohydrate side chains, comprised a heterogeneous population ranging in size from 1 to about 12 sugars. When most of these chains that are shorter than hexasaccharides are fully characterized it may be possible to reconstruct the overall structure of the carbohydrate moiety of the blood group substances and account for their biological activities.     

Seed Germination Biology of the Narrowly Endemic Species Lesquerella stonensis (Brassicaceae)

Abstract Lesquerella stonensis (Brassicaceae) is an obligate winter annual endemic to a small portion of Rutherford County in the Central Basin of Tennessee, where it grows in disturbed habitats. This species forms a persistent seed bank, and seeds remain viable in the soil for at least 6 years. Seeds are dormant at maturity in May and are dispersed as soon as they ripen. Some of the seeds produced in the current year, as well as some of those in the persistent seed bank, afterripen during late spring and summer; others do not afterripen and thus remain dormant. Seeds require actual or simulated spring/summer temperatures to come out of dormancy. Germination occurs in September and October. Fully afterripened seeds germinate over a wide range of thermoperiods (15/6–35/20°C) and to a much higher percentage in light (14 h photoperiod) than in darkness. The optimum daily thermoperiod for germination was 30/15°C. Nondormant seeds that do not germinate in autumn are induced back into dormancy (secondary dormancy) by low temperatures (e.g., 5°C) during winter, and those that are dormant do not afterripen; thus seeds cannot germinate in spring. These seed dormancy/ germination characteristics of L. stonensis do not differ from those reported for some geographically widespread, weedy species of winter annuals and thus do not help account for the narrow endemism of this species.     

Skeletal adaptations and phylogeny of the oldest mole Eotalpa (Talpidae,Lipotyphla, Mammalia) from the UK Eocene: the beginning of fossoriality in moles

The oldest talpid, Eotalpa, was previously known only from isolated cheek teeth from the European late Middle Eocene to earliest Oligocene. Screenwashing of Late Eocene sediments of the Hampshire Basin, UK, has yielded cranial and postcranial elements: maxilla, dentary, ulna, metacarpals, distal tibia, astragalus, calcaneum, metatarsals and phalanges. In addition to M_1–2 myotodonty, typical talpid features are as follows: ulna with long medially curved olecranon and deep abductor fossa and astragalar body with lateral process. However, Eotalpa retains certain soricid‐like primitive states (M¹ preparacrista, P⁴ with prominent mesiolingual protocone lobe, strongly angled astragalar neck and calcaneum with no space for a cuboid medial process) not found in modern talpids. Eotalpa is more derived than the most primitive living talpid Uropsilus in having lost the M^1–2 talon shelf, developed a convex radial facet on the ulna, an incipient proximal olecranon crest, relatively shorter metapodials and depressed manual unguals. Its astragalus with medial trochlear ridge taller than the lateral one and massive medial plantar process is typical of the Lipotyphla. Eotalpa lacks synostosis of tibia and fibula, found in other Talpidae, Soricidae and Erinaceidae, suggesting that synostosis in these groups has been independently acquired. Cladistic analysis places Eotalpa as stem member of the Talpidae and shows that much homoplasy arose during the early evolution of the family. Ground dwelling in Eotalpa is indicated by the following: astragalus with a medially dipping head, curved in a single plane; calcaneum with distal peroneal process and strongly overlapping ectal and sustentacular facets; and matching sized ectal and sustentacular facets on calcaneum and astragalus. These features would have restricted ankle mobility. Ungual and metatarsal shape and ulnar structure suggest a primitive stage in fossorial evolution and argue against a semiaquatic precursor stage in talpid fossoriality. Shrew‐moles may represent a reversal to surface foraging rather than an intermediate stage in fossoriality.     

Calcium mobilization and muscle contraction induced by acetylcholine in swine trachealis

The roles of Ca²⁺ mobilization in development of tension induced by acetylcholine (ACh, 0.1–100 µM) in swine tracheal smooth muscle strips were studied. Under control conditions, ACh induced a transient increase in free cytosolic calcium concentration ([Ca²⁺]_i) that declined to a steady-state level. The peak increase in [Ca²⁺]_i correlated with the magnitude of tension at each [ACh] after a single exposure to ACh, while the steady-state [Ca²⁺]_i did not. Removal of extracellular Ca²⁺ had little effect on peak [Ca²⁺]_i but greatly reduced steady-state increases in [Ca²⁺]_i and tension. Verapamil inhibited steady-state [Ca²⁺]_i only at [ACh]<1 µM. After depletion of internal Ca²⁺ stores by 10 min exposure to ACh in Ca²⁺-free solution and then washout of ACh for 5 min in Ca²⁺-free solution, simultaneous re-exposure to ACh in the presence of 2.5 mM Ca²⁺ increased [Ca²⁺]_i to the control steady-state level without overshoot. The tension attained was the same as control for each [ACh] used. Continuous exposure to successively increasing [ACh] (0.1–100 µM) also reduced the overshoot of [Ca²⁺]_i at 10 and 100 µM ACh, yet tension reached control levels at each [ACh] used. We conclude that the steady-state increase in [Ca²⁺]_i is necessary for tension maintenance and is dependent on Ca²⁺ influx through voltage-gated calcium channels at 0.1 µM ACh and through a verapamil-insensitive pathway at 10 and 100 µM. The initial transient increase in calcium arises from intracellular stores and is correlated with the magnitude of tension only in muscles that have completely recovered from previous exposure to agonists.     

GERMINATION ECOPHYSIOLOGY OF HYDROPHYLLUM APPENDICULATUM,A MESIC FOREST BIENNIAL

In north central Kentucky, seeds of the mesic forest biennial Hydrophyllum appendiculatum Michx., are innately dormant at maturity in June. Under natural and simulated seasonal temperature changes, dormancy break occurred in two stages. Root dormancy was broken by high summer temperatures, and shoot dormancy was broken by low winter temperatures. Consequently, roots emerged from seeds during autumn, and cotyledons emerged the following spring. A 90-day warm (30/15 C) stratification treatment broke root dormancy, but the roots emerged only after transfer to lower temperatures. After the warm stratification treatment, roots emerged from 93, 73, 6 and 9% of the seeds incubated at 5, 15/6, 20/10 and 30/15 C (12/12 hr), respectively. Zero, 28, 56 and 84 days of cold (5 C) stratification of seeds with emerged roots resulted in 9, 21, 49 and 82% cotyledon emergence, respectively, at 20/10 C. Thus, H. appendiculatum exhibits a type of morpho-physiological dormancy known as epicotyl dormancy. Although many seeds germinate the first year, others remain dormant and germinate in successive years until the fourth season after ripening.     

A simple,rapid assay for cysteamine and other thiols

Cysteamine is under investigation as an aid in radiation therapy and as a treatment for the inherited disorder cystinosis. An assay is presented for its measurement in biological fluids. The specific reaction of thiosulfonates with sulfhydryl compounds is employed to form a radiolabeled derivative of cysteamine which is then isolated by high-voltage electrophoresis on paper. Cysteamine can be measured in aqueous solutions, plasma, and urine with this method.     

A second generation of double mutant cholera toxin adjuvants: enhanced immunity without intracellular trafficking

Nasal application of native cholera toxin (nCT) as a mucosal adjuvant has potential toxicity for the CNS through binding to GM1 gangliosides in the olfactory nerves. Although mutants of cholera toxin (mCTs) have been developed that show mucosal adjuvant activity without toxicity, it still remains unclear whether these mCTs will induce CNS damage. To help overcome these concerns, in this study we created new double mutant CTs (dmCTs) that have two amino acid substitutions in the ADP-ribosyltransferase active center (E112K) and COOH-terminal KDEL (E112K/KDEV or E112K/KDGL). Confocal microscopic analysis showed that intracellular localization of dmCTs differed from that of mCTs and nCTs in intestinal epithelial T84 cells. Furthermore, both dmCTs exhibited very low toxicity in the Y1 cell assay and mouse ileal loop tests. When mucosal adjuvanticity was examined, both dmCTs induced enhanced OVA-specific immune responses in both mucosal and systemic lymphoid tissues. Interestingly, although both dmCT E112K/KDEV and dmCT E112K/KDGL showed high Th2-type and significant Th1-type cytokine responses by OVA-specific CD4+ T cells, dmCT E112K/KDEV exhibited significantly lower Th1-type cytokine responses than did nCT and dmCT E112K/KDGL. These results show that newly developed dmCTs retain strong biological adjuvant activity without CNS toxicity.