A consensus genome of sika deer (Cervus nippon) and transcriptome analysis provided novel insights on the regulation mechanism of transcript factor in antler development

BACKGROUND

Sika deer (Cervus nippon) holds significance among cervids, with three genomes recently published. However, these genomes still contain hundreds of gaps and display significant discrepancies in continuity and accuracy. This poses challenges to functional genomics research and the selection of an appropriate reference genome. Thus, obtaining a high-quality reference genome is imperative to delve into functional genomics effectively.

FINDINGS

Here we report a high-quality consensus genome of male sika deer. All 34 chromosomes are assembled into single-contig pseudomolecules without any gaps, which is the most complete assembly. The genome size is 2.7G with 23,284 protein-coding genes. Comparative genomics analysis found that the genomes of sika deer and red deer are highly conserved, an approximately 2.4G collinear regions with up to 99% sequence similarity. Meanwhile, we observed the fusion of red deer's Chr23 and Chr4 during evolution, forming sika deer's Chr1. Additionally, we identified 607 transcription factors (TFs) that are involved in the regulation of antler development, including RUNX2, SOX6, SOX8, SOX9, PAX8, SIX2, SIX4, SIX6, SPI1, NFAC1, KLHL8, ZN710, JDP2, and TWST2, based on this consensus reference genome.

CONCLUSIONS

Our results indicated that we acquired a high-quality consensus reference genome. That provided valuable resources for understanding functional genomics. In addition, discovered the genetic basis of sika-red hybrid fertility and identified 607 significant TFs that impact antler development.

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.

De novo reconstruction of satellite repeat units from sequence data

Satellite DNA are long tandemly repeating sequences in a genome and may be organized as high-order repeats (HORs). They are enriched in centromeres and are challenging to assemble. Existing algorithms for identifying satellite repeats either require the complete assembly of satellites or only work for simple repeat structures without HORs. Here we describe Satellite Repeat Finder (SRF), a new algorithm for reconstructing satellite repeat units and HORs from accurate reads or assemblies without prior knowledge on repeat structures. Applying SRF to real sequence data, we show that SRF could reconstruct known satellites in human and well-studied model organisms. We also find satellite repeats are pervasive in various other species, accounting for up to 12% of their genome contents but are often underrepresented in assemblies. With the rapid progress in genome sequencing, SRF will help the annotation of new genomes and the study of satellite DNA evolution even if such repeats are not fully assembled.

Comparative Analysis of Microsatellite and SNP Markers for Genetic Management of Red Deer

The analysis of population genetic structure and individual multilocus heterozygosity are crucial for wildlife management and conservation. Microsatellite markers have traditionally been used to assess these genetic parameters. However, single-nucleotide polymorphisms (SNPs) are becoming increasingly popular. Our goal here was to determine to what extent SNPs can provide better insights than microsatellites into the overall genetic status and population genetic processes in the species. To this end, we genotyped 210 red deer (Cervus elaphus) in the Spanish wild population with both 11 microsatellites and 31,712 SNPs. We compared parameters related to population genetic structure and individual multilocus heterozygosity obtained with both types of markers. Our results showed correlations between parameters measured using both microsatellites and SNPs, particularly those related to the level of genetic diversity and genetic differentiation. However, we found notably lower precision of microsatellites in measuring the distribution of genetic diversity among individuals. We conclude that microsatellites can be used to monitor the overall genetic status and detect broad patterns in red deer populations. Nevertheless, the greater precision of SNPs in inferring genetic structure and multilocus heterozygosity leads us to encourage scientists and wildlife managers to prioritize their use whenever possible.

Comparative studies of X chromosomes in Cervidae family

The family Cervidae is the second most diverse in the infraorder Pecora and is characterized by variability in the diploid chromosome numbers among species. X chromosomes in Cervidae evolved through complex chromosomal rearrangements of conserved segments within the chromosome, changes in centromere position, heterochromatic variation, and X-autosomal translocations. The family Cervidae consists of two subfamilies: Cervinae and Capreolinae. Here we build a detailed X chromosome map with 29 cattle bacterial artificial chromosomes of representatives of both subfamilies: reindeer (Rangifer tarandus), gray brocket deer (Mazama gouazoubira), Chinese water deer (Hydropotes inermis) (Capreolinae); black muntjac (Muntiacus crinifrons), tufted deer (Elaphodus cephalophus), sika deer (Cervus nippon) and red deer (Cervus elaphus) (Cervinae). To track chromosomal rearrangements during Cervidae evolution, we summarized new data, and compared them with available X chromosomal maps and chromosome level assemblies of other species. We demonstrate the types of rearrangements that may have underlined the variability of Cervidae X chromosomes. We detected two types of cervine X chromosome-acrocentric and submetacentric. The acrocentric type is found in three independent deer lineages (subfamily Cervinae and in two Capreolinae tribes-Odocoileini and Capreolini). We show that chromosomal rearrangements on the X-chromosome in Cervidae occur at a higher frequency than in the entire Ruminantia lineage: the rate of rearrangements is 2 per 10 million years.

Selection, recombination and population history effects on runs of homozygosity (ROH) in wild red deer (Cervus elaphus)

The distribution of runs of homozygosity (ROH) may be shaped by a number of interacting processes such as selection, recombination and population history, but little is known about the importance of these mechanisms in shaping ROH in wild populations. We combined an empirical dataset of >3000 red deer genotyped at >35,000 genome-wide autosomal SNPs and evolutionary simulations to investigate the influence of each of these factors on ROH. We assessed ROH in a focal and comparison population to investigate the effect of population history. We investigated the role of recombination using both a physical map and a genetic linkage map to search for ROH. We found differences in ROH distribution between both populations and map types indicating that population history and local recombination rate have an effect on ROH. Finally, we ran forward genetic simulations with varying population histories, recombination rates and levels of selection, allowing us to further interpret our empirical data. These simulations showed that population history has a greater effect on ROH distribution than either recombination or selection. We further show that selection can cause genomic regions where ROH is common, only when the effective population size (Ne) is large or selection is particularly strong. In populations having undergone a population bottleneck, genetic drift can outweigh the effect of selection. Overall, we conclude that in this population, genetic drift resulting from a historical population bottleneck is most likely to have resulted in the observed ROH distribution, with selection possibly playing a minor role.

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.