DNA methylation and expression of estrogen receptor alpha in fathead minnows exposed to 17α-ethynylestradiol.

The gene expression response thought to underlie the negative apical effects resulting from estrogen exposure have been thoroughly described in fish. Although epigenetics are believed to play a critical role translating environmental exposures into the development of adverse apical effects, they remain poorly characterized in fish species. This study investigated alterations of DNA methylation of estrogen receptor alpha (esr1) in brain and liver tissues from 8 to 10 month old male fathead minnows (Pimephales promelas) after a 2d exposure to either 2.5 ng/L or 10 ng/L 17α-ethynylestradiol (EE2). Changes in the patterns of methylation were evaluated using targeted deep sequencing of bisulfite treated DNA in the 5' region of esr1. Methylation and gene expression were assessed at 2d of exposure and after a 7 and 14d depuration period. After 2d EE2 exposure, males exhibited significant demethylation in the 5' upstream region of esr1 in liver tissue, which was inversely correlated to gene expression. This methylation pattern reflected what was seen in females. No gene body methylation (GBM) was observed for liver of exposed males. Differential methylation was observed for a single upstream CpG site in the liver after the 14d depuration. A less pronounced methylation response was observed in the upstream region in brain tissue, however, several CpGs were necessarily excluded from the analysis. In contrast to the liver, a significant GBM response was observed across the entire gene body, which was sustained until at least 7d post-exposure. No differential expression was observed in the brain, limiting functional interpretation of methylation changes. The identification of EE2-dependent changes in methylation levels strongly suggests the importance of epigenetic mechanisms as a mediator of the organismal response to environmental exposures and the need for further characterization of the epigenome. Further, differential methylation following depuration indicates estrogenic effects persist well after the active exposure, which has implications for the risk posed by repeated exposures..

Phylogeny of rodentia (Mammalia) inferred from the nuclear-encoded gene IRBP.

The order Rodentia includes nearly half of all living mammalian species. Phylogenetic relationships among 22 species of rodents were investigated by use of a 1.2-kb region from exon 1 of the single-copy nuclear gene IRBP. IRBP has been extensively used for study of interordinal phylogeny in mammals, which allowed inclusion of 50 outgroup species, representing every eutherian order plus seven marsupials. Several clades were strongly supported, regardless of analytical method or inclusion/exclusion of data. These include a monophyletic Muroidea, with a clade including Spalax and Rhizomys as the first divergence; a clade uniting Zapus with Dipus, but excluding Sicista; a monophyletic Myodonta (Muroidea plus Dipodidae); and a clade including Aplodontidae as sister to Sciuridae. One bipartition, separating Hystricognathi and Geomyoidea from the remaining rodents, is strongly supported in all analyses that include third-position sites but almost completely absent from analyses that exclude third-position sites. A combination of nonstationary nucleotide composition and branch length effects may be causing all methods examined (including those using the LogDet distance) to support an incorrect conclusion when third-position sites are analyzed together with first- and second-position sites.

Parallel adaptive radiations in two major clades of placental mammals.

Higher level relationships among placental mammals, as well as the historical biogeography and morphological diversification of this group, remain unclear. Here we analyse independent molecular data sets, having aligned lengths of DNA of 5,708 and 2,947 base pairs, respectively, for all orders of placental mammals. Phylogenetic analyses resolve placental orders into four groups: Xenarthra, Afrotheria, Laurasiatheria, and Euarchonta plus Glires. The first three groups are consistently monophyletic with different methods of analysis. Euarchonta plus Glires is monophyletic or paraphyletic depending on the phylogenetic method. A unique nine-base-pair deletion in exon 11 of the BRCA1 gene provides additional support for the monophyly of Afrotheria, which includes proboscideans, sirenians, hyracoids, tubulidentates, macroscelideans, chrysochlorids and tenrecids. Laurasiatheria contains cetartiodactyls, perissodactyls, carnivores, pangolins, bats and eulipotyphlan insectivores. Parallel adaptive radiations have occurred within Laurasiatheria and Afrotheria. In each group, there are aquatic, ungulate and insectivore-like forms.

Mitochondrial versus nuclear gene sequences in deep-level mammalian phylogeny reconstruction.

Both mitochondrial and nuclear gene sequences have been employed in efforts to reconstruct deep-level phylogenetic relationships. A fundamental question in molecular systematics concerns the efficacy of different types of sequences in recovering clades at different taxonomic levels. We compared the performance of four mitochondrial data sets (cytochrome b, cytochrome oxidase II, NADH dehydrogenase subunit I, 12S rRNA-tRNA-16S rRNA) and eight nuclear data sets (exonic regions of alpha-2B adrenergic receptor, aquaporin, ss-casein, gamma-fibrinogen, interphotoreceptor retinoid binding protein, kappa-casein, protamine, von Willebrand Factor) in recovering deep-level mammalian clades. We employed parsimony and minimum-evolution with a variety of distance corrections for superimposed substitutions. In 32 different pairwise comparisons between these mitochondrial and nuclear data sets, we used the maximum set of overlapping taxa. In each case, the variable-length bootstrap was used to resample at the size of the smaller data set. The nuclear exons consistently performed better than mitochondrial protein and rRNA-tRNA coding genes on a per-residue basis in recovering benchmark clades. We also concatenated nuclear genes for overlapping taxa and made comparisons with concatenated mitochondrial protein-coding genes from complete mitochondrial genomes. The variable-length bootstrap was used to score the recovery of benchmark clades as a function of the number of resampled base pairs. In every case, the nuclear concatenations were more efficient than the mitochondrial concatenations in recovering benchmark clades. Among genes included in our study, the nuclear genes were much less affected by superimposed substitutions. Nuclear genes having appropriate rates of substitution should receive strong consideration in efforts to reconstruct deep-level phylogenetic relationships.

Maximum likelihood analysis of gene-based and structure-based process partitions, using mammalian mitochondrial genomes.

Aligned protein-coding genes from 19 completely sequenced mammalian mitochondrial genomes were examined by parsimony and maximum likelihood analyses. Particular attention is given to a comparison between gene-based and structure-based data partitions. Because actual structures are not known for most of the mitochondrially encoded proteins, three different surrogate partitioning schemes were examined, each based on the identity of the consensus amino acid at a specific homologous position. One of the amino-acid-based partitioning schemes gave the highest likelihood, but that scheme was based on concordance with a well-corroborated phylogeny from an earlier parsimony analysis. The gene-based partitioning scheme gave a significantly higher likelihood compared to the only structure-based scheme examined that could be generated without prior assumptions about the phylogeny. Two contrasting phylogenetic inferences were supported by the analyses. Both unpartitioned analyses and analyses in which all partitions were constrained to have identical patterns of branch lengths supported ((Artiodactyla, Cetacea) (Perissodactyla, Carnivora)), whereas all analyses with that constraint relaxed supported (((Artiodactyla, Cetacea) Carnivora) Perissodactyla).

Comparative analysis of evolution in a rodent histone H2a pseudogene.

Sequences were obtained from five species of rodents that are orthologous to an H2a histone pseudogene from Mus musculus. The pseudogene is part of the cluster of replication-dependent histone genes found on Mus musculus chromosome 13. Comparative analysis of these five sequences together with the previously published sequence from M. musculus shows that this gene has likely been a pseudogene throughout the evolution of the genus Mus, while the gene from Rattus norvegicus is likely functional. Three large (> 20 bp) deletions were found among the Mus pseudogenes, a feature that is very unusual compared to surveys of processed pseudogenes. In addition, there are two single-base deletions and one 4-bp insertion among the Mus pseudogenes. The species distributions of one of the large deletions and the 4-bp insertion require either independent insertions of an identical sequence, independent deletions with identical boundaries, or a deletion followed by precise reintegration of the original sequence. The evidence favors the hypothesis of multiple deletions with identical boundaries. The "coding" regions of the Mus pseudogenes show a much reduced level of among-species variability in the 3' half of the pseudogene, compared both to the 5' half and to flanking sequences. This supports a hypothesis that the 3' end of the pseudogene is the target of frequent gene conversion by functional H2a genes.

Characterization of the 55-kb mouse histone gene cluster on chromosome 3.

The histone gene cluster on mouse chromosome 3 has been isolated as a series of overlapping P1 clones, covering 110-120 kb, by probing with the histone H3-614 gene that had been mapped previously to mouse chromosome 3. There are genes for 10 core histone proteins present in a 55-kb cluster within this contig. There are three histone H3 genes, two of which are identical; four histone H2a genes, two of which are identical, one histone H4 gene; and two histone H2b genes. These histone H3 and H2a genes encode approximately 40% of the total H3 and H2a mRNA, whereas the histone H4 and histone H2b genes encode < 10% of the total H4 and H2b mRNA. There are no histone H1 genes present in this cluster. All of the histone H2a genes encode histone H2a.2 proteins (or variants of H2a.2), and account for all the H2a.2 genes in the mouse genome. All three histone H3 genes encode the histone H3.2 protein. A 21-kb region containing the adjacent H3-614 and H2a-614 genes has been duplicated and is present in an inverted repeat separated by 4.5 kb. The other two H2a genes are adjacent, with the 3' ends of their mRNAs separated by only 49 nucleotides in the DNA and the U7 snRNP binding sites separated by only 20 nucleotides. One of the histone H2b genes has lost the stem-loop sequence characteristic of the replication-dependent histone mRNAs and encodes only polyadenylated mRNAs.

Selection on silent sites in the rodent H3 histone gene family.

Selection promoting differential use of synonymous codons has been shown for several unicellular organisms and for Drosophila, but not for mammals. Selection coefficients operating on synonymous codons are likely to be extremely small, so that a very large effective population size is required for selection to overcome the effects of drift. In mammals, codon-usage bias is believed to be determined exclusively by mutation pressure, with differences between genes due to large-scale variation in base composition around the genome. The replication-dependent histone genes are expressed at extremely high levels during periods of DNA synthesis, and thus are among the most likely mammalian genes to be affected by selection on synonymous codon usage. We suggest that the extremely biased pattern of codon usage in the H3 genes is determined in part by selection. Silent site G + C content is much higher than expected based on flanking sequence G + C content, compared to other rodent genes with similar silent site base composition but lower levels of expression. Dinucleotide-mediated mutation bias does affect codon usage, but the affect is limited to the choice between G and C in some fourfold degenerate codons. Gene conversion between the two clusters of histone genes has not been an important force in the evolution of the H3 genes, but gene conversion appears to have had some effect within the cluster on chromosome 13.

The consistency of several phylogeny-inference methods under varying evolutionary rates.

A phylogenetic method is a consistent estimator of phylogeny if and only if it is guaranteed to give the correct tree, given that sufficient (possibly infinite) independent data are examined. The following methods are examined for consistency: UPGMA (unweighted pair-group method, averages), NJ (neighbor joining), MF (modified Farris), and P (parsimony). A two-parameter model of nucleotide sequence substitution is used, and the expected distribution of character states is calculated. Without perfect correction for superimposed substitutions, all four methods may be inconsistent if there is but one branch evolving at a faster rate than the other branches. Partial correction of observed distances improves the robustness of the NJ method to rate variation, and perfect correction makes the NJ method a consistent estimator for all combinations of rates that were examined. The sensitivity of all the methods to unequal rates varies over a wide range, so relative-rate tests are unlikely to be a reliable guide for accepting or rejecting phylogenies based on parsimony analysis.