ZBED6, a novel transcription factor derived from a domesticated DNA transposon regulates IGF2 expression and muscle growth

Circular RNA transcriptome across various development periods of Paralichthys olivaceus reveal skeletal muscle-specific circchd6 regulating myogenesis

The Japanese flounder (Paralichthys olivaceus) is greatly influenced in terms of muscle quality and quantity by the development of skeletal muscle. While the mechanisms underlying skeletal muscle development are well-studied, the role of non-coding RNAs, particularly circRNAs, in the skeletal muscle development of Japanese flounder remains unclear. To investigate the expression patterns of circRNAs during different developmental stages (JP1: 7 days, JP2: 90 days, JP3: 24 months (female), JP4: 24 months (male)) in Japanese flounder, we performed transcriptome sequencing analysis. We identified a total of 3523 circRNAs, of which 10.19 % were differentially expressed. These differentially expressed (DE) circRNAs were studied, and their impacts on muscle development were analyzed. The RNA interaction network revealed that skeletal muscle-specific circchd6 targeted novel-miR-508 and further regulated dual specificity tyrosine-phosphorylation regulated kinase 2 (dyrk2). Functional analysis showed that overexpressed circchd6 and dyrk2 promoted myoblast proliferation and differentiation, while novel-miR-508 inhibited both. Our study identified the circchd6-novel-miR-508-dyrk2 axis as a regulatory mechanism and provided new evidence for the use of epigenetic approaches in genetic breeding.

Spatial 3D genome organization reveals intratumor heterogeneity in primary glioblastoma samples

Glioblastoma (GBM) is the most prevalent malignant brain tumor with poor prognosis. Although chromatin intratumoral heterogeneity is a characteristic feature of GBM, most current studies are conducted at a single tumor site. To investigate the GBM-specific 3D genome organization and its heterogeneity, we conducted Hi-C experiments in 21 GBM samples from nine patients, along with three normal brain samples. We identified genome subcompartmentalization and chromatin interactions specific to GBM, as well as extensive intertumoral and intratumoral heterogeneity at these levels. We identified copy number variants (CNVs) and structural variations (SVs) and demonstrated how they disrupted 3D genome structures. SVs could not only induce enhancer hijacking but also cause the loss of enhancers to the same gene, both of which contributed to gene dysregulation. Our findings provide insights into the GBM-specific 3D genome organization and the intratumoral heterogeneity of this organization and open avenues for understanding this devastating disease.

Targeting ZC3H11A elicits immunogenic cancer cell death through augmentation of antigen presentation and interferon response

Zinc finger CCCH containing 11A (ZC3H11A) is a stress-induced protein that is upregulated in various conditions such as heat shock and virus infection. It has also been reported to be upregulated in certain cancers. The aim of this study was to evaluate the feasibility of targeting ZC3H11A as a therapeutic approach for cancer treatment, using nuclease-resistant, affinity-enhanced antisense oligonucleotide (ASO). An ASO targeting ZC3H11A was validated and evaluated in vitro and in the B16 melanoma model in vivo. Antigen presentation, interferon response, cell proliferation, and apoptosis were transcriptionally affected. These findings were validated on the protein level by the upregulation of major histocompatibility complex class I (MHC class I), an increased secretion of interferon-β (IFN-β), and induction of apoptosis observed as upregulation of caspases and annexin V. Immunogenic features of the induced apoptosis were evidenced by the surface exposure of calreticulin (CRT) and the secretion of ATP leading to enhanced dendritic cell (DC) phagocytosis, maturation, and activation. Treatment with the ZC3H11A-targeted ASO had limited efficacy in vivo, while constitutive lentiviral shRNA knockdown of ZC3H11A in murine B16 melanoma cells and human HeLa cells led to reduced tumor growth with prolonged survival of mice, validating ZC3H11A as a relevant target for cancer therapy.

Tumor-suppressive activities for pogo transposable element derived with KRAB domain via ribosome biogenesis restriction

Transposable elements (TEs) are indispensable for human development, with critical functions in pluripotency and embryogenesis. TE sequences also contribute to human pathologies, especially cancer, with documented activities as cis/trans transcriptional regulators, as sources of non-coding RNAs, and as mutagens that disrupt tumor suppressors. Despite this knowledge, little is known regarding the involvement of TE-derived genes (TEGs) in tumor pathogenesis. Here, systematic analyses of TEG expression across human cancer reveal a prominent role for pogo TE derived with KRAB domain (POGK). We show that POGK acts as a tumor suppressor in triple-negative breast cancer (TNBC) cells and that it couples with the co-repressor TRIM28 to directly block the transcription of ribosomal genes RPS16 and RPS29, in turn causing widespread inhibition of ribosomal biogenesis. We report that POGK undergoes deactivation by isoform switching in clinical TNBC, altogether revealing its exapted activities in tumor growth control.

Genome-wide identification and characterization of SLEEPER, a transposon-derived gene family and their expression pattern in Brassica napus L

BACKGROUND

The transposons of the hAT superfamily are the most widespread transposons ever known. SLEEPER genes encode domesticated transposases from the hAT superfamily, which may have lost their transposable functions during long-term evolution and transformed into host proteins that regulate plant growth and development.

RESULTS

This study identified 162 members of the SLEEPER gene family from Brassica napus. These members are widely distributed on 19 chromosomes, mainly in the Cn subgenome, and have promoters with various cis-acting elements related to hormone regulation, abiotic stress, and growth and development regulation. Most of the genes in this family contain similar conserved domains and motifs, and the closer the genes are distributed on evolutionary branches, the more similar their structures are. Transcriptome sequencing performed on tissues at different growth stages from B. napus line 3529 indicated that these genes had different expression patterns, and nearly half of the genes were not detectably expressed in all samples.

CONCLUSIONS

This study investigated the gene structure, expression patterns, evolutionary features, and gene localization of the SLEEPER family members to confirm the significance of these genes in the growth of B. napus, providing a reference for the study of transposon domestication and outstanding genetic resources for the genetic improvement of B. napus.

The hAT family hopper transposon exists as highly similar yet discontinuous elements in the Bactrocera tephritid fly genus

The hAT family transposable element, hopper, was originally discovered as a defective 3120-bp full-length element in a wild-type strain of the oriental fruit fly, Bactrocera dorsalis (Hendel) (Diptera: Tephritidae), and subsequently a functional 3131-bp element, hopperBdwe, was isolated from a white eye mutant strain. The latter study showed that closely related elements exist in melonfly, Zeugodacus cucurbitae (Coquillett) (Diptera: Tephritidae), a closely related subgenus, suggesting that hopper could have a widespread presence in the Bactrocera genus. To further understand the distribution of hopper within and beyond the B. dorsalis species complex, primer pairs from hopperBdwe and its adjacent genomic insertion site were used to survey the presence and relatedness of hopper in five species within the complex and four species beyond the complex. Based on sequence identity of a 1.94 kb internal nucleotide sequence, the closest relationships were with mutated elements from B. dorsalis s.s. and species synonymized with B. dorsalis including B. papayae, B. philippinensis and B. invadens, ranging in identity between 88.4% and 99.5%. Notably, Bactrocera carambolae (Drew & Hancock) (Diptera: Tephritidae), which is most closely related to B. dorsalis beyond the synonymized species, shared hopper identities of 97.3%-99.5%. Beyond the B. dorsalis complex, Z. cucurbitae, Bactrocera tryoni (Froggatt) (Diptera: Tephritidae) and Bactrocera zonata (Saunders) (Diptera: Tephritidae) shared identities of 83.1%-97.1%, while hopper was absent from the Bactrocera oleae (Gmelin) (Diptera: Tephritidae) strain tested. While the functional autonomous hopperBdwe element was not detected in these species, another closely related hopper element isolated from a B. dorsalis genetic sexing strain has an uninterrupted transposase open reading frame. The discontinuous presence of hopper in the Bactrocera genus has implications for its use for genomic manipulation and understanding the phylogenetic relationship of these species.

Identification and functional analysis of circpdlim5a generated from pdlim5a gene splicing in the skeletal muscle of Japanese flounder (Paralichthys olivaceus)

Circular RNAs (circRNAs) are non-coding RNAs with endogenous regulatory functions, including regulating skeletal muscle development. However, its role in the development of skeletal muscle in Japanese flounder (Paralichthys olivaceus) is not clear. Therefore we screened a candidate circpdlim5a, which is derived from the gene pdlim5a, from the skeletal muscle transcriptome of Japanese flounder. We characterized circpdlim5a, which was more stable compared to the linear RNA pdlim5a. Distributional characterization of circpdlim5a showed that circpdlim5a was predominantly distributed in the nucleus and was highly expressed in the skeletal muscle of adult Japanese flounder (24 months). When we further studied the circpdlim5a function, we found that it inhibited the expression of proliferation and differentiation genes according to the over-expression experiment of circpdlim5a in myoblasts. We concluded that circpdlim5a may inhibit the proliferation and differentiation of myoblasts and thereby inhibit skeletal muscle development in Japanese flounder. This experiment provides information for the study of circRNAs by identifying circpdlim5a and exploring its function, and offers clues for molecular breeding from an epigenetic perspective.

Three Copies of zbed1 Specific in Chromosome W Are Essential for Female-Biased Sexual Size Dimorphism in Cynoglossus semilaevis

The sex chromosome, especially specific in one sex, generally determines sexual size dimorphism (SSD), a phenomenon with dimorphic sexual difference in the body size. For Cynoglossus semilaevis, a flatfish in China, although the importance of chromosome W and its specific gene zbed1 in female-biased SSD have been suggested, its family members and regulation information are still unknown. At present, three zbed1 copies gene were identified on chromosome W, with no gametologs. Phylogenetic analysis for the ZBED family revealed an existence of ZBED9 in the fish. Nine members were uncovered from C. semilaevis, clustering into three kinds, ZBED1, ZBED4 and ZBEDX, which is less than the eleven kinds of ZBED members in mammals. The predominant expression of zbed1 in the female brain and pituitary tissues was further verified by qPCR. Transcription factor c/ebpα could significantly enhance the transcriptional activity of zbed1 promoter, which is opposite to its effect on the male determinant factor-dmrt1. When zbed1 was interfered with, piwil1, esr2 and wnt7b were up-regulated, while cell-cycle-related genes, including cdk4 and ccng1, were down-regulated. Thus, zbed1 is involved in cell proliferation by regulating esr2, piwil1, cell cycle and the Wnt pathway. Further research on their interactions would be helpful to understand fish SSD.

Livestock species as emerging models for genomic imprinting

Genomic imprinting is an epigenetically-regulated process of central importance in mammalian development and evolution. It involves multiple levels of regulation, with spatio-temporal heterogeneity, leading to the context-dependent and parent-of-origin specific expression of a small fraction of the genome. Genomic imprinting studies have therefore been essential to increase basic knowledge in functional genomics, evolution biology and developmental biology, as well as with regard to potential clinical and agrigenomic perspectives. Here we offer an overview on the contribution of livestock research, which features attractive resources in several respects, for better understanding genomic imprinting and its functional impacts. Given the related broad implications and complexity, we promote the use of such resources for studying genomic imprinting in a holistic and integrative view. We hope this mini-review will draw attention to the relevance of livestock genomic imprinting studies and stimulate research in this area.

The polymorphism of the ovine insulin like growth factor-2 (IGF2) gene and their associations with growth related traits in Tibetan sheep

In the current study, the role of the ovine IGF2 as a potential candidate gene was investigated as though marker-assisted selection in Chinese Tibetan sheep. The Sanger DNA sequencing method explored five single nucleotide polymorphisms (SNPs) in 5'UTR of the ovine IGF2 gene (C15640T, G15801A, G15870A, C15982G and G15991A) in Chinese Tibetan sheep. The frequencies of four SNPs were within the Hardy-Weinberg Equilibrium (chi-square test) except C15982G. The statistical analysis indicated that the C15640T and G15801A were significantly associated with body height, body length, chest circumference, and body weight (P < 0.05 or P < 0.01). Furthermore, C15982G variant exhibited significant correlation with the body weight (P < 0.01). These findings suggests that the promoter variants of IGF2 gene could be used as a candidate gene through marker-assisted selection for the body weight and body measurement traits in Tibetan sheep breeding program.

Single-base editing in IGF2 improves meat production and intramuscular fat deposition in Liang Guang Small Spotted pigs

BACKGROUND

Chinese indigenous pigs are popular with consumers for their juiciness, flavour and meat quality, but they have lower meat production. Insulin-like growth factor 2 (IGF2) is a maternally imprinted growth factor that promotes skeletal muscle growth by regulating cell proliferation and differentiation. A single nucleotide polymorphism (SNP) within intron 3 of porcine IGF2 disrupts a binding site for the repressor, zinc finger BED-type containing 6 (ZBED6), leading to up-regulation of IGF2 and causing major effects on muscle growth, heart size, and backfat thickness. This favorable mutation is common in Western commercial pig populations, but absent in most Chinese indigenous pig breeds. To improve meat production of Chinese indigenous pigs, we used cytosine base editor 3 (CBE3) to introduce IGF2-intron3-C3071T mutation into porcine embryonic fibroblasts (PEFs) isolated from a male Liang Guang Small Spotted pig (LGSS), and single-cell clones harboring the desired mutation were selected for somatic cell nuclear transfer (SCNT) to generate the founder line of IGF2T/T pigs.

RESULTS

We found the heterozygous progeny IGF2C/T pigs exhibited enhanced expression of IGF2, increased lean meat by 18%-36%, enlarged loin muscle area by 3%-17%, improved intramuscular fat (IMF) content by 18%-39%, marbling score by 0.75-1, meat color score by 0.53-1.25, and reduced backfat thickness by 5%-16%. The enhanced accumulation of intramuscular fat in IGF2C/T pigs was identified to be regulated by the PI3K-AKT/AMPK pathway, which activated SREBP1 to promote adipogenesis.

CONCLUSIONS

We demonstrated the introduction of IGF2-intron3-C3071T in Chinese LGSS can improve both meat production and quality, and first identified the regulation of IMF deposition by IGF2 through SREBP1 via the PI3K-AKT/AMPK signaling pathways. Our study provides a further understanding of the biological functions of IGF2 and an example for improving porcine economic traits through precise base editing.

A Functional Variant Alters the Binding of Bone morphogenetic protein 2 to the Transcription Factor NF-κB to Regulate Bone morphogenetic protein 2 Gene Expression and Chicken Abdominal Fat Deposition

In this study, we employed a dual-luciferase reporter assay and electrophoretic mobility shift analysis (EMSA) in vitro to explore whether a 12-base pair (bp) insertion/deletion (InDel) variant (namely g.14798187_14798188insTCCCTGCCCCCT) within intron 2 of the chicken BMP2 gene, which was significantly associated with chicken abdominal fat weight and abdominal fat percentage, is a functional marker and its potential regulatory mechanism. The reporter analysis demonstrated that the luciferase activity of the deletion allele was extremely significantly higher than that of the insertion allele (p < 0.01). A bioinformatics analysis revealed that compared to the deletion allele, the insertion allele created a transcription factor binding site of nuclear factor-kappa B (NF-κB), which exhibited an inhibitory effect on fat deposition. A dual-luciferase reporter assay demonstrated that the inhibitory effect of NF-κB on the deletion allele was stronger than that on the insertion allele. EMSA indicated that the binding affinity of NF-κB for the insertion allele was stronger than that for the deletion allele. In conclusion, the 12-bp InDel chicken BMP2 gene variant is a functional variant affecting fat deposition in chickens, which may partially regulate BMP2 gene expression by affecting the binding of transcription factor NF-κB to the BMP2 gene.

An organism-wide ATAC-seq peak catalog for the bovine and its use to identify regulatory variants

We report the generation of an organism-wide catalog of 976,813 cis-acting regulatory elements for the bovine detected by the assay for transposase accessible chromatin using sequencing (ATAC-seq). We regroup these regulatory elements in 16 components by nonnegative matrix factorization. Correlation between the genome-wide density of peaks and transcription start sites, correlation between peak accessibility and expression of neighboring genes, and enrichment in transcription factor binding motifs support their regulatory potential. Using a previously established catalog of 12,736,643 variants, we show that the proportion of single-nucleotide polymorphisms mapping to ATAC-seq peaks is higher than expected and that this is owing to an approximately 1.3-fold higher mutation rate within peaks. Their site frequency spectrum indicates that variants in ATAC-seq peaks are subject to purifying selection. We generate eQTL data sets for liver and blood and show that variants that drive eQTL fall into liver- and blood-specific ATAC-seq peaks more often than expected by chance. We combine ATAC-seq and eQTL data to estimate that the proportion of regulatory variants mapping to ATAC-seq peaks is approximately one in three and that the proportion of variants mapping to ATAC-seq peaks that are regulatory is approximately one in 25. We discuss the implication of these findings on the utility of ATAC-seq information to improve the accuracy of genomic selection.

Loss of ZBED6 Protects Against Sepsis-Induced Muscle Atrophy by Upregulating DOCK3-Mediated RAC1/PI3K/AKT Signaling Pathway in Pigs

Sepsis-induced muscle atrophy often increases morbidity and mortality in intensive care unit (ICU) patients, yet neither therapeutic target nor optimal animal model is available for this disease. Here, by modifying the surgical strategy of cecal ligation and puncture (CLP), a novel sepsis pig model is created that for the first time recapitulates the whole course of sepsis in humans. With this model and sepsis patients, increased levels of the transcription factor zinc finger BED-type containing 6 (ZBED6) in skeletal muscle are shown. Protection against sepsis-induced muscle wasting in ZBED6-deficient pigs is further demonstrated. Mechanistically, integrated analysis of RNA-seq and ChIP-seq reveals dedicator of cytokinesis 3 (DOCK3) as the direct target of ZBED6. In septic ZBED6-deficient pigs, DOCK3 expression is increased in skeletal muscle and myocytes, activating the RAC1/PI3K/AKT pathway and protecting against sepsis-induced muscle wasting. Conversely, opposite gene expression patterns and exacerbated muscle wasting are observed in septic ZBED6-overexpressing myotubes. Notably, sepsis patients show increased ZBED6 expression along with reduced DOCK3 and downregulated RAC1/PI3K/AKT pathway. These findings suggest that ZBED6 is a potential therapeutic target for sepsis-induced muscle atrophy, and the established sepsis pig model is a valuable tool for understanding sepsis pathogenesis and developing its therapeutics.

Transcriptomic analysis of glutamate-induced HT22 neurotoxicity as a model for screening anti-Alzheimer's drugs

Glutamate-induced neurotoxicity in the HT22 mouse hippocampal neuronal cell line has been recognized as a valuable cell model for the study of neurotoxicity associated with neurodegenerative diseases including Alzheimer's disease (AD). However, the relevance of this cell model for AD pathogenesis and preclinical drug screening remains to be more elucidated. While there is increasing use of this cell model in a number of studies, relatively little is known about its underlying molecular signatures in relation to AD. Here, our RNA sequencing study provides the first transcriptomic and network analyses of HT22 cells following glutamate exposure. Several differentially expressed genes (DEGs) and their relationships specific to AD were identified. Additionally, the usefulness of this cell model as a drug screening system was assessed by determining the expression of those AD-associated DEGs in response to two medicinal plant extracts, Acanthus ebracteatus and Streblus asper, that have been previously shown to be protective in this cell model. In summary, the present study reports newly identified AD-specific molecular signatures in glutamate-injured HT22 cells, suggesting that this cell can be a valuable model system for the screening and evaluation of new anti-AD agents, particularly from natural products.

Effect of Zbed6 Single-Allele Knockout on the Growth and Development of Skeletal Muscle in Mice

ZBED6, a key transcription factor, plays an important role in skeletal muscle and organ growth. ZBED6 knockout (ZBED6^-/-) leads to the upregulation of IGF2 in pig and mice muscle, thereby increasing muscle mass. However, the effects and mechanism of Zbed6 single-allele knockout (Zbed6^+/-) on mice muscle remain unknown. Here, we reported that Zbed6^+/- promotes muscle growth by a new potential target gene rather than Igf2 in mice muscle. Zbed6^+/- mice showed markedly higher muscle mass (25%) and a markedly higher muscle weight ratio (18%) than wild-type (WT) mice, coinciding with a larger muscle fiber area (28%). Despite a significant increase in muscle growth, Zbed6^+/- mice showed similar Igf2 expression with WT mice, indicating that a ZBED6-Igf2-independent regulatory pathway exists in Zbed6^+/- mice muscle. RNA-seq of muscle between the Zbed6^+/- and WT mice revealed two terms related to muscle growth. Overlapping the DEGs and C2C12 Chip-seq data of ZBED6 screened out a potential ZBED6 target gene Barx2, which may regulate muscle growth in Zbed6^+/- mice. These results may open new research directions leading to a better understanding of the integral functions of ZBED6 and provide evidence of Zbed6^+/- promoting muscle growth by regulating Barx2 in mice.

Genomic Structure, Protein Character, Phylogenic Implication, and Embryonic Expression Pattern of a Zebrafish New Member of Zinc Finger BED-Type Gene Family

We reported a new member of the C₂H₂-zinc-finger BED-type (ZBED) protein family found in zebrafish (Danio rerio). It was previously assigned as an uncharacterized protein LOC569044 encoded by the Zgc:161969 gene, the transcripts of which were highly expressed in the CNS after the spinal cord injury of zebrafish. As such, this novel gene deserves a more detailed investigation. The 2.79-kb Zgc:161969 gene contains one intron located on Chromosome 6 at 16,468,776-16,475,879 in the zebrafish genome encoding a 630-aa protein LOC569044. This protein is composed of a DNA-binding BED domain, which is highly conserved among the ZBED protein family, and a catalytic domain consisting of an α-helix structure and an hAT dimerization region. Phylogenetic analysis revealed the LOC569044 protein to be clustered into the monophyletic clade of the ZBED protein family of golden fish. Specifically, the LOC569044 protein was classified as closely related to the monophyletic clades of zebrafish ZBED4-like isoforms and ZBED isoform 2. Furthermore, Zgc:161969 transcripts represented maternal inheritance, expressed in the brain and eyes at early developmental stages and in the telencephalon ventricular zone at late developmental stages. After characterizing the LOC569044 protein encoded by the Zgc:161969 gene, it was identified as a new member of the zebrafish ZBED protein family, named the ZBEDX protein.

Transposon-derived transcription factors across metazoans

Transposable elements (TE) could serve as sources of new transcription factors (TFs) in plants and some other model species, but such evidence is lacking for most animal lineages. Here, we discovered multiple independent co-options of TEs to generate 788 TFs across Metazoa, including all early-branching animal lineages. Six of ten superfamilies of DNA transposon-derived conserved TF families (ZBED, CENPB, FHY3, HTH-Psq, THAP, and FLYWCH) were identified across nine phyla encompassing the entire metazoan phylogeny. The most extensive convergent domestication of potentially TE-derived TFs occurred in the hydroid polyps, polychaete worms, cephalopods, oysters, and sea slugs. Phylogenetic reconstructions showed species-specific clustering and lineage-specific expansion; none of the identified TE-derived TFs revealed homologs in their closest neighbors. Together, our study established a framework for categorizing TE-derived TFs and informing the origins of novel genes across phyla.

Indel mutations of sheep PLAG1 gene and their associations with growth traits

Pleiomorphic adenoma gene 1 (PLAG1) is mainly expressed in embryonic development, and it is reported to take an effect on the growth performance of mice, cattle, pigs, and sheep. To explore how conservative the PLAG1 is in different sheep breeds, the effects of the two indel variants on the growth traits of the Chinese Luxi blackhead (LXBH) sheep were firstly detected. The P2-del 30 bp and P4-del 45 bp indel loci of the sheep PLAG1 gene were significantly related to 15 growth traits (P < 0.05). Genotype ID and genotype II were dominant for the P2-del 30 bp and P4-del 45 bp loci, respectively. The above findings indicated that the two indel mutations in the ovine PLAG1 gene were suggested to become the molecular markers for the selection of economic traits in sheep.

Identification and Validation of Ferroptosis-Related Genes in Sevoflurane-Induced Hippocampal Neurotoxicity

Background

Sevoflurane is one of the most popular inhalational anesthetics during perioperative period but presenting neurotoxicity among pediatric and aged populations. Recent experiments in vivo and in vitro have indicated that ferroptosis may contribute to the neurotoxicity of sevoflurane anesthesia. However, the exact mechanism is still unclear.

Methods

In current study, we explored the differential expressed genes (DEGs) in HT-22 mouse hippocampal neuronal cells after sevoflurane anesthesia using RNA-seq. Differential expressed ferroptosis-related genes (DEFRGs) were screened and analyzed by Gene Ontology (GO) and pathway enrichment analysis. Protein-to-protein interaction (PPI) network was constructed by the Search Tool for the Retrieval of Interacting Genes (STRING). Significant modules and the hub genes were identified by using Cytoscape. The Connectivity Map (cMAP) was used for screening drug candidates targeting the identified DEFRGs. Potential TF-gene network and drug-gene pairs were established towards the hub genes. In final, we validated these results in experiments.

Results

A total of 37 ferroptosis-related genes (18 upregulated and 19 downregulated) after sevoflurane exposure in hippocampal neuronal cells were finally identified. These differentially expressed genes were mainly involved into the biological processes of cellular response to oxidative stress. Pathway analysis indicated that these genes were involved in ferroptosis, mTOR signaling pathway, and longevity-regulating pathway. PPI network was constructed. 10 hub genes including Prkaa2, Chac1, Arntl, Tfrc, Slc7a11, Atf4, Mgst1, Lpin1, Atf3, and Sesn2 were found. Top 10 drug candidates, gene-drug networks, and TFs targeting these genes were finally identified. These results were validated in experiments.

Conclusion

Our results suggested that ferroptosis-related genes play roles in sevoflurane anesthesia-related hippocampal neuron injury and offered the hub genes and potential therapeutic agents for investigating and treatment of this neurotoxicity after sevoflurane exposure. Finally, therapeutic effect of these drug candidates and function of potential ferroptosis targets should be further investigated for treatment and clarifying mechanisms of sevoflurane anesthesia-induced neuron injury in future research.

A Comprehensive Analysis of the Myocardial Transcriptome in ZBED6-Knockout Bama Xiang Pigs

The ZBED6 gene is a transcription factor that regulates the expression of IGF2 and affects muscle growth and development. However, its effect on the growth and development of the heart is still unknown. Emerging evidence suggests that long noncoding RNAs (lncRNAs) can regulate genes at the epigenetic, transcriptional, and posttranscriptional levels and play an important role in the development of eukaryotes. To investigate the function of ZBED6 in the cardiac development of pigs, we constructed the expression profiles of mRNAs and lncRNAs in myocardial tissue obtained from Bama Xiang pigs in the ZBED6 knockout group (ZBED6-KO) and the wild-type group (ZBED6-WT). A total of 248 differentially expressed genes (DEGs) and 209 differentially expressed lncRNAs (DELs) were detected, and 105 potential cis target genes of DELs were identified. The functional annotation analysis based on the Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) databases revealed two GO items related to muscle development by the cis target genes of DELs. Moreover, IGF2 was the direct target gene of ZBED6 by ChIP-PCR experiment. Our results explored the mechanism and expression profile of mRNAs and lncRNAs of ZBED6 gene knockout on myocardium tissue development, mining the key candidate genes in that process like IGF2.

Hsa-miR-422a Originated from Short Interspersed Nuclear Element Increases ARID5B Expression by Collaborating with NF-E2

MicroRNAs (miRNAs) are a class of small non-coding RNAs that regulate the expression of target messenger RNA (mRNA) complementary to the 3' untranslated region (UTR) at the post-transcriptional level. Hsa-miR-422a, which is commonly known as miRNA derived from transposable element (MDTE), was derived from short interspersed nuclear element (SINE). Through expression analysis, hsa-miR-422a was found to be highly expressed in both the small intestine and liver of crab-eating monkey. AT-Rich Interaction Domain 5 B (ARID5B) was selected as the target gene of hsa-miR-422a, which has two binding sites in both the exon and 3'UTR of ARID5B. To identify the interaction between hsa-miR-422a and ARID5B, a dual luciferase assay was conducted in HepG2 cell line. The luciferase activity of cells treated with the hsa-miR-422a mimic was upregulated and inversely downregulated when both the hsa-miR-422a mimic and inhibitor were administered. Nuclear factor erythroid-2 (NF-E2) was selected as the core transcription factor (TF) via feed forward loop analysis. The luciferase expression was downregulated when both the hsa-miR-422a mimic and siRNA of NF-E2 were treated, compared to the treatment of the hsa-miR-422a mimic alone. The present study suggests that hsa-miR-422a derived from SINE could bind to the exon region as well as the 3'UTR of ARID5B. Additionally, hsa-miR-422a was found to share binding sites in ARID5Bwith several TFs, including NF-E2. The hsa-miR-422a might thus interact with TF to regulate the expression of ARID5B, as demonstrated experimentally. Altogether, hsa-miR-422a acts as a super enhancer miRNA of ARID5Bby collaborating with TF and NF-E2.

Molecular genetic variation of animals and plants under domestication

Domesticated plants and animals played crucial roles as models for evolutionary change by means of natural selection and for establishing the rules of inheritance, originally proposed by Charles Darwin and Gregor Mendel, respectively. Here, we review progress that has been made during the last 35 y in unraveling the molecular genetic variation underlying the stunning phenotypic diversity in crops and domesticated animals that inspired Mendel and Darwin. We notice that numerous domestication genes, crucial for the domestication process, have been identified in plants, whereas animal domestication appears to have a polygenic background with no obvious "domestication genes" involved. Although model organisms, such as Drosophila and Arabidopsis, have replaced domesticated species as models for basic research, the latter are still outstanding models for evolutionary research because phenotypic change in these species represents an evolutionary process over thousands of years. A consequence of this is that some alleles contributing to phenotypic diversity have evolved by accumulating multiple changes in the same gene. The continued molecular characterization of crops and farm animals with ever sharper tools is essential for future food security.

Efficient Editing of the ZBED6-Binding Site in Intron 3 of IGF2 in a Bovine Model Using the CRISPR/Cas9 System

Background: Insulin-like growth factor 2 is a growth-promoting factor that plays an important role in the growth and development of mammals. A nucleotide substitution in intron 3 of IGF2—which disrupts the ZBED6-binding site—affects muscle mass, organ size, and fat deposition in pigs. The ZBED6-binding site is also conserved in cattle. Methods: In the present study, we introduced mutations in the ZBED6-binding site in intron3 of IGF2 in bovine fetal fibroblasts using the CRISPR/Cas9 system, and investigated the effect of disruption of ZBED6 binding on IGF2 expression. Results: Eleven biallelic-mutant single-cell clones were established, three of which contained no foreign DNA residues. Single-cell clones 93 and 135 were used to produce cloned embryos. Dual-luciferase reporter assay in C2C12 cells demonstrated that the mutation in the ZBED6-binding site increases the promoter 3 activity of bovine IGF2. A total of 49 mutant cloned embryos were transplanted into surrogate cows. Unfortunately, all cloned embryos died before birth. IGF2 was found to be hypomethylated in the only fetus born (stillborn), which may have been due to the incomplete reprogramming. Conclusions: We efficiently constructed IGF2-edited cell lines and cloned embryos, which provided a theoretical basis and experimental materials for beef cattle breeding.

Porcine ZBED6 regulates growth of skeletal muscle and internal organs via multiple targets

ZBED6 (zinc finger BED domain containing protein 6) is a transcription factor unique to placental mammals and its interaction with the IGF2 (insulin-like growth factor 2) locus plays a prominent role in the regulation of postnatal skeletal muscle growth. Here, we generated lean Bama miniature pigs by generating ZBED6-knockout (ZBED6^−/−) and investigated the mechanism underlying ZBED6 in growth of muscle and internal organs of placental mammals. ZBED6^−/− pigs show markedly higher lean mass, lean mass rate, larger muscle fiber area and heavier internal organs (heart and liver) than wild-type (WT) pigs. The striking phenotypic changes of ZBED6^-/- pigs coincided with remarkable upregulation of IGF2 mRNA and protein expression across three tissues (gastrocnemius muscle, longissimus dorsi, heart). Despite a significant increase in liver weight, ZBED6^-/- pigs show comparable levels of IGF2 expression to those of WT controls. A mechanistic study revealed that elevated methylation in the liver abrogates ZBED6 binding at the IGF2 locus, explaining the unaltered hepatic IGF2 expression in ZBED6^-/- pigs. These results indicate that a ZBED6-IGF2-independent regulatory pathway exists in the liver. Transcriptome analysis and ChIP-PCR revealed new ZBED6 target genes other than IGF2, including cyclin dependent kinase inhibitor 1A (CDKN1A) and tsukushi, small leucine rich proteoglycan (TSKU), that regulates growth of muscle and liver, respectively.

The lean meat rate is an important economic trait for the swine industry and it is determined by muscle growth and development. A single base change in intron 3 of the insulin-like growth factor 2 (IGF2) gene increases meat production in pigs by disrupting a binding site for zinc finger BED domain containing protein 6 (ZBED6). Chinese indigenous pig breeds carrying the homozygous IGF2 wild-type allele produce low lean meat. We thus generate a lean pig model in Chinese Bama pig by knocking out ZBED6. In this model, we demonstrate that ZBED6 KO increases muscle and internal organ growth through ZBED6-IGF2 axis and other target genes. These results not only open new strategies for lean meat breeding in Chinese indigenous pigs, but also provide new insights to the global function of ZBED6 in organ growth and development.

ZBED6 counteracts high-fat diet-induced glucose intolerance by maintaining beta cell area and reducing excess mitochondrial activation

AIMS/HYPOTHESIS

ZBED6 (zinc finger, BED-type containing 6) is known to regulate muscle mass by suppression of Igf2 gene transcription. In insulin-producing cell lines, ZBED6 maintains proliferative capacity at the expense of differentiation and beta cell function. The aim was to study the impact of Zbed6 knockout on beta cell function and glucose tolerance in C57BL/6 mice.

METHODS

Beta cell area and proliferation were determined in Zbed6 knockout mice using immunohistochemical analysis. Muscle and fat distribution were assessed using micro-computed tomography. Islet gene expression was assessed by RNA sequencing. Effects of a high-fat diet were analysed by glucose tolerance and insulin tolerance tests. ZBED6 was overexpressed in EndoC-βH1 cells and human islet cells using an adenoviral vector. Beta cell cell-cycle analysis, insulin release and mitochondrial function were studied in vitro using propidium iodide staining and flow cytometry, ELISA, the Seahorse technique, and the fluorescent probes JC-1 and MitoSox.

RESULTS

Islets from Zbed6 knockout mice showed lowered expression of the cell cycle gene Pttg1, decreased beta cell proliferation and decreased beta cell area, which occurred independently from ZBED6 effects on Igf2 gene expression. Zbed6 knockout mice, but not wild-type mice, developed glucose intolerance when given a high-fat diet. The high-fat diet Zbed6 knockout islets displayed upregulated expression of oxidative phosphorylation genes and genes associated with beta cell differentiation. In vitro, ZBED6 overexpression resulted in increased EndoC-βH1 cell proliferation and a reduced glucose-stimulated insulin release in human islets. ZBED6 also reduced mitochondrial JC-1 J-aggregate formation, mitochondrial oxygen consumption rates (OCR) and mitochondrial reactive oxygen species (ROS) production, both at basal and palmitate + high glucose-stimulated conditions. ZBED6-induced inhibition of OCR was not rescued by IGF2 addition. ZBED6 reduced levels of the mitochondrial regulator PPAR-γ related coactivator 1 protein (PRC) and bound its promoter/enhancer region. Knockdown of PRC resulted in a lowered OCR.

CONCLUSIONS/INTERPRETATION

It is concluded that ZBED6 is required for normal beta cell replication and also limits excessive beta cell mitochondrial activation in response to an increased functional demand. ZBED6 may act, at least in part, by restricting PRC-mediated mitochondrial activation/ROS production, which may lead to protection against beta cell dysfunction and glucose intolerance in vivo.

ZBED6 regulates Igf2 expression partially through its regulation of miR483 expression

The expression of Igf2 in mammals shows a complex regulation involving multiple promoters and epigenetic mechanisms. We previously identified a novel regulatory mechanism based on the interaction between the transcriptional factor ZBED6 and Igf2 intron. Disruption of the ZBED6-Igf2 interaction leads to a dramatic up-regulation of IGF2 expression postnatally. In the current study we characterize an additional layer of regulation involving miR483 encoded by another Igf2 intron. We found a highly significant up-regulation of miR483 expression when the ZBED6-Igf2 axis is disrupted in transgenic mice. Furthermore, CRISPR/Cas9 mediated knock-out of miR483 in C2C12 myoblast cells, both wild-type and cells with disrupted ZBED6-Igf2 axis (Igf2dGGCT), resulted in down-regulation of Igf2 expression and a reduced proliferation rate. This was further validated using miR483 mimics and inhibitors. RNA-seq analysis revealed a significant enrichment of genes involved in the PI3K-Akt signaling pathway among genes down-regulated in miR483-/- cells, including Igf2 down-regulation. The opposite pattern was observed in Igf2dGGCT cells, where Igf2 is up-regulated. Our data suggest a positive feedback between miR483 and Igf2 promoter activity, strongly affecting how ZBED6 controls Igf2 expression in various cell types.

Evidence for and localization of proposed causative variants in cattle and pig genomes

BACKGROUND

This paper reviews the localization of published potential causative variants in contemporary pig and cattle reference genomes, and the evidence for their causality. In spite of the difficulties inherent to the identification of causative variants from genetic mapping and genome-wide association studies, researchers in animal genetics have proposed putative causative variants for several traits relevant to livestock breeding.

RESULTS

For this review, we read the literature that supports potential causative variants in 13 genes (ABCG2, DGAT1, GHR, IGF2, MC4R, MSTN, NR6A1, PHGK1, PRKAG3, PLRL, RYR1, SYNGR2 and VRTN) in cattle and pigs, and localized them in contemporary reference genomes. We review the evidence for their causality, by aiming to separate the evidence for the locus, the proposed causative gene and the proposed causative variant, and report the bioinformatic searches and tactics needed to localize the sequence variants in the cattle or pig genome.

CONCLUSIONS

Taken together, there is usually good evidence for the association at the locus level, some evidence for a specific causative gene at eight of the loci, and some experimental evidence for a specific causative variant at six of the loci. We recommend that researchers who report new potential causative variants use referenced coordinate systems, show local sequence context, and submit variants to repositories.

Comprehensive analysis of lncRNAs involved in skeletal muscle development in ZBED6-knockout Bama pigs

BACKGROUND

The mutation of insulin-like growth factor 2 (IGF2 mutation) that a single-nucleotide substitution (G→A) in the third intron of IGF2 abrogates the interaction with zinc finger BED-type containing 6 (ZBED6) and leads to increased muscle mass in pigs. IGF2 mutation knock-in (IGF2 KI) and ZBED6 knockout (ZBED6 KO) lead to changes in IGF2 expression and increase muscle mass in mice and pigs. Long noncoding RNAs (lncRNAs) may participate in numerous biological processes, including skeletal muscle development. However, the role of the ZBED6-lncRNA axis in skeletal muscle development is poorly characterized.

RESULTS

In this study, we assembled transcriptomes using RNA-seq data published in previous studies by our group and identified 11,408 known lncRNAs and 2269 potential lncRNAs in seven tissues, heart, longissimus dorsi, gastrocnemius muscle, liver, spleen, lung and kidney, of ZBED6 KO (lean mass model) and WT Bama pigs. ZBED6 affected the expression of 1570 lncRNAs (differentially expressed lncRNAs [DE-lncRNAs]; log2-fold change ≥ 1, nominal p-value ≤ 0.05) in the seven examined tissues. The expressed lncRNAs (FPKM > 0.1) exhibited tissue-specific patterns in WT pigs. Specifically, 3410 lncRNAs were expressed exclusively in only one tissue. Potential functions of lncRNAs were indirectly predicted by searching their target cis- and trans-regulated protein-coding genes. LncRNAs with tissue-specific expression influence numerous genes related to tissue functions. Weighted gene coexpression network analysis (WGCNA) of 1570 DE-lncRNAs between WT and ZBED6 KO pigs was used to define the following six lncRNA modules specific to different tissues: skeletal muscle, heart, lung, spleen, kidney and liver modules. Furthermore, by conjoint analysis of longissimus dorsi data (tissue-specific expression, muscle module and DE-lncRNAs) and ChIP-PCR revealed NONSUSG002145.1 (adjusted p-values = 0.044), which is coexpressed with the IGF2 gene and binding with ZBED6, may play important roles in ZBED6 KO pig skeletal muscle development.

CONCLUSIONS

These findings indicate that the identified lncRNAs may play essential roles in tissue function and regulate the mechanism of ZBED6 action in skeletal muscle development in pigs. To our knowledge, this is the first study describing lncRNAs in ZBED6 KO pigs. These results may open new research directions leading to a better understanding of the global functions of ZBED6 and of lncRNA functions in skeletal muscle development in pigs.

A High-Resolution Melting Analysis with an Unlabeled Probe for CRISPR/Cas9-Induced ZBED6 Knockout Pigs Detection

BACKGROUND

The zinc finger BED-type containing six knockout (ZBED6-KO) pigs were created to improve economic traits by increasing the expression of insulin-like growth factor 2. They were generated by CRISPR/CRISPR-associated protein 9 (Cas9) technology and a single-base deletion of ZBED6 was found. An efficient and rapid method was needed to detect this type of pig.

OBJECTIVE

This study aimed to develop a high-resolution melting (HRM) method to detect ZBED6-KO pigs.

METHODS

An unlabeled probe and two primers were designed to develop the HRM method. The limit of detection, specificity, and accuracy of the established method were tested by the constructed plasmid and DNA extracts of tissue specimens.

RESULTS

The limit of detection by the established method was 102 copies/µL. The HRM method with an unlabeled probe showed good specificity and high accuracy.

CONCLUSIONS

The established HRM analysis with an unlabeled probe showed it to be a highly effective, rapid, and reliable method to distinguish ZBED6-KO pigs from wild-type pigs.

HIGHLIGHTS

It is the first time that HRM analysis with an unlabeled probe has been used in the detection of genome editing pigs by the CRISPR/Cas9 technology.

Transposable element-derived sequences in vertebrate development

Transposable elements (TEs) are major components of all vertebrate genomes that can cause deleterious insertions and genomic instability. However, depending on the specific genomic context of their insertion site, TE sequences can sometimes get positively selected, leading to what are called "exaptation" events. TE sequence exaptation constitutes an important source of novelties for gene, genome and organism evolution, giving rise to new regulatory sequences, protein-coding exons/genes and non-coding RNAs, which can play various roles beneficial to the host. In this review, we focus on the development of vertebrates, which present many derived traits such as bones, adaptive immunity and a complex brain. We illustrate how TE-derived sequences have given rise to developmental innovations in vertebrates and how they thereby contributed to the evolutionary success of this lineage.

Liver Expression of IGF2 and Related Proteins in ZBED6 Gene-Edited Pig by RNA-Seq

Simple Summary

Zinc finger BED-type containing 6 (ZBED6), as a regulatory factor, has different regulatory mechanisms in animal development. The intron of insulin-like growth factor 2 (IGF2) regulates the development of animal muscle and adipose by combining with the binding site of ZBED6. As a member of the insulin-like growth factor family, IGF2 plays an important role in embryonic growth and development, cell proliferation, muscle growth and genome imprinting. In order to further study the regulatory mechanism of ZBED6 on IGF2, we detected the expression of IGF2 and related genes in ZBED6 single allele knockout (ZBED6-SKO) pig tissues and analyzed differently expressed genes of the transcriptome of ZBED6-SKO pig liver. The results showed that the partial knockout of ZBED6 could affect the secretion of IGF2 in pig liver but had no significant difference at the protein level. This research provides a new idea for the interaction between IGF2 and ZBED6.

Abstract

Zinc finger BED-type containing 6 (ZBED6), a highly conservative transcription factor of placental mammals, has conservative interaction of insulin-like growth factor 2 (IGF2) based on the 16 bp binding sites of ZBED6 on the IGF2 sequence. IGF2 is related to embryo growth and cell proliferation. At the same time, its functions in muscle and adipose in mammals have been widely mentioned in recent studies. To further investigate the mechanism of ZBED6 on IGF2, we detected the expression of IGF2 and related genes in ZBED6 single allele knockout (ZBED6-SKO) pig tissues and analyzed the transcriptome of ZBED6-SKO pig liver. Through RNA-seq, we captured nine up-regulated genes and eight down-regulated genes which related to lipid metabolism. The results showed that the mRNA of IGF2 had an upward trend after the partial knockout of ZBED6 in liver and had no significant difference in protein expression of IGF2. In summary, ZBED6-SKO could affect the secretion of IGF2 in pig liver and its own lipid metabolism. Our research has provided basic information for revealing the regulatory mechanism of the interaction between ZBED6 and IGF2 in mammals.

ZBED6 binding motifs correlate with endogenous retroviruses and Syncytin genes

Retroviruses have infiltrated vertebrate germlines for millions of years as inherited endogenous retroviruses (ERVs). Mammalian genomes host large numbers of ERVs and transposable elements (TEs), including retrotransposons and DNA transposons, that contribute to genomic innovation and evolution as coopted genes and regulators of diverse functions. To explore features distinguishing coopted ERVs and TEs from other integrations, we focus on the potential role of ZBED6 and repeated ERV domestication as repurposed Syncytin genes. The placental mammal-specific ZBED6 is a DNA transposon-derived transcription regulator and we demonstrate that its binding motifs are associated with distinct Syncytins and that ZBED6 binding motifs are 2- to 3-fold more frequent in ERVs than in flanking DNA. Our observations suggest that ZBED6 could contribute an extended regulatory role of genomic expression, utilizing ERVs as platforms for genomic innovation and evolution.

Goat sperm associated antigen 17 protein gene (SPAG17): Small and large fragment genetic variation detection, association analysis, and mRNA expression in gonads

Transcriptome sequencing analyses have suggested that sperm associated antigen 17 protein gene (SPAG17) may play important regulating roles in litter size. In this study, the expression profiles and genetic variations of the SPAG17 were studied in Shaanbei White Cashmere (SBWC) goats (n=1567). SPAG17 was highly expressed in testis and ovary of SBWC goats. At different developmental stages, it also continued to be highly expressed in testis. In addition, two variations of SPAG17, one indel locus and one copy number variation locus, were significantly associated with first-born litter size. Joint analysis results suggested that two polymorphic loci of the SPAG17 gene may regulate host gene expression in goat ovary and testis. Overall, the results indicated the important role of SPAG17 in the reproductive process of goats.

ZBED1/DREF: A transcription factor that regulates cell proliferation

Maintenance of genomic diversity is critically dependent on gene regulation at the transcriptional level. This occurs via the interaction of regulatory DNA sequence motifs with DNA-binding transcription factors. The zinc finger, BED-type (ZBED) gene family contains major DNA-binding motifs present in human transcriptional factors. It encodes proteins that present markedly diverse regulatory functions. ZBED1 has similar structural and functional properties to its Drosophila homolog DNA replication-related element-binding factor (DREF) and plays a critical role in the regulation of transcription. ZBED1 regulates the expression of several genes associated with cell proliferation, including cell cycle regulation, chromatin remodeling and protein metabolism, and some genes associated with apoptosis and differentiation. In the present review, the origin, structure and functional role of ZBED1 were comprehensively assessed. In addition, the similarities and differences between ZBED1 and its Drosophila homolog DREF were highlighted, and future research directions, particularly in the area of clinical cancer, were discussed.

The importance of the ZBED6-IGF2 axis for metabolic regulation in mouse myoblast cells

The transcription factor ZBED6 acts as a repressor of Igf2 and affects directly or indirectly the transcriptional regulation of thousands of genes. Here, we use gene editing in mouse C2C12 myoblasts and show that ZBED6 regulates Igf2 exclusively through its binding site 5'-GGCTCG-3' in intron 1 of Igf2. Deletion of this motif (Igf2^ΔGGCT ) or complete ablation of Zbed6 leads to ~20-fold upregulation of the IGF2 protein. Quantitative proteomics revealed an activation of Ras signaling pathway in both Zbed6^-/- and Igf2^ΔGGCT myoblasts, and a significant enrichment of mitochondrial membrane proteins among proteins showing altered expression in Zbed6^-/- myoblasts. Both Zbed6^-/- and Igf2^ΔGGCT myoblasts showed a faster growth rate and developed myotube hypertrophy. These cells exhibited an increased O₂ consumption rate, due to IGF2 upregulation. Transcriptome analysis revealed ~30% overlap between differentially expressed genes in Zbed6^-/- and Igf2^ΔGGCT myotubes, with an enrichment of upregulated genes involved in muscle development. In contrast, ZBED6-overexpression in myoblasts led to cell apoptosis, cell cycle arrest, reduced mitochondrial activities, and ceased myoblast differentiation. The similarities in growth and differentiation phenotypes observed in Zbed6^-/- and Igf2^ΔGGCT myoblasts demonstrates that ZBED6 affects mitochondrial activity and myogenesis largely through its regulation of IGF2 expression. This study adds new insights how the ZBED6-Igf2 axis affects muscle metabolism.

Identification of eQTLs associated with lipid metabolism in Longissimus dorsi muscle of pigs with different genetic backgrounds

Intramuscular fat content and its fatty acid composition affect porcine meat quality and its nutritional value. The present work aimed to identify genomic variants regulating the expression in the porcine muscle (Longissimus dorsi) of 45 candidate genes for lipid metabolism and fatty acid composition in three experimental backcrosses based on the Iberian breed. Expression genome-wide association studies (eGWAS) were performed between the muscle gene expression values, measured by real-time quantitative PCR, and the genotypes of 38,426 SNPs distributed along all chromosomes. The eGWAS identified 186 eSNPs located in ten Sus scrofa regions and associated with the expression of ACSM5, ACSS2, ATF3, DGAT2, FOS and IGF2 (FDR < 0.05) genes. Two expression quantitative trait loci (eQTLs) for IGF2 and ACSM5 were classified as cis-acting eQTLs, suggesting a mutation in the same gene affecting its expression. Conversely, ten eQTLs showed trans-regulatory effects on gene expression. When the eGWAS was performed for each backcross independently, only three common trans-eQTL regions were observed, indicating different regulatory mechanisms or allelic frequencies among the breeds. In addition, hotspot regions regulating the expression of several genes were detected. Our results provide new data to better understand the functional regulatory mechanisms of lipid metabolism genes in muscle.

ZBED2 is an antagonist of interferon regulatory factor 1 and modifies cell identity in pancreatic cancer

Lineage plasticity is a prominent feature of pancreatic ductal adenocarcinoma (PDA) cells, which can occur via deregulation of lineage-specifying transcription factors. Here, we show that the zinc finger protein ZBED2 is aberrantly expressed in PDA and alters tumor cell identity in this disease. Unexpectedly, our epigenomic experiments reveal that ZBED2 is a sequence-specific transcriptional repressor of IFN-stimulated genes, which occurs through antagonism of IFN regulatory factor 1 (IRF1)-mediated transcriptional activation at cooccupied promoter elements. Consequently, ZBED2 attenuates the transcriptional output and growth arrest phenotypes downstream of IFN signaling in multiple PDA cell line models. We also found that ZBED2 is preferentially expressed in the squamous molecular subtype of human PDA, in association with inferior patient survival outcomes. Consistent with this observation, we show that ZBED2 can repress the pancreatic progenitor transcriptional program, enhance motility, and promote invasion in PDA cells. Collectively, our findings suggest that high ZBED2 expression is acquired during PDA progression to suppress the IFN response pathway and to promote lineage plasticity in this disease.

An 11-bp Indel Polymorphism within the CSN1S1 Gene Is Associated with Milk Performance and Body Measurement Traits in Chinese Goats

The casein alpha s1 (CSN1S1) gene encodes α-s1 casein, one of the proteins constituting milk, which affects milk performance, as well as improving the absorption of calcium and bone development in mammals. A previous study found that an 11-bp insertion/deletion (indel) of this gene strongly affected litter size in goats. However, to our knowledge, the relationships between this polymorphism and the milk performance and body measurement traits of goats have not been reported. In this paper, the previously identified indel has been recognized in three Chinese goat breeds, namely the Guanzhong dairy goat (GZDG; n = 235), Shaanbei white cashmere goat (SBWC; n = 1092), and Hainan black goat (HNBG; n = 278), and the following three genotypes have been studied for all of the breeds: insertion/insertion (II), deletion/deletion (DD), and insertion/deletion (ID). The allele frequencies analyzed signified that the frequencies of the "D" allele were higher (47.8%-65.5%), similar to the previous report, which indicates that this polymorphism is genetically stable in different goat breeds. Further analysis showed that this indel was markedly associated with milk fat content, total solids content, solids-not-fat content, freezing point depression, and acidity in GZDG (p < 0.05), and also affected different body measurement traits in all three breeds (p < 0.05). The goats with II genotypes had superior milk performance, compared with the others; however, goats with DD genotypes had better body measurement sizes. Hence, it may be necessary to select goats with an II or DD genotype, based on the desired traits, while breeding. Our study provides information on the potential impact of the 11-bp indel polymorphism of the CSN1S1 gene for improving the milk performance and body measurement traits in goats.

Pig Hsd17b3: Alternative splice variants expression, insertion/deletion (indel) in promoter region and their associations with male reproductive traits

Hydroxysteroid 17-Beta Dehydrogenase 3 (Hsd17b3), primarily expressed in Leydig cells (LCs) of the mammalian testes, is essential for testosterone biosynthesis and male fertility. The aim of our study was to profile the expression, splice variants (SV) and novel insertion/deletion (indel) of Hsd17b3 in boars. Quantitative analysis showed that the expression level of Hsd17b3 in the testis was significantly highest. Among different testicular cell types, the Hsd17b3 mRNA expression level of LCs was significantly higher than that of SSCs (spermatogonial stem cells) and SCs (Sertoli cells). Furthermore, the SV was firstly identified in pigs and it was highly expressed in LCs comparing with SSCs and SCs. In addition, two mutations were identified in pig Hsd17b3 gene promotor and intron, respectively, which were associated with male reproductive traits (P <  0.05). In conclusion, both transcripts of Hsd17b3 gene were highly expressed in pig testes and LCs; the two novel indel variants of Hsd17b3 gene can be used as potential DNA makers for the marker-assisted selection in pigs. All these findings would enrich the study of Hsd17b3 gene in pig genetic breeding.

A novel mutation in the promoter region of RPL8 regulates milk fat traits in dairy cattle by binding transcription factor Pax6

Ribosomal protein L8 (RPL8) was considered as a promising candidate gene for the milk fat percentage trait in dairy cattle in our previous genome-wide association studies, but the mechanism remains to be determined. Here we investigated the molecular mechanism underlying the effect of bovine RPL8 on milk fat percentage. We demonstrated that RPL8 silencing in bovine mammary epithelial cells affected the expression of genes encoding fat-related enzymes (ACACA, FASN, ACSS1, FABP3, SREBP-1, DGAT1, GPAM, PLIN2, PLIN5 and CIDEA). Furthermore, we showed here that a single nucleotide polymorphism, g.-931G > T (chr14:1508300, UMD3.1) in the putative RPL8 promoter region significantly reduced its promoter activity. Interestingly, this decrease in activity was paralleled by lower RPL8 expression in mammary gland tissues of dairy cattle with the homozygous TT genotype compared to that of cattle with the wild-type homozygous GG genotype. Importantly, we found g.-931G > T added a paired box 6 (Pax6)-binding site and this mutation located in the presumed Pax6-binding site. EMSA and co-immunoprecipitation (Co-IP) assays confirmed the interaction between RPL8 and Pax6 and the T allele exhibited a higher affinity of DNA/protein interactions than G allele, suggesting that Pax6 is an important transcription factor for RPL8 expression. In addition, lactating cows with the GG and GT genotypes presented a significant decrease in milk fat percentage compared to cows with TT genotypes. Altogether, our study indicated that g.-931G > T at RPL8 promoter altered its expression by affecting the interplay between Pax6 and RPL8, which may account for the association with milk fat traits. Findings herein first elucidated the biological function of RPL8 gene in milk fat and the identified SNP g.-931G > T may be considered as genetic makers for breeding in dairy cattle.

Analysis of porcine IGF2 gene expression in adipose tissue and its effect on fatty acid composition

IGF2:g.3072G>A polymorphism has been described as the causal mutation of a maternally imprinted QTL for muscle growth and fat deposition in pigs. The objective of the current work was to study the association between the IGF2:g.3072G>A polymorphism and the IGF2 gene expression and its effect on fatty acid composition in adipose tissue in different pig genetic backgrounds. A cis-eQTL region associated with the IGF2 mRNA expression in adipose tissue was identified in an eGWAS with 355 animals. The IGF2 gene was located in this genomic interval and IGF2g.3072G>A was the most significant SNP, explaining a 25% of the gene expression variance. Significant associations between IGF2:g.3072G>A polymorphism and oleic (C18:1(n-9); p-value = 4.18x10^-07), hexadecanoic (C16:1(n-9); p-value = 4.04x10^-07), linoleic (C18:2(n-6); p-value = 6.44x10^-09), α-linoleic (C18:3(n-3); p-value = 3.30x10^-06), arachidonic (C20:4(n-6); p-value = 9.82x10^-08) FAs and the MUFA/PUFA ratio (p-value = 2.51x10^-9) measured in backfat were identified. Animals carrying the A allele showed an increase in IGF2 gene expression and higher PUFA and lower MUFA content. However, in additional studies was observed that there could be other proximal genetic variants affecting FA composition in adipose tissue. Finally, no differences in the IGF2 gene expression in adipose tissue were found between heterozygous animals classified according to the IGF2:g.3072G>A allele inherited from the father (A^PG^M or A^MG^P). However, pyrosequencing analysis revealed that there is imprinting of the IGF2 gene in muscle and adipose tissues, with stronger differences among the paternally and maternally inherited alleles in muscle. Our results suggested that IGF2:g.3072G>A polymorphism plays an important role in the regulation of IGF2 gene expression and can be involved in the fatty acid composition in adipose tissue. In both cases, further studies are still needed to deepen the mechanism of regulation of IGF2 gene expression in adipose tissue and the IGF2 role in FA composition.

The insulin-like growth factor 2 gene in mammals: Organizational complexity within a conserved locus

The secreted protein, insulin-like growth factor 2 (IGF2), plays a central role in fetal and prenatal growth and development, and is regulated at the genetic level by parental imprinting, being expressed predominantly from the paternally derived chromosome in mice and humans. Here, IGF2/Igf2 and its locus has been examined in 19 mammals from 13 orders spanning ~166 million years of evolutionary development. By using human or mouse DNA segments as queries in genome analyses, and by assessing gene expression using RNA-sequencing libraries, more complexity was identified within IGF2/Igf2 than was annotated previously. Multiple potential 5’ non-coding exons were mapped in most mammals and are presumably linked to distinct IGF2/Igf2 promoters, as shown for several species by interrogating RNA-sequencing libraries. DNA similarity was highest in IGF2/Igf2 coding exons; yet, even though the mature IGF2 protein was conserved, versions of 67 or 70 residues are produced secondary to species-specific maintenance of alternative RNA splicing at a variable intron-exon junction. Adjacent H19 was more divergent than IGF2/Igf2, as expected in a gene for a noncoding RNA, and was identified in only 10/19 species. These results show that common features, including those defining IGF2/Igf2 coding and several non-coding exons, were likely present at the onset of the mammalian radiation, but that others, such as a putative imprinting control region 5’ to H19 and potential enhancer elements 3’ to H19, diversified with speciation. This study also demonstrates that careful analysis of genomic and gene expression repositories can provide new insights into gene structure and regulation.

sgRNA-shRNA Structure Mediated SNP Site Editing on Porcine IGF2 Gene by CRISPR/StCas9

The SNP within intron 3 of the porcine IGF2 gene (G3072A) plays an important role for muscle growth and fat deposition in pigs. In this study, the StCas9 derived from Streptococcus thermophilus together with the Drosha-mediated sgRNA-shRNA structure were combined to boost the G to A base editing on the IGF2 SNP site, which we called “SNP editing.” The codon-humanized StCas9 as we previously reported was firstly compared with the prevalently used SpCas9 derived from Streptococcus pyogenes using our idiomatic surrogate report assay, and the StCas9 demonstrated a comparable targeting activity. On the other hand, by combining shRNA with sgRNA, simultaneous gene silencing and genome targeting can be achieved. Thus, the novel IGF2.sgRNA-LIG4.shRNA-IGF2.sgRNA structure was constructed to enhance the sgRNA/Cas9-mediated HDR-based IGF2 SNP editing by silencing the LIG4 gene, which is a key molecule of the HDR’s competitive NHEJ pathway. The sgRNA-shRNA/StCas9 all-in-one expression vector and the IGF2.sgRNA/StCas9 as control were separately used to transfect porcine PK15 cells together with an ssODNs donor for the IGF2 SNP editing. The editing events were detected by the RFLP assay, Sanger sequencing as well as Deep-sequencing, and the Deep-sequencing results finally demonstrated a significant higher HDR-based editing efficiency (16.38%) for our sgRNA-shRNA/StCas9 strategy. In short, we achieved effective IGF2 SNP editing by using the combined sgRNA-shRNA/StCas9 strategy, which will facilitate the further production of base-edited animals and perhaps extend for the gene therapy for the base correction of some genetic diseases.

IMP2 Increases Mouse Skeletal Muscle Mass and Voluntary Activity by Enhancing Autocrine Insulin-Like Growth Factor 2 Production and Optimizing Muscle Metabolism

Insulin-like growth factor 2 (IGF2) mRNA binding protein 2 (IMP2) was selectively deleted from adult mouse muscle; two phenotypes were observed: decreased accrual of skeletal muscle mass after weaning and reduced wheel-running activity but normal forced treadmill performance. Reduced wheel running occurs when mice are fed a high-fat diet but is normalized when mice consume standard chow. The two phenotypes are due to altered output from different IMP2 client mRNAs. The reduced fiber size of IMP2-deficient muscle is attributable, in part, to diminished autocrine Igf2 production; basal tyrosine phosphorylation of the insulin and IGF1 receptors is diminished, and Akt1 activation is selectively reduced. Gsk3α is disinhibited, and S536-phosphorylated ε subunit of eukaryotic initiation factor 2B [eIF2Bε(S536)] is hyperphosphorylated. Protein synthesis is reduced despite unaltered mTOR complex 1 activity. The diet-dependent reduction in voluntary exercise is likely due to altered muscle metabolism, as contractile function is normal. IMP2-deficient muscle exhibits reduced fatty acid oxidation, due to a reduced abundance of mRNA of peroxisome proliferator-activated receptor α (PPARα), an IMP2 client, and PPARα protein. IMP2-deficient muscle fibers treated with a mitochondrial uncoupler to increase electron flux, as occurs with exercise, exhibit reduced oxygen consumption from fatty acids, with higher oxygen consumption from glucose. The greater dependence on muscle glucose metabolism during increased oxygen demand may promote central fatigue and thereby diminish voluntary activity.

Comparative transcriptome analysis reveals regulators mediating breast muscle growth and development in three chicken breeds

Objective: The goal of this study was to investigate the mechanisms of muscle growth and development of three chicken breeds. Participants: Eighteen chickens, including three different breeds with different growth speeds (White Broiler, Daheng, and Commercial Layers of Roman), were used. Methods: Total RNA from breast muscle of these chickens was subjected to a gene expression microarray. Differentially expressed genes (DEGs) were screened and functional enrichment analysis was performed using DAVID. Seven DEGs were confirmed by quantitative reverse transcription PCR. Results: Overall, 8,398 DEGs were found among the different lines. The DEGs between each two lines that were unique for a developmental stage were greater than those that were common during all stages. Functional analysis revealed that DEGs across the entire developmental process were primarily involved in positive cell proliferation, growth, cell differentiation, and developmental processes. Genes involved in muscle regulation, muscle construction, and muscle cell differentiation were upregulated in the faster-growing breed compared to the slower-growing breed. DEGs including myosin heavy chain 15 (MYH15), myozenin 2 (MYOZ2), myosin-binding protein C (MYBPC3), insulin-like growth factor 2 (IGF2), apoptosis regulator (BCL-2), AP-1 transcription factor subunit (JUN), and AP-1 transcription factor subunit (FOS) directly regulated muscle growth or were in the center of the protein-protein interaction network. Pathways, including the extracellular matrix (ECM)-receptor interaction, mitogen-activated protein kinase (MAPK) signaling pathway, and focal adhesion, were the most enriched DEGs between lines or within lines under different developmental stages. Conclusions: Genes involved in muscle construction and cell differentiation were differentially expressed among the three breeds.

Editing porcine IGF2 regulatory element improved meat production in Chinese Bama pigs

Insulin-like growth factor 2 (IGF2) is an important growth factor, which promotes growth and development in mammals during fetal and postnatal stages. Using CRISPR-Cas9 system, we generated multiple founder pigs containing 12 different mutant alleles around a regulatory element within the intron 3 of IGF2 gene. Crossing two male founders passed four mutant alleles onto F1 generation, and these mutations abolished repressor ZBED6 binding and rendered this regulatory element nonfunctional. Both founders and F1 animals showed significantly faster growth, without affecting meat quality. These results indicated that editing IGF2 intron 3-3072 site using CRISPR-Cas9 technology improved meat production in Bama pigs. This is the first demonstration that editing non-coding region can improve economic traits in livestock.

Disruption of the ZBED6 binding site in intron 3 of IGF2 by CRISPR/Cas9 leads to enhanced muscle development in Liang Guang Small Spotted pigs

Insulin-like growth factor 2 (IGF2) plays an important role in the development of the foetus and in post-natal growth and development. A SNP within intron 3 of porcine IGF2 disrupts a binding site for the repressor, zinc finger BED-type containing 6 (ZBED6), leading to up-regulation of IGF2 in skeletal muscle and major effects on muscle growth, heart size, and fat deposition. This favourable mutation is common in Western commercial pig populations, but is not present in most indigenous Chinese pig breeds. Here, we described the efficient disruption of the ZBED6 binding site motif in intron 3 of IGF2 by CRISPR/Cas9 in porcine embryonic fibroblasts (PEFs) from the indigenous Chinese pig breed, Liang Guang Small Spotted pig. Disruption of the binding motif led to a drastic up-regulation of IGF2 expression in PEFs and enhanced myogenic potential and cell proliferation of PEFs. IGF2-edited pigs were then generated using somatic cell nuclear transfer. Enhanced muscle development was evident in one pig with biallelic deletion of the ZBED6 binding site motif, implying that the release of ZBED6 repression has a major effect on porcine muscle development. Our study confirmed the important effect of a mutation in the ZBED6 binding site motif on IGF2 expression and myogenesis, thus providing the basis for breeding a new line of Liang Guang Small Spotted pigs with improved lean meat percentage, a trait of great commercial value to pig producers.

The insulin-like growth factor 2 gene and locus in nonmammalian vertebrates: Organizational simplicity with duplication but limited divergence in fish

The small, secreted peptide, insulin-like growth factor 2 (IGF2), is essential for fetal and prenatal growth in humans and other mammals. Human IGF2 and mouse Igf2 genes are located within a conserved linkage group and are regulated by parental imprinting, with IGF2/Igf2 being expressed from the paternally derived chromosome, and H19 from the maternal chromosome. Here, data retrieved from genomic and gene expression repositories were used to examine the Igf2 gene and locus in 8 terrestrial vertebrates, 11 ray-finned fish, and 1 lobe-finned fish representing >500 million years of evolutionary diversification. The analysis revealed that vertebrate Igf2 genes are simpler than their mammalian counterparts, having fewer exons and lacking multiple gene promoters. Igf2 genes are conserved among these species, especially in protein-coding regions, and IGF2 proteins also are conserved, although less so in fish than in terrestrial vertebrates. The Igf2 locus in terrestrial vertebrates shares additional genes with its mammalian counterparts, including tyrosine hydroxylase (Th), insulin (Ins), mitochondrial ribosomal protein L23 (Mrpl23), and troponin T3, fast skeletal type (Tnnt3), and both Th and Mrpl23 are present in the Igf2 locus in fish. Taken together, these observations support the idea that a recognizable Igf2 was present in the earliest vertebrate ancestors, but that other features developed and diversified in the gene and locus with speciation, especially in mammals. This study also highlights the need for correcting inaccuracies in genome databases to maximize our ability to accurately assess contributions of individual genes and multigene families toward evolution, physiology, and disease.