Journeying towards best practice data management in biodiversity genomics

Advances in sequencing technologies and declining costs are increasing the accessibility of large-scale biodiversity genomic datasets. To maximize the impact of these data, a careful, considered approach to data management is essential. However, challenges associated with the management of such datasets remain, exacerbated by uncertainty among the research community as to what constitutes best practices. As an interdisciplinary team with diverse data management experience, we recognize the growing need for guidance on comprehensive data management practices that minimize the risks of data loss, maximize efficiency for stand-alone projects, enhance opportunities for data reuse, facilitate Indigenous data sovereignty and uphold the FAIR and CARE Guiding Principles. Here, we describe four fictional personas reflecting differing user experiences with data management to identify data management challenges across the biodiversity genomics research ecosystem. We then use these personas to demonstrate realistic considerations, compromises and actions for biodiversity genomic data management. We also launch the Biodiversity Genomics Data Management Hub (https://genomicsaotearoa.github.io/data-management-resources/), containing tips, tricks and resources to support biodiversity genomics researchers, especially those new to data management, in their journey towards best practice. The Hub also provides an opportunity for those biodiversity researchers whose expertise lies beyond genomics and are keen to advance their data management journey. We aim to support the biodiversity genomics community in embedding data management throughout the research lifecycle to maximize research impact and outcomes.

Species-wide genomics of kākāpō provides tools to accelerate recovery

The kākāpō is a critically endangered, intensively managed, long-lived nocturnal parrot endemic to Aotearoa New Zealand. We generated and analysed whole-genome sequence data for nearly all individuals living in early 2018 (169 individuals) to generate a high-quality species-wide genetic variant callset. We leverage extensive long-term metadata to quantify genome-wide diversity of the species over time and present new approaches using probabilistic programming, combined with a phenotype dataset spanning five decades, to disentangle phenotypic variance into environmental and genetic effects while quantifying uncertainty in small populations. We find associations for growth, disease susceptibility, clutch size and egg fertility within genic regions previously shown to influence these traits in other species. Finally, we generate breeding values to predict phenotype and illustrate that active management over the past 45 years has maintained both genome-wide diversity and diversity in breeding values and, hence, evolutionary potential. We provide new pathways for informing future conservation management decisions for kākāpō, including prioritizing individuals for translocation and monitoring individuals with poor growth or high disease risk. Overall, by explicitly addressing the challenge of the small sample size, we provide a template for the inclusion of genomic data that will be transformational for species recovery efforts around the globe.

Fish germ cell cryobanking and transplanting for conservation

The unprecedented loss of global biodiversity is linked to multiple anthropogenic stressors. New conservation technologies are urgently needed to mitigate this loss. The rights, knowledge and perspectives of Indigenous peoples in biodiversity conservation-including the development and application of new technologies-are increasingly recognised. Advances in germplasm cryopreservation and germ cell transplantation (termed 'broodstock surrogacy') techniques offer exciting tools to preserve biodiversity, but their application has been underappreciated. Here, we use teleost fishes as an exemplar group to outline (1) the power of these techniques to preserve genome-wide genetic diversity, (2) the need to apply a conservation genomic lens when selecting individuals for germplasm cryobanking and broodstock surrogacy and (3) the value of considering the cultural significance of these genomic resources. We conclude by discussing the opportunities and challenges of these techniques for conserving biodiversity in threatened teleost fish and beyond.

The promise and challenges of characterizing genome-wide structural variants: A case study in a critically endangered parrot

There is growing interest in the role of structural variants (SVs) as drivers of local adaptation and speciation. From a biodiversity genomics perspective, the characterization of genome-wide SVs provides an exciting opportunity to complement single nucleotide polymorphisms (SNPs). However, little is known about the impacts of SV discovery and genotyping strategies on the characterization of genome-wide SV diversity within and among populations. Here, we explore a near whole-species resequence data set, and long-read sequence data for a subset of highly represented individuals in the critically endangered kākāpō (Strigops habroptilus). We demonstrate that even when using a highly contiguous reference genome, different discovery and genotyping strategies can significantly impact the type, size and location of SVs characterized genome-wide. Further, we found that the mean number of SVs in each of two kākāpō lineages differed both within and across generations. These combined results suggest that genome-wide characterization of SVs remains challenging at the population-scale. We are optimistic that increased accessibility to long-read sequencing and advancements in bioinformatic approaches including multireference approaches like genome graphs will alleviate at least some of the challenges associated with resolving SV characteristics below the species level. In the meantime, we address caveats, highlight considerations, and provide recommendations for the characterization of genome-wide SVs in biodiversity genomic research.

Hidden impacts of conservation management on fertility of the critically endangered kākāpō

Background

Animal conservation often requires intensive management actions to improve reproductive output, yet any adverse effects of these may not be immediately apparent, particularly in threatened species with small populations and long lifespans. Hand-rearing is an example of a conservation management strategy which, while boosting populations, can cause long-term demographic and behavioural problems. It is used in the recovery of the critically endangered kākāpō (Strigops habroptilus), a flightless parrot endemic to New Zealand, to improve the slow population growth that is due to infrequent breeding, low fertility and low hatching success.

Methods

We applied Bayesian mixed models to examine whether hand-rearing and other factors were associated with clutch fertility in kākāpō. We used projection predictive variable selection to compare the relative contributions to fertility from the parents' rearing environment, their age and previous copulation experience, the parental kinship, and the number of mates and copulations for each clutch. We also explored how the incidence of repeated copulations and multiple mates varied with kākāpō density.

Results

The rearing status of the clutch father and the number of mates and copulations of the clutch mother were the dominant factors in predicting fertility. Clutches were less likely to be fertile if the father was hand-reared compared to wild-reared, but there was no similar effect for mothers. Clutches produced by females copulating with different males were more likely to be fertile than those from repeated copulations with one male, which in turn had a higher probability of fertility than those from a single copulation. The likelihood of multiple copulations and mates increased with female:male adult sex ratio, perhaps as a result of mate guarding by females. Parental kinship, copulation experience and age all had negligible associations with clutch fertility.

Conclusions

These results provide a rare assessment of factors affecting fertility in a wild threatened bird species, with implications for conservation management. The increased fertility due to multiple mates and copulations, combined with the evidence for mate guarding and previous results of kākāpō sperm morphology, suggests that an evolutionary mechanism exists to optimise fertility through sperm competition in kākāpō. The high frequency of clutches produced from single copulations in the contemporary population may therefore represent an unnatural state, perhaps due to too few females. This suggests that opportunity for sperm competition should be maximised by increasing population densities, optimising sex ratios, and using artificial insemination. The lower fertility of hand-reared males may result from behavioural defects due to lack of exposure to conspecifics at critical development stages, as seen in other taxa. This potential negative impact of hand-rearing must be balanced against the short-term benefits it provides.

Leveraging an existing whole genome resequencing population dataset to characterize toll‐like receptor gene diversity in a threatened bird

Species recovery programs are increasingly using genomic data to measure neutral genetic diversity and calculate metrics like relatedness. While these measures can inform conservation management, determining the mechanisms underlying inbreeding depression requires information about functional genes associated with adaptive or maladaptive traits. Toll‐like receptors (TLRs) are one family of functional genes, which play a crucial role in recognition of pathogens and activation of the immune system. Previously, these genes have been analysed using species‐specific primers and PCR. Here, we leverage an existing short‐read reference genome, whole‐genome resequencing population data set, and bioinformatic tools to characterize TLR gene diversity in captive and wild tchūriwat’/tūturuatu/shore plover (Thinornis novaeseelandiae), a threatened bird endemic to Aotearoa New Zealand. Our results show that TLR gene diversity in tchūriwat’/tūturuatu is low, and forms two distinct captive and wild genetic clusters. The bioinformatic approach presented here has broad applicability to other threatened species with existing genomic resources in Aotearoa New Zealand and beyond.

Comparing genome-based estimates of relatedness for use in pedigree-based conservation management

Researchers have long debated which estimator of relatedness best captures the degree of relationship between two individuals. In the genomics era, this debate continues, with relatedness estimates being sensitive to the methods used to generate markers, marker quality, and levels of diversity in sampled individuals. Here, we compare six commonly used genome-based relatedness estimators (kinship genetic distance (KGD), Wang Maximum Likelihood (TrioML), Queller and Goodnight (R_xy ), Kinship INference for Genome-wide association studies (KING-robust), and Pairwise Relatedness (R_AB ), allele-sharing co-ancestry (AS)) across five species bred in captivity-including three birds and two mammals-with varying degrees of reliable pedigree data, using reduced-representation and whole genome resequencing data. Genome-based relatedness estimates varied widely across estimators, sequencing methods, and species, yet the most consistent results for known first order relationships were found using R_xy , R_AB , and AS. However, AS was found to be less consistently correlated with known pedigree relatedness than either R_xy or R_AB . Our combined results indicate there is not a single genome-based estimator that is ideal across different species and data types. To determine the most appropriate genome-based relatedness estimator for each new dataset, we recommend assessing the relative: (1) correlation of candidate estimators with known relationships in the pedigree and (2) precision of candidate estimators with known first-order relationships. These recommendations are broadly applicable to conservation breeding programs, particularly where genome-based estimates of relatedness can complement and complete poorly pedigreed populations. Given a growing interest in the application of wild pedigrees, our results are also applicable to in-situ wildlife management.

Weaving place‐based knowledge for culturally significant species in the age of genomics: Looking to the past to navigate the future

Relationships with place provide critical context for characterizing biocultural diversity. Yet, genetic and genomic studies are rarely informed by Indigenous or local knowledge, processes, and practices, including the movement of culturally significant species. Here, we show how place‐based knowledge can better reveal the biocultural complexities of genetic or genomic data derived from culturally significant species. As a case study, we focus on culturally significant southern freshwater kōura (crayfish) in Aotearoa me Te Waipounamu (New Zealand, herein Aotearoa NZ). Our results, based on genotyping‐by‐sequencing markers, reveal strong population genetic structure along with signatures of population admixture in 19 genetically depauperate populations across the east coast of Te Waipounamu. Environment association and differentiation analyses for local adaptation also indicate a role for hydroclimatic variables—including temperature, precipitation, and water flow regimes—in shaping local adaptation in kōura. Through trusted partnerships between community and researchers, weaving genomic markers with place‐based knowledge has both provided invaluable context for the interpretation of data and created opportunities to reconnect people and place. We envisage such trusted partnerships guiding future genomic research for culturally significant species in Aotearoa NZ and beyond.

The relevance of pedigrees in the conservation genomics era

Over the past 50 years conservation genetics has developed a substantive toolbox to inform species management. One of the most long‐standing tools available to manage genetics—the pedigree—has been widely used to characterize diversity and maximize evolutionary potential in threatened populations. Now, with the ability to use high throughput sequencing to estimate relatedness, inbreeding, and genome‐wide functional diversity, some have asked whether it is warranted for conservation biologists to continue collecting and collating pedigrees for species management. In this perspective, we argue that pedigrees remain a relevant tool, and when combined with genomic data, create an invaluable resource for conservation genomic management. Genomic data can address pedigree pitfalls (e.g., founder relatedness, missing data, uncertainty), and in return robust pedigrees allow for more nuanced research design, including well‐informed sampling strategies and quantitative analyses (e.g., heritability, linkage) to better inform genomic inquiry. We further contend that building and maintaining pedigrees provides an opportunity to strengthen trusted relationships among conservation researchers, practitioners, Indigenous Peoples, and Local Communities.

Expanding the conservation genomics toolbox: Incorporating structural variants to enhance genomic studies for species of conservation concern

Structural variants (SVs) are large rearrangements (>50 bp) within the genome that impact gene function and the content and structure of chromosomes. As a result, SVs are a significant source of functional genomic variation, that is, variation at genomic regions underpinning phenotype differences, that can have large effects on individual and population fitness. While there are increasing opportunities to investigate functional genomic variation in threatened species via single nucleotide polymorphism (SNP) data sets, SVs remain understudied despite their potential influence on fitness traits of conservation interest. In this future-focused Opinion, we contend that characterizing SVs offers the conservation genomics community an exciting opportunity to complement SNP-based approaches to enhance species recovery. We also leverage the existing literature-predominantly in human health, agriculture and ecoevolutionary biology-to identify approaches for readily characterizing SVs and consider how integrating these into the conservation genomics toolbox may transform the way we manage some of the world's most threatened species.

Designing monitoring protocols to measure population trends of threatened insects: A case study of the cryptic, flightless grasshopper Brachaspis robustus

Statistically robust monitoring of threatened populations is essential for effective conservation management because the population trend data that monitoring generates is often used to make decisions about when and how to take action. Despite representing the highest proportion of threatened animals globally, the development of best practice methods for monitoring populations of threatened insects is relatively uncommon. Traditionally, population trend data for the Nationally Endangered New Zealand grasshopper Brachaspis robustus has been determined by counting all adults and nymphs seen on a single ~1.5 km transect searched once annually. This method lacks spatial and temporal replication, both of which are essential to overcome detection errors in highly cryptic species like B. robustus. It also provides no information about changes in the grasshopper’s distribution throughout its range. Here, we design and test new population density and site occupancy monitoring protocols by comparing a) comprehensive plot and transect searches at one site and b) transect searches at two sites representing two different habitats (gravel road and natural riverbed) occupied by the species across its remaining range. Using power analyses, we determined a) the number of transects, b) the number of repeated visits and c) the grasshopper demographic to count to accurately detect long term change in relative population density. To inform a monitoring protocol design to track trends in grasshopper distribution, we estimated the probability of detecting an individual with respect to a) search area, b) weather and c) the grasshopper demographic counted at each of the two sites. Density estimates from plots and transects did not differ significantly. Population density monitoring was found to be most informative when large adult females present in early summer were used to index population size. To detect a significant change in relative density with power > 0.8 at the gravel road habitat, at least seventeen spatial replicates (transects) and four temporal replicates (visits) were required. Density estimates at the natural braided river site performed poorly and likely require a much higher survey effort. Detection of grasshopper presence was highest (p_g > 0.6) using a 100 m x 1 m transect at both sites in February under optimal (no cloud) conditions. At least three visits to a transect should be conducted per season for distribution monitoring. Monitoring protocols that inform the management of threatened species are crucial for better understanding and mitigation of the current global trends of insect decline. This study provides an exemplar of how appropriate monitoring protocols can be developed for threatened insect species.

A comparison of pedigree, genetic and genomic estimates of relatedness for informing pairing decisions in two critically endangered birds: Implications for conservation breeding programmes worldwide

Conservation management strategies for many highly threatened species include conservation breeding to prevent extinction and enhance recovery. Pairing decisions for these conservation breeding programmes can be informed by pedigree data to minimize relatedness between individuals in an effort to avoid inbreeding, maximize diversity and maintain evolutionary potential. However, conservation breeding programmes struggle to use this approach when pedigrees are shallow or incomplete. While genetic data (i.e., microsatellites) can be used to estimate relatedness to inform pairing decisions, emerging evidence indicates this approach may lack precision in genetically depauperate species, and more effective estimates will likely be obtained from genomic data (i.e., thousands of genome-wide single nucleotide polymorphisms, or SNPs). Here, we compare relatedness estimates and subsequent pairing decisions using pedigrees, microsatellites and SNPs from whole-genome resequencing approaches in two critically endangered birds endemic to New Zealand: kakī/black stilt (Himantopus novaezelandiae) and kākāriki karaka/orange-fronted parakeet (Cyanoramphus malherbi). Our findings indicate that SNPs provide more precise estimates of relatedness than microsatellites when assessing empirical parent-offspring and full sibling relationships. Further, our results show that relatedness estimates and subsequent pairing recommendations using PMx are most similar between pedigree- and SNP-based approaches. These combined results indicate that in lieu of robust pedigrees, SNPs are an effective tool for informing pairing decisions, which has important implications for many poorly pedigreed conservation breeding programmes worldwide.

Sperm competition risk drives rapid ejaculate adjustments mediated by seminal fluid

In many species, males can make rapid adjustments to ejaculate performance in response to sperm competition risk; however, the mechanisms behind these changes are not understood. Here, we manipulate male social status in an externally fertilising fish, chinook salmon (Oncorhynchus tshawytscha), and find that in less than 48 hr, males can upregulate sperm velocity when faced with an increased risk of sperm competition. Using a series of in vitro sperm manipulation and competition experiments, we show that rapid changes in sperm velocity are mediated by seminal fluid and the effect of seminal fluid on sperm velocity directly impacts paternity share and therefore reproductive success. These combined findings, completely consistent with sperm competition theory, provide unequivocal evidence that sperm competition risk drives plastic adjustment of ejaculate quality, that seminal fluid harbours the mechanism for the rapid adjustment of sperm velocity and that fitness benefits accrue to males from such adjustment.

Building strong relationships between conservation genetics and primary industry leads to mutually beneficial genomic advances

Several reviews in the past decade have heralded the benefits of embracing high-throughput sequencing technologies to inform conservation policy and the management of threatened species, but few have offered practical advice on how to expedite the transition from conservation genetics to conservation genomics. Here, we argue that an effective and efficient way to navigate this transition is to capitalize on emerging synergies between conservation genetics and primary industry (e.g., agriculture, fisheries, forestry and horticulture). Here, we demonstrate how building strong relationships between conservation geneticists and primary industry scientists is leading to mutually-beneficial outcomes for both disciplines. Based on our collective experience as collaborative New Zealand-based scientists, we also provide insight for forging these cross-sector relationships.

Conservation and losses of non-coding RNAs in avian genomes

Here we present the results of a large-scale bioinformatics annotation of non-coding RNA loci in 48 avian genomes. Our approach uses probabilistic models of hand-curated families from the Rfam database to infer conserved RNA families within each avian genome. We supplement these annotations with predictions from the tRNA annotation tool, tRNAscan-SE and microRNAs from miRBase. We identify 34 lncRNA-associated loci that are conserved between birds and mammals and validate 12 of these in chicken. We report several intriguing cases where a reported mammalian lncRNA, but not its function, is conserved. We also demonstrate extensive conservation of classical ncRNAs (e.g., tRNAs) and more recently discovered ncRNAs (e.g., snoRNAs and miRNAs) in birds. Furthermore, we describe numerous “losses” of several RNA families, and attribute these to either genuine loss, divergence or missing data. In particular, we show that many of these losses are due to the challenges associated with assembling avian microchromosomes. These combined results illustrate the utility of applying homology-based methods for annotating novel vertebrate genomes.

Development of polymorphic microsatellite markers for the New Zealand black stilt (Himantopus novaezelandiae) and cross-amplification in the pied stilt (Himantopus himantopus leucocephalus)

Eight polymorphic microsatellite primer pairs were developed for the critically endangered New Zealand black stilt, Himantopus novaezelandiae, representing the first microsatellite markers available for birds in the family Recurvirostridae. The number of alleles ranged from two to four per locus. Observed and expected heterozygosities ranged from 0.30 to 0.80 and from 0.37 to 0.70, respectively. All eight loci were polymorphic in the related species Himantopus himantopus leucocephalus, indicating these primer pairs may be useful for additional taxa in the globally distributed genus Himantopus.

Sampling for microsatellite-based population genetic studies: 25 to 30 individuals per population is enough to accurately estimate allele frequencies

One of the most common questions asked before starting a new population genetic study using microsatellite allele frequencies is "how many individuals do I need to sample from each population?" This question has previously been answered by addressing how many individuals are needed to detect all of the alleles present in a population (i.e. rarefaction based analyses). However, we argue that obtaining accurate allele frequencies and accurate estimates of diversity are much more important than detecting all of the alleles, given that very rare alleles (i.e. new mutations) are not very informative for assessing genetic diversity within a population or genetic structure among populations. Here we present a comparison of allele frequencies, expected heterozygosities and genetic distances between real and simulated populations by randomly subsampling 5-100 individuals from four empirical microsatellite genotype datasets (Formica lugubris, Sciurus vulgaris, Thalassarche melanophris, and Himantopus novaezelandia) to create 100 replicate datasets at each sample size. Despite differences in taxon (two birds, one mammal, one insect), population size, number of loci and polymorphism across loci, the degree of differences between simulated and empirical dataset allele frequencies, expected heterozygosities and pairwise F(ST) values were almost identical among the four datasets at each sample size. Variability in allele frequency and expected heterozygosity among replicates decreased with increasing sample size, but these decreases were minimal above sample sizes of 25 to 30. Therefore, there appears to be little benefit in sampling more than 25 to 30 individuals per population for population genetic studies based on microsatellite allele frequencies.

Genetic analyses reveal hybridization but no hybrid swarm in one of the world's rarest birds

Hybridization facilitated by human activities has dramatically altered the evolutionary trajectories of threatened taxa around the globe. Whereas introduced mammalian predators and widespread habitat loss and degradation clearly imperil the recovery and survival of the New Zealand endemic black stilt or kakī (Himantopus novaezelandiae), the risk associated with hybridization between this critically endangered endemic and its self-introduced congener, the pied stilt or poaka (Himantopus himantopus leucocephalus) is less clear. Here, we combine Bayesian admixture analyses of microsatellite data with mitochondrial DNA sequence data to assess the levels of hybridization and introgression between kakī and poaka. We show that birds classified as hybrids on the basis of adult plumage are indeed of hybrid origin and that hybridization between kakī and poaka is both extensive and bidirectional. Despite this, we found almost no evidence for introgression from poaka to kakī, thus negating the popular belief that kakī represent a hybrid swarm. To our knowledge, ours represents the first comprehensive study to document a lack of widespread introgression for a species at risk despite a recent history of extensive bidirectional human-induced hybridization. We attribute this rather surprising result, in part, to reduced reproductive success in female hybrids combined with a transient male-biased kakī sex ratio. To maximize the evolutionary potential of kakī, we use these data to recommend conservation management activities aimed to maintain the genetic integrity and to maximize the genetic diversity of this iconic rare bird.

Comparative phylogeography of brown (Sula leucogaster) and red-footed boobies (S. sula): the influence of physical barriers and habitat preference on gene flow in pelagic seabirds

To test the hypothesis that both physical and ecological barriers to gene flow drive population differentiation in tropical seabirds, we surveyed mitochondrial control region variation in 242 brown boobies (Sula leucogaster), which prefer inshore habitat, and 271 red-footed boobies (S. sula), which prefer pelagic habitat. To determine the relative influence of isolation and gene flow on population structure, we used both traditional methods and a recently developed statistical method based on coalescent theory and Bayesian inference (Isolation with Migration). We found that global population genetic structure was high in both species, and that female-mediated gene flow among ocean basins apparently has been restricted by major physical barriers including the Isthmus of Panama, and the periodic emergence of the Sunda and Sahul Shelves in Southeast Asia. In contrast, the evolutionary history of populations within ocean basins differed markedly between the two species. In brown boobies, we found high levels of population genetic differentiation and limited gene flow among colonies, even at spatial scales as small as 500km. Although red-footed booby colonies were also genetically differentiated within ocean basins, coalescent analyses indicated that populations have either diverged in the face of ongoing gene flow, or diverged without gene flow but recently made secondary contact. Regardless, gene flow among red-footed booby populations was higher than among brown booby populations. We suggest that these contrasting patterns of gene flow within ocean basins may be explained by the different habitat preferences of brown and red-footed boobies.

Merging ancient and modern DNA: extinct seabird taxon rediscovered in the North Tasman Sea

Ancient DNA has revolutionized the way in which evolutionary biologists research both extinct and extant taxa, from the inference of evolutionary history to the resolution of taxonomy. Here, we present, to our knowledge, the first study to report the rediscovery of an 'extinct' avian taxon, the Tasman booby (Sula tasmani), using classical palaeontological data combined with ancient and modern DNA data. Contrary to earlier work, we show an overlap in size between fossil and modern birds in the North Tasman Sea (classified currently as S. tasmani and Sula dactylatra fullagari, respectively). In addition, we show that Holocene fossil birds have mitochondrial control region sequences that are identical to those found in modern birds. These results indicate that the Tasman booby is not an extinct taxon: S. dactylatra fullagari O'Brien & Davies, 1990 is therefore a junior synonym of Sula tasmani van Tets, Meredith, Fullagar & Davidson, 1988 and all North Tasman Sea boobies should be known as S. d. tasmani. In addition to reporting the rediscovery of an extinct avian taxon, our study highlights the need for researchers to be cognizant of multidisciplinary approaches to understanding taxonomy and past biodiversity.

A role for nonphysical barriers to gene flow in the diversification of a highly vagile seabird, the masked booby (Sula dactylatra)

To test the hypothesis that nonphysical barriers to gene flow play a role in the divergence of low-latitude seabird populations, we applied phylogeographic methods to mitochondrial control region sequence variation in a global sample of masked boobies (Sula dactylatra). In accord with previous studies, we found that Indo-Pacific and Atlantic haplotypes form two divergent lineages, excluding one haplotype previously attributed to secondary contact between the Indian Ocean and the Caribbean Sea. Within the Indo-Pacific and the Atlantic, we found a relatively large number of haplotypes, many of which were unique to a single population. Although haplotypes from most populations were found in more than one higher-level clade, nested clade analysis revealed a significant association between clades and geography for the majority of higher-level clades, most often interpreted as a consequence of isolation by distance. We found low levels of gene flow within Indo-Pacific and Atlantic populations, and a significant correlation between gene flow and geographical distance among Indo-Pacific populations. We estimate that Indo-Pacific masked boobies experienced rapid population growth approximately 180,000 years ago and that the majority of Indo-Pacific and Atlantic populations diverged within the last approximately 115,000 years. These combined data suggest that the predominant pattern between Indo-Pacific and Atlantic populations is long-term isolation by physical barriers to gene flow. In contrast, populations within these regions appear to have diverged despite few obvious physical barriers to gene flow, perhaps as a consequence of limited natal dispersal combined with local adaptation and/or genetic drift.

The Isthmus of Panama: a major physical barrier to gene flow in a highly mobile pantropical seabird

To further test the hypothesis that the Isthmus of Panama is a major barrier to gene flow in pantropical seabirds, we applied phylogeographic methods to mitochondrial control sequence variation in masked booby (Sula dactylatra) populations on either side of the Isthmus of Panama and the southern tip of Africa. In accord with Steeves et al. (2003), we found that all Caribbean masked boobies with the 'secondary contact' cytochrome b haplotype (m-B) shared a control region haplotype (Sd_100), which grouped with Indian-Pacific haplotypes and not Caribbean-Atlantic haplotypes. In addition, Sd_100 was more closely related to control region haplotypes in the Indian Ocean than in the Pacific. We also found that the 'secondary contact' birds diverged more recently from extant populations in the Indian Ocean than in the Pacific. Thus, it appears that these masked boobies did not breach the Isthmus of Panama. Rather, birds likely dispersed around the southern tip of Africa during favourable oceanographic conditions in the Pleistocene.

In vitro development and nutrient uptake by embryos derived from oocytes of pre-pubertal and adult cows

The present study compared the developmental potential and uptake of nutrients by embryos from pre-pubertal and adult cows. Oocytes retrieved from ovaries of 5 to 7 month old calves and adult cows were matured and fertilized in vitro. Embryos were cultured in SOFaa to the blastocyst stage (7 days post-insemination). At successive stages of development, rates of glucose and pyruvate uptake were measured non-invasively by microfluorescence for individual embryos. Fertilization was equivalent in embryos from pre-pubertal and adult cows (P > 0.05), however development to blastocyst was significantly lower in embryos from pre-pubertal cows (9.8% versus 33.7%, respectively; P < 0.05). Total blastocyst cell number was not different between pre-pubertal and adult material (P > 0.05). Glucose uptake was exponential (pre-pubertal, r = 0.82; adult, r = 0. 82; P < 0.05), with an increase in uptake beyond the 8- to 16-cell stage. Glucose uptake was significantly lower in embryos from pre-pubertal cows at the 2- to 4-cell stages (1.5 versus 3.0 pmoles/embryo/hr; P < 0.05), but was equivalent to the adult cow at all other stages of development (P > 0.05). Pyruvate uptake was low until the blastocyst stage. Pyruvate uptake by embryos from pre-pubertal cows was significantly different to adult cows at the 1-cell stage (2.7 versus 4.6 pmoles/embryo/hr, respectively; P < 0. 05) and 2- to 4-cell stages (4.9 versus 3.6 pmoles/embryo/hr, respectively; P < 0.05). Pyruvate uptake was equivalent in the two groups in the later stages of development (P > 0.05). Perturbations in the uptake of nutrients by embryos from pre-pubertal cows were most likely due to the presence of a high proportion of developmentally incompetent embryos. Further, embryos from pre-pubertal cows that did develop to the blastocyst were as viable as blastocysts from adult cows with respect to nutrient uptakes and total cell number.

Metabolism of glucose, pyruvate, and glutamine during the maturation of oocytes derived from pre-pubertal and adult cows

The aim of the study was to compare the energy metabolism of oocytes from pre-pubertal (2 to 3 months) and adult cows during maturation, to identify the cause of poor developmental potential in many pre-pubertal oocytes. The metabolism of [5-(3)H] glucose, [2-(14)C] pyruvate, and [G-(3)H] glutamine was measured at 0 hr, 12 hr, and 24 hr maturation. Oxidative metabolism was important during maturation of oocytes from both pre-pubertal and adult cows, with pyruvate metabolism peaking at 12 hr and glutamine metabolism increasing linearly and peaking at 24 hr. Peak oxidative metabolism was significantly lower in oocytes from pre-pubertal animals, for both pyruvate and glutamine (P < 0.05). Glucose metabolism increased significantly during oocyte maturation in both groups (0hr to 24 hr). Glucose metabolism was significantly lower in oocytes from pre-pubertal cows at 12 hr (P < 0.05). Oocytes from pre-pubertal animals were significantly smaller than oocytes from adult cows at 0 hr, 12 hr, and 24 hr maturation (P < 0.05). When metabolic rates were corrected for oocyte volume, there were no significant differences in substrate metabolism between oocytes from pre-pubertal and adult cows. There was however, a delay in the increase in glucose metabolism in pre-pubertal oocytes 0 hr to 12 hr maturation. Germinal vesicle breakdown was slower in oocytes from pre-pubertal animals with more oocytes still at the germinal vesicle stage approximately 5 hr post-aspiration, compared to oocytes from adult cows (P < 0.05). By 24 hr, development to metaphase II was equivalent for pre-pubertal and adult oocytes. This study identified differences in energy metabolism, oocyte size, and meiotic progression between the oocytes from pre-pubertal and adult cows that may account for the poor developmental potential of many pre-pubertal oocytes.

Temporal and differential effects of amino acids on bovine embryo development in culture

The aim of the study was to determine the amino acid requirements of the in vitro-produced bovine embryo as it develops from the zygote to the blastocyst, using a two-step culture system. When added to synthetic oviduct fluid (SOF) for the first 72-h culture, Eagle's nonessential amino acids and glutamine (NeGln) significantly increased development to the 8- to 16-cell stage (Day 4 postinsemination [pi]) and subsequent blastocyst development (Day 7 pi). Glutamine alone during the first 72-h culture did not stimulate development to the 8- to 16-cell stage (p > 0.05); however, the removal of glutamine from NeGln reduced the stimulatory effects of the nonessential amino acids. Replacing glutamine with betaine (an organic osmolyte) in NeGln did not stimulate development to the 8- to 16-cell stage compared to culture in SOF, but it did improve subsequent blastocyst development, indicating an osmolytic function of glutamine during the first 72-h culture. The addition of Eagle's essential amino acids and glutamine to SOF, or to medium already containing nonessential amino acids and glutamine for the first 72-h culture, did not affect cleavage to the 8- to 16-cell stage or subsequent blastocyst development (p > 0.05). Beyond Day 4 pi, culture with 20aa (nonessential and essential amino acids and glutamine) increased blastocyst development, total cell number, and the number of cells in both the trophectoderm and inner cell mass, compared to culture with other groups of amino acids (p < 0.05). Substituting betaine for glutamine in 20aa reduced blastocyst formation, indicating a non-osmolytic function of glutamine during the second 72-h culture. Further, there was a significant negative correlation between the concentration of essential amino acids (quarter, half, or single strength) and embryo development during both the first 72-h and second 72-h culture (p < 0.01), indicating that the concentration of essential amino acids was too high during culture of the bovine embryo. This study identified the temporal and differential effects of amino acids during development of the bovine embryo from the zygote to the blastocyst.