Genomic analysis of the domestication and post-Spanish conquest evolution of the llama and alpaca

Genetic legacy and adaptive signatures: investigating the history, diversity, and selection signatures in Rendena cattle resilient to eighteenth century rinderpest epidemics

BACKGROUND

Rendena is a dual-purpose cattle breed, which is primarily found in the Italian Alps and the eastern areas of the Po valley, and recognized for its longevity, fertility, disease resistance and adaptability to steep Alpine pastures. It is categorized as 'vulnerable to extinction' with only 6057 registered animals in 2022, yet no comprehensive analyses of its molecular diversity have been performed to date. The aim of this study was to analyse the origin, genetic diversity, and genomic signatures of selection in Rendena cattle using data from samples collected in 2000 and 2018, and shed light on the breed's evolution and conservation needs.

RESULTS

Genetic analysis revealed that the Rendena breed shares genetic components with various Alpine and Po valley breeds, with a marked genetic proximity to the Original Braunvieh breed, reflecting historical restocking efforts across the region. The breed shows signatures of selection related to both milk and meat production, environmental adaptation and immune response, the latter being possibly the result of multiple rinderpest epidemics that swept across the Alps in the eighteenth century. An analysis of the Rendena cattle population spanning 18 years showed an increase in the mean level of inbreeding over time, which is confirmed by the mean number of runs of homozygosity per individual, which was larger in the 2018 sample.

CONCLUSIONS

The Rendena breed, while sharing a common origin with Brown Swiss, has developed distinct traits that enable it to thrive in the Alpine environment and make it highly valued by local farmers. Preserving these adaptive features is essential, not only for maintaining genetic diversity and enhancing the ability of this traditional animal husbandry to adapt to changing environments, but also for guaranteeing the resilience and sustainability of both this livestock system and the livelihoods within the Rendena valley.

History of Diversification and Adaptation from North to South Revealed by Genomic Data: Guanacos from the Desert to Sub-Antarctica

The increased availability of quality genomic data has greatly improved the scope and resolution of our understanding of the recent evolutionary history of wild species adapted to extreme environments and their susceptibility to anthropogenic impacts. The guanaco (Lama guanicoe), the largest wild ungulate in South America, is a good example. The guanaco is well adapted to a wide range of habitats, including the Sechura Desert, the high Andes Mountains to the north, and the extreme temperatures and conditions of Navarino Island to the south. Guanacos also have a long history of overexploitation by humans. To assess the evolutionary impact of these challenging habitats on the genomic diversity, we analyzed 38 genomes (∼10 to 16×) throughout their extensive latitudinal distribution from the Sechura and Atacama Desert to southward into Tierra del Fuego Island. These included analyses of patterns of unique differentiation in the north and geographic region further south with admixture among L. g. cacsilensis and L. g. guanicoe. Our findings provide new insights on the divergence of the subspecies ∼800,000 yr BP and document two divergent demographic trajectories and to the initial expansion of guanaco into the more southern portions of the Atacama Desert. Patagonian guanacos have experienced contemporary reductions in effective population sizes, likely the consequence of anthropogenic impacts. The lowest levels of genetic diversity corresponded to their northern and western limits of distribution and some varying degrees of genetic differentiation. Adaptive genomic diversity was strongly linked with environmental variables and was linked with colonization toward the south followed by adaptation.

Conflict between Farmers and Guanacos (Lama guanicoe cacsilensis): Field Surveys, Remote Sensing, and Interviews Provide Information for Conservation of a Critically Endangered Species in Southern Peru

The Peruvian guanaco (Lama guanicoe cacsilensis) is classified as being "in critical danger of extinction" by the government. In this study, we evaluate how the conflict between farmers and guanacos in the Susapaya and Estique Districts, Tacna Department (Southern Peru) may represent a threat to their survival. To evaluate the situation, we 1. Conducted field surveys to monitor guanaco presence, 2. Used available remote sensing data to map guanaco movement, and 3. interviewed the impacted farmers concerning their losses. Remote sensing data showed that sedentary guanaco family groups located in prime steppe vegetation habitats never entered agricultural areas, while field surveys showed that bachelor bands and solitary individuals did, perhaps seeking forage due to growing population pressure. Interview data found that 90% of community farmers felt that guanacos were a problem best resolved by better fencing (45%), hunting (19%), or increased security (16%), and 92% saw no value in the conservation of the species. Hunting is illegal, given the critically endangered status of guanacos in Peru, so additional efforts are needed to both educate those who feel guanacos are a menace and involve them in efforts to preserve the species.

Whole genome sequencing analysis of alpaca suggests TRPV3 as a candidate gene for the suri phenotype

BACKGROUND

Alpaca is a domestic South American camelid probably arising from the domestication of two wild camelids, the vicugna and the guanaco. Two phenotypes are described for alpaca, known as huacaya and suri. Huacaya fleece is characterized by compact, soft, and highly crimped fibers, while suri fleece is longer, straight, less crimped, and lustrous. The gene variants determining these phenotypes are still unknown, although previous studies suggested a dominant inheritance of the suri. Based on that, the aim of this study was the identification of the gene variants determining alpaca coat phenotypes through whole genome sequencing (WGS) analysis.

RESULTS

The sample used includes two test-cross alpaca families, suri × huacaya, which produced two offspring, one with the suri phenotype and one with the huacaya phenotype. The analyzed sample was expanded through the addition of WGS data from six vicugnas and six guanacos; this because we assumed the absence of the gene variants linked to the suri phenotype in these wild species. The analysis of gene variant segregation with the suri phenotype, coupled with the filtering of gene variants present in the wild species, disclosed the presence in all the suri samples of a premature termination codon (PTC) in TRPV3 (transient receptor potential cation channel subfamily V member 3), a gene known to be involved in hair growth and cycling, thermal sensation, cold tolerance and adaptation in several species. Mutations in TRPV3 were previously associated with the alteration of hair structure leading to an impaired formation of the hair canal and the hair shaft in mouse. This PTC in TRPV3, due to a G > T substitution (p.Glu475*), results in a loss of 290 amino acids from the canonical translated protein, plausibly leading to a physiological dysfunction.

CONCLUSION

The present results suggest that the suri phenotype may arise from a TRPV3 gene variant which may explain some of the suri features such as its longer hair fibre with lower number of cuticular scales compared to huacaya.

Genetic variability among and within domestic Old and New World camels at the α-lactalbumin gene (LALBA) reveals new alleles and polymorphisms responsible for differential expression

α-Lactalbumin (α-LA), which is encoded by the LALBA gene, is a major whey protein that binds to Ca2+ and facilitates lactose synthesis as a regulatory subunit of the synthase enzyme complex. In addition, it has been shown to play central roles in immune modulation, cell-growth regulation, and antimicrobial activity. In this study, a multitechnical approach was used to fully characterize the LALBA gene and its variants in both coding and regulatory regions for domestic camelids (dromedary, Bactrian camel, alpaca, and llama). The gene analysis revealed a conserved structure among the camelids, but a slight difference in size (2,012 bp on average) due to intronic variations. Promoters were characterized for the transcription factor binding sites (11 found in total). Intraspecies sequence comparison showed 36 SNPs in total (2 in the dromedary, none in the Bactrian camel, 22 in the alpaca, and 12 in the llama), whereas interspecies comparison showed 86 additional polymorphic sites. Eight SNPs were identified as trans-specific polymorphisms, and 2 of them (g.112A>G and g.1229A>G) were particularly interesting in the New World camels. The first creates a new binding site for transcription factor SP1. An enhancing effect of the g.112G variant on the expression was demonstrated by 3 independent pGL3 gene reporter assays. The latter is responsible for the p.78Ile>Val AA replacement and represents novel allelic variants (named LALBA A and B). A link to protein variants has been established by isoelectric focusing (IEF), and bioinformatics analysis revealed that carriers of valine (g.1229G) have a higher glycosylation rate. Genotyping methods based on restriction fragment length polymorphism (PCR-RFLP) were set up for both SNPs. Overall, adenine was more frequent (0.54 and 0.76) at both loci. Four haplotypes were found, and the AA and GA were the most common with a frequency of 0.403 and 0.365, respectively. Conversely, a putative biological gain characterizes the haplotype GG. Therefore, opportunities for rapid directional selection can be realized if this haplotype is associated with favorable milk protein properties. This study adds knowledge at the gene and protein level for α-LA (LALBA) in camelids and importantly contributes to a relatively unexplored research area in these species.

Genome-wide scan for runs of homozygosity in South American Camelids

BACKGROUND

Alpaca (Vicugna pacos), llama (Lama glama), vicugna (Vicugna vicugna) and guanaco (Lama guanicoe), are the camelid species distributed over the Andean high-altitude grasslands, the Altiplano, and the Patagonian arid steppes. Despite the wide interest on these animals, most of the loci under selection are still unknown. Using whole-genome sequencing (WGS) data we investigated the occurrence and the distribution of Runs Of Homozygosity (ROHs) across the South American Camelids (SACs) genome to identify the genetic relationship between the four species and the potential signatures of selection.

RESULTS

A total of 37 WGS samples covering the four species was included in the final analysis. The multi-dimensional scaling approach showed a clear separation between the four species; however, admixture analysis suggested a strong genetic introgression from vicugna and llama to alpaca. Conversely, very low genetic admixture of the guanaco with the other SACs was found. The four species did not show significant differences in the number, length of ROHs (100-500 kb) and genomic inbreeding values. Longer ROHs (> 500 kb) were found almost exclusively in alpaca. Seven overlapping ROHs were shared by alpacas, encompassing nine loci (FGF5, LOC107034918, PRDM8, ANTXR2, LOC102534792, BSN, LOC116284892, DAG1 and RIC8B) while nine overlapping ROHs were found in llama with twenty-five loci annotated (ERC2, FZD9, BAZ1B, BCL7B, LOC116284208, TBL2, MLXIPL, PHF20, TRNAD-AUC, LOC116284365, RBM39, ARFGEF2, DCAF5, EXD2, HSPB11, LRRC42, LDLRAD1, TMEM59, LOC107033213, TCEANC2, LOC102545169, LOC116278408, SMIM15, NDUFAF2 and RCOR1). Four overlapping ROHs, with three annotated loci (DLG1, KAT6B and PDE4D) and three overlapping ROHs, with seven annotated genes (ATP6V1E1, BCL2L13, LOC116276952, BID, KAT6B, LOC116282667 and LOC107034552), were detected for vicugna and guanaco, respectively.

CONCLUSIONS

The signatures of selection revealed genomic areas potentially selected for production traits as well as for natural adaptation to harsh environment. Alpaca and llama hint a selection driven by environment as well as by farming purpose while vicugna and guanaco showed selection signals for adaptation to harsh environment. Interesting, signatures of selection on KAT6B gene were identified for both vicugna and guanaco, suggesting a positive effect on wild populations fitness. Such information may be of interest to further ecological and animal production studies.

Analysis of the genetic loci of pigment pattern evolution in vertebrates

Vertebrate pigmentation patterns are amongst the best characterised model systems for studying the genetic basis of adaptive evolution. The wealth of available data on the genetic basis for pigmentation evolution allows for analysis of trends and quantitative testing of evolutionary hypotheses. We employed Gephebase, a database of genetic variants associated with natural and domesticated trait variation, to examine trends in how cis-regulatory and coding mutations contribute to vertebrate pigmentation phenotypes, as well as factors that favour one mutation type over the other. We found that studies with lower ascertainment bias identified higher proportions of cis-regulatory mutations, and that cis-regulatory mutations were more common amongst animals harbouring a higher number of pigment cell classes. We classified pigmentation traits firstly according to their physiological basis and secondly according to whether they affect colour or pattern, and identified that carotenoid-based pigmentation and variation in pattern boundaries are preferentially associated with cis-regulatory change. We also classified genes according to their developmental, cellular, and molecular functions. We found a greater proportion of cis-regulatory mutations in genes implicated in upstream developmental processes compared to those involved in downstream cellular functions, and that ligands were associated with a higher proportion of cis-regulatory mutations than their respective receptors. Based on these trends, we discuss future directions for research in vertebrate pigmentation evolution.

A look into the wild. Pathological analysis of a modern collection of guanacos from the Dry Chaco and its implications for South American camelid paleopathological studies

OBJECTIVE

This article evaluates the prevalence of lesions in a modern osteological collection of guanacos (Lama guanicoe) and discusses the potential of paleopathological data to assess human intervention and environmental stress.

MATERIALS

A modern osteological collection of guanacos (NISP = 862) from north-western Córdoba, Central Argentina.

METHODS

The prevalence of pathological specimens per skeletal element and the pathological index (Bartosiewicz et al., 1997) was used. The prevalence of arthropathies, trauma and infections was quantified. Additionally, thorn lesions in the autopodium were recorded.

RESULTS

11.03 % of the specimens presented pathological changes and the mean pathological index was 0.01. Degenerative lesions were the most prevalent type (10.34 %), followed by traumatic (0.81 %) and infectious pathologies (0.12 %). Thorn lesions (2.55 %) were recorded especially in metapodials.

CONCLUSIONS

Guanacos are exposed to the development of degenerative lesions, mainly in the autopodium and vertebrae. These lesions are probably common in camelids and should not be used to argue human management. Traumatic and infectious lesions are less frequent.

SIGNIFICANCE

This work provides baseline information for the paleopathological study of South American camelids and contributes to the characterization of a regionally endangered species.

LIMITATIONS

The nature of the faunal assemblage did not allow for direct correlations between pathologies and individual variables such as sex or age.

SUGGESTIONS FOR FURTHER RESEARCH

The comparison of our results with other wild and domesticated modern populations would be valuable to expand the baseline information for paleopathological studies. The use of quantitative methods is encouraged for future comparative and diachronic studies.

Signatures of purifying selection and site-specific positive selection on the mitochondrial DNA of dromedary camels (Camelus dromedarius)

The two species of the Old World Camelini tribe, dromedary and Bactrian camels, show superior adaptability to the different environmental conditions they populate, e.g. desert, mountains and coastal areas, which might be associated with adaptive variations on their mitochondrial DNA. Here, we investigate signatures of natural selection in the 13-mitochondrial protein-coding genes of different dromedary camel populations from the Arabian Peninsula, Africa and southwest Asia. The full mitogenome sequences of 42 dromedaries, 38 domestic Bactrian, 29 wild Bactrian camels and 31 samples representing the New World Lamini tribe reveal species-wise genetic distinction among Camelidae family species, with no evidence of geographic distinction among dromedary camels. We observe gene-wide signals of adaptive divergence between the Old World and New World camels, with evidence of purifying selection among Old World camel species. Upon comparing the different Camelidae tribes, 27 amino acid substitutions across ten mtDNA protein-coding genes were found to be under positive selection, in which, 24 codons were defined to be under positive adaptive divergence between Old World and New World camels. Seven codons belonging to three genes demonstrated positive selection in dromedary lineage. A total of 89 codons were found to be under positive selection in Camelidae family based on investigating the impact of amino acid replacement on the physiochemical properties of proteins, including equilibrium constant and surrounding hydrophobicity. These mtDNA variants under positive selection in the Camelidae family might be associated with their adaptation to their contrasting environments.

Arctic introgression and chromatin regulation facilitated rapid Qinghai-Tibet Plateau colonization by an avian predator

The Qinghai-Tibet Plateau (QTP), possesses a climate as cold as that of the Arctic, and also presents uniquely low oxygen concentrations and intense ultraviolet (UV) radiation. QTP animals have adapted to these extreme conditions, but whether they obtained genetic variations from the Arctic during cold adaptation, and how genomic mutations in non-coding regions regulate gene expression under hypoxia and intense UV environment, remain largely unknown. Here, we assemble a high-quality saker falcon genome and resequence populations across Eurasia. We identify female-biased hybridization with Arctic gyrfalcons in the last glacial maximum, that endowed eastern sakers with alleles conveying larger body size and changes in fat metabolism, predisposing their QTP cold adaptation. We discover that QTP hypoxia and UV adaptations mainly involve independent changes in non-coding genomic variants. Our study highlights key roles of gene flow from Arctic relatives during QTP hypothermia adaptation, and cis-regulatory elements during hypoxic response and UV protection.

Genotyping-by-sequencing (GBS) as a tool for interspecies hybrid detection

Genotyping-by-sequencing (GBS) is an extremely useful, modern and relatively inexpensive approach to discovering high-quality single-nucleotide polymorphisms (SNPs), which seem to be the most promising markers for identifying hybrid individuals between different species, especially those that can create backcrosses. In addition, GBS could become an invaluable tool in finding backcrosses, even several generations back. Its potential for the use of restriction enzymes and species is almost unlimited. It can also be successfully applied to species for which a reference genome is not established. In this paper, we describe the GBS technique, its main advantages and disadvantages, and the research carried out using this method concerning interspecies hybridisation and the identification of fertile hybrids. We also present future approaches that could be of interest in the context of the GBS method.

Changes in Neuropeptide Prohormone Genes among Cetartiodactyla Livestock and Wild Species Associated with Evolution and Domestication

The impact of evolution and domestication processes on the sequences of neuropeptide prohormone genes that participate in cell–cell signaling influences multiple biological process that involve neuropeptide signaling. This information is important to understand the physiological differences between Cetartiodactyla domesticated species such as cow, pig, and llama and wild species such as hippopotamus, giraffes, and whales. Systematic analysis of changes associated with evolutionary and domestication forces in neuropeptide prohormone protein sequences that are processed into neuropeptides was undertaken. The genomes from 118 Cetartiodactyla genomes representing 22 families were mined for 98 neuropeptide prohormone genes. Compared to other Cetartiodactyla suborders, Ruminantia preserved PYY2 and lost RLN1. Changes in GNRH2, IAPP, INSL6, POMC, PRLH, and TAC4 protein sequences could result in the loss of some bioactive neuropeptides in some families. An evolutionary model suggested that most neuropeptide prohormone genes disfavor sequence changes that incorporate large and hydrophobic amino acids. A compelling finding was that differences between domestic and wild species are associated with the molecular system underlying ‘fight or flight’ responses. Overall, the results demonstrate the importance of simultaneously comparing the neuropeptide prohormone gene complement from close and distant-related species. These findings broaden the foundation for empirical studies about the function of the neuropeptidome associated with health, behavior, and food production.

Occurrence of sarcoptic mange in free-ranging vicuñas (Vicugna vicugna) of the Andean high plateau region of Argentina

Free-ranging vicuñas (Vicugna vicugna) are handled in some areas of the Andean high plateau region following an ancestral practice known as chaku, which consists in their transient capture and shearing of their fiber for commercialization. In this study, 807 vicuñas captured during 12 chaku events that took place in 2019 in the province of Jujuy, Argentina, were examined for typical mange skin lesions. Twenty-eight of the examined vicuñas presented alopecia with erythema, exudation, hyperkeratosis, and/or bleeding scarred lesions, mostly in the chest, rear and front legs, and inguinal zone. Most of the cases (82%) appeared in Laguna Cucho at 4900 masl, where 23% of the animals presented these skin reactions. Microscopic evaluation of skin scrapings revealed the presence of a great number of 0.1- to 0.4-mm-long mites of different life cycle stages, morphologically compatible with the species Sarcoptes scabiei. This etiological agent was confirmed by PCR amplification and sequencing of a cox-1 species-specific segment. Histopathological examination of skin biopsies showed extensive infiltration of the dermis with lymphocytes, neutrophils and eosinophils, hyperplasia at different stages, epidermis degeneration, and hyperkeratosis. This is the first characterization of sarcoptic mange in free-ranging vicuñas by clinical examination, mite morphology, histopathological studies, and molecular confirmation in the region. Mange hampers the welfare of vicuñas and the economy of the local communities that organize chaku events since infested vicuñas cannot be sheared. Its long-term effects are unknown but it might affect the fitness and survival of this iconic South American camelid.

Whole-genome sequence analysis unveils different origins of European and Asiatic mouflon and domestication-related genes in sheep

The domestication and subsequent development of sheep are crucial events in the history of human civilization and the agricultural revolution. However, the impact of interspecific introgression on the genomic regions under domestication and subsequent selection remains unclear. Here, we analyze the whole genomes of domestic sheep and their wild relative species. We found introgression from wild sheep such as the snow sheep and its American relatives (bighorn and thinhorn sheep) into urial, Asiatic and European mouflons. We observed independent events of adaptive introgression from wild sheep into the Asiatic and European mouflons, as well as shared introgressed regions from both snow sheep and argali into Asiatic mouflon before or during the domestication process. We revealed European mouflons might arise through hybridization events between a now extinct sheep in Europe and feral domesticated sheep around 6000-5000 years BP. We also unveiled later introgressions from wild sheep to their sympatric domestic sheep after domestication. Several of the introgression events contain loci with candidate domestication genes (e.g., PAPPA2, NR6A1, SH3GL3, RFX3 and CAMK4), associated with morphological, immune, reproduction or production traits (wool/meat/milk). We also detected introgression events that introduced genes related to nervous response (NEURL1), neurogenesis (PRUNE2), hearing ability (USH2A), and placental viability (PAG11 and PAG3) into domestic sheep and their ancestral wild species from other wild species.

The Quest for Genes Involved in Adaptation to Climate Change in Ruminant Livestock

Livestock radiated out from domestication centres to most regions of the world, gradually adapting to diverse environments, from very hot to sub-zero temperatures and from wet and humid conditions to deserts. The climate is changing; generally global temperature is increasing, although there are also more extreme cold periods, storms, and higher solar radiation. These changes impact livestock welfare and productivity. This review describes advances in the methodology for studying livestock genomes and the impact of the environment on animal production, giving examples of discoveries made. Sequencing livestock genomes has facilitated genome-wide association studies to localize genes controlling many traits, and population genetics has identified genomic regions under selection or introgressed from one breed into another to improve production or facilitate adaptation. Landscape genomics, which combines global positioning and genomics, has identified genomic features that enable animals to adapt to local environments. Combining the advances in genomics and methods for predicting changes in climate is generating an explosion of data which calls for innovations in the way big data sets are treated. Artificial intelligence and machine learning are now being used to study the interactions between the genome and the environment to identify historic effects on the genome and to model future scenarios.

Genomics of Adaptations in Ungulates

Ungulates are a group of hoofed animals that have long interacted with humans as essential sources of food, labor, clothing, and transportation. These consist of domesticated, feral, and wild species raised in a wide range of habitats and biomes. Given the diverse and extreme environments inhabited by ungulates, unique adaptive traits are fundamental for fitness. The documentation of genes that underlie their genomic signatures of selection is crucial in this regard. The increasing availability of advanced sequencing technologies has seen the rapid growth of ungulate genomic resources, which offers an exceptional opportunity to understand their adaptive evolution. Here, we summarize the current knowledge on evolutionary genetic signatures underlying the adaptations of ungulates to different habitats.

Ancient DNA reveals the lost domestication history of South American camelids in Northern Chile and across the Andes

The study of South American camelids and their domestication is a highly debated topic in zooarchaeology. Identifying the domestic species (alpaca and llama) in archaeological sites based solely on morphological data is challenging due to their similarity with respect to their wild ancestors. Using genetic methods also presents challenges due to the hybridization history of the domestic species, which are thought to have extensively hybridized following the Spanish conquest of South America that resulted in camelids slaughtered en masse. In this study, we generated mitochondrial genomes for 61 ancient South American camelids dated between 3,500 and 2,400 years before the present (Early Formative period) from two archaeological sites in Northern Chile (Tulán-54 and Tulán-85), as well as 66 modern camelid mitogenomes and 815 modern mitochondrial control region sequences from across South America. In addition, we performed osteometric analyses to differentiate big and small body size camelids. A comparative analysis of these data suggests that a substantial proportion of the ancient vicuña genetic variation has been lost since the Early Formative period, as it is not present in modern specimens. Moreover, we propose a domestication hypothesis that includes an ancient guanaco population that no longer exists. Finally, we find evidence that interbreeding practices were widespread during the domestication process by the early camelid herders in the Atacama during the Early Formative period and predating the Spanish conquest.

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.