Development and application of camelid molecular cytogenetic tools

Ovarian Dynamics and Pathological Conditions in Camelids

Domestic camelids comprise Dromedary, Bactrian (old-world camelids, OWC), alpaca, and Lama (new-world camelids, NWC). They possess unique reproductive physiology distinct from other livestock. A key feature is being induced ovulators, triggered by beta nerve growth factor (βNGF) present in the seminal plasma. Proper ovarian function is crucial for reproductive management and biotechnologies. In the absence of ovulation, camelids display overlapping follicular waves, with puberty and sexual receptivity influenced by factors like weight and age. Follicular growth patterns vary among these species and are better monitored via transrectal ultrasonography. Physiological variations in ovarian dynamics and hormonal activity differ among members of the Camelidae. Understanding these differences is essential for production and conservation efforts, particularly in wild species. Ovarian disorders, including failure to ovulate or respond to therapy, frequently cause subfertility or infertility in affected females. Superstimulatory and synchronisation protocols that rely on exogenous progestin, eCG, and FSH supplementation have been developed to address some causes of irregular ovarian dynamics. This paper reviews the current knowledge on follicular wave dynamics, their manipulation, and the common ovarian disorders affecting reproduction in camelids.

Observations on chromosome-specific sequencing for the construction of cross-species chromosome homology maps and its resolution of human:alpaca homology

Background

The history of comparative chromosome mapping is briefly reviewed, with discussion about the problem that occurs in chromosome painting when size heteromorphisms between homologues cause contamination in chromosomes sorted by flow cytometry that are used in the preparation of chromosome-specific DNA probes.

Main body

As an example, we show in the alpaca (Vicagna pacos) that sequencing of contaminated chromosome sorts can reveal chromosome homologies from alignment with human and mouse genome reference sequences. The procedure identifies syntenic blocks of DNA separated in the human karyotype that are associated in the closely related alpaca and dromedary (Camelus dromedarius) karyotypes. This example provides proof of principal for the validity of the method for comparative chromosome mapping.

Conclusion

It is suggested that the approach presented here may have application in the construction of comparative chromosome maps between distantly related taxa, such as monotremes and mammals.

Clinical and Histological Features of Ovarian Hypoplasia/Dysgenesis in Alpacas

Alpacas have a high incidence of congenital reproductive tract abnormalities, including ovarian hypoplasia/dysgenesis. Diagnosis of this condition is often challenging. The present study describes the clinical, ultrasonographic, and histologic features of ovarian hypoplasia/dysgenesis syndrome in 5 female alpacas. Additionally, serum AMH levels were compared between female alpacas diagnosed with ovarian hypoplasia/dysgenesis and a group of reproductively sound females (n = 11). The syndrome was suspected based on the presence of an infantile uterus and lack of ovaries by ultrasonography and laparoscopy. All females had normal female karyotype (n = 74 XX), but one presented a minute chromosome. The ovaries from these cases showed 3 main histological classifications: hypoplasia (n = 2), dysgenesis (n = 2), and dysplasia (n = 1). Serum AMH levels in affected females were significantly lower (P < 0.05) than those of reproductively sound control females. In conclusion, Serum AMH level may be helpful in the rapid diagnosis of ovarian hypoplasia/dysgenesis syndrome in alpacas. Furthermore, this syndrome in alpacas presents a variety of histological features. Different mechanisms may be involved in the derangement of ovarian differentiation. Further studies are needed to elucidate the causes of the syndrome.

State of the art on the physical mapping of the Y-chromosome in the <i>Bovidae</i> and comparison with other species

The next generation sequencing has significantly contributed to clarify the genome structure of many species of zootechnical interest. However, to date, some portions of the genome, especially those linked to a heterogametic nature such as the Y chromosome, are difficult to assemble and many gaps are still present. It is well known that the fluorescence in situ hybridization (FISH) is an excellent tool for identifying genes unequivocably mapped on chromosomes. Therefore, FISH can contribute to the localization of unplaced genome sequences, as well as to correct assembly errors generated by comparative bioinformatics. To this end, it is necessary to have starting points; therefore, in this study, we reviewed the physically mapped genes on the Y chromosome of cattle, buffalo, sheep, goats, pigs, horses and alpacas. A total of 208 loci were currently mapped by FISH. 89 were located in the male-specific region of the Y chromosome (MSY) and 119 were identified in the pseudoautosomal region (PAR). The loci reported in MSY and PAR were respectively: 18 and 25 in Bos taurus, 5 and 7 in Bubalus bubalis, 5 and 24 in Ovis aries, 5 and 19 in Capra hircus, 10 and 16 in Sus scrofa, 46 and 18 in Equus caballus. While in Vicugna pacos only 10 loci are reported in the PAR region. The correct knowledge and assembly of all genome sequences, including those of genes mapped on the Y chromosome, will help to elucidate their biological processes, as well as to discover and exploit potentially epistasis effects useful for selection breeding programs.

Investigation of male and female infertility in llamas and alpacas

Llamas and alpacas are important production animals in South America, with increasing interest in other parts of the world. Poor reproductive efficiency combined with several unique anatomical and physiological reproductive features offer challenges in the diagnosis and treatment of infertility in camelids. This review presents an approach to the clinical investigation and common causes of infertility and subfertility in the male and female. The selection of males for breeding should be made based on complete evaluation to eliminate congenital and possibly hereditary disorders. Common disorders of the male reproductive system include testicular hypoplasia, testicular and epididymal cysts and testicular degeneration. Semen evaluation presents some challenges owing to the viscous nature of the ejaculate in these species. Females should be screened for congenital genital defects before breeding. Causes of subfertility in the female are dominated by ovarian and uterine disorders. A systematic clinical approach and the use of endometrial biopsy and advanced techniques, such as laparoscopy, allow early identification of these disorders. Further research is needed for continued understanding of the reproductive pathological processes in these species.

Cytogenetic Mapping of 35 New Markers in the Alpaca (Vicugna pacos)

Alpaca is a camelid species of broad economic, biological and biomedical interest, and an essential part of the cultural and historical heritage of Peru. Recently, efforts have been made to improve knowledge of the alpaca genome, and its genetics and cytogenetics, to develop molecular tools for selection and breeding. Here, we report cytogenetic mapping of 35 new markers to 19 alpaca autosomes and the X chromosome. Twenty-eight markers represent alpaca SNPs, of which 17 are located inside or near protein-coding genes, two are in ncRNA genes and nine are intergenic. The remaining seven markers correspond to candidate genes for fiber characteristics (BMP4, COL1A2, GLI1, SFRP4), coat color (TYR) and development (CHD7, PAX7). The results take the tally of cytogenetically mapped markers in alpaca to 281, covering all 36 autosomes and the sex chromosomes. The new map assignments overall agree with human–camelid conserved synteny data, except for mapping BMP4 to VPA3, suggesting a hitherto unknown homology with HSA14. The findings validate, refine and correct the current alpaca assembly VicPac3.1 by anchoring unassigned sequence scaffolds, and ordering and orienting assigned scaffolds. The study contributes to the improvement in the alpaca reference genome and advances camelid molecular cytogenetics.

Casein Gene Cluster in Camelids: Comparative Genome Analysis and New Findings on Haplotype Variability and Physical Mapping

The structure of casein genes has been fully understood in llamas, whereas in other camelids, this information is still incomplete. In fact, structure and polymorphisms have been identified in three (CSN1S1, αs1-CN; CSN2, β-CN; CSN3, κ-CN) out of four casein genes, whereas controversial information is available for the CSN1S2 (αs2-CN) in terms of structure and genetic diversity. Data from the genome analysis, whose assembly is available for feral camel, Bactrian, dromedary, and alpaca, can contribute to a better knowledge. However, a majority of the scaffolds available in GenBank are still unplaced, and the comparative annotation is often inaccurate or lacking.Therefore, the aims of this study are 1) to perform a comparative genome analysis and synthesize the literature data on camelids casein cluster; 2) to analyze the casein variability in two dromedary populations (Sudanese and Nigerian) using polymorphisms at CSN1S1 (c.150G > T), CSN2 (g.2126A > G), and CSN3 (g.1029T > C); and 3) to physically map the casein cluster in alpaca. Exon structures, gene and intergenic distances, large insertion/deletion events, SNPs, and microsatellites were annotated. In all camelids, the CSN1S2 consists of 17 exons, confirming the structure of llama CSN1S2 gene. The comparative analysis of the complete casein cluster (∼190kb) shows 12,818 polymorphisms. The most polymorphic gene is the CSN1S1 (99 SNPs in Bactrian vs. 248 in dromedary vs. 626 in alpaca). The less polymorphic is the CSN3 in the Bactrian (22 SNPs) and alpaca (301 SNPs), whereas it is the CSN1S2 in dromedary (79 SNPs). In the two investigated dromedary populations, the allele frequencies for the three markers are slightly different: the allele C at CSN1S1 is very rare in Nigerian (0.054) and Sudanese dromedaries (0.094), whereas the frequency of the allele G at CSN2 is almost inverted. Haplotype analysis evidenced GAC as the most frequent (0.288) and TGC as the rarest (0.005). The analysis of R-banding metaphases hybridized with specific probes mapped the casein genes on chromosome 2q21 in alpaca. These data deepen the information on the structure of the casein cluster in camelids and add knowledge on the cytogenetic map and haplotype variability.

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.

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.

The major histocompatibility complex in Old World camelids and low polymorphism of its class II genes

BACKGROUND

The Major Histocompatibility Complex (MHC) is a genomic region containing genes with crucial roles in immune responses. MHC class I and class II genes encode antigen-presenting molecules expressed on the cell surface. To counteract the high variability of pathogens, the MHC evolved into a region of considerable heterogeneity in its organization, number and extent of polymorphism. Studies of MHCs in different model species contribute to our understanding of mechanisms of immunity, diseases and their evolution. Camels are economically important domestic animals and interesting biomodels. Three species of Old World camels have been recognized: the dromedary (Camelus dromedarius), Bactrian camel (Camelus bactrianus) and the wild camel (Camelus ferus). Despite their importance, little is known about the MHC genomic region, its organization and diversity in camels. The objectives of this study were to identify, map and characterize the MHC region of Old World camelids, with special attention to genetic variation at selected class MHC II loci.

RESULTS

Physical mapping located the MHC region to the chromosome 20 in Camelus dromedarius. Cytogenetic and comparative analyses of whole genome sequences showed that the order of the three major sub-regions is "Centromere - Class II - Class III - Class I". DRA, DRB, DQA and DQB exon 2 sequences encoding the antigen binding site of the corresponding class II antigen presenting molecules showed high degree of sequence similarity and extensive allele sharing across the three species. Unexpectedly low extent of polymorphism with low numbers of alleles and haplotypes was observed in all species, despite different geographic origins of the camels analyzed. The DRA locus was found to be polymorphic, with three alleles shared by all three species. DRA and DQA sequences retrieved from ancient DNA samples of Camelus dromedarius suggested that additional polymorphism might exist.

CONCLUSIONS

This study provided evidence that camels possess an MHC comparable to other mammalian species in terms of its genomic localization, organization and sequence similarity. We described ancient variation at the DRA locus, monomorphic in most species. The extent of molecular diversity of MHC class II genes seems to be substantially lower in Old World camels than in other mammalian species.

Skeletal variation in Tennessee Walking Horses maps to the LCORL/NCAPG gene region

Conformation has long been a driving force in horse selection and breed creation as a predictor for performance. The Tennessee Walking Horse (TWH) ranges in size from 1.5 to 1.7 m and is often used as a trail, show, and pleasure horse. To investigate the contribution of genetics to body conformation in the TWH, we collected DNA samples, body measurements, and gait/training information from 282 individuals. We analyzed the 32 body measures with a principal component analysis. Principal component (PC)1 captured 28.5% of the trait variance, while PC2 comprised just 9.5% and PC3 6.4% of trait variance. All 32 measures correlated positively with PC1, indicating that PC1 describes overall body size. We genotyped 109 horses using the EquineSNP70 bead chip and marker association assessed the data using PC1 scores as a phenotype. Mixed-model linear analysis (EMMAX) revealed a well-documented candidate locus on ECA3 (raw P = 3.86 × 10(-9)) near the LCORL gene. A custom genotyping panel enabled fine-mapping of the PC1 body-size trait to the 3'-end of the LCORL gene (P = 7.09 × 10(-10)). This position differs from other reports suggesting single nucleotide polymorphisms (SNPs) upstream of the LCORL coding sequence regulate expression of the gene and, therefore, body size in horses. Fluorescent in situ hybridization analysis defined the position of a highly homologous 5 kb retrogene copy of LCORL (assigned to unplaced contigs of the EquCab 2.0 assembly) at ECA9 q12-q13. This is the first study to identify putative causative SNPs within the LCORL transcript itself, which are associated with skeletal size variation in horses.

A comprehensive whole-genome integrated cytogenetic map for the alpaca (Lama pacos)

Genome analysis of the alpaca (Lama pacos, LPA) has progressed slowly compared to other domestic species. Here, we report the development of the first comprehensive whole-genome integrated cytogenetic map for the alpaca using fluorescence in situ hybridization (FISH) and CHORI-246 BAC library clones. The map is comprised of 230 linearly ordered markers distributed among all 36 alpaca autosomes and the sex chromosomes. For the first time, markers were assigned to LPA14, 21, 22, 28, and 36. Additionally, 86 genes from 15 alpaca chromosomes were mapped in the dromedary camel (Camelus dromedarius, CDR), demonstrating exceptional synteny and linkage conservation between the 2 camelid genomes. Cytogenetic mapping of 191 protein-coding genes improved and refined the known Zoo-FISH homologies between camelids and humans: we discovered new homologous synteny blocks (HSBs) corresponding to HSA1-LPA/CDR11, HSA4-LPA/CDR31 and HSA7-LPA/CDR36, and revised the location of breakpoints for others. Overall, gene mapping was in good agreement with the Zoo-FISH and revealed remarkable evolutionary conservation of gene order within many human-camelid HSBs. Most importantly, 91 FISH-mapped markers effectively integrated the alpaca whole-genome sequence and the radiation hybrid maps with physical chromosomes, thus facilitating the improvement of the sequence assembly and the discovery of genes of biological importance.

A cytogenetic and comparative map of camelid chromosome 36 and the minute in alpacas

Recent advances in camelid genomics have provided draft sequence assemblies and the first comparative and gene maps for the dromedary (CDR) and the alpaca (LPA). However, no map information is currently available for the smallest camelid autosome-chr36. The chromosome is also of clinical interest because of its involvement in the minute chromosome syndrome (MCS) in infertile alpacas. Here, we developed molecular markers for camelid chr36 by direct sequencing CDR36 and LPA minute and by bioinformatics analysis of alpaca unplaced sequence scaffolds. We constructed a cytogenetic map for chr36 in the alpaca, llama, and dromedary and showed its homology to human chromosome 7 (HSA7) at 49.8-55.5 Mb. The chr36 map comprised seven markers, including two genes-ZPBP and WVC2. Comparative status of HSA7 was further refined by cytogenetic mapping of 16 HSA7 orthologs in camelid chromosomes 7 and 18 and by the analysis of HSA7-conserved synteny blocks across 11 vertebrate species. Finally, mapping chr36 markers in infertile alpacas confirmed that the minute chromosome was a derivative of chr36, but the small size was not a result of a large deletion or a translocation. Instead, cytogenetic mapping of 5.8S, 18S, and 28S rRNA genes (nucleolus organizer region (NOR)) revealed that the size difference between chr36 homologs in infertile alpacas was due to a heterozygous presence of NOR, whereas chr36 in fertile alpacas had no NOR. We theorized that the heterozygous NOR might affect chr36 pairing, recombination, and segregation in meiosis and, thus fertility.

Comparative organization and gene expression profiles of the porcine pseudoautosomal region

The pseudoautosomal region (PAR) has important biological functions in spermatogenesis, male fertility and early development. Even though pig (Sus scrofa, SSC) is an agriculturally and biomedically important species, and its genome is sequenced, current knowledge about the porcine PAR is sparse. Here we defined the PAR in SSCXp/Yp by demarcating the sequence of the pseudoautosomal boundary at X:6,743,567 bp in intron 3-4 of SHROOM2 and showed that SHROOM2 is truncated in SSCY. Cytogenetic mapping of 20 BAC clones containing 15 PAR and X-specific genes revealed that the pig PAR is largely collinear with other mammalian PARs or Xp terminal regions. The results improved the current SSCX sequence assembly and facilitated distinction between the PAR and X-specific genes to study their expression in adult and embryonic tissues. A pilot analysis showed that the PAR genes are expressed at higher levels than X-specific genes during early development, whereas the expression of PAR genes was higher at day 60 compared to day 26, and higher in embryonic tissues compared to placenta. The findings advance the knowledge about the comparative organization of the PAR in mammals and suggest that the region might have important functions in early development in pigs.