Neurosurgery at Baylor Scott & White Health: Recent History and a Decade of Neurosurgery Residency Training

Neurosurgery at Baylor Scott & White Memorial Hospital in Temple, Texas began as a division in the Department of Surgery many decades ago. The hospital has long served as the flagship tertiary referral center for the Baylor Scott & White healthcare system, which merged in 2013 with Baylor University Medical Center, a hospital system based in Dallas. It is now the largest non-profit hospital system as well as the most awarded hospital system by the US News and World Report within the state of Texas. The Department of Neurosurgery was established at Baylor Scott & White Memorial Hospital in the 2006-2007 academic year. Between then and 2014, four neurosurgeons served as department chair or interim chair: Dr. Robert Buchanan, Dr. Gerhard Friehs, Dr. Ibrahim El Nihum, and Dr. David Garrett Jr. In 2014, Dr. Jason Huang was appointed chairman after a national search and established the neurosurgery residency program in 2015. The department has undergone tremendous growth under the leadership of Dr. Huang, and the residency program is a priority of the department. Surgical excellence is honed at primarily three campuses: Baylor Scott & White Memorial Hospital, Baylor Scott & White McLane Children's Medical Center, and Baylor Scott & White Medical Center - Hillcrest. In this editorial, we provide a brief history of the institution, a recent history of the neurosurgical presence at Baylor Scott & White Memorial Hospital in Temple, Texas, and briefly describe the program's future directions under the continued leadership of Dr. Jason Huang.

Supporting the democratization of science during a pandemic: genomics Course-based Undergraduate Research Experiences (CUREs) as an effective remote learning strategy

The initial phase of the COVID-19 pandemic changed the nature of course delivery from largely in-person to exclusively remote, thus disrupting the well-established pedagogy of the Genomics Education Partnership (GEP; https://www.thegep.org). However, our web-based research adapted well to the remote learning environment. As usual, students who engaged in the GEP's Course-based Undergraduate Research Experience (CURE) received digital projects based on genetic information within assembled Drosophila genomes. Adaptations for remote implementation included moving new member faculty training and peer Teaching Assistant office hours from in-person to online. Surprisingly, our faculty membership significantly increased and, hence, the number of supported students. Furthermore, despite the mostly virtual instruction of the 2020-2021 academic year, there was no significant decline in student learning nor attitudes. Based on successfully expanding the GEP CURE within a virtual learning environment, we provide four strategic lessons we infer toward democratizing science education. First, it appears that increasing access to scientific research and professional development opportunities by supporting virtual, cost-free attendance at national conferences attracts more faculty members to educational initiatives. Second, we observed that transitioning new member training to an online platform removed geographical barriers, reducing time and travel demands, and increased access for diverse faculty to join. Third, developing a Virtual Teaching Assistant program increased the availability of peer support, thereby improving the opportunities for student success. Finally, increasing access to web-based technology is critical for providing equitable opportunities for marginalized students to fully participate in research courses. Online CUREs have great potential for democratizing science education.

Investigating the phylogenetic history of toxin tolerance in mushroom-feeding Drosophila

Understanding how and when key novel adaptations evolved is a central goal of evolutionary biology. Within the immigrans-tripunctata radiation of Drosophila, many mushroom-feeding species are tolerant of host toxins, such as cyclopeptides, that are lethal to nearly all other eukaryotes. In this study, we used phylogenetic and functional approaches to investigate the evolution of cyclopeptide tolerance in the immigrans-tripunctata radiation of Drosophila. First, we inferred the evolutionary relationships among 48 species in this radiation using 978 single copy orthologs. Our results resolved previous incongruities within species groups across the phylogeny. Second, we expanded on previous studies of toxin tolerance by assaying 16 of these species for tolerance to α-amanitin and found that six of them could develop on diet with toxin. Finally, we asked how α-amanitin tolerance might have evolved across the immigrans-tripunctata radiation, and inferred that toxin tolerance was ancestral in mushroom-feeding Drosophila and subsequently lost multiple times. Our findings expand our understanding of toxin tolerance across the immigrans-tripunctata radiation and emphasize the uniqueness of toxin tolerance in this adaptive radiation and the complexity of biochemical adaptations.

Manual annotation of Drosophila genes: a Genomics Education Partnership protocol

Annotating the genomes of multiple species allows us to analyze the evolution of their genes. While many eukaryotic genome assemblies already include computational gene predictions, these predictions can benefit from review and refinement through manual gene annotation. The Genomics Education Partnership (GEP; https://thegep.org/) developed a structural annotation protocol for protein-coding genes that enables undergraduate student and faculty researchers to create high-quality gene annotations that can be utilized in subsequent scientific investigations. For example, this protocol has been utilized by the GEP faculty to engage undergraduate students in the comparative annotation of genes involved in the insulin signaling pathway in 27 Drosophila species, using D. melanogaster as the reference genome. Students construct gene models using multiple lines of computational and empirical evidence including expression data (e.g., RNA-Seq), sequence similarity (e.g., BLAST and multiple sequence alignment), and computational gene predictions. Quality control measures require each gene be annotated by at least two students working independently, followed by reconciliation of the submitted gene models by a more experienced student. This article provides an overview of the annotation protocol and describes how discrepancies in student submitted gene models are resolved to produce a final, high-quality gene set suitable for subsequent analyses. The protocol can be adapted to other scientific questions (e.g., expansion of the Drosophila Muller F element) and species (e.g., parasitoid wasps) to provide additional opportunities for undergraduate students to participate in genomics research. These student annotation efforts can substantially improve the quality of gene annotations in publicly available genomic databases.

Long-read genome assemblies for the study of chromosome expansion: Drosophila kikkawai, Drosophila takahashii, Drosophila bipectinata, and Drosophila ananassae

Flow cytometry estimates of genome sizes among species of Drosophila show a 3-fold variation, ranging from ∼127 Mb in Drosophila mercatorum to ∼400 Mb in Drosophila cyrtoloma. However, the assembled portion of the Muller F element (orthologous to the fourth chromosome in Drosophila melanogaster) shows a nearly 14-fold variation in size, ranging from ∼1.3 Mb to >18 Mb. Here, we present chromosome-level long-read genome assemblies for 4 Drosophila species with expanded F elements ranging in size from 2.3 to 20.5 Mb. Each Muller element is present as a single scaffold in each assembly. These assemblies will enable new insights into the evolutionary causes and consequences of chromosome size expansion.

Investigating the phylogenetic history of toxin tolerance in mushroom-feeding Drosophila

Understanding how and when key novel adaptations evolved is a central goal of evolutionary biology. Within the immigrans-tripunctata radiation of Drosophila , many mushroom-feeding species are tolerant of host toxins, such as cyclopeptides, that are lethal to nearly all other eukaryotes. In this study, we used phylogenetic and functional approaches to investigate the evolution of cyclopeptide tolerance in the immigrans-tripunctata radiation of Drosophila . We first inferred the evolutionary relationships among 48 species in this radiation using 978 single copy orthologs. Our results resolved previous incongruities within species groups across the phylogeny. Second, we expanded on previous studies of toxin tolerance by assaying 16 of these species for tolerance to α-amanitin and found that six of these species could develop on diet with toxin. Third, we examined fly development on a diet containing a natural mix of toxins extracted from the Death Cap Amanita phalloides mushroom. Both tolerant and susceptible species developed on diet with this mix, though tolerant species survived at significantly higher concentrations. Finally, we asked how cyclopeptide tolerance might have evolved across the immigrans-tripunctata radiation and inferred that toxin tolerance was ancestral and subsequently lost multiple times. Our results suggest the evolutionary history of cyclopeptide tolerance is complex, and simply describing this trait as present or absent does not fully capture the occurrence or impact on this adaptive radiation. More broadly, the evolution of novelty can be more complex than previously thought, and that accurate descriptions of such novelties are critical in studies examining their evolution.

Quamdiu? Time to proficiency in endoscope-assisted minimally invasive clot evacuation

INTRODUCTION

Spontaneous intraparenchymal brain hemorrhages are a devastating disease associated with significant disability or death. Minimally invasive clot evacuation (MICE) techniques can reduce mortality. We reviewed our experience with learning endoscope-assisted MICE to determine whether adequate results could be obtained in less than 10 cases.

METHODS

We performed a retrospective chart review of patients undergoing endoscope-assisted MICE at a single institution by a single surgeon from January 1, 2018 to January 1, 2023 using a neuro-endoscope, a commercial clot evacuation device, and frameless stereotaxis. Demographic data was collected along with surgical results and complications. Image analysis using software determined the degree of clot removal. Hospital length of stay and functional outcomes were assessed using the Glasgow Coma Scale score (GCS) and Glasgow Outcome Score (extended) (GOS-E).

RESULTS

Eleven patients were identified: average age 60.82 years old, 64 % male, all had hypertension. There was a clear improvement in IPH evacuation over the series. By case #7, greater than 80 % of clot volume was evacuated consistently. All patients remained neurologically stable or improved following surgery. In long-term follow-up, four patients (36.4 %) had good outcomes (GOS-E ≥ 6) and 2 patients had fair outcomes (GOS-E = 4) (18 %). There were no surgical mortalities, re-hemorrhages, or infections.

CONCLUSIONS

With an experience of less than 10 cases, it is possible to obtain results comparable to most published series of endoscope-assisted MICE. Benchmarks such as greater than 80 % volume removal, less than 15 mL residual, and 40 % good functional outcomes can be obtained.

Long-read genome assemblies for the study of chromosome expansion: Drosophila kikkawai , Drosophila takahashii , Drosophila bipectinata , and Drosophila ananassae

Flow cytometry estimates of genome sizes among species of Drosophila show a 3-fold variation, ranging from ∼127 Mb in Drosophila mercatorum to ∼400 Mb in Drosophila cyrtoloma . However, the assembled portion of the Muller F Element (orthologous to the fourth chromosome in Drosophila melanogaster ) shows a nearly 14-fold variation in size, ranging from ∼1.3 Mb to > 18 Mb. Here, we present chromosome-level long read genome assemblies for four Drosophila species with expanded F Elements ranging in size from 2.3 Mb to 20.5 Mb. Each Muller Element is present as a single scaffold in each assembly. These assemblies will enable new insights into the evolutionary causes and consequences of chromosome size expansion.

Early Surgery for Moyamoya Is Not Associated with Worse Outcomes: Analysis of a North American Adult Cohort

BACKGROUND

Patients with symptomatic moyamoya disease (MMD) or moyamoya syndrome (MMS) are at high risk for recurrent stroke. Surgical revascularization with either direct or indirect superficial temporal artery to the middle cerebral artery (STA-MCA) bypass is a well-accepted treatment. The optimal timing and surgical technique for adult patients with MMD or MMS remains unknown.

METHODS

We performed a retrospective chart review of patients who underwent STA-MCA bypass for MMD or MMS from January 1, 2017, to January 1, 2022. Data collected included demographics, comorbidities, complications, and angiographic/clinical outcomes. Early surgery occurred within two weeks of the last stroke; delayed after two weeks. Statistical analysis compared early vs. delayed surgery and direct vs. indirect bypass.

RESULTS

Nineteen patients underwent bypass surgery on 24 hemispheres. There were 10 early and 14 delayed bypasses. Seventeen were direct, and seven were indirect. There was no statistical difference in total complications in the early (3/10, 30%) vs. the delayed groups (3/14, 21%) (P = 0.67). There were five total complications in the direct group (5/17, 29%) and one in the indirect group (1/7, 14%) (P=0.63). There were no mortalities related to surgery. Angiographic follow-up showed more extensive revascularization with early direct bypass than delayed indirect bypass.

CONCLUSIONS

In this population of North American adults undergoing surgical revascularization for MMD or MMS, early surgery within two weeks of the last stroke did not differ from delayed surgery in terms of complications or outcomes. Early direct bypass showed more revascularization on angiography than delayed indirect surgery.

The Genomics Education Partnership: First findings on genomics research in community colleges

The Genomics Education Partnership (GEP), a consortium of diverse colleges/universities, provides support for integrating genomics research into undergraduate curricula. To increase research opportunities for underrepresented students, GEP is expanding to more community colleges (CC). Genomics research, requiring only a computer with internet access, may be particularly accessible for 2-year institutions with limited research capacity and significant budget constraints. To understand how GEP supports student research at CCs, we analyzed student knowledge and self-reported outcomes. We found that CC student gains are comparable to non-CC student gains, with improvements in attitudes toward science and thriving in science. Our early findings suggest that the GEP model of centralized support with flexible CURE implementation benefits CC students and may help mitigate barriers to implementing research at CCs.

Safety of Contemporary Carotid Artery Stenting and Associated Risk Factors for Treatment of Carotid Artery Stenosis

BACKGROUND

The prevalence of carotid artery atherosclerosis and stroke is increasing. Carotid artery stenting (CAS), which treats carotid artery stenosis and prevents strokes, has benefited from advances in techniques and technology. In this study, we aimed to identify significant risk factors for major complications following contemporary CAS.

METHODS

Utilizing the National Surgical Quality Improvement Program registry, we investigated individuals who underwent CAS between the years 2015-2018. Multivariable logistic regression was performed to identify risk factors associated with stroke, myocardial infarction (MI), and death within 30 days of CAS.

RESULTS

Of the 770 patients within the National Surgical Quality Improvement Program registry who met study criteria, 3.2% experienced strokes, 1.6% experienced MIs, and 2.2% died within 30 days of CAS. A total of 6.1% of patients experienced any of these major complications. Of the 395 asymptomatic patients, 3.5% suffered a stroke, 1.5% experienced an MI, and 1.0% died for a total of 5.6% patients suffering from any event. Of the 375 symptomatic patients, there were 2.9% strokes, 1.6% MIs, and 3.5% deaths with 6.7% suffering any event. Age, race, and hyponatremia were significantly associated with stroke, and hyponatremia was associated with MI. Age, anemia, hypoalbuminemia, international normalized ratio, and emergency case were associated with death.

CONCLUSIONS

These results suggest that preoperative risk stratification of CAS patients may not be straightforward as each adverse event has its own risk factors, precluding a generalized CAS risk stratification scoring system from being created. This paper can aid in patient selection for CAS by highlighting risk factors in modern CAS.

Student Attitudes Contribute to the Effectiveness of a Genomics CURE

The Genomics Education Partnership (GEP) engages students in a course-based undergraduate research experience (CURE). To better understand the student attributes that support success in this CURE, we asked students about their attitudes using previously published scales that measure epistemic beliefs about work and science, interest in science, and grit. We found, in general, that the attitudes students bring with them into the classroom contribute to two outcome measures, namely, learning as assessed by a pre- and postquiz and perceived self-reported benefits. While the GEP CURE produces positive outcomes overall, the students with more positive attitudes toward science, particularly with respect to epistemic beliefs, showed greater gains. The findings indicate the importance of a student's epistemic beliefs to achieving positive learning outcomes.

Drosophila eugracilis - Akt

Gene Model for Akt in the D. eugracilis (DeugGB2) assembly (GCA_000236325.2).

Neurovascular Syndromes

Patients with cerebrovascular syndromes are at risk for additional concerns associated with their syndrome. A wide variety of syndromes are associated with cerebrovascular diseases. Multidisciplinary care is helpful to ensure comprehensive evaluation and management. Precise diagnosis and appreciation for the underlying syndrome is critical for effective cerebrovascular and broader care. This text focuses on these conditions with a focus on underlying pathophysiology and associated genetics, presentation, diagnosis, and management of each disease.

Drosophila yakuba - Tsc1

Gene Model for the ortholog of Tsc1 in the Drosophila yakuba DyakCAF1 assembly (GCA_000005975.1).

PulseRider Treated Aneurysm with Significant Artifact on Postoperative Magnetic Resonance Angiography: A Case Report and Literature Review

The PulseRider is a neuroendovascular adjunct for wide-necked intracranial aneurysms. The decreased metal burden of the PulseRider theoretically reduces artifact on radiologic imaging. However, we report here on a case of a patient who underwent PulseRider-assisted stent-coiling of a basilar tip aneurysm. He returned 19 months later for intermittent diplopia and darkening of vision but was neurologically intact on exam. Both contrast-enhanced and time-of-flight magnetic resonance angiography (MRA) demonstrated absence of signal in the basilar artery in the proximal anchors of the PulseRider. Given his lack of reproducible symptoms and high functional status, it is presumed that the imaging reflected artifact and not thrombosis/stenosis. Although the PulseRider is a useful treatment option for wide-necked intracranial aneurysms, the clinician should be aware that even contrast-enhanced MRA can produce artifact that resembles thrombosis/stenosis. Non-angiogram radiologic imaging modalities may be appropriate for evaluation for residual aneurysm but not patency of the parent artery.

Drosophila yakuba - Tsc1

Gene Model for Tsc1 in the Drosophila yakuba's DyakCAF1 assembly (GCA_000005975.1).

Effect of gut microbiota on α-amanitin tolerance in Drosophila tripunctata

The bacterial gut microbiota of many animals is known to be important for many physiological functions including detoxification. The selective pressures imposed on insects by exposure to toxins may also be selective pressures on their symbiotic bacteria, who thus may contribute to the mechanism of toxin tolerance for the insect. Amatoxins are a class of cyclopeptide mushroom toxins that primarily act by binding to RNA polymerase II and inhibiting transcription. Several species of mycophagous Drosophila are tolerant to amatoxins found in mushrooms of the genus Amanita, despite these toxins being lethal to most other known eukaryotes. These species can tolerate amatoxins in natural concentrations to utilize toxic mushrooms as larval hosts, but the mechanism by which these species are tolerant remains unknown. Previous data have shown that a local population of D. tripunctata exhibits significant genetic variation in toxin tolerance. This study assesses the potential role of the microbiome in α-amanitin tolerance in six wild-derived strains of Drosophila tripunctata. Normal and antibiotic-treated samples of six strains were reared on diets with and without α-amanitin, and then scored for survival from the larval stage to adulthood and for development time to pupation. Our results show that a substantial reduction in bacterial load does not influence toxin tolerance in this system, while confirming genotype and toxin-specific effects on survival are independent of the microbiome composition. Thus, we conclude that this adaptation to exploit toxic mushrooms as a host is likely intrinsic to the fly's genome and not a property of their microbiome.

Exhaustive extraction of cyclopeptides from Amanita phalloides: Guidelines for working with complex mixtures of secondary metabolites

Understanding plant-insect interactions is an active area of research in both ecology and evolution. Much attention has been focused on the impact of secondary metabolites in the host plant or fungi on these interactions. Plants and fungi contain a variety of biologically active compounds, and the secondary metabolite profile can vary significantly between individual samples. However, many experiments characterize the biological effects of only a single secondary metabolite or a subset of these compounds.Here, we develop an exhaustive extraction protocol using an accelerated solvent extraction protocol to recover the complete suite of cyclopeptides and other secondary metabolites found in Amanita phalloides (death cap mushrooms) and compare its efficacy to the "Classic" extraction method used in earlier works.We demonstrate that our extraction protocol recovers the full suite of cyclopeptides and other secondary metabolites in A. phalloides unlike the "Classic" method that favors polar cyclopeptides.Based on these findings, we provide recommendations for how to optimize protocols to ensure exhaustive extracts and also the best practices when using natural extracts in ecological experiments.

Successful endovascular removal of a chondroblastic osteosarcoma embolus within the left middle cerebral artery

A 57-year-old woman with metastatic chondroblastic osteosarcoma developed sudden-onset right-sided weakness and aphasia. She was found to have an occluded left middle cerebral artery (M1 segment) and underwent endovascular stroke intervention with return to a normal neurological exam within 24 h. Histologic analysis of the embolus was consistent with chondroblastic osteosarcoma.

Facilitating Growth through Frustration: Using Genomics Research in a Course-Based Undergraduate Research Experience

A hallmark of the research experience is encountering difficulty and working through those challenges to achieve success. This ability is essential to being a successful scientist, but replicating such challenges in a teaching setting can be difficult. The Genomics Education Partnership (GEP) is a consortium of faculty who engage their students in a genomics Course-Based Undergraduate Research Experience (CURE). Students participate in genome annotation, generating gene models using multiple lines of experimental evidence. Our observations suggested that the students' learning experience is continuous and recursive, frequently beginning with frustration but eventually leading to success as they come up with defendable gene models. In order to explore our "formative frustration" hypothesis, we gathered data from faculty via a survey, and from students via both a general survey and a set of student focus groups. Upon analyzing these data, we found that all three datasets mentioned frustration and struggle, as well as learning and better understanding of the scientific process. Bioinformatics projects are particularly well suited to the process of iteration and refinement because iterations can be performed quickly and are inexpensive in both time and money. Based on these findings, we suggest that a dynamic of "formative frustration" is an important aspect for a successful CURE.

Enhanced responses to somatostatin interneuron activation in developmentally malformed cortex

Intractable epilepsy is commonly associated with developmental cortical malformations. Using the rodent freeze lesion model, we have sought the underlying circuit abnormalities contributing to the epileptiform activity that occurs in association with the structural pathology of four-layered microgyria. We showed previously that within the epileptogenic paramicrogyral region (PMR) surrounding the malformation, non-fast-spiking neurons commonly containing somatostatin (SSt) exhibit alterations, including having a greater maximum firing rate. Here we examined the output of SSt interneurons with optogenetics, using SSt-Cre mice mated to mice with floxed channelrhodopsin-2. Voltage clamp recordings from layer V pyramidal neurons in ex vivo slices had significantly enhanced SSt-evoked inhibitory postsynaptic currents in PMR cortex compared to control. In addition, under conditions of low-Mg2+ artificial cerebral spinal fluid (aCSF), light activation of the SSt neurons evoked field potential epileptiform activity in the PMR cortex, but not in control. These data suggest that within the PMR cortex, SSts have a significantly larger effect on excitatory neurons. Surprisingly, the network effect of this enhanced inhibition is hyperexcitability with propagating epileptiform activity, perhaps due to disinhibition of other interneuron cell types or to enhanced synchrony of excitatory cortical elements. This identification creates a new locus for potential modulation of epileptiform activity associated with cortical malformation.

Random Forest Analysis of Untargeted Metabolomics Data Suggests Increased Use of Omega Fatty Acid Oxidation Pathway in Drosophila Melanogaster Larvae Fed a Medium Chain Fatty Acid Rich High-Fat Diet

Obesity is a complex disease, shaped by both genetic and environmental factors such as diet. In this study, we use untargeted metabolomics and Drosophila melanogaster to model how diet and genotype shape the metabolome of obese phenotypes. We used 16 distinct outbred genotypes of Drosophila larvae raised on normal (ND) and high-fat (HFD) diets, to produce three distinct phenotypic classes; genotypes that stored more triglycerides on a ND relative to the HFD, genotypes that stored more triglycerides on a HFD relative to ND, and genotypes that showed no change in triglyceride storage on either of the two diets. Using untargeted metabolomics we characterized 350 metabolites: 270 with definitive chemical IDs and 80 that were chemically unidentified. Using random forests, we determined metabolites that were important in discriminating between the HFD and ND larvae as well as between the triglyceride phenotypic classes. We found that flies fed on a HFD showed evidence of an increased use of omega fatty acid oxidation pathway, an alternative to the more commonly used beta fatty acid oxidation pathway. Additionally, we observed no correlation between the triglyceride storage phenotype and free fatty acid levels (laurate, caprate, caprylate, caproate), indicating that the distinct metabolic profile of fatty acids in high-fat diet fed Drosophila larvae does not propagate into triglyceride storage differences. However, dipeptides did show moderate differences between the phenotypic classes. We fit Gaussian graphical models (GGMs) of the metabolic profiles for HFD and ND flies to characterize changes in metabolic network structure between the two diets, finding the HFD to have a greater number of edges indicating that metabolome varies more across samples on a HFD. Taken together, these results show that, in the context of obesity, metabolomic profiles under distinct dietary conditions may not be reliable predictors of phenotypic outcomes in a genetically diverse population.

A phylogenetic examination of host use evolution in the quinaria and testacea groups of Drosophila

Adaptive radiations provide an opportunity to examine complex evolutionary processes such as ecological specialization and speciation. While a well-resolved phylogenetic hypothesis is critical to completing such studies, the rapid rates of evolution in these groups can impede phylogenetic studies. Here we study the quinaria and testacea species groups of the immigrans-tripunctata radiation of Drosophila, which represent a recent adaptive radiation and are a developing model system for ecological genetics. We were especially interested in understanding host use evolution in these species. In order to infer a phylogenetic hypothesis for this group we sampled loci from both the nuclear genome and the mitochondrial DNA to develop a dataset of 43 protein-coding loci for these two groups along with their close relatives in the immigrans-tripunctata radiation. We used this dataset to examine their evolutionary relationships along with the evolution of feeding behavior. Our analysis recovers strong support for the monophyly of the testacea but not the quinaria group. Results from our ancestral state reconstruction analysis suggests that the ancestor of the testacea and quinaria groups exhibited mushroom-feeding. Within the quinaria group, we infer that transition to vegetative feeding occurred twice, and that this transition did not coincide with a genome-wide change in the rate of protein evolution.

The TreadWheel: Interval Training Protocol for Gently Induced Exercise in Drosophila melanogaster

The incidence of complex metabolic diseases has increased as a result of a widespread transition towards lifestyles of increased caloric intake and lowered activity levels. These multifactorial diseases arise from a combination of genetic, environmental, and behavioral factors. One such complex disease is Metabolic Syndrome (MetS), which is a cluster of metabolic disorders, including hypertension, hyperglycemia, and abdominal obesity. Exercise and dietary intervention are the primary treatments recommended by doctors to mitigate obesity and its subsequent metabolic diseases. Exercise intervention, in particular aerobic interval training, stimulates favorable changes in the common risk factors for Type 2 Diabetes Mellitus (T2DM), Cardiovascular Disease (CVD), and other conditions. With the influx of evidence describing the therapeutic effect exercise has on metabolic health, establishing a system that models exercise in a controlled setting provides a valuable tool for assessing the effects of exercise in an experimental context. Drosophila melanogaster is a great tool for investigating the physiological and molecular changes that result from exercise intervention. The flies have short lifespans and similar mechanisms of metabolizing nutrients when compared to humans. To induce exercise in Drosophila, we developed a machine called the TreadWheel, which utilizes the fly's innate, negative geotaxis tendency to gently induce climbing. This enables researchers to perform experiments on large cohorts of genetically diverse flies to better understand the genotype-by-environment interactions underlying the effects of exercise on metabolic health.

Automated assembly of species metabolomes through data submission into a public repository

Following similar global efforts to exchange genomic and other biomedical data, global databases in metabolomics have now been established. MetaboLights, the first general purpose, publically available, cross-species, cross-application database in metabolomics, has become the fastest growing data repository at the European Bioinformatics Institute in terms of data volume. Here we present the automated assembly of species metabolomes in MetaboLights, a crucial reference for chemical biology, which is growing through user submissions.

Responsible RAD: Striving for best practices in population genomic studies of adaptation

Two recent articles were written in response to our paper "Breaking RAD: An evaluation of the utility of restriction site associated DNA sequencing scans of adaptation." While we agree with some of the comments made by the authors of these two response papers, we still believe caution should be employed in RADseq studies that aim to detect loci that contribute to adaptation. In this rebuttal, we evaluate the key points made in these papers, attempt to identify a middle ground and make suggestions for responsibly conducting future studies to understand the genomewide mechanisms of adaptation.

Breaking RAD: an evaluation of the utility of restriction site-associated DNA sequencing for genome scans of adaptation

Understanding how and why populations evolve is of fundamental importance to molecular ecology. Restriction site-associated DNA sequencing (RADseq), a popular reduced representation method, has ushered in a new era of genome-scale research for assessing population structure, hybridization, demographic history, phylogeography and migration. RADseq has also been widely used to conduct genome scans to detect loci involved in adaptive divergence among natural populations. Here, we examine the capacity of those RADseq-based genome scan studies to detect loci involved in local adaptation. To understand what proportion of the genome is missed by RADseq studies, we developed a simple model using different numbers of RAD-tags, genome sizes and extents of linkage disequilibrium (length of haplotype blocks). Under the best-case modelling scenario, we found that RADseq using six- or eight-base pair cutting restriction enzymes would fail to sample many regions of the genome, especially for species with short linkage disequilibrium. We then surveyed recent studies that have used RADseq for genome scans and found that the median density of markers across these studies was 4.08 RAD-tag markers per megabase (one marker per 245 kb). The length of linkage disequilibrium for many species is one to three orders of magnitude less than density of the typical recent RADseq study. Thus, we conclude that genome scans based on RADseq data alone, while useful for studies of neutral genetic variation and genetic population structure, will likely miss many loci under selection in studies of local adaptation.

Eigenvector metabolite analysis reveals dietary effects on the association among metabolite correlation patterns, gene expression, and phenotypes

INTRODUCTION

'Multi-omics' datasets obtained from an organism of interest reared under different environmental treatments are increasingly common. Identifying the links among metabolites and transcripts can help to elucidate our understanding of the impact of environment at different levels within the organism. However, many methods for characterizing physiological connections cannot address unidentified metabolites.

OBJECTIVES

Here, we use Eigenvector Metabolite Analysis (EvMA) to examine links between metabolomic, transcriptomic, and phenotypic variation data and to assess the impact of environmental factors on these associations. Unlike other methods, EvMA can be used to analyze datasets that include unidentified metabolites and unannotated transcripts.

METHODS

To demonstrate the utility of EvMA, we analyzed metabolomic, transcriptomic, and phenotypic datasets produced from 20 Drosophila melanogaster genotypes reared on four dietary treatments. We used a hierarchical distance-based method to cluster the metabolites. The links between metabolite clusters, gene expression, and overt phenotypes were characterized using the eigenmetabolite (first principal component) of each cluster.

RESULTS

EvMA recovered chemically related groups of metabolites within the clusters. Using the eigenmetabolite, we identified genes and phenotypes that significantly correlated with each cluster. EvMA identifies new connections between the phenotypes, metabolites, and gene transcripts. Conclusion EvMA provides a simple method to identify correlations between metabolites, gene expression, and phenotypes, which can allow us to partition multivariate datasets into meaningful biological modules and identify under-studied metabolites and unannotated gene transcripts that may be central to important biological processes. This can be used to inform our understanding of the effect of environmental mechanisms underlying physiological states of interest.

Finding the Genomic Basis of Local Adaptation: Pitfalls, Practical Solutions, and Future Directions

Uncovering the genetic and evolutionary basis of local adaptation is a major focus of evolutionary biology. The recent development of cost-effective methods for obtaining high-quality genome-scale data makes it possible to identify some of the loci responsible for adaptive differences among populations. Two basic approaches for identifying putatively locally adaptive loci have been developed and are broadly used: one that identifies loci with unusually high genetic differentiation among populations (differentiation outlier methods) and one that searches for correlations between local population allele frequencies and local environments (genetic-environment association methods). Here, we review the promises and challenges of these genome scan methods, including correcting for the confounding influence of a species' demographic history, biases caused by missing aspects of the genome, matching scales of environmental data with population structure, and other statistical considerations. In each case, we make suggestions for best practices for maximizing the accuracy and efficiency of genome scans to detect the underlying genetic basis of local adaptation. With attention to their current limitations, genome scan methods can be an important tool in finding the genetic basis of adaptive evolutionary change.

The changing biodiversity of Alabama Drosophila: important impacts of seasonal variation, urbanization, and invasive species

Global warming and anthropogenic disturbances significantly influence the biosphere, tremendously increasing species extinction rates. In Central Alabama, we analyzed Drosophilidae species composition change nearly 100 years after the previous survey. We found ten Drosophilid species that were not reported during the last major biodiversity studies, two of which are invasive pests. In addition, we analyzed the influence of seasonal environmental variables characteristic of the subtropical climate zone on Drosophila abundance and biodiversity. We found a significant correlation between temperature and abundance of total Drosophila as well as for six of the seven most represented species individually, with a maximum abundance at intermediate temperatures (18–26°C). In addition, temperature was positively correlated with biodiversity of Drosophila. Precipitation produced a significant effect on the abundance of five species of Drosophila, with different optima for each species, but did not affect overall biodiversity. We concluded that in the subtropical climate zone of Central Alabama, seasonal temperature and precipitation changes produce a significant effect on Drosophila abundance and biodiversity, while local land use also impacts fly abundance, contributing to an apparent shift in species composition over the last century. We expect global climate change and other anthropogenic factors to further impact Drosophila species composition in the subtropical climate zone into the future.

Genetic and Sex-Specific Transgenerational Effects of a High Fat Diet in Drosophila melanogaster

An organism's phenotype is the product of its environment and genotype, but an ancestor's environment can also be a contributing factor. The recent increase in caloric intake and decrease in physical activity of developed nations' populations is contributing to deteriorating health and making the study of the longer term impacts of a changing lifestyle a priority. The dietary habits of ancestors have been shown to affect phenotype in several organisms, including humans, mice, and the fruit fly. Whether the ancestral dietary effect is purely environmental or if there is a genetic interaction with the environment passed down for multiple generations, has not been determined previously. Here we used the fruit fly, Drosophila melanogaster, to investigate the genetic, sex-specific, and environmental effects of a high fat diet for three generations' on pupal body weights across ten genotypes. We also tested for genotype-specific transgenerational effects on metabolic pools and egg size across three genotypes. We showed that there were substantial differences in transgenerational responses to ancestral diet between genotypes and sexes through both first and second descendant generations. Additionally, there were differences in phenotypes between maternally and paternally inherited dietary effects. We also found a treated organism's reaction to a high fat diet was not a consistent predictor of its untreated descendants' phenotype. The implication of these results is that, given our interest in understanding and preventing metabolic diseases like obesity, we need to consider the contribution of ancestral environmental experiences. However, we need to be cautious when drawing population-level generalization from small studies because transgenerational effects are likely to exhibit substantial sex and genotype specificity.

Evolutionary Divergence of Circadian and Photoperiodic Phenotypes in the Pitcher-Plant Mosquito, Wyeomyia smithii

For decades, chronobiologists have investigated the relationship between the circadian clock that mediates daily activities and the photoperiodic timer that mediates seasonal activities. The main experiment used to infer a circadian basis for photoperiodic time measurement is the Nanda-Hamner protocol (NH). Herein, the authors compare additive and nonadditive (dominance and epistasis) genetic effects that lead to the divergence of populations of the pitcher-plant mosquito, Wyeomyia smithii, for critical photoperiod (CPP) and amplitude of the rhythmic response to NH for 3 temporal-geographic scales: 1) Over geological time between populations in northern and southern clades, 2) over millennial time between populations within the northern clade, and 3) over generational time between lines selected for long and short CPP from within a single population. The authors show that the pattern of additive, dominance, and epistatic effects depends on the time scale over which populations or lines have diverged. Patterns for genetic differences between populations for CPP and response to NH reveal similarities over geological and millennial time scales but differences over shorter periods of evolution. These results, and the observation that neither the period nor amplitude of the NH rhythm are significantly correlated with CPP among populations, lead the authors to conclude that the rhythmic response to NH has evolved independently of photoperiodic response in populations of W. smithii. The implication is that in this species, genetic modification of the circadian clock has not been the basis for the adaptive modification of photoperiodic time measurement over the climatic gradient of North America.

The Time Is Right to Focus on Model Organism Metabolomes

Model organisms are an essential component of biological and biomedical research that can be used to study specific biological processes. These organisms are in part selected for facile experimental study. However, just as importantly, intensive study of a small number of model organisms yields important synergies as discoveries in one area of science for a given organism shed light on biological processes in other areas, even for other organisms. Furthermore, the extensive knowledge bases compiled for each model organism enable systems-level understandings of these species, which enhance the overall biological and biomedical knowledge for all organisms, including humans. Building upon extensive genomics research, we argue that the time is now right to focus intensively on model organism metabolomes. We propose a grand challenge for metabolomics studies of model organisms: to identify and map all metabolites onto metabolic pathways, to develop quantitative metabolic models for model organisms, and to relate organism metabolic pathways within the context of evolutionary metabolomics, i.e., phylometabolomics. These efforts should focus on a series of established model organisms in microbial, animal and plant research.

A course-based research experience: how benefits change with increased investment in instructional time

There is widespread agreement that science, technology, engineering, and mathematics programs should provide undergraduates with research experience. Practical issues and limited resources, however, make this a challenge. We have developed a bioinformatics project that provides a course-based research experience for students at a diverse group of schools and offers the opportunity to tailor this experience to local curriculum and institution-specific student needs. We assessed both attitude and knowledge gains, looking for insights into how students respond given this wide range of curricular and institutional variables. While different approaches all appear to result in learning gains, we find that a significant investment of course time is required to enable students to show gains commensurate to a summer research experience. An alumni survey revealed that time spent on a research project is also a significant factor in the value former students assign to the experience one or more years later. We conclude: 1) implementation of a bioinformatics project within the biology curriculum provides a mechanism for successfully engaging large numbers of students in undergraduate research; 2) benefits to students are achievable at a wide variety of academic institutions; and 3) successful implementation of course-based research experiences requires significant investment of instructional time for students to gain full benefit.

Genetic architecture of hybrid male sterility in Drosophila: analysis of intraspecies variation for interspecies isolation

BACKGROUND

The genetic basis of postzygotic isolation is a central puzzle in evolutionary biology. Evolutionary forces causing hybrid sterility or inviability act on the responsible genes while they still are polymorphic, thus we have to study these traits as they arise, before isolation is complete.

METHODOLOGY/PRINCIPAL FINDINGS

Isofemale strains of D. mojavensis vary significantly in their production of sterile F(1) sons when females are crossed to D. arizonae males. We took advantage of the intraspecific polymorphism, in a novel design, to perform quantitative trait locus (QTL) mapping analyses directly on F(1) hybrid male sterility itself. We found that the genetic architecture of the polymorphism for hybrid male sterility (HMS) in the F(1) is complex, involving multiple QTL, epistasis, and cytoplasmic effects.

CONCLUSIONS/SIGNIFICANCE

The role of extensive intraspecific polymorphism, multiple QTL, and epistatic interactions in HMS in this young species pair shows that HMS is arising as a complex trait in this system. Directional selection alone would be unlikely to maintain polymorphism at multiple loci, thus we hypothesize that directional selection is unlikely to be the only evolutionary force influencing postzygotic isolation.

Polytene chromosomal maps of 11 Drosophila species: the order of genomic scaffolds inferred from genetic and physical maps

The sequencing of the 12 genomes of members of the genus Drosophila was taken as an opportunity to reevaluate the genetic and physical maps for 11 of the species, in part to aid in the mapping of assembled scaffolds. Here, we present an overview of the importance of cytogenetic maps to Drosophila biology and to the concepts of chromosomal evolution. Physical and genetic markers were used to anchor the genome assembly scaffolds to the polytene chromosomal maps for each species. In addition, a computational approach was used to anchor smaller scaffolds on the basis of the analysis of syntenic blocks. We present the chromosomal map data from each of the 11 sequenced non-Drosophila melanogaster species as a series of sections. Each section reviews the history of the polytene chromosome maps for each species, presents the new polytene chromosome maps, and anchors the genomic scaffolds to the cytological maps using genetic and physical markers. The mapping data agree with Muller's idea that the majority of Drosophila genes are syntenic. Despite the conservation of genes within homologous chromosome arms across species, the karyotypes of these species have changed through the fusion of chromosomal arms followed by subsequent rearrangement events.

Multilocus nuclear sequences reveal intra- and interspecific relationships among chromosomally polymorphic species of cactophilic Drosophila

Drosophila mojavensis and Drosophila arizonae, a pair of sibling species endemic to North America, constitute an important model system to study ecological genetics and the evolution of reproductive isolation. This species pair can produce fertile hybrids in some crosses and are sympatric in a large part of their ranges. Despite the potential for hybridization in nature, however, evidence of introgression has not been rigorously sought. Further, the evolutionary relationships within and among the geographically distant populations of the two species have not been characterized in detail using high-resolution molecular studies. Both species have six chromosomes: five large acrocentrics and one 'dot' chromosome. Fixed inversion differences between the species exist in three chromosomes (X, 2 and 3) while three are colinear (4, 5 and 6), suggesting that were introgression to occur, it would be most likely in the colinear chromosomes. We utilized nucleotide sequence variation at multiple loci on five chromosomes to test for evidence of introgression, and to test various scenarios for the evolutionary relationships of these two species and their populations. While we do not find evidence of recent introgression, loci in the colinear chromosomes appear to have participated in exchange in the past. We also found considerable population structure within both species. The level of differentiation discovered among D. arizonae populations was unexpectedly high and suggests that its populations, as well as those of D. mojavensis, may be themselves undergoing incipient speciation and merit further attention.

Evolutionary relationships of Drosophila mojavensis geographic host races and their sister species Drosophila arizonae

The cactophilic Drosophila mojavensis species group living in the deserts and dry tropical forests of the southwestern United States and Mexico provides a valuable system for studies in diversification and speciation. Rigorous studies of the relationships between host races of D. mojavensis and the relationships among the members of the species group (D. mojavensis, Drosophila arizona, and Drosophila navojoa) are lacking. We used mitochondrial CO1 sequence data to address the phylogenetics and population genetics of this species group. In this study we have found that the sister species D. mojavensis and D. arizonae share no mitochondrial haplotypes and thus show no evidence for recent introgression. We estimate the divergence time between D. mojavensis and D. arizonae to be between 1.91 and 2.97 million years ago. D. arizonae shows little structure in our population genetic analyses but there is phylogenetic differentiation between southeastern and northern populations of D. arizonae. Drosophila mojavensis shows significant population and phylogenetic structure across the four geographic regions of its distribution. The mitochondrial data support an origin of D. mojavensis on the mainland with early differentiation into the populations now found in the Mojave Desert and the Mainland Sonoran Desert and later colonization of the Baja Peninsula, in contrast to previous models. Also, the sister clade to D. mojavensis/D. arizonae includes D. navojoa and Drosophila huaylasi. By defining the genetic relationships among these populations, we provide a foundation for more sophisticated hypothesis testing regarding the timing of early speciation events and host switches in this species group.