Efficient correction of a deleterious point mutation in primary horse fibroblasts with CRISPR-Cas9

Phenotypic selection during animal domestication has resulted in unwanted incorporation of deleterious mutations. In horses, the autosomal recessive condition known as Glycogen Branching Enzyme Deficiency (GBED) is the result of one of these deleterious mutations (102C > A), in the first exon of the GBE1 gene (GBE1102C>A). With recent advances in genome editing, this type of genetic mutation can be precisely repaired. In this study, we used the RNA-guided nuclease CRISPR-Cas9 (clustered regularly-interspaced short palindromic repeats/CRISPR-associated protein 9) to correct the GBE1102C>A mutation in a primary fibroblast cell line derived from a high genetic merit heterozygous stallion. To correct this mutation by homologous recombination (HR), we designed a series of single guide RNAs (sgRNAs) flanking the mutation and provided different single-stranded donor DNA templates. The distance between the Cas9-mediated double-stranded break (DSB) to the mutation site, rather than DSB efficiency, was the primary determinant for successful HR. This framework can be used for targeting other harmful diseases in animal populations.

The Indian cobra reference genome and transcriptome enables comprehensive identification of venom toxins

Snakebite envenoming is a serious and neglected tropical disease that kills ~100,000 people annually. High-quality, genome-enabled comprehensive characterization of toxin genes will facilitate development of effective humanized recombinant antivenom. We report a de novo near-chromosomal genome assembly of Naja naja, the Indian cobra, a highly venomous, medically important snake. Our assembly has a scaffold N50 of 223.35 Mb, with 19 scaffolds containing 95% of the genome. Of the 23,248 predicted protein-coding genes, 12,346 venom-gland-expressed genes constitute the ‘venom-ome’ and this included 139 genes from 33 toxin families. Among the 139 toxin genes were 19 ‘venom-ome-specific toxins’ (VSTs) that showed venom-gland-specific expression, and these probably encode the minimal core venom effector proteins. Synthetic venom reconstituted through recombinant VST expression will aid in the rapid development of safe and effective synthetic antivenom. Additionally, our genome could serve as a reference for snake genomes, support evolutionary studies and enable venom-driven drug discovery.

Analysis of a near-chromosomal genome assembly and transcriptome profiling of the Indian cobra identifies genes expressed in the venom glands. These data should help develop a new antivenom.

Ten years of the horse reference genome: insights into equine biology, domestication and population dynamics in the post-genome era

The horse reference genome from the Thoroughbred mare Twilight has been available for a decade and, together with advances in genomics technologies, has led to unparalleled developments in equine genomics. At the core of this progress is the continuing improvement of the quality, contiguity and completeness of the reference genome, and its functional annotation. Recent achievements include the release of the next version of the reference genome (EquCab3.0) and generation of a reference sequence for the Y chromosome. Horse satellite‐free centromeres provide unique models for mammalian centromere research. Despite extremely low genetic diversity of the Y chromosome, it has been possible to trace patrilines of breeds and pedigrees and show that Y variation was lost in the past approximately 2300 years owing to selective breeding. The high‐quality reference genome has led to the development of three different SNP arrays and WGSs of almost 2000 modern individual horses. The collection of WGS of hundreds of ancient horses is unique and not available for any other domestic species. These tools and resources have led to global population studies dissecting the natural history of the species and genetic makeup and ancestry of modern breeds. Most importantly, the available tools and resources, together with the discovery of functional elements, are dissecting molecular causes of a growing number of Mendelian and complex traits. The improved understanding of molecular underpinnings of various traits continues to benefit the health and performance of the horse whereas also serving as a model for complex disease across species.

Chromosome-Level Alpaca Reference Genome VicPac3.1 Improves Genomic Insight Into the Biology of New World Camelids

The development of high-quality chromosomally assigned reference genomes constitutes a key feature for understanding genome architecture of a species and is critical for the discovery of the genetic blueprints of traits of biological significance. South American camelids serve people in extreme environments and are important fiber and companion animals worldwide. Despite this, the alpaca reference genome lags far behind those available for other domestic species. Here we produced a chromosome-level improved reference assembly for the alpaca genome using the DNA of the same female Huacaya alpaca as in previous assemblies. We generated 190X Illumina short-read, 8X Pacific Biosciences long-read and 60X Dovetail Chicago^® chromatin interaction scaffolding data for the assembly, used testis and skin RNAseq data for annotation, and cytogenetic map data for chromosomal assignments. The new assembly VicPac3.1 contains 90% of the alpaca genome in just 103 scaffolds and 76% of all scaffolds are mapped to the 36 pairs of the alpaca autosomes and the X chromosome. Preliminary annotation of the assembly predicted 22,462 coding genes and 29,337 isoforms. Comparative analysis of selected regions of the alpaca genome, such as the major histocompatibility complex (MHC), the region involved in the Minute Chromosome Syndrome (MCS) and candidate genes for high-altitude adaptations, reveal unique features of the alpaca genome. The alpaca reference genome VicPac3.1 presents a significant improvement in completeness, contiguity and accuracy over VicPac2 and is an important tool for the advancement of genomics research in all New World camelids.

Chromosomal Localization of Candidate Genes for Fiber Growth and Color in Alpaca (Vicugna pacos)

The alpaca (Vicugna pacos) is an economically important and cultural signature species in Peru. Thus, molecular genomic information about the genes underlying the traits of interest, such as fiber properties and color, is critical for improved breeding and management schemes. Current knowledge about the alpaca genome, particularly the chromosomal location of such genes of interest is limited and lags far behind other livestock species. The main objective of this work was to localize alpaca candidate genes for fiber growth and color using fluorescence in situ hybridization (FISH). We report the mapping of candidate genes for fiber growth COL1A1, CTNNB1, DAB2IP, KRT15, KRTAP13-1, and TNFSF12 to chromosomes 16, 17, 4, 16, 1, and 16, respectively. Likewise, we report the mapping of candidate genes for fiber color ALX3, NCOA6, SOX9, ZIC1, and ZIC5 to chromosomes 9, 19, 16, 1, and 14, respectively. In addition, since KRT15 clusters with five other keratin genes (KRT31, KRT13, KRT9, KRT14, and KRT16) in scaffold 450 (Vic.Pac 2.0.2), the entire gene cluster was assigned to chromosome 16. Similarly, mapping NCOA6 to chromosome 19, anchored scaffold 34 with 8 genes, viz., RALY, EIF2S2, XPOTP1, ASIP, AHCY, ITCH, PIGU, and GGT7 to chromosome 19. These results are concordant with known conserved synteny blocks between camelids and humans, cattle and pigs.

Comparative FISH-Mapping of MC1R, ASIP, and TYRP1 in New and Old World Camelids and Association Analysis With Coat Color Phenotypes in the Dromedary (Camelus dromedarius)

Melanocortin 1 receptor (MC1R), the agouti signaling protein (ASIP), and tyrosinase related protein 1 (TYRP1) are among the major regulators of pigmentation in mammals. Recently, MC1R and ASIP sequence variants were associated with white and black/dark brown coat colors, respectively, in the dromedary. Here we confirmed this association by independent sequencing and mutation discovery of MC1R and ASIP coding regions and by TaqMan genotyping in 188 dromedaries from Saudi Arabia and United States, including 38 black, 53 white, and 97 beige/brown/red animals. We showed that heterozygosity for a missense mutation c.901C > T in MC1R is sufficient for the white coat color suggesting a possible dominant negative effect. Likewise, we confirmed that the majority of black dromedaries were homozygous for a frameshift mutation in ASIP exon 2, except for 4 animals, which were heterozygous. In search for additional mutations underlying the black color, we identified another frameshift mutation in ASIP exon 4 and 6 new variants in MC1R including a significantly associated SNP in 3'UTR. In pursuit of sequence variants that may modify dromedary wild-type color from dark-reddish brown to light beige, we identified 4 SNPs and one insertion in TYRP1 non-coding regions. However, none of these were associated with variations in wild-type colors. Finally, the three genes were cytogenetically mapped in New World (alpaca) and Old World (dromedary and Bactrian camel) camelids. The MC1R was assigned to chr21, ASIP to chr19 and TYRP1 to chr4 in all 3 species confirming extensive conservation of camelid karyotypes. Notably, while the locations of ASIP and TYRP1 were in agreement with human-camelid comparative map, mapping MC1R identified a new evolutionary conserved synteny segment between camelid chromosome 21 and HSA16. The findings contribute to coat color genomics and the development of molecular tests in camelids and toward the chromosome level reference assemblies of camelid genomes.

An Autosomal Translocation 73,XY,t(12;20)(q11;q11) in an Infertile Male Llama (Lama glama) With Teratozoospermia

Structural chromosome abnormalities, such as translocations and inversions occasionally occur in all livestock species and are typically associated with reproductive and developmental disorders. Curiously, only a few structural chromosome aberrations have been reported in camelids, and most involved sex chromosomes. This can be attributed to a high diploid number (2n = 74) and complex chromosome morphology, which makes unambiguous identification of camelid chromosomes difficult. Additionally, molecular tools for camelid cytogenetics are sparse and have become available only recently. Here we present a case report about an infertile male llama with teratozoospermia and abnormal chromosome number 2n = 73,XY. This llama carries an autosomal translocation of chromosomes 12 and 20, which is the likely cause of defective spermatogenesis and infertility in this individual. Our analysis underlines the power of molecular cytogenetics methods over conventional banding-based chromosome analysis for explicit identification of normal and aberrant chromosomes in camelid karyotypes. This is the first case of a translocation and the first autosomal aberration reported in any camelid species. It is proof of principle that, like in other mammalian species, structural chromosome abnormalities contribute to reproductive disorders in camelids.

The horse Y chromosome as an informative marker for tracing sire lines

Analysis of the Y chromosome is the best-established way to reconstruct paternal family history in humans. Here, we applied fine-scaled Y-chromosomal haplotyping in horses with biallelic markers and demonstrate the potential of our approach to address the ancestry of sire lines. We de novo assembled a draft reference of the male-specific region of the Y chromosome from Illumina short reads and then screened 5.8 million basepairs for variants in 130 specimens from intensively selected and rural breeds and nine Przewalski's horses. Among domestic horses we confirmed the predominance of a young'crown haplogroup' in Central European and North American breeds. Within the crown, we distinguished 58 haplotypes based on 211 variants, forming three major haplogroups. In addition to two previously characterised haplogroups, one observed in Arabian/Coldblooded and the other in Turkoman/Thoroughbred horses, we uncovered a third haplogroup containing Iberian lines and a North African Barb Horse. In a genealogical showcase, we distinguished the patrilines of the three English Thoroughbred founder stallions and resolved a historic controversy over the parentage of the horse 'Galopin', born in 1872. We observed two nearly instantaneous radiations in the history of Central and Northern European Y-chromosomal lineages that both occurred after domestication 5,500 years ago.

Evolutionary conservation of Y Chromosome ampliconic gene families despite extensive structural variation

Despite claims that the mammalian Y Chromosome is on a path to extinction, comparative sequence analysis of primate Y Chromosomes has shown the decay of the ancestral single-copy genes has all but ceased in this eutherian lineage. The suite of single-copy Y-linked genes is highly conserved among the majority of eutherian Y Chromosomes due to strong purifying selection to retain dosage-sensitive genes. In contrast, the ampliconic regions of the Y Chromosome, which contain testis-specific genes that encode the majority of the transcripts on eutherian Y Chromosomes, are rapidly evolving and are thought to undergo species-specific turnover. However, ampliconic genes are known from only a handful of species, limiting insights into their long-term evolutionary dynamics. We used a clone-based sequencing approach employing both long- and short-read sequencing technologies to assemble ∼2.4 Mb of representative ampliconic sequence dispersed across the domestic cat Y Chromosome, and identified the major ampliconic gene families and repeat units. We analyzed fluorescence in situ hybridization, qPCR, and whole-genome sequence data from 20 cat species and revealed that ampliconic gene families are conserved across the cat family Felidae but show high transcript diversity, copy number variation, and structural rearrangement. Our analysis of ampliconic gene evolution unveils a complex pattern of long-term gene content stability despite extensive structural variation on a nonrecombining background.

Generation of an equine biobank to be used for Functional Annotation of Animal Genomes project

The Functional Annotation of Animal Genomes (FAANG) project aims to identify genomic regulatory elements in both sexes across multiple stages of development in domesticated animals. This study represents the first stage of the FAANG project for the horse, Equus caballus. A biobank of 80 tissue samples, two cell lines and six body fluids was created from two adult Thoroughbred mares. Ante-mortem assessments included full physical examinations, lameness, ophthalmologic and neurologic evaluations. Complete blood counts and serum biochemistries were also performed. At necropsy, in addition to tissue samples, aliquots of serum, ethylenediaminetetraacetic acid (EDTA) plasma, heparinized plasma, cerebrospinal fluid, synovial fluid, urine and microbiome samples from all regions of the gastrointestinal and urogenital tracts were collected. Epidermal keratinocytes and dermal fibroblasts were cultured from skin samples. All tissues were grossly and histologically evaluated by a board-certified veterinary pathologist. The results of the clinical and pathological evaluations identified subclinical eosinophilic and lymphocytic infiltration throughout the length of the gastrointestinal tract as well as a mild clinical lameness in both animals. Each sample was cryo-preserved in multiple ways, and nuclei were extracted from selected tissues. These samples represent the first published systemically healthy equine-specific biobank with extensive clinical phenotyping ante- and post-mortem. The tissues in the biobank are intended for community-wide use in the functional annotation of the equine genome. The use of the biobank will improve the quality of the reference annotation and allow all equine researchers to elucidate unknown genomic and epigenomic causes of disease.

Comparative Chromosome Painting in Genets (Carnivora, Viverridae, Genetta), the Only Known Feliforms with a Highly Rearranged Karyotype

Mammalian carnivores have been extensively studied by cross-species chromosome painting, which indicated a high degree of karyotypic conservatism in the cat-like suborder Feliformia relative to the ancestral carnivore karyotype (ACK). The first exception to this high degree of karyotypic conservation in feliforms was recently confirmed in genets, mesocarnivores belonging to the basal family Viverridae. Here, we present a comparative analysis of the chromosome rearrangements among 2 subspecies of the small-spotted genet Genetta genetta (the Iberian nominate and the Arabian grantii) and the panther genet G. maculata, the 2 most common and widespread genets, using whole-chromosome paints from the domestic cat (Felis catus). The chromosome homology maps and the presence of numerous interstitial telomeric sites in both genet species strengthen the hypothesis that a highly rearranged karyotype compared to the ACK may occur throughout Genetta. The karyotype of G. maculata appears to have undergone more rearrangements than that of G. genetta, which is an older lineage. Notably, we identified a tandem fusion distinguishing G. g. genetta and G. g.grantii. As G. g. grantii is morphologically and genetically distinctive, and tandem fusions have been associated with substantial postzygotic isolation in mammals, this cytogenetic finding flags the subspecies for future taxonomic investigations.

Horse Y chromosome assembly displays unique evolutionary features and putative stallion fertility genes

Dynamic evolutionary processes and complex structure make the Y chromosome among the most diverse and least understood regions in mammalian genomes. Here, we present an annotated assembly of the male specific region of the horse Y chromosome (eMSY), representing the first comprehensive Y assembly in odd-toed ungulates. The eMSY comprises single-copy, equine specific multi-copy, PAR transposed, and novel ampliconic sequence classes. The eMSY gene density approaches that of autosomes with the highest number of retained X-Y gametologs recorded in eutherians, in addition to novel Y-born and transposed genes. Horse, donkey and mule testis RNAseq reveals several candidate genes for stallion fertility. A novel testis-expressed XY ampliconic sequence class, ETSTY7, is shared with the parasite Parascaris genome, providing evidence for eukaryotic horizontal transfer and inter-chromosomal mobility. Our study highlights the dynamic nature of the Y and provides a reference sequence for improved understanding of equine male development and fertility.

Long-term expansion of primary equine keratinocytes that maintain the ability to differentiate into stratified epidermis

BACKGROUND

Skin injuries in horses frequently lead to chronic wounds that lack a keratinocyte cover essential for healing. The limited proliferation of equine keratinocytes using current protocols has limited their use for regenerative medicine. Previously, equine induced pluripotent stem cells (eiPSCs) have been produced, and eiPSCs could be differentiated into equine keratinocytes suitable for stem cell-based skin constructs. However, the procedure is technically challenging and time-consuming. The present study was designed to evaluate whether conditional reprogramming (CR) could expand primary equine keratinocytes rapidly in an undifferentiated state but retain their ability to differentiate normally and form stratified epithelium.

METHODS

Conditional reprogramming was used to isolate and propagate two equine keratinocyte cultures. PCR and FISH were employed to evaluate the equine origin of the cells and karyotyping to perform a chromosomal count. FACS analysis and immunofluorescence were used to determine the purity of equine keratinocytes and their proliferative state. Three-dimensional air-liquid interphase method was used to test the ability of cells to differentiate and form stratified squamous epithelium.

RESULTS

Conditional reprogramming was an efficient method to isolate and propagate two equine keratinocyte cultures. Cells were propagated at the rate of 2.39 days/doubling for more than 40 population doublings. A feeder-free culture method was also developed for long-term expansion. Rock-inhibitor is critical for both feeder and feeder-free conditions and for maintaining the proliferating cells in a stem-like state. PCR and FISH validated equine-specific markers in the cultures. Karyotyping showed normal equine 64, XY chromosomes. FACS using pan-cytokeratin antibodies showed a pure population of keratinocytes. When ROCK inhibitor was withdrawn and the cells were transferred to a three-dimensional air-liquid culture, they formed a well-differentiated stratified squamous epithelium, which was positive for terminal differentiation markers.

CONCLUSIONS

Our results prove that conditional reprogramming is the first method that allows for the rapid and continued in vitro propagation of primary equine keratinocytes. These unlimited supplies of autologous cells could be used to generate transplants without the risk of immune rejection. This offers the opportunity for treating recalcitrant horse wounds using autologous transplantation.

Construction of two whole genome radiation hybrid panels for dromedary (Camelus dromedarius): 5000RAD and 15000RAD

The availability of genomic resources including linkage information for camelids has been very limited. Here, we describe the construction of a set of two radiation hybrid (RH) panels (5000RAD and 15000RAD) for the dromedary (Camelus dromedarius) as a permanent genetic resource for camel genome researchers worldwide. For the 5000RAD panel, a total of 245 female camel-hamster radiation hybrid clones were collected, of which 186 were screened with 44 custom designed marker loci distributed throughout camel genome. The overall mean retention frequency (RF) of the final set of 93 hybrids was 47.7%. For the 15000RAD panel, 238 male dromedary-hamster radiation hybrid clones were collected, of which 93 were tested using 44 PCR markers. The final set of 90 clones had a mean RF of 39.9%. This 15000RAD panel is an important high-resolution complement to the main 5000RAD panel and an indispensable tool for resolving complex genomic regions. This valuable genetic resource of dromedary RH panels is expected to be instrumental for constructing a high resolution camel genome map. Construction of the set of RH panels is essential step toward chromosome level reference quality genome assembly that is critical for advancing camelid genomics and the development of custom genomic tools.

Chromosome Aberrations and Fertility Disorders in Domestic Animals

The association between chromosomal abnormalities and reduced fertility in domestic animals is well recorded and has been studied for decades. Chromosome aberrations directly affect meiosis, gametogenesis, and the viability of zygotes and embryos. In some instances, balanced structural rearrangements can be transmitted, causing fertility problems in subsequent generations. Here, we aim to give a comprehensive overview of the current status and future prospects of clinical cytogenetics of animal reproduction by focusing on the advances in molecular cytogenetics during the genomics era. We describe how advancing knowledge about animal genomes has improved our understanding of connections between gross structural or molecular chromosome variations and reproductive disorders. Further, we expand on a key area of reproduction genetics: cytogenetics of animal gametes and embryos. Finally, we describe how traditional cytogenetics is interfacing with advanced genomics approaches, such as array technologies and next-generation sequencing, and speculate about the future prospects.

Fertility Assessment in Sorraia Stallions by Sperm-Fish and Fkbp6 Genotyping

The Sorraia, a critically endangered indigenous Iberian horse breed, is characterized by low genetic variability, high rate of inbreeding, bad sperm quality and subfertility. Here, we studied 11 phenotypically normal but subfertile Sorraia stallions by karyotyping, sex chromosome sperm-FISH and molecular analysis of FKBP6 - a susceptibility locus for impaired acrosome reaction (IAR). The stallions had normal sperm concentration (>300 million cells/ml), but the numbers of progressively motile sperm (21%) and morphologically normal sperm (28%) were invariably low. All stallions had a normal 64,XY karyotype. The majority of sperm (89%) had normal haploid sex chromosome content, although 11% of sperm carried various sex chromosome aneuploidies. No correlation was found between the percentage of sperm sex chromosome abnormalities and inbreeding, sperm morphology or stallion age. Direct sequencing of FKBP6 exon 4 for SNPs g.11040315G>A and g.11040379C>A revealed that none of the stallions had the susceptibility genotype (A/A-A/A) for IAR. Instead, all animals had a G/G-A/A genotype - a testimony of low genetic variability. The findings ruled out chromosomal abnormalities and genetic predisposition for IAR as contributing factors for subfertility. However, low fertility of the Sorraia stallions could be partly attributed to relatively higher rate of sex chromosome aneuploidies in the sperm.

Analysis of genomic copy number variation in equine recurrent airway obstruction (heaves)

We explored the involvement of genomic copy number variants (CNVs) in susceptibility to recurrent airway obstruction (RAO), or heaves-an asthmalike inflammatory disease in horses. Analysis of 16 RAO-susceptible (cases) and six RAO-resistant (control) horses on a custom-made whole-genome 400K equine tiling array identified 245 CNV regions (CNVRs), 197 previously known and 48 new, distributed on all horse autosomes and the X chromosome. Among the new CNVRs, 30 were exclusively found in RAO cases and were further analyzed by quantitative PCR, including additional cases and controls. Suggestive association (P = 0.03; corrected P = 0.06) was found between RAO and a loss on chromosome 5 involving NME7, a gene necessary for ciliary functions in lungs and involved in primary ciliary dyskinesia in humans. The CNVR could be a potential marker for RAO susceptibility but needs further study in additional RAO cohorts. Other CNVRs were not associated with RAO, although several involved genes of interest, such as SPI2/SERPINA1 from the serpin gene family, which are associated with chronic obstructive pulmonary disease and asthma in humans. The SPI2/SERPINA1 CNVR showed striking variation among horses, but it was not significantly different between RAO cases and controls. The findings provide baseline information on the relationship between CNVs and RAO susceptibility. Discovery of new CNVs and the use of a larger population of RAO-affected and control horses are needed to shed more light on their significance in modulating this complex and heterogeneous disease.