Genetic analysis of key agronomic traits of local sheep breeds in Xinjiang, China.

The formation of sheep (Ovis aries) breeds is influenced by different ecological environments and populations with different living habits, resulting in the development of germplasm resources with stable genetic key agronomic traits. Thus, investigating the genetic mechanisms behind various agronomic traits can enhance the conservation and utilization of diverse sheep breeds. Here, we explored the sheep variome and selection signatures using the Ovine Infinium HD SNP BeadChip (600 K SNPs) from 23 sheep breeds, comprising a total of 1215 sheep. The genetic mechanisms of wool quality and tail morphology were analyzed by selective sweep and genome-wide association study. Based on the results of within-population selective sweep analysis, we performed gene network analysis and divided them into 6 gene communities. We identified genetic regions containing genes linked to sheep wool and tail, which have been and may continue to be important targets for breeding and selection. Furthermore, our results revealed the expression profiles of genes in these regions across different biological systems. Our study provides insights into categorizing sheep breeds into distinct gene communities, as well as references for constructing genetic network pathways related to key agronomic traits in sheep and other domestic animals.

Landscape genomics reveals adaptive divergence of indigenous sheep in different ecological environments of Xinjiang, China.

Sheep are important livestock animals that have evolved under various ecological pressures. Xinjiang is a region with diverse and harsh environments that have shaped many local sheep breeds with unique characteristics and environmental adaptability. However, these breeds are losing ecological flexibility due to the promotion of intensive farming practices. Here we sequenced 14 local sheep breeds from Xinjiang and analyzed their genetic structure and gene flow with other sheep breeds from neighboring regions. The Tibetan Plateau was the geographic origin of Xinjiang native sheep evolution. We performed genome-environment association analysis and identified Bio9: Mean Temperature of Driest Quarter and Bio15: Precipitation Seasonality as the key environmental factors affecting Xinjiang local sheep and the key genes involved in their survival and adaptation. We classified Xinjiang native sheep breeds into six groups based on their differential genes by pairwise selective sweep analysis and Community Network Analysis. We analyzed transcriptome expression data of 832 sheep tissues and detected tissue-specific enrichment of six group-specific genes in different biological systems. Our results revealed the genetic basis of year-round estrus, drought tolerance, hypoxia resistance, and cold tolerance traits of Xinjiang sheep breeds. Moreover, we proposed conservation strategies for Xinjiang local sheep breeds and provided theoretical guidance for breeding new sheep breeds under global extreme environments.

Population structure, genetic diversity and prolificacy in pishan red sheep under an extreme desert environment.

Extreme environmental conditions are a major challenge for livestock production. Changes in climate conditions, especially those that lead to extreme weather, can reduce livestock production. The screening of genes and molecular markers is of great significance to explore the genetic mechanism of sheep prolificacy traits in Taklimakan Desert environment. We selected healthy adult Pishan Red Sheep (PRS) and Qira Black Sheep (QR) which live in Taklimakan Desert environment, collected blood from jugular vein, extracted DNA, and prepared Illumina Ovine SNP50 chip. For PRS, linkage disequilibrium (LD) was calculated using the ovine SNP50 Beadchip and the effective population size (Ne) was estimated using SMC++. The genetic characteristics of PRS were analyzed by integrated haplotype score (iHS) and fixation index (F ST ). The result showed that r 2 of PRS was 0.233 ± 0.280 in the range of 0-10 Kb and decreased with increasing distances. SMC++ tested that the Ne of PRS remained at 236.99 in recent generations. 184 genes were screened out under iHS 1% threshold, and 1148 genes were screened out with F ST under the 5% threshold, and 29 genes were obtained from the intersection of the two gene sets. In this study, the genetic characteristics of PRS and QR were compared by ovine genome chip, and the related excellent genes were searched, providing reference for the protection of sheep germplasm resources and molecular breeding in a desert environment.

Exploring the growth trait molecular markers in two sheep breeds based on Genome-wide association analysis.

Growth traits are quantitative traits controlled by multiple micro-effect genes. we identified molecular markers related to sheep growth traits, which formed the basis of molecular breeding. In this study, we randomly selected 100 Qira Black sheep and 84 German Merino sheep for the blood collection the jugular vein to genotype by using the Illumina Ovine SNP 50K Bead Chip. quality control criteria for statistical analysis were: rejection detection rate < 90% and minimum allele frequency (MAF) < 5%. Then, we performed Genome-wide association studies (GWAS) on sheep body weight, body height, body length, and chest circumference using mixed linear models. After getting 55 SNPs with significant correlation, they were annotated by reference genome of Ovis aries genome (Oar_v4.0) and We obtained a total of 84 candidate genes associated with production traits (BMPR1B, HSD17B3, TMEM63C, etc.). We selected BMPR1B for population validation and found a correlation between the FecB locus and body weight traits. Therefore, this study not only supplements the existing knowledge of molecular markers of sheep growth traits, but also has important theoretical significance and reference value for the mining of functional genes of sheep growth traits.

Analysis of Methylation and mRNA Expression of Lin28B Gene Promoter Region in the Hypothalamus of Dolang Sheep During Pubertal Initiation.

Lin28B plays an important role in puberty initiation in sheep. This study aimed to discuss the correlation between different growth periods and the methylation status of cytosine-guanine dinucleotide (CpG) islands in the promoter region of the Lin28B gene in the Dolang sheep's hypothalamus. In this study, the sequence of the Lin28B gene promoter region in Dolang sheep was obtained by cloning and sequencing, and methyl groups of the CpG island of the Lin28B gene promoter in the hypothalamus were detected by bisulfite sequencing PCR during the three periods of prepuberty, adolescence, and postpuberty in Dolang sheep. Lin28B expression in the Dolang sheep's hypothalamus was detected by fluorescence quantitative PCR at three stages: prepuberty, puberty, and postpuberty. In this experiment, the 2993-bp Lin28B promoter region was obtained, and it was predicted that there was a CpG island containing 15 transcription factor binding sites and 12 CpG sites, which may play a role in gene expression regulation. Overall, methylation levels increased from prepuberty to postpuberty, while Lin28B expression levels decreased, indicating that Lin28B expression was negatively correlated with promoter methylation levels. Variance analysis showed significant differences in the methylation status of CpG5, CpG7, and CpG9 between pre- and postpuberty (p < 0.05). Our data show that Lin28B expression is increased by demethylation of promoter CpG islands, with CpG5, CpG7, and CpG9 implicated as critical regulatory sites.

Identification of Signatures of Selection for Litter Size and Pubertal Initiation in Two Sheep Populations.

Fecundity is an important economic trait in sheep that directly affects their economic and productive efficiency. Our study aimed to identify SNP loci associated with sheep puberty or litter size which could be used in future breeding programs to improve fertility. Genomic DNA was obtained from Hetian and Cele Black sheep breeds and used for reduced-representation genome sequencing to identify SNP loci associated with pubertal initiation and litter size. Selective signatures analysis was performed based on the fixation index and nucleotide diversity, followed by pathway analysis of the genes contained in the selected regions. The selected SNP loci in the genes associated with pubertal initiation and litter size were validated using both sheep breeds. In total, 384,718 high quality SNPs were obtained and 376 genes were selected. Functional annotation of genes and enrichment analysis identified 12 genes associated with pubertal initiation and 11 genes associated with litter size. SNP locus validation showed that two SNP on PAK1 and four on ADCY1 may be associated with pubertal initiation, and one SNP on GNAQ gene (NC_040253.1: g.62677376G > A) was associated with litter size in Cele Black sheep. Our results provide new theoretical support for sheep breeding.

Expression and functional analysis of GnRH at the onset of puberty in sheep.

Gonadotropin-releasing hormone (GnRH) is a key factor at the onset of puberty. This decapeptide has been found in mammalian ovaries, but its regulatory mechanism in the ovary of sheep at the onset of puberty is not clear. This study investigated the coding sequence (CDS) of the GnRH gene in the ovary of Duolang sheep and the expression of GnRH mRNA in different tissues at the onset of puberty, and analyzed the effect of GnRH on ovarian granulosa cells (GCs) of Duolang sheep. The results showed that the GnRH CDS of sheep was cloned, the full length of the GnRH CDS in sheep ovary was 279 bp, and the nucleotide sequence was completely homologous to that in the hypothalamus. The expression of GnRH mRNA was highest in the hypothalamus and ovary. The expression of related hormones and receptors in GCs of Duolang sheep treated with different concentrations of GnRH for 24 h was affected. GnRH significantly inhibited LH synthesis and LHR expression in GCs. Low concentration (100 ng mL - 1 ) had the most obvious therapeutic effect on follicle-stimulating hormone (FSH) and FSHR. Higher concentration (250 ng mL - 1 ) significantly promoted estradiol and ER ß mRNA. These findings provide strong evidence that ovarian GnRH is an important regulatory factor at the onset of puberty in sheep.

Analysis on the desert adaptability of indigenous sheep in the southern edge of Taklimakan Desert.

The southern margin of the Taklimakan Desert is characterized by low rainfall, heavy sandstorms, sparse vegetation and harsh ecological environment. The indigenous sheep in this area are rich in resources, with the advantages of perennial estrus and good resistance to stress in most sheep. Exploring the molecular markers of livestock adaptability in this environment will provide the molecular basis for breeding research to cope with extreme future changes in the desert environment. In this study, we analyzed the population genetic structure and linkage imbalance of five sheep breeds with three different agricultural geographic characteristics using four complementary genomic selection signals: fixation index (FST), cross-population extended haplotype homozygosity (xp-EHH), Rsb (extended haplotype homozygosity between-populations) and iHS (integrated haplotype homozygosity score). We used Illumina Ovine SNP 50K Genotyping BeadChip Array, and gene annotation and enrichment analysis were performed on selected regions of the obtained genome. The ovary of Qira Black sheep (Follicular phase, Luteal phase, 30th day of pregnancy, 45th day of pregnancy) was collected, and the differentially expressed genes were screened by transcriptomic sequencing. Genome-wide selective sweep results and transcriptome data were combined for association analysis to obtain candidate genes associated with perennial estrus and stable reproduction. In order to verify the significance of the results, 15 resulting genes were randomly selected for fluorescence quantitative analysis. The results showed that Dolang sheep and Qira Black sheep evolved from Kazak sheep. Linkage disequilibrium analysis showed that the decay rate of sheep breeds in the Taklimakan Desert was higher than that in Yili grassland. The signals of FST, xp-EHH, Rsb and iHS detected 526, 332, 308 and 408 genes, respectively, under the threshold of 1% and 17 overlapping genes under the threshold of 5%. A total of 29 genes were detected in association analysis of whole-genome and transcriptome data. This study reveals the genetic mechanism of perennial estrus and environmental adaptability of indigenous sheep breeds in the Taklimakan Desert. It provides a theoretical basis for the conservation and exploitation of genetic resources of indigenous sheep breeds in extreme desert environment. This provides a new perspective for the quick adaptation of sheep and other mammals to extreme environments and future climate changes.

Transcriptome Sequencing-Based Mining of Genes Associated With Pubertal Initiation in Dolang Sheep.

Improving the fertility of sheep is an important goal in sheep breeding as it greatly increases the productivity. Dolang sheep is a typical representative breed of lamb in Xinjiang and is the main local sheep breed and meat source in the region. To explore the genes associated with the initiation of puberty in Dolang sheep, the hypothalamic tissues of Dolang sheep prepubertal, pubertal, and postpubertal periods were collected for RNA-seq analysis on the Illumina platform, generating 64.08 Gb clean reads. A total of 575, 166, and 648 differentially expressed genes (DEGs) were detected in prepuberty_vs._puberty, postpuberty_vs._prepuberty, and postpuberty_vs._puberty analyses, respectively. Based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, the related genes involved in the initiation of puberty in Dolang sheep were mined. Ten genes that have direct or indirect functions in the initiation of puberty in Dolang sheep were screened using the GO and KEGG results. Additionally, quantitative real-time PCR was used to verify the reliability of the RNA-Seq data. This study provided a new approach for revealing the mechanism of puberty initiation in sheep and provided a theoretical basis and candidate genes for the breeding of early-pubertal sheep by molecular techniques, and at the same time, it is also beneficial for the protection, development, and utilization of the fine genetic resources of Xinjiang local sheep.

Discovery and Molecular Basis of a Diverse Set of Polycomb Repressive Complex 2 Inhibitors Recognition by EED.

Polycomb repressive complex 2 (PRC2), a histone H3 lysine 27 methyltransferase, plays a key role in gene regulation and is a known epigenetics drug target for cancer therapy. The WD40 domain-containing protein EED is the regulatory subunit of PRC2. It binds to the tri-methylated lysine 27 of the histone H3 (H3K27me3), and through which stimulates the activity of PRC2 allosterically. Recently, we disclosed a novel PRC2 inhibitor EED226 which binds to the K27me3-pocket on EED and showed strong antitumor activity in xenograft mice model. Here, we further report the identification and validation of four other EED binders along with EED162, the parental compound of EED226. The crystal structures for all these five compounds in complex with EED revealed a common deep pocket induced by the binding of this diverse set of compounds. This pocket was created after significant conformational rearrangement of the aromatic cage residues (Y365, Y148 and F97) in the H3K27me3 binding pocket of EED, the width of which was delineated by the side chains of these rearranged residues. In addition, all five compounds interact with the Arg367 at the bottom of the pocket. Each compound also displays unique features in its interaction with EED, suggesting the dynamics of the H3K27me3 pocket in accommodating the binding of different compounds. Our results provide structural insights for rational design of novel EED binder for the inhibition of PRC2 complex activity.

Discovery and Characterization of Allosteric WNK Kinase Inhibitors.

Protein kinases are known for their highly conserved adenosine triphosphate (ATP)-binding site, rendering the discovery of selective inhibitors a major challenge. In theory, allosteric inhibitors can achieve high selectivity by targeting less conserved regions of the kinases, often with an added benefit of retaining efficacy under high physiological ATP concentration. Although often overlooked in favor of ATP-site directed approaches, performing a screen at high ATP concentration or stringent hit triaging with high ATP concentration offers conceptually simple methods of identifying inhibitors that bind outside the ATP pocket. Here, we applied the latter approach to the With-No-Lysine (K) (WNK) kinases to discover lead molecules for a next-generation antihypertensive that requires a stringent safety profile. This strategy yielded several ATP noncompetitive WNK1-4 kinase inhibitors, the optimization of which enabled cocrystallization with WNK1, revealing an allosteric binding mode consistent with the observed exquisite specificity for WNK1-4 kinases. The optimized compound inhibited rubidium uptake by sodium chloride cotransporter 1 (NKCC1) in HT29 cells, consistent with the reported physiology of WNK kinases in renal electrolyte handling.

Matrix-Based Activity Pattern Classification as a Novel Method for the Characterization of Enzyme Inhibitors Derived from High-Throughput Screening.

One of the central questions in the characterization of enzyme inhibitors is determining the mode of inhibition (MOI). Classically, this is done with a number of low-throughput methods in which inhibition models are fitted to the data. The ability to rapidly characterize the MOI for inhibitors arising from high-throughput screening in which hundreds to thousands of primary inhibitors may need to be characterized would greatly help in lead selection efforts. Here we describe a novel method for determining the MOI of a compound without the need for curve fitting of the enzyme inhibition data. We provide experimental data to demonstrate the utility of this new high-throughput MOI classification method based on nonparametric analysis of the activity derived from a small matrix of substrate and inhibitor concentrations (e.g., from a 4S × 4I matrix). Lists of inhibitors from four different enzyme assays are studied, and the results are compared with the previously described IC50-shift method for MOI classification. The MOI results from this method are in good agreement with the known MOI and compare favorably with those from the IC50-shift method. In addition, we discuss some advantages and limitations of the method and provide recommendations for utilization of this MOI classification method.

Allosteric inhibition of SHP2 phosphatase inhibits cancers driven by receptor tyrosine kinases.

The non-receptor protein tyrosine phosphatase SHP2, encoded by PTPN11, has an important role in signal transduction downstream of growth factor receptor signalling and was the first reported oncogenic tyrosine phosphatase. Activating mutations of SHP2 have been associated with developmental pathologies such as Noonan syndrome and are found in multiple cancer types, including leukaemia, lung and breast cancer and neuroblastoma. SHP2 is ubiquitously expressed and regulates cell survival and proliferation primarily through activation of the RAS­ERK signalling pathway. It is also a key mediator of the programmed cell death 1 (PD-1) and B- and T-lymphocyte attenuator (BTLA) immune checkpoint pathways. Reduction of SHP2 activity suppresses tumour cell growth and is a potential target of cancer therapy. Here we report the discovery of a highly potent (IC50 = 0.071 µM), selective and orally bioavailable small-molecule SHP2 inhibitor, SHP099, that stabilizes SHP2 in an auto-inhibited conformation. SHP099 concurrently binds to the interface of the N-terminal SH2, C-terminal SH2, and protein tyrosine phosphatase domains, thus inhibiting SHP2 activity through an allosteric mechanism. SHP099 suppresses RAS­ERK signalling to inhibit the proliferation of receptor-tyrosine-kinase-driven human cancer cells in vitro and is efficacious in mouse tumour xenograft models. Together, these data demonstrate that pharmacological inhibition of SHP2 is a valid therapeutic approach for the treatment of cancers.

Allosteric Inhibition of SHP2: Identification of a Potent, Selective, and Orally Efficacious Phosphatase Inhibitor.

SHP2 is a nonreceptor protein tyrosine phosphatase (PTP) encoded by the PTPN11 gene involved in cell growth and differentiation via the MAPK signaling pathway. SHP2 also purportedly plays an important role in the programmed cell death pathway (PD-1/PD-L1). Because it is an oncoprotein associated with multiple cancer-related diseases, as well as a potential immunomodulator, controlling SHP2 activity is of significant therapeutic interest. Recently in our laboratories, a small molecule inhibitor of SHP2 was identified as an allosteric modulator that stabilizes the autoinhibited conformation of SHP2. A high throughput screen was performed to identify progressable chemical matter, and X-ray crystallography revealed the location of binding in a previously undisclosed allosteric binding pocket. Structure-based drug design was employed to optimize for SHP2 inhibition, and several new protein-ligand interactions were characterized. These studies culminated in the discovery of 6-(4-amino-4-methylpiperidin-1-yl)-3-(2,3-dichlorophenyl)pyrazin-2-amine (SHP099, 1), a potent, selective, orally bioavailable, and efficacious SHP2 inhibitor.

Application of Titration-Based Screening for the Rapid Pilot Testing of High-Throughput Assays.

Pilot testing of an assay intended for high-throughput screening (HTS) with small compound sets is a necessary but often time-consuming step in the validation of an assay protocol. When the initial testing concentration is less than optimal, this can involve iterative testing at different concentrations to further evaluate the pilot outcome, which can be even more time-consuming. Quantitative HTS (qHTS) enables flexible and rapid collection of assay performance statistics, hits at different concentrations, and concentration-response curves in a single experiment. Here we describe the qHTS process for pilot testing in which eight-point concentration-response curves are produced using an interplate asymmetric dilution protocol in which the first four concentrations are used to represent the range of typical HTS screening concentrations and the last four concentrations are added for robust curve fitting to determine potency/efficacy values. We also describe how these data can be analyzed to predict the frequency of false-positives, false-negatives, hit rates, and confirmation rates for the HTS process as a function of screening concentration. By taking into account the compound pharmacology, this pilot-testing paradigm enables rapid assessment of the assay performance and choosing the optimal concentration for the large-scale HTS in one experiment.

Selective inhibition of Ezh2 by a small molecule inhibitor blocks tumor cells proliferation.

Ezh2 (Enhancer of zeste homolog 2) protein is the enzymatic component of the Polycomb repressive complex 2 (PRC2), which represses gene expression by methylating lysine 27 of histone H3 (H3K27) and regulates cell proliferation and differentiation during embryonic development. Recently, hot-spot mutations of Ezh2 were identified in diffused large B-cell lymphomas and follicular lymphomas. To investigate if tumor growth is dependent on the enzymatic activity of Ezh2, we developed a potent and selective small molecule inhibitor, EI1, which inhibits the enzymatic activity of Ezh2 through direct binding to the enzyme and competing with the methyl group donor S-Adenosyl methionine. EI1-treated cells exhibit genome-wide loss of H3K27 methylation and activation of PRC2 target genes. Furthermore, inhibition of Ezh2 by EI1 in diffused large B-cell lymphomas cells carrying the Y641 mutations results in decreased proliferation, cell cycle arrest, and apoptosis. These results provide strong validation of Ezh2 as a potential therapeutic target for the treatment of cancer.

A sensitive, homogeneous, and high-throughput assay for lysine-specific histone demethylases at the H3K4 site.

Histone methylation is a regulated feature of nucleosomes that can have an impact on gene expression. The methylation state of histone residues has also been found in recent years to be associated with various disorders. Tools for detecting methylation state changes are very useful for dissecting the function of these epigenetic marks. In this work, a sensitive homogeneous assay for histone demethylase activity at the H3K4 site has been developed in a time-resolved fluorescent resonance energy transfer assay format. The assay is based on the detection of the unmethylated H3 peptide by a fluorescent europium-chelate labeled monoclonal antibody binding specifically to the H3K4 site. The assay was validated for histone lysine-specific demethylase 1 and was demonstrated to be a suitable assay for inhibitor profiling and high-throughput screening.

3,5-diarylazoles as novel and selective inhibitors of protein kinase D.

The synthesis and preliminary studies of the SAR of novel 3,5-diarylazole inhibitors of Protein Kinase D (PKD) are reported. Notably, optimized compounds in this class have been found to be active in cellular assays of phosphorylation-dependant HDAC5 nuclear export, orally bioavailable, and highly selective versus a panel of additional putative histone deacetylase (HDAC) kinases. Therefore these compounds could provide attractive tools for the further study of PKD/HDAC5 signaling.

Compound set enrichment: a novel approach to analysis of primary HTS data.

The main goal of high-throughput screening (HTS) is to identify active chemical series rather than just individual active compounds. In light of this goal, a new method (called compound set enrichment) to identify active chemical series from primary screening data is proposed. The method employs the scaffold tree compound classification in conjunction with the Kolmogorov-Smirnov statistic to assess the overall activity of a compound scaffold. The application of this method to seven PubChem data sets (containing between 9389 and 263679 molecules) is presented, and the ability of this method to identify compound classes with only weakly active compounds (potentially latent hits) is demonstrated. The analysis presented here shows how methods based on an activity cutoff can distort activity information, leading to the incorrect activity assignment of compound series. These results suggest that this method might have utility in the rational selection of active classes of compounds (and not just individual active compounds) for followup and validation.

Identification of orally available naphthyridine protein kinase D inhibitors.

A novel 2,6-naphthyridine was identified by high throughput screen (HTS) as a dual protein kinase C/D (PKC/PKD) inhibitor. PKD inhibition in the heart was proposed as a potential antihypertrophic mechanism with application as a heart failure therapy. As PKC was previously identified as the immediate upstream activator of PKD, PKD vs PKC selectivity was essential to understand the effect of PKD inhibition in models of cardiac hypertrophy and heart failure. The present study describes the modification of the HTS hit to a series of prototype pan-PKD inhibitors with routine 1000-fold PKD vs PKC selectivity. Example compounds inhibited PKD activity in vitro, in cells, and in vivo following oral administration. Their effects on heart morphology and function are discussed herein.