<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.1d3 20150301//EN" "http://jats.nlm.nih.gov/publishing/1.1d3/JATS-journalpublishing1.dtd">
<article article-type="research-article" dtd-version="1.1d3" xml:lang="en" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="nlm-ta">PLoS ONE</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plosone</journal-id>
<journal-title-group>
<journal-title>PLOS ONE</journal-title>
</journal-title-group>
<issn pub-type="epub">1932-6203</issn>
<publisher>
<publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, CA USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.1371/journal.pone.0164823</article-id>
<article-id pub-id-type="publisher-id">PONE-D-16-25972</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Research Article</subject>
</subj-group>
<subj-group subj-group-type="Discipline-v3"><subject>Biology and life sciences</subject><subj-group><subject>Organisms</subject><subj-group><subject>Animals</subject><subj-group><subject>Vertebrates</subject><subj-group><subject>Amniotes</subject><subj-group><subject>Birds</subject><subj-group><subject>Ostriches</subject></subj-group></subj-group></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and life sciences</subject><subj-group><subject>Molecular biology</subject><subj-group><subject>Molecular biology techniques</subject><subj-group><subject>Artificial gene amplification and extension</subject><subj-group><subject>Polymerase chain reaction</subject></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Research and analysis methods</subject><subj-group><subject>Molecular biology techniques</subject><subj-group><subject>Artificial gene amplification and extension</subject><subj-group><subject>Polymerase chain reaction</subject></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and life sciences</subject><subj-group><subject>Genetics</subject><subj-group><subject>DNA</subject><subj-group><subject>Ancient DNA</subject></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and life sciences</subject><subj-group><subject>Biochemistry</subject><subj-group><subject>Nucleic acids</subject><subj-group><subject>DNA</subject><subj-group><subject>Ancient DNA</subject></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and life sciences</subject><subj-group><subject>Paleontology</subject><subj-group><subject>Paleogenetics</subject><subj-group><subject>Ancient DNA</subject></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Earth sciences</subject><subj-group><subject>Paleontology</subject><subj-group><subject>Paleogenetics</subject><subj-group><subject>Ancient DNA</subject></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and life sciences</subject><subj-group><subject>Paleontology</subject><subj-group><subject>Fossils</subject></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Earth sciences</subject><subj-group><subject>Paleontology</subject><subj-group><subject>Fossils</subject></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and life sciences</subject><subj-group><subject>Paleontology</subject><subj-group><subject>Paleogenetics</subject></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Earth sciences</subject><subj-group><subject>Paleontology</subject><subj-group><subject>Paleogenetics</subject></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and life sciences</subject><subj-group><subject>Organisms</subject><subj-group><subject>Animals</subject><subj-group><subject>Vertebrates</subject><subj-group><subject>Amniotes</subject><subj-group><subject>Birds</subject></subj-group></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Research and analysis methods</subject><subj-group><subject>Database and informatics methods</subject><subj-group><subject>Bioinformatics</subject><subj-group><subject>Sequence analysis</subject><subj-group><subject>DNA sequence analysis</subject></subj-group></subj-group></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Biology and life sciences</subject><subj-group><subject>Paleontology</subject><subj-group><subject>Paleobiology</subject></subj-group></subj-group></subj-group><subj-group subj-group-type="Discipline-v3"><subject>Earth sciences</subject><subj-group><subject>Paleontology</subject><subj-group><subject>Paleobiology</subject></subj-group></subj-group></subj-group></article-categories>
<title-group>
<article-title>Ancient DNA Reveals Late Pleistocene Existence of Ostriches in Indian Sub-Continent</article-title>
<alt-title alt-title-type="running-head">Ancient DNA Reveals Ostriches' Existence in India</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes" xlink:type="simple">
<contrib-id authenticated="true" contrib-id-type="orcid">http://orcid.org/0000-0003-0439-6797</contrib-id>
<name name-style="western">
<surname>Jain</surname>
<given-names>Sonal</given-names>
</name>
<xref ref-type="aff" rid="aff001"><sup>1</sup></xref>
<xref ref-type="corresp" rid="cor001">*</xref>
</contrib>
<contrib contrib-type="author" xlink:type="simple">
<name name-style="western">
<surname>Rai</surname>
<given-names>Niraj</given-names>
</name>
<xref ref-type="aff" rid="aff002"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author" xlink:type="simple">
<name name-style="western">
<surname>Kumar</surname>
<given-names>Giriraj</given-names>
</name>
<xref ref-type="aff" rid="aff003"><sup>3</sup></xref>
</contrib>
<contrib contrib-type="author" xlink:type="simple">
<name name-style="western">
<surname>Pruthi</surname>
<given-names>Parul Aggarwal</given-names>
</name>
<xref ref-type="aff" rid="aff001"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author" xlink:type="simple">
<name name-style="western">
<surname>Thangaraj</surname>
<given-names>Kumarasamy</given-names>
</name>
<xref ref-type="aff" rid="aff002"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author" xlink:type="simple">
<name name-style="western">
<surname>Bajpai</surname>
<given-names>Sunil</given-names>
</name>
<xref ref-type="aff" rid="aff004"><sup>4</sup></xref>
<xref ref-type="fn" rid="currentaff001"><sup>¤</sup></xref>
</contrib>
<contrib contrib-type="author" xlink:type="simple">
<name name-style="western">
<surname>Pruthi</surname>
<given-names>Vikas</given-names>
</name>
<xref ref-type="aff" rid="aff001"><sup>1</sup></xref>
</contrib>
</contrib-group>
<aff id="aff001"><label>1</label> <addr-line>Department of Biotechnology, Indian Institute of Technology, Roorkee, Uttarakhand, India</addr-line></aff>
<aff id="aff002"><label>2</label> <addr-line>CSIR-Centre for Cellular and Molecular Biology, Hyderabad, India</addr-line></aff>
<aff id="aff003"><label>3</label> <addr-line>Dayalbagh Educational Institute, Dayalbagh, Agra, UP, India</addr-line></aff>
<aff id="aff004"><label>4</label> <addr-line>Department of Earth Sciences, Indian Institute of Technology, Roorkee, Uttarakhand, India</addr-line></aff>
<contrib-group>
<contrib contrib-type="editor" xlink:type="simple">
<name name-style="western">
<surname>Calafell</surname>
<given-names>Francesc</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"/>
</contrib>
</contrib-group>
<aff id="edit1"><addr-line>Universitat Pompeu Fabra, SPAIN</addr-line></aff>
<author-notes>
<fn fn-type="conflict" id="coi001">
<p>The authors have declared that no competing interests exist.</p>
</fn>
<fn fn-type="con">
<p><list list-type="simple">
<list-item><p><bold>Conceptualization:</bold> SJ SB VP.</p></list-item>
<list-item><p><bold>Data curation:</bold> SJ NR PAP.</p></list-item>
<list-item><p><bold>Formal analysis:</bold> SJ NR PAP VP.</p></list-item>
<list-item><p><bold>Funding acquisition:</bold> SJ NR KT.</p></list-item>
<list-item><p><bold>Investigation:</bold> SJ NR.</p></list-item>
<list-item><p><bold>Methodology:</bold> SJ NR KT VP.</p></list-item>
<list-item><p><bold>Project administration:</bold> KT SB VP.</p></list-item>
<list-item><p><bold>Resources:</bold> SJ NR KT SB GK VP.</p></list-item>
<list-item><p><bold>Software:</bold> SJ NR PAP.</p></list-item>
<list-item><p><bold>Supervision:</bold> KT PAP SB VP.</p></list-item>
<list-item><p><bold>Validation:</bold> SJ NR PAP SB KT VP.</p></list-item>
<list-item><p><bold>Visualization:</bold> SJ NR PAP.</p></list-item>
<list-item><p><bold>Writing – original draft:</bold> SJ NR PAP.</p></list-item>
<list-item><p><bold>Writing – review &amp; editing:</bold> SJ NR PAP KT SB VP.</p></list-item></list></p>
</fn>
<fn fn-type="current-aff" id="currentaff001">
<label>¤</label><p>Current address: Birbal Sahni Institute of Palaeosciences, Lucknow, U.P., India</p>
</fn>
<corresp id="cor001">* E-mail: <email xlink:type="simple">jainsonal20@gmail.com</email></corresp>
</author-notes>
<pub-date pub-type="epub">
<day>8</day>
<month>3</month>
<year>2017</year>
</pub-date>
<pub-date pub-type="collection">
<year>2017</year>
</pub-date>
<volume>12</volume>
<issue>3</issue>
<elocation-id>e0164823</elocation-id>
<history>
<date date-type="received">
<day>29</day>
<month>6</month>
<year>2016</year>
</date>
<date date-type="accepted">
<day>1</day>
<month>10</month>
<year>2016</year>
</date>
</history>
<permissions>
<copyright-year>2017</copyright-year>
<copyright-holder>Jain et al</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/" xlink:type="simple">
<license-p>This is an open access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/" xlink:type="simple">Creative Commons Attribution License</ext-link>, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p>
</license>
</permissions>
<self-uri content-type="pdf" xlink:href="info:doi/10.1371/journal.pone.0164823"/>
<abstract>
<p>Ancient DNA (aDNA) analysis of extinct ratite species is of considerable interest as it provides important insights into their origin, evolution, paleogeographical distribution and vicariant speciation in congruence with continental drift theory. In this study, DNA hotspots were detected in fossilized eggshell fragments of ratites (dated ≥25000 years B.P. by radiocarbon dating) using confocal laser scanning microscopy (CLSM). DNA was isolated from five eggshell fragments and a 43 base pair (bp) sequence of a 16S rRNA mitochondrial-conserved region was successfully amplified and sequenced from one of the samples. Phylogenetic analysis of the DNA sequence revealed a 92% identity of the fossil eggshells to <italic>Struthio camelus</italic> and their position basal to other palaeognaths, consistent with the vicariant speciation model. Our study provides the first molecular evidence for the presence of ostriches in India, complementing the continental drift theory of biogeographical movement of ostriches in India, and opening up a new window into the evolutionary history of ratites.</p>
</abstract>
<funding-group>
<award-group id="award001">
<funding-source>
<institution-wrap>
<institution-id institution-id-type="funder-id">http://dx.doi.org/10.13039/501100001412</institution-id>
<institution>Council of Scientific and Industrial Research</institution>
</institution-wrap>
</funding-source>
<principal-award-recipient>
<name name-style="western">
<surname>Rai</surname>
<given-names>Niraj</given-names>
</name>
</principal-award-recipient>
</award-group>
<award-group id="award002">
<funding-source>
<institution-wrap>
<institution-id institution-id-type="funder-id">http://dx.doi.org/10.13039/501100001501</institution-id>
<institution>University Grants Commission</institution>
</institution-wrap>
</funding-source>
<principal-award-recipient>
<contrib-id authenticated="true" contrib-id-type="orcid">http://orcid.org/0000-0003-0439-6797</contrib-id>
<name name-style="western">
<surname>Jain</surname>
<given-names>Sonal</given-names>
</name>
</principal-award-recipient>
</award-group>
<funding-statement>For the financial support, we would like to acknowledge the Council for Scientific and Industrial Research (CSIR) and University Grants Commission (UGC), New Delhi.</funding-statement>
</funding-group>
<counts>
<fig-count count="5"/>
<table-count count="3"/>
<page-count count="14"/>
</counts>
<custom-meta-group>
<custom-meta id="data-availability">
<meta-name>Data Availability</meta-name>
<meta-value>All relevant data are within the paper and its Supporting Information file.</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body>
<sec id="sec001" sec-type="intro">
<title>Introduction</title>
<p>The origin and evolution of ratites, the flightless birds, is widely attributed to the continental drifting of Gondwanaland [<xref ref-type="bibr" rid="pone.0164823.ref001">1</xref>]. One hundred and fifty million years ago, Gondwanaland was a consortium of South America, Arabia, Africa, Australia, Antarctica, India and Madagascar [<xref ref-type="bibr" rid="pone.0164823.ref002">2</xref>]. Initial breakup of this supercontinent during the Early Cretaceous, ~130 to 100 million years ago (Ma), separated Africa and Indo-Madagascar (<xref ref-type="fig" rid="pone.0164823.g001">Fig 1</xref>). Later, the remaining Gondwana landmass further split into South America and Australia/Antarctica [<xref ref-type="bibr" rid="pone.0164823.ref001">1</xref>]. During the Late Cretaceous (~80 Ma), Australia/Antarctica and Indo-Madagascar may have been connected through Kerguelen Plateau until India drifted further northwards following its break up from Madagascar [<xref ref-type="bibr" rid="pone.0164823.ref003">3</xref>]. This continental drift theory is complemented by the vicariant speciation that led to the evolution of ratites [<xref ref-type="bibr" rid="pone.0164823.ref001">1</xref>]. In the fractured Gondwana, rhea and emu migrated to South America and Australia, respectively, while elephant birds and ostrich existed together on Indo-Madagascar [<xref ref-type="bibr" rid="pone.0164823.ref002">2</xref>]. Continued northerly movement of India split the Indo-Madagascar landmass at ~88 Ma, with elephant birds remaining in Madagascar and ostriches in India. This biogeographical dispersion eventually led to hopping of ostriches in Africa through Eurasia via a land route around ~20 Ma [<xref ref-type="bibr" rid="pone.0164823.ref003">3</xref>]. The above studies impliy that elephant birds and ostriches should be the oldest ratite lineages. However, recent investigations by Mitchell et al, 2014 suggested that kiwis and elephant birds form the closest group, in contradiction with the continental drift theory [<xref ref-type="bibr" rid="pone.0164823.ref001">1</xref>].</p>
<fig id="pone.0164823.g001" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0164823.g001</object-id>
<label>Fig 1</label>
<caption>
<title>Continental drifting and vicariant speciation of ratites.</title>
<p>Breakup of supercontinent Gondwanaland during Early Cretaceous and biogeographical movement of ratites to different continents.</p>
</caption>
<graphic mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0164823.g001" xlink:type="simple"/>
</fig>
<p>Known occurrences of fossil eggshells of ratites in the Indian subcontinent range in age from the Middle Miocene to Late Pleistocene epochs [<xref ref-type="bibr" rid="pone.0164823.ref004">4</xref>–<xref ref-type="bibr" rid="pone.0164823.ref005">5</xref>]. Aigner reported ratite eggshells in 1935 from the Late Miocene Siwalik rocks of Dhokpathan near Hasnot in Punjab, Pakistan [<xref ref-type="bibr" rid="pone.0164823.ref006">6</xref>]. Evidence from the Siwalik formations revealed the presence of Struthionidae (ostrich family) from the Middle Miocene to Middle Pleistocene. No later records of the family are known from the Himalayan foreland but it is known to have existed in Central India [<xref ref-type="bibr" rid="pone.0164823.ref007">7</xref>]. Skeletal remains of the Siwalik ratites were identified as <italic>Struthio asiaticus</italic> (1871), whose eggshells differ in size from those of <italic>Struthio camelus</italic> [<xref ref-type="bibr" rid="pone.0164823.ref008">8</xref>–<xref ref-type="bibr" rid="pone.0164823.ref009">9</xref>]. Recently, AMS dating of eggshells confirmed the earliest evidence of ostriches on the Indian subcontinent and their expansion into Central and Western India, prior to 60 thousand years ago (ka) [<xref ref-type="bibr" rid="pone.0164823.ref005">5</xref>]. The eggshells can be identified on the basis of morphological characters such as eggshell thickness, pore complex shape, density and diameter [<xref ref-type="bibr" rid="pone.0164823.ref010">10</xref>–<xref ref-type="bibr" rid="pone.0164823.ref013">13</xref>]. Eggshells with struthinoid type pore patterns were found in younger Siwalik deposits ranging in age from 2.24 Ma to 0.5 Ma [<xref ref-type="bibr" rid="pone.0164823.ref014">14</xref>]. However, from Indian sites (Rajasthan, Uttar Pradesh, M.P., Maharashtra), fossil eggshells dating between &gt;60000–18000 years ago have been identified as <italic>Struthio</italic> sp. [<xref ref-type="bibr" rid="pone.0164823.ref006">6</xref>,<xref ref-type="bibr" rid="pone.0164823.ref007">7</xref>,<xref ref-type="bibr" rid="pone.0164823.ref015">15</xref>]. This taxonomic identification is based on the pattern of decreasing pore diameter (1 mm), shell thickness (2.5 cm) and increasing pore density through time in comparison to other ostrich eggshells classified as <italic>Struthio karingarbensis</italic>, <italic>Struthio daberasensis</italic>, <italic>Struthio kakesiensis</italic> and <italic>S</italic>. <italic>camelus</italic> from Late Miocene and Pliocene of Namibia and East Africa [<xref ref-type="bibr" rid="pone.0164823.ref005">5</xref>]. Ostrich eggshells with a mean shell thickness of 2.54mm, identified as <italic>S</italic>. <italic>asiaticus</italic>, were found in Ahl al Oughlam, Morocco from Late Pliocene deposits. This shell thickness overlaps with the range observed in <italic>S</italic>. <italic>daberasensis</italic>, but the latter has a greater pore density comparable to <italic>Struthio camelus</italic>, thus following the general African trend [<xref ref-type="bibr" rid="pone.0164823.ref016">16</xref>]. The Late Pleistocene ostrich populations in India also show decreased shell thickness and high pore density similar to the African taxa. However, assignment of the Indian material to the African species based on morphological similarity alone may be problematic as this similarity may be a result of similar paleoenvironmental conditions [<xref ref-type="bibr" rid="pone.0164823.ref017">17</xref>–<xref ref-type="bibr" rid="pone.0164823.ref020">20</xref>]. DNA based species identification offers an alternative tool for authenticating the presence of ostriches in Indian peninsula.</p>
<p>Recently, avian fossil eggshells have been characterized as a potent source of preserved aDNA biomolecules [<xref ref-type="bibr" rid="pone.0164823.ref021">21</xref>]. The ultrastructure of eggshells shows a characteristic morphology consisting of several layers (mammillary layer, palisade layer and cuticle layer) beginning from the inner to outer surface (<xref ref-type="fig" rid="pone.0164823.g002">Fig 2</xref>). The matrix of the avian eggshell is porous, an ordered and heterogeneous complex, which is coupled with an extracellular structure composed of mineralized and non-mineralized regions having calcium carbonate (CaCO<sub>3</sub>) (calcite 97%) and 3.5% organic matrix [<xref ref-type="bibr" rid="pone.0164823.ref018">18</xref>,<xref ref-type="bibr" rid="pone.0164823.ref022">22</xref>–<xref ref-type="bibr" rid="pone.0164823.ref023">23</xref>]. The intracrystalline organic matrix of avian eggshells prevents diffusion losses and isotopic exchange of organic constituents [<xref ref-type="bibr" rid="pone.0164823.ref019">19</xref>,<xref ref-type="bibr" rid="pone.0164823.ref021">21</xref>,<xref ref-type="bibr" rid="pone.0164823.ref024">24</xref>]. This structural hallmark provides an external skeletal support to the developing avian embryo and controls gases, water exchange as well as guards against microbial interventions and physical stresses [<xref ref-type="bibr" rid="pone.0164823.ref021">21</xref>, <xref ref-type="bibr" rid="pone.0164823.ref025">25</xref>–<xref ref-type="bibr" rid="pone.0164823.ref027">27</xref>]. This structure also imparts stability as well as excellent mechanical strength to the eggshell, besides providing antibacterial and antifouling properties thus enhancing the preservation potential of biomolecules in eggshells [<xref ref-type="bibr" rid="pone.0164823.ref025">25</xref>,<xref ref-type="bibr" rid="pone.0164823.ref028">28</xref>].</p>
<fig id="pone.0164823.g002" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0164823.g002</object-id>
<label>Fig 2</label>
<caption>
<title>Ostrich Eggshell.</title>
<p>Pictorial view and radial cross section showing various layers, shell thickness and pore diameter.</p>
</caption>
<graphic mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0164823.g002" xlink:type="simple"/>
</fig>
<p>In this study, we examined the extent of aDNA preservation in eleven fossilized avian eggshell samples collected from eight archaeological sites in central and western parts of India (<xref ref-type="supplementary-material" rid="pone.0164823.s001">S1 Fig</xref>). The eggshells used in this investigation have survived harsh environmental conditions [<xref ref-type="bibr" rid="pone.0164823.ref005">5</xref>] and have been dated to Upper Palaeolithic (25,000 to 40,000 years B.P) based on the stone tool technology [<xref ref-type="bibr" rid="pone.0164823.ref006">6</xref>,<xref ref-type="bibr" rid="pone.0164823.ref029">29</xref>–<xref ref-type="bibr" rid="pone.0164823.ref030">30</xref>] and radiocarbon dating at Groningen University Lab, Germany [<xref ref-type="bibr" rid="pone.0164823.ref031">31</xref>,<xref ref-type="bibr" rid="pone.0164823.ref032">32</xref>]. Some of these eggshells are engraved into finished or unfinished beads [<xref ref-type="bibr" rid="pone.0164823.ref033">33</xref>–<xref ref-type="bibr" rid="pone.0164823.ref035">35</xref>]. The perforated beads belong to archaeological assemblages that correspond to the same time period when modern human populations would have dispersed out of Africa [<xref ref-type="bibr" rid="pone.0164823.ref036">36</xref>]. The presence of DNA hotspots in the mammillary layer was detected with DNA binding fluorescent dye Hoechst 33342 and imaged through Confocal Laser Scanning Microscopy. DNA was isolated, quantitized using Bioanalyzer 2100 (Agilent) and a 43 bp sequence of the mitochondrial coding region (16S rRNA) was amplified and sequenced from one of the samples using species-specific primers. The DNA sequence have been submitted to NCBI GenBank (accession number KU251475).</p>
</sec>
<sec id="sec002" sec-type="materials|methods">
<title>Materials and Methods</title>
<sec id="sec003" sec-type="materials|methods">
<title>Materials</title>
<p>Fossilized eggshell samples used in this study came from the excavation sites and personal museum collections of our team members. These eggshell fragments were recovered at excavation sites at Bundi (Rajasthan state, western India) and from the Pleistocene surface deposits of Gujarat, western India. Eggshell fragments, collected from eight different sites in India were analyzed (<xref ref-type="table" rid="pone.0164823.t001">Table 1</xref>). These samples have been assigned to the genus <italic>Struthio</italic> on the basis of morphological characteristics, especially thickness and pore pattern on the eggshells [<xref ref-type="bibr" rid="pone.0164823.ref006">6</xref>, <xref ref-type="bibr" rid="pone.0164823.ref029">29</xref>]</p>
<table-wrap id="pone.0164823.t001" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0164823.t001</object-id>
<label>Table 1</label> <caption><title>Sample collection sites and their respective states.</title></caption>
<alternatives>
<graphic id="pone.0164823.t001g" mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0164823.t001" xlink:type="simple"/>
<table>
<colgroup>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
</colgroup>
<thead>
<tr>
<th align="center" style="background-color:#FFFFFF">S.No.</th>
<th align="center" style="background-color:#FFFFFF">Site</th>
<th align="center" style="background-color:#FFFFFF">Area</th>
<th align="center" style="background-color:#FFFFFF">State &amp; GPS Co-ordinates</th>
<th align="center" style="background-color:#FFFFFF">Catalogue No.</th>
<th align="center" style="background-color:#FFFFFF">Dates</th>
<th align="center" style="background-color:#FFFFFF">Method</th>
<th align="center" style="background-color:#FFFFFF">Reference</th>
</tr>
</thead>
<tbody>
<tr>
<td align="center" rowspan="2" style="background-color:#FFFFFF"><bold>1</bold></td>
<td align="center" rowspan="2" style="background-color:#FFFFFF">Chavni Baroda</td>
<td align="center" rowspan="2" style="background-color:#FFFFFF">Bundi</td>
<td align="center" rowspan="2" style="background-color:#FFFFFF"><bold>Rajasthan</bold> (25° 26′ 24″ N, 75° 38′ 24″ E)</td>
<td align="center" style="background-color:#FFFFFF">SJ/CB/001</td>
<td align="center" rowspan="2" style="background-color:#FFFFFF">-</td>
<td align="center" rowspan="2" style="background-color:#FFFFFF">-</td>
<td align="center" rowspan="2" style="background-color:#FFFFFF">-</td>
</tr>
<tr>
<td align="center" style="background-color:#FFFFFF">SJ/CB/002</td>
</tr>
<tr>
<td align="center" style="background-color:#FFFFFF"><bold>2</bold></td>
<td align="center" style="background-color:#FFFFFF">Anjar</td>
<td align="center" style="background-color:#FFFFFF">Kachchh</td>
<td align="center" style="background-color:#FFFFFF"><bold>Gujarat</bold> (23° 6′ 48.32″ N, 70° 1′ 39.88″ E)</td>
<td align="center" style="background-color:#FFFFFF">SB/AN/003</td>
<td align="center" style="background-color:#FFFFFF">&gt;25000</td>
<td align="center" style="background-color:#FFFFFF"><bold>C</bold><sup><bold>14</bold></sup></td>
<td align="center" style="background-color:#FFFFFF"><bold>Sahni et al. (1989)</bold></td>
</tr>
<tr>
<td align="center" rowspan="2" style="background-color:#FFFFFF"><bold>3</bold></td>
<td align="center" rowspan="2" style="background-color:#FFFFFF">Chandresal- 1</td>
<td align="center" rowspan="2" style="background-color:#FFFFFF">Kota</td>
<td align="center" rowspan="2" style="background-color:#FFFFFF"><bold>Rajasthan</bold> (25° 10′ 48″ N, 75° 49′ 48″ E)</td>
<td align="center" style="background-color:#FFFFFF">GK/CH1/004</td>
<td align="center" rowspan="2" style="background-color:#FFFFFF">38900±<bold>750</bold></td>
<td align="center" rowspan="2" style="background-color:#FFFFFF"><bold>C</bold><sup><bold>14</bold></sup></td>
<td align="center" rowspan="2" style="background-color:#FFFFFF"><bold>Mishra (1995)</bold></td>
</tr>
<tr>
<td align="center" style="background-color:#FFFFFF">GK/CH1/005</td>
</tr>
<tr>
<td align="center" rowspan="2" style="background-color:#FFFFFF"><bold>4</bold></td>
<td align="center" rowspan="2" style="background-color:#FFFFFF">Chandresal -2</td>
<td align="center" rowspan="2" style="background-color:#FFFFFF">Kota</td>
<td align="center" rowspan="2" style="background-color:#FFFFFF"><bold>Rajasthan</bold> (25°13'34.1"N 75°55'34.7"E)</td>
<td align="center" style="background-color:#FFFFFF">GK/CH2/006</td>
<td align="center" rowspan="2" style="background-color:#FFFFFF">36500±<bold>600</bold></td>
<td align="center" rowspan="2" style="background-color:#FFFFFF"><bold>C</bold><sup><bold>14</bold></sup></td>
<td align="center" rowspan="2" style="background-color:#FFFFFF"><bold>Agarwal et al. (1991)</bold></td>
</tr>
<tr>
<td align="center" style="background-color:#FFFFFF">GK/CH2/007</td>
</tr>
<tr>
<td align="center" style="background-color:#FFFFFF"><bold>5</bold></td>
<td align="center" style="background-color:#FFFFFF">Nagda</td>
<td align="center" style="background-color:#FFFFFF">Chambal</td>
<td align="center" style="background-color:#FFFFFF"><bold>Rajasthan</bold> (25°10'53.3"N 75°48'45.5"E)</td>
<td align="center" style="background-color:#FFFFFF">GK/NA/008</td>
<td align="center" style="background-color:#FFFFFF">&gt;31000</td>
<td align="center" style="background-color:#FFFFFF"><bold>C</bold><sup><bold>14</bold></sup></td>
<td align="center" style="background-color:#FFFFFF"><bold>Mishra (1995)</bold></td>
</tr>
<tr>
<td align="center" style="background-color:#FFFFFF"><bold>6</bold></td>
<td align="center" style="background-color:#FFFFFF">Runija</td>
<td align="center" style="background-color:#FFFFFF">Ujjain</td>
<td align="center" style="background-color:#FFFFFF"><bold>Madhya Pradesh</bold> (23°09'39.4"N 75°16'02.6"E)</td>
<td align="center" style="background-color:#FFFFFF">GK/RU/009</td>
<td align="center" style="background-color:#FFFFFF">&gt;25000</td>
<td align="center" style="background-color:#FFFFFF"><bold>C</bold><sup><bold>14</bold></sup></td>
<td align="center" style="background-color:#FFFFFF"><bold>Mishra (1995)</bold></td>
</tr>
<tr>
<td align="center" style="background-color:#FFFFFF"><bold>7</bold></td>
<td align="center" style="background-color:#FFFFFF">Khajurna</td>
<td align="center" style="background-color:#FFFFFF">Bijora</td>
<td align="center" style="background-color:#FFFFFF"><bold>Madhya Pradesh</bold> (26°44'33.6"N 78°47'45.6"E)</td>
<td align="center" style="background-color:#FFFFFF">GK/KH/010</td>
<td align="center" style="background-color:#FFFFFF">&gt;25000</td>
<td align="center" style="background-color:#FFFFFF"><bold>C</bold><sup><bold>14</bold></sup></td>
<td align="center" style="background-color:#FFFFFF"><bold>Mishra (1995)</bold></td>
</tr>
<tr>
<td align="center" style="background-color:#FFFFFF"><bold>8</bold></td>
<td align="center" style="background-color:#FFFFFF">Ravishankar Nagar</td>
<td align="center" style="background-color:#FFFFFF">Bhopal</td>
<td align="center" style="background-color:#FFFFFF"><bold>Madhya Pradesh</bold> (23°15'34.8"N 77°24'24.9"E)</td>
<td align="center" style="background-color:#FFFFFF">GK/RN/011</td>
<td align="center" style="background-color:#FFFFFF">&gt;25000</td>
<td align="center" style="background-color:#FFFFFF"><bold>C</bold><sup><bold>14</bold></sup></td>
<td align="center" style="background-color:#FFFFFF"><bold>Mishra (1995)</bold></td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
</sec>
<sec id="sec004">
<title>Confocal Laser Scanning Microscopy</title>
<p>To determine the presence and location of DNA preserved in ratites eggshells, we used a confocal scanning microscope. The samples were incubated in square chambered plates with DNA binding fluorescent dye, Hoechst 33342 (1mg/ml) (Sigma, USA) at room temperature for 15 minutes and imaged under Leica TCS SP5 AOPS Confocal system. Images were captured (405 nm) using 10X and 40X (1.25 oil immersion) objectives, while 3-D sections snapshots were taken using Galvo-Z stage for Z sectioning. Cross sectional images of eggshells obtained were analysed using LASAF LITE software. Unstained eggshell samples were used as control.</p>
</sec>
<sec id="sec005">
<title>Quality control</title>
<p>A major concern throughout the present molecular analysis was contamination by any exogenous DNA. The DNA extraction was performed in a dedicated ancient DNA facility at CCMB, Hyderabad, India. This lab has already carried out investigations on human aDNA [<xref ref-type="bibr" rid="pone.0164823.ref037">37</xref>]. Disposable gloves, full body suits and face masks were used throughout the sample processing. The ancient DNA lab functions under positive pressure created through 5μM HEPA filters and is UV irradiated every night.</p>
<p>All samples and reagents were exposed to UV light. All equipments and benchtops were regularly cleaned with 20% bleach solution. Disposable laboratory wares, gloves and laboratory coats were used during the experiments. For each sample two independent extractions were prepared and done by two independent researchers. Every day, a single sample was extracted using controls. Controls consisted of PCR reaction mixtures without the sample. All the parameters were followed in order to avoid the risk of contamination [<xref ref-type="bibr" rid="pone.0164823.ref038">38</xref>–<xref ref-type="bibr" rid="pone.0164823.ref040">40</xref>].</p>
</sec>
<sec id="sec006">
<title>Sample preparation and DNA isolation</title>
<p>Eggshell samples were wiped with 5% hypochlorite solution followed by surface cleaning using 70% alcohol. A small piece of eggshell was cut from the source eggshell using Dremel MultiPro tools and around 700-800mg of powder was collected per sample. DNA extraction was performed on the powdered eggshell samples using a published protocol with slight modification [<xref ref-type="bibr" rid="pone.0164823.ref037">37</xref>,<xref ref-type="bibr" rid="pone.0164823.ref041">41</xref>]. Briefly, one gram of powder was incubated in a digestion buffer (700 mL per sample) containing final volumes of 0.47 M EDTA (pH 8.0), 20 mM Tris (pH 8.0), 1% Triton X-100, 10 mM Dithiothreitol (DTT) and 1 mg/mL proteinase K for up to 24 hours, followed by a final heating at 92°C. The heating aids in solubilisation of calcite and in releasing the DNA from crystalline matrix. The heating step was modified to prevent the complete denaturation of DNA. The solution was then concentrated with 30,000 KDa MWCO columns (Millipore) and purified using commercial silica spin columns (Qiagen) with slight change in manufacturer’s protocol i.e incubating the elution step for 5 minutes at 55 degree centigrade.</p>
</sec>
<sec id="sec007">
<title>Ancient DNA quantity</title>
<p>Eluted DNA samples were quantified using the high sensitivity chip on the Bioanalyzer 2100 (Agilent) and prepared for PCR amplification and sequencing. The quantity of total DNA was estimated to be close to the lowest detection limit of the Bioanalyzer.</p>
</sec>
<sec id="sec008">
<title>PCR amplification and sequencing</title>
<p>For amplification of 16S rRNA region and control region of mitochondrial DNA, primers were designed as shown in <xref ref-type="table" rid="pone.0164823.t002">Table 2</xref>. Same fragments from each sample was amplified in three DNA extracts to validate results. Annealing temperature was changed and the annealing time was extended. Primers were designed using Primer 3 plus and using mtDNA reference sequence of <italic>Struthio camelus</italic>. As DNA quantity was extremely low, Titanium Taq DNA Polymerase (Roche) was used, which is mostly used for High-Throughput genotyping using MassArray System (Agena Biosciences).</p>
<table-wrap id="pone.0164823.t002" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0164823.t002</object-id>
<label>Table 2</label> <caption><title>Primer sequences used for PCR and sequencing of 16S and Hypervariable regions of mitochondrial DNA.</title></caption>
<alternatives>
<graphic id="pone.0164823.t002g" mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0164823.t002" xlink:type="simple"/>
<table>
<colgroup>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
</colgroup>
<thead>
<tr>
<th align="center">Primer name</th>
<th align="center">Sequences</th>
<th align="center">Gene</th>
</tr>
</thead>
<tbody>
<tr>
<td align="center">2F</td>
<td align="center"><monospace>5'GTACCGCAAGGGAAAGATGA3'</monospace></td>
<td align="center" rowspan="2">16S rRNA</td>
</tr>
<tr>
<td align="center">2R</td>
<td align="center"><monospace>5'GGATGGCAAGCTTAAATTCG3'</monospace></td>
</tr>
<tr>
<td align="center">4F</td>
<td align="center"><monospace>5'GGGGGTTACCCTCTAATGGA3'</monospace></td>
<td align="center" rowspan="2">16S rRNA</td>
</tr>
<tr>
<td align="center">4R</td>
<td align="center"><monospace>5'AGCTGGTTGCCTGTGAAAAG3'</monospace></td>
</tr>
<tr>
<td align="center">7F</td>
<td align="center"><monospace>5'CGACTCAGGAGCGCCTATTA3'</monospace></td>
<td align="center" rowspan="2">16S rRNA</td>
</tr>
<tr>
<td align="center">7R</td>
<td align="center"><monospace>5'TGGCTGAAGGCTATGTTTTTG3'</monospace></td>
</tr>
<tr>
<td align="center">1CF</td>
<td align="center"><monospace>5'TCGCGATTAAGAGGGACAAT3'</monospace></td>
<td align="center" rowspan="2">Hypervariable Region</td>
</tr>
<tr>
<td align="center">1CR</td>
<td align="center"><monospace>5'TCCATTCACGTTCCCCTTTA3'</monospace></td>
</tr>
</tbody>
</table>
</alternatives>
</table-wrap>
<p>PCR reaction was set up in 5 μl total volume per reaction. The reaction mixture contained 2 μl of eluted DNA, reaction buffer supplied by the manufacturer (Agena Biosciences), dNTPs at 200 μM each, Taq polymerase (Roche) at 0.1 unit (0.02 μl at 5 U/μl, primers at 200 nM and MgCl<sub>2</sub> at the final concentration as mentioned in the protocol. PCR cycling conditions for Titanium Taq DNA polymerase (Invitrogen) included: initial heating to 95°C, followed by 45 cycles of 94°C for 20 sec, 56°C for 30 sec, and 72°C for 1 min, and a final step of 72°C for 3 min. After PCR, the products were checked on 3% agarose gel (Sigma Aldrich) using 50bp DNA ladders (NEB) and the expected product sizes were confirmed. PCR products were purified using Shrimp Alkaline Phosphatase (SAP) from Agena Biosciences as per the manufacturer protocol.</p>
</sec>
<sec id="sec009">
<title>Cloning</title>
<p>PCR products were cloned using the pMOS Blunt-ended PCR cloning kit (GE Healthcare, USA). PCR products were ligated to pMossBlue dephosphorylated blunt vector as per the manufacturer protocol. 2 μl of ligated product and DH5α competent cells were used for transformation. Sample to vector ratio was adjusted to optimize insert-specific transformed cells. Transformed cells were spread on LB agar media and plate was incubated at 37°C overnight. For screening, white blue selection was made, Almost half of the colonies were directly used for colony PCR and for the remaining half plasmid DNA was isolated using plasmid isolation kit (Qiagen).</p>
</sec>
<sec id="sec010">
<title>DNA sequencing</title>
<p>The purified plasmids were sequenced using M13 forward primer as well as insert specific forward and reverse primers using Sanger’s dideoxy chain terminator cycle sequencing method. In case of colony PCR, PCR products were purified by treating with Exonuclease I and Shrimp Alkaline Phosphatase (ExoSAP-IT®; USB Corporation, Cleveland, Ohio, USA) and incubated at 37°C for 15 min. then at 80°C for 15 min. Purified PCR products (1.0 μl) were subjected to sequencing reaction, by adding 0.65–2.0 pmoles of the primer and 3.2 μl of BigDye<sup>TM</sup> (Applied Biosystems, Foster City, CA, USA; containing fluorescently labeled ddNTPs and unlabeled dNTP). The sequencing PCR conditions include: 30 cycles of 94°C for 10 sec, 55°C for 5 sec and 60°C for 4 min. The PCR product was precipitated using 3 M sodium acetate (pH 5.2) and ethanol, and centrifuged at 4,000 rpm at 4°C for 15 min. The pellet was washed in 70% alcohol, air dried and suspended in 10 μl of HiDi<sup>TM</sup> formamide (Applied Biosystems, Foster city, USA) and loaded on to ABI 3730 Automated DNA analyzer, after denaturation at 94°C for 10 min. The samples were run using POP-7<sup>TM</sup> polymer and analyzed using ‘Sequencing Analysis’ software (Applied Biosystems, Foster city, USA) and MEGA 6.</p>
</sec>
<sec id="sec011">
<title>Negative Control Amplification</title>
<p>Deionized water was used as negative controls for the PCR amplification instead of the DNA extracts, following the same protocols and using the same reagents. No amplification of the negative control reactions confirmed the absence of lab-based contamination.</p>
</sec>
</sec>
<sec id="sec012" sec-type="results">
<title>Results</title>
<p>Eleven fossil samples from India, collected from eight archaeological sites (<xref ref-type="table" rid="pone.0164823.t001">Table 1</xref>) and from the personal collections of palaeontologists, were investigated to identify the extent of DNA preservation in these fossil avian eggshells and for DNA based species identification. To maximize the DNA recovery from fossil eggshells, it was important to determine the physical location of DNA and to microscopically identify its presence in the inner, outer or calcified layers of the eggshell.</p>
<sec id="sec013">
<title>CLSM of eggshells</title>
<p>DNA in ratite eggshells is distributed uniformly throughout the eggshell matrix but may be concentrated around the periphery of the mammillary cones [<xref ref-type="bibr" rid="pone.0164823.ref021">21</xref>]. Confocal imaging of the avian eggshells studied here demonstrated the distribution of DNA throughout the eggshell matrix, as evident from the DNA hotspots observed after staining with Hoechst 33342, a fluorescent dye used to stain DNA (<xref ref-type="fig" rid="pone.0164823.g003">Fig 3A and 3B</xref>). The control image of unstained eggshells did not show any hotspots, confirming that the fluorescence observed in sample images was due to the binding of dye with DNA. The inner surface images of the eggshell fragments showed the concentration of DNA hotspots on the periphery of mammillary cones (<xref ref-type="fig" rid="pone.0164823.g003">Fig 3C</xref>). Three dimensional images captured using Z stage at different angles clearly validate the presence of DNA hotspots distributed uniformly throughout the matrix (<xref ref-type="fig" rid="pone.0164823.g003">Fig 3D and 3E</xref>).</p>
<fig id="pone.0164823.g003" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0164823.g003</object-id>
<label>Fig 3</label>
<caption>
<title/>
<p><bold>Confocal images of ratite eggshell</bold> (A) Control sample, without staining showing no fluoroscence (B) Sample stained with Hoechst 33342, displaying the DNA distributed throughout the matrix. Arrow indicates the presence of hotspots in fossilized eggshells (C) Fluorescently labelled DNA can be observed on the periphery of mammillary cones. Scale Bar 50 μm. (D) Three-dimensional imaging of inner layers of eggshell shows DNA distributed throughout the matrix as observed through fluorescent hotspots in horizontal view (E) Vertical view [Inset- Picture of eggshell. Marked area was observed under CLSM].</p>
</caption>
<graphic mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0164823.g003" xlink:type="simple"/>
</fig>
</sec>
<sec id="sec014">
<title>DNA extraction and species determination</title>
<p>Five samples (GK/RN/011, GK/RU/009, GK/CH1/004, GK/KH/010 and GK/CH2/007), from five different sites, yielded aDNA (<xref ref-type="table" rid="pone.0164823.t003">Table 3</xref>). Since the copy number of DNA is traditionally low in ancient samples and aDNA is fragmented in nature, we optimized the amplification protocol targeting the mitochondrial hypervariable region (HVR) and 16S rRNA gene by using PCR primers designed to amplify short fragments, using multiple aliquots of the DNA extracts as template for the PCR reactions. The extractions and PCR set up were done at CCMB, Hyderabad, India in facilities dedicated to the analysis of low copy number template. The DNA extracts were visualized on the Bioanalyzer 2100 (Agilent). Bioanalyzer revealed a concentration of &lt;100pg/ml in the five samples GK/RN/011 (123.6 pg/ml), GK/KH/010 (566.32 pg/ml) GK/CH1/004 (749.06 pg/ml), GK/RU/009 (316.94 pg/ml) and GK/CH2/007 (223.84 pg/ml). Quantification of DNA extracted from the fossil eggshell samples showed a strong band around 45 bp in GK/RN/011 and 65bp in GK/KH/010 (<xref ref-type="fig" rid="pone.0164823.g004">Fig 4</xref>). GK/RN/011 was further used for phylogentic analysis. Bioanalyzer images of other samples showed multiple bands and smears which could be due to degraded nature of aDNA. The extracted and amplified DNA were further cloned using pMOS blunt ended PCR cloning kit and multiple clones were sequenced.</p>
<fig id="pone.0164823.g004" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0164823.g004</object-id>
<label>Fig 4</label>
<caption>
<title>Quantification of fossilized DNA.</title>
<p>Bioanalyzer image showing bands of 65 bp and 43 bp in GK/KH/010 and GK/RN/011 respectively. Concentration of DNA in pg/ml of five samples, calculated by bioanalyzer.</p>
</caption>
<graphic mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0164823.g004" xlink:type="simple"/>
</fig>
<table-wrap id="pone.0164823.t003" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0164823.t003</object-id>
<label>Table 3</label> <caption><title>List of eggshell specimens from India with associated genes and predicted taxon.</title></caption>
<alternatives>
<graphic id="pone.0164823.t003g" mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0164823.t003" xlink:type="simple"/>
<table>
<colgroup>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
<col align="left" valign="middle"/>
</colgroup>
<thead>
<tr>
<th align="left">Site</th>
<th align="left">Specimen ID</th>
<th align="left">Dates</th>
<th align="left">Predicted taxa</th>
<th align="left">Gene</th>
<th align="left">Amplicon site</th>
<th align="left">Fragment length (bp)</th>
<th align="left">E value (NCBI)</th>
<th align="left">Identity</th>
<th align="left">Taxon</th>
<th align="left">Common name</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left">Ravishankar Nagar</td>
<td align="left">GK/RN/011</td>
<td align="left">&gt;25000</td>
<td align="left"><italic>Struthio Camelus</italic></td>
<td align="left">16s RNA</td>
<td align="left">1255–1298</td>
<td align="left">43</td>
<td align="left">2e-07</td>
<td align="left">92%</td>
<td align="left"><italic>S</italic>. <italic>Camelus</italic></td>
<td align="left">Ostrich</td>
</tr>
<tr>
<td align="left">Runija</td>
<td align="left">GK/RU/009</td>
<td align="left">&gt;25000</td>
<td align="left"><italic>S</italic>. <italic>Camelus</italic></td>
<td align="left">16s RNA</td>
<td align="left">1931–1954</td>
<td align="left">23</td>
<td align="left">6e-09</td>
<td align="left">100%</td>
<td align="left"><italic>S</italic>. <italic>Camelus</italic></td>
<td align="left">Ostrich</td>
</tr>
<tr>
<td align="left">Chandresal-1</td>
<td align="left">GK/CH1/004</td>
<td align="left">38,900 ± 750 BP</td>
<td align="left"><italic>S</italic>. <italic>Camelus</italic></td>
<td align="left">16s RNA</td>
<td align="left">4859–4882</td>
<td align="left">23</td>
<td align="left">2e-07</td>
<td align="left">96%</td>
<td align="left"><italic>Apteryx mantellis</italic></td>
<td align="left">Kiwi</td>
</tr>
<tr>
<td align="left">Khajurna</td>
<td align="left">GK/KH/010</td>
<td align="left">&gt;25000</td>
<td align="left"><italic>S</italic>. <italic>Camelus</italic></td>
<td align="left">HVR</td>
<td align="left">15586–15606</td>
<td align="left">20</td>
<td align="left">0.019</td>
<td align="left">89%</td>
<td align="left"><italic>S</italic>. <italic>Camelus</italic></td>
<td align="left">Ostrich</td>
</tr>
<tr>
<td align="left">Chandresal-2</td>
<td align="left">GK/CH2/007</td>
<td align="left">36,550 ± 600 BP</td>
<td align="left"><italic>S</italic>. <italic>Camelus</italic></td>
<td align="left">HVR</td>
<td align="left">15707–15724</td>
<td align="left">17</td>
<td align="left">1e-05</td>
<td align="left">87%</td>
<td align="left"><italic>S</italic>. <italic>Camelus</italic></td>
<td align="left">Ostrich</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot>
<fn id="t003fn001"><p>Dates according to Kumar et al, 1988; 1990</p></fn>
<fn id="t003fn002"><p>Taxon based on the closest GenBank BLAST match</p></fn>
<fn id="t003fn003"><p>Predicted taxa on the basis of eggshell morphology and location</p></fn>
<fn id="t003fn004"><p>HVR Hypervariable Region</p></fn>
</table-wrap-foot>
</table-wrap>
<p>A fragment length of 43 bp and 23 bp, corresponding to GK/RN/011 and GK/RU/009 respectively, was obtained from 25,000 years old samples collected from Ravishankar Nagar (GK/RN/011) and Runija (GK/RU/009), respectively, which showed 92% and 100% similarity with <italic>Struthio camelus</italic>. Alignments of partial 16S rRNA gene sequence of ratites (obtained from NCBI) and samples GK/RN/011 and GK/RU/009 were done (<xref ref-type="fig" rid="pone.0164823.g005">Fig 5</xref>). Other samples GK/CH1/004, GK/KH/010 and GK/CH2/007 collected from Chandresal 1, Khajurna and Chandresal 2 also attained a fragment length of 23 bp, 20 bp and 17 bp, respectively. Obtained sequences were compared with published partial mitochondrial sequences on NCBI. BLAST (NCBI) and multiple sequence alignment showed 92% similarity with <italic>Struthio</italic> species.</p>
<fig id="pone.0164823.g005" position="float">
<object-id pub-id-type="doi">10.1371/journal.pone.0164823.g005</object-id>
<label>Fig 5</label>
<caption>
<title>Alignment of partial 16s rRNA gene sequences of samples and ratites and evolutionary relationships of taxa.</title>
<p>(A) and (B) Alignment of ratites sequences and sample GK/RN/011 and GK/RU/009. Dots signify identity to the <italic>Struthio camelus</italic> sequence. (C) The evolutionary history was inferred using the Neighbor-Joining method [<xref ref-type="bibr" rid="pone.0164823.ref042">42</xref>]. The bootstrap consensus tree inferred from 1000 replicates [<xref ref-type="bibr" rid="pone.0164823.ref043">43</xref>] is taken to represent the evolutionary history of the taxa analyzed [<xref ref-type="bibr" rid="pone.0164823.ref043">43</xref>]. The evolutionary distances were computed using the Maximum Composite Likelihood method [<xref ref-type="bibr" rid="pone.0164823.ref044">44</xref>] and are in the units of the number of base substitutions per site. The analysis involved 7 nucleotide sequences. The codon positions which were included: 1st+2nd+3rd+Noncoding. All positions containing gaps and missing data were eliminated. There were a total of 87 positions in the final dataset. Evolutionary analyses were conducted in MEGA6 [<xref ref-type="bibr" rid="pone.0164823.ref045">45</xref>].</p>
</caption>
<graphic mimetype="image" position="float" xlink:href="info:doi/10.1371/journal.pone.0164823.g005" xlink:type="simple"/>
</fig>
</sec>
<sec id="sec015">
<title>Phylogenetic analysis</title>
<p>To assess the phylogenetic relationship of our DNA sequences, a comparative study with published avian mtDNA sequences was done and the phylogenetic tree was reconstructed (<xref ref-type="fig" rid="pone.0164823.g005">Fig 5</xref>) using MEGA 6.06. The phylogenetic analysis was based on maximum parsimony, neighbour joining and maximum likelihood [<xref ref-type="bibr" rid="pone.0164823.ref042">42</xref>–<xref ref-type="bibr" rid="pone.0164823.ref045">45</xref>]. These three methods are based on different assumptions and, thus, can be used complementarily to confirm the obtained phylogeny. The neighbour joining method was further extended and bootstrap consensus tree was inferred from 1000 replicates. The sample sequence shows phylogenetical similarity to the ratites family, forming a sister group with S. <italic>camelus</italic> (<xref ref-type="fig" rid="pone.0164823.g005">Fig 5</xref>). Phylogenetic tree showed the ostrich clade as basal to studied palaeognaths, which is consistent with the vicariant speciation model and other molecular phylogeny studies on ratites [<xref ref-type="bibr" rid="pone.0164823.ref001">1</xref>,<xref ref-type="bibr" rid="pone.0164823.ref003">3</xref>]. The estimates of evolutionary divergence between sequences were conducted in MEGA 6 (<xref ref-type="supplementary-material" rid="pone.0164823.s002">S1 Table</xref>) using the Maximum Composite Likelihood model.</p>
</sec>
</sec>
<sec id="sec016" sec-type="conclusions">
<title>Discussion</title>
<p>The breakup of the supercontinent Gondwanaland during the Early Cretaceous resulted in the vicariant speciation of ratites, the flightless birds, to various countries that we know today as South America, Australia, New Zealand, India and Madagascar. The continental drift theory implies that ostriches migrated to India from other landmasses, consistent with the geographical distribution of fossil records of ratites. The fossil eggshells used in this study date to 25–40 ka and morphologically characterized as belonging to <italic>Struthio camelus</italic> [<xref ref-type="bibr" rid="pone.0164823.ref006">6</xref>,<xref ref-type="bibr" rid="pone.0164823.ref030">30</xref>]. Ostriches were believed to be extinct in India during the Late Pleistocene but our investigation reports the first genetic evidence of their presence on the subcontinent. Fossil eggshells have been reported as being more conducive to biomolecular preservation compared to fossil bones due to their intracrystalline structure, which minimizes microbial contamination [<xref ref-type="bibr" rid="pone.0164823.ref025">25</xref>,<xref ref-type="bibr" rid="pone.0164823.ref028">28</xref>]. This study is significant as it is the first report, to the best of our knowledge, of long-term DNA preservation in fossil eggshells collected from tropical environments with heterogeneous climatic conditions as those encountered in India. aDNA survival in fossils is usually determined by the thermal age theory[<xref ref-type="bibr" rid="pone.0164823.ref046">46</xref>]. This theory requires specific parameters including temperature, pressure and climatic conditions of fossil locations. These parameters are unfortunately not available for our samples, which are from open sites. To evaluate the presence of DNA, Hoechst 33342 stained fossil eggshell samples were visualized under CLSM. DNA hotspots observed on the periphery of mammillary layer are consistent with shape and size of epithelial fossil cells instead of random shapeless patches, eliminating the possibility of these hotspots being of microbial origin. Modifications in the aDNA isolation procedure of Oskam and Bunce [<xref ref-type="bibr" rid="pone.0164823.ref041">41</xref>] were made at the heating step to prevent complete denaturation of aDNA. Species-specific primers based on <italic>S</italic>. <italic>camelus</italic> were designed to target short fragments. Amplified fragments were further cloned and sequenced. The multiple alignment of sequences confirms the relationship of our samples with other avian species and their identity with <italic>Stuthio camelus</italic>. Molecular phylogenetic tree further revealed the placement of our samples in this clade and a basal position relative to other studied paleognaths, which supports the vicariant speciation theory and other ratites phylogenetic studies.</p>
<p>Morphological studies conducted earlier on these eggshells [<xref ref-type="bibr" rid="pone.0164823.ref006">6</xref>,<xref ref-type="bibr" rid="pone.0164823.ref017">17</xref>,<xref ref-type="bibr" rid="pone.0164823.ref030">30</xref>] indicated their close proximity to <italic>Struthio</italic> species. The present genetic study performed using species-specific primers confirms that the eggshell fragments belong to the <italic>Struthio</italic> species and that these ratites were present in India during the Late Pleistocene. Paleoenvironmental studies indicate favorable environment for survival and expansion of this species in the Indian subcontinent during Late Pleistocene [<xref ref-type="bibr" rid="pone.0164823.ref005">5</xref>]. Currently, ostriches are extinct in India and their sudden extinction coincides with the period of expansion of humans from Africa towards South Asia [<xref ref-type="bibr" rid="pone.0164823.ref005">5</xref>]. There is evidence in the archaeological context for the occurrence of ostrich eggshell beads and crisscross motif patterns on these eggshells with striking similarities to material culture of Africa [<xref ref-type="bibr" rid="pone.0164823.ref036">36</xref>]. These records suggest a connection between the earliest modern humans in southern Asia and their probable ancestors in eastern and southern Africa [<xref ref-type="bibr" rid="pone.0164823.ref036">36</xref>]. However, no conclusive evidence in this regard is available as yet.</p>
<p>The archaeological records of eggshells fragments, beads and rock shelter paintings discovered at Indian sites are consistent with our genetic study, which is the first molecular evidence indicating the presence of ostriches in India. These evidences further endorse the biogeographical movement of ostriches in India when Indo-Madagascar broke apart [<xref ref-type="bibr" rid="pone.0164823.ref003">3</xref>]. Furthermore, this genetic characterization of fossil eggshells is important in studying aspects related to biology, ecology and extinction of this species.</p>
</sec>
<sec id="sec017">
<title>Supporting Information</title>
<supplementary-material id="pone.0164823.s001" mimetype="image/tiff" position="float" xlink:href="info:doi/10.1371/journal.pone.0164823.s001" xlink:type="simple">
<label>S1 Fig</label>
<caption>
<title>Sample collection locations shown on map.</title>
<p>(A) Bundi (B) Anjar (C) Chandresal-1 (D) Nagda (E) Runija (F) Khajurna (G) Ravishankar nagar. (Source- Maps at the CIA (public domain): <ext-link ext-link-type="uri" xlink:href="https://www.cia.gov/library/publications/the-world-factbook/index.html" xlink:type="simple">https://www.cia.gov/library/publications/the-world-factbook/index.html</ext-link>).</p>
<p>(TIFF)</p>
</caption>
</supplementary-material>
<supplementary-material id="pone.0164823.s002" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" position="float" xlink:href="info:doi/10.1371/journal.pone.0164823.s002" xlink:type="simple">
<label>S1 Table</label>
<caption>
<title>Estimates of Evolutionary Divergence between Sequences.</title>
<p>The number of base substitutions per site from between sequences are shown. Analyses were conducted using the Maximum Composite Likelihood model. The analysis involved 7 nucleotide sequences. Codon positions included were 1st+2nd+3rd+Noncoding. All positions containing gaps and missing data were eliminated. There were a total of 87 positions in the final dataset. Evolutionary analyses were conducted in MEGA6.</p>
<p>(DOCX)</p>
</caption>
</supplementary-material>
</sec>
</body>
<back>
<ack>
<p>We acknowledge the much-needed help received from the All India Institute of Medical Sciences (AIIMS), New Delhi in conducting Confocal Laser Scanning Microscopy imaging. We are thankful to CCMB Director Ch. Mohan Rao, to allow us to perform experiments in his Ancient DNA laboratory. For the financial support, we would like to acknowledge the Council for Scientific and Industrial Research and University Grants Commission (UGC), New Delhi. We would like to thank Mr. Om Prakash Sharma for helping in procuring samples from Bundi site.</p>
</ack>
<ref-list>
<title>References</title>
<ref id="pone.0164823.ref001"><label>1</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Mitchell</surname> <given-names>KJ</given-names></name>, <name name-style="western"><surname>Llamas</surname> <given-names>B</given-names></name>, <name name-style="western"><surname>Soubrier</surname> <given-names>J</given-names></name>, <name name-style="western"><surname>Rawlence</surname> <given-names>NJ</given-names></name>, <name name-style="western"><surname>Worthy</surname> <given-names>TH</given-names></name>, <name name-style="western"><surname>Wood</surname> <given-names>J</given-names></name>, <etal>et al</etal>. <article-title>Ancient DNA reveals elephant birds and kiwi are sister taxa and clarifies ratite bird evolution</article-title>. <source>Science</source>. <year>2014</year> <month>May</month> <day>23</day>;<volume>344</volume>(<issue>6186</issue>):<fpage>898</fpage>–<lpage>900</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1126/science.1251981" xlink:type="simple">10.1126/science.1251981</ext-link></comment> <object-id pub-id-type="pmid">24855267</object-id></mixed-citation></ref>
<ref id="pone.0164823.ref002"><label>2</label><mixed-citation publication-type="book" xlink:type="simple"><name name-style="western"><surname>Dawkins</surname> <given-names>R</given-names></name>. <chapter-title>The ancestor's tale: a pilgrimage to the dawn of evolution</chapter-title>. <source>New York Times Book Review</source>. <volume>688</volume>, <publisher-name>Houghton Mifflin</publisher-name>, <year>2004</year>.</mixed-citation></ref>
<ref id="pone.0164823.ref003"><label>3</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Cooper</surname> <given-names>A</given-names></name>, <name name-style="western"><surname>Lalueza-Fox</surname> <given-names>C</given-names></name>, <name name-style="western"><surname>Anderson</surname> <given-names>S</given-names></name>, <name name-style="western"><surname>Rambaut</surname> <given-names>A</given-names></name>, <name name-style="western"><surname>Austin</surname> <given-names>J</given-names></name>, <name name-style="western"><surname>Ward</surname> <given-names>R</given-names></name>. <article-title>Complete mitochondrial genome sequences of two extinct moas clarify ratite evolution</article-title>. <source>Nature</source>. <year>2001</year> <month>Feb</month> <day>8</day>;<volume>409</volume>(<issue>6821</issue>):<fpage>704</fpage>–<lpage>7</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1038/35055536" xlink:type="simple">10.1038/35055536</ext-link></comment> <object-id pub-id-type="pmid">11217857</object-id></mixed-citation></ref>
<ref id="pone.0164823.ref004"><label>4</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Roberts</surname> <given-names>P</given-names></name>, <name name-style="western"><surname>Delson</surname> <given-names>E</given-names></name>, <name name-style="western"><surname>Miracle</surname> <given-names>P</given-names></name>, <name name-style="western"><surname>Ditchfield</surname> <given-names>P</given-names></name>, <name name-style="western"><surname>Roberts</surname> <given-names>RG</given-names></name>, <name name-style="western"><surname>Jacobs</surname> <given-names>Z</given-names></name>, <etal>et al</etal>. <article-title>Continuity of mammalian fauna over the last 200,000 y in the Indian subcontinent</article-title>. <source>Proceedings of the National Academy of Sciences</source>. <year>2014</year> <month>Apr</month> <day>22</day>;<volume>111</volume>(<issue>16</issue>):<fpage>5848</fpage>–<lpage>53</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref005"><label>5</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Blinkhorn</surname> <given-names>J</given-names></name>, <name name-style="western"><surname>Achyuthan</surname> <given-names>H</given-names></name>, <name name-style="western"><surname>Petraglia</surname> <given-names>MD</given-names></name>. <article-title>Ostrich expansion into India during the Late Pleistocene: Implications for continental dispersal corridors</article-title>. <source>Palaeogeography, Palaeoclimatology, Palaeoecology</source>. <year>2015</year> <month>Jan</month> <day>1</day>;<volume>417</volume>:<fpage>80</fpage>–<lpage>90</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref006"><label>6</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Sahni</surname> <given-names>A</given-names></name>, <name name-style="western"><surname>Kumar</surname> <given-names>G</given-names></name>, <name name-style="western"><surname>Bajpai</surname> <given-names>S</given-names></name>, <name name-style="western"><surname>Srinivasan</surname> <given-names>S</given-names></name>. <article-title>A review of Late Pleistocene ostriches (Struthio sp.) in India</article-title>. <source>Man and Environment</source>. <year>1990</year>;<volume>15</volume>(<issue>1</issue>):<fpage>41</fpage>–<lpage>7</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref007"><label>7</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Badam</surname> <given-names>GL</given-names></name>. <article-title>A note on the ostrich in India since the Miocene</article-title>. <source>Man and Environment</source>. <year>2005</year>;<volume>30</volume>(<issue>2</issue>):<fpage>97</fpage>–<lpage>104</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref008"><label>8</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Davies</surname> <given-names>W</given-names></name>. <article-title>III.—On some Fossil Bird-Remains from the Siwalik Hills in the British Museum</article-title>. <source>Geological Magazine (Decade II)</source>. <year>1880</year>;<volume>7</volume>(<issue>01</issue>):<fpage>18</fpage>–<lpage>27</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref009"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Manegold A, Louchart A, Carrier J, Elzanowski A. The early Pliocene avifauna of Langebaanweg (South Africa): a review and update. In Paleontological Research 2013—Proceedings of the 8th International Meeting of the Society of Avian Palaeontology and Evolution 2013 (pp. 135–152).</mixed-citation></ref>
<ref id="pone.0164823.ref010"><label>10</label><mixed-citation publication-type="book" xlink:type="simple"><name name-style="western"><surname>Leakey</surname> <given-names>MG</given-names></name>, <name name-style="western"><surname>Harris</surname> <given-names>JM</given-names></name>. <chapter-title>Lothagam: the dawn of humanity in eastern Africa</chapter-title>. <publisher-name>Columbia University Press</publisher-name>; <year>2003</year>.</mixed-citation></ref>
<ref id="pone.0164823.ref011"><label>11</label><mixed-citation publication-type="book" xlink:type="simple"><name name-style="western"><surname>Pickford</surname> <given-names>M</given-names></name>, <name name-style="western"><surname>Senut</surname> <given-names>B</given-names></name>. <source>Geology and Palaeobiology of the Central and Southern Namib Desert, Southwestern Africa: Geology and History of Study</source>. <publisher-name>Geological survey of Namibia</publisher-name>; <year>2000</year>.</mixed-citation></ref>
<ref id="pone.0164823.ref012"><label>12</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Sauer</surname> <given-names>EG</given-names></name>, <name name-style="western"><surname>Sauer</surname> <given-names>EM</given-names></name>. <article-title>Ratite Eggshell Fragments from Mio-Pliocene-Continental Sediments in the District of Ouarzazate, Morocco</article-title>. <source>Palaeontographica Abteilung A</source>. <year>1978</year> <month>Jan</month> <volume>1</volume>:<fpage>1</fpage>–<lpage>54</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref013"><label>13</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Stidham</surname> <given-names>TA</given-names></name>. <article-title>Extinct ostrich eggshell (Aves: Struthionidae) from the Pliocene Chiwondo Beds, Malawi: implications for the potential biostratigraphic correlation of African Neogene deposits</article-title>. <source>Journal of human evolution</source>. <year>2004</year> <month>Apr</month> <day>30</day>;<volume>46</volume>(<issue>4</issue>):<fpage>489</fpage>–<lpage>96</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1016/j.jhevol.2004.02.002" xlink:type="simple">10.1016/j.jhevol.2004.02.002</ext-link></comment> <object-id pub-id-type="pmid">15066381</object-id></mixed-citation></ref>
<ref id="pone.0164823.ref014"><label>14</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Stern</surname> <given-names>LA</given-names></name>, <name name-style="western"><surname>Johnson</surname> <given-names>GD</given-names></name>, <name name-style="western"><surname>Chamberlain</surname> <given-names>CP</given-names></name>. <article-title>Carbon isotope signature of environmental change found in fossil ratite eggshells from a South Asian Neogene sequence</article-title>. <source>Geology</source>. <year>1994</year> <month>May</month> <day>1</day>;<volume>22</volume>(<issue>5</issue>):<fpage>419</fpage>–<lpage>22</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref015"><label>15</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Harrison</surname> <given-names>T</given-names></name>, <name name-style="western"><surname>Msuya</surname> <given-names>CP</given-names></name>. <article-title>Fossil struthionid eggshells from Laetoli, Tanzania: Taxonomic and biostratigraphic significance</article-title>. <source>Journal of African Earth Sciences</source>. <year>2005</year> <month>Apr</month> <day>30</day>;<volume>41</volume>(<issue>4</issue>):<fpage>303</fpage>–<lpage>15</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref016"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Mourer-Chauviré C, Geraads D. The Struthionidae and Pelagornithidae (Aves: Struthioniformes, Odontopterygiformes) from the late Pliocene of Ahl al Oughlam, Morocco. In6ème Symposium international de la Society for Avian Paleontology and Evolution (SAPE) 2005 (Vol. 7, pp. 169–194).</mixed-citation></ref>
<ref id="pone.0164823.ref017"><label>17</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Patnaik</surname> <given-names>R</given-names></name>, <name name-style="western"><surname>Sahni</surname> <given-names>A</given-names></name>, <name name-style="western"><surname>Cameron</surname> <given-names>D</given-names></name>, <name name-style="western"><surname>Pillans</surname> <given-names>B</given-names></name>, <name name-style="western"><surname>Chatrath</surname> <given-names>P</given-names></name>, <name name-style="western"><surname>Simons</surname> <given-names>E</given-names></name>, <etal>et al</etal>. <article-title>Ostrich-like eggshells from a 10.1 million-yr-old Miocene ape locality, Haritalyangar, Himachal Pradesh, India</article-title>. <source>Current Science</source>. <year>2009</year> <month>Jun</month> <day>10</day>;<volume>96</volume>(<issue>11</issue>):<fpage>1485</fpage>–<lpage>95</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref018"><label>18</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Hincke</surname> <given-names>MT</given-names></name>, <name name-style="western"><surname>Gautron</surname> <given-names>J</given-names></name>, <name name-style="western"><surname>Panheleux</surname> <given-names>M</given-names></name>, <name name-style="western"><surname>Garcia-Ruiz</surname> <given-names>J</given-names></name>, <name name-style="western"><surname>McKee</surname> <given-names>MD</given-names></name>, <name name-style="western"><surname>Nys</surname> <given-names>Y</given-names></name>. <article-title>Identification and localization of lysozyme as a component of eggshell membranes and eggshell matrix</article-title>. <source>Matrix Biology</source>. <year>2000</year> <month>Sep</month> <day>30</day>;<volume>19</volume>(<issue>5</issue>):<fpage>443</fpage>–<lpage>53</lpage>. <object-id pub-id-type="pmid">10980420</object-id></mixed-citation></ref>
<ref id="pone.0164823.ref019"><label>19</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Johnson</surname> <given-names>BJ</given-names></name>, <name name-style="western"><surname>Miller</surname> <given-names>GH</given-names></name>, <name name-style="western"><surname>Fogel</surname> <given-names>ML</given-names></name>, <name name-style="western"><surname>Beaumont</surname> <given-names>PB</given-names></name>. <article-title>The determination of late Quaternary paleoenvironments at Equus Cave, South Africa, using stable isotopes and amino acid racemization in ostrich eggshell</article-title>. <source>Palaeogeography, Palaeoclimatology, Palaeoecology</source>. <year>1997</year> <month>Dec</month> <day>15</day>;<volume>136</volume>(<issue>1</issue>):<fpage>121</fpage>–<lpage>37</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref020"><label>20</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Miller</surname> <given-names>GH</given-names></name>, <name name-style="western"><surname>Hart</surname> <given-names>CP</given-names></name>, <name name-style="western"><surname>Roark</surname> <given-names>EB</given-names></name>, <name name-style="western"><surname>Johnson</surname> <given-names>BJ</given-names></name>. <article-title>Isoleucine epimerization in eggshells of the flightless Australian birds Genyomis and Dromaius</article-title>. <source>Perspectives in Amino Acid and Protein Geochemistry</source>. <year>2000</year>:<volume>161</volume>.</mixed-citation></ref>
<ref id="pone.0164823.ref021"><label>21</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Oskam</surname> <given-names>CL</given-names></name>, <name name-style="western"><surname>Haile</surname> <given-names>J</given-names></name>, <name name-style="western"><surname>McLay</surname> <given-names>E</given-names></name>, <name name-style="western"><surname>Rigby</surname> <given-names>P</given-names></name>, <name name-style="western"><surname>Allentoft</surname> <given-names>ME</given-names></name>, <name name-style="western"><surname>Olsen</surname> <given-names>ME</given-names></name>, <etal>et al</etal>. <article-title>Fossil avian eggshell preserves ancient DNA</article-title>. <source>Proceedings of the Royal Society of London B: Biological Sciences</source>. <year>2010</year> <month>Jul</month> <day>7</day>;<volume>277</volume>(<issue>1690</issue>):<fpage>1991</fpage>–<lpage>2000</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref022"><label>22</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Carrino</surname> <given-names>DA</given-names></name>, <name name-style="western"><surname>Dennis</surname> <given-names>JE</given-names></name>, <name name-style="western"><surname>Wu</surname> <given-names>TM</given-names></name>, <name name-style="western"><surname>Arias</surname> <given-names>JL</given-names></name>, <name name-style="western"><surname>Fernandez</surname> <given-names>MS</given-names></name>, <name name-style="western"><surname>Rodriguez</surname> <given-names>JP</given-names></name>, <etal>et al</etal>. <article-title>The avian eggshell extracellular matrix as a model for biomineralization</article-title>. <source>Connective tissue research</source>. <year>1996</year> <month>Jan</month> <day>1</day>;<volume>35</volume>(<issue>1–4</issue>):<fpage>325</fpage>–<lpage>8</lpage>. <object-id pub-id-type="pmid">9084671</object-id></mixed-citation></ref>
<ref id="pone.0164823.ref023"><label>23</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Feng</surname> <given-names>QL</given-names></name>, <name name-style="western"><surname>Zhu</surname> <given-names>X</given-names></name>, <name name-style="western"><surname>Li</surname> <given-names>HD</given-names></name>, <name name-style="western"><surname>Kim</surname> <given-names>TN</given-names></name>. <article-title>Crystal orientation regulation in ostrich eggshells</article-title>. <source>Journal of crystal growth</source>. <year>2001</year> <month>Dec</month> <day>31</day>;<volume>233</volume>(<issue>3</issue>):<fpage>548</fpage>–<lpage>54</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref024"><label>24</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Miller</surname> <given-names>GH</given-names></name>, <name name-style="western"><surname>Beaumont</surname> <given-names>PB</given-names></name>, <name name-style="western"><surname>Jull</surname> <given-names>AJ</given-names></name>, <name name-style="western"><surname>Johnson</surname> <given-names>B</given-names></name>. <article-title>Pleistocene geochronology and palaeothermometry from protein diagenesis in ostrich eggshells: implications for the evolution of modern humans</article-title>. <source>Philosophical Transactions of the Royal Society of London B: Biological Sciences</source>. <year>1992</year> <month>Aug</month> <day>29</day>;<volume>337</volume>(<issue>1280</issue>):<fpage>149</fpage>–<lpage>57</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1098/rstb.1992.0092" xlink:type="simple">10.1098/rstb.1992.0092</ext-link></comment> <object-id pub-id-type="pmid">1357689</object-id></mixed-citation></ref>
<ref id="pone.0164823.ref025"><label>25</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Nys</surname> <given-names>Y</given-names></name>, <name name-style="western"><surname>Gautron</surname> <given-names>J</given-names></name>, <name name-style="western"><surname>Garcia-Ruiz</surname> <given-names>JM</given-names></name>, <name name-style="western"><surname>Hincke</surname> <given-names>MT</given-names></name>. <article-title>Avian eggshell mineralization: biochemical and functional characterization of matrix proteins</article-title>. <source>Comptes Rendus Palevol</source>. <year>2004</year> <month>Oct</month> <day>31</day>;<volume>3</volume>(<issue>6</issue>):<fpage>549</fpage>–<lpage>62</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref026"><label>26</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Wellman-Labadie</surname> <given-names>O</given-names></name>, <name name-style="western"><surname>Lakshminarayanan</surname> <given-names>R</given-names></name>, <name name-style="western"><surname>Hincke</surname> <given-names>MT</given-names></name>. <article-title>Antimicrobial properties of avian eggshell-specific Ctype lectin-like proteins</article-title>. <source>FEBS letters</source>. <year>2008</year> <month>Mar</month> <day>5</day>;<volume>582</volume>(<issue>5</issue>):<fpage>699</fpage>–<lpage>704</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1016/j.febslet.2008.01.043" xlink:type="simple">10.1016/j.febslet.2008.01.043</ext-link></comment> <object-id pub-id-type="pmid">18258195</object-id></mixed-citation></ref>
<ref id="pone.0164823.ref027"><label>27</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Wellman-Labadie</surname> <given-names>O</given-names></name>, <name name-style="western"><surname>Picman</surname> <given-names>J</given-names></name>, <name name-style="western"><surname>Hincke</surname> <given-names>MT</given-names></name>. <article-title>Antimicrobial activity of cuticle and outer eggshell protein extracts from three species of domestic birds</article-title>. <source>British poultry science</source>. <year>2008</year> <month>Mar</month> <day>1</day>;<volume>49</volume> (<issue>2</issue>):<fpage>133</fpage>–<lpage>43</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1080/00071660802001722" xlink:type="simple">10.1080/00071660802001722</ext-link></comment> <object-id pub-id-type="pmid">18409087</object-id></mixed-citation></ref>
<ref id="pone.0164823.ref028"><label>28</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Fink</surname> <given-names>D</given-names></name>, <name name-style="western"><surname>Rojas-Chapana</surname> <given-names>J</given-names></name>, <name name-style="western"><surname>Petrov</surname> <given-names>A</given-names></name>, <name name-style="western"><surname>Tributsch</surname> <given-names>H</given-names></name>, <name name-style="western"><surname>Friedrich</surname> <given-names>D</given-names></name>, <name name-style="western"><surname>Küppers</surname> <given-names>U</given-names></name>, <etal>et al</etal>. <article-title>The “artificial ostrich eggshell” project: Sterilizing polymer foils for food industry and medicine</article-title>. <source>Solar energy materials and solar cells</source>. <year>2006</year> <month>Jun</month> <day>15</day>;<volume>90</volume> (<issue>10</issue>):<fpage>1458</fpage>–<lpage>70</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref029"><label>29</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Kumar</surname> <given-names>G</given-names></name>, <name name-style="western"><surname>Narvare</surname> <given-names>G</given-names></name>, <name name-style="western"><surname>Pancholi</surname> <given-names>RK</given-names></name>. <article-title>Engraved ostrich eggshell objects: new evidence of Upper Palaeolithic art in India</article-title>. <source>Rock Art Research</source>. <year>1988</year>;<volume>5</volume>(<issue>1</issue>):<fpage>43</fpage>–<lpage>53</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref030"><label>30</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Sahni</surname> <given-names>AS</given-names></name>, <name name-style="western"><surname>Kumar</surname> <given-names>GI</given-names></name>, <name name-style="western"><surname>Bajpai</surname> <given-names>SU</given-names></name>, <name name-style="western"><surname>Srinivasan</surname> <given-names>S</given-names></name>. <article-title>Ultrastructure and taxonomy of ostrich eggshells from Upper Palaeolithic sites of India</article-title>. <source>J Palaeontol Soc (India)</source>. <year>1989</year>;<volume>34</volume>:<fpage>91</fpage>–<lpage>8</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref031"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">Kumar G. Chronology: Upper Palaeolithic Culture. Archaeology of northwestern Malwa (Ph.D. Thesis),. Ujjain: Vikram University. 227–229, 1983.</mixed-citation></ref>
<ref id="pone.0164823.ref032"><label>32</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Mohabey</surname> <given-names>DM</given-names></name>. <article-title>Avian Egg Shells from Pleistocene of Kutch</article-title>. <source>Geological Society of India</source>. <year>1989</year> <month>May</month> <day>1</day>;<volume>33</volume>(<issue>5</issue>):<fpage>477</fpage>–<lpage>81</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref033"><label>33</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Bednarik</surname> <given-names>RG</given-names></name>. <article-title>About Palaeolithic ostrich eggshell in India</article-title>. <source>Bulletin of the Indo-Pacific Prehistory Association</source>. <year>1993</year> <month>Jan</month> <day>24</day>;<volume>13</volume>:<fpage>34</fpage>–<lpage>43</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref034"><label>34</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Bednarik</surname> <given-names>RG</given-names></name>. <article-title>Palaeolithic art in India</article-title>. <source>Man and Environment</source>. <year>1993</year>;<volume>18</volume>(<issue>2</issue>):<fpage>33</fpage>–<lpage>40</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref035"><label>35</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Bednarik</surname> <given-names>RG</given-names></name>. <article-title>A taphonomy of palaeoart</article-title>. <source>Antiquity</source>. <year>1994</year> <month>Mar</month> <day>1</day>;<volume>68</volume>(<issue>258</issue>):<fpage>68</fpage>–<lpage>74</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref036"><label>36</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Mellars</surname> <given-names>P</given-names></name>. <article-title>Why did modern human populations disperse from Africa ca. 60,000 years ago? A new model</article-title>. <source>Proceedings of the National Academy of Sciences</source>. <year>2006</year> <month>Jun</month> <day>20</day>;<volume>103</volume>(<issue>25</issue>):<fpage>9381</fpage>–<lpage>6</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref037"><label>37</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Rai</surname> <given-names>N</given-names></name>, <name name-style="western"><surname>Taher</surname> <given-names>N</given-names></name>, <name name-style="western"><surname>Singh</surname> <given-names>M</given-names></name>, <name name-style="western"><surname>Chaubey</surname> <given-names>G</given-names></name>, <name name-style="western"><surname>Jha</surname> <given-names>AN</given-names></name>, <name name-style="western"><surname>Singh</surname> <given-names>L</given-names></name>, <etal>et al</etal>. <article-title>Relic excavated in western India is probably of Georgian Queen Ketevan</article-title>. <source>Mitochondrion</source>. <year>2014</year> <month>Jan</month> <day>31</day>;<volume>14</volume>:<fpage>1</fpage>–<lpage>6</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1016/j.mito.2013.12.002" xlink:type="simple">10.1016/j.mito.2013.12.002</ext-link></comment> <object-id pub-id-type="pmid">24355295</object-id></mixed-citation></ref>
<ref id="pone.0164823.ref038"><label>38</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Poinar</surname> <given-names>HN</given-names></name>, <name name-style="western"><surname>Cooper</surname> <given-names>A</given-names></name>. <article-title>Ancient DNA: do it right or not at all</article-title>. <source>Science</source>. <year>2000</year>;<volume>5482</volume>:<fpage>1139</fpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref039"><label>39</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Haak</surname> <given-names>W</given-names></name>, <name name-style="western"><surname>Balanovsky</surname> <given-names>O</given-names></name>, <name name-style="western"><surname>Sanchez</surname> <given-names>JJ</given-names></name>, <name name-style="western"><surname>Koshel</surname> <given-names>S</given-names></name>, <name name-style="western"><surname>Zaporozhchenko</surname> <given-names>V</given-names></name>, <name name-style="western"><surname>Adler</surname> <given-names>CJ</given-names></name>, <etal>et al</etal>. <article-title>Ancient DNA from European early neolithic farmers reveals their near eastern affinities</article-title>. <source>PLoS Biol</source>. <year>2010</year> <month>Nov</month> <day>9</day>;<volume>8</volume>(<issue>11</issue>):<fpage>e1000536</fpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1371/journal.pbio.1000536" xlink:type="simple">10.1371/journal.pbio.1000536</ext-link></comment> <object-id pub-id-type="pmid">21085689</object-id></mixed-citation></ref>
<ref id="pone.0164823.ref040"><label>40</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Pääbo</surname> <given-names>S</given-names></name>, <name name-style="western"><surname>Poinar</surname> <given-names>H</given-names></name>, <name name-style="western"><surname>Serre</surname> <given-names>D</given-names></name>, <name name-style="western"><surname>Jaenicke-Després</surname> <given-names>V</given-names></name>, <name name-style="western"><surname>Hebler</surname> <given-names>J</given-names></name>, <name name-style="western"><surname>Rohland</surname> <given-names>N</given-names></name>, <etal>et al</etal>. <article-title>Genetic analyses from ancient DNA</article-title>. <source>Annu. Rev. Genet.</source> <year>2004</year> <month>Dec</month> <day>15</day>;<volume>38</volume>:<fpage>645</fpage>–<lpage>79</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1146/annurev.genet.37.110801.143214" xlink:type="simple">10.1146/annurev.genet.37.110801.143214</ext-link></comment> <object-id pub-id-type="pmid">15568989</object-id></mixed-citation></ref>
<ref id="pone.0164823.ref041"><label>41</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Oskam</surname> <given-names>CL</given-names></name>, <name name-style="western"><surname>Bunce</surname> <given-names>M</given-names></name>. <article-title>DNA extraction from fossil eggshell</article-title>. <source>Ancient DNA: Methods and Protocols</source>. <year>2012</year>:<fpage>65</fpage>–<lpage>70</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref042"><label>42</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Saitou</surname> <given-names>N</given-names></name>, <name name-style="western"><surname>Nei</surname> <given-names>M</given-names></name>. <article-title>The neighbor-joining method: a new method for reconstructing phylogenetic trees</article-title>. <source>Molecular biology and evolution</source>. <year>1987</year> <month>Jul</month> <day>1</day>;<volume>4</volume>(<issue>4</issue>):<fpage>406</fpage>–<lpage>25</lpage>. <object-id pub-id-type="pmid">3447015</object-id></mixed-citation></ref>
<ref id="pone.0164823.ref043"><label>43</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Felsenstein</surname> <given-names>J</given-names></name>. <article-title>Confidence limits on phylogenies: an approach using the bootstrap</article-title>. <source>Evolution</source>. <year>1985</year> <month>Jul</month> <volume>1</volume>:<fpage>783</fpage>–<lpage>91</lpage>.</mixed-citation></ref>
<ref id="pone.0164823.ref044"><label>44</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Tamura</surname> <given-names>K</given-names></name>, <name name-style="western"><surname>Nei</surname> <given-names>M</given-names></name>, <name name-style="western"><surname>Kumar</surname> <given-names>S</given-names></name>. <article-title>Prospects for inferring very large phylogenies by using the neighbor-joining method</article-title>. <source>Proceedings of the National Academy of Sciences of the United States of America</source>. <year>2004</year> <month>Jul</month> <day>27</day>;<volume>101</volume>(<issue>30</issue>):<fpage>11030</fpage>–<lpage>5</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1073/pnas.0404206101" xlink:type="simple">10.1073/pnas.0404206101</ext-link></comment> <object-id pub-id-type="pmid">15258291</object-id></mixed-citation></ref>
<ref id="pone.0164823.ref045"><label>45</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Tamura</surname> <given-names>K</given-names></name>, <name name-style="western"><surname>Stecher</surname> <given-names>G</given-names></name>, <name name-style="western"><surname>Peterson</surname> <given-names>D</given-names></name>, <name name-style="western"><surname>Filipski</surname> <given-names>A</given-names></name>, <name name-style="western"><surname>Kumar</surname> <given-names>S</given-names></name>. <article-title>MEGA6: molecular evolutionary genetics analysis version 6.0</article-title>. <source>Molecular biology and evolution</source>. <year>2013</year> <month>Dec</month> <day>1</day>;<volume>30</volume>(<issue>12</issue>):<fpage>2725</fpage>–<lpage>9</lpage>. <comment>doi: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.1093/molbev/mst197" xlink:type="simple">10.1093/molbev/mst197</ext-link></comment> <object-id pub-id-type="pmid">24132122</object-id></mixed-citation></ref>
<ref id="pone.0164823.ref046"><label>46</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Smith</surname> <given-names>CI</given-names></name>, <name name-style="western"><surname>Chamberlain</surname> <given-names>AT</given-names></name>, <name name-style="western"><surname>Riley</surname> <given-names>MS</given-names></name>, <name name-style="western"><surname>Stringer</surname> <given-names>C</given-names></name>, <name name-style="western"><surname>Collins</surname> <given-names>MJ</given-names></name>. <article-title>The thermal history of human fossils and the likelihood of successful DNA amplification</article-title>. <source>Journal of Human Evolution</source>. <year>2003</year> <month>Sep</month> <day>30</day>;<volume>45</volume>(<issue>3</issue>):<fpage>203</fpage>–<lpage>17</lpage>. <object-id pub-id-type="pmid">14580590</object-id></mixed-citation></ref>
</ref-list>
</back>
</article>