Optimizing Bangkaew dog breed identification using DNA technology

BACKGROUND

The Bangkaew dog is an indigenous dog breed in the Phitsanulok province of Thailand. This breed is recognized by the Fédération Cynologique Internationale (FCI), a global canine organization. The unique traits of the Bangkaew breed lead to purebred selection for breeding, while only their traits and pedigree from parental history are recorded. Determination of the risk of inbreeding depression and the origin of unknown DNA profiles is essential due to the challenges in predicting puppy characteristics, which are crucial for breed management and conservation.

OBJECTIVE

This study aimed to emphasize that current allelic frequency data for the Bangkaew dog breed must be considered for precise individual identification.

METHODS

Approximately 82 Bangkaew dogs from various Thai localities were studied using 15 microsatellite markers for genotypic monitoring and individual identification. Maternal genetic inheritance was assessed via mtDNA D-loop analysis.

RESULTS

The results revealed high genetic diversity in the Bangkaew breed, indicating low potential for inbreeding. We also found that using a 15 loci microsatellite panel was effective for the identification of Bangkaew dogs. The optimized 10 loci microsatellite genotyping panel developed in this study presents improved identification testing efficiency, promoting both time- and cost-effectiveness.

CONCLUSION

Analysis of microsatellite DNA markers in Bangkaew dogs using an optimized panel of 10 loci selected from 15 loci effectively facilitated individual identification. This approach not only enhances time and cost efficiency, but also provides accurate allelic frequency estimates, which are crucial for the realistic evaluation of DNA evidence.

Genome-wide sequence divergence of satellite DNA could underlie meiotic failure in male hybrids of bighead catfish and North African catfish (Clarias, Clariidae)

Hybrid sterility, a hallmark of postzygotic isolation, arises from parental genome divergence disrupting meiosis. While chromosomal incompatibility is often implicated, the underlying mechanisms remain unclear. This study investigated meiotic behavior and genome-wide divergence in bighead catfish (C. macrocephalus), North African catfish (C. gariepinus), and their sterile male hybrids (important in aquaculture). Repetitive DNA analysis using bioinformatics and cytogenetics revealed significant divergence in satellite DNA (satDNA) families between parental species. Notably, one hybrid exhibited successful meiosis and spermatozoa production, suggesting potential variation in sterility expression. Our findings suggest that genome-wide satDNA divergence, rather than chromosome number differences, likely contributes to meiotic failure and male sterility in these catfish hybrids.

Terrestrial wildlife as indicators of microplastic pollution in western Thailand

Plastic pollution in terrestrial wildlife represents a new conservation challenge, with research in this area, especially within protected areas (PAs), being scant. This study documents the accumulation of microplastics (MPs) in terrestrial wildlife both inside and outside PAs in western Thailand. Carcasses of road-killed vertebrates in good condition, as well as live tadpoles, were collected to examine their exposure to plastic pollution. The digestive tracts of the vertebrate carcasses and the entire bodies of tadpoles were analyzed for MPs, which were identified if they measured over 50 µm. A total of 136 individuals from 48 vertebrate species were examined. The sample comprised snakes (44.12%), birds (11.03%), lizards (5.15%), tadpoles (32.25%), amphibians (5.88%), and mammals (1.47%). In total, 387 MPs were found in 44 species (91.67%), with an average occurrence of 3.25 ± 3.63 MPs per individual or 0.05 ± 0.08 MPs per gram of body weight. The quantities of MPs significantly varied among the animal groups, both in terms of number per individual (p < 0.05) and number per gram of body weight (p < 0.01). Furthermore, a significant difference in MP quantities was observed between specimens collected inside and outside PAs on an individual basis (p < 0.05), but not on a body weight basis (p = 0.07). Most MPs were fibers (77%), followed by fragments (22.22%), with only a minimal presence of film (0.52%) and foam (0.26%). Of all the MPs identified, 36.84% were confirmed as plastics or fibers made from natural materials, and 31.58% were plastics, including Polyethylene (PE), polyethylene terephthalate (PET), polypropylene (PP), Polyvinylidene chloride (PVDC), and polyester (PES). Additionally, fibers made of cotton, and those containing polyurethane (PU), rayon, PES, and combinations of rayon and PU, were identified. The quantities of MPs were significantly influenced by animal body weight, factors associated with human settlement/activity, and land use types. Our findings highlight the prevalence of plastic pollution in terrestrial vertebrates within Thai PAs. Further toxicological studies are required to establish plastic pollution standards. It is proposed that snakes, obtained from road kills, could serve as a non-invasive method for monitoring plastic pollution, thus acting as an indicator of the pollution threat to species within terrestrial ecosystems. There is an urgent need for the standardization of solid waste management at garbage dump sites in remote areas, especially within PAs. Conservation education focusing on MP occurrence, potential sources, and impacts could enhance awareness, thereby influencing changes in behaviors and attitudes toward plastic waste management at the household level.

Weak purifying selection in allelic diversity of the ADSL gene in indigenous and local chicken breeds and red junglefowl in Thailand

High levels of purine and uric acid, which are associated with health issues such as gout and cardiovascular disease, are found in the meat of fast-growing broiler chickens, which raises concerns about the quality of chicken meat and the health of the consumers who consume it. High genetic homogeneity and uniformity, particularly in genes involved in the synthesis of inosine monophosphate (IMP) and subsequent process of purine synthesis, which are associated with the meat quality, are exhibited in commercial broiler chickens owing to intensive inbreeding programs. Adenosine succinate lyase (ADSL) is a key enzyme involved in de novo purine biosynthetic pathway and its genetic polymorphisms affect IMP metabolism and purine content. In this study, we investigated the polymorphism of the ADSL gene in indigenous and local chicken breeds and red junglefowl in Thailand, using metabarcoding and genetic diversity analyses. Five alleles with 73 single nucleotide polymorphisms in exon 2, including missense and silent mutations, which may act on the synthesis efficiency of IMP and purine. Their protein structures revealed changes in amino acid composition that may affect ADSL enzyme activity. Weak purifying selection in these ADSL alleles was observed in the chicken population studied, implying that the variants have minor fitness impacts and a greater probability of fixation of beneficial mutations than strong purifying selection. A potential selective sweep was observed in Mae Hong Son chickens, whose purine content was lower than that in other breeds. This suggests a potential correlation between variations of the ADSL gene and reduced purine content and an impact of ADSL expression on the quality of chicken meat. However, further studies are required to validate its potential availability as a genetic marker for selecting useful traits that are beneficial to human health and well-being.

Questioning inbreeding: Could outbreeding affect productivity in the North African catfish in Thailand?

The North African catfish (Clarias gariepinus) is a significant species in aquaculture, which is crucial for ensuring food and nutrition security. Their high adaptability to diverse environments has led to an increase in the number of farms that are available for their production. However, long-term closed breeding adversely affects their reproductive performance, leading to a decrease in production efficiency. This is possibly caused by inbreeding depression. To investigate the root cause of this issue, the genetic diversity of captive North African catfish populations was assessed in this study. Microsatellite genotyping and mitochondrial DNA D-loop sequencing were applied to 136 catfish specimens, collected from three populations captured for breeding in Thailand. Interestingly, extremely low inbreeding coefficients were obtained within each population, and distinct genetic diversity was observed among the three populations, indicating that their genetic origins are markedly different. This suggests that outbreeding depression by genetic admixture among currently captured populations of different origins may account for the low productivity of the North African catfish in Thailand. Genetic improvement of the North African catfish populations is required by introducing new populations whose origins are clearly known. This strategy should be systematically integrated into breeding programs to establish an ideal founder stock for selective breeding.

Research Note: Possible influence of thermal selection on patterns of HSP70 and HSP90 gene polymorphisms in Thai indigenous and local chicken breeds and red junglefowls

The thermal stress caused by global climate change adversely affects the welfare, productivity, and reproductive performance of farm animals, including chickens, and causes substantial economic losses. However, the understanding of the genetic basis of the indigenous chicken adaptation to high ambient temperatures is limited. Hence, to reveal the genetic basis of thermal stress adaptation in chickens, this study investigated polymorphisms in the heat shock protein 70 (HSP70) and HSP90 genes, known mechanisms of cellular defense against thermal stress in indigenous and local chicken breeds and red junglefowls in Thailand. The result revealed seven alleles of the HSP70 gene. One allele exhibited a missense mutation, where an amino acid changed from Asn to His in the substrate-binding and peptide-binding domains, which is exclusive to the Lao Pa Koi chicken breed. Twenty new alleles with silent mutations in the HSP90 gene highlighted its greater complexity. Despite this diversity, distinct population structures were not found for either HSP70 or HSP90, which suggests incomplete impact on the domestication process and selection. The low genetic diversity, shown by the sharing of alleles between red junglefowls and Thai indigenous and local chicken breeds, aligns with the hypothesis that these alleles have undergone selection in tropical regions, such as Thailand. Selection signature analysis suggests the purifying selection of HSP70 for thermotolerance. This study provides valuable insights for enhancing the conservation of genetic resources with thermotolerant traits, which are essential for developing breeding programs to increase poultry production in the context of global climate change.

Disclosing the hidden nucleotide sequences: a journey into DNA barcoding of raptor species in public repositories

BACKGROUND

In nucleotide public repositories, studies discovered data errors which resulted in incorrect species identification of several accipitrid raptors considered for conservation. Mislabeling, particularly in cases of cryptic species complexes and closely related species, which were identified based on morphological characteristics, was discovered. Prioritizing accurate species labeling, morphological taxonomy, and voucher documentation is crucial to rectify spurious data.

OBJECTIVE

Our study aimed to identify an effective DNA barcoding tool that accurately reflects the efficiency status of barcodes in raptor species (Accipitridae).

METHODS

Barcode sequences, including 889 sequences from the mitochondrial cytochrome c oxidase I (COI) gene and 1052 sequences from cytochrome b (Cytb), from 150 raptor species within the Accipitridae family were analyzed.

RESULTS

The highest percentage of intraspecific nearest neighbors from the nearest neighbor test was 88.05% for COI and 95.00% for Cytb, suggesting that the Cytb gene is a more suitable marker for accurately identifying raptor species and can serve as a standard region for DNA barcoding. In both datasets, a positive barcoding gap representing the difference between inter-and intra-specific sequence divergences was observed. For COI and Cytb, the cut-off score sequence divergences for species identification were 4.00% and 3.00%, respectively.

CONCLUSION

Greater accuracy was demonstrated for the Cytb gene, making it the preferred primary DNA barcoding marker for raptors.

Genetic insights: mapping sex-specific loci in Siamese cobra (Naja kaouthia) sheds light on the putative sex determining region

The location of female-specific/linked loci identified in Siamese cobra (Naja kaouthia) previously has been determined through in silico chromosome mapping of the Indian cobra genome (N. naja) as a reference genome. In the present study, we used in silico chromosome mapping to identify sex-specific and linked loci in Siamese cobra. Many sex-specific and sex-linked loci were successfully mapped on the Z sex chromosome, with 227 of the 475 specific loci frequently mapped in a region covering 57 Mb and positioned at 38,992,675-95,561,177 bp of the Indian cobra genome (N. naja). This suggested the existence of a putative sex-determining region (SDR), with one specific locus (PA100000600) homologous to the TOPBP1 gene. The involvement of TOPBP1 gene may lead to abnormal synaptonemal complexes and meiotic chromosomal defects, resulting in male infertility. These findings offer valuable insights into the genetic basis and functional aspects of sex-specific traits in the Siamese cobra, which will contribute to our understanding of snake genetics and evolutionary biology.

Genome-wide SNP analysis provides insights into the XX/XY sex-determination system in silver barb (Barbonymus gonionotus)

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.

Hi-C sequencing unravels dynamic three-dimensional chromatin interactions in muntjac lineage: insights from chromosome fusions in Fea's muntjac genome

Eukaryotes have varying numbers and structures of characteristic chromosomes across lineages or species. The evolutionary trajectory of species may have been affected by spontaneous genome rearrangements. Chromosome fusion drastically alters karyotypes. However, the mechanisms and consequences of chromosome fusions, particularly in muntjac species, are poorly understood. Recent research-based advancements in three-dimensional (3D) genomics, particularly high-throughput chromatin conformation capture (Hi-C) sequencing, have allowed for the identification of chromosome fusions and provided mechanistic insights into three muntjac species: Muntiacus muntjak, M. reevesi, and M. crinifrons. This study aimed to uncover potential genome rearrangement patterns in the threatened species Fea's muntjac (Muntiacus feae), which have not been previously examined for such characteristics. Deep Hi-C sequencing (31.42 × coverage) was performed to reveal the 3D chromatin architecture of the Fea's muntjac genome. Patterns of repeated chromosome fusions that were potentially mediated by high-abundance transposable elements were identified. Comparative Hi-C maps demonstrated linkage homology between the sex chromosomes in Fea's muntjac and autosomes in M. reevesi, indicating that fusions may have played a crucial role in the evolution of the sex chromosomes of the lineage. The species-level dynamics of topologically associated domains (TADs) suggest that TAD organization could be altered by differential chromosome interactions owing to repeated chromosome fusions. However, research on the effect of TADs on muntjac genome evolution is insufficient. This study generated Hi-C data for the Fea's muntjac, providing a genomic resource for future investigations of the evolutionary patterns of chromatin conformation at the chromosomal level.

Preserving Pure Siamese Crocodile Populations: A Comprehensive Approach Using Multi-Genetic Tools

Hybrids between the critically endangered Siamese crocodile (Crocodylus siamensis) and least-concern saltwater crocodile (C. porosus) in captive populations represent a serious challenge for conservation and reintroduction programs due to the impact of anthropogenic activities. A previous study used microsatellite and mitochondrial DNA data to establish the criteria for identifying species and their hybrids; however, the results may have been influenced by biased allelic frequencies and genetic drift within the examined population. To overcome these limitations and identify the true signals of selection, alternative DNA markers and a diverse set of populations should be employed. Therefore, this study used DArT sequencing to identify genome-wide single nucleotide polymorphisms (SNPs) in both species and confirm the genetic scenario of the parental species and their hybrids. A population of saltwater crocodiles from Australia was used to compare the distribution of species-diagnostic SNPs. Different analytical approaches were compared to diagnose the level of hybridization when an admixture was present, wherein three individuals had potential backcrossing. Approximately 17.00-26.00% of loci were conserved between the Siamese and saltwater crocodile genomes. Species-diagnostic SNP loci for Siamese and saltwater crocodiles were identified as 8051 loci and 1288 loci, respectively. To validate the species-diagnostic SNP loci, a PCR-based approach was used by selecting 20 SNP loci for PCR primer design, among which 3 loci were successfully able to differentiate the actual species and different hybridization levels. Mitochondrial and nuclear genetic information, including microsatellite genotyping and species-diagnostic DNA markers, were combined as a novel method that can compensate for the limitations of each method. This method enables conservation prioritization before release into the wild, thereby ensuring sustainable genetic integrity for long-term species survival through reintroduction and management programs.

Genetic admixture and diversity in Thai domestic chickens revealed through analysis of Lao Pa Koi fighting cocks

Lao Pa Koi (LPK) chicken is a popular fighting breed in Thailand, prized for (its unique characteristics acquired by selective breeding), and a valuable model for exploring the genetic diversity and admixture of red junglefowls and domestic chickens. In this study, genetic structure and diversity of LPK chicken were assessed using 28 microsatellite markers and mitochondrial DNA (mtDNA) D-loop sequences, and the findings were compared to a gene pool library from "The Siam Chicken Bioresource Project". High genetic variability was observed in LPK chickens using mtDNA D-loop haplotype analysis, and six haplotypes were identified. Microsatellite data revealed 182 alleles, with an average of 6.5 alleles per locus. These results confirmed the occurrence of genetic admixture of red junglefowl and Thai domestic chickens in LPK chicken breed. A maximum entropy modeling approach was used to analyze the spatial suitability and to assess the adaptive evolution of LPK chickens in diverse local environments. The model identified 82.52% of the area studied as unsuitable, and 9.34%, 7.11%, and 2.02% of the area indicated moderate, low, and high suitability, respectively. The highest contribution rate to land suitability for LPK chickens was found at an elevation of 100-250 m, suggesting the importance of elevation for their potential distribution. The results of this study provide valuable insights into the genetic origin of LPK chicken breed and identify resources for future genetic improvement.

MicrosatNavigator: exploring nonrandom distribution and lineage-specificity of microsatellite repeat motifs on vertebrate sex chromosomes across 186 whole genomes

Microsatellites are short tandem DNA repeats, ubiquitous in genomes. They are believed to be under selection pressure, considering their high distribution and abundance beyond chance or random accumulation. However, limited analysis of microsatellites in single taxonomic groups makes it challenging to understand their evolutionary significance across taxonomic boundaries. Despite abundant genomic information, microsatellites have been studied in limited contexts and within a few species, warranting an unbiased examination of their genome-wide distribution in distinct versus closely related-clades. Large-scale comparisons have revealed relevant trends, especially in vertebrates. Here, "MicrosatNavigator", a new tool that allows quick and reliable investigation of perfect microsatellites in DNA sequences, was developed. This tool can identify microsatellites across the entire genome sequences. Using this tool, microsatellite repeat motifs were identified in the genome sequences of 186 vertebrates. A significant positive correlation was noted between the abundance, density, length, and GC bias of microsatellites and specific lineages. The (AC)n motif is the most prevalent in vertebrate genomes, showing distinct patterns in closely related species. Longer microsatellites were observed on sex chromosomes in birds and mammals but not on autosomes. Microsatellites on sex chromosomes of non-fish vertebrates have the lowest GC content, whereas high-GC microsatellites (≥ 50 M% GC) are preferred in bony and cartilaginous fishes. Thus, similar selective forces and mutational processes may constrain GC-rich microsatellites to different clades. These findings should facilitate investigations into the roles of microsatellites in sex chromosome differentiation and provide candidate microsatellites for functional analysis across the vertebrate evolutionary spectrum.

Optimizing Microsatellite Marker Panels for Genetic Diversity and Population Genetic Studies: An Ant Colony Algorithm Approach with Polymorphic Information Content

Microsatellites are polymorphic and cost-effective. Optimizing reduced microsatellite panels using heuristic algorithms eases budget constraints in genetic diversity and population genetic assessments. Microsatellite marker efficiency is strongly associated with its polymorphism and is quantified as the polymorphic information content (PIC). Nevertheless, marker selection cannot rely solely on PIC. In this study, the ant colony optimization (ACO) algorithm, a widely recognized optimization method, was adopted to create an enhanced selection scheme for refining microsatellite marker panels, called the PIC-ACO selection scheme. The algorithm was fine-tuned and validated using extensive datasets of chicken (Gallus gallus) and Chinese gorals (Naemorhedus griseus) from our previous studies. In contrast to basic optimization algorithms that stochastically initialize potential outputs, our selection algorithm utilizes the PIC values of markers to prime the ACO process. This increases the global solution discovery speed while reducing the likelihood of becoming trapped in local solutions. This process facilitated the acquisition of a cost-efficient and optimized microsatellite marker panel for studying genetic diversity and population genetic datasets. The established microsatellite efficiency metrics such as PIC, allele richness, and heterozygosity were correlated with the actual effectiveness of the microsatellite marker panel. This approach could substantially reduce budgetary barriers to population genetic assessments, breeding, and conservation programs.

Overcoming taxonomic challenges in DNA barcoding for improvement of identification and preservation of clariid catfish species

DNA barcoding without assessing reliability and validity causes taxonomic errors of species identification, which is responsible for disruptions of their conservation and aquaculture industry. Although DNA barcoding facilitates molecular identification and phylogenetic analysis of species, its availability in clariid catfish lineage remains uncertain. In this study, DNA barcoding was developed and validated for clariid catfish. 2,970 barcode sequences from mitochondrial cytochrome c oxidase I (COI) and cytochrome b (Cytb) genes and D-loop sequences were analyzed for 37 clariid catfish species. The highest intraspecific nearest neighbor distances were 85.47%, 98.03%, and 89.10% for COI, Cytb, and D-loop sequences, respectively. This suggests that the Cytb gene is the most appropriate for identifying clariid catfish and can serve as a standard region for DNA barcoding. A positive barcoding gap between interspecific and intraspecific sequence divergence was observed in the Cytb dataset but not in the COI and D-loop datasets. Intraspecific variation was typically less than 4.4%, whereas interspecific variation was generally more than 66.9%. However, a species complex was detected in walking catfish and significant intraspecific sequence divergence was observed in North African catfish. These findings suggest the need to focus on developing a DNA barcoding system for classifying clariid catfish properly and to validate its efficacy for a wider range of clariid catfish. With an enriched database of multiple sequences from a target species and its genus, species identification can be more accurate and biodiversity assessment of the species can be facilitated.

Environmental and Socio-Cultural Factors Impacting the Unique Gene Pool Pattern of Mae Hong-Son Chicken

Understanding the genetic diversity of domestic chicken breeds under the impact of socio-cultural and ecological dynamics is vital for the conservation of natural resources. Mae Hong Son chicken is a local breed of North Thai domestic chicken widely distributed in Mae Hong Son Province, Thailand; however, its genetic characterization, origin, and diversity remain poorly understood. Here, we studied the socio-cultural, environmental, and genetic aspects of the Mae Hong Son chicken breed and investigated its diversity and allelic gene pool. We genotyped 28 microsatellite markers and analyzed mitochondrial D-loop sequencing data to evaluate genetic diversity and assessed spatial habitat suitability using maximum entropy modeling. Sequence diversity analysis revealed a total of 188 genotyped alleles, with overall nucleotide diversity of 0.014 ± 0.007, indicating that the Mae Hong Son chicken population is genetically highly diverse, with 35 (M1-M35) haplotypes clustered into haplogroups A, B, E, and F, mostly in the North ecotype. Allelic gene pool patterns showed a unique DNA fingerprint of the Mae Hong Son chicken, as compared to other breeds and red junglefowl. A genetic introgression of some parts of the gene pool of red junglefowl and other indigenous breeds was identified in the Mae Hong Son chicken, supporting the hypothesis of the origin of the Mae Hong Son chicken. During domestication in the past 200-300 years after the crossing of indigenous chickens and red junglefowl, the Mae Hong Son chicken has adapted to the highland environment and played a significant socio-cultural role in the Northern Thai community. The unique genetic fingerprint of the Mae Hong Son chicken, retaining a high level of genetic variability that includes a dynamic demographic and domestication history, as well as a range of ecological factors, might reshape the adaptation of this breed under selective pressure.

The identity of Nanorana yunnanensis (Anderson, 1879) from Thailand, with a description of its tadpole (Anura, Dicroglossidae)

Nanorana yunnanensis has been reported in Thailand for many years; however, no specimens have been verifiably described. To verify the identity of this species in Thailand, we observed Thai specimens formerly assigned to the species N. bourreti and N. yunnanensis. The 12S rRNA, 16S rRNA, and ND2 mitochondrial genes were used to clarify the taxonomic statuses of the Thai specimens. The results showed that our specimens formed a clade with N. yunnanensis, containing specimens from southern Yunnan and Vietnam. The respective mean sequence divergences of the 12S and 16S rRNA genes between populations from Thailand versus China were 0.59% and 1.09%, for populations from Thailand versus Vietnam were 0.57% and 0.61%, and for populations from China versus Vietnam were 0.59% and 1.09%. We provide descriptions of the tadpole and adult specimens of N. yunnanensis from Thailand.

Should the Identification Guidelines for Siamese Crocodiles Be Revised? Differing Post-Occipital Scute Scale Numbers Show Phenotypic Variation Does Not Result from Hybridization with Saltwater Crocodiles

Populations of Siamese crocodiles (Crocodylus siamensis) have severely declined because of hunting and habitat fragmentation, necessitating a reintroduction plan involving commercial captive-bred populations. However, hybridization between Siamese and saltwater crocodiles (C. porosus) has occurred in captivity. Siamese crocodiles commonly have post-occipital scutes (P.O.) with 4–6 scales, but 2–6 P.O. scales were found in captives on Thai farms. Here, the genetic diversity and population structure of Siamese crocodiles with large P.O. variations and saltwater crocodiles were analyzed using mitochondrial DNA D-loop and microsatellite genotyping. Possible crocodile hybrids or phenotypic variations were ascertained by comparison with our previous library from the Siam Crocodile Bioresource Project. Siamese crocodiles with <4 P.O. scales in a row exhibit normal species-level phenotypic variation. This evidence encourages the revised description of Siamese crocodiles. Moreover, the STRUCTURE plot revealed large distinct gene pools, suggesting crocodiles in each farm were derived from distinct lineages. However, combining both genetic approaches provides evidence of introgression for several individual crocodiles, suggesting possible hybridization between Siamese and saltwater crocodiles. We proposed a schematic protocol with patterns observed in phenotypic and molecular data to screen hybrids. Identifying non-hybrid and hybrid individuals is important for long-term in situ/ex situ conservation.

Quality control of fighting fish nucleotide sequences in public repositories reveals a dark matter of systematic taxonomic implication

BACKGROUND

The number of nucleotide sequences in public repositories has exploded recently. However, the data contain errors, leading to incorrect species identification. Several fighting fish (Betta spp.) are poorly described, with unresolved cryptic species complexes masking undescribed species. Here, DNA barcoding was used to detect erroneous sequences in public repositories.

OBJECTIVE

This study reflects the current quantitative and qualitative status of DNA barcoding in fighting fish and provides a rapid and reliable identification tool.

METHODS

A total of 1034 barcode sequences were analyzed from mitochondrial cytochrome c oxidase I (COI) and cytochrome b (Cytb) genes from 71 fighting fish species.

RESULTS

The nearest neighbor test showed the highest percentage of intraspecific nearest neighbors at 93.41% for COI and 91.67% for Cytb, which can be used as reference barcodes for certain taxa. Intraspecific variation was usually less than 13%, while most species differed by more than 54%. The barcoding gap, calculated from the difference between inter- and intraspecific sequence divergences, was negative in the COI data set indicating overlapping intra- and interspecific sequence divergence. Sequence saturation was observed in the Cytb data set but not in the COI data set.

CONCLUSION

The COI gene should thus be used as the main barcoding marker for fighting fish.

Abundance and physical factors affecting the appearance of selected terrestrial birds in Khao Yai National Park using camera trapping

Abstract. Kanka P, Sukmasuang R, Duengkae P, Siripattaranugul K. 2023. Abundance and physical factors affecting the appearance of selected terrestrial birds in Khao Yai National Park using camera trapping. Biodiversitas 24: 222 232. This study used results from camera traps set for Dhole (Cuon alpinus) and their prey during 2017-2020 in Khao Yai National Park to assess the species richness of ground-dwelling birds and identify abundance and also temporal appearance of the species in the area. We selected only the Pheasant species that typically forage on the ground to investigate the physical factors affecting the appearance in Khao Yai National Park. Based on camera trapping data with a total of 4139 trap days, 115 camera locations showed 36 species of birds belonging to 21 families, and 13 orders of ground-dwelling bird species were recorded. Siamese Fireback (Lophura diardi), Red Junglefowl (Gallus gallus), and Silver Pheasant (Lophura nycthemera) were detected with an abundance index of 10.65%, 7.27%, and 1.84%, respectively. The ground-dwelling bird abundance index from the camera traps accounted for 26.85% of the total fauna. Although the Phasianidae were found to be diurnal species, this study also shows a statistically significant difference in the average time between Red Junglefowl and Siamese Fireback. Meanwhile, there was no statistically significant difference between Red Junglefowl and Silver Pheasant during their active time. The water sources showed a significant effect on the presence of the Siamese Fireback in data analysis; however, no physical environmental factors significantly affected the distribution of Red Junglefowl and Silver Pheasant. This study found that surrounding creeks and streams significantly impacted the appearance of the Siamese Fireback. As a result, area management needs to prevent and reduce disturbances along a creek in order to effectively manage the habitat of the Siamese Fireback and other wildlife.

Species diversity, abundance, and movement of small mammals in the dry evergreen forest at Khao Yai National Park, Thailand

Abstract. Chanachai Y, Nathalang A, Duengkae P, Sukmasuang R. 2022. Species diversity, abundance, and movement of small mammals in the dry evergreen forest at Khao Yai National Park, Thailand. Biodiversitas 23: 5892-5901. Small mammals play an integral role in the forest ecosystem. This research was conducted on The Mo Singto Forest Dynamics Plot, Khao Yai National Park. Trapping of small mammals was carried out in June-September 2019 and from January-April 2020 based on a total of 64 Tomahawk for small mammal live traps were systematically placed in an 8×8 grid of 20 m intervals, ripe bananas and palm oil as bait was used. Trap cages were opened for 3 consecutive months per season in each study period covering a total study period of 1 year. The result showed that a total of 4015 individuals from 9 species, 9 genera, 4 families and 3 orders were captured with a sampling effort of 32,555 trap nights. The most common species were Red Spiny Rat (Maxomys surifer), followed by Northern Tree Shrew (Tupaia belangeri), Long-tailed Giant Rat (Leopoldamys sabanus), Indochinese Ground Squirrel (Menetes berdmorei), Asian House Rat (Rattus tanezumi), Indomalayan Niviventer (Niviventer fulvescens), Savile's Bandicoot Rat (Bandicota savilei), (Tamiops mcclellandii), and Short-tailed Gymnure (Hylomys suillus) respectively. Considering the number of small mammals caught, Maxomys surifer was the most abundant, accounting for 61.69% of all mammals caught. The sex ratio was found to somewhat favor females in both seasons. This study indicated that the small mammal diversity was high in the primary dry evergreen forest. These results have important implications not only for conservation but are also useful for further investigation if there is any disturbance or change in the area or any potential disasters that may occur in the future.

Effect of enviromental factors at multiple landscape scales on bird community in riparian ecosystem at Mun-Chi River confluence, Thailand

Abstract. Chaleekarn W, Duengkae P, Pongcharoen C, Sutummawong N, Nakmuenwai P, Siripin S, Chirachitmichi C, Kummoo W, Paansri P, Suksavate W. 2022. Effect of environmental factors at multiple landscape scales on bird community in riparian ecosystem at Mun-Chi River confluence, Thailand. Biodiversitas 23: 5194-5204. Wetland and riparian ecosystem is an important migratory stopover for land and water birds in the East Asian - Australasian Flyway. Understanding relationship patterns between bird communities and environmental factors at multi-spatial scales within a landscape context could contribute to the conservation and management of bird biodiversity in wetland ecosystems. The landscape metrics index is critical in revealing the relationship between the composition of bird communities and habitats at both local and landscape scales. This study aims to determine the effect of the environmental factors at different designated spatial scales on the composition of local bird communities in terms of species and feeding guilds. Our study conducted a bird survey using 227-point transects along 40 tracks across different land cover types surrounding the Mun-Chi River confluence. Canonical correspondence analysis (CCA) was used to quantify the association between bird communities, represented by species and feeding guilds, and environmental factors with the integration of multilevel habitat metrics. From the results, the CCA showed patterns of the community-environmental association at multiple scales of patch, class, and landscape characteristics with the proportional explanation of 54% and 61.82% for the composition of species and feeding guilds, respectively. The results indicated the premise that the majority of bird species respond to the habitat at the local scale. Large forest patches can maintain migratory and resident bird species. Moreover, most avian groups were arranged primarily in a large forest core area, forest area, and Shrubland PA. The results confirmed existing information on feeding guilds. The prediction map of the principal component of avian species composition was created from the association with the drivers of land use, including crops, perennial farmland, and water body on the edge of forests. Therefore, wetland management must be done at both local and landscape scales to preserve suitable avian habitats.

Monetary evaluation of supporting ecosystem services as a habitat provider for birds in Thailand urban park

Abstract. Yarnvudhi A, Leksungnoen N, Siri S, Ponpithuk Y, Sukmasuang R, Duengkae P, Pongcharoen C, Sutummawong N, Marod D, Wachrinrat C, Premashthira A, Tor-ngern P, Poungcharean S, Hermuk S, Kachina P. 2022. Monetary evaluation of supporting ecosystem services as a habitat provider for birds in Thailand urban park. Biodiversitas 23: 4747-4758. Supporting services are important services to maintain ecosystems by providing habitats for organisms and genetic diversity. In this study, the monetary value of supporting services for bird habitats in urban parks was evaluated based on the price for nursing each bird species and market value. Bird diversity was conducted using point count observation. We found total of 53 bird species, 27 families and 8 orders in the park. Total monetary benefit for bird habitats in this park was estimated to be around USD60,354.12 per individual bird with an average value of USD1,138.76 per bird per species. The top three species with the highest monetary value were the Painted Stork (Mycteria leucocephala), Asian Openbill (Anastomus oscitans), and Chinese Pond Heron (Ardeola bacchus). The species diversity index of the bird community (H') was 2.73 and the most abundant bird species were Eastern Spotted Dove (Spilopelia chinensis), Eurasian tree sparrow (Passer montanus), Eastern jungle crow (Corvus levaillantii), Oriental magpie robin (Copsychus saularis), and Coppersmith Barbet (Psilopogon haemacephalus). Among 4 microhabitats, birds were found the most in the trees, followed by lawn, grassland, and wetland areas. The small urban parks should be designed with diverse microhabitats to provide various ecological functions to attract and ensure adequate resources for organisms.

Coincidence of low genetic diversity and increasing population size in wild gaur populations in the Khao Phaeng Ma Non-Hunting Area, Thailand: A challenge for conservation management under human-wildlife conflict

The gaur (Bos gaurus) is found throughout mainland South and Southeast Asia but is listed as an endangered species in Thailand with a decreasing population size and a reduction in suitable habitat. While gaur have shown a population recovery from 35 to 300 individuals within 30 years in the Khao Phaeng Ma (KPM) Non-Hunting Area, this has caused conflict with villagers along the border of the protected area. At the same time, the ecotourism potential of watching gaurs has boosted the local economy. In this study, 13 mitochondrial displacement-loop sequence samples taken from gaur with GPS collars were analyzed. Three haplotypes identified in the population were defined by only two parsimony informative sites (from 9 mutational steps of nucleotide difference). One haplotype was shared among eleven individuals located in different subpopulations/herds, suggesting very low genetic diversity with few maternal lineages in the founder population. Based on the current small number of sequences, neutrality and demographic expansion test results also showed that the population was likely to contract in the near future. These findings provide insight into the genetic diversity and demography of the wild gaur population in the KPM protected area that can inform long-term sustainable management action plans.

Mapping Risk of Nipah Virus Transmission from Bats to Humans in Thailand

Nipah virus (NiV) is a zoonotic virus that can pose a serious threat to human and livestock health. Old-world fruit bats (Pteropus spp.) are the natural reservoir hosts for NiV, and Pteropus lylei, Lyle's flying fox, is an important host of NiV in mainland Southeast Asia. NiV can be transmitted from bats to humans directly via bat-contaminated foods (i.e., date palm sap or fruit) or indirectly via livestock or other intermediate animal hosts. Here we construct risk maps for NiV spillover and transmission by combining ecological niche models for the P. lylei bat reservoir with other spatial data related to direct or indirect NiV transmission (livestock density, foodborne sources including fruit production, and human population). We predict the current and future (2050 and 2070) distribution of P. lylei across Thailand, Cambodia, and Vietnam. Our best-fit model predicted that central and western regions of Thailand and small areas in Cambodia are currently the most suitable habitats for P. lylei. However, due to climate change, the species range is predicted to expand to include lower northern, northeastern, eastern, and upper southern Thailand and almost all of Cambodia and lower southern Vietnam. This expansion will create additional risk areas for human infection from P. lylei in Thailand. Our combined predictive risk maps showed that central Thailand, inhabited by 2.3 million people, is considered highly suitable for the zoonotic transmission of NiV from P. lylei. These current and future NiV transmission risk maps can be used to prioritize sites for active virus surveillance and developing awareness and prevention programs to reduce the risk of NiV spillover and spread in Thailand.

Something Fishy about Siamese Fighting Fish (Betta splendens) Sex: Polygenic Sex Determination or a Newly Emerged Sex-Determining Region?

Fishes provide a unique and intriguing model system for studying the genomic origin and evolutionary mechanisms underlying sex determination and high sex-chromosome turnover. In this study, the mode of sex determination was investigated in Siamese fighting fish, a species of commercial importance. Genome-wide SNP analyses were performed on 75 individuals (40 males and 35 females) across commercial populations to determine candidate sex-specific/sex-linked loci. In total, 73 male-specific loci were identified and mapped to a 5.6 kb region on chromosome 9, suggesting a putative male-determining region (pMDR) containing localized dmrt1 and znrf3 functional sex developmental genes. Repeat annotations of the pMDR revealed an abundance of transposable elements, particularly Ty3/Gypsy and novel repeats. Remarkably, two out of the 73 male-specific loci were located on chromosomes 7 and 19, implying the existence of polygenic sex determination. Besides male-specific loci, five female-specific loci on chromosome 9 were also observed in certain populations, indicating the possibility of a female-determining region and the polygenic nature of sex determination. An alternative explanation is that male-specific loci derived from other chromosomes or female-specific loci in Siamese fighting fish recently emerged as new sex-determining loci during domestication and repeated hybridization.

Genome-Wide SNP Analysis of Hybrid Clariid Fish Reflects the Existence of Polygenic Sex-Determination in the Lineage

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.

Remnant of Unrelated Amniote Sex Chromosomal Linkage Sharing on the Same Chromosome in House Gecko Lizards, Providing a Better Understanding of the Ancestral Super-Sex Chromosome

Comparative chromosome maps investigating sex chromosomal linkage groups in amniotes and microsatellite repeat motifs of a male house gecko lizard (Hemidactylus frenatus, HFR) and a flat-tailed house gecko lizard (H. platyurus, HPL) of unknown sex were examined using 75 bacterial artificial chromosomes (BACs) from chicken and zebra finch genomes. No massive accumulations of microsatellite repeat motifs were found in either of the gecko lizards, but 10 out of 13 BACs mapped on HPL chromosomes were associated with other amniote sex chromosomes. Hybridization of the same BACs onto multiple different chromosome pairs suggested transitions to sex chromosomes across amniotes. No BAC hybridization signals were found on HFR chromosomes. However, HFR diverged from HPL about 30 million years ago, possibly due to intrachromosomal rearrangements occurring in the HFR lineage. By contrast, heterochromatin likely reshuffled patterns between HPL and HFR, as observed from C-positive heterochromatin distribution. Six out of ten BACs showed partial homology with squamate reptile chromosome 2 (SR2) and snake Z and/or W sex chromosomes. The gecko lizard showed shared unrelated sex chromosomal linkages-the remnants of a super-sex chromosome. A large ancestral super-sex chromosome showed a correlation between SR2 and snake W sex chromosomes.

Reduced genetic variability in a captive-bred population of the endangered Hume's pheasant (Syrmaticus humiae, Hume 1881) revealed by microsatellite genotyping and D-loop sequencing

Captive breeding programs are crucial to ensure the survival of endangered species and ultimately to reintroduce individuals into the wild. However, captive-bred populations can also deteriorate due to inbreeding depression and reduction of genetic variability. We genotyped a captive population of 82 individuals of the endangered Hume's pheasant (Syrmaticus humiae, Hume 1881) at the Doi Tung Wildlife Breeding Center to assess the genetic consequences associated with captive breeding. Analysis of microsatellite loci and mitochondrial D-loop sequences reveal significantly reduced genetic differentiation and a shallow population structure. Despite the low genetic variability, no bottleneck was observed but 12 microsatellite loci were informative in reflecting probable inbreeding. These findings provide a valuable source of knowledge to maximize genetic variability and enhance the success of future conservation plans for captive and wild populations of Hume's pheasant.

Mitochondrial genome of bronze-winged jacana (Metopidius indicus, Latham 1790)

We reported the mitochondrial genome (mitogenome) of bronze-winged jacana (Metopidius indicus, Latham 1790). The circular mitogenome was 17,208 base pairs (bp) in length, containing 13 protein-coding genes, two rRNAs, 22 tRNAs, and a non-coding control region. A DNA spacer 109 bp long was also detected between ND5 and Cytb. Phylogenetic analysis indicated that M. indicus was more closely related with the genera Himantopus, Jacana and Hydrophasianus. This annotated mitogenome reference can be utilized as a data resource for comparative mitogenomics of waders or shorebirds, with possible use in ecological and evolutionary studies.

Two decades of one health surveillance of Nipah virus in Thailand

BACKGROUND

Nipah virus (NiV) infection causes encephalitis and has > 75% mortality rate, making it a WHO priority pathogen due to its pandemic potential. There have been NiV outbreak(s) in Malaysia, India, Bangladesh, and southern Philippines. NiV naturally circulates among fruit bats of the genus Pteropus and has been detected widely across Southeast and South Asia. Both Malaysian and Bangladeshi NiV strains have been found in fruit bats in Thailand. This study summarizes 20 years of pre-emptive One Health surveillance of NiV in Thailand, including triangulated surveillance of bats, and humans and pigs in the vicinity of roosts inhabited by NiV-infected bats.

METHODS

Samples were collected periodically and tested for NiV from bats, pigs and healthy human volunteers from Wat Luang village, Chonburi province, home to the biggest P. lylei roosts in Thailand, and other provinces since 2001. Archived cerebrospinal fluid specimens from encephalitis patients between 2001 and 2012 were also tested for NiV. NiV RNA was detected using nested reverse transcription polymerase chain reaction (RT-PCR). NiV antibodies were detected using enzyme-linked immunosorbent assay or multiplex microsphere immunoassay.

RESULTS

NiV RNA (mainly Bangladesh strain) was detected every year in fruit bats by RT-PCR from 2002 to 2020. The whole genome sequence of NiV directly sequenced from bat urine in 2017 shared 99.17% identity to NiV from a Bangladeshi patient in 2004. No NiV-specific IgG antibodies or RNA have been found in healthy volunteers, encephalitis patients, or pigs to date. During the sample collection trips, 100 community members were trained on how to live safely with bats.

CONCLUSIONS

High identity shared between the NiV genome from Thai bats and the Bangladeshi patient highlights the outbreak potential of NiV in Thailand. Results from NiV cross-sectoral surveillance were conveyed to national authorities and villagers which led to preventive control measures, increased surveillance of pigs and humans in vicinity of known NiV-infected roosts, and increased vigilance and reduced risk behaviors at the community level. This proactive One Health approach to NiV surveillance is a success story; that increased collaboration between the human, animal, and wildlife sectors is imperative to staying ahead of a zoonotic disease outbreak.

High-Level Gene Flow Restricts Genetic Differentiation in Dairy Cattle Populations in Thailand: Insights from Large-Scale Mt D-Loop Sequencing

Domestication and artificial selection lead to the development of genetically divergent cattle breeds or hybrids that exhibit specific patterns of genetic diversity and population structure. Recently developed mitochondrial markers have allowed investigation of cattle diversity worldwide; however, an extensive study on the population-level genetic diversity and demography of dairy cattle in Thailand is still needed. Mitochondrial D-loop sequences were obtained from 179 individuals (hybrids of Bos taurus and B. indicus) sampled from nine different provinces. Fifty-one haplotypes, of which most were classified in haplogroup "I", were found across all nine populations. All sampled populations showed severely reduced degrees of genetic differentiation, and low nucleotide diversity was observed in populations from central Thailand. Populations that originated from adjacent geographical areas tended to show high gene flow, as revealed by patterns of weak network structuring. Mismatch distribution analysis was suggestive of a stable population, with the recent occurrence of a slight expansion event. The results provide insights into the origins and the genetic relationships among local Thai cattle breeds and will be useful for guiding management of cattle breeding in Thailand.

Complete mitochondrial genome of Mekong fighting fish, Betta smaragdina (Teleostei: Osphronemidae)

Mekong fighting fish (Betta smaragdina) are found in Northeast Thailand. A complete mitochondrial genome (mitogenome) of B. smaragdina was assembled and annotated. Mitogenome sequences were 16,372 bp in length, with slight AT bias (59.8%), containing 37 genes with identical order to most teleost mitogenomes. Phylogenetic analysis of B. smaragdina showed closer relationship with B. splendens and B. mahachaiensis as the bubble-nesting group, compared to the mouthbrooder group (B. apollon, B. simplex, and B. pi). Results will allow the creation of a reference annotated genome that can be utilized to sustain biodiversity and eco-management of betta bioresources to improve conservation programs.

Evidence for SARS-CoV-2 related coronaviruses circulating in bats and pangolins in Southeast Asia

Among the many questions unanswered for the COVID-19 pandemic are the origin of SARS-CoV-2 and the potential role of intermediate animal host(s) in the early animal-to-human transmission. The discovery of RaTG13 bat coronavirus in China suggested a high probability of a bat origin. Here we report molecular and serological evidence of SARS-CoV-2 related coronaviruses (SC2r-CoVs) actively circulating in bats in Southeast Asia. Whole genome sequences were obtained from five independent bats (Rhinolophus acuminatus) in a Thai cave yielding a single isolate (named RacCS203) which is most related to the RmYN02 isolate found in Rhinolophus malayanus in Yunnan, China. SARS-CoV-2 neutralizing antibodies were also detected in bats of the same colony and in a pangolin at a wildlife checkpoint in Southern Thailand. Antisera raised against the receptor binding domain (RBD) of RmYN02 was able to cross-neutralize SARS-CoV-2 despite the fact that the RBD of RacCS203 or RmYN02 failed to bind ACE2. Although the origin of the virus remains unresolved, our study extended the geographic distribution of genetically diverse SC2r-CoVs from Japan and China to Thailand over a 4800-km range. Cross-border surveillance is urgently needed to find the immediate progenitor virus of SARS-CoV-2.

A bat origin for SARS-CoV-2 has been proposed. Here, by sampling wild Rhinolophus acuminatus bats from Thailand, the authors identified a SARS-CoV-2-related coronavirus (SC2r-CoV), designated as RacCS203, with 91.5% genome similarity to SARS-CoV-2, and show that sera obtained from bats and Malayan pangolin neutralize SARS-CoV-2.

Origin and evolutionary history of domestic chickens inferred from a large population study of Thai red junglefowl and indigenous chickens

In this study, we aimed to elucidate the origin of domestic chickens and their evolutionary history over the course of their domestication. We conducted a large-scale genetic study using mitochondrial DNA D-loop sequences and 28 microsatellite DNA markers to investigate the diversity of 298 wild progenitor red junglefowl (Gallus gallus) across two subspecies (G. g. gallus and G. g. spadiceus) from 12 populations and 138 chickens from 10 chicken breeds indigenous to Thailand. Twenty-nine D-loop sequence haplotypes were newly identified: 14 and 17 for Thai indigenous chickens and red junglefowl, respectively. Bayesian clustering analysis with microsatellite markers also revealed high genetic diversity in the red junglefowl populations. These results suggest that the ancestral populations of Thai indigenous chickens were large, and that a part of the red junglefowl population gene pool was not involved in the domestication process. In addition, some haplogroups that are distributed in other countries of Southeast Asia were not observed in either the red junglefowls or the indigenous chickens examined in the present study, suggesting that chicken domestication occurred independently across multiple regions in Southeast Asia.

Next-generation sequencing yields complete mitochondrial genome assembly of peaceful betta fish, Betta imbellis (Teleostei: Osphronemidae)

The complete mitochondrial genome (mitogenome) of the peaceful betta (Betta imbellis) was obtained using next-generation sequencing. The sample of B. imbellis was collected from its native habitat in Southern Thailand. The mitogenome sequence was 16,897 bp in length, containing 37 genes with identical order to most teleost mitogenomes. Overall nucleotide base composition of the complete mitogenome was determined as AT bias. Phylogenetic analysis of B. imbellis showed a closer relationship with bubble-nesting fighting fish. This annotated mitogenome reference can be utilized as a bioresource for phylogenetic studies to support betta conservation programs.

Genome Complexity Reduction High-Throughput Genome Sequencing of Green Iguana (Iguana iguana) Reveal a Paradigm Shift in Understanding Sex-Chromosomal Linkages on Homomorphic X and Y Sex Chromosomes

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.

Partial Amniote Sex Chromosomal Linkage Homologies Shared on Snake W Sex Chromosomes Support the Ancestral Super-Sex Chromosome Evolution in Amniotes

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.

Complete mitochondrial genome of Mahachai betta, Betta mahachaiensis (Teleostei: Osphronemidae)

Mahachai bettas (Betta mahachaiensis) are distributed in areas of brackish water with Nipa Palms in Samut Sakhon, Thailand but urbanization is restricting their biodiversity. A complete mitochondrial genome (mitogenome) of B. mahachaiensis was determined to support conservation programs. Mitogenome sequences were 16,980 bp in length with slight AT bias (61.91%), containing 37 genes with identical order to most teleost mitogenomes. Phylogenetic analysis of B. mahachaiensis showed a closer relationship with B. splendens. Results will allow the creation of a reference annotated genome that can be utilized to sustain biodiversity and eco-management of the betta to improve conservation programs.

Predictive genetic plan for a captive population of the Chinese goral (Naemorhedus griseus) and prescriptive action for ex situ and in situ conservation management in Thailand

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, H_e = 0.455 ± 0.219; H_e < H_o) 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.

Take one step backward to move forward: Assessment of genetic diversity and population structure of captive Asian woolly-necked storks (Ciconia episcopus)

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 (N_e) 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 N_e 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.

The complete mitochondrial genome of Fea's muntjac (Muntiacus feae Thomas and Doria, 1889) with phylogenetic analysis†

The complete mitochondrial genome (mitogenome) of the rare Fea's muntjac (Muntiacus feae) was sequenced (GenBank accession nos. MG857662-MG857664). The mitogenome was found to be 16,355 bp in length with base compositions of 33.16% A, 24.59% C, 13.46% G, and 28.78% T and a GC content of 38.06%. The genome is comprised of 13 protein-coding genes (PCGs), 22 transfer RNA genes, 2 ribosomal RNA genes, and a control region (D-loop). Phylogenetic analysis revealed that Fea's muntjac is more closely related to Black muntjac (M. cronifrons) than to Red muntjac (M. muntjak). These data will be useful for further studies on the genetic diversity and molecular phylogenetic relationship of the genus Muntiacus.

Longitudinal study of age-specific pattern of coronavirus infection in Lyle's flying fox (Pteropus lylei) in Thailand

Background

Bats are natural reservoirs for several highly pathogenic and novel viruses including coronaviruses (CoVs) (mainly Alphacoronavirus and Betacoronavirus). Lyle’s flying fox (Pteropus lylei)‘s roosts and foraging sites are usually in the proximity to humans and animals. Knowledge about age-specific pattern of CoV infection in P. lylei, prevalence, and viral shedding at roosts and foraging sites may have an impact on infection-age-structure model to control CoV outbreak.

Methods

P. lylei bats were captured monthly during January–December 2012 for detection of CoV at three areas in Chonburi province; two human dwellings, S1 and S2, where few fruit trees were located with an open pig farm, 0.6 km and 5.5 km away from the bat roost, S3. Nested RT-PCR of RNA-dependent RNA polymerase (RdRp) gene from rectal swabs was used for CoV detection. The strain of CoV was confirmed by sequencing and phylogenetic analysis.

Results

CoV infection was found in both juveniles and adult bats between May and October (January, in adults only and April, in juveniles only). Of total rectal swab positives (68/367, 18.5%), ratio was higher in bats captured at S1 (11/44, 25.0%) and S2 (35/99, 35.4%) foraging sites than at roost (S3) (22/224, 9.8%). Juveniles (forearm length ≤ 136 mm) were found with more CoV infection than adults at all three sites; S1 (9/24, 37.5% vs 2/20, 10%), S2 (22/49, 44.9% vs 13/50, 26.0%), and S3 (10/30, 33.3% vs 12/194, 6.2%). The average BCI of CoV infected bats was significantly lower than uninfected bats. No gender difference related to infection was found at the sites. Phylogenetic analysis of conserved RdRp gene revealed that the detected CoVs belonged to group D betacoronavirus (n = 64) and alphacoronavirus (n = 4).

Conclusions

The fact that CoV infection and shedding was found in more juvenile than adult bats may suggest transmission from mother during peripartum period. Whether viral reactivation during parturition period or stress is responsible in maintaining transmission in the bat colony needs to be explored.