Three decades of rat genomics: approaching the finish(ed) line.

The rat, Rattus norvegicus, has provided an important model for investigation of a range of characteristics of biomedical importance. Here we survey the origins of this species, its introduction into laboratory research and the emergence of genetic and genomic methods that utilize this model organism. Genomic studies have yielded important progress and provided new insight into several biologically important traits. However, some studies have been impeded by the lack of a complete and accurate reference genome for this species. New sequencing and genome assembly methods applied to the rat have resulted in a new reference genome assembly, GRCr8, which is a near telomere-to-telomere assembly of high base level accuracy that incorporates several elements not captured in prior assemblies. As genome assembly methods continue to advance and production costs become a less significant obstacle, genome assemblies for multiple inbred rat strains are emerging. These assemblies will allow a rat pangenome assembly to be constructed which captures all the genetic variation in strains selected for their utility in research and will overcome reference bias, a limitation associated with reliance on a single reference assembly. By this means, the full utility of this model organism to genomic studies will begin to be revealed.

mRatBN7.2: familiar and unfamiliar features of a new rat genome reference assembly.

Rat genomic tools have been slower to emerge than for those of humans and mice and have remained less thorough and comprehensive. The arrival of a new and improved rat reference genome, mRatBN7.2, in late 2020 is a welcome event. This assembly, like predecessor rat reference assemblies, is derived from an inbred Brown Norway rat. In this "user" survey we hope to provide other users of this assembly some insight into its characteristics and some assessment of its improvements as well as a few caveats that arise from the unique aspects of this assembly. mRatBN7.2 was generated by the Wellcome Sanger Institute as part of the large Vertebrate Genomes Project. This rat assembly has now joined human, mouse, chicken, and zebrafish in the National Center for Biotechnology Information (NCBI)'s Genome Reference Consortium, which provides ongoing curation of the assembly. Here we examine the technical procedures by which the assembly was created and assess how this assembly constitutes an improvement over its predecessor. We also indicate the technical limitations affecting the assembly, providing illustrations of how these limitations arise and the impact that results for this reference assembly.

Equine arteritis virus long-term persistence is orchestrated by CD8+ T lymphocyte transcription factors, inhibitory receptors, and the CXCL16/CXCR6 axis.

Equine arteritis virus (EAV) has the unique ability to establish long-term persistent infection in the reproductive tract of stallions and be sexually transmitted. Previous studies showed that long-term persistent infection is associated with a specific allele of the CXCL16 gene (CXCL16S) and that persistence is maintained despite the presence of local inflammatory and humoral and mucosal antibody responses. Here, we performed transcriptomic analysis of the ampullae, the primary site of EAV persistence in long-term EAV carrier stallions, to understand the molecular signatures of viral persistence. We demonstrated that the local CD8+ T lymphocyte response is predominantly orchestrated by the transcription factors eomesodermin (EOMES) and nuclear factor of activated T-cells cytoplasmic 2 (NFATC2), which is likely modulated by the upregulation of inhibitory receptors. Most importantly, EAV persistence is associated with an enhanced expression of CXCL16 and CXCR6 by infiltrating lymphocytes, providing evidence of the implication of this chemokine axis in the pathogenesis of persistent EAV infection in the stallion reproductive tract. Furthermore, we have established a link between the CXCL16 genotype and the gene expression profile in the ampullae of the stallion reproductive tract. Specifically, CXCL16 acts as a "hub" gene likely driving a specific transcriptional network. The findings herein are novel and strongly suggest that RNA viruses such as EAV could exploit the CXCL16/CXCR6 axis in order to modulate local inflammatory and immune responses in the male reproductive tract by inducing a dysfunctional CD8+ T lymphocyte response and unique lymphocyte homing in the reproductive tract.

Inflamed synovial fluid induces a homeostatic response in bone marrow mononuclear cells in vitro: Implications for joint therapy.

Synovial inflammation is a central feature of osteoarthritis (OA), elicited when local regulatory macrophages (M2-like) become overwhelmed, activating an inflammatory response (M1-like). Bone marrow mononuclear cells (BMNC) are a source of naïve macrophages capable of reducing joint inflammation and producing molecules essential for cartilage metabolism. This study investigated the response of BMNC to normal (SF) and inflamed synovial fluid (ISF). Equine BMNC cultured in autologous SF or ISF (n = 8 horses) developed into macrophage-rich cultures with phenotypes similar to cells native to normal SF and became more confluent in ISF (~100%) than SF (~25%). BMNC cultured in SF or ISF were neither M1- nor M2-like, but exhibited aspects of both phenotypes and a regulatory immune response, characterized by increasing counts of IL-10+ macrophages, decreasing IL-1ß concentrations and progressively increasing IL-10 and IGF-1 concentrations. Changes were more marked in ISF and suggest that homeostatic mechanisms were preserved over time and were potentially favored by progressive cell proliferation. Collectively, our data suggest that intra-articular BMNC could increase synovial macrophage counts, potentiating the macrophage- and IL-10-associated mechanisms of joint homeostasis lost during the progression of OA, preserving the production of cytokines involved in tissue repair (PGE2 , IL-10) generally impaired by frequently used corticosteroids.

Transcriptomic analysis of equine chorioallantois reveals immune networks and molecular mechanisms involved in nocardioform placentitis.

Nocardioform placentitis (NP) continues to result in episodic outbreaks of abortion and preterm birth in mares and remains a poorly understood disease. The objective of this study was to characterize the transcriptome of the chorioallantois (CA) of mares with NP. The CA were collected from mares with confirmed NP based upon histopathology, microbiological culture and PCR for Amycolatopsis spp. Samples were collected from the margin of the NP lesion (NPL, n = 4) and grossly normal region (NPN, n = 4). Additionally, CA samples were collected from normal postpartum mares (Control; CRL, n = 4). Transcriptome analysis identified 2892 differentially expressed genes (DEGs) in NPL vs. CRL and 2450 DEGs in NPL vs. NPN. Functional genomics analysis elucidated that inflammatory signaling, toll-like receptor signaling, inflammasome activation, chemotaxis, and apoptosis pathways are involved in NP. The increased leukocytic infiltration in NPL was associated with the upregulation of matrix metalloproteinase (MMP1, MMP3, and MMP8) and apoptosis-related genes, such as caspases (CASP3 and CASP7), which could explain placental separation associated with NP. Also, NP was associated with downregulation of several placenta-regulatory genes (ABCG2, GCM1, EPAS1, and NR3C1), angiogenesis-related genes (VEGFA, FLT1, KDR, and ANGPT2), and glucose transporter coding genes (GLUT1, GLUT10, and GLUT12), as well as upregulation of hypoxia-related genes (HIF1A and EGLN3), which could elucidate placental insufficiency accompanying NP. In conclusion, our findings revealed for the first time, the key regulators and mechanisms underlying placental inflammation, separation, and insufficiency during NP, which might lead to the development of efficacious therapies or diagnostic aids by targeting the key molecular pathways.

A Reference Genome Assembly of American Bison, Bison bison bison.

Bison are an icon of the American West and an ecologically, commercially, and culturally important species. Despite numbering in the hundreds of thousands today, conservation concerns remain for the species, including the impact on genetic diversity of a severe bottleneck around the turn of the 20th century and genetic introgression from domestic cattle. Genetic diversity and admixture are best evaluated at genome-wide scale, for which a high-quality reference is necessary. Here, we use trio binning of long reads from a bison-Simmental cattle (Bos taurus taurus) male F1 hybrid to sequence and assemble the genome of the American plains bison (Bison bison bison). The male haplotype genome is chromosome-scale, with a total length of 2.65 Gb across 775 scaffolds (839 contigs) and a scaffold N50 of 87.8 Mb. Our bison genome is ~13× more contiguous overall and ~3400× more contiguous at the contig level than the current bison reference genome. The bison genome sequence presented here (ARS-UCSC_bison1.0) will enable new research into the evolutionary history of this iconic megafauna species and provide a new tool for the management of bison populations in federal and commercial herds.

A Reference Genome Assembly of Simmental Cattle, Bos taurus taurus.

Genomics research has relied principally on the establishment and curation of a reference genome for the species. However, it is increasingly recognized that a single reference genome cannot fully describe the extent of genetic variation within many widely distributed species. Pangenome representations are based on high-quality genome assemblies of multiple individuals and intended to represent the broadest possible diversity within a species. A Bovine Pangenome Consortium (BPC) has recently been established to begin assembling genomes from more than 600 recognized breeds of cattle, together with other related species to provide information on ancestral alleles and haplotypes. Previously reported de novo genome assemblies for Angus, Brahman, Hereford, and Highland breeds of cattle are part of the initial BPC effort. The present report describes a complete single haplotype assembly at chromosome-scale for a fullblood Simmental cow from an F1 bison-cattle hybrid fetus by trio binning. Simmental cattle, also known as Fleckvieh due to their red and white spots, originated in central Europe in the 1830s as a triple-purpose breed selected for draught, meat, and dairy production. There are over 50 million Simmental cattle in the world, known today for their fast growth and beef yields. This assembly (ARS_Simm1.0) is similar in length to the other bovine assemblies at 2.86 Gb, with a scaffold N50 of 102 Mb (max scaffold 156.8 Mb) and meets or exceeds the continuity of the best Bos taurus reference assemblies to date.

Dynamic alteration in miRNA and mRNA expression profiles at different stages of chronic arsenic exposure-induced carcinogenesis in a human cell culture model of skin cancer.

Chronic arsenic exposure causes skin cancer, although the underlying molecular mechanisms are not well defined. Altered microRNA and mRNA expression likely play a pivotal role in carcinogenesis. Changes in genome-wide differential expression of miRNA and mRNA at 3 strategic time points upon chronic sodium arsenite (As3+) exposure were investigated in a well-validated HaCaT cell line model of arsenic-induced cutaneous squamous cell carcinoma (cSCC). Quadruplicate independent HaCaT cell cultures were exposed to 0 or 100 nM As3+ for up to 28-weeks (wk). Cell growth was monitored throughout the course of exposure and epithelial-mesenchymal transition (EMT) was examined employing immunoblot. Differentially expressed miRNA and mRNA profiles were generated at 7, 19, and 28-wk by RNA-seq, followed by identification of differentially expressed mRNA targets of differentially expressed miRNAs through expression pairing at each time point. Pathway analyses were performed for total differentially expressed mRNAs and for the miRNA targeted mRNAs at each time point. RNA-seq predictions were validated by immunoblot of selected target proteins. While the As3+-exposed cells grew slower initially, growth was equal to that of unexposed cells by 19-wk (transformation initiation), and exposed cells subsequently grew faster than passage-matched unexposed cells. As3+-exposed cells had undergone EMT at 28-wk. Pathway analyses demonstrate dysregulation of carcinogenesis-related pathways and networks in a complex coordinated manner at each time point. Immunoblot data largely corroborate RNA-seq predictions in the endoplasmic reticulum stress (ER stress) pathway. This study provides a detailed molecular picture of changes occurring during the arsenic-induced transformation of human keratinocytes.

MC1R and KIT Haplotypes Associate With Pigmentation Phenotypes of North American Yak (Bos grunniens).

Small numbers of domestic yak (Bos grunniens) were imported to North America in the late 19th century indirectly from the Qinghai-Tibetan Plateau. Coat color of yak is of interest for fiber production, aesthetics, and as a potential indicator of recent hybridization with cattle. North American yak are classified into 3 major coat color patterns depending upon the presence and extent of white markings. They are further classified by nose pigmentation (black or gray). The aim of this study was to identify loci involved in white patterning and nose pigmentation of North American yak. Genotyping by mass spectrometry of markers identified through Sanger and whole-genome sequencing revealed a 388 kb haplotype of KIT associated in a semi-dominant manner with white coloration in this population of yak. This KIT haplotype is similar to both a haplotype found in white-faced Chinese yak and to haplotypes found in cattle but is divergent from other Bos species such as bison, gaur, and banteng. Melanocortin 1 receptor (MC1R) was implicated as a dominant determinant of black nose color with a single haplotype containing 2 missense mutations perfectly associated with the phenotype. The MC1R haplotype associated with black nose pigment is also similar to cattle haplotypes. No cattle studied, however, shared either of the 2 haplotypes associated with color in yak, suggesting these alleles were introgressed into yak before they were imported to North America. These results provide molecular insight into the history of North American yak and information from which breeders can determine possible color outcomes of matings.

Genetics, Genomics, and Emergent Precision Medicine 12 Years After the Equine Reference Genome Was Published.

The first equine reference genome was completed in 2007 and published in 2009. This major accomplishment has enabled equine science to advance in ways that broadly parallel the transformative impact that genomics has had on many animal species including humans. A conceptual overview of reference genomes, genome annotation, and the major implications for equine science is presented. The relationship between genomic sequencing and the accelerating application of precision P4 medicine is discussed in the context of human and equine patients. Emergent technologies built on the foundation of genomic sequencing and rapidly gaining traction in research and clinical settings are introduced.

Downregulation of MicroRNA eca-mir-128 in Seminal Exosomes and Enhanced Expression of CXCL16 in the Stallion Reproductive Tract Are Associated with Long-Term Persistence of Equine Arteritis Virus.

Equine arteritis virus (EAV) can establish long-term persistent infection in the reproductive tract of stallions and is shed in the semen. Previous studies showed that long-term persistence is associated with a specific allele of the CXCL16 gene (CXCL16S) and that persistent infection is maintained despite the presence of a local inflammatory and humoral and mucosal antibody responses. In this study, we demonstrated that equine seminal exosomes (SEs) are enriched in a small subset of microRNAs (miRNAs). Most importantly, we demonstrated that long-term EAV persistence is associated with the downregulation of an SE-associated miRNA (eca-mir-128) and with an enhanced expression of CXCL16 in the reproductive tract, a putative target of eca-mir-128. The findings presented here suggest that SE eca-mir-128 is implicated in the regulation of the CXCL16/CXCR6 axis in the reproductive tract of persistently infected stallions, a chemokine axis strongly implicated in EAV persistence. This is a novel finding and warrants further investigation to identify its specific mechanism in modulating the CXCL16/CXCR6 axis in the reproductive tract of the EAV long-term carrier stallion.IMPORTANCE Equine arteritis virus (EAV) has the ability to establish long-term persistent infection in the stallion reproductive tract and to be shed in semen, which jeopardizes its worldwide control. Currently, the molecular mechanisms of viral persistence are being unraveled, and these are essential for the development of effective therapeutics to eliminate persistent infection. Recently, it has been determined that long-term persistence is associated with a specific allele of the CXCL16 gene (CXCL16S) and is maintained despite induction of local inflammatory, humoral, and mucosal antibody responses. This study demonstrated that long-term persistence is associated with the downregulation of seminal exosome miRNA eca-mir-128 and enhanced expression of its putative target, CXCL16, in the reproductive tract. For the first time, this study suggests complex interactions between eca-mir-128 and cellular elements at the site of EAV persistence and implicates this miRNA in the regulation of the CXCL16/CXCR6 axis in the reproductive tract during long-term persistence.

Equine Arteritis Virus Uses Equine CXCL16 as an Entry Receptor.

UNLABELLED: Previous studies in our laboratory have identified equine CXCL16 (EqCXCL16) to be a candidate molecule and possible cell entry receptor for equine arteritis virus (EAV). In horses, the CXCL16 gene is located on equine chromosome 11 (ECA11) and encodes a glycosylated, type I transmembrane protein with 247 amino acids. Stable transfection of HEK-293T cells with plasmid DNA carrying EqCXCL16 (HEK-EqCXCL16 cells) increased the proportion of the cell population permissive to EAV infection from <3% to almost 100%. The increase in permissiveness was blocked either by transfection of HEK-EqCXCL16 cells with small interfering RNAs (siRNAs) directed against EqCXCL16 or by pretreatment with guinea pig polyclonal antibody against EqCXCL16 protein (Gp anti-EqCXCL16 pAb). Furthermore, using a virus overlay protein-binding assay (VOPBA) in combination with far-Western blotting, gradient-purified EAV particles were shown to bind directly to the EqCXCL16 protein in vitro. The binding of biotinylated virulent EAV strain Bucyrus at 4°C was significantly higher in HEK-EqCXCL16 cells than nontransfected HEK-293T cells. Finally, the results demonstrated that EAV preferentially infects subpopulations of horse CD14(+) monocytes expressing EqCXCL16 and that infection of these cells is significantly reduced by pretreatment with Gp anti-EqCXCL16 pAb. The collective data from this study provide confirmatory evidence that the transmembrane form of EqCXCL16 likely plays a major role in EAV host cell entry processes, possibly acting as a primary receptor molecule for this virus. IMPORTANCE: Outbreaks of EVA can be a source of significant economic loss for the equine industry from high rates of abortion in pregnant mares, death in young foals, establishment of the carrier state in stallions, and trade restrictions imposed by various countries. Similar to other arteriviruses, EAV primarily targets cells of the monocyte/macrophage lineage, which, when infected, are believed to play a critical role in EVA pathogenesis. To this point, however, the host-specified molecules involved in EAV binding and entry into monocytes/macrophages have not been identified. Identification of the cellular receptors for EAV may provide insights to design antivirals and better prophylactic reagents. In this study, we have demonstrated that EqCXCL16 acts as an EAV entry receptor in EAV-susceptible cells, equine monocytes. These findings represent a significant advance in our understanding of the fundamental mechanisms associated with the entry of EAV into susceptible cells.

Genomic signatures of Mannheimia haemolytica that associate with the lungs of cattle with respiratory disease, an integrative conjugative element, and antibiotic resistance genes.

BACKGROUND: Mannheimia haemolytica typically resides in cattle as a commensal member of the upper respiratory tract microbiome. However, some strains can invade their lungs and cause respiratory disease and death, including those with multi-drug resistance. A nucleotide polymorphism typing system was developed for M. haemolytica from the genome sequences of 1133 North American isolates, and used to identify genetic differences between isolates from the lungs and upper respiratory tract of cattle with and without clinical signs of respiratory disease. RESULTS: A total of 26,081 nucleotide polymorphisms were characterized after quality control filtering of 48,403 putative polymorphisms. Phylogenetic analyses of nucleotide polymorphism genotypes split M. haemolytica into two major genotypes (1 and 2) that each were further divided into multiple subtypes. Multiple polymorphisms were identified with alleles that tagged genotypes 1 or 2, and their respective subtypes. Only genotype 2 M. haemolytica associated with the lungs of diseased cattle and the sequence of a particular integrative and conjugative element (ICE). Additionally, isolates belonging to one subtype of genotype 2 (2b), had the majority of antibiotic resistance genes detected in this study, which were assorted into seven combinations that ranged from 1 to 12 resistance genes. CONCLUSIONS: Typing of diverse M. haemolytica by nucleotide polymorphism genotypes successfully identified associations with diseased cattle lungs, ICE sequence, and antibiotic resistance genes. Management of cattle by their carriage of M. haemolytica could be an effective intervention strategy to reduce the prevalence of respiratory disease and supplemental needs for antibiotic treatments in North American herds.

Interaction between smoking history and gene expression levels impacts survival of breast cancer patients.

In contrast to studies focused on cigarette smoking and risk of breast cancer occurrence, this study explored the influence of smoking on breast cancer recurrence and progression. The goal was to evaluate the interaction between smoking history and gene expression levels on recurrence and overall survival of breast cancer patients. Multivariable Cox proportional hazards models were fitted for 48 cigarette smokers, 50 non-smokers, and the total population separately to determine which gene expressions and gene expression/cigarette usage interaction terms were significant in predicting overall and disease-free survival in breast cancer patients. Using methods similar to Andres et al. (BMC Cancer 13:326, 2013a; Horm Cancer 4:208-221, 2013b), multivariable analyses revealed CENPN, CETN1, CYP1A1, IRF2, LECT2, and NCOA1 to be important predictors for both breast carcinoma recurrence and mortality among smokers. Additionally, COMT was important for recurrence, and NAT1 and RIPK1 were important for mortality. In contrast, only IRF2, CETN1, and CYP1A1 were significant for disease recurrence and mortality among non-smokers, with NAT2 additionally significant for survival. Analysis of interaction between smoking status and gene expression values using the combined samples revealed significant interactions between smoking status and CYP1A1, LECT2, and CETN1. Signatures consisting of 7-8 genes were highly predictive for breast cancer recurrence and overall survival among smokers, with median C-index values of 0.8 and 0.73 for overall survival and recurrence, respectively. In contrast, median C-index values for non-smokers was only 0.59. Hence, significant interactions between gene expression and smoking status can play a key role in predicting breast cancer patient outcomes.

Reduced lentivirus susceptibility in sheep with TMEM154 mutations.

Visna/Maedi, or ovine progressive pneumonia (OPP) as it is known in the United States, is an incurable slow-acting disease of sheep caused by persistent lentivirus infection. This disease affects multiple tissues, including those of the respiratory and central nervous systems. Our aim was to identify ovine genetic risk factors for lentivirus infection. Sixty-nine matched pairs of infected cases and uninfected controls were identified among 736 naturally exposed sheep older than five years of age. These pairs were used in a genome-wide association study with 50,614 markers. A single SNP was identified in the ovine transmembrane protein (TMEM154) that exceeded genome-wide significance (unadjusted p-value 3×10(-9)). Sanger sequencing of the ovine TMEM154 coding region identified six missense and two frameshift deletion mutations in the predicted signal peptide and extracellular domain. Two TMEM154 haplotypes encoding glutamate (E) at position 35 were associated with infection while a third haplotype with lysine (K) at position 35 was not. Haplotypes encoding full-length E35 isoforms were analyzed together as genetic risk factors in a multi-breed, matched case-control design, with 61 pairs of 4-year-old ewes. The odds of infection for ewes with one copy of a full-length TMEM154 E35 allele were 28 times greater than the odds for those without (p-value<0.0001, 95% CI 5-1,100). In a combined analysis of nine cohorts with 2,705 sheep from Nebraska, Idaho, and Iowa, the relative risk of infection was 2.85 times greater for sheep with a full-length TMEM154 E35 allele (p-value<0.0001, 95% CI 2.36-3.43). Although rare, some sheep were homozygous for TMEM154 deletion mutations and remained uninfected despite a lifetime of significant exposure. Together, these findings indicate that TMEM154 may play a central role in ovine lentivirus infection and removing sheep with the most susceptible genotypes may help eradicate OPP and protect flocks from reinfection.

Thermodynamic Details of Pinholin S2168 Activation Revealed Using Alchemical Free Energy Simulations.

Pinholin S2168 is a viral integral membrane protein whose function is to form nanoscopic "pinholes" in bacterial cell membranes to induce cell lysis as part of the viral replication cycle. Pinholin can transition from an inactive to an active conformation by exposing a transmembrane domain (TMD1) to the extracellular fluid. Upon activation, several copies of the protein assemble via interactions among a second transmembrane domain (TMD2) to form a single pore, thus hastening cell lysis and viral escape. The following experiments provide conformational descriptors of pinholin in active and inactive states and elucidate the molecular driving forces that control pinholin activity. In the present study, molecular dynamics (MD) simulations have been used to refine experimentally derived conformational descriptors into an atomistically detailed model of irsS2168, an antiholin mutant. To provide additional details about the thermodynamics of pinholin activation and to overcome large intrinsic kinetic barriers to activation, alchemical free energy simulations have been conducted. Alchemical mutations reveal the change in folding free energy upon mutation. The results suggest that alchemical mutations are an effective tool to rationalize experimental observations and predict the effects of site mutations on conformational states for proteins integrated into lipid bilayers. S16F, A17Q, A17Q+G21Q, and A17Q+G21Q+G14Q mutants reveal how changes in hydrophilicity and disruption of the glycine zipper motif influence pinholin's thermodynamic equilibrium, favoring the active conformation. These findings align with experimental observations from DEER spectroscopy, demonstrating that mutations increasing the hydrophilicity of TMD1 promote activation by making TMD1 more likely to exit the membrane and enter the extracellular fluid.

The Ruminant Telomere-to-Telomere (RT2T) Consortium.

Telomere-to-telomere (T2T) assemblies reveal new insights into the structure and function of the previously 'invisible' parts of the genome and allow comparative analyses of complete genomes across entire clades. We present here an open collaborative effort, termed the 'Ruminant T2T Consortium' (RT2T), that aims to generate complete diploid assemblies for numerous species of the Artiodactyla suborder Ruminantia to examine chromosomal evolution in the context of natural selection and domestication of species used as livestock.

Trio-binning of a hinny refines the comparative organization of the horse and donkey X chromosomes and reveals novel species-specific features.

We generated single haplotype assemblies from a hinny hybrid which significantly improved the gapless contiguity for horse and donkey autosomal genomes and the X chromosomes. We added over 15 Mb of missing sequence to both X chromosomes, 60 Mb to donkey autosomes and corrected numerous errors in donkey and some in horse reference genomes. We resolved functionally important X-linked repeats: the DXZ4 macrosatellite and ampliconic Equine Testis Specific Transcript Y7 (ETSTY7). We pinpointed the location of the pseudoautosomal boundaries (PAB) and determined the size of the horse (1.8 Mb) and donkey (1.88 Mb) pseudoautosomal regions (PARs). We discovered distinct differences in horse and donkey PABs: a testis-expressed gene, XKR3Y, spans horse PAB with exons1-2 located in Y and exon3 in the X-Y PAR, whereas the donkey XKR3Y is Y-specific. DXZ4 had a similar ~ 8 kb monomer in both species with 10 copies in horse and 20 in donkey. We assigned hundreds of copies of ETSTY7, a sequence horizontally transferred from Parascaris and massively amplified in equids, to horse and donkey X chromosomes and three autosomes. The findings and products contribute to molecular studies of equid biology and advance research on X-linked conditions, sex chromosome regulation and evolution in equids.

First gene-edited calf with reduced susceptibility to a major viral pathogen.

Bovine viral diarrhea virus (BVDV) is one of the most important viruses affecting the health and well-being of bovine species throughout the world. Here, we used CRISPR-mediated homology-directed repair and somatic cell nuclear transfer to produce a live calf with a six amino acid substitution in the BVDV binding domain of bovine CD46. The result was a gene-edited calf with dramatically reduced susceptibility to infection as measured by reduced clinical signs and the lack of viral infection in white blood cells. The edited calf has no off-target edits and appears normal and healthy at 20 months of age without obvious adverse effects from the on-target edit. This precision bred, proof-of-concept animal provides the first evidence that intentional genome alterations in the CD46 gene may reduce the burden of BVDV-associated diseases in cattle and is consistent with our stepwise, in vitro and ex vivo experiments with cell lines and matched fetal clones.