Inherent differences of small airway contraction and Ca2+ oscillations in airway smooth muscle cells between BALB/c and C57BL/6 mouse strains.

BALB/c and C57BL/6 mouse strains are widely used as animal model in studies of respiratory diseases, such as asthma. Asthma is characterized by airway hyperresponsiveness, which is eventually resulted from the excessive airway smooth muscle (ASM) contraction mediated by Ca2+ oscillations in ASM cells. It is reported that BALB/c mice have inherently higher airway responsiveness, but show no different contractive response of tracheal ring as compared to C57BL/6 mice. However, whether the different airway responsiveness is due to the different extents of small airway contraction, and what's underlying mechanism remains unknown. Here, we assess agonist-induced small airway contraction and Ca2+ oscillations in ASM cells between BALB/c and C57BL/6 mice by using precision-cut lung slices (PCLS). We found that BALB/c mice showed an intrinsically stronger extent of small airway narrowing and faster Ca2+ oscillations in ASM cells in response to agonists. These differences were associated with a higher magnitude of Ca2+ influx via store-operated Ca2+ entry (SOCE), as a result of increased expression of SOCE components (STIM1, Orai1) in the ASM cells of small airway of BALB/c mice. An established mathematical model and experimental results suggested that the increased SOC current could result in increased agonist-induced Ca2+ oscillations. Therefore, the inherently higher SOC underlies the increased Ca2+ oscillation frequency in ASM cells and stronger small airway contraction in BALB/c mice, thus higher airway responsiveness in BALB/c than C57BL/6 mouse strain.

Origin and Diversification of the Saguaro Cactus (Carnegiea gigantea): A Within-Species Phylogenomic Analysis.

Reconstructing accurate historical relationships within a species poses numerous challenges, not least in many plant groups in which gene flow is high enough to extend well beyond species boundaries. Nonetheless, the extent of tree-like history within a species is an empirical question on which it is now possible to bring large amounts of genome sequence to bear. We assess phylogenetic structure across the geographic range of the saguaro cactus, an emblematic member of Cactaceae, a clade known for extensive hybridization and porous species boundaries. Using 200 Gb of whole genome resequencing data from 20 individuals sampled from 10 localities, we assembled two data sets comprising 150,000 biallelic single nucleotide polymorphisms (SNPs) from protein coding sequences. From these, we inferred within-species trees and evaluated their significance and robustness using five qualitatively different inference methods. Despite the low sequence diversity, large census population sizes, and presence of wide-ranging pollen and seed dispersal agents, phylogenetic trees were well resolved and highly consistent across both data sets and all methods. We inferred that the most likely root, based on marginal likelihood comparisons, is to the east and south of the region of highest genetic diversity, which lies along the coast of the Gulf of California in Sonora, Mexico. Together with striking decreases in marginal likelihood found to the north, this supports hypotheses that saguaro's current range reflects postglacial expansion from the refugia in the south of its range. We conclude with observations about practical and theoretical issues raised by phylogenomic data sets within species, in which SNP-based methods must be used rather than gene tree methods that are widely used when sequence divergence is higher. These include computational scalability, inference of gene flow, and proper assessment of statistical support in the presence of linkage effects. [Phylogenomics; phylogeography; rooting; Sonoran Desert.].

Accuracy in Near-Perfect Virus Phylogenies.

Phylogenetic trees from real-world data often include short edges with very few substitutions per site, which can lead to partially resolved trees and poor accuracy. Theory indicates that the number of sites needed to accurately reconstruct a fully resolved tree grows at a rate proportional to the inverse square of the length of the shortest edge. However, when inferred trees are partially resolved due to short edges, "accuracy" should be defined as the rate of discovering false splits (clades on a rooted tree) relative to the actual number found. Thus, accuracy can be high even if short edges are common. Specifically, in a "near-perfect" parameter space in which trees are large, the tree length $\xi$ (the sum of all edge lengths) is small, and rate variation is minimal, the expected false positive rate is less than $\xi/3$; the exact value depends on tree shape and sequence length. This expected false positive rate is far below the false negative rate for small $\xi$ and often well below 5% even when some assumptions are relaxed. We show this result analytically for maximum parsimony and explore its extension to maximum likelihood using theory and simulations. For hypothesis testing, we show that measures of split "support" that rely on bootstrap resampling consistently imply weaker support than that implied by the false positive rates in near-perfect trees. The near-perfect parameter space closely fits several empirical studies of human virus diversification during outbreaks and epidemics, including Ebolavirus, Zika virus, and SARS-CoV-2, reflecting low substitution rates relative to high transmission/sampling rates in these viruses.[Ebolavirus; epidemic; HIV; homoplasy; mumps virus; perfect phylogeny; SARS-CoV-2; virus; West Nile virus; Yule-Harding model; Zika virus.].

Plastome Structural Evolution and Homoplastic Inversions in Neo-Astragalus (Fabaceae).

The plastid genomes of photosynthetic green plants have largely maintained conserved gene content and order as well as structure over hundreds of millions of years of evolution. Several plant lineages, however, have departed from this conservation and contain many plastome structural rearrangements, which have been associated with an abundance of repeated sequences both overall and near rearrangement endpoints. We sequenced the plastomes of 25 taxa of Astragalus L. (Fabaceae), a large genus in the inverted repeat-lacking clade of legumes, to gain a greater understanding of the connection between repeats and plastome inversions. We found plastome repeat structure has a strong phylogenetic signal among these closely related taxa mostly in the New World clade of Astragalus called Neo-Astragalus. Taxa without inversions also do not differ substantially in their overall repeat structure from four taxa each with one large-scale inversion. For two taxa with inversion endpoints between the same pairs of genes, differences in their exact endpoints indicate the inversions occurred independently. Our proposed mechanism for inversion formation suggests the short inverted repeats now found near the endpoints of the four inversions may be there as a result of these inversions rather than their cause. The longer inverted repeats now near endpoints may have allowed the inversions first mediated by shorter microhomologous sequences to propagate, something that should be considered in explaining how any plastome rearrangement becomes fixed regardless of the mechanism of initial formation.

Evolution of l-DOPA 4,5-dioxygenase activity allows for recurrent specialisation to betalain pigmentation in Caryophyllales.

The evolution of l-DOPA 4,5-dioxygenase activity, encoded by the gene DODA, was a key step in the origin of betalain biosynthesis in Caryophyllales. We previously proposed that l-DOPA 4,5-dioxygenase activity evolved via a single Caryophyllales-specific neofunctionalisation event within the DODA gene lineage. However, this neofunctionalisation event has not been confirmed and the DODA gene lineage exhibits numerous gene duplication events, whose evolutionary significance is unclear. To address this, we functionally characterised 23 distinct DODA proteins for l-DOPA 4,5-dioxygenase activity, from four betalain-pigmented and five anthocyanin-pigmented species, representing key evolutionary transitions across Caryophyllales. By mapping these functional data to an updated DODA phylogeny, we then explored the evolution of l-DOPA 4,5-dioxygenase activity. We find that low l-DOPA 4,5-dioxygenase activity is distributed across the DODA gene lineage. In this context, repeated gene duplication events within the DODA gene lineage give rise to polyphyletic occurrences of elevated l-DOPA 4,5-dioxygenase activity, accompanied by convergent shifts in key functional residues and distinct genomic patterns of micro-synteny. In the context of an updated organismal phylogeny and newly inferred pigment reconstructions, we argue that repeated convergent acquisition of elevated l-DOPA 4,5-dioxygenase activity is consistent with recurrent specialisation to betalain synthesis in Caryophyllales.

The Implications of Lineage-Specific Rates for Divergence Time Estimation.

Rate variation adds considerable complexity to divergence time estimation in molecular phylogenies. Here, we evaluate the impact of lineage-specific rates-which we define as among-branch-rate-variation that acts consistently across the entire genome. We compare its impact to residual rates-defined as among-branch-rate-variation that shows a different pattern of rate variation at each sampled locus, and gene-specific rates-defined as variation in the average rate across all branches at each sampled locus. We show that lineage-specific rates lead to erroneous divergence time estimates, regardless of how many loci are sampled. Further, we show that stronger lineage-specific rates lead to increasing error. This contrasts to residual rates and gene-specific rates, where sampling more loci significantly reduces error. If divergence times are inferred in a Bayesian framework, we highlight that error caused by lineage-specific rates significantly reduces the probability that the 95% highest posterior density includes the correct value, and leads to sensitivity to the prior. Use of a more complex rate prior-which has recently been proposed to model rate variation more accurately-does not affect these conclusions. Finally, we show that the scale of lineage-specific rates used in our simulation experiments is comparable to that of an empirical data set for the angiosperm genus Ipomoea. Taken together, our findings demonstrate that lineage-specific rates cause error in divergence time estimates, and that this error is not overcome by analyzing genomic scale multilocus data sets. [Divergence time estimation; error; rate variation.].

The flagellar protein Enkurin is required for mouse sperm motility and for transport through the female reproductive tract.

Enkurin was identified initially in mouse sperm where it was suggested to act as an intracellular adaptor protein linking membrane calcium influx to intracellular signaling pathways. In order to examine the function of this protein, a targeted mutation was introduced into the mouse Enkurin gene. Males that were homozygous for this mutated allele were subfertile. This was associated with lower rates of sperm transport in the female reproductive tract, including reduced entry into the oviduct and slower migration to the site of fertilization in the distal oviduct, and with poor progressive motility in vitro. Flagella from wild-type animals exhibited symmetrical bending and progressive motility in culture medium, and demembranated flagella exhibited the "curlicue" response to Ca2+ in vitro. In contrast, flagella of mice homozygous for the mutated allele displayed only asymmetric bending, nonprogressive motility, and a loss of Ca2+-responsiveness following demembrantion. We propose that Enkurin is part of a flagellar Ca2+-sensor that regulates bending and that the motility defects following mutation of the locus are the proximate cause of subfertility.

The prevalence of terraced treescapes in analyses of phylogenetic data sets.

BACKGROUND: The pattern of data availability in a phylogenetic data set may lead to the formation of terraces, collections of equally optimal trees. Terraces can arise in tree space if trees are scored with parsimony or with partitioned, edge-unlinked maximum likelihood. Theory predicts that terraces can be large, but their prevalence in contemporary data sets has never been surveyed. We selected 26 data sets and phylogenetic trees reported in recent literature and investigated the terraces to which the trees would belong, under a common set of inference assumptions. We examined terrace size as a function of the sampling properties of the data sets, including taxon coverage density (the proportion of taxon-by-gene positions with any data present) and a measure of gene sampling "sufficiency". We evaluated each data set in relation to the theoretical minimum gene sampling depth needed to reduce terrace size to a single tree, and explored the impact of the terraces found in replicate trees in bootstrap methods. RESULTS: Terraces were identified in nearly all data sets with taxon coverage densities < 0.90. They were not found, however, in high-coverage-density (i.e., ≥ 0.94) transcriptomic and genomic data sets. The terraces could be very large, and size varied inversely with taxon coverage density and with gene sampling sufficiency. Few data sets achieved a theoretical minimum gene sampling depth needed to reduce terrace size to a single tree. Terraces found during bootstrap resampling reduced overall support. CONCLUSIONS: If certain inference assumptions apply, trees estimated from empirical data sets often belong to large terraces of equally optimal trees. Terrace size correlates to data set sampling properties. Data sets seldom include enough genes to reduce terrace size to one tree. When bootstrap replicate trees lie on a terrace, statistical support for phylogenetic hypotheses may be reduced. Although some of the published analyses surveyed were conducted with edge-linked inference models (which do not induce terraces), unlinked models have been used and advocated. The present study describes the potential impact of that inference assumption on phylogenetic inference in the context of the kinds of multigene data sets now widely assembled for large-scale tree construction.

Orosomucoid-like 3 (ORMDL3) upregulates airway smooth muscle proliferation, contraction, and Ca2+ oscillations in asthma.

BACKGROUND: Airway hyperresponsiveness is a major feature of asthma attributed predominantly to an extrinsic immune/inflammatory response increasing airway smooth muscle (ASM) contractility. OBJECTIVE: We investigated whether increased ASM expression of orosomucoid-like 3 (ORMDL3), a gene on chromosome 17q21 highly linked to asthma, induced increased ASM proliferation and contractility in vitro and influenced airway contractility and calcium flux in ASM in precision-cut lung slices (PCLSs) from wild-type and hORMDL3Zp3-Cre mice (which express increased levels of human ORMDL3 [hORMDL3]). METHODS: Levels of ASM proliferation and contraction were assessed in ASM cells transfected with ORMDL3 in vitro. In addition, airway contractility and calcium oscillations were quantitated in ASM cells in PCLSs derived from naive wild-type and naive hORMDL3Zp3-Cre mice, which do not have a blood supply. RESULTS: Increased ASM expression of ORMDL3 in vitro resulted in increased ASM proliferation and contractility. PCLSs derived from naive hORMDL3Zp3-Cre mice, which do not have airway inflammation, exhibit increased airway contractility with increased calcium oscillations in ASM cells. Increased ASM ORMDL3 expression increases levels of ASM sarcoplasmic reticulum Ca2+ ATPase 2b (SERCA2b), which increases ASM proliferation and contractility. CONCLUSION: Overall, these studies provide evidence that an intrinsic increase in ORMDL3 expression in ASM can induce increased ASM proliferation and contractility, which might contribute to increased airway hyperresponsiveness in the absence of airway inflammation in asthmatic patients.

Extensive gene tree discordance and hemiplasy shaped the genomes of North American columnar cacti.

Few clades of plants have proven as difficult to classify as cacti. One explanation may be an unusually high level of convergent and parallel evolution (homoplasy). To evaluate support for this phylogenetic hypothesis at the molecular level, we sequenced the genomes of four cacti in the especially problematic tribe Pachycereeae, which contains most of the large columnar cacti of Mexico and adjacent areas, including the iconic saguaro cactus (Carnegiea gigantea) of the Sonoran Desert. We assembled a high-coverage draft genome for saguaro and lower coverage genomes for three other genera of tribe Pachycereeae (Pachycereus, Lophocereus, and Stenocereus) and a more distant outgroup cactus, Pereskia We used these to construct 4,436 orthologous gene alignments. Species tree inference consistently returned the same phylogeny, but gene tree discordance was high: 37% of gene trees having at least 90% bootstrap support conflicted with the species tree. Evidently, discordance is a product of long generation times and moderately large effective population sizes, leading to extensive incomplete lineage sorting (ILS). In the best supported gene trees, 58% of apparent homoplasy at amino sites in the species tree is due to gene tree-species tree discordance rather than parallel substitutions in the gene trees themselves, a phenomenon termed "hemiplasy." The high rate of genomic hemiplasy may contribute to apparent parallelisms in phenotypic traits, which could confound understanding of species relationships and character evolution in cacti.

Role of transient receptor potential vanilloid 1 in the modulation of airway smooth muscle tone and calcium handling.

Asthma is a common disorder characterized, in part, by airway smooth muscle (ASM) hyperresponsiveness. Transient receptor potential vanilloid 1 (TRPV1) is a nonselective cation channel expressed on airway nerve fibers that modulates afferent signals, resulting in cough, and potentially bronchoconstriction. In the present study, the TRPV1 transcript was detected by RT-PCR in primary cultured human ASM cells, and the TRPV1 protein was detected in ASM of human trachea by immunohistochemistry. Proximity ligation assays suggest that TRPV1 is expressed in the sarcoplasmic reticulum membrane of human ASM cells in close association with sarco/endoplasmic reticulum Ca2+-ATPase-2. In guinea pig tracheal ring organ bath experiments, the TRPV1 agonist capsaicin led to ASM contraction, but this contraction was significantly attenuated by the sodium channel inhibitor bupivacaine (n = 4, P < 0.05) and the neurokinin-2 receptor antagonist GR-159897 (n = 4, P < 0.05), suggesting that this contraction is neutrally mediated. However, pretreatment of guinea pig and human ASM in organ bath experiments with the TRPV1 antagonist capsazepine inhibited the maintenance phase of an acetylcholine-induced contraction (n = 4, P < 0.01 for both species). Similarly, capsazepine inhibited methacholine-induced contraction of peripheral airways in mouse precision-cut lung slice (PCLS) experiments (n = 4-5, P < 0.05). Although capsazepine did not inhibit store-operated calcium entry in mouse ASM cells in PCLS (n = 4-7, P = nonsignificant), it did inhibit calcium oscillations (n = 3, P < 0.001). These studies suggest that TRPV1 is expressed on ASM, including the SR, but that ASM TRPV1 activation does not play a significant role in initiation of ASM contraction. However, capsazepine does inhibit maintenance of contraction, likely by inhibiting calcium oscillations.

On the dynamical structure of calcium oscillations.

Oscillations in the concentration of free cytosolic Ca2+ are an important and ubiquitous control mechanism in many cell types. It is thus correspondingly important to understand the mechanisms that underlie the control of these oscillations and how their period is determined. We show that Class I Ca2+ oscillations (i.e., oscillations that can occur at a constant concentration of inositol trisphosphate) have a common dynamical structure, irrespective of the oscillation period. This commonality allows the construction of a simple canonical model that incorporates this underlying dynamical behavior. Predictions from the model are tested, and confirmed, in three different cell types, with oscillation periods ranging over an order of magnitude. The model also predicts that Ca2+ oscillation period can be controlled by modulation of the rate of activation by Ca2+ of the inositol trisphosphate receptor. Preliminary experimental evidence consistent with this hypothesis is presented. Our canonical model has a structure similar to, but not identical to, the classic FitzHugh-Nagumo model. The characterization of variables by speed of evolution, as either fast or slow variables, changes over the course of a typical oscillation, leading to a model without globally defined fast and slow variables.

Store-operated calcium entry is required for sustained contraction and Ca2+ oscillations of airway smooth muscle.

KEY POINTS: Airway hyper-responsiveness in asthma is driven by excessive contraction of airway smooth muscle cells (ASMCs). Agonist-induced Ca2+ oscillations underlie this contraction of ASMCs and the magnitude of this contraction is proportional to the Ca2+ oscillation frequency. Sustained contraction and Ca2+ oscillations require an influx of extracellular Ca2+ , although the mechanisms and pathways mediating this Ca2+ influx during agonist-induced ASMC contraction are not well defined. By inhibiting store-operated calcium entry (SOCE) or voltage-gated Ca2+ channels (VGCCs), we show that SOCE, rather than Ca2+ influx via VGCCs, provides the major Ca2+ entry pathway into ASMCs to sustain ASMCs contraction and Ca2+ oscillations. SOCE may therefore serve as a potential target for new bronchodilators to reduce airway hyper-responsiveness in asthma. ABSTRACT: Asthma is characterized by airway hyper-responsiveness: the excessive contraction of airway smooth muscle. The extent of this airway contraction is proportional to the frequency of Ca2+ oscillations within airway smooth muscle cells (ASMCs). Sustained Ca2+ oscillations require a Ca2+ influx to replenish Ca2+ losses across the plasma membrane. Our previous studies implied store-operated calcium entry (SOCE) as the major pathway for this Ca2+ influx. In the present study, we explore this hypothesis, by examining the effects of SOCE inhibitors (GSK7975A and GSK5498A) as well as L-type voltage-gated Ca2+ channel inhibitors (nifedipine and nimodipine) on airway contraction and Ca2+ oscillations and SOCE-mediated Ca2+ influx in ASMCs within mouse precision-cut lung slices. We found that both GSK7975A and GSK5498A were able to fully relax methacholine-induced airway contraction by abolishing the Ca2+ oscillations, in a manner similar to that observed in zero extracellular Ca2+ ([Ca2+ ]e ). In addition, GSK7975A and GSK5498A inhibited increases in intracellular Ca2+ ([Ca2+ ]i ) in ASMCs with depleted Ca2+ -stores in response to increased [Ca2+ ]e , demonstrating a response consistent with the inhibition of SOCE. However, GSK7975A and GSK5498A did not reduce Ca2+ release via IP3 receptors stimulated with IP3 released from caged-IP3 . By contrast, nifedipine and nimodipine only partially reduced airway contraction, Ca2+ oscillation frequency and SOCE-mediated Ca2+ influx. These data suggest that SOCE is the major Ca2+ influx pathway for ASMCs with respect to sustaining agonist-induced airway contraction and the underlying Ca2+ oscillations. The mechanisms of SOCE may therefore form novel targets for new bronchodilators.

Toward a Self-Updating Platform for Estimating Rates of Speciation and Migration, Ages, and Relationships of Taxa.

Rapidly growing biological data-including molecular sequences and fossils-hold an unprecedented potential to reveal how evolutionary processes generate and maintain biodiversity. However, researchers often have to develop their own idiosyncratic workflows to integrate and analyze these data for reconstructing time-calibrated phylogenies. In addition, divergence times estimated under different methods and assumptions, and based on data of various quality and reliability, should not be combined without proper correction. Here we introduce a modular framework termed SUPERSMART (Self-Updating Platform for Estimating Rates of Speciation and Migration, Ages, and Relationships of Taxa), and provide a proof of concept for dealing with the moving targets of evolutionary and biogeographical research. This framework assembles comprehensive data sets of molecular and fossil data for any taxa and infers dated phylogenies using robust species tree methods, also allowing for the inclusion of genomic data produced through next-generation sequencing techniques. We exemplify the application of our method by presenting phylogenetic and dating analyses for the mammal order Primates and for the plant family Arecaceae (palms). We believe that this framework will provide a valuable tool for a wide range of hypothesis-driven research questions in systematics, biogeography, and evolution. SUPERSMART will also accelerate the inference of a "Dated Tree of Life" where all node ages are directly comparable. [Bayesian phylogenetics; data mining; divide-and-conquer methods; GenBank; multilocus multispecies coalescent; next-generation sequencing; palms; primates; tree calibration.].

The relative contributions of store-operated and voltage-gated Ca2+ channels to the control of Ca2+ oscillations in airway smooth muscle.

KEY POINTS: Agonist-dependent oscillations in the concentration of free cytosolic calcium are a vital mechanism for the control of airway smooth muscle contraction and thus are a critical factor in airway hyper-responsiveness. Using a mathematical model, closely tied to experimental work, we show that the oscillations in membrane potential accompanying the calcium oscillations have no significant effect on the properties of the calcium oscillations. In addition, the model shows that calcium entry through store-operated calcium channels is critical for calcium oscillations, but calcium entry through voltage-gated channels has much less effect. The model predicts that voltage-gated channels are less important than store-operated channels in the control of airway smooth muscle tone. ABSTRACT: Airway smooth muscle contraction is typically the key mechanism underlying airway hyper-responsiveness, and the strength of muscle contraction is determined by the frequency of oscillations of intracellular calcium (Ca2+ ) concentration. In airway smooth muscle cells, these Ca2+ oscillations are caused by cyclic Ca2+ release from the sarcoplasmic reticulum, although Ca2+ influx via plasma membrane channels is also necessary to sustain the oscillations over longer times. To assess the relative contributions of store-operated and voltage-gated Ca2+ channels to this Ca2+ influx, we generated a comprehensive mathematical model, based on experimental Ca2+ measurements in mouse precision-cut lung slices, to simulate Ca2+ oscillations and changes in membrane potential. Agonist-induced Ca2+ oscillations are accompanied by oscillations in membrane potential, although the membrane potential oscillations are too small to generate large Ca2+ currents through voltage-gated Ca2+ channels, and thus have little effect on the Ca2+ oscillations. Ca2+ entry through voltage-gated channels only becomes important when the cell is depolarized (e.g. by a high external K+ concentration). As a result, agonist-induced Ca2+ oscillations are critically dependent on Ca2+ entry through store-operated channels but do not depend strongly on Ca2+ entry though voltage-gated channels.

Cryopreserved Human Precision-Cut Lung Slices as a Bioassay for Live Tissue Banking. A Viability Study of Bronchodilation with Bitter-Taste Receptor Agonists.

Human precision-cut lung slices (hPCLSs) provide a unique ex vivo model for translational research. However, the limited and unpredictable availability of human lung tissue greatly impedes their use. Here, we demonstrate that cryopreservation of hPCLSs facilitates banking of live human lung tissue for routine use. Our results show that cryopreservation had little effect on overall cell viability and vital functions of immune cells, including phagocytes and T lymphocytes. In addition, airway contraction and relaxation in response to specific agonists and antagonists, respectively, were unchanged after cryopreservation. At the subcellular level, cryopreserved hPCLSs maintained Ca(2+)-dependent regulatory mechanisms for the control of airway smooth muscle cell contractility. To exemplify the use of cryopreserved hPCLSs in smooth muscle research, we provide evidence that bitter-taste receptor (TAS2R) agonists relax airways by blocking Ca(2+) oscillations in airway smooth muscle cells. In conclusion, the banking of cryopreserved hPCLSs provides a robust bioassay for translational research of lung physiology and disease.

Tradeoff between robustness and elaboration in carotenoid networks produces cycles of avian color diversification.

BACKGROUND: Resolution of the link between micro- and macroevolution calls for comparing both processes on the same deterministic landscape, such as genomic, metabolic or fitness networks. We apply this perspective to the evolution of carotenoid pigmentation that produces spectacular diversity in avian colors and show that basic structural properties of the underlying carotenoid metabolic network are reflected in global patterns of elaboration and diversification in color displays. Birds color themselves by consuming and metabolizing several dietary carotenoids from the environment. Such fundamental dependency on the most upstream external compounds should intrinsically constrain sustained evolutionary elongation of multi-step metabolic pathways needed for color elaboration unless the metabolic network gains robustness - the ability to synthesize the same carotenoid from an additional dietary starting point. RESULTS: We found that gains and losses of metabolic robustness were associated with evolutionary cycles of elaboration and stasis in expressed carotenoids in birds. Lack of metabolic robustness constrained lineage's metabolic explorations to the immediate biochemical vicinity of their ecologically distinct dietary carotenoids, whereas gains of robustness repeatedly resulted in sustained elongation of metabolic pathways on evolutionary time scales and corresponding color elaboration. CONCLUSIONS: The structural link between length and robustness in metabolic pathways may explain periodic convergence of phylogenetically distant and ecologically distinct species in expressed carotenoid pigmentation; account for stasis in carotenoid colors in some ecological lineages; and show how the connectivity of the underlying metabolic network provides a mechanistic link between microevolutionary elaboration and macroevolutionary diversification.

Exceptional reduction of the plastid genome of saguaro cactus (Carnegiea gigantea): Loss of the ndh gene suite and inverted repeat.

UNLABELLED: • PREMISE OF THE STUDY: Land-plant plastid genomes have only rarely undergone significant changes in gene content and order. Thus, discovery of additional examples adds power to tests for causes of such genome-scale structural changes.• METHODS: Using next-generation sequence data, we assembled the plastid genome of saguaro cactus and probed the nuclear genome for transferred plastid genes and functionally related nuclear genes. We combined these results with available data across Cactaceae and seed plants more broadly to infer the history of gene loss and to assess the strength of phylogenetic association between gene loss and loss of the inverted repeat (IR).• KEY RESULTS: The saguaro plastid genome is the smallest known for an obligately photosynthetic angiosperm (∼113 kb), having lost the IR and plastid ndh genes. This loss supports a statistically strong association across seed plants between the loss of ndh genes and the loss of the IR. Many nonplastid copies of plastid ndh genes were found in the nuclear genome, but none had intact reading frames; nor did three related nuclear-encoded subunits. However, nuclear pgr5, which functions in a partially redundant pathway, was intact.• CONCLUSIONS: The existence of an alternative pathway redundant with the function of the plastid NADH dehydrogenase-like complex (NDH) complex may permit loss of the plastid ndh gene suite in photoautotrophs like saguaro. Loss of these genes may be a recurring mechanism for overall plastid genome size reduction, especially in combination with loss of the IR.

DRC3 connects the N-DRC to dynein g to regulate flagellar waveform.

The nexin-dynein regulatory complex (N-DRC), which is a major hub for the control of flagellar motility, contains at least 11 different subunits. A major challenge is to determine the location and function of each of these subunits within the N-DRC. We characterized a Chlamydomonas mutant defective in the N-DRC subunit DRC3. Of the known N-DRC subunits, the drc3 mutant is missing only DRC3. Like other N-DRC mutants, the drc3 mutant has a defect in flagellar motility. However, in contrast to other mutations affecting the N-DRC, drc3 does not suppress flagellar paralysis caused by loss of radial spokes. Cryo-electron tomography revealed that the drc3 mutant lacks a portion of the N-DRC linker domain, including the L1 protrusion, part of the distal lobe, and the connection between these two structures, thus localizing DRC3 to this part of the N-DRC. This and additional considerations enable us to assign DRC3 to the L1 protrusion. Because the L1 protrusion is the only non-dynein structure in contact with the dynein g motor domain in wild-type axonemes and this is the only N-DRC-dynein connection missing in the drc3 mutant, we conclude that DRC3 interacts with dynein g to regulate flagellar waveform.