Levels Of Nucleotide Diversity Research Articles

Species divergence and environmental adaptation of Picea asperata complex at the whole genome level.

To study the interspecific differentiation characteristics of species originating from recent radiation, the genotyping-by-sequencing (GBS) technique was used to explore the kinship, population structure, gene flow, genetic variability, genotype-environment association and selective sweeps of Picea asperata complex with similar phenotypes from a genome-wide perspective. The following results were obtained: 14 populations of P. asperata complex could be divided into 5 clades; P. wilsonii and P. neoveitchii diverged earlier and were more distantly related to the remaining 6 spruce species. Various geological events have promoted the species differentiation of P. asperata complex. There were four instances of gene flow among P. koraiensis, P. meyeri, P. asperata, P. crassifolia and P. mongolica. The population of P. mongolica had the highest level of nucleotide diversity, and P. neoveitchii may have experienced a bottleneck recently. Genotype-environment association found that a total of 20,808 genes were related to the environmental variables, which enhanced the adaptability of spruce in different environments. Genes that were selectively swept in the P. asperata complex were primarily associated with plant stress resistance. Among them were some genes involved in plant growth and development, heat stress, circadian rhythms and flowering. In addition to the commonly selected genes, different spruce species also displayed unique genes subjected to selective sweeps that improved their adaptability to different habitats. Understanding the interspecific gene flow and adaptive evolution of Picea species is beneficial to further understanding the species relationships of spruce and can provide a basis for studying spruce introgression and functional genomics.

The genome-wide signature of short-term temporal selection

Despite evolutionary biology's obsession with natural selection, few studies have evaluated multigenerational series of patterns of selection on a genome-wide scale in natural populations. Here, we report on a 10-y population-genomic survey of the microcrustacean Daphnia pulex. The genome sequences of [Formula: see text]800 isolates provide insights into patterns of selection that cannot be obtained from long-term molecular-evolution studies, including the following: the pervasiveness of near quasi-neutrality across the genome (mean net selection coefficients near zero, but with significant temporal variance about the mean, and little evidence of positive covariance of selection across time intervals); the preponderance of weak positive selection operating on minor alleles; and a genome-wide distribution of numerous small linkage islands of observable selection influencing levels of nucleotide diversity. These results suggest that interannual fluctuating selection is a major determinant of standing levels of variation in natural populations, challenge the conventional paradigm for interpreting patterns of nucleotide diversity and divergence, and motivate the need for the further development of theoretical expressions for the interpretation of population-genomic data.

Plasmodium vivax HAP2 genetic diversity and population structure from worldwide clinical samples. A potential Transmission-Blocking malaria vaccine candidate

BackgroundMalaria infections due to Plasmodium vivax are increasingly becoming a significant public health problem and are directly impacting malaria elimination efforts in many countries. Recent functional studies have shown that a pre-fertilising antigen; gamete membrane fusion protein HAP2/Generative Cell specific 1 (GCS1) domain significantly reduces the formation of P. vivax oocysts in mosquitos, implying that it could be developed as a transmission-blocking vaccine. MethodsIn this study, the genetic diversity, natural selection, haplotype network analysis and population structure of PvHAP2/GCS1 full-length gene using worldwide samples from 16 countries are determined and the implications for its use as a potential vaccine candidate against P. vivax are discussed. ResultsThe study revealed low levels of nucleotide diversity (π ∼ 0.003, SNPs = 95) and purifying selection for the full-length PvHAP2 gene indicating functional constraints. Among the 95 SNPs (65 synonymous and 30 non-synonymous), 15 non-synonymous substitutions were localised to the C-terminal domain of the protein (674a.a − 862a.a) majority of these non-synonymous substitutions were from two countries of Central Africa i.e. Cameroon and Gabon. Phylogeographic haplotype network analysis indicated that a major haplotype (Hap_2, n = 34) is shared among countries from Asian countries and South American countries. Population structure analysis using a robust Bayesian algorithm indicated the existence of two sub-populations of PvHAP2/GCS1 among the worldwide samples. ConclusionThis study suggests that, due to low genetic diversity across 15 countries, and shared haplotype identified, the HAP2 gene can be a good target for vaccine development as a transmission-blocking vaccine with the potential to provide strain-transcending immunity. The C-terminal region from two countries from Central Africa (Cameroon and Gabon) identified high number of non-synonymous substitutions, which may indicate new adaption of parasites from these countries. However, more detailed field based studies with a larger sample size from these countries are needed.

Population genetics of the critically endangered three-striped turtle, Batagur dhongoka, from the Ganga river system using mitochondrial DNA and microsatellite analysis

The three-striped roofed (Batagur dhongoka) is a semi-aquatic turtle that belongs to family Geoemydidae. Due to anthropogenic pressure, it has been facing an intense decline of over 80% in its distribution range in the past 50 years. It is considered as 'Critically Endangered' so effective conservation strategies are needed to protect the species by determining their genetic diversity and population genetic structure. This study investigates the genetic diversity, population structure and demographic pattern of B. dhongoka from two Turtle Rescue and Rehabilitation Centre established near Ganga river using mitochondrial cytochrome b (Cyt b: 1140 bp) ; control region (CR: 451 bp) and ten nuclear microsatellite loci. mtDNA results show low levels of nucleotide diversity (π = 0.0022) in B. dhongoka haplotypes and provide evidence for a low substitution rate. The demographic pattern estimated by the Bayesian skyline plot (BSP) analysis indicates historical stability followed by growth in the effective population size, with a recent reduction in population size from ~ 2 thousand years ago. The microsatellite findings show a moderate level of observed heterozygosity (Ho: 0.49). Bayesian-based clustering analysis revealed weak genetic structures in B. dhongoka and presence of admixed assignations suggesting close genetic relationships. These findings shed light on B. dhongoka's genetic status and underline the necessity of comprehensive rehabilitation and relocation programs and conservation and management techniques to ensure the species' long-term survival. In order to ensure the effective protection and conservation of B. dhongoka, the Government of India has taken a proactive measure by incorporating it into Schedule I of the Wildlife (Protection) Act, 1972, as amended in 2022.

New reference genomes to distinguish the sympatric malaria parasites, Plasmodium ovale curtisi and Plasmodium ovale wallikeri

Despite Plasmodium ovale curtisi (Poc) and wallikeri (Pow) being important human-infecting malaria parasites that are widespread across Africa and Asia, little is known about their genome diversity. Morphologically identical, Poc and Pow are indistinguishable and commonly misidentified. Recent rises in the incidence of Poc/Pow infections have renewed efforts to address fundamental knowledge gaps in their biology, and to develop diagnostic tools to understand their epidemiological dynamics and malaria burden. A major roadblock has been the incompleteness of available reference assemblies (PocGH01, PowCR01; ~ 33.5 Mbp). Here, we applied multiple sequencing platforms and advanced bioinformatics tools to generate new reference genomes, Poc221 (South Sudan; 36.0 Mbp) and Pow222 (Nigeria; 34.3 Mbp), with improved nuclear genome contiguity (> 4.2 Mbp), annotation and completeness (> 99% Plasmodium spp., single copy orthologs). Subsequent sequencing of 6 Poc and 15 Pow isolates from Africa revealed a total of 22,517 and 43,855 high-quality core genome SNPs, respectively. Genome-wide levels of nucleotide diversity were determined to be 2.98 × 10–4 (Poc) and 3.43 × 10–4 (Pow), comparable to estimates for other Plasmodium species. Overall, the new reference genomes provide a robust foundation for dissecting the biology of Poc/Pow, their population structure and evolution, and will contribute to uncovering the recombination barrier separating these species.

Demographic dynamics and molecular evolution of the rare and endangered subsect. Gerardianae of Pinus: insights from chloroplast genomes and mitochondrial DNA markers.

The divergence of subsect. Gerardianae was likely triggered by the uplift of the Qinghai-Tibetan Plateau and adjacent mountains. Pinus bungeana might have probably experienced expansion since Last Interglacial period. Historical geological and climatic oscillations have profoundly affected patterns of nucleotide variability, evolutionary history, and species divergence in numerous plants of the Northern Hemisphere. However, how long-lived conifers responded to geological and climatic fluctuations in East Asia remain poorly understood. Here, based on paternally inherited chloroplast genomes and maternally inherited mitochondrial DNA markers, we investigated the population demographic history and molecular evolution of subsect. Gerardianae (only including three species, Pinus bungeana, P. gerardiana, and P. squamata) of Pinus. A low level of nucleotide diversity was found in P. bungeana (π was 0.00016 in chloroplast DNA sequences, and 0.00304 in mitochondrial DNAs). The haplotype-based phylogenetic topology and unimodal distributions of demographic analysis suggested that P. bungeana probably originated in the southern Qinling Mountains and experienced rapid population expansion since Last Interglacial period. Phylogenetic analysis revealed that P. gerardiana and P. squamata had closer genetic relationship. The species divergence of subsect. Gerardianae occurred about 27.18 million years ago (Mya) during the middle to late Oligocene, which was significantly associated with the uplift of the Qinghai-Tibetan Plateau and adjacent mountains from the Eocene to the mid-Pliocene. The molecular evolutionary analysis showed that two chloroplast genes (psaI and ycf1) were under positive selection, the genetic lineages of P. bungeana exhibited higher transition and nonsynonymous mutations, which were involved with the strongly environmental adaptation. These findings shed light on the population evolutionary history of white pine species and provide striking insights for comprehension of their species divergence and molecular evolution.

High-Quality Genome Assemblies Reveal Evolutionary Dynamics of Repetitive DNA and Structural Rearrangements in the Drosophila virilis Subgroup.

High-quality genome assemblies across a range of nontraditional model organisms can accelerate the discovery of novel aspects of genome evolution. The Drosophila virilis group has several attributes that distinguish it from more highly studied species in the Drosophila genus, such as an unusual abundance of repetitive elements and extensive karyotype evolution, in addition to being an attractive model for speciation genetics. Here, we used long-read sequencing to assemble five genomes of three virilis group species and characterized sequence and structural divergence and repetitive DNA evolution. We find that our contiguous genome assemblies allow characterization of chromosomal arrangements with ease and can facilitate analysis of inversion breakpoints. We also leverage a small panel of resequenced strains to explore the genomic pattern of divergence and polymorphism in this species and show that known demographic histories largely predicts the extent of genome-wide segregating polymorphism. We further find that a neo-X chromosome in Drosophila americana displays X-like levels of nucleotide diversity. We also found that unusual repetitive elements were responsible for much of the divergence in genome composition among species. Helitron-derived tandem repeats tripled in abundance on the Y chromosome in D. americana compared to Drosophila novamexicana, accounting for most of the difference in repeat content between these sister species. Repeats with characteristics of both transposable elements and satellite DNAs expanded by 3-fold, mostly in euchromatin, in both D. americana and D. novamexicana compared to D. virilis. Our results represent a major advance in our understanding of genome biology in this emerging model clade.

History cooling events caused migrations of endangered Cathaya argyrophylla into low latitude refuges

Cathaya argyrophylla, referred as “the giant panda of plants”, is an endangered relict conifer species endemic to subtropical mountainous in China. The remanent populations of C. argyrophylla only survive in four refuges, while how the species migrated into these refuges remains unknown. To investigate the genetic structure and history of C. argyrophylla, here we genotyped four natural populations of C. argyrophylla using restriction site-related DNA sequencing (RAD-Seq). The four populations exhibited a relatively high level of nucleotide diversity (π = 0.29–0.33) and low differentiation (Fst = 0.04–0.12), with relatively high gene flow frequency (Nm = 1.86–6.00). We identified two historical migration events from high latitude populations to low latitude ones. The divergence times were corresponded to the starting time of population decline, which coincided with Naynayxungla Glaciation (∼0.8–0.5 Mya), the last glacial period (LGP, ∼115–11.7 kya) and the Holocene (∼10 kya). These results indicate that historical cooling events drove the endemic status of C. argyrphylla and its migration into low latitude refuges. Our findings provide insights into the evolutionary history and genetic diversity of this endangered species.

Population Genomics of the Facultatively Sexual Liverwort Marchantia polymorpha.

The population genomics of facultatively sexual organisms are understudied compared with their abundance across the tree of life. We explore patterns of genetic diversity in two subspecies of the facultatively sexual liverwort Marchantia polymorpha using samples from across Southern Ontario, Canada. Despite the ease with which M. polymorpha should be able to propagate asexually, we find no evidence of strictly clonal descent among our samples and little to no signal of isolation by distance. Patterns of identity-by-descent tract sharing further showed evidence of recent recombination and close relatedness between geographically distant isolates, suggesting long distance gene flow and at least a modest frequency of sexual reproduction. However, the M. polymorpha genome contains overall very low levels of nucleotide diversity and signs of inefficient selection evidenced by a relatively high fraction of segregating deleterious variants. We interpret these patterns as possible evidence of the action of linked selection and a small effective population size due to past generations of asexual propagation. Overall, the M. polymorpha genome harbors signals of a complex history of both sexual and asexual reproduction.

Genomic diversity and natural recombination of equid gammaherpesvirus 5 isolates

BackgroundEquid gammaherpesvirus 5 (EHV5) is closely related to equid gammaherpesvirus 2 (EHV2). Detection of EHV5 is frequent in horse populations worldwide, but it is often without a clear and significant clinical impact. Infection in horses can often present as subclinical disease; however, it has been associated with respiratory disease, including equine multinodular pulmonary fibrosis (EMPF). Genetic heterogeneity within small regions of the EHV5 glycoprotein B (gB) sequences have been reported and multiple genotypes of this virus have been identified within individual horses, but full genome sequence data for these viruses is limited. The primary focus of this study was to assess the genomic diversity and natural recombination among EHV5 isolates. ResultsThe genome size of EHV5 prototype strain and the five EHV5 isolates cultured for this study, including four isolates from the same horse, ranged from 181,929 to 183,428 base pairs (bp), with the sizes of terminal repeat regions varying from 0 to 10 bp. The nucleotide sequence identity between the six EHV5 genomes ranged from 95.5 to 99.1%, and the estimated average nucleotide diversity between isolates was 1%. Individual genes displayed varying levels of nucleotide diversity that ranged from 0 to 19%. The analysis of nonsynonymous substitution (Ka > 0.025) revealed high diversity in eight genes. Genome analysis using RDP4 and SplitsTree programs detected evidence of past recombination events between EHV5 isolates. ConclusionGenomic diversity and recombination hotspots were identified among EHV5 strains. Recombination can drive genetic diversity, particularly in viruses that have a low rate of nucleotide substitutions. Therefore, the results from this study suggest that recombination is an important contributing factor to EHV5 genomic diversity. The findings from this study provide additional insights into the genetic heterogeneity of the EHV5 genome.

Weaker selection on genes with treatment-specific expression consistent with a limit on plasticity evolution in Arabidopsis thaliana.

Differential gene expression between environments often underlies phenotypic plasticity. However, environment-specific expression patterns are hypothesized to relax selection on genes, and thus limit plasticity evolution. We collated over 27 terabases of RNA-sequencing data on Arabidopsis thaliana from over 300 peer-reviewed studies and 200 treatment conditions to investigate this hypothesis. Consistent with relaxed selection, genes with more treatment-specific expression have higher levels of nucleotide diversity and divergence at nonsynonymous sites but lack stronger signals of positive selection. This result persisted even after controlling for expression level, gene length, GC content, the tissue specificity of expression, and technical variation between studies. Overall, our investigation supports the existence of a hypothesized trade-off between the environment specificity of a gene's expression and the strength of selection on said gene in A. thaliana. Future studies should leverage multiple genome-scale datasets to tease apart the contributions of many variables in limiting plasticity evolution.

Assessing diversity and phylogeny of Indonesian breadfruit (Artocarpus spp.) using internal transcribed spacer (ITS) region and leaf morphology

BackgroundBreadfruit (Artocarpus spp.) is the main genus of Moraceae with multipurpose benefits, both ecologically and economically important, e.g., food ingredients, building materials, traditional medicine, and natural insecticides. However, most endemic Artocarpus have been threatened due to natural disasters and habitat degradation. The objective of our study was to determine the genetic diversity and relationships of endemic Artocarpus from South Borneo, Indonesia, using an internal transcribed spacer (ITS) region and leaf morphology. ResultsMorphologically, endemic Artocarpus endemic to South Borneo, Indonesia, has a different leaf shape, i.e., narrow-obovate to broad-elliptic, from simple to deeply dissected. Following the ITS region, this germplasm has a moderate level of nucleotide diversity (0.069). The phylogenetic analysis revealed Artocarpus into four (4) main clades, where the nearest is shown by the ‘Puyian’ (Artocarpus rigidus) and ‘Binturung’ (Artocarpus odoratissimus) at a coefficient divergence of 0.027, whereas the furthest by ‘Kulur’ (A. camansi) and ‘Tiwadak’ (A. integer) at a coefficient of 0.132. ConclusionIn brief, although an endemic Artocarpus of South Borneo, Indonesia, has a moderate level of nucleotide diversity, this germplasm also shows a unique phylogenetic relationship. Thus, this information is urgent in supporting the future Artocarpus breeding and preservation programs, mainly to save this germplasm from being threatened.

Demography and linked selection interact to shape the genomic landscape of codistributed woodpeckers during the Ice Age.

The influence of genetic drift on population dynamics during Pleistocene glacial cycles is well understood, but the role of selection in shaping patterns of genomic variation during these events is less explored. We resequenced whole genomes to investigate how demography and natural selection interact to generate the genomic landscapes of Downy and Hairy Woodpecker, species codistributed in previously glaciated North America. First, we explored the spatial and temporal patterns of genomic diversity produced by neutral evolution. Next, we tested (i) whether levels of nucleotide diversity along the genome are correlated with intrinsic genomic properties, such as recombination rate and gene density, and (ii) whether different demographic trajectories impacted the efficacy of selection. Our results revealed cycles of bottleneck and expansion, and genetic structure associated with glacial refugia. Nucleotide diversity varied widely along the genome, but this variation was highly correlated between the species, suggesting the presence of conserved genomic features. In both taxa, nucleotide diversity was positively correlated with recombination rate and negatively correlated with gene density, suggesting that linked selection played a role in reducing diversity. Despite strong fluctuations in effective population size, the maintenance of relatively large populations during glaciations may have facilitated selection. Under these conditions, we found evidence that the individual demographic trajectory of populations modulated linked selection, with purifying selection being more efficient in removing deleterious alleles in large populations. These results highlight that while genome-wide variation reflects the expected signature of demographic change during climatic perturbations, the interaction of multiple processes produces a predictable and highly heterogeneous genomic landscape.

Survey and Molecular Characterization of Echinococcus granulosus sensu stricto from Livestock and Humans in the Altai Region of Xinjiang, China

Cystic echinococcosis (CE), caused by the metacestode Echinococcus granulosus sensu stricto (s.s.), is an important zoonotic parasite, endemic in the Altai region of China. It is a serious human health risk and causes livestock losses. To evaluate the prevalence, genetic variation, and population structure of CE, 2898 sheep and 703 cattle were examined from October 2019 to mid-February 2020 in the Altai region (Altai, Habahe, Fuhai, and Buerjin). Sheep had an infection rate of 4.52% (131/2898) and cattle had an infection rate of 4.84% (34/703). In total, 180 cyst isolates were obtained, including 131 sheep, 34 cattle, and 15 from CE human patients. The cysts were investigated using mitochondrial cytochrome C oxidase subunit 1 (cox1). Polymerase Chain Reaction (PCR) results showed that, among the two genotypes of E. granulosus s.s., there were 22 different haplotypes (Haps). Phylogenetic analysis and parsimony network indicated that seventeen (77.27%) Haps belonged to the sheep strain (G1 genotype) and five Haps (22.73%) belonged to the buffalo strain (G3 genotype). Hap3 was the most common haplotype (65.00%, 112/180), which belongs to the G1 genotype. Hap18−Hap22 were found in human samples, indicating that sheep and cattle reservoirs of human CE. Molecular diversity indices revealed the high levels of haplotype diversity and relatively low levels of nucleotide diversity. Tajima’s D and Fu’s Fs tests displayed that the Altai population had a significant deviation from neutrality. Based on pairwise fixation index (Fst) values, a low level of genetic differentiation was found between the populations of E. granulosus s.s. isolated from different regions. The present survey findings represent an epidemiological survey of CE in the Altai region where there were two genotypes simultaneously and will provide more information on the genetic structure of E. granulosus s.s. within this region.

Whole-genome resequencing reveals genetic diversity and selection characteristics of dairy goat.

The dairy goat is one of the earliest dairy livestock species, which plays an important role in the economic development, especially for developing countries. With the development of agricultural civilization, dairy goats have been widely distributed across the world. However, few studies have been conducted on the specific characteristics of dairy goat. In this study, we collected the whole-genome data of 89 goat individuals by sequencing 48 goats and employing 41 publicly available goats, including five dairy goat breeds (Saanen, Nubian, Alpine, Toggenburg, and Guanzhong dairy goat; n = 24, 15, 11, 6, 6), and three goat breeds (Guishan goat, Longlin goat, Yunshang Black goat; n = 6, 15, 6). Through compared the genomes of dairy goat and non-dairy goat to analyze genetic diversity and selection characteristics of dairy goat. The results show that the eight goats could be divided into three subgroups of European, African, and Chinese indigenous goat populations, and we also found that Australian Nubian, Toggenburg, and Australian Alpine had the highest linkage disequilibrium, the lowest level of nucleotide diversity, and a higher inbreeding coefficient, indicating that they were strongly artificially selected. In addition, we identified several candidate genes related to the specificity of dairy goat, particularly genes associated with milk production traits (GHR, DGAT2, ELF5, GLYCAM1, ACSBG2, ACSS2), reproduction traits (TSHR, TSHB, PTGS2, ESR2), immunity traits (JAK1, POU2F2, LRRC66). Our results provide not only insights into the evolutionary history and breed characteristics of dairy goat, but also valuable information for the implementation and improvement of dairy goat cross breeding program.

The complete chloroplast genome sequence of Rhododendron fortunei: Structural comparative and phylogenetic analysis in the ericaceae family

Rhododendron fortunei (Ericaceae) possesses valuable horticultural and medicinal values. However, the genomic information on R. fortunei is very limited. In this study, the complete chloroplast genome (cp) of R. fortunei was assembled and annotated, SSR loci were characterised, comparative genomic analysis was carried out, and phylogenetic research was also performed. The results showed that the R. fortunei cp genome was of a typical quadripartite structure (200,997 bp). The lengths of the large single copy region (LSC), the inverted repeat regions (IR), and the small single copy region (SSC) were 109,151 bp, 2,604 bp, and 44,619 bp, respectively. A total of 147 unique genes were identified, including 99 protein-coding genes, 42 tRNA genes, and 6 rRNA genes, respectively. Leucine (11.51%) and cysteine (1.15%) were the highest and lowest representative amino acids, respectively. The total of 30 codons with obvious codon usage bias were all A/U-ending codons. Among the 77 simple sequence repeats, the majority were mononucleotide A/T repeats located in the intergenic spacer region. Five gene regions showed high levels of nucleotide diversity (Pi > 0.03). The comparative genome analysis revealed 7 hotspot intergenic regions (trnI-rpoB, trnTrpl16, rpoA-psbJ, rps7-rrn16, ndhI-rps16, rps16-rps19, and rrn16-trnI), showing great potential as molecular makers for species authentication. Expansion and contraction were detected in the IR region of the R. fortunei cp genome. In the phylogenetic tree, R. fortunei was closely related to R. platypodum. This research will be beneficial for evolutionary and genetic diversity studies of R. fortunei and related species among the Ericaceae family.

The demographic history of house mice (Mus musculus domesticus) in eastern North America.

The Western European house mouse (Mus musculus domesticus) is a widespread human commensal that has recently been introduced to North America. Its introduction to the Americas is thought to have resulted from the transatlantic movements of Europeans that began in the early 16th century. To study the details of this colonization history, we examine population structure, explore relevant demographic models, and infer the timing of divergence among house mouse populations in the eastern United States using published exome sequences from five North American populations and two European populations. For North American populations of house mice, levels of nucleotide variation were lower, and low-frequency alleles were less common than for European populations. These patterns provide evidence of a mild bottleneck associated with the movement of house mice into North America. Several analyses revealed that one North American population is genetically admixed, which indicates at least two source populations from Europe were independently introduced to eastern North America. Estimated divergence times between North American and German populations ranged between ∼1,000 and 7,000 years ago and overlapped with the estimated divergence time between populations from Germany and France. Demographic models comparing different North American populations revealed that these populations diverged from each other mostly within the last 500 years, consistent with the timing of the arrival of Western European settlers to North America. Together, these results support a recent introduction of Western European house mice to eastern North America, highlighting the effects of human migration and colonization on the spread of an invasive human commensal.

Genetic Diversity, Haplotype Relationships, and kdr Mutation of Malaria Anopheles Vectors in the Most Plasmodium knowlesi-Endemic Area of Thailand.

Plasmodium knowlesi, a malaria parasite that occurs naturally in long-tailed macaques, pig-tailed macaques, and banded leaf monkeys, is currently regarded as the fifth of the human malaria parasites. We aimed to investigate genetic diversity based on the cytochrome c oxidase subunit I (COI) gene, detect Plasmodium parasites, and screen for the voltage-gated sodium channel (VGSC)-mutation-mediated knockdown resistance (kdr) of Anopheles mosquitoes in Ranong province, which is the most P. knowlesi-endemic area in Thailand. One hundred and fourteen Anopheles females belonging to eight species, including An. baimaii (21.05%), An. minimus s.s. (20.17%), An. epiroticus (19.30%), An. jamesii (19.30%), An. maculatus s.s. (13.16%), An. barbirostris A3 (5.26%), An. sawadwongporni (0.88%), and An. aconitus (0.88%), were caught in three geographical regions of Ranong province. None of the Anopheles mosquitoes sampled in this study were infected with Plasmodium parasites. Based on the sequence analysis of COI sequences, An. epiroticus had the highest level of nucleotide diversity (0.012), followed by An. minimus (0.011). In contrast, An. maculatus (0.002) had the lowest level of nucleotide diversity. The Fu's Fs and Tajima's D values of the Anopheles species in Ranong were all negative, except the Tajima's D values of An. minimus (0.077). Screening of VGSC sequences showed no presence of the kdr mutation of Anopheles mosquitoes. Our results could be used to further select effective techniques for controlling Anopheles populations in Thailand's most P. knowlesi-endemic area.

Origin of domesticated Qingtian paddy-field carp and its genetic differentiation from wild common carp populations

Common carp is the first domesticated fish species in the world, which is the main fish candidate for being cocultured with rice in paddy fields. There are several strains of rice field-domesticated common carp in China, such as Qingtian (paddy-field) carp, Congjiang Carp, Wuyi Carp etc. Unfortunately, their geographical origins were unclear and the adaptability to rice field of these common carp varieties have not yet been well studied. Here we used 2b-RAD simplified genome sequencing to obtain 108,813 SNP markers from two paddy-field carp populations (Qingtian and Wuyi Mountain) and eight wild carp populations to explore the phylogenetic relationship and differences between the Qingtian paddy field carp and the other neighboring wild and domesticated carp populations. Compared with the wild carp populations, the levels of nucleotide diversity of the two paddy-field carp populations were lower. The neighbor-joining (NJ) tree, principal component analysis (PCA), and population structure analysis all showed three clusters of the carp populations: Qingtian paddy-field carps as one of the clusters, paddy-field carps from Wuyi Mountain areas the second and all the wild common carp populations the third. Qingtian paddy-field carps were clearly separated from other populations, which was consistent with the population differentiation index (Fst) analysis. All these analyses indicated that Qingtian paddy-field carps showed the closest relationship with the Oujiang wild carp population (OJ), and less close to all other populations, providing clear evidence of the geographical origin of Qingtian paddy-field carp from Oujiang River. With wild carps as the control group and paddy-field carps as the selection group, 180 genes with a selection signature were obtained. The gene ontology (GO) classification and enrichment analysis showed that the genes with a selection signature were enriched in the pathway related to signal transduction function. We can thus conclude that Qingtian paddy field carp of the Asia first GIAHS rice-fish coculture system were originated from the domestication from the local Oujiang wild common carp population, which showed a clear genetic adaptive evolution in genes with function in signal transduction. These findings will help us to further disclose the domestication effects on fish that leading us better understanding the adaptive mechanism of fish to the rice field environment.

Genetic population structure and demographic history of the endemic fish Paralichthys olivaceus of the Northwest Pacific Ocean.

The Northwest Pacific marginal waters comprising the South China Sea, East China Sea, Yellow Sea, and the Sea of Japan have unique geomorphic features. The Japanese flounder Paralichthys olivaceus, which is endemic to the Northwest Pacific, has high nutritional, economic, and ecological value. To allow the examination of the demographic history and population structure of the most common P. olivaceus species range over the five marginal seas (East China Sea, Yellow Sea, Bohai Sea, Northwest Pacific Ocean, and the Sea of Japan), the mitochondrial DNA control region of 91 individuals from six populations in China was sequenced. These sequences were combined with 233 sequences from four populations distributed in the Northwestern Pacific Ocean for analysis. Higher levels of nucleotide diversity (0.032 ± 0.016) and haplotype diversity (0.996 ± 0.001) were observed. The peripheral Fuqing population in the East China Sea had the relatively lowest genetic diversity and highest differentiation. Furthermore, when the results of the isolation by distance test, spatial analysis of molecular variation and geographic barrier analysis are also considered, there is a clear need to prioritize resource conservation and enhancement measures in this area. The phylogenetic trees, structure assignment test, and haplotypes network revealed no significant differences in the genealogical structure among ten populations. Mismatch distribution analysis, Bayesian skyline plots, and neutrality tests suggested that P. olivaceus experienced population expansion during the Pleistocene. Ocean currents and climate change play important roles in shaping the geographical distribution and genetic population structure of P. olivaceus.