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 3 m to >100 m. 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 âˆ¼30 mya 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" (>100 m) 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.

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.

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 FST = 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.

Development and validation of a RAD-Seq target-capture based genotyping assay for routine application in advanced black tiger shrimp (Penaeus monodon) breeding programs.

BACKGROUND: The development of genome-wide genotyping resources has provided terrestrial livestock and crop industries with the unique ability to accurately assess genomic relationships between individuals, uncover the genetic architecture of commercial traits, as well as identify superior individuals for selection based on their specific genetic profile. Utilising recent advancements in de-novo genome-wide genotyping technologies, it is now possible to provide aquaculture industries with these same important genotyping resources, even in the absence of existing genome assemblies. Here, we present the development of a genome-wide SNP assay for the Black Tiger shrimp (Penaeus monodon) through utilisation of a reduced-representation whole-genome genotyping approach (DArTseq). RESULTS: Based on a single reduced-representation library, 31,262 polymorphic SNPs were identified across 650 individuals obtained from Australian wild stocks and commercial aquaculture populations. After filtering to remove SNPs with low read depth, low MAF, low call rate, deviation from HWE, and non-Mendelian inheritance, 7542 high-quality SNPs were retained. From these, 4236 high-quality genome-wide loci were selected for baits-probe development and 4194 SNPs were included within a finalized target-capture genotype-by-sequence assay (DArTcap). This assay was designed for routine and cost effective commercial application in large scale breeding programs, and demonstrates higher confidence in genotype calls through increased call rate (from 80.2 ± 14.7 to 93.0% ± 3.5%), increased read depth (from 20.4 ± 15.6 to 80.0 ± 88.7), as well as a 3-fold reduction in cost over traditional genotype-by-sequencing approaches. CONCLUSION: Importantly, this assay equips the P. monodon industry with the ability to simultaneously assign parentage of communally reared animals, undertake genomic relationship analysis, manage mate pairings between cryptic family lines, as well as undertake advance studies of genome and trait architecture. Critically this assay can be cost effectively applied as P. monodon breeding programs transition to undertaking genomic selection.

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.

Infection dynamics, dispersal, and adaptation: understanding the lack of recovery in a remnant frog population following a disease outbreak.

Emerging infectious diseases can cause dramatic declines in wildlife populations. Sometimes, these declines are followed by recovery, but many populations do not recover. Studying differential recovery patterns may yield important information for managing disease-afflicted populations and facilitating population recoveries. In the late 1980s, a chytridiomycosis outbreak caused multiple frog species in Australia's Wet Tropics to decline. Populations of some species (e.g., Litoria nannotis) subsequently recovered, while others (e.g., Litoria dayi) did not. We examined the population genetics and current infection status of L. dayi, to test several hypotheses regarding the failure of its populations to recover: (1) a lack of individual dispersal abilities has prevented recolonization of previously occupied locations, (2) a loss of genetic variation has resulted in limited adaptive potential, and (3) L. dayi is currently adapting to chytridiomycosis. We found moderate-to-high levels of gene flow and diversity (Fst range: <0.01-0.15; minor allele frequency (MAF): 0.192-0.245), which were similar to previously published levels for recovered L. nannotis populations. This suggests that dispersal ability and genetic diversity do not limit the ability of L. dayi to recolonize upland sites. Further, infection intensity and prevalence increased with elevation, suggesting that chytridiomycosis is still limiting the elevational range of L. dayi. Outlier tests comparing infected and uninfected individuals consistently identified 18 markers as putatively under selection, and several of those markers matched genes that were previously implicated in infection. This suggests that L. dayi has genetic variation for genes that affect infection dynamics and may be undergoing adaptation.

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 Ne 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%; Ne  = 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 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; AvHecorr=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.

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.

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.

Reversing the decline of threatened koala (Phascolarctos cinereus) populations in New South Wales: Using genomics to enhance conservation outcomes.

Genetic management is a critical component of threatened species conservation. Understanding spatial patterns of genetic diversity is essential for evaluating the resilience of fragmented populations to accelerating anthropogenic threats. Nowhere is this more relevant than on the Australian continent, which is experiencing an ongoing loss of biodiversity that exceeds any other developed nation. Using a proprietary genome complexity reduction-based method (DArTSeq), we generated a data set of 3239 high quality Single Nucleotide Polymorphisms (SNPs) to investigate spatial patterns and indices of genetic diversity in the koala (Phascolarctos cinereus), a highly specialised folivorous marsupial that is experiencing rapid and widespread population declines across much of its former range. Our findings demonstrate that current management divisions across the state of New South Wales (NSW) do not fully represent the distribution of genetic diversity among extant koala populations, and that care must be taken to ensure that translocation paradigms based on these frameworks do not inadvertently restrict gene flow between populations and regions that were historically interconnected. We also recommend that koala populations should be prioritised for conservation action based on the scale and severity of the threatening processes that they are currently faced with, rather than placing too much emphasis on their perceived value (e.g., as reservoirs of potentially adaptive alleles), as our data indicate that existing genetic variation in koalas is primarily partitioned among individual animals. As such, the extirpation of koalas from any part of their range represents a potentially critical reduction of genetic diversity for this iconic Australian species.

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.

The El Niño Southern Oscillation drives multidirectional inter-reef larval connectivity in the Great Barrier Reef.

The El Niño Southern Oscillation (ENSO) is the strongest source of interannual global climate variability, and extreme ENSO events are projected to increase in frequency under climate change. Interannual variability in the Coral Sea circulation has been associated with ENSO, although uncertainty remains regarding ENSO's influence on hydrodynamics and larval dispersal in the adjacent Great Barrier Reef (GBR). We investigated larval connectivity during ENSO events from 2010 to 2017 throughout the GBR, based on biophysical modelling of a widespread predatory reef fish, Lutjanus carponotatus. Our results indicate a well-connected system over the study period with high interannual variability in inter-reef connectivity associated with ENSO. Larval connectivity patterns were highly correlated to variations in the Southern Oscillation Index (SOI). During El Niño conditions and periods of weak SOI, larval dispersal patterns were predominantly poleward in the central and southern regions, reversing to a predominant equatorward flow during very strong SOI and extreme La Niña conditions. These ENSO-linked connectivity patterns were associated with positive connectivity anomalies among reefs. Our findings identify ENSO as an important source of variation in larval dispersal and connectivity patterns in the GBR, which can influence the stability of population dynamics and patterns of biodiversity in the region.

Genetic approaches for increasing fitness in endangered species.

The global rate of wildlife extinctions is accelerating, and the persistence of many species requires conservation breeding programs. A central paradigm of these programs is to preserve the genetic diversity of the founder populations. However, this may preserve original characteristics that make them vulnerable to extinction. We introduce targeted genetic intervention (TGI) as an alternative approach that promotes traits that enable species to persist in the face of threats by changing the incidence of alleles that impact on fitness. The TGI toolkit includes methods with established efficacy in model organisms and agriculture but are largely untried for conservation, such as synthetic biology and artificial selection. We explore TGI approaches as a species-restoration tool for intractable threats including infectious disease and climate change.

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.

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/.