Combining genotypic and phenotypic variation in a geospatial framework to identify sources of mussels in northern New Zealand

The New Zealand green-lipped mussel aquaculture industry is largely dependent on the supply of young mussels that wash up on Ninety Mile Beach (so-called Kaitaia spat), which are collected and trucked to aquaculture farms. The locations of source populations of Kaitaia spat are unknown and this lack of knowledge represents a major problem because spat supply may be irregular. We combined genotypic (microsatellite) and phenotypic (shell geochemistry) data in a geospatial framework to determine if this new approach can help identify source populations of mussels collected from two spat-collecting and four non-spat-collecting sites further south. Genetic analyses resolved differentiated clusters (mostly three clusters), but no obvious source populations. Shell geochemistry analyses resolved six differentiated clusters, as did the combined genotypic and phenotypic data. Analyses revealed high levels of spatial and temporal variability in the geochemistry signal. Whilst we have not been able to identify the source site(s) of Kaitaia spat our analyses indicate that geospatial testing using combined genotypic and phenotypic data is a powerful approach. Next steps should employ analyses of single nucleotide polymorphism markers with shell geochemistry and in conjunction with high resolution physical oceanographic modelling to resolve the longstanding question of the origin of Kaitaia spat

Population structure, genetic connectivity, and signatures of local adaptation of the giant black tiger shrimp (Penaeus monodon) throughout the indo-pacific region

The giant black tiger shrimp (Penaeus monodon) is native to the Indo-Pacific and is the second most farmed penaeid shrimp species globally. Understanding genetic structure, connectivity, and local adaptation among Indo-Pacific black tiger shrimp populations is important for informing sustainable fisheries management and aquaculture breeding programs. Population genetic and outlier detection analyses were undertaken using 10,593 genome-wide single nucleotide polymorphisms (SNPs) from 16 geographically disparate Indo-Pacific P. monodon populations. Levels of genetic diversity were highest for Southeast Asian populations and were lowest for Western Indian Ocean (WIO) populations. Both neutral (n = 9,930) and outlier (n = 663) loci datasets revealed a pattern of strong genetic structure of P. monodon corresponding with broad geographical regions and clear genetic breaks among samples within regions. Neutral loci revealed seven genetic clusters and the separation of Fiji and WIO clusters from all other clusters, whereas outlier loci revealed six genetic clusters and high genetic differentiation among populations. The neutral loci dataset estimated five migration events that indicated migration to Southeast Asia from the WIO, with partial connectivity to populations in both oceans. We also identified 26 putatively adaptive SNPs that exhibited significant Pearson correlation (P < 0.05) between minor allele frequency and maximum or minimum sea surface temperature. Matched transcriptome contig annotations suggest putatively adaptive SNPs involvement in cellular and metabolic processes, pigmentation, immune response, and currently unknown functions. This study provides novel genome-level insights that have direct implications for P. monodon aquaculture and fishery management practices

Fine-scale population structure and evidence for local adaptation in Australian giant black tiger shrimp (Penaeus monodon) using SNP analysis

Background: Restrictions to gene flow, genetic drift, and divergent selection associated with different environments are significant drivers of genetic differentiation. The black tiger shrimp (Penaeus monodon), is widely distributed throughout the Indian and Pacific Oceans including along the western, northern and eastern coastline of Australia, where it is an important aquaculture and fishery species. Understanding the genetic structure and the influence of environmental factors leading to adaptive differences among populations of this species is important for farm genetic improvement programs and sustainable fisheries management. Results: Based on 278 individuals obtained from seven geographically disparate Australian locations, 10,624 high-quality SNP loci were used to characterize genetic diversity, population structure, genetic connectivity, and adaptive divergence. Significant population structure and differentiation were revealed among wild populations (average FST = 0.001–0.107; p <  0.05). Eighty-nine putatively outlier SNPs were identified to be potentially associated with environmental variables by using both population differentiation (BayeScan and PCAdapt) and environmental association (redundancy analysis and latent factor mixed model) analysis methods. Clear population structure with similar spatial patterns were observed in both neutral and outlier markers with three genetically distinct groups identified (north Queensland, Northern Territory, and Western Australia). Redundancy, partial redundancy, and multiple regression on distance matrices analyses revealed that both geographical distance and environmental factors interact to generate the structure observed across Australian P. monodon populations. Conclusion: This study provides new insights on genetic population structure of Australian P. monodon in the face of environmental changes, which can be used to advance sustainable fisheries management and aquaculture breeding programs

Search for lepton-flavor violation at HERA

A search for lepton-flavor-violating interactions $e p \to \mu X$ and $e p\to \tau X$ has been performed with the ZEUS detector using the entire HERA I data sample, corresponding to an integrated luminosity of 130 pb^{-1}. The data were taken at center-of-mass energies, $\sqrt{s}$, of 300 and 318 GeV. No evidence of lepton-flavor violation was found, and constraints were derived on leptoquarks (LQs) that could mediate such interactions. For LQ masses below $\sqrt{s}$, limits were set on $\lambda_{eq_1} \sqrt{\beta_{\ell q}}$, where $\lambda_{eq_1}$ is the coupling of the LQ to an electron and a first-generation quark $q_1$, and $\beta_{\ell q}$ is the branching ratio of the LQ to the final-state lepton $\ell$ ($\mu$ or $\tau$) and a quark $q$. For LQ masses much larger than $\sqrt{s}$, limits were set on the four-fermion interaction term $\lambda_{e q_\alpha} \lambda_{\ell q_\beta} / M_{\mathrm{LQ}}^2$ for LQs that couple to an electron and a quark $q_\alpha$ and to a lepton $\ell$ and a quark $q_\beta$, where $\alpha$ and $\beta$ are quark generation indices. Some of the limits are also applicable to lepton-flavor-violating processes mediated by squarks in $R$-Parity-violating supersymmetric models. In some cases, especially when a higher-generation quark is involved and for the process $e p\to \tau X$, the ZEUS limits are the most stringent to date.Comment: 37 pages, 10 figures, Accepted by EPJC. References and 1 figure (Fig. 6) adde

Multijet production in neutral current deep inelastic scattering at HERA and determination of alpha_s

Multijet production rates in neutral current deep inelastic scattering have been measured in the range of exchanged boson virtualities 10 < Q2 < 5000 GeV2. The data were taken at the ep collider HERA with centre-of-mass energy sqrt(s) = 318 GeV using the ZEUS detector and correspond to an integrated luminosity of 82.2 pb-1. Jets were identified in the Breit frame using the k_T cluster algorithm in the longitudinally invariant inclusive mode. Measurements of differential dijet and trijet cross sections are presented as functions of jet transverse energy E_{T,B}{jet}, pseudorapidity eta_{LAB}{jet} and Q2 with E_{T,B}{jet} > 5 GeV and -1 < eta_{LAB}{jet} < 2.5. Next-to-leading-order QCD calculations describe the data well. The value of the strong coupling constant alpha_s(M_Z), determined from the ratio of the trijet to dijet cross sections, is alpha_s(M_Z) = 0.1179 pm 0.0013(stat.) {+0.0028}_{-0.0046}(exp.) {+0.0064}_{-0.0046}(th.)Comment: 22 pages, 5 figure

Modeling the Subsurface Structure of Sunspots

While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat.Comment: 73 pages, 19 figures, accepted by Solar Physic

Measurement of (anti)deuteron and (anti)proton production in DIS at HERA

The first observation of (anti)deuterons in deep inelastic scattering at HERA has been made with the ZEUS detector at a centre-of-mass energy of 300--318 GeV using an integrated luminosity of 120 pb-1. The measurement was performed in the central rapidity region for transverse momentum per unit of mass in the range 0.3<p_T/M<0.7. The particle rates have been extracted and interpreted in terms of the coalescence model. The (anti)deuteron production yield is smaller than the (anti)proton yield by approximately three orders of magnitude, consistent with the world measurements.Comment: 26 pages, 9 figures, 5 tables, submitted to Nucl. Phys.

Photoproduction of D∗±D^{*\pm} mesons associated with a leading neutron

The photoproduction of $D^{*\pm} (2010)$ mesons associated with a leading neutron has been observed with the ZEUS detector in $ep$ collisions at HERA using an integrated luminosity of 80 pb$^{-1}$. The neutron carries a large fraction, {$x_L>0.2$}, of the incoming proton beam energy and is detected at very small production angles, {$\theta_n<0.8$ mrad}, an indication of peripheral scattering. The $D^*$ meson is centrally produced with pseudorapidity {$|\eta| 1.9$ GeV}, which is large compared to the average transverse momentum of the neutron of 0.22 GeV. The ratio of neutron-tagged to inclusive $D^*$ production is $8.85\pm 0.93({\rm stat.})^{+0.48}_{-0.61}({\rm syst.})\%$ in the photon-proton center-of-mass energy range {$130 <W<280$ GeV}. The data suggest that the presence of a hard scale enhances the fraction of events with a leading neutron in the final state.Comment: 28 pages, 4 figures, 2 table

Measurement of the open-charm contribution to the diffractive proton structure function

Production of D*+/-(2010) mesons in diffractive deep inelastic scattering has been measured with the ZEUS detector at HERA using an integrated luminosity of 82 pb^{-1}. Diffractive events were identified by the presence of a large rapidity gap in the final state. Differential cross sections have been measured in the kinematic region 1.5 < Q^2 < 200 GeV^2, 0.02 < y < 0.7, x_{IP} < 0.035, beta 1.5 GeV and |\eta(D*+/-)| < 1.5. The measured cross sections are compared to theoretical predictions. The results are presented in terms of the open-charm contribution to the diffractive proton structure function. The data demonstrate a strong sensitivity to the diffractive parton densities.Comment: 35 pages, 11 figures, 6 table

Measurement of beauty production in deep inelastic scattering at HERA

The beauty production cross section for deep inelastic scattering events with at least one hard jet in the Breit frame together with a muon has been measured, for photon virtualities Q^2 > 2 GeV^2, with the ZEUS detector at HERA using integrated luminosity of 72 pb^-1. The total visible cross section is sigma_b-bbar (ep -> e jet mu X) = 40.9 +- 5.7 (stat.) +6.0 -4.4 (syst.) pb. The next-to-leading order QCD prediction lies about 2.5 standard deviations below the data. The differential cross sections are in general consistent with the NLO QCD predictions; however at low values of Q^2, Bjorken x, and muon transverse momentum, and high values of jet transverse energy and muon pseudorapidity, the prediction is about two standard deviations below the data.Comment: 18 pages, 4 figure