Yulin Chen

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.

Exploring Differentially Expressed Genes by RNA-Seq in Cashmere Goat (Capra hircus) Skin during Hair Follicle Development and Cycling

Cashmere goat (Capra hircus) hair follicle development and cycling can be divided into three stages: anagen, catagen and telogen. To elucidate the genes involved in hair follicle development and cycling in cashmere goats, transcriptome profiling of skin was carried out by analysing samples from three hair follicle developmental stages using RNA-Seq. The RNA-Seq analysis generated 8487344, 8142514 and 7345335 clean reads in anagen, catagen and telogen stages, respectively, which provided abundant data for further analysis. A total of 1332 differentially expressed genes (DEGs) were identified, providing evidence that the development of hair follicles among the three distinct stages changed considerably. A total of 683 genes with significant differential expression were detected between anagen and catagen, 530 DEGs were identified between anagen and telogen, and 119 DEGs were identified between catagen and telogen. A large number of DEGs were predominantly related to cellular process, cell & cell part, binding, biological regulation and metabolic process among the different stages of hair follicle development. In addition, the Wnt, Shh, TGF-β and Notch signaling pathways may be involved in hair follicle development and the identified DEGs may play important roles in these signaling pathways. These results will expand our understanding of the complex molecular mechanisms of hair follicle development and cycling in cashmere goats and provide a foundation for future studies.

Transcriptome profile analysis of adipose tissues from fat and short-tailed sheep

Recent studies in domestic animals have used RNA-seq to explore the transcriptome of different tissues in a limited number of individuals. In the present study, de novo transcriptome sequencing was used to compare sheep adipose tissue transcriptome profiles between a fat-tailed breed (Kazak sheep; KS) and a short-tailed (Tibetan sheep; TS). The RNA-seq data from these two groups revealed that 646 genes were differentially expressed between the KS and TS groups, including 280 up-regulated and 366 down-regulated genes. We identified genes relevant to fat metabolism in adipose tissues, including two top genes with the largest fold change (NELL1 and FMO3). Pathway analysis revealed that the differentially expressed genes between the KS and TS breeds belong to fatty acid metabolism relevant pathways (e.g. fat digestion and absorption, glycine, serine, and threonine metabolism) and cell junction-related pathways (e.g. cell adhesion molecules) which contribute to fat deposition. This work highlighted potential genes and gene networks that affect fat deposition and meat quality in sheep.

Generation of an artificial double promoter for protein expression in Bacillus subtilis through a promoter trap system.

Bacillus subtilis is an attractive host for production of recombinant proteins. Promoters and expression plasmid backbones have direct impacts on the efficiency of gene expression. To screen and isolate strong promoters, a promoter trap vector pShuttleF was developed in this study. Using the vector, approximately 1000 colonies containing likely promoters from Bacillus licheniformis genomic DNA were obtained. Amongst them, pShuttle-09 exhibited the highest β-Gal activities in both Escherichia coli and B. subtilis. The activity of pShuttle-09 in B. subtilis was eight times of that of the P43 promoter, a commonly used strong promoter for B. subtilis. A sequence analysis showed that pShuttle-09 contained PluxS and truncated luxS in-frame fused with the reporter gene as well as another fragment upstream of PluxS containing a putative promoter. This putative promoter was a hybrid promoter and its β-Gal activity was higher than PluxS. Reconstructing the hybrid promoter from pShuttle-09 to PlapS further improved the β-Gal production by 60%. The usefulness of our promoter trap system is likely due to random shuffling and recombination of DNA fragments and adoption of a rapid and high-throughput screening. Thus, our data provide additional evidence to support the concept of using a promoter trap system to create new promoters.

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.

Multiplex gene editing via CRISPR/Cas9 exhibits desirable muscle hypertrophy without detectable off-target effects in sheep

The CRISPR/Cas9 system provides a flexible approach for genome engineering of genetic loci. Here, we successfully achieved precise gene targeting in sheep by co-injecting one-cell-stage embryos with Cas9 mRNA and RNA guides targeting three genes (MSTN, ASIP, and BCO2). We carefully examined the sgRNAs:Cas9-mediated targeting effects in injected embryos, somatic tissues, as well as gonads via cloning and sequencing. The targeting efficiencies in these three genes were within the range of 27–33% in generated lambs, and that of simultaneously targeting the three genes was 5.6%, which demonstrated that micro-injection of zygotes is an efficient approach for generating gene-modified sheep. Interestingly, we observed that disruption of the MSTN gene resulted in the desired muscle hypertrophy that is characterized by enlarged myofibers, thereby providing the first detailed evidence supporting that gene modifications had occurred at both the genetic and morphological levels. In addition, prescreening for the off-target effect of sgRNAs was performed on fibroblasts before microinjection, to ensure that no detectable off-target mutations from founder animals existed. Our findings suggested that the CRISPR/Cas9 method can be exploited as a powerful tool for livestock improvement by simultaneously targeting multiple genes that are responsible for economically significant traits.

Whole-genome sequencing of eight goat populations for the detection of selection signatures underlying production and adaptive traits.

The goat (Capra hircus) is one of the first farm animals that have undergone domestication and extensive natural and artificial selection by adapting to various environments, which in turn has resulted in its high level of phenotypic diversity. Here, we generated medium-coverage (9–13×) sequences from eight domesticated goat breeds, representing morphologically or geographically specific populations, to identify genomic regions representing selection signatures. We discovered ~10 million single nucleotide polymorphisms (SNPs) for each breed. By combining two approaches, ZHp and di values, we identified 22 genomic regions that may have contributed to the phenotypes in coat color patterns, body size, cashmere traits, as well as high altitude adaptation in goat populations. Candidate genes underlying strong selection signatures including coloration (ASIP, KITLG, HTT, GNA11, and OSTM1), body size (TBX15, DGCR8, CDC25A, and RDH16), cashmere traits (LHX2, FGF9, and WNT2), and hypoxia adaptation (CDK2, SOCS2, NOXA1, and ENPEP) were identified. We also identified candidate functional SNPs within selected genes that may be important for each trait. Our results demonstrated the potential of using sequence data in identifying genomic regions that are responsible for agriculturally significant phenotypes in goats, which in turn can be used in the selection of goat breeds for environmental adaptation and domestication.

Comparative Transcriptome Analysis of Fetal Skin Reveals Key Genes Related to Hair Follicle Morphogenesis in Cashmere Goats.

Cashmere goat skin contains two types of hair follicles (HF): primary hair follicles (PHF) and secondary hair follicles (SHF). Although multiple genetic determinants associated with HF formation have been identified, the molecules that determine the independent morphogenesis of HF in cashmere goats remain elusive. The growth and development of SHF directly influence the quantity and quality of cashmere production. Here, we report the transcriptome profiling analysis of nine skin samples from cashmere goats using 60- and 120-day-old embryos (E60 and E120, respectively), as well as newborns (NB), through RNA-sequencing (RNA-seq). HF morphological changes indicated that PHF were initiated at E60, with maturation from E120, while differentiation of SHF was identified at E120 until formation of cashmere occurred after birth (NB). The RNA-sequencing analysis generated over 20.6 million clean reads from each mRNA library. The number of differentially expressed genes (DEGs) in E60 vs. E120, E120 vs. NB, and E60 vs. NB were 1,024, 0 and 1,801, respectively, indicating that no significant differences were found at transcriptomic levels between E120 and NB. Key genes including B4GALT4, TNC, a-integrin, and FGFR1, were up-regulated and expressed in HF initiation from E60 to E120, while regulatory genes such as GPRC5D, PAD3, HOXC13, PRR9, VSIG8, LRRC15, LHX2, MSX-2, and FOXN1 were up-regulated and expressed in HF keratinisation and hair shaft differentiation from E120 and NB to E60. Several genes belonging to the KRT and KRTAP gene families were detected throughout the three HF developmental stages. The transcriptional trajectory analyses of all DEGs indicated that immune privilege, glycosaminoglycan biosynthesis, extracellular matrix receptor interaction, and growth factor receptors all played dominant roles in the epithelial-mesenchymal interface and HF formation. We found that the Wnt, transforming growth factor-beta/bone morphogenetic protein, and Notch family members played vital roles in HF differentiation and maturation. The DEGs we found could be attributed to the generation and development of HF, and thus will be critically important for improving the quantity and quality of fleece production in animals for fibres.

Fibroblast growth factor 5-short (FGF5s) inhibits the activity of FGF5 in primary and secondary hair follicle dermal papilla cells of cashmere goats.

To determine the relationship between fibroblast growth factor 5 (FGF5) and FGF5-short (FGF5s) in dermal papilla cells of cashmere goat primary and secondary hair follicles. We isolated dermal papilla cells from primary hair follicle (PHF) and secondary hair follicle (SHF) of cashmere goat, and found that the FGF5 receptor, fibroblast growth factor receptor 1 (FGFR1), was expressed in these two types of dermal papilla cells. Moreover, adenovirus-mediated overexpression of FGF5 could upregulate the mRNA expression of insulin-like growth factor-1 (IGF-1), versican and noggin that were important for follicle growth maintenance, whereas downregulate the expression of anagen chalone bone morphogenetic protein 4 (BMP4) in dermal papilla cells. However, these alterations were partly reversed by FGF5s overexpression. In conclusion, our results demonstrated that FGF5s acted as an inhibitor of FGF5 in the regulation of anagen-catagen transition of cashmere goat dermal papilla cells.