Estimating F-statistics

Selection signatures for high altitude adaptation in livestock: A review

High altitude adapted livestock species (cattle, yak, goat, sheep, and horse) has critical role in the human socioeconomic sphere and acts as good source of animal source products including milk, meat, and leather, among other things. These species sustain production and reproduction even in harsh environments on account of adaptation resulting from continued evolution of beneficial traits. Selection pressure leads to various adaptive strategies in livestock whose footprints are evident at the different genomic sites as the "Selection Signature". Scrutiny of these signatures provides us crucial insight into the evolutionary process and domestication of livestock adapted to diverse climatic conditions. These signatures have the potential to change the sphere of animal breeding and further usher the selection programmes in right direction. Technological revolution and recent strides made in genomic studies has opened the routes for the identification of selection signatures. Numerous statistical approaches and bioinformatics tools have been developed to detect the selection signature. Consequently, studies across years have identified candidate genes under selection region found associated with numerous traits which have a say in adaptation to high-altitude environment. This makes it pertinent to have a better understanding about the selection signature, the ways to identify and how to utilize them for betterment of livestock populations as well as farmers. This review takes a closer look into the general concept, various methodologies, and bioinformatics tools commonly employed in selection signature studies and summarize the results of recent selection signature studies related to high-altitude adaptation in various livestock species. This review will serve as an informative and useful insight for researchers and students in the field of animal breeding and evolutionary biology.

Linking the spatial and genomic structure of adaptive potential for conservation management: a review

We unified the recent literature with the goal to contribute to the discussion on how genetic diversity might best be conserved. We argue that this decision will be guided by how genomic variation is distributed among manageable populations (i.e., its spatial structure), the degree to which adaptive potential is best predicted by variation across the entire genome or the subset of that variation that is identified as putatively adaptive (i.e., its genomic structure), and whether we are managing species as single entities or as collections of diversifying lineages. The distribution of genetic variation and our ultimate goal will have practical implications for on-the-ground management. If adaptive variation is largely polygenic or responsive to change, its spatial structure might be broadly governed by the forces determining genome-wide variation (linked selection, drift, and gene flow), making measurement and prioritization straightforward. If we are managing species as single entities, then population-level prioritization schemes are possible so as to maximize future pooled genetic variation. We outline one such scheme based on the popular Shapley value from cooperative game theory that considers the relative genetic contribution of a population to an unknown future collection of populations.

MeHA: A Computational Framework in Revealing the Genetic Basis of Animal Mental Health Traits Under an Intensive Farming System-A Case Study in Pigs

Intensively farmed animals such as pigs inevitably experience a certain degree of psychological stress, which leads to a reduction in production performance. Mental health traits are currently difficult to measure, resulting in a gap in understanding their genetic basis. To address this challenge, we propose a computational framework called mental health of animals (MeHA), capable of revealing genes related to animal mental health traits. Using MeHA, we identified 109 candidate genes associated with pig mental health and discovered their intricate connections with critical functions, such as memory, cognition, and neural development, which are essential components of mental health and cognitive performance. Importantly, our findings provide evidence of the potential impact of these genes on economically important traits, including meat quality and piglet survival. This research underscores the importance of genetic studies in enhancing our understanding of animal behavior and cognition, as well as promoting agricultural practices. By applying our approach to study the genetic basis of mental health in pigs as a case, we confirmed that our framework is an effective way to reveal genetic factors affecting animal mental health traits, which contributes to animal welfare and has potential implications for understanding human mental disorders.

Notes from the taxonomic disaster zone: Evolutionary drivers of intractable species boundaries in an Australian lizard clade (Scincidae: Ctenotus)

Genomic-scale datasets, sophisticated analytical techniques, and conceptual advances have disproportionately failed to resolve species boundaries in some groups relative to others. To understand the processes that underlie taxonomic intractability, we dissect the speciation history of an Australian lizard clade that arguably represents a "worst-case" scenario for species delimitation within vertebrates: the Ctenotus inornatus species group, a clade beset with decoupled genetic and phenotypic breaks, uncertain geographic ranges, and parallelism in purportedly diagnostic morphological characters. We sampled hundreds of localities to generate a genomic perspective on population divergence, structure, and admixture. Our results revealed rampant paraphyly of nominate taxa in the group, with lineages that are either morphologically cryptic or polytypic. Isolation-by-distance patterns reflect spatially continuous differentiation among certain pairs of putative species, yet genetic and geographic distances are decoupled in other pairs. Comparisons of mitochondrial and nuclear gene trees, tests of nuclear introgression, and historical demographic modelling identified gene flow between divergent candidate species. Levels of admixture are decoupled from phylogenetic relatedness; gene flow is often higher between sympatric species than between parapatric populations of the same species. Such idiosyncratic patterns of introgression contribute to species boundaries that are fuzzy while also varying in fuzziness. Our results suggest that "taxonomic disaster zones" like the C. inornatus species group result from spatial variation in the porosity of species boundaries and the resulting patterns of genetic and phenotypic variation. This study raises questions about the origin and persistence of hybridizing species and highlights the unique insights provided by taxa that have long eluded straightforward taxonomic categorization.

Disentangling effects of dispersal, environment and anthropogenic barriers on functional connectivity in aquatic systems

Disentangling the roles of structural landscape factors and animal movement behaviour can present challenges for practitioners managing landscapes to maintain functional connectivity and achieve conservation goals. We used a landscape genetics approach to combine robust demographic, behavioural and genetic datasets with spatially explicit simulations to evaluate the effects of anthropogenic barriers (dams, culverts) and natural landscape resistance (gradient, elevation) affecting dispersal behaviour, genetic connectivity and genetic structure in a resident population of Westslope Cutthroat Trout (Oncorhynchus clarkii lewisi). Analyses based on 10 years of sampling effort revealed a pattern of restricted dispersal, and population genetics identified discrete population clusters between distal tributaries and the mainstem stream and no structure within the mainstem stream. Demogenetic simulations demonstrated that, for this population, the effects of existing anthropogenic barriers on population structure are redundant with effects of restricted dispersal associated with the underlying environmental resistance. Our approach provides an example of how extensive field sampling combined with landscape genetics can be incorporated into spatially explicit simulation modelling to explore how, together, movement ecology and landscape resistance can be used to inform decisions around restoration and connectivity.

Association of Genetic Profile with Muscle Mass Gain and Muscle Injury Prevention in Professional Football Players after Creatine Supplementation

BACKGROUND

In recent years, the study of creatine supplementation in professional athletes has been of great interest. However, the genetics involved in response to supplementation is unknown. The aim of this study was to analyse, for the first time, the relationship between muscle performance-related genes and the risk of an increased body mass index (BMI) and muscle mass and a decrease in fat mass in professional football players after creatine supplementation.

METHODS

For this longitudinal study, one hundred and sixty-one men's professional football players were recruited. The polymorphisms ACE I/D, ACTN3 c.1729C>T, AMPD1 c.34C>T, CKM c.*800A>G, and MLCK (c.49C>T and c.37885C>A) were genotyped using Single-Nucleotide Primer Extension (SNPE). To assess the combined impact of these six polymorphisms, a total genotype score (TGS) was calculated. The creatine supplementation protocol consisted of 20 g/day of creatine monohydrate for 5 days (loading dose) and 3-5 g/day for 7 weeks (maintenance dose). Anthropometric characteristics (body mass index (BMI), fat, and muscle mass) were recorded before and after the creatine supplementation protocol. Characteristics of non-contact muscle injuries during the 2022/2023 season were classified according to a consensus statement for injury recording. The results showed that the allelic frequencies of ACE and AMPD1 differed between responders and non-responders in muscle mass increase (all p < 0.05). Players with a TGS exceeding 54.16 a.u. had an odds ratio (OR) of 2.985 (95%CI: 1.560-5.711; p = 0.001) for muscle mass increase. By contrast, those with a TGS below 54.16 a.u. had an OR of 9.385 (95%CI: 4.535-19.425; p < 0.001) for suffering non-contact muscle injuries during the season.

CONCLUSIONS

The increase in BMI and muscle mass in response to creatine supplementation in professional football players was influenced by a TGS derived from the combination of favourable genotypes linked to muscle performance. The CC genotype and C allele of AMPD1 were particularly associated with a higher likelihood of muscle mass increase under creatine supplementation in this group of professional football players.

The Association of Genetic Markers Involved in Muscle Performance Responding to Lactate Levels during Physical Exercise Therapy by Nordic Walking in Patients with Long COVID Syndrome: A Nonrandomized Controlled Pilot Study

Several genetic markers have shown associations with muscle performance and physical abilities, but the response to exercise therapy is still unknown. The aim of this study was to test the response of patients with long COVID through an aerobic physical therapy strategy by the Nordic walking program and how several genetic polymorphisms involved in muscle performance influence physical capabilities. Using a nonrandomized controlled pilot study, 29 patients who previously suffered from COVID-19 (long COVID = 13, COVID-19 = 16) performed a Nordic walking exercise therapy program for 12 sessions. The influence of the ACE (rs4646994), ACTN3 (rs1815739), AMPD1 (rs17602729), CKM (rs8111989), and MLCK (rs2849757 and rs2700352) polymorphisms, genotyped by using single nucleotide primer extension (SNPE) in lactic acid concentration was established with a three-way ANOVA (group × genotype × sessions). For ACE polymorphism, the main effect was lactic acid (p = 0.019). In ACTN3 polymorphism, there were no main effects of lactic acid, group, or genotype. However, the posthoc analysis revealed that, in comparison with nonlong COVID, long COVID increased lactic acid concentrations in Nordic walking sessions in CT and TT genotypes (all p < 0.05). For AMPD1 polymorphism, there were main effects of lactic acid, group, or genotype and lactic acid × genotype or lactic acid × group × genotype interactions (all p < 0.05). The posthoc analysis revealed that, in comparison with nonlong COVID, long COVID increased lactic acid concentrations in Nordic walking sessions in CC and CT genotypes (all p < 0.05). Physical therapy strategy through Nordic walking enhanced physical capabilities during aerobic exercise in post-COVID19 patients with different genotypes in ACTN3 c.1729C>T and AMPD1 c.34C>T polymorphisms. These findings suggest that individuals who reported long COVID who presumably exercised less beforehand appeared to be less able to exercise, based on lactate levels, and the effect of aerobic physical exercise enhanced physical capabilities conditioned by several genetic markers in long COVID patients.

Association between glioma and neurodegenerative diseases risk: a two-sample bi-directional Mendelian randomization analysis

Background

Earlier observational studies have demonstrated a correlation between glioma and the risk of neurodegenerative diseases (NDs), but the causality and direction of their associations remain unclear. The objective of this study was to ascertain the causal link between glioma and NDs using Mendelian randomization (MR) methodology.

Methods

Genome-wide association study (GWAS) data were used in a two-sample bi-directional MR analysis. From the largest meta-analysis GWAS, encompassing 18,169 controls and 12,488 cases, summary statistics data on gliomas was extracted. Summarized statistics for NDs, including Alzheimer's disease (AD), multiple sclerosis (MS), amyotrophic lateral sclerosis (ALS) and Parkinson's disease (PD) were obtained from the GWAS of European ancestry. Inverse variance weighted (IVW) method was elected as the core MR approach with weighted median (WM) method and MR-Egger method as complementary methods. In addition, sensitivity analyses were performed. A Bonferroni correction was used to correct the results.

Results

Genetically predicted glioma had been related to decreased risk of AD. Specifically, for all glioma (IVW: OR = 0.93, 95% CI = 0.90-0.96, p = 4.88 × 10-6) and glioblastoma (GBM) (IVW: OR = 0.93, 95% CI = 0.91-0.95, p = 5.11 × 10-9). We also found that genetically predicted all glioma has a suggestive causative association with MS (IVW: OR = 0.90, 95% CI = 0.81-1.00, p = 0.045). There was no evidence of causal association between glioma and ALS or PD. According to the results of reverse MR analysis, no discernible causal connection of NDs was found on glioma. Sensitivity analyses validated the robustness of the above associations.

Conclusion

We report evidence in support of potential causal associations of different glioma subtypes with AD and MS. More studies are required to uncover the underlying mechanisms of these findings.

Weakened tanning ability is an important mechanism for evolutionary skin lightening in East Asians

The evolution of light-skin pigmentation among Eurasians is considered as an adaptation to the high-latitude environments. East Asians are ideal populations for studying skin color evolution because of the complex environment of East Asia. Here, we report a strong selection signal for the pigmentation gene phenylalanine hydroxylase (PAH) in light-skinned Han Chinese individuals. The intron mutation rs10778203 in PAH is enriched in East Asians and is significantly associated with skin color of the back of the hand in Han Chinese males (P < 0.05). In vitro luciferase and transcription factor binding assays show that the ancestral allele of rs10778203 could bind to SMAD2 and has a significant enhancer activity for PAH. However, the derived T allele (the major allele in East Asians) of rs10778203 decreases the binding activity of transcription factors and enhancer activity. Meanwhile, the derived T allele of rs10778203 shows a weaker ultraviolet radiation response in A375 cells and zebrafish embryos. Furthermore, rs10778203 decreases melanin production in transgenic zebrafish embryos after ultraviolet B (UVB) treatment. Collectively, PAH is a potential pigmentation gene that regulates skin tanning ability. Natural selection has enriched the adaptive allele, resulting in weakened tanning ability in East Asians, suggesting a unique genetic mechanism for evolutionary skin lightening in East Asians.

Multiomics analyses of Jining Grey goat and Boer goat reveal genomic regions associated with fatty acid and amino acid metabolism and muscle development

OBJECTIVE

Jining Grey goat is a local Chinese goat breed that is well known for its high fertility and excellent meat quality but shows low meat production performance. Numerous studies have focused on revealing the genetic mechanism of its high fertility, but its highlighting meat quality and muscle growth mechanism still need to be studied.

METHODS

In this research, an integrative analysis of the genomics and transcriptomics of Jining Grey goats compared with Boer goats was performed to identify candidate genes and pathways related to the mechanisms of meat quality and muscle development.

RESULTS

Our results overlap among five genes (ABHD2, FN1, PGM2L1, PRKAG3, RAVER2) and detected a set of candidate genes associated with fatty acid metabolism (PRKAG3, HADHB, FASN, ACADM), amino acid metabolism (KMT2C, PLOD3, NSD2, SETDB1, STT3B, MAN1A2, BCKDHB, NAT8L, P4HA3) and muscle development (MSTN, PPARGC1A, ANKRD2). Several pathways have also been detected, such as the FoxO signaling pathway and Apelin signaling pathway that play roles in lipid metabolism, lysine degradation, N-glycan biosynthesis, valine, leucine and isoleucine degradation that involving with amino acid metabolism.

CONCLUSION

The comparative genomic and transcriptomic analysis of Jining Grey goat and Boer goat revealed the mechanisms underlying the meat quality and meat productive performance of goats. These results provide valuable information for future breeding of goats.

Low frequency of the wild-type freezing-tolerance LsCBF7 allele among lettuce population suggests a negative selection during domestication and breeding

KEY MESSAGE

Sustainable winter production in lettuce requires freezing tolerant varieties. This study identified a wild-type allele of LsCBF7 that could contribute to freezing tolerance improvement in lettuce. Lettuce is one of the most consumed vegetables globally. While ideally grown in 13-21 °C, its cultivation extends into winter in milder climates. However, occasional freezing temperatures can significantly reduce yields. Therefore, the development of freezing-tolerant lettuce varieties has become a long-term goal of lettuce breeding programs. Despite its significance, our understanding of freezing tolerance in lettuce remains limited. Plants have evolved a coping mechanism against freezing, known as cold acclimation, whereby they can increase freezing tolerance when pre-exposed to low nonfreezing temperatures. The CBF pathway is well-known for its central role in cold acclimation. Previously, we identified 14 CBF genes in lettuce and discovered that one of them, LsCBF7, had a loss-of-function mutation. In this study, we uncovered that accessions from colder regions carried the wild-type allele of LsCBF7 and this allele likely contributed to increased freezing tolerance, with 14% of the lettuce population carrying this allele. Interestingly, in wild lettuce (L. serriola) that is considered a progenitor of cultivated lettuce, this wild-type allele was much more common, with a frequency of 90%. This finding suggests that this wild-type allele may have undergone negative selection during the domestication or breeding of lettuce. Our data strongly indicate that this allele could be linked to early bolting, an undesirable trait in lettuce, which may have driven the negative selection. While this wild-type allele shows promise for improving freezing tolerance in lettuce, it is crucial to decouple it from the early bolting trait to fully harness its potential in lettuce breeding.

Reconstructing the ancestral gene pool to uncover the origins and genetic links of Hmong-Mien speakers

BACKGROUND

Hmong-Mien (HM) speakers are linguistically related and live primarily in China, but little is known about their ancestral origins or the evolutionary mechanism shaping their genomic diversity. In particular, the lack of whole-genome sequencing data on the Yao population has prevented a full investigation of the origins and evolutionary history of HM speakers. As such, their origins are debatable.

RESULTS

Here, we made a deep sequencing effort of 80 Yao genomes, and our analysis together with 28 East Asian populations and 968 ancient Asian genomes suggested that there is a strong genetic basis for the formation of the HM language family. We estimated that the most recent common ancestor dates to 5800 years ago, while the genetic divergence between the HM and Tai-Kadai speakers was estimated to be 8200 years ago. We proposed that HM speakers originated from the Yangtze River Basin and spread with agricultural civilization. We identified highly differentiated variants between HM and Han Chinese, in particular, a deafness-related missense variant (rs72474224) in the GJB2 gene is in a higher frequency in HM speakers than in others.

CONCLUSIONS

Our results indicated complex gene flow and medically relevant variants involved in the HM speakers' evolution history.

Sequence-based population structure, relatedness, and inbreeding estimates for forensic autosomal STR markers

Population data have become available for sequence data to aid forensic investigations and prepare the forensic community in the move towards implementing NGS methods. This comes with a need for updated population genetic parameters estimates to allow DNA evidence evaluations using sequence data. Initial work has been done on a small sample and here we expand this work by providing estimates of population structure and relatedness for autosomal STR data generated by sequencing technologies. We also discuss the effect of inbreeding on forensic calculations and discuss why the use of genotypic-based estimates may be preferred over allelic-based estimates.

Assessing the genetic integrity of sugarcane germplasm in the USDA-ARS National Plant Germplasm System collection using single-dose SNP markers

The World Collection of Sugarcane and Related Grasses, maintained at the USDA-ARS in Miami, FL, is one of the largest sugarcane germplasm repositories in the world. However, the genetic integrity of the Saccharum spp. germplasm in this collection has not been fully analyzed. In this study, we employed a single-dose SNP panel to genotype 901 sugarcane accessions, representing six Saccharum species and various hybrids. Our analysis uncovered a high rate of clone mislabeling in the collection. Specifically, we identified 86 groups of duplicates, characterized by identical SNP genotypes, which encompassed 211 accessions (23% of the total clones), while 135 groups, constituting 471 clones (52% of the total), exhibited near-identical genotypes. In addition, twenty-seven homonymous groups were detected, which shared the same clone name but differed in SNP genotypes. Hierarchical analysis of population structure partitioned the Saccharum germplasm into five clusters, corresponding to S. barberi, S. sinense, S. officinarum, S. spontaneum and S. robustum/S. edule. An assignment test, based on the five Saccharum species, enabled correcting 141 instances of mislabeled species memberships and inaccuracies. Moreover, we clarified the species membership and parentage of 298 clones that had ambiguous passport records (e.g., 'Saccharum spp', 'unknown', and 'hybrid'). Population structure and genetic diversity in these five species were further supported by Principal Coordinate Analysis and neighbor-joining clustering analysis. Analysis of Molecular Variance revealed that within-species genetic variations accounted for 85% of the total molecular variance, with the remaining 15% attributed to among-species genetic variations. The single-dose SNP markers developed in this study offer a robust tool for characterizing sugarcane germplasm worldwide. These findings have important implications for sugarcane genebank management, germplasm exchange, and crop genetic improvement.

Animal board invited review: Practical applications of genomic information in livestock

Access to high-dimensional genomic information in many livestock species is accelerating. This has been greatly aided not only by continual reductions in genotyping costs but also an expansion in the services available that leverage genomic information to create a greater return-on-investment. Genomic information on individual animals has many uses including (1) parentage verification and discovery, (2) traceability, (3) karyotyping, (4) sex determination, (5) reporting and monitoring of mutations conferring major effects or congenital defects, (6) better estimating inbreeding of individuals and coancestry among individuals, (7) mating advice, (8) determining breed composition, (9) enabling precision management, and (10) genomic evaluations; genomic evaluations exploit genome-wide genotype information to improve the accuracy of predicting an animal's (and by extension its progeny's) genetic merit. Genomic data also provide a huge resource for research, albeit the outcome from this research, if successful, should eventually be realised through one of the ten applications already mentioned. The process for generating a genotype all the way from sample procurement to identifying erroneous genotypes is described, as are the steps that should be considered when developing a bespoke genotyping panel for practical application.

An allele-sharing, moment-based estimator of global, population-specific and population-pair FST under a general model of population structure

Being able to properly quantify genetic differentiation is key to understanding the evolutionary potential of a species. One central parameter in this context is FST, the mean coancestry within populations relative to the mean coancestry between populations. Researchers have been estimating FST globally or between pairs of populations for a long time. More recently, it has been proposed to estimate population-specific FST values, and population-pair mean relative coancestry. Here, we review the several definitions and estimation methods of FST, and stress that they provide values relative to a reference population. We show the good statistical properties of an allele-sharing, method of moments based estimator of FST (global, population-specific and population-pair) under a very general model of population structure. We point to the limitation of existing likelihood and Bayesian estimators when the populations are not independent. Last, we show that recent attempts to estimate absolute, rather than relative, mean coancestry fail to do so.

Patterns of genetic differentiation imply distinct phylogeographic history of the mosquito species Anopheles messeae and Anopheles daciae in Eurasia

Detailed knowledge of phylogeography is important for control of mosquito species involved in the transmission of human infectious diseases. Anopheles messeae is a geographically widespread and genetically diverse dominant vector of malaria in Eurasia. A closely related species, An. daciae, was originally distinguished from An. messeae based on five nucleotide substitutions in its ribosomal DNA (rDNA). However, the patterns of phylogeographic history of these species in Eurasia remain poorly understood. Here, using internal transcribed spacer 2 (ITS2) of rDNA and karyotyping for the species identification we determined the composition of five Anopheles species in 28 locations in Eurasia. Based on the frequencies of 11 polymorphic chromosomal inversions used as genetic markers, a large-scale population genetics analysis was performed of 1932 mosquitoes identified as An. messeae, An. daciae and their hybrids. The largest genetic differences between the species were detected in the X sex chromosome suggesting a potential involvement of this chromosome in speciation. The frequencies of autosomal inversions in the same locations differed by 13%-45% between the species demonstrating a restricted gene flow between the species. Overall, An. messeae was identified as a diverse species with a more complex population structure than An. daciae. The clinal gradients in frequencies of chromosomal inversions were determined in both species implicating their possible involvement in climate adaptations. The frequencies of hybrids were low ~1% in northern Europe but high up to 50% in south-eastern populations. Thus, our study revealed critical differences in patterns of phylogeographic history between An. messeae and An. daciae in Eurasia. This knowledge will help to predict the potential of the malaria transmission in the northern territories of the continent.

Genetic connectivity and population expansion inferred from multilocus analysis in Lutjanus alexandrei (Lutjanidae-Perciformes), an endemic snapper from Northeastern Brazilian coast

Previous studies about the genetic diversity, connectivity and demographic history in Lutjanidae fishes have reported a common pattern of genetic homogeneity and expansion in populations from Western South Atlantic. In the present work, we inferred the population structure, the levels of genetic diversity and the demographic history of the Brazilian snapper Lutjanus alexandrei, a recently described and endemic species from Northeastern coast of Brazil. Five different fragments, including mitochondrial DNA (Control Region, Cyt b and ND4) and nuclear DNA (Myostatin and S7) regions were analyzed in 120 specimens of L. alexandrei from four localities in Northeastern Brazil, representing the first study of population genetics in this species. High levels of genetic diversity were observed following a panmictic pattern, probably related to the larval dispersal by the current tides along the Brazilian coast. In addition, both demographic history and neutrality tests indicated that L. alexandrei has undergone population expansion during Pleistocene. In this sense, the sea level variation from this period could have increased the available resources and suitable habitats for the Brazilian snapper.

Multi-locus genomic signatures of local adaptation to snow across the landscape in California populations of a willow leaf beetle

Organisms living in mountains contend with extreme climatic conditions, including short growing seasons and long winters with extensive snow cover. Anthropogenic climate change is driving unprecedented, rapid warming of montane regions across the globe, resulting in reduced winter snowpack. Loss of snow as a thermal buffer may have serious consequences for animals overwintering in soil, yet little is known about how variability in snowpack acts as a selective agent in montane ecosystems. Here, we examine genomic variation in California populations of the leaf beetle Chrysomela aeneicollis, an emerging natural model system for understanding how organisms respond to climate change. We used a genotype-environment association approach to identify genomic signatures of local adaptation to microclimate in populations from three montane regions with variable snowpack and a coastal region with no snow. We found that both winter-associated environmental variation and geographical distance contribute to overall genomic variation across the landscape. We identified non-synonymous variation in novel candidate loci associated with cytoskeletal function, ion transport and membrane stability, cellular processes associated with cold tolerance in other insects. These findings provide intriguing evidence that variation in snowpack imposes selective gradients in montane ecosystems.

Does subnetting and port hardening influence human adversarial decisions? An investigation via a HackIT tool

Prior research in cyber deception has investigated the effectiveness of the timing of deception on human decisions using simulation tools. However, there exists a gap in the literature on how the availability of subnets and port-hardening influence human decisions to attack a system. We tested the influence of subnets and port-hardening on human attack decisions in a simulated environment using the HackIT tool. Availability of subnets (present/absent) within a network and port-hardening (easy-to-attack/hard-to-attack) were varied across four between-subject conditions (N = 30 in each condition): with-subnet with easy-to-attack, with-subnet with hard-to-attack, without-subnet with easy-to-attack, and without-subnet with hard-to-attack. In with-subnet conditions, 40 systems were connected in a hybrid topology network with ten subnets connected linearly, and each subnet contained four connected systems. In without-subnet conditions, all 40 systems were connected in a bus topology. In hard-to-attack (easy-to-attack) conditions, the probabilities of successfully attacking real systems and honeypots were kept low (high) and high (low), respectively. In an experiment, human participants were randomly assigned to one of the four conditions to attack as many real systems as possible and steal credit card information. Results revealed a significant decrease in the proportion of real system attacks in the availability of subnetting and port hardening within the network. Also, more honeypots were attacked in with-subnet conditions than without-subnet conditions. Moreover, a significantly lower proportion of real systems were attacked in the port-hardened condition. This research highlights the implications of subnetting and port-hardening with honeypots to reduce real system attacks. These findings are relevant in developing advanced intrusion detection systems trained on hackers' behavior.

Fast and accurate joint inference of coancestry parameters for populations and/or individuals

We introduce a fast, new algorithm for inferring from allele count data the FST parameters describing genetic distances among a set of populations and/or unrelated diploid individuals, and a tree with branch lengths corresponding to FST values. The tree can reflect historical processes of splitting and divergence, but seeks to represent the actual genetic variance as accurately as possible with a tree structure. We generalise two major approaches to defining FST, via correlations and mismatch probabilities of sampled allele pairs, which measure shared and non-shared components of genetic variance. A diploid individual can be treated as a population of two gametes, which allows inference of coancestry coefficients for individuals as well as for populations, or a combination of the two. A simulation study illustrates that our fast method-of-moments estimation of FST values, simultaneously for multiple populations/individuals, gains statistical efficiency over pairwise approaches when the population structure is close to tree-like. We apply our approach to genome-wide genotypes from the 26 worldwide human populations of the 1000 Genomes Project. We first analyse at the population level, then a subset of individuals and in a final analysis we pool individuals from the more homogeneous populations. This flexible analysis approach gives advantages over traditional approaches to population structure/coancestry, including visual and quantitative assessments of long-standing questions about the relative magnitudes of within- and between-population genetic differences.

Genetic Diversity, Population Structure and Subset Development in a Sesbania sesban Collection

Sesbania sesban (L.) Merr. is a multipurpose legume tree grown primarily for fodder and forage in the tropical and subtropical world. In this study, the Sesbania sesban collection maintained in the International Livestock Research Institute (ILRI) forage Genebank was studied using genome-wide markers generated on the DArTseq platform. Genotyping produced 84,673 and 60,626 SNP and SilicoDArT markers with a mean polymorphic information content of 0.153 and 0.123, respectively. From the generated markers, 7587 and 15,031 highly informative SNP and SilicoDArT markers, respectively, were filtered and used for genetic diversity analysis and subset development. Analysis of molecular variance (AMOVA) revealed higher variability 'within' (52.73% for SNP markers and 67.36% for SilicoDArT markers) than 'between' accessions. Hierarchical cluster analysis showed the presence of four main clusters in the collection. Mantel correlation analysis showed a lack of relationship between genetic variation of the germplasm and their geographical origin. A representative subset of 34 accessions containing germplasm from diverse origins and agro-ecologies was developed using SNP markers. The genetic diversity information generated in this study could be used for marker-assisted screening for stress tolerance, gap analysis and identification and acquisition of new distinct genotype(s) to broaden the genetic basis of the collection for future improvement programs to develop high-yielding, stress-tolerant varieties for enhancing food and environmental security in crop-livestock-based production systems.

Russian honey bee genotype identification through enhanced marker panel set

Russian honey bees (RHB) are a breeding population developed by USDA-ARS as an effort to provide Varroa-resistant honey bees to beekeepers. The selection strategy for this breeding population was the first in honey bees to incorporate genetic stock identification (GSI). The original GSI approach has been in use for over a decade, and though effective, novel technologies and analytical approaches recently developed provide an opportunity for improvement. Here we outline a novel genotyping assay that capitalizes on the markers used in the GSI as well as new loci recently identified in a whole genome pooled study of commercial honey bee stocks. Our approach utilizes a microfluidic platform and machine learning analyses to arrive at an accurate, high throughput assay. This novel approach provides an improved tool that can be readily incorporated into breeding decisions towards healthier more productive bees.

Genetic profile in genes associated with muscle injuries and injury etiology in professional soccer players

Many causes define injuries in professional soccer players. In recent years, the study of genetics in association with injuries has been of great interest. The purpose of this study was to examine the relationship between muscle injury-related genes, injury risk and injury etiology in professional soccer players. In a cross-sectional cohort study, one hundred and twenty-two male professional football players were recruited. AMPD1 (rs17602729), ACE (rs4646994), ACTN3 (rs1815739), CKM (rs8111989) and MLCK (rs2849757 and rs2700352) polymorphisms were genotyped by using Single Nucleotide Primer Extension (SNPE). The combined influence of the six polymorphisms studied was calculated using a total genotype score (TGS). A genotype score (GS) of 2 was assigned to the “protective” genotype for injuries, a GS of 1 was assigned to the heterozygous genotype while a GS of 0 was assigned to the “worst” genotype. Injury characteristics and etiology during the 2021/2022 season were classified following a Consensus Statement for injuries recording. The distribution of allelic frequencies in the AMPD1 and MLCK c.37885C&gt;A polymorphisms were different between non-injured and injured soccer players (p &lt; 0.001 and p = 0.003, respectively). The mean total genotype score (TGS) in non-injured soccer players (57.18 ± 14.43 arbitrary units [a.u.]) was different from that of injured soccer players (51.71 ± 12.82 a.u., p = 0.034). There was a TGS cut-off point (45.83 a.u.) to discriminate non-injured from injured soccer players. Players with a TGS beyond this cut-off had an odds ratio of 1.91 (95%CI: 1.14–2.91; p = 0.022) to suffer an injury when compared with players with lower TGS. In conclusion, TGS analysis in muscle injury-related genes presented a relationship with professional soccer players at increased risk of injury. Future studies will help to develop this TGS as a potential tool to predict injury risk and perform prevention methodology in this cohort of football players.

Genetic adaptation of skin pigmentation in highland Tibetans

Strong ultraviolet (UV) radiation at high altitude imposes a serious selective pressure, which may induce skin pigmentation adaptation of indigenous populations. We conducted skin pigmentation phenotyping and genome-wide analysis of Tibetans in order to understand the underlying mechanism of adaptation to UV radiation. We observe that Tibetans have darker baseline skin color compared with lowland Han Chinese, as well as an improved tanning ability, suggesting a two-level adaptation to boost their melanin production. A genome-wide search for the responsible genes identifies GNPAT showing strong signals of positive selection in Tibetans. An enhancer mutation (rs75356281) located in GNPAT intron 2 is enriched in Tibetans (58%) but rare in other world populations (0 to 18%). The adaptive allele of rs75356281 is associated with darker skin in Tibetans and, under UVB treatment, it displays higher enhancer activities compared with the wild-type allele in in vitro luciferase assays. Transcriptome analyses of gene-edited cells clearly show that with UVB treatment, the adaptive variant of GNPAT promotes melanin synthesis, likely through the interactions of CAT and ACAA1 in peroxisomes with other pigmentation genes, and they act synergistically, leading to an improved tanning ability in Tibetans for UV protection.

Plant‐associate interactions and diversification across trophic levels

Interactions between species are widely understood to have promoted the diversification of life on Earth, but how interactions spur the formation of new species remains unclear. Interacting species often become locally adapted to each other, but they may also be subject to shared dispersal limitations and environmental conditions. Moreover, theory predicts that different kinds of interactions have different effects on diversification. To better understand how species interactions promote diversification, we compiled population genetic studies of host plants and intimately associated herbivores, parasites, and mutualists. We used Bayesian multiple regressions and the BEDASSLE modeling framework to test whether host and associate population structures were correlated over and above the potentially confounding effects of geography and shared environmental variation. We found that associates' population structure often paralleled their hosts' population structure, and that this effect is robust to accounting for geographic distance and climate. Associate genetic structure was significantly explained by plant genetic structure somewhat more often in antagonistic interactions than in mutualistic ones. This aligns with a key prediction of coevolutionary theory that antagonistic interactions promote diversity through local adaptation of antagonists to hosts, while mutualistic interactions more often promote diversity via the effect of hosts' geographic distribution on mutualists' dispersal.

Whole-genome resequencing reveals the origin of tea in Lincang

Phylogeographic, population genetics and diversity analysis are crucial for local tea resource conservation and breeding programs. Lincang in Yunnan has been known as the possible place of domestication for tea worldwide, yet, its genetic makeup and unique Lincang origin are little understood. Here, we reported a large-scale whole-genome resequencing based population genomic analysis in eight main tea-producing areas of Lincang in Yunnan (1,350 accessions), and the first comprehensive map of tea genome variation in Lincang was constructed. Based on the population structure, tea sample in Lincang was divided into three subgroups, and inferred Xigui and Nahan Tea Mountain in Linxiang, Baiying Mountain Ancient Tea Garden in Yun, and Jinxiu Village of Xiaowan Town in Fengqing, which belong to the birthplace of the three subgroups, were all likely to be the origin center of Lincang tea. Meanwhile, the history population sizes analysis show that similar evolutionary patterns were observed for the three subgroups of Lincang. It also was observed that the hybrid among eight areas of Lincang was noticeable, resulting in insignificant genetic differentiation between geographical populations and low genetic diversity. The findings of this study clarified the genetic make-up and evolutionary traits of the local population of tea, which gave some insight into the development of Lincang tea.

Evolutionary Genomics of a Subdivided Species

The ways in which genetic variation is distributed within and among populations is a key determinant of the evolutionary features of a species. However, most comprehensive studies of these features have been restricted to studies of subdivision in settings known to have been driven by local adaptation, leaving our understanding of the natural dispersion of allelic variation less than ideal. Here, we present a geographic population-genomic analysis of 10 populations of the freshwater microcrustacean Daphnia pulex, an emerging model system in evolutionary genomics. These populations exhibit a pattern of moderate isolation-by-distance, with an average migration rate of 0.6 individuals per generation, and average effective population sizes of ∼650,000 individuals. Most populations contain numerous private alleles, and genomic scans highlight the presence of islands of excessively high population subdivision for more common alleles. A large fraction of such islands of population divergence likely reflect historical neutral changes, including rare stochastic migration and hybridization events. The data do point to local adaptive divergence, although the precise nature of the relevant variation is diffuse and cannot be associated with particular loci, despite the very large sample sizes involved in this study. In contrast, an analysis of between-species divergence highlights positive selection operating on a large set of genes with functions nearly nonoverlapping with those involved in local adaptation, in particular ribosome structure, mitochondrial bioenergetics, light reception and response, detoxification, and gene regulation. These results set the stage for using D. pulex as a model for understanding the relationship between molecular and cellular evolution in the context of natural environments.

A gene-based capture assay for surveying patterns of genetic diversity and insecticide resistance in a worldwide group of invasive mosquitoes

Understanding patterns of diversification, genetic exchange, and pesticide resistance in arthropod disease vectors is necessary for effective population management. With the availability of next-generation sequencing technologies, one of the best approaches for surveying such patterns involves the simultaneous genotyping of many samples for a large number of genetic markers. To this end, the targeting of gene sequences of known function can be a cost-effective strategy. One insect group of substantial health concern are the mosquito taxa that make up the Culex pipiens complex. Members of this complex transmit damaging arboviruses and filariae worms to humans, as well as other pathogens such as avian malaria parasites that are detrimental to birds. Here we describe the development of a targeted, gene-based assay for surveying genetic diversity and population structure in this mosquito complex. To test the utility of this assay, we sequenced samples from several members of the complex, as well as from distinct populations of the relatively under-studied Culex quinquefasciatus. The data generated was then used to examine taxonomic divergence and population clustering between and within these mosquitoes. We also used this data to investigate genetic variants present in our samples that had previously been shown to correlate with insecticide-resistance. Broadly, our gene capture approach successfully enriched the genomic regions of interest, and proved effective for facilitating examinations of taxonomic divergence and geographic clustering within the Cx. pipiens complex. It also allowed us to successfully survey genetic variation associated with insecticide resistance in Culex mosquitoes. This enrichment protocol will be useful for future studies that aim to understand the genetic mechanisms underlying the evolution of these ubiquitous and increasingly damaging disease vectors.

The mosquito taxa that make up the Culex pipiens complex are important vectors of the agents of several human diseases such as West Nile and St. Louis encephalitides, and lymphatic filariasis. They are also important vectors of avian malaria, which impacts livestock and wildlife. The development of effective strategies for the control of these mosquitoes requires knowledge of their origins, distribution, dispersal patterns, and the extent to which discreet taxonomic entities within the complex interbreed. To achieve these objectives, it is necessary to compare patterns of genetic diversity across many mosquito samples, which can be cost-prohibitive. To address this limitation, we developed a targeted, gene-based assay that allowed us to cost-effectively genotype a large number of genetic variants from a representative global sampling of individual Cx. pipiens complex mosquitoes. We show that this assay is a powerful tool for examining genetic structure and hybridization among populations. We also explore its utility for surveying alleles previously shown to be associated with insecticide resistance. Future use of this enrichment assay and the bioinformatics methods described here will allow researchers to study evolutionary patterns across the Cx. pipiens complex as well as monitor the presence of genetic variation that could affect control efforts.

Population Structure of German Cockroaches (Blattodea: Ectobiidae) in an Urban Environment Based on Single Nucleotide Polymorphisms

German cockroaches (Blattella germanica L.) harbor and disperse medically important pathogens and are a source of allergens that impact human health and wellbeing. Management of this pest requires an understanding of their distribution and dispersal. In this study, we collected German cockroaches from three apartment buildings in New Jersey, USA. We identified single-nucleotide polymorphisms (SNPs) from DNA extractions using next generation sequencing. We analyzed the SNPs and characterized cockroach population genetic structure using Fst, principal component, phylogenetic, and STRUCTURE analyses. We found significant differences in German cockroach population structure among the buildings. Within buildings, we found variable population structure that may be evidence for multiple colonization events. This study shows that SNPs derived from next generation sequencing provide a powerful tool for analyzing the genetic population structure of these medically important pests.

Maltese Allelic Variants in Corneal Dystrophy Genes in a Worldwide Setting

INTRODUCTION

This study aimed to establish which worldwide population cohorts have a genetic make-up closest to that of a large sample of the Maltese population with regard to corneal dystrophy (CD) genes.

METHODS

Single nucleotide polymorphisms (SNPs) in the Maltese cohort were compared with worldwide cohorts. Fixation index (F_ST) values were calculated to evaluate population differentiation. The genetic prevalence of CD subtypes in worldwide and Maltese cohorts were calculated, and single nucleotide missense mutations present in the Maltese cohort were evaluated for potential pathogenicity.

RESULTS

F_ST values showed that CD-related genes differ substantially among the studied cohorts. F_ST values for each SNP showed greatest differentiation between the Maltese and African cohorts and least differentiation with the Puerto Rican, Mexican, and Colombian cohorts. One TGFBI casual CD mutation, 502V, which causes a Bowman's layer CD/atypical Thiel-Behnke CD was identified in the Maltese cohort. The KRT3 NC_000012.11:g.53186088G>C mutation was potentially deleterious.

CONCLUSION

Identifying populations with least genetic differentiation can facilitate and help guide future diagnostic and treatment strategies for Maltese individuals with CDs in the absence of comparable Maltese data. Analysing the previously unknown CD genetic pool present in a large Maltese cohort adds to the global genetic bank that researchers rely on for medical progress.

Genetic Diversity and Population Structure of Psidium Species from Restinga: A Coastal and Disturbed Ecosystem of the Brazilian Atlantic Forest

The Atlantic Forest is one of the most threatened biomes in the world. Here, we use a common set of microsatellite markers to assess the genetic diversity and population structure of three species from the genus Psidium (P. guajava, P. macahense, and P. guineense), located in a disturbed environment of the Atlantic Forest, the restinga, in Espírito Santo, Brazil. Psidium guajava populations presented the highest number of alleles (95) followed by P. guineense (81) and P. macahense (68). The genetic variability was high (P. guajava = 0.71; P. guineense = 0.74; P. macahense = 0.63), with greater variation within populations (72 to 84%) than among populations (15 to 27%), reflecting elevated values of genetic differentiation (P. guajava, F_ST: 0.15; P. macahense, Ø_ST: 0.27; P. guineense, Ø_ST: 0.21). The populations were clustered into two main groups and considered moderately structured. This is the first report of genetic studies and evidence of polyploidy to P. macahense. Our results may provide information that can be used in management and conservation strategies, to preserve the diversity of Psidium populations.

EDClust: an EM-MM hybrid method for cell clustering in multiple-subject single-cell RNA sequencing

MOTIVATION

Single-cell RNA sequencing (scRNA-seq) has revolutionized biological research by enabling the measurement of transcriptomic profiles at the single-cell level. With the increasing application of scRNA-seq in larger-scale studies, the problem of appropriately clustering cells emerges when the scRNA-seq data are from multiple subjects. One challenge is the subject-specific variation; systematic heterogeneity from multiple subjects may have a significant impact on clustering accuracy. Existing methods seeking to address such effects suffer from several limitations.

RESULTS

We develop a novel statistical method, EDClust, for multi-subject scRNA-seq cell clustering. EDClust models the sequence read counts by a mixture of Dirichlet-multinomial distributions and explicitly accounts for cell-type heterogeneity, subject heterogeneity and clustering uncertainty. An EM-MM hybrid algorithm is derived for maximizing the data likelihood and clustering the cells. We perform a series of simulation studies to evaluate the proposed method and demonstrate the outstanding performance of EDClust. Comprehensive benchmarking on four real scRNA-seq datasets with various tissue types and species demonstrates the substantial accuracy improvement of EDClust compared to existing methods.

AVAILABILITY AND IMPLEMENTATION

The R package is freely available at https://github.com/weix21/EDClust.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Whole genome sequencing of Luxi Black Head sheep for screening selection signatures associated with important traits

Objective

Luxi Black Head sheep (LBH) is the first crossbreed specialized for meat production and was developed by crossbreeding Black Head Dorper sheep (DP) and Small Tailed Han sheep (STH) in the farming areas of northern China. Research on the genomic variations and selection signatures of LBH caused by continuous artificial selection is of great significance for identifying the genetic mechanisms of important traits of sheep and for the continuous breeding of LBH.

Methods

We explored the genetic relationships of LBH, DP, and several Mongolian sheep breeds by constructing phylogenetic tree, principal component analysis and linkage disequilibrium analysis. In addition, we analysed 29 whole genomes of sheep. The genome-wide selection signatures have been scanned with four methods: heterozygosity (H_P), fixation index (F_ST), cross-population extended haplotype homozygosity (XP-EHH) and the nucleotide diversity (θ_π) ratio.

Results

The genetic relationships analysis showed that LBH appeared to be an independent cluster closer to DP. The candidate signatures of positive selection in sheep genome revealed candidate genes for developmental process (HoxA gene cluster, BCL2L11, TSHR), immunity (CXCL6, CXCL1, SKAP2, PTK6, MST1R), growth (PDGFD, FGF18, SRF, SOCS2), and reproduction (BCAS3, TRIM24, ASTL, FNDC3A). Moreover, two signalling pathways closely related to reproduction, the thyroid hormone signalling pathway and the oxytocin signalling pathway, were detected.

Conclusion

The selective sweep analysis of LBH genome revealed candidate genes and signalling pathways associated with developmental process, immunity, growth, and reproduction. Our findings provide a valuable resource for sheep breeding and insight into the mechanisms of artificial selection.

Genome-Wide Selection Sweep between Wild and Local Pigs from Europe for the Investigation of the Hereditary Characteristics of Domestication in Sus Scrofa

The phenotypic characteristics of existing domestic pigs (DPs) greatly differ from those of wild boar (WB) populations thousands of years ago. After thousands of years of human domestication, WB and DP have diverged greatly in terms of genetics. Theoretically, worldwide local pigs have independent contributions from their local WBs at the beginning of Sus scrofa domestication. The investigation of the vicissitude of the heredity material between domestic populations and their wild ancestors will help in further understanding the domestication history of domestic animals. In the present study, we performed a genome-wide association scan (GWSA) and phylogeny estimation with a total of 1098 public European Illumina 60K single nucleotide polymorphism data, which included 650 local DPs and 448 WBs. The results revealed that the phylogenetic relationship of WBs corresponds to their geographical distribution and carries large divergence with DPs, and all WB breeds (e.g., HRWB, SBWB, and TIWB) presents a closely linkage with the middle WB (e.g., HRWB, and PLWB). In addition, 64 selected candidate genes (e.g., IDH2, PIP5K1B, SMARCA2, KIF5C, and TJP2) were identified from GWSA. A total of 63 known multiple biological functional pathways were annotated by 22 genes, and ubiquinone and other terpenoid-quinone biosynthesis pathways that belong to the metabolism of cofactors and vitamins were significantly enriched (p < 0.05). The most frequent (28.57%) pathways were classified under metabolism. We confirmed that the middle European WB has made an important genetic contribution to the entire European WB populations. A series of selected genes discovered from this study provides the scientific community with a deeper understanding of the heredity performance of metabolism and emotion and the real purpose behind domestication.

Detection of QTLs for Plant Height Architecture Traits in Rice (Oryza sativa L.) by Association Mapping and the RSTEP-LRT Method

Plant height (PH) and its component traits are critical determinants of lodging resistance and strongly influence yield in rice. The genetic architecture of PH and its component traits were mined in two mapping populations. In the natural population composed of 504 accessions, a total of forty simple sequence repeat (SSR) markers associated with PH and its component traits were detected across two environments via association mapping. Allele RM305-210 bp on chromosome 5 for PH had the largest phenotypic effect value (PEV) (−51.42 cm) with a reducing effect. Allele RM3533-220 bp on chromosome 9 for panicle length and allele RM264-120 bp on chromosome 8 for the length of upper first elongated internode (1IN) showed the highest positive PEV. Among the elongated internodes with negative effects being desirable, the allele RM348-130 bp showed the largest PEV (−7.48 cm) for the length of upper second elongated internode. In the chromosome segment substitution line population consisting of 53 lines, a total of nine QTLs were detected across two environments, with the phenotypic variance explained (PVE) ranging 10.07–28.42%. Among the detected QTLs, q1IN-7 explained the largest PVE (28.42%) for the 1IN, with an additive of 5.31 cm. The favorable allele RM257-125 bp on chromosome 9 for the 1IN increasing was detected in both populations. The favorable alleles provided here could be used to shape PH architecture against lodging.

Genomic Screening to Identify Food Trees Potentially Dispersed by Precolonial Indigenous Peoples

Over millennia, Indigenous peoples have dispersed the propagules of non-crop plants through trade, seasonal migration or attending ceremonies; and potentially increased the geographic range or abundance of many food species around the world. Genomic data can be used to reconstruct these histories. However, it can be difficult to disentangle anthropogenic from non-anthropogenic dispersal in long-lived non-crop species. We developed a genomic workflow that can be used to screen out species that show patterns consistent with faunal dispersal or long-term isolation, and identify species that carry dispersal signals of putative human influence. We used genotyping-by-sequencing (DArTseq) and whole-plastid sequencing (SKIMseq) to identify nuclear and chloroplast Single Nucleotide Polymorphisms in east Australian rainforest trees (4 families, 7 genera, 15 species) with large (>30 mm) or small (<30 mm) edible fruit, either with or without a known history of use by Indigenous peoples. We employed standard population genetic analyses to test for four signals of dispersal using a limited and opportunistically acquired sample scheme. We expected different patterns for species that fall into one of three broadly described dispersal histories: (1) ongoing faunal dispersal, (2) post-megafauna isolation and (3) post-megafauna isolation followed by dispersal of putative human influence. We identified five large-fruited species that displayed strong population structure combined with signals of dispersal. We propose coalescent methods to investigate whether these genomic signals can be attributed to post-megafauna isolation and dispersal by Indigenous peoples.

Plant collections for conservation and restoration: can they be adapted and adaptable?

Plant collections are important for the conservation of threatened species, and can provide material for ecological restoration. Typically we want collections to have high genetic diversity so populations founded from it are adaptable to future challenges. Sometimes, we have additional objectives for collections, such as enrichment for desirable traits controlled by adaptive alleles. We used landscape genomic datasets for two plants, Westringia fruticosa and Wilkiea huegeliana, to design collections that are genetically diverse, and that are adapted to warming climates. We characterized temperature adaptation by: (i) using the mean annual temperature of the sites of origin of the plants, and (ii) using the representation of alleles that are associated with warm temperatures. In Westringia fruticosa, there was a negative correlation, or tradeoff, between designing a collection that was both genetically diverse and adapted to warm temperatures. This tradeoff was weaker in Wilkiea huegeliana. We hypothesized this was because neutral genetic variation was strongly correlated with temperature in Westringia fruticosa, and not in Wilkiea huegeliana. Accordingly, when we shuffled the temperature data, breaking up the covariance between Westringia fruticosa genetic variation and temperature, there was a relaxation of the observed tradeoff. In sum, we explore tradeoffs between promoting genetic diversity and selecting for a specific trait in plant collections, and show that the strength of this tradeoff varies between two species. This represents a useful step towards understanding when selection will have a large cost in genetic diversity, and when it will be possible to design a collection that is both adapted and adaptable.

Genetic isolation by distance underlies colour pattern divergence in red-eyed treefrogs (Agalychnis callidryas)

Investigating the spatial distribution of genetic and phenotypic variation can provide insights into the evolutionary processes that shape diversity in natural systems. We characterized patterns of genetic and phenotypic diversity to learn about drivers of colour-pattern diversification in red-eyed treefrogs (Agalychnis callidryas) in Costa Rica. Along the Pacific coast, red-eyed treefrogs have conspicuous leg colour patterning that transitions from orange in the north to purple in the south. We measured phenotypic variation of frogs, with increased sampling at sites where the orange-to-purple transition occurs. At the transition zone, we discovered the co-occurrence of multiple colour-pattern morphs. To explore possible causes of this variation, we generated a single nucleotide polymorphism data set to analyse population genetic structure, measure genetic diversity and infer the processes that mediate genotype-phenotype dynamics. We investigated how patterns of genetic relatedness correspond to individual measures of colour pattern along the coast, including testing for the role of hybridization in geographic regions where orange and purple phenotypic groups co-occur. We found no evidence that colour-pattern polymorphism in the transition zone arose through recent hybridization. Instead, a strong pattern of genetic isolation by distance indicates that colour-pattern variation was either retained through other processes such as ancestral colour polymorphisms or ancient secondary contact, or else it was generated by novel mutations. We found that phenotype changes along the Pacific coast more than would be expected based on genetic divergence and geographic distance alone. Combined, our results suggest the possibility of selective pressures acting on colour pattern at a small geographic scale.

Haplotype-resolved de novo assembly of a Tujia genome suggests the necessity for high-quality population-specific genome references

Even though the human reference genome assembly is continually being improved, it remains debatable whether a population-specific reference is necessary for every ethnic group. Here, we de novo assembled an individual genome (TJ1) from the Tujia population, an ethnic minority group most closely related to the Han Chinese. TJ1 provided a high-quality haplotype-resolved assembly of chromosome-scale with a scaffold N50 size >78 Mb. Compared with GRCh38 and other de novo assemblies, TJ1 improved short-read mapping, enhanced calling precision for structural variants, and detected rare and low-frequency variants. This revealed fine-scale differences between the closely related Han and Tujia populations, such as population-stratified variants of LCT and UBXN8, and improved screening for ancestry informative markers. We demonstrated that TJ1 could reduce false positives in clinical diagnosis and analyzed the PRSS1-PRSS2 locus as a test case. Our results suggest that population-specific assemblies are necessary for genetic and medical analysis, especially when closely related populations are studied. A record of this paper's transparent peer review process is included in the supplemental information.

Worldwide diversity, association potential, and natural selection in the superimposed taste genes, CD36 and GNAT3

CD36 and GNAT3 mediate taste responses, with CD36 acting as a lipid detector and GNAT3 acting as the α subunit of gustducin, a G protein governing sweet, savory, and bitter transduction. Strikingly, the genes encoding CD36 and GNAT3 are genomically superimposed, with CD36 completely encompassing GNAT3. To characterize genetic variation across the CD36-GNAT3 region, its implications for phenotypic diversity, and its recent evolution, we analyzed from ~2,500 worldwide subjects sequenced by the 1000 Genomes Project (1000GP). CD36-GNAT3 harbored extensive diversity including 8,688 single-nucleotide polymorphisms (SNPs), 414 indels, and other complex variants. Sliding window analyses revealed that nucleotide diversity and population differentiation across CD36-GNAT3 were consistent with genome-wide trends in the 1000GP (π = 0.10%, P = 0.64; FST = 9.0%, P = 0.57). In addition, functional predictions using SIFT and PolyPhen-2 identified 60 variants likely to alter protein function, and they were in weak linkage disequilibrium (r2 < 0.17), suggesting their effects are largely independent. However, the frequencies of predicted functional variants were low (P¯ = 0.0013), indicating their contributions to phenotypic variance on population scales are limited. Tests using Tajima's D statistic revealed that pressures from natural selection have been relaxed across most of CD36-GNAT3 during its recent history (0.39 < P < 0.67). However, CD36 exons showed signs of local adaptation consistent with prior reports (P < 0.035). Thus, CD36 and GNAT3 harbor numerous variants predicted to affect taste sensitivity, but most are rare and phenotypic variance on a population level is likely mediated by a small number of sites.

Rank-invariant estimation of inbreeding coefficients

The two alleles an individual carries at a locus are identical by descent (ibd) if they have descended from a single ancestral allele in a reference population, and the probability of such identity is the inbreeding coefficient of the individual. Inbreeding coefficients can be predicted from pedigrees with founders constituting the reference population, but estimation from genetic data is not possible without data from the reference population. Most inbreeding estimators that make explicit use of sample allele frequencies as estimates of allele probabilities in the reference population are confounded by average kinships with other individuals. This means that the ranking of those estimates depends on the scope of the study sample and we show the variation in rankings for common estimators applied to different subdivisions of 1000 Genomes data. Allele-sharing estimators of within-population inbreeding relative to average kinship in a study sample, however, do have invariant rankings across all studies including those individuals. They are unbiased with a large number of SNPs. We discuss how allele sharing estimates are the relevant quantities for a range of empirical applications.

Evolution of flowering time in a selfing annual plant: Roles of adaptation and genetic drift

Resurrection studies are a useful tool to measure how phenotypic traits have changed in populations through time. If these trait modifications correlate with the environmental changes that occurred during the time period, it suggests that the phenotypic changes could be a response to selection. Selfing, through its reduction of effective size, could challenge the ability of a population to adapt to environmental changes. Here, we used a resurrection study to test for adaptation in a selfing population of Medicago truncatula, by comparing the genetic composition and flowering times across 22 generations. We found evidence for evolution toward earlier flowering times by about two days and a peculiar genetic structure, typical of highly selfing populations, where some multilocus genotypes (MLGs) are persistent through time. We used the change in frequency of the MLGs through time as a multilocus fitness measure and built a selection gradient that suggests evolution toward earlier flowering times. Yet, a simulation model revealed that the observed change in flowering time could be explained by drift alone, provided the effective size of the population is small enough (<150). These analyses suffer from the difficulty to estimate the effective size in a highly selfing population, where effective recombination is severely reduced.

The correlation of substitution effects across populations and generations in the presence of nonadditive functional gene action

Allele substitution effects at quantitative trait loci (QTL) are part of the basis of quantitative genetics theory and applications such as association analysis and genomic prediction. In the presence of nonadditive functional gene action, substitution effects are not constant across populations. We develop an original approach to model the difference in substitution effects across populations as a first order Taylor series expansion from a "focal" population. This expansion involves the difference in allele frequencies and second-order statistical effects (additive by additive and dominance). The change in allele frequencies is a function of relationships (or genetic distances) across populations. As a result, it is possible to estimate the correlation of substitution effects across two populations using three elements: magnitudes of additive, dominance, and additive by additive variances; relationships (Nei's minimum distances or Fst indexes); and assumed heterozygosities. Similarly, the theory applies as well to distinct generations in a population, in which case the distance across generations is a function of increase of inbreeding. Simulation results confirmed our derivations. Slight biases were observed, depending on the nonadditive mechanism and the reference allele. Our derivations are useful to understand and forecast the possibility of prediction across populations and the similarity of GWAS effects.

Waste not, want not: Microsatellites remain an economical and informative technology for conservation genetics

Comparisons of microsatellites and single-nucleotide polymorphisms (SNPs) have found that SNPs outperform microsatellites in population genetic analyses, questioning the continued utility of microsatellites in population and landscape genetics. Yet, highly polymorphic markers may be of value in species that have reduced genetic variation. This study repeated previous analyses that used microsatellites with SNPs developed from ddRAD sequencing in the black-capped vireo source-sink system. SNPs provided greater resolution of genetic diversity, population differentiation, and migrant detection but could not reconstruct parentage relationships due to insufficient heterozygosities. The biological inferences made by both sets of markers were similar: asymmetrical gene flow from source sites to the remaining sink sites. With the landscape genetic analyses, we found different results between the two molecular markers, but associations of the top environmental features (riparian, open habitat, agriculture, and human development) with dispersal estimates were shared between marker types. Despite the higher precision of SNPs, we find that microsatellites effectively uncover population processes and patterns and are superior for parentage analyses in this species with reduced genetic diversity. This study illustrates the continued applicability and relevance of microsatellites in population genetic research.

New Cytogenetic Photomap and Molecular Diagnostics for the Cryptic Species of the Malaria Mosquitoes Anopheles messeae and Anopheles daciae from Eurasia

The Eurasian malaria vector Anopheles messeae is a widely spread and genetically diverse species. Five widespread polymorphic chromosomal inversions were found in natural populations of this mosquito. A cryptic species, Anopheles daciae, was differentiated from An. messeae by the presence of several nucleotide substitutions in the Internal Transcribed Spacer 2 (ITS2) region of ribosomal DNA. However, because of the absence of a high-quality reference cytogenetic map, the inversion polymorphisms in An. daciae and An. messeae remain poorly understood. Moreover, a recently determined heterogeneity in ITS2 in An. daciae questioned the accuracy of the previously used Restriction Fragment Length Polymorphism (RFLP) assay for species diagnostics. In this study, a standard-universal cytogenetic map was constructed based on orcein stained images of chromosomes from salivary glands for population studies of the chromosomal inversions that can be used for both An. messeae and An. daciae. In addition, a new ITS2-RFLP approach for species diagnostics was developed. Both methods were applied to characterize inversion polymorphism in populations of An. messeae and An. daciae from a single location in Western Siberia in Russia. The analysis demonstrates that cryptic species are remarkably different in their frequencies of chromosomal inversion variants. Our study supports previous observations that An. messeae has higher inversion polymorphism in all autosomes than the cryptic species An. daciae.