diabetes mellitus is a chronic disease that has become a major health and economic problem nowadays .
387 million people in the world suffer from diabetes and it is forecast that the number will be 205 million more in 2035 .
the state expenses incurred on the disease include not only medications , hospital inpatient care , and medical visits but also indirect costs like absenteeism , reduced productivity , and inability to work . because of that
, more attention is paid to prophylaxis of obesity as well as to optimization and individualization of diabetes treatment .
type 2 diabetes mellitus ( t2 dm ) is determined by genetic and environmental causes . the risk factors contributing to the disease are , inter alia , age , body mass index , physical inactivity , and sedentary lifestyle [ 4 , 5 ] . as far as genetic factors predisposing to t2 dm
are concerned , genetic variants modifying the risk of development of the disease have been identified ; however , altogether , they explain only a small proportion of the genetic background of t2 dm [ 6 , 7 ] .
among possible causes of  missing heritability  gene , environment interactions are mentioned . in recent years
, there have been many studies on the relationship between the composition of the gastrointestinal bacterial flora , t2 dm , and obesity [ 911 ] .
these analyses indicate that the etiology of the various forms of diabetes , especially those associated with obesity and insulin resistance , such as t2 dm , may be related to intestinal microflora [ 9 , 10 ] . on the other hand , monogenic forms of the disease , like maturity onset diabetes of the young ( mody ) , are determined mostly by genetic factors , with minor effect of the environment .
hnf1a - mody is one of the most common forms of subtypes of mody [ 12 , 13 ] .
it is caused by mutation in hnf1a ( hepatic nuclear factor 1 alpha ) gene .
hnf1a - mody is characterized by progressive hyperglycemia , with the onset of diabetes in youth .
the differences in the gut microbiota between patients with t2 dm and hnf1a - mody in relation to healthy individuals could be a new argument for an influence of microbiota as an environmental risk factor for some forms of diabetes .
the gut microbiota perform multiple functions , produce vitamins , inhibit the growth of potential pathogenic bacteria , and lower the cholesterol level in blood .
the bacterial metabolite indole can modulate the secretion of glucagon - like peptide-1 ( glp-1 ) by l - cells .
it is hypothesized that remission of diabetes after bariatric surgery is connected with neurohormonal signaling or changes in the ecosystem of microbiota of the gastrointestinal tract [ 18 , 19 ] . the intestinal microbiota play a role in the circulation of bile acids .
it is suggested that the increased levels of bile acids after bariatric operation stimulate glp-1 secretion and lead to improvement in glucose metabolism .
the nuclear farnesoid x receptor ( fxr ) is also the molecular target for metabolic condition after bariatric surgery .
the new approach is to assess the bacterial flora through the use of the noncultivated methods allowed for identifying bacteria that are difficult to culture .
these include targeted gene sequencing , genome sequencing , and shotgun metagenomic sequencing . 
 the aim .
the aim of the investigation was to compare the qualitative parameters of the colonic flora in patients with hnf1a - mody diabetes , matched controls , and patients with t2 dm in a single stool sample .
the study was designed in collaboration with the department of metabolic diseases , jagiellonian university medical college , krakw , poland ; university hospital , krakw , poland ; and center for medical genomics omicron and department of microbiology , jagiellonian university medical college , krakw , poland .
the project has obtained the consent of the ethics committee of the jagiellonian university number kbet/272/b/2013 .
participants were the patients of the department of metabolic diseases and volunteers , all of whom consented to participate in the research .
we collected stool samples from ten patients with hnf1a - mody and twenty - three patients with t2 dm .
all participants declared not using antibiotics 4 weeks prior to stool sample delivery and they denied having met the exclusion criteria ( see the following list ) . 

exclusion criteria 
 patients who did not agree to participate in the study and who withdrew during the investigation  patients below 18 years of age and over 65 years of age  patients taking antibiotics for up to 30 days before giving the sample of stool  patients using probiotics  patients with confirmed infection of the gastrointestinal tract  patients with chronic inflammatory bowel disease of unknown etiology  patients with active cancer ( especially of the gastrointestinal tract )  patients with immunodeficiency  patients with features of liver damage ( with the exception of nonalcoholic fatty liver transaminase levels less than three times the upper limit of the normal level ) 
  patients who did not agree to participate in the study and who withdrew during the investigation patients below 18 years of age and over 65 years of age patients taking antibiotics for up to 30 days before giving the sample of stool patients using probiotics patients with confirmed infection of the gastrointestinal tract patients with chronic inflammatory bowel disease of unknown etiology patients with active cancer ( especially of the gastrointestinal tract ) patients with immunodeficiency patients with features of liver damage ( with the exception of nonalcoholic fatty liver transaminase levels less than three times the upper limit of the normal level )   in most patients , blinded continuous glucose monitoring ( cgm , ipro2 , medtronic ) was implemented .
in addition , we utilized results of stool sample analysis of healthy individuals done previously ( project supported by the national science centre in poland , number dec-2011/03/d / nz5/00551 ) .
one control group of healthy individuals was matched with mody patients with respect to age and bmi ( body mass index ) , and the second one was matched to t2 dm patients with respect to age .
the bacterial dna was isolated using genomic mini ax stool spin ( a&a biotechnology ) modified to include enzymatic treatment ( lysozyme , lysostaphin , and lyticase ) and bead - beating step .
libraries were prepared according to illumina 16s metagenomic sequencing library preparation protocol ( https://support.illumina.com/downloads/16s_metagenomic_sequencing_library_preparation.html ) . in short ,
after the polymerase chain reaction ( pcr ) clean - up , samples were indexed , cleaned , and pooled .
10  pm libraries with 10% phix spike - in were sequenced on illumina miseq using v3 sequencing kit ( 300 
resulting short reads were joined on overlapping regions using fastq - join tool from ea - utils package .
otus ( operational taxonomic units ) were picked using open - reference protocol . in the first step , closed reference otu
subsequent reads that failed to hit the reference database were filtered out and used to perform de novo otu picking .
relative otu abundance was calculated using qiime ( http://qiime.org/ ) . to estimate alpha diversity chao 1 , the observed otus and phylogenetic
distance metrics were calculated . both weighted and unweighted unifrac distances were calculated to analyze beta diversity .
frequency of otu across sample groups was compared using nonparametric anova ( kruskal - wallis test ) and mann - whitney u test .
post hoc tests using statistica were performed if there were significant differences between groups . to correlate clinical parameters with relative otu abundance , spearman 's rho test was applied .
to compare the clinical parameters , the adequate parametric ( t - test ) and nonparametric ( mann - whitney u test , kruskal - wallis test ) tests were used . as diet
is concerned at disease diagnosis , all mody and t2 dm patients were educated to follow the diet recommended in diabetes with 4050% of calories coming from carbohydrates , 2030% from fat , and 20% from protein .
there were 23 t2 dm patients ( 13 women and 10 men ) and 10 hnf1a - mody individuals ( 5 women and 5 men ) included in the study .
the median age in the t2 dm group was 60 ( 5662 ) years and in the hnf1a - mody group was 36.5 ( 3056 ) , p = 0.0133 .
the mean duration of diabetes in t2 dm group was 4.46  2.95 years versus 19.1  13.37 years , p < 0.0001 .
the median hba1c ( glycosylated haemoglobin ) level in t2 dm group was 8.12% ( 7.11%8.78% ) versus 7.205% ( 5.4%7.86% ) , p = 0.0479 , respectively .
the first control group , matched for age ( p = 1.0 ) and bmi ( p = 0.26 ) with hnf1a - mody group , consisted of 16 healthy individuals ( 11 women and 5 men ) with median age of 39.5 ( 31.549.5 ) years and median bmi of 23.77 ( 22.8524.95 )  kg / m .
the median bmi in hnf1a - mody group was 25.74 ( 24.2229 )  kg / m and in t2 dm group was 30.25 ( 27.6833.25 )  kg / m , p = 0.0290 .
the t2 dm group differed in age ( p = 0.0006 ) and bmi ( p < 0.0001 ) in comparison to the first control group ( figure 1 ) .
the second control group consisted of 10 people ( 7 women and 3 men ) .
the median bmi in the second control group was 25.57 ( 24.3227.38 )  kg / m , and the median age was 56 ( 53.558 ) years .
there was no difference in age between this control group and t2 dm group ( p = 0.3059 ) , but the groups differ in bmi ( p = 0.0012 ) .
based on the questionnaire data , we found no major differences in the dietary profile between t2 dm and mody patients .
the questionnaire included four possible answers regarding following the prescribed diet : always , most of the time , rarely , or never .
the majority ( 16 out of 23 patients with t2 dm and 8 out of 10 individuals with mody ) declared following always or most of the time the diet recommended in diabetes with respect to the proportion of calories coming from carbohydrates , fat , and protein and with respect to food quality . unfortunately , we had no reliable data concerning diet in control group .
the patients ' treatment included metformin , sulfonylurea , acarbose , and insulin ( supplementary material available online at http://dx.doi.org/10.1155/2016/3876764 ) .
the analysis of the samples by sequencing obtained an average of 180353.52 paired reads per sample , with median 140780 and standard deviation 152430.2 .
after preprocessing ( removing adapters , pcr primers , filtration of low quality bases , pair joining , and removal of chimeric reads ) , we received an average of 115782.34 readings per sample , with median 87994 and standard deviation 104992.4 .
the bacterial composition was analyzed at the phylum , class , order , family , and genus level .
the bacteria and archaea profiles for the first control , hnf1a - mody , and t2 dm at the phylum level were as follows : euryarchaeota 0.12% versus 0.11% versus 0.13% ; actinobacteria 9.53% versus 9.31% versus 11.70 ; bacteroidetes 6.39% versus 1.09% versus 1.30% ; firmicutes 77.02% versus 87.36% versus 80.09% ; proteobacteria 0.36% versus 0.35% versus 2.82% ; synergistetes 0.00% versus 0.00% versus 0.10% ; tenericutes 0.05% versus 0.005% versus 0.001% ; verrucomicrobia 6.52% versus 1.78% versus 3.85% ( figure 2 ) .
there were significant differences between t2 dm and the first control groups at the phylum level ( kruskal - wallis test with a false discovery rate ( fdr ) correction ) . in t2
dm group , in comparison to the first control group , the amount of proteobacteria was higher ( p = 0.0006 ) and the amount of bacteroidetes was lower ( p = 0.0005 )
. in hnf1a - mody group , the frequency of bacteroidetes was lower than in the control group ( p = 0.0143 ) .
there were no significant differences between t2 dm and hnf1a - mody groups at the phylum level .
the firmicutes / bacteroidetes ( f / b ) ratio was significantly higher in both t2 dm and hnf1a - mody groups than in the first control group ( figure 3 ) .
the firmicutes / bacteroidetes ratio was correlated positively with the hba1c level in hnf1a - mody group ; we did not observe such phenomenon in t2 dm .
there were significant differences at the class level : the proportion of flavobacteria and bacteroidia was higher in the first control group than in t2 dm and hnf1a - mody groups .
the amount of gammaproteobacteria was higher in t2 dm group than in the first control group ( p = 0.0002 ) .
the amount of erysipelotrichi was higher in hnf1a - mody group than in the first control group and t2 dm group ( supplementary data ) .
there were significant differences at the order level : flavobacteriales and bacteroidales were higher in the first control group than in hnf1a - mody and t2 dm groups .
actinomycetales was lower in the first control group than in hnf1a - mody and t2 dm groups .
turicibacterales was higher in the first control group than in t2 dm group ( p = 0.0073 ) , and also turicibacterales was higher in hnf1a - mody group than in t2 dm group .
erysipelotrichales was higher in hnf1a - mody group than in the first control group and in t2 dm group ( supplementary data ) .
there were significant differences at the family level : flavobacteriaceae and bacteroidaceae were higher and promicromonosporaceae and the other family of actinomycetales were lower in the first control group than in both t2 dm and hnf1a - mody groups .
enterobacteriaceae was higher and porphyromonadaceae was lower in t2 dm group than in the first control group ( p = 0.004 ) .
the frequency of clostridiaceae and turicibacteraceae was lower in t2 dm group than in the first control group and hnf1a - mody group .
the amount of erysipelotrichaceae was higher in hnf1a - mody group than in the first control group and t2 dm group ( supplementary data ) .
there were significant differences at the genus level : bacteroides and the unnamed genus of flavobacteriaceae were higher and cellulosimicrobium and the other genus of the other family of actinomycetales were lower in the first control group than in t2 dm and hnf1a - mody groups .
the amount of parabacteroides , unnamed genus of clostridiaceae , turicibacter , lachnospira , and anaerostipes was higher , and the unnamed genus of enterobacteriaceae was lower in the first control group than in t2 dm group .
the frequency of the other genus of clostridiaceae and smb53 was lower in t2 dm group than in the first control group and hnf1a - mody group ( supplementary data ) .
comparing the microbiome between t2 dm group and the second control group ( matched with respect to age ) revealed the differences in the following otus : the frequency of proteobacteria ( phyla level ) , actinomycetales ( order level ) , the other family of actinomycetales , promicromonosporaceae ( family level ) , the other genus of the other family of actinomycetales , cellulosimicrobium , bulleidia , and eubacterium ( genus level ) was higher and the frequency of clostridiaceae , ruminococcaceae ( family level ) , and the unnamed and the other genus of clostridiaceae and smb53 ( genus level ) was lower in t2 dm group than in the second control group ( supplementary data ) .
the observed otus alpha diversity ( p = 0.027 ) but not chao 1 ( p = 0.063 ) was significantly lower in t2 dm than in the control group ; there was no difference in alpha diversity between control and hnf1a - mody , hnf1a , and t2 dm groups ( figure 4 ) .
the beta diversity analysis revealed the differences in communities structure between t2 dm and the second control group .
we obtained data of continuous glucose monitoring from 21 patients with t2 dm and 9 patients with hnf1a - mody .
the mean glucose level in t2 dm group was 167.33  35.34  mg / dl , and the mean glucose level in hnf1a - mody group was 143.11  38.07 
mg / dl ; the groups did not differ in mean glucose level ( t - test , p = 0.1037 ) .
the percentage of time with glucose level above 140  mg / dl was higher in t2 dm group ( mann - whitney u test , p = 0.0235 ) ; the percentage of time with glucose level within 70140  mg / dl was higher in hnf1a - mody group ( mann - whitney u test ; 29% ( 14%54% ) versus 49% ( 32%81% ) , p = 0.0352 ) ; the percentage of time with glucose level below 70  mg / dl was higher in hnf1a - mody group ( mann - whitney u test ; 2% ( 0%6% ) versus 0% ( 0%-0% ) , p = 0.0127 ) .
this study for the first time compares the bacterial flora of patients with hnf1a - mody with the bacterial flora of individuals with t2 dm and control group .
herein , we have shown for the first time significant differences in the bacteria profile of stool samples between hnf1a - mody and both control and t2 dm groups .
it is true that the diabetes subpopulations differed at the baseline in age , duration of diabetes in years , bmi , and hba1c .
it is also possible that other confounders like diet   or medications [ 2527 ] could influence the results .
some differences between the control group and diabetes cohorts can result from medications including antidiabetic drugs [ 26 , 27 ] .
it was demonstrated that the f / b ratio depends on age . in infants and elders , it is significantly lower than in adults .
it has also been presented that the obese people have increased f / b ratio .
t2 dm and hnf1a - mody groups in comparison with control group share some between - group differences but the study also revealed some type of diabetes - related differences . between diabetes groups in otu taxa
, we observed that the frequency of turicibacterales was significantly higher in hnf1a - mody group than in t2 dm group and in the first control group than in t2 dm group .
in the study results , the amount of enterobacteriales was statistically higher in t2 dm group than in the first control group .
we can not forget that only hnf1a - mody group was matched with respect to age and bmi with the first control group .
the baseline differences between t2 dm and hnf1a - mody group could be prone to and mask some differences in the microbiome profile between diabetes groups . nevertheless , the study showed differences in the bacterial profile between hnf1a - mody and control groups
. this can be an argument for the thesis that also hnf1a - mody is associated with changes in colonic bacterial flora .
the differences in microbiome profiles between the first control group and t2 dm group were not the same as between the second control group and t2 dm group .
this could result from the differences in age ; only the second control group was matched with respect to age with the t2 dm group .
still , the second control group was not matched with respect to bmi with the t2 dm group . in recent studies
, it has been revealed that the alpha diversity ( within samples diversity ) is decreased in obese individuals , but no significant differences were found in alpha diversity between healthy people and diabetes patients [ 30 , 31 ] . in our study , we observed lower otus alpha diversity in t2 dm group than in the control group , but the groups differ in bmi .
we did not observe differences in alpha diversity between hnf1a - mody and both control and t2 dm groups .
nevertheless , the alpha diversity in hnf1a - mody group was lower than in the control group and higher than in the t2 dm group , but the differences were not significant .
one could suggest that the colonic microbiome plays a different role in the metabolic control in t2 dm in comparison to hnf1a - mody .
one could speculate that manipulation in the structure of gut microbiota , in order to improve glucose metabolism in patients with diabetes , could be more effective in some types of diabetes ( like t2 dm ) than in others and should be individualized to the type or subtype of diabetes . some bacteria that belong to the firmicutes phylum are butyrate producers .
they play a key role in human gut energy supply by producing butyrate , the main source of energy for the colonic epithelium .
furthermore , butyrate improves insulin sensitivity   and is able to trigger the secretion of glp-1 from l - cells .
demonstrated that short fatty acids stimulate the secretion of glp-1 in vitro . also , some suggest that there is a relationship between the short fatty acid as histone deacetylases inhibitors and the role of epigenetics in the development of diabetes .
it has to be emphasized that , in recent studies , it is reported that the amount of butyrate producers such as roseburia species and faecalibacterium prausnitzii is lower in patients affected with t2 dm [ 11 , 36 ] . in our data ,
the frequency of the genus faecalibacterium was 0.5% in control , 0.2% in hnf1a - mody , and 0.5% in t2 dm group but the differences did not approach a level of significance ( p > 0.05 ) .
the genus roseburia was abundant in 0.5% in control , in 1.1% in hnf1a - mody , and in 0.3% in t2 dm ( p > 0.05 for each comparison ) . in our study , the relative abundance of akkermansia tended to be higher in the first control group ( 6.5% ) than in hnf1a - mody ( 1.8% ) and t2 dm ( 3.9% ) groups .
it has been shown that in humans the percentage of akkermansia muciniphila in the healthy population is 14% [ 37 , 38 ] .
the abundance of the bacteria was reported to be decreased in genetically ( t2 dm ) and diet - induced obesity mice model . on the other hand
, it has been shown that in t2 dm microbiota genes belonging to akkermansia are increased . of note ,
four individuals in the t2 dm group had very high relative abundance of akkermansia : 31.4% , 32.1% , 21.1% , and 20.6% .
three controls from the first control group ( 29.2% , 17.8% , and 14.3% ) and two hnf1a - mody patients ( 9.5% and 13.7% ) were also characterized by high frequency of akkermansia .
it was demonstrated that higher amount of akkermansia muciniphila in overweight and obese people resulted in healthier metabolic status and better glucose homeostasis after introducing calorie restriction diet .
definitely , one of the limitations of our study is due to the age differences between the t2 dm and mody individuals .
the difference is due to the fact that mody individuals are usually diagnosed at much younger age in comparison to patients with type 2 diabetes .
the individuals from older generations representing mody families were usually not available for genetic testing ( either death or no willingness for testing ) . on the other hand , hardly ever do we diagnose type 2 diabetes in patients below 50 years of age .
this is why matching t2 dm and mody patients for the age was not possible . to address these differences
we failed , however , to match controls with t2 dm as far as bmi is concerned ( nondiabetic obese individuals were not willing to follow the study protocol ) .
the estimation was based on a self - reported questionnaire and was not verified ; thus , we can not fully exclude the notion that the results of our study to some degree were affected by differences in dietary patterns . however , the major limitation of our study was the relatively small sample size , which was partially due to protocol requirements .
comparison of bacterial profiles between patients with hnf1a - mody and individuals with mutation in hnf1a gene without clinical diabetes could be significant .
it could be suggested that all of the particular genus of bacteria , proportion of bacteria ( including shifts and the proportion within the same phyla , class , order , or family ) , and percentage of bacteria with the same function should be taken into account in consideration of the influence of microbiota on the development of different types of diabetes .
the differences in microbiota between patients with different types of diabetes could be also a novel strategy in distinguishing subtypes of diabetes . in conclusion
, it appears that there are differences in gut microbiome composition between patients with hnf1a - mody , control individuals , and patients with type 2 diabetes .
 
 background . 
 type 2 diabetes mellitus ( t2 dm ) is determined by genetic and environmental factors . there have been many studies on the relationship between the composition of the gastrointestinal bacterial flora , t2 dm , and obesity . 
 there are no data , however , on the gut microbiome structure in monogenic forms of the disease including maturity onset diabetes of the young ( mody ) 
 . methods . 
 the aim of the investigation was to compare the qualitative parameters of the colonic flora in patients with hnf1a - mody and t2 dm and healthy individuals . 
 16s sequencing of bacterial dna isolated from the collected fecal samples using the miseq platform was performed . 
 results . 
 there were significant between - group differences in the bacterial profile . at the phylum level , the amount of proteobacteria was higher ( p = 0.0006 ) and the amount of bacteroidetes was lower ( p = 0.0005 ) in t2 dm group in comparison to the control group . in hnf1a - mody group , 
 the frequency of bacteroidetes was lower than in the control group ( p = 0.0143 ) . at the order level , turicibacterales was more abundant in hnf1a - mody group than in t2 dm group . 
 conclusions . 
 it appears that there are differences in the gut microbiome composition between patients with hnf1a - mody and type 2 diabetes 
 . further investigation on this matter should be conducted . 
1. Introduction
2. Materials and Methods
3. Results
4. Discussion