Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Protein sequence 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 361 records · Page 20

Prochlorococcus phage ferredoxin: structural characterization and electron transfer to cyanobacterial sulfite reductases

Marine cyanobacteria are infected by phages whose genomes encode ferredoxin (Fd) electron carriers. These Fds are thought to redirect the energy harvested from light to phage-encoded oxidoreductases that enhance viral fitness, but it is unclear how the biophysical properties and partner specificities of phage Fds relate to those of photosynthetic organisms. Here, results of a bioinformatics analysis using a sequence similarity network revealed that phage Fds are most closely related to cyanobacterial Fds that transfer electrons from photosystems to oxidoreductases involved in nutrient assimilation. Structural analysis of myovirus P-SSM2 Fd (pssm2-Fd), which infects the cyanobacterium Prochlorococcus marinus, revealed high levels of similarity to cyanobacterial Fds (root mean square deviations of ≤0.5 Å). Additionally, pssm2-Fd exhibited a low midpoint reduction potential (–336 mV versus a standard hydrogen electrode), similar to other photosynthetic Fds, although it had lower thermostability (Tm = 28 °C) than did many other Fds. When expressed in an Escherichia coli strain deficient in sulfite assimilation, pssm2-Fd complemented bacterial growth when coexpressed with a P. marinus sulfite reductase, revealing that pssm2-Fd can transfer electrons to a host protein involved in nutrient assimilation. The high levels of structural similarity with cyanobacterial Fds and reactivity with a host sulfite reductase suggest that phage Fds evolved to transfer electrons to cyanobacterially encoded oxidoreductases.

cyanobacteria↗

Characterization of Heparin’s Conformational Ensemble by Molecular Dynamics Simulations and Nuclear Magnetic Resonance Spectroscopy

Heparin is a highly charged, polysulfated polysaccharide and serves as an anticoagulant. Heparin binds to multiple proteins throughout the body, suggesting a large range of potential therapeutic applications. Although its function has been characterized in multiple physiological contexts, heparin’s solution conformational dynamics and structure–function relationships are not fully understood. Molecular dynamics (MD) simulations facilitate the analysis of a molecule’s underlying conformational ensemble, which then provides important information necessary for understanding structure–function relationships. However, for MD simulations to afford meaningful results, they must both provide adequate sampling and accurately represent the energy properties of a molecule. The aim of this study is to compare heparin’s conformational ensemble using two well-developed force fields for carbohydrates, known as GLYCAM06 and CHARMM36, using replica exchange molecular dynamics (REMD) simulations, and to validate these results with NMR experiments. In this work, the anticoagulant sequence, an ultra-low-molecular-weight heparin, known as Arixtra (fondaparinux, sodium), was simulated with both parameter sets. The results suggest that GLYCAM06 matches experimental nuclear magnetic resonance three-bond J-coupling values measured for Arixtra better than CHARMM36. In addition, NOESY and ROESY experiments suggest that Arixtra is very flexible in the sub-millisecond time scale and does not adopt a unique structure at 25 C. Moreover, GLYCAM06 affords a much more dynamic conformational ensemble for Arixtra than CHARMM36.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ARG1 (altered response to gravity) encodes a DnaJ-like protein that potentially interacts with the cytoskeleton

Gravitropism allows plant organs to direct their growth at a specific angle from the gravity vector, promoting upward growth for shoots and downward growth for roots. Little is known about the mechanisms underlying gravitropic signal transduction. We found that mutations in the ARG1 locus of Arabidopsis thaliana alter root and hypocotyl gravitropism without affecting phototropism, root growth responses to phytohormones or inhibitors of auxin transport, or starch accumulation. The positional cloning of ARG1 revealed a DnaJ-like protein containing a coiled-coil region homologous to coiled coils found in cytoskeleton-interacting proteins. These data suggest that ARG1 participates in a gravity-signaling process involving the cytoskeleton. A combination of Northern blot studies and analysis of ARG1-GUS fusion-reporter expression in transgenic plants demonstrated that ARG1 is expressed in all organs. Ubiquitous ARG1 expression in Arabidopsis and the identification of an ortholog in Caenorhabditis elegans suggest that ARG1 is involved in other essential processes.

Non-NASA Center↗

Microarray analysis identifies Salmonella genes belonging to the low-shear modeled microgravity regulon

The low-shear environment of optimized rotation suspension culture allows both eukaryotic and prokaryotic cells to assume physiologically relevant phenotypes that have led to significant advances in fundamental investigations of medical and biological importance. This culture environment has also been used to model microgravity for ground-based studies regarding the impact of space flight on eukaryotic and prokaryotic physiology. We have previously demonstrated that low-shear modeled microgravity (LSMMG) under optimized rotation suspension culture is a novel environmental signal that regulates the virulence, stress resistance, and protein expression levels of Salmonella enterica serovar Typhimurium. However, the mechanisms used by the cells of any species, including Salmonella, to sense and respond to LSMMG and identities of the genes involved are unknown. In this study, we used DNA microarrays to elucidate the global transcriptional response of Salmonella to LSMMG. When compared with identical growth conditions under normal gravity (1 x g), LSMMG differentially regulated the expression of 163 genes distributed throughout the chromosome, representing functionally diverse groups including transcriptional regulators, virulence factors, lipopolysaccharide biosynthetic enzymes, iron-utilization enzymes, and proteins of unknown function. Many of the LSMMG-regulated genes were organized in clusters or operons. The microarray results were further validated by RT-PCR and phenotypic analyses, and they indicate that the ferric uptake regulator is involved in the LSMMG response. The results provide important insight about the Salmonella LSMMG response and could provide clues for the functioning of known Salmonella virulence systems or the identification of uncharacterized bacterial virulence strategies.

Non-NASA Center↗

Putting Humpty Dumpty Back Together Again: What Does Protein Quantification Mean in Bottom-Up Proteomics?

Bottom-up proteomics provides peptide measurements and has been invaluable for moving proteomics into large-scale analyses. Commonly, a single quantitative value is reported for each protein-coding gene by aggregating peptide quantities into protein groups following protein inference or parsimony. However, given the complexity of both RNA splicing and post-translational protein modification, it is overly simplistic to assume that all peptides that map to a singular protein-coding gene will demonstrate the same quantitative response. Here, by assuming that all peptides from a protein-coding sequence are representative of the same protein, we may miss the discovery of important biological differences. To capture the contributions of existing proteoforms, we need to reconsider the practice of aggregating protein values to a single quantity per protein-coding gene.

59 BASIC BIOLOGICAL SCIENCES↗

RNA Splicing Events in Circulation Distinguish Individuals With and Without New-onset Type 1 Diabetes

Context: Alterations in RNA splicing may influence protein isoform diversity that contributes to or reflects the pathophysiology of certain diseases. Whereas specific RNA splicing events in pancreatic islets have been investigated in models of inflammation in vitro, how RNA splicing in the circulation correlates with or is reflective of type 1 diabetes (T1D) disease pathophysiology in humans remains unexplored. Objective: To use machine learning to investigate if alternative RNA splicing events differ between individuals with and without new-onset T1D and to determine if these splicing events provide insight into T1D pathophysiology. Methods: RNA deep sequencing was performed on whole blood samples from 2 independent cohorts: a training cohort consisting of 12 individuals with new-onset T1D and 12 age- and sex-matched nondiabetic controls and a validation cohort of the same size and demographics. Machine learning analysis was used to identify specific isoforms that could distinguish individuals with T1D from controls. Results: Distinct patterns of RNA splicing differentiated participants with T1D from unaffected controls. Notably, certain splicing events, particularly involving retained introns, showed significant association with T1D. Machine learning analysis using these splicing events as features from the training cohort demonstrated high accuracy in distinguishing between T1D subjects and controls in the validation cohort. Gene Ontology pathway enrichment analysis of the retained intron category showed evidence for a systemic viral response in T1D subjects. Conclusion: Alternative RNA splicing events in whole blood are significantly enriched in individuals with new-onset T1D and can effectively distinguish these individuals from unaffected controls. Further, our findings also suggest that RNA splicing profiles offer the potential to provide insights into disease pathogenesis.

60 APPLIED LIFE SCIENCES↗

Genome Analysis of a Verrucomicrobial Endosymbiont With a Tiny Genome Discovered in an Antarctic Lake

Organic Lake in Antarctica is a marine-derived, cold (–13°C), stratified (oxic-anoxic), hypersaline (>200 gl –1 ) system with unusual chemistry (very high levels of dimethylsulfide) that supports the growth of phylogenetically and metabolically diverse microorganisms. Symbionts are not well characterized in Antarctica. However, unicellular eukaryotes are often present in Antarctic lakes and theoretically could harbor endosymbionts. Here, we describe Candidatus Organicella extenuata, a member of the Verrucomicrobia with a highly reduced genome, recovered as a metagenome-assembled genome with genetic code 4 (UGA-to-Trp recoding) from Organic Lake. It is closely related to Candidatus Pinguicocccus supinus (163,218 bp, 205 genes), a newly described cytoplasmic endosymbiont of the freshwater ciliate Euplotes vanleeuwenhoeki (Serra et al., 2020). At 158,228 bp (encoding 194 genes), the genome of Ca. Organicella extenuata is among the smallest known bacterial genomes and similar to the genome of Ca. Pinguicoccus supinus (163,218 bp, 205 genes). Ca. Organicella extenuata retains a capacity for replication, transcription, translation, and protein-folding while lacking any capacity for the biosynthesis of amino acids or vitamins. Notably, the endosymbiont retains a capacity for fatty acid synthesis (type II) and iron–sulfur (Fe-S) cluster assembly. Metagenomic analysis of 150 new metagenomes from Organic Lake and more than 70 other Antarctic aquatic locations revealed a strong correlation in abundance between Ca. Organicella extenuata and a novel ciliate of the genus Euplotes. Like Ca. Pinguicoccus supinus, we infer that Ca. Organicella extenuata is an endosymbiont of Euplotes and hypothesize that both Ca. Organicella extenuata and Ca. Pinguicocccus supinus provide fatty acids and Fe-S clusters to their Euplotes host as the foundation of a mutualistic symbiosis. The discovery of Ca. Organicella extenuata as possessing genetic code 4 illustrates that in addition to identifying endosymbionts by sequencing known symbiotic communities and searching metagenome data using reference endosymbiont genomes, the potential exists to identify novel endosymbionts by searching for unusual coding parameters.

59 BASIC BIOLOGICAL SCIENCES↗

Alteration of gene expression profiles in skeletal muscle of rats exposed to microgravity during a spaceflight

To clarify the mechanism of skeletal muscle wasting during spaceflights, we investigated whether intramuscular gene expression profiles are affected, by using DNA microarray methods. Male rats sent on the 17-day NASA STS-90 Neurolab spaceflight were sacrificed 24 hours after return to earth (MG group). Ground control rats were maintained for 17 days in flight-simulated cages (CS group). Spaceflight induced a 19% and 23% loss of tibialis anterior and gastrocnemius muscle mass, respectively, as compared to ground controls. Muscle RNA was analyzed by the Clontech Atlas DNA expression array in four rats, with two MG/ CS pairs for the tibialis anterior, and one pair for the gastrocnemius. Alterations in gene expression were verified for selected genes by reverse-transcription PCR. In both muscles of MG rats, mRNAs for 12 genes were up-regulated by over 2-fold, and 38 were down-regulated compared to controls. There was inhibition of genes for cell proliferation and growth factor cascades, including cell cycle genes and signal transduction proteins, such as p21 Cip1, retinoblastoma (Rb), cyclins G1/S, -E and -D3, MAP kinase 3, MAD3, and ras related protein RAB2. These data indicate that following exposure to microgravity, there is downregulation of genes involved in regulation of muscle satellite cell replication.

NASA Program Fundamental Space Biology↗

Methylation of p16(INK4a) promoters occurs in vivo in histologically normal human mammary epithelia

Cultures of human mammary epithelial cells (HMECs) contain a subpopulation of variant cells with the capacity to propagate beyond an in vitro proliferation barrier. These variant HMECs, which contain hypermethylated and silenced p16(INK4a) (p16) promoters, eventually accumulate multiple chromosomal changes, many of which are similar to those detected in premalignant and malignant lesions of breast cancer. To determine the origin of these variant HMECs in culture, we used Luria-Delbruck fluctuation analysis and found that variant HMECs exist within the population before the proliferation barrier, thereby raising the possibility that variant HMECs exist in vivo before cultivation. To test this hypothesis, we examined mammary tissue from normal women for evidence of p16 promoter hypermethylation. Here we show that epithelial cells with methylation of p16 promoter sequences occur in focal patches of histologically normal mammary tissue of a substantial fraction of healthy, cancer-free women.

NASA Discipline Radiation Health↗

Ancestral archaea expanded the genetic code with pyrrolysine

The pyrrolysyl-tRNA synthetase (PylRS) facilitates the cotranslational installation of the 22nd amino acid pyrrolysine. Owing to its tolerance for diverse amino acid substrates, and its orthogonality in multiple organisms, PylRS has emerged as a major route to install noncanonical amino acids into proteins in living cells. Recently, a novel class of PylRS enzymes was identified in a subset of methanogenic archaea. Enzymes within this class (ΔPylSn) lack the N-terminal tRNA-binding domain that is widely conserved amongst PylRS enzymes, yet remain active and orthogonal in bacteria and eukaryotes. In this study, we use biochemical and in vivo UAG-readthrough assays to characterize the aminoacylation efficiency and substrate spectrum of a ΔPylSn class PylRS from the archaeon Candidatus Methanomethylophilus alvus. We show that, compared with the full-length enzyme from Methanosarcina mazei, the Ca. M. alvus PylRS displays reduced aminoacylation efficiency but an expanded amino acid substrate spectrum. To gain insight into the evolution of ΔPylSn enzymes, we performed molecular phylogeny using 156 PylRS and 105 pyrrolysine tRNA (tRNA Pyl ) sequences from diverse archaea and bacteria. This analysis suggests that the PylRS•tRNA Pyl pair diverged before the evolution of the three domains of life, placing an early limit on the evolution of the Pyl-decoding trait. Furthermore, our results document the coevolutionary history of PylRS and tRNA Pyl and reveal the emergence of tRNA Pyl sequences with unique A73 and U73 discriminator bases. The orthogonality of these tRNA Pyl species with the more common G73-containing tRNA Pyl will enable future efforts to engineer PylRS systems for further genetic code expansion.

59 BASIC BIOLOGICAL SCIENCES↗

Developmental regulation of the gene for chimeric calcium/calmodulin-dependent protein kinase in anthers

Chimeric Ca(2+)/calmodulin-dependent protein kinase (CCaMK) was cloned from developing anthers of lily (Lilium longiflorum Thumb. cv. Nellie White) and tobacco (Nicotiana tabacum L. cv. Xanthi). Previous biochemical characterization and structure/function studies had revealed that CCaMK has dual modes of regulation by Ca(2+) and Ca(2+)/calmodulin. The unique structural features of CCaMK include a catalytic domain, a calmodulin-binding domain, and a neural visinin-like Ca(2+)-binding domain. The existence of these three features in a single polypeptide distinguishes it from other kinases. Western analysis revealed that CCaMK is expressed in a stage-specific manner in developing anthers. Expression of CCaMK was first detected in pollen mother cells and continued to increase, reaching a peak around the tetrad stage of meiosis. Following microsporogenesis, CCaMK expression rapidly decreased and at later stages of microspore development, no expression was detected. A tobacco genomic clone of CCaMK was isolated and transgenic tobacco plants were produced carrying the CCaMK promoter fused to the beta-glucuronidase reporter gene. Both CCaMK mRNA and protein were detected in the pollen sac and their localizations were restricted to the pollen mother cells and tapetal cells. Consistent results showing a stage-specific expression pattern were obtained by beta-glucuronidase analysis, in-situ hybridization and immunolocalization. The stage- and tissue-specific appearance of CCaMK in anthers suggests that it could play a role in sensing transient changes in free Ca(2+) concentration in target cells, thereby controlling developmental events in the anther.

NASA Discipline Plant Biology↗

Rev1 overexpression accelerates N -methyl- N -nitrosourea (MNU)-induced thymic lymphoma by increasing mutagenesis

Rev1 has two important functions in the translesion synthesis pathway, including dCMP transferase activity, and acts as a scaffolding protein for other polymerases involved in translesion synthesis. However, the role of Rev1 in mutagenesis and tumorigenesis in vivo remains unclear. We previously generated Rev1-overexpressing (Rev1-Tg) mice and reported that they exhibited a significantly increased incidence of intestinal adenoma and thymic lymphoma (TL) after N-methyl-N-nitrosourea (MNU) treatment. In this study, we investigated mutagenesis of MNU-induced TL tumorigenesis in wild-type (WT) and Rev1-Tg mice using diverse approaches, including whole-exome sequencing (WES). In Rev1-Tg TLs, the mutation frequency was higher than that in WT TL in most cases. However, no difference in the number of nonsynonymous mutations in the Catalogue of Somatic Mutations in Cancer (COSMIC) genes was observed, and mutations involved in Notch1 and MAPK signaling were similarly detected in both TLs. Mutational signature analysis of WT and Rev1-Tg TLs revealed cosine similarity with COSMIC mutational SBS5 (aging-related) and SBS11 (alkylation-related). Interestingly, the total number of mutations, but not the genotypes of WT and Rev1-Tg, was positively correlated with the relative contribution of SBS5 in individual TLs, suggesting that genetic instability could be accelerated in Rev1-Tg TLs. Finally, we demonstrated that preleukemic cells could be detected earlier in Rev1-Tg mice than in WT mice, following MNU treatment. In conclusion, Rev1 overexpression accelerates mutagenesis and increases the incidence of MNU-induced TL by shortening the latency period, which may be associated with more frequent DNA damage-induced genetic instability.

60 APPLIED LIFE SCIENCES↗

DeepComplex: A Web Server of Predicting Protein Complex Structures by Deep Learning Inter-chain Contact Prediction and Distance-Based Modelling

Proteins interact to form complexes. Predicting the quaternary structure of protein complexes is useful for protein function analysis, protein engineering, and drug design. However, few user-friendly tools leveraging the latest deep learning technology for inter-chain contact prediction and the distance-based modelling to predict protein quaternary structures are available. To address this gap, we develop DeepComplex, a web server for predicting structures of dimeric protein complexes. It uses deep learning to predict inter-chain contacts in a homodimer or heterodimer. The predicted contacts are then used to construct a quaternary structure of the dimer by the distance-based modelling, which can be interactively viewed and analysed. The web server is freely accessible and requires no registration. It can be easily used by providing a job name and an email address along with the tertiary structure for one chain of a homodimer or two chains of a heterodimer. The output webpage provides the multiple sequence alignment, predicted inter-chain residue-residue contact map, and predicted quaternary structure of the dimer.

59 BASIC BIOLOGICAL SCIENCES↗

Immunological characterization of plant ornithine transcarbamylases

Pea (Pisum sativum L.) ornithine transcarbamylase (OTC) antisera were used to investigate the immunological relatedness of several plant and animal OTC enzymes. The antisera immunoprecipitated OTC activity in all monocot and dicot species tested, and sodium dodecyl sulfate polyacrylamide gel electrophoresis analysis of immunoprecipitated protein revealed monomeric proteins ranging from 35,200 to 36,800 daltons in size. Pea OTC antisera did not recognize mammalian OTC protein. OTC activity and protein levels detected on sodium dodecyl sulfate polyacrylamide gel electrophoresis immunoblots from homogenates of green leaf, etiolated epicotyl and cotyledon, and root tissues of pea were poorly correlated. This might result from differences in amounts of enzymatically active OTC protein in the homogenates. Alternatively, the antisera may fail to recognize different isozyme forms of OTC, which have been reported for some plant species. A putative cytosolic precursor OTC (pOTC) polypeptide exhibiting and Mr = 39,500 to 40,000 daltons was immunoprecipitated from in vitro translation mixtures of total pea leaf poly(A)+ RNA. The size of the pOTC polypeptide, as compared with mature OTC monomer (36,000 daltons), suggests that a 4 kilodalton N-terminal leader sequence, like that responsible for mitochondrial targeting of the mammalian enzyme, may be involved in organellar import of the plant enzyme.

NASA Discipline Plant Biology↗

Over-Expression, Purification and Crystallization of Human Dihydrolipoamide Dehydrogenase

Dehydrolipoamide dehydrogenase (E3; dihydrolipoan-tide:NAD+ oxidoreductase, EC 1.8.1.4) is a common catalytic component found in pyruvate dehydrogenase complex, alpha-ketoglutarate dehydrogenase complex, and branched-chain cc-keto acid dehydrogenase complex. E3 is also a component (referred to as L protein) of the glycine cleavage system in bacterial metabolism (2). Active E3 forms a homodimer with four distinctive subdomain structures (FAD binding, NAD+ binding, central and interface domains) with non-covalently but tightly bound FAD in the holoenzyme. Deduced amino acids from cloned full-length human E3 gene showed a total of 509 amino acids with a leader sequence (N-terminal 35 amino acids) that is excised (mature form) during transportation of expressed E3 into mitochondria membrane. So far, three-dimensional structure of human E3 has not been reported. Our effort to achieve the elucidation of the X-ray crystal structure of human E3 will be presented. Recombinant pPROEX-1 expression vector (from GIBCO BRL Life Technologies) having the human E3 gene without leader sequence was constructed by Polymerase Chain Reaction (PCR) and subsequent ligation, and cloned in E.coli XL1-Blue by transformation. Since pPROEX-1 vector has an internal His-tag (six histidine peptide) located at the upstream region of a multicloning site, one-step affinity purification of E3 using nickelnitriloacetic acid (Ni-NTA) agarose resin, which has a strong affinity to His-tag, was feasible. Also a seven-amino-acid spacer peptide and a recombinant tobacco etch virus protease recognition site (seven amino acids peptide) found between His-tag and first amino acid of expressed E3 facilitated the cleavage of His-tag from E3 after the affinity purification. By IPTG induction, ca. 15 mg of human E3 (mature form) was obtained from 1L LB culture with overnight incubation at 25C. Over 98% of purity of E3 from one-step Ni-NTA agarose affinity purification was confirmed by SDS-PAGE analysis. For crystallization, E3 samples were prepared with and without His-tag. To minimize the aggregation of E3, apo- and holo- forms of E3s were tested, as well as a mutated E3. Dynamic light scattering measurements revealed that the E3 preparations without His-tag and substrate are highly monodispersive with regard to homodimers. Consequent crystallization trials of this E3 preparation led to single crystals of E3 grown by the vapor diffusion method. Crystals were obtained within a few days from solution containing poly (ethylene glycol) monomethyl ether 5000 as a precipitant. Autoindexing and integration of the X-ray diffraction data showed that E3 crystals belong to an orthorhombic system with unit cell parameters a-- 123. 1, b= 165.3 and c=214.3A. Further optimization of protein preparation and crystallization experiments for the structural determination will be discussed.

Hong, Y. S.↗

Standardized Residue Numbering and Secondary Structure Nomenclature in the Class D β-Lactamases

Over 1370 class D β-lactamases are currently known, and they pose a serious threat to the effective treatment of many infectious diseases, particularly in some pathogenic bacteria where evolving carbapenemase activity has been reported. Detailed understanding of their molecular biology, enzymology, and structural biology are critically important, but the lack of a standardized residue numbering scheme and inconsistent secondary structure annotation has made comparative analyses sometimes difficult and cumbersome. Compounding this, in the post-AlphaFold world where we currently find ourselves, an extraordinary wealth of detailed structural information on these enzymes is literally at our fingertips; therefore it is vitally important that a standard numbering system is in place to facilitate the accurate and straightforward analysis of their structures. In conclusion, here we present a residue numbering and secondary structure scheme for the class D enzymes based on the sequence and structure of OXA-48 and apply it to test targets to demonstrate the ease with which it can be used.

59 BASIC BIOLOGICAL SCIENCES↗

Transcriptomic Analysis of the CAM Species Kalanchoë fedtschenkoi Under Low- and High-Temperature Regimes

Temperature stress is one of the major limiting environmental factors that negatively impact global crop yields. Kalanchoë fedtschenkoi is an obligate crassulacean acid metabolism (CAM) plant species, exhibiting much higher water-use efficiency and tolerance to drought and heat stresses than C 3 or C 4 plant species. Previous studies on gene expression responses to low- or high-temperature stress have been focused on C 3 and C 4 plants. There is a lack of information about the regulation of gene expression by low and high temperatures in CAM plants. To address this knowledge gap, we performed transcriptome sequencing (RNA-Seq) of leaf and root tissues of K. fedtschenkoi under cold (8 °C), normal (25 °C), and heat (37 °C) conditions at dawn (i.e., 2 h before the light period) and dusk (i.e., 2 h before the dark period). Our analysis revealed differentially expressed genes (DEGs) under cold or heat treatment in comparison to normal conditions in leaf or root tissue at each of the two time points. In particular, DEGs exhibiting either the same or opposite direction of expression change (either up-regulated or down-regulated) under cold and heat treatments were identified. In addition, we analyzed gene co-expression modules regulated by cold or heat treatment, and we performed in-depth analyses of expression regulation by temperature stresses for selected gene categories, including CAM-related genes, genes encoding heat shock factors and heat shock proteins, circadian rhythm genes, and stomatal movement genes. Our study highlights both the common and distinct molecular strategies employed by CAM and C 3 /C 4 plants in adapting to extreme temperatures, providing new insights into the molecular mechanisms underlying temperature stress responses in CAM species.

59 BASIC BIOLOGICAL SCIENCES↗

Kingdom-Wide Analysis of Fungal Protein-Coding and tRNA Genes Reveals Conserved Patterns of Adaptive Evolution

Protein-coding genes evolved codon usage bias due to the combined but uneven effects of adaptive and nonadaptive influences. Studies in model fungi agree on codon usage bias as an adaptation for fine-tuning gene expression levels; however, such knowledge is lacking for most other fungi. Our comparative genomics analysis of over 450 species supports codon usage and transfer RNAs (tRNAs) as coadapted for translation speed and this is most likely a realization of convergent evolution. Rather than drift, phylogenetic reconstruction inferred adaptive radiation as the best explanation for the variation of interspecific codon usage bias. Although the phylogenetic signals for individual codon and tRNAs frequencies are lower than expected by genetic drift, we found remarkable conservation of highly expressed genes being codon optimized for translation by the most abundant tRNAs, especially by inosine-modified tRNAs. As an application, we present a sequence-to-expression neural network that uses codons to reliably predict highly expressed transcripts. The kingdom Fungi, with over a million species, includes many key players in various ecosystems and good targets for biotechnology. Collectively, our results have implications for better understanding the evolutionary success of fungi, as well as informing the biosynthetic manipulation of fungal genes.

59 BASIC BIOLOGICAL SCIENCES↗