Conserved vertebrate chromosome segments in the large salamander genome

Urodele amphibians (salamanders) are important models for embryological, physiological, and natural history research and are also a biomedically important group because they are the only vertebrates capable of regenerating entire organ systems. To enhance the utility of salamanders for biomedical research and for understanding genome evolution, genetic linkage analysis was used to identify chromosome segments that are homologous between ambystomatid salamanders and distantly related vertebrate model organisms. A total of 347 loci (AFLPs, RAPDs, and protein-coding loci) were mapped using an interspecific meiotic mapping panel (Ambystoma mexicanum and A. tigrinum tigrinum; family Ambystomatidae). Genome size in Ambystoma was estimated to be 7291 cM, the largest linkage map estimate reported for any organism. However, the relatively large size of the salamander genome did not hinder efforts to map and identify conserved syntenies from a small sample of 24 protein-coding loci. Chromosomal segments that are conserved between fishes and mammals are also conserved in these salamanders. Thus, comparative gene mapping appears to be an efficient strategy for identifying orthologous loci between ambystomatid salamanders and genomically well-characterized vertebrate model organisms.     

Rates of Nucleotide Substitution in Angiosperm Mitochondrial DNA Sequences and Dates of Divergence Between Brassica and Other Angiosperm Lineages

We obtained 16 nucleotide sequences (∼1400 bp each) of the first intron of the mitochondrial (mt) gene for NADH subunit 4 (nad4) from 10 species of Brassicaceae. Using these new sequences and five published sequences from GenBank, we constructed a phylogenetic tree of the Brassicaceae species under study and showed that the rate of nucleotide substitution in the first intron of nad4 is very low, about 0.16–0.23 × 10⁻⁹ substitution per site per year, which is about half of the silent rate in exons of nad4. The ratios of substitution rates in this intron, ITS, and IGS are approximately 1:23:73, where ITS is the nuclear intergenic spacer between 18S and 25S rRNA genes and IGS is the intergenic spacer of 5S rRNA genes. A segment (335 bp) in the first intron of nad4 in Brassicaceae species that is absent in wheat was considered as a nonfunctional sequence and used to estimate the neutral rate (the rate of mutation) in mtDNA to be 0.5–0.7 × 10⁻⁹ substitution per site per year, which is about three times higher than the substitution rate in the rest of the first intron of nad4. We estimated that the dates of divergence are 170–235 million years (Myr) for the monocot–dicot split, 112–156 Myr for the Brassicaceae–Lettuce split, 14.5–20.4 Myr for the Brassica–Arabidopsis split, and 14.5–20.4 Myr for the Arabidopsis–Arabideae split. Received: 14 July 1998 / Accepted: 1 October 1998     

Unisexual salamanders (genus Ambystoma) present a new reproductive mode for eukaryotes.

To persist, unisexual and asexual eukaryotes must have reproductive modes that circumvent normal bisexual reproduction. Parthenogenesis, gynogenesis, and hybridogenesis are the modes that have generally been ascribed to various unisexuals. Unisexual Ambystoma are abundant around the Great Lakes region of North America, and have variously been described as having all 3 reproductive modes. Diploid and polyploid unisexuals have nuclear genomes that combine the haploid genomes of 2 to 4 distinct sexual species, but the mtDNA is unlike any of those 4 species and is similar to another species, Ambystoma barbouri. To obtain better resolution of the reproductive mode used by unisexual Ambystoma and to explore the relationship of A. barbouri to the unisexuals, we sequenced the mitochondrial control and highly variable intergenic spacer region of 48 ambystomatids, which included 28 unisexuals, representatives of the 4 sexual species and A. barbouri. The unisexuals have similar sequences over most of their range, and form a close sister group to A. barbouri, with an estimated time of divergence of 2.4-3.9 million years ago. Individuals from the Lake Erie Islands (Kelleys, Pelee, North Bass) have a haplotype that demonstrates an isolation event. We examined highly variable microsatellite loci, and found that the genetic makeup of the unisexuals is highly variable and that unisexual individuals share microsatellite alleles with sexual individuals within populations. Although many progeny from the same female had the same genotype for 5 microsatellite DNA loci, there was no indication that any particular genome is consistently inherited in a clonal fashion in a population. The reproductive mode used by unisexual Ambystoma appears to be unique; we suggest kleptogenesis as a new unisexual reproductive mode that is used by these salamanders.     

A comprehensive expressed sequence tag linkage map for tiger salamander and Mexican axolotl: enabling gene mapping and comparative genomics in Ambystoma

Expressed sequence tag (EST) markers were developed for Ambystoma tigrinum tigrinum (Eastern tiger salamander) and for A. mexicanum (Mexican axolotl) to generate the first comprehensive linkage map for these model amphibians. We identified 14 large linkage groups (125.5-836.7 cM) that presumably correspond to the 14 haploid chromosomes in the Ambystoma genome. The extent of genome coverage for these linkage groups is apparently high because the total map size (5251 cM) falls within the range of theoretical estimates and is consistent with independent empirical estimates. Unlike most vertebrate species, linkage map size in Ambystoma is not strongly correlated with chromosome arm number. Presumably, the large physical genome size ( approximately 30 Gbp) is a major determinant of map size in Ambystoma. To demonstrate the utility of this resource, we mapped the position of two historically significant A. mexicanum mutants, white and melanoid, and also met, a quantitative trait locus (QTL) that contributes to variation in metamorphic timing. This new collection of EST-based PCR markers will better enable the Ambystoma system by facilitating development of new molecular probes, and the linkage map will allow comparative studies of this important vertebrate group.     

Biogeography of the southeastern United States: a comparison of salamander phylogeographic studies

Most phylogeographic studies of species from the southeastern United States have shown a simple east-west division of mtDNA variation. However, a study of the salamander Ambystoma maculatum resulted in a more complex pattern that includes a close affinity between populations from the Central Highlands of Missouri and Arkansas and the Coastal Plain separated by a genetically distinct central group of populations. We test the generality of this observation by surveying mitochondrial DNA (mtDNA) variation in the closely related species A. talpoideum. An Ambystoma-specific intergenic spacer was amplified and sequenced. The 26 resulting haplotypes varied from 380 to 800 base pairs, and alignments, including the outgroup, required 101 insertions/deletions. Sequence divergence among haplotypes ranged from 0.001 to 0.758. Population subdivision was extensive (theta = 0.64). Phylogenetic analysis of A. talpoideum mtDNA sequence reveals a close relationship between the populations from the Central Highlands and the Coastal Plain. This result is similar to that obtained for A. maculatum, although the A. talpoideum clade is not as well differentiated from its sister clades. We discuss the differences and similarities between the two Ambystoma species and previous studies and call for increased focus on multiple species with similar ecologies as a way to detect subtle biogeographic events.     

Specific dynamic action of ambystomatid salamanders and the effects of meal size, meal type, and body temperature

The past decade has witnessed a dramatic increase in studies of amphibian and reptile specific dynamic action (SDA). These studies have demonstrated that SDA, the summed energy expended on meal digestion and assimilation, is affected significantly by meal size, meal type, and body size and to some extent by body temperature. While much of this attention has been directed at anuran and reptile SDA, we investigated the effects of meal size, meal type, and body temperature on the postprandial metabolic responses and the SDA of the tiger salamander (Ambystoma tigrinum tigrinum). We also compared the SDA responses among six species of Ambystoma salamanders representing the breadth of Ambystoma phylogeny. Postprandial peaks in VO(2) and VO(2), duration of elevated metabolism, and SDA of tiger salamanders increased with the size of cricket meals (2.5%-12.5% of body mass). For A. tigrinum, as for other ectotherms, a doubling of meal size results in an approximate doubling of SDA, a function of equal increases in peak VO(2) and duration. For nine meal types of equivalent size (5% of body mass), the digestion of hard-bodied prey (crickets, superworms, mealworms, beetles) generated larger SDA responses than the digestion of soft-bodied prey (redworms, beetle larvae). Body temperature affected the profile of postprandial metabolism, increasing the peak and shortening the duration of the profile as body temperature increased. SDA was equivalent among three body temperatures (20 degrees, 25 degrees, and 30 degrees C) but decreased significantly at 15 degrees C. Comparatively, the postprandial metabolic responses and SDA of Ambystoma jeffersonianum, Ambystoma maculatum, Ambystoma opacum, Ambystoma talpoideum, Ambystoma texanum, and the conspecific Ambystoma tigrinum mavortium digesting cricket meals that were 5% of their body mass were similar (independent of body mass) to those of A. t. tigrinum. Among the six species, standard metabolic rate, peak postprandial VO(2), and SDA scaled with body mass with mass exponents of 0.72, 0.78, and 1.05, respectively.     

Multiple nuclear gene sequences identify phylogenetic species boundaries in the rapidly radiating clade of Mexican ambystomatid salamanders

Delimiting the boundaries of species involved in radiations is critical to understanding the tempo and mode of lineage formation. Single locus gene trees may or may not reflect the underlying pattern of population divergence and lineage formation, yet they constitute the vast majority of the empirical data in species radiations. In this study we make use of an expressed sequence tag (EST) database to perform nuclear (nDNA) and mitochondrial (mtDNA) genealogical tests of species boundaries in Ambystoma ordinarium, a member of an adaptive radiation of metamorphic and paedomorphic salamanders (the Ambystoma tigrinum complex) that have diversified across terrestrial and aquatic environments. Gene tree comparisons demonstrate extensive nonmonophyly in the mtDNA genealogy of A. ordinarium, while seven of eight independent nuclear loci resolve the species as monophyletic or nearly so, and diagnose it as a well-resolved genealogical species. A differential introgression hypothesis is supported by the observation that western A. ordinarium localities contain mtDNA haplotypes that are identical or minimally diverged from haplotypes sampled from a nearby paedomorphic species, Ambystoma dumerilii, while most nDNA trees place these species in distant phylogenetic positions. These results provide a strong example of how historical introgression can lead to radical differences between gene trees and species histories, even among currently allopatric species with divergent life history adaptations and morphologies. They also demonstrate how EST-based nuclear resources can be used to more fully resolve the phylogenetic history of species radiations.     

Sequence Characterization of a Unique Intergenic Spacer in Gadiformes Mitochondrial DNA

The nucleotide sequences of intergenic spacers located between the tRNA^Thr and tRNA^Pro genes in mitochondrial DNA of cod fishes (order Godiformes) were determined. Spacers from eight species representing two families of cod fishes were analyzed and found to vary in size from 25 to 99 bp. Each spacer sequence contains one or two copies of a conserved 17-bp motif. Four to five central nucleotides of this motif constitute a substitutional hot spot as observed from interspecific and intraspecific comparisons. The substitution rate of the spacer is approximately twice that of the variable part I of the mitochondrial DNA control region, making this sequence region interesting as a molecular marker in population studies or stock assessments of cod fishes. We propose that the spacer originated in a duplication event and evolved into a functional domain, perhaps by binding regulatory proteins. Accepted February 26, 1999     

丝带凤蝶线粒体基因组全序列及其系统学意义(英文)

锯凤蝶类与凤蝶科其他类群的系统发生关系及其分类学地位一直存在争议。本研究采用PCR和long PCR技术测定了属于锯凤蝶类的丝带凤蝶Sericinus montelus线粒体基因组全序列; 结合已有的其他凤蝶科物种的相应序列数据, 基于13个蛋白质编码基因重建了凤蝶科主要类群的系统发生树, 探讨了它们之间的系统发生关系。基因组分析结果表明： 丝带凤蝶线粒体基因组全长15 242 bp, 包括13个编码蛋白基因(ATP6, ATP8, COⅠ-Ⅲ, ND1-6, ND4L和Cytb)、 22个tRNA基因、 16S和12S rRNA基因以及非编码的控制区; 基因组A, T, G和C含量分别为40.1%, 40.8%, 7.4%和11.7%, 表现出明显的AT偏倚。所有的蛋白质编码基因都使用标准的起始密码子(ATN); 除ND4 和 ND4L基因使用单个的T作为终止密码子外, 其余蛋白编码基因都使用了标准的终止密码子(TAA)。除丝氨酸 tRNA的二氢尿苷突环缺失外, 所有tRNA基因都形成典型的三叶草型结构。基因组中共存在12个大小介于2~65 bp之间的基因间隔区以及15个大小介于1~8 bp之间的基因重叠区, 其中, 存在于COⅡ和tRNA^Lys之间的24 bp的间隔区在其他鳞翅目昆虫中未曾见到。以邻接法和最大简约法并基于13个蛋白质编码基因序列对凤蝶科进行了系统发生分析。结果显示, 丝带凤蝶和中华虎凤蝶Luehdorfia chinensis先构成一个支系, 再和冰清绢蝶Parnassius bremeri构成姊妹群; 表明锯凤蝶类应作为族级分类单元归于凤蝶科下的绢蝶亚科。     

中国淋巴囊肿病毒胸苷酸合酶基因结构特点及分析

胸苷酸合酶(Thymidylate synthase,TS)是进行DNA合成所必需的酶类,与细胞分化及肿瘤发生密切相关。淋巴囊肿病毒属于虹彩病毒科成员,是能引起百余种淡、海水鱼感染,并产生肿瘤的病毒病原。在已完成中国淋巴囊肿病毒株(Lymphocystis disease virus-China,LCDV-C)基因组序列测定的基础上,本文对位于LCDV-C基因组开放阅读框ORF011L的TS基因结构、及其推定蛋白结构进行了分析。该基因全长858bp,GC含量为28．2％,编码一个长为286aa、分子量为32．7kD、等电点为7．1的推定蛋白,称之为中国淋巴囊肿病毒胸苷酸合酶(LCDV-CTS)。该酶所具有的叶酸结合区在第44和70位氨基酸之间,dUMP结合区在第163和206位氨基酸之间,24个必需氨基酸具有高度保守性。二级结构预测结果显示LCDV-C TS含8个a螺旋,6个β折叠和23个环,表明其具备酶活性分子所有的柔性和可变性结构特征。这是迄今所知在脊椎动物虹彩病毒中唯一含两个完整结合区的TS。对来自包括LCDV在内二十个物种的TS结构进行同源性分析,显示LCDV-C TS被单独分为一枝。进一步对LCDV-C TS可能的起源途径及与宿主的作用等进行了探讨。     

A DNA-based assay identifies Batrachochytrium dendrobatidis in amphibians

Chytridiomycosis caused by Batrachochytrium dendrobatidis (Chytridiomycota) has been implicated in declines of amphibian populations on four continents. We have developed a sensitive and specific polymerase chain reaction-based assay to detect this pathogen. We isolated B. dendrobatidis from captive and wild amphibians collected across North America and sequenced the internal transcribed spacer regions of the rDNA cassette of multiple isolates. We identified two primers (Bd1a and Bd2a) that are specific to B. dendrobatidis under amplification conditions described in this study. DNA amplification with Bd1a/Bd2a primers produced a fragment of approximately 300 bp from B. dendrobatidis DNA but not from DNA of other species of chytrids or common soil fungi. The assay detected 10 zoospores or 10 pg of DNA from B. dendrobatidis and detected infections in skin samples from a tiger salamander (Ambystoma tigrinum), boreal toads (Bufo boreas), Wyoming toads (Bufo baxteri), and smooth-sided toads (Bufo guttatus). This assay required only small samples of skin and can be used to process a large number of samples.     

Evolution of the A+T-rich region of mitochondrial DNA in the melanogaster species subgroup of Drosophila

We determined the nucleotide sequences of two regions in the A+T-rich region of mitochondrial DNA (mtDNA) in the siI and siII types of D. simulans, the maII type of D. mauritiana, and D. sechellia. The sequences were aligned with those of the corresponding regions of siIII of D. simulans and maI of D. mauritiana, D. melanogaster, and D. yakuba. The type I and type II elements and the T-stretches were detected in all eight of the mtDNA types compared, indicating that the three elements are essential in the A+T-rich region of this species subgroup. The alignment revealed several short repetitive sequences and relatively large deletions in the central portions of the region. In the highly conserved sequence elements in the type II elements, the substitution rates were not uniform among lineages and acceleration in the substitution rate might have been due to loss of functional constraint in the stem–loop-forming sequences predicted in the type II elements. Patterns of nucleotide substitutions observed in the A+T-rich region were further compared with those in the coding regions and in the intergenic regions of mtDNA. Substitutions between A and T were particularly repressed in the highly conserved sequence elements and in the intergenic regions compared with those in the A+T-rich region excluding the highly conserved sequence elements and in the fourfold degenerate sites in the coding regions. The functional and structural characteristics of the A+T-rich region that might be involved in this substitutional bias are discussed.     

残锷线蛱蝶线粒体基因组全序列及其系统学意义(英文)

对残锷线蛱蝶(Parathymasulpitia)(鳞翅目:蛱蝶科)线粒体基因组全序列进行了测定。结果表明:残锷线蛱蝶线粒体基因组全序列全长为15268bp,除了在trnS1(AGN)和trnE基因之间有一段121bp长的基因间隔外,其基因的排列顺序及排列方向与大多数已测鳞翅目物种基本一致。在蛋白质编码基因中,除cox1以CGA作为其起始密码子之外,其余12个蛋白质编码基因都以标准的ATN作为起始密码子。此外,除nad4基因以单独的T为终止密码子,其余12个蛋白质编码基因都以TAA结尾。除trnS1(AGN)缺少DHU臂之外,22个tRNA基因都显示典型的三叶草形二级结构。除A+T富集区外的非编码序列中,线粒体基因组共含有11个基因间隔区。其中,最长的一个121bp的基因间隔区位于trnS1(AGN)和trnE之间,其A+T含量高达100%。另外,和其他鳞翅目物种一样,在其A+T富集区的3’端有一段长达18bp的poly-T结构。A+T富集区内部没有明显的小卫星样多拷贝重复序列,而含有一些微卫星样的重复结构。本研究基于13种蛋白编码基因序列的组合数据,用最大似然法和贝叶斯法对蛱蝶科几个主要亚科间共9个代表物种间的系统发生关系进行了分析。结果表明,本研究的结果与前人的分子系统学研究结论基本吻合(其中,线蛱蝶亚科和釉蛱蝶亚科互为姐妹群),而与形态学的研究结论不一致。     

A study of the phylogeny of Brassica rapa,B. nigra,Raphanus sativus,and their related genera using noncoding regions of chloroplast DNA

There are two evolutionary lineages in the genus Brassica: the rapa/oleracea lineage and the nigra lineage. Using nuclear DNA sequences such as the intergenic spacer between 5S rRNA genes and the internal transcribed spacer between 18S and 25S rRNA genes, we and others had previously demonstrated that Raphanus sativus is closely related to the nigra lineage. In the present study, we sequenced the chloroplast noncoding region between trnT and trnF and that between trnD and trnT in seven species and showed that R. sativus is more closely related to the rapa/oleracea lineage than to the nigra lineage. The conflicting results from nuclear DNA and chloroplast DNA support the hypothesis that Raphanus was derived from a hybridization between the rapa/oleracea and the nigra lineages. We estimated the date of this hybridization event to be 60% of the divergence time between the two Brassica lineages. In addition, the pattern and rate of nucleotide substitution were studied. There were more transversions than transitions in these noncoding regions, which have a high AT content. Furthermore, the proportion of transversions among the substitutions at a site increases with increasing A + T content of its two adjacent nucleotides. An influence of immediate 5(') pyrimidine on substitution pattern is also observed when both adjacent bases in the two DNA strands are A or T. The rate of nucleotide substitution in the trnL group I intron is only about one third of the rate in the nearby intergenic spacers in the trnT-trnF fragment. The rate of nucleotide substitution in the rapa/oleracea lineage is at least 1.5 times that in the nigra lineage.     

Low mitochondrial DNA variation among American alligators and a novel non-coding region in crocodilians

We analyzed 1317-1823 base pairs (bp) of mitochondrial DNA sequence beginning in the 5' end of cytochrome b (cyt b) and ending in the central domain of the control region for 25 American alligators (Alligator mississippiensis) and compared these to a homologous sequence from a Chinese alligator (A. sinensis). Both species share a non-coding spacer between cyt b and tRNA(Thr). Chinese alligator cyt b differs from that of the American alligator by 17.5% at the nucleotide level and 13.8% for inferred amino acids, which is consistent with their presumed ancient divergence. Only two cyt b haplotypes were detected among the 25 American alligators (693-1199 bp surveyed), with one haplotype shared among 24 individuals. One alligator from Mississippi differed from all other alligators by a single silent substitution. The control region contained only slightly more variation among the 25 American alligators, with two variable positions (624 bp surveyed), yielding three haplotypes with 22, two, and one individuals in each of these groups. Previous genetic studies examining allozymes and the proportion of variable microsatellite DNA loci also found low levels of genetic diversity in American alligators. However, in contrast with allozymes, microsatellites, and morphology, the mtDNA data shows no evidence of differentiation among populations from the extremes of the species range. These results suggest that American alligators underwent a severe population bottleneck in the late Pleistocene, resulting in nearly homogenous mtDNA among all American alligators today.     

Intraspecific DNA sequence variation of the mitochondrial control region of white sturgeon (Acipenser transmontanus)

Intraspecific sequence variation in the D-loop region of mtDNA in white sturgeon (Acipenser transmontanus), a relict North American fish species, was examined in 27 individuals from populations of the Columbia and Fraser rivers. Thirty-three varied nucleotide positions were present in a 462-nucleotide D-loop sequence, amplified using the polymerase chain reaction. Bootstrapped neighbor-joining and maximum- parsimony trees of sequences from 19 haplotypes suggest that the two populations have recently diverged. This is consistent with the hypothesis that the Columbia River, a Pleistocene refugium habitat, was the source of founders for the Fraser River after the last glacial recession. On the basis of a divergence time of 10-12 thousand years ago, the estimated substitution rate of the white sturgeon D-loop region is 1.1-1.3 x 10(-7) nucleotides/site/year, which is comparable to rates for hypervariable sequences in the human D-loop region. Furthermore, the ratio of mean percent nucleotide differences in the D- loop (2.27%) to that in whole mtDNA (0.54%, as estimated from restriction-enzyme data) is 4.3, which is similar to the fourfold-to- fivefold-higher substitution rate estimated for the human D-loop. The high nucleotide substitution rate of the hypervariable region indicates that the vertebrate D-loop has potential as a genetic marker in molecular population studies.      

High mutation rates in the mitochondrial genomes of Daphnia pulex

Despite the great utility of mitochondrial DNA (mtDNA) sequence data in population genetics and phylogenetics, key parameters describing the process of mitochondrial mutation (e.g., the rate and spectrum of mutational change) are based on few direct estimates. Furthermore, the variation in the mtDNA mutation process within species or between lineages with contrasting reproductive strategies remains poorly understood. In this study, we directly estimate the mtDNA mutation rate and spectrum using Daphnia pulex mutation-accumulation (MA) lines derived from sexual (cyclically parthenogenetic) and asexual (obligately parthenogenetic) lineages. The nearly complete mitochondrial genome sequences of 82 sexual and 47 asexual MA lines reveal high mtDNA mutation rate of 1.37 × 10(-7) and 1.73 × 10(-7) per nucleotide per generation, respectively. The Daphnia mtDNA mutation rate is among the highest in eukaryotes, and its spectrum is dominated by insertions and deletions (70%), largely due to the presence of mutational hotspots at homopolymeric nucleotide stretches. Maximum likelihood estimates of the Daphnia mitochondrial effective population size reveal that between five and ten copies of mitochondrial genomes are transmitted per female per generation. Comparison between sexual and asexual lineages reveals no statistically different mutation rates and highly similar mutation spectra.