Search

Silver barb (Barbonymus gonionotus) is among the most economically important freshwater fish species in Thailand. It ranks fourth in economic value and third in production weight for fisheries and culture in Thailand. An XX/XY sex-determination system based on gynogenesis was previously reported for this fish. In this study, the molecular basis underlying the sex-determination system was further investigated. Genome-wide single-nucleotide polymorphism data were generated for 32 captive-bred silver barb individuals, previously scored by phenotypic sex, to identify sex-linked regions associated with sex determination. Sixty-three male-linked loci, indicating putative XY chromosomes, were identified. Male-specific loci were not observed, which indicates that the putative Y chromosome is young and the sex determination region is cryptic. A homology search revealed that most male-linked loci were homologous to the Mariner/Tc1 and Gypsy transposable elements and are probably the remnants of an initial accumulation of repeats on the Y chromosome from the early stages of sex chromosome differentiation. This research provides convincing insights into the mechanism of sex determination and reveals the potential sex determination regions in silver barb. The study provides the basic data necessary for increasing the commercial value of silver barbs through genetic improvements.

The African catfish ( Clarias gariepinus ) may exhibit the co-existence of XX/XY and ZZ/ZW sex-determination systems (SDSs). However, the SDS of African catfish might be influenced by a polygenic sex-determination (PSD) system, comprising multiple independently segregating sex "switch" loci to determine sex within a species. Here, we aimed to detect the existence of PSD using hybrid. The hybrid produced by crossing male African catfish with female bighead catfish ( C. macrocephalus , XX/XY) is a good animal model to study SDSs. Determining the SDS of hybrid catfish can help in understanding the interactions between these two complex SDS systems. Using the genotyping-by-sequencing "DART-seq" approach, we detected seven moderately male-linked loci and seventeen female-linked loci across all the examined hybrid specimens. Most of these loci were not sex-linked in the parental species, suggesting that the hybrid exhibits a combination of different alleles. Annotation of the identified sex-linked loci revealed the presence of one female-linked locus homologous with the B4GALNT1 gene, which is involved in the spermatogenesis pathway and hatchability. However, this locus was not sex-linked in the parental species, and the African catfish might also exhibit PSD., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Nguyen, Ponjarat, Laopichienpong, Panthum, Singchat, Ahmad, Kraichak, Muangmai, Duengkae, Peyachoknagul, Na-Nakorn and Srikulnath.)

An investigation of sex-specific loci may provide important insights into fish sex determination strategies. This may be useful for biotechnological purposes, for example, to produce all-male or all-female fish for commercial breeding. The North African catfish species, Clarias gariepinus , has been widely adopted for aquaculture because its superior growth and disease resistance render the species suitable for hybridization with other catfish to improve the productivity and quality of fish meat. This species has either a ZZ/ZW or XX/XY sex determination system. Here, we investigate and characterize these systems using high-throughput genome complexity reduction sequencing as Diversity Arrays Technology. This approach was effective in identifying moderately sex-linked loci with both single-nucleotide polymorphisms (SNPs) and restriction fragment presence/absence (PA) markers in 30 perfectly sexed individuals of C. gariepinus . However, SNPs based markers were not found in this study. In total, 41 loci met the criteria for being moderately male-linked (with male vs. female ratios 80:20 and 70:30), while 25 loci were found to be moderately linked to female sex. No strictly male- or female-linked loci were detected. Seven moderately male-linked loci were partially homologous to some classes of transposable elements and three moderately male-linked loci were partially homologous to functional genes. Our data showed that the male heterogametic XX/XY sex determination system should co-exist with the ZZ/ZW system in C. gariepinus . Our finding of the co-existence of XX/XY and ZZ/ZW systems can be applied to benefit commercial breeding of this species in Thailand. This approach using moderately sex-linked loci provides a solid baseline for revealing sex determination mechanisms and identify potential sex determination regions in catfish, allowing further investigation of genetic improvements in breeding programs., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Nguyen, Panthum, Ponjarat, Laopichienpong, Kraichak, Singchat, Ahmad, Muangmai, Peyachoknagul, Na-Nakorn and Srikulnath.)

Elucidation of the process of sex chromosome differentiation is necessary to understand the dynamics of evolutionary mechanisms in organisms. The W sex chromosome of the Siamese cobra (Naja kaouthia) contains a large number of repeats and shares amniote sex chromosomal linkages. Diversity Arrays Technology provides an effective approach to identify sex-specific loci that are epoch-making, to understand the dynamics of molecular transitions between the Z and W sex chromosomes in a snake lineage. From a total of 543 sex-specific loci, 90 showed partial homology with sex chromosomes of several amniotes and 89 loci were homologous to transposable elements. Two loci were confirmed as W-specific nucleotides after PCR amplification. These loci might result from a sex chromosome differentiation process and involve putative sex-determination regions in the Siamese cobra. Sex-specific loci shared linkage homologies among amniote sex chromosomes, supporting an ancestral super-sex chromosome., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Transposable elements (TEs) are dynamic elements present in all eukaryotic genomes. They can "jump" and amplify within the genome and promote segmental genome rearrangements on both autosomes and sex chromosomes by disruption of gene structures. The Bovine-B long interspersed nuclear element (Bov-B LINE) is among the most abundant TE-retrotransposon families in vertebrates due to horizontal transfer (HT) among vertebrate lineages. Recent studies have shown multiple HTs or the presence of diverse Bov-B LINE groups in the snake lineage. It is hypothesized that Bov-B LINEs are highly dynamic and that the diversity reflects multiple HTs in snake lineages. Partial sequences of Bov-B LINE from 23 snake species were characterized. Phylogenetic analysis resolved at least two Bov-B LINE groups that might correspond to henophidian and caenophidian snakes; however, the tree topology differed from that based on functional nuclear and mitochondrial gene sequences. Several Bov-B LINEs of snakes showed greater than 80% similarity to sequences obtained from insects, whereas the two Bov-B LINE groups as well as sequences from the same snake species classified in different Bov-B LINE groups showed sequence similarities of less than 80%. Calculation of estimated divergence time and pairwise divergence between all individual Bov-B LINE copies suggest invasion times ranging from 79.19 to 98.8 million years ago in snakes. Accumulation of elements in a lineage-specific fashion ranged from 9 × 10 -6 % to 5.63 × 10 -2 % per genome. The genomic proportion of Bov-B LINEs varied among snake species but was not directly associated with genome size or invasion time. No differentiation in Bov-B LINE copy number between males and females was observed in any of the snake species examined. Incongruence in tree topology between Bov-B LINEs and other snake phylogenies may reflect past HT events. Sequence divergence of Bov-B LINEs between copies suggests that recent multiple HTs occurred within the same evolutionary timeframe in the snake lineage. The proportion of Bov-B LINEs varies among species, reflecting species specificity in TE invasion. The rapid speciation of snakes, coinciding with Bov-B LINE invasion in snake genomes, leads us to better understand the effect of Bov-B LINEs on snake genome evolution., Competing Interests: The authors declare no conflicts of interest.

The majority of lizards classified in the superfamily Iguanoidea have an XX/XY sex-determination system in which sex-chromosomal linkage shows homology with chicken ( Gallus gallus ) chromosome 15 (GGA15). However, the genomics of sex chromosomes remain largely unexplored owing to the presence of homomorphic sex chromosomes in majority of the species. Recent advances in high-throughput genome complexity reduction sequencing provide an effective approach to the identification of sex-specific loci with both single-nucleotide polymorphisms (SNPs) and restriction fragment presence/absence (PA), and a better understanding of sex chromosome dynamics in Iguanoidea. In this study, we applied Diversity Arrays Technology (DArTseq TM ) in 29 phenotypic sex assignments (14 males and 15 females) of green iguana ( Iguana iguana ). We confirmed a male heterogametic (XX/XY) sex determination mode in this species, identifying 29 perfectly sex-linked SNP/PA loci and 164 moderately sex-linked SNP/PA loci, providing evidence probably indicative of XY recombination. Three loci from among the perfectly sex-linked SNP/PA loci showed partial homology with several amniote sex chromosomal linkages. The results support the hypothesis of an ancestral super-sex chromosome with overlaps of partial sex-chromosomal linkages. However, only one locus among the moderately sex-linked loci showed homology with GGA15, which suggests that the specific region homologous to GGA15 was located outside the non-recombination region but in close proximity to this region of the sex chromosome in green iguana. Therefore, the location of GGA15 might be further from the putative sex-determination locus in green iguana. This is a paradigm shift in understanding linkages on homomorphic X and Y sex chromosomes. The DArTseq platform provides an easy-to-use strategy for future research on the evolution of sex chromosomes in Iguanoidea, particularly for non-model species with homomorphic or highly cryptic sex chromosomes., (Copyright © 2020 Koomgun, Laopichienpong, Singchat, Panthum, Phatcharakullawarawat, Kraichak, Sillapaprayoon, Ahmad, Muangmai, Peyachoknagul, Duengkae, Ezaz and Srikulnath.)

Centromeric satellite DNA (cen-satDNA) sequences of the Asian swamp eel (Monopterus albus) were characterized. Three GC-rich cen-satDNA sequences were detected as a 233 bp MALREP-A and a 293 bp MALREP-B localized to all chromosomes, and a 293 bp MALREP-C distributed on eight chromosome pairs. Sequence lengths of MALREP-B and MALREP-C were 60 bp larger than that of MALREP-A, showing partial homology with core sequences (233 bp). Size differences between MALREP-A and MALREP-B/C suggest the possible occurrence of two satDNA families. The presence of an additional 60 bp in MALREP-B/C resulted from an ancient dimer of 233 bp monomers and subsequent mutation and homogenization between the two monomers. All MALREPs showed partial homology with transposable elements (TEs), suggesting that the MALREPs originated from the TEs. The MALREPs might have been acquired in the Asian swamp eel, thereby promoting fixation in the species., Competing Interests: Declaration of Competing Interest The authors declare no conflict of interest., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Squamate reptile chromosome 2 (SR2) is thought to be an important remnant of an ancestral amniote super-sex chromosome, but a recent study showed that the Siamese cobra W sex chromosome is also a part of this larger ancestral chromosome. To confirm the existence of an ancestral amniote super-sex chromosome and understand the mechanisms of amniote sex chromosome evolution, chromosome maps of two snake species [Russell's viper: Daboia russelii (DRU) and the common tiger snake: Notechis scutatus (NSC)] were constructed using bacterial artificial chromosomes (BACs) derived from chicken and zebra finch libraries containing amniote sex chromosomal linkages. Sixteen BACs were mapped on the W sex chromosome of DRU and/or NSC, suggesting that these BACs contained a common genomic region shared with the W sex chromosome of these snakes. Two of the sixteen BACs were co-localized to DRU2 and NSC2, corresponding to SR2. Prediction of genomic content from all BACs mapped on snake W sex chromosomes revealed a large proportion of long interspersed nuclear element (LINE) and short interspersed nuclear element (SINE) retrotransposons. These results led us to predict that amplification of LINE and SINE may have occurred on snake W chromosomes during evolution. Genome compartmentalization, such as transposon amplification, might be the key factor influencing chromosome structure and differentiation. Multiple sequence alignments of all BACs mapped on snake W sex chromosomes did not reveal common sequences. Our findings indicate that the SR2 and snake W sex chromosomes may have been part of a larger ancestral amniote super-sex chromosome, and support the view of sex chromosome evolution as a colorful myriad of situations and trajectories in which many diverse processes are in action., (Copyright © 2020 Singchat, Ahmad, Sillapaprayoon, Muangmai, Duengkae, Peyachoknagul, O’Connor, Griffin and Srikulnath.)

Captive breeding programs for endangered species can increase population numbers for eventual reintroduction to the wild. Captive populations are typically small and isolated, which results in inbreeding and reduction of genetic variability, and may lead to an increased risk of extinction. The Omkoi Wildlife Breeding Center maintains the only Thai captive Chinese goral (Naemorhedus griseus) population, and has plans to reintroduce individuals into natural isolated populations. Genetic variability was assessed within the captive population using microsatellite data. Although no bottleneck was observed, genetic variability was low (allelic richness = 7.091 ± 0.756, He = 0.455 ± 0.219; He < Ho) and 11 microsatellite loci were informative that likely reflect inbreeding. Estimates of small effective population size and limited numbers of founders, combined with wild-born individuals within subpopulations, tend to cause reduction of genetic variability over time in captive programs. This leads to low reproductive fitness and limited ability to adapt to environmental change, thereby increasing the risk of extinction. Management of captive populations as evolutionarily significant units with diverse genetic backgrounds offers an effective strategy for population recovery. Relocation of individuals among subpopulations, or introduction of newly captured wild individuals into the captive program will help to ensure the future security of Chinese goral. Implications for future conservation actions for the species are discussed herein., Competing Interests: The authors have declared that no competing interests exist.

Sex chromosomes in some amniotes share linkage homologies with distantly related taxa in regions orthologous to squamate reptile chromosome 2 (SR2) and the snake W sex chromosome. Thus, the SR2 and W chromosomes may formerly have been part of a larger ancestral amniote super-sex chromosome. Comparison of various sex chromosomal linkage homologies in Toxicofera with those in other amniotes offers an excellent model to assess key cytological differences, to understand the mechanisms of amniote sex chromosome evolution in each lineage and the existence of an ancestral amniote super-sex chromosome. Chromosome maps of four species of Toxicofera were constructed using bacterial artificial chromosomes (BACs) derived from chicken and zebra finch libraries containing amniote sex chromosomal linkages. Different macrochromosome linkage homologies were highly conserved among Toxicofera, and at least two BACs (CH261-125F1 and CH261-40D6) showed partial homology with sex chromosomes of amniotes associated with SR2, which supports the hypothesis of an ancestral super-sex chromosome with overlaps of partial linkage homologies. The present data also suggest a possible multiple fission mechanism of an ancestral super-sex chromosome, which resulted in further development of various sex chromosomal linkages of Toxicofera based on particular properties that favored the role of sex chromosomes.

To better understand PBI-DdeI satellite DNA located in the centromeric region of python, molecular evolution analysis was conducted on 40 snake species. A ladder-like pattern of DNA bands with repetition of the 194-210 bp monomer was observed in 15 species using PCR. Molecular cloning was performed to obtain 97 AT-rich monomer sequences. Phylogenetic and network analyses showed three PBI-DdeI subfamilies with sequences grouped in species-specific clusters, suggesting rapid evolution. Slow evolution was found in eight species with shared PBI-DdeI sequences, suggesting recent species diversification, allowing PBI-DdeI no time to diverge, with limited homogenization and fixation processes. Quantitative real-time PCR showed large differences in copy number between Python bivittatus and other snakes, consistent with repeat scanning of whole genome sequences. Copy numbers were significantly higher in female Naja kaouthia than in males, concurring with chromosomal distribution of PBI-DdeI specifically localized to female W chromosomes. PBI-DdeI might act as an evolutionary driver with several repeats to promote W chromosome differentiation and heterochromatinization in N. kaouthia. Analysis revealed PBI-DdeI with a reduced copy number, compared to P. bivittatus, in most snakes studied, and it is possible that it subsequently dispersed and amplified on W chromosomes with different functional roles in N. kaouthia.

The fragmentation of habitats and hunting have impacted the Asian woolly-necked stork (Ciconia episcopus), leading to a serious risk of extinction in Thailand. Programs of active captive breeding, together with careful genetic monitoring, can play an important role in facilitating the creation of source populations with genetic variability to aid the recovery of endangered species. Here, the genetic diversity and population structure of 86 Asian woolly-necked storks from three captive breeding programs [Khao Kheow Open Zoo (KKOZ) comprising 68 individuals, Nakhon Ratchasima Zoo (NRZ) comprising 16 individuals, and Dusit Zoo (DSZ) comprising 2 individuals] were analyzed using 13 microsatellite loci, to aid effective conservation management. Inbreeding and an extremely low effective population size (Ne) were found in the KKOZ population, suggesting that deleterious genetic issues had resulted from multiple generations held in captivity. By contrast, a recent demographic bottleneck was observed in the population at NRZ, where the ratio of Ne to abundance (N) was greater than 1. Clustering analysis also showed that one subdivision of the KKOZ population shared allelic variability with the NRZ population. This suggests that genetic drift, with a possible recent and mixed origin, occurred in the initial NRZ population, indicating historical transfer between captivities. These captive stork populations require improved genetic variability and a greater population size, which could be achieved by choosing low-related individuals for future transfers to increase the adaptive potential of reintroduced populations. Forward-in-time simulations such as those described herein constitute the first step in establishing an appropriate source population using a scientifically managed perspective for an in situ and ex situ conservation program in Thailand., Competing Interests: The authors have declared that no competing interests exist.

Telomeres comprise tandem repeated DNA sequences that protect the ends of chromosomes from deterioration or fusion with neighboring chromosomes, and their lengths might vary with sex and age. Here, age- and sex-related telomere lengths in male and female captive Siamese cobras ( Naja kaouthia ) were investigated using quantitative real-time polymerase chain reaction based on cross-sectional data. A negative correlation was shown between telomere length and body size in males but not in females. Age-related sex differences were also recorded. Juvenile female snakes have shorter telomeres relative to males at up to 5 years of age, while body size also rapidly increases during this period. This suggests that an accelerated increase in telomere length of female cobra results from sex hormone stimulation to telomerase activity, reflecting sexually dimorphic phenotypic traits. This might also result from amplification of telomeric repeats on sex chromosomes. By contrast, female Siamese cobras older than 5 years had longer telomeres than males. Diverse sex hormone levels and oxidative stress parameters between sexes may affect telomere length., Competing Interests: None declared.

Background: Unlike the chromosome constitution of most snakes (2n=36), the cobra karyotype shows a diploid chromosome number of 38 with a highly heterochromatic W chromosome and a large morphologically different chromosome 2. To investigate the process of sex chromosome differentiation and evolution between cobras, most snakes, and other amniotes, we constructed a chromosome map of the Siamese cobra (Naja kaouthia) with 43 bacterial artificial chromosomes (BACs) derived from the chicken and zebra finch libraries using the fluorescence in situ hybridization (FISH) technique, and compared it with those of the chicken, the zebra finch, and other amniotes., Results: We produced a detailed chromosome map of the Siamese cobra genome, focusing on chromosome 2 and sex chromosomes. Synteny of the Siamese cobra chromosome 2 (NKA2) and NKAZ were highly conserved among snakes and other squamate reptiles, except for intrachromosomal rearrangements occurring in NKA2. Interestingly, twelve BACs that had partial homology with sex chromosomes of several amniotes were mapped on the heterochromatic NKAW as hybridization signals such as repeat sequences. Sequence analysis showed that most of these BACs contained high proportions of transposable elements. In addition, hybridization signals of telomeric repeat (TTAGGG) n and six microsatellite repeat motifs ((AAGG) 8 , (AGAT) 8 , (AAAC) 8 , (ACAG) 8 , (AATC) 8 , and (AAAAT) 6 ) were observed on NKAW, and most of these were also found on other amniote sex chromosomes., Conclusions: The frequent amplification of repeats might involve heterochromatinization and promote sex chromosome differentiation in the Siamese cobra W sex chromosome. Repeat sequences are also shared among amniote sex chromosomes, which supports the hypothesis of an ancestral super-sex chromosome with overlaps of partial syntenies. Alternatively, amplification of microsatellite repeat motifs could have occurred independently in each lineage, representing convergent sex chromosomal differentiation among amniote sex chromosomes.

Mitochondrial genomes (mitogenomes) of five Cyrtodactylus were determined. Their compositions and structures were similar to most of the available gecko lizard mitogenomes as 13 protein-coding, two rRNA and 22 tRNA genes. The non-coding control region (CR) of almost all Cyrtodactylus mitogenome structures contained a repeated sequence named the 75-bp box family, except for C. auribalteatus which contained the 225-bp box. Sequence similarities indicated that the 225-bp box resulted from the duplication event of 75-bp boxes, followed by homogenization and fixation in C. auribalteatus . The 75-bp box family was found in most gecko lizards with high conservation (55-75% similarities) and could form secondary structures, suggesting that this repeated sequence family played an important role under selective pressure and might involve mitogenome replication and the likelihood of rearrangements in CR. The 75-bp box family was acquired in the common ancestral genome of the gecko lizard, evolving gradually through each lineage by independent nucleotide mutation. Comparison of gecko lizard mitogenomes revealed low structural diversity with at least six types of mitochondrial gene rearrangements. Cyrtodactylus mitogenome structure showed the same gene rearrangement as found in most gecko lizards. Advanced mitogenome information will enable a better understanding of structure evolution mechanisms., Competing Interests: The authors declare there are no competing interests.

Coconuts (Cocos nucifera L.) are divided by the height into tall and dwarf types. In many plants the short phenotype was emerged by mutation of the GA20ox gene encoding the enzyme involved in gibberellin (GA) biosynthesis. Two CnGA20ox genes, CnGA20ox1 and CnGA20ox2, were cloned from tall and dwarf types coconut. The sequences, gene structures and expressions were compared. The structure of each gene comprised three exons and two introns. The CnGA20ox1 and CnGA20ox2 genes consisted of the coding region of 1110 and 1131 bp, encoding proteins of 369 and 376 amino acids, respectively. Their amino acid sequences are highly homologous to GA20ox1 and GA20ox2 genes of Elaeis guineensis, but only 57% homologous to each other. However, the characteristic amino acids two histidines and one aspartic acid which are the two iron (Fe 2+ ) binding residues, and arginine and serine which are the substrate binding residues of the dioxygenase enzyme in the 20G-FeII&#95;Oxy domain involved in GA biosynthesis, were found in the active site of both enzymes. The evolutionary relationship of their proteins revealed three clusters in vascular plants, with two subgroups in dicots and three subgroups in monocots. This result confirmed that CnGA20ox was present as multi-copy genes, and at least two groups CnGA20ox1 and CnGA20ox2 were found in coconut. The nucleotide sequences of CnGA20ox1 gene in both coconut types were identical but its expression was about three folds higher in the leaves of tall coconut than in those of dwarf type which was in good agreement with their height. In contrast, the nucleotide sequences of CnGA20ox2 gene in the two coconut types were different, but the expression of CnGA20ox2 gene could not be detected in either coconut type. The promoter region of CnGA20ox1 gene was cloned, and the core promoter sequences and various cis-elements were found. The CnGA20ox1 gene should be responsible for the height in coconut, which is different from other plants because no mutation was present in CnGA20ox1 gene of dwarf type coconut.

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., Competing Interests: The authors report no conflict of interest and are responsible for the content and writing of this article. Animal care and all experimental procedures were approved by the Animal Experiment Committee, Kasetsart University, Thailand (approval no. ACKU01157), and conducted in accordance with the Regulations on Animal Experiments at Kasetsart University., (© 2017 The Author(s). Published by Informa UK Limited, trading as Taylor & Francis Group.)

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.

Background: Satellite 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., Results: Results 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., Conclusions: The 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.

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 CTNNB1 W 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.

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.

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., (Copyright © 2016 Elsevier B.V. All rights reserved.)

Genomic alteration in head and neck squamous cell carcinoma (HNSCC) was studied in two cell line pairs (HN30-HN31 and HN4-HN12) using conventional C-banding, multiplex fluorescence in situ hybridization (M-FISH), and array comparative genomic hybridization (array CGH). HN30 and HN4 were derived from primary lesions in the pharynx and base of tongue, respectively, and HN31 and HN12 were derived from lymph-node metastatic lesions belonging to the same patients. Gain of chromosome 1, 7, and 11 were shared in almost all cell lines. Hierarchical clustering revealed that HN31 was closely related to HN4, which shared eight chromosome alteration cases. Large C-positive heterochromatins were found in the centromeric region of chromosome 9 in HN31 and HN4, which suggests complex structural amplification of the repetitive sequence. Array CGH revealed amplification of 7p22.3p11.2, 8q11.23q12.1, and 14q32.33 in all cell lines involved with tumorigenesis and inflammation genes. The amplification of 2p21 (SIX3), 11p15.5 (H19), and 11q21q22.3 (MAML2, PGR, TRPC6, and MMP family) regions, and deletion of 9p23 (PTPRD) and 16q23.1 (WWOX) regions were identified in HN31 and HN12. Interestingly, partial loss of PTPRD (9p23) and WWOX (16q23.1) genes was identified in HN31 and HN12, and the level of gene expression tended to be the down-regulation of PTPRD, with no detectable expression of the WWOX gene. This suggests that the scarcity of PTPRD and WWOX genes might have played an important role in progression of HNSCC, and could be considered as a target for cancer therapy or a biomarker in molecular pathology.

We analyzed the structure and the expression of Kunitz chymotrypsin inhibitor (WCI) genes in winged bean. WCI was encoded by a multigene family which comprised at least seven members. From their primary structures, four genes (WCI-2, WCI-3a, WCI-3b, and WCI-x) were expected to be functional ones and the other three (WCI-P1, WCI-P2, and WCI-P3) to be pseudogenes. The nucleotide sequences of the WCI-3a and WCI-3b genes were nearly identical, and they encoded the WCI-3 protein, the major chymotrypsin inhibitor in seeds. The WCI-2 gene also encoded the chymotrypsin inhibitor found in seeds and the WCI-x gene was expected to encode an unidentified chymotrypsin inhibitor. WCI messenger RNA and protein accumulated mainly in developing seeds and tuberous roots, small amounts of WCI mRNA being present in stems and pericarps. In seeds, transient accumulation of WCI mRNA was observed during the seed maturation period. These results suggest that the expression of WCI genes is regulated organ-specifically and developmentally in winged bean.

The accumulation of the Kunitz-type chymotrypsin inhibitor WCI-3 in winged bean seeds is controlled developmentally. In vitro translation experiments showed that the WCI-3 mRNA was present in 35- and 40-day-old immature seeds after flowering. The size of the in vitro translation product is about 2 000 Da larger than that of the mature WCI-3 protein. The WCI-3 cDNA clones were isolated from a λgtll cDNA library of 35-day-old immature seeds by immunoscreening. A nearly full-length cDNA clone was obtained containing an open reading frame of 207 amino acid residues. The deduced sequence of the 183 carboxy terminal amino acids coincides precisely with the amino acid sequence determined for purified WCI-3. The amino terminal extension of 24 residues has the characteristics of a signal peptide. Northern hybridization analysis of total poly(A)(+) RNA showed that the WCI-3 mRNA is approximately 900 nucleotides long and accumulates in 35- and 40-day-old but not in 30-day-old immature seeds.