ALBERT

Description: The predicted effect of effective population size on the distribution of fitness effects and substitution rate is critically dependent on the relationship between sequence and fitness. This highlights the importance of using models that are informed by the molecular biology, biochemistry, and biophysics of the evolving systems. We describe a computational model based on fundamental aspects of biophysics, the requirement for (most) proteins to be thermodynamically stable. Using this model, we find that differences in population size have minimal impact on the distribution of population-scaled fitness effects, as well as on the rate of molecular evolution. This is because larger populations result in selection for more stable proteins that are less affected by mutations. This reduction in the magnitude of the fitness effects almost exactly cancels the greater selective pressure resulting from the larger population size. Conversely, changes in the population size in either direction cause transient increases in the substitution rate. As differences in population size often correspond to changes in population size, this makes comparisons of substitution rates in different lineages difficult to interpret.

Description: Populations of widely distributed species encounter and must adapt to local environmental conditions. However, comprehensive characterization of the genetic basis of adaptation is demanding, requiring genome-wide genotype data, multiple sampled populations, and an understanding of population structure and potential selection pressures. Here, we used single-nucleotide polymorphism genotyping and data on numerous environmental variables to describe the genetic basis of local adaptation in 21 populations of teosinte, the wild ancestor of maize. We found complex hierarchical genetic structure created by altitude, dispersal events, and admixture among subspecies, which complicated identification of locally beneficial alleles. Patterns of linkage disequilibrium revealed four large putative inversion polymorphisms showing clinal patterns of frequency. Population differentiation and environmental correlations suggest that both inversions and intergenic polymorphisms are involved in local adaptation.

Description: Wolbachia , endosymbiotic bacteria of the order Rickettsiales, are widespread in arthropods but also present in nematodes. In arthropods, A and B supergroup Wolbachia are generally associated with distortion of host reproduction. In filarial nematodes, including some human parasites, multiple lines of experimental evidence indicate that C and D supergroup Wolbachia are essential for the survival of the host, and here the symbiotic relationship is considered mutualistic. The origin of this mutualistic endosymbiosis is of interest for both basic and applied reasons: How does a parasite become a mutualist? Could intervention in the mutualism aid in treatment of human disease? Correct rooting and high-quality resolution of Wolbachia relationships are required to resolve this question. However, because of the large genetic distance between Wolbachia and the nearest outgroups, and the limited number of genomes so far available for large-scale analyses, current phylogenies do not provide robust answers. We therefore sequenced the genome of the D supergroup Wolbachia endosymbiont of Litomosoides sigmodontis , revisited the selection of loci for phylogenomic analyses, and performed a phylogenomic analysis including available complete genomes (from isolates in supergroups A, B, C, and D). Using 90 orthologous genes with reliable phylogenetic signals, we obtained a robust phylogenetic reconstruction, including a highly supported root to the Wolbachia phylogeny between a (A + B) clade and a (C + D) clade. Although we currently lack data from several Wolbachia supergroups, notably F, our analysis supports a model wherein the putatively mutualist endosymbiotic relationship between Wolbachia and nematodes originated from a single transition event.

Description: Many insects rely on bacterial symbionts with tiny genomes specialized for provisioning nutrients lacking in host diets. Xylem sap and phloem sap are both deficient as insect diets, but differ dramatically in nutrient content, potentially affecting symbiont genome evolution. For sap-feeding insects, sequenced symbiont genomes are available only for phloem-feeding examples from the suborder Sternorrhyncha and xylem-feeding examples from the suborder Auchenorrhyncha, confounding comparisons. We sequenced genomes of the obligate symbionts, Sulcia muelleri and Nasuia deltocephalinicola , of the phloem-feeding pest insect, Macrosteles quadrilineatus (Auchenorrhyncha: Cicadellidae). Our results reveal that Nasuia- ALF has the smallest bacterial genome yet sequenced (112 kb), and that the Sulcia- ALF genome (190 kb) is smaller than that of Sulcia in other insect lineages. Together, these symbionts retain the capability to synthesize the 10 essential amino acids, as observed for several symbiont pairs from xylem-feeding Auchenorrhyncha. Nasuia retains genes enabling synthesis of two amino acids, DNA replication, transcription, and translation. Both symbionts have lost genes underlying ATP synthesis through oxidative phosphorylation, possibly as a consequence of the enriched sugar content of phloem. Shared genomic features, including reassignment of the UGA codon from Stop to tryptophan, and phylogenetic results suggest that Nasuia -ALF is most closely related to Zinderia , the betaproteobacterial symbiont of spittlebugs. Thus, Nasuia / Zinderia and Sulcia likely represent ancient associates that have co-resided in hosts since the divergence of leafhoppers and spittlebugs 〉200 Ma, and possibly since the origin of the Auchenorrhyncha, 〉260 Ma.

Description: Thanks to the microarray technology, our understanding of transcriptome evolution at the genome level has been considerably advanced in the past decade. Yet, further investigation was challenged by several technical limitations of this technology. Recent innovation of next-generation sequencing, particularly the invention of RNA-seq technology, has shed insightful lights on resolving this problem. Though a number of statistical and computational methods have been developed to analyze RNA-seq data, the analytical framework specifically designed for evolutionary genomics remains an open question. In this article we develop a new method for estimating the genome expression distance from the RNA-seq data, which has explicit interpretations under the model of gene expression evolution. Moreover, this distance measure takes the data overdispersion, gene length variation, and sequencing depth variation into account so that it can be applied to multiple genomes from different species. Using mammalian RNA-seq data as example, we demonstrated that this expression distance is useful in phylogenomic analysis.

Description: The control of RNA splicing is often modulated by exonic motifs near splice sites. Chief among these are exonic splice enhancers (ESEs). Well-described ESEs in mammals are purine rich and cause predictable skews in codon and amino acid usage toward exonic ends. Looking across species, those with relatively abundant intronic sequence are those with the more profound end of exon skews, indicative of exonization of splice site recognition. To date, the only intron-rich species that have been analyzed are mammals, precluding any conclusions about the likely ancestral condition. Here, we examine the patterns of codon and amino acid usage in the vicinity of exon–intron junctions in the brown alga Ectocarpus siliculosus , a species with abundant large introns, known SR proteins, and classical splice sites. We find that amino acids and codons preferred/avoided at both 3' and 5' ends in Ectocarpus , of which there are many, tend, on average, to also be preferred/avoided at the same exon ends in humans. Moreover, the preferences observed at the 5' ends of exons are largely the same as those at the 3' ends, a symmetry trend only previously observed in animals. We predict putative hexameric ESEs in Ectocarpus and show that these are purine rich and that there are many more of these identified as functional ESEs in humans than expected by chance. These results are consistent with deep phylogenetic conservation of SR protein binding motifs. Assuming codons preferred near boundaries are "splice optimal" codons, in Ectocarpus , unlike Drosophila, splice optimal and translationally optimal codons are not mutually exclusive. The exclusivity of translationally optimal and splice optimal codon sets is thus not universal.

Description: Autonomous retrotransposons lacking long terminal repeats (LTR) account for much of the variation in genome size and structure among vertebrates. Mammalian genomes contain hundreds of thousands of non-LTR retrotransposon copies, mostly resulting from the amplification of a single clade known as L1. The genomes of teleost fish and squamate reptiles contain a much more diverse array of non-LTR retrotransposon families, whereas copy number is relatively low. The majority of non-LTR retrotransposon insertions in nonmammalian vertebrates also appear to be very recent, suggesting strong purifying selection limits the accumulation of non-LTR retrotransposon copies. It is however unclear whether this turnover model, originally proposed in Drosophila , applies to nonmammalian vertebrates. Here, we studied the population dynamics of L1 in the green anole lizard ( Anolis carolinensis ). We found that although most L1 elements are recent in this genome, truncated insertions accumulate readily, and many are fixed at both the population and species level. In contrast, full-length L1 insertions are found at lower population frequencies, suggesting that the turnover model only applies to longer L1 elements in Anolis . We also found that full-length L1 inserts are more likely to be fixed in populations of small effective size, suggesting that the strength of purifying selection against deleterious alleles is highly dependent on host demographic history. Similar mechanisms seem to be controlling the fate of non-LTR retrotransposons in both Anolis and teleostean fish, which suggests that mammals have considerably diverged from the ancestral vertebrate in terms of how they interact with their intragenomic parasites.

Description: Increasing evidence from sequence data from various environments, including the human gut, suggests the existence of a previously unknown putative seventh order of methanogens. The first genomic data from members of this lineage, Methanomassiliicoccus luminyensis and " Candidatus Methanomethylophilus alvus," provide insights into its evolutionary history and metabolic features. Phylogenetic analysis of ribosomal proteins robustly indicates a monophyletic group independent of any previously known methanogenic order, which shares ancestry with the Marine Benthic Group D, the Marine Group II, the DHVE2 group, and the Thermoplasmatales. This phylogenetic position, along with the analysis of enzymes involved in core methanogenesis, strengthens a single ancient origin of methanogenesis in the Euryarchaeota and indicates further multiple independent losses of this metabolism in nonmethanogenic lineages than previously suggested. Genomic analysis revealed an unprecedented loss of the genes coding for the first six steps of methanogenesis from H 2 /CO 2 and the oxidative part of methylotrophic methanogenesis, consistent with the fact that M. luminyensis and " Ca. M. alvus" are obligate H 2 -dependent methylotrophic methanogens. Genomic data also suggest that these methanogens may use a large panel of methylated compounds. Phylogenetic analysis including homologs retrieved from environmental samples indicates that methylotrophic methanogenesis (regardless of dependency on H 2 ) is not restricted to gut representatives but may be an ancestral characteristic of the whole order, and possibly also of ancient origin in the Euryarchaeota. 16S rRNA and McrA trees show that this new order of methanogens is very diverse and occupies environments highly relevant for methane production, therefore representing a key lineage to fully understand the diversity and evolution of methanogenesis.

Description: Streptococcus dysgalactiae subsp. equisimilis (SDSE) is an emerging human pathogen that causes life-threatening invasive infections such as streptococcal toxic shock syndrome. Recent epidemiological studies reveal that invasive SDSE infections have been increasing in Asia, Europe, and the United States. Almost all SDSE carry Lancefield group G or C antigen. We have determined the complete genome sequence of a human group C SDSE 167 strain. A comparison of its sequence with that of four SDSE strains, three in Lancefield group G and one in Lancefield group A, showed approximately 90% coverage. Most regions showing little or no homology were located in the prophages. There was no evidence of massive rearrangement in the genome of SDSE 167. Bayesian phylogeny using entire genome sequences showed that the most recent common ancestor of the five SDSE strains appeared 446 years ago. Interestingly, we found that SDSE 167 harbors sugar metabolizing enzymes in a unique region and streptodornase in the phage region, which presumably contribute to the degradation of host tissues and the prompted covRS mutation, respectively. A comparison of these five SDSE strains, which differ in Lancefield group antigens, revealed a gene cluster presumably responsible for the synthesis of the antigenic determinant. These results may provide the basis for molecular epidemiological research of SDSE.

Description: Mitochondria are intracellular organelles where oxidative phosphorylation is carried out to complete ATP synthesis. Mitochondria have their own genome; in metazoans, this is a small, circular molecule encoding 13 electron transport proteins, 22 tRNAs, and 2 rRNAs. In invertebrates, mitochondrial gene rearrangement is common, and it is correlated with increased substitution rates. In vertebrates, mitochondrial gene rearrangement is rare, and its relationship to substitution rate remains unexplored. Mitochondrial genes can also show spatial variation in substitution rates around the genome due to the mechanism of mtDNA replication, which produces a mutation gradient. To date, however, the strength of the mutation gradient and whether movement along the gradient in rearranged (or otherwise modified) genomes impacts genic substitution rates remain unexplored in the majority of vertebrates. Salamanders include both normal mitochondrial genomes and independently derived rearrangements and expansions, providing a rare opportunity to test the effects of large-scale changes to genome architecture on vertebrate mitochondrial gene sequence evolution. We show that: 1) rearranged/expanded genomes have higher substitution rates; 2) most genes in rearranged/expanded genomes maintain their position along the mutation gradient, substitution rates of the genes that do move are unaffected by their new position, and the gradient in salamanders is weak; and 3) genomic rearrangements/expansions occur independent of levels of selective constraint on genes. Together, our results demonstrate that large-scale changes to genome architecture impact mitochondrial gene evolution in predictable ways; however, despite these impacts, the same functional constraints act on mitochondrial protein-coding genes in both modified and normal genomes.

Description: Very little is known about genetic factors that regulate life history transitions during ontogeny. Closely related tiger salamanders ( Ambystoma species complex) show extreme variation in metamorphic timing, with some species foregoing metamorphosis altogether, an adaptive trait called paedomorphosis. Previous studies identified a major effect quantitative trait locus ( met1 ) for metamorphic timing and expression of paedomorphosis in hybrid crosses between the biphasic Eastern tiger salamander ( Ambystoma tigrinum tigrinum ) and the paedomorphic Mexican axolotl ( Ambystoma mexicanum ). We used existing hybrid mapping panels and a newly created hybrid cross to map the met1 genomic region and determine the effect of met1 on larval growth, metamorphic timing, and gene expression in the brain. We show that met1 maps to the position of a urodele-specific chromosome rearrangement on linkage group 2 that uniquely brought functionally associated genes into linkage. Furthermore, we found that more than 200 genes were differentially expressed during larval development as a function of met1 genotype. This list of differentially expressed genes is enriched for proteins that function in the mitochondria, providing evidence of a link between met1 , thyroid hormone signaling, and mitochondrial energetics associated with metamorphosis. Finally, we found that met1 significantly affected metamorphic timing in hybrids, but not early larval growth rate. Collectively, our results show that met1 regulates species and morph-specific patterns of brain transcription and life history variation .

Description: Oxidative phosphorylation (OXPHOS), the major energy-producing pathway in aerobic organisms, includes protein subunits encoded by both mitochondrial (mt) and nuclear (nu) genomes. How these independent genomes have coevolved is a long-standing question in evolutionary biology. Although mt genes evolve faster than most nu genes, maintenance of OXPHOS structural stability and functional efficiency may involve correlated evolution of mt and nu OXPHOS genes. The nu OXPHOS genes might be predicted to exhibit accelerated evolutionary rates to accommodate the elevated substitution rates of mt OXPHOS subunits with which they interact. Evolutionary rates of nu OXPHOS genes should, therefore, be higher than that of nu genes that are not involved in OXPHOS (nu non-OXPHOS). We tested the compensatory evolution hypothesis by comparing the evolutionary rates (synonymous substitution rate d S and nonsynonymous substitution rate d N ) among 13 mt OXPHOS genes, 60 nu OXPHOS genes, and 77 nu non-OXPHOS genes in vertebrates (7 fish and 40 mammal species). The results from a combined analysis of all OXPHOS subunits fit the predictions of the hypothesis. However, results from two OXPHOS complexes did not fit this pattern when analyzed separately. We found that the d N of nu OXPHOS genes for "core" subunits (those involved in the major catalytic activity) was lower than that of "noncore" subunits, whereas there was no significant difference in d N between genes for nu non-OXPHOS and core subunits. This latter finding suggests that compensatory changes play a minor role in the evolution of OXPHOS genes and that the observed accelerated nu substitution rates are due largely to reduced functional constraint on noncore subunits.

Description: Recently, we found dramatic mitochondrial DNA divergence of Israeli Chamaeleo chamaeleon populations into two geographically distinct groups. We aimed to examine whether the same pattern of divergence could be found in nuclear genes. However, no genomic resource is available for any chameleon species. Here we present the first chameleon transcriptome, obtained using deep sequencing (SOLiD). Our analysis identified 164,000 sequence contigs of which 19,000 yielded unique BlastX hits. To test the efficacy of our sequencing effort, we examined whether the chameleon and other available reptilian transcriptomes harbored complete sets of genes comprising known biochemical pathways, focusing on the nDNA-encoded oxidative phosphorylation (OXPHOS) genes as a model. As a reference for the screen, we used the human 86 (including isoforms) known structural nDNA-encoded OXPHOS subunits. Analysis of 34 publicly available vertebrate transcriptomes revealed orthologs for most human OXPHOS genes. However, OXPHOS subunit COX8 (Cytochrome C oxidase subunit 8), including all its known isoforms, was consistently absent in transcriptomes of iguanian lizards, implying loss of this subunit during the radiation of this suborder. The lack of COX8 in the suborder Iguania is intriguing, since it is important for cellular respiration and ATP production. Our sequencing effort added a new resource for comparative genomic studies, and shed new light on the evolutionary dynamics of the OXPHOS system.

Description: Human hemoglobins, the oxygen carriers in the blood, are composed by two α-like and two β-like globin monomers. The β-globin gene cluster located at 11p15.5 comprises one pseudogene and five genes whose expression undergoes two critical switches: the embryonic-to-fetal and fetal-to-adult transition. HBD encodes the -globin chain of the minor adult hemoglobin (HbA 2 ), which is assumed to be physiologically irrelevant. Paradoxically, reduced diversity levels have been reported for this gene. In this study, we sought a detailed portrait of the genetic variation within the β-globin cluster in a large human population panel from different geographic backgrounds. We resequenced the coding and noncoding regions of the two adult β-globin genes ( HBD and HBB ) in European and African populations, and analyzed the data from the β-globin cluster ( HBE , HBG2 , HBG1 , HBBP1 , HBD, and HBB ) in 1,092 individuals representing 14 populations sequenced as part of the 1000 Genomes Project. Additionally, we assessed the diversity levels in nonhuman primates using chimpanzee sequence data provided by the PanMap Project. Comprehensive analyses, based on classic neutrality tests, empirical and haplotype-based studies, revealed that HBD and its neighbor pseudogene HBBP1 have mainly evolved under purifying selection, suggesting that their roles are essential and nonredundant. Moreover, in the light of recent studies on the chromatin conformation of the β-globin cluster, we present evidence sustaining that the strong functional constraints underlying the decreased contemporary diversity at these two regions were not driven by protein function but instead are likely due to a regulatory role in ontogenic switches of gene expression.

Description: A tandem repeat’s (TR) propensity to mutate increases with repeat number, and can become very pronounced beyond a critical boundary, transforming it into a microsatellite (MS). However, a clear understanding of the mutational behavior of different TR classes and motifs and related mechanisms is lacking, as is a consensus on the existence of a boundary separating short TRs (STRs) from MSs. This hinders our understanding of MSs’ mutational properties and their effective use as genetic markers. Using indel calls for 179 individuals from 1000 Genomes Pilot-1 Project, we determined polymorphism incidence for four major TR classes, and formalized its varying relationship with repeat number using segmented regression. We observed a biphasic regime with a transition from a faster to a slower exponential growth at 9, 5, 4, and 4 repeats for mono-, di-, tri-, and tetranucleotide TRs, respectively. We used an in vitro mutagenesis assay to evaluate the contribution of strand slippage errors to mutability. STRs and MSs differ in their absolute polymorphism levels, but more importantly in their rates of mutability growth. Although strand slippage is a major factor driving mononucleotide polymorphism incidence, dinucleotide polymorphism incidence is greater than that expected due to strand slippage alone, indicating that additional cellular factors might be driving dinucleotide mutability in the human genome. Leveraging on hundreds of human genomes, we present the first comprehensive, genome-wide analysis of TR mutational behavior, encompassing several motif sizes and compositions.

Description: Eukaryotic genome sequencing projects often yield bacterial DNA sequences, data typically considered as microbial contamination. However, these sequences may also indicate either symbiont genes or lateral gene transfer (LGT) to host genomes. These bacterial sequences can provide clues about eukaryote–microbe interactions. Here, we used the genome of the primitive animal Trichoplax adhaerens (Metazoa: Placozoa), which is known to harbor an uncharacterized Gram-negative endosymbiont, to search for the presence of bacterial DNA sequences. Bioinformatic and phylogenomic analyses of extracted data from the genome assembly (181 bacterial coding sequences [CDS]) and trace read archive (16S rDNA) revealed a dominant proteobacterial profile strongly skewed to Rickettsiales ( Alphaproteobacteria ) genomes. By way of phylogenetic analysis of 16S rDNA and 113 proteins conserved across proteobacterial genomes, as well as identification of 27 rickettsial signature genes, we propose a Rickettsiales endosymbiont of T. adhaerens (RETA). The majority (93%) of the identified bacterial CDS belongs to small scaffolds containing prokaryotic-like genes; however, 12 CDS were identified on large scaffolds comprised of eukaryotic-like genes, suggesting that T . adhaerens might have recently acquired bacterial genes. These putative LGTs may coincide with the placozoan’s aquatic niche and symbiosis with RETA. This work underscores the rich, and relatively untapped, resource of eukaryotic genome projects for harboring data pertinent to host–microbial interactions. The nature of unknown (or poorly characterized) bacterial species may only emerge via analysis of host genome sequencing projects, particularly if these species are resistant to cell culturing, as are many obligate intracellular microbes. Our work provides methodological insight for such an approach.

Description: Following polyploidy, duplicate genes are often deleted, and if they are not, then duplicate regulatory regions are sometimes lost. By what mechanism is this loss and what is the chance that such a loss removes function? To explore these questions, we followed individual Arabidopsis thaliana–A. thaliana conserved noncoding sequences (CNSs) into the Brassica ancestor, through a paleohexaploidy and into Brassica rapa . Thus, a single Brassicaceae CNS has six potential orthologous positions in B. rapa ; a single Arabidopsis CNS has three potential homeologous positions. We reasoned that a CNS, if present on a singlet Brassica gene, would be unlikely to lose function compared with a more redundant CNS, and this is the case. Redundant CNSs go nondetectable often. Using this logic, each mechanism of CNS loss was assigned a metric of functionality. By definition, proved deletions do not function as sequence. Our results indicated that CNSs that go nondetectable by base substitution or large insertion are almost certainly still functional (redundancy does not matter much to their detectability frequency), whereas those lost by inferred deletion or indels are approximately 75% likely to be nonfunctional. Overall, an average nondetectable, once-redundant CNS more than 30 bp in length has a 72% chance of being nonfunctional, and that makes sense because 97% of them sort to a molecular mechanism with "deletion" in its description, but base substitutions do cause loss. Similarly, proved-functional G-boxes go undetectable by deletion 82% of the time. Fractionation mutagenesis is a procedure that uses polyploidy as a mutagenic agent to genetically alter RNA expression profiles, and then to construct testable hypotheses as to the function of the lost regulatory site. We show fractionation mutagenesis to be a "deletion machine" in the Brassica lineage.

Description: Understanding the molecular basis of within and between species phenotypic variation is one of the main goals of Biology. In Drosophila , most of the work regarding this issue has been performed in D. melanogaster , but other distantly related species must also be studied to verify the generality of the findings obtained for this species. Here, we make the case for D. americana , a species of the virilis group of Drosophila that has been diverging from the model species, D. melanogaster , for approximately 40 Myr. To determine the suitability of this species for such studies, polymorphism and recombination estimates are presented for D. americana based on the largest nucleotide sequence polymorphism data set so far analyzed (more than 100 data sets) for this species. The polymorphism estimates are also compared with those obtained from the comparison of the genome assembly of two D. americana strains (H5 and W11) here reported. As an example of the general utility of these resources, we perform a preliminary study on the molecular basis of lifespan differences in D. americana . First, we show that there are lifespan differences between D. americana populations from different regions of the distribution range. Then, we perform five F2 association experiments using markers for 21 candidate genes previously identified in D. melanogaster . Significant associations are found between polymorphism at two genes ( hep and Lim3 ) and lifespan. For the F2 association study involving the two sequenced strains (H5 and W11), we identify amino acid differences at Lim3 and Hep that could be responsible for the observed changes in lifespan. For both genes, no large gene expression differences were observed between the two strains.

Description: Comparative mitochondrial genomics of arbuscular mycorrhizal fungi (AMF) provide new avenues to overcome long-lasting obstacles that have hampered studies aimed at understanding the community structure, diversity, and evolution of these multinucleated and genetically polymorphic organisms. AMF mitochondrial (mt) genomes are homogeneous within isolates, and their intergenic regions harbor numerous mobile elements that have rapidly diverged, including homing endonuclease genes, small inverted repeats, and plasmid-related DNA polymerase genes ( dpo ), making them suitable targets for the development of reliable strain-specific markers. However, these elements may also lead to genome rearrangements through homologous recombination, although this has never previously been reported in this group of obligate symbiotic fungi. To investigate whether such rearrangements are present and caused by mobile elements in AMF, the mitochondrial genomes from two Glomeraceae members (i.e., Glomus cerebriforme and Glomus sp. ) with substantial mtDNA synteny divergence, were sequenced and compared with available glomeromycotan mitochondrial genomes . We used an extensive nucleotide/protein similarity network-based approach to investigate dpo diversity in AMF as well as in other organisms for which sequences are publicly available. We provide strong evidence of dpo -induced inter-haplotype recombination, leading to a reshuffled mitochondrial genome in Glomus sp. These findings raise questions as to whether AMF single spore cultivations artificially underestimate mtDNA genetic diversity. We assessed potential dpo dispersal mechanisms in AMF and inferred a robust phylogenetic relationship with plant mitochondrial plasmids. Along with other indirect evidence, our analyses indicate that members of the Glomeromycota phylum are potential donors of mitochondrial plasmids to plants.

Description: The study of genetic and phenotypic variation is fundamental for understanding the dynamics of bacterial genome evolution and untangling the evolution and epidemiology of bacterial pathogens. Neisseria meningitidis ( Nm ) is among the most intriguing bacterial pathogens in genomic studies due to its dynamic population structure and complex forms of pathogenicity. Extensive genomic variation within identical clonal complexes (CCs) in Nm has been recently reported and suggested to be the result of homologous recombination, but the extent to which recombination contributes to genomic variation within identical CCs has remained unclear. In this study, we sequenced two Nm strains of identical serogroup (C) and multi-locus sequence type (ST60), and conducted a systematic analysis with an additional 34 Nm genomes. Our results revealed that all gene content variation between the two ST60 genomes was introduced by homologous recombination at the conserved flanking genes, and 94.25% or more of sequence divergence was caused by homologous recombination. Recombination was found in genes associated with virulence factors, antigenic outer membrane proteins, and vaccine targets, suggesting an important role of homologous recombination in rapidly altering the pathogenicity and antigenicity of Nm . Recombination was also evident in genes of the restriction and modification systems, which may undermine barriers to DNA exchange. In conclusion, homologous recombination can drive both gene content variation and sequence divergence in Nm . These findings shed new light on the understanding of the rapid pathoadaptive evolution of Nm and other recombinogenic bacterial pathogens.

Description: Organelle DNA is no stranger to palindromic repeats. But never has a mitochondrial or plastid genome been described in which every coding region is part of a distinct palindromic unit. While sequencing the mitochondrial DNA of the nonphotosynthetic green alga Polytomella magna , we uncovered precisely this type of genic arrangement. The P. magna mitochondrial genome is linear and made up entirely of palindromes, each containing 1–7 unique coding regions. Consequently, every gene in the genome is duplicated and in an inverted orientation relative to its partner. And when these palindromic genes are folded into putative stem-loops, their predicted translational start sites are often positioned in the apex of the loop. Gel electrophoresis results support the linear, 28-kb monomeric conformation of the P. magna mitochondrial genome. Analyses of other Polytomella taxa suggest that palindromic mitochondrial genes were present in the ancestor of the Polytomella lineage and lost or retained to various degrees in extant species. The possible origins and consequences of this bizarre genomic architecture are discussed.

Description: Expression quantitative trait loci (eQTLs) have been found to be enriched in trait-associated single-nucleotide polymorphisms (SNPs). However, whether eQTLs are adaptive to different environmental factors and its relative evolutionary significance compared with nonsynonymous SNPs (NS SNPs) are still elusive. Compiling environmental correlation data from three studies for more than 500,000 SNPs and 42 environmental factors, including climate, subsistence, pathogens, and dietary patterns, we performed a systematic examination of the adaptive patterns of eQTLs to local environment. Compared with intergenic SNPs, eQTLs are significantly enriched in the lower tail of a transformed rank statistic in the environmental correlation analysis, indicating possible adaptation of eQTLs to the majority of 42 environmental factors. The mean enrichment of eQTLs across 42 environmental factors is as great as, if not greater than, that of NS SNPs. The enrichment of eQTLs, although significant across all levels of recombination rate, is inversely correlated with recombination rate, suggesting the presence of selective sweep or background selection. Further pathway enrichment analysis identified a number of pathways with possible environmental adaption from eQTLs. These pathways are mostly related with immune function and metabolism. Our results indicate that eQTLs might have played an important role in recent and ongoing human adaptation and are of special importance for some environmental factors and biological pathways.

Description: Environmental or geological changes can create new niches that drive ecological species divergence without the immediate cessation of gene flow. However, few such cases have been characterized. On a recently formed volcano, Mt. Etna, Senecio aethnensis and S. chrysanthemifolius inhabit contrasting environments of high and low altitude, respectively. They have very distinct phenotypes, despite hybridizing promiscuously, and thus may represent an important example of ecological speciation "in action," possibly as a response to the rapid geological changes that Mt. Etna has recently undergone. To elucidate the species’ evolutionary history, and help establish the species as a study system for speciation genomics, we sequenced the transcriptomes of the two Etnean species, and the outgroup, S. vernalis , using Illumina sequencing. Despite the species’ substantial phenotypic divergence, synonymous divergence between the high- and low-altitude species was low (d S = 0.016 ± 0.017 [SD]). A comparison of species divergence models with and without gene flow provided unequivocal support in favor of the former and demonstrated a recent time of species divergence (153,080 ya ± 11,470 [SE]) that coincides with the growth of Mt. Etna to the altitudes that separate the species today. Analysis of d N /d S revealed wide variation in selective constraint between genes, and evidence that highly expressed genes, more "multifunctional" genes, and those with more paralogs were under elevated purifying selection. Taken together, these results are consistent with a model of ecological speciation, potentially as a response to the emergence of a new, high-altitude niche as the volcano grew.

Description: Admixture mapping has been enormously resourceful in identifying genetic variations linked to phenotypes, adaptation, and diseases. In this study through analysis of copy number variable regions (CNVRs), we report extensive restructuring in the genomes of the recently admixed African-Indian population (OG-W-IP) that inhabits a highly saline environment in Western India. The study included subjects from OG-W-IP (OG), five different Indian and three HapMap populations that were genotyped using Affymetrix version 6.0 arrays. Copy number variations (CNVs) detected using Birdsuite were used to define CNVRs. Population structure with respect to CNVRs was delineated using random forest approach. OG genomes have a surprising excess of CNVs in comparison to other studied populations. Individual ancestry proportions computed using STRUCTURE also reveals a unique genetic component in OGs. Population structure analysis with CNV genotypes indicates OG to be distant from both the African and Indian ancestral populations. Interestingly, it shows genetic proximity with respect to CNVs to only one Indian population IE-W-LP4, which also happens to reside in the same geographical region. We also observe a significant enrichment of molecular processes related to ion binding and receptor activity in genes encompassing OG-specific CNVRs. Our results suggest that retention of CNVRs from ancestral natives and de novo acquisition of CNVRs could accelerate the process of adaptation especially in an extreme environment. Additionally, this population would be enormously useful for dissecting genes and delineating the involvement of CNVs in salt adaptation.

Description: As an important subtype of structural variations, chromosomal translocation is associated with various diseases, especially cancers, by disrupting gene structures and functions. Traditional methods for identifying translocations are time consuming and have limited resolutions. Recently, a few studies have employed next-generation sequencing (NGS) technology for characterizing chromosomal translocations on human genome, obtaining high-throughput results with high resolutions. However, these studies are mainly focused on mechanism-specific or site-specific translocation mapping. In this study, we conducted a comprehensive genome-wide analysis on the characterization of human chromosomal material exchange with regard to the chromosome translocations. Using NGS data of 1,481 subjects from the 1000 Genomes Project, we identified 15,349,092 translocated DNA fragment pairs, ranging from 65 to 1,886 bp and with an average size of approximately 102 bp. On average, each individual genome carried about 10,364 pairs, covering approximately 0.069% of the genome. We identified 16 translocation hot regions, among which two regions did not contain repetitive fragments. Results of our study overlapped with a majority of previous results, containing approximately 79% of approximately 2,340 translocations characterized in three available translocation databases. In addition, our study identified five novel potential recurrent chromosomal material exchange regions with greater than 20% detection rates. Our results will be helpful for an accurate characterization of translocations in human genomes, and contribute as a resource for future studies of the roles of translocations in human disease etiology and mechanisms.

Description: Angiosperm mitochondrial genomes exhibit many unusual properties, including heterogeneous nucleotide composition and exceptionally large and variable genome sizes. Determining the role of nonadaptive mechanisms such as mutation bias in shaping the molecular evolution of these unique genomes has proven challenging because their dynamic structures generally prevent identification of homologous intergenic sequences for comparative analyses. Here, we report an analysis of angiosperm mitochondrial DNA sequences that are derived from inserted plastid DNA ( mtpts ). The availability of numerous completely sequenced plastid genomes allows us to infer the evolutionary history of these insertions, including the specific nucleotide substitutions and indels that have occurred because their incorporation into the mitochondrial genome. Our analysis confirmed that many mtpts have a complex history, including frequent gene conversion and multiple examples of horizontal transfer between divergent angiosperm lineages. Nevertheless, it is clear that the majority of extant mtpt sequence in angiosperms is the product of recent transfer (or gene conversion) and is subject to rapid loss/deterioration, suggesting that most mtpts are evolving relatively free from functional constraint. The evolution of mtpt sequences reveals a pattern of biased mutational input in angiosperm mitochondrial genomes, including an excess of small deletions over insertions and a skew toward nucleotide substitutions that increase AT content. However, these mutation biases are far weaker than have been observed in many other cellular genomes, providing insight into some of the notable features of angiosperm mitochondrial architecture, including the retention of large intergenic regions and the relatively neutral GC content found in these regions.

Description: Neotropical primates (NP) are presently distributed in the New World from Mexico to northern Argentina, comprising three large families, Cebidae, Atelidae, and Pitheciidae, consequently to their diversification following their separation from Old World anthropoids near the Eocene/Oligocene boundary, some 40 Ma. The evolution of NP has been intensively investigated in the last decade by studies focusing on their phylogeny and timescale. However, despite major efforts, the phylogenetic relationship between these three major clades and the age of their last common ancestor are still controversial because these inferences were based on limited numbers of loci and dating analyses that did not consider the evolutionary variation associated with the distribution of gene trees within the proposed phylogenies. We show, by multispecies coalescent analyses of selected genome segments, spanning along 92,496,904 bp that the early diversification of extant NP was marked by a 2-fold increase of their effective population size and that Atelids and Cebids are more closely related respective to Pitheciids. The molecular phylogeny of NP has been difficult to solve because of population-level phenomena at the early evolution of the lineage. The association of evolutionary variation with the distribution of gene trees within proposed phylogenies is crucial for distinguishing the mean genetic divergence between species (the mean coalescent time between loci) from speciation time. This approach, based on extensive genomic data provided by new generation DNA sequencing, provides more accurate reconstructions of phylogenies and timescales for all organisms.

Description: Dogs shared a much closer relationship with humans than any other domesticated animals, probably due to their unique social cognitive capabilities, which were hypothesized to be a by-product of selection for tameness toward humans. Here, we demonstrate that genes involved in glutamate metabolism, which account partially for fear response, indeed show the greatest population differentiation by whole-genome comparison of dogs and wolves. However, the changing direction of their expression supports a role in increasing excitatory synaptic plasticity in dogs rather than reducing fear response. Because synaptic plasticity are widely believed to be cellular correlates of learning and memory, this change may alter the learning and memory abilities of ancient scavenging wolves, weaken the fear reaction toward humans, and prompt the initial interspecific contact.

Description: Hirsutella minnesotensis [Ophiocordycipitaceae (Hypocreales, Ascomycota)] is a dominant endoparasitic fungus by using conidia that adhere to and penetrate the secondary stage juveniles of soybean cyst nematode. Its genome was de novo sequenced and compared with five entomopathogenic fungi in the Hypocreales and three nematode-trapping fungi in the Orbiliales (Ascomycota). The genome of H. minnesotensis is 51.4 Mb and encodes 12,702 genes enriched with transposable elements up to 32%. Phylogenomic analysis revealed that H. minnesotensis was diverged from entomopathogenic fungi in Hypocreales. Genome of H. minnesotensis is similar to those of entomopathogenic fungi to have fewer genes encoding lectins for adhesion and glycoside hydrolases for cellulose degradation, but is different from those of nematode-trapping fungi to possess more genes for protein degradation, signal transduction, and secondary metabolism. Those results indicate that H. minnesotensis has evolved different mechanism for nematode endoparasitism compared with nematode-trapping fungi. Transcriptomics analyses for the time-scale parasitism revealed the upregulations of lectins, secreted proteases and the genes for biosynthesis of secondary metabolites that could be putatively involved in host surface adhesion, cuticle degradation, and host manipulation. Genome and transcriptome analyses provided comprehensive understanding of the evolution and lifestyle of nematode endoparasitism.

Description: Anopheles gambiae is a major mosquito vector of malaria in Africa. Although increased use of insecticide-based vector control tools has decreased malaria transmission, elimination is likely to require novel genetic control strategies. It can be argued that the absence of an A. gambiae inbred line has slowed progress toward genetic vector control . In order to empower genetic studies and enable precise and reproducible experimentation, we set out to create an inbred line of this species. We found that amenability to inbreeding varied between populations of A. gambiae . After full-sib inbreeding for ten generations, we genotyped 112 individuals—56 saved prior to inbreeding and 56 collected after inbreeding—at a genome-wide panel of single nucleotide polymorphisms (SNPs). Although inbreeding dramatically reduced diversity across much of the genome, we discovered numerous, discrete genomic blocks that maintained high heterozygosity. For one large genomic region, we were able to definitively show that high diversity is due to the persistent polymorphism of a chromosomal inversion. Inbred lines in other eukaryotes often exhibit a qualitatively similar retention of polymorphism when typed at a small number of markers. Our whole-genome SNP data provide the first strong, empirical evidence supporting associative overdominance as the mechanism maintaining higher than expected diversity in inbred lines. Although creation of A. gambiae lines devoid of nearly all polymorphism may not be feasible, our results provide critical insights into how more fully isogenic lines can be created.

Description: Greenstreet, S. P. R., Fraser, H. M., Rogers, S. I., Trenkel, V. M., Simpson, S. D., and Pinnegar, J. K. 2012. Redundancy in metrics describing the composition, structure, and functioning of the North Sea demersal fish community. – ICES Journal of Marine Science, 69: 8–22. Broader ecosystem management objectives for North Sea demersal fish currently focus on restoring community size structure. However, most policy drivers explicitly concentrate on restoring and conserving biodiversity, and it has not yet been established that simply restoring demersal fish size composition will be sufficient to reverse declines in biodiversity and ensure a generally healthy community. If different aspects of community composition, structure, and function vary independently, then to monitor all aspects of community general health will require application of a suite of metrics. This assumes low redundancy among the metrics used in any such suite and implies that addressing biodiversity issues specifically will require explicit management objectives for particular biodiversity metrics. This issue of metric redundancy is addressed, and 15 metrics covering five main attributes of community composition, structure, and function are applied to groundfish survey data. Factor analysis suggested a new interpretation of the metric information and indicated that a minimum suite of seven metrics was necessary to ensure that all changes in the general health of the North Sea demersal fish community were monitored properly. Covariance among size-based and species-diversity metrics was low, implying that restoration of community size structure would not necessarily reverse declines in species diversity.

Description: Dankel, D. J., Aps, R., Padda, G., Röckmann, C., van der Sluijs, J. P., Wilson, D. C., and Degnbol, P. 2012. Advice under uncertainty in the marine system. – ICES Journal of Marine Science, 69: 3–7. There is some uncertainty in the fisheries science–policy interface. Although progress has been made towards more transparency and participation in fisheries science in ICES Areas, routine use of state-of-the-art quantitative and qualitative tools to address uncertainty systematically is still lacking. Fisheries science that gives advice to policy-making is plagued by uncertainties; the stakes of the policies are high and value-laden and need therefore to be treated as an example of "post-normal science" (PNS). To achieve robust governance, understanding of the characteristics and implications of the scientific uncertainties for management strategies need to come to the centre of the table. This can be achieved using state-of-the-art tools such as pedigree matrices and uncertainty matrices, as developed by PNS scholars and used in similar science–policy arenas on other complex issues. An explicit extension of the peer community within maritime systems will be required to put these new tools in place. These new competences become even more important as many countries within the ICES Area are now embarking on new policies.

Description: Browman, H.I. 2012. Quo Vadimus. – ICES Journal of Marine Science, 69: 1–2.

Description: Hüssy, K., Coad, J. O., Farrell, E. D., Clausen, L. A. W., and Clarke, M. W. 2012. Age verification of boarfish ( Capros aper ) in the Northeast Atlantic. – ICES Journal of Marine Science, 69: 34–40. The boarfish ( Capros aper ) is a pelagic species of recent interest to the fishing industry, with landings increasing by 〉500% over the past 3 years. The objective of the study was to provide a method for age determination based on whole sagittal otoliths, with the results to be used in stock assessment. Translucent zones laid down at regular intervals are identified by marginal increment analysis as seasonally recurring. Translucent zones are formed between September/October and March/April, regardless of fish age. The occurrence of the first annulus is validated by analysis of presumed daily growth increments. Subsequent annulus deposition is homogenous between individuals and allows general guidelines to be derived for interpreting the age of boarfish using their otoliths.

Description: Kim, D., Yang, E. J., Kim, K. H., Shin, C-W., Park, J., Yoo, S., and Hyun, J-H. 2012. Impact of an anticyclonic eddy on the summer nutrient and chlorophyll a distributions in the Ulleung Basin, East Sea (Japan Sea). – ICES Journal of Marine Science, 69: 23–29. The impact of the anticyclonic Ulleung Warm Eddy (UWE) on the vertical distributions of nutrient and chlorophyll a (Chl a ) concentrations in the Ulleung Basin (UB) was investigated during the contrasting summers of 2005 and 2007. The physical structure of the water column was characterized by an intrathermocline eddy (ITE) in 2005, whereas the UWE remained distant from the sampling transect in 2007. Water column structures appeared to be highly stratified, and nutrients in the surface waters were totally depleted at all stations. In 2005, an exceptionally high concentration of Chl a (5.5 mg m –3 ) was measured below the surface mixed layer in the eddy core (station D3), and values of ~2.5 mg m –3 were observed at the eddy edge (stations D2 and D4). Formation of an ITE efficiently mixed surface and deep-ocean waters, the latter supplying sufficient nutrients to generate an extremely high concentration of Chl a at the base of the subsurface layer. Overall, the results indicated that the anticyclonic UWE plays a key ecological role in supporting substantial phytoplankton biomass in the nutrient-depleted surface waters in summer and maintaining high benthic mineralization in the deep-sea sediments of the UB.

Description: Franchini, P., Sola, L., Crosetti, D., Milana, V., and Rossi, A. R. 2012. Low levels of population genetic structure in the gilthead sea bream, Sparus aurata , along the coast of Italy. – ICES Journal of Marine Science, 69: 41–50. The gilthead sea bream, Sparus aurata , is a coastal, commercially important fish. Contrasting results concerning the genetic structure of the species at different geographic scales have been reported. Here, an investigation is made into the population genetic structure of S. aurata along the coast of Italy, using samples analysed previously and material from new sampling sites (12) and using different microsatellite loci (10). One sample from the eastern Atlantic and three temporal replicates from one site were also included. The presence of a weak (overall F ST = 0.0072), but significant, genetic population subdivision was detected by F -statistics. Temporal replicates indicate genetic data consistency over time. Isolation by distance between the Atlantic and the coast of Italy is suggested by a Mantel test. The distributional pattern of genetic variance obtained by analysis of molecular variation reflects the geographic sampling areas, but is only partially congruent with the results obtained with fewer sites and loci. The dispersal of passive eggs/larvae by the main currents appears to contribute to shaping the gene flow. Given the intensity of sea bream aquaculture activities in Italy, the possibility that aquaculture may have partially contributed to the population genetic pattern detected cannot be excluded.

Description: Charrier, F., Mazel, V., Caraguel, J-M., Abdallah, Y., Le Gurun, L. L., Legault, A., and Laffaille, P. 2012. Escapement of silver-phase European eels, Anguilla anguilla , determined from fishing activities in a Mediterranean lagoon (Or, France). – ICES Journal of Marine Science, 69: 30–33. Escapement of silver eels from a Mediterranean lagoon was estimated by a capture–tagging–recapture and automated tag-reading study. The population of silver-phase eels in the lagoon was estimated to be 13.2 kg ha –1 , with an escapement rate from the commercial fishery of 76.8%.

Description: Jennings, S., and Lee, J. 2012. Defining fishing grounds with vessel monitoring system data. – ICES Journal of Marine Science, 69: 51–63. Methods for defining fishing grounds to support marine spatial planning and management are developed, applied, and compared. The methods are broadly applicable and repeatable because they use vessel monitoring system (VMS) data that are archived and increasingly accessible. For several fleets at regional and national scales, an attempt is made to assess how the choice of criteria for defining grounds influences (i) size, shape, and location, (ii) overlap among grounds, and (iii) the extent to which annual and multi-annual patterns of fishing activity describe grounds used seasonally or by individual vessels. The results show that grounds defined by excluding infrequently fished margins (areas with 〈10% of total fishing activity) are typically 50% smaller than total fished area. However, landings weight or value (LWV) per unit activity can be higher at the margins, with 10% of activity usually accounting for 10–20% of LWV. The removal of fishing activity in the margins, as a consequence of regulation or fleet behaviour, would lead to disproportionately greater reductions in interactions with other fisheries, sectors, and the environment. Accessible high-resolution information on the "anatomy" of all fishing grounds would better inform debates on the allocation and the use of marine space and the integration of fisheries and environmental management.

Description: Bentley, N., Kendrick, T. H., Starr, P. J., and Breen, P. A. 2012. Influence plots and metrics: tools for better understanding fisheries catch-per-unit-effort standardizations. – ICES Journal of Marine Science, 69: 84–88. Standardization of catch per unit effort using generalized linear models (GLMs) is a common procedure that attempts to remove the confounding effects of variables other than abundance. Simple plots and metrics are described to assist understanding the standardization effects of explanatory variables included in GLMs, illustrated with an example based on New Zealand trevally ( Caranx lutescens ) data.

Description: Aldrin, M., Mortensen, B., Storvik, G., Nedreaas, K., Aglen, A., and Aanes, S. 2012. Improving management decisions by predicting fish bycatch in the Barents Sea shrimp fishery. – ICES Journal of Marine Science, 69: 64–74. When the bycatch of juvenile fish within the Barents Sea shrimp fishery is too large, the area is closed to fishing for a certain period. Bycatch is estimated from sampled trawl hauls, for which the shrimp yield is recorded, along with the total number of various bycatch fish species. At present, bycatch estimation is based on a simple estimator, the sum of the number of fish caught within the area of interest within a small time window, divided by the corresponding shrimp yield (in weight). No historical data are used. A model-based estimation is proposed in which spatio-temporal models are constructed for the variation in both the yield of shrimp and the amount of bycatch in space and time. The main effects are described through generalized additive models, and local dependence structures are specified through correlated random effects. Model estimation includes historical and recent data. Experiments with both simulated and real data show that the model-based estimator outperforms the present simple estimator when a low or moderate number of samples (e.g. 〈20) is available, whereas the two estimators are equally good when the number of samples is high.

Description: Jiao, Y., Smith, E. P., O'Reilly, R., and Orth, D. J. 2012. Modelling non-stationary natural mortality in catch-at-age models. – ICES Journal of Marine Science, 69: 105–118. Non-stationarity in the natural mortality of many species has been discussed among research scientists, but no generally applicable models/methods have been developed where a statistical catch-at-age (SCA) model framework is used. Using the Atlantic weakfish ( Cynoscion regalis ) fishery as an example, several SCA models are developed to assess the population dynamics, then compared. Models used included (i) an SCA with constant natural mortality, (ii) an SCA with unknown natural mortality, but with a hierarchical prior distribution from a mixture of distributions based on life-history information, (iii) an SCA with age-specific unknown natural mortality, (iv) an SCA with time-varying natural mortality, following a random-walk process, and (v) an SCA with age-specific time-varying natural mortality. The last two models imply that the population dynamics are not stationary. A Bayesian approach was used to estimate parameters, and performance of the models was compared through goodness-of-fit and the retrospective patterns of the models. A simulation study was then conducted to test the uncertainty resulting from model selection, the efficiency of using the best model selected based on deviance information criterion, and whether changes in natural mortality could be detected. An SCA with time-varying natural mortality, following a random-walk process, is proposed for the example fishery here. The estimated non-stationary temporal patterns in natural mortality were linked further to climate–ocean oscillation indices, to diagnose possible mechanisms/linkages to these temporal variations in natural mortality.

Description: Juntunen, T., Vanhatalo, J., Peltonen, H., and Mäntyniemi, S. 2012. Bayesian spatial multispecies modelling to assess pelagic fish stocks from acoustic- and trawl-survey data. – ICES Journal of Marine Science, 69: 95–104. A Bayesian spatial model was constructed to estimate the abundance of multiple fish species in a pelagic environment. Acoustic- and trawl-survey data were combined with environmental data to predict the spatial distribution of (i) the acoustic backscattering of fish, (ii) the relative proportion of each species, and (iii) their mean length in the Gulf of Finland in the northeastern Baltic Sea. By combining the three spatial model layers, the spatial distribution of the biomass of each species was estimated. The model consists of a linear predictor on environmental variables and a spatial random effect given by a Gaussian process. A Bayesian approach is a natural choice for the task because it provides a theoretically justified means of summarizing the uncertainties from various model layers. In the study area, three species dominate pelagic waters: sprat ( Sprattus sprattus ), herring ( Clupea harengus ), and three-spined stickleback ( Gasterosteus aculeatus ). Results are presented for each model layer and for estimated total biomass for each species in 2 x 2 km lattices. The posterior mean and central 95% credible intervals of total biomass were sprat 45.7 kt (27.7–71.6), herring 24.6 kt (9.7–41.3), and three-spined stickleback 1.9 kt (0.9–3.2).

Description: Smith, W. K., and Solow, A. R. 2012. Missing and presumed lost: extinction in the ocean and its inference. – ICES Journal of Marine Science, 69: 89–94. The number of modern extinctions in the ocean is unknown. The actual demise of the last individual of a species is essentially unobservable, so extinction can only be inferred. Statistical methods are described for inferring extinction from sighting records, species–area considerations, and taxonomic samples collected at two different times. The methods are illustrated using a variety of real datasets, including a sighting record of the Caribbean monk seal and results from three surveys of benthic invertebrates.

Description: Delavenne, J., Metcalfe, K., Smith, R. J., Vaz, S., Martin, C. S., Dupuis, L., Coppin, F., and Carpentier, A. 2012. Systematic conservation planning in the eastern English Channel: comparing the Marxan and Zonation decision-support tools. – ICES Journal of Marine Science, 69: 75–83. The systematic conservation approach is now commonly used for the design of efficient marine protected area (MPA) networks, and identifying these priority areas often involves using specific conservation-planning software. Several such software programmes have been developed in recent years, each differing in the underlying algorithms used. Here, an investigation is made into whether the choice of software influences the location of priority areas by comparing outputs from Marxan and Zonation, two widely used conservation-planning, decision-support tools. Using biological and socio-economic data from the eastern English Channel, outputs are compared and it is shown that the two software packages identified similar sets of priority areas, although the relatively wide distribution of habitat types and species considered offered much flexibility. Moreover, the similarity increased with increasing spatial constraint, especially when using real-world cost data, suggesting that the choice of cost metric has a greater influence on conservation-planning analyses than the choice of software. However, Marxan generally produced more efficient results and Zonation produced results with greater connectivity, so the most appropriate software package will depend on the overall goals of the MPA planning process.

Description: Andersen, B. S., Ulrich, C., Eigaard, O. R., and Christensen, A-S. 2012. Short-term choice behaviour in a mixed fishery: investigating métier selection in the Danish gillnet fishery. – ICES Journal of Marine Science, 69: 131–143. The study presents a short-term effort allocation modelling approach based on a discrete choice random utility model combined with a survey questionnaire to examine the selection of métiers (a combination of fishing area and target species) in the Danish North Sea gillnet fishery. Key decision variables were identified from the survey questionnaire, and relevant proxies for the decision function were identified based on available landings and effort information. Additional variables from the survey questionnaire were further used to validate and verify the outcome of the choice model. Commercial fishers in a mixed fishery make use of a number of decision variables used previously in the literature, but also a number of decision parameters rarely explicitly accounted for, such as price, weather, and management regulation. The seasonal availability of individual target species and within-year changes in monthly catch ration were the main explanatory drivers, but gillnetters were also responsive to information on the whole fishery, fish prices, and distance travelled to fishing grounds. Heterogeneous responses were evident from geographic differences in home harbour, which underpins the need to understand alternative fishing strategies among individual gillnetters better.

Description: Calise, L., and Knutsen, T. 2012. Multifrequency target strength of northern krill ( Meganyctiphanes norvegica ) swimming horizontally. – ICES Journal of Marine Science, 69: 119–130. Multifrequency acoustic measurements on ex situ horizontally swimming krill were made in a novel experimental setting. An ensemble of northern krill ( Meganyctiphanes norvegica ) was introduced to a large enclosure (a mesocosm), and acoustic backscatter was sampled using a multifrequency (70, 120, and 200 kHz) echosounder (Simrad EK60). Two submerged lamps were placed at opposite sides of the mesocosm and switched on and off to induce the krill, by light attraction, to swim horizontally through the acoustic beams. By tracking echoes, animal displacement, swimming speed, and target strength ( TS ) by frequency were estimated. The dominant and secondary modes of the total-length distribution were 21.8 ± 3.0 and 27.8 ± 2.7 mm, respectively. Although krill orientation was assumed stable and the ping rate was high, the range and inter-ping variability of the average TS values were large, decreasing and increasing with frequency, respectively. The overall TS frequency response observed and concurrent measurements at 120 and 200 kHz confirm the theoretical expectation that the acoustic backscatter from the investigated organisms were confined to the Rayleigh and Geometric scattering regions, a finding that might both aid acoustic identification and size-group separation of in situ northern krill.

Description: Rosen, S., Engås, A., Fernö, A., and Jörgensen, T. 2012. The reactions of shoaling adult cod to a pelagic trawl: implications for commercial trawling. – ICES Journal of Marine Science, 69: 303–312. The reactions of shoaling adult Atlantic cod to a pelagic trawl were measured during fishing off the north coast of Norway. Cod remaining in the trawl track dived at rates as fast as 0.35 m s –1 following vessel passage and swam away from the vessel, in the direction of the approaching trawl, at an average rate of 0.6 m s –1 . They did not attempt to swim ahead of the trawl as documented previously, but passed into the lower half of the trawl entrance and swam slowly in the direction of trawling at a rate of 0.2–0.5 m s –1 as the trawl's greater speed through water carried them deeper into the trawl. Shoals compressed vertically once inside the trawl, suggesting that packing density increased at least fourfold. Fish remained in the lower part of the trawl as they moved through its tapered section towards the codend, with little to no clearance above the bottom panel, but significant clearance beneath the top panel. Catches were sufficient to support commercial harvest, and the behaviour observed suggests that changes in trawl design and fishing strategy might improve fuel economy and species selectivity.

Description: Shackell, N. L., Bundy, A., Nye, J. A., and Link, J. S. 2012. Common large-scale responses to climate and fishing across Northwest Atlantic ecosystems. – ICES Journal of Marine Science, 69: 151–152. Investigating whether there were common biological responses to climate and fishing across seven Northwest Atlantic ecosystems, a minimum/maximum autocorrelation factor analysis of biological indicators for each region revealed a common primary multivariate trend of a rapid change during the 1980s and early 1990s. There was a strong common pattern in the biological indicators responsible for the primary multivariate temporal trend in the five more northerly regions: an increase in the abundance of phytoplankton, an increase in biomass at mid-trophic levels, and a decline in predatory groundfish size. The common associations between patterns and drivers were fishing indices and the Atlantic Multidecadal Oscillation, but all associations weakened when co-varying drivers were held constant. The results are consistent with known long-term effects of intense fishing, such as a decline in average fish size and changes in species composition. Less fishing pressure has allowed some regions to recover to former predatory biomass levels since the late 1990s, but the bulk of the biomass consists of fewer species. However, fishing was not the only driver, and a more mechanistic understanding of how the climate affects lower trophic levels is needed to contextualize climate effects in heavily fished ecosystems.

Description: Deporte, N., Ulrich, C., Mahévas, S., Demanèche, S., and Bastardie, F. 2012. Regional métier definition: a comparative investigation of statistical methods using a workflow applied to international otter trawl fisheries in the North Sea. – ICES Journal of Marine Science, 69: 331–342. The European Common Fisheries Policy recognizes the importance of accounting for heterogeneity in fishing practices, and métier-based sampling is now at the core of the EU Data Collection Framework. The implementation of such an approach would require Member States to agree on the standard regional métier definitions and on practical rules to categorize logbook records into métiers. Several alternative approaches have been used in the past to categorize landings profiles, but no consensus has yet emerged. A generic open-source workflow is developed to test and compare a selection of methods, including principal components analysis (PCA), hierarchical agglomerative clustering (HAC), K -means, and Clustering LARge Applications (CLARA), and to provide simple allocation rules. This workflow is applied to a unique regional dataset consisting of bottom-trawl logbooks of five North Sea countries. No method proved to be infallible, but combining PCA with either CLARA or HAC performed best. For 2008, a hierarchical classification with 14 species assemblages is proposed. Discriminant analysis proved more robust than simple ordination methods for allocating a new logbook record into an existing métier. The whole approach is directly operational and could contribute to defining more objective and consistent métiers across European fisheries.

Description: Greenstreet, S. P. R., Rogers, S. I., Rice, J. C., Piet, G. J., Guirey, E. J., Fraser, H. M., and Fryer, R. J. 2012. A reassessment of trends in the North Sea Large Fish Indicator and a re-evaluation of earlier conclusions. – ICES Journal of Marine Science, 69: 343–345. Previous analysis of the Large Fish Indicator, the basis for the North Sea "fish communities" EcoQO, suggested a strong recovery between 2001 and 2008. Discovery of a mistake in this earlier analysis now suggests that this recovery is not as strong as first thought, but reanalysis of the corrected data still supports the majority of conclusions drawn in the original paper.

Description: Xie, Y., Michielsens, C. G. J., and Martens, F. J. 2012. Classification of fish and non-fish acoustic tracks using discriminant function analysis. – ICES Journal of Marine Science, 69: 313–322. Hydroacoustic data acquired for estimating fish populations contain information on both fish and non-fish targets, so sonar technicians traditionally rely on their knowledge of fish behaviour and experience with hydroacoustics to remove non-fish targets from the hydroacoustic data. This process is often labour-intensive and time-consuming, making real-time assessment of fish populations difficult. Simple solutions are not always available for all circumstances. However, the split-beam sonar data collected in the lower Fraser River, British Columbia, Canada, showed distinct signatures between actively swimming fish and non-fish objects such as drifting debris, surface bubbles, and stationary objects in the water column and off the river bottom. Acoustic tracks of fish and non-fish targets were characterized by differentiable statistical patterns that were amenable to discriminant function analysis (DFA). An application of DFA to segregate fish and non-fish targets detected by a split-beam sonar system in the lower Fraser River is presented, characteristics of user-identified fish and non-fish acoustic tracks being utilized as learning samples for the DFA. Also, a method to rank the discriminating power of individual variables is presented, providing guidance for constructing efficient and effective discriminant functions with variables that offer high discriminating power. The DFA yielded classification accuracies of 96% for fish and 91% for non-fish tracks and reduced the manual sorting time by 50–75%.

Description: Rey, F. 2012. Declining silicate concentrations in the Norwegian and Barents Seas. – ICES Journal of Marine Science, 69: 208–212. Since 1990, a decline in silicate concentrations together with increasing salinities has been observed in the Atlantic water of the Norwegian and Barents Seas. This decline in silicate has been found to be related to the relative proportion in which eastern and western source water masses from the northeastern North Atlantic enter the Norwegian Sea.

Description: Zhou, S., Milton, D. A., and Fry, G. C. 2012. Integrated risk analysis for rare marine species impacted by fishing: sustainability assessment and population trend modelling. – ICES Journal of Marine Science, 69: 271–280. Risk assessment and conservation management of rare species are challenging due to a lack of data. We developed an integrated risk assessment approach to assess human impact on population sustainability of rare species. The approach involved two components: a quantitative sustainability assessment coupled with modelling trends in relative abundance. Both components took nil catches into account through zero-inflated statistical distributions that simultaneously modelled the zero and non-zero catches separately in submodels. The sustainability assessment used detection–non-detection data for population estimation and linked sustainability to easily collected life-history traits. This component provides an assessment of population sustainability at one point in time. The trend modelling applied zero-inflated negative binomial models to temporal trends in density and dispersion of species. It provided a complement to the static sustainability assessment. We applied this integrated approach to assess the risk to 14 species of rare, protected sea snakes incidentally caught in the Australian Northern Prawn Fishery. This approach can be applicable for risk assessment of many species with limited abundance data, a large number of absences and some presence–absence information only.

Description: Frid, C. L. J., and Paramor, O. A. L. 2012. Feeding the world: what role for fisheries? – ICES Journal of Marine Science, 69: 145–150. Fisheries (wild capture and aquaculture) deliver more than 110 million tonnes of food and around 15% of the dietary protein to the 7 billion people currently living on the planet. With the global population expected to peak at 9 billion by 2050, and 〉80% of global fish stocks currently fully or overexploited (and aquaculture is at least in part dependent on capture fisheries), the contribution of fisheries looks set to decline. The challenge is therefore determining how better management, an ecosystem perspective, and more efficient utilization of fisheries waste can support fisheries products continuing to contribute significantly to "feeding the world" up to and beyond the population peak.

Description: Boström, M. K., Östman, Ö., Bergenius, M. A. J., and Lunneryd, S-G. 2012. Cormorant diet in relation to temporal changes in fish communities. – ICES Journal of Marine Science, 69: 175–183. The stomach contents of 229 great cormorants ( Phalacrocorax carbo sinensis ) shot between March and October 2009 along the Swedish east coast were analysed for differences in diet between gender, age, and breeding phase. Sticklebacks ( Gasterosteus aculeatus and Pungitius pungitius ) were the most common prey, followed by eelpout ( Zoarces viviparus ), herring ( Clupea harengus membras ), and cyprinids (Cyprinidae). Diet did not differ between age and gender, but changed over the breeding season. The different phases explained around 10% of the total variation in stomach content between cormorants, suggesting no major shift in diet over the breeding season. The diet of cormorants in 2009 was compared with the results of a study conducted in the same area in 1992. There were evident changes in the diet between 1992 and 2009, with less perch ( Perca fluviatilis ) and cyprinids and more eelpout and herring in 2009. This change in diet could partly be related to changes in the fish community. The seasonal changes in diet composition of whole stomachs were less notable than in many previous studies, but long-term changes in the fish communities may induce changes in cormorant diet. It is clearly important to use stomach contents in areas with many small fish species for a comprehensive assessment of cormorant diet.

Description: Cañás, L., Stransky, C., Schlickeisen, J., Sampedro, M. P., and Fariña, A. C. 2012. Use of the otolith shape analysis in stock identification of anglerfish ( Lophius piscatorius ) in the Northeast Atlantic. – ICES Journal of Marine Science, 69: 250–256. The geographical variability in the anglerfish Lophius piscatorius otolith shape was investigated to examine biological evidence for the current separation into Northern and Southern stocks of L. piscatorius established in the Southwest European waters (Southern Shelf area) by the International Council for the Exploration of the Sea (ICES). In all, 629 otoliths of L. piscatorius (12–90 cm total length) from commercial fishing and research surveys were collected in two consecutive years, 2008 and 2009. The otolith shape was compared between stocks for differences in univariate descriptors using non-parametric analysis of covariance, as well as elliptical Fourier shape descriptors, through a metric multidimensional scaling and a linear discriminant analysis (LDA). The differences found in the otolith shape between Northern and Southern stocks were not large enough to hold the current status of L. piscatorius regarding stock definition in the Southern Shelf area, although homogeneous mixing between both stocks is not expected as there were significant differences between three univariate descriptors and the LDA was significant. These results are in agreement with the previous L. piscatorius stock identification studies (genetics, body morphology, and migratory behaviour), which do not support the current separation of these stocks but found substructures within the area.

Description: The kallikrein ( KLK ) gene family comprises the largest uninterrupted locus of serine proteases in the human genome and represents a notable case for studying the evolutionary fate of duplicated genes. In primates, a recent duplication event gave rise to KLK2 and KLK3 , both encoding essential proteins for the cascade of seminal plasma liquefaction. We reconstructed the evolutionary history of KLK2 and KLK3 by comparative analysis of the orthologous sequences from 22 primate species, calculated d N / d S ratios, and addressed the hypothesis of coevolution with their substrates, the semenogelins (SEMG1 and SEMG2). Our findings support the placement of the KLK2–KLK3 duplication in the Catarrhini ancestor and unveil the frequent loss of KLK2 throughout primate evolution by different genomic mechanisms, including unequal crossing-over, deletions, and pseudogenization. We provide evidences for an adaptive evolution of KLK3 toward an expanded enzymatic spectrum, with an effect on the hydrolysis of semen coagulum. Furthermore, we found associations between mating system, the number of SEMG repeat units, and the number of functional KLK2 and KLK3 , suggesting complex evolutionary dynamics shaped by reproductive biology.

Description: The question of Jewish ancestry has been the subject of controversy for over two centuries and has yet to be resolved. The "Rhineland hypothesis" depicts Eastern European Jews as a "population isolate" that emerged from a small group of German Jews who migrated eastward and expanded rapidly. Alternatively, the "Khazarian hypothesis" suggests that Eastern European Jews descended from the Khazars, an amalgam of Turkic clans that settled the Caucasus in the early centuries CE and converted to Judaism in the 8th century. Mesopotamian and Greco–Roman Jews continuously reinforced the Judaized empire until the 13th century. Following the collapse of their empire, the Judeo–Khazars fled to Eastern Europe. The rise of European Jewry is therefore explained by the contribution of the Judeo–Khazars. Thus far, however, the Khazars’ contribution has been estimated only empirically, as the absence of genome-wide data from Caucasus populations precluded testing the Khazarian hypothesis. Recent sequencing of modern Caucasus populations prompted us to revisit the Khazarian hypothesis and compare it with the Rhineland hypothesis. We applied a wide range of population genetic analyses to compare these two hypotheses. Our findings support the Khazarian hypothesis and portray the European Jewish genome as a mosaic of Near Eastern-Caucasus, European, and Semitic ancestries, thereby consolidating previous contradictory reports of Jewish ancestry. We further describe a major difference among Caucasus populations explained by the early presence of Judeans in the Southern and Central Caucasus. Our results have important implications for the demographic forces that shaped the genetic diversity in the Caucasus and for medical studies.

Description: The impact of transposable elements (TEs) on genome structure, plasticity, and evolution is still not well understood. The recent availability of complete genome sequences makes it possible to get new insights on the evolutionary dynamics of TEs from the phylogenetic analysis of their multiple copies in a wide range of species. However, this source of information is not always fully exploited. Here, we show how the history of transposition activity may be qualitatively and quantitatively reconstructed by considering the distribution of transposition events in the phylogenetic tree, along with the tree topology. Using statistical models developed to infer speciation and extinction rates in species phylogenies, we demonstrate that it is possible to estimate the past transposition rate of a TE family, as well as how this rate varies with time. This methodological framework may not only facilitate the interpretation of genomic data, but also serve as a basis to develop new theoretical and statistical models.

Description: We generated a genome-wide replication profile in the genome of Lachancea kluyveri and assessed the relationship between replication and base composition . This species diverged from Saccharomyces cerevisiae before the ancestral whole genome duplication . The genome comprises eight chromosomes among which a chromosomal arm of 1 Mb has a G + C-content much higher than the rest of the genome. We identified 252 active replication origins in L. kluyveri and found considerable divergence in origin location with S. cerevisiae and with Lachancea waltii. Although some global features of S. cerevisiae replication are conserved: Centromeres replicate early, whereas telomeres replicate late, we found that replication origins both in L. kluyveri and L. waltii do not behave as evolutionary fragile sites. In L. kluyveri, replication timing along chromosomes alternates between regions of early and late activating origins, except for the 1 Mb GC-rich chromosomal arm. This chromosomal arm contains an origin consensus motif different from other chromosomes and is replicated early during S-phase. We showed that precocious replication results from the specific absence of late firing origins in this chromosomal arm. In addition, we found a correlation between GC-content and distance from replication origins as well as a lack of replication-associated compositional skew between leading and lagging strands specifically in this GC-rich chromosomal arm. These findings suggest that the unusual base composition in the genome of L. kluyveri could be linked to replication.

Description: The vomeronasal organ (VNO) is an olfactory structure that detects pheromones and environmental cues. It consists of sensory neurons that express evolutionary unrelated groups of transmembrane chemoreceptors. The predominant V1R and V2R receptor repertoires are believed to detect airborne and water-soluble molecules, respectively. It has been suggested that the shift in habitat of early tetrapods from water to land is reflected by an increase in the ratio of V1R/V2R genes. Snakes, which have a very large VNO associated with a sophisticated tongue delivery system, are missing from this analysis. Here, we use RNA-seq and RNA in situ hybridization to study the diversity, evolution, and expression pattern of the corn snake vomeronasal receptor repertoires. Our analyses indicate that snakes and lizards retain an extremely limited number of V1R genes but exhibit a large number of V2R genes, including multiple lineages of reptile-specific and snake-specific expansions. We finally show that the peculiar bigenic pattern of V2R vomeronasal receptor gene transcription observed in mammals is conserved in squamate reptiles, hinting at an important but unknown functional role played by this expression strategy. Our results do not support the hypothesis that the shift to a vomeronasal receptor repertoire dominated by V1Rs in mammals reflects the evolutionary transition of early tetrapods from water to land. This study sheds light on the evolutionary dynamics of the vomeronasal receptor families in vertebrates and reveals how mammals and squamates differentially adapted the same ancestral vomeronasal repertoire to succeed in a terrestrial environment.

Description: Bachiller, E., Cotano, U., Boyra, G., and Irigoien, X. 2013. Spatial distribution of the stomach weights of juvenile anchovy ( Engraulis encrasicolus L.) in the Bay of Biscay. – ICES Journal of Marine Science, 70: 362–378. Previous studies have shown that the survival of larval and juvenile anchovy off the continental shelf in the Bay of Biscay is not significantly different from that observed on the shelf, even though the food concentration is significantly higher on the shelf. In this paper we investigate the causes for the lack of relationship between food and survival for anchovy juvenile through analyses of feeding activity versus zooplankton distribution (in the transition from summer to fall from 2003 to 2010). The spatial distribution of the stomach weights of juvenile anchovy in relation to body size on and off the shelf revealed higher stomach fullness in areas off the shelf, where less zooplankton biomass was available. This result indicates that the food concentration is not always the main factor determining ingestion in fish. A situation of comparatively lesser prey abundance may provide better feeding conditions when combined with lower predation risk and lower light attenuation. In addition, the relatively heavier stomachs found in anchovies caught in years of higher recruitment indices suggest a large stomach content may be a symptom of good biological condition that may favour the winter survival of a larger fraction of the juvenile stock.

Description: Van Hoof, L. 2013. Design or pragmatic evolution: applying ITQs in EU fisheries management. – ICES Journal of Marine Science, 70: 462–470. Among the proposals for the 2012 revision of the EU Common Fisheries Policy, a strong case is made for the introduction of a system of rights-based management. The EU perceives individual fishing concessions as an important instrument for capacity management. We will use the introduction of individual tradable quotas in the management of the Dutch North Sea beam trawl fisheries as a case for exploring the effect of the introduction of such an instrument. The effect will be assessed in terms of reduction of fishing capacity in the Dutch beam trawl fleet and its economic and social impact. These Dutch experiences will be translated to the current debate on the reform of the EU Common Fisheries Policy. Especially, we will focus on the issues of "relative stability", the concentration of rights, and the effects on the small-scale fisheries sector. Some of the negative effects associated with individual tradable rights can be addressed through design. However, trying to maintain stability and counter perceived negative impacts on fishing communities will modify the effect of introducing individual fishing concessions.

Description: Spitz, J., Chouvelon, T., Cardinaud, M., Kostecki, C., and Lorance, P. 2013. Prey preferences of adult sea bass Dicentrarchus labrax in the northeastern Atlantic: implications for bycatch of common dolphin Delphinus delphis – ICES Journal of Marine Science, 70: 452–461. In the northeastern Atlantic, adult sea bass ( Dicentrarchus labrax ) is one of largest fish living on the shelf, and this species has important commercial value. However, pelagic trawl fisheries that target sea bass have negative operational interactions with common dolphins ( Delphinus delphis ). Our goal was to determine the diet of adult sea bass in the Bay of Biscay from stomach-content and stable-isotope analyses, and explore the dietary overlap between sea bass and common dolphins. We found that sea bass primarily target small pelagic fish, most notably mackerel ( Scomber scombrus ), scads ( Trachurus spp.), anchovy ( Engraulis encrasicolus ), and sardine ( Sardina pilchardus ). These four species also dominated the diets of common dolphins. This overlap in feeding preferences could increase the risk of dolphins being caught by trawl fisheries while feeding among sea bass, and may be an underlying mechanism to explain the high rate of common dolphin bycatch observed in the pelagic trawl fishery for sea bass in the Bay of Biscay. Understanding the foraging ecology and trophic interactions of predator species is an essential step for identifying and resolving management issues in the northeastern Atlantic and other marine ecosystems.

Description: Lassen, H., Pedersen, S. A., Frost, H., and Hoff. A. 2013. Fishery management advice with ecosystem considerations. – ICES Journal of Marine Science, 70: 471–479. The European Union Marine Strategy Framework Directive (MSFD) includes four descriptors of Good Environmental Status (GES) which are affected by fishing activity. These descriptors are: biodiversity, fish stocks, foodweb, and seabed integrity. This paper shows how these descriptors can be related to variables within an ecological model and how an ecological model can be used to analyse whether the fishing pressure that is estimated based on bioeconomic criteria is within general sustainable limits. The paper presents an example of such an analysis of the Eastern Baltic cod fishery using two models: a bioeconomic model and an ecological model. The models are calibrated based on historic data. The mapping between the descriptors specified by MSFD and variables available for analysis in the models is incomplete, e.g. genetics and spatial structures are not included in the models. The models can be used strategically, providing a qualitative understanding of the anticipated relative changes.

Description: Finley, C. and Oreskes, N. 2013. Maximum sustained yield: a policy disguised as science. – ICES Journal of Marine Science, 70: 245–250. Overfishing is most commonly explained as an example of the tragedy of the commons, where individuals are unable to control their activities, leading to the destruction of the resource they are dependent on. The historical record suggests otherwise. Between1949 and 1958, the US State Department used fisheries science, and especially the concept of maximum sustained yield (MSY) as a political tool to achieve its foreign policy objectives. During the Cold War, the Department thought that if countries were allowed to restrict fishing in their waters, it might lead to restrictions on passage of military vessels. While there has been much criticism of MSY and its failure to conserve fish stocks, there has been little attention paid to the political context in which MSY was adopted.

Description: The sharing of the Northeast Atlantic mackerel ( Scomber scombrus ) stock is analysed as a game between four parties: the European Union, Norway, the Faroe Islands, and Iceland. Consideration is given to how the outcome depends on the nature of the stock's migrations. Two types of migrations are considered: (i) density-dependent, where the mackerel migrates into the Icelandic economic zone only if it exceeds 3.5 million t, and (ii) stochastic migrations, where the said migrations are stochastic. It is determined that the Faroe Islands would never accept a cooperative solution wherein they can only fish with the globally optimal fishing mortality within their own zone. This is also true for Iceland when the migrations into her zone are stochastic, but not if they are density-dependent. In the latter case, the other players have incentives to retaliate to Icelandic overfishing by fishing harder, which greatly reduces the number of years when mackerel are available in the Icelandic zone. It is assumed that the objective is maximization of the catch volume over a time-horizon of 50 years.

Description: Harris, P. T., Bridge, T. C. L., Beaman, R. J., Webster, J. M., Nichol, S. L., and Brooke, B. P. 2013. Submerged banks in the Great Barrier Reef, Australia, greatly increase available coral reef habitat. – ICES Journal of Marine Science, 70: 284–293. Anthropogenic global ocean warming is predicted to cause bleaching of many near-sea-surface (NSS) coral reefs, placing increased importance on deeper reef habitats to maintain coral reef biodiversity and ecosystem function. However, the location and spatial extent of many deep reef habitats is poorly known. The question arises: how common are deep reef habitats in comparison with NSS reefs? We used a dataset from the Great Barrier Reef (GBR) to show that only about 39% of available seabed on submerged banks is capped by NSS coral reefs (16 110 km 2 ); the other 61% of bank area (25 600 km 2 ) is submerged at a mean depth of around 27 m and represents potential deep reef habitat that is spatially distributed along the GBR continental shelf in the same latitudinal distribution as NSS reefs. Out of 25 600 km 2 of submerged bank area, predictive habitat modelling indicates that more than half (around 14 000 km 2 ) is suitable habitat for coral communities.

Description: Dutertre, M., Hamon D., Chevalier C., and Ehrhold, A. 2013. The use of the relationships between environmental factors and benthic macrofaunal distribution in the establishment of a baseline for coastal management – ICES Journal of Marine Science, 70: 294–308. Relationships between benthic macrofauna and natural abiotic factors were studied along the coastal fringe of South Brittany, situated north of the Gulf of Biscay on the French Atlantic continental shelf. Within the framework of the REBENT network, sediment characteristics, depth, and macrofaunal abundance were determined for 95 stations spread over five subtidal sectors, using a combination of seabed acoustic remote sensing systems and grab sampling. The physico-chemical properties of the water column and the hydrodynamic conditions were generated by validated three-dimensional environmental models which take into account variations over shorter temporal scales. Multivariate analyses ranked 16 natural abiotic variables according to the significance of their influence on the macrofauna. Together these variables explained 51% of spatial variation in the macrofauna, with morpho-sedimentological and hydrological factors contributing 22% and 26%, respectively. The outputs from validated three-dimensional environmental models appear to be useful interpretational tools for benthic ecology studies, especially in estuarine and coastal ecosystems with high environmental variability due to regular freshwater inputs. Ten major species assemblages were identified using biological and physical characteristics. The results provide important baseline knowledge for future ecosystem and resource management.

Description: Prigge, E., Marohn, L., Oeberst R., and Hanel, R. 2013. Model prediction vs. reality—testing the predictions of a European eel ( Anguilla anguilla ) stock dynamics model against the in situ observation of silver eel escapement in compliance with the European eel regulation – ICES Journal of Marine Science, 70: 309–318. A direct monitoring of European silver eel ( Anguilla anguilla , L) escapement from rivers and estuaries has been proven to be challenging, and a Europe-wide documentation of escaping silver eel numbers therefore hardly seems realistic. To reinforce management decisions, policy-makers are thus widely reliant on the accuracy of escapement models. A 3-year programme of silver eel escapement monitoring was undertaken to compile model input data and revise an eel population model (German Eel Model II; GEM II) already used in the decision-making process of management authorities. By compiling necessary input data and analysing vital system-specific population characteristics, it was possible to compare the documented silver eel escapement with the modelled potential silver eel escapement. Resulting model predictions were close to actually monitored escapement numbers, which were distinctly lower than reference escapement values for the same freshwater system given in the implementation report of the German Eel Management Plans. Applying different commercial and recreational catch scenarios revealed the sensitivity of the model. The results show the potential of the GEM II and highlight the importance of high-quality input data to use model predictions as the basis for management measures.

Description: Sumaila, U. R. 2013. How to make progress in disciplining overfishing subsidies. ICES Journal of Marine Science, 70: 251–258. The World Trade Organization (WTO) has been working for more than seven years now to discipline overfishing subsidies, as mandated by the global community, without success. I argue that this failure is partly because WTO negotiators aim for an all-inclusive deal, i.e. negotiations are conducted as a "single undertaking", whereby results must be achieved in all areas. Negotiators are required to broker an all-inclusive deal for all maritime WTO member countries and for all fisheries, whether domestic or international; small or large scale; developing or developed country fisheries. It is argued here that this commitment to a "single undertaking" does not align the incentive to remove subsidies with national interests, and therefore needs to be changed by splitting the world's fisheries into domestic and international fisheries. In this way, the battle for eliminating overfishing subsidies for some stocks would shift to home countries, and for others this would still rest with the international community. This split, it is argued, would align the incentives and improve the chances of eliminating overfishing subsidies.

Description: Morgan, M. J., Garabana, D., Rideout, R. M., Román, E., Pérez-Rodriguez, A. and Saborido-Rey, F. 2013. Changes in distribution of Greenland halibut in a varying environment. – ICES Journal of Marine Science, 70: 352–361. Fish are expected to respond to changing oceanographic temperature by altering their distribution. Off the coast of Newfoundland, Canada, there have been major changes in oceanographic temperature over the last several decades, with both record cold and record warm years being observed. Greenland halibut is an important flatfish species in the area, and is distributed in deep waters over a very wide geographic range. Thus, it might be buffered from temperature change in the overall area by reduced temperature variation at depth, and the diverse temperature conditions over its wide range. We examined intrapopulation variation in temperature and depth distribution, and the biological changes in relation to changes in available temperature. On the Flemish Cap, variation in available temperature was limited, and changes in depth were related to changing age composition and the differential depth distribution with age/size. In other areas there was a larger decline in available temperature, and associated with this, Greenland halibut moved to deeper waters and occupied warmer temperatures than they had previously. Concurrently, growth declined and condition increased. This study shows that shifts in distribution may not result in maintenance of homogeneous environmental conditions, and that resulting biological changes will be difficult to predict.

Description: Hufnagl, M., Huebert, K. B., and Temming, A. 2013. How does seasonal variability in growth, recruitment, and mortality affect the performance of length-based mortality and asymptotic length estimates in aquatic resources? – ICES Journal of Marine Science, 70: 329–341. We tested the sensitivity of eight methods for estimating total mortality from size frequencies (modified Wetherall; Powell; Beverton and Holt; Jones and van Zalinge; Hoenig; Ssentongo and Larkin; seasonal and non-seasonal Length Converted Catch Curve) to violations of basic assumptions, such as seasonal growth, mortality, recruitment and variable asymptotic length L or growth parameter K . For each method, bias was estimated by simulating length frequency distributions with different combinations of known L , Z and K values, calculating ( Z / K ) and L estimates, and comparing the true input with the estimated output values. Input mortality was generally underestimated by all methods and in 27% of all simulations no method provided estimates within ± 1. Spring recruitment especially negatively influenced the mortality estimate. A decision tree was developed that provides general guidance in selecting appropriate methods despite violated assumptions, but species-specific case studies are recommended. An example of a species-specific study is provided for the brown shrimp, Crangon crangon. Despite inherent limitations for all methods, the results illustrate that estimates of and Z for brown shrimp can be improved substantially by selecting suitable methods and correcting for observed bias.

Description: Brunel, T., and Piet, G. J. 2013. Is age structure a relevant criterion for the health of fish stocks? – ICES Journal of Marine Science, 70: 270–283. The age and size structure of exploited fish stocks is one of the criteria for Good Environmental Status of commercial fish. However, two underlying assumptions to this criterion remain to be tested: first, that a well-balanced age structure is indeed indicative of a "healthier" stock, and second, that managers can exert a control on the age structure, independently from the regulation of stock abundance. This study investigates these questions using simulations from a population model in which recruitment is based on egg production, which is more sensitive to age structure variations than the traditionally used spawning stock biomass (SSB) and that was parameterized to represent the population dynamics of North Sea cod, plaice, and herring. Our results show that (i) the age structure is highly dependent on the selection pattern, as well as on the level of fishing mortality; (ii) the selection pattern determines the ability of fish stocks to withstand, and recover from, external perturbation; and (iii) the selection pattern determines the output of the fishery providing the management option to balance stable but relatively low yields vs strongly fluctuating high yields. Therefore, we propose to make the selection pattern, for which clear management targets can be set, a policy goal instead of the age structure that is currently in place.

Description: Bréchon, A. L., Coombs S. H., Sims D. W., and Griffiths A. M. 2013. Development of a rapid genetic technique for the identification of clupeid larvae in the Western English Channel and investigation of mislabelling in processed fish products – ICES Journal of Marine Science, 70: 399–407. A novel genetic technique, involving real-time polymerase chain reaction (PCR) and high-resolution melt (HRM) analysis was developed, which successfully identified four species of clupeid larvae from the Western English Channel and revealed the presence of herring ( Clupea harengus ). This was unexpected in a context of global warming as herring are at the southern limit of their range off southwest England, where the larvae are generally very rare. The application of the method in cases of degraded DNA, such as in processed fishery products, was also assessed, revealing a low level of misidentification or equivocal identity (8.3%). Eight-year-old formalin-preserved samples were also analysed successfully. This highlights the potentially broad applications for HRM-based approaches to species identification in marine science.

Description: Quetglas, A., Ordines, F., Hidalgo, M., Monserrat, S., Ruiz, S., Amores, Á., Moranta, J., and Massutí, E. 2013. Synchronous combined effects of fishing and climate within a demersal community. – ICES Journal of Marine Science, 70: 319–328. Accumulating evidence shows that fishing exploitation and environmental variables can synergistically affect the population dynamics of exploited populations. Here, we document an interaction between fishing impact and climate variability that triggered a synchronic response in the population fluctuations of six exploited species in the Mediterranean from 1965–2008. Throughout this period, the fishing activity experienced a sharp increase in fishing effort, which caused all stocks to shift from an early period of underexploitation to a later period of overexploitation. This change altered the population resilience of the stocks and brought about an increase in the sensitivity of its dynamics to climate variability. Landings increased exponentially when underexploited but displayed an oscillatory behaviour once overexploited. Climatic indices, related to the Mediterranean mesoscale hydrography and large-scale north Atlantic climatic variability, seemed to affect the species with broader age structure and longer lifespan, while the global-scale El Niño Southern Oscillation index (ENSO) positively influenced the population abundances of species with a narrow age structure and short lifespan. The species affected by ENSO preferentially inhabit the continental shelf, suggesting that Mediterranean shelf ecosystems are sensitive to the hydroclimatic variability linked to global climate.

Description: Papetti, C., Di Franco, A., Zane, L., Guidetti, P., De Simone, V., Spizzotin, M., Zorica, B., Cikes Kec, V. and Mazzoldi, C. 2013. Single population and common natal origin for Adriatic Scomber scombrus stocks: evidence from an integrated approach – ICES Journal of Marine Science, 70: 387–398. In order to implement proper fishery management strategies aimed at avoiding stock declines, information about connectivity among stocks and populations is critically required. In this perspective, the present study investigated population structure of the Atlantic mackerel Scomber scombrus in the northern–central Adriatic Sea by integrating multiple approaches (analysis of fisheries data, population genetics, and otolith chemistry). Monthly data of fishery landings indicate a latitudinal trend along the western Adriatic coast, with Atlantic mackerel disappearing from the northern waters in winter, corresponding to the reproductive season. Population genetic analyses by genotyping of eight microsatellites clearly point to the presence of a single panmictic population in the northern–central Adriatic Sea. Otolith cores of samples from the northern–central Adriatic were chemically homogeneous, suggesting a common spawning ground. These results strongly suggest that Atlantic mackerel perform an autumn–winter migration in the northern–central Adriatic Sea, from the northern to the central sector, to reach a single spawning ground, and that a single population is present in this area. Considering that S. scombrus has shown a marked decline in the last 40 years in the Adriatic, this study highlights a potential high vulnerability to collapse by overfishing for the Atlantic mackerel stocks in this geographic area.

Description: Uusitalo, L., Fleming-Lehtinen, V., Hällfors, H., Jaanus A., Hällfors, S., and London, L. 2013. A novel approach for estimating phytoplankton biodiversity – ICES Journal of Marine Science, 70: 408–417. Maintaining biodiversity is one of the main priorities in environmental protection. The biodiversity of phytoplankton, the key primary producers in the marine ecosystem, is, however, often very difficult to estimate, since the phytoplankton assemblage includes a vast number of taxa, many of which occur in such small quantities that they may not be recorded in routine sampling. Moreover, many taxa cannot be identified to species level by routine methods such as light microscopy of preserved samples, even by a skilled taxonomist. This means that, in general, we cannot assume to have a complete list of species in the ecosystem at any given point in time. In this paper, we present an approach for evaluating phytoplankton biodiversity in spite of this challenge. Since eutrophication, which increases phytoplankton biomass, has been identified as the most important factor causing degradation of the Baltic Sea ecosystem, the proposed approach was evaluated against total phytoplankton biomass. When phytoplankton biomass was low, both low and high biodiversity values were observed, and, as the phytoplankton biomass increased, the high biodiversity values disappeared. These results were consistent both using data based on individual samples and using yearly sampling station averages.

Description: Frie, A. K., Hammill, M. O., Hauksson, E., Lind, Y., Lockyer, C., Stenman, O., and Svetocheva, O. 2013. Error patterns in age estimation and tooth readability assignment of grey seals ( Halichoerus grypus ): results from a transatlantic, image-based, blind-reading study using known-age animals – ICES Journal of Marine Science, 70: 418–430. We analysed error patterns in a first interlaboratory grey seal ( Halichoerus grypus ) age-reading experiment. The experiment involved ten readers, who estimated age using images of cementum growth layers from teeth of 68 known-age seals (0–22 years). The percentages of correct estimates ranged from 32.4% to 60.3% among readers, and 89.3% of all errors were by ±1–2 years. Six readers showed increasing underageing with increasing seal age. An elevated risk of underestimation by 1 year occurred in teeth collected 0–5 months after breeding and was attributed to more frequent absence of a distinct growth layer for the new year and lack of information on months between the last birthday and the date of sample collection (plusmonths). For plusmonths 6–11, positive bias was predominant, suggesting that overestimation is the more common error when plusmonth information is available. Readers assigned readability scores to the tooth sections, and 79.1% of all ageing errors occurred in sections of low or intermediate readability. Excluding these sections would, however, also exclude 43.0% of all correct estimates. Neither levels of age estimation error nor predictive values of readability assignments were associated with reader experience levels. Analyses of image markings identified common errors in delineations of annual increment layers.

Description: Domínguez-Petit, R., Ouellet, P., and Lambert, Y. 2013. Reproductive strategy, egg characteristics and embryonic development of Greenland halibut ( Reinhardtius hippoglossoides ). – ICES Journal of Marine Science, 70: 342–351. Despite the commercial importance of Greenland halibut (GH), important gaps exist in our knowledge of the reproductive and early life stage biology for this species. The present study examined through laboratory experiments the spawning strategy, realized fecundity, egg characteristics, biochemical composition, and embryonic development of GH. The results confirmed the hypothesis that GH is a single-batch spawner producing large eggs, resulting in low realized fecundity. Embryonic development and hatching time are highly dependent on incubation temperature; 50% hatching occurred after 46, 30, and 24 days at 2, 4, and 6°C, respectively. Few changes in the biochemical composition of the eggs are observed during embryonic development. Newly hatched larvae are not well developed, having a large yolk sac, no pigmentation and incomplete development of the jaws. Egg specific density confirmed the mesopelagic distribution of the eggs at sea. However, important buoyancy changes occurring in the last 3–4 days before hatching indicate that larvae hatch higher in the water column. These results are important for understanding advection and dispersion processes of GH eggs and larvae and the connectivity between spawning grounds and nursery areas.

Description: Johnsen, E., and Harbitz, A. 2013. Small-scale spatial structuring of burrowed sandeels and the catching properties of the dredge – ICES Journal of Marine Science, 70: 379–386. The lesser sandeel ( Ammodytes marinus ) remains burrowed for long periods in the seabed in sandy areas, and its habitat preference is reflected in a large-scale (10 km) patchy distribution in the North Sea. The small-scale spatial structure of sandeels in the seabed is less known. In this study, it was found that sandeel catches in grabs and dredges fit gamma distributions with a common shape parameter on both a 10 m and a 100 m scale. These spatial models were used to examine the catching properties of the dredge, which, in the new ICES sandeel assessment, is used to estimate recruitment and maturity curves. Analyses of the empirical data and simulations of the dredge and grab catch rates show very low dredge catchability (~5%), but no length-dependent catchability in the dredge was found. However, analyses of the catch rates of parallel dredge hauls indicated increasing dredge catchability with sandeel density. The observed small-scale patchiness and the low and seemingly density-dependent dredge catchability result in biased abundance estimates with low precision. Therefore, we recommend a joint international effort to develop a more catch-efficient sampling dredge, where the effect of density on catchability is further examined.

Description: Rousi, H., Laine, A. O., Peltonen, H., Kangas, P., Andersin, A-B., Rissanen, J., Sandberg-Kilpi, E., and Bonsdorff, E. 2013. Long-term changes in coastal zoobenthos in the northern Baltic Sea: the role of abiotic environmental factors – ICES Journal of Marine Science, 70: 440–451. We investigated site-specific changes in a coastal zoobenthic community during 1964–2007 in the northern Baltic Sea. Multivariate analysis indicated that the community structure had changed. The amphipods Monoporeia affinis and Pontoporeia femorata decreased drastically in the late 1970s–early 1980s, and by the early 1990s P. femorata vanished and M. affinis abundance was low. The decline of M . affinis and P . femorata was followed by an increase of the bivalve Macoma balthica and the arrival of the polychaete genus Marenzelleria in the 1990s. Trend analyses showed a rising trend for near-bottom temperature from the late 1960s. There was a significant decline in salinity during the early 1980s, which stabilized during the early 1990s. A negative trend was observed for oxygen concentration during the entire study period. There were interannual variations in the phosphorus and nitrogen concentrations, and a step-like increase in phosphorus concentration at the beginning of the 21st century. The correlations of temperature, salinity, oxygen, total phosphorus and nitrogen with zoobenthic communities were examined using Constrained Correspondence Analysis. Temperature was the most highly correlated explanatory variable for the benthic species. The study highlights the importance of long-term data sets in assessing the state and ecological processes of zoobenthic systems.

Description: O'Driscoll, R. L., Oeffner, J., and Dunford, A. J. 2013. In situ target strength estimates of optically verified southern blue whiting ( Micromesistius australis ) – ICES Journal of Marine Science, 70: 431–439. Estimates of the acoustic target strength (TS) of southern blue whiting ( Micromesistius australis ) at 38 kHz were obtained using an autonomous acoustic–optical system (AOS) mounted on a demersal trawl. Data were collected from aggregations of spawning adult [mean fork length (FL) 34.4 cm] and immature (mean FL 24.6 cm) southern blue whiting south of New Zealand. Mean TS was estimated from 162 tracks containing 695 echoes from targets identified from video as southern blue whiting. The mean TS was –37.9 dB with a 95% confidence interval ( CI ) of –39.7 to –36.6 dB for 21 immature fish and –34.6 dB (95% CI –35.4 to –34.0 dB) for 141 adults. A logarithmic fit through the mean TS values produced a TS–fork length (FL) relationship from optically verified targets of TS = 22.06 log 10 FL – 68.54. This new relationship gives TS values within 1 dB of those estimated using the relationship recently adopted by ICES for blue whiting ( Micromesistius poutassou ) of TS = 20 log 10 TL – 65.2 (where TL is total length) obtained from in situ measurements, but higher values than those estimated from the previous relationship for southern blue whiting of TS = 38 log 10 FL – 97, which was based on swimbladder modelling.

Description: De Robertis, A., Wilson, C. D., Furnish, S. R., and Dahl, P. H. 2013. Underwater radiated noise measurements of a noise-reduced fisheries research vessel. – ICES Journal of Marine Science, 70: 480–484. Vessel-radiated noise is traditionally measured at naval acoustic ranges, but lower-cost options are desirable for routine monitoring of research vessels. Measurements of a noise-reduced research vessel made at a naval noise range are compared to those made using an experimental mooring equipped with commercially available instrumentation. The measurements from the mooring were precise and within 2.5 dB of those from the noise range at third-octave bands 〈500 Hz. At higher frequencies, direct comparisons were precluded by an intermittent shaft-related noise present only during the mooring measurements, but previously observed at the navy range. The agreement between the two methods suggests that simplified, field-deployable hydrophone systems can be used to accurately characterize the noise signatures of research vessels.

Description: Orphan genes are defined as genes that lack detectable similarity to genes in other species and therefore no clear signals of common descent (i.e., homology) can be inferred. Orphans are an enigmatic portion of the genome because their origin and function are mostly unknown and they typically make up 10% to 30% of all genes in a genome. Several case studies demonstrated that orphans can contribute to lineage-specific adaptation. Here, we study orphan genes by comparing 30 arthropod genomes, focusing in particular on seven recently sequenced ant genomes. This setup allows analyzing a major metazoan taxon and a comparison between social Hymenoptera (ants and bees) and nonsocial Diptera (flies and mosquitoes). First, we find that recently split lineages undergo accelerated genomic reorganization, including the rapid gain of many orphan genes. Second, between the two insect orders Hymenoptera and Diptera, orphan genes are more abundant and emerge more rapidly in Hymenoptera, in particular, in leaf-cutter ants. With respect to intragenomic localization, we find that ant orphan genes show little clustering, which suggests that orphan genes in ants are scattered uniformly over the genome and between nonorphan genes. Finally, our results indicate that the genetic mechanisms creating orphan genes—such as gene duplication, frame-shift fixation, creation of overlapping genes, horizontal gene transfer, and exaptation of transposable elements—act at different rates in insects, primates, and plants. In Formicidae, the majority of orphan genes has their origin in intergenic regions, pointing to a high rate of de novo gene formation or generalized gene loss, and support a recently proposed dynamic model of frequent gene birth and death.

Description: The most bacteria-like mitochondrial genome known is that of the jakobid flagellate Reclinomonas americana NZ. This genome also encodes the largest known gene set among mitochondrial DNAs (mtDNAs), including the RNA subunit of RNase P (transfer RNA processing), a reduced form of transfer–messenger RNA (translational control), and a four-subunit bacteria-like RNA polymerase, which in other eukaryotes is substituted by a nucleus-encoded, single-subunit, phage-like enzyme. Further, protein-coding genes are preceded by potential Shine–Dalgarno translation initiation motifs. Whether similarly ancestral mitochondrial characters also exist in relatives of R. americana NZ is unknown. Here, we report a comparative analysis of nine mtDNAs from five distant jakobid genera: Andalucia, Histiona, Jakoba, Reclinomonas , and Seculamonas . We find that Andalucia godoyi has an even larger mtDNA gene complement than R. americana NZ. The extra genes are rpl35 (a large subunit mitoribosomal protein) and cox15 (involved in cytochrome oxidase assembly), which are nucleus encoded throughout other eukaryotes. Andalucia cox15 is strikingly similar to its homolog in the free-living α-proteobacterium Tistrella mobilis . Similarly, a long, highly conserved gene cluster in jakobid mtDNAs, which is a clear vestige of prokaryotic operons, displays a gene order more closely resembling that in free-living α-proteobacteria than in Rickettsiales species. Although jakobid mtDNAs, overall, are characterized by bacteria-like features, they also display a few remarkably divergent characters, such as 3'-tRNA editing in Seculamonas ecuadoriensis and genome linearization in Jakoba libera . Phylogenetic analysis with mtDNA-encoded proteins strongly supports monophyly of jakobids with Andalucia as the deepest divergence. However, it remains unclear which α-proteobacterial group is the closest mitochondrial relative.

Description: Large-scale evolutionary studies often require the automated construction of alignments of a large number of homologous gene families. The majority of eukaryotic genes can produce different transcripts due to alternative splicing or transcription initiation, and many such transcripts encode different protein isoforms. As analyses tend to be gene centered, one single-protein isoform per gene is selected for the alignment, with the de facto approach being to use the longest protein isoform per gene (Longest), presumably to avoid including partial sequences and to maximize sequence information. Here, we show that this approach is problematic because it increases the number of indels in the alignments due to the inclusion of nonhomologous regions, such as those derived from species-specific exons, increasing the number of misaligned positions. With the aim of ameliorating this problem, we have developed a novel heuristic, Protein ALignment Optimizer (PALO), which, for each gene family, selects the combination of protein isoforms that are most similar in length. We examine several evolutionary parameters inferred from alignments in which the only difference is the method used to select the protein isoform combination: Longest, PALO, the combination that results in the highest sequence conservation, and a randomly selected combination. We observe that Longest tends to overestimate both nonsynonymous and synonymous substitution rates when compared with PALO, which is most likely due to an excess of misaligned positions. The estimation of the fraction of genes that have experienced positive selection by maximum likelihood is very sensitive to the method of isoform selection employed, both when alignments are constructed with MAFFT and with Prank +F . Longest performs better than a random combination but still estimates up to 3 times more positively selected genes than the combination showing the highest conservation, indicating the presence of many false positives. We show that PALO can eliminate the majority of such false positives and thus that it is a more appropriate approach for large-scale analyses than Longest. A web server has been set up to facilitate the use of PALO given a user-defined set of gene families; it is available at http://evolutionarygenomics.imim.es/palo .

Description: Corten, A. 2013. Recruitment depressions in North Sea herring. – ICES Journal of Marine Science, 70:1–15. Two periods of sustained low recruitment have been observed in North Sea herring in recent history; one from 1971 to 1979 and one that started in 2002 and continued up until the most recent year for which information is available (2010). This paper compares both periods of recruitment depression and considers the possible causes for each of them. The first depression (1971–1979) has been commonly ascribed to insufficient egg production by the depleted parent stock. However, recruitment to the central and northern populations was probably also affected by an environmental factor. There are indications of a reduced Atlantic inflow into the northwestern North Sea at the time of the depression. This could have affected survival of the larvae by slowing down their transport to the nursery areas in the southeastern North Sea. For the second period (2002–2010), the low recruitment has been ascribed to an increase in temperature and/or a regime shift in the ecosystem. An alternative explanation could be predation by the large stock of adult herring on its own larvae.

Description: Punt, A. E., Huang, T., and Maunder, M. N. 2013. Review of integrated size-structured models for stock assessment of hard-to-age crustacean and mollusc species. – ICES Journal of Marine Science, 70:16–33. Crustaceans and molluscs such as crabs, rock lobsters, prawns, abalone, and oysters constitute large and valuable fisheries. However, assessments of these species are hampered because they cannot be production aged, in contrast to many teleosts. The major data sources for these species, in addition to catch and abundance index data, are the size compositions of the catches and of any fishery-independent indices. Assessments of such species have been conducted using age-based methods of stock assessment, as well as surplus production models. However, size-structured methods are now preferred because they can make full use of size-composition data, are able to integrate multiple sources of data, and produce the types of outputs which are needed for management purposes. An advantage of size-based models over age-based models is that all processes can be size-based, and these processes can modify the (unmodelled) size-at-age distribution. We review these methods, highlighting the choices that need to be made when developing integrated size-structured stock assessments, the data sources which are typically available and how they are used for parameter estimation, and contrast a number of such assessments worldwide.

Description: De Robertis, A. and Handegard, N. O. 2013. Fish avoidance of research vessels and the efficacy of noise-reduced vessels: a review. – ICES Journal of Marine Science, 70:34–45. It has long been recognized that fish can avoid approaching vessels and that these behaviours can bias fishery surveys. Underwater noise is considered the primary stimulus, and standards for research vessel noise have been established to minimize fish reactions. We review the literature on fish reactions to vessels appearing since these recommendations were made, focusing on acoustic surveys, and compare how fish react to noise-reduced and conventional vessels. Reactions to approaching vessels are variable and difficult to predict. However, the behaviour can bias acoustic abundance measurements, and should be considered when performing acoustic surveys. The few comparisons of acoustic abundance measurements from noise-reduced and conventional vessels are contradictory, but demonstrate that the sound pressure level, on which the noise-reduction criterion is based, is insufficient to explain how fish react to survey vessels. Further research is needed to identify the stimuli fish perceive from approaching vessels and the factors affecting whether fish perceiving these stimuli will react before further recommendations to reduce vessel-avoidance reactions can be made. In the interim, measurement of the biases introduced by fish avoidance reactions during surveys, and timing of surveys when fish are in a less reactive state, may reduce errors introduced by vessel avoidance.

Description: Xavier, J. C., Cherel, Y., Roberts, J., and Piatkowski, U. 2013. How do cephalopods become available to seabirds: can fish gut contents from tuna fishing vessels be a major food source of deep-dwelling cephalopods? – ICES Journal of Marine Science, 70:46–49. Cephalopods are important prey for numerous seabird species. However, the physical mechanisms by which cephalopods (particularly species considered as deep-dwelling) become available to seabirds are poorly understood, and it has recently been suggested that the discarded stomachs of gutted fish captured by tuna longliners can be a major source of deep-dwelling species. Here, we identify some deep-dwelling cephalopods that appear in the diet of seabirds, review the current knowledge of their vertical distribution, and compare the stomach contents of commercially captured tuna with those of seabirds foraging in the same area. The limited available information leads us to conclude that tuna longliners are unlikely to be a major source of deep-dwelling cephalopods for seabirds. However, much more information is required on the ecology of seabird prey, particularly commercially unexploited cephalopod species, which may be obtained from scientific cruises devoted to cephalopod biological research. In addition multispecies/foodweb modelling studies may be required to explore potential interactions between seabirds, their predators and prey, and commercial fishing operations.

Description: Karlson, S., Michalsen, K., and Folkvord, A. 2013. Age determination of Atlantic halibut ( Hippoglossus hippoglossus L.) along the coast of Norway: status and improvements. – ICES Journal of Marine Science, 70:50–55. This study re-evaluates the current ageing methodology for the Atlantic halibut, Hippoglossus hippoglossus . The traditional method is through surface readings of otoliths, but, based on new experiments with different preparation treatments and techniques, a more accurate and cost-efficient procedure for the age determination of Atlantic halibut is proposed.

Description: In vitro studies of the haloarchaeal genus Haloferax have demonstrated their ability to frequently exchange DNA between species, whereas rates of homologous recombination estimated from natural populations in the genus Halorubrum are high enough to maintain random association of alleles between five loci. To quantify the effects of gene transfer and recombination of commonly held (relaxed core) genes during the evolution of the class Halobacteria (haloarchaea), we reconstructed the history of 21 genomes representing all major groups. Using a novel algorithm and a concatenated ribosomal protein phylogeny as a reference, we created a directed horizontal genetic transfer (HGT) network of contemporary and ancestral genomes. Gene order analysis revealed that 90% of testable HGTs were by direct homologous replacement, rather than nonhomologous integration followed by a loss. Network analysis revealed an inverse log-linear relationship between HGT frequency and ribosomal protein evolutionary distance that is maintained across the deepest divergences in Halobacteria. We use this mathematical relationship to estimate the total transfers and amino acid substitutions delivered by HGTs in each genome, providing a measure of chimerism. For the relaxed core genes of each genome, we conservatively estimate that 11–20% of their evolution occurred in other haloarchaea. Our findings are unexpected, because the transfer and homologous recombination of relaxed core genes between members of the class Halobacteria disrupts the coevolution of genes; however, the generation of new combinations of divergent but functionally related genes may lead to adaptive phenotypes not available through cumulative mutations and recombination within a single population.

Description: Waldo, S., and Paulrud, A. 2013. ITQs in Swedish demersal fisheries. – ICES Journal of Marine Science, 70:68–77. Individual transferable quota (ITQ) is a management measure that is widely discussed, not least in the reform of the European Common Fisheries Policy. While the system is expected to reduce overcapacity and improve economic performance, questions are raised concerning the future of small-scale fisheries. This paper uses a model for Swedish fisheries (the Swedish Resource Rent Model for the Commercial Fisheries, SRRMCF) where the economics and fleet structure in a potential Swedish ITQ-system are analysed. The model is an optimization model based on linear programming and data from the European Union's data collection framework. The modelling approach can readily be used by other member states. Three main conclusions can be drawn on how ITQs will affect fisheries: The fishing fleet measured in number of vessels will decrease by approximately 30–50%. Profitability will increase so that the fishing industry will be able to provide competitive wages and make a financial contribution to fisheries management. The system can be designed so that small-scale fisheries are not disadvantaged.

Description: Tzanatos, E., Castro, J., Forcada, A., Matic-Skoko, S., Gaspar, M., and Koutsikopoulos, C. 2013. A Métier-Sustainability-Index (MSI25) to evaluate fisheries components: assessment of cases from data-poor fisheries from southern Europe. – ICES Journal of Marine Science, 70:78–98. The present work presents an integrated comparative assessment of the sustainability of fishing tactics or métiers (combinations of area, season, fishing gear, and target species). An index, which includes biological, conservation, socioeconomic, and management criteria, was created for this purpose: the "Métier Sustainability Index" (MSI25). It was designed using the traffic lights approach and comprises 25 indicators that refer to the "health" of the fisheries. Additionally, the MSI25 provides an indication of the "certainty" of that information with regard to the data type/heterogeneity. Subsequently, the index was tested using data from six case studies (81 métiers) from the Mediterranean and the Atlantic. The index showed that 63.0% of the métiers were "sustainable", 21.0% were "unsustainable" and 16.0% were "under pressure"; 59.3% of characterizations were considered as "certain". The index evaluation showed significant differences between areas and gears. This paper considers the characteristics of the datasets used, the case management status, and the applicability of the index in improving the management of specific fisheries. Overall, the MSI25 was found to be a useful tool for the evaluation and management of highly heterogeneous fisheries and data from various sources and types.