Patterns of population differentiation in annual killifishes from the Paraná–Uruguay–La Plata Basin: the role of vicariance and dispersal

Aim To elucidate the role of vicariance versus dispersal at the microevolutionary scale in annual killifish populations belonging to the Austrolebias bellottii species complex (Rivulidae). Within this complex, A. bellottii and A. apaii have low vagility and occur widely within the study area, making them excellent models for testing biogeographic hypotheses of differentiation. Location South America, in the Paraná–Uruguay–La Plata river basin. Methods Molecular data and morphometric analyses were used to reconstruct the phylogeographic history and morphological variation of 24 populations of two taxa of the A. bellottii species complex. Phylogenetic analyses using maximum likelihood (ML) and Bayesian inference (BI) model‐based methods, estimates of clade divergence times implemented in beast , non‐metric multidimensional scaling, analysis of molecular variance results, and morphological analyses elucidated the role of vicariance versus dispersal hypotheses in population differentiation in the aforementioned river basin. Results In the A. bellottii species complex from the Paraná–Uruguay–La Plata river basin, past allopatric fragmentation from vicariance events seems to be the most plausible scenario for diversification since the Late Miocene and more recently since the Plio‐Pleistocene. The Plio‐Pleistocene vicariance produced the differentiation of three major clades in A. bellottii populations. One clade from the eastern Uruguay River drainage was separated from another in western Uruguay and the Paraná–La Plata River drainages. A later vicariance event split populations to the south (lower Paraná–La Plata Basin) and north (middle Paraná) of the western Paraná River drainage. However, our results do not exclude the possibility of dispersal events among A. bellottii populations from both the Uruguay and Paraná river drainages, which could occur in these river basins during hypothesized connectivity cycles of the Late Pliocene and Pleistocene. Main conclusions Past allopatric fragmentation caused by different vicariance events seems to be the main driver of diversification in the A. bellottii species complex since the Plio‐Pleistocene. However, the current molecular data suggest that populations from both drainages of the Paraná–Uruguay rivers may have experienced cycles of connectivity during the Pleistocene, perhaps including multiple vicariance or dispersal events from populations located in the western lower Uruguay River drainage, which encompassed climatic and geological changes in the Paraná–Uruguay–La Plata Basin.     

Dispersal vs. vicariance in the Mediterranean: historical biogeography of the Palearctic Pachydeminae (Coleoptera, Scarabaeoidea)

Aim The geological evolution of the Mediterranean region is largely the result of the Tertiary collision of the African and Eurasian Plates, but also a mosaic of migrating island arcs, fragmenting tectonic belts, and extending back‐arc basins. Such complex paleogeography has resulted in a ‘reticulate’ biogeographical history, in which Mediterranean biotas repeatedly fragmented and merged as dispersal barriers appeared and disappeared through time. In this study, dispersal‐vicariance analysis (DIVA) is used to assess the relative role played by dispersal and vicariance in shaping distribution patterns in the beetle subfamily Pachydeminae Reitter, 1902 (Scarabaeoidea), an example of east–west Mediterranean disjunction. Location The Mediterranean region, including North Africa, the western Mediterranean, Balkans–Anatolia, Middle East, Caucasus, the Iranian Plateau, and Central Asia. Methods A phylogenetic hypothesis of the Palearctic genera of Pachydeminae in conjunction with distributional data was analysed using DIVA. This method reconstructs the ancestral distribution in a given phylogeny based on the vicariance model, while allowing dispersal and extinction to occur. Unlike other methods, DIVA does not enforce area relationships to conform to a hierarchical ‘area cladogram’, so it can be used to reconstruct ‘reticulate’ biogeographical scenarios. Results Optimal reconstructions, requiring 23 dispersal events, suggest that the ancestor of Pachydeminae was originally present in the south‐east Mediterranean region. Basal splitting within the subfamily was caused by vicariance events related to the late Tertiary collision of the African microplates Apulia and Arabia with Eurasia, and the resultant arise of successive dispersal barriers (e.g. the Red Sea, the Zagros Mountains). Subsequent diversification in Pachydeminae involved multiple speciation events within the Middle East and Iran–Afghanistan regions, which gave rise to the least speciose genera of Pachydeminae (e.g. Otoclinius Brenske, 1896). Finally, the presence of Pachydeminae in the western Mediterranean region seems to be the result of a recent dispersal event. The ancestor of the Iberian genera Ceramida Baraud, 1987 and Elaphocera Gené, 1836 probably dispersed from the Middle East to the Iberian Peninsula across North Africa and the Gibraltar Strait during the ‘Messinian salinity crisis’ at the end of the Miocene. Main conclusions Although the basal diversification of Pachydeminae around the Mediterranean appears to be related to vicariance events linked to the geological formation of the Mediterranean Basin, dispersal has also played a very important role. Nearly 38% of the speciation events in the phylogeny resulted from dispersal to a new area followed by allopatric speciation between lineages. Relationships between western and eastern Mediterranean disjuncts are usually explained by dispersal through Central Europe. The biogeographical history of the Pachydeminae corroborates other biogeographical studies that consider North Africa to be an alternative dispersal route by which Mediterranean taxa could have achieved circum‐Mediterranean distributions.     

An evaluation of new parsimony‐based versus parametric inference methods in biogeography: a case study using the globally distributed plant family Sapindaceae

Aim Recently developed parametric methods in historical biogeography allow researchers to integrate temporal and palaeogeographical information into the reconstruction of biogeographical scenarios, thus overcoming a known bias of parsimony‐based approaches. Here, we compare a parametric method, dispersal–extinction–cladogenesis (DEC), against a parsimony‐based method, dispersal–vicariance analysis (DIVA), which does not incorporate branch lengths but accounts for phylogenetic uncertainty through a Bayesian empirical approach (Bayes‐DIVA). We analyse the benefits and limitations of each method using the cosmopolitan plant family Sapindaceae as a case study. Location World‐wide. Methods Phylogenetic relationships were estimated by Bayesian inference on a large dataset representing generic diversity within Sapindaceae. Lineage divergence times were estimated by penalized likelihood over a sample of trees from the posterior distribution of the phylogeny to account for dating uncertainty in biogeographical reconstructions. We compared biogeographical scenarios between Bayes‐DIVA and two different DEC models: one with no geological constraints and another that employed a stratified palaeogeographical model in which dispersal rates were scaled according to area connectivity across four time slices, reflecting the changing continental configuration over the last 110 million years. Results Despite differences in the underlying biogeographical model, Bayes‐DIVA and DEC inferred similar biogeographical scenarios. The main differences were: (1) in the timing of dispersal events – which in Bayes‐DIVA sometimes conflicts with palaeogeographical information, and (2) in the lower frequency of terminal dispersal events inferred by DEC. Uncertainty in divergence time estimations influenced both the inference of ancestral ranges and the decisiveness with which an area can be assigned to a node. Main conclusions By considering lineage divergence times, the DEC method gives more accurate reconstructions that are in agreement with palaeogeographical evidence. In contrast, Bayes‐DIVA showed the highest decisiveness in unequivocally reconstructing ancestral ranges, probably reflecting its ability to integrate phylogenetic uncertainty. Care should be taken in defining the palaeogeographical model in DEC because of the possibility of overestimating the frequency of extinction events, or of inferring ancestral ranges that are outside the extant species ranges, owing to dispersal constraints enforced by the model. The wide‐spanning spatial and temporal model proposed here could prove useful for testing large‐scale biogeographical patterns in plants.     

Are the Northern Andes a species pump for Neotropical birds? Phylogenetics and biogeography of a clade of Neotropical tanagers (Aves: Thraupini)

Aim We used mitochondrial DNA sequence data to reconstruct the phylogeny of a large clade of tanagers (Aves: Thraupini). We used the phylogeny of this Neotropical bird group to identify areas of vicariance, reconstruct ancestral zoogeographical areas and elevational distributions, and to investigate the correspondence of geological events to speciation events. Location The species investigated are found in 18 of the 22 zoogeographical regions of South America, Central America and the Caribbean islands; therefore, we were able to use the phylogeny to address the biogeographical history of the entire region. Methods Molecular sequence data were gathered from two mitochondrial markers (cytochrome b and ND2) and analysed using Bayesian and maximum‐likelihood approaches. Dispersal–vicariance analysis (DIVA) was used to reconstruct zoogeographical areas and elevational distributions. A Bayesian framework was also used to address changes in elevation during the evolutionary history of the group. Results Our phylogeny was similar to previous tanager phylogenies constructed using fewer species; however, we identified three genera that are not monophyletic and uncovered high levels of sequence divergence within some species. DIVA identified early diverging nodes as having a Northern Andean distribution, and the most recent common ancestor of the species included in this study occurred at high elevations. Most speciation events occurred either within highland areas or within lowland areas, with few exchanges occurring between the highlands and lowlands. The Northern Andes has been a source for lineages in other regions, with more dispersals out of this area relative to dispersals into this area. Most of the dispersals out of the Northern Andes were dispersals into the Central Andes; however, a few key dispersal events were identified out of the Andes and into other zoogeographical regions. Main conclusions The timing of diversification of these tanagers correlates well with the main uplift of the Northern Andes, with the highest rate of speciation occurring during this timeframe. Central American tanagers included in this study originated from South American lineages, and the timing of their dispersal into Central America coincides with or post‐dates the completion of the Panamanian isthmus.     

Divergence time uncertainty and historical biogeography reconstruction – an example from Urophylleae (Rubiaceae)

Aim When hypotheses of historical biogeography are evaluated, age estimates of individual nodes in a phylogeny often have a direct impact on what explanation is concluded to be most likely. Confidence intervals of estimated divergence times obtained in molecular dating analyses are usually very large, but the uncertainty is rarely incorporated in biogeographical analyses. The aim of this study is to use the group Urophylleae, which has a disjunct pantropical distribution, to explore how the uncertainty in estimated divergence times affects conclusions in biogeographical analysis. Two hypotheses are evaluated: (1) long‐distance dispersal from Africa to Asia and the Neotropics, and (2) a continuous distribution in the boreotropics, probably involving migration across the North Atlantic Land Bridge, followed by isolation in equatorial refugia. Location Tropical and subtropical Asia, tropical Africa, and central and southern tropical America. Methods This study uses parsimony and Bayesian phylogenetic analyses of chloroplast DNA and nuclear ribosomal DNA data from 56 ingroup species, beast molecular dating and a Bayesian approach to dispersal–vicariance analysis (Bayes‐DIVA) to reconstruct the ancestral area of the group, and the dispersal–extinction–cladogenesis method to test biogeographical hypotheses. Results When the two models of geographic range evolution were compared using the maximum likelihood (ML) tree with mean estimates of divergence times, boreotropical migration was indicated to be much more likely than long‐distance dispersal. Analyses of a large sample of dated phylogenies did, however, show that this result was not consistent. The age estimate of one specific node had a major impact on likelihood values and on which model performed best. The results show that boreotropical migration provides a slightly better explanation of the geographical distribution patterns of extant Urophylleae than long‐distance dispersal. Main conclusions This study shows that results from biogeographical analyses based on single phylogenetic trees, such as a ML or consensus tree, can be misleading, and that it may be very important to take the uncertainty in age estimates into account. Methods that account for the uncertainty in topology, branch lengths and estimated divergence times are not commonly used in biogeographical inference today but should definitely be preferred in order to avoid unwarranted conclusions.     

Molecular systematics and biogeography of Wisteria inferred from nucleotide sequences of nuclear and plastid genes

Previous molecular phylogenetic studies of Fabaceae indicated that species of Wisteria, an intercontinental disjunct genus between eastern Asia and eastern North America, formed a clade derived from within Callerya. However, interspecific relationships were not well resolved or supported. In this study, we used sequences of the nuclear ribosomal DNA internal transcribed spacer region and the chloroplast gene matK to examine interspecific relationships and explore implications of the phylogeny for the systematics and biogeography of Wisteria. Our results showed that Wisteria with deciduous leaves and racemose inflorescences formed a strongly supported clade derived from within the paraphyletic Callerya. Afgekia was also found to be included within Callerya. Therefore, our data support the merger ofAfgekia, Callerya, and Wisteria. The phylogenetic pattern suggested that the deciduousness in Wisteria may be a derived trait likely in response to temperate climate, and the racemose inflorescences in the Afgekia–Callerya–Wisteria clade may have evolved from panicles. Our study also provided strong support for the sister relationship of the North American and eastern Asian species of Wisteria. In the Asian clade, Wisteria brachybotrys Siebold & Zucc. of Japan was sister to the clade containing W. floribunda (Willd.) DC of Japan and Korea, and W. sinensis (Sims) Sweet of China. However, our data offered weak support for the sister relationship ofW. floribunda and W. sinensis. Our divergence time and biogeographic analyses suggested that the eastern Asian–North American disjunction in Wisteria may have occurred through a dispersal event in the middle Miocene (13.4 Mya) from the Old World to the New World across the Bering land bridge followed by vicariance in the late Miocene (6.8 Mya). This study added another example to the “out of Asia” migration for the eastern Asian–eastern North American disjunction.     

Diversification of subgenus Calathus (Coleoptera: Carabidae) in the Mediterranean region – glacial refugia and taxon pulses

Aim To investigate the effects of Pleistocene climatic variations on the diversification rate of the subgenus Calathus (Coleoptera: Carabidae), and to estimate the role of vicariance and dispersal for explaining current distributional patterns. Location Western Palaearctic Region, particularly the Mediterranean Basin. Methods Fragments of the mitochondrial cox1–cox2 and the nuclear 28S and EF1α genes were analysed by Bayesian inference. Lineage divergence times were estimated using a Bayesian relaxed molecular clock. Three diversification rate analyses were conducted, namely gamma (γ)‐statistic, birth–death likelihood (BDL) test and survival analyses, in order to test departures from a constant rate model of diversification. A Bayesian approach to dispersal–vicariance analysis was developed to reconstruct the most probable ancestral area of subgenus Calathus and subsequent events of dispersal and colonization. Results A constant rate of speciation events from the late Miocene onwards was found for the subgenus Calathus, whereas recent Pleistocene climatic oscillations played an important role only in shaping intraspecific diversity. Overall diversification patterns for the subgenus are best explained by at least four westward dispersal events from the eastern Mediterranean Basin. Three distinct phylogroups were found for the widely distributed Calathus fuscipes. Incongruence between mitochondrial and nuclear loci was found for a number of species. Main conclusions Diversification analyses suggest either a constant rate of diversification (BDL analysis) or a decrease in diversification rates for the subgenus (survival or γ‐statistics analyses), but not an increase related to the effects of glaciation cycles. Diversification patterns in the subgenus Calathus agree with predictions of the taxon pulse model. From the middle Miocene onwards the Anatolian Peninsula was possibly the main centre of diversification, with successive dispersal events towards the western Mediterranean Basin. Range expansion and secondary contact zones are postulated between members of different phylogroups in C. fuscipes.     

Historical biogeography of the coffee family (Rubiaceae,Gentianales) in Madagascar: case studies from the tribes Knoxieae,Naucleeae, Paederieae and Vanguerieae

Aim In Madagascar the family Rubiaceae includes an estimated 650 species representing 95 genera. As many as 98% of the species and 30% of the genera are endemic. Several factors make the Rubiaceae a model system for developing an understanding of the origins of the Malagasy flora. Ancestral area distributions are explicitly reconstructed for four tribes (Knoxieae, Naucleeae, Paederieae and Vanguerieae) with the aim of understanding how many times, and from where, these groups have originated in Madagascar. Location Indian Ocean Basin, with a focus on Madagascar. Methods Bayesian phylogenetic analyses are conducted on the four tribes. The results are used for reconstructing ancestral areas using dispersal–vicariance analyses. Phylogenetic uncertainties in the reconstructions are accounted for by conducting all analyses on the posterior distribution from the analyses. Results Altogether, 11 arrivals in Madagascar (one in Paederieae, five in Knoxieae, three in Vanguerieae, and two in Naucleeae) are reconstructed. The most common pattern is a dispersal event (followed by vicariance) from Eastern Tropical Africa. The Naucleeae and Paederieae in Madagascar differ and originate from Asia. Numerous out‐of‐Madagascar dispersals, mainly in the dioecious Vanguerieae, are reconstructed. Main conclusions The four tribes arrived several times in Madagascar via dispersal events from Eastern Tropical Africa, Southern Africa and Tropical Asia. The presence of monophyletic groups that include a number of species only found in Madagascar indicates that much endemism in the tribes results from speciation events occurring well after their arrival in Madagascar. Madagascar is the source of origin for almost all Rubiaceae found on the neighbouring islands of the Comoros, Mascarenes and Seychelles.     

Evolution of the Madrean–Tethyan disjunctions and the North and South American amphitropical disjunctions in plants

Abstract The present paper reviews advances in the study of two major intercontinental disjunct biogeographic patterns: (i) between Eurasian and western North American deserts with the Mediterranean climate (the Madrean–Tethyan disjunctions); and (ii) between the temperate regions of North and South America (the amphitropical disjunctions). Both disjunct patterns have multiple times of origin. The amphitropical disjunctions have largely resulted from long‐distance dispersal, primarily from the Miocene to the Holocene, with available data indicating that most lineages dispersed from North to South America. Results of recent studies on the Mediterranean disjuncts between the deserts of Eurasia and western North America support the multiple modes of origin and are mostly consistent with hypotheses of long‐distance dispersal and the North Atlantic migration. Axelrod's Madrean–Tethyan hypothesis, which implies vicariance between the two regions in the early Tertiary, has been favored by a few studies. The Beringian migration corridor for semiarid taxa is also supported in some cases.     

Phylogenomics,biogeography, and evolution of morphology and ecological niche of the eastern Asian–eastern North American Nyssa (Nyssaceae)

Nyssa (Nyssaceae, Cornales) represents a classical example of the well‐known eastern Asian–eastern North American floristic disjunction. The genus consists of three species in eastern Asia, four species in eastern North America, and one species in Central America. Species of the genus are ecologically important trees in eastern North American and eastern Asian forests. The distribution of living species and a rich fossil record of the genus make it an excellent model for understanding the origin and evolution of the eastern Asian–eastern North American floristic disjunction. However, despite the small number of species, relationships within the genus have remained unclear and have not been elucidated using a molecular approach. Here, we integrate data from 48 nuclear genes, fossils, morphology, and ecological niche to resolve species relationships, elucidate its biogeographical history, and investigate the evolution of morphology and ecological niches, aiming at a better understanding of the well‐known EA–ENA floristic disjunction. Results showed that the Central American (CAM) Nyssa talamancana was sister to the remaining species, which were divided among three, rapidly diversified subclades. Estimated divergence times and biogeographical history suggested that Nyssa had an ancestral range in Eurasia and western North America in the late Paleocene. The rapid diversification occurred in the early Eocene, followed by multiple dispersals between and within the Erasian and North American continents. The genus experienced two major episodes of extinction in the early Oligocene and end of Neogene, respectively. The Central American N. talamancana represents a relic lineage of the boreotropical flora in the Paleocene/Eocene boundary that once diversified in western North America. The results supported the importance of both the North Atlantic land bridge and the Bering land bridge (BLB) for the Paleogene dispersals of Nyssa and the Neogene dispersals, respectively, as well as the role of Central America as refugia of the Paleogene flora. The total‐evidence‐based dated phylogeny suggested that the pattern of macroevolution of Nyssa coincided with paleoclimatic changes. We found a number of evolutionary changes in morphology (including wood anatomy and leaf traits) and ecological niches (precipitation and temperature) between the EA–ENA disjunct, supporting the ecological selection driving trait evolutions after geographic isolation. We also demonstrated challenges in phylogenomic studies of lineages with rapid diversification histories. The concatenation of gene data can lead to inference of strongly supported relationships incongruent with the species tree. However, conflicts in gene genealogies did not seem to impose a strong effect on divergence time dating in our case. Furthermore, we demonstrated that rapid diversification events may not be recovered in the divergence time dating analysis using BEAST if critical fossil constraints of the relevant nodes are not available. Our study provides an example of complex bidirectional exchanges of plants between Eurasia and North America in the Paleogene, but “out of Asia” migrations in the Neogene, to explain the present disjunct distribution of Nyssa in EA and ENA.     

Electrochemical analysis of copper ion using a Gly–Gly–His tripeptide modified poly(3-thiopheneacetic acid) biosensor

A novel biosensor harnessing a conducting polymer functionalized with a copper ion specific peptide proved to be highly effective for electrochemical analysis of copper ions. The developed sensor comprised a transducer based on a conducting polymer (poly(3-thiopheneacetic acid)) electrode and a probe (tripeptide, Gly–Gly–His) selectively cognitive of copper ions. For functionalization of the electrode, the carboxylic group of the polymer was covalently coupled with the amine group of the tripeptide, and its structural features were confirmed by X-ray photoelectron spectroscopy (XPS) and attenuated total reflection infrared (ATR-IR) spectroscopy. The peptide modified polythiophene biosensor was used for the electrochemical analysis of various trace metal ions by square wave voltammetry. The electrode was found to be highly sensitive and selective to Cu²⁺ in the range of 0.02–20 μM with almost no cross binding to other metal ions such as Ni²⁺ and Pb²⁺. Furthermore, the developed sensor exhibited a high stability and reproducibility despite the repeated use of the sensor electrode and probe. With the advent of more diverse affinity bioprobes specific towards a broad range of analytes, the demonstrated strategy harnessing peptide modified polythiophene biosensor is likely to provide an excellent platform for the selective determination of trace amount of analytes whose detection is otherwise cumbersome.     

In silico structure–function analysis of E. cloacae nitroreductase

Reduction, catalyzed by the bacterial nitroreductases, is the quintessential first step in the biodegradation of a variety of nitroaromatic compounds from contaminated waters and soil. The Enterobacter cloacae nitroreductase (EcNR) enzyme is considered as a prospective biotechnological tool for bioremediation of hazardous nitroaromatic compounds. Using diverse computational methods, we obtain insights into the structural basis of activity and mechanism of its function. We have performed molecular dynamics simulation of EcNR in three different states (free EcNR in oxidized form, fully reduced EcNR with benzoate inhibitor and fully reduced EcNR with nitrobenzene) in explicit solvent and with full electrostatics. Principal Component Analysis (PCA) of the variance‐covariance matrix showed that the complexed nitroreductase becomes more flexible overall upon complexation, particularly helix H6, in the vicinity of the binding site. A multiple sequence alignment was also constructed in order to examine positional constraints on substitution in EcNR. Five regions which are highly conserved within the flavin mononucleotide (FMN) binding site were identified. Obtained results and their implications for EcNR functioning are discussed, and new plausible mechanism has been proposed. Proteins 2012; © 2012 Wiley Periodicals, Inc.     

Determination of clomipramine by flow‐injection analysis with acidic potassium permanganate–formic acid chemiluminescence detection

A sensitive and simple chemiluminescent (CL) method for the determination of clomipramine has been developed by combining the flow‐injection analysis (FIA) technique, which is based on the CL intensity generated from the redox reaction of potassium permanganate (KMnO₄)–formic acid in sulphuric acid (H₂SO₄) medium. Under the optimum conditions, the linear range for the determination of clomipramine was 0.04–4 µg/mL, with a correlation coefficient of 0.9988 (n = 10) and a detection limit of 0.008 µg/mL (3σ), and the relative standard deviation (RSD) for 2.0 µg/mL clomipramine (n = 11) is 1.26%. The proposed method has been successfully applied to the determination of the studied clomipramine in pharmaceutical preparations. The possible reaction mechanism is discussed. Copyright © 2011 John Wiley & Sons, Ltd.     

Solid‐phase extraction with metal–organic frameworks for the analysis of chiral compounds

Metal–organic frameworks (MOFs) are excellent porous materials with nanoscale cavities and high surface areas, which make them promising as novel adsorbents in solid‐phase extraction (SPE). In this article we report a new application of the chiral MOF [Zn₂(D‐Cam)₂(4,4′‐bpy)]_n in SPE used for the measurement of the enantiomeric excess (ee) of (±)‐1,1′‐bi‐2‐naphthol. Several important experimental parameters that may influence the extraction efficiency were investigated and optimized. Under the optimum conditions, a good linearity (R² > 0.999) was found between the ee value and the reciprocal of the peak areas. When compared with the actual ee measured using chiral HPLC, the SPE‐based assay also showed good accuracy and precision. The results showed that SPE based on chiral MOFs as adsorbents is a simple and effective method for the determination of the ee values of chiral compounds.     

A cost–benefit analysis of the environmental taxation policy in China: A frontier analysis‐based environmentally extended input–output optimization method

China's high‐speed economic development and reliance on overconsumption of natural resources have led to serious environmental pollution. Environmental taxation is seen as an effective economic tool to help mitigate air pollution. In order to assess the effects of different scenarios of environmental taxation policies, we propose a frontier‐based environmentally extended input–output optimization model with explicit emission abatement sectors to reflect the inputs and benefits of abatement. Frontier analysis ensures policy scenarios are assessed under the same technical efficiency benchmark, while input–output analysis depicts the wide range of economic transactions among sectors of an economy. Four scenarios are considered in this study, which are increasing specific tax rates of SO₂, NO_x, and soot and dust separately and increasing all three tax rates simultaneously. Our estimation results show that: raising tax rates of SO₂, NO_x, and soot and dust simultaneously would have the highest emission reduction effects, with the SO₂ tax rate making the greatest contribution to emission reduction. Raising the soot and dust tax rate is the most environmentally friendly strategy due to its highest abatement to welfare through avoided health costs. The combination of frontier analysis and input–output analysis provides policy makers a comprehensive and sectoral approach to assess costs and benefits of environmental taxation.     

Meta‐analysis on blood transcriptomic studies identifies consistently coexpressed protein–protein interaction modules as robust markers of human aging

The bodily decline that occurs with advancing age strongly impacts on the prospects for future health and life expectancy. Despite the profound role of age in disease etiology, knowledge about the molecular mechanisms driving the process of aging in humans is limited. Here, we used an integrative network-based approach for combining multiple large-scale expression studies in blood (2539 individuals) with protein–protein Interaction (PPI) data for the detection of consistently coexpressed PPI modules that may reflect key processes that change throughout the course of normative aging. Module detection followed by a meta-analysis on chronological age identified fifteen consistently coexpressed PPI modules associated with chronological age, including a highly significant module (P = 3.5 × 10⁻³⁸) enriched for ‘T-cell activation’ marking age-associated shifts in lymphocyte blood cell counts (R² = 0.603; P = 1.9 × 10⁻¹⁰). Adjusting the analysis in the compendium for the ‘T-cell activation’ module showed five consistently coexpressed PPI modules that robustly associated with chronological age and included modules enriched for ‘Translational elongation’, ‘Cytolysis’ and ‘DNA metabolic process’. In an independent study of 3535 individuals, four of five modules consistently associated with chronological age, underpinning the robustness of the approach. We found three of five modules to be significantly enriched with aging-related genes, as defined by the GenAge database, and association with prospective survival at high ages for one of the modules including ASF1A. The hereby-detected age-associated and consistently coexpressed PPI modules therefore may provide a molecular basis for future research into mechanisms underlying human aging.