Search NASA⌕ Search

SEARCH · Search NASA

Results for “gene expression analysis”

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 577 records · Page 32

Quantitative trait locus (QTL) mapping and transcriptome profiling identify QTLs and candidate genes associated with heat stress response during reproductive development in Camelina sativa

Camelina sativa (L.) Crantz is a low-input oilseed crop that has great potential in providing sustainable feedstock for biofuels and bioproducts. Climate change is threatening production of camelina with rising global temperatures. Elucidating the genetic response to high temperatures is essential for successful breeding of heat-tolerant camelina varieties. Here, we report a combinatorial approach to identifying candidate genes associated with heat stress by quantitative trait locus (QTL) mapping and comparative transcriptome profiling. A population of recombinant inbred lines (RILs) was grown in a controlled growth chamber under the high-temperature regimes for 14 days beginning at the onset of the reproductive stage. Several traits related to seed production were evaluated at maturity. The QTL analysis identified several regions with co-located traits on chromosomes 8, 10, and 12. Two RILs with contrasting phenotypic responses to heat stress were chosen for gene expression profiling via RNA sequencing. Multiple pathways and genes were found to be strongly affected by heat stress, and many genes expressed differently between the two RILs. Several genes identified within the QTL regions were considered strong candidates that may control heat tolerance during reproduction in camelina. These studies provide resources for future studies that may assist in improving the heat tolerance of camelina.

60 APPLIED LIFE SCIENCES↗

CRISPR/Cas9 Directed Reprogramming of iPSC for Accelerated Motor Neuron Differentiation Leads to Dysregulation of Neuronal Fate Patterning and Function

Neurodegeneration causes a significant disease burden and there are few therapeutic interventions available for reversing or slowing the disease progression. Induced pluripotent stem cells (iPSCs) hold significant potential since they are sourced from adult tissue and have the capacity to be differentiated into numerous cell lineages, including motor neurons. This differentiation process traditionally relies on cell lineage patterning factors to be supplied in the differentiation media. Genetic engineering of iPSC with the introduction of recombinant master regulators of motor neuron (MN) differentiation has the potential to shorten and streamline cell developmental programs. We have established stable iPSC cell lines with transient induction of exogenous LHX3 and ISL1 from the Tet-activator regulatory region and have demonstrated that induction of the transgenes is not sufficient for the development of mature MNs in the absence of neuron patterning factors. Comparative global transcriptome analysis of MN development from native and Lhx-ISL1 modified iPSC cultures demonstrated that the genetic manipulation helped to streamline the neuronal patterning process. However, leaky gene expression of the exogenous MN master regulators in iPSC resulted in the premature activation of genetic pathways characteristic of the mature MN function. Dysregulation of metabolic and regulatory pathways within the developmental process affected the MN electrophysiological responses.

59 BASIC BIOLOGICAL SCIENCES↗

Site-specific acetylation of polynucleotide kinase 3'-phosphatase regulates its distinct role in DNA repair pathways

Mammalian polynucleotide kinase 3'-phosphatase (PNKP), a DNA end-processing enzyme with 3'-phosphatase and 5'-kinase activities, is involved in multiple DNA repair pathways, including base excision (BER), single-strand break (SSBR), and double-strand break repair (DSBR). However, little is known as to how PNKP functions in such diverse repair processes. Here we report that PNKP is acetylated at K142 (AcK142) by p300 constitutively but at K226 (AcK226) by CBP, only after DSB induction. Co-immunoprecipitation analysis using AcK142 or AcK226 PNKP-specific antibodies showed that AcK142-PNKP associates only with BER/SSBR, and AcK226 PNKP with DSBR proteins. Despite the modest effect of acetylation on PNKP’s enzymatic activity in vitro, cells expressing non-acetylable PNKP (K142R or K226R) accumulated DNA damage in transcribed genes. Intriguingly, in striatal neuronal cells of a Huntington's Disease (HD)-based mouse model, K142, but not K226, was acetylated. This is consistent with the reported degradation of CBP, but not p300, in HD cells. Moreover, transcribed genomes of HD cells progressively accumulated DSBs. Chromatin-immunoprecipitation analysis demonstrated the association of Ac-PNKP with the transcribed genes, consistent with PNKP’s role in transcription-coupled repair. Thus, our findings demonstrate that acetylation at two lysine residues, located in different domains of PNKP, regulates its distinct role in BER/SSBR versus DSBR.

59 BASIC BIOLOGICAL SCIENCES↗

Systems Analysis of NADH Dehydrogenase Mutants Reveals Flexibility and Limits of Pseudomonas taiwanensis VLB120’s Metabolism

Obligate aerobic organisms rely on a functional electron transport chain for energy conservation and NADH oxidation. Because of this essential requirement, the genes of this pathway are likely constitutively and highly expressed to avoid a cofactor imbalance and energy shortage under fluctuating environmental conditions. In this work, we investigated the essentiality of the three NADH dehydrogenases of the respiratory chain of the obligate aerobe Pseudomonas taiwanensis VLB120 and the impact of the knockouts of corresponding genes on its physiology and metabolism. While a mutant lacking all three NADH dehydrogenases seemed to be nonviable, the single or double knockout mutant strains displayed no, or only a weak, phenotype. Only the mutant deficient in both type 2 dehydrogenases showed a clear phenotype with biphasic growth behavior and a strongly reduced growth rate in the second phase. In-depth analyses of the metabolism of the generated mutants, including quantitative physiological experiments, transcript analysis, proteomics, and enzyme activity assays revealed distinct responses to type 2 and type 1 dehydrogenase deletions. An overall high metabolic flexibility enables P. taiwanensis to cope with the introduced genetic perturbations and maintain stable phenotypes, likely by rerouting of metabolic fluxes. This metabolic adaptability has implications for biotechnological applications. While the phenotypic robustness is favorable in large-scale applications with inhomogeneous conditions, the possible versatile redirecting of carbon fluxes upon genetic interventions can thwart metabolic engineering efforts. While Pseudomonas has the capability for high metabolic activity and the provision of reduced redox cofactors important for biocatalytic applications, exploitation of this characteristic might be hindered by high, constitutive activity of and, consequently, competition with the NADH dehydrogenases of the respiratory chain. The in-depth analysis of NADH dehydrogenase mutants of Pseudomonas taiwanensis VLB120 presented here provides insight into the phenotypic and metabolic response of this strain to these redox metabolism perturbations. This high degree of metabolic flexibility needs to be taken into account for rational engineering of this promising biotechnological workhorse toward a host with a controlled and efficient supply of redox cofactors for product synthesis.

36 MATERIALS SCIENCE↗

Dual reporter transgene driven by 2.3Col1a1 promoter is active in differentiated osteoblasts

AIM: As quantitative and spatial analyses of promoter reporter constructs are not easily performed in intact bone, we designed a reporter gene specific to bone, which could be analyzed both visually and quantitatively by using chloramphenicol acetyltransferase (CAT) and a cyan version of green fluorescent protein (GFPcyan), driven by a 2.3-kb fragment of the rat collagen promoter (Col2.3). METHODS: The construct Col2.3CATiresGFPcyan was used for generating transgenic mice. Quantitative measurement of promoter activity was performed by CAT analysis of different tissues derived from transgenic animals; localization was performed by visualized GFP in frozen bone sections. To assess transgene expression during in vitro differentiation, marrow stromal cell and neonatal calvarial osteoblast cultures were analyzed for CAT and GFP activity. RESULTS: In mice, CAT activity was detected in the calvaria, long bone, teeth, and tendon, whereas histology showed that GFP expression was limited to osteoblasts and osteocytes. In cell culture, increased activity of CAT correlated with increased differentiation, and GFP activity was restricted to mineralized nodules. CONCLUSION: The concept of a dual reporter allows a simultaneous visual and quantitative analysis of transgene activity in bone.

Non-NASA Center↗

Identification of candidate host-specificity genes in Exserohilum turcicum using comparative genomics and transcriptomics

Abstract Exserohilum turcicum causes northern corn leaf blight and sorghum leaf blight. While the same species cause disease in both crops, the strains are host-specific. Here, we report the sequence and de novo annotated assemblies of one sorghum- and one maize-specific E. turcicum strain. The strains were sequenced using the PacBio Sequel II system. The total genome length for both assemblies was between 44 and 45 Mb with N50 of ∼2.5 Mb. Ninety-eight percent of the Benchmarking Universal Single-Copy Orthologs (BUSCO) for both assemblies had complete status. The estimated number of genes was 11,762 and 12,029 in the sorghum- and maize-specific isolates, respectively. Funannotate, EffectorP, SignalP, and transcriptome data were used to create functional annotation of each genome. The whole-genome comparison identified ten large-scale inversions and three translocations between the maize- and sorghum-specific strains, along with homologous genes and gene duplications. RNA was sequenced from the maize- and sorghum-specific isolate 10 days post-inoculation in maize and sorghum and from axenic cultures. Gene expression data from planta and axenic growth experiments were compared for each strain. Candidate host-specificity genes were identified by combining results from whole-genome comparison, synteny analysis, gene annotations, and transcriptome data. Overall, this study identified several candidate host-specificity genes that provide insights into E. turcicum interaction with its hosts.

Krone, Mara J. (ORCID:0000000159006624)↗

COMPILE: a GWAS computational pipeline for gene discovery in complex genomes

Abstract Background Genome-Wide Association Studies (GWAS) are used to identify genes and alleles that contribute to quantitative traits in large and genetically diverse populations. However, traits with complex genetic architectures create an enormous computational load for discovery of candidate genes with acceptable statistical certainty. We developed a streamlined computational pipeline for GWAS (COMPILE) to accelerate identification and annotation of candidate maize genes associated with a quantitative trait, and then matches maize genes to their closest rice and Arabidopsis homologs by sequence similarity. Results COMPILE executed GWAS using a Mixed Linear Model that incorporated, without compression, recent advancements in population structure control, then linked significant Quantitative Trait Loci (QTL) to candidate genes and RNA regulatory elements contained in any genome. COMPILE was validated using published data to identify QTL associated with the traits of α-tocopherol biosynthesis and flowering time, and identified published candidate genes as well as additional genes and non-coding RNAs. We then applied COMPILE to 274 genotypes of the maize Goodman Association Panel to identify candidate loci contributing to resistance of maize stems to penetration by larvae of the European Corn Borer ( Ostrinia nubilalis ). Candidate genes included those that encode a gene of unknown function, WRKY and MYB-like transcriptional factors, receptor-kinase signaling, riboflavin synthesis, nucleotide-sugar interconversion, and prolyl hydroxylation. Expression of the gene of unknown function has been associated with pathogen stress in maize and in rice homologs closest in sequence identity. Conclusions The relative speed of data analysis using COMPILE allowed comparison of population size and compression. Limitations in population size and diversity are major constraints for a trait and are not overcome by increasing marker density. COMPILE is customizable and is readily adaptable for application to species with robust genomic and proteome databases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Disruption of histone acetylation homeostasis reveals multilayered chromatin regulation for transcriptional resiliency

Background Epigenetic modifications, nucleosome occupancy, and three-dimensional chromatin architecture collectively create a multi-layered, highly interactive regulatory system for controlling genomic functionality. Dysregulation of epigenetic processes leads to a plethora of abnormalities including disease states. Therapies focused on epigenetic modulation can alter gene expression to correct dysfunction, though the perpetuation of these states and the relationships among chromatin regulatory layers is not well understood. Results Here, we investigated global and local chromatin structural and functional responses after acute histone deacetylase inhibitor treatment (suberoylanilide hydroxamic acid) in lung cancer cells across time. Treatment substantially increased global histone acetylation resulting in a pervasive but not distinctive signature. The spread of acetylation did not significantly impact global chromatin accessibility, and nucleosome remodeling largely occurred at finer scales in functionally relevant genomic regions. Indeed, both H 3 K 4 trimethylation, a mark of active transcription, and gene expression changes were altered in a controlled locus-specific manner, suggesting aberrant acetylation indirectly leads to balanced and bidirectional gene expression profiles from tighter regulation of other chromatin features. HDACi treatment induced (13%) genomic rearrangement in chromatin compartmentalization and moderate weakening of topologically associating domains. Conclusions Continuous wavelet analysis of these features demonstrates that scale-dependent, locus-specific factors influence the relationship between chromatin architecture and functional output, suggesting that regulation of transcription and nucleosome remodeling is not entirely (nor linearly) dependent upon large scale compartment exchange. Structural and functional responses are most pronounced early after treatment with partial persistence of differential local chromatin features and expression later in time; this highlights the plasticity of chromatin regulation, which may have implications for the efficacy of epigenetic treatments. These results demonstrate the effectiveness of multi-layered regulation of transcription: in resilient systems, disruption of one chromatin feature does not distort the regulation of other features in supporting a transcriptional program that allows for survival.

59 BASIC BIOLOGICAL SCIENCES↗

Alicyclobacillus mali FL18 as a Novel Source of Glycosyl Hydrolases: Characterization of a New Thermophilic β-Xylosidase Tolerant to Monosaccharides

A thermo-acidophilic bacterium, Alicyclobacillus mali FL18, was isolated from a hot spring of Pisciarelli, near Naples, Italy; following genome analysis, a novel putative β-xylosidase, AmβXyl, belonging to the glycosyl hydrolase (GH) family 3 was identified. A synthetic gene was produced, cloned in pET-30a(+), and expressed in Escherichia coli BL21 (DE3) RIL. The purified recombinant protein, which showed a dimeric structure, had optimal catalytic activity at 80 °C and pH 5.6, exhibiting 60% of its activity after 2 h at 50 °C and displaying high stability (more than 80%) at pH 5.0–8.0 after 16 h. AmβXyl is mainly active on both para-nitrophenyl-β-D-xylopyranoside (K M 0.52 mM, k cat 1606 s -1 , and k cat /K M 3088.46 mM -1 ·s -1 ) and para-nitrophenyl-α-L-arabinofuranoside (K M 10.56 mM, k cat 2395.8 s -1 , and k cat /K M 226.87 mM -1 ·s -1 ). Thin-layer chromatography showed its ability to convert xylooligomers (xylobiose and xylotriose) into xylose, confirming that AmβXyl is a true β-xylosidase. Furthermore, no inhibitory effect on enzymatic activity by metal ions, detergents, or EDTA was observed except for 5 mM Cu 2+ . AmβXyl showed an excellent tolerance to organic solvents; in particular, the enzyme increased its activity at high concentrations (30%) of organic solvents such as ethanol, methanol, and DMSO. Lastly, the enzyme showed not only a good tolerance to inhibition by xylose, arabinose, and glucose, but was activated by 0.75 M xylose and up to 1.5 M by both arabinose and glucose. The high tolerance to organic solvents and monosaccharides together with other characteristics reported above suggests that AmβXyl may have several applications in many industrial fields.

59 BASIC BIOLOGICAL SCIENCES↗

Biological Insights at the Interface of Multiple Arabidopsis Legacy Datasets

The NASA GeneLab database includes an open-access collection of datasets yielded by space biology experiments. Six Gene Lab Data Sets (GLDS’s) performed in Arabidopsis were selected for analysis (7/17/44/121/205/213), all of which included transcriptome data from spaceflight and ground control environments. Hardware, ecotype, environmental conditions, and other experimental conditions varied, allowing the observations of overarching gene expression impacts of microgravity on plant life without focusing on effects of specific experimental conditions. Using GeneLab pre-processed datasets as the basis for the study, RNA microarray data were analyzed to identify genes that showed altered expression in microgravity when compared to control samples for each individual GLDS. All differentially expressed genes were compared to locate differentially expressed genes common between spaceflight experiments. The most noteworthy result is that not one gene shared differential expression among the six GLDS’s. However, gene expression was not influenced randomly by the microgravity environment, as there were several gene ontology terms that were significantly enriched across all experiments. These included 20 significantly enriched biological processes, and although the genes which enriched each term varied, there were many cases of specific genes common to clusters of multiple GLDS’s. Gene expression such as NAC92 and ERF011 or membrane structural element FFP6 provide insight and direction toward understanding the plant response to spaceflight. Characterizing these common processes and the shared differentially expressed genes has demonstrated potential targets for further study to understand and modulate the biological response of plants in microgravity. Life on Earth has never been subjected to the absence of gravity as a selective pressure, so observing how life forms react to a microgravity environment could provide insight to our shared fundamental biological processes. It is also feasible that the genetic modification of specific genes linked to the microgravity response could improve health and yield of space crops.

Joseph Emhof↗

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↗

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↗

Two cis elements collaborate to spatially repress transcription from a sea urchin promoter

The expression pattern of many territory-specific genes in metazoan embryos is maintained by an active process of negative spatial regulation. However, the mechanism of this strategy of gene regulation is not well understood in any system. Here we show that reporter constructs containing regulatory sequence for the SM30-alpha gene of Stronglyocentrotus purpuratus are expressed in a pattern congruent with that of the endogenous SM30 gene(s), largely as a result of active transcriptional repression in cell lineages in which the gene is not normally expressed. Chloramphenicol acetyl transferase assays of deletion constructs from the 2600-bp upstream region showed that repressive elements were present in the region from -1628 to -300. In situ hybridization analysis showed that the spatial fidelity of expression was severely compromised when the region from -1628 to -300 was deleted. Two highly repetitive sequence motifs, (G/A/C)CCCCT and (T/C)(T/A/C)CTTTT(T/A/C), are present in the -1628 to -300 region. Representatives of these elements were analyzed by gel mobility shift experiments and were found to interact specifically with protein in crude nuclear extracts. When oligonucleotides containing either sequence element were co-injected with a correctly regulated reporter as potential competitors, the reporter was expressed in inappropriate cells. When composite oligonucleotides, containing both sequence elements, were fused to a misregulated reporter, the expression of the reporter in inappropriate cells was suppressed. Comparison of composite oligonucleotides with oligonucleotides containing single constituent elements show that both sequence elements are required for effective spatial regulation. Thus, both individual elements are required, but only a composite element containing both elements is sufficient to function as a tissue-specific repressive element.

Non-NASA Center↗

A multi-plex protein expression system for production of complex enzyme formulations in Trichoderma reesei

Abstract Heterologous protein production has been challenging in the hyper-cellulolytic fungus, Trichoderma reesei as the species is known for poor transformation efficiency, low homologous recombination frequency, and marginal screening systems for the identification of successful transformants. We have applied the 2A-peptide multi-gene expression system to co-express four proteins, which include three cellulases: a cellobiohydrolase (CBH1), an endoglucanase (EG1), and a β-D-glucosidase (BGL1), as well as the enhanced green fluorescent protein (eGFP) marker protein. We designed a new chassis vector, pTrEno-4X-2A, for this work. Expression of these cellulase enzymes was confirmed by real-time quantitative reverse transcription PCR and immunoblot analysis. The activity of each cellulase was assessed using chromogenic substrates, which confirmed the functionality of the enzymes. Expression and activity of these enzymes were proportional to the level of eGFP fluorescence, thereby validating the reliability of this screening technique. An 18-fold differencein protein expression was observed between the first and third genes within the 2A-peptide construct. The availability of this new multi-gene expression and screening tool is expected to greatly impact multi-enzyme applications, such as the production of complex commercial enzyme formulations and metabolic pathway enzymes, especially those destined for cell-free applications.

09 BIOMASS FUELS↗

Biallelic variants in SNUPN cause a limb girdle muscular dystrophy with myofibrillar-like features

Abstract Alterations in RNA-splicing are a molecular hallmark of several neurological diseases, including muscular dystrophies, where mutations in genes involved in RNA metabolism or characterized by alterations in RNA splicing have been described. Here, we present five patients from two unrelated families with a limb-girdle muscular dystrophy (LGMD) phenotype carrying a biallelic variant in SNUPN gene. Snurportin-1, the protein encoded by SNUPN, plays an important role in the nuclear transport of small nuclear ribonucleoproteins (snRNPs), essential components of the spliceosome. We combine deep phenotyping, including clinical features, histopathology and muscle MRI, with functional studies in patient-derived cells and muscle biopsies to demonstrate that variants in SNUPN are the cause of a new type of LGMD according to current definition. Moreover, an in vivo model in Drosophila melanogaster further supports the relevance of Snurportin-1 in muscle. SNUPN patients show a similar phenotype characterized by proximal weakness starting in childhood, restrictive respiratory dysfunction and prominent contractures, although inter-individual variability in terms of severity even in individuals from the same family was found. Muscle biopsy showed myofibrillar-like features consisting of myotilin deposits and Z-disc disorganization. MRI showed predominant impairment of paravertebral, vasti, sartorius, gracilis, peroneal and medial gastrocnemius muscles. Conservation and structural analyses of Snurportin-1 p.Ile309Ser variant suggest an effect in nuclear-cytosol snRNP trafficking. In patient-derived fibroblasts and muscle, cytoplasmic accumulation of snRNP components is observed, while total expression of Snurportin-1 and snRNPs remains unchanged, which demonstrates a functional impact of SNUPN variant in snRNP metabolism. Furthermore, RNA-splicing analysis in patients’ muscle showed widespread splicing deregulation, in particular in genes relevant for muscle development and splicing factors that participate in the early steps of spliceosome assembly. In conclusion, we report that SNUPN variants are a new cause of limb girdle muscular dystrophy with specific clinical, histopathological and imaging features, supporting SNUPN as a new gene to be included in genetic testing of myopathies. These results further support the relevance of splicing-related proteins in muscle disorders.

Iruzubieta, Pablo (ORCID:0000000303316222)↗

Information-incorporated gene network construction with FDR control

Abstract Motivation Large-scale gene expression studies allow gene network construction to uncover associations among genes. To study direct associations among genes, partial correlation-based networks are preferred over marginal correlations. However, FDR control for partial correlation-based network construction is not well-studied. In addition, currently available partial correlation-based methods cannot take existing biological knowledge to help network construction while controlling FDR. Results In this paper, we propose a method called Partial Correlation Graph with Information Incorporation (PCGII). PCGII estimates partial correlations between each pair of genes by regularized node-wise regression that can incorporate prior knowledge while controlling the effects of all other genes. It handles high-dimensional data where the number of genes can be much larger than the sample size and controls FDR at the same time. We compare PCGII with several existing approaches through extensive simulation studies and demonstrate that PCGII has better FDR control and higher power. We apply PCGII to a plant gene expression dataset where it recovers confirmed regulatory relationships and a hub node, as well as several direct associations that shed light on potential functional relationships in the system. We also introduce a method to supplement observed data with a pseudogene to apply PCGII when no prior information is available, which also allows checking FDR control and power for real data analysis. Availability and implementation R package is freely available for download at https://cran.r-project.org/package=PCGII.

59 BASIC BIOLOGICAL SCIENCES↗

Repeated exposure to eucalyptus wood smoke alters pulmonary gene and metabolic profiles in male Long-Evans rats

Abstract Exposure to wildfire smoke is associated with both acute and chronic cardiopulmonary illnesses, which are of special concern for wildland firefighters who experience repeated exposure to wood smoke. It is necessary to better understand the underlying pathophysiology by which wood smoke exposure increases pulmonary disease burdens in this population. We hypothesize that wood smoke exposure produces pulmonary dysfunction, lung inflammation, and gene expression profiles associated with future pulmonary complications. Male Long-Evans rats were intermittently exposed to smoldering eucalyptus wood smoke at 2 concentrations, low (11.0 ± 1.89 mg/m3) and high (23.7 ± 0.077 mg/m3), over a 2-week period. Whole-body plethysmography was measured intermittently throughout. Lung tissue and lavage fluid were collected 24 h after the final exposure for transcriptomics and metabolomics. Increasing smoke exposure upregulated neutrophils and select cytokines in the bronchoalveolar lavage fluid. In total, 3446 genes were differentially expressed in the lungs of rats in the high smoke exposure and only 1 gene in the low smoke exposure (Cd151). Genes altered in the high smoke group reflected changes to the Eukaryotic Initiation Factor 2 stress and oxidative stress responses, which mirrored metabolomics analyses. xMWAS-integrated analysis revealed that smoke exposure significantly altered pathways associated with oxidative stress, lung morphogenesis, and tumor proliferation pathways. These results indicate that intermittent, 2-week exposure to eucalyptus wood smoke leads to transcriptomic and metabolic changes in the lung that may predict future lung disease development. Collectively, these findings provide insight into cellular signaling pathways that may contribute to the chronic pulmonary conditions observed in wildland firefighters.

Toxicology↗

Integrated Proteomic and Glycoproteomic Characterization of Human High-Grade Serous Ovarian Carcinoma

Many gene products exhibit great structural heterogeneity because of an array of modifications. These modifications are not directly encoded in the genomic template but often affect the functionality of proteins. Protein glycosylation plays a vital role in proper protein functions. However, the analysis of glycoproteins has been challenging compared with other protein modifications, such as phosphorylation. Here, we perform an integrated proteomic and glycoproteomic analysis of 83 prospectively collected high-grade serous ovarian carcinoma (HGSC) and 23 non-tumor tissues. Integration of the expression data from global proteomics and glycoproteomics reveals tumor-specific glycosylation, uncovers different glycosylation associated with three tumor clusters, and identifies glycosylation enzymes that were correlated with the altered glycosylation. In addition to providing a valuable resource, these results provide insights into the potentialroles of glycosylation in the pathogenesis of HGSC, with the possibility of distinguishing pathological outcomes of ovarian tumors from non-tumors, as well as classifying tumor clusters.

Hu, Yingwei↗