Integrating Bayesian genomic cline analyses and association mapping of morphological and ecological traits to dissect reproductive isolation and introgression in a Louisiana Iris hybrid zone

Hybrid zones provide unique opportunities to examine reproductive isolation and introgression in nature. We utilized 45,384 single nucleotide polymorphism (SNP) loci to perform association mapping of 14 floral, vegetative and ecological traits that differ between Iris hexagona and Iris fulva, and to investigate, using a Bayesian genomic cline (BGC) framework, patterns of genomic introgression in a large and phenotypically diverse hybrid zone in southern Louisiana. Many loci of small effect size were consistently found to be associated with phenotypic variation across all traits, and several individual loci were revealed to influence phenotypic variation across multiple traits. Patterns of genomic introgression were quite heterogeneous throughout the Louisiana Iris genome, with I. hexagona alleles tending to be favoured over those of I. fulva. Loci that were found to have exceptional patterns of introgression were also found to be significantly associated with phenotypic variation in a small number of morphological traits. However, this was the exception rather than the rule, as most loci that were associated with morphological trait variation were not significantly associated with excess ancestry. These findings provide insights into the complexity of the genomic architecture of phenotypic differences and are a first step towards identifying loci that are associated with both trait variation and reproductive isolation in nature.     

Methylation of lysine residues of histone fractions in synchronized mommalian cells

Synchronized cultures of mammalian cells were labeled with ¹⁴C-methyl methionine. Labeled methionine methyl groups were incorporated into certain histone fractions, forming methyl lysine. Incorporation of labeled methyl group into histone fractions as ¹⁴C-methyl lysine was followed through the cell cycle from late G₁ into early M. The ¹⁴C-methyl lysine contents of fractions F2a and F3 began to rise in S and reached maxima after termination of DNA and histone synthesis, coincident with the beginning of mitosis, and began to fall by mid-M. The ¹⁴C-methyl lysine content of fraction F2b rose to a maximum early in S, coincident with initiation of DNA synthesis, and rapidly decreased to its original unmethylated level by late S. Fraction F1 remained unmethylated during the period G₁-M. Evidence is presented to demonstrate differential methylation of histone fractions and to substantiate differential temporal coupling of the methylation of specific histone fractions with histone and DNA biosynthesis.     

Primary sequence and structural analysis of sterol carrier protein 2 from rat liver: homology with immunoglobulins

Sterol carrier protein 2 (SCP2) is involved in the later steps of cholesterol biosynthesis and in the intracellular transport of cholesterol. In the present investigation, the amino acid sequence of SCP2 from rat liver has been determined. It is a single polypeptide chain with 122 amino acid residues. Secondary structure prediction indicates an amphipathic alpha-helix region for residues 21-34 and antiparallel beta-sheet structure for residues 35-95. A major finding is the significant homology which exists over approximately 80 residues between SCP2 and the variable domains of the heavy chain of immunoglobulin G.     

Recovery of dynamic muscular endurance

Recovery of the rate of dynamic muscular endurance was measured in two groups of college-aged males. Subjects were required to perform elbow flexion (between the angles of 70 and 170 degrees) for as long as possible at the rate of 38 contractions/min while loaded with 1/6 of their maximum isometric strength (MVC). The task was terminated when the subject fell four contractions behind the required cadence or failed to complete two successive contractions. Subsequent to the task the subject was given a predetermined rest period after which a second fatigue bout to failure was performed. The rest intervals for Gp I (n = 22) were 5, 15, 45, 135, 405, and 1215 seconds, while the rest intervals for Gp II (n = 17) were 10, 30, 90, 270, 810, and 2550 s. Each subject completed six recovery intervals with the order of administration assigned at random. The percentage of recovery was calculated by dividing the exercise time of the first bout into the time of the second bout. These normalized data for the two groups were combined for analysis providing a 12 point recovery curve. The percentage of recovery ranged from 15.4% after 5 s to 91.8% after 2550 s. Analysis of the data revealed that the recovery pattern of dynamic muscular endurance progressed very rapidly initially, reached 50% at approximately 2 min and 15 s and was slightly less than 90% complete at 20 min. Exponential analysis of these data yielded a three-component curve.     

THE GENETIC ARCHITECTURE OF REPRODUCTIVE ISOLATION IN LOUISIANA IRISES: HYBRID FITNESS IN NATURE

Negative epistasis in hybrid genomes commonly results in postzygotic isolation between divergent lineages. However, some genomic regions may be selectively neutral or adaptive in hybrids and thus may potentially cross species barriers. We examined postzygotic isolation between ecologically similar species of Louisiana Iris: Iris brevicaulis and I. fulva to determine the potential for adaptive introgression in nature. Line-cross analyses allowed us a general overview of the gene action responsible for fitness-related traits. We then used a QTL mapping approach to detect genomic regions responsible for variation in these traits. Although hybrid classes suffered reduced fitness for many traits, hybrid means were equivalent to at least one of the parental species in overall estimates of maternal and paternal fitness during the two years of the field study. The genetic architecture underlying the fitness-related traits varied across field site and year of the study, thus emphasizing the importance of the environment in determining the degree of postzygotic isolation and potential for introgression across natural hybrid zones.     

Folding, stability, and physical properties of the alpha subunit of bacterial luciferase

Bacterial luciferase is a heterodimeric (alphabeta) enzyme composed of homologous subunits. When the Vibrio harveyi luxA gene is expressed in Escherichia coli, the alpha subunit accumulates to high levels. The alpha subunit has a well-defined near-UV circular dichroism spectrum and a higher intrinsic fluorescence than the heterodimer, demonstrating fluorescence quenching in the enzyme which is reduced in the free subunit [Sinclair, J. F., Waddle, J. J., Waddill, W. F., and Baldwin, T. O. (1993) Biochemistry 32, 5036-5044]. Analytical ultracentrifugation of the alpha subunit has revealed a reversible monomer to dimer equilibrium with a dissociation constant of 14.9 +/- 4.0 microM at 18 degrees C in 50 mM phosphate and 100 mM NaCl, pH 7.0. The alpha subunit unfolded and refolded reversibly in urea-containing buffers by a three-state mechanism. The first transition occurred over the range of 0-2 M urea with an associated free-energy change of 2.24 +/- 0.25 kcal/mol at 18 degrees C in 50 mM phosphate buffer, pH 7.0. The second, occurring between 2.5 and 3.5 M urea, comprised a cooperative transition with a free-energy change of 6.50 +/- 0.75 kcal/mol. The intermediate species, populated maximally at ca. 2 M urea, has defined near-UV circular dichroism spectral properties distinct from either the native or the denatured states. The intrinsic fluorescence of the intermediate suggested that, although the quantum yield had decreased, the tryptophanyl residues remained largely buried. The far-UV circular dichroism spectrum of the intermediate indicated that it had lost ca. 40% of its native secondary structure. N-Terminal sequencing of the products of limited proteolysis of the intermediate showed that the C-terminal region of the alpha subunit became protease labile over the urea concentration range at which the intermediate was maximally populated. These observations have led us to propose an unfolding model in which the first transition is the unfolding of a C-terminal subdomain and the second transition represents the unfolding of a more stable N-terminal subdomain. Comparison of the structural properties of the unfolding intermediate using spectroscopic probes and limited proteolysis of the alpha subunit with those of the alphabeta heterodimer suggested that the unfolding pathway of the alpha subunit is the same, whether it is in the form of the free subunit or in the heterodimer.     

Demonstration of two independently folding domains in the alpha subunit of bacterial luciferase by preferential ligand binding-induced stabilization

The alpha subunit of bacterial luciferase unfolds and refolds reversibly by a three-state mechanism in urea-containing buffer. It has been proposed that the three-state unfolding of the alpha subunit arises from a stepwise unfolding of a C-terminal folding domain at lower concentrations of urea, followed by unfolding of the N-terminal domain at higher concentrations of urea (Noland, B. W., Dangott, L. J., and Baldwin, T. O. (1999) Biochemistry 38, 16136-16145). The location of an anion binding site in the proposed N-terminal folding domain allowed the folding mechanism to be probed in the context of the intact polypeptide. Anions preferentially stabilized the N-terminal domain in a concentration-dependent manner. The polyvalent anions sulfate and phosphate were found to be more stabilizing than monovalent chloride ion. Cations did not show a similar stabilizing effect, demonstrating that the stabilization was due to the anions alone. The purified N-terminal domain prepared by limited proteolysis and anion exchange chromatography was found to refold cooperatively with a midpoint approximately that of the second unfolding transition of the alpha subunit. Phosphate ion stabilized this fragment to roughly the same extent as it did the alpha subunit. The results presented are consistent with the proposed two-domain folding model and demonstrate that anion binding to the N-terminal folding domain stabilizes the alpha subunit of bacterial luciferase.     

Activation of cyclic AMP-dependent protein kinase and stimulation of protein phosphorylation in response to adenosine in C-1300 murine neuroblastoma

DEAE-cellulose chromatography of the 20,000g supernatant fraction of homogenates of C-1300 murine neuroblastoma (clone N2a) yields one major and two minor peaks of cyclic AMP-dependent protein kinase activity. Assessment of the endogenous activation state of the enzyme(s) reveals that the enzyme is fully activated by the treatment of whole cells with adenosine (10 μM) in the presence of the phosphodiesterase inhibitor Ro 20 1724 (0.7 mM). This treatment produces a large elevation in the cyclic AMP content of the cells. The treatment of whole cells with adenosine alone (1–100 μM) or Ro 20 1724 alone (0.1–0.7 mM) produces minimal elevations in cyclic AMP but nevertheless causes significant activations of cyclic AMP-dependent protein kinase. The autophosphorylation of whole homogenates of treated and untreated cells was studied using [γ-³²P] ATP, sodium dodecyl sulfate-polyacrylamide gel electrophoresis and autoradiography. Treatments which activate cyclic AMP-dependent protein kinase selectively stimulate the incorporation of ³²P into several proteins. This stimulation is most prominent in the 15,000-dalton protein band. The addition of cyclic AMP to phosphorylation reactions containing homogenate of untreated cells stimulates the phosphorylation of the same protein bands. These results indicate that adenosine may have regulatory functions through its effect on the cyclic AMP: cyclic AMP-dependent protein kinase system.