Evidence for autotrophic growth of purple sulfur bacteria using pyrite as electron and sulfur source.

Purple sulfur bacteria (PSB) are capable of anoxygenic photosynthesis via oxidizing reduced sulfur compounds and are considered key drivers of the sulfur cycle in a range of anoxic environments. In this study, we show that Allochromatium vinosum (a PSB species) is capable of autotrophic growth using pyrite as the electron and sulfur source. Comparative growth profile, substrate characterization, and transcriptomic sequencing data provided valuable insight into the molecular mechanisms underlying the bacterial utilization of pyrite and autotrophic growth. Specifically, the pyrite-supported cell cultures ("py"') demonstrated robust but much slower growth rates and distinct patterns from their sodium sulfide-amended positive controls. Up to ~200-fold upregulation of genes encoding various c- and b-type cytochromes was observed in "py," pointing to the high relevance of these molecules in scavenging and relaying electrons from pyrite to cytoplasmic metabolisms. Conversely, extensive downregulation of genes related to LH and RC complex components indicates that the electron source may have direct control over the bacterial cells' photosynthetic activity. In terms of sulfur metabolism, genes encoding periplasmic or membrane-bound proteins (e.g., FccAB and SoxYZ) were largely upregulated, whereas those encoding cytoplasmic proteins (e.g., Dsr and Apr groups) are extensively suppressed. Other notable differentially expressed genes are related to flagella/fimbriae/pilin(+), metal efflux(+), ferrienterochelin(-), and [NiFe] hydrogenases(+). Characterization of the biologically reacted pyrite indicates the presence of polymeric sulfur. These results have, for the first time, put the interplay of PSB and transition metal sulfide chemistry under the spotlight, with the potential to advance multiple fields, including metal and sulfur biogeochemistry, bacterial extracellular electron transfer, and artificial photosynthesis. IMPORTANCE: Microbial utilization of solid-phase substrates constitutes a critical area of focus in environmental microbiology, offering valuable insights into microbial metabolic processes and adaptability. Recent advancements in this field have profoundly deepened our knowledge of microbial physiology pertinent to these scenarios and spurred innovations in biosynthesis and energy production. Furthermore, research into interactions between microbes and solid-phase substrates has directly linked microbial activities to the surrounding mineralogical environments, thereby enhancing our understanding of the relevant biogeochemical cycles. Our study represents a significant step forward in this field by demonstrating, for the first time, the autotrophic growth of purple sulfur bacteria using insoluble pyrite (FeS2) as both the electron and sulfur source. The presented comparative growth profiles, substrate characterizations, and transcriptomic sequencing data shed light on the relationships between electron donor types, photosynthetic reaction center activities, and potential extracellular electron transfer in these organisms capable of anoxygenic photosynthesis. Furthermore, the findings of our study may provide new insights into early-Earth biogeochemical evolutions, offering valuable constraints for understanding the environmental conditions and microbial processes that shaped our planet's history.

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.

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.

ISOGlyP: de novo prediction of isoform-specific mucin-type O-glycosylation.

Mucin-type O-glycosylation is one of the most common posttranslational modifications of proteins. The abnormal expression of various polypeptide GalNAc-transferases (GalNAc-Ts) which initiate and define sites of O-glycosylation are linked to many cancers and other diseases. Current O-glycosyation prediction programs utilize O-glycoproteomics data obtained without regard to the transferase isoform (s) responsible for the glycosylation. With 20 different GalNAc-Ts in humans, having an ability to predict and interpret O-glycosylation sites in terms of specific GalNAc-T isoforms is invaluable. To fill this gap, ISOGlyP (Isoform-Specific O-Glycosylation Prediction) has been developed. Using position-specific enhancement values generated based on GalNAc-T isoform-specific amino acid preferences, ISOGlyP predicts the propensity that a site would be glycosylated by a specific transferase. ISOGlyP gave an overall prediction accuracy of 70% against in vivo data, which is comparable to that of the NetOGlyc4.0 predictor. Additionally, ISOGlyP can identify the known effects of long- and short-range prior glycosylation and can generate potential peptide sequences selectively glycosylated by specific isoforms. ISOGlyP is freely available for use at ISOGlyP.utep.edu. The code is also available on GitHub (https://github.com/jonmohl/ISOGlyP).

Schlafen 11 Restricts Flavivirus Replication.

Schlafen 11 (Slfn11) is an interferon-stimulated gene that controls the synthesis of proteins by regulating tRNA abundance. Likely through this mechanism, Slfn11 has previously been shown to impair human immunodeficiency virus type 1 (HIV-1) infection and the expression of codon-biased open reading frames. Because replication of positive-sense single-stranded RNA [(+)ssRNA] viruses requires the immediate translation of the incoming viral genome, whereas negative-sense single-stranded RNA [(-)ssRNA] viruses carry at infection an RNA replicase that makes multiple translation-competent copies of the incoming viral genome, we reasoned that (+)ssRNA viruses will be more sensitive to the effect of Slfn11 on protein synthesis than (-)ssRNA viruses. To evaluate this hypothesis, we tested the effects of Slfn11 on the replication of a panel of ssRNA viruses in the human glioblastoma cell line A172, which naturally expresses Slfn11. Depletion of Slfn11 significantly increased the replication of (+)ssRNA viruses from the Flavivirus genus, including West Nile virus (WNV), dengue virus (DENV), and Zika virus (ZIKV), but had no significant effect on the replication of the (-)ssRNA viruses vesicular stomatitis virus (VSV) (Rhabdoviridae family) and Rift Valley fever virus (RVFV) (Phenuiviridae family). Quantification of the ratio of genome-containing viral particles to PFU indicated that Slfn11 impairs WNV infectivity. Intriguingly, Slfn11 prevented WNV-induced downregulation of a subset of tRNAs implicated in the translation of 11.8% of the viral polyprotein. Low-abundance tRNAs might promote optimal protein folding and enhance viral infectivity, as previously reported. In summary, this study demonstrates that Slfn11 restricts flavivirus replication by impairing viral infectivity.IMPORTANCE We provide evidence that the cellular protein Schlafen 11 (Slfn11) impairs replication of flaviviruses, including West Nile virus (WNV), dengue virus (DENV), and Zika virus (ZIKV). However, replication of single-stranded negative RNA viruses was not affected. Specifically, Slfn11 decreases the infectivity of WNV potentially by preventing virus-induced modifications of the host tRNA repertoire that could lead to enhanced viral protein folding. Furthermore, we demonstrate that Slfn11 is not the limiting factor of this novel broad antiviral pathway.

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.

Predicting mucin-type O-Glycosylation using enhancement value products from derived protein features.

Mucin-type O-glycosylation is one of the most common post-translational modifications of proteins. This glycosylation is initiated in the Golgi by the addition of the sugar N-acetylgalactosamine (GalNAc) onto protein Ser and Thr residues by a family of polypeptide GalNAc transferases. In humans there are 20 isoforms that are differentially expressed across tissues that serve multiple important biological roles. Using random peptide substrates, isoform specific amino acid preferences have been obtained in the form of enhancement values (EV). These EVs alone have previously been used to predict O-glycosylation sites via the web based ISOGlyP (Isoform Specific O-Glycosylation Prediction) tool. Here we explore additional protein features to determine whether these can complement the random peptide derived enhancement values and increase the predictive power of ISOGlyP. The inclusion of additional protein substrate features (such as secondary structure and surface accessibility) was found to increase sensitivity with minimal loss of specificity, when tested with three different published in vivo O-glycoproteomics data sets, thus increasing the overall accuracy of the ISOGlyP predictions.

ISOGlyP: O-Glycosylation Site Prediction Using Peptide Sequences and GALNTs.

Mucin-type O-glycosylation is one of the most common posttranslational modifications of proteins. The abnormal expression of various polypeptide GalNAc-transferases (GALNTs) which initiate and define sites of O-glycosylation is linked to many cancers and other diseases. Many current O-glycosylation prediction programs utilize O-glycoproteomics data obtained without using the transferase isoform(s) responsible for the glycosylation. With 20 different GALNTs in humans, having the ability to predict and interpret O-glycosylation sites in terms of specific GALNT isoforms is invaluable.To fill this gap, ISOGlyP (isoform-specific O-glycosylation prediction) has been developed. Using position-specific enhancement values generated based on GalNAc-T isoform-specific amino acid preferences, ISOGlyP predicts the propensity that a site would be glycosylated by a specific transferase. ISOGlyP gave an overall prediction accuracy of 70% against in vivo data, which is comparable to that of the NetOGlyc4.0 predictor. Additionally, ISOGlyP can identify the known effects of long- and short-range prior glycosylation and can generate potential peptide sequences selectively glycosylated by specific isoforms. ISOGlyP is freely available for use at https://ISOGlyP.utep.edu . The code is also available on GitHub ( https://github.com/jonmohl/ISOGlyP ).

Chromosomal-level reference genome of a wild North American mallard (Anas platyrhynchos).

The mallard (Anas platyrhynchos) is one of the most common, economically, and socially important birds around the world. Mallards were not only an important food source for early humans but eventually becoming intimately linked with people as they were domesticated over the last 2,000 years. To date, mallard genomes are largely reconstructed from samples of domestic or unknown genetic heritage. Here, we report the first high-quality genome assembly and annotation of a genetically vetted wild mallard from North America (NAwild_v1.0). The genome was assembled using a combination of shotgun libraries, proximity ligation Chicago, and Dovetail Hi-C libraries. The final assembly is ∼1.04 Gb in size, with 98.3% of the sequence located in 30 full or nearly full chromosome-level scaffolds, and with a N50/L50 of 79.1 Mb/4 scaffolds. We used a combination of gene prediction and similarity approaches to annotate a total of 23,584 functional genes, of which 19,242 were associated to GO terms. The genome assembly and the set of annotated genes yielded a 95.4% completeness score when compared with the BUSCO aves_odb10 dataset. Next, we aligned 3 previously published mallard genomes to ours, and demonstrate how runs of homozygosity and nucleotide diversity are substantially higher and lower, respectively, to ours and how these artificially changed genomes resulted in profoundly different and biased demographic histories. Our wild mallard assembly not only provides a valuable resource to shed light onto genome evolution, speciation, and other adaptive processes, but also helping with identifying functional genes that have been significantly altered during the domestication process.

The meaning of wild: Genetic and adaptive consequences from large-scale releases of domestic mallards.

The translocation of individuals around the world is leading to rising incidences of anthropogenic hybridization, particularly between domestic and wild congeners. We apply a landscape genomics approach for thousands of mallard (Anas platyrhynchos) samples across continental and island populations to determine the result of over a century of supplementation practices. We establish that a single domestic game-farm mallard breed is the source for contemporary release programs in Eurasia and North America, as well as for established feral populations in New Zealand and Hawaii. In particular, we identify central Europe and eastern North America as epicenters of ongoing anthropogenic hybridization, and conclude that the release of game-farm mallards continues to affect the genetic integrity of wild mallards. Conversely, self-sustaining feral populations in New Zealand and Hawaii not only show strong differentiation from their original stock, but also signatures of local adaptation occurring in less than a half-century since game-farm mallard releases have ceased. We conclude that 'wild' is not singular, and that even feral populations are capable of responding to natural processes. Although considered paradoxical to biological conservation, understanding the capacity for wildness among feral and feral admixed populations in human landscapes is critical as such interactions increase in the Anthropocene.

Rapid Shift from SARS-CoV-2 Delta to Omicron Sub-Variants within a Dynamic Southern U.S. Borderplex.

COVID-19, caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), remains an ongoing global health challenge. This study analyzed 3641 SARS-CoV-2 positive samples from the El Paso, Texas, community and hospitalized patients over 48 weeks from Fall 2021 to Summer 2022. The binational community along the U.S. southern border was predominantly SARS-CoV-2 Delta variant (B.1.617.2) positive for a 5-week period from September 2021 to January 2022 and quickly transitioned to the Omicron variant (B.1.1.529), which was first detected at the end of December 2021. Omicron replaced Delta as the predominant detectable variant in the community and was associated with a sharp increase in COVID-19 positivity rate, related hospitalizations, and newly reported cases. In this study, Omicron BA.1, BA.4, and BA.5 variants were overwhelmingly associated with S-gene dropout by qRT-PCR analysis unlike the Delta and Omicron BA.2 variants. The study reveals that a dominant variant, like Delta, can be rapidly replaced by a more transmissible variant, like Omicron, within a dynamic metropolitan border city, necessitating enhanced monitoring, readiness, and response from public health officials and healthcare workers.

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.

Landscape Genomics Provides Evidence of Ecotypic Adaptation and a Barrier to Gene Flow at Treeline for the Arctic Foundation Species Eriophorum vaginatum.

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.

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.

GPCR-PEnDB: a database of protein sequences and derived features to facilitate prediction and classification of G protein-coupled receptors.

G protein-coupled receptors (GPCRs) constitute the largest group of membrane receptor proteins in eukaryotes. Due to their significant roles in various physiological processes such as vision, smell and inflammation, GPCRs are the targets of many prescription drugs. However, the functional and sequence diversity of GPCRs has kept their prediction and classification based on amino acid sequence data as a challenging bioinformatics problem. There are existing computational approaches, mainly using machine learning and statistical methods, to predict and classify GPCRs based on amino acid sequence and sequence derived features. In this paper, we describe a searchable MySQL database, named GPCR-PEnDB (GPCR Prediction Ensemble Database), of confirmed GPCRs and non-GPCRs. It was constructed with the goal of allowing users to conveniently access useful information of GPCRs in a wide range of organisms and to compile reliable training and testing datasets for different combinations of computational tools. This database currently contains 3129 confirmed GPCR and 3575 non-GPCR sequences collected from the UniProtKB/Swiss-Prot protein database, encompassing over 1200 species. The non-GPCR entries include transmembrane proteins for evaluating various prediction programs' abilities to distinguish GPCRs from other transmembrane proteins. Each protein is linked to information about its source organism, classification, sequence lengths and composition, and other derived sequence features. We present examples of using this database along with its graphical user interface, to query for GPCRs with specific sequence properties and to compare the accuracies of five tools for GPCR prediction. This initial version of GPCR-PEnDB will provide a framework for future extensions to include additional sequence and feature data to facilitate the design and assessment of software tools and experimental studies to help understand the functional roles of GPCRs. Database URL: gpcr.utep.edu/database.

Thermal Bioprinting Causes Ample Alterations of Expression of LUCAT1, IL6, CCL26, and NRN1L Genes and Massive Phosphorylation of Critical Oncogenic Drug Resistance Pathways in Breast Cancer Cells.

Bioprinting technology merges engineering and biological fields and together, they possess a great translational potential, which can tremendously impact the future of regenerative medicine and drug discovery. However, the molecular effects elicited by thermal inkjet bioprinting in breast cancer cells remains elusive. Previous studies have suggested that bioprinting can be used to model tissues for drug discovery and pharmacology. We report viability, apoptosis, phosphorylation, and RNA sequence analysis of bioprinted MCF7 breast cancer cells at separate timepoints post-bioprinting. An Annexin A5-FITC apoptosis stain was used in combination with flow cytometry at 2 and 24 h post-bioprinting. Antibody arrays using a Human phospho-MAPK array kit was performed 24 h post-bioprinting. RNA sequence analysis was conducted in samples collected at 2, 7, and 24 h post-bioprinting. The post-bioprinting cell viability averages were 77 and 76% at 24 h and 48 h, with 31 and 64% apoptotic cells at 2 and 24 h after bioprinting. A total of 21 kinases were phosphorylated in the bioprinted cells and 9 were phosphorylated in the manually seeded controls. The RNA seq analysis in the bioprinted cells identified a total of 12,235 genes, of which 9.7% were significantly differentially expressed. Using a ±2-fold change as the cutoff, 266 upregulated and 206 downregulated genes were observed in the bioprinted cells, with the following 5 genes uniquely expressed NRN1L, LUCAT1, IL6, CCL26, and LOC401585. This suggests that thermal inkjet bioprinting is stimulating large scale gene alterations that could potentially be utilized for drug discovery. Moreover, bioprinting activates key pathways implicated in drug resistance, cell motility, proliferation, survival, and differentiation.

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

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.