The N-terminal intrinsically disordered region of Ncb5or docks with the cytochrome b5 core to form a helical motif that is of ancient origin.

NADH cytochrome b5 oxidoreductase (Ncb5or) is a cytosolic ferric reductase implicated in diabetes and neurological conditions. Ncb5or comprises cytochrome b5 (b5 ) and cytochrome b5 reductase (b5 R) domains separated by a CHORD-Sgt1 (CS) linker domain. Ncb5or redox activity depends on proper inter-domain interactions to mediate electron transfer from NADH or NADPH via FAD to heme. While full-length human Ncb5or has proven resistant to crystallization, we have succeeded in obtaining high-resolution atomic structures of the b5 domain and a construct containing the CS and b5 R domains (CS/b5 R). Ncb5or also contains an N-terminal intrinsically disordered region of 50 residues that has no homologs in other protein families in animals but features a distinctive, conserved L34 MDWIRL40 motif also present in reduced lateral root formation (RLF) protein in rice and increased recombination center 21 in baker's yeast, all attaching to a b5 domain. After unsuccessful attempts at crystallizing a human Ncb5or construct comprising the N-terminal region naturally fused to the b5 domain, we were able to obtain a high-resolution atomic structure of a recombinant rice RLF construct corresponding to residues 25-129 of human Ncb5or (52% sequence identity; 74% similarity). The structure reveals Trp120 (corresponding to invariant Trp37 in Ncb5or) to be part of an 11-residue α-helix (S116 QMDWLKLTRT126 ) packing against two of the four helices in the b5 domain that surround heme (α2 and α5). The Trp120 side chain forms a network of interactions with the side chains of four highly conserved residues corresponding to Tyr85 and Tyr88 (α2), Cys124 (α5), and Leu47 in Ncb5or. Circular dichroism measurements of human Ncb5or fragments further support a key role of Trp37 in nucleating the formation of the N-terminal helix, whose location in the N/b5 module suggests a role in regulating the function of this multi-domain redox enzyme. This study revealed for the first time an ancient origin of a helical motif in the N/b5 module as reflected by its existence in a class of cytochrome b5 proteins from three kingdoms among eukaryotes.

Role of global public sector research in discovering new drugs and vaccines.

Analysis of international public-sector contributions to Food and Drug Administration (FDA)-approved drugs and vaccines allows for a more thorough examination of the global biomedical innovation ecosystem by institution of origin. Using new and existing methods, we have identified 364 FDA-approved drugs and vaccines approved from 1973 to 2016 discovered in whole or in part by Public Sector Research Institutions (PSRIs) worldwide. We identified product-specific intellectual property contributions to FDA-approved small molecule and biologic drugs and vaccines from the FDA Orange Book, our peer network, published studies, and three new sources: reports of medical product manufacturers' payments to physicians and teaching hospitals under The Sunshine Act of 2010, a paper by Kneller and 64 royalty monetization transactions by academic institutions and/or their faculty that one of us (AS) maintains. We include a total of 293 drugs discovered either wholly by a US PSRI or jointly by a U.S. and a non-U.S. PSRI. 119 FDA-approved drugs and vaccines were discovered by PSRIs outside the U.S. Of these, 71 were solely discovered outside the US, while 48 also involved intellectual property contributions by US PSRIs. In the context of the global public sector landscape, the US dominates drug discovery, accounting for two-thirds of these drugs and many of the important, innovative vaccines introduced over the past 30 years. Contributions by Canada, UK, Germany, Belgium, Japan, and others each amount to 5.4% or less of the total. Supplementary Information: The online version contains supplementary material available at 10.1007/s10961-023-10007-z.

A Review on Electric and Fuel Cell Vehicle Anatomy, Technology Evolution and Policy Drivers towards EVs and FCEVs Market Propagation.

The transportation sector is the largest consumer of fossil fuels; making it a major producer of greenhouse gases. Due to declining fossil fuel reserves and increasingly stringent vehicle emission regulations globally, it is essential to shift to alternative energy sources. Economic and eco-friendly fuel-efficient hybrid, electric, and fuel cell vehicles are regarded as one of the best alternative solutions to cope with the government policies and to reduce the rise in global temperature caused by the automotive sector. Technological advancements in fuel cells, batteries, and chargers have further supported the development of electric vehicles. The major challenges of range and charging time in electric vehicles can be countered by range extension technology and developing all-electric hybrid vehicles. In this review, a comprehensive study of different type of vehicles and their architectures are presented. Insights on energy storage devices and converters of electric vehicles currently in use were also provided. Furthermore, various fuel cell advancements and the technical challenges faced during the commercialization of fuel cell vehicles were highlighted.

Fast-changing life-styles and ecotoxicity of hair dyes drive the emergence of hidden toxicants threatening environmental sustainability in Asia.

The practice of hair dyeing is a rapidly expanding industry on a global scale; however, it has become a major concern for Asian countries because they have been undergoing rapid transformations of their environment and lifestyles. While the socio-economic benefits and impacts of this globalization trend are widely understood, the environmental effects are largely unknown. In particular, commonly available oxidative dyes potentially pose specific environmental risks due to their use of a toxic aromatic amine p-Phenylenediamine (PPD). In investigating the environmental impacts of PPD chemicals, we first provide context to the study by setting out the socio-psychological drivers to industrial expansion in Asian countries along with an overview of research into its effects, to show that its environmental impacts are under-researched. We then investigate the environmental toxicity of PPD by focusing on the role of microbes in metabolizing waste products. Results show that Acinetobacter baumannii EB1 isolated from dye effluent prevents autoxidation of PPD under oxygen-enriched (shaking) or oxygen-deficient (static) conditions representing different environmental settings. Microbes transformed PPD into more toxic metabolites, which then significantly reduced plant growth, thereby having a direct bearing on ecosystem services. Based on the findings, we argue that stricter regulatory controls on hair dye wastewater are necessary, particularly in newly industrialising Asian countries where the expansion of commercial practice is most prevalent.

19th International Meeting on Frankia and Actinorhizal Plants.

It has been 40 years since the first meeting dedicated to Frankia and actinorhizal plants, which was held at Petersham, Massachusetts (reported in Torrey and Tjepkema, 1979). Since then biennial meetings have been organised and held in different venues around the globe (Table 1). The most recent meeting, the "19th International Meeting on Frankia and Actinorhizal Plants", organised in Hammamet, Tunisia from 17th to 19th of March, 2018, gathered scientists from Algeria, Argentina, Belgium, China, Egypt, France, India, Portugal, Senegal, Sweden, UK, USA and Tunisia. The event was a stimulating opportunity for active researchers to share many advances since the previous meeting held in Montpellier, France (Franche et al. 2016) and to discuss new perspectives in this research field.

Clinical and financial implications of emergency department visits for synthetic marijuana.

BACKGROUND: Many users believe that synthetic cannabinoids offer a safe and legal means of getting high. However, spikes in emergency department visits have been associated with use of synthetic cannabinoids. The purpose of the current study was to document emergency department visits from three large hospitals in one metropolitan area over a two month period. METHOD: This was a retrospective chart review examining 218 patients presenting to three inner city emergency departments between March and April 2014. Data collected included demographic information, information regarding ED diagnosis and treatment, signs and symptoms, ancillary testing, ED disposition, and cost of the medical treatment. RESULTS: The majority of patients (75.7%) were discharged after ED workup, but 12.4% were admitted for medical treatment and 11.5% were admitted for psychiatric treatment. Ten patients (4.6%) were admitted to the ICU. Symptoms experienced most frequently include: hypertension, tachycardia, agitation, drowsiness, nausea, and confusion. Cluster analysis revealed four symptom clusters of individuals presenting after using synthetic cannabinoids: 1) confusion, hostility, agitation, 2) nausea, vomiting, abdominal pain, 3) drowsiness, and 4) the absence of these symptoms. CONCLUSION: This study has three important findings. First, significant ED resources are being used to treat individuals presenting due to effects of synthetic cannabis. Second, synthetic cannabis is not a benign substance. Third, while the hostile and agitated user is generally presented in the media, this study finds significant heterogeneity in presentation. Further research is needed to fully understand the implications of synthetic cannabinoid use.

Automated Estimation of Aquifer Parameters from Arbitrary-Rate Pumping Tests in Python and MATLAB.

Inspired by the analysis by Mishra et al. (2012) of variable pumping rate tests using piecewise-linear reconstructions of the pumping history, this article contains a derivation of the convolutional form of pumping tests in which the pumping history may take any possible form. The solution is very similar to the classical Theis (1935) equation but uses the Green's function for a pumped aquifer given by taking the time derivative of the well function W ( u ( t ) ) . This eliminates one integration inside another and renders the convolution including the pumping history about as computationally demanding as calculating the well function alone, so that the convolution can be completed using handy mathematical software. It also allows nonlinear well losses, and because an easily-computed deterministic model exists for all data points and pumping history, an objective function may include all data, so that errors are reduced in calculating any nonlinear-well losses. In addition, data from multiple observation wells may be used simultaneously in the inversion. We provide codes in MATLAB and Python to solve for drawdown resulting from an arbitrary pumping history and compute the optimal aquifer parameters to fit the data. We find that the subtleties in parameter dependencies and constructing an appropriate objective function have a substantial effect on the interpreted parameters. Furthermore, the optimization from step-drawdown tests is typically nonunique and strongly suggests that a Bayesian inversion should be used to fully estimate the joint probability density of the parameter vector.

Simulating Metal-Imidazole Complexes.

One commonly observed binding motif in metalloproteins involves the interaction between a metal ion and histidine's imidazole side chains. Although previous imidazole-M(II) parameters established the flexibility and reliability of the 12-6-4 Lennard-Jones (LJ)-type nonbonded model by simply tuning the ligating atom's polarizability, they have not been applied to multiple-imidazole complexes. To fill this gap, we systematically simulate multiple-imidazole complexes (ranging from one to six) for five metal ions (Co(II), Cu(II), Mn(II), Ni(II), and Zn(II)) which commonly appear in metalloproteins. Using extensive (40 ns per PMF window) sampling to assemble free energy association profiles (using OPC water and standard HID imidazole charge models from AMBER) and comparing the equilibrium distances to DFT calculations, a new set of parameters was developed to focus on energetic and geometric features of multiple-imidazole complexes. The obtained free energy profiles agree with the experimental binding free energy and DFT calculated distances. To validate our model, we show that we can close the thermodynamic cycle for metal-imidazole complexes with up to six imidazole molecules in the first solvation shell. Given the success in closing the thermodynamic cycles, we then used the same extended sampling method for six other metal ions (Ag(I), Ca(II), Cd(II), Cu(I), Fe(II), and Mg(II)) to obtain new parameters. Since these new parameters can reproduce the one-imidazole geometry and energy accurately, we hypothesize that they will reasonably predict the binding free energy of higher-level coordination numbers. Hence, we did not extend the analysis of these ions up to six imidazole complexes. Overall, the results shed light on metal-protein interactions by emphasizing the importance of ligand-ligand interaction and metal-π-stacking within metalloproteins.

Communication: A full solution of the annihilation reaction A + B → [empty-set] based on time-subordination.

The connection between the governing equations of chemical reaction and the underlying stochastic processes of particle collision and transformation have been developed previously along two end-member conditions: perfectly mixed and maximally diffusion-limited. The complete governing equation recognizes that in the perfectly mixed case, the particle (i.e., molecular or macro-particle) number state evolution is markovian, but that spatial self-organization of reactants decreases the probability of reactant pairs finding themselves co-located. This decreased probability manifests itself as a subordination of the clock time: as reactant concentrations become spatially variable (unmixed), the time required for reactants to find each other increases and the random operational time that particles spend in the active reaction process is less than the clock time. For example, in the system A + B → [empty-set], a simple approximate calculation for the return time of a brownian motion to a moving boundary allows a calculation of the operational time density, and the total solution is a subordination integral of the perfectly-mixed solution with a modified inverse gaussian subordinator. The system transitions from the well-mixed solution to the asymptotic diffusion-limited solution that decays as t(-d∕4) in d-dimensions.

Hydraulic Conductivity Fields: Gaussian or Not?

Hydraulic conductivity (K) fields are used to parameterize groundwater flow and transport models. Numerical simulations require a detailed representation of the K field, synthesized to interpolate between available data. Several recent studies introduced high resolution K data (HRK) at the Macro Dispersion Experiment (MADE) site, and used ground-penetrating radar (GPR) to delineate the main structural features of the aquifer. This paper describes a statistical analysis of these data, and the implications for K field modeling in alluvial aquifers. Two striking observations have emerged from this analysis. The first is that a simple fractional difference filter can have a profound effect on data histograms, organizing non-Gaussian ln K data into a coherent distribution. The second is that using GPR facies allows us to reproduce the significantly non-Gaussian shape seen in real HRK data profiles, using a simulated Gaussian ln K field in each facies. This illuminates a current controversy in the literature, between those who favor Gaussian ln K models, and those who observe non-Gaussian ln K fields. Both camps are correct, but at different scales.

Fractional calculus in hydrologic modeling: A numerical perspective.

Fractional derivatives can be viewed either as handy extensions of classical calculus or, more deeply, as mathematical operators defined by natural phenomena. This follows the view that the diffusion equation is defined as the governing equation of a Brownian motion. In this paper, we emphasize that fractional derivatives come from the governing equations of stable Lévy motion, and that fractional integration is the corresponding inverse operator. Fractional integration, and its multi-dimensional extensions derived in this way, are intimately tied to fractional Brownian (and Lévy) motions and noises. By following these general principles, we discuss the Eulerian and Lagrangian numerical solutions to fractional partial differential equations, and Eulerian methods for stochastic integrals. These numerical approximations illuminate the essential nature of the fractional calculus.

Mixing-Driven Equilibrium Reactions in Multidimensional Fractional Advection Dispersion Systems.

We study instantaneous, mixing-driven, bimolecular equilibrium reactions in a system where transport is governed by a multidimensional space fractional dispersion equation. The superdiffusive, nonlocal nature of the system causes the location and magnitude of reactions that take place to change significantly from a classical Fickian diffusion model. In particular, regions where reaction rates would be zero for the Fickian case become regions where the maximum reaction rate occurs when anomalous dispersion operates. We also study a global metric of mixing in the system, the scalar dissipation rate and compute its asymptotic scaling rates analytically. The scalar dissipation rate scales asymptotically as t-(d+α)/α , where d is the number of spatial dimensions and α is the fractional derivative exponent.

Growth and development of Frankia spp. strain CcI3 at the single-hypha level in liquid culture.

Filamentous actinobacteria from the genus Frankia grow by hyphal tip extension and branching. The growth kinetics and branching pattern of Frankia are not well studied, especially at the early stages of mycelial development. Here, we compare the growth of Frankia sp. strain CcI3 in liquid cultures with and without proteose peptone #3 (PP3) using time-lapse photomicrography and image analysis. Individual hyphae showed a pseudolinear increase in length at early stages of development, whereas at the mycelial level, the aggregate length of hyphae described an exponential rate before slowing. Growth based on optical density or microscopic observations was similar in medium with or without PP3. However, PP3 altered the pattern of mycelial development by increasing branching. Distances between the hyphal apex and first branches were on average shorter in PP3-containing media. The final interbranch distances were also shorter in PP3 medium indicating that hyphae tended to branch earlier and more often when supplemented with PP3 to give a more compact mycelium. Vesicle development in nitrogen-fixing cultures limited cell expansion as a result of vesicles truncating growth on new branches. The results provide some explanation for the growth kinetics of Frankia and some indication of how growth rates may be improved.

A microarray biosensor for multiplexed detection of microbes using grating-coupled surface plasmon resonance imaging.

Grating-coupled surface plasmon resonance imaging (GCSPRI) utilizes an optical diffraction grating embossed on a gold-coated sensor chip to couple collimated incident light into surface plasmons. The angle at which this coupling occurs is sensitive to the capture of analyte at the chip surface. This approach permits the use of disposable biosensor chips that can be mass-produced at low cost and spotted in microarray format to greatly increase multiplexing capabilities. The current GCSPRI instrument has the capacity to simultaneously measure binding at over 1000 unique, discrete regions of interest (ROIs) by utilizing a compact microarray of antibodies or other specific capture molecules immobilized on the sensor chip. In this report, we describe the use of GCSPRI to directly detect multiple analytes over a large dynamic range, including soluble protein toxins, bacterial cells, and viruses, in near real-time. GCSPRI was used to detect a variety of agents that would be useful for diagnostic and environmental sensing purposes, including macromolecular antigens, a nontoxic form of Pseudomonas aeruginosa exotoxin A (ntPE), Bacillus globigii, Mycoplasma hyopneumoniae, Listeria monocytogenes, Escherichia coli, and M13 bacteriophage. These studies indicate that GCSPRI can be used to simultaneously assess the presence of toxins and pathogens, as well as quantify specific antibodies to environmental agents, in a rapid, label-free, and highly multiplexed assay requiring nanoliter amounts of capture reagents.

Effect of actinorhizal root exudates on the proteomes of Frankia soli NRRL B-16219, a strain colonizing the root tissues of its actinorhizal host via intercellular pathway.

Frankia and actinorhizal plants exchange signals in the rhizosphere leading to specific mutual recognition of partners and nitrogen-fixing nodule organogenesis. Frankia soli strain NRRL B-16219, from the Elaeagnus specificity group, colonizes the root tissues of its actinorhizal host through direct intercellular penetration of root epidermis cells and cortex. Here, we studied the early proteogenomic response of strain NRRL B-16219 to treatment with root exudates from compatible Elaeagnus angustifolia, and incompatible Ceanothus thyrsiflorus and Coriaria myrtifolia, host plants grown in nitrogen depleted hydroponic medium. Next-generation proteomics was used to identify the main Frankia proteins differentially expressed in response to the root exudates. No products of the nod genes present in B-16219 were detected. Proteins specifically upregulated in presence of E. angustifolia root exudates include those connected to nitrogen fixation and assimilation (glutamate synthetase, hydrogenase and squalene synthesis), respiration (oxidative phosphorylation and citric acid cycle pathways), oxidative stress (catalase, superoxide dismutase, and peroxidase), proteolysis (proteasome, protease, and peptidase) and plant cell wall degrading proteins involved in the depolymerization of celluloses (endoglucanase, glycosyltransferase, beta-mannanases, glycoside hydrolase and glycosyl hydrolase). Proteomic data obtained in this study will help link signaling molecules/factors to their biosynthetic pathways once those factors have been fully characterized.

Genome sequence of "Candidatus Frankia datiscae" Dg1, the uncultured microsymbiont from nitrogen-fixing root nodules of the dicot Datisca glomerata.

Members of the noncultured clade of Frankia enter into root nodule symbioses with actinorhizal species from the orders Cucurbitales and Rosales. We report the genome sequence of a member of this clade originally from Pakistan but obtained from root nodules of the American plant Datisca glomerata without isolation in culture.

Accommodating a nonconservative internal mutation by water-mediated hydrogen bonding between ß-sheet strands: a comparison of human and rat type B (mitochondrial) cytochrome b5.

Mammalian type B (mitochondrial) b(5) cytochromes exhibit greater amino acid sequence diversity than their type A (microsomal) counterparts, as exemplified by the type B proteins from human (hCYB5B) and rat (rCYB5B). The comparison of X-ray crystal structures of hCYB5B and rCYB5B reported herein reveals a striking difference in packing involving the five-strand ß-sheet, which can be attributed to fully buried residue 21 in strand ß4. The greater bulk of Leu21 in hCYB5B in comparison to that of Thr21 in rCYB5B results in a substantial displacement of the first two residues in ß5, and consequent loss of two of the three hydrogen bonds between ß5 and ß4. Hydrogen bonding between the residues is instead mediated by two well-ordered, fully buried water molecules. In a 10 ns molecular dynamics simulation, one of the buried water molecules in the hCYB5B structure exchanged readily with solvent via intermediates having three water molecules sandwiched between ß4 and ß5. When the buried water molecules were removed prior to a second 10 ns simulation, ß4 and ß5 formed persistent hydrogen bonds identical to those in rCYB5B, but the Leu21 side chain was forced to adopt a rarely observed conformation. Despite the apparently greater ease of access of water to the interior of hCYB5B than of rCYB5B suggested by these observations, the two proteins exhibit virtually identical stability, dynamic, and redox properties. The results provide new insight into the factors stabilizing the cytochrome b(5) fold.

Study of the individual cytochrome b5 and cytochrome b5 reductase domains of Ncb5or reveals a unique heme pocket and a possible role of the CS domain.

NADH cytochrome b(5) oxidoreductase (Ncb5or) is found in animals and contains three domains similar to cytochrome b(5) (b(5)), CHORD-SGT1 (CS), and cytochrome b(5) reductase (b(5)R). Ncb5or has an important function, as suggested by the diabetes and lipoatrophy phenotypes in Ncb5or null mice. To elucidate the structural and functional properties of human Ncb5or, we generated its individual b(5) and b(5)R domains (Ncb5or-b(5) and Ncb5or-b(5)R, respectively) and compared them with human microsomal b(5) (Cyb5A) and b(5)R (Cyb5R3). A 1.25 Å x-ray crystal structure of Ncb5or-b(5) reveals nearly orthogonal planes of the imidazolyl rings of heme-ligating residues His(89) and His(112), consistent with a highly anisotropic low spin EPR spectrum. Ncb5or is the first member of the cytochrome b(5) family shown to have such a heme environment. Like other b(5) family members, Ncb5or-b(5) has two helix-loop-helix motifs surrounding heme. However, Ncb5or-b(5) differs from Cyb5A with respect to location of the second heme ligand (His(112)) and of polypeptide conformation in its vicinity. Electron transfer from Ncb5or-b(5)R to Ncb5or-b(5) is much less efficient than from Cyb5R3 to Cyb5A, possibly as a consequence of weaker electrostatic interactions. The CS linkage probably obviates the need for strong interactions between b(5) and b(5)R domains in Ncb5or. Studies with a construct combining the Ncb5or CS and b(5)R domains suggest that the CS domain facilitates docking of the b(5) and b(5)R domains. Trp(114) is an invariant surface residue in all known Ncb5or orthologs but appears not to contribute to electron transfer from the b(5)R domain to the b(5) domain.

The biology of Frankia sp. strains in the post-genome era.

Progress in understanding symbiotic determinants involved in the N(2)-fixing actinorhizal plant symbioses has been slow but steady. Problems persist with studying the bacterial contributions to the symbiosis using traditional microbiological techniques. However, recent years have seen the emergence of several genomes from Frankia sp. strains and the development of techniques for manipulating plant gene expression. Approaches to understanding the bacterial side of the symbiosis have employed a range of techniques that reveal the proteomes and transcriptomes from both cultured and symbiotic frankiae. The picture beginning to emerge provides some perspective on the heterogeneity of frankial populations in both conditions. In general, frankial populations in root nodules seem to maintain a rather robust metabolism that includes nitrogen fixation and substantial biosynthesis and energy-generating pathways, along with a modified ammonium assimilation program. To date, particular bacterial genes have not been implicated in root nodule formation but some hypotheses are emerging with regard to how the plant and microorganism manage to coexist. In particular, frankiae seem to present a nonpathogenic presence to the plant that may have the effect of minimizing some plant defense responses. Future studies using high-throughput approaches will likely clarify the range of bacterial responses to symbiosis that will need to be understood in light of the more rapidly advancing work on the plant host.

Transcriptomes of Frankia sp. strain CcI3 in growth transitions.

BACKGROUND: Frankia sp. strains are actinobacteria that form N2-fixing root nodules on angiosperms. Several reference genome sequences are available enabling transcriptome studies in Frankia sp. Genomes from Frankia sp. strains differ markedly in size, a consequence proposed to be associated with a high number of indigenous transposases, more than 200 of which are found in Frankia sp. strain CcI3 used in this study. Because Frankia exhibits a high degree of cell heterogeneity as a consequence of its mycelial growth pattern, its transcriptome is likely to be quite sensitive to culture age. This study focuses on the behavior of the Frankia sp. strain CcI3 transcriptome as a function of nitrogen source and culture age. RESULTS: To study global transcription in Frankia sp. CcI3 grown under different conditions, complete transcriptomes were determined using high throughput RNA deep sequencing. Samples varied by time (five days vs. three days) and by culture conditions (NH4+ added vs. N2 fixing). Assembly of millions of reads revealed more diversity of gene expression between five-day and three-day old cultures than between three day old cultures differing in nitrogen sources. Heat map analysis organized genes into groups that were expressed or repressed under the various conditions compared to median expression values. Twenty-one SNPs common to all three transcriptome samples were detected indicating culture heterogeneity in this slow-growing organism. Significantly higher expression of transposase ORFs was found in the five-day and N2-fixing cultures, suggesting that N starvation and culture aging provide conditions for on-going genome modification. Transposases have previously been proposed to participate in the creating the large number of gene duplication or deletion in host strains. Subsequent RT-qPCR experiments confirmed predicted elevated transposase expression levels indicated by the mRNA-seq data. CONCLUSIONS: The overall pattern of gene expression in aging cultures of CcI3 suggests significant cell heterogeneity even during normal growth on ammonia. The detection of abundant transcription of nif (nitrogen fixation) genes likely reflects the presence of anaerobic, N-depleted microsites in the growing mycelium of the culture, and the presence of significantly elevated transposase transcription during starvation indicates the continuing evolution of the Frankia sp. strain CcI3 genome, even in culture, especially under stressed conditions. These studies also sound a cautionary note when comparing the transcriptomes of Frankia grown in root nodules, where cell heterogeneity would be expected to be quite high.