Short communication: Identification of the pseudoautosomal region in the Hereford bovine reference genome assembly ARS-UCD1.2

Genomic regions and biological pathways associated with sex-limited reproductive traits in bovine species

Many animal species exhibit sex-limited traits, where certain phenotypes are exclusively expressed in one sex. Yet, the genomic regions that contribute to these sex-limited traits in males and females remain a subject of debate. Reproductive traits are ideal phenotypes to study sexual differences since they are mostly expressed in a sex-limited way. Therefore, this study aims to use local correlation analyses to identify genomic regions and biological pathways significantly associated with male and female sex-limited traits in two distinct cattle breeds (Brahman [BB] and Tropical Composite [TC]). We used the Correlation Scan method to perform local correlation analysis on 42 trait pairs consisting of six female and seven male reproductive traits recorded on ~1,000 animals for each sex in each breed. To pinpoint a specific region associated with these sex-limited reproductive traits, we investigated the genomic region(s) consistently identified as significant across the 42 trait pairs in each breed. The genes found in the identified regions were subjected to Quantitative Trait Loci (QTL) colocalization, QTL enrichment analyses, and functional analyses to gain biological insight into sexual differences. We found that the genomic regions associated with the sex-limited reproductive phenotypes are widely distributed across all the chromosomes. However, no single region across the genome was associated with all the 42 reproductive trait pairs in the two breeds. Nevertheless, we found a region on the X-chromosome to be most significant for 80% to 90% (BB: 33 and TC: 38) of the total 42 trait pairs. A considerable number of the genes in this region were regulatory genes. By considering only genomic regions that were significant for at least 50% of the 42 trait pairs, we observed more regions spread across the autosomes and the X-chromosome. All genomic regions identified were highly enriched for trait-specific QTL linked to sex-limited traits (percentage of normal sperm, metabolic weight, average daily gain, carcass weight, age at puberty, etc.). The gene list created from these identified regions was enriched for biological pathways that contribute to the observed differences between sexes. Our results demonstrate that genomic regions associated with male and female sex-limited reproductive traits are distributed across the genome. Yet, chromosome X seems to exert a relatively larger effect on the phenotypic variation observed between the sexes.

Identification of candidate novel production variants on the Bos taurus chromosome X

Chromosome X is often excluded from bovine genetic studies due to complications caused by the sex specific nature of the chromosome. As chromosome X is the second largest cattle chromosome and makes up approximately 6% of the female genome, finding ways to include chromosome X in dairy genetic studies is important. Using female animals and treating chromosome X as an autosome, we performed X chromosome inclusive genome-wide association studies in the selective breeding environment of the New Zealand dairy industry, aiming to identify chromosome X variants associated with milk production traits. We report on the findings of these genome-wide association studies and their potential effect within the dairy industry. We identify missense mutations in the MOSPD1 and CCDC160 genes that are associated with decreased milk volume and protein production and increased fat production. Both of these mutations are exonic SNP that are more prevalent in the Jersey breed than in Holstein-Friesians. Of the 2 candidates proposed it is likely that only one is causal, though we have not been able to identify which is more likely.

X-linked genes influence various complex traits in dairy cattle

BACKGROUND

The search for quantitative trait loci (QTL) affecting traits of interest in mammals is frequently limited to autosomes, with the X chromosome excluded because of its hemizygosity in males. This study aimed to assess the importance of the X chromosome in the genetic determinism of 11 complex traits related to milk production, milk composition, mastitis resistance, fertility, and stature in 236,496 cows from three major French dairy breeds (Holstein, Montbéliarde, and Normande) and three breeds of regional importance (Abondance, Tarentaise, and Vosgienne).

RESULTS

Estimates of the proportions of heritability due to autosomes and X chromosome (h²_X) were consistent among breeds. On average over the 11 traits, h²_X=0.008 and the X chromosome explained ~ 3.5% of total genetic variance. GWAS was performed within-breed at the sequence level (~ 200,000 genetic variants) and then combined in a meta-analysis. QTL were identified for most breeds and traits analyzed, with the exception of Tarentaise and Vosgienne and two fertility traits. Overall, 3, 74, 59, and 71 QTL were identified in Abondance, Montbéliarde, Normande, and Holstein, respectively, and most were associated with the most-heritable traits (milk traits and stature). The meta-analyses, which assessed a total of 157 QTL for the different traits, highlighted new QTL and refined the positions of some QTL found in the within-breed analyses. Altogether, our analyses identified a number of functional candidate genes, with the most notable being GPC3, MBNL3, HS6ST2, and DMD for dairy traits; TMEM164, ACSL4, ENOX2, HTR2C, AMOT, and IRAK1 for udder health; MAMLD1 and COL4A6 for fertility; and NRK, ESX1, GPR50, GPC3, and GPC4 for stature.

CONCLUSIONS

This study demonstrates the importance of the X chromosome in the genetic determinism of complex traits in dairy cattle and highlights new functional candidate genes and variants for these traits. These results could potentially be extended to other species as many X-linked genes are shared among mammals.

A potent allele marker related to low bull conception rate in Japanese Black bulls

Over the years, there has been considerable variation in the bull conception rate (BCR) of Japanese Black cattle; moreover, several Japanese Black bulls with a low BCR of ≤10% have been identified. However, the alleles responsible for the low BCR are not determined yet. Therefore, in this study, we aimed to identify single-nucleotide polymorphisms (SNPs) for predicting low BCR. To this end, the genome of Japanese Black bulls was comprehensively examined by a genome-wide association study with whole-exome sequencing (WES), and the effect of the identified marker regions on BCR was determined. The WES analysis of six sub-fertile bulls with a BCR of ≤10% and 73 normal bulls with a BCR of ≥40% identified a homozygous genotype for low BCR in Bos taurus autosome 5 in the region between 116.2 and 117.9 Mb. The g.116408653G > A SNP in this region had the most significant effect on the BCR (P-value = 1.0 × 10^-23), and the GG (55.4 ± 11.2%) and AG (54.4 ± 9.4%) genotypes in the SNP had a higher phenotype than the AA (9.5 ± 6.1%) genotype for the BCR. The mixed model analysis revealed that g.116408653G > A was related to approximately 43% of the total genetic variance. In conclusion, the AA genotype of g.116408653G > A is a useful index for identifying sub-fertile Japanese Black bulls. Some positive and negative effects of SNP on the BCR were presumed to identify the causative mutations, which can help evaluate bull fertility.

Deviations from Mendelian Inheritance on Bovine X-Chromosome Revealing Recombination, Sex-of-Offspring Effects and Fertility-Related Candidate Genes

Transmission ratio distortion (TRD), or significant deviations from Mendelian inheritance, is a well-studied phenomenon on autosomal chromosomes, but has not yet received attention on sex chromosomes. TRD was analyzed on 3832 heterosomal single nucleotide polymorphisms (SNPs) and 400 pseudoautosomal SNPs spanning the length of the X-chromosome using 436,651 genotyped Holstein cattle. On the pseudoautosomal region, an opposite sire-TRD pattern between male and female offspring was identified for 149 SNPs. This finding revealed unique SNPs linked to a specific-sex (Y- or X-) chromosome and describes the accumulation of recombination events across the pseudoautosomal region. On the heterosomal region, 13 SNPs and 69 haplotype windows were identified with dam-TRD. Functional analyses for TRD regions highlighted relevant biological functions responsible to regulate spermatogenesis, development of Sertoli cells, homeostasis of endometrium tissue and embryonic development. This study uncovered the prevalence of different TRD patterns across both heterosomal and pseudoautosomal regions of the X-chromosome and revealed functional candidate genes for bovine reproduction.

Genetic Modeling and Genomic Analyses of Yearling Temperament in American Angus Cattle and Its Relationship With Productive Efficiency and Resilience Traits

Cattle temperament has been considered by farmers as a key breeding goal due to its relevance for cattlemen's safety, animal welfare, resilience, and longevity and its association with many economically important traits (e.g., production and meat quality). The definition of proper statistical models, accurate variance component estimates, and knowledge on the genetic background of the indicator trait evaluated are of great importance for accurately predicting the genetic merit of breeding animals. Therefore, 266,029 American Angus cattle with yearling temperament records (1-6 score) were used to evaluate statistical models and estimate variance components; investigate the association of sex and farm management with temperament; assess the weighted correlation of estimated breeding values for temperament and productive, reproductive efficiency and resilience traits; and perform a weighted single-step genome-wide association analysis using 69,559 animals genotyped for 54,609 single-nucleotide polymorphisms. Sex and extrinsic factors were significantly associated with temperament, including conception type, age of dam, birth season, and additional animal-human interactions. Similar results were observed among models including only the direct additive genetic effect and when adding other maternal effects. Estimated heritability of temperament was equal to 0.39 on the liability scale. Favorable genetic correlations were observed between temperament and other relevant traits, including growth, feed efficiency, meat quality, and reproductive traits. The highest approximated genetic correlations were observed between temperament and growth traits (weaning weight, 0.28; yearling weight, 0.28). Altogether, we identified 11 genomic regions, located across nine chromosomes including BTAX, explaining 3.33% of the total additive genetic variance. The candidate genes identified were enriched in pathways related to vision, which could be associated with reception of stimulus and/or cognitive abilities. This study encompasses large and diverse phenotypic, genomic, and pedigree datasets of US Angus cattle. Yearling temperament is a highly heritable and polygenic trait that can be improved through genetic selection. Direct selection for temperament is not expected to result in unfavorable responses on other relevant traits due to the favorable or low genetic correlations observed. In summary, this study contributes to a better understanding of the impact of maternal effects, extrinsic factors, and various genomic regions associated with yearling temperament in North American Angus cattle.

Genome-wide detection of non-additive quantitative trait loci for semen production traits in beef and dairy bulls

Semen production traits are important aspects of bull fertility, because semen quantity leads to direct profits for artificial insemination centres, and semen quality is associated with the probability of achieving a pregnancy. Most genome-wide association studies (GWASs) for semen production traits have assumed that each quantitative trait locus (QTL) has an additive effect. However, GWASs that account for non-additive effects are also important in fitness traits, such as bull fertility. Here, we performed a GWAS using models that accounted for additive and non-additive effects to evaluate the importance of non-additive effects on five semen production traits in beef and dairy bulls. A total of 65 463 records for 615 Japanese Black bulls (JB) and 50 734 records for 873 Holstein bulls (HOL), which were previously genotyped using the Illumina BovineSNP50 BeadChip, were used to estimate genetic parameters and perform GWAS. The heritability estimates were low (ranged from 0.11 to 0.23), and the repeatability estimates were low to moderate (ranged from 0.28 to 0.45) in both breeds. The estimated repeatability was approximately twice as high as the estimated heritability for all traits. In this study, only one significant region with an additive effect was detected in each breed, but multiple significant regions with non-additive effects were detected for each breed. In particular, the region at approximately 64 Mbp on Bos taurus autosome 17 had the highest significant non-additive effect on four semen production traits in HOL. The rs41843851 single nucleotide polymorphism (SNP) in the region had a much lower P-value for the non-additive effect (P-value = 1.1 × 10^-31) than for the additive effect (P-value = 1.1 × 10^-8) in sperm motility. The AA and AB genotypes on the SNP had a higher phenotype than the BB genotype in HOL, and there was no bull with the BB genotype in JB. Our results showed that non-additive QTLs affect semen production traits, and a novel QTL accounting for non-additive effects could be detected by GWAS. This study provides new insights into non-additive QTLs that affect fitness traits, such as semen production traits in beef and dairy bulls.

Haplotype-Based Single-Step GWAS for Yearling Temperament in American Angus Cattle

Behavior is a complex trait and, therefore, understanding its genetic architecture is paramount for the development of effective breeding strategies. The objective of this study was to perform traditional and weighted single-step genome-wide association studies (ssGWAS and WssGWAS, respectively) for yearling temperament (YT) in North American Angus cattle using haplotypes. Approximately 266 K YT records and 70 K animals genotyped using a 50 K single nucleotide polymorphisms (SNP) panel were used. Linkage disequilibrium thresholds (LD) of 0.15, 0.50, and 0.80 were used to create the haploblocks, and the inclusion of non-LD-clustered SNPs (NCSNP) with the haplotypes in the genomic models was also evaluated. WssGWAS did not perform better than ssGWAS. Cattle YT was found to be a highly polygenic trait, with genes and quantitative trait loci (QTL) broadly distributed across the whole genome. Association studies using LD-based haplotypes should include NCSNPs and different LD thresholds to increase the likelihood of finding the relevant genomic regions affecting the trait of interest. The main candidate genes identified, i.e., ATXN10, ADAM10, VAX2, ATP6V1B1, CRISPLD1, CAPRIN1, FA2H, SPEF2, PLXNA1, and CACNA2D3, are involved in important biological processes and metabolic pathways related to behavioral traits, social interactions, and aggressiveness in cattle. Future studies should further investigate the role of these candidate genes.

Design and performance of a bovine 200 k SNP chip developed for endangered German Black Pied cattle (DSN)

BACKGROUND

German Black Pied cattle (DSN) are an endangered dual-purpose breed which was largely replaced by Holstein cattle due to their lower milk yield. DSN cattle are kept as a genetic reserve with a current herd size of around 2500 animals. The ability to track sequence variants specific to DSN could help to support the conservation of DSN's genetic diversity and to provide avenues for genetic improvement.

RESULTS

Whole-genome sequencing data of 304 DSN cattle were used to design a customized DSN200k SNP chip harboring 182,154 variants (173,569 SNPs and 8585 indels) based on ten selection categories. We included variants of interest to DSN such as DSN unique variants and variants from previous association studies in DSN, but also variants of general interest such as variants with predicted consequences of high, moderate, or low impact on the transcripts and SNPs from the Illumina BovineSNP50 BeadChip. Further, the selection of variants based on haplotype blocks ensured that the whole-genome was uniformly covered with an average variant distance of 14.4 kb on autosomes. Using 300 DSN and 162 animals from other cattle breeds including Holstein, endangered local cattle populations, and also a Bos indicus breed, performance of the SNP chip was evaluated. Altogether, 171,978 (94.31%) of the variants were successfully called in at least one of the analyzed breeds. In DSN, the number of successfully called variants was 166,563 (91.44%) while 156,684 (86.02%) were segregating at a minor allele frequency > 1%. The concordance rate between technical replicates was 99.83 ± 0.19%.

CONCLUSION

The DSN200k SNP chip was proved useful for DSN and other Bos taurus as well as one Bos indicus breed. It is suitable for genetic diversity management and marker-assisted selection of DSN animals. Moreover, variants that were segregating in other breeds can be used for the design of breed-specific customized SNP chips. This will be of great value in the application of conservation programs for endangered local populations in the future.

Genetic and Genomic Analyses of Service Sire Effect on Female Reproductive Traits in Holstein Cattle

Fertility and reproductive performance are key drivers of dairy farm profitability. Hence, reproduction traits have been included in a large majority of worldwide dairy cattle selection indexes. The reproductive traits are lowly heritable but can be improved through direct genetic selection. However, most scientific studies and dairy cattle breeding programs have focused solely on the genetic effects of the dam (GED) on reproductive performance and, therefore, ignored the contribution of the service sire in the phenotypic outcomes. This study aimed to investigate the service sire effects on female reproductive traits in Holstein cattle from a genomic perspective. Genetic parameter estimation and genome-wide association studies (GWAS) were performed for the genetic effect of service sire (GESS) on conception rate (CR), 56-day non-return rate (NRR56), calving ease (CE), stillbirth (SB), and gestation length (GL). Our findings indicate that the additive genetic effects of both sire and dam contribute to the phenotypic variance of reproductive traits measured in females (0.0196 vs. 0.0109, 0.0237 vs. 0.0133, 0.0040 vs. 0.0289, 0.0782 vs. 0.0083, and 0.1024 vs. 0.1020 for GESS and GED heritability estimates for CR, NRR56, CE, SB, and GL, respectively), and these two genetic effects are positively correlated for SB (0.1394) and GL (0.7871). Interestingly, the breeding values for GESS on insemination success traits (CR and NRR56) are unfavorably and significantly correlated with some production, health, and type breeding values (ranging from -0.449 to 0.274), while the GESS values on calving traits (CE, SB, and GL) are usually favorably associated with those traits (ranging from -0.493 to 0.313). One hundred sixty-two significant single-nucleotide polymorphisms (SNPs) and their surrounding protein-coding genes were identified as significantly associated with GESS and GED, respectively. Six genes overlapped between GESS and GED for calving traits and 10 genes overlapped between GESS for success traits and calving traits. Our findings indicate the importance of considering the GESS when genetically evaluating the female reproductive traits in Holstein cattle.

Effect of genomic X-chromosome regions on Nelore bull fertility

Scrotal circumference (SC) is a commonly used trait related to sexual precocity in bulls. Genome-wide association studies have uncovered a lot of genes related to this trait, however, only those present on autosomes. The inclusion of the second biggest chromosome (BTAX) can improve the knowledge of the genetic architecture of this trait. In this study, we performed a weighted, single-step, genome-wide association study using a 777 k BovineHD BeadChip (IllumHD) to analyze the association between SNPs and SC in Brazilian Nelore cattle. Phenotypes from 79,300 males and 3263 genotypes (2017 from females and 1246 from males)-(39,367 SNPs markers located at ChrX) were used. We identified eight regions on chromosome X that displayed important associations with SC. The results showed that together the genomic windows explained 28.52% of the genetic variance for the examined trait. Genes with potential functions in reproduction and fertility regulation were highlighted as candidates for sexual precocity rates in Nelore cattle (AFF2 and PJA1). Moreover, we found 10 genes that had not previously been identified as being associated with sexual precocity traits in cattle. These findings will further advance our understanding of the genetic architecture, considering mainly the presence of the chromosome X, for indicine cattle reproductive traits, being useful in the context of genomic prediction in beef cattle.

Genetic parameter estimates for male and female fertility traits using genomic data to improve fertility in Australian beef cattle

Context Studies have shown that favourable genetic correlations exist between female and male fertility traits. However, investigations regarding these correlations in Australian tropical beef cattle are limited to either pedigree or single-breed analysis. Aim The study aims to use genomic information to estimate genetic parameters of six female and seven male fertility traits measured during the first 2 years of life, in two tropical breeds. Methods Single-, bivariate and multi-trait models were used to analyse fertility data from Brahman (BB; 996 cows and 1022 bulls); and Tropical Composite (TC; 1091 cows and 998 bulls) cattle genotyped with high-density single-nucleotide polymorphism chip assay. Key results Heritability estimates in BB cows ranged from low (0.07 ± 0.04) for days to calving at the first calving opportunity (DC1, days) to high (0.57 ± 0.08) for age at first corpus luteum (AGECL, days). In BB bulls, estimates varied from low (0.09 ± 0.05) for sperm motility (score 1–5) to high (0.64 ± 0.06) for scrotal circumference (SC) measured at 24 months (SC24, cm). Similarly, heritability estimates in TC cows were low (0.04 ± 0.03) for DC1 and high (0.69 ± 0.02) for AGECL. In TC bulls, the heritability was low (0.09 ± 0.05) for sperm motility and high (0.69 ± 0.07) for SC24. Within-sex for both breeds, blood concentrations of insulin growth-factor 1 (IGF1) measured in cows at 18 months (IGF1c) were negatively correlated with female fertility phenotypes. In BB, across-sex, bulls’ blood concentration of IGF1 measured at 6 months (IGF1b) was a good indicator trait for the following four female traits: AGECL, the first postpartum anoestrus interval, age at first calving and DC1. In TC, IGF1b and percentage normal sperm were good predictors of female fertility phenotypes. Conclusions The heritability estimates and genomic correlations from the present study generally support and confirmed the earlier estimates from pedigree analyses. The findings suggest that selection for female fertility traits will benefit male fertility, and vice versa. Implications Heritability estimates and genomic correlations suggest that we can select for fertility traits measured early in life, with benefits within and across sex. Using traits available through veterinary assessment of bull fertility as selection indicators will enhance bull and cow fertility, which can lead to better breeding rates in tropical herds.

Theoretical and empirical comparisons of expected and realized relationships for the X-chromosome

Background

X-chromosomal loci present different inheritance patterns compared to autosomal loci and must be modeled accordingly. Sexual chromosomes are not systematically considered in whole-genome relationship matrices although rules based on genealogical or marker information have been derived. Loci on the X-chromosome could have a significant contribution to the additive genetic variance, in particular for some traits such as those related to reproduction. Thus, accounting for the X-chromosome relationship matrix might be informative to better understand the architecture of complex traits (e.g., by estimating the variance associated to this chromosome) and to improve their genomic prediction. For such applications, previous studies have shown the benefits of combining information from genotyped and ungenotyped individuals.

Results

In this paper, we start by presenting rules to compute a genomic relationship matrix (GRM) for the X-chromosome (G^X) without making any assumption on dosage compensation, and based on coding of gene content with 0/1 for males and 0/1/2 for females. This coding adjusts naturally to previously derived pedigree-based relationships (S) for the X-chromosome. When needed, we propose to accommodate and estimate dosage compensation and genetic heterogeneity across sexes via multiple trait models. Using a Holstein dairy cattle dataset, including males and females, we then empirically illustrate that realized relationships (G^X) matches expectations (S). However, G^X presents high deviations from S. G^X has also a lower dimensionality compared to the autosomal GRM. In particular, individuals are frequently identical along the entire chromosome. Finally, we confirm that the heritability of gene content for markers on the X-chromosome that are estimated by using S is 1, further demonstrating that S and G^X can be combined. For the pseudo-autosomal region, we demonstrate that the expected relationships vary according to position because of the sex-gradient. We end by presenting the rules to construct the 'H matrix’ by combining both relationship matrices.

Conclusions

This work shows theoretically and empirically that a pedigree-based relationship matrix built with rules specifically developed for the X-chromosome (S) matches the realized GRM for the X-chromosome. Therefore, applications that combine expected relationships and genotypes for markers on the X-chromosome should use S and G^X.

X chromosome variants are associated with male fertility traits in two bovine populations

BACKGROUND

Twenty-five phenotypes were measured as indicators of bull fertility (1099 Brahman and 1719 Tropical Composite bulls). Measurements included sperm morphology, scrotal circumference, and sperm chromatin phenotypes such as DNA fragmentation and protamine deficiency. We estimated the heritability of these phenotypes and carried out genome-wide association studies (GWAS) within breed, using the bovine high-density chip, to detect quantitative trait loci (QTL).

RESULTS

Our analyses suggested that both sperm DNA fragmentation and sperm protamine deficiency are heritable (h2 from 0.10 to 0.22). To confirm these first estimates of heritability, further studies on sperm chromatin traits, with larger datasets are necessary. Our GWAS identified 12 QTL for bull fertility traits, based on at least five polymorphisms (P < 10-8) for each QTL. Five QTL were identified in Brahman and another seven in Tropical Composite bulls. Most of the significant polymorphisms detected in both breeds and nine of the 12 QTL were on chromosome X. The QTL were breed-specific, but for some traits, a closer inspection of the GWAS results revealed suggestive single nucleotide polymorphism (SNP) associations (P < 10-7) in both breeds. For example, the QTL for inhibin level in Braham could be relevant to Tropical Composites too (many polymorphisms reached P < 10-7 in the same region). The QTL for sperm midpiece morphological abnormalities on chromosome X (QTL peak at 4.92 Mb, P < 10-17) is an example of a breed-specific QTL, supported by 143 significant SNPs (P < 10-8) in Brahman, but absent in Tropical Composites. Our GWAS results add evidence to the mammalian specialization of the X chromosome, which during evolution has accumulated genes linked to spermatogenesis. Some of the polymorphisms on chromosome X were associated to more than one genetically correlated trait (correlations ranged from 0.33 to 0.51). Correlations and shared polymorphism associations support the hypothesis that these phenotypes share the same underlying cause, i.e. defective spermatogenesis.

CONCLUSIONS

Genetic improvement for bull fertility is possible through genomic selection, which is likely more accurate if the QTL on chromosome X are considered in the predictions. Polymorphisms associated with male fertility accumulate on this chromosome in cattle, as in humans and mice, suggesting its specialization.

Genetic architecture of individual variation in recombination rate on the X chromosome in cattle

Meiotic recombination is an essential biological process that ensures proper chromosome segregation and creates genetic diversity. Individual variation in global recombination rates has been shown to be heritable in several species, and variants significantly associated with this trait have been identified. Recombination on the sex chromosome has often been ignored in these studies although this trait may be particularly interesting as it may correspond to a biological process distinct from that on autosomes. For instance, recombination in males is restricted to the pseudo-autosomal region (PAR). We herein used a large cattle pedigree with more than 100,000 genotyped animals to improve the genetic map of the X chromosome and to study the genetic architecture of individual variation in recombination rate on the sex chromosome (XRR). The length of the genetic map was 46.4 and 121.2 cM in males and females, respectively, but the recombination rate in the PAR was six times higher in males. The heritability of CO counts on the X chromosome was comparable to that of autosomes in males (0.011) but larger than that of autosomes in females (0.024). XRR was highly correlated (0.76) with global recombination rate (GRR) in females, suggesting that both traits might be governed by shared variants. In agreement, a set of eleven previously identified variants associated with GRR had correlated effects on female XRR (0.86). In males, XRR and GRR appeared to be distinct traits, although more accurate CO counts on the PAR would be valuable to confirm these results.

Symposium review: Advances in sequencing technology herald a new frontier in cattle genomics and genome-enabled selection

The cattle reference genome assembly has underpinned major innovations in beef and dairy genetics through genome-enabled selection, including removal of deleterious recessive variants and selection for favorable alleles affecting quantitative production traits. The initial reference assemblies, up to and including UMD3.1 and Btau4.1, were based on a combination of clone-by-clone sequencing of bacterial artificial chromosome clones generated from blood DNA of a Hereford bull and whole-genome shotgun sequencing of blood DNA from his inbred daughter/granddaughter named L1 Dominette 01449 (Dominette). The approach introduced assembly gaps, misassemblies, and errors, and it limited the ability to assemble regions that undergo rearrangement in blood cells, such as immune gene clusters. Nonetheless, the reference supported the creation of genotyping tools and provided a basis for many studies of gene expression. Recently, long-read sequencing technologies have emerged that facilitated a re-assembly of the reference genome, using lung tissue from Dominette to resolve many of the problems and providing a bridge to place historical studies in common context. The new reference, ARS-UCD1.2, successfully assembled germline immune gene clusters and improved overall continuity (i.e., reduction of gaps and inversions) by over 250-fold. This reference properly places nearly all of the legacy genetic markers used for over a decade in the industry. In this review, we discuss the improvements made to the cattle reference; remaining issues present in the assembly; tools developed to support genome-based studies in beef and dairy cattle; and the emergence of newer genome assembly methods that are producing even higher-quality assemblies for other breeds of cattle at a fraction of the cost. The new frontier for cattle genomics research will likely include a transition from the individual Hereford reference genome, to a "pan-genome" reference, representing all the DNA segments existing in commonly used cattle breeds, bringing the cattle reference into line with the current direction of human genome research.

New insights into mammalian sex chromosome structure and evolution using high-quality sequences from bovine X and Y chromosomes

BACKGROUND

Mammalian X chromosomes are mainly euchromatic with a similar size and structure among species whereas Y chromosomes are smaller, have undergone substantial evolutionary changes and accumulated male specific genes and genes involved in sex determination. The pseudoautosomal region (PAR) is conserved on the X and Y and pair during meiosis. The structure, evolution and function of mammalian sex chromosomes, particularly the Y chromsome, is still poorly understood because few species have high quality sex chromosome assemblies.

RESULTS

Here we report the first bovine sex chromosome assemblies that include the complete PAR spanning 6.84 Mb and three Y chromosome X-degenerate (X-d) regions. The PAR comprises 31 genes, including genes that are missing from the X chromosome in current cattle, sheep and goat reference genomes. Twenty-nine PAR genes are single-copy genes and two are multi-copy gene families, OBP, which has 3 copies and BDA20, which has 4 copies. The Y chromosome X-d1, 2a and 2b regions contain 11, 2 and 2 gametologs, respectively.

CONCLUSIONS

The ruminant PAR comprises 31 genes and is similar to the PAR of pig and dog but extends further than those of human and horse. Differences in the pseudoautosomal boundaries are consistent with evolutionary divergence times. A bovidae-specific expansion of members of the lipocalin gene family in the PAR reported here, may affect immune-modulation and anti-inflammatory responses in ruminants. Comparison of the X-d regions of Y chromosomes across species revealed that five of the X-Y gametologs, which are known to be global regulators of gene activity and candidate sexual dimorphism genes, are conserved.