Search

We study two fringe subtree counting statistics, the number of cherries and that of pitchforks for Ford's $\alpha$ model, a one-parameter family of random phylogenetic tree models that includes the uniform and the Yule models, two tree models commonly used in phylogenetics. Based on a nonuniform version of the extended P\'olya urn models in which negative entries are permitted for their replacement matrices, we obtain the strong law of large numbers and the central limit theorem for the joint distribution of these two count statistics for the Ford model. Furthermore, we derive a recursive formula for computing the exact joint distribution of these two statistics. This leads to exact formulas for their means and higher order asymptotic expansions of their second moments, which allows us to identify a critical parameter value for the correlation between these two statistics. That is, when $n$ is sufficiently large, they are negatively correlated for $0\le \alpha \le 1/2$ and positively correlated for $1/2, Comment: 23 pages, 2 figures

RNA molecules are single-stranded analogues of DNA that can fold into various structures which influence their biological function within the cell. RNA structures can be modelled combinatorially in terms of a certain type of graph called an RNA diagram. In this paper we introduce a new poset of RNA diagrams $\mathcal{B}^r_{f,k}$, $r\ge 0$, $k \ge 1$ and $f \ge 3$, which we call the Penner-Waterman poset, and, using results from the theory of multitriangulations, we show that this is a pure poset of rank $k(2f-2k+1)+r-f-1$, whose geometric realization is the join of a simplicial sphere of dimension $k(f-2k)-1$ and an $\left((f+1)(k-1)-1\right)$-simplex in case $r=0$. As a corollary for the special case $k=1$, we obtain a result due to Penner and Waterman concerning the topology of the space of RNA secondary structures. These results could eventually lead to new ways to investigate landscapes of RNA $k$-noncrossing structures.

SUMMARYCalcium signalling is central to many plant processes, with families of calcium decoder proteins having expanded across the green lineage and redundancy existing between decoders. The liverwort Marchantia polymorpha has fast become a new model plant, but it is unclear what calcium decoders exist in this species.We have performed phylogenetic analyses to identify the Calcineurin B-Like (CBL) and CBL-Interacting Protein Kinase (CIPK) network of M. polymorpha. We analysed CBL-CIPK expression during salt stress, and determined protein-protein interactions using yeast two-hybrid and bimolecular fluorescence complementation. We also created genetic knockouts using CRISPR/Cas9.We confirm that M. polymorpha has two CIPKs and three CBLs. Both CIPKs and only one CBL show salt-responsive transcriptional changes. All M. polymorpha CBL-CIPKs interact with each other in planta. Knocking out CIPK-B causes increased sensitivity to salt suggesting that this CIPK is involved in salt signalling.We have identified CBL-CIPKs that form part of a salt tolerance pathway in M. polymorpha. Phylogeny and interaction studies imply that these CBL-CIPKs form an evolutionarily conserved Salt Overly Sensitive (SOS) pathway. Hence, salt responses may be some of the early functions of CBL-CIPK networks and increased abiotic stress tolerance required for land plant emergence.

Recent advances in single-cell RNA sequencing (scRNA-seq) provide exciting opportunities for transcriptome analysis at single-cell resolution. Clustering individual cells is a key step to reveal cell subtypes and infer cell lineage in scRNA-seq analysis. Although many dedicated algorithms have been proposed, clustering quality remains a computational challenge for scRNA-seq data, which is exacerbated by inflated zero counts due to various technical noise. To address this challenge, we assess the combinations of nine popular dropout imputation methods and eight clustering methods on a collection of 10 well-annotated scRNA-seq datasets with different sample sizes. Our results show that (i) imputation algorithms do typically improve the performance of clustering methods, and the quality of data visualization using t-Distributed Stochastic Neighbor Embedding; and (ii) the performance of a particular combination of imputation and clustering methods varies with dataset size. For example, the combination of single-cell analysis via expression recovery and Sparse Subspace Clustering (SSC) methods usually works well on smaller datasets, while the combination of adaptively-thresholded low-rank approximation and single-cell interpretation via multikernel learning (SIMLR) usually achieves the best performance on larger datasets.

A rooted phylogenetic network is a directed acyclic graph with a single root, whose sinks correspond to a set of species. As such networks are useful for representing the evolution of species that have undergone reticulate evolution, there has been great interest in developing the theory behind and algorithms for constructing them. However, unlike evolutionary trees, these networks can be highly non-planar, which can make them difficult to visualise and interpret. Here we investigate properties of planar rooted phylogenetic networks and algorithms for deciding whether or not rooted networks have certain special planarity properties. In particular, we introduce three natural subclasses of planar rooted phylogenetic networks and show that they form a hierarchy. In addition, for the well-known level-k networks, we show that level-1, -2, -3 networks are always outer, terminal, and upward planar, respectively, and that level-4 networks are not necessarily planar. Finally, we show that a regular network is terminal planar if and only if it is pyramidal. Our results make use of the highly developed field of planar digraphs, and we believe that the link between phylogenetic networks and planar graphs should prove useful in future for developing new approaches to both construct and visualise phylogenetic networks., 11 pages

In order to reduce the pollutant emission and alleviate the pressure of petroleum resources shortage and greenhouse gas emission at the same time, the use of clean and renewable alternative fuel for marine engines is a promising option. In this study, a marine diesel engine, which was modified to run in diesel methanol compound combustion (DMCC) mode, was investigated. After the diesel injection parameters were calibrated, and combined with a sample after-treatment device DOC (diesel oxidation catalyst), the engine could meet the requirements of China II legislation. The overall MSP (methanol substitute percent) reached 54.1%. The value of each pollutant emission was much lower than that in China II emission legislation, and there was almost no methanol and formaldehyde emissions. When methanol was injected into the inlet manifold, the intake air temperature decreased a lot, as well as the exhaust gas temperature, which were beneficial to increase engine thermal efficiency and improve engine room environment. Compared with the engine running in pure diesel mode, when the engine ran in diesel/methanol dual fuel mode, the combustion phase was advanced, and the combustion duration became shorter. Therefore, the engine thermal efficiency increased, and fuel consumption decreased significantly.

We developed a method, ArchaicSeeker 2.0, to identify introgressed hominin sequences and model multiple-wave admixture. The new method enabled us to discern two waves of introgression from both Denisovan-like and Neanderthal-like hominins in present-day Eurasian populations and an ancient Siberian individual. We estimated that an early Denisovan-like introgression occurred in Eurasia around 118.8–94.0 thousand years ago (kya). In contrast, we detected only one single episode of Denisovan-like admixture in indigenous peoples eastern to the Wallace-Line. Modeling ancient admixtures suggested an early dispersal of modern humans throughout Asia before the Toba volcanic super-eruption 74 kya, predating the initial peopling of Asia as proposed by the traditional Out-of-Africa model. Survived archaic sequences are involved in various phenotypes including immune and body mass (e.g., ZNF169), cardiovascular and lung function (e.g., HHAT), UV response and carbohydrate metabolism (e.g., HYAL1/HYAL2/HYAL3), while “archaic deserts” are enriched with genes associated with skin development and keratinization., Existing methods to identify the presence of DNA from other hominin species can be limited in the ability to accurately estimate introgression waves, or can only be applied to specific populations. Here, the authors have developed a generalizable method to identify introgression in multi-wave situations.

Phylogenetic networks are a generalization of evolutionary or phylogenetic trees that are commonly used to represent the evolution of species which cross with one another. A special type of phylogenetic network is an {\em $X$-cactus}, which is essentially a cactus graph in which all vertices with degree less than three are labelled by at least one element from a set $X$ of species. In this paper, we present a way to {\em encode} $X$-cactuses in terms of certain collections of partitions of $X$ that naturally arise from $X$-cactuses. Using this encoding, we also introduce a partial order on the set of $X$-cactuses (up to isomorphism), and derive some structural properties of the resulting partially ordered set. This includes an analysis of some properties of its least upper and greatest lower bounds. Our results not only extend some fundamental properties of phylogenetic trees to $X$-cactuses, but also provides a new approach to solving topical problems in phylogenetic network theory such as deriving consensus networks.

Particulate matter (PM) emissions have become an increasingly noticeable concern for gasoline vehicles, especially for gasoline direct injection (GDI) one. The PM index (PMI) model developed by Aikiwa et al. provides an effective model linking the fuel properties to the PM emissions from gasoline vehicles. In this study, a more practical reduced PM index (RPMI) was developed based on the detailed analysis of 22 fuels with various physical properties and chemical compositions. Two approaches of statistical mathematics of correlation analysis and multivariate linear regression (MLR) were adopted in the modeling process. The RPMI involving two global fuel properties of T90 (distillation temperature of 90% by volume) and T70 (distillation temperature of 70% by volume) was tested to be statistically valid. In addition, the model was verified through the engine bench tests and the vehicle emissions tests. The results reveal that the RPMI showed good correlations to the engine-out particle number (PN) emissions under all of the typical test conditions, including PFI (port fuel injection) mode, GDI (gasoline direct injection) mode and compound injection (PFI + GDI) mode. In addition, the RPMI demonstrated a significantly high correlation to vehicle-out PN emissions over the New European Driving Cycle (NEDC) with the determination coefficient R2 = 0.947. Moreover, a moderate correlation (R2 = 0.801) of the index to the filter PM mass emissions was observed. In the light of no additional components analysis is needed except for the legitimate tests of fuels, the RPMI has the potential to be a useful tool for original equipment manufacturers (OEMs) and environmental certification departments to expediently evaluate the fuel effect on PM emissions.

Phylogenetic trees and networks are leaf-labelled graphs used to model evolution. Display graphs are created by identifying common leaf labels in two or more phylogenetic trees or networks. The treewidth of such graphs is bounded as a function of many common dissimilarity measures between phylogenetic trees and this has been leveraged in fixed param- eter tractability results. Here we further elucidate the properties of dis- play graphs and their interaction with treewidth. We show that it is NP-hard to recognize display graphs, but that display graphs of bounded treewidth can be recognized in linear time. Next we show that if a phyloge- netic network displays (i.e. topologically embeds) a phylogenetic tree, the treewidth of their display graph is bounded by a function of the treewidth of the original network (and also by various other parameters). In fact, us- ing a bramble argument we show that this treewidth bound is sharp up to an additive term of 1. We leverage this bound to give an FPT algorithm, parameterized by treewidth, for determining whether a network displays a tree, which is an intensively-studied problem in the field. We conclude with a discussion on the future use of display graphs and treewidth in phylogenetics.

The late-injection strategy (close after main-injection) of common-rail diesel engine is capable of enhancing combustion turbulence and reducing particulate matter (PM) emissions. In this work, experimental study was performed to investigate the effect of diesel late-injection on combustion and emissions characteristics of diesel/methanol dual fuel (DMDF) engine. The experiments were carried out at a constant engine speed of 1340 rpm and a medium load of 1.0 MPa brake mean effective pressure (BMEP) with various methanol substitution ratio (MSR) on a common rail DMDF engine. The results reveal that higher MSR caused simultaneous decrease of nitric oxides (NOx) and accumulation mode PM emissions in spite of late-injection strategy. In particular, an augment of up to 12.8% in nucleation mode particle number (PN) was observed as MSR increased from 15% to 50%. A small late-injection of diesel ( 1200 μs) induced a simultaneous increase of particles with different size.

Methanol is difficult to be ignited by compression due to its low cetane number. The strategy of diesel–methanol compound combustion (DMCC) was proposed, which starts up with pure diesel fuel and switches to the mode of diesel–methanol dual fuel after the engine fully warms up. The strategy of DMCC will substitute methanol to diesel fuel up to 40% for new engine and 30% for retrofitting engine when it is applied to heavy-duty vehicle and constructive machine as well as marine power unit. The engines equipped with DMCC system are able to reduce NOx without assistance of selective catalyst reduction (SCR) urea system, as well as PM emission at the same time, which can meet the requirements of China V emission legislation and have the potential to meet the demand of China VI in future. The chapter introduces the strategy of DMCC, the main components dealing with methanol while applied to those engines mentioned above, and real examples on different power units. Additionally, the characteristics of diesel–methanol dual fuel (DMDF) combustion mode, emission control strategy, and system related to catalyst, retrofitting, and modifying engine from pure diesel to running DMCC mode are also illustrated. Again, it provides with the research results about the characteristics of DMDF combustion mode to assist the readers to understand the dual fuel combustion in depth as well as the examples applied to various power units, such as heavy-duty trucks, constructive machine, marine propeller, and locomotive engines.

SUMMARYLymphocyte biology research commonly involves purification of lymphocyte subpopulations by fluorescence-activated cell sorting (FACS) or immunomagnetic separation (IMS), both of which typically rely on antibody labeling of validated cell markers. Methods enabling label-free segregation of lymphocyte subpopulations would be invaluable with regard to less-perturbation, simplicity and cost-effectiveness. Here, we introduce TRuST, a label-free approach forTcellreconstitutionusingside-scatter (SSC). TRuST-sorted SSClowcells enrich for CD4+T and naïve T cells, while SSChighcells enrich for CD8+T, NK and differentiated T cells. Enrichment purity can be improved by computational gate design. SSClowcells have superior expansion capacity and generate more central memory precursors with naïve-resembling cytokine responses. Moreover, we find that both T cell differentiation status and CD4/CD8 T ratio in the starting cellular material are critical attributes predicting T cell product quality and quantity. TRuST presents an effective and reliable technique for label-free lymphocytes selection and reconstitution.

Rearrangements are discrete processes whereby discrete segments of DNA are deleted, replicated and inserted into novel positions. A sequence of such configurations, termed a rearrangement evolution, results in jumbled DNA arrangements, frequently observed in cancer genomes. We introduce a method that allows us to precisely count these different evolutions for a range of processes including breakage-fusion-bridge-cycles, tandem-duplications, inverted-duplications, reversals, transpositions and deletions, showing that the space of rearrangement evolution is super-exponential in size. These counts assume the infinite sites model of unique breakpoint usage.

Tree shape statistics are important for investigating evolutionary mechanisms mediating phylogenetic trees. As a step towards bridging shape statistics between rooted and unrooted trees, we present a comparison study on two subtree statistics known as numbers of cherries and pitchforks for the proportional to distinguishable arrangements (PDA) and the Yule-Harding-Kingman (YHK) models. Based on recursive formulas on the joint distribution of the number of cherries and that of pitchforks, it is shown that cherry distributions are log-concave for both rooted and unrooted trees under these two models. Furthermore, the mean number of cherries and that of pitchforks for unrooted trees converge respectively to those for rooted trees under the YHK model while there exists a limiting gap of 1/4 for the PDA model. Finally, the total variation distances between the cherry distributions of rooted and those of unrooted trees converge for both models. Our results indicate that caution is required for conducting statistical analysis for tree shapes involving both rooted and unrooted trees., Comment: 26 pages

In this paper, an experimental investigation on achieving ultra-low raw emissions in a diesel/methanol dual fuel engine based on dual-loop EGR was carried out. The effect of dual-loop EGR (High pressure EGR and low pressure EGR) on the combustion, performance and emissions of methanol engine has been studied comprehensively. The results show that ultra-low NOx (

UPGMA (Unweighted Pair Group Method with Arithmetic Mean) is a widely used clustering method. Here we show that UPGMA is a greedy heuristic for the normalized equidistant minimum evolution (NEME) problem, that is, finding a rooted tree that minimizes the minimum evolution score relative to the dissimilarity matrix among all rooted trees with the same leaf-set in which all leaves have the same distance to the root. We prove that the NEME problem is NP-hard. In addition, we present some heuristic and approximation algorithms for solving the NEME problem, including a polynomial time algorithm that yields a binary, rooted tree whose NEME score is within O(log^2 n) of the optimum. We expect that these results to eventually provide further insights into the behavior of the UPGMA algorithm., Comment: 29 pages, 8 figures

This paper explored the reason that the amount of methanol consumed is much lower than that of its theoretical value in terms of calorific value in diesel methanol dual fuel (DMDF) combustion by the means of engine bench test and modeling analysis. The results from experiments show that, the DMDF has much higher combustion rate than that of the correspondent diesel (D) mode and the accelerating effect caused by methanol changes with engine load, methanol-air mixture concentration and intake air temperature. At low load condition, the ignition delay caused by methanol is dominated while the accelerating effect is weak. However, this situation just turned upside down at high loads, to bring forth great improvement in isochoric degree. When methanol conducts quickly burning, both the combustion efficiency and the conversion efficiency from heat to work are improved, and the replacement ratio S R becomes much lower than its theoretical value 2.16. In order to reveal the mechanism of high efficiency DMDF at various running conditions, a 3-dimensional CFD model to simulate and analysis the process of methanol burning together with diesel was built. By fixing constant operating parameters and boundary conditions, both the putting off and accelerating effects caused by methanol are enhanced with increasing methanol concentration. The putting off effect weakens with rising air temperature, while the accelerating effect becomes more effective in this process. Hence, the isochoric degree for DMDF deteriorates with rising air temperature at lean methanol mixture but improves as the concentration increases. Finally, a relational graph among isochoric degree, heat release concentration degree, methanol-air mixture concentration and temperature are derived, from which we get the understanding of controllable methanol quick combustion (CMQC) to further improve the DMDF thermal efficiency.

Based on European Steady-state Cycle (ESC) test procedure, an engine test was carried out to meet the demand of Euro V legislation with diesel methanol compound combustion (DMCC) technology on a common-rail, turbocharged, inter-cooling heavy-duty and high pressure exhaust gas recirculation (HP EGR) diesel engine with meeting the requirements of Euro IV legislation. The engine was just modified with DMCC technology and combined with after-treatment device diesel oxidation catalyst coupled with particulate oxidation catalyst (DPOC). The test results show that it can fully meet the demand of Euro V legislation without the help of selected catalyst reduction (SCR). Additionally, the comparative analysis of the dynamic and economic performances between the DMCC Euro V engine and the baseline Euro IV engine was presented in this study. The regulated emissions (included HC, CO, NOX, PM (particulate matter)) and unregulated emissions (included formaldehyde (HCHO) and unburned methanol (CH3OH)) before and after DPOC from DMCC Euro V engine were also investigated. The dynamic performance of DMCC Euro V engine is basically consistent with the baseline Euro IV engine, while there is a slight deterioration in economic performance at certain operation conditions in terms of calorific value. Compared with the baseline engine before DPOC, the DMCC engine causes an increase in HC, CO, HCHO and CH3OH emissions, but the NOX and PM emissions are much better. The DPOC efficiency in the conversion rate of HC, CO, HCHO and CH3OH is over 98% under the whole ESC test procedure, and for PM, the conversion efficiency is only around 21%. However, DPOC has no obvious influence on the conversion of NOX emissions.

The increasing of fuel temperature can reduce the droplet size and have an advantage of improving spray atomization, while investigations of the effect of temperature on the methanol injector internal flow and external spray is rare. Firstly, a detailed three dimensional numerical simulations of nozzle internal flow have been conducted to probe into the cavitation in methanol injector nozzles, and then an experimental study has been carried out to investigate the droplet size and velocity of methanol spray at various temperatures using the Phase Doppler Particle Analyzer (PDPA) detecting system. And results show that the region of cavitations in nozzle orifice enlarges as methanol temperature and injection pressure increases, and the temperature for 'super-cavitation' occurring decreases gradually with the increasing of injection pressure. Moreover, the nozzle exit velocity, discharge coefficient and cavitations number were also analyzed. However, the discharge coefficient reduces nearly equal under various pressure when the methanol temperature is higher than 60 °C. In addition, the Sauter Mean Diameter (SMD) and velocity of methanol droplet were also analyzed, and found that the droplet velocity reaches the maximum value when the methanol temperature is 60 °C.

A diesel engine with 6-cylinder, turbocharged intercooling and common rail system was modified to operate on diesel/methanol dual fuel (DMDF). To optimize combustion process at high methanol substitution ratio (MSR) and further decrease exhaust emissions of DMDF engine at low load condition, pilot injection strategy was applied on it and the combustion, regulated and unregulated gaseous emission (included NOx, CO, THC, carbon dioxide (CO 2 ), nitrous oxide (N 2 O), unburned methanol (CH 3 OH), formaldehyde (HCHO), formic acid (HCOOH), 1,3-butadiene (1,3-C 4 H 6 ), benzene (C 6 H 6 ) and toluene (C 7 H 8 )) characteristics were experimentally investigated. Experimental results reveal that the application of pilot injection could improve combustion stability and fuel economy at high MSR, and it can also reduce regulated emissions CO, THC except NOx, and unregulated emissions tested in this study except CO 2 on M0 and M10 mode and toluene on M50 mode when compared with single injection cases. Increase pilot injection quantity and advance pilot injection timing both cause an increase in in-cylinder temperature and pressure before main combustion. The variations of pilot injection quantity and timing have less effect on regulated and unregulated gaseous emissions except NOx on M0 and M10 mode than higher MSR mode. With the rise of pilot quantity and the advance of pilot injection timing, HC and CO emissions decrease gradually but NOx emissions ascend; almost all unregulated emissions tested in this study on M30 mode reduce gradually; but when the MSR continues to rise to 50%, increasing emissions in unburned methanol, 1,3-butadiene and benzene are observed.

The reconstruction of the evolutionary history of a set of species is an important problem in classification and phylogenetics. Phylogenetic networks are a generalization of evolutionary trees that are used to represent histories for species that have undergone reticulate evolution, an important evolutionary force for many organisms (e.g. plants or viruses). In this paper, we present a novel approach to understanding the structure of networks that are not necessarily binary. More specifically, we define the concept of a closed set and show that the collection of closed sets of a network forms a hierarchy, and that this hierarchy can be deduced from either the subtrees or subnetworks on all 3-subsets. This allows us to also show that closed sets generalize the concept of the SN-sets of a binary network, sets which have proven very useful in elucidating the structure of binary networks. We also characterize the minimal closed sets (under set inclusion) for a special class of networks (2-terminal networks). Taken together, we anticipate that our results should be useful for the development of new phylogenetic network reconstruction algorithms.

Diesel/natural gas dual fuel (DNDF) and diesel/methanol dual fuel (DMDF) both have very good potential for use in alternative energy technologies for internal combustion engine. In this work, a comparative study was done to better understand the ignition characteristics of dual fuel combustion. First, their apparent ignition characteristics were studied experimentally on a constant volume chamber and apparent ignition delay correlations of DNDF and DMDF were built; second, chemical kinetic analysis was used to explain the phenomenon observed; lastly, an improved KIVA-3 V code was applied to study the importance of the chemical interaction between diesel and methane or methanol on diesel ignition timing in an engine. Experimental results showed that both methanol and methane addition retarded the ignition of diesel and methane retarded less, which was due to the dehydrogenation of methanol and methane, reaction CH3 + HO2 = CH3O + OH and the formation of H2O2. CFD numerical results showed that, in the DMDF ...

The methanol fuel used in diesel engine can reduce the emissions effectively, while the performance of engine keeps unchanged or even improved, the study of methanol used in diesel engine is paid more and more attention. Experiments were conducted in a 4-cylinder, turbocharged direct injection diesel engine fueled with diesel/methanol dual fuel (DMDF). In the experiment, methanol was injected into inlet duct by two low-pressure injectors, and directly mixed with intake air before entering the engine. The diesel fuel was directly injected in cylinder. The methanol substitution percent (MSP), intake air temperature and engine speed were changed to investigate the unevenness degree (UED) cylinder-to-cylinder. The experimental results show that the UED of engine increases with the increase of MSP under various engine loads, the COV pp increases with MSP under low engine load, and the COV pp curves vary very little under low MSP under 75% and 100% engine load, but a significantly rise appeared when MSP further increases. The UED gradually increases with the increase of intake air temperature, and it is usually less than 8% when intake air temperatures lower than 80 °C under 40% MSP. It suggests it is helpful to decrease the UED when the intake air temperature is in the range of 35–80 °C as MSP less than 40%. In addition, the difference of pressure between cylinder-to-cylinder decreases with the increase of engine speed, it can be considered raising MSP at high engine speed in order to gain better engine economy.

This paper presented the analytical results for the energy balance of a diesel methanol dual-fuel (DMDF) engine. A methanol injection system was fixed on a diesel engine and let the engine to alternate to run in either the pure diesel (D) mode or DMDF mode. Primary application of the thermal balance analysis was to investigate the cause of methanol replacement ratio SR changing with engine loads. Calculations were conducted within the control volume by comparing each of the energy terms. The results show that the reduction in cooling loss in DMDF mode is the dominant factor for low SR at high loads. While, a substantial increase in incomplete combustion loss in DMDF mode is the most important reason for high SR at low loads. Another application of the thermal balance analysis was to research the cooling loss in DMDF mode which implied unique characteristics. Therefore, the changing methanol replacement rate SP experiments were carried out in various engine loads, giving the results that the cooling loss in DMDF mode is always lower than that in D mode, and it goes even lower as SP rises. Finally, the effect of methanol temperature on SR was also investigated from the thermal balance point of view. The results reveal that the incomplete combustion loss decreasing in higher methanol temperature conditions will lead to a certain extent of SR drop. All these results are related to the energy flow distribution within the engine.

An experimental study was conducted to investigate into the cyclic variations of diesel methanol dual fuel (DMDF) combustion on a methanol fumigated diesel engine. The effects of engine load, methanol substitution percent (MSP), injection timing and intake temperature on cyclic combustion variation were studied experimentally. In-cylinder pressure of 100 consecutive combustion cycles were recorded for each test case. The coefficient of variation (COV) of maximum in-cylinder pressure ( p max ), maximum mass-averaged temperature ( T max ) and the indicated mean effective pressure (IMEP) were used to evaluate the cyclic variations of DMDF combustion. The results showed that the stability of DMDF combustion at high loads is comparable with that of neat diesel combustion. But the magnitude of cyclic DMDF combustion variations increases at light loads and more sensitive to the quantity of fumigated methanol. At light loads, COV p max and COV T max increase significantly and the crank angles corresponding to p max and T max scatter more around with the increment of MSP. IMEP values show obvious fluctuations and are distributed in a wider range with an increase in MSP. DMDF combustion stability is very sensitive to intake temperature. At light loads the COVs of all combustion parameters decrease with the increment of intake temperature. In contrast, at high loads, the effect of intake temperature on the cyclic variations is quite small while there is a significant increase in the COVs when the intake temperature exceeds over 75 °C. Overall, injection timings does no effect on cyclic variations of DMDF combustion.

In order to reveal the mechanism of high thermal efficiency of diesel methanol dual fuel (DMDF) engine, the exergy loss of DMDF engine under different operating parameters is investigated by using experimental and calculative methods based on the second law of thermodynamics. It has been found that DMDF has lower heat transfer exergy loss under all engine loads due to the high latent heat of vaporization of methanol which reduces the temperature in-cylinder. The exhaust exergy loss of DMDF under high load is relatively low compared with that of pure diesel fuel, while it is opposite to that under medium-low load. This is due to the influence of different temperatures on the incomplete combustion of methanol, which also indicates the main cause for the high thermal efficiency of DMDF under high load. Then studies on exergy analysis of DMDF under different temperatures were carried out. The studies covers the effect of intake air temperature on exergy efficiency, the heat transfer exergy loss and the exhaust physical exergy loss increase with rising intake air temperature as well. The results shows that the exhaust chemical exergy loss decreases significantly with the increase of intake air temperature and higher intake air temperature helps to improve exergy efficiency. Additionally, it is concluded that higher methanol temperature and cooling water temperature can improve exergy efficiency by decreasing exhaust exergy loss and combustion irreversible exergy loss. Finally, some feasible measures are proposed to further improve DMDF engine’s performance.

An important problem in phylogenetics is the construction of phylogenetic trees. One way to approach this problem, known as the supertree method, involves inferring a phylogenetic tree with leaves consisting of a set $X$ of species from a collection of trees, each having leaf-set some subset of $X$. In the 1980's characterizations, certain inference rules were given for when a collection of 4-leaved trees, one for each 4-element subset of $X$, can all be simultaneously displayed by a single supertree with leaf-set $X$. Recently, it has become of interest to extend such results to phylogenetic networks. These are a generalization of phylogenetic trees which can be used to represent reticulate evolution (where species can come together to form a new species). It has been shown that a certain type of phylogenetic network, called a level-1 network, can essentially be constructed from 4-leaved trees. However, the problem of providing appropriate inference rules for such networks remains unresolved. Here we show that by considering 4-leaved networks, called quarnets, as opposed to 4-leaved trees, it is possible to provide such rules. In particular, we show that these rules can be used to characterize when a collection of quarnets, one for each 4-element subset of $X$, can all be simultaneously displayed by a level-1 network with leaf-set $X$. The rules are an intriguing mixture of tree inference rules, and an inference rule for building up a cyclic ordering of $X$ from orderings on subsets of $X$ of size 4. This opens up several new directions of research for inferring phylogenetic networks from smaller ones, which could yield new algorithms for solving the supernetwork problem in phylogenetics., 17 pages

SummarySplit-networks are a generalization of phylogenetic trees that have proven to be a powerful tool in phylogenetics. Various ways have been developed for computing such networks, including split-decomposition, NeighborNet, QNet and FlatNJ. Some of these approaches are implemented in the user-friendly SplitsTree software package. However, to give the user the option to adjust and extend these approaches and to facilitate their integration into analysis pipelines, there is a need for robust, open-source implementations of associated data structures and algorithms. Here we present SPECTRE, a readily available, open-source library of data structures written in Java, that comes complete with new implementations of several pre-published algorithms and a basic interactive graphical interface for visualizing planar split networks. SPECTRE also supports the use of longer running algorithms by providing command line interfaces, which can be executed on servers or in High Performance Computing (HPC) environments.AvailabilityFull source code is available under the GPLv3 license at: https://github.com/maplesond/SPECTRESPECTRE’s core library is available from Maven Central at: https://mvnrepository.com/artifactuk.ac.uea.cmp.spectre/coreDocumentation is available at: http://spectre-suite-of-phylogenetic-tools-for-reticulate-evolution.readthedocs.io/en/latest/Contactsarah.bastkowski@earlham.ac.ukSupplementary Information (SI)Supplementary information is available at Bioinformatics online.

Diesel–methanol dual fuel (DMDF) engines at light loads suffer from low thermal efficiency and high unburned percentages of fuel. Pilot fuel injection timing and intake temperature are two important parameters which affect the combustion process in DMDF engines. In present experimental work, the combined effects of intake temperature and injection timing on the performance of a DMDF engine have been studied. The experiments were conducted on a methanol-fumigated diesel engine at 25% of full load and the results concerning performance, combustion characteristics and emissions were analyzed. Results show that the low efficiency at light loads can be improved significantly by raising the intake temperature and advancing the injection timing of direct-injected diesel. Increasing the intake temperature also significantly decreases the heat release rate of premixed combustion and increases the combustion rate of methanol burned by flame propagation. Flame propagation of the methanol–air mixture disappears gradually and DMDF combustion transforms into single stage combustion as the injection timing is retarded. When injection timing is retarded after 4.6° crank angle, misfire occurs at higher methanol substitute percent (MSP) and lower intake temperature, while the auto-ignition of methanol occurs at lower MSP and higher intake temperature. Under DMDF operation, soot and nitrogen oxides trade-off dilemma is completely broken at lower intake temperature and higher MSP.

Tandem duplication is a rearrangement process whereby a segment of DNA is replicated and proximally inserted. A sequence of these events is termed an evolution. Many different configurations can arise from such evolutions, generating some interesting combinatorial properties. Firstly, new DNA connections arising in an evolution can be algebraically represented with a word producing automaton. The number of words arising from n tandem duplications can then be recursively derived. Secondly, many distinct evolutions result in the same sequence of words. With the aid of a bi-colored 2d-tree, a Hasse diagram corresponding to a partially ordered set is constructed, for which the number of linear extensions equates to the number of evolutions generating a given word sequence. Thirdly, we implement some subtree prune and graft operations on this structure to show that the total number of possible evolutions arising from n tandem duplications is ? k = 1 n ( 4 k - ( 2 k + 1 ) ) . The space of structures arising from tandem duplication thus grows at a super-exponential rate with leading order term O ( 4 1 2 n 2 ) .

In phylogenetics, a common strategy used to construct an evolutionary tree for a set of species $$X$$ is to search in the space of all such trees for one that optimizes some given score function (such as the minimum evolution, parsimony or likelihood score). As this can be computationally intensive, it was recently proposed to restrict such searches to the set of all those trees that are compatible with some circular ordering of the set $$X$$ . To inform the design of efficient algorithms to perform such searches, it is therefore of interest to find bounds for the number of trees compatible with a fixed ordering in the neighborhood of a tree that is determined by certain tree operations commonly used to search for trees: the nearest neighbor interchange (nni), the subtree prune and regraft (spr) and the tree bisection and reconnection (tbr) operations. We show that the size of such a neighborhood of a binary tree associated with the nni operation is independent of the tree’s topology, but that this is not the case for the spr and tbr operations. We also give tight upper and lower bounds for the size of the neighborhood of a binary tree for the spr and tbr operations and characterize those trees for which these bounds are attained.

Phylogenetic networks are a generalization of evolutionary trees that are used by biologists to represent the evolution of organisms which have undergone reticulate evolution. Essentially, a phylogenetic network is a directed acyclic graph having a unique root in which the leaves are labelled by a given set of species. Recently, some approaches have been developed to construct phylogenetic networks from collections of networks on 2- and 3-leaved networks, which are known as binets and trinets, respectively. Here we study in more depth properties of collections of binets, one of the simplest possible types of networks into which a phylogenetic network can be decomposed. More specifically, we show that if a collection of level-1 binets is compatible with some binary network, then it is also compatible with a binary level-1 network. Our proofs are based on useful structural results concerning lowest stable ancestors in networks. In addition, we show that, although the binets do not determine the topology of the network, they do determine the number of reticulations in the network, which is one of its most important parameters. We also consider algorithmic questions concerning binets. We show that deciding whether an arbitrary set of binets is compatible with some network is at least as hard as the well-known graph isomorphism problem. However, if we restrict to level-1 binets, it is possible to decide in polynomial time whether there exists a binary network that displays all the binets. We also show that to find a network that displays a maximum number of the binets is NP-hard, but that there exists a simple polynomial-time 1/3-approximation algorithm for this problem. It is hoped that these results will eventually assist in the development of new methods for constructing phylogenetic networks from collections of smaller networks.

The Southern Ocean houses a diverse and productive community of organisms. Unicellular eukaryotic diatoms are the main primary producers in this environment, where photosynthesis is limited by low concentrations of dissolved iron and large seasonal fluctuations in light, temperature and the extent of sea ice. How diatoms have adapted to this extreme environment is largely unknown. Here we present insights into the genome evolution of a cold-Adapted diatom from the Southern Ocean, Fragilariopsis cylindrus, based on a comparison with temperate diatoms. We find that approximately 24.7 per cent of the diploid F. cylindrus genome consists of genetic loci with alleles that are highly divergent (15.1 megabases of the total genome size of 61.1 megabases). These divergent alleles were differentially expressed across environmental conditions, including darkness, low iron, freezing, elevated temperature and increased CO 2 . Alleles with the largest ratio of non-synonymous to synonymous nucleotide substitutions also show the most pronounced condition-dependent expression, suggesting a correlation between diversifying selection and allelic differentiation. Divergent alleles may be involved in adaptation to environmental fluctuations in the Southern Ocean. © 2017 Macmillan Publishers Limited, part of Springer Nature. All rights reserved.

Phylogenetic networks are a generalization of phylogenetic trees that allow for representation of reticulate evolution. Recently, a space of unrooted phylogenetic networks was introduced, where such a network is a connected graph in which every vertex has degree 1 or 3 and whose leaf-set is a fixed set $X$ of taxa. This space, denoted $\mathcal{N}(X)$, is defined in terms of two operations on networks -- the nearest neighbor interchange and triangle operations -- which can be used to transform any network with leaf set $X$ into any other network with that leaf set. In particular, it gives rise to a metric $d$ on $\mathcal N(X)$ which is given by the smallest number of operations required to transform one network in $\mathcal N(X)$ into another in $\mathcal N(X)$. The metric generalizes the well-known NNI-metric on phylogenetic trees which has been intensively studied in the literature. In this paper, we derive a bound for the metric $d$ as well as a related metric $d_{N\!N\!I}$ which arises when restricting $d$ to the subset of $\mathcal{N}(X)$ consisting of all networks with $2(|X|-1+i)$ vertices, $i \ge 1$. We also introduce two new metrics on networks -- the SPR and TBR metrics -- which generalize the metrics on phylogenetic trees with the same name and give bounds for these new metrics. We expect our results to eventually have applications to the development and understanding of network search algorithms., Comment: 17 pages, 5 figures. This version has a new figure to illustrate Lemma 3.2, and a new Corollary 5.7 that bounds the distance between networks in different tiers

In evolutionary biology, phylogenetic networks are constructed to represent the evolution of species in which reticulate events are thought to have occurred, such as recombination and hybridization. It is therefore useful to have efficiently computable metrics with which to systematically compare such networks. Through developing an optimal algorithm to enumerate all trinets displayed by a level-1 network (a type of network that is slightly more general than an evolutionary tree), here we propose a cubic-time algorithm to compute the trinet distance between two level-1 networks. Employing simulations, we also present a comparison between the trinet metric and the so-called Robinson-Foulds phylogenetic network metric restricted to level-1 networks. The algorithms described in this paper have been implemented in JAVA and are freely available at https://www.uea.ac.uk/computing/TriLoNet., Comment: 11pages, 5 figures

International audience; Binets and trinets are phylogenetic networks with two and three leaves, respectively. Here we consider the problem of deciding if there exists a binary level-1 phylogenetic network displaying a given set T of binary binets or trinets over a taxon set X , and constructing such a network whenever it exists. We show that this is NP-hard for trinets but polynomial-time solvable for binets. Moreover, we show that the problem is still polynomial-time solvable for inputs consisting of binets and trinets as long as the cycles in the trinets have size three. Finally, we present an O(3 |X | poly(|X |)) time algorithm for general sets of binets and trinets. The latter two algorithms generalise to instances containing level-1 networks with arbitrarily many leaves, and thus provide some of the first supernetwork algorithms for computing networks from a set of rooted phylogenetic networks. B Leo van Iersel

In this article we study the treewidth of the \emph{display graph}, an auxiliary graph structure obtained from the fusion of phylogenetic (i.e., evolutionary) trees at their leaves. Earlier work has shown that the treewidth of the display graph is bounded if the trees are in some formal sense topologically similar. Here we further expand upon this relationship. We analyse a number of reduction rules which are commonly used in the phylogenetics literature to obtain fixed parameter tractable algorithms. In some cases (the \emph{subtree} reduction) the reduction rules behave similarly with respect to treewidth, while others (the \emph{cluster} reduction) behave very differently, and the behaviour of the \emph{chain reduction} is particularly intriguing because of its link with graph separators and forbidden minors. We also show that the gap between treewidth and Tree Bisection and Reconnect (TBR) distance can be infinitely large, and that unlike, for example, planar graphs the treewidth of the display graph can be as much as linear in its number of vertices. On a slightly different note we show that if a display graph is formed from the fusion of a phylogenetic network and a tree, rather than from two trees, the treewidth of the display graph is bounded whenever the tree can be topologically embedded ("displayed") within the network. This opens the door to the formulation of the display problem in Monadic Second Order Logic (MSOL). A number of other auxiliary results are given. We conclude with a discussion and list a number of open problems., Comment: 29 pages, 12 figures