Genetic analysis of populations of the threatened bat Pteropus mariannus

Thirty years of conservation breeding: Assessing the genetic diversity of captive Livingstone's fruit bats

Fruit bats (genus Pteropus) are typically island-endemic species important in seed dispersal and reforestation that are vulnerable to increased extinction risk. An effective method of reducing extinction risk in vulnerable species that cannot be conserved in their native habitat is establishing an ex-situ captive breeding programme. Due to anthropogenic threats and low population numbers, in the early 1990s, a captive breeding programme was established at Jersey Zoo, British Isles, for Critically Endangered Livingstone's fruit bats (Pteropus livingstonii). Here we use six polymorphic microsatellite loci to assess genetic diversity in the captive breeding population of Livingstone's fruit bats (P. livingstonii), 30 years after the programme's establishment, investigating change over generations and comparing our findings with published data from the wild population. We found no significant difference between the genetic diversity in the captive and wild populations of Livingstone's fruit bats (P. livingstonii), in both expected heterozygosity and allelic richness. The captive population has retained a comparable level of genetic diversity to that documented in the wild, and there has been no significant decline in genetic diversity over the last 30 years. We advise that a full pedigree of the paternal lineage is created to improve the management of the captive breeding programme and further reduce the possibility of inbreeding. However, it appears that the captive breeding programme is currently effective at maintaining genetic diversity at levels comparable to those seen in the wild population, which suggests reintroductions could be viable if genetic diversity remains stable in captivity.

The taxonomic basis of subspecies listed as threatened and endangered under the endangered species act

More than 170 subspecies are listed as threatened or endangered under the US Endangered Species Act. Most of these subspecies were described decades ago on the basis of geographical variation in morphology using relatively primitive taxonomic methods. The US Fish and Wildlife Service defaults to subspecies descriptions by taxonomists working with specific groups of organisms, but there is no single definition of subspecies across plants and animals. Valid tests today usually entail molecular analyses of variation within and among populations, although there is no reason that behavioral, ecological or molecular characters could not be used, and include tests for significant differences between samples of the putative endangered subspecies and its nearest geographic relatives. We evaluated data gathered since subspecies listed under the ESA were described finding about one-third are valid (distinct evolutionary taxa), one-third are not, and one-third have not been tested. Therefore, it should not be assumed that because a subspecies occurs in a checklist, it is taxonomically valid. If the US Fish and Wildlife Service intends to continue listing subspecies, we suggest that they convene taxonomic experts representing various groups of organisms to provide a minimal set of criteria for a subspecies to be listed under the ESA.

Genetic Differentiation and Demographic Trajectory of the Insular Formosan and Orii's Flying Foxes

Insular flying foxes are keystone species in island ecosystems due to their critical roles in plant pollination and seed dispersal. These species are vulnerable to population decline because of their small populations and low reproductive rates. The Formosan flying fox (Pteropus dasymallus formosus) is one of the 5 subspecies of the Ryukyu flying fox. Pteropus dasymallus formosus has suffered from a severe decline and is currently recognized as a critically endangered population in Taiwan. On the contrary, the Orii's flying fox (Pteropus dasymallus inopinatus) is a relatively stable population inhabiting Okinawa Island. Here, we applied a genomic approach called double digest restriction-site associated DNA sequencing to study these 2 subspecies for a total of 7 individuals. We detected significant genetic structure between the 2 populations. Despite their contrasting contemporary population sizes, both populations harbor very low degrees of genetic diversity. We further inferred their demographic history based on the joint folded site frequency spectrum and revealed that both P. d. formosus and P. d. inopinatus had maintained small population sizes for a long period of time after their divergence. Recently, these populations experienced distinct trajectories of demographic changes. While P. d. formosus suffered from a drastic ~10-fold population decline not long ago, P. d. inopinatus underwent a ~4.5-fold population expansion. Our results suggest separate conservation management for the 2 populations-population recovery is urgently needed for P. d. formosus while long-term monitoring for adverse genetic effects should be considered for P. d. inopinatus.

Genetic diversity and population structure in the Ryukyu flying fox inferred from remote sampling in the Yaeyama archipelago

The Ryukyu flying fox (Pteropus dasymallus) is distributed throughout the island chain spanning across southern Japan, Taiwan, and possibly the Philippines. Although P. dasymallus is listed as VU (vulnerable) in the IUCN Red List, only few genetic works have been conducted to support its conservation. In this study we analyzed 19 markers (mtDNA haplotypes and 18 microsatellite markers) to evaluate genetic diversity and investigate the genetic structure of this species. mtDNA analysis was conducted with 142 DNA remote samples, mostly from faeces, and wing tissues collected on eight islands (Miyako, Ishigaki, Kohama, Kuroshima, Hateruma, Taketomi, Iriomote, Yonaguni). 39 haplotypes were identified in 526bp of the control region, and haplotype network showed no clear genetic structure. Microsatellite analysis was also conducted with 155 samples collected on six islands (Miyako, Ishigaki, Kohama, Taketomi, Iriomote, Yonaguni). It showed that the Yonaguni population exhibits low genetic diversity, high inbreeding, and clear genetic differentiation from other populations. Gene flow between Ishigaki and Miyako through small stepstone islands might be preventing inbreeding of the Miyako population. We provide for the first time indirect proof of long-distance inter-island dispersal in the Ryukyu flying fox and revealed genetic diversity, gene flow and genetic differentiation among the archipelago’s populations. These results will be useful for delineating conservation units and designing specific conservation policies for each island based on metapopulation genetic structure.

Genetic diversity of toll-like receptor genes in the vulnerable Chinese egret (Egretta eulophotes)

Toll-like receptor (TLR) genes have recently been employed to assess genetic diversity, as they can be used to infer both demographic history and adaptation to environments with different pathogen pressure. Here, we sampled 120 individuals of the Chinese egret (Egretta eulophotes), a globally vulnerable species, from four breeding populations across China. We assessed the levels of genetic diversity, selection pressure, and population differentiation at seven TLR loci (TLR1LB, TLR2A, TLR3, TLR4, TLR5, TLR7, and TLR15). Using a variety of metrics (SNPs, heterozygosity, nucleotides, haplotypes), our analyses showed that genetic diversity was lower at 4 of the 7 TLR loci in the vulnerable Chinese egret compared to the more common little egret (Egretta garzetta). The selection test indicated TLRs, except for TLR5, were under purifying selection in TLR evolution, suggesting that low TLR genetic diversity in the Chinese egret may be caused by purifying selection. Moreover, analysis of molecular variance indicated low but significant population differentiation among four populations at all of the TLR loci in this egret. However, some comparisons based on fixation index analyses did not show significant population differentiation, and Bayesian clustering showed admixture. Our finding suggested that these four populations of the Chinese egret in China may be considered a single unit for conservation planning. These results, the new report of TLR genetic diversity in a long-distance migratory vulnerable Ardeid species, will provide fundamental TLR information for further studies on the conservation genetics of the Chinese egret and other Ardeids.

Dispersal out of Wallacea spurs diversification of Pteropus flying foxes, the world's largest bats (Mammalia: Chiroptera)

AIM

Islands provide opportunities for isolation and speciation. Many landmasses in the Indo-Australian Archipelago (IAA) are oceanic islands, and founder-event speciation is expected to be the predominant form of speciation of volant taxa on these islands. We studied the biogeographic history of flying foxes, a group with many endemic species and a predilection for islands, to test this hypothesis and infer the biogeographic origin of the group.

LOCATION

Australasia, Indo-Australian Archipelago, Madagascar, Pacific Islands.

TAXON

Pteropus (Pteropodidae).

METHODS

To infer the biogeographic history of Pteropus, we sequenced up to 6169 bp of genetic data from 10 markers and reconstructed a multilocus species tree of 34 currently recognized Pteropus species and subspecies with 3 Acerodon outgroups using BEAST and subsequently estimated ancestral areas using models implemented in BioGeoBEARS.

RESULTS

Species-level resolution was occasionally low because of slow rates of molecular evolution and/or recent divergences. Older divergences, however, were more strongly supported and allow the evolutionary history of the group to be inferred. The genus diverged in Wallacea from its common ancestor with Acerodon; founder-event speciation out of Wallacea was a common inference. Pteropus species in Micronesia and the western Indian Ocean were also inferred to result from founder-event speciation.

MAIN CONCLUSIONS

Dispersal between regions of the IAA and the islands found therein fostered diversification of Pteropus throughout the IAA and beyond. Dispersal in Pteropus is far higher than in most other volant taxa studied to date, highlighting the importance of inter-island movement in the biogeographic history of this large clade of large bats.

Population genetics of fruit bat reservoir informs the dynamics, distribution and diversity of Nipah virus

The structure and connectivity of wildlife host populations may influence zoonotic disease dynamics, evolution and therefore spillover risk to people. Fruit bats in the genus Pteropus, or flying foxes, are the primary natural reservoir for henipaviruses-a group of emerging paramyxoviruses that threaten livestock and public health. In Bangladesh, Pteropus medius is the reservoir for Nipah virus-and viral spillover has led to human fatalities nearly every year since 2001. Here, we use mitochondrial DNA and nuclear microsatellite markers to measure the population structure, demographic history and phylogeography of P. medius in Bangladesh. We combine this with a phylogeographic analysis of all known Nipah virus sequences and strains currently available to better inform the dynamics, distribution and evolutionary history of Nipah virus. We show that P. medius is primarily panmictic, but combined analysis of microsatellite and morphological data shows evidence for differentiation of two populations in eastern Bangladesh, corresponding to a divergent strain of Nipah virus also found in bats from eastern Bangladesh. Our demographic analyses indicate that a large, expanding population of flying foxes has existed in Bangladesh since the Late Pleistocene, coinciding with human population expansion in South Asia, suggesting repeated historical spillover of Nipah virus likely occurred. We present the first evidence of mitochondrial introgression, or hybridization, between P. medius and flying fox species found in South-East Asia (P. vampyrus and P. hypomelanus), which may help to explain the distribution of Nipah virus strains across the region.

Characterization of the mitochondrial genome of an ancient amphipod Halice sp. MT-2017 (Pardaliscidae) from 10,908 m in the Mariana Trench

Small amphipods (Halice sp. MT-2017) with body length <1 cm were collected from the Challenger Deep (~10,920 m below sea level). The divergence time of their lineage was approximately 109 Mya, making this group ancient compared to others under study. The mitochondrial genome of Halice sp. shared the usual gene components of metazoans, comprising 13 protein coding genes (PCGs), 22 transfer RNAs (tRNAs), and 2 ribosomal RNAs (rRNAs). The arrangement of these genes, however, differed greatly from that of other amphipods. Of the 15 genes that were rearranged with respect to the pancrustacean gene pattern, 12 genes (2 PCGs, 2 rRNAs, and 8 tRNAs) were both translocated and strand-reversed. In contrast, the mitochondrial genomes in other amphipods never show so many reordered genes, and in most instances, only tRNAs were involved in strand-reversion-coupled translocation. Other characteristics, including reversed strand nucleotide composition bias, relatively higher composition of non-polar amino acids, and lower evolutionary rate, were also identified. Interestingly, the latter two features were shared with another hadal amphipod, Hirondellea gigas, suggesting their possible associations with the adaptation to deep-sea extreme habitats. Overall, our data provided a useful resource for future studies on the evolutionary and adaptive mechanisms of hadal faunas.

Unmasking the complexity of species identification in Australasian flying-foxes

Pteropus (flying-foxes) are a speciose group of non-echolocating large bats, with five extant Australian species and 24 additional species distributed amongst the Pacific Islands. In 2015, an injured flying-fox with unusual facial markings was found in Sydney, Australia, following severe and widespread storms. Based on an initial assessment, the individual belonged to Pteropus but could not be readily identified to species. As a consequence, four hypotheses for its identification/origin were posited: the specimen represented (1) an undescribed Australian species; or (2) a morphological variant of a recognised Australian species; or (3) a hybrid individual; or (4) a vagrant from the nearby Southwest Pacific Islands. We used a combination of morphological and both mitochondrial- and nuclear DNA-based identification methods to assess these hypotheses. Based on the results, we propose that this morphologically unique Pteropus most likely represents an unusual P. alecto (black flying-fox) potentially resulting from introgression from another Pteropus species. Unexpectedly, this individual, and the addition of reference sequence data from newly vouchered specimens, revealed a previously unreported P. alecto mitochondrial DNA lineage. This lineage was distinct from currently available haplotypes. It also suggests long-term hybridisation commonly occurs between P. alecto and P. conspicillatus (spectacled flying-fox). This highlights the importance of extensive reference data, and the inclusion of multiple vouchered specimens for each species to encompass both intraspecific and interspecific variation to provide accurate and robust species identification. Moreover, our additional reference data further demonstrates the complexity of Pteropus species relationships, including hybridisation, and potential intraspecific biogeographical structure that may impact on their management and conservation.

The mitochondrial genome of a sea anemone Bolocera sp. exhibits novel genetic structures potentially involved in adaptation to the deep-sea environment

The deep sea is one of the most extensive ecosystems on earth. Organisms living there survive in an extremely harsh environment, and their mitochondrial energy metabolism might be a result of evolution. As one of the most important organelles, mitochondria generate energy through energy metabolism and play an important role in almost all biological activities. In this study, the mitogenome of a deep‐sea sea anemone (Bolocera sp.) was sequenced and characterized. Like other metazoans, it contained 13 energy pathway protein‐coding genes and two ribosomal RNAs. However, it also exhibited some unique features: just two transfer RNA genes, two group I introns, two transposon‐like noncanonical open reading frames (ORFs), and a control region‐like (CR‐like) element. All of the mitochondrial genes were coded by the same strand (the H‐strand). The genetic order and orientation were identical to those of most sequenced actiniarians. Phylogenetic analyses showed that this species was closely related to Bolocera tuediae. Positive selection analysis showed that three residues (31 L and 42 N in ATP6, 570 S in ND5) of Bolocera sp. were positively selected sites. By comparing these features with those of shallow sea anemone species, we deduced that these novel gene features may influence the activity of mitochondrial genes. This study may provide some clues regarding the adaptation of Bolocera sp. to the deep‐sea environment.

Taxonomy, distribution, and natural history of flying foxes (Chiroptera, Pteropodidae) in the Mortlock Islands and Chuuk State, Caroline Islands

The taxonomy, biology, and population status of flying foxes (Pteropus spp.) remain little investigated in the Caroline Islands, Micronesia, where multiple endemic taxa occur. Our study evaluated the taxonomic relationships between the flying foxes of the Mortlock Islands (a subgroup of the Carolines) and two closely related taxa from elsewhere in the region, and involved the first ever field study of the Mortlock population. Through a review of historical literature, the name Pteropus pelagicus Kittlitz, 1836 is resurrected to replace the prevailing but younger name Pteropus phaeocephalus Thomas, 1882 for the flying fox of the Mortlocks. On the basis of cranial and external morphological comparisons, Pteropus pelagicus is united taxonomically with Pteropus insularis "Hombron and Jacquinot, 1842" (with authority herein emended to Jacquinot and Pucheran 1853), and the two formerly monotypic species are now treated as subspecies - Pteropus pelagicus pelagicus in the Mortlocks, and Pteropus phaeocephalus insularis on the islands of Chuuk Lagoon and Namonuito Atoll. The closest relative of Pteropus pelagicus is Pteropus tokudae Tate, 1934, of Guam, which is best regarded as a distinct species. Pteropus pelagicus pelagicus is the only known resident bat in the Mortlock Islands, a chain of more than 100 atoll islands with a total land area of <12 km(2). Based on field observations in 2004, we estimated a population size of 925-1,200 bats, most of which occurred on Satawan and Lukunor Atolls, the two largest and southernmost atolls in the chain. Bats were absent on Nama Island and possibly extirpated from Losap Atoll in the northern Mortlocks. Resident Mortlockese indicated bats were more common in the past, but that the population generally has remained stable in recent years. Most Pteropus phaeocephalus pelagicus roosted alone or in groups of 5-10 bats; a roost of 27 was the largest noted. Diet is comprised of at least eight plant species, with breadfruit (Artocarpus spp.) being a preferred food. Records of females with young (April, July) and pregnant females (July) suggest an extended breeding season. Pteropus pelagicus pelagicus appears most threatened by the prospect of sea level rise associated with global climate change, which has the potential to submerge or reduce the size of atolls in the Mortlocks. Occasional severe typhoons probably temporarily reduce populations on heavily damaged atolls, but hunting and ongoing habitat loss are not current problems for the subspecies.

The complete mitochondrial genome sequence and gene organization of the mud crab (Scylla paramamosain) with phylogenetic consideration

The complete mitochondrial genome is of great importance for better understanding the genome-level characteristics and phylogenetic relationships among related species. In the present study, we determined the complete mitochondrial genome DNA sequence of the mud crab (Scylla paramamosain) by 454 deep sequencing and Sanger sequencing approaches. The complete genome DNA was 15,824 bp in length and contained a typical set of 13 protein-coding genes, 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes and a putative control region (CR). Of 37 genes, twenty-three were encoded by the heavy strand (H-strand), while the other ones were encoded by light strand (L-strand). The gene order in the mitochondrial genome was largely identical to those obtained in most arthropods, although the relative position of gene tRNA(His) differed from other arthropods. Among 13 protein-coding genes, three (ATPase subunit 6 (ATP6), NADH dehydrogenase subunits 1 (ND1) and ND3) started with a rare start codon ATT, whereas, one gene cytochrome c oxidase subunit I (COI) ended with the incomplete stop codon TA. All 22 tRNAs could fold into a typical clover-leaf secondary structure, with the gene sizes ranging from 63 to 73 bp. The phylogenetic analysis based on 12 concatenated protein-coding genes showed that the molecular genetic relationship of 19 species of 11 genera was identical to the traditional taxonomy.