Evolution of gilled mushrooms and puffballs inferred from ribosomal DNA sequences

Homobasidiomycete fungi display many complex fruiting body morphologies, including mushrooms and puffballs, but their anatomical simplicity has confounded efforts to understand the evolution of these forms. We performed a comprehensive phylogenetic analysis of homobasidiomycetes, using sequences from nuclear and mitochondrial ribosomal DNA, with an emphasis on understanding evolutionary relationships of gilled mushrooms and puffballs. Parsimony-based optimization of character states on our phylogenetic trees suggested that strikingly similar gilled mushrooms evolved at least six times, from morphologically diverse precursors. Approximately 87% of gilled mushrooms are in a single lineage, which we call the “euagarics.” Recently discovered 90 million-year-old fossil mushrooms are probably euagarics, suggesting that (i) the origin of this clade must have occurred no later than the mid-Cretaceous and (ii) the gilled mushroom morphology has been maintained in certain lineages for tens of millions of years. Puffballs and other forms with enclosed spore-bearing structures (Gasteromycetes) evolved at least four times. Derivation of Gasteromycetes from forms with exposed spore-bearing structures (Hymenomycetes) is correlated with repeated loss of forcible spore discharge (ballistospory). Diverse fruiting body forms and spore dispersal mechanisms have evolved among Gasteromycetes. Nevertheless, it appears that Hymenomycetes have never been secondarily derived from Gasteromycetes, which suggests that the loss of ballistospory has constrained evolution in these lineages.

Instability of repetitive DNA sequences: The role of replication in multiple mechanisms

Rearrangements between tandem sequence homologies of various lengths are a major source of genomic change and can be deleterious to the organism. These rearrangements can result in either deletion or duplication of genetic material flanked by direct sequence repeats. Molecular genetic analysis of repetitive sequence instability in Escherichia coli has provided several clues to the underlying mechanisms of these rearrangements. We present evidence for three mechanisms of RecA-independent sequence rearrangements: simple replication slippage, sister-chromosome exchange-associated slippage, and single-strand annealing. We discuss the constraints of these mechanisms and contrast their properties with RecA-dependent homologous recombination. Replication plays a critical role in the two slipped misalignment mechanisms, and difficulties in replication appear to trigger rearrangements via all these mechanisms.

Differences in water release for the binding of EcoRI to specific and nonspecific DNA sequences.

The free energy difference between complexes of the restriction nuclease EcoRI with nonspecific DNA and with the enzyme's recognition sequence is linearly dependent on the water chemical potential of the solution, set using several very different solutes, ranging from glycine and glycerol to triethylene glycol and sucrose. This osmotic dependence indicates that the nonspecific complex sequesters some 110 waters more than the specific complex with the recognition sequence. The insensitivity of the difference in number of waters released to the solute identity further indicates that this water is sequestered in a space that is sterically inaccessible to solutes, most likely at the protein-DNA interface of the nonspecific complex. Calculations based on the structure of the specific complex suggest that the apposing DNA and protein surfaces in the nonspecific complex retain approximately a full hydration layer of water.

Does behavior reflect phylogeny in swiftlets (Aves: Apodidae)? A test using cytochrome b mitochondrial DNA sequences.

Swiftlets are small insectivorous birds, many of which nest in caves and are known to echolocate. Due to a lack of distinguishing morphological characters, the taxonomy of swiftlets is primarily based on the presence or absence of echolocating ability, together with nest characters. To test the reliability of these behavioral characters, we constructed an independent phylogeny using cytochrome b mitochondrial DNA sequences from swiftlets and their relatives. This phylogeny is broadly consistent with the higher classification of swifts but does not support the monophyly of swiftlets. Echolocating swiftlets (Aerodramus) and the nonecholocating "giant swiftlet" (Hydrochous gigas) group together, but the remaining nonecholocating swiftlets belonging to Collocalia are not sister taxa to these swiftlets. While echolocation may be a synapomorphy of Aerodramus (perhaps secondarily lost in Hydrochous), no character of Aerodramus nests showed a statistically significant fit to the molecular phylogeny, indicating that nest characters are not phylogenetically reliable in this group.

DNA sequences of Alu elements indicate a recent replacement of the human autosomal genetic complement.

DNA sequences of neutral nuclear autosomal loci, compared across diverse human populations, provide a previously untapped perspective into the mode and tempo of the emergence of modern humans and a critical comparison with published clonally inherited mitochondrial DNA and Y chromosome measurements of human diversity. We obtained over 55 kilobases of sequence from three autosomal loci encompassing Alu repeats for representatives of diverse human populations as well as orthologous sequences for other hominoid species at one of these loci. Nucleotide diversity was exceedingly low. Most individuals and populations were identical. Only a single nucleotide difference distinguished presumed ancestral alleles from descendants. These results differ from those expected if alleles from divergent archaic populations were maintained through multiregional continuity. The observed virtual lack of sequence polymorphism is the signature of a recent single origin for modern humans, with general replacement of archaic populations.

Inferring phylogenies from DNA sequences of unequal base compositions.

A new method for computing evolutionary distances between DNA sequences is proposed. Contrasting with classical methods, the underlying model does not assume that sequence base compositions (A, C, G, and T contents) are at equilibrium, thus allowing unequal base compositions among compared sequences. This makes the method more efficient than the usual ones in recovering phylogenetic trees from sequence data when base composition is heterogeneous within the data set, as we show by using both simulated and empirical data. When applied to small-subunit ribosomal RNA sequences from several prokaryotic or eukaryotic organisms, this method provides evidence for an early divergence of the microsporidian Vairimorpha necatrix in the eukaryotic lineage.

Coiled bodies contain U7 small nuclear RNA and associate with specific DNA sequences in interphase human cells.

Coiled bodies (CBs) are nuclear organelles whose structures appear to be highly conserved in evolution. In rapidly cycling cells, they are typically located in the nucleoplasm but are often found in contact with the nucleolus. The CBs in human cells contain a unique protein, called p80-coilin. Studies on amphibian oocyte nuclei have revealed a protein within the "sphere" organelle that shares significant structural similarity to p80-coilin. Spheres and CBs are also highly enriched in small nuclear ribonucleoproteins and other RNA-processing components. We present evidence that, like spheres, CBs contain U7 small nuclear RNA (snRNA) and associate with specific chromosomal loci. Using biotinylated 2'-O-methyl oligonucleotides complementary to the 5' end of U7 snRNA and fluorescence in situ hybridization, we show that U7 is distributed throughout the nucleoplasm, excluding nucleoli, and is concentrated in CBs. Interestingly, we found that CBs often associate with subsets of the histone, U1, and U2 snRNA gene loci in interphase HeLa-ATCC and HEp-2 monolayer cells. However, in a strain of suspension-grown HeLa cells, called HeLa-JS1000, we found a much lower rate of association between CBs and snRNA genes. Possible roles for CBs in the metabolism of these various histone and snRNAs are discussed.

Specific molecular recognition of mixed nucleic acid sequences: An aromatic dication that binds in the DNA minor groove as a dimer

Phenylamidine cationic groups linked by a furan ring (furamidine) and related compounds bind as monomers to AT sequences of DNA. An unsymmetric derivative (DB293) with one of the phenyl rings of furamidine replaced with a benzimidazole has been found by quantitative footprinting analyses to bind to GC-containing sites on DNA more strongly than to pure AT sequences. NMR structural analysis and surface plasmon resonance binding results clearly demonstrate that DB293 binds in the minor groove at specific GC-containing sequences of DNA in a highly cooperative manner as a stacked dimer. Neither the symmetric bisphenyl nor bisbenzimidazole analogs of DB293 bind significantly to the GC containing sequences. DB293 provides a paradigm for design of compounds for specific recognition of mixed DNA sequences and extends the boundaries for small molecule-DNA recognition.

[Leu5]enkephalin-encoding sequences are targets for a specific DNA-binding factor.

A DNA-binding factor with high affinity and specificity for the [Leu5]enkephalin-encoding sequences in the prodynorphin and proenkephalin genes has been characterized. The factor has the highest affinity for the [Leu5]-enkephalin-encoding sequence in the dynorphin B-encoding region of the prodynorphin gene, has relatively high affinity for other [Leu5]enkephalin-encoding sequences in the prodynorphin and proenkephalin genes, but has no apparent affinity for similar DNA sequences coding for [Met5]-enkephalin in the prodynorphin or proopiomelanocortin genes. The factor has been named [Leu5]enkephalin-encoding sequence DNA-binding factor (LEF). LEF has a nuclear localization and is composed of three subunits of about 60, 70, and 95 kDa, respectively. The highest levels were observed in rat testis, cerebellum, and spleen and were generally higher in late embryonal compared to newborn or adult animals. LEF activity was also recorded in human clonal tumor cell lines. LEF inhibited the transcription of reporter genes in artificial gene constructs where a [Leu5]enkephalin-encoding DNA fragment had been inserted between the transcription initiation site and the coding region of the reporter genes. These observations suggest that the [Leu5]enkephalin-encoding sequences in the prodynorphin and proenkephalin genes also have regulatory functions realized through interaction with a specific DNA-binding factor.

Two sterol regulatory element-like sequences mediate up-regulation of caveolin gene transcription in response to low density lipoprotein free cholesterol

Caveolae form the terminus for a major pathway of intracellular free cholesterol (FC) transport. Caveolin mRNA levels in confluent human skin fibroblasts were up-regulated following increased uptake of low density lipoprotein (LDL) FC. The increase induced by FC was not associated with detectable change in mRNA stability, indicating that caveolin mRNA levels were mediated at the level of gene transcription. A total of 924 bp of 5′ flanking region of the caveolin gene were cloned and sequenced. The promoter sequence included three G+C-rich potential sterol regulatory elements (SREs), a CAAT sequence and a Sp1 consensus sequence. Deletional mutagenesis of individual SRE-like sequences indicated that of these two (at −646 and −395 bp) were essential for the increased transcription rates mediated by LDL-FC, whereas the third was inconsequential. Gel shift analysis of protein binding from nuclear extracts to these caveolin promoter DNA sequences, together with DNase I footprinting, confirmed nucleoprotein binding to the SRE-like elements as part of the transcriptional response to LDL-FC. A supershift obtained with antibody to SRE-binding protein 1 (SPEBP-1) indicated that this protein binds at −395 bp. There was no reaction at −395 bp with anti-Sp1 antibody nor with either antibody at −646 bp. The cysteine protease inhibitor N-acetyl-leu-leu-norleucinal (ALLN), which inhibits SREBP catabolism, superinhibited caveolin mRNA levels regardless of LDL-FC. This finding suggests that SREBP inhibits caveolin gene transcription in contrast to its stimulating effect on other promoters. The findings of this study are consistent with the postulated role for caveolin as a regulator of cellular FC homeostasis in quiescent peripheral cells, and the coordinate regulation by SREBP of FC influx and efflux.

Inhibition of RNA polymerase II transcription in human cells by synthetic DNA-binding ligands

Sequence-specific DNA-binding small molecules that can permeate human cells potentially could regulate transcription of specific genes. Multiple cellular DNA-binding transcription factors are required by HIV type 1 for RNA synthesis. Two pyrrole–imidazole polyamides were designed to bind DNA sequences immediately adjacent to binding sites for the transcription factors Ets-1, lymphoid-enhancer binding factor 1, and TATA-box binding protein. These synthetic ligands specifically inhibit DNA-binding of each transcription factor and HIV type 1 transcription in cell-free assays. When used in combination, the polyamides inhibit virus replication by >99% in isolated human peripheral blood lymphocytes, with no detectable cell toxicity. The ability of small molecules to target predetermined DNA sequences located within RNA polymerase II promoters suggests a general approach for regulation of gene expression, as well as a mechanism for the inhibition of viral replication.

A zinc finger directory for high-affinity DNA recognition

We have used two monovalent phage display libraries containing variants of the Zif268 DNA-binding domain to obtain families of zinc fingers that bind to alterations in the last 4 bp of the DNA sequence of the Zif268 consensus operator, GCG TGGGCG. Affinity selection was performed by altering the Zif268 operator three base pairs at a time, and simultaneously selecting for sets of 16 related DNA sequences. In this way, only four experiments were required to select for all possible 64 combinations of DNA triplet sequences. The results show that (i) for high-affinity DNA binding in the range observed for the Zif268 wild-type complex (Kd = 0.5–5 nM), finger 1 specifically requires the arginine at the carboxy terminus of its recognition helix that forms a bidentate hydrogen-bond with the guanine base (G) in the crystal structure of Zif268 complexed to its DNA operator sequence GCG TGG GCG; (ii) when the guanine base (G) is replaced by A, C, or T, a lower-affinity family (Kd ⩾ 50 nM) can be detected that shows an overall tendency to bind G-rich DNA; (iii) the residues at position 2 on the finger 2 recognition helix do not appear to interact strongly with the complementary 5′ base in the finger 1 binding site; and (iv) unexpected substitutions at the amino terminus of finger 1 can occasionally result in specificity for the 3′ base in the finger 1 binding site. A DNA recognition directory was constructed for high-affinity zinc fingers that recognize all three bases in a DNA triplet for seven sequences of the type GNN. Similar approaches may be applied to other zinc fingers to broaden the scope of the directory.

Conserved characteristics of heterochromatin-forming DNA at the 15q11-q13 imprinting center

Nuclear matrix binding assays (NMBAs) define certain DNA sequences as matrix attachment regions (MARs), which often have cis-acting epigenetic regulatory functions. We used NMBAs to analyze the functionally important 15q11-q13 imprinting center (IC). We find that the IC is composed of an unusually high density of MARs, located in close proximity to the germ line elements that are proposed to direct imprint switching in this region. Moreover, we find that the organization of MARs is the same at the homologous mouse locus, despite extensive divergence of DNA sequence. MARs of this size are not usually associated with genes but rather with heterochromatin-forming areas of the genome. In contrast, the 15q11-q13 region contains multiple transcribed genes and is unusual for being subject to genomic imprinting, causing the maternal chromosome to be more transcriptionally silent, methylated, and late replicating than the paternal chromosome. We suggest that the extensive MAR sequences at the IC are organized as heterochromatin during oogenesis, an organization disrupted during spermatogenesis. Consistent with this model, multicolor fluorescence in situ hybridization to halo nuclei demonstrates a strong matrix association of the maternal IC, whereas the paternal IC is more decondensed, extending into the nuclear halo. This model also provides a mechanism for spreading of the imprinting signal, because heterochromatin at the IC on the maternal chromosome may exert a suppressive position effect in cis. We propose that the germ line elements at the 15q11-q13 IC mediate their effects through the candidate heterochromatin-forming DNA identified in this study.

Analysis of microsatellite mutations in the mitochondrial DNA of Saccharomyces cerevisiae

In the nuclear genome of Saccharomyces cerevisiae, simple, repetitive DNA sequences (microsatellites) mutate at rates much higher than nonrepetitive sequences. Most of these mutations are deletions or additions of repeat units. The yeast mitochondrial genome also contains many microsatellites. To examine the stability of these sequences, we constructed a reporter gene (arg8m) containing out-of-frame insertions of either poly(AT) or poly(GT) tracts within the coding sequence. Yeast strains with this reporter gene inserted within the mitochondrial genome were constructed. Using these strains, we showed that poly(GT) tracts were considerably less stable than poly(AT) tracts and that alterations usually involved deletions rather than additions of repeat units. In contrast, in the nuclear genome, poly(GT) and poly(AT) tracts had similar stabilities, and alterations usually involved additions rather than deletions. Poly(GT) tracts were more stable in the mitochondria of diploid cells than in haploids. In addition, an msh1 mutation destabilized poly(GT) tracts in the mitochondrial genome.

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.