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

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.

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.

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.