Influence of historical landscapes, drainage evolution and ecological traits on patterns of genetic diversity in Southeast Asian freshwater snakehead fishes

Snakehead fishes in the family Channidae are obligate freshwater fishes represented by two extant genera, the African Parachannna and the Asian Channa. These species prefer still or slow flowing water bodies, where they are top predators that exercise high levels of parental care, have the ability to breathe air, can tolerate poor water quality, and interestingly, can aestivate or traverse terrestrial habitat in response to seasonal changes in freshwater habitat availability. These attributes suggest that snakehead fishes may possess high dispersal potential, irrespective of the terrestrial barriers that would otherwise constrain the distribution of most freshwater fishes. A number of biogeographical hypotheses have been developed to account for the modern distributions of snakehead fishes across two continents, including ancient vicariance during Gondwanan break-up, or recent colonisation tracking the formation of suitable climatic conditions. Taxonomic uncertainty also surrounds some members of the Channa genus, as geographical distributions for some taxa across southern and Southeast (SE) Asia are very large, and in one case is highly disjunct. The current study adopted a molecular genetics approach to gain an understanding of the evolution of this group of fishes, and in particular how the phylogeography of two Asian species may have been influenced by contemporary versus historical levels of dispersal and vicariance. First, a molecular phylogeny was constructed based on multiple DNA loci and calibrated with fossil evidence to provide a dated chronology of divergence events among extant species, and also within species with widespread geographical distributions. The data provide strong evidence that trans-continental distribution of the Channidae arose as a result of dispersal out of Asia and into Africa in the mid–Eocene. Among Asian Channa, deep divergence among lineages indicates that the Oligocene-Miocene boundary was a time of significant species radiation, potentially associated with historical changes in climate and drainage geomorphology. Mid-Miocene divergence among lineages suggests that a taxonomic revision is warranted for two taxa. Deep intra-specific divergence (~8Mya) was also detected between C. striata lineages that occur sympatrically in the Mekong River Basin. The study then examined the phylogeography and population structure of two major taxa, Channa striata (the chevron snakehead) and the C. micropeltes (the giant snakehead), across SE Asia. Species specific microsatellite loci were developed and used in addition to a mitochondrial DNA marker (Cyt b) to screen neutral genetic variation within and among wild populations. C. striata individuals were sampled across SE Asia (n=988), with the major focus being the Mekong Basin, which is the largest drainage basin in the region. The distributions of two divergent lineages were identified and admixture analysis showed that where they co-occur they are interbreeding, indicating that after long periods of evolution in isolation, divergence has not resulted in reproductive isolation. One lineage is predominantly confined to upland areas of northern Lao PDR to the north of the Khorat Plateau, while the other, which is more closely related to individuals from southern India, has a widespread distribution across mainland SE Asian and Sumatra. The phylogeographical pattern recovered is associated with past river networks, and high diversity and divergence among all populations sampled reveal that contemporary dispersal is very low for this taxon, even where populations occur in contiguous freshwater habitats. C. micropeltes (n=280) were also sampled from across the Mekong River Basin, focusing on the lower basin where it constitutes an important wild fishery resource. In comparison with C. striata, allelic diversity and genetic divergence among populations were extremely low, suggesting very recent colonisation of the greater Mekong region. Populations were significantly structured into at least three discrete populations in the lower Mekong. Results of this study have implications for establishing effective conservation plans for managing both species, that represent economically important wild fishery resources for the region. For C. micropeltes, it is likely that a single fisheries stock in the Tonle Sap Great Lake is being exploited by multiple fisheries operations, and future management initiatives for this species in this region will need to account for this. For C. striata, conservation of natural levels of genetic variation will require management initiatives designed to promote population persistence at very localised spatial scales, as the high level of population structuring uncovered for this species indicates that significant unique diversity is present at this fine spatial scale.

An assessment of the geographical scale of recurrent gene flow in wild populations of two species of Mekong River carps (Henicorhynchus spp.)

The Mekong is the most productive river fishery in the world, and such as, the Mekong River Basin (MRB) is very important to very large human populations across the region as a source of revenue (through fishing and marketing of aquatic resources products) and as the major source for local animal protein. Threats to biodiversity in the MRB, either to the fishery sector itself or to other sectors are a major concern, even though currently, fisheries across this region are still very productive. If not managed properly however, fish population declines will cause significant economic impact and affect livelihoods of local people and will have a major impact on food security and nutrition. Biodiversity declines will undoubtedly affect food security, income and socio-economic status of people in the MRB that depend on aquatic resources. This is an indicator of unsustainable development and hence should be avoided. Genetic diversity (biodiversity) that can be measured using techniques based on DNA markers; refers to variation within and among populations within the same species or reproductive units. In a population, new genetic variation is generated by sexual recombination contributed by individuals with mutations in genes and chromosomes. Over time, populations of a species that are not reproducing together will diverge as differential impacts of selection and genetic drift change their genetic attributes. For mud carp (Henicorhynchus spp.), understanding the status of breeding units in the MRB will be important for their long term persistence, sustainability and for implementing effective management strategies. Earlier analysis of stock structure in two economically important mud carp species (Henicorhynchus siamensis and H. lobatus) in the MRB completed with mtDNA markers identified a number of populations of both species where gene flow had apparently been interrupted or reduced but applying these data directly to management unit identification is potentially compromised because information was only available about female dispersal patterns. The current study aimed to address this problem and to fully assess the extent of current gene flow (nDNA) and reproductive exchange among selected wild populations of two species of carp (Henicorhynchus spp.) of high economic importance in the MRB using combined mtDNA and nDNA markers. In combination, the data can be used to define effective management units for each species. In general, nDNA diversity for H. lobatus (with average allelic richness (A) 7.56 and average heterozygosity (Ho) 0.61) was very similar to that identified for H. siamensis (A = 6.81 and Ho = 0.75). Both mud carp species show significant but low FST estimates among populations as a result of lower genetic diversity among sampled populations compared with genetic diversity within populations that may potentially mask any 'real' population structure. Overall, population genetic structure patterns from mtDNA and nDNA in both Henicorhynchus species were largely congruent. Different population structures however, were identified for the two Henicorhynchus species across the same geographical area. Apparent co-similarity in morphology and co-distribution of these two relatively closely related species does not apparently imply parallel evolutionary histories. Differences in each species population structure likely reflect historical drainage rearrangement of the Mekong River. The data indicate that H. siamensis is likely to have occupied the Mekong system for much longer than has H. lobatus in the past. Two divergent stocks were identified for H. lobatus in the MRB below the Khone Falls while a single stock had been evident in the earlier mtDNA study. This suggests that the two Henicorhynchus species may possess different life history traits and that different patterns of gene flow has likely influenced modern genetic structure in these close congeners. In combination, results of the earlier mtDNA and the current study have implications for effective management of both Henicorhynchus species across the MRB. Currently, both species are essentially treated as a single management unit in this region. This strategy may be appropriate for H. lobatus as a single stock was evident in the main stream of the MRB, but may not be appropriate for H. siamensis as more than a single stock was identified across the same range for this species. Management strategies should consider this difference to conserve overall biodiversity (local discrete populations) and this will include maintaining natural habitat and migration pathways, provision of fish sanctuaries (refuges) and may also require close monitoring of any stock declines, a signal that may require effective recovery strategies.

Genetic variability in spotted seatrout (Cynoscion nebulosus), determined with microsatellite DNA markers

Variation in the allele frequencies of five microsatellite loci was surveyed in 1256 individual spotted seatrout (Cynoscion nebulosus) obtained from 12 bays and estuaries from Laguna Madre, Texas, to Charlotte Harbor, Florida, to St. John’s River on the Florida Atlantic Coast. Texas and Louisiana collection sites were resampled each year for two to four years (1998−2001). Genetic differentiation was observed. Spotted seatrout from Florida waters were strongly differentiated from spotted seatrout collected in Louisiana and Texas. The greatest genetic discontinuity was observed between Tampa Bay and Charlotte Harbor, and Charlotte Harbor seatrout were most similar to Atlantic Coast spotted seatrout. Texas and Louisiana samples were not strongly structured within the northwestern Gulf of Mexico and there was little evidence of temporal differentiation within bays. These findings are contrary to those of earlier analyses with allozymes and mitochondrial DNA (mtDNA) where evidence of spatial differentiation was found for spotted seatrout resident on the Texas coast. The differences in genetic structure observed among these markers may reflect differences in response to selective pressure, or may be due to differences in underlying genetic processes.

Origin of rabbit (Oryctolagus cuniculus) in China: evidence from mitochondrial DNA control region sequence analysis

A fragment of mitochondrial DNA (mtDNA) control region (similar to700 bp) was sequenced in 104 individuals from 20 breeds (three Chinese domestic breeds, five recently derived breeds and 12 introduced breeds) of domestic rabbits, Oryctolagus cuniculus . Nineteen sites were polymorphic, with 18 transitions and one insertion/deletion, and eight haplotypes (A1, A2, A3, A4, A5, A6, A7 and A8) were identified. Haplotype A1 was the most common and occurred in 89 individuals. In the 25 Chinese rabbits, only haplotype A1 was observed, while four haplotypes (A1, A3, A5 and A6) were found in 26 recently derived individuals. Haplotype A2 was shared by seven individuals among three introduced strains. The other six haplotypes accounted for 0. 96-1. 92% of the animals. Combined with the published sequences of European rabbits, a reduced median-joining network was constructed. The Chinese rabbit mtDNAs were scattered into two clusters of European rabbits. These results suggest that the (so-called) Chinese rabbits were introduced from Europe. Genetic diversity in Chinese rabbits was very low.

Identification, cloning, and characterization of microsatellite DNA in Euglena gracilis

Microsatellite DNA has been developed into one of the most popular genetic markers. We have identified and cloned microsatellite loci in the genome of a free-living protozoan Euglena gracilis FACHB-848, using the random amplified microsatellites method (RAMS). The digoxigenin-labelled oligonucleotides(CT)(10) and (GT)(10) served as probes to detect complementary sequences in the randomly amplified polymorphic DNA (RAPD) fingerprints produced by means of Southern blotting. Subsequently, positive RAPD fragments were cloned. From a total of 31 RAPD primer profiles, eight microsatellite loci of E. gracilis were detected and characterized. Further, six sites (i.e. EGMS1, EGMS3, EGMS4, EGMS5, EGMS6, and EGMS7) showed polymorphisms. We found a GT or CT microsatellite every 10.5 kb in the genome of E. gracilis, and similar to animal genomes, the (GT)(n) motif was much more abundant than the (CT)(n) motif. These polymorphic microsatellite DNA will serve as advantageous molecular markers for studying the genetic diversity and molecular ecology of Euglena.

Genetic maternity and paternity in a local population of armadillos assessed by microsatellite DNA markers and field data.

Genetic data from polymorphic microsatellite loci were employed to estimate paternity and maternity in a local population of nine-banded armadillos (Dasypus novemcinctus) in northern Florida. The parentage assessments took advantage of maximum likelihood procedures developed expressly for situations when individuals of neither gender can be excluded a priori as candidate parents. The molecular data for 290 individuals, interpreted alone and in conjunction with detailed biological and spatial information for the population, demonstrate high exclusion probabilities and reasonably strong likelihoods of genetic parentage assignment in many cases; low mean probabilities of successful reproductive contribution to the local population by individual armadillo adults in a given year; and statistically significant microspatial associations of parents and their offspring. Results suggest that molecular assays of highly polymorphic genetic systems can add considerable power to assessments of biological parentage in natural populations even when neither parent is otherwise known.

Lack of genetic divergence found with microsatellite DNA markers in the tarakihi Nemadactylus macropterus

Three classes of molecular markers are commonly employed during population genetic studies of marine taxa: allozymes, mitochondrial DNA (mtDNA), and microsatellite DNA. These markers differ in their levels of polymorphism, and the ease and cost of their application. Nemadactylus macropterus is a commercially important marine fish from New Zealand and southern Australia that has been the subject of genetic (allozyme, mtDNA) and non-genetic (otolith microchemistry, larval advection) studies of stock structure. We collected microsatellite DNA data from this species to compare the utility of these molecular markers with those genetic methods previously applied to N. macropterus. Microsatellites did not indicate significant divergence among Australian samples, or between Australian and New Zealand samples. The latter is incongruent with the allozyme and mtDNA studies, and it is suggested that allelic homoplasy has hindered the resolution of population structure when using microsatellites.

Microsatellite DNA analysis of southeast Australian Haliotis laevigata (Donovan) populations - implications for ranching in Port Phillip Bay

The genetic composition of greenlip abalone (Haliotis laevigata) from Point Cook in Port Phillip Bay was examined prior to the aggregation of individuals from this site for ranching. The very thinly distributed natural population at Point Cook was believed to be of low genetic diversity, because the animals all originated from a single spawning event 5 y previously. Animals from Point Cook were compared with other H. laevigata from two sampling sites within Port Phillip Bay, and two sites outside the Bay in Bass Strait, to examine their genetic diversity and origin. Variation was assessed at five microsatellite loci. Deviations from Hardy-Weinberg equilibrium (HWE) were observed at some loci in various populations, but the Point Cook population was in HWE at all five loci. Mean heterozygosity and number of alleles was similar in all populations. Hierarchical analysis of molecular variance indicated significant genetic variation among populations, but did not differentiate Port Phillip Bay from Bass Strait populations. Pairwise comparisons of multilocus F_STand R_ST indicated significant genetic differences between Point Cook and some populations, as well as between other populations, but no consistent spatial pattern of differentiation was observed. There was no significant correlation between genetic and geographic distance. The level of genetic variation observed in the Point Cook individuals was similar to that in individuals from the other four sites, and sufficient to support a ranching program. However, this variation should be monitored to maximize genetic potential, and avoid commercially undesirable effects of inbreeding. Implications of this study in relation to the management of a ranching population in Port Phillip Bay are discussed.

Population structure in the highly fragmented range of Tor douronensis (Cyprinidae) in Sarawak, Malaysia revealed by microsatellite DNA markers

1. Understanding the structure of threatened populations, particularly those that exist in degraded or fragmented habitats is crucial for their effective management and conservation. Recently developed methods of individual-based analysis of genetic data provide an unprecedented opportunity to understand the relationships amongst fragmented populations.

2. In the present study, population structure of an important cyprinid species (Tor douronensis), which is indigenous to Sarawak, Malaysia, is investigated as part of an ongoing conservation effort to restore threatened wild populations of the species. The population structure inferred using data from seven autosomal microsatellite loci was generally consistent with geography and habitat fragmentation.

3. The results indicate that there are two well-defined clusters of T. douronensis in Sarawak, namely the 'northeastern' and the 'southwestern' clusters. In addition, a further subdivision was observed in each of the clusters distributed between river systems. Low levels of gene flow were also observed and migrants between habitat fragments were identified, possibly resulting from human-mediated translocations.

4. Implications of the findings for management and conservation of T. douronensis are discussed.

THEORETICAL STUDIES ON THE METHODS OF RECONSTRUCTING PHYLOGENETIC TREES FROM DNA SEQUENCE DATA (MOLECULAR EVOLUTION, HOMINOID EVOLUTION, COMPUTER SIMULATION)

(1) A mathematical theory for computing the probabilities of various nucleotide configurations is developed, and the probability of obtaining the correct phylogenetic tree (model tree) from sequence data is evaluated for six phylogenetic tree-making methods (UPGMA, distance Wagner method, transformed distance method, Fitch-Margoliash's method, maximum parsimony method, and compatibility method). The number of nucleotides (m*) necessary to obtain the correct tree with a probability of 95% is estimated with special reference to the human, chimpanzee, and gorilla divergence. m* is at least 4,200, but the availability of outgroup species greatly reduces m* for all methods except UPGMA. m* increases if transitions occur more frequently than transversions as in the case of mitochondrial DNA. (2) A new tree-making method called the neighbor-joining method is proposed. This method is applicable either for distance data or character state data. Computer simulation has shown that the neighbor-joining method is generally better than UPGMA, Farris' method, Li's method, and modified Farris method on recovering the true topology when distance data are used. A related method, the simultaneous partitioning method, is also discussed. (3) The maximum likelihood (ML) method for phylogeny reconstruction under the assumption of both constant and varying evolutionary rates is studied, and a new algorithm for obtaining the ML tree is presented. This method gives a tree similar to that obtained by UPGMA when constant evolutionary rate is assumed, whereas it gives a tree similar to that obtained by the maximum parsimony tree and the neighbor-joining method when varying evolutionary rate is assumed. ^

Sedimentological and palynological study of sediment cores from the northern Mekong River Delta

A sedimentological and palynological study of three sediment cores from the northern Mekong River Delta shows the regional sedimentary and environmental development since the mid-Holocene sea level highstand. A sub- to intertidal flat deposit of mid-Holocene age is recorded in the northernmost core. Shoreline deposits in all three cores show descending ages from N to S documenting 1) the early stages of the late Holocene regression and 2) the subsequent delta progradation. The delta plain successions vary from floodplain deposits with swamp-like elements to natural levee sediments. The uppermost sediments in all cores show human disturbance to varying degrees. The most intense alteration is recorded in the northernmost core where the palynological signal together with a charcoal peak indicates the profound change of the environment during the modern land reclamation. The sediments from at least one of the three presented cores do not show a "true" delta facies succession, but rather estuary-like features, as also observed in records from southern Cambodia. This absence is probably due to lack of accommodation space during the initial phase of rapid delta progradation which impeded the development of "true" delta successions as shown in cores from the southern Mekong River Delta.

Sedimentological and palynological analysis of sediment cores from the northeastern Vietnamese Mekong River Delta

Three radiocarbon-dated sediment cores from the northeastern Vietnamese Mekong River Delta have been analysed with a multiproxy approach (grain size, pollen and spores, macro-charcoal, carbon content) to unravel the palaeoenvironmental history of the region since the mid Holocene. During the mid-Holocene sea-level highstand a diverse, zoned and widespread mangrove belt (dominated by Rhizophora) covered the extended tidal flats. The subsequent regression and coeval delta progradation led to the rapid development of a back-mangrove community dominated by Ceriops and Bruguiera but also represented locally by e.g. Kandelia, Excoecaria and Phoenix. Along rivers this community seems to have endured even when the adjoining floodplain had already shifted to freshwater vegetation. Generally this freshwater vegetation has a strong swamp signature but locally Arecaceae, Fabaceae, Moraceae/Urticaceae and Myrsinaceae are important and mirror the geomorphological diversity of the delta plain. The macro-charcoal record implies that natural burning of vegetation occurred throughout the records, however, the occurrence of the highest amounts of macro-charcoal particles is linked with modern human activity.