Ancient and Modern Genomes Reveal Microsatellites Maintain a Dynamic Equilibrium Through Deep Time.

Microsatellites are widely used in population genetics, but their evolutionary dynamics remain poorly understood. It is unclear whether microsatellite loci drift in length over time. This is important because the mutation processes that underlie these important genetic markers are central to the evolutionary models that employ microsatellites. We identify more than 27 million microsatellites using a novel and unique dataset of modern and ancient Adélie penguin genomes along with data from 63 published chordate genomes. We investigate microsatellite evolutionary dynamics over 2 timescales: one based on Adélie penguin samples dating to ∼46.5 ka and the other dating to the diversification of chordates aged more than 500 Ma. We show that the process of microsatellite allele length evolution is at dynamic equilibrium; while there is length polymorphism among individuals, the length distribution for a given locus remains stable. Many microsatellites persist over very long timescales, particularly in exons and regulatory sequences. These often retain length variability, suggesting that they may play a role in maintaining phenotypic variation within populations.

A Unique Mitochondrial Gene Block Inversion in Antarctic Trematomin Fishes: A Cautionary Tale.

Many Antarctic notothenioid fishes have major rearrangements in their mitochondrial (mt) genomes. Here, we report the complete mt genomes of 3 trematomin notothenioids: the bald notothen (Trematomus (Pagothenia) borchgrevinki), the spotted notothen (T. nicolai), and the emerald notothen (T. bernacchii). The 3 mt genomes were sequenced using next-generation Illumina technology, and the assemblies verified by Sanger sequencing. When compared with the canonical mt gene order of the Antarctic silverfish (Pleuragramma antarctica), we found a large gene inversion in the 3 trematomin mt genomes that included tRNAIle, ND1, tRNALeu2, 16S, tRNAVal, 12S, tRNAPhe, and the control region. The trematomin mt genomes contained 3 intergenic spacers, which are thought to be the remnants of previous gene and control region duplications. All control regions included the characteristic conserved regulatory sequence motifs. Although short-read next-generation DNA sequencing technology has allowed the rapid and cost-effective sequencing of a large number of complete mt genomes, it is essential in all cases to verify the assembly in order to prevent the publication and use of erroneous data.

The design and application of a 50 K SNP chip for a threatened Aotearoa New Zealand passerine, the hihi.

Next-generation sequencing has transformed the fields of ecological and evolutionary genetics by allowing for cost-effective identification of genome-wide variation. Single nucleotide polymorphism (SNP) arrays, or "SNP chips", enable very large numbers of individuals to be consistently genotyped at a selected set of these identified markers, and also offer the advantage of being able to analyse samples of variable DNA quality. We used reduced representation restriction-aided digest sequencing (RAD-seq) of 31 birds of the threatened hihi (Notiomystis cincta; stitchbird) and low-coverage whole genome sequencing (WGS) of 10 of these birds to develop an Affymetrix 50 K SNP chip. We overcame the limitations of having no hihi reference genome and a low quantity of sequence data by separate and pooled de novo assembly of each of the 10 WGS birds. Reads from all individuals were mapped back to these de novo assemblies to identify SNPs. A subset of RAD-seq and WGS SNPs were selected for inclusion on the chip, prioritising SNPs with the highest quality scores whose flanking sequence uniquely aligned to the zebra finch (Taeniopygia guttata) genome. Of the 58,466 SNPs manufactured on the chip, 72% passed filtering metrics and were polymorphic. By genotyping 1,536 hihi on the array, we found that SNPs detected in multiple assemblies were more likely to successfully genotype, representing a cost-effective approach to identify SNPs for genotyping. Here, we demonstrate the utility of the SNP chip by describing the high rates of linkage disequilibrium in the hihi genome, reflecting the history of population bottlenecks in the species.

Sex-specific foraging of an apex predator puts females at risk of human-wildlife conflict.

Urbanisation and anthropogenic alteration of ecosystems has led to conflict between humans and wildlife. Such conflict is often observed in apex predators. Although human-wildlife conflict has been extensively studied, male/female differences in behaviour are rarely considered. We investigated male/female differences in foraging behaviour of the predatory/scavenging brown skua Catharacta antarctica lonnbergi breeding on a New Zealand island nature reserve in proximity to farmland. These skuas are subject to culling, when perceived as a threat to livestock. As part of a long-term ecological study, we used high-resolution Global Positioning System (GPS) devices to characterise the space-use of foraging brown skuas. We also analysed stable isotopes of carbon (δ13 C) and nitrogen (δ15 N) from modern and archived blood samples to investigate possible changes in diet over the past ~30 years. Analysis of 100 GPS tracks collected from 2014 to 2016 demonstrated that males and females consistently visited different habitats. Males spent most of their time close to their breeding territory on the island nature reserve and females frequently visited a farmed island approximately two kilometres away. Consistent with this finding, we show that male and female skuas also differed markedly in their diets: males specialised on burrow-nesting white-faced storm petrels Pelagodroma marina (80%) with only a small proportion of sheep remains Ovis aries (<6%) contributing to their diet. In contrast, female diet comprised 27% white-faced storm petrels, other seabirds (18%) and a relatively large proportion of sheep remains (47%). Further, our data (186 blood samples from 122 individuals) show that this male/female difference in diet has persisted at least since 1987. Because females fed disproportionally on sheep remains, they may be more vulnerable to being culled by farmers. Importantly, our case study suggests that intersexual differences in diet and foraging patterns can have major implications for the reproduction and survival of apex predators that interact with farming. We strongly suggest that intersexual differences in behaviour should be considered when investigating human-wildlife conflicts.

Mitogenomic diversity in Sacred Ibis Mummies sheds light on early Egyptian practices.

The ancient catacombs of Egypt harbor millions of well-preserved mummified Sacred Ibis (Threskiornis aethiopicus) dating from ~600BC. Although it is known that a very large number of these 'votive' mummies were sacrificed to the Egyptian God Thoth, how the ancient Egyptians obtained millions of these birds for mummification remains unresolved. Ancient Egyptian textual evidences suggest they may have been raised in dedicated large-scale farms. To investigate the most likely method used by the priests to secure birds for mummification, we report the first study of complete mitochondrial genomes of 14 Sacred Ibis mummies interred ~2500 years ago. We analysed and compared the mitogenomic diversity among Sacred Ibis mummies to that found in modern Sacred Ibis populations from throughout Africa. The ancient birds show a high level of genetic variation comparable to that identified in modern African populations, contrary to the suggestion in ancient hieroglyphics (or ancient writings) of centralized industrial scale farming of sacrificial birds. This suggests a sustained short-term taming of the wild migratory Sacred Ibis for the ritual yearly demand.

DNA barcoding a unique avifauna: an important tool for evolution, systematics and conservation.

BACKGROUND: DNA barcoding utilises a standardised region of the cytochrome c oxidase I (COI) gene to identify specimens to the species level. It has proven to be an effective tool for identification of avian samples. The unique island avifauna of New Zealand is taxonomically and evolutionarily distinct. We analysed COI sequence data in order to determine if DNA barcoding could accurately identify New Zealand birds. RESULTS: We sequenced 928 specimens from 180 species. Additional Genbank sequences expanded the dataset to 1416 sequences from 211 of the estimated 236 New Zealand species. Furthermore, to improve the assessment of genetic variation in non-endemic species, and to assess the overall accuracy of our approach, sequences from 404 specimens collected outside of New Zealand were also included in our analyses. Of the 191 species represented by multiple sequences, 88.5% could be successfully identified by their DNA barcodes. This is likely a conservative estimate of the power of DNA barcoding in New Zealand, given our extensive geographic sampling. The majority of the 13 groups that could not be distinguished contain recently diverged taxa, indicating incomplete lineage sorting and in some cases hybridisation. In contrast, 16 species showed evidence of distinct intra-species lineages, some of these corresponding to recognised subspecies. For species identification purposes a character-based method was more successful than distance and phylogenetic tree-based methods. CONCLUSIONS: DNA barcodes accurately identify most New Zealand bird species. However, low levels of COI sequence divergence in some recently diverged taxa limit the identification power of DNA barcoding. A small number of currently recognised species would benefit from further systematic investigations. The reference database and analysis presented will provide valuable insights into the evolution, systematics and conservation of New Zealand birds.

Ancient nuclear genomes enable repatriation of Indigenous human remains.

After European colonization, the ancestral remains of Indigenous people were often collected for scientific research or display in museum collections. For many decades, Indigenous people, including Native Americans and Aboriginal Australians, have fought for their return. However, many of these remains have no recorded provenance, making their repatriation very difficult or impossible. To determine whether DNA-based methods could resolve this important problem, we sequenced 10 nuclear genomes and 27 mitogenomes from ancient pre-European Aboriginal Australians (up to 1540 years before the present) of known provenance and compared them to 100 high-coverage contemporary Aboriginal Australian genomes, also of known provenance. We report substantial ancient population structure showing strong genetic affinities between ancient and contemporary Aboriginal Australian individuals from the same geographic location. Our findings demonstrate the feasibility of successfully identifying the origins of unprovenanced ancestral remains using genomic methods.

The Sacred Ibis debate: The first test of evolution.

In 1798, Napoleon Bonaparte's army invaded Egypt, returning with many treasures including large numbers of Sacred Ibis mummies. The ancient Egyptians revered the ibis and mummified literally millions of them. The French naturalist Georges Cuvier used these mummies to challenge an emerging idea of the time, namely Jean-Baptiste Lamarck's theory of evolution. Cuvier detected no measurable differences between mummified Sacred Ibis and contemporary specimens of the same species. Consequently, he argued that this was evidence for the "fixity of species." The "Sacred Ibis debate" predates the so-called "Great Debate" between Cuvier and Geoffroy Saint-Hilaire and the publication of Darwin's On the Origin of Species five decades later. Cuvier's views and his study had a profound influence on the scientific and public perception of evolution, setting back the acceptance of evolutionary theory in Europe for decades.

Genetic structure of the grey side-gilled sea slug (Pleurobranchaea maculata) in coastal waters of New Zealand.

Pleurobranchaea maculata is a rarely studied species of the Heterobranchia found throughout the south and western Pacific-and recently recorded in Argentina-whose population genetic structure is unknown. Interest in the species was sparked in New Zealand following a series of dog deaths caused by ingestions of slugs containing high levels of the neurotoxin tetrodotoxin. Here we describe the genetic structure and demographic history of P. maculata populations from five principle locations in New Zealand based on extensive analyses of 12 microsatellite loci and the COI and CytB regions of mitochondrial DNA (mtDNA). Microsatellite data showed significant differentiation between northern and southern populations with population structure being associated with previously described regional variations in tetrodotoxin concentrations. However, mtDNA sequence data did not support such structure, revealing a star-shaped haplotype network with estimates of expansion time suggesting a population expansion in the Pleistocene era. Inclusion of publicly available mtDNA sequence sea slugs from Argentina did not alter the star-shaped network. We interpret our data as indicative of a single founding population that fragmented following geographical changes that brought about the present day north-south divide in New Zealand waters. Lack of evidence of cryptic species supports data indicating that differences in toxicity of individuals among regions are a consequence of differences in diet.

Genetic Kinship Analyses Reveal That Gray's Beaked Whales Strand in Unrelated Groups.

Some marine mammals are so rarely seen that their life history and social structure remain a mystery. Around New Zealand, Gray's beaked whales (Mesoplodon grayi) are almost never seen alive, yet they are a commonly stranded species. Gray's are unique among the beaked whales in that they frequently strand in groups, providing an opportunity to investigate their social organization. We examined group composition and genetic kinship in 113 Gray's beaked whales with samples collected over a 20-year period. Fifty-six individuals stranded in 19 groups (2 or more individuals), and 57 whales stranded individually. Mitochondrial control region haplotypes and microsatellite genotypes (16 loci) were obtained for 103 whales. We estimated pairwise relatedness between all pairs of individuals and average relatedness within, and between, groups. We identified 6 mother-calf pairs and 2 half-siblings, including 2 whales in different strandings 17 years and 1500 km apart. Surprisingly, none of the adults stranding together were related suggesting that groups are not formed through the retention of kin. These data suggest that both sexes may disperse from their mothers, and groups consisting of unrelated subadults are common. We also found no instances of paternity within the groups. Our results provide the first insights into dispersal, social organization, and the mating system in this rarely sighted species. Why whales strand is still unknown but, in Gray's beaked whales, the dead can tell us much about the living.

Ancient mtDNA sequences from the First Australians revisited.

The publication in 2001 by Adcock et al. [Adcock GJ, et al. (2001) Proc Natl Acad Sci USA 98(2):537-542] in PNAS reported the recovery of short mtDNA sequences from ancient Australians, including the 42,000-y-old Mungo Man [Willandra Lakes Hominid (WLH3)]. This landmark study in human ancient DNA suggested that an early modern human mitochondrial lineage emerged in Asia and that the theory of modern human origins could no longer be considered solely through the lens of the "Out of Africa" model. To evaluate these claims, we used second generation DNA sequencing and capture methods as well as PCR-based and single-primer extension (SPEX) approaches to reexamine the same four Willandra Lakes and Kow Swamp 8 (KS8) remains studied in the work by Adcock et al. Two of the remains sampled contained no identifiable human DNA (WLH15 and WLH55), whereas the Mungo Man (WLH3) sample contained no Aboriginal Australian DNA. KS8 reveals human mitochondrial sequences that differ from the previously inferred sequence. Instead, we recover a total of five modern European contaminants from Mungo Man (WLH3). We show that the remaining sample (WLH4) contains ∼1.4% human DNA, from which we assembled two complete mitochondrial genomes. One of these was a previously unidentified Aboriginal Australian haplotype belonging to haplogroup S2 that we sequenced to a high coverage. The other was a contaminating modern European mitochondrial haplotype. Although none of the sequences that we recovered matched those reported by Adcock et al., except a contaminant, these findings show the feasibility of obtaining important information from ancient Aboriginal Australian remains.

The effects of the general anaesthetic isoflurane on the honey bee (Apis mellifera) circadian clock.

General anaesthesia administered during the day has previously been shown to phase shift the honey bee clock. We describe a phase response curve for honey bees (n=105) to six hour isoflurane anaesthesia. The honey bee isoflurane PRC is "weak" with a delay portion (maximum shift of -1.88 hours, circadian time 0 - 3) but no advance zone. The isoflurane-induced shifts observed here are in direct opposition to those of light. Furthermore, concurrent administration of light and isoflurane abolishes the shifts that occur with isoflurane alone. Light may thus provide a means of reducing isoflurane-induced phase shifts.

High coverage of the complete mitochondrial genome of the rare Gray's beaked whale (Mesoplodon grayi) using Illumina next generation sequencing.

Using an Illumina platform, we shot-gun sequenced the complete mitochondrial genome of Gray's beaked whale (Mesoplodon grayi) to an average coverage of 152X. We performed a de novo assembly using SOAPdenovo2 and determined the total mitogenome length to be 16,347 bp. The nucleotide composition was asymmetric (33.3% A, 24.6% C, 12.6% G, 29.5% T) with an overall GC content of 37.2%. The gene organization was similar to that of other cetaceans with 13 protein-coding genes, 2 rRNAs (12S and 16S), 22 predicted tRNAs and 1 control region or D-loop. We found no evidence of heteroplasmy or nuclear copies of mitochondrial DNA in this individual. Beaked whales within the genus Mesoplodon are rarely seen at sea and their basic biology is poorly understood. These data will contribute to resolving the phylogeography and population ecology of this speciose group.

Distance-dependent patterns of molecular divergences in Tuatara mitogenomes.

Population genetic models predict that populations that are geographically close to each other are expected to be genetically more similar to each other compared to those that are widely separate. However the patterns of relationships between geographic distance and molecular divergences at neutral and constrained regions of the genome are unclear. We attempted to clarify this relationship by sequencing complete mitochondrial genomes of the relic species Tuatara (Sphenodon punctatus) from ten offshore islands of New Zealand. We observed a positive relationship that showed a proportional increase in the neutral diversity at synonymous sites (dS), with increasing geographical distance. In contrast we showed that diversity at evolutionarily constrained sites (dC) was elevated in the case of comparisons involving closely located populations. Conversely diversity was reduced in the case of comparisons between distantly located populations. These patterns were confirmed by a significant negative relationship between the ratio of dC/dS and geographic distance. The observed high dC/dS could be explained by the abundance of deleterious mutations in comparisons involving closely located populations, due to the recent population divergence times. Since distantly related populations were separated over long periods of time, deleterious mutations might have been removed by purifying selection.

Complete mitochondrial genomes of Tuatara endemic to different islands of New Zealand.

Tuatara are the sister taxon to the Squamata (including lizards and snakes) and are regarded as the most distinctive surviving reptilian genus. They are currently inhabits on offshore islands around New Zealand and have been recognized as a species in need of active conservation management. In this study, we report a total number of five nearly complete mitochondrial genomes, which were sequenced by Sanger and Next Generation DNA sequencing methods. Our phylogenomic analysis revealed distinct clustering of tuatara populations from the north and south islands of New Zealand.

Reconstruction and in vivo analysis of the extinct tbx5 gene from ancient wingless moa (Aves: Dinornithiformes).

BACKGROUND: The forelimb-specific gene tbx5 is highly conserved and essential for the development of forelimbs in zebrafish, mice, and humans. Amongst birds, a single order, Dinornithiformes, comprising the extinct wingless moa of New Zealand, are unique in having no skeletal evidence of forelimb-like structures. RESULTS: To determine the sequence of tbx5 in moa, we used a range of PCR-based techniques on ancient DNA to retrieve all nine tbx5 exons and splice sites from the giant moa, Dinornis. Moa Tbx5 is identical to chicken Tbx5 in being able to activate the downstream promotors of fgf10 and ANF. In addition we show that missexpression of moa tbx5 in the hindlimb of chicken embryos results in the formation of forelimb features, suggesting that Tbx5 was fully functional in wingless moa. An alternatively spliced exon 1 for tbx5 that is expressed specifically in the forelimb region was shown to be almost identical between moa and ostrich, suggesting that, as well as being fully functional, tbx5 is likely to have been expressed normally in moa since divergence from their flighted ancestors, approximately 60 mya. CONCLUSIONS: The results suggests that, as in mice, moa tbx5 is necessary for the induction of forelimbs, but is not sufficient for their outgrowth. Moa Tbx5 may have played an important role in the development of moa's remnant forelimb girdle, and may be required for the formation of this structure. Our results further show that genetic changes affecting genes other than tbx5 must be responsible for the complete loss of forelimbs in moa.

Way-finding in displaced clock-shifted bees proves bees use a cognitive map.

Mammals navigate by means of a metric cognitive map. Insects, most notably bees and ants, are also impressive navigators. The question whether they, too, have a metric cognitive map is important to cognitive science and neuroscience. Experimentally captured and displaced bees often depart from the release site in the compass direction they were bent on before their capture, even though this no longer heads them toward their goal. When they discover their error, however, the bees set off more or less directly toward their goal. This ability to orient toward a goal from an arbitrary point in the familiar environment is evidence that they have an integrated metric map of the experienced environment. We report a test of an alternative hypothesis, which is that all the bees have in memory is a collection of snapshots that enable them to recognize different landmarks and, associated with each such snapshot, a sun-compass-referenced home vector derived from dead reckoning done before and after previous visits to the landmark. We show that a large shift in the sun-compass rapidly induced by general anesthesia does not alter the accuracy or speed of the homeward-oriented flight made after the bees discover the error in their initial postrelease flight. This result rules out the sun-referenced home-vector hypothesis, further strengthening the now extensive evidence for a metric cognitive map in bees.