Pedigree reconstruction from SNP data: parentage assignment, sibship clustering and beyond

MalKinID: A Likelihood-Based Model for Identifying Malaria Parasite Genealogical Relationships Using Identity-by-Descent

Pathogen genomics is a powerful tool for tracking infectious disease transmission. In malaria, identity-by-descent (IBD) is used to assess the genetic relatedness between parasites and has been used to study transmission and importation. In theory, IBD can be used to distinguish genealogical relationships to reconstruct transmission history or identify parasites for genotype-to-phenotype quantitative-trait-locus experiments. MalKinID (Malaria Kinship Identifier) is a new likelihood-based classification model designed to identify genealogical relationships among malaria parasites based on genome-wide IBD proportions and IBD segment distributions. MalKinID was calibrated to the genomic data from three laboratory-based genetic crosses (yielding 440 parent-child and 9060 full-sibling comparisons). MalKinID identified lab generated F1 progeny with >80% sensitivity and showed that 0.39 (95% CI 0.28, 0.49) of the second-generation progeny of a NF54 and NHP4026 cross were F1s and 0.56 (0.45, 0.67) were backcrosses of an F1 with the parental NF54 strain. In simulated outcrossed importations, MalKinID accurately reconstructs genealogy history with high precision and sensitivity, with F1-scores exceeding 0.84. However, when importation involves inbreeding, such as during serial co-transmission, the precision and sensitivity of MalKinID declined, with F1-scores of 0.76 (0.56, 0.92) and 0.23 (0.0, 0.4) for PC and FS and <0.05 for second-degree and third-degree relatives. Genealogical inference is most powered 1) when outcrossing is the norm or 2) when multi-sample comparisons based on a predefined pedigree are used. MalKinID lays the foundations for using IBD to track parasite transmission history and for separating progeny for quantitative-trait-locus experiments.

Patterns of reproduction and autozygosity distinguish the breeding from nonbreeding gray wolves of Yellowstone National Park

For species of management concern, accurate estimates of inbreeding and associated consequences on reproduction are crucial for predicting their future viability. However, few studies have partitioned this aspect of genetic viability with respect to reproduction in a group-living social mammal. We investigated the contributions of foundation stock lineages, putative fitness consequences of inbreeding, and genetic diversity of the breeding versus nonreproductive segment of the Yellowstone National Park gray wolf population. Our dataset spans 25 years and seven generations since reintroduction, encompassing 152 nuclear families and 329 litters. We found more than 87% of the pedigree foundation genomes persisted and report influxes of allelic diversity from two translocated wolves from a divergent source in Montana. As expected for group-living species, mean kinship significantly increased over time but with minimal loss of observed heterozygosity. Strikingly, the reproductive portion of the population carried a significantly lower genome-wide inbreeding coefficients, autozygosity, and more rapid decay for linkage disequilibrium relative to the nonbreeding population. Breeding wolves had significantly longer lifespans and lower inbreeding coefficients than nonbreeding wolves. Our model revealed that the number of litters was negatively significantly associated with heterozygosity (R = -0.11). Our findings highlight genetic contributions to fitness, and the importance of the reproductively active individuals in a population to counteract loss of genetic variation in a wild, free-ranging social carnivore. It is crucial for managers to mitigate factors that significantly reduce effective population size and genetic connectivity, which supports the dispersion of genetic variation that aids in rapid evolutionary responses to environmental challenges.

Group vocal composition and decision-making during roost finding in Spix's disk-winged bats

Theoretical work suggests that having many informed individuals within social groups can promote efficient resource location. However, it may also give rise to group fragmentation if members fail to reach consensus on their direction of movement. In this study, we investigate whether the number of informed individuals, exemplified by bats emitting calls from different roosts, influences group cohesion in Spix's disk-winged bats (Thyroptera tricolor). Additionally, we explore the role of signal reliability, quantified through signalling rates, in group consensus on where to roost. These bats use contact calls to announce the location of a roost site and recruit conspecifics. The groups they form exhibit high levels of cohesion and consist of both vocal and non-vocal bats, with vocal behaviour being consistent over time. Our findings revealed that an increase in the number of roosts broadcasting calls is strongly associated with the likelihood of groups fragmenting among multiple roosts. Additionally, we found that a majority of group members enter the roost with higher calling rates. This phenomenon can mitigate the risk of group fragmentation, as bats emitting more calls may contribute to greater group consensus on roosting locations, thereby reducing the likelihood of individuals separating and enhancing overall group cohesion. Our results highlight the potential costs of having too many information producers for group coordination, despite their established role in finding critical resources. This article is part of the theme issue 'The power of sound: unravelling how acoustic communication shapes group dynamics'.

Calling to the collective: contact calling rates within groups of disc-winged bats do not vary by kinship or association

Many group-living animals coordinate social behaviours using contact calls, which can be produced for all group members or targeted at specific individuals. In the disc-winged bat, Thyroptera tricolor, group members use 'inquiry' and 'response' calls to coordinate daily movements into new roosts (furled leaves). Rates of both calls show consistent among-individual variation, but causes of within-individual variation remain unknown. Here, we tested whether disc-winged bats produce more contact calls towards group members with higher kinship or association. In 446 experimental trials, we recorded 139 random within-group pairs of one flying bat (producing inquiry calls for roost searching) and one roosting bat (producing response calls for roost advertising). Using generalized linear mixed-effect models (GLMM), we assessed how response and inquiry calling rates varied by sender, receiver, genetic kinship and co-roosting association rate. Calling rates varied consistently across senders but not by receiver. Response calling was influenced by inquiry calling rates, but neither calling rate was higher when the interacting pair had higher kinship or association. Rather than dyadic calling rates indicating within-group relationships, our findings are consistent with the hypothesis that bats produce contact calls to maintain contact with any or all individuals within a group while collectively searching for a new roost site. This article is part of the theme issue 'The power of sound: unravelling how acoustic communication shapes group dynamics'.

Sex-specific overdominance at the maturation vgll3 gene for reproductive fitness in wild Atlantic salmon

Linking reproductive fitness with adaptive traits at the genomic level can shed light on the mechanisms that produce and maintain sex-specific selection. Here, we construct a multigenerational pedigree to investigate sex-specific selection on a maturation gene, vgll3, in a wild Atlantic salmon population. The vgll3 locus is responsible for ~40% of the variation in maturation (sea age at first reproduction). Genetic parentage analysis was conducted on 18,265 juveniles (parr) and 685 adults collected at the same spawning ground over eight consecutive years. A high proportion of females (26%) were iteroparous and reproduced two to four times in their lifetime. A smaller proportion of males (9%) spawned at least twice in their lifetime. Sex-specific patterns of reproductive fitness were related to vgll3 genotype. Females showed a pattern of overdominance where vgll3*EL genotypes had three-fold more total offspring than homozygous females. In contrast, males demonstrated that late-maturing vgll3*LL individuals had two-fold more offspring than either vgll3*EE or vgll3*EL males. Taken together, these data suggest that balancing selection in females contributes to the maintenance of variation at this locus via increased fitness of iteroparous vgll3*EL females. This study demonstrates the utility of multigenerational pedigrees for uncovering complex patterns of reproduction, sex-specific selection and the maintenance of genetic variation.

Translating genomic advances into biodiversity conservation

A key action of the new Global Biodiversity Framework is the maintenance of genetic diversity in all species to safeguard their adaptive potential. To achieve this goal, a translational mindset, which aims to convert results of basic research into direct practical benefits, needs to be applied to biodiversity conservation. Despite much discussion on the value of genomics to conservation, a disconnect between those generating genomic resources and those applying it to biodiversity management remains. As global efforts to generate reference genomes for non-model species increase, investment into practical biodiversity applications is critically important. Applications such as understanding population and multispecies diversity and longitudinal monitoring need support alongside education for policymakers on integrating the data into evidence-based decisions. Without such investment, the opportunity to revolutionize global biodiversity conservation using genomics will not be fully realized.

Genetic architecture of inbreeding depression may explain its persistence in a population of wild red deer

Inbreeding depression is of major concern in declining populations, but relatively little is known about its genetic architecture in wild populations, such as the degree to which it is composed of large or small effect loci and their distribution throughout the genome. Here, we combine fitness and genomic data from a wild population of red deer to investigate the genomic distribution of inbreeding effects. Based on the runs of homozygosity (ROH)-based inbreeding coefficient, FROH, we use chromosome-specific inbreeding coefficients (FROHChr) to explore whether the effect of inbreeding varies between chromosomes. Under the assumption that within an individual the probability of being identical-by-descent is equal across all chromosomes, we used a multi-membership model to estimate the deviation of FROHChr from the average inbreeding effect. This novel approach ensures effect sizes are not overestimated whilst maximising the power of our available dataset of >3000 individuals genotyped on >35,000 autosomal SNPs. We find that most chromosomes confer a minor reduction in fitness-related traits, which when these effects are summed, results in the observed inbreeding depression in birth weight, survival and lifetime breeding success. However, no chromosomes had a significant detrimental effect compared to the overall effect of inbreeding, indicating no major effect loci. We conclude that in this population, inbreeding depression is likely the result of multiple mildly or moderately deleterious mutations spread across all chromosomes, which are difficult to detect with statistical confidence. Such mutations will be inefficiently purged, which may explain the persistence of inbreeding depression in this population.

Collecting baleen whale blow samples by drone: A minimally intrusive tool for conservation genetics

In coastal British Columbia, Canada, marine megafauna such as humpback whales (Megaptera novaeangliae) and fin whales (Balaenoptera physalus velifera) have been subject to a history of exploitation and near extirpation. While their populations have been in recovery, significant threats are posed to these vulnerable species by proposed natural resource ventures in this region, in addition to the compounding effects of anthropogenic climate change. Genetic tools play a vital role in informing conservation efforts, but the associated collection of tissue biopsy samples can be challenging for the investigators and disruptive to the ongoing behaviour of the targeted whales. Here, we evaluate a minimally intrusive approach based on collecting exhaled breath condensate, or respiratory 'blow' samples, from baleen whales using an unoccupied aerial system (UAS), within Gitga'at First Nation territory for conservation genetics. Minimal behavioural responses to the sampling technique were observed, with no response detected 87% of the time (of 112 UAS deployments). DNA from whale blow (n = 88 samples) was extracted, and DNA profiles consisting of 10 nuclear microsatellite loci, sex identification and mitochondrial (mt) DNA haplotypes were constructed. An average of 7.5 microsatellite loci per individual were successfully genotyped. The success rates for mtDNA and sex assignment were 80% and 89% respectively. Thus, this minimally intrusive sampling method can be used to describe genetic diversity and generate genetic profiles for individual identification. The results of this research demonstrate the potential of UAS-collected whale blow for conservation genetics from a remote location.

The origin, connectivity, and individual specialization of island wolves after deer extirpation

Wolves are assumed to be ungulate obligates, however, a recently described pack on Pleasant Island, Alaska USA, is persisting on sea otters and other marine resources without ungulate prey, violating this long-held assumption. We address questions about these wolves regarding their origin and fate, degree of isolation, risk of inbreeding depression, and diet specialization by individual and sex. We applied DNA metabarcoding and genotyping by amplicon sequencing using 957 scats collected from 2016 to 2022, and reduced representation sequencing of tissue samples to establish a detailed understanding of Pleasant Island wolf ecology and compare them with adjacent mainland wolves. Dietary overlap was higher among individual wolves on Pleasant Island (Pianka's index mean 0.95 ± 0.03) compared to mainland wolves (0.70 ± 0.21). The individual diets of island wolves were dominated by sea otter, ranging from 40.6% to 63.2% weighted percent of occurrence (wPOO) (mean 55.5 ± 8.7). In contrast, individual mainland wolves primarily fed on ungulates (42.2 ± 21.3) or voles during a population outbreak (31.2 ± 23.2). We traced the origin of the Pleasant Island pack to a mainland pair that colonized around 2013 and produced several litters. After this breeding pair was killed, their female offspring and an immigrant male became the new breeders in 2019. We detected 20 individuals of which 8 (40%) were trapped and killed while two died of natural causes during the 6-year study. Except for the new breeding male, the pedigree analysis and genotype results showed no additional movement to or from the island, indicating limited dispersal but no evidence of inbreeding. Our findings suggest wolves exhibit more flexible foraging behavior than previously believed, and hunting strategies can substantially differ between individuals within or between packs. Nevertheless, anthropogenic and natural mortality combined with limited connectivity to the mainland may inhibit the continued persistence of Pleasant Island wolves.

Investigating pedigree- and SNP-associated components of heritability in a wild population of Soay sheep

Estimates of narrow sense heritability derived from genomic data that contain related individuals may be biased due to the within-family effects such as dominance, epistasis and common environmental factors. However, for many wild populations, removal of related individuals from the data would result in small sample sizes. In 2013, Zaitlen et al. proposed a method to estimate heritability in populations that include close relatives by simultaneously fitting an identity-by-state (IBS) genomic relatedness matrix (GRM) and an identity-by-descent (IBD) GRM. The IBD GRM is identical to the IBS GRM, except relatedness estimates below a specified threshold are set to 0. We applied this method to a sample of 8557 wild Soay sheep from St. Kilda, with genotypic information for 419,281 single nucleotide polymorphisms. We aimed to see how this method would partition heritability into population-level (IBS) and family-associated (IBD) variance for a range of genetic architectures, and so we focused on a mixture of polygenic and monogenic traits. We also implemented a variant of the model in which the IBD GRM was replaced by a GRM constructed from SNPs with low minor allele frequency to examine whether any additive genetic variance is captured by rare alleles. Whilst the inclusion of the IBD GRM did not significantly improve the fit of the model for the monogenic traits, it improved the fit for some of the polygenic traits, suggesting that dominance, epistasis and/or common environment not already captured by the non-genetic random effects fitted in our models may influence these traits.

Ultimate drivers of forced extra-pair copulations in birds lacking a penis: jackdaws as a case-study

Forced copulation is common, presumably because it can increase male reproductive success. Forced extra-pair copulation (FEPC) occurs in birds, even though most species lack a penis and are widely thought to require female cooperation for fertilization. How FEPC persists, despite a presumed lack of siring success and likely non-negligible costs to the male, is unknown. Using the jackdaw (Corvus monedula) as a case study, we use SNPs to quantify the extra-pair paternity rate through FEPC and evaluate explanations for the persistence of FEPC in species without a penis. We then collate evidence for FEPC across penis-lacking birds. Combining genetic and behavioural analyses, our study suggests that the most likely explanations for the maintenance of FEPC in jackdaws are that it provides a selective advantage to males or it is a relic. Our literature review shows that across birds lacking a penis, FEPC is taxonomically widespread, and yet, little is known about its evolution. A broader implementation of the approach used here, combining both genetic and behavioural data, may shed light on why this widespread sexual behaviour persists. Additional work is necessary to understand whether a penis is needed for paternity through forced copulation and to quantify the costs of FEPC.

How Can Genomics Help or Hinder Wildlife Conservation?

Genomic data are becoming increasingly affordable and easy to collect, and new tools for their analysis are appearing rapidly. Conservation biologists are interested in using this information to assist in management and planning but are typically limited financially and by the lack of genomic resources available for non-model taxa. It is therefore important to be aware of the pitfalls as well as the benefits of applying genomic approaches. Here, we highlight recent methods aimed at standardizing population assessments of genetic variation, inbreeding, and forms of genetic load and methods that help identify past and ongoing patterns of genetic interchange between populations, including those subjected to recent disturbance. We emphasize challenges in applying some of these methods and the need for adequate bioinformatic support. We also consider the promises and challenges of applying genomic approaches to understand adaptive changes in natural populations to predict their future adaptive capacity.

PMSeeker: A Scheme Based on the Greedy Algorithm and the Exhaustive Algorithm to Screen Low-Redundancy Marker Sets for Large-Scale Parentage Assignment with Full Parental Genotyping

Parentage assignment is a genetic test that utilizes genetic characteristics, such as molecular markers, to identify the parental relationships within populations, which, in commercial fish farming, are almost always large and where full information on potential parents is known. To accurately find the true parents, the genotypes of all loci in the parentage marker set (PMS) are required for each individual being tested. With the same accuracy, a PMS containing a smaller number of markers will undoubtedly save experimental costs. Thus, this study established a scheme to screen low-redundancy PMSs using the exhaustive algorithm and greedy algorithm. When screening PMSs, the greedy algorithm selects markers based on the parental dispersity index (PDI), a uniquely defined metric that outperforms the probability of exclusion (PE). With the conjunctive use of the two algorithms, non-redundant PMSs were found for more than 99.7% of solvable cases in three groups of random sample experiments in this study. Then, a low-redundancy PMS can be composed using two or more of these non-redundant PMSs. This scheme effectively reduces the number of markers in PMSs, thus conserving human and experimental resources and laying the groundwork for the widespread implementation of parentage assignment technology in economic species breeding.

Heritability and developmental plasticity of growth in an oviparous lizard

Selective processes act on phenotypic variation although the evolutionary potential of a trait relies on the underlying heritable variation. Developmental plasticity is an important source of phenotypic variation, but it can also promote changes in genetic variation, yet we have a limited understanding of how they are both impacted. Here, we quantified the influence of developmental temperature on growth in delicate skinks (Lampropholis delicata) and partitioned total phenotypic variance using an animal model fitted with a genomic relatedness matrix. We measured mass for 261 individuals (nhot = 125, ncold = 136) over 16 months (nobservations = 3002) and estimated heritability and maternal effects over time. Our results show that lizards reared in cold developmental temperatures had consistently higher mass across development compared to lizards that were reared in hot developmental temperatures. However, developmental temperature did not impact the rate of growth. On average, additive genetic variance, maternal effects and heritability were higher in the hot developmental temperature treatment; however, these differences were not statistically significant. Heritability increased with age, whereas maternal effects decreased upon hatching but increased again at a later age, which could be driven by social competition or intrinsic changes in the expression of variation as an individual's growth. Our work suggests that the evolutionary potential of growth is complex, age-dependent and not overtly affected by extremes in natural nest temperatures.

Development and evaluation of a novel single nucleotide polymorphism panel for North American bison

Genome-wide single nucleotide polymorphism (SNP) genotyping platforms have become increasingly popular in characterizing livestock and wildlife populations, replacing traditional methods such as microsatellite fragment analysis. Herein, we report the development and evaluation of a novel bison SNP panel for population management and conservation. Initially, 2474 autosomal SNPs were selected from existing bison whole-genome sequences and variable sites among bison on the GGSP bovine 50K Chip, based on minor allele frequency, data completeness, and chromosome location. Additionally, 20 mitochondrial SNPs were chosen to identify known mitochondrial haplotypes in bison according to previous research. The SNPs were further evaluated using genotyping-by-sequencing with 190 bison, representing the historical lineages that survived the major population crash of the late 1800s. Variants with high potential for genotyping error were filtered out, and the remaining SNPs were placed on a custom Illumina™ array. The final panel consisting of 798 autosomal and 13 mitochondrial SNPs was used to establish baseline genetic parameters, compare populations, and assign mitochondrial haplotypes in 995 bison across ten populations. These SNPs were also found to be highly informative for individual animal identification and parentage assignment. This SNP panel provides a powerful new method to establish a baseline for estimating genetic health of bison populations and a new tool for bison managers to make informed management decisions based on genetic information specific to their populations.

The estimation of additive genetic variance of body size in a wild passerine is sensitive to the method used to estimate relatedness among the individuals

Assessing additive genetic variance is a crucial step in predicting the evolutionary response of a target trait. However, the estimated genetic variance may be sensitive to the methodology used, e.g., the way relatedness is assessed among the individuals, especially in wild populations where social pedigrees can be inaccurate. To investigate this possibility, we investigated the additive genetic variance in tarsus length, a major proxy of skeletal body size in birds. The model species was the collared flycatcher (Ficedula albicollis), a socially monogamous but genetically polygamous migratory passerine. We used two relatedness matrices to estimate the genetic variance: (1) based solely on social links and (2) a genetic similarity matrix based on a large array of single-nucleotide polymorphisms (SNPs). Depending on the relatedness matrix considered, we found moderate to high additive genetic variance and heritability estimates for tarsus length. In particular, the heritability estimates were higher when obtained with the genetic similarity matrix instead of the social pedigree. Our results confirm the potential for this crucial trait to respond to selection and highlight methodological concerns when calculating additive genetic variance and heritability in phenotypic traits. We conclude that using a social pedigree instead of a genetic similarity matrix to estimate relatedness among individuals in a genetically polygamous wild population may significantly deflate the estimates of additive genetic variation.

Genetic diversity and breed-informative SNPs identification in domestic pig populations using coding SNPs

Background: The use of breed-informative genetic markers, specifically coding Single Nucleotide Polymorphisms (SNPs), is crucial for breed traceability, authentication of meat and dairy products, and the preservation and improvement of pig breeds. By identifying breed informative markers, we aimed to gain insights into the genetic mechanisms that influence production traits, enabling informed decisions in animal management and promoting sustainable pig production to meet the growing demand for animal products. Methods: Our dataset consists of 300 coding SNPs genotyped from three Italian commercial pig populations: Landrace, Yorkshire, and Duroc. Firstly, we analyzed the genetic diversity among the populations. Then, we applied a discriminant analysis of principal components to identify the most informative SNPs for discriminating between these populations. Lastly, we conducted a functional enrichment analysis to identify the most enriched pathways related to the genetic variation observed in the pig populations. Results: The alpha diversity indexes revealed a high genetic diversity within the three breeds. The higher proportion of observed heterozygosity than expected revealed an excess of heterozygotes in the populations that was supported by negative values of the fixation index (FIS) and deviations from the Hardy-Weinberg equilibrium. The Euclidean distance, the pairwise FST, and the pairwise Nei's GST genetic distances revealed that Yorkshire and Landrace breeds are genetically the closest, with distance values of 2.242, 0.029, and 0.033, respectively. Conversely, Landrace and Duroc breeds showed the highest genetic divergence, with distance values of 2.815, 0.048, and 0.052, respectively. We identified 28 significant SNPs that are related to phenotypic traits and these SNPs were able to differentiate between the pig breeds with high accuracy. The Functional Enrichment Analysis of the informative SNPs highlighted biological functions related to DNA packaging, chromatin integrity, and the preparation of DNA into higher-order structures. Conclusion: Our study sheds light on the genetic underpinnings of phenotypic variation among three Italian pig breeds, offering potential insights into the mechanisms driving breed differentiation. By prioritizing breed-specific coding SNPs, our approach enables a more focused analysis of specific genomic regions relevant to the research question compared to analyzing the entire genome.

Adaptive significance of affiliative behaviour differs between sexes in a wild reptile population

In recent years, we have begun to appreciate that social behaviours might exhibit repeatable among-individual variation. Such behavioural traits may even covary and have critical evolutionary implications. Importantly, some social behaviours such as aggressiveness have been shown to provide fitness benefits, including higher reproductive success and survival. However, fitness consequences of affiliative behaviours, especially between or among sexes, can be more challenging to establish. Using a longitudinal behavioural dataset (2014-2021) collected on eastern water dragons (Intellagama lesueurii), we investigated whether various aspects of affiliative behaviour (i) were repeatable across years, (ii) covaried with each other at the among-individual level, and (iii) influenced individuals' fitness. In particular, we considered affiliative behaviours towards opposite-sex and same-sex conspecifics separately. We found that social traits were repeatable and covaried with each other similarly for both sexes. More notably, we found that male reproductive success was positively correlated with the number of female associates and the proportion of time spent with females, while females' reproductive success was not correlated with any of the measured social behaviour metrics. Overall, these findings suggest that selection may be acting differently on social behaviour of male and female eastern water dragons.

Genomics in animal breeding from the perspectives of matrices and molecules

BACKGROUND

This paper describes genomics from two perspectives that are in use in animal breeding and genetics: a statistical perspective concentrating on models for estimating breeding values, and a sequence perspective concentrating on the function of DNA molecules.

MAIN BODY

This paper reviews the development of genomics in animal breeding and speculates on its future from these two perspectives. From the statistical perspective, genomic data are large sets of markers of ancestry; animal breeding makes use of them while remaining agnostic about their function. From the sequence perspective, genomic data are a source of causative variants; what animal breeding needs is to identify and make use of them.

CONCLUSION

The statistical perspective, in the form of genomic selection, is the more applicable in contemporary breeding. Animal genomics researchers using from the sequence perspective are still working towards this the isolation of causative variants, equipped with new technologies but continuing a decades-long line of research.

Correlates of early reproduction and apparent fitness consequences in male Soay sheep

Life history trade-offs are ubiquitous across species and place constraints on the timing of life history events, including the optimal age at first reproduction. However, studies on lifetime breeding success of male mammals are rare due to sex-biased dispersal and the requirement for genetic paternity inferences. We studied the correlates and apparent fitness consequences of early life reproduction among males in a free-living population of Soay sheep (Ovis aries) on St Kilda, Scotland. We investigated the factors associated with early breeding success and the apparent consequences of early success for survival and future reproduction. We used genetic paternity inferences, population data, and individual morphology measurements collected over 30 years. We found that individuals born in years with low-density population size had the highest early life breeding success and singletons were more likely to be successful than twins. Individuals that bred successfully at 7 months were more likely to survive their first winter. For individuals that survived their first winter, early breeding success was not associated with later breeding success. As individual heterogeneity affects breeding success, we believe that variation in individual quality masks the costs of early reproduction in this population. Our findings provide no evidence for selection for delayed age at reproduction in male Soay sheep.

Effects of an IgE receptor polymorphism acting on immunity, susceptibility to infection, and reproduction in a wild rodent

The genotype of an individual is an important predictor of their immune function, and subsequently, their ability to control or avoid infection and ultimately contribute offspring to the next generation. However, the same genotype, subjected to different intrinsic and/or extrinsic environments, can also result in different phenotypic outcomes, which can be missed in controlled laboratory studies. Natural wildlife populations, which capture both genotypic and environmental variability, provide an opportunity to more fully understand the phenotypic expression of genetic variation. We identified a synonymous polymorphism in the high-affinity Immunoglobulin E (IgE) receptor (GC and non-GC haplotypes) that has sex-dependent effects on immune gene expression, susceptibility to infection, and reproductive success of individuals in a natural population of field voles (Microtus agrestis). We found that the effect of the GC haplotype on the expression of immune genes differed between sexes. Regardless of sex, both pro-inflammatory and anti-inflammatory genes were more highly relatively expressed in individuals with the GC haplotype than individuals without the haplotype. However, males with the GC haplotype showed a stronger signal for pro-inflammatory genes, while females showed a stronger signal for anti-inflammatory genes. Furthermore, we found an effect of the GC haplotype on the probability of infection with a common microparasite, Babesia microti, in females - with females carrying the GC haplotype being more likely to be infected. Finally, we found an effect of the GC haplotype on reproductive success in males - with males carrying the GC haplotype having a lower reproductive success. This is a rare example of a polymorphism whose consequences we are able to follow across immunity, infection, and reproduction for both males and females in a natural population.

Pedigree reconstruction and population structure using SNP markers in Gir cattle

Our objective was to establish a SNPs panel for pedigree reconstruction using microarrays of different densities and evaluate the genomic relationship coefficient of the inferred pedigree, in addition to analyzing the population structure based on genomic analyses in Gir cattle. For parentage analysis and genomic relationship, 16,205 genotyped Gir animals (14,458 females and 1747 males) and 1810 common markers to the four SNP microarrays were used. For population structure analyses, including linkage disequilibrium, effective population size, and runs of homozygosity (ROH), genotypes from 21,656 animals were imputed. Likelihood ratio (LR) approach was used to reconstruct the pedigree, deepening the pedigree and showing it is well established in terms of recent information. Coefficients for each relationship category of the inferred pedigree were adequate. Linkage disequilibrium showed rapid decay. We detected a decrease in the effective population size over the last 50 generations, with the average generation interval around 9.08 years. Higher ROH-based inbreeding coefficient in a class of short ROH segments, with moderate to high values, was also detected, suggesting bottlenecks in the Gir genome. Breeding strategies to minimize inbreeding and avoid massive use of few proven sires with high genetic value are suggested to maintain genetic variability in future generations. In addition, we recommend reducing the generation interval to maximize genetic progress and increase effective population size.

Wild snapdragon plant pedigree sheds light on limited connectivity enhanced by higher migrant reproductive success in a fragmented landscape

Background: In contrast with historical knowledge, a recent view posits that a non-negligible proportion of populations thrive in a fragmented landscape. One underlying mechanism is the maintenance of functional connectivity, i.e., the net flow of individuals or their genes moving among suitable habitat patches. Alternatively, functional connectivity might be typically limited but enhanced by a higher reproductive success of migrants. Methods: We tested for this hypothesis in wild snapdragon plants inhabiting six patches separated by seawater in a fragmented Mediterranean scrubland landscape. We reconstructed their pedigree by using a parentage assignment method based on microsatellite genetic markers. We then estimated functional connectivity and the reproductive success of plants resulting from between-patch dispersal events. Results: We found that wild snapdragon plants thrived in this fragmented landscape, although functional connectivity between habitat patches was low (i.e. 2.9%). The progeny resulting from between-patch dispersal events had a higher reproductive success than residents. Conclusion: Our findings imply that low functional connectivity in a fragmented landscapes may have been enhanced by higher reproductive success after migration. This original mechanisms might be partly compensating the negative impact of fragmentation.

Taking quantitative genomics into the wild

We organized this special issue to highlight new work and review recent advances at the cutting edge of 'wild quantitative genomics'. In this editorial, we will present some history of wild quantitative genetic and genomic studies, before discussing the main themes in the papers published in this special issue and highlighting the future outlook of this dynamic field.

Genetic diversity and connectivity of chemosynthetic cold seep mussels from the U.S. Atlantic margin

Background

Deep-sea mussels in the subfamily Bathymodiolinae have unique adaptations to colonize hydrothermal-vent and cold-seep environments throughout the world ocean. These invertebrates function as important ecosystem engineers, creating heterogeneous habitat and promoting biodiversity in the deep sea. Despite their ecological significance, efforts to assess the diversity and connectivity of this group are extremely limited. Here, we present the first genomic-scale diversity assessments of the recently discovered bathymodioline cold-seep communities along the U.S. Atlantic margin, dominated by Gigantidas childressi and Bathymodiolus heckerae.

Results

A Restriction-site Associated DNA Sequencing (RADSeq) approach was used on 177 bathymodiolines to examine genetic diversity and population structure within and between seep sites. Assessments of genetic differentiation using single-nucleotide polymorphism (SNP) data revealed high gene flow among sites, with the shallower and more northern sites serving as source populations for deeper occurring G. childressi. No evidence was found for genetic diversification across depth in G. childressi, likely due to their high dispersal capabilities. Kinship analyses indicated a high degree of relatedness among individuals, and at least 10–20% of local recruits within a particular site. We also discovered candidate adaptive loci in G. childressi and B. heckerae that suggest differences in developmental processes and depth-related and metabolic adaptations to chemosynthetic environments.

Conclusions

These results highlight putative source communities for an important ecosystem engineer in the deep sea that may be considered in future conservation efforts. Our results also provide clues into species-specific adaptations that enable survival and potential speciation within chemosynthetic ecosystems.

The Unusual Value of Long-Term Studies of Individuals: The Example of the Isle of Rum Red Deer Project

Long-term studies of individuals enable incisive investigations of questions across ecology and evolution. Here, we illustrate this claim by reference to our long-term study of red deer on the Isle of Rum, Scotland. This project has established many of the characteristics of social organization, selection, and population ecology typical of large, polygynous, seasonally breeding mammals, with wider implications for our understanding of sexual selection and the evolution of sex differences, as well as for their population dynamics and population management. As molecular genetic techniques have developed, the project has pivoted to investigate evolutionary genetic questions, also breaking new ground in this field. With ongoing advances in genomics and statistical approaches and the development of increasingly sophisticated ways to assay new phenotypic traits, the questions that long-term studies such as the red deer study can answer become both broader and ever more sophisticated. They also offer powerful means of understanding the effects of ongoing climate change on wild populations.

Maternal effects do not resolve the paradox of stasis in birth weight in a wild red deer populaton

In natural populations, quantitative traits seldom show short-term evolution at the rate predicted by evolutionary models. Resolving this "paradox of stasis" is a key goal in evolutionary biology, as it directly challenges our capacity to predict evolutionary change. One particularly promising hypothesis to explain the lack of evolutionary responses in a key offspring trait, body weight, is that positive selection on juveniles is counterbalanced by selection against maternal investment in offspring growth, given that reproduction is costly for the mothers. Here, we used data from one of the longest individual-based studies of a wild mammal population to test this hypothesis. We first showed that despite positive directional selection on birth weight, and heritable variation for this trait, no genetic change has been observed for birth weight over the past 47 years in the study population. Contrarily to our expectation, we found no evidence of selection against maternal investment in birth weight-if anything, selection favors mothers that produce large calves. Accordingly, we show that genetic change in birth weight over the study period is actually lower than that predicted from models including selection on maternal performance; ultimately our analysis here only deepens rather than resolves the paradox of stasis.

Design and validation of a 63K genome-wide SNP-genotyping platform for caribou/reindeer (Rangifer tarandus)

BACKGROUND

Development of large single nucleotide polymorphism (SNP) arrays can make genomic data promptly available for conservation problematic. Medium and high-density panels can be designed with sufficient coverage to offer a genome-wide perspective and the generated genotypes can be used to assess different genetic metrics related to population structure, relatedness, or inbreeding. SNP genotyping could also permit sexing samples with unknown associated metadata as it is often the case when using non-invasive sampling methods favored for endangered species. Genome sequencing of wild species provides the necessary information to design such SNP arrays. We report here the development of a SNP-array for endangered Rangifer tarandus using a multi-platform sequencing approach from animals found in diverse populations representing the entire circumpolar distribution of the species.

RESULTS

From a very large comprehensive catalog of SNPs detected over the entire sample set (N = 894), a total of 63,336 SNPs were selected. SNP selection accounted for SNPs evenly distributed across the entire genome (~ every 50Kb) with known minor alleles across populations world-wide. In addition, a subset of SNPs was selected to represent rare and local alleles found in Eastern Canada which could be used for ecotype and population assignments - information urgently needed for conservation planning. In addition, heterozygosity from SNPs located in the X-chromosome and genotyping call-rate of SNPs located into the SRY gene of the Y-chromosome yielded an accurate and robust sexing assessment. All SNPs were validated using a high-throughput SNP-genotyping chip.

CONCLUSION

This design is now integrated into the first genome-wide commercially available genotyping platform for Rangifer tarandus. This platform would pave the way to future genomic investigation of populations for this endangered species, including estimation of genetic diversity parameters, population assignments, as well as animal sexing from genetic SNP data for non-invasive samples.

Genome-wide analysis of Schistosoma mansoni reveals limited population structure and possible praziquantel drug selection pressure within Ugandan hot-spot communities

Populations within schistosomiasis control areas, especially those in Africa, are recommended to receive regular mass drug administration (MDA) with praziquantel (PZQ) as the main strategy for controlling the disease. The impact of PZQ treatment on schistosome genetics remains poorly understood, and is limited by a lack of high-resolution genetic data on the population structure of parasites within these control areas. We generated whole-genome sequence data from 174 individual miracidia collected from both children and adults from fishing communities on islands in Lake Victoria in Uganda that had received either annual or quarterly MDA with PZQ over four years, including samples collected immediately before and four weeks after treatment. Genome variation within and between samples was characterised and we investigated genomic signatures of natural selection acting on these populations that could be due to PZQ treatment. The parasite population on these islands was more diverse than found in nearby villages on the lake shore. We saw little or no genetic differentiation between villages, or between the groups of villages with different treatment intensity, but slightly higher genetic diversity within the pre-treatment compared to post-treatment parasite populations. We identified classes of genes significantly enriched within regions of the genome with evidence of recent positive selection among post-treatment and intensively treated parasite populations. The differential selection observed in post-treatment and pre-treatment parasite populations could be linked to any reduced susceptibility of parasites to praziquantel treatment.

Integer programming for selecting set of informative markers in paternity inference

Background

Parentage information is fundamental to various life sciences. Recent advances in sequencing technologies have made it possible to accurately infer parentage even in non-model species. The optimization of sets of genome-wide markers is valuable for cost-effective applications but requires extremely large amounts of computation, which presses for the development of new efficient algorithms.

Results

Here, for a closed half-sib population, we generalized the process of marker loci selection as a binary integer programming problem. The proposed systematic formulation considered marker localization and the family structure of the potential parental population, resulting in an accurate assignment with a small set of markers. We also proposed an efficient heuristic approach, which effectively improved the number of markers, localization, and tolerance to missing data of the set. Applying this method to the actual genotypes of apple (Malus × domestica) germplasm, we identified a set of 34 SNP markers that distinguished 300 potential parents crossed to a particular cultivar with a greater than 99% accuracy.

Conclusions

We present a novel approach for selecting informative markers based on binary integer programming. Since the data generated by high-throughput sequencing technology far exceeds the requirement for parentage assignment, a combination of the systematic marker selection with targeted SNP genotyping, such as KASP, allows flexibly enlarging the analysis up to a scale that has been unrealistic in various species. The method developed in this study can be directly applied to unsolved large-scale problems in breeding, reproduction, and ecological research, and is expected to lead to novel knowledge in various biological fields. The implementation is available at https://github.com/SoNishiyama/IP-SIMPAT.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04801-z.

Genetic variance in fitness indicates rapid contemporary adaptive evolution in wild animals

The rate of adaptive evolution, the contribution of selection to genetic changes that increase mean fitness, is determined by the additive genetic variance in individual relative fitness. To date, there are few robust estimates of this parameter for natural populations, and it is therefore unclear whether adaptive evolution can play a meaningful role in short-term population dynamics. We developed and applied quantitative genetic methods to long-term datasets from 19 wild bird and mammal populations and found that, while estimates vary between populations, additive genetic variance in relative fitness is often substantial and, on average, twice that of previous estimates. We show that these rates of contemporary adaptive evolution can affect population dynamics and hence that natural selection has the potential to partly mitigate effects of current environmental change.

Diverse parentage relationships in paternal mouthbrooding fishes

While mouthbrooding is not an uncommon parental care strategy in fishes, paternal mouthbrooding only occurs in eight fish families and is little studied. The high cost of paternal mouthbrooding to the male implies a low risk of investment in another male's offspring but genetic parentage patterns are poorly known for paternal mouthbrooders. Here, we used single-nucleotide polymorphism genetic data to investigate parentage relationships of broods of two mouthbrooders of northern Australian rivers, mouth almighty Glossamia aprion and blue catfish Neoarius graeffei. For N. graeffei, we found that the parentage pattern was largely monogamous with the brooder male as the sire. For G. aprion, the parentage pattern was more heterogeneous including observations of monogamous broods with the brooder male as the sire (73%), polygyny (13%), cuckoldry (6%) and a brood genetically unrelated to the brooder male (6%). Findings demonstrate the potential for complex interrelationships of male care, paternity confidence and mating behaviour in mouthbrooding fishes.

Single‐parentage assignments reveal negative‐assortative mating in an endangered salmonid

Understanding reproductive patterns in endangered species is critical for supporting their recovery efforts. In this study we use a combination of paired‐parent and single‐parent assignments to examine the reproductive patterns in an endangered population of sockeye salmon (Oncorhynchus nerka) that uses Redfish Lake in central Idaho as a spawning and nursery lake. Recovery efforts include the release of maturing adults into the lake for volitional spawning. The lake is also inhabited by a population of resident O. nerka that is genetically indistinguishable, but phenotypically smaller, to the maturing adults released into the lake. The resident population is difficult to sample and the reproductive patterns between the two groups are unknown. We used results of paired‐ and single‐parentage assignments to specifically examine the reproductive patterns of male fish released into the lake under an equal sex ratio and a male‐biased sex ratio. Assignment results of offspring leaving the lake indicated a reproductive shift by males under the two scenarios. Males displayed an assortative mating pattern under an equal sex ratio and spawned almost exclusively with the released females. Under a male‐biased sex ratio most males shifted to a negative‐assortative mating pattern and spawned with smaller females from the resident population. These males were younger and smaller than males that spawned with released females suggesting they were unable to compete with larger males for spawning opportunities with the larger, released females. The results provided insights into the reproductive behavior of this endangered population and has implications for recovery efforts.

Unraveling the complex genetic basis of growth in New Zealand silver trevally (Pseudocaranx georgianus)

Growth directly influences production rate and therefore is one of the most important and well-studied traits in animal breeding. However, understanding the genetic basis of growth has been hindered by its typically complex polygenic architecture. Here, we performed quantitative trait locus mapping and genome-wide association studies for 10 growth traits that were observed over 2 years in 1,100 F1 captive-bred trevally (Pseudocaranx georgianus). We constructed the first high-density linkage map for trevally, which included 19,861 single nucleotide polymorphism markers, and discovered 8 quantitative trait loci for height, length, and weight on linkage groups 3, 14, and 18. Using genome-wide association studies, we further identified 113 single nucleotide polymorphism-trait associations, uncovering 10 genetic hot spots involved in growth. Two of the markers found in the genome-wide association studies colocated with the quantitative trait loci previously mentioned, demonstrating that combining quantitative trait locus mapping and genome-wide association studies represents a powerful approach for the identification and validation of loci controlling complex traits. This is the first study of its kind for trevally. Our findings provide important insights into the genetic architecture of growth in this species and supply a basis for fine mapping quantitative trait loci, genomic selection, and further detailed functional analysis of the genes underlying growth in trevally.

Rescued back from extinction in the wild: past, present and future of the genetics of the Arabian oryx in Oman

The Arabian oryx was the first species to be rescued from extinction in the wild by the concerted efforts of captive programmes in zoos and private collections around the world. Reintroduction efforts have used two main sources: the 'World Herd', established at the Phoenix Zoo, and private collections in Saudi Arabia. The breeding programme at the Al-Wusta Wildlife Reserve (WWR) in Oman has played a central role in the rescue of the oryx. Individuals from the 'World Herd' and the United Arab Emirates have been the main source for the WWR programme. However, no breeding strategies accounting for genetic diversity have been implemented. To address this, we investigated the diversity of the WWR population and historical samples using mitochondrial DNA (mtDNA) and single nucleotide polymorphisms (SNPs). We found individuals at WWR contain 58% of the total mtDNA diversity observed globally. Inference of ancestry and spatial patterns of SNP variation shows the presence of three ancestral sources and three different groups of individuals. Similar levels of diversity and low inbreeding were observed between groups. We identified individuals and groups that could most effectively contribute to maximizing genetic diversity. Our results will be valuable to guide breeding and reintroduction programmes at WWR.

Quantitative multidimensional phenotypes improve genetic analysis of laterality traits

Handedness is the most commonly investigated lateralised phenotype and is usually measured as a binary left/right category. Its links with psychiatric and neurodevelopmental disorders prompted studies aimed at understanding the underlying genetics, while other measures and side preferences have been less studied. We investigated the heritability of hand, as well as foot, and eye preference by assessing parental effects (n ≤ 5028 family trios) and SNP-based heritability (SNP-h², n ≤ 5931 children) in the Avon Longitudinal Study of Parents and Children (ALSPAC). An independent twin cohort from Hong Kong (n = 358) was used to replicate results from structural equation modelling (SEM). Parental left-side preference increased the chance of an individual to be left-sided for the same trait, with stronger maternal than paternal effects for footedness. By regressing out the effects of sex, age, and ancestry, we transformed laterality categories into quantitative measures. The SNP-h² for quantitative handedness and footedness was 0.21 and 0.23, respectively, which is higher than the SNP-h² reported in larger genetic studies using binary handedness measures. The heritability of the quantitative measure of handedness increased (0.45) compared to a binary measure for writing hand (0.27) in the Hong Kong twins. Genomic and behavioural SEM identified a shared genetic factor contributing to handedness, footedness, and eyedness, but no independent effects on individual phenotypes. Our analysis demonstrates how quantitative multidimensional laterality phenotypes are better suited to capture the underlying genetics than binary traits.

Genomic Predictions of Phenotypes and Pseudo-Phenotypes for Viral Nervous Necrosis Resistance, Cortisol Concentration, Antibody Titer and Body Weight in European Sea Bass

Simple Summary

Selective breeding programs based on genomic data are still not a common practice in aquaculture, although genomic selection has been widely demonstrated to be advantageous when trait phenotyping is a difficult task. In this study, we investigated the accuracy of predicting the phenotype and the estimated breeding value (EBV) of three Bayesian models and a Random Forest algorithm exploiting the information of a genome-wide SNP panel for European sea bass. The genomic predictions were developed for mortality caused by viral nervous necrosis, post-stress cortisol concentration, antibody titer against nervous necrosis virus and body weight. Selective breeding based on genomic data is a possible option for improving these traits while overcoming difficulties related to individual phenotyping of the investigated traits. Our results evidenced that the EBV used as a pseudo-phenotype enhances the predictive performances of genomic models, and that EBV can be predicted with satisfactory accuracy. The genomic prediction of the EBV for mortality might also be used to classify the phenotype for the same trait.

Abstract

In European sea bass (Dicentrarchus labrax L.), the viral nervous necrosis mortality (MORT), post-stress cortisol concentration (HC), antibody titer (AT) against nervous necrosis virus and body weight (BW) show significant heritability, which makes selective breeding a possible option for their improvement. An experimental population (N = 650) generated by a commercial broodstock was phenotyped for the aforementioned traits and genotyped with a genome-wide SNP panel (16,075 markers). We compared the predictive accuracies of three Bayesian models (Bayes B, Bayes C and Bayesian Ridge Regression) and a machine-learning method (Random Forest). The prediction accuracy of the EBV for MORT was approximately 0.90, whereas the prediction accuracies of the EBV and the phenotype were 0.86 and 0.21 for HC, 0.79 and 0.26 for AT and 0.71 and 0.38 for BW. The genomic prediction of the EBV for MORT used to classify the phenotype for the same trait showed moderate classification performance. Genome-wide association studies confirmed the polygenic nature of MORT and demonstrated a complex genetic structure for HC and AT. Genomic predictions of the EBV for MORT could potentially be used to classify the phenotype of the same trait, though further investigations on a larger experimental population are needed.

Threatened Species Initiative: Empowering conservation action using genomic resources

Globally, 15,521 animal species are listed as threatened by the International Union for the Conservation of Nature, and of these less than 3% have genomic resources that can inform conservation management. To combat this, global genome initiatives are developing genomic resources, yet production of a reference genome alone does not conserve a species. The reference genome allows us to develop a suite of tools to understand both genome-wide and functional diversity within and between species. Conservation practitioners can use these tools to inform their decision-making. But, at present there is an implementation gap between the release of genome information and the use of genomic data in applied conservation by conservation practitioners. In May 2020, we launched the Threatened Species Initiative and brought a consortium of genome biologists, population biologists, bioinformaticians, population geneticists, and ecologists together with conservation agencies across Australia, including government, zoos, and nongovernment organizations. Our objective is to create a foundation of genomic data to advance our understanding of key Australian threatened species, and ultimately empower conservation practitioners to access and apply genomic data to their decision-making processes through a web-based portal. Currently, we are developing genomic resources for 61 threatened species from a range of taxa, across Australia, with more than 130 collaborators from government, academia, and conservation organizations. Developed in direct consultation with government threatened-species managers and other conservation practitioners, herein we present our framework for meeting their needs and our systematic approach to integrating genomics into threatened species recovery.

Contemporary selection on MHC genes in a free-living ruminant population

Genes within the major histocompatibility complex (MHC) are the most variable identified in vertebrates. Pathogen-mediated selection is believed to be the main force maintaining MHC diversity. However, relatively few studies have demonstrated contemporary selection on MHC genes. Here, we examine associations between MHC variation and several fitness measurements including total fitness and five fitness components, in 3400 wild Soay sheep (Ovis aries) monitored between 1989 and 2012. In terms of total fitness, measured as lifetime breeding success of all individuals born, we found haplotypes named C and D were associated with decreased and increased male total fitness respectively. In terms of fitness components, juvenile survival was associated with haplotype divergence while individual haplotypes (C, D and F) were associated with adult fitness components. Consistent with the increased male total fitness, the rarest haplotype D has increased in frequency throughout the study period more than expected under neutral expectations. Our results demonstrate contemporary natural selection is acting on MHC class II genes in Soay sheep and the mode of selection on specific fitness components can be different mode from selection on total fitness.

How much is enough? Sampling intensity influences estimates of reproductive variance in an introduced population

Conservation introductions to islands and fenced enclosures are increasing as in situ mitigations fail to keep pace with population declines. Few studies consider the potential loss of genetic diversity and increased inbreeding if released individuals breed disproportionately. As funding is limited and post-release monitoring expensive for conservation programs, understanding how sampling effort influences estimates of reproductive variance is useful. To investigate this relationship, we used a well-studied population of Tasmanian devils (Sarcophilus harrisii) introduced to Maria Island, Tasmania, Australia. Pedigree reconstruction based on molecular data revealed high variance in number of offspring per breeder and high proportions of unsuccessful individuals. Computational subsampling of 20%, 40%, 60%, and 80% of observed offspring resulted in inaccurate estimates of reproductive variance compared to the pedigree reconstructed with all sampled individuals. With decreased sampling effort, the proportion of inferred unsuccessful individuals was overestimated and the variance in number of offspring per breeder was underestimated. To accurately estimate reproductive variance, we recommend sampling as many individuals as logistically possible during the early stages of population establishment. Further, we recommend careful selection of colonizing individuals as they may be disproportionately represented in subsequent generations. Within the conservation management context, our results highlight important considerations for sample collection and post-release monitoring during population establishment.

The genome sequence of the red deer, Cervus elaphus Linnaeus 1758

We present a genome assembly from an individual female Cervus elaphus (the red deer; Chordata; Mammalia; Artiodactyla; Cervidae). The genome sequence is 2,887 megabases in span. The majority of the assembly is scaffolded into 34 chromosomal pseudomolecules, with the X sex chromosome assembled.

Genomic analysis reveals a polygenic architecture of antler morphology in wild red deer (Cervus elaphus)

Sexually selected traits show large variation and rapid evolution across the animal kingdom, yet genetic variation often persists within populations despite apparent directional selection. A key step in solving this long-standing paradox is to determine the genetic architecture of sexually selected traits to understand evolutionary drivers and constraints at the genomic level. Antlers are a form of sexual weaponry in male red deer (Cervus elaphus). On the island of Rum, Scotland, males with larger antlers have increased breeding success, yet there has been no evidence of any response to selection at the genetic level. To try and understand the mechanisms underlying this observation, we investigate the genetic architecture of ten antler traits and their principal components using genomic data from >38,000 SNPs. We estimate the heritabilities and genetic correlations of the antler traits using a genomic relatedness approach. We then use genome-wide association and haplotype-based regional heritability to identify regions of the genome underlying antler morphology, and an empirical Bayes approach to estimate the underlying distributions of allele effect sizes. We show that antler morphology is highly repeatable over an individual's lifetime, heritable and has a polygenic architecture and that almost all antler traits are positively genetically correlated with some loci identified as having pleiotropic effects. Our findings suggest that a large mutational target and genetic covariances among antler traits, in part maintained by pleiotropy, are likely to contribute to the maintenance of genetic variation in antler morphology in this population.

Bonsai: An efficient method for inferring large human pedigrees from genotype data

Pedigree inference from genotype data is a challenging problem, particularly when pedigrees are sparsely sampled and individuals may be distantly related to their closest genotyped relatives. We present a method that infers small pedigrees of close relatives and then assembles them into larger pedigrees. To assemble large pedigrees, we introduce several formulas and tools including a likelihood for the degree separating two small pedigrees, a generalization of the fast DRUID point estimate of the degree separating two pedigrees, a method for detecting individuals who share background identity-by-descent (IBD) that does not reflect recent common ancestry, and a method for identifying the ancestral branches through which distant relatives are connected. Our method also takes several approaches that help to improve the accuracy and efficiency of pedigree inference. In particular, we incorporate age information directly into the likelihood rather than using ages only for consistency checks and we employ a heuristic branch-and-bound-like approach to more efficiently explore the space of possible pedigrees. Together, these approaches make it possible to construct large pedigrees that are challenging or intractable for current inference methods.

A Semi-Automated SNP-Based Approach for Contaminant Identification in Biparental Polyploid Populations of Tropical Forage Grasses

Artificial hybridization plays a fundamental role in plant breeding programs since it generates new genotypic combinations that can result in desirable phenotypes. Depending on the species and mode of reproduction, controlled crosses may be challenging, and contaminating individuals can be introduced accidentally. In this context, the identification of such contaminants is important to avoid compromising further selection cycles, as well as genetic and genomic studies. The main objective of this work was to propose an automated multivariate methodology for the detection and classification of putative contaminants, including apomictic clones (ACs), self-fertilized individuals, half-siblings (HSs), and full contaminants (FCs), in biparental polyploid progenies of tropical forage grasses. We established a pipeline to identify contaminants in genotyping-by-sequencing (GBS) data encoded as allele dosages of single nucleotide polymorphism (SNP) markers by integrating principal component analysis (PCA), genotypic analysis (GA) measures based on Mendelian segregation, and clustering analysis (CA). The combination of these methods allowed for the correct identification of all contaminants in all simulated progenies and the detection of putative contaminants in three real progenies of tropical forage grasses, providing an easy and promising methodology for the identification of contaminants in biparental progenies of tetraploid and hexaploid species. The proposed pipeline was made available through the polyCID Shiny app and can be easily coupled with traditional genetic approaches, such as linkage map construction, thereby increasing the efficiency of breeding programs.