Search

Abstract The infant gut microbiome is increasingly recognized as a reservoir of antibiotic resistance genes, yet the assembly of gut resistome in infants and its influencing factors remain largely unknown. We characterized resistome in 4132 metagenomes from 963 infants in six countries and 4285 resistance genes were observed. The inherent resistome pattern of healthy infants (N = 272) could be distinguished by two stages: a multicompound resistance phase (Months 0–7) and a tetracycline‐mupirocin‐β‐lactam‐dominant phase (Months 8–14). Microbial taxonomy explained 40.7% of the gut resistome of healthy infants, with Escherichia (25.5%) harboring the most resistance genes. In a further analysis with all available infants (N = 963), we found age was the strongest influencer on the resistome and was negatively correlated with the overall resistance during the first 3 years (p

Abstract The carbohydrate fraction of most mammalian milks contains a variety of oligosaccharides that encompass a range of structures and monosaccharide compositions. Human milk oligosaccharides have received considerable attention due to their biological roles in neonatal gut microbiota, immunomodulation, and brain development. However, a major challenge in understanding the biology of milk oligosaccharides across other mammals is that reports span more than 5 decades of publications with varying data reporting methods. In the present study, publications on milk oligosaccharide profiles were identified and harmonized into a standardized format to create a comprehensive, machine-readable database of milk oligosaccharides across mammalian species. The resulting database, MilkOligoDB, includes 3193 entries for 783 unique oligosaccharide structures from the milk of 77 different species harvested from 113 publications. Cross-species and cross-publication comparisons of milk oligosaccharide profiles reveal common structural motifs within mammalian orders. Of the species studied, only chimpanzees, bonobos, and Asian elephants share the specific combination of fucosylation, sialylation, and core structures that are characteristic of human milk oligosaccharides. However, agriculturally important species do produce diverse oligosaccharides that may be valuable for human supplementation. Overall, MilkOligoDB facilitates cross-species and cross-publication comparisons of milk oligosaccharide profiles and the generation of new data-driven hypotheses for future research.

ABSTRACT Milk oligosaccharides (MOs) can be prebiotic and antiadhesive, while fatty acids (MFAs) can be antimicrobial. Both have been associated with milk microbes or mammary gland inflammation in humans. Relationships between these milk components and milk microbes or inflammation have not been determined for cows and could help elucidate a novel approach for the dairy industry to promote desired milk microbial composition for improvement of milk quality and reduction of milk waste. We aimed to determine relationships among milk microbiota, MFAs, MOs, lactose, and somatic cell counts (SCC) from Holstein cows, using our previously published data. Raw milk samples were collected at three time points, ranging from early to late lactation. Data were analyzed using linear mixed-effects modeling and repeated-measures correlation. Unsaturated MFA and short-chain MFA had mostly negative relationships with potentially pathogenic genera, including Corynebacterium, Pseudomonas, and an unknown Enterobacteriaceae genus but numerous positive relationships with symbionts Bifidobacterium and Bacteroides. Conversely, many MOs were positively correlated with potentially pathogenic genera (e.g., Corynebacterium, Enterococcus, and Pseudomonas), and numerous MOs were negatively correlated with the symbiont Bifidobacterium. The neutral, nonfucosylated MO composed of eight hexoses had a positive relationship with SCC, while lactose had a negative relationship with SCC. One interpretation of these trends might be that in milk, MFAs disrupt primarily pathogenic bacterial cells, causing a relative increase in abundance of beneficial microbial taxa, while MOs respond to and act on pathogenic taxa primarily through antiadhesive methods. Further research is needed to confirm the potential mechanisms driving these correlations. IMPORTANCE Bovine milk can harbor microbes that cause mastitis, milk spoilage, and foodborne illness. Fatty acids found in milk can be antimicrobial and milk oligosaccharides can have antiadhesive, prebiotic, and immune-modulatory effects. Relationships among milk microbes, fatty acids, oligosaccharides, and inflammation have been reported for humans. To our knowledge, associations among the milk microbial composition, fatty acids, oligosaccharides, and lactose have not been reported for healthy lactating cows. Identifying these potential relationships in bovine milk will inform future efforts to characterize direct and indirect interactions of the milk components with the milk microbiota. Since many milk components are associated with herd management practices, determining if these milk components impact milk microbes may provide valuable information for dairy cow management and breeding practices aimed at minimizing harmful and spoilage-causing microbes in raw milk.

Photo-based dietary assessment is becoming more feasible as artificial intelligence methods improve. However, advancement of these methods for dietary assessment in research settings has been hindered by the lack of an appropriate dataset against which to benchmark algorithm performance. We conducted the Surveying Nutrient Assessment with Photographs of Meals (SNAPMe) study (ClinicalTrials ID: NCT05008653) to pair meal photographs with traditional food records. Participants were recruited nationally, and 110 enrollment meetings were completed via web-based video conferencing. Participants uploaded and annotated their meal photos using a mobile phone app called Bitesnap and completed food records using the Automated Self-Administered 24-h Dietary Assessment Tool (ASA24®) version 2020. Participants included photos before and after eating non-packaged and multi-serving packaged meals, as well as photos of the front and ingredient labels for single-serving packaged foods. The SNAPMe Database (DB) contains 3311 unique food photos linked with 275 ASA24 food records from 95 participants who photographed all foods consumed and recorded food records in parallel for up to 3 study days each. The use of the SNAPMe DB to evaluate ingredient prediction demonstrated that the publicly available algorithms FB Inverse Cooking and Im2Recipe performed poorly, especially for single-ingredient foods and beverages. Correlations between nutrient estimates common to the Bitesnap and ASA24 dietary assessment tools indicated a range in predictive capacity across nutrients (cholesterol, adjusted R2 = 0.85, p < 0.0001; food folate, adjusted R2 = 0.21, p < 0.05). SNAPMe DB is a publicly available benchmark for photo-based dietary assessment in nutrition research. Its demonstrated utility suggested areas of needed improvement, especially the prediction of single-ingredient foods and beverages.

Dairy products are a good source of essential nutrients and past reviews have shown associations of dairy consumption with decreased systemic inflammation. Links between dairy intake and gastrointestinal (GI) inflammation are under-investigated. Therefore, we examined associations between reported dairy intake and markers of GI inflammation in healthy adults in a cross-sectional observational study, hypothesizing a negative association with yogurt intake, suggesting a protective effect, and no associations with total dairy, fluid milk, and cheese intake. Participants completed 24-h dietary recalls and a food frequency questionnaire (FFQ) to assess recent and habitual intake, respectively. Those who also provided a stool sample (n = 295), and plasma sample (n = 348) were included in analysis. Inflammation markers from stool, including calprotectin, neopterin, and myeloperoxidase, were measured along with LPS-binding protein (LBP) from plasma. Regression models tested associations between dairy intake variables and inflammation markers with covariates: age, sex, and body mass index (BMI). As yogurt is episodically consumed, we examined differences in inflammation levels between consumers (>0 cup equivalents/day reported in recalls) and non-consumers. We found no significant associations between dairy intake and markers of GI inflammation. In this cohort of healthy adults, dairy intake was not associated with GI inflammation.

ABSTRACT Antimicrobial resistance (AMR) represents a significant source of morbidity and mortality worldwide, with expectations that AMR-associated consequences will continue to worsen throughout the coming decades. Since resistance to antibiotics is encoded in the microbiome, interventions aimed at altering the taxonomic composition of the gut might allow us to prophylactically engineer microbiomes that harbor fewer antibiotic resistant genes (ARGs). Diet is one method of intervention, and yet little is known about the association between diet and antimicrobial resistance. To address this knowledge gap, we examined diet using the food frequency questionnaire (FFQ; habitual diet) and 24-h dietary recalls (Automated Self-Administered 24-h [ASA24®] tool) coupled with an analysis of the microbiome using shotgun metagenome sequencing in 290 healthy adult participants of the United States Department of Agriculture (USDA) Nutritional Phenotyping Study. We found that aminoglycosides were the most abundant and prevalent mechanism of AMR in these healthy adults and that aminoglycoside-O-phosphotransferases (aph3-dprime) correlated negatively with total calories and soluble fiber intake. Individuals in the lowest quartile of ARGs (low-ARG) consumed significantly more fiber in their diets than medium- and high-ARG individuals, which was concomitant with increased abundances of obligate anaerobes, especially from the family Clostridiaceae, in their gut microbiota. Finally, we applied machine learning to examine 387 dietary, physiological, and lifestyle features for associations with antimicrobial resistance, finding that increased phylogenetic diversity of diet was associated with low-ARG individuals. These data suggest diet may be a potential method for reducing the burden of AMR. IMPORTANCE Antimicrobial resistance (AMR) represents a considerable burden to health care systems, with the public health community largely in consensus that AMR will be a major cause of death worldwide in the coming decades. Humans carry antibiotic resistance in the microbes that live in and on us, collectively known as the human microbiome. Diet is a powerful method for shaping the human gut microbiome and may be a tractable method for lessening antibiotic resistance, and yet little is known about the relationship between diet and AMR. We examined this relationship in healthy individuals who contained various abundances of antibiotic resistance genes and found that individuals who consumed diverse diets that were high in fiber and low in animal protein had fewer antibiotic resistance genes. Dietary interventions may be useful for lessening the burden of antimicrobial resistance and might ultimately motivate dietary guidelines which will consider how nutrition can reduce the impact of infectious disease.

Abstract Background It has been estimated that at least 3% of the USA population consumes unpasteurized (raw) milk from animal sources, and the demand to legalize raw milk sales continues to increase. However, consumption of raw milk can cause foodborne illness and be a source of bacteria containing transferrable antimicrobial resistance genes (ARGs). To obtain a comprehensive understanding of the microbiome and antibiotic resistome in both raw and processed milk, we systematically analyzed 2034 retail milk samples including unpasteurized milk and pasteurized milk via vat pasteurization, high-temperature-short-time pasteurization, and ultra-pasteurization from the United States using complementary culture-based, 16S rRNA gene, and metagenomic sequencing techniques. Results Raw milk samples had the highest prevalence of viable bacteria which were measured as all aerobic bacteria, coliform, and Escherichia coli counts, and their microbiota was distinct from other types of milk. 16S rRNA gene sequencing revealed that Pseudomonadaceae dominated raw milk with limited levels of lactic acid bacteria. Among all milk samples, the microbiota remained stable with constant bacterial populations when stored at 4 °C. In contrast, storage at room temperature dramatically enriched the bacterial populations present in raw milk samples and, in parallel, significantly increased the richness and abundance of ARGs. Metagenomic sequencing indicated raw milk possessed dramatically more ARGs than pasteurized milk, and a conjugation assay documented the active transfer of bla CMY-2, one ceftazidime resistance gene present in raw milk-borne E. coli, across bacterial species. The room temperature-enriched resistome differed in raw milk from distinct geographic locations, a difference likely associated with regionally distinct milk microbiota. Conclusion Despite advertised “probiotic” effects, our results indicate that raw milk microbiota has minimal lactic acid bacteria. In addition, retail raw milk serves as a reservoir of ARGs, populations of which are readily amplified by spontaneous fermentation. There is an increased need to understand potential food safety risks from improper transportation and storage of raw milk with regard to ARGs. Video Abstract

Abstract Background Shotgun metagenomes are often assembled prior to annotation of genes which biases the functional capacity of a community towards its most abundant members. For an unbiased assessment of community function, short reads need to be mapped directly to a gene or protein database. The ability to detect genes in short read sequences is dependent on pre- and post-sequencing decisions. The objective of the current study was to determine how library size selection, read length and format, protein database, e-value threshold, and sequencing depth impact gene-centric analysis of human fecal microbiomes when using DIAMOND, an alignment tool that is up to 20,000 times faster than BLASTX. Results Using metagenomes simulated from a database of experimentally verified protein sequences, we find that read length, e-value threshold, and the choice of protein database dramatically impact detection of a known target, with best performance achieved with longer reads, stricter e-value thresholds, and a custom database. Using publicly available metagenomes, we evaluated library size selection, paired end read strategy, and sequencing depth. Longer read lengths were acheivable by merging paired ends when the sequencing library was size-selected to enable overlaps. When paired ends could not be merged, a congruent strategy in which both ends are independently mapped was acceptable. Sequencing depths of 5 million merged reads minimized the error of abundance estimates of specific target genes, including an antimicrobial resistance gene. Conclusions Shotgun metagenomes of DNA extracted from human fecal samples sequenced using the Illumina platform should be size-selected to enable merging of paired end reads and should be sequenced in the PE150 format with a minimum sequencing depth of 5 million merge-able reads to enable detection of specific target genes. Expecting the merged reads to be 180-250 bp in length, the appropriate e-value threshold for DIAMOND would then need to be more strict than the default. Accurate and interpretable results for specific hypotheses will be best obtained using small databases customized for the research question.

Livestock play a role in maintaining antibiotic resistance reservoirs. Here, Liu et al. use metagenomics to study the dynamics of antibiotic resistance genes (ARGs) in the gut of dairy cattle during early life, showing that colostrum is a potential source of ARGs and ARG abundance declines during nursing covarying with dietary transition.

Abstract Background Idiopathic chronic diarrhea (ICD) is a common cause of morbidity and mortality among juvenile rhesus macaques. Characterized by chronic inflammation of the colon and repeated bouts of diarrhea, ICD is largely unresponsive to medical interventions, including corticosteroid, antiparasitic, and antibiotic treatments. Although ICD is accompanied by large disruptions in the composition of the commensal gut microbiome, no single pathogen has been concretely identified as responsible for the onset and continuation of the disease. Results Fecal samples were collected from 12 ICD-diagnosed macaques and 12 age- and sex-matched controls. RNA was extracted for metatranscriptomic analysis of organisms and functional annotations associated with the gut microbiome. Bacterial, fungal, archaeal, protozoan, and macaque (host) transcripts were simultaneously assessed. ICD-afflicted animals were characterized by increased expression of host-derived genes involved in inflammation and increased transcripts from bacterial pathogens such as Campylobacter and Helicobacter and the protozoan Trichomonas. Transcripts associated with known mucin-degrading organisms and mucin-degrading enzymes were elevated in the fecal microbiomes of ICD-afflicted animals. Assessment of colon sections using immunohistochemistry and of the host transcriptome suggests differential fucosylation of mucins between control and ICD-afflicted animals. Interrogation of the metatranscriptome for fucose utilization genes reveals possible mechanisms by which opportunists persist in ICD. Bacteroides sp. potentially cross-fed fucose to Haemophilus whereas Campylobacter expressed a mucosa-associated transcriptome with increased expression of adherence genes. Conclusions The simultaneous profiling of bacterial, fungal, archaeal, protozoan, and macaque transcripts from stool samples reveals that ICD of rhesus macaques is associated with increased gene expression by pathogens, increased mucin degradation, and altered fucose utilization. The data suggest that the ICD-afflicted host produces fucosylated mucins that are leveraged by potentially pathogenic microbes as a carbon source or as adhesion sites.

The molecular complexity of the carbohydrates consumed by humans has been deceptively oversimplified due to a lack of analytical methods that possess the throughput, sensitivity, and resolution required to provide quantitative structural information. However, such information is becoming an integral part of understanding how specific glycan structures impact health through their interaction with the gut microbiome and host physiology. This work presents a detailed catalogue of the glycans present in complementary foods commonly consumed by toddlers during weaning and foods commonly consumed by American adults. The monosaccharide compositions of over 800 foods from diverse food groups including Fruits, Vegetables, Grain Products, Beans, Peas, Other Legumes, Nuts, Seeds; Sugars, Sweets and Beverages; Animal Products, and more were obtained and used to construct the “Davis Food Glycopedia” (DFG), an open-access database that provides quantitative structural information on the carbohydrates in food. While many foods within the same group possessed similar compositions, hierarchical clustering analysis revealed similarities between different groups as well. Such a Glycopedia can be used to formulate diets rich in specific monosaccharide residues to provide a more targeted modulation of the gut microbiome, thereby opening the door for a new class of prophylactic or therapeutic diets.

Among people of European descent, the ability to digest lactose into adulthood arose via strong positive selection of a highly advantageous allele encompassing the lactase gene. Lactose-tolerant and intolerant individuals may have different disease risks due to the shared genetics of their haplotype block. Therefore, the overall objective of the study was to assess the genetic association of the lactase persistence haplotype to disease risk. Using data from the 1000Genomes project, we estimated the size of the lactase persistence haplotype block to be 1.9 Mbp containing up to 9 protein-coding genes and a microRNA. Based on the function of the genes and microRNA, we studied health phenotypes likely to be impacted by the lactase persistence allele: prostate cancer status, cardiovascular disease status, and bone mineral density. We used summary statistics from large genome-wide metanalyses—32,965 bone mineral density, 140,306 prostate cancer and 184,305 coronary artery disease subjects—to evaluate whether the lactase persistence allele was associated with these disease phenotypes. Despite the fact that previous work demonstrated that the lactase persistence haplotype block harbors increased deleterious mutations, these results suggest little effect on the studied disease phenotypes.

Abstract Background Complex microbial communities are an area of growing interest in biology. Metatranscriptomics allows researchers to quantify microbial gene expression in an environmental sample via high-throughput sequencing. Metatranscriptomic experiments are computationally intensive because the experiments generate a large volume of sequence data and each sequence must be compared with reference sequences from thousands of organisms. Results SAMSA2 is an upgrade to the original Simple Annotation of Metatranscriptomes by Sequence Analysis (SAMSA) pipeline that has been redesigned for standalone use on a supercomputing cluster. SAMSA2 is faster due to the use of the DIAMOND aligner, and more flexible and reproducible because it uses local databases. SAMSA2 is available with detailed documentation, and example input and output files along with examples of master scripts for full pipeline execution. Conclusions SAMSA2 is a rapid and efficient metatranscriptome pipeline for analyzing large RNA-seq datasets in a supercomputing cluster environment. SAMSA2 provides simplified output that can be examined directly or used for further analyses, and its reference databases may be upgraded, altered or customized to fit the needs of any experiment.

Mechanistic research suggests a unique evolutionary relationship between complex milk oligosaccharides and cognate bifidobacteria enriched in breast-fed infants. Bovine milk oligosaccharides (BMO) were recently identified as structurally and functionally similar to human milk oligosaccharides. The present single-blind three-way crossover study is the first to determine the safety and tolerability of BMO consumption by healthy human participants (n 12) and its effects on faecal microbiota and microbial metabolism. Participants consumed each supplement (placebo-control; low- and high-BMO doses) for eleven consecutive days, followed by a 2-week washout period prior to initiating the next supplement arm. Low and high BMO doses were consumed as 25 and 35 % of each individual's daily fibre intake, respectively. Safety and tolerability were measured using standardised questionnaires on gut and stomach discomfort and stool consistency. Faecal extracts were profiled for bacterial populations by next-generation sequencing (NGS) and bifidobacteria presence was confirmed using quantitative PCR. Urine was analysed for changes in microbial metabolism using nuclear magnetic resonance spectroscopy (1H-NMR). Consumption of both the low and high BMO doses was well tolerated and did not change stool consistency from baseline. Multivariate analysis of the NGS results demonstrated no change in faecal microbiota phyla among the placebo-control and BMO supplement groups. In conclusion, BMO supplementation was well tolerated in healthy adults and has the potential to shift faecal microbiota toward beneficial strains as part of a synbiotic treatment with probiotic cultures that selectively metabolise oligosaccharides.

Lactase persistence (LP) is a trait in which lactose can be digested throughout adulthood, while lactase non-persistence (LNP) can cause lactose intolerance and influence dairy consumption. One single nucleotide polymorphism (SNP ID: rs4988235) is often used as a predictor for dairy intake, since it is responsible for LP in people in European descent, and can occur in other ethnic groups. The objective of this study was to determine whether rs4988235 genotypes and ethnicity influence reported dairy consumption in the United States (U.S.). A food frequency questionnaire (FFQ) and multiple Automated Self-Administered 24-h recalls (ASA24®) were used to measure habitual and recent intake, respectively, of total dairy, cheese, cow’s milk, plant-based alternative milk, and yogurt in a multi-ethnic U.S. cohort genotyped for rs4988235. Within Caucasian subjects, LP individuals reported consuming more recent total dairy and habitual total cow’s milk intake. For subjects of all ethnicities, LP individuals consumed more cheese (FFQ p = 0.043, ASA24 p = 0.012) and recent total dairy (ASA24 p = 0.005). For both dietary assessments, Caucasians consumed more cheese than all non-Caucasians (FFQ p = 0.036, ASA24 p = 0.002) independent of genotype, as well as more recent intake of yogurt (ASA24 p = 0.042). LP subjects consumed more total cow’s milk than LNP, but only when accounting for whether subjects were Caucasian or not (FFQ p = 0.015). Fluid milk and alternative plant-based milk consumption were not associated with genotypes or ethnicity. Our results show that both LP genotype and ethnicity influence the intake of some dairy products in a multi-ethnic U.S. cohort, but the ability of rs4988235 genotypes to predict intake may depend on ethnic background, the specific dairy product, and whether intake is reported on a habitual or recent basis. Therefore, ethnicity and the dietary assessment method should also be considered when determining the suitability of rs4988235 as a proxy for dairy intake.

Peroxisome proliferator-activated receptor (PPAR) agonists are currently used therapeutically in humans, even though many of their direct gene targets are unknown. Because PPARs can directly regulate gene expression through peroxisome proliferator response elements (PPREs), we pursued the computational prediction of PPREs on a genome-wide scale. Contrary to current hypotheses, PPREs are not isotype-specific, nor do flanking nucleotides confer additional information. However, a position weight matrix-based search for PPREs within upstream conserved elements yielded sufficient selectivity for a genome-wide search. Additionally, a novel motif occurring with greater prevalence than PPREs was revealed. Microarray and gene ontology analyses further validated our search technique and provided new functional clusters of genes that were not previously known to be directly regulated by PPARs (e.g., chromatin remodeling, DNA damage response, Wnt, and mitogen-activated protein kinase signaling). This first genome-wide library of high-confidence predicted PPAR target genes will be a valuable resource to PPAR biologists.

Over half of all children with autism spectrum disorders (ASD) have gastrointestinal (GI) co-morbidities including chronic constipation, diarrhea, and irritable bowel syndrome. The severity of these symptoms has been correlated with the degree of GI microbial dysbiosis. The study objective was to assess tolerability of a probiotic (Bifidobacterium infantis) in combination with a bovine colostrum product (BCP) as a source of prebiotic oligosaccharides and to evaluate GI, microbiome and immune factors in children with ASD and GI co-morbidities. This pilot study is a randomized, double blind, controlled trial of combination treatment (BCP + B. infantis) vs. BCP alone in a cross-over study in children ages 2-11 with ASD and GI co-morbidities (n = 8). This 12-week study included 5 weeks of probiotic-prebiotic supplementation, followed by a two-week washout period, and 5 weeks of prebiotic only supplementation. The primary outcome of tolerability was assessed using validated questionnaires of GI function and atypical behaviors, along with side effects. Results suggest that the combination treatment is well-tolerated in this cohort. The most common side effect was mild gassiness. Some participants on both treatments saw a reduction in the frequency of certain GI symptoms, as well as reduced occurrence of particular aberrant behaviors. Improvement may be explained by a reduction in IL-13 and TNF-α production in some participants. Although limited conclusions can be drawn from this small pilot study, the results support the need for further research into the efficacy of these treatments.

Collaboration among scientists has a major influence on scientific progress. Such collaboration often results from scientific meetings, where scientists gather to present and discuss their research and to meet potential collaborators. However, most scientific meetings have inherent biases, such as the availability of research funding or the selection bias of professional societies that make it difficult to study the effect of the meeting per se on scientific productivity. To evaluate the effects of scientific meetings on collaboration and progress independent of these biases, we conducted a study of the annual symposia held by the International Milk Genomics Consortium (IMGC) over a 12-year period. In our study, we conducted permutation testing to analyze the effectiveness of the IMGC in facilitating collaboration and productivity in a community of milk scientists who were meeting attendees relative to non-attendees. Using the number of co-authorships on published papers as a measure of collaboration, our analysis revealed that scientists who attended the symposium were associated with more collaboration than were scientists who did not attend. Furthermore, we evaluated the scientific progress of consortium attendees by analyzing publication rate and article impact. We found that IMGC attendees, in addition to being more collaborative, were also more productive and influential than were non-attendees who published in the same field. The results of our study suggest that the annual symposium encouraged interactions among disparate scientists and increased research productivity, exemplifying the positive effect of scientific meetings on both collaboration and progress.

Unlike other tissues, development and differentiation of the mammary gland occur mostly after birth. The roles of systemic hormones and local growth factors important for this development and functional differentiation are well-studied. In other tissues, it has been shown that chromatin organization plays a key role in transcriptional regulation and underlies epigenetic regulation during development and differentiation. However, the role of chromatin organization in mammary gland development and differentiation is less well-defined. Here, we have studied the changes in chromatin organization at the milk protein gene loci (casein, whey acidic protein, and others) in the mouse mammary gland before and after functional differentiation.Distal regulatory elements within the casein gene cluster and whey acidic protein gene region have an open chromatin organization after pubertal development, while proximal promoters only gain open-chromatin marks during pregnancy in conjunction with the major induction of their expression. In contrast, other milk protein genes, such as alpha-lactalbumin, already have an open chromatin organization in the mature virgin gland. Changes in chromatin organization in the casein gene cluster region that are present after puberty persisted after lactation has ceased, while the changes which occurred during pregnancy at the gene promoters were not maintained. In general, mammary gland expressed genes and their regulatory elements exhibit developmental stage- and tissue-specific chromatin organization.A progressive gain of epigenetic marks indicative of open/active chromatin on genes marking functional differentiation accompanies the development of the mammary gland. These results support a model in which a chromatin organization is established during pubertal development that is then poised to respond to the systemic hormonal signals of pregnancy and lactation to achieve the full functional capacity of the mammary gland.

Even in genomes lacking operons, a gene's position in the genome influences its potential for expression. The mechanisms by which adjacent genes are co-expressed are still not completely understood. Using lactation and the mammary gland as a model system, we explore the hypothesis that chromatin state contributes to the co-regulation of gene neighborhoods. The mammary gland represents a unique evolutionary model, due to its recent appearance, in the context of vertebrate genomes. An understanding of how the mammary gland is regulated to produce milk is also of biomedical and agricultural importance for human lactation and dairying. Here, we integrate epigenomic and transcriptomic data to develop a comprehensive regulatory model. Neighborhoods of mammary-expressed genes were determined using expression data derived from pregnant and lactating mice and a neighborhood scoring tool, G-NEST. Regions of open and closed chromatin were identified by ChIP-Seq of histone modifications H3K36me3, H3K4me2, and H3K27me3 in the mouse mammary gland and liver tissue during lactation. We found that neighborhoods of genes in regions of uniquely active chromatin in the lactating mammary gland, compared with liver tissue, were extremely rare. Rather, genes in most neighborhoods were suppressed during lactation as reflected in their expression levels and their location in regions of silenced chromatin. Chromatin silencing was largely shared between the liver and mammary gland during lactation, and what distinguished the mammary gland was mainly a small tissue-specific repertoire of isolated, expressed genes. These findings suggest that an advantage of the neighborhood organization is in the collective repression of groups of genes via a shared mechanism of chromatin repression. Genes essential to the mammary gland's uniqueness are isolated from neighbors, and likely have less tolerance for variation in expression, properties they share with genes responsible for an organism's survival.

Aware of the important benefits of human milk, most U.S. women initiate breastfeeding but difficulties with milk supply lead some to quit earlier than intended. Yet, the contribution of maternal physiology to lactation difficulties remains poorly understood. Human milk fat globules, by enveloping cell contents during their secretion into milk, are a rich source of mammary cell RNA. Here, we pair this non-invasive mRNA source with RNA-sequencing to probe the milk fat layer transcriptome during three stages of lactation: colostral, transitional, and mature milk production. The resulting transcriptomes paint an exquisite portrait of human lactation. The resulting transcriptional profiles cluster not by postpartum day, but by milk Na:K ratio, indicating that women sampled during similar postpartum time frames could be at markedly different stages of gene expression. Each stage of lactation is characterized by a dynamic range (10(5)-fold) in transcript abundances not previously observed with microarray technology. We discovered that transcripts for isoferritins and cathepsins are strikingly abundant during colostrum production, highlighting the potential importance of these proteins for neonatal health. Two transcripts, encoding β-casein (CSN2) and α-lactalbumin (LALBA), make up 45% of the total pool of mRNA in mature lactation. Genes significantly expressed across all stages of lactation are associated with making, modifying, transporting, and packaging milk proteins. Stage-specific transcripts are associated with immune defense during the colostral stage, up-regulation of the machinery needed for milk protein synthesis during the transitional stage, and the production of lipids during mature lactation. We observed strong modulation of key genes involved in lactose synthesis and insulin signaling. In particular, protein tyrosine phosphatase, receptor type, F (PTPRF) may serve as a biomarker linking insulin resistance with insufficient milk supply. This study provides the methodology and reference data set to enable future targeted research on the physiological contributors of sub-optimal lactation in humans.