The first linkage map for Australo-Papuan Treefrogs (family: Pelodryadidae) reveals the sex-determination system of the Green-eyed Treefrog (Litoria serrata)

Amphibians represent a useful taxon to study the evolution of sex determination because of their highly variable sex-determination systems. However, the sex-determination system for many amphibian families remains unknown, in part because of a lack of genomic resources. Here, using an F1 family of Green-eyed Treefrogs (Litoria serrata), we produce the first genetic linkage map for any Australo-Papuan Treefrogs (family: Pelodryadidae). The resulting linkage map contains 8662 SNPs across 13 linkage groups. Using an independent set of sexed adults, we identify a small region in linkage group 6 matching an XY sex-determination system. These results suggest Litoria serrata possesses a male heterogametic system, with a candidate sex-determination locus on linkage group 6. Furthermore, this linkage map represents the first genomic resource for Australo-Papuan Treefrogs, an ecologically diverse family of over 220 species.

SNP data reveals the complex and diverse evolutionary history of the blue-ringed octopus genus (Octopodidae: Hapalochlaena) in the Asia-Pacific

The blue-ringed octopus species complex (Hapalochlaena spp.), known to occur from Southern Australia to Japan, currently contains four formally described species (Hapalochlaena maculosa, Hapalochlaena fasciata, Hapalochlaena lunulata and Hapalochlaena nierstraszi). These species are distinguished based on morphological characters (iridescent blue rings and/or lines) along with reproductive strategies. However, the observation of greater morphological diversity than previously captured by the current taxonomic framework indicates that a revision is required. To examine species boundaries within the genus we used mitochondrial (12S rRNA, 16S rRNA, cytochrome c oxidase subunit 1 [COI], cytochrome c oxidase subunit 3 [COIII] and cytochrome b [Cytb]) and genome-wide SNP data (DaRT seq) from specimens collected across its geographic range including variations in depth from 3m to >100m. This investigation indicates substantially greater species diversity present within the genus Hapalochlaena than is currently described. We identified 10,346 SNPs across all locations, which when analysed support a minimum of 11 distinct clades. Bayesian phylogenetic analysis of the mitochondrial COI gene on a more limited sample set dates the diversification of the genus to ∼30mya and corroborates eight of the lineages indicated by the SNP analyses. Furthermore, we demonstrate that the diagnostic lined patterning of H. fasciata found in North Pacific waters and NSW, Australia is polyphyletic and therefore likely the result of convergent evolution. Several "deep water" (> 100m) lineages were also identified in this study with genetic convergence likely to be driven by external selective pressures. Examination of morphological traits, currently being undertaken in a parallel morphological study, is required to describe additional species within the complex.

Genome assembly of the Australian black tiger shrimp (Penaeus monodon) reveals a novel fragmented IHHNV EVE sequence

Shrimp are a valuable aquaculture species globally; however, disease remains a major hindrance to shrimp aquaculture sustainability and growth. Mechanisms mediated by endogenous viral elements have been proposed as a means by which shrimp that encounter a new virus start to accommodate rather than succumb to infection over time. However, evidence on the nature of such endogenous viral elements and how they mediate viral accommodation is limited. More extensive genomic data on Penaeid shrimp from different geographical locations should assist in exposing the diversity of endogenous viral elements. In this context, reported here is a PacBio Sequel-based draft genome assembly of an Australian black tiger shrimp (Penaeus monodon) inbred for 1 generation. The 1.89 Gbp draft genome is comprised of 31,922 scaffolds (N50: 496,398 bp) covering 85.9% of the projected genome size. The genome repeat content (61.8% with 30% representing simple sequence repeats) is almost the highest identified for any species. The functional annotation identified 35,517 gene models, of which 25,809 were protein-coding and 17,158 were annotated using interproscan. Scaffold scanning for specific endogenous viral elements identified an element comprised of a 9,045-bp stretch of repeated, inverted, and jumbled genome fragments of infectious hypodermal and hematopoietic necrosis virus bounded by a repeated 591/590 bp host sequence. As only near complete linear ∼4 kb infectious hypodermal and hematopoietic necrosis virus genomes have been found integrated in the genome of P. monodon previously, its discovery has implications regarding the validity of PCR tests designed to specifically detect such linear endogenous viral element types. The existence of joined inverted infectious hypodermal and hematopoietic necrosis virus genome fragments also provides a means by which hairpin double-stranded RNA could be expressed and processed by the shrimp RNA interference machinery.

Population Structure, Genetic Connectivity, and Signatures of Local Adaptation of the Giant Black Tiger Shrimp (Penaeus monodon) throughout the Indo-Pacific Region

The giant black tiger shrimp (Penaeus monodon) is native to the Indo-Pacific and is the second most farmed penaeid shrimp species globally. Understanding genetic structure, connectivity, and local adaptation among Indo-Pacific black tiger shrimp populations is important for informing sustainable fisheries management and aquaculture breeding programs. Population genetic and outlier detection analyses were undertaken using 10,593 genome-wide single nucleotide polymorphisms (SNPs) from 16 geographically disparate Indo-Pacific P. monodon populations. Levels of genetic diversity were highest for Southeast Asian populations and were lowest for Western Indian Ocean (WIO) populations. Both neutral (n = 9,930) and outlier (n = 663) loci datasets revealed a pattern of strong genetic structure of P. monodon corresponding with broad geographical regions and clear genetic breaks among samples within regions. Neutral loci revealed seven genetic clusters and the separation of Fiji and WIO clusters from all other clusters, whereas outlier loci revealed six genetic clusters and high genetic differentiation among populations. The neutral loci dataset estimated five migration events that indicated migration to Southeast Asia from the WIO, with partial connectivity to populations in both oceans. We also identified 26 putatively adaptive SNPs that exhibited significant Pearson correlation (P < 0.05) between minor allele frequency and maximum or minimum sea surface temperature. Matched transcriptome contig annotations suggest putatively adaptive SNPs involvement in cellular and metabolic processes, pigmentation, immune response, and currently unknown functions. This study provides novel genome-level insights that have direct implications for P. monodon aquaculture and fishery management practices.

Comparing Genomic Signatures of Selection Between the Abbassa Strain and Eight Wild Populations of Nile Tilapia (Oreochromis niloticus) in Egypt

Domestication to captive rearing conditions, along with targeted selective breeding have genetic consequences that vary from those in wild environments. Nile tilapia (Oreochromis niloticus) is one of the most translocated and farmed aquaculture species globally, farmed throughout Asia, North and South America, and its African native range. In Egypt, a breeding program established the Abbassa Strain of Nile tilapia (AS) in 2002 based on local broodstock sourced from the Nile River. The AS has been intensively selected for growth and has gone through genetic bottlenecks which have likely shifted levels and composition of genetic diversity within the strain. Consequently, there are questions on the possible genetic impact AS escapees may have on endemic populations of Nile tilapia. However, to date there have been no genetic studies comparing genetic changes in the domesticated AS to local wild populations. This study used 9,827 genome-wide SNPs to investigate population genetic structure and signatures of selection in the AS (generations 9–11) and eight wild Nile tilapia populations from Egypt. SNP analyses identified two major genetic clusters (captive and wild populations), with wild populations showing evidence of isolation-by-distance among the Nile Delta and upstream riverine populations. Between genetic clusters, approximately 6.9% of SNPs were identified as outliers with outliers identified on all 22 O. niloticus chromosomes. A lack of localized outlier clustering on the genome suggests that no genes of major effect were presently detected. The AS has retained high levels of genetic diversity (H_{o_All} = 0.21 ± 0.01; H_{e_All} = 0.23 ± 0.01) when compared to wild populations (H_{o_All} = 0.18 ± 0.01; H_{e_All} = 0.17 ± 0.01) after 11 years of domestication and selective breeding. Additionally, 565 SNPs were unique within the AS line. While these private SNPs may be due to domestication signals or founder effects, it is suspected that introgression with blue tilapia (Oreochromis aureus) has occurred. This study highlights the importance of understanding the effects of domestication in addition to wild population structure to inform future management and dissemination decisions. Furthermore, by conducting a baseline genetic study of wild populations prior to the dissemination of a domestic line, the effects of aquaculture on these populations can be monitored over time.

Fine-scale population structure and evidence for local adaptation in Australian giant black tiger shrimp (Penaeus monodon) using SNP analysis

BACKGROUND

Restrictions to gene flow, genetic drift, and divergent selection associated with different environments are significant drivers of genetic differentiation. The black tiger shrimp (Penaeus monodon), is widely distributed throughout the Indian and Pacific Oceans including along the western, northern and eastern coastline of Australia, where it is an important aquaculture and fishery species. Understanding the genetic structure and the influence of environmental factors leading to adaptive differences among populations of this species is important for farm genetic improvement programs and sustainable fisheries management.

RESULTS

Based on 278 individuals obtained from seven geographically disparate Australian locations, 10,624 high-quality SNP loci were used to characterize genetic diversity, population structure, genetic connectivity, and adaptive divergence. Significant population structure and differentiation were revealed among wild populations (average F_ST = 0.001-0.107; p <  0.05). Eighty-nine putatively outlier SNPs were identified to be potentially associated with environmental variables by using both population differentiation (BayeScan and PCAdapt) and environmental association (redundancy analysis and latent factor mixed model) analysis methods. Clear population structure with similar spatial patterns were observed in both neutral and outlier markers with three genetically distinct groups identified (north Queensland, Northern Territory, and Western Australia). Redundancy, partial redundancy, and multiple regression on distance matrices analyses revealed that both geographical distance and environmental factors interact to generate the structure observed across Australian P. monodon populations.

CONCLUSION

This study provides new insights on genetic population structure of Australian P. monodon in the face of environmental changes, which can be used to advance sustainable fisheries management and aquaculture breeding programs.

Microbiome diversity and composition varies across body areas in a freshwater turtle

There is increasing recognition that microbiomes are important for host health and ecology, and understanding host microbiomes is important for planning appropriate conservation strategies. However, microbiome data are lacking for many taxa, including turtles. To further our understanding of the interactions between aquatic microbiomes and their hosts, we used next generation sequencing technology to examine the microbiomes of the Krefft's river turtle (Emydura macquarii krefftii). We examined the microbiomes of the buccal (oral) cavity, skin on the head, parts of the shell with macroalgae and parts of the shell without macroalgae. Bacteria in the phyla Proteobacteria and Bacteroidetes were the most common in most samples (particularly buccal samples), but Cyanobacteria, Deinococcus-thermus and Chloroflexi were also common (particularly in external microbiomes). We found significant differences in community composition among each body area, as well as significant differences among individuals. The buccal cavity had lower bacterial richness and evenness than any of the external microbiomes, and it had many amplicon sequence variants (ASVs) with a low relative abundance compared to other body areas. Nevertheless, the buccal cavity also had the most unique ASVs. Parts of the shell with and without algae also had different microbiomes, with particularly obvious differences in the relative abundances of the families Methylomonaceae, Saprospiraceae and Nostocaceae. This study provides novel, baseline information about the external microbiomes of turtles and is a first step in understanding their ecological roles.

The return of the frogs: The importance of habitat refugia in maintaining diversity during a disease outbreak

Recent decades have seen the emergence and spread of numerous infectious diseases, often with severe negative consequences for wildlife populations. Nevertheless, many populations survive the initial outbreaks, and even undergo recoveries. Unfortunately, the long-term effects of these outbreaks on host population genetics are poorly understood; to increase this understanding, we examined the population genetics of two species of rainforest frogs (Litoria nannotis and Litoria serrata) that have largely recovered from a chytridiomycosis outbreak at two national parks in the Wet Tropics of northern Australia. At the wetter, northern park there was little evidence of decreased genetic diversity in either species, and all of the sampled sites had high minor allele frequencies (mean MAF = 0.230-0.235), high heterozygosity (0.318-0.325), and few monomorphic markers (1.4%-4.0%); however, some recovered L. nannotis populations had low N_e values (59.3-683.8) compared to populations that did not decline during the outbreak (1,537.4-1,756.5). At the drier, southern park, both species exhibited lower diversity (mean MAF = 0.084-0.180; heterozygosity = 0.126-0.257; monomorphic markers = 3.7%-43.5%; N_e  = 18.4-676.1). The diversity patterns in this park matched habitat patterns, with both species having higher diversity levels and fewer closely related individuals at sites with higher quality habitat. These patterns were more pronounced for L. nannotis, which has lower dispersal rates than L. serrata. These results suggest that refugia with high quality habitat are important for retaining genetic diversity during disease outbreaks, and that gene flow following disease outbreaks is important for re-establishing diversity in populations where it was reduced.

Genomic Selection in Aquaculture: Application, Limitations and Opportunities With Special Reference to Marine Shrimp and Pearl Oysters

Within aquaculture industries, selection based on genomic information (genomic selection) has the profound potential to change genetic improvement programs and production systems. Genomic selection exploits the use of realized genomic relationships among individuals and information from genome-wide markers in close linkage disequilibrium with genes of biological and economic importance. We discuss the technical advances, practical requirements, and commercial applications that have made genomic selection feasible in a range of aquaculture industries, with a particular focus on molluscs (pearl oysters, Pinctada maxima) and marine shrimp (Litopenaeus vannamei and Penaeus monodon). The use of low-cost genome sequencing has enabled cost-effective genotyping on a large scale and is of particular value for species without a reference genome or access to commercial genotyping arrays. We highlight the pitfalls and offer the solutions to the genotyping by sequencing approach and the building of appropriate genetic resources to undertake genomic selection from first-hand experience. We describe the potential to capture large-scale commercial phenotypes based on image analysis and artificial intelligence through machine learning, as inputs for calculation of genomic breeding values. The application of genomic selection over traditional aquatic breeding programs offers significant advantages through being able to accurately predict complex polygenic traits including disease resistance; increasing rates of genetic gain; minimizing inbreeding; and negating potential limiting effects of genotype by environment interactions. Further practical advantages of genomic selection through the use of large-scale communal mating and rearing systems are highlighted, as well as presenting rate-limiting steps that impact on attaining maximum benefits from adopting genomic selection. Genomic selection is now at the tipping point where commercial applications can be readily adopted and offer significant short- and long-term solutions to sustainable and profitable aquaculture industries.

Dwarf minke whales from the South Pacific share a matrilineal lineage distinct from Balaenoptera acutorostrata acutorostrata and B. a. scammoni

Dwarf minke whales are regarded as an undescribed subspecies of common minke whales (Balaenoptera acutorostrata), but appropriate conservation action requires taxonomic confirmation. The relationship of the Australian Great Barrier Reef (GBR) dwarf minke whale aggregation to other minke whales is unknown. This study aims to clarify the phylogenetic relationship of GBR dwarf minke whales, using partial mitochondrial DNA sequence data from 23 GBR dwarf minke whales, compared with other available minke whale sequences. GBR dwarf minke whales share haplotypes with other West South Pacific (WSP) dwarf minke whales. Satellite tagging studies corroborate these findings, indicating that GBR dwarf minke whales migrate south along the east Australian coast towards the Southern Ocean. Despite nuclear data not being available, GBR and WSP dwarf minke whales share a distinctive mitochondrial lineage compared with other common minke whales and should be managed independently of North Pacific and Atlantic Ocean common minke whale populations.

Genome-wide SNP analyses reveal high gene flow and signatures of local adaptation among the scalloped spiny lobster (Panulirus homarus) along the Omani coastline

BACKGROUND

The scalloped spiny lobster (Panulirus homarus) is a popular seafood commodity worldwide and an important export item from Oman. Annual catches in commercial fisheries are in serious decline, which has resulted in calls for the development of an integrated stock management approach. In Oman, the scalloped spiny lobster is currently treated as a single management unit (MU) or stock and there is an absence of information on the genetic population structure of the species that can inform management decisions, particularly at a fine-scale level. This work is the first to identify genome-wide single nucleotide polymorphisms (SNPs) for P. homarus using Diversity Arrays Technology sequencing (DArT-seq) and to elucidate any stock structure in the species.

RESULTS

After stringent filtering, 7988 high utility SNPs were discovered and used to assess the genetic diversity, connectivity and structure of P. homarus populations from Al Ashkharah, Masirah Island, Duqm, Ras Madrakah, Haitam, Ashuwaymiyah, Mirbat and Dhalkut landing sites. Pairwise FST estimates revealed low differentiation among populations (pairwise FST range = - 0.0008 - 0.0021). Analysis of genetic variation using putatively directional FST outliers (504 SNPs) revealed higher and significant pairwise differentiation (p < 0.01) for all locations, with Ashuwaymiyah being the most diverged population (Ashuwaymiyah pairwise FST range = 0.0288-0.0736). Analysis of population structure using Discriminant Analysis of Principal Components (DAPC) revealed a broad admixture among P. homarus, however, Ashuwaymiyah stock appeared to be potentially under local adaptive pressures. Fine scale analysis using Netview R provided further support for the general admixture of P. homarus.

CONCLUSIONS

Findings here suggested that stocks of P. homarus along the Omani coastline are admixed. Yet, fishery managers need to treat the lobster stock from Ashuwaymiyah with caution as it might be subject to local adaptive pressures. We emphasize further study with larger number of samples to confirm the genetic status of the Ashuwaymiyah stock. The approach utilised in this study has high transferability in conservation and management of other marine stocks with similar biological and ecological attributes.

Sex-specific dmrt1 and cyp19a1 methylation and alternative splicing in gonads of the protandrous hermaphrodite barramundi

Epigenetics is involved in sex differentiation of gonochoristic and hermaphroditic fish species, whereby two genes dmrt1 (pro-male) and cyp19a1 (pro-female) are known to play major roles. Barramundi, Lates calcarifer, is an important tropical aquaculture species that undergo natural and permanent male to female sex change, a process for which the exact underlying molecular mechanisms are still unknown. To elucidate whether DNA methylation is involved in sex control of barramundi, a next-generation bisulfite amplicon sequencing approach was used to target 146 CpG sites within proximal promoters and first exons of seven sex-related genes (dmrt1, cyp19a1, amh, foxl2, nr5a2, sox8 and sox9) of 24 testis and 18 ovaries of captive and wild adult barramundi. Moreover, comparative expression profiles of the key dmrt1 and cyp19a1 genes were further investigated using RT-qPCR and Sanger sequencing approaches, whereas expression levels of remaining targeted genes were based on available literature for the species. Results showed that cyp19a1 and amh were more methylated in males, whereas dmrt1 and nr5a2 were more methylated in females (P < 0.001), with no gender differences found for foxl2, sox8 or sox9 genes (P > 0.05). Sex-biased promoter DNA methylation was inversely related to gene expression only for dmrt1 and nr5a2, and directly related to amh expression, whereas no differences in cyp19a1 expression were found between testes and ovaries. Notably, unique sex-specific alternative splicing of dmrt1 and cyp19a1 were discovered, whereby males lacked the full-length aromatase coding cyp19a1 mRNA due to partial or total exon splicing, and females lacked the dmrt1 exon containing the DM-domain sequence. This study advances the current knowledge aiming to elucidate the genetic mechanisms within male and female gonads of this large protandrous hermaphrodite by providing the first evidence of epigenetics and alternative splicing simultaneously affecting key genes (cyp19a1 and dmrt1) central to sex differentiation pathways.

Genomic comparisons reveal biogeographic and anthropogenic impacts in the koala (Phascolarctos cinereus): a dietary-specialist species distributed across heterogeneous environments

The Australian koala is an iconic marsupial with highly specific dietary requirements distributed across heterogeneous environments, over a large geographic range. The distribution and genetic structure of koala populations has been heavily influenced by human actions, specifically habitat modification, hunting and translocation of koalas. There is currently limited information on population diversity and gene flow at a species-wide scale, or with consideration to the potential impacts of local adaptation. Using species-wide sampling across heterogeneous environments, and high-density genome-wide markers (SNPs and PAVs), we show that most koala populations display levels of diversity comparable to other outbred species, except for those populations impacted by population reductions. Genetic clustering analysis and phylogenetic reconstruction reveals a lack of support for current taxonomic classification of three koala subspecies, with only a single evolutionary significant unit supported. Furthermore, ~70% of genetic variance is accounted for at the individual level. The Sydney Basin region is highlighted as a unique reservoir of genetic diversity, having higher diversity levels (i.e., Blue Mountains region; AvHe^corr=0.20, PL% = 68.6). Broad-scale population differentiation is primarily driven by an isolation by distance genetic structure model (49% of genetic variance), with clinal local adaptation corresponding to habitat bioregions. Signatures of selection were detected between bioregions, with no single region returning evidence of strong selection. The results of this study show that although the koala is widely considered to be a dietary-specialist species, this apparent specialisation has not limited the koala’s ability to maintain gene flow and adapt across divergent environments as long as the required food source is available.

The State of "Omics" Research for Farmed Penaeids: Advances in Research and Impediments to Industry Utilization

Elucidating the underlying genetic drivers of production traits in agricultural and aquaculture species is critical to efforts to maximize farming efficiency. "Omics" based methods (i.e., transcriptomics, genomics, proteomics, and metabolomics) are increasingly being applied to gain unprecedented insight into the biology of many aquaculture species. While the culture of penaeid shrimp has increased markedly, the industry continues to be impeded in many regards by disease, reproductive dysfunction, and a poor understanding of production traits. Extensive effort has been, and continues to be, applied to develop critical genomic resources for many commercially important penaeids. However, the industry application of these genomic resources, and the translation of the knowledge derived from "omics" studies has not yet been completely realized. Integration between the multiple "omics" resources now available (i.e., genome assemblies, transcriptomes, linkage maps, optical maps, and proteomes) will prove critical to unlocking the full utility of these otherwise independently developed and isolated resources. Furthermore, emerging "omics" based techniques are now available to address longstanding issues with completing keystone genome assemblies (e.g., through long-read sequencing), and can provide cost-effective industrial scale genotyping tools (e.g., through low density SNP chips and genotype-by-sequencing) to undertake advanced selective breeding programs (i.e., genomic selection) and powerful genome-wide association studies. In particular, this review highlights the status, utility and suggested path forward for continued development, and improved use of "omics" resources in penaeid aquaculture.

Gene flow and genetic structure in Nile perch, Lates niloticus, from African freshwater rivers and lakes

BACKGROUND

Geological evolution of the African continent has been subject to complex processes including uplift, volcanism, desert formation and tectonic rifting. This complex geology has created substantial biogeographical barriers, and coupled with anthropogenic introductions of freshwater fishes, has influenced the genetic diversity, connectivity and sub-structuring of the teleost fauna. Nile perch, Lates niloticus, is an iconic fish in Africa and is of high commercial importance, both in the species' native range and where it has been translocated. However, the species is in decline and there is a need to understand its population genetic structure to facilitate sustainable management of the fishery and aquaculture development.

METHODOLOGY

Nile perch tissue samples were acquired from two West and four East (Lakes; Albert, Kyoga, Victoria and Turkana) African locations. Nineteen polymorphic microsatellite loci were used to study the genetic variation among populations across regions (West and East Africa), as well as between native and introduced environments within East Africa.

PRINCIPAL FINDINGS AND THEIR SIGNIFICANCE

Results revealed strong and significant genetic structuring among populations across the sampled distribution (divergence across regions, FCT = 0.26, P = 0.000). STRUCTURE analysis at a broad scale revealed K = 2 clusters, the West African individuals were assigned to one cluster, while all individuals from the East African region, regardless of whether native or introduced, were assigned to another cluster. The distinct genetic clusters identified in the current study between the West and East African Nile perch, appear to have been maintained by presence of biogeographic barriers and restricted gene flow between the two regions. Therefore, any translocations of Nile perch should be carefully considered across the regions of West and East Africa. Further analysis at a regional scale revealed further structuring of up to K = 3 genetic clusters in East African Nile perch. Significantly (P < 0.05) lower genetic diversity based on analysis of allelic richness (AR) was obtained for the two translocated populations of Lake Kyoga (AR = 3.61) and Lake Victoria (AR = 3.52), compared to Nile perch populations from their putative origins of Lakes Albert (AR = 4.12) and Turkana (AR = 4.43). The lower genetic diversity in the translocated populations may be an indication of previous bottlenecks and may also indicate a difficulty for these populations to persist and adapt to climatic changes and anthropogenic pressures that are currently present in the East African region.

Genome-wide comparisons reveal a clinal species pattern within a holobenthic octopod-the Australian Southern blue-ringed octopus, Hapalochlaena maculosa (Cephalopoda: Octopodidae)

The southern blue-ringed octopus, Hapalochlaena maculosa (Hoyle, 1883) lacks a planktonic dispersal phase, yet ranges across Australia's southern coastline. This species' brief and holobenthic life history suggests gene flow might be limited, leaving distant populations prone to strong genetic divergence. This study used 17,523 genome-wide SNP loci to investigate genetic structuring and local adaptation patterns of H. maculosa among eight sampling sites along its reported range. Within sites, interrelatedness was very high, consistent with the limited dispersal of this taxon. However, inbreeding coefficients were proportionally lower among sites where substructuring was not detected, suggesting H. maculosa might possess a mechanism for inbreeding avoidance. Genetic divergence was extremely high among all sites, with the greatest divergence observed between both ends of the distribution, Fremantle, WA, and Stanley, TAS. Genetic distances closely followed an isolation by geographic distance pattern. Outlier analyses revealed distinct selection signatures at all sites, with the strongest divergence reported between Fremantle and the other Western Australian sites. Phylogenetic reconstructions using the described sister taxon H. fasciata (Hoyle, 1886) further supported that the genetic divergence between distal H. maculosa sites in this study was equivalent to that of between established heterospecifics within this genus. However, it is advocated that taxonomic delineations within this species should be made with caution. These data indicate that H. maculosa forms a clinal species pattern across its geographic range, with gene flow present through allele sharing between adjacent populations. Morphological investigations are recommended for a robust resolution of the taxonomic identity and ecotype boundaries of this species.

Extensive genetic differentiation detected within a model marsupial, the tammar wallaby (Notamacropus eugenii)

The tammar wallaby (Notamacropus eugenii) is one of the most intensively studied of all macropodids and was the first Australasian marsupial to have its genome sequenced. However, comparatively little is known about genetic diversity and differentiation amongst the morphologically distinct allopatric populations of tammar wallabies found in Western (WA) and South Australia (SA). Here we compare autosomal and Y-linked microsatellite genotypes, as well as sequence data (~600 bp) from the mitochondrial DNA (mtDNA) control region (CR) in tammar wallabies from across its distribution. Levels of diversity at autosomal microsatellite loci were typically high in the WA mainland and Kangaroo Island (SA) populations (A = 8.9-10.6; He = 0.77-0.78) but significantly reduced in other endemic island populations (A = 3.8-4.1; He = 0.41-0.48). Autosomal and Y-linked microsatellite loci revealed a pattern of significant differentiation amongst populations, especially between SA and WA. The Kangaroo Island and introduced New Zealand population showed limited differentiation. Multiple divergent mtDNA CR haplotypes were identified within both SA and WA populations. The CR haplotypes of tammar wallabies from SA and WA show reciprocal monophyly and are highly divergent (14.5%), with levels of sequence divergence more typical of different species. Within WA tammar wallabies, island populations each have unique clusters of highly related CR haplotypes and each is most closely related to different WA mainland haplotypes. Y-linked microsatellite haplotypes show a similar pattern of divergence although levels of diversity are lower. In light of these differences, we suggest that two subspecies of tammar wallaby be recognized; Notamacropus eugenii eugenii in SA and N. eugenii derbianus in WA. The extensive neutral genetic diversity and inter-population differentiation identified within tammar wallabies should further increase the species value and usefulness as a model organism.

Swept away: ocean currents and seascape features influence genetic structure across the 18,000 Km Indo-Pacific distribution of a marine invertebrate, the black-lip pearl oyster Pinctada margaritifera

BACKGROUND

Genetic structure in many widely-distributed broadcast spawning marine invertebrates remains poorly understood, posing substantial challenges for their fishery management, conservation and aquaculture. Under the Core-Periphery Hypothesis (CPH), genetic diversity is expected to be highest at the centre of a species' distribution, progressively decreasing with increased differentiation towards outer range limits, as populations become increasingly isolated, fragmented and locally adapted. The unique life history characteristics of many marine invertebrates such as high dispersal rates, stochastic survival and variable recruitment are also likely to influence how populations are organised. To examine the microevolutionary forces influencing population structure, connectivity and adaptive variation in a highly-dispersive bivalve, populations of the black-lip pearl oyster Pinctada margaritifera were examined across its ~18,000 km Indo-Pacific distribution.

RESULTS

Analyses utilising 9,624 genome-wide SNPs and 580 oysters, discovered differing patterns of significant and substantial broad-scale genetic structure between the Indian and Pacific Ocean basins. Indian Ocean populations were markedly divergent (F st = 0.2534-0.4177, p < 0.001), compared to Pacific Ocean oysters, where basin-wide gene flow was much higher (F st = 0.0007-0.1090, p < 0.001). Partitioning of genetic diversity (hierarchical AMOVA) attributed 18.1% of variance between ocean basins, whereas greater proportions were resolved within samples and populations (45.8% and 35.7% respectively). Visualisation of population structure at selectively neutral loci resolved three and five discrete genetic clusters for the Indian and Pacific Oceans respectively. Evaluation of genetic structure at adaptive loci for Pacific populations (89 SNPs under directional selection; F st = 0.1012-0.4371, FDR = 0.05), revealed five clusters identical to those detected at neutral SNPs, suggesting environmental heterogeneity within the Pacific. Patterns of structure and connectivity were supported by Mantel tests of isolation by distance (IBD) and independent hydrodynamic particle dispersal simulations.

CONCLUSIONS

It is evident that genetic structure and connectivity across the natural range of P. margaritifera is highly complex, and produced by the interaction of ocean currents, IBD and seascape features at a broad scale, together with habitat geomorphology and local adaptation at regional levels. Overall population organisation is far more elaborate than generalised CPH predictions, however valuable insights for regional fishery management, and a greater understanding of range-wide genetic structure in a highly-dispersive marine invertebrate have been gained.

A Parallel Population Genomic and Hydrodynamic Approach to Fishery Management of Highly-Dispersive Marine Invertebrates: The Case of the Fijian Black-Lip Pearl Oyster Pinctada margaritifera

Fishery management and conservation of marine species increasingly relies on genetic data to delineate biologically relevant stock boundaries. Unfortunately for high gene flow species which may display low, but statistically significant population structure, there is no clear consensus on the level of differentiation required to resolve distinct stocks. The use of fine-scale neutral and adaptive variation, considered together with environmental data can offer additional insights to this problem. Genome-wide genetic data (4,123 SNPs), together with an independent hydrodynamic particle dispersal model were used to inform farm and fishery management in the Fijian black-lip pearl oyster Pinctada margaritifera, where comprehensive fishery management is lacking, and the sustainability of exploitation uncertain. Weak fine-scale patterns of population structure were detected, indicative of broad-scale panmixia among wild oysters, while a hatchery-sourced farmed population exhibited a higher degree of genetic divergence (Fst = 0.0850-0.102). This hatchery-produced population had also experienced a bottleneck (NeLD = 5.1; 95% C.I. = [5.1-5.3]); compared to infinite NeLD estimates for all wild oysters. Simulation of larval transport pathways confirmed the existence of broad-scale mixture by surface ocean currents, correlating well with fine-scale patterns of population structuring. Fst outlier tests failed to detect large numbers of loci supportive of selection, with 2-5 directional outlier SNPs identified (average Fst = 0.116). The lack of biologically significant population genetic structure, absence of evidence for local adaptation and larval dispersal simulation, all indicate the existence of a single genetic stock of P. margaritifera in the Fiji Islands. This approach using independent genomic and oceanographic tools has allowed fundamental insights into stock structure in this species, with transferability to other highly-dispersive marine taxa for their conservation and management.

Unravelling the effects of gene flow and selection in highly connected populations of the silver-lip pearl oyster (Pinctada maxima)

Many marine organisms often display weak levels of population genetic structuring as a result of both environmental characteristics (e.g., ocean currents) and life history traits (e.g., widely dispersed planktonic larval stages) maintaining high levels of gene flow. This can lead to the assumption that these organisms can be managed as a single stock based on high levels of population connectivity. However, this neglects to account for other micro-evolutionary forces such as selection, which also shape these populations. This study utilizes 1130 genome-wide SNP loci to unravel the effects of gene flow and selection shaping three highly connected populations of the silver-lip pearl oyster (Pinctada maxima) in the ecologically and economically important Indo-Pacific region (Aru, Bali, and West Papua). Twenty-two loci under directional selection were identified amongst the populations, providing further supporting evidence of strong local adaptation (i.e., G×E effects) among populations in this region. Global Fst values for directional outliers (0.348) were up to eight times greater than for neutral markers (0.043). Pairwise Fst comparisons between Aru and Bali revealed the largest directional differences (0.488), while Bali and West Papua had the least (0.062). Unrooted neighbour-joining (NJ) distance trees and genetic diversity indices of directional outliers revealed that individuals from Bali and West Papua had reduced allelic variation (MAFavg=0.144, Ho=0.238 and MAFavg=0.232, Ho=0.369, respectively) compared to Aru (MAFavg=0.292, Ho=0.412). This indicates that directional selection is most likely acting upon genetic variation within the Bali and West Papua populations. NJ distance trees, discriminant analysis of principal components, and Fst analyses of directional outliers revealed two divergent groups ("Bali/West Papua"; "Aru") that had previously gone unrecognized. This study not only illustrates that relatively strong local adaptive forces are occurring despite high gene flow, but identifies the populations that are most likely experiencing selection. Additionally, this study highlights the need to understand all micro-evolutionary forces acting on populations when resolving stock structure.

Fishing for divergence in a sea of connectivity: The utility of ddRADseq genotyping in a marine invertebrate, the black-lip pearl oyster Pinctada margaritifera

Population genomic investigations on highly dispersive marine organisms typically require thousands of genome-wide SNP loci to resolve fine-scale population structure and detect signatures of selection. This information is important for species conservation efforts and stock management in both wild and captive populations, as well as genome mapping and genome wide association studies. Double digest Restriction site-Associated DNA Sequencing (ddRADseq) is a recent tool for delivering genome wide SNPs for non-model organisms. However, its application to marine invertebrate taxa has been limited, particularly given the complex and highly repetitive nature of many of these organisms' genomes. This study develops and evaluates an optimised ddRADseq technique together with associated analyses for generating genome-wide SNP data, and performs population genomic analyses to inform aquaculture and fishery management of a marine bivalve, the black-lip pearl oyster Pinctada margaritifera. A total of 5243 high-quality genome-wide SNP markers were detected, and used to assess population structure, genome diversity, detect Fst outliers and perform association testing in 156 individuals belonging to three wild and one hatchery produced populations from the Fiji Islands. Shallow but significant population structure was revealed among all wild populations (average pairwise Fst=0.046) when visualised with DAPC and an individual network analysis (NetView P), with clear evidence of a genetic bottleneck in the hatchery population (NeLD=6.1), compared to wild populations (NeLD>192.5). Fst outlier detection revealed 42-62 highly differentiated SNPs (p<0.02), while case-control association discovered up to 152 SNPs (p<0.001). Both analyses were able to successfully differentiate individuals between the orange and black tissue colour morphotypes characteristic of this species. BLAST searches revealed that five of these SNPs were associated with a melanin biosynthesis pathway, demonstrating their biological relevance. This study has produced highly informative SNP and population genomic data in P. margaritifera, and using the same approach promises to be of substantial value to a range of other non-model, broadcast-spawning or marine invertebrate taxa.

netview p: a network visualization tool to unravel complex population structure using genome-wide SNPs

Network-based approaches are emerging as valuable tools for the analysis of complex genetic structure in wild and captive populations. netview p combines data quality control with the construction of population networks through mutual k-nearest neighbours thresholds applied to genome-wide SNPs. The program is cross-platform compatible, open-source and efficiently operates on data ranging from hundreds to hundreds of thousands of SNPs. The pipeline was used for the analysis of pedigree data from simulated (n = 750, SNPs = 1279) and captive silver-lipped pearl oysters (n = 415, SNPs = 1107), wild populations of the European hake from the Atlantic and Mediterranean (n = 834, SNPs = 380) and grey wolves from North America (n = 239, SNPs = 78 255). The population networks effectively visualize large- and fine-scale genetic structure within and between populations, including family-level structure and relationships. netview p comprises a network-based addition to other population analysis tools and provides user-friendly access to a complex network analysis pipeline through implementation in python.

Nocturnal mating behaviour and dynamic male investment of copulation time in the southern blue-ringed octopus, Hapalochlaena maculosa (Cephalopoda: Octopodidae)

The southern blue-ringed octopus,Hapalochlaena maculosaHoyle (1883), is a nocturnal species that exhibits a mating system in which females hold sperm from multiple males over a one to two month breeding window before laying a single egg clutch. Contrary to most studied animal mating systems where anisogamy exists, gamete package production is limited for both males and females of this species (approx. 50 spermatophores/eggs). This presents an animal model for studying aspects of sperm competition and dynamic mate choice behaviours. The present study reports on the mating behaviour ofH. maculosaobserved under laboratory conditions using infrared closed-circuit television video footage. Rates of male copulation attempts increased with male size, while female receptivity to mating attempts increased with female size, resulting in larger animals of both sexes gaining more copulations and spending more time per day in copulation. There was some evidence of female preference of larger males, but no male preference of females based on measured morphological traits. Both sexes terminated copulations in equal frequencies but male-terminated copulations were significantly shorter in duration. Males were more likely to terminate copulation early with females they had previously mated with, however were less likely to do so if the female had recently mated with a different male. Among male-terminated copulations, males mated for longer with females that had previously mated with other males in the trial. Male–male mounts were as common as male–female mounts, suggesting that maleH. maculosaare not able to discriminate the sex of conspecifics. These findings suggest male strategic allocation of spermatophores based female mating history is an important factor influencing mating behaviours of this species.

Effects of growth phase and nitrogen starvation on expression of fatty acid desaturases and fatty acid composition of Isochrysis aff. galbana (TISO)

Very long-chain polyunsaturated fatty acids (VLC-PUFAs) are important dietary requirements for maintaining human health. Many marine microalgae are naturally high in ω-3 VLC-PUFAs, however, the molecular mechanisms underpinning fatty acid (FA) desaturation and elongation in algae are poorly understood. An advanced molecular understanding would facilitate improvements of this nascent industry. We aimed to investigate expression responses of four front-end fatty acid desaturase genes and downstream effects on FA profiles to nitrogen limitation and cultivation growth stage in Isochrysis aff. galbana (TISO). Cultures were grown in nitrogen-replete and -deplete medium; samples were harvested during logarithmic, late logarithmic and stationary growth phases to analyse FA content/composition and gene expression of ∆(6)-, ∆(8)-, ∆(5)- and ∆(4)-desaturases (d6FAD (putative), d8FAD, d5FAD and d4FAD, respectively). d6FAD (putative) exhibited no differential expression, while d8FAD, d5FAD and d4FAD were significantly upregulated during logarithmic growth of nutrient-replete cultures, coinciding with rapid cell division. In conclusion, it is demonstrated that expression of some FADs in I. aff. galbana varies with culture age and nitrogen status which has downstream consequences on FA desaturation levels. This has implications for the commercial production of VLC-PUFAs where a trade-off between total lipid yield and VLC-PUFAs has to be made.

Genome-wide SNP validation and mantle tissue transcriptome analysis in the silver-lipped pearl oyster, Pinctada maxima

Pearl oysters are not only farmed for their gemstone quality pearls worldwide, but they are also becoming important model organisms for investigating genetic mechanisms of biomineralisation. Despite their economic and scientific significance, limited genomic resources are available for this important group of bivalves, hampering investigations into identifying genes that regulate important pearl quality traits and unique biological characteristics (i.e. biomineralisation). The silver-lipped pearl oyster, Pinctada maxima, is one species where there is interest in understanding genes that regulate commercially important pearl traits, but presently, there is a dearth of genomic information. The objective of this study was to develop and validate a large number of type I genome-wide single nucleotide polymorphisms (SNPs) for P. maxima suitable for high-throughput genotyping. In addition, sequence annotations and Gene Ontology terms were assigned to a large mantle tissue 454 expressed sequence tag assembly (96,794 contigs) and information on known bivalve biomineralisation genes was incorporated into SNP discovery. The SNP discovery effort resulted in the de novo identification of 172,625 SNPs, of which 9,108 were identified as high value [minor allele frequency (MAF)≥ 0.15, read depth  ≥ 8]. Validation of 2,782 of these SNPs using Illumina iSelect Infinium genotyping technology returned some of the highest assay conversion (86.6 %) and validation (59.9 %; mean MAF 0.28) rates observed in aquaculture species to date. Genomic resources presented here will be pivotal to future research investigating the biological mechanisms behind biomineralisation and will form a strong foundation for genetic selective breeding programs in the P. maxima pearling industry.

A high-density SNP genetic linkage map for the silver-lipped pearl oyster, Pinctada maxima: a valuable resource for gene localisation and marker-assisted selection

Background

The silver-lipped pearl oyster, Pinctada maxima, is an important tropical aquaculture species extensively farmed for the highly sought "South Sea" pearls. Traditional breeding programs have been initiated for this species in order to select for improved pearl quality, but many economic traits under selection are complex, polygenic and confounded with environmental factors, limiting the accuracy of selection. The incorporation of a marker-assisted selection (MAS) breeding approach would greatly benefit pearl breeding programs by allowing the direct selection of genes responsible for pearl quality. However, before MAS can be incorporated, substantial genomic resources such as genetic linkage maps need to be generated. The construction of a high-density genetic linkage map for P. maxima is not only essential for unravelling the genomic architecture of complex pearl quality traits, but also provides indispensable information on the genome structure of pearl oysters.

Results

A total of 1,189 informative genome-wide single nucleotide polymorphisms (SNPs) were incorporated into linkage map construction. The final linkage map consisted of 887 SNPs in 14 linkage groups, spans a total genetic distance of 831.7 centimorgans (cM), and covers an estimated 96% of the P. maxima genome. Assessment of sex-specific recombination across all linkage groups revealed limited overall heterochiasmy between the sexes (i.e. 1.15:1 F/M map length ratio). However, there were pronounced localised differences throughout the linkage groups, whereby male recombination was suppressed near the centromeres compared to female recombination, but inflated towards telomeric regions. Mean values of LD for adjacent SNP pairs suggest that a higher density of markers will be required for powerful genome-wide association studies. Finally, numerous nacre biomineralization genes were localised providing novel positional information for these genes.

Conclusions

This high-density SNP genetic map is the first comprehensive linkage map for any pearl oyster species. It provides an essential genomic tool facilitating studies investigating the genomic architecture of complex trait variation and identifying quantitative trait loci for economically important traits useful in genetic selection programs within the P. maxima pearling industry. Furthermore, this map provides a foundation for further research aiming to improve our understanding of the dynamic process of biomineralization, and pearl oyster evolution and synteny.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-810) contains supplementary material, which is available to authorized users.

Next-generation transcriptome profiling reveals insights into genetic factors contributing to growth differences and temperature adaptation in Australian populations of barramundi (Lates calcarifer)

Identification of genetically-regulated adaptation in fish is a precursor to understanding how populations will respond to future climate induced stressors like temperature. Australian populations of barramundi (Lates calcarifer) show strong evidence of local adaptation to temperature. However, the phenotypic consequences of this adaptation are unknown and the genetic mechanisms underlying this response are poorly understood. In the current study, two populations of barramundi from temperature extremes of the species Australian distribution were communally reared at cool (22°C), control (28°C) and hot (36°C) water temperatures for 3.5months. Southern populations of barramundi originating from a cooler environment grew significantly faster at 22°C than northern populations of warm adapted barramundi. However, no difference in population growth was present at either 28°C or 36°C. The underlying transcriptome profile of barramundi was examined via Illumina mRNA deep sequencing to determine the major contributing gene categories giving rise to phenotypic differences in barramundi population growth. Gene ontology (GO) analysis revealed enrichment in categories relating to the regulation of peptidase activity as well as microtubule, cytoplasmic and cellular metabolic based processes. Further analysis of the GO category "microtubule based process" with associated genes from the "response to stress" category revealed an apparent re-organization of cytoskeletal elements in response to an induced cold stress in northern barramundi reared at 22°C, when compared with northern barramundi reared at 36°C. Between southern barramundi and northern barramundi reared at 36°C, an analysis of the "endopeptidase inhibitor activity" GO category in conjunction with stress genes indicated a suppression of the complement system in southern barramundi along with an increase in the cellular stress response. The results of the present study show that southern populations of barramundi exhibit underlying molecular adaptation to cooler water temperatures, but still retain a tolerance for warm water temperatures. Furthermore, GO profiling has revealed groups of genes that underlie population differences in temperature tolerance as a means to prioritize the analysis of differential gene expression in studies of local adaptation in the future.

Paternally inherited genetic markers reveal new insights into genetic structuring within Macropus fuliginosus and hybridisation with sympatric Macropus giganteus

There are several aspects of biology in which the contribution of males and females is unequal. In these instances the examination of Y chromosome markers may be used to elucidate male-specific attributes. Here, male dispersal patterns and genetic structuring were examined using four Y-microsatellite loci in 186 male western grey kangaroos, Macropus fuliginosus, from throughout the species’ trans-continental distribution. In addition, 52 male grey kangaroos were examined to investigate hybridisation between M. fuliginosus and the eastern grey kangaroo, Macropus giganteus, in their region of sympatry in eastern Australia. Detected Y chromosome diversity was low, resulting from low effective male population size due to skewed sex ratios and a polygynous mating system. As expected, male dispersal was high across the range. However, the Lake Torrens–Flinders Ranges region appears to have significantly restricted male movement between eastern and central/western Australia. There was little evidence to suggest that other barriers (Nullarbor Plain and Swan River Valley) previously identified by nuclear and mitochondrial DNA marker studies restrict male movement. Hence, the admixture events previously identified may be associated with high male dispersal. Within the region of sympatry between M. fuliginosus and M. giganteus in eastern Australia, four M. giganteus individuals were found to possess M. fuliginosus Y-haplotypes. These results confirm the occurrence of hybridisation between male M. fuliginosus and female M. giganteus. Additionally, the introgression of M. fuliginosus Y-haplotypes into M. giganteus populations indicates that at least some male hybrids are fertile, despite evidence to the contrary from captive studies. This study has provided insights into the male contribution to population history, structure and hybridisation in M. fuliginosus, which were not predicted by comparisons between biparentally and maternally inherited markers. This highlights the importance of direct examination of the Y chromosome to provide novel insights into male-mediated processes, especially where the contribution of the sexes may differ.

Genetic analysis reveals a distinct and highly diverse koala (Phascolarctos cinereus) population in South Gippsland, Victoria, Australia

Population genetics can reveal otherwise hidden information involving a species’ history in a given region. Koalas were thought to have been virtually exterminated from the Australian state of Victoria during the koala fur trade of the late 1800s. Koalas in the South Gippsland region of Victoria were examined using microsatellite markers to infer population structure and gene flow and to locate a possible remnant gene pool. The results indicate that the South Gippsland koala population had higher genetic diversity (A = 5.97, HO = 0.564) than other published Victorian populations, and was genetically distinct from other koala populations examined. South Gippsland koalas, therefore, may have survived the population reductions of the koala fur trade and now represent a remnant Victorian gene pool that has been largely lost from the remainder of Victoria. This paper illustrates that historic anthropogenic impacts have had little effect on reducing the genetic diversity of a population in the South Gippsland region. However, the South Gippsland population is now subject to threats such as logging and loss of habitat from housing and agriculture expansion. Our results suggest that the South Gippsland koalas require an alternative conservation management program.

Genome sequence of an Australian kangaroo, Macropus eugenii, provides insight into the evolution of mammalian reproduction and development

BACKGROUND

We present the genome sequence of the tammar wallaby, Macropus eugenii, which is a member of the kangaroo family and the first representative of the iconic hopping mammals that symbolize Australia to be sequenced. The tammar has many unusual biological characteristics, including the longest period of embryonic diapause of any mammal, extremely synchronized seasonal breeding and prolonged and sophisticated lactation within a well-defined pouch. Like other marsupials, it gives birth to highly altricial young, and has a small number of very large chromosomes, making it a valuable model for genomics, reproduction and development.

RESULTS

The genome has been sequenced to 2 × coverage using Sanger sequencing, enhanced with additional next generation sequencing and the integration of extensive physical and linkage maps to build the genome assembly. We also sequenced the tammar transcriptome across many tissues and developmental time points. Our analyses of these data shed light on mammalian reproduction, development and genome evolution: there is innovation in reproductive and lactational genes, rapid evolution of germ cell genes, and incomplete, locus-specific X inactivation. We also observe novel retrotransposons and a highly rearranged major histocompatibility complex, with many class I genes located outside the complex. Novel microRNAs in the tammar HOX clusters uncover new potential mammalian HOX regulatory elements.

CONCLUSIONS

Analyses of these resources enhance our understanding of marsupial gene evolution, identify marsupial-specific conserved non-coding elements and critical genes across a range of biological systems, including reproduction, development and immunity, and provide new insight into marsupial and mammalian biology and genome evolution.

A first-generation integrated tammar wallaby map and its use in creating a tammar wallaby first-generation virtual genome map

BACKGROUND

The limited (2X) coverage of the tammar wallaby (Macropus eugenii) genome sequence dataset currently presents a challenge for assembly and anchoring onto chromosomes. To provide a framework for this assembly, it would be a great advantage to have a dense map of the tammar wallaby genome. However, only limited mapping data are available for this non-model species, comprising a physical map and a linkage map.

RESULTS

We combined all available tammar wallaby mapping data to create a tammar wallaby integrated map, using the Location DataBase (LDB) strategy. This first-generation integrated map combines all available information from the second-generation tammar wallaby linkage map with 148 loci, and extensive FISH mapping data for 492 loci, especially for genes likely to be located at the ends of wallaby chromosomes or at evolutionary breakpoints inferred from comparative information. For loci whose positions are only approximately known, their location in the integrated map was refined on the basis of comparative information from opossum (Monodelphis domestica) and human. Interpolation of segments from the opossum and human assemblies into the integrated map enabled the subsequent construction of a tammar wallaby first-generation virtual genome map, which comprises 14336 markers, including 13783 genes recruited from opossum and human assemblies. Both maps are freely available at http://compldb.angis.org.au.

CONCLUSIONS

The first-generation integrated map and the first-generation virtual genome map provide a backbone for the chromosome assembly of the tammar wallaby genome sequence. For example, 78% of the 10257 gene-scaffolds in the Ensembl annotation of the tammar wallaby genome sequence (including 10522 protein-coding genes) can now be given a chromosome location in the tammar wallaby virtual genome map.

A second-generation anchored genetic linkage map of the tammar wallaby (Macropus eugenii)

Background

The tammar wallaby, Macropus eugenii, a small kangaroo used for decades for studies of reproduction and metabolism, is the model Australian marsupial for genome sequencing and genetic investigations. The production of a more comprehensive cytogenetically-anchored genetic linkage map will significantly contribute to the deciphering of the tammar wallaby genome. It has great value as a resource to identify novel genes and for comparative studies, and is vital for the ongoing genome sequence assembly and gene ordering in this species.

Results

A second-generation anchored tammar wallaby genetic linkage map has been constructed based on a total of 148 loci. The linkage map contains the original 64 loci included in the first-generation map, plus an additional 84 microsatellite loci that were chosen specifically to increase coverage and assist with the anchoring and orientation of linkage groups to chromosomes. These additional loci were derived from (a) sequenced BAC clones that had been previously mapped to tammar wallaby chromosomes by fluorescence in situ hybridization (FISH), (b) End sequence from BACs subsequently FISH-mapped to tammar wallaby chromosomes, and (c) tammar wallaby genes orthologous to opossum genes predicted to fill gaps in the tammar wallaby linkage map as well as three X-linked markers from a published study. Based on these 148 loci, eight linkage groups were formed. These linkage groups were assigned (via FISH-mapped markers) to all seven autosomes and the X chromosome. The sex-pooled map size is 1402.4 cM, which is estimated to provide 82.6% total coverage of the genome, with an average interval distance of 10.9 cM between adjacent markers. The overall ratio of female/male map length is 0.84, which is comparable to the ratio of 0.78 obtained for the first-generation map.

Conclusions

Construction of this second-generation genetic linkage map is a significant step towards complete coverage of the tammar wallaby genome and considerably extends that of the first-generation map. It will be a valuable resource for ongoing tammar wallaby genetic research and assembling the genome sequence. The sex-pooled map is available online at http://compldb.angis.org.au/.

In silico whole-genome EST analysis reveals 2322 novel microsatellites for the silver-lipped pearl oyster, Pinctada maxima

Molecular stock improvement techniques such as marker assisted selection have great potential in accelerating selective breeding programmes for animal production industries. However, the discovery and application of trait/marker associations usually requires a large number of genome-wide polymorphic loci. Here, we present 2322 unique microsatellites for the silver-lipped pearl oyster, Pinctada maxima, a species of aquaculture importance throughout the Indo-Australian Archipelago for production of the highly valued South Sea pearl. More than 1.2 million Roche 454 expressed sequence tag (EST) reads were screened for microsatellite repeat motifs. A total of 12,604 sequences contained either a di, tri, tetra, penta or hexa microsatellite repeat motif (n ≥ 6), with 6435 of these sequences having sufficient flanking regions for primer development. All identified microsatellites with designed primers were condensed into 2322 unique clusters (i.e., unique loci) of which 360 were shown to be polymorphic based on multiple sequence reads with different repeat motifs. Genotyping of five microsatellite loci demonstrated that in silico evaluation of polymorphism levels was a very useful method for identification of polymorphic loci, with the variation uncovered being a lower bound. Gene Ontology annotations of sequences containing microsatellites suggest that most are derived from a diverse array of unique genes. This EST derived microsatellite database will be a valuable resource for future studies in genetic map construction, diversity analysis, quantitative trait loci analysis, association mapping and marker assisted selection, not only for P. maxima, but also closely related species within the genus Pinctada.

Defining spatial genetic structure and management units for vulnerable koala (Phascolarctos cinereus) populations in the Sydney region, Australia

Context. Mammal populations around the world are increasingly threatened with population fragmentation because of loss of habitat or barriers to gene flow. The investigation of koala populations in the Sydney region not only provides valuable information about this vulnerable species, but also serves as a model for other species that have suffered major rapid declines in population size, and are now recovering in fragmented habitat. The peri-urban study region allows investigation of the impact of landscape features such as major roads and housing developments on koala gene flow. Aims. Animals originating from four geographic sampling areas around Sydney, New South Wales, Australia, were examined to determine population structure and gene flow and to identify barriers to gene flow and management units. Methods. The present study examined 12 microsatellite loci and used Bayesian assignment methods and genic frequency analysis methods to identify demographically separate populations and barriers to gene flow between those populations. Key results. Three discrete populations were resolved, with all displaying moderate to high levels of genetic differentiation among them (θ = 0.141–0.224). The allelic richness and heterozygosity of the Blue Mountains population (A = 6.46, HO = 0.66) is comparable to the highest diversity found in any koala population previously investigated. However, considerably lower genetic diversity was found in the Campbelltown population (A = 3.17, HO = 0.49), which also displayed evidence of a recent population bottleneck (effective population size estimated at 16–21). Conclusions. Animals separated by a military reserve were identified as one population, suggesting that the reserve maintains gene flow within this population. By contrast, strong differentiation of two geographically close populations separated by several potential barriers to gene flow suggested these land-use features pose barriers to gene flow. Implications. Implications of these findings for management of koala populations in the Greater Sydney region are discussed. In particular, the need to carefully consider the future of a military reserve is highlighted, along with possible solutions to allow gene flow across the proposed barrier regions. Because these are demographically separate populations, specific management plans tailored to the needs of each population will need to be formulated.

Landscape discontinuities influence gene flow and genetic structure in a large, vagile Australian mammal, Macropus fuliginosus

Large vagile mammals typically exhibit little genetic structuring across their range, particularly when their habitat is essentially continuous. We investigated the population genetic structure of a large vagile Australian macropodid, Macropus fuliginosus, which is continuously distributed across most of southern Australia, using nine highly polymorphic nuclear microsatellite loci. Five distinct genetic units were identified across the range, four on the mainland and one on Kangaroo Island. In addition to the predicted historic Nullarbor Plain Barrier, two unexpected mainland barriers to gene flow were identified. Both were associated with landscape discontinuities (Swan River, Flinders Ranges), which appear within the dispersal capabilities of M. fuliginosus. Typical of large vagile mammals, M. fuliginosus displays high genetic diversity (with the exception of an insular population) and weak genetic structuring (within genetic units). However, the expansion of M. fuliginosus from southwestern Australia during the Pleistocene has resulted in significantly reduced genetic diversity in eastern populations. No significant sex-biased dispersal was detected, although differences in habitat, densities and climatic conditions between the eastern and western regions of the range appear to influence dispersal with the effects of isolation by distance only evident in the west. These results suggest that the biogeography of southern Australia is more complex than previously thought and reveal that seemingly minor landscape features can significantly impact genetic structuring in large vagile mammals.