Lei Qu

Generation of gene-modified goats targeting MSTN and FGF5 via zygote injection of CRISPR/Cas9 system

Recent advances in the study of the CRISPR/Cas9 system have provided a precise and versatile approach for genome editing in various species. However, the applicability and efficiency of this method in large animal models, such as the goat, have not been extensively studied. Here, by co-injection of one-cell stage embryos with Cas9 mRNA and sgRNAs targeting two functional genes (MSTN and FGF5), we successfully produced gene-modified goats with either one or both genes disrupted. The targeting efficiency of MSTN and FGF5 in cultured primary fibroblasts was as high as 60%, while the efficiency of disrupting MSTN and FGF5 in 98 tested animals was 15% and 21% respectively, and 10% for double gene modifications. The on- and off-target mutations of the target genes in fibroblasts, as well as in somatic tissues and testis of founder and dead animals, were carefully analyzed. The results showed that simultaneous editing of several sites was achieved in large animals, demonstrating that the CRISPR/Cas9 system has the potential to become a robust and efficient gene engineering tool in farm animals, and therefore will be critically important and applicable for breeding.

Discovery of cashmere goat (Capra hircus) microRNAs in skin and hair follicles by Solexa sequencing

BackgroundMicroRNAs (miRNAs) are a large family of endogenous, non-coding RNAs, about 22 nucleotides long, which regulate gene expression through sequence-specific base pairing with target mRNAs. Extensive studies have shown that miRNA expression in the skin changes remarkably during distinct stages of the hair cycle in humans, mice, goats and sheep.ResultsIn this study, the skin tissues were harvested from the three stages of hair follicle cycling (anagen, catagen and telogen) in a fibre-producing goat breed. In total, 63,109,004 raw reads were obtained by Solexa sequencing and 61,125,752 clean reads remained for the small RNA digitalisation analysis. This resulted in the identification of 399 conserved miRNAs; among these, 326 miRNAs were expressed in all three follicular cycling stages, whereas 3, 12 and 11 miRNAs were specifically expressed in anagen, catagen, and telogen, respectively. We also identified 172 potential novel miRNAs by Mireap, 36 miRNAs were expressed in all three cycling stages, whereas 23, 29 and 44 miRNAs were specifically expressed in anagen, catagen, and telogen, respectively. The expression level of five arbitrarily selected miRNAs was analyzed by quantitative PCR, and the results indicated that the expression patterns were consistent with the Solexa sequencing results. Gene Ontology and KEGG pathway analyses indicated that five major biological pathways (Metabolic pathways, Pathways in cancer, MAPK signalling pathway, Endocytosis and Focal adhesion) accounted for 23.08% of target genes among 278 biological functions, indicating that these pathways are likely to play significant roles during hair cycling.ConclusionsDuring all hair cycle stages of cashmere goats, a large number of conserved and novel miRNAs were identified through a high-throughput sequencing approach. This study enriches the Capra hircus miRNA databases and provides a comprehensive miRNA transcriptome profile in the skin of goats during the hair follicle cycle.

Rumen bacterial diversity of 80 to 110-day-old goats using 16S rRNA sequencing.

The ability of rumen microorganisms to use fibrous plant matter plays an important role in ruminant animals; however, little information about rumen colonization by microbial populations after weaning has been reported. In this study, high-throughput sequencing was used to investigate the establishment of this microbial population in 80 to 110-day-old goats. Illumina sequencing of goat rumen samples yielded 101,356,610 nucleotides that were assembled into 256,868 reads with an average read length of 394 nucleotides. Taxonomic analysis of metagenomic reads indicated that the predominant phyla were distinct at different growth stages. The phyla Firmicutes and Synergistetes were predominant in samples taken from 80 to 100-day-old goats, but Bacteroidetes and Firmicutes became the most abundant phyla in samples from 110-day-old animals. There was a remarkable variation in the microbial populations with age; Firmicutes and Synergistetes decreased after weaning, but Bacteroidetes and Proteobacteria increased from 80 to 110 day of age. These findings suggested that colonization of the rumen by microorganisms is related to their function in the rumen digestive system. These results give a better understanding of the role of rumen microbes and the establishment of the microbial population, which help to maintain the host’s health and improve animal performance.

Insertion/Deletion Within the KDM6A Gene Is Significantly Associated With Litter Size in Goat

A previous whole-genome association analysis identified lysine demethylase 6A (KDM6A), which encodes a type of histone demethylase, as a candidate gene associated to goat fecundity. KDM6A gene knockout mouse disrupts gametophyte development, suggesting that it has a critical role in reproduction. In this study, goat KDM6A mRNA expression profiles were determined, insertion/deletion (indel) variants in the gene identified, indel variants effect on KDM6A gene expression assessed, and their association with first-born litter size analyzed in 2326 healthy female Shaanbei white cashmere goats. KDM6A mRNA was expressed in all tissues tested (heart, liver, spleen, lung, kidney, muscle, brain, skin and testis); the expression levels in testes at different developmental stages [1-week-old (wk), 2, 3 wk, 1-month-old (mo), 1.5 and 2 mo] indicated a potential association with the mitosis-to-meiosis transition, implying that KDM6A may have an essential role in goat fertility. Meanwhile, two novel intronic indels of 16 bp and 5 bp were identified. Statistical analysis revealed that only the 16 bp indel was associated with first-born litter size (P < 0.01), and the average first-born litter size of individuals with an insertion/insertion genotype higher than that of those with the deletion/deletion genotype (P < 0.05). There was also a significant difference in genotype distributions of the 16 bp indel between mothers of single-lamb and multi-lamb litters in the studied goat population (P = 0.001). Consistently, the 16 bp indel also had a significant effect on KDM6A gene expression. Additionally, there was no significant linkage disequilibrium (LD) between these two indel loci, consistent with the association analysis results. Together, these findings suggest that the 16 bp indel in KDM6A may be useful for marker-assisted selection (MAS) of goats.

Use of green fluorescent protein to monitor Lactobacillus plantarum in the gastrointestinal tract of goats.

The experiment aimed to specifically monitor the passage of lactobacilli in vivo after oral administration. The green fluorescent protein (GFP) gene was cloned downstream from the constitutive p32 promoter from L. lactis subsp. cremoris Wg2. The recombinant expression vector, pLEM415-gfp-p32, was electroporated into Lactobacillus plantarum (L. plantarum) isolated from goat. Green fluorescent protein (GFP) was successfully expressed in L. plantarum. After 2 h post-administration, transformed Lactobacillus could be detectable in all luminal contents. In the rumen, bacteria concentration initially decreased, reached the minimum at 42 h post-oral administration and then increased. However, this concentration decreased constantly in the duodenum. This result indicated that L. plantarum could colonize in the rumen but not in the duodenum.

CRISPR/Cas9-mediated MSTN disruption and heritable mutagenesis in goats causes increased body mass

Genetic engineering in livestock has been greatly enhanced through the use of artificial programmed nucleases such as the recently emerged clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated 9 (Cas9) system. We recently reported our successful application of the CRISPR/Cas9 system to engineer the goat genome through micro-injection of Cas9 mRNA and sgRNAs targeting MSTN and FGF5 in goat embryos. The phenotypes induced by edited loss-of-function mutations of MSTN remain to be evaluated extensively. We demonstrate the utility of this approach by disrupting MSTN, resulting in enhanced body weight and larger muscle fiber size in Cas9-mediated gene-modified goats. The effects of genome modifications were further characterized by H&E staining, quantitative PCR, Western blotting and immunofluorescence staining. Morphological and genetic analyses indicated the occurrence of phenotypic and genotypic modifications. We further provide sufficient evidence, including breeding data, to demonstrate the transmission of the knockout alleles through the germline. By phenotypic and genotypic characterization, we demonstrated the merit of using the CRISPR/Cas9 approach for establishing genetically modified livestock with an enhanced production trait.

Tβ4-overexpression based on the piggyBac transposon system in cashmere goats alters hair fiber characteristics

AbstractIncreasing cashmere yield is one of the vital aims of cashmere goats breeding. Compared to traditional breeding methods, transgenic technology is more efficient and the piggyBac (PB) transposon system has been widely applied to generate transgenic animals. For the present study, donor fibroblasts were stably transfected via a PB donor vector containing the coding sequence of cashmere goat thymosin beta-4 (Tβ4) and driven by a hair follicle-specific promoter, the keratin-associated protein 6.1 (KAP6.1) promoter. To obtain genetically modified cells as nuclear donors, we co-transfected donor vectors into fetal fibroblasts of cashmere goats. Five transgenic cashmere goats were generated following somatic cell nuclear transfer (SCNT). Via determination of the copy numbers and integration sites, the Tβ4 gene was successfully inserted into the goat genome. Histological examination of skin tissue revealed that Tβ4-overexpressing, transgenic goats had a higher secondary to primary hair follicle (S/P) ratio compared to wild type goats. This indicates that Tβ4-overexpressing goats possess increased numbers of secondary hair follicles (SHF). Our results indicate that Tβ4-overexpression in cashmere goats could be a feasible strategy to increase cashmere yield.

RNA-seq reveals transcriptome changes in goats following myostatin gene knockout

Myostatin (MSTN) is a powerful negative regulator of skeletal muscle mass in mammalian species that is primarily expressed in skeletal muscles, and mutations of its encoding gene can result in the double-muscling trait. In this study, the CRISPR/Cas9 technique was used to edit MSTN in Shaanbei Cashmere goats and generate knockout animals. RNA sequencing was used to determine and compare the transcriptome profiles of the muscles from three wild-type (WT) goats, three fibroblast growth factor 5 (FGF5) knockout goats (FGF5+/- group) and three goats with disrupted expression of both the FGF5 and MSTN genes (FM+/- group). The sequence reads were obtained using the Illumina HiSeq 2000 system and mapped to the Capra hircus reference genome using TopHat (v2.0.9). In total, 68.93, 62.04 and 66.26 million clean sequencing reads were obtained from the WT, FM+/- and FGF5+/- groups, respectively. There were 201 differentially expressed genes (DEGs) between the WT and FGF5+/- groups, with 86 down- and 115 up-regulated genes in the FGF5+/- group. Between the WT and FM+/- groups, 121 DEGs were identified, including 81 down- and 40 up-regulated genes in the FM+/- group. A total of 198 DEGs were detected between the FGF5+/- group and FM+/- group, with 128 down- and 70 up-regulated genes in the FM+/- group. At the transcriptome level, we found substantial changes in genes involved in fatty acid metabolism and the biosynthesis of unsaturated fatty acids, such as stearoyl-CoA dehydrogenase, 3-hydroxyacyl-CoA dehydratase 2, ELOVL fatty acid elongase 6 and fatty acid synthase, suggesting that the expression levels of these genes may be directly regulated by MSTN and that these genes are likely downstream targets of MSTN with potential roles in lipid metabolism in goats. Moreover, five randomly selected DEGs were further validated with qRT-PCR, and the results were consistent with the transcriptome analysis. The present study provides insight into the unique transcriptome profile of the MSTN knockout goat, which is a valuable resource for studying goat genomics.

Trio-Based Deep Sequencing Reveals a Low Incidence of Off-Target Mutations in the Offspring of Genetically Edited Goats

Unintended off-target mutations induced by CRISPR/Cas9 nucleases may result in unwanted consequences, which will impede the efficient applicability of this technology for genetic improvement. We have recently edited the goat genome through CRISPR/Cas9 by targeting MSTN and FGF5, which increased muscle fiber diameter and hair fiber length, respectively. Using family trio-based sequencing that allow better discrimination of variant origins, we herein generated offspring from edited goats, and sequenced the members of four family trios (gene-edited goats and their offspring) to an average of ∼36.8× coverage. This data was to systematically examined for mutation profiles using a stringent pipeline that comprehensively analyzed the sequence data for de novo single nucleotide variants, indels, and structural variants from the genome. Our results revealed that the incidence of de novo mutations in the offspring was equivalent to normal populations. We further conducted RNA sequencing using muscle and skin tissues from the offspring and control animals, the differentially expressed genes (DEGs) were related to muscle fiber development in muscles, skin development, and immune responses in skin tissues. Furthermore, in contrast to recently reports of Cas9 triggered p53 expression alterations in cultured cells, we provide primary evidence to show that Cas9-mediated genetic modification does not induce apparent p53 expression changes in animal tissues. This work provides adequate molecular evidence to support the reliability of conducting Cas9-mediated genome editing in large animal models for biomedicine and agriculture.