The nucleotide changes governing cuticular hydrocarbon variation and their evolution in Drosophila melanogaster

The cuticular hydrocarbon (CH) pheromones in Drosophila melanogaster exhibit strong geographic variation. African and Caribbean populations have a high ratio of 5,9 heptacosadiene/7,11 heptacosadiene (the “High” CH type), whereas populations from all other areas have a low ratio (“Low” CH type). Based on previous genetic mapping, DNA markers were developed that localized the genetic basis of this CH polymorphism to within a 13-kb region. We then carried out a hierarchical search for diagnostic nucleotide sites starting with four lines, and increasing to 24 and 43 lines from a worldwide collection. Within the 13-kb region, only one variable site shows a complete concordance with the CH phenotype. This is a 16-bp deletion in the 5′ region of a desaturase gene (desat2) that was recently suggested to be responsible for the CH polymorphism on the basis of its expression [Dallerac, R., Labeur, C., Jallon, J.-M., Knipple, D. C., Roelofs, W. L. & Wicker-Thomas, C. (2000) Proc. Natl. Acad. Sci. 97, 9449–9454]. The cosmopolitan Low type is derived from the ancestral High type, and DNA sequence variations suggest that the former spread worldwide with the aid of positive selection. Whether this CH variation could be a component of the sexual isolation between Zimbabwe and other cosmopolitan populations remains an interesting and unresolved question.

Coordinating DNA replication initiation with cell growth: differential roles for DnaA and SeqA proteins.

We describe here the development of a new approach to the analysis of Escherichia coli replication control. Cells were grown at low growth rates, in which case the bacterial cell cycle approximates that of eukaryotic cells with G1, S, and G2 phases: cell division is followed sequentially by a gap period without DNA replication, replication of the single chromosome, another gap period, and finally the next cell division. Flow cytometry of such slowly growing cells reveals the timing of replication initiation as a function of cell mass. The data show that initiation is normally coupled to cell physiology extremely tightly: the distribution of individual cell masses at the time of initiation in wild-type cells is very narrow, with a coefficient of variation of less than 9%. Furthermore, a comparison between wild-type and seqA mutant cells shows that initiation occurs at a 10-20% lower mass in the seqA mutant, providing direct evidence that SeqA is a bona fide negative regulator of replication initiation. In dnaA (Ts) mutants the opposite is found: the mass at initiation is dramatically increased and the variability in cell mass at initiation is much higher than that for wild-type cells. In contrast to wild-type and dnaA(Ts) cells, seqA mutant cells frequently go through two initiation events per cell division cycle, and all the origins present in each cell are not initiated in synchrony. The implications for the complex interplay amongst growth, cell division, and DNA replication are discussed.

DNA repeats identify novel virulence genes in Haemophilus influenzae.

The whole genome sequence (1.83 Mbp) of Haemophilus influenzae strain Rd was searched to identify tandem oligonucleotide repeat sequences. Loss or gain of one or more nucleotide repeats through a recombination-independent slippage mechanism is known to mediate phase variation of surface molecules of pathogenic bacteria, including H. influenzae. This facilitates evasion of host defenses and adaptation to the varying microenvironments of the host. We reasoned that iterative nucleotides could identify novel genes relevant to microbe-host interactions. Our search of the Rd genome sequence identified 9 novel loci with multiple (range 6-36, mean 22) tandem tetranucleotide repeats. All were found to be located within putative open reading frames and included homologues of hemoglobin-binding proteins of Neisseria, a glycosyltransferase (IgtC gene product) of Neisseria, and an adhesin of Yersinia. These tetranucleotide repeat sequences were also shown to be present in two other epidemiologically different H. influenzae type b strains, although the number and distribution of repeats was different. Further characterization of the IgtC gene showed that it was involved in phenotypic switching of a lipopolysaccharide epitope and that this variable expression was associated with changes in the number of tetranucleotide repeats. Mutation of IgtC resulted in attenuated virulence of H. influenzae in an infant rat model of invasive infection. These data indicate the rapidity, economy, and completeness with which whole genome sequences can be used to investigate the biology of pathogenic bacteria.

Fos and Jun bend the AP-1 site: effects of probe geometry on the detection of protein-induced DNA bending.

The effect of Fos and Jun binding on the structure of the AP-1 recognition site is controversial. Results from phasing analysis and phase-sensitive detection studies of DNA bending by Fos and Jun have led to opposite conclusions. The differences between these assays, the length of the spacer between two bends and the length of the sequences flanking the bends, are investigated here using intrinsic DNA bend standards. Both an increase in the spacer length as well as a decrease in the length of flanking sequences resulted in a reduction in the phase-dependent variation in electrophoretic mobilities. Probes with a wide separation between the bends and short flanking sequences, such as those used in the phase-sensitive detection studies, displayed no phase-dependent mobility variation. This shape-dependent variation in electrophoretic mobilities was reproduced by complexes formed by truncated Fos and Jun. Results from ligase-catalyzed cyclization experiments have been interpreted to indicate the absence of DNA bending in the Fos-Jun-AP-1 complex. However, truncated Fos and Jun can alter the relative rates of inter- and intramolecular ligation through mechanisms unrelated to DNA bending, confounding the interpretation of cyclization data. The analogous phase- and shape-dependence of the electrophoretic mobilities of the Fos-Jun-AP-1 complex and an intrinsic DNA bend confirm that Fos and Jun bend DNA, which may contribute to their functions in transcription regulation.

Expression of a recoded nuclear gene inserted into yeast mitochondrial DNA is limited by mRNA-specific translational activation.

Genetic code differences prevent expression of nuclear genes within Saccharomyces cerevisiae mitochondria. To bridge this gap a synthetic gene, ARG8m, designed to specify an arginine biosynthetic enzyme when expressed inside mitochondria, has been inserted into yeast mtDNA in place of the COX3 structural gene. This mitochondrial cox3::ARG8m gene fully complements a nuclear arg8 deletion at the level of cell growth, and it is dependent for expression upon nuclear genes that encode subunits of the COX3 mRNA-specific translational activator. Thus, cox3::ARG8m serves as a mitochondrial reporter gene. Measurement of cox3::ARG8m expression at the levels of steady-state protein and enzymatic activity reveals that glucose repression operates within mitochondria. The levels of this reporter vary among strains whose nuclear genotypes lead to under- and overexpression of translational activator subunits, in particular Pet494p, indicating that mRNA-specific translational activation is a rate-limiting step in this organellar system. Whereas the steady-state level of cox3::ARG8m mRNA was also glucose repressed in an otherwise wild-type strain, absence of translational activation led to essentially repressed mRNA levels even under derepressing growth conditions. Thus, the mRNA is stabilized by translational activation, and variation in its level may be largely due to modulation of translation.

A general additive distance with time-reversibility and rate variation among nucleotide sites.

As additivity is a very useful property for a distance measure, a general additive distance is proposed under the stationary time-reversible (SR) model of nucleotide substitution or, more generally, under the stationary, time-reversible, and rate variable (SRV) model, which allows rate variation among nucleotide sites. A method for estimating the mean distance and the sampling variance is developed. In addition, a method is developed for estimating the variance-covariance matrix of distances, which is useful for the statistical test of phylogenies and molecular clocks. Computer simulation shows (i) if the sequences are longer than, say, 1000 bp, the SR method is preferable to simpler methods; (ii) the SR method is robust against deviations from time-reversibility; (iii) when the rate varies among sites, the SRV method is much better than the SR method because the distance is seriously underestimated by the SR method; and (iv) our method for estimating the sampling variance is accurate for sequences longer than 500 bp. Finally, a test is constructed for testing whether DNA evolution follows a general Markovian model.

Changing nucleosome positions through modification of the DNA rotational information.

The effects of the rotational information of DNA in determining the in vitro localization of nucleosomal core particles (ncps) have been studied in the Saccharomyces cerevisiae 5S rRNA repeat gene. We have altered the distribution of the phased series of flexibility signals present on this DNA by inserting a 25-bp tract, and we have analyzed the effects of this mutation on the distribution and on the frequencies of ncps, as compared with the wild type and a reference 21-bp insertion mutant. The variation of the standard free energy of nucleosome reconstitution was determined. The results show that the DNA rotational information is a major determinant of ncps positioning, define how many rotationally phased signals are required for the formation of a stable particle, and teach how to modify their distribution through the alteration of the rotational signals.

Antisense DNA downregulation of the ERBB2 oncogene measured by a flow cytometric assay.

A causal role has been inferred for ERBB2 overexpression in the etiology of breast cancer and other epithelial malignancies. The development of therapeutics that inhibit this tyrosine kinase cell surface receptor remains a high priority. This report describes the specific downregulation of ERBB2 protein and mRNA in the breast cancer cell line SK-BR-3 by using antisense DNA phosphorothioates. An approach was developed to examine antisense effects which allows simultaneous measurements of antisense dose and gene specific regulation on a per cell basis. A fluorescein isothiocyanate end-labeled tracer oligonucleotide was codelivered with antisense DNA followed by immunofluorescent staining for ERBB2 protein expression. Two-color flow cytometry measured the amount of both intracellular oligonucleotide and ERBB2 protein. In addition, populations of cells that received various doses of nucleic acids were physically separated and studied. In any given transfection, a 100-fold variation in oligonucleotide dosage was found. ERBB2 protein expression was reduced greater than 50%, but only in cells within a relatively narrow uptake range. Steady-state ERBB2 mRNA levels were selectively diminished, indicating a specific antisense effect. Cells receiving the optimal antisense dose were sorted and analyzed for cell cycle changes. After 2 days of ERBB2 suppression, breast cancer cells showed an accumulation in the G1 phase of the cell cycle.

Complete congruence between gene diversity estimates derived from genotypic data at enzyme and random amplified polymorphic DNA loci in black spruce.

Controversy still exists over the adaptive nature of variation of enzyme loci. In conifers, random amplified polymorphic DNAs (RAPDs) represent a class of marker loci that is unlikely to fall within or be strongly linked to coding DNA. We have compared the genetic diversity in natural populations of black spruce [Picea mariana (Mill.) B.S.P.] using genotypic data at allozyme loci and RAPD loci as well as phenotypic data from inferred RAPD fingerprints. The genotypic data for both allozymes and RAPDs were obtained from at least six haploid megagametophytes for each of 75 sexually mature individuals distributed in five populations. Heterozygosities and population fixation indices were in complete agreement between allozyme loci and RAPD loci. In black spruce, it is more likely that the similar levels of variation detected at both enzyme and RAPD loci are due to such evolutionary forces as migration and the mating system, rather than to balancing selection and overdominance. Furthermore, we show that biased estimates of expected heterozygosity and among-population differentiation are obtained when using allele frequencies derived from dominant RAPD phenotypes.

On the magnitude of the electrostatic contribution to ligand-DNA interactions.

A model based on the nonlinear Poisson-Boltzmann equation is used to study the electrostatic contribution to the binding free energy of a simple intercalating ligand, 3,8-diamino-6-phenylphenanthridine, to DNA. We find that the nonlinear Poisson-Boltzmann model accurately describes both the absolute magnitude of the pKa shift of 3,8-diamino-6-phenylphenanthridine observed upon intercalation and its variation with bulk salt concentration. Since the pKa shift is directly related to the total electrostatic binding free energy of the charged and neutral forms of the ligand, the accuracy of the calculations implies that the electrostatic contributions to binding are accurately predicted as well. Based on our results, we have developed a general physical description of the electrostatic contribution to ligand-DNA binding in which the electrostatic binding free energy is described as a balance between the coulombic attraction of a ligand to DNA and the disruption of solvent upon binding. Long-range coulombic forces associated with highly charged nucleic acids provide a strong driving force for the interaction of cationic ligands with DNA. These favorable electrostatic interactions are, however, largely compensated for by unfavorable changes in the solvation of both the ligand and the DNA upon binding. The formation of a ligand-DNA complex removes both charged and polar groups at the binding interface from pure solvent while it displaces salt from around the nucleic acid. As a result, the total electrostatic binding free energy is quite small. Consequently, nonpolar interactions, such as tight packing and hydrophobic forces, must play a significant role in ligand-DNA stability.