Forensic and genetic characterizations of diverse southern Thai populations based on 15 autosomal STRs

Southern Thailand is home to various populations; the Moklen, Moken and Urak Lawoi’ sea nomads and Maniq negrito are the minority, while the southern Thai groups (Buddhist and Muslim) are the majority. Although previous studies have generated forensic STR dataset for major groups, such data of the southern Thai minority have not been included; here we generated a regional forensic database of southern Thailand. We newly genotyped common 15 autosomal STRs in 184 unrelated southern Thais, including all minorities and majorities. When combined with previously published data of major southern Thais, this provides a total of 334 southern Thai samples. The forensic parameter results show appropriate values for personal identification and paternity testing; the probability of excluding paternity is 0.99999622, and the combined discrimination power is 0.999999999999999. Probably driven by genetic drift and/or isolation with small census size, we found genetic distinction of the Maniq and sea nomads from the major groups, which were closer to the Malay and central Thais than the other Thai groups. The allelic frequency results can strength the regional forensic database in southern Thailand and also provide useful information for anthropological perspective.

Demographic history and conservation genomics of caribou (Rangifer tarandus) in Québec

The loss of genetic diversity is a challenge many species are facing, and genomics is a potential tool that can inform and prioritize decision making. Caribou populations have experienced significant recent declines throughout Québec, Canada, and some are considered threatened or endangered. We calculated the ancestral and contemporary patterns of genomic diversity of five caribou populations and applied a comparative framework to assess the interplay between demography and genomic diversity. We calculated a caribou specific mutation rate, μ, by extracting orthologous genes from related ungulates. Whole genome re-sequencing was completed on 67 caribou and genotype likelihoods were estimated. We calculated nucleotide diversity, θπ and estimated the coalescent or ancestral Ne, which ranged from 12,030 to 15,513. When compared to the census size, NC, the endangered Gaspésie Mountain caribou population had the highest Ne:NC ratio which is consistent with recent work suggesting high ancestral Ne:NC is of conservation concern. These ratios were highly correlated with genomic signatures (i.e. Tajima’s D) and explicit demographic model parameters. Values of contemporary Ne, estimated from linkage-disequilibrium, ranged from 11 to 162, with Gaspésie having among the highest contemporary Ne:NC ratio. Importantly, conservation genetics theory would predict this population to be of less concern based on this ratio. Of note, F varied only slightly between populations, and runs of homozygosity were not abundant in the genome. Our study highlights how genomic patterns are nuanced and misleading if viewed only through a contemporary lens; a holistic view should integrate ancestral Ne and Tajima’s D into conservation decisions.

PSXII-8 Mitochondrial DNA diversity within domestic reindeer in Russia

Domestic reindeer in Russia are a valuable resource of vital importance to the physical and cultural survival of the Northern indigenous minority. During the last decades, the mitochondrial (mt) genetic markers have been widely used as a molecular tool to investigate genetic structure and diversity of livestock species. Here we aimed at the assessing the mtDNA diversity of the domestic reindeer inhabiting the area from the Kola Peninsula in the west to the Chukotka region in the east. A complete cytochrome b (cytb) sequences (1,140 bp) from representatives of six populations, including Nenets (NEN, n = 16), Evenk (EVK, n = 12), Even (EVN, n = 6), Chukotka (CHU, n = 6), Chukotka-Khargin (CHUKH, n = 6) and Tuva (TUVA, n = 6) were obtained. Sequences’ alignment was conducted using MUSCLE algorithm in R package msa. In total, 34 haplotypes were identified. Median-joining network, constructed in PopART 1.7, revealed three major groups of haplotypes: the first one joined the samples of all the populations, the second one included NEN, EVN and CHUKH, and the third group was presented by the one sample of CHU. AMOVA, calculated in Arlequin 3.5.2.2, showed that only 9.58% of molecular variance could be explained by the differences between populations and 90.42% - within populations. Genetic diversity parameters calculated in DnaSP 6.12.03, demonstrated that average number of nucleotide differences (K) was highest in CHUKH (28.333) and EVN (27.409) and lowest in TUVA (4.533) and EVK (5.400). Nucleotide diversity (Pi) was 0.01238±0.00559, 0.00474±0.00091, 0.02404±0.00453, 0.01281±0.00464, 0.02485±0.00744, and 0.00398±0.00110 for NEN, EVK, EVN, CHU, CHUKH and TUVA, respectively. Our study demonstrated the lack of clear genetic structure of the studied reindeer populations in relation to cytb sequence. The level of genetic diversity was associated with census size and was lowest in the smallest Tuva population. This study was supported by RSF-21-16-00071 and Russian Ministry of Science and Higher Education-0445-2019-0024.

PSXIV-26 Selection footprints in Russian red cattle identified by linkage disequilibrium blocks based on SNP data

Among the variety of cattle breeds in Russia, the Russian red dual-purpose cattle breeds have great importance because of their ability to produce high milk yields as well as to provide excellent milk quality. The low census size of the Russian red cattle breeds requires development of programs for conservation of their biodiversity. Our study aimed to investigate selection footprints in Russian red cattle breeds, using high values of linkage disequilibrium (LD) in SNP haplotype blocks as indicators. For finding such LD blocks, we used the genotypes (≈35K SNPs) of Red Gorbatov (RGB, n = 26), Bestuzhev (BST,n = 27), and Suksun (SKS,n = 17) breeds, as well as Red Holstein (RH,n = 16) as an outgroup. Quality control and LD calculations for different distances were performed in Plink 1.90. Top 0.01% SNP pairs by LD value (0.9≤r2< 1.0) were selected for further analysis. The effective population size derived from LD patterns was estimated using SNeP tool. Comparison of LD values for 70 kb interval between breeds and chromosomes by MANOVA pairwise testing significantly distinguished RH/RGB and BST/SKS breeds (P < 0.05-0.001). LD values among chromosomes were 0.195–0.287 for RH, 0.194–0.272 for RGB, 0.172–0.237 for BST, and 0.157–0.217 for SKS. The SKS and BST breeds had higher Ne values (84 and 113, respectively) compared to RH (63) and RGB (79). Selection footprints by LD blocks in Russian red cattle genome covered several relevant genes on BTA1 (EPHA6,DGKG), BTA2 (LRP1B,THSD7B,STAT1), BTA5 (CPM,BAIAP2L2), BTA9 (TRDN,UTRN), BTA10 (KCNN2,CAPN3), BTA11 (SH3RF3,RABGAP1,RALGPS1), BTA14 (ZNF16,ARHGAP39,TOX,DGAT1), and BTA19 (MYH10), BTA22 (FHIT). Detected genes were found to be responsible for milk fat and protein contents, fatty acid composition, somatic cells score, fertility, feet and legs, and udder conformation traits. Our results can be useful for developing the breeding and conservation programs of the Russian red cattle genetic resources. The study was funded by RFBR within project No. 20-516-00020

Estimating movement rates between Eurasian and North American birds that are vectors of avian influenza (AI)

Avian influenza (AI) is an emerging zoonotic disease that will likely be involved in future pandemics. Because waterbird movements are difficult to quantify, determining the host-specific risk of Eurasian-origin AI movements into North America is challenging. We estimated relative rates of movements, based on long-term evolutionary averages of gene flow, between Eurasian and North American waterbird populations to obtain bidirectional baseline rates of the intercontinental movements of these AI hosts. We used population genomics and coalescent-based demographic models to obtain these gene-flow-based movement estimates. Inferred rates of movement between these populations varies greatly among species. Within dabbling ducks, gene flow, relative to effective population size, varies from ~3-24 individuals/generation between Eurasian and American wigeons (Mareca penelope — M. americana) to ~100-300 individuals/generation between continental populations of northern pintails (Anas acuta). These are evolutionary long-term averages and provide a solid foundation for understanding the relative risks of each of these host species in potential intercontinental AI movements. We scale these values to census size for evaluation in that context. In addition to being AI hosts, many of these species are also important in the subsistence diets of Alaskans, increasing the risk of direct bird-to-human exposure to Eurasian-origin AI virus. We contrast species-specific rates of intercontinental movements with the importance of each species in Alaskan diets to understand the relative risk of these taxa to humans. Greater scaup (Aythya marila), mallard (Anas platyrhynchos), and northern pintail (Anas acuta) were the top three species presenting the highest risks for intercontinental AI movement both within the natural system and through exposure to subsistence hunters. These directly comparable, species-based intercontinental movement rates and relative risk rankings should help in modeling, monitoring, and mitigating the impacts of intercontinental host and AI movements.

Quantifying the relationship between genetic diversity and population size suggests natural selection cannot explain Lewontin's paradox

Neutral theory predicts that genetic diversity increases with population size, yet observed levels of diversity across metazoans vary only two orders of magnitude while population sizes vary over several. This unexpectedly narrow range of diversity is known as Lewontin’s Paradox of Variation (1974). While some have suggested selection constrains diversity, tests of this hypothesis seem to fall short. Here, I revisit Lewontin’s Paradox to assess whether current models of linked selection are capable of reducing diversity to this extent. To quantify the discrepancy between pairwise diversity and census population sizes across species, I combine previously-published estimates of pairwise diversity from 172 metazoan taxa with newly derived estimates of census sizes. Using phylogenetic comparative methods, I show this relationship is significant accounting for phylogeny, but with high phylogenetic signal and evidence that some lineages experience shifts in the evolutionary rate of diversity deep in the past. Additionally, I find a negative relationship between recombination map length and census size, suggesting abundant species have less recombination and experience greater reductions in diversity due to linked selection. However, I show that even assuming strong and abundant selection, models of linked selection are unlikely to explain the observed relationship between diversity and census sizes across species.

Trunk perimeter correlates with genetic bottleneck intensity and the level of genetic diversity in populations of Taxus baccata L

Key message Taxus baccata remnants established recently tend to contribute less to the species’ overall genetic variation than historical populations because they are subjected to a greater impact of the founder effect and genetic isolation. As tree trunk perimeter is a rough indicator of genetic variation in a population, this measure should be considered in conservation programs. Context Genetic variation within Taxus baccata (L.) populations is not associated with the current census size but correlates well with the effective size, suggesting that genetic drift intensity reflects variation in demographic histories. Aims We hypothesize that recently established populations are subjected to greater bottleneck than old remnants. Using the mean trunk perimeter as a surrogate of tree age, we test whether the demographic history and genetic variation are associated with the mean tree age. Methods Using 18 microsatellite markers, we analyze the genetic diversity and demographic history of 11 yew populations in Poland to assess the relationship between the mean trunk perimeter and the inferred genetic parameters. Results Populations reveal significant differences in levels of genetic variation and in the intensity and time of genetic bottleneck. After excluding an apparent outlier, the genetic variation is significantly greater while the bottleneck intensity lower in populations with a greater perimeter. Conclusion Due to continuous species decline and increasing fragmentation, the non-uniform contribution of yew remnants to the overall genetic variation tends to decrease together with the mean tree age. Germplasm collections for the species should take into account tree perimeter as a rough indicator of the genetic variation of a population.

Allee Effects Determine Pollinator Services at Local and Geographic Scales

Climatic factors have attracted much attention in the study of species’ distributions, while little is known about the role of biotic interactions. Here, we tested for variation in pollinator service across the distribution of a plant species, and evaluated the driving mechanisms. We monitored insect pollinators using time-lapse cameras in populations of North American Arabidopsis lyrata from the southern to the northern range limit. We spotted 67 pollinating insect taxa, indicating that this plant-pollinator network is a generalist system. Pollinator service increased with latitude. Higher pollinator visitation was correlated with the richness of other flowering plants and with plant census size, which was largest in northern populations. Furthermore, pollinator service reached a maximum at intermediate local flower density. Synthesis: This study indicates that pollination service underlies Allee effects on a local and species’ range scale, and that plant populations at range limits receive only marginal pollination service if they are small.