Genetic variation and bill size dimorphism in a passerine bird, the reed bunting Emberiza schoeniclus

In passerine birds morphological differentiation in bill size within species is not commonly observed. Bill size is usually associated with a trophic niche, and strong differences in it may reflect the process of genetic differentiation and, possibly, speciation. We used both mitochondrial DNA (mtDNA) and nuclear microsatellites to study genetic variation between two subspecies of reed bunting, Emberiza schoeniclus schoeniclus and E.s. intermedia, along their distributional boundary in western Europe. These two subspecies are characterized by a high dimorphism in bill size and, although breeding populations of the two subspecies are found very close to each other in northern Italy, apparently no interbreeding occurs. The observed morphological pattern between the two subspecies may be maintained by geographically varying selective forces or, alternatively, may be the result of a long geographical separation followed by a secondary contact. MtDNA sequences of cytochrome b and ND5 (515 bp) showed little variation and did not discriminate between the two subspecies, indicating a divergence time of less than 500 000 years. The analysis of four microsatellite loci suggested a clear, although weak, degree of genetic differentiation in the large- and small-billed populations, as indicated by FST and RST values and genetic distances. The correlation between bill size and genetic distance between populations remained significant after accounting for the geographical distances between sampling localities. Altogether, these results indicate a very recent genetic differentiation between the two bill morphs and suggest that a strong selection for large bills in the southern part of the breeding range is probably involved in maintaining the geographical differentiation of this species.     

The molecular evolution of sperm bindin in six species of sea urchins (Echinoida: Strongylocentrotidae)

The acrosomal protein bindin attaches sperm to eggs during sea urchin fertilization. Complementary to ongoing functional biochemical studies, I take a comparative approach to explore the molecular evolution of bindin in a group of closely related free-spawning echinoid species. Two alleles of the mature bindin gene were sequenced for each of six species in the sea urchin family Strongylocentrotidae. The nucleotide sequences diverged by at least 1% per Myr at both silent and replacement sites. Two short sections flanking the conserved block show an excess of nonsynonymous substitutions. Each is homologous to a region that had been identified as a target of selection in other sea urchin comparisons. A large proportion of the bindin-coding sequence consists of a highly variable repeat region. Bindin sequences, even including the large intron, could not resolve the branching order among five of the species.     

Determinants of postprandial appetite sensations: macronutrient intake and glucose metabolism

Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis was used to characterize mitochondrial DNA (mtDNA) variation in screwworms, Cochliomyia hominivorax, and secondary screwworm, C.macellaria, from the Caribbean, North America and South America. Four amplicons, totaling 7.1 kb, were analysed with sixteen restriction enzymes. A total of 133 restriction sites was observed in the two species, 104 in C.hominivorax, of which nineteen were variable, and ninety-five in C.macellaria, none of which was variable. Fourteen mtDNA haplotypes were observed among eighteen C.hominivorax examined. Mean divergence between C.hominivorax haplotypes (d) was 0.0064 substitutions per base-pair and genotypic diversity (G) was 0.97. Mean divergence between C.hominivorax and C.macellaria was 0.0824. Cochliomyia hominivorax haplotypes could be divided into three assemblages representing North America, South America and Jamaica, based on UPGMA clustering with d values. The assemblages did not exhibit complete geographic fidelity. These data were discordant with previously published allozyme data indicating little differentiation between screwworm populations. A scenario invoking historically isolated populations coming into contact with the introduction and movement of European livestock is proposed to explain the observed population structure of screwworm.     

Contrasting population structure from nuclear intron sequences and mtDNA of humpback whales

Powerful analyses of population structure require information from multiple genetic loci. To help develop a molecular toolbox for obtaining this information, we have designed universal oligonucleotide primers that span conserved intron-exon junctions in a wide variety of animal phyla. We test the utility of exon-primed, intron-crossing amplifications by analyzing the variability of actin intron sequences from humpback, blue, and bowhead whales and comparing the results with mitochondrial DNA (mtDNA) haplotype data. Humpback actin introns fall into two major clades that exist in different frequencies in different oceanic populations. It is surprising that Hawaii and California populations, which are very distinct in mtDNAs, are similar in actin intron alleles. This discrepancy between mtDNA and nuclear DNA results may be due either to differences in genetic drift in mitochondrial and nuclear genes or to preferential movement of males, which do not transmit mtDNA to offspring, between separate breeding grounds. Opposing mtDNA and nuclear DNA results can help clarify otherwise hidden patterns of structure in natural populations.     

Molecular phylogeny and ecological diversification in a clade of New World songbirds (genus Vireo)

We constructed a molecular phylogeny for a clade of eye-ringed vireos (Vireo flavifrons and the V. solitarius complex) to examine existing hypotheses of speciation and ecological diversification. Complete sequences of the mtDNA cytochrome b gene were obtained from 47 individuals of this group plus four vireonid outgroups. Mean levels of sequence divergence in the clade varied from 0.29% to 5.7%. Differences were greatest between V. flavifrons and four taxa of 'V. solitarius'. The latter separated into three taxonomic, geographical and ecological groups: V. plumbeus plumbeus, V. cassinii cassinii, and V. solitarius solitarius plus V. solitarius alticola. These differed by an average of 2.6-3.2%. Populations within each group revealed low levels of sequence variation (x = 0.20%) and little geographical structuring. The mtDNA data generally corroborate results from allozymes. V. plumbeus shows a loss of yellow-green carotenoid pigmentation from the ancestral condition. The occupancy of relatively dry habitats by this species and V. cassinii represents a derived ecological shift from more-humid environments occupied by other species of vireonids. Ecological divergence in this clade occurred in allopatry and is associated with generic-level stability in morphometrics and foraging styles. Migratory behaviour and seasonal habitat shifts apparently evolved multiple times in vireos breeding in temperate environments. Present geographical and ecological distributions, and low levels of intrataxon genetic divergence, are hypothesized to be the result of postglacial regionalization of climate-plant associations and rapid northward expansion of breeding ranges.     

Phylogeography of Bufo marinus from its natural and introduced ranges

The marine toad, Bufo marinus, has a broad natural distribution extending from the south-west of the USA to southern Peru and the central Amazon. It was introduced to several localities in the Caribbean and Pacific Oceans to control sugar cane pests. We sequenced 468 bp of mitochondrial DNA (mtDNA) containing the ND3 gene, and flanking tRNA genes from toads spanning the broad natural and introduced ranges. Consistent with the known history of introductions and expected effects of serial bottlenecks, mtDNA within introduced populations in Hawaii and Australia was uniform and most closely related to samples from eastern Venezuela and French Guiana. However, mtDNA nucleotide diversity in the geographic region spanning the source areas is also relative low (0.18-0.46%) and the absence of variation in the introduced populations precludes quantitative assessment of the reduction in genetic diversity. Unexpectedly, there was a large phylogeographic break (5.4% sequence divergence) within the natural range separating populations east and west of the Venezuelan Andes. We hypothesize that the two major lineages of B. marinus were isolated by the uplift of the eastern Andean cordillera which was completed approximately 2.7 Ma. Another species of the marinus group, B. paracnemis, had mtDNA paraphyletic, with marinus, being nested within the eastern lineage. Thus, at least one speciation event within the marinus group postdates the split within marinus. These findings suggest that the taxonomy of B. marinus should be re-evaluated and that the search for pathogens to control Australian populations should be conducted in populations from both lineages in the natural range.     

Cytochrome b gene haplotypes characterize chromosomal lineages of anoa, the Sulawesi dwarf buffalo (Bovidae: Bubalus sp.)

Partial mitochondrial cytochrome b gene sequences reveal two deeply differentiated mtDNA lineages in anoa dwarf buffaloes (Bubalus depressicornis) from the studbook herd in European zoos. Three matrilinear lineages of lowland anoas (depressicornis type) contributed three rather similar sequence haplotypes, but one remarkably distinct haplotype was observed exclusively in mountain anoas (quarlesi type) descended from one founder female. The carriers of the distinctive mtDNA haplotype were also distinguished by several chromosomal and phenotypic peculiarities too. The differentiation between the mtDNA lineages of anoa approached or even surpassed the genetic divergence between some uncontested species of wild cattle. The depth of this haplotype divergence in anoas is discussed against the background of the phylogenetic age of these paleoendemic inhabitants of a predator-free island refugium, Sulawesi, who are among the most plesiomorphic living bovines. The studbook breeding of captive anoas as a safeguard against extinction might profit from such population genetic markers. These cytochrome b gene sequences were unable to resolve the phylogeny of nine bovine taxa robustly, except the divergence of Bubalus, Synceros, Bison, and Bos (sensu lato) genera.     

Habitat islands, genetic diversity, and gene flow in a Patagonian rodent

The effects of terrestrial habitat islands on gene flow and genetic diversity in animal populations have been predicted and discussed in theoretical terms, but empirical data are needed to test these predictions and provide an understanding of the relationships of life-history characteristics to genetics of insular species. We studied saxicolous mice (Phyllotis xanthopygus) in Patagonia to explore genetic structure, phylogeography, and gene flow in a species inhabiting natural habitat islands. Phylogeographic analyses based on mtDNA sequences revealed two haplotype clades, which presumably reflect early Pleistocene factors that temporarily separated the mice into two geographically isolated groups. The Río Chubut, which lies within a glacial drainage basin bisecting northern Patagonia, might have affected gene flow in the species. Although we anticipated isolation by distance and founder phenomena associated with habitat islands, in some habitat patches we found evidence of high local genetic diversity. The amount of divergence in the mitochondrial cytochrome b gene (approximately 3.4%) in animals at a single locality could best be explained through a combination of historical factors and metapopulation source-sink theory. Demographic shifts, dispersal, and episodic recolonization are important in the life history and genetic population structure of P. xanthopygus.     

Population genetic structure of North Atlantic, Mediterranean Sea and Sea of Cortez fin whales, Balaenoptera physalus (Linnaeus 1758): analysis of mitochondrial and nuclear loci

Samples were collected from 407 fin whales, Balaenoptera physalus, at four North Atlantic and one Mediterranean Sea summer feeding area as well as the Sea of Cortez in the Pacific Ocean. For each sample, the sex, the sequence of the first 288 nucleotides of the mitochondrial (mt) control region and the genotype at six microsatellite loci were determined. A significant degree of divergence was detected at all nuclear and mt loci between North Atlantic/Mediterranean Sea and the Sea of Cortez. However, the divergence time estimated from the mt sequences was substantially lower than the time elapsed since the rise of the Panama Isthmus, suggesting occasional gene flow between the North Pacific and North Atlantic ocean after the separation of the two oceans. Within the North Atlantic and Mediterranean Sea, significant levels of heterogeneity were observed in the mtDNA between the Mediterranean Sea, the eastern (Spain) and the western (the Gulf of Maine and the Gulf of St Lawrence) North Atlantic. Samples collected off West Greenland and Iceland could not be unequivocally assigned to either of the two areas. The homogeneity tests performed using the nuclear data revealed significant levels of divergence only between the Mediterranean Sea and the Gulf of St Lawrence or West Greenland. In conclusion, our results suggest the existence of several recently diverged populations in the North Atlantic and Mediterranean Sea, possibly with some limited gene flow between adjacent populations, a population structure which is consistent with earlier population models proposed by Kellogg, Ingebrigtsen, and Sergeant.     

The effect of mating system differences on nucleotide diversity at the phosphoglucose isomerase locus in the plant genus Leavenworthia

To test the theoretical prediction that highly inbreeding populations should have low neutral genetic diversity relative to closely related outcrossing populations, we sequenced portions of the cytosolic phosphoglucose isomerase (PgiC) gene in the plant genus Leavenworthia, which includes both self-incompatible and inbreeding taxa. On the basis of sequences of intron 12 of this gene, the expected low diversity was seen in both populations of the selfers Leavenworthia uniflora and L. torulosa and in three highly inbreeding populations of L. crassa, while high diversity was found in self-incompatible L. stylosa, and moderate diversity in L. crassa populations with partial or complete self-incompatibility. In L. stylosa, the nucleotide diversity was strongly structured into three haplotypic classes, differing by several insertion/deletion sequences, with linkage disequilibrium between sequences of the three types in intron 12, but not in the adjacent regions. Differences between the three kinds of haplotypes are larger than between sequences of this gene region from different species. The haplotype divergence suggests the presence of a balanced polymorphism at this locus, possibly predating the split between L. stylosa and its two inbreeding sister taxa, L. uniflora and L. torulosa. It is therefore difficult to distinguish between different potential causes of the much lower sequence diversity at this locus in inbreeding than outcrossing populations. Selective sweeps during the evolution of these populations are possible, or background selection, or merely loss of a balanced polymorphism maintained by overdominance in the populations that evolved high selfing rates.     

Transposition of group II intron aI1 in yeast and invasion of mitochondrial genes at new locations

Intron mobility at the RNA level by splicing reversal at allelic (homing) and non-allelic locations (transposition) has been reported in vitro. In the living cell, however, only intron homing by unidirectional gene conversion has been described. Supposing that intron insertions at non-allelic sites might occur in vivo, we speculated that group II splice-site-associated macro-deletions in fungal mitochondrial DNA might result from group II intron transposition to new locations followed by recombination. We used polymerase chain reaction techniques to detect this critical, infrequent intermediate in mtDNA populations. Here we report on group II intron aI1 transposition to non-allelic, splicing-compatible locations within the cox1 gene of yeast mtDNA. The identified integration sites are preceded by motifs similar to the upstream exon A1. Sequences flanking intron aI1 are not co-converted to the insertion sites and cis- and trans-acting mutations within aI1 reduce intron mobility below detection levels. These findings suggest the involvement of an RNA intermediate in group II intron transposition.     

Specific binding of aminoglycoside antibiotics to RNA

BACKGROUND: Aminoglycoside antibiotics interfere with ribosomal protein synthesis and with intron splicing. Various lines of evidence suggest that RNA is the molecular target for aminoglycosides, but little is known about the recognition process. Is recognition of a particular aminoglycoside specific for certain RNA structures? If so, what are the rules for recognition? We have begun to investigate this problem by in vitro selection of RNA molecules that can specifically bind to the aminoglycoside antibiotic tobramycin. RESULTS: An RNA diversity library was used to select for sequences capable of binding to the aminoglycoside antibiotic tobramycin. After six cycles of selection, 82% of the RNA bound to tobramycin specifically. The selected RNA was reverse-transcribed into DNA, which was then cloned. At low selection stringency, an extremely large number of clones, on the order of 10(7), produced RNAs capable of binding tobramycin with Kds in the microM range (values similar to that observed for the binding of tobramycin to Escherichia coli ribosomes). Sequencing of 18 of the clones revealed no obvious consensus sequence. At higher selection stringencies (Kds in the nM range) only two consensus sequences for binding were observed. CONCLUSIONS: We have shown that RNA molecules can be readily selected that bind the aminoglycoside tobramycin. The RNAs that bind tobramycin with high affinity contain consensus binding regions that may be confined to predicted stem-loop structures. These studies open the way for understanding the basis of RNA-aminoglycoside recognition.     

Function of proteoglycans in the extracellular matrix

Estimating phylogenies from DNA sequence data has become the major methodology of molecular phylogenetics. To date, molecular phylogenetics of the vertebrates has been very dependent on mtDNA, but studies involving mtDNA are limited because the several genes comprising the mt-genome are inherited as a single linkage group. The only apparent solution to this problem is to sequence additional genes, each representing a distinct linkage group, so that the resultant gene trees provide independent estimates of the species tree. There exists the need to find novel gene sequences which contain enough phylogenetic information to resolve relationships between closely related species. A possible source is the nuclear-encoded introns, because they evolve more rapidly than exons. We designed primers to amplify and sequence the 7 intron from the beta-fibrinogen gene for a recently evolved group, the woodpeckers. We sequenced the entire intron for 10 specimens representing five species. Nucleotide substitutions are randomly distributed along the length of the intron, suggesting selective neutrality. A preliminary analysis indicates that the phylogenetic signal in the intron is as strong as that in the mitochondrial encoded cytochrome b (cyt b) gene. The topology of the beta-fibrinogen tree is identical to that of the cyt b tree. This analysis demonstrates the ability of the 7 intron of beta-fibrinogen to provide well resolved, independent gene trees for recently evolved groups and establishes it as a source of sequences to be used in other phylogenetic studies.     

Concerted evolution at the population level: pupfish HindIII satellite DNA sequences

The canonical monomers (approximately 170 bp) of an abundant (1.9 x 10(6) copies per diploid genome) satellite DNA sequence family in the genome of Cyprinodon variegatus, a "pupfish" that ranges along the Atlantic coast from Cape Cod to central Mexico, are divergent in base sequence in 10 of 12 samples collected from natural populations. The divergence involves substitutions, deletions, and insertions, is marked in scope (mean pairwise sequence similarity = 61.6%; range = 35-95.9%), is largely confined to the 3' half of the monomer, and is not correlated with the distance among collecting sites. Repetitive cloning and direct genomic sequencing experiments failed to detect intrapopulation and intraindividual variation, suggesting high levels of sequence homogeneity within populations. The satellite sequence has therefore undergone "concerted evolution," at the level of the local population. Concerted evolution has previously almost always been discussed in terms of the divergence of species or higher taxa; its intraspecific occurrence apparently has not been reported previously. The generality of the observation is difficult to evaluate, for although satellite DNAs from a large number of organisms have been studied in detail, there appear to be little or no other data on their sequence variation in natural populations. The relationship (if any) between concerted, population level, satellite DNA divergence and the extent of gene flow/genetic isolation among conspecific natural populations remains to be established.     

Different forces drive the evolution of the Acp26Aa and Acp26Ab accessory gland genes in the Drosophila melanogaster species complex

The Acp26Aa and Acp26Ab genes that code for male accessory gland proteins are tandemly arranged in the species of the Drosophila melanogaster complex. An approximately 1.6-kb region encompassing both genes has been sequenced in 10, 24, and 18 lines from Spain, Ivory Coast, and Malawi, respectively; the previously studied 10 lines from North Carolina have also been included in the analyses. A total of 110 nucleotide and 4 length polymorphisms were detected. Silent variation for the whole Acp26A region was slightly higher in African than in non-African populations, while for both genes nonsynonymous variation was similar in all populations studied. Based on Fst estimates no major genetic differentiation was detected between East and West Africa, while in general non-African populations were strongly differentiated from both African populations. Comparison of polymorphism and divergence at synonymous and nonsynonymous sites revealed that directional selection acting on amino acid replacement changes has driven the evolution of the Acp26Aa protein in the last 2.5 myr.     

Structural and evolutionary studies on sterol 14-demethylase P450 (CYP51), the most conserved P450 monooxygenase: II. Evolutionary analysis of protein and gene structures

Phylogenetic analyses based on protein sequence data indicated that sterol 14-demethylase P450 (CYP51) and bacterial CYP51-like protein were joined into a distinctive evolutionary cluster, CYP51 cluster, within the CYP protein superfamily. The most probable branch topology of the CYP51 phylogenetic tree was (bacteria, (plants, (fungi, mammals))), which is comparable to the phylogeny of major kingdoms of living matter, suggesting that CYP51 has been conserved from the era of prokaryotic evolution. This may be strong evidence supporting the prokaryotic origin of P450. Structure of flanking regions and the number and insertion sites of introns are quite different between mammalian and fungal CYP51s. This fact indicates that different mechanisms are operative in evolution of protein sequences and gene structures. CYP51 is the first example violating the well-documented rule that the basic structure of a gene, including intron insertion sites, is well conserved in each P450 family. One CYP51 processed a pseudogene was found in rat genome. Nonsynonymous nucleotide divergence observed between the pseudogene and CYP51 cDNA was less than one-fifth of the synonymous divergence. This unusually low rate of nonsynonymous nucleotide changes in the pseudogene suggests that it may be derived from another CYP51, which might have been active for a significant duration in the past.     

Molecular evolution of the Amy multigenes in the subgenus Sophophora of Drosophila

Molecular evolution of the Amy multigenes in Drosophila was investigated using PCR amplification. Twenty-five partial Amy sequences from 13 species belonging mainly to the subgenus Sophophora were determine, and a molecular phylogeny of the Amy genes in Drosophila was constructed, together with published Amy sequences. Clustering of species are mostly consistent with the traditional classification and that inferred from other genes. From sequence divergence between PCR products, several species, including D. elegans and D. fuyamai, were suggested to have multiple copies of the Amy genes. The loss of an intron took place at least three times after the Sophophora radiation. In order to investigate the mechanism of sequence evolution, the numbers of amino acid replacement and synonymous substitutions in five lineages were estimated. The heterogeneity in the relative numbers of synonymous and replacement substitutions among the lineages was found. Possible roles of selection in the sequence evolution of the Amy gene are discussed.     

Genetic variability in mitochondrial DNA of the screwworm, Cochliomyia hominivorax (Diptera: Calliphoridae), from Brazil

Restriction fragment length polymorphism (RFLP) analysis of mitochondrial DNA (mtDNA) was used to examine genetic variation and population structure of screwworm flies in four populations from S?o Paulo State, Brazil. The total DNA of 405 individuals was digested with 15 restriction endonucleases and probed with five cloned HindIII fragments representing the entire mitochondrial genome of Cochliomyia hominivorax. The survey revealed that four enzymes (HaeIII, HindIII, MspI, and PvuII) were suitable to detect mtDNA variation among all populations. Based on the fragment patterns obtained for these four enzymes, a total of 15 haplotypes in combination was detected. Heteroplasmic individuals for the PvuII pattern were obtained in one of the populations. The estimated average for nucleotide sequence divergence (delta) was 0.92%. The cladogram of the geographical distribution among the observed haplotypes suggests that the sampled screwworms probably belong to a single evolutionary lineage with populations interconnected by reduced gene flow.