Molecular characterisation and additional morphological descriptions of Eimeria spp. (Apicomplexa: Eimeriidae) from brown kiwi (Apteryx mantelli Bartlett).

Brown kiwi (Apteryx mantelli Bartlett), a ratite endemic to New Zealand, is currently listed as "Vulnerable" under the IUCN classification system due to predation by introduced mammals. Operation Nest Egg (ONE) raises chicks and juveniles in predator-proof enclosures until they are large enough to defend themselves. These facilities experience an environmental accumulation of coccidial oocysts, which leads to severe morbidity and mortality of these kiwi. Four species of coccidia have been morphologically described from sporulated oocysts with additional opportunistic descriptions of endogenous stages. This research continues the morphological descriptions of these species of Eimeria with an additional novel morphotype also morphologically described. It also provides the first genetic characterisation targeting the mitochondrial cytochrome c oxidase I (COI) gene. Based on these findings, it was determined there are at least five morphotypes of Eimeria that infect brown kiwi and co-infections are common at the ONE facilities surveyed. The COI amplicon targeted for this study was sufficient to provide differentiation from other members of this genus. Sanger sequencing yielded ambiguous bases, indicating the need for more in-depth sequencing.

Symmetry breaking in the embryonic skin triggers directional and sequential plumage patterning.

The development of an organism involves the formation of patterns from initially homogeneous surfaces in a reproducible manner. Simulations of various theoretical models recapitulate final states of natural patterns, yet drawing testable hypotheses from those often remains difficult. Consequently, little is known about pattern-forming events. Here, we surveyed plumage patterns and their emergence in Galliformes, ratites, passerines, and penguins, together representing the three major taxa of the avian phylogeny, and built a unified model that not only reproduces final patterns but also intrinsically generates shared and varying directionality, sequence, and duration of patterning. We used in vivo and ex vivo experiments to test its parameter-based predictions. We showed that directional and sequential pattern progression depends on a species-specific prepattern: an initial break in surface symmetry launches a travelling front of sharply defined, oriented domains with self-organising capacity. This front propagates through the timely transfer of increased cell density mediated by cell proliferation, which controls overall patterning duration. These results show that universal mechanisms combining prepatterning and self-organisation govern the timely emergence of the plumage pattern in birds.

Sequence comparisons of cytochrome P450 aromatases from Australian animals predict differences in enzymatic activity and/or efficiency†.

Aromatase (P450arom, CYP19A1) is the terminal enzyme in the synthesis of the steroid hormone family of estrogens. Not surprisingly, this enzyme has structural similarities between the limited number of species studied thus far. This study examined the structure of aromatases from four diverse Australian species including a marsupial (tammar wallaby; Macropus eugenii), monotreme (platypus; Ornithorhynchus anatinus), ratite (emu; Dromaius novaehollandiae) and lizard (bearded dragon; Pogona vitticeps). We successfully built homology models for each species, using the only crystallographically determined structure available, human aromatase. The amino acid sequences showed high amino acid sequence identity to the human aromatase: wallaby 81%, platypus 73%, emu 75% and bearded dragon at 74%. The overall structure was highly conserved among the five species, although there were non-secondary structures (loops and bends) that were variable and flexible that may result in some differences in catalytic activity. At the N-terminal regions, there were deletions and variations that suggest that functional distinctions may be found. We found that the active sites of all these proteins were identical, except for a slight variation in the emu. The electrostatic potential across the surfaces of these aromatases highlighted likely variations to the protein-protein interactions of these enzymes with both redox partner cytochrome P450 reductase and possibly homodimerization in the case of the platypus, which has been postulated for the human aromatase enzyme. Given the high natural selection pressures on reproductive strategies, the relatively high degree of conservation of aromatase sequence and structure across species suggests that there is biochemically very little scope for changes to have evolved without the loss of enzyme activity.

Whole-Genome Analyses Resolve the Phylogeny of Flightless Birds (Palaeognathae) in the Presence of an Empirical Anomaly Zone.

Palaeognathae represent one of the two basal lineages in modern birds, and comprise the volant (flighted) tinamous and the flightless ratites. Resolving palaeognath phylogenetic relationships has historically proved difficult, and short internal branches separating major palaeognath lineages in previous molecular phylogenies suggest that extensive incomplete lineage sorting (ILS) might have accompanied a rapid ancient divergence. Here, we investigate palaeognath relationships using genome-wide data sets of three types of noncoding nuclear markers, together totaling 20,850 loci and over 41 million base pairs of aligned sequence data. We recover a fully resolved topology placing rheas as the sister to kiwi and emu + cassowary that is congruent across marker types for two species tree methods (MP-EST and ASTRAL-II). This topology is corroborated by patterns of insertions for 4274 CR1 retroelements identified from multispecies whole-genome screening, and is robustly supported by phylogenomic subsampling analyses, with MP-EST demonstrating particularly consistent performance across subsampling replicates as compared to ASTRAL. In contrast, analyses of concatenated data supermatrices recover rheas as the sister to all other nonostrich palaeognaths, an alternative that lacks retroelement support and shows inconsistent behavior under subsampling approaches. While statistically supporting the species tree topology, conflicting patterns of retroelement insertions also occur and imply high amounts of ILS across short successive internal branches, consistent with observed patterns of gene tree heterogeneity. Coalescent simulations and topology tests indicate that the majority of observed topological incongruence among gene trees is consistent with coalescent variation rather than arising from gene tree estimation error alone, and estimated branch lengths for short successive internodes in the inferred species tree fall within the theoretical range encompassing the anomaly zone. Distributions of empirical gene trees confirm that the most common gene tree topology for each marker type differs from the species tree, signifying the existence of an empirical anomaly zone in palaeognaths.

Genetic parameters for tonic immobility, body weight, and morphological traits of the red-winged tinamou (Rhynchotus rufescens).

This study was carried out to estimate genetic parameters for morphology, body weight, and tonic immobility traits in the red-winged tinamou (Rhynchotus rufescens). Information on 690 birds was used and genetic parameters were estimated using Bayesian methods under a multi-trait animal model. The following traits were considered in this study: tarsal length (TL), bill length (BL), wing length (WL), head width (HW), bill width (BW), mature weight (MW), weight at 90 days (W90), and tonic immobility (TI). The heritability showed estimates between 0.15 for wing length and 0.56 for bill length. Positive and negative genetic correlations were estimated, ranging from - 0.33 to 0.81. All the morphological, production, and behavioral traits studied will have moderate to high response to selection. The body weight at 90 days is a better alternative for use in breeding programs and its selection would not lead to an increase in the time of tonic immobility. Both the selection for weight gain and for reduction of tonic immobility time would lead to an increase in the size of the legs of the red-winged tinamou, which could be advantageous for thermal control of these birds in tropical systems.

The phylogenetic significance of the morphology of the syrinx, hyoid and larynx, of the southern cassowary, Casuarius casuarius (Aves, Palaeognathae).

BACKGROUND: Palaeognathae is a basal clade within Aves and include the large and flightless ratites and the smaller, volant tinamous. Although much research has been conducted on various aspects of palaeognath morphology, ecology, and evolutionary history, there are still areas which require investigation. This study aimed to fill gaps in our knowledge of the Southern Cassowary, Casuarius casuarius, for which information on the skeletal systems of the syrinx, hyoid and larynx is lacking - despite these structures having been recognised as performing key functional roles associated with vocalisation, respiration and feeding. Previous research into the syrinx and hyoid have also indicated these structures to be valuable for determining evolutionary relationships among neognath taxa, and thus suggest they would also be informative for palaeognath phylogenetic analyses, which still exhibits strong conflict between morphological and molecular trees. RESULTS: The morphology of the syrinx, hyoid and larynx of C. casuarius is described from CT scans. The syrinx is of the simple tracheo-bronchial syrinx type, lacking specialised elements such as the pessulus; the hyoid is relatively short with longer ceratobranchials compared to epibranchials; and the larynx is comprised of entirely cartilaginous, standard avian anatomical elements including a concave, basin-like cricoid and fused cricoid wings. As in the larynx, both the syrinx and hyoid lack ossification and all three structures were most similar to Dromaius. We documented substantial variation across palaeognaths in the skeletal character states of the syrinx, hyoid, and larynx, using both the literature and novel observations (e.g. of C. casuarius). Notably, new synapomorphies linking Dinornithiformes and Tinamidae are identified, consistent with the molecular evidence for this clade. These shared morphological character traits include the ossification of the cricoid and arytenoid cartilages, and an additional cranial character, the articulation between the maxillary process of the nasal and the maxilla. CONCLUSION: Syrinx, hyoid and larynx characters of palaeognaths display greater concordance with molecular trees than do other morphological traits. These structures might therefore be less prone to homoplasy related to flightlessness and gigantism, compared to typical morphological traits emphasised in previous phylogenetic studies.

Explosive ice age diversification of kiwi.

Molecular dating largely overturned the paradigm that global cooling during recent Pleistocene glacial cycles resulted in a burst of species diversification although some evidence exists that speciation was commonly promoted in habitats near the expanding and retracting ice sheets. Here, we used a genome-wide dataset of more than half a million base pairs of DNA to test for a glacially induced burst of diversification in kiwi, an avian family distributed within several hundred kilometers of the expanding and retracting glaciers of the Southern Alps of New Zealand. By sampling across the geographic range of the five kiwi species, we discovered many cryptic lineages, bringing the total number of kiwi taxa that currently exist to 11 and the number that existed just before human arrival to 16 or 17. We found that 80% of kiwi diversification events date to the major glacial advances of the Middle and Late Pleistocene. During this period, New Zealand was repeatedly fragmented by glaciers into a series of refugia, with the tiny geographic ranges of many kiwi lineages currently distributed in areas adjacent to these refugia. Estimates of effective population size through time show a dramatic bottleneck during the last glacial cycle in all but one kiwi lineage, as expected if kiwi were isolated in glacially induced refugia. Our results support a fivefold increase in diversification rates during key glacial periods, comparable with levels observed in classic adaptive radiations, and confirm that at least some lineages distributed near glaciated regions underwent rapid ice age diversification.

Cryptic inbreeding depression in a growing population of a long-lived species.

Genetic effects are often overlooked in endangered species monitoring, and populations showing positive growth are often assumed to be secure. However, the continued reproductive success of a few individuals may mask issues such as inbreeding depression, especially in long-lived species. Here, we test for inbreeding depression in little spotted kiwi (Apteryx owenii) by comparing a population founded with two birds to one founded with 40 birds, both from the same source population and both showing positive population growth. We used a combination of microsatellite genotypes, nest observations and modelling to examine the consequences of assessing population viability exclusively via population growth. We demonstrate (i) significantly lower hatching success despite significantly higher reproductive effort in the population with two founders; (ii) positive growth in the population with two founders is mainly driven by ongoing chick production of the founding pair; and (iii) a substantial genetic load in the population founded with two birds (10-15 diploid lethal equivalents). Our results illustrate that substantial, cryptic inbreeding depression may still be present when a population is growing, especially in long-lived species with overlapping generations.

Sixteen kiwi (Apteryx spp) transcriptomes provide a wealth of genetic markers and insight into sex chromosome evolution in birds.

BACKGROUND: Kiwi represent the most basal extant avian lineage (paleognaths) and exhibit biological attributes that are unusual or extreme among living birds, such as large egg size, strong olfaction, nocturnality, flightlessness and long lifespan. Despite intense interest in their evolution and their threatened status, genomic resources for kiwi were virtually non-existent until the recent publication of a single genome. Here we present the most comprehensive kiwi transcriptomes to date, obtained via Illumina sequencing of whole blood and de novo assembly of mRNA sequences of eight individuals from each of the two rarest kiwi species, little spotted kiwi (LSK; Apteryx owenii) and rowi (A. rowi). RESULTS: Sequences obtained were orthologous with a wide diversity of functional genes despite the sequencing of a single tissue type. Individual and composite assemblies contain more than 7900 unique protein coding transcripts in each of LSK and rowi that show strong homology with chicken (Gallus gallus), including those associated with growth, development, disease resistance, reproduction and behavior. The assemblies also contain 66,909 SNPs that distinguish between LSK and rowi, 12,384 SNPs among LSK (associated with 3088 genes), and 29,313 SNPs among rowi (associated with 4953 genes). We found 3084 transcripts differentially expressed between LSK and rowi and 150 transcripts differentially expressed between the sexes. Of the latter, 83 could be mapped to chicken chromosomes with 95% syntenic with chromosome Z. CONCLUSIONS: Our study has simultaneously sequenced multiple species, sexes, and individual kiwi at thousands of genes, and thus represents a significant leap forward in genomic resources available for kiwi. The expression pattern we observed among chromosome Z related genes in kiwi is similar to that observed in ostriches and emu, suggesting a common and ancestral pattern of sex chromosome homomorphy, recombination, and gene dosage among living paleognaths. The transcriptome assemblies described here will provide a rich resource for polymorphic marker development and studies of adaptation of these highly unusual and endangered birds.

Genomic support for a moa-tinamou clade and adaptive morphological convergence in flightless ratites.

One of the most startling discoveries in avian molecular phylogenetics is that the volant tinamous are embedded in the flightless ratites, but this topology remains controversial because recent morphological phylogenies place tinamous as the closest relative of a monophyletic ratite clade. Here, we integrate new phylogenomic sequences from 1,448 nuclear DNA loci totaling almost 1 million bp from the extinct little bush moa, Chilean tinamou, and emu with available sequences from ostrich, elegant crested tinamou, four neognaths, and the green anole. Phylogenetic analysis using standard homogeneous models and heterogeneous models robust to common topological artifacts recovered compelling support for ratite paraphyly with the little bush moa closest to tinamous within ratites. Ratite paraphyly was further corroborated by eight independent CR1 retroposon insertions. Analysis of morphological characters reinterpreted on a 27-gene paleognath topology indicates that many characters are convergent in the ratites, probably as the result of adaptation to a cursorial life style.

Ancient DNA microsatellite analyses of the extinct New Zealand giant moa (Dinornis robustus) identify relatives within a single fossil site.

By analysing ancient DNA (aDNA) from 74 (14)C-dated individuals of the extinct South Island giant moa (Dinornis robustus) of New Zealand, we identified four dyads of closely related adult females. Although our total sample included bones from four fossil deposits located within a 10 km radius, these eight individuals had all been excavated from the same locality. Indications of kinship were based on high pairwise genetic relatedness (rXY) in six microsatellite markers genotyped from aDNA, coupled with overlapping radiocarbon ages. The observed rXY values in the four dyads exceeded a conservative cutoff value for potential relatives obtained from simulated data. In three of the four dyads, the kinship was further supported by observing shared and rare mitochondrial haplotypes. Simulations demonstrated that the proportion of observed dyads above the cutoff value was at least 20 times higher than expected in a randomly mating population with temporal sampling, also when introducing population structure in the simulations. We conclude that the results must reflect social structure in the moa population and we discuss the implications for future aDNA research.

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.

Genomic structures of the kW1 loci on the Z and W chromosomes in ratite birds: structural changes at an early stage of W chromosome differentiation.

The W chromosome of ratite birds shows minimal morphological differentiation and retains homology of genetic linkage and gene order with a substantial stretch of the Z chromosome; however, the molecular structure in the differentiated region is still not well known. The kW1 sequence was isolated from the kiwi as a W-specific DNA marker for PCR-based molecular sexing of ratite birds. In ratite W chromosomes, this sequence commonly contains a ∼200-bp deletion. To characterize the very early event of avian sex chromosome differentiation, we performed molecular cytogenetic analyses of kW1 and its flanking sequences in paleognathous and neognathous birds and reptiles. Female-specific repeats were found in the kW1-flanking sequence of the cassowary (Casuarius casuarius), and the repeats have been amplified in the pericentromeric region of the W chromosomes of ratites, which may have resulted from the cessation of meiotic recombination between the Z and W chromosomes at an early stage of sex chromosome differentiation. The presence of the kW1 sequence in neognathous birds and a crocodilian species suggests that the kW1 sequence was present in the ancestral genome of Archosauria; however, it disappeared in other reptilian taxa and several lineages of neognathous birds.

Ratite nonmonophyly: independent evidence from 40 novel Loci.

Large-scale multilocus studies have become common in molecular phylogenetics, but the best way to interpret these studies when their results strongly conflict with prior information about phylogeny remains unclear. An example of such a conflict is provided by the ratites (the large flightless birds of southern land masses, including ostriches, emus, and rheas). Ratite monophyly is strongly supported by both morphological data and many earlier molecular studies and is used as a textbook example of vicariance biogeography. However, recent studies have indicated that ratites are not monophyletic; instead, the volant tinamous nest inside the ratites rather than forming their sister group within the avian superorder Palaeognathae. Large-scale studies can exhibit biases that reflect a number of factors, including limitations in the fit of the evolutionary models used for analyses and problems with sequence alignment, so the unexpected conclusion that ratites are not monophyletic needs to be rigorously evaluated. A rigorous approach to testing novel hypotheses generated by large-scale studies is to collect independent evidence (i.e., excluding the loci and/or traits used to generate the hypotheses). We used 40 nuclear loci not used in previous studies that investigated the relationship among ratites and tinamous. Our results strongly support the recent molecular studies, revealing that the deepest branch within Palaeognathae separates the ostrich from other members of the clade, rather than the traditional hypothesis that separates the tinamous from the ratites. To ensure these results reflected evolutionary history, we examined potential biases in types of loci used, heterotachy, alignment biases, and discordance between gene trees and the species tree. All analyses consistently supported nonmonophyly of the ratites and no confounding biases were identified. This confirmation that ratites are not monophyletic using independent evidence will hopefully stimulate further comparative research on paleognath development and genetics that might reveal the basis of the morphological convergence in these large, flightless birds.

Evolution of bird sex chromosomes: a cytogenomic approach in Palaeognathae species.

BACKGROUND: Different patterns of sex chromosome differentiation are seen in Palaeognathae birds, a lineage that includes the ratites (Struthioniformes, Rheiformes, Apterygiformes, Casuariiformes, and the sister group Tinamiformes). While some Tinamiform species have well-differentiated W chromosomes, both Z and W of all the flightless ratites are still morphologically undifferentiated. Here, we conducted a comprehensive analysis of the ZW differentiation in birds using a combination of cytogenetic, genomic, and bioinformatic approaches. The whole set of satDNAs from the emu (Dromaius novaehollandiae) was described and characterized. Furthermore, we examined the in situ locations of these satDNAs alongside several microsatellite repeats and carried out Comparative Genomic Hybridizations in two related species: the greater rhea (Rhea americana) and the tataupa tinamou (Crypturellus tataupa). RESULTS: From the 24 satDNA families identified (which represent the greatest diversity of satDNAs ever uncovered in any bird species), only three of them were found to accumulate on the emu's sex chromosomes, with no discernible accumulation observed on the W chromosome. The W chromosomes of both the greater rhea and the emu did not exhibit a significant buildup of either C-positive heterochromatin or repetitive DNAs, indicating their large undifferentiation both at morphological and molecular levels. In contrast, the tataupa tinamou has a highly differentiated W chromosome that accumulates several DNA repeats. CONCLUSION: The findings provide new information on the architecture of the avian genome and an inside look at the starting points of sex chromosome differentiation in birds.

Genetic consequences of a century of protection: serial founder events and survival of the little spotted kiwi (Apteryx owenii).

We present the outcome of a century of post-bottleneck isolation of a long-lived species, the little spotted kiwi (Apteryx owenii, LSK) and demonstrate that profound genetic consequences can result from protecting few individuals in isolation. LSK were saved from extinction by translocation of five birds from South Island, New Zealand to Kapiti Island 100 years ago. The Kapiti population now numbers some 1200 birds and provides founders for new populations. We used 15 microsatellite loci to compare genetic variation among Kapiti LSK and the populations of Red Mercury, Tiritiri Matangi and Long Islands that were founded with birds from Kapiti. Two LSK native to D'Urville Island were also placed on Long Island. We found extremely low genetic variation and signatures of acute and recent genetic bottleneck effects in all four populations, indicating that LSK have survived multiple genetic bottlenecks. The Long Island population appears to have arisen from a single mating pair from Kapiti, suggesting there is no genetic contribution from D'Urville birds among extant LSK. The Ne/NC ratio of Kapiti Island LSK (0.03) is exceptionally low for terrestrial vertebrates and suggests that genetic diversity might still be eroding in this population, despite its large census size.

The evolution of sexual dimorphism in New Zealand giant moa (Dinornis) and other ratites.

The extinct giant moa Dinornis is one of the most remarkable known examples of reversed sexual size dimorphism (RSD), with males weighing 34-85 kg, but females weighing up to 240 kg. However, there has been little consideration of the evolutionary mechanism that produced this level of dimorphism, and most living palaeognaths also exhibit varying levels of RSD. Using male and female body mass data for extant ratites and tinamous and four extinct moa genera, and tests of phylogenetic dependence (λ) of body size evolution among these species, we investigated whether Dinornis was truly unusual with respect to RSD relative to other palaeognaths, which sex was under greater pressure to change in size over evolutionary time, and which candidate hypotheses explaining the presence and variability of RSD in the genus are most plausible. We demonstrate that the extreme level of RSD exhibited by Dinornis represents a straightforward consequence of positive allometric scaling of body size. However, Dinornis females have undergone more evolutionary change than males, and larger females from high-productivity environments are associated with greater differentiation, possibly driven by intraspecific competition and female-biased selection for increased offspring investment.

Simulating retention of rare alleles in small populations to assess management options for species with different life histories.

Preserving allelic diversity is important because it provides the capacity for adaptation and thus enables long-term population viability. Allele retention is difficult to predict in animals with overlapping generations, so we used a new computer model to simulate retention of rare alleles in small populations of 3 species with contrasting life-history traits: North Island Brown Kiwi (Apteryx mantelli; monogamous, long-lived), North Island Robins (Petroica longipes; monogamous, short-lived), and red deer (Cervus elaphus; polygynous, moderate lifespan). We simulated closed populations under various demographic scenarios and assessed the amounts of artificial immigration needed to achieve a goal of retaining 90% of selectively neutral rare alleles (frequency in the source population = 0.05) after 10 generations. The number of immigrants per generation required to meet the genetic goal ranged from 11 to 30, and there were key similarities and differences among species. None of the species met the genetic goal without immigration, and red deer lost the most allelic diversity due to reproductive skew among polygynous males. However, red deer required only a moderate rate of immigration relative to the other species to meet the genetic goal because nonterritorial breeders had a high turnover. Conversely, North Island Brown Kiwi needed the most immigration because the long lifespan of locally produced territorial breeders prevented a large proportion of immigrants from recruiting. In all species, the amount of immigration needed generally decreased with an increase in carrying capacity, survival, or reproductive output and increased as individual variation in reproductive success increased, indicating the importance of accurately quantifying these parameters to predict the effects of management. Overall, retaining rare alleles in a small, isolated population requires substantial investment of management effort. Use of simulations to explore strategies optimized for the populations in question will help maximize the value of this effort..

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.

Effect of Different Types of Sequence Data on Palaeognath Phylogeny.

Palaeognathae consists of five groups of extant species: flighted tinamous (1) and four flightless groups: kiwi (2), cassowaries and emu (3), rheas (4), and ostriches (5). Molecular studies supported the groupings of extinct moas with tinamous and elephant birds with kiwi as well as ostriches as the group that diverged first among the five groups. However, phylogenetic relationships among the five groups are still controversial. Previous studies showed extensive heterogeneity in estimated gene tree topologies from conserved nonexonic elements, introns, and ultraconserved elements. Using the noncoding loci together with protein-coding loci, this study investigated the factors that affected gene tree estimation error and the relationships among the five groups. Using closely related ostrich rather than distantly related chicken as the outgroup, concatenated and gene tree-based approaches supported rheas as the group that diverged first among groups (1)-(4). Whereas gene tree estimation error increased using loci with low sequence divergence and short length, topological bias in estimated trees occurred using loci with high sequence divergence and/or nucleotide composition bias and heterogeneity, which more occurred in trees estimated from coding loci than noncoding loci. Regarding the relationships of (1)-(4), the site patterns by parsimony criterion appeared less susceptible to the bias than tree construction assuming stationary time-homogeneous model and suggested the clustering of kiwi and cassowaries and emu the most likely with ∼40% support rather than the clustering of kiwi and rheas and that of kiwi and tinamous with 30% support each.