Search

BackgroundThe accurate identification of the functional elements in the bovine genome is a fundamental requirement for high-quality analysis of data informing both genome biology and genomic selection. Functional annotation of the bovine genome was performed to identify a more complete catalog of transcript isoforms across bovine tissues.ResultsA total of 160,820 unique transcripts (50% protein coding) representing 34,882 unique genes (60% protein coding) were identified across tissues. Among them, 118,563 transcripts (73% of the total) were structurally validated by independent datasets (PacBio isoform sequencing data, Oxford Nanopore Technologies sequencing data, de novo assembled transcripts from RNA sequencing data) and comparison with Ensembl and NCBI gene sets. In addition, all transcripts were supported by extensive data from different technologies such as whole transcriptome termini site sequencing, RNA Annotation and Mapping of Promoters for the Analysis of Gene Expression, chromatin immunoprecipitation sequencing, and assay for transposase-accessible chromatin using sequencing. A large proportion of identified transcripts (69%) were unannotated, of which 86% were produced by annotated genes and 14% by unannotated genes. A median of two 5' untranslated regions were expressed per gene. Around 50% of protein-coding genes in each tissue were bifunctional and transcribed both coding and noncoding isoforms. Furthermore, we identified 3,744 genes that functioned as noncoding genes in fetal tissues but as protein-coding genes in adult tissues. Our new bovine genome annotation extended more than 11,000 annotated gene borders compared to Ensembl or NCBI annotations. The resulting bovine transcriptome was integrated with publicly available quantitative trait loci data to study tissue-tissue interconnection involved in different traits and construct the first bovine trait similarity network.ConclusionsThese validated results show significant improvement over current bovine genome annotations.

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.

Abstract Background In previous studies, gene neighborhoods—spatial clusters of co-expressed genes in the genome—have been defined using arbitrary rules such as requiring adjacency, a minimum number of genes, a fixed window size, or a minimum expression level. In the current study, we developed a Gene Neighborhood Scoring Tool (G-NEST) which combines genomic location, gene expression, and evolutionary sequence conservation data to score putative gene neighborhoods across all possible window sizes simultaneously. Results Using G-NEST on atlases of mouse and human tissue expression data, we found that large neighborhoods of ten or more genes are extremely rare in mammalian genomes. When they do occur, neighborhoods are typically composed of families of related genes. Both the highest scoring and the largest neighborhoods in mammalian genomes are formed by tandem gene duplication. Mammalian gene neighborhoods contain highly and variably expressed genes. Co-localized noisy gene pairs exhibit lower evolutionary conservation of their adjacent genome locations, suggesting that their shared transcriptional background may be disadvantageous. Genes that are essential to mammalian survival and reproduction are less likely to occur in neighborhoods, although neighborhoods are enriched with genes that function in mitosis. We also found that gene orientation and protein-protein interactions are partially responsible for maintenance of gene neighborhoods. Conclusions Our experiments using G-NEST confirm that tandem gene duplication is the primary driver of non-random gene order in mammalian genomes. Non-essentiality, co-functionality, gene orientation, and protein-protein interactions are additional forces that maintain gene neighborhoods, especially those formed by tandem duplicates. We expect G-NEST to be useful for other applications such as the identification of core regulatory modules, common transcriptional backgrounds, and chromatin domains. The software is available at http://docpollard.org/software.html

This paper resulted from a conference entitled "Lactation and Milk: Defining and refining the critical questions" held at the University of Colorado School of Medicine from January 18-20, 2012. The mission of the conference was to identify unresolved questions and set future goals for research into human milk composition, mammary development and lactation. We first outline the unanswered questions regarding the composition of human milk (Section I) and the mechanisms by which milk components affect neonatal development, growth and health and recommend models for future research. Emerging questions about how milk components affect cognitive development and behavioral phenotype of the offspring are presented in Section II. In Section III we outline the important unanswered questions about regulation of mammary gland development, the heritability of defects, the effects of maternal nutrition, disease, metabolic status, and therapeutic drugs upon the subsequent lactation. Questions surrounding breastfeeding practice are also highlighted. In Section IV we describe the specific nutritional challenges faced by three different populations, namely preterm infants, infants born to obese mothers who may or may not have gestational diabetes, and infants born to undernourished mothers. The recognition that multidisciplinary training is critical to advancing the field led us to formulate specific training recommendations in Section V. Our recommendations for research emphasis are summarized in Section VI. In sum, we present a roadmap for multidisciplinary research into all aspects of human lactation, milk and its role in infant nutrition for the next decade and beyond.

BackgroundThe newly assembled Bos taurus genome sequence enables the linkage of bovine milk and lactation data with other mammalian genomes.ResultsUsing publicly available milk proteome data and mammary expressed sequence tags, 197 milk protein genes and over 6,000 mammary genes were identified in the bovine genome. Intersection of these genes with 238 milk production quantitative trait loci curated from the literature decreased the search space for milk trait effectors by more than an order of magnitude. Genome location analysis revealed a tendency for milk protein genes to be clustered with other mammary genes. Using the genomes of a monotreme (platypus), a marsupial (opossum), and five placental mammals (bovine, human, dog, mice, rat), gene loss and duplication, phylogeny, sequence conservation, and evolution were examined. Compared with other genes in the bovine genome, milk and mammary genes are: more likely to be present in all mammals; more likely to be duplicated in therians; more highly conserved across Mammalia; and evolving more slowly along the bovine lineage. The most divergent proteins in milk were associated with nutritional and immunological components of milk, whereas highly conserved proteins were associated with secretory processes.ConclusionsAlthough both copy number and sequence variation contribute to the diversity of milk protein composition across species, our results suggest that this diversity is primarily due to other mechanisms. Our findings support the essentiality of milk to the survival of mammalian neonates and the establishment of milk secretory mechanisms more than 160 million years ago.

Mitochondria operate as a central hub for many metabolic processes by sensing and responding to the cellular environment. Developmental cues from the environment have been implicated in selective autophagy, or mitophagy, of mitochondria during cell differentiation and tissue development. Mitophagy occurring in this context, termed programmed mitophagy, responds to cell state rather than mitochondrial damage and is often accompanied by a metabolic transition. However, little is known about the mechanisms that engage and execute mitophagy under physiological or developmental conditions. As the mammary gland undergoes post-natal development and lactation challenges mitochondrial homeostasis, we investigated the contribution of mitochondria to differentiation of mammary epithelial cells (MECs). Using lactogenic differentiation of the HC11 mouse MEC line, we demonstrated that HC11 cells transition to a highly energetic state during differentiation by engaging both oxidative phosphorylation and glycolysis. Interestingly, this transition was lost when autophagy was inhibited with bafilomycin A1 or knockdown of Atg7 (autophagy related 7). To evaluate the specific targeting of mitochondria, we traced mitochondrial oxidation and turnover in vitro with the fluorescent probe, pMitoTimer. Indeed, we found that differentiation engaged mitophagy. To further evaluate the requirement of mitophagy during differentiation, we knocked down the expression of Prkn/parkin in HC11 cells. We found that MEC differentiation was impaired in shPrkn cells, implying that PRKN is required for MEC differentiation. These studies suggest a novel regulation of MEC differentiation through programmed mitophagy and provide a foundation for future studies of development and disease associated with mitochondrial function in the mammary gland. Abbreviations: AA: antimycin A; ATG5: autophagy related 5; BAF: bafilomycin A1; BNIP3: BCL2 interacting protein 3; BNIP3L/NIX: BCL2 interacting protein 3 like; COX8A: cytochrome c oxidase subunit 8A; CQ: chloroquine; CSN2: casein beta; ECAR: extracellular acidification rate; FCCP: trifluoromethoxy carbonylcyanide phenylhydrazone; FUNDC1: FUN14 domain containing 1; HIF1A: hypoxia inducible factor 1 subunit alpha; L1: lactation day 1; MAP1LC3B: microtubule associated protein 1 light chain 3 beta; MEC: mammary epithelial cell; mitoQ: mitoquinol; mROS: mitochondrial reactive oxygen species; OCR: oxygen consumption rate; P: priming; P16: pregnancy day 16; PARP1: poly(ADP-ribose) polymerase 1; PINK1: PTEN induced kinase 1; PPARGC1A: PPARG coactivator 1 alpha; PRKN: parkin RBR E3 ubiquitin protein ligase; shNT: short hairpin non-targeting control; SQSTM1: sequestosome 1; STAT3: signal transducer and activator of transcription 3; TEM: transmission electron microscopy; TFAM: transcription factor A, mitochondrial; U: undifferentiated. [ABSTRACT FROM AUTHOR]

The molecular clock plays key roles in daily physiological functions, development and cancer. Period 2 (PER2) is a repressive element, which inhibits transcription activated by positive clock elements, resulting in diurnal cycling of genes. However, there are gaps in our understanding of the role of the clock in normal development outside of its time-keeping function. Here, we show that PER2 has a noncircadian function that is crucial to mammalian mammary gland development. Virgin Per2-deficient mice, Per2-/-, have underdeveloped glands, containing fewer bifurcations and terminal ducts than glands of wildtype mice. Using a transplantation model, we show that these changes are intrinsic to the gland and further identify changes in cell fate commitment. Per2-/- mouse mammary glands have a dual luminal/basal phenotypic character in cells of the ductal epithelium. We identified colocalization of E-cadherin and keratin 14 in luminal cells. Similar results were demonstrated using MCF10A and shPER2 MCF10A human cell lines. Collectively this study reveals a crucial noncircadian function of PER2 in mammalian mammary gland development, validates the Per2-/- model, and describes a potential role for PER2 in breast cancer. [ABSTRACT FROM AUTHOR]