Studies on genetic diversity, gene flow and landscape genetic in Avicennia marina: Spatial PCA, Random Forest, and phylogeography approaches

Mangrove forests grow in coastal areas, lagoons, estuaries, and deltas and form the main vegetation in tidal and saline wetlands. Due to the mankind activities and also changes in climate, these forests face degradations and probably extinction in some areas. Avicennia marina is one of the most distributed mangrove species throughout the world. The populations of A. marina occur in a limited region in southern parts of Iran. Very few genetic and spatial analyses are available on these plants from our country. Therefore, the present study was planned to provide detailed information on Avicennia marina populations with regard to genetic diversity, gene flow versus genetic isolation, effects of spatial variables on connectivity and structuring the genetic content of trees populations and also identifying adaptive genetic regions in respond too spatial variables. We used SCoT molecular markers for genetic analyses and utilized different computational approaches for population genetics and landscapes analyses. The results of present study showed a low to moderate genetic diversity in the studied populations and presence of significant Fst values among them. Genetic fragmentation was also observed within each province studied. A limited gene flow was noticed among neighboring populations within a particular province. One population was almost completely isolated from the gene flow with other populations and had peculiar genetic content.Spatial PCA analysis revealed both significant global and local genetic structuring in the studied populations. Spatial variables like humidity, longitude and altitude were the most important spatial features affecting genetic structure in these populations.

Mitochondrial sequence data reveal population structure within Pustulosa pustulosa

Unionid mussels are among the most imperiled group of organisms in North America, and Pustulosa pustulosa is a freshwater species with a relatively wide latitudinal distribution that extends from southern Ontario, Canada, to Texas, USA. Considerable morphological and geographic variation in the genus Pustulosa (formerly Cyclonaias) has led to uncertainty over species boundaries, and recent studies have suggested revisions to species-level classifications by synonymizing C. aurea, C. houstonensis, C. mortoni, and C. refulgens with C. pustulosa (currently P. pustulosa). Owing to its wide range and shallow phylogenetic differentiation, we analyzed individuals of P. pustulosa using mitochondrial DNA sequence data under a population genetics framework. We included 496 individuals, which were comprised of 166 samples collected during this study and 330 additional sequences retrieved from GenBank. Pairwise ΦST measures based on ND1 data suggested there may be up to five major geographic groups present within P. pustulosa. Genetic differentiation between regions within Texas was higher compared to populations from the Mississippi and Great Lakes populations, which may reflect differences in historical connectivity. Mitochondrial sequence data also revealed varying demographic histories for each major group suggesting each geographic region has also experienced differential population dynamics in the past. Future surveys should consider exploring variation within species after phylogeographic delimitation has been performed. In this study, we begin to address this need for freshwater mussels via the P. pustulosa system.

Genomic insights into selection for heterozygous alleles and woody traits in Populus tomentosa

Heterozygous alleles are widespread in outcrossing and clonally propagated woody plants. The variation in heterozygosity that underlies population adaptive evolution and phenotypic variation, however, remains largely unknown. Here, we describe a de novo chromosome-level genome assembly of Populus tomentosa, an economic and ecologically important native tree in northern China. By resequencing 302 natural accessions, we determined that the South subpopulation (Pop_S) encompasses the ancestral strains of P. tomentosa, while the Northwest subpopulation (Pop_NW) and Northeast subpopulation (Pop_NE) experienced different selection pressures during population evolution, resulting in significant population differentiation and a decrease in the extent of heterozygosity. Analysis of heterozygous selective sweep regions (HSSR) suggested that selection for lower heterozygosity contributed to the local adaptation of P. tomentosa by dwindling gene expression and genetic load in the Pop_NW and Pop_NE subpopulations. Genome-wide association studies (GWAS) revealed that 88 single nucleotide polymorphisms (SNPs) within 63 genes are associated with nine wood composition traits. Among them, the selection for the homozygous AA allele in PtoARF8 is associated with reductions in cellulose and hemicellulose contents by attenuating PtoARF8 expression, and the increase in lignin content is attributable to the selection for decreases in exon heterozygosity in PtoLOX3 during adaptive evolution of natural populations. This study provides novel insights into allelic variations in heterozygosity associated with adaptive evolution of P. tomentosa in response to the local environment and identifies a series of key genes for wood component traits, thereby facilitating genomic-based breeding of important traits in perennial woody plants.

Molecular variability of the Ancylostoma secreted Protein-2 (Aca-asp-2) gene from Ancylostoma caninum contributes to expand information on population genetic studies of hookworms

Hookworm infection is a major public health problem in many regions of the world. Given the high levels of host morbidity and even mortality of the host caused by these infections, it is crucial to understand the genetic structure of hookworm populations. This understanding can provide insights into the ecology, transmission patterns, mechanisms of drug resistance, and the development of vaccines and immunotherapeutic strategies. Previously, we examined presumably neutral molecular markers, such as microsatellites and COI (Cytochrome C oxidase subunit 1) in Brazilian populations of Ancylostoma caninum. Here we analyze the molecular variability of a genomic fragment of the Aca-asp-2 (Ancylostoma secreted protein-2) gene from Ancylostoma caninum. This gene is a highly expressed and activated following the infection of the L3 larvae in the host. We obtained individuals of A. caninum from five different geographic locations in Brazil, sequenced and analyzed parts of the gene. The results revealed extensive polymorphism at this fragment, especially in the intronic region, indicating low selective pressure acting on these sequences. However, we also observed irregular distributions of nucleotides and polymorphisms in the coding region of this gene, resulting in the identification of 27 alleles. The data presented here contribute to expanding the understanding of population genetic studies of hookworms.

Comprehensive analyses of genetic diversities and population structure of the Guizhou Dong group based on 44 Y-markers

Background

The non-recombining region of the human Y chromosome (NRY) is a strictly paternally inherited genetic marker and the best material to trace the paternal lineages of populations. Y chromosomal short tandem repeat (Y-STR) is characterized by high polymorphism and paternal inheritance pattern, so it has been widely used in forensic medicine and population genetic research. This study aims to understand the genetic distribution of Y-STRs in the Guizhou Dong population, provide reference data for forensic application, and explore the phylogenetic relationships between the Guizhou Dong population and other comparison populations.

Methods

Based on the allele profile of 44 Y-markers in the Guizhou Dong group, we estimate their allele frequencies and haplotype frequencies. In addition, we also compare the forensic application efficiency of different Y-STR sets in the Guizhou Dong group. Finally, genetic relationships among Guizhou Dong and other reference populations are dissected by the multi-dimensional scaling and the phylogenetic tree.

Results

A total of 393 alleles are observed in 312 Guizhou Dong individuals for these Y-markers, with allele frequencies ranging from 0.0032 to 0.9679. The haplotype diversity and discriminatory capacity for these Y-markers in the Guizhou Dong population are 0.99984 and 0.97440, respectively. The population genetic analyses of the Guizhou Dong group and other reference populations show that the Guizhou Dong group has the closest genetic relationship with the Hunan Dong population, and followed by the Guizhou Tujia population.

Conclusions

In conclusion, these 44 Y-markers can be used as an effective tool for male differentiation in the Guizhou Dong group. The haplotype data in this study not only enrich the Y-STR data of different ethnic groups in China, but also have important significance for population genetics and forensic research.

Population genetics of marmosets in Asian primate research centers and loci associated with epileptic risk revealed by whole-genome sequencing

The common marmoset ( Callithrix jacchus) has emerged as a valuable nonhuman primate model in biomedical research with the recent release of high-quality reference genome assemblies. Epileptic marmosets have been independently reported in two Asian primate research centers. Nevertheless, the population genetics within these primate centers and the specific genetic variants associated with epilepsy in marmosets have not yet been elucidated. Here, we characterized the genetic relationships and risk variants for epilepsy in 41 samples from two epileptic marmoset pedigrees using whole-genome sequencing. We identified 14 558 184 single nucleotide polymorphisms (SNPs) from the 41 samples and found higher chimerism levels in blood samples than in fingernail samples. Genetic analysis showed fourth-degree of relatedness among marmosets at the primate centers. In addition, SNP and copy number variation (CNV) analyses suggested that the WW domain-containing oxidoreductase ( WWOX) and Tyrosine-protein phosphatase nonreceptor type 21 ( PTPN21) genes may be associated with epilepsy in marmosets. Notably, KCTD18-like gene deletion was more common in epileptic marmosets than control marmosets. This study provides valuable population genomic resources for marmosets in two Asian primate centers. Genetic analyses identified a reasonable breeding strategy for genetic diversity maintenance in the two centers, while the case-control study revealed potential risk genes/variants associated with epilepsy in marmosets.

Genetic modification of inflammation- and clonal hematopoiesis-associated cardiovascular risk

Clonal hematopoiesis of indeterminate potential (CHIP) is associated with an increased risk of cardiovascular diseases (CVDs), putatively via inflammasome activation. We pursued an inflammatory gene modifier scan for CHIP-associated CVD risk among 424,651 UK Biobank participants. We identified CHIP using whole-exome sequencing data of blood DNA and modeled as a composite, considering all driver genes together, as well as separately for common drivers (DNMT3A, TET2, ASXL1, and JAK2). We developed predicted gene expression scores for 26 inflammasome-related genes and assessed how they modify CHIP-associated CVD risk. We identified IL1RAP as a potential key molecule for CHIP-associated CVD risk across genes and increased AIM2 gene expression leading to heightened JAK2- and ASXL1-associated CVD risk. We show that CRISPR-induced Asxl1-mutated murine macrophages had a particularly heightened inflammatory response to AIM2 agonism, associated with an increased DNA damage response, as well as increased IL-10 secretion, mirroring a CVD-protective effect of IL10 expression in ASXL1 CHIP. Our study supports the role of inflammasomes in CHIP-associated CVD and provides evidence to support gene-specific strategies to address CHIP-associated CVD risk.

Development and characterization of microsatellite markers for the tsetse species Glossina brevipalpis and preliminary population genetics analyses

Tsetse flies, the vectors of African trypanosomes are of key medical and economic importance and one of the constraints for the development of Africa. Tsetse fly control is one of the most effective and sustainable strategies used for controlling the disease. Knowledge about population structure and level of gene flow between neighbouring populations of the target vector is of high importance to develop appropriate strategies for implementing effective management programmes. Microsatellites are commonly used to identify population structure and assess dispersal of the target populations and have been developed for several tsetse species but were lacking for Glossina brevipalpis. In this study, we screened the genome of G. brevipalpis to search for suitable microsatellite markers and nine were found to be efficient enough to distinguish between different tsetse populations. The availability of these novel microsatellite loci will help to better understand the population biology of G. brevipalpis and to assess the level of gene flow between different populations. Such information will help with the development of appropriate strategies to implement the sterile insect technique (SIT) in the framework of an area-wide integrated pest management (AW-IPM) approach to manage tsetse populations and ultimately address the trypanosomoses problem in these targeted areas.

Benchmarking of human Y-chromosomal haplogroup classifiers with whole-genome and whole-exome sequence data

In anthropological, medical, and forensic studies, the nonrecombinant region of the human Y chromosome (NRY) enables accurate reconstruction of pedigree relationships and retrieval of ancestral information. Using high-throughput sequencing (HTS) data, we present a benchmarking analysis of command-line tools for NRY haplogroup classification. The evaluation was performed using paired Illumina data from whole-genome sequencing (WGS) and whole-exome sequencing (WES) experiments from 50 unrelated donors. Additionally, as a validation, we also used paired WGS/WES datasets of 54 individuals from the 1000 Genomes Project. Finally, we evaluated the tools on data from third-generation HTS obtained from a subset of donors and one reference sample. Our results show that WES, despite typically offering less genealogical resolution than WGS, is an effective method for determining the NRY haplogroup. Y-LineageTracker and Yleaf showed the highest accuracy for WGS data, classifying precisely 98% and 96% of the samples, respectively. Yleaf outperforms all benchmarked tools in the WES data, classifying approximately 90% of the samples. Yleaf, Y-LineageTracker, and pathPhynder can correctly classify most samples (88%) sequenced with third-generation HTS. As a result, Yleaf provides the best performance for applications that use WGS and WES. Overall, our study offers researchers with a guide that allows them to select the most appropriate tool to analyze the NRY region using both second- and third-generation HTS data.

A scaffolded and annotated reference genome of giant kelp (Macrocystis pyrifera)

Macrocystis pyrifera (giant kelp), is a brown macroalga of great ecological importance as a primary producer and structure-forming foundational species that provides habitat for hundreds of species. It has many commercial uses (e.g. source of alginate, fertilizer, cosmetics, feedstock). One of the limitations to exploiting giant kelp's economic potential and assisting in giant kelp conservation efforts is a lack of genomic tools like a high quality, contiguous reference genome with accurate gene annotations. Reference genomes attempt to capture the complete genomic sequence of an individual or species, and importantly provide a universal structure for comparison across a multitude of genetic experiments, both within and between species. We assembled the giant kelp genome of a haploid female gametophyte de novo using PacBio reads, then ordered contigs into chromosome level scaffolds using Hi-C. We found the giant kelp genome to be 537 MB, with a total of 35 scaffolds and 188 contigs. The assembly N50 is 13,669,674 with GC content of 50.37%. We assessed the genome completeness using BUSCO, and found giant kelp contained 94% of the BUSCO genes from the stramenopile clade. Annotation of the giant kelp genome revealed 25,919 genes. Additionally, we present genetic variation data based on 48 diploid giant kelp sporophytes from three different Southern California populations that confirms the population structure found in other studies of these populations. This work resulted in a high-quality giant kelp genome that greatly increases the genetic knowledge of this ecologically and economically vital species.

The genetic variability of grapevine Pinot gris virus (GPGV) in Australia

Grapevine Pinot gris virus (GPGV; genus Trichovirus in the family Betaflexiviridae) was detected in Australia in 2016, but its impact on the production of nursery material and fruit in Australia is still currently unknown. This study investigated the prevalence and genetic diversity of GPGV in Australia. GPGV was detected by reverse transcription-polymerase chain reaction (RT-PCR) in a range of rootstock, table and wine grape varieties from New South Wales, South Australia, and Victoria, with 473/2171 (21.8%) samples found to be infected. Genomes of 32 Australian GPGV isolates were sequenced and many of the isolates shared high nucleotide homology. Phylogenetic and haplotype analyses demonstrated that there were four distinct clades amongst the 32 Australian GPGV isolates and that there were likely to have been at least five separate introductions of the virus into Australia. Recombination and haplotype analysis indicate the emergence of new GPGV strains after introduction into Australia. When compared with 168 overseas GPGV isolates, the analyses suggest that the most likely origin of Australian GPGV isolates is from Europe. There was no correlation between specific GPGV genotypes and symptoms such as leaf mottling, leaf deformation, and shoot stunting, which were observed in some vineyards, and the virus was frequently found in symptomless grapevines.

Holistic view of the seascape dynamics and environment impact on macro-scale genetic connectivity of marine plankton populations

BACKGROUND

Plankton seascape genomics studies have revealed different trends from large-scale weak differentiation to microscale structures. Previous studies have underlined the influence of the environment and seascape on species differentiation and adaptation. However, these studies have generally focused on a few single species, sparse molecular markers, or local scales. Here, we investigated the genomic differentiation of plankton at the macro-scale in a holistic approach using Tara Oceans metagenomic data together with a reference-free computational method.

RESULTS

We reconstructed the FST-based genomic differentiation of 113 marine planktonic taxa occurring in the North and South Atlantic Oceans, Southern Ocean, and Mediterranean Sea. These taxa belong to various taxonomic clades spanning Metazoa, Chromista, Chlorophyta, Bacteria, and viruses. Globally, population genetic connectivity was significantly higher within oceanic basins and lower in bacteria and unicellular eukaryotes than in zooplankton. Using mixed linear models, we tested six abiotic factors influencing connectivity, including Lagrangian travel time, as proxies of oceanic current effects. We found that oceanic currents were the main population genetic connectivity drivers, together with temperature and salinity. Finally, we classified the 113 taxa into parameter-driven groups and showed that plankton taxa belonging to the same taxonomic rank such as phylum, class or order presented genomic differentiation driven by different environmental factors.

CONCLUSION

Our results validate the isolation-by-current hypothesis for a non-negligible proportion of taxa and highlight the role of other physicochemical parameters in large-scale plankton genetic connectivity. The reference-free approach used in this study offers a new systematic framework to analyse the population genomics of non-model and undocumented marine organisms from a large-scale and holistic point of view.

Absence of known knockdown resistance mutations but fixation of CYP337B3 was detected in field populations of Helicoverpa armigera across China

The cotton bollworm (Helicoverpa armigera) is a worldwide agricultural pest that infests many important crops. Pyrethroids targeting the voltage-gated sodium channel (VGSC) have been long used in the control of the cotton bollworm. Two amino acid substitutions (D1561V and E1565G) in H. armigera VGSC (HaVGSC) and the presence of a chimeric P450 gene (CYP337B3) have been documented to be associated with pyrethroid resistance. To understand the current occurrence of kdr mutations and the CYP337B3 gene in Chinese H. armigera populations, high-throughput amplicon sequencing was adopted to detect potential nucleotide variations in three fragments of the VGSC gene that cover 10 reported knockdown resistance (kdr) sites in insects, and gene-specific PCR was performed to examine the presence of CYP337B3 gene in H. armigera samples collected across China. The nucleotide variation analysis revealed a wealth of nucleotide variations in not only exons but also introns in the VGSC gene in Chinese H. armigera populations. However, neither previously reported kdr-conferring amino acid replacements nor other non-synonymous mutations were observed in a total of 1439 examined individuals. Population genetic analysis suggested that the H. armigera population in Nanchang, Jiangxi Province (JNC) had a moderate genetic differentiation from other populations, while no significant divergence was observed in other populations in northern and northwestern China. The CYP337B3 was present in all the examined individuals, indicating that CYP337B3 is extensively fixed in H. armigera populations across China. These results support that point mutations in VGSC are not a major factor involved in the current pyrethroid resistance in H. armigera. Instead, CYP337B3 plays a prevalent role in the development of resistance to pyrethroids in H. armigera.

Population genetic differentiation of the ubiquitous brooding coral Pocillopora acuta along Phuket Island reefs in the Andaman Sea, Thailand

BACKGROUND

The widespread Indo-Pacific coral species Pocillopora acuta Lamarck, 1816 displays varying levels of asexual versus sexual reproduction, with strong repercussions on genetic diversity, connectivity and genetic structuring within and among populations. For many geographic regions, baseline information on genetic diversity is still lacking, particularly in the Andaman Sea. The region suffered a massive heat-induced bleaching event in 2010 with high coral cover loss of branching coral species such as P. acuta. A subsequent bleaching in 2016, however, revealed a mild bleaching response in pocilloporids compared to other coral taxa in the region, suggesting that rare, heat tolerant genotypes had been selected by the 2010 bleaching event. In order to test whether this potential 'evolutionary rescue' event has led to a low genetic diversity, we conducted a population genetic survey covering a total of nine different P. acuta populations (336 individuals) along a 50 km coastal stretch around Phuket Island, Thailand. We used six microsatellite markers to assess genotypic diversity and to determine the prevalent mode of reproduction (i.e. sexual or asexual recruitment).

RESULTS

In contrast to other Indian Ocean P. acuta populations, the majority of corals in this study adopted a sexual reproduction mode (75% across all populations). At the same time, substantial regional gene flow was observed around Phuket Island with strong genetic differentiation as indicated by three genetic clusters that were separated by only a few kilometers. Patterns of isolation by distance over 0.7 - 40 km suggest small-scale genetic barriers, such as changing currents throughout each monsoonal season, potentially contributing to locally restricted dispersal of P. acuta larvae.

CONCLUSIONS

The occurrence of distinct genetic clusters within short coastal stretches suggests that the 2010 bleaching event has not led to extreme genetic impoverishment. While more in-depth genomic analyses are necessary to investigate changes in genetic diversity following extreme bleaching events, our results will help guide conservation efforts to maintain genetic diversity of a coral species that likely will be dominant in future, warmer Andaman Sea reefs.

Population genetics analysis of SLC3A1 and SLC7A9 revealed the etiology of cystine stone may be more than what our current genetic knowledge can explain

BACKGROUND

Cystine stone is a Mendelian genetic disease caused by SLC3A1 or SLC7A9. In this study, we aimed to estimate the genetic prevalence of cystine stones and compare it with the clinical prevalence to better understand the disease etiology.

METHODS

We analyzed genetic variants in the general population using the 1000 Genomes project and the Human Gene Mutation Database to extract all SLC3A1 and SLC7A9 pathogenic variants. All variants procured from both databases were intersected. Pathogenic allele frequency, carrier rate, and affected rate were calculated and estimated based on Hardy-Weinberg equilibrium.

RESULTS

We found that 9 unique SLC3A1 pathogenic variants were carried by 26 people and 5 unique SLC7A9 pathogenic variants were carried by 12 people, all of whom were heterozygote carriers. No homozygote, compoun d heterozygote, or double heterozygote was identified in the 1000 Genome database. Based on the Hardy-Weinberg equilibrium, the calculated genetic prevalence of cystine stone disease is 1 in 30,585.

CONCLUSION

The clinical prevalence of cystine stone has been previously reported as 1 in 7,000, a notably higher figure than the genetic prevalence of 1 in 30,585 calculated in this study. This suggests that the etiology of cystine stone is more complex than what our current genetic knowledge can explain. Possible factors that may contribute to this difference include novel causal genes, undiscovered pathogenic variants, alternative inheritance models, founder effects, epigenetic modifications, environmental factors, or other modifying factors. Further investigation is needed to fully understand the etiology of cystine stone.

Advances in the genetic characterization of the malaria vector, Anopheles funestus, and implications for improved surveillance and control

Anopheles mosquitoes present a major public health challenge in sub-Saharan Africa; notably, as vectors of malaria that kill over half a million people annually. In parts of the east and southern Africa region, one species in the Funestus group, Anopheles funestus, has established itself as an exceptionally dominant vector in some areas, it is responsible for more than 90% of all malaria transmission events. However, compared to other malaria vectors, the species is far less studied, partly due to difficulties in laboratory colonization and the unresolved aspects of its taxonomy and systematics. Control of An. funestus is also increasingly difficult because it has developed widespread resistance to public health insecticides. Fortunately, recent advances in molecular techniques are enabling greater insights into species identity, gene flow patterns, population structure, and the spread of resistance in mosquitoes. These advances and their potential applications are reviewed with a focus on four research themes relevant to the biology and control of An. funestus in Africa, namely: (i) the taxonomic characterization of different vector species within the Funestus group and their role in malaria transmission; (ii) insecticide resistance profile; (iii) population genetic diversity and gene flow, and (iv) applications of genetic technologies for surveillance and control. The research gaps and opportunities identified in this review will provide a basis for improving the surveillance and control of An. funestus and malaria transmission in Africa.

Lower viral evolutionary pressure under stable versus fluctuating conditions in subzero Arctic brines

BACKGROUND

Climate change threatens Earth's ice-based ecosystems which currently offer archives and eco-evolutionary experiments in the extreme. Arctic cryopeg brine (marine-derived, within permafrost) and sea ice brine, similar in subzero temperature and high salinity but different in temporal stability, are inhabited by microbes adapted to these extreme conditions. However, little is known about their viruses (community composition, diversity, interaction with hosts, or evolution) or how they might respond to geologically stable cryopeg versus fluctuating sea ice conditions.

RESULTS

We used long- and short-read viromics and metatranscriptomics to study viruses in Arctic cryopeg brine, sea ice brine, and underlying seawater, recovering 11,088 vOTUs (~species-level taxonomic unit), a 4.4-fold increase of known viruses in these brines. More specifically, the long-read-powered viromes doubled the number of longer (≥25 kb) vOTUs generated and recovered more hypervariable regions by >5-fold compared to short-read viromes. Distribution assessment, by comparing to known viruses in public databases, supported that cryopeg brine viruses were of marine origin yet distinct from either sea ice brine or seawater viruses, while 94% of sea ice brine viruses were also present in seawater. A virus-encoded, ecologically important exopolysaccharide biosynthesis gene was identified, and many viruses (~half of metatranscriptome-inferred "active" vOTUs) were predicted as actively infecting the dominant microbial genera Marinobacter and Polaribacter in cryopeg and sea ice brines, respectively. Evolutionarily, microdiversity (intra-species genetic variations) analyses suggested that viruses within the stable cryopeg brine were under significantly lower evolutionary pressures than those in the fluctuating sea ice environment, while many sea ice brine virus-tail genes were under positive selection, indicating virus-host co-evolutionary arms races.

CONCLUSIONS

Our results confirmed the benefits of long-read-powered viromics in understanding the environmental virosphere through significantly improved genomic recovery, expanding viral discovery and the potential for biological inference. Evidence of viruses actively infecting the dominant microbes in subzero brines and modulating host metabolism underscored the potential impact of viruses on these remote and underexplored extreme ecosystems. Microdiversity results shed light on different strategies viruses use to evolve and adapt when extreme conditions are stable versus fluctuating. Together, these findings verify the value of long-read-powered viromics and provide foundational data on viral evolution and virus-microbe interactions in Earth's destabilized and rapidly disappearing cryosphere. Video Abstract.

Isolation and characterization of 17 polymorphic microsatellite loci with tri- and tetra-nucleotide repeat motifs for Ayu (Plecoglossus altivelis) using next-generation sequencing

BACKGROUND

Ayu or sweetfish, Plecoglossus altivelis, an amphidromous fish ranging in the northwestern Pacific, is economically important inland fisheries and aquaculture resources. Genetic characterization of wild Ayu and derived culture seeds with competent molecular genetic markers is still insufficient for their sustainable use. Microsatellite DNA markers with larger repeat motifs (e.g. tri- and tetra-nucleotide motifs) are convenient and accurate compared with those having mono- and di-nucleotide motifs, but the latter motifs characterized most Ayu microsatellite markers developed previously.

METHODS AND RESULTS

Here, we isolated and characterized 17 polymorphic microsatellite DNA markers with tri- and tetra-nucleotide repeat motif using next-generation sequencing. Alleles per locus varied from 6 to 23. The observed and expected heterozygosities ranged from 0.542 to 1.000 and 0.709 to 0.951, respectively. Polymorphic information content (PIC) of 15 out of the 17 loci were high (≧ 0.700), suggesting them to be highly informative. Twelve of the 17 loci were used for preliminary assignment test among three collections, and successfully allocated the examined fish to the original populations.

CONCLUSION

The novel polymorphic microsatellite markers developed herein will be useful to examine the genetic diversity and population structure of wild Ayu and the effect of seed transplantation on native populations, providing a tool for conservation and sustainable adaptive management of this species.

The first maternal genetic study of hunter-gatherers from Vietnam

The current limitation of ancient DNA data from Vietnam led to the controversy surrounding the prehistory of people in this region. The combination of high heat and humidity damaged ancient bones that challenged the study of human evolution, especially when using DNA as study materials. So far, only 4 k years of history have been recorded despite the 65 k years of history of anatomically modern human occupations in Vietnam. Here we report, to our knowledge, the oldest mitogenomes of two hunter-gatherers from Vietnam. We extracted DNA from the femurs of two individuals aged 6.2 k cal BP from the Con Co Ngua (CCN) site in Thanh Hoa, Vietnam. This archeological site is the largest cemetery of the hunter-gatherer population in Southeast Asia (SEA) that was discovered, but their genetics have not been explored until the present. We indicated that the CCN haplotype belongs to a rare haplogroup that was not detected in any present-day Vietnamese individuals. Further matrilineal analysis on CCN mitogenomes showed a close relationship with ancient farmers and present-day populations in SEA. The mitogenomes of hunter-gatherers from Vietnam debate the "two layers" model of peopling history in SEA and provide an alternative solution for studying challenging ancient human samples from Vietnam.

INTERPRETING GENERATIVE ADVERSARIAL NETWORKS TO INFER NATURAL SELECTION FROM GENETIC DATA

Understanding natural selection in humans and other species is a major focus for the use of machine learning in population genetics. Existing methods rely on computationally intensive simulated training data. Unlike efficient neutral coalescent simulations for demographic inference, realistic simulations of selection typically requires slow forward simulations. Because there are many possible modes of selection, a high dimensional parameter space must be explored, with no guarantee that the simulated models are close to the real processes. Mismatches between simulated training data and real test data can lead to incorrect inference. Finally, it is difficult to interpret trained neural networks, leading to a lack of understanding about what features contribute to classification. Here we develop a new approach to detect selection that requires relatively few selection simulations during training. We use a Generative Adversarial Network (GAN) trained to simulate realistic neutral data. The resulting GAN consists of a generator (fitted demographic model) and a discriminator (convolutional neural network). For a genomic region, the discriminator predicts whether it is "real" or "fake" in the sense that it could have been simulated by the generator. As the "real" training data includes regions that experienced selection and the generator cannot produce such regions, regions with a high probability of being real are likely to have experienced selection. To further incentivize this behavior, we "fine-tune" the discriminator with a small number of selection simulations. We show that this approach has high power to detect selection in simulations, and that it finds regions under selection identified by state-of-the art population genetic methods in three human populations. Finally, we show how to interpret the trained networks by clustering hidden units of the discriminator based on their correlation patterns with known summary statistics. In summary, our approach is a novel, efficient, and powerful way to use machine learning to detect natural selection.

Development of microsatellite markers for the invasive mosquito Aedes koreicus

BACKGROUND

Aedes koreicus is a mosquito species native to East Asia which has recently invaded several countries in Europe. In Italy, this mosquito was first detected in the North-East in 2011 and is now widely distributed in the entire northern part of the country. The development of specific genetic markers, such as microsatellites, is necessary to uncover the dispersal routes of this mosquito from its native areas and, eventually, to plan future control interventions.

METHODS

Available raw sequences of genomic DNA of Ae. koreicus were screened in silico using BLASTn to identify possible microsatellite-containing sequences. Specific primer pairs were then designed, and their efficiency was determined through polymerase chain reaction (PCR) on 32 individuals of Ae. koreicus collected in Italy. PCR conditions were optimised in three multiplex reactions. Genotyping of individual mosquitoes was performed on both single and multiplex PCR reactions. Finally, analysis of intra-population variation was performed to assess the level of polymorphism of the markers.

RESULTS

Mosquito genotyping provided consistent results in both single and multiplex reactions. Out of the 31 microsatellite markers identified in the Ae. koreicus genome raw sequences, 11 were polymorphic in the examined mosquito samples.

CONCLUSIONS

The results show that the 11 microsatellite markers developed here hold potential for investigating the genetic structure of Ae. koreicus populations. These markers could thus represent a novel and useful tool to infer the routes of invasion of this mosquito species into Europe and other non-native areas.

Genome resequencing and transcriptome profiling reveal molecular evidence of tolerance to water deficit in barley

INTRODUCTION

Frequent climate change-induced drought events are detrimental environmental stresses affecting global crop production and ecosystem health. Several efforts have facilitated crop breeding for resilient varieties to counteract stress. However, progress is hampered due to the complexity of drought tolerance; a greater variety of novel genes are required across varying environments. Tibetan annual wild barley is a unique and precious germplasm that is well adapted to abiotic stress and can provide elite genes for crop improvement in drought tolerance.

OBJECTIVES

To identify the genetic basis and unique mechanisms for drought tolerance in Tibetan wild barley.

METHODS

Whole genome resequencing and comparative RNA-seq approaches were performed to identify candidate genes associated with drought tolerance via investigating the genetic diversity and transcriptional variation between cultivated and Tibetan wild barley. Bioinformatics, population genetics, and gene silencing were conducted to obtain insights into ecological adaptation in barley and functions of key genes.

RESULTS

Over 20 million genetic variants and a total of 15,361 significantly affected genes were identified in our dataset. Combined genomic, transcriptomic, evolutionary, and experimental analyses revealed 26 water deficit resilience-associated genes in the drought-tolerant wild barley XZ5 with unique genetic variants and expression patterns. Functional prediction revealed Tibetan wild barley employs effective regulators to activate various responsive pathways with novel genes, such as Zinc-Induced Facilitator-Like 2 (HvZIFL2) and Peroxidase 11 (HvPOD11), to adapt to water deficit conditions. Gene silencing and drought tolerance evaluation in a natural barley population demonstrated that HvZIFL2 and HvPOD11 positively regulate drought tolerance in barley.

CONCLUSION

Our findings reveal functional genes that have been selected across barley's complex history of domestication to thrive in water deficit environments. This will be useful for molecular breeding and provide new insights into drought-tolerance mechanisms in wild relatives of major cereal crops.

Blood Group Serology and "Race": Looking Back to Move Forward

Less than a decade after the discovery of the ABO antigens as a Mendelian inherited trait, blood group antigen frequencies were first used to define racial groups. This approach, known as seroanthropology, was the basis for collecting large amounts of blood group frequency data in different populations and was also sometimes used for racist purposes. Ultimately, population geneticists used these data to disprove race as a biological construct. Through understanding the history of seroanthropology, and recognizing the harms of its lingering presence, healthcare providers can better practice race-conscious, as opposed to race-based, transfusion medicine.

Interspecific variation and phylogenetic relationship between mangrove and non-mangrove species of a same family (Meliaceae)-insights from comparative analysis of complete chloroplast genome

The mahogany family, Meliaceae, contains 58 genera with only one mangrove genus: Xylocarpus. Two of the three species of the genus Xylocarpus are true mangroves (X. granatum and X. moluccensis), and one is a non-mangrove (X. rumphii). In order to resolve the phylogenetic relationship between the mangrove and non-mangrove species, we sequenced chloroplast genomes of these Xylocarpus species along with two non-mangrove species of the Meliaceae family (Carapa guianensis and Swietenia macrophylla) and compared the genome features and variations across the five species. The five Meliaceae species shared 130 genes (85 protein-coding genes, 37 tRNA, and eight rRNA) with identical direction and order, with a few variations in genes and intergenic spacers. The repetitive sequences identified in the rpl22 gene region only occurred in Xylocarpus, while the repetitive sequences in accD were found in X. moluccensis and X. rumphii. The TrnH-GUG and rpl32 gene regions and four non-coding gene regions showed high variabilities between X. granatum and the two non-mangrove species (S. macrophylla and C. guianensis). In addition, among the Xylocarpus species, only two genes (accD and clpP) showed positive selection. Carapa guianensis and S. macrophylla owned unique RNA editing sites. The above genes played an important role in acclimation to different stress factors like heat, low temperature, high UV light, and high salinity. Phylogenetic analysis with 22 species in the order Sapindales supported previous studies, which revealed that the non-mangrove species X. rumphii is closer to X. moluccensis than X. granatum. Overall, our results provided important insights into the variation of genetic structure and adaptation mechanism at interspecific (three Xylocarpus species) and intergeneric (mangrove and non-mangrove genera) levels.

Spatial and temporal population genetic analysis of Semaprochilodus insignis (Prochilodontidae), an overexploited fish from the Amazon basin

Background

Semaprochilodus insignis is a migratory fish of commercial and subsistence importance to communities in the Amazon. Despite the high intensity of exploitation, recent studies have not been carried out to assess the genetic status of its stocks.

Methods

This study is the first to estimate genetic diversity and to test the existence of spatial and temporal structuring of S. insignis through sequencing of the mtDNA control region (n = 241) and eight microsatellite loci (n = 180) of individuals sampled at 11 sites distributed in the Brazilian Amazon basin.

Results

Results for both markers were congruent, revealing a homogeneous genetic diversity in all the sampled locations, in addition to the absence of spatial and temporal genetic structure, indicating that the species forms a large panmictic population in the Brazilian Amazon.

Discussion

Although overfishing does not yet appear to have affected the levels of genetic variability of S. insignis, signals of reduction of the effective population size and a bottleneck provide an early alert to the effects of overfishing. Thus, the ever-decreasing populations may threaten S. insignis in the future. Therefore, it is hoped that the results of this study may contribute to the elaboration of management plans or any other measures that aim at the management and conservation of this species of great importance for the Amazon basin.

Lack of genetic structure in euryhaline Chelon labrosus from the estuaries under anthropic pressure in the Southern Bay of Biscay to the coastal waters of the Mediterranean Sea

Over the last decade, xenoestrogenic effects have been reported in populations of thicklip grey mullet Chelon labrosus from contaminated estuaries in the Bay of Biscay, resulting in intersex condition. To understand the level of gene flow in individuals of different Basque estuaries microsatellite markers were used to evaluate the population structure and connectivity of C. labrosus from estuaries of the Basque coast. 46 microsatellites were tested and 10 validated for the analysis of 204 individuals collected from 5 selected Basque estuaries and 2 outgroups in the Bay of Cadiz and Thermaic Gulf. The polymorphic microsatellites revealed 74 total alleles, 2-19 alleles per locus. The mean observed heterozygosity (0.49 ± 0.02) was lower than the expected one (0.53 ± 0.01). There was no evidence of genetic differentiation (F_ST = 0.0098, P = 0.0000) among individuals or sites. Bayesian clustering analysis revealed a single population in all sampled locations. The results of this study indicate widespread genetic homogeneity and panmixia of C. labrosus across the current sampling areas spanning the Atlantic and Mediterranean basins. The hypothesis of panmixia could therefore be well supported so individuals inhabiting estuaries with high prevalence of intersex condition should be considered as members of the same single genetic group as those inhabiting adjacent estuaries without incidence of xenoestrogenicity.

Comprehensive analysis of geographic and breed-purpose influences on genetic diversity and inherited disease risk in the Doberman dog breed

BACKGROUND

Publicly available phenotype data and genotyping array data from two citizen science projects: "Doberman Health Surveys" and "The Doberman Diversity Project" were analyzed to explore relative homozygosity, diversity, and disorder risk according to geographical locale and breeding purpose in the Doberman.

RESULTS

From the phenotypic data cohort, life expectancy of a Doberman at birth is 9.1 years. The leading causes of death were heart disease (accounting for 28% of deaths) and cancers (collectively accounting for 14% of deaths). By genotyping, the world Doberman population exists as four major cohorts (European exhibition-bred, Americas exhibition-bred, European work, Americas pet/informal). Considering the entire Doberman population, four genomic regions longer than 500 Kb are fixed in 90% or more of 3,226 dogs included in this study. The four fixed regions reside on two autosomal chromosomes: CFA3:0.8-2.3 Mb (1.55 Mb); CFA3: 57.9-59.8 Mb (1.8 Mb); CFA31:0-1.2 Mb (1.2 Mb); and CFA31:4.80-6.47 Mb (1.67 Mb). Using public variant call files including variants for eight Doberman pinschers, we observed 30 potentially functional alternate variants that were evolutionarily diverged relative to the wider sequenced dog population within the four strongly homozygous chromosomal regions. Effective population size (Ne) is a statistical measure of breed diversity at the time of sampling that approximates the number of unique individuals. The major identified sub-populations of Dobermans demonstrated Ne in the range 70-236. The mean level of inbreeding in the Doberman breed is 40% as calculated by the number of array variants in runs of homozygosity divided by the assayed genome size (excluding the X chromosome). The lowest observed level of inbreeding in the Dobermans assayed was 15% in animals that were first generation mixes of European and USA bred Dobermans. Array variant analysis shows that inter-crossing between European and USA-bred Dobermans has capacity to re-introduce variation at many loci that are strongly homozygous.

CONCLUSIONS

We conclude that efforts to improve breed diversity first should focus on regions with the highest fixation levels, but managers must ensure that mutation loads are not worsened by increasing the frequencies of rarer haplotypes in the identified regions. The analysis of global data identified regions of strong fixation that might impact known disorder risks in the breed. Plausible gene candidates for future analysis of the genetic basis of cardiac disease and cancer were identified in the analysis.

Development of a multiplex panel with 36 insertion/deletion markers (InDel) for individual identification

In recent years, the insertion/deletion (InDel) polymorphism has become a preferred genetic marker in forensic genetics due to its low mutation rates and small amplicon sizes. In this study, a 36-InDelplex identification panel, consisting of autosomal 34 InDel loci, 1 Y InDel locus, and amelogenin, was developed, and gene frequencies in the Turkish population were determined. The loci of the InDel panel with global minimum allele frequencies (MAF) ≥ 0.4 were selected from the 1000 Genomes Project Phase 3 data. The amplicon sizes of the loci were designed in the range of 69-252 bp. In the validation study of the developed panel, analysis threshold, dynamic range, sensitivity, stochastic threshold, inhibitor tolerance, and reproducibility parameters were studied by following the Scientific Working Group on DNA Analysis Methods (SWGDAM) guidelines. The sensitivity studies indicated that complete and reliable InDel profiles could be obtained with 0.25 ng of DNA. A population study was evaluated using 250 samples from Turkey. The mean observed heterozygosity ratio (Ho) of all loci was 0.48. The combined discrimination power (CPD) is 0.999999999990867 and the combined exclusion probability (CPE) was 0.9930. The population comparison was also made using Turkish and the five major populations from the 1000 Genomes Phase 3 populations' data (Africa, Europe, East Asia, South Asia, and America). In conclusion, the results showed that the 36-InDelplex panel is a reliable, sensitive, and accurate system that is suitable for human identification and population genetics purposes.

[Preliminary study on the mechanism underlying the ecological isolation of Oncomelania hupensis populations in Changde City]

OBJECTIVE

To investigate ecological isolation between Oncomelania hupensis snail populations in hilly regions and marshland and lake regions in Yuanjiang valley, Changde City, Hunan Province, and to unravel its underlying mechanisms.

METHODS

Taoyuan County, Shimen County, Linli County and Lixian County in Changde City were selected as snail sampling sites in hilly regions, and Lixian County, Jinshi City, West Lake Administration District, Hanshou County and Dingcheng District were selected as snail sampling sites in marshland and lake areas. Cytochrome C oxidase 1 (cox 1) gene was amplified in snail samples and sequenced. The genetic sequences of O. hupensis snails were aligned using the software MEGA 11, and the haplotypes of O. hupensis snails were determined using the software DNASP 5.10.01. The phylogenetic tree was generated using Bayesian inference with the software MrBayes 3.2, and analysis of molecular variance (AMOVA) was performed to analyze the source of genetic divergence and estimate the genetic divergence index (F_ST) among snail populations with the software Arlequin 3.5.2.2. The genetic barrier among 11 O. hupensis snail populations was estimated using the Monmonier algorithm of adegenet toolkit in R package. The settings with "land in winter and water in summer" in the Yuanjian River section were divided into two categories according to the upstream and downstream, and the areas with "land in winter and water in summer" in the upstream and downstream were transformed into raster data, and then loaded into the software Fragstats 4 for analysis of landscape indicators. The trends in changes of digital elevation were extracted from the Yuanjiang River section based on the digital elevation model, and made three-dimensional visualization using the R package.

RESULTS

The mitochondrial cox 1 gene were amplified in 165 O. hupensis snais from 11 sampling sites and sequenced, and a total of 152 valid gene sequences were obtained, with 46 haplotypes or 9 populations determined. No haplotype was shared in snails between Taoyuan County and Dingcheng District and Hanshou County along the downstream of the Yuanjiang River. The total area of settings with "land in winter and water in summer" was 617.66 hm² in the upsteram of the Yuanjiang River, which consisted of 473 patches, with each patch measuring 1.31 hm², the largest area index of 0.735 2, the landscape division index of 0.999 9, and the landscape shape index of 45.293 7. The total area of settings with "land in winter and water in summer" was 9 956.92 hm² in the downstream of the Yuanjiang River, which consisted of 771 patches, with each patch measuring 12.91 hm², the largest area index of 97.839 9, the landscape division index of 0.042 7, and the landscape shape index of 7.249 6. The area of settings with "land in winter and water in summer" was much larger in the downstream than that in the upstream of the Yuanjiang River, and the stronger landscape connectivity and non-remarkable alteration of riverbed elevation provided suitable habitats for snail breeding.

CONCLUSIONS

The hydrological and environmental characteristics of the upstream of the Yuanjiang River restrain the breeding and spread of O. hupensis, resulting in ecological isolation between Oncomelania hupensis in Taoyuan County and those in the downstream of Yuanjiang River.

Counterintuitive properties of evolutionary measures: A stochastic process study in cyclic population structures with periodic environments

Local environmental interactions are a major factor in determining the success of a new mutant in structured populations. Spatial variations in the concentration of genotype-specific resources change the fitness of competing strategies locally and thus can drastically change the outcome of evolutionary processes in unintuitive ways. The question is how such local environmental variations in network population structures change the condition for selection and fixation probability of an advantageous (or deleterious) mutant. We consider linear graph structures and focus on the case where resources have a spatial periodic pattern. This is the simplest model with two parameters, length scale and fitness scales, representing heterogeneity. We calculate fixation probability and fixation times for a constant population birth-death process as fitness heterogeneity and period vary. Fixation probability is affected by not only the level of fitness heterogeneity but also spatial scale of resources variations set by period of distribution T. We identify conditions for which a previously a deleterious mutant (in a uniform environment) becomes beneficial as fitness heterogeneity is increased. We observe cases where the fixation probability of both mutant and resident types are more than their neutral value, 1/N, simultaneously. This coincides with exponential increase in time to fixation which points to potential coexistence of resident and mutant types. Finally, we discuss the effect of the 'fitness shift' where the fitness function of two types has a phase difference. We observe significant increases (or decreases) in the fixation probability of the mutant as a result of such phase shift.

Resequencing of Rosa rugosa accessions revealed the history of population dynamics, breed origin, and domestication pathways

BACKGROUND

Rosa rugosa is a shrub that originated in China and has economic and ecological value. However, during the development of R. rugosa, the genetic background was chaotic, and the genetic structure among different wild populations was unclear, as well as wild and cultivated accessions. Here, we report whole-genome resequencing of wild and cultivated R. rugosa accessions.

RESULTS

A total of 19,041,284 SNPs were identified in 188 R. rugosa accessions and 3 R. chinensis accessions by resequencing. Population genetic analysis revealed that cultivated and wild groups were separated very early. All R. rugosa accessions were divided into 8 categories based on genetic structure: (1) Weihai, Yantai, and Liaoning category, (2) Jilin category, and (3) Hammonasset category (above three are wild); (4) traditional varieties, (5) hybrids between R. rugosa and R. chinensis, (6) Zizhi Rose, (7) Kushui Rose, (8) hybrids between R. rugosa and R. multiflora. We found that the heterozygosity and genetic diversity of wild accessions were generally lower than those of cultivated individuals. The genes that were selected during cultivation were identified, and it was found that these genes were mainly related to environmental adaptation and growth.

CONCLUSIONS

The Jilin population was the oldest population and later migrated to Liaoning and then migrated to Yantai and Weihai by sea regression in the Bohai Basin. The Hammonasset naturalized population probably originated from the Jilin population and then experienced separate differentiation. The long-term asexual reproduction pattern of R. rugosa decreased genetic diversity in the wild population. During R. rugosa cultivation, the ancestors of the Jilin population were involved in breeding traditional varieties, after which almost no wild individuals were engaged in breeding. However, in recent decades, cross breeding of R. rugosa started the utilization of wild germplasms. In comparison, some other species play important roles in variety formation. Few genes related to economic traits were selected, suggesting no directional domestication in the R. rugosa cultivation process.

Phylogenomic analyses reveal incongruences between divergence times and fossil records of freshwater snails in East Asia

Fossils provide important insight into our understanding of phylogenetic history by serving as calibration points for divergence time estimation. However, uncertainties in the fossil record due to parallel evolution and convergent evolution can critically affect estimates of node ages. Here, we compare and contrast estimates of phylogenetic divergence with geologic and fossil history for two freshwater snail genera of the family Viviparidae in East Asia (Cipangopaludina and Margarya). Cipangopaludina species are commonly widely distributed species in East Asia, but extant Margarya species are endemic to the ancient lakes in Yunnan, China. According to some previous studies, parallel evolution or convergent evolution of shell morphology has occurred in the family several times which may affect divergence time estimation using fossil records. In this study, we used SNP data derived from ddRAD-seq loci to investigate population demographic history of both genera. Our results show a common pattern of lake endemic lineages diversifying from widely distributed lineages in the Miocene, and multiple colonization to a single ancient lake occurred in the Pleistocene. Our results indicate substantial incongruence among estimated phylogenomic divergence times, some fossil records, and formation ages of ancient lakes. These findings suggest some fossil records may be misidentified in these groups and highlight the need to carefully evaluate geological evidence and fossil records when using these for divergence time estimation.

Priors, population sizes, and power in genome-wide hypothesis tests

BACKGROUND

Genome-wide tests, including genome-wide association studies (GWAS) of germ-line genetic variants, driver tests of cancer somatic mutations, and transcriptome-wide association tests of RNAseq data, carry a high multiple testing burden. This burden can be overcome by enrolling larger cohorts or alleviated by using prior biological knowledge to favor some hypotheses over others. Here we compare these two methods in terms of their abilities to boost the power of hypothesis testing.

RESULTS

We provide a quantitative estimate for progress in cohort sizes and present a theoretical analysis of the power of oracular hard priors: priors that select a subset of hypotheses for testing, with an oracular guarantee that all true positives are within the tested subset. This theory demonstrates that for GWAS, strong priors that limit testing to 100-1000 genes provide less power than typical annual 20-40% increases in cohort sizes. Furthermore, non-oracular priors that exclude even a small fraction of true positives from the tested set can perform worse than not using a prior at all.

CONCLUSION

Our results provide a theoretical explanation for the continued dominance of simple, unbiased univariate hypothesis tests for GWAS: if a statistical question can be answered by larger cohort sizes, it should be answered by larger cohort sizes rather than by more complicated biased methods involving priors. We suggest that priors are better suited for non-statistical aspects of biology, such as pathway structure and causality, that are not yet easily captured by standard hypothesis tests.

High-resolution genotyping of 58 STRs in 635 Northern Han Chinese with MiSeq FGx ® Forensic Genomics System

Sequence polymorphisms were characterized at 27 autosomal STRs (A-STRs), 7 X chromosomal STRs (X-STRs), and 24 Y chromosomal STRs (Y-STRs) in 635 Northern Han Chinese with the ForenSeq DNA Signature Prep Kit on the MiSeq FGx Forensic Genomics System. Since repeat region (RR) and flanking region (FR) variation can be detected by massively parallel sequencing (MPS), the increase in the number of unique alleles and the average of gene diversity was 78.18% and 3.51% between sequence and length, respectively. A total of 74 novel RR variants were identified at 33 STRs compared with STRSeq and previous studies, and 13 FR variants (rs1770275883, rs2053373277, rs2082557941, rs1925525766, rs1926380862, rs1569322793, rs2051848492, rs2051848696, rs2016239814, rs2053269960, rs2044518192, rs2044536444, and rs2089968964) were first submitted to dbSNP. Also, 99.94% of alleles were concordant between the ForenSeq DNA Signature Prep Kit and commercial CE kits. Discordance resulted from the low performance at D22S1045 and occasionally at DYS392, flanking region deletions at D7S820 and DXS10074, and the strict alignment algorithm at DXS7132. Null alleles at DYS505 and DYS448 and multialleles at DYS387S1a/b, DYS385a/b, DYS448, DYS505, DXS7132, and HPRTB were validated with other MPS and CE kits. Thus, a high-resolution sequence-based (SB) and length-based (LB) allele frequencies dataset from Northern Han Chinese has been established already. As expected, forensic parameters increased significantly on combined power of discrimination (PD) and combined power of exclusion (PE) at A-STRs, mildly on combined PD and combined mean exclusion chance (MEC) at X-STRs, and barely on discrimination capacity (DC) at Y-STRs. Additionally, MiSeq FGx quality metrics and MPS performance were evaluated in this study, which presented the high-quality of the dataset at 20 consecutive runs, such as ≥ 60% bases with a quality score of 20 or higher (%≥ Q20), > 60% of effective reads, > 2000 × of depth of coverage (DoC), ≥ 60% of allele coverage ratio (ACR) or heterozygote balance, ≥ 70% of inter-locus balance, and ≤ 0.4 of the absolute value of observed minus expected heterozygosity (|H_exp - H_obs|). In conclusion, MiSeq FGx can help us generate a high-resolution and high-quality dataset for human identification and population genetic studies.

Exploring the forensic effectiveness and population genetic differentiation by self-constructed 41 multi-InDel panel in Yunnan Zhuang group

Yunnan is one of the main residences of the Zhuang group which is one of the 55 ethnic minorities in China. At present, there are relatively few researches on population genetics and forensic science of the Yunnan Zhuang group. Therefore, this study used a self-constructed panel containing 41 multi-InDel markers to analyze the genetic polymorphisms of 173 individuals from Yunnan Zhuang group. The results indicated that these 41 multi-InDels in Yunnan Zhuang group were highly polymorphic markers expect for three markers. The cumulative match probability and combined exclusion probability values of the 40 multi-InDels (MI38 marker was excluded) were 8.0671E-26 and 0.9999995959, respectively. In addition, population genetic analyses were performed on genotyping data of 41 multi-InDel markers among the Yunnan Zhuang and 26 reference populations, revealing that the Yunnan Zhuang group was genetically close to the five populations in East Asia. According to the STRUCTURE analysis, the Yunnan Zhuang group presented similar ancestral compositions to the five populations from East Asia, and when the K value was three, the five intercontinental populations showed their different genetic structures. In conclusion, the 41 multi-InDel markers could be used as an effective tool for individual identification and paternity testing of the Zhuang group in Yunnan province, as well as for their ancestry information inference studies.

Isolation and characterization of novel microsatellite loci for the Eastern Pacific marine sponge Mycale cecilia by Illumina MiSeq sequencing

BACKGROUND

Mycale cecilia is an abundant Eastern Tropical Pacific sponge living in a wide variety of habitats, including coral reefs where it may directly interact with corals. It is also known to possess secondary metabolites of pharmacological value. These aspects highlight the importance of having a better understanding of its biology, and genetic and population diversity.

METHODS AND RESULTS

In the present study, we isolated and characterized twelve novel microsatellite loci by Illumina MiSeq sequencing. The loci were tested in 30 specimens collected from two coral reef localities (La Paz, Baja California Sur and Isabel Island, Nayarit) from the Mexican Pacific using M13(-21) labeling. All loci were polymorphic, with two to nine alleles per locus. Expected heterozygosities varied from 0.616 to 0.901. Eleven loci were tested and successfully amplified in M. microsigmatosa from the Gulf of Mexico.

CONCLUSION

Here we report the first microsatellite loci developed for a sponge species from the Eastern Pacific coast. These molecular markers will be used for population genetic studies of M. cecilia, and potentially in other congeneric species; particularly in vulnerable marine areas that require protection, such as coral reefs.

Genetic architecture of human brain evolution

Comparative studies of hominids have long sought to identify mutational events that shaped the evolution of the human nervous system. However, functional genetic differences are outnumbered by millions of nearly neutral mutations, and the developmental mechanisms underlying human nervous system specializations are difficult to model and incompletely understood. Candidate-gene studies have attempted to map select human-specific genetic differences to neurodevelopmental functions, but it remains unclear how to contextualize the relative effects of genes that are investigated independently. Considering these limitations, we discuss scalable approaches for probing the functional contributions of human-specific genetic differences. We propose that a systems-level view will enable a more quantitative and integrative understanding of the genetic, molecular and cellular underpinnings of human nervous system evolution.

Development of twenty-one novel microsatellite loci for Gila topminnow, Poeciliopsis occidentalis occidentalis

BACKGROUND

Gila topminnow (Poeciliopsis occidentalis occidentalis) was once highly abundant throughout the Lower Colorado River Basin of the southwestern United States. However, this Sonoran Desert endemic suffered extreme population declines over the past century because of habitat degradation and nonnative species introductions. Much of the prior conservation genetic work conducted on the species relied upon a small number of microsatellite loci; many exhibiting low variability in extant populations. Consequently, there was a need for additional microsatellite loci to provide high-resolution delimitation of populations for conservation purposes.

METHODS AND RESULTS

Paired-end Illumina sequencing was utilized to screen the Gila topminnow genome for novel microsatellite loci. We identified 21 novel loci that exhibited no deviations from expectations of genetic equilibrium, and cross-amplified in Yaqui topminnow (P. o. sonoriensis). These loci were amplified from 401 samples representing eight populations of Gila topminnow and Yaqui topminnow. Although diversity was low for all populations (observed heterozygosity = 0.12 to 0.45), these novel markers provided ample power to identify population of origin for each individual in Bayesian assignment tests.

CONCLUSIONS

This novel set of microsatellite loci provide a useful genetic tool to assess population genetic parameters of the endangered Gila topminnow and delineate populations for identifying conservation priorities. The cross-amplification of these loci in Yaqui topminnow shows promise for application to other Poeciliopsis species of Mexico and Central America.

Comprehensive molecular characterization of complete mitogenome assemblies of 33 Eimeria isolates infecting domestic chickens

BACKGROUND

Coccidiosis caused by Eimeria is one of the most severe chicken diseases and poses a great economic threat to the poultry industry. Understanding the evolutionary biology of chicken Eimeria parasites underpins development of new interactions toward the improved prevention and control of this poultry disease.

METHODS

We presented an evolutionary blueprint of chicken coccidia by genetically characterizing complete mitogenome assemblies of 33 isolates representing all seven known Eimeria species infecting chickens in China. Further genome- and gene-level phylogenies were also achieved to better understand the evolutionary relationships of these chicken Eimeria at the species level.

RESULTS

33 mitogenomes of chicken eimerian parasites ranged from 6148 bp to 6480 bp in size and encoded typical mitochondrial compositions of apicomplexan parasites including three protein-coding genes (PCGs), seven fragmented small subunit (SSU) and 12/13 fragmented large subunit (LSU) rRNAs. Comparative genomics provided an evolutionary scenario for the genetic diversity of PCGs-cytochrome c oxidase subunits 1 and 3 (cox1 and cox3) and cytochrome b (cytb); all were under purifying selection with cox1 and cox3 being the lowest and highest evolutionary rates, respectively. Genome-wide phylogenies classified the 33 Eimeria isolates into seven subgroups, and furthermore Eimeria tenella and Eimeria necatrix were determined to be more closely related to each other than to the other eight congenic species. Single/concatenated mitochondrial protein gene-based phylogenies supported cox1 as the genetic marker for evolutionary and phylogenetic studies for avain coccidia.

CONCLUSIONS

To our knowledge, these are the first population-level mitogenomic data on the genus Eimeria, and its comprehensive molecular characterization provides valuable resources for systematic, population genetic and evolutionary biological studies of apicomplexan parasites in poultry.

Population genomics of helminth parasites

Next generation sequencing technologies have facilitated a shift from a few targeted loci in population genetic studies to whole genome approaches. Here, we review the types of questions and inferences regarding the population biology and evolution of parasitic helminths being addressed within the field of population genomics. Topics include parabiome, hybridization, population structure, loci under selection and linkage mapping. We highlight various advances, and note the current trends in the field, particularly a focus on human-related parasites despite the inherent biodiversity of helminth species. We conclude by advocating for a broader application of population genomics to reflect the taxonomic and life history breadth displayed by helminth parasites. As such, our basic knowledge about helminth population biology and evolution would be enhanced while the diversity of helminths in itself would facilitate population genomic comparative studies to address broader ecological and evolutionary concepts.

Fine-scale spatial distribution of deltamethrin resistance and population structure of Anopheles funestus and Anopheles arabiensis populations in Southern Mozambique

BACKGROUND

Insecticide resistance in malaria vectors can be spatially highly heterogeneous, yet population structure analyses frequently find relatively high levels of gene flow among mosquito populations. Few studies have contemporaneously assessed phenotypic, genotypic and population structure analysis on mosquito populations and none at fine geographical scales. In this study, genetic diversity, population structure, and insecticide resistance profiles of Anopheles funestus and Anopheles arabiensis were examined across mosquito populations from and within neighbouring villages.

METHODS

Mosquitoes were collected from 11 towns in southern Mozambique, as well as from different neighbourhoods within the town of Palmeira, during the peak malaria transmission season in 2016. CDC bottle bioassay and PCR assays were performed with Anopheles mosquitoes at each site to determine phenotypic and molecular insecticide resistance profiles, respectively. Microsatellite analysis was conducted on a subsample of mosquitoes to estimate genetic diversity and population structure.

RESULTS

Phenotypic insecticide resistance to deltamethrin was observed in An. funestus sensu stricto (s.s.) throughout the area, though a high level of mortality variation was seen. However, 98% of An. funestus s.s. were CYP6P9a homozygous resistant. An. arabiensis was phenotypically susceptible to deltamethrin and 99% were kdr homozygous susceptible. Both Anopheles species exhibited high allelic richness and heterozygosity. Significant deviations from Hardy-Weinberg equilibrium were observed, and high linkage disequilibrium was seen for An. funestus s.s., supporting population subdivision. However, the F_ST values were low for both anophelines (- 0.00457 to 0.04213), N_m values were high (9.4-71.8 migrants per generation), AMOVA results showed almost 100% genetic variation among and within individuals, and Structure analysis showed no clustering of An. funestus s.s. and An. arabiensis populations. These results suggest high gene flow among mosquito populations.

CONCLUSION

Despite a relatively high level of phenotypic variation in the An. funestus population, molecular analysis shows the population is admixed. These data indicate that CYP6P9a resistance markers do not capture all phenotypic variation in the area, but also that resistance genes of high impact are likely to easily spread in the area. Conversely, other strategies, such as transgenic mosquito release programmes will likely not face challenges in this locality.

Allele frequencies and selection coefficients in locally adapted populations

Understanding the role of natural selection in driving evolutionary change requires accurate estimates of the strength of selection acting at the genetic level in the wild. This is challenging to achieve but may be easier in the case of populations in migration-selection balance. When two populations are at equilibrium under migration-selection balance, there exist loci whose alleles are selected different ways in the two populations. Such loci can be identified from genome sequencing by their high values of F_ST. This raises the question of what is the strength of selection on locally-adaptive alleles. To answer this question we analyse a 1-locus 2-allele model of a population distributed between two niches. We show by simulation of selected cases that the outputs from finite-population models are essentially the same as those from deterministic infinite-population models. We then derive theory for the infinite-population model showing the dependence of selection coefficients on equilibrium allele frequencies, migration rates, dominance and relative population sizes in the two niches. An Excel spreadsheet is provided for the calculation of selection coefficients and their approximate standard errors from observed values of population parameters. We illustrate our results with a worked example, with graphs showing the dependence of selection coefficients on equilibrium allele frequencies, and graphs showing how F_ST depends on the selection coefficients acting on the alleles at a locus. Given the extent of recent progress in ecological genomics, we hope our methods may help those studying migration-selection balance to quantify the advantages conferred by adaptive genes.

Development and evaluation of a novel panel containing 188 microhaplotypes for 2nd-degree kinship testing in the Hebei Han population

Distant kinship identification is one of the critical problems in forensic genetics. As a new type of genetic marker defined and discussed in the last decade, the microhaplotype (MH) has drawn much attention in such identification owing to its specific advantages to traditional short tandem repeat (STR) or single nucleotide polymorphism (SNP) markers. In this study, MH markers were screened step by step from the 1000 Genomes Project database, and a novel multiplex panel containing 188 MHs (in which 181 are reported the first time, while 1 was reported in a previous study and the other 6 have partial overlaps with known markers) was constructed for application in 2nd- and 3rd-degree kinship identification. Along with the construction, a novel MH nomenclature was proposed, in which the SNP position information they contained was taken into account to eliminate the possibility that the same locus was named differently interlaboratory. After a series of evaluations, the panel was shown to have good sequencing accuracy, high sensitivity, species specificity, and resistance to anti-PCR inhibitors or degradation. Population data of the 188 MHs were calculated based on the genetic information of 221 unrelated Hebei Han individuals, and the effective number of alleles (Ae) ranged from 2.0925 to 8.2634 (with an average of 2.9267). For the whole system, the cumulative matching probability (CMP), the cumulative power of exclusion in paternity testing of duos (CPEduo) and that of trios (CPEtrio) reached 2.8422 × 10^-137, 1-1.3109 × 10^-21, and 1-2.8975 × 10^-39, respectively, indicating that this panel was satisfactory for individual identification and paternity testing. Then, the efficiency of the 188 MHs in 2nd- and 3rd-degree kinship testing was studied based on 30 extended families consisting of 179 2nd-degree and 121 3rd-degree relatives, as well as simulations of 0.5 million pairs of those two kinships. The results showed that clear opinions would be given in 83.36% of 2nd-degree identifications with a false rate less than 10^-5, when the confirming and excluding thresholds of cumulative likelihood ratio (CLR) were set as 10⁴ and 10^-4, respectively. This panel is still not sufficient to solve the problem of 3rd-degree kinship identification alone, and approximately 300 or 870 MH loci would be needed in 2nd- or 3rd-degree kinship identification, respectively, to achieve a system efficiency not less than 0.99 with such a threshold set; such necessary numbers would be used only as a reference in further research.

FST estimates of 94 populations in China based on STR markers

The proper assessment of DNA evidence in cases of personal identification is a recurring theme in forensics. It is common practice to evaluate the strength of DNA evidence using the likelihood ratio (LR). The accurate use of population allele frequencies is a crucial problem in LR calculation. Allele frequency differences among different populations could be estimated by the F_ST values. Thus, F_ST would also affect LR values by correcting the allele frequencies. In this study, Chinese population allele frequency data were selected from population reports published in Chinese and English journals. The population-specific F_ST values of each population, the overall F_ST values of each province, each region, and the whole country, and the locus-specific F_ST values of different loci were calculated. The LRs using different allele frequencies and different F_ST values were compared based on the combination of simulated genotypes. As a result, the F_ST values of 94 populations, 19 provinces, 7 regions, and the whole country were obtained. The LR was overestimated using allele frequencies of the combined population containing multiple populations rather than using allele frequencies of a population, and the LRs after F_ST correction were lower than those without correction. Conclusively, the correction in conjunction with corresponding F_ST values can make the LRs more accurate and reasonable.

Genetic Structure of Aedes (Stegomyia) albopictus Populations in Russia

Background

Aedes (Stegomyia) albopictus was found for the first time in 2011 on the Black Sea coast in Russia, and during 2011-2019, the species expanded over two climate zones Cfa and Csa.

Methods

Here, we studied the sequence diversity of the mitochondrial cytochrome c oxidase I (COI) gene, 1317-1433bp in length. In total, 131 specimens of Ae. albopictus sampled from 21 locations in Russia and Abkhazia were examined.

Results

Two of the six identified mitochondrial haplotypes were detected for the first time. Four COI haplotypes were shared by at least two studied local populations. The most prevalent H1 and H2 haplotypes dominated in all the sampled localities in the Cfa zone. The H3 haplotype was prevalent in the Csa zone. Other haplotypes were rare. Phylogenetic analyses, spatial isolation and limited gene flow revealed that the samples from the Csa zone differed significantly from those from the Cfa zone.

Conclusion

Two spatially isolated genetic lineages exist in Ae. albopictus population in southern region of Russia. One lineage obtained on the seacoast and inland (in valleys of the Caucasus Mountains and steppe zone) is widely distributed worldwide including Mediterranean populations. This confirms the hypothesis that the emergence of Ae. albopictus population in southern region of Russia may be associated with the terrestrial spread of mosquitoes from the well-established European population due to human activity. The other lineage, discovered in Novorossiysk, a maritime port, is similar to Ae. albopictus from the USA and Japan, suggesting the independent introduction of these mosquitoes.

The race to understand immunopathology in COVID-19: Perspectives on the impact of quantitative approaches to understand within-host interactions

The COVID-19 pandemic has revealed the need for the increased integration of modelling and data analysis to public health, experimental, and clinical studies. Throughout the first two years of the pandemic, there has been a concerted effort to improve our understanding of the within-host immune response to the SARS-CoV-2 virus to provide better predictions of COVID-19 severity, treatment and vaccine development questions, and insights into viral evolution and the impacts of variants on immunopathology. Here we provide perspectives on what has been accomplished using quantitative methods, including predictive modelling, population genetics, machine learning, and dimensionality reduction techniques, in the first 26 months of the COVID-19 pandemic approaches, and where we go from here to improve our responses to this and future pandemics.

Phylogeography and population structure of the global, wide host-range hybrid pathogen Phytophthora × cambivora

Invasive, exotic plant pathogens pose a major threat to native and agricultural ecosystems. Phytophthora × cambivora is an invasive, destructive pathogen of forest and fruit trees causing severe damage worldwide to chestnuts (Castanea), apricots, peaches, plums, almonds and cherries (Prunus), apples (Malus), oaks (Quercus), and beech (Fagus). It was one of the first damaging invasive Phytophthora species to be introduced to Europe and North America, although its origin is unknown. We determined its population genetic history in Europe, North and South America, Australia and East Asia (mainly Japan) using genotyping-by-sequencing. Populations in Europe and Australia appear clonal, those in North America are highly clonal yet show some degree of sexual reproduction, and those in East Asia are partially sexual. Two clonal lineages, each of opposite mating type, and a hybrid lineage derived from these two lineages, dominated the populations in Europe and were predominantly found on fagaceous forest hosts (Castanea, Quercus, Fagus). Isolates from fruit trees (Prunus and Malus) belonged to a separate lineage found in Australia, North America, Europe and East Asia, indicating the disease on fruit trees could be caused by a distinct lineage of P. × cambivora, which may potentially be a separate sister species and has likely been moved with live plants. The highest genetic diversity was found in Japan, suggesting that East Asia is the centre of origin of the pathogen. Further surveys in unsampled, temperate regions of East Asia are needed to more precisely identify the location and range of the centre of diversity.

Genomic-based microsatellite development for Ternstroemia (Pentaphylacaceae) and transferability to other Ericales

BACKGROUND

The genus Ternstroemia is associated with the vulnerable tropical montane cloud forest in Mexico and with other relevant vegetation types worldwide. It contains threatened and pharmacologically important species and has taxonomic issues regarding its species limits. This study describes 38 microsatellite markers generated using a genomic-based approach.

METHODS AND RESULTS

We tested 23 of these markers in a natural population of Ternstroemia lineata. These markers are highly polymorphic (all loci polymorphic with 3-14 alleles per locus and expected heterozygosity between 0.202 and 0.908), most of them (19 out of 23) are in Hardy-Weinberg Equilibrium and free of null alleles (18 out of 23). Also we found no evidence of linkage among them. Finally, we tested the transferability to six other American species of Ternstroemia, two other Pentaphylacaceae species, and four species from different families within the order Ericales.

CONCLUSIONS

These molecular resources are promising tools to investigate genetic diversity loss and as barcodes for ethnopharmacological applications and species delimitation in the family Pentaphylacaceae and some Ericales, among other applications.

Population-based genomic study of Plasmodium vivax malaria in seven Brazilian states and across South America

Background

Brazil is a unique and understudied setting for malaria, with complex foci of transmission associated with human and environmental conditions. An understanding of the population genomic diversity of P. vivax parasites across Brazil can support malaria control strategies.

Methods

Through whole genome sequencing of P. vivax isolates across 7 Brazilian states, we use population genomic approaches to compare genetic diversity within country (n = 123), continent (6 countries, n = 315) and globally (26 countries, n = 885).

Findings

We confirm that South American isolates are distinct, have more ancestral populations than the other global regions, with differentiating mutations in genes under selective pressure linked to antimalarial drugs (pvmdr1, pvdhfr-ts) and mosquito vectors (pvcrmp3, pvP45/48, pvP47). We demonstrate Brazil as a distinct parasite population, with signals of selection including ABC transporter (PvABCI3) and PHIST exported proteins.

Interpretation

Brazil has a complex population structure, with evidence of P. simium infections and Amazonian parasites separating into multiple clusters. Overall, our work provides the first Brazil-wide analysis of P. vivax population structure and identifies important mutations, which can inform future research and control measures.

Funding

AI is funded by an MRC LiD PhD studentship. TGC is funded by the Medical Research Council (Grant no. MR/M01360X/1, MR/N010469/1, MR/R025576/1, MR/R020973/1 and MR/X005895/1). SC is funded by Medical Research Council UK grants (MR/M01360X/1, MR/R025576/1, MR/R020973/1 and MR/X005895/1) and Bloomsbury SET (ref. CCF17-7779). FN is funded by The Shloklo Malaria Research Unit - part of the Mahidol Oxford Research Unit, supported by the Wellcome Trust (Grant no. 220211). ARSB is funded by São Paulo Research Foundation - FAPESP (Grant no. 2002/09546-1). RLDM is funded by Brazilian National Council for Scientific and Technological Development - CNPq (Grant no. 302353/2003-8 and 471605/2011-5); CRFM is funded by FAPESP (Grant no. 2020/06747-4) and CNPq (Grant no. 302917/2019-5 and 408636/2018-1); JGD is funded by FAPESP fellowships (2016/13465-0 and 2019/12068-5) and CNPq (Grant no. 409216/2018-6).