The genetic background of the cynomolgus macaque (Macaca fascicularis) is made complex by the high genetic diversity, population structure, and gene introgression from the closely related rhesus macaque (Macaca mulatta). Herein we report the whole-genome sequence of a Malaysian cynomolgus macaque male with more than 40-fold coverage, which was determined using a resequencing method based on the Indian rhesus macaque genome.
Reassortment of genetic segments between and within influenza B lineages (Victoria and Yamagata) has been shown to generate novel reassortants with unique genetic characteristics. Based on hemagglutinin (HA) and neuraminidase (NA) genes, recent surveillance study has identified reassortment properties in B/Phuket/3073/2013-like virus, which is currently used in the WHO-recommended influenza vaccine. To understand the potential reassortment patterns for all gene segments, four B/Phuket/3073/2013-like viruses and two unique reassortants (one each from Yamagata and Victoria) detected in Malaysia from 2012-2014 were subjected to whole-genome sequencing. Each gene was phylogenetically classified into lineages, clades and sub-clades. Three B/Phuket/3073/2013-like viruses from Yamagata lineage were found to be intra-clade reassortants, possessing PA and NA genes derived from Stockholm/12-like sub-clade, while the remaining genes from Wisconsin/01-like sub-clade (both sub-clades were within Yamagata Clade 3/Yam-3). However, the other B/Phuket/3073/2013-like virus had NS gene that derived from Stockholm/12-like sub-clade instead of Wisconsin/01-like sub-clade. One inter-clade reassortant had Yamagata Clade 2/Yam-2-derived HA and NP, and its remaining genes were Yam-3-derived. Within Victoria Clade 1/Vic-1 in Victoria lineage, one virus had intra-clade reassortment properties: HA and PB2 from Vic-1B sub-clade, MP and NS from a unique sub-clade "Vic-1C", and the remaining genes from Vic-1A sub-clade. Although random reassortment event may generate unique reassortants, detailed phylogenetic classification of gene segments showed possible genetic linkage between PA and NA genes in B/Phuket/3073/2013-like viruses, which requires further investigation. Understanding on reassortment patterns in influenza B evolution may contribute to future vaccine design.
Here we report the draft genomes and annotation of four N-acyl homoserine lactone (AHL)-producing members from the family Sphingomonadaceae. Comparative genomic analyses of 62 Sphingomonadaceae genomes were performed to gain insights into the distribution of the canonical luxI/R-type quorum sensing (QS) network within this family. Forty genomes contained at least one luxR homolog while the genome of Sphingobium yanoikuyae B1 contained seven Open Reading Frames (ORFs) that have significant homology to that of luxR. Thirty-three genomes contained at least one luxI homolog while the genomes of Sphingobium sp. SYK6, Sphingobium japonicum, and Sphingobium lactosutens contained four luxI. Using phylogenetic analysis, the sphingomonad LuxR homologs formed five distinct clades with two minor clades located near the plant associated bacteria (PAB) LuxR solo clade. This work for the first time shows that 13 Sphingobium and one Sphingomonas genome(s) contain three convergently oriented genes composed of two tandem luxR genes proximal to one luxI (luxR-luxR-luxI). Interestingly, luxI solos were identified in two Sphingobium species and may represent species that contribute to AHL-based QS system by contributing AHL molecules but are unable to perceive AHLs as signals. This work provides the most comprehensive description of the luxI/R circuitry and genome-based taxonomical description of the available sphingomonad genomes to date indicating that the presence of luxR solos and luxI solos are not an uncommon feature in members of the Sphingomonadaceae family.
A report on the de novo Whole Genome Sequence (WGS) of Marinobacter litoralis SW-45, a moderately salt-tolerant bacterium isolated from the seawater in Malaysia is presented. The strain has a genome size of 3.45 Mb and is capable of producing halophilic lipase, protease and esterase enzymes. Computational prediction of non-coding RNA (ncRNA) genes in M. litoralis SW-45 was performed using standalone software known as the non-coding RNA characterization (nocoRNAc). In addition, a phylogenetic tree showing the evolutionary relationship between the strain and other members of the genus Marinobacter was constructed using 16SrRNA sequence information. A total of 385 ncRNA transcripts, 1124 terminator region, and 2350 Stress Induced Duplex Destabilization sites were predicted. The current WGS shotgun project has provided the relevant genetic information that may be useful for the strain's improvement studies. This manuscript gives the first description of M. litoralis with a complete genome.
Cases of autochthonous infections of dengue virus type 1 (DENV-1) were detected in Japan after a 70-year period devoid of dengue outbreaks. We previously showed that E gene sequences are identical in 11 of the 12 DENV-1 strains autochthonous to Japan. However, the E sequence represents only 14% of the DENV-1 genome. In the present study, we have sequenced the entire genome of 6 autochthonous DENV-1 strains that were isolated from patients during the 2014 outbreak. Sequencing of 5 Yoyogi group strains with identical E sequences and 1 Shizuoka strain with a different E sequence revealed that the first Yoyogi group strain differed from the Shizuoka strain by 18 amino acid residues. Furthermore, 2 Yoyogi group strains had different genomic sequences while the other 3 had identical genomes. Phylogenetic analyses indicated that the Hyogo strain, a Yoyogi group strain, was the first to diverge from the other 4 Yoyogi group strains. The E gene sequence of the Yoyogi group strains exhibits the highest homology to those of the strains isolated in Malaysia and Singapore between 2013 and 2014. The patient infected with the Hyogo strain visited Malaysia before the onset of dengue fever, suggesting that this was a case of dengue infection imported from Malaysia.
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.
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.
Epidemics caused by fungal plant pathogens pose a major threat to agro-ecosystems and impact global food security. High-throughput sequencing enabled major advances in understanding how pathogens cause disease on crops. Hundreds of fungal genomes are now available and analyzing these genomes highlighted the key role of effector genes in disease. Effectors are small secreted proteins that enhance infection by manipulating host metabolism. Fungal genomes carry 100s of putative effector genes, but the lack of homology among effector genes, even for closely related species, challenges evolutionary and functional analyses. Furthermore, effector genes are often found in rapidly evolving chromosome compartments which are difficult to assemble. We review how population and comparative genomics toolsets can be combined to address these challenges. We highlight studies that associated genome-scale polymorphisms with pathogen lifestyles and adaptation to different environments. We show how genome-wide association studies can be used to identify effectors and other pathogenicity-related genes underlying rapid adaptation. We also discuss how the compartmentalization of fungal genomes into core and accessory regions shapes the evolution of effector genes. We argue that an understanding of genome evolution provides important insight into the trajectory of host-pathogen co-evolution.
The suppression of recombination is considered a hallmark of sex chromosome evolution. However, previous research has identified undifferentiated sex chromosomes and sex determination by single SNP in the greater amberjack (Seriola dumerili). We observed the same phenomena in the golden pompano (Trachinotus ovatus) of the same family Carangidae and discovered a different sex-determining SNP within the same gene Hsd17b1. We propose an evolutionary model elucidating the turnover of sex-determining mutations by highlighting the contrasting dynamics between purifying selection, responsible for maintaining W-linked Hsd17b1, and neutral evolution, which drives Z-linked Hsd17b1. Additionally, sporadic loss-of-function mutations in W-linked Hsd17b1 contribute to the conversion of W chromosomes into Z chromosomes. This model was directly supported by simulations, closely related species, and indirectly by zebrafish mutants. These findings shed new light on the early stages of sex chromosome evolution.
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.
In infections by apicomplexan parasites including Plasmodium, Toxoplasma gondii, and Eimeria, host interactions are mediated by proteins including families of membrane-anchored cysteine-rich surface antigens (SAGs) and SAG-related sequences (SRS). Eimeria tenella causes caecal coccidiosis in chickens and has a SAG family with over 80 members making up 1% of the proteome. We have solved the structure of a representative E. tenella SAG, EtSAG19, revealing that, despite a low level of sequence similarity, the entire Eimeria SAG family is unified by its three-layer αβα fold which is related to that of the CAP superfamily. Furthermore, sequence comparisons show that the Eimeria SAG fold is conserved in surface antigens of the human coccidial parasite Cyclospora cayetanensis but this fold is unrelated to that of the SAGs/SRS proteins expressed in other apicomplexans including Plasmodium species and the cyst-forming coccidia Toxoplasma gondii, Neospora caninum and Besnoitia besnoiti. However, despite having very different structures, Consurf analysis showed that Eimeria SAG and Toxoplasma SRS families each exhibit marked hotspots of sequence hypervariability that map to their surfaces distal to the membrane anchor. This suggests that the primary and convergent purpose of the different structures is to provide a platform onto which sequence variability can be imposed.
Peninsular Malaysia is a strategic region which might have played an important role in the initial peopling and subsequent human migrations in Asia. However, the genetic diversity and history of human populations--especially indigenous populations--inhabiting this area remain poorly understood. Here, we conducted a genome-wide study using over 900,000 single nucleotide polymorphisms (SNPs) in four major Malaysian ethnic groups (MEGs; Malay, Proto-Malay, Senoi and Negrito), and made comparisons of 17 world-wide populations. Our data revealed that Peninsular Malaysia has greater genetic diversity corresponding to its role as a contact zone of both early and recent human migrations in Asia. However, each single Orang Asli (indigenous) group was less diverse with a smaller effective population size (N(e)) than a European or an East Asian population, indicating a substantial isolation of some duration for these groups. All four MEGs were genetically more similar to Asian populations than to other continental groups, and the divergence time between MEGs and East Asian populations (12,000--6,000 years ago) was also much shorter than that between East Asians and Europeans. Thus, Malaysian Orang Asli groups, despite their significantly different features, may share a common origin with the other Asian groups. Nevertheless, we identified traces of recent gene flow from non-Asians to MEGs. Finally, natural selection signatures were detected in a batch of genes associated with immune response, human height, skin pigmentation, hair and facial morphology and blood pressure in MEGs. Notable examples include SYN3 which is associated with human height in all Orang Asli groups, a height-related gene (PNPT1) and two blood pressure-related genes (CDH13 and PAX5) in Negritos. We conclude that a long isolation period, subsequent gene flow and local adaptations have jointly shaped the genetic architectures of MEGs, and this study provides insight into the peopling and human migration history in Southeast Asia.
Salmonellae display intricate evolutionary patterns comprising over 2500 serovars having diverse pathogenic profiles. The acquisition and/or exchange of various virulence factors influences the evolutionary framework. To gain insights into evolution of Salmonella in association with the CRISPR-Cas genes we performed phylogenetic surveillance across strains of 22 Salmonella serovars. The strains differed in their CRISPR1-leader and cas operon features assorting into two main clades, CRISPR1-STY/cas-STY and CRISPR1-STM/cas-STM, comprising majorly typhoidal and non-typhoidal Salmonella serovars respectively. Serovars of these two clades displayed better relatedness, concerning CRISPR1-leader and cas operon, across genera than between themselves. This signifies the acquisition of CRISPR1/Cas region could be through a horizontal gene transfer event owing to the presence of mobile genetic elements flanking CRISPR1 array. Comparison of CRISPR and cas phenograms with that of multilocus sequence typing (MLST) suggests differential evolution of CRISPR/Cas system. As opposed to broad-host-range, the host-specific serovars harbor fewer spacers. Mapping of protospacer sources suggested a partial correlation of spacer content with habitat diversity of the serovars. Some serovars like serovar Enteritidis and Typhimurium that inhabit similar environment/infect similar hosts hardly shared their protospacer sources.
DARPP-32 (dopamine and adenosine 3', 5'-monophosphate-regulated phosphoprotein of 32 kDa), which belongs to PPP1R1 gene family, is known to act as an important integrator in dopamine-mediated neurotransmission via the inhibition of protein phosphatase-1 (PP1). Besides its neuronal roles, this protein also behaves as a key player in pathological and pharmacological aspects. Use of bioinformatics and phylogenetics approaches to further characterize the molecular features of DARPP-32 can guide future works. Predicted phosphorylation sites on DARPP-32 show conservation across vertebrates. Phylogenetics analysis indicates evolutionary strata of phosphorylation site acquisition at the C-terminus, suggesting functional expansion of DARPP-32, where more diverse signalling cues may involve in regulating DARPP-32 in inhibiting PP1 activity. Moreover, both phylogenetics and synteny analyses suggest de novo origination of PPP1R1 gene family via chromosomal rearrangement and exonization.
Southeast Asian Ovalocytosis (SAO) is a common red blood cell disorder that is maintained as a balanced polymorphism in human populations. In individuals heterozygous for the SAO-causing mutation there are minimal detrimental effects and well-documented protection from severe malaria caused by Plasmodium vivax and Plasmodium falciparum; however, the SAO-causing mutation is fully lethal in utero when homozygous. The present-day high frequency of SAO in Island Southeast Asia indicates the trait is maintained by strong heterozygote advantage. Our study elucidates the evolutionary origin of SAO by characterizing DNA sequence variation in a 9.5 kilobase region surrounding the causal mutation in the SLC4A1 gene. We find substantial haplotype diversity among SAO chromosomes and estimate the age of the trait to be approximately 10,005 years (95% CI: 4930-23,200 years). This date is far older than any other human malaria-resistance trait examined previously in Southeast Asia, and considerably pre-dates the widespread adoption of agriculture associated with the spread of speakers of Austronesian languages some 4000 years ago. Using a genealogy-based method we find no evidence of historical positive selection acting on SAO (s=0.0, 95% CI: 0.0-0.03), in sharp contrast to the strong present-day selection coefficient (e.g., 0.09) estimated from the frequency of this recessively lethal trait. This discrepancy may be due to a recent increase in malaria-driven selection pressure following the spread of agriculture, with SAO targeted as a standing variant by positive selection in malarial populations.
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.
Plant pathogens secrete effector proteins to manipulate the host and facilitate infection. Cognate hosts trigger strong defence responses upon detection of these effectors. Consequently, pathogens and hosts undergo rapid coevolutionary arms races driven by adaptive evolution of effectors and receptors. Because of their high rate of turnover, most effectors are thought to be species-specific and the evolutionary trajectories are poorly understood. Here, we investigate the necrosis-inducing protein 1 (NIP1) effector in the multihost pathogen genus Rhynchosporium. We retraced the evolutionary history of the NIP1 locus using whole-genome assemblies of 146 strains covering four closely related species. NIP1 orthologues were present in all species but the locus consistently segregated presence-absence polymorphisms suggesting long-term balancing selection. We also identified previously unknown paralogues of NIP1 that were shared among multiple species and showed substantial copy-number variation within R. commune. The NIP1A paralogue was under significant positive selection suggesting that NIP1A is the dominant effector variant coevolving with host immune receptors. Consistent with this prediction, we found that copy number variation at NIP1A had a stronger effect on virulence than NIP1B. Our analyses unravelled the origins and diversification mechanisms of a pathogen effector family shedding light on how pathogens gain adaptive genetic variation.
Weedy rice (Oryza spp.) is a major nuisance to rice farmers from all over the world. Although the emergence of weedy rice in East Malaysia on the island of Borneo is very recent, the threat to rice yield has reached an alarming stage. Using 47,027 genotyping-by-sequencing (GBS)-derived SNPs and candidate gene analysis of the plant architecture domestication gene TAC1, we assessed the genetic variations and evolutionary origin of weedy rice in East Malaysia. Our findings revealed two major evolutionary paths for genetically distinct weedy rice types. Whilst the cultivar-like weedy rice are very likely to be the weedy descendant of local coexisting cultivars, the wild-like weedy rice appeared to have arisen through two possible routes: (i) accidental introduction from Peninsular Malaysia weedy rice populations, and (ii) weedy descendants of coexisting cultivars. The outcome of our genetic analyses supports the notion that Sabah cultivars and Peninsular Malaysia weedy rice are the potential progenitors of Sabah weedy rice. Similar TAC1 haplotypes were shared between Malaysian cultivated and weedy rice populations, which further supported the findings of our GBS-SNP analyses. These different strains of weedy rice have convergently evolved shared traits, such as seeds shattering and open tillers. A comparison with our previous simple-sequence repeat-based population genetic analyses highlights the strength of genome-wide SNPs, including detection of admixtures and low-level introgression events. These findings could inform better strategic management for controlling the spread of weedy rice in the region.
Zanthoxylum armatum and Zanthoxylum bungeanum, known as 'Chinese pepper', are distinguished by their extraordinary complex genomes, phenotypic innovation of adaptive evolution and species-special metabolites. Here, we report reference-grade genomes of Z. armatum and Z. bungeanum. Using high coverage sequence data and comprehensive assembly strategies, we derived 66 pseudochromosomes comprising 33 homologous phased groups of two subgenomes, including autotetraploid Z. armatum. The genomic rearrangements and two whole-genome duplications created large (~4.5 Gb) complex genomes with a high ratio of repetitive sequences (>82%) and high chromosome number (2n = 4x = 132). Further analysis of the high-quality genomes shed lights on the genomic basis of involutional reproduction, allomones biosynthesis and adaptive evolution in Chinese pepper, revealing a high consistent relationship between genomic evolution, environmental factors and phenotypic innovation. Our study provides genomic resources and new insights for investigating diversification and phenotypic innovation in Chinese pepper, with broader implications for the protection of plants under severe environmental changes.
The mitochondrial genome sequence of the porcellanid crab, Petrolisthes haswelli is provided, making it the second for the family Porcellanidae and the third for the superfamily Galatheoidea. Petrolisthes haswelli has a mitogenome of 15,348 bp consisting of 13 protein-coding genes, two ribosomal subunit genes, 22 transfer RNAs and a non-coding AT-rich region. The base composition of the P. haswelli mitogenome is 35.66% for T, 18.65% for C, 34.35% for A and 11.34% for G, with an AT bias of 70.01%. The mitogenome gene order is identical to the mitogenome of Neopetrolisthes maculatus, the only other species of the family with a sequenced mitogenome.