Using flow cytometry to estimate pollen DNA content: improved methodology and applications

Background and Aims

Flow cytometry has been used to measure nuclear DNA content in pollen, mostly to understand pollen development and detect unreduced gametes. Published data have not always met the high-quality standards required for some applications, in part due to difficulties inherent in the extraction of nuclei. Here we describe a simple and relatively novel method for extracting pollen nuclei, involving the bursting of pollen through a nylon mesh, compare it with other methods and demonstrate its broad applicability and utility.

Methods

The method was tested across 80 species, 64 genera and 33 families, and the data were evaluated using established criteria for estimating genome size and analysing cell cycle. Filter bursting was directly compared with chopping in five species, yields were compared with published values for sonicated samples, and the method was applied by comparing genome size estimates for leaf and pollen nuclei in six species.

Key Results

Data quality met generally applied standards for estimating genome size in 81 % of species and the higher best practice standards for cell cycle analysis in 51 %. In 41 % of species we met the most stringent criterion of screening 10 000 pollen grains per sample. In direct comparison with two chopping techniques, our method produced better quality histograms with consistently higher nuclei yields, and yields were higher than previously published results for sonication. In three binucleate and three trinucleate species we found that pollen-based genome size estimates differed from leaf tissue estimates by 1·5 % or less when 1C pollen nuclei were used, while estimates from 2C generative nuclei differed from leaf estimates by up to 2·5 %.

Conclusions

The high success rate, ease of use and wide applicability of the filter bursting method show that this method can facilitate the use of pollen for estimating genome size and dramatically improve unreduced pollen production estimation with flow cytometry.     

Properties of crowding indices and statistical tools to analyze parasite crowding data

Crowding, i.e., the size of the infrapopulation inhabiting an individual host, is a major component of parasites' environment, which often influences both morphological and life-history characters (the so-called density-dependent characters) in different parasite taxa. Although crowding equals intensity in case of a single parasite individual, mean intensity of the host population does not define mean crowding of the parasite population. Crowding indices are notoriously hard to handle statistically because of the inherently large number of nonindependent values in data. In this study, we aim to investigate the apparently paradox features of crowding indices and to make some proposals and also to introduce statistical methods to calculate confidence intervals and 1-sample and 2-sample tests for mean crowding. All methods described in this study are supported by the freely distributed statistical software Quantitative Parasitology.     

An improved statistical method for detecting heterotachy in nucleotide sequences

The principle of heterotachy states that the substitution rate of sites in a gene can change through time. In this article, we propose a powerful statistical test to detect sites that evolve according to the process of heterotachy. We apply this test to an alignment of 1289 eukaryotic rRNA molecules to 1) determine how widespread the phenomenon of heterotachy is in ribosomal RNA, 2) to test whether these heterotachous sites are nonrandomly distributed, that is, linked to secondary structure features of ribosomal RNA, and 3) to determine the impact of heterotachous sites on the bootstrap support of monophyletic groupings. Our study revealed that with 21 monophyletic taxa, approximately two-thirds of the sites in the considered set of sequences is heterotachous. Although the detected heterotachous sites do not appear bound to specific structural features of the small subunit rRNA, their presence is shown to have a large beneficial influence on the bootstrap support of monophyletic groups. Using extensive testing, we show that this may not be due to heterotachy itself but merely due to the increased substitution rate at the detected heterotachous sites.     

Respiratory transfer impedances with pressure input at the mouth and chest

Two methods of measuring respiratory transfer impedance (Ztr) were compared in 14 normal subjects, from 4 to 30 Hz, 1) studying the relationship between transrespiratory pressure (Prs) and flow at the chest when varying pressure at the mouth (Ztrm) and 2) studying the relationship between Prs and flow at the mouth when varying pressure around the chest wall (Ztrw). The similarity of the two relationships was expected on the basis of a T-network model. Almost identical phase responses were obtained from the two methods. Pressure-flow ratios were slightly larger for Ztrw than for Ztrm, but differences did not exceed 2% on average in 11 of 14 subjects. When the data were analyzed with the six-coefficient model proposed by DuBois et al. (J. Appl. Physiol. 8: 587-594, 1956), similar values were found for tissue compliance and tissue inertance but slightly different values for gaseous inertance in the airways (1.97 +/- 0.35 X 10(-2) cmH2O X l-1 X s2 for Ztrw vs. 1.73 +/- 0.26 for Ztrm; P less than 0.01). Similar results were also found for total respiratory resistance but with a slightly larger contribution of airway resistance for Ztrw (64 +/- 14 vs. 57 +/- 10%; P less than 0.05). As a practical conclusion it is recommended to measure Ztrw, which is technically much easier.     

An improved mathematical method to estimate DNA synthesis time of bromodeoxyuridine-labelled cells, using FCM-derived data

Abstract. Chinese hamster ovary cells in vitro were pulse-labelled with bromodeoxyuridine (BrdUrd and were then allowed to progress through the cell cycle. Every half hour after labelling, cells were harvested and prepared for simultaneous flow cytometric determination of DNA content and incorporated BrdUrd, with the intercalating dye propidium iodide and with a monoclonal antibody against incorporated BrdUrd, respectively. The relative movement (RM), i.e. the relative mean DNA content of the moving cohort of BrdUrd-labelled cells in relation to that of G₁ and G₂ cells, was calculated. RM was then used to calculate DNA synthesis time (T_S), at all post-labelling times (t). Since labelled cells in G₂ and mitosis (M) in addition to S phase cells, are included in the cohort of moving labelled cells, and since the time of G₂ and M (T_g2+M) phases is finite, a non-linear relationship exists between RM and post-labelling time. Because of this, the use of a linear formula in the calculation of T_S yields results that are affected by t. We found that RM data can be corrected with regard to T_G2+M resulting in the derivation of a non-linear T_S formula. This non-linear T_S formula gave results that were nearly independent of t. Moreover, windows were set in the mid DNA distributions for G₁, S and G₂+ M cells in the bivariate DNA v. BrdUrd cytograms, to estimate the fraction of BrdUrd-labelled cells in each window at every post-labelling time. Plots of the fraction of BrdUrd-labelled cells v. post-labelling time were then made for each window. T_S obtained in this way was in agreement with T_S obtained with the corrected RM method. In conclusion, we present a method to calculate T_s which theoretically first makes the determination of RM independent of T_G2+M, and secondly compensates for the non-linear function of RM with post-labelling time caused by accumulation of BrdUrd-labelled cells in G₂+ M.     

Respiratory impedances and acinar gas transfer in a canine model for emphysema

Barnas, George M., Paul A. Delaney, Ileana Gheorghiu,Srinivas Mandava, Robert G. Russell, Renée Kahn, and Colin F. Mackenzie. Respiratory impedances and acinar gas transfer in acanine model for emphysema. J. Appl.Physiol. 83(1): 179-188, 1997.We examined howthe changes in the acini caused by emphysema affected gas transfer outof the acinus (T_aci) and lungand chest wall mechanical properties. Measurements were taken from fivedogs before and 3 mo after induction of severe bilateral emphysema byexposure to papain aerosol (170-350 mg/dose) for 4 consecutive wk.With the dogs anesthetized, paralyzed, and mechanically ventilated at0.2 Hz and 20 ml/kg, we measuredT_aci by the rate of washout of¹³³Xe from an area of the lungwith occluded blood flow. Measurements were repeated at positiveend-expiratory pressures (PEEP) of 10, 5, 15, 0, and 20 cmH₂O. We also measured dynamicelastances and resistances of the lungs(EL andRL, respectively) and chest wall at the different PEEP and during sinusoidal forcing in the normal rangeof breathing frequency and tidal volume. After final measurements, tissue sections from five randomly selected areas of the lung eachshowed indications of emphysema.T_aci during emphysema was similarto that in control dogs. ELdecreased by ~50% during emphysema (P < 0.05) but did not change itsdependence on frequency or tidal volume.RL did not change(P > 0.05) at the lowest frequencystudied (0.2 Hz), but in some dogs it increased compared with control at the higher frequencies. Chest wall properties were not changed byemphysema (P > 0.05). We suggestthat although large changes in acinar structure andEL occur during uncomplicatedbilateral emphysema, secondary complications must be present to causeseveral of the characteristic dysfunctions seen in patients withemphysema.

     

An improved immunoelectrophoretic methodology for the evaluation of protein caloric malnutrition.

An improved technique of two dimensional immunoelectrophoresis of serum proteins on gelatinized cellulose acetate for routine evaluation of protein nutritional state is described. The method, which allows the estimation of more than 30 protein fractions, proves extremely useful in monitoring the effectiveness of total parenteral nutrition and protein metabolic conditions. Sera samples from healthy subjects and cachectical neoplastic patients have been examined and preliminary results are reported.     

An improved model of heat transfer through penguin feathers and down

Penguins, mostly live in the extremely cold Antarctic, are known to have feathers and down, which are light weight, compact and extremely efficient in preventing heat loss. Nevertheless, the mechanisms of heat transfer through the penguin feathers and down, and how the unique characteristics of penguin feathers and down make them such good thermal insulators are not fully understood. In this paper, an integrated model of heat transfer through the penguin feathers and down is developed and computed using finite volume method, with the geometrical structure of the barbules being considered. Monte-Carlo method is adopted to determine the radiative absorption and emission constant in the integrated model. The effective thermal conductance of penguin feathers and down computed from our model compared well with the experimentally measured value reported in the literature. Three models (penguin model, random fibre model (fibre radius=3microm) and random fibre model (fibre radius=10microm)) are further simulated and compared. Results showed that the relative small radius of the barbules of penguin feather and their geometrical structure are responsible for the reduction of heat loss in cold environment.     

Correlation functions quantify super-resolution images and estimate apparent clustering due to over-counting

We present an analytical method using correlation functions to quantify clustering in super-resolution fluorescence localization images and electron microscopy images of static surfaces in two dimensions. We use this method to quantify how over-counting of labeled molecules contributes to apparent self-clustering and to calculate the effective lateral resolution of an image. This treatment applies to distributions of proteins and lipids in cell membranes, where there is significant interest in using electron microscopy and super-resolution fluorescence localization techniques to probe membrane heterogeneity. When images are quantified using pair auto-correlation functions, the magnitude of apparent clustering arising from over-counting varies inversely with the surface density of labeled molecules and does not depend on the number of times an average molecule is counted. In contrast, we demonstrate that over-counting does not give rise to apparent co-clustering in double label experiments when pair cross-correlation functions are measured. We apply our analytical method to quantify the distribution of the IgE receptor (FcεRI) on the plasma membranes of chemically fixed RBL-2H3 mast cells from images acquired using stochastic optical reconstruction microscopy (STORM/dSTORM) and scanning electron microscopy (SEM). We find that apparent clustering of FcεRI-bound IgE is dominated by over-counting labels on individual complexes when IgE is directly conjugated to organic fluorophores. We verify this observation by measuring pair cross-correlation functions between two distinguishably labeled pools of IgE-FcεRI on the cell surface using both imaging methods. After correcting for over-counting, we observe weak but significant self-clustering of IgE-FcεRI in fluorescence localization measurements, and no residual self-clustering as detected with SEM. We also apply this method to quantify IgE-FcεRI redistribution after deliberate clustering by crosslinking with two distinct trivalent ligands of defined architectures, and we evaluate contributions from both over-counting of labels and redistribution of proteins.     

Merging fluorescence resonance energy transfer and expressed protein ligation to analyze protein-protein interactions

Determination of protein oligomerization state can be technically challenging. We have combined the methods of expressed protein ligation (EPL) and fluorescence resonance energy transfer (FRET) for the analysis of protein homo-oligomerization states. We have attached fluorescein (donor) and rhodamine (acceptor) chromophores via dipeptide linkages to the C-termini of three recombinant proteins and examined the potential for FRET between mixtures of these semisynthetic proteins. The known protein dimer (glutathione S-transferase) showed evidence of FRET and the known protein monomer (SH2 domain phosphatase-1) did not display FRET. Using this method, the previously uncharacterized circadian rhythm enzyme, serotonin N-acetyltransferase, displayed significant FRET, indicating its likely propensity for dimerization or more complex oligomerization. These results establish the potential of the union of EPL and FRET in the analysis of protein-protein interactions and provide insight into the unusual enzymatic behavior of a key circadian rhythm enzyme.     

Effects of posture and bronchoconstriction on low-frequency input and transfer impedances in humans.

We simultaneously evaluated the mechanical response of the total respiratory system, lung, and chest wall to changes in posture and to bronchoconstriction. We synthesized the optimal ventilation waveform (OVW) approach, which simultaneously provides ventilation and multifrequency forcing, with optoelectronic plethysmography (OEP) to measure chest wall flow globally and locally. We applied an OVW containing six frequencies from 0.156 to 4.6 Hz to the mouth of six healthy men in the seated and supine positions, before and after methacholine challenge. We measured mouth, esophageal, and transpulmonary pressures, airway flow by pneumotachometry, and total chest wall, pulmonary rib cage, and abdominal volumes by OEP. We computed total respiratory, lung, and chest wall input impedances and the total and regional transfer impedances (Ztr). These data were appropriately sensitive to changes in posture, showing added resistance in supine vs. seated position. The Ztr were also highly sensitive to lung constriction, more so than input impedance, as the former is minimally distorted by shunting of flow into alveolar gas compression and airway walls. Local impedances show that, during bronchoconstriction and at typical breathing frequencies, the contribution of the abdomen becomes amplified relative to the rib cage. A similar redistribution occurs when passing from seated to supine. These data suggest that the OEP-OVW approach for measuring Ztr could noninvasively track important lung and respiratory conditions, even in subjects who cannot cooperate. Applications might range from routine evaluation of airway hyperreactivity in asthmatic subjects to critical conditions in the supine position during mechanical ventilation.     

DNA bar coding and pyrosequencing to analyze adverse events in therapeutic gene transfer

Gene transfer has been used to correct inherited immunodeficiencies, but in several patients integration of therapeutic retroviral vectors activated proto-oncogenes and caused leukemia. Here, we describe improved methods for characterizing integration site populations from gene transfer studies using DNA bar coding and pyrosequencing. We characterized 160 232 integration site sequences in 28 tissue samples from eight mice, where Rag1 or Artemis deficiencies were corrected by introducing the missing gene with gamma-retroviral or lentiviral vectors. The integration sites were characterized for their genomic distributions, including proximity to proto-oncogenes. Several mice harbored abnormal lymphoproliferations following therapy—in these cases, comparison of the location and frequency of isolation of integration sites across multiple tissues helped clarify the contribution of specific proviruses to the adverse events. We also took advantage of the large number of pyrosequencing reads to show that recovery of integration sites can be highly biased by the use of restriction enzyme cleavage of genomic DNA, which is a limitation in all widely used methods, but describe improved approaches that take advantage of the power of pyrosequencing to overcome this problem. The methods described here should allow integration site populations from human gene therapy to be deeply characterized with spatial and temporal resolution.     

Cross-intensity functions and the estimate of spike-time jitter

Correlation measures are important tools for the analysis of simultaneously recorded spike trains. A well-known measure with probabilistic interpretation is the cross-intensity function (CIF), which is an estimate of the conditional probability that a neuron spikes as a function of the time lag to spikes in another neuron. The non-commutative nature of the CIF is particularly useful when different neuron classes are studied that can be distinguished based on their anatomy or physiology. Here we explore the utility of the CIF for estimating spike-time jitter in synaptic interactions between neuron pairs of connected classes. When applied to spike train pairs from sleeping songbirds, we are able to distinguish fast synaptic interactions mediated primarily by AMPA receptors from slower interactions mediated by NMDA receptors. We also find that spike jitter increases with the time lag between spikes, reflecting the accumulation of noise in neural activity sequences, such as in synfire chains. In conclusion, we demonstrate some new utility of the CIF as a spike-train measure.     

An improved multi-joint EMG-assisted optimization approach to estimate joint and muscle forces in a musculoskeletal model of the lumbar spine

Muscle force partitioning methods and musculoskeletal system simplifications are key modeling issues that can alter outcomes, and thus change conclusions and recommendations addressed to health and safety professionals. A critical modeling concern is the use of single-joint equilibrium to estimate muscle forces and joint loads in a multi-joint system, an unjustified simplification made by most lumbar spine biomechanical models. In the context of common occupational tasks, an EMG-assisted optimization method (EMGAO) is modified in this study to simultaneously account for the equilibrium at all lumbar joints (M-EMGAO). The results of this improved approach were compared to those of its conventional single-joint equivalent (S-EMGAO) counterpart, the latter method being applied to the same lumbar joints but one at a time. Despite identical geometrical configurations and passive contributions used in both models, computed outcomes clearly differed between single- and multi-joint methods, especially at larger trunk flexed postures and during asymmetric lifting. Moreover, muscle forces predicted by L5-S1 single-joint analyses do not maintain mechanical equilibrium at other spine joints crossed by the same muscles. Assuming that the central nervous system does not attempt to balance the external moments one joint at a time and that a given muscle cannot exert different forces at different joints, the proposed multi-joint method represents a substantial improvement over its single-joint counterpart. This improved approach, hence, resolves trunk muscle forces with biological integrity but without compromising mechanical equilibrium at the lumbar joints.