Tom Druet

The RZooRoH R package was implemented to characterize individual inbreeding levels. It identifies DNA segments inherited twice from a common ancestor through different paths, which are known as homozygous‐by‐descent (HBD) segments. The package accepts different data formats and provides multiple outputs: HBD segments, inbreeding rates and genome‐wide and locus‐specific HBD probabilities. In addition, it partitions HBD levels into multiple HBD classes. The length distribution varies between these classes, which therefore correspond to distinct groups of ancestors that can be traced back to different generations in the past. This provides information about mating structure and recent demographic history. The computational performance of the package has been substantially improved, enabling, for example, computing times to be reduced when working with whole‐genome sequence data and more HBD classes to be fitted. It is now possible to fit one class per past generation, which facilitates interpretation of the results. Since we have previously demonstrated that the ZooRoH model can be used to characterize identity‐by‐descent (IBD) between haploid individuals or phased haplotypes, this option has been included in the new package version. Estimating kinship by characterizing IBD levels between the four possible pairs of haplotypes from two individuals is another feature we added to the package. Finally, new options allow models to be refined, for instance by defining HBD classes as intervals or constant inbreeding rates for neighbouring classes. Overall, the new version of the package offers improved computational efficiency and interpretability when characterizing inbreeding, IBD and relatedness levels.

Increasing agricultural production while reducing the reliance on synthetic inputs such as antibiotics is an important challenge. For cattle breeding, this implies better understanding the genetics underlying meat production and the immune response. Here, we use systems immunology to investigate the genetic and environmental drivers of immune variation in Belgian White Blue male cattle, a breed historically bred for meat production. While seasonality and other non-genetic factors account for much of the immune variation observed, genome-wide association studies identify loci with major effects on specific immunophenotypes. Genetics also emerges as the primary driver of cytokine production. Finally, we develop a predictive model linking genetic data to cytokine responses. Our findings support the selection of cattle with improved immunity and advance our understanding of mammalian immune variation.

The RZooRoH R package was implemented to characterize individual inbreeding levels. It identifies DNA segments inherited twice from a common ancestor through different paths, which are known as homozygous-by-descent (HBD) segments. The package accepts different data formats and provides multiple outputs: HBD segments, inbreeding rates and genome-wide and locus-specific HBD probabilities. In addition, it partitions HBD levels into multiple HBD classes. The length distribution varies between these classes, which therefore correspond to distinct groups of ancestors that can be traced back to different generations in the past. This provides information about mating structure and recent demographic history. The computational performance of the package has been substantially improved, enabling, for example, computing times to be reduced when working with whole-genome sequence data and more HBD classes to be fitted. It is now possible to fit one class per past generation, which facilitates interpretation of the results. Since we have previously demonstrated that the ZooRoH model can be used to characterize identity-by-descent (IBD) between haploid individuals or phased haplotypes, this option has been included in the new package version. Estimating kinship by characterizing IBD levels between the four possible pairs of haplotypes from two individuals is another feature we added to the package. Finally, new options allow models to be refined, for instance by defining HBD classes as intervals or constant inbreeding rates for neighboring classes. Overall, the new version of the package offers improved computational efficiency and interpretability when characterizing inbreeding, IBD and relatedness levels.

In populations of small effective size (N_e), such as those in conservation programmes, companion animals or livestock species, inbreeding control is essential. Homozygosity‐by‐descent (HBD) segments provide relevant information in that context, as they allow accurate estimation of the inbreeding coefficient, provide locus‐specific information and their length is informative about the “age” of inbreeding. Our objective was to evaluate tools for predicting HBD in future offspring based on parental genotypes, a problem equivalent to identifying segments identical‐by‐descent (IBD) among the four parental chromosomes. In total, we reviewed and evaluated 16 approaches using simulated and real data from populations with small N_e. The methods included model‐based approaches as well as more computationally efficient rule‐based approaches. The accuracy of the methods was then evaluated, including with low‐density marker panels, genotyping‐by‐sequencing data and small groups of individuals, typical features of such populations. Two model‐based approaches performed consistently well, while some rule‐based approaches proved accurate for genome‐wide predictions. The model‐based approaches were particularly efficient when genomic information was sparse or degraded. Methods using phased data proved to be more accurate, while some approaches relying on unphased genotype data were sensitive to the assumed allele frequencies. In some settings, pedigree‐based predictions ranked high for recent inbreeding levels. Finally, we showed that our evaluation is also informative about the accuracy of the methods for estimating relatedness and identifying IBD segments between pairs of present‐day individuals. This study shows that future inbreeding can be accurately predicted, including at specific loci, but not all methods perform equally well.

The feral cattle of the subantarctic island of Amsterdam provide an outstanding case study of a large mammalian population that was established by a handful of founders and thrived within a few generations in a seemingly inhospitable environment. Here, we investigated the genetic history and composition of this population using genotyping and sequencing data. Our inference showed an intense but brief founding bottleneck around the late 19th century and revealed contributions from European taurine and Indian Ocean Zebu in the founder ancestry. Comparative analysis of whole-genome sequences further revealed a moderate reduction in genetic diversity despite high levels of inbreeding. The brief and intense bottleneck was associated with high levels of drift, a flattening of the site frequency spectrum and a slight relaxation of purifying selection on mildly deleterious variants. Unlike some populations that have experienced prolonged reductions in effective population size, we did not observe any significant purging of highly deleterious variants. Interestingly, the population’s success in the harsh environment can be attributed to preadaptation from their European taurine ancestry, suggesting no strong bioclimatic challenge, and also contradicting evidence for insular dwarfism. Genome scan for footprints of selection uncovered a majority of candidate genes related to nervous system function, likely reflecting rapid feralization driven by behavioral changes and complex social restructuring. The Amsterdam Island cattle offers valuable insights into rapid population establishment, feralization, and genetic adaptation in challenging environments. It also sheds light on the unique genetic legacies of feral populations, raising ethical questions according to conservation efforts.

In populations of small effective size (N _e ), such as those in conservation programs, companion animals or livestock species, management of diversity and inbreeding is essential. Homozygosity-by-descent (HBD) segments provide relevant information in that context, as they allow efficient estimation of the inbreeding coefficient, provide locus-specific information and their length is informative about the “age” of inbreeding. Therefore, our objective was to evaluate tools for predicting HBD in future offspring based on parental genotypes, a problem equivalent to identifying segments identical-by-descent (IBD) among the four parental chromosomes. In total, we reviewed and evaluated 16 approaches using simulated and real data with small N _e . The methods included model-based approaches, mostly hidden Markov models (HMMs), which considered up to 15 IBD configurations among the four parental chromosomes, as well as more computationally efficient rule-based approaches. The accuracy of the methods was then evaluated, including with low-density marker panels, genotyping-by-sequencing data and small groups of individuals, typical features in such populations. Two HMMs performed consistently well, while two rule-based approaches proved efficient for genome-wide predictions. The model-based approaches were particularly efficient when information was reduced (low marker density, locus-specific estimation). Methods using phased data proved to be more efficient, while some approaches relying on unphased genotype data proved to be sensitive to the allele frequencies used. In some settings, pedigree information was competitive in predicting recent inbreeding levels. Finally, we showed that our evaluation is also informative about the accuracy of the methods for estimating relatedness and identifying IBD segments between pairs of individuals.

The feral cattle of the subantarctic island of Amsterdam provide an outstanding case study of a large mammalian population that was established by a handful of founders and thrived within a few generations in a seemingly inhospitable environment. Here, we investigated the genetic history and composition of this population using genotyping and sequencing data. Our inference showed an intense but brief founding bottleneck around the late 19 ^th century and revealed contributions from European taurine and Indian Ocean zebu in the founder ancestry. Comparative analysis of whole genome sequences further revealed a moderate reduction in genetic diversity despite high levels of inbreeding. The brief and intense bottleneck was associated with high levels of drift, a flattening of the site frequency spectrum and a slight relaxation of purifying selection on mildly deleterious variants. However, we did not observe any significant purge of highly deleterious variants. Interestingly, the population’s success in the harsh environment can be attributed to pre-adaptation from their European taurine ancestry, suggesting no strong bioclimatic challenge, and also contradicting evidence for insular dwarfism. Genome scan for footprints of selection uncovered a majority of candidate genes related to nervous system function, likely reflecting rapid feralization driven by behavioral changes and complex social restructuring. This unprecedented case study provides valuable insights into rapid population establishment, feralization, and genetic adaptation in challenging environments. It highlights the importance of preserving the unique genetic legacies of feral populations and raises ethical questions in the eyes of conservation efforts.

We report the generation of an organism-wide catalog of 976,813cis-acting regulatory elements for the bovine detected by the assay for transposase accessible chromatin using sequencing (ATAC-seq). We regroup these regulatory elements in 16 components by nonnegative matrix factorization. Correlation between the genome-wide density of peaks and transcription start sites, correlation between peak accessibility and expression of neighboring genes, and enrichment in transcription factor binding motifs support their regulatory potential. Using a previously established catalog of 12,736,643 variants, we show that the proportion of single-nucleotide polymorphisms mapping to ATAC-seq peaks is higher than expected and that this is owing to an approximately 1.3-fold higher mutation rate within peaks. Their site frequency spectrum indicates that variants in ATAC-seq peaks are subject to purifying selection. We generate eQTL data sets for liver and blood and show that variants that drive eQTL fall into liver- and blood-specific ATAC-seq peaks more often than expected by chance. We combine ATAC-seq and eQTL data to estimate that the proportion of regulatory variants mapping to ATAC-seq peaks is approximately one in three and that the proportion of variants mapping to ATAC-seq peaks that are regulatory is approximately one in 25. We discuss the implication of these findings on the utility of ATAC-seq information to improve the accuracy of genomic selection.

Cohorts of individuals that have been genotyped and phenotyped for genomic selection programs offer the opportunity to better understand genetic variation associated with complex traits. Here, we perform an association study for traits related to body size and muscular development in intensively selected beef cattle. We leveraged multiple trait information to refine and interpret the significant associations. After a multiple-step genotype imputation to the sequence-level for 14,762 Belgian Blue beef (BBB) cattle cows, we performed a GWAS for 11 traits related to muscular development and body size. The 37 identified genome-wide significant QTL could be condensed in 11 unique QTL regions based on their position. There was evidence for pleiotropic effects in most of these regions (e.g., correlated association signals, overlap between credible sets of candidate variants – CSCV). We consequently applied a multiple-trait approach to combine information from different traits to refine the CSCV. In several QTL regions, we identified strong candidate genes known to be related to growth and height in other species such as LCORL-NCAPG or CCND2 . For some of these genes, relevant candidate variants were identified in the CSCV, including three new missense variants in EZH2 , PAPPA2 and ADAM12 , possibly two additional coding variants in LCORL , and candidate regulatory variants linked to CCND2 and ARMC12 . Strikingly, four other QTL regions were related to five (recessive) deleterious coding variants previously identified. Heterozygotes for several of these mutations have favorable effects for muscular development traits. Our study further supports that a set of common genes controls body size across mammalian species. In particular, we added new genes to the list of those associated with height in both human and cattle. We also identified new strong candidate causing variants in some of those genes, strengthening the evidence of the causality of these genes. Several breed-specific recessive deleterious variants were identified in our QTL regions, probably as a result of the extreme selection for muscular development in BBB cattle.

Natural variation among individuals and populations exists in all species, playing key roles in response to environmental stress and adaptation. Micro- and macro-nutrients have a wide range of functions in photosynthetic organisms and mineral nutrition plays thus a sizable role in biomass production. To maintain nutrient concentrations inside the cell within physiological limits and prevent the detrimental effects of deficiency or excess, complex homeostatic networks have evolved in photosynthetic cells. The microalga Chlamydomonas reinhardtii (Chlamydomonas) is a unicellular eukaryotic model for studying such mechanisms. In this work, twenty-four Chlamydomonas strains, comprising field isolates and laboratory strains, were examined for intraspecific differences in nutrient homeostasis. Growth and mineral content were quantified in mixotrophy, as full nutrition control, and compared to autotrophy and 9 deficiency conditions for macronutrients (-Ca, -Mg, -N, -P, -S) and micronutrients (-Cu, -Fe, -Mn, -Zn). Growth differences among strains were relatively limited. However, similar growth was accompanied by highly divergent mineral accumulation among strains. The expression of nutrient status marker genes and photosynthesis were scored in pairs of contrasting field strains, revealing distinct transcriptional regulation and nutrient requirements. Leveraging this natural variation should enable a better understanding of nutrient homeostasis in Chlamydomonas.