Antioxidant activity and metabolic regulation of sodium salicylate on goat sperm at low temperature

OBJECTIVE

The purpose of this study was to explore the effect of sodium salicylate (SS) on semen preservation and metabolic regulation in goats.

METHODS

Under the condition of low temperature, SS was added to goat semen diluent to detect goat sperm motility, plasma membrane, acrosome, antioxidant capacity, mitochondrial membrane potential (MMP) and metabonomics.

RESULTS

The results show that at the 8th day of low-temperature storage, the sperm motility of the 20 μM SS group was 66.64%, and the integrity rates of the plasma membrane and acrosome were both above 60%, significantly higher than those of the other groups. The activities of catalase and superoxide dismutase in the sperm of the 20 μM SS group were significantly higher than those of the control group, the contents of reactive oxygen species and malondialdehyde were significantly lower than those in the control group, the MMP was significantly higher than that in the control group, and the contents of Ca2+ and total cholesterol were significantly higher than those in the control group. Through metabonomics analysis, there were significant metabolic differences between the control group and the 20 μM SS group. Twenty of the most significant metabolic markers were screened, mainly involving five metabolic pathways, of which nicotinic acid and nicotinamide metabolic pathways were the most significant.

CONCLUSION

The results indicate that SS can effectively improve the low-temperature preservation quality of goat sperm.

Multi-omics analysis on the mechanism of the effect of Isatis leaf on the growth performance of fattening sheep

Introduction

This study evaluated the effects of Isatis Leaf (ISL) on the growth performance, gastrointestinal tissue morphology, rumen and intestinal microbiota, rumen, serum and urine metabolites, and rumen epithelial tissue transcriptome of fattening sheep.

Methods

Twelve 3.5-month-old healthy fattening sheep were randomly divided into two groups, each with 6 replicates, and fed with basal diet (CON) and basal diet supplemented with 80 g/kg ISL for 2.5 months. Gastrointestinal tract was collected for histological analysis, rumen fluid and feces were subjected to metagenomic analysis, rumen fluid, serum, and urine for metabolomics analysis, and rumen epithelial tissue for transcriptomics analysis.

Results

The results showed that in the ISL group, the average daily gain and average daily feed intake of fattening sheep were significantly lower than those of the CON group (P < 0.05), and the rumen ammonia nitrogen level was significantly higher than that of the CON group (P < 0.01). The thickness of the reticulum and abomasum muscle layer was significantly increased (P < 0.05). At the genus level, the addition of ISL modified the composition of rumen and fecal microorganisms, and the relative abundance of Methanobrevibacter and Centipeda was significantly upregulated in rumen microorganisms, The relative abundance of Butyrivibrio, Saccharofermentans, Mogibacterium, and Pirellula was significantly downregulated (P < 0.05). In fecal microorganisms, the relative abundance of Papillibacter, Pseudoflavonifractor, Butyricicoccus, Anaerovorax, and Methanocorpusculum was significantly upregulated, while the relative abundance of Roseburia, Coprococcus, Clostridium XVIII, Butyrivibrio, Parasutterella, Macellibacteroides, and Porphyromonas was significantly downregulated (P < 0.05). There were 164, 107, and 77 different metabolites in the rumen, serum, and urine between the ISL and CON groups (P < 0.05). The differential metabolic pathways mainly included thiamine metabolism, niacin and nicotinamide metabolism, vitamin B6 metabolism, taurine and taurine metabolism, beta-Alanine metabolism and riboflavin metabolism. These metabolic pathways were mainly involved in the regulation of energy metabolism and immune function in fattening sheep. Transcriptome sequencing showed that differentially expressed genes were mainly enriched in cellular physiological processes, development, and immune regulation.

Conclusion

In summary, the addition of ISL to the diet had the effect of increasing rumen ammonia nitrogen levels, regulating gastrointestinal microbiota, promoting body fat metabolism, and enhancing immunity in fattening sheep.

Roles of Y-27632 on sheep sperm metabolism

To investigate the effect of Y-27632 on low-temperature metabolism of sheep sperm, different concentrations of Y-27632 were added to sheep semen at 4 °C in this experiment to detect indicators such as sperm motility, plasma membrane, acrosome, antioxidant performance, mitochondrial membrane potential (MMP), and metabolomics. The results showed that the addition of 20 µM Y-27632 significantly increased sperm motility, plasma membrane integrity rate, acrosome integrity rate, antioxidant capacity, MMP level, significantly increased sperm adenosine triphosphate (ATP) and total cholesterol content, and significantly reduced sperm Ca2+ content. In metabolomics analysis, compared with the control group, the 20 µM Y-27632 group screened 20 differential metabolites, mainly involved in five metabolic pathways, with the most significant difference in Histidine metabolism (P = 0.001). The results confirmed that Y-27632 significantly improved the quality of sheep sperm preservation under low-temperature conditions.

Comprehensive analyses of 435 goat transcriptomes provides insight into male reproduction

A systematic analysis of genes related to reproduction is crucial for obtaining a comprehensive understanding of the molecular mechanisms that underlie male reproductive traits in mammals. Here, we utilized 435 goat transcriptome datasets to unveil the testicular tissue-specific genes (TSGs), allele-specific expression (ASE) genes and their uncharacterized transcriptional features related to male goat reproduction. Results showed a total of 1790 TSGs were identified in goat testis, which was the most among all tissues. GO enrichment analyses suggested that testicular TSGs were mainly involved in spermatogenesis, multicellular organism development, spermatid development, and flagellated sperm motility. Subsequently, a total of 95 highly conserved TSGs (HCTSGs), 508 middle conserved TSGs (MCTSGs) and 42 no conserved TSGs (NCTSGs) were identified in goat testis. GO enrichment analyses suggested that the HCTSGs and MCTSGs has a more important association with male reproduction than NCTSGs. Additionally, we identified 644 ASE genes, including 88 tissue-specific ASE (TS-ASE) genes (e.g., FSIP2, TDRD9). GO enrichment analyses indicated that both ASE genes and TS-ASE genes were associated with goat male reproduction. Overall, this study revealed an extensive gene set involved in the regulation of male goat reproduction and their dynamic transcription patterns. Data reported here provide valuable insights for a further improvement of the economic benefits of goats as well as future treatments for male infertility.

Metabolomics Analysis of Sodium Salicylate Improving the Preservation Quality of Ram Sperm

The aim of this study was to investigate the effects of sodium salicylate (SS) on the preservation and metabolic regulation of sheep sperm. Under 4 °C low-temperature conditions, SS (at 10 µM, 20 µM, 30 µM, and 50 µM) was added to the semen diluent to detect sperm motility, plasma membrane, and acrosome integrity. Based on the selected optimal concentration of SS (20 µM), the effects of 20 µM of SS on sperms' antioxidant capacity and mitochondrial membrane potential (MMP) were evaluated, and metabolomics analysis was conducted. The results showed that on the 20th day of low-temperature storage, the sperm motility of the 20 µM SS group was 62.80%, and the activities of catalase (CAT) and superoxide dismutase (SOD) were significantly higher than those of the control group (p < 0.01). The content of Ca2+, reactive oxygen species (ROS), and malondialdehyde (MDA) were significantly lower than those of the control group (p < 0.01), and the total antioxidant capacity (T-AOC) was significantly higher than that of the control group (p < 0.05); mitochondrial activity and the total cholesterol (TC) content were significantly higher than those in the control group (p < 0.01). An ultrastructural examination showed that in the SS group, the sperm plasma membrane and acrosome were intact, the fibrous sheath and axoneme morphology of the outer dense fibers were normal, and the mitochondria were arranged neatly. In the control group, there was significant swelling of the sperm plasma membrane, rupture of the acrosome, and vacuolization of mitochondria. Using metabolomics analysis, 20 of the most significant differential metabolic markers were screened, mainly involving 6 metabolic pathways, with the amino acid biosynthesis pathway being the most abundant. In summary, 20 µM of SS significantly improved the preservation quality of sheep sperm under low-temperature conditions of 4 °C.

Multi-Omics Analysis of the Mechanism of Mentha Haplocalyx Briq on the Growth and Metabolic Regulation of Fattening Sheep

Mentha haplocalyx Briq (MHB) and its components have been proven to improve the growth performance of livestock and poultry. The aim of this experiment was to investigate the effects of MHB addition on growth performance, rumen and fecal microbiota, rumen fluid, serum and urine metabolism, and transcriptomics of rumen epithelial cells in meat sheep. Twelve Hu sheep were selected for the experiment and fed with basic diet (CON) and a basal diet supplemented with 80 g/kg DM of Mentha haplocalyx Briq (MHB). The experimental period was 10 weeks with the first 2 weeks as the pre-trial period. The results showed that compared with the CON group, the average daily weight gain of meat sheep in the MHB group increased by 20.1%; the total volatile fatty acid (VFA) concentration significantly increased (p < 0.05); The thickness of the cecal mucosal layer was significantly reduced (p < 0.01), while the thickness of the colonic mucosal layer was significantly increased (p < 0.05), the length of ileal villi significantly increased (p < 0.01), the thickness of colonic mucosal layer and rectal mucosal muscle layer significantly increased (p < 0.05), and the thickness of cecal mucosal layer significantly decreased (p < 0.05); The serum antioxidant capacity has increased. At the genus level, the addition of MHB changed the composition of rumen and fecal microbiota, increased the relative abundance of Paraprevotella, Alloprevotella, Marinilabilia, Saccharibacteria_genera_incertae_sedis, Subdivision5_genera_incertae_sedis and Ornatilinea in rumen microbiota, and decreased the relative abundance of Blautia (p < 0.05). The relative abundance of Prevotella, Clostridium XlVb and Parasutterella increased in fecal microbiota, while the relative abundance of Blautia and Coprococcus decreased (p < 0.05). There were significant differences in the concentrations of 105, 163, and 54 metabolites in the rumen, serum, and urine between the MHB group and the CON group (p < 0.05). The main metabolic pathways of the differences were pyrimidine metabolism, taurine and taurine metabolism, glyceride metabolism, and pentose phosphate pathway (p < 0.05), which had a significant impact on protein synthesis and energy metabolism. The transcriptome sequencing results showed that differentially expressed genes were mainly enriched in immune regulation, energy metabolism, and protein modification. Therefore, adding MHB improved the growth performance of lambs by altering rumen and intestinal microbiota, rumen, serum and urine metabolomics, and transcriptome.

MiR-188-5p regulates the proliferation and differentiation of goat skeletal muscle satellite cells by targeting calcium/calmodulin dependent protein kinase II beta

OBJECTIVE

The aim of this study was to reveal the role and regulatory mechanism of miR-188-5p in the proliferation and differentiation of goat muscle satellite cells.

METHODS

Goat skeletal muscle satellite cells isolated in the pre-laboratory were used as the test material. First, the expression of miR-188-5p in goat muscle tissues at different developmental stages was detected by quantitative reverse transcription polymerase chain reaction (qRT-PCR). In addition, miR-188-5p was transfected into goat skeletal muscle satellite cells by constructing mimics and inhibitors of miR-188-5p, respectively. The changes of differentiation marker gene expression were detected by qPCR method.

RESULTS

It was highly expressed in adult goat latissimus dorsi and leg muscles, goat fetal skeletal muscle, and at the differentiation stage of muscle satellite cells. Overexpression and interference of miR-188-5p showed that miR-188-5p inhibited the proliferation and promoted the differentiation of goat muscle satellite cells. Target gene prediction and dual luciferase assays showed that miR-188-5p could target the 3'untranslated region of the calcium/calmodulin dependent protein kinase II beta (CAMK2B) gene and inhibit luciferase activity. Further functional studies revealed that CAMK2B promoted the proliferation and inhibited the differentiation of goat muscle satellite cells, whereas si-CAMK2B restored the function of miR-188-5p inhibitor.

CONCLUSION

These results suggest that miR-188-5p inhibits the proliferation and promotes the differentiation of goat muscle satellite cells by targeting CAMK2B. This study will provide a theoretical reference for future studies on the molecular mechanisms of skeletal muscle development in goats.

Integrative proteomic and phosphoproteomic analysis in the female goat hypothalamus to study the onset of puberty

BACKGROUND

Puberty marks the end of childhood and achieve sexual maturation and fertility. The role of hypothalamic proteins in regulating puberty onset is unclear. We performed a comprehensive differential proteomics and phosphoproteomics analysis in prepubertal and pubertal goats to determine the roles of hypothalamic proteins and phosphoproteins during the onset of puberty.

RESULTS

We used peptide and posttranslational modifications peptide quantification and statistical analyses, and identified 69 differentially expressed proteins from 5,057 proteins and 576 differentially expressed phosphopeptides from 1574 phosphorylated proteins. Combined proteomic and phosphoproteomics, 759 correlated proteins were identified, of which 5 were differentially expressed only at the protein level, and 201 were only differentially expressed at the phosphoprotein level. Pathway enrichment analyses revealed that the majority of correlated proteins were associated with glycolysis/gluconeogenesis, Fc gamma R-mediated phagocytosis, focal adhesion, GABAergic synapse, and Rap1 signaling pathway. These pathways are related to cell proliferation, neurocyte migration, and promoting the release of gonadotropin-releasing hormone in the hypothalamus. CTNNB1 occupied important locations in the protein-protein interaction network and is involved in focal adhesion.

CONCLUSION

The results demonstrate that the proteins differentially expression only at the protein level or only differentially expressed at the phosphoprotein level and their related signalling pathways are crucial in regulating puberty in goats. These differentially expressed proteins and phosphorylated proteins may constitute the proteomic backgrounds between the two different stages.

CRISPR-based m6A modification and its potential applications in telomerase regulation

Telomerase determines cell lifespan by controlling chromosome stability and cell viability, m⁶A epigenetic modification plays an important role in the regulation of telomerase activity. Using CRISPR epigenome editing to analyze specific m⁶A modification sites in telomerase will provide an important tool for analyzing the molecular mechanism of m⁶A modification regulating telomerase activity. In this review, we clarified the relevant applications of CRISPR system, paid special attention to the regulation of m⁶A modification in stem cells and cancer cells based on CRISPR system, emphasized the regulation of m⁶A modification on telomerase activity, pointed out that m⁶A modification sites regulate telomerase activity, and discussed strategies based on telomerase activity and disease treatment, which are helpful to promote the research of anti-aging and tumor related diseases.

Analysis of serum reproductive hormones and ovarian genes in pubertal female goats

BACKGROUND

Age at puberty is an important factor affecting goat fertility, with endocrine and genetic factors playing a crucial role in the onset of puberty. To better understand the relationship between endocrine and genetic factors and mechanisms underlying puberty onset in goats, reproductive hormone levels were analyzed by ELISA and ultraperformance liquid chromatography-multiple reaction monitoring-multistage/mass spectrometry and RNA sequencing was performed to analyze ovarian genes.

RESULTS

Serum follicle stimulating hormone, luteinizing hormone, estradiol, 11-deoxycortisol, 11-deoxycorticosterone, corticosterone, cortisone, and cortisol levels were found to be higher but progesterone were lower in pubertal goats as compared to those in prepubertal goats (P < 0.05). A total of 18,139 genes were identified in cDNA libraries, and 75 differentially expressed genes (DEGs) were identified (|log₂ fold change|≥ 1, P ≤ 0.05), of which 32 were significantly up- and 43 were down-regulated in pubertal goats. Gene ontology enrichment analyses indicated that DEGs were mainly involved in "metabolic process," "signaling," "reproduction," and "growth." Further, DEGs were significantly enriched in 91 Kyoto Encyclopedia of Genes and Genomes pathways, including estrogen signaling pathway, steroid hormone biosynthesis, and cAMP signaling pathway. Bioinformatics analysis showed that PRLR and THBS1 were highly expressed in pubertal ovaries, and ZP3, ZP4, and ASTL showed low expression, suggesting their involvement in follicular development and lutealization.

CONCLUSIONS

To summarize, serum hormone changes and ovarian DEGs expression were investigated in our study. Further studies are warranted to comprehensively explore the functions of DEGs in goat puberty.

Effect of GnRH immunocastration on immune function in male rats

The present study aimed to reveal the effects of immunocastration on the development of the immune system in rats. Seventy rats were randomly assigned into two groups: Control (n = 35) and immunized (n = 35). Twenty-day-old rats were immunized with gonadotropin-releasing hormone (GnRH) and booster immunization was administered every two weeks (three immunizations in total). From 20-day-old rats, we collected samples every two weeks, including five immunized rats and five control rats (seven collections in total). We collected blood samples, testicles, thymuses, and spleens. The results showed that GnRH immunization increased the GnRH antibody titers and reduced the testosterone concentration (both P < 0.05). Compared with the control group, the number of CD4+CD8- cells, CD4-CD8+ cells, and CD4+CD8+ cells increased (P < 0.05) whereas the number of CD4-CD8- cells and CD4+CD25+ cells reduced in the immunized group (P < 0.05) over time. GnRH immunization also increased the relative weights of thymus and spleen (P < 0.05), serum concentrations of interleukin (IL)-2, IL-4, IL-6, IL-10, IL-17 and Interferon-γ (IFN-γ) over time (P < 0.05), and changed the mRNA levels of IL-2, IL-4, IL-6. IL-10, IL-17, IFN-γ, CD4, D8, CD19 GnRH, and GnRH receptor (GnRH-R) in thymus and spleen. Thus, GnRH immunization enhanced the immune markers in thymus, spleen, and blood immune cytokines in rats.

Quantitative proteomics analysis to assess protein expression levels in the ovaries of pubescent goats

Background

Changes in the abundance of ovarian proteins play a key role in the regulation of reproduction. However, to date, no studies have investigated such changes in pubescent goats. Herein we applied isobaric tags for relative and absolute quantitation (iTRAQ) and liquid chromatography–tandem mass spectrometry to analyze the expression levels of ovarian proteins in pre-pubertal (n = 3) and pubertal (n = 3) goats.

Results

Overall, 7,550 proteins were recognized; 301 (176 up- and 125 downregulated) were identified as differentially abundant proteins (DAPs). Five DAPs were randomly selected for expression level validation by Western blotting; the results of Western blotting and iTRAQ analysis were consistent. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis indicated that DAPs were enriched in olfactory transduction, glutathione metabolism, and calcium signaling pathways. Besides, gene ontology functional enrichment analysis revealed that several DAPs enriched in biological processes were associated with cellular process, biological regulation, metabolic process, and response to stimulus. Protein–protein interaction network showed that proteins interacting with CDK1, HSPA1A, and UCK2 were the most abundant.

Conclusions

We identified 301 DAPs, which were enriched in olfactory transduction, glutathione metabolism, and calcium signaling pathways, suggesting the involvement of these processes in the onset of puberty. Further studies are warranted to more comprehensively explore the function of the identified DAPs and aforementioned signaling pathways to gain novel, deeper insights into the mechanisms underlying the onset of puberty.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08699-y.

Knockdown of Ptprn-2 delays the onset of puberty in female rats

In the present study, we evaluated how Ptprn-2 (encoding tyrosine phosphatase, receptor type, N2 polypeptide protein) affects the onset of puberty in female rats. We evaluated the expression of Ptprn-2 mRNA and protein in the hypothalamus-pituitary-ovary axis of female rats using real-time quantitative reverse transcription polymerase chain reaction (RT-qPCR) and immunofluorescence at infancy, prepuberty, puberty, peripuberty, and adulthood. We evaluated the effects of Ptprn-2 gene knockdown on different aspects of reproduction-related biology in female rats, including the expression levels of puberty-related genes in vivo and in vitro, the time to onset of puberty, the concentration of serum reproductive hormones, the morphology of ovaries, and the ultrastructure of pituitary gonadotropin cells. Our results demonstrated that PTPRN-2 was primarily distributed in the arcuate nucleus (ARC), periventricular nucleus (PeN), adenohypophysis, and the ovarian follicular theca, stroma, and granulosa cells of female rats at various stages. Ptprn-2 mRNA levels significantly varied between peripuberty and puberty (P < 0.05) in the hypothalamus and pituitary gland. In hypothalamic cells, Ptprn-2 knockdown decreased the expression of Ptprn-2 and Rfrp-3 mRNA (P < 0.05) and increased the levels of Gnrh and Kiss-1 mRNA (P < 0.05). Ptprn-2 knockdown in the hypothalamus resulted in delayed vaginal opening compared to the control group (n = 12, P < 0.01), and Ptprn-2, Gnrh, and Kiss-1 mRNA levels (P < 0.05) all decreased, while the expression of Igf-1 (P < 0.05) and Rfrp-3 mRNA (P < 0.01) increased. The concentrations of FSH and P₄ in the serum of Ptprn-2 knockdown rats were lower than in control animals (P < 0.05). Large transverse perimeters and longitudinal perimeters (P < 0.05) were found in the ovaries of Ptprn-2 knockdown rats. There were fewer large secretory particles from gonadotropin cells in adenohypophysis tissue of the Ptprn-2 knockdown group compared to the control group. This indicates that Ptprn-2 knockdown can regulate levels of Gnrh, Kiss-1, and Rfrp-3 mRNA in the hypothalamus, regulate the concentration of serum FSH and P₄, and alter the morphology of ovarian and gonadotropin cells, delaying the onset of puberty in female rats.

Dynamic changes of miRNAs in skeletal muscle development at New Zealand rabbits

BACKGROUND

miRNA is one of the crucial roles in the complex and dynamic network that regulates the development of skeletal muscle. The landscape of skeletal muscle miRNAs from fetus to adult in New Zealand rabbits has not been revealed yet.

RESULTS

In this study, nine RNA-seq libraries of fetus, child and adult rabbits' leg muscles were constructed. A total of 278 differentially expressed miRNAs (DEmiRNAs) were identified. In the fetus vs. child group, the main functional enrichments were involved in membrane and transport. Pathway enriched terms of up-regulated DEmiRNAs were connected with the differentiation and hypertrophy of skeletal muscle, and down-regulated ones were related to muscle structure and metabolic capacity. In the child vs. adult group, functions were associated to positioning and transportation, and pathways were relevant to ECM, muscle structure and hypertrophy. Finally, ocu-miR-185-3p and ocu-miR-370-3p, which had the most target genes, were identified as hub-miRNAs in these two groups.

CONCLUSIONS

In short, we summarized the highly expressed and uniquely expressed DEmiRNAs of fetus, child and adult rabbits' leg muscles. Besides, the potential functional changes of miRNAs in two consecutive stages have been explored. Among them, the ocu-miR-185-3p and ocu-miR-370-3p with the most target genes were selected as hub-miRNAs. These data improved the understanding of the regulatory molecules of meat rabbit development, and provided a novel perspective for molecular breeding of meat rabbits.

Morphological changes and functional circRNAs screening of rabbit skeletal muscle development

BACKGROUND

The temporal expression pattern of circular RNAs (circRNAs) across developmental stages is essential for skeletal muscle growth and functional analysis. However, there are few analyses on the potential functions of circRNAs in rabbit skeletal muscle development.

RESULTS

Initially, the paraffin sections showed extremely significant differences in the diameter, number, area and density of skeletal muscle fibers of the fetus, child, adult rabbit hind legs (P < 0.01). Then, RNA-seq libraries of these three stages were constructed. A total of 481 differentially expressed circRNAs (DE-circRNAs) and 5,658 differentially expressed genes (DEGs) were identified. Subsequently, DE-circRNAs, whose host genes were DEGs or non-DEGs, were analyzed by GO respectively. In the fetus vs. child group, up-regulated DE-circRNAs (whose host genes were DEGs) were related to muscle fiber structure, and down-regulated ones were related to mitosis. The up-regulated DE-circRNAs (whose host genes were non-DEGs) were involved in enzyme activity, methylation and glycosylation, and the down-regulated ones were involved in mitosis and catabolism. In the fetus vs. adult group, the up-regulated DE-circRNAs (whose host genes were DEGs) were related to skeletal muscle basic structure, and the down-regulated ones were also associated with cell proliferation. But the up-regulated DE-circRNAs (whose host genes were non-DEGs) were connected with regulation of histone ubiquitination, chromatin and organelles. The down-regulated DE-circRNAs were connected with the catabolism processes. In addition, novel_circ_0022663 and novel_circ_0005489, which might have coding potential, and novel_circ_0004210 and novel_circ_0001669, which might have miRNA sponge capability, were screened out.

CONCLUSIONS

In this study, hind leg muscles of fetus, child and adult rabbits were collected for paraffin section and RNA-seq to observe the structural changes of skeletal muscle and obtain circRNA expression profiles at different stages. These data provided a catalog of circRNAs related to muscle development in New Zealand rabbits, allowing us to better understand the functional transitions in mammalian muscle development.

miR-130a-3p regulates steroid hormone synthesis in goat ovarian granulosa cells by targeting the PMEPA1 gene

MicroRNAs (miRNAs) are key epigenomic regulators of proliferation, differentiation, and secretion in cells involved in follicular development. We here studied the functional role of one such molecule, miR-130a-3p, in goat ovarian granulosa cells (GCs). High expression of this miRNA was evident in goat GCs by fluorescence in situ hybridization and suppressed estradiol and progesterone secretion from these cells, as determined by ELISA. miR-130a-3p was predicted to have a binding site for the 3' UTR of the prostate transmembrane protein androgen induced 1 gene (PMEPA1), and this was verified by a dual-luciferase reporter assay. PMEPA1 mRNA and protein expression were both found to be regulated by miR-130a-3p in GCs. Moreover, the overexpression or knockdown of PMEPA1 enhanced or suppressed estradiol and progesterone secretion from these cells, respectively. Furthermore, the secretion of estradiol and progesterone did not change significantly after the offsetting of PMEPA1 overexpression in GCs by miR-130a-3p. In summary, our present data indicate that miR-130a-3p inhibits the secretion of estradiol and progesterone in GCs by targeting PMEPA1. Our study thus provides seminal data and important new insights into the regulation of reproductive mechanisms in the nanny goat and other female mammals.

Transcription profiles of oocytes during maturation and embryos during preimplantation development in vivo in the goat

RNA sequencing performed on goat matured oocytes and preimplantation embryos generated invivo enabled us to define the transcriptome for goat preimplantation embryo development. The largest proportion of changes in gene expression in goat was found at the 16-cell stage, not as previously defined at the 8-cell stage, and is later than in other mammalian species. In all, 6482 genes were identified to be significantly differentially expressed across all consecutive developmental stage comparisons, and the important signalling pathways involved in each development transition were determined. In addition, we identified genes that appear to be transcribed only at a specific stage of development. Using weighted gene coexpression network analysis, we found nine stage-specific modules of coexpressed genes that represent the corresponding stage of development. Furthermore, we identified conserved key members (or hub genes) of the goat transcriptional networks. Their association with other embryo genes suggests that they may have important regulatory roles in embryo development. Our cross-mammalian species transcriptomic comparisons demonstrate both conserved and goat-specific features of preimplantation development.

RNA-Seq Reveals miRNA Role Shifts in Seven Stages of Skeletal Muscles in Goat Fetuses and Kids

MicroRNAs (miRNAs) are indispensable for the regulation of skeletal muscle. We performed RNA sequencing (RNA-seq) to establish a comprehensive miRNA profiling of goats in seven stages, namely, 45- (F45), 65- (F65), 90- (F90), 120- (F120), and 135-day (F135) fetuses, newborn (B1), and 90-day-old (B90) kids. In total, 421 known miRNAs and 228 goat novel miRNAs were identified in the data, and the average abundance of 19 miRNAs in seven stages exceeds 10,000 reads per million. Furthermore, 420 differentially expressed miRNAs (DEmiRNAs) were identified in all comparison group at seven stages, 80 of which were uniquely differentially expressed in the B1 and B90 comparison groups. Pathway analysis indicated that this group was associated with the release of muscle hypertrophy and regulation of myoblast proliferation. Besides, 305 DEmiRNAs were clustered into three significantly enriched profiles (profiles 11, 16, and 19). Function analysis revealed that profile 16 was related to muscle hypertrophy and differentiation. Profile 11 was involved in multiple enzyme activities and metabolic processes in muscle cells. And profile 19 was involved in material transport and structural stability. Two highly expressed miRNAs and three key miRNAs (chi-miR-328-3p, chi-miR-767, and chi-miR-150) of these profiles were verified to be consistent with the data by quantitative real-time PCR. These results provided a catalog of goat muscle-associated miRNAs, allowing us to better understand the transformation of miRNA roles during mammalian muscle development.

Effect of GABA-T on Reproductive Function in Female Rats

This study explored the role of γ-aminobutyric acid transaminase (GABA-T) in the puberty and reproductive performance of female rats. Immunofluorescence technique, quantitative real-time PCR (RT-qPCR) and enzyme-linked immunosorbent assay (ELISA) were used to detect the distribution of GABA-T and the expression of genes and hormones in female rats, respectively. The results showed that GABA-T was mainly distributed in the arcuate nucleus (ARC), paraventricular nucleus (PVN) and periventricular nucleus (PeN) of the hypothalamus, and in the adenohypophysis, ovarian granulosa cells and oocytes. Abat mRNA level at 28 d was lowest in the hypothalamus and the pituitary; at puberty, it was lowest in the ovary. Abat mRNA level was highest in adults in the hypothalamus; at infancy and puberty, it was highest in the pituitary; and at 21 d it was highest in the ovary. After vigabatrin (GABA-T irreversible inhibitor) was added to hypothalamus cells, the levels of Abat mRNA and Rfrp-3 mRNA were significantly reduced, but Gnrh mRNA increased at the dose of 25 and 50 μg/mL; Kiss1 mRNA was significantly increased but Gabbr1 mRNA was reduced at the 50 μg/mL dose. In prepubertal rats injected with vigabatrin, puberty onset was delayed. Abat mRNA, Kiss1 mRNA and Gnrh mRNA levels were significantly reduced, but Rfrp-3 mRNA level increased in the hypothalamus. Vigabatrin reduced the concentrations of GABA-T, luteinizing hormone (LH) and progesterone (P₄), and the ovarian index. Lactation performance was reduced in adult rats with vigabatrin treatment. Four hours after vigabatrin injection, the concentrations of GABA-T and LH were significantly reduced in adult and 25 d rats, but follicle-stimulating hormone (FSH) increased in 25 d rats. In conclusion, GABA-T affects the reproductive function of female rats by regulating the levels of Gnrh, Kiss1 and Rfrp-3 in the hypothalamus as well as the concentrations of LH and P₄.

Trend analysis of the role of circular RNA in goat skeletal muscle development

Background

Circular RNA (circRNA) is produced during the splicing of mRNA (in addition to linear splicing) and is part of the gene regulatory network. The temporal expression patterns the different developmental stages were inseparable from these molecules’ function.

Results

Skeletal muscles of Anhui white goat (AWG) across seven fetal to postnatal development stages were sequenced and 21 RNA sequencing libraries were constructed. We thereby identified 9090 circRNAs and analyzed their molecular properties, temporal expression patterns, and potential functions at the different stages. CircRNAs showed complexities and diversity of formation as the same host gene produces multiple isoforms of these nucleic acids with different expression profiles. The differential expression of 2881 circRNAs (DECs, P < 0.05) was identified and four were randomly selected and validated by qPCR. Moreover, 1118 DECs under strict selected (SDECs, |log2^FC| > 2 and P-adj value < 0.01) showed 4 expression trends (Clusters 0, 19, 16 and 18). Cluster 0 molecules had increasing expression at all stages with effects on muscle through metabolism, regulation of enzyme activity, and biosynthesis. Cluster 16 circRNAs had high expression in the early and late stages and are involved in “Wnt signaling pathway”, “AMPK signaling pathway” and others. Cluster 18 molecules were mainly expressed at F120 and participate in “cytoskeletal protein binding”, “Notch signaling pathway” and so on. Cluster 19 circRNAs were down-regulated at all stages and related to muscle structure and development. Lastly, the SDECs divided the period of skeletal muscle development into three transitional stages: stage 1 (F45 to F90), which related to muscle satellite cell proliferation and muscle fiber structure; stage 2 (F90 to B1), in which the attachment of the cytoplasmic surface to the actin cytoskeleton initiates; and stage 3, which involved the “cGMP-PKG signaling pathway”. Moreover, the paraffin sections messages also validated that there are three transitional stages of skeletal muscle development.

Conclusion

Our current study provides a catalog of goat muscle-related circRNAs that can stratify skeletal muscle development fetus 45 days to newborn 90 days into three developmental stages. These findings better our understanding of functional transitions during mammalian muscle development.

Local origin or external input: modern horse origin in East Asia

BACKGROUND

Despite decades of research, the horse domestication scenario in East Asia remains poorly understood.

RESULTS

The study identified 16 haplogroups with fine-scale phylogenetic resolution using mitochondrial genomes of 317 horse samples. The time to the most recent common ancestor of the 16 haplogroups ranges from [0.8-3.1] thousand years ago (KYA) to [7.9-27.1] KYA. With combined analyses of the mitochondrial control region for 35 extant Przewalski's horses, 3544 modern and 203 ancient horses across the world, researchers provide evidence for that East Asian prevalent haplogroups Q and R were indigenously domesticated or they were involved in numerous distinct genetic components from wild horses in the southern part of East Asia. These events of haplotypes Q and R occurred during 4.7 to 16.3 KYA and 2.1 to 11.5 KYA, respectively. The diffusion of preponderant European haplogroups L from west to East Asia is consistent with the external gene input. Furthermore, genetic differences were detected between northern East Asia and southern East Asia cohorts by Principal Component Analysis, Analysis of Molecular Variance test, the χ² test and phylogeographic analyses.

CONCLUSIONS

All results suggest a complex picture of horse domestication, as well as geographic pattern in East Asia. Both local origin and external input occurred in East Asia horse populations. And besides, there are at least two different domestication or hybridization centers in East Asia.

Comprehensive Analysis of LncRNA Reveals the Temporal-Specific Module of Goat Skeletal Muscle Development

A series of complex processes regulate muscle development, and lncRNAs play essential roles in the regulation of skeletal myogenesis. Using RNA sequencing, we profiled the lncRNA expression during goat (Capra hircus) skeletal muscle development, which included seven stages across fetal 45 (F45), 65 (F65), 90 (F90), 120 (F120), 135 (F135) days, born for 24 h (B1) and 90 (B90) days. A total of 15,079 lncRNAs were identified in the seven stages, and they were less conservative with other species (human, cow, and mouse). Among them, 547 were differentially expressed, and they divided the seven stages into three functional transition periods. Following weighted gene co-expression network analysis (WGCNA), five lncRNA modules specific for developmental stages were defined as three types: 'Early modules', 'late modules', and 'individual-stage-specific modules'. The enrichment content showed that 'early modules' were related to muscle structure formation, 'late modules' participated in the 'p53 signaling pathway' and other pathways, the F90-highly related module was involved in the 'MAPK signaling pathway', and other pathways. Furthermore, we identified hub-lncRNA in three types of modules, and LNC_011371, LNC_ 007561, and LNC_001728 may play important roles in goat skeletal muscle. These data will facilitate further exploration of skeletal muscle lncRNA functions at different developmental stages in goats.

The expression and regulation of miR-1 in goat skeletal muscle and satellite cell during muscle growth and development

MicroRNA-1 (miR-1) has been shown to play an important role in muscle growth and development, however, it was mainly discovered in model animals. To explore the function and mechanism of miR-1 in goat, we firstly explored the expression profile of miR-1 in goat tissues and cells. Furthermore, the target gene of miR-1 was predicted, and the relationship between miR-1 and one of its target genes, histone deacetylase 4 (HDAC4), was analyzed through double luciferase reporter assay, real-time PCR, and western blot. It was found that the miR-1 is most abundantly expressed in goat heart and skeletal muscle tissue. Meanwhile, the expression of miR-1 showed an increasing tendency from new-born goats to the 7-month-old goats, and then its expression decreases as the goats mature further. In addition, the expression levels of miR-1 decreased in goat skeletal muscle satellite cells with the algebraic increasing of cells. At last, the results showed that HDAC4 is a target gene of miR-1 in goat, and miR-1 can inhibit the post-transcriptional expression of HDAC4, but had no significant influence on the mRNA level of HDAC4. It was hypothesized that miR-1 promotes muscle development by inhibiting the post-transcriptional expression of HDAC4 in goat.

Comparative analysis of differentially expressed genes between the ovaries from pregnant and nonpregnant goats using RNA-Seq

Background

A multitude of genes tightly regulate ovarian follicular development and hormone secretion. These complex and coordinated biological processes are altered during pregnancy. In order to further understand the regulatory role of these genes during pregnancy, it is important to screen the differentially expressed genes (DEGs) in the ovaries of pregnant and nonpregnant mammals. To detect the genes associated with the development of pregnancy in goats, DEGs from the ovaries from pregnant and nonpregnant Anhui white goats (pAWGs and nAWGs, respectively) were analyzed using RNA sequencing technology (RNA-Seq).

Results

In this study, 13,676,394 and 13,549,560 clean reads were generated from pAWGs and nAWGs, respectively, and 1724 DEGs were identified between the two libraries. Compared with nAWGs, 1033 genes were upregulated and 691 genes were downregulated in pAWGs, including PGR, PRLR, STAR and CYP19A1, which play important roles in goat reproduction. Gene Ontology analysis showed that the DEGs were enriched for 49 functional GO terms. Kyoto Encyclopedia of Genes and Genomes analysis revealed that 397 DEGs were significantly enriched in 13 pathways, including “cell cycle”, “cytokine–cytokine receptor interaction” and “steroid biosynthesis”, suggesting that the genes may be associated with cell cycle regulation, follicular development and hormone secretion to regulate the reproduction process. Additionally, quantitative real-time PCR was used to verify the reliability of the RNA-Seq data.

Conclusions

The data obtained in this work enrich the genetic resources of goat and provide a further understanding of the complex molecular regulatory mechanisms occurring during the development of pregnancy and reproduction in goats. The DEGs screened in this study may play an important role in follicular development and hormone secretion and they would provide scientific basis for related research in the future.

Electronic supplementary material

The online version of this article (10.1186/s40709-019-0095-9) contains supplementary material, which is available to authorized users.

Filtered reproductive long non-coding RNAs by genome-wide analyses of goat ovary at different estrus periods

Background

The goat is an important farm animal. Reproduction is an important process of goat farming. The ovary is the most important reproductive organ for goats. In recent years, an increasing number of long non-coding RNAs (lncRNAs) have been implicated in the regulation of mammal reproduction. However, there are few studies on the function of lncRNAs in reproduction, particularly lncRNAs in the ovary.

Results

The sequencing of goat ovaries generated 1,122,014,112 clean reads, and 4926 lncRNAs and 1454 TUCPs (transcripts of uncertain coding potential) were identified for further analysis by using the coding potential analysis software, CNCI, CPC and Pfam-sca. There were 115 /22 differential lncRNAs /TUCPs transcripts between the ovaries of the luteal phase and the follicular phase. We predicted the related genes of lncRNA /TUCP based on co-expression and co-localization methods. In total, 2584 /904 genes were predicted by co-expression, and 326/73 genes were predicted by co-localization. The functions of these genes were further analyzed with GO and KEGG analysis. The results showed that lncRNAs /TUCPs, which are highly expressed in goat ovaries in the luteal phase, are mainly associated with the synthesis of progesterone, and we filtered the lncRNAs /TUCPs, such as XR_001918177.1 and TUCP_001362, which may regulate the synthesis of progesterone; lncRNAs /TUCPs, which are highly expressed in goat ovaries in the follicular phase, are mainly associated with oogenesis and the maturation of oocytes, and we filtered the lncRNAs /TUCPs that may regulate the oogenesis and maturation of oocyte, such as XR_001917388.1 and TUCP_000849.

Conclusion

The present study provided the genome expression profile of lncRNAs /TUCPs in goat ovaries at different estrus periods and filtered the potential lncRNAs /TUCPs associated with goat reproduction. These results are helpful to further study the molecular mechanisms of goat reproduction.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5268-7) contains supplementary material, which is available to authorized users.

The interaction between DNA methylation and long non-coding RNA during the onset of puberty in goats

Epigenetics plays an important role in controlling female puberty. Both DNA methylation and long non-coding RNAs (lncRNA) regulate the initiation of puberty by affecting the expression of genes related to puberty. While recent studies have indicated that DNA methylation of lncRNA represses the expression of lncRNA, its role in regulating puberty remains unclear. To explore the mechanism between DNA methylation and lncRNAs during puberty onset, we performed whole-genome bisulphite sequencing (WGBS) and RNA-sequencing (RNA-seq). We found that DNA methylation was inversely correlated to gene expression levels during puberty. Methylation levels gradually decreased near the transcription initiation site and were present at high levels in the exon, intron and 3' untranslated regions. In the promoter, lncRNA expression was negatively related to DNA methylation. We reported hypermethylation in the gene body and downstream of the lncRNA compared with upstream regions. In GO and KEGG analyses, we found enriched target genes of lncRNA, XLOC_960044 and XLOC_767346. During puberty, methylation of these genes increased while expression decreased. Our study indicates that DNA methylation of the promoter is negatively correlated with lncRNA during puberty onset, and methylation regulates the initiation of puberty via lncRNA, which provides new insight into the epigenetic mechanism of puberty onset.

RNA-seq analysis of lncRNA-controlled developmental gene expression during puberty in goat & rat

Background

Puberty is a pivotal stage in female animal development, and marks the onset of reproductive capability. However, little is known about the function of lncRNAs (long noncoding RNAs) in puberty. Therefore, RNA-seq analysis were performed between goats and rats to clarify the roles of lncRNAs and mRNAs in the onset of puberty.

Results

In the present study, the length of lncRNAs, the length of the open reading frame and the exon count were compared between the two species. Furthermore, functional annotation analysis based on Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analysis of lncRNAs target genes and differentially expressed mRNA demonstrated the significantly enriched terms, such as AMPK signaling pathway, oxytocin signaling pathway, insulin secretion as well as pheromone receptor activity, and some other signaling pathways which were involved in the regulation of female puberty. Moreover, our results of siRNA interference in vitro showed the candidate lncRNA XLOC_446331 may play a crucial role in regulating female puberty.

Conclusion

In conclusion, the RNA-seq analysis between goat and rat provide novel candidate regulators for genetic and molecular studies on female puberty.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0608-9) contains supplementary material, which is available to authorized users.

Identification and analysis of differentially expressed long non-coding RNAs between multiparous and uniparous goat (Capra hircus) ovaries

Long non-coding RNAs (lncRNAs) play important roles in almost all biological processes. However, there is little information on the effects of lncRNAs on ovulation and lambing rates. In the present study, we used high-throughput RNA sequencing to identify differentially expressed lncRNAs between the ovaries of multiparous (Mul) and uniparous (Uni) Anhui White goats. Among the 107,255,422 clean reads, 183,754 lncRNAs were significantly differentially expressed between the Uni and Mul. Among them, 455 lncRNAs were co-expressed between the two samples, whereas, 157,523 lncRNAs were uniquely expressed in the Uni, and 25,776 uniquely lncRNAs were expressed in the Mul. Through Cis role analysis, 24 lncRNAs were predicted to overlap with cis-regulatory elements, which involved in Progesterone-mediated oocyte maturation, Steroid biosynthesis, Oocyte meiosis, and gonadotropin-releasing hormone (GnRH) signaling pathway. These 4 pathways were related to ovulation, and the KEGG pathway analysis on target genes of the differentially expressed lncRNAs confirmed this results. In addition, 10 lncRNAs harbored precursors of 40 miRNAs, such as TCONS_00320849 related to a mature miRNA sequence, miR-365a, which was reported to be related to proliferation, were annotated in the precursor analysis of miRNAs. The present expand the understanding of lncRNA biology and contribute to the annotation of the goat genome. The study will provide a resource for lncRNA studies of ovulation and lambing.

Screening and evaluating of long noncoding RNAs in the puberty of goats

Background

Long noncoding RNAs (lncRNAs) are involved in regulating animal development, however, their function in the onset of puberty in goats remain largely unexplored. To identify the genes controlling the regulation of puberty in goats, we measured lncRNA and mRNA expression levels from the hypothalamus.

Results

We applied RNA sequencing analysis to examine the hypothalamus of pubertal (case; n = 3) and prepubertal (control; n = 3) goats. Our results showed 2943 predicted lncRNAs, including 2012 differentially expressed lncRNAs, which corresponded to 5412 target genes. We also investigated the role of lncRNAs that act cis and trans to the target genes and found a number of lncRNAs involved in the regulation of puberty and reproduction, as well as several pathways related to these processes. For example, oxytocin signaling pathway, sterol biosynthetic process, and pheromone receptor activity signaling pathway were enriched as Kyoto Encyclopedia of Genes and Genomes (KEGG) or gene ontology (GO) analyses showed.

Conclusion

Our results clearly demonstrate that lncRNAs play an important role in regulating puberty in goats. However, further research is needed to explore the functions of lncRNAs and their predicted targets to provide a detailed expression profile of lncRNAs on goat puberty.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3578-9) contains supplementary material, which is available to authorized users.

DNA Methylation Patterns in the Hypothalamus of Female Pubertal Goats

Female pubertal development is tightly controlled by complex mechanisms, including neuroendocrine and epigenetic regulatory pathways. Specific gene expression patterns can be influenced by DNA methylation changes in the hypothalamus, which can in turn regulate timing of puberty onset. In order to understand the relationship between DNA methylation changes and gene expression patterns in the hypothalamus of pubertal goats, whole-genome bisulfite sequencing and RNA-sequencing analyses were carried out. There was a decline in DNA methylation levels in the hypothalamus during puberty and 268 differentially methylated regions (DMR) in the genome, with differential patterns in different gene regions. There were 1049 genes identified with distinct expression patterns. High levels of DNA methylation were detected in promoters, introns and 3′-untranslated regions (UTRs). Levels of methylation decreased gradually from promoters to 5′-UTRs and increased from 5′-UTRs to introns. Methylation density analysis demonstrated that methylation level variation was consistent with the density in the promoter, exon, intron, 5′-UTRs and 3′-UTRs. Analyses of CpG island (CGI) sites showed that the enriched gene contents were gene bodies, intergenic regions and introns, and these CGI sites were hypermethylated. Our study demonstrated that DNA methylation changes may influence gene expression profiles in the hypothalamus of goats during the onset of puberty, which may provide new insights into the mechanisms involved in pubertal onset.

Comparative profiling of differentially expressed microRNAs between the follicular and luteal phases ovaries of goats

To explore if the regulation at post-transcriptional level of follicular phase (Fols) to luteal phase (Luts) transition occurs in the ovaries of Anhuai goats, the differentially expressed microRNAs (miRNAs) of ovaries in the Fols and Luts were analyzed using Solexa sequencing in the study. In total, 320 known miRNAs were co-expressed in the two phases, 339 and 353 known miRNAs were expressed in the ovary in the Fols and Luts, respectively. In addition, 45 novel miRNAs were co-expressed in the two phases, 70 and 94 novel miRNAs were expressed in the ovary in the Fols and Luts, respectively. Let-7f was the highest expressed significantly different known miRNA in the two phases, and mir-159 was the highest expressed significantly different novel miRNA in the two phases, which may participate in the follicular-luteal transition of Anhuai goats. GO annotation and KEGG pathway analysis were applied to analyze the target genes of differentially expressed miRNAs detected in the two phases. The results will help to further understand the role of miRNAs in the regulation of follicular to luteal transition in goat ovaries.

Expression profiles of differentially expressed genes affecting fecundity in goat ovarian tissues

Although RNA-Seq is an effective method for identifying and exploring novel functional genes in mammals, it has rarely been applied to study fertility-related genes in the goat. In this study, RNA-Seq was used to screen the estrus ovaries of uniparous and multiparous Anhui white goats (AWGs). In total, 15,890 genes were identified and 2201 of these were found to be differentially expressed between the genetic libraries from uniparous and multiparous goats. Compared to the uniparous library, 1583 genes were up-regulated and 618 genes were down-regulated in the multiparous library. The FER1L4 gene showed the level of highest up-regulation in the multiparous library, while SRD5A2 expression showed the greatest down-regulation. In order to determine the functions of FER1L4 and SRD5A2 in goats, the expression profiles of the two genes in different tissues from AWGs and Boer goats at diestrus were analyzed by quantitative PCR. FER1L4 and SRD5A2 showed tissue specific expression patterns and were highly expressed in ovaries from both AWGs and Boer goats. FER1L4 was more highly expressed in ovaries from multiparous than uniparous AWGs. In contrast, SRD5A2 was expressed at a lower level in multiparous AWGs. These results indicated that FER1L4 and SRD5A2 may be associated with the high fecundity of AWGs.

Genome-Wide Dynamic Profiling of Histone Methylation during Nuclear Transfer-Mediated Porcine Somatic Cell Reprogramming

The low full-term developmental efficiency of porcine somatic cell nuclear transfer (SCNT) embryos is mainly attributed to imperfect epigenetic reprogramming in the early embryos. However, dynamic expression patterns of histone methylation involved in epigenetic reprogramming progression during porcine SCNT embryo early development remain to be unknown. In this study, we characterized and compared the expression patterns of multiple histone methylation markers including transcriptionally repressive (H3K9me2, H3K9me3, H3K27me2, H3K27me3, H4K20me2 and H4K20me3) and active modifications (H3K4me2, H3K4me3, H3K36me2, H3K36me3, H3K79me2 and H3K79me3) in SCNT early embryos from different developmental stages with that from in vitro fertilization (IVF) counterparts. We found that the expression level of H3K9me2, H3K9me3 and H4K20me3 of SCNT embryos from 1-cell to 4-cell stages was significantly higher than that in the IVF embryos. We also detected a symmetric distribution pattern of H3K9me2 between inner cell mass (ICM) and trophectoderm (TE) in SCNT blastocysts. The expression level of H3K9me2 in both lineages from SCNT expanded blastocyst onwards was significantly higher than that in IVF counterparts. The expression level of H4K20me2 was significantly lower in SCNT embryos from morula to blastocyst stage compared with IVF embryos. However, no aberrant dynamic reprogramming of H3K27me2/3 occurred during early developmental stages of SCNT embryos. The expression of H3K4me3 was higher in SCNT embryos at 4-cell stage than that of IVF embryos. H3K4me2 expression in SCNT embryos from 8-cell stage to blastocyst stage was lower than that in the IVF embryos. Dynamic patterns of other active histone methylation markers were similar between SCNT and IVF embryos. Taken together, histone methylation exhibited developmentally stage-specific abnormal expression patterns in porcine SCNT early embryos.

Inflammasome activity is essential for one kidney/deoxycorticosterone acetate/salt-induced hypertension in mice

BACKGROUND AND PURPOSE

Inflammasomes are multimeric complexes that facilitate caspase-1-mediated processing of the pro-inflammatory cytokines IL-1β and IL-18. Clinical hypertension is associated with renal inflammation and elevated circulating levels of IL-1β and IL-18. Therefore, we investigated whether hypertension in mice is associated with increased expression and/or activation of the inflammasome in the kidney, and if inhibition of inflammasome activity reduces BP, markers of renal inflammation and fibrosis.

EXPERIMENTAL APPROACH

Wild-type and inflammasome-deficient ASC(-/-) mice were uninephrectomized and received deoxycorticosterone acetate and saline to drink (1K/DOCA/salt). Control mice were uninephrectomized but received a placebo pellet and water. BP was measured by tail cuff; renal expression of inflammasome subunits and inflammatory markers was measured by real-time PCR and immunoblotting; macrophage and collagen accumulation was assessed by immunohistochemistry.

KEY RESULTS

1K/DOCA/salt-induced hypertension in mice was associated with increased renal mRNA expression of inflammasome subunits NLRP3, ASC and pro-caspase-1, and the cytokine, pro-IL-1β, as well as protein levels of active caspase-1 and mature IL-1β. Following treatment with 1K/DOCA/salt, ASC(-/-) mice displayed blunted pressor responses and were also protected from increases in renal expression of IL-6, IL-17A, CCL2, ICAM-1 and VCAM-1, and accumulation of macrophages and collagen. Finally, treatment with a novel inflammasome inhibitor, MCC950, reversed hypertension in 1K/DOCA/salt-treated mice.

CONCLUSIONS AND IMPLICATIONS

Renal inflammation, fibrosis and elevated BP induced by 1K/DOCA/salt treatment are dependent on inflammasome activity, highlighting the inflammasome/IL-1β pathway as a potential therapeutic target in hypertension.

Characterization of microRNAs from goat (Capra hircus) by Solexa deep-sequencing technology

MicroRNAs (miRNAs) are an important class of small noncoding RNAs that are highly conserved in plants and animals. Many miRNAs are known to mediate a myriad of cell processes, including proliferation and differentiation, via the regulation of some transcription and signaling factors, which are closely related to muscle development and disease. In this study, small RNA cDNA libraries of Boer goats were constructed. In addition, we obtained the goat muscle miRNAs by using Solexa deep-sequencing technology and analyzed these miRNA characteristics by combining it with the bioinformatics technology. Based on Solexa sequencing and bioinformatics analysis, 562 species-conserved and 5 goat genome-specific miRNAs were identified, 322 of which exceeded 100 in the expression levels. The results of real-time quantitative polymerase chain reaction from 8 randomly selected miRNAs showed that the 8 miRNAs were expressed in goat muscle, and the expression patterns were consistent with the Solexa sequencing results. The identification and characterization of miRNAs in goat muscle provide important information on the role of miRNA regulation in muscle growth and development. These data will help to facilitate studies on the regulatory roles played by miRNAs during goat growth and development.

Temporal and spatial expression of podocyte-associated molecules are accompanied by proteinuria in IgA nephropathy rat model

We used a rat model to assess the role of nephrin, podocin, and desmin in the pathogenesis of IgA nephropathy (IgAN). A rat IgAN model was established by administration of BSA, CCl(4), and lipopolysaccharide (LPS) and compared with healthy control rats. Urinary protein, urine red blood cells, and biochemical parameters were measured for 12 weeks. Renal morphology and ultrastructure were examined by light and electron microscopy. Immunofluorescence was used to assess IgA deposition in the glomeruli and to measure expression of nephrin, podocin, and desmin. Real-time quantitative PCR was used to measure expression of nephrin, podocin, and desmin mRNAs. IgAN rats developed proteinuria at week-6 and this worsened over time. Pathological changes were evident under light microscopy at week-8 and under electron microscopy at week-4. Immunofluorescence analysis showed deposition of IgA in the kidneys of IgAN rats, but not control rats. IgAN rats had increased expression of glomerular podocin, nephrin, and desmin mRNAs and proteins at week-4. The expression of nephrin, podocin and desmin proteins and the expression of podocin and desmin mRNAs preceded the increase in urinary protein. Taken together, our study of a rat model of IgAN indicates that changes in the expression and distribution of nephrin, podocin, and desmin precede and may cause foot process fusion and proteinuria.

Genetic differentiation of chinese indigenous meat goats ascertained using microsatellite information

To investigate the genetic diversity of seven Chinese indigenous meat goat breeds (Tibet goat, Guizhou white goat, Shannan white goat, Yichang white goat, Matou goat, Changjiangsanjiaozhou white goat and Anhui white goat), explain their genetic relationship and assess their integrity and degree of admixture, 302 individuals from these breeds and 42 Boer goats introduced from Africa as reference samples were genotyped for 11 microsatellite markers. Results indicated that the genetic diversity of Chinese indigenous meat goats was rich. The mean heterozygosity and the mean allelic richness (AR) for the 8 goat breeds varied from 0.697 to 0.738 and 6.21 to 7.35, respectively. Structure analysis showed that Tibet goat breed was genetically distinct and was the first to separate and the other Chinese goats were then divided into two sub-clusters: Shannan white goat and Yichang white goat in one cluster; and Guizhou white goat, Matou goat, Changjiangsanjiaozhou white goat and Anhui white goat in the other cluster. This grouping pattern was further supported by clustering analysis and Principal component analysis. These results may provide a scientific basis for the characteristization, conservation and utilization of Chinese meat goats.

Evaluation of the genetic diversity and population structure of Chinese indigenous horse breeds using 27 microsatellite markers

We determined the genetic diversity and evolutionary relationships among 26 Chinese indigenous horse breeds and two introduced horse breeds by genotyping these animals for 27 microsatellite loci. The 26 Chinese horse breeds come from 12 different provinces. Two introduced horse breeds were the Mongolia B Horse from Mongolia and the Thoroughbred Horse from the UK. A total of 330 alleles were detected, and the expected heterozygosity ranged from 0.719 (Elenchuns) to 0.780 (Dali). The mean number of alleles among the horse breeds ranged from 6.74 (Hequ) to 8.81 (Debao). Although there were abundant genetic variations found, the genetic differentiation was low between the Chinese horses, which displayed only 2.4% of the total genetic variance among the different breeds. However, genetic differentiation (pairwise FST) among Chinese horses, although moderate, was still apparent and varied from 0.001 for the Guizou-Luoping pair to 0.064 for the Jingjiang-Elenchuns pair. The genetic differentiation patterns and genetic relationships among Chinese horse breeds were also consistent with their geographical distribution. The Thoroughbred and Mongolia B breeds could be discerned as two distinct breeds, but the Mongolia B horse in particular suffered genetic admixture with Chinese horses. The Chinese breeds could be divided into five major groups, i.e. the south or along the Yangtze river group (Bose, Debao, Wenshan, Lichuan, Jianchang, Guizhou, Luoping, Jinjiang and Dali), the Qinghai-Tibet Plateau group (Chaidamu, Hequ, Datong, Yushu, Tibet Grassland and Tibet Valley), the Northeast of China group (Elenchuns, Jilin and Heihe), the Northwest of China group (Kazakh, Yili and Yanqi) and the Inner Mongolia group (Mongolia A, Sanhe, Xinihe,Wuzhumuqin and Sengeng). This grouping pattern was further supported by principal component analysis and structure analysis.

Identification of Y chromosome genetic variations in Chinese indigenous horse breeds

Y chromosome acts as a single nonrecombining unit that is male specific and in effect haploid, thus ensuring the preservation of mutational events as a single haplotype via male lines. In this study, 6 Y chromosome-specific microsatellites (SSR) were tested for the patrilineal genetic variations of 573 male samples from Chinese domestic horse (30 breeds), Przewalski's horse, and donkey. All the 6 loci appeared as a haplotype block in Przewalski's horse and the domestic donkey. There were notable differences, however, at Y chromosome markers between horse and donkey. There were 2 haplotypes of Eca.YA16 in the domestic horse breeds, Haplotype A (Allele A: 156 bp) and Haplotype B (Allele B: 152 bp). Allele A was the common allele among 30 horse breeds, and Allele B was found in 11 horse breeds. This is the first description of a Y chromosome variant for horses. The 2 haplotypes of Y chromosome discovered in the domestic horse breeds in China could be helpful in unveiling their intricate genetic genealogy.

Major late toxicities after conformal radiotherapy for nasopharyngeal carcinoma-patient- and treatment-related risk factors

PURPOSE

To retrospectively analyze the factors affecting late toxicity for nasopharyngeal carcinoma.

METHODS AND MATERIALS

Between 1998 and 2003, 422 patients were treated with a conformal technique with 2-Gy daily fractions to a total dose of 70 Gy. Conventional fractionation (5 fractions weekly) was used in 232 patients and accelerated fractionation (6 fractions weekly) in 190 patients. One hundred seventy-one patients were treated with the basic radiotherapy course alone (Group 1), 55 patients had an additional boost of 5 Gy in 2 fractions (Group 2), and 196 patients underwent concurrent cisplatin-based chemotherapy (Group 3).

RESULTS

The 5-year overall toxicity rate was significantly greater in Group 3 than in Group 1 (37% vs. 27%, p = 0.009). Although the overall rate in Group 2 was not elevated (28% vs. 27%, p = 0.697), a significant increase in temporal lobe necrosis was observed (4.8% vs. 0%, p = 0.015). Multivariate analyses showed that age and concurrent chemotherapy were significant factors. The hazard ratio of overall toxicity attributed to chemotherapy was 1.99 (95% confidence interval, 1.32-2.99, p = 0.001). The mean radiation dose to the cochlea was another significant factor affecting deafness, with a hazard ratio of 1.03 (95% confidence interval, 1.01-1.05, p = 0.005) per 1-Gy increase. The cochlea that received >50 Gy had a significantly greater deaf rate (Group 1, 18% vs. 7%; and Group 3, 22% vs. 14%).

CONCLUSION

The therapeutic margin for nasopharyngeal carcinoma is extremely narrow, and a significant increase in brain necrosis could result from dose escalation. The significant factors affecting the risk of deafness included age, concurrent chemoradiotherapy, and greater radiation dose to the cochlea.

Sensitivity to topoisomerase I inhibitors and cisplatin is associated with epidermal growth factor receptor expression in human cervical squamous carcinoma ME180 sublines

The relationship between expression and function of the epidermal growth factor (EGF) family of receptors and chemosensitivity remains controversial. We studied the chemosensitivity to various anticancer agents of human cervical squamous carcinoma ME180 cells, and two resistant subclones, ME180/TNF and ME180/Pt, which also differ in their EGF receptor (EGFR) expression. Compared with ME180 cells, EGFR is overexpressed sixfold in ME180/TNF cells and is barely detectable in ME 180/Pt cells. Cell cycle analysis by flow cytometry and BrdU incorporation into DNA showed a correlation between EGFR expression and percentage of cells in S phase and active DNA replication (35% in high EGFR-expressing ME180/TNF cells, 19% in non-EGFR-expressing ME180/Pt cells and 23% in parental, intermediate-level EGFR-expressing ME 180 cells). By MTT assay and compared with parental, intermediate-level EGFR-expressing ME180 cells, high EGFR-expressing ME180/TNF cells had a three- to fourfold increased sensitivity to cisplatin, camptothecin (CPT), and topotecan, and low EGFR-expressing ME180/Pt cells had a five- to ninefold reduced sensitivity to the same agents. In contrast, the degree of cross-resistance with the topoisomerase II inhibitors doxorubicin and etoposide was minimal and the pattern of sensitivity to the anti-microtubulin agents vinblastine and paclitaxel was different, with a two- to fourfold decreased sensitivity in the high EGFR-expressing ME180/TNF cells and only a 1.5-fold decreased sensitivity in the low EGFR-expressing ME180/Pt cells. Neither alterations in intracellular CPT levels nor changes in topoisomerase I expression or activity, measured as ability to form DNA-protein complexes, were found to explain the differences in sensitivity to CPT among the three cell lines. Co-treatment with CP358774, a specific EGFR tyrosine kinase inhibitor, reduced the enhanced sensitivity of high EGFR-expressing ME180/TNF cells to the values observed in intermediate EGFR-expressing ME180 cells, but only reduced modestly the sensitivity of intermediate expressing ME180 cells. As a result, the resistance index of low EGFR-expressing ME180/Pt cells compared with intermediate EGFR-expressing ME180 cells was reduced only from five- to fourfold for cisplatin and from seven- to fourfold for CPT when ME180 cells were exposed to CP358774. CP358774 did not affect the sensitivity to either agent in low EGFR-expressing ME180/Pt cells. These results provide evidence that changes in EGFR expression or function may play a role in determining chemosensitivity to platinum and topoisomerase I poisons in some human tumor systems, and that the EGFR-related changes in chemosensitivity may vary depending on the level of EGFR expression and/or function.

3-m-bromoacetylamino benzoic acid ethyl ester: a new cancericidal agent that activates the apoptotic pathway through caspase-9

The mechanism underlying the cancericidal activity of 3-m-bromoacetylamino benzoic acid ethyl ester (3-BAABE) was investigated. 3-BAABE exerted a strong cancericidal effect on human leukemia and lymphoma cells (IC(50) < 0.2 microgram/mL) and on cell lines of prostate, colon, ductal, and kidney cancer (IC(50) 0.8 to 0.88 microgram/mL). Multiple drug resistance (MDR) had no effect on the susceptibility of human lymphoma cells to 3-BAABE, since Daudi/MDR(20) and wild-type Daudi cells had a similar susceptibility to the cytotoxic effect of 3-BAABE. The cancericidal effect of 3-BAABE, which was not associated with changes in the cell cycle, was mediated by apoptosis. Thus, cells exposed to 3-BAABE displayed the DNA fragmentation ladder characteristic for apoptosis, associated with a marked increase of the activity of apoptosis effector caspases-3 and -6, which was followed by proteolytic cleavage of DNA fragmentation factor (DFF) and poly(ADP-ribose) polymerase (PARP). Exposure of tumor cells to 3-BAABE increased the activity of apical caspase-9, but had no effect on caspase-8. Complete inhibition of 3-BAABE-induced apoptosis was exerted by LEHD-FMK, a caspase-9 inhibitor. DEVD-FMK, a caspase-3 inhibitor, and VEID-FMK, a caspase-6 inhibitor, partially inhibited 3-BAABE-induced apoptosis, whereas exposure to IETD-FMK, a caspase-8 inhibitor, had no effect. The fragmentation and elevated activity of caspase-9 in 3-BAABE-treated cells and the fact that only an inhibitor of caspase-9 abrogated 3-BAABE-induced apoptosis indicate that 3-BAABE is a distinctive compound that elicits apoptosis through a pathway that is limited specifically to activation of apical caspase-9.

Response and determinants of sensitivity to paclitaxel in human non-small cell lung cancer tumors heterotransplanted in nude mice

The lack of tumor models that can reliably predict for response to anticancer agents remains a major deficiency in the field of experimental cancer therapy. Although heterotransplants of certain human solid tumors can be successfully grown in nude mice, they have never been appropriately explored for prediction of in vivo chemosensitivity to anticancer agents. We determined the tumor response rate and studied the influence of several biological and molecular tumor parameters on the in vivo sensitivity to paclitaxel in a series of heterotransplanted human non-small cell lung cancer (NSCLC) tumors. One hundred consecutive resected NSCLC tumors were heterotransplanted s.c. in nude mice. The in vivo sensitivity to i.v. paclitaxel (60 mg/kg every 3 weeks) was studied in 34 successfully grown heterotransplants. Treatment started when the tumors reached a size of 5 mm in diameter, and strict standard clinical criteria (>50% shrinkage in tumor weight or cross-sectional surface) were used to define tumor response. Baseline multidrug resistance protein (MRP), Her-2/neu, and epidermal growth factor receptor (EGFR) expression, and pre- and posttherapy bax and bcl-2 expression were determined by Western blot analysis. p53 status was determined by sequencing. The overall take rate was 46% (95% confidence interval, 36-56%) and was significantly higher (P < 0.05) for squamous carcinoma tumors (75%) than for adenocarcinoma tumors (30%) and bronchoalveolar tumors (23%). The heterotransplants were morphologically very similar to the original tumors. The response rate to paclitaxel was 21% (95% confidence interval, 9-38%). Baseline tumor parameters associated with response were no Her-2/neu expression (none of the responding tumors expressed Her-2/neu versus 48% of the nonresponding tumors, P = 0.05) and baseline bcl-2 expression (all responding tumors expressed bcl-2 versus only 43% of the nonresponding tumors, P = 0.02). There was a trend toward a higher response rate in bax-positive tumors, and MRP- and EGFR-negative tumors, but it was not statistically significant. The response was independent of baseline p53 status and baseline mitotic index. Responding tumors had a higher bax/bcl-2 ratio 24 h after therapy, but the difference was only marginally significant (2.8 for responding tumors versus 1.1 for nonresponding tumors, P = 0.07). The extent of mitotic arrest at 24 h after therapy was not associated with response. Human NSCLC heterotransplants are morphologically identical to the original tumors and have a response rate to paclitaxel that is equivalent to that reported in Phase II studies in patients with advanced NSCLC treated with single-agent paclitaxel. NSCLC heterotransplants deserve to be explored to evaluate new agents for lung cancer and to predict clinical response on an individual basis in selected groups of patients.

Development of cationic liposome formulations for intratracheal gene therapy of early lung cancer

Regional (intratracheal or aerosol) delivery of cationic liposome-DNA complexes for gene therapy of lung disease offers distinct advantages over systemic (intravenous) administration. However, optimal formulations for early lung cancer treatment have not been established. Therefore, we investigated >50 different liposome and micelle formulations for factors that may affect their transcription efficiency and tested the ideal formulations in an in vivo mouse model. Our data showed that cationic liposomes were generally more effective at transfecting genes than were micelles of the same lipid composition, thus suggesting a role for the bilayer structure in facilitating transfection. In addition, the transfection efficiency of liposome-delivered genes was highly dependent upon the lipid composition, lipid/DNA ratio, particle size of the liposome-DNA complex, and cell lines used. By optimizing these factors in vitro and in vivo, we developed a novel liposome formulation (DP3) suitable for intratracheal administration. Using G67 liposome as control, we found that DP3 was more effective than G67 in vitro and as effective as G67 at both preventing lung tumor growth and prolonging survival in our lung cancer mouse model. We observed a positive correlation between the in vitro p53 function and the in vivo antitumoral activities of liposome-p53 formulations, which had not been reported previously in studies of an intravenous liposome gene delivery system. This correlation may facilitate the development and optimization of a liposome-p53 formulation for aerosol use in lung cancer patients.

[Staining of complete denture: a preliminary clinical study]

OBJECTIVE:To study the relation between staining of complete denture and drinking tea, coffee and smoking, and to evaluate the result of denture cleaning agent on removing stain of complete denture.METHODS:A survey of 176 patients with complete dentures were carried out, which included the history of drinking tea, coffee and smoking, the use of denture-cleaning agent and the times of daily use. Statistical analysis was performed to determine the relationship between staining of denture and the above factors. RESULTS:The results showed (1)The longer the complete denture was weared, the more serious the staining. (2)Drinking tea and smoking were the main causes for denture staining. (3)Daily use of cleaning agent can effictively remove stain of denture. CONCLUSION:Change of life habit(drinking less tea and stop of smoking) and daily use of cleaning agent can reduce denture staining.

Induction of senescence-like phenotype and loss of paclitaxel sensitivity after wild-type p53 gene transfection of p53-null human non-small cell lung cancer H358 cells

The p53 gene plays an important role in the regulation of cell-cycle progression and apoptosis. Recent studies have implicated p53 in determining cell fate, and shown that p53 status is associated with cellular sensitivity to anticancer agents. However, the role of p53 in paclitaxel-induced cytotoxicity remains unclear. Here we show that the induction of exogenous wild-type (wt) p53 genes in p53-null human NSCLC H358 cells via transient gene transfection with cationic liposome-wt p53 complexes resulted in a typical senescence-like phenotype. In short, cell growth was reduced, homeostasis occurred, cell morphology became enlarged and flat, the cell cycle was arrested at G1 phase, cyclin B1 and cdc2 expression was down-regulated, and DNA synthesis was suppressed. The sensitivity of wt p53-transfected cells (H358/p53) to paclitaxel was approximately 3-fold lower than that of H358 cells. Paclitaxel treatment gradually and significantly blocked cell-cycle progression at G2/M phase and increased the accumulation of cyclin B1 and cdc2 in H358 cells. In contrast, the same treatment slightly arrested the cell cycle at G2/M phase and slightly elevated cyclin B1 expression in H358/p53 cells. The rate of uptake and efflux of paclitaxel was not significantly different between H358 and H358/p53 cells, indicating that the reduction in cellular sensitivity caused by p53 transfection was not due to alterasion in intracellular drug concentration. Together, our findings suggest that the induction of exogenous wt p53 gene expression in cells lacking p53 function can trigger the senescence program and that loss of sensitivity to paclitaxel by p53-transfected cells may be associated, at least in part, with the induction of a senescence-like phenotype.

EGF receptor and p21WAF1 expression are reciprocally altered as ME-180 cervical carcinoma cells progress from high to low cisplatin sensitivity

Cell cycle regulators and signal transduction pathways can influence apoptotic sensitivity of tumor cells, and we previously described an association between EGFr overexpression, reduced DNA repair activity, and increased apoptotic sensitivity of ME-180 cervical carcinoma cells toward cis-diammedichloroplatinum (cDDP; K. Nishikawa, et al., Cancer Res., 52: 4758-4765, 1992). In the present study, the characteristics of ME-180 cells selected for high or low apoptotic sensitivity to cDDP (or camptothecin) were examined and compared to determine whether signal transduction components and cell cycle regulation were distinct in these isogenic drug response variant populations. As ME-180 cells progressed from high to low cDDP sensitivity [IC50 approximately 80 ng/ml in cDDP sensitive (PT-S) to approximately 2000 ng/ml in cDDP-resistant (Pt-R) cells], there was a significant decrease in EGFr expression that paralleled the relative reduction in cDDP apoptotic responsiveness (approximately 30-fold). cDDP-resistant cells had the slowest rate of growth and more effectively reduced DNA adduct levels following cDDP exposure than parental cells. Cellular levels of the cell cycle inhibitor p21WAF1 inversely correlated with cDDP responsiveness with high levels of p21WAF1 expressed in drug-resistant Pt-R cells in the absence of elevated p53. cDDP stimulated a 2-fold increase in p53 levels in both drug-sensitive and drug-resistant cells but caused a delayed reduction in p21WAF1 levels, suggesting p53-independent regulation of p21WAF1 in ME-180 cells. Activation of EGFr in Pt-R cells stimulated cell cycle progression (2-fold), reduced p21WAF1 levels (>2-fold), and increased sensitivity to cDDP (3-fold), suggesting that receptor signaling enhanced the efficacy of cDDP to induce cell death by relieving cell cycle restriction. These results demonstrate that the transition of ME-180 cells from a drug-sensitive to drug-resistant phenotype correlates with reciprocal changes in EGFr and p21WAF1 expression and provides additional evidence that the pathways controlled by these proteins may contribute to some forms of drug resistance.

Novel suicide ligands of tubulin arrest cancer cells in S-phase

It is presently accepted that the mechanism of action for all anti-tumor tubulin ligands involves the perturbation of microtubule dynamics during the G2/M phase of cell division and subsequent entry into apoptosis [1]. In this report, we challenge the established dogma by describing a unique mechanism of action caused by a novel series of tubulin ligands, halogenated derivatives of acetamido benzoyl ethyl ester. We have developed a suicide ligand for tubulin, which covalently attaches to the target and shows potent cancericidal activity in tissue culture assays and in animal tumor models. These compounds target early S-phase at the G1/S transition rather than the G2/M phase and mitotic arrest. Bcl-2 phosphorylation, a marker of mitotic microtubule inhibition by other tubulin ligands was dramatically altered, phosphorylation was rapid and biphasic rather than a slow linear event. The halogenated ethyl ester series of derivatives thus constitute a unique set of tubulin ligands which induce a novel mechanism of apoptosis.