As anomalous heat waves are causing the widespread decline of coral reefs worldwide, there is an urgent need to identify coral populations tolerant to thermal stress. Heat stress adaptive potential is the degree of tolerance expected from evolutionary processes and, for a given reef, depends on the arrival of propagules from reefs exposed to recurrent thermal stress. For this reason, assessing spatial patterns of thermal adaptation and reef connectivity is of paramount importance to inform conservation strategies. In this work, we applied a seascape genomics framework to characterize the spatial patterns of thermal adaptation and connectivity for coral reefs of New Caledonia (Southern Pacific). In this approach, remote sensing of seascape conditions was combined with genomic data from three coral species. For every reef of the region, we computed a probability of heat stress adaptation, and two indices forecasting inbound and outbound connectivity. We then compared our indicators to field survey data, and observed that decrease of coral cover after heat stress was lower at reefs predicted with high probability of adaptation and inbound connectivity. Last, we discussed how these indicators can be used to inform local conservation strategies and preserve the adaptive potential of New Caledonian reefs.
The majority of lizards classified in the superfamily Iguanoidea have an XX/XY sex-determination system in which sex-chromosomal linkage shows homology with chicken (Gallus gallus) chromosome 15 (GGA15). However, the genomics of sex chromosomes remain largely unexplored owing to the presence of homomorphic sex chromosomes in majority of the species. Recent advances in high-throughput genome complexity reduction sequencing provide an effective approach to the identification of sex-specific loci with both single-nucleotide polymorphisms (SNPs) and restriction fragment presence/absence (PA), and a better understanding of sex chromosome dynamics in Iguanoidea. In this study, we applied Diversity Arrays Technology (DArTseqTM) in 29 phenotypic sex assignments (14 males and 15 females) of green iguana (Iguana iguana). We confirmed a male heterogametic (XX/XY) sex determination mode in this species, identifying 29 perfectly sex-linked SNP/PA loci and 164 moderately sex-linked SNP/PA loci, providing evidence probably indicative of XY recombination. Three loci from among the perfectly sex-linked SNP/PA loci showed partial homology with several amniote sex chromosomal linkages. The results support the hypothesis of an ancestral super-sex chromosome with overlaps of partial sex-chromosomal linkages. However, only one locus among the moderately sex-linked loci showed homology with GGA15, which suggests that the specific region homologous to GGA15 was located outside the non-recombination region but in close proximity to this region of the sex chromosome in green iguana. Therefore, the location of GGA15 might be further from the putative sex-determination locus in green iguana. This is a paradigm shift in understanding linkages on homomorphic X and Y sex chromosomes. The DArTseq platform provides an easy-to-use strategy for future research on the evolution of sex chromosomes in Iguanoidea, particularly for non-model species with homomorphic or highly cryptic sex chromosomes.
Conservation genomics research often relies on accurate sex information to make inferences about species demography, dispersal, and population structure. However, field determined sex data are not always available and can be subject to human error, while laboratory sex determination is costly and often challenging for non-model species. Conservation genomics programs increasingly use reduced-representation genome sequencing to assess neutral and functional genetic diversity, population structure, gene flow and pedigrees in threatened species. Here we demonstrate that sex can be determined from reduced-representation sequencing data produced by the increasingly popular Diversity Arrays Technology sequencing workflow (DArT-seq) using a program originally designed for application to shotgun data. This program – sexassign – compares the “dosage” of sequencing reads mapping to autosomes versus the X chromosome. In the present study, sexassign accurately determined the sex of 60 field-collected Greater Stick-Nest Rat (Leporillus conditor) samples, despite the absence of an annotated reference genome for the species. This “read-dosage” approach is not only more accurate and affordable than traditional sex determination methods, but can be applied to any diploid organism with a heterogametic sex determination system – including non-model and understudied species of conservation importance – by using FASTQs generated by DArT.
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.
An assessment of genetic diversity of marine populations is critical not only for the understanding and preserving natural biodiversity but also for its commercial potential. As commercial demand rises for marine resources, it is critical to generate baseline information for monitoring wild populations. Furthermore, anthropogenic stressors on the coastal environment, such as warming sea temperatures and overharvesting of wild populations, are leading to the destruction of keystone marine species such as kelps. In this study, we conducted a fine-scale genetic analysis using genome-wide high-density markers on Northwest Atlantic sugar kelp. The population structure for a total of 149 samples from the Gulf of Maine (GOM) and Southern New England (SNE) was investigated using AMOVA, FST, admixture, and PCoA. Genome-wide association analyses were conducted for six morphological traits, and the extended Lewontin and Krakauer (FLK) test was used to detect selection signatures. Our results indicate that the GOM region is more heterogeneous than SNE. These two regions have large genetic difference (between-location FST ranged from 0.21 to 0.32) and were separated by Cape Cod, which is known to be the biogeographic barrier for other taxa. We detected one significant SNP (P = 2.03 × 10–7) associated with stipe length, and 248 SNPs with higher-than-neutral differentiation. The findings of this study provide baseline knowledge on sugar kelp population genetics for future monitoring, managing and potentially restoring wild populations, as well as assisting in selective breeding to improve desirable traits for future commercialization opportunities.
Ecological restoration requires balancing levels of genetic diversity to achieve present-day establishment as well as long-term sustainability. Assumptions based on distributional, taxonomic or functional generalizations are often made when deciding how to source plant material for restoration. We investigate this assumption and ask whether species-specific data is required to optimize provenancing strategies. We use population genetic and environmental data from five congeneric and largely co-distributed species of Acacia to specifically ask how different species-specific genetic provenancing strategies are based on empirical data and how well a simple, standardized collection strategy would work when applied to the same species. We find substantial variability in terms of patterns of genetic diversity and differentiation across the landscape among these five co-distributed Acacia species. This variation translates into substantial differences in genetic provenancing recommendations among species (ranging from 100% to less than 1% of observed genetic variation across species) that could not have been accurately predicted a priori based on simple observation or overall distributional patterns. Furthermore, when a common provenancing strategy was applied to each species, the recommended collection areas and the evolutionary representativeness of such artificially standardized areas were substantially different (smaller) from those identified based on environmental and genetic data. We recommend the implementation of the increasingly accessible array of evolutionary-based methodologies and information to optimize restoration efforts
We characterize geographical variation in the genetics, ecology, morphology and plumage of the Hooded Plover Charadrius cucullatus, a threatened shorebird widespread on coasts (south-east) and saltlakes (south-west) across the southern mainland of Australia. We confirm the distinctness of populations on either side of the Nullabor Plain/Great Australian Bight indicating that this apparent biogeographic barrier to terrestrial flora and fauna also influenced the phylogeographic pattern of a resident shorebird which is dependent on access to saline/hypersaline water. Western birds are genetically distinct, larger and have darker plumage than those in the east. Distinct ecological differences also exist, with eastern birds being resident, usually occurring on ocean beaches, with generalised diets and seasonal, solitary breeding. Western birds are nomadic, breeding when suitable conditions prevail, and exhibit more specialized diets, and may at times breed semi-colonially. Within the extensive eastern population, genetic analyses revealed an isolation by distance pattern, as well as genetic evidence of a long distance dispersal event. South Australian birds are distinct from Victorian birds, but cluster closely with Victorian birds in all genetic analyses. The western population features little genetic structuring, consistent with nomadism and congregatory behavior associated with drought/flood events. The distinctness of western populations found using multiple lines of evidence and the lack of evidence for any recent migration events across the Nullabor Plain, suggests that previously described subspecies warrant reinstatement, and that these evolutionary units require conservation assessment, as has been recognized by national conservation authorities. Conservation management of the eastern population should occur across its distribution, taking care to capture the genetic variability which is evident within this population. Further exploration of genetic variation between Eastern states is suggested.
The Neotropical region exhibits the greatest worldwide diversity and the diversification history of several clades is related to the puzzling geomorphologic and climatic history of this region. The freshwater Amazon ecoregion contains the main hydrographic basins of the Neotropical region that are highly dendritic and ecologically diverse. It contains a rich and endemic fish fauna, including one of its most iconic and economically important representatives, the bony-tongue Arapaima gigas (Teleostei, Osteoglossiformes). Here, we evaluated the projected distribution of the genus in different historical periods (Present, Last Glacial Maximum, Last Interglacial Maximum and Near Future) and interpreted these results in light of the genomic diversity and modeled historical demography. For that, we combined species distribution models, population genetic analysis using SNPs and deep learning model selection. We analyzed a representative sample of the genus from the two basins where it naturally occurs, four localities in the Amazon (Am) and three in the Tocantins-Araguaia (To-Ar) basin, as well as individuals from three fish farms. We inferred a potentially smaller distribution in the glacial period, with a possible refuge in central Am. Our genetic data agrees with this result, suggesting a higher level of genetic diversity in the Am basin, compared to that observed in To-Ar. Our deep learning model comparison indicated that the To-Ar basin was colonized by the population from the Am basin. Considering a global warming scenario in the near future, A. gigas could reach an even larger range, especially if anthropogenic related dispersal occurs, potentially invading new areas and impacting their communities.
With recent advances in sequencing technology, genomic data are changing how important conservation management decisions are made. Applications such as Close-Kin Mark-Recapture demand large amounts of data to estimate population size and structure, and their full potential can only be realised through ongoing improvements in genotyping strategies. Here we introduce DArTcap, a cost-efficient method that combines DArTseq and sequence capture, and illustrate its use in a high resolution population analysis of Glyphis garricki, a rare, poorly known and threatened euryhaline shark. Clustering analyses and spatial distribution of kin pairs from four different regions across northern Australia and one in Papua New Guinea, representing its entire known range, revealed that each region hosts at least one distinct population. Further structuring is likely within Van Diemen Gulf, the region that included the most rivers sampled, suggesting additional population structuring would be found if other rivers were sampled. Coalescent analyses and spatially explicit modelling suggest that G. garricki experienced a recent range expansion during the opening of the Gulf of Carpentaria following the conclusion of the Last Glacial Maximum. The low migration rates between neighbouring populations of a species that is found only in restricted coastal and riverine habitats show the importance of managing each population separately, including careful monitoring of local and remote anthropogenic activities that may affect their environments. Overall we demonstrated how a carefully chosen SNP panel combined with DArTcap can provide highly accurate kinship inference and also support population structure and historical demography analyses, therefore maximising cost-effectiveness.
Several studies suggest the relation of DNA methylation to diseases in humans and important phenotypes in plants drawing attention to this epigenetic mark as an important source of variability. In the last decades, several methodologies were developed to assess the methylation state of a genome. However, there is still a lack of affordable and precise methods for genome wide analysis in large sample size studies. Methyl sensitive double digestion MS-DArT sequencing method emerges as a promising alternative for methylation profiling. We developed a computational pipeline for the identification of DNA methylation using MS-DArT-seq data and carried out a pilot study using the Eucalyptus grandis tree sequenced for the species reference genome. Using a statistic framework as in differential expression analysis, 72,515 genomic sites were investigated and 5,846 methylated sites identified, several tissue specific, distributed along the species 11 chromosomes. We highlight a bias towards identification of DNA methylation in genic regions and the identification of 2,783 genes and 842 transposons containing methylated sites. Comparison with WGBS, DNA sequencing after treatment with bisulfite, data demonstrated a precision rate higher than 95% for our approach. The availability of a reference genome is useful for determining the genomic context of methylated sites but not imperative, making this approach suitable for any species. Our approach provides a cost effective, broad and reliable examination of DNA methylation profile on MspI/HpaII restriction sites, is fully reproducible and the source code is available on GitHub (https://github.com/wendelljpereira/ms-dart-seq).
