Goat Genome Research Articles

Comparative genomic approach reveals novel conserved microRNAs in Inner Mongolia cashmere goat skin and longissimus dorsi.

MicroRNA (miRNA) is a large class of non-coding RNA which usually acts a fine-tuned regulation in repressing gene expression on transcriptional and post-transcriptional level during ontogenetic development, metabolism, and occurrence of disease. Due to the lateness of goat genome investigation, registered goat microRNAs are little known and the function of it is poorly understood. In this study, we identified 5 novel miRNAs express in cashmere goat skin and longissimus dorsi using comparative genomic approach combined with expression profile analysis. Further qPCR and sequencing validation demonstrated that the novel miRNAs we identified expressed in goat skin and longissimus dorsi, three of them with the highest expression in February and two in October, indicating these novel microRNAs were involved in growth and cessation of goat hair production. Subsequently, potential target genes were explored via base-pairing with complementary sequences miRNA/mRNA interaction searching. WDR12 and CSNK1A1 involved in Notch/Wnt signaltransduction pathway were finally identified. Collectively, the approach of comparative genomics combining expression profile analysis is a powerful tool to identify novel miRNA in goat and the 5 identified novel miRNAs(chi-miR-2284n, chi-miR-421*, chi-miR-421, chi-miR-1839 and chi-miR-374; Accession number: JQ002550-JQ002554) play a essential role in goat hair production of cashmere goat, both in entering growth and cessation phase. This study is also a meaningful complement to the present goat miRNA database for further understanding function of miRNAs in regulation of goat cashmere production.

An update of the goat genome assembly using dense radiation hybrid maps allows detailed analysis of evolutionary rearrangements in Bovidae.

BackgroundThe domestic goat (Capra hircus), an important livestock species, belongs to a clade of Ruminantia, Bovidae, together with cattle, buffalo and sheep. The history of genome evolution and chromosomal rearrangements on a small scale in ruminants remain speculative. Recently completed goat genome sequence was released but is still in a draft stage. The draft sequence used a variety of assembly packages, as well as a radiation hybrid (RH) map of chromosome 1 as part of its validation.ResultsUsing an improved RH mapping pipeline, whole-genome dense maps of 45,953 SNP markers were constructed with statistical confidence measures and the saturated maps provided a fine map resolution of approximate 65 kb. Linking RH maps to the goat sequences showed that the assemblies of scaffolds/super-scaffolds were globally accurate. However, we observed certain flaws linked to the process of anchoring chromosome using conserved synteny with cattle. Chromosome assignments, long-range order, and orientation of the scaffolds were reassessed in an updated genome sequence version. We also present new results exploiting the updated goat genome sequence to understand genomic rearrangements and chromosome evolution between mammals during species radiations. The sequence architecture of rearrangement sites between the goat and cattle genomes presented abundant segmental duplication on regions of goat chromosome 9 and 14, as well as new insertions in homologous cattle genome regions. This complex interplay between duplicated sequences and Robertsonian translocations highlights the rearrangement mechanism of centromeric nonallelic homologous recombination (NAHR) in mammals. We observed that species-specific shifts in ANKRD26 gene duplication are coincident with breakpoint reuse in divergent lineages and this gene family may play a role in chromosome stabilization in chromosome evolution.ConclusionsWe generated dense maps of the complete whole goat genome. The chromosomal maps allowed us to anchor and orientate assembled genome scaffolds along the chromosomes, annotate chromosome rearrangements and thereby get a better understanding of the genome evolution of ruminants and other mammals.Electronic supplementary materialThe online version of this article (doi:10.1186/1471-2164-15-625) contains supplementary material, which is available to authorized users.

Design and Characterization of a 52K SNP Chip for Goats

The success of Genome Wide Association Studies in the discovery of sequence variation linked to complex traits in humans has increased interest in high throughput SNP genotyping assays in livestock species. Primary goals are QTL detection and genomic selection. The purpose here was design of a 50–60,000 SNP chip for goats. The success of a moderate density SNP assay depends on reliable bioinformatic SNP detection procedures, the technological success rate of the SNP design, even spacing of SNPs on the genome and selection of Minor Allele Frequencies (MAF) suitable to use in diverse breeds. Through the federation of three SNP discovery projects consolidated as the International Goat Genome Consortium, we have identified approximately twelve million high quality SNP variants in the goat genome stored in a database together with their biological and technical characteristics. These SNPs were identified within and between six breeds (meat, milk and mixed): Alpine, Boer, Creole, Katjang, Saanen and Savanna, comprising a total of 97 animals. Whole genome and Reduced Representation Library sequences were aligned on >10 kb scaffolds of the de novo goat genome assembly. The 60,000 selected SNPs, evenly spaced on the goat genome, were submitted for oligo manufacturing (Illumina, Inc) and published in dbSNP along with flanking sequences and map position on goat assemblies (i.e. scaffolds and pseudo-chromosomes), sheep genome V2 and cattle UMD3.1 assembly. Ten breeds were then used to validate the SNP content and 52,295 loci could be successfully genotyped and used to generate a final cluster file. The combined strategy of using mainly whole genome Next Generation Sequencing and mapping on a contig genome assembly, complemented with Illumina design tools proved to be efficient in producing this GoatSNP50 chip. Advances in use of molecular markers are expected to accelerate goat genomic studies in coming years.

Seasonal variation of urinary microRNA expression in male goats (Capra hircus) as assessed by next generation sequencing

Testosterone plays a key role in preparation of a male domesticated goat (Capra hircus) to breeding season including changes in the urogenital tract of a male goat (buck). microRNAs are important regulators of cellular metabolism, differentiation and function. They are powerful intermediaries of hormonal activity in the body, including the urogenital tract. We investigated seasonal changes in expression of microRNAs in goat buck urine and their potential consequences using next generation sequencing (microRNA-Seq). We determined the location of each microRNA gene in the goat genome. Testosterone was measured by radioimmunoassay and the androgen receptor binding sites (ARBS) in the promoters of the microRNA genes were determined by MatInspector. The overall impact of regulated microRNAs on cellular physiology was assessed by mirPath. We observed high testosterone levels during the breeding season and changes in the expression of forty microRNAs. Nineteen microRNAs were upregulated, while twenty-one were downregulated. We identified several ARBS in the promoters of regulated microRNAs. Notably, the mostly inhibited microRNA, miR-1246, has a unique set of several gene copy variants associated with a cluster of androgen receptor binding sites. Concomitant changes in regulated microRNA expression could promote transcription, proliferation and differentiation of urogenital tract cells. Together, these findings indicate that in a domesticated goat (Capra hircus), there are specific changes in the microRNA expression profile in buck urine during breeding season, which could be attributable to high testosterone levels during breeding, and could help in preparation of the urogenital tract for high metabolic demands of that season.

Cloning and transcriptional activity of ovarian-specific promoter in dairy goats

Ovarian-specific promoter (OSP) is a tissue-specific promoter only expressed in the ovary to regulate its development. To investigate the activity of OSP in diverse dairy goat cells, OSP was cloned from dairy goat genome and used to construct the luciferase reporter vector pGL3-OSP to examine the transcriptional activity of OSP. Consequently, a 461-bp OSP fragment was obtained. Bioinformatics analysis indicated that the fragment contained multiple transcription factor binding sites, such as SRY, HSF, AP-1, C/EBP, CAAT, TATA, and so on, and the nucleotide sequence of OSP shared 99% and 97% similarity with those of OSP-1 and OSP-2 in rat. In addition, two variants were identified in the OSP (g.164T was deletion and g.375T<C). In diverse cells, the expression of luciferase in ovarian granulosa cells (GCs) and stromal cells of ovarian cortex was extremely significantly higher than that in fetal fibroblasts and oviduct epithelial cells. Meanwhile, the expression of luciferase in ovarian GCs was extremely significantly higher than that in the stromal cells of ovarian cortex. In summary, we successfully extracted dairy goat OSP and identified its transcriptional activity, which in part shed light on the study of ovarian-specific transcription and provided further insights into the development and regulation of dairy goat's ovary as well.

A tailing genome walking method suitable for genomes with high local GC content

The tailing genome walking strategies are simple and efficient. However, they sometimes can be restricted due to the low stringency of homo-oligomeric primers. Here we modified their conventional tailing step by adding polythymidine and polyguanine to the target single-stranded DNA (ssDNA). The tailed ssDNA was then amplified exponentially with a specific primer in the known region and a primer comprising 5′ polycytosine and 3′ polyadenosine. The successful application of this novel method for identifying integration sites mediated by φC31 integrase in goat genome indicates that the method is more suitable for genomes with high complexity and local GC content.

Targeted editing of goat genome with modular-assembly zinc finger nucleases based on activity prediction by computational molecular modeling

Zinc finger nuclease (ZFN) technology can mediate targeted genome modification to produce transgenic animals in a high-efficient and biological-safe way. Modular assembly is a rapid, convenient and open-source method for the synthesis of ZFNs. However, this biotechnology is hampered by multistep construction, low-efficiency editing and off-target cleavage. Here we synthesized and tested six pairs of three- or four-finger ZFNs to target one site in goat beta-lactoglobulin (BLG, a dominant allergen in goat milk) gene. Homology modeling was applied to build the structure model of ZFNs to predict their editing activities targeting at goat BLG gene. Goat fibroblast cells were transfected with plasmids that encoded ZFN pairs, and genomic DNA was isolated 72 h later for genome editing efficiency assay. The results of editing efficiency assay demonstrated that ZFNs with optimal interaction modes can edit goat BLG gene more efficiently, whereas ZFNs with unexpected interaction modes showed lower activities in editing BLG gene. We concluded that modular-assembly ZFNs can provide a rapid, public-available, and easy-to-practice platform for transgenic animal research and molecular modeling would help as a useful tool for ZFNs activity prediction.

Goat genome sequence by optical mapping

Goat genome sequence by optical mapping

Genetic diversity and investigation of polledness in divergent goat populations using 52 088 SNPs

The recent availability of a genome-wide SNP array for the goat genome dramatically increases the power to investigate aspects of genetic diversity and to conduct genome-wide association studies in this important domestic species. We collected and analysed genotypes from 52 088 SNPs in Boer, Cashmere and Rangeland goats that had both polled and horned individuals. Principal components analysis revealed a clear genetic division between animals for each population, and model-based clustering successfully detected evidence of admixture that matched aspects of their recorded history. For example, shared co-ancestry was detected, suggesting Boer goats have been introgressed into the Rangeland population. Further, allele frequency data successfully tracked the altered genetic profile that has taken place after 40 years of breeding Australian Cashmere goats using the Rangeland animals as the founding population. Genome-wide association mapping of the POLL locus revealed a strong signal on goat chromosome 1. The 769-kb critical interval contained the polled intersex syndrome locus, confirming the genetic basis in non-European animals is the same as identified previously in Saanen goats. Interestingly, analysis of the haplotypes carried by a small set of sex-reversed animals, known to be associated with polledness, revealed some animals carried the wild-type chromosome associated with the presence of horns. This suggests a more complex basis for the relationship between polledness and the intersex condition than initially thought while validating the application of the goat SNP50 BeadChip for fine-mapping traits in goat.

A note on the organization and expression of β-defensin genes in Polish goats

Nucleotide (cds) and amino acids sequences of the caprine β2-defensin genes were in silico compared to search for the sequence variation and for the LAP gene sequences in the goat genome and for the presence of LAP gene transcripts in goat tissues. The comparison of the exon sequences revealed that the first 64 amino acids are identical in both LAP and β1-defensin. However, the GBD-1 prepropeptide is shorter by 18 amino acids due to the presence of the stop codon UAA at position 209–211 in GBD-1 mRNA. The LAP gene, which was found, so far, only in Indian goat breeds, is absent in the genome of Polish dairy goats. The introns of the caprine β1- and β2-defensin genes were, for the first time, sequenced; their sequences showed 99.6 % identity, differing in six nucleotide positions.

PiggyBac Transposon-Mediated Gene Transfer in Cashmere Goat Fetal Fibroblast Cells

PiggyBac (PB) has recently been found to be functional in various organisms. To verify and exploit its application in the cashmere goat, a PB transposon system including donor and helper vector of was developed, in which the EGFP gene in donor of vector was used as reporter. Cashmere goat fetal fibroblasts cells (GFFs) were transfected with the PB transposon system and the efficiency of gene transfer was determined. Compared with random integration, PB-mediated EGFP expression levels increased 7.78-fold in the GFFs, confirming that the PB transposon system constructed successfully mediated efficient foreign gene integration in the GFFs. To further investigate the characteristics of PB-mediated integration instance, PB integration site distribution in the goat genome was examined. The results showed that PB had a preference for AT rich regions of the goat genome. Thus this study confirms the function of PB transposon in GFFs and provides a potential genetic tool for producing transgenic goats.

A whole-genome radiation hybrid panel for goat

A 5000rad goat–hamster panel of 121 whole-genome radiation hybrids was generated and preliminarily characterized. A normal diploid fibroblast culture from a male Boer goat was fused with a recipient thymidine kinase-deficient hamster cell line. The generated 121 radiation hybrids were grown and produced an average of 8.4mg of DNA per radiation hybrid. A SINE-PCR test showed that almost all radiation hybrids retained goat DNA. The retention frequencies of the 121 hybrids were preliminarily estimated using a collection of 42 unlinked molecular markers. An optimized panel of therein 90 radiation hybrids (CHRH5000) with an average retention frequency of 34.2% was screened on the basis of interspersed repetitive DNA content and chromosome retention frequency. The development of this WG-RH panel will provide a fundamental tool for advanced goat genome mapping studies and for mammalian comparative mapping.

Transgenesis applied to goat: current applications and ongoing research

Compared to experiments involving pigs, cows and/or sheep, transgenesis applied to goats is probably less advertised. However, recent successes and increasing amount of dedicated research make this species of special interest for ongoing biological and physiological questions on genome engineering in large animals. This short review aims at highlighting the current applications and limitations of the goat genome manipulation.

Optimization of Saanen sperm genes amplification: evaluation of standardized protocols in genetically uncharacterized rural goats reared under a subtropical environment

The purpose of this research is to optimize quantitatively the amplification of specific sperm genes in reference genomically characterized Saanen goat and to evaluate the standardized protocols applicability on sperms of uncharacterized genome of rural goats reared under subtropical environment for inclusion in future selection programs. The optimization of the protocols in Saanen sperms included three production genes (growth hormone (GH) exons 2, 3, and 4, αS1-casein (CSN1S1), and α-lactalbumin) and two health genes (MHC class II DRB and prion (PrP)). The optimization was based on varying the primers concentrations and the inclusion of a PCR cosolvent (Triton X). The impact of the studied variables on statistically significant increase in the yield of amplicons was noticed in four out of five (80%) optimized protocols, namely in those related to GH, CSN1S1, α-lactalbumin, and PrP genes (P < 0.05). There was no significant difference in the yield of amplicons related to MHC class II DRB gene, regardless of the variables used (P > 0.05). The applicability of the optimized protocols of Saanen sperm genes on amplification of uncharacterized rural goat sperms revealed a 100% success in tested individuals for amplification of GH, CSN1S1, α-lactalbumin, and MHC class II DRB genes and a 75% success for the PrP gene. The significant success in applicability of the Saanen quantitatively optimized protocols to other uncharacterized genome of rural goats allows for their inclusion in future selection, targeting the sustainability of this farming system in a subtropical environment and the improvement of the farmers livelihood.

Stable silencing of β-lactoglobulin (BLG) gene by lentivirus-mediated RNAi in goat fetal fibroblasts

β-lactoglobulin (BLG), a dominant allergen in goat milk, is difficult to remove by traditional biochemical methods. Its elimination from goat milk by genetic modification therefore poses a major challenge for modern goat breeders. A shRNA targeting BLG mRNA with high interference efficiency was identified, with which lentiviral vectors were used for mediating stable shRNA interference in goat-fetal fibroblast cells. Apart from high efficiency in the knockdown of BLG expression in these cells, lentivector-mediated RNAi manifested stable integration into the goat genome itself. Consequently, an in vitro model for goat BLG-content control was compiled, and a goat-cell line for accompanying transgenetic goat production created.

RsaI repetitive DNA in Buffalo Bubalus bubalis representing retrotransposons, conserved in bovids, are part of the functional genes

BackgroundRepetitive sequences are the major components of the eukaryotic genomes. Association of these repeats with transcribing sequences and their regulation in buffalo Bubalus bubalis has remained largely unresolved.ResultsWe cloned and sequenced RsaI repeat fragments pDp1, pDp2, pDp3, pDp4 of 1331, 651, 603 and 339 base pairs, respectively from the buffalo, Bubalus bubalis. Upon characterization, these fragments were found to represent retrotransposons and part of some functional genes. The resultant clones showed cross hybridization only with buffalo, cattle, goat and sheep genomic DNA. Real Time PCR, detected ~2 × 104 copies of pDp1, ~ 3000 copies of pDp2 and pDp3 and ~ 1000 of pDp4 in buffalo, cattle, goat and sheep genomes, respectively. RsaI repeats are transcriptionally active in somatic tissues and spermatozoa. Accordingly, pDp1 showed maximum expression in lung, pDp2 and pDp3 both in Kidney, and pDp4 in ovary. Fluorescence in situ hybridization showed repeats to be distributed all across the chromosomes.ConclusionsThe data suggest that RsaI repeats have been incorporated into the exonic regions of various transcribing genes, possibly contributing towards the architecture and evolution of the buffalo and related genomes. Prospects of our present work in the context of comparative and functional genomics are highlighted.

A first comparative map of copy number variations in the sheep genome

We carried out a cross species cattle–sheep array comparative genome hybridization experiment to identify copy number variations (CNVs) in the sheep genome analysing ewes of Italian dairy or dual-purpose breeds (Bagnolese, Comisana, Laticauda, Massese, Sarda, and Valle del Belice) using a tiling oligonucleotide array with ~ 385,000 probes designed on the bovine genome. We identified 135 CNV regions (CNVRs; 24 reported in more than one animal) covering ~ 10.5 Mb of the virtual sheep genome referred to the bovine genome (0.398%) with a mean and a median equal to 77.6 and 55.9 kb, respectively. A comparative analysis between the identified sheep CNVRs and those reported in cattle and goat genomes indicated that overlaps between sheep and both other species CNVRs are highly significant (P < 0.0001), suggesting that several chromosome regions might contain recurrent interspecies CNVRs. Many sheep CNVRs include genes with important biological functions. Further studies are needed to evaluate their functional relevance.

An initial comparative map of copy number variations in the goat (Capra hircus) genome

BackgroundThe goat (Capra hircus) represents one of the most important farm animal species. It is reared in all continents with an estimated world population of about 800 million of animals. Despite its importance, studies on the goat genome are still in their infancy compared to those in other farm animal species. Comparative mapping between cattle and goat showed only a few rearrangements in agreement with the similarity of chromosome banding. We carried out a cross species cattle-goat array comparative genome hybridization (aCGH) experiment in order to identify copy number variations (CNVs) in the goat genome analysing animals of different breeds (Saanen, Camosciata delle Alpi, Girgentana, and Murciano-Granadina) using a tiling oligonucleotide array with ~385,000 probes designed on the bovine genome.ResultsWe identified a total of 161 CNVs (an average of 17.9 CNVs per goat), with the largest number in the Saanen breed and the lowest in the Camosciata delle Alpi goat. By aggregating overlapping CNVs identified in different animals we determined CNV regions (CNVRs): on the whole, we identified 127 CNVRs covering about 11.47 Mb of the virtual goat genome referred to the bovine genome (0.435% of the latter genome). These 127 CNVRs included 86 loss and 41 gain and ranged from about 24 kb to about 1.07 Mb with a mean and median equal to 90,292 bp and 49,530 bp, respectively. To evaluate whether the identified goat CNVRs overlap with those reported in the cattle genome, we compared our results with those obtained in four independent cattle experiments. Overlapping between goat and cattle CNVRs was highly significant (P < 0.0001) suggesting that several chromosome regions might contain recurrent interspecies CNVRs. Genes with environmental functions were over-represented in goat CNVRs as reported in other mammals.ConclusionsWe describe a first map of goat CNVRs. This provides information on a comparative basis with the cattle genome by identifying putative recurrent interspecies CNVs between these two ruminant species. Several goat CNVs affect genes with important biological functions. Further studies are needed to evaluate the functional relevance of these CNVs and their effects on behavior, production, and disease resistance traits in goats.

Capturing the Perfect Reference Genome

Capturing the Perfect Reference Genome

Production of transgenic piglets using ICSI–sperm-mediated gene transfer in combination with recombinase RecA

Sperm-mediated gene transfer (SMGT) is a method for the production of transgenic animals based on the intrinsic ability of sperm cells to bind and internalize exogenous DNA molecules and to transfer them into the oocyte at fertilization. Recombinase-A (RecA) protein-coated exogenous DNA has been used previously in pronuclear injection systems increasing integration into goat and pig genomes. However, there are no data regarding transgene expression after ICSI. Here, we set out to investigate whether the expression of transgenic DNA in porcine embryos is improved by recombinase-mediated DNA transfer and if it is possible to generate transgenic animals using this methodology. Different factors which could affect the performance of this transgenic methodology were analyzed by studying 1) the effect of the presence of exogenous DNA and RecA protein on boar sperm functionality; 2) the effect of recombinase RecA on in vitro enhanced green fluorescent protein (EGFP)-expressing embryos produced by ICSI or IVF; and 3) the efficiency of generation of transgenic piglets by RecA-mediated ICSI. Our results suggested that 1) the presence of exogenous DNA and RecA-DNA complexes at 5 microg/ml did not affect sperm functionality in terms of motility, viability, membrane lipid disorder, or reactive oxygen species generation; 2) EGFP-expressing embryos were obtained with a high efficiency using the SMGT-ICSI technique in combination with recombinase; however, the use of IVF system did not result in any fluorescent embryos; and 3) transgenic piglets were produced by this methodology. To our knowledge, this is the first time that transgenic pigs have been produced by ICSI-SGMT and a recombinase.