Biodiversity of Costa Rican salamanders: Implications of high levels of genetic differentiation and phylogeographic structure for species formation

Although salamanders are characteristic amphibians in Holarctic temperate habitats, in tropical regions they have diversified evolutionarily only in tropical America. An adaptive radiation centered in Middle America occurred late in the history of a single clade, the supergenus Bolitoglossa (Plethodontidae), and large numbers of species now occur in diverse habitats. Sublineages within this clade decrease in number from the northern to southern parts of Middle America, and in Costa Rica, there are but three. Despite this phylogenetic constraint, Costa Rica has many species; the number of salamander species on one local elevational transect in the Cordillera de Talamanca may be the largest for any such transect in the world. Extraordinary variation in sequences of the mitochondrial gene cytochrome b within a clade of the genus Bolitoglossa in Costa Rica reveals strong phylogeographic structure within a single species, Bolitoglossa pesrubra. Allozymic variation in 19 proteins reveals a pattern largely concordant with the mitochondrial DNA phylogeography. More species exist than are currently recognized. Diversification occurs in restricted geographic areas and involves sharp geographic and elevational differentiation and zonation. In their degree of genetic differentiation at a local scale, these species of the deep tropics exceed the known variation of extratropical salamanders, which also differ in being less restricted in elevational range. Salamanders display “tropicality” in that although speciose, they are usually local in distribution and rare. They display strong ecological and physiological differentiation that may contribute importantly to morphological divergence and species formation.

Population genetics of Ice Age brown bears

The Pleistocene was a dynamic period for Holarctic mammal species, complicated by episodes of glaciation, local extinctions, and intercontinental migration. The genetic consequences of these events are difficult to resolve from the study of present-day populations. To provide a direct view of population genetics in the late Pleistocene, we measured mitochondrial DNA sequence variation in seven permafrost-preserved brown bear (Ursus arctos) specimens, dated from 14,000 to 42,000 years ago. Approximately 36,000 years ago, the Beringian brown bear population had a higher genetic diversity than any extant North American population, but by 15,000 years ago genetic diversity appears similar to the modern day. The older, genetically diverse, Beringian population contained sequences from three clades now restricted to local regions within North America, indicating that current phylogeographic patterns may provide misleading data for evolutionary studies and conservation management. The late Pleistocene phylogeographic data also indicate possible colonization routes to areas south of the Cordilleran ice sheet.

The ATPase and protease domains of yeast mitochondrial Lon: Roles in proteolysis and respiration-dependent growth

The ATP-dependent Lon protease of Saccharomyces cerevisiae mitochondria is required for selective proteolysis in the matrix, maintenance of mitochondrial DNA, and respiration-dependent growth. Lon may also possess a chaperone-like function that facilitates protein degradation and protein-complex assembly. To understand the influence of Lon’s ATPase and protease activities on these functions, we examined several Lon mutants for their ability to complement defects of Lon-deleted yeast cells. We also developed a rapid procedure for purifying yeast Lon to homogeneity to study the enzyme’s activities and oligomeric state. A point mutation in either the ATPase or the protease site strongly inhibited the corresponding activity of the pure protein but did not alter the protein’s oligomerization; when expressed at normal low levels neither of these mutant enzymes supported respiration-dependent growth of Lon-deleted cells. When the ATPase- or the protease-containing regions of Lon were expressed as separate truncated proteins, neither could support respiration-dependent growth of Lon-deleted cells; however, coexpression of these two separated regions sustained wild-type growth. These results suggest that yeast Lon contains two catalytic domains that can interact with one another even as separate proteins, and that both are essential for the different functions of Lon.

A knob-associated tandem repeat in maize capable of forming fold-back DNA segments: Are chromosome knobs megatransposons?

A class of tandemly repeated DNA sequences (TR-1) of 350-bp unit length was isolated from the knob DNA of chromosome 9 of Zea mays L. Comparative fluorescence in situ hybridization revealed that TR-1 elements are also present in cytologically detectable knobs on other maize chromosomes in different proportions relative to the previously described 180-bp repeats. At least one knob on chromosome 4 is composed predominantly of the TR-1 repeat. In addition, several small clusters of the TR-1 and 180-bp repeats have been found in different chromosomes, some not located in obvious knob heterochromatin. Variation in restriction fragment fingerprints and copy number of the TR-1 elements was found among maize lines and among maize chromosomes. TR-1 tandem arrays up to 70 kilobases in length can be interspersed with stretches of 180-bp tandem repeat arrays. DNA sequence analysis and restriction mapping of one particular stretch of tandemly arranged TR-1 units indicate that these elements may be organized in the form of fold-back DNA segments. The TR-1 repeat shares two short segments of homology with the 180-bp repeat. The longest of these segments (31 bp; 64% identity) corresponds to the conserved region among 180-bp repeats. The polymorphism and complex structure of knob DNA suggest that, similar to the fold-back DNA-containing giant transposons in Drosophila, maize knob DNA may have some properties of transposable elements.

Ash1p is a site-specific DNA-binding protein that actively represses transcription

ASH1 encodes a protein that is localized specifically to the daughter cell nucleus, where it has been proposed to repress transcription of the HO gene. Using Ash1p purified from baculovirus-infected insect cells, we have shown that Ash1p binds specific DNA sequences in the HO promoter. DNase I protection analyses showed that Ash1p recognizes a consensus sequence, YTGAT. Mutation of this consensus abolishes Ash1p DNA binding in vitro. We have shown that Ash1p requires an intact zinc-binding domain in its C terminus for repression of HO in vivo and that this domain may be involved in DNA binding. A heterologous DNA-binding domain fused to an N-terminal segment of Ash1p functions as an active repressor of transcription. Our studies indicate that Ash1p is a DNA-binding protein of the GATA family with a separable transcriptional repression domain.

Target DNA chromatinization modulates nicking by L1 endonuclease

L1 elements are human transposons which replicate via an RNA intermediate. At least 15% of the human genome is composed of L1 sequence. An important initial step in the transposition reaction is nicking of the genomic DNA by L1 endonuclease (L1 EN). In vivo much of the genome exists in the form of chromatin or is undergoing biochemical transactions such as transcription, replication or repair, which may alter the accessibility of the L1 transposition machinery to DNA. To investigate this possibility we have examined the effect of substrate chromatinization on the ability of L1 EN to nick DNA. We find that DNA incorporated into nucleosomes is generally refractory to nicking by L1 EN. Interestingly, nicking of a minority of DNA sequences is enhanced when included in chromatin. Thus, dynamic epigenetic factors such as chromatinization are likely to influence the relatively permanent placement of L1 and other retroelements in the human genome.

The human mitochondrial deoxynucleotide carrier and its role in the toxicity of nucleoside antivirals

The synthesis of DNA in mitochondria requires the uptake of deoxynucleotides into the matrix of the organelle. We have characterized a human cDNA encoding a member of the family of mitochondrial carriers. The protein has been overexpressed in bacteria and reconstituted into phospholipid vesicles where it catalyzed the transport of all four deoxy (d) NDPs, and, less efficiently, the corresponding dNTPs, in exchange for dNDPs, ADP, or ATP. It did not transport dNMPs, NMPs, deoxynucleosides, nucleosides, purines, or pyrimidines. The physiological role of this deoxynucleotide carrier is probably to supply deoxynucleotides to the mitochondrial matrix for conversion to triphosphates and incorporation into mitochondrial DNA. The protein is expressed in all human tissues that were examined except for placenta, in accord with such a central role. The deoxynucleotide carrier also transports dideoxynucleotides efficiently. It is likely to be medically important by providing the means of uptake into mitochondria of nucleoside analogs, leading to the mitochondrial impairment that underlies the toxic side effects of such drugs in the treatment of viral illnesses, including AIDS, and in cancer therapy.

Induced fit DNA recognition by a minor groove binding analogue of Hoechst 33258: fluctuations in DNA A tract structure investigated by NMR and molecular dynamics simulations

NMR analysis and molecular dynamics simulations of d(GGTAATTACC)2 and its complex with a tetrahydropyrimidinium analogue of Hoechst 33258 suggest that DNA minor groove recognition in solution involves a combination of conformational selection and induced fit, rather than binding to a preorganised site. Analysis of structural fluctuations in the bound and unbound states suggests that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the minor groove resulting in groove narrowing through: (i) changes in base step parameters, including increased helical twist and propeller twist; (ii) changes to the sugar–phosphate backbone conformation to engulf the bound ligand; (iii) suppression of bending modes at the TpA steps. In contrast, the geometrical arrangement of hydrogen bond acceptors on the groove floor appears to be relatively insensitive to DNA conformation (helical twist and propeller twist). We suggest that effective recognition of DNA sequences (in this case an A tract structure) appears to depend to a significant extent on the sequence being flexible enough to be able to adopt the geometrically optimal conformation compatible with the various binding interactions, rather than involving ‘lock and key’ recognition.

Visualization of oligonucleotide probes and point mutations in interphase nuclei and DNA fibers using rolling circle DNA amplification

Rolling circle amplification (RCA) is a surface-anchored DNA replication reaction that can be exploited to visualize single molecular recognition events. Here we report the use of RCA to visualize target DNA sequences as small as 50 nts in peripheral blood lymphocytes or in stretched DNA fibers. Three unique target sequences within the cystic fibrosis transmembrane conductance regulator gene could be detected simultaneously in interphase nuclei, and could be ordered in a linear map in stretched DNA. Allele-discriminating oligonucleotide probes in conjunction with RCA also were used to discriminate wild-type and mutant alleles in the cystic fibrosis transmembrane conductance regulator, p53, BRCA-1, and Gorlin syndrome genes in the nuclei of cultured cells or in DNA fibers. These observations demonstrate that signal amplification by RCA can be coupled to nucleic acid hybridization and multicolor fluorescence imaging to detect single nucleotide changes in DNA within a cytological context or in single DNA molecules. This provides a means for direct physical haplotyping and the analysis of somatic mutations on a cell-by-cell basis.

Detection of mitochondrial single nucleotide polymorphisms using a primer elongation reaction on oligonucleotide microarrays

We have developed a novel allele-specific primer elongation protocol using a DNA polymerase on oligonucleotide chips. Oligonucleotide primers carrying polymorphic sites at their free 3́end were covalently bound to glass slides. The generation of single-stranded targets of genomic DNA containing single nuclotide polymorphisms (SNPs) to be typed was achieved by an asymmetric PCR reaction or exonuclease treatment of phosphothioate (PTO)-modified PCR products. In the presence of DNA polymerase and all four dNTPs, with Cy3-dUTP replacing dTTP, allele-specific extension of the immobilized primers took place along a stretch of target DNA sequence. The yield of elongated products was increased by repeated reaction cycles. We performed multiplexed assays with many small DNA targets, or used single targets of up to 4.4 kb mitochondrial DNA (mtDNA) sequence to detect multiple SNPs in one reaction. The latter approach greatly simplifies preamplification of SNP-containing regions, thereby providing a framework for typing hundreds of mtDNA polymorphisms.

The evolution of coloration and toxicity in the poison frog family (Dendrobatidae)

The poison frogs (family Dendrobatidae) are terrestrial anuran amphibians displaying a wide range of coloration and toxicity. These frogs generally have been considered to be aposematic, but relatively little research has been carried out to test the predictions of this hypothesis. Here we use a comparative approach to test one prediction of the hypothesis of aposematism: that coloration will evolve in tandem with toxicity. Recently, we developed a phylogenetic hypothesis of the evolutionary relationships among representative species of poison frogs, using sequences from three regions of mitochondrial DNA. In our analysis, we use that DNA-based phylogeny and comparative analysis of independent contrasts to investigate the correlation between coloration and toxicity in the poison frog family (Dendrobatidae). Information on the toxicity of different species was obtained from the literature. Two different measures of the brightness and extent of coloration were used. (i) Twenty-four human observers were asked to rank different photos of each different species in the analysis in terms of contrast to a leaf-littered background. (ii) Color photos of each species were scanned into a computer and a computer program was used to obtain a measure of the contrast of the colors of each species relative to a leaf-littered background. Comparative analyses of the results were carried out with two different models of character evolution: gradual change, with branch lengths proportional to the amount of genetic change, and punctuational change, with all change being associated with speciation events. Comparative analysis using either method or model indicated a significant correlation between the evolution of toxicity and coloration across this family. These results are consistent with the hypothesis that coloration in this group is aposematic.

Exploring the DNA-binding specificities of zinc fingers with DNA microarrays

A key step in the regulation of networks that control gene expression is the sequence-specific binding of transcription factors to their DNA recognition sites. A more complete understanding of these DNA–protein interactions will permit a more comprehensive and quantitative mapping of the regulatory pathways within cells, as well as a deeper understanding of the potential functions of individual genes regulated by newly identified DNA-binding sites. Here we describe a DNA microarray-based method to characterize sequence-specific DNA recognition by zinc-finger proteins. A phage display library, prepared by randomizing critical amino acid residues in the second of three fingers of the mouse Zif268 domain, provided a rich source of zinc-finger proteins with variant DNA-binding specificities. Microarrays containing all possible 3-bp binding sites for the variable zinc fingers permitted the quantitation of the binding site preferences of the entire library, pools of zinc fingers corresponding to different rounds of selection from this library, as well as individual Zif268 variants that were isolated from the library by using specific DNA sequences. The results demonstrate the feasibility of using DNA microarrays for genome-wide identification of putative transcription factor-binding sites.

Multiple DNA and protein sequence alignment based on segment-to-segment comparison.

In this paper, a new way to think about, and to construct, pairwise as well as multiple alignments of DNA and protein sequences is proposed. Rather than forcing alignments to either align single residues or to introduce gaps by defining an alignment as a path running right from the source up to the sink in the associated dot-matrix diagram, we propose to consider alignments as consistent equivalence relations defined on the set of all positions occurring in all sequences under consideration. We also propose constructing alignments from whole segments exhibiting highly significant overall similarity rather than by aligning individual residues. Consequently, we present an alignment algorithm that (i) is based on segment-to-segment comparison instead of the commonly used residue-to-residue comparison and which (ii) avoids the well-known difficulties concerning the choice of appropriate gap penalties: gaps are not treated explicity, but remain as those parts of the sequences that do not belong to any of the aligned segments. Finally, we discuss the application of our algorithm to two test examples and compare it with commonly used alignment methods. As a first example, we aligned a set of 11 DNA sequences coding for functional helix-loop-helix proteins. Though the sequences show only low overall similarity, our program correctly aligned all of the 11 functional sites, which was a unique result among the methods tested. As a by-product, the reading frames of the sequences were identified. Next, we aligned a set of ribonuclease H proteins and compared our results with alignments produced by other programs as reported by McClure et al. [McClure, M. A., Vasi, T. K. & Fitch, W. M. (1994) Mol. Biol. Evol. 11, 571-592]. Our program was one of the best scoring programs. However, in contrast to other methods, our protein alignments are independent of user-defined parameters.

A cytoplasmic male sterility-associated mitochondrial protein causes pollen disruption in transgenic tobacco.

In higher plants, dominant mitochondrial mutations are associated with pollen sterility. This phenomenon is known as cytoplasmic male sterility (CMS). It is thought that the disruption in pollen development is a consequence of mitochondrial dysfunction. To provide definitive evidence that expression of an abnormal mitochondrial gene can interrupt pollen development, a CMS-associated mitochondrial DNA sequence from common bean, orf239, was introduced into the tobacco nuclear genome. Several transformants containing the orf239 gene constructs, with or without a mitochondrial targeting sequence, exhibited a semi sterile or male-sterile phenotype. Expression of the gene fusions in transformed anthers was confirmed using RNA gel blotting, ELISA, and light and electron microscopic immunocytochemistry. Immunocytological analysis showed that the ORF239 protein could associate with the cell wall of aberrant developing microspores. This pattern of extracellular localization was earlier observed in the CMS common bean line containing orf239 in the mitochondrial genome. Results presented here demonstrate that ORF239 causes pollen disruption in transgenic tobacco plants and may do so without targeting of the protein to the mitochondrion.

Efficient in vivo manipulation of mouse genomic sequences at the zygote stage.

We describe a transgenic mouse line carrying the cre transgene under the control of the adenovirus EIIa promoter that targets expression of the Cre recombinase to the early mouse embryo. To assess the ability of this recombinase to excise loxP-flanked DNA sequences at early stages of development, we bred EIIa-cre transgenic mice to two different mouse lines carrying loxP-flanked target sequences: (i) a strain with a single gene-targeted neomycin resistance gene flanked by 1oxP sites and (ii) a transgenic line carrying multiple transgene copies with internal loxP sites. Mating either of these loxP-carrying mouse lines to EIIa-cre mice resulted in first generation progeny in which the loxP-flanked sequences had been efficiently deleted from all tissues tested, including the germ cells. Interbreeding of these first generation progeny resulted in efficient germ-line transmission of the deletion to subsequent generations. These results demonstrate a method by which loxP-flanked DNA sequences can be efficiently deleted in the early mouse embryo. Potential applications of this approach are discussed, including reduction of multicopy transgene loci to produce single-copy transgenic lines and introduction of a variety of subtle mutations into the line.