Lions are one of the world's most iconic megafauna, yet little is known about their temporal and spatial demographic history and population differentiation. We analyzed a genomic dataset of 20 specimens: two ca. 30,000-y-old cave lions (Panthera leo spelaea), 12 historic lions (Panthera leo leo/Panthera leo melanochaita) that lived between the 15th and 20th centuries outside the current geographic distribution of lions, and 6 present-day lions from Africa and India. We found that cave and modern lions shared an ancestor ca. 500,000 y ago and that the 2 lineages likely did not hybridize following their divergence. Within modern lions, we found 2 main lineages that diverged ca. 70,000 y ago, with clear evidence of subsequent gene flow. Our data also reveal a nearly complete absence of genetic diversity within Indian lions, probably due to well-documented extremely low effective population sizes in the recent past. Our results contribute toward the understanding of the evolutionary history of lions and complement conservation efforts to protect the diversity of this vulnerable species.
Two Malaysian infectious bronchitis virus isolates, MH5365/95 and V9/04 were characterized based on sequence and phylogenetic analyses of S1, S2, M, and N genes. Nucleotide sequence alignments revealed many point mutations, short deletions, and insertions in S1 region of both IBV isolates. Phylogenetic analysis of S1 gene and sequences analysis of M gene indicated that MH5365/95 and V9/04 belong to non-Massachusetts strain. However, both isolates share only 77% identity. Analysis based on S1 gene showed that MH5365/95 shared more than 87% identity to several Chinese strains. Meanwhile, V9/04 showed only 67-77% identity to all the previously studied IBV strains included in this study suggesting it is a variant of IBV isolate that is unique to Malaysia. Phylogenetic analysis suggests, although both isolates were isolated 10 years apart from different states in Malaysia, they shared a common origin. Analysis based on S2 and N genes indicated that both strains are highly related to each other, and there are fewer mutations which occurred in the respective genes.
Genetic structure and biodiversity of the medicinal plant Ficus deltoidea have rarely been scrutinized. To fill these lacunae, five varieties, consisting of 30 F. deltoidea accessions were collected across the country and studied on the basis of molecular and morphological data. Molecular analysis of the accessions was performed using nine Inter Simple Sequence Repeat (ISSR) markers, seven of which were detected as polymorphic markers. ISSR-based clustering generated four clusters supporting the geographical distribution of the accessions to some extent. The Jaccard's similarity coefficient implied the existence of low diversity (0.50-0.75) in the studied population. STRUCTURE analysis showed a low differentiation among the sampling sites, while a moderate varietal differentiation was unveiled with two main populations of F. deltoidea. Our observations confirmed the occurrence of gene flow among the accessions; however, the highest degree of this genetic interference was related to the three accessions of FDDJ10, FDTT16 and FDKT25. These three accessions may be the genetic intervarietal fusion points of the plant's population. Principal Components Analysis (PCA) relying on quantitative morphological characteristics resulted in two principal components with Eigenvalue >1 which made up 89.96% of the total variation. The cluster analysis performed by the eight quantitative characteristics led to grouping the accessions into four clusters with a Euclidean distance ranged between 0.06 and 1.10. Similarly, a four-cluster dendrogram was generated using qualitative traits. The qualitative characteristics were found to be more discriminating in the cluster and PCA analyses, while ISSRs were more informative on the evolution and genetic structure of the population.
Several anuran groups of Laurasian origin are each co-distributed in four isolated regions of the Northern Hemisphere: central/southern Europe and adjacent areas, Korean Peninsula and adjacent areas, Indo-Malaya, and southern North America. Similar distribution patterns have been observed in diverse animal and plant groups. Savage [Savage, J.M., 1973. The geographic distribution of frogs: patterns and predictions. In: Vial, J.L. (Ed.), Evolutionary Biology of the Anurans. University of Missouri Press, Columbia, pp. 351-445] hypothesized that the Miocene global cooling and increasing aridities in interiors of Eurasia and North America caused a southward displacement and range contraction of Laurasian frogs (and other groups). We use the frog genus Bombina to test Savage's biogeographical hypothesis. A phylogeny of Bombina is reconstructed based on three mitochondrial and two nuclear gene fragments. The genus is divided into three major clades: an Indo-Malaya clade includes B. fortinuptialis, B. lichuanensis, B. maxima, and B. microdeladigitora; a European clade includes B. bombina, B. pachypus, and B. variegata; and a Korean clade contains B. orientalis. The European and Korean clades form sister-group relationship. Molecular dating of the phylogenetic tree using the penalized likelihood and Bayesian analyses suggests that the divergence between the Indo-Malaya clade and other Bombina species occurred 5.9-28.6 million years ago. The split time between the European clade and the Korean clade is estimated at 5.1-20.9 million years ago. The divergence times of these clades are not significantly later than the timing of Miocene cooling and drying, and therefore can not reject Savage's hypothesis. Some other aspects of biogeography of Bombina also are discussed. The Korean Peninsula and the Shandong Peninsula might have supplied distinct southern refugia for B. orientalis during the Pleistocene glacial maxima. In the Indo-Malaya clade, the uplift of the Tibetan Plateau might have promoted the split between B. maxima and the other species.
Bamboo and rattan are widely grown for manufacturing, horticulture, and agroforestry. Bamboo and rattan production might help reduce poverty, boost economic growth, mitigate climate change, and protect the natural environment. Despite progress in research, sufficient molecular and genomic resources to study these species are lacking. We launched the Genome Atlas of Bamboo and Rattan (GABR) project, a comprehensive, coordinated international effort to accelerate understanding of bamboo and rattan genetics through genome analysis. GABR includes 2 core subprojects: Bamboo-T1K (Transcriptomes of 1000 Bamboos) and Rattan-G5 (Genomes of 5 Rattans), and several other subprojects. Here we describe the organization, directions, and status of GABR.
Data from three published studies of genetic variation at 18 microsatellite loci in water buffalo populations in China (18 swamp type, two river type), Nepal (one wild, one domestic river, one hybrid) and south-east Asia (eight swamp, three river) were combined so as to gain a broader understanding of genetic relationships among the populations and their demographic history. Mean numbers of alleles and expected heterozygosities were significantly different among populations. Estimates of θ (a measure of population differentiation) were significant among the swamp populations for all loci and among the river populations for most loci. Differentiation among the Chinese swamp populations (which was due primarily to just one population) was much less than among the south-east Asian. The Nepal wild animals, phenotypically swamp type but genetically like river type, are significantly different from all the domestic river populations and presumably represent the ancestral Bubalus arnee (possibly with some river-type introgression). Relationships among the swamp populations (D(A) genetic distances, principal component analysis and structure analyses) show the south-east Asian populations separated into two groups by the Chinese populations. Given these relationships and the patterns of genetic variability, we postulate that the swamp buffalo was domesticated in the region of the far south of China, northern Thailand and Indochina. Following domestication, it spread south through peninsular Malaysia to Sumatra, Java and Sulawesi, and north through China, and then to Taiwan, the Philippines and Borneo.
Despite the tremendous growth of the DNA sequencing data in the last decade, our understanding of the human genome is still in its infancy. To understand the implications of genetic variants in the light of population genetics and molecular evolution, we developed a database, PGG.SNV ( https://www.pggsnv.org ), which gives much higher weight to previously under-investigated indigenous populations in Asia. PGG.SNV archives 265 million SNVs across 220,147 present-day genomes and 1018 ancient genomes, including 1009 newly sequenced genomes, representing 977 global populations. Moreover, estimation of population genetic diversity and evolutionary parameters is available in PGG.SNV, a unique feature compared with other databases.
Disjunctive distributions across paleotropical regions in the Indian Ocean Basin (IOB) often invoke dispersal/vicariance debates. Exacum (Gentianaceae, tribe Exaceae) species are spread around the IOB, in Africa, Madagascar, Socotra, the Arabian peninsula, Sri Lanka, India, the Himalayas, mainland Southeast Asia including southern China and Malaysia, and northern Australia. The distribution of this genus was suggested to be a typical example of vicariance resulting from the breakup of the Gondwanan supercontinent. The molecular phylogeny of Exacum is in principle congruent with morphological conclusions and shows a pattern that resembles a vicariance scenario with rapid divergence among lineages, but our molecular dating analysis demonstrates that the radiation is too recent to be associated with the Gondwanan continental breakup. We used our dating analysis to test the results of DIVA and found that the program predicted impossible vicariance events. Ancestral area reconstruction suggests that Exacum originated in Madagascar, and divergence dating suggests its origin was not before the Eocene. The Madagascan progenitor, the most recent common ancestor of Exacum, colonized Sri Lanka and southern India via long-distance dispersals. This colonizer underwent an extensive range expansion and spread to Socotra-Arabia, northern India, and mainland Southeast Asia in the northern IOB when it was warm and humid in these regions. This widespread common ancestor retreated subsequently from most parts of these regions and survived in isolation in Socotra-Arabia, southern India-Sri Lanka, and perhaps mainland Southeast Asia, possibly as a consequence of drastic climatic changes, particularly the spreading drought during the Neogene. Secondary diversification from these surviving centers and Madagascar resulted in the extant main lineages of the genus. The vicariance-like pattern shown by the phylogeny appears to have resulted from long-distance dispersals followed by extensive range expansion and subsequent fragmentation. The extant African species E. oldenlandioides is confirmed to be recently dispersed from Madagascar.
he Chinese tiger frog Hoplobatrachus rugulosus is widely distributed in southern China, Malaysia, Myanmar, Thailand, and Vietnam. It is listed in Appendix II of CITES as the only Class II nationally-protected frog in China. The bred tiger frog known as the Thailand tiger frog, is also identified as H. rugulosus. Our analysis of the Cyt b gene showed high genetic divergence (13.8%) between wild and bred samples of tiger frog. Unexpected genetic divergence of the complete mt genome (14.0%) was also observed between wild and bred samples of tiger frog. Yet, the nuclear genes (NCX1, Rag1, Rhod, Tyr) showed little divergence between them. Despite this and their very similar morphology, the features of the mitochondrial genome including genetic divergence of other genes, different three-dimensional structures of ND5 proteins, and gene rearrangements indicate that H. rugulosus may be a cryptic species complex. Using Bayesian inference, maximum likelihood, and maximum parsimony analyses, Hoplobatrachus was resolved as a sister clade to Euphlyctis, and H. rugulosus (BT) as a sister clade to H. rugulosus (WT). We suggest that we should prevent Thailand tiger frogs (bred type) from escaping into wild environments lest they produce hybrids with Chinese tiger frogs (wild type).
Dragonflies of the genus Orthetrum are members of the suborder Anisoptera, family Libellulidae. There are species pairs whose members are not easily separated from each other by morphological characters. In the present study, the DNA nucleotide sequences of mitochondrial and nuclear genes were employed to elucidate the phylogeny and systematics of Orthetrum dragonflies. Phylogenetic analyses could not resolve the various subfamilies of the family Libellulidae unequivocally. The nuclear 28S rRNA gene is highly conserved and could not resolve congeneric species of Orthetrum. Individual mitochondrial genes (COI, COII, and 16S rRNA) and combination of these genes as well as the nuclear ITS1&2 genes clearly differentiate morphologically similar species, such as the reddish species pairs O. chrysis and O. testaceum, and the bluish-coloured species O. glaucum and O. luzonicum. This study also reveals distinct genetic lineages between O. pruinosum schneideri (occurring in Malaysia) and O. pruinosum neglectum (occurring north of Peninsular Malaysia from India to Japan), indicating these taxa are cryptic species.
Originating in Africa, the Zika virus (ZIKV) has spread to Asia, Pacific Islands and now to the Americas and beyond. Since the first isolation in 1947, ZIKV strains have been sampled at various times in the last 69 years, but this history has not been reflected in studying the patterns of mutation accumulation in their genomes. Implementing the viral history, we show that the ZIKV ancestor appeared sometime in 1930-1945 and, at that point, its mutation rate was probably less than 0.2 × 10-3/nucleotide site/year and subsequently increased significantly in most of its descendants. Sustaining a high mutation rate of 4 × 10-3/site/year throughout its evolution, the Ancestral Asian strain, which was sampled from a mosquito in Malaysia, accumulated 13 mutations in the 3'-untranslated region of RNA stem-loops prior to 1963, seven of which generate more stable stem-loop structures and are likely to inhibit cellular antiviral activities, including immune and RNA interference (RNAi) pathways. The seven mutations have been maintained in all Asian and American strains and may be responsible for serious medical problems we are facing today and offer testable hypotheses to examine their roles in molecular interactions during brain development.
Typhoid fever is an infectious disease of global importance that is caused by Salmonella enterica subsp. enterica serovar Typhi (S. Typhi). This disease causes an estimated 200,000 deaths per year and remains a serious global health threat. S. Typhi is strictly a human pathogen, and some recovered individuals become long-term carriers who continue to shed the bacteria in their faeces, thus becoming main reservoirs of infection.
The current outbreak of coronavirus disease-2019 (COVID-19) poses unprecedented challenges to global health1. The new coronavirus responsible for this outbreak-severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-shares high sequence identity to SARS-CoV and a bat coronavirus, RaTG132. Although bats may be the reservoir host for a variety of coronaviruses3,4, it remains unknown whether SARS-CoV-2 has additional host species. Here we show that a coronavirus, which we name pangolin-CoV, isolated from a Malayan pangolin has 100%, 98.6%, 97.8% and 90.7% amino acid identity with SARS-CoV-2 in the E, M, N and S proteins, respectively. In particular, the receptor-binding domain of the S protein of pangolin-CoV is almost identical to that of SARS-CoV-2, with one difference in a noncritical amino acid. Our comparative genomic analysis suggests that SARS-CoV-2 may have originated in the recombination of a virus similar to pangolin-CoV with one similar to RaTG13. Pangolin-CoV was detected in 17 out of the 25 Malayan pangolins that we analysed. Infected pangolins showed clinical signs and histological changes, and circulating antibodies against pangolin-CoV reacted with the S protein of SARS-CoV-2. The isolation of a coronavirus from pangolins that is closely related to SARS-CoV-2 suggests that these animals have the potential to act as an intermediate host of SARS-CoV-2. This newly identified coronavirus from pangolins-the most-trafficked mammal in the illegal wildlife trade-could represent a future threat to public health if wildlife trade is not effectively controlled.
As one of the great survivors of the plant kingdom, barnyard grasses (Echinochloa spp.) are the most noxious and common weeds in paddy ecosystems. Meanwhile, at least two Echinochloa species have been domesticated and cultivated as millets. In order to better understand the genomic forces driving the evolution of Echinochloa species toward weed and crop characteristics, we assemble genomes of three Echinochloa species (allohexaploid E. crus-galli and E. colona, and allotetraploid E. oryzicola) and re-sequence 737 accessions of barnyard grasses and millets from 16 rice-producing countries. Phylogenomic and comparative genomic analyses reveal the complex and reticulate evolution in the speciation of Echinochloa polyploids and provide evidence of constrained disease-related gene copy numbers in Echinochloa. A population-level investigation uncovers deep population differentiation for local adaptation, multiple target-site herbicide resistance mutations of barnyard grasses, and limited domestication of barnyard millets. Our results provide genomic insights into the dual roles of Echinochloa species as weeds and crops as well as essential resources for studying plant polyploidization, adaptation, precision weed control and millet improvements.
Bats are important flagship species for biodiversity research; however, diversity in Southeast Asia is considerably underestimated in the current checklists and field guides. Incorporation of DNA barcoding into surveys has revealed numerous species-level taxa overlooked by conventional methods. Inclusion of these taxa in inventories provides a more informative record of diversity, but is problematic as these species lack formal description. We investigated how frequently documented, but undescribed, bat taxa are encountered in Peninsular Malaysia. We discuss whether a barcode library provides a means of recognizing and recording these taxa across biodiversity inventories. Tissue was sampled from bats trapped at Pasir Raja, Dungun Terengganu, Peninsular Malaysia. The DNA was extracted and the COI barcode region amplified and sequenced. We identified 9 species-level taxa within our samples, based on analysis of the DNA barcodes. Six specimens matched to four previously documented taxa considered candidate species but currently lacking formal taxonomic status. This study confirms the high diversity of bats within Peninsular Malaysia (9 species in 13 samples) and demonstrates how DNA barcoding allows for inventory and documentation of known taxa lacking formal taxonomic status.
Mycobacteria a genus of Actinobacteria are widespread in nature ranging from soil-dwelling saprophytes to human and animal pathogens. The rate of growth has been a classifying factor for the Mycobacterium spp., dividing them into the rapid growers and the slow growers. Here we have performed a comparative genome study of mycobacterial species in order to get better understanding of their evolution, particularly to understand the distinction between the rapid and slow growers. Our study shows that the slow growers had generally gained and lost more genes compared to the rapid growers. The slow growers might haved eventually lost genes (LivFGMH operon, shaACDEFG genes and MspA porin) that could contribute to the slow growth rate of the slow growers. The genes gained and lost in mycobacteria had eventually helped these bacteria to adapt to different environments and have led to the evolution of the present day rapid and slow growers. Our results also show high number of Mycobacterium abscessus specific genes (811 genes) and some of them are associated with the known bacterial quorum sensing genes that might be important for Mycobacterium abscessus to adapt and survive in variety of unfavorable environments. Mycobacterium abscessus also does not contains genes involved in the bacterial defense system and together with the quorum sensing genes may have contributed to the high gene gain rate of Mycobacterium abscessus.
The aim of the present study was to analyze the sequence of the VP1 gene in enterovirus 71 (EV71) isolates and to explore their genetic evolution, so as to provide a scientific basis for the clinical prevention and treatment of hand, foot and mouth disease. The fecal samples of 590 patients with suspected hand, foot and mouth disease treated at Yan'an Hospital (Kunming, China) between January 2015 and December 2016 were collected and EV71 nucleic acid was detected by fluorescence PCR. The viral RNA of EV71-positive samples was extracted, the VP1 gene was amplified by PCR and the products were sequenced. The VP1 gene sequence was analyzed using DNAMAN and MEGA (version 4.0) software and homologous modeling was performed using Pymol software. A total of 50 EV71-positive samples were identified and the detection rate was 8.47% (50/590 cases). All of the 50 EV71 strains were of the C4 subtype. The genetic distance between the strains detected in the present study and EV71 strains detected in Beijing, Anhui and Malaysia was 0.01-0.03, while that between the strains detected in the present study and Australian strains was 2.11. Homologous modeling indicated that the amino acid sequence of the VP1 gene of the detected strains had a H144Y mutation. There was no significant genetic variation in the EV71 strain within the 2-year period. In conclusion, the EV71 strains detected in the present study was similar to that detected in Beijing, Anhui and Malaysia but different to that from Australia. A point mutation was present in the amino acid sequence of the VP1 gene.
We previously described three new Malaysian orthoreoviruses designated Pulau virus, Melaka virus and Kampar virus. Melaka and Kampar viruses were shown to cause respiratory disease in humans. These viruses, together with Nelson Bay virus, isolated from Australian bats, are tentatively classified as different strains within the species Pteropine orthoreovirus (PRV), formerly known as Nelson Bay orthoreovirus, based on the small (S) genome segments. Here we report the sequences of the large (L) and medium (M) segments, thus completing the whole-genome characterization of the four PRVs. All L and M segments were highly conserved in size and sequence. Conserved functional motifs previously identified in other orthoreovirus gene products were also found in the deduced proteins encoded by the cognate segments of these viruses. Detailed sequence analysis identified two genetic lineages divided into the Australian and Malaysian PRVs, and potential genetic reassortment among the M and S segments of the three Malaysian viruses.
The causative agents of leptospirosis are responsible for an emerging zoonotic disease worldwide. One of the major routes of transmission for leptospirosis is the natural environment contaminated with the urine of a wide range of reservoir animals. Soils and surface waters also host a high diversity of non-pathogenic Leptospira and species for which the virulence status is not clearly established. The genus Leptospira is currently divided into 35 species classified into three phylogenetic clusters, which supposedly correlate with the virulence of the bacteria. In this study, a total of 90 Leptospira strains isolated from different environments worldwide including Japan, Malaysia, New Caledonia, Algeria, mainland France, and the island of Mayotte in the Indian Ocean were sequenced. A comparison of average nucleotide identity (ANI) values of genomes of the 90 isolates and representative genomes of known species revealed 30 new Leptospira species. These data also supported the existence of two clades and 4 subclades. To avoid classification that strongly implies assumption on the virulence status of the lineages, we called them P1, P2, S1, S2. One of these subclades has not yet been described and is composed of Leptospira idonii and 4 novel species that are phylogenetically related to the saprophytes. We then investigated genome diversity and evolutionary relationships among members of the genus Leptospira by studying the pangenome and core gene sets. Our data enable the identification of genome features, genes and domains that are important for each subclade, thereby laying the foundation for refining the classification of this complex bacterial genus. We also shed light on atypical genomic features of a group of species that includes the species often associated with human infection, suggesting a specific and ongoing evolution of this group of species that will require more attention. In conclusion, we have uncovered a massive species diversity and revealed a novel subclade in environmental samples collected worldwide and we have redefined the classification of species in the genus. The implication of several new potentially infectious Leptospira species for human and animal health remains to be determined but our data also provide new insights into the emergence of virulence in the pathogenic species.
In the industrial processes, lipases are expected to operate at temperatures above 45 °C and could retain activity in organic solvents. Hence, a C-terminal truncated lipase from Staphylococcus epidermis AT2 (rT-M386) was engineered by directed evolution. A mutant with glycine-to-cysteine substitution (G210C) demonstrated a remarkable improvement of thermostability, whereby the mutation enhanced the activity five-fold when compared to the rT-M386 at 50 °C. The rT-M386 and G210C lipases were purified concurrently using GST-affinity chromatography. The biochemical and biophysical properties of both enzymes were investigated. The G210C lipase showed a higher optimum temperature (45 °C) and displayed a more prolonged half-life in the range of 40-60 °C as compared to rT-M386. Both lipases exhibited optimal activity and stability at pH 8. The G210C showed the highest stability in the presence of polar organic solvents at 50 °C compared to the rT-M386. Denatured protein analysis presented a significant change in the molecular ellipticity value above 60 °C, which verified the experimental result on the temperature and thermostability profile of G210C.