( ( ( ( ( ( ( (! ! 1!
|
|
- Marilynn Lambert
- 5 years ago
- Views:
Transcription
1 ( ( ( ( ( ( ( ( SocialBehaviorandGutMicrobiotainWild Red7belliedLemurs(Eulemur(rubriventer) (In SearchoftheRoleofImmunityintheEvolutionof Sociality( ProGradu AuraRaulo DepartmentofBiosciences DivisionofEcologyandEvolutionaryBiology UniversityofHelsinki May2015 1
2 Table&of&Contents& 1. INTRODUCTION GUTMICROBIOTAANDTHEHOLOBIOMEAPPROACH SOCIALBEHAVIORANDMICROBETRANSMISSIONINASOCIALGROUP 5 SOCIALBEHAVIORANDDIVERSITYOFMICROBIOTA 5 SOCIALBEHAVIORANDSIMILARITYOFMICROBIOTA 6 SOCIALBEHAVIORANDPRESENCEOFCERTAINBENEFICIALMICROBES STRESSANDMICROBIOTA AIMS STUDYSPECIESANDSTUDYSITE MATERIALSANDMETHODS COLLECTIONOFBEHAVIORALDATA 12 CONSTRUCTINGSOCIALITYINDICESANDINTERACTIONMATRICES FECALSAMPLECOLLECTIONANDANALYSISFORCORTISOLMETABOLITES FECALSAMPLECOLLECTIONANDGUTMICROBIALDNAANALYSIS STATISTICALANALYSES 16 DESCRIBINGMICROBIALDATA 16 CLUSTERINGOFMICROBIALDATA 17 ANALYSESUSINGDISSIMILARITYINDICES 18 ANALYSESUSINGDIVERSITYANDRICHNESS 18 WHICHBACTERIALTAXAAREBEHINDTHETREND? RESULTS ENVIRONMENTAL,DEMOGRAPHICANDGENETICFACTORSAFFECTINGMICROBIOTA SEASONALCHANGEINGUTMICROBIOTA GUTMICROBIALCOMPOSITIONANDINTERACTIONMATRICES GUTMICROBIALCOMPOSITIONANDSOCIALITYINDICES GUTMICROBIALCOMPOSITIONANDSTRESSHORMONES DISCUSSION ENVIRONMENTALANDDEMOGRAPHICVARIATIONINGUTMICROBIOTA GROUPASPECIFICMICROBIOTA MICROBIALTRANSMISSIONINASOCIALNETWORK GUTMICROBIOTA,HPAAAXISANDINDIVIDUALDIFFERENCESINSOCIALBEHAVIOR SEASONALCHANGEINGUTMICROBIOTA FUTUREDIRECTIONSANDAFRAMEWORKFORSOCIALTRANSMISSIONOFBENEFICIAL GUTMICROBIOTA CONCLUSIONS ACKNOWLEDGEMENTS 36 2
3 1. Introduction& 1.1. Gutµbiota&and&the&holobiome&approach&& Microbialcommunitieslivinginthegut,skinandglandsofvertebratesareofgrowing scientificinterest,duetotheiremergingroleasafunctionallinkbetweentheindividual anditssurroundingecosystem.recentdevelopmentsinmolecularmethodshave inspiredagreatamountofresearchonmutualisticmicrobiotainalltaxa.infields rangingfromendocrinologyandphysiologytobehavioralandevolutionaryecology,this new holobiomeapproach (Rosenbergetal.2007)focusesonthereciprocaland dynamicrelationshipsofmutualisticmicrobiotaandtheirhost. Mostofourcurrentknowledgeongutmicrobiotacomesfromhumanmedical researchfocusingonpathogens,whilethemutualisticmicrobiotaofwildvertebrate specieshasreceivedmuchlessattention.furthermore,whiletheeffectsofmicrobiota onhealthareincreasinglyacknowledged,currentresearchisforthemostpartlacking evolutionaryperspectiveonthistopic(butseestillingetal.2014).itisknown, however,thatvertebrategutmicrobiotahaveanimportantinfluenceonindividual immunity(reviewedbyhooperetal.2012)andisneededforprovidingenzymesfor digestion(mackie2000).forexampleruminantsareknownfortheirbacteria7rich rumenthatslowsdownthepassageoffoodmaterialtoenableefficientfermentationof cellulosebymutualisticmicrobes(mackie2000).newevidencealsosuggeststhat consistentwiththe"fermentationhypothesis"(alboneetal.1974,albone1990), microbiota,especiallythatlivinginmammalianscentglands,affectssocialrecognition bysynthesizingpheromones(sinetal.2012,theisetal.2012,2013,douglas&dobson 2013).Inadditiontothispotentialtoalterbehaviorindirectly,anorganism s microbiotaisknowntohavedirecteffectsonitsmoodandbehavior(bravoetal.2011, Ezenwaetal.2012).Thedirecteffectsofgutmicrobesonbehaviorlikelyhappen throughmicrobe7gut7brain(mgb)7axis,aphysiologicalpathwayinvolvingthevagus7 nerve(bravo2011,montiel7castroetal.2013). Mutualisticmicrobesplayanotableroleinananimal sphenotype,butthe determinantsofindividual smicrobiotaarescantlyknown.studiesofwildspecies, 3
4 laboratoryrodentsandhumansyieldinconsistentresultsregardingtherelative importanceofgeneticbackground,maternaltransmission,andsubsequent environmentalcontactsindetermininggutmicrobialcompositionsinceallthesefactors canhaveamajorimpactonthemicrobiota(pendersetal.2006,lanyonetal.2007, Khachatryanetal.2008,Turnbaughetal.2009a,b,Bensonetal.2010,Degnanetal. 2012,Ubedaetal.2012,Stevensonetal.2014,Leclaireetal.2014).Firstly,itisknown thatthecompositionofmicrobiotacanhaveageneticbasis.forexample,immunity7 relatedgenesareresponsibleforselectivelyenablingtheestablishmentofcertain microbesinthebody,andthusformthephysiologicalbasisformicrobiotainmice (Lanyonetal.2007,Zomeretal.2009,Bensonetal.2010).Ontheotherhand, relatednesswithinhumanfamiliesdoesnotpredictsimilarityofgutmicrobiota (Turnbaughetal.2009b).Secondly,themicrobesthatcolonizeananimalbodycome primarilyfromthemother.manygutmicrobescannotliveoutsidethehostbody,and mustrelyentirelyonmaternaltransmission(reviewedbyfunkhouser&bordenstein 2013).Thus,aswithmanyphenotypictraits,earlylifeisimportantindetermining microbiota.forexamplebirth(mandaretal.1996)andearlypostnatalexposureby nursing(jinetal.2011,fernandezetal.2013)orevenbyeatingsoilorfeces(troyer 1984b,Soave&Brand1991)havebeenconsideredthemostimportantperiodsof transmissiondeterminingindividualmicrobiota.inaddition,recentevidencesuggests thatpartofthemicrobiotacanbeacquiredpriortobirth(funkhouser&bordenstein 2013).Lastly,gutmicrobialcompositionisneverthelessknowntofluctuatewith changingdiet,hormonalprofiles,seasonandenvironment,hencehighlightingthe environmentalimpactonmicrobiota.whilesomeoftheseeffectsareduetoachanging dietandtheabioticenvironment,socialenvironmentislikelyanimportantfactoras well.thusamixtureofvertical(genetic,maternal)andhorizontal(social, environmental)transmissionmostlikelydetermineindividualmicrobiota.forexample, bothspottedhyenas(crocuta(crocuta)andchimpanzees(pan(troglodytes)havebeen foundtohavemoresimilarmicrobiotawithingroupsduetosharedenvironmentand socialcontact,butimmigrantindividualsstillretainedaslightlydifferentmicrobiota afteryearsofco7living(degnanetal.2012,theisetal.2012) 4
5 1.2. Social&behavior&andµbe&transmission&in&a&social&group& Microbialtransmissionhasbeenmainlystudiedinthecontextofpathogensand parasites.itisknownthatoneimportantfactoraffectingpathogenorparasite transmissionissocialcontactandthisisconsideredoneofthemaincostsofsociallife (Alexanderetal.1974,Mölleretal.1993,Altizeretal.2003).Forexamplevulnerability toinfectionsincreaseswithincreasinggroupsizeininsects(turnbulletal.2011). However,itislessoftensuggestedthatinadditiontoentailingimmunologicalchallenge, socialitymightactuallyalsoenhanceimmunitythroughtransmissionofbeneficial microbes(troyeretal.1984a,lombardo2008,archie&theis2011,gilbert2015).this isanimportantandunderstudiedpointofviewwhenbalancingtheprosandconsof sociallifestyletounderstandtheevolutionofsociality.tosketchtheroleofgut microbiotaintheevolutionandsocialdynamicsofsocialspecies,oneneedstopay attentiontothreeimportantaspectsofthemicrobiotainasocialcontext:diversity (affectsimmunity),similaritybetweenindividuals(affectsrecognitionandimmunity) andthepresenceofcertainbacteriawithknownbehavioralimpacts. Social'behavior'and'diversity'of'microbiota' Manysocialbehaviorsincludephysicalcontactandcanfunctionaspotentialpathways formicrobialtransmission.forexample,kulkarni&heeb(2007)experimentally infectedzebrafinches(taeniopygia(guttata)withbacillus(licheniformis(bacteria(and followedhowtheinfectionspreadintheirsocialgroup;preeningandsexualbehaviors transmittedthebacteriatoallgroupmembersinlessthanaday.similarly,frequent intimatekissingenhancedmutualtransmissionofmouthmicrobiotainhumans(kortet al.2014).behaviorcanaffectmicrobialtransmissionandrecentfindingssuggestthat microbiotacanaffecthostbehavioraswell(ezenwaetal.2012,montiel7castroetal. 2014).Forexamplegerm7freemicewerefoundtobemoreactiveandlessanxiousthan micewithnormalmicrobiota(heijtzaetal.2011).similarly,autism7relatedlackofpro7 socialityisassociatedwithalteredmicrobiota(mingetal2012,caoetal.2013) Manyauthorshavesuggestedthatpro7socialandaffiliativebehaviors,suchas grooming,licking,preening,kissingorsocio7sexualbehaviors,mighthavepartly evolvedtoservebeneficialmicrobialtransmission(troyer1984,lombardoetal.2008, Ezenwaetal.2012,Montiel7Castroetal.2013).Insocialgroups,exchangingmutualistic 5
6 microbesbetweenmultipleindividualsallowsmorediversemicrobiotatobepresent. Diversemicrobiotainturnhasbeenlongsuggestedtobearequisiteforaresilient immunity(hooperetal.2012,lozuponeetal.2012),inthesamewaybiodiversity makesecosystemsmoreresilienttochange(levineetal.1999,gunderson2000).for example,sociallytransmittedmicrobiotawasfoundtoprotectbumblebees(bombus( terrestris)fromacommonparasite(koch&schmid7hempel2011).whilesociality bringsimmunechallengesintheformofincreasedexposure(infectionfromothers)and evensusceptibility(immunosuppressivesocialstress)tosimilarlytransmittable parasites(alexander1974,mölleretal.1993,altizeretal.2003,turnbulletal.2011), transmissionofbeneficialbacteriamightcompensatethis.thisisanimportantissuefor thetheoriesoftheevolutionofsociality.becausemutualisticmicrobescanaffect behaviorandsocialbehaviorservesasatransmissionrouteforthesemicrobes,a positivefeedbackloopofincreasedsocialityandmicrobialtransmissionisapossible mechanismintheevolutionofsociality.thiscouldhelpexplaintheassociationbetween socialityandhealthfoundinhumansandotheranimals(houseetal.1988,cohenetal. 2003,Sapolskyetal.2004,Holt7Lunstadetal.2010,Archieetal.2014) Social'behavior'and'similarity'of'microbiota' Manygroup7livinganimalsusefriendlyandaffiliativebehaviorstoenhancesocial bondsorgroupcohesionandtoavoidorreconcileconflictsbetweengroup7members (Vasey1995,Aurelietal.2002,Wittigetal.2008,Hohmannetal.2009,Radford2011). Whenmicrobesaretransmittedthroughsocialbehaviorinasocialnetwork,individuals mostcloselysociallylinkedcanbeexpectedtosharethesamemicrobes.mutualistic microbiotainsocialsystemshavereceivedrelativelylittleattention,whereas transmissionpatternsofparasitesinsocialnetworksarewellstudied(dreweetal. 2010,Godfreyetal.2010,Griffin&Nunn2011,Fenneretal.2011,Bulletal.2012, MacIntoshetal2012,Rimbachetal.2015).Forexample,MacIntoshetal.(2012)found thatfemalejapanesemacaques(macaca(fuscata(yakui)withmorecentralpositionina socialnetwork,andhencemoregroomingpartners,hadmoreendoparasites.consistent withthis,malesocialrolepredictedgutmicrobialcompositioninmeerkats(suricata( suricatta,leclaireetal.2014).whensocialrelationshipsaresymmetric,transmission leadstomoresimilarmicrobiotabetweensocialpartners.forexamplehumanfamily memberslivingtogethersharemoresimilargutmicrobiotawitheachotherthanwith 6
7 peopleoutsidethefamily,regardlessofrelatedness(turnbaughetal.2009b,songetal. 2013).Actually,evendogssharemoresimilargutmicrobiotawiththeirhumanowners thanotherdogsoutsidethefamily(songetal.2013).suchgroup7specificmicrobiota hasbeenfoundinyellowbaboons(tungetal.2015),gorillas(gomezetal.2015), hyenas(theisetal.2012),meerkats(leclaireetal.2014)andchimpanzees(degnanet al.2012).similarmicrobialcompositionsinsocialgroupshavebeensuggestedtoplaya roleinchemicalcommunication,sincemutualisticbacteriaareresponsiblefor synthesizingmostoftheodorsandpheromonesusedinrecognitionbyavarietyof species.forexampleinspottedhyenas,groupmemberssharea scentmarksignature, asimilarscentmarkchemicalcompositioncreatedbyasimilarscentglandmicrobial community(theisetal.2012). Thefactthatindividualsdifferintheirsocialpersonality(Englishetal.2010,Koski 2011,Seyfarthetal.2012,Neumann2013,Carteretal.2014)andpossiblyunderlying physiology,speaksforafluctuatingbalancebetweentheimmunologicalbenefitsand challengesofsociallife.thisbalancedeterminesforexampletherelativevaluesof microbialdiversityandsimilarityunderselection.eventhoughgainingadiverse microbialcommunitythroughsocialcontactmightbebeneficial,contactwith individualswithtoodissimilarmicrobiotamightbeharmful,becauseindividuals becomeco7adaptedtotheirmicrobesandunfamiliarbacteriamightbepathogenicina non7adaptedindividual(fengetal.2011,stillingetal.2014).inlinewiththis, Barribeauetal.(2012)concludedthattadpolesmostlikelypreferthecompanyof individualswithsimilarmhc7genotypeaswellassimilarmicrobiota.furthermore,such avoidanceofopportunisticpathogensmightbethereasonbehindmatechoicedrivenby similarmicrobiota.forexample,drosophila(melanogastermaleswereobservedtomate morereadilywithfemaleswhohadhadthesamedietandsimilargutmicrobiota (Sharonetal.2010). Social'behavior'and'presence'of'certain'beneficial'microbes' Gutmicrobiotaconsistsofavarietyofbacteria,ofwhichsomearemorerelevantfor socialbehaviorandtransmissionthanothers.microbesthatvarylittleinabundance andaresharedmoreorlessbyallindividualsofagivenspeciesarethoughttoformthe CoreMicrobiota (Turnbaughetal.2007,butseealsoLozuponeetal.2012).Moreover, coremicrobiotacanbethoughttoconsistofthematernallytransmittedpartof 7
8 microbiota.inaddition,somemicrobesaremoreeasilysociallytransmittableand highlyvariableinabundance(heretermedas PeripheralMicrobiota ).Forexampleina recentstudy,tungetal.(2015)foundthataspecificsetof sociallystructuredbacteria weretransmittedthroughsocialinteraction,leadingmoresimilarmicrobiotabetween socialpartners.oneimportantgenusamongthesesociallystructuredbacteriawas Bifidobacterium,whichisassociatedwithmanyhealthbenefitsinhumans(Turronietal. 2008).Anotherbacteriallineagewithemergingimportanceinsocialbehavioris Lactobacillus.RecentstudieshaverevealedthatLactobacillusbacteriaareassociated withdecreasedanxietyinmice(bravoetal.2011),possiblyduetotheircapabilityto promoteoxytocinsecretion(poutahidisetal.2013) Stress&andµbiota& Socialbehaviorandmicrobiotaseemtobemoreintimatelylinkedthanpreviously thought.onefactorentwinedwiththisrelationshipisstressphysiologysincestress affectsbothsocialbehaviorandmicrobiota.stressisknowntoaffectphysiological homeostasisandtoresultinmoreorlesssocialbehaviordependingonthecontext. Specifically,theeffectsofsocialstresscandifferfromtheeffectsofnon7socialstress. Theevolutionofsociallifehascreatedawholenewnicheforthehormonalstress response,mediatedbythehypothalamic7pituitary7adrenalaxis(hpa).whereas environmentalstresscanreducesocialproximity(barbosa&dasilva7mota2009, Overdorff&Tecot2006,Tecot2013),socialstresscaninfactenhanceaffiliationaimed atreducingstressandavoidingconflicts(devriesetal.1996,carter1998,enghetal. 2006,Stöweetal.2008).Itisclearthatstresscanaffectmicrobialtransmissionthrough itseffectsonsocialbehavior,butglucocorticoidstresshormonesalsohavedirecteffects onmicrobiota.hormonalprofilesareassociatedwithmicrobiotaingeneral,whichis evidentinmicrobiotadifferencesbetweenthesexes(alexyetal.2003,saagetal.2011, Markleetal.2013,Leclaireetal.2014),butglucocorticoidstresshormoneshave especiallyfastandprofoundeffectsongutmicrobes.inmiceandprimates,stressis knowntoreducegutmicrobialdiversity(bailey&coe1999,baileyetal.2010),reduce therelativeamountoflactobacilli(bailey&coe1999,galleyetal.2014,galley&bailey 2014)andchangerelativeabundancesofresidentbacteriaingeneral(Parketal.2003). Inhumans,stressfultimeperiodswithhighercortisollevelsareassociatedwith changesingutmicrobiota,mostimportantlyreducedlactobacilli(lizkoetal.1987, 8
9 Knowlesetal.2008).However,stresseffectsonmicribiotaarenotunidirectional; throughthevagusnerveofthegut7brain7axis(bravoetal.2011,montiel7castroetal. 2013)microbiotaisknowntoaffectHPA7axisfunctionaswell.Forexample,lackof functionalmicrobiotasensitizesmicetostressorsandincreasedhpa7axisreactivityand anxiety7likebehavior(sudoetal.2004,crumeyrolle7ariasetal.2014). Figure 1 9
10 Figure1depictstherelationshipsbetweenmicrobiota,immunity,socialbehavior andstressphysiology.thisholobiomicmaplinksanindividualtoitssurrounding ecosystem,blurringtheboundariesofindividualsasfunctionalunitsofnature.some interactions,suchasthosebetweenimmunityandmicrobiotaorstressphysiologyand socialbehavior,areextensivelystudied,butotherinteractions,suchasthosebetween socialbehaviororstressphysiologyandmicrobiotaorstressphysiologyand microbiota,callformoreresearchbeforewecandrawanextensivemapofthese interactions Aims& Toshedlightontheinteractionsofgutmicrobiota,socialbehaviorandHPA7axis,I collecteddataaboutsocialbehavior,taxonomiccompositionofgutmicrobiotaandfecal cortisolmetabolitesfromawildpopulationofred7belliedlemurs(eulemur(rubriventer)( intheranomafananationalpark,easternmadagascar.iaskwhetherpatternsofsocial contact(group7membership,groupsize,socialnetwork)areassociatedwithpatternsof gutmicrobialcomposition(similarity)andwhetherindividualdifferencesinsociality arerelatedtospecificaspectsofgutmicrobiota(diversity,(presence(of(certain(microbes).( Furthermore,Iaddressthequestionofgeneticvs.environmentaldeterminantsof microbiotabycomparinggutmicrobialcompositionoffamilymemberswithdifferent relatedness.lastly,iexaminecorrelationsbetweenfecalcortisollevelsandgut microbialcompositiontotakethefirststeptowardsrelatinghpa7activityandgut microbialcompositioninthisspeciesandtoraisefurtherquestionsabouttheoriginof possiblecovariance. Red7belliedlemursarecooperativeandhighlysocialprimates.Individualsvaryin howmuchtheycontributetosocialactionsandwithwhom.forexampleindividual significantlydifferinhowmuchallomaternalcaretheyprovidetoinfants(tecot& Baden,unpublisheddata).Byusingsocialnetworkanalysisandbyconstructing individualsocialityindicesfromtheamountofsocialbehaviorsanindividualexhibits,i examinetheintra7groupdifferencesinsocialbehavior.myoverallhypothesisisthatthe amountanddirectionofclosesocialinteractionispositivelycorrelatedwithgut microbialdiversityandsimilarity,respectively.thispatternshouldbevisiblewhen comparingi)individualswithdifferentlevelsofsocialinteractionsinthesamegroup,ii) 10
11 individualsbetweengroups(nointeraction),andiii)individualsfromgroupsof differentsizes.inaddition,theassociationbetweengutmicrobiotaandsocialbehavior couldpartlybemediatedbyinteractionwiththestresshormonesofthehpa7axis.to testthishypothesis,severalspecificpredictionscanbemade: 1. Gutmicrobialcompositionshouldbemoresimilarbetweengroup7membersthan betweenindividualsfromdifferentgroups 2. Gutmicrobialdiversityshouldbegreaterinlargergroups 3. Moresocialindividualsshouldhavemorediversegutmicrobiotathanlesssocial individuals. 4. Pairswithshortdistanceinasocialnetwork(frequentmutualinteraction) shouldhavemoresimilargutmicrobialcompositionthanpairswithagreater distance(rarelyinteractwitheachother) 5. Gutmicrobialcompositionandfecalcortisolmetabolitesshouldcovary. Specifically,highfecalcortisollevelsshouldbeassociatedwithlowgutmicrobial diversity Inadditiontosociallife,Iinvestigatewhethergeneticrelatednessandearly environmentwithingrouparedeterminantsofgutmicrobiotainthisspecies.totest thishypothesisipredictthefollowing:ifgenesaremoreimportantthanthesocial environmentindetermininggutmicrobiota,adultsshouldhavemoresimilargut microbialcompositionwiththeiroffspringthanwitheachother Study&species&and&study&site& IcarriedoutthebehavioralobservationsandfecalsamplecollectionintheRanomafana NationalPark,EasternMadagascar,fromAugust2013toFebruary2014.Theareaisa rainforestwithtopographicvariation,montanecloudforestgrowingonridgesand lowlandrainforestinvalleys(wright1999).myspecificstudysite,vatoharanana,isa highaltituderegionlocated5kmsouthofthecentrevalbioresearchstation,bythe edgeoftheforest.thereareatleast12lemurspecieslivinginthearea. Withinthisstudyarea,thefocalstudysubjectswereeightfamilygroupsofred7bellied lemurs(eulemur(rubriventer).additionalfecalsampleswerecollectedfromother groupswhenencountered(total36individuals).red7belliedlemursarefrugivorous strepssirrhines,livinginsmall,relativelystable,monogamousfamilyunitsof276 11
12 individualswithonebreedingpairandtheiroffspringofdifferentages(overdorff 1996b,OverdorffandTecot2006,Tecot2008).Theyanegalitarianspeciesthatexhibit allomaternalcareandhaveaseasonalbreedingcycle(tecot2008,2010).somefathers andjuvenileshelptoraiseinfantsandallfamilymembersparticipateinsocial interactions,suchasgrooming,huddlingorplayingwitheachother.groupsofred7 belliedlemurshavefixedterritoriesthatvaryintheirfruitavailabilityduringtheyear (Overdorff1996a,Tecot2008).Mydatawerecollectedduringthetimeoftheyearwhen infantswerebornandfruitavailabilitywasgenerallylow(tecot2008),thoughwithin thisseasonthedietoftheanimalschangedmanytimesdependingonwhichparticular treeswerefruiting.furthermore,becauseofthegeneralfruitingpatternsofthefood plants,allstudyindividualshadarelativelysimilardietatanygiventime. 2. Materials&and&Methods& 2.1. Collection&of&Behavioral&Data&& Therewerefourpeoplecollectingbehavioraldatainourfieldteam.Twoofuswere localfieldtechnicianswithextensiveexpertiseinlemurbehavioralandhormonal research.wefollowedtheeightfamilygroupsofred7belliedlemursinthe Vatoharananaareaopportunistically:Wheneveragroupwasfound,itwasfollowedfor therestoftheday.ikepttrackofgroupsizeandcompositionandofanychanges therein.continuousbehavioraldatawerecollectedfrom27focalindividualsineight groups,byfollowingeachgrouponewholedayopportunisticallyeverytimetheywere found(approximatelysevenhoursperday,everyonortwoweeks).behavioraldata wereobtainedbyusingscansamplingmethodwithfive7minuteintervals;everyfive minutes,eachgroupmemberwasclassifiedtoonebehavioralcategoryofthered7 Bellied7LemurAllomaternalCareProjectethogram(Table1;Tecot&Baden, unpublishedfieldmanual).forsocialbehaviors(mutualgrooming,huddling,playing),i recordedwithwhomthebehaviorwasexhibited.becausefourdifferentpeople participatedinthebehavioralclassification,inter7observer7reliabilitywastested repeatedly(gwet2001).weaccumulated157150hoursofbehavioraldatafromeach 12
13 group(table2).onlythegroupswithmorethantwoindividualsandmorethan50 hoursofbehavioraldatawereusedtoinvestigatequestionsrelatedtosocialbehavior. Table 1. Ethogram, Tecot & Baden, unpublished field manual Mark Behavior Description M( Move Movinginsideasingletree(lessthan3meters) T( Travel Directedmovementmorethan3meters R( Rest Animalismotionless(eyesopenorclosed) HUD( Huddle Restinginphysicalcontactwithothers,tailswrapped aroundeachother PPLY( PartnerPlay Disorientednon7continuousplayfulbehaviorbetween morethanoneindividuals SPLY( Self7play Disorientednon7continuousplayfulbehaviornot directedtootherindividuals SGR( Self7groom Groomingone sownfur AGR( Allogroom Groominganotherindividual RGR( ReceivingGroom Receivinggroomfromanotherindividual MGR( MutualGroom Twoindividualsgroomingeachothersimultaneously FD( Feed Handling,sniffingandeatingfooditems AGG( Aggression Fastandloudchasingorcuffingotherindividual AGM( Scentmarking Ano7genitalscentmarkingbehavior ' Table 2. Amount of behavioral data per group and the criteria used to choose groups to study social behavior. All groups were used in other analyses described in the methods Group GroupSize Hoursofbehavioraldata Usedtostudybehavior X X X X
14 Constructing'sociality'indices'and'interaction'matrices' Individualsocialityindiceswereconstructedbasedonthebehavioralscansampledata. Theywerecalculatedastheproportionsoftimeeachsubjectspentinparticularsocial behaviors.specifically,iconstructedtwoindicesfortwodifferentsocialbehaviors, respectively: 1. Huddleindex:proportionoftimespenthuddlingwithothersoutofthetotal amountofresttime. 2. Socialgroomindex:proportionoftimespentsociallygrooming(mutualgroom, allogroom,receivinggroom)outofthetotalamountoftimeobserved. Toconstructinteractionmatrices,Iusedastandardmethodforsocialnetworkanalyses (Whitehead2008).Pairwisecomparisonsofmutualinteractionsyieldedsocialdistance valuesforeachpairofindividuals.distancevalueswerecalculatedforapairof individualsastheproportionoftimespentinmutualinteractionoutoftotaltimespent insocialinteractionofthetwoindividualsineachpair.interactionmatriceswere constructedfortwosocialbehaviors: 1. Huddlematrix:adistancematrixofhowmucheachpairspenthuddling together 2. Groommatrix:adistancematrixofhowmucheachpairspentgroomingeach other Socialityindicesandinteractionmatriceswereconstructedseparatelyfortwotime periods:1)beforeinfantswereborn(september7octoberin2013),and2)andafter infantswereborn(december2013untiljanuary2014). 2.2&& Fecal&Sample&Collection&and&Analysis&for&Cortisol&Metabolites& Icollectedfecalsamplesforhormonalassaysbywrappingtheminfoilanddryingthem byfire(tecot2001,2008,2013).sampleswerekeptdrybystoringtheminabagfullof silicapearls.hormonalsampleswerecollectedfromindividualsofthefocalgroups wheneverpossibletodetectchangesinhormonalstatus.sampleswerekeptdryand broughttothelaboratoryfortheevolutionaryendocrinologyofprimates(leep)in UniversityofArizona.Fecalsampleswereground,siftedand0.1gwasextractedfrom eachsampleintoasolutionofhalfethanolhalfwater.sampleswerethenmixed 14
15 (vortexed),centrifugedandthesupernatantwaspouredoff.onemlofeachsamplewas takenintonewtubesand4mlofethylacetatewasaddedtofreeconjugantsteroids. Samplesweremixedandcentrifugedagainandthetoplayerwasaspiratedintoculture tubes,evaporatedinwaterbathafterwhich1mlof30%methanolwasadded.samples wereconditionedwith1mlofmethanoland1mlofwater(distilled).1mlofsample wasaddedandwashedwith5%methanol.steroidswereelutedfromthecolumnwith 2mlofmethanol,whichwasthenevaporatedandreplacedwith1mlofethanol. Cortisolassaysweredoneusinganenzymeimmunoassay(EIA)protocolby CoralieMunro:Enzyme7labeledantigen(horseradishperoxidasecortisol)wasdiluted ineia7buffer(1:100000)andmixed.0.5mlofsampleswerefirstevaporatedand0.3 mloftheantigendilution(t:hrp7eia)wasaddedtoallsamplesandstandards.samples andstandardswerethenmixedandtransferredtoassaytubes.forcontrol,0.3mlpure EIA7bufferwasaddedtoNSB7tubesand0.3mlT:HRP7EIAwasaddedtoB07assaytubes. Platescoatedwithsteroidconjugatewerewashedand0.1mlofeachsampleand standardmixturewereplatedinduplicate,shakenandincubatedfor2hoursina humiditychambertoenablebinding.plateswerethenwashedtoseparateboundand freecortisol.1mlsubstratewasadded,theplatewasmixedforfiveminutesand absorbancewascalculatedwithabiotekplatereader.theslopeofseriallydilutedfecal extractswasparalleltothatofthecortisolstandard,f=0.12(11,12),p=0.730.fecal cortisolassaysweredeterminedtobeaccurate(96.68±1.79%standarderrorofthe mean) 2.3& Fecal&sample&collection&and&gutµbial&DNA&analysis&& IcollectedfecalsamplesformicrobialDNAanalysisinEppendorftubesfilledwith RNAlater.Minimumofonesamplewascollectedfromallfocalindividuals.Second samplesforseasonalcomparisonwerecollectedfrom15individualsinfivegroups.to increasesamplesize,icollectedadditionalsamplesopportunisticallywhenever encounteringindividualsfromnon7focalgroups.eventhoughididnothavematching behavioraldataforthesesamples,theywereusedintheanalysisofdemographic correlates,suchaswithage,sexorgroupsize.altogether110microbialfecalsamples werecollected. 15
16 Eppendorftubeswerekeptinafreezer(780C)for6monthsbeforesendingthem totheknightlaboratoryintheuniversityofcoloradotobeanalyzed.dna7extraction wasdoneusingmobiopowersoil kits.pcrwasrunwithprimerstargetedatv17v3 regionofthe16srrnagene.theresultingampliconsweresequencedusingillumina MiSeqplatform.IusedQIIMEsoftwaretomatchbarcodeswithsamples,toeliminate badqualitydnaandtoflipreversesequences(caporasoetal.2010).here,quality filteringwasdonewithdefaultsettingsofthe split_libraries_fastq.py command(q>20, 1%riskmarginalinbases).Sequenceswereclusteredandmatchedwithknown taxonomicdatabase(greengeens)usingopenreferenceotupicking(pick_openm reference_otus.py)whileenablingreversestrandmatching.thiswaythelistofthe knownotus(operationaltaxonomicunit,comparableto species )thatmatchedwith thedatabasecouldbeenrichedwiththedenovoclusteredunknownotus.thisis important,becauseasignificantproportionofgutmicrobiota,especiallythatoflemurs, belongstotheseunknowntaxa(katherineamato,personalcommunication).herefor example,openreferenceotupickingyielded1,131,454reads(sequences)whereas closedreferenceotupicking(withoutdenovoclustering)yieldedonly304,163reads. Theamountofsamplesthatweresuccessfullysequencedandhadaminimumof3,000 readswas92.thedatawasrarefiedto3,000readspersample Statistical&analyses& Describingmicrobialdata IperformedstatisticalanalysesusingR(Rversion3.1.2, PumpkinHelmet,RCore Team2014).Themicrobialdatawerefirststandardizedasrelativeabundances(the frequencyofcountsofeachoturelativetothetotalcount)andconvertedinto dissimilarityindices.theuseddissimilarityindicesdescribetherelativedifferencesof microbialspeciescommunitiesbetweenindividualsbasedonpresence absence information(sørensendissimilarityindex)orrelativeabundance(bray Curtis dissimilarityindex)(legendre&legendre2012).thebray Curtisdissimilarityindexis anindexusedtocalculatethedissimilarityofspeciescompositionbetweentwosites (here,betweensamples).itiscalculatedasfollows: " " = 1 ", 16
17 where " isthesumofthelessercountsofeachotupresentinbothsamples"#, and "# arethetotalnumberofotuspresentinbothsamples.thesørensen dissimilarityindexisasimilarindex,butunlikebray Curtis,itonlyappliesfor presence absencedata.itiscalculatedasfollows: " =, wherecistheabsolutenumberofotussharedbythetwosamples,andaandbarethe totalnumberofotusineachsamples,respectively.iusedthesetwodifferentindicesto getanoverviewofgutmicrobialdissimilaritybasedonboth(a)havingdifferentotus presentorabsent(sørensen),and(b)havingdifferentrelativeabundancesofotus (Bray Curtis).SørensenindiceswerealsousedtoreinforceBray7Curtismetricsof dissimilarity,becausedifferencesinthepresence absenceofcertainbacteriaarea morereliablemeasureofdissimilarityofmicrobialcommunitiesthandifferencesin relativeabundances.inadditiontodissimilarityindices,iusedshannondiversityindex andrichnesstodescribethegutmicrobialdiversity. Clusteringofmicrobialdata Iuseddissimilarityindicestoconstructdissimilaritymatricesofallindividuals,and thesewereusedtoclusterdatatorevealthedistributionofgutmicrobiotaacrossthe studypopulation.clusteringfromthedissimilaritymatriceswasdoneusingthek7 meansmethod,whichisanalgorithmthatassignsdatapointstoclustersbyminimizing thesumofsquareddistancesofpointstothecenterofeachcluster(legendre& Legendre2012).Eightclusterswerefoundusingthismethod,whichcoincideswiththe numberofstudygroups. Sincethemonthlychangesinthegutmicrobiotaweresmall(controlledfor individualandpregnancy,r 2 =0.02,p=0.012),Ipooledthesamplesfromtwomonths beforeandafterinfantswereborn(september OctoberandDecember January)to formtwocomparabletimeperiods(seasons).clusteringwasthendoneonasubsetof thedatacontainingonlyonesampleperindividual,alltakenduringthefirstseason (Sep Oct).Here,clusteringwasdonewithadefinednumberofeightclusters,derived fromtheobservedclusteringpatternofthewholedata.thiswasfurthervisualized usingprincipalcoordinateanalysis(pcoa),whichprovidesagoodoverviewofthedata 17
18 distributionbeforeanyfurtheranalyses.pcoaaimsatfindingtheaxisalongthe multidimensionaldistributionthatexhibitsthemostvariation(legendre&legendre 2012);thisaxisiscalledthefirstprincipalcoordinate. Analysesusingdissimilarityindices Totesttheeffectoftime(monthandseason)andpregnancyonthemicrobiota composition,iusedpermutationalmultivariateanalysisofvarianceusingdissimilarity indicesasdistancematrices.thisanalysiswasperformedwiththeadonisfunctionas implementedintheveganpackageinr(oksanenetal.2014).permutational multivariateanalysiswaschosenbecausetheelementsinadissimilaritymatrixarenot independent.thistestisbasedoncomparingtheobservedassociationbetween dissimilarityindicesandothervariableswiththosearisingwhentheindexvaluesare repeatedlyrandomlypermuted.thisallowsreliableestimationofstatistical significance.thesametestwasusedforthesmallersubsetofdatawithonesampleper individualcollectedinthefirstseasontotestfortheeffectsofage,sex,groupandgroup size.similarly,icomparedsamplesofasetoffemalesandtheirmatesbeforeandafter pregnancytotestforthepossibleeffectsofpregnancy,andusedprincipalcoordinate analysistovisualizethis. Analysesusingdiversityandrichness Inadditiontothese,Iusedgeneralizedlinearmodelstotestfortheeffectofseasonand socialityindicesongutmicrobialdiversityandrichness.thiswasdonebyfittingaline bylinearregressionthedata(minimizingthesumofsquaredresiduals)andcalculating whetherthisregressionexplainssignificantlymoreoftheobservedvariationthat wouldbeexpectedbychance,assumingadistributionofexponentialfamily(grafen& Hails2002).Furthermore,totestcorrelationsbetweeninteractionmatrices(converted fromsocialproximitytosocialdistancematrices)andmicrobialdissimilaritymatrices,i usedthemanteltestofmatrixcorrelation(legendre&legendre2012).themanteltest isapermutationaltestofcorrelation,wherestatisticalsignificanceistestedby repeatedlyrandomlypermutingtherowsandcolumnsofoneofthematrices,andthe correlationcoefficientsarecomparedwiththatobservedfortheoriginalmatrices.p7 valueiscalculatedasaproportionofpermutationsthatleadtohigherthanthe observedcorrelationcoefficient.fortestsinvolvingsocialityindicesandinteraction 18
19 matricesiuseddataonlyfromgroups,forwhichtherewassufficientbehavioraldatato constructtheseindicesandmatricesforbothseasons(beforeandafterinfantswere born).thissubsetofdataconsistedof26samplesencompassing13individualsfrom fourgroups. Whichbacterialtaxaarebehindthetrend? Tofindoutwhichbacterialtaxaareresponsibleforagiventrend,IusedDufrene LegendreindicatoranalysistocalculateindicatorvaluesforeachOTUinthegroupsof samplesundercomparison(dufrene&legendre1997).thedufrene Legendre indicatorvalueiscalculatedasfollows: "#$%& " = " " 100,where " isthenumberofotuiinclassj(suchasthefamilygroup)dividedbythe totalnumberofotui,and " isthenumberofsamples(individuals)inclassjthat containotuidividedbytotalamountofsamplesinclassj.thestatisticalsignificanceof theresultingindicatorvaluesarethentestedbyrandompermutation.here,p7valueof 0.01(except0.05whenanalyzingtheeffectofpregnancywithasmalldataset)wasused asathresholdforsignificance,andofeachclass,tenotuswiththehighestindicator valueswerechosen. 3. Results& 3.1. Environmental,&demographic&and&genetic&factors&affectingµbiota&& Thegutmicrobiotaofred7belliedlemurswasdominatedbythephylaBacterioidetes, ProteobacteriaandFirmicutes(Fig.5).Withinoneseason,gutmicrobiotaclustered stronglyaccordingtothefamilygroupandthegroupsformedthreebiggerclusters partiallyparallelingthelocationoftheirterritories(fig.2and3).thatis,thegroupsin thetwosmallerclustershaveneighboringoroverlappingterritories,butnotallgroups inthebiggestclusterseemtobesocloselylocated.individualdifferencesingut microbialcomposition(bybothdissimilarityindices)werebestexplainedbygroup identity(seetable3).inaddition,within7groupdifferencesweresignificantlyexplained 19
20 bypregnancy(seetable3).additionalanalysesrevealedthatmicrobiotaduringand afterpregnancychangedsignificantlyinadultfemales(r 2 =0.16,p=0.02withboth dissimilarityindices;seefig.4),butnotinacontrolgroupofthecorrespondingmales (p=0.5).thisdifferencewasnotinthegutmicrobialdiversity,butratherindifferences inthespecies(otus)ofthebacterialcommunity. Indicatoranalysisrevealedthatdifferencesbetweenthemicrobiotaofthesame femalesduringandafterpregnancyweredrivenbyphylum7levelchanges.specifically, pregnantfemales gutmicrobiotawascharacterizedbybacteriabelongingtophyla ProteobacteriaandVerrucomicrobia,whereaspost7pregnantfemales gutmicrobiota wascharacterizedbyabundanceofbacteriabelongingtophylumfirmicutes.moreover, thebestindicatorspeciesforgroup7specificdifferencesingutmicrobialcomposition werebacteriabelongingtotheorderclostridiales(inphylumfirmicutes).figure5 summarizesgroupdifferencesingutmicrobialcompositionatthephylum7level.group size,age,sexandage7sexinteractioncontrollingforpregnancyhadnosignificanteffect ongutmicrobialcomposition(seetable3). Table 3. Correlations of demographic factors (age, sex, group, group size) with the gut microbial composition (Sørensen dissimilarity index). All data are from the same season. Df SS MS F7value R 2 p( Pregnancy <0.01 Age Sex Group <0.01 GroupSize Age:Sex Residuals Total
21 Figure 2. PCA-ordination of the filtered data. Gut microbes cluster as family groups (numbered colored areas) and these groups form three bigger clusters, partially paralleling their territory location (See Fig. 3). Figure 3. Geographical territories of the study groups. Group colors are the same as in Fig
22 pca[,2] DURING AFTER pca[,1] Figure 4. Breeding females gut microbial compositions in relation to each other (lines) and during and after pregnancy (colors) visualized with principal component analysis Bacteroidetes Proteobacteria Verrucomicrobia Firmicutes Euryarchaeota Acidobacteria Planctomycetes Figure 5 Gut microbial composition of the eight study groups, allocated by phyla. White area covers unknown species 22
23 3.2. Seasonal&Change&in&Gut&Microbiota& Seasonalchangeinthegutmicrobialcompositionwasinvestigatedwitha datasetwherebreedingfemaleshadbeenremovedtoeliminatethestrongeffect ofpregnancy.therelativeabundancesofgutmicrobes(bray7curtisdissimilarity index)variedbetweenthetwoseasonssignificantly,thoughtheeffectwasquite small(table4).analysisusingpresence7absencedataformicrobialcomposition (Sørensendissimilarityindex)showednosignificanteffectoftheseason(Table 5).Gutmicrobialdiversity(Shannon)andrichnessdecreasedsignificantlyfrom beforetoafterinfantswereborn(p=0.02forrichness,p=0.04fordiversity,fig. 5).Asimilartrendindiversitywasnotdetectedinbreedingfemaleswhenthis wasanalyzedseparately.indicatoranalysisrevealedthatthetrendofdecreasing gutmicrobialdiversitywasassociatedwithincreasingabundanceofbacteria belongingtothephylumbacteroidetes,mostimportantlygenusprevotella. Table 4. The effect of season on gut microbial composition using relative abundances (Bray- Curtis dissimilarity index). Because data were unevenly distributed in time (Different groups had different amount of data in different months), Month, Group and their interaction were added to the analysis as control variables before Season. Df SS MS F7value R 2 p Group Month <0.01 Season <0.01 Group:Month Residuals Total Table 5. The effect of season on gut microbial composition using relative abundances (Sørensen dissimilarity index). Because data was unevenly distributed in time (Different groups had different amount of data in different months), Month, Group and their interaction were added to the analysis as control variables before Season. Df SS MS F7value& R 2 p Group Month <0.01 Season Group:Month Residuals Total
24 A (Before) B (After) A (Before) B (After) Figure 6 Gut microbial diversity (left) and richness (right) are higher before than after infants are born 3.3. Gut&Microbial&Composition&and&Interaction&Matrices& Therewasnosignificantcorrelationbetweengutmicrobialdissimilarity(by Bray7CurtisorSørensendissimilarityindices)andsocialconnectednesswithin groupwithanyoftheinteractionmatrices(huddlematrixandgroommatrix beforeandafterinfantswereborn) Gut&Microbial&Composition&and&Sociality&Indices& Individualsvariedintheamountofrestingtimetheyspenthuddlingtogether andtheamountoftimetheyspentgroomingothers(fig.7).groomindiceswere significantlynegativelycorrelatedwithgutmicrobialdiversitybetween individualsthroughtime(p=0.01).afteraddingtimetotheanalysis,negative correlationswereapparentinbothindices(table6),butonlybeforeinfants wereborn(fig.7.a,b).richnessfollowedasimilartrend(fig.7.c,d),thoughit wasnotsignificantlycorrelatedwiththehuddleindexatanytime.overall, huddleindexdecreased(p<0.01)andgroomindexdidnotchangefrombeforeto afterinfantswereborn. 24
25 Table 6. Results of generalized linear model testing for the effects of both sociality indices on gut microbial diversity, after the effect of time (Pre/Post infant) is taken into account. Estimate( Std.Error( tvalue( p (Intercept) <0.001 Pre/PostInfant GroomIndex Estimate( Std.Error( tvalue( p (Intercept) <0.001 Pre/PostInfant HuddleIndex Figure7Gutmicrobialdiversityandrichnessinrelationtosocialityindices.Beforeinfantsare born(inred),diversityisnegativelycorrelatedwithgroomindex(a)andhuddleindex(b),but correlationisnotapparentafterinfantsareborn(inblack).richnessfollowssimilartrends(c,d), butcorrelationwithhuddleindex(d)isnotsignificant 25
26 3.5. Gutµbial&composition&and&stress&hormones& TocontrolfortheeffectsofpregnancyonHPA7axisactivity,breedingfemales wereexcludedfromtheanalysisofcorrelationsbetweencortisolprofilesandgut microbialcomposition.fecalcortisollevelswerepositivelycorrelatedwithgut microbialcomposition(r 2 =0.14,p=0.05withSørensendissimilarityindex,R 2 = 0.14,p=0.04withBray7Curtisdissimilarityindex).Dissimilarityindices correlatedwitheachother(p<0.01),indicatingthattheyspeakforthesame thing.afteraddinggroupandsex7ageclassintotheanalysisascontrolvariables, theresultwasevenmoresignificant(table7).theeffectofcortisolongut microbiotawasnotrelatedtodiversityorrichness,butrathertodifferencesin thespeciesofthebacterialcommunity,specificallyspeciesbelongingtophylum Firmicutes. Table 7. Correlation of Cortisol levels (CORT) with gut microbial composition using relative abundances (Bray-Curtis dissimilarity index). Group and sex-age class are control variables. Df SS MS F7value R 2 p Group <0.01 Sex7AgeClass CORT Residuals Total Discussion& 4.1. Environmental&and&Demographic&Variation&in&Gut&Microbiota& Groupidentity,age,andseasonalchangewerethemainfactorsexplaining variationinthegutmicrobiotaofred7belliedlemurs,whichwasverysimilar betweenallindividualswithinonegroup.noeffectofrelatedness(differencesin microbiotabetweenparent7offspring,parent7parentandoffspring7offspring pairs)wasdetected.themainfactorexplainingwithin7groupdifferencesingut microbiotawaspregnancy.inparticularfemalesgutmicrobiotadifferedfrom theothergroupmember smicrobiotaandbreedingfemalesalsoexperienceda changeinthegutmicrobialcompositionatthetimeofgivingbirth.similarlyto 26
27 theseresults,themicrobiotaofpregnantwomen(colladoetal.2008,korenetal. 2012,Romeroetal.2014)aswellaswellbreedingfemaleEuropeanbadgers (Meles(meles)werefoundtodifferfromnon7breedingfemales'microbiota. Moreover,parallelingthetrendfoundbyKorenetal.(2012)inpregnantwomen, ourresultsshowthatpregnantfemalesgutisabundantinbacteriabelongingto thephylumproteobacteria,whilepost7pregnantfemales gutsareabundantin bacteriabelongingtothephylumfirmicutes.sincebirthandearlyenvironment areimportantinthedevelopmentofahealthymicrobiota,pregnancy7related modificationofgutmicrobiotamightbeanadaptivewayofreadyingthe microbialenvironmentfortheinfant(korenetal.2012).inadditionto pregnancy,noothersexoragedifferenceswerefound.thisisquitesurprising, giventhenumberofstudiesindicatingthatmicrobiotacandifferbetweensexes (Alexyetal.2003,Saagetal.2011,Markleetal.2013,Leclaireetal.2014)and tendstochangeduringdevelopmentuntilindividualsreachmaturity(palmeret al.2007,mariatetal.2009,sinetal.2012,leclaireetal.2014). Thedissimilarityofgutmicrobiotabetweenfamilygroupsmightbe partiallyduetodifferencesinthefoodsourcesandmicrohabitatsofdifferent territories.geographicalrangehasbeenpreviouslylinkedwithdifferencesingut microbiotainprimates(degnanetal.2012,moelleretal.2013,gomezetal. 2015)However,sincetheterritoriesofourstudygroupshaveoverlap, environmentaldifferencesseemaninadequateexplanationfortheobserved dissimilarityofmicrobiota.furthermore,red7belliedlemursarefrugivoresand allindividualswithinthesameareaexperiencemoreorlesssynchronized changesindietbecauseofthefruitingpatternsoftheirfoodplants.despitethis, theseasonalchangeingutmicrobiotawasnotnearlyasrobustasdifferences betweengroups,suggestingthatspatialvariationindietorenvironmental conditionscannotaloneexplainthedifferencesbetweengroups. However,group7specificmicrobiotasformedthreebiggerclustersthat mightbeexplainedbythelocationofgroupterritories(fig3).indeed,thegroups inthetwosmallerclusters(8+11and5+10)haveneighboringoroverlapping territories,butnotallgroupsinthebiggestcluster( )seemtobeso closelylocated.thisterritoryinformationisnotbasedongps7coordinates,buta trail7locatingsystemandmightbeslightlymisleading.morespecificinformation 27
28 ofterritoriesfromvariousdistancesofeachotherisneededbeforetheeffectof territorydistanceongutmicrobialcompositioncanbedescribed GroupKspecificµbiota& Group7specificgutmicrobiotareflectsthesociallifeofthisspecies,withhigh within7groupcohesionandextremelylowinteractionbetweenmembersof differentgroups.theseresultsaddtotheincreasingevidenceofgroup7specific microbiotainhighlysocialspecies(turnbaughetal.2009b,songetal.2013, Degnanetal.2012,Theisetal.2012,Leclaireetal.2014,Tungetal.2015,Gomez etal.2015).red7belliedlemurs,justlikecooperativelybreedingmeerkatsand hyenas,provideallomaternalcareandrelyonin7groupcooperationincaringfor theyoung.inthiscontext,sharedmicrobiotamightbecrucialforgroup7specific scentmarks,recognitionandbonding.moreimportantly,however,shared microbiotainasmallgroupcanbeseenasamechanismofimmunity synchronization;individualsbecomeaccustomedtotheirgutmicrobesand indeedthesamegutbacteriamightbemutualisticinanaccustomedindividual butpathogenicasanovelinfection(barribeauetal.2011,fengetal.2011, Stillingetal.2014).Sharingmicrobiotamakesitcertainthateverybodyis accustomedtothesamebacterialallies,andwillnotinfecttheirgroupmembers withpotentialpathogens.insupportofthishypothesis,cooperativebreeders andatleastonespecieswithsharedinfantcarearegenerallyknowntopossess synchronizedphysiologicalprocesseswithinthegroup.forexample,inmany cooperativelybreedingprimates,othergroupmemberscanreflectthehormonal statusofthepregnantfemale(storeyetal.,2000,nunesetal.,2001,ziegleretal., 2004,Tecot2008)possiblyreadyingthemtotheirtaskinprovidingallomaternal care(ziegler2000).furthermore,cooperativelybreedingmarmosets(callithrix( jacchus)werefoundtohavemoresimilarpersonalities,andthusmostlikely moresimilarhormonalprofiles(koolhaasetal.,1999,carereetal.2010, Koolhaasetal.,2010Sih2011)thanchimpanzeesrelyingondyadicratherthan group7widecooperation(koski&burkart2015).toexplorethisfurther, synchronizedhormonalpatternsincooperativegroupsneedstobestudied more. 28
29 4.3. Microbial&transmission&in&a&social&network& Ifoundnocorrelationbetweensocialcontactandsimilarityofmicrobiotawithin groups,butthisislikelybecausefamilygroupsofthisspeciesaresmall(276 individuals)andsocialconnectednesswithingroupisrelativelyhighbetweenall individuals.however,group7specificgutmicrobiotaspeakfortheimportanceof socialcontactindeterminingthemicrobiotainthisspecies.furthermore,even thebreedingpairsharedsimilargutmicrobiota,whichismostlikelyformedonly throughsocialorsexualcontact,sincetheyarenotrelatedanddonotshareearly environment(i.e.theyderivefromdifferentnatalgroups).thisisinteresting, becauseifsocialcontactcanfunctionasatransmissionrouteforgutmicrobes (Kulkarni&Heeb2007,Kortetal.2014,Tungetal.2015),thesemicrobescan carryinformationofsocialcontactaswell.inwildanimalresearch,thereisa strongbiasinthestudyofsocialbehaviortowardsdiurnalandhighlysocial species.gutmicrobiotacouldbeusedtodrawmapsofsocialnetworksalsoin asocial,nocturnalorotherwisecrypticspecies,whose(socialorsexual)behavior ishardtostudyinwild.parallellingthisidea,zohdyetal.(2012)managedto capturethepatternofsocialcontactsinwildnocturnalmouselemurs (Microcebus(rufus)bymarkingandrecapturingtheirparasites(lice).Whereas socialmicrobetransmissioninasocialspecieswouldmainlymeansexualor agonisticcontact,insocialspeciessimilarityofmicrobiotamightbeagoodproxy forthestrengthofasocialbond(montiel7castroetal.2013). Inadditiontosimilaritywithingroups,sociallytransmittedgutmicrobes mightenhancegutmicrobialdiversityasgroupsgetbigger.however,ifoundno correlationbetweengroupsizeandthegutmicrobialdiversity.thisislikelyto beduethesmallvariationingroupsize(276).inaddition,montiel7castroetal. (2013)pointedoutthatitisnotgroupsizeperse,butrathermodulargroup organizationthatisassociatedwiththebenefitsofsociallytransmittedmicrobes. Unlikeourpair7living,cooperativestudyspecies,biggerprimategroupstend havehighmodularity,leadingtopartialsocialisolationofdifferentsocial moduleswithinonegroup.thediversityofmicrobiotamighthavebeenone factordrivingtheevolutionoflargegroupsize,butselective(modular)bonding 29
30 patternsmayhaverisenwithaneedtobalancethetransmissionofbeneficial microbesandpathogens(montiel7castroetal.2013).indeed,socialmodularity isknowntoreduceparasitetransmission(griffin&nunn2011). Sincebondingbehaviorinprimatesservesasagoodwayformicrobial transmission,oneneedstodecidecarefullywithwhomtoaffiliate.affiliating withasmallsetofclosekinorfriendsisknowntoreducestress(wittigetal. 2008).Inadditiontothis,itmightalsoreducetheriskofinfection.Thus,asocial bondortheassociatedbondingbehaviorinhighlysocialspeciescanbe additionallyseenasamanifestationofimmunologicaltrust.inspecieswith allomaternalcare,likered7belliedlemursallindividualsaretightlybondedand seemtosharethemicrobiota,hencenomodularityisneeded Gut&Microbiota,&HPAKaxis&and&Individual&Differences&in&Social& Behavior& Contrarytomyprediction,bothsocialityindices(HuddleindexandGroom index)werenegativelycorrelatedwithgutmicrobialdiversitybetween individuals.thenegativecorrelationbetweensocialityandgutmicrobial diversitywassignificantandstrongbeforeinfantswereborn,notaftertheywere born.whereasgutmicrobialdiversitydecreasedfrombeforetoafterinfants wereborn,individualhuddleindicesdecreasedbutgroomindicesstayedthe same.thefactthatthenegativetrenddisappearsafterinfantsareborncouldbe explainedbydecreasingphysiologicaldissimilaritybetweenindividualsatthe timewheninfantsareborn.familymembersmostlikelyhavemoredifferent physiologicalstateswhilefemalesarepregnant.inaddition,individualscould havemoresimilarphysiologicalstatesafterinfantsarebornwhenallindividuals mightexperiencestressorothersynchronizedhormonalchangesrelatedto allomaternalcare(ziegler2000). Beforeinfantswereborn,individualsthatparticipatedmoreingrooming, groomedothersmoreandsleptinclosecontactwithothers,hadthelowestgut microbialdiversity.thisisoppositetowhatwaspredicted,assumingthat groomingfunctionsasawayofgutmicrobialtransmission.increasedmicrobial transmission,however,mightnotalwaysleadtohighermicrobialdiversity. 30
Supplementary Materials for
advances.sciencemag.org/cgi/content/full/2/1/e1500997/dc1 Supplementary Materials for Social behavior shapes the chimpanzee pan-microbiome Andrew H. Moeller, Steffen Foerster, Michael L. Wilson, Anne E.
More informationSoil microbial responses to climate perturbations Sarah E. Evans Kellogg Biological Station Michigan State University
Soil microbial responses to climate perturbations Sarah E. Evans Kellogg Biological Station Michigan State University The soil environment Harsh environment, many abiotic stressors Highly diverse, many
More informationRoot- associated Bacterial Communities in Flax (Linum usitatissimum) and Their Response to Arbuscular Mycorrhizal (AM) Inoculation
Root- associated Bacterial Communities in Flax (Linum usitatissimum) and Their Response to Arbuscular Mycorrhizal (AM) Inoculation M.A.K. Wijesinghe, J.J. Germida & F.L. Walley Department of Soil Sciences
More informationDETAILED RESULTS ON FUNGAL AND BACTERIAL COMMUNITIES COLONIZING THE COMPOST-INCUBATED POLYMER CARD SURFACE
DETAILED RESULTS ON FUNGAL AND BACTERIAL COMMUNITIES COLONIZING THE COMPOST-INCUBATED POLYMER CARD SURFACE PART 1: FUNGAL COMMUNITIES COLONIZING THE COMPOST-INCUBATED POLYMER CARD SURFACE Fungal OTUs richness,
More informationLecture 2: Diversity, Distances, adonis. Lecture 2: Diversity, Distances, adonis. Alpha- Diversity. Alpha diversity definition(s)
Lecture 2: Diversity, Distances, adonis Lecture 2: Diversity, Distances, adonis Diversity - alpha, beta (, gamma) Beta- Diversity in practice: Ecological Distances Unsupervised Learning: Clustering, etc
More informationMicrobiota: Its Evolution and Essence. Hsin-Jung Joyce Wu "Microbiota and man: the story about us
Microbiota: Its Evolution and Essence Overview q Define microbiota q Learn the tool q Ecological and evolutionary forces in shaping gut microbiota q Gut microbiota versus free-living microbe communities
More informationVariations in pelagic bacterial communities in the North Atlantic Ocean coincide with water bodies
The following supplement accompanies the article Variations in pelagic bacterial communities in the North Atlantic Ocean coincide with water bodies Richard L. Hahnke 1, Christina Probian 1, Bernhard M.
More informationKatrina Agno Lutap Iowa State University. Iowa State University Capstones, Theses and Dissertations. Graduate Theses and Dissertations
Graduate Theses and Dissertations Iowa State University Capstones, Theses and Dissertations 2015 Microbial community diversity in freshwater sponges from Lake Baikal/Protein size in three domains of life
More informationUsing Topological Data Analysis to find discrimination between microbial states in human microbiome data
Using Topological Data Analysis to find discrimination between microbial states in human microbiome data Mehrdad Yazdani 1,2, Larry Smarr 1,3 and Rob Knight 4 1 California Institute for Telecommunications
More informationMultivariate Statistics 101. Ordination (PCA, NMDS, CA) Cluster Analysis (UPGMA, Ward s) Canonical Correspondence Analysis
Multivariate Statistics 101 Ordination (PCA, NMDS, CA) Cluster Analysis (UPGMA, Ward s) Canonical Correspondence Analysis Multivariate Statistics 101 Copy of slides and exercises PAST software download
More informationH. Pieter J. van Veelen *, Joana Falcao Salles and B. Irene Tieleman
van Veelen et al. Microbiome (2017) 5:156 DOI 10.1186/s40168-017-0371-6 RESEARCH Open Access Multi-level comparisons of cloacal, skin, feather and nest-associated microbiota suggest considerable influence
More informationGenomics and Bioinformatics
Genomics and Bioinformatics Associate Professor Bharat Patel Genomes in terms of earth s history: Earth s environment & cellular evolution Genomes in terms of natural relationships: (Technology driven
More informationThe Thirteen Colonies
The Thirteen Colonies Objective: Students will create a 3 panel foldable that describes the geography, economy and social make up of the 3 regions of the English colonies For each region you must include
More informationSTAT Chapter 10: Analysis of Variance
STAT 515 -- Chapter 10: Analysis of Variance Designed Experiment A study in which the researcher controls the levels of one or more variables to determine their effect on the variable of interest (called
More informationBIO 682 Multivariate Statistics Spring 2008
BIO 682 Multivariate Statistics Spring 2008 Steve Shuster http://www4.nau.edu/shustercourses/bio682/index.htm Lecture 11 Properties of Community Data Gauch 1982, Causton 1988, Jongman 1995 a. Qualitative:
More informationAmplicon Sequencing. Dr. Orla O Sullivan SIRG Research Fellow Teagasc
Amplicon Sequencing Dr. Orla O Sullivan SIRG Research Fellow Teagasc What is Amplicon Sequencing? Sequencing of target genes (are regions of ) obtained by PCR using gene specific primers. Why do we do
More informationPCA Advanced Examples & Applications
PCA Advanced Examples & Applications Objectives: Showcase advanced PCA analysis: - Addressing the assumptions - Improving the signal / decreasing the noise Principal Components (PCA) Paper II Example:
More informationFIG S1: Rarefaction analysis of observed richness within Drosophila. All calculations were
Page 1 of 14 FIG S1: Rarefaction analysis of observed richness within Drosophila. All calculations were performed using mothur (2). OTUs were defined at the 3% divergence threshold using the average neighbor
More informationThe microbiome and its influence on plant growth and health. Friederike Trognitz
The microbiome and its influence on plant growth and health Friederike Trognitz Friederike.trognitz@ait.ac.at Fact sheet: 82% of variation in carbon cycling are explained by microbiological soil indicators
More informationOutline Classes of diversity measures. Species Divergence and the Measurement of Microbial Diversity. How do we describe and compare diversity?
Species Divergence and the Measurement of Microbial Diversity Cathy Lozupone University of Colorado, Boulder. Washington University, St Louis. Outline Classes of diversity measures α vs β diversity Quantitative
More informationCharacterizing and predicting cyanobacterial blooms in an 8-year
1 2 3 4 5 Characterizing and predicting cyanobacterial blooms in an 8-year amplicon sequencing time-course Authors Nicolas Tromas 1*, Nathalie Fortin 2, Larbi Bedrani 1, Yves Terrat 1, Pedro Cardoso 4,
More informationTrites and Larkin, 1996). The dashed line shows the division between the declining
Fig. 1. Locations of major geographic features cited in the text. The inserted graph shows estimated numbers of Steller sea lions (all ages) in Alaska from 1956 to 2000 (based on Trites and Larkin, 1996).
More informationInnovation and Regional Growth in the European Union
Riccardo Crescenzi Andres Rodriguez-Pose Innovation and Regional Growth in the European Union Springer Contents 1 Introduction 1 2 Theoretical Framework: A Spatial Perspective On Innovation and the Genesis
More informationDissimilarity and transformations. Pierre Legendre Département de sciences biologiques Université de Montréal
and transformations Pierre Legendre Département de sciences biologiques Université de Montréal http://www.numericalecology.com/ Pierre Legendre 2017 Definitions An association coefficient is a function
More informationA meta-analysis of changes in bacterial and archaeal communities with time
OPEN (2013), 1 14 & 2013 International Society for Microbial Ecology All rights reserved 1751-7362/13 www.nature.com/ismej ORIGINAL ARTICLE A meta-analysis of changes in bacterial and archaeal communities
More informationANOVA approach. Investigates interaction terms. Disadvantages: Requires careful sampling design with replication
ANOVA approach Advantages: Ideal for evaluating hypotheses Ideal to quantify effect size (e.g., differences between groups) Address multiple factors at once Investigates interaction terms Disadvantages:
More informationSeasonal, spatial, and maternal effects on gut microbiome in wild red squirrels
Ren et al. Microbiome (2017) 5:163 DOI 10.1186/s40168-017-0382-3 RESEARCH Seasonal, spatial, and maternal effects on gut microbiome in wild red squirrels Open Access Tiantian Ren 1, Stan Boutin 2, Murray
More information1Department of Demography and Organization Studies, University of Texas at San Antonio, One UTSA Circle, San Antonio, TX
Well, it depends on where you're born: A practical application of geographically weighted regression to the study of infant mortality in the U.S. P. Johnelle Sparks and Corey S. Sparks 1 Introduction Infant
More information19-1 Notes Bacteria. Named after the Greek word Little stick because many bacteria have a stick-like shape when viewed under a microscope
19-1 Notes Bacteria Bacteria are found virtually everywhere! Named after the Greek word Little stick because many bacteria have a stick-like shape when viewed under a microscope We know that bacteria are
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION ARTICLE NUMBER: 16161 DOI: 10.1038/NMICROBIOL.2016.161 A reference gene catalogue of the pig gut microbiome Liang Xiao 1, Jordi Estellé 2, Pia Kiilerich 3, Yuliaxis Ramayo-Caldas
More informationChemometrics. Matti Hotokka Physical chemistry Åbo Akademi University
Chemometrics Matti Hotokka Physical chemistry Åbo Akademi University Linear regression Experiment Consider spectrophotometry as an example Beer-Lamberts law: A = cå Experiment Make three known references
More informationOther resources. Greengenes (bacterial) Silva (bacteria, archaeal and eukarya)
General QIIME resources http://qiime.org/ Blog (news, updates): http://qiime.wordpress.com/ Support/forum: https://groups.google.com/forum/#!forum/qiimeforum Citing QIIME: Caporaso, J.G. et al., QIIME
More informationBivariate data analysis
Bivariate data analysis Categorical data - creating data set Upload the following data set to R Commander sex female male male male male female female male female female eye black black blue green green
More informationPractical Statistics
Practical Statistics Lecture 1 (Nov. 9): - Correlation - Hypothesis Testing Lecture 2 (Nov. 16): - Error Estimation - Bayesian Analysis - Rejecting Outliers Lecture 3 (Nov. 18) - Monte Carlo Modeling -
More information16S Metagenomics Report
Sample: BL-09 Report Date: 05/11/2017 18:40:20 Sample Information Sample ID: Sample Name: Run Folder: Taxonomy File: BL_09 BL_09 D:\Illumina\MiSeqAnalysis\170509_M01675_0059_000000000-AW8GB gg_13_5_species_32bp.dat
More informationBacteria Outline. 1. Overview. 2. Structural & Functional Features. 3. Taxonomy. 4. Communities
Bacteria Outline 1. Overview 2. Structural & Functional Features 3. Taxonomy 4. Communities Bacteria - Taxonomy PHYLUM CLASS ORDER FAMILY GENUS SPECIES SUB-SPECIES & STRAINS Bacteria - Phyla Firmicutes
More informationExploring Microbes in the Sea. Alma Parada Postdoctoral Scholar Stanford University
Exploring Microbes in the Sea Alma Parada Postdoctoral Scholar Stanford University Cruising the ocean to get us some microbes It s all about the Microbe! Microbes = microorganisms an organism that requires
More informationRevision: Chapter 1-6. Applied Multivariate Statistics Spring 2012
Revision: Chapter 1-6 Applied Multivariate Statistics Spring 2012 Overview Cov, Cor, Mahalanobis, MV normal distribution Visualization: Stars plot, mosaic plot with shading Outlier: chisq.plot Missing
More informationSection 5: Dummy Variables and Interactions
Section 5: Dummy Variables and Interactions Carlos M. Carvalho The University of Texas at Austin McCombs School of Business http://faculty.mccombs.utexas.edu/carlos.carvalho/teaching/ 1 Example: Detecting
More informationChad Burrus April 6, 2010
Chad Burrus April 6, 2010 1 Background What is UniFrac? Materials and Methods Results Discussion Questions 2 The vast majority of microbes cannot be cultured with current methods Only half (26) out of
More informationLab 11. Multilevel Models. Description of Data
Lab 11 Multilevel Models Henian Chen, M.D., Ph.D. Description of Data MULTILEVEL.TXT is clustered data for 386 women distributed across 40 groups. ID: 386 women, id from 1 to 386, individual level (level
More informationDr. Emily A. Fogarty, Coordinator History, Politics and Geography Dept School of Arts & Sciences
Geography Dr. Emily A. Fogarty, Coordinator History, Politics and Geography Dept. Emily.Fogarty@farmingdale.edu 631-420-2739 School of Arts & Sciences The History, Politics and Geography Department offers
More informationNiche Modeling. STAMPS - MBL Course Woods Hole, MA - August 9, 2016
Niche Modeling Katie Pollard & Josh Ladau Gladstone Institutes UCSF Division of Biostatistics, Institute for Human Genetics and Institute for Computational Health Science STAMPS - MBL Course Woods Hole,
More informationBacterial Community Variation in Human Body Habitats Across Space and Time
Bacterial Community Variation in Human Body Habitats Across Space and Time Elizabeth K. Costello, 1 Christian L. Lauber, 2 Micah Hamady, 3 Noah Fierer, 2,4 Jeffrey I. Gordon, 5 Rob Knight 1* 1 Department
More informationStructure, function and host control of rhizosphere microbiome
Structure, function and host control of rhizosphere microbiome Davide Bulgarelli PhD, Microbiome AgBioTech Europe London, September 21st, 2017 Outline The rhizosphere microbiome Barley as a model to study
More informationMultivariate Analysis of Ecological Data
Multivariate Analysis of Ecological Data MICHAEL GREENACRE Professor of Statistics at the Pompeu Fabra University in Barcelona, Spain RAUL PRIMICERIO Associate Professor of Ecology, Evolutionary Biology
More informationDETECTING BIOLOGICAL AND ENVIRONMENTAL CHANGES: DESIGN AND ANALYSIS OF MONITORING AND EXPERIMENTS (University of Bologna, 3-14 March 2008)
Dipartimento di Biologia Evoluzionistica Sperimentale Centro Interdipartimentale di Ricerca per le Scienze Ambientali in Ravenna INTERNATIONAL WINTER SCHOOL UNIVERSITY OF BOLOGNA DETECTING BIOLOGICAL AND
More informationDOWNLOAD OR READ : TRAIT BASED ECOLOGY FROM STRUCTURE TO FUNCTION PDF EBOOK EPUB MOBI
DOWNLOAD OR READ : TRAIT BASED ECOLOGY FROM STRUCTURE TO FUNCTION PDF EBOOK EPUB MOBI Page 1 Page 2 trait based ecology from structure to function trait based ecology from pdf trait based ecology from
More informationMultiple Linear Regression
Multiple Linear Regression ST 430/514 Recall: a regression model describes how a dependent variable (or response) Y is affected, on average, by one or more independent variables (or factors, or covariates).
More informationWhy Spatial Data Mining?
Intelligent Data Analysis for Spatial Data Mining Applications Wei Ding Knowledge Discovery Lab Department of Computer Science University of Massachusetts Boston Why Spatial Data Mining? Spatial Data mining
More informationThis gives us an upper and lower bound that capture our population mean.
Confidence Intervals Critical Values Practice Problems 1 Estimation 1.1 Confidence Intervals Definition 1.1 Margin of error. The margin of error of a distribution is the amount of error we predict when
More informationAnalysis of Variance: Part 1
Analysis of Variance: Part 1 Oneway ANOVA When there are more than two means Each time two means are compared the probability (Type I error) =α. When there are more than two means Each time two means are
More informationFecal Bacterial Communities as an Indicator of Trophic Interactions Among Anuran Larvae
The University of Southern Mississippi The Aquila Digital Community Master's Theses Fall 12-2016 Fecal Bacterial Communities as an Indicator of Trophic Interactions Among Anuran Larvae Steven Jacob Everman
More informationSTATISTICAL LEARNING OF INTEGRATIVE ANALYSIS. Meilei Jiang
STATISTICAL LEARNING OF INTEGRATIVE ANALYSIS Meilei Jiang A dissertation submitted to the faculty of the University of North Carolina at Chapel Hill in partial fulfillment of the requirements for the degree
More informationFigure 1. Locations of major geographic features cited in the text. The inserted graph
Figure 1. Locations of major geographic features cited in the text. The inserted graph shows estimated numbers of Steller sea lions (all ages) in Alaska from 1956 to 2000 (based on Trites and Larkin, 1996).
More informationAnnouncements. Lecture 1 - Data and Data Summaries. Data. Numerical Data. all variables. continuous discrete. Homework 1 - Out 1/15, due 1/22
Announcements Announcements Lecture 1 - Data and Data Summaries Statistics 102 Colin Rundel January 13, 2013 Homework 1 - Out 1/15, due 1/22 Lab 1 - Tomorrow RStudio accounts created this evening Try logging
More informationClimate controls prokaryotic community composition in desert soils of the southwestern United States
FEMS Microbiology Ecology, 93, 2017, fix116 doi: 10.1093/femsec/fix116 Advance Access Publication Date: 11 September 2017 Research Article RESEARCH ARTICLE Climate controls prokaryotic community composition
More informationEdwin A. Hernández-Delgado*
Long-term Coral Reef Ecological Change Monitoring Program of the Luis Peña Channel Marine Fishery Reserve, Culebra Island, Puerto Rico: I. Status of the coral reef epibenthic communities (1997-2002). Edwin
More informationCAP. Canonical Analysis of Principal coordinates. A computer program by Marti J. Anderson. Department of Statistics University of Auckland (2002)
CAP Canonical Analysis of Principal coordinates A computer program by Marti J. Anderson Department of Statistics University of Auckland (2002) 2 DISCLAIMER This FORTRAN program is provided without any
More informationCommunity Health Needs Assessment through Spatial Regression Modeling
Community Health Needs Assessment through Spatial Regression Modeling Glen D. Johnson, PhD CUNY School of Public Health glen.johnson@lehman.cuny.edu Objectives: Assess community needs with respect to particular
More information2/19/2018. Dataset: 85,122 islands 19,392 > 1km 2 17,883 with data
The group numbers are arbitrary. Remember that you can rotate dendrograms around any node and not change the meaning. So, the order of the clusters is not meaningful. Taking a subset of the data changes
More informationECOLOGY OF CALLIPHORID LARVAL MASSES AND POSTMORTEM COLONIZATION ESTIMATE VARIABILITY. Courtney R. Weatherbee
ECOLOGY OF CALLIPHORID LARVAL MASSES AND POSTMORTEM COLONIZATION ESTIMATE VARIABILITY By Courtney R. Weatherbee A THESIS Submitted to Michigan State University in partial fulfillment of the requirements
More informationIstituto di Microbiologia. Università Cattolica del Sacro Cuore, Roma. Gut Microbiota assessment and the Meta-HIT program.
Istituto di Microbiologia Università Cattolica del Sacro Cuore, Roma Gut Microbiota assessment and the Meta-HIT program Giovanni Delogu 1 Most of the bacteria species living in the gut cannot be cultivated
More informationMICROBIOLOGY ECOLOGY. Chronic N-amended soils exhibit an altered bacterial community structure in Harvard Forest, MA, USA.
RESEARCH ARTICLE Chronic N-amended soils exhibit an altered bacterial community structure in Harvard Forest, MA, USA Swathi A. Turlapati 1, Rakesh Minocha 2, Premsai S. Bhiravarasa 1, Louis S. Tisa 3,
More informationG E INTERACTION USING JMP: AN OVERVIEW
G E INTERACTION USING JMP: AN OVERVIEW Sukanta Dash I.A.S.R.I., Library Avenue, New Delhi-110012 sukanta@iasri.res.in 1. Introduction Genotype Environment interaction (G E) is a common phenomenon in agricultural
More informationSpatial Variation in Local Road Pedestrian and Bicycle Crashes
2015 Esri International User Conference July 20 24, 2015 San Diego, California Spatial Variation in Local Road Pedestrian and Bicycle Crashes Musinguzi, Abram, Graduate Research Assistant Chimba,Deo, PhD.,
More informationChapter 44. Table of Contents. Section 1 Development of Behavior. Section 2 Types of Animal Behavior. Animal Behavior
Animal Behavior Table of Contents Section 1 Development of Behavior Section 2 Types of Animal Behavior Section 1 Development of Behavior Objectives Identify four questions asked by biologists who study
More informationAn Analysis of College Algebra Exam Scores December 14, James D Jones Math Section 01
An Analysis of College Algebra Exam s December, 000 James D Jones Math - Section 0 An Analysis of College Algebra Exam s Introduction Students often complain about a test being too difficult. Are there
More informationGlobal Clinical Data Classification: A Discriminate Analysis
Global Clinical Data Classification: A Discriminate Analysis Amurthur Ramamurthy, Gordon Kapke and Jodi Yoder Covance Central Laboratories, Indianapolis, IN How different is Clinical Laboratory data sets
More informationINTRODUCTION TO MULTIVARIATE ANALYSIS OF ECOLOGICAL DATA
INTRODUCTION TO MULTIVARIATE ANALYSIS OF ECOLOGICAL DATA David Zelený & Ching-Feng Li INTRODUCTION TO MULTIVARIATE ANALYSIS Ecologial similarity similarity and distance indices Gradient analysis regression,
More informationTaxonomy and Clustering of SSU rrna Tags. Susan Huse Josephine Bay Paul Center August 5, 2013
Taxonomy and Clustering of SSU rrna Tags Susan Huse Josephine Bay Paul Center August 5, 2013 Primary Methods of Taxonomic Assignment Bayesian Kmer Matching RDP http://rdp.cme.msu.edu Wang, et al (2007)
More informationMultivariate Data Analysis a survey of data reduction and data association techniques: Principal Components Analysis
Multivariate Data Analysis a survey of data reduction and data association techniques: Principal Components Analysis For example Data reduction approaches Cluster analysis Principal components analysis
More informationIII Introduction to Populations III Introduction to Populations A. Definitions A population is (Krebs 2001:116) a group of organisms same species
III Introduction to s III Introduction to s A. Definitions B. characteristics, processes, and environment C. Uses of dynamics D. Limits of a A. Definitions What is a? A is (Krebs 2001:116) a group of organisms
More informationMicrobial diversity pa1erns in rela3on to hydrocarbon seepage
Microbial diversity pa1erns in rela3on to hydrocarbon seepage I. Microbial mat pushcores II. Sediment gravity cores The MC118 subsurface microbial community microbial players within and outside of subsurface
More informationPubH 7405: REGRESSION ANALYSIS. MLR: INFERENCES, Part I
PubH 7405: REGRESSION ANALYSIS MLR: INFERENCES, Part I TESTING HYPOTHESES Once we have fitted a multiple linear regression model and obtained estimates for the various parameters of interest, we want to
More informationClusters. Unsupervised Learning. Luc Anselin. Copyright 2017 by Luc Anselin, All Rights Reserved
Clusters Unsupervised Learning Luc Anselin http://spatial.uchicago.edu 1 curse of dimensionality principal components multidimensional scaling classical clustering methods 2 Curse of Dimensionality 3 Curse
More informationAnalysis of Bivariate Data
Analysis of Bivariate Data Data Two Quantitative variables GPA and GAES Interest rates and indices Tax and fund allocation Population size and prison population Bivariate data (x,y) Case corr® 2 Independent
More informationUser Guide. Affirmatively Furthering Fair Housing Data and Mapping Tool. U.S. Department of Housing and Urban Development
User Guide Affirmatively Furthering Fair Housing Data and Mapping Tool U.S. Department of Housing and Urban Development December, 2015 1 Table of Contents 1. Getting Started... 5 1.1 Software Version...
More informationAssessment Schedule 2016 Biology: Demonstrate understanding of the responses of plants and animals to their external environment (91603)
NCEA Level 3 Biology (91603) 2016 page 1 of 6 Assessment Schedule 2016 Biology: Demonstrate understanding of the responses of plants and animals to their external environment (91603) Evidence Statement
More informationGeographical knowledge and understanding scope and sequence: Foundation to Year 10
Geographical knowledge and understanding scope and sequence: Foundation to Year 10 Foundation Year 1 Year 2 Year 3 Year 4 Year 5 Year 6 Year level focus People live in places Places have distinctive features
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature25973 Power Simulations We performed extensive power simulations to demonstrate that the analyses carried out in our study are well powered. Our simulations indicate very high power for
More informationATINER's Conference Paper Series PLA
ATINER CONFERENCE PAPER SERIES No: LNG2014-1176 Athens Institute for Education and Research ATINER ATINER's Conference Paper Series PLA2015-1709 A Tool for Quantifying the Urban Gradient Agnes Lipovits
More informationSupplementary Information
Supplementary Information Table S1. Per-sample sequences, observed OTUs, richness estimates, diversity indices and coverage. Samples codes as follows: YED (Young leaves Endophytes), MED (Mature leaves
More informationCHAPTER 2: KEY ISSUE 1 Where Is the World s Population Distributed? p
CHAPTER 2: KEY ISSUE 1 Where Is the World s Population Distributed? p. 45-49 Always keep your vocabulary packet out whenever you take notes. As the term comes up in the text, add to your examples for the
More informationThe Scope and Growth of Spatial Analysis in the Social Sciences
context. 2 We applied these search terms to six online bibliographic indexes of social science Completed as part of the CSISS literature search initiative on November 18, 2003 The Scope and Growth of Spatial
More informationCatalonia is a small region, managed by an autonomous government (depending from Spain), and placed in NE. Spain, next to Mediterranean sea.
Characterization of the river basin according to the WFD The Catalan Basins, a case of small Mediterranean water district Planning Department c/ Provença, 204-208 08036-Barcelona Introduction to Catalonia
More informationExample. Multiple Regression. Review of ANOVA & Simple Regression /749 Experimental Design for Behavioral and Social Sciences
36-309/749 Experimental Design for Behavioral and Social Sciences Sep. 29, 2015 Lecture 5: Multiple Regression Review of ANOVA & Simple Regression Both Quantitative outcome Independent, Gaussian errors
More informationAppendix S1 Replacement, richness difference and nestedness indices
Appendix to: Legendre, P. (2014) Interpreting the replacement and richness difference components of beta diversity. Global Ecology and Biogeography, 23, 1324-1334. Appendix S1 Replacement, richness difference
More informationIntro to Linear Regression
Intro to Linear Regression Introduction to Regression Regression is a statistical procedure for modeling the relationship among variables to predict the value of a dependent variable from one or more predictor
More informationWorkshop 9.3a: Randomized block designs
-1- Workshop 93a: Randomized block designs Murray Logan November 23, 16 Table of contents 1 Randomized Block (RCB) designs 1 2 Worked Examples 12 1 Randomized Block (RCB) designs 11 RCB design Simple Randomized
More informationMicrobiome: 16S rrna Sequencing 3/30/2018
Microbiome: 16S rrna Sequencing 3/30/2018 Skills from Previous Lectures Central Dogma of Biology Lecture 3: Genetics and Genomics Lecture 4: Microarrays Lecture 12: ChIP-Seq Phylogenetics Lecture 13: Phylogenetics
More informationdiversity(datamatrix, index= shannon, base=exp(1))
Tutorial 11: Diversity, Indicator Species Analysis, Cluster Analysis Calculating Diversity Indices The vegan package contains the command diversity() for calculating Shannon and Simpson diversity indices.
More informationAdvanced Mantel Test
Advanced Mantel Test Objectives: Illustrate Flexibility of Simple Mantel Test Discuss the Need and Rationale for the Partial Mantel Test Illustrate the use of the Partial Mantel Test Summary Mantel Test
More informationMicrobes and mountains: metagenetics on Mount Fuji, Japan. Jonathan Adams, Biology Department, SNU, Korea
Microbes and mountains: metagenetics on Mount Fuji, Japan Jonathan Adams, Biology Department, SNU, Korea jonadams@snu.ac.kr Until about a decade ago, culturing could only yield 8,000 described species
More informationSupplementary Material
Supplementary Material The impact of logging and forest conversion to oil palm on soil bacterial communities in Borneo Larisa Lee-Cruz 1, David P. Edwards 2,3, Binu Tripathi 1, Jonathan M. Adams 1* 1 Department
More informationDepartment of Geography: Vivekananda College for Women. Barisha, Kolkata-8. Syllabus of Post graduate Course in Geography
India: Regional Problems and Resource management Module 11 (Full Marks 50) Unit I: Region and Regionalisation 1.1 Various bases of regionalisation of India; problems of identification and delineation.
More informationunadjusted model for baseline cholesterol 22:31 Monday, April 19,
unadjusted model for baseline cholesterol 22:31 Monday, April 19, 2004 1 Class Level Information Class Levels Values TRETGRP 3 3 4 5 SEX 2 0 1 Number of observations 916 unadjusted model for baseline cholesterol
More informationSimple Linear Regression
Simple Linear Regression ST 430/514 Recall: A regression model describes how a dependent variable (or response) Y is affected, on average, by one or more independent variables (or factors, or covariates)
More informationUnit 7: Multiple linear regression 1. Introduction to multiple linear regression
Announcements Unit 7: Multiple linear regression 1. Introduction to multiple linear regression Sta 101 - Fall 2017 Duke University, Department of Statistical Science Work on your project! Due date- Sunday
More informationMission Geography and Missouri Show-Me Standards Connecting Mission Geography to State Standards
Module 1: Volcanoes local hazard, global issue Mission Geography and Missouri Show-Me Standards Connecting Mission Geography to State Standards Grades 5-8 Inv Geography for Life State Standard(s) Connection
More information