Cross-species chromosome painting among camel, cattle, pig and human: further insights into the putative Cetartiodactyla ancestral karyotype

Genetic determinants of micronucleus formation in vivo

Genomic instability arising from defective responses to DNA damage1 or mitotic chromosomal imbalances2 can lead to the sequestration of DNA in aberrant extranuclear structures called micronuclei (MN). Although MN are a hallmark of ageing and diseases associated with genomic instability, the catalogue of genetic players that regulate the generation of MN remains to be determined. Here we analyse 997 mouse mutant lines, revealing 145 genes whose loss significantly increases (n = 71) or decreases (n = 74) MN formation, including many genes whose orthologues are linked to human disease. We found that mice null for Dscc1, which showed the most significant increase in MN, also displayed a range of phenotypes characteristic of patients with cohesinopathy disorders. After validating the DSCC1-associated MN instability phenotype in human cells, we used genome-wide CRISPR-Cas9 screening to define synthetic lethal and synthetic rescue interactors. We found that the loss of SIRT1 can rescue phenotypes associated with DSCC1 loss in a manner paralleling restoration of protein acetylation of SMC3. Our study reveals factors involved in maintaining genomic stability and shows how this information can be used to identify mechanisms that are relevant to human disease biology1.

Genome-Wide Association Study of Fiber Diameter in Alpacas

The aim of this study was the identification of candidate genomic regions associated with fiber diameter in alpacas. DNA samples were collected from 1011 female Huacaya alpacas from two geographical Andean regions in Peru (Pasco and Puno), and three alpaca farms within each region. The samples were genotyped using an Affymetrix Custom Alpaca genotyping array containing 76,508 SNPs. After the quality controls, 960 samples and 51,742 SNPs were retained. Three association study methodologies were performed. The GWAS based on a linear model allowed us to identify 11 and 35 SNPs (-log10(p-values) > 4) using information on all alpacas and alpacas with extreme values of fiber diameter, respectively. The haplotype and marker analysis method allowed us to identify nine haplotypes with standardized haplotype heritability higher than six standard deviations. The selection signatures based on cross-population extended haplotype homozygosity (XP-EHH) allowed us to identify 180 SNPs with XP-EHH values greater than |3|. Four candidate regions with adjacent SNPs identified via two association methods of analysis are located on VPA6, VPA9, VPA29 and one chromosomally unassigned scaffold. This study represents the first analysis of alpaca whole genome association with fiber diameter, using a recently assembled alpaca SNP microarray.

Construction and Characterization of Immortalized Fibroblast Cell Line from Bactrian Camel

Immortalized cell lines with many advantages are widely used in various experimental contexts by many different labs. However, the absence of available cell lines poses difficulties for research in some species, such as camels. To establish an immortalized Bactrian camel fibroblast (iBCF) cell line and understand its biological characteristics, primary fibroblast cells from Bactrian camels were isolated and purified using enzymatic digestion in this study, and telomerase reverse transcriptase (hTERT) vectors were introduced into primary BCF (pBCF) for continuous passage to 80 generations after screening with G418. The cell morphology of different generations was examined under a microscope. Cell cycle and viability were evaluated by flow cytometry and CCK-8 assay, respectively. Cellular genes expression was monitored by qPCR, immunofluorescence, and Western blot, respectively. Chromosomes were determined by karyotyping. The results showed that like most other cells, both pBCF and iBCF were sensitive to nutrient concentrations and adapted to culture in the medium with 4.5 g/L glucose and 10% fetal bovine serum (FBS) concentration. hTERT gene was introduced and stably expressed in iBCF cells, which promoted BCF cell immortalization. The fibroblast specific marker vimentin (VIM) is expressed in both pBCF and iBCF, but epithelial marker cytokeratin18 (CK18) expression is weak in BCF cells. Proliferation and viability detection showed that hTERT-induced iBCF exhibits faster growth rates and higher viability than pBCF. Karyotyping showed that iBCF maintained the same number and morphology of chromosomes as the pBCF. This study demonstrated that we have successfully constructed an immortalized Bactrian camel fibroblast cell line, which was named BCF23. The establishment of the BCF23 cell line provides a foundation for expanding camel-related research.

Chromosome-level assembly of the Rangifer tarandus genome and validation of cervid and bovid evolution insights

BACKGROUND

Genome assembly into chromosomes facilitates several analyses including cytogenetics, genomics and phylogenetics. Despite rapid development in bioinformatics, however, assembly beyond scaffolds remains challenging, especially in species without closely related well-assembled and available reference genomes. So far, four draft genomes of Rangifer tarandus (caribou or reindeer, a circumpolar distributed cervid species) have been published, but none with chromosome-level assembly. This emblematic northern species is of high interest in ecological studies and conservation since most populations are declining.

RESULTS

We have designed specific probes based on Oligopaint FISH technology to upgrade the latest published reindeer and caribou chromosome-level genomes. Using this oligonucleotide-based method, we found six mis-assembled scaffolds and physically mapped 68 of the largest scaffolds representing 78% of the most recent R. tarandus genome assembly. Combining physical mapping and comparative genomics, it was possible to document chromosomal evolution among Cervidae and closely related bovids.

CONCLUSIONS

Our results provide validation for the current chromosome-level genome assembly as well as resources to use chromosome banding in studies of Rangifer tarandus.

Establishment of Bactrian Camel Induced Pluripotent Stem Cells and Prediction of Their Unique Pluripotency Genes

Induced pluripotent stem cells (iPSCs) can differentiate into all types of cells and can be used in livestock for research on biological development, genetic breeding, and in vitro genetic resource conservation. The Bactrian camel is a large domestic animal that inhabits extreme environments and holds value in the treatment of various diseases and the development of the local economy. Therefore, we transferred four mouse genes (Oct4, Sox2, Klf4, and c-Myc) into Bactrian camel fetal fibroblasts (BCFFs) using retroviruses with a large host range to obtain Bactrian camel induced pluripotent stem cells (bciPSCs). They were comprehensively identified based on cell morphology, pluripotency gene and marker expression, chromosome number, transcriptome sequencing, and differentiation potential. The results showed the pluripotency of bciPSCs. However, unlike stem cells of other species, late formation of stem cell clones was observed; moreover, the immunofluorescence of SSEA1, SSEA3, and SSEA4 were positive, and teratoma formation took four months. These findings may be related to the extremely long gestation period and species specificity of Bactrian camels. By mining RNA sequence data, 85 potential unique pluripotent genes of Bactrian camels were predicted, which could be used as candidate genes for the production of bciPSC in the future. Among them, ASF1B, DTL, CDCA5, PROM1, CYTL1, NUP210, Epha3, and SYT13 are more attractive. In conclusion, we generated bciPSCs for the first time and obtained their transcriptome information, expanding the iPSC genetic information database and exploring the applicability of iPSCs in livestock. Our results can provide an experimental basis for Bactrian camel ESC establishment, developmental research, and genetic resource conservation.

Comparative Fluorescence In Situ Hybridization (FISH) Mapping of Twenty-Three Endogenous Jaagsiekte Sheep Retrovirus (enJSRVs) in Sheep (Ovis aries) and River Buffalo (Bubalus bubalis) Chromosomes

Endogenous retroviruses (ERVs) are the remnants of ancient infections of host germline cells, thus representing key tools to study host and viral evolution. Homologous ERV sequences often map at the same genomic locus of different species, indicating that retroviral integration occurred in the genomes of the common ancestors of those species. The genome of domestic sheep (Ovis aries) harbors at least twenty-seven copies of ERVs related to the exogenous and pathogenic Jaagsiekte sheep retrovirus (JSRVs), thus referred to as enJSRVs. Some of these loci are unequally distributed between breeds and individuals of the host species due to polymorphic insertions, thereby representing invaluable tools to trace the evolutionary dynamics of virus populations within their hosts. In this study, we extend the cytogenetic physical maps of sheep and river buffalo by performing fluorescent in situ hybridization (FISH) mapping of twenty-three genetically characterized enJSRVs. Additionally, we report the first comparative FISH mapping of enJSRVs in domestic sheep (2n = 54) and river buffalo (Bubalus bubalis, 2n = 50). Finally, we demonstrate that enJSRV loci are conserved in the homologous chromosomes and chromosome bands of both species. Altogether, our results support the hypothesis that enJSRVs were present in the genomes of both species before they differentiated within the Bovidae family.

Aspects of Molecular Genetics in Dromedary Camel

Dromedary camels are unique in their morphological and physiological characteristics and are capable of providing milk and meat even under extreme environmental conditions. Like other species, the dromedary camel has also benefitted from the development of the molecular genetics to increase the knowledge about different aspect in camel genetics (genetic variation, molecular marker, parentage control, gene of interest, whole genome, dating…etc.). In this paper we review the different molecular genetic technics used in this particular species and future prospects. Dromedary genetic studies started in the end of the 1980s with phenotypic evaluation and the attempts to highlight the protein and biochemical diversity. In the 2000s, with the development of molecular markers such as microsatellites, genetic diversity of different types in several countries were estimated and microsatellites were also used for parentage control. In terms of genetic characterization, microsatellites revealed a defined global structure, differentiating East African and South Arabian dromedaries from North African, North Arabian, and South Asian individuals, respectively. Also, mitochondrialDNA sequence analysis of ancient DNA proved to be crucial in resolving domestication processes in dromedaries. Ancient and modern DNA revealed dynamics of domestication and cross-continental dispersion of the dromedary. Nuclear SNPs, single nucleotide polymorphisms changes that occur approximately each 1000 bps in the mammalian genome were also applied in some studies in dromedary. These markers are a very useful alternative to microsatellites and have been employed in some studies on genetic diversity and relevant phenotypic traits in livestock. Finally, thanks to the use of Next Generation Sequencing (NGS) the whole-genome assemblies of the dromedary (Camelus dromedarius) and a work to establish the organization of the dromedary genome at chromosome level were recently published.

Development of a 76k Alpaca (Vicugna pacos) Single Nucleotide Polymorphisms (SNPs) Microarray

Small farm producers’ sustenance depends on their alpaca herds and the production of fiber. Genetic improvement of fiber characteristics would increase their economic benefits and quality of life. The incorporation of molecular marker technology could overcome current limitations for the implementation of genetic improvement programs. Hence, the aim of this project was the generation of an alpaca single nucleotide polymorphism (SNP) microarray. A sample of 150 Huacaya alpacas from four farms, two each in Puno and Cerro de Pasco were used for SNP discovery by genotyping by sequencing (GBS). Reduced representation libraries, two per animal, were produced after DNA digestion with ApeK1 and double digestion with Pst1-Msp1. Ten alpaca genomes, sequenced at depths between 12× to 30×, and the VicPac3.1 reference genome were used for read alignments. Bioinformatics analysis discovered 76,508 SNPs included in the microarray. Candidate genes SNPs (302) for fiber quality and color are also included. The microarray SNPs cover 90.5% of the genome length with a density of about 39 ± 2.51 SNPs/Mb of DNA at an average interval of 26.45 ± 18.57 kbp. The performance was evaluated by genotyping 30 family trios and comparing them to their pedigrees, as well as comparing microarray to GBS genotypes. Concordance values of 0.93 and 0.94 for ApeK1 and Pst1-Msp1 generated SNPs were observed. Similarly, 290 fiber quality and color candidate gene SNPs were validated. Availability of this microarray will facilitate genome-wide association studies, marker-assisted selection and, in time, genomic selection.

An 8.22 Mb Assembly and Annotation of the Alpaca (Vicugna pacos) Y Chromosome

The unique evolutionary dynamics and complex structure make the Y chromosome the most diverse and least understood region in the mammalian genome, despite its undisputable role in sex determination, development, and male fertility. Here we present the first contig-level annotated draft assembly for the alpaca (Vicugna pacos) Y chromosome based on hybrid assembly of short- and long-read sequence data of flow-sorted Y. The latter was also used for cDNA selection providing Y-enriched testis transcriptome for annotation. The final assembly of 8.22 Mb comprised 4.5 Mb of male specific Y (MSY) and 3.7 Mb of the pseudoautosomal region. In MSY, we annotated 15 X-degenerate genes and two novel transcripts, but no transposed sequences. Two MSY genes, HSFY and RBMY, are multicopy. The pseudoautosomal boundary is located between SHROOM2 and HSFY. Comparative analysis shows that the small and cytogenetically distinct alpaca Y shares most of MSY sequences with the larger dromedary and Bactrian camel Y chromosomes. Most of alpaca X-degenerate genes are also shared with other mammalian MSYs, though WWC3Y is Y-specific only in alpaca/camels and the horse. The partial alpaca Y assembly is a starting point for further expansion and will have applications in the study of camelid populations and male biology.

An Indo-Pacific Humpback Dolphin Genome Reveals Insights into Chromosome Evolution and the Demography of a Vulnerable Species

The Indo-Pacific humpback dolphin (Sousa chinensis) is a small inshore species of odontocete cetacean listed as Vulnerable on the IUCN Red List. Here, we report on the evolution of S. chinensis chromosomes from its cetruminant ancestor and elucidate the evolutionary history and population genetics of two neighboring S. chinensis populations. We found that breakpoints in ancestral chromosomes leading to S. chinensis could have affected the function of genes related to kidney filtration, body development, and immunity. Resequencing of individuals from two neighboring populations in the northwestern South China Sea, Leizhou Bay and Sanniang Bay, revealed genetic differentiation, low diversity, and small contemporary effective population sizes. Demographic analyses showed a marked decrease in the population size of the two investigated populations over the last ~4,000 years, possibly related to climatic oscillations. This study implies a high risk of extinction and strong conservation requirement for the Indo-Pacific humpback dolphin.

•

Deducing chromosome evolution from ancestral Cetruminantia and ancestral Odontoceti

•

Reconstructing the demographic history of Sousa chinensis

•

Implying high risk of extinction and strong conservation requirement for S. chinensis

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.

Evolution of the Human Chromosome 13 Synteny: Evolutionary Rearrangements, Plasticity, Human Disease Genes and Cancer Breakpoints

The history of each human chromosome can be studied through comparative cytogenetic approaches in mammals which permit the identification of human chromosomal homologies and rearrangements between species. Comparative banding, chromosome painting, Bacterial Artificial Chromosome (BAC) mapping and genome data permit researchers to formulate hypotheses about ancestral chromosome forms. Human chromosome 13 has been previously shown to be conserved as a single syntenic element in the Ancestral Primate Karyotype; in this context, in order to study and verify the conservation of primate chromosomes homologous to human chromosome 13, we mapped a selected set of BAC probes in three platyrrhine species, characterised by a high level of rearrangements, using fluorescence in situ hybridisation (FISH). Our mapping data on Saguinus oedipus, Callithrix argentata and Alouatta belzebul provide insight into synteny of human chromosome 13 evolution in a comparative perspective among primate species, showing rearrangements across taxa. Furthermore, in a wider perspective, we have revised previous cytogenomic literature data on chromosome 13 evolution in eutherian mammals, showing a complex origin of the eutherian mammal ancestral karyotype which has still not been completely clarified. Moreover, we analysed biomedical aspects (the OMIM and Mitelman databases) regarding human chromosome 13, showing that this autosome is characterised by a certain level of plasticity that has been implicated in many human cancers and diseases.

Chromosome-level assembly of wild Bactrian camel genome reveals organization of immune gene loci

Camelids are characterized by their unique adaptive immune system that exhibits the generation of homodimeric heavy-chain immunoglobulins, somatic hypermutation of T-cell receptors, and low genetic diversity of major histocompatibility complex (MHC) genes. However, short-read assemblies are typically highly fragmented in these gene loci owing to their repetitive and polymorphic nature. Here, we constructed a chromosome-level assembly of wild Bactrian camel genome based on high-coverage long-read sequencing and chromatin interaction mapping. The assembly with a contig N50 of 5.37 Mb and a scaffold N50 of 76.03 Mb, represents the most contiguous camelid genome to date. The genomic organization of immunoglobulin heavy-chain locus was similar between the wild Bactrian camel and alpaca, and genes encoding for conventional and heavy-chain antibodies were intermixed. The organizations of two immunoglobulin light-chain loci and four T cell receptor loci were also fully deciphered using the new assembly. Additionally, the complete classical MHC region was resolved into a single contig. The high-quality assembly presented here provides an essential reference for future investigations examining the camelid immune system.

Comparative Chromosome Mapping of Musk Ox and the X Chromosome among Some Bovidae Species

: Bovidae, the largest family in Pecora infraorder, are characterized by a striking variability in diploid number of chromosomes between species and among individuals within a species. The bovid X chromosome is also remarkably variable, with several morphological types in the family. Here we built a detailed chromosome map of musk ox (Ovibos moschatus), a relic species originating from Pleistocene megafauna, with dromedary and human probes using chromosome painting. We trace chromosomal rearrangements during Bovidae evolution by comparing species already studied by chromosome painting. The musk ox karyotype differs from the ancestral pecoran karyotype by six fusions, one fission, and three inversions. We discuss changes in pecoran ancestral karyotype in the light of new painting data. Variations in the X chromosome structure of four bovid species nilgai bull (Boselaphus tragocamelus), saola (Pseudoryx nghetinhensis), gaur (Bos gaurus), and Kirk's Dikdik (Madoqua kirkii) were further analyzed using 26 cattle BAC-clones. We found the duplication on the X in saola. We show main rearrangements leading to the formation of four types of bovid X: Bovinae type with derived cattle subtype formed by centromere reposition and Antilopinae type with Caprini subtype formed by inversion in XSB3.

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.

Improving Illumina assemblies with Hi-C and long reads: An example with the North African dromedary

Researchers have assembled thousands of eukaryotic genomes using Illumina reads, but traditional mate‐pair libraries cannot span all repetitive elements, resulting in highly fragmented assemblies. However, both chromosome conformation capture techniques, such as Hi‐C and Dovetail Genomics Chicago libraries and long‐read sequencing, such as Pacific Biosciences and Oxford Nanopore, help span and resolve repetitive regions and therefore improve genome assemblies. One important livestock species of arid regions that does not have a high‐quality contiguous reference genome is the dromedary (Camelus dromedarius). Draft genomes exist but are highly fragmented, and a high‐quality reference genome is needed to understand adaptation to desert environments and artificial selection during domestication. Dromedaries are among the last livestock species to have been domesticated, and together with wild and domestic Bactrian camels, they are the only representatives of the Camelini tribe, which highlights their evolutionary significance. Here we describe our efforts to improve the North African dromedary genome. We used Chicago and Hi‐C sequencing libraries from Dovetail Genomics to resolve the order of previously assembled contigs, producing almost chromosome‐level scaffolds. Remaining gaps were filled with Pacific Biosciences long reads, and then scaffolds were comparatively mapped to chromosomes. Long reads added 99.32 Mbp to the total length of the new assembly. Dovetail Chicago and Hi‐C libraries increased the longest scaffold over 12‐fold, from 9.71 Mbp to 124.99 Mbp and the scaffold N50 over 50‐fold, from 1.48 Mbp to 75.02 Mbp. We demonstrate that Illumina de novo assemblies can be substantially upgraded by combining chromosome conformation capture and long‐read sequencing.

Evaluation of SNP Genotyping in Alpacas Using the Bovine HD Genotyping Beadchip

Alpacas are one of four South American Camelid species living in the highlands of the Andes. Production of alpaca fiber contributes to the economy of the region and the livelihood of many rural families. Fiber quantity and quality are important and in need of a modern breeding program based on genomic selection to accelerate genetic gain. To achieve this is necessary to discover enough molecular markers, single nucleotide polymorphisms (SNPs) in particular, to provide genome coverage and facilitate genome wide association studies to fiber production characteristics. The aim of this study was to discover alpaca SNPs by genotyping forty alpaca DNA samples using the BovineHD Genotyping Beadchip. Data analysis was performed with GenomeStudio (Illumina) software. Because different filters and thresholds are reported in the literature we investigated the effects of no-call threshold (≥0.05, ≥0.15, and ≥0.25) and call frequency (≥0.9 and =1.0) in identifying positive SNPs. Average GC Scores, calculated as the average of the 10% and 50% GenCall scores for each SNP (≥0.70) and the GenTrain score ≥ 0.25 parameters were applied to all comparisons. SNPs with minor allele frequency (MAF) ≥ 0.05 or ≥ 0.01 were retained. Since detection of SNPs is based on the stable binding of oligonucleotide probes to the target DNA immediately adjacent to the variant nucleotide, all positive SNP flanking sequences showing perfect alignments between the bovine and alpaca genomes for the first 21 or 26 nucleotides flanking the variant nucleotide at either side were selected. Only SNPs localized in one scaffold were assumed unique. Unique SNPs identified in both reference genomes were kept and mapped on the Vicugna_pacos 2.0.2 genome. The effects of the no-call threshold ≥ 0.25, call frequency = 1 and average GC ≥ 0.7 were meaningful and identified 6756 SNPs of which 400 were unique and polymorphic (MAF ≥ 0.01). Assignment to alpaca chromosomes was possible for 292 SNPs. Likewise, 209 SNPs were localized in 202 alpaca gene loci and 29 of these share the same loci with the dromedary. Interestingly, 69 of 400 alpaca SNPs have 100% similarity with dromedary.

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.

Evolution of gene regulation in ruminants differs between evolutionary breakpoint regions and homologous synteny blocks

The role of chromosome rearrangements in driving evolution has been a long-standing question of evolutionary biology. Here we focused on ruminants as a model to assess how rearrangements may have contributed to the evolution of gene regulation. Using reconstructed ancestral karyotypes of Cetartiodactyls, Ruminants, Pecorans, and Bovids, we traced patterns of gross chromosome changes. We found that the lineage leading to the ruminant ancestor after the split from other cetartiodactyls was characterized by mostly intrachromosomal changes, whereas the lineage leading to the pecoran ancestor (including all livestock ruminants) included multiple interchromosomal changes. We observed that the liver cell putative enhancers in the ruminant evolutionary breakpoint regions are highly enriched for DNA sequences under selective constraint acting on lineage-specific transposable elements (TEs) and a set of 25 specific transcription factor (TF) binding motifs associated with recently active TEs. Coupled with gene expression data, we found that genes near ruminant breakpoint regions exhibit more divergent expression profiles among species, particularly in cattle, which is consistent with the phylogenetic origin of these breakpoint regions. This divergence was significantly greater in genes with enhancers that contain at least one of the 25 specific TF binding motifs and located near bovidae-to-cattle lineage breakpoint regions. Taken together, by combining ancestral karyotype reconstructions with analysis of cis regulatory element and gene expression evolution, our work demonstrated that lineage-specific regulatory elements colocalized with gross chromosome rearrangements may have provided valuable functional modifications that helped to shape ruminant evolution.

A Near Chromosome Assembly of the Dromedary Camel Genome

The dromedary camel is an economically and socially important species of livestock in many parts of the world, being used for transport and the production of milk and meat. Much like cattle and horses, the camel may be found in industrial farming conditions as well as used in sporting. Camel racing is a multi-million dollar industry, with some specimens being valued at upward of 9.5 million USD. Despite its apparent value to humans, the dromedary camel is a neglected species in genomics. While cattle and other domesticated species have had much attention in terms of genome assembly, the camel has only been assembled to scaffold level, which does not give a clear indication of the order or chromosomal location of sequenced fragments. In this study, the Reference Assistant Chromosome Assembly (RACA) algorithm was implemented to use read-pair information of camel scaffolds, aligned with the cattle and human genomes in order to organize and orient these scaffolds in a near-chromosome level assembly. This method generated 72 large size fragments (N50 54.36 Mb). These predicted chromosome fragments (PCFs) were then compared with comparative maps of camel and cytogenetic map of alpaca chromosomes, allowing us to further upgrade the assembly. This dromedary camel assembly will be an invaluable tool to verify future camel assemblies generated with chromatin conformation or/and long read technologies. This study provides the first near-chromosome assembly of the dromedary camel, thus adding this economically important species to a growing pool of knowledge regarding the genome structure of domesticated livestock.

Emerging patterns of genome organization in Notopteridae species (Teleostei, Osteoglossiformes) as revealed by Zoo-FISH and Comparative Genomic Hybridization (CGH)

Notopteridae (Teleostei, Osteoglossiformes) represents an old fish lineage with ten currently recognized species distributed in African and Southeastern Asian rivers. Their karyotype structures and diploid numbers remained conserved over long evolutionary periods, since African and Asian lineages diverged approximately 120 Mya. However, a significant genetic diversity was already identified for these species using molecular data. Thus, why the evolutionary relationships within Notopteridae are so diverse at the genomic level but so conserved in terms of their karyotypes? In an attempt to develop a more comprehensive picture of the karyotype and genome evolution in Notopteridae, we performed comparative genomic hybridization (CGH) and cross-species (Zoo-FISH) whole chromosome painting experiments to explore chromosome-scale intergenomic divergence among seven notopterid species, collected in different African and Southeast Asian river basins. CGH demonstrated an advanced stage of sequence divergence among the species and Zoo-FISH experiments showed diffuse and limited homology on inter-generic level, showing a temporal reduction of evolutionarily conserved syntenic regions. The sharing of a conserved chromosomal region revealed by Zoo-FISH in these species provides perspectives that several other homologous syntenic regions have remained conserved among their genomes despite long temporal isolation. In summary, Notopteridae is an interesting model for tracking the chromosome evolution as it is (i) ancestral vertebrate group with Gondwanan distribution and (ii) an example of animal group exhibiting karyotype stasis. The present study brings new insights into degree of genome divergence vs. conservation at chromosomal and sub-chromosomal level in representative sampling of this group.

ATM orchestrates the DNA-damage response to counter toxic non-homologous end-joining at broken replication forks

Mutations in the ATM tumor suppressor gene confer hypersensitivity to DNA-damaging chemotherapeutic agents. To explore genetic resistance mechanisms, we performed genome-wide CRISPR-Cas9 screens in cells treated with the DNA topoisomerase I inhibitor topotecan. Thus, we here establish that inactivating terminal components of the non-homologous end-joining (NHEJ) machinery or of the BRCA1-A complex specifically confer topotecan resistance to ATM-deficient cells. We show that hypersensitivity of ATM-mutant cells to topotecan or the poly-(ADP-ribose) polymerase (PARP) inhibitor olaparib reflects delayed engagement of homologous recombination at DNA-replication-fork associated single-ended double-strand breaks (DSBs), allowing some to be subject to toxic NHEJ. Preventing DSB ligation by NHEJ, or enhancing homologous recombination by BRCA1-A complex disruption, suppresses this toxicity, highlighting a crucial role for ATM in preventing toxic LIG4-mediated chromosome fusions. Notably, suppressor mutations in ATM-mutant backgrounds are different to those in BRCA1-mutant scenarios, suggesting new opportunities for patient stratification and additional therapeutic vulnerabilities for clinical exploitation.

High-quality whole-genome sequence of an abundant Holarctic odontocete, the harbour porpoise (Phocoena phocoena)

The harbour porpoise (Phocoena phocoena) is a highly mobile cetacean found across the Northern hemisphere. It occurs in coastal waters and inhabits basins that vary broadly in salinity, temperature and food availability. These diverse habitats could drive subtle differentiation among populations, but examination of this would be best conducted with a robust reference genome. Here, we report the first harbour porpoise genome, assembled de novo from an individual originating in the Kattegat Sea (Sweden). The genome is one of the most complete cetacean genomes currently available, with a total size of 2.39 Gb and 50% of the total length found in just 34 scaffolds. Using 122 of the longest scaffolds, we were able to show high levels of synteny with the genome of the domestic cattle (Bos taurus). Our draft annotation comprises 22,154 predicted genes, which we further annotated through matches to the NCBI nucleotide database, GO categorization and motif prediction. Within the predicted genes, we have confirmed the presence of >20 genes or gene families that have been associated with adaptive evolution in other cetaceans. Overall, this genome assembly and draft annotation represent a crucial addition to the genomic resources currently available for the study of porpoises and Phocoenidae evolution, phylogeny and conservation.

Genetic characterization of the oxytocin-neurophysin I gene (OXT) and its regulatory regions analysis in domestic Old and New World camelids

Oxytocin is a neurohypophysial peptide linked to a wide range of biological functions, including milk ejection, temperament and reproduction. Aims of the present study were a) the characterization of the OXT (Oxytocin-neurophysin I) gene and its regulatory regions in Old and New world camelids; b) the investigation of the genetic diversity and the discovery of markers potentially affecting the gene regulation. On average, the gene extends over 814 bp, ranging between 825 bp in dromedary, 811 bp in Bactrian and 810 bp in llama and alpaca. Such difference in size is due to a duplication event of 21 bp in dromedary. The main regulatory elements, including the composite hormone response elements (CHREs), were identified in the promoter, whereas the presence of mature microRNAs binding sequences in the 3'UTR improves the knowledge on the factors putatively involved in the OXT gene regulation, although their specific biological effect needs to be still elucidated. The sequencing of genomic DNA allowed the identification of 17 intraspecific polymorphisms and 69 nucleotide differences among the four species. One of these (MF464535:g.622C>G) is responsible, in alpaca, for the loss of a consensus sequence for the transcription factor SP1. Furthermore, the same SNP falls within a CpG island and it creates a new methylation site, thus opening future possibilities of investigation to verify the influence of the novel allelic variant in the OXT gene regulation. A PCR-RFLP method was setup for the genotyping and the frequency of the allele C was 0.93 in a population of 71 alpacas. The obtained data clarify the structure of OXT gene in domestic camelids and add knowledge to the genetic variability of a genomic region, which has received little investigation so far. These findings open the opportunity for new investigations, including association studies with productive and reproductive traits.

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.

First Chromosomal Analysis in Hepsetidae (Actinopterygii, Characiformes): Insights into Relationship between African and Neotropical Fish Groups

Hepsetidae is a small fish family with only the genus Hepsetus, with six described species distributed throughout the South, Central and Western regions of Africa, showing a close relationship with the Alestidae and some Neotropical fish families. However, no cytogenetic information is available for both Hepsetidae and Alestidae species, thus preventing any evolutionary comparative studies at the chromosomal level. In the present study, we are providing new cytogenetic data for Hepsetus odoe, including the standard karyotype, C-banding, repetitive DNAs mapping, comparative genomic hybridization (CGH) and whole chromosome painting (WCP), providing chromosomal patterns and subsidies for comparative cytogenetics with other characiform families. Both males and females H. odoe have 2n = 58 chromosomes (10m + 28sm + 20st/a), with most of the C-band positive heterochromatin localized in the centromeric and subtelomeric regions. Only one pair of chromosomes bears proximal 5S rDNA sites in the short arms, contrasting with the 18S rDNA sequences which are located in the terminal regions of four chromosome pairs. Clear interstitial hybridization signals are evidenced for the U1 and U2 snDNA probes, but in only one and two chromosome pairs, respectively. Microsatellite motifs are widely distributed in the karyotype, with exception for the (CGG)₁₀, (GAA)₁₀ and (GAG)₁₀ probes, which highlight conspicuous interstitial signals on an unique pair of chromosomes. Comparative data from conventional and molecular cytogenetics, including CGH and WCP experiments, indicate that H. odoe and some Erythrinidae species, particularly Erythrinus erythrinus, share similar chromosomal sequences suggesting some relatedness among them, although bearing genomic specificities in view of their divergent evolutionary histories.

X Chromosome Evolution in Cetartiodactyla

The phenomenon of a remarkable conservation of the X chromosome in eutherian mammals has been first described by Susumu Ohno in 1964. A notable exception is the cetartiodactyl X chromosome, which varies widely in morphology and G-banding pattern between species. It is hypothesized that this sex chromosome has undergone multiple rearrangements that changed the centromere position and the order of syntenic segments over the last 80 million years of Cetartiodactyla speciation. To investigate its evolution we have selected 26 evolutionarily conserved bacterial artificial chromosome (BAC) clones from the cattle CHORI-240 library evenly distributed along the cattle X chromosome. High-resolution BAC maps of the X chromosome on a representative range of cetartiodactyl species from different branches: pig (Suidae), alpaca (Camelidae), gray whale (Cetacea), hippopotamus (Hippopotamidae), Java mouse-deer (Tragulidae), pronghorn (Antilocapridae), Siberian musk deer (Moschidae), and giraffe (Giraffidae) were obtained by fluorescent in situ hybridization. To trace the X chromosome evolution during fast radiation in specious families, we performed mapping in several cervids (moose, Siberian roe deer, fallow deer, and Pere David's deer) and bovid (muskox, goat, sheep, sable antelope, and cattle) species. We have identified three major conserved synteny blocks and rearrangements in different cetartiodactyl lineages and found that the recently described phenomenon of the evolutionary new centromere emergence has taken place in the X chromosome evolution of Cetartiodactyla at least five times. We propose the structure of the putative ancestral cetartiodactyl X chromosome by reconstructing the order of syntenic segments and centromere position for key groups.

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.

Comparative Chromosome Map and Heterochromatin Features of the Gray Whale Karyotype (Cetacea)

Cetacean karyotypes possess exceptionally stable diploid numbers and highly conserved chromosomes. To date, only toothed whales (Odontoceti) have been analyzed by comparative chromosome painting. Here, we studied the karyotype of a representative of baleen whales, the gray whale (Eschrichtius robustus, Mysticeti), by Zoo-FISH with dromedary camel and human chromosome-specific probes. We confirmed a high degree of karyotype conservation and found an identical order of syntenic segments in both branches of cetaceans. Yet, whale chromosomes harbor variable heterochromatic regions constituting up to a third of the genome due to the presence of several types of repeats. To investigate the cause of this variability, several classes of repeated DNA sequences were mapped onto chromosomes of whale species from both Mysticeti and Odontoceti. We uncovered extensive intrapopulation variability in the size of heterochromatic blocks present in homologous chromosomes among 3 individuals of the gray whale by 2-step differential chromosome staining. We show that some of the heteromorphisms observed in the gray whale karyotype are due to distinct amplification of a complex of common cetacean repeat and heavy satellite repeat on homologous autosomes. Furthermore, we demonstrate localization of the telomeric repeat in the heterochromatin of both gray and pilot whale (Globicephala melas, Odontoceti). Heterochromatic blocks in the pilot whale represent a composite of telomeric and common repeats, while heavy satellite repeat is lacking in the toothed whale consistent with previous studies.

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.

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.

Sequential cross-species chromosome painting among river buffalo, cattle, sheep and goat: a useful tool for chromosome abnormalities diagnosis within the family Bovidae

The main goal of this study was to develop a comparative multi-colour Zoo-FISH on domestic ruminants metaphases using a combination of whole chromosome and sub-chromosomal painting probes obtained from the river buffalo species (Bubalus bubalis, 2n = 50,XY). A total of 13 DNA probes were obtained through chromosome microdissection and DOP-PCR amplification, labelled with two fluorochromes and sequentially hybridized on river buffalo, cattle (Bos taurus, 2n = 60,XY), sheep (Ovis aries, 2n = 54,XY) and goat (Capra hircus, 2n = 60,XY) metaphases. The same set of paintings were then hybridized on bovine secondary oocytes to test their potential use for aneuploidy detection during in vitro maturation. FISH showed excellent specificity on metaphases and interphase nuclei of all the investigated species. Eight pairs of chromosomes were simultaneously identified in buffalo, whereas the same set of probes covered 13 out 30 chromosome pairs in the bovine and goat karyotypes and 40% of the sheep karyotype (11 out of 27 chromosome pairs). This result allowed development of the first comparative M-FISH karyotype within the domestic ruminants. The molecular resolution of complex karyotypes by FISH is particularly useful for the small chromosomes, whose similarity in the banding patterns makes their identification very difficult. The M-FISH karyotype also represents a practical tool for structural and numerical chromosome abnormalities diagnosis. In this regard, the successful hybridization on bovine secondary oocytes confirmed the potential use of this set of probes for the simultaneous identification on the same germ cell of 12 chromosome aneuploidies. This is a fundamental result for monitoring the reproductive health of the domestic animals in relation to management errors and/or environmental hazards.

Comparative chromosome painting of pronghorn (Antilocapra americana) and saola (Pseudoryx nghetinhensis) karyotypes with human and dromedary camel probes

BACKGROUND

Pronghorn (Antilocapridae, 2n = 58) and saola (Bovidae, 2n = 50) are members of Pecora, a highly diversified group of even-toed hoofed mammals. Karyotypes of these species were not involved in chromosome painting studies despite their intriguing phylogenetic positions in Pecora.

RESULTS

To trace the chromosome evolution during very fast radiation of main families from the common Pecoran ancestor, high-resolution comparative chromosome maps of pronghorn and saola with human (HSA) and dromedary camel (CDR) painting probes were established. The human and dromedary camel painting probes revealed 50 and 64 conserved segments respectively in the pronghorn genome, while 51 and 63 conserved segments respectively in the saola genome. Integrative analysis with published comparative maps showed that inversions in chromosomes homologous to CDR19/35/19 (HSA 10/20/10), CDR12/34/12 (HSA12/22/12/22), CDR10/33/10 (HSA 11) are present in representatives of all five living Pecoran families. The pronghorn karyotype could have formed from a putative 2n = 58 Pecoran ancestral karyotype by one fission and one fusion and that the saola karyotype differs from the presumed 2n = 60 bovid ancestral karyotype (2n = 60) by five fusions.

CONCLUSION

The establishment of high-resolution comparative maps for pronghorn and saola has shed some new insights into the putative ancestral karyotype, chromosomal evolution and phylogenic relationships in Pecora. No cytogenetic signature rearrangements were found that could unite the Antilocapridae with Giraffidae or with any other Pecoran families. Our data on the saola support a separate position of Pseudorigyna subtribe rather than its affinity to either Bovina or Bubalina, but the saola phylogenetic position within Bovidae remains unresolved.

The utility of chromosome microdissection in clinical cytogenetics: a new reciprocal translocation in sheep

Local sheep breeders and scientists in Italy cooperate and conduct research on the genetic improvement of autochthonous genetic types (AGTs) by various approaches, including a cytogenetic breeding selection since 2011. The Laticauda sheep (Ovis aries, 2n = 54) breed is one of the AGTs reared in the Campania region (southern Italy). Performing cytogenetic analyses, we have detected and described a novel reciprocal translocation in a Laticauda sheep identified as 54,XX t(18;23)(q14;q26). Our data support recurring appeals that suggest the regular performance of cytogenetic analyses for monitoring genetic health of livestock species. In total, 5 cases of reciprocal translocations in sheep are known, including the new case. None of them has any phenotypic effect on the living offspring. However, affected animals are characterized by sterility or have a low fertility which can have an effect on breeding success and on economical balance. Presence and kind of the described novel chromosomal aberration were detected by performing CBA-banding and FISH mapping with telomeric probes. RBA-banding allowed the karyotyping of sheep chromosomes and the identification of aberrant chromosomes and regions involved in the new reciprocal translocation. Whole chromosome painting (WCP) probes received from equivalent chromosomes in cattle and the derivative sheep chromosome 18 confirmed the cytogenetic data. This way, our study underlined both the importance of WCP probes by chromosome microdissection and a new way to use WCP probes directly generated from derivative chromosomes.

Expression of pregnancy-associated glycoprotein family in the epitheliochorial placenta of two Camelidae species (C. dromedarius and C. bactrianus)

This study describes placental morphology and immunolocalization of the placental pregnancy associated glycoprotein-like family (PAGs) identified in two selected taxa of Old-World camels of the Camelidae family: Camelus dromedarius (Cd) and Camelus bactrianus (Cb). Placental tissues of Cd from days 140-293 post-coitum (dpc), term (404 dpc); and of Cb from term (440 dpc) were examined. Histological staining (hematoxylin/eosin and propidium iodine) revealed the development of the placental structure, while chorionic folding increased the feto-placental surface during the progress of pregnancy. The camelid placenta during early pregnancy is similar to the diffuse epitheliochorial type, and during later stages of pregnancy resembles the synepitheliochorial (cotyledonary) type. Placental expression of the PAGs was detected (Alexa 488 - green) within camelid trophectoderm cells (TRD - chorionic epithelium as outer layer of embryonic cells) among all placental cells with nuclei stained by propidium iodide (red). The PAGs, identified in both Camelidae taxa, were named CbPAGs and CdPAGs. Placental CbPAG and CdPAG expression is restricted to the TRD cells, which are differentially developed throughout gestation. Cross-reactivity of polyvalent anti-pPAG polyclonals with the CbPAGs and CdPAGs revealed high structural similarities of the PAG-like epitopes in pigs and camels. This is the first study identifying PAG expression in chorionic cells of the camel placenta.

Development and application of camelid molecular cytogenetic tools

Cytogenetic chromosome maps offer molecular tools for genome analysis and clinical cytogenetics and are of particular importance for species with difficult karyotypes, such as camelids (2n = 74). Building on the available human-camel zoo-fluorescence in situ hybridization (FISH) data, we developed the first cytogenetic map for the alpaca (Lama pacos, LPA) genome by isolating and identifying 151 alpaca bacterial artificial chromosome (BAC) clones corresponding to 44 specific genes. The genes were mapped by FISH to 31 alpaca autosomes and the sex chromosomes; 11 chromosomes had 2 markers, which were ordered by dual-color FISH. The STS gene mapped to Xpter/Ypter, demarcating the pseudoautosomal region, whereas no markers were assigned to chromosomes 14, 21, 22, 28, and 36. The chromosome-specific markers were applied in clinical cytogenetics to identify LPA20, the major histocompatibility complex (MHC)-carrying chromosome, as a part of an autosomal translocation in a sterile male llama (Lama glama, LGL; 2n = 73,XY). FISH with LPAX BACs and LPA36 paints, as well as comparative genomic hybridization, were also used to investigate the origin of the minute chromosome, an abnormally small LPA36 in infertile female alpacas. This collection of cytogenetically mapped markers represents a new tool for camelid clinical cytogenetics and has applications for the improvement of the alpaca genome map and sequence assembly.

Ovarian dysgenesis in an alpaca with a minute chromosome 36

A 4-year-old female alpaca (Lama pacos [LPA]) was presented to the Oregon State Veterinary Teaching Hospital for failure to display receptive behavior to males. Although no abnormalities were found on physical examination, transrectal ultrasonographic examination of the reproductive tract revealed uterine hypoplasia and ovarian dysgenesis. Cytogenetic analysis demonstrated a normal female 74,XX karyotype with 1 exceptionally small (minute) homologue of autosome LPA36. Chromosome analysis by Giemsa staining and DAPI- and C-banding revealed that the minute LPA36 was submetacentric, AT-rich, and largely heterochromatic. Because of the small size and lack of molecular markers, it was not possible to identify the origin of the minute. There is a need to improve molecular cytogenetic tools to further study the phenomenon of this minute chromosome and its relation to female reproduction in alpacas and llamas.

Chromosome evolution in Perissodactyla

Comparative painting has provided a wealth of useful information and helped to reconstruct the pathways of karyotype evolution within major eutherian phylogenetic clades. New data have come from gene localizations, BAC mapping and high throughout sequencing projects that enrich and provide new details of genome evolution. Extensive research on perissodactyl genomes has revealed not only increased rates of chromosomal rearrangements, but also an exceptionally high number of centromere repositioning events in equids. Here were combined new physical mapping, comparative painting and genome sequencing data to refine the putative ancestral karyotype maps and to revise the previously proposed scenario of perissodactyl karyotype evolution.

Chromosome evolution in bats as revealed by FISH: the ongoing search for the ancestral chiropteran karyotype

Chiroptera, the second largest order of mammals, comprises more than 1,000 species in 18 highly morphologically diverse families. Chromosome painting with human probes has been applied to 10 bat species from 8 families. Except for the combination 10/12pq/22q, all syntenic segmental associations proposed for the mammalian ancestor have been found in Chiroptera. Bat-specific painting probes, established from 4 species of 3 families, have been used in whole chromosome painting experiments in 29 species from 8 families. The results show that the prevailing mode of chromosomal evolution in bats is Robertsonian translocation with a large number of convergent events. Given our present knowledge of chiropteran karyotypes, only a few elements of the ancestral chiropteran karyotype can be reconstructed with confidence.

Comparative molecular cytogenetics in Cetartiodactyla

Cetartiodactyla comprises Artiodactyla (even-toed ungulates) and Cetacea (whales, dolphins and porpoises). Artiodactyla is a large taxon represented by about 200 living species ranked in 10 families. Cetacea are classified into 13 families with almost 80 species. Many publications concerning karyotypic relationships in Cetartiodactyla have been published in previous decades. Formerly, the karyotypes of closely related species were compared by chromosome banding. Introduction of molecular cytogenetic methods facilitated comparative mapping between species with highly rearranged karyotypes and distantly related species. Such information is a prerequisite for the understanding of karyotypic phylogeny and the reconstruction of the karyotypes of common ancestors. This study summarizes the data on chromosome evolution in Cetartiodactyla, mainly derived from molecular cytogenetic studies. Traditionally, phylogenetic relationships of most groups have been estimated using morphological data. However, the results of some molecular studies of mammalian phylogeny are discordant with traditional conceptions of phylogeny. Cetartiodactyls provide several examples of incongruence between traditional morphological and molecular data. Such cases of conflict include the relationships of the major clades of artiodactyls, the relationships among the extant families of the suborder Ruminantia or the phylogeny of the family Bovidae. The most unexpected aspect of the molecular phylogeny was the recognition that Cetacea is a deeply nested member of Artiodactyla. The largest living order of terrestrial hoofed mammals is the even-toed hoofed mammals, or Artiodactyla. The artiodactyls are composed of over 190 living species including pigs, peccaries, hippos, camels, llamas, deer, pronghorns, giraffes, sheep, goats, cattle and antelopes. Cetacea is an order of wholly aquatic mammals, which include whales, dolphins and porpoises. Cetartiodactyla has become the generally accepted name for the clade containing both of these orders.

Chromosome painting in Tragulidae facilitates the reconstruction of Ruminantia ancestral karyotype

Although Tragulidae, as the basal family in Ruminantia phylogenetic tree, is the key taxon for understanding the early chromosome evolution of extant ruminants, comparative molecular cytogenetic data on the tragulids are scarce. Here, we present the first genome-wide comparative map of the Java mouse deer (Tragulus javanicus, Tragulidae) revealed by chromosome painting with human and dromedary probes. Together with the published comparative maps of major representative cetartiodactyl species established with the same set of probes, our results allowed us to reconstruct a 2n = 48 Ruminantia ancestral karyotype, which is similar to the cetartiodactyl ancestral karyotype. The karyotype evolution of T. javanicus has involved multiple rearrangements, most of which appear to be apomorphic and have not found in karyotype evolution of pecoran species (i.e., Ruminantia excluding Tragulidae). The rate of chromosome evolution of the mouse deer was rather low-0.4 R/Ma, while the estimated tempo of chromosome changes on the lineages leading from Cetartiodactyla ancestor to Ruminantia and from Ruminantia to Pecora were roughly the same (about 1.2 R/Ma).

Restriction of porcine endogenous retrovirus by porcine APOBEC3 cytidine deaminases

Xenotransplantation of porcine cells, tissues, and organs shows promise to surmount the shortage of human donor materials. Among the barriers to pig-to-human xenotransplantation are porcine endogenous retroviruses (PERV) since functional representatives of the two polytropic classes, PERV-A and PERV-B, are able to infect human embryonic kidney cells in vitro, suggesting that a xenozoonosis in vivo could occur. To assess the capacity of human and porcine cells to counteract PERV infections, we analyzed human and porcine APOBEC3 (A3) proteins. This multigene family of cytidine deaminases contributes to the cellular intrinsic immunity and act as potent inhibitors of retroviruses and retrotransposons. Our data show that the porcine A3 gene locus on chromosome 5 consists of the two single-domain genes A3Z2 and A3Z3. The evolutionary relationships of the A3Z3 genes reflect the evolutionary history of mammals. The two A3 genes encode at least four different mRNAs: A3Z2, A3Z3, A3Z2-Z3, and A3Z2-Z3 splice variant A (SVA). Porcine and human A3s have been tested toward their antiretroviral activity against PERV and murine leukemia virus (MuLV) using novel single-round reporter viruses. The porcine A3Z2, A3Z3 and A3Z2-Z3 were packaged into PERV particles and inhibited PERV replication in a dose-dependent manner. The antiretroviral effect correlated with editing by the porcine A3s with a trinucleotide preference for 5' TGC for A3Z2 and A3Z2-Z3 and 5' CAC for A3Z3. These results strongly imply that human and porcine A3s could inhibit PERV replication in vivo, thereby reducing the risk of infection of human cells by PERV in the context of pig-to-human xenotransplantation.