Population genomics of a predatory mammal reveals patterns of decline and impacts of exposure to toxic toads.

Mammal declines across northern Australia are one of the major biodiversity loss events occurring globally. There has been no regional assessment of the implications of these species declines for genomic diversity. To address this, we conducted a species-wide assessment of genomic diversity in the northern quoll (Dasyurus hallucatus), an Endangered marsupial carnivore. We used next generation sequencing methods to genotype 10,191 single nucleotide polymorphisms (SNPs) in 352 individuals from across a 3220-km length of the continent, investigating patterns of population genomic structure and diversity, and identifying loci showing signals of putative selection. We found strong heterogeneity in the distribution of genomic diversity across the continent, characterized by (i) biogeographical barriers driving hierarchical population structure through long-term isolation, and (ii) severe reductions in diversity resulting from population declines, exacerbated by the spread of introduced toxic cane toads (Rhinella marina). These results warn of a large ongoing loss of genomic diversity and associated adaptive capacity as mammals decline across northern Australia. Encouragingly, populations of the northern quoll established on toad-free islands by translocations appear to have maintained most of the initial genomic diversity after 16 years. By mapping patterns of genomic diversity within and among populations, and investigating these patterns in the context of population declines, we can provide conservation managers with data critical to informed decision-making. This includes the identification of populations that are candidates for genetic management, the importance of remnant island and insurance/translocated populations for the conservation of genetic diversity, and the characterization of putative evolutionarily significant units.

An assembly-free method of phylogeny reconstruction using short-read sequences from pooled samples without barcodes.

A current strategy for obtaining haplotype information from several individuals involves short-read sequencing of pooled amplicons, where fragments from each individual is identified by a unique DNA barcode. In this paper, we report a new method to recover the phylogeny of haplotypes from short-read sequences obtained using pooled amplicons from a mixture of individuals, without barcoding. The method, AFPhyloMix, accepts an alignment of the mixture of reads against a reference sequence, obtains the single-nucleotide-polymorphisms (SNP) patterns along the alignment, and constructs the phylogenetic tree according to the SNP patterns. AFPhyloMix adopts a Bayesian inference model to estimate the phylogeny of the haplotypes and their relative abundances, given that the number of haplotypes is known. In our simulations, AFPhyloMix achieved at least 80% accuracy at recovering the phylogenies and relative abundances of the constituent haplotypes, for mixtures with up to 15 haplotypes. AFPhyloMix also worked well on a real data set of kangaroo mitochondrial DNA sequences.

Loss of vocal culture and fitness costs in a critically endangered songbird.

Cultures in humans and other species are maintained through interactions among conspecifics. Declines in population density could be exacerbated by culture loss, thereby linking culture to conservation. We combined historical recordings, citizen science and breeding data to assess the impact of severe population decline on song culture, song complexity and individual fitness in critically endangered regent honeyeaters (Anthochaera phrygia). Song production in the remaining wild males varied dramatically, with 27% singing songs that differed from the regional cultural norm. Twelve per cent of males, occurring in areas of particularly low population density, completely failed to sing any species-specific songs and instead sang other species' songs. Atypical song production was associated with reduced individual fitness, as males singing atypical songs were less likely to pair or nest than males that sang the regional cultural norm. Songs of captive-bred birds differed from those of all wild birds. The complexity of regent honeyeater songs has also declined over recent decades. We therefore provide rare evidence that a severe decline in population density is associated with the loss of vocal culture in a wild animal, with concomitant fitness costs for remaining individuals. The loss of culture may be a precursor to extinction in declining populations that learn selected behaviours from conspecifics, and therefore provides a useful conservation indicator.

Correction to: Effective machine-learning assembly for next-generation amplicon sequencing with very low coverage.

Following publication of the original article [1], the author reported that there are several errors in the original article.

Effective machine-learning assembly for next-generation amplicon sequencing with very low coverage.

BACKGROUND: In short-read DNA sequencing experiments, the read coverage is a key parameter to successfully assemble the reads and reconstruct the sequence of the input DNA. When coverage is very low, the original sequence reconstruction from the reads can be difficult because of the occurrence of uncovered gaps. Reference guided assembly can then improve these assemblies. However, when the available reference is phylogenetically distant from the sequencing reads, the mapping rate of the reads can be extremely low. Some recent improvements in read mapping approaches aim at modifying the reference according to the reads dynamically. Such approaches can significantly improve the alignment rate of the reads onto distant references but the processing of insertions and deletions remains challenging. RESULTS: Here, we introduce a new algorithm to update the reference sequence according to previously aligned reads. Substitutions, insertions and deletions are performed in the reference sequence dynamically. We evaluate this approach to assemble a western-grey kangaroo mitochondrial amplicon. Our results show that more reads can be aligned and that this method produces assemblies of length comparable to the truth while limiting error rate when classic approaches fail to recover the correct length. Finally, we discuss how the core algorithm of this method could be improved and combined with other approaches to analyse larger genomic sequences. CONCLUSIONS: We introduced an algorithm to perform dynamic alignment of reads on a distant reference. We showed that such approach can improve the reconstruction of an amplicon compared to classically used bioinformatic pipelines. Although not portable to genomic scale in the current form, we suggested several improvements to be investigated to make this method more flexible and allow dynamic alignment to be used for large genome assemblies.

HaploJuice : accurate haplotype assembly from a pool of sequences with known relative concentrations.

BACKGROUND: Pooling techniques, where multiple sub-samples are mixed in a single sample, are widely used to take full advantage of high-throughput DNA sequencing. Recently, Ranjard et al. (PLoS ONE 13:0195090, 2018) proposed a pooling strategy without the use of barcodes. Three sub-samples were mixed in different known proportions (i.e. 62.5%, 25% and 12.5%), and a method was developed to use these proportions to reconstruct the three haplotypes effectively. RESULTS: HaploJuice provides an alternative haplotype reconstruction algorithm for Ranjard et al.'s pooling strategy. HaploJuice significantly increases the accuracy by first identifying the empirical proportions of the three mixed sub-samples and then assembling the haplotypes using a dynamic programming approach. HaploJuice was evaluated against five different assembly algorithms, Hmmfreq (Ranjard et al., PLoS ONE 13:0195090, 2018), ShoRAH (Zagordi et al., BMC Bioinformatics 12:119, 2011), SAVAGE (Baaijens et al., Genome Res 27:835-848, 2017), PredictHaplo (Prabhakaran et al., IEEE/ACM Trans Comput Biol Bioinform 11:182-91, 2014) and QuRe (Prosperi and Salemi, Bioinformatics 28:132-3, 2012). Using simulated and real data sets, HaploJuice reconstructed the true sequences with the highest coverage and the lowest error rate. CONCLUSION: HaploJuice provides high accuracy in haplotype reconstruction, making Ranjard et al.'s pooling strategy more efficient, feasible, and applicable, with the benefit of reducing the sequencing cost.

Modelling competition and dispersal in a statistical phylogeographic framework.

Competition between organisms influences the processes governing the colonization of new habitats. As a consequence, species or populations arriving first at a suitable location may prevent secondary colonization. Although adaptation to environmental variables (e.g., temperature, altitude, etc.) is essential, the presence or absence of certain species at a particular location often depends on whether or not competing species co-occur. For example, competition is thought to play an important role in structuring mammalian communities assembly. It can also explain spatial patterns of low genetic diversity following rapid colonization events or the "progression rule" displayed by phylogenies of species found on archipelagos. Despite the potential of competition to maintain populations in isolation, past quantitative analyses have largely ignored it because of the difficulty in designing adequate methods for assessing its impact. We present here a new model that integrates competition and dispersal into a Bayesian phylogeographic framework. Extensive simulations and analysis of real data show that our approach clearly outperforms the traditional Mantel test for detecting correlation between genetic and geographic distances. But most importantly, we demonstrate that competition can be detected with high sensitivity and specificity from the phylogenetic analysis of genetic variation in space.

Integration over song classification replicates: song variant analysis in the hihi.

Human expert analyses are commonly used in bioacoustic studies and can potentially limit the reproducibility of these results. In this paper, a machine learning method is presented to statistically classify avian vocalizations. Automated approaches were applied to isolate bird songs from long field recordings, assess song similarities, and classify songs into distinct variants. Because no positive controls were available to assess the true classification of variants, multiple replicates of automatic classification of song variants were analyzed to investigate clustering uncertainty. The automatic classifications were more similar to the expert classifications than expected by chance. Application of these methods demonstrated the presence of discrete song variants in an island population of the New Zealand hihi (Notiomystis cincta). The geographic patterns of song variation were then revealed by integrating over classification replicates. Because this automated approach considers variation in song variant classification, it reduces potential human bias and facilitates the reproducibility of the results.

Soil bacterial community composition is more stable in kiwifruit orchards relative to phyllosphere communities over time.

BACKGROUND: Soil and phyllosphere (leaves and fruit) microbes play critical roles in the productivity and health of crops. However, microbial community dynamics are currently understudied in orchards, with a limited number incorporating temporal monitoring. We used 16S rRNA gene amplicon sequencing to investigate bacterial community temporal dynamics and community assembly processes on the leaves and fruit, and in the soil of 12 kiwifruit orchards across a cropping season in New Zealand. RESULTS: Community composition significantly differed (P < 0.001) among the three sample types. However, the communities in the phyllosphere substrates more closely resembled each other, relative to the communities in the soil. There was more temporal stability in the soil bacterial community composition, relative to the communities residing on the leaves and fruit, and low similarity between the belowground and aboveground communities. Bacteria in the soil were more influenced by deterministic processes, while stochastic processes were more important for community assembly in the phyllosphere. CONCLUSIONS: The higher temporal variability and the stochastic nature of the community assembly processes observed in the phyllosphere communities highlights why predicting the responsiveness of phyllosphere communities to environmental change, or the likelihood of pathogen invasion, can be challenging. The relative temporal stability and the influence of deterministic selection on soil microbial communities suggests a greater potential for their prediction and reliable manipulation.

pgHMA: Application of the heteroduplex mobility assay analysis in phylogenetics and population genetics.

The heteroduplex mobility assay (HMA) has proven to be a robust tool for the detection of genetic variation. Here, we describe a simple and rapid application of the HMA by microfluidic capillary electrophoresis, for phylogenetics and population genetic analyses (pgHMA). We show how commonly applied techniques in phylogenetics and population genetics have equivalents with pgHMA: phylogenetic reconstruction with bootstrapping, skyline plots, and mismatch distribution analysis. We assess the performance and accuracy of pgHMA by comparing the results obtained against those obtained using standard methods of analyses applied to sequencing data. The resulting comparisons demonstrate that: (a) there is a significant linear relationship (R2  = .992) between heteroduplex mobility and genetic distance, (b) phylogenetic trees obtained by HMA and nucleotide sequences present nearly identical topologies, (c) clades with high pgHMA parametric bootstrap support also have high bootstrap support on nucleotide phylogenies, (d) skyline plots estimated from the UPGMA trees of HMA and Bayesian trees of nucleotide data reveal similar trends, especially for the median trend estimate of effective population size, and (e) optimized mismatch distributions of HMA are closely fitted to the mismatch distributions of nucleotide sequences. In summary, pgHMA is an easily-applied method for approximating phylogenetic diversity and population trends.

Assemblage of Focal Species Recognizers-AFSR: A technique for decreasing false indications of presence from acoustic automatic identification in a multiple species context.

Passive acoustic monitoring (PAM) coupled with automated species identification is a promising tool for species monitoring and conservation worldwide. However, high false indications of presence are still an important limitation and a crucial factor for acceptance of these techniques in wildlife surveys. Here we present the Assemblage of Focal Species Recognizers-AFSR, a novel approach for decreasing false positives and increasing models' precision in multispecies contexts. AFSR focusses on decreasing false positives by excluding unreliable sound file segments that are prone to misidentification. We used MatlabHTK, a hidden Markov models interface for bioacoustics analyses, for illustrating AFSR technique by comparing two approaches, 1) a multispecies recognizer where all species are identified simultaneously, and 2) an assemblage of focal species recognizers (AFSR), where several recognizers that each prioritise a single focal species are then summarised into a single output, according to a set of rules designed to exclude unreliable segments. Both approaches (the multispecies recognizer and AFSR) used the same sound files training dataset, but different processing workflow. We applied these recognisers to PAM recordings from a remote island colony with five seabird species and compared their outputs with manual species identifications. False positives and precision improved for all the five species when using AFSR, achieving remarkable 0% false positives and 100% precision for three of five seabird species, and < 6% false positives, and >90% precision for the other two species. AFSR' output was also used to generate daily calling activity patterns for each species. Instead of attempting to withdraw useful information from every fragment in a sound recording, AFSR prioritises more trustworthy information from sections with better quality data. AFSR can be applied to automated species identification from multispecies PAM recordings worldwide.

Genomic impact of severe population decline in a nomadic songbird.

Uncovering the population genetic histories of non-model organisms is increasingly possible through advances in next generation sequencing and DNA sampling of museum specimens. This new information can inform conservation of threatened species, particularly those for which historical and contemporary population data are unavailable or challenging to obtain. The critically endangered, nomadic regent honeyeater Anthochaera phrygia was abundant and widespread throughout south-eastern Australia prior to a rapid population decline and range contraction since the 1970s. A current estimated population of 250-400 individuals is distributed sparsely across 600,000 km2 from northern Victoria to southern Queensland. Using hybridization RAD (hyRAD) techniques, we obtained a SNP dataset from 64 museum specimens (date 1879-1960), 102 'recent' (1989-2012) and 52 'current' (2015-2016) wild birds sampled throughout the historical and contemporary range. We aimed to estimate population genetic structure, genetic diversity and population size of the regent honeyeater prior to its rapid decline. We then assessed the impact of the decline on recent and current population size, structure and genetic diversity. Museum sampling showed population structure in regent honeyeaters was historically low, which remains the case despite a severe fragmentation of the breeding range. Population decline has led to minimal loss of genetic diversity since the 1980's. Capacity to quantify the overall magnitude of both genetic diversity loss and population decline was limited by the poorer quality of genomic data derived from museum specimens. A rapid population decline, coupled with the regent honeyeater's high mobility, means a detectable genomic impact of this decline has not yet manifested. Extinction may occur in this nomadic species before a detectable genomic impact of small population size is realised. We discuss the implications for genetic management of endangered mobile species and enhancing the value of museum specimens in population genomic studies.

Unsupervised bird song syllable classification using evolving neural networks.

Evolution of bird vocalizations is subjected to selection pressure related to their functions. Passerine bird songs are also under a neutral model of evolution because of the learning process supporting their transmission; thus they contain signals of individual, population, and species relationships. In order to retrieve this information, large amounts of data need to be processed. From vocalization recordings, songs are extracted and encoded as sequences of syllables before being compared. Encoding songs in such a way can be done either by ear and spectrogram visual analysis or by specific algorithms permitting reproducible studies. Here, a specific automatic method is presented to compute a syllable distance measure allowing an unsupervised classification of song syllables. Results obtained from the encoding of White-crowned Sparrow (Zonotrichia leucophrys pugetensis) songs are compared to human-based analysis.

Reassembling haplotypes in a mixture of pooled amplicons when the relative concentrations are known: A proof-of-concept study on the efficient design of next-generation sequencing strategies.

Next-generation sequencing can be costly and labour intensive. Usually, the sequencing cost per sample is reduced by pooling amplified DNA = amplicons) derived from different individuals on the same sequencing lane. Barcodes unique to each amplicon permit short-read sequences to be assigned appropriately. However, the cost of the library preparation increases with the number of barcodes used. We propose an alternative to barcoding: by using different known proportions of individually-derived amplicons in a pooled sample, each is characterised a priori by an expected depth of coverage. We have developed a Hidden Markov Model that uses these expected proportions to reconstruct the input sequences. We apply this method to pools of mitochondrial DNA amplicons extracted from kangaroo meat, genus Macropus. Our experiments indicate that the sequence coverage can be efficiently used to index the short-reads and that we can reassemble the input haplotypes when secondary factors impacting the coverage are controlled. We therefore demonstrate that, by combining our approach with standard barcoding, the cost of the library preparation is reduced to a third.

Complete mitochondrial genome of the green-lipped mussel, Perna canaliculus (Mollusca: Mytiloidea), from long nanopore sequencing reads.

We describe here the first complete genome assembly of the New Zealand green-lipped mussel, Perna canaliculus, mitochondrion. The assembly was performed de novo from a mix of long nanopore sequencing reads and short sequencing reads. The genome is 16,005 bp long. Comparison to other Mytiloidea mitochondrial genomes indicates important gene rearrangements in this family.

Visualization and Quantification of Mesenchymal Cell Adipogenic Differentiation Potential with a Lineage Specific Marker.

Several dyes are currently available for use in detecting differentiation of mesenchymal cells into adipocytes. Dyes, such as Oil Red O, are cheap, easy to use and widely utilized by laboratories analyzing the adipogenic potential of mesenchymal cells. However, they are not specific to changes in gene transcription. We have developed a gene-specific differentiation assay to analyze when a mesenchymal cell has switched its fate to an adipogenic lineage. Immuno-labelling against fatty acid binding protein-4 (FABP4), a lineage-specific marker of adipogenic differentiation, enabled visualization and quantification of differentiated cells. The ability to quantify adipogenic differentiation potential of mesenchymal cells in a 96 well microplate format has promising implications for a number of applications. Hundreds of clinical trials involve the use of adult mesenchymal stromal cells and it is currently difficult to correlate therapeutic outcomes within and especially between such clinical trials. This simple high-throughput FABP4 assay provides a quantitative assay for assessing the differentiation potential of patient-derived cells and is a robust tool for comparing different isolation and expansion methods. This is particularly important given the increasing recognition of the heterogeneity of the cells being administered to patients in mesenchymal cell products. The assay also has potential utility in high throughput drug screening, particularly in obesity and pre-diabetes research.

Cell cycle gene networks are associated with melanoma prognosis.

BACKGROUND: Our understanding of the molecular pathways that underlie melanoma remains incomplete. Although several published microarray studies of clinical melanomas have provided valuable information, we found only limited concordance between these studies. Therefore, we took an in vitro functional genomics approach to understand melanoma molecular pathways. METHODOLOGY/PRINCIPAL FINDINGS: Affymetrix microarray data were generated from A375 melanoma cells treated in vitro with siRNAs against 45 transcription factors and signaling molecules. Analysis of this data using unsupervised hierarchical clustering and Bayesian gene networks identified proliferation-association RNA clusters, which were co-ordinately expressed across the A375 cells and also across melanomas from patients. The abundance in metastatic melanomas of these cellular proliferation clusters and their putative upstream regulators was significantly associated with patient prognosis. An 8-gene classifier derived from gene network hub genes correctly classified the prognosis of 23/26 metastatic melanoma patients in a cross-validation study. Unlike the RNA clusters associated with cellular proliferation described above, co-ordinately expressed RNA clusters associated with immune response were clearly identified across melanoma tumours from patients but not across the siRNA-treated A375 cells, in which immune responses are not active. Three uncharacterised genes, which the gene networks predicted to be upstream of apoptosis- or cellular proliferation-associated RNAs, were found to significantly alter apoptosis and cell number when over-expressed in vitro. CONCLUSIONS/SIGNIFICANCE: This analysis identified co-expression of RNAs that encode functionally-related proteins, in particular, proliferation-associated RNA clusters that are linked to melanoma patient prognosis. Our analysis suggests that A375 cells in vitro may be valid models in which to study the gene expression modules that underlie some melanoma biological processes (e.g., proliferation) but not others (e.g., immune response). The gene expression modules identified here, and the RNAs predicted by Bayesian network inference to be upstream of these modules, are potential prognostic biomarkers and drug targets.

Contemporary and historical separation of transequatorial migration between genetically distinct seabird populations.

Pelagic seabirds are highly mobile, reducing the likelihood of allopatric speciation where disruption of gene flow between populations is caused by physically insurmountable, extrinsic barriers. Spatial segregation during the non-breeding season appears to provide an intrinsic barrier to gene flow among seabird populations that otherwise occupy nearby or overlapping regions during breeding, but how this is achieved remains unclear. Here we show that the two genetically distinct populations of Cook's petrel (Pterodroma cookii) exhibit transequatorial separation of non-breeding ranges at contemporary (ca. 2-3 yrs) and historical (ca. 100 yrs) time scales. Segregation during the non-breeding season per se appears as an unlikely barrier to gene flow. Instead we provide evidence that habitat specialization during the non-breeding season is associated with breeding asynchrony which, in conjunction with philopatry, restricts gene flow. Habitat specialization during breeding and non-breeding likely promotes evolutionary divergence between these two populations via local adaptation.