Mitonuclear Coevolution,but not Nuclear Compensation,Drives Evolution of OXPHOS Complexes in Bivalves

In Metazoa, four out of five complexes involved in oxidative phosphorylation (OXPHOS) are formed by subunits encoded by both the mitochondrial (mtDNA) and nuclear (nuDNA) genomes, leading to the expectation of mitonuclear coevolution. Previous studies have supported coadaptation of mitochondria-encoded (mtOXPHOS) and nuclear-encoded OXPHOS (nuOXPHOS) subunits, often specifically interpreted with regard to the “nuclear compensation hypothesis,” a specific form of mitonuclear coevolution where nuclear genes compensate for deleterious mitochondrial mutations due to less efficient mitochondrial selection. In this study, we analyzed patterns of sequence evolution of 79 OXPHOS subunits in 31 bivalve species, a taxon showing extraordinary mtDNA variability and including species with “doubly uniparental” mtDNA inheritance. Our data showed strong and clear signals of mitonuclear coevolution. NuOXPHOS subunits had concordant topologies with mtOXPHOS subunits, contrary to previous phylogenies based on nuclear genes lacking mt interactions. Evolutionary rates between mt and nuOXPHOS subunits were also highly correlated compared with non-OXPHO-interacting nuclear genes. Nuclear subunits of chimeric OXPHOS complexes (I, III, IV, and V) also had higher dN/dS ratios than Complex II, which is formed exclusively by nuDNA-encoded subunits. However, we did not find evidence of nuclear compensation: mitochondria-encoded subunits showed similar dN/dS ratios compared with nuclear-encoded subunits, contrary to most previously studied bilaterian animals. Moreover, no site-specific signals of compensatory positive selection were detected in nuOXPHOS genes. Our analyses extend the evidence for mitonuclear coevolution to a new taxonomic group, but we propose a reconsideration of the nuclear compensation hypothesis.     

Ecological and geological processes impacting speciation modes drive the formation of wide-range disjunctions within tribe Putorieae (Rubiaceae)

Wide-range geographically discontinuous distributions have long intrigued scientists. We explore the role of ecology, geology, and dispersal in the formation of these large-scale disjunctions, using the angiosperm tribe Putorieae (Rubiaceae) as a case study. From DNA sequences of nuclear ITS and six plastid markers, we inferred a phylogeny with 65% of all known Putorieae species. Divergence times, ancestral ranges, and diversification rate shifts were then estimated using Bayesian inference. We further explored species climatic tolerances and performed ancestral niche reconstruction to discriminate among alternative speciation modes, including geographical and ecological vicariance, and ecogeographical, ecological, and dispersal-mediated speciation. As a result, we identified seven major clades in Putorieae, some of which exhibit striking geographical disjunctions, matching the Rand Flora pattern, with sister species in the Canary Islands andeastern and southern Africa. Initial diversification within the tribe occurred in the early Miocene, coincident with a period of climate warming; however, most clades diverged within the last 10 Myr. Aridification and high extinction rates, coupled with ecological vicariance, explain the oldest disjunctions. Adaptation to new environmental conditions, after allopatry, is observed in several clades. Dispersal, either long-distance or via corridors made available by mountain uplift, is behind the most recent disjunctions. Some of these events were followed by ecological speciation and rapid diversification, with species becoming adapted to xeric or increasingly colder continental climates. We show that an integrative approach may help discriminate among speciation modes invoked to explain disjunctions at macroevolutionary time scales, even when extinction has erased the signature of past events.     

Better alone? A demographic case study of the hemiparasite Castilleja tenuiflora (Orobanchaceae): A first approximation

Castilleja tenuiflora is a facultative root hemiparasitic plant that has colonized a disturbed lava field in central Mexico. To determine the effects of hemiparasitism on the population dynamics of the parasite, we identified a set of potential hosts and quantified their effects on the vital rates of C. tenuiflora during 2016–2018. Connections between the roots of the hemiparasite and the hosts were confirmed with a scanning electron microscope. Annual matrices considering two conditions (with and without potential hosts) were built based on vital rates for each year, and annual stochastic finite rate growth rates (λ_s) were calculated. Plants produced more reproductive structures with hosts than without hosts. A Life Table Response Experiment (LTRE) was performed to compare the contributions of vital rates between conditions. We identified 19 species of potential hosts for this generalist hemiparasite. Stochastic lambda with hosts λ_s = 1.02 (CI = 0.9999, 1.1) tended to be higher than without them λ_s = 0.9503 (CI = 0.9055, 0.9981). The highest elasticity values correspond to survival. LTRE indicated that the most important parameters are survival and fecundity; the total contribution of fecundity (0.0192) to the difference in growth was three times lower than that of survival (0.0603). Piqueria trinervia was the most abundant host, and C. tenuiflora had a higher lambda with it than with other species. Individuals can grow alone, but hosts can have a positive effect on the vital parameters of C. tenuiflora and on λ.     

Photon beam energy dependent single-arc volumetric modulated arc optimization

PurposeThe purpose of this work was to present a new single-arc mixed photon (6&18MV) VMAT (SAMP) optimization framework that concurrently optimizes for two photon energies with corresponding partial arc lengths.Methods and materialsOwing to simultaneous optimization of energy dependent intensity maps and corresponding arc locations, the proposed model poses nonlinearity. Unique relaxation constraints based on McCormick approximations were introduced for linearization. Energy dependent intensity maps were then decomposed to generate apertures. Feasibility of the proposed framework was tested on a sample of ten prostate cancer cases with lateral separation ranging from 34 cm (case no.1) to 52 cm (case no.6). The SAMP plans were compared against single energy (6MV) VMAT (SE) plans through dose volume histograms (DVHs) and radiobiological parameters including normal tissue complication probability (NTCP) and equivalent uniform dose (EUD).ResultsThe contribution of higher energy photon beam optimized by the algorithm demonstrated an increase for cases with a lateral separation >40 cm. SAMP–VMAT notably improved bladder and rectum sparing in large size cases. Compared to single energy, SAMP–VMAT plans reduced bladder and rectum NTCP in cases with large lateral separation. With the exception of one case, SAMP–VMAT either improved or maintained femoral heads compared to SE–VMAT. SAMP–VMAT reduced the nontarget tissue integral dose in all ten cases.ConclusionsA single-arc VMAT optimization framework comprising mixed photon energy partial arcs was presented. Overall results underline the feasibility and potential of the proposed approach for improving OAR sparing in large size patients without compromising the target homogeneity and coverage.     

Viscoelastic creep induced by repetitive spine flexion and its relationship to dynamic spine stability

Repetitive trunk flexion elicits passive tissue creep, which has been hypothesized to compromise spine stability. The current investigation determined if increased spine flexion angle at the onset of flexion relaxation (FR) in the lumbar extensor musculature was associated with altered dynamic stability of spine kinematics. Twelve male participants performed 125 consecutive cycles of full forward trunk flexion. Spine kinematics and lumbar erector spinae (LES) electromyographic (EMG) activity were obtained throughout the repetitive trunk flexion trial. Dynamic stability was evaluated with maximum finite-time Lyapunov exponents over five sequential blocks of 25 cycles. Spine flexion angle at FR onset, and peak LES EMG activity were determined at baseline and every 25th cycle. Spine flexion angle at FR increased on average by 1.7° after baseline with significant increases of 1.7° and 2.4° at the 50th and 100th cycles. Maximum finite-time Lyapunov exponents demonstrated a transient, non-statistically significant, increase between cycles 26 and 50 followed by a recovery to baseline over the remainder of the repetitive trunk flexion cycles. Recovery of dynamic stability may be the consequence of increased active spine stiffness demonstrated by the non-significant increase in peak LES EMG that occurred as the repetitive trunk flexion progressed.     

Desulfurization of dibenzothiophene using the 4S enzymatic route: Influence of operational conditions on initial reaction rates

The influence of operational conditions (pH, temperature and oxygen transfer rate) on the initial reaction rates of the four reactions involved in the 4S biodesulfurization route of dibenzothiophenes (DBT) has been studied. The bioprocess was carried out using a genetically modified organism, Pseudomonas putida CECT 5279. The rates of the four reactions were calculated from the rates of production of different compounds involved in the 4S pathway, by matrix manipulation. The initial (zero time) reaction rates showed a slight dependence on oxygen transfer rate. Temperature and pH were optimal at 30°C and 9, respectively, temperature being the most important variable. This study also identifies the last reaction as the limiting step in the pathway.     

A comparison of cartilage stress-relaxation models in unconfined compression: QLV and stretched exponential in combination with fluid flow

Cartilage exhibits nonlinear viscoelastic behaviour. Various models have been proposed to explain cartilage stress relaxation, but it is unclear whether explicit modelling of fluid flow in unconfined compression is needed. This study compared Fung's quasi-linear viscoelastic (QLV) model with a stretched-exponential model of cartilage stress relaxation and examined each of these models both alone and in combination with a fluid-flow model in unconfined compression. Cartilage explants were harvested from bovine calf patellofemoral joints and equilibrated in tissue culture for 5 days before stress-relaxation testing in unconfined compression at 5% nominal strain. The stretched exponential models fit as well as the QLV models. Furthermore, the average stretched exponential relaxation time determined by this model lies within the range of experimentally measured relaxation times for extracted proteoglycan aggregates, consistent with the hypothesis that the stretched exponential model represents polymeric mechanisms of cartilage viscoelasticity.