BACKGROUND
micrornas  are small noncoding rna molecules that modulate gene expression by targeting mrna transcripts for degradation or interfering with mrna translation  <cit> .
these 21- to 22-nucleotide molecules are derived from the processing of longer mrnas transcribed from intergenic regions or introns of protein-coding genes  <cit> .
approximately  <dig> mirna genes have been recognized in the human and mouse genomes, many of which are expressed in a tissue-dependent manner  <cit> .
bioinformatic predictions estimate that 30–80% of mammalian mrnas are targeted by mirnas, and a given mrna may be targeted by multiple mirnas  <cit> .
in addition, a single mirna can regulate multiple mrna transcripts, potentially orchestrating coordinate regulation of several genes within a metabolic pathway  <cit> .
it is generally thought that mirnas act to fine-tune mrna levels, but there is also evidence of full repression of specific protein production in some cases .
mirnas have important roles in the regulation of metabolic homeostasis.
for example, mir-33a and mir-33b, which are embedded within genes for sterol regulatory element-binding proteins srebf <dig> and srebf <dig>  respectively, have key roles in the modulation of cholesterol homeostasis  <cit> .
the discovery of the cellular roles of mir-33a/b in the repression of genes involved in cellular cholesterol export has revealed a new therapeutic target.
indeed, studies in mice and non-human primates have shown that antagonism of mir-33a/b reduces plasma lipid levels and atherosclerosis  <cit> .
many other examples of important roles for mirnas in metabolism have emerged recently .
among these are roles for mirnas in the regulation of adipogenesis .
specific mirnas that enhance adipocyte differentiation  or inhibit adipocyte differentiation  have been identified.
some mirnas that are induced during adipocyte development are dysregulated in obese mice  <cit> .
while a few specific mirnas have been well characterized, the roles of the majority of mirnas expressed in adipose tissue are unknown.
many properties of mammalian adipose tissue accumulation, distribution, and metabolism differ between males and females  <cit> .
to understand the basis for sex differences in fat storage, previous studies have assessed differences in mrna expression levels in male vs. females fat depots  <cit> .
it is likely that sex differences also occur in mirna expression levels in adipose tissue, and that these differences may contribute to sex differences in mrna levels.
consistent with this, sex differences in mirna expression levels have been reported in brain  <cit> , lung  <cit> , and liver  <cit>  and may influence sex-specific disease development and pathogenesis  <cit> .
however, the effect of sex on mirna expression in adipose tissue has not been investigated.
the mechanisms underlying sex differences in mirna levels have not been studied.
sex differences in metabolism can be attributed to both hormonal and genetic factors  <cit> .
gonadal hormones have been considered to be the primary drivers of sex differences, and some mirna levels in a variety of tissues are responsive to estradiols  <cit> .
importantly, however, the sex chromosome complement also plays a major role in determining adiposity  <cit> , and it is possible that this is mediated in part by mirnas.
to dissect the contributions of gonadal secretions and sex chromosome complement to sex differences, we have used the four core genotypes  mouse model.
the fcg model generates mice with four combinations of gonads and sex chromosomes: xx female, xx male, xy female, and xy male mice  <cit> .
fcg mice have a y chromosome that carries a mutation  in the testis-determining gene, sry, such that xy− mice develop ovaries rather than testes.
in addition, fcg mice carry an autosomal sry transgene that independently segregates from the y− chromosome.
mice that inherit the sry transgene develop testes.
a cross between xx females and xy−
sry males produces four combinations of gonads and sex chromosomes .
using this model, we previously demonstrated that the presence of two x chromosomes leads to increased adiposity compared to xy mice, independently from effects of ovaries or testes  <cit> .
here, we performed rna sequencing  of small rnas present in gonadal fat from fcg mice to identify sex differences in mirna expression levels in adipose tissue.
we determined that sex hormones and sex chromosomes each influence the mirna expression profile.
in addition, comparison of mice fed chow vs. high fat diets revealed sex-specific changes in mirna profiles in response to diet-induced obesity.
our findings have implications for understanding sex differences in adipose tissue metabolism and the development of diet-induced obesity.
methods
animals
four core genotypes  c57bl/ <dig> mice were bred and genotyped as described previously  <cit> .
briefly, xx female mice were mated with xy− male mice to generate xx, xx, xy−, and xy− offspring, and genotyping was performed by pcr to detect the sry transgene and a y-chromosome–specific sequence.
where indicated, gonadectomy was performed at 75 days of age as described previously  <cit> .
during surgery, gonads were removed while leaving surrounding adipose tissue in place.
gonadal males and females were housed in separate cages and maintained at 23 °c with a 12:12 h light:dark cycle.
gonadally intact females were analyzed without estrous cycle synchronization, such that gene expression values represent an average over the estrous phases.
all mice were initially fed purina mouse chow diet containing 5% fat .
where specified, mice were fed a high fat diet  for 16 weeks beginning at  <dig>  months of age .
adipose tissue was harvested from all mice at  <dig>  months of age.
mouse studies were conducted after approval by the institutional animal research committee of the university of california, los angeles.
rna extraction and quality control
at the time of dissection, gonadal fat tissue was flash frozen in liquid nitrogen and stored at −80 °c.
small rnas were isolated from 100 mg tissue samples using qiazol and qiagen’s mirneasy mini kit .
after homogenization, samples were centrifuged at  <dig>  × g for 10 min to separate the transparent lipid layer from the pink organic layer.
only the organic layer was used in chloroform extraction.
all subsequent steps followed the qiagen protocol.
rna samples were submitted to agilent bioanalyzer eukaryote total nano-rna chip analysis, yielding rna integrity numbers of  <dig>  or greater.
mirna library preparation
the feasibility of using pooled sequencing libraries was assessed by sequencing mirnas from individual samples separately and after pooling.
indexed libraries were generated from adipose tissue of three xx females fed a high fat diet.
average counts of mapped mirnas from the individual libraries were compared to mirna counts of the pooled library using pearson’s product-moment correlation.
for the remaining conditions, three samples of each genotype were pooled into equimolar amounts for library preparation.
in total, twelve mirna libraries were made: libraries for each of the four genotypes in chow-fed, gonadally intact mice, chow-fed, gonadectomized mice, and gonadectomized mice fed a high fat diet.
mirna libraries were processed individually using a standard protocol from illumina truseq small rna kit, with indices 1– <dig>  and gel purified according to manufacturer’s instructions.
final sequencing library concentration  was determined using kapa library quantification qpcr kit .
sequencing was performed at the broad stem cell research center core facility at ucla, on illumina hiseq  <dig>
reference sequence determination
mirna gene expression is typically quantified by counting reads that map to the mirna genes.
however, in some families of mirnas, several genes give rise to identical mature sequences, thus it is impossible to distinguish which mirna gene gave rise to the mature sequence based on sequencing alone.
we performed reference preprocessing to compile a list of uniquely expressed mouse mirna sequences, regardless of their gene of origin, so that quantification was done at the level of mature mirna sequence, rather than at the gene level.
a reference sequence was compiled based on all mature and precursor mirna sequences available from the main repository of mirna studies, mirbase version  <dig> , and included all mature and precursor sequences.
the respective mature sequences within precursor sequences were masked to “n”, to prevent mature sequences mapping both to mature and precursor.
thus, reads mapping to precursor sequences mapped to precursors exclusively.
in cases where several genes gave rise to the same mature mirna sequence, that sequence was represented only once in the reference under one name, with an additional column listing all matching genes as potential targets in the results.
in total,  <dig> mature mirna sequences were unique , while an additional  <dig> mirnas grouped into  <dig> unique sequences .
in  <dig> cases, the mirnas were not completely equivalent, with one mirna being  <dig> base shorter than the other.
in such cases, the shorter sequence would map to both mature sequences, while the longer one would map uniquely to the longer sequence.
for these pairs, the counts of the shorter sequence were determined as the total counts mapping to longer sequence minus the unique mappings to longer sequence .
sequence read processing
reads were demultiplexed in.qseq format based on perfect match to barcode sequence, and converted to.fastq format using in-house perl scripts .
read quality was assessed with fastqc  and adaptor sequence was clipped using cutadapt tool .
clean reads 18–30 bp long were retained for subsequent mapping and analysis.
identical reads were collapsed using collapse.pl script from mirdeep package .
reads were then mapped to the reference with bwa aligner allowing for up to  <dig> mismatch and no gaps.
using these criteria, 80–90% of the reads mapped to the reference sequence.
read counting and result table
reads that uniquely mapped to one target were counted towards that target.
reads that had both perfect and imperfect  matches, were counted only towards the perfect match target.
results were summarized in a count table that listed the counts for each target in each sample .
precursor mirna counts are composed of reads mapping exclusively to precursor sequences, and not to any mature mirna sequence in the reference.
for each reference sequence, the total read count and the number of reads that originate from 1-mismatch reads were presented separately.
the number of mismatched reads is a subset of total reads mapping to the reference.
percent abundance of mirnas was calculated by normalizing the reads from one mirna to the total number of reads mapped for all mirnas.
the mean of the four genotypes across all diet/gonadectomy conditions is represented in fig.
2b.
quantitative pcr
select mirnas were validated with quantitative pcr, using individual samples that comprised the pooled samples that were sequenced.
taqman primers  were used according to manufacturer’s protocols.
briefly, 10 ng rna was reverse transcribed using taqman microrna reverse transcription kit .
quantitative real-time pcr was performed on bio-rad cfx connect real-time pcr detection system using ssofast probes supermix .
mirnas of interest were normalized to two housekeeping genes: u <dig> small nuclear rna and small nucleolar rna   <dig>  expression was quantified using Δct for mir- <dig> or standard curve for mir-133a, mir- <dig>  and mir- <dig>  standard curves were generated using  <dig> serial dilution points of cdna combined from all biological samples.
mirna target prediction
mirnas with exclusively male–female or xx–xy log <dig> fold changes greater than  <dig>  were selected for target prediction.
targets for mirnas were predicted by mirdb.org, version  <dig>   <cit> .
all targets with a prediction score of  <dig> or above were included in subsequent functional analysis using david, version  <dig>   <cit> .
the functional annotation tool was used to identify enriched annotation terms in the kegg pathway.
unique kegg pathways with a bonferroni-corrected p-value <  <dig>  were included in additional file 5: table s <dig>
statistics
pearson’s product-moment correlation, exact binomial tests, wilcoxon rank sum tests, and principal component analyses were performed using r and visualized with the r packages ggplot <dig> and ggbiplot .
because each genotype and diet condition was measured in one pooled sample from each genotype/diet/sex condition, the present analysis did not have adequate power to assess the significance of treatment and genotypic effects on expression levels of individual mirnas.
instead, we viewed each mirna as an independent contributor to the overall mirna expression profile of each genotype under each experimental condition, and we report patterns of expression among groups.
in qpcr experiments, groups were compared using two-way anova , with gonadal sex and sex chromosome complement as independent factors.
values for mir- <dig> and mir- <dig> were log10-transformed before statistical analysis.
RESULTS
study design to identify sex differences in adipose tissue mirnas under three physiological conditions
to assess the effect of gonadal hormones and sex chromosomes on mirna levels in adipose tissue, we sequenced mirnas in the gonadal fat depot of fcg mice.
the four genotypes of this model allow analysis of sex differences based on the presence of male vs. female gonads , or the presence of xx vs. xy chromosome complement  .
additionally, we analyzed the effect of circulating gonadal hormones by comparing changes in mirna profiles in mice with intact gonads vs. mice after gonadectomy .
finally, we assessed the effect on adipose tissue mirna profile of high fat diet, and its interaction with sex, by comparing mice on chow and high fat diets .
adipose tissue samples from the gonadal fat depot were collected from all cohorts at the same age .
sequenced and mapped mirnas were counted to establish mirna expression profiles.fig.
<dig> study design to identify sex and dietary effects on mirna expression.
a mirnas were sequenced in adipose tissue from the four core genotypes mouse model .
a comparison of mirna levels in female xx and xy adipose tissue with those in male xx and xy tissue allows detection of effects due to gonadal type.
a comparison of mirna levels in female and male xx adipose tissue with those on female and male xy tissue allows detection of effects due to sex chromosome type.
b to determine the effects of acute gonadal hormones of diet, mirna sequencing was performed in four core genotypes cohorts that were fed a chow diet in the gonadally intact state, fed a chow diet and gonadectomized as adults to remove the acute effects of gonadal secretions, and in mice that were gonadectomized as adults and fed a high fat diet for 16 weeks
sex differences in adipose tissue mirna profiles
we sequenced and mapped  <dig>  mature and precursor mirnas in gonadal adipose tissue of fcg mice.
we identified  <dig> mirnas  that were expressed at substantial levels in gonadal fat  .
three mirna species—mir-10b-5p, mir-143-3p, and mir-22-3p—accounted for nearly 50% of all mirna reads .
an additional  <dig> mirna species each accounted for at least  <dig> % of the gonadal adipose tissue mirna pool.
<dig> of the  <dig>  most prevalent mirnas in mouse adipose tissue have previously been shown to be abundant in human subcutaneous white adipose tissue  <cit> .fig.
<dig> highly expressed mirnas in gonadal fat.
a out of  <dig> mapped mirnas,  <dig> were expressed at over  <dig> counts per million reads in at least one of the twelve libraries.
these mirnas are termed “highly expressed” mirnas.
b percent abundance of mirnas in gonadal fat of fcg mice.
values represent mean ± sd of the twelve sequencing libraries.
c, d three individual mirna sequencing libraries representing biological replicates of xx females fed a high fat diet are correlated with the pooled sequencing library composed of those same mice.
highly expressed gonadal fat mirnas  are more correlated between the individual and pooled libraries than mirnas that are expressed at a lower level .
values represent mean ± sd.
r <dig> and p-values were calculated using pearson product–moment correlation
to streamline our sequencing process for the twelve conditions shown in fig.
<dig>  we assessed the feasibility of pooling three biological replicates for mirna sequencing for each condition.
we sequenced samples from three individual mice  separately and after pooling together.
we found that the average mirna reads of the three individual samples were highly correlated with the values from the pooled sample .
based on these findings, we concluded that pooled libraries were a satisfactory representation of individual biological samples.
for analyses described below, we focused on the  <dig> most abundantly expressed adipose tissue mirnas, as these are most likely to have physiological significance.
to establish a baseline profile of sex effects on mirna expression in adipose tissue, we sequenced mirnas from mice fed a standard chow diet.
we first asked whether there was a male or female bias in overall mirna expression by examining the ratio of female-to-male mirna levels .
female-to-male ratios of each mirna were log2-transformed and plotted as a histogram .
adipose tissue mirnas had a significant bias toward higher expression in male compared to female gonadal adipose tissue .
in the chow-fed cohort, there was no significant bias in mirna distribution between xx and xy mice, although we cannot exclude the possibility that the sex chromosome complement may alter expression levels of individual mirna species.
thus, in the basal physiological condition of chow-fed animals, the presence of male vs. female gonads led to a significant difference in adipose tissue mirna expression profile.fig.
<dig> distribution of mirnas in chow-fed gonadally intact mice.
distribution of female-to-male ratios  and xx-to-xy ratios  in highly expressed mirnas.
ratios were calculated for each mirna, log2-transformed, and binned in  <dig>  increments.
values less than − <dig> were grouped into the first bin, and values greater than  <dig> were grouped into the last bin.
red bars refer to male -biased mirnas, black bars refer to female -biased mirnas, blue bars refer to xy-biased mirnas, and white bars refer to xx-biased mirnas.
p-values were calculated using the exact binomial test
adipose tissue mirnas regulated by gonadal hormones
we hypothesized that male–female differences in mirna expression profile were caused by differences in gonadal hormones.
gonadal hormones cause sex differences in two primary ways.
first, they impose permanent, or “organizational,” effects during development, such as the effect of testosterone on male genital differentiation.
second, male and female gonadal hormones in the circulation have acute/reversible effects on processes such as gene expression.
to investigate the effects of acute gonadal hormone action on mirna expression, we gonadectomized adult fcg mice  and assessed mirna profiles 5 months later, during which time circulating gonadal hormones were absent.
with this experimental design, male–female differences that disappear after gonadectomy can be attributed to acute hormone effects.
in gonadectomized mice, the mirna profile expressed as female-to-male ratio exhibited a small bias toward higher mirna expression levels in females compared to males .
comparison of female–male ratio distributions between gonadectomized mice and gonadally intact mice revealed a significant difference .
this represented a shift from a profile with male > female expression ratios in gonadally intact mice to a profile with a slight bias toward female > male ratios in gonadectomized mice .
the mirna ratio distribution of xx vs. xy mice was similar in intact and gonadectomized mice .
taken together, these data suggest that acute effects of gonadal secretions, but not sex chromosome complement, produce sex biases in overall mirna levels in mice fed a chow diet.
we note, however, that there are subsets of mirnas with extreme values in female-to-male ratios and in xx-to-xy ratios in gonadectomized mice, suggesting that, in the absence of circulating gonadal hormones, sex and sex chromosome differences remain in individual mirnas.fig.
<dig> gonadectomy alters distribution of mirnas in chow-fed mice.
distribution of female-to-male ratios  and xx-to-xy ratios  in highly expressed mirnas.
ratios were calculated for each mirna and log2-transformed in gonadally intact  and gonadectomized  fcg mice fed a chow diet.
data for gonadally intact  mice are identical to fig.
3a–b, but binned in  <dig>  increments and compared to data distribution for gdx mice.
each dot represents one of the  <dig> highly expressed mirnas.
black bars represent median values.
distributions were significantly different for female-to-male ratios, but not for xx-to-xy ratios.
p-values were calculated using the wilcoxon rank sum test.
f, female; m, male
adipose tissue mirna sex differences influenced by high fat diet or obesity
several studies have shown that mirna expression levels are regulated during adipogenesis or are altered in obesity  <cit> .
however, it is unknown whether sex influences mirna expression levels in obese adipose tissue.
to identify sex and/or sex chromosome factors that alter mirna expression in obesity apart from the acute effect of gonadal hormones, we performed mirna-seq on gonadal fat from fcg mice that were gonadectomized  and then fed a high fat diet  beginning 1 month later for 4 months.
the distribution of female-to-male mirna read count ratios was significantly male-biased .
this male bias in gonadectomized high fat diet-fed mice  represents a significant shift in mirna expression compared to the distribution in chow-fed mice that had been gonadectomized at the same age and therefore represents an effect of diet .
interestingly, mirna expression levels were also biased toward the xy sex chromosome complement in mice fed a high fat diet .
this indicates that the sex chromosome complement, in addition to gonadal hormones, influenced mirna expression in mice fed a high fat diet.
comparison of the distributions of xx-to-xy ratios indicated a shift from a slight xx bias in chow-fed mice  to an xy bias in fat-fed mice .
in addition, there are a number of mirnas that shift to the extreme ends of the sex chromosome distribution in mice fed a high fat diet.
these data suggest that, in gonadectomized mice, high fat diet and obesity modulate mirna expression in a sex- and sex chromosome-specific manner.fig.
<dig> high fat diet alters mirna expression levels in a sex-dependent and sex chromosome complement-dependent manner.
distribution of female-to-male ratios  and xx-to-xy ratios  in highly expressed mirnas.
ratios were calculated for each mirna and log2-transformed in gonadectomized  fcg mice fed a high fat diet .
data for chow-fed gonadectomized mice are identical to fig.
4a–b.
distributions were significantly different for female-to-male ratios and for xx-to-xy ratios.
p-values were calculated using the wilcoxon rank sum test.
f, female; m, male
we followed up on select mirnas that had distinct sex-influenced expression patterns and/or were previously implicated in adipocyte biology.
quantitative real-time pcr  of mir-133a-3p, previously shown to inhibit brown adipocyte differentiation  <cit> , was upregulated in xx compared to xy mice in the chow-fed, gonadally intact group .
this difference was abolished after gonadectomy, suggesting an interplay between gonadal hormones and sex chromosome complement.
mir-192-5p levels were higher in males compared to females in gonad-intact mice, but this sex difference is lost after gonadectomy, and an xx > xy difference emerges .
conversely, mir-205-5p showed a sex-chromosome bias  in gonadally intact mice and a sex bias  in gonadectomized mice fed a chow diet .
however, when mice were fed a high fat diet, mir-192-5p and mir-205-5p showed no sex or sex chromosome differences.
mir-221-3p, which is located on the x chromosome and has been previously associated with obesity and adipogenesis  <cit> , showed no sex or sex chromosome difference in mice fed a chow diet .
however, mice fed high fat diet showed increased levels of mir-221-3p in xx male and female mice compared to xy mice.
these mirna species illustrate several distinct patterns of expression that are dependent on gonadal type, gonadal hormones, diet, and sex chromosome complement, and emphasize that sex differences exist in mirnas that have been studied previously in the context of adipose tissue and obesity.
to examine the relationship between diet, acute hormone effects, gonadal type, and sex chromosome complement effects on mirna expression, we performed principal component analysis.
diet distinguished mirna expression levels to the greatest degree, explaining  <dig> % of mirna covariance .
gonadal state  and gonadal type  led to a milder but detectable separation .
the sex chromosome complement did not lead to divergence of mirna covariance.
because diet had such a pronounced effect on mirna covariance, we performed principal component analysis on mice fed only the chow diet to determine if sex-biasing factors still influenced mirna covariance.
the first two components revealed deviation of gonadally intact mice from gonadectomized mice .
in addition, when examining the third and fourth principal components, females clustered separately from males, and xx mice exhibited a slight degree of separation from xy mice .
we note that the effect of the sex chromosome complement was detectable after removing the overwhelming effect of diet.
taken together, these data suggest that diet is a major factor in mirna expression levels, and that sex factors, such as acute gonadal secretions, male–female gonads, and xx–xy chromosome complement, also play important roles in modulating mirna levels.fig.
<dig> principal component analysis  reveals factors influencing mirna correlation covariance.
pca of all twelve sequencing libraries  or chow-fed samples only .
each dot represents a sequencing library and consists of normalized read counts for each of the  <dig> highly expressed mirnas.
dots are colored and encircled according to diet , gonadal state , gonadal sex , or sex chromosome complement .
hfd, high fat diet; int, gonadally intact; gdx, gonadectomized; f, female; m, male
discussion
mirnas are key regulators of gene expression through their effects on the degradation or translation of mrna transcripts.
in metabolic tissues, mirnas influence gene expression related to adipogenesis, cholesterol homeostasis, and glucose homeostasis  <cit> .
these metabolic processes are also affected by sex.
while there is evidence for sexual dimorphism in mirna expression in brain, skeletal muscle, and immune cells  <cit> , sex effects on mirnas in adipose tissue have not been characterized.
in this study, we used the fcg mouse model to identify sex chromosome complement and gonadal hormone effects on mirna abundance in gonadal fat.
we sequenced mirnas in three distinct cohorts of fcg mice in different hormonal and dietary conditions to determine the effects of circulating gonadal hormones and high fat diet on sex differences in gonadal fat mirnas.
we selected the gonadal fat depot for our studies because it has been widely studied in the mouse obesity field.
however, it is clear that differences exist between gonadal fat and other visceral fat depots, as well as subcutaneous fat depots, and it would ultimately be of interest to delineate sex differences across these distinct depots.
our analyses focused on the most abundant mirnas, as we suspect that sex differences in their levels are the most likely to have significant physiological effects.
by comparing overall mirna profiles across diet, gonadal state, gonadal type, and sex chromosome complement, we determined that each of these factors influence mirna expression in white adipose tissue.
in gonadally intact mice, we identified a male bias in mirna expression levels and no significant bias between xx and xy mice.
as expected, the distribution of male–female differences in the levels of mirnas was altered by gonadectomy, suggesting that circulating gonadal hormones regulate mirna levels and are responsible for some male–female sex differences.
gonadectomy also influences adiposity  <cit> , which may impact mirna expression.
we have previously shown that gonadectomized xx mice have increased body weight and adiposity compared to gonadectomized xy mice  <cit> .
because we studied mice 5 months after gonadectomy, we cannot distinguish whether changes in mirna expression preceded increases in adiposity or resulted from them.
however, male–female differences in mirna levels do not correlate with xx–xy differences in body weight, suggesting that sex differences in adiposity do not determine sex differences in overall mirna levels.
studies measuring mirna levels after gonadectomy but before weight gain will be required to address these questions.
high fat diet feeding alters sex and sex chromosome differences in mirna levels, suggesting a diet–sex interaction in mirna expression.
it is notable that the observed sex and sex chromosome biases result from changes in overall mirna distributions and not from changes in a few specific mirnas.
high fat feeding revealed a sex chromosome bias in mirna expression ratios that was not apparent in chow-fed animals.
this raises the possibility that differential mirna expression levels may contribute to the dramatic difference in adiposity of xx compared to xy mice in response to high fat diet  <cit> .
thus, mice with two x chromosomes gain more fat than xy mice, and this difference is enhanced upon feeding a high fat diet  <cit> .
further studies will be required to elucidate the physiological roles of individual sex-biased mirnas in adipogenesis and diet-induced obesity.
some mirnas with distinct sex-biased patterns in expression, such as mir-133a and mir- <dig>  have been implicated in adipogenesis  <cit> .
other mirnas may have yet unidentified roles in lipid metabolism.
an initial analysis of potential target genes for sex-biased mirnas that we identified revealed enrichment in specific metabolic pathways .
these predictions suggest that pathways such as mtor signaling and phosphatidylinositol signaling may be important in mediating sex differences in healthy and/or diseased adipocytes.
studies to determine the impact of these predicted target genes and corresponding pathways on sex differences in obesity provide a unique area for future research.
CONCLUSIONS
this study represents a first look at mechanisms underlying sex differences in mirna profiles in adipose tissue.
we demonstrate that mirnas in adipose tissue are influenced by diet, gonadal state, gonadal type, and sex chromosome complement.
principal component analysis revealed that diet is the leading factor in mirna expression covariance, while gonadal state, gonadal type, and sex chromosome complement were milder but important regulators.
by recognizing the innate sex differences in mirna levels, we can better understand sex differences in fat accumulation, distribution, and adipose pathophysiology.
additional files
additional file 1:  <dig> mature mirna sequences that mapped uniquely to the mirbase  <dig>
additional file 2:  <dig> mirnas that mapped to  <dig> unique sequences.
additional file 3:  <dig> cases of mirnas that were of unequal length, and thus read counts had to be adjusted for non-unique mapping.
additional file 4: count table summarizing read counts for  <dig> mature and precursor mirna sequences in fifteen sequencing libraries.
perfectly matched counts are normalized to total library reads per million.
additional file 5: figure s <dig>  expression of select mirnas with distinct sex- and sex-chromosome-biased patterns, measured by taqman qpcr.
table s <dig>  enrichment of specific cellular pathways in predicted targets of sex-influenced mirnas.
abbreviations
fcgfour core genotypes
mir, mirnamicrorna
mrnamessenger rna

 micrornas are small noncoding rna molecules that modulate gene expression by targeting mrna transcripts for degradation or interfering with mrna translation. these 21- to 22-nucleotide molecules are derived from the processing of longer mrnas transcribed from intergenic regions or introns of protein-coding genes. approximately mirna genes have been recognized in the human and mouse genomes, many of which are expressed in a tissue-dependent manner. bioinformatic predictions estimate that 30–80% of mammalian mrnas are targeted by mirnas, and a given mrna may be targeted by multiple mirnas. in addition, a single mirna can regulate multiple mrna transcripts, potentially orchestrating coordinate regulation of several genes within a metabolic pathway. it is generally thought that mirnas act to fine-tune mrna levels, but there is also evidence of full repression of specific protein production in some cases. mirnas have important roles in the regulation of metabolic homeostasis. for example, mir-33a and mir-33b, which are embedded within genes for sterol regulatory element-binding proteins srebf and srebf respectively, have key roles in the modulation of cholesterol homeostasis. the discovery of the cellular roles of mir-33a/b in the repression of genes involved in cellular cholesterol export has revealed a new therapeutic target. indeed, studies in mice and non-human primates have shown that antagonism of mir-33a/b reduces plasma lipid levels and atherosclerosis. many other examples of important roles for mirnas in metabolism have emerged recently. among these are roles for mirnas in the regulation of adipogenesis. specific mirnas that enhance adipocyte differentiation or inhibit adipocyte differentiation have been identified. some mirnas that are induced during adipocyte development are dysregulated in obese mice. while a few specific mirnas have been well characterized, the roles of the majority of mirnas expressed in adipose tissue are unknown. many properties of mammalian adipose tissue accumulation, distribution, and metabolism differ between males and females. to understand the basis for sex differences in fat storage, previous studies have assessed differences in mrna expression levels in male vs. females fat depots. it is likely that sex differences also occur in mirna expression levels in adipose tissue, and that these differences may contribute to sex differences in mrna levels. consistent with this, sex differences in mirna expression levels have been reported in brain, lung, and liver and may influence sex-specific disease development and pathogenesis. however, the effect of sex on mirna expression in adipose tissue has not been investigated. the mechanisms underlying sex differences in mirna levels have not been studied. sex differences in metabolism can be attributed to both hormonal and genetic factors. gonadal hormones have been considered to be the primary drivers of sex differences, and some mirna levels in a variety of tissues are responsive to estradiols. importantly, however, the sex chromosome complement also plays a major role in determining adiposity, and it is possible that this is mediated in part by mirnas. to dissect the contributions of gonadal secretions and sex chromosome complement to sex differences, we have used the four core genotypes mouse model. the fcg model generates mice with four combinations of gonads and sex chromosomes: xx female, xx male, xy female, and xy male mice. fcg mice have a y chromosome that carries a mutation in the testis-determining gene, sry, such that xy− mice develop ovaries rather than testes. in addition, fcg mice carry an autosomal sry transgene that independently segregates from the y− chromosome. mice that inherit the sry transgene develop testes. a cross between xx females and xy− sry males produces four combinations of gonads and sex chromosomes. using this model, we previously demonstrated that the presence of two x chromosomes leads to increased adiposity compared to xy mice, independently from effects of ovaries or testes. here, we performed rna sequencing of small rnas present in gonadal fat from fcg mice to identify sex differences in mirna expression levels in adipose tissue. we determined that sex hormones and sex chromosomes each influence the mirna expression profile. in addition, comparison of mice fed chow vs. high fat diets revealed sex-specific changes in mirna profiles in response to diet-induced obesity. our findings have implications for understanding sex differences in adipose tissue metabolism and the development of diet-induced obesity. methods animals four core genotypes c57bl/ mice were bred and genotyped as described previously. briefly, xx female mice were mated with xy− male mice to generate xx, xx, xy−, and xy− offspring, and genotyping was performed by pcr to detect the sry transgene and a y-chromosome–specific sequence. where indicated, gonadectomy was performed at 75 days of age as described previously. during surgery, gonads were removed while leaving surrounding adipose tissue in place. gonadal males and females were housed in separate cages and maintained at 23 °c with a 12:12 h light:dark cycle. gonadally intact females were analyzed without estrous cycle synchronization, such that gene expression values represent an average over the estrous phases. all mice were initially fed purina mouse chow diet containing 5% fat. where specified, mice were fed a high fat diet for 16 weeks beginning at  months of age. adipose tissue was harvested from all mice at  months of age. mouse studies were conducted after approval by the institutional animal research committee of the university of california, los angeles. rna extraction and quality control at the time of dissection, gonadal fat tissue was flash frozen in liquid nitrogen and stored at −80 °c. small rnas were isolated from 100 mg tissue samples using qiazol and qiagen’s mirneasy mini kit. after homogenization, samples were centrifuged at  × g for 10 min to separate the transparent lipid layer from the pink organic layer. only the organic layer was used in chloroform extraction. all subsequent steps followed the qiagen protocol. rna samples were submitted to agilent bioanalyzer eukaryote total nano-rna chip analysis, yielding rna integrity numbers of or greater. mirna library preparation the feasibility of using pooled sequencing libraries was assessed by sequencing mirnas from individual samples separately and after pooling. indexed libraries were generated from adipose tissue of three xx females fed a high fat diet. average counts of mapped mirnas from the individual libraries were compared to mirna counts of the pooled library using pearson’s product-moment correlation. for the remaining conditions, three samples of each genotype were pooled into equimolar amounts for library preparation. in total, twelve mirna libraries were made: libraries for each of the four genotypes in chow-fed, gonadally intact mice, chow-fed, gonadectomized mice, and gonadectomized mice fed a high fat diet. mirna libraries were processed individually using a standard protocol from illumina truseq small rna kit, with indices 1– and gel purified according to manufacturer’s instructions. final sequencing library concentration was determined using kapa library quantification qpcr kit. sequencing was performed at the broad stem cell research center core facility at ucla, on illumina hiseq reference sequence determination mirna gene expression is typically quantified by counting reads that map to the mirna genes. however, in some families of mirnas, several genes give rise to identical mature sequences, thus it is impossible to distinguish which mirna gene gave rise to the mature sequence based on sequencing alone. we performed reference preprocessing to compile a list of uniquely expressed mouse mirna sequences, regardless of their gene of origin, so that quantification was done at the level of mature mirna sequence, rather than at the gene level. a reference sequence was compiled based on all mature and precursor mirna sequences available from the main repository of mirna studies, mirbase version, and included all mature and precursor sequences. the respective mature sequences within precursor sequences were masked to “n”, to prevent mature sequences mapping both to mature and precursor. thus, reads mapping to precursor sequences mapped to precursors exclusively. in cases where several genes gave rise to the same mature mirna sequence, that sequence was represented only once in the reference under one name, with an additional column listing all matching genes as potential targets in the results. in total, mature mirna sequences were unique, while an additional mirnas grouped into unique sequences. in cases, the mirnas were not completely equivalent, with one mirna being base shorter than the other. in such cases, the shorter sequence would map to both mature sequences, while the longer one would map uniquely to the longer sequence. for these pairs, the counts of the shorter sequence were determined as the total counts mapping to longer sequence minus the unique mappings to longer sequence. sequence read processing reads were demultiplexed in.qseq format based on perfect match to barcode sequence, and converted to.fastq format using in-house perl scripts. read quality was assessed with fastqc and adaptor sequence was clipped using cutadapt tool. clean reads 18–30 bp long were retained for subsequent mapping and analysis. identical reads were collapsed using collapse.pl script from mirdeep package. reads were then mapped to the reference with bwa aligner allowing for up to mismatch and no gaps. using these criteria, 80–90% of the reads mapped to the reference sequence. read counting and result table reads that uniquely mapped to one target were counted towards that target. reads that had both perfect and imperfect matches, were counted only towards the perfect match target. results were summarized in a count table that listed the counts for each target in each sample. precursor mirna counts are composed of reads mapping exclusively to precursor sequences, and not to any mature mirna sequence in the reference. for each reference sequence, the total read count and the number of reads that originate from 1-mismatch reads were presented separately. the number of mismatched reads is a subset of total reads mapping to the reference. percent abundance of mirnas was calculated by normalizing the reads from one mirna to the total number of reads mapped for all mirnas. the mean of the four genotypes across all diet/gonadectomy conditions is represented in fig. 2b. quantitative pcr select mirnas were validated with quantitative pcr, using individual samples that comprised the pooled samples that were sequenced. taqman primers were used according to manufacturer’s protocols. briefly, 10 ng rna was reverse transcribed using taqman microrna reverse transcription kit. quantitative real-time pcr was performed on bio-rad cfx connect real-time pcr detection system using ssofast probes supermix. mirnas of interest were normalized to two housekeeping genes: u small nuclear rna and small nucleolar rna expression was quantified using Δct for mir- or standard curve for mir-133a, mir- and mir- standard curves were generated using serial dilution points of cdna combined from all biological samples. mirna target prediction mirnas with exclusively male–female or xx–xy log fold changes greater than were selected for target prediction. targets for mirnas were predicted by mirdb.org, version. all targets with a prediction score of or above were included in subsequent functional analysis using david, version. the functional annotation tool was used to identify enriched annotation terms in the kegg pathway. unique kegg pathways with a bonferroni-corrected p-value <  were included in additional file 5: table s statistics pearson’s product-moment correlation, exact binomial tests, wilcoxon rank sum tests, and principal component analyses were performed using r and visualized with the r packages ggplot and ggbiplot. because each genotype and diet condition was measured in one pooled sample from each genotype/diet/sex condition, the present analysis did not have adequate power to assess the significance of treatment and genotypic effects on expression levels of individual mirnas. instead, we viewed each mirna as an independent contributor to the overall mirna expression profile of each genotype under each experimental condition, and we report patterns of expression among groups. in qpcr experiments, groups were compared using two-way anova, with gonadal sex and sex chromosome complement as independent factors. values for mir- and mir- were log10-transformed before statistical analysis.  study design to identify sex differences in adipose tissue mirnas under three physiological conditions to assess the effect of gonadal hormones and sex chromosomes on mirna levels in adipose tissue, we sequenced mirnas in the gonadal fat depot of fcg mice. the four genotypes of this model allow analysis of sex differences based on the presence of male vs. female gonads, or the presence of xx vs. xy chromosome complement. additionally, we analyzed the effect of circulating gonadal hormones by comparing changes in mirna profiles in mice with intact gonads vs. mice after gonadectomy. finally, we assessed the effect on adipose tissue mirna profile of high fat diet, and its interaction with sex, by comparing mice on chow and high fat diets. adipose tissue samples from the gonadal fat depot were collected from all cohorts at the same age. sequenced and mapped mirnas were counted to establish mirna expression profiles.fig. study design to identify sex and dietary effects on mirna expression. a mirnas were sequenced in adipose tissue from the four core genotypes mouse model. a comparison of mirna levels in female xx and xy adipose tissue with those in male xx and xy tissue allows detection of effects due to gonadal type. a comparison of mirna levels in female and male xx adipose tissue with those on female and male xy tissue allows detection of effects due to sex chromosome type. b to determine the effects of acute gonadal hormones of diet, mirna sequencing was performed in four core genotypes cohorts that were fed a chow diet in the gonadally intact state, fed a chow diet and gonadectomized as adults to remove the acute effects of gonadal secretions, and in mice that were gonadectomized as adults and fed a high fat diet for 16 weeks sex differences in adipose tissue mirna profiles we sequenced and mapped mature and precursor mirnas in gonadal adipose tissue of fcg mice. we identified mirnas that were expressed at substantial levels in gonadal fat. three mirna species—mir-10b-5p, mir-143-3p, and mir-22-3p—accounted for nearly 50% of all mirna reads. an additional mirna species each accounted for at least % of the gonadal adipose tissue mirna pool. of the most prevalent mirnas in mouse adipose tissue have previously been shown to be abundant in human subcutaneous white adipose tissue.fig. highly expressed mirnas in gonadal fat. a out of mapped mirnas, were expressed at over counts per million reads in at least one of the twelve libraries. these mirnas are termed “highly expressed” mirnas. b percent abundance of mirnas in gonadal fat of fcg mice. values represent mean ± sd of the twelve sequencing libraries. c, d three individual mirna sequencing libraries representing biological replicates of xx females fed a high fat diet are correlated with the pooled sequencing library composed of those same mice. highly expressed gonadal fat mirnas are more correlated between the individual and pooled libraries than mirnas that are expressed at a lower level. values represent mean ± sd. r and p-values were calculated using pearson product–moment correlation to streamline our sequencing process for the twelve conditions shown in fig.  we assessed the feasibility of pooling three biological replicates for mirna sequencing for each condition. we sequenced samples from three individual mice separately and after pooling together. we found that the average mirna reads of the three individual samples were highly correlated with the values from the pooled sample. based on these findings, we concluded that pooled libraries were a satisfactory representation of individual biological samples. for analyses described below, we focused on the most abundantly expressed adipose tissue mirnas, as these are most likely to have physiological significance. to establish a baseline profile of sex effects on mirna expression in adipose tissue, we sequenced mirnas from mice fed a standard chow diet. we first asked whether there was a male or female bias in overall mirna expression by examining the ratio of female-to-male mirna levels. female-to-male ratios of each mirna were log2-transformed and plotted as a histogram. adipose tissue mirnas had a significant bias toward higher expression in male compared to female gonadal adipose tissue. in the chow-fed cohort, there was no significant bias in mirna distribution between xx and xy mice, although we cannot exclude the possibility that the sex chromosome complement may alter expression levels of individual mirna species. thus, in the basal physiological condition of chow-fed animals, the presence of male vs. female gonads led to a significant difference in adipose tissue mirna expression profile.fig. distribution of mirnas in chow-fed gonadally intact mice. distribution of female-to-male ratios and xx-to-xy ratios in highly expressed mirnas. ratios were calculated for each mirna, log2-transformed, and binned in increments. values less than − were grouped into the first bin, and values greater than were grouped into the last bin. red bars refer to male -biased mirnas, black bars refer to female -biased mirnas, blue bars refer to xy-biased mirnas, and white bars refer to xx-biased mirnas. p-values were calculated using the exact binomial test adipose tissue mirnas regulated by gonadal hormones we hypothesized that male–female differences in mirna expression profile were caused by differences in gonadal hormones. gonadal hormones cause sex differences in two primary ways. first, they impose permanent, or “organizational,” effects during development, such as the effect of testosterone on male genital differentiation. second, male and female gonadal hormones in the circulation have acute/reversible effects on processes such as gene expression. to investigate the effects of acute gonadal hormone action on mirna expression, we gonadectomized adult fcg mice and assessed mirna profiles 5 months later, during which time circulating gonadal hormones were absent. with this experimental design, male–female differences that disappear after gonadectomy can be attributed to acute hormone effects. in gonadectomized mice, the mirna profile expressed as female-to-male ratio exhibited a small bias toward higher mirna expression levels in females compared to males. comparison of female–male ratio distributions between gonadectomized mice and gonadally intact mice revealed a significant difference. this represented a shift from a profile with male > female expression ratios in gonadally intact mice to a profile with a slight bias toward female > male ratios in gonadectomized mice. the mirna ratio distribution of xx vs. xy mice was similar in intact and gonadectomized mice. taken together, these data suggest that acute effects of gonadal secretions, but not sex chromosome complement, produce sex biases in overall mirna levels in mice fed a chow diet. we note, however, that there are subsets of mirnas with extreme values in female-to-male ratios and in xx-to-xy ratios in gonadectomized mice, suggesting that, in the absence of circulating gonadal hormones, sex and sex chromosome differences remain in individual mirnas.fig. gonadectomy alters distribution of mirnas in chow-fed mice. distribution of female-to-male ratios and xx-to-xy ratios in highly expressed mirnas. ratios were calculated for each mirna and log2-transformed in gonadally intact and gonadectomized fcg mice fed a chow diet. data for gonadally intact mice are identical to fig. 3a–b, but binned in increments and compared to data distribution for gdx mice. each dot represents one of the highly expressed mirnas. black bars represent median values. distributions were significantly different for female-to-male ratios, but not for xx-to-xy ratios. p-values were calculated using the wilcoxon rank sum test. f, female; m, male adipose tissue mirna sex differences influenced by high fat diet or obesity several studies have shown that mirna expression levels are regulated during adipogenesis or are altered in obesity. however, it is unknown whether sex influences mirna expression levels in obese adipose tissue. to identify sex and/or sex chromosome factors that alter mirna expression in obesity apart from the acute effect of gonadal hormones, we performed mirna-seq on gonadal fat from fcg mice that were gonadectomized and then fed a high fat diet beginning 1 month later for 4 months. the distribution of female-to-male mirna read count ratios was significantly male-biased. this male bias in gonadectomized high fat diet-fed mice represents a significant shift in mirna expression compared to the distribution in chow-fed mice that had been gonadectomized at the same age and therefore represents an effect of diet. interestingly, mirna expression levels were also biased toward the xy sex chromosome complement in mice fed a high fat diet. this indicates that the sex chromosome complement, in addition to gonadal hormones, influenced mirna expression in mice fed a high fat diet. comparison of the distributions of xx-to-xy ratios indicated a shift from a slight xx bias in chow-fed mice to an xy bias in fat-fed mice. in addition, there are a number of mirnas that shift to the extreme ends of the sex chromosome distribution in mice fed a high fat diet. these data suggest that, in gonadectomized mice, high fat diet and obesity modulate mirna expression in a sex- and sex chromosome-specific manner.fig. high fat diet alters mirna expression levels in a sex-dependent and sex chromosome complement-dependent manner. distribution of female-to-male ratios and xx-to-xy ratios in highly expressed mirnas. ratios were calculated for each mirna and log2-transformed in gonadectomized fcg mice fed a high fat diet. data for chow-fed gonadectomized mice are identical to fig. 4a–b. distributions were significantly different for female-to-male ratios and for xx-to-xy ratios. p-values were calculated using the wilcoxon rank sum test. f, female; m, male we followed up on select mirnas that had distinct sex-influenced expression patterns and/or were previously implicated in adipocyte biology. quantitative real-time pcr of mir-133a-3p, previously shown to inhibit brown adipocyte differentiation, was upregulated in xx compared to xy mice in the chow-fed, gonadally intact group. this difference was abolished after gonadectomy, suggesting an interplay between gonadal hormones and sex chromosome complement. mir-192-5p levels were higher in males compared to females in gonad-intact mice, but this sex difference is lost after gonadectomy, and an xx > xy difference emerges. conversely, mir-205-5p showed a sex-chromosome bias in gonadally intact mice and a sex bias in gonadectomized mice fed a chow diet. however, when mice were fed a high fat diet, mir-192-5p and mir-205-5p showed no sex or sex chromosome differences. mir-221-3p, which is located on the x chromosome and has been previously associated with obesity and adipogenesis, showed no sex or sex chromosome difference in mice fed a chow diet. however, mice fed high fat diet showed increased levels of mir-221-3p in xx male and female mice compared to xy mice. these mirna species illustrate several distinct patterns of expression that are dependent on gonadal type, gonadal hormones, diet, and sex chromosome complement, and emphasize that sex differences exist in mirnas that have been studied previously in the context of adipose tissue and obesity. to examine the relationship between diet, acute hormone effects, gonadal type, and sex chromosome complement effects on mirna expression, we performed principal component analysis. diet distinguished mirna expression levels to the greatest degree, explaining % of mirna covariance. gonadal state and gonadal type led to a milder but detectable separation. the sex chromosome complement did not lead to divergence of mirna covariance. because diet had such a pronounced effect on mirna covariance, we performed principal component analysis on mice fed only the chow diet to determine if sex-biasing factors still influenced mirna covariance. the first two components revealed deviation of gonadally intact mice from gonadectomized mice. in addition, when examining the third and fourth principal components, females clustered separately from males, and xx mice exhibited a slight degree of separation from xy mice. we note that the effect of the sex chromosome complement was detectable after removing the overwhelming effect of diet. taken together, these data suggest that diet is a major factor in mirna expression levels, and that sex factors, such as acute gonadal secretions, male–female gonads, and xx–xy chromosome complement, also play important roles in modulating mirna levels.fig. principal component analysis reveals factors influencing mirna correlation covariance. pca of all twelve sequencing libraries or chow-fed samples only. each dot represents a sequencing library and consists of normalized read counts for each of the highly expressed mirnas. dots are colored and encircled according to diet, gonadal state, gonadal sex, or sex chromosome complement. hfd, high fat diet; int, gonadally intact; gdx, gonadectomized; f, female; m, male discussion mirnas are key regulators of gene expression through their effects on the degradation or translation of mrna transcripts. in metabolic tissues, mirnas influence gene expression related to adipogenesis, cholesterol homeostasis, and glucose homeostasis. these metabolic processes are also affected by sex. while there is evidence for sexual dimorphism in mirna expression in brain, skeletal muscle, and immune cells, sex effects on mirnas in adipose tissue have not been characterized. in this study, we used the fcg mouse model to identify sex chromosome complement and gonadal hormone effects on mirna abundance in gonadal fat. we sequenced mirnas in three distinct cohorts of fcg mice in different hormonal and dietary conditions to determine the effects of circulating gonadal hormones and high fat diet on sex differences in gonadal fat mirnas. we selected the gonadal fat depot for our studies because it has been widely studied in the mouse obesity field. however, it is clear that differences exist between gonadal fat and other visceral fat depots, as well as subcutaneous fat depots, and it would ultimately be of interest to delineate sex differences across these distinct depots. our analyses focused on the most abundant mirnas, as we suspect that sex differences in their levels are the most likely to have significant physiological effects. by comparing overall mirna profiles across diet, gonadal state, gonadal type, and sex chromosome complement, we determined that each of these factors influence mirna expression in white adipose tissue. in gonadally intact mice, we identified a male bias in mirna expression levels and no significant bias between xx and xy mice. as expected, the distribution of male–female differences in the levels of mirnas was altered by gonadectomy, suggesting that circulating gonadal hormones regulate mirna levels and are responsible for some male–female sex differences. gonadectomy also influences adiposity, which may impact mirna expression. we have previously shown that gonadectomized xx mice have increased body weight and adiposity compared to gonadectomized xy mice. because we studied mice 5 months after gonadectomy, we cannot distinguish whether changes in mirna expression preceded increases in adiposity or resulted from them. however, male–female differences in mirna levels do not correlate with xx–xy differences in body weight, suggesting that sex differences in adiposity do not determine sex differences in overall mirna levels. studies measuring mirna levels after gonadectomy but before weight gain will be required to address these questions. high fat diet feeding alters sex and sex chromosome differences in mirna levels, suggesting a diet–sex interaction in mirna expression. it is notable that the observed sex and sex chromosome biases result from changes in overall mirna distributions and not from changes in a few specific mirnas. high fat feeding revealed a sex chromosome bias in mirna expression ratios that was not apparent in chow-fed animals. this raises the possibility that differential mirna expression levels may contribute to the dramatic difference in adiposity of xx compared to xy mice in response to high fat diet. thus, mice with two x chromosomes gain more fat than xy mice, and this difference is enhanced upon feeding a high fat diet. further studies will be required to elucidate the physiological roles of individual sex-biased mirnas in adipogenesis and diet-induced obesity. some mirnas with distinct sex-biased patterns in expression, such as mir-133a and mir- have been implicated in adipogenesis. other mirnas may have yet unidentified roles in lipid metabolism. an initial analysis of potential target genes for sex-biased mirnas that we identified revealed enrichment in specific metabolic pathways. these predictions suggest that pathways such as mtor signaling and phosphatidylinositol signaling may be important in mediating sex differences in healthy and/or diseased adipocytes. studies to determine the impact of these predicted target genes and corresponding pathways on sex differences in obesity provide a unique area for future research.  this study represents a first look at mechanisms underlying sex differences in mirna profiles in adipose tissue. we demonstrate that mirnas in adipose tissue are influenced by diet, gonadal state, gonadal type, and sex chromosome complement. principal component analysis revealed that diet is the leading factor in mirna expression covariance, while gonadal state, gonadal type, and sex chromosome complement were milder but important regulators. by recognizing the innate sex differences in mirna levels, we can better understand sex differences in fat accumulation, distribution, and adipose pathophysiology. additional files additional file 1: mature mirna sequences that mapped uniquely to the mirbase additional file 2: mirnas that mapped to unique sequences. additional file 3: cases of mirnas that were of unequal length, and thus read counts had to be adjusted for non-unique mapping. additional file 4: count table summarizing read counts for mature and precursor mirna sequences in fifteen sequencing libraries. perfectly matched counts are normalized to total library reads per million. additional file 5: figure s expression of select mirnas with distinct sex- and sex-chromosome-biased patterns, measured by taqman qpcr. table s enrichment of specific cellular pathways in predicted targets of sex-influenced mirnas. abbreviations fcgfour core genotypes mir, mirnamicrorna mrnamessenger rna 
this model, which consists of xx female, xx male, xy female, and xy male mice, allowed us to assess independent effects of gonadal type and sex chromosome complement on mirna expression profiles. we identified a male–female effect on the overall mirna expression profile in mice fed a chow diet, with a bias toward higher expression in male compared to female gonadal adipose tissue. this sex bias disappeared after gonadectomy, suggesting that circulating levels of gonadal secretions modulate the mirna expression profile. our results demonstrate that the overall mirna expression profile in adipose tissue is influenced by gonadal hormones and the sex chromosome complement, and that expression profiles change in response to gonadectomy and high fat diet. differential mirna expression profiles may contribute to sex differences in adipose tissue gene expression, adipose tissue development, and diet-induced obesity. principal component analysis revealed that high fat diet has a substantial effect on mirna profile variance, while gonadal secretions and sex chromosome complement each have milder effects. specific mirnas are key regulators of adipogenesis, and are expressed at different levels in adipose tissue from lean and obese mice. to identify sex differences in adipose tissue mirna expression and to understand the underlying mechanisms, we performed high-throughput mirna sequencing in gonadal fat depots of the four core genotypes mouse model. micrornas are small non-coding rna molecules that regulate gene expression by targeting specific mrna species for degradation or interfering with translation. the degree of lipid accumulation and distribution of white adipose tissue differs between males and females, and it is unknown whether sex differences in adipose tissue-specific mirna expression may contribute to this dimorphism. we have also assessed the effects of a high fat diet on sex differences in adipose tissue mirna profiles. after 16 weeks of high fat diet, the mirna expression distribution was shifted toward higher expression in xy vs. xx adipose tissue. micrornaadipose tissuesex differencessex chromosomehttp://dx.doi.org/ /100006955office of extramural research, national institutes of healthgm007185link jenny c. http://dx.doi.org/ /100000062national institute of diabetes and digestive and kidney diseasesr56dk083561r56dk083561arnold arthur p. reue karen http://dx.doi.org/ /100000050national heart, lung, and blood institutep01hl28481p01hl28481hasin-brumshtein yehudit lusis aldons j. issue-copyright-statement© the author 2017 however, the sex chromosome complement is also a determinant of sex differences and may regulate mirna expression in adipocytes. typically, sex differences are attributed to hormones secreted from ovaries or testes.
 micrornas are small non-coding rna molecules that regulate gene expression by targeting specific mrna species for degradation or interfering with translation. specific mirnas are key regulators of adipogenesis, and are expressed at different levels in adipose tissue from lean and obese mice. the degree of lipid accumulation and distribution of white adipose tissue differs between males and females, and it is unknown whether sex differences in adipose tissue-specific mirna expression may contribute to this dimorphism. typically, sex differences are attributed to hormones secreted from ovaries or testes. however, the sex chromosome complement is also a determinant of sex differences and may regulate mirna expression in adipocytes.  to identify sex differences in adipose tissue mirna expression and to understand the underlying mechanisms, we performed high-throughput mirna sequencing in gonadal fat depots of the four core genotypes mouse model. this model, which consists of xx female, xx male, xy female, and xy male mice, allowed us to assess independent effects of gonadal type and sex chromosome complement on mirna expression profiles. we have also assessed the effects of a high fat diet on sex differences in adipose tissue mirna profiles. we identified a male–female effect on the overall mirna expression profile in mice fed a chow diet, with a bias toward higher expression in male compared to female gonadal adipose tissue. this sex bias disappeared after gonadectomy, suggesting that circulating levels of gonadal secretions modulate the mirna expression profile. after 16 weeks of high fat diet, the mirna expression distribution was shifted toward higher expression in xy vs. xx adipose tissue. principal component analysis revealed that high fat diet has a substantial effect on mirna profile variance, while gonadal secretions and sex chromosome complement each have milder effects.  our results demonstrate that the overall mirna expression profile in adipose tissue is influenced by gonadal hormones and the sex chromosome complement, and that expression profiles change in response to gonadectomy and high fat diet. differential mirna expression profiles may contribute to sex differences in adipose tissue gene expression, adipose tissue development, and diet-induced obesity. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords micrornaadipose tissuesex differencessex chromosomehttp://dx.doi.org/ /100006955office of extramural research, national institutes of healthgm007185link jenny c. http://dx.doi.org/ /100000062national institute of diabetes and digestive and kidney diseasesr56dk083561r56dk083561arnold arthur p. reue karen http://dx.doi.org/ /100000050national heart, lung, and blood institutep01hl28481p01hl28481hasin-brumshtein yehudit lusis aldons j. issue-copyright-statement© the author 2017 