TREEasy: An automated workflow to infer gene trees,species trees,and phylogenetic networks from multilocus data

Multilocus genomic data sets can be used to infer a rich set of information about the evolutionary history of a lineage, including gene trees, species trees, and phylogenetic networks. However, user‐friendly tools to run such integrated analyses are lacking, and workflows often require tedious reformatting and handling time to shepherd data through a series of individual programs. Here, we present a tool written in Python—TREEasy—that performs automated sequence alignment (with MAFFT), gene tree inference (with IQ‐Tree), species inference from concatenated data (with IQ‐Tree and RaxML‐NG), species tree inference from gene trees (with ASTRAL, MP‐EST, and STELLS2), and phylogenetic network inference (with SNaQ and PhyloNet). The tool only requires FASTA files and nine parameters as inputs. The tool can be run as command line or through a Graphical User Interface (GUI). As examples, we reproduced a recent analysis of staghorn coral evolution, and performed a new analysis on the evolution of the “WGD clade” of yeast. The latter revealed novel patterns that were not identified by previous analyses. TREEasy represents a reliable and simple tool to accelerate research in systematic biology ( https://github.com/MaoYafei/TREEasy ).     

Genomic and phenomic analysis of island ant community assembly

Island biodiversity has long fascinated biologists as it typically presents tractable systems for unpicking the eco‐evolutionary processes driving community assembly. In general, two recurring themes are of central theoretical interest. First, immigration, diversification, and extinction typically depend on island geographical properties (e.g., area, isolation, and age). Second, predictable ecological and evolutionary trajectories readily occur after colonization, such as the evolution of adaptive trait syndromes, trends toward specialization, adaptive radiation, and eventual ecological decline. Hypotheses such as the taxon cycle draw on several of these themes to posit particular constraints on colonization and subsequent eco‐evolutionary dynamics. However, it has been challenging to examine these integrated dynamics with traditional methods. Here, we combine phylogenomics, population genomics and phenomics, to unravel community assembly dynamics among Pheidole (Hymenoptera, Formicidae) ants in the isolated Fijian archipelago. We uphold basic island biogeographic predictions that isolated islands accumulate diversity primarily through in situ evolution rather than dispersal, and population genomic support for taxon cycle predictions that endemic species have decreased dispersal ability and demography relative to regionally widespread taxa. However, rather than trending toward island syndromes, ecomorphological diversification in Fiji was intense, filling much of the genus‐level global morphospace. Furthermore, while most endemic species exhibit demographic decline and reduced dispersal, we show that the archipelago is not an evolutionary dead‐end. Rather, several endemic species show signatures of population and range expansion, including a successful colonization to the Cook islands. These results shed light on the processes shaping island biotas and refine our understanding of island biogeographic theory.     

The influence of experimentally induced polyploidy on the relationships between endopolyploidy and plant function in Arabidopsis thaliana

Whole genome duplication, leading to polyploidy and endopolyploidy, occurs in all domains and kingdoms and is especially prevalent in vascular plants. Both polyploidy and endopolyploidy increase cell size, but it is unclear whether both processes have similar effects on plant morphology and function, or whether polyploidy influences the magnitude of endopolyploidy. To address these gaps in knowledge, fifty‐five geographically separated diploid accessions of Arabidopsis thaliana that span a gradient of endopolyploidy were experimentally manipulated to induce polyploidy. Both the diploids and artificially induced tetraploids were grown in a common greenhouse environment and evaluated with respect to nine reproductive and vegetative characteristics. Induced polyploidy decreased leaf endopolyploidy and stem endopolyploidy along with specific leaf area and stem height, but increased days to bolting, leaf size, leaf dry mass, and leaf water content. Phenotypic responses to induced polyploidy varied significantly among accessions but this did not affect the relationship between phenotypic traits and endopolyploidy. Our results provide experimental support for a trade‐off between induced polyploidy and endopolyploidy, which caused induced polyploids to have lower endopolyploidy than diploids. Though polyploidy did not influence the relationship between endopolyploidy and plant traits, phenotypic responses to experimental genome duplication could not be easily predicted because of strong cytotype by accession interactions.     

Approaches to regulating recreational fisheries: balancing biology with angler satisfaction

Reviews in Fish Biology and Fisheries - Recreational fishing is practiced by?~?350 million people globally, and while it historically has been thought to have minimal ecological impact...     

Primate Location Preference in a Double-Tier Cage: The Effects of Illumination and Cage Height

Nonhuman primates are frequently housed in double-tier arrangements with significant differences between the environments of the upper and lower-row cages. Although several studies have investigated whether this arrangement alters monkeys' behavior, no studies have addressed the two most notable differences, light and height, individually to determine their relative importance. This experiment examined how rhesus and long-tailed macaques allocated their time between the upper and lower-row cages of a 1-over-1 apartment module under different lighting conditions. In Condition A, monkeys' baseline degree of preference for the upper- and lower-row was tested. In Condition B, the lighting environment was reversed by limiting illumination in the upper-row cage and increasing illumination in the lower-row cage. In both conditions, monkeys spent more time in the upper-row cage, thus indicating a strong preference for elevation regardless of illumination. The amount of time that monkeys spent in the lower-row cage increased by 7% under reversed lighting, but this trend was not significant. These results corroborate the importance of providing captive primates with access to elevated areas.     

The response of soil organic carbon of a rich fen peatland in interior Alaska to projected climate change

It is important to understand the fate of carbon in boreal peatland soils in response to climate change because a substantial change in release of this carbon as CO₂ and CH₄ could influence the climate system. The goal of this research was to synthesize the results of a field water table manipulation experiment conducted in a boreal rich fen into a process‐based model to understand how soil organic carbon (SOC) of the rich fen might respond to projected climate change. This model, the peatland version of the dynamic organic soil Terrestrial Ecosystem Model (peatland DOS‐TEM), was calibrated with data collected during 2005–2011 from the control treatment of a boreal rich fen in the Alaska Peatland Experiment (APEX). The performance of the model was validated with the experimental data measured from the raised and lowered water‐table treatments of APEX during the same period. The model was then applied to simulate future SOC dynamics of the rich fen control site under various CO₂ emission scenarios. The results across these emissions scenarios suggest that the rate of SOC sequestration in the rich fen will increase between year 2012 and 2061 because the effects of warming increase heterotrophic respiration less than they increase carbon inputs via production. However, after 2061, the rate of SOC sequestration will be weakened and, as a result, the rich fen will likely become a carbon source to the atmosphere between 2062 and 2099. During this period, the effects of projected warming increase respiration so that it is greater than carbon inputs via production. Although changes in precipitation alone had relatively little effect on the dynamics of SOC, changes in precipitation did interact with warming to influence SOC dynamics for some climate scenarios.     

Proteomic analysis of ERK1/2-mediated human sickle red blood cell membrane protein phosphorylation

Background

In sickle cell disease (SCD), the mitogen-activated protein kinase (MAPK) ERK1/2 is constitutively active and can be inducible by agonist-stimulation only in sickle but not in normal human red blood cells (RBCs). ERK1/2 is involved in activation of ICAM-4-mediated sickle RBC adhesion to the endothelium. However, other effects of the ERK1/2 activation in sickle RBCs leading to the complex SCD pathophysiology, such as alteration of RBC hemorheology are unknown.

Results

To further characterize global ERK1/2-induced changes in membrane protein phosphorylation within human RBCs, a label-free quantitative phosphoproteomic analysis was applied to sickle and normal RBC membrane ghosts pre-treated with U0126, a specific inhibitor of MEK1/2, the upstream kinase of ERK1/2, in the presence or absence of recombinant active ERK2. Across eight unique treatment groups, 375 phosphopeptides from 155 phosphoproteins were quantified with an average technical coefficient of variation in peak intensity of 19.8%. Sickle RBC treatment with U0126 decreased thirty-six phosphopeptides from twenty-one phosphoproteins involved in regulation of not only RBC shape, flexibility, cell morphology maintenance and adhesion, but also glucose and glutamate transport, cAMP production, degradation of misfolded proteins and receptor ubiquitination. Glycophorin A was the most affected protein in sickle RBCs by this ERK1/2 pathway, which contained 12 unique phosphorylated peptides, suggesting that in addition to its effect on sickle RBC adhesion, increased glycophorin A phosphorylation via the ERK1/2 pathway may also affect glycophorin A interactions with band 3, which could result in decreases in both anion transport by band 3 and band 3 trafficking. The abundance of twelve of the thirty-six phosphopeptides were subsequently increased in normal RBCs co-incubated with recombinant ERK2 and therefore represent specific MEK1/2 phospho-inhibitory targets mediated via ERK2.

Conclusions

These findings expand upon the current model for the involvement of ERK1/2 signaling in RBCs. These findings also identify additional protein targets of this pathway other than the RBC adhesion molecule ICAM-4 and enhance the understanding of the mechanism of small molecule inhibitors of MEK/1/2/ERK1/2, which could be effective in ameliorating RBC hemorheology and adhesion, the hallmarks of SCD.     

Perturbation Biology: Inferring Signaling Networks in Cellular Systems

We present a powerful experimental-computational technology for inferring network models that predict the response of cells to perturbations, and that may be useful in the design of combinatorial therapy against cancer. The experiments are systematic series of perturbations of cancer cell lines by targeted drugs, singly or in combination. The response to perturbation is quantified in terms of relative changes in the measured levels of proteins, phospho-proteins and cellular phenotypes such as viability. Computational network models are derived de novo, i.e., without prior knowledge of signaling pathways, and are based on simple non-linear differential equations. The prohibitively large solution space of all possible network models is explored efficiently using a probabilistic algorithm, Belief Propagation (BP), which is three orders of magnitude faster than standard Monte Carlo methods. Explicit executable models are derived for a set of perturbation experiments in SKMEL-133 melanoma cell lines, which are resistant to the therapeutically important inhibitor of RAF kinase. The resulting network models reproduce and extend known pathway biology. They empower potential discoveries of new molecular interactions and predict efficacious novel drug perturbations, such as the inhibition of PLK1, which is verified experimentally. This technology is suitable for application to larger systems in diverse areas of molecular biology.