Dividing population genetic distance data with the software Partitioning Optimization with Restricted Growth Strings (PORGS): an application for Chinook salmon (Oncorhynchus tshawytscha), Vancouver Island, British Columbia

A new method of finding the optimal group membership and number of groupings to partition population genetic distance data is presented. The software program Partitioning Optimization with Restricted Growth Strings (PORGS), visits all possible set partitions and deems acceptable partitions to be those that reduce mean intracluster distance. The optimal number of groups is determined with the gap statistic which compares PORGS results with a reference distribution. The PORGS method was validated by a simulated data set with a known distribution. For efficiency, where values of n were larger, restricted growth strings (RGS) were used to bipartition populations during a nested search (bi-PORGS). Bi-PORGS was applied to a set of genetic data from 18 Chinook salmon (Oncorhynchus tshawytscha) populations from the west coast of Vancouver Island. The optimal grouping of these populations corresponded to four geographic locations: 1) Quatsino Sound, 2) Nootka Sound, 3) Clayoquot +Barkley sounds, and 4) southwest Vancouver Island. However, assignment of populations to groups did not strictly reflect the geographical divisions; fish of Barkley Sound origin that had strayed into the Gold River and close genetic similarity between transferred and donor populations meant groupings crossed geographic boundaries. Overall, stock structure determined by this partitioning method was similar to that determined by the unweighted pair-group method with arithmetic averages (UPGMA), an agglomerative clustering algorithm.

The use of otolith morphology to indicate the stock structure of common coral trout (Plectropomus leopardus) on the Great Barrier Reef, Australia

We investigated the use of otolith morphology to indicate the stock structure of an exploited serranid coral reef fish, Plectropomus leopardus, on the Great Barrier Reef (GBR), Australia. Otoliths were measured by traditional one-and two-dimensional measures (otolith length, width, area, perimeter, circularity, and rectangularity), as well as by Fourier analysis to capture the finer details of otolith shape. Variables were compared among four regions of the GBR separated by hundreds of kilometers, as well as among three reefs within each region, hundreds of meters to tens of kilometers apart. The temporal stability in otolith structure was examined by comparing two cohorts of fully recruited four-year-old P. leopardus collected two years before and two years after a signif icant disturbance in the southern parts of the GBR caused by a large tropical cyclone in March 1997. Results indicated the presence of at least two stocks of P. leopardus, although the structure of each stock varied depending on the cohort considered. The results highlight the importance of incorporating data from several years in studies using otolith morphology to discriminate temporary and possibly misleading signals from those that indicate persistent spatial structure in stocks. We conclude that otolith morphology can be used as an initial step to direct further research on groups of P. leopardus that have lived at least a part of their life in different environments.

A Bayesian method for identification of stock mixtures from molecular marker data

Molecular markers have been demonstrated to be useful for the estimation of stock mixture proportions where the origin of individuals is determined from baseline samples. Bayesian statistical methods are widely recognized as providing a preferable strategy for such analyses. In general, Bayesian estimation is based on standard latent class models using data augmentation through Markov chain Monte Carlo techniques. In this study, we introduce a novel approach based on recent developments in the estimation of genetic population structure. Our strategy combines analytical integration with stochastic optimization to identify stock mixtures. An important enhancement over previous methods is the possibility of appropriately handling data where only partial baseline sample information is available. We address the potential use of nonmolecular, auxiliary biological information in our Bayesian model.

Inferring Population Genetic Structure in Widely and Continuously Distributed Carnivores: The Stone Marten (Martes foina) as a Case Study

The stone marten is a widely distributed mustelid in the Palaearctic region that exhibits variable habitat preferences in different parts of its range. The species is a Holocene immigrant from southwest Asia which, according to fossil remains, followed the expansion of the Neolithic farming cultures into Europe and possibly colonized the Iberian Peninsula during the Early Neolithic (ca. 7,000 years BP). However, the population genetic structure and historical biogeography of this generalist carnivore remains essentially unknown. In this study we have combined mitochondrial DNA (mtDNA) sequencing (621 bp) and microsatellite genotyping (23 polymorphic markers) to infer the population genetic structure of the stone marten within the Iberian Peninsula. The mtDNA data revealed low haplotype and nucleotide diversities and a lack of phylogeographic structure, most likely due to a recent colonization of the Iberian Peninsula by a few mtDNA lineages during the Early Neolithic. The microsatellite data set was analysed with a) spatial and non-spatial Bayesian individual-based clustering (IBC) approaches (STRUCTURE, TESS, BAPS and GENELAND), and b) multivariate methods [discriminant analysis of principal components (DAPC) and spatial principal component analysis (sPCA)]. Additionally, because isolation by distance (IBD) is a common spatial genetic pattern in mobile and continuously distributed species and it may represent a challenge to the performance of the above methods, the microsatellite data set was tested for its presence. Overall, the genetic structure of the stone marten in the Iberian Peninsula was characterized by a NE-SW spatial pattern of IBD, and this may explain the observed disagreement between clustering solutions obtained by the different IBC methods. However, there was significant indication for contemporary genetic structuring, albeit weak, into at least three different subpopulations. The detected subdivision could be attributed to the influence of the rivers Ebro, Tagus and Guadiana, suggesting that main watercourses in the Iberian Peninsula may act as semi-permeable barriers to gene flow in stone martens. To our knowledge, this is the first phylogeographic and population genetic study of the species at a broad regional scale. We also wanted to make the case for the importance and benefits of using and comparing multiple different clustering and multivariate methods in spatial genetic analyses of mobile and continuously distributed species.

New Nuclear SNP Markers Unravel the Genetic Structure and Effective Population Size of Albacore Tuna (Thunnus alalunga)

In the present study we have investigated the population genetic structure of albacore (Thunnus alalunga, Bonnaterre 1788) and assessed the loss of genetic diversity, likely due to overfishing, of albacore population in the North Atlantic Ocean. For this purpose, 1,331 individuals from 26 worldwide locations were analyzed by genotyping 75 novel nuclear SNPs. Our results indicated the existence of four genetically homogeneous populations delimited within the Mediterranean Sea, the Atlantic Ocean, the Indian Ocean and the Pacific Ocean. Current definition of stocks allows the sustainable management of albacore since no stock includes more than one genetic entity. In addition, short-and long-term effective population sizes were estimated for the North Atlantic Ocean albacore population, and results showed no historical decline for this population. Therefore, the genetic diversity and, consequently, the adaptive potential of this population have not been significantly affected by overfishing.

A comparison of two fishery-independent survey programs used to define the population structure of American lobster (Homarus americanus) in the Gulf of Maine

The population structure and abundance of the American lobster (Homarus americanus) stock in the Gulf of Maine are defined by data derived from a fishery-independent trawl survey program conducted by the National Marine Fisheries Service (NMFS). Few sampling stations in the survey area are located inshore, in particular along coastal Maine. According to statistics, however, more than two thirds of the lobster landings come from inshore waters within three miles off the coast of Maine. In order to include an inshore survey program, complementary to the NMFS survey, the Maine Department of Marine Resources (DMR) initialized an inshore survey program in 2000. The survey was modeled on the NMFS survey program, making these two survey programs comparable. Using data from both survey programs, we evaluated the population structure of the American lobster in the Gulf of Maine. Our findings indicate that lobsters in the Gulf of Maine tend to have a size-dependent inshore-off-shore distribution; smaller lobsters are more likely to stay inshore and larger lobsters are more likely to stay offshore. The DMR inshore and NMFS survey programs focused on different areas in the Gulf of Maine and likely targeted different segments of the stock. We suggest that data from both survey programs be used to assess the lobster stock and to describe the dynamics of the stock in the Gulf of Maine.

Population structure and variance effective size of red snapper (Lutjanus campechanus) in the northern Gulf of Mexico*

We assayed allelic variation at 19 nuclear-encoded microsatellites among 1622 Gulf red snapper (Lutjanus campechanus) sampled from the 1995 and 1997 cohorts at each of three offshore localities in the northern Gulf of Mexico (Gulf). Localities represented western, central, and eastern subregions within the northern Gulf. Number of alleles per microsatellite per sample ranged from four to 23, and gene diversity ranged from 0.170 to 0.917. Tests of conformity to Hardy-Weinberg equilibrium expectations and of genotypic equilibrium between pairs of micro-satellites were generally nonsignificant following Bonferroni correction. Significant genic or genotypic heterogeneity (or both) among samples was detected at four microsatellites and over all microsatellites. Levels of divergence among samples were low (FST ≤0.001). Pairwise exact tests revealed that six of seven “significant” comparisons involved temporal rather than spatial heterogeneity. Contemporaneous or variance effective size (NeV) was estimated from the temporal variance in allele frequencies by using a maximum-likelihood method. Estimates of NeV ranged between 1098 and >75,000 and differed significantly among localities; the NeV estimate for the sample from the northcentral Gulf was >60 times as large as the estimates for the other two localities. The differences in variance effective size could ref lect differences in number of individuals successfully reproducing, differences in patterns and intensity of immigration, or both, and are consistent with the hypothesis, supported by life-history data, that different “demographic stocks” of red snapper are found in the northern Gulf. Estimates of NeV for red snapper in the northern Gulf were at least three orders of magnitude lower than current estimates of census size (N). The ratio of effective to census size (Ne/N) is far below that expected in an ideal population and may reflect high variance in individual reproductive success, high temporal and spatial variance in productivity among subregions or a combination of the two.

Structure-function analysis of USP1: insights into the role of Ser313 phosphorylation site and the effect of cancer-associated mutations on autocleavage

Background: In complex with its cofactor UAF1, the USP1 deubiquitinase plays an important role in cellular processes related to cancer, including the response to DNA damage. The USP1/UAF1 complex is emerging as a novel target in cancer therapy, but several aspects of its function and regulation remain to be further clarified. These include the role of the serine 313 phosphorylation site, the relative contribution of different USP1 sequence motifs to UAF1 binding, and the potential effect of cancer-associated mutations on USP1 regulation by autocleavage. Methods: We have generated a large set of USP1 structural variants, including a catalytically inactive form (C90S), non-phosphorylatable (S313A) and phosphomimetic (S313D) mutants, deletion mutants lacking potential UAF1 binding sites, a mutant (GG/AA) unable to undergo autocleavage at the well-characterized G670/G671 diglycine motif, and four USP1 mutants identified in tumor samples that cluster around this cleavage site (G667A, L669P, K673T and A676T). Using cell-based assays, we have determined the ability of these mutants to bind UAF1, to reverse DNA damage-induced monoubiquitination of PCNA, and to undergo autocleavage. Results: A non-phosphorylatable S313A mutant of USP1 retained the ability to bind UAF1 and to reverse PCNA ubiquitination in cell-based assays. Regardless of the presence of a phosphomimetic S313D mutation, deletion of USP1 fragment 420-520 disrupted UAF1 binding, as determined using a nuclear relocation assay. The UAF1 binding site in a second UAF1-interacting DUB, USP46, was mapped to a region homologous to USP1(420-520). Regarding USP1 autocleavage, co-expression of the C90S and GG/AA mutants did not result in cleavage, while the cancer-associated mutation L669P was found to reduce cleavage efficiency. Conclusions: USP1 phosphorylation at S313 is not critical for PCNA deubiquitination, neither for binding to UAF1 in a cellular environment. In this context, USP1 amino acid motif 420-520 is necessary and sufficient for UAF1 binding. This motif, and a homologous amino acid segment that mediates USP46 binding to UAF1, map to the Fingers sub-domain of these DUBs. On the other hand, our results support the view that USP1 autocleavage may occur in cis, and can be altered by a cancer-associated mutation.

Energy analysis of the ring seine operations, off Cochin, Kerala

Ring seines are lightly constructed purse seines adapted for operation in the traditional sector. Fish production and energy requirement in the ring seine operations, off Cochin, Kerala, India are discussed in this paper, based on data collected during 1997- 1998. The results reflect the Gross Energy Requirement (GER) situation that existed during 1997-1998. Mean catch per ring seiner per year worked out to be 211.9 t of which sardines (Sardinella spp.) constituted 44.3%, followed by Indian mackerel (Rastrelliger kanagurta) 29.7%, carangids 11.4%, penaeid prawns 2.2%, pomfrets 1.1% and miscellaneous fish 11.3%. Total energy inputs into the ring seine operations were estimated to be 1300.8 GJ. Output by way of fish production was determined to be 931.85 GJ. GER is the sum of all non-renewable energy resources consumed in making available a product or service and is a measure of intensity of non-renewable resource use. GER per tonne of fish landed by ring seiners was estimated to be 6.14. Among the operational inputs, kerosene constituted 73.4% of the GER, followed by petrol (12.7%), diesel (6.7%) and lubricating oil (2.4%). Fishing gear contributed 3.8%, engine 0.8% and fishing craft 0.3% of the GER. Energy ratio for ring seining was 0.72 and energy intensity 1.40.

Geomática aplicada às medidas preventivas e corretivas para combater danos ambientais: o caso dos incidentes rodoviários com produtos perigosos.

Este trabalho aborda a aplicação da Geomática no apoio às ações de rápida resposta às emergências ambientais passíveis de ocorrência, nas rodovias do País, por incidentes com produtos perigosos. A cadeia logística nacional de infraestrutura de transporte de carga está apoiada, principalmente, no modal rodoviário, que é responsável por cerca de 60% do volume total transportado. Seguindo essa tendência, a maioria dos produtos perigosos também é transportada através de estradas. Nessas operações, apesar de todas as medidas de segurança empregadas, tanto na embalagem quanto no manuseio, há sempre a possibilidade da ocorrência de incidentes com esse tipo de carga, no deslocamento entre o ponto de partida e o destino final. Diante dessas potenciais ameaças, ficam evidentes os riscos à integridade física e patrimonial das populações lindeiras, bem como ao meio ambiente, no tocante à qualidade das águas superficiais e subterrâneas, do solo e do ar. Este trabalho sugere, com base na integração dos recursos tecnológicos do Sensoriamento Remoto, do Sistema de Posicionamento Global-GPS e do Geoprocessamento, implantar um sistema de informações georreferenciadas que proporcione os elementos necessários às ações de rápida resposta a esses eventos. A estrutura proposta, denominada Bureau de Informações Territoriais, está moldada para proporcionar uma solução para os problemas relacionados com o posicionamento na superfície terrestre, e que engloba os conceitos da Geomática, através da aplicação das suas quatro etapas básicas: coleta de dados, análise de dados, distribuição da informação e uso da informação. Assim, atendendo ao princípio da mobilidade, o projeto visa, ainda, a disponibilizar na tela de um laptop as informações necessárias ao planejamento e à execução das ações, no menor tempo possível. Por outro lado, considerando que o Bureau também proporciona às organizações que participam das ações de combate aos incidentes rodoviários com produtos perigosos as informações necessárias, através de uma conexão à Internet, espera-se a maximização dos benefícios do planejamento das ações preventivas e corretivas de resposta aos sinistros, assim como a minimização dos efeitos desse tipo de evento sobre a população, as propriedades e o meio ambiente.

Calculation of the electronic structure of carbon films using electron energy loss spectroscopy.

The detailed understanding of the electronic properties of carbon-based materials requires the determination of their electronic structure and more precisely the calculation of their joint density of states (JDOS) and dielectric constant. Low electron energy loss spectroscopy (EELS) provides a continuous spectrum which represents all the excitations of the electrons within the material with energies ranging between zero and about 100 eV. Therefore, EELS is potentially more powerful than conventional optical spectroscopy which has an intrinsic upper information limit of about 6 eV due to absorption of light from the optical components of the system or the ambient. However, when analysing EELS data, the extraction of the single scattered data needed for Kramers Kronig calculations is subject to the deconvolution of the zero loss peak from the raw data. This procedure is particularly critical when attempting to study the near-bandgap region of materials with a bandgap below 1.5 eV. In this paper, we have calculated the electronic properties of three widely studied carbon materials; namely amorphous carbon (a-C), tetrahedral amorphous carbon (ta-C) and C60 fullerite crystal. The JDOS curve starts from zero for energy values below the bandgap and then starts to rise with a rate depending on whether the material has a direct or an indirect bandgap. Extrapolating a fit to the data immediately above the bandgap in the stronger energy loss region was used to get an accurate value for the bandgap energy and to determine whether the bandgap is direct or indirect in character. Particular problems relating to the extraction of the single scattered data for these materials are also addressed. The ta-C and C60 fullerite materials are found to be direct bandgap-like semiconductors having a bandgaps of 2.63 and 1.59eV, respectively. On the other hand, the electronic structure of a-C was unobtainable because it had such a small bandgap that most of the information is contained in the first 1.2 eV of the spectrum, which is a region removed during the zero loss deconvolution.

Structure learning in a sensorimotor association task.

Learning is often understood as an organism's gradual acquisition of the association between a given sensory stimulus and the correct motor response. Mathematically, this corresponds to regressing a mapping between the set of observations and the set of actions. Recently, however, it has been shown both in cognitive and motor neuroscience that humans are not only able to learn particular stimulus-response mappings, but are also able to extract abstract structural invariants that facilitate generalization to novel tasks. Here we show how such structure learning can enhance facilitation in a sensorimotor association task performed by human subjects. Using regression and reinforcement learning models we show that the observed facilitation cannot be explained by these basic models of learning stimulus-response associations. We show, however, that the observed data can be explained by a hierarchical Bayesian model that performs structure learning. In line with previous results from cognitive tasks, this suggests that hierarchical Bayesian inference might provide a common framework to explain both the learning of specific stimulus-response associations and the learning of abstract structures that are shared by different task environments.

Discovering transcriptional modules by Bayesian data integration.

MOTIVATION: We present a method for directly inferring transcriptional modules (TMs) by integrating gene expression and transcription factor binding (ChIP-chip) data. Our model extends a hierarchical Dirichlet process mixture model to allow data fusion on a gene-by-gene basis. This encodes the intuition that co-expression and co-regulation are not necessarily equivalent and hence we do not expect all genes to group similarly in both datasets. In particular, it allows us to identify the subset of genes that share the same structure of transcriptional modules in both datasets. RESULTS: We find that by working on a gene-by-gene basis, our model is able to extract clusters with greater functional coherence than existing methods. By combining gene expression and transcription factor binding (ChIP-chip) data in this way, we are better able to determine the groups of genes that are most likely to represent underlying TMs. AVAILABILITY: If interested in the code for the work presented in this article, please contact the authors. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Integration of fisheries acoustics surveys and bathymetric mapping to characterize midwater-seafloor habitats of US Virgin Islands and Puerto Rico (2008-2010)

NOAA’s Coral Reef Conservation program (CRCP) develops coral reef management priorities by bringing together various partners to better understand threats to coral reef ecosystems with the goal of conserving, protecting and restoring these resources. Place-based and ecosystem-based management approaches employed by CRCP require that spatially explicit information about benthic habitats and fish utilization are available to characterize coral reef ecosystems and set conservation priorities. To accomplish this, seafloor habitat mapping of coral reefs around the U.S. Virgin Islands (USVI) and Puerto Rico has been ongoing since 2004. In 2008, fishery acoustics surveys were added to NOAA survey missions in the USVI and Puerto Rico to assess fish distribution and abundance in relation to benthic habitats in high priority conservation areas. NOAA’s National Centers for Coastal Ocean Science (NCCOS) have developed fisheries acoustics survey capabilities onboard the NOAA ship Nancy Foster to complement the CRCP seafloor habitat mapping effort spearheaded by the Center for Coastal Monitoring and Assessment Biogeography Branch (CCMA-BB). The integration of these activities has evolved on the Nancy Foster over the three years summarized in this report. A strategy for improved operations and products has emerged over that time. Not only has the concurrent operation of multibeam and fisheries acoustics surveys been beneficial in terms of optimizing ship time and resources, this joint effort has advanced an integrated approach to characterizing bottom and mid-water habitats and the fishes associated with them. CCMA conducts multibeam surveys to systematically map and characterize coral reef ecosystems, resulting in products such as high resolution bathymetric maps, backscatter information, and benthic habitat classification maps. These products focus on benthic features and live bottom habitats associated with them. NCCOS Centers (the Center for Coastal Fisheries and Habitat Research and the Center for Coastal Environmental Health and Biomolecular Research) characterize coral reef ecosystems by using fisheries acoustics methods to capture biological information through the entire water column. Spatially-explicit information on marine resources derived from fisheries acoustics surveys, such as maps of fish density, supports marine spatial planning strategies and decision making by providing a biological metric for evaluating coral reef ecosystems and assessing impacts from pollution, fishing pressure, and climate change. Data from fisheries acoustics surveys address management needs by providing a measure of biomass in management areas, detecting spatial and temporal responses in distribution relative to natural and anthropogenic impacts, and identifying hotspots that support high fish abundance or fish aggregations. Fisheries acoustics surveys conducted alongside multibeam mapping efforts inherently couple water column data with information on benthic habitats and provide information on the heterogeneity of both benthic habitats and biota in the water column. Building on this information serves to inform resource managers regarding how fishes are organized around habitat structure and the scale at which these relationships are important. Where resource managers require place-based assessments regarding the location of critical habitats along with high abundances of fish, concurrent multibeam and fisheries acoustics surveys serve as an important tool for characterizing and prioritizing coral reef ecosystems. This report summarizes the evolution of fisheries acoustics surveys onboard the NOAA ship Nancy Foster from 2008 to 2010, in conjunction with multibeam data collection, aimed at characterizing benthic and mid-water habitats in high priority conservation areas around the USVI and Puerto Rico. It also serves as a resource for the continued development of consistent data products derived from acoustic surveys. By focusing on the activities of 2010, this report highlights the progress made to date and illustrates the potential application of fisheries data derived from acoustic surveys to the management of coral reef ecosystems.

Red snapper (Lutjanus campechanus) demographic structure in the northern Gulf of Mexico based on spatial patterns in growth rates and morphometrics

Red snapper (Lutjanus campechanus) in the United States waters of the Gulf of Mexico (GOM) has been considered a single unit stock since management of the species began in 1991. The validity of this assumption is essential to management decisions because measures of growth can differ for nonmixing populations. We examined growth rates, size-at-age, and length and weight information of red snapper collected from the recreational harvests of Alabama (n=2010), Louisiana (n=1905), and Texas (n =1277) from 1999 to 2001. Ages were obtained from 5035 otolith sections and ranged from one to 45 years. Fork length, total weight, and age-frequency distributions differed significantly among all states; Texas, however, had a much higher proportion of smaller, younger fish. All red snapper showed rapid growth until about age 10 years, after which growth slowed considerably. Von Bertalanffy growth models of both mean fork length and mean total weight-at-age predicted significantly smaller fish at age from Texas, whereas no differences were found between Alabama and Louisiana models. Texas red snapper were also shown to differ significantly from both Alabama and Louisiana red snapper in regressions of mean weight at age. Demographic variation in growth rates may indicate the existence of separate management units of red snapper in the GOM. Our data indicate that the red snapper inhabiting the waters off Texas are reaching smaller maximum sizes at a faster rate and have a consistently smaller total weight at age than those collected from Louisiana and Alabama waters. Whether these differences are environmentally induced or are the result of genetic divergence remains to be determined, but they should be considered for future management regulations.