Identifying causal variants by fine mapping across multiple studies

Increasingly large Genome-Wide Association Studies (GWAS) have yielded numerous variants associated with many complex traits, motivating the development of “fine mapping” methods to identify which of the associated variants are causal. Additionally, GWAS of the same trait for different populations are increasingly available, raising the possibility of refining fine mapping results further by leveraging different linkage disequilibrium (LD) structures across studies. Here, we introduce multiple study causal variants identification in associated regions (MsCAVIAR), a method that extends the popular CAVIAR fine mapping framework to a multiple study setting using a random effects model. MsCAVIAR only requires summary statistics and LD as input, accounts for uncertainty in association statistics using a multivariate normal model, allows for multiple causal variants at a locus, and explicitly models the possibility of different SNP effect sizes in different populations. We demonstrate the efficacy of MsCAVIAR in both a simulation study and a trans-ethnic, trans-biobank fine mapping analysis of High Density Lipoprotein (HDL).     

Local Recurrence of Soft Tissue Sarcoma Revisited: Is there a Role for “Selective” Radiation?

BackgroundRadiation therapy (RT) is often utilized in cases of high-grade soft tissue sarcoma (STS), but there remain situations where treatment is with surgical excision alone. Our goals were to determine (1) the local recurrence (LR) rate with and without perioperative RT and (2) associations between local recurrence, patient, tumor, and treatment variables.MethodsWe performed a retrospective review of 165 consecutive STS patients. A Cox proportional hazards model was used to investigate variables associated with local recurrence.ResultsLR occurred in 15/78 (19%) without RT, 4/29 (14%) with postoperative RT, and 0/58 with preoperative RT (p=0.002). We found increased rates of local recurrence at 24 months for myxofibrosarcoma (p=0.001) and no-RT (p=0.003). Myxofibrosarcoma accounted for 33 (20%) of the study patients and 12 (63%) of the local recurrences.ConclusionThe LR rate in patients treated with surgery alone was disproportionately attributable to myxofibrosarcoma (11/23 cases, 48%). Other subtypes demonstrated a lower rate of LR in the absence of RT (4/55 cases, 7%), and no LR occurred when final margins were >2 mm. In certain circumstances treatment with a negative margin surgical resection followed by close observation is justifiable. RT is effective and should continue to be considered routinely in myxofibrosarcoma or when surgical margins are inadequate. Level of Evidence: III     

Intrathecal Morphine Use in Adolescent Idiopathic Scoliosis Surgery is Associated with Decreased Opioid Use and Decreased Length of Stay

BackgroundLength of stay (LOS) in the hospital following posterior spinal fusion (PSF) for adolescent idiopathic scoliosis (AIS) has decreased over the past decade due to well-defined postoperative clinical pathways, earlier mobilization, and improved pain control methods. Historically, liberal use of parenteral and oral opioids for pain control caused side effects, resulting in delayed discharge. Intraoperative intrathecal morphine (ITM) has been posited to reduce the need for postoperative opioids and to expedite the discharge process. This study examines the relationship between the use of ITM with average required postoperative opioid usage and with average LOS.MethodsThis IRB-approved retrospective cohort study examined 105 patients with AIS who received PSF with instrumentation split into two cohorts. One cohort underwent PSF via standard surgical protocol (n=40) while the other cohort received intraoperative ITM with the standard surgical protocol (n=65). Power analysis demonstrated a study power of 0.8. LOS and total postoperative opioid analgesic medication (morphine milligram equivalent, MME) data were collected. Age at surgery, gender, number of spinal levels fused, estimated intraoperative blood loss (EBL), preoperative Cobb angle, and any complications related to the use of ITM were also recorded. Continuous variables were analyzed with Student’s t-test and categorical variables were analyzed with chi-square independent-sample tests using SAS 9.4 (α = 0.05).ResultsPatients who were treated with ITM displayed shorter LOS (p<0.0001) and reduced postoperative analgesic requirement (p<0.0001). Patients who received ITM spent an average of 1.8 fewer midnights in the hospital and received an average of 221.2 MME less than patients who received standard protocol (57% decrease). There were no significant differences between the two groups for any other variable.ConclusionIntraoperative ITM is a simple and effective treatment for scoliosis surgeons to better control postoperative pain in patients, reduce the risk of dependency, and achieve earlier discharge from the hospital. Shortened LOS reduces the overall cost of care, benefitting patients, hospitals, and insurance companies. Based on the results of this study and several earlier studies, the authors recommended that scoliosis surgeons consider incorporating use of ITM into their standard operative protocols. Level of Evidence: IV     

The Giardia ventrolateral flange is a lamellar membrane protrusion that supports attachment

Attachment to the intestinal epithelium is critical to the lifestyle of the ubiquitous parasite Giardia lamblia. The ventrolateral flange is a sheet-like membrane protrusion at the interface between parasites and attached surfaces. This structure has been implicated in attachment, but its role has been poorly defined. Here, we identified a novel actin associated protein with putative WH2-like actin binding domains we named Flangin. Flangin complexes with Giardia actin (GlActin) and is enriched in the ventrolateral flange making it a valuable marker for studying the flanges’ role in Giardia biology. Live imaging revealed that the flange grows to around 1 μm in width after cytokinesis, then remains uniform in size during interphase, grows in mitosis, and is resorbed during cytokinesis. A flangin truncation mutant stabilizes the flange and blocks cytokinesis, indicating that flange disassembly is necessary for rapid myosin-independent cytokinesis in Giardia. Rho family GTPases are important regulators of membrane protrusions and GlRac, the sole Rho family GTPase in Giardia, was localized to the flange. Knockdown of Flangin, GlActin, and GlRac result in flange formation defects. This indicates a conserved role for GlRac and GlActin in forming membrane protrusions, despite the absence of canonical actin binding proteins that link Rho GTPase signaling to lamellipodia formation. Flangin-depleted parasites had reduced surface contact and when challenged with fluid shear force in flow chambers they had a reduced ability to remain attached, confirming a role for the flange in attachment. This secondary attachment mechanism complements the microtubule based adhesive ventral disc, a feature that may be particularly important during mitosis when the parental ventral disc disassembles in preparation for cytokinesis. This work supports the emerging view that Giardia’s unconventional actin cytoskeleton has an important role in supporting parasite attachment.     

Diverse lineages of pathogenic Leptospira species are widespread in the environment in Puerto Rico,USA

BackgroundLeptospirosis, caused by Leptospira bacteria, is a common zoonosis worldwide, especially in the tropics. Reservoir species and risk factors have been identified but surveys for environmental sources are rare. Furthermore, understanding of environmental Leptospira containing virulence associated genes and possibly capable of causing disease is incomplete, which may convolute leptospirosis diagnosis, prevention, and epidemiology.Methodology/Principal findingsWe collected environmental samples from 22 sites in Puerto Rico during three sampling periods over 14-months (Dec 2018-Feb 2020); 10 water and 10 soil samples were collected at each site. Samples were screened for DNA from potentially pathogenic Leptospira using the lipL32 PCR assay and positive samples were sequenced to assess genetic diversity. One urban site in San Juan was sampled three times over 14 months to assess persistence in soil; live leptospires were obtained during the last sampling period. Isolates were whole genome sequenced and LipL32 expression was assessed in vitro.We detected pathogenic Leptospira DNA at 15/22 sites; both soil and water were positive at 5/15 sites. We recovered lipL32 sequences from 83/86 positive samples (15/15 positive sites) and secY sequences from 32/86 (10/15 sites); multiple genotypes were identified at 12 sites. These sequences revealed significant diversity across samples, including four novel lipL32 phylogenetic clades within the pathogenic P1 group. Most samples from the serially sampled site were lipL32 positive at each time point. We sequenced the genomes of six saprophytic and two pathogenic Leptospira isolates; the latter represent a novel pathogenic Leptospira species likely belonging to a new serogroup.Conclusions/SignificanceDiverse and novel pathogenic Leptospira are widespread in the environment in Puerto Rico. The disease potential of these lineages is unknown but several were consistently detected for >1 year in soil, which could contaminate water. This work increases understanding of environmental Leptospira diversity and should improve leptospirosis surveillance and diagnostics.     

Moths complement bumblebee pollination of red clover: a case for day-and-night insect surveillance

Recent decades have seen a surge in awareness about insect pollinator declines. Social bees receive the most attention, but most flower-visiting species are lesser known, non-bee insects. Nocturnal flower visitors, e.g. moths, are especially difficult to observe and largely ignored in pollination studies. Clearly, achieving balanced monitoring of all pollinator taxa represents a major scientific challenge. Here, we use time-lapse cameras for season-wide, day-and-night pollinator surveillance of Trifolium pratense (L.; red clover) in an alpine grassland. We reveal the first evidence to suggest that moths, mainly Noctua pronuba (L.; large yellow underwing), pollinate this important wildflower and forage crop, providing 34% of visits (bumblebees: 61%). This is a remarkable finding; moths have received no recognition throughout a century of T. pratense pollinator research. We conclude that despite a non-negligible frequency and duration of nocturnal flower visits, nocturnal pollinators of T. pratense have been systematically overlooked. We further show how the relationship between visitation and seed set may only become clear after accounting for moth visits. As such, population trends in moths, as well as bees, could profoundly affect T. pratense seed yield. Ultimately, camera surveillance gives fair representation to non-bee pollinators and lays a foundation for automated monitoring of species interactions in future.     

Host plant preference based on salicylate chemistry in a willow leaf beetle (Chrysomela aeneicollis)

Summary Chrysomela aeneicollis (Coleoptera: Chrysomelidae) uses salicin from its host plant (Salix spp.) to produce a defensive secretion, salicylaldehyde. Because it requires salicin for this secretion, I predicted that C. aeneicollis should be attracted to willows which possess salicin and other salicylates. To test this prediction, I determined the host-plant preferences of C. aeneicollis among four potential hosts which occur in the Sierra Nevada range of eastern California. These species have very different salicylate chemistries but do not differ in nutritional quality for C. aeneicollis. In oviposition-preference tests, gravid females showed no preference between a salicylate-poor species, S. lutea, and a salicylate-rich species, S. orestera. However in feeding-choice tests, both larvae and adults preferred S. orestera over S. lutea. This preference was not affected by the species on which the larvae were reared. In other feeding tests, adults preferred S. orestera over two medium-salicylate species, S. boothi and S. geyeriana, regardless of which host species they had been feeding on in nature. In a final feeding test, adults were stimulated to feed by salicin itself. In nature, the relative abundances of C. aeneicollis adults and egg clutches among these species correspond to the adult feeding preference in the laboratory. Additionally, multiple regression analyses showed that adult abundance was not related to among-clone differences in leaf toughness or nutritional quality, but rather to salicin content and plant size. Thus for C. aeneicollis, both laboratory and field results demonstrate a preference for salicylate-rich willows which is partly responsible for the increased level of attack on them.     

Complementary genetic and genomic approaches help characterize the linkage group I seed protein QTL in soybean

Background

Cotton fibers (produced by Gossypium species) are the premier natural fibers for textile production. The two tetraploid species, G. barbadense (Gb) and G. hirsutum (Gh), differ significantly in their fiber properties, the former having much longer, finer and stronger fibers that are highly prized. A better understanding of the genetics and underlying biological causes of these differences will aid further improvement of cotton quality through breeding and biotechnology. We evaluated an inter-specific Gh × Gb recombinant inbred line (RIL) population for fiber characteristics in 11 independent experiments under field and glasshouse conditions. Sites were located on 4 continents and 5 countries and some locations were analyzed over multiple years.

Results

The RIL population displayed a large variability for all major fiber traits. QTL analyses were performed on a per-site basis by composite interval mapping. Among the 651 putative QTLs (LOD > 2), 167 had a LOD exceeding permutation based thresholds. Coincidence in QTL location across data sets was assessed for the fiber trait categories strength, elongation, length, length uniformity, fineness/maturity, and color. A meta-analysis of more than a thousand putative QTLs was conducted with MetaQTL software to integrate QTL data from the RIL and 3 backcross populations (from the same parents) and to compare them with the literature. Although the global level of congruence across experiments and populations was generally moderate, the QTL clustering was possible for 30 trait x chromosome combinations (5 traits in 19 different chromosomes) where an effective co-localization of unidirectional (similar sign of additivity) QTLs from at least 5 different data sets was observed. Most consistent meta-clusters were identified for fiber color on chromosomes c6, c8 and c25, fineness on c15, and fiber length on c3.

Conclusions

Meta-analysis provided a reliable means of integrating phenotypic and genetic mapping data across multiple populations and environments for complex fiber traits. The consistent chromosomal regions contributing to fiber quality traits constitute good candidates for the further dissection of the genetic and genomic factors underlying important fiber characteristics, and for marker-assisted selection.     

Evidence for widespread genomic methylation in the migratory locust, Locusta migratoria (Orthoptera: Acrididae)

The importance of DNA methylation in mammalian and plant systems is well established. In recent years there has been renewed interest in DNA methylation in insects. Accumulating evidence, both from mammals and insects, points towards an emerging role for DNA methylation in the regulation of phenotypic plasticity. The migratory locust (Locusta migratoria) is a model organism for the study of phenotypic plasticity. Despite this, there is little information available about the degree to which the genome is methylated in this species and genes encoding methylation machinery have not been previously identified. We therefore undertook an initial investigation to establish the presence of a functional DNA methylation system in L. migratoria. We found that the migratory locust possesses genes that putatively encode methylation machinery (DNA methyltransferases and a methyl-binding domain protein) and exhibits genomic methylation, some of which appears to be localised to repetitive regions of the genome. We have also identified a distinct group of genes within the L. migratoria genome that appear to have been historically methylated and show some possible functional differentiation. These results will facilitate more detailed research into the functional significance of DNA methylation in locusts.     

Exome sequencing and arrayCGH detection of gene sequence and copy number variation between ILS and ISS mouse strains

It has been well documented that genetic factors can influence predisposition to develop alcoholism. While the underlying genomic changes may be of several types, two of the most common and disease associated are copy number variations (CNVs) and sequence alterations of protein coding regions. The goal of this study was to identify CNVs and single-nucleotide polymorphisms that occur in gene coding regions that may play a role in influencing the risk of an individual developing alcoholism. Toward this end, two mouse strains were used that have been selectively bred based on their differential sensitivity to alcohol: the Inbred long sleep (ILS) and Inbred short sleep (ISS) mouse strains. Differences in initial response to alcohol have been linked to risk for alcoholism, and the ILS/ISS strains are used to investigate the genetics of initial sensitivity to alcohol. Array comparative genomic hybridization (arrayCGH) and exome sequencing were conducted to identify CNVs and gene coding sequence differences, respectively, between ILS and ISS mice. Mouse arrayCGH was performed using catalog Agilent 1 × 244 k mouse arrays. Subsequently, exome sequencing was carried out using an Illumina HiSeq 2000 instrument. ArrayCGH detected 74 CNVs that were strain-specific (38 ILS/36 ISS), including several ISS-specific deletions that contained genes implicated in brain function and neurotransmitter release. Among several interesting coding variations detected by exome sequencing was the gain of a premature stop codon in the alpha-amylase 2B (AMY2B) gene specifically in the ILS strain. In total, exome sequencing detected 2,597 and 1,768 strain-specific exonic gene variants in the ILS and ISS mice, respectively. This study represents the most comprehensive and detailed genomic comparison of ILS and ISS mouse strains to date. The two complementary genome-wide approaches identified strain-specific CNVs and gene coding sequence variations that should provide strong candidates to contribute to the alcohol-related phenotypic differences associated with these strains.