Identification of gene sets and pathways associated with lactation performance in mice

Single-Cell Transcription Mapping of Murine and Human Mammary Organoids Responses to Female Hormones

During female adolescence and pregnancy, rising levels of hormones result in a cyclic source of signals that control the development of mammary tissue. While such alterations are well understood from a whole-gland perspective, the alterations that such hormones bring to organoid cultures derived from mammary glands have yet to be fully mapped. This is of special importance given that organoids are considered suitable systems to understand cross species breast development. Here we utilized single-cell transcriptional profiling to delineate responses of murine and human normal breast organoid systems to female hormones across evolutionary distinct species. Collectively, our study represents a molecular atlas of epithelial dynamics in response to estrogen and pregnancy hormones.

Single-cell transcription mapping of murine and human mammary organoids responses to female hormones

During female adolescence and pregnancy, rising levels of hormones result in a cyclic source of signals that control the development of mammary tissue. While such alterations are well understood from a whole-gland perspective, the alterations that such hormones bring to organoid cultures derived from mammary glands have yet to be fully mapped. This is of special importance given that organoids are considered suitable systems to understand cross species breast development. Here we utilized single-cell transcriptional profiling to delineate responses of murine and human normal breast organoid systems to female hormones across evolutionary distinct species. Collectively, our study represents a molecular atlas of epithelial dynamics in response to estrogen and pregnancy hormones.

Fur removal promotes an earlier expression of involution-related genes in mammary gland of lactating mice

Peak lactation occurs when milk production is at its highest. The factors limiting peak lactation performance have been subject of intense debate. Milk production at peak lactation appears limited by the capacity of lactating females to dissipate body heat generated as a by-product of processing food and producing milk. As a result, manipulations that enhance capacity to dissipate body heat (such as fur removal) increase peak milk production. We investigated the potential correlates of shaving-induced increases in peak milk production in laboratory mice. By transcriptomic profiling of the mammary gland, we searched for the mechanisms underlying experimentally increased milk production and its consequences for mother-young conflict over weaning, manifested by advanced or delayed involution of mammary gland. We demonstrated that shaving-induced increases in milk production were paradoxically linked to reduced expression of some milk synthesis-related genes. Moreover, the mammary glands of shaved mice had a gene expression profile indicative of earlier involution relative to unshaved mice. Once provided with enhanced capacity to dissipate body heat, shaved mice were likely to rear their young to independence faster than unshaved mothers.

Comparative transcriptomic analysis for identification of candidate sex-related genes and pathways in Crimson seabream (Parargyrops edita)

Teleost fishes display the largest array of sex-determining systems among animals, resulting in various reproductive strategies. Research on sex-related genes in teleosts will broaden our understanding of the process, and provide important insight into the plasticity of the sex determination process in vertebrates in general. Crimson seabream (Parargyrops edita Tanaka, 1916) is one of the most valuable and abundant fish resources throughout Asia. However, little genomic information on P. edita is available. In the present study, the transcriptomes of male and female P. edita were sequenced with RNA-seq technology. A total of 388,683,472 reads were generated from the libraries. After filtering and assembling, a total of 79,775 non redundant unigenes were obtained with an N50 of 2,921 bp. The unigenes were annotated with multiple public databases, including NT (53,556, 67.13%), NR (54,092, 67.81%), Swiss-Prot (45,265, 56.74%), KOG (41,274, 51.74%), KEGG (46,302, 58.04%), and GO (11,056, 13.86%) databases. Comparison of the unigenes of different sexes of P. edita revealed that 11,676 unigenes (9,335 in females, 2,341 in males) were differentially expressed between males and females. Of these, 5,463 were specifically expressed in females, and 1,134 were specifically expressed in males. In addition, the expression levels of ten unigenes were confirmed to validate the transcriptomic data by qRT-PCR. Moreover, 34,473 simple sequence repeats (SSRs) were identified in SSR-containing sequences, and 50 loci were randomly selected for primer development. Of these, 36 loci were successfully amplified, and 19 loci were polymorphic. Finally, our comparative analysis identified many sex-related genes (zps, amh, gsdf, sox4, cyp19a, etc.) and pathways (MAPK signaling pathway, p53 signaling pathway, etc.) of P. edita. This informative transcriptomic analysis provides valuable data to increase genomic resources of P. edita. The results will be useful for clarifying the molecular mechanism of sex determination and for future functional analyses of sex-associated genes.

GWAS-Based Identification of New Loci for Milk Yield, Fat, and Protein in Holstein Cattle

High-yield and high-quality of milk are the primary goals of dairy production. Understanding the genetic architecture underlying these milk-related traits is beneficial so that genetic variants can be targeted toward the genetic improvement. In this study, we measured five milk production and quality traits in Holstein cattle population from China. These traits included milk yield, fat, and protein. We used the estimated breeding values as dependent variables to conduct the genome-wide association studies (GWAS). Breeding values were estimated through pedigree relationships by using a linear mixed model. Genotyping was carried out on the individuals with phenotypes by using the Illumina BovineSNP150 BeadChip. The association analyses were conducted by using the fixed and random model Circulating Probability Unification (FarmCPU) method. A total of ten single-nucleotide polymorphisms (SNPs) were detected above the genome-wide significant threshold (p < 4.0 × 10-7), including six located in previously reported quantitative traits locus (QTL) regions. We found eight candidate genes within distances of 120 kb upstream or downstream to the associated SNPs. The study not only identified the effect of DGAT1 gene on milk fat and protein, but also discovered novel genetic loci and candidate genes related to milk traits. These novel genetic loci would be an important basis for molecular breeding in dairy cattle.

Genome-wide association and pathway-based analysis using latent variables related to milk protein composition and cheesemaking traits in dairy cattle

The aim of this study was to perform genome-wide associations (GWAS) and gene-set enrichment analyses with protein composition and cheesemaking-related latent variables (factors; F) in a cohort of 1,011 Italian Brown Swiss cows. Factor analysis was applied to identify latent structures of 26 phenotypes related to bovine milk quantity and quality, protein fractions [α_S1-, α_S2-, β-, and κ-casein (CN), β-lactoglobulin, and α-lactalbumin (α-LA)], coagulation and curd firming at time t (CF_t) measures, and cheese properties [cheese yield (%CY) and nutrients recovery in the curd] of individual cows. Ten orthogonal F were extracted, explaining 74% of the original variability. Factor 1_%CY underlined the %CY characteristics, F2_CFt was related to the CF_t process parameters, F3_Yield was considered as descriptor of milk and solids yield, whereas F4_{Cheese N} underscored the presence of nitrogenous compounds (N) into the cheese. Four more F were related to the milk caseins (F5_αS1-β-CN, F7_β-κ-CN, F8_αS2-CN, and F9_αS1-CN-Ph) and 1 F was linked to the whey protein (F10_α-LA); 1 F underlined the udder health status (F6_{Udder health}). All cows were genotyped with the Illumina BovineSNP50 Bead Chip v.2 (Illumina Inc., San Diego, CA). Single marker regression GWAS were fitted. Gene-set enrichment analysis was run on GWAS results, using the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway databases, to reveal ontologies or pathways associated with the F. All F but F3_Yield showed significance in GWAS. Signals in 10 Bos taurus autosomes (BTA) were detected. High peaks on BTA6 (∼87 Mbp) were found for F6_β-κ-CN, F5_αS1-β-CN, and at the tail of BTA11 (∼104 Mbp) for F4_{Cheese N}. Gene-set enrichment analyses showed significant results (false discovery rate at 5%) for F8_αS2-CN, F1_%CY, F4_{Cheese N}, and F10_α-LA. For F8_αS2-CN, 33 Gene Ontology terms and 3 Kyoto Encyclopedia of Genes and Genomes categories were enriched, including terms related to ion transport and homeostasis, neuron function or part, and GnRH signaling pathway. Our results support the feasibility of factor analysis as a dimension reduction technique in genomic studies and evidenced a potential key role of α_S2-CN in milk quality and composition.

De novo transcriptome sequencing and analysis of male, pseudo-male and female yellow perch, Perca flavescens

Transcriptome sequencing could facilitate discovery of sex-biased genes, biological pathways and molecular markers, which could help clarify the molecular mechanism of sex determination and sexual dimorphism, and assist with selective breeding in aquaculture. Yellow perch has unique gonad system and sexual dimorphism and is an alternative model to study mechanism of sex determination, sexual dimorphism and sexual selection. In this study, we performed the de novo assembly of yellow perch gonads and muscle transcriptomes by high throughput Illumina sequencing. A total of 212,180 contigs were obtained, ranging from 127 to 64,876 bp, and N50 of 1,066 bp. The assembly RNA-Seq contigs (≥200bp) were then used for subsequent analyses, including annotation, pathway analysis, and microsatellites discovery. No female- and pseudo-male-biased genes were involved in any pathways while male-biased genes were involved in 29 pathways, and neuroactive ligand receptor interaction and enzyme of trypsin (enzyme code, EC: 3.4.21.4) was highly involved. Pyruvate kinase (enzyme code, EC: 2.7.1.40), which plays important roles in cell proliferation, was highly expressed in muscles. In addition, a total of 183,939 SNPs, 11,286 InDels and 41,479 microsatellites were identified. This study is the first report on transcriptome information in Percids, and provides rich resources for conducting further studies on understanding the molecular basis of sex determinations, sexual dimorphism, and sexual selection in fish, and for population studies and marker-assisted selection in Percids.

Pathway-based genome-wide association analysis of milk coagulation properties, curd firmness, cheese yield, and curd nutrient recovery in dairy cattle

It is becoming common to complement genome-wide association studies (GWAS) with gene-set enrichment analysis to deepen the understanding of the biological pathways affecting quantitative traits. Our objective was to conduct a gene ontology and pathway-based analysis to identify possible biological mechanisms involved in the regulation of bovine milk technological traits: coagulation properties, curd firmness modeling, individual cheese yield (CY), and milk nutrient recovery into the curd (REC) or whey loss traits. Results from 2 previous GWAS studies using 1,011 cows genotyped for 50k single nucleotide polymorphisms were used. Overall, the phenotypes analyzed consisted of 3 traditional milk coagulation property measures [RCT: rennet coagulation time defined as the time (min) from addition of enzyme to the beginning of coagulation; k₂₀: the interval (min) from RCT to the time at which a curd firmness of 20 mm is attained; a₃₀: a measure of the extent of curd firmness (mm) 30 min after coagulant addition], 6 curd firmness modeling traits [RCT_eq: RCT estimated through the CF equation (min); CF_P: potential asymptotic curd firmness (mm); k_CF: curd-firming rate constant (% × min^-1); k_SR: syneresis rate constant (% × min^-1); CF_max: maximum curd firmness (mm); and t_max: time to CF_max (min)], 3 individual CY-related traits expressing the weight of fresh curd (%CY_CURD), curd solids (%CY_SOLIDS), and curd moisture (%CY_WATER) as a percentage of weight of milk processed and 4 milk nutrient and energy recoveries in the curd (REC_FAT, REC_PROTEIN, REC_SOLIDS, and REC_ENERGY calculated as the % ratio between the nutrient in curd and the corresponding nutrient in processed milk), milk pH, and protein percentage. Each trait was analyzed separately. In total, 13,269 annotated genes were used in the analysis. The Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway databases were queried for enrichment analyses. Overall, 21 Gene Ontology and 17 Kyoto Encyclopedia of Genes and Genomes categories were significantly associated (false discovery rate at 5%) with 7 traits (RCT, RCT_eq, k_CF, %CY_SOLIDS, REC_FAT, REC_SOLIDS, and REC_ENERGY), with some being in common between traits. The significantly enriched categories included calcium signaling pathway, salivary secretion, metabolic pathways, carbohydrate digestion and absorption, the tight junction and the phosphatidylinositol pathways, as well as pathways related to the bovine mammary gland health status, and contained a total of 150 genes spanning all chromosomes but 9, 20, and 27. This study provided new insights into the regulation of bovine milk coagulation and cheese ability that were not captured by the GWAS.

Anesthetic Propofol-Induced Gene Expression Changes in Patients Undergoing Coronary Artery Bypass Graft Surgery Based on Dynamical Differential Coexpression Network Analysis

We aimed to determine the influence of anesthetic propofol on gene expression in patients treated by coronary artery bypass graft (CABG) surgery based on differential coexpression network (DCN) and to further reveal the novel mechanisms of the cardioprotective effects of propofol. Firstly, we constructed the DCN for disease condition based on Pearson correlation coefficient (PCC) and weight value. Secondly, the inference of modules was applied to search modules from DCN with same members but varied connectivity. Furthermore, we measured the statistical significance of the modules for selecting differential modules (DMs). Finally, attract method was used for DMs analysis to select key modules. Based on the δ value, 11928 edges and 2956 nodes were chosen to construct DCNs. A total of 29 seed genes were selected. Moreover, by quantifying connectivity changes in shared gene modules across different conditions, 8 DMs with higher connectivity dynamics were identified. Then, we extracted key modules using attract method, there were 8 key modules, and the top 3 modules were module 1, 2, and 3. Furthermore, GCG, PPY, and PON1 were initial seed genes of these 3 key modules, respectively. Accordingly, GCG and PON1 might exert important roles in the cardioprotective effects of propofol during CABG.

Targeted Analysis Reveals an Important Role of JAK-STAT-SOCS Genes for Milk Production Traits in Australian Dairy Cattle

The Janus kinase and signal transducer and activator of transcription (JAK-STAT) pathway genes along with suppressors of cytokine signalling (SOCS) family genes play a crucial role in controlling cytokine signals in the mammary gland and thus mammary gland development. Mammary gene expression studies showed differential expression patterns for all the JAK-STAT pathway genes. Gene expression studies using qRT-PCR revealed differential expression of SOCS2, SOCS4, and SOCS5 genes across the lactation cycle in dairy cows. Using genotypes from 1,546 Australian Holstein-Friesian bulls, a statistical model for an association analysis based on SNPs within 500 kb of JAK-STAT pathway genes, and SOCS genes alone was constructed. The analysis suggested that these genes and pathways make a significant contribution to the Australian milk production traits. There were 24 SNPs close to SOCS1, SOCS3, SOCS5, SOCS7, and CISH genes that were significantly associated with Australian Profit Ranking (APR), Australian Selection Index (ASI), and protein yield (PY). This study supports the view that there may be some merit in choosing SNPs around functionally relevant genes for the selection and genetic improvement schemes for dairy production traits.

An integrated genomics approach to identify genetic regions associated with neonatal growth trait in mice

Neonatal growth during the early post-partum period is closely associated with lactation performance. Neonatal growth reflects milk output and is a complex variable trait among inbred mouse strains, but few studies have compared this trait systematically across more than a few strains. In the present study, 11 inbred strains of mice were measured for a neonatal growth phenotype during the first eight days of lactation. Significant differences in neonatal growth trait were observed with QSi5 (3.71±0.05 g) and DBA/1J (2.67±0.06 g) strains defining the two extremes of the phenotype. In silico association analysis was performed for trait variability using the high density SNP information on inbred strains of mice. We found strong evidence to refine a previously identified large neonatal growth QTL on mouse chromosome 9, Neogq1. When an integrated strategy that combined fine mapping and analysis of mammary transcriptome expression profiles of lactating mice with divergent phenotypes was applied, we identified neogenin (Neo1), a gene important for mammary gland morphogenesis, as a likely quantitative trait gene (QTG) underlying the Neogq1 QTL in mice.

DanQi Pill protects against heart failure through the arachidonic acid metabolism pathway by attenuating different cyclooxygenases and leukotrienes B4

Background

Chinese herbal formulae are composed of complex components and produce comprehensive pharmacological effects. Unlike chemical drugs that have only one clear single target, the components of Chinese herbal formulae have multiple channels and targets. How to discover the pharmacological targets of Chinese herbal formulae and their underlying molecular mechanism are still under investigation.

Methods

DanQi pill (DQP), which is one of the widely prescribed traditional Chinese medicines, is applied as an example drug. In this study, we used the drug target prediction model (DrugCIPHER-CS) to examine the underlying molecular mechanism of DQP, followed by experimental validation.

Results

A novel therapeutic effect pattern of DQP was identified. After determining the compounds in DQP, we used DrugCIPHER-CS to predict their potential targets. These potential targets were significantly enriched in well-known cardiovascular disease-related pathways. For example, the biological processes of neuroactive ligand–receptor interaction, calcium-signaling pathway, and aminoacyl–tRNA biosynthesis were involved. A new and significant pathway, arachidonic acid (AA) metabolism, was also identified in this study. This predicted pathway alteration was validated with an animal model of heart failure (HF). Results show that DQP had effect both on thromboxane B2 (TXB2) and Prostaglandin I2 (PGI2) in different patterns. It can down-regulate the TXB2 and up-regulate the PGI2 in diverse way. Remarkably, it also had effect on cyclooxygenase (COX)-1 and COX2 by suppressing their levels, which may be the critical and novel mechanism of cardiacprotective efficacy for DQP. Furthermore, leukotrienes B4 (LTB4) receptor, another key molecule of AA metabolism which finally mediated gastrotoxic leukotrienes, was also reduced by DQP.

Conclusions

The combination of drug target prediction and experimental validation provides new insights into the complicated mechanism of DQP.

Sequencing the transcriptome of milk production: milk trumps mammary tissue

BACKGROUND

Studies of normal human mammary gland development and function have mostly relied on cell culture, limited surgical specimens, and rodent models. Although RNA extracted from human milk has been used to assay the mammary transcriptome non-invasively, this assay has not been adequately validated in primates. Thus, the objectives of the current study were to assess the suitability of lactating rhesus macaques as a model for lactating humans and to determine whether RNA extracted from milk fractions is representative of RNA extracted from mammary tissue for the purpose of studying the transcriptome of milk-producing cells.

RESULTS

We confirmed that macaque milk contains cytoplasmic crescents and that ample high-quality RNA can be obtained for sequencing. Using RNA sequencing, RNA extracted from macaque milk fat and milk cell fractions more accurately represented RNA from mammary epithelial cells (cells that produce milk) than did RNA from whole mammary tissue. Mammary epithelium-specific transcripts were more abundant in macaque milk fat, whereas adipose or stroma-specific transcripts were more abundant in mammary tissue. Functional analyses confirmed the validity of milk as a source of RNA from milk-producing mammary epithelial cells.

CONCLUSIONS

RNA extracted from the milk fat during lactation accurately portrayed the RNA profile of milk-producing mammary epithelial cells in a non-human primate. However, this sample type clearly requires protocols that minimize RNA degradation. Overall, we validated the use of RNA extracted from human and macaque milk and provided evidence to support the use of lactating macaques as a model for human lactation.