Muscle-specific mutations accumulate with aging in critical human mtDNA control sites for replication

The recently discovered aging-dependent large accumulation of point mutations in the human fibroblast mtDNA control region raised the question of their occurrence in postmitotic tissues. In the present work, analysis of biopsied or autopsied human skeletal muscle revealed the absence or only minimal presence of those mutations. By contrast, surprisingly, most of 26 individuals 53 to 92 years old, without a known history of neuromuscular disease, exhibited at mtDNA replication control sites in muscle an accumulation of two new point mutations, i.e., A189G and T408A, which were absent or marginally present in 19 individuals younger than 34 years. These two mutations were not found in fibroblasts from 22 subjects 64 to 101 years of age (T408A), or were present only in three subjects in very low amounts (A189G). Furthermore, in several older individuals exhibiting an accumulation in muscle of one or both of these mutations, they were nearly absent in other tissues, whereas the most frequent fibroblast-specific mutation (T414G) was present in skin, but not in muscle. Among eight additional individuals exhibiting partial denervation of their biopsied muscle, four subjects >80 years old had accumulated the two muscle-specific point mutations, which were, conversely, present at only very low levels in four subjects ≤40 years old. The striking tissue specificity of the muscle mtDNA mutations detected here and their mapping at critical sites for mtDNA replication strongly point to the involvement of a specific mutagenic machinery and to the functional relevance of these mutations.

CONTRASTING PHYLOGEOGRAPHIC PATTERNS IN MITOCHONDRIAL DNA AND MICROSATELLITES: EVIDENCE OF FEMALE PHILOPATRY AND MALE-BIASED GENE FLOW AMONG REGIONAL POPULATIONS OF THE BLUE-AND-YELLOW MACAW (PSITTACIFORMES: ARA ARARAUNA) IN BRAZIL

Comparing the patterns of population differentiation among genetic markers with different modes of inheritance call provide insights into patterns of sex-biased dispersal and gene flow. The blue-and-yellow Macaw (Ara ararauna) is a Neotropical parrot with a broad geographic distribution ill South America. However, little is known about the natural history and current status Of remaining wild populations, including levels of genetic variability. The progressive decline and possible fragmentation of populations may endanger this species in the near future. We analyzed mitochondrial DNA (mtDNA) control-region sequences and six microsatellite 106 Of Blue-and-yellow Macaws sampled throughout their geographic range ill Brazil to describe population genetic Structure, to make inferences about historical demography and dispersal behavior, and to provide insight for conservation efforts. Analyses of population genetic structure based on mtDNA showed evidence of two major populations ill western and eastern Brazil that share a few low-frequency haplotypes. This phylogeographic pattern seems to have originated by the historical isolation of Blue-and-yellow Macaw populations similar to 374,000 years ago and has been maintained by restricted gene flow and female philopatry. By contrast, variation ill biparentally inherited microsatellites was not structured geographically, Male-biased dispersal and female philopatry best explain the different patterns observed in these two markers. Because females disperse less than males, the two regional populations with well-differentiated mtDNA haplogroups should be considered two different management units for conservation purposes. Received 4 November 2007 accepted 10 December 2008.

Population subdivision in marine environments: the contributions of biogeography, geographic distance, and discontinuous habitat to genetic differentiation in a blennioid

The relative importance of factors that may promote genetic differentiation in marine organisms is largely unknown. Here, contributions to population structure from biogeography, habitat distribution, and isolation by distance were investigated in Axoclinus nigricaudus, a small subtidal rock reef fish, throughout its range in the Gulf of California. A 408 basepair fragment of the mitochondrial control region was sequenced from 105 individuals. Variation was significantly partitioned between many pairs of populations. Phylogenetic analyses, hierarchical analyses of variance, and general linear models substantiated a major break between two putative biogeographic regions. This genetic discontinuity coincides with an abrupt change in ecological characteristics (including temperature and salinity) but does not coincide with known oceanographic circulation patterns. Geographic distance and the nature of habitat separating populations (continuous habitat along a shoreline, discontinuous habitat along a shoreline, and open water) also contributed to population structure in general linear model analyses. To verify that local populations are genetically stable over time, one population was resampled on four occasions over eighteen months; it showed no evidence of a temporal component to diversity. These results indicate that having a planktonic life stage does not preclude geographically partitioned genetic variation over relatively small geographic distances in marine environments. Moreover, levels of genetic differentiation among populations of Axoclinus nigricaudus cannot be explained by a single factor, but are due to the combined influences of a biogeographic boundary, habitat, and geographic distance.

Concordance between dispersal and mitochondrial gene flow: isolation by distance in a tropical teleost, Lates calcarifer (Australian barramundi)

Patterns of population subdivision and the relationship between gene flow and geographical distance in the tropical estuarine fish Lares calcarifer (Centropomidae) were investigated using mtDNA control region sequences. Sixty-three putative haplotypes were resolved from a total of 270 individuals from nine localities within three geographical regions spanning the north Australian coastline. Despite a continuous estuarine distribution throughout the sampled range, no haplotypes were shared among regions. However, within regions, common haplotypes were often shared among localities. Both sequence-based (average Phi(ST)=0.328) and haplotype-based (average Phi(ST)=0.182) population subdivision analyses indicated strong geographical structuring. Depending on the method of calculation, geographical distance explained either 79 per cent (sequence-based) or 23 per cent (haplotype-based) of the variation in mitochondrial gene flow. Such relationships suggest that genetic differentiation of L. calcarifer has been generated via isolation-by-distance, possibly in a stepping-stone fashion. This pattern of genetic structure is concordant with expectations based on the life history of L. calcarifer and direct studies of its dispersal patterns. Mitochondrial DNA variation, although generally in agreement with patterns of allozyme variation, detected population subdivision at smaller spatial scales. Our analysis of mtDNA variation in L. calcarifer confirms that population genetic models can detect population structure of not only evolutionary significance but also of demographic significance. Further, it demonstrates the power of inferring such structure from hypervariable markers, which correspond to small effective population sizes.

When oceans meet: a teleost shows secondary intergradation at an Indian-Pacific interface

The Indo-West Pacific is characterized by extraordinary marine species diversity. The evolutionary mechanisms responsible for generating this diversity remain puzzling, but are often linked to Pleistocene sea level fluctuations. The impact of these sea level changes on the population genetic architecture of the estuarine fish Lates calcarifer are investigated via a natural experiment in a region of the Indo-West Pacific known to have undergone considerable change during the Pleistocene. L. calcarifer, a coastline-restricted catadromous teleost, provides an excellent model for studying the effects of sea level change as its habitat requirements potentially make it sensitive to the region's physical history. Evidence was found for a large phylogenetic break (4% mtDNA control region; 0.47% ATPase 6 and 8) either side of the Torres Strait, which separates the Western Pacific and Indian Oceans, although some mixing of the clades was evident. This suggests clinal secondary introgression of the clades via contemporary gene flow. Further, populations on Australia's east coast appear to have passed through a bottleneck. This was linked to the historical drying of the Great Barrier Reef coastal lagoon, which resulted in a significant loss of habitat and forced retreat into isolated refugia. These results suggest that historical eustatic changes have left a significant imprint on the molecular diversity within marine species as well as among them in the Indo-West Pacific.

Strong genetic structuring in a habitat specialist, the Oxleyan Pygmy Perch Nannoperca oxleyana

This study used allozyme and mitochondrial DNA variation to examine genetic structure in the Oxleyan Pygmy Perch Nannoperca oxleyana. This small-bodied freshwater fish has a very restricted distribution occurring only in some small coastal streams in south-east Queensland and northern New South Wales. It was expected that subpopulations may contain little genetic variation and be highly differentiated from one another. The results, based on allozyme and mitochondrial DNA control region variation were in agreement with these expectations. Allozyme variation was very low overall, with only one locus showing variation at most sites. The high differentiation was because a different locus tended to be polymorphic at each site. Mitochondrial variation within sites was also low, but some sites had unique haplotypes. The patterns of similarity among mitochondrial DNA haplotypes were not as expected from geographical proximity alone. In particular, although some northern sites had unique haplotypes, four sites spread along 200 km of coastline were remarkably similar, sharing the same common haplotype at similar frequencies. We suggest that these four streams may have had a confluence relatively recently, possibly when sea levels were lower, 8000-10 000 BP.

Genetic population structure of the catadromous perciform: Macquaria novemaculeata (Percichthyidae)

Genetic population structure in the catadromous Australian bass Macquaria novemaculeata was investigated using samples from four locations spanning 600 km along the eastern Australian coastline. Both allozymes and mtDNA control region sequences were examined. Population subdivision estimates based on allozymes revealed low levels of population structuring (G(st)=0.043, P<0.05). However, mtDNA indicated moderate levels of geographic population structure (G(st)=0.146, P<0.01). Phylogenetic analysis of mtDNA control region sequences (mean sequence divergence 1.9%) indicated little phylogeographic structuring. Results suggested that genotypic variation within each river population, while bring affected primarily by genetic drift, was also prevented from more significant divergence by homogenizing levels of gene flow-synonymous with a one-dimensional stepping-stone model of population structure. The catadromous life history of Macquaria novemaculeata was considered to br influential on the pattern of population structure displayed. Results were compared to the few population genetic studies involving catadromous fishes, indicating that catadromy alone is unlikely to be a good predictor of population structure. A more comprehensive suite of biological characteristics than simple life-history traits must be considered fully to allow reliable predictive models of population structure to be formulated. (C) 1997 The Fisheries Society of the British Isles.

The GHEP-EMPOP collaboration on mtDNA population data-A new resource for forensic casework

Mitochondrial DNA (mtDNA) population data for forensic purposes are still scarce for some populations, which may limit the evaluation of forensic evidence especially when the rarity of a haplotype needs to be determined in a database search. In order to improve the collection of mtDNA lineages from the Iberian and South American subcontinents, we here report the results of a collaborative study involving nine laboratories from the Spanish and Portuguese Speaking Working Group of the International Society for Forensic Genetics (GHEP-ISFG) and EMPOP. The individual laboratories contributed population data that were generated throughout the past 10 years, but in the majority of cases have not been made available to the scientific community. A total of 1019 haplotypes from Iberia (Basque Country, 2 general Spanish populations, 2 North and 1 Central Portugal populations), and Latin America (3 populations from Sao Paulo) were collected, reviewed and harmonized according to defined EMPOP criteria. The majority of data ambiguities that were found during the reviewing process (41 in total) were transcription errors confirming that the documentation process is still the most error-prone stage in reporting mtDNA population data, especially when performed manually. This GHEP-EMPOP collaboration has significantly improved the quality of the individual mtDNA datasets and adds mtDNA population data as valuable resource to the EMPOP database (www.empop.org). (C) 2010 Elsevier Ireland Ltd. All rights reserved.

Numerous gene rearrangements in the mitochondrial genome of the wallaby louse, Heterodoxus macropus (Phthiraptera)

The complete arrangement of genes in the mitochondrial (mt) genome is known for 12 species of insects, and part of the gene arrangement in the mt genome is known for over 300 other species of insects. The arrangement of genes in the mt genome is very conserved in insects studied, since all of the protein-coding and rRNA genes and most of the tRNA genes are arranged in the same way. We sequenced the entire mt genome of the wallaby louse, Heterodoxus macropus, which is 14,670 bp long and has the 37 genes typical of animals and some noncoding regions. The largest noncoding region is 73 bp long (93% A+T), and the second largest is 47 bp long (92% AST). Both of these noncoding regions seem to be able to form stem-loop structures. The arrangement of genes in the mt genome of this louse is unlike that of any other animal studied. All tRNA genes have moved and/or inverted relative to the ancestral gene arrangement of insects, which is present in the fruit fly Drosophila yakuba. At least nine protein-coding genes (atp6, atp8, cox2, cob, nad1-nad3, nad5, and nad6) have moved; moreover, four of these genes (atp6, atp8, nad1, and nad3) have inverted. The large number of gene rearrangements in the mt genome of H. macropus is unprecedented for an arthropod.

Complete DNA sequence and gene organization of the mitochondrial genome of the liverfluke, Fasciola hepatica L. (Platyhelminthes; Trematoda)

The complete nucleotide sequence of the mitochondrial (mt) DNA molecule of the liverfluke, Fasciola hepatica (phylum Platyhelminthes, class Trematoda, family Fasciolidae), was determined, It comprises 14462 bp, contains 12 protein-encoding, 2 ribosomal and 22 transfer RNA genes, and is the second complete flatworm (and the first trematode) mitochondrial sequence to be described in detail. All of the genes are transcribed from the same strand. Of the genes typically found in mitochondrial genomes of eumetazoans, only atp8 is absent. The nad4L and nad4 genes overlap by 40 nt. Most intergenic sequences are very short. Two larger non-coding regions are present. The longer one (817 nt) is located between trnG and cox3 and consists of 8 identical tandem repeats of 85 nt, rich in G and C, followed by 1 imperfect repeat. The shorter non-coding region (187 nt) exhibits no special features and is separated from the longer region by trnG. The gene arrangement resembles that of some other trematodes including the eastern Asian Schistosoma species (and cyclophyllidean cestode species) but it is strikingly different from that of the African schistosomes, represented by Schistosoma mansoni. The genetic code is as inferred previously for flatworms. Transfer RNA genes range in length from 58 to 70 nt, their products producing characteristic 'clover leaf' structures, except for tRNA(S-VON) and tRNA(S-AGN) lacking the DHU arm.

Oceanic interchange and nonequilibrium population structure in the estuarine dependent Indo-Pacific tasselfish, Polynemus sheridani

We assayed mtDNA haplotype [300 base pairs (bp) control region] geography and genealogy in the Indo-Pacific tasselfish, Polynemus sheridani from its contiguous estuarine distribution across northern Australia (n = 169). Eight estuaries were sampled from three oceanographic regions (Timor Sea, Gulf of Carpentaria and the Coral Sea) to assess the impact of Pleistocene sea level changes on the historical connectivity among P. sheridani populations. Specifically, we investigated the genetic consequences of disruption to Indian-Pacific Ocean connectivity brought about by the closure of the Torres Strait. Overall there was significant population subdivision among estuaries (F-ST = 0.161, (Phi(ST) = 0.187). Despite a linear distribution, P. sheridani did not show isolation by distance over the entire sampled range because of genetic similarity of estuaries greater than 3000 km apart. However, significant isolation by distance was detected between estuaries separated by less than 3000 km of coastline. Unlike many genetic studies of Indo-Pacific marine species, there was no evidence for an historical division between eastern and western populations. Instead, phylogeographical patterns were dominated by a starlike intraspecific phylogeny coupled with evidence for population expansion in both the Gulf of Carpentaria and the Coral Sea but not the Timor Sea. This was interpreted as evidence for recent west to east recolonization across of northern Australia following the last postglacial marine advance. We argue that although sufficient time has elapsed postcolonization for populations to approach gene flow-drift equilibrium over smaller spatial scales (< 3000 km), the signal of historical colonization persists to obscure the expected equilibrium pattern of isolation by distance over large spatial scales (> 3000 km).

The value of idiosyncratic markers and changes to conserved tRNA sequences from the mitochondrial genome of hard ticks (Acari: Ixodida: Ixodidae) for phylogenetic inference

Idiosyncratic markers are features of genes and genomes that are so unusual that it is unlikely that they evolved more than once in a lineage of organisms. Here we explore further the potential of idiosyncratic markers and changes to typically conserved tRNA sequences for phylogenetic inference. Hard ticks were chosen as the model group because their phylogeny has been studied extensively. Fifty-eight candidate markers from hard ticks ( family Ixodidae) and 22 markers from the subfamily Rhipicephalinae sensu lato were mapped onto phylogenies of these groups. Two of the most interesting markers, features of the secondary structure of two different tRNAs, gave strong support to the hypothesis that species of the Prostriata ( Ixodes spp.) are monophyletic. Previous analyses of genes and morphology did not strongly support this relationship, instead suggesting that the Prostriata is paraphyletic with respect to the Metastriata ( the rest of the hard ticks). Parallel or convergent evolution was not found in the arrangements of mitochondrial genes in ticks nor were there any reversals to the ancestral arthropod character state. Many of the markers identified were phylogenetically informative, whereas others should be informative with study of additional taxa. Idiosyncratic markers and changes to typically conserved nucleotides in tRNAs that are phylogenetically informative were common in this data set, and thus these types of markers might be found in other organisms.

The highly rearranged mitochondrial genome of the plague thrips, Thrips imaginis (Insecta : thysanoptera): Convergence of two novel gene boundaries and an extraordinary arrangement of rRNA genes

To help understand the mechanisms of gene rearrangement in the mitochondrial (mt) genomes of hemipteroid insects, we sequenced the mt genome of the plague thrips, Thrips imaginis (Thysanoptera). This genome is circular, 15,407 by long, and has many unusual features, including (1) rRNA genes inverted and distant from one another, (2) an extra gene for tRNA-Ser, (3) a tRNA-Val lacking a D-arm, (4) two pseudo-tRNA genes, (5) duplicate control regions, and (6) translocations and/or inversions of 24 of the 37 genes. The mechanism of rRNA gene transcription in T. imaginis may be different from that of other arthropods since the two rRNA genes have inverted and are distant from one another. Further, the rRNA genes are not adjacent or even close to either of the two control regions. Tandem duplication and deletion is a plausible model for the evolution of duplicate control regions and for the gene translocations, but intramitochondrial recombination may account for the gene inversions in T. imaginis. All the 18 genes between control regions #1 and #2 have translocated and/or inverted, whereas only six of the 20 genes outside this region have translocated and/or inverted. Moreover, the extra tRNA gene and the two pseudo-tRNA genes are either in this region or immediately adjacent to one of the control regions. These observations suggest that tandem duplication and deletion may be facilitated by the duplicate control regions and may have occurred a number of times in the lineage leading to T. imaginis. T. imaginis shares two novel gene boundaries with a lepidopsocid species from another order of hemipteroid insects, the Psocoptera. The evidence available suggests that these shared gene boundaries evolved by convergence and thus are not informative for the interordinal phylogeny of hemipteroid insects. We discuss the potential of hemipteroid insects as a model system for studies of the evolution of animal rut genomes and outline some fundamental questions that may be addressed with this system.

Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects

A number of studies indicated that lineages of animals with high rates of mitochondrial (mt) gene rearrangement might have high rates of mt nucleotide substitution. We chose the hemipteroid assemblage and the Insecta to test the idea that rates of mt gene rearrangement and mt nucleotide substitution are correlated. For this purpose, we sequenced the mt genome of a lepidopsocid from the Psocoptera, the only order of hemipteroid insects for which an entire mtDNA sequence is not available. The mt genome of this lepidopsocid is circular, 16,924 bp long, and contains 37 genes and a putative control region; seven tRNA genes and a protein-coding gene in this genome have changed positions relative to the ancestral arrangement of mt genes of insects. We then compared the relative rates of nucleotide substitution among species from each of the four orders of hemipteroid insects and among the 20 insects whose mt genomes have been sequenced entirely. All comparisons among the hernipteroid insects showed that species with higher rates of gene rearrangement also had significantly higher rates of nucleotide substitution statistically than did species with lower rates of gene rearrangement. In comparisons among the 20 insects, where the mt genomes of the two species differed by more than five breakpoints, the more rearranged species always had a significantly higher rate of nucleotide substitution than the less rearranged species. However, in comparisons where the mt genomes of two species differed by five or less breakpoints, the more rearranged species did not always have a significantly higher rate of nucleotide substitution than the less rearranged species. We tested the statistical significance of the correlation between the rates of mt gene rearrangement and mt nucleotide substitution with nine pairs of insects that were phylogenetically independent from one 2 another. We found that the correlation was positive and statistically significant (R-2 = 0.73, P = 0.01; R-s = 0.67, P < 0.05). We propose that increased rates of nucleotide substitution may lead to increased rates of gene rearrangement in the mt genomes of insects.

Faeces as a source of DNA for molecular studies in a threatened population of great bustards

With recent advances in molecular biology, it is now possible to use the trace amounts of DNA in faeces to non-invasively sample endangered species for genetic studies. A highly vulnerable population of approximately 100 great bustards (Otis tarda) exists in Morocco necessitating the use of non-invasive protocols to study their genetic structure. Here we report a reliable silica-based method to extract DNA from great bustard faeces. We found that successful extraction and amplification correlated strongly with faeces freshness and composition. We could not extract amplifiable DNA from 30% of our samples as they were dry or contained insect material. However 100% of our fresh faecal samples containing no obvious insect material worked, allowing us to assess the levels of genetic variation among 25 individuals using a 542 bp control region sequence. We were able to extract DNA from four out of five other avian species, demonstrating that faeces represents a suitable source of DNA for population genetics studies in a broad range of species.