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 tRNA^Ile, ND1, tRNA^Leu2, 16S, tRNA^Val, 12S, tRNA^Phe, 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 (δ¹³ C) and nitrogen (δ¹⁵ 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.

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.

Bucking the trend: genetic analysis reveals high diversity, large population size and low differentiation in a deep ocean cetacean

Understanding the genetic structure of a population is essential to its conservation and management. We report the level of genetic diversity and determine the population structure of a cryptic deep ocean cetacean, the Gray's beaked whale (Mesoplodon grayi). We analysed 530 bp of mitochondrial control region and 12 microsatellite loci from 94 individuals stranded around New Zealand and Australia. The samples cover a large area of the species distribution (~6000 km) and were collected over a 22-year period. We show high genetic diversity (h=0.933-0.987, π=0.763-0.996% and Rs=4.22-4.37, He=0.624-0.675), and, in contrast to other cetaceans, we found a complete lack of genetic structure in both maternally and biparentally inherited markers. The oceanic habitats around New Zealand are diverse with extremely deep waters, seamounts and submarine canyons that are suitable for Gray's beaked whales and their prey. We propose that the abundance of this rich habitat has promoted genetic homogeneity in this species. Furthermore, it has been suggested that the lack of beaked whale sightings is the result of their low abundance, but this is in contrast to our estimates of female effective population size based on mitochondrial data. In conclusion, the high diversity and lack of genetic structure can be explained by a historically large population size, in combination with no known exploitation, few apparent behavioural barriers and abundant habitat.

Ancient population genomics and the study of evolution

Recently, the study of ancient DNA (aDNA) has been greatly enhanced by the development of second-generation DNA sequencing technologies and targeted enrichment strategies. These developments have allowed the recovery of several complete ancient genomes, a result that would have been considered virtually impossible only a decade ago. Prior to these developments, aDNA research was largely focused on the recovery of short DNA sequences and their use in the study of phylogenetic relationships, molecular rates, species identification and population structure. However, it is now possible to sequence a large number of modern and ancient complete genomes from a single species and thereby study the genomic patterns of evolutionary change over time. Such a study would herald the beginnings of ancient population genomics and its use in the study of evolution. Species that are amenable to such large-scale studies warrant increased research effort. We report here progress on a population genomic study of the Adélie penguin (Pygoscelis adeliae). This species is ideally suited to ancient population genomic research because both modern and ancient samples are abundant in the permafrost conditions of Antarctica. This species will enable us to directly address many of the fundamental questions in ecology and evolution.

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.

Two Antarctic penguin genomes reveal insights into their evolutionary history and molecular changes related to the Antarctic environment

BACKGROUND

Penguins are flightless aquatic birds widely distributed in the Southern Hemisphere. The distinctive morphological and physiological features of penguins allow them to live an aquatic life, and some of them have successfully adapted to the hostile environments in Antarctica. To study the phylogenetic and population history of penguins and the molecular basis of their adaptations to Antarctica, we sequenced the genomes of the two Antarctic dwelling penguin species, the Adélie penguin [Pygoscelis adeliae] and emperor penguin [Aptenodytes forsteri].

RESULTS

Phylogenetic dating suggests that early penguins arose ~60 million years ago, coinciding with a period of global warming. Analysis of effective population sizes reveals that the two penguin species experienced population expansions from ~1 million years ago to ~100 thousand years ago, but responded differently to the climatic cooling of the last glacial period. Comparative genomic analyses with other available avian genomes identified molecular changes in genes related to epidermal structure, phototransduction, lipid metabolism, and forelimb morphology.

CONCLUSIONS

Our sequencing and initial analyses of the first two penguin genomes provide insights into the timing of penguin origin, fluctuations in effective population sizes of the two penguin species over the past 10 million years, and the potential associations between these biological patterns and global climate change. The molecular changes compared with other avian genomes reflect both shared and diverse adaptations of the two penguin species to the Antarctic environment.

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.

DNA fingerprinting in zoology: past, present, future

In 1962, Thomas Kuhn famously argued that the progress of scientific knowledge results from periodic 'paradigm shifts' during a period of crisis in which new ideas dramatically change the status quo. Although this is generally true, Alec Jeffreys' identification of hypervariable repeat motifs in the human beta-globin gene, and the subsequent development of a technology known now as 'DNA fingerprinting', also resulted in a dramatic shift in the life sciences, particularly in ecology, evolutionary biology, and forensics. The variation Jeffreys recognized has been used to identify individuals from tissue samples of not just humans, but also of many animal species. In addition, the technology has been used to determine the sex of individuals, as well as paternity/maternity and close kinship. We review a broad range of such studies involving a wide diversity of animal species. For individual researchers, Jeffreys' invention resulted in many ecologists and evolutionary biologists being given the opportunity to develop skills in molecular biology to augment their whole organism focus. Few developments in science, even among the subsequent genome discoveries of the 21st century, have the same wide-reaching significance. Even the later development of PCR-based genotyping of individuals using microsatellite repeats sequences, and their use in determining multiple paternity, is conceptually rooted in Alec Jeffreys' pioneering work.

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.

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.

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.

Adélie penguins and temperature changes in Antarctica: a long-term view

During the summer months, Adélie penguins represent the dominant biomass of terrestrial Antarctica. Literally millions of individuals nest in ice-free areas around the coast of the continent. Hence, these modern populations of Adélie penguins have often been championed as an ideal biological indicator of ecological and environmental changes that we currently face. In addition, Adélie penguins show an extraordinary record of sub-fossil remains, dating back to the late Pleistocene. At this time, temperatures were much lower than now. Hence, this species offers unique long-term information, at both the genomic and ecological levels, about how a species has responded to climate change over more than 40 000 years.

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.

A honey bee (Apis mellifera) light phase response curve

The authors report a phase response curve (PRC) for individual honey bees (Apis mellifera) to single 1-h light pulses (1000 lux) using an Aschoff Type 1 protocol (n = 134). The bee PRC is a weak (Type 1) PRC with a maximum advance of 1.5 h between circadian time (CT) 18 and 3 and a maximum delay of 1.5 h between CT 12 and 18. This is the first published honey bee light PRC and provides an important resource for chronobiologists and honey bee researchers. It may also have practical applications for what is an economically important species frequently transported across different time zones.

Birdstrikes and barcoding: can DNA methods help make the airways safer?

While flying remains one of the safest means of travel, reported birdstrikes on aircraft have risen. This is a result of increased aircraft flight movements, changes in agricultural methods and greater environmental awareness contributing to growing populations of hazardous bird species, as well as more diligent reporting of incidents. Measures to mitigate this hazard require accurate data about the species involved; however, the remains of birds from these incidents are often not easy to identify. Reported birdstrikes include a substantial number where the species cannot be determined from morphology alone. DNA barcoding offers a reliable method of identifying species from very small amounts of organic material such as blood, muscle and feathers. We compare species identification based on morphological criteria and identifications based on mitochondrial cytochrome c oxidase subunit I DNA barcoding methods for New Zealand species. Our data suggest that DNA-based identification can substantially add to the accuracy of species identifications, and these methods represent an important addition to existing procedures to improve air safety. In addition, we outline simple and effective protocols for the recovery and processing of samples for DNA barcoding.

Ancient DNA recovers the origins of Māori feather cloaks

Feather cloaks ("kakahu"), particularly those adorned with kiwi feathers, are treasured items or "taonga" to the Māori people of "Aotearoa"/New Zealand. They are considered iconic expression of Māori culture. Despite their status, much of our knowledge of the materials used to construct cloaks, the provenance of cloaks, and the origins of cloak making itself, has been lost. We used ancient DNA methods to recover mitochondrial DNA sequences from 849 feather samples taken from 109 cloaks. We show that almost all (>99%) of the cloaks were constructed using feathers from North Island brown kiwi. Molecular sexing of nuclear DNA recovered from 92 feather cloak samples also revealed that the sex ratio of birds deviated from a ratio of 1:1 observed in reference populations. Additionally, we constructed a database of 185 mitochondrial control region DNA sequences of kiwi feathers comprising samples collected from 26 North Island locations together with data available from the literature. Genetic subdivision (G(ST)), nucleotide subdivision (N(ST)) and Spatial Analysis of Molecular Variants (SAMOVA) analyses revealed high levels of genetic structuring in North Island brown kiwi. Together with sequence data from previously studied ancient and modern kiwi samples, we were able to determine the geographic provenance of 847 cloak feathers from 108 cloaks. A surprising proportion (15%) of cloaks were found to contain feathers from different geographic locations, providing evidence of kiwi trading among Māori tribes or organized hunting trips into other tribal areas. Our data also suggest that the east of the North Island of New Zealand was the most prolific of all kiwi cloak making areas, with over 50% of all cloaks analyzed originating from this region. Similar molecular approaches have the potential to discover a wealth of lost information from artifacts of endemic cultures worldwide.

Next generation sequencing and analysis of a conserved transcriptome of New Zealand's kiwi

BACKGROUND

Kiwi is a highly distinctive, flightless and endangered ratite bird endemic to New Zealand. To understand the patterns of molecular evolution of the nuclear protein-coding genes in brown kiwi (Apteryx australis mantelli) and to determine the timescale of avian history we sequenced a transcriptome obtained from a kiwi embryo using next generation sequencing methods. We then assembled the conserved protein-coding regions using the chicken proteome as a scaffold.

RESULTS

Using 1,543 conserved protein coding genes we estimated the neutral evolutionary divergence between the kiwi and chicken to be ~45%, which is approximately equal to the divergence computed for the human-mouse pair using the same set of genes. A large fraction of genes was found to be under high selective constraint, as most of the expressed genes appeared to be involved in developmental gene regulation. Our study suggests a significant relationship between gene expression levels and protein evolution. Using sequences from over 700 nuclear genes we estimated the divergence between the two basal avian groups, Palaeognathae and Neognathae to be 132 million years, which is consistent with previous studies using mitochondrial genes.

CONCLUSIONS

The results of this investigation revealed patterns of mutation and purifying selection in conserved protein coding regions in birds. Furthermore this study suggests a relatively cost-effective way of obtaining a glimpse into the fundamental molecular evolutionary attributes of a genome, particularly when no closely related genomic sequence is available.

Size dimorphism and avian-perceived sexual dichromatism in a New Zealand endemic bird, the whitehead Mohoua albicilla

Sex differences in behavior, morphology, and physiology are common in animals. In many bird species, differences in the feather colors of the sexes are apparent when judged by human observers and using physical measures of plumage reflectance, cryptic (to human) plumage dichromatism has also been detected in several additional avian lineages. However, it remains to be confirmed in almost all species whether sexual dichromatism is perceivable by individuals of the studied species. This latter step is essential because it allows the evaluation of alternative hypotheses regarding the signaling and communication functions of plumage variation. We applied perceptual modeling of the avian visual system for the first time to an endemic New Zealand bird to provide evidence of subtle but consistent sexual dichromatism in the whitehead, Mohoua albicilla. Molecular sexing techniques were also used in this species to confirm the extent of the sexual size dimorphism in plumage and body mass. Despite the small sample sizes, we now validate previous reports based on human perception that in male whiteheads head and chest feathers are physically brighter than in females. We further suggest that the extent of sexual plumage dichromatism is pronounced and can be perceived by these birds. In contrast, although sexual dimorphism was also detectable in the mass among the DNA-sexed individuals, it was found to be less extensive than previously thought. Sexual size dimorphism and intraspecifically perceivable plumage dichromatism represent reliable traits that differ between female and male whiteheads. These traits, in turn, may contribute to honest communication displays within the complex social recognition systems of communally breeding whitehead and other group-breeding taxa.

Why does insect RNA look degraded?

The integrity of extracted ribonucleic acid (RNA) is commonly assessed by gel electrophoresis and subsequent analysis of the ribosomal RNA (rRNA) bands. Using the honey bee, Apis mellifera (Hymenoptera: Apidae), as an example, the electrophoretic rRNA profile of insects is explained. This profile differs significantly from the standard benchmark since the 28S rRNA of most insects contains an endogenous "hidden break." Upon denaturation, the masking hydrogen bonds are disrupted, releasing two similar sized fragments that both migrate closely with 18S rRNA. The resulting rRNA profile thus reflects the endogenous composition of insect rRNA and should not be misinterpreted as degradation.

The molecular ecology of the extinct New Zealand Huia

The extinct Huia (Heteralocha acutirostris) of New Zealand represents the most extreme example of beak dimorphism known in birds. We used a combination of nuclear genotyping methods, molecular sexing, and morphometric analyses of museum specimens collected in the late 19(th) and early 20(th) centuries to quantify the sexual dimorphism and population structure of this extraordinary species. We report that the classical description of Huia as having distinctive sex-linked morphologies is not universally correct. Four Huia, sexed as females had short beaks and, on this basis, were indistinguishable from males. Hence, we suggest it is likely that Huia males and females were indistinguishable as juveniles and that the well-known beak dimorphism is the result of differential beak growth rates in males and females. Furthermore, we tested the prediction that the social organisation and limited powers of flight of Huia resulted in high levels of population genetic structure. Using a suite of microsatellite DNA loci, we report high levels of genetic diversity in Huia, and we detected no significant population genetic structure. In addition, using mitochondrial hypervariable region sequences, and likely mutation rates and generation times, we estimated that the census population size of Huia was moderately high. We conclude that the social organization and limited powers of flight did not result in a highly structured population.

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.

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.

New genetic approach to detecting individuals of rare and endangered species

Many rare and endangered species are difficult to locate, observe, and study. Consequently, many individuals, breeding pairs, and even populations of such species could remain undetected. Genetic markers can potentially be used to detect the existence of undiscovered individuals and populations, and we propose a method to do so that requires 3 conditions. First, sampling of the known population(s) of the target species must be comprehensive. Second, the species must display a reasonable level of philopatry and genetic structuring. Third, individuals must be able to be caught outside of breeding locations (e.g., at courtship or feeding areas, in flight), and the level of recapture must be reasonably high. We applied our method to the Chatham Island Taiko (Pterodroma magentae), one of the world's most endangered seabirds. We sequenced the Taiko mitochondrial cytochrome b gene and both copies of a fragment of the duplicated domain I of the control region. Twenty-one haplotypes were revealed, including 4 (19%) not found in birds at known burrows. These results suggest there are more burrow groups yet to be located. The species is a pelagic gadfly petrel that inhabits land only in the breeding season during which it is nocturnal and nests in burrows. Taiko burrows are situated in dense forest in a remote area of Chatham Island, and are consequently difficult to locate and study. It is important that all Taiko burrows be discovered to enable monitoring and protection of the birds from exotic predators.

New developments in ancient genomics

Ancient DNA research is on the crest of a 'third wave' of progress due to the introduction of a new generation of DNA sequencing technologies. Here we review the advantages and disadvantages of the four new DNA sequencers that are becoming available to researchers. These machines now allow the recovery of orders of magnitude more DNA sequence data, albeit as short sequence reads. Hence, the potential reassembly of complete ancient genomes seems imminent, and when used to screen libraries of ancient sequences, these methods are cost effective. This new wealth of data is also likely to herald investigations into the functional properties of extinct genes and gene complexes and will improve our understanding of the biological basis of extinct phenotypes.

Evaluating molecular and behavioural sexing methods for the Australasian gannet (Morus serrator)

The availability of molecular methods for avian sex identification has revolutionised the study of sexual differences in behaviour, morphology, life-history traits and conservation management. We implemented the recommendations of a recent review of DNA-based sex-identification by (1) verifying the sex-specificity and (2) estimating the accuracy of different sex-assignment methods in an apparently monomorphic seabird, the Australasian gannet (Morus serrator). The polymerase chain reaction (PCR) method based on the amplification of the sex-linked chromodomain-helicase-DNA binding gene (CHD) repeatedly assigned the same sex in 96% (n = 27 replicates) and correctly sexed all individuals with known gonadal anatomy (n = 6). PCR and sex-specific restriction fragment length polymorphism (RFLPs) showed agreement for 99.5% of individuals (n = 201). DNA-sexed pairs known to be social mates consisted of a male and a female in 96% of pairs sexed by PCR (n = 77) and 98% of pairs sexed by RFLP (n = 65). DNA-sexed females were in the bottom and males in the top copulatory position in 86% of observed copulations (n = 43 individuals). These results validate assumptions that both membership in social pairs and different copulatory positions can serve as reliable behavioural proxies for field-based sex identification in this colonial and obligately biparental seabird.

Microevolution and mega-icebergs in the Antarctic

Microevolution is regarded as changes in the frequencies of genes in populations over time. Ancient DNA technology now provides an opportunity to demonstrate evolution over a geological time frame and to possibly identify the causal factors in any such evolutionary event. Using nine nuclear microsatellite DNA loci, we genotyped an ancient population of Adélie penguins (Pygoscelis adeliae) aged approximately 6,000 years B.P. Subfossil bones from this population were excavated by using an accurate stratigraphic method that allowed the identification of individuals even within the same layer. We compared the allele frequencies in the ancient population with those recorded from the modern population at the same site in Antarctica. We report significant changes in the frequencies of alleles between these two time points, hence demonstrating microevolutionary change. This study demonstrates a nuclear gene-frequency change over such a geological time frame. We discuss the possible causes of such a change, including the role of mutation, genetic drift, and the effects of gene mixing among different penguin populations. The latter is likely to be precipitated by mega-icebergs that act to promote migration among penguin colonies that typically show strong natal return.

Reconstructing the tempo and mode of evolution in an extinct clade of birds with ancient DNA: the giant moas of New Zealand

The tempo and mode of evolution of the extinct giant moas of New Zealand remain obscure because the number of lineages and their divergence times cannot be estimated reliably by using fossil bone characters only. We therefore extracted ancient DNA from 125 specimens and genetically typed them for a 658-bp mtDNA control region sequence. The sequences detected 14 monophyletic lineages, 9 of which correspond to currently recognized species. One of the newly detected lineages was a genetically divergent form of Megalapteryx originally described as a separate species, two more were lineages of Pachyornis in southern and northeastern New Zealand, and two were basal lineages of South Island Dinornis. When results from genetic typing and previous molecular sexing were combined, at least 33.6% of the specimens were incorrectly classified. We used longer sequences of the control region and nine other mtDNA genes totaling 2,814 base pairs to derive a strongly supported phylogeny of the 14 moa lineages. Molecular dating estimated the most recent common ancestor of moas existed after the Oligocene drowning of New Zealand. However, a cycle of lineage-splitting occurred approximately 4-10 million years ago, when the landmass was fragmented by tectonic and mountain-building events and general cooling of the climate. These events resulted in the geographic isolation of lineages and ecological specialization. The spectacular radiation of moa lineages involved significant changes in body size, shape, and mass and provides another example of the general influence of large-scale paleoenvironmental changes on vertebrate evolutionary history.

Is a large-scale DNA-based inventory of ancient life possible?

A complete DNA-based inventory of the Earth's present biota using large-scale high-throughput DNA sequencing of signature region(s) (DNA barcoding) is an ambitious proposal rivaling the Human Genome Project. We examine whether this approach will also enable us to assess the past diversity of the earth's biota. To test this, we sequenced the 5' terminus of the mitochondrial cytochrome c oxidase I (COI) gene of individuals belonging to a group of extinct ratite birds, the moa of New Zealand. Moa comprised a large number of taxa that radiated in isolation on this oceanic landmass. Using a phylogenetic approach based on a large data set including protein coding and 12S DNA sequences as well as morphology, we now have precise information about the number of moa species that once existed. We show that each of the moa species detected using this extensive data set has a unique COI barcode(s) and that they all show low levels of within-species COI variation. Consequently, we conclude that COI sequences accurately identify the species discovered using the larger data set. Hence, more generally, this study suggests that DNA barcoding might also help us detect other extinct animal species and that a large-scale inventory of ancient life is possible.

Ancient DNA enables timing of the pleistocene origin and holocene expansion of two adélie penguin lineages in antarctica

The timing of divergent events in history is one of the central goals of contemporary evolutionary biology. Such studies are however dependent on accurate evolutionary rates. Recent developments in ancient DNA analysis enable the estimation of more accurate evolutionary rates and therefore more accurate timing of divergence events. Consequently, this leads to a better understanding of changes in populations through time. We use an evolutionary rate calculated from ancient DNA of Adélie penguins (Pygoscelis adeliae) to time divergent events in their history. We report the presence of two distinct and highly variable mitochondrial DNA lineages and track changes in these lineages through space and time. When the ancient DNA and the phylogenetic rates are used to estimate the time of origin of the lineages, two very different estimates resulted. In addition, these same rates provide very different estimates of the time of expansion of these lineages. We suggest that the rate calculated from ancient DNA is more consistent with the glacial history of Antarctica and requires fewer assumptions than does a narrative based on the phylogenetic rate. Finally, we suggest that our study indicates an important new role for ancient DNA studies in the timing of divergent events in history.

Nuclear DNA sequences detect species limits in ancient moa

Ancient DNA studies have typically used multi-copy mitochondrial DNA sequences. This is largely because single-locus nuclear genes have been difficult to recover from sub-fossil material, restricting the scope of ancient DNA research. Here, we have isolated single-locus nuclear DNA markers to assign the sex of 115 extinct moa and, in combination with a mitochondrial DNA phylogeny, tested competing hypotheses about the specific status of moa taxa. Moa were large ratite birds that showed extreme size variation both within and among species. For some taxa, this large variation was hypothesized to represent sexual dimorphism, while for others it was argued to reflect the existence of different species. Our results show that moa were characterized by extreme reverse sexual dimorphism and as a result we have been able to clarify the number of moa species. For example, we show that the three recognized 'species' of Dinornis comprised only two monophyletic groups and that two of these 'species' comprised individuals of one sex only. This study also illustrates that single-locus nuclear DNA sequences can be consistently recovered from ancient material.

A DNA test to sex ratite birds

DNA-based sex tests now exist for many avian species. However, none of these tests are widely applicable to ratites. We present DNA sequence data for a locus that is W chromosome-linked in the kiwi, ostrich, cassowary, rhea, and emu. At the amino acid level, this sequence has significant homology to X-linked genes in platyfish and Caenorhabditis elegans. Polymerase chain reaction (PCR) primers designed to this locus allow the assignment of sex in all species of living ratites.

Rates of evolution in ancient DNA from Adélie penguins

Well-preserved subfossil bones of Adélie penguins, Pygoscelis adeliae, underlie existing and abandoned nesting colonies in Antarctica. These bones, dating back to more than 7000 years before the present, harbor some of the best-preserved ancient DNA yet discovered. From 96 radiocarbon-aged bones, we report large numbers of mitochondrial haplotypes, some of which appear to be extinct, given the 380 living birds sampled. We demonstrate DNA sequence evolution through time and estimate the rate of evolution of the hypervariable region I using a Markov chain Monte Carlo integration and a least-squares regression analysis. Our calculated rates of evolution are approximately two to seven times higher than previous indirect phylogenetic estimates.

Molecular sexing of the communally breeding pukeko: an important ecological tool

A central biological parameter in the study of any animal population is the accurate assignment of sex. Indeed any ecological study of a population requires information on sex composition in relation to such biological factors as behaviour, movement, mortality and birth rate. However, our ability to assign the sex of adults of many avian species is poor and the sexing of young is universally difficult. We report here the successful application of a molecular technique for the assignment of sex in the communally breeding pukeko or purple swamphen (Porphyrio porphyrio melanotus). W- and Z-linked chromosome fragments in digested genomic DNA of pukeko were detected with the DNA probe pMg1. We consequently show that this species breeds in polyandrous, polygynous and polygynandrous groups. Finally we discuss why recent molecular methods represent important new tools in ecology.

DNA science and conservation

A wide array of DNA-based genetic techniques are now available for the study of many problems in conservation biology. Either directly or indirectly, such techniques are becoming increasingly available to scientists and managers alike. Although these technologies are generally known to conservationists, there is a need to clearly outline the principal characteristics of such genetic tools and to detail how they can most appropriately be used in the management of wildlife species. The essential characteristics of mitochondrial and chloroplast restriction fragment analyses are detailed, together with discussions of single locus nuclear restriction fragment length polymorphisms (RFLPs), multilocus DNA fingerprinting, microsatellite DNA, randomly amplified polymorphic DNAs (RAPDs), and DNA sequence variation. We also give relevant information about the development times for these techniques, their relative costs, and the quality of tissue required. In addition, we discuss which conservation problems are appropriate to each of these methods and give examples of their application and potential use in relation to New Zealand organisms. The following problems are considered: sex assignment; parentage and kinship; migration; species, population and strain identification; forensic applications; genetic effects of population bottlenecks; disease identification; feeding preferences; philopatry; pest control; and understanding population extinction. Finally, we suggest that both microsatellite and minisatellite DNA techniques have particular advantages over many other currently-available techniques and conclude that these two approaches are applicable to a wide range of the conservation problems.

Single- and multilocus DNA fingerprinting of communally breeding pukeko: do copulations or dominance ensure reproductive success?

In behavioral and ecological studies the "function" of dominance hierarchies is thought to be related to reproductive success. In particular, dominant males are regarded as likely to gain a reproductive advantage due to enhanced "access" to females. We compare the dominance status of adults with the frequency of copulations and the patterns of parentage in communally breeding pukeko or purple swamphen (Porphyrio porphyrio melanotus). This avian species has an unusual social system, often breeding in polygynandrous groups in which there is a strong dominance hierarchy. Typically, during the breeding season, there is considerable sexual activity, with heterosexual and homosexual copulations between adults being commonplace. Hae III-digested DNA from individuals belonging to breeding groups was hybridized to the minisatellite DNA probe YNH24, revealing putative single-locus profiles, while hybridization of the same DNA to the minisatellite probes pV47-2, 3'HVR, and per revealed typical multilocus profiles. The numbers of unattributable restriction fragments allowed the maternity and paternity of all individuals to be conclusively determined, despite high levels of band sharing among individuals within breeding groups. These close genetic similarities are a likely consequence of strong philopatry and inbreeding. We report instances of males which are high on the dominance hierarchy but have only a limited reproductive output in comparison with others and males which are subordinate but achieve a significant proportion of fertilizations. Generally these data reveal no consistent relationship between dominance, the frequency of copulations, and parentage among males. We conclude that pukeko highlight some difficulties with conventional explanations for the "function" of dominance.