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.

Major Histocompatibility Complex Genes Map to Two Chromosomes in an Evolutionarily Ancient Reptile, the Tuatara Sphenodon punctatus.

Major histocompatibility complex (MHC) genes are a central component of the vertebrate immune system and usually exist in a single genomic region. However, considerable differences in MHC organization and size exist between different vertebrate lineages. Reptiles occupy a key evolutionary position for understanding how variation in MHC structure evolved in vertebrates, but information on the structure of the MHC region in reptiles is limited. In this study, we investigate the organization and cytogenetic location of MHC genes in the tuatara (Sphenodon punctatus), the sole extant representative of the early-diverging reptilian order Rhynchocephalia. Sequencing and mapping of 12 clones containing class I and II MHC genes from a bacterial artificial chromosome library indicated that the core MHC region is located on chromosome 13q. However, duplication and translocation of MHC genes outside of the core region was evident, because additional class I MHC genes were located on chromosome 4p. We found a total of seven class I sequences and 11 class II ß sequences, with evidence for duplication and pseudogenization of genes within the tuatara lineage. The tuatara MHC is characterized by high repeat content and low gene density compared with other species and we found no antigen processing or MHC framework genes on the MHC gene-containing clones. Our findings indicate substantial differences in MHC organization in tuatara compared with mammalian and avian MHCs and highlight the dynamic nature of the MHC. Further sequencing and annotation of tuatara and other reptile MHCs will determine if the tuatara MHC is representative of nonavian reptiles in general.

De novo sequence assembly and characterisation of a partial transcriptome for an evolutionarily distinct reptile, the tuatara (Sphenodon punctatus).

BACKGROUND: The tuatara (Sphenodon punctatus) is a species of extraordinary zoological interest, being the only surviving member of an entire order of reptiles which diverged early in amniote evolution. In addition to their unique phylogenetic placement, many aspects of tuatara biology, including temperature-dependent sex determination, cold adaptation and extreme longevity have the potential to inform studies of genome evolution and development. Despite increasing interest in the tuatara genome, genomic resources for the species are still very limited. We aimed to address this by assembling a transcriptome for tuatara from an early-stage embryo, which will provide a resource for genome annotation, molecular marker development and studies of development and adaptation in tuatara. RESULTS: We obtained 30 million paired-end 50 bp reads from an Illumina Genome Analyzer and assembled them with Velvet and Oases using a range of kmers. After removing redundancy and filtering out low quality transcripts, our transcriptome dataset contained 32911 transcripts, with an N50 of 675 and a mean length of 451 bp. Almost 50% (15965) of these transcripts could be annotated by comparison with the NCBI non-redundant (NR) protein database or the chicken, green anole and zebrafish UniGene sets. A scan of candidate genes and repetitive elements revealed genes involved in immune function, sex differentiation and temperature-sensitivity, as well as over 200 microsatellite markers. CONCLUSIONS: This dataset represents a major increase in genomic resources for the tuatara, increasing the number of annotated gene sequences from just 60 to almost 16,000. This will facilitate future research in sex determination, genome evolution, local adaptation and population genetics of tuatara, as well as inform studies on amniote evolution.

Securing the demographic and genetic future of tuatara through assisted colonization.

Climate change poses a particular threat to species with fragmented distributions and little or no capacity to migrate. Assisted colonization, moving species into regions where they have not previously occurred, aims to establish populations where they are expected to survive as climatic envelopes shift. However, adaptation to the source environment may affect whether species successfully establish in new regions. Assisted colonization has spurred debate among conservation biologists and ecologists over whether the potential benefits to the threatened species outweigh the potential disruption to recipient communities. In our opinion, the debate has been distracted by controversial examples, rather than cases where assisted colonization may be a viable strategy. We present a strategic plan for the assisted migration of tuatara (Sphenodon punctatus), an endemic New Zealand reptile. The plan includes use of extant populations as reference points for comparisons with assisted-colonization populations with respect to demography, phenotypic plasticity, and phenology; optimization of genetic variation; research to fill knowledge gaps; consideration of host and recipient communities; and inclusion of stakeholders in the planning stage. When strategically planned and monitored, assisted colonization could meet conservation and research goals and ultimately result in the establishment of long-term sustainable populations capable of persisting during rapid changes in climate.

Characterisation of class II B MHC genes from a ratite bird, the little spotted kiwi (Apteryx owenii).

Major histocompatibility complex (MHC) genes are important for vertebrate immune response and typically display high levels of diversity due to balancing selection from exposure to diverse pathogens. An understanding of the structure of the MHC region and diversity among functional MHC genes is critical to understanding the evolution of the MHC and species resilience to disease exposure. In this study, we characterise the structure and diversity of class II MHC genes in little spotted kiwi Apteryx owenii, a ratite bird representing the basal avian lineage (paleognaths). Results indicate that little spotted kiwi have a more complex MHC structure than that of other non-passerine birds, with at least five class II MHC genes, three of which are expressed and likely to be functional. Levels of MHC variation among little spotted kiwi are extremely low, with 13 birds assayed having nearly identical MHC genotypes (only two genotypes containing four alleles, three of which are fixed). These results suggest that recent genetic drift due to a species-wide bottleneck of at most seven birds has overwhelmed past selection for high MHC diversity in little spotted kiwi, potentially leaving the species highly susceptible to disease.

Genetic diversity and differentiation at MHC genes in island populations of tuatara (Sphenodon spp.).

Neutral genetic markers are commonly used to understand the effects of fragmentation and population bottlenecks on genetic variation in threatened species. Although neutral markers are useful for inferring population history, the analysis of functional genes is required to determine the significance of any observed geographical differences in variation. The genes of the major histocompatibility complex (MHC) are well-known examples of genes of adaptive significance and are particularly relevant to conservation because of their role in pathogen resistance. In this study, we survey diversity at MHC class I loci across a range of tuatara populations. We compare the levels of MHC variation with that observed at neutral microsatellite markers to determine the relative roles of balancing selection, diversifying selection and genetic drift in shaping patterns of MHC variation in isolated populations. In general, levels of MHC variation within tuatara populations are concordant with microsatellite variation. Tuatara populations are highly differentiated at MHC genes, particularly between the northern and Cook Strait regions, and a trend towards diversifying selection across populations was observed. However, overall our results indicate that population bottlenecks and isolation have a larger influence on patterns of MHC variation in tuatara populations than selection.

Influence of major histocompatibility complex genotype on mating success in a free-ranging reptile population.

Major histocompatibility complex (MHC) genes are highly polymorphic components of the vertebrate immune system, which play a key role in pathogen resistance. MHC genes may also function as odour-related cues for mate choice, thus ensuring optimal MHC diversity in offspring. MHC-associated mate choice has been demonstrated in some fish, bird and mammal species but it is not known whether this is a general vertebrate phenomenon. We investigated whether MHC-associated mate choice occurs in a wild population of tuatara (Sphenodon punctatus), a territorial and sexually dimorphic reptile. We found weak evidence for MHC-disassortative mating, based on amino acid genotypic distance between pairs, when mated pairs were directly compared with potential pairs in close spatial proximity. No significant association was found between male mating success, number of MHC sequences, microsatellite heterozygosity or MHC lineage. The major determinant of mating success in tuatara was male body size, which was not related to MHC lineage or microsatellite heterozygosity. Our results suggest that male competitive ability is the primary driver of mating success in tuatara. However, MHC-associated preferences also appear to play a role, possibly as a kin avoidance mechanism during territory formation.

Two patterns of variation among MHC class I loci in Tuatara (Sphenodon punctatus).

The genes of the major histocompatibility complex (MHC) are a central component of the immune system in vertebrates and have become important markers of functional, fitness-related genetic variation. We have investigated the evolutionary processes that generate diversity at MHC class I genes in a large population of an archaic reptile species, the tuatara (Sphenodon punctatus), found on Stephens Island, Cook Strait, New Zealand. We identified at least 2 highly polymorphic (UA type) loci and one locus (UZ) exhibiting low polymorphism. The UZ locus is characterized by low nucleotide diversity and weak balancing selection and may be either a nonclassical class I gene or a pseudogene. In contrast, the UA-type alleles have high nucleotide diversity and show evidence of balancing selection at putative peptide-binding sites. Twenty-one different UA-type genotypes were identified among 26 individuals, suggesting that the Stephens Island population has high levels of MHC class I variation. UA-type allelic diversity is generated by a mixture of point mutation and gene conversion. As has been found in birds and fish, gene conversion obscures the genealogical relationships among alleles and prevents the assignment of alleles to loci. Our results suggest that the molecular mechanisms that underpin MHC evolution in nonmammals make locus-specific amplification impossible in some species.

A molecular phylogeny of New Zealand's Petroica (Aves: Petroicidae) species based on mitochondrial DNA sequences.

The New Zealand robin (Petroica australis), tomtit (P. macrocephala), and Chatham Island black robin (P. traversi) are members of the Petroicidae family of Australo-Papuan robins, found throughout Australasia and the western Pacific. In the nearly 200 years since the New Zealand members of Petroicidae were first described, the division of species, subspecies, and even genera has undergone many changes. In this study, we investigate whether molecular phylogenies based on mitochondrial DNA sequences support current taxonomic classifications based on morphology. Petroica traversi, P. australis, and P. macrocephala form distinct clades in phylogenetic trees constructed from Cytochrome b and control region sequences, however the position of the black robin is at odds with the morphological and behavioral data. The black robin does not appear to be a derivative of the New Zealand robin, instead it groups strongly with the tomtit, indicating that lineage sorting and/or introgressive hybridization may have occurred. There is some evidence to support the hypothesis that two invasions of Petroica from Australia have occurred, however additional data from Australian Petroica taxa are required to confirm this. Control region sequences confirm a deep split between the North and South Island P. australis lineages, but suggest a recent radiation of P. macrocephala.

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. MHC Class I genes in the Tuatara (Sphenodon spp.): evolution of the MHC in an ancient reptilian order.

The major histocompatibility complex (MHC) is an extremely dynamic region of the genome, characterized by high polymorphism and frequent gene duplications and rearrangements. This has resulted in considerable differences in MHC organization and evolution among vertebrate lineages, particularly between birds and mammals. As nonavian reptiles are ancestral to both mammals and birds, they occupy an important phylogenetic position for understanding these differences. However, little is known about reptile MHC genes. To address this, we have characterized MHC class I sequences from the tuatara (Sphenodon spp.), the last survivor of an ancient order of reptiles, Sphenodontia. We isolated two different class I cDNA sequences, which share 93% sequence similarity with each other but are highly divergent from other vertebrate MHC genes. Southern blotting and polymerase chain reaction amplification of class I sequences from seven adult tuatara plus a family group indicate that these sequences represent at least two to three loci. Preliminary analysis of variation among individuals from an island population of tuatara indicates that these loci are highly polymorphic. Maximum likelihood analysis of reptile MHC class I sequences indicates that gene duplication has occurred within reptilian orders. However, the evolutionary relationships among sequences from different reptilian orders cannot be resolved, reflecting the antiquity of the major reptile lineages.

Characterization of MHC class II genes from an ancient reptile lineage, Sphenodon (tuatara).

The organization and evolution of major histocompatibility complex (MHC) genes vary considerably among vertebrate lineages. MHC genes have been well characterized in mammals, birds, amphibians and fish, but little is known about their organization in reptiles, despite the fact that reptiles occupy an important phylogenetic position for understanding the evolutionary history of both mammalian and avian MHC genes. Here we describe the characterization of the first MHC class II B cDNA sequences from a non-avian reptile, the tuatara (Sphenodon spp.). Three class II B sequences were isolated from a tuatara cDNA library, and four additional partial sequences were isolated by reverse transcriptase-polymerase chain reaction. Six of these sequences appear to belong to the same gene family, which we have named SppuDAB. The remaining sequence (named SppuDBB) shares only 43.9% amino acid similarity with SppuDAB and thus appears to represent a separate gene family. SppuDBB may be a non-classical locus as it does not contain all the conserved residues expected of a classical MHC class II gene. Southern blot analysis indicates that only a single copy of SppuDBB exists in tuatara, but that multiple loci related to SppuDAB are present. The SppuDAB sequences have the highest amino acid similarity (57.2-62.4%) with class II B sequences from the spectacled caiman, but only 26.4-48.7% similarity with sequences from other vertebrates. The tuatara sequences do not strongly group with other reptile sequences on a phylogenetic tree, reflecting the antiquity of the Sphenodon lineage and the lack of closely related sequences for comparison.

Genetic drift outweighs balancing selection in shaping post-bottleneck major histocompatibility complex variation in New Zealand robins (Petroicidae).

The Chatham Island black robin, Petroica traversi, is a highly inbred, endangered passerine with extremely low levels of variation at hypervariable neutral DNA markers. In this study we investigated variation in major histocompatibility complex (MHC) class II genes in both the black robin and its nonendangered relative, the South Island robin Petroica australis australis. Previous studies have shown that Petroica have at least four expressed class II B MHC genes. In this study, the sequences of introns flanking exon 2 of these loci were characterized to design primers for peptide-binding region (PBR) sequence analysis. Intron sequences were comprised of varying numbers of repeated units, with highly conserved regions immediately flanking exon 2. Polymerase chain reaction primers designed to this region amplified three or four sequences per black robin individual, and eight to 14 sequences per South Island robin individual. MHC genes are fitness-related genes thought to be under balancing selection, so they may be more likely to retain variation in bottlenecked populations. To test this, we compared MHC variation in the black robin with artificially bottlenecked populations of South Island robin, and with their respective source populations, using restriction fragment length polymorphism analyses and DNA sequencing of the PBR. Our results indicate that the black robin is monomorphic at class II B MHC loci, while both source and bottlenecked populations of South Island robin have retained moderate levels of variation. Comparison of MHC variation with minisatellite DNA variation indicates that genetic drift outweighs balancing selection in determining MHC diversity in the bottlenecked populations. However, balancing selection appears to influence MHC diversity over evolutionary timescales, and the effects of gene conversion are evident.

Gene duplication and gene conversion in class II MHC genes of New Zealand robins (Petroicidae).

In contrast to mammals, the evolution of MHC genes in birds appears to be characterized by high rates of gene duplication and concerted evolution. To further our understanding of the evolution of passerine MHC genes, we have isolated class II B sequences from two species of New Zealand robins, the South Island robin (Petroica australis australis), and the endangered Chatham Island black robin (Petroica traversi). Using an RT-PCR based approach we isolated four transcribed class II B MHC sequences from the black robin, and eight sequences from the South Island robin. RFLP analysis indicated that all class II B loci were contained within a single linkage group. Analysis of 3'-untranslated region sequences enabled putative orthologous loci to be identified in the two species, and indicated that multiple rounds of gene duplication have occurred within the MHC of New Zealand robins. The orthologous relationships are not retained within the coding region of the gene, instead the sequences group within species. A number of putative gene conversion events were identified across the length of our sequences that may account for this. Exon 2 sequences are highly diverse and appear to have diverged under balancing selection. It is also possible that gene conversion involving short stretches of sequence within exon 2 adds to this diversity. Our study is the first report of putative orthologous MHC loci in passerines, and provides further evidence for the importance of gene duplication and gene conversion in the evolution of the passerine MHC.