Theoretical analysis of the footprinting experiment

In the footprinting experiment, an end-radiolabeled DNA restriction fragment is subjected to digest by an endonuclease in the presence and absence of a ligand which alters the endonuclease cleavage rate at sites of ligand-DNA contact. The location of these sites, and the strength of the ligand binding, are then deduced from the measured concentrations of the different oligonucleotides produced by the digest. We analyze the experiment in terms of coupled kinetic equations which take into account the cutting rates of endonuclease for sites with ligand present and absent, and the rates of binding and dissociation of the ligand to a site. As long as the ligand concentration remains essentially constant (which occurs, for example, if digest is terminated early enough to assure that all fragments result from single cuts by the endonuclease), the oligonucleotide concentrations reflect only the ligand binding equilibrium constant (ratio of rate constants) and the cutting rates in the presence and absence of ligand. We also show how the measured oligonucleotide concentrations (from, e.g. an autoradiogram) can be used to deduce the ligand equilibrium binding constants for the various sites on the polymer.     

Theoretical Analysis of the Footprinting Experiment

Abstract

In the footprinting experiment, an end-radiolabeled DNA restriction fragment is subjected to digest by an endonuclease in the presence and absence of a ligand which alters the endonuclease cleavage rate at sites of ligand-DNA contact. The location of these sites, and the strength of the ligand binding, are then deduced from the measured concentrations of the different oligonucleotides produced by the digest. We analyze the experiment in terms of coupled kinetic equations which take into account the cutting rates of endonuclease for sites with ligand present and absent, and the rates of binding and dissociation of the ligand to a site. As long as the ligand concentration remains essentially constant (which occurs, for example, if digest is terminated early enough to assure that all fragments result from single cuts by the endonuclease), the oligonucleotide concentrations reflect only the ligand binding equilibrium constant (ratio of rate constants) and the cutting rates in the presence and absence of ligand. We also show how the measured oligonucleotide concentrations (from, e.g. an autoradiogram) can be used to deduce the ligand equilibrium binding constants for the various sites on the polymer.     

The significance of multiple mating in the social wasp Vespula maculifrons

The evolution of the complex societies displayed by social insects depended partly on high relatedness among interacting group members. Therefore, behaviors that depress group relatedness, such as multiple mating by reproductive females (polyandry), are unexpected in social insects. Nevertheless, the queens of several social insect species mate multiply, suggesting that polyandry provides some benefits that counteract the costs. However, few studies have obtained evidence for links between rates of polyandry and fitness in naturally occurring social insect populations. We investigated if polyandry was beneficial in the social wasp Vespula maculifrons. We used genetic markers to estimate queen mate number in V. maculifrons colonies and assessed colony fitness by counting the number of cells that colonies produced. Our results indicated that queen mate number was directly, strongly, and significantly correlated with the number of queen cells produced by colonies. Because V. maculifrons queens are necessarily reared in queen cells, our results demonstrate that high levels of polyandry are associated with colonies capable of producing many new queens. These data are consistent with the explanation that polyandry is adaptive in V. maculifrons because it provides a fitness advantage to queens. Our research may provide a rare example of an association between polyandry and fitness in a natural social insect population and help explain why queens in this taxon mate multiply.     

Gene expression and the evolution of phenotypic diversity in social wasps

Background  

Organisms are capable of developing different phenotypes by altering the genes they express. This phenotypic plasticity provides a means for species to respond effectively to environmental conditions. One of the most dramatic examples of phenotypic plasticity occurs in the highly social hymenopteran insects (ants, social bees, and social wasps), where distinct castes and sexes all arise from the same genes. To elucidate how variation in patterns of gene expression affects phenotypic variation, we conducted a study to simultaneously address the influence of developmental stage, sex, and caste on patterns of gene expression in Vespula wasps. Furthermore, we compared the patterns found in this species to those found in other taxa in order to investigate how variation in gene expression leads to phenotypic evolution.     

Nestmate relatedness and population genetic structure of the Australian social crab spider Diaea ergandros (Araneae: Thomisidae)

We characterized the population genetic structure of the Australian social spider Diaea ergandros using polymorphic allozyme markers. Our main objectives were to understand the social organization of D. ergandros and discern patterns of gene flow across distantly separated geographical areas. Spiders were sampled from nests located within 100 m wide locales, which were distributed within larger 50 km wide regions. Our results indicated that nestmates could have been produced by a single mother and father in 88.9% of D. ergandros nests. The remainder of nests contained spiders that were probably produced by polyandrous females or were immigrants from foreign nests. Nestmate relatedness was relatively high (r = 0.44) and did not differ significantly between the sexes or among juvenile, subadult and adult life stages. We also discovered that D. ergandros populations were highly structured, with significant differentiation detected among locales (FLR = 0.23) and regions (FRT = 0.081). Spiders within locales were also substantially inbred (FIL = 0.15). Overall, our data show that significant population subdivision exists in D. ergandros populations, and we suggest that the poor dispersal ability of Diaea spiders can account for the observed genetic structure.     

Footprinting and circular dichroism studies on paromomycin binding to the packaging region of human immunodeficiency virus type-1

We have studied the interaction of the aminoglycoside drug, paromomycin, with a 171-mer from the packaging region of HIV-1 (psi-RNA), using quantitative footprinting and circular dichroism spectroscopy. The footprinting autoradiographic data were obtained by cutting end-labeled RNA with RNase I or RNase T1 in the presence of varying paromomycin concentrations. Scanning the autoradiograms produced footprinting plots showing cleavage intensities for specific sites on the psi-RNA as functions of drug concentration. Footprinting plots showing binding were analyzed using a two-state model to give apparent binding constants for specific sites of the psi-RNA. These plots show that the highest-affinity paromomycin binding site involves nucleotides near bulges in the main stem and SL-1, and other nucleotides in SL-4 of the psi-RNA. RNase I gives an apparent value of K for this drug site of approximately 1.7 x 10(5) M(-1) while RNase T1 reports a value of K of approximately 8 x 10(4) M(-1) (10 mM Tris HCl, pH 7). Footprinting shows that loading the highest affinity site with paromomycin causes structural changes in the single-stranded linker regions, between the stem-loops and main stem and the loops of SL-1 and SL-3. Drug-induced structural changes also affect the intensity of the 208 nm band in the circular dichroism spectrum of the psi-RNA. Fitting the changes in CD band intensity to a two-state model yielded a binding constant for the highest-affinity drug site of 6 x 10(6) M(-1). Thus, the binding constants from footprinting are lower than those obtained for the highest-affinity site from the circular dichroism spectrum, and lower than those earlier obtained using absorption spectroscopy (Sullivan, J. M.; Goodisman, J.; Dabrowiak, J. C., Bioorg. Med. Chem. Lett. 2002, 12, 615). The discrepancy may be due to competitive binding between drug and cleavage agent in the footprinting experiments, but other explanations are discussed. In addition to revealing sites of binding and regions of drug-induced structural change, footprinting showed that the loop regions of SL-1, SL-3 and SL-4 are exposed in the RNA, whereas the linker region between SL-1 and SL-2 is 'buried' and not accessible to cutting by RNase I or RNase T1.     

CYTONUCLEAR THEORY FOR HAPLODIPLOID SPECIES AND X-LINKED GENES. II. STEPPING-STONE MODELS OF GENE FLOW AND APPLICATION TO A FIRE ANT HYBRID ZONE

We develop cytonuclear, hybrid zone models for haplodiploid species or X-linked genes in diploid species using a stepping-stone framework of migration, in which migration rates vary with both direction and sex. The equilibrium clines for the allele frequencies, cytonuclear disequilibria, and frequencies of pure parental types are examined for species with diagnostic markers, under four important migration schemes: uniform migration of both sexes in both directions, greater migration of both sexes from one direction, greater migration of females, and greater migration of males. Of the three cytonuclear variables examined, the allele frequency clines are the most informative in differentiating among the various migration patterns. The cytonuclear disequilibria and the frequency of the pure parental types tend to be useful only in revealing directional asymmetries in migration. The extent of hybrid zone subdivision has quantitative but not qualitative effects on the distribution of cytonuclear variables, in that the allele frequency clines become more gradual, the cytonuclear disequilibria decrease in magnitude, and the frequencies of pure parentals decline with increasing subpopulation number. Also, the only major difference between the X-linked and haplodiploid frameworks is that a higher frequency of pure parentals is found when considering haplodiploids, in which male production does not require mating. The final important theoretical result is that censusing after migration yields greater disequilibria and parental frequencies than censusing after mating. We analyzed cytonuclear data from two transects from a naturally occurring hybrid zone between two haplodiploid fire ant species, Solenopsis invicta and S. richteri, using our stepping-stone framework. The frequency of S. invicta mtDNA exceeds the frequency of the S. invicta nuclear markers through much of this hybrid zone, indicating that sex differences in migration or selection may be occurring. Maximum-likelihood estimates for the migration rates are very high, due to an unexpectedly large number of pure parental types in the hybrid zone, and differ substantially between the two transects. Overall, our model does not provide a good fit, in part because the S. invicta–S. richteri hybrid zone has not yet reached equilibrium.     

Site-specific binding constants for actinomycin D on DNA determined from footprinting studies.

We report site-specific binding constants for the intercalating anticancer drug actinomycin D (Act-D), binding to a 139-base-pair restriction fragment from pBR 322 DNA. The binding constants are derived from analysis of footprinting experiments, in which the radiolabeled 139-mer is cleaved using DNase I, the cleavage products undergo gel electrophoresis, and, from the gel autoradiogram, spot intensities, proportional to amounts of cleaved fragments, are measured. A bound drug prevents DNase I from cleaving at approximately 7 bonds, leading to decreased amounts of corresponding fragments. With the radiolabel on the 3' end of the noncoding strand (A-label), we measured relative amounts of 54 cleavage products at 25 Act-D concentrations. For cleavage of the 139-mer with the label on the 3' end of the coding strand (G-label), relative amounts of 43 cleavage products at 11 Act-D concentrations were measured. These measurements give information about approximately 120 base pairs of the restriction fragment (approximately 12 turns of the DNA helix); in this region, 14 strong and weak Act-D binding sites were identified. The model used to interpret the footprinting plots is derived in detail. Binding constants for 14 sites on the fragment are obtained simultaneously. It is important to take into account the effect of drug binding at its various sites on the local concentration of probe elsewhere. It is also necessary to include in the model weak as well as strong Act-D sites on the carrier DNA which is present, since the carrier DNA controls the free-drug concentration. As expected, the strongest sites are those with the sequence (all sequences are 5'----3') GC, with TGCT having the highest binding constant, 6.4 x 10(6) M-1. Sites having the sequence GC preceded by G are weak binding sites, having binding constants approximately 1 order of magnitude lower than those of the strong sites. Also, the non-GC-containing sequences CCG and CCC bind Act-D with a binding constant comparable to those of the weak GGC sites. The analysis may reveal drug-induced structural changes on the DNA, which are discussed in terms of the mechanism of Act-D binding.     

Epigenetic inheritance and genome regulation: is DNA methylation linked to ploidy in haplodiploid insects?

Organisms show great variation in ploidy level. For example, chromosome copy number varies among cells, individuals and species. One particularly widespread example of ploidy variation is found in haplodiploid taxa, wherein males are typically haploid and females are typically diploid. Despite the prevalence of haplodiploidy, the regulatory consequences of having separate haploid and diploid genomes are poorly understood. In particular, it remains unknown whether epigenetic mechanisms contribute to regulatory compensation for genome dosage. To gain greater insights into the importance of epigenetic information to ploidy compensation, we examined DNA methylation differences among diploid queen, diploid worker, haploid male and diploid male Solenopsis invicta fire ants. Surprisingly, we found that morphologically dissimilar diploid males, queens and workers were more similar to one another in terms of DNA methylation than were morphologically similar haploid and diploid males. Moreover, methylation level was positively associated with gene expression for genes that were differentially methylated in haploid and diploid castes. These data demonstrate that intragenic DNA methylation levels differ among individuals of distinct ploidy and are positively associated with levels of gene expression. Thus, these results suggest that epigenetic information may be linked to ploidy compensation in haplodiploid insects. Overall, this study suggests that epigenetic mechanisms may be important to maintaining appropriate patterns of gene regulation in biological systems that differ in genome copy number.