Search

In local chickens targeted for niche markets, genotyping costs are relatively high due to the small population size and diverse breeding goals. The single-step genomic best linear unbiased prediction (ssGBLUP) model, which combines pedigree and genomic information, has been introduced to increase the accuracy of genomic estimated breeding value (GEBV). Therefore, this model may be more beneficial than the genomic BLUP (GBLUP) model for genomic selection in local chickens. Additionally, the single-step genome-wide association study (ssGWAS) can be used to extend the ssGBLUP model results to animals with available phenotypic information but without genotypic data. In this study, we compared the accuracy of (G)EBVs using the pedigree-based BLUP (PBLUP), GBLUP, and ssGBLUP models. Moreover, we conducted single-SNP GWAS (SNP-GWAS), GBLUP-GWAS, and ssGWAS methods to identify genes associated with egg production traits in the NCHU-G101 chicken to understand the feasibility of using genomic selection in a small population. The average prediction accuracy of (G)EBV for egg production traits using the PBLUP, GBLUP, and ssGBLUP models is 0.536, 0.531, and 0.555, respectively. In total, 22 suggestive- and 5% Bonferroni genome-wide significant-level SNPs for total egg number (EN), average laying rate (LR), average clutch length, and total clutch number are detected using 3 GWAS methods. These SNPs are mapped onto Gallus gallus chromosomes (GGA) 4, 6, 10, 18, and 25 in NCHU-G101 chicken. Furthermore, through SNP-GWAS and ssGWAS methods, we identify 2 genes on GGA4 associated with EN and LR: ENSGALG00000023172 and PPARGC1A. In conclusion, the ssGBLUP model demonstrates superior prediction accuracy, performing on average 3.41% than the PBLUP model. The implications of our gene results may guide future selection strategies for Taiwan Country chickens. Our results highlight the applicability of the ssGBLUP model for egg production traits selection in a small population, specifically NCHU-G101 chicken in Taiwan., Competing Interests: DISCLOSURES The authors declare no conflicts of interest., (Copyright © 2024. Published by Elsevier Inc.)

Objective: Phenotypic misclassification in genetic association analyses can impact the accuracy of PRS-based prediction models. The bias reduction method proposed by Tong et al. (2019) has demonstrated its efficacy in reducing the effects of bias on the estimation of association parameters between genotype and phenotype while minimizing variance by employing chart reviews on a subset of the data for validating phenotypes, however its improvement of subsequent PRS prediction accuracy remains unclear. Our study aims to fill this gap by assessing the performance of simulated PRS models and estimating the optimal number of chart reviews needed for validation., Methods: To comprehensively assess the efficacy of the bias reduction method proposed by Tong et al. in enhancing the accuracy of PRS-based prediction models, we simulated each phenotype under different correlation structures (an independent model, a weakly correlated model, a strongly correlated model) and introduced error-prone phenotypes using two distinct error mechanisms (differential and non-differential phenotyping errors). To facilitate this, we used genotype and phenotype data from 12 case-control datasets in the Alzheimer's Disease Genetics Consortium (ADGC) to produce simulated phenotypes. The evaluation included analyses across various misclassification rates of original phenotypes as well as quantities of validation set. Additionally, we determined the median threshold, identifying the minimal validation size required for a meaningful improvement in the accuracy of PRS-based predictions across a broad spectrum., Results: This simulation study demonstrated that incorporating chart review does not universally guarantee enhanced performance of PRS-based prediction models. Specifically, in scenarios with minimal misclassification rates and limited validation sizes, PRS models utilizing debiased regression coefficients demonstrated inferior predictive capabilities compared to models using error-prone phenotypes. Put differently, the effectiveness of the bias reduction method is contingent upon the misclassification rates of phenotypes and the size of the validation set employed during chart reviews. Notably, when dealing with datasets featuring higher misclassification rates, the advantages of utilizing this bias reduction method become more evident, requiring a smaller validation set to achieve better performance., Conclusion: This study highlights the importance of choosing an appropriate validation set size to balance between the efforts of chart review and the gain in PRS prediction accuracy. Consequently, our study establishes a valuable guidance for validation planning, across a diverse array of sensitivity and specificity combinations., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Prevalent Gene Regulatory Network (GRN) construction methods rely on generalized correlation analysis. However, in biological systems, regulation is essentially a causal relationship that cannot be adequately captured solely through correlation. Therefore, it is more reasonable to infer GRNs from a causal perspective. Existing causal discovery algorithms typically rely on Directed Acyclic Graphs (DAGs) to model causal relationships, but it often requires traversing the entire network, which result in computational demands skyrocketing as the number of nodes grows and make causal discovery algorithms only suitable for small networks with one or two hundred nodes or fewer. In this study, we propose the SLIVER (cauSaL dIscovery Via dimEnsionality Reduction) algorithm which integrates causal structural equation model and graph decomposition. SLIVER introduces a set of factor nodes, serving as abstractions of different functional modules to integrate the regulatory relationships between genes based on their respective functions or pathways, thus reducing the GRN to the product of two low-dimensional matrices. Subsequently, we employ the structural causal model (SCM) to learn the GRN within the gene node space, enforce the DAG constraint in the low-dimensional space, and guide each factor to aggregate various functions through cosine similarity. We evaluate the performance of the SLIVER algorithm on 12 real single cell transcriptomic datasets, and demonstrate it outperforms other 12 widely used methods both in GRN inference performance and computational resource usage. The analysis of the gene information integrated by factor nodes also demonstrate the biological explanation of factor nodes in GRNs. We apply it to scRNA-seq of Type 2 diabetes mellitus to capture the transcriptional regulatory structural changes of β cells under high insulin demand., Competing Interests: Declaration of competing interest None Declared., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Along the last decades, the genetic evaluation methodology has evolved, improving breeding value estimates. Many breeding programmes have historical phenotypic records and large number of generations, but to make use of them could result in more inconveniences than benefits. In this study, the prediction ability of genotyped young animals was assessed by simultaneously evaluating the removal of historical data, two pedigree deepness and two methodologies (traditional BLUP and single-step genomic BLUP or ssGBLUP), using milk yield records of 40 years of three Latxa dairy sheep populations. The linear regression method was used to compare predictions of young rams before and after progeny testing, with six cut-off points, by intervals of 4 years (from 1992 to 2012), and statistics of ratio of accuracies, bias, and dispersion were calculated. The prediction accuracy of selection candidates, when genomic information was included, was the highest in all Latxa populations (between 0.54 and 0.69 with full data set). Nevertheless, the deletion of historical phenotypic data resulted on moderate accuracy gain in the bigger data size populations (mean gain 2.5%), and the smaller population took advantage of a moderate data deletion (2.7% gain by removing data until 2004), reducing accuracy when more records were removed. The bias of validation individuals was lower when the breeding value was predicted based on genomic information (between 2.1 and 13.9), being lower when the biggest amount of data was deleted in the bigger data size populations (5.2% reduction), and the smaller population was benefited from data deletion between 1996 and 2008 (3.8% bias reduction). Meanwhile, the slope of estimated genetic trend was lower when less data were included, and an overestimation of the unknown parent group estimates was observed. The results indicated that ssGBLUP evaluations were outstanding, compared with traditional BLUP evaluations, while the depth of pedigree had a very small influence, and deletion of historical phenotypic data was beneficial. Thus, Latxa routine genetic evaluations would benefit from truncating phenotypic records between 2000 and 2004, the use of two pedigree generations and the implementation of ssGBLUP methodology., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

In beef cattle, dams play a crucial role in shaping the pre- and postnatal environment for the growth of their offspring. Acknowledging the substantial impact of maternal influence on the early development of calves, researchers utilize maternal animal models. These models take into account both maternal genetic and permanent environmental effects, operating under the assumption that these influences remain constant throughout the productive life of the cow. Nevertheless, it cannot be ruled out that these genetic and environmental effects may evolve throughout the lifespan of the cows. Therefore, this study aims to describe the changes in genetic and environmental maternal effects over the productive lifespan of cows. To accomplish this goal, we utilized random regression models, incorporating the age of the dam effect, maternal genetic effects, and environmental permanent effects using Legendre orthogonal polynomials. Additionally, the analytical model incorporated a covariate to adjust for the calf's age at recording, a two-level sex effect, a random herd-year-season effect, and an additive direct genetic effect associated linked to the calf. The dataset comprised information from dams aged between 2 and 16 years, resulting in a final database that comprised weight records of 58 332 calves from 21 673 dams. The average weight at 90 days was 135.0 ± 39.3 kg, and the mean age of the dam at calving was 7.03 ± 3.41 years. We evaluated models incorporating 2, 3, 4, 5, and 6 orthogonal polynomials alongside the standard maternal animal model. Afterward, we selected the model with five orthogonal polynomials based on the Akaike Information Criteria. The Restricted Maximum Likelihood estimates within this model indicated a direct heritability of around 0.50, and a maternal heritability ranging between 0.15 and 0.25, exhibiting a consistent increase between 4, 5 to 13 years. The genetic correlation estimates between direct and maternal genetic effects remained stable at approximately -0.55 across the lifespan of the cows. Furthermore, maternal genetic correlations between different ages of the dam decreased to around 0.7 for more distant age points. The maternal permanent correlations were notably lower, occasionally even reaching negative values, suggesting variability in environmental influence on maternal effects over the productive lifespan of the cow. Finally, the model enables the prediction of breeding values for the maternal genetic effects of the cow across its lifespan, providing opportunities for innovative selection strategies on the maternal side., (Copyright © 2024 The Author(s). Published by Elsevier B.V. All rights reserved.)

In the realm of animal phenotyping, manual measurements are frequently utilised. While machine-generated data show potential for enhancing high-throughput breeding, additional research and validation are imperative before incorporating them into genetic evaluation processes. This research presents a method for managing meat sheep and collecting data, utilising the Sheep Data Recorder system for data input and the Sheep Body Size Collector system for image capture. The study aimed to investigate the genetic parameter changes of growth traits in Ujumqin sheep by comparing machine-generated measurements with manual measurements. The dataset consisted of 552 data points from the offspring of 75 breeding rams and 399 breeding ewes. Six distinct random regression models were assessed to pinpoint the most suitable model for estimating genetic parameters linked to growth traits. These models were distinguished based on the inclusion or exclusion of maternal genetic effects, maternal permanent environmental effects, and covariance between maternal and direct genetic effects. Fixed factors such as individual age, individual sex, and ewe age were taken into account in the analysis. The genetic parameters for the yearling growth traits of Ujumqin sheep were calculated using ASReml software. The Akaike information criterion, the Bayesian information criterion, and fivefold cross-validation were employed to identify the optimal model. Research findings indicate that the most accurate models for manually measured data revealed heritability estimates of 0.12 ± 0.15 for BW, 0.05 ± 0.07 for body slanting length, 0.03 ± 0.07 for withers height, 0.15 ± 0.12 for hip height, 0.11 ± 0.11 for chest depth, 0.13 ± 0.13 for shoulder width, and 0.53 ± 0.15 for chest circumference. The optimal models for machine-predicted data showed heritability estimates of 0.1 ± 0.09 for body slanting length, 0.14 ± 0.12 for withers height, 0.55 ± 0.15 for hip height, 0.34 ± 0.15 for chest depth, 0.26 ± 0.15 for shoulder width, and 0.47 ± 0.16 for chest circumference. In manually measured data, genetic correlations ranged from 0.35 to 0.99, while phenotypic correlations ranged from 0.07 to 0.90. In machine data, genetic correlations ranged from -0.05 to 0.99, while phenotypic correlations ranged from 0.03 to 0.84. The results suggest that machine-based estimations may lead to an overestimation of heritability, but this discrepancy does not impact the selection of breeding models., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

Tree shape statistics based on peripheral structures have been utilized to study evolutionary mechanisms and inference methods. Partially motivated by a recent study by Pouryahya and Sankoff on modeling the accumulation of subgenomes in the evolution of polyploids, we present the distribution of subtree patterns with four or fewer leaves for the unrooted Proportional to Distinguishable Arrangements (PDA) model. We derive a recursive formula for computing the joint distributions, as well as a Strong Law of Large Numbers and a Central Limit Theorem for the joint distributions. This enables us to confirm several conjectures proposed by Pouryahya and Sankoff, as well as provide some theoretical insights into their observations. Based on their empirical datasets, we demonstrate that the statistical test based on the joint distribution could be more sensitive than those based on one individual subtree pattern to detect the existence of evolutionary forces such as whole genome duplication., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

The traditional genetic evaluation methods generally consider additive genetic effects only and often ignore non-additive (dominance and epistasis) effects that may have contributed to genetic variation of complex traits of livestock species. The available dense single nucleotide polymorphisms (SNPs) panels offer to investigate the potential benefits of including non-additive genetic effects in the genomic evaluation models. Data from 16 971 genotyped (Illumina Bovine 50 K SNP chip) Korean Hanwoo cattle were used to estimate genetic variance components and prediction accuracy of genomic breeding values (GEBVs) for four carcass and meat quality traits: carcass weight (CWT), eye muscle area (EMA), back fat thickness (BFT) and marbling score (MS). Five different genetic models were evaluated through including additive, dominance and epistatic interactions (additive by additive, A × A; additive by dominance, A × D and dominance by dominance, D × D) successively in the models. The estimates of additive genetic variances and narrow sense heritabilities (h a 2 ) were found similar across the evaluated models and traits except when additive interaction (A × A) was included. The dominance variance estimates relative to phenotypic variance ranged from 1.7-3.4% for CWT and MS traits, whereas, they were close to zero for EMA and BFT traits. The magnitude of A × A epistatic heritability (h aa 2 ) ranged between 14.8 and 27.7% in all traits. However, heritability estimates for A × D and D × D epistatic interactions (h ad 2 and h dd 2 ) were quite low compared to h aa 2 and were contributed only 0.0-9.7% of the total phenotypic variation. In general, broad sense heritability (h G 2 ) estimates were almost twice (ranging between 0.54 and 0.68) the h a 2 for all of the investigated traits. The inclusion of dominance effects did not improve the prediction accuracy of GEBV but improved 2.0-3.0% when epistatic effects were included in the model. More importantly, rank correlation revealed that partitioning of variance components considering dominance and epistatic effects in the model would enable to re-rank of top animals with better prediction of GEBV. The present result suggests that dominance and epistatic effects could be included in the genomic evaluation model for better estimates of variance components and more accurate prediction of GEBV for carcass and meat quality traits in Korean Hanwoo cattle., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

Studying gene regulatory networks associated with cancer provides valuable insights for therapeutic purposes, given that cancer is fundamentally a genetic disease. However, as the number of genes in the system increases, the complexity arising from the interconnections between network components grows exponentially. In this study, using Boolean logic to adjust the existing relationships between network components has facilitated simplifying the modeling process, enabling the generation of attractors that represent cell phenotypes based on breast cancer RNA-seq data. A key therapeutic objective is to guide cells, through targeted interventions, to transition from the current cancer attractor to a physiologically distinct attractor unrelated to cancer. To achieve this, we developed a computational method that identifies network nodes whose inhibition can facilitate the desired transition from one tumor attractor to another associated with apoptosis, leveraging transcriptomic data from cell lines. To validate the model, we utilized previously published in vitro experiments where the downregulation of specific proteins resulted in cell growth arrest and death of a breast cancer cell line. The method proposed in this manuscript combines diverse data sources, conducts structural network analysis, and incorporates relevant biological knowledge on apoptosis in cancer cells. This comprehensive approach aims to identify potential targets of significance for personalized medicine., Competing Interests: Declaration of Competing Interest None Declared., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

[S U M M A R Y] Many miRNA-disease association prediction models incorporate Gaussian interaction profile kernel similarity (GIPS). However, the GIPS fails to consider the specificity of the miRNA-disease association matrix, where matrix elements with a value of 0 represent miRNA and disease relationships that have not been discovered yet. To address this issue and better account for the impact of known and unknown miRNA-disease associations on similarity, we propose a method called vector projection similarity-based method for miRNA-disease association prediction (VPSMDA). In VPSMDA, we introduce three projection rules and combined with logistic functions for the miRNA-disease association matrix and propose a vector projection similarity measure for miRNAs and diseases. By integrating the vector projection similarity matrix with the original one, we obtain the improved miRNA and disease similarity matrix. Additionally, we construct a weight matrix using different numbers of neighbors to reduce the noise in the similarity matrix. In performance evaluation, both LOOCV and 5-fold CV experiments demonstrate that VPSMDA outperforms seven other state-of-the-art methods in AUC. Furthermore, in a case study, VPSMDA successfully predicted 10, 9, and 10 out of the top 10 associations for three important human diseases, respectively, and these predictions were confirmed by recent biomedical resources., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Inc. All rights reserved.)

Nowadays, several countries are developing or adopting genomic selection in the dairy goat sector. The most used method to estimate breeding values is Single-Step Genomic Best Linear Unbiased Prediction (ssGBLUP) which offers several advantages in terms of computational process and accuracy of the estimated breeding values (EBVs). Saanen and Alpine are the predominant dairy goat breeds in Italy, and both have similar breeding programs where EBVs for productive traits are currently calculated using BLUP. This work describes the implementation of genomic selection for these two breeds in Italy, aligning with the selection practices already carried out in the international landscape. The available dataset included 3 611 genotyped animals, 11 470 lactation records, five traits (milk, protein and fat yields, and fat and protein percentages), and three-generation pedigrees. EBVs were estimated using BLUP, GBLUP, and ssGBLUP both with single and multiple trait approaches. The methods were compared in terms of correlation between EBVs and genetic trends. Results were also validated with the linear regression method excluding part of the phenotypic data. In both breeds, EBVs and GEBVs were strongly correlated and the trend of each trait was similar comparing the three methods. The average increase in accuracy across traits and methods amounted to +13 and +10% from BLUP to ssGBLUP for Alpine and Saanen breeds, respectively. Results indicated higher prediction accuracy and correlation for GBLUP and ssGBLUP compared to BLUP, implying that the use of genotypes increases the accuracy of EBVs, particularly in the absence of phenotypic data. Therefore, ssGBLUP is likely to be the most effective method to enhance genetic gain in Italian Saanen and Alpine goats., (Copyright © 2024 The Author(s). Published by Elsevier B.V. All rights reserved.)

Boolean modeling is a mathematical modeling framework used for defining and studying gene-regulatory networks (GRNs). It serves as a means to develop mechanistic models, offering insights into the trajectories and dynamic properties of GRNs. In this review, I delve into seminal papers published in the Journal of Theoretical Biology that have spearheaded this field. Additionally, I explore the application of these modeling methods in the current era of data-intensive science., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

The size of the reference group is among the most critical determinants of genomic estimated breeding values (GEBVs) accuracy. However, small- and medium-sized pig farms often need help accumulating adequate reference data, posing significant challenges to breeding programs. To solve this problem, exploring the potential benefits of combining reference groups of different sizes is necessary to improve GEBV accuracy. The primary objective of this investigation was to assess a more effective statistical model for combined multi-populations and its potential to enhance the accuracy of GEBVs for small and medium populations. Three populations were simulated using the QMSim software, each consisting of different sizes (300, 600, and 1 500, respectively). To assess the impact of heritability on the accuracy of GEBVs, four different levels of heritability (0.05, 0.15, 0.35, and 0.5) were simulated. Simultaneously, to investigate the impact of kinship on multi-populations, the study created four distinct scenarios for the three sizes of populations. These scenarios included: (1) the three groups are all independent, (2) the large group and the small group with a familial connection (n = 1 800), a middle group (n = 600) acting independently with no kinship, (3) the large group with a familial connection to the middle group (n = 2 100) but no connection to the small group (n = 300), and (4) the small group with a familial connection to the middle group (n = 900), while the large group (n = 1 500) acted independently with no kinship. This study evaluates and compares the accuracy of predicting breeding values using four different methods, including genomic best linear unbiased prediction (GBLUP), single-stepGBLUP (ssGBLUP), and two Bayesian models (Bayes A and Bayes B), with varying sizes of reference groups. In each scenario, three different prediction strategies were compared: (1) Merging all three different sizes of populations for predicting, (2) predicting each independent population separately, and (3) the other two populations predict the population. Our findings reveal that combining populations enhances the Bayesian models, with Bayes B yielding the highest accuracy. In independent populations, the best linear unbiased prediction (BLUP) models demonstrated the highest accuracy. However, in cases where populations were related and the heritability was high, the Bayes B model exhibited the highest overall accuracy (slightly higher than BLUP models) in the independent population. Our results underscore the importance of considering population combinations when using genetic models to predict breeding values, particularly for pig farmers with limited resources., (Copyright © 2023 The Author(s). Published by Elsevier B.V. All rights reserved.)

Australian beef cattle experience variable conditions, which may give rise to genotype-by-environment interactions depending on the genotypes' macro- and/or micro-genetic environmental sensitivity (GES). Macro-GES gives rise to genotype-by-environment interactions across definable and shared environments, while micro-GES causes heritable variation of phenotypes, e.g., the performance of progeny from one sire may be more variable than other sires. Yearling weight (YW) is a key trait in Australian Angus cattle that may be impacted by both macro- and micro-GES. Current models for genetic evaluation of YW attempt to account for macro-GES by fitting sire-by-herd interactions (S × H). Variation in micro-GES had not yet been estimated for YW in Australian Angus. The aim of this study was to estimate genetic variation due to macro- and micro-GES in YW of Australian Angus cattle. A reaction norm with contemporary group effects as the environmental covariate was fitted either as an alternative to or in combination with a random S × H effect to account for macro-GES. Double hierarchical generalised linear models (DHGLM), fitted as sire models, were used to estimate the genetic variance of the dispersion as a measure of micro-GES. Variation due to both macro- and micro-GES were found in YW. The variance of the slope of the reaction norm was 0.02-0.03 (SEs 0.00), while the S × H variance accounted for 7% of the phenotypic variance in all models. Results showed that both a random S × H effect and a reaction norm should be included to account for both macro-GES and the additional variation captured by an S × H effect. The heritability of the dispersion on the measurement scale ranged from 0.06 to 0.10 (SEs 0.00) depending on which model was used. It should therefore be possible to alter both macro- and micro-GES of YW in Australian Angus through selection. However, care should be taken to ensure an appropriate data structure when including sire-by-herd interactions in the mean part of a DHGLM; otherwise, it can cause biased estimates of micro-GES., (Copyright © 2023. Published by Elsevier B.V.)

Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

The purebred-crossbred genetic correlation (r pc ) is a key parameter to determine whether the optimal selection of purebred animals to improve crossbred performance should rely on crossbred phenotypes, purebred phenotypes, or both. We reviewed published estimates of the r pc in poultry. In total, 19 studies were included, of which four were on broilers and 15 on laying hens, with 150 r pc estimates for nine different trait categories. Average reported r pc estimates were highest for egg weight, egg quality and egg colour (0.74-0.82), intermediate for BW, maturity and mortality (0.61-0.70) and egg number (0.58), and low for resilience (0.40) and body conformation (0.14). Most studies were based on measuring purebred and crossbred phenotypes in the same environment and thus did not capture the contribution of genotype by environment interactions to the r pc , suggesting that the presented average estimates may be higher than values that apply in practice. Nearly all studies were based on two-way crossbred animals. We hypothesised that r pc values for a two-way cross are good proxies for r pc of a four-way cross. Only eight out of 19 studies were published in the last 25 years, and only two of those used genomic data. We expect that more studies using genomic data may be published in the coming years, as the required data may be generated when implementing genomic selection for crossbred performance, which will lead to more accurate r pc estimates. Future studies that aim to estimate r pc are encouraged to capture the genotype by environment interaction component by housing purebred and crossbred animals differently as is done in practice. Moreover, there is a need for further studies that enable to explicitly estimate the magnitude of genotype by environment versus genotype by genotype interactions for multiple trait categories. Further, studies are advised to report: the specific housing conditions of the animals, any differences between measurements of purebred versus crossbred performance, and the heritabilities of purebred and crossbred performance., (Copyright © 2023 The Author(s). Published by Elsevier B.V. All rights reserved.)

Genomic prediction (GP) has greatly advanced animal and plant breeding over the past two decades. GP in joint populations is a feasible method to improve the accuracy of genomic estimated breeding values in small populations. However, there is still a need to understand the factors that influence GP in joint populations. This study used simulated data and real data from Duroc pig populations to examine the impact of linkage disequilibrium (LD), causal variants effect sizes (CVESs), and minor allele frequencies (MAF) of SNPs on the accuracy of genomic prediction in joint populations. Three prediction methods were used: genomic best linear unbiased prediction (GBLUP), single-step GBLUP and multi-trait GBLUP. Results from the simulated datasets showed that the accuracies of GP in joint populations were always higher than those in a single population when only LD inconsistencies existed. However, single-step GBLUP accuracy in joint populations decreased as the correlation of MAF between populations decreased, while the accuracy of GBLUP is consistently higher in joint populations than in a single population. As the correlation of CVES between populations decreased, the accuracy of both GBLUP and single-step GBLUP in joint populations declined. Analysis of real Duroc populations showed low genetic correlation, similar to the simulated relationship between the most distant populations. In most cases in Duroc populations, GP have higher accuracies in joint populations than in individual population. In conclusion, the consistency of CVES plays a more important role in multi-population GP. The genetic relatedness of the Duroc populations is so weak that the prediction accuracy of GP in joint populations is reduced in some traits. Multi-trait GBLUP is a competitive method for the joint breeding evaluation., (Copyright © 2023 The Author(s). Published by Elsevier B.V. All rights reserved.)

Intrinsic hybrid incompatibility is one of the important isolating barriers between species. In organisms with sex chromosomes, intrinsic hybrid incompatibility often follows two rules: Haldane's rule and large-X effects. One explanation for these two rules is that sex chromosomes are hotspots for meiotic drivers that can cause intrinsic hybrid incompatibility between geographically isolated populations. Although this hypothesis seems plausible and several empirical data are consistent with it, we are still unsure whether such mechanisms occur in nature, particularly with respect to speciation with gene flow. Here, we review empirical studies that have investigated the roles of meiotic drive in sex-chromosome evolution and speciation and propose future studies necessary for testing this hypothesis., Competing Interests: Declaration of Competing Interest The authors declare no conflict of interest., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

The interphase genome of vertebrates contains roughly 100 000 dynamic loops formed by cohesins. These loops are thought to play important roles in many functions, but their exact contribution in each case remains hotly disputed. The key challenge in studying these loops is the lack of a single experimental technique that could reliably and comprehensively visualize their locations and dynamics. Yet, we can infer them using theoretical models that integrate complementary experimental observations. Modeling proved instrumental in showing that cohesins form loops via extrusion. The loop extrusion model made numerous successful qualitative and quantitative predictions and inspired many experiments. However, it also demonstrated limited accuracy in predicting contact maps. Recent research suggests that the original model did not fully account for the intricate details of the mechanism of loop extrusion and its complex regulation. Here, we review the progress in visualizing extrusion and characterizing the cohesin cofactors. These discoveries can be summarized as 'rules' of cohesin movement along chromosomes and incorporated into the next generation of models. Such improved models will enable more accurate inferences of positions and dynamics of cohesin loops and generate better predictions for designing experiments., Competing Interests: Declaration of Competing Interest The authors have no conflicts of interest to declare., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Understanding the role of natural selection in driving evolutionary change requires accurate estimates of the strength of selection acting at the genetic level in the wild. This is challenging to achieve but may be easier in the case of populations in migration-selection balance. When two populations are at equilibrium under migration-selection balance, there exist loci whose alleles are selected different ways in the two populations. Such loci can be identified from genome sequencing by their high values of F ST . This raises the question of what is the strength of selection on locally-adaptive alleles. To answer this question we analyse a 1-locus 2-allele model of a population distributed between two niches. We show by simulation of selected cases that the outputs from finite-population models are essentially the same as those from deterministic infinite-population models. We then derive theory for the infinite-population model showing the dependence of selection coefficients on equilibrium allele frequencies, migration rates, dominance and relative population sizes in the two niches. An Excel spreadsheet is provided for the calculation of selection coefficients and their approximate standard errors from observed values of population parameters. We illustrate our results with a worked example, with graphs showing the dependence of selection coefficients on equilibrium allele frequencies, and graphs showing how F ST depends on the selection coefficients acting on the alleles at a locus. Given the extent of recent progress in ecological genomics, we hope our methods may help those studying migration-selection balance to quantify the advantages conferred by adaptive genes., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Genomic selection using single nucleotide polymorphism (SNP) markers is now intensively investigated in breeding and has been widely utilized for genetic improvement. Currently, several studies have used haplotype (consisting of multiallelic SNPs) for genomic prediction and revealed its performance advantage. In this study, we comprehensively evaluated the performance of haplotype models for genomic prediction in 15 traits, including 6 growth, 5 carcass, and 4 feeding traits in a Chinese yellow-feathered chicken population. We adopted 3 methods to define haplotypes from high-density SNP panels, and our strategy included combining Kyoto Encyclopedia of Genes and Genomes pathway information and considering linkage disequilibrium (LD) information. Our results showed an increase in prediction accuracy due to haplotypes ranging from -0.04∼27.16% in all traits, where the significant improvements were found in 12 traits. The estimates of haplotype epistasis heritability were strongly correlated with the accuracy increase by haplotype models. In addition, incorporating genomic annotation information could further increase the accuracy of the haplotype model, where the further increase in accuracy is significantly relative to the increase of relative haplotype epistasis heritability. The genomic prediction using LD information for constructing haplotypes has the best prediction performance among the 4 traits. These results uncovered that haplotype methods were beneficial for genomic prediction, and the accuracy could be further increased by incorporating genomic annotation information. Moreover, using LD information would potentially improve the performance of genomic prediction., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Random regression modelling has been used across multiple animal species to model longitudinal data. The random regression model for growth accounts for the genetic correlation between measures of the same trait over time and the wide environmental variability in growth, but this requires adequate weight records across the age range. However, contemporary management practices in sheep in the United Kingdom generally focus on growing lambs and neglect mature weight recordings. This study examined modelling strategies for growth data in Suffolk and Charollais sheep, provided by the Agriculture and Horticulture Development Board, with polynomial random regression modelling with many early life weight recordings but limited weight recordings in mature animals. Two methods were employed to model the data. In Method A, missing mature weight records were predicted for those animals that did not have a recorded mature weight. The animals were sorted into groups based on the identity of their sires and the year in which the animal was born. Mature weight values were predicted within each group with a multiple regression model. The dataset, including predicted values, was analysed with random regression models using polynomials and simple linear regression for animal and permanent environmental (PE) effects. In Method B, the dataset with missing mature weight records was analysed using a random linear regression animal model with random animal and PE effects. Due to problems of convergence because the parameters were close to the boundary space, fixing the correlation between the intercept and slope of the Legendre polynomial at different levels was investigated. The heritability values resulting from the model with a fixed correlation between intercept and slope parameters at 0.5 for the Suffolk dataset resulted in heritability values ranging from 0.2 to 0.5 from 1 to 619 days of age. Corresponding estimates for the Charollais dataset ranged from 0.18 to 0.49 from 1 to 640 days of age. For the Suffolk data, the genetic correlations ranged from 1.00 to 0.08 between weight at day 1 to weight at day 619, while for the Charollais, the correlations ranged from 1.00 to 0.05 from 1 to 640 days of age. Validation procedures were undertaken using a multitrait approach to examine the estimated breeding values when the correlation between the intercept and slope are fixed at different levels. The results indicated that fixing the correlation at 0.5 gave the most appropriate estimates for the Suffolk and Charollais datasets., (Copyright © 2023 The Authors. Published by Elsevier B.V. All rights reserved.)

Neutral theory predicts that the genetic diversity within a population is proportional to the effective population size. In contrast, the observed genetic diversity for various species vary in a narrow range for several orders of magnitude change in the population sizes (Lewontin's paradox). The selective sweeps and background selection, reduce the genetic variation at the linked neutral sites and have been studied considering the environment to be selectively constant. However, in a natural population, the selective environment varies with time. Here, we investigate the impact of selective sweeps and background selection on neutral genetic diversity when the selection coefficient changes periodically over time. The reduction in genetic variation due to selective sweeps is known to depend on the conditional fixation time. Here, we find that the effect of changing environment on conditional mean fixation time is most substantial for the randomly mating population than the inbreeding population with arbitrary inbreeding coefficient. We also study the effect of background selection on neutral sites when the selection coefficient of linked deleterious mutation change periodically in time. In the slowly changing environment, we find that neutral heterozygosity is significantly different, and the site frequency spectrum has a different shape than that in the static environment., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Nowadays, in some populations, the number of genotyped animals is too large to obtain the inverse of the genomic relationship matrix. The algorithm for proven and young animals (APY) can be used to overcome this problem. In the present work, different strategies for defining core animals in APY were tested using either simulated or real data. In particular, core definitions based on random choice or on the contribution to the genomic relationship matrix (G CONTR ) calculated using Principal Component Analysis were tested. Core sizes able to explain 90, 95, 98, and 99% of the total variance of the genomic relationship matrix (G) were used. Analyzed phenotypes were three simulated traits for 3 000 individuals, and milkability records for 136 406 Italian Simmental cows. The number of genotypes was 4 100 for the simulated dataset, and 11 636 for the Simmental data, respectively. The G CONTR values in Simmental dataset were moderately correlated with the analyzed phenotype, and they showed a decreasing trend according to the year of birth of genotyped animals. The accuracy increased as the size of the core increased in both datasets. The inclusion in the core of animals with largest G CONTR values led to the lowest accuracies (0.50 and 0.71 for the simulated and Simmental datasets, respectively; average across traits and core sizes). On the contrary, the selection of animals with the lowest rank according to their contribution to the G provided slightly higher accuracies, especially in the simulated dataset (0.68 for the simulated dataset, and 0.76 for the Simmental data; average across traits and core sizes). In real data, particularly for larger sizes of core animals, the criteria of choice appear less important, confirming the results of earlier studies. Anyway, the inclusion in the core of animals with the lowest values of G CONTR led to increases in accuracy. These are preliminary results based on a small sample size that need to be confirmed on a larger number of genotypes., (Copyright © 2023 The Author(s). Published by Elsevier B.V. All rights reserved.)

The moral rule "Risk your life to save your family members" is, at the same time, a biological phenomenon. The prominent population geneticist, J.B.S. Haldane told his friends that he would risk his life to save two drowning brothers, but not one - so the story goes. In biological terms, Haldane's arithmetic claims that sib altruism is evolutionarily rational, whenever by "self-sacrifice" an altruistic gene "rescues", on average, more than one copy of itself in its lineage. Here, we derive conditions for evolutionary stability of sib altruism, using population genetic models for three mating systems (monogamy, promiscuity and polygyny) with linear and non-linear group effect on the siblings' survival rate. We show that for all considered selection situations, the condition of evolutionary stability is equivalent to Haldane's arithmetic. The condition for evolutionary stability is formulated in terms of genetic relatedness and the group effect on the survival probability, similarly to the classical Hamilton's rule. We can set up a "scale of mating systems", since in pairwise interactions the chance of evolutionary stability of sib altruism decreases in this order: monogamy, polygyny and promiscuity. Practice of marrying and siblings' solidarity are moral rules in a secular world and in various religious traditions. These moral rules are not evolutionarily independent, in the sense that the subsistence of sib altruism is more likely in a monogamous population., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

Pedigree likelihood ratios are widely used in forensic genetics. Linkage of genetic markers affects the pedigree likelihood ratio. In this study, a series of novel formulae was established based on simplification due to founder and founder couple symmetry of the Lander-Green algorithm for dealing with linkage. These new formulae are based on the inheritance vectors and differ from existing formulae that are based on the number of alleles shared identical-by-descent. In contrast to previous formulae, the formulae of this study expand the calculation to more general relationships, which enables specifying the number of typed individuals in pedigrees, and the new formulae are no longer limited to calculation of pairs of linked markers. Alternatively, pedigree likelihood values incorporating multiple linked markers could be obtained by iterating the vector transition probability between pairs of linkage markers and the probability of genotypes of typed individuals given inheritance vectors; the vector transition probability formulae are provided in a table form. To demonstrate utilization of the proposed formulae system, a series of cases were designed based on the relationships to validate GeneVisa software (www.genevisa.net) with frequency data of three linked markers, SE33-D6S1043-D6S474, on chromosome 6. The calculation results were verified with the FamLink software. The results demonstrated that the proposed formulae can correctly obtain the pedigree likelihood ratio, demonstrating the application potential of these formulae to verify effectiveness of DNA testing software., (Copyright © 2022 Elsevier B.V. All rights reserved.)

The fluctuation assay remains an important tool for analyzing the levels of mutagenesis in microbial populations. The mutant counts originating from some average number of mutations are usually assumed to obey the Luria-Delbrück distribution. While several tools for estimating mutation rates are available, they sometimes lack accuracy or versatility under non-standard conditions. In this work, extensions to the Luria-Delbrück protocol to account for phenotypic lag and cellular death with either perfect or partial plating were developed. Hence, the novel MLE MUtation Rate calculator, or mlemur, is the first tool that provides a user-friendly graphical interface allowing the researchers to model their data with consideration for partial plating, differential growth of mutants and non-mutants, phenotypic lag, cellular death, variability of the final number of cells, post-exponential-phase mutations, and the size of the inoculum. Additionally, mlemur allows the users to incorporate most of these special conditions at the same time to obtain highly accurate estimates of mutation rates and P values, confidence intervals for an arbitrary function of data (such as fold), and perform power analysis and sample size determination for the likelihood ratio test. The accuracy of point and interval estimates produced by mlemur against historical and simulated fluctuation experiments are assessed. Both mlemur and the analyses in this work might be of great help when evaluating fluctuation experiments and increase the awareness of the limitations of the widely-used Lea-Coulson formulation of the Luria-Delbrück distribution in the more realistic biological contexts., Competing Interests: Declaration of Competing Interest The author declares no conflict of interest., (Copyright © 2023 The Authors. Published by Elsevier B.V. All rights reserved.)

The origin of the genetic code is an abiding mystery in biology. Hints of a 'code within the codons' suggest biophysical interactions, but these patterns have resisted interpretation. Here, we present a new framework, grounded in the autotrophic growth of protocells from CO 2 and H 2 . Recent work suggests that the universal core of metabolism recapitulates a thermodynamically favoured protometabolism right up to nucleotide synthesis. Considering the genetic code in relation to an extended protometabolism allows us to predict most codon assignments. We show that the first letter of the codon corresponds to the distance from CO 2 fixation, with amino acids encoded by the purines (G followed by A) being closest to CO 2 fixation. These associations suggest a purine-rich early metabolism with a restricted pool of amino acids. The second position of the anticodon corresponds to the hydrophobicity of the amino acid encoded. We combine multiple measures of hydrophobicity to show that this correlation holds strongly for early amino acids but is weaker for later species. Finally, we demonstrate that redundancy at the third position is not randomly distributed around the code: non-redundant amino acids can be assigned based on size, specifically length. We attribute this to additional stereochemical interactions at the anticodon. These rules imply an iterative expansion of the genetic code over time with codon assignments depending on both distance from CO 2 and biophysical interactions between nucleotide sequences and amino acids. In this way the earliest RNA polymers could produce non-random peptide sequences with selectable functions in autotrophic protocells., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Authors. Published by Elsevier B.V. All rights reserved.)

The Wright-Fisher binomial model of allele frequency change is often approximated by a scaling limit in which selection, mutation and drift all decrease at the same 1/N rate. This construction restricts the applicability of the resulting 'Wright-Fisher diffusion equation' to the weak selection, weak mutation regime of evolution. We argue that diffusion approximations of the Wright-Fisher model can be used more generally, for instance in cases where genetic drift is much weaker than selection. One important example of this regime is Muller's ratchet phenomenon, whereby deleterious mutations slowly but irreversibly accumulate through rare stochastic fluctuations. Using a modified diffusion equation we derive improved analytical estimates for the mean click time of the ratchet., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

The coalescent model represents how individuals sampled from a population may have originated from a last common ancestor. The bounded coalescent model is obtained by conditioning the coalescent model such that the last common ancestor must have existed after a certain date. This conditioned model arises in a variety of applications, such as speciation, horizontal gene transfer or transmission analysis, and yet the bounded coalescent model has not been previously analysed in detail. Here we describe a new algorithm to simulate from this model directly, without resorting to rejection sampling. We show that this direct simulation algorithm is more computationally efficient than the rejection sampling approach. We also show how to calculate the probability of the last common ancestor occurring after a given date, which is required to compute the probability density of realisations under the bounded coalescent model. Our results are applicable in both the isochronous (when all samples have the same date) and heterochronous (where samples can have different dates) settings. We explore the effect of setting a bound on the date of the last common ancestor, and show that it affects a number of properties of the resulting phylogenies. All our methods are implemented in a new R package called BoundedCoalescent which is freely available online., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Piezo1 and Piezo2 are recently discovered mechanosensory ion channels. Piezo channels transduce mechanical stimulation into cellular signaling in a variety of tissues and organ systems. The functional roles of Piezo1 and Piezo2 have been revealed in both developmental and physiological scenarios by using mouse genetic models. Mechanotransduction by Piezo1 channels regulates osteoblast/osteocyte activity and, thus, strengthens the skeleton enabling it to adapt to a wide range of mechanical loadings. Deletion of the Piezo1 gene in the developing skeleton causes bone malformations that lead to spontaneous bone fractures, while inactivity of Piezo1 in adulthood results in osteoporosis. Furthermore, Piezo2 channels in sensory neurons might provide another route of skeletal regulation. Piezo channels also regulate the proliferation and differentiation of various types of stem cells. PIEZO1 and PIEZO2 mutations and channel malfunctions have been implicated in an increasing number of human diseases, and PIEZO channels are currently emerging as potential targets for disease treatment. This review summarizes the important findings of Piezo channels for skeletal development and homeostasis using the mouse genetic model system., (Copyright © 2022 International Society of Differentiation. Published by Elsevier B.V. All rights reserved.)

The first step in statistical inference of the evolutionary histories of species is developing a probability model that describes the mutation process as accurately and realistically as possible. A major complication of this inference is that different loci on the genome can have histories that diverge from the common species history and each other. The multispecies coalescent process is commonly used to model one source of this divergence, incomplete lineage sorting, or ILS. In Chifman and Kubatko (2015), the authors computed the site pattern probabilities for four taxa under a full probability model based on the Jukes-Cantor substitution model when the molecular clock holds. This paper generalizes that work to a relaxed clock model, allowing for mutation rates to differ among species. This will enable better phylogentic inference in cases where the molecular clock does not hold., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022. Published by Elsevier Ltd.)