Results for “comparative transcriptomics” — Page 6

DOE OSTI2021

Effect of Cultivation Conditions on the Bioconversion of 4-Hydroxybenzoic Acid in Two White-Rot Fungi

White-rot fungi (WRF) are the most efficient organisms for lignin degradation in nature. A recent study has also shown the ability of WRF in funneling lignin-derived aromatics, including 4-hydroxybenzoic acid (4HBA) and vanillic acid, to central metabolism. However, it is unknown how the cultivation conditions affect the conversion of lignin-derived aromatic compounds. To address this, we performed multi-omic analyses and tracked the conversion of 4HBA in different cultivation conditions in two white-rot fungi: Trametes versicolor and Ceriporiopsis subvermispora. Specifically, we evaluated the effect of static and agitation cultivation conditions in the absence and the presence of antioxidants on fungal performance. We found that the comparison between static and agitation has the higher number of differentially expressed genes or protein abundances at the transcriptomic and proteomic levels, respectively, compared to the presence or the absence of antioxidants. Gene mapping on the enzymatic steps of the proposed 4HBA conversion pathway revealed that agitation causes upregulation of several oxidase genes such as aldehyde dehydrogenases, aldehyde oxidases, and hydroxylases in both fungi. Time-course experiments including 4HBA and cellobiose as carbon sources showed that T. versicolor exhibited higher 4HBA conversion levels in agitation (~110% more conversion, compared to static cultivation conditions) and in the presence of antioxidants (~10% more conversion, compared to the treatments without antioxidants). Interestingly, these trends were not found in C. subvermispora, suggesting different carbon metabolism in both species. Our findings are fundamental to better understand carbon sequestration from lignin in WRF, which will be key for future applications in the lignin valorization field.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

NASA NTRSDate not supplied

Translational Line of Sight: a multi-omics longitudinal study of the murine retinal response to spaceflight hazard analogs

The space environment includes unique hazards like radiation and microgravity which adversely affect physiology and behavior of humans and rodent models. To better characterize the retinal response to spaceflight, we assessed a multi-omics NASA GeneLab dataset where 6-month-old female mice were gamma irradiated and/or hindlimb unloaded for 21 days followed by whole transcriptome shotgun sequencing (RNA-Seq) and reduced representation bisulfite sequencing (RRBS) of retina samples collected at 7 days, 1 month or 4 months post-exposure. We compared time-matched epigenomic and transcriptomic retinal profiles revealing a total of 4,178 differentially methylated loci or regions, and 457 differentially expressed genes. Highest correlation in methylation differences was seen across different conditions at the same time point (e.g., between radiation exposure and hindlimb unloaded at 7 days). Biological processes related to nucleotide metabolism were enriched in all groups with activation at 1 month and suppression at 7 days and 4 months. Genes and processes related to Notch and Wnt signaling showed alterations 4 months post-exposure. Interestingly, Notch3 and Lrg1 showed differential patterns in the NASA Twins Study in-flight samples and in response to stressors in the murine retina in the current study. A total of 23 genes were both differentially methylated and expressed, including genes involved in retinal disease or cataract development (Crybb3, Fgfr1, Pitpnm3, Sipa1l3, Sox9) and inflammatory response (B4galt6, Ppm1a, Sphk1). To our knowledge, the current multi-omics analysis is the first multi-omics study to interrogate the epigenomic and transcriptomic impacts of radiation and hindlimb unloading on the retina in isolation and in combination. The results provide an insight into the retinal response to individual spaceflight hazard analogs and their interplay at different post-exposure stages and contributes towards a mechanistic understanding of spaceflight-induced vision impairment using ground-based models.

Prachi Kothiyal↗

NASA NTRSDate not supplied

A Multi-omics Longitudinal Study of the Murine Retinal Response to Chronic Low-dose Irradiation and/or Simulated Microgravity

The space environment includes unique hazards like radiation and microgravity which adversely affect physiology and behavior of humans and rodent models. To better characterize the retinal response to spaceflight, we assessed a multi-omics NASA GeneLab dataset where 6-month-old female mice were gamma irradiated and/or hindlimb unloaded for 21 days followed by whole transcriptome shotgun sequencing (RNA-Seq) and reduced representation bisulfite sequencing (RRBS) of retina samples collected at 7 days, 1 month or 4 months post-exposure. We compared time-matched epigenomic and transcriptomic retinal profiles revealing a total of 4,178 differentially methylated loci or regions, and 457 differentially expressed genes. Highest correlation in methylation differences was seen across different conditions at the same time point (e.g., between radiation exposure and hindlimb unloaded at 7 days). Biological processes related to nucleotide metabolism were enriched in all groups with activation at 1 month and suppression at 7 days and 4 months. Genes and processes related to Notch and Wnt signaling showed alterations 4 months post-exposure. Interestingly, Notch3 and Lrg1 showed differential patterns in the NASA Twins Study in-flight samples and in response to stressors in the murine retina in the current study. A total of 23 genes were both differentially methylated and expressed, including genes involved in retinal disease or cataract development (Crybb3, Fgfr1, Pitpnm3, Sipa1l3, Sox9) and inflammatory response (B4galt6, Ppm1a, Sphk1). To our knowledge, the current multi-omics analysis is the first multi-omics study to interrogate the epigenomic and transcriptomic impacts of radiation and hindlimb unloading on the retina in isolation and in combination. The results provide an insight into the retinal response to individual spaceflight hazard analogs and their interplay at different post-exposure stages and contributes towards a mechanistic understanding of spaceflight-induced vision impairment using ground-based models.

Prachi Kothiyal↗

DOE OSTI2024

Exploring genetic diversity, population structure, and subgenome differences in the allopolyploid Camelina sativa : implications for future breeding and research studies

Abstract Camelina (Camelina sativa), an allohexaploid species, is an emerging aviation biofuel crop that has been the focus of resurgent interest in recent decades. To guide future breeding and crop improvement efforts, the community requires a deeper comprehension of subgenome dominance, often noted in allopolyploid species, “alongside an understanding of the genetic diversity” and population structure of material present within breeding programs. We conducted population genetic analyses of a C. sativa diversity panel, leveraging a new genome, to estimate nucleotide diversity and population structure, and analyzed for patterns of subgenome expression dominance among different organs. Our analyses confirm that C. sativa has relatively low genetic diversity and show that the SG3 subgenome has substantially lower genetic diversity compared to the other two subgenomes. Despite the low genetic diversity, our analyses identified 13 distinct subpopulations including two distinct wild populations and others putatively representing founders in existing breeding populations. When analyzing for subgenome composition of long non-coding RNAs, which are known to play important roles in (a)biotic stress tolerance, we found that the SG3 subgenome contained significantly more lincRNAs compared to other subgenomes. Similarly, transcriptome analyses revealed that expression dominance of SG3 is not as strong as previously reported and may not be universal across all organ types. From a global analysis, SG3 “was only significant higher expressed” in flower, flower bud, and fruit organs, which is an important discovery given that the crop yield is associated with these organs. Collectively, these results will be valuable for guiding future breeding efforts in camelina.

Agriculture↗

DOE OSTI2024

Entire expressed peripheral blood transcriptome in pediatric severe malarial anemia

Abstract This study on severe malarial anemia (SMA: Hb < 6.0 g/dL), a leading global cause of childhood morbidity and mortality, compares the entire expressed whole blood host transcriptome between Kenyan children (3-48 mos.) with non-SMA (Hb ≥ 6.0 g/dL, n = 39) and SMA ( n = 18). Differential expression analyses reveal 1403 up-regulated and 279 down-regulated transcripts in SMA, signifying impairments in host inflammasome activation, cell death, and innate immune and cellular stress responses. Immune cell profiling shows decreased memory responses, antigen presentation, and immediate pathogen clearance, suggesting an immature/improperly regulated immune response in SMA. Module repertoire analysis of blood-specific gene signatures identifies up-regulation of erythroid genes, enhanced neutrophil activation, and impaired inflammatory responses in SMA. Enrichment analyses converge on disruptions in cellular homeostasis and regulatory pathways for the ubiquitin-proteasome system, autophagy, and heme metabolism. Pathway analyses highlight activation in response to hypoxic conditions [Hypoxia Inducible Factor (HIF)−1 target and Reactive Oxygen Species (ROS) signaling] as a central theme in SMA. These signaling pathways are also top-ranking in protein abundance measures and a Ugandan SMA cohort with available transcriptomic data. Targeted RNA-Seq validation shows strong concordance with our entire expressed transcriptome data. These findings identify key molecular themes in SMA pathogenesis, offering potential targets for new malaria therapies.

60 APPLIED LIFE SCIENCES↗

NASA NTRS2022

A Multi-omics Longitudinal Study of the Murine Retinal Response to Chronic Low-dose Irradiation and Simulated Microgravity

The space environment includes unique hazards like radiation and microgravity which can adversely affect biological systems. We assessed a multi-omics NASA GeneLab dataset where mice were hindlimb unloaded and/or gamma irradiated for 21 days followed by retinal analysis at 7 days, 1 month or 4 months post-exposure. We compared time-matched epigenomic and transcriptomic retinal profiles resulting in a total of 4,178 differentially methylated loci or regions, and 457 differentially expressed genes. Highest correlation in methylation difference was seen across different conditions at the same time point. Nucleotide metabolism biological processes were enriched in all groups with activation at 1 month and suppression at 7 days and 4 months. Genes and processes related to Notch and Wnt signaling showed alterations 4 months post-exposure. A total of 23 genes showed significant changes in methylation and expression compared to unexposed controls, including genes involved in retinal function and inflammatory response. This multi-omics analysis interrogates the epigenomic and transcriptomic impacts of radiation and hindlimb unloading on the retina in isolation and in combination and highlights important molecular mechanisms at different post-exposure stages.

Prachi Kothiyal↗

DOE OSTI2026

Integrative analysis of CAM photosynthesis reveals its impact on primary metabolism in Yucca

Crassulacean Acid Metabolism (CAM) is an adaptation that temporally separates carbon uptake at night from photosynthesis during the day. CAM has evolved repeatedly across vascular plants, as its emergence may depend on simple regulatory changes to deeply conserved metabolic pathways. Modern CAM research relies heavily on interpretation of transcriptomic data, though regulation occurs at multiple levels following transcription. Additionally, while most research to date has focused on a handful of genes and metabolites in the core CAM pathway, the co-option of conserved regulatory and functional genes is bound to have wide-ranging effects on other aspects of primary metabolism. In this study, we integrate transcriptomic, proteomic, and metabolomic data to compare primary metabolism between the CAM species Yucca aloifolia and closely related C 3 species, Y. filamentosa. We observe minimal correlation between protein abundance and mRNA expression, suggesting significant post-transcriptional regulation in CAM species. We also find evidence of shifts in gene expression and metabolite accumulation outside of the central CAM pathway, suggesting that the shift to CAM has cascading effects across primary metabolism, especially nitrogen metabolism. Our findings provide insights into the metabolic shifts associated with CAM evolution, and highlight the complexity of its regulation at multiple biological levels.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Mef2c regulates bone mass through Sost -dependent and -independent mechanisms

We report Mef2c is a transcription factor that mediates key cellular behaviors that promote endochondral ossification and bone formation. Previously, Mef2c has been shown to regulate Sost transcription via its osteocyte-specific enhancer, ECR5, and conditional deletions of Mef2 cfl/fl with either Col1-Cre or Dmp1-Cre produced generalized high bone mass (HBM) consistent with Van Buchem Disease phenotypes. However, Sost -/- ; Mef2 cfl/fl ; Dmp1-Cre mice produced a significantly higher bone mass phenotype that Sost-/- alone suggesting that Mef2c modulates bone mass through additional mechanisms, independent of Sost. To identify new Mef2c transcriptional targets important in bone metabolism, we profiled gene expression by single-cell RNA sequencing in subpopulations of cells isolated from Mef2 cfl/fl ; Dmp1-Cre and Mef2 cfl/fl ; Bglap-Cre femurs, both strains exhibiting similar high bone mass phenotypes. However, we found Mef2cfl/fl; Bglap-Cre to also display a growth plate defect characterized by an expansion of several osteoprogenitor subpopulations. Differential gene expression analysis identified a total of 96 up- and 2434 down- regulated genes in Mef2 cfl/fl ; Bglap-Cre and 176 up- and 1041 down- regulated genes in Mef2 cfl/fl ; Dmp1-Cre bone cell subpopulations compared to wildtype mice. Mef2c deletion affected the transcriptomes across several cell types including mesenchymal progenitors (MP), osteoprogenitors (OSP), osteoblast (OB), and osteocyte (OCY) subpopulations. Several energy metabolism genes such as Uqcrb, Ndufv2, Ndufs3, Ndufa13, Ndufb9, Ndufb5, Cox6a1, Cox5a, Atp5o, Atp5g2, Atp5b, Atp5 were significantly down regulated in Mef2c-deficient OBs and OCYs, in both strains. Binding motif analysis of promoter regions of differentially expressed genes identified Mef2c binding in Bone Sialoprotein (BSP/Ibsp), a gene known to cause increased trabecular BV/TV in the femurs of Ibsp -/- mice. Immunohistochemical analysis confirmed the absence of Ibsp protein in OBs and OCYs. These findings suggests that the HBM in Sost -/- ; Mef2 cfl/fl ; Dmp1-Cre is caused by a multitude of transcriptional changes in genes that regulate bone formation, two of which are Sost and Ibsp.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Prediction of plant complex traits via integration of multi-omics data

The formation of complex traits is the consequence of genotype and activities at multiple molecular levels. However, connecting genotypes and these activities to complex traits remains challenging. Here, we investigate whether integrating genomic, transcriptomic, and methylomic data can improve prediction for six Arabidopsis traits. We find that transcriptome- and methylome-based models have performances comparable to those of genome-based models. However, models built for flowering time using different omics data identify different benchmark genes. Nine additional genes identified as important for flowering time from our models are experimentally validated as regulating flowering. Gene contributions to flowering time prediction are accession-dependent and distinct genes contribute to trait prediction in different genotypes. Models integrating multi-omics data perform best and reveal known and additional gene interactions, extending knowledge about existing regulatory networks underlying flowering time determination. These results demonstrate the feasibility of revealing molecular mechanisms underlying complex traits through multi-omics data integration.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

High light and temperature reduce photosynthetic efficiency through different mechanisms in the C 4 model Setaria viridis

C 4 plants frequently experience high light and high temperature conditions in the field, which reduce growth and yield. However, the mechanisms underlying these stress responses in C 4 plants have been under-explored, especially the coordination between mesophyll (M) and bundle sheath (BS) cells. We investigated how the C4 model plant Setaria viridis responded to a four-hour high light or high temperature treatment at photosynthetic, transcriptomic, and ultrastructural levels. Although we observed a comparable reduction of photosynthetic efficiency in high light or high temperature treated leaves, detailed analysis of multi-level responses revealed important differences in key pathways and M/BS specificity responding to high light and high temperature. We provide a systematic analysis of high light and high temperature responses in S. viridis, reveal different acclimation strategies to these two stresses in C 4 plants, discover unique light/temperature responses in C 4 plants in comparison to C 3 plants, and identify potential targets to improve abiotic stress tolerance in C 4 crops.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Discovering Innovations in Stress Tolerance through Comparative Gene Regulatory Network Analysis and Cell-Type Specific Expression Maps (Final Technical Report with Cover Page)

Through this grant, we developed a comparative framework to elucidate the mechanisms behind variations in environmental stress responses among a diverse group of species within the Brassicaceae family. Our focus was on the differences in physiological and transcriptomic responses to abscisic acid (ABA), a hormone associated with water stress. We examined the differential growth responses of four Brassicaceae species, finding that most exhibited reduced root growth correlated with smaller meristem size. In contrast, Schrenkiella parvula showed accelerated growth due to increased root cell elongation. We employed RNA sequencing to analyze the transcriptional responses to ABA across these species, and innovative bioinformatics techniques were used to pinpoint biological pathways with significant divergence. Additionally, we utilized DAP-seq to map the gene regulatory networks associated with ABAresponsive transcription factors, revealing that variations in the regulation of growth hormone biosynthesis play a critical role in the distinct ABA effects on root growth among the species. This research sets a new standard for comparative physiology by integrating comparative genomics and transcriptomics to uncover pathway divergences.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRSDate not supplied

Oxygen Deficiency in Spaceflight & its Impact on Plants’ Adaptive Changes

The goal of this study was to investigate the effects of hypoxic conditions in spaceflight. The distribution of genes involved with hypoxia in Arabidopsis thaliana and Brassica rapa were analyzed with the results from past spaceflight experiments to evaluate genes for future studies. Transcriptomes data of two different spaceflight studies of Arabidopsis thaliana from the NASA GeneLab database, GLDS-7 and GLDS-17, were compared. DNA microarrays were utilized for transcription profiling to conduct these studies. For GLDS-7, the response in spaceflight was studied with approaches that collected gene expression data. Leaves, hypocotyls, and root tissues were compared to the whole plant. For GLDS-17, seedlings and undifferentiated cultured cells were placed in the Biological Research in Canisters (BRIC), specifically BRIC-16. The genes related to hypoxia in Arabidopsis thaliana from these two studies were compared to genes in Brassica rapa with the TOAST database to evaluate similarities. When transcriptomes were analyzed for GLDS-7 and 17, genes that were considered significant had p-values ≤ 0.05 and log fold change values ≤ -1 or ≥1. Sixteen genes fulfilled the criteria. The genes related to hypoxia were alcohol dehydrogenase, elongation factor, ethylene-responsive factor, GUS, heat-shock proteins, NAP, RAP2.12, and RD20. The genes most impacted by spaceflight were heat-shock proteins. These genes were compared with Brassica rapa through Arabidopsis Ensemble Orthology from the TOAST Database. Similarities were seen in alcohol dehydrogenase, elongation factor, ethylene-responsive factor, heat-shock proteins, NAP, and RAP2.12. Overall, transcription profiling indicates that plants’ survival in spaceflight is dependent on adaptive changes with gene expression. This study also indicates that there are similarities in gene expression between Arabidopsis thaliana and Brassica rapa with comparable gene expression. Future studies could include analyzing additional species to understand which genes could be modified to ensure better yield of space crops amid hypoxic conditions.

hypoxia↗

NASA NTRSDate not supplied

A Gene Expression and Histologic Approach to Study Production and Outflow of Cerebrospinal Fluid in Hindlimb Suspended Rats

INTRODUCTION: The Spaceflight Associated Neuro-ocular Syndrome (SANS) is thought to be associated with weightlessness-induced cephalad fluid shift, possibly associated with a chronic elevation of intracranial pressure (ICP) during long-duration ISS missions. Changes in cerebrospinal fluid (CSF) dynamics and cranial compliance might be involved in the ICP increase. It is not known whether CSF production and/or outflow are altered in microgravity, but changes at the molecular and cellular level in the structures that produce and regulate the transcellular and paracellular secretion and reabsorption of CSF may be relevant. In this study, we used the rat hindlimb suspension (HS) model to examine the relationship between intracranial pressure (ICP) and the cellular responses to the prolonged change in body posture elicited specifically in the choroid plexus (CP). This was evaluated by transcriptomics, histopathology and ultrastructure of the CP and arachnoid villi (AV). METHODS: ICP was measured by telemetry. The morphology, ultrastructure, and gene expression profile of the CP was examined using male 9-month-old Long Evans rats subjected to HS for 14 and 90 days. A subset of animals completing 90-day HS returned to normal posture for 14 and 90 additional days of recovery. All HS rats had age-matched cage controls maintained in normal posture. A group of animals was maintained in conditions of 1% CO2 throughout the entire protocol. The rat brains were carefully removed and preserved for various analyses, including transmission electron microscopy (TEM), immunohistochemical analysis of specific targets involved in CSF regulation, and RNA sequencing analysis of laser capture micro-dissected CP tissue from the lateral ventricles. SUMMARY OF RESULTS: The ICP record was limited, which impacted any conclusions derived from the results. However, the ICP of rats in the CO2-enriched atmosphere appeared increased compared to those in a normal air composition. HS of rats resulted in an altered transcriptomics profile in the choroid plexus compared to animals maintained at normal posture. This was observed during both the HS period and during normal posture recovery period following 90 days of suspension. The greatest number of differentially expressed genes was observed at 90 days HS. Elevated CO2 also led to a different transcriptomics profile in the CP. The histology and ultrastructure results should be considered preliminary due to the limited number of samples. Examination of the CP by TEM showed blood vessel congestion and microvilli swelling in the CP, as well as some subpial and periventricular gliosis, but no associations were observed with treatment. There was no evidence of a reduced clearance of b-amyloid in the periventricular, perivascular and subpial regions of the brain in HS animals versus those in normal posture. Immunohistochemical staining of aquaporin 4, showed the localization to the subpial region and ependyma with increased immunoreactivity in the brains of CO2 exposed rats. This work was supported by awards NNX15AW48G to S.Z. Animal tissue and live data was provided under a tissue/data sharing agreement with Dr. Charles Fuller, UC Davis.

S B Zanello↗

NASA NTRSDate not supplied

A Gene Expression and Histologic Approach to Study Production and Outflow of Cerebrospinal Fluid in Hindlimb Suspended RATS

INTRODUCTION: The Spaceflight Associated Neuro-ocular Syndrome (SANS) is thought to be associated with weightlessness-induced cephalad fluid shift, possibly associated with a chronic elevation of intracranial pressure (ICP) during long-duration ISS missions. Changes in cerebrospinal fluid (CSF) dynamics and cranial compliance might be involved in the ICP increase. It is not known whether CSF production and/or outflow are altered in microgravity, but changes at the molecular and cellular level in the structures that produce and regulate the transcellular and paracellular secretion and reabsorption of CSF may be relevant. In this study, we used the rat hindlimb suspension (HS) model to examine the relationship between intracranial pressure (ICP) and the cellular responses to the prolonged change in body posture elicited specifically in the choroid plexus (CP). This was evaluated by transcriptomics, histopathology and ultrastructure of the CP and arachnoid villi (AV). METHODS: ICP was measured by telemetry. The morphology, ultrastructure, and gene expression profile of the CP was examined using male 9-month-old Long Evans rats subjected to HS for 14 and 90 days. A subset of animals completing 90-day HS returned to normal posture for 14 and 90 additional days of recovery. All HS rats had age-matched cage controls maintained in normal posture. A group of animals was maintained in conditions of 1% CO2 throughout the entire protocol. The rat brains were carefully removed and preserved for various analyses, including transmission electron microscopy (TEM), immunohistochemical analysis of specific targets involved in CSF regulation, and RNA sequencing analysis of laser capture micro-dissected CP tissue from the lateral ventricles. SUMMARY OF RESULTS: The ICP record was limited, which impacted any conclusions derived from the results. However, the ICP of rats in the CO2-enriched atmosphere appeared increased compared to those in a normal air composition. HS of rats resulted in an altered transcriptomics profile in the choroid plexus compared to animals maintained at normal posture. This was observed during both the HS period and during normal posture recovery period following 90 days of suspension. The greatest number of differentially expressed genes was observed at 90 days HS. Elevated CO2 also led to a different transcriptomics profile in the CP. The histology and ultrastructure results should be considered preliminary due to the limited number of samples. Examination of the CP by TEM showed blood vessel congestion and microvilli swelling in the CP, as well as some subpial and periventricular gliosis, but no associations were observed with treatment. There was no evidence of a reduced clearance of b-amyloid in the periventricular, perivascular and subpial regions of the brain in HS animals versus those in normal posture. Immunohistochemical staining of aquaporin 4, showed the localization to the subpial region and ependyma with increased immunoreactivity in the brains of CO2 exposed rats. This work was supported by awards NNX15AW48G to S.Z. Animal tissue and live data was provided under a tissue/data sharing agreement with Dr. Charles Fuller, UC Davis.

S B Zanello↗

DOE OSTI2022

A transcriptome software comparison for the analyses of treatments expected to give subtle gene expression responses

Background: In this comparative study we evaluate the performance of four software tools: DNAstar-D (DESeq2), DNAstar-E (edgeR), CLC Genomics and Partek Flow for identification of differentially expressed genes (DEGs) using a transcriptome of E. coli. The RNA-seq data are from the effect of below-background radiation 5.5 nGy total dose (0.2nGy/hr) on E. coli grown shielded from natural radiation 655 m below ground in a pre-World War II steel vault. The gene expression response to three supplemented sources of radiation designed to mimic natural background, 1952 – 5720 nGy in total dose (71–208 nGy/hr), are compared to this “radiation-deprived” treatment. In addition, RNA-seq data of Caenorhabditis elegans nematode from similar radiation treatments was analyzed by three of the software packages. Results: In E. coli, the four software programs identified one of the supplementary sources of radiation (KCl) to evoke about 5 times more transcribed genes than the minus-radiation treatment (69–114 differentially expressed genes, DEGs), and so the rest of the analyses used this KCl vs “Minus” comparison. After imposing a 30-read minimum cutoff, one of the DNAStar options shared two of the three steps (mapping, normalization, and statistic) with Partek Flow (they both used median of ratios to normalize and the DESeq2 statistical package), and these two programs identified the highest number of DEGs in common with each other (53). In contrast, when the programs used different approaches in each of the three steps, between 31 and 40 DEGs were found in common. Regarding the extent of expression differences, three of the four programs gave high fold-change results (15–178 fold), but one (DNAstar’s DESeq2) resulted in more conservative fold-changes (1.5–3.5). In a parallel study comparing three qPCR commercial validation software programs, these programs also gave variable results as to which genes were significantly regulated. Similarly, the C. elegans analysis showed exaggerated fold-changes in CLC and DNAstar’s edgeR while DNAstar-D was more conservative. Conclusions: Regarding the extent of expression (fold-change), and considering the subtlety of the very low level radiation treatments, in E. coli three of the four programs gave what we consider exaggerated fold-change results (15 – 178 fold), but one (DNAstar’s DESeq2) gave more realistic fold-changes (1.5–3.5). When RT-qPCR validation comparisons to transcriptome results were carried out, they supported the more conservative DNAstar-D’s expression results. When another model organism’s (nematode) response to these radiation differences was similarly analyzed, DNAstar-D also resulted in the most conservative expression patterns. Therefore, we would propose DESeq2 (“DNAstar-D”) as an appropriate software tool for differential gene expression studies for treatments expected to give subtle transcriptome responses.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Transcriptional networks underlying a primary ovarian insufficiency disorder in alligators naturally exposed to EDCs

Interactions between the endocrine system and environmental contaminants are responsible for impairing reproductive development and function. Despite the taxonomic diversity of affected species and attendant complexity inherent to natural systems, the underlying signaling pathways and cellular consequences are mostly studied in lab models. To resolve the genetic and endocrine pathways that mediate affected ovarian function in organisms exposed to endocrine disrupting contaminants in their natural environments, we assessed broad-scale transcriptional and steroidogenic responses to exogenous gonadotropin stimulation in juvenile alligators (Alligator missippiensis) originating from a lake with well-documented pollution (Lake Apopka, FL) and a nearby reference site (Lake Woodruff, FL). Here, we found that individuals from Lake Apopka are characterized by hyperandrogenism and display hyper-sensitive transcriptional responses to gonadotropin stimulation when compared to individuals from Lake Woodruff. Site-specific transcriptomic divergence appears to be driven by wholly distinct subsets of transcriptional regulators, indicating alterations to fundamental genetic pathways governing ovarian function. Consistent with broad-scale transcriptional differences, ovaries of Lake Apopka alligators displayed impediments to folliculogenesis, with larger germinal beds and decreased numbers of late-stage follicles. After resolving the ovarian transcriptome into clusters of co-expressed genes, most site-associated modules were correlated to ovarian follicule phenotypes across individuals. However, expression of two site-specific clusters were independent of ovarian cellular architecture and are hypothesized to represent alterations to cell-autonomous transcriptional programs. Collectively, our findings provide high resolution mapping of transcriptional patterns to specific reproductive function and advance our mechanistic understanding regarding impaired reproductive health in an established model of environmental endocrine disruption.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Phenotypic and functional characterization of corneal endothelial cells during in vitro expansion

The advent of cell culture-based methods for the establishment and expansion of human corneal endothelial cells (CEnC) has provided a source of transplantable corneal endothelium, with a significant potential to challenge the one donor-one recipient paradigm. However, concerns over cell identity remain, and a comprehensive characterization of the cultured CEnC across serial passages has not been performed. To this end, we compared two established CEnC culture methods by assessing the transcriptomic changes that occur during in vitro expansion. In confluent monolayers, low mitogenic culture conditions preserved corneal endothelial cell state identity better than culture in high mitogenic conditions. Expansion by continuous passaging induced replicative cell senescence. Transcriptomic analysis of the senescent phenotype identified a cell senescence signature distinct for CEnC. We identified activation of both classic and new cell signaling pathways that may be targeted to prevent senescence, a significant barrier to realizing the potential clinical utility of in vitro expansion.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

The genome of Shorea leprosula (Dipterocarpaceae) highlights the ecological relevance of drought in aseasonal tropical rainforests

Hyperdiverse tropical rainforests, such as the aseasonal forests in Southeast Asia, are supported by high annual rainfall. Its canopy is dominated by the species-rich tree family of Dipterocarpaceae (Asian dipterocarps), which has both ecological (e.g., supports flora and fauna) and economical (e.g., timber production) importance. Recent ecological studies suggested that rare irregular drought events may be an environmental stress and signal for the tropical trees. We assembled the genome of a widespread but near threatened dipterocarp, Shorea leprosula, and analyzed the transcriptome sequences of ten dipterocarp species representing seven genera. Comparative genomic and molecular dating analyses suggested a whole-genome duplication close to the Cretaceous-Paleogene extinction event followed by the diversification of major dipterocarp lineages (i.e. Dipterocarpoideae). Interestingly, the retained duplicated genes were enriched for genes upregulated by no-irrigation treatment. These findings provide molecular support for the relevance of drought for tropical trees despite the lack of an annual dry season.

59 BASIC BIOLOGICAL SCIENCES↗