Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Expression”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Concentration-response gene expression analysis in zebrafish reveals phenotypically-anchored transcriptional responses to retene

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous environmental contaminants and are associated with human disease. Canonically, many PAHs induce toxicity via activation of the aryl hydrocarbon receptor (AHR) pathway. While the interaction between PAHs and the AHR is well-established, understanding which AHR-regulated transcriptional effects directly result in observable phenotypes and which are adaptive or benign is important to better understand PAH toxicity. Retene is a frequently detected PAH in environmental sampling and has been associated with AHR2-dependent developmental toxicity in zebrafish, though its mechanism of toxicity has not been fully elucidated. To interrogate transcriptional changes causally associated with retene toxicity, we conducted whole-animal RNA sequencing at 48 hours post-fertilization after exposure to eight retene concentrations. The concentrations were selected to produce effects ranging from no phenotype to mortality and malformations in 100% of animals at 5 days post-fertilization. We identified a concentration-response relationship between retene teratogenicity and differential gene expression in both number of DEGs and magnitude of expression change. Elevated expression of cyp1a at retene concentrations below the threshold for teratogenicity suggested that while cyp1a expression is a sensitive biomarker of AHR activation, it may be too sensitive to serve as a biomarker of AHR-dependent teratogenicity. Genes differentially expressed at only non-teratogenic concentrations were enriched for transforming growth factor-ß (TGF-ß) signaling pathway disruption while DEGs identified at only teratogenic concentrations were significantly enriched for response to xenobiotic stimulus and reduction-oxidation reaction activity. DEGs which spanned both non-teratogenic and teratogenic concentrations showed similar disrupted biological processes to those unique to teratogenic concentrations, indicating these processes were disrupted at low exposure concentrations. Gene co-expression network analysis identified several gene modules, including those associated with PAHs and AHR2 activation. One, Module 7, was strongly enriched for AHR2-associated genes and contained the strongest responses to retene. Benchmark concentration (BMC) of Module 7 genes identified a median BMC of 7.5 µM, nearly the highest retene concentration with no associated teratogenicity, supporting the hypothesis that Module 7 genes are largely responsible for retene toxicity.

Toxin, Zebrafish, Retene, Transcriptomics, network↗

Data for FUN-PROSE: A Deep Learning Approach to Predict Condition-Specific Gene Expression in Fungi

mRNA levels of all genes in a genome is a critical piece of information defining the overall state of the cell in a given environmental condition. Being able to reconstruct such condition-specific expression in fungal genomes is particularly important to metabolically engineer these organisms to produce desired chemicals in industrially scalable conditions. Most previous deep learning approaches focused on predicting the average expression levels of a gene based on its promoter sequence, ignoring its variation across different conditions. Here we present FUN-PROSE—a deep learning model trained to predict differential expression of individual genes across various conditions using their promoter sequences and expression levels of all transcription factors. We train and test our model on three fungal species and get the correlation between predicted and observed condition-specific gene expression as high as 0.85. We then interpret our model to extract promoter sequence motifs responsible for variable expression of individual genes. We also carried out input feature importance analysis to connect individual transcription factors to their gene targets. A sizeable fraction of both sequence motifs and TF-gene interactions learned by our model agree with previously known biological information, while the rest corresponds to either novel biological facts or indirect correlations.

Genomics↗

Data for A Fluorescence-Based Transient Expression Assay for the Analysis of Upstream Open Reading Frames in Plants

Scripts for the manuscript "A fluorescence-based transient expression assay for the analysis of upstream open reading frames in plant" by Haas et al. Upstream open reading frames (uORFs) are regulatory elements present in the 5′ leaders of mRNA that can significantly impact downstream gene expression in eukaryotes. In crop engineering, editing of uORFs can provide an avenue to upregulate expression of native genes without the need to add persistent transgenic copies. Even with genome- wide methods to identify translated uORFs such as ribosome profiling, their functional characterization depends on validation through reporter gene assays and mutagenesis studies. Current screening methods for plants use luciferases or protoplasts to measure differential gene expression between wild- type and mutated transcript leaders, which requires tissue processing and/or substrate addition. Here, we present a time- and cost- efficient alternative to investigate transcript leaders by co- expression of two fluorescent proteins in Nicotiana benthamiana leaf tissue and test our assay on genes involved in photoprotection, editing of which could provide a pathway to increase CO2 assimilation during sun–shade transitions.

Gene Editing↗

Data for Expression of a Bacterial Trehalose 6-Phosphate Synthase Gene otsA in Camelina sativa Seeds Promotes the Channelling of Carbon Towards Oil Accumulation

Improving seed oil yield is essential for developing Camelina sativa as a sustainable biofuel crop. Fatty acid synthesis depends on the production of acetyl-CoA from photosynthetically derived sugars. Trehalose 6-phosphate (T6P), a proxy for sucrose availability, can link sugar status to plant growth and development. Synthesised by trehalose 6-phosphate synthase (TPS) from UDP-glucose and glucose-6-phosphate, T6P plays a regulatory role in metabolism. Our previous studies on Arabidopsis transgenic lines constitutively expressing the E. coli otsA (encoding TPS) showed increased T6P levels and seed triacylglycerol, along with stunted growth. In the present study we express otsA in camelina under the control of a seed-specific Phaseolin promoter. Seeds of the resulting transgenic lines accumulated high levels of T6P, and a 15%–20% increase in total fatty acids and triacylglycerol compared to wild-type. Molecular analysis showed the transgenic seeds had reduced SnRK1 activity, elevated WRI1 protein levels, and increased the levels of WRI1 and its target genes, along with enhanced rates of fatty acid synthesis that increased seed weights relative to wild type. Notably, the increase in oil did not affect seed protein levels but did reduce the soluble metabolite fraction. Crucially, seed-specific expression of otsA mitigated the growth defects associated with constitutive otsA expression, and the transgenic lines showed normal seed development and germination. These findings demonstrate that targeted T6P modulation via seed-specific otsA expression is an effective metabolic engineering strategy to boost oil production in camelina and potentially in other oilseed crops and bioenergy crops such as energycane, sorghum and miscanthus.

Lipids↗

Data for A Role for Differential Rubisco Activase Isoform Expression in C4 Bioenergy Grasses at High Temperature

Rubisco activase (Rca) facilitates the release of sugar-phosphate inhibitors at Rubisco catalytic sites during CO2 fixation. Most plant species express two Rca isoforms, the larger Rca-α and the shorter Rca-β, either by alternative splicing from a single gene or expression from separate genes. The mechanism of Rubisco activation by Rca isoforms has been intensively studied in C3 plants. However, the functional role of Rca in C4 plants where Rubisco and Rca are located in a much higher [CO2] compartment is less clear. In this study, we selected four C4 bioenergy grasses and the model C4 grass setaria ( Setaria viridis ) to investigate the role of Rca in C4 photosynthesis. All five C4 grass species contained two Rca genes, one encoding Rca-α and the other Rca-β, which were positioned closely together in the genomes. A variety of abiotic stress-related motifs were identified in the Rca-α promoter of each grass, and while the Rca-β gene was constantly highly expressed at ambient temperature, Rca-α isoforms were expressed only at high temperature but never surpassed 30% of Rca-β content. The pattern of Rca-α induction on transition to high temperature and reduction on return to ambient temperature was the same in all five C4 grasses. In sorghum ( Sorghum bicolor ), sugarcane ( Saccharum officinarum ), and setaria, the induction rate of Rca-α was similar to the recovery rate of photosynthesis and Rubisco activation at high temperature. This association between Rca-α isoform expression and maintenance of Rubisco activation at high temperature suggests that Rca-α has a functional thermo-protective role in carbon fixation in C4 grasses by sustaining Rubisco activation at high temperature.

Genomics↗

Multiple highly expressed phospho enol pyruvate carboxylase genes have divergent enzyme kinetic properties in two C4 grasses

Abstract Background and Aims Phosphoenolpyruvate (PEP) carboxylase (PEPC) catalyses the irreversible carboxylation of PEP with bicarbonate to produce oxaloacetate. This reaction powers the carbon-concentrating mechanism (CCM) in plants that perform C4 photosynthesis. This CCM is generally driven by a single PEPC gene product that is highly expressed in the cytosol of mesophyll cells. We found two C4 grasses, Panicum miliaceum and Echinochloa colona, that each have two highly expressed PEPC genes. We characterized the kinetic properties of the two most abundant PEPCs in E. colona and P. miliaceum to better understand how the enzyme’s amino acid structure influences its function. Methods Coding sequences of the two most abundant PEPC proteins in E. colona and P. miliaceum were synthesized by GenScript and were inserted into bacteria expression plasmids. Point mutations resulting in substitutions at conserved amino acid residues (e.g. N-terminal serine and residue 890) were created via site-directed PCR mutagenesis. The kinetic properties of semi-purified plant PEPCs from Escherichia coli were analysed using membrane-inlet mass spectrometry and a spectrophotometric enzyme-coupled reaction. Key Results The two most abundant P. miliaceum PEPCs (PmPPC1 and PmPPC2) have similar sequence identities (>95 %), and as a result had similar kinetic properties. The two most abundant E. colona PEPCs (EcPPC1 and EcPPC2) had identities of ~78 % and had significantly different kinetic properties. The PmPPCs and EcPPCs had different responses to allosteric inhibitors and activators, and substitutions at the conserved N-terminal serine and residue 890 resulted in significantly altered responses to allosteric regulators. Conclusions The two, significantly expressed C4Ppc genes in P. miliaceum were probably the result of genomes combining from two closely related C4Panicum species. We found natural variation in PEPC’s sensitivity to allosteric inhibition that seems to bypass the conserved 890 residue, suggesting alternative evolutionary pathways for increased malate tolerance and other kinetic properties.

DiMario, Robert J. (ORCID:0000000250566868)↗

Comparing four heat-inducible promoters in stably transformed sugarcane regarding spatial and temporal control of transgene expression reveals candidates to drive stem-preferred transgene expression

Small heat shock protein (sHSP) promoters contain cis-regulatory elements that facilitate transcription in response to heat stress, making them valuable tools for functional studies through controlled gene expression and the precise regulation of gene-editing tools or morphogenic regulators. To evaluate their utility, GUS reporter gene expression driven by four plant-sourced HSP promoters (p GmHSP17.5 , p HvHSP17 , p ZmHSP17.7 , and p ZmHSP26 ) was compared across various tissues of stably transformed sugarcane before and after heat treatment. At 22°C, all promoters showed minimal activity in leaves and roots, although p ZmHSP17.7 and p HvHSP17 displayed moderate expression in stems. Following heat treatment, all promoters exhibited their highest activity in stems, followed by leaves and roots. In stem tissues, p GmHSP17.5 displayed heat-induced uidA expression comparable to the constitutive p ZmUbi promoter. Notably, heat-induced reporter gene activity in stem middle sections of single-copy transgenic lines containing p ZmHSP17.7 , p HvHSP17 , or p ZmHSP26 exceeded p ZmUbi -derived uidA activity by 9.7-fold, 3.8-fold, and 3.0-fold, respectively, with 346- to 3,672-fold induction compared to control conditions. Most promoters showed peak expression in the middle sections of the stem, while p HvHSP17 was the most active in the stem apices. Histochemical analysis revealed that p ZmHSP17.7 and p HvHSP17 were active in both parenchyma cells and vascular bundles within sugarcane stems. Among leaf tissues, mature leaves exhibited greater expression than senescing or immature leaves, while root activity remained consistently minimal across all promoters. Temperature-course experiments identified distinct activation thresholds: 34°C–36°C for p ZmHSP17.7 , 36°C for p ZmHSP26 , 36°C–38°C for p HvHSP17 , and 40°C–42°C for p GmHSP17.5 . Drought stress also induced reporter gene transcription in stems under HSP promoters, although with lower fold induction than heat treatment. These findings provide valuable tools for gene function studies and biotechnology applications, including heat stress tolerance research, controlled transgene expression in metabolic engineering, precision gene editing, and developmental biology studies.

60 APPLIED LIFE SCIENCES↗

Nonadditive gene expression is correlated with nonadditive phenotypic expression in interspecific triploid hybrids of willow ( Salix spp.)

Many studies have highlighted the complex and diverse basis for heterosis in inbred crops. Despite the lack of a consensus model, it is vital that we turn our attention to understanding heterosis in undomesticated, heterozygous, and polyploid species, such as willow (Salix spp.). Shrub willow is a dedicated energy crop bred to be fast-growing and high yielding on marginal land without competing with food crops. A trend in willow breeding is the consistent pattern of heterosis in triploids produced from crosses between diploid and tetraploid species. Here, we test whether differentially expressed genes are associated with heterosis in triploid families derived from diploid Salix purpurea, diploid Salix viminalis, and tetraploid Salix miyabeana parents. Three biological replicates of shoot tips from all family progeny and parents were collected after 12 weeks in the greenhouse and RNA extracted for RNA-Seq analysis. This study provides evidence that nonadditive patterns of gene expression are correlated with nonadditive phenotypic expression in interspecific triploid hybrids of willow. Expression-level dominance was most correlated with heterosis for biomass yield traits and was highly enriched for processes involved in starch and sucrose metabolism. In addition, there was a global dosage effect of parent alleles in triploid hybrids, with expression proportional to copy number variation. Importantly, differentially expressed genes between family parents were most predictive of heterosis for both field and greenhouse collected traits. Altogether, these data will be used to progress models of heterosis to complement the growing genomic resources available for the improvement of heterozygous perennial bioenergy crops.

59 BASIC BIOLOGICAL SCIENCES↗

Expression of blue pigment synthetase a from Streptomyces lavenduale reveals insights on the effects of refactoring biosynthetic megasynthases for heterologous expression in Escherichia coli .

High GC bacteria from the genus Streptomyces harbor expansive secondary metabolism. The expression of biosynthetic proteins and the characterization and identification of biological "parts" for synthetic biology purposes from such pathways are of interest. However, the high GC content of proteins from actinomycetes in addition to the large size and multi-domain architecture of many biosynthetic proteins (such as non-ribosomal peptide synthetases; NRPSs, and polyketide synthases; PKSs often called "megasynthases") often presents issues with full-length translation and folding. Here we evaluate a non-ribosomal peptide synthetase (NRPS) from Streptomyces lavenduale, a multidomain "megasynthase" gene that comes from a high GC (72.5%) genome. While a preliminary step in revealing differences, to our knowledge this presents the first head-to-head comparison of codon-optimized sequences versus a native sequence of proteins of streptomycete origin heterologously expressed in E. coli. We found that any disruption in co-translational folding from codon mismatch that reduces the titer of indigoidine is explainable via the formation of more inclusion bodies as opposed to compromising folding or posttranslational modification in the soluble fraction. In conclusion, this result supports that one could apply any refactoring strategies that improve soluble expression in E. coli without concern that the protein that reaches the soluble fraction is differentially folded.

59 BASIC BIOLOGICAL SCIENCES↗

Expression of V-nitrogenase and Fe-nitrogenase in Methanosarcina acetivorans is controlled by molybdenum, fixed nitrogen, and the expression of Mo-nitrogenase

ABSTRACT All nitrogen-fixing bacteria and archaea (diazotrophs) use molybdenum (Mo) nitrogenase to reduce dinitrogen (N 2 ) to ammonia, with some also containing vanadium (V) and iron-only (Fe) nitrogenases that lack Mo. Among diazotrophs, the regulation and usage of the alternative V-nitrogenase and Fe-nitrogenase in methanogens are largely unknown. Methanosarcina acetivorans contains nif , vnf , and anf gene clusters encoding putative Mo-nitrogenase, V-nitrogenase, and Fe-nitrogenase, respectively. This study investigated nitrogenase expression and growth by M. acetivorans in response to fixed nitrogen, Mo/V availability, and CRISPRi repression of the nif , vnf , and/or anf gene clusters. The availability of Mo and V significantly affected growth of M. acetivorans with N 2 but not with NH 4 Cl. M. acetivorans exhibited the fastest growth rate and highest cell yield during growth with N 2 in medium containing Mo, and the slowest growth in medium lacking Mo and V. qPCR analysis revealed the transcription of the nif operon is only moderately affected by depletion of fixed nitrogen and Mo, whereas vnf and anf transcription increased significantly when fixed nitrogen and Mo were depleted, with removal of Mo being key. Immunoblot analysis revealed Mo-nitrogenase is detected when fixed nitrogen is depleted regardless of Mo availability, while V-nitrogenase and Fe-nitrogenase are detected only in the absence of fixed nitrogen and Mo. CRISPRi repression studies revealed that V-nitrogenase and/or Fe-nitrogenase are required for Mo-independent diazotrophy, and unexpectedly that the expression of Mo-nitrogenase is also required. These results reveal that alternative nitrogenase production in M. acetivorans is tightly controlled and dependent on Mo-nitrogenase expression. IMPORTANCE Methanogens and closely related methanotrophs are the only archaea known or predicted to possess nitrogenase. Methanogens play critical roles in both the global biological nitrogen and carbon cycles. Moreover, methanogens are an ancient microbial lineage and nitrogenase likely originated in methanogens. An understanding of the usage and properties of nitrogenases in methanogens can provide new insight into the evolution of nitrogen fixation and aid in the development nitrogenase-based biotechnology. This study provides the first evidence that a methanogen can produce all three forms of nitrogenases, including simultaneously. The results reveal components of Mo-nitrogenase regulate or are needed to produce V-nitrogenase and Fe-nitrogenase in methanogens, a result not seen in bacteria. Overall, this study provides a foundation to understand the assembly, regulation, and activity of the alternative nitrogenases in methanogens.

Biotechnology & Applied Microbiology↗

Mechanically induced localisation of SECONDARY WALL INTERACTING bZIP is associated with thigmomorphogenic and secondary cell wall gene expression

Plant growth requires the integration of internal and external cues, perceived and transduced into a developmental programme of cell division, elongation and wall thickening. Mechanical forces contribute to this regulation, and thigmomorphogenesis typically includes reducing stem height, increasing stem diameter, and a canonical transcriptomic response. We present data on a bZIP transcription factor involved in this process in grasses. Brachypodium distachyon SECONDARY WALL INTERACTING bZIP (SWIZ) protein translocated into the nucleus following mechanostimulation. Classical touch-responsive genes were upregulated in B. distachyon roots following touch, including significant induction of the glycoside hydrolase 17 family, which may be unique to grass thigmomorphogenesis. SWIZ protein binding to an E-box variant in exons and introns was associated with immediate activation followed by repression of gene expression. SWIZ overexpression resulted in plants with reduced stem and root elongation. These data further define plant touch-responsive transcriptomics and physiology, offering insights into grass mechanotranduction dynamics.

Coomey, Joshua↗

Comparison of gene expression in the skin tissue of gray, humpback, and fin whales

Analyses of gene expression in the skin of several species of whales identified genes that are differentially expressed in association with environmental factors, suggesting that skin transcriptomics may provide a valuable tool for assessing physiological responses in marine mammals. Previous work exploring differing levels of gene expression has focused on odontocetes with comparatively limited investigation of skin gene expression has been explored in mysticetes. Here, we describe the identity of genes expressed in skin tissue of three species of baleen whales to establish a baseline of gene expression and compare gene identity and expression patterns across species. We also evaluate sex-specific differences in skin gene expression through a comparison of expression levels between males and females in gray and humpback whales. A total of 16 skin tissue samples were collected from free-ranging gray, humpback and fin whales off the central Oregon coast in the eastern North Pacific. Comparison of the expressed genes in the humpback and gray whale skin tissue to the blue whale reference database identified enriched gene ontology terms in the skin tissue of each species, suggesting genes over-represented in the whale skin related to cell epithelial development, regulation of gene expression and cell maintenance . Comparison of gene expression between male and female samples revealed sex-specific differences in gray and humpback whales. A differential gene expression analysis identified several x-linked genes that have been previously identified and show gene expression differences in male and female cetaceans, such as ZFX, DDX3X and USP9X. Establishing baseline skin gene expression profiles for these three baleen whale species sampled off the Oregon coast provides a foundation for linking transcriptome variation with physiological condition and environment.

Sremba, Angela↗

Small Signaling Peptides in Sorghum bicolor : Integrating Phylogeny and Gene Expression to Characterize Roles in Stem Development

Small signaling peptides (SSPs) are important regulators of plant growth, development, and responses to biotic and abiotic stress, yet their role in the C4 grass Sorghum bicolor is largely uncharacterized. To help fill this knowledge gap, 219 sorghum genes that encode SSPs were identified based on SSP sequences previously identified in Arabidopsis thaliana, Zea mays, Oryza sativa, Triticum aestivum , and Brachypodium distachyon . The 219 sorghum SSP-encoding genes were assigned to 19 gene families, analyzed for the presence of motifs, and aligned with genes that encode SSPs in other plants using phylogenetic analysis. Sorghum genes in 12 of the 19 SSP gene families had not been previously characterized. Expression of the 219 SSP-encoding genes in sorghum organs, during stem development, and in stem tissues and cell types revealed distinct spatial, temporal, and developmental patterns of expression. Genes associated with the SbCEP and SbRGF families were preferentially expressed in roots, whereas SbEPF genes were expressed in stem epidermal and pith parenchyma cells and panicles. The expression of genes during bioenergy sorghum stem growth and development was investigated because stems account for ~80% of harvested biomass and serve as conduits for water and nutrient transport between leaves and roots. During stem development, 28 SSP genes in several families ( CLE, EPF, CEP, GASS, PSY, ES, PSK, CAPE, POE ) were expressed at higher levels in zones of cell proliferation. For example, the TDIF homologs SbCLE41 and SbCLE42 were expressed at high levels in nascent stem nodes where they may regulate vascular bundle cambial activity and cell differentiation. A different set of 15 genes in the CIF, POE, CAPE, PSY, CEP, RALF , and CLE families were expressed at higher levels in zones of stem tissue differentiation highlighted by elevated expression of five SbRALFR s in the stem nodal plexus. Cell type–specific expression of many sorghum genes that encode SSPs was observed in fully elongated internodes indicating gene expression is regulated with high spatial resolution. Overall, the results provide a foundation of information for analysis of SSP function in sorghum that can be integrated with knowledge of sorghum gene regulatory networks to modulate traits important for production of sorghum crops.

bioenergy sorghum↗

Novel cis-regulatory elements as synthetic promoters to drive recombinant protein expression from the Chlamydomonas reinhardtii nuclear genome

Eukaryotic green microalgae represent a sustainable, photosynthetic biotechnology platform for generating high-value products. The model green alga Chlamydomonas reinhardtii has already been used to generate high value bioproducts such as recombinant proteins and terpenoids. However, low, unstable, and variable nuclear transgene expression has limited the ease and speed of metabolic engineering and recombinant protein expression in this system. Here, novel genetic devices for transgene expression in C. reinhardtii have been developed by identifying cis-regulatory DNA elements capable of driving high transgene expression in C. reinhardtii promoters using de novo motif discovery informatics approaches. Thirteen putative motifs were synthesized as concatemers, linked to a common minimal basal promoter, and assayed for their activity to drive expression of a yellow fluorescent protein reporter gene. Following transformation of the vectors into C. reinhardtii by electroporation, in vivo measurements of yellow fluorescent protein expression by flow cytometry revealed that five of the DNA motifs analyzed displayed significantly higher reporter expression compared to the basal promoter control. Two of the concatemerized motifs, despite being much smaller minimal cis-regulatory elements, drove reporter expression at levels approaching that of the conventionally-used AR1 promoter. This analysis provides insight into C. reinhardtii promoter structure and gene regulation, and provides a new toolbox of cis-regulatory elements that can be used to drive transgene expression at a variety of expression levels.

59 BASIC BIOLOGICAL SCIENCES↗

Evolutionary Analyses of Gene Expression Divergence in Panicum hallii : Exploring Constitutive and Plastic Responses Using Reciprocal Transplants

Abstract The evolution of gene expression is thought to be an important mechanism of local adaptation and ecological speciation. Gene expression divergence occurs through the evolution of cis- polymorphisms and through more widespread effects driven by trans-regulatory factors. Here, we explore expression and sequence divergence in a large sample of Panicum hallii accessions encompassing the species range using a reciprocal transplantation experiment. We observed widespread genotype and transplant site drivers of expression divergence, with a limited number of genes exhibiting genotype-by-site interactions. We used a modified FST–QST outlier approach (QPC analysis) to detect local adaptation. We identified 514 genes with constitutive expression divergence above and beyond the levels expected under neutral processes. However, no plastic expression responses met our multiple testing correction as QPC outliers. Constitutive QPC outlier genes were involved in a number of developmental processes and responses to abiotic environments. Leveraging earlier expression quantitative trait loci results, we found a strong enrichment of expression divergence, including for QPC outliers, in genes previously identified with cis and cis–environment interactions but found no patterns related to trans-factors. Population genetic analyses detected elevated sequence divergence of promoters and coding sequence of constitutive expression outliers but little evidence for positive selection on these proteins. Our results are consistent with a hypothesis of cis-regulatory divergence as a primary driver of expression divergence in P. hallii.

3′ TagSeq↗