Conservation of lipid functions in cytochrome bc complexes

Lipid binding sites and properties are compared in two sub-families of hetero-oligomeric membrane protein complexes known to have similar functions in order to gain further understanding of the role of lipid in the function, dynamics, and assembly of these complexes. Using the crystal structure information for both complexes, we compared the lipid binding properties of the cytochrome b₆f and bc₁ complexes that function in photosynthetic and respiratory membrane energy transduction. Comparison of lipid and detergent binding sites in the b₆f complex with those in bc₁ shows significant conservation of lipid positions. Seven lipid binding sites in the cyanobacterial b₆f complex overlap three natural sites in the Chlamydomonas reinhardtii algal complex and four sites in the yeast mitochondrial bc₁ complex. The specific identity of lipids is different in b₆f and bc₁ complexes: b₆f contains sulfoquinovosyldiacylglycerol, phosphatidylglycerol, phosphatidylcholine, monogalactosyldiacylglycerol, and digalactosyldiacylglycerol, whereas cardiolipin, phosphatidylethanolamine, and phosphatidic acid are present in the yeast bc₁ complex. The lipidic chlorophyll a and β-carotene (β-car) in cyanobacterial b₆f, as well as eicosane in C. reinhardtii, are unique to the b₆f complex. Inferences of lipid binding sites and functions were supported by sequence, interatomic distance, and B-factor information on interacting lipid groups and coordinating amino acid residues. The lipid functions inferred in the b₆f complex are as follows: (i) substitution of a transmembrane helix by a lipid and chlorin ring, (ii) lipid and β-car connection of peripheral and core domains, (iii) stabilization of the iron-sulfur protein transmembrane helix, (iv) n-side charge and polarity compensation, and (v) β-car-mediated super-complex with the photosystem I complex.     

Association of P2Y2 receptor SNPs with bone mineral density and osteoporosis risk in a cohort of Dutch fracture patients

The P2Y₂ receptor is a G-protein-coupled receptor with adenosine 5′-triphosphate (and UTP) as natural ligands. It is thought to be involved in bone physiology in an anti-osteogenic manner. As several non-synonymous single nucleotide polymorphisms (SNPs) have been identified within the P2Y₂ receptor gene in humans, we examined associations between genetic variations in the P2Y₂ receptor gene and bone mineral density (BMD) (i.e., osteoporosis risk), in a cohort of fracture patients. Six hundred and ninety women and 231 men aged ≥50 years, visiting an osteoporosis outpatient clinic at Maastricht University Medical Centre for standard medical follow-up after a recent fracture, were genotyped for three non-synonymous P2Y₂ receptor gene SNPs. BMD was measured at three locations (total hip, lumbar spine, and femoral neck) using dual-energy X-ray absorptiometry. Differences in BMD between different genotypes were tested using analysis of covariance. In women, BMD values at all sites were significantly different between the genotypes for the Leu46Pro polymorphism, with women homozygous for the variant allele showing the highest BMD values (0.05 > p > 0.01). The Arg312Ser and Arg334Cys polymorphisms showed no differences in BMD values between the different genotypes. This is the first report that describes the association between the Leu46Pro polymorphism of the human P2Y₂ receptor and the risk of osteoporosis.     

Ascochlorin is a novel, specific inhibitor of the mitochondrial cytochrome bc1 complex

Ascochlorin is an isoprenoid antibiotic that is produced by the phytopathogenic fungus Ascochyta viciae. Similar to ascofuranone, which specifically inhibits trypanosome alternative oxidase by acting at the ubiquinol binding domain, ascochlorin is also structurally related to ubiquinol. When added to the mitochondrial preparations isolated from rat liver, or the yeast Pichia (Hansenula) anomala, ascochlorin inhibited the electron transport via CoQ in a fashion comparable to antimycin A and stigmatellin, indicating that this antibiotic acted on the cytochrome bc₁ complex. In contrast to ascochlorin, ascofuranone had much less inhibition on the same activities. On the one hand, like the Q_i site inhibitors antimycin A and funiculosin, ascochlorin induced in H. anomala the expression of nuclear-encoded alternative oxidase gene much more strongly than the Q_o site inhibitors tested. On the other hand, it suppressed the reduction of cytochrome b and the generation of superoxide anion in the presence of antimycin A₃ in a fashion similar to the Q_o site inhibitor myxothiazol. These results suggested that ascochlorin might act at both the Q_i and the Q_o sites of the fungal cytochrome bc₁ complex. Indeed, the altered electron paramagnetic resonance (EPR) lineshape of the Rieske iron-sulfur protein, and the light-induced, time-resolved cytochrome b and c reduction kinetics of Rhodobacter capsulatus cytochrome bc₁ complex in the presence of ascochlorin demonstrated that this inhibitor can bind to both the Q_o and Q_i sites of the bacterial enzyme. Additional experiments using purified bovine cytochrome bc₁ complex showed that ascochlorin inhibits reduction of cytochrome b by ubiquinone through both Q_i and Q_o sites. Moreover, crystal structure of chicken cytochrome bc₁ complex treated with excess ascochlorin revealed clear electron densities that could be attributed to ascochlorin bound at both the Q_i and Q_o sites. Overall findings clearly show that ascochlorin is an unusual cytochrome bc₁ inhibitor that acts at both of the active sites of this enzyme.     

Evolution of I-SceI homing endonucleases with increased DNA recognition site specificity

Elucidating how homing endonucleases undergo changes in recognition site specificity will facilitate efforts to engineer proteins for gene therapy applications. I-SceI is a monomeric homing endonuclease that recognizes and cleaves within an 18-bp target. It tolerates limited degeneracy in its target sequence, including substitution of a C:G₊₄ base pair for the wild-type A:T₊₄ base pair. Libraries encoding randomized amino acids at I-SceI residue positions that contact or are proximal to A:T₊₄ were used in conjunction with a bacterial one-hybrid system to select I-SceI derivatives that bind to recognition sites containing either the A:T₊₄ or the C:G₊₄ base pairs. As expected, isolates encoding wild-type residues at the randomized positions were selected using either target sequence. All I-SceI proteins isolated using the C:G₊₄ recognition site included small side-chain substitutions at G100 and either contained (K86R/G100T, K86R/G100S and K86R/G100C) or lacked (G100A, G100T) a K86R substitution. Interestingly, the binding affinities of the selected variants for the wild-type A:T₊₄ target are 4- to 11-fold lower than that of wild-type I-SceI, whereas those for the C:G₊₄ target are similar. The increased specificity of the mutant proteins is also evident in binding experiments in vivo. These differences in binding affinities account for the observed ∼36-fold difference in target preference between the K86R/G100T and wild-type proteins in DNA cleavage assays. An X-ray crystal structure of the K86R/G100T mutant protein bound to a DNA duplex containing the C:G₊₄ substitution suggests how sequence specificity of a homing enzyme can increase. This biochemical and structural analysis defines one pathway by which site specificity is augmented for a homing endonuclease.     

Structure-Function Relationships in Fungal Large-Subunit Catalases

Neurospora crassa has two large-subunit catalases, CAT-1 and CAT-3. CAT-1 is associated with non-growing cells and accumulates particularly in asexual spores; CAT-3 is associated with growing cells and is induced under different stress conditions. It is our interest to elucidate the structure-function relationships in large-subunit catalases. Here we have determined the CAT-3 crystal structure and compared it with the previously determined CAT-1 structure. Similar to CAT-1, CAT-3 hydrogen peroxide (H₂O₂) saturation kinetics exhibited two components, consistent with the existence of two active sites: one saturated in the millimolar range and the other in the molar range. In the CAT-1 structure, we found three interesting features related to its unusual kinetics: (a) a constriction in the channel that conveys H₂O₂ to the active site; (b) a covalent bond between the tyrosine, which forms the fifth coordination bound to the iron of the heme, and a vicinal cysteine; (c) oxidation of the pyrrole ring III to form a cis-hydroxyl group in C5 and a cis-γ-spirolactone in C6. The site of heme oxidation marks the starts of the central channel that communicates to the central cavity and the shortest way products can exit the active site. CAT-3 has a similar constriction in its major channel, which could function as a gating system regulated by the H₂O₂ concentration before the gate. CAT-3 functional tyrosine is not covalently bonded, but has instead the electron relay mechanism described for the human catalase to divert electrons from it. Pyrrole ring III in CAT-3 is not oxidized as it is in other large-subunit catalases whose structure has been determined. Different in CAT-3 from these enzymes is an occupied central cavity. Results presented here indicate that CAT-3 and CAT-1 enzymes represent a functional group of catalases with distinctive structural characteristics that determine similar kinetics.     

Lysophosphatidic acid receptor-5 negatively regulates cellular responses in mouse fibroblast 3T3 cells

Lysophosphatidic acid (LPA) signaling via G protein-coupled LPA receptors (LPA₁–LPA₆) mediates a variety of biological functions, including cell migration. Recently, we have reported that LPA₁ inhibited the cell motile activities of mouse fibroblast 3T3 cells. In the present study, to evaluate a role of LPA₅ in cellular responses, Lpar5 knockdown (3T3-L5) cells were generated from 3T3 cells. In cell proliferation assays, LPA markedly stimulated the cell proliferation activities of 3T3-L5 cells, compared with control cells. In cell motility assays with Cell Culture Inserts, the cell motile activities of 3T3-L5 cells were significantly higher than those of control cells. The activity levels of matrix metalloproteinases (MMPs) were measured by gelatin zymography. 3T3-L5 cells stimulated the activation of Mmp-2, correlating with the expression levels of Mmp-2 gene. Moreover, to assess the co-effects of LPA₁ and LPA₅ on cell motile activities, Lpar5 knockdown (3T3a1-L5) cells were also established from Lpar1 over-expressing (3T3a1) cells. 3T3a1-L5 cells increased the cell motile activities of 3T3a1 cells, while the cell motile activities of 3T3a1 cells were significantly lower than those of control cells. These results suggest that LPA₅ may act as a negative regulator of cellular responses in mouse fibroblast 3T3 cells, similar to the case for LPA₁.     

Morphological and Molecular Differentiation of Sporidiobolus johnsonii ATCC 20490 and Its Coenzyme Q10 Overproducing Mutant Strain UF16

Coenzyme Q₁₀ (CoQ₁₀) is an industrially important molecule having nutraceutical and cosmeceutical applications. CoQ₁₀ is mainly produced by microbial fermentation and the process demands the use of strains with high productivity and yields of CoQ₁₀. During strain improvement program consisting of sequential induced mutagenesis, rational selection and screening process, a mutant strain UF16 was generated from Sporidiobolus johnsonii ATCC 20490 with 2.3-fold improvements in CoQ₁₀ content. EMS and UV rays were used as mutagenic agents for generating UF16 and it was rationally selected based on atorvastatin resistance as well as survival at free radicals exposure. We investigated the genotypic and phenotypic changes in UF16 in order to differentiate it from wild type strain. Morphologically it was distinct due to reduced pigmentation of colony, reduced cell size and significant reduction in mycelial growth forms with abundance of yeast forms. At molecular level, UF16 was differentiated based on PCR fingerprinting method of RAPD as well as large and small-subunit rRNA gene sequences. Rapid molecular technique of RAPD analysis using six primers showed 34 % polymorphic fragments with mean genetic distance of 0.235. The partial sequences of rRNA-gene revealed few mutation sites on nucleotide base pairs. However, the mutations detected on rRNA gene of UF16 were less than 1 % of total base pairs and its sequence showed 99 % homology with the wild type strain. These mutations in UF16 could not be linked to phenotypic or genotypic changes on CoQ₁₀ biosynthetic pathway that resulted in improved yield. Hence, investigating the mutations responsible for deregulation of CoQ₁₀ pathway is essential to understand the cause of overproduction in UF16. Phylogenetic analysis based on RAPD bands and rRNA gene sequences coupled with morphological variations, exhibited the novelty of mutant UF16 having potential for improved CoQ₁₀ production.     

C-terminal truncation of a bovine B12 trafficking chaperone enhances the sensitivity of the glutathione-regulated thermostability

The human B₁₂ trafficking chaperone hCblC is well conserved in mammals and non-mammalian eukaryotes. However, the C-terminal ∼40 amino acids of hCblC vary significantly and are predicted to be deleted by alternative splicing of the encoding gene. In this study, we examined the thermostability of the bovine CblC truncated at the C-terminal variable region (t-bCblC) and its regulation by glutathione. t-bCblC is highly thermolabile (T_m = ∼42℃) similar to the full-length protein (f-bCblC). However, t-bCblC is stabilized to a greater extent than f-bCblC by binding of reduced glutathione (GSH) with increased sensitivity to GSH. In addition, binding of oxidized glutathione (GSSG) destabilizes t-bCblC to a greater extent and with increased sensitivity as compared to f-bCblC. These results indicate that t-bCblC is a more sensitive form to be regulated by glutathione than the full-length form of the protein. [BMB Reports 2013; 46(3): 169-174]     

Influence of Histidine-198 of the D1 subunit on the properties of the primary electron donor, P680, of photosystem II in Thermosynechococcus elongatus

The influence of the histidine axial ligand to the P_D1 chlorophyll of photosystem II on the redox potential and spectroscopic properties of the primary electron donor, P_680, was investigated in mutant oxygen-evolving photosystem II (PSII) complexes purified from the thermophilic cyanobacterium Thermosynechococcus elongatus. To achieve this aim, a mutagenesis system was developed in which the psbA₁ and psbA₂ genes encoding D1 were deleted from a His-tagged CP43 strain (to generate strain WT^?) and mutations D1-H198A and D1-H198Q were introduced into the remaining psbA₃ gene. The O₂-evolving activity of His-tagged PSII isolated from WT^? was found to be significantly higher than that measured from His-tagged PSII isolated from WT in which psbA₁ is expected to be the dominantly expressed form. PSII purified from both the D1-H198A and D1-H198Q mutants exhibited oxygen-evolving activity as high as that from WT^?. Surprisingly, a variety of kinetic and spectroscopic measurements revealed that the D1-H198A and D1-H198Q mutations had little effect on the redox and spectroscopic properties of P₆₈₀, in contrast to the earlier results from the analysis of the equivalent mutants constructed in Synechocystis sp. PCC 6803 [B.A. Diner, E. Schlodder, P.J. Nixon, W.J. Coleman, F. Rappaport, J. Lavergne, W.F. Vermaas, D.A. Chisholm, Site-directed mutations at D1-His198 and D2-His197 of photosystem II in Synechocystis PCC 6803: sites of primary charge separation and cation and triplet stabilization, Biochemistry 40 (2001) 9265-9281]. We conclude that the nature of the axial ligand to P_D1 is not an important determinant of the redox and spectroscopic properties of P₆₈₀ in T. elongatus.     

Crystal Structure of the Electron Carrier Domain of the Reaction Center Cytochrome cz Subunit from Green Photosynthetic Bacterium Chlorobium tepidum

In green sulfur photosynthetic bacteria, the cytochrome c_z (cyt c_z) subunit in the reaction center complex mediates electron transfer mainly from menaquinol/cytochrome c oxidoreductase to the special pair (P840) of the reaction center. The cyt c_z subunit consists of an N-terminal transmembrane domain and a C-terminal soluble domain that binds a single heme group. The periplasmic soluble domain has been proposed to be highly mobile and to fluctuate between oxidoreductase and P840 during photosynthetic electron transfer. We have determined the crystal structure of the oxidized form of the C-terminal functional domain of the cyt c_z subunit (C-cyt c_z) from thermophilic green sulfur bacterium Chlorobium tepidum at 1.3-Å resolution. The overall fold of C-cyt c_z consists of four α-helices and is similar to that of class I cytochrome c proteins despite the low similarity in their amino acid sequences. The N-terminal structure of C-cyt c_z supports the swinging mechanism previously proposed in relation with electron transfer, and the surface properties provide useful information on possible interaction sites with its electron transfer partners. Several characteristic features are observed for the heme environment: These include orientation of the axial ligands with respect to the heme plane, surface-exposed area of the heme, positions of water molecules, and hydrogen-bond network involving heme propionate groups. These structural features are essential for elucidating the mechanism for regulating the redox state of cyt c_z.     

Association of 25-Hydroxyvitamin D status and genetic variation in the vitamin D metabolic pathway with FEV1 in the Framingham Heart Study

Background

Vitamin D is associated with lung function in cross-sectional studies, and vitamin D inadequacy is hypothesized to play a role in the pathogenesis of chronic obstructive pulmonary disease. Further data are needed to clarify the relation between vitamin D status, genetic variation in vitamin D metabolic genes, and cross-sectional and longitudinal changes in lung function in healthy adults.

Methods

We estimated the association between serum 25-hydroxyvitamin D [25(OH)D] and cross-sectional forced expiratory volume in the first second (FEV₁) in Framingham Heart Study (FHS) Offspring and Third Generation participants and the association between serum 25(OH)D and longitudinal change in FEV₁ in Third Generation participants using linear mixed-effects models. Using a gene-based approach, we investigated the association between 241 SNPs in 6 select vitamin D metabolic genes in relation to longitudinal change in FEV₁ in Offspring participants and pursued replication of these findings in a meta-analyzed set of 4 independent cohorts.

Results

We found a positive cross-sectional association between 25(OH)D and FEV₁ in FHS Offspring and Third Generation participants (P = 0.004). There was little or no association between 25(OH)D and longitudinal change in FEV₁ in Third Generation participants (P = 0.97). In Offspring participants, the CYP2R1 gene, hypothesized to influence usual serum 25(OH)D status, was associated with longitudinal change in FEV₁ (gene-based P < 0.05). The most significantly associated SNP from CYP2R1 had a consistent direction of association with FEV₁ in the meta-analyzed set of replication cohorts, but the association did not reach statistical significance thresholds (P = 0.09).

Conclusions

Serum 25(OH)D status was associated with cross-sectional FEV₁, but not longitudinal change in FEV₁. The inconsistent associations may be driven by differences in the groups studied. CYP2R1 demonstrated a gene-based association with longitudinal change in FEV₁ and is a promising candidate gene for further studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12931-015-0238-y) contains supplementary material, which is available to authorized users.     

Growth of Arabidopsis seedlings on high fungal doses of Piriformospora indica has little effect on plant performance,stress, and defense gene expression in spite of elevated jasmonic acid and jasmonic acid-isoleucine levels in the roots

The endophytic fungus Piriformospora indica colonizes the roots of many plant species including Arabidopsis and promotes their performance, biomass, and seed production as well as resistance against biotic and abiotic stress. Imbalances in the symbiotic interaction such as uncontrolled fungal growth result in the loss of benefits for the plants and activation of defense responses against the microbe. We exposed Arabidopsis seedlings to a dense hyphal lawn of P. indica. The seedlings continue to grow, accumulate normal amounts of chlorophyll, and the photosynthetic parameters demonstrate that they perform well. In spite of high fungal doses around the roots, the fungal material inside the roots was not significantly higher when compared with roots that live in a beneficial symbiosis with P. indica. Fifteen defense- and stress-related genes including PR2, PR3, PAL2, and ERF1 are only moderately upregulated in the roots on the fungal lawn, and the seedlings did not accumulate H₂O₂/radical oxygen species. However, accumulation of anthocyanin in P. indica-exposed seedlings indicates stress symptoms. Furthermore, the jasmonic acid (JA) and jasmonic acid-isoleucine (JA-Ile) levels were increased in the roots, and consequently PDF1.2 and a newly characterized gene for a 2-oxoglurate and Fe²⁺-dependent oxygenase were upregulated more than 7-fold on the dense fungal lawn, in a JAR1- and EIN3-dependent manner. We conclude that growth of A. thaliana seedlings on high fungal doses of P. indica has little effect on the overall performance of the plants although elevated JA and JA-Ile levels in the roots induce a mild stress or defense response.     

Rapid In Situ Hybridization using Oligonucleotide Probes on Paraformaldehyde-prefixed Brain of Rats with Serotonin Syndrome

3,4-Methylenedioxymethamphetamine (MDMA; ecstasy) toxicity may cause region-specific changes in serotonergic mRNA expression due to acute serotonin (5-hydroxytryptamine; 5-HT) syndrome. This hypothesis can be tested using in situ hybridization to detect the serotonin 5-HT_2A receptor gene htr2a. In the past, such procedures, utilizing radioactive riboprobe, were difficult because of the complicated workflow that needs several days to perform and the added difficulty that the technique required the use of fresh frozen tissues maintained in an RNase-free environment. Recently, the development of short oligonucleotide probes has simplified in situ hybridization procedures and allowed the use of paraformaldehyde-prefixed brain sections, which are more widely available in laboratories. Here, we describe a detailed protocol using non-radioactive oligonucleotide probes on the prefixed brain tissues. Hybridization probes used for this study include dapB (a bacterial gene coding for dihydrodipicolinate reductase), ppiB (a housekeeping gene coding for peptidylprolyl isomerase B), and htr2a (a serotonin gene coding for 5-HT_2A receptors). This method is relatively simply, cheap, reproducible and requires less than two days to complete.     

Moving Iron through ferritin protein nanocages depends on residues throughout each four α-helix bundle subunit

Eukaryotic H ferritins move iron through protein cages to form biologically required, iron mineral concentrates. The biominerals are synthesized during protein-based Fe²⁺/O₂ oxidoreduction and formation of [Fe³⁺O]_n multimers within the protein cage, en route to the cavity, at sites distributed over ∼50 Å. Recent NMR and Co²⁺-protein x-ray diffraction (XRD) studies identified the entire iron path and new metal-protein interactions: (i) lines of metal ions in 8 Fe²⁺ ion entry channels with three-way metal distribution points at channel exits and (ii) interior Fe³⁺O nucleation channels. To obtain functional information on the newly identified metal-protein interactions, we analyzed effects of amino acid substitution on formation of the earliest catalytic intermediate (diferric peroxo-A_{650 nm}) and on mineral growth (Fe³⁺O-A_{350 nm}), in A26S, V42G, D127A, E130A, and T149C. The results show that all of the residues influenced catalysis significantly (p < 0.01), with effects on four functions: (i) Fe²⁺ access/selectivity to the active sites (Glu¹³⁰), (ii) distribution of Fe²⁺ to each of the three active sites near each ion channel (Asp¹²⁷), (iii) product (diferric oxo) release into the Fe³⁺O nucleation channels (Ala²⁶), and (iv) [Fe³⁺O]_n transit through subunits (Val⁴², Thr¹⁴⁹). Synthesis of ferritin biominerals depends on residues along the entire length of H subunits from Fe²⁺ substrate entry at 3-fold cage axes at one subunit end through active sites and nucleation channels, at the other subunit end, inside the cage at 4-fold cage axes. Ferritin subunit-subunit geometry contributes to mineral order and explains the physiological impact of ferritin H and L subunits.     

Synergistic effect of an Escherichia coli mutator gene on mutagenesis by ultraviolet radiation and by alkylating agents

E. coli strains differing in a gene responsible for high spontaneous mutability (mut HI) were compared for their mutability by UV radiation and by the alkylating agents ethyl methanesulfonate and methyl methanesulfonate. All three exogenous mutagenic agents induced significantly higher frequencies of mutants with impaired carbohydrate-fermenting ability when the mutator allele rather than the wild-type allele was present. Thus the mut HI gene product possibly increases the probability of replication error due to alterations in the structure of the template strand of DNA. An attempt to detect an synergistic effect for UV-induced suppressor mutations was unsuccessful. The failure may have been due to the particular method used for scoring this type of mutation.     

A new positive relationship between pCO2 and stomatal frequency in Quercus guyavifolia (Fagaceae): a potential proxy for palaeo-CO2 levels

Background and Aims The inverse relationship between atmospheric CO₂ partial pressure (pCO₂) and stomatal frequency in many species of plants has been widely used to estimate palaeoatmospheric CO₂ (palaeo-CO₂) levels; however, the results obtained have been quite variable. This study attempts to find a potential new proxy for palaeo-CO₂ levels by analysing stomatal frequency in Quercus guyavifolia (Q. guajavifolia, Fagaceae), an extant dominant species of sclerophyllous forests in the Himalayas with abundant fossil relatives.Methods Stomatal frequency was analysed for extant samples of Q. guyavifolia collected from17 field sites at altitudes ranging between 2493 and 4497 m. Herbarium specimens collected between 1926 and 2011 were also examined. Correlations of pCO₂–stomatal frequency were determined using samples from both sources, and these were then applied to Q. preguyavaefolia fossils in order to estimate palaeo-CO₂ concentrations for two late-Pliocene floras in south-western China.Key Results In contrast to the negative correlations detected for most other species that have been studied, a positive correlation between pCO₂ and stomatal frequency was determined in Q. guyavifolia sampled from both extant field collections and historical herbarium specimens. Palaeo-CO₂ concentrations were estimated to be approx. 180–240 ppm in the late Pliocene, which is consistent with most other previous estimates.Conclusions A new positive relationship between pCO₂ and stomatal frequency in Q. guyavifolia is presented, which can be applied to the fossils closely related to this species that are widely distributed in the late-Cenozoic strata in order to estimate palaeo-CO₂ concentrations. The results show that it is valid to use a positive relationship to estimate palaeo-CO₂ concentrations, and the study adds to the variety of stomatal density/index relationships that available for estimating pCO₂. The physiological mechanisms underlying this positive response are unclear, however, and require further research.