Comparative transcriptomics of an arctic foundation species, tussock cottongrass (Eriophorum vaginatum), during an extreme heat event

© The Author(s), 2020. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Mohl, J. E., Fetcher, N., Stunz, E., Tang, J., & Moody, M. L. Comparative transcriptomics of an arctic foundation species, tussock cottongrass (Eriophorum vaginatum), during an extreme heat event. Scientific Reports, 10(1), (2020): 8990, doi:10.1038/s41598-020-65693-8.Tussock cottongrass (Eriophorum vaginatum) is a foundation species for much of the arctic moist acidic tundra, which is currently experiencing extreme effects of climate change. The Arctic is facing higher summer temperatures and extreme weather events are becoming more common. We used Illumina RNA-Seq to analyse cDNA libraries for differential expression of genes from leaves of ecologically well-characterized ecotypes of tussock cottongrass found along a latitudinal gradient in the Alaskan Arctic and transplanted into a common garden. Plant sampling was performed on a typical summer day and during an extreme heat event. We obtained a de novo assembly that contained 423,353 unigenes. There were 363 unigenes up-regulated and 1,117 down-regulated among all ecotypes examined during the extreme heat event. Of these, 26 HSP unigenes had >log2-fold up-regulation. Several TFs associated with heat stress in previous studies were identified that had >log2-fold up- or down-regulation during the extreme heat event (e.g., DREB, NAC). There was consistent variation in DEGs among ecotypes, but not specifically related to whether plants originated from taiga or tundra ecosystems. As the climate changes it is essential to determine ecotypic diversity at the genomic level, especially for widespread species that impact ecosystem function.We thank Thomas Parker for providing crucial logistical support at Toolik Field station and Darrel Dech, Stephen Turner, and Mayra Melendez for assistance in field sampling. Funding for this research was provided through the National Science Foundation (NSF/PLR 1418010 to NF, NSF/PLR 1417645 to MLM, NSF/PLR 1417763 to JT) and JEM received funding in part from NIH Grant #5G12RR007592 from the National Center for Research Resources (NCRR)/NIH to UTEP’s Border Biomedical Research Center. Significant logistic support came from Toolik Field Station and the Arctic LTER (NSF/PLR 1637459)

Landscape genomics provides evidence of ecotypic adaptation and a barrier to gene flow at treeline for the arctic foundation species Eriophorum vaginatum0

© The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Stunz, E., Fetcher, N., Lavretsky, P., Mohl, J., Tang, J., & Moody, M. Landscape genomics provides evidence of ecotypic adaptation and a barrier to gene flow at treeline for the arctic foundation species Eriophorum vaginatum. Frontiers in Plant Science, 13, (2022): 860439, https://doi.org/10.3389/fpls.2022.860439.Global climate change has resulted in geographic range shifts of flora and fauna at a global scale. Extreme environments, like the Arctic, are seeing some of the most pronounced changes. This region covers 14% of the Earth’s land area, and while many arctic species are widespread, understanding ecotypic variation at the genomic level will be important for elucidating how range shifts will affect ecological processes. Tussock cottongrass (Eriophorum vaginatum L.) is a foundation species of the moist acidic tundra, whose potential decline due to competition from shrubs may affect ecosystem stability in the Arctic. We used double-digest Restriction Site-Associated DNA sequencing to identify genomic variation in 273 individuals of E. vaginatum from 17 sites along a latitudinal gradient in north central Alaska. These sites have been part of 30 + years of ecological research and are inclusive of a region that was part of the Beringian refugium. The data analyses included genomic population structure, demographic models, and genotype by environment association. Genome-wide SNP investigation revealed environmentally associated variation and population structure across the sampled range of E. vaginatum, including a genetic break between populations north and south of treeline. This structure is likely the result of subrefugial isolation, contemporary isolation by resistance, and adaptation. Forty-five candidate loci were identified with genotype-environment association (GEA) analyses, with most identified genes related to abiotic stress. Our results support a hypothesis of limited gene flow based on spatial and environmental factors for E. vaginatum, which in combination with life history traits could limit range expansion of southern ecotypes northward as the tundra warms. This has implications for lower competitive attributes of northern plants of this foundation species likely resulting in changes in ecosystem productivity.This research was made possible by funding provided by NSF/PLR-1417645 to MM. The Botanical Society of America Graduate Student Research Award and the Dodson Research Grant from the Graduate School of the University of Texas at El Paso provided assistance to ES. The grant 5U54MD007592 from the National Institute on Minority Health and Health Disparities (NIMHD), a component of the National Institutes of Health (NIH) provided bioinformatics resources and support of JM

What Makes GPCRs from Different Families Bind to the Same Ligand?

G protein-coupled receptors (GPCRs) are the largest class of cell-surface receptor proteins with important functions in signal transduction and often serve as therapeutic drug targets. With the rapidly growing public data on three dimensional (3D) structures of GPCRs and GPCR-ligand interactions, computational prediction of GPCR ligand binding becomes a convincing option to high throughput screening and other experimental approaches during the beginning phases of ligand discovery. In this work, we set out to computationally uncover and understand the binding of a single ligand to GPCRs from several different families. Three-dimensional structural comparisons of the GPCRs that bind to the same ligand revealed local 3D structural similarities and often these regions overlap with locations of binding pockets. These pockets were found to be similar (based on backbone geometry and side-chain orientation using APoc), and they correlate positively with electrostatic properties of the pockets. Moreover, the more similar the pockets, the more likely a ligand binding to the pockets will interact with similar residues, have similar conformations, and produce similar binding affinities across the pockets. These findings can be exploited to improve protein function inference, drug repurposing and drug toxicity prediction, and accelerate the development of new drugs

RNA Post-transcriptional Modifications of an Early-Stage Large-Subunit Ribosomal Intermediate

Protein production by ribosomes is fundamental to life, and proper assembly of the ribosome is required for protein production. The RNA, which is post-transcriptionally modified, provides the platform for ribosome assembly. Thus, a complete understanding of ribosome assembly requires the determination of the RNA post-transcriptional modifications in all of the ribosome assembly intermediates and on each pathway. There are 26 RNA post-transcriptional modifications in 23S RNA of the mature Escherichia coli (E. coli) large ribosomal subunit. The levels of these modifications have been investigated extensively only for a small number of large subunit intermediates and under a limited number of cellular and environmental conditions. In this study, we determined the level of incorporations of 2-methyl adenosine, 3-methyl pseudouridine, 5-hydroxycytosine, and seven pseudouridines in an early-stage E. coli large-subunit assembly intermediate with a sedimentation coefficient of 27S. The 27S intermediate is one of three large subunit intermediates accumulated in E. coli cells lacking the DEAD-box RNA helicase DbpA and expressing the helicase inactive R331A DbpA construct. The majority of the investigated modifications are incorporated into the 27S large subunit intermediate to similar levels to those in the mature 50S large subunit, indicating that these early modifications or the enzymes that incorporate them play important roles in the initial events of large subunit ribosome assembly

Assessing changes in genomic divergence following a century of human-mediated secondary contact among wild and captive-bred ducks

Along with manipulating habitat, the direct release of domesticated individuals into the wild is a practice used worldwide to augment wildlife populations. We test between possible outcomes of human-mediated secondary contact using genomic techniques at both historical and contemporary timescales for two iconic duck species. First, we sequence several thousand ddRAD-seq loci for contemporary mallards (Anas platyrhynchos) throughout North America and two domestic mallard types (i.e., known game-farm mallards and feral Khaki Campbell's). We show that North American mallards may well be becoming a hybrid swarm due to interbreeding with domesticated game-farm mallards released for hunting. Next, to attain a historical perspective, we applied a bait-capture array targeting thousands of loci in century-old (1842-1915) and contemporary (2009-2010) mallard and American black duck (Anas rubripes) specimens. We conclude that American black ducks and mallards have always been closely related, with a divergence time of ~600,000 years before present, and likely evolved through prolonged isolation followed by limited bouts of gene flow (i.e., secondary contact). They continue to maintain genetic separation, a finding that overturns decades of prior research and speculation suggesting the genetic extinction of the American black duck due to contemporary interbreeding with mallards. Thus, despite having high rates of hybridization, actual gene flow is limited between mallards and American black ducks. Conversely, our historical and contemporary data confirm that the intensive stocking of game-farm mallards during the last ~100 years has fundamentally changed the genetic integrity of North America's wild mallard population, especially in the east. It thus becomes of great interest to ask whether the iconic North American mallard is declining in the wild due to introgression of maladaptive traits from domesticated forms. Moreover, we hypothesize that differential gene flow from domestic game-farm mallards into the wild mallard population may explain the overall temporal increase in differentiation between wild black ducks and mallards, as well as the uncoupling of genetic diversity and effective population size estimates across time in our results. Finally, our findings highlight how genomic methods can recover complex population histories by capturing DNA preserved in traditional museum specimens

JAK1 Pseudokinase V666G Mutant Dominantly Impairs JAK3 Phosphorylation and IL-2 Signaling

Overactive Janus kinases (JAKs) are known to drive leukemia, making them well-suited targets for treatment. We sought to identify new JAK-activating mutations and instead found a JAK1-inactivating pseudokinase mutation, V666G. In contrast to other pseudokinase mutations that canonically lead to an active kinase, the JAK1 V666G mutation led to under-activation seen by reduced phosphorylation. To understand the functional role of JAK1 V666G in modifying kinase activity we investigated its influence on other JAK kinases and within the Interleukin-2 pathway. JAK1 V666G not only inhibited its own activity, but its presence could inhibit other JAK kinases. These findings provide new insights into the potential of JAK1 pseudokinase to modulate its own activity, as well as of other JAK kinases. Thus, the features of the JAK1 V666 region in modifying JAK kinases can be exploited to allosterically inhibit overactive JAKs

Identification of a Potent Cytotoxic Pyrazole with Anti-Breast Cancer Activity That Alters Multiple Pathways

In this study, we identified a novel pyrazole-based derivative (P3C) that displayed potent cytotoxicity against 27 human cancer cell lines derived from different tissue origins with 50% cytotoxic concentrations (CC50) in the low micromolar and nanomolar range, particularly in two triple-negative breast cancer (TNBC) cell lines (from 0.25 to 0.49 µM). In vitro assays revealed that P3C induces reactive oxygen species (ROS) accumulation leading to mitochondrial depolarization and caspase-3/7 and -8 activation, suggesting the participation of both the intrinsic and extrinsic apoptotic pathways. P3C caused microtubule disruption, phosphatidylserine externalization, PARP cleavage, DNA fragmentation, and cell cycle arrest on TNBC cells. In addition, P3C triggered dephosphorylation of CREB, p38, ERK, STAT3, and Fyn, and hyperphosphorylation of JNK and NF-kB in TNBC cells, indicating the inactivation of both p38MAPK/STAT3 and ERK1/2/CREB signaling pathways. In support of our in vitro assays, transcriptome analyses of two distinct TNBC cell lines (MDA-MB-231 and MDA-MB-468 cells) treated with P3C revealed 28 genes similarly affected by the treatment implicated in apoptosis, oxidative stress, protein kinase modulation, and microtubule stability