Effective Population Research Articles

Benchmarking and Optimization of Methods for the Detection of Identity-By-Descent in High-Recombining Plasmodium falciparum Genomes.

Genomic surveillance is crucial for identifying at-risk populations for targeted malaria control and elimination. Identity-by-descent (IBD) is increasingly being used in Plasmodium population genomics to estimate genetic relatedness, effective population size (N e ), population structure, and signals of positive selection. Despite its potential, a thorough evaluation of IBD segment detection tools for species with high recombination rates, such as P. falciparum, remains absent. Here, we perform comprehensive benchmarking of IBD callers - probabilistic (hmmIBD, isoRelate), identity-by-state-based (hap-IBD, phased IBD) and others (Refined IBD) - using population genetic simulations tailored for high recombination, and IBD quality metrics at both the IBD segment level and the IBD-based downstream inference level. Our results demonstrate that low marker density per genetic unit, related to high recombination relative to mutation, significantly compromises the accuracy of detected IBD segments. In genomes with high recombination rates resembling P. falciparum, most IBD callers exhibit high false negative rates for shorter IBD segments, which can be partially mitigated through optimization of IBD caller parameters, especially those related to marker density. Notably, IBD detected with optimized parameters allows for more accurate capture of selection signals and population structure; IBD-based N e inference is very sensitive to IBD detection errors, with IBD called from hmmIBD uniquely providing less biased estimates of N e in this context. Validation with empirical data from the MalariaGEN Pf 7 database, representing different transmission settings, corroborates these findings. We conclude that context-specific evaluation and parameter optimization are essential for accurate IBD detection in high-recombining species and recommend hmmIBD for quality-sensitive analysis, such as estimation of N e in these species. Our optimization and high-level benchmarking methods not only improve IBD segment detection in high-recombining genomes but also enhance overall genomic analysis, paving the way for more accurate genomic surveillance and targeted intervention strategies for malaria.

Evaluation of heritability partitioning approaches in livestock populations

BackgroundHeritability partitioning approaches estimate the contribution of different functional classes, such as coding or regulatory variants, to the genetic variance. This information allows a better understanding of the genetic architecture of complex traits, including complex diseases, but can also help improve the accuracy of genomic selection in livestock species. However, methods have mainly been tested on human genomic data, whereas livestock populations have specific characteristics, such as high levels of relatedness, small effective population size or long-range levels of linkage disequilibrium.ResultsHere, we used data from 14,762 cows, imputed at the whole-genome sequence level for 11,537,240 variants, to simulate traits in a typical livestock population and evaluate the accuracy of two state-of-the-art heritability partitioning methods, GREML and a Bayesian mixture model. In simulations where a single functional class had increased contribution to heritability, we observed that the estimators were unbiased but had low precision. When causal variants were enriched in variants with low (< 0.05) or high (> 0.20) minor allele frequency or low (below 1st quartile) or high (above 3rd quartile) linkage disequilibrium scores, it was necessary to partition the genetic variance into multiple classes defined on the basis of allele frequencies or LD scores to obtain unbiased results. When multiple functional classes had variable contributions to heritability, estimators showed higher levels of variation and confounding between certain categories was observed. In addition, estimators from small categories were particularly imprecise. However, the estimates and their ranking were still informative about the contribution of the classes. We also demonstrated that using methods that estimate the contribution of a single category at a time, a commonly used approach, results in an overestimation. Finally, we applied the methods to phenotypes for muscular development and height and estimated that, on average, variants in open chromatin regions had a higher contribution to the genetic variance (> 45%), while variants in coding regions had the strongest individual effects (> 25-fold enrichment on average). Conversely, variants in intergenic or intronic regions showed lower levels of enrichment (0.2 and 0.6-fold on average, respectively).ConclusionsHeritability partitioning approaches should be used cautiously in livestock populations, in particular for small categories. Two-component approaches that fit only one functional category at a time lead to biased estimators and should not be used.

Epidemiology of carbapenem-resistant Klebsiella pneumoniae in China and the evolving trends of predominant clone ST11: a multicentre, genome-based study.

Carbapenem-resistant Klebsiella pneumoniae (CRKP) is a major nosocomial infectious pathogen with rapidly increasing prevalence. The genomic epidemiological characteristics of CRKP nationwide, especially the evolving trends within the predominant clones, should be evaluated clearly. We collected 3415 K. pneumoniae strains from 28 hospitals across China. Antimicrobial susceptibility testing and WGS were performed. Subsequent genomic analyses, including sequence typing, K-locus (KL) identification, antimicrobial resistance gene screening, and virulence score assessment were performed. The phylogenetic relationship of clonal group 11 was determined based on core-genome analysis, and the presence of the pLVPK-like virulence plasmid in ST11 isolates was confirmed using plasmid core-gene analysis. Additionally, the trends of the ST11 lineage with different KL types on a global scale were investigated using Beast2. Of the K. pneumoniae strains, 708 were identified as CRKP isolates (20.7%), of which 97.7% were MDR. ST11 was the predominant clone, and KPC-2 was the prevalent carbapenemase in China, although the prevalence of specific clones and carbapenemases varied by geographic region. Among ST11 isolates, KL47 and KL64 were the predominant KL types, and KL64 gradually replaced KL47, with a higher percentage of KL64 isolates harbouring the pLVPK-like plasmid. Global genome data showed a significant increase in the effective population size of KL64 over the last 5 years. The prevalence of CRKP was very high in certain regions in China. The increasing convergence of virulence and resistance, particularly in ST11-KL64 isolates, should be given more attention and further investigation.

ONeSAMP 3.0: estimation of effective population size via single nucleotide polymorphism data from one population.

The genetic effective size (Ne) is arguably one of the most important characteristics of a population as it impacts the rate of loss of genetic diversity. Methods that estimate Ne are important in population and conservation genetic studies as they quantify the risk of a population being inbred or lacking genetic diversity. Yet there are very few methods that can estimate the Ne from data from a single population and without extensive information about the genetics of the population, such as a linkage map, or a reference genome of the species of interest. We present ONeSAMP 3.0, an algorithm for estimating Ne from single nucleotide polymorphism data collected from a single population sample using approximate Bayesian computation and local linear regression. We demonstrate the utility of this approach using simulated Wright-Fisher populations, and empirical data from five endangered Channel Island fox (Urocyon littoralis) populations to evaluate the performance of ONeSAMP 3.0 compared to a commonly used Ne estimator. Our results show that ONeSAMP 3.0 is broadly applicable to natural populations and is flexible enough that future versions could easily include summary statistics appropriate for a suite of biological and sampling conditions. ONeSAMP 3.0 is publicly available under the GNU General Public License at https://github.com/AaronHong1024/ONeSAMP_3.

Illuminating Genetic Diversity and Selection Signatures in Matou Goats through Whole-Genome Sequencing Analysis.

(1) Background: Matou goats, native to Hunan and Hubei provinces in China, are renowned for their exceptional meat and skin quality. However, a comprehensive whole-genome-based exploration of the genetic architecture of this breed is scant in the literature. (2) Methods: To address this substantial gap, we used whole-genome sequences of 20 Matou goats and compared them with published genomic data of 133 goats of different breeds across China. This comprehensive investigation sought to assess genetic diversity, population structure, and the presence of genomic selection signals. (3) Results: The whole genome of Matou goat populations yielded a substantial catalog of over 19 million single nucleotide polymorphisms (SNPs), primarily distributed within intergenic and intron regions. The phylogenetic tree analysis revealed distinct clades corresponding to each goat population within the dataset. Notably, this analysis positioned Matou goats in a closer genetic affinity with Guizhou White goats, compared to other recognized goat breeds. This observation was corroborated by principal component analysis (PCA) and admixture analysis. Remarkably, Matou goats exhibited diminished genetic diversity and a notable degree of inbreeding, signifying a reduced effective population size. Moreover, the study employed five selective sweep detection methods (including PI, CLR, PI-Ratio, Fst, and XP-EHH) to screen top signal genes associated with critical biological functions, encompassing cardiomyocytes, immunity, coat color, and meat quality. (4) Conclusions: In conclusion, this study significantly advances our understanding of the current genetic landscape and evolutionary dynamics of Matou goats. These findings underscore the importance of concerted efforts in resource conservation and genetic enhancement for this invaluable breed.

Mitonuclear compatibility is maintained despite relaxed selection on male mitochondrial DNA in bivalves with doubly uniparental inheritance.

Mitonuclear coevolution is common in eukaryotes, but bivalve lineages that have doubly uniparental inheritance (DUI) of mitochondria may be an interesting example. In this system, females transmit mtDNA (F mtDNA) to all offspring, while males transmit a different mtDNA (M mtDNA) solely to their sons. Molecular evolution and functional data suggest oxidative phosphorylation (OXPHOS) genes encoded in M mtDNA evolve under relaxed selection due to their function being limited to sperm only (vs. all other tissues for F mtDNA). This has led to the hypothesis that mitonuclear coevolution is less important for M mtDNA. Here, we use comparative phylogenetics, transcriptomics, and proteomics to understand mitonuclear interactions in DUI bivalves. We found nuclear OXPHOS proteins coevolve and maintain compatibility similarly with both F and M mtDNA OXPHOS proteins. Mitochondrial recombination did not influence mitonuclear compatibility and nuclear-encoded OXPHOS genes were not upregulated in tissues with M mtDNA to offset dysfunction. Our results support that selection maintains mitonuclear compatibility with F and M mtDNA despite relaxed selection on M mtDNA. Strict sperm transmission, lower effective population size, and higher mutation rates may explain the evolution of M mtDNA. Our study highlights that mitonuclear coevolution and compatibility may be broad features of eukaryotes.

Repeated horizontal acquisition of lagriamide-producing symbionts in Lagriinae beetles.

Microbial symbionts associate with multicellular organisms on a continuum from facultative associations to mutual codependency. In some of the oldest intracellular symbioses there is exclusive vertical symbiont transmission, and co-diversification of symbiotic partners over millions of years. Such symbionts often undergo genome reduction due to low effective population sizes, frequent population bottlenecks, and reduced purifying selection. Here, we describe multiple independent acquisition events of closely related defensive symbionts followed by genome erosion in a group of Lagriinae beetles. Previous work in Lagria villosa revealed the dominant genome-eroded symbiont of the genus Burkholderia produces the antifungal compound lagriamide and protects the beetle's eggs and larvae from antagonistic fungi. Here, we use metagenomics to assemble 11 additional genomes of lagriamide-producing symbionts from seven different host species within Lagriinae from five countries, to unravel the evolutionary history of this symbiotic relationship. In each host species, we detected one dominant genome-eroded Burkholderia symbiont encoding the lagriamide biosynthetic gene cluster (BGC). Surprisingly, however, we did not find evidence for host-symbiont co-diversification, or for a monophyly of the lagriamide-producing symbionts. Instead, our analyses support at least four independent acquisition events of lagriamide-encoding symbionts and subsequent genome erosion in each of these lineages. By contrast, a clade of plant-associated relatives retained large genomes but secondarily lost the lagriamide BGC. In conclusion, our results reveal a dynamic evolutionary history with multiple independent symbiont acquisitions characterized by high degree of specificity. They highlight the importance of the specialized metabolite lagriamide for the establishment and maintenance of this defensive symbiosis.

Molecular Assessment of Genetic Diversity and Bottleneck in Hair Goat Reared in Türkiye

Objective: Being the most preferred and geographically distributed in Türkiye, the Hair goat was screened at a molecular level to evaluate genetic diversity and population structure via microsatellite DNA markers. This paper also aimed to investigate the effects of genetic bottleneck to evaluate whether the Hair goat has maintained its effective population size in recent past. Material and Methods: A total of 411 Hair goats were sampled from farms participating in the "Hair Goat Breeding" project, initiated by the General Directorate of Agricultural Research and Policies in Aydın and Denizli provinces. Sampled animals were genotyped with 18 microsatellite loci to assess genetic diversity, population structure, and genetic bottleneck. Results: A total of 341 different alleles were observed across 18 microsatellite loci in which the highest number of alleles (26) and effective alleles (10.18) were detected in INRA005 and HSC loci, respectively. The average observed heterozygosity (0.73) was lower than the expected value (0.83), whereas all loci turned out to be highly informative (PIC&amp;gt;0.50). Factorial correspondence analysis separated animals into two groups, while a genetic admixture was detected between these groups. STRUCTURE analysis, on the other hand, confirmed that 411 animals were derived from three ancestral populations in which the third group is drawn due to admixed individuals. The Wilcoxon test and mode-shift indicator detected a lack of genetic bottleneck indicating that Hair goats reared in Türkiye have maintained their effective population size in recent past. Conclusion: This study validates that used microsatellite markers are highly polymorphic and could be utilized for revealing genetic diversity in different local goat breeds. The findings recovered in this study could be integrated into breeding and conservation programs, while further studies should adopt SNP array technologies and next-generation sequencing platforms to reveal deeper knowledge about the genetic diversity and population structure of Anatolian goat breeds.

Drivers of genomic diversity and phenotypic development in early phases of domestication in Hermetia illucens.

The black soldier fly (BSF), Hermetia illucens, has the ability to efficiently bioremediate organic waste into usable bio-compounds. Understanding the impact of domestication and mass rearing on fitness and production traits is therefore important for sustainable production. This study aimed to assess patterns of genomic diversity and its association to phenotypic development across early generations of mass rearing under two selection strategies: selection for greater larval mass (SEL lines) and no direct artificial selection (NS lines). Genome-wide single nucleotide polymorphism (SNP) data were generated using 2bRAD sequencing, while phenotypic traits relating to production and population fitness were measured. Declining patterns of genomic diversity were observed across three generations of captive breeding, with the lowest diversity recorded for the F3 generation of both selection lines, most likely due to founder effects. The SEL cohort displayed statistically significantly greater larval weight com the NS lines with pronounced genetic and phenotypic directional changes across generations. Furthermore, lower genetic and phenotypic diversity, particularly for fitness traits, were evident for SEL lines, illustrating the trade-off between selecting for mass and the resulting decline in population fitness. SNP-based heritability was significant for growth, but was low or non-significant for fitness traits. Genotype-phenotype correlations were observed for traits, but individual locus effect sizes where small and very few of these loci demonstrated a signature for selection. Pronounced genetic drift, due to small effective population sizes, is likely overshadowing the impacts of selection on genomic diversity and consequently phenotypic development. The results hold particular relevance for genetic management and selective breeding for BSF in future.

Ecological genetics of isolated loach populations indicate compromised adaptive potential

Many endangered species live in fragmented and isolated populations with low genetic variability, signs of inbreeding, and small effective population sizes - all features elevating their extinction risk. The flat-headed loach (Oreonectes platycephalus), a small noemacheilid fish, is widely across southern China, but only in the headwaters of hillstreams; as a result, they are spatially isolated from conspecific populations. We surveyed single nucleotide polymorphisms in 16 Hong Kong populations of O. platycephalus to determine whether loach populations from different streams were genetically isolated from each other, showed low levels of genetic diversity, signs of inbreeding, and had small contemporary effective population sizes. Estimates of average observed heterozygosity (HO = 0.0473), average weighted nucleotide diversity (πw = 0.0546) and contemporary effective population sizes (Ne = 10.2 ~ 129.8) were very low, and several populations showed clear signs of inbreeding as judged from relatedness estimates. The degree of genetic differentiation among populations was very high (average FST = 0.668), even over short geographic distances (<1.5 km), with clear patterns of isolation by distance. These results suggest that Hong Kong populations of O. platycephalus have experienced strong genetic drift and loss of genetic variability because sea-level rise after the last glaciation reduced connectedness among paleodrainages, isolating populations in headwaters. All this, together with the fact that the levels of genetic diversity and contemporary effective population sizes within O. platycephalus populations are lower than most other freshwater fishes, suggests that they face high local extinction risk and have limited capacity for future adaptation.

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

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

Genetic consequences in the southern African endemic seabream Polysteganus undulosus (Sparidae) after eight decades of overfishing

The Critically Endangered seventy-four seabream Polysteganus undulosus, a slow-growing sparid that forms spawning aggregations off South Africa, faced heavy exploitation from 1910 until a fishery moratorium was put in place in 1998. Utilising temporal samples from 1962/1963 (mid-collapse) and 2005/2006 (post-collapse), we assessed genetic diversity at six microsatellite loci. Amplification success for archived samples was low (43%), necessitating a rarefaction approach, revealing a 40% decrease in allelic diversity. Significant genetic differences between recent and archived samples confirmed the impact of overfishing. Simulation studies indicated that missing genotypes did not affect these tests, validating the genetic differences. Using a coalescent-based approach, a 10-fold decrease in effective population size (Ne) was estimated over this 43-year period (Ne = 43.88 to 4.76). Simulations provided corrected Ne of 36.86, accounting for missing genotypes that were responsible for inflated values. Assuming 25% of pristine levels in the 1960s, pristine Ne values ranged from 147 to 176 in 1910, suggesting a 96.8–98.3% genetic decline over 88 years of exploitation. This study emphasises severe genetic consequences of overfishing on P. undulosus diversity and effective population size. The results provide vital genetic baseline data for recovery assessment and conservation efforts, as well as an analytical framework to evaluate the stock decline using partial genotyping data gathered from degraded archived samples.

Estimation of effective number of breeders and effective population size in an abundant and heavily exploited marine teleost.

Obtaining reliable estimates of the effective number of breeders (N b) and generational effective population size (N e) for fishery-important species is challenging because they are often iteroparous and highly abundant, which can lead to bias and imprecision. However, recent advances in understanding of these parameters, as well as the development of bias correction methods, have improved the capacity to generate reliable estimates. We utilized samples of both single-cohort young of the year and mixed-age adults from two geographically and genetically isolated stocks of the Australasian snapper (Chrysophrys auratus) to investigate the feasibility of generating reliable N b and N e estimates for a fishery species. Snapper is an abundant, iteroparous broadcast spawning teleost that is heavily exploited by recreational and commercial fisheries. Employing neutral genome-wide SNPs and the linkage-disequilibrium method, we determined that the most reliable N b and N e estimates could be derived by genotyping at least 200 individuals from a single cohort. Although our estimates made from the mixed-age adult samples were generally lower and less precise than those based on a single cohort, they still proved useful for understanding relative differences in genetic effective size between stocks. The correction formulas applied to adjust for biases due to physical linkage of loci and age structure resulted in substantial upward modifications of our estimates, demonstrating the importance of applying these bias corrections. Our findings provide important guidelines for estimating N b and N e for iteroparous species with large populations. This work also highlights the utility of samples originally collected for stock structure and stock assessment work for investigating genetic effective size in fishery-important species.

Pitfalls and windfalls of detecting demographic declines using population genetics in long-lived species.

Detecting recent demographic changes is a crucial component of species conservation and management, as many natural populations face declines due to anthropogenic habitat alteration and climate change. Genetic methods allow researchers to detect changes in effective population size (Ne) from sampling at a single timepoint. However, in species with long lifespans, there is a lag between the start of a decline in a population and the resulting decrease in genetic diversity. This lag slows the rate at which diversity is lost, and therefore makes it difficult to detect recent declines using genetic data. However, the genomes of old individuals can provide a window into the past, and can be compared to those of younger individuals, a contrast that may help reveal recent demographic declines. To test whether comparing the genomes of young and old individuals can help infer recent demographic bottlenecks, we use forward-time, individual-based simulations with varying mean individual lifespans and extents of generational overlap. We find that age information can be used to aid in the detection of demographic declines when the decline has been severe. When average lifespan is long, comparing young and old individuals from a single timepoint has greater power to detect a recent (within the last 50 years) bottleneck event than comparing individuals sampled at different points in time. Our results demonstrate how longevity and generational overlap can be both a hindrance and a boon to detecting recent demographic declines from population genomic data.

Genetic diversity and shallow genetic differentiation of the endangered scaly-sided merganser Mergus squamatus.

Examining patterns of genetic diversity are crucial for conservation planning on endangered species, while inferring the underlying process of recent anthropogenic habitat modifications in the context potential long-term demographic changes remains challenging. The globally endangered scaly-sided merganser (SSME), Mergus squamatus, is endemic to a narrow range in Northeast Asia, and its population has recently been contracted into two main breeding areas. Although low genetic diversity has been suggested in the Russian population, the genetic status and demographic history of these individuals have not been fully elucidated. We therefore examined the genetic diversity and structure of the breeding populations of the SSME and investigated the relative importance of historical and recent demographic changes to the present-day pattern of genetic diversity. Using 10 nuclear microsatellite (SSR) markers and mitochondrial DNA (mtDNA) control region sequences, we found limited female-inherited genetic diversity and a high level of nuclear genetic diversity. In addition, analysis of both markers consistently revealed significant but weak divergence between the breeding populations. Inconsistent demographic history parameters calculated from mtDNA and bottleneck analysis results based on SSR suggested a stable historical effective population size. By applying approximate Bayesian computation, it was estimated that populations started to genetically diverge from each other due to recent fragmentation events caused by anthropogenic effects rather than isolation during Last Glacial Maximum (LGM) and post-LGM recolonization. These results suggest that limited historical population size and shallow evolutionary history may be potential factors contributing to the contemporary genetic diversity pattern of breeding SSME populations. Conservation efforts should focus on protecting the current breeding habitats from further destruction, with priority given to both the Russian and Chinese population, as well as restoring the connected suitable breeding grounds.

Genetic diversity and population structure of wild and artificially selected populations of the ridgetail white prawn (Exopalaemon carinicauda) using consolidated mitochondrial cytochrome b and microsatellite markers

The ridgetail white prawn (Exopalaemon carinicauda) belongs to the family Palaemonidae and genus Exopalaemon. Owing to its significant economic value and excellent nutritional properties, it is an important species for the development of saltwater aquaculture models in China. In this study, 160 individuals were collected from four wild populations and one artificially selected population using consolidated mitochondrial cytochrome b (mt-cyt b), and microsatellite molecular markers. The results showed that the hatcheries had lower genetic diversity than the wild samples, both at the levels of mt-cyt b (haplotype diversity: 0.123 VS 0.572; nucleotide diversity: 0.0003 VS 0.0014) and microsatellite (expected heterozygosity: 0.672 VS 0.752). Among the four wild populations, the Ningbo population sampled from the East China Sea showed higher genetic diversity in both mt-cyt b and microsatellites. A pairwise comparison of genetic differentiation coefficients revealed that the artificially selected population had not genetically diverged from the Cangzhou and Rizhao populations. However, the Ningbo population exhibited moderate to high genetic differentiation (compared to other populations), whereas the Zhanjiang population showed a high level of genetic differentiation, indicating a distinct geographical distribution pattern. Structural analysis revealed the presence of three lineage clusters in all current populations, and the artificially selected population predominantly consisted of lineage 3, accounting for 71.68 % of the population. Currently, the effective population size of hatcheries was <50, rendering them susceptible to the genetic threat of inbreeding depression. These findings provide a scientific basis for the conservation of wild germplasm resources of E. carinicauda and inform future breeding and management strategies for improved varieties.

Genomic Characterization of Local Croatian Sheep Breeds-Effective Population Size, Inbreeding & Signatures of Selection.

The Istrian (IS) and the Pag sheep (PS) are local Croatian breeds which provide significant income for the regional economy and have a cultural and traditional importance for the inhabitants. The aim of this study was to estimate some important population specific genetic parameters in IS (N = 1293) and PS (N = 2637) based on genome wide SNPs. Estimates of linkage disequilibrium effective population size (Ne) evidenced more genetic variability in PS (Ne = 838) compared to IS (Ne = 197), regardless of historical time (both recent and ancient genetic variability). The discrepancy in the recent genetic variability between these breeds was additionally confirmed by the estimates of genomic inbreeding (FROH), which was estimated to be notably higher in IS (FROH>2 = 0.062) than in PS (FROH>2 = 0.029). The average FROH2-4, FROH4-8, FROH8-16, and FROH>16 were 0.26, 1.65, 2.14, and 3.72 for IS and 0.22, 0.61, 0.75, and 1.58 for PS, thus evidencing a high contribution of recent inbreeding in the overall inbreeding. One ROH island with > 30% of SNP incidence in ROHs was detected in IS (OAR6; 34,253,440-38,238,124 bp) while there was no ROH islands detected in PS. Seven genes (CCSER1, HERC3, LCORL, NAP1L5, PKD2, PYURF, and SPP1) involved in growth, feed intake, milk production, immune responses, and resistance were associated with the found autozygosity. The results of this study represent the first comprehensive insight into genomic variability of these two Croatian local sheep breeds and will serve as a baseline for setting up the most promising strategy of genomic Optimum Contribution Selection.

Continent-wide genomic analysis of the African buffalo (Syncerus caffer)

The African buffalo (Syncerus caffer) is a wild bovid with a historical distribution across much of sub-Saharan Africa. Genomic analysis can provide insights into the evolutionary history of the species, and the key selective pressures shaping populations, including assessment of population level differentiation, population fragmentation, and population genetic structure. In this study we generated the highest quality de novo genome assembly (2.65 Gb, scaffold N50 69.17 Mb) of African buffalo to date, and sequenced a further 195 genomes from across the species distribution. Principal component and admixture analyses provided little support for the currently described four subspecies. Estimating Effective Migration Surfaces analysis suggested that geographical barriers have played a significant role in shaping gene flow and the population structure. Estimated effective population sizes indicated a substantial drop occurring in all populations 5-10,000 years ago, coinciding with the increase in human populations. Finally, signatures of selection were enriched for key genes associated with the immune response, suggesting infectious disease exert a substantial selective pressure upon the African buffalo. These findings have important implications for understanding bovid evolution, buffalo conservation and population management.

Habitat quality or quantity? Niche marginality across 21 plants and animals suggests differential responses between highland and lowland species to past climatic changes

Climatic changes can affect species distributions, population abundance, and evolution. Such organismal responses could be determined by the amount and quality of available habitats, which can vary independently. In this study, we assessed changes in habitat quantity and quality independently to generate explicit predictions of the species' responses to climatic changes between Last Glacial Maximum (LGM) and present day. We built ecological niche models for genetic groups within 21 reptile, mammal, and plant taxa from the Baja California peninsula inhabiting lowland or highland environments. Significant niche divergence was detected for all clades within species, along with significant differences in the niche breadth and area of distribution between northern and southern clades. We quantified habitat quantity from the distribution models, and most clades showed a reduction in distribution area towards LGM. Further, niche marginality (used as a measure of habitat quality) was higher during LGM for most clades, except for northern highland species. Our results suggest that changes in habitat quantity and quality can affect organismal responses independently. This allows the prediction of genomic signatures associated with changes in effective population size and selection pressure that could be explicitly tested from our models.

Chromosome-level genome assembly of a cliff plant Taihangia rupestris var. ciliata provides insights into its adaptation and demographic history

BackgroundCliffs are recognized as one of the most challenging environments for plants, characterized by harsh conditions such as drought, infertile soil, and steep terrain. However, they surprisingly host ancient and diverse plant communities and play a crucial role in protecting biodiversity. The Taihang Mountains, which act as a natural boundary in eastern China, support a rich variety of plant species, including many unique to cliff habitats. However, it is little known how cliff plants adapt to harsh habitats and the demographic history in this region.ResultsTo better understand the demographic history and adaptation of cliff plants in this area, we analyzed the chromosome-level genome of a representative cliff plant, T. rupestris var. ciliata, which has a genome size of 769.5 Mb, with a scaffold N50 of 104.92 Mb. The rapid expansion of transposable elements may have contributed to the increasing genome and its ability to adapt to unique and challenging cliff habitats. Comparative analysis of the genome evolution between Taihangia and non-cliff plants in Rosaceae revealed a significant expansion of gene families associated with oxidative phosphorylation, which is likely a response to the abiotic stresses faced by cliff plants. This expansion may explain the long-term adaptation of Taihangia to harsh cliff environments. The effective population size of the two varieties has continuously decreased due to climatic fluctuations during the Quaternary period. Furthermore, significant differences in gene expression between the two varieties may explain the varied leaf phenotypes and adaptations to harsh conditions in different natural distributions.ConclusionOur study highlights the extraordinary adaptation of T. rupestris var. ciliata, shedding light on the evolution of cliff plants worldwide.