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 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.

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.

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 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.

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.

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.

General anesthesia alters time perception by phase shifting the circadian clock.

Following general anesthesia, people are often confused about the time of day and experience sleep disruption and fatigue. It has been hypothesized that these symptoms may be caused by general anesthesia affecting the circadian clock. The circadian clock is fundamental to our well-being because it regulates almost all aspects of our daily biochemistry, physiology, and behavior. Here, we investigated the effects of the most common general anesthetic, isoflurane, on time perception and the circadian clock using the honeybee (Apis mellifera) as a model. A 6-h daytime anesthetic systematically altered the time-compensated sun compass orientation of the bees, with a mean anticlockwise shift in vanishing bearing of 87° in the Southern Hemisphere and a clockwise shift in flight direction of 58° in the Northern Hemisphere. Using the same 6-h anesthetic treatment, time-trained bees showed a delay in the start of foraging of 3.3 h, and whole-hive locomotor-activity rhythms were delayed by an average of 4.3 h. We show that these effects are all attributable to a phase delay in the core molecular clockwork. mRNA oscillations of the central clock genes cryptochrome-m and period were delayed by 4.9 and 4.3 h, respectively. However, this effect is dependent on the time of day of administration, as is common for clock effects, and nighttime anesthesia did not shift the clock. Taken together, our results suggest that general anesthesia during the day causes a persistent and marked shift of the clock effectively inducing "jet lag" and causing impaired time perception. Managing this effect in humans is likely to help expedite postoperative recovery.

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.

Resurrecting ancient animal genomes: the extinct moa and more.

Recently two developments have had a major impact on the field of ancient DNA (aDNA). First, new advances in DNA sequencing, in combination with improved capture/enrichment methods, have resulted in the recovery of orders of magnitude more DNA sequence data from ancient animals. Second, there has been an increase in the range of tissue types employed in aDNA. Hair in particular has proven to be very successful as a source of DNA because of its low levels of contamination and high level of ancient endogenous DNA. These developments have resulted in significant advances in our understanding of recently extinct animals: namely their evolutionary relationships, physiology, and even behaviour. Hair has been used to recover the first complete ancient nuclear genome, that of the extinct woolly mammoth, which then facilitated the expression and functional analysis of haemoglobins. Finally, we speculate on the consequences of these developments for the possibility of recreating extinct animals.

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.

Evidence for a recent origin of penguins.

Penguins are a remarkable group of birds, with the 18 extant species living in diverse climatic zones from the tropics to Antarctica. The timing of the origin of these extant penguins remains controversial. Previous studies based on DNA sequences and fossil records have suggested widely differing times for the origin of the group. This has given rise to widely differing biogeographic narratives about their evolution. To resolve this problem, we sequenced five introns from 11 species representing all genera of living penguins. Using these data and other available DNA sequences, together with the ages of multiple penguin fossils to calibrate the molecular clock, we estimated the age of the most recent common ancestor of extant penguins to be 20.4 Myr (17.0-23.8 Myr). This time is half of the previous estimates based on molecular sequence data. Our results suggest that most of the major groups of extant penguins diverged 11-16 Ma. This overlaps with the sharp decline in Antarctic temperatures that began approximately 12 Ma, suggesting a possible relationship between climate change and penguin evolution.

Ancient DNA reveals extreme egg morphology and nesting behavior in New Zealand's extinct moa.

New Zealand's extinct flightless moa radiated rapidly into a large number of morphologically diverse species, which produced an equally large range of egg morphologies. The exact number of moa species, as well as the characteristics of the eggs they laid, remains contentious. Moreover, like most extinct species, we understand little about their nesting and incubation habits. We used a modified ancient DNA extraction procedure to recover exogenous mitochondrial and nuclear DNA from the inside and outside surfaces of moa eggs. We used sequences from the inside of 69 eggshells to directly assign these remains to seven of the 10 currently recognized moa species. In addition we were able to assign, to the species level, six of the rare reconstructed "whole" eggs. These molecular results enabled us to identify two distinct lineages within the genus Euryapteryx. Members of these lineages differed in eggshell thickness, with one lineage being characterized by a relatively thin eggshell. Unexpectedly, several thin-shelled eggs were also shown to belong to the heaviest moa of the genera Dinornis, Euryapteryx and Emeus, making these, to our knowledge, the most fragile of all avian eggs measured to date. Moreover, sex-specific DNA recovered from the outer surfaces of eggshells belonging to species of Dinornis and Euryapteryx suggest that these very thin eggs were likely to have been incubated by the lighter males. The thin nature of the eggshells of these larger species of moa, even if incubated by the male, suggests that egg breakage in these species would have been common if the typical contact method of avian egg incubation was used.

High mitogenomic evolutionary rates and time dependency.

Using entire modern and ancient mitochondrial genomes of Adélie penguins (Pygoscelis adeliae) that are up to 44000 years old, we show that the rates of evolution of the mitochondrial genome are two to six times greater than those estimated from phylogenetic comparisons. Although the rate of evolution at constrained sites, including nonsynonymous positions and RNAs, varies more than twofold with time (between shallow and deep nodes), the rate of evolution at synonymous sites remains the same. The time-independent neutral evolutionary rates reported here would be useful for the study of recent evolutionary events.

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.

Rapid molecular evolution in a living fossil.

The tuatara of New Zealand is a unique reptile that coexisted with dinosaurs and has changed little morphologically from its Cretaceous relatives. Tuatara have very slow metabolic and growth rates, long generation times and slow rates of reproduction. This suggests that the species is likely to exhibit a very slow rate of molecular evolution. Our analysis of ancient and modern tuatara DNA shows that, surprisingly, tuatara have the highest rate of molecular change recorded in vertebrates. Our work also suggests that rates of neutral molecular and phenotypic evolution are decoupled.

Mutation and evolutionary rates in adélie penguins from the antarctic.

Precise estimations of molecular rates are fundamental to our understanding of the processes of evolution. In principle, mutation and evolutionary rates for neutral regions of the same species are expected to be equal. However, a number of recent studies have shown that mutation rates estimated from pedigree material are much faster than evolutionary rates measured over longer time periods. To resolve this apparent contradiction, we have examined the hypervariable region (HVR I) of the mitochondrial genome using families of Adélie penguins (Pygoscelis adeliae) from the Antarctic. We sequenced 344 bps of the HVR I from penguins comprising 508 families with 915 chicks, together with both their parents. All of the 62 germline heteroplasmies that we detected in mothers were also detected in their offspring, consistent with maternal inheritance. These data give an estimated mutation rate (micro) of 0.55 mutations/site/Myrs (HPD 95% confidence interval of 0.29-0.88 mutations/site/Myrs) after accounting for the persistence of these heteroplasmies and the sensitivity of current detection methods. In comparison, the rate of evolution (k) of the same HVR I region, determined using DNA sequences from 162 known age sub-fossil bones spanning a 37,000-year period, was 0.86 substitutions/site/Myrs (HPD 95% confidence interval of 0.53 and 1.17). Importantly, the latter rate is not statistically different from our estimate of the mutation rate. These results are in contrast to the view that molecular rates are time dependent.