Genetics of fibre production and fleece characteristics in small ruminants, Angora rabbit and South American camelids

Fiber quality and fertility in male alpacas in the Cusco region of Peru

Introduction

High testosterone levels might be associated with coarser fiber in alpacas, affecting fiber quality. In the husbandry systems employed in the Andes, males with higher libido might have an increased frequency of successful matings than males with lower libido. This study aimed to determine some of the factors affecting fiber quality in Peruvian alpacas and to evaluate the association between fiber quality and male mating behavior.

Methods

The study population consisted of 189 adult male alpacas at La Raya, Cusco, Peru, at 4,400 m above sea level, belonging to the National University of San Marcos, Lima, and the National University of San Antonio Abad of Cusco. Fiber samples were collected from male alpacas in September 2015; body condition score (BCS) was evaluated. After washing and drying, the fiber samples were analyzed using an Optical Fiber Diameter Analyser 2000; linear regression analysis between predictor variables (location of alpacas, age, BCS, and fiber color) and fiber quality outcomes was conducted. Fertility data were available only for some males belonging to San Marcos University, consisting of hand-written records from the breeding season January to April 2015; individual fertility quotients were calculated for each male.

Results

Age was associated with fiber quality, young alpacas having the finest fibers (p < 0.05). An increased BCS was associated with increased fiber length, fiber diameter and spinning fineness, but decreased comfort factor (p < 0.05). White fibers were thinner (p = 0.05) than colored fibers, with lower comfort factor and spinning fineness. A significant association between the fertility quotient and fiber curvature was observed (p = 0.018).

Discussion

These results suggest that careful selection of breeding individuals and attention to husbandry could result in improved fiber quality among alpaca herds in Peru. However, it would be advisable to increase the number of males studied, using more reliable methods for evaluating male fertility and pregnancy diagnosis than were available for this study.

Transcriptomic Analysis Reveals Candidate Ligand-Receptor Pairs and Signaling Networks Mediating Intercellular Communication between Hair Matrix Cells and Dermal Papilla Cells from Cashmere Goats

Hair fiber growth is determined by the spatiotemporally controlled proliferation, differentiation, and apoptosis of hair matrix cells (HMCs) inside the hair follicle (HF); however, dermal papilla cells (DPCs), the cell population surrounded by HMCs, manipulate the above processes via intercellular crosstalk with HMCs. Therefore, exploring how the mutual commutations between the cells are molecularly achieved is vital to understanding the mechanisms underlying hair growth. Here, based on our previous successes in cultivating HMCs and DPCs from cashmere goats, we combined a series of techniques, including in vitro cell coculture, transcriptome sequencing, and bioinformatic analysis, to uncover ligand-receptor pairs and signaling networks mediating intercellular crosstalk. Firstly, we found that direct cellular interaction significantly alters cell cycle distribution patterns and changes the gene expression profiles of both cells at the global level. Next, we constructed the networks of ligand-receptor pairs mediating intercellular autocrine or paracrine crosstalk between the cells. A few pairs, such as LEP-LEPR, IL6-EGFR, RSPO1-LRP6, and ADM-CALCRL, are found to have known or potential roles in hair growth by acting as bridges linking cells. Further, we inferred the signaling axis connecting the cells from transcriptomic data with the advantage of CCCExplorer. Certain pathways, including INHBA-ACVR2A/ACVR2B-ACVR1/ACVR1B-SMAD3, were predicted as the axis mediating the promotive effect of INHBA on hair growth via paracrine crosstalk between DPCs and HMCs. Finally, we verified that LEP-LEPR and IL1A-IL1R1 are pivotal ligand-receptor pairs involved in autocrine and paracrine communication of DPCs and HMCs to DPCs, respectively. Our study provides a comprehensive landscape of intercellular crosstalk between key cell types inside HF at the molecular level, which is helpful for an in-depth understanding of the mechanisms related to hair growth.

Assessment of hair and cashmere properties and their genetic background of several goat breeds in Southwest China

The aim of this study was to determine the properties and quality characteristics of hair and cashmere fibres of three goat breeds raised in Southwest China, namely; Dazu black goat (DBG, n = 203; ♂99, ♀104), Inner Mongolia cashmere goat (IMCG, n = 65; 21♂, 44♀) and their first cross (F₁, n = 79; 39♂, 40♀). Totals of 5219, 2130 and 2981 fibre samples, from the three breeds respectively, were taken prior to shearing at 32.5 ± 01.25 months of age from four body sites; shoulder, side-portion, abdomen and leg. Breed effect was significant (P < 0.01) for most hair and cashmere properties. IMCG and F₁ hair lengths were longer (P < 0.001) with less variable lengths than DBG. Shoulder hair diameters of the three breeds were not different (P > 0.05) but biggest of the side-portion and abdomen sites of DBG were bigger (P > 0.01), however, the biggest (P < 0.001) hair diameter was recorded for the leg site of F₁ and the smallest (P > 0.01) for IMCG. IMCG recorded the longest value for cashmere lengths followed by DBG, while F₁ recorded the lowest (P = 0.001), whilst F₁ recorded the biggest (P = 0.001) diameter whereas no differences existed between parents' breeds. The cortical cell lengths of IMCG and DBG were 94.57 and 86.85 μm without significant difference detected between breeds. Differences between hair length and diameter for body sites of the studied goat breeds were significant (P < 0.01) but between whiteness, cashmere diameter and diameter of cortical cells were not. Sex had no significant effect on all hair/cashmere properties. Quality characteristics of cashmere fibres from IMCG and F₁ were better (P < 0.001) than from DBG. Leg hair diameter, curl recovery rate and cashmere diameter were superior in the crossbred F₁ compared to pure breed parents, and DBG was superior to IMCG for fibre elasticity and intensity traits. FGF-5 gene was detected as a candidate gene for hair and cashmere traits in IMCG breed. Whilst, KIT gene was found to be associated with coat colour in the studied breeds. Extra investigations to examine more cashmere goat breeds and crosses are needed to discover genetic variability in cashmere production locally and worldwide.

Establishment of Angora rabbits' whisker hair follicle model and optimal culture conditions

To establish the model of whisker hair follicle culture in vitro and explore the best culture conditions, the whisker hair follicles of Angora rabbits were separated with stereomicroscope and cultured in William's E, DMEM, MEM media. The surface of the cultured whisker hair follicles was not damaged due to manual operation, resulting in the hair shaft's growth. This indicated the success of the in vitro whisker hair follicle model. The hair shaft grew at the fastest rate in the William's E culture (p < 0.05), which was significantly higher than that in the DMEM and MEM media. The hematoxylin-eosin results showed that compared to the William's E group, the atrophy of whisker hair follicles in the DMEM and MEM media was evident, especially in the MEM medium. PCNA immunofluorescence staining was employed to detect the expression of whisker hair follicles. The results showed that the PCNA positive expression of the William's E group was significantly stronger than that of the DMEM and MEM groups. Furthermore, CCK-8 and Annexin V-FITC/PI methods were used to detect the proliferation and apoptosis of the dermal papilla cells (DPCs). The results of this study provide a model for studying the hair growth of fur animals.

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.

Identification of potential key genes and pathways associated with the Pashmina fiber initiation using RNA-Seq and integrated bioinformatics analysis

Pashmina goat (Capra hircus) is an economically important livestock species, which habitats the cold arid desert of the Ladakh region (India), and produces a princely animal fiber called Pashmina. The Pashmina goat has a double coat fleece as an adaptation to the very harsh cold winters the outer long coarse hair (guard hair) produced from primary hair follicles and the inner fine Pashmina fiber produced from secondary hair follicles. Pashmina fiber undergoes a circannual and synchronized growth cycle. In the present study, we analyzed transcriptome profiles from 10 different Pashmina goats during anagen and telogen to delineate genes and signaling pathways regulating active (anagen) and regressive (telogen) phases of the follicle growth. During anagen, 150 genes were expressed at significantly higher levels with log (FC) > 2 and p_adj < 0.05. The RNA seq results were subjected to qRT-PCR validation. Among the nine genes selected, the expression of HAS1, TRIB2, P2RX1. PRG4, CNR2, and MMP25 were significantly higher (p < 0.05) in the anagen phase, whereas MC4R, GIPC2, and CDO1 were significantly expressed (p < 0.05) in the telogen phase which supports and validates the gene expression pattern from the RNA-sequencing. Differentially expressed genes revealed that Pashmina fiber initiation is largely controlled by signaling pathways like Wnt, NF-Kappa, JAK-STAT, Hippo, MAPK, Calcium, and PI3K-Akt. Expression of genes from the Integrin family, Cell adhesion molecules, and ECM-receptors were observed to be at much higher levels during anagen. We identified key genes (IL36RN, IGF2, ITGAV, ITGA5, ITCCR7, CXCL5, C3, CCL19, and CXCR3) and a collagen cluster which might be tightly correlated with anagen-induction. The regulatory network suggests the potential role of RUNX3, NR2F1/2, and GATA family transcription factors in anagen-initiation and maintaining fiber quality in Pashmina goats.

Rabbits - their domestication and molecular genetics of hair coat development and quality

The European rabbit (Oryctolagus cuniculus) is the only representative of its genus living in present-day Europe and North Africa, and all domestic rabbits are descendants of this one species, which is native to the Iberian Peninsula. There are over 300 breeds of rabbits that differ in size, coat color, length of ears and type of fur. Rabbits are bred for various reasons, such as for laboratory animals and a source of meat, wool and fur, as well as for pets and exhibition animals. The hair coat is a important economic trait of rabbits. Its development and quality are influenced by various factors, both environmental and genetic. The genetic mechanisms underlying its development have not been thoroughly researched. The aim of this review is to discuss the domestication of rabbits and the different aspects of rabbit genetics. A brief review of the properties of rabbit hair coat, hair coat development and hair cycle will be provided, followed by discussion of the factors regulating hair coat development, molecular control of hair coat development and the role of non-coding RNAs in the regulation of gene expression in the hair follicles of rabbits. Information about genetic regulation of pathways could provide useful tools for improving hair coat quality and be of practical use in rabbit breeding.

Shedding light on cashmere goat hair follicle biology: from morphology analyses to transcriptomic landascape

BACKGROUND

Cashmere goat is known for its precious undercoat. Being photoperiod-dictated, cashmere growth has been studied focusing mainly on hair follicle cycle phases (anagen, catagen and telogen). An accurate molecular knowledge of the goat hair follicle cycle, disentangling gene expression changes during phases and recognizing timing boundaries, could be useful to improve cashmere goat management and ultimately cashmere production.

RESULTS

To better describe goat's hair follicle transcriptome we applied RNA-sequencing to isolated hair follicles from five Italian cashmere goats, during the anagen and catagen phase, identifying total of 214 differentially expressed genes (DEGs): 97 were up-regulated while 117 were down-regulated in catagen with respect to anagen. Gene Ontology and pathway analysis were performed. We detected 144 significant pathways spanning from estrogen, pluripotency of stem cells, thermogenesis and fatty acid metabolism that were strongly expressed during the hair follicle phases analysed. Finally, we validated promising DEGs by RT-qPCR in the same set of samples as well as in hair follicles and entire skin biopsies of another cashmere goats cohort accounting for early anagen, anagen, early catagen, and catagen phases.

CONCLUSIONS

As in the isolated hair follicles, some target genes were homogenously modulated during the four hair follicle phases. Ceruloplasmin (CP) and Keratin 4 (K4), confirmed their clear cut expression between growing and resting phase. In fact, K4 was almost absent in catagen phases while CP was barely expressed in anagen phases. In particular, the strong expression of K4 in early anagen makes it an eligible marker to track the beginning of a new hair cycle, and therefore defining the optimum time for cashmere harvesting.

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.

Optimum wool harvest interval of angora rabbits under organised farm conditions in East China

The present study was conducted to evaluate the commercial lifespan and optimum wool harvest interval of Angora rabbits. One hundred shorn Angora rabbits were housed in an organised farm to describe the wool production curve. It showed that the optimum wool harvest interval was 75 d, when fibre length reached 55.0 mm. Wool production was lower from the 3^rd to the 6^th mo (young stage) than from 7^th to 28^th mo and rapidly decreased from 28^th to 31^st mo and was the lowest from 31^st to 33^rd mo of age. Feed intake-to-wool production ratio was higher from 3^rd to 4^th and from 7^th to 9^th mo of age than during the adult stage, and increased from 31^st to 33^rd mo of age. Daily weight gain was significantly higher from 3^rd to 4^th mo of age than in any other periods of the adult stage, and was negative from 23^rd to 33^rd mo of age. Therefore, the study reveals that the commercial lifespan of Angora rabbits was approximately 28 mo. Furthermore, wool production was higher in spring and winter than in autumn, and was the lowest in summer. Concomitantly, feed intake-to-wool production ratio was lower in spring than in autumn and winter, and was the highest in summer. Finally, daily weight gain was higher in spring and autumn than in winter, and was the lowest in summer. This indicates that wool production was depending on the season, and decreased significantly in summer. Moreover, the spring provided the best conditions for Angora rabbits.

RAPID COMMUNICATION: Generation of FGF5 knockout sheep via the CRISPR/Cas9 system

Sheep are an important source of fiber production. Fibroblast growth factor 5 (FGF5) is a dominant inhibitor of length of the anagen phase of the hair cycle. Knockout or silencing of the gene results in a wooly coat in mice, donkeys, dogs, and rabbits. In sheep breeding, wool length is one of the most important wool quality traits. However, traditional breeding cannot accurately and efficiently mediate an advanced genotype into the sheep genome. In this study, we generated 3 knockout sheep via the 1-step clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9 system. Sequencing analysis confirmed that mutations in the gene existed in all germ lines of 3 founders: besides the intact sequence, 3 kinds of deletions in the gene (including 5, 13, and 33 bp) were detected. The changes in the primary and senior structure of the FGF5 protein due to the 3 deletions in founders suggested that the FGF5 protein was dysfunctional. In addition, the expression level of intact mRNA in heterozygous individuals decreased compared with the wild types ( < 0.01). Functionally, we discovered that wool length in founders was significantly longer than in wild types ( < 0.05). Collectively, the knockout sheep with the longer wool length phenotype will provide an efficient way for fast genetic improvement of sheep breeding and promote the development of wool industry.

Competition for Materno-Fetal Resource Partitioning in a Rabbit Model of Undernourished Pregnancy

The major goal of animal production is to obtain abundant and healthy meat for consumers. Maternal food restriction (MFR) is often applied in farms to reduce production costs. However, the suitability of MFR in livestock animals is questionable, as this management may compromise maternal fitness due to a severe negative energetic balance and can induce Intrauterine Growth Restriction (IUGR) and prenatal programming in the offspring. Here, we sought to determine, using pregnant rabbits, the consequences of MFR on maternal endocrine and metabolic status and conceptus development. Pregnant dams were distributed into three groups: CONTROL (ad libitum feeding throughout the entire pregnancy; mean pregnancy length being around 31 days), UNDERFED (50% MFR during the entire pregnancy) and EARLY-UNDERFED (50% MFR only during the preimplantation period, Days 0–7). Maternal leptin concentrations and glycemic and lipid profiles were determined throughout pregnancy, whilst conceptus development was assessed ex-vivo at Day 28. Placental parameters were determined by macroscopic and histological evaluations and apoptotic assessments (TUNEL and Caspase-3). The main results of the study showed that, despite MFR altered maternal plasma lipid concentration (P<0.05), there were no effects on maternal bodyweight, plasma leptin concentration or glycemic profile. Fetal crown-rump lengths were reduced in both undernourished groups (P<0.001), but a significant reduction in fetal weight was only observed in the UNDERFED group (P<0.001). Growth in both undernourished groups was asymmetrical, with reduced liver weight (P<0.001) and significantly increased brain: fetal weight-ratio (P<0.001) and brain: liver weight-ratio (P<0.001) when compared to the CONTROL group. A significant reduction in placental weight was only observed in the UNDERFED group (P<0.001), despite both undernourished groups showing higher apoptotic rates at decidua and labyrinth zone (P<0.05) than the CONTROL group. Thus, these groups evidenced signs of placental degeneration, necrosis and stromal collapse. In summary, MFR may encourage the mother to make strategic decisions to safeguard her metabolic status and fitness at the expense of growth reduction in the litter, resulting in enhanced apoptotic and pathological processes at placental level and IUGR.

Rare phenotypes in domestic animals: unique resources for multiple applications

Preservation of specific and inheritable phenotypes of current or potential future importance is one of the main purposes of conservation of animal genetic resources. In this review, we investigate the issues behind the characterisation, utilisation and conservation of rare phenotypes, considering their multiple paths of relevance, variable levels of complexity and mode of inheritance. Accurately assessing the rarity of a given phenotype, especially a complex one, is not a simple task, because it requires the phenotypic and genetic characterisation of a large number of animals and populations and remains dependent of the scale of the study. Once characterised, specific phenotypes may contribute to various purposes (adaptedness, production, biological model, aesthetics, etc.) with adequate introgression programmes, which justifies the consideration of (real or potential) existence of such characteristics in in situ or ex situ conservation strategies. Recent biotechnological developments (genomic and genetic engineering) will undoubtedly bring important changes to the way phenotypes are characterised, introgressed and managed.

Molecular Characterization of the Llamas (Lama glama) Casein Cluster Genes Transcripts (CSN1S1, CSN2, CSN1S2, CSN3) and Regulatory Regions

In the present paper, we report for the first time the characterization of llama (Lama glama) caseins at transcriptomic and genetic level. A total of 288 casein clones transcripts were analysed from two lactating llamas. The most represented mRNA populations were those correctly assembled (85.07%) and they encoded for mature proteins of 215, 217, 187 and 162 amino acids respectively for the CSN1S1, CSN2, CSN1S2 and CSN3 genes. The exonic subdivision evidenced a structure made of 21, 9, 17 and 6 exons for the αs1-, β-, αs2- and κ-casein genes respectively. Exon skipping and duplication events were evidenced. Two variants A and B were identified in the αs1-casein gene as result of the alternative out-splicing of the exon 18. An additional exon coding for a novel esapeptide was found to be cryptic in the κ-casein gene, whereas one extra exon was found in the αs2-casein gene by the comparison with the Camelus dromedaries sequence. A total of 28 putative phosphorylated motifs highlighted a complex heterogeneity and a potential variable degree of post-translational modifications. Ninety-six polymorphic sites were found through the comparison of the lama casein cDNAs with the homologous camel sequences, whereas the first description and characterization of the 5'- and 3'-regulatory regions allowed to identify the main putative consensus sequences involved in the casein genes expression, thus opening the way to new investigations -so far- never achieved in this species.

Weighting fibre and morphological traits in a genetic index for an alpaca breeding programme

Nowadays, the fibre diameter (FD) is considered the main selection objective in alpaca populations all over the world. International Committee for Animal Recording recommendations define the FD and its CV as the first two traits to be considered in breeding programmes for this specie. In addition to these main criteria, other selection criteria of economic value used are comfort factor (CF) or standard deviation (s.d.); also other less important traits being used as selection objectives are these morphological traits: density (DE), crimp (CR) or lock structure (LS) for, respectively, Huacaya (HU) and Suri (SU) ecotypes, head (HE), coverage (CO) and balance (BA). The goal of this study was to establish how to implement a combined selection index starting from genetic parameters and to study the expected correlation between genetic trends by considering different alternative procedures of weighting all the involved traits, and the consequences of a wrongly proceeding way. Heritabilities and genetic and phenotypic correlations were estimated from the data set belonging to the PACOMARCA experimental farm for SU and HU. Two approaches were used to check the consequences of a set of subjective weights essayed. The coefficients of selection indexes were obtained for two sets of reference weights. In addition, equivalent weights were drawn if applied those reference values as coefficients of hypothetical selection indexes directly on phenotypes; relative expected genetic responses were computed in different cases. Results showed that almost in all cases for both ecotypes, the weight applied to CF should be surprisingly negative. Concerning genetic responses, only CO was compromised in some cases for the HU ecotype. The essayed methodology allowed explaining the differences between ecotypes in the genetic trends. The proposed methodology was shown to be effective to study the relative importance of the traits granted by the manager of a breeding scheme.

De Novo Transcriptome Assembly and Differential Gene Expression Profiling of Three Capra hircus Skin Types during Anagen of the Hair Growth Cycle

Despite that goat is one of the best nonmodel systems for villus growth studies and hair biology, limited gene resources associated with skin or hair follicles are available. In the present study, using Illumina/Solexa sequencing technology, we de novo assembled 130 million mRNA-Seq reads into a total of 49,115 contigs. Searching public databases revealed that about 45% of the total contigs can be annotated as known proteins, indicating that some of the assembled contigs may have previously uncharacterized functions. Functional classification by KOG and GO showed that activities associated with metabolism are predominant in goat skin during anagen phase. Many signaling pathways was also created based on the mapping of assembled contigs to the KEGG pathway database, some of which have been previously demonstrated to have diverse roles in hair follicle and hair shaft formation. Furthermore, gene expression profiling of three skin types identified ~6,300 transcript-derived contigs that are differentially expressed. These genes mainly enriched in the functional cluster associated with cell cycle and cell division. The large contig catalogue as well as the genes which were differentially expressed in different skin types provide valuable candidates for further characterization of gene functions.

Differential gene expression analysis between anagen and telogen of Capra hircus skin based on the de novo assembled transcriptome sequence

Capra hircus, an economically important livestock, plays an indispensable role in the world animal fiber industry. In the present study, using Illumina/Solexa high throughput sequencing technology, we sequenced and de novo assembled the goat skin transcriptome corresponding to the anagen and telogen of the hair growth cycle. Approximately 53Mb of transcriptome sequences consisting of 57,040 high quality contigs was obtained. More than 8300 contigs were predicted to contain a full length coding sequence. Approximately 43% of the total contigs were identified as harboring homologs of sequences from other organisms in the public database. Based on the assembled transcript-derived contigs, we identified about 7000 transcripts that were differentially expressed between the anagen and telogen libraries. These differentially expressed genes were mainly enriched in signal transduction mechanisms, extracellular structures and cytoskeleton from the KOG database and in ECM receptor interaction, focal adhesion and gap junction from the KEGG pathway database, indicating the essential roles of these genes may play in cell-to-cell and cell-to-matrix communications during the active hair growth phase. In addition, many signaling pathway associated ligands and/or receptors were also identified as up-regulated genes during the anagen phase compared with the telogen stage, suggesting that enhanced cross-talk among signaling transduction pathways may be required for anagen of the hair cycle. These differentially expressed genes, especially those that were over-represented in each of the functional clusters and biochemical pathways, provide valuable resources and opportunities for characterizing the gene functions associated with hair fiber growth as well as for breeding elite Cashmere goat species.

Fleece variation in alpaca (Vicugna pacos): a two-locus model for the Suri/Huacaya phenotype

Background

Genetic improvement of fibre-producing animal species has often induced transition from double coated to single coated fleece, accompanied by dramatic changes in skin follicles and hair composition, likely implying variation at multiple loci. Huacaya, the more common fleece phenotype in alpaca (Vicugna pacos), is characterized by a thick dense coat growing perpendicularly from the body, whereas the alternative rare and more prized single-coated Suri phenotype is distinguished by long silky fibre that grows parallel to the body and hangs in separate, distinctive pencil locks. A single-locus genetic model has been proposed for the Suri-Huacaya phenotype, where Huacaya is recessive.

Results

Two reciprocal experimental test-crosses (Suri × Huacaya) were carried out, involving a total of 17 unrelated males and 149 unrelated females. An additional dataset of 587 offspring of Suri × Suri crosses was analyzed. Segregation ratios, population genotype frequencies, and/or recombination fraction under different genetic models were estimated by maximum likelihood. The single locus model for the Suri/Huacaya phenotype was rejected. In addition, we present two unexpected observations: 1) a large proportion (about 3/4) of the Suri animals are segregating (with at least one Huacaya offspring), even in breeding conditions where the Huacaya trait would have been almost eliminated; 2) a model with two different values of the segregation ratio fit the data significantly better than a model with a single parameter.

Conclusions

The data support a genetic model in which two linked loci must simultaneously be homozygous for recessive alleles in order to produce the Huacaya phenotype. The estimated recombination rate between these loci was 0.099 (95% C.L. = 0.029-0.204). Our genetic analysis may be useful for other species whose breeding system produces mainly half-sib families.