The enormous repetitive Antarctic krill genome reveals environmental adaptations and population insights.

Antarctic krill (Euphausia superba) is Earth's most abundant wild animal, and its enormous biomass is vital to the Southern Ocean ecosystem. Here, we report a 48.01-Gb chromosome-level Antarctic krill genome, whose large genome size appears to have resulted from inter-genic transposable element expansions. Our assembly reveals the molecular architecture of the Antarctic krill circadian clock and uncovers expanded gene families associated with molting and energy metabolism, providing insights into adaptations to the cold and highly seasonal Antarctic environment. Population-level genome re-sequencing from four geographical sites around the Antarctic continent reveals no clear population structure but highlights natural selection associated with environmental variables. An apparent drastic reduction in krill population size 10 mya and a subsequent rebound 100 thousand years ago coincides with climate change events. Our findings uncover the genomic basis of Antarctic krill adaptations to the Southern Ocean and provide valuable resources for future Antarctic research.

Genome scans for selection signatures identify candidate virulence genes for adaptation of the soybean cyst nematode to host resistance.

Plant pathogens are constantly under selection pressure for host resistance adaptation. Soybean cyst nematode (SCN, Heterodera glycines) is a major pest of soybean primarily managed through resistant cultivars; however, SCN populations have evolved virulence in response to selection pressures driven by repeated monoculture of the same genetic resistance. Resistance to SCN is mediated by multiple epistatic interactions between Rhg (for resistance to H. glycines) genes. However, the identity of SCN virulence genes that confer the ability to overcome resistance remains unknown. To identify candidate genomic regions showing signatures of selection for increased virulence, we conducted whole genome resequencing of pooled individuals (Pool-Seq) from two pairs of SCN populations adapted on soybeans with Peking-type (rhg1-a, rhg2, and Rhg4) resistance. Population differentiation and principal component analysis-based approaches identified approximately 0.72-0.79 million SNPs, the frequency of which showed potential selection signatures across multiple genomic regions. Chromosomes 3 and 6 between population pairs showed the greatest density of outlier SNPs with high population differentiation. Conducting multiple outlier detection tests to identify overlapping SNPs resulted in a total of 966 significantly differentiated SNPs, of which 285 exon SNPs were mapped to 97 genes. Of these, six genes encoded members of known stylet-secreted effector protein families potentially involved in host defence modulation including venom-allergen-like, annexin, glutathione synthetase, SPRYSEC, chitinase, and CLE effector proteins. Further functional analysis of identified candidate genes will provide new insights into the genetic mechanisms by which SCN overcomes soybean resistance and inform the development of molecular markers for rapidly screening the virulence profile of an SCN-infested field.

Genetic diversity and positive signatures of selection in indigenous cattle breeds of Iran.

Animal domestication, climate changes over time, and artificial selection have played significant roles in shaping the genome structure of various animal species, including cattle. These processes have led to the emergence of several indigenous cattle breeds with distinct genetic characteristics. This study focused on unraveling the genetic diversity and identifying candidate genomic regions in eight indigenous cattle breeds of Iran. The data consisted of ∼777 962 single nucleotide polymorphisms (SNPs) of 89 animals from Iranian indigenous cattle scattered throughout the country. We employed various methods, including integrated haplotype score, FST, and cross-population composite likelihood ratio, to conduct a genome scan for detecting selection signals within and between cattle populations. Average observed heterozygosity across the populations was 0.36, with a range of 0.32-0.40. In addition, negative and low rates of inbreeding (FIS) in the populations were observed. The genome-wide analysis revealed several genomic regions that harbored candidate genes associated with production traits (e.g., MFSD1, TYW5, ADRB2, BLK, and CRTC3), adaptation to local environmental constraints (CACNA2D1, CXCL3, and GRO1), and coat color (DYM). Finally, the study of the reported quantitative trait loci (QTL) regions in the cattle genome demonstrated that the identified regions were associated with QTL related to important traits such as milk composition, body weight, daily gain, feed conversion, and residual feed intake. Overall, this study contributes to a better understanding of the genetic diversity and potential candidate genes underlying important traits in Iranian indigenous cattle breeds, which can inform future breeding and conservation efforts.

Population genetic differentiation and phenotypic plasticity of Ambrosia artemisiifolia under different nitrogen levels.

Rapid adaptive evolution and phenotypic plasticity are two mechanisms that often underlie invasiveness of alien plant species, but whether they can co-occur within invasive plant populations under altered environmental conditions such as nitrogen (N) enrichment has seldom been explored. Latitudinal clines in plant trait responses to variation in environmental factors may provide evidence of local adaptation. Here, we inferred the relative contributions of phenotypic plasticity and local adaptation to the performance of the invasive plant Ambrosia artemisiifolia under different soil N levels, using a common garden approach. We grew A. artemisiifolia individuals raised from seeds that were sampled from six invasive populations along a wide latitudinal cline in China (23°42' N to 45°43' N) under three N (0, 5, and 10 g N m-2 ) levels in a common garden. Results show significant interpopulation genetic differentiation in plant height, number of branches, total biomass, and transpiration rate of the invader A. artemisiifolia across the N treatments. The populations also expressed genetic differentiation in basal diameter, growth rate, leaf area, seed width, root biomass, aboveground biomass, stomatal conductance, and intercellular CO2 concentration regardless of N treatments. Moreover, plants from different populations of the invader displayed plastic responses in time to first flower, hundred-grain weight, net photosynthetic rate, and relative biomass allocation to roots and shoots and seed length under different N treatments. Additionally, individuals of A. artemisiifolia from higher latitudes grew shorter and allocated less biomass to the roots regardless of N treatment, while latitudinal cline (or lack thereof) in other traits depended on the level of N in which the plants were grown. Overall, these results suggest that rapid adaptive evolution and phenotypic plasticity in the various traits that we quantified may jointly contribute to invasiveness of A. artemisiifolia under different levels of N availability. More broadly, the results support the idea that phenotypic plasticity and rapid adaptive evolution can jointly enable invasive plants to colonize a wide range of environmental conditions.

Genotyping-by-Sequencing Reveals Population Differentiation and Linkage Disequilibrium in Alternaria linariae from Tomato.

Alternaria linariae is an economically important foliar pathogen that causes early blight disease in tomatoes. Understanding genetic diversity, population genetic structure, and evolutionary potential is crucial to contemplating effective disease management strategies. We leveraged genotyping-by-sequencing (GBS) technology to compare genome-wide variation in 124 isolates of Alternaria spp. (A. alternata, A. linariae, and A. solani) for comparative genome analysis and to test the hypotheses of genetic differentiation and linkage disequilibrium (LD) in A. linariae collected from tomatoes in western North Carolina. We performed a pangenome-aware variant calling and filtering with GBSapp and identified 53,238 variants conserved across the reference genomes of three Alternaria spp. The highest marker density was observed on chromosome 1 (7 Mb). Both discriminant analysis of principal components and Bayesian model-based STRUCTURE analysis of A. linariae isolates revealed three subpopulations with minimal admixture. The genetic differentiation coefficients (FST) within A. linariae subpopulations were similar and high (0.86), indicating that alleles in the subpopulations are fixed and the genetic structure is likely due to restricted recombination. Analysis of molecular variance indicated higher variation among populations (89%) than within the population (11%). We found long-range LD between pairs of loci in A. linariae, supporting the hypothesis of low recombination expected for a fungal pathogen with limited sexual reproduction. Our findings provide evidence of a high level of population genetic differentiation in A. linariae, which reinforces the importance of developing tomato varieties with broad-spectrum resistance to various isolates of A. linariae.

A comprehensive analysis of copy number variations in diverse apple populations.

BACKGROUND: As an important source of genetic variation, copy number variation (CNV) can alter the dosage of DNA segments, which in turn may affect gene expression level and phenotype. However, our knowledge of CNV in apple is still limited. Here, we obtained high-confidence CNVs and investigated their functional impact based on genome resequencing data of two apple populations, cultivars and wild relatives. RESULTS: In this study, we identified 914,610 CNVs comprising 14,839 CNV regions (CNVRs) from 346 apple accessions, including 289 cultivars and 57 wild relatives. CNVRs summed to 71.19 Mb, accounting for 10.03% of the apple genome. Under the low linkage disequilibrium (LD) with nearby SNPs, they could also accurately reflect the population structure of apple independent of SNPs. Furthermore, A total of 3,621 genes were covered by CNVRs and functionally involved in biological processes such as defense response, reproduction and metabolic processes. In addition, the population differentiation index ([Formula: see text]) analysis between cultivars and wild relatives revealed 127 CN-differentiated genes, which may contribute to trait differences in these two populations. CONCLUSIONS: This study was based on identification of CNVs from 346 diverse apple accessions, which to our knowledge was the largest dataset for CNV analysis in apple. Our work presented the first comprehensive CNV map and provided valuable resources for understanding genomic variations in apple.

Perishing rich, expanding poor: Demography and population genetic patterns in two congeneric butterflies.

In human-altered landscapes, specialist butterflies typically form spatially restricted populations, genetically differentiated due to dispersal restrictions. Generalists, in contrast, display minimum differentiation but high genetic diversity. While local-level actions suffice to conserve specialists and landscape-level actions are necessary for generalists, minimum information exists regarding conservation of species with intermediate features. We targeted two congeneric butterflies, the recently re-expanding Argynnis adippe and the strongly declining A. niobe, co-occurring in the pastoral landscape of the Carpathian Mountains, Czech Republic. We integrated species distribution models, mark-recapture and microsatellite analysis to compare their habitat requirements, adult demography, dispersal and genetic patterns, and expanded the genetic analysis across the Carpathian Arc and beyond to delimit spatial conservation units. In two mountain valleys, both species formed interconnected populations numbering thousands of individuals. Mobility patterns suggested the populations' interconnection across the Czech Carpathians. Genetic diversity was extremely poor in the nonthreatened A. adippe and moderate in the declining A. niobe. No population differentiation was detected within the Czech Carpathians (~1500 km2 ). Low genetic diversity and no differentiation was preserved in A. adippe across East Central Europe, whereas in A. niobe, populations from Serbia were differentiated from the Carpathian Arc + Alps. The high adult mobility linked to low differentiation probably reflects the distribution of larval resources, historically widespread but sparse and currently declining for A. niobe (grazing-disturbed grounds), while currently increasing for A. adippe (abandonment scrub, disturbed woodlands). Units as large as entire mountain systems define population boundaries, and hence conservation management units, for both species.

Cultural evolution of conformity and anticonformity.

Conformist bias occurs when the probability of adopting a more common cultural variant in a population exceeds its frequency, and anticonformist bias occurs when the reverse is true. Conformist and anticonformist bias have been widely documented in humans, and conformist bias has also been observed in many nonhuman animals. Boyd and Richerson used models of conformist and anticonformist bias to explain the evolution of large-scale cooperation, and subsequent research has extended these models. We revisit Boyd and Richerson's original analysis and show that, with conformity based on more than three role models, the evolutionary dynamics can be more complex than previously assumed. For example, we show the presence of stable cycles and chaos under strong anticonformity and the presence of new equilibria when both conformity and anticonformity act at different variant frequencies, with and without selection. We also investigate the case of population subdivision with migration and find that the common claim that conformity can maintain between-group differences is not always true. Therefore, the effect of conformity on the evolution of cooperation by group selection may be more complicated than previously stated. Finally, using Feldman and Liberman's modifier approach, we investigate the conditions under which a rare modifier of the extent of conformity or the number of role models can invade a population. Understanding the dynamics of conformist- and anticonformist-biased transmission may have implications for research on human and nonhuman animal behavior, the evolution of cooperation, and frequency-dependent transmission in general.

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

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

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

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

Assessment of genetic diversity of the fat-tailed Dumba sheep of India by mitochondrial and microsatellite markers.

India has a centuries-old tradition of sheep production and breeding that accomplish economic, agricultural, and religious roles. In addition to the 44 registered sheep breeds, there is a fat-tailed sheep population referred to as Dumba. This study evaluated genetic variation in Dumba sheep and its differentiation from other Indian sheep breeds using mitochondrial DNA and genomic microsatellite loci. Haplotype and nucleotide diversity based on mitochondrial DNA analysis revealed substantially high maternal genetic diversity in Dumba sheep. Major ovine haplogroups A and B observed in sheep populations across the globe registered their presence in the Dumba sheep. The molecular genetic analysis using microsatellite markers also showed high measures of allele (10.125 ± 0.762) and gene diversity (0.749 ± 0.029). Results correspond to the non-bottleneck population that is near mutation-drift equilibrium despite some deficiency in the number of heterozygotes (FIS = 0.043 ± 0.059). Phylogenetic clustering confirmed Dumba to be a distinct population. Results of this study endow authorities with critical information imperative for sustainable utilization and conservation of Indian fat-tailed sheep, which is considered to be an untapped genetic resource contributing to the food security, livelihood, and economic sustainability of rural households in marginal areas of the country.

Whole-genome resequencing of Dorper and Hu sheep to reveal selection signatures associated with important traits.

Dorper and Hu sheep exhibit different characteristics in terms of reproduction, growth, and meat quality. Comparison of the genomes of two breeds help to reveal important genomic information. In this study, whole genome resequencing of 30 individuals (Dorper, DB and Hu sheep, HY) identified 15,108,125 single nucleotide polymorphisms (SNPs). Population differentiation (Fst) and cross population extended haplotype homozygosity (XP-EHH) were performed for selective signal analysis. In total, 106 and 515 overlapped genes were present in both the Fst results and XP-EHH results in HY vs DB and in DB vs HY, respectively. In HY vs DB, 106 genes were enriched in 12 GO terms and 83 KEGG pathways, such as ATP binding (GO:0005524) and PI3K-Akt signaling pathway (oas04151). In DB vs HY, 515 genes were enriched in 109 GO terms and 215 KEGG pathways, such as skeletal muscle cell differentiation (GO:0035914) and MAPK signaling pathway (oas04010). According to the annotation results, we identified a series of candidate genes associated with reproduction (UNC5C, BMPR1B, and GLIS1), meat quality (MECOM, MEF2C, and MYF6), and immunity (GMDS, GALK1, and ITGB4). Our investigation has uncovered genomic information for important traits in sheep and provided a basis for subsequent studies of related traits.

Discovery of SNP Molecular Markers and Candidate Genes Associated with Sacbrood Virus Resistance in Apis cerana cerana Larvae by Whole-Genome Resequencing.

Sacbrood virus (SBV) is a significant problem that impedes brood development in both eastern and western honeybees. Whole-genome sequencing has become an important tool in researching population genetic variations. Numerous studies have been conducted using multiple techniques to suppress SBV infection in honeybees, but the genetic markers and molecular mechanisms underlying SBV resistance have not been identified. To explore single nucleotide polymorphisms (SNPs), insertions, deletions (Indels), and genes at the DNA level related to SBV resistance, we conducted whole-genome resequencing on 90 Apis cerana cerana larvae raised in vitro and challenged with SBV. After filtering, a total of 337.47 gigabytes of clean data and 31,000,613 high-quality SNP loci were detected in three populations. We used ten databases to annotate 9359 predicted genes. By combining population differentiation index (FST) and nucleotide polymorphisms (π), we examined genome variants between resistant (R) and susceptible (S) larvae, focusing on site integrity (INT < 0.5) and minor allele frequency (MAF < 0.05). A selective sweep analysis with the top 1% and top 5% was used to identify significant regions. Two SNPs on the 15th chromosome with GenBank KZ288474.1_322717 (Guanine > Cytosine) and KZ288479.1_95621 (Cytosine > Thiamine) were found to be significantly associated with SBV resistance based on their associated allele frequencies after SNP validation. Each SNP was authenticated in 926 and 1022 samples, respectively. The enrichment and functional annotation pathways from significantly predicted genes to SBV resistance revealed immune response processes, signal transduction mechanisms, endocytosis, peroxisomes, phagosomes, and regulation of autophagy, which may be significant in SBV resistance. This study presents novel and useful SNP molecular markers that can be utilized as assisted molecular markers to select honeybees resistant to SBV for breeding and that can be used as a biocontrol technique to protect honeybees from SBV.

Population Structure Limits Parallel Evolution in Sticklebacks.

Population genetic theory predicts that small effective population sizes (Ne) and restricted gene flow limit the potential for local adaptation. In particular, the probability of evolving similar phenotypes based on shared genetic mechanisms (i.e., parallel evolution), is expected to be reduced. We tested these predictions in a comparative genomic study of two ecologically similar and geographically codistributed stickleback species (viz. Gasterosteus aculeatus and Pungitius pungitius). We found that P. pungitius harbors less genetic diversity and exhibits higher levels of genetic differentiation and isolation-by-distance than G. aculeatus. Conversely, G. aculeatus exhibits a stronger degree of genetic parallelism across freshwater populations than P. pungitius: 2,996 versus 379 single nucleotide polymorphisms located within 26 versus 9 genomic regions show evidence of selection in multiple freshwater populations of G. aculeatus and P. pungitius, respectively. Most regions involved in parallel evolution in G. aculeatus showed increased levels of divergence, suggestive of selection on ancient haplotypes. In contrast, haplotypes involved in freshwater adaptation in P. pungitius were younger. In accordance with theory, the results suggest that connectivity and genetic drift play crucial roles in determining the levels and geographic distribution of standing genetic variation, providing evidence that population subdivision limits local adaptation and therefore also the likelihood of parallel evolution.

Increased Floral Rewards due to Local Adaptation Drives Plant Ecological Speciation via Learned Preferences of Pollinators.

AbstractIn animal-pollinated plants, the growth environment and pollination environment are two important agents of natural selection. However, their simultaneous effects on plant speciation remain underexplored. Here, we report a theoretical finding that if plants' local adaptation to the growth environment increases their floral rewards for pollinators, it can strongly facilitate ecological speciation in plants. We consider two evolving plant traits, vegetative and floral signal traits, in a population genetic model for two plant populations under divergent selection from different growth environments. The vegetative trait determines plants' local adaptation. Locally adapted plants reward pollinators better than maladapted plants. By associative learning, pollinators acquire learned preferences for floral signal traits expressed by better-rewarding plants. If pollinators' learned preferences become divergent between populations, floral signal divergence occurs and plants develop genetic associations between vegetative and floral signal traits, leading to ecological speciation via a two-allele mechanism. Interestingly, speciation is contingent on whether novel floral signal variants arise before or after plant populations become locally adapted to the growth environment. Our results suggest that simultaneous selection from growth and pollination environments might be important for the ecological speciation of animal-pollinated plants.

More opportunities more species: Pleistocene differentiation and northward expansion of an evergreen broad-leaved tree species Machilus thunbergii (Lauraceae) in Southeast China.

BACKGROUND: The broad continuum between tropical and temperate floras in Eastern Asia (EAS) are thought to be one of the main factors responsible for a prominent species diversity anomaly of temperate plants between EAS and eastern North America (ENS). However, how the broad continuum and niche evolution between tropical and temperate floras in EAS contributes to lineage divergence and species diversity remains largely unknown. RESULTS: Population genetic structure, demography, and determinants of genetic structure [i.e., isolation-by-distance (IBD), isolation-by-resistance (IBR), and isolation-by-environment (IBE)] of Machilus thunbergii Sieb. et Zucc. (Lauraceae) were evaluated by examining sequence variation of ten low-copy nuclear genes across 43 populations in southeast China. Climatic niche difference and potential distributions across four periods (Current, mid-Holocene, the last glacial maximum, the last interglacial) of two genetic clusters were determined by niche modelling. North and south clusters of populations in M. thunbergii were revealed and their demarcation line corresponds well with the northern boundary of tropical zone in China of Zhu & Wan. The divergence time between the clusters and demographic expansion of M. thunbergii occurred after the mid-Pleistocene climate transition (MPT, 0.8-1.2 Ma). Migration rates between clusters were asymmetrical, being much greater from north to south than the reverse. Significant effects of IBE, but non-significant effects of IBD and IBR on population genetic divergence were detected. The two clusters have different ecological niches and require different temperature regimes. CONCLUSIONS: The north-south genetic differentiation may be common across the temperate-tropical boundary in southeast China. Divergent selection under different temperature regimes (possibly above and below freezing temperature in winter) could account for this divergence pattern. The broad continuum between tropical and temperate floras in EAS may have provided ample opportunities for tropical plant lineages to acquire freezing tolerance and to colonize the temperate regions during the late-Cenozoic global cooling. Our findings shed deeper insights into the high temperate plant species diversity in EAS.

Population structure in Neotropical plants: Integrating pollination biology, topography and climatic niches.

Animal pollinators mediate gene flow among plant populations, but in contrast to well-studied topographic and (Pleistocene) environmental isolating barriers, their impact on population genetic differentiation remains largely unexplored. Comparing how these multifarious factors drive microevolutionary histories is, however, crucial for better resolving macroevolutionary patterns of plant diversification. Here we combined genomic analyses with landscape genetics and niche modelling across six related Neotropical plant species (424 individuals across 33 localities) differing in pollination strategy to test the hypothesis that highly mobile (vertebrate) pollinators more effectively link isolated localities than less mobile (bee) pollinators. We found consistently higher genetic differentiation (FST ) among localities of bee- than vertebrate-pollinated species with increasing geographical distance, topographic barriers and historical climatic instability. High admixture among montane populations further suggested relative climatic stability of Neotropical montane forests during the Pleistocene. Overall, our results indicate that pollinators may differentially impact the potential for allopatric speciation, thereby critically influencing diversification histories at macroevolutionary scales.

Geological and climatic influences on population differentiation of the Phrynocephalus vlangalii species complex (Sauria: Agamidae) in the northern Qinghai-Tibet Plateau.

Extremely heterogeneous topography and complex paleoclimatic history of the Qinghai-Tibet Plateau (QTP) have a key role in promoting genetic divergence among populations and lineage/species formation. Here, we sequenced one nuclear and three mitochondrial markers of 532 individuals from the entire range of the Phrynocephalus vlangalii species complex including two species, P. putjatai and P. vlangalii, endemic to the northern QTP. We integrated multilocus phylogeny, demographic analysis and geographic barrier detection to evaluate the population structure and dynamics. We found a new mitochondrial clade (PV-I) in the Gonghe County population of P. vlangalii, partial mitochondrial DNA replacement within P. vlangalii and complete mitochondrial DNA replacement between P. putjatai and P. vlangalii. Neutrality test, mismatch distribution analysis and Extended Bayesian Skyline Plot (EBSP) analysis all supported a significant expansion of the Qaidam Basin population of P. vlangalii (PV-II-2) from 0.091 to 0.026 Ma after Penultimate Glaciation. The uplift of the Arjin and Anyemanqen Mountains during the Kunhuang Movement (∼1.2 Ma) split populations of P. vlangalii in Akesai, Qaidam Basin and source of the Yellow River. The uplift of the Elashan Mountains during the second phase of the Qingzang Movement (∼2.5 Ma) contributed to the divergence of the Gonghe County population of P. vlangalii from other conspecific populations. The third phase of the Qingzang Movement (∼1.7 Ma) contributed to the divergence of the Xinghai population of P. vlangalii from P. putjatai and to the divergence of the northern populations of P. putjatai from the southern conspecific populations. Our data support the idea that the geological and climatic changes following the orogeny of the QTP may have promoted population differentiation and shaped the current population patterns of the P. vlangalii species complex in the northeastern QTP.

Genomic diversity and breed composition of Vietnamese smallholder dairy cows.

Vietnamese smallholder dairy cows (VDC) are the result of crossbreeding between different zebu (ZEB) and taurine dairy breeds through many undefined generations. Thus, the predominant breed composition of VDC is currently unknown. This study aimed to evaluate the level of genetic diversity and breed composition of VDC. The SNP data of 344 animals from 32 farms located across four dairy regions of Vietnam were collected and merged with genomic reference data, which included three ZEB breeds: Red Sindhi, Sahiwal and Brahman, three taurine breeds: Holstein (HOL), Jersey (JER) and Brown Swiss (BSW), and a composite breed: Chinese Yellow cattle. Diversity and admixture analyses were applied to the merged data set. The VDC were not excessively inbred, as indicated by very low inbreeding coefficients (Wright's FIS ranged from -0.017 to 0.003). The genetic fractions in the test herds suggested that the VDC are primarily composed of HOL (85.0%); however, JER (6.0%), BSW 5.3%) and ZEB (4.5%) had also contributed. Furthermore, major genotype groupings in the test herds were pure HOL (48%), B3:15/16HOL_1/16ZEB (22%) and B2:7/8HOL_1/8ZEB (12%). The genetic makeup of the VDC is mainly components of various dairy breeds but also has a small percentage of ZEB; thus, the VDC could be a good genetic base for selecting high milk-producing cows with some degree of adaptation to tropical conditions.

Prioritising positively selected variants in whole-genome sequencing data using FineMAV.

BACKGROUND: In population genomics, polymorphisms that are highly differentiated between geographically separated populations are often suggestive of Darwinian positive selection. Genomic scans have highlighted several such regions in African and non-African populations, but only a handful of these have functional data that clearly associates candidate variations driving the selection process. Fine-Mapping of Adaptive Variation (FineMAV) was developed to address this in a high-throughput manner using population based whole-genome sequences generated by the 1000 Genomes Project. It pinpoints positively selected genetic variants in sequencing data by prioritizing high frequency, population-specific and functional derived alleles. RESULTS: We developed a stand-alone software that implements the FineMAV statistic. To graphically visualise the FineMAV scores, it outputs the statistics as bigWig files, which is a common file format supported by many genome browsers. It is available as a command-line and graphical user interface. The software was tested by replicating the FineMAV scores obtained using 1000 Genomes Project African, European, East and South Asian populations and subsequently applied to whole-genome sequencing datasets from Singapore and China to highlight population specific variants that can be subsequently modelled. The software tool is publicly available at https://github.com/fadilla-wahyudi/finemav . CONCLUSIONS: The software tool described here determines genome-wide FineMAV scores, using low or high-coverage whole-genome sequencing datasets, that can be used to prioritize a list of population specific, highly differentiated candidate variants for in vitro or in vivo functional screens. The tool displays these scores on the human genome browsers for easy visualisation, annotation and comparison between different genomic regions in worldwide human populations.