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.

To understand the ecosystem dynamics that underpin the year-round presence of a large generalist consumer, the Bryde’s whale (Balaenoptera edeni brydei), we use a DNA metabarcoding approach and systematic zooplankton surveys to investigate seasonal and regional changes in zooplankton communities and if whale diet reflects such changes. Twenty-four zooplankton community samples were collected from three regions throughout the Hauraki Gulf, New Zealand, over two temperature regimes (warm and cool seasons), as well as 20 samples of opportunistically collected Bryde’s whale scat. Multi-locus DNA barcode libraries were constructed from 18S and COI gene fragments, representing a trade-off between identification and resolution of metazoan taxa. Zooplankton community OTU occurrence and relative read abundance showed regional and seasonal differences based on permutational analyses of variance in both DNA barcodes, with significant changes in biodiversity indices linked to season in COI only. In contrast, we did not find evidence that Bryde’s whale diet shows seasonal or regional trends, but instead indicated clear prey preferences for krill-like crustaceans, copepods, salps and ray-finned fishes independent of prey availability. The year-round presence of Bryde’s whales in the Hauraki Gulf is likely associated with the patterns of distribution and abundance of these key prey items.

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.

BackgroundOur 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 FindingsAffymetrix 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/SignificanceThis 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.

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.