ALBERT

Beschreibung: Pre-main sequence (PMS) stars evolve into main sequence (MS) phase over a period of time. Interestingly, we found a scarcity of studies in existing literature that examines and attempts to better understand the stars in PMS to MS transition phase. The purpose of the present study is to detect such rare stars, which we named as ‘Transition Phase’ (TP) candidates – stars evolving from the PMS to the MS phase. We identified 98 TP candidates using photometric analysis of a sample of 2167 classical Be (CBe) and 225 Herbig Ae/Be (HAeBe) stars. This identification is done by analyzing the near- and mid-infrared excess and their location in the optical colour-magnitude diagram. The age and mass of 58 of these TP candidates are determined to be between 0.1–5 Myr and 2–10.5 M⊙, respectively. The TP candidates are found to possess rotational velocity and colour excess values in between CBe and HAeBe stars, which is reconfirmed by generating a set of synthetic samples using the machine learning approach.

Beschreibung: Recent searches for the hosts of z ∼ 4 damped Lyα absorbers (DLAs) have detected bright galaxies at distances of tens of kpc from the DLA. Using the FIRE-2 cosmological zoom simulations, we argue that these relatively large distances are due to a predominantly cool and neutral inner circumgalactic medium (CGM) surrounding high-redshift galaxies. The inner CGM is cool because of the short cooling time of hot gas in ≲ 1012 M⊙ haloes, which implies that accretion and feedback energy are radiated quickly, while it is neutral due to high volume densities and column densities at high redshift which shield cool gas from photoionization. Our analysis predicts large DLA covering factors ($gtrsim 50{{ m per cent}}$) out to impact parameters ∼0.3((1 + z)/5)3/2Rvir from the central galaxies at z ≳ 1, equivalent to a proper distance of $sim 21, M_{12}^{1/3} left(left(1+z ight)/5 ight)^{1/2}, { m kpc}$ (Rvir and M12 are the halo virial radius and mass in units of 1012 M⊙, respectively). This implies that DLA covering factors at z ∼ 4 may be comparable to unity out to a distance ∼10 times larger than stellar half-mass radii. A predominantly neutral inner CGM in the early universe suggests that its mass and metallicity can be directly constrained by absorption surveys, without resorting to the large ionization corrections as required for ionized CGM.

Beschreibung: Motivation Accurate automatic annotation of protein function relies on both innovative models and robust data sets. Due to their importance in biological processes, the identification of DNA-binding proteins directly from protein sequence has been the focus of many studies. However, the data sets used to train and evaluate these methods have suffered from substantial flaws. We describe some of the weaknesses of the data sets used in previous DNA-binding protein literature and provide several new data sets addressing these problems. We suggest new evaluative benchmark tasks that more realistically assess real-world performance for protein annotation models. We propose a simple new model for the prediction of DNA-binding proteins and compare its performance on the improved data sets to two previously published models. Additionally, we provide extensive tests showing how the best models predict across taxonomies. Results Our new gradient boosting model, which uses features derived from a published protein language model, outperforms the earlier models. Perhaps surprisingly, so does a baseline nearest neighbor model using BLAST percent identity. We evaluate the sensitivity of these models to perturbations of DNA-binding regions and control regions of protein sequences. The successful data-driven models learn to focus on DNA-binding regions. When predicting across taxonomies, the best models are highly accurate across species in the same kingdom and can provide some information when predicting across kingdoms. Code and Data Availability The data and results for this paper can be found at https://doi.org/10.5281/zenodo.5153906. The code for this paper can be found at https://doi.org/10.5281/zenodo.5153683. The code, data and results can also be found at https://github.com/AZaitzeff/tools_for_dna_binding_proteins.

Beschreibung: The Internal Linear Combination (ILC) method has been extensively used to extract the cosmic microwave background (CMB) anisotropy map from foreground contaminated multi-frequency maps. However, the performance of simple ILC is limited and can be significantly improved by heavily constraint equations, dubbed constrained ILC (cILC). The standard ILC and cILC works on spin-0 fields. Recently, a generalised version of ILC has been developed, named Polarization ILC (PILC), in which Q ± iU at multiple frequencies are combined using complex coefficients to estimate Stokes Q and U maps. A statistical moment expansion method has recently been developed for high precision modelling of the Galactic foregrounds. This paper develops a semi-blind component separation method combining the moment approach of foreground modelling with a generalised version of the PILC method for heavily constraint equations. The algorithm is developed in pixel space over a spin-2 field. We demonstrate the performance of the method on three sets of absolutely calibrated simulated maps at WMAP and Planck frequencies with varying foreground models. We apply this component separation technique in simultaneous estimation of Stokes Q and U maps of the thermal dust at 353 GHz and synchrotron at 30 GHz. We also recover both dust and synchrotron maps at 100 and 143 GHz, where separating two components is challenging.

Beschreibung: Luminous active galactic nuclei and X-ray binaries often contain geometrically thin, radiatively cooled accretion discs. According to theory, these are – in many cases – initially highly misaligned with the black hole equator. In this work, we present the first general relativistic magnetohydrodynamic simulations of very thin (h/r ∼ 0.015–0.05) accretion discs around rapidly spinning (a ∼ 0.9) black holes and tilted by 45°–65°. We show that the inner regions of the discs with h/r ≲ 0.03 align with the black hole equator, though out to smaller radii than predicted by analytic work. The inner aligned and outer misaligned disc regions are separated by a sharp break in tilt angle accompanied by a sharp drop in density. We find that frame dragging by the spinning black hole overpowers the disc viscosity, which is self-consistently produced by magnetized turbulence, tearing the disc apart and forming a rapidly precessing inner sub-disc surrounded by a slowly precessing outer sub-disc. We find that the system produces a pair of relativistic jets for all initial tilt values. At small distances, the black hole launched jets precess rapidly together with the inner sub-disc, whereas at large distances they partially align with the outer sub-disc and precess more slowly. If the tearing radius can be modeled accurately in future work, emission model independent measurements of black hole spin based on precession-driven quasi-periodic oscillations may become possible.

Beschreibung: Atlantic salmon (Salmo salar) populations throughout the North Atlantic have declined in recent decades largely due to reduced marine survival, yet our understanding of marine distribution patterns and migratory routes remains limited. Here, we assigned archived individual samples (n = 3891) collected over a half century (1968–2018) throughout the North Atlantic to region of origin using range-wide genetic assignment. In the Northwest Atlantic, the distribution of assignments reinforced the importance of the Labrador Sea as an aggregation area, with 73% of all reporting groups detected. Moreover, individuals from six European reporting groups were identified in the Northwest Atlantic, and detections decreased with decreasing latitude spanning an area from Greenland to southern Newfoundland. In the Northeast Atlantic, six North American reporting groups were detected in samples from around the Faroe Islands. Based on the distribution of samples, estimates of trans-Atlantic migration distance averaged 3861 and 2889 km for North American and European salmon respectively. Our analysis highlights the widespread importance of the Labrador Sea and Faroe Islands to the species marine distribution patterns, and the prevalence of long-distance trans-Atlantic migration. Ultimately, the results suggest that environmental conditions experienced by many Atlantic salmon populations span much of the North Atlantic Ocean.

Beschreibung: Application of machine and deep learning methods in drug discovery and cancer research has gained a considerable amount of attention in the past years. As the field grows, it becomes crucial to systematically evaluate the performance of novel computational solutions in relation to established techniques. To this end, we compare rule-based and data-driven molecular representations in prediction of drug combination sensitivity and drug synergy scores using standardized results of 14 high-throughput screening studies, comprising 64 200 unique combinations of 4153 molecules tested in 112 cancer cell lines. We evaluate the clustering performance of molecular representations and quantify their similarity by adapting the Centered Kernel Alignment metric. Our work demonstrates that to identify an optimal molecular representation type, it is necessary to supplement quantitative benchmark results with qualitative considerations, such as model interpretability and robustness, which may vary between and throughout preclinical drug development projects.

Beschreibung: Motivation The investigation of quantitative trait loci (QTL) is an essential component in our understanding of how organisms vary phenotypically. However, many important crop species are polyploid (carrying more than two copies of each chromosome), requiring specialized tools for such analyses. Moreover, deciphering meiotic processes at higher ploidy levels is not straightforward, but is necessary to understand the reproductive dynamics of these species, or uncover potential barriers to their genetic improvement. Results Here, we present polyqtlR, a novel software tool to facilitate such analyses in (auto)polyploid crops. It performs QTL interval mapping in F1 populations of outcrossing polyploids of any ploidy level using identity-by-descent probabilities. The allelic composition of discovered QTL can be explored, enabling favourable alleles to be identified and tracked in the population. Visualization tools within the package facilitate this process, and options to include genetic co-factors and experimental factors are included. Detailed information on polyploid meiosis including prediction of multivalent pairing structures, detection of preferential chromosomal pairing and location of double reduction events can be performed. Availabilityand implementation polyqtlR is freely available from the Comprehensive R Archive Network (CRAN) at http://cran.r-project.org/package=polyqtlR. Supplementary information Supplementary data are available at Bioinformatics online.

Beschreibung: Circular RNAs (circRNAs) are widely expressed in highly diverged eukaryotes. Although circRNAs have been known for many years, their function remains unclear. Interaction with RNA-binding protein (RBP) to influence post-transcriptional regulation is considered to be an important pathway for circRNA function, such as acting as an oncogenic RBP sponge to inhibit cancer. In this study, we design a deep learning framework, CRPBsites, to predict the binding sites of RBPs on circRNAs. In this model, the sequences of variable-length binding sites are transformed into embedding vectors by word2vec model. Bidirectional LSTM is used to encode the embedding vectors of binding sites, and then they are fed into another LSTM decoder for decoding and classification tasks. To train and test the model, we construct four datasets that contain sequences of variable-length binding sites on circRNAs, and each set corresponds to an RBP, which is overexpressed in bladder cancer tissues. Experimental results on four datasets and comparison with other existing models show that CRPBsites has superior performance. Afterwards, we found that there were highly similar binding motifs in the four binding site datasets. Finally, we applied well-trained CRPBsites to identify the binding sites of IGF2BP1 on circCDYL, and the results proved the effectiveness of this method. In conclusion, CRPBsites is an effective prediction model for circRNA-RBP interaction site identification. We hope that CRPBsites can provide valuable guidance for experimental studies on the influence of circRNA on post-transcriptional regulation.

Beschreibung: Intratumoral heterogeneity is a well-documented feature of human cancers and is associated with outcome and treatment resistance. However, a heterogeneous tumor transcriptome contributes an unknown level of variability to analyses of differentially expressed genes (DEGs) that may contribute to phenotypes of interest, including treatment response. Although current clinical practice and the vast majority of research studies use a single sample from each patient, decreasing costs of sequencing technologies and computing power have made repeated-measures analyses increasingly economical. Repeatedly sampling the same tumor increases the statistical power of DEG analysis, which is indispensable toward downstream analysis and also increases one’s understanding of within-tumor variance, which may affect conclusions. Here, we compared five different methods for analyzing gene expression profiles derived from repeated sampling of human prostate tumors in two separate cohorts of patients. We also benchmarked the sensitivity of generalized linear models to linear mixed models for identifying DEGs contributing to relevant prostate cancer pathways based on a ground-truth model.

Beschreibung: Efforts to elucidate protein–DNA interactions at the molecular level rely in part on accurate predictions of DNA-binding residues in protein sequences. While there are over a dozen computational predictors of the DNA-binding residues, they are DNA-type agnostic and significantly cross-predict residues that interact with other ligands as DNA binding. We leverage a custom-designed machine learning architecture to introduce DNAgenie, first-of-its-kind predictor of residues that interact with A-DNA, B-DNA and single-stranded DNA. DNAgenie uses a comprehensive physiochemical profile extracted from an input protein sequence and implements a two-step refinement process to provide accurate predictions and to minimize the cross-predictions. Comparative tests on an independent test dataset demonstrate that DNAgenie outperforms the current methods that we adapt to predict residue-level interactions with the three DNA types. Further analysis finds that the use of the second (refinement) step leads to a substantial reduction in the cross predictions. Empirical tests show that DNAgenie’s outputs that are converted to coarse-grained protein-level predictions compare favorably against recent tools that predict which DNA-binding proteins interact with double-stranded versus single-stranded DNAs. Moreover, predictions from the sequences of the whole human proteome reveal that the results produced by DNAgenie substantially overlap with the known DNA-binding proteins while also including promising leads for several hundred previously unknown putative DNA binders. These results suggest that DNAgenie is a valuable tool for the sequence-based characterization of protein functions. The DNAgenie’s webserver is available at http://biomine.cs.vcu.edu/servers/DNAgenie/.

Beschreibung: Accurate prediction of immunogenic peptide recognized by T cell receptor (TCR) can greatly benefit vaccine development and cancer immunotherapy. However, identifying immunogenic peptides accurately is still a huge challenge. Most of the antigen peptides predicted in silico fail to elicit immune responses in vivo without considering TCR as a key factor. This inevitably causes costly and time-consuming experimental validation test for predicted antigens. Therefore, it is necessary to develop novel computational methods for precisely and effectively predicting immunogenic peptide recognized by TCR. Here, we described DLpTCR, a multimodal ensemble deep learning framework for predicting the likelihood of interaction between single/paired chain(s) of TCR and peptide presented by major histocompatibility complex molecules. To investigate the generality and robustness of the proposed model, COVID-19 data and IEDB data were constructed for independent evaluation. The DLpTCR model exhibited high predictive power with area under the curve up to 0.91 on COVID-19 data while predicting the interaction between peptide and single TCR chain. Additionally, the DLpTCR model achieved the overall accuracy of 81.03% on IEDB data while predicting the interaction between peptide and paired TCR chains. The results demonstrate that DLpTCR has the ability to learn general interaction rules and generalize to antigen peptide recognition by TCR. A user-friendly webserver is available at http://jianglab.org.cn/DLpTCR/. Additionally, a stand-alone software package that can be downloaded from https://github.com/jiangBiolab/DLpTCR.

Beschreibung: We present the second part of results of the ongoing project aimed at searching for and studying eXtremely Metal-Poor (XMP) – adopted as those with Zgas ≲ Z⊙/30, or with 12+log (O/H) ≲ 7.21 dex − very gas-rich blue dwarfs in voids. They were first identified in the course of the ‘unbiased’ study of the galaxy population in the nearby Lynx–Cancer void. These very rare and unusual galaxies seem to be the best proxies of so-called Very Young Galaxies (VYGs) defined recently in model simulations by Tweed et al. To date, for 16 pre-selected void XMP candidates, using the Big Telescope Alt-azimuth (BTA), the SAO 6-m telescope, we have obtained spectra suitable for the determination of O/H. For majority of the observed galaxies, the principal line [O iii] λ4363 , used for the direct classical Te method of O/H determination, is undetected. Therefore, to estimate O/H, we use a new ‘strong-lines’ method by Izotov et al. This appears to be the most accurate empirical O/H estimator for the range of 12+log (O/H) ≲ 7.4–7.5. For objects with higher O/H, we use the semi-empirical method by Izotov & Thuan with our modification accounting for variance of the excitation parameter O32. Six of those 16 candidates are found, with confidence, to be XMP dwarfs. In addition, eight studied galaxies are less metal-poor, with 12+log (O/H) = 7.24–7.33, and these can also fall into the category of VYG candidates. Taking into account our recently published work and the previously known (nine prototype galaxies) XMP gas-rich void objects, the new findings increase the number of this type of galaxy known to date to a total of 19.

Beschreibung: Heterotrophic protists are essential components of the marine ecosystem, yet they are often excluded from monitoring programmes. With limited resources, monitoring strategies need to be optimised considering both scientific knowledge and available resources. In doing so, it is crucial to understand how sampling frequency affects the value of the data. We analysed 11 years of weekly heterotrophic protist time-series data from Station L4 in the Western English Channel to explore how different sampling intervals impact data quality. In the L4 dataset, comprising 55 protist taxa, the reduction of sampling frequency from weekly to four times a year at specific seasons decreased the number of taxa encountered by 38% for ciliates and 29% for heterotrophic dinoflagellates while the mean annual biomass or its mean variation were not affected. Furthermore, when samples were taken only four times a year, biomass peaks of the ten most important taxa were often missed. The primary motivator for this study was furthering the development of the heterotrophic protist monitoring in temperate and subarctic marine areas, e.g. the Baltic Sea. Based on our findings, we give recommendations on sampling frequency to optimise the value of heterotrophic protist monitoring.

Beschreibung: Motivation The investigation of the structure of biological systems at the molecular level gives insights about their functions and dynamics. Shape and surface of biomolecules are fundamental to molecular recognition events. Characterizing their geometry can lead to more adequate predictions of their interactions. In the present work, we assess the performance of reference shape retrieval methods from the computer vision community on protein shapes. Results Shape retrieval methods are efficient in identifying orthologous proteins and tracking large conformational changes. This work illustrates the interest for the protein surface shape as a higher-level representation of the protein structure that (i) abstracts the underlying protein sequence, structure or fold, (ii) allows the use of shape retrieval methods to screen large databases of protein structures to identify surficial homologs and possible interacting partners and (iii) opens an extension of the protein structure–function paradigm toward a protein structure-surface(s)-function paradigm. Availabilityand implementation All data are available online at http://datasetmachat.drugdesign.fr. Supplementary information Supplementary data are available at Bioinformatics online.

Beschreibung: Fisheries management is usually supported by technical and financial measurements (i.e. logbooks and market data), which are helpful for ecological or economic assessments. Yet this information is not able to address social heterogeneity and fisher motivations, which are key to understanding fisher behaviour. This case study of the demersal segment in the Netherlands shows that combining quantitative analysis of logbooks with qualitative data collected by engaging with fishers can capture both fishing activity and its motivations, generating a more social understanding of fisher behaviour. A métier analysis of logbook data describes five dominant fishing practices among the selected segment. Twenty-five in-depth interviews with fishers along with focus groups including other experts identify three social factors that influence fisher behaviour in the Dutch demersal fleet: business structure, working rhythm, and polyvalence. The results show that motivations for fisher behaviour are more complex than complying with regulations or seeking profit: social factors also influence fishing activity. Furthermore, these social factors have real implications for the impacts of management measures on both the fishing communities and the environment, especially in times of change. These results are useful for management strategy development or evaluation because they are feasibly observable through existing data collection protocols.

Beschreibung: Harbour seals are surveyed aerially when they haul-out to moult in August. The proportion of the population hauled out throughout the year is related to temporal, environmental, and meteorological variables. Thus, monitoring is conducted under predefined ranges of conditions. Effects of variation within these ranges are rarely reviewed. We used linear models to assess effects of time, date and weather on the difference between counts predicted by a population growth model and observed counts, based on a 30-year time-series. Our top-ranked model explained 34.4% of the variance. Survey date and its interaction with survey year were the most important variables, with higher counts earlier in August, particularly early in the time series, where surveys may not have been timed optimally to capture the peak in the moult. Cloud cover, wind speed, temperature, and interactions between these were of lesser importance; there were fewer seals on land during cloudy, windy days and on clear, warm days. These effects of weather are likely related to temperature regulation. Power analyses suggested that correction for survey conditions would allow detection of a one percentage point annual change in population growth rate with 80% power 4 years sooner than without taking survey conditions into account.

Beschreibung: Sargassum forests play an important role in coastal waters as habitats for marine organisms, including commercial species. However, human activities have negatively affected their distribution causing a worldwide decline of Sargassum forests. Mapping and monitoring the distribution and biomass of these habitats using acoustic remote sensing techniques is key for their conservation. Nonetheless, most researches based on acoustic remote sensing methods focus on estimations of macrophyte area and its canopy height, and less researches reporting 3D visualization of these habitats. This study demonstrates the use of high-resolution multibeam echosounder (MBES) bathymetric data to visualize the 3D structure of Sargassum forests. Comparing acoustic data and underwater camera photos collected in field surveys, we identified Sargassum individuals as vertical clusters of contiguous sounding points with a base close to the sea bottom in the sounding data of the MBES. Using this criterion, we could distinguish Sargassum echoes, visualize the 3D structure of Sargassum forests and estimate the number of Sargassum individuals in the survey area. Using the relation between thallus length and dry weight of sampled Sargassum plants, standing stock and biomass could be estimated assuming the thallus length was the height of Sargassum plants identified with the MBES.

Beschreibung: This paper reviews application of the Daily Egg Production Method (DEPM) to sardine (Sardinops sagax) off southern Australia between 1995 and 2019. Coefficients of variation (CVs) of estimates of spawning biomass (SB) were reduced from 23–59% to 8–12% by: (i) estimating mean daily egg production (P0), spawning fraction (S), and sex ratio (R) from all historical data rather than annually; and (ii) combining batch fecundity (F) and female weight (W) into a single parameter, relative fecundity (F′ = F^/W). Total daily egg production was estimated most precisely from annual estimates of spawning area (A) and estimates of P0 obtained from historical data. Both S and R were estimated most precisely from historical data. Estimating W and F from historical data did not increase precision. F′ had lower CVs than both W and F, and was stable across years and a wide range of W. Findings demonstrate that A can be converted into a precise estimates of SB using estimates of P0, S, R, and F′ obtained from historical data. However, the possibility that DEPM parameters may change in the future cannot be discounted. Future monitoring should include annual estimation of P0 and periodic (e.g. 3–5 years) re-estimation of adult parameters.

Beschreibung: In his 1953 paper, Harald Sverdrup argued that the development of a spring bloom in the ocean depends on the juxtaposition of two depth horizons: the mixed-layer depth and the critical depth. Mixed-layer depth shallower than the critical depth favours phytoplankton growth in the layer and vice versa. However, mathematically, Sverdrup left the problem unsolved in the form of a transcendental equation. In spite of the high number of citations that this paper has garnered, the solution to this equation has not been found, until now. In this work, we present an analytical solution for the critical depth, as originally defined by Sverdrup. The paper opens with the definition of the critical depth and the description of the Lambert W function. The analytical solution for critical depth follows. Sverdrup’s original model is extended to include the effect of light attenuation by phytoplankton and the analytical solution for steady-state biomass in the mixed layer is derived. The expression for mixed-layer production at steady state is also presented. Two novel variants of the critical depth are defined: the optically uncoupled critical depth and the optically coupled critical depth. It is demonstrated that at steady state the optically coupled critical depth equals the mixed-layer depth and that the irradiance at the base of the mixed layer equals the irradiance at the optically uncoupled critical depth. Competitive exclusion is demonstrated to hold and the optically uncoupled critical depth is linked to the critical light intensity in multi-species competition. Finally, a conservation principle for the critical depth is found.

Beschreibung: Fishery bycatch poses a serious threat to seabird populations globally. Traditional haul-only post-capture observations are inadequate and inefficient to document seabird bycatch due to the substantial bycatch loss known to occur. Pre-capture observations offer an alternative by documenting seabird interactions leading up to captures. Based on the long-term large-scale dedicated field observations, this study revealed significant risk factors for the pre-capture stages of the seabird bycatch process in pelagic longline fisheries using Bayesian methods. Rough sea conditions were found to correlate with more seabirds following fishing vessels. Species identity, density effect, inter-species effect, and sea condition were found to significantly affect how frequently seabirds aggregated around a fishing vessel engage in bait-taking interactions. Intra-species competition was found to be the dominant type of density effect. Moreover, a web of inter-species interactions was identified to facilitate the bait-taking of superior competitors at the expense of inferior ones. The findings of this study are relevant to fishery managers in updating current data collection protocols to alleviate data issues caused by bycatch loss, to conservation biologists in quantifying bycatch risks for susceptible seabird populations, and in aiding the design and evaluation of bycatch mitigation measures.

Beschreibung: Fishing and translocation of marine species for use in aquaculture is widespread. Corkwing, goldsinny, and ballan wrasse (Symphodus melops, Ctenolabrus rupestris, and Labrus bergylta) are fished on the Swedish west coast for use as cleaner-fish in Norwegian salmon farms. Here, we aim to provide knowledge and recommendations to support ecosystem-based management for wrasse fisheries in Sweden. We compared fished and non-fished areas to test if current fishery levels have led to stock depletion. To gain insight on the role of wrasse in the algal belt trophic chain, we analysed the gut contents of goldsinny and corkwing using metabarcoding. Finally, we analysed the trophic interactions of wrasse and potential prey in a mesocosm study. We could not detect any signs of stock depletion or altered size structure in fished areas compared to the protected control area. Gut analyses confirmed both goldsinny and corkwing as non-specialized, omnivorous opportunists and revealed, with 189 prey taxa detected, a broader spectrum of prey than previously known. Common prey items included mesoherbivores such as small gastropods and crustaceans, but also insects and algae. We conclude that there are no visible signs of stock depletion at the current removal level of wrasses by the fishery. However, this emerging fishery should be closely monitored for potential cascading effects on the algal belt ecosystem, and our study could provide a baseline for future monitoring.

Beschreibung: This paper explores methods for constructing low multipole temperature and polarization likelihoods from maps of the cosmic microwave background anisotropies that have complex noise properties and partial sky coverage. We use Planck 2018 High Frequency Instrument (HFI) and updated SRoll2 temperature and polarization maps to test our methods. We present three likelihood approximations based on quadratic cross spectrum estimators: (i) a variant of the simulation-based likelihood (SimBaL) techniques used in the Planck legacy papers to produce a low multipole EE likelihood; (ii) a semi-analytical likelihood approximation (momento) based on the principle of maximum entropy; (iii) a density-estimation ‘likelihood-free’ scheme (delfi). Approaches (ii) and (iii) can be generalized to produce low multipole joint temperature-polarization (TTTEEE) likelihoods. We present extensive tests of these methods on simulations with realistic correlated noise. We then analyse the Planck data and confirm the robustness of our method and likelihoods on multiple inter- and intra-frequency detector set combinations of SRoll2 maps. The three likelihood techniques give consistent results and support a low value of the optical depth to reoinization, τ, from the HFI. Our best estimate of τ comes from combining the low multipole SRoll2momento (TTTEEE) likelihood with the CamSpec high multipole likelihood and is $au = 0.0627^{+0.0050}_{-0.0058}$. This is consistent with the SRoll2 team’s determination of τ, though slightly higher by ∼0.5σ, mainly because of our joint treatment of temperature and polarization.

Beschreibung: X-ray flares in gamma-ray bursts (GRBs) are believed to be generated by the late activities of central engine, and thus provide an useful tool to diagnose the properties of central objects. In this paper, we work on a GRB X-ray flare sample whose bulk Lorentz factors are constrained by two different methods and the jet opening angles are determined by the jet breaks in afterglow light curves. Considering a hyperaccreting stellar-mass black hole (BH) as the central engine of GRBs and the Blandford & Znajek process (BZ) as the jet production mechanism, we constrain the parameters of central engine by using the X-ray flare data. We find that the BZ mechanism is so powerful making it possible to interpret both GRB prompt emissions and bright X-ray flares. The wind parameter (p) and accreted mass (Md) fall into reasonable ranges. Our result is also applied to GRB 170817A. The late X-ray flare in GRB 170817A, if it is true, might not be a BH origin.

Beschreibung: Young neutron stars (NSs) have magnetic fields in the range 1012–1015 G, believed to be generated by dynamo action at birth. We argue that such a dynamo is actually too inefficient to explain the strongest of these fields. Dynamo action in the mature star is also unlikely. Instead we propose a promising new precession-driven dynamo and examine its basic properties, as well as arguing for a revised mean-field approach to NS dynamos. The precession-driven dynamo could also play a role in field generation in main-sequence stars.

Beschreibung: We report the discovery of a rare close binary system, SMSS J160639.78−100010.7, comprising a magnetic white dwarf with a field of about 30 MG and a brown dwarf. We measured an orbital period of 92 min which is consistent with the photometric period. Minimum and maximum light occur at the orbital quadratures Φ = 0.25 and 0.75, respectively, and cannot be caused by reflection on the brown dwarf, but, instead, by a spot on the synchronously rotating magnetic white dwarf. The brown dwarf does not fill its Roche lobe and the system may be in a low-accretion state or, more likely, in a detached state following episodes of mass transfer. SMSS J160639.78−100010.7 is the nearest known magnetic white dwarf plus brown dwarf system.

Beschreibung: In this work, we expand and test the capabilities of our recently developed superresolution (SR) model to generate high-resolution (HR) realizations of the full phase-space matter distribution, including both displacement and velocity, from computationally cheap low-resolution (LR) cosmological N-body simulations. The SR model enhances the simulation resolution by generating 512 times more tracer particles, extending into the deeply nonlinear regime where complex structure formation processes take place. We validate the SR model by deploying the model in 10 test simulations of box size 100 h−1 Mpc, and examine the matter power spectra, bispectra, and two-dimensional power spectra in redshift space. We find the generated SR field matches the true HR result at per cent level down to scales of k ∼ 10 h  Mpc−1. We also identify and inspect dark matter haloes and their substructures. Our SR model generates visually authentic small-scale structures that cannot be resolved by the LR input, and are in good statistical agreement with the real HR results. The SR model performs satisfactorily on the halo occupation distribution, halo correlations in both real and redshift space, and the pairwise velocity distribution, matching the HR results with comparable scatter, thus demonstrating its potential in making mock halo catalogues. The SR technique can be a powerful and promising tool for modelling small-scale galaxy formation physics in large cosmological volumes.

Beschreibung: Over the past decade, genome-wide assays for chromatin interactions in single cells have enabled the study of individual nuclei at unprecedented resolution and throughput. Current chromosome conformation capture techniques survey contacts for up to tens of thousands of individual cells, improving our understanding of genome function in 3D. However, these methods recover a small fraction of all contacts in single cells, requiring specialised processing of sparse interactome data. In this review, we highlight recent advances in methods for the interpretation of single-cell genomic contacts. After discussing the strengths and limitations of these methods, we outline frontiers for future development in this rapidly moving field.

Beschreibung: Good knowledge of a peptide’s tertiary structure is important for understanding its function and its interactions with its biological targets. APPTEST is a novel computational protocol that employs a neural network architecture and simulated annealing methods for the prediction of peptide tertiary structure from the primary sequence. APPTEST works for both linear and cyclic peptides of 5–40 natural amino acids. APPTEST is computationally efficient, returning predicted structures within a number of minutes. APPTEST performance was evaluated on a set of 356 test peptides; the best structure predicted for each peptide deviated by an average of 1.9Å from its experimentally determined backbone conformation, and a native or near-native structure was predicted for 97% of the target sequences. A comparison of APPTEST performance with PEP-FOLD, PEPstrMOD and PepLook across benchmark datasets of short, long and cyclic peptides shows that on average APPTEST produces structures more native than the existing methods in all three categories. This innovative, cutting-edge peptide structure prediction method is available as an online web server at https://research.timmons.eu/apptest, facilitating in silico study and design of peptides by the wider research community.

Beschreibung: Target strength model inputs including morphometry, material properties, lipid composition, and in situ orientations were measured for sub-Arctic krill (Euphausia pacifica, Thysanoessa spinifera, T. inermis, and T. raschii) in the eastern Bering Sea (EBS, 2016) and Gulf of Alaska (GOA, 2017). Inter-species and -regional animal lengths were significantly different (F1,680 = 114.10, p 

Beschreibung: Atlantic cod (Gadus morhua) is one of the most commercially important fish species in the North Atlantic. Environmental factors, such as water temperatures, influence growth of individuals over time, thus forming population-specific growth patterns across climatic regions. Here we develop an integrative approach to investigate the role of temperature in shaping geographic differences of cod growth in the Celtic Sea, North Sea, Iceland, and Barents Sea. We combine a physiology-based growth model and 50-years observational temperature data of 0.5 × 0.5° spatial resolution to simulate continuous growth of cod. The model generated weight-at-age data for the period 1959–2007 which we compared to observational data from fishery-independent scientific surveys. In the Celtic and the northern North Sea, simulated growth matches well observational data. We also show that relatively warm temperatures in the Celtic Sea facilitate maximum growth rates; future warming is likely to have a negative impact on growth of these cod stocks. Growth simulations in Icelandic waters and the Barents Sea are less consistent with local observational data. More complex growth patterns in these regions are probably shaped by ontogenetic shifts in temperature regimes, feeding conditions and physiological adaptations. These findings should stimulate further research on critical processes to be considered in population-specific projections of growth of cod and productivity.

Beschreibung: In designing and performing surveys of animal abundance, monitoring programs often struggle to determine the sampling intensity and design required to achieve their objectives, and this problem greatly increases in complexity for multispecies surveys with inherent trade-offs among species. To address these issues, we conducted a multispecies stratified random survey design optimization using a spatiotemporal operating model and a genetic algorithm that optimizes both the stratification (defined by depth and longitude) and the minimum optimal allocation of samples across strata subject to prespecified precision limits. Surveys were then simulated under those optimized designs and performance was evaluated by calculating the precision and accuracy of a resulting design-based abundance index. We applied this framework to a multispecies fishery-independent bottom trawl survey in the Gulf of Alaska, USA. Incorporating only spatial variation in the optimization failed to produce population estimates within the prespecified precision constraints, whereas including additional spatiotemporal variation ensured that estimates were both unbiased and within prespecified precision constraints. In general, results were not sensitive to the number of strata in the optimized solutions. This optimization approach provides an objective quantitative framework for designing new, or improving existing, survey designs for many different ecosystems.

Beschreibung: Changes in animal movement and behaviour at fine scales (tens of metres) in immediate proximity to tidal stream turbine structures are largely unknown and have implications for risks of animal collision with turbine blades. This study used upward-facing multibeam echosounder data to detect and track animal movement comprising fish, diving seabirds, and marine mammals. Measurements over spring-neap tidal cycles at a turbine structure (no blades present) are compared to a neighbouring reference area with no structure and comparable conditions, with measurements consecutive in time to maximize comparability. The majority of tracked animals (93.4% around turbine structure and 99.1% without turbine structure) were observed swimming against the flow, with 87.5% and 97.8%, respectively, making ground and showing capability of manoeuvring in tidal stream flow speeds. Track tortuosity increased around the turbine structure compared to the reference site, particularly in the wake and at low flow speeds, indicating animal station-holding or milling behaviour. These data also evidence the benefits of multibeam echosounders to measure animal movement through larger measurement volumes rather than relying on single-beam echosounders to measure animal presence alone, including to avoid large biases overestimating the size of schools swimming against the flow measured by time-in-beam.

Beschreibung: The Antarctic Peninsula is one of the most rapidly warming regions on earth, and it is likely that the abundance and distribution of marine predators will change as a result.Procellariiform seabirds are highly mobile predators, which target specific habitat characteristics associated with underlying distributions of prey and areas of increased prey availability. We use ship surveys and hurdle models, to estimate the summer distribution and relative density of 11 seabird species within the northern Antarctic Peninsula marine ecosystem. Models differed among species; however, sea surface temperature and depth were frequently associated with seabird occurrence and had the greatest explanatory power across many species. Null models based on observation data were better at predicting seabird density than models that included environmental covariates. This suggests that the main driver of distribution patterns is the broad-scale habitat features, and fine-scale aggregations within these ranges are harder to predict. Our seabird distribution models reflect known habitat associations, species hotspots, and community organization relative to oceanic and coastal marine processes. Application of species distribution models will benefit the assessments of critical habitat and potential responses to climate change and anthropogenic disturbance, which will provide insight into how species may change in polar ecosystems.

Beschreibung: Many data-limited fish stocks worldwide require management advice. Simple empirical management procedures have been used to manage data-limited fisheries but do not necessarily ensure compliance with maximum sustainable yield objectives and precautionary principles. Genetic algorithms are efficient optimization procedures for which the objectives are formalized as a fitness function. This optimization can be included when testing management procedures in a management strategy evaluation. This study explored the application of a genetic algorithm to an empirical catch rule and found that this approach could substantially improve the performance of the catch rule. The optimized parameterization and the magnitude of the improvement were dependent on the specific stock, stock status, and definition of the fitness function. The genetic algorithm proved to be an efficient and automated method for tuning the catch rule and removed the need for manual intervention during the optimization process. Therefore, we conclude that the approach could also be applied to other management procedures, case-specific tuning, and even data-rich stocks. Finally, we recommend the phasing out of the current generic ICES “2 over 3” advice rule in favour of case-specific catch rules of the form tested here, although we caution that neither works well for fast-growing stocks.

Beschreibung: Ocean grabbing occurs when traditional users, such as small-scale fishers, are pushed aside by new development activities. This grabbing must be prevented to avoid sea uses that maintain or increase social inequity. In this paper, we show that in tropical Atlantic countries, such as Brazil and Senegal, examples of ocean grabbing already occur. In this context, we analyse if Maritime Spatial Planning (MSP) may be an opportunity to limit ocean grabbing or, to the contrary, poses a risk to increase it. MSP calls for an ecosystem approach that requires integrated coastal and marine management and involves stakeholders in developing a shared vision of the future, where society and environment are preserved. However, recent studies have shown that MSP is a process to be used cautiously to ensure equitable decisions. Meanwhile, the concept is spreading worldwide including in tropical Atlantic countries. We highlight that context matters and the specificities of the tropical Atlantic must be taken into account when deploying MSP processes. In the tropical Atlantic context, there is increased imbalances of stakeholder power, traps from decision support tools, and a need for adaptive management. These specific features must be addressed when deploying MSP in a way to avoid ocean grabbing.

Beschreibung: The influence of oceanographic variables on assemblages of meso- and bathypelagic fish was investigated along a Northeastern Atlantic Ocean transect (Cape Verde to the Bay of Biscay) during May 2019. Fish were collected using a mrozooplankton trawl during daylight hours at ten stations. Along the transect, 17 hydrographic stations were also performed with a CTD (Conductivity, Temperature, and Depth). A total of 130 fish taxa were identified. The dominant family was Gonostomatidae, with four species (Cyclothone braueri, Cyclothone microdon, Cyclothone pseudopallida, and Cyclothone pallida) being responsible of more than 78% of the total density. The most frequent species that appeared to be ubiquitous were C. braueri and C. pseudopallida, while Myctophidae was the most diverse family. Multivariate analyses revealed two clusters related with the latitudinal gradient. The fish community in the southern stations (25–37°N) was more diverse than in the northern stations (42–48°N). Temperature from 300 to 700 m depth explained 65% of variation in terms of density and 58% in terms of biomass, both statistically significant. The investigated variation in the deep-pelagic ecosystems on a large spatial scale gives essential information to ecosystem management approhes and marine spatial planning.

Beschreibung: In the European Union (EU), the ecosystem approach to fisheries is implemented through several directives and polices, which are overarched by the Marine Strategy Framework Directive (MSFD). The MSFD requires the assessment of the environmental status of exploited fish and shellfish stocks, among others, to advise the EU Common Fisheries Policy (CFP) on sustainable catch options. The stock assessments for the CFP and the MSFD are supposed to be congruent, yet they differ in several substantial subjects. While the annual CFP assessment is based on two indicators, fishing mortality (F) and spawning stock biomass (SSB), in its Descriptor 3 (D3) the MSFD requires the assessment of three criteria (F, SSB, and age or size structure) within a six year period. Assessing exploited fish and shellfish stocks under the MSFD by using existing CFP assessments has therefore been a half-completed task, which had to be fulfilled by the member states of the EU. This paper suggests six easy steps, referred to as “quick wins” (QW), which are based on existing information from CFP stock assessments. The implementation of these six QW would allow for assessments of exploited fish and shellfish stocks that are compliant to D3. These QW are to (i) assess length/age structure within a stock, (ii) analyse and assess selectivity by fisheries, (iii) use all available information to assess a stock, (iv) use response indicators to assess environmental targets, (v) provide integrated stock-specific advice, and (vi) provide assessments with a mid-term perspective over a-period of six years. International Council for the Exploration of the Sea (ICES) combines the infrastructure and expertise to produce stock-specific assessment products tailored to the requirements of MSFD D3. Thereby, ICES could provide a valuable service to EU member states in the north-east Atlantic region by providing scientifically validated, quality-assured, and MSFD-compliant single-stock assessment products. This would be a big advancement towards implementing the ecosystem approach to fisheries management within Europe.

Beschreibung: Deep generative models have been an upsurge in the deep learning community since they were proposed. These models are designed for generating new synthetic data including images, videos and texts by fitting the data approximate distributions. In the last few years, deep generative models have shown superior performance in drug discovery especially de novo molecular design. In this study, deep generative models are reviewed to witness the recent advances of de novo molecular design for drug discovery. In addition, we divide those models into two categories based on molecular representations in silico. Then these two classical types of models are reported in detail and discussed about both pros and cons. We also indicate the current challenges in deep generative models for de novo molecular design. De novo molecular design automatically is promising but a long road to be explored.

Beschreibung: DNA methylation may be regulated by genetic variants within a genomic region, referred to as methylation quantitative trait loci (mQTLs). The changes of methylation levels can further lead to alterations of gene expression, and influence the risk of various complex human diseases. Detecting mQTLs may provide insights into the underlying mechanism of how genotypic variations may influence the disease risk. In this article, we propose a methylation random field (MRF) method to detect mQTLs by testing the association between the methylation level of a CpG site and a set of genetic variants within a genomic region. The proposed MRF has two major advantages over existing approaches. First, it uses a beta distribution to characterize the bimodal and interval properties of the methylation trait at a CpG site. Second, it considers multiple common and rare genetic variants within a genomic region to identify mQTLs. Through simulations, we demonstrated that the MRF had improved power over other existing methods in detecting rare variants of relatively large effect, especially when the sample size is small. We further applied our method to a study of congenital heart defects with 83 cardiac tissue samples and identified two mQTL regions, MRPS10 and PSORS1C1, which were colocalized with expression QTL in cardiac tissue. In conclusion, the proposed MRF is a useful tool to identify novel mQTLs, especially for studies with limited sample sizes.

Beschreibung: Hydrogen cyanide (HCN) and its isomer hydrogen isocyanide (HNC) play an important role in molecular cloud chemistry and the formation of more complex molecules. We investigate here the impact of protostellar shocks on the HCN and HNC abundances from high-sensitivity IRAM 30 m observations of the prototypical shock region L1157-B1 and the envelope of the associated Class 0 protostar, as a proxy for the pre-shock gas. The isotopologues H12CN, HN12C, H13CN, HN13C, HC15N, H15NC, DCN, and DNC were all detected towards both regions. Abundances and excitation conditions were obtained from radiative transfer analysis of molecular line emission under the assumption of local thermodynamical equilibrium. In the pre-shock gas, the abundances of the HCN and HNC isotopologues are similar to those encountered in dark clouds, with an HCN/HNC abundance ratio ≈1 for all isotopologues. A strong D-enrichment (D/H ≈ 0.06) is measured in the pre-shock gas. There is no evidence of 15N fractionation neither in the quiescent nor in the shocked gas. At the passage of the shock, the HCN and HNC abundances increase in the gas phase in different manners so that the HCN/HNC relative abundance ratio increases by a factor 20. The gas-grain chemical and shock model uclchem allows us to reproduce the observed trends for a C-type shock with pre-shock density n(H) = $10^5hbox{cm$^{-3}$}$ and shock velocity $V_mathrm{ s}= 40hbox{kms$^{-1}$}$. We conclude that the HCN/HNC variations across the shock are mainly caused by the sputtering of the grain mantle material in relation with the history of the grain ices.

Beschreibung: Motivation The mathematically optimal solution in computational protein folding simulations does not always correspond to the native structure, due to the imperfection of the energy force fields. There is therefore a need to search for more diverse suboptimal solutions in order to identify the states close to the native. We propose a novel multimodal optimization protocol to improve the conformation sampling efficiency and modeling accuracy of de novo protein structure folding simulations. Results A distance-assisted multimodal optimization sampling algorithm, MMpred, is proposed for de novo protein structure prediction. The protocol consists of three stages: The first is a modal exploration stage, in which a structural similarity evaluation model DMscore is designed to control the diversity of conformations, generating a population of diverse structures in different low-energy basins. The second is a modal maintaining stage, where an adaptive clustering algorithm MNDcluster is proposed to divide the populations and merge the modal by adjusting the annealing temperature to locate the promising basins. In the last stage of modal exploitation, a greedy search strategy is used to accelerate the convergence of the modal. Distance constraint information is used to construct the conformation scoring model to guide sampling. MMpred is tested on a large set of 320 non-redundant proteins, where MMpred obtains models with TM-score≥0.5 on 291 cases, which is 28% higher than that of Rosetta guided with the same set of distance constraints. In addition, on 320 benchmark proteins, the enhanced version of MMpred (E-MMpred) has 167 targets better than trRosetta when the best of five models are evaluated. The average TM-score of the best model of E-MMpred is 0.732, which is comparable to trRosetta (0.730). Availability and implementation The source code and executable are freely available at https://github.com/iobio-zjut/MMpred. Supplementary information Supplementary data are available at Bioinformatics online.

Beschreibung: We present results from a polarization study of the radio-intermediate quasar, III Zw 2, at a redshift of 0.089, with the upgraded Giant Metrewave Radio Telescope (uGMRT) at 685 MHz and the Karl G. Jansky Very Large Array (VLA) at 5 and 34 GHz. We detect a kpc-scale outflow, exhibiting transverse magnetic (B-) fields. The curved jet terminates in a bow-shock-like radio structure with inferred B-fields aligned with the lobe edges. We suggest that the radio outflow in III Zw 2 is a combination of a collimated jet along with a wind-like component. This ‘wind’ component could be a magnetized accretion disc wind or the outer layers of a broadened jet or a combination of both. The current data cannot differentiate between these possibilities. We also detect kpc-scale lobe emission that is misaligned with the primary lobes in the uGMRT images. The spectral indices and the electron lifetimes in the misaligned lobe are similar to the primary lobe, suggesting that the misaligned lobe is not a relic. We propose that changing spectral states of the accretion disc, and the subsequent intermittent behaviour of the outflow, along with the close interplay between the jet and ‘wind’ could explain the radio-intermediate nature of III Zw 2. Our study shows that radio-intermediate quasars are promising sources for understanding the role of jets and winds in galaxy evolution and demonstrates the power of radio polarization studies towards achieving this.

Beschreibung: Managers, stakeholders, and scientists recognize the need for collaborative, transparent, integrated approaches to complex resource management issues, and frameworks to address these complex issues are developing. Through the course of 2019, the Mid-Atlantic Fishery Management Council developed a conceptual model of ecosystem linkages and risks for summer flounder, a species of recreational and commercial fisheries importance. The proximal aim of the model was to develop a list of integrated management questions that could be refined and addressed through a future quantitative management strategy evaluation. As such, this conceptual model served as a scoping tool. However, the true value of the conceptual model lays elsewhere: familiarizing resource managers historically focused on single-species management with the potential utility of an ecosystem approach to management. This paper details the goals and development of the conceptual model and situates this process in the broader context of best practices for collaborative open science and scientific reproducibility. Further, it highlights a successful path by which the shift towards ecosystem-based management can be actuated.

Beschreibung: Understanding the competitive interactions of ecological similar species is essential to determine their role and niche in the ecosystem. Using both conventional dietary methods and stable isotope analysis, we examined the feeding ecology, trophic position and possible niche partitioning of Northeast Atlantic mackerel (Scomber scombrus), Icelandic summer spawning (ISS) herring (Clupea harengus), and Norwegian spring spawning (NSS) herring from Icelandic waters during the spring and summer 2012 and 2014. The stomach analysis showed differences in prey preferences among the species during summer, where mackerel diets were almost exclusively calanoid copepods, while herring ate larger zooplankton, i.e. euphausiids and amphipods. Analysis of isotopic diet contribution of mackerel and herring using Bayesian mixing models, representing mainly the spring consumption, revealed that euphausiids were the main dietary contributor of both mackerel and herring. Fish prey were also an important dietary source for ISS herring, and amphipods were important to the diet of NSS herring. In addition, mackerel encompassed a broad isotopic niche, whereas the two herring stocks had narrower and separated niches from each other. The results from this analysis reveal new and holistic information into the diet of these species, which have various ecological implications.

Beschreibung: Theoretical size at the onset of maturity (TSOM) for female Norway lobster was estimated by a new methodology based on the probability distributions of mature individuals built on physiological maturity measures. Onset of maturity using TSOM varied from 18.4- to 33.7-mm carapace length for the Irish functional management units (FUs). These estimates showed a significant negative linear relationship (R2 = 0.60) with population density and a significant positive linear relationship with average size in females (R2 = 0.84). The size class at which 50% of the females are sexually mature (L50) was linked to the new TSOM metric by a significant positive linear relationship (R2 = 0.40). This set of linear relationships ultimately allowed TSOM and L50 to be estimated without a requirement for maturity stages to be distinguished. As well as contributing to the stock assessment and management of Nephrops (e.g. in data-limited FUs) and its potential for the calibration of more routinely used estimates, TSOM might be applied in new species and meta-analyses where size of maturity data are scarce. This new metric also better-defines the maturity process since, taken together, TSOM, L50, and smallest berried female represent sequential maturity events: (i) onset of maturity, (ii) 50% mature (from gonad staging), and (iii) berried females.

Beschreibung: Although gelatinous zooplankton are an important component of marine ecosystems, gelatinous mesozooplankton that are

Beschreibung: The carapace length (CL) at which American lobster (Homarus americanus) females reach maturity can be used to evaluate egg production, growth patterns, and the overall health of lobster stocks. The female maturity datasets used to represent Gulf of Maine (GOM) lobsters in the 2015 Atlantic States Marine Fisheries Commission American Lobster Stock Assessment were collected in the 1990s by the Maine Department of Marine Resources at two coastal sites. Many studies have demonstrated an inverse relationship between temperature and the size at maturity in female lobsters, and GOM waters have warmed significantly over this period. To update these GOM maturity datasets, we used ovarian staging to determine the maturity status of over 1200 females from fives sites over 3 years. Broad application of this methodology in tandem with key growth measurements on females 50–120 mm CL allowed us to characterize reproductive development and generate maturity ogives (proportion mature at a given CL). We observed a latitudinal gradient in the size at maturity across this coastal region of the GOM and quantified a decrease in this size over 25 years. These findings have implications for future stock assessment approaches and management measures implemented to sustain this valuable fishery.

Beschreibung: The carbonate chemistry in coastal waters is more variable compared with that of open oceans, both in magnitude and time scale of its fluctuations. However, knowledge of the responses of coastal phytoplankton to dynamic changes in pH/pCO2 has been scarcely documented. Hence, we investigated the physiological performance of a coastal isolate of the coccolithophore Emiliania huxleyi (PML B92/11) under fluctuating and stable pCO2 regimes (steady ambient pCO2, 400 μatm; steady elevated pCO2, 1200 μatm; diurnally fluctuating elevated pCO2, 600–1800 μatm). Elevated pCO2 inhibited the calcification rate in both the steady and fluctuating regimes. However, higher specific growth rates and lower ratios of calcification to photosynthesis were detected in the cells grown under diurnally fluctuating elevated pCO2 conditions. The fluctuating pCO2 regime alleviated the negative effects of elevated pCO2 on effective photochemical quantum yield and relative photosynthetic electron transport rate compared with the steady elevated pCO2 treatment. Our results suggest that growth of E. huxleyi could benefit from diel fluctuations of pH/pCO2 under future-projected ocean acidification, but its calcification was reduced by the fluctuation and the increased concentration of CO2, reflecting a necessity to consider the influences of dynamic pH fluctuations on coastal carbon cycles associated with ocean global changes.

Beschreibung: Millions of people’s livelihoods rely on artisanal fisheries. However, in many regions fishers are increasingly facing ciguatera poisoning, a seafood-borne illness. The toxin, produced by benthic dinoflagellates, can spread through marine food webs and to humans by direct consumption. Ciguatera risk can play a major role in fisher’s activities but has never been considered in any marine spatial plans thus far. To fill this gap, we examined if integrating ciguatera in systematic conservation plans could affect these decisions. We developed through map-based interviews, a novel seven-step framework to collect and map local knowledge on ciguatera risk and fisheries activities with two innovations: (i) better mapping of fishing grounds by combining geomorphological habitat and fishing gear information, and (ii) integrating ciguatera risk directly into systematic spatial planning designs and scenarios conceived to maximize benthic habitat conservation while minimizing impacts to fishers. The approach is illustrated for Raivavae Island, in French Polynesia, Pacific Ocean. We found that integrating ciguatera significantly improved prioritization solutions with a 24–38% decrease of costs to fishers compared with scenarios based solely on fishery data. This framework was designed for scientists and managers to optimize the implementation of conservation plans and could be generalized to ciguatera-prone areas.

Beschreibung: Deep-sea marine fishes support important fisheries but estimates of their distributions are often incomplete as the data behind them may reflect fishing practices, access rights, or political boundaries, rather than actual geographic distributions. We use a simple suitable habitat model based on bottom depth, temperature, and salinity to estimate the potential distribution of Greenland halibut (Reinhardtius hippoglossoides). A large presence-only dataset is examined using multivariate kernel densities to define environmental envelopes, which we link to spatial distribution using a pan-Arctic oceanographic model. Occurrences generally fit the model well, although there were gaps in the predicted circum-Arctic distribution likely due to limited survey activity in many of the ice-covered seas around the Arctic Ocean. Bottom temperature and depth were major factors defining model fit to observations, but other factors, such as ecosystem interactions and larval drift could also influence distribution. Model predictions can be tested by increasing sampling effort in poorly explored regions and by studying the connectivity of putative populations. While abundances of Greenland halibut in the High Arctic are currently low, some areas are predicted to be suitable habitat for this species, suggesting that on-going sea-ice melt may lead to fisheries expansion into new areas.

Beschreibung: The peripheral regions of globular clusters (GCs) are extremely challenging to study due to their low surface brightness nature and the dominance of Milky Way contaminant populations along their sightlines. We have developed a probabilistic approach to this problem through utilizing a mixture model in spatial and proper motion space that separately models the cluster, extra-tidal, and contaminant stellar populations. We demonstrate the efficacy of our method through application to Gaia third early data release photometry and astrometry in the direction of NGC 5139 (ω Centauri, hereafter ω Cen), a highly challenging target on account of its Galactic latitude (b ≈ 15○) and low proper motion contrast with the surrounding field. We recover the spectacular tidal extensions, spanning the 10○ on the sky explored here, seen in earlier work and quantify the star count profile and ellipticity of the system out to a cluster-centric radius of 4○. We show that both RR Lyrae and blue horizontal branch stars consistent with belonging to ω Cen are found in the tidal tails, and calculate that these extensions contain at least ≈0.1 per cent of the total stellar mass in the system. Our high-probability members provide prime targets for future spectroscopic studies of ω Cen out to unprecedented radii.

Beschreibung: Antimicrobial resistance (AMR) poses a threat to global public health. To mitigate the impacts of AMR, it is important to identify the molecular mechanisms of AMR and thereby determine optimal therapy as early as possible. Conventional machine learning-based drug-resistance analyses assume genetic variations to be homogeneous, thus not distinguishing between coding and intergenic sequences. In this study, we represent genetic data from Mycobacterium tuberculosis as a graph, and then adopt a deep graph learning method—heterogeneous graph attention network (‘HGAT–AMR’)—to predict anti-tuberculosis (TB) drug resistance. The HGAT–AMR model is able to accommodate incomplete phenotypic profiles, as well as provide ‘attention scores’ of genes and single nucleotide polymorphisms (SNPs) both at a population level and for individual samples. These scores encode the inputs, which the model is ‘paying attention to’ in making its drug resistance predictions. The results show that the proposed model generated the best area under the receiver operating characteristic (AUROC) for isoniazid and rifampicin (98.53 and 99.10%), the best sensitivity for three first-line drugs (94.91% for isoniazid, 96.60% for ethambutol and 90.63% for pyrazinamide), and maintained performance when the data were associated with incomplete phenotypes (i.e. for those isolates for which phenotypic data for some drugs were missing). We also demonstrate that the model successfully identifies genes and SNPs associated with drug resistance, mitigating the impact of resistance profile while considering particular drug resistance, which is consistent with domain knowledge.

Beschreibung: Many astrophysical hydrodynamics simulations must account for gravity, and evaluating the gravitational field at the positions of all resolution elements can incur significant cost. Typical algorithms update the gravitational field at the position of each resolution element every time the element is updated hydrodynamically, but the actual required update frequencies for hydrodynamics and gravity can be different in general. We show that the gravity calculation in hydrodynamics simulations can be optimized by only updating gravity on a time-scale dictated by the already determined maximum time-step for accurate gravity integration Δtgrav, while staying well within the typical error budget of hydro schemes and gravity solvers. Our implementation in the gizmo code uses the time-scale derived from the tidal tensor $t_{ m tidal} = Vert mathbf {T}Vert ^{-1/2}$ to determine Δtgrav and the force update frequency in turn, and uses the rate of change of acceleration evaluated by the gravity solver to construct a predictor of the acceleration for use between updates. We test the scheme on standard self-gravitating hydrodynamics test problems, finding solutions very close to the standard scheme while evaluating far fewer gravity forces, optimizing the simulations. We also demonstrate a $sim 70{{ m per cent}}$ speed-up in an example simulation of a giant molecular cloud. In general, this scheme introduces a new tunable parameter for obtaining an optimal compromise between accuracy and computational cost, in conjunction with, e.g. time-step tolerance, numerical resolution, and gravity solver tolerance.

Beschreibung: The rotation powered pulsar loses angular momentum at a rate of the rotation power divided by the angular velocity Ω*. This means that the length of the lever arm of the angular momentum extracted by the photons, relativistic particles, and wind must be on average c/Ω*, which is known as the light cylinder radius. Therefore, any deposition of the rotation power within the light cylinder causes insufficient loss of angular momentum. In this paper, we investigate two cases of this type of energy release: polar cap acceleration and Ohmic heating in the magnetospheric current inside the star. As for the first case, the outer magnetosphere beyond the light cylinder is found to compensate the insufficient loss of the angular momentum. We argue that the energy flux coming from the sub-rotating magnetic field lines must be larger than the solid-angle average value, and as a result, an enhanced energy flux emanating beyond the light cylinder is observed in different phases in the light curve from those of emission inside the light cylinder. As for the second case, the stellar surface rotates more slowly than the stellar interior. We find that the way the magnetospheric current closes inside the star is linked to how the angular momentum is transferred inside the star. We obtain numerical solutions that shows that the magnetospheric current inside the star spreads over the polar cap magnetic flux embedded in the star in such a way that electromotive force is gained efficiently.

Beschreibung: The structure and evolution of protoplanetary discs (PPDs) are largely governed by disc angular momentum transport, mediated by magnetic fields. In the most observable outer disc, ambipolar diffusion is the primary non-ideal magnetohydrodynamic (MHD) effect. In this work, we study the gas dynamics in outer PPDs by conducting a series of global three-dimensional non-ideal MHD simulations with ambipolar diffusion and net poloidal magnetic flux, using the Athena++ MHD code, with resolution comparable to local simulations. Our simulations demonstrate the coexistence of magnetized disc winds and turbulence driven by the magneto-rotational instability (MRI). While MHD winds dominate disc angular momentum transport, the MRI turbulence also contributes significantly. We observe that magnetic flux spontaneously concentrates into axisymmetric flux sheets, leading to radial variations in turbulence levels, stresses, and accretion rates. Annular substructures arise as a natural consequence of magnetic flux concentration. The flux concentration phenomena show diverse properties with different levels of disc magnetization and ambipolar diffusion. The disc generally loses magnetic flux over time, though flux sheets could prevent the leak of magnetic flux in some cases. Our results demonstrate the ubiquity of disc annular substructures in weakly MRI turbulent outer PPDs and imply a stochastic nature of disc evolution.

Beschreibung: Motivated by recent theoretical work on tidal disruption events and other peculiar transients, we present moving-mesh radiation-hydrodynamic simulations of radiative luminosity emitted by a central source being reprocessed by a wind-like outflow. We couple the moving-mesh hydrodynamic code jet with our newly developed radiation module based on mixed-frame grey flux-limited diffusion with implicit timestep update. This allows us to study the self-consistent multidimensional radiation-hydrodynamic evolution over more than 10 orders of magnitude in both space and time in a single run. We simulate an optically thick spherical wind with constant or evolving mass-loss rate, which is irradiated by a central isotropic or angularly dependent radiation source. Our spherically symmetric simulations confirm previous analytic results by identifying different stages of radiation reprocessing: radiation trapped in the wind, diffusing out through the wind, and reaching constant maximum attenuation. We find that confining the central radiation source in a cone with moderate opening angles decrease up to one order of magnitude the early flux along sightlines oriented away from the direction of radiation injection but that the reprocessed radiation becomes isotropic roughly after one lateral diffusion time through the ejecta. We discuss further applications and guidelines for the use of our novel radiation-hydrodynamics tool in the context of transient modelling.

Beschreibung: Protein engineering and design principles employing the 20 standard amino acids have been extensively used to achieve stable protein scaffolds and deliver their specific activities. Although this confers some advantages, it often restricts the sequence, chemical space, and ultimately the functional diversity of proteins. Moreover, although site-specific incorporation of non-natural amino acids (nnAAs) has been proven to be a valuable strategy in protein engineering and therapeutics development, its utility in the affinity-maturation of nanobodies is not fully explored. Besides, current experimental methods do not routinely employ nnAAs due to their enormous library size and infinite combinations. To address this, we have developed an integrated computational pipeline employing structure-based protein design methodologies, molecular dynamics simulations and free energy calculations, for the binding affinity prediction of an nnAA-incorporated nanobody toward its target and selection of potent binders. We show that by incorporating halogenated tyrosines, the affinity of 9G8 nanobody can be improved toward epidermal growth factor receptor (EGFR), a crucial cancer target. Surface plasmon resonance (SPR) assays showed that the binding of several 3-chloro-l-tyrosine (3MY)-incorporated nanobodies were improved up to 6-fold into a picomolar range, and the computationally estimated binding affinities shared a Pearson’s r of 0.87 with SPR results. The improved affinity was found to be due to enhanced van der Waals interactions of key 3MY-proximate nanobody residues with EGFR, and an overall increase in the nanobody’s structural stability. In conclusion, we show that our method can facilitate screening large libraries and predict potent site-specific nnAA-incorporated nanobody binders against crucial disease-targets.

Beschreibung: Over the past few years, meta-analysis has become popular among biomedical researchers for detecting biomarkers across multiple cohort studies with increased predictive power. Combining datasets from different sources increases sample size, thus overcoming the issue related to limited sample size from each individual study and boosting the predictive power. This leads to an increased likelihood of more accurately predicting differentially expressed genes/proteins or significant biomarkers underlying the biological condition of interest. Currently, several meta-analysis methods and tools exist, each having its own strengths and limitations. In this paper, we survey existing meta-analysis methods, and assess the performance of different methods based on results from different datasets as well as assessment from prior knowledge of each method. This provides a reference summary of meta-analysis models and tools, which helps to guide end-users on the choice of appropriate models or tools for given types of datasets and enables developers to consider current advances when planning the development of new meta-analysis models and more practical integrative tools.

Beschreibung: Since the detection of X-ray pulses from ultraluminous X-ray sources (ULXs) in 2014, neutron stars have been considered as their central objects. However, it remains unclear how neutron stars can be brighter than the Eddington luminosity, and no unified view exists on the magnetic field of neutron stars and the degree of beaming. Recent observations suggest that some X-ray pulsating ULXs have Be-type donors, and some of them occupy the same region as Be-type high-mass X-ray binaries (Be-HMXBs) on the Corbet diagram, which reveals the relation between spin and orbital periods. This suggests that at least some ULXs are special cases of Be-HMXBs. In this study, we use the framework of mass accretion models for Be-HMXBs to investigate the conditions under which neutron stars achieve mass accretion rates beyond the Eddington limit and become observable as ULXs. We show that a Be-HMXB may become a ULX if the magnetic field of the neutron star and the density of the Be disc meet certain conditions. We also show that, although a stronger magnetic field increases the brightness of a neutron star ULX with a Be donor, its brightness cannot exceed the Eddington limit by a more than a factor of ≈50. Finally, we propose a scenario whereby some normal Be-HMXBs may evolve into ULXs as the donor evolves into a giant.

Beschreibung: We present an analysis of the intracluster light (ICL) in the Frontier Field Cluster MACS J1149.5+2223 (z = 0.544), which combines new and archival Hubble WFC3/IR imaging to provide continuous radial coverage out to 2.8 Mpc from the brightest cluster galaxy (BCG). Employing careful treatment of potential systematic biases and using data at the largest radii to determine the background sky level, we reconstruct the surface brightness profile out to a radius of 2 Mpc. This radius is the largest to which the ICL has been measured for an individual cluster. Within this radius, we measure a total luminosity of 1.5 × 1013 L⊙ for the BCG plus ICL. From the profile and its logarithmic slope, we identify the transition from the BCG to ICL at r ∼ 70 kpc. Remarkably, we also detect an apparent inflection in the profile centred in the 1.2–1.7 Mpc (0.37–0.52 r200m) radial bin, a signature of an infall caustic in the stellar distribution. Based upon the shape and strength of the feature, we interpret it as potentially being at the splashback radius, although the radius is smaller than theoretical predictions. If this is the splashback radius, then it is the first such detection in the ICL and the first detection of the splashback radius for an individual cluster. Similar analyses should be possible with the other Frontier Field clusters, and eventually with clusters observed by the Euclid and Roman missions.

Beschreibung: On the continental shelf, New Zealand bryozoans dominate ecologically-important three-dimensional benthic habitat providing structural complexity which hosts a wide variety of fauna, including economically valuable species (e.g. oysters and blue cod). The association between these species and bryozoan-dominated biogenic habitat commonly results in trawling damage to the benthos; eliminating pressure from destructive fishing practices could support bryozoan regrowth. In 2002, a voluntary fishing restriction was designated over part of one such bryozoan-dominated biogenic habitat, but the efficacy of this ban has not been assessed. Statutory marine protected areas have been proposed in the same area; they aim to achieve bryozoan protection and recovery. A comparison of the bryozoan fauna using benthic images from 2003 and 2019 from the Otago shelf, South Island, New Zealand/Te Waipounamu, Aotearoa shows that a 17-year voluntary fishing ban has been ineffective at increasing the cover and richness of bryozoan-dominated biogenic habitat on the Otago shelf. This study provides a baseline for future monitoring and develops useful metrics for recovery, while calling for further research on the biology and ecology of habitat-forming bryozoans.

Beschreibung: Ecosystem-based management is generally viewed as one of the most promising avenues for addressing the various anthropogenic pressures facing the world’s marine ecosystems. These approaches have been developed to varying degrees by individual countries or international organisations, but there remain a large proportion of marine ecosystems, particularly in developing regions, that have not yet been the subject of such research. In these areas, lack of effective regulation and the often high importance of the marine environment in providing food and economic opportunities, together create conditions where marine resources and habitats come under unsustainable levels of pressure. Here, we present a data-limited assessment approach to discern marine ecological patterns, in this case for the exclusive economic zone of Vietnam. By combining data from environmental and biological surveys from the Vietnamese national survey dataset and local oceanographic models, we have identified a series of 12 candidate ecological production units, delineated by their environmental characteristics, and the key commercial species that exist within them. These units are suggested as a possible foundation for a spatial management structure in the Vietnamese exclusive economic zone including considerations such as placement of marine protected areas, or ecological boundaries of key areas of socio-economic importance.

Beschreibung: This paper explores the possibility of using the ensemble modelling paradigm to fully capture assessment uncertainty and improve the robustness of advice provision. We identify and discuss advantages and challenges of ensemble modelling approaches in the context of scientific advice. There are uncertainties associated with every phase in the stock assessment process: data collection, assessment model choice, model assumptions, interpretation of risk, up to the implementation of management advice. Additionally, the dynamics of fish populations are complex, and our incomplete understanding of those dynamics and limited observations of important mechanisms, necessitate that models are simpler than nature. The aim is for the model to capture enough of the dynamics to accurately estimate trends and abundance, and provide the basis for robust advice about sustainable harvests. The status quo approach to assessment modelling has been to identify the “best” model and generate advice from that model, mostly ignoring advice from other model configurations regardless of how closely they performed relative to the chosen model. We discuss and make suggestions about the utility of ensemble models, including revisions to the formal process of providing advice to management bodies, and recommend further research to evaluate potential gains in modelling and advice performance.

Beschreibung: Motivation Co-evolution analysis can be used to accurately predict residue–residue contacts from multiple sequence alignments. The introduction of machine-learning techniques has enabled substantial improvements in precision and a shift from predicting binary contacts to predict distances between pairs of residues. These developments have significantly improved the accuracy of de novo prediction of static protein structures. With AlphaFold2 lifting the accuracy of some predicted protein models close to experimental levels, structure prediction research will move on to other challenges. One of those areas is the prediction of more than one conformation of a protein. Here, we examine the potential of residue–residue distance predictions to be informative of protein flexibility rather than simply static structure. Results We used DMPfold to predict distance distributions for every residue pair in a set of proteins that showed both rigid and flexible behaviour. Residue pairs that were in contact in at least one reference structure were classified as rigid, flexible or neither. The predicted distance distribution of each residue pair was analysed for local maxima of probability indicating the most likely distance or distances between a pair of residues. We found that rigid residue pairs tended to have only a single local maximum in their predicted distance distributions while flexible residue pairs more often had multiple local maxima. These results suggest that the shape of predicted distance distributions contains information on the rigidity or flexibility of a protein and its constituent residues. Supplementary information Supplementary data are available at Bioinformatics online.

Beschreibung: Climate effects on marine fish depend on life stage, particularly when life stages differ in habitat utilization. In the present study, we investigated life stage-dependent responses of lesser sandeel (Ammodytes marinus) to temperature at contrasting geographical scales. We related population density and individual growth to temperature and found different temperature responses between the first and the second years of life. During the first year of life, fish size was the single most important factor influencing sandeel abundances, indicating a positive relationship between growth and survival. In contrast, during the second year of life, autumn bottom temperature was negatively correlated with sandeel abundance, suggesting elevated mortality in warm years. Southerly areas, experiencing higher temperatures in general, were also the areas showing the strongest response to temperature. This study sheds light on how warming impacts population dynamics of one of the most important forage fishes in the North Sea and supports the discussion of underlying mechanisms.

Beschreibung: Motivation The well-known fact that protein structures are more conserved than their sequences forms the basis of several areas of computational structural biology. Methods based on the structure analysis provide more complete information on residue conservation in evolutionary processes. This is crucial for the determination of evolutionary relationships between proteins and for the identification of recurrent structural patterns present in biomolecules involved in similar functions. However, algorithmic structural alignment is much more difficult than multiple sequence alignment. This study is devoted to the development and applications of DAMA—a novel effective environment capable to compute and analyze multiple structure alignments. Results DAMA is based on local structural similarities, using local 3D structure descriptors and thus accounts for nearest-neighbor molecular environments of aligned residues. It is constrained neither by protein topology nor by its global structure. DAMA is an extension of our previous study (DEDAL) which demonstrated the applicability of local descriptors to pairwise alignment problems. Since the multiple alignment problem is NP-complete, an effective heuristic approach has been developed without imposing any artificial constraints. The alignment algorithm searches for the largest, consistent ensemble of similar descriptors. The new method is capable to capture most of the biologically significant similarities present in canonical test sets and is discriminatory enough to prevent the emergence of larger, but meaningless, solutions. Tests performed on the test sets, including protein kinases, demonstrate DAMA’s capability of identifying equivalent residues, which should be very useful in discovering the biological nature of proteins similarity. Performance profiles show the advantage of DAMA over other methods, in particular when using a strict similarity measure QC, which is the ratio of correctly aligned columns, and when applying the methods to more difficult cases. Availability and implementation DAMA is available online at http://dworkowa.imdik.pan.pl/EP/DAMA. Linux binaries of the software are available upon request. Supplementary information Supplementary data are available at Bioinformatics online.

Beschreibung: Hi intensity mapping will provide maps of the large-scale distribution of neutral hydrogen (Hi) in the universe. These are prime candidates to be used to constrain primordial non-Gaussianity using the Large Scale Structure of the Universe as well as to provide further tests of Einstein’s theory of Gravity (GR). But Hi maps are contaminated by foregrounds, which can be several orders of magnitude above the cosmological signal. Here we quantify how degenerated are the large-scale effects (fNL and GR effects) with the residual foregrounds. We conclude that a joint analysis does not provide a catastrophic degradation of constraints and provides a framework to determine the marginal errors of large scale-effects in the presence of foregrounds. Similarly, we conclude that the macroscopical properties of the foregrounds can be measured with high precision. Notwithstanding, such results are highly dependent on accurate forward modelling of the foregrounds, which incorrectly done catastrophically bias the best fit values of cosmological parameters, foreground parameterisations, and large-scale effects.

Beschreibung: Motivation Assigning new sequences to known protein families and subfamilies is a prerequisite for many functional, comparative and evolutionary genomics analyses. Such assignment is commonly achieved by looking for the closest sequence in a reference database, using a method such as BLAST. However, ignoring the gene phylogeny can be misleading because a query sequence does not necessarily belong to the same subfamily as its closest sequence. For example, a hemoglobin which branched out prior to the hemoglobin alpha/beta duplication could be closest to a hemoglobin alpha or beta sequence, whereas it is neither. To overcome this problem, phylogeny-driven tools have emerged but rely on gene trees, whose inference is computationally expensive. Results Here, we first show that in multiple animal and plant datasets, 18 to 62% of assignments by closest sequence are misassigned, typically to an over-specific subfamily. Then, we introduce OMAmer, a novel alignment-free protein subfamily assignment method, which limits over-specific subfamily assignments and is suited to phylogenomic databases with thousands of genomes. OMAmer is based on an innovative method using evolutionarily-informed k-mers for alignment-free mapping to ancestral protein subfamilies. Whilst able to reject non-homologous family-level assignments, we show that OMAmer provides better and quicker subfamily-level assignments than approaches relying on the closest sequence, whether inferred exactly by Smith-Waterman or by the fast heuristic DIAMOND. Availability OMAmer is available from the Python Package Index (as omamer), with the source code and a precomputed database available at https://github.com/DessimozLab/omamer. Supplementary information Supplementary data are available at Bioinformatics online.

Beschreibung: The [O iii] 88 μm fine structure emission line has been detected into the Epoch of Reionization (EoR) from star-forming galaxies at redshifts 6 〈 z ≲ 9 with ALMA. These measurements provide valuable information regarding the properties of the interstellar medium (ISM) in the highest redshift galaxies discovered thus far. The [O iii] 88 μm line observations leave, however, a degeneracy between the gas density and metallicity in these systems. Here we quantify the prospects for breaking this degeneracy using future ALMA observations of the [O iii] 52 μm line. Among the current set of ten [O iii] 88 μm emitters at 6 〈 z ≲ 9, we forecast 52 μm detections (at 6-σ) in SXDF-NB1006-2, B14-6566, J0217-0208, and J1211-0118 within on-source observing times of 2-10 hours, provided their gas densities are larger than about nH ≳ 102 − 103 cm−3. Other targets generally require much longer integration times for a 6-σ detection. Either successful detections of the 52 μm line, or reliable upper limits, will lead to significantly tighter constraints on ISM parameters. The forecasted improvements are as large as ∼3 dex in gas density and ∼1 dex in metallicity for some regions of parameter space. We suggest SXDF-NB1006-2 as a promising first target for 52 μm line measurements. We discuss how such measurements will help in understanding the mass metallicity relationship during the EoR.

Beschreibung: Misalignments within protoplanetary discs are now commonly observed, and features such as shadows in scattered light images indicate departure from a co-planar geometry. VLT/SPHERE observations of the disc around HD 143006 show a large-scale asymmetry, and two narrow dark lanes which are indicative of shadowing. ALMA observations also reveal the presence of rings and gaps in the disc, along with a bright arc at large radii. We present new hydrodynamic simulations of HD 143006, and show that a configuration with both a strongly inclined binary and an outer planetary companion is the most plausible to explain the observed morphological features. We compute synthetic observations from our simulations, and successfully reproduce both the narrow shadows and the brightness asymmetry seen in IR scattered light. Additionally, we reproduce the large dust observed in the mm continuum, due to a 10 Jupiter mass planet detected in the CO kinematics. Our simulations also show the formation of a circumplanetary disc, which is misaligned with respect to the outer disc. The narrow shadows cast by the inner disc and the planet-induced “kink” in the disc kinematics are both expected to move on a time-scale of ∼5–10 years, presenting a potentially observable test of our model. If confirmed, HD 143006 would be the first known example of a circumbinary planet on a strongly misaligned orbit.

Beschreibung: Short-period super-Earth-sized planets are common. Explaining how they form near their present orbits requires understanding the structure of the inner regions of protoplanetary discs. Previous studies have argued that the hot inner protoplanetary disc is unstable to the magneto-rotational instability (MRI) due to thermal ionization of potassium, and that a local gas pressure maximum forms at the outer edge of this MRI-active zone. Here we present a steady-state model for inner discs accreting viscously, primarily due to the MRI. The structure and MRI-viscosity of the inner disc are fully coupled in our model; moreover, we account for many processes omitted in previous such models, including disc heating by both accretion and stellar irradiation, vertical energy transport, realistic dust opacities, dust effects on disc ionization and non-thermal sources of ionization. For a disc around a solar-mass star with a standard gas accretion rate ($dot{M}$∼10−8 M⊙ yr−1) and small dust grains, we find that the inner disc is optically thick, and the accretion heat is primarily released near the midplane. As a result, both the disc midplane temperature and the location of the pressure maximum are only marginally affected by stellar irradiation, and the inner disc is also convectively unstable. As previously suggested, the inner disc is primarily ionized through thermionic and potassium ion emission from dust grains, which, at high temperatures, counteract adsorption of free charges onto grains. Our results show that the location of the pressure maximum is determined by the threshold temperature above which thermionic and ion emission become efficient.

Beschreibung: Summary VCF files with results of sequencing projects take a lot of space. We propose the VCFShark, which is able to compress VCF files up to an order of magnitude better than the de facto standards (gzipped VCF and BCF). The advantage over competitors is the greatest when compressing VCF files containing large amounts of genotype data. The processing speeds up to 100 MB/s and main memory requirements lower than 30 GB allow to use our tool at typical workstations even for large datasets. Availability and Implementation https://github.com/refresh-bio/vcfshark Supplementary information Supplementary data are available at publisher’s Web site.

Beschreibung: Local star-forming galaxies show properties that are thought to differ from galaxies in the early Universe. Among them, the ionizing stellar populations and the gas geometry make the recipes designed to derive chemical abundances from nebular emission lines to differ from those calibrated in the Local Universe. A sample of 1969 Extreme Emission Line Galaxies (EELGs) at a redshift 0 ≲ z ≲ 0.49, selected from the Sloan Digital Sky Survey (SDSS) to be local analogues of high-redshift galaxies, was used to analyse their most prominent emission lines and to derive total oxygen abundances and nitrogen-to-oxygen ratios following the direct method in the ranges 7.7 〈 12+log(O/H) 〈 8.6 and -1.8 〈 log(N/O) 〈 -0.8. They allow us to obtain new empirically calibrated strong-line methods and to evaluate other recipes based on photoionization models that can be later used for a chemical analysis of actively star-forming galaxies in very early stages of galaxy evolution. Our new relations are in agreement with others found for smaller samples of objects at higher redshifts. When compared with other relations calibrated in the local Universe, they differ when the employed strong-line ratio depends on the hardness of the ionizing radiation, such as O32 or Ne3O2, but they do not when the main dependence is on the ionization parameter, such as S23. In the case of strong-line ratios depending on [N ii] lines, the derivation of O/H becomes very uncertain due to the very high N/O values derived in this sample, above all in the low-metallicity regime. Finally, we adapt the bayesian-like code HII-Chi-mistry for the conditions found in this kind of galaxies and we prove that it can be used to derive within errors both O/H and N/O, in consistency with the direct method .

Beschreibung: Stellar tidal streams are sensitive tracers of the properties of the gravitational potential in which they orbit and detailed observations of their density structure can be used to place stringent constraints on fluctuations in the potential caused by, e.g., the expected populations of dark matter subhalos in the standard cold dark matter paradigm (CDM). Simulations of the evolution of stellar streams in live N-body halos without low-mass dark-matter subhalos, however, indicate that streams exhibit significant perturbations on small scales even in the absence of substructure. Here we demonstrate, using high-resolution N-body simulations combined with sophisticated semi-analytic and simple analytic models, that the mass resolutions of 104–$10^5, m {M}_{odot }$ commonly used to perform such simulations cause spurious stream density variations with a similar magnitude on large scales as those expected from a CDM-like subhalo population and an order of magnitude larger on small, yet observable, scales. We estimate that mass resolutions of $approx 100, m {M}_{odot }$ ($approx 1, m {M}_{odot }$) are necessary for spurious, numerical density variations to be well below the CDM subhalo expectation on large (small) scales. That streams are sensitive to a simulation’s particle mass down to such small masses indicates that streams are sensitive to dark matter clustering down to these low masses if a significant fraction of the dark matter is clustered or concentrated in this way, for example, in MACHO models with masses of 10–$100, m {M}_{odot }$.

Beschreibung: We present TES, a new n-body integration code for the accurate and rapid propagation of planetary systems in the presence of close encounters. TES builds upon the classic Encke method and integrates only the perturbations to Keplerian trajectories to reduce both the error and runtime of simulations. Variable step size is used throughout to enable close encounters to be precisely handled. A suite of numerical improvements are presented that together make TES optimal in terms of energy error. Lower runtimes are found in the majority of test problems considered when compared to direct integration using IAS15. TES is freely available.

Beschreibung: As deep spectroscopic campaigns extend to higher redshifts and lower stellar masses, the interpretation of galaxy spectra depends increasingly upon models for very young stellar populations. Here we present new HST/COS ultraviolet spectroscopy of seven nearby (

Beschreibung: Supernova Ia magnitude surveys measure the dimensionless luminosity distance H0DL. However, from the distances alone one cannot obtain quantities like H(z) or the dark energy equation of state, unless further cosmological assumptions are imposed. Here we show that by measuring the power spectrum of density contrast and of peculiar velocities of supernovae one can estimate also H(z)/H0 regardless of background or linearly perturbed cosmology and of galaxy-matter bias. This method, dubbed Clustering of Standard Candles (CSC) also yields the redshift distortion parameter β(k, z) and the biased matter power spectrum in a model-independent way. We forecast that an optimistic (pessimistic) LSST may be able to constrain H(z)/H0 to 5–13% (9–40%) in redshift bins of Δz = 0.1 up to at least z = 0.6.

Beschreibung: We identify 228 periodic variables in the region of young open cluster NGC 281 using time series photometry carried out from 1 m class ARIES telescopes, Nainital. The cluster membership of these identified variables is determined on the basis colour-magnitude, two colour diagrams and kinematic data. Eighty one variable stars are found to be members of the cluster NGC 281. Of 81 variables, 30 and 51 are probable main-sequence and pre-main-sequence members, respectively. Present study classifies main-sequence variable stars into different variability types according to their periods of variability, shape of light curves and location in the Hertzsprung-Russell diagram. These identified main-sequence variables could be β Cep, δ Scuti, slowly pulsating B type and new class variables. Among 51 pre-main-sequence variables, majority of them are weak line T Tauri stars. The remaining 147 variables could belong to the field population. The variability characteristics of the field population indicate that these variables could be RR Lyrae, δ Scuti and binaries type variables.

Beschreibung: ΛCDM cosmological models with Early Dark Energy (EDE) have been proposed to resolve tensions between the Hubble constant H0 = 100 h kmṡ−1Ṁpc−1 measured locally, giving h ≈ 0.73, and H0 deduced from Planck cosmic microwave background (CMB) and other early universe measurements plus ΛCDM, giving h ≈ 0.67. EDE models do this by adding a scalar field that temporarily adds dark energy equal to about 10% of the cosmological energy density at the end of the radiation-dominated era at redshift z ∼ 3500. Here we compare linear and nonlinear predictions of a Planck-normalized ΛCDM model including EDE giving h = 0.728 with those of standard Planck-normalized ΛCDM with h = 0.678. We find that nonlinear evolution reduces the differences between power spectra of fluctuations at low redshifts. As a result, at z = 0 the halo mass functions on galactic scales are nearly the same, with differences only 1-2%. However, the differences dramatically increase at high redshifts. The EDE model predicts 50% more massive clusters at z = 1 and twice more galaxy-mass haloes at z = 4. Even greater increases in abundances of galaxy-mass haloes at higher redshifts may make it easier to reionize the universe with EDE. Predicted galaxy abundances and clustering will soon be tested by the James Webb Space Telescope (JWST) observations. Positions of baryonic acoustic oscillations (BAOs) and correlation functions differ by about 2% between the models – an effect that is not washed out by nonlinearities. Both standard ΛCDM and the EDE model studied here agree well with presently available acoustic-scale observations, but the Dark Energy Spectroscopic Instrument (DESI) and Euclid measurements will provide stringent new tests.

Beschreibung: Motivation As the generation of complex single-cell RNA sequencing datasets becomes more commonplace it is the responsibility of researchers to provide access to these data in a way that can be easily explored and shared. Whilst it is often the case that data is deposited for future bioinformatic analysis many studies do not release their data in a way that is easy to explore by non-computational researchers. Results In order to help address this we have developed ShinyCell, an R package that converts single-cell RNA sequencing datasets into explorable and shareable interactive interfaces. These interfaces can be easily customised in order to maximise their usability and can be easily uploaded to online platforms to facilitate wider access to published data. Availability ShinyCell is available at https://github.com/SGDDNB/ShinyCell and https://figshare.com/projects/ShinyCell/100439.

Beschreibung: Motivation Genomic selection (GS) is currently deemed the most effective approach to speed up breeding of agricultural varieties. It has been recognized that consideration of multiple traits in GS can improve accuracy of prediction for traits of low heritability. However, since GS forgoes statistical testing with the idea of improving predictions, it does not facilitate mechanistic understanding of the contribution of particular single nucleotide polymorphisms (SNP). Results Here we propose a L2,1-norm regularized multivariate regression model and devise a fast and efficient iterative optimization algorithm, called L2,1-joint, applicable in multi-trait GS. The usage of the L2,1-norm facilitates variable selection in a penalized multivariate regression that considers the relation between individuals, when the number of SNPs is much larger than the number of individuals. The capacity for variable selection allows us to define master regulators that can be used in a multi-trait GS setting to dissect the genetic architecture of the analyzed traits. Our comparative analyses demonstrate that the proposed model is a favorable candidate compared to existing state-of-the-art approaches. Prediction and variable selection with data sets from Brassica napus, wheat and Arabidopsis thaliana diversity panels are conducted to further showcase the performance of the proposed model. Availability and implementation The model is implemented using R programming language and the code is freely available from https://github.com/alainmbebi/L21-norm-GS. Supplementary information Supplementary data are available at Bioinformatics online.

Beschreibung: Summary Finding informative predictive features in high dimensional biological case-control datasets is challenging. The Extreme Pseudo-Sampling (EPS) algorithm offers a solution to the challenge of feature selection via a combination of deep learning and linear regression models. First, using a variational autoencoder, it generates complex latent representations for the samples. Second, it classifies the latent representations of cases and controls via logistic regression. Third, it generates new samples (pseudo-samples) around the extreme cases and controls in the regression model. Finally, it trains a new regression model over the upsampled space. The most significant variables in this regression are selected. We present an open-source implementation of the algorithm that is easy to set up, use, and customize. Our package enhances the original algorithm by providing new features and customizability for data preparation, model training and classification functionalities. We believe the new features will enable the adoption of the algorithm for a diverse range of datasets. Availability The software package for Python is available online at https://github.com/roohy/eps

Beschreibung: The detection of the binary neutron star (BNS) merger GW170817 and the associated electromagnetic (EM) counterpart, the “kilonova” (kN) AT2017gfo, opened a new era in multi-messenger astronomy. However, despite many efforts, it has been proven very difficult to find additional kilonovae, even though LIGO/Virgo has reported at least one BNS event during their latest run, O3. The focus of this work is the exploration of the sensitivity of the adopted optical surveys searching for kNe during O3. We propose ways to optimize the choices of filters and survey depth to boost the detection efficiency for these faint and fast-evolving transients in the future. In particular, we use kN models to explore the dependence on ejecta mass, geometry, viewing angle, wavelength coverage and source distance. We find that the kN detection efficiency has a strong viewing-angle dependency, especially for filters blueward of i-band. This loss of sensitivity can be mitigated by early, deep, observations. Efficient gri counterpart searches for kNe at ∼200 Mpc would require reaching a limiting magnitude mlim = 23 mag, to ensure good sensitivity over a wide range of the model phase-space. We conclude that kN searches during O3 were generally too shallow to detect BNS optical counterparts, even under optimistic assumptions.

Beschreibung: Theoretical models are vital for exploring the galaxy merger process, which plays a crucial role in the evolution of galaxies. Recent advances in modelling have placed tight constraints on the buildup of stellar material in galaxies across cosmic time. Despite these successes, extracting the merger rates from observable data remains a challenge. Differences in modelling techniques, combined with limited observational data, drive conflicting conclusions on the merging timescales of close pairs. We employ an empirical model for galaxy formation that links galaxy properties to the growth of simulated dark matter haloes, along with mock lightcone galaxy catalogues, to probe the dependencies of pair merging probabilities and merging timescales. In this work, we demonstrate that the pair merging probabilities are best described by a logistic function and that mean merging timescales can be well approximated by a linear relation in the projected separation and line of sight velocity difference in observed pairs. Together, our fitting formulae can accurately predict merger rates from galaxy pairs to at least z ∼ 4 under a wide variety of pair selection criteria. Additionally, we show that some commonly used pair selection criteria may not represent a suitable sample of galaxies to reproduce underlying merger rates. Finally, we conclude from our analysis that observation timescales are primarily driven by dynamics and are not strongly impacted by the star formation properties of the component galaxies.

Beschreibung: Sue Bowler remembers the Guernsey-born astronomer

Beschreibung: Jacqueline Mitton examines one strand in the life and research of the pioneering astronomer

Beschreibung: Robert Massey ponders the potential post-pandemic world for those working in astronomy and geophysics, following an online meeting in January

Beschreibung: A digest of some of the big stories from the past two months