Narongrit Muangmai

Evolutionary Dynamics of the Gametologous CTNNB1 Gene on the Z and W Chromosomes of Snakes

Snakes exhibit genotypic sex determination with female heterogamety (ZZ males and ZW females), and the state of sex chromosome differentiation also varies among lineages. To investigate the evolutionary history of homologous genes located in the nonrecombining region of differentiated sex chromosomes in snakes, partial sequences of the gametologous CTNNB1 gene were analyzed for 12 species belonging to henophid (Cylindrophiidae, Xenopeltidae, and Pythonidae) and caenophid snakes (Viperidae, Elapidae, and Colubridae). Nonsynonymous/synonymous substitution ratios (Ka/Ks) in coding sequences were low (Ka/Ks < 1) between CTNNB1Z and CTNNB1W, suggesting that these 2 genes may have similar functional properties. However, frequencies of intron sequence substitutions and insertion–deletions were higher in CTNNB1Z than CTNNB1W, suggesting that Z-linked sequences evolved faster than W-linked sequences. Molecular phylogeny based on both intron and exon sequences showed the presence of 2 major clades: 1) Z-linked sequences of Caenophidia and 2) W-linked sequences of Caenophidia clustered with Z-linked sequences of Henophidia, which suggests that the sequence divergence between CTNNB1Z and CTNNB1W in Caenophidia may have occurred by the cessation of recombination after the split from Henophidia.

High genetic diversity and demographic history of captive Siamese and Saltwater crocodiles suggest the first step toward the establishment of a breeding and reintroduction program in Thailand

The Siamese crocodile (Crocodylus siamensis) and Saltwater crocodile (C. porosus) are two of the most endangered animals in Thailand. Their numbers have been reduced severely by hunting and habitat fragmentation. A reintroduction plan involving captive-bred populations that are used commercially is important and necessary as a conservation strategy to aid in the recovery of wild populations. Here, the genetic diversity and population structure of 69 individual crocodiles, mostly members of captive populations, were analyzed using both mitochondrial D-loop DNA and microsatellite markers. The overall haplotype diversity was 0.924–0.971 and the mean expected heterozygosity across 22 microsatellite loci was 0.578–0.701 for the two species. This agreed with the star-like shaped topology of the haplotype network, which suggests a high level of genetic diversity. The mean ratio of the number of alleles to the allelic range (M ratio) for the populations of both species was considerably lower than the threshold of 0.68, which was interpreted as indicative of a historical genetic bottleneck. Microsatellite markers provided evidence of introgression for three individual crocodiles, which suggest that hybridization might have occurred between C. siamensis and C. porosus. D-loop sequence analysis detected bi-directional hybridization between male and female individuals of the parent species. Therefore, identification of genetically non-hybrid and hybrid individuals is important for long-term conservation management. Relatedness values were low within the captive populations, which supported their genetic integrity and the viability of a breeding and reintroduction management plan. This work constitutes the first step in establishing an appropriate source population from a scientifically managed perspective for an in situ/ex situ conservation program and reintroduction of crocodile individuals to the wild in Thailand.

Role of Chromosome Changes in Crocodylus Evolution and Diversity

The karyotypes of most species of crocodilians were studied using conventional and molecular cytogenetics. These provided an important contribution of chromosomal rearrangements for the evolutionary processes of Crocodylia and Sauropsida (birds and reptiles). The karyotypic features of crocodilians contain small diploid chromosome numbers (30~42), with little interspecific variation of the chromosome arm number (fundamental number) among crocodiles (56~60). This suggested that centric fusion and/or fission events occurred in the lineage, leading to crocodilian evolution and diversity. The chromosome numbers of Alligator, Caiman, Melanosuchus, Paleosuchus, Gavialis, Tomistoma, Mecistops, and Osteolaemus were stable within each genus, whereas those of Crocodylus (crocodylians) varied within the taxa. This agreed with molecular phylogeny that suggested a highly recent radiation of Crocodylus species. Karyotype analysis also suggests the direction of molecular phylogenetic placement among Crocodylus species and their migration from the Indo-Pacific to Africa and The New World. Crocodylus species originated from an ancestor in the Indo-Pacific around 9~16 million years ago (MYA) in the mid-Miocene, with a rapid radiation and dispersion into Africa 8~12 MYA. This was followed by a trans-Atlantic dispersion to the New World between 4~8 MYA in the Pliocene. The chromosomes provided a better understanding of crocodilian evolution and diversity, which will be useful for further study of the genome evolution in Crocodylia.

Assessment of snake DNA barcodes based on mitochondrial COI and Cytb genes revealed multiple putative cryptic species in Thailand

DNA barcodes of mitochondrial cytochrome c oxidase I (COI), cytochrome b (Cytb) genes, and their combined data sets were constructed from 35 snake species in Thailand. No barcoding gap was detected in either of the two genes from the observed intra- and interspecific sequence divergences. Intra- and interspecific sequence divergences of the COI gene differed 14 times, with barcode cut-off scores ranging over 2%-4% for threshold values differentiated among most of the different species; the Cytb gene differed 6 times with cut-off scores ranging over 2%-6%. Thirty-five specific nucleotide mutations were also found at interspecific level in the COI gene, identifying 18 snake species, but no specific nucleotide mutation was observed for Cytb in any single species. This suggests that COI barcoding was a better marker than Cytb. Phylogenetic clustering analysis indicated that most species were represented by monophyletic clusters, suggesting that these snake species could be clearly differentiated using COI barcodes. However, the two-marker combination of both COI and Cytb was more effective, differentiating snake species by over 2%-4%, and reducing species numbers in the overlap value between intra- and interspecific divergences. Three species delimitation algorithms (general mixed Yule-coalescent, automatic barcoding gap detection, and statistical parsimony network analysis) were extensively applied to a wide range of snakes based on both barcodes. This revealed cryptic diversity for eleven snake species in Thailand. In addition, eleven accessions from the database previously grouped under the same species were represented at different species level, suggesting either high genetic diversity, or the misidentification of these sequences in the database as a consequence of cryptic species.

A novel method of caenophidian snake sex identification using molecular markers based on two gametologous genes

Abstract Sex identification provides important information for ecological and evolutionary studies, as well as benefiting snake conservation management. Traditional methods such as cloacal probing or cloacal popping are counterproductive for sex identification concerning very small species, resulting in difficulties in the management of their breeding programs. In this study, the nucleotide sequences of gametologous genes (CTNNB1 and WAC genes) were used for the development of molecular sexing markers in caenophidian snakes. Two candidate markers were developed with the two primer sets, and successfully amplified by a single band on the agarose gel in male (ZZ) and two bands, differing in fragment sizes, in female (ZW) of 16 caenophidian snakes for CTNNB1 and 12 caenophidian snakes for WAC. Another candidate marker was developed with the primer set to amplify the specific sequence for CTNNB1W homolog, and the PCR products were successfully obtained in a female‐specific 250‐bp DNA bands. The three candidate PCR sexing markers provide a simple sex identification method based on the amplification of gametologous genes, and they can be used to facilitate effective caenophidian snake conservation and management programs.

Lack of satellite DNA species-specific homogenization and relationship to chromosomal rearrangements in monitor lizards (Varanidae, Squamata)

BackgroundSatellite DNAs (stDNAs) are highly repeated sequences that constitute large portions of any genome. The evolutionary dynamics of stDNA (e.g. copy number, nucleotide sequence, location) can, therefore, provide an insight into genome organization and evolution. We investigated the evolutionary origin of VSAREP stDNA in 17 monitor lizards (seven Asian, five Australian, and five African) at molecular and cytogenetic level.ResultsResults revealed that VSAREP is conserved in the genome of Asian and Australian varanids, but not in African varanids, suggesting that these sequences are either differentiated or lost in the African varanids. Phylogenetic and arrangement network analyses revealed the existence of at least four VSAREP subfamilies. The similarity of each sequence unit within the same VSAREP subfamily from different species was higher than those of other VSAREP subfamilies belonging to the same species. Additionally, all VSAREP subfamilies isolated from the three Australian species (Varanus rosenbergi, V. gouldii, and V. acanthurus) were co-localized near the centromeric or pericentromeric regions of the macrochromosomes, except for chromosomes 3 and 4 in each Australian varanid. However, their chromosomal arrangements were different among species.ConclusionsThe VSAREP stDNA family lack homogenized species-specific nucleotide positions in varanid lineage. Most VSAREP sequences were shared among varanids within the four VSAREP subfamilies. This suggests that nucleotide substitutions in each varanid species accumulated more slowly than homogenization rates in each VSAREP subfamily, resulting in non-species-specific evolution of stDNA profiles. Moreover, changes in location of VSAREP stDNA in each Australian varanid suggests a correlation with chromosomal rearrangements, leading to karyotypic differences among these species.

Complete mitochondrial genome of mouthbrooding fighting fish (Betta pi) compared with bubble nesting fighting fish (B. splendens)

Abstract Betta pi is the largest species of mouthbrooding fighting fish, while B. splendens is a globally ornamental bubble nesting fish. Complete mitochondrial genomes (mitogenomes) of wild individuals of B. pi and B. splendens were determined. The mitogenome sequences were 16,521 and 16,980 base pair in length, containing 37 genes with gene order identical to most teleost mitogenomes. Overall A + T content was 57.72% for B. pi and 61.92% for B. splendens. Phylogenetic analysis showed that B. pi and B. splendens were highly supported monophyletic clades. Our results will facilitate further genetic studies, including mitochondrial variations and population structure of fighting fishes.