Current status of community resources and priorities for weed genomics research.

Weeds are attractive models for basic and applied research due to their impacts on agricultural systems and capacity to swiftly adapt in response to anthropogenic selection pressures. Currently, a lack of genomic information precludes research to elucidate the genetic basis of rapid adaptation for important traits like herbicide resistance and stress tolerance and the effect of evolutionary mechanisms on wild populations. The International Weed Genomics Consortium is a collaborative group of scientists focused on developing genomic resources to impact research into sustainable, effective weed control methods and to provide insights about stress tolerance and adaptation to assist crop breeding.

Polygenic contributions to performance on the Balloon Analogue Risk Task.

Risky decision-making is a common, heritable endophenotype seen across many psychiatric disorders. Its underlying genetic architecture is incompletely explored. We examined behavior in the Balloon Analogue Risk Task (BART), which tests risky decision-making, in two independent samples of European ancestry. One sample (n = 1138) comprised healthy participants and some psychiatric patients (53 schizophrenia, 42 bipolar disorder, 47 ADHD); the other (n = 911) excluded for recent treatment of various psychiatric disorders but not ADHD. Participants provided DNA and performed the BART, indexed by mean adjusted pumps. We constructed a polygenic risk score (PRS) for discovery in each dataset and tested it in the other as replication. Subsequently, a genome-wide MEGA-analysis, combining both samples, tested genetic correlation with risk-taking self-report in the UK Biobank sample and psychiatric phenotypes characterized by risk-taking (ADHD, Bipolar Disorder, Alcohol Use Disorder, prior cannabis use) in the Psychiatric Genomics Consortium. The PRS for BART performance in one dataset predicted task performance in the replication sample (r = 0.13, p = 0.000012, pFDR = 0.000052), as did the reciprocal analysis (r = 0.09, p = 0.0083, pFDR=0.04). Excluding participants with psychiatric diagnoses produced similar results. The MEGA-GWAS identified a single SNP (rs12023073; p = 3.24 × 10-8) near IGSF21, a protein involved in inhibitory brain synapses; replication samples are needed to validate this result. A PRS for self-reported cannabis use (p = 0.00047, pFDR = 0.0053), but not self-reported risk-taking or psychiatric disorder status, predicted behavior on the BART in our MEGA-GWAS sample. The findings reveal polygenic architecture of risky decision-making as measured by the BART and highlight its overlap with cannabis use.

Metabolic Exploration for Cyhalofop-Butyl Antagonism in Barnyardgrass [Echinochloa crus-galli (L.) P. Beauv.] Following Pretreatment of Malathion.

The present study investigated the effects of broad-spectrum metabolic inhibitors malathion (cytochrome P450 inhibitor) and/or 4-chloro-7-nitrobenzofurazan (NBD-Cl; glutathione S-transferase inhibitor) on the metabolism of cyhalofop-butyl (CyB) in barnyardgrass [Echinochloa crus-galli (L.) P. Beauv.] biotypes confirmed previously with multiple resistance to two herbicides CyB and florpyrauxifen-benzyl. The metabolic inhibitors were not effective at recovering the sensitivity of resistant barnyardgrass biotypes to CyB treated at the labeled rate (313 g ai ha-1). Rather, treatment with malathion followed by CyB caused antagonism, reducing the efficacy of CyB and promoting the growth of resistant biotypes. Pretreatment with malathion did not influence absorption/translocation of the applied form CyB and its conversion to the active herbicide form cyhalofop-acid (CyA), in both susceptible and resistant biotypes. In contrast, metabolism of the applied form (CyB) decreased 1.5 to 10.5 times by the malathion pretreatment. Taken together, the maintained CyA production against the reduced CyB metabolism could be the mechanism to account for the cause of CyB antagonism observed in barnyardgrass following malathion pretreatment. Additionally, the evolution of CyB resistance in barnyardgrass might be associated with reduced production of CyA in resistant biotypes, independent of activities of cytochrome P450 or GST enzymes.

Exploratory Analysis on Herbicide Metabolism and Very-Long-Chain Fatty Acid Production in Metolachlor-Resistant Palmer Amaranth (Amaranthus palmeri S. Wats.).

A Palmer amaranth (Amaranthus palmeri S. Wats.) biotype resistant to S-metolachlor was confirmed from crop fields in Arkansas, USA. This study investigated the metabolic effects of malathion (cytochrome P450 inhibitor) and 4-chloro-7-nitrobenzofurazan [NBD-Cl; glutathione S-transferase inhibitor] on the S-metolachlor-resistant A. palmeri biotype. Root elongation of the resistant biotype was 20% more inhibited by treatment of NBD-Cl (50 nM) and S-metolachlor (2 µM) in mixture than by treatment of S-metolachlor alone. Metabolites of S-metolachlor were 1.4-12.1 times greater produced in the resistant biotype for 7 d than in the susceptible standard. Production of cerotic acid, one of the very-long-chain fatty acids containing 26 carbons, was more reduced in the susceptible standard (3.8-fold) than in the resistant biotype (1.8-fold) by S-metolachlor treatment. Conclusively, evolution of S-metolachlor resistance observed in this study was likely associated with improved activity of glutathione S-transferases. Further studies are needed to genetically evaluate plant endogenous enzymes involving cerotic acid production.

Detection of subgenome bias using an anchored syntenic approach in Eleusine coracana (finger millet).

BACKGROUND: Finger millet (Eleusine coracana 2n = 4x = 36) is a hardy, nutraceutical, climate change tolerant, orphan crop that is consumed throughout eastern Africa and India. Its genome has been sequenced multiple times, but A and B subgenomes could not be separated because no published genome for E. indica existed. The classification of A and B subgenomes is important for understanding the evolution of this crop and provide a means to improve current and future breeding programs. RESULTS: We produced subgenome calls for 704 syntenic blocks and inferred A or B subgenomic identity for 59,377 genes 81% of the annotated genes. Phylogenetic analysis of a super matrix containing 455 genes shows high support for A and B divergence within the Eleusine genus. Synonymous substitution rates between A and B genes support A and B calls. The repetitive content on highly supported B contigs is higher than that on similar A contigs. Analysis of syntenic singletons showed evidence of biased fractionation showed a pattern of A genome dominance, with 61% A, 37% B and 1% unassigned, and was further supported by the pattern of loss observed among cyto-nuclear interacting genes. CONCLUSION: The evidence of individual gene calls within each syntenic block, provides a powerful tool for inference for subgenome classification. Our results show the utility of a draft genome in resolving A and B subgenomes calls, primarily it allows for the proper polarization of A and B syntenic blocks. There have been multiple calls for the use of phylogenetic inference in subgenome classification, our use of synteny is a practical application in a system that has only one parental genome available.

A novel mutation A212T in chloroplast Protoporphyrinogen oxidase (PPO1) confers resistance to PPO inhibitor Oxadiazon in Eleusine indica.

BACKGROUND: Protoporphyrinogen oxidase (PPO) with two isoforms, chloroplast-targeted (PPO1) and mitochondrial-targeted (PPO2), catalyzes a step in the biosynthesis of chlorophyll and heme. PPO1 and PPO2 are herbicide target sites of PPO-inhibiting herbicides. Target-site mutations conferring resistance to PPO inhibitors have all thus far been in PPO2. Oxadiazon is a unique PPO inhibitor utilized for preemergence Eleusine indica control. In this research, we evaluated the response of two previously confirmed oxadiazon-resistant and susceptible E. indica biotypes to other PPO inhibitors and identified the resistance mechanism in two oxadiazon-resistant E. indica biotypes. RESULTS: Two E. indica biotypes were resistant to oxadiazon, but not to other structurally unrelated PPO inhibitors, such as lactofen, flumioxazin and sulfentrazone. A novel mutation A212T was identified in the chloroplast-targeted PPO1, conferring resistance to oxadiazon in a heterologous expression system. Computational structural modeling provided a mechanistic explanation for reduced herbicide binding to the variant protein: the presence of a methyl group of threonine 212 changes the PPO1 active site and produces repulsive electrostatic interactions that repel oxadiazon from the binding pocket. CONCLUSION: The novel A212T mutation in PPO1 conferring resistance specifically to PPO inhibitor oxadiazon was characterized. This is the first evidence of the direct role of PPO1 in the PPO mode of action, and the first evidence of evolved resistance in PPO1. © 2019 Society of Chemical Industry.

Transcriptome Analysis Reveals Unique Relationships Among Eleusine Species and Heritage of Eleusine coracana.

Relationships in the genus Eleusine were obtained through transcriptome analysis. Eleusine coracana (E. coracana ssp. coracana), also known as finger millet, is an allotetraploid minor crop primarily grown in East Africa and India. Domesticated E. coracana evolved from wild E. africana (E. coracana ssp. africana) with the maternal genome donor largely supported to be E. indica; however, the paternal genome donor remains elusive. We developed transcriptomes for six Eleusine species from fully developed seedlings using Illumina technology and three de novo assemblers (Trinity, Velvet, and SOAPdenovo2) with the redundancy-reducing EvidentialGene pipeline. Mapping E. coracana reads to the chloroplast genes of all Eleusine species detected fewer variants between E. coracana and E. indica compared to all other species. Phylogenetic analysis further supports E. indica as the maternal parent of E. coracana and E. africana, in addition to a close relationship between E. indica and E. tristachya, and between E. floccifolia and E. multiflora, and E. intermedia as a separate group. A close relationship between E. floccifolia and E. multiflora was unexpected considering they are reported to have distinct nuclear genomes, BB and CC, respectively. Further, it was expected that E. intermedia and E. floccifolia would have a closer relationship considering they have similar nuclear genomes, AB and BB, respectively. A rethinking of the labeling of ancestral genomes of E. floccifolia, E. multiflora, and E. intermedia is maybe needed based on this data.

Rapid identification of a candidate nicosulfuron sensitivity gene (Nss) in maize (Zea mays L.) via combining bulked segregant analysis and RNA-seq.

KEY MESSAGE: A candidate nicosulfuron sensitivity gene Nss was identified by combining bulked segregant analysis and RNA-seq. Multiple mutations of this gene were discovered in nicosulfuron-sensitive maize compared with the tolerant. It has been demonstrated that variabilities exist in maize response to nicosulfuron. Two nicosulfuron-sensitive inbred lines (HB39, HB41) and two tolerant inbred lines (HB05, HB09) were identified via greenhouse and field trials. Genetic analysis indicated that the sensitivity to nicosulfuron in maize was controlled by a single, recessive gene. To precisely and rapidly map the nicosulfuron sensitivity gene (Nss), two independent F2 segregating populations, Population A (HB41 × HB09) and Population B (HB39 × HB05), were constructed. By applying bulked segregant RNA-Seq (BSR-Seq), the Nss gene was, respectively, mapped on the short arm of chromosome 5 (chr5: 1.1-15.3 Mb) and (chr5: 0.5-18.2 Mb) using two populations, with 14.2 Mb region in common. Further analysis revealed that there were 43 and 119 differentially expressed genes in the mapping intervals, with 18 genes in common. Gene annotation results showed that a cytochrome P450 gene (CYP81A9) appeared to be the candidate gene of Nss associated with nicosulfuron sensitivity in maize. Sequence analysis demonstrated that two common deletion mutations existed in the sensitive maize, which might lead to the nicosulfuron sensitivity in maize. The results will make valuable contributions to the understanding of molecular mechanism of herbicide sensitivity in maize.

The power and potential of genomics in weed biology and management.

There have been previous calls for, and efforts focused on, realizing the power and potential of weed genomics for better understanding of weeds. Sustained advances in genome sequencing and assembly technologies now make it possible for individual research groups to generate reference genomes for multiple weed species at reasonable costs. Here, we present the outcomes from several meetings, discussions, and workshops focused on establishing an International Weed Genomics Consortium (IWGC) for a coordinated international effort in weed genomics. We review the 'state of the art' in genomics and weed genomics, including technologies, applications, and on-going weed genome projects. We also report the outcomes from a workshop and a global survey of the weed science community to identify priority species, key biological questions, and weed management applications that can be addressed through greater availability of, and access to, genomic resources. Major focus areas include the evolution of herbicide resistance and weedy traits, the development of molecular diagnostics, and the identification of novel targets and approaches for weed management. There is increasing interest in, and need for, weed genomics, and the establishment of the IWGC will provide the necessary global platform for communication and coordination of weed genomics research. © 2018 Society of Chemical Industry.

Fast and reliable method for As speciation in urine samples containing low levels of As by LC-ICP-MS: Focus on epidemiological studies.

The speciation analysis of As in urine samples can provide important information for epidemiological studies. Considering that these studies involve hundreds or thousands of samples, a fast and reliable method using a simple LC system with short-length mixed bed ion exchange chromatographic column coupled to ICP-MS for As speciation in human urine samples was developed in this work. Separation of AB+TMAO, DMA, AC, MMA and AsIII+AsV was accomplished within 5min with good resolution when ammonium carbonate solutions were used as mobile phases and H2O2 was added to samples to quantitatively convert AsIII-AsV. Repeatability studies yielded RSD values from 2.0% to 4.8% for all species evaluated. Limits of detection (LOD) for As species ranged from 0.003 to 0.051ngg-1. Application of the method to human urine samples from a non-contaminated area showed that the sum of species measured corresponded to 62-125% of the total As in the sample. The recovery values for these species in urine SRM 2669 were in the range of 89-112% and demonstrated the suitability of the proposed method for epidemiological studies.

EMR documentation of physician-patient communication following genomic counseling for actionable complex disease and pharmacogenomic results.

Genomic risk information for potentially actionable complex diseases and pharmacogenomics communicated through genomic counseling (GC) may motivate physicians and patients to take preventive actions. The Ohio State University-Coriell Personalized Medicine Collaborative is a randomized trial to measure the effects of in-person GC on chronic disease patients provided with multiplex results. Nine personalized genomic risk reports were provided to patients through a web portal, and to physicians via electronic medical record (EMR). Active arm participants (98, 39% female) received GC within 1 month of report viewing; control arm subjects (101, 54% female) could access counseling 3-months post-report viewing. We examined whether GC affected documentation of physician-patient communication by reviewing the first clinical note following the patient's GC visit or report upload to the EMR. Multivariable logistic regression modeling estimated the independent effect of GC on physician-patient communication, as intention to treat (ITT) and per protocol (PP), adjusted for physician educational intervention. Counselees in the active arm had more physician-patient communications than control subjects [ITT, odds ratio (OR): 3.76 (95% confidence interval (CI): 1.38-10.22, p < 0.0094); PP, OR: 5.53 (95% CI: 2.20-13.90, p = 0.0017). In conclusion, GC appreciably affected physician-patient communication following receipt of potentially actionable genomic risk information.

Complete plastid genome sequence of goosegrass (Eleusine indica) and comparison with other Poaceae.

Eleusine indica, also known as goosegrass, is a serious weed in at least 42 countries. In this paper we report the complete plastid genome sequence of goosegrass obtained by de novo assembly of paired-end and mate-paired reads generated by Illumina sequencing of total genomic DNA. The goosegrass plastome is a circular molecule of 135,151bp in length, consisting of two single-copy regions separated by a pair of inverted repeats (IRs) of 20,919 bases. The large (LSC) and the small (SSC) single-copy regions span 80,667 bases and 12,646 bases, respectively. The plastome of goosegrass has 38.19% GC content and includes 108 unique genes, of which 76 are protein-coding, 28 are transfer RNA, and 4 are ribosomal RNA. The goosegrass plastome sequence was compared to eight other species of Poaceae. Although generally conserved with respect to Poaceae, this genomic resource will be useful for evolutionary studies within this weed species and the genus Eleusine.

Transcriptome Assembly and Comparison of an Allotetraploid Weed Species, Annual Bluegrass, with its Two Diploid Progenitor Species, Schrad and Kunth.

Annual bluegrass ( L.) is one of the most widespread weed species in this world. As a young allotetraploid, has occupied diverse environments from Antarctic area to subtropical regions. To unveil the evolutionary mystery behind 's wide distribution, extensive adaptability and phenotypic plasticity needs collaboration from multiple research scopes from ecology and plant physiology to population genetics and molecular biology. However, the lack of omic data and reference has greatly hampered the study. This is the first comprehensive transcriptome study on species. Total RNA was extracted from and its two proposed diploid parents, Schrad and Kunth, and sequenced in Illumina Hiseq2000. Optimized, nonredundant transcriptome references were generated for each species using four de novo assemblers (Trinity, Velvet, SOAPdenovo, and CLC Genomics Workbench) and a redundancy-reducing pipeline (CD-HIT-EST and EvidentialGene tr2aacds). Using the constructed transcriptomes together with sequencing reads, we found high similarity in nucleotide sequences and homeologous polymorphisms between and the two proposed parents. Comparison of chloroplast and mitochondrion genes further confirmed as the maternal parent. Less nucleotide percentage differences were observed between and homeologs than between and homeologs, indicating a higher nucleotide substitution rates in homeologs than in homeologs. Gene ontology (GO) enrichment analysis suggested the more compatible cytoplasmic environment and cellular apparatus for homeologs as the major cause for this phenomenon.

Optimizing Transcriptome Assemblies for Eleusine indica Leaf and Seedling by Combining Multiple Assemblies from Three De Novo Assemblers.

Due to rapid advances in sequencing technology, increasing amounts of genomic and transcriptomic data are available for plant species, presenting enormous challenges for biocomputing analysis. A crucial first step for a successful transcriptomics-based study is the building of a high-quality assembly. Here, we utilized three different de novo assemblers (Trinity, Velvet, and CLC) and the EvidentialGene pipeline tr2aacds to assemble two optimized transcript sets for the notorious weed species, Eleusine indica. Two RNA sequencing (RNA-seq) datasets from leaf and aboveground seedlings were processed using three assemblers, which resulted in 20 assemblies for each dataset. The contig numbers and N50 values of each assembly were compared to study the effect of read number, k-mer size, and in silico normalization on assembly output. The 20 assemblies were then processed through the tr2aacds pipeline to remove redundant transcripts and to select the transcript set with the best coding potential. Each assembly contributed a considerable proportion to the final transcript combination with the exception of the CLC-k14. Thus each assembler and parameter set did assemble better contigs for certain transcripts. The redundancy, total contig number, N50, fully assembled contig number, and transcripts related to target-site herbicide resistance were evaluated for the EvidentialGene and Trinity assemblies. Comparing the EvidentialGene set with the Trinity assembly revealed improved quality and reduced redundancy in both leaf and seedling EvidentialGene sets. The optimized transcriptome references will be useful for studying herbicide resistance in E. indica and the evolutionary process in the three allotetraploid E. indica offspring.

Utilizing next-generation sequencing to study homeologous polymorphisms and herbicide-resistance-endowing mutations in Poa annua acetolactate synthase genes.

BACKGROUND: Detection of single nucleotide polymorphisms (SNPs) related to herbicide resistance in non-model polyploid weed species is fraught with difficulty owing to the gene duplication and lack of reference sequences. Our research seeks to overcome these obstacles by Illumina HiSeq read mapping, SNP calling and allele frequency determinations. Our focus is on the acetolactate synthase (ALS) gene, the target site of ALS-inhibiting herbicides, in Poa annua, an allotetraploid weed species originating from two diploid parents, P. supina and P. infirma. RESULTS: ALS contigs with complete coding regions of P. supina, P. infirma and P. annua were assembled and compared with ALS genes from other plant species. The ALS infirma-homeolog of P. annua showed higher levels of nucleotide sequence variability than the supina-homeolog. Comparisons of read mappings of P. annua and a simulated P. supina × P. infirma hybrid showed high resemblance. Two homeolog-specific primer pairs were designed and used to amplify a 1860 bp region covering all resistance-conferring codons in the ALS gene. Four P. annua populations, GN, RB, GW and LG, showed high resistance to two ALS inhibitors, bispyribac-sodium and foramsulfuron, and two populations, HD and RS, showed lower resistance in the rate-response trial. Mutations conferring Trp-574-Leu substitution were observed in the infirma-homeolog of GN and RB and in the supina-homeolog of GW and LG, but no resistance-conferring mutation was observed in the two populations of lower resistance, HD and RS. CONCLUSION: In this study we have demonstrated the use of NGS data to study homeologous polymorphisms, parentage and herbicide resistance in an allotetraploid weed species, P. annua. Complete coding sequences of the ALS gene were assembled for P. infirma, P. supina, infirma-homeolog and supina-homeolog in P. annua. A pipeline consisting of read mapping, SNP calling and allele frequency calculation was developed to study the parentage of P. annua, which provided a new perspective on this topic besides the views of morphology, karyotype and phylogeny. Our two homeolog-specific primer pairs can be utilized in future research to separate the homeologs of the ALS gene in P. annua and cover all the codons that have been reported to confer herbicide resistance.

Effects of elevated CO2 on biomass and fungi associated with two ecotypes of ragweed (Ambrosia artemisiifolia L.).

Herbicide resistant weed populations have developed due to the repeated application of herbicides. Elevated concentrations of atmospheric CO2 can have positive effects on weed growth, but how rising CO2 might affect herbicide resistant weeds is not known. Ragweed (Ambrosia artemisiifolia L.) ecotypes known to be resistant or susceptible to glyphosate herbicide were exposed to either ambient or elevated (ambient +200 µ mol mol(-1)) concentrations of CO2 in open top chambers. Plants were harvested following 8 weeks of CO2 exposure; at this time, they had begun to exhibit disease symptoms including spots on leaves and stems. Elevated CO2 significantly increased top, root, and total plant biomass. Also, glyphosate resistant plants had significantly greater top, root, and total biomass than plants susceptible to the herbicide. There were no significant CO2 by ecotype interactions. Fungi from 13 genera were associated with ragweed, several of which can be either pathogens (i.e., Alternaria, Fusarium, Rhizoctonia), aiding the decline in health of the ragweed plants, or saprophytes existing on dead plant tissues. The common foliar disease powdery mildew was significantly higher on susceptible compared with resistant ragweed. Susceptible plants also showed an increased frequency of Rhizoctonia on leaves and Alternaria on stems; however, Fusarium occurred more frequently on resistant ragweed leaves. Fungi were not affected by CO2 concentration or its interaction with ecotype. This study reports the first information on the effects of elevated CO2 on growth of herbicide resistant weeds. This is also the first study examining the impact of herbicide resistance and elevated CO2 on fungi associated with weeds. What effects herbicide resistance might have on plant diseases and how rising atmospheric CO2 might impact these effects needs to be addressed, not only with important weeds but also with crops.

Analysis of the Citrullus colocynthis transcriptome during water deficit stress.

Citrullus colocynthis is a very drought tolerant species, closely related to watermelon (C. lanatus var. lanatus), an economically important cucurbit crop. Drought is a threat to plant growth and development, and the discovery of drought inducible genes with various functions is of great importance. We used high throughput mRNA Illumina sequencing technology and bioinformatic strategies to analyze the C. colocynthis leaf transcriptome under drought treatment. Leaf samples at four different time points (0, 24, 36, or 48 hours of withholding water) were used for RNA extraction and Illumina sequencing. qRT-PCR of several drought responsive genes was performed to confirm the accuracy of RNA sequencing. Leaf transcriptome analysis provided the first glimpse of the drought responsive transcriptome of this unique cucurbit species. A total of 5038 full-length cDNAs were detected, with 2545 genes showing significant changes during drought stress. Principle component analysis indicated that drought was the major contributing factor regulating transcriptome changes. Up regulation of many transcription factors, stress signaling factors, detoxification genes, and genes involved in phytohormone signaling and citrulline metabolism occurred under the water deficit conditions. The C. colocynthis transcriptome data highlight the activation of a large set of drought related genes in this species, thus providing a valuable resource for future functional analysis of candidate genes in defense of drought stress.

Copy number variation in pediatric multiple sclerosis.

BACKGROUND: Pediatric onset multiple sclerosis (MS) accounts for 2-4% of all MS. It is unknown whether the disease shares the same underlying pathophysiology found in adult patients or an extreme early onset phenotype triggered by distinct biological mechanisms. It has been hypothesized that copy number variations (CNVs) may result in extreme early onset diseases because CNVs can have major effects on many genes in large genomic regions. OBJECTIVES AND METHODS: The objective of the current research was to identify CNVs, with a specific focus on de novo CNVs, potentially causing early onset MS by competitively hybridizing 30 white non-Hispanic pediatric MS patients with each of their parents via comparative genomic hybridization (CGH) analysis on the Agilent 1M CGH array. RESULTS AND DISCUSSION: We identified 10 CNVs not overlapping with any CNV regions currently reported in the Database of Genomic Variants (DGV). Fifty-five putatively de novo CNVs were also identified: all but one common in the DGV. We found the single rare CNV was a private variation harboring the SACS gene. SACS mutations cause autosomal-recessive spastic ataxia of Charlevoix-Saguenay (ARSACS) disease. Additional clinical review revealed that the patient with the SACS gene CNV shared some features of both MS and ARSACS. CONCLUSIONS: This is the first reported study analyzing pediatric MS CNVs. While not yielding causal variation in our initial pediatric dataset, our approach confirmed diagnosis of an ARSACS-like disease in addition to MS in the affected individual, which led to a more complete understanding of the patient's disease course and prognosis.