The estimation of population differentiation with microsatellite markers

Genetic analysis of the circumsporozoite gene in Plasmodium falciparum isolates from Cameroon: Implications for efficacy and deployment of RTS,S/AS01 vaccine

Current understanding of genetic polymorphisms and natural selection in Plasmodium falciparum circumsporozoite (PfCSP), the leading malaria vaccine, is crucial for the development of next-generation vaccines, and such data is lacking in Africa. Blood samples were collected among Plasmodium-infected individuals living in four Cameroonian areas (Douala, Maroua, Mayo-Oulo, Pette). DNA samples were amplified using nested PCR protocols, sequenced, and BLASTed. Single nucleotide polymorphisms (SNPs) were analysed in each PfCSP region, and their impact on PfCSP function/structure was predicted in silico. The N-terminal region showed a limited polymorphism with four haplotypes, and three novel SNPs (N68Y, R87W, K93E) were found. Thirty-five haplotypes were identified in the central region, with several variants (e.g., NVNP and KANP). The C-terminal region was also highly diverse, with 25 haplotypes and eight novel SNPs (N290D, N308I, S312G, K317A, V344I, D356E, E357L, D359Y). Most polymorphic codon sites were mainly observed in the Th2R subregion in isolates from Douala and Pette. The codon site 321 was under episodic positive selection. One novel (E357L) and three known (K322I, G349D, D359Y) SNPs show an impact on function/structure. This study showed extensive genetic diversity with geographical patterns and evidence of the selection of Cameroonian PfCSP central and C-terminal regions.

Wolbachia-Based Emerging Strategies for Control of Vector-Transmitted Disease

Dengue fever is a mosquito-transmitted disease of great public health importance. Dengue lacks adequate vaccine protection and insecticide-based methods of mosquito control are proving increasingly ineffective. Here we review the emerging use of mosquitoes transinfected with the obligate intracellular bacterium Wolbachia pipientis for vector control. Wolbachia often induces cytoplasmic incompatibility in its mosquito hosts, resulting in infertile progeny between an infected male and an uninfected female. Wolbachia infection also suppresses the replication of pathogens in the mosquito, a process known as "pathogen blocking". Two strategies have emerged. The first one releases Wolbachia-carriers (both male and female) to replace the wild mosquito population, a process driven by cytoplasmic incompatibility and that becomes irreversible once a threshold is reached. This suppresses disease transmission mainly by pathogen blocking and frequently requires a single intervention. The second strategy floods the field population with an exclusively male population of Wolbachia-carrying mosquitoes to generate infertile hybrid progeny. In this case, transmission suppression depends largely on decreasing the population density of mosquitoes caused by infertility and requires continued mosquito release. The efficacy of both Wolbachia-based approaches has been conclusively demonstrated by randomized and non-randomized studies of deployments across the world. However, results conducted in one setting cannot be directly or easily extrapolated to other settings because dengue incidence is highly affected by the conditions into which the mosquitoes are released. Compared to traditional methods, Wolbachia-based approaches are much more environmentally friendly and can be effective in the medium/long term. On the flip side, they are much more complex and cost-intensive operations, requiring a substantial investment, infrastructure, trained personnel, coordination between agencies, and community engagement. Finally, we discuss recent evidence suggesting that transinfected Wolbachia in released mosquitoes has a moderate potential risk of spreading genes in the environment.

Spatio-temporal analysis of genetic diversity of merozoite surface protein-3 alpha in Myanmar Plasmodium vivax isolates

Myanmar aims to eliminate malaria by 2030. However, recent increase of malaria incidence is a great challenge to archive that goal. Increasing prevalence of Plasmodium vivax also hinders this endeavor. Monitoring genetic structure of the parasite is necessary to understand genetic nature and evolutionary aspect of P. vivax population in Myanmar. Partial fragment flanking blocks I and II of merozoite surface protein-3 alpha of P. vivax (pvmsp-3α) was amplified from P. vivax isolates collected in Pyin Oo Lwin, Mandalay Region, Myanmar in 2013-2015. Sequence analysis of pvmsp-3α was performed to determine genetic diversity and natural selection of this gene. Spatio-temporal genetic changes of pvmsp-3α in Myanmar P. vivax population were also investigated via comparative analysis of gene sequences obtained in this study and previously reported Myanmar pvmsp-3α sequences. Genetic diversity of Myanmar pvmsp-3α was detected in P. vivax isolates analyzed. Size polymorphisms in block I and amino acid changes and recombination events in block II were main factors contributing to the genetic diversity of pvmsp-3α. Comparative spatio-temporal analysis with previously reported Myanmar pvmsp-3α populations revealed the presence of genetic differences by population with moderate genetic differentiation between populations. Similar pattern of natural selection was also detected in Myanmar pvmsp-3α populations. These suggested that enough size of the P. vivax population sufficient to generate or maintain the genetic diversity remains in the population. Thus, continuous molecular surveillance of genetic structure of Myanmar P. vivax is necessary.

Across a phylogeographic break in the Qinling Mountains-Huaihe River Line: Quaternary evolutionary history of a medicinal and edible homologous plant (Allium macrostemon) in China

Biogeographic barriers to gene flow are central to studies of plant phylogeography. There are many physical and geographic barriers in China, but few studies have used molecular ecological evidence to investigate the natural geographic isolation barrier of the Qinling Mountains-Huaihe River Line (QHL). Allium macrostemon is a precious Chinese perennial herb belonging to the Amaryllidaceae family. It is used as a food and medicine, with a variety of health and healing properties. Five SSR markers, three chloroplast DNA (cpDNA) markers (psbA-trnH, rps16 and trnL-F), one nuclear ribosomal DNA (nrDNA) marker (ITS), and simplified genome GBS sequencing were used to analyse the genetic diversity and structure of A. macrostemon. Combining SSR, cpDNA, nrDNA ITS data and GBS analysis results, we divided A. macrostemon populations into northern and southern groups, with the southern group further divided into southwestern and central-southeastern groups. Niche simulation results reveal that the distribution area of A. macrostemon will reach its maximum in the future. These data indicate that the regional separation of A. macrostemon has been maintained by the combined influence of a geographical barrier and Quaternary climate, and that the back-and-forth fluctuations of QHL and Quaternary climate have played an important role in this process. QHL acts as a north-south dividing line in phylogeography and population genetic structure, promoting physical geographic isolation. This study provides a theoretical basis for the conservation, development, and utilization of A. macrostemon resources. It further provides a reference for understanding the systematic geographical pattern of the large-scale spatial distribution of plants in China and enriches our understanding of Quaternary plant evolution in areas with complex terrain.

Genetic structure of apical membrane antigen-1 in Plasmodium falciparum isolates from Pakistan

Plasmodium falciparum apical membrane antigen-1 (PfAMA-1) is a major candidate for the blood-stage malaria vaccine. Genetic polymorphisms of global pfama-1suggest that the genetic diversity of the gene can disturb effective vaccine development targeting this antigen. This study was conducted to explore the genetic diversity and gene structure of pfama-1 among P. falciparum isolates collected in the Khyber Pakhtunkhwa (KP) province of Pakistan. A total of 19 full-length pfama-1 sequences were obtained from KP-Pakistan P. falciparum isolates, and genetic polymorphism and natural selection were investigated. KP-Pakistan pfama-1 exhibited genetic diversity, wherein 58 amino acid changes were identified, most of which were located in ectodomains, and domains I, II, and III. The amino acid changes commonly found in the ectodomain of global pfama-1 were also detected in KP-Pakistan pfama-1. Interestingly, 13 novel amino acid changes not reported in the global population were identified in KP-Pakistan pfama-1. KP-Pakistan pfama-1 shared similar levels of genetic diversity with global pfama-1. Evidence of natural selection and recombination events were also detected in KP-Pakistan pfama-1.

Genetic Diversity and Natural Selection of Plasmodium vivax Merozoite Surface Protein 8 in Global Populations

Plasmodium vivax Merozoite Surface Protein 8 (PvMSP8) is a promising candidate target for the development of multi-component vaccines. Therefore, determining the genetic variation pattern of Pvmsp8 is essential in providing a reference for the rational design of the P. vivax malaria vaccines. This study delves into the genetic characteristics of the Pvmsp8 gene, specifically focusing on samples from the China-Myanmar border (CMB) region, and contrasts these findings with broader global patterns. The study uncovers that Pvmsp8 exhibits a notable level of conservation across different populations, with limited polymorphisms and relatively low nucleotide diversity (0.00023-0.00120). This conservation contrasts starkly with the high polymorphisms found in other P. vivax antigens such as Pvmsp1. A total of 25 haplotypes and 14 amino acid mutation sites were identified in the global populations, and all mutation sites were confined to non-functional regions. The study also notes that most CMB Pvmsp8 haplotypes are shared among Burmese, Cambodian, Thai, and Vietnamese populations, indicating less geographical variance, but differ notably from those found in Pacific island regions or the Panama. The findings underscore the importance of considering regional genetic diversity in P. vivax when developing targeted malaria vaccines. Non departure from neutral evolution were found by Tajima's D test, however, statistically significant differences were observed between the kn/ks rates. The study's findings are crucial in understanding the evolution and population structure of the Pvmsp8 gene, particularly during regional malaria elimination efforts. The highly conserved nature of Pvmsp8, combined with the lack of mutations in its functional domain, presents it as a promising candidate for developing a broad and effective P. vivax vaccine. This research thus lays a foundation for the rational development of multivalent malaria vaccines targeting this genetically stable antigen.

Plasmodium vivax populations in the western Greater Mekong Subregion evaluated using a genetic barcode

An improved understanding of the Plasmodium vivax populations in the Great Mekong Subregion (GMS) is needed to monitor the progress of malaria elimination. This study aimed to use a P. vivax single nucleotide polymorphism (SNP) barcode to evaluate the population dynamics and explore the gene flow among P. vivax parasite populations in the western GMS (China, Myanmar and Thailand). A total of 315 P. vivax patient samples collected in 2011 and 2018 from four regions of the western GMS were genotyped for 42 SNPs using the high-throughput MassARRAY SNP genotyping technology. Population genetic analysis was conducted to estimate the genetic diversity, effective population size, and population structure among the P. vivax populations. Overall, 291 samples were successfully genotyped at 39 SNPs. A significant difference was observed in the proportion of polyclonal infections among the five P. vivax populations (P = 0.0012, Pearson Chi-square test, χ2 = 18.1), with western Myanmar having the highest proportion (96.2%, 50/52) in 2018. Likewise, the average complexity of infection was also highest in western Myanmar (1.31) and lowest in northeast Myanmar (1.01) in 2018. The older samples from western China in 2011 had the highest pairwise nucleotide diversity (π, 0.388 ± 0.046), expected heterozygosity (He, 0.363 ± 0.02), and the largest effective population size. In comparison, in the neighboring northeast Myanmar, the more recent samples in 2018 showed the lowest values (π, 0.224 ± 0.036; He, 0.220 ± 0.026). Furthermore, the 2018 northeast Myanmar parasites showed high and moderate genetic differentiation from other populations with FST values of 0.162-0.252, whereas genetic differentiation among other populations was relatively low (FST ≤ 0.059). Principal component analysis, phylogeny, and STRUCTURE analysis showed that the P. vivax population in northeast Myanmar in 2018 substantially diverged from other populations. Although the 42 SNP barcode is a valuable tool for tracking parasite origins of worldwide parasite populations, a more extended barcode with additional SNPs is needed to distinguish the more related parasite populations in the western GMS.

Declining Genetic Polymorphisms of the C-terminus Merozoite Surface Protein-1 Amidst Increased Plasmodium knowlesi Transmission in Thailand

Recent reports from Thailand reveal a substantial surge in Plasmodium knowlesi cases over the past decade, with a more than eightfold increase in incidence by 2023 compared to 2018. This study investigates temporal changes in genetic polymorphism associated with the escalating transmission of P. knowlesi malaria in Thailand over the past two decades. Twenty-five P. knowlesi samples collected in 2018-2023 were sequenced for the 42-kDa region of pkmsp1 and compared with 24 samples collected in 2000-2009, focusing on nucleotide diversity, natural selection, recombination rate, and population differentiation. Seven unique haplotypes were identified in recent samples, compared to 15 in earlier samples. Nucleotide and haplotype diversities were lower in recent samples (π = 0.016, Hd = 0.817) than in earlier samples (π = 0.018, Hd = 0.942). Significantly higher synonymous substitution rates were observed in both sample sets (dS - dN = 2.77 and 2.43, p < 0.05), indicating purifying selection and reduced genetic diversity over time. Additionally, 8 out of 17 mutation points were located on B-cell epitopes, suggesting an adaptive response by the parasites to evade immune recognition. Population differentiation analysis using the fixation index (Fst) revealed high genetic differentiation between parasite populations in central and southern Thailand or Malaysia. Conversely, the relatively lower Fst value between southern Thailand and Malaysia suggests a closer genetic relationship, possibly reflecting historical gene flow. In conclusion, our findings highlight a decline in genetic diversity and evidence of purifying selection associated with the recently increased incidence of P. knowlesi malaria in Thailand. The minor genetic differentiation between P. knowlesi populations from southern Thailand and Malaysia suggests a shared recent ancestry of these parasites and underscores the need for coordinated efforts between the two countries for the elimination of P. knowlesi.

Molecular genetic identification of Wuzhishan ant chicken, a newly discovered resource in China

Introduction

The Wuzhishan ant (MY) chicken exhibits significant differences from other chicken breeds. However, the molecular genetic relationship between the MY breed and other chicken breeds has not been assessed.

Methods

Whole-genome resequencing was used to compare genetic diversity, nucleotide diversity, the fixation index, the linkage disequilibrium coefficient, and phylogenetic tree relationships between the MY breed and the Wenchang (WC), Danzhou (DZ), Bawangling (BW), and Longsheng Feng (LF) breeds.

Results

A total of 21,586,378 singlenucleotide polymorphisms and 4,253,341 insertions/deletions were screened out among the five breeds. The MY breed had the second highest genomic genetic diversity and nucleotide diversity and the lowest LD coefficient among the five breeds. Moreover, the phylogenetic tree analysis showed that individual birds of each breed clustered together with those of their respective breeds.

Discussion

Our data indicated that the MY breed is distinct from the other breeds and can be considered a new genetic resource.

The Variability of Phytophthora infestans Isolates Collected from Estonian Islands in the Baltic Sea

Knowledge of a pathogen's genetic variability and population structure is of great importance to effective disease management. In this study, 193 isolates of Phytophthora infestans collected from three Estonian islands were characterized over 3 years using simple sequence repeat (SSR) marker data complemented by information on their mating type and resistance to metalaxyl. In combination with SSR marker data from samples in the neighboring Pskov region of Northwest Russia, the impact of regional and landscape structure on the level of genetic exchange was also examined. Among the 111 P. infestans isolates from Estonian islands, 49 alleles were detected among 12 SSR loci, and 59 SSR multilocus genotypes were found, of which 64% were unique. The genetic variation was higher among years than that among islands, as revealed by the analysis of molecular variance. The frequency of metalaxyl-resistant isolates increased from 9% in 2012 to 30% in 2014, and metalaxyl resistance was most frequent among A1 isolates. The test for isolation by distance among the studied regions was not significant, and coupled with the absence of genetic differentiation, the result revealed gene flow and the absence of local adaptation. The data are consistent with a sexual population in which diversity is driven by an annual germination of soilborne oospores. The absence of shared genotypes over the years has important implications when it comes to the management of diseases. Such population diversity can make it difficult to predict the nature of the outbreak in the coming year as the genetic makeup is different for each year.

The Poultry Red Mite, Dermanyssus gallinae, travels far but not frequently, and takes up permanent residence on farms

Management of Dermanyssus gallinae, a cosmopolitan hematophagous mite responsible for damage in layer poultry farming, is hampered by a lack of knowledge of its spatio-temporal population dynamics. Previous studies have shown that the circulation of this pest between farms is of strictly anthropogenic origin, that a mitochondrial haplogroup has been expanding on European farms since the beginning of the 21st century and that its local population growth may be particularly rapid. To refine our understanding of how D. gallinae spreads within and among farms, we characterized the genetic structure of mite populations at different spatial scales and sought to identify the main factors interrupting gene flow between poultry houses and between mitochondrial haplogroups. To this end, we selected and validated the first set of nuclear microsatellite markers for D. gallinae and sequenced a region of the CO1-encoding mitochondrial gene in a subsample of microsatellite-genotyped mites. We also tested certain conditions required for effective contamination of a poultry house through field experimentation, and conducted a survey of practices during poultry transfers. Our results confirm the role of poultry transport in the dissemination of mite populations, but the frequency of effective contamination after the introduction of contaminated material into poultry houses seems lower than expected. The high persistence of mites on farms, even during periods when poultry houses are empty and cleaned, and the very large number of nodes in the logistic network (large number of companies supplying pullets or transporting animals) undoubtedly explain the very high prevalence on farms. Substantial genetic diversity was measured in farm populations, probably as a result of the mite's known haplodiploid mode of sexual reproduction, coupled with the dense logistic network. The possibility of the occasional occurrence of asexual reproduction in this sexually reproducing mite was also revealed in our analyses, which could explain the extreme aggressiveness of its demographic dynamics under certain conditions.

Population genetic structure of Pomacea canaliculata in China based on the COI and ITS1 genes

Comprehending the phylogeography of invasive organisms enhances our insight into their distribution dynamics, which is instrumental for the development of effective prevention and management strategies. In China, Pomacea canaliculata and Pomacea maculata are the two most widespread and damaging species of the non-native Pomacea spp.. Given this species' rapid spread throughout country, it is urgent to investigate the genetic diversity and structure of its different geographic populations, a task undertaken in the current study using the COI and ITS1 mitochondrial and ribosomal DNA genes, respectively. The result of this study, based on a nationwide systematic survey, a collection of Pomacea spp., and the identification of cryptic species, showed that there is a degree of genetic diversity and differentiation in P. canaliculata, and that all of its variations are mainly due to differences between individuals within different geographical populations. Indeed, this species contains multiple haplotypes, but none of them form a systematic geographical population structure. Furthermore, the COI gene exhibits higher genetic diversity than the ITS1 gene. Our study further clarifies the invasive pathways and dispersal patterns of P. canaliculata in China to provide a theoretical basis.

A method for noninvasive individual genotyping of black-footed cat (Felis nigripes)

The black-footed cat (Felis nigripes) is endemic to the arid regions of southern Africa. One of the world's smallest wild felids, the species occurs at low densities and is secretive and elusive, which makes ecological studies difficult. Genetic data could provide key information such as estimates on population size, sex ratios, and genetic diversity. In this study, we test if microsatellite loci can be successfully amplified from scat samples that could be noninvasively collected from the field. Using 21 blood and scat samples collected from the same individuals, we statistically tested whether nine microsatellites previously designed for use in domestic cats can be used to identify individual black-footed cats. Genotypes recovered from blood and scat samples were compared to assess loss of heterozygosity, allele dropout, and false alleles resulting from DNA degradation or PCR inhibitors present in scat samples. The microsatellite markers were also used to identify individuals from scats collected in the field that were not linked to any blood samples. All nine microsatellites used in this study were amplified successfully and were polymorphic. Microsatellite loci were found to have sufficient discriminatory power to distinguish individuals and identify clones. In conclusion, these molecular markers can be used to monitor populations of wild black-footed cats noninvasively. The genetic data will be able to contribute important information that may be used to guide future conservation initiatives.

Genetic differentiation of Plasmodium vivax duffy binding protein in Ethiopia and comparison with other geographical isolates

BACKGROUND

Plasmodium vivax Duffy binding protein (PvDBP) is a merozoite surface protein located in the micronemes of P. vivax. The invasion of human reticulocytes by P. vivax merozoites depends on the parasite DBP binding domain engaging Duffy Antigen Receptor for Chemokine (DARC) on these red blood cells (RBCs). PvDBPII shows high genetic diversity which is a major challenge to its use in the development of a vaccine against vivax malaria.

METHODS

A cross-sectional study was conducted from February 2021 to September 2022 in five study sites across Ethiopia. A total of 58 blood samples confirmed positive for P. vivax by polymerase chain reaction (PCR) were included in the study to determine PvDBPII genetic diversity. PvDBPII were amplified using primers designed from reference sequence of P. vivax Sal I strain. Assembling of sequences was done using Geneious Prime version 2023.2.1. Alignment and phylogenetic tree constructions using MEGA version 10.1.1. Nucleotide diversity and haplotype diversity were analysed using DnaSP version 6.12.03, and haplotype network was generated with PopART version 1.7.

RESULTS

The mean age of the participants was 25 years, 5 (8.6%) participants were Duffy negatives. From the 58 PvDBPII sequences, seven haplotypes based on nucleotide differences at 8 positions were identified. Nucleotide diversity and haplotype diversity were 0.00267 ± 0.00023 and 0.731 ± 0.036, respectively. Among the five study sites, the highest numbers of haplotypes were identified in Arbaminch with six different haplotypes while only two haplotypes were identified in Gambella. The phylogenetic tree based on PvDBPII revealed that parasites of different study sites shared similar genetic clusters with few exceptions. Globally, a total of 39 haplotypes were identified from 223 PvDBPII sequences representing different geographical isolates obtained from NCBI archive. The nucleotide and haplotype diversity were 0.00373 and 0.845 ± 0.015, respectively. The haplotype prevalence ranged from 0.45% to 27.3%. Two haplotypes were shared among isolates from all geographical areas of the globe.

CONCLUSIONS

PvDBPII of the Ethiopian P. vivax isolates showed low nucleotide but high haplotype diversity, this pattern of genetic variability suggests that the population may have undergone a recent expansion. Among the Ethiopian P. vivax isolates, almost half of the sequences were identical to the Sal-I reference sequence. However, there were unique haplotypes observed in the Ethiopian isolates, which does not share with isolates from other geographical areas. There were two haplotypes that were common among populations across the globe. Categorizing population haplotype frequency can help to determine common haplotypes for designing an effective blood-stage vaccine which will have a significant role for the control and elimination of P. vivax.

Study on the characteristics of genetic diversity of different populations of Guizhou endemic plant Rhododendron pudingense based on microsatellite markers

BACKGROUND

Rhododendron pudingense, firstly discovered in Puding county of Guizhou province in 2020, have adapted to living in rocky fissure habitat, which has important ornamental and economic values. However, the genetic diversity and population structure of this species have been rarely described, which seriously affects the collection and protection of wild germplasm resources.

RESULTS

In the present study, 13 pairs of primers for polymorphic microsatellite were used to investigate the genetic diversity of 65 R. pudingense accessions from six different geographic populations. A total of 254 alleles (Na) were obtained with an average of 19.5 alleles per locus. The average values of polymorphic information content (PIC), observed heterozygosity (Ho), and expected heterozygosity (He) were 0.8826, 0.4501, and 0.8993, respectively, These results indicate that the microsatellite primers adopted demonstrate good polymorphism, and the R. pudingense exhibits a high level of genetic diversity at the species level. The average genetic differentiation coefficient (Fst) was 0.1325, suggested that moderate divergence occurred in R. pudingense populations. The average values of genetic differentiation coefficient and gene flow among populations were 0.1165 and 3.1281, respectively. The analysis of molecular variance (AMOVA) indicated that most of the population differences (88%) were attributed to within-population variation. The PCoA results are consistent with the findings of the UPGMA clustering analysis, supporting the conclusion that the six populations of R. pudingense can be clearly grouped into two separate clusters. Based on Mantel analysis, we speculate that the PD population may have migrated from WM-1 and WM-2. Therefore, it is advised to protect the natural habitat of R. pudingense in situ as much as possible, in order to maximize the preservation of its genetic diversity.

CONCLUSIONS

This is the first comprehensive analysis of genetic diversity and population structure of R. pudingense in Guizhou province. The research results revealed the high genetic diversity and moderate population diferentiation in this horticulture plant. This study provide a theoretical basis for the conservation of wild resources of the R. pudingense and lay the foundation for the breeding or cultivation of this new species.

Morphometric characteristics and microsatellite markers based diversity and differentiation recognizes the first prospective cattle breed from the Jharkhand state of India

India is bestowed with immense cattle biodiversity with 50 registered breeds. However, the majority (59.3%) is yet not characterized. Identification and characterization are the gateways to the management of prized indigenous resources. Present research described a unique cattle population of Jharkhand state, managed under a traditional low-input, low-output system. It was characterized by morphological traits, performance parameters, and management practices. Animals have the characteristic pre-scapular location of the hump. Genetic variation within this population and its differentiation with the six closely distributed cattle breeds were evaluated using FAO recommended microsatellite markers. Jharkhandi cattle have substantial genetic variation based on gene diversity (>0.6) and the average number of alleles per locus (>8). The population did not suffer from a genetic bottleneck in the recent past. Pairwise Nei's genetic distance, phylogenetic relationship, population differentiation, and the correct assignment of all the animals to self group substantiated its separate genetic identity. Since gene flow (Nm = 2.8-7.32) was identified and admixture was indicated by the Bayesian analysis there is a pressing need for scientific management of this population. Results endow authorities with critical information for registering a new Indian cattle breed (Medini) that contributes to the food security, livelihood, and economic sustainability of rural tribal households.

Genetic diversity and differentiation of Thutho cattle from northeast India using microsatellite markers

Cattle are losing maximum breeds among the world's livestock. Genetic variability data is essentially required for conservation decision-making. Thutho is a recently registered Indian cattle breed (INDIA_CATTLE_1400_THUTHO_03047) from the northeast region (NE), a biodiversity hotspot. Genetic diversity in the Thutho population and its differentiation from the only other cattle breed of NE (Siri) and cattle (Bachaur) of the neighboring region was established using highly polymorphic, FAO-recommended microsatellite markers. Numerous alleles (253) were detected across the 25 loci. The mean observed and expected numbers of alleles in the population were 10.12 ± 0.5 and 4.5 ± 0.37, respectively. The observed heterozygosity (0.67 ± 0.04) was lower than the expected heterozygosity (0.73 ± 0.03) which indicated a departure from the Hardy-Weinberg equilibrium. A positive FIS value (0.097) confirmed the heterozygote deficiency in the Thutho population. Genetic distance, phylogenetic relationships, differentiation parameters, population assignment, and Bayesian analysis explicitly ascertained the unique genetic identity of the Thutho cattle. The population did not suffer any bottlenecks in the past. Thutho has minimum diversity among the three populations; hence, its scientific management needs to be initiated immediately. Interestingly, genetic variation is enough for formulating breeding programs for managing, improving, and conserving this precious indigenous cattle germplasm.

Ancestral Information Analysis of Chinese Korean Ethnic Group via a Novel Multiplex DIP System

Deletion/insertion polymorphism (DIP) is one of the more promising genetic markers in the field of forensic genetics for personal identification and biogeographic ancestry inference. In this research, we used an in-house developed ancestry-informative marker-DIP system, including 56 autosomal diallelic DIPs, three Y-chromosomal DIPs, and an Amelogenin gene, to analyze the genetic polymorphism and ancestral composition of the Chinese Korean group, as well as to explore its genetic relationships with the 26 reference populations. The results showed that this novel panel exhibited high genetic polymorphism in the studied Korean group and could be effectively applied for forensic individual identification in the Korean group. In addition, the results of multiple population genetic analyses indicated that the ancestral component of the Korean group was dominated by northern East Asia. Moreover, the Korean group was more closely related to the East Asian populations, especially to the Japanese population in Tokyo. This study enriched the genetic data of the Korean ethnic group in China and provided information on the ancestry of the Korean group from the perspective of population genetics.

Plasmodium falciparum population structure and genetic diversity of cell traversal protein for ookinetes and sporozoites (CelTOS) during malaria resurgences in Dielmo, Senegal

The ability to accurately measure the intensity of malaria transmission in areas with low transmission is extremely important to guide elimination efforts. Plasmodium falciparum Cell-traversal protein for ookinetes and sporozoites (PfCelTOS) is an important conserved sporozoite antigen reported as one of the promising malaria vaccine candidates, and could be used to estimate malaria transmission intensity. This study aimed at determining whether the diversity of PfCelTOS gene reflects the changes in malaria transmission that occurred between 2007 and 2014 in Dielmo, a Senegalese village, before and after the implementation of insecticide treated bed nets (ITNs). Of the 109 samples positive for PfCelTOS PCR, 96 (88%) were successfully sequenced and analysed for polymorphisms and population diversity. The number of segregating sites was higher during the pre-intervention period (13) and the malaria resurgences (11) than during the intervention period (5). Similarly, the number and diversity of haplotypes were higher during the pre-intervention period (16 and 0.914, respectively) and the malaria resurgences (6 and 0.821, respectively) than during the intervention period (4 and 0.758, respectively). Moreover, the average number of nucleotide differences was higher during the pre-intervention (3.792) and during malaria resurgences (3.467) than during the intervention period (2.189). The 3D7 KSSFNEP haplotype was only observed during the intervention period. Only two haplotypes were shared in both the pre-intervention and intervention periods while four haplotypes were shared between the pre-intervention and the malaria resurgences. The Fst values indicate moderate differentiation between pre-intervention and intervention periods (0.17433), and between intervention and malaria resurgences period (0.19198) as well as between pre-intervention and malaria resurgences periods (0.06607). PfCelTOS genetic diversity reflected changes of malaria transmission, with higher polymorphisms recorded before the large-scale implementation of ITNs and during the malaria resurgences. PfCelTOS is also a candidate vaccine; mapping its diversity across multiple endemic environments will facilitate the design and optimisation of a broad and efficacious vaccine.

Genetic characteristics of P. falciparum parasites collected from 2012 to 2016 and anti-malaria resistance along the China-Myanmar border

BACKGROUNDS

The therapeutic efficacy studies of DHA-PIP for uncomplicated Plasmodium falciparum patients were implemented from 2012 to 2016 along China (Yunnan province)-Myanmar border, which verified the high efficacy of DHA-PIP. With the samples collected in these studies, the genetic characteristics of P. falciparum parasites based on in vivo parasite clearance time (PCT) was investigated to explore if these parasites had developed resistance to DHA and PIP at molecular level.

METHODS

The genetic characteristics were investigated based on K13 genotypes, copy numbers of genes pfpm2 and pfmdr1, and nine microsatellite loci (Short Tandem Repeats, STR) flanking the K13 gene on chromosome 13. The PCT 50s were compared based on different K13 genotypes, sites, periods and copy numbers.

RESULTS

In the NW (North-West Yunnan province bordering with Myanmar) region, F446I was the main K13 genotype. No significant differences for PCT 50s presented among three K13 genotypes. In SW (South-West Yunnan province bordering with Myanmar) region, only wild K13 genotype was detected in all parasite isolates whose PCT 50s was significantly longer than those in NW region. For the copy numbers of genes, parasite isolates containing multiple copies of pfmdr1 gene were found in both regions, but only single copy of pfpm2 gene was detected. Though the prevalence of parasite isolates with multiple copies of pfmdr1 gene in SW region was higher than that in NW region, no difference in PCT 50s were presented between isolates with single and multiple copies of pfmdr1 gene. The median He values of F446I group and Others (Non-F446I K13 mutation) group were 0.08 and 0.41 respectively. The mean He values of ML group (Menglian County in SW) and W (wild K13 genotype in NW) group were 0 and 0.69 respectively. The mean Fst values between ML and W groups were significantly higher than the other two K13 groups.

CONCLUSIONS

P. falciparum isolates in NW and SW regions had very different genetic characteristics. The F446I was hypothesized to have independently appeared and spread in NW region from 2012 and 2016. The high susceptibility of PIP had ensured the efficacy of DHA-PIP in vivo. Multiple copy numbers of pfmdr1 gene might be a potential cause of prolonged clearance time of ACTs drugs along China-Myanmar border.

TRIAL REGISTRATION

Trial registration: ISRCTN, ISRCTN 11775446. Registered 17 April 2020-Retrospectively registered, the registered name was Investigating resistance to DHA-PIP for the treatment of Plasmodium falciparum malaria and chloroquine for the treatment of Plasmodium vivax malaria in Yunnan, China. http://www.isrctn.com/ISRCTN11775446.

Introduction and effect of natural selection analysis at common mutations of SARS-CoV-2 spike gene in Iran

The epidemic of coronavirus disease 2019 (COVID-19) was caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The spike (S) protein of SARS-Cov-2 is composed of two subunits, S1 and S2. This study aimed to describe SARS-CoV-2 haplotypes in Iranians based on the S gene, which plays a key role in the receptor recognition and cell membrane fusion proses. 95 positive saliva samples for SARS-CoV-2 were amplified and sequenced for the S gene. The sequences were classified into 35 haplotypes, which 11 haplotypes were new (H1, H2, H3, H4, H6, H7, H11, H13, H15, H16, H25) and have not been reported so far. Amino acid substitutions were found at 40 positions that 23 were located at S1 subunit and 16 were at S2 subunit and one was at cleavage loop (P681H/R), thus polymorphisms at S1 subunit were found to be higher than S2. The neutrality index (NI) analyses showed a negative departure from the neutral substitution patterns (NI > 1) for S1 and S2 subunit in the studied sequences. The co-occurrence of B-cell epitopes and mutation sites were found in seven positions with more probably to be exposed the immune system pressure. In conclusion, the results provide the significant data to design an effective vaccine based on this protein.

Species Delimitation and Genetic Relationship of Castanopsis hainanensis and Castanopsis wenchangensis (Fagaceae)

Castanopsis is one of the most common genus of trees in subtropical evergreen broad-leaved forests and tropical monsoon rainforests in China. Castanopsis hainanensis and Castanopsis wenchangensis are endemic to Hainan Island, but they were once confused as the same species due to very similar morphologies. In this study, nuclear microsatellite markers and chloroplast genomes were used to delimit C. hainanensis and C. wenchangensis. The allelic variations of nuclear microsatellites revealed that C. hainanensis and C. wenchangensis were highly genetically differentiated with very limited gene admixture. Both showed higher genetic diversity within populations and lower genetic diversity among populations, and neither had further population genetic structure. Furthermore, C. wenchangensis and C. hainanensis had very different chloroplast genomes. The independent genetic units, very limited gene admixture, different distribution ranges, and distinct habitats all suggest that C. wenchangensis and C. hainanensis are independent species, thus they should be treated as distinct conservation units.

Genetic diversity and molecular evolution of Plasmodium vivax Duffy Binding Protein and Merozoite Surface Protein-1 in northwestern Thailand

The local diversity and population structure of malaria parasites vary across different regions of the world, reflecting variations in transmission intensity, host immunity, and vector species. This study aimed to use amplicon sequencing to investigate the genotypic patterns and population structure of P. vivax isolates from a highly endemic province of Thailand in recent years. Amplicon deep sequencing was performed on 70 samples for the 42-kDa region of pvmsp1 and domain II of pvdbp. Unique haplotypes were identified and a network constructed to illustrate genetic relatedness in northwestern Thailand. Based on this dataset of 70 samples collected between 2015 and 2021, 16 and 40 unique haplotypes were identified in pvdbpII and pvmsp142kDa, respectively. Nucleotide diversity was higher in pvmsp142kDa than in pvdbpII (π = 0.027 and 0.012), as was haplotype diversity (Hd = 0.962 and 0.849). pvmsp142kDa also showed a higher recombination rate and higher levels of genetic differentiation (Fst) in northwestern Thailand versus other regions (0.2761-0.4881). These data together suggested that the genetic diversity of P. vivax in northwestern Thailand at these two studied loci evolved under a balancing selection, most likely host immunity. The lower genetic diversity of pvdbpII may reflect its stronger functional constrain. In addition, despite the balancing selection, a decrease in genetic diversity was observed. Hd of pvdbpII decreased from 0.874 in 2015-2016 to 0.778 in 2018-2021; π of pvmsp142kDa decreased from 0.030 to 0.022 over the same period. Thus, the control activities must have had a strong impact on the parasite population size. The findings from this study provide an understanding of P. vivax population structure and the evolutionary force on vaccine candidates. They also established a new baseline for tracking future changes in P. vivax diversity in the most malarious area of Thailand.

Associations between genetic variants of HSD17B13 and fasting plasma glucose in Chinese children

BACKGROUND AND AIMS

Genetic variants in 17-β hydroxysteroid dehydrogenase 13 (HSD17B13) were demonstrated to protect against NAFLD, which is highly related with insulin resistance and dyslipidemia. However, the effects of NAFLD associated HSD17B13 variants on circulating glucose and lipids have not been adequately investigated in children. This study aimed to investigate associations between single nucleotide polymorphisms (SNPs) of HSD17B13 and NAFLD or its related phenotypes, such as blood glucose and serum lipids in Chinese children.

METHODS AND RESULTS

We studied 1027 Chinese Han children aged 7-18 years old, which included 162 NAFLD children and 865 controls without NAFLD. Three SNPs (rs13112695, rs7692397, rs6834314) in HSD17B13 were genotyped. The multivariable logistic and linear regression models were applied to detect the associations between three SNPs and NAFLD or its related phenotypes [alanine transaminase (ALT), fasting plasma glucose (FPG) and serum lipids]. The effect allele A of rs7692397 was negatively associated with FPG [β (SE) = -0.088 (0.027) mmol/L, P = 0.001], whereas the effect allele G of rs6834314 was positively associated with FPG (β (SE) = 0.060 (0.019) mmol/L, P = 0.002). After Bonferroni correction, the significant associations still remained (both P < 0.0024). No significant associations were found for NAFLD or serum lipids.

CONCLUSION

The study firstly revealed the association between two HSD17B13 variants and FPG in Chinese children, providing evidence for HSD17B13 variants and abnormal glucose metabolism.

Genetic diversity and conservation of Siberian apricot (Prunus sibirica L.) based on microsatellite markers

Siberian apricot (Prunus sibirica L.) is a woody tree species of ecological, economic, and social importance. To evaluate the genetic diversity, differentiation, and structure of P. sibirica, we analyzed 176 individuals from 10 natural populations using 14 microsatellite markers. These markers generated 194 alleles in total. The mean number of alleles (13.8571) was higher than the mean number of effective alleles (6.4822). The average expected heterozygosity (0.8292) was higher than the average observed heterozygosity (0.3178). Shannon information index and polymorphism information content were separately 2.0610 and 0.8093, demonstrating the rich genetic diversity of P. sibirica. Analysis of molecular variance revealed that 85% of the genetic variation occurred within populations, with only 15% among them. The genetic differentiation coefficient and gene flow were separately 0.151 and 1.401, indicating a high degree of genetic differentiation. Clustering results showed that a genetic distance coefficient of 0.6 divided the 10 natural populations into two subgroups (subgroups A and B). STRUCTURE and principal coordinate analysis divided the 176 individuals into two subgroups (clusters 1 and 2). Mantel tests revealed that genetic distance was correlated with geographical distance and elevation differences. These findings can contribute to the effective conservation and management of P. sibirica resources.

Genome-wide scans for selection signatures in indigenous pigs revealed candidate genes relating to heat tolerance

Heat stress is a major problem that constrains pig productivity. Understanding and identifying adaptation to heat stress has been the focus of recent studies, and the identification of genome-wide selection signatures can provide insights into the mechanisms of environmental adaptation. Here, we generated whole-genome re-sequencing data from six Chinese indigenous pig populations to identify genomic regions with selection signatures related to heat tolerance using multiple methods: three methods for intra-population analyses (Integrated Haplotype Score, Runs of Homozygosity and Nucleotide diversity Analysis) and three methods for inter-population analyses (Fixation index (F_ST), Cross-population Composite Likelihood Ratio and Cross-population Extended Haplotype Homozygosity). In total, 1 966 796 single nucleotide polymorphisms were identified in this study. Genetic structure analyses and F_ST indicated differentiation among these breeds. Based on information on the location environment, the six breeds were divided into heat and cold groups. By combining two or more approaches for selection signatures, outlier signals in overlapping regions were identified as candidate selection regions. A total of 163 candidate genes were identified, of which, 29 were associated with heat stress injury and anti-inflammatory effects. These candidate genes were further associated with 78 Gene Ontology functional terms and 30 Kyoto Encyclopedia of Genes and Genomes pathways in enrichment analysis (P < 0.05). Some of these have clear relevance to heat resistance, such as the AMPK signalling pathway and the mTOR signalling pathway. The results improve our understanding of the selection mechanisms responsible for heat resistance in pigs and provide new insights of introgression in heat adaptation.

The Effect of Landscape Environmental Factors on Gene Flow of Red Deer (Cervus canadensis xanthopygus) in the Southern of the Greater Khingan Mountains, China

Red deer (Cervus canadensis xanthopygus) living in the north of China are restricted and threatened due to human activities and the changes in the natural environment, which influence the dispersal and effective gene flow between different groups of red deer. Effective gene flow plays an important role in maintaining genetic diversity and structure and ensuring population health. In order to evaluate the genetic diversity level and understand the gene flow between different red deer groups, 231 fresh fecal samples were collected from the southern part of the Greater Khingan Mountains, China. A microsatellite marker was used for genetic analysis. The results showed that the genetic diversity of red deer was intermediate in this region. Significant genetic differentiation among different groups was found in the main distribution area (p < 0.01) using F-statistics and the program STRUCTURE. Different degrees of gene flow existed in red deer groups, and the roads (importance = 40.9), elevation (importance = 38.6), and settlements (importance = 14.1) exerted main effects on gene flow between red deer groups. Human-made factors should be noticed and strictly supervised in this region to avoid excessive disturbance to the normal movement of the red deer. Further conservation and management of red deer should reduce the intensity of vehicular traffic in the concentrated distribution areas of red deer, especially during the heat season. This research helps us better understand the genetic level and health status of red deer in the southern part of the Greater Khingan Mountains and provides theoretical references for protecting and restoring the red deer populations in China.

Aspidistradaibuensis var. longkiauensis, a new variety of Aspidistra (Asparagaceae) from Taiwan, identified through morphological and genetic analyses

Aspidistra Ker Gawl. is one of the the most diverse and fastest-growing genera of angiosperm. Most Aspidistra species have been discovered in a limited area or a single site through morphological comparison. Because of the lack of population studies, morphological variation within species and the boundaries of some species remain unclear. In recent years, combining genetic and morphological markers has become a powerful approach for species delimitation. In this study, we performed population sampling and integrated morphometrics and microsatellite genetic diversity analyses to determine the species diversity of Aspidistra in Taiwan. We identified three species, namely Aspidistraattenuata Hayata; A.daibuensisHayatavar.daibuensis; A.mushaensisHayatavar.mushaensis; and reduced A.longiconnectiva C.T.Lu, K.C.Chuang & J.C.Wang to the variety level, and described a new variety, A.daibuensisHayatavar.longkiauensis. The description, diagnosis, distribution, and photographs of this new variety as well as a key to the known Taiwanese Aspidistra are provided.

Hawkmoth and bee pollinators impact pollen dispersal at the landscape but not local scales in two species of Oenothera

PREMISE

Animal pollinators play an important role in pollen dispersal. Here, we assessed differences in pollen and seed dispersal and the role of pollinator functional groups with different foraging behaviors in generating patterns of genetic diversity over similar geographic ranges for two closely related taxa. We focused on two members of Oenothera section Calylophus (Onagraceae) that co-occur on gypsum outcrops throughout the northern Chihuahuan Desert but differ in floral phenotype and primary pollinator: Oenothera gayleana (bee) and O. hartwegii subsp. filifolia (hawkmoth).

METHODS

We measured breeding system and floral traits and studied gene flow and population differentiation at the local (<13 km; four populations) and landscape (60-440 km; five populations) scales using 10-11 nuclear (pollen dispersal) and three plastid (seed dispersal) microsatellite markers.

RESULTS

Both taxa were self-incompatible and floral traits were consistent with expectations for different pollinators. Seed and pollen dispersal patterns were distinctly different for both species. We found no evidence of genetic structure at the local scale but did at the landscape scale; O. gayleana showed greater differentiation and significant isolation by distance than in O. hartwegii subsp. filifolia. The plastid data were consistent with gravity dispersal of seeds and suggest that pollen dispersal is the principal driver of genetic structure in both species.

CONCLUSIONS

We demonstrated that pollinator functional groups can impact genetic differentiation in different and predictable ways. Hawkmoths, with larger foraging distances, can maintain gene flow across greater spatial scales than bees.

Phylogeography of the Plateau Pika (Ochotona curzoniae) in Response to the Uplift of the Qinghai-Tibet Plateau

The evolution and current distribution of species on the Qinghai-Tibet Plateau have been significantly impacted by historical occurrences, including the uplift of the plateau and the Quaternary climate upheaval. As a remnant species, the plateau pika (Ochotona curzoniae) is a great model for researching historical events. In this study, 302 samples from 42 sample sites were utilized to analyze the impact of historical events on the evolution and distribution pattern of plateau pikas. The genetic diversity, patterns of differentiation, and historical dynamics of the plateau pika were investigated using molecular markers that included four mitochondrial genes (COI, D-loop, Cytb, and 12S rRNA) and three nuclear genes (GHR, IRBP, and RAG1). The results showed that: (1) The genetic diversity of the plateau pika was high in the Tibetan Plateau (Hd = 0.9997, π = 0.01205), and the plateau pika evolved into five lineages that occupied different geographical areas, with lineage 1 (Group 1) in the south of the Yarlung Zangbo River, lineage 2 (Group 2) in the hinterland of the plateau, lineage 3 (Group 3) in the northeastern part of the plateau, lineage 4 (Group 4) in the Hengduan Mountains, and lineage 5 (Group 5) in the eastern part of the plateau. (2) The gene flow among the five lineages was low, and the differentiation level was high (Nm < 0.25; Fst > 0.25), indicating that the geographical barriers between the five lineages, such as the Yarlung Zangbo River, the Qaidam-Ghuong-Guide Basin, and the Lancang River, effectively promoted the population differentiation of the plateau pika. (3) The plateau pika first spread from the Hengduan Mountains to the entire Qinghai-Tibet Plateau and then conducted small-scale migration and dispersal in several refuges across the plateau in response to climate changes during the glacial and interglacial periods. (4) Except for Group 1 and Group 4, all the other populations exhibited a rapid expansion between 0.06 and 0.01 Mya, but the expansion was considerably delayed or halted by the effects of climate change during the last glacial maximum (0.02 Mya). Overall, the plateau pika on the Qinghai-Tibet Plateau exhibits high genetic diversity, and topographic obstacles, including mountains, valleys, and basins, created by the uplift of the plateau and climatic changes since the Quaternary period have played an important role in the differentiation and historical dynamics of the plateau pika population.

Genetic Diversity and Population Structure of Cowpea (Vigna unguiculata (L.) Walp.) Landraces from Portugal and Mozambique

Cowpea (Vigna unguiculata (L.) Walp.) is currently a legume crop of minor importance in Europe but a highly relevant staple crop in Africa and the second most cultivated legume in Mozambique. In Portugal and Mozambique, cowpea's phenotypic and genetic variation has been maintained locally by farmers in some areas. We used the molecular markers SSR, SilicoDArT and SNP to analyze the genetic diversity and population structure of 97 cowpea accessions, mainly from Portugal (Southern Europe) and Mozambique (Southern Africa). As far as we know, this is the first time that the genetic variation and the relationship between cowpea landraces collected in Portugal with those originated in Mozambique is reported. Despite the shared historical past, the Portuguese landraces did not share a common genetic background with those from Mozambique, and two different gene pools were revealed. Knowledge of the genetic structure of cowpea landraces offers an opportunity for individual selection within landraces adapted to particular eco-physiological conditions and suggests the existence of a valuable gene pool for exploitation in future Portugal-PALOP (Portuguese-speaking African countries) cowpea breeding programs.

Association of BMAL1 clock gene polymorphisms with fasting glucose in children

BACKGROUND

The brain and muscle Arnt-like protein-1 (BMAL1) gene is an important circadian clock gene and previous studies have found that certain polymorphisms are associated with type 2 diabetes in adults. However, it remains unknown if such polymorphisms can affect fasting glucose in children and if other factors modify the associations.

METHODS

A school-based cross-sectional study with 947 Chinese children was conducted. A multivariable linear regression model was used to analyze the association between BMAL1 gene polymorphisms and fasting glucose level.

RESULTS

After adjusting for age, sex, body mass index (BMI), physical activity, and unhealthy diet, GG genotype carriers of BMAL1 rs3789327 had higher fasting glucose than AA/GA genotype carriers (b = 0.101, SE = 0.050, P = 0.045). Adjusting for the same confounders, rs3816358 was shown to be significantly associated with fasting glucose (b = 0.060, SE = 0.028, P = 0.032). Furthermore, a significant interaction between rs3789327 and nutritional status on fasting glucose was identified (P_interaction = 0.009); rs3789327 was associated with fasting glucose in the overweight/obese subgroup (b = 0.353, SE = 0.126, P = 0.006), but not in non-overweight/non-obese children.

CONCLUSIONS

BMAL1 polymorphisms were significantly associated with the fasting glucose level in children. Additionally, the observed interaction between nutritional status and BMAL1 supports promoting an optimal BMI in children genetically predisposed to higher glucose level.

IMPACT

Polymorphisms in the essential circadian clock gene BMAL1 were associated with fasting blood glucose levels in children. Additionally, there was a significant interaction between nutritional status and BMAL1 affecting fasting glucose levels. BMAL1 rs3789327 was associated with fasting glucose only in overweight/obese children. This finding could bring novel insights into mechanisms by which nutritional status influences fasting glucose in children.

Genome-wide re-sequencing data reveals the genetic diversity and population structure of Wenchang chicken in China

Wenchang (WC) chicken, the only indigenous chicken breed listed in Chinese genetic resources in Hainan province, is well known for its excellent meat quality and is sold all over southeast Asia. In recent years, the number of WC has decreased sharply with considerable variability in the quality at market. To explore the genetic diversity and population structure of WC chickens, the whole-genome data of 235 WC individuals from three conservation farms were obtained using the Illumina 150 bp paired-end platform and used in conjunction with the sequencing data from 123 individuals from other chicken breeds (including eight Chinese indigenous chicken breeds and three foreign or commercial breeds) downloaded from a public database. A total of 12 111 532 SNPs were identified, of which 11 541 878 SNPs were identified in WC. The results of gene enrichment analyses revealed that the SNPs harbored in WC genomes are mainly related to environmental adaptation, disease resistance and meat quality traits. Genetic diversity statistics, quantified by expected heterozygosity, observed heterozygosity, linkage disequilibrium, nucleotide diversity and fixation statistics, indicated that WC displays high genetic diversity compared with other Chinese indigenous chicken breeds. Genetic structure analyses showed that each population displayed great differentiation between WC and the other breeds, indicating the uniqueness of WC. In conclusion, the results of our study provide the first genomic overview of genetic variants, genetic diversity and population structure of WC from three conservation farms. This information will be valuable for the future breeding and conservation of WC and other surveyed populations.

Outlier analyses and genome-wide association study identify glgC and ERD6-like 4 as candidate genes for foliar water-soluble carbohydrate accumulation in Trifolium repens

Increasing water-soluble carbohydrate (WSC) content in white clover is important for improving nutritional quality and reducing environmental impacts from pastoral agriculture. Elucidation of genes responsible for foliar WSC variation would enhance genetic improvement by enabling molecular breeding approaches. The aim of the present study was to identify single nucleotide polymorphisms (SNPs) associated with variation in foliar WSC in white clover. A set of 935 white clover individuals, randomly sampled from five breeding pools selectively bred for divergent (low or high) WSC content, were assessed with 14,743 genotyping-by-sequencing SNPs, using three outlier detection methods: PCAdapt, BayeScan and KGD-F_ST. These analyses identified 33 SNPs as discriminating between high and low WSC populations and putatively under selection. One SNP was located in the intron of ERD6-like 4, a gene coding for a sugar transporter located on the vacuole membrane. A genome-wide association study using a subset of 605 white clover individuals and 5,757 SNPs, identified a further 12 SNPs, one of which was associated with a starch biosynthesis gene, glucose-1-phosphate adenylyltransferase, glgC. Our results provide insight into genomic regions underlying WSC accumulation in white clover, identify candidate genomic regions for further functional validation studies, and reveal valuable information for marker-assisted or genomic selection in white clover.

Genetic structure of Spirometra mansoni (Cestoda: Diphyllobothriidae) populations in China revealed by a Target SSR-seq method

BACKGROUND

In China, the plerocercoid of the cestode Spirometra mansoni is the main causative agent of human and animal sparganosis. However, the population genetic structure of this parasite remains unclear. In this study, we genotyped S. mansoni isolates with the aim to improve current knowledge on the evolution and population diversity of this cestode.

METHODS

We first screened 34 perfect simple sequence repeats (SSRs) using all available omic data and then constructed target sequencing technology (Target SSR-seq) based on the Illumina NovaSeq platform. Next, a series of STRUCTURE. clustering, principal component, analysis of molecular variance and TreeMix analyses were performed on 362 worm samples isolated from 12 different hosts in 16 geographical populations of China to identify the genetic structure.

RESULTS

A total of 170 alleles were detected. The whole population could be organized and was found to be derived from the admixture of two ancestral clusters. TreeMix analysis hinted that possible gene flow occurred from Guizhou (GZ) to Sichuan (SC), SC to Jaingxi (JX), SC to Hubei (HB), GZ to Yunnan (YN) and GZ to Jiangsu (JS). Both neighbor-joining clustering and principal coordinate analysis showed that isolates from intermediate hosts tend to cluster together, while parasites from definitive hosts revealed greater genetic differences. Generally, a S. mansoni population was observed to harbor high genetic diversity, moderate genetic differentiation and a little genetic exchange among geographical populations.

CONCLUSIONS

A Target SSR-seq genotyping method was successfully developed, and an in-depth view of genetic diversity and genetic relationship will have important implications for the prevention and control of sparganosis.

Current Situation and Utilization of Velvet Deer Germplasm Resources in China

Velvet deer are not only a representative special economic animal but also an important part of livestock. With the increasing awareness of international competition for germplasm resources in China, more and more attention has been paid to the protection and utilization of germplasm resources. However, there is poor understanding about velvet deer resources. Therefore, we are providing a comprehensive introduction of Chinese velvet deer germplasm resources from the aspects of ecological distribution, domestication and breeding.

Spatial and temporal variation of genetic diversity and genetic differentiation in Daphnia galeata populations in four large reservoirs in southern China

Daphnia galeata is a common and dominant species in warmer waters, and has a strong top-down effect on both phytoplankton and bacteria. The knowledge of its temporal and spatial patterns of genetic diversity is fundamental in understanding its population dynamics and potential ecological function in ecosystems. Its population genetics have been investigated at regional scales but few within regions or at smaller spatial scales. Here, we examined the fine-scale spatial genetic variation of D. galeata within four large, deep reservoirs in wet and dry seasons and the six-year variation of genetic diversity in one of the reservoirs by using cytochrome c oxidase subunit I and microsatellites (simple sequence repeat). Our study shows that fine-scale spatial genetic variation commonly occurred within the reservoirs, indicating strong environmental selection at least in the two of reservoirs with strong longitudinal gradients. Since the environmental gradients established in the dry season was largely reduced in the wet season, the fine-scale spatial genetic variation was much higher in the dry season. The dynamics of local genetic diversity did not follow the theoretical pattern of rapid erosion but peaked in mid or mid-late growth season. The local genetic diversity of D. galeata appears to be shaped and maintained not only by recruitment from resting egg banks but also by gene flow within reservoirs. The temporal and fine-scale genetic variation within a water body suggests that it is necessary to pay attention to sampling periods and locations of a given water body in regional studies.

Microsatellite analysis and polymorphic marker development based on the full-length transcriptome of Camellia chekiangoleosa

Camellia chekiangoleosa is a popular variety of Oil-camellia that has high oil production and ornamental value. Microsatellite (SSR) markers are the preferred tool for the molecular marker-assisted breeding of C. chekiangoleosa. By focusing on the problems of the low development efficiency of polymorphic SSR markers and the lack of available functional markers in Oil-camellia, we identified 97,510 SSR loci based on the full-length transcriptome sequence of C. chekiangoleosa. An analysis of SSR characteristics showed that mononucleotide (51.29%) and dinucleotide (34.36%) SSRs were the main repeat types. The main SSR distribution areas based on proportion covered were ordered as follows: 5'UTR > 3'UTR > CDS. By comparing our data with those in databases such as GO and KEGG, we obtained functional annotations of unigene sequences containing SSR sites. The data showed that the amplification efficiency of the SSR primers was 51.72%, and the development efficiency of polymorphic SSR primers was 26.72%. Experiments verified that dinucleotide and pentanucleotide SSRs located in UTR regions could produce more polymorphic markers. An investigation into the genetic diversity of several C. chekiangoleosa populations also suggested that the developed SSR markers had higher levels of polymorphism. This study will provide a reference and high-quality markers for the large-scale development of functional SSR markers and genetic research in Oil-camellia.

Genetic diversity in the transmission-blocking vaccine candidate Plasmodium vivax gametocyte protein Pvs230 from the China-Myanmar border area and central Myanmar

BACKGROUND

Sexual stage surface antigens are potential targets of transmission-blocking vaccines (TBVs). The gametocyte and gamete surface antigen P230, a leading TBV candidate, is critical for red blood cell binding during exflagellation and subsequent oocyst development. Here, the genetic diversity of Pvs230 was studied in Plasmodium vivax parasite isolates from the China-Myanmar border (CMB) and central Myanmar.

METHODS

Plasmodium vivax isolates were collected in clinics from malaria-endemic areas of the CMB (143 samples) and Myanmar (23 samples). The interspecies variable part (IVP, nucleotides 1-807) and interspecies conserved part (ICP, 808-2862) of Pvs230 were amplified by PCR and sequenced. Molecular evolution studies were conducted to evaluate the genetic diversity, signature of selection, population differentiation, haplotype network, and population structure of the study parasite populations and publicly available Pvs230 sequences from six global P. vivax populations.

RESULTS

Limited genetic diversity was observed for the CMB (π = 0.002) and Myanmar (π = 0.001) isolates. Most amino acid substitutions were located in the IVP and cysteine-rich domain of Pvs230. Evidence of positive selection was observed for IVP and purifying selection for ICP. Codon-based tests identified specific codons under natural selection in both IVP and ICP. The fixation index (F_ST) showed low genetic differentiation between East and Southeast Asian populations, with F_ST ranging from 0.018 to 0.119. The highest F_ST value (F_ST = 0.503) was detected between the Turkey and Papua New Guinea populations. A total of 92 haplotypes were identified in global isolates, with the major haplotypes 2 and 9 being the most abundant and circulating in East and Southeast Asia populations. Several detected non-synonymous substitutions were mapped in the predicted structure and B-cell epitopes of Pvs230.

CONCLUSIONS

We detected low levels of genetic diversity of Pvs230 in global P. vivax populations. Geographically specific haplotypes were identified for Pvs230. Some mutations are located within a potential B-cell epitope region and need to be considered in future TBV designs.

Population genetic structure of Wikstroemia monnula highlights the necessity and feasibility of hierarchical analysis for a highly differentiated species

Population genetic structure can provide valuable insights for conserving genetic resources and understanding population evolution, but it is often underestimated when using the most popular method and software, STRUCTURE and delta K, to assess. Although the hierarchical STRUCTURE analysis has been proposed early to overcome the above potential problems, this method was just utilized in a few studies and its reliability needs to be further tested. In this study, the genetic structure of populations of Wikstroemia monnula was evaluated by sequencing 12 nuclear microsatellite loci of 905 individuals from 38 populations. The STRUCTURE analysis suggested the most likely number of clusters was two, but using multi-hierarchical structure analysis, almost every population was determined with an endemic genetic component. The latter result is consistent with the extremely low gene flow among populations and a large number of unique cpDNA haplotypes in this species, indicating one level of structure analysis would extremely underestimate its genetic component. The simulation analysis shows the number of populations and the genetic dispersion among populations are two key factors to affect the estimation of K value using the above tools. When the number of populations is more than a certain amount, K always is equal to 2, and when a simulation only includes few populations, the underestimation of K value also may occur only if these populations consist of two main types of significantly differentiated genetic components. Our results strongly support that the hierarchical STRUCTURE analysis is necessary and practicable for the species with lots of subdivisions.

The genetic diversity and structure in the European polecat were not affected by the introduction of the American mink in Poland

The introduction and expansion of an invasive non-native species could have important consequences for the genetic patterns and processes of native species, moreover if the new arrival competes strongly for resources and space. This may result in the demographic decline of the native species. Knowing the effects on the levels of genetic diversity and structure in native species is key in terms of their conservation. We analysed temporal (over 50 years) genetic variation of the population of the European polecat (Mustela putorius), a species under threat in several European countries, in the Białowieża Primeval Forest (BPF), Poland, before and after the invasion of the American mink (Neovison vison). Using 11 microsatellite loci and a fragment of the mitochondrial control region we show that levels of diversity changed in the polecat population over 53 generations (over the period 1959–2012) and after the invasion of mink. When compared with other threatened European polecat populations, high levels of diversity are observed in the population in BPF in both periods, as well as in other areas in Poland. Our data shows that genetic structure was not present either before or after the mink invasion in BPF. This would suggest that the polecat population in Poland was not affected by invasive species and other negative factors and would be a potential good source of individuals for captive breeding or genetic rescue conservation management actions in areas where such actions are needed, for example the UK.

Fine-scale genetic structure and wolbachia infection of aedes albopictus (Diptera: Culicidae) in Nanjing city, China

Background:Aedes albopictus is an indigenous primary vector of dengue and Zika viruses in China. Wolbachia is a gram-negative and common intracellular bacteria, which is maternally inherited endosymbionts and could expand their propagation in host populations by means of various manipulations. Compared with research on the dispersion of Ae. albopictus at the macrospatial level (mainly at the country or continent level), little is known about its variation and Wolbachia infection at the microspatial level, which is essential for its management. Meanwhile, no local cases of dengue fever have been recorded in the history of Nanjing, which implies that few adulticides have been applied in the city. Thus, the present study examines how the Ae. albopictus population varies and the Wolbachia infection status of each population among microspatial regions of Nanjing City.

Methods: The genetic structure of 17 Aedes albopictus populations collected from urban, urban fringe, and rural regions of Nanjing City was investigated based on 9 microsatellite loci and the mitochondrial coxI gene. The Wolbachia infection status of each population was also assessed with Wolbachia A- and Wolbachia B-specific primers.

Results: Nine out of 58 tested pairs of microsatellite markers were highly polymorphic, with a mean PIC value of 0.560, and these markers were therefore chosen for microsatellite genotyping analysis. The Na value of each Ae. albopictus population was very high, and the urban area populations (7.353 ± 4.975) showed a lower mean value than the urban fringe region populations (7.866 ± 5.010). A total of 19 coxI haplotypes were observed among 329 Ae. albopictus individuals via haplotype genotyping, with the highest diversity observed among the urban fringe Ae. albopictus populations (Hd = 0.456) and the lowest among the urban populations (Hd = 0.277). Each Ae. albopictus population showed significant departure from HWE, and significant population expansion was observed in only three populations from the urban (ZSL), urban fringe (HAJY), and rural areas (HSZY) (p < 0.05). Combined with DAPC analysis, all the Ae. albopictus populations were adequately allocated to two clades with significant genetic differences according to population structure analysis, and the best K value was equal to two. AMOVA results showed that most (96.18%) of the genetic variation detected in Ae. albopictus occurred within individuals (F_IT = 0.22238, p < 0.0001), while no significant positive correlation was observed via isolation by distance (IBD) analysis (R² = 0.03262, p = 0.584). The TCS network of all haplotypes showed that haplotype 1 (H1) and haplotype 4 (H4) were the most frequent haplotypes among all populations, and the haplotype frequency significantly increased from urban regions (36.84%) to rural regions (68.42%). Frequent migration was observed among Ae. albopictus populations from rural to urban regions via the urban fringe region, with four direct migration routes between rural and urban regions. Furthermore, Wolbachia genotyping results showed that most of the individuals of each population were coinfected with Wolbachia A and Wolbachia B. The independent infection rate of Wolbachia A was slightly higher than that of Wolbachia B, and no significant differences were observed among different regions.

Conclusion: In the microspatial environment of Nanjing City, the urban fringe region is an important region for the dispersion of Ae. albopictus populations between rural and urban areas, and Wolbachia A and Wolbachia B coinfection is the most common Wolbachia infection status in all Ae. albopictus populations among different regions.

Tunisian camel casein gene characterization reveals similarities and differences with Sudanese and Nigerian populations

Milk is a primary protein source that has always played a role in mammalian health. Despite the intensification of research projects on dromedary and the knowledge of the genetic diversity at the casein loci, the genetic structure of the Tunisian camel population still needs exploration. This study sought to determine the genetic diversity of 3 casein gene variants in 5 Tunisian camel ecotypes: c.150G>T at CSN1S1 (α_S1-casein), g.2126A>G at CSN2 (β-casein), and g.1029T>C at CSN3 (κ-casein). The obtained results were compared with data published on Sudanese and Nigerian camels to establish the level of differentiation within and between populations. A total of 159 blood samples were collected from 5 Tunisian camel ecotypes and the extracted DNA was genotyped by PCR-RFLP. A streamlined genotyping protocol was also developed for CSN3. Results indicated that allele T was quite rare (0.06) at CSN1S1 for all ecotypes. Minor allele frequency was found for G (0.462) in CSN2 except for Ardhaoui Medenine ecotype who deviated from the average CSN2 allele frequency of the total population. Allele C showed minor allele frequency of 0.384 in CSN3. Among the Tunisian population, GAT (0.343) was the most represented haplotype in all ecotypes except for Ardhaoui Medenine, where GGC (0.322) was the most frequent one. Significant differences in heterozygosity and local inbreeding were observed across the Tunisian, Sudanese, and Nigerian populations, although the global fixation index indicated that only 2.2% of the genetic variance is related to ecotype differences. Instead, phylogenetic analysis revealed a closer link between the Tunisian and Sudanese populations through a clade subdivision with 3 main branches among the ecotypes. This study represents the first attempt to understand casein gene variability in Tunisian camels; with further study, milk traits and genetic differentiation among populations can be associated with the history of camel domestication.

Microsatellite Characteristics of Silver Carp (Hypophthalmichthysmolitrix) Genome and Genetic Diversity Analysis in Four Cultured Populations

Hypophthalmichthys molitrix is one of the four most important fish in China and has high breeding potential. However, simple sequence repeat (SSR) markers developed on H. molitrix genome level for genetic diversity analysis are limited. In this study, the distribution characteristics of SSRs in the assembled H. molitrix genome were analyzed, and new markers were developed to preliminarily evaluate the genetic diversity of the four breeding populations. A total of 368,572 SSRs were identified from the H. molitrix genome. The total length of SSRs was 6,492,076 bp, accounting for 0.77% of the total length of the genome sequence. The total frequency and total density were 437.73 loci/Mb and 7713.16 bp/Mb, respectively. Among the 2–6 different nucleotide repeat types, SSRs were dominated by di-nucleotide repeats (204,873, 55.59%), and AC/GT was the most abundant motif. The number of SSRs on each chromosome was positively correlated with the length. The 13 pairs of markers developed were used to analyze the genetic diversity of four cultivated populations in Hubei Province. The results showed that the genetic diversity of the four populations was low, and the ranges of alleles (Na), effective alleles (Ne), observed heterozygosity (Ho), and Shannon’s index information (I) were 3.538–4.462, 2.045–2.461, 0.392–0.450, and 0.879–0.954, respectively. Genetic variation occurs mainly among individuals within populations (95.35%). UPGMA tree and Bayesian analysis showed that four populations could be divided into two different branches. Therefore, the genome-wide SSRs were effectively in genetic diversity analysis on H. molitrix.

Genetic diversity, gene flow, and differentiation among wild, semiwild, and landrace chile pepper (Capsicum annuum) populations in Oaxaca, Mexico

PREMISE

Capsicum annuum (Solanaceae) was originally domesticated in Mexico, where wild (C. annuum var. glabriusculum) and cultivated (C. annuum var. annuum) chile pepper populations (>60 landraces) are common, and wild-resembling individuals (hereafter semiwild) grow spontaneously in anthropogenic environments. Here we analyze the role of elevation and domestication gradients in shaping the genetic diversity in C. annuum from the state of Oaxaca, Mexico.

METHODS

We collected samples of 341 individuals from 28 populations, corresponding to wild, semiwild (C. annuum var. glabriusculum) and cultivated C. annuum, and closely related species Capsicum frutescens and C. chinense. From the genetic variation of 10 simple sequence repeat (SSR) loci, we assessed the population genetic structure, inbreeding, and gene flow through variance distribution analyses, genetic clustering, and connectivity estimations.

RESULTS

Genetic diversity (H_E ) did not differ across domestication levels. However, inbreeding coefficients were higher in semiwild and cultivated chiles than in wild populations. We found evidence for gene flow between wild populations and cultivated landraces along the coast. Genetic structure analysis revealed strong differentiation between most highland and lowland landraces.

CONCLUSIONS

Gene flow between wild and domesticated populations may be mediated by backyards and smallholder farms, while mating systems may facilitate gene flow between landraces and semiwild populations. Domestication and elevation may overlap in their influence on genetic differentiation. Lowland Gui'ña dani clustered with highland landraces perhaps due to the social history of the Zapotec peoples. In situ conservation may play an important role in preserving semiwild populations and private alleles found in landraces.

Genetic Variations of the DPYD Gene and Its Relationship with Ancestry Proportions in Different Ecuadorian Trihybrid Populations

Dihydropyrimidine dehydrogenase is one of the main pharmacological metabolizers of fluoropyrimidines, a group of drugs widely used in clinical oncology. Around 20 to 30% of patients treated with fluoropyrimidines experience severe toxicity caused by a partial or total decrease in enzymatic activity. This decrease is due to molecular variants in the DPYD gene. Their prevalence and allelic frequencies vary considerably worldwide, so their description in heterogeneous groups such as the Ecuadorian population will allow for the description of pharmacogenetic variants and proper characterization of this population. Thus, we genotyped all the molecular variants with a predictive value for DPYD in a total of 410 Ecuadorian individuals belonging to Mestizo, Afro-Ecuadorian, and Indigenous ethnic groups. Moreover, we developed a genetic ancestry analysis using 46 autosomal ancestry informative markers. We determined 20 genetic variations in 5 amplified regions, including 3 novel single nucleotide variants. The allele frequencies for DPYD variants c.1627G>A (*5, rs1801159), c.1129-15T>C (rs56293913), c.1218G>A (rs61622928), rs1337752, rs141050810, rs2786783, rs2811178, and g.97450142G>A (chr1, GRCh38.p13) are significantly related to Native American and African ancestry proportions. In addition, the FST calculated from these variants demonstrates the closeness between Indigenous and Mestizo populations, and evidences genetic divergence between Afro-Ecuadorian groups when compared with Mestizo and Indigenous ethnic groups. In conclusion, the genetic variability in the DPYD gene is related to the genetic component of ancestral populations in different Ecuadorian ethnic groups. The absence and low frequency of variants with predictive value for fluoropyrimidine toxicity such as DPYD *2A, HapB3, and c.2846A>T (prevalent in populations with European ancestry) is consistent with the genetic background found.

Genetic diversity of Prunus armeniaca L. var. ansu Maxim. germplasm revealed by simple sequence repeat (SSR) markers

The genetic diversity and genetic structure of P. armeniaca var. ansu were analyzed based on SSR markers. The aim was to provide scientific basis for conservation, efficient utilization, molecular marker assisted breeding and improved variety selection of P. armeniaca var. ansu germplasm resources. The results showed that the level of genetic diversity within the population was high. Among the 30 SSR markers, the mean number of observed alleles was 11.433, the mean number of effective alleles was 4.433, the mean of Shannon information index was 1.670, and the mean of polymorphic information content was 0.670. Among the eight provenances, Tuanjie Township, Xinyuan County, Xinjiang had the highest genetic diversity. The observed alleles, effective alleles, Shannon information index and Nei’s gene diversity index among provenances were higher than those within provenances. Based on Bayesian mathematical modeling and UPGMA cluster analysis, 86 P. armeniaca var. ansu accessions were divided into three subpopulations and four groups, which reflected individual differences in provenances. Subpopulations classified by Bayesian mathematical modeling and groups classified by UPGMA cluster analysis were significantly correlated with geographical provenance (Sig<0.01) and the provenances significantly impacted classification of groups. The provenances played an important role in classification of groups. The genetic distance between Tuanjie Township of Xinyuan County and Alemale Township of Xinyuan County was the smallest, while the genetic relationship between them was the closest and the degree of genetic differentiation was small.