Pluripotency of a founding field: rebranding developmental biology.

The field of developmental biology has declined in prominence in recent decades, with off-shoots from the field becoming more fashionable and highly funded. This has created inequity in discovery and opportunity, partly due to the perception that the field is antiquated or not cutting edge. A 'think tank' of scientists from multiple developmental biology-related disciplines came together to define specific challenges in the field that may have inhibited innovation, and to provide tangible solutions to some of the issues facing developmental biology. The community suggestions include a call to the community to help 'rebrand' the field, alongside proposals for additional funding apparatuses, frameworks for interdisciplinary innovative collaborations, pedagogical access, improved science communication, increased diversity and inclusion, and equity of resources to provide maximal impact to the community.

A long noncoding RNA at the cortex locus controls adaptive coloration in butterflies.

Evolutionary variation in the wing pigmentation of butterflies and moths offers striking examples of adaptation by crypsis and mimicry. The cortex locus has been independently mapped as the locus controlling color polymorphisms in 15 lepidopteran species, suggesting that it acts as a genomic hotspot for the diversification of wing patterns, but functional validation through protein-coding knockouts has proven difficult to obtain. Our study unveils the role of a long noncoding RNA (lncRNA) which we name ivory, transcribed from the cortex locus, in modulating color patterning in butterflies. Strikingly, ivory expression prefigures most melanic patterns during pupal development, suggesting an early developmental role in specifying scale identity. To test this, we generated CRISPR mosaic knock-outs in five nymphalid butterfly species and show that ivory mutagenesis yields transformations of dark pigmented scales into white or light-colored scales. Genotyping of Vanessa cardui germline mutants associates these phenotypes to small on-target deletions at the conserved first exon of ivory. In contrast, cortex germline mutant butterflies with confirmed null alleles lack any wing phenotype and exclude a color patterning role for this adjacent gene. Overall, these results show that a lncRNA gene acts as a master switch of color pattern specification and played key roles in the adaptive diversification of wing patterns in butterflies.

The ivory lncRNA regulates seasonal color patterns in buckeye butterflies.

Long noncoding RNAs (lncRNAs) are transcribed elements increasingly recognized for their roles in regulating gene expression. Thus far, however, we have little understanding of how lncRNAs contribute to evolution and adaptation. Here, we show that a conserved lncRNA, ivory, is an important color patterning gene in the buckeye butterfly Junonia coenia. ivory overlaps with cortex, a locus linked to multiple cases of crypsis and mimicry in Lepidoptera. Along with a companion paper by Livraghi et al., we argue that ivory, not cortex, is the color pattern gene of interest at this locus. In J. coenia, a cluster of cis-regulatory elements (CREs) in the first intron of ivory are genetically associated with natural variation in seasonal color pattern plasticity, and targeted deletions of these CREs phenocopy seasonal phenotypes. Deletions of different ivory CREs produce other distinct phenotypes as well, including loss of melanic eyespot rings, and positive and negative changes in overall wing pigmentation. We show that the color pattern transcription factors Spineless, Bric-a-brac, and Ftz-f1 bind to the ivory promoter during wing pattern development, suggesting that they directly regulate ivory. This case study demonstrates how cis-regulation of a single noncoding RNA can exert diverse and nuanced effects on the evolution and development of color patterns, including modulating seasonally plastic color patterns.

Patterns of selection across gene regulatory networks.

Gene regulatory networks (GRNs) are the core engine of organismal development. If we would like to understand the origin and diversification of phenotypes, it is necessary to consider the structure of GRNs in order to reconstruct the links between genetic mutations and phenotypic change. Much of the progress in evolutionary developmental biology, however, has occurred without a nuanced consideration of the evolution of functional relationships between genes, especially in the context of their broader network interactions. Characterizing and comparing GRNs across traits and species in a more detailed way will allow us to determine how network position influences what genes drive adaptive evolution. In this perspective paper, we consider the architecture of developmental GRNs and how positive selection strength may vary across a GRN. We then propose several testable models for these patterns of selection and experimental approaches to test these models.

Mimicking Clinical Rejection Patterns in a Rat Osteomyocutaneous Flap Model of Vascularized Composite Allotransplantation.

INTRODUCTION: Graft loss in vascularized composite allotransplantation (VCA) is more often associated with vasculopathy and chronic rejection (CR) than acute cellular rejection (ACR). We present a rat osteomyocutaneous flap model using titrated tacrolimus administration that mimics the graft rejection patterns in our clinical hand transplant program. Comparison of outcomes in these models support a role for ischemia reperfusion injury (IRI) and microvascular changes in CR of skin and large-vessel vasculopathy. The potential of the surgical models for investigating mechanisms of rejection and vasculopathy in VCA and treatment interventions is presented. MATERIALS AND METHODS: Four rodent groups were evaluated: syngeneic controls (Group 1), allogeneic transient immunosuppression (Group 2), allogeneic suboptimal immunosuppression (Group 3), and allogeneic standard immunosuppression (Group 4). Animals were monitored for ACR, vasculopathy, and CR of the skin. RESULTS: Transient immunosuppression resulted in severe ACR within 2 wk of tacrolimus discontinuation. Standard immunosuppression resulted in minimal rejection but subclinical microvascular changes, including capillary thrombosis and luminal narrowing in arterioles in the donor skin. Further reduction in tacrolimus dose led to femoral vasculopathy and CR of the skin. Surprisingly, femoral vasculopathy was also observed in the syngeneic control group. CONCLUSIONS: Titration of tacrolimus in the allogeneic VCA model resulted in presentations of rejection and vasculopathy similar to those in patients and suggests vasculopathy starts at the microvascular level. This adjustable experimental model will allow the study of variables and interventions, such as external trauma or complement blockade, that may initiate or mitigate vasculopathy and CR in VCA.

Successful lung transplantation using an allograft from a COVID-19-recovered donor: a potential role for subgenomic RNA to guide organ utilization.

Although the risk of SARS-CoV-2 transmission through lung transplantation from acutely infected donors is high, the risks of virus transmission and long-term lung allograft outcomes are not as well described when using pulmonary organs from COVID-19-recovered donors. We describe successful lung transplantation for a COVID-19-related lung injury using lungs from a COVID-19-recovered donor who was retrospectively found to have detectable genomic SARS-CoV-2 RNA in the lung tissue by multiple highly sensitive assays. However, SARS-CoV-2 subgenomic RNA (sgRNA), a marker of viral replication, was not detectable in the donor respiratory tissues. One year after lung transplantation, the recipient has a good functional status, walking 1 mile several times per week without the need for supplemental oxygen and without any evidence of donor-derived SARS-CoV-2 transmission. Our findings highlight the limitations of current clinical laboratory diagnostic assays in detecting the persistence of SARS-CoV-2 RNA in the lung tissue. The persistence of SARS-CoV-2 RNA in the donor tissue did not appear to represent active viral replication via sgRNA testing and, most importantly, did not negatively impact the allograft outcome in the first year after lung transplantation. sgRNA is easily performed and may be a useful assay for assessing viral infectivity in organs from donors with a recent infection.

Chromosome Fusion Affects Genetic Diversity and Evolutionary Turnover of Functional Loci but Consistently Depends on Chromosome Size.

Major changes in chromosome number and structure are linked to a series of evolutionary phenomena, including intrinsic barriers to gene flow or suppression of recombination due to chromosomal rearrangements. However, chromosome rearrangements can also affect the fundamental dynamics of molecular evolution within populations by changing relationships between linked loci and altering rates of recombination. Here, we build chromosome-level assembly Eueides isabella and, together with a recent chromosome-level assembly of Dryas iulia, examine the evolutionary consequences of multiple chromosome fusions in Heliconius butterflies. These assemblies pinpoint fusion points on 10 of the 20 autosomal chromosomes and reveal striking differences in the characteristics of fused and unfused chromosomes. The ten smallest autosomes in D. iulia and E. isabella, which have each fused to a longer chromosome in Heliconius, have higher repeat and GC content, and longer introns than predicted by their chromosome length. When fused, these characteristics change to become more in line with chromosome length. The fusions also led to reduced diversity, which likely reflects increased background selection and selection against introgression between diverging populations, following a reduction in per-base recombination rate. We further show that chromosome size and fusion impact turnover rates of functional loci at a macroevolutionary scale. Together these results provide further evidence that chromosome fusion in Heliconius likely had dramatic effects on population level processes shaping rates of neutral and adaptive divergence. These effects may have impacted patterns of diversification in Heliconius, a classic example of an adaptive radiation.

Metoprolol for the Prevention of Acute Exacerbations of COPD.

BACKGROUND: Observational studies suggest that beta-blockers may reduce the risk of exacerbations and death in patients with moderate or severe chronic obstructive pulmonary disease (COPD), but these findings have not been confirmed in randomized trials. METHODS: In this prospective, randomized trial, we assigned patients between the ages of 40 and 85 years who had COPD to receive either a beta-blocker (extended-release metoprolol) or placebo. All the patients had a clinical history of COPD, along with moderate airflow limitation and an increased risk of exacerbations, as evidenced by a history of exacerbations during the previous year or the prescribed use of supplemental oxygen. We excluded patients who were already taking a beta-blocker or who had an established indication for the use of such drugs. The primary end point was the time until the first exacerbation of COPD during the treatment period, which ranged from 336 to 350 days, depending on the adjusted dose of metoprolol. RESULTS: A total of 532 patients underwent randomization. The mean (±SD) age of the patients was 65.0±7.8 years; the mean forced expiratory volume in 1 second (FEV1) was 41.1±16.3% of the predicted value. The trial was stopped early because of futility with respect to the primary end point and safety concerns. There was no significant between-group difference in the median time until the first exacerbation, which was 202 days in the metoprolol group and 222 days in the placebo group (hazard ratio for metoprolol vs. placebo, 1.05; 95% confidence interval [CI], 0.84 to 1.32; P = 0.66). Metoprolol was associated with a higher risk of exacerbation leading to hospitalization (hazard ratio, 1.91; 95% CI, 1.29 to 2.83). The frequency of side effects that were possibly related to metoprolol was similar in the two groups, as was the overall rate of nonrespiratory serious adverse events. During the treatment period, there were 11 deaths in the metoprolol group and 5 in the placebo group. CONCLUSIONS: Among patients with moderate or severe COPD who did not have an established indication for beta-blocker use, the time until the first COPD exacerbation was similar in the metoprolol group and the placebo group. Hospitalization for exacerbation was more common among the patients treated with metoprolol. (Funded by the Department of Defense; BLOCK COPD ClinicalTrials.gov number, NCT02587351.).

Nusinersen by subcutaneous intrathecal catheter for symptomatic spinal muscular atrophy patients with complex spine anatomy.

INTRODUCTION/AIMS: Intrathecal administration of nusinersen is challenging in patients with spinal muscular atrophy (SMA) who have spine deformities or fusions. We prospectively studied the safety and efficacy of nusinersen administration via an indwelling subcutaneous intrathecal catheter (SIC) for SMA patients with advanced disease. METHODS: Seventeen participants commenced nusinersen therapy between 2.7 and 31.5 years of age and received 9 to 12 doses via SIC. Safety was assessed in all participants. A separate efficacy analysis comprised 11 nonambulatory, treatment-naive SMA patients (18.1 ± 6.8 years) with three SMN2 copies and complex spine anatomy. RESULTS: In the safety analysis, 14 treatment-related adverse events (AEs) occurred among 12 (71%) participants; all were related to the SIC and not nusinersen. Device-related AEs interfered with 2.5% of nusinersen doses. Four SICs (24%) required surgical revision due to mechanical malfunction with or without cerebrospinal fluid leak (n = 2), and one (6%) was removed due to Staphylococcus epidermidis meningitis. In the efficacy analysis, mean performance on the nine-hole peg test improved in dominant (15.9%, P = 0.012) and nondominant (19.0%, P = 0.008) hands and grip strength increased by 44.9% (P = 0.031). We observed no significant changes in motor scales, muscle force, pulmonary function, or SMA biomarkers. All participants in the efficacy cohort reported one or more subjective improvement(s) in endurance, purposeful hand use, arm strength, head control, and/or speech. DISCUSSION: For SMA patients with complex spine anatomy, the SIC allows for reliable outpatient administration of nusinersen that results in meaningful improvements in upper limb function, but introduces risks of technical malfunction and iatrogenic infection.

Classic indicators and diel dissolved oxygen versus trend analysis in assessing eutrophication of potable-water reservoirs.

Potable source-water reservoirs are the main water supplies in many urbanizing regions, yet their long-term responses to cultural eutrophication are poorly documented in comparison with natural lakes, creating major management uncertainties. Here, long-term discrete data (June 2006-June 2018) for classical eutrophication water quality indicators, continuous depth-profile data for dissolved oxygen (DO), and an enhanced hybrid statistical trend analysis model were used to evaluate the eutrophication status of a potable source-water reservoir. Based on classical indicators (nitrogen, N and phosphorus, P concentrations and ratios; phytoplankton biomass as chlorophyll a, chl a; and trophic state indices), the reservoir was eutrophic to hypereutrophic and stoichiometrically imbalanced. Anoxia/hypoxia occurred for 7-8 months annually systemwide, even throughout the water column for days to weeks in some years; and elevated total ammonia (up to ~900 µg tNH3 L-1 ) in surface waters from late summer/fall through late winter/early spring suggested substantial internal legacy nutrient loading. These surprising DO and tNH3 phenomena may characterize many reservoirs in urbanizing areas, and the associated cascade of negative impacts may increasingly affect them under global warming. Total organic carbon (TOC), seasonally influenced by phytoplankton biomass, commonly exceeded 6 mg L-1 , which is problematic for potable-water treatment, and significantly trended up over time. Wet-year inflow dilution influenced an apparent decreasing trend in nutrients within the hypereutrophic upper reservoir, which receives most tributary inputs. Nevertheless, significant reservoirwide trends (increasing total phosphorus [TP], phytoplankton chl a, TOC) and mid- and/or lower region trends (increasing total nitrogen [TN], tNH3 , decreasing TN:TP ratios) suggest that water quality degradation from eutrophication has worsened over time. These findings support broadly applicable recommendations to strengthen protection of potable source-water reservoirs in urbanizing watersheds: (1) protective numeric water quality criteria are needed for TOC as well as TN, TP, and chl a; (2) continuous diel data capture more realistic DO conditions than traditional sampling, and can provide important insights for water treatment managers; and (3) assessment of reservoir eutrophication status to track management progress over time should emphasize classic indicators equally as statistical trends, which are highly sensitive to short-term meteorological forcing.

Parallel evolution of ancient, pleiotropic enhancers underlies butterfly wing pattern mimicry.

Color pattern mimicry in Heliconius butterflies is a classic case study of complex trait adaptation via selection on a few large effect genes. Association studies have linked color pattern variation to a handful of noncoding regions, yet the presumptive cis-regulatory elements (CREs) that control color patterning remain unknown. Here we combine chromatin assays, DNA sequence associations, and genome editing to functionally characterize 5 cis-regulatory elements of the color pattern gene optix We were surprised to find that the cis-regulatory architecture of optix is characterized by pleiotropy and regulatory fragility, where deletion of individual cis-regulatory elements has broad effects on both color pattern and wing vein development. Remarkably, we found orthologous cis-regulatory elements associate with wing pattern convergence of distantly related comimics, suggesting that parallel coevolution of ancestral elements facilitated pattern mimicry. Our results support a model of color pattern evolution in Heliconius where changes to ancient, multifunctional cis-regulatory elements underlie adaptive radiation.

Multiethnic genome-wide and HLA association study of total serum IgE level.

BACKGROUND: Total serum IgE (tIgE) is an important intermediate phenotype of allergic disease. Whole genome genetic association studies across ancestries may identify important determinants of IgE. OBJECTIVE: We aimed to increase understanding of genetic variants affecting tIgE production across the ancestry and allergic disease spectrum by leveraging data from the National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine program; the Consortium on Asthma among African-ancestry Populations in the Americas (CAAPA); and the Atopic Dermatitis Research Network (N = 21,901). METHODS: We performed genome-wide association within strata of study, disease, and ancestry groups, and we combined results via a meta-regression approach that models heterogeneity attributable to ancestry. We also tested for association between HLA alleles called from whole genome sequence data and tIgE, assessing replication of associations in HLA alleles called from genotype array data. RESULTS: We identified 6 loci at genome-wide significance (P < 5 × 10-9), including 4 loci previously reported as genome-wide significant for tIgE, as well as new regions in chr11q13.5 and chr15q22.2, which were also identified in prior genome-wide association studies of atopic dermatitis and asthma. In the HLA allele association study, HLA-A∗02:01 was associated with decreased tIgE level (Pdiscovery = 2 × 10-4; Preplication = 5 × 10-4; Pdiscovery+replication = 4 × 10-7), and HLA-DQB1∗03:02 was strongly associated with decreased tIgE level in Hispanic/Latino ancestry populations (PHispanic/Latino discovery+replication = 8 × 10-8). CONCLUSION: We performed the largest genome-wide association study and HLA association study of tIgE focused on ancestrally diverse populations and found several known tIgE and allergic disease loci that are relevant in non-European ancestry populations.

Epigenetic marker of telomeric age is associated with exacerbations and hospitalizations in chronic obstructive pulmonary disease.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) is an age-related condition that has been associated with early telomere attrition; the clinical implications of telomere shortening in COPD are not well known. In this study we aimed to determine the relationship of the epigenetic regulation of telomeric length in peripheral blood with the risk of exacerbations and hospitalization in patients with COPD. METHODS: Blood DNA methylation profiles were obtained from 292 patients with COPD enrolled in the placebo arm of the Macrolide Azithromycin to Prevent Rapid Worsening of Symptoms Associated with Chronic Obstructive Pulmonary Disease (MACRO) Study and who were followed for 1-year. We calculated telomere length based on DNA methylation markers (DNAmTL) and related this biomarker to the risk of exacerbation and hospitalization and health status (St. George Respiratory Questionnaire [SGRQ]) score over time using a Cox proportional hazards model. We also used linear models to investigate the associations of DNAmTL with the rates of exacerbation and hospitalization (adjusted for chronological age, lung function, race, sex, smoking, body mass index and cell composition). RESULTS: Participants with short DNAmTL demonstrated increased risk of exacerbation (P = 0.02) and hospitalization (P = 0.03) compared to those with longer DNAmTL. DNAmTL age acceleration was associated with higher rates of exacerbation (P = 1.35 × 10-04) and hospitalization (P = 5.21 × 10-03) and poor health status (lower SGRQ scores) independent of chronological age (P = 0.03). CONCLUSION: Telomeric age based on blood DNA methylation is associated with COPD exacerbation and hospitalization and thus a promising biomarker for poor outcomes in COPD.

Locally abundant, endangered Mariana swiftlets impact the abundance, behavior, and body condition of an invasive predator.

Invasive predators are known to have negative consumptive and non-consumptive effects on native species, but few examples show how the abundance of native prey may influence an established invasive predator. We compared invasive brown treesnakes (Boiga irregularis; BTS) found in caves occupied by endangered Mariana swiftlets (Aerodramus bartschi) to snakes found in nearby forests and caves without birds to quantify how the abundance of native avian prey impacts BTS abundance and behavior on Guam. From 2011 to 2017 we removed 151 BTS in caves occupied by swiftlets and never observed BTS in caves without birds. Notable locations included snakes foraging near swiftlets and in holes that allowed cave access and escape from capture. Of 43 BTS with gut contents, 27 (63%) contained swiftlets. BTS in swiftlet-occupied caves had greater fat mass compared to forests, indicating access to swiftlets may increase body condition and promote reproduction. Number of ovarian follicles was significantly greater in female snakes from swiftlet-occupied caves compared to those from ravine, but not limestone forests; evidence of male BTS being more capable of reproduction was limited (i.e., fewer non-discernible but not significantly larger testes in snakes from caves). Assuming other limiting factors are considered, altering the functional response of predators through the modification of caves or interdiction lures to exclude or hinder the largest BTS could bolster swiftlet populations by increasing nesting refugia in currently-occupied caves and facilitate recolonization of historical caves.

Necroptosis and Apoptosis Contribute to Cisplatin and Aminoglycoside Ototoxicity.

Ototoxic side effects of cisplatin and aminoglycosides have been extensively studied, but no therapy is available to date. Sensory hair cells, upon exposure to cisplatin or aminoglycosides, undergo apoptotic and necrotic cell death. Blocking these cell death pathways has therapeutic potential in theory, but incomplete protection and lack of therapeutic targets in the case of necrosis, has hampered the development of clinically applicable drugs. Over the past decade, a novel form of necrosis, termed necroptosis, was established as an alternative cell death pathway. Necroptosis is distinguished from passive necrotic cell death, in that it follows a cellular program, involving the receptor-interacting protein kinase (RIPK) 1 and RIPK3. In this study, we used pharmacological and genetic interventions in the mouse to test the relative contributions of necroptosis and caspase-8-mediated apoptosis toward cisplatin and aminoglycoside ototoxicity. We find that ex vivo, only apoptosis contributes to cisplatin and aminoglycoside ototoxicity, while in vivo, necroptosis as well as apoptosis are involved in both sexes. Inhibition of necroptosis and apoptosis using pharmacological compounds is thus a viable strategy to ameliorate aminoglycoside and cisplatin ototoxicity.SIGNIFICANCE STATEMENT The clinical application of cisplatin and aminoglycosides is limited due to ototoxic side effects. Here, using pharmaceutical and genetic intervention, we present evidence that two types of programmed cell death, apoptosis and necroptosis, contribute to aminoglycoside and cisplatin ototoxicity. Key molecular factors mediating necroptosis are well characterized and druggable, presenting new avenues for pharmaceutical intervention.

Genome-Wide Regulatory Adaptation Shapes Population-Level Genomic Landscapes in Heliconius.

Cis-regulatory evolution is an important engine of organismal diversification. Although recent studies have looked at genomic patterns of regulatory evolution between species, we still have a poor understanding of the magnitude and nature of regulatory variation within species. Here, we examine the evolution of regulatory element activity over wing development in three Heliconius erato butterfly populations to determine how regulatory variation is associated with population structure. We show that intraspecific divergence in chromatin accessibility and regulatory activity is abundant, and that regulatory variants are spatially clustered in the genome. Regions with strong population structure are highly enriched for regulatory variants, and enrichment patterns are associated with developmental stage and gene expression. We also found that variable regulatory elements are particularly enriched in species-specific genomic regions and long interspersed nuclear elements. Our findings suggest that genome-wide selection on chromatin accessibility and regulatory activity is an important force driving patterns of genomic divergence within Heliconius species. This work also provides a resource for the study of gene regulatory evolution in H. erato and other heliconiine butterflies.

Multiple roles forlaccase2 in butterfly wing pigmentation, scale development, and cuticle tanning.

Lepidopteran wing scales play important roles in a number of functions including color patterning and thermoregulation. Despite the importance of wing scales, however, we still have a limited understanding of the genetic mechanisms that underlie scale patterning, development, and coloration. Here, we explore the function of the phenoloxidase-encoding gene laccase2 in wing and scale development in the nymphalid butterfly Vanessa cardui. Somatic deletion mosaics of laccase2 generated by CRISPR/Cas9 genome editing presented several distinct mutant phenotypes. Consistent with the work in other nonlepidopteran insect groups, we observed reductions in melanin pigmentation and defects in cuticle formation. We were also surprised, however, to see distinct effects on scale development including complete loss of wing scales. This study highlights laccase2 as a gene that plays multiple roles in wing and scale development and provides new insight into the evolution of lepidopteran wing coloration.

Long-term exposure to particulate air pollution and brachial artery flow-mediated dilation in the Old Order Amish.

BACKGROUND: Atmospheric particulate matter (PM) has been associated with endothelial dysfunction, an early marker of cardiovascular risk. Our aim was to extend this research to a genetically homogenous, geographically stable rural population using location-specific moving-average air pollution exposure estimates indexed to the date of endothelial function measurement. METHODS: We measured endothelial function using brachial artery flow-mediated dilation (FMD) in 615 community-dwelling healthy Amish participants. Exposures to PM < 2.5 µm (PM2.5) and PM < 10 µm (PM10) were estimated at participants' residential addresses using previously developed geographic information system-based spatio-temporal models and normalized. Associations between PM exposures and FMD were evaluated using linear mixed-effects regression models, and polynomial distributed lag (PDL) models followed by Bayesian model averaging (BMA) were used to assess response to delayed effects occurring across multiple months. RESULTS: Exposure to PM10 was consistently inversely associated with FMD, with the strongest (most negative) association for a 12-month moving average (- 0.09; 95% CI: - 0.15, - 0.03). Associations with PM2.5 were also strongest for a 12-month moving average but were weaker than for PM10 (- 0.07; 95% CI: - 0.13, - 0.09). Associations of PM2.5 and PM10 with FMD were somewhat stronger in men than in women, particularly for PM10. CONCLUSIONS: Using location-specific moving-average air pollution exposure estimates, we have shown that 12-month moving-average estimates of PM2.5 and PM10 exposure are associated with impaired endothelial function in a rural population.

Single master regulatory gene coordinates the evolution and development of butterfly color and iridescence.

The optix gene has been implicated in butterfly wing pattern adaptation by genetic association, mapping, and expression studies. The actual developmental function of this gene has remained unclear, however. Here we used CRISPR/Cas9 genome editing to show that optix plays a fundamental role in nymphalid butterfly wing pattern development, where it is required for determination of all chromatic coloration. optix knockouts in four species show complete replacement of color pigments with melanins, with corresponding changes in pigment-related gene expression, resulting in black and gray butterflies. We also show that optix simultaneously acts as a switch gene for blue structural iridescence in some butterflies, demonstrating simple regulatory coordination of structural and pigmentary coloration. Remarkably, these optix knockouts phenocopy the recurring "black and blue" wing pattern archetype that has arisen on many independent occasions in butterflies. Here we demonstrate a simple genetic basis for structural coloration, and show that optix plays a deeply conserved role in butterfly wing pattern development.

Macroevolutionary shifts of WntA function potentiate butterfly wing-pattern diversity.

Butterfly wing patterns provide a rich comparative framework to study how morphological complexity develops and evolves. Here we used CRISPR/Cas9 somatic mutagenesis to test a patterning role for WntA, a signaling ligand gene previously identified as a hotspot of shape-tuning alleles involved in wing mimicry. We show that WntA loss-of-function causes multiple modifications of pattern elements in seven nymphalid butterfly species. In three butterflies with a conserved wing-pattern arrangement, WntA is necessary for the induction of stripe-like patterns known as symmetry systems and acquired a novel eyespot activator role specific to Vanessa forewings. In two Heliconius species, WntA specifies the boundaries between melanic fields and the light-color patterns that they contour. In the passionvine butterfly Agraulis, WntA removal shows opposite effects on adjacent pattern elements, revealing a dual role across the wing field. Finally, WntA acquired a divergent role in the patterning of interveinous patterns in the monarch, a basal nymphalid butterfly that lacks stripe-like symmetry systems. These results identify WntA as an instructive signal for the prepatterning of a biological system of exuberant diversity and illustrate how shifts in the deployment and effects of a single developmental gene underlie morphological change.