Next‐generation sequencing workflow for assembly of nonmodel mitogenomes exemplified with North Pacific albatrosses (Phoebastria spp.)

Use of complete mitochondrial genomes (mitogenomes) can greatly increase the resolution achievable in phylogeographic and historical demographic studies. Using next‐generation sequencing methods, it is now feasible to efficiently sequence mitogenomes of large numbers of individuals once a reference mitogenome is available. However, assembling the initial mitogenomes of nonmodel organisms can present challenges, for example, in birds, where mtDNA is often subject to gene rearrangements and duplications. We developed a workflow based on Illumina paired‐end, whole‐genome shotgun sequencing, which we used to generate complete 19‐kilobase mitogenomes for each of three species of North Pacific albatross, a group of birds known to carry a tandem duplication. Although this duplication had been described previously, our procedure did not depend on this prior knowledge, nor did it require a closely related reference mitogenome (e.g. a mammalian mitogenome was sufficient). We employed an iterative process including de novo assembly, reference‐guided assembly and gap closing, which enabled us to detect duplications, determine gene order and identify sequence for primer positioning to resolve any mitogenome ambiguity (via minimal targeted Sanger sequencing). We present full mtDNA annotations, including 22 tRNAs, 2 rRNAs, 13 protein‐coding genes, a control region and a duplicated feature for all three species. Pairwise comparisons supported previous hypotheses regarding the phylogenetic relationships within this group and occurrence of a shared tandem duplication. The resulting mitogenome sequences will enable rapid, high‐throughput NGS mitogenome sequencing of North Pacific albatrosses via direct reference‐guided assembly. Moreover, our approach to assembling mitogenomes should be applicable to any taxon.     

A high‐quality genome of taro (Colocasia esculenta (L.) Schott), one of the world's oldest crops

Taro (Colocasia esculenta (L.), Schott), from the Araceae family, is one of the oldest crops with important edible, medicinal, nutritional and economic value. Taro is a highly polymorphic species including diverse genotypes adapted to a broad range of environments, but the taro genome has rarely been investigated. Here, a high‐quality chromosome‐level genome of C. esculenta was assembled using data sequenced by Illumina, PacBio and Nanopore platforms. The assembled genome size was 2,405 Mb with a contig N50 of 400.0 kb and a scaffold N50 of 159.4 Mb. In total, 2,311 Mb (96.09%) of the contig sequences was anchored onto 14 chromosomes to form pseudomolecules, and 2,126 Mb (88.43%) was annotated as repetitive sequences. Of the 28,695 predicted protein‐coding genes, 26,215 genes (91.4%) could be functionally annotated. On the basis of phylogenetic analysis using 769 genes, C. esculenta and Spirodela polyrhiza were placed on one branch of the tree that diverged approximately 73.23 million years ago. The synteny analyses showed that there have been two whole‐genome duplication events in C. esculenta separated by a relatively short gap. According to comparative genome analysis, a larger number (1,189) of distinct gene families and long terminal repeats were enriched in C. esculenta. Our high‐quality taro genome will provide valuable resources for further genetic, ecological and evolutionary analyses of taro or other species in the Araceae.     

Effects of invasion of fire‐free arid shrublands by a fire‐promoting invasive alien grass (Pennisetum setaceum) in South Africa

Arid shrublands in the Karoo (South Africa) seldom accumulate sufficient combustible fuel to support fire. However, as a result of invasion by an alien perennial grass (Pennisetum setaceum), they could become flammable. This paper reports on an experiment to assess the effects of fire following invasion by P. setaceum. We established 10 plots (5 × 10 m) separated by 2.5 m, and added grass fuel to five plots (5 and 10 tons ha^?1 to alternate halves of the plot) leaving the remaining five plots as interspersed controls. Plots with fuel added were burnt, and fire behaviour was measured during the burns. Rates of fire spread were generally low (0.01–0.07 m s^?1) and did not differ significantly between burn treatments. Mean fireline intensities were higher in the high compared with the low fuel treatments (894 and 427 kW m^?1, respectively). We recorded plant species and their cover before and after burning on each of the plots. After 15 months of follow‐up monitoring in the burn plots, only two species, the dwarf shrub (Tripteris sinuata) and the perennial herb (Gazania krebsiana) resprouted. Most individuals of other species were killed and did not reseed during the 15‐month study. The mass of added fuel load (high or low) did not influence vegetation recovery rates after fire. Should future invasions by P. setaceum lead to similar fuel loads in these shrublands, inevitable fires could change the vegetation and may favour spread of the flammable grass. Our results have important implications for predicting the effects of invasive alien plants (especially grasses) on fire‐free ecosystems elsewhere. The predicted impacts of fire may alter species composition, ultimately affecting core natural resources that support the Karoo economy.     

High‐Efficiency Thermoelectric Power Generation Enabled by Homogeneous Incorporation of MXene in (Bi,Sb)2Te3 Matrix

The (Bi,Sb)₂Te₃ (BST) compounds have long been considered as the benchmark of thermoelectric (TE) materials near room temperature especially for refrigeration. However, their unsatisfactory TE performances in wide‐temperature range severely restrict the large‐scale applications for power generation. Here, using a self‐assembly protocol to deliver a homogeneous dispersion of 2D inclusion in matrix, the first evidence is shown that incorporation of MXene (Ti₃C₂T_x) into BST can simultaneously achieve the improved power factor and greatly reduced thermal conductivity. The oxygen‐terminated Ti₃C₂T_x with proper work function leads to highly increased electrical conductivity via hole injection and retained Seebeck coefficient due to the energy barrier scattering. Meanwhile, the alignment of Ti₃C₂T_x with the layered structure significantly suppresses the phonon transport, resulting in higher interfacial thermal resistance. Accordingly, a peak ZT of up to 1.3 and an average ZT value of 1.23 from 300 to 475 K are realized for the 1 vol% Ti₃C₂T_x/BST composite. Combined with the high‐performance composite and rational device design, a record‐high thermoelectric conversion efficiency of up to 7.8% is obtained under a temperature gradient of 237 K. These findings provide a robust and scalable protocol to incorporate MXene as a versatile 2D inclusion for improving the overall performance of TE materials toward high energy‐conversion efficiency.     

Species‐level phylogeographical history of the endemic species Calligonum roborovskii and its close relatives in Calligonum section Medusa (Polygonaceae) in arid north‐western China

Quaternary climatic oscillations appear to have influenced the genetic diversity and evolutionary history of arid‐adapted plants. To understand the processes involved and reveal evolutionary relationships, haplotypes were examined from Calligonum roborovskii, an endemic species occurring in the arid zones across the desert regions of north‐western China, and seven other species also from Calligonum section Medusa, including C. gobicum, C. mongolicum and the narrow endemic species C. ebi‐nuricum, C. pumilum, C. taklimakanense, C. trifarium and C. yengisaricum. Forty‐three haplotypes were identified in 422 individuals from 51 natural populations, from variation of two plastid DNA intergenic spacers (rpl32‐trnL and ycf6‐psbM). A high level of total genetic diversity was found across species for which more than two populations were examined, including C. gobicum, C. mongolicum, C. pumilum and C. roborovskii. A distinct isolation‐by‐distance pattern in each of these species was suggested by the Mantel test, indicating that restricted gene flow caused high genetic differentiation among populations. Three haplotypes were shared by two or three species each, but the other 40 haplotypes were species‐specific. The 43 haplotypes split into three major clades, but not species‐specific lineages; most of the Calligonum species were not reciprocally monophyletic, probably due to incomplete lineage sorting or introgression. The identified haplotypes were dated to 1.97 Mya (95% highest posterior density: 2.95–0.99 Mya) and diverged until the late Pleistocene, possibly linked to aridification and enlargement of deserts caused by climate changes. Variation of desert habitats during the Pleistocene might play a key role in causing the divergence.     

Range expansion during the Pleistocene drove morphological radiation of the fir genus (Abies,Pinaceae) in the Qinghai‐Tibet Plateau and Himalayas

Range expansion caused by climate oscillations in the past probably promoted morphological radiation in a few plant groups. In this study, we aim to test this hypothesis through phylogeographical analysis of the cold‐tolerant fir genus (Abies) in the Qinghai‐Tibet Plateau (QTP) and Himalayas, where it comprises 12 described species. We examined sequence variation in two maternally inherited mitochondrial (mt) DNA fragments (nad5‐4 and nad7‐1) and two paternally inherited plastid DNA fragments (trnS‐G and trnL‐F) for 733 individuals from 75 populations of the species in a monophyletic group. Only six mtDNA haplotypes were recovered, but five were shared between multiple species and one occurred at a high frequency, providing strong evidence of range expansion. Forty‐three plastid DNA haplotypes were detected, 19 of which were shared between species and three occurred at high frequency. Network, mismatch and Bayesian skyline plot analyses of all plastid DNA haplotypes from this clade clearly suggested range expansion. This expansion was dated as having occurred during the longest and most extensive glaciation in the Pleistocene. Our results therefore supported the range expansion hypothesis for this clade of Abies during the Pleistocene; expansion probably drove the morphological radiation of the clade in the QTP and Himalayas, although it remains unclear whether the different morphotypes should be acknowledged as independent, reproductively isolated species. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 179 , 444–453.     

Mosaics in the mangroves: allopatric diversification of tree‐climbing mudwhelks (Gastropoda: Potamididae: Cerithidea) in the Indo‐West Pacific

The Indo‐Australian Archipelago (IAA) is the richest area of biodiversity in the marine realm, yet the processes that generate and maintain this diversity are poorly understood and have hardly been studied in the mangrove biotope. Cerithidea is a genus of marine and brackish‐water snails restricted to mangrove habitats in the Indo‐West Pacific, and its species are believed to have a short pelagic larval life. Using molecular and morphological techniques, we demonstrate the existence of 15 species, reconstruct their phylogeny and plot their geographical ranges. Sister species show a pattern of narrowly allopatric ranges across the IAA, with overlap only between clades that show evidence of ecological differentiation. These allopatric mosaic distributions suggest that speciation may have been driven by isolation during low sea‐level stands, during episodes preceding the Plio‐Pleistocene glaciations. The Makassar Strait forms a biogeographical barrier hindering eastward dispersal, corresponding to part of Wallace's Line in the terrestrial realm. Areas of maximum diversity of mangrove plants and their associated molluscs do not coincide closely. © 2013 The Natural History Museum. Biological Journal of the Linnean Society © 2013 The Linnean Society of London, 2013, 110 , 564–580.     

Crystal structure of a trimeric form of the KV7.1 (KCNQ1) A‐domain tail coiled‐coil reveals structural plasticity and context dependent changes in a putative coiled‐coil trimerization motif

Coiled-coils are widespread protein–protein interaction motifs typified by the heptad repeat (abcdefg)_n in which “a” and “d” positions are hydrophobic residues. Although identification of likely coiled-coil sequences is robust, prediction of strand order remains elusive. We present the X-ray crystal structure of a short form (residues 583–611), “Q1-short,” of the coiled-coil assembly specificity domain from the voltage-gated potassium channel Kv7.1 (KCNQ1) determined at 1.7 Å resolution. Q1-short lacks one and half heptads present in a previously studied tetrameric coiled-coil construct, Kv7.1 585–621, “Q1-long.” Surprisingly, Q1-short crystallizes as a trimer. In solution, Q1-short self-assembles more poorly than Q1-long and depends on an R-h-x-x-h-E motif common to trimeric coiled-coils. Addition of native sequences that include “a” and “d” positions C-terminal to Q1-short overrides the R-h-x-x-h-E motif influence and changes assembly state from a weakly associated trimer to a strongly associated tetramer. These data provide a striking example of a naturally occurring amino sequence that exhibits context-dependent folding into different oligomerization states, a three-stranded versus a four-stranded coiled-coil. The results emphasize the degenerate nature of coiled-coil energy landscapes in which small changes can have drastic effects on oligomerization. Discovery of these properties in an ion channel assembly domain and prevalence of the R-h-x-x-h-E motif in coiled-coil assembly domains of a number of different channels that are thought to function as tetrameric assemblies raises the possibility that such sequence features may be important for facilitating the assembly of intermediates en route to the final native state.     

Cross‐strand disulfides in the hydrogen bonding site of antiparallel β‐sheet (aCSDhs): Forbidden disulfides that are highly strained,easily broken

Some disulfide bonds perform important structural roles in proteins, but another group has functional roles via redox reactions. Forbidden disulfides are stressed disulfides found in recognizable protein contexts, which currently constitute more than 10% of all disulfides in the PDB. They likely have functional redox roles and constitute a major subset of all redox‐active disulfides. The torsional strain of forbidden disulfides is typically higher than for structural disulfides, but not so high as to render them immediately susceptible to reduction under physionormal conditions. Previously we characterized the most abundant forbidden disulfide in the Protein Data Bank, the aCSDn: a canonical motif in which disulfide‐bonded cysteine residues are positioned directly opposite each other on adjacent anti‐parallel β‐strands such that the backbone hydrogen‐bonded moieties are directed away from each other. Here we perform a similar analysis for the aCSDh, a less common motif in which the opposed cysteine residues are backbone hydrogen bonded. Oxidation of two Cys in this context places significant strain on the protein system, with the β‐chains tilting toward each other to allow disulfide formation. Only left‐handed aCSDh conformations are compatible with the inherent right‐handed twist of β‐sheets. aCSDhs tend to be more highly strained than aCSDns, particularly when both hydrogen bonds are formed. We discuss characterized roles of aCSDh motifs in proteins of the dataset, which include catalytic disulfides in ribonucleotide reductase and ahpC peroxidase as well as a redox‐active disulfide in P1 lysozyme, involved in a major conformation change. The dataset also includes many binding proteins.     

Phylogeography of the Indo‐West Pacific maskrays (Dasyatidae,Neotrygon): a complex example of chondrichthyan radiation in the Cenozoic

Maskrays of the genus Neotrygon (Dasyatidae) have dispersed widely in the Indo‐West Pacific being represented largely by an assemblage of narrow‐ranging coastal endemics. Phylogenetic reconstruction methods reproduced nearly identical and statistically robust topologies supporting the monophyly of the genus Neotrygon within the family Dasyatidae, the genus Taeniura being consistently basal to Neotrygon, and Dasyatis being polyphyletic to the genera Taeniurops and Pteroplatytrygon. The Neotrygon kuhlii complex, once considered to be an assemblage of color variants of the same biological species, is the most derived and widely dispersed subgroup of the genus. Mitochondrial (COI, 16S) and nuclear (RAG1) phylogenies used in synergy with molecular dating identified paleoclimatic fluctuations responsible for periods of vicariance and dispersal promoting population fragmentation and speciation in Neotrygon. Signatures of population differentiation exist in N. ningalooensis and N. annotata, yet a large‐scale geological event, such as the collision between the Australian and Eurasian Plates, coupled with subsequent sea‐level falls, appears to have separated a once homogeneous population of the ancestral form of N. kuhlii into southern Indian Ocean and northern Pacific taxa some 4–16 million years ago. Repeated climatic oscillations, and the subsequent establishment of land and shallow sea connections within and between Australia and parts of the Indo‐Malay Archipelago, have both promoted speciation and established zones of secondary contact within the Indian and Pacific Ocean basins.     

Transient formation of the oxo–iron(IV) porphyrin radical cation during the reaction of iron(III) tetrakis‐5,10,15,20‐(N‐methyl‐4‐pyridyl)porphyrin with hydrogen peroxide in aqueous solution

The reaction of iron(III) tetrakis‐5,10,15,20‐(N‐methyl‐4‐pyridyl)porphyrin (Fe(III)TMPyP) with hydrogen peroxide (H₂O₂) and the catalytic activity of the reaction intermediates on the luminescent peroxidation of luminol in aqueous solution were studied by using a double‐mixing stopped‐flow system. The observed luminescence intensities showed biphasic decay depending on the conditions. The initial flashlight decayed within <1 s followed by a sustained emission for more than 30 s. Computer deconvolution of the time‐resolved absorption spectra under the same conditions revealed that the initial flashlight appeared during the formation of the oxo–iron(IV) porphyrin, TMPyPFe(IV) = O, which is responsible for the sustained emission. The absorption spectra 0.0–0.5 s did not reproduce well by a simple combination of the two spectra of Fe(III)TMPyP and TMPyPFe(IV) = O, indicating that transient species was formed at the initial stage. Addition of uric acid (UA) caused a significant delay in the initiation of the luminol emission as well as in the formation of the TMPyPFe(IV) = O. Both of them were completely diminished in the presence of UA equimolar with H₂O₂, while mannitol had no effect at all. The delay of the light emission as well as the appearance of TMPyPFe(IV) = O was directly proportional to the [UA]₀ but other kinetic profiles were not changed significantly. Based on these observations and the kinetic analysis, we confirmed the involvement of the oxo–iron(IV) porphyrin radical cation, (TMPyP)^·+Fe(IV) = O, as an obligatory intermediate in the rate‐determining step of the overall reaction, Fe(III)TMPyP + H₂O₂ → TMPyPFe(IV) = O, with a rate constant of k = 4.3 × 10⁴/mol/L/s. The rate constants for the reaction between the (TMPyP)^·+Fe(IV) = O and luminol, and between the TMPyPFe(IV) = O and luminol were estimated to be 3.6 × 10⁶/mol/L/s and 1.31 × 10⁴/mol/L/s, respectively. Copyright © 2003 John Wiley & Sons, Ltd.     

Crossing the uncrossable: novel trans‐valley biogeographic patterns revealed in the genetic history of low‐dispersal mygalomorph spiders (Antrodiaetidae,Antrodiaetus) from California

Antrodiaetus riversi is a dispersal‐limited, habitat‐specialized mygalomorph spider species endemic to mesic woodlands of northern and central California. Here, we build upon prior phylogeographic research using a much larger geographic sample and include additional nuclear genes, providing more detailed biogeographic insights throughout the range of this complex. Of particular interest is the uncovering of unexpected and replicated trans‐valley biogeographic patterns, where in two separate genetic clades western haplotypes in the California south Coast Ranges are phylogenetically closely related to eastern haplotypes from central and northern Sierran foothills. In both instances, these trans‐valley phylogenetic patterns are strongly supported by multiple genes. These western and eastern populations are currently separated by the Central Valley, a well‐recognized modern‐day and historical biogeographic barrier in California. For one clade, the directionality is clearly northeast to southwest, and all available evidence is consistent with a jump dispersal event estimated at 1.2–1.3 Ma. During this time period, paleogeographic data indicate that northern Sierran rivers emptied to the ocean in the south Coast Ranges, rather than at the San Francisco Bay. For the other trans‐valley clade genetic evidence is less conclusive regarding the mechanism and directionality of biogeographic exchange, although the estimated timeframe is similar (approximately 1.8 Ma). Despite the large number of biogeographic studies previously conducted in central California, to the best of our knowledge no prior studies have discussed or revealed a northern Sierran to south Coast Range biogeographic connection. This uniqueness may reflect the low‐dispersal biology of mygalomorph spiders, where ‘post‐event’ gene exchange rarely erases historical biogeographic signal.     

The geography of speciation in narrow‐range endemics of the ‘Haenydra’ lineage (Coleoptera,Hydraenidae, Hydraena)

Aim We test whether species of western Mediterranean aquatic Coleoptera of the ‘Haenydra’ lineage (Hydraenidae, Hydraena) originated through: (1) successive periods of dispersal and speciation, (2) range fragmentation by random vicariance, or (3) range fragmentation by geographic isolation owing to a general reduction of population density. Location Europe. Methods To discriminate between scenarios we use contrasting predictions of the relationship between phylogenetic and geographic distance. The phylogeny was based on 3 kb of four mitochondrial and two nuclear gene fragments of about half of the known species of ‘Haenydra’, including most western Mediterranean taxa. Divergences were estimated using a molecular clock. The relationship between phylogenetic and geographic distance was tested using bivariate plots, Mantel tests and comparison of the observed phylogeny with the one minimizing geographic distances between species, as measured using Euclidean minimum spanning trees (EMSTs). Results The monophyly of ‘Haenydra’ was strongly supported, although its phylogenetic placement was not resolved. ‘Haenydra’ was estimated to be of late Miocene age, with most species originating during the Pleistocene. In two clades (Hydraena tatii and Hydraena emarginata clades) there was a significant association between geographic and phylogenetic distance, and the reconstructed phylogeny was identical to that obtained through the EMST, demonstrating a strong non‐randomness of the geographic distribution of the species. In two other clades (Hydraena iberica and Hydraena bitruncata clades) there was no association between geographic and phylogenetic distance, and the observed phylogeny was not the one minimizing geographic distances. In one of the clades this seems to be due to a secondary, recent range expansion of one species (H. iberica), which erased the geographic signal of their distributions. Main conclusions We show that it is possible to obtain strong evidence of stasis of the geographic ranges of narrow‐range endemic species through the study of their phylogenetic relationships and current distributions. In at least two of the studied clades, current species seem to have originated through the fragmentation of a more widely distributed species, without further range movements. A process of range expansion and fragmentation may have occurred repeatedly within the ‘Haenydra’ lineage, contributing to the accumulation of narrow‐range endemics in Mediterranean Pleistocene refugia.     

Validated spectrofluorimetric method for the determination of carbamazepine in pharmaceutical dosage forms after reaction with 4‐chloro‐7‐‐nitrobenzo‐2‐oxa‐1,3‐diazole (NBD‐Cl)

A sensitive and simple spectrofluorimetric method has been developed and validated for the determination of the anti‐epileptic drug carbamazepine (CBZ) in its dosage forms. The method was based on a nucleophilic substitution reaction of CBZ with 4‐chloro‐7‐nitrobenzo‐2‐ oxa‐1,3‐diazole (NBD‐Cl) in borate buffer (pH 9) to form a highly fluorescent derivative that was measured at 530 nm after excitation at 460 nm. Factors affecting the formation of the reaction product were studied and optimized, and the reaction mechanism was postulated. The fluorescence–concentration plot is rectilinear over the range of 0.6–8 µg/mL with limit of detection of 0.06 µg/mL and limit of quantitation of 0.19 µg/mL. The method was applied to the analysis of commercial tablets and the results were in good agreement with those obtained using the reference method. Validation of the analytical procedures was evaluated according to ICH guidelines. Copyright © 2015 John Wiley & Sons, Ltd.     

Testing and estimation of proportion (or risk) ratio under the matched‐pair design with multiple binary endpoints

The proportion ratio (PR) of responses between an experimental treatment and a control treatment is one of the most commonly used indices to measure the relative treatment effect in a randomized clinical trial. We develop asymptotic and permutation‐based procedures for testing equality of treatment effects as well as derive confidence intervals of PRs for multivariate binary matched‐pair data under a mixed‐effects exponential risk model. To evaluate and compare the performance of these test procedures and interval estimators, we employ Monte Carlo simulation. When the number of matched pairs is large, we find that all test procedures presented here can perform well with respect to Type I error. When the number of matched pairs is small, the permutation‐based test procedures developed in this paper is of use. Furthermore, using test procedures (or interval estimators) based on a weighted linear average estimator of treatment effects can improve power (or gain precision) when the treatment effects on all response variables of interest are known to fall in the same direction. Finally, we apply the data taken from a crossover clinical trial that monitored several adverse events of an antidepressive drug to illustrate the practical use of test procedures and interval estimators considered here.     

Diploid hybrid origin of Ostryopsis intermedia (Betulaceae) in the Qinghai‐Tibet Plateau triggered by Quaternary climate change

Despite the well‐known effects that Quaternary climate oscillations had on shaping intraspecific diversity, their role in driving homoploid hybrid speciation is less clear. Here, we examine their importance in the putative homoploid hybrid origin and evolution of Ostryopsis intermedia, a diploid species occurring in the Qinghai‐Tibet Plateau (QTP), a biodiversity hotspot. We investigated interspecific relationships between this species and its only other congeners, O. davidiana and O. nobilis, based on four sets of nuclear and chloroplast population genetic data and tested alternative speciation hypotheses. All nuclear data distinguished the three species clearly and supported a close relationship between O. intermedia and the disjunctly distributed O. davidiana. Chloroplast DNA sequence variation identified two tentative lineages, which distinguished O. intermedia from O. davidiana; however, both were present in O. nobilis. Admixture analyses of genetic polymorphisms at 20 SSR loci and sequence variation at 11 nuclear loci and approximate Bayesian computation (ABC) tests supported the hypothesis that O. intermedia originated by homoploid hybrid speciation from O. davidiana and O. nobilis. We further estimated that O. davidiana and O. nobilis diverged 6–11 Ma, while O. intermedia originated 0.5–1.2 Ma when O. davidiana is believed to have migrated southward, contacted and hybridized with O. nobilis possibly during the largest Quaternary glaciation that occurred in this region. Our findings highlight the importance of Quaternary climate change in the QTP in causing hybrid speciation in this important biodiversity hotspot.